Front. Endocrinol.Frontiers in EndocrinologyFront. Endocrinol.1664-2392Frontiers Media S.A.10.3389/fendo.2022.860286EndocrinologyOriginal ResearchVitamin D Deficiency During Development Permanently Alters Liver Cell Composition and FunctionLundyKassidy1GreallyJohn F.2Essilfie-BondzieGrace1OlivierJosephine B.1Doña-TermineReanna1GreallyJohn M.1SuzukiMasako1*1Department of Genetics, Albert Einstein College of Medicine, Bronx, NY, United States2Retired, Galway, Ireland
Edited by: Takahiro Nemoto, Nippon Medical School, Japan
Reviewed by: Nobuaki Shiraki, Tokyo Institute of Technology, Japan; Satoru Ikenoue, Keio University School of Medicine, Japan
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.
Vitamin D, a fat-soluble vitamin, plays a critical role in calcium homeostasis, the immune system, and normal development. Many epidemiological cohort studies globally have found high prevalence rates of vitamin D deficiency and insufficiency, recognized as an important health issue that needs to be solved. In particular, reproductive age and pregnant women low in vitamin D status may confer risks of diseases like obesity on their offspring. While observational studies have suggested associations between prenatal vitamin D deficiency and metabolic phenotypes in offspring, not yet determined is whether prenatal vitamin D deficiency permanently alters the development of the liver, a major metabolic organ. We tested the histopathology and the transcriptomic profiles of livers from male C57BL/6J mice exposed to prenatal vitamin D deficiency through a maternal dietary intervention model. We found that prenatal vitamin D deficiency increases the prevalence of histopathological changes in the liver, and alters its gene expression profile. Cell subtype proportion analysis showed that the liver of prenatal vitamin D deficiency alters non-parenchymal cells of the liver, specifically macrophages, a subset of endothelial cells, and dendritic cells. Our results indicate the long-term memory of prenatal vitamin D deficiency exposure in the adult liver, a potential contributor to offspring health risks.
vitamin D deficiencyprenatal environmentlivertranscriptional alterationsDOHaD (developmental origins of health and disease)R01HL145302National Institutes of Health10.13039/100000002Introduction
Vitamin D is a crucial micronutrient that plays many physiologic functions in mammals. Humans create vitamin D from the action of sunlight on the skin, with the consequence that limiting sun exposure, covering skin with clothes, or darker skin pigmentation reduces vitamin D production and may cause its deficiency. The serum 25-hydroxyvitamin D concentration is the marker of vitamin D nutritional adequacy. Serum 25-hydroxyvitamin D concentrations in the range of 50-125 nmol/L indicate adequate vitamin D intake, 30-49 nmol/L is inadequate, and less than 30 nmol/L is classified as deficient, while greater than 125 nmol/L represents excess (1, 2). Many countries provide recommended daily vitamin D intake values and fortified foods such as milk, margarine, flours, cereals, and juices to reduce the prevalence rate of vitamin D deficiency (VDD), but VDD remains highly prevalent worldwide (3–6). In particular, the prevalence of VDD in reproductive-age women and pregnant women is recognized to be an urgent issue that needs to be addressed in low- and mid-income countries, and in high-income countries within the specific race and ethnicity groups (6–8).
While a major role of vitamin D is in regulating calcium homeostasis [reviewed in (9)], it has broader physiological effects. The bioactive metabolite of vitamin D, 1,25(OH)2D3 (calcitriol), is a steroid hormone, influencing transcriptional regulation through interaction with the high-affinity nuclear vitamin D receptor (VDR), a member of the nuclear receptor superfamily of ligand-activated transcription factors. The VDR pathway regulates the transcription of many genes, including genes with essential roles in fetal development and ensures the fetal supply of calcium for bone development (10–13). Intriguingly, epidemiological studies exploring the effects of maternal vitamin D deficiency revealed that offspring exposed to vitamin D deficiency in the uterus increased not only the risk of infantile rickets (4) but also the risk of childhood obesity (14–17) later in adulthood, suggesting long-term memories of the deficiency while developing in utero.
A mouse dietary intervention study revealed that mice born to mothers fed a vitamin D deficient diet had rapid weight gain after weaning and were more susceptible to high-fat diet-induced adipocyte hypertrophy than those born to mothers fed a vitamin D sufficient diet (18). In addition, a study of mice exposed to prenatal VDD found that prenatal VDD alters the DNA methylation status of the liver and sperm of two generations of offspring (19). These findings suggest that prenatal vitamin D deficiency predisposes offspring to long-term metabolic phenotypic alterations. However, little is known if prenatal vitamin D deficiency permanently alters the physiology of the liver, a major metabolic organ. To fill this knowledge gap, we assessed the effects of prenatal vitamin D deficiency on liver physiology and gene expression at the adult stage.
Materials And MethodsMaternal Vitamin D Deficiency Mouse Model
Five week old C57Bl/6J (Strain #000664) female mice (F0) purchased from the Jackson Laboratory were fed vitamin D deficient (VDD) or sufficient (VDS) diets for five weeks before mating to the control diet-fed male mice. The assigned diets were maintained throughout the subsequent pregnancy. VDD (0.0 IU/g vitamin D) and nutrient-matched VDS (1.0 IU/g vitamin D, D12450J) diets were obtained from Research Diets Inc. (10 kcal%fat, 20 kcal%Protein, and 70 kcal%Carbonate). VDD-fed females were supplemented with 1.5% calcium gluconate to maintain calcium homeostasis (20). Vitamin D deficiency of the F0 animals was confirmed by measuring serum 25-OH vitamin D levels using the Mouse Rat 25-OG Vitamin D ELISA kit (Eagle BIOSCIENCES) before mating and postnatal day 1 offspring (Supplementary Figure 1A). After delivery, all F0 mice were fed VDS diets. After weaning, the offspring (F1) were fed a VDS diet that was maintained until animal sacrifice. All animal studies were approved by the Institutional Animal Care and Use Committee at the Albert Einstein College of Medicine.
Histopathological Study
After dissection, the liver samples were fixed in 4% paraformaldehyde solution at 4°C for overnight and processed for paraffin embedding. The paraffin-embedded tissue slides were stained with hematoxylin and eosin (HE), trichrome, and Gomori Reticulin stainings. We assessed 10 VDS-F1 and 19 VDD-F1 male mice.
Transcription Analysis
We used male offspring for gene expression studies (n=6 VDD-F1 and n=6 VDS-F1). Total RNA was extracted from mouse livers at 16 weeks of age using the AllPrep DNA/RNA Micro Kits (QIAGEN). The RNA-seq libraries were generated using KAPA RNA HyperPrep with RiboErase kit (Roche) and sequenced on an Illumina NovaSeq sequencer (Novogene Co., Ltd., USA). After removing reads that failed the quality check and trimming adapter sequences, the resulting sequences were aligned to the mouse mm10 reference genome, quantifying gene expression by transcript counting (GENECODE release M15, GRCm38) using the STAR aligner (21). Differential expression as a result of dietary manipulation was determined using the Bioconductor package DESeq2 (22). We eliminated genes with fewer than the mean read counts per sample per gene is less than 12 from the analysis. Significantly differentially expressed genes (DEGs) were determined based on a log2-fold change of less than -1 or greater than 1 and a false discovery rate (FDR) adjusted p-value less than 0.05. DEGs were further assessed for their biological properties by Gene Ontology (GO) enrichment analysis using a Bioconductor package, clusterProfiler (23, 24). Quantitative reverse transcription PCR (qRT-PCR) was performed on the same RNA samples to verify the results. We used SuperScript III with random hexamer priming to synthesize cDNA and LightCycler FastStart Universal SYBR Green Master Mix for quantitative PCR. The measurement was performed using the LightCycler 480 system. The primers used in this study are listed in Supplementary Table 1.
Cell Subtype Deconvolution
We used CIBERSORTx, a web interface-based tool, for quantifying the proportion of cell types in the RNA-seq data (25). A publicly available mouse liver (parenchymal and non-parenchymal liver cells isolated by a two-step collagenase perfusion) single-cell RNA-seq data (GSE134134) (26, 27) was used to generate a custom signature from single-cell liver gene expression data that represents the gene expression signatures of liver cell and immune cell subtypes. We combined four datasets of single-cell RNA seq libraries listed under GSE134134, 1000 (GSM3937757), 2000 (GSM3937758), 5000 (GSM3937759), and 10000 (GSM3937760) (26, 27). We used the Seurat R package for the single-cell RNA-seq analysis (28). We excluded cells with the detected number of genes were less than 1000 and reads from mitochondrial were greater than 10% before the analysis. We performed normalization and variance stabilization using the sctransform algorism (29), and the sctransformed datasets were merged using integration functions and identified clusters based on the gene expression profiles (29). We identified cell subtype clusters based on the expression status of maker genes (Supplementary Figure 4).
Statistical Analysis
All statistical analyses were performed using R version 4.1.1 (https://www.r-project.org/). The Student’s t-test was used to analyze continuous variables between VDS and VDD unless otherwise stated, and Fisher’s exact test to analyze categorical variables. P-values of <0.05 or FDR-adjusted p-values in multiple testing analyses were considered significant.
ResultsPrenatal Vitamin D Deficiency Is Associated With Histopathological Changes of the Livers of Offspring
The growth trajectory of offspring was comparable in VDD-F1 and VDS-F1 (Supplementary Figure 1B). We dissected F1 male mice at 16 weeks of age to assess the long-term effects of prenatal vitamin D deficiency. We found histopathological alterations of the livers including hydropic degenerations, steatosis, and fibrosis in VDD-F1 offspring (Figure 1). While the prevalence rates were not statistically significant due to the small sample numbers, steatosis and fibrosis were only observed in the VDD-F1 offspring. The sizes of lipid droplets were varied between the individuals, but the morphology of the nuclei remained unaltered. We summarize the histopathological alterations in Table 1.
Representative histopathological alterations in the liver of VDD-F1 offspring. Top, hematoxylin-eosin stain; bottom, Gomori’s reticulin stain. The arrows indicate the described changes.
A summary of histopathological alterations in liver.
VDD-F1 (n=19)
VDS-F1 (n=10)
p-values*
Hydropic degeneration
11
4
0.675
Steatosis
1
0
0.063
Steatosis and Fibrosis
5
0
–
No alterations observed
8
6
0.562
*Fisher’s exact test.
Prenatal Vitamin D Deficiency Is Associated With Permanently Altered Gene Expression of the Liver
The histopathological analyses suggested that prenatal vitamin D deficiency alters liver physiology. We, therefore, conducted a transcriptional analysis using bulk RNA-seq to assess the gene expression alterations linked to prenatal vitamin D deficiency. We extracted total RNA from the left lateral lobes of the liver and prepared libraries (n=6 per group, Supplementary Figure 2). We sequenced libraries at the depth of at least 26 million paired reads per library. After the quality check of the sequences and removal of low-quality reads, we performed a principal components analysis to compare the expression pattern between the individual F1 mice (Figure 2A, Supplementary Table 2). We observed three large clusters that consist of only VDS samples, VDD samples, and both VDS and VDD. There were no significant histopathological differences between VDD samples in different clusters. Next, we performed a differential expression analysis and identified 281 significant differentially expressed genes (DEGs) when comparing samples from mice who were exposed to prenatal VDD versus those exposed to prenatal VDS. 249 of DEGs were upregulated and 32 of DEGs were downregulated (Figure 2B, Supplementary Table 3). We verified the results using the same RNA samples and confirmed significant alterations or the same trends of DEGs by quantitative RT-PCR (Figure 2C, Supplementary Figure 3).
Prenatal vitamin D deficiency leads to altered gene expression in adult mouse liver. (A) PCA demonstrates variation in expression patterns between VDD and VDS offspring. (B) A volcano plot shows differentially expressed genes (DEGs) found in the analysis comparing mouse offspring exposed to a prenatal vitamin D deficiency versus mice exposed to a vitamin D sufficient diet. Genes are significantly differentially expressed if they have a p adjusted value less than 0.05 and -2 < fold change < 2. (C) qRT-PCR using cDNA generated from offspring liver samples validates the upregulation of the collagen, type III, alpha I (Col3a1), and Cidec genes.
The DEGs Are Enriched in Pathways for the Extracellular Structural Organization and Lipid Catabolism
We performed gene ontology (GO) enrichment analyses on upregulated and downregulated DEGs respectively to test if the DEGs are enriched in specific gene pathways. We found that lipid metabolism-related pathways such as lipid catabolic process, xenobiotic metabolism, and linoleic acid metabolic process are enriched in down-regulated DEGs. 249 upregulated DEGs were enriched in extracellular structure organization, nuclear division, and mitosis-related pathways (Figure 3, Supplementary Table 4).
Results from GO enrichment analysis (biological processes) using a Bioconductor package, clusterProfiler, for up and down-regulated DEGs. The top 15 up (right) and downregulated (left) terms were. The size of the dots represents the number of the genes, and the color indicates the significance of the enrichment.
Cell Subtype Proportion Analysis of Bulk RNA-Seq Data Revealed That Prenatal Vitamin D Deficiency Alters Cell Subtype Composition of Non-Parenchymal Liver Cells
We next tested if the proportions of cell subtypes are affected based on prenatal exposure to VDS or VDD. To perform cell subtype proportion analysis, we generated reference expression signatures of each liver cell subtype using publicly available adult liver single-cell RNA-seq data. We segregated 21 clusters of the adult liver cells, as shown in Figure 4A and Supplementary Figure 4. Based on the relative expression of well-known marker genes of liver cell subtypes identified by the mouse cell atlas, we identified the cell subtypes of each cluster. The cluster information was then used to generate a custom signature of gene expression in liver cells using CIBERSORTx. The generated reference expression signature was used to perform deconvolution of bulk RNA-seq data of each sample (Supplementary Table 5). As we expected, the hepatocyte proportion was the highest among the cell subtypes (80.6% VDD and 81.4% VDS, p=0.39) (Supplementary Table 6). While the proportions of hepatocytes were not significantly altered, we observed significant alterations in non-parenchymal liver cell subtypes, macrophages, endothelial cells, and dendritic cells (Figure 4B). Interestingly, while Cd163-positive macrophages were reduced in prenatal VDD offspring liver, Ccr2-positive macrophages were increased compared to that of VDS offspring. We previously reported that cell subtype composition significantly contributes to the gene expression profile (30). We therefore tested if the DEGs we identified were attributable to the cell subtype proportion differences between samples. We did not observe significant contributions of cell subtype proportion variations to the gene expression profile (Supplementary Figure 5). This result indicates that the identified DEGs were independent of the cell subtype proportion variations.
Liver cell subtypes were found in different cell subtype proportions in prenatal VDD compared to VDS offspring. (A) Publicly available single-cell RNA-seq data was processed and clustered as shown in the UMAP plot. Gene expression patterns of clusters were used to define cell identities. (B) The proportions of Cd163+ macrophages and Ccr2+macrophages were inversely related between prenatal VDD and VDS offspring. Dendritic cells were found in a significantly lower proportion in VDD samples while a subset of endothelial cells was found in a higher proportion in VDD samples.
Discussion
In this study, we found that prenatal vitamin D deficiency alters liver histology and transcription profiles in offspring adulthood. We observed livers with varying ranges of cytoplasmic clearing, hydropic changes, and fibrosis in histopathological analysis. Transcriptional analysis showed that lipid metabolism-related genes, mitosis-related genes, and extracellular structure organization-related genes were enriched in DEGs. Moreover, cell subtype deconvolution analysis indicates alterations of non-parenchymal liver cell proportions. These findings suggest that prenatal vitamin D deficiency permanently alters the physiological and metabolic properties of offspring livers.
Herrick et al. reported that more than 23% of Americans equal or older than 1 year are at risk of vitamin D inadequacy or deficiency and the prevalence of the risk of deficiency is at its highest among adults among 20-39 years old using 2011-2014 National Health and Nutrition Examination Survey (NHANES) data sets (1). They also reported that the risk of deficiency was higher among non-Hispanic Black than non-Hispanic Asian, non-Hispanic White, and Hispanic participants. Moreover, the high prevalence of vitamin D deficiency and inadequacy is associated with obesity (5, 31, 32). These findings suggest that reproductive-age women with obesity have a higher risk of vitamin D deficiency during pregnancies, specifically women of the vitamin D deficiency vulnerable ethnic/race groups.
The associations between maternal vitamin D deficiency and adverse maternal and fetal outcomes have been well reported (33–37). A systematic review of 3357 studies revealed that maternal vitamin D deficiency increases the risks of gestational diabetes, preeclampsia, and small for gestational age infants (35). In addition, prior studies have implicated associations between maternal vitamin D deficiency and offspring obesity later in life in both human observation studies (14–17, 38) and dietary manipulation studies in rodents (18, 39). We did not observe an obese phenotype; however, we found that mice exposed to prenatal vitamin D deficiency showed transcriptional alterations in lipid metabolism-related pathways. Moreover, the histopathological analysis suggested that exposure to a prenatal vitamin D deficiency affects the physiology of the liver at the adult stage. Our results suggest the existence of long-term memory of prenatal vitamin D deficiency exposure in the liver.
The association between vitamin D deficiency in adulthood and chronic liver diseases (CLD) (40), including cirrhosis (41), non-alcoholic fatty liver disease (NAFLD, recently been proposed to rename as metabolic-associated fatty liver disease [MAFLD) (42)] (43–45), and non-alcoholic steatohepatitis (NASH) (44, 45), have been identified in many cross-sectional studies. Stokes et al. reviewed the associations between CLD and vitamin D deficiency (46). We fed a vitamin D sufficient diet to all offspring after weaning, finding that the serum vitamin D concentration of mice exposed to prenatal VDD became normal by 5 weeks of age. Therefore, the pathogenesis of the liver phenotype we observed in the liver of prenatal VDD offspring was not by deficiency at the time, likely by a long-term memory of VDD exposure. However, there are several similarities that exist. For instance, the infiltration and activation of macrophages is a key feature of endoplasmic reticulum stress-induced inflammation in liver injury, with VDR signaling regulating the inflammatory response in macrophages (47). We found that two distinct macrophage proportions, Cd163-positive, the anti-inflammatory (M2 phenotype) macrophages (48), and Ccr2-positive macrophages had distinct patterns of alteration in prenatal VDD offspring. Associations between the Ccr2-positive Kupffer cells (hepatic residential macrophage) and choline-deficient amino acid-defined diet-induced hepatic steatosis and fibrosis have been identified (49, 50). Several studies reported that a high-fat diet-induced hepatic steatosis increases bone-marrow (BM)-derived macrophages, which predominantly express Ccr2 (51, 52). The cell subtype reference expression signatures of both Cd163-positive and Ccr2-positive macrophages in this study showed higher expression of Emr1 (F4/80) but low expression of Itgam (Cd11b) (52), indicating that they may be tissue-resident macrophages. Ccr2 expression in Kupffer cells contributes to BM-derived macrophage recruitment (49). The resident macrophages are derived from multiple anatomical locations during development suggesting developmental origin cell subtype proportion alterations may exist in the liver that is associated with a predisposing condition of prenatal VDD to the liver phenotype at the adult stage. The GO enrichment analyses on DEGs in this study showed enrichment of biological processes GO terms including extracellular matrix organization, regulation of vascularized development, and regulation of angiogenesis, processes which occur with liver fibrosis. Our study provides new insight into how prenatal VDD affects the liver phenotype and physiology as well as how this results in changes in the proportions of liver-resident immune cell subtypes in adult offspring.
There are several limitations of the current study that we would like to address in the future. First, we performed dietary manipulation on F0 female mice for five weeks prior to mating until the delivery. Belenchia et al. reported that it takes at least two weeks for serum vitamin D levels to return to normal after switching the diet to VDS from long-term VDD feeding (53). Thus, the F1 offspring were exposed to vitamin D deficient and insufficient conditions for the first two weeks of life. Using foster mothers would help us dissociate the effects between prenatal and neonatal periods. Second, while we observed cell subtype proportion changes in the prenatal VDD offspring, we cannot conclude if this alteration contributes to prenatal VDD exposure or results from liver damage. We believe testing cell subtype proportions in the liver before liver injury started would resolve this question.
Moreover, we estimated cell subtype proportion changes by the bulk RNA-seq results; therefore, the estimated proportion may be affected by the expression changes in specific cell subtypes. Testing the expression profiles at the single-cell levels will test if the alterations we observed in this study are attributed to the gene expression alterations between the same cell subtypes. Lastly, Wang et al. estimated the relationship between mice and humans at stages of growth and reported that the age range of 10 to 64.29 weeks of age in mice corresponds to the adult stage in humans (20 to 51 years old) (54), suggesting 16 weeks of age is at the early adult stage. The prevalence of NAFLD increases after age 30 (55) and age older than 40 is a risk factor of NAFLD (56). This study could be improved by maintaining the mice for longer than 16 weeks of age following prenatal VDD in order to observe more changes related to prenatal VDD and at the adulthood stage.
In summary, we found that prenatal vitamin D deficiency alters gene expression profiles of the liver and cell subtype proportions of non-parenchymal cells at the adult stage. The prevalence of at risk of vitamin D deficiency and insufficiency is high among reproductive-age women, and a significant association of obesity to low vitamin D status has been implicated. Childhood obesity in the world is rising in frequency, with the prevalence of overweight/obesity among children aged 5-18 in the world rising from 4% in 1974 to 18% in 2016 (57). Epidemiological studies have identified significant associations between excess weight during childhood and the risk of adulthood obesity (58–62). Moreover, the incidence of NASH and NAFLD/MAFLD in young people has increased to 1.52 times between 1990 and 2017, and the prevalence of NAFLD/MAFLD in children varies by ethnicity (56). Prenatal and neonatal vitamin D deficiency may be contributing to this increase in prevalence rates. Reducing maternal vitamin D deficiency could reduce the adverse health outcome of the next generation, while understanding mechanisms of how our body remembers prenatal vitamin D status over the lifespan should help us to find interventions and treatments to prevent developmental origins of health and diseases in the liver.
Data Availability Statement
The original contributions presented in the study are publicly available. This data can be found here: GEO, GSE194417.
Ethics Statement
The animal study was reviewed and approved by Institutional Animal Care and Use Committee at the Albert Einstein College of Medicine.
Author Contributions
Conceptualization, MS. Methodology, KL, RD-T, and MS. Formal analysis, KL and GE-B. Histopathological analysis support, JO. Histopathological analysis, JFG. Data curation, KL and MS. Writing—original draft preparation, KL and MS. Writing—review and editing, KL, JFG, RD-T, JMG, and MS. Visualization, KL, GE-B, and MS. Supervision, JFG, JMG, and MS. Funding acquisition, MS. All authors have read and agreed to the published version of the manuscript.
Funding
This work was supported by the Human Genome Project Pilot Grant of the Department of Genetics at Albert Einstein College of Medicine and partially supported by the National Institutes of Health under award number R01HL145302 (MS).
Author Disclaimer
The content is solely the responsibility of the authors and does not necessarily represent the official views of the National Institutes of Health.
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.
Acknowledgments
We thank to Analytical Imaging Facility (AIF) of Albert Einstein College of Medicine for their support. All histology slides were scanned with a P250 slide scanner supported by Shared Instrumentation Grant, 1S10OD019961-01. The AIF is partially funded by Cancer Center Support Grant of National Cancer Institutes (P30CA013330).
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fendo.2022.860286/full#supplementary-material
ReferencesHerrickKAStorandtRJAffulJPfeifferCMSchleicherRLGahcheJJ. Vitamin D Status in the United States, 2011-2014. Am J Clin Nutr (2019) 110:150–7. doi: 10.1093/ajcn/nqz037RossACTaylorCLYaktineALDel ValleHB, eds. Institute of Medicine (US) Committee to Review Dietary Reference Intakes for Vitamin D and Calcium. In: Dietary Reference Intakes for Calcium and Vitamin D. Washington, DC): National Academies Press (US. doi: 10.17226/13050LipsP. Worldwide Status of Vitamin D Nutrition. J Steroid Biochem Mol Biol (2010) 121:297–300. doi: 10.1016/j.jsbmb.2010.02.021ÖzdemirAAErcan GündemirYKüçükMYıldıran SarıcıDElgörmüşYÇağY. Vitamin D Deficiency in Pregnant Women and Their Infants. J Clin Res Pediatr Endocrinol (2018) 10:44–50. doi: 10.4274/jcrpe.4706LiuXBaylinALevyPD. Vitamin D Deficiency and Insufficiency Among US Adults: Prevalence, Predictors and Clinical Implications. Br J Nutr (2018) 119:928–36. doi: 10.1017/S0007114518000491RothDEAbramsSAAloiaJBergeronGBourassaMWBrownKH. Global Prevalence and Disease Burden of Vitamin D Deficiency: A Roadmap for Action in Low- and Middle-Income Countries. Ann N Y Acad Sci (2018) 1430:44–79. doi: 10.1111/nyas.13968KanataniKTNakayamaTAdachiYHamazakiKOnishiKKonishiY. High Frequency of Vitamin D Deficiency in Current Pregnant Japanese Women Associated With UV Avoidance and Hypo-Vitamin D Diet. PloS One (2019) 14:e0213264. doi: 10.1371/journal.pone.0213264JohnsonDDWagnerCLHulseyTCMcNeilRBEbelingMHollisBW. Vitamin D Deficiency and Insufficiency is Common During Pregnancy. Am J Perinatol (2011) 28:7–12. doi: 10.1055/s-0030-1262505DeLucaHF. Evolution of Our Understanding of Vitamin D. Nutr Rev (2008) 66:S73–87. doi: 10.1111/j.1753-4887.2008.00105.xYuSBruceDFroicuMWeaverVCantornaMT. Failure of T Cell Homing, Reduced CD4/CD8alphaalpha Intraepithelial Lymphocytes, and Inflammation in the Gut of Vitamin D Receptor KO Mice. Proc Natl Acad Sci USA (2008) 105:20834–9. doi: 10.1073/pnas.0808700106YoshizawaTHandaYUematsuYTakedaSSekineKYoshiharaY. Mice Lacking the Vitamin D Receptor Exhibit Impaired Bone Formation, Uterine Hypoplasia and Growth Retardation After Weaning. Nat Genet (1997) 16:391–6. doi: 10.1038/ng0897-391LiYCPirroAEAmlingMDellingGBaronRBronsonR. Targeted Ablation of the Vitamin D Receptor: An Animal Model of Vitamin D-Dependent Rickets Type II With Alopecia. Proc Natl Acad Sci USA (1997) 94:9831–5. doi: 10.1073/pnas.94.18.9831DardenneOPrud’hommeJGlorieuxFHSt-ArnaudR. Rescue of the Phenotype of CYP27B1 (1alpha-Hydroxylase)-Deficient Mice. J Steroid Biochem Mol Biol (2004) 89–90:327–30. doi: 10.1016/j.jsbmb.2004.03.026MoralesERodriguezAValviDIñiguezCEspluguesAVioqueJ. Deficit of Vitamin D in Pregnancy and Growth and Overweight in the Offspring. Int J Obes (Lond) (2015) 39:61–8. doi: 10.1038/ijo.2014.165GaleCRRobinsonSMHarveyNCJavaidMKJiangBMartynCN. Maternal Vitamin D Status During Pregnancy and Child Outcomes. Eur J Clin Nutr (2008) 62:68–77. doi: 10.1038/sj.ejcn.1602680DarakiVRoumeliotakiTChalkiadakiGKatrinakiMKarachaliouMLeventakouV. Low Maternal Vitamin D Status in Pregnancy Increases the Risk of Childhood Obesity. Pediatr Obes (2018) 13:467–75. doi: 10.1111/ijpo.12267TintMTChongMFArisIMGodfreyKMQuahPLKapurJ. Association Between Maternal Mid-Gestation Vitamin D Status and Neonatal Abdominal Adiposity. Int J Obes (Lond) (2018) 42:1296–305. doi: 10.1038/s41366-018-0032-2BelenchiaAMJohnsonSAEllersieckMRRosenfeldCSPetersonCA. In Utero Vitamin D Deficiency Predisposes Offspring to Long-Term Adverse Adipose Tissue Effects. J Endocrinol (2017) 234:301–13. doi: 10.1530/JOE-17-0015XueJSchoenrockSAValdarWTarantinoLMIderaabdullahFY. Maternal Vitamin D Depletion Alters DNA Methylation at Imprinted Loci in Multiple Generations. Clin Epigenet (2016) 8:107. doi: 10.1186/s13148-016-0276-4NicholasCDavisJFisherTSegalTPettiMSunY. Maternal Vitamin D Deficiency Programs Reproductive Dysfunction in Female Mice Offspring Through Adverse Effects on the Neuroendocrine Axis. Endocrinology (2016) 157:1535–45. doi: 10.1210/en.2015-1638DobinADavisCASchlesingerFDrenkowJZaleskiCJhaS. STAR: Ultrafast Universal RNA-Seq Aligner. Bioinformatics (2013) 29:15–21. doi: 10.1093/bioinformatics/bts635LoveMIHuberWAndersS. Moderated Estimation of Fold Change and Dispersion for RNA-Seq Data With Deseq2. Genome Biol (2014) 15:550. doi: 10.1186/s13059-014-0550-8YuGWangL-GHanYHeQ-Y. Clusterprofiler: An R Package for Comparing Biological Themes Among Gene Clusters. OMICS (2012) 16:284–7. doi: 10.1089/omi.2011.0118WuTHuEXuSChenMGuoPDaiZ. Clusterprofiler 4.0: A Universal Enrichment Tool for Interpreting Omics Data. Innovation (NY) (2021) 2:100141. doi: 10.1016/j.xinn.2021.100141NewmanAMSteenCBLiuCLGentlesAJChaudhuriAASchererF. Determining Cell Type Abundance and Expression From Bulk Tissues With Digital Cytometry. Nat Biotechnol (2019) 37:773–82. doi: 10.1038/s41587-019-0114-2BaierFASánchez-TaltavullDYarahmadovTCastellàCGJebbawiFKeoghA. Loss of Claudin-3 Impairs Hepatic Metabolism, Biliary Barrier Function, and Cell Proliferation in the Murine Liver. Cell Mol Gastroenterol Hepatol (2021) 12:745–67. doi: 10.1016/j.jcmgh.2021.04.003Sanchez-TaltavullDPerkinsTJDommannNMelinNKeoghACandinasD. Bayesian Correlation Is a Robust Gene Similarity Measure for Single-Cell RNA-Seq Data. NAR Genom Bioinform (2020) 2:lqaa002. doi: 10.1093/nargab/lqaa002SatijaRFarrellJAGennertDSchierAFRegevA. Spatial Reconstruction of Single-Cell Gene Expression Data. Nat Biotechnol (2015) 33:495–502. doi: 10.1038/nbt.3192HafemeisterCSatijaR. Normalization and Variance Stabilization of Single-Cell RNA-Seq Data Using Regularized Negative Binomial Regression. Genome Biol (2019) 20:296. doi: 10.1186/s13059-019-1874-1KongYRastogiDSeoigheCGreallyJMSuzukiM. Insights From Deconvolution of Cell Subtype Proportions Enhance the Interpretation of Functional Genomic Data. PloS One (2019) 14:e0215987. doi: 10.1371/journal.pone.0215987BurrisHHThomasAZeraCAMcElrathTF. Prenatal Vitamin Use and Vitamin D Status During Pregnancy, Differences by Race and Overweight Status. J Perinatol (2015) 35:241–5. doi: 10.1038/jp.2014.198TurerCBLinHFloresG. Prevalence of Vitamin D Deficiency Among Overweight and Obese US Children. Pediatrics (2013) 131:e152–61. doi: 10.1542/peds.2012-1711BodnarLMCatovJMSimhanHNHolickMFPowersRWRobertsJM. Maternal Vitamin D Deficiency Increases the Risk of Preeclampsia. J Clin Endocrinol Metab (2007) 92:3517–22. doi: 10.1210/jc.2007-0718WeiSQAudibertFHidiroglouNSarafinKJulienPWuY. Longitudinal Vitamin D Status in Pregnancy and the Risk of Pre-Eclampsia. BJOG (2012) 119:832–9. doi: 10.1111/j.1471-0528.2012.03307.xAghajafariFNagulesapillaiTRonksleyPEToughSCO’BeirneMRabiDM. Association Between Maternal Serum 25-Hydroxyvitamin D Level and Pregnancy and Neonatal Outcomes: Systematic Review and Meta-Analysis of Observational Studies. BMJ (2013) 346:f1169. doi: 10.1136/bmj.f1169van der PligtPWillcoxJSzymlek-GayEAMurrayEWorsleyADalyRM. Associations of Maternal Vitamin D Deficiency With Pregnancy and Neonatal Complications in Developing Countries: A Systematic Review. Nutrients (2018) 10(5):640. doi: 10.3390/nu10050640GilaniSJanssenP. Maternal Vitamin D Levels During Pregnancy and Their Effects on Maternal-Fetal Outcomes: A Systematic Review. J Obstet Gynaecol Can (2019) 42(9):1129–37. doi: 10.1016/j.jogc.2019.09.013BoyleVTThorstensenEBThompsonJMDMcCowanLMEMitchellEAGodfreyKM. The Relationship Between Maternal 25-Hydroxyvitamin D Status in Pregnancy and Childhood Adiposity and Allergy: An Observational Study. Int J Obes (Lond) (2017) 41:1755–60. doi: 10.1038/ijo.2017.182WenJHongQWangXZhuLWuTXuP. The Effect of Maternal Vitamin D Deficiency During Pregnancy on Body Fat and Adipogenesis in Rat Offspring. Sci Rep (2018) 8:365. doi: 10.1038/s41598-017-18770-4HeXXuCLuZ-HFangX-ZTanJSongY. Low Serum 25-Hydroxyvitamin D Levels are Associated With Liver Injury Markers in the US Adult Population. Public Health Nutr (2020) 23:2915–22. doi: 10.1017/S1368980020000348KubeschAQuenstedtLSalehMRüschenbaumSSchwarzkopfKMartinezY. Vitamin D Deficiency Is Associated With Hepatic Decompensation and Inflammation in Patients With Liver Cirrhosis: A Prospective Cohort Study. PloS One (2018) 13:e0207162. doi: 10.1371/journal.pone.0207162EslamMSanyalAJGeorgeJ.International Consensus Panel. MAFLD: A Consensus-Driven Proposed Nomenclature for Metabolic Associated Fatty Liver Disease. Gastroenterology (2020) 158:1999–2014.e1. doi: 10.1053/j.gastro.2019.11.312ParkEParkEY. Inverse Association Between Serum 25-Hydroxyvitamin D Levels and Risk of Suspected Non-Alcoholic Fatty Liver Disease in Obese Population. Int J Environ Res Public Health (2021) 18(16):8682. doi: 10.3390/ijerph18168682LivadariuRTimofteDTrifanADanilaRIonescuLSîngeapAM. Vitamin D Deficiency, a Noninvasive Marker of Steatohepatitis in Patients With Obesity and Biopsy Proven Nonalcoholic Fatty Liver Disease. Acta Endocrinol (Buchar) (2018) 14:76–84. doi: 10.4183/aeb.2018.76TargherGBertoliniLScalaLCigoliniMZenariLFalezzaG. Associations Between Serum 25-Hydroxyvitamin D3 Concentrations and Liver Histology in Patients With Non-Alcoholic Fatty Liver Disease. Nutr Metab Cardiovasc Dis (2007) 17:517–24. doi: 10.1016/j.numecd.2006.04.002StokesCSVolmerDAGrünhageFLammertF. Vitamin D in Chronic Liver Disease. Liver Int (2013) 33:338–52. doi: 10.1111/liv.12106ZhouYDongBKimKHChoiSSunZWuN. Vitamin D Receptor Activation in Liver Macrophages Protects Against Hepatic Endoplasmic Reticulum Stress in Mice. Hepatology (2020) 71:1453–66. doi: 10.1002/hep.30887WangHZhangC-SFangB-BHouJLiW-DLiZ-D. Dual Role of Hepatic Macrophages in the Establishment of the Echinococcus Multilocularis Metacestode in Mice. Front Immunol (2020) 11:600635. doi: 10.3389/fimmu.2020.600635MiuraKYangLvan RooijenNOhnishiHSekiE. Hepatic Recruitment of Macrophages Promotes Nonalcoholic Steatohepatitis Through CCR2. Am J Physiol Gastrointest Liver Physiol (2012) 302:G1310–21. doi: 10.1152/ajpgi.00365.2011SekiEde MinicisSInokuchiSTauraKMiyaiKvan RooijenN. CCR2 Promotes Hepatic Fibrosis in Mice. Hepatology (2009) 50:185–97. doi: 10.1002/hep.22952KleinICornejoJCPolakosNKJohnBWuenschSATophamDJ. Kupffer Cell Heterogeneity: Functional Properties of Bone Marrow Derived and Sessile Hepatic Macrophages. Blood (2007) 110:4077–85. doi: 10.1182/blood-2007-02-073841ObstfeldAESugaruEThearleMFranciscoA-MGayetCGinsbergHN. C-C Chemokine Receptor 2 (CCR2) Regulates the Hepatic Recruitment of Myeloid Cells That Promote Obesity-Induced Hepatic Steatosis. Diabetes (2010) 59:916–25. doi: 10.2337/db09-1403BelenchiaAMJohnsonSAKieschnickACRosenfeldCSPetersonCA. Time Course of Vitamin D Depletion and Repletion in Reproductive-Age Female C57BL/6 Mice. Comp Med (2017) 67:483–90.WangSLaiXDengYSongY. Correlation Between Mouse Age and Human Age in Anti-Tumor Research: Significance and Method Establishment. Life Sci (2020) 242:117242. doi: 10.1016/j.lfs.2019.117242EstesCRazaviHLoombaRYounossiZSanyalAJ. Modeling the Epidemic of Nonalcoholic Fatty Liver Disease Demonstrates an Exponential Increase in Burden of Disease. Hepatology (2018) 67:123–33. doi: 10.1002/hep.29466AmarapurkarDKamaniPPatelNGuptePKumarPAgalS. Prevalence of Non-Alcoholic Fatty Liver Disease: Population Based Study. Ann Hepatol (2007) 6:161–3. doi: 10.1016/S1665-2681(19)31922-2Le GarfSNègreVAntyRGualP. Metabolic Fatty Liver Disease in Children: A Growing Public Health Problem. Biomedicines (2021) 9(12):1915. doi: 10.3390/biomedicines9121915WhitakerRCWrightJAPepeMSSeidelKDDietzWH. Predicting Obesity in Young Adulthood From Childhood and Parental Obesity. N Engl J Med (1997) 337:869–73. doi: 10.1056/NEJM199709253371301VennAJThomsonRJSchmidtMDClelandVJCurryBAGennatHC. Overweight and Obesity From Childhood to Adulthood: A Follow-Up of Participants in the 1985 Australian Schools Health and Fitness Survey. Med J Aust (2007) 186:458–60. doi: 10.5694/j.1326-5377.2007.tb00997.xEvensenEWilsgaardTFurbergA-SSkeieG. Tracking of Overweight and Obesity From Early Childhood to Adolescence in a Population-Based Cohort - The Tromsø Study, Fit Futures. BMC Pediatr (2016) 16:64. doi: 10.1186/s12887-016-0599-5WangLYChyenDLeeSLowryR. The Association Between Body Mass Index in Adolescence and Obesity in Adulthood. J Adolesc Health (2008) 42:512–8. doi: 10.1016/j.jadohealth.2007.10.010LlewellynASimmondsMOwenCGWoolacottN. Childhood Obesity as a Predictor of Morbidity in Adulthood: A Systematic Review and Meta-Analysis. Obes Rev (2016) 17:56–67. doi: 10.1111/obr.12316