ARTICLE   Open Access    

The Association between Vitamin D and Hashimoto Thyroiditis: An Up-to-date Systematic Review and Meta-analysis

More Information
  • The objective of the surrent study was to summarize the up-to-date studies to investigate the relationship between vitamin D and Hashimoto thyroiditis (HT). An online search of English and Chinese databases was performed. The studies concerned the investigation of the relationship between vitamin D and HT including meta-analysis, meanwhile the heterogeneities were revealed by subgroup analysis. Fourty six elated studies containing 15,336 participants (HT: 6,138 versus control: 9,198) were included. HT patients had lower levels of 25(OH)D3 (standardised mean difference, −1.09; 95%CI: [−1.42, −0.75]; P < 0.01), and were more likely to be deficient in 25(OH)D3 (OR, 2.77; 95%CI, [1.88, 3.91]; P < 0.05). Obvious heterogeneities in the results of meta-analysis were down to the difference of detection methods and criteria of vitamin D insufficiency among studies. Vitamin D deficiency was colncluded to have a significant relation with HT.
  • The Lonicera Linn. genus is a constituent member of the Caprifoliaceae family[1]. It is the largest genus in this family and comprises at least 200 species with a notable presence in North Africa, North America, Asia, and Europe[1]. Members of the Lonicera genus possess a wide range of economic benefits from their use as ornamental plants to food and as plants credited with numerous health benefits. Conspicuous among the numerous members of this genus with known medicinal uses are L. japonica, L. macranthoides, L. hypoglauca, L. confusa, and L. fulvotomentosa[2]. Though these species feature prominently in the Chinese Pharmacopoeia, other species such as L. acuminata, L. buchananii, and L. similis are recognized as medicinal resources in certain parts of China[1]. Among the aforementioned species, L. japonica takes precedence over the rest due to its high medicinal and nutritional value[3,4]. For instance, the microRNA MIR2911, an isolate from L. japonica, has been reported to inhibit the replication of viruses[57]. Also, the water extract of L. japonica has been used to produce various beverages and health products[8]. The Lonicera genus therefore possesses huge prospects in the pharmaceutical, food, and cosmetic industries as an invaluable raw material[9].

    The main active constituents of the Lonicera genus include organic acids, flavonoids, iridoids, and triterpene saponins. Chlorogenic acids, iridoids, and flavones are mainly credited with the anti-inflammatory, antiviral, anticancer, and antioxidant effects of the various Lonicera species[1013]. Their hepatoprotective, immune modulatory, anti-tumor and anti-Alzheimer’s effects are for the most part ascribed to the triterpene saponins[1416]. As stated in the Chinese Pharmacopoeia and backed by the findings of diverse research groups, the plants of the Lonicera genus are known to possess high amounts of organic acids (specifically chlorogenic acid) and pentacyclic triterpenoid saponins[2,1719]. The flower and flower bud have traditionally served as the main medicinal parts of the Lonicera genus even though there is ample evidence that the leaves possess the same chemical composition[20]. A perusal of the current scientific literature reveals the fact that little attention has been devoted to exploring the biosynthesis of the chemical constituents of the Lonicera genus with the view to finding alternative means of obtaining higher yields. It is therefore imperative that priority is given to the exploration of the biosynthesis of these bioactive compounds as a possible means of resource protection. There is also the need for further research on ways to fully tap the medicinal benefits of other plant parts in the Lonicera genus.

    Here, we provide a comprehensive review of relevant scientific literature covering the structure, pharmacology, multi-omics analyses, phylogenetic analysis, biosynthesis, and metabolic engineering of the main bioactive constituents of the Lonicera genus. Finally, we proffer suggestions on the prospects of fully exploiting and utilizing plants of the Lonicera genus as useful medicinal plant resources.

    A total of at least 400 secondary metabolites have been reported for the Lonicera genus. These metabolites are categorized into four main groups (Fig. 1a), including not less than 50 organic acids, 80 flavonoids, 80 iridoids, and 80 triterpene saponins[2123]. Organic acids are mainly derivatives of p-hydroxycinnamic acid and quinic acid. Among the organic acids, chlorogenic acids are reported to be the main bioactive compounds in L. japonica[2426]. The organic acids are most abundant in the leaves, while the least amounts are found in the stem of L. japonica. The flowers of the plant are known to contain moderately high amounts of organic acids[27]. The basic core structure of the flavonoids is 2-phenylchromogen. Luteolin and its glycoside which are characteristic flavonoids of the Lonicera genus are most abundant in L. japonica[28]. On the whole, the flavonoid contents in L. japonica are also highest in the leaves, available in moderate amounts in the flowers, and in lowest amounts in the stem[21]. The core structures of the iridoids are iridoid alcohols, the chemical properties of which are similar to hemiacetal. The iridoids often exist in the form of iridoid glycosides in plants. Secoiridoids glycosides are predominant in the Lonicera genus[25]. In L. japonica, the contents of the iridoids are most abundant in the flowers, moderate in leaves, and lowest in the stem[21]. The characteristic saponins of the Lonicera genus are mainly pentacyclic triterpenoids, including the hederin-, oleanane-, lupane-, ursulane- and fernane-types, etc[22]. The hederin-type saponins are reported in the highest amounts in L. macranthoides[17] (Fig. 1b).

    Figure 1.  Core structures of main secondary metabolites and their distribution in five species of Lonicera. (a) 1 and 2, the main core structures of organic acids; 3, the main core structures of flavonoids; 4, the main core structures of iridoids; 5, the main core structures of triterpene saponins. (b) Comparison of dry weight of four kinds of substances in five species of Lonicera[17,28].

    The similarities between chlorogenic acid (CGA) and flavonoids can be traced back to their biosynthesis since p-coumaroyl CoA serves as the common precursor for these compounds[29]. p-coumaroyl CoA is obtained through sequential catalysis of phenylalanine and its biosynthetic intermediates by phenylalanine-ammonia-lyase (PAL), cinnamate 4-hydroxylase (C4H) and 4-coumarate CoA ligase (4CL)[3033].

    CGA is a phenolic acid composed of caffeic acid and quinic acid and is the most important bioactive compound among the organic acids. Its biosynthesis has been relatively well-established; three main biosynthetic routes have been propounded (Fig. 2a). One route relates to the catalysis of caffeoyl-CoA and quinic acid by hydroxycinnamoyl-CoA quinate transferase (HQT)/hydroxycinnamoyl CoA shikimate/quinate hydroxycinnamoyl transferase (HCT) to produce CGA[3437]. The HQT-mediated pathway has been deemed the major route for CGA synthesis in in different plant species[38,39]. The second biosynthetic route stems from the biosynthesis of p-coumaroyl quinate through the catalytic effect of HCT/HQT and subsequent hydroxylation of p-coumaroyl quinate under the catalysis of p-coumarate 3'-hydroxylase (C3’H)[34,36,37]. For the third route, caffeoyl glucoside serves as the intermediate to form CGA, a process that is catalyzed by hydroxycinnamyl D-glucose: quinic acid hydroxycinnamyl transferase (HCGQT)[40,41].

    Figure 2.  Biosynthetic pathways of main bioactive constituents of Lonicera. (a) Biosynthetic pathways of chlorogenic acid. (b) Biosynthetic pathways of luteoloside. (c) Biosynthetic pathways of secologanin. (d) Biosynthetic pathways of hederin-type triterpene saponins. PAL, phenylalanine ammonia-lyase; C4H, cinnamate 4-hydroxylase; 4CL, 4-hydroxycinnamoyl CoA ligase; HCT, hydroxycinnamoyl CoA shikimate/quinate hydroxycinnamoyl transferase; C3’H, p-coumaroyl 3-hydroxylase; HQT, hydroxycinnamoyl-CoA quinate transferase; UGCT, UDP glucose: cinnamate glucosyl transferase; CGH, p-coumaroyl-D-glucose hydroxylase; HCGQT, hydroxycinnamoyl D-glucose: quinate hydroxycinnamoyl transferase; CHS, Chalcone synthase; CHI, Chalcone isomerase; FNS, Flavone synthase; F3H, flavonoid 30-monooxygenase/flavonoid 30-hydroxylase; UF7GT, flavone 7-O-β-glucosyltransferase; GPS, Geranyl pyrophosphatase; GES, geraniol synthase; G8O, geraniol 10-hydroxylase/8-oxidase; 8HO, 8-hydroxygeraniol oxidoreductase; IS, iridoid synthase; IO, iridoid oxidase; 7DLGT, 7-deoxyloganetic acid glucosyltransferase; 7DLH, 7-deoxyloganic acid hydroxylase; LAMT, loganic acid O-methyltransferase; SLS, secologanin synthase; FPS, farnesyl pyrophosphate synthase; SS, squalene synthase; SE, squalene epoxidase; β-AS, β-amyrin synthase; OAS, oleanolic acid synthetase.

    The key enzymes in the biosynthesis of p-coumaroyl CoA, and invariably CGA, thus, PAL, C4H, and 4CL have been established in diverse studies such as enzyme gene overexpression/knockdown[42], enzyme activity analysis[33] and transcriptomics[18]. However, the centrality of HQT in the biosynthesis of CGA remains disputable. While some studies have reported a strong correlation between HQT expression level with CGA content and distribution[18,35,39,43,44], others found no such link[45], bringing into question the role of HQT as a key enzyme in CGA biosynthesis.

    Few studies have been conducted on the regulation of CGA biosynthesis in the Lonicera genus. It was found that overexpression of the transcription factor, LmMYB15 in Nicotiana benthamiana can promote CGA accumulation by directly activating 4CL or indirectly binding to MYB3 and MYB4 promoters[46]. LjbZIP8 can specifically bind to PAL2 and act as a transcriptional repressor to reduce PAL2 expression levels and CGA content[47]. Under NaCl stress, increased PAL expression promoted the accumulation of phenolic substances in leaves without oxidative damage, a condition that was conducive to the accumulation of the bioactive compounds in leaves[48].

    Luteolin and its derivative luteolin 7-O- glucoside (luteoloside) are representative flavonoids of the Lonicera genus. Similar to CGA, luteolin is biosynthesized from p-coumaroyl CoA but via a different route. The transition from p-coumaroyl CoA to luteolin is underpinned by sequential catalysis by chalcone synthetase (CHS), chalcone isomerase (CHI), flavonoid synthetase (FNS), and flavonoid 3'-monooxygenase/flavonoid 3'-hydroxylase (F3'H)[45,49,50] (Fig. 2b). Luteoloside is synthesized from luteolin by UDP glucose-flavonoid 7-O-β-glucosyltransferase (UF7GT)[51]. Similar to CGA biosynthesis, the key enzymes of luteolin synthesis include PAL, C4H, and 4CL in addition to FNS[33,45,52]. The content of luteoloside was found to be highly abundant in senescing leaves relative to other tissues such as stem, flowers, and even young leaves[52]. Through transcriptomic analysis, luteoloside biosynthesis-related differentially expressed unigenes (DEGs) such as PAL2, C4H2, flavone 7-O-β-glucosyltransferase (UFGT), 4CL, C4H, chalcone synthase 2 and flavonoid 3'-monooxygenase (F3'H) genes were found to be upregulated in the senescing leaves. Biosynthesis-related transcription factors such as MYB, bHLH, and WD40 were also differentially expressed during leaf senescence[52], while bHLH, ERF, MYB, bZIP, and NAC were differentially expressed during flower growth[53]. Further analysis of the transcription factors revealed that MYB12, MYB44, MYB75, MYB114, MYC12, bHLH113, and TTG1 are crucial in luteoloside biosynthesis[52,53].

    The biosynthesis of terpenoids mainly involves three stages; formation of intermediates, formation of basic structural skeleton, and modification of basic skeleton[54]. The intermediates of terpenoids are mainly formed through the mevalonate (MVA) and methylerythritol phosphate (MEP) pathways, and eventually converted to the universal isoprenoid precursors, isopentenyl pyrophosphate (IPP) and its isomer dimethylallyl pyrophosphate (DMAPP) through a series of enzyme-catalyzed reactions. Under the catalysis of geranyl pyrophosphatase (GPS), IPP is then converted to geranyl pyrophosphate (GPP). Different terpenoids are subsequently derived from GPP as the intermediate product. For instance, in the formation of secoiridoid, GPP first removes the phosphoric acid group to obtain geraniol, second through a series of reactions such as oxidation and cyclization, the skeleton of iridoid, namely iridodial, can be obtained. Finally, through a series of reactions, the basic carbon skeleton of the secoiridoid, namely secologanin, is obtained[5561] (Fig. 2c). In the formation of triterpene saponins, the key step lies in the formation of the precursor, 2,3-oxidosqualene, a reaction that is catalyzed by squalene epoxidase (SE). There are many pentacyclic triterpenes in the Lonicera genus, the most important type being the hederin-type saponins with hederagenin as aglycones. Hederin-type saponins are produced after the synthesis of oleanolic acid from β-amyrin and catalyzed by β-starch synthetase (β-AS) and Oleanolic acid synthase (OAS)[62,63]. The skeletal modification of the triterpenoid saponins is mainly achieved via the activities of the CYP450 enzymes and UDP-glycosyltransferase (UGT). Hence, the corresponding aglycones are first obtained via oxidation by the CYP450 enzymes (e.g., CYP72A), and further subjected to glycosylation by the appropriate UGT enzyme[6365] (Fig. 2d). Skeletal formations of the iridoids and triterpene saponins in general have been well researched, but the same cannot be said about the enzymes involved in biosynthesis of these groups of compounds in the Lonicera genus. To fully utilize the iridoids and triterpene saponins in the Lonicera genus, it is necessary to further explore their biosyntheses with the view to enhancing and optimizing the process.

    Given the importance of the bioactive compounds in the Lonicera genus, continual isolation of these compounds using the traditional methods are not only tedious and time-consuming, but also unsustainable. With the development and application of microbial metabolic engineering, different strategies have been introduced to produce these bioactive compounds by heterologous synthesis (Table 1).

    Table 1.  Biosynthesis of Lonicera-specialized metabolites using metabolic engineering.
    Engineering bacteriaOperational methodsProductsYieldRefs
    S. cerevisiaeEliminate the tyrosine-induced feedback inhibition, delete genes involved in competing pathways and overexpress rate-limiting enzymesCaffeic acid569.0 mg/L[69]
    S. cerevisiaeEmploye a heterologous tyrosine ammonia lyase and a 4HPA3H complex composed of HpaB and HpaC derived from different speciesCaffeic acid289.4 mg/L[73]
    S. cerevisiaeSupply and recycle of three cofactors: FADH2, S-adenosyl-L-methion, NADPHCaffeic acid
    Ferulic acid
    Caffeic acid: 5.5 g/L;
    Ferulic acid: 3.8 g/L
    [117]
    E. coliKnocking out competing pathwaysCaffeic acid7,922 mg/L[118]
    E. coliArtificial microbial community, a polyculture of three recombinant Escherichia coli strainsChlorogenic acid250 μM[68]
    Cell-free biosynthesisExtract and purify spy-cyclized enzymes (CFBS-mixture)Chlorogenic acid711.26 mg/L[70]
    S. cerevisiaeThree metabolic engineering modules were systematically optimized: shikimate pathway and carbon distribution, branch pathways, CGA pathway genesChlorogenic acidFlask fermentation: 234.8 mg/L;
    Fed-batch fermentation:
    806.8 mg/L
    [119]
    E. coliUsing modular coculture engineering: construction of the defective strain improves the production and utilization of precursor substancesChlorogenic acid131.31 mg/L[122]
    E. coliIntroduce heterologous UDP-glucose biosynthetic genesLuteolin34 mg/L[120]
    Y. lipolyticaOverexpression of the key genes involved in the mevalonate pathway, the gene encoding cytochrome P450 (CYP716A12) to that encoding NADPH-P450 reductaseOleanolic acid129.9 mg/L[85]
    S. cerevisiaeImprove the pairing efficiency between Cytochrome P450 monooxygenase and reductase and the expression level of key genesOleanolic acid606.9 mg/L[121]
    S. cerevisiaeHeterologous expression and optimization of CrAS, CrAO, and AtCPR1, and regulation of ERG1 and NADPH regeneration systemOleanolic acid433.9 mg/L[123]
     | Show Table
    DownLoad: CSV

    Due to the demand for CGA in the food, pharmaceutical, chemical, and cosmetic industries, the traditional means of obtaining the same requires a relatively longer period for plant maturation to obtain low yields of the desired product. This therefore brings into question the sustainability and efficiency of this approach. The alternative and sustainable approach has been to produce CGA using synthetic biology and metabolic engineering.

    Current research has sought to utilize Escherichia coli (and its mutant strain) and Saccharomyces cerevisiae to synthetically generate CGA and other flavonoids[6673]. For instance, Cha et al. employed two strains of E. coli to produce a relatively good yield of CGA (78 mg/L). Their approach was based on the ability of one strain to generate caffeic acid from glucose and the other strain to use the caffeic acid produced and quinic acid as starting materials to synthesize CGA[66]. Using a bioengineered mutant of E. coli (aroD mutant), Kim et al. increased the yield of CGA to as high as 450 mg/L[67]. Others have sought to increase the yield of CGA by employing a polyculture of three E. coli strains that act as specific modules for the de novo biosynthesis of caffeic acid, quinic acid and CGA. This strategy eliminates the competition posed by the precursor of CGA (i.e., caffeic acid and quinic acid) and generally results in improved production of CGA[68]. Saccharomyces cerevisiae is a chassis widely used for the production of natural substances from plants with an intimal structure that can be used for the expression of cytochrome P450 enzymes that cannot be expressed in E. coli. Researchers have used yeast to increase the production of organic acids[69]. A de novo biosynthetic pathway for the construction of CGA in yeast has been reported new cell-free biosynthetic system based on a mixture of chassis cell extracts and purified Spy cyclized enzymes were adopted by Niu et al. to a produce the highest yield of CGA reported so far up to 711.26 ± 15.63 mg/L[70].

    There are many studies on the metabolic engineering for the synthesis of flavonoids, but few on luteolin and its glycosides. Strains of E. coli have been engineered with specific uridine diphosphate (UDP)-dependent glycosyltransferase (UGT) to synthesize three novel flavonoid glycosides. These glycosides were quercetin 3-O-(N-acetyl) quinovosamine (158.3 mg/L), luteolin 7-O-(N-acetyl) glucosaminuronic acid (172.5 mg/L) and quercetin 3-O-(N-acetyl)-xylosamine (160.8 mg/L)[71]. Since most of the flavonoid glycosides synthesized in E. coli are glucosylated, Kim et al. in their bid to synthesize luteolin-7-O-glucuronide, deleted the araA gene that encodes UDP-4-deoxy-4-formamido-L-arabinose formyltransferase/UDP-glucuronic acid C-4'' decarboxylase in E.coli and were able to obtain a yield of 300 mg/L of the desired product[72].

    Terpenoidal saponins are mostly derived from slow-growing plants and usually possess multiple chiral centers[74]. Traditional isolation and even chemical synthesis of the terpenoidal saponins are both tedious and uneconomical for large-scale production. Therefore, it is necessary to find other ways to synthesize these compounds known to have diverse pharmacological functions.

    Heterologous synthesis has become an important way to improve the target products. With the development of synthetic biology, heterologous synthesis of triterpene saponins involves chassis of both plant and microbial origin. In this regard, Nicotiana benthamiana is a model plant species for the reconstruction of the biosynthetic pathways of different bioactive compounds including monoterpenes, hemiterpenes, and diterpenes[59,7577]. Aside from Nicotiana benthamiana, other plants have also been used as heterologous hosts[78]. Heterologous synthesis using microbial hosts mainly involves Saccharomyces cerevisiae and Escherichia coli[7981], and other microorganisms[82,83]. Comparatively, plants as biosynthetic hosts have the advantages of an established photosynthetic system, abundant supply of relevant enzymes, and presence of cell compartments, etc. They are however not as fast growing as the microorganisms, and it is also difficult to extract and separate the desired synthesized compounds from them as hosts.

    Although heterologous synthesis has many advantages, the premise of successful construction of synthetic pathway in host is to elucidate the unique structure of the compound and the key enzyme reaction mechanism in the biosynthetic pathway. There is little research on metabolic engineering of the hederin-type pentacyclic triterpene saponins in Lonicera, but there are studies on the heterologous synthesis of its aglycone precursor, oleanolic acid[84,85]. There is a dearth of scientific literature on key enzymes in the biosynthesis of pentacyclic triterpenoid saponins in the Lonicera genus.

    Scientific evidence by diverse research groups has linked members of the Lonicera genus to a wide range of pharmacological effects (Fig. 3). These pharmacological effects are elicited by different chemical constituents, much of the underlying mechanisms of which have been elucidated by the omics techniques. Here, we summarize the pharmacological effects and pharmacodynamics of the Lonicera genus in the last 6 years.

    Figure 3.  Schematic summary of four main pharmacological effects (anti-inflammatory, antimicrobial, anti-oxidative and hepatoprotective effects) of the Lonicera genus and the underlying mechanisms of actions.

    Bioactive compounds of plants in the Lonicera genus have demonstrated varying degrees of anti-inflammatory actions. In a recent study, Lv et al. showed that lonicerin inhibits the activation of NOD-like receptor thermal protein domain associated protein 3 (NLRP3) through regulating EZH2/AtG5-mediated autophagy in bone marrow-derived macrophages of C57BL/6 mice[86]. The polysaccharide extract of L. japonica reduces atopic dermatitis in mice by promoting Nrf2 activation and NLRP3 degradation through p62[87]. Several products of Lonicera have been reported to have ameliorative effects on DSS-induced colitis. Among them, flavonoids of L. rupicola can improve the ulcerative colitis of C57BL/6 mice by inhibiting PI3K/AKT, and pomace of L. japonica can improve the ulcerative colitis of C57BL/6 mice by improving the intestinal barrier and intestinal flora[88,89]. The flavonoids can also ameliorate ulcerative colitis induced by local enema of 2,4,6-trinitrobenzene sulfonic acid (TNBS) in Wistar rats by inhibiting NF-κB pathway[90]. Ethanol extract from L. Japonica has demonstrated the potential to inhibit the expressions of inflammatory cytokines in serum and macrophages of LPS-induced ICR mice[91]. The water extract of L. japonica and luteolin were found to exhibit their anti-inflammatory effects via the inhibition of the JAK/STAT1/3-dependent NF-κB pathway and induction of HO-1 expression in RAW263.7 cells induced by pseudorabies virus (PRV)[92].

    Existing scientific evidence indicates that the extracts of plants in the Lonicera genus exhibit strong inhibition against different pathogenic microorganisms. Phenolic compounds from L. japonica demonstrated a particularly significant inhibitory effect against Staphylococcus aureus and Escherichia coli, in vitro, making these compounds potential food preservatives[93]. Influenza A virus is a serious threat to human health. Recent research has found the ethanol extract of L. japonica to possess a strong inhibitory effect against H1N1 influenza virus-infected MDCK cells and ICR mice[94]. The incidence of the COVID-19 pandemic called to action various scientists in a bid to find safe and efficacious treatment[95]. Traditional Chinese medicines became an attractive alternative in this search. The water extract of the flower bud of L. japonica which has traditionally served as a good antipyretic and antitussive agent attracted the attention of researchers. Scientific evidences have confirmed that the water extract of L. japonica can induce let-7a expression in human rhabdomyosarcoma cells or neuronal cells and blood of lactating mice, inhibiting the entry and replication of the virus in vitro and in vivo[96]. In addition, the water extract of L. japonica also inhibits the fusion of human lung cancer cells Calu-3 expressing ACE2 receptor and BGK-21 cells transfected with SARS-CoV-2 spike protein, and up-regulates the expression of miR-148b and miR-146a[97].

    Oxidative stress has been implicated in the pathophysiology of many diseases, hence, amelioration of the same could be a good therapeutic approach[98,99]. In keeping with this therapeutic strategy, various compounds from the Lonicera genus have demonstrated the ability to relieve oxidative stress due to their pronounced antioxidant effects. For instance, the polyphenolic extract of L. caerulea berry was found to activate the expression of AMPK-PGC1α-NRF1-TFAM proteins in the skeletal muscle mitochondria, improve the activity of SOD, CAT and GSH-Px enzymes in blood and skeletal muscle, relieve exercise fatigue in mice by reducing oxidative stress in skeletal muscle, and enhance mitochondrial biosynthesis and cell proliferation[100]. The diverse health benefits of the anthocyanins from L. japonica have been mainly credited to their antioxidant and anti-inflammatory effects. The anthocyanin and cyanidin-3-o-glucoside have been reported to possess the potential to prolong life and delay senescence of Drosophila through the activation of the KEAP1/NRF2 signaling pathway[101].

    The liver is an essential organ that contributes to food digestion and detoxification of the body. These functions expose the liver to diverse toxins and metabolites. The Lonicera genus is rich in phytochemicals that confer protection on the liver against various toxins. The phenolic compound, 4, 5-di-O-Caffeoylquinic acid methyl ester was shown to be able to improve H2O2-induced liver oxidative damage in HepG2 cells by targeting the Keap1/Nrf2 pathway[102]. Hepatic fibrosis is a complex dynamic process, with the propensity to progress to liver cancer in severe cases. The L. japonicae flos water extract solution increased the cell viability of FL83B cells treated with thioacetamide (TAA), decreased the levels of serum alanine aminotransferase (ALT) and alkaline phosphatase (ALP), inhibited the transformation growth factor β1 (TGF-β1) and liver collagen deposition[103]. Sweroside, a secoiridoid glucoside isolate of L. japonica is known to protect the C57BL/6 mice liver from hepatic fibrosis by up-regulating miR-29a and inhibiting COL1 and TIMP1[104].

    Aside from the aforementioned, other pharmacological effects have been ascribed to the Lonicera genus. The ethanolic extract of L. caerulea has been reported to inhibit the proliferation of SMMC-7721 and H22 hepatoma cells, while its anthocyanins induced the apoptosis of tumor cells via the release of cytochrome C and activation of caspase[105]. AMPK/PPARα axes play an important role in lipid metabolism. A chlorogenic acid-rich extract of L. Japonica was found to significantly decrease the early onset of high-fat diet-induced diabetes in Sprague-Dawley rats via the CTRPs-AdipoRs-AMPK/PPARα axes[106]. In a high-fat diet-induced non-alcoholic fatty liver disease in C57BL/6 mice, treatment with L. caerulea polyphenol extract decreased serum inflammatory factors and endotoxin levels and the Firmicutes/Bacteroidetes ratio, an indication of its modulatory effect on the gut microbiota[107]. The iridoid-anthocyanin extract from L. caerulea berry contributed to alleviating the symptoms of intestinal infection with spirochaeta in mice[108].

    The traditional classification of the Lonicera genus based on the morphology of member plants is further categorized into two subgenera, Chamaecerasus and Periclymenum. The Chamaecerasus includes four categories, Coeloxylosteum, Isika, Isoxylosteum and Nintooa. The Periclymenum includes two categories, Subsect. Lonicera and Subsect. Phenianthi (Supplemental Table S1).

    High-throughput chloroplast genome sequencing of L. japonica found its length to be 155078 bp, which is similar to the structure of the typical angiosperm chloroplast genome. It contains a pair of inverted repeat regions (IRa and IRb, 23774 bp), a large single copy region (LSC, 88858 bp) and a small single copy area (SSC, 18672 bp)[109,110]. However, compared with chloroplast genomes of other plants, the chloroplast genome of L. japonica has a unique rearrangement between trnI-CAU and trnN-GUU[110]. Based on the phylogenetic analysis of the plastid genomes of seven plants in the Lonicera genus, 16 diverging hot spots were identified as potential molecular markers for the development of the Lonicera plants[111]. The phylogeny of Lonicera is rarely researched at the molecular level and the pattern of repetitive variation and adaptive evolution of the genome sequence is still unknown. Chloroplast genome sequences are highly conserved, but insertions and deletions, inversions, substitutions, genome rearrangements, and translocations also occur and have become powerful tools for studying plant phylogeny[112,113].

    We present here the phylogenetic tree of the Lonicera genus based on the published complete chloroplast genome sequences downloaded from the National Center for Biotechnology Information (NCBI) database using the Maximum likelihood method (Fig. 4). Based on our chloroplast phylogenies, we propose to merge L. harae into Sect. Isika and L. insularis into Chamaecerasus, but whether L. insularis belongs to Sect. Isika or Sect. Coeloxylosteum is uncertain. Based on protein-coding regions (CDS) of the chloroplast genome or complete chloroplast genomes, Liu et al. and Chen et al. supported the classification of the two subgenera in Lonicera[111,114]. Sun et al. and Srivastav et al. demonstrated a classification between the two subgenera with more species by using sequences of nuclear loci generated, chloroplast genome, and restriction site-associated DNA sequencing (RADSeq)[115,116]. However, our phylogenetic analysis and that of Sun et al. show relations within the subgenus Chamaecerasus are tanglesome in some respects[116]. Plant traits are affected by the environment to varying degrees. Since evidence of plant speciation is implicit in its genome sequence, comparative analysis at the molecular level provides a relatively accurate depiction of inherent changes that might have occurred over time. These findings suggest the need for more species of the Lonicera genus to be sequenced to provide a more accurate theoretical basis for the evolution of the Lonicera plants and a more effective revision in the classification of the Lonicera genus.

    Figure 4.  Phylogenetic tree of 42 species of the Lonicera genus based on complete chloroplast genome sequence data. The phylogenetic tree was constructed by the maximum likelihood method. Coeloxylosteum, Isika, Isoxylosteum, and Nintooa belong to Chamaecerasus and Subsect. Lonicera belongs to the Periclymenum. Chamaecerasus and Periclymenum are the two subgenera of Lonicera. 'Not retrieved' indicates that the species failed to retrieve a subordinate taxon in the Lonicera.

    The Lonicera genus is rich in diverse bioactive compounds with immeasurable prospects in many fields. Members of this genus have been used for thousands of years in traditional Chinese medicine for heat-clearing and detoxification. These plants generally have a good taste and form part of the ingredients of various fruit juices. In cosmetics, they are known to possess anti-aging and moisturizing functions. Plants of the Lonicera genus are also known for their good ecological adaptability and can be used to improve soil and ecological environment. Based on the value of the Lonicera genus, besides researching their use through molecular biological means, their efficient utilization can also be promoted in the following ways: (1) The stems and leaves of the plants could be developed for consumption and use since the chemical profiles of these parts do not differ significantly from the flowers. This way, the wastage of this scarce resource could be minimized or avoided. (2) Most of the Lonicera plants are vines or shrubs and their natural regeneration speed is slow, so the introduction and domestication of species could be strengthened to avoid overexploitation of wild resources.

    At present, only the research on the biosynthesis and efficacy of chlorogenic acid is quite comprehensive and has been used widely in various fields. There is limited research on various aspects of other bioactive compounds and should therefore be given priority in future research goals. Currently, the multi-omics analytical approach has gradually evolved as a reliable and helpful analytical platform. Hence, multi-omics research on the Lonicera genus could lead to discoveries in drug discovery and human health.

    The authors confirm contribution to the paper as follows: study conception and design, draft manuscript preparation: Yin X, Chen X, Li W, Tran LSP, Lu X; manuscript revision: Yin X, Chen X, Li W, Tran LSP, Lu X, Chen X, Yin X, Alolga RN; data/literature collection: Chen X, Yin X; figure preparation: Chen X, Yin X; figure revision: Alolga RN, Yin X, Chen X, Li W, Tran LSP, Lu X. All authors reviewed the results and approved the final version of the manuscript.

    All data generated or analyzed during this study are included in this published article and its supplementary information file.

    This work was partially supported by the National Natural Science Foundation of China (NSFC, Nos 82173918 and 82373983).

  • The authors declare that they have no conflict of interest. Xiaojian Yin is the Editorial Board member of Medicinal Plant Biology who was blinded from reviewing or making decisions on the manuscript. The article was subject to the journal's standard procedures, with peer-review handled independently of this Editorial Board member and the research groups.

  • [1]

    Antonelli A, Ferrari SM, Corrado A, Di Domenicantonio A, Fallahi P. 2015. Autoimmune thyroid disorders. Autoimmunity Reviews 14:174−80

    doi: 10.1016/j.autrev.2014.10.016

    CrossRef   Google Scholar

    [2]

    von Hafe M, Neves JS, Vale C, Borges-Canha M, Leite-Moreira A. 2019. The impact of thyroid hormone dysfunction on ischemic heart disease. Endocrine Connections 8:R76−R90

    doi: 10.1530/EC-19-0096

    CrossRef   Google Scholar

    [3]

    Delitala AP, Scuteri A, Doria C. 2020. Thyroid hormone diseases and osteoporosis. Journal of Clinical Medicine 9:1034

    doi: 10.3390/jcm9041034

    CrossRef   Google Scholar

    [4]

    Biondi B, Kahaly GJ, Robertson RP. 2019. Thyroid dysfunction and diabetes mellitus: Two Closely associated disorders. Endocrine Reviews 40:789−824

    doi: 10.1210/er.2018-00163

    CrossRef   Google Scholar

    [5]

    Chen WH, Chen YK, Lin CL, Yeh JH, Kao CH. 2015. Hashimoto's thyroiditis, risk of coronary heart disease, and L-thyroxine treatment: a nationwide cohort study. The Journal of Clinical Endocrinology and Metabolism 100:109−14

    doi: 10.1210/jc.2014-2990

    CrossRef   Google Scholar

    [6]

    Liu CL, Cheng SP, Lin HW, Lai YL. 2014. Risk of thyroid cancer in patients with thyroiditis: a population-based cohort study. Annals of Surgical Oncology 21:843−49

    doi: 10.1245/s10434-013-3363-1

    CrossRef   Google Scholar

    [7]

    Hussein O, Abdelwahab K, Hamdy O, Awny S, Megahed NA, et al. 2020. Thyroid cancer associated with Hashimoto thyroiditis: similarities and differences in an endemic area. Journal of the Egyptian National Cancer Institute 32:7

    doi: 10.1186/s43046-020-0017-9

    CrossRef   Google Scholar

    [8]

    Vukovic R, Zeljkovic A, Bufan B, Spasojevic-Kalimanovska V, Milenkovic T, et al. 2019. Hashimoto thyroiditis and dyslipidemia in childhood: A review. Frontiers in Endocrinology 10:868

    doi: 10.3389/fendo.2019.00868

    CrossRef   Google Scholar

    [9]

    Verstuyf A, Carmeliet G, Bouillon R, Mathieu C. 2010. Vitamin D: a pleiotropic hormone. Kidney International 78:140−45

    doi: 10.1038/ki.2010.17

    CrossRef   Google Scholar

    [10]

    Krysiak R, Szkróbka W, Okopień B. 2019. The relationship between statin action on thyroid autoimmunity and vitamin D status: A pilot study. Experimental and Clinical Endocrinology & Diabetes 127:23−28

    doi: 10.1055/a-0669-9309

    CrossRef   Google Scholar

    [11]

    Nodehi M, Ajami A, Izad M, Asgarian Omran H, Chahardoli R, et al. 2019. Effects of vitamin D supplements on frequency of CD4+ T-cell subsets in women with Hashimoto's thyroiditis: a double-blind placebo-controlled study. European Journal of Clinical Nutrition 73:1236−43

    doi: 10.1038/s41430-019-0395-z

    CrossRef   Google Scholar

    [12]

    Jamka M, Ruchala M, Walkowiak J. 2019. Vitamin D and Hashimoto's disease. Polski merkuriusz lekarski 47:111−13

    Google Scholar

    [13]

    Effraimidis G, Badenhoop K, Tijssen JGP, Wiersinga WM. 2012. Vitamin D deficiency is not associated with early stages of thyroid autoimmunity. European Journal of Endocrinology 167:43−48

    doi: 10.1530/EJE-12-0048

    CrossRef   Google Scholar

    [14]

    Moher D, Liberati A, Tetzlaff J, Altman DG. 2009. Preferred reporting items for systematic reviews and meta-analyses: the PRISMA statement. BMJ 339:b2535

    doi: 10.1136/bmj.b2535

    CrossRef   Google Scholar

    [15]

    Meads C, Sutton A, Małysiak S, Kowalska M, Zapalska A, et al. 2013. Sentinel lymph node status in vulval cancer: systematic reviews of test accuracy and decision-analytic model-based economic evaluation. Health Technology Assessment 17:1−216

    doi: 10.3310/hta17600

    CrossRef   Google Scholar

    [16]

    Botelho IMB, Moura Neto A, Silva CA, Tambascia MA, Alegre SM, et al. 2018. Vitamin D in Hashimoto's thyroiditis and its relationship with thyroid function and inflammatory status. Endocrine Journal 65:1029−37

    doi: 10.1507/endocrj.EJ18-0166

    CrossRef   Google Scholar

    [17]

    Bozkurt NC, Karbek B, Ucan B, Sahin M, Cakal E, et al. 2013. The association between severity of vitamin D deficiency and Hashimoto's thyroiditis. Endocrine Practice 19:479−84

    doi: 10.4158/EP12376.OR

    CrossRef   Google Scholar

    [18]

    De Pergola G, Triggiani V, Bartolomeo N, Giagulli VA, Anelli M, et al. 2018. Low 25 hydroxyvitamin D levels are independently associated with autoimmune thyroiditis in a cohort of apparently healthy overweight and obese subjects. Endocrine, Metabolic & Immune Disorders Drug Targets 18:646−52

    doi: 10.2174/1871530318666180406163426

    CrossRef   Google Scholar

    [19]

    Ke W, Sun T, Zhang Y, He L, Wu Q, et al. 2017. 25-Hydroxyvitamin D serum level in Hashimoto's thyroiditis, but not Graves' disease is relatively deficient. Endocrine Journal 64:581−87

    doi: 10.1507/endocrj.EJ16-0547

    CrossRef   Google Scholar

    [20]

    Kim D. 2016. Low vitamin D status is associated with hypothyroid Hashimoto's thyroiditis. Hormones 15:385−93

    doi: 10.14310/horm.2002.1681

    CrossRef   Google Scholar

    [21]

    Kivity S, Agmon-Levin N, Zisappl M, Shapira Y, Nagy EV, et al. 2011. Vitamin D and autoimmune thyroid diseases. Cellular & Molecular Immunology 8:243−47

    doi: 10.1038/cmi.2010.73

    CrossRef   Google Scholar

    [22]

    Ma J, Wu D, Li C, Fan C, Chao N, et al. 2015. Lower serum 25-hydroxyvitamin D level is associated with 3 types of autoimmune thyroid diseases. Medicine 94:e1639

    doi: 10.1097/MD.0000000000001639

    CrossRef   Google Scholar

    [23]

    Maciejewski A, Wójcicka M, Roszak M, Losy J, Łącka K. 2015. Assessment of Vitamin D Level in Autoimmune Thyroiditis Patients and a Control Group in the Polish Population. Advances in Clinical and Experimental Medicine 24:801−6

    doi: 10.17219/acem/29183

    CrossRef   Google Scholar

    [24]

    Mansournia N, Mansournia MA, Saeedi S, Dehghan J. 2014. The association between serum 25OHD levels and hypothyroid Hashimoto’s thyroiditis. Journal of Endocrinological Investigation 37:473−76

    doi: 10.1007/s40618-014-0064-y

    CrossRef   Google Scholar

    [25]

    Sönmezgöz E, Ozer S, Yilmaz R, Önder Y, Bütün I, et al. 2016. Hypovitaminosis D in children with Hashimoto's thyroiditis. Revista Medica De Chile 144:611−16

    doi: 10.4067/S0034-98872016000500009

    CrossRef   Google Scholar

    [26]

    Tamer G, Arik S, Tamer I, Coksert D. 2011. Relative Vitamin D Insufficiency in Hashimoto's Thyroiditis. Thyroid 21:891−96

    doi: 10.1089/thy.2009.0200

    CrossRef   Google Scholar

    [27]

    Ucan B, Sahin M, Sayki Arslan M, Colak Bozkurt N, Kizilgul M, et al. 2016. Vitamin D treatment in patients with Hashimoto's thyroiditis may decrease the development of hypothyroidism. International Journal for Vitamin and Nutrition Research 86:9−17

    doi: 10.1024/0300-9831/a000269

    CrossRef   Google Scholar

    [28]

    Xu J, Zhu X, Sun H, Xu X, Xu S, et al. 2018. Low vitamin D levels are associated with cognitive impairment in patients with Hashimoto thyroiditis. BMC Endocrine Disorders 18:87

    doi: 10.1186/s12902-018-0314-7

    CrossRef   Google Scholar

    [29]

    Yasmeh J, Farpour F, Rizzo V, Kheradnam S, Sachmechi I. 2016. Hashimoto thyroiditis not associated with vitamin D deficiency. Endocrine Practice 22:809−13

    doi: 10.4158/EP15934.OR

    CrossRef   Google Scholar

    [30]

    Unal AD, Tarcin O, Parildar H, Cigerli O, Eroglu H, et al. 2014. Vitamin D deficiency is related to thyroid antibodies in autoimmune thyroiditis. Central-European Journal of Immunology 39:493−97

    doi: 10.5114/ceji.2014.47735

    CrossRef   Google Scholar

    [31]

    Evliyaoğlu O, Acar M, Özcabı B, Erginöz E, Bucak F, et al. 2015. Vitamin D deficiency and Hashimoto's thyroiditis in children and adolescents: a critical vitamin D level for this association? Journal of Clinical Research in Pediatric Endocrinology 7:128−33

    doi: 10.4274/jcrpe.2011

    CrossRef   Google Scholar

    [32]

    Chaudhary N, Kumar R, Sachdeva N, Dayal D. 2018. Vitamin D levels in children with Hashimoto's thyroiditis: Before and after L-thyroxine therapy. Thyroid Research and Practice 15:23

    doi: 10.4103/trp.trp_45_17

    CrossRef   Google Scholar

    [33]

    Camurdan OM, Döğer E, Bideci A, Celik N, Cinaz P. 2012. Vitamin D status in children with Hashimoto thyroiditis. Journal of Pediatric Endocrinology & Metabolism 25:467−70

    Google Scholar

    [34]

    Dellal FD, Niyazoglu M, Ademoglu E, Gorar S, Candan Z, et al. 2013. Evaluation of serum trace elements and vitamin levels in Hashimoto’s thyroiditis: single centre experience from Turkey. Open Journal of Endocrine and Metabolic Diseases 3:236−40

    doi: 10.4236/ojemd.2013.34031

    CrossRef   Google Scholar

    [35]

    Şıklar Z, Karataş D, Doğu F, Hacıhamdioğlu B, İkincioğulları A, et al. 2016. Regulatory T cells and vitamin D status in children with chronic autoimmune thyroiditis. Journal of Clinical Research in Pediatric Endocrinology 8:276−81

    doi: 10.4274/jcrpe.2766

    CrossRef   Google Scholar

    [36]

    Nalbant A, Gökosmanoğlu F, Cinemre H, Varım C, Kaya T, et al. 2017. The relation between serum vitamin D levels and Hashimoto thyroiditis in women. The Journal of the Kuwait Medical Association 49:223−26

    Google Scholar

    [37]

    Giovinazzo S, Vicchio TM, Certo R, Alibrandi A, Palmieri O, et al. 2017. Vitamin D receptor gene polymorphisms/haplotypes and serum 25(OH)D3 levels in Hashimoto's thyroiditis. Endocrine 55:599−606

    doi: 10.1007/s12020-016-0942-5

    CrossRef   Google Scholar

    [38]

    Guleryuz B, Akin F, Tunc Ata M, Mergen Dalyanoglu M, Turgut S. 2016. Vitamin-D Receptor (VDR) gene polymorphisms (TaqI, FokI) in Turkish patients with Hashimoto's thyroiditis: relationship to the levels of Vit-D and cytokines. Endocrine, Metabolic & Immune Disorders - Drug Targets 16:131−39

    doi: 10.2174/1871530316666160728092613

    CrossRef   Google Scholar

    [39]

    Perga S, Martire S, Montarolo F, Giordani I, Spadaro M, et al. 2018. The footprints of poly-autoimmunity: evidence for common biological factors involved in multiple sclerosis and Hashimoto's thyroiditis. Frontiers in Immunology 9:311

    doi: 10.3389/fimmu.2018.00311

    CrossRef   Google Scholar

    [40]

    Rezaee H, Najafipour F, Ranjdoust F, Sadra V, Hamishekar H, et al. 2017. Evaluation of serum vitamin D levels in patients with Hashimoto's thyroiditis. International Journal of Advanced Biotechnology and Research 2:890−95

    doi: 10.4274/turkderm.galenos.2020.66664

    CrossRef   Google Scholar

    [41]

    Prasad I, Kumari R, A S. 2016. Vitamin D evaluation in autoimmune thyroid diseases. International Journal of Contemporary Medical Research 12:3415−18

    Google Scholar

    [42]

    Chao G, Zhu Y, Fang L. 2020. Correlation between Hashimoto's thyroiditis-related thyroid hormone levels and 25-hydroxyvitamin D. Frontiers in Endocrinology 11:4

    doi: 10.3389/fendo.2020.00004

    CrossRef   Google Scholar

    [43]

    Ahi S, Dehdar MR, Hatami N. 2020. Vitamin D deficiency in non-autoimmune hypothyroidism: a case-control study. BMC Endocrine Disorders 20:41

    doi: 10.1186/s12902-020-0522-9

    CrossRef   Google Scholar

    [44]

    Feng Y, Qiu T, Chen H, Wei Y, Jiang X. et al. 2020. Association of serum IL-21 and vitamin D concentrations in Chinese children with autoimmune thyroid disease. Clin Chim Acta 507:194−98

    Google Scholar

    [45]

    Chen P, Chen G, Xu S, Li W, Liu C. 2015. Detection and significance of serum 25-hydroxyvitamin D in patients with Hashimoto thyroiditis. Jiangsu Medical Journal 41:675−77

    Google Scholar

    [46]

    Liu X, Zhang H. 2012. Correlation between serum level of 25-(OH)D3 and humoral immunity in patients with autoimmune thyroid diseases. The Journal of Practical Medicine 28:1455−57

    Google Scholar

    [47]

    Xiang Q, Zhang X, Li W, Guo Y, Bai Y, et al. 2017. Research on the relationship between 25-hydroxyvitamin D level and hashimoto's thyroiditis. Modern Journal of Integrated Traditional Chinese and Western Medicine 26:3877−79+95

    Google Scholar

    [48]

    Zhang W, An M. 2015. The association of serum 25(OH)D levels with thyroid hormones and autoantibodies in AITD patients. Zhejiang Clinical Medical Journal 2015:1067−68,69

    Google Scholar

    [49]

    Li N, Jiang Y, Yan S, Fang F, Gu L, et al. 2015. Correlation analysis between serum vitamin D and Hashimoto's thyroiditis. World Clinical Drugs 36:591−95

    Google Scholar

    [50]

    Cvek M, Kaličanin D, Barić A, Vuletić M, Gunjača I et al. 2021. Vitamin D and Hashimoto’s thyroiditis: observations from CROHT biobank. Nutrients 13:2793

    doi: 10.3390/nu13082793

    CrossRef   Google Scholar

    [51]

    Salem TM, Abdelmonem E, Fayad A. 2021. Hashimoto's thyroiditis, iron, and vitamin D deficiency among Egyptian female patients: associations and possible causalities. Hormones 20:833−36

    doi: 10.1007/s42000-021-00297-z

    CrossRef   Google Scholar

    [52]

    Hanna HWZ, Rizzo C, Abdel Halim RM, El Haddad HE, Salam R, et al. 2021. Vitamin D status in Hashimoto's thyroiditis and its association with vitamin D receptor genetic variants. The Journal of Steroid Biochemistry and Molecular Biology 212:105922

    doi: 10.1016/j.jsbmb.2021.105922

    CrossRef   Google Scholar

    [53]

    Głowinska-Olszewska B, Borysewicz-Sańczyk H, Sawicka B, Klonowska B, Charemska D, et al. 2020. Does Hashimoto's Thyroiditis increase the risk of cardiovascular disease in young type 1 diabetic patients? Frontiers in Endocrinology 11:431

    doi: 10.3389/fendo.2020.00431

    CrossRef   Google Scholar

    [54]

    Ahi S, Dehdar MR, Hatami N. 2020. Vitamin D deficiency in non-autoimmune hypothyroidism: a case-control study. BMC Endocrine Disorders 20:41

    doi: 10.1186/s12902-020-0522-9

    CrossRef   Google Scholar

    [55]

    Ren X, Wu H, Zhang T. 2021. Study on the changes of serum 25(OH)D concentration in patients with Type 2 Diabetes complicated with Hashimoto's thyroiditis and its correlation with clinical indicators. Journal of Medical Research 50:71−74

    Google Scholar

    [56]

    Huang S, Xue X, Li Y. 2018. Relevance of 25-hydroxyvitamin D and autoimmune thyroid diseases in children and adolescents. Modern Hospital 18:1792−94

    Google Scholar

    [57]

    Chi M, Liu J. 2020. Relationship analysis between 25-hydroxyvitamin D3 level and autoimmune thyroid disease in children and adolescents. Smart Healthcare 6:86−87

    Google Scholar

    [58]

    Yang Y, Qiu W, Zou H. 2021. Expression and correlation of CD4+ T lymphocytes and 25(OH)D in peripheral blood in patients with Hashimoto's thyroiditis. Chinese Journal of Practical Medicine 48:63−66

    Google Scholar

    [59]

    Ke W, Gu Y, Liu J, Jin X. 2021. The correlation between serum IL-17 and IL-35 and vitamin D in Hashimoto's thyroiditis. Lab Med 36:181−84

    Google Scholar

    [60]

    Wu H. 2019. Effects and significance of vitamin D on Hashimoto thyroiditis. Journal of New Medicine 29:279−280+283

    Google Scholar

    [61]

    Fu J, Li X, Zhang H, Xu Y, Zhu Y. 2021. Study on the relationship between vitamin D deficiency and Hashimoto's thyroiditis. Journal of Internal Medicine Concepts & Practice 16:27−31

    Google Scholar

    [62]

    Pearce EN, Farwell AP, Braverman LE. 2003. Thyroiditis. New England Journal of Medicine 348:2646−55

    doi: 10.1056/NEJMra021194

    CrossRef   Google Scholar

    [63]

    Dankers W, Colin EM, van Hamburg JP, Lubberts E. 2017. Vitamin D in autoimmunity: molecular mechanisms and therapeutic potential. Frontiers in Immunology 7:697

    doi: 10.3389/fimmu.2016.00697

    CrossRef   Google Scholar

    [64]

    Wang J, Lv S, Chen G, Gao C, He J, et al. 2015. Meta-analysis of the association between vitamin D and autoimmune thyroid disease. Nutrients 7:2485−98

    doi: 10.3390/nu7042485

    CrossRef   Google Scholar

    [65]

    Štefanić M, Tokić S. 2020. Serum 25-hydoxyvitamin D concentrations in relation to Hashimoto's thyroiditis: a systematic review, meta-analysis and meta-regression of observational studies. European Journal of Nutrition 59:859−72

    doi: 10.1007/s00394-019-01991-w

    CrossRef   Google Scholar

    [66]

    Boonstra A, Barrat FJ, Crain C, Heath VL, Savelkoul HF, et al. 2001. 1α,25-Dihydroxyvitamin d3 has a direct effect on naive CD4+ T cells to enhance the development of Th2 cells. Journal of Immunology 167:4974−80

    doi: 10.4049/jimmunol.167.9.4974

    CrossRef   Google Scholar

    [67]

    Baeke F, Takiishi T, Korf H, Gysemans C, Mathieu C. 2010. Vitamin D: modulator of the immune system. Current Opinion in Pharmacology 10:482−96

    doi: 10.1016/j.coph.2010.04.001

    CrossRef   Google Scholar

    [68]

    Fisher SA, Rahimzadeh M, Brierley C, Gration B, Doree C, et al. 2019. The role of vitamin D in increasing circulating T regulatory cell numbers and modulating T regulatory cell phenotypes in patients with inflammatory disease or in healthy volunteers: A systematic review. PLoS One 14:e0222313

    doi: 10.1371/journal.pone.0222313

    CrossRef   Google Scholar

    [69]

    Chahardoli R, Saboor-Yaraghi AA, Amouzegar A, Khalili D, Vakili AZ, et al. 2019. Can supplementation with vitamin D modify thyroid autoantibodies (Anti-TPO Ab, Anti-Tg Ab) and Thyroid profile (T3, T4, TSH) in Hashimoto's thyroiditis? A double blind, randomized clinical trial HormMetab Res 51:296−301

    doi: 10.1055/a-0856-1044

    CrossRef   Google Scholar

    [70]

    Vahabi AP, Aminorroaya A, Amini M, Momeni F, Feizi A, et al. 2017. Effect of Vitamin D deficiency treatment on thyroid function and autoimmunity markers in Hashimoto's thyroiditis: A double-blind randomized placebo-controlled clinical trial. Journal of Research in Medical Sciences 22:103

    doi: 10.4103/jrms.JRMS_1048_16

    CrossRef   Google Scholar

    [71]

    Kmiec P, Minkiewicz I, Rola R, Sworczak K, Zmijewski MA, et al. 2018. Vitamin D status including 3-epi-25(OH)D3 among adult patients with thyroid disorders during summer months. Endokrynologia Polska 69:653−60

    doi: 10.5603/EP.a2018.0065

    CrossRef   Google Scholar

  • Cite this article

    Liu Z, Feng L, He Y, Yuan S, Xu C. 2022. The Association between Vitamin D and Hashimoto Thyroiditis: An Up-to-date Systematic Review and Meta-analysis. Food Materials Research 2:9 doi: 10.48130/FMR-2022-0009
    Liu Z, Feng L, He Y, Yuan S, Xu C. 2022. The Association between Vitamin D and Hashimoto Thyroiditis: An Up-to-date Systematic Review and Meta-analysis. Food Materials Research 2:9 doi: 10.48130/FMR-2022-0009

Figures(8)  /  Tables(1)

Article Metrics

Article views(5624) PDF downloads(1068)

Other Articles By Authors

ARTICLE   Open Access    

The Association between Vitamin D and Hashimoto Thyroiditis: An Up-to-date Systematic Review and Meta-analysis

Food Materials Research  2 Article number: 9  (2022)  |  Cite this article

Abstract: The objective of the surrent study was to summarize the up-to-date studies to investigate the relationship between vitamin D and Hashimoto thyroiditis (HT). An online search of English and Chinese databases was performed. The studies concerned the investigation of the relationship between vitamin D and HT including meta-analysis, meanwhile the heterogeneities were revealed by subgroup analysis. Fourty six elated studies containing 15,336 participants (HT: 6,138 versus control: 9,198) were included. HT patients had lower levels of 25(OH)D3 (standardised mean difference, −1.09; 95%CI: [−1.42, −0.75]; P < 0.01), and were more likely to be deficient in 25(OH)D3 (OR, 2.77; 95%CI, [1.88, 3.91]; P < 0.05). Obvious heterogeneities in the results of meta-analysis were down to the difference of detection methods and criteria of vitamin D insufficiency among studies. Vitamin D deficiency was colncluded to have a significant relation with HT.

    • T and B lymphocyte mediated immune tolerance disorder led to abnormal auto-antibodies invasion in thyroid glands, which was the main mechanism of Hashimoto thyroiditis (HT). The increase of HT continued[1], and gradually became the most common cause of thyroid hormone insufficiency. Besides, it was suggested a link to diseases including ischemic heart disease, osteoporosis, diabetes, cardiovascular disease and cancer[26]. Papillary thyroid carcinoma was the most common form of cancer associated with HT[7]. Even the health of children was threatened by HT with higher risk of dyslipidemia and cardiovascular disease[8]. Along with the progress of HT, the secretion of thyroid hormone was significantly insufficient, and patients had to replenish with Euthyroxin throughout their life, but still suffered from higher risk of HT related complications. Therefore, effective methods for HT prevention were considered to be urgent research.

      Since Vitamin D receptors were shown to be present in the thyroid[9], Vitamin D was considered as HT prevention in resent studies. Krysiak et al.[10] suggested that Vitamin D combined with atorvastatin could improve thyroid autoimmunity. In addition, supplementation with cholecalciferol for HT patients was indicated to twist the balance of CD4+ T-cell subsets toward ameliorative composition[11]. Although the relationship between Vitamin D and HT seemed clear, it still needed further investigation[12], for there existed study reported they had no evident correlation[13]. These contradictory conclusions maybe due to the small sample size of the local population, but it was very difficult to conduct a high-quality epidemiological study with a large sample size. Therefore, to indicate whether vitamin D deficiency was really correlated to HT, which may affect clinical strategy whether we should use Vitamin D to prevent HT, evidence-based medical research tools such as systematic review and meta-analysis were necessary for a more reliable conclusion.

      We conducted a systematic review and meta-analysis of the domestic and foreign studies that investigated the relationship between Vitamin D and HT. Compatible data were pooled into meta-analysis to provide rigorous evidence-based medical reference for the prevention of HT.

    • The present systematic review and meta-analysis were performed under PRISMA guidelines[14]. An online bibliographic search was performed in Pubmed (for English), and China National Knowledge Infrastructure (CNKI) and Wanfang Databases (both for Chinese) by two investigators with key words of 'Vitamin D' and 'Hashimoto thyroiditis'. Studies were considered if they were in English or Chinese, and updated to 27 September 2021. If the abstract displayed investigation about the relationship between Vitamin D and HT, the full text was read in detail by at least two authors (Fig. 1).

      Figure 1. 

      Forest plot of 25(OH)D3 level (random model).

    • Studies were finally included if they fulfilled the following criteria: (1) included the comparison of HT patient group with a healthy control group; (2) serum levels of 25(OH)D3 level and/or the quantity of patients with 25(OH)D3 insufficiency were reported; (3) written in English or Chinese; (4) with a quality score above or equal to 6 according to the coding manual for case-control studies[15] assessed by two authors respectively.

    • The following information was extracted from each study by two investigators independently: (1) the first author; (2) year of publication; (3) region; (4) sample size; (5) serum levels of 25(OH)D3; (6) cut-off of serum 25(OH)D3 insufficiency; (7) the quantity of patients with 25(OH)D3 insufficiency; and (8) quality score. After that, the extracted information was summarized and checked by another two authors.

    • RevMan 5.3 (the Cochrane Collaboration) was used to perform a meta-analysis on the data obtained. Firstly, Weighted mean differences (WMD) for continuous variable and Odds Ratios (OR) for binary variables were calculated. Subsequently, statistical heterogeneity was assessed with I2 test, and the main source of heterogeneity was revealed by subgroup analysis. Lastly, publication bias was evaluated by funnel plots. For statistical analysis above, 'P < 0.05' was considered a significant difference between groups.

    • Online search obtained 336 studies, in which 280 were excluded by abstract screening. Then, of the 56 remaining, 10 were excluded by full text in-detail evaluation; finally 46 studies with 15,336 individuals in total (6,138 HT patients and 9,198 healthy controls) were included into the present study for the systematic review[1661]. Characteristics of included studies are summarized in Table 1.

      Table 1.  Characteristics of included studies.

      Study
      (published year)
      RegionSample size (HT:C)25(OH)D3
      Assay method
      Serum 25(OH)D3 level
      (HT vs C)
      (ng/mL)
      Serum 25(OH)D3 insufficiency cut off (ng/mL)Number of 25(OH)D3
      insufficiency (HT:C)
      Quality
      score
      Maciejewski et al. 2015[23]Poland62/32ELISA
      8.00 ± 5.06 vs
      12.12 ± 7.80
      < 3061/277
      Ucan et al. 2016[27]Turkey
      75/43RIA
      9.37 ± 0.69 vs
      11.9 ± 1.01
      < 2075/369
      Bozkurt et al. 2013[12-17]]Turkey
      360/180CLS12.2 ± 5.6 vs
      15.4 ± 6.8
      < 10150/378
      Kim 2016[20]Korea221/555CLS36.84 ± 22.96 vs
      39.84 ± 21.48
      < 30108/2068
      Sonmezgoz et al. 2016[25]Turkey
      68/68CLS16.8 ± 9.2 vs
      24.1 ± 9.4
      < 3061/548
      De Pergola et al. 2018[18]Italy
      45/216CLS< 2031/1138
      Botelho et al. 2018[16]Brazil
      88/71CLS26.4 (7.6–48.2) vs
      28.6 (13–51.2)
      < 3061/397
      Ma et al. 2015[22]China70/70ELISA
      12.40 ± 4.46 vs
      16.53 ± 5.79
      < 3070/677
      Yasmeh et al. 2016[29]America97/88CLS24.5 ± 6.42 vs
      20.6 ± 6.5
      < 3066/747
      Xu et al. 2018[28]China194/200CPBA16.16 (13.72–18.76) vs
      23.32 (20.84–25.92)
      7
      Kivity et al. 2011[21]Israel
      28/98CLS< 1022/308
      Mansournia et al. 2014[24]Iran
      41/45SC15.9 ± 1.21 vs
      24.4 ± 1.73
      < 2034/248
      Tamer et al. 2011[26]Turkey
      161/162RIA
      16.3 ± 10.4 vs
      29.6 ± 2.55
      < 30148/1028
      Chaudhary et al. 2018[32]India
      35/50HPLC13.39 ± 6.8 vs
      26.16 ± 12.28
      < 2031/388
      Evliyaoğlu et al. 2015[31]Turkey
      90/79HPLC16.67 ± 11.65 vs
      20.99 ± 9.86
      < 2080/698
      Unal et al. 2014[30]Turkey
      254/124CLS17.05 (5.4−80) vs
      19.9 (9−122.7)
      < 20160/-7
      Ke et al. 2017[19]China
      61/51EBL22.10 ± 1.52 vs
      33.40 ± 1.56
      < 2034/127
      Camurdan et al. 2012[33]Turkey
      78/74HPLC31.2 ± 11.5 vs
      57.9 ± 19.7
      < 2069/247
      Dellal et al 2013[34]Turkey51/27RIA
      17.3 ± 8.0 vs
      21.8 ± 15.2
      6
      Siklar et al. 2016[35]Turkey32/24HPLC16.02 ± 9.84 vs
      21.91 ± 7.68
      < 2022/107
      Nalbant et al. 2017[36]Turkey253/200CLS33 ± 29.6 vs
      43.7 ± 26.2
      < 20161/1118
      Giovinazzo et al. 2017[37]Italy
      100/100HPLC21.2 ± 12.9 vs
      35.7 ± 16.7
      < 2070/187
      Guleryuz et al. 2016[38]Turkey
      136/50HPLC14.88 ± 8.23 vs
      15.52 ± 1.34
      6
      Perga et al. 2018[39]Italy
      55/59CLS< 2037/42
      Yavuzer et al. 2017Turkey
      49/34ELISA19.5 ± 15 vs
      23.8 ± 19
      6
      Priya et al. 2016India25/27ELISA14.3 (12.65−17.90)
      vs 26.2 (21.00−32.8)
      6
      Chao et al. 2020[42]China373/4889RIA
      16.66 ± 6.51 vs
      15.81 ± 6.42
      < 20363/47389
      Feng et al. 2020[44]China36/30ELISA17.39 ± 8.49 vs
      35.15 ± 14.16
      6
      Ahi et al. 2020[43]Iran633/200CLS13.22 (8.1−24.27) vs
      20.4 (11.2−29.6)
      7
      Liu and Zhang. 2012[46]China30/20RIA
      16.48 ± 6.25 vs
      24.31 ± 7.88
      7
      Xiang et al. 2017[47]China41/106CLS19.71 ± 8.43 vs
      20.56 ± 11.64
      < 3038/906
      Zhang et al. 2015[48]China31/19HPLC17 ± 6 vs
      24 ± 7
      6
      Chen et al. 2015[45]China34/52CLS14.4 ± 5.6 vs
      17.4 ± 5.6
      < 2029/377
      Li et al. 2015[49]China50/5621.19 (18.40−25.28) vs
      24.06 (18.94−33.90)
      < 3044/376
      Cvek et al. 2021[50]Croatian461/176CLS19.7 (14.4−25.2) vs
      17.3 (13.2−22.7)
      < 20127/657
      Salem et al. 2021[51]Egypt120/120ELISA7.6 ± 4.4 vs 20.6 ± 5.5< 10120/1127
      Hana et al. 2021[52]Egypt112/48HPLC10.1 (8.7−11.7) vs 12.0 (9.3−15.6)< 30101/406
      Olszewska et al. 2020[53]Italy30/2017.9 ± 7.9 vs 18.5 ± 8.16
      Rezaee et al. 2017[40]Iran51/45CLS6
      Ren et al. 2021[55]China62/8013.49 ± 4.32 vs 15.75 ± 5.85< 3060/766
      Huang et al. 2018[56]China61/50CLS16.27 ± 6.99 vs 29.01 ± 9.72< 206
      Chi et al. 2020[57]China32/30CLS15.27 ± 5.98 vs 28.89 ± 9.586
      Yang et al. 2021[58]China88/6013.37 ± 3.49 vs 17.58 ± 5.636
      Ke et al. 2021[59]China152/50CLS20.56 ± 1.4 vs 33.4 ± 6.5< 2090/67
      Wang et al. 2015[64]China31/30ELISA10.08 ± 0.44 vs 14.32 ± 3.746
      Fu et al. 2021[61]China334/30016.84 (11.81, 23.39) vs 16.66 (11.98, 22.13)< 30214/2097
      H: hashimoto thyroiditis group; C: Healthy control group; ELISA: Enzyme Linked Immunosorbent Assay; RIA: Radioimmunoassay; CLS: Chemiluminesent lmmunoassay Assay; CPBA: competitive protein binding assay; SC: Solid Chromatography, HPLC: High Performance Liquid Chromatography, EBL: Euglobulin lysis method, −: Non reported.
    • Meta-analysis included 33 studys with 3,161 patients in HT group and 7,488 healthy individuals in the control group for comparison. Random model indicated 25(OH)D3 levels of HT group were significantly lower than the control group (WMD: −7.44, 95%CI [−9.29, −5.60], P < 0.01). I2 test (98%) suggested significant heterogeneity in the meta-analysis (Fig. 1). The subgroup meta-analysis basing on 25(OH)D3 assays in a fixed model revealed similar results (WMDs: −0.55; 95%CI [−0.60, −0.49], P < 0.01), and its significant heterogeneity among subgroups represented by I2 = 98.3% suggested the difference of 25(OH)D3 assays was the main source of heterogeneity (Fig. 2,). Lastly, we separated the Chinese studies with 6,639 individuals (HT: 1,102 vs C: 5537) to perform another particle meta-analysis in random model. Result showed that, in the Chinese population, serum 25(OH)D3 level of HT patients was significantly lower than that of healthy individuals (WMD: −7.04, 95%CI [−10.37, −3.71], P < 0.01 ). Meanwhile, I2 = 98.0% also suggested a significant heterogeneity (Fig. 3).

      Figure 2. 

      Subgroup forest plot of 25(OH)D3 level(fixed model).

      Figure 3. 

      Forest plot of prevalence of Vitamin D insufficiency (random model).

    • A total of 29 studies comprising 11,795 individuals (HT: 3,709 vs C: 8,086) were pooled for OR of 25(OH)D3 insufficiency. Random model indicated HT patients had higher prevalence of Vitamin D insufficiency compared to healthy individuals (OR: 2.54, 95%CI [1.77, 3.63], P < 0.01). I2 test (86%) suggested significant heterogeneity in meta-analysis (Fig. 4). Subgroup meta-analysis in a fixed model based on different 25(OH)D3 insufficiency cut-off also revealed similar results as above (OR: 1.84; 95%CI [1.64, 2.07], P < 0.01). Meanwhile, I2 equaled to 93% suggested the main source of heterogeneity was from the different cut-off of 25(OH)D3 insufficiency (Fig. 5). Chinese studies with 1,373 individuals (HT: 700 vs C: 673) were separated to perform another particle meta-analysis in random model. Results displayed a trend that the HT population had a higher prevalence of 25(OH)D3 insufficiency compared to healthy individuals, but it was not statistically significant (P > 0.05), and meanwhile significant heterogeneity was indicated by I2 equal to 91% (Fig. 6).

      Figure 4. 

      Subgroup forest plot of prevalence of Vitamin D insufficiency (fixed model).

      Figure 5. 

      Forest plot of 25(OH)D3 level (random model, Chinese studies).

      Figure 6. 

      Forest plot of prevalence of Vitamin D insufficiency (random model, Chinese studies).

    • A funnel plot of serum 25(OH)D3 level in subgroup analysis exhibited that the included studies accumulated at the top of the funnel, which suggested that publication bias may exert little adverse effect on the confidence in the meta-analysis (Fig. 7). Similarly, results of the funnel plot suggested low risk of publication bias in prevalence of 25(OH)D3 insufficiency comparisons (Fig. 8).

      Figure 7. 

      Funnel plot of 25(OH)D3 level.

      Figure 8. 

      Funnel plot of quantity of individuals with 25(OH)D3 insufficiency.

    • The present study reinforced the close relationship between Vitamin D insufficiency and HT with methods of systematic review and meta-analysis. To our knowledge, the present systematic review summarized the most related studies to date; among them, Chinese studies, which may be ignored by other foreign researchers, were also included. Hence, we believe our conclusion produce more confident evidence for a relationship between Vitamin D insufficiency and HT.

      Although a series of related factors of HT have been revealed, the real etiology has so far not been clearly understood[62]. Vitamin D has been proved to closely relate to HT, for it plays a vital role in regulating inflammatory response and maintaining immune balance[63]. Multiple epidemiological studies suggested a close relationship between Vitamin D and HT; however, differences in quality, region and population may affect the conclusions. Therefore, high quality systematic review or meta-analysis was still needed to acquire more reliable evidence. Wang et al.[64] published a meta-analysis in 2015 to indicate the relationship between Vitamin D insufficiency and HT. Štefanić et al.[65] subsequently included more recent studies for meta-analysis and drew a similar conclusion, but their results seemed to weaken the relationship of Vitamin D insufficiency and HT compared with Wang et al.[64]. However, these two studied did not include enough recent Chinese studies which should not be ignored. This may not only decrease the confidence of the conclusions, but also weaken the reliability for Chinese researchers. To fulfill this deficiency, we included Chinese studies into the present systematic review and meta-analysis. As anticipated, the general results including serum 25(OH)D3 level and quantity of individuals with Vitamin D insufficiency displayed similar results to Wang et al.[64] and Štefanić et al.[65], which meant the relationship of Vitamin D insufficiency and HT also existed in the Chinese population. We next separated the Chinese studies to perform a particle meta-analysis. The particle result of serum 25(OH)D3 levels of Chinese HT patients were significantly lower than healthy individuals generally, but its difference was less (−7.05 vs −7.44). However, concerning the prevalence of Vitamin D insufficiency, Chinese HT patients were not likely to have more Vitamin D insufficiency cases compared to healthy individuals, suggesting the relationship between Vitamin D insufficiency and HT in the Chinese population may not be as strong as in the global population. Note that, the Chinese population in the present study was only a small part of the total, the negative result maybe due to the small sample size. Further studies with larger sample sizes and high quality investigating the prevalence of Vitamin D insufficiency in the HT population are necessary in China in the future.

      With regards to the studies included in the systematic review, most studies reported lower serum 25(OH)D3 levels and higher prevalence of Vitamin D insufficiency in HT patients compared to healthy individuals. The present study had drawn a similar conclusion, but the results contained significant heterogeneity. According to the systematic review, this heterogeneity may be due to the difference of 25(OH)D3 assays, thus we performed an analysis which separated studies with the same assay into several sub-groups. Results similarly indicated lower serum 25(OH)D3 levels in HT patients, and the most significant heterogeneity among sub-groups (I2 = 99.5%), which hinted that the heterogeneity was mainly caused by the difference in 25(OH)D3 assays. In parallel with lower serum 25(OH)D3 levels, HT patients were at higher risk of 25(OH)D3 insufficiency, indicated by our meta-analysis. Meanwhile, significant heterogeneity was indicated owing to the difference of serum 25(OH)D3 insufficiency criteria. With regard to the publication bias, we determined that heterogeneity would bring significant publication bias displayed by the funnel plot, but the results showed that the included studies accumulating at the top of the funnel; this suggested the publication bias exerted little influence on our results. However, our study also had limitations, as 25(OH)D3 level in a population could be affected by many other factors such as sunshine duration, season, area, economy, and education, which was not considered in our study. These factors may affect the conclusions of epidemiological research, and bring bias to the meta-analysis. Therefore, we can only indicate that vitamin D insufficiency was related to HT. Whether vitamin D insufficiency could lead to HT should be further investigated by biological research in the future.

      Until recently, further investigations focussed on the HT mechanism in which vitamin D was involved. As is well understood, T lymphocytes including Th1, Th2 and Th17 cells infiltrate the thyroid gland due to immunological disorders in HT patients. Vitamin D can inhibit the differentiation of Th1 cells, and the production of inflammatory cytokines such as TNF-α, INF-γ. It could also suppress inflammatory Th1, but induce anti-inflammatory Th2 which produced anti-inflammatory cytokines such as IL-4 and IL-5[66]. Furthermore, the Th17 cells with their production of IL-17A could also be inhibited by Vitamin D at the transcriptional level[67]. On the other hand, vitamin D could increase the proportion of Treg cells to exert immune regulation[68]. Taken together, vitamin D may have the potential to prevent HT. However, clinical studies have shown contradictory results: Chahardoli et al.[69] reported activated vitamin D supplementation can decrease TSH and TG-Ab antibodies levels in HT patients, but another study showed that activated vitamin D supplementation had no effect on improving HT[70]. To our knowledge, the studies mentioned above may ignore the vitamin D receptors polymorphism. Vitamin D receptors in the thyroid gland have single nucleotide polymorphism and most typical Apal, Bsml, Fokl and Taql single nucleotide variations have been shown to be closely related to autoimmune diseases[71]. Therefore, more prospective studies are needed to confirm the preventive effect of vitamin D on HT.

    • In conclusion, the present systematic review and meta-analysis strengthened the relationship between vitamin D insufficiency and HT. HT patients potentially had higher propensity for having lower serum 25(OH)D3 levels compared to healthy individuals. Clinical staff may have to carefully consider the possibility of vitamin D insufficiency in HT patients.

      • The authors declare that they have no conflict of interest.

      • Copyright: © 2023 by the author(s). Published by Maximum Academic Press on behalf of Nanjing Agricultural University. This article is an open access article distributed under Creative Commons Attribution License (CC BY 4.0), visit https://creativecommons.org/licenses/by/4.0/.
    Figure (8)  Table (1) References (71)
  • About this article
    Cite this article
    Liu Z, Feng L, He Y, Yuan S, Xu C. 2022. The Association between Vitamin D and Hashimoto Thyroiditis: An Up-to-date Systematic Review and Meta-analysis. Food Materials Research 2:9 doi: 10.48130/FMR-2022-0009
    Liu Z, Feng L, He Y, Yuan S, Xu C. 2022. The Association between Vitamin D and Hashimoto Thyroiditis: An Up-to-date Systematic Review and Meta-analysis. Food Materials Research 2:9 doi: 10.48130/FMR-2022-0009

Catalog

  • About this article

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return