Skip to main content

ORIGINAL RESEARCH article

Front. Mol. Biosci., 02 February 2024
Sec. Molecular Diagnostics and Therapeutics

Circulating miR-206, miR-181b, and miR-21 as promising biomarkers in hypothyroidism and their relationship to related hyperlipidemia and hepatic steatosis

Asmaa MohammedAsmaa Mohammed1Olfat G. ShakerOlfat G. Shaker2Mahmoud A. F. Khalil
Mahmoud A. F. Khalil3*Abeer K. Abu-El-AzayemAbeer K. Abu-El-Azayem4Amira SamyAmira Samy5Shaimaa A. FathyShaimaa A. Fathy6Mohamed M. K. AbdElguaadMohamed M. K. AbdElguaad7Fatma A. M. MahmoudFatma A. M. Mahmoud8Randa Erfan
Randa Erfan2*
  • 1Department of Medical Biochemistry and Molecular Biology, Faculty of Medicine, Fayoum University, Fayoum, Egypt
  • 2Department of Medical Biochemistry and Molecular Biology, Faculty of Medicine, Cairo University, Cairo, Egypt
  • 3Department of Microbiology and Immunology, Faculty of Pharmacy, Fayoum University, Fayoum, Egypt
  • 4Department of Medical Microbiology and Immunology, Faculty of Medicine, Cairo University, Cairo, Egypt
  • 5Department of Clinical and Chemical Pathology, Faculty of Medicine, Cairo University, Cairo, Egypt
  • 6Department of Internal Medicine, Diabetes and Endocrinology, Faculty of Medicine, Cairo University, Cairo, Egypt
  • 7Department of Physiology, Faculty of Medicine, Fayoum University, Fayoum, Egypt
  • 8Department of Tropical Medicine, Faculty of Medicine, Fayoum University, Fayoum, Egypt

Background: Thyroid hormones (THs) signaling has profound effects on many physiological processes. The regulation of THs signaling in various tissues involves the action of microRNAs (miRNAs) on thyroid deiodinases and receptors. THs regulate the expression of certain miRNAs and their target messenger RNAs (mRNAs) in various tissues and cells. The modulation of miRNA levels by THs affects their functions in processes such as liver lipid metabolism, skin physiology, and muscle and heart performance.

Aim: This research aimed to investigate miR-181b, miR-206, and miR-21 in the serum of patients with subclinical and overt hypothyroidism to determine their possible role in the diagnosis of the disease and their relationship to clinical disorders related to hypothyroidism.

Methods: This study included ninety participants, divided evenly into three groups as follows: patients with overt hypothyroidism diagnosed clinically, radiologically, and by investigation, subclinical hypothyroid patients, and healthy volunteers. The patients had a thorough medical history and underwent a clinical examination. Laboratory tests included plasma cholesterol, LDL, HDL, TGs, liver and renal function tests, CBC, fasting insulin, HOMA-IR, HbA1c, TSH, and free T4. The serum levels of miR-21, miR-206, and miR-181b were measured using qRT-PCR.

Results: miR-206 and miR-181b levels were higher in the subclinical group, followed by the hypothyroid and control groups. For miR-21, there was a significantly lower mean value in both the hypothyroid and subclinical groups than in the control group, with no difference between the two groups. Both miR-206 and miR-181b showed a significant negative association with albumin and free T4 levels and a significant direct association with GGT, ALT, AST, creatinine, uric acid, TGs, TC, LDL, TSH, thyroid volume, and CAP score. The same correlation pattern was observed for miR-181b, except that it was not significantly correlated with the TGs. For miR-21 levels, there was a significant positive correlation with albumin, free T4 level, and kPa score and a negative correlation with GGT, ALT, AST, creatinine, uric acid, HOMA-IR, HbA1c, TC, LDL, TSH, and CAP score. Cases with F1 kPa score and S2 CAP scores had significantly higher averages for miR-206 and miR-181b, with a p-value of 0.05. Moreover, miR-21 levels were significantly lower in the S2 CAP score group.

Conclusion: These miRNAs (miR-206, miR-181b, and miR-21) may be used as diagnostic biomarkers for hypothyroidism. They may be used as therapeutic targets to control dyslipidemia and hepatic steatosis during hypothyroid disease.

1 Introduction

Primary hypothyroidism is a prevalent condition that depends on the dietary iodine intake, age, and sex of the population. The high sensitivity of the hypothalamic-pituitary axis to changes in blood thyroid hormone (TH) levels makes the serum thyroid-stimulating hormone (TSH) level a reliable indicator for diagnosing this condition (Chaker et al., 2017). Subclinical hypothyroidism is an early stage of thyroid failure characterized by elevated TSH and normal TH levels. This stage was initially defined as being asymptomatic. However, when TH levels drop below normal, overt hypothyroidism develops along with severe symptoms (Hegedus et al., 2022). Various clinical conditions such as abnormal blood lipid levels (van Vliet et al., 2021), non-alcoholic fatty liver (NAFLD) (Teimouri et al., 2022), increased insulin resistance (assessed by HOMA-IR) (Shehata et al., 2022), type 2 diabetes mellitus (T2DM) (Han et al., 2022), and impaired kidney function (Ammar et al., 2021) are also associated with thyroid dysfunction.

MicroRNAs (miRNAs) play crucial roles in gene regulation. These are short RNAs of approximately 23 nucleotides, derived from a precursor with a double-stranded or hairpin structure. miRNA biogenesis involves RNA polymerase II (Landthaler et al., 2004). miRNAs can inhibit gene expression and translation via RISC, which binds to the 3′-UTRs of target mRNAs (Bartel, 2009). miRNAs are associated with various metabolic disorders, such as obesity and diabetes (Ji and Guo, 2019). Hypothyroidism, for example, is influenced by several miRNAs, such as miR-29a-3p (Tokic et al., 2018). miR-224-5p may also affect the development of hypothyroidism by targeting deiodinases, which convert T4 to rT3 (Wang et al., 2020).

A previous study demonstrated that hepatic expression of miR-206, miR-133a, and miR-133b is modulated by TH status in mice (Visser et al., 2009). It was also shown that serum levels of miR-206 are reduced in hyperthyroid patients and that this miRNA participates in the T3-mediated regulation of lipid metabolism in HepG2 cells, a human liver cancer cell line. This study suggested a novel molecular mechanism for TH-mediated lipid metabolism (Zheng et al., 2018). MiR-181b is another miRNA overexpressed in papillary thyroid cancer (PTC) cells (Nikiforova et al., 2008). Li et al. demonstrated that miR-181b downregulation inhibited cell proliferation and induced apoptosis by increasing the expression of CYLD, a lysine 63 deubiquitinase.

Additionally, miR-181b expression was approximately eight-fold higher in tumor tissues than in normal tissues (Li et al., 2014). Moreover, miR-181b overexpression has been linked to an increased potential for cancer recurrence and lymph node metastasis (Pamedytyte et al., 2019). MiR-21 is a key miRNA involved in PTC pathogenesis. Its expression is altered in tumor tissues (Samsonov et al., 2016). Ortiz et al. revealed that hypomethylation of DNA leads to overexpression of miR-141b and miR-21, which results in reduced transcription of their target genes (Ortiz et al., 2018). Sondermann et al. revealed a positive association between low miR-21 expression and an increased PTC recurrence rate. The authors also identified miR-21 and miR-9 as strong predictors of PTC recurrence (Sondermann et al., 2015). In our study, we chose three miRNAs that are regulated by TH (miR-206, miR-181b, and miR-21) to examine their levels in patients with hypothyroidism and subclinical hypothyroidism. This is the first study to examine the fold changes of miR-206, miR-181b, and miR-21 in the serum of patients with subclinical and overt hypothyroidism to detect their possible role in the diagnosis of the disease and their relationship with the development of associated clinical disorders.

2 Methods

2.1 Ethical statement

This study followed the ethical principles established by the pertinent national and institutional committees on human experimentation in line with the revised 2008 Helsinki Declaration. Approval for all procedures involving human subjects or patients (code: R384) was obtained from the Research Ethics Committee of the General Organization for Teaching Hospitals and Institutes (GOTHI). All the patients provided written informed consent before participating in the study. The manuscript was prepared according to the Strengthening the Reporting of Observational Studies in Epidemiology (STROBE) guidelines (Moher et al., 2012).

2.2 Study design and patient selection

This research involved three groups. The first group included 30 patients with overt hypothyroidism diagnosed clinically, radiologically, or by investigation. The second group comprised 30 subclinical hypothyroid patients, and the last group consisted of 30 healthy volunteers.

The following definitions were used to define overt hypothyroidism and subclinical hypothyroidism patients:

The term “subclinical hypothyroidism” is used to define that grade of primary hypothyroidism in which there is an elevated TSH concentration in the presence of normal serum-free thyroxine (T4) and triiodothyronine (T3) concentrations. Subclinical hypothyroidism may progress to overt hypothyroidism in approximately 2%–5% of cases annually. At the same time, overt hypothyroidism is clear hypothyroidism characterized by an increased TSH and a decreased T4 level (Khandelwal and Tandon, 2012).

All patients were subjected to full history and clinical examination. Full history taking covered asking about symptoms suggesting hypothyroidism, including Fatigue, loss of energy, lethargy, weight gain, decreased appetite, cold intolerance, dry skin, hair loss, sleepiness, muscle pain, joint pain, weakness in the extremities, depression, emotional lability, mental impairment, forgetfulness, impaired memory, inability to concentrate, constipation, menstrual disturbances, impaired fertility, decreased perspiration, paresthesias, nerve entrapment syndromes, blurred vision, decreased hearing, fullness in the throat and hoarseness (Spanou et al., 2019).

Regarding physical examination, vital signs, including blood pressure, pulse, respiratory rate, and temperature, were assessed. A detailed general examination, head and neck, chest, heart, abdominal and neurological examination were performed to look for the following signs: weight gain, slowed speech and movements, dry skin (rarely, yellow-hued from carotene), Jaundice, pallor, coarse, brittle, straw-like hair, loss of scalp hair, axillary hair, pubic hair, or a combination, dull facial expression, coarse facial features, periorbital puffiness, macroglossia, goiter (simple or nodular), hoarseness, decreased systolic blood pressure and increased diastolic blood pressure, bradycardia, pericardial effusion, abdominal distention, ascites (uncommon), hypothermia (only in severe hypothyroid states), nonpitting edema (myxedema), pitting edema of lower extremities. Hyporeflexia with delayed relaxation, ataxia, or both (Piantanida et al., 2016).

2.3 Diagnosis of liver steatosis and fibrosis

Transient elastography (TE) using a FibroScan® device (Echosens, Paris) was carried out to investigate liver stiffness (LS). TE was performed with a standard M probe, and an XL probe (for obese patients) (Myers et al., 2012). Investigations were performed through the intercostal spaces, where patients were lying in the dorsal decubitus position with the right arm in maximal abduction. Measurements were taken following an overnight fast. The tip of the probe is placed in contact with the intercostal skin through a coupling gel in the ninth to 11th intercostal space at the level where a liver biopsy would be conducted. After locating a section of the liver that is at least 6 cm thick and devoid of significant vascular structures, the operator initiates shots by pressing the probe button.

Measurement depth was between 25 and 65 mm below the skin surface. The success or failure of each measurement is determined by the software. The subsequent parameters are applied for validating successful measurements:

1) Number of shots ≥10.

2) Ratio of valid shots to the total number of shots ≥60%.

3) Interquartile range (IQR) less than 30% of the median liver stiffness measurement (LSM) value (IQR/LSM ≤30%).

The results are described in kilopascal (kPa). The median value of the successful measurements was representative of LS (Castéra et al., 2010).

In NAFLD population:

F0, no fibrosis à <5.5 kPa.

F1, mild fibrosis à 5.5–8.0 kPa.

F2, moderate fibrosis à 8.0–10.0 kPa.

F3, severe fibrosis à 11–16 kPa.

F4, cirrhosis à >16 kPa (Kamali et al., 2019a).

CAP score:

A FibroScan® device (Echosens, Paris) was performed to measure LS with the assessment of hepatic steatosis using CAP. TE was carried out with a standard M probe. CAP measurement depends on the attenuation of US propagation in liver tissue due to the presence of fat and is calculated via a sophisticated process. CAP is measured only on validated measurements of TE according to the same criteria used for LSM and on the same signals. This ensures that the operator obtains a liver ultrasonic attenuation simultaneously and in the same volume of liver parenchyma as the LSM. The cut-off points used to stage steatosis were S1 ≥ 222 dB/m, S2 ≥233 dB/m, and S3 ≥290 dB/m (Sasso et al., 2012).

In NAFLD population:

S0, no steatosis à <237 dB/m.

S1, mild steatosis à 237–259 dB/m.

S2, moderate steatosis à 259–291 dB/m.

S3, severe steatosis à 291–400 dB/m (Kamali et al., 2019b).

2.4 Laboratory investigation

Liver functions included alanine aminotransferase (ALT), total bilirubin (diazo), albumin (BCP), gamma-glutamyl transferase (GGT) (enzymatic), and aspartate aminotransferase (AST), according to the methodology recommended by the International Federation of Clinical Chemistry (IFCC). Renal functions included Urea (enzymatic), Creatinine (Jaffé), and uric acid (uricase) on the Beckman Coulter AU analyzer (Beckman Coulter Inc., CA 92821, United States). The assessment of cholesterol levels followed the accuracy criteria of the National Cholesterol Education Program (NCEP), and enzymatic techniques were used on the same instrument to analyze low-density lipoprotein (LDL), triglycerides, and high-density lipoprotein (HDL).

Fasting insulin using a two-site sandwich immunochemiluminometric assay (ICMA) on ADVIA Centaur (Siemens Healthcare, 91052 Erlangen, Germany) and fasting blood glucose (hexokinase) were used to calculate the homeostatic model assessment for insulin resistance (HOMA-IR) formula. Hba1c (hemoglobin a1c) measurement was performed using high-performance liquid chromatography (Ion Exchange HPLC) on a Tosoh G8 Analyzer (TOSOH Bioscience, US), which is the gold standard method for Hba1c measurement.

A direct chemiluminescence TSH assay using third-generation acridinium ester technology (detection limit less than 0.01 mU/L) was used to measure TSH concentrations. Serum TSH and free T4 (fT4) levels were analyzed to assess thyroid function using a Siemens ADVIA Centaur. The TSH assay is a two-site sandwich method, and the FT4 test is a competitive immunoassay. A 10-mL blood sample was obtained from each patient using a vacutainer device. The blood samples were stored in tubes containing separator gels. The serum layer (top) and packed cells were separated using gels. After 15 min, the tubes were centrifuged at 4,000 g for 10 min. The serum was stored at −80 °C after separating it from the clotted blood and then used for RNA extraction (Bush et al., 2001).

2.5 RNA extraction

RNA extraction was conducted following the protocol provided by the miRNeasy extraction kit (Qiagen, Valencia, CA, USA) with a total sample volume of 100 µL of serum. Initially, the reaction mixture was treated with the QIAzol lysis reagent (500 µL) and allowed to stand for 5 min at RT. Subsequently, 100 µL chloroform was added, followed by vortexing for 15 s and incubation for 2–3 min at RT. Centrifugation was performed at 4 °C for 15 min at 120,00× g. Upon removing the top aqueous phase, 1.5 times the volume of 100% ethanol was added to the remaining mixture. Each 700 µL of this solution was then centrifuged for 15 s at 8,000× g and RT on an RNeasy Mini spin column and placed in a 2 mL collection tube. The mixture was eluted through each column, and 700 µL of buffer RW1 was added. Centrifugation was performed at RT for 15 s at a speed of 8,000× g. The column received 500 µL of buffer RPE and then underwent centrifugation at RT at 8,000× g for 15 s. To further purify the mixture, the same amount of buffer RPE was added, and the column was centrifuged for 2 min at RT at the same speed. Using a pipette, 50 µL of RNase-free water was applied directly to the column to elute RNA. Centrifugation was performed for 1 minute at 8,000× g. A DNase Max Kit (Qiagen, Valencia, CA, United States) was used to perform DNase treatment to remove any DNA before RNA was converted into cDNA by reverse transcription. A NanoDrop 2000 spectrophotometer (Thermo Scientific, Waltham, MA, United States) was used to measure the concentration and purity of the RNA at 260/280 nm.

2.6 Reverse transcription

A high-capacity cDNA reverse transcription kit (Applied Biosystems, Foster City, CA, USA) was used to convert 1 µg of RNA into cDNA following the manufacturer’s instructions. The reaction was incubated at 37 °C for 60 min, denatured at 95 °C for 5 min, and stored at 4 °C.

2.7 miRNA expression by real-time quantitative PCR

qRT-PCR was used to measure serum levels of miR181b-5p and miR21-5p. The miRCURY LNA SYBR Green PCR Kit (Qiagen, Valencia, CA, United States), miRCURY LNA miRNA PCR Assay (Qiagen, Valencia, CA, United States), and Rotor-Gene Q RT-PCR system (Qiagen, Germany) were used for analysis. Predesigned primers from Qiagen (Valencia, CA, United States) were used for target miRNAs and SNORD 68 (reference gene) (Ayeldeen et al., 2018), and primer assays for miR206 (catalog no. MS00003787), miR181b (catalog no. YP00204530), miR21 (Catalog no. YP00204230), and SNORD68 (catalog no. MS00033712). A modified thermal cycle was performed, with an initial denaturation step of 2 min at 95°C, followed by 40cycles. Each cycle consisted of 10 s denaturation at 95°C and 60 s of annealing/extension at 56°C. The 2−ΔΔCT equation was used to calculate the serum fold-change of miR181b-5p and miR21-5p for each sample (Livak and Schmittgen, 2001).

Ct, ΔCt, ΔΔCt, and Fold Change (FC) have the following mathematical relationships:

ΔCt (patients) = Ct (miRNA)–Ct (housekeeping gene).

ΔCt (control) = Ct (miRNA)–Ct (housekeeping gene).

ΔΔCt (patients) = ΔCt (patients)–ΔCt (Control).

FC = 2−ΔΔCTmiRNA overexpression was indicated when the calculated fold change was greater than one, whereas a fold change of less than one signified miRNA downregulation. In this context, because the -ΔΔCt value for control individuals is equal to zero and 2° equals one, the control value was assumed to be 1 (Santoro et al., 2016).

3 Results

There were no statistically significant differences (p > 0.05) in the demographic characteristics of the study groups (Table 1).

TABLE 1
www.frontiersin.org

TABLE 1. Comparison of different study groups based on demographic characteristics.

3.1 Laboratory investigations

The CBC investigations did not reveal any statistically significant differences (p > 0.05) among the study groups. Creatinine and uric acid levels were significantly higher in hypothyroidism and subclinical cases than in controls, with no differences between the two patient groups. Mean ALT and HOMA-IR levels were significantly higher in the hypothyroidism group than in the subclinical and control groups. GGT and HbA1c levels differed among the three study groups, with a higher mean in the hypothyroid group and a lower mean among the controls. AST showed a significant difference between the hypothyroid group and the controls, with a higher mean in the hypothyroid group. The hypothyroid and subclinical groups of cases had significantly higher mean TGs, TC, and LDL levels than the controls, with no difference between case subgroups. In addition, there was a significantly lower mean HDL level in the hypothyroid group than in the control group; however, there was no difference between the subclinical group and the controls. The hypothyroid and subclinical groups showed a significantly higher mean TSH level than the controls, with no difference between the case subgroups. In addition, there was a significant difference between all groups regarding the level of free T4, with a lower mean among the hypothyroids, followed by the subclinical group, compared with the controls. Regarding thyroid volume, there was a significant difference between all three groups, with a higher mean among the hypothyroids, followed by the subclinical group (Table 2).

TABLE 2
www.frontiersin.org

TABLE 2. Laboratory investigations comparing the study groups.

3.2 Radiological findings

The hypothyroid and subclinical groups of cases showed a significantly higher mean CAP score and a higher percentage of F1 levels by FibroScan than controls, with no difference between case subgroups. In addition, there was a significant difference between all three groups, with a higher percentage of S3 CAP scores among patients with hypothyroidism and S2 and S3 CAP scores among the subclinical group compared to controls who showed the S0 CAP scores (Table 3).

TABLE 3
www.frontiersin.org

TABLE 3. Comparison of study groups based on radiological findings.

3.3 Serum levels of biomarkers

Our results demonstrated a significant difference (p < 0.05) between the three study groups concerning miR-206 and miR-181b, with a higher level among the subclinical group, followed by the hypothyroid group, and then the control group. For miR-21, there was a significantly lower mean value among both the hypothyroid and subclinical groups than in the control group, with no difference between the two groups (Table 4; Figures 1A–C).

TABLE 4
www.frontiersin.org

TABLE 4. Comparison of markers.

FIGURE 1
www.frontiersin.org

FIGURE 1. (A) Serum miR-206 levels in hypothyroid, subclinical, and control groups. (B) Serum miR-181b levels in hypothyroid, subclinical, and control groups. (C) Serum levels of miR-21 in the hypothyroid, subclinical, and control groups.

3.4 Comparison of biomarker levels in different demographic and radiological data among the patient group

For cases with an F1 kPa score and S2 CAP score, the patient group had a substantially greater mean value of miR-206 and miR-181b (p = 0.05). Additionally, miR-21 levels were significantly lower in the S2 CAP scores. However, there were no sex differences for any of these markers (Table 5).

TABLE 5
www.frontiersin.org

TABLE 5. Comparison of markers in different demographic and radiological data among the patient group.

3.5 Correlation between biomarker levels and study variables

The patient group demonstrated a significant negative association between albumin, miR-206, and free T4 levels and a significant positive association between miR-206 and several parameters, such as GGT, ALT, AST, creatinine, uric acid, HOMA-IR, HbA1c, TGs, TC, LDL, TSH, thyroid volume, and CAP score. Similarly, miR-181b showed a significantly negative association with albumin and free T4 levels and a significantly positive association with various parameters, such as GGT, ALT, AST, creatinine, uric acid, HbA1c, TGs, TC, LDL, TSH, thyroid volume, and cap score. For miR-21 level, there was a significant positive correlation with albumin, free T4 level, and kPa score and a negative correlation with GGT, ALT, AST, creatinine, uric acid, HOMA-IR, HbA1c, TC, LDL, TSH, and CAP score.

The expression levels of miR-206 and miR-181b showed a strong positive association, while miR-21 exhibited a strong negative association with miR-181b and miR-21, as indicated by the correlation coefficients in Table 6.

TABLE 6
www.frontiersin.org

TABLE 6. Relationship between study variables and biomarker levels in the patient group.

3.6 Detection of predictive power of miR-206, miR-181b, and miR-21 in hypothyroidism and subclinical cases

The sensitivity and specificity test for the diagnosis of hypothyroidism showed higher sensitivity and specificity among the miR-206, miR-181b, and miR-21 marker levels, with a sensitivity of 96.7%, 93.3%, and 83.3%, respectively, and a specificity of 96.7%, 96.7%, and 90% at a cut-off value of 1.55, 1.45, and 0.8 respectively.

For subclinical cases, the sensitivity test for miR-206, miR-181b, and miR-21 showed a sensitivity of 100% and 83.3% and specificity of 96.7% and 100% at a cut-off value of 1.55 and 0.77, respectively (Table 7; Figures 2A–D).

TABLE 7
www.frontiersin.org

TABLE 7. Diagnostic accuracy of marker levels for hypothyroidism and its subclinical cases measured by sensitivity and specificity.

FIGURE 2
www.frontiersin.org

FIGURE 2. (Continued).

4 Discussion

THs modulate the expression levels of certain miRNAs and their target mRNAs in different cell types and organs. Many miRNAs contain thyroid hormone elements (TREs) in their regulatory regions, which are transcriptionally regulated by hormone-bound receptors. THs exert significant effects on processes, including lipid and liver metabolism, muscle and heart function, and skin physiology (Aranda, 2021), by altering the levels of miRNAs. Among these miRNAs, we selected miR-206, miR-181b, and miR-21, which have been reported to affect lipid metabolism in multiple disorders (Wang et al., 2017; Zheng et al., 2018; Zhao et al., 2021).

Our findings demonstrated that miRNAs in patients with thyroid illness were dysregulated compared to those in healthy controls. Regarding miR-206, a significant difference was found between the study groups (p < 0.05), with the subclinical group having the highest levels, followed by the hypothyroid and control groups. miR-206 was found to have a significantly higher mean among cases with F1 kPa score and S2 CAP score. Moreover, a significant inverse association between miR-206 and albumin and free T4 levels was observed in the patient group, as was a positive correlation between miR-206 and several biochemical markers (GGT, ALT, AST, creatinine, uric acid, HOMA-IR, HbA1c, TGs, TC, LDL, TSH, thyroid volume, and CAP score). Recent research has reported the contribution of thyroid hormones to the regulation of miR-206 in the mouse liver and corresponding cell line models (Dong et al., 2010) and in human skeletal muscle in hypothyroid patients (Visser et al., 2009). Thyroid hormones and their analogs reduce miR-206 levels, whereas hypothyroidism enhances its expression (Visser et al., 2009; Dong et al., 2010). Therefore, miR-206 is also suggested to act as a key mediator of thyroid hormone function since some of its targets (e.g., Mup1 and Gpd2) were shown to be modulated by thyroid hormones (Zhao et al., 2021). Zheng and colleagues showed the involvement of miR-206 in lipid metabolism in hyperthyroidism and hypothyroidism by finding that serum miR-206 expression decreases in hyperthyroid patients. They suggested that miR-206 was involved in regulating lipid metabolism by T3 in HepG2 cells, implying that miR-206 may play a role in thyroid hormone-related liver lipid metabolism disorders.

Additionally, miR-206 overexpression in HepG2 cells increased lipid accumulation and attenuated the T3-mediated decrease in TGs, and hepatic steatosis was associated with TH deficiency (Zheng et al., 2018). This previous study is consistent with ours, as we observed a positive correlation between miR-206 and GGT, ALT, AST, TGs, TC, and LDL levels. In addition, our findings confirmed the role of this miRNA in hepatic steatosis during hypothyroidism, as we found a positive correlation between miR-206 and the CAP score, which measures the percentage of fatty change in the liver. Additionally, we found a significantly higher mean miR-206 level in patients with an F1 kPa score than in those with an F0 kPa score, which measures the degree of liver stiffness.

Our findings revealed an inverse association between miR-206 and free T4 levels and a positive association between miR-206 and TSH levels. These results support those of previous studies, suggesting the influence of thyroid hormones on regulating miR-206 expression. In agreement with Mugdha et al., our results align with their finding that there is no relationship between the transcripts miR-206 and glucokinase in human islets, contrasting with the situation observed in mouse islets (Joglekar et al., 2018). The results reported by Vinod et al., who suggested that miR-206 regulated glucokinase mRNA post-transcriptionally in mouse pancreatic islets, indicating a possible role of miR-206 in glucose homeostasis regulation, possibly through glucokinase (Vinod et al., 2016), are inconsistent with our findings. In our study, a positive association was found between miR-206 and HbA1c as well as HOMA-IR, which further supports the potential involvement of miR-206 in metabolic processes.

Regarding miR-181b, we found that subclinical hypothyroid patients had the highest levels, followed by overt hypothyroid patients. Similar to our findings, Huang et al. discovered that thyroid hormones downregulate several miRNAs, including miR-200, miR-92a, and miR-181b (Huang et al., 2019).

The role of miR-181b in NAFLD and its direct interaction with NAD-dependent deacetylase sirtuin1 (SIRT1) were examined by Wang et al. They found that human patients with NAFLD and NAFLD mouse models induced by high-fat diet (HFD) had higher miR-181b levels and lower SIRT1 levels. They also showed that miR-181b directly targeted SIRT1 in HepG2 cells and that miR-181b inhibition reduced hepatic steatosis both in vitro and in vivo. Furthermore, SIRT1 overexpression reversed the effects of miR-181b on steatosis. Their data suggested that miR-181b may be involved in lipid metabolism dysregulation in NAFLD (Wang et al., 2017). Moreover, they discovered that the SIRT1, a direct target of miR-181b, modulated TGs production in human liver cells and influenced lipid metabolism. This previous study may explain our results, as we demonstrated a significant positive correlation between miR-181b and GGT, ALT, AST, HbA1c, TGs, TC, LDL, and CAP scores. In addition, miR-181b had a significantly higher mean in cases with F1 kPa score and S2 CAP score.

For miR-21, there was a significantly lower mean value in both the hypothyroid and subclinical groups than in the control group, with no difference between the two groups. For miR-21 levels, there was a significant positive correlation with albumin, free T4 level, and kPa score and a negative correlation with GGT, ALT, AST, creatinine, uric acid, HOMA-IR, HbA1c, TC, LDL, TSH, and CAP score. The present study agrees with the results of Huang et al., who reported that miR-21 expression is regulated by TH in vivo and in vitro. They showed that miR-21 levels were lower in the livers of rats with hypothyroidism and higher in rats with hyperthyroidism. They also demonstrated that tumor growth and metastasis of a TRα-overexpressing cell line were influenced by the thyroid status of the mice, with hyperthyroid mice having higher miR-21 and lower TIAM Rac1-associated GEF 1 (TIAM1) expression than euthyroid or hypothyroid mice. This confirmed that TH modulates miR-21 expression and its downstream target, TIAM1 (Huang et al., 2013).

Moreover, miR-21 has been widely investigated in thyroid cancers and has been identified as one of the miRNAs that could improve the diagnostic accuracy of PTC and is associated with TNM stage, metastasis, and tumor size. It can be used as a diagnostic marker to distinguish PTC from benign lesions (Wang et al., 2019; Samsonov et al., 2016; Zhang et al., 2017). MiR-21 was negatively correlated with GGT, ALT, AST, TC, LDL, and CAP scores. This finding aligns with prior research indicating downregulation of miR-21 expression in individuals affected by non-alcoholic fatty liver disease or in mice subjected to a HFD. Additionally, ablation of miR-21 has been observed to induce hepatic steatosis, expedited atherosclerosis, plaque necrosis, and vascular inflammation. These observations are consistent with and reinforce the conclusions of previous studies (Sun et al., 2015). Zhou et al. recently showed that miR-21 can be used in treating hyperlipidemia, as exercise also affects the expression of miR-21 and has beneficial effects on hyperlipidemia (Zhao et al., 2021). The sensitivity and specificity tests for diagnosing hypothyroidism showed a higher level of sensitivity and specificity among the miR-206 marker levels, with a sensitivity of 96.7% and a specificity of 96.7% at a cut-off value of 1.55. For subclinical disease, the sensitivity test for miR-206 and miR-181b showed the same level (100%) of sensitivity and specificity (96.7%) at a cut-off value of 1.55.

5 Conclusion

Finally, it was concluded that these miRNAs (miR-206, miR-181b, and miR-21) can be used as diagnostic biomarkers for hypothyroidism. TH is a significant regulator of these biomarkers, which influence related disorders such as dyslipidemia and hepatic steatosis; therefore, they may be used as therapeutic targets for controlling the previously mentioned disorders during hypothyroid disease. Further studies are required to validate these results.

Data availability statement

The original contributions presented in the study are included in the article/Supplementary material, further inquiries can be directed to the corresponding authors.

Ethics statement

The studies involving humans were approved by the Institutional Review Board of the Faculty of Medicine, Fayoum University Ethical Committee. The studies were conducted in accordance with the local legislation and institutional requirements. The participants provided their written informed consent to participate in this study.

Author contributions

AM: Conceptualization, Formal Analysis, Methodology, Writing–original draft. OS: Conceptualization, Investigation, Methodology, Supervision, Writing–original draft. MK: Conceptualization, Writing–original draft, Writing–review and editing. AA-E-A: Supervision, Writing–review and editing, Formal Analysis, Methodology, Investigation. AS: Formal Analysis, Methodology, Data curation, Investigation, Writing–original draft. SF: Methodology, Validation, Visualization, Writing–original draft. MA: Methodology, Formal Analysis, Investigation, Writing–review and editing. FM: Formal Analysis, Methodology, Writing–review and editing. RE: Methodology, Conceptualization, Supervision, Writing–original draft.

Funding

The author(s) declare that no financial support was received for the research, authorship, and/or publication of this article.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

References

Ammar, A., Mir, K. B., Batool, S., Marwat, N., Saeed, M., Saeed, A., et al. (2021). Assessment of decrease in renal function associated with hypothyroidism. World J. Adv. Res. Rev. 12 (1), 275–286. doi:10.30574/wjarr.2021.12.1.0477

CrossRef Full Text | Google Scholar

Aranda, A. (2021). MicroRNAs and thyroid hormone action. Mol. Cell Endocrinol. 525, 111175. doi:10.1016/j.mce.2021.111175

PubMed Abstract | CrossRef Full Text | Google Scholar

Ayeldeen, G., Nassar, Y., Ahmed, H., Shaker, O., and Gheita, T. (2018). Possible use of miRNAs-146a and -499 expression and their polymorphisms as diagnostic markers for rheumatoid arthritis. Mol. Cell Biochem. 449 (1-2), 145–156. doi:10.1007/s11010-018-3351-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Bartel, D. P. (2009). MicroRNAs: target recognition and regulatory functions. Cell 136 (2), 215–233. doi:10.1016/j.cell.2009.01.002

PubMed Abstract | CrossRef Full Text | Google Scholar

Bush, V. J., Janu, M. R., Bathur, F., Wells, A., and Dasgupta, A. (2001). Comparison of BD Vacutainer SST™ plus tubes with BD SST™ II plus tubes for common analytes. Clin. Chim. Acta. 306 (1-2), 139–143. doi:10.1016/s0009-8981(01)00396-5

PubMed Abstract | CrossRef Full Text | Google Scholar

Castéra, L., Foucher, J., Bernard, P. H., Carvalho, F., Allaix, D., Merrouche, W., et al. (2010). Pitfalls of liver stiffness measurement: a 5-year prospective study of 13,369 examinations. Hepatology 51 (3), 828–835. doi:10.1002/hep.23425

PubMed Abstract | CrossRef Full Text | Google Scholar

Chaker, L., Bianco, A. C., Jonklaas, J., and Peeters, R. P. (2017). Hypothyroidism. Lancet 390 (10101), 1550–1562. doi:10.1016/S0140-6736(17)30703-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Dong, H., Paquette, M., Williams, A., Zoeller, R. T., Wade, M., and Yauk, C. (2010). Thyroid hormone may regulate mRNA abundance in liver by acting on microRNAs. PLoS One 5 (8), e12136. doi:10.1371/journal.pone.0012136

PubMed Abstract | CrossRef Full Text | Google Scholar

Han, M., Wu, H., Yang, W., and Chen, J. (2022). Analysis of risk factors for the development of type 2 diabetes mellitus complicated with Hashimoto's thyroiditis. BMC Endocr. Disord. 22 (1), 173. doi:10.1186/s12902-022-01092-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Hegedus, L., Bianco, A. C., Jonklaas, J., Pearce, S. H., Weetman, A. P., and Perros, P. (2022). Primary hypothyroidism and quality of life. Nat. Rev. Endocrinol. 18 (4), 230–242. doi:10.1038/s41574-021-00625-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Huang, P. S., Wang, C. S., Yeh, C. T., and Lin, K. H. (2019). Roles of thyroid hormone-associated microRNAs affecting oxidative stress in human hepatocellular carcinoma. Int. J. Mol. Sci. 20 (20), 5220. doi:10.3390/ijms20205220

PubMed Abstract | CrossRef Full Text | Google Scholar

Huang, Y. H., Lin, Y. H., Chi, H. C., Liao, C. H., Liao, C. J., Wu, S. M., et al. (2013). Thyroid hormone regulation of miR-21 enhances migration and invasion of hepatoma. Cancer Res. 73 (8), 2505–2517. doi:10.1158/0008-5472.CAN-12-2218

PubMed Abstract | CrossRef Full Text | Google Scholar

Ji, C., and Guo, X. (2019). The clinical potential of circulating microRNAs in obesity. Nat. Rev. Endocrinol. 15 (12), 731–743. doi:10.1038/s41574-019-0260-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Joglekar, M. V., Wong, W. K. M., Maynard, C. L., Umrani, M. R., Martin, D., Loudovaris, T., et al. (2018). Expression of miR-206 in human islets and its role in glucokinase regulation. Am. J. Physiol. Endocrinol. Metab. 315 (4), E634–E637. doi:10.1152/ajpendo.00152.2018

PubMed Abstract | CrossRef Full Text | Google Scholar

Kamali, F., Zamanlou, M., Ghanbari, A., Alipour, A., and Bervis, S. (2019b). Comparison of manipulation and stabilization exercises in patients with sacroiliac joint dysfunction patients: a randomized clinical trial. J. Bodyw. Mov. Ther. 23 (1), 177–182. doi:10.1016/j.jbmt.2018.01.014

PubMed Abstract | CrossRef Full Text | Google Scholar

Kamali, L., Adibi, A., Ebrahimian, S., Jafari, F., and Sharifi, M. (2019a). Diagnostic performance of ultrasonography in detecting fatty liver disease in comparison with fibroscan in people suspected of fatty liver. Adv. Biomed. Res. 8, 69. doi:10.4103/abr.abr_114_19

PubMed Abstract | CrossRef Full Text | Google Scholar

Khandelwal, D., and Tandon, N. (2012). Overt and subclinical hypothyroidism: who to treat and how. Drugs 72 (1), 17–33. doi:10.2165/11598070-000000000-00000

PubMed Abstract | CrossRef Full Text | Google Scholar

Landthaler, M., Yalcin, A., and Tuschl, T. (2004). The human DiGeorge syndrome critical region gene 8 and its D. melanogaster homolog are required for miRNA biogenesis. Curr. Biol. 14 (23), 2162–2167. doi:10.1016/j.cub.2004.11.001

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, D., Jian, W., Wei, C., Song, H., Gu, Y., Luo, Y., et al. (2014). Down-regulation of miR-181b promotes apoptosis by targeting CYLD in thyroid papillary cancer. Int. J. Clin. Exp. Pathol. 7 (11), 7672–7680.

PubMed Abstract | Google Scholar

Livak, K. J., and Schmittgen, T. D. (2001). Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) Method. Methods 25 (4), 402–408. doi:10.1006/meth.2001.1262

PubMed Abstract | CrossRef Full Text | Google Scholar

Moher, D., Hopewell, S., Schulz, K. F., Montori, V., Gotzsche, P. C., Devereaux, P. J., et al. (2012). CONSORT 2010 explanation and elaboration: updated guidelines for reporting parallel group randomised trials. Int. J. Surg. 10 (1), 28–55. doi:10.1016/j.ijsu.2011.10.001

PubMed Abstract | CrossRef Full Text | Google Scholar

Myers, R. P., Pomier-Layrargues, G., Kirsch, R., Pollett, A., Duarte-Rojo, A., Wong, D., et al. (2012). Feasibility and diagnostic performance of the FibroScan XL probe for liver stiffness measurement in overweight and obese patients. Hepatology 55 (1), 199–208. doi:10.1002/hep.24624

PubMed Abstract | CrossRef Full Text | Google Scholar

Nikiforova, M. N., Tseng, G. C., Steward, D., Diorio, D., and Nikiforov, Y. E. (2008). MicroRNA expression profiling of thyroid tumors: biological significance and diagnostic utility. J. Clin. Endocrinol. Metab. 93 (5), 1600–1608. doi:10.1210/jc.2007-2696

PubMed Abstract | CrossRef Full Text | Google Scholar

Ortiz, I., Barros-Filho, M. C., Dos Reis, M. B., Beltrami, C. M., Marchi, F. A., Kuasne, H., et al. (2018). Loss of DNA methylation is related to increased expression of miR-21 and miR-146b in papillary thyroid carcinoma. Clin. Epigenetics 10 (1), 144. doi:10.1186/s13148-018-0579-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Pamedytyte, D., Simanaviciene, V., Dauksiene, D., Leipute, E., Zvirbliene, A., Kazokaite, M., et al. (2019). SAT-LB096 relationship between MicroRNA expression levels and clinicopathological parameters and recurrence of papillary thyroid carcinoma. J. Endocr. Soc. 3 (1), SAT-LB096. doi:10.1210/js.2019-SAT-LB096

CrossRef Full Text | Google Scholar

Piantanida, E., Gallo, D., Veronesi, G., Pariani, N., Masiello, E., Premoli, P., et al. (2016). Masked hypertension in newly diagnosed hypothyroidism: a pilot study. J. Endocrinol. Invest. 39 (10), 1131–1138. doi:10.1007/s40618-016-0488-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Samsonov, R., Burdakov, V., Shtam, T., Radzhabovsmall a, C. Z., Vasilyev, D., Tsyrlina, E., et al. (2016). Plasma exosomal miR-21 and miR-181a differentiates follicular from papillary thyroid cancer. Tumour Biol. 37 (9), 12011–12021. doi:10.1007/s13277-016-5065-3

PubMed Abstract | CrossRef Full Text | Google Scholar

Santoro, M., Nociti, V., Lucchini, M., De Fino, C., Losavio, F. A., and Mirabella, M. (2016). Expression profile of long non-coding RNAs in serum of patients with multiple sclerosis. J. Mol. Neurosci. 59 (1), 18–23. doi:10.1007/s12031-016-0741-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Sasso, M., Tengher-Barna, I., Ziol, M., Miette, V., Fournier, C., Sandrin, L., et al. (2012). Novel controlled attenuation parameter for noninvasive assessment of steatosis using Fibroscan(®): validation in chronic hepatitis C. J. Viral Hepat. 19 (4), 244–253. doi:10.1111/j.1365-2893.2011.01534.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Shehata, A. E., Elfayoumy, K. N., Emran, T. M., and Eladl, H. E. (2022). Study of insulin resistance in patients with subclinical hypothyroidism. J Int. J. Med. Arts 4 (11), 2794–2800. doi:10.21608/ijma.2022.137559.1453

CrossRef Full Text | Google Scholar

Sondermann, A., Andreghetto, F. M., Moulatlet, A. C., da Silva Victor, E., de Castro, M. G., Nunes, F. D., et al. (2015). MiR-9 and miR-21 as prognostic biomarkers for recurrence in papillary thyroid cancer. Clin. Exp. Metastasis 32 (6), 521–530. doi:10.1007/s10585-015-9724-3

PubMed Abstract | CrossRef Full Text | Google Scholar

Spanou, I., Bougea, A., Liakakis, G., Rizonaki, K., Anagnostou, E., Duntas, L., et al. (2019). Relationship of migraine and tension-type headache with hypothyroidism: a literature review. Headache 59 (8), 1174–1186. doi:10.1111/head.13600

PubMed Abstract | CrossRef Full Text | Google Scholar

Sun, C., Huang, F., Liu, X., Xiao, X., Yang, M., Hu, G., et al. (2015). miR-21 regulates triglyceride and cholesterol metabolism in non-alcoholic fatty liver disease by targeting HMGCR. Int. J. Mol. Med. 35 (3), 847–853. doi:10.3892/ijmm.2015.2076

PubMed Abstract | CrossRef Full Text | Google Scholar

Teimouri, K., Pakravan, S., and Azadbakht, K. (2022). Examining the relationship between non-alcoholic fatty liver disease and hypothyroidism. J. Parathyr. Dis. 10 (1), e9148. doi:10.34172/jpd.2022.9148

CrossRef Full Text | Google Scholar

Tokic, S., Stefanic, M., Glavas-Obrovac, L., Kishore, A., Navratilova, Z., and Petrek, M. (2018). miR-29a-3p/T-bet regulatory circuit is altered in T cells of patients with hashimoto's thyroiditis. Front. Endocrinol. (Lausanne) 9, 264. doi:10.3389/fendo.2018.00264

PubMed Abstract | CrossRef Full Text | Google Scholar

van Vliet, N. A., Bos, M. M., Thesing, C. S., Chaker, L., Pietzner, M., Houtman, E., et al. (2021). Higher thyrotropin leads to unfavorable lipid profile and somewhat higher cardiovascular disease risk: evidence from multi-cohort Mendelian randomization and metabolomic profiling. BMC Med. 19 (1), 266. doi:10.1186/s12916-021-02130-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Vinod, M., Patankar, J. V., Sachdev, V., Frank, S., Graier, W. F., Kratky, D., et al. (2016). MiR-206 is expressed in pancreatic islets and regulates glucokinase activity. Am. J. Physiol. Endocrinol. Metab. 311 (1), E175–E185. doi:10.1152/ajpendo.00510.2015

PubMed Abstract | CrossRef Full Text | Google Scholar

Visser, W. E., Heemstra, K. A., Swagemakers, S. M., Ozgur, Z., Corssmit, E. P., Burggraaf, J., et al. (2009). Physiological thyroid hormone levels regulate numerous skeletal muscle transcripts. J. Clin. Endocrinol. Metab. 94 (9), 3487–3496. doi:10.1210/jc.2009-0782

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, C., Zhu, J., Zhang, Z., Chen, H., Ji, M., Chen, C., et al. (2020). Rno-miR-224-5p contributes to 2,2',4,4'-tetrabromodiphenyl ether-induced low triiodothyronine in rats by targeting deiodinases. Chemosphere 246, 125774. doi:10.1016/j.chemosphere.2019.125774

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, L., Duan, Y., Peng, W., Qu, C., Lin, J., Deng, Z., et al. (2019). miR-21 facilitates the diagnostic value of miR-138 for papillary thyroid cancer in formalin-fixed paraffin-embedded tissues. Transl. Cancer Res. 8 (5), 1718–1726. doi:10.21037/tcr.2019.08.14

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, Y., Zhu, K., Yu, W., Wang, H., Liu, L., Wu, Q., et al. (2017). MiR-181b regulates steatosis in non-alcoholic fatty liver disease via targeting SIRT1. Biochem. Biophys. Res. Commun. 493 (1), 227–232. doi:10.1016/j.bbrc.2017.09.042

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, Y., Xu, D., Pan, J., Yang, Z., Chen, M., Han, J., et al. (2017). Dynamic monitoring of circulating microRNAs as a predictive biomarker for the diagnosis and recurrence of papillary thyroid carcinoma. Oncol. Lett. 13 (6), 4252–4266. doi:10.3892/ol.2017.6028

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhao, J., Song, Y., Zeng, Y., Chen, L., Yan, F., Chen, A., et al. (2021). Improvement of hyperlipidemia by aerobic exercise in mice through a regulatory effect of miR-21a-5p on its target genes. Sci. Rep. 11 (1), 11966. doi:10.1038/s41598-021-91583-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Zheng, Y., Zhao, C., Zhang, N., Kang, W., Lu, R., Wu, H., et al. (2018). Serum microRNA miR-206 is decreased in hyperthyroidism and mediates thyroid hormone regulation of lipid metabolism in HepG2 human hepatoblastoma cells. Mol. Med. Rep. 17 (4), 5635–5641. doi:10.3892/mmr.2018.8633

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: thyroid hormone, hepatic steatosis, microRNA, miR-206, MiR-181b, miR-21

Citation: Mohammed A, Shaker OG, Khalil MAF, Abu-El-Azayem AK, Samy A, Fathy SA, AbdElguaad MMK, Mahmoud FAM and Erfan R (2024) Circulating miR-206, miR-181b, and miR-21 as promising biomarkers in hypothyroidism and their relationship to related hyperlipidemia and hepatic steatosis. Front. Mol. Biosci. 11:1307512. doi: 10.3389/fmolb.2024.1307512

Received: 04 October 2023; Accepted: 08 January 2024;
Published: 02 February 2024.

Edited by:

Federico Ávila-Moreno, National Autonomous University of Mexico, Mexico

Reviewed by:

Ali Sahin Kucukaslan, Beykent University, Türkiye
Tapahsama Banerjee, University of Colorado, United States

Copyright © 2024 Mohammed, Shaker, Khalil, Abu-El-Azayem, Samy, Fathy, AbdElguaad, Mahmoud and Erfan. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Mahmoud A. F. Khalil, maf04@fayoum.edu.eg; Randa Erfan, randa.erfan@cu.edu.eg

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.