Multi-omics Analysis of Umbilical Cord Hematopoietic Stem Cells from a Multi-ethnic Cohort of Hawaii Reveals the Transgenerational Effect of Maternal Pre-Pregnancy Obesity
===========================================================================================================================================================================

* Yuheng Du
* Paula A. Benny
* Ryan J. Schlueter
* Alexandra Gurary
* Annette Lum-Jones
* Cameron B Lassiter
* Fadhl M. AlAkwaa
* Maarit Tiirikainen
* Dena Towner
* W. Steven Ward
* Lana X Garmire

## Abstract

**Background** Maternal obesity is a health concern that may predispose newborns to a high risk of medical problems later in life. To understand the transgenerational effect of maternal obesity, we conducted a multi-omics study, using DNA methylation and gene expression in the CD34+/CD38-/Lin-umbilical cord blood hematopoietic stem cells (uHSCs) and metabolomics of the cord blood, all from a multi-ethnic cohort (n=72) from Kapiolani Medical Center for Women and Children in Honolulu, Hawaii (collected between 2016 and 2018).

**Results** Differential methylation (DM) analysis unveiled a global hypermethylation pattern in the maternal pre-pregnancy obese group (BH adjusted p<0.05), after adjusting for major clinical confounders. Comprehensive functional analysis showed hypermethylation in promoters of genes involved in cell cycle, protein synthesis, immune signaling, and lipid metabolism. Utilizing Shannon entropy on uHSCs methylation, we discerned notably higher quiescence of uHSCs impacted by maternal obesity. Additionally, the integration of multi-omics data-including methylation, gene expression, and metabolomics-provided further evidence of dysfunctions in adipogenesis, erythropoietin production, cell differentiation, and DNA repair, aligning with the findings at the epigenetic level.

**Conclusions** This study reveals the significant correlation between pre-pregnancy maternal obesity and multi-omics level molecular changes in the uHSCs of offspring, particularly in DNA methylation.

Keywords
*   obesity
*   Native Hawaiian
*   Hematopoietic stem cells
*   Multi-omics
*   cord blood
*   methylation
*   pregnancy

## Introduction

Maternal obesity has emerged as a primary health concern during pregnancy, with its prevalence alarmingly increasing. According to a study by the Centers for Disease Control and Prevention, the percentage of women experiencing pre-pregnancy obesity in the United States escalated from 26% to 29% between 2016 and 2019 1. Born to mothers with obesity, higher birth weight is associated with a higher incidence of childhood cancers such as leukemia and neuroblastoma 2,3, as well as greater risks of prostate and testicular cancers in men 4–6 and breast cancer in women 7. Moreover, maternal obesity may have a transgenerational effect and set the stage for increased chronic disease susceptibility later in the adulthood of offspring 8,9. The hypothesis of the utero origin of diseases proposes that numerous chronic diseases have their origins in the fetal stage, the earliest phase of human development 10,11. Some researchers have speculated higher stem cell burdens in newborn babies born from obese mothers 12. In particular, a study showed increases in cord blood CD34+CD38-stem cell and CD34+ progenitor cell concentrations with maternal obesity 13, suggesting that the higher proportions of stem cells in cord blood may make the babies more susceptible to obesity and cancer risks. However, so far little work provides the direct molecular links as to how maternal obesity affects the cellular function and increases the disease risk in offspring.

To seek answers in this area, we conducted an epigenome centered multi-omics study to directly pinpoint the effect of maternal obesity in umbilical cord blood hematopoietic stem cells (uHSCs). Epigenetics is chosen as the center of multi-omics integration, as it is both inheritable and susceptible to modification by diseases. Thus it may serve as a plausible mediator in the transmission of the effects of maternal obesity to offspring. We incorporate gene expression for cord blood stem cells and metabolomics data from the cord blood serum as the downstream readout of epigenetics changes. By elucidating these molecular connections, we provide a systematic understanding of how maternal obesity during pregnancy can influence the multiple types of molecular profiles of newborns. Such knowledge may ultimately help develop early therapeutic interventions at the molecular level to mitigate these transgenerational health risks due to maternal obesity.

## Methods

### Overview of the maternal pre-pregnancy cohort with baby cord blood

In this study, baby cord blood samples from 72 pregnant women (34 obese; 38 non-obese) who delivered at Kapiolani Medical Center for Women and Children in Honolulu, Hawaii (2015-2018) were collected. The study was approved by the Western IRB (WIRB Protocol #20151223). Patients meeting the inclusion criteria were identified from pre-admission medical records with pre-pregnancy BMI > 30.0 (maternal obesity) or 18.5-25.0 (normal weight). All participants involved in this study provided written informed consent before the collection of cord blood samples. Pregnant women undergoing elected C-sections at >37 weeks gestation were included, in order to minimize confounding events during the labor. Patient exclusion criteria included pregnant women with preterm rupture of membranes, labor, multiple gestations, pregestational diabetes, hypertensive disorders, cigarette smokers, infection of human immunodeficiency virus or hepatitis B virus, and chronic drug use. Demographic and phenotypic information was recorded, including maternal and paternal age, ethnicity, gestational weight gain, gestational age, parity, and gravidity. For newborns, Apgar scores were documented at 1 minute and 5 minutes post-delivery. The Apgar score serves as a comprehensive assessment of a newborn’s health, with a normal range considered to be between 7 to 10 14.

### Sample preparation and methylation profiling

The baby cord blood sample was collected in the operating room under sterile conditions at the time of the C-section (Pall Medical Cord Blood Collection Kit containing 25ml citrate phosphate dextrose). The umbilical cord was first cleansed with chlorhexidine swabs before cord blood collection. The total volume of collected blood was measured and recorded before aliquoting to conical tubes for centrifugation. The tubes were centrifuged at 200g for 10 min, and plasma was collected. The plasma volume was replaced with 2% FBS/PBS. Negative selection reagents were added to the blood and incubated for 20 min (Miltenyi Biotec, Auburn, CA). The cord blood was diluted with an equal volume of 2% FBS/PBS. A 20ml aliquot of the diluted blood was layered over a density gradient (15ml Lymphoprep) and centrifuged at 1200g for 20 min. The top layer containing an enriched population of stem cells was collected, centrifuged at 300g for 8 min, and then washed in 2% FBS/PBS. Red blood cells were lysed using ammonium chloride (9:1) with incubation on ice for 10 min, washed twice, and then resuspended in 100µl of 2% FBS/PBS. Cells were stained with Lineage FITC and CD34 APC for 45 min on ice, washed twice, and then sorted using the BD FacsAria III. Hematopoietic stem cells (CD34+/CD38-/Lin-) were collected and stored at -80°C until DNA/RNA extraction.

DNA and RNA were extracted simultaneously using the AllPrep DNA/RNA extraction kit (Qiagen). DNA purity and concentration were quantified in Nanodrop 2000 and Picogreen assay. Bisulfite conversion of 500 ng DNA was performed using the EZ DNA Methylation kit (Zymo), followed by sample processing for Infinium HumanMethylation450 bead chips (Illumina) according to the manufacturer’s instructions. Bead chips were analyzed at the Genomics Shared Resource at the University of Hawaii Cancer Center.

### Bulk RNA sequencing

A total of 50 RNA samples were prepared for bulk RNA Sequencing. RNA concentration and RIN score were assayed using Nanodrop 2000 and Agilent Bioanalyzer. A total of 200 ng of high-quality RNA (RIN>7) was subjected to library construction (polyA) and sequenced on HS4000 (2x100) at the Yale Center for Genome Analysis, Connecticut to obtain at least 25M paired reads per sample.

### Methylation data pre-processing

R version 3.6.3 was used for all downstream analyses. Raw intensity data (.idat) were extracted using the ‘ChAMP’ package (version 2.16.2) in R 15–18. Quality control and processing were performed following the ChAMP pipeline. Background controls were subtracted from the data, and raw data that did not pass detection P-value of 0.05 were removed. For each CpG site, the methylation score was initially calculated as the beta value, a fluorescence intensity ratio between 0 and 1. CpG sites whose probes had known underlying SNPs and association with XY chromosomes were removed from analysis due to potential confounding. After BMIQ normalization 19, the batch effect due to non-biological technical variation caused by experiment handling was removed using the ComBat function in the ChAMP package, confirmed by the singular value decomposition (SVD) plot. The M-values for differential analysis were transformed from beta-values using lumi (ver 3.1.4) in R 20–23. A total of 410,765 CpG sites remained for downstream analysis after probe filtering, normalization, and batch removal.

### Source of variation analysis and confounding adjustment

To eliminate potential confounding factors of pre-pregnant maternal obesity among the 13 clinical factors, we conducted a source of variation analysis to identify the clinical factors that significantly contribute to the methylation level variation, as done before 24,25. The variables with F statistics greater than 1 (the error value) were determined as confounders and subjected to confounding adjustment. These factors include the baby’s sex, net weight gain during the pregnancy, maternal age, maternal ethnicity, paternal ethnicity, gravidity, and gestational age. To adjust for confounding, a linear regression model is built using the ‘limma’ package to fit methylation M values of each CpG site, using the confounding factors above. The remaining residuals on the M values were considered to be confounding-adjusted for the subsequent bioinformatics analysis of DNA methylation.

### Bioinformatics analysis of differential methylation (DM)

A moderated t-test from the ‘limma’ R package (version 3.42.2) 26 was used for detecting DM CpG sites between healthy controls and cases with M values. The p-values were adjusted for multiple hypotheses testing using Benjamini-Hochberg FDR. CpG sites with FDR <0.05 were considered statistically significant. To minimize the effect of the gestational age, CpG sites located within the gestational-age-related differentially methylated regions (DMR) were removed. A total of 130 DMRs related to gestational age were identified using linear regression analysis performed with bumphunter 27 across eight public datasets including a total of 248 patients.: GSE31781 28, GSE36829 29, GSE59274 28,30, GSE44667 31, GSE74738 32, GSE49343 33, GSE69502 34, and GSE98224 35,36. The complete list of DMRs was included in **Supplementary Table 1**. Hypermethylation and hypomethylation states were defined by the values of log2 Fold Change (log2FC) of M values in cases compared to controls: hypermethylation if bigger than 0, and hypomethylation if less than 0. Corresponding genes and feature locations of these differential CpG sites were annotated using IlluminaHumanMethylation450kanno.ilmn12.hg19 (ver 0.6.0) 37 .

### KEGG pathway enrichment analysis

‘gometh’ function from R package “missMethyl” (version 1.26.1) 38–41 was used for KEGG pathway enrichment 42–44 with DNA methylation data. DM sites were used for pathway enrichment within five supergroups from KEGG pathways: Metabolism, Genetic Information Processing, Environmental Information Processing, Cellular Processes, and Organismal Systems. Pathways with adjusted p-values less than 0.05 were considered significant. Pathway scores for protein pathways (KEGG: Transcription, Translation, Folding, sorting and degradation) and immune pathways (KEGG: Immune system) were calculated with averaged beta values from the promoter region CpG sites.

### Weighted co-expression network analysis

Firstly, we adjusted all beta values with clinical confounders, then summarized the DM CpG sites at the gene level by averaging the beta values in the promoter regions (those in the TSS200 and TSS1500 promoter regions). Next, we transformed adjusted beta values to adjusted M values for the downstream adjacency matrix construction. We used adjusted M values for the weighted gene co-expression network analysis (WGCNA) with R package ‘WGCNA’ (version 1.70-3) 45,46. The soft threshold for the weighted adjacency matrix with an adjusted R2>0.8 was 7. The topological overlap matrix was constructed for hierarchical clustering. Modules were identified by the dynamic tree-cut algorithm. The networks were exported to Cytoscape with an edge weight greater than 0.03 in each module. The genes with the highest betweenness and degree in the WGCNA network were identified as the hub genes for different modules.

### Protein-protein interaction network analysis

For the protein-protein interaction (PPI) network analysis, DM genes are used as the inputs and were mapped on the STRING database (version 10) 47. Significantly functionally associated protein pairs were identified using PANDA (Preferential Attachment based common Neighbor Distribution derived Associations) (version 0.9.9) 48. KEGG pathways associated with these protein pairs (in terms of genes) were found using PANDA. The bipartite network graph was visualized using Cytoscape (version 3.8.1) 49.

### Stemness score computation

The stemness score was based on Shannon entropy and scaled plasticity, as proposed previously 50. Shannon entropy has been widely applied in developmental biology, particularly in stem cell research 51–53. The formulas are shown below: ![Formula][1]</img>  N is the total number of CpG sites. CpG is represented by the beta value on each CpG probe. The stemness score was calculated for all samples using all remaining 410,765 CpG sites after the preprocessing. A Wilcoxon rank test was performed between the stemness scores of the healthy and maternally obese groups.

### Bulk RNA-seq data processing

The Illumina universal adapter regions of raw RNA-seq data were first trimmed using BBMap (version 38.91) 54. All raw sequences passed the quality control using fastqc (version 0.11.8) 55. The trimmed .fastq files were aligned by STAR (version 2.7.0f) 56 to the human Ensembl genome (Homo_sapiens.GRCh38.dna.primary_assembly.fa) and Ensembl annotation (Homo_sapiens.GRCh38.94.gtf). The gene expression counts were calculated using featureCount57 from Subread (ver 1.6.4) 58.

### Differential expression (DE) of RNA-Seq data

The statistically significant DE genes between healthy controls and maternally obese cases were found adjusting for the same clinical confounders for methylation analysis using the ‘DESeq2’ (version 1.26.0) 59 and ‘limma-voom’ function from ‘limma’ package 26. The p-values were adjusted for multiple hypotheses testing using BH adjustment. Genes with adjusted p-values less than 0.05 were considered statistically significant.

### Correlation analysis between bulk RNA-seq and methylation data

A subset of 47 patients have done both methylation and RNA-seq assyas. Differential expression analysis was done on the bulk RNA-seq data using the ‘limma’ package, with a threshold of BH adjusted p<0.05 to be differential genes. Pearson correlation coefficients (PCC) were calculated between gene expression and methylation beta values from the promoter regions, among the same patients. As mostly a negative correlation between gene expression and DNA methylation in the promoter region is expected, genes with a high negative correlation (PCC<-0.2) were used for pathway enrichment using TOPPFUN 60–62. Top genes of interest were selected with the absolute value Fold Change>1.5 in gene expression and gene-methyl correlation <-0.3 for hyper- and hypo-methylated CpG sites.

### Metabolomics analysis

Metabolomics data were acquired from a previously published study involving 87 patients in the same cohort from three batches (metabolomics workbench study ID ST001114) 63. Targeted metabolites were generated with LC-MS, and untargeted metabolites were generated with GC-MS. After removal of compounds missing in more than 10% samples, a total of 185 metabolites remained, including 10 amino acids (AA), 40 acylcarnitines (AC), 35 acyl/acyl phosphatidylcholines (PC aa), 38 acyl/alkyl phosphatidylcholines (PC ae) and 62 untargeted metabolites. The raw metabolite data were log transformed, standardized, normalized using variance stabilization normalization (VSN), and batch corrected with ComBat function in sva pacakge64. Differential metabolites were identified by limma, with clinical confounders adjustment.

### Multi-omics integration on metabolomics, epigenomics, and transcriptomics

A subset of 42 patients have the matched methylation, gene expression, and metabolomics data. We applied multi-omics integration with Data Integration Analysis for Biomarker discovery using Latent cOmponents (DIABLO) implemented in the mixomics package 65. DIABLO finds the correlated consensus latent variables among different omics in the supervised manner. Top DIABLO features for each omic were selected based on the loading values. We integrated the pathway level methylation, gene, and metabolite interaction using pathview 66.

## Results

### Overview of study design and cohort characteristics

This study aims to investigate the transgenerational effect of pre-pregnancy maternal obesity on offspring. A total of 72 patients who elected to deliver full-term babies through C-sections were recruited from Kapiolani Medical Center for Women and Children in Honolulu, Hawaii from 2016 to 2018. Among them, 38 deliveries are in the healthy control group and 34 are cases with pre-pregnancy maternal obesity. We excluded natural virginal births, to avoid its potential confounding effect on multi-omics profiles. We also carried out stringent recruitment selection criteria, including matching the mothers’ ages as much as possible, as well as similar net gestational weight gain to minimize its confounding effect over maternal pre-pregnancy maternal obesity. The overall study design is shown in **Figure 1**. Briefly, upon collecting the blood samples, umbilical cord blood hematopoietic stem cells (uHSCs) were enriched by FACS sorting with CD34+CD35-LIN-markers (see **Methods**). We extracted DNA and RNA from these uHSCs for Illumina 450k array based DNA methylation and bulk RNA-Seq sequencing respectively. The plasma from these cord blood samples was subject to untargeted metabolomics assays using GC-MS and targeted metabolomics assays using LC–MS 63. Given the rationale that DNA methylation could be the mediator for exerting the transgenerational effect of maternal obesity 67,68, we carried out multi-omics data integration analysis in the DNA methylation-centric manner.

![Figure 1.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2024/07/27/2024.07.27.24310936/F1.medium.gif)

[Figure 1.](http://medrxiv.org/content/early/2024/07/27/2024.07.27.24310936/F1)

Figure 1. Overview of the study design and analysis.
In the preparation step, cord blood plasma samples are collected for metabolome profiling and stem cell sorting. DNA and RNA extraction assays are performed on the enriched stem cells for the methylation and RNA-seq analyses. Downstream analyses are mainly focused on the methylation data. Bulk RNA-seq data were used for validations for methylation discoveries. (Created with BioRender.com)

The demographic details and clinical information of these patients are summarized in **Table 1**. The distributions of the most representative variables are shown in **Figure 2**. Among categorical demographic variables, the distribution of baby sex had no statistical difference between obese and health groups, whereas the ethnicity distributions among mothers and fathers, parity and gravidity are statistically different (P<0.05) between the two groups (**Figure 2A-2E**). Besides maternal pre-pregnancy BMI, other maternal parameters such as maternal age, gestational week, net weight gain and hemoglobin are also not statistically significantly different between the two groups per study design (**Figure 2F-2I**, **Table 1**). While mothers of Asian ethnicity are the majority in the control group, NHPIs account for the majority of the maternal-obese group, revealing the health disparity issue known in the state of Hawaii 69. Moreover, the control group has lower parities and gravidities, compared to the cases. Babies born to obese mothers show significantly higher (P<0.05) body weights compared to the control group, as expected 70. Other parameters including the baby gender, head circumference, body length, and APGAR score at 5 min after birth are not statistically significantly different between case and control groups (**Figure 2J-2M**).

![Figure 2.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2024/07/27/2024.07.27.24310936/F2.medium.gif)

[Figure 2.](http://medrxiv.org/content/early/2024/07/27/2024.07.27.24310936/F2)

Figure 2. Mother and newborns statistics of the multi-ethnic cohort from Hawaii.
(**A-E**) Categorical variables including baby sex, maternal ethnicity, paternal ethnicity, parity and gravidity between control and obese groups are shown in the barplots. P-values using Chi-square test are annotated comparing control and obese groups. **(F-I**) The distributions of maternal age, gestation age, maternal net weight gain during pregnancy, and maternal hemoglobin between control and obese groups are compared. Mean and standard deviation are shown in boxplot. P-values using t-test are annotated. (**J-M**) The distributions of baby weight, baby head circumference, baby length, and APGAR score after 5 minutes of delivery between control and obese groups are compared. Mean and standard deviation are shown in boxplot. P-values using t-test are annotated.

View this table:
[Table 1.](http://medrxiv.org/content/early/2024/07/27/2024.07.27.24310936/T1)

Table 1. Summary statistics of the study cohort.

### Global hypermethylation pattern revealed by CpG level methylation analysis

For scientific rigor, it is critical to adjust for confounding in DNA methylation association analysis 71. Thus we performed the source of variance (SOV) analysis on the beta values of the DNA methylation with respect to physiological and phenotypic information, in order to assess potential confounding factors systematically 24,25,71. As shown in **Figure 3A**, marginal F-statistics in the SOV analysis show that the dominating contribution to DNA methylation variation is maternal pre-pregnancy obesity status, confirming the quality of the study design which aimed to minimize other confounders’ effect. The other minor confounding factors include baby sex, maternal age, maternal ethnicity, net weight gain during pregnancy, paternal ethnicity, gravidity, and gestational age (F-statistics>1). After adjusting these factors by linear regression, all have reduced F-statistics of less than 0.5 (**Figure 3B**) except maternal pre-pregnancy obesity, confirming the success of confounding removal.

![Figure 3.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2024/07/27/2024.07.27.24310936/F3.medium.gif)

[Figure 3.](http://medrxiv.org/content/early/2024/07/27/2024.07.27.24310936/F3)

Figure 3. DNA methylation analysis on uHSCs.
(**A-B**) Source of variance plot before and after confounding adjustment. F-statistics are reported for each clinical factor. F statistics greater than 1 are considered to have confounding effects in addition to the case/control difference due to pre-pregnancy maternal obesity. (**C**) Volcano plot of -log(BH adjusted p-values) against logFC. The cutoff line for adjusted p-value < 0.05 is shown as the red horizontal line. The hyper/hypo threshold is shown as a blue vertical line where logFC=0. Non-significant methylation CpG sites after the differential analysis were shown in gray. Significant CpG sites are colored. After the removal of gestational age-related CpG sites, 10,254 CpG sites are hypermethylated and 5,394 sites are hypomethylated. (**D-E**) Normalized location distribution of differentially methylated CpG sites according to their CpG features in terms of isle regions and gene regions based on the chip annotation. Isle regions include shelf, shore, island, and open-sea. Gene regions include gene body, intergenic region (IGR), TSS200, TSS1500, 5’UTR, 3’UTR, and 1st Exon.

Next, we conducted differential methylation (DE) analysis on the confounding adjusted DNA methylation data (**Methods**). We observed a global hypermethylation pattern in pre-pregnancy obese samples, with 10,254 hypermethylated vs. 5394 hypomethylated CpG sites (**Figure 3C**). The top 20 differentially hypermethylated and hypomethylated CpG sites are reported in **Table 2**, respectively. These CpG sites are related to a wide variety of biological functions, including inflammation (CD69, ADAM12), transcription factors (ZNF222, HMGN4, LHX6, TAF3), proliferation and apoptosis (HDAC4, DHRS4, LRCH3, SAFB2, CRADD, EBF3, PRKAR1B). Some top DM CpG sites are associated with obesity directly, including HDAC4 72 and PLEC173.

View this table:
[Table 2.](http://medrxiv.org/content/early/2024/07/27/2024.07.27.24310936/T2)

Table 2. Top 20 hypermethylated CpG sites and top 20 hypomethylated CpG sites.

We further examined the distributions of these differentially methylated sites, relative to the CpG island regions and promoter proximity. (**Figure 3D-E**). A big fraction (38.7%) of the DM sites are located in CpG islands 74,75, significantly higher than that from the Illumina 450K annotation (P<2E-16). CpG islands are more frequent in the hypermethylated sites (41.3%) as compared to these in the hypomethylated sites (33.6%), consistent with the global hypermethylation pattern. Relative to gene localization, DM sites are most frequent (32.1%) in the promoter regions (including 16.8% and 15.3% in TSS200 and TSS1500 respectively) as expected.

### Functional analyses reveal the association between maternal obesity and cell cycle, immune function and metabolic changes in the cord blood of offspring

To investigate the biological functions related to the epigenome alternation, we conducted systematic analysis of DM sites employing multiple methods: KEGG pathway enrichment analysis, Weighted Gene Co-expression Network Analysis (WGCNA), and Protein-Protein Interaction (PPI) network analysis.

KEGG pathway enrichment analysis on hypermethylated CpG sites identified five significant pathways with FDR<0.05 (**Figure 4A)**, including the cell cycle, ribosome, nucleocytoplasmic transport, ribosome biogenesis in eukaryotes, and mTOR signaling pathway. Cell cycle, ribosome, and nucleocytoplasmic transport pathways are essential to normal cell functioning. mTOR signaling pathway coordinates the nutrient-mediated metabolism, immune responses and cell cycle progression, and dysregulation of mTOR could lead to various diseases such as cancer and obesity 76. There was no significantly enriched pathway emerging from hypomethylated CpG sites. The maternally obese group shows significantly higher methylation levels in KEGG protein synthesis and immune system pathway collections compared to the control group, indicating repression in immune response as well as translation and protein synthesis (**Figure 4B-C**). Similarly, we further explored the differential potential, or stemness, of umbilical cord Hematopoietic Stem Cells (uHSCs). We first confirmed the homogeneity of uHSCs by single-cell RNA sequencing UMAP plot (**Supplementary Figure 2**). We calculated the cell stemness scores using the DNA methylation beta values similar to others 77. uHSCs derived from the maternally obese group exhibit significantly elevated stemness scores (P<0.01) in comparison to the control group **(Figure 4D)**, confirming the results in KEGG pathway analysis.

![Figure 4.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2024/07/27/2024.07.27.24310936/F4.medium.gif)

[Figure 4.](http://medrxiv.org/content/early/2024/07/27/2024.07.27.24310936/F4)

Figure 4. Pathway and network analysis.
(**A**) KEGG pathway enrichment for hypermethylated CpG sites from promotor region. Enriched KEGG pathway names, adjusted p-values (-log10 transformed), and the size of enriched gene list are reported for CpG sites from TSS200+TSS1500 regions. The red dotted line shows the threshold cutoff for FDR at -log10(0.05). (**B-C**) Boxplots of averaged beta values for KEGG protein pathway collection and immune pathway collection with Wilcoxon P-values. (**D**) Violin plots of cell entropy scores between control and obese groups with Wilcoxon P-values. (**E-F**) WGCNA network analysis results. WGCNA modules are shown for both the control and the obese group. The top two modules with largest degrees are turquoise and brown modules. Each node represents a gene. Genes co-expressed in each module are annotated. (**G**)Protein-protein interaction (PPI) network. Bipartite graphs represent enriched KEGG pathways and associated genes with significant PPIs. Red nodes represent genes with hypermethylated CpG sites. Blue nodes represent genes with hypomethylated CpG sites. Yellow nodes represented the enriched KEGG pathways. Number of inter-pathway PPIs are annotated in the rectangular boxes.

Next, we applied WGCNA to cluster co-regulation of gene-level methylation, by averaging CpG sites to affiliated genes (see **Methods**). Five co-expression modules are identified, using the M-values adjusted for clinical confounders (**Supplementary Figure 1A)**, and all modules show positive correlations with maternal obesity except one. The largest turquoise module (**Figure 4E**) is related to cell cycle, protein synthesis, and transport and vesicle trafficking pathways through pathway enrichment analysis. Some hub genes in this module are identified, including INTU, ANAPC7, and AGBL5. These genes were reported essential for maintaining cell polarity (INTU)78, proliferation (ANAPC7) 79 and glycemic control (AGBL5) 80. The brown module (**Figure 4F**) is enriched with immune response pathways, in which TLR6 is identified as a hub gene. The other yellow module is related to ion homeostasis, and the gray module is related to the p53 pathway, apoptosis, cell senescence, and ER stress (**Supplementary Figure 1B)**. The only negatively correlated blue module is associated with axon guidance and VEGF signaling pathway (**Supplementary Figure 1B**).

Furthermore, we examined the PPI network, using the gene-level DNA methylation as surrogates (**Figure 4G**). The PPI analysis identifies 14 unique pathways (FDR < 0.05) predominantly associated with hypermethylated CpG sites in the TSS200 and TSS1500 regions. The top five largest pathways included ribosome, proteasome, cell cycle, axon guidance, RNA polymerase, and neuroactive ligand-receptor interaction. Taken all three types of systematic analyses together, cell cycle, immune function and protein synthesis are ubiquitously highlighted, suggesting that these biological functions in cord blood stem cells are negatively impacted by maternal obesity.

### Multi-omics analysis reveals disruptions in cell cycle and metabolic pathways

To systematically investigate the epigenetic, transcriptomic, and metabolomic alterations induced by maternal obesity, we performed multi-omics integration analysis on this cohort. We employed DIABLO, a supervised integration method that extracts features associated with maternal obesity, based on the correlations in the embedding space 65. **Figure 5A-C** shows that methylation data provide the clearest separation between obese and control groups, confirming the value of the earlier DNA methylation centered analysis.

![Figure 5.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2024/07/27/2024.07.27.24310936/F5.medium.gif)

[Figure 5.](http://medrxiv.org/content/early/2024/07/27/2024.07.27.24310936/F5)

Figure 5. Multi-omics integration analysis
(**A-C**) Omics-specific sample plots from DIABLO showing the separation of obese and control samples in methylation data, gene expression data, and metabolomics data respectively. (**D-F**) Importance plot of top 25 features in methylation, gene expression and metabolomics modalities with the highest loadings extracted from the embedding space. The color represents the condition which features contribute the most.

The top 25 features from each omic with the highest feature weights (loadings) following integrated canonical correlation analysis are demonstrated in **Figure 5D-F**. The methylation features with the highest weights related to maternal obesity include CpG sites involved in cell-cycle control, glucose metabolism, and adipogenesis (FOXO181), DNA repair (LIG3, SMUG1), erythropoietin pathway and differentiation (EPO, CSNK2A1, CSF1), which are hypermethylated in the obese group. Hypomethylation of LEP (encoding leptin) was also observed as a top feature, aligning with prior findings that maternal obesity is associated with elevated maternal leptin levels, a known marker of adipose tissue 82. These featured CpG sites indicate repression in fat metabolism and DNA repair and reduced differentiation potential. In the transcriptomic space, many genes related to mRNA splicing (SRRM1, IGF2BP1, IGF2BP2, CNOT4) have increased expression levels due to maternal obesity. Among the metabolite features, essential sugars (glucose, xylose), poly-unsatuated fatty acids (oleic acids, DHA, arachidonic acid), and phosphatidylcholine (PCs) are mostly decreased in the obese group; whereas most acylcarnitines (C) are elevated. The metabolic changes show an overall accumulation of saturated fatty acid, but repression on fat breakdown, glucose and unsaturated fatty acid generation. As poly-unsatuated fatty acids (eg. arachidonic acid) have important anti-inflammatory effects, the results indicate a pro-inflammatory environment in offspring born of pre-pregnant obese mothers.

## Discussion

Maternal obesity is one of the most urgent health concerns worldwide. Pre-pregnancy maternal obesity could cause various pregnancy-related complications and predispose offspring to cardiometabolic complications and chronic diseases in the long term 9. Multiple cross-continental large cohort meta-analyses have shown that maternal obesity is directly associated with offspring’s risk of obesity, coronary heart disease, insulin resistance, and adverse neurodevelopmental outcome based on longitudinal observational studies 9,83,84. To directly ping-point the molecular level changes in offspring by maternal pre-pregnancy obesity, we used cord blood stem cells as the studying material, which serve as a great surrogate revealing newborn’s metabolism and immune system changes at the time of birth 85. The current study expands previous effects and investigates the direct impact of maternal obesity on uHSCs programming, the progenitor of immune cell population, using a multi-omics (epigenetic, gene expression, and metabolite) analysis approach from a unique multi-ethnic cohort.

Centered around methylation changes, three complimentary functional analysis approaches (KEGG, WGCNA, and PPI) consistently demonstrated that maternal obesity impacts multiple bilogical functions including hypermethylation in promoters of genes involved in cell cycle, ribosome biogenesis, and mTOR signaling pathways. Moreover, mTOR signaling pathway also plays a crucial role in metabolism and cell cycle regulation, disruption in this pathway leads to insulin resistance and long-term diseases 86. We observed a significant increase in stemness scores among uHSCs affected by maternal obesity, aligning with expected downregulation in the cell cycle gene expression due to observed hypermethylation in the promoters of these genes.

Higher stembess scores indicate enhanced quiescence, shift the balance between stem cell maintenance and differentiation towards the former. Unlike adult HSCs, fetal/neonatal HSCs typically exhibit higher proliferation and self-renewal capabilities, crucial for blood cell regeneration and innate immune system development 87. Our findings provide strong epigenetic evidence that maternal obesity compromises the maturation processes in neonatal uHSCs, which may predispose newborns to immunological disorders.

The subsequent multi-omics integration analysis expanded conclusions from methylation analysis to additional metabolomics readouts that are also linked to biological functions eg. cell cycle and inflammatory pathway. We thus propose the conceptual model to illustrate the effect of maternal pre-pregnancy obesity (**Figure 6**). Maternal obesity leads to nutrient deficiency with lower levels of essential amino acids and fatty acids in the newborn blood and disrupts the lipid metabolism homeosis in offspring. These metabolite changes further induce cell membrane instability and repress cell cycle progression, cell proliferation88, enhancing the dysregulation of these functions preexisting at the methylation level. Lipid dysregulation may also enhance the pro-inflammatory environment, which in turn induces complications in offspring later in life, such as cardiovascular diseases. Such proposed model is also consistent with and further strengthens previous studies at the metabolomics or epigenome levels. For example, previous metabolomics studies of cord blood showed metabolic derangement predisposes newborns to cardiometabolic and endocrine diseases, and disrupt the normal hormone function and neonatal adiposity 82,89. Previous epigenome-wide association study (EWAS) with cord blood found a strong association between DNA methylation pattern and postnatal BMI trajectory until adolescent 90.

![Figure 6.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2024/07/27/2024.07.27.24310936/F6.medium.gif)

[Figure 6.](http://medrxiv.org/content/early/2024/07/27/2024.07.27.24310936/F6)

Figure 6. A proposed model of maternal obesity’s impact on neonatal development.

In summary, this newborn study demonstrates the direct impact of maternal pre-pregnancy obesity and on newborn blood at the multi-omics level, which includes increased cell cycle arrest, impairment in the uHSC differentiation capacity, more inflammation and disruption in lipid metabolism.

### Declaration of generative AI and AI-assisted technologies in the writing process

During the preparation of this work the author(s) used GPT-4.0 in order to improve the readability. After using this tool/service, the author(s) reviewed and edited the content thoroughly and take(s) full responsibility for the content of the publication.

#### Data availability statement

Data generated in this study have been submitted and will be available through the National Institutes of Health Gene Expression Omnibus (accession number pending). Code to produce the analyses in this manuscript are available through GitHub ([https://github.com/lanagarmire/](https://github.com/lanagarmire/))

## Supporting information

supplemental figures and tables [[supplements/310936_file09.pdf]](pending:yes)

## Data Availability

All data produced in the present work are contained in the manuscript

## Author contributions

LG envisioned this project, obtained the funding, supervised the study and revised the manuscript. YD performed the data analysis, generated the figures, and wrote the initial manuscript. RS consented patients and obtained the samples from the hospital, with coordination from PB. DT and SW coordinated with the patient recruitment and study. PB coordinated all the multi-omics assays. PB, CL, and FA designed the DNA methylation assays. AG performed FACS sorting of cord blood cells. ALJ performed the Illumina Meth 450 assay, MT supervised the Genomics Shared Resource analyses and provided a critical review of the manuscript. All authors have read the manuscript.

## Conflicts of interest

None

## Acknowledgments

This research was supported by grants R01 LM012373 and LM012907 awarded by NLM, and R01 HD084633 awarded by NICHD to L.X. Garmire, as well as in part by the NCI Cancer Center Support Grant (CCSG) number P30 CA071789 awarded to Genomics and Bioinformatics Shared Resource (RRID:SCR_019085). This research was supported in part by training funding provided by the NIH grant T32 GM141746 and Advanced Proteogenomics of Cancer (T32 CA140044).

## Abbreviation

BH
:   Benjamini-Hochberg
BMI
:   Body mass index
C
:   Acylcarnitines
DE
:   Differential expression
DIABLO
:   Data Integration Analysis for Biomarker discovery using Latent cOmponents
DM
:   Differential methylation
DMR
:   Differentially methylated regions
DOHaD
:   Developmental Origins of Health and Disease
EWAS
:   Epigenome-wide association studies
FC
:   Fold Change
FDR
:   False positive results
KEGG
:   Kyoto Encyclopedia of Genes and Genomes
LASSO
:   Least absolute shrinkage and selection operator
LDA
:   Linear Discriminant Analysis
LOG
:   Logistic regression
PAM
:   Partition Around Medoids
PANDA
:   Preferential Attachment-based common Neighbor Distribution derived Associations
PC aa
:   Diacyl phosphatidylcholines
PC ae
:   Acyl-alkylphosphatidylcholines
PCC
:   Pearson correlation coefficients
PPI
:   Protein-Protein Interaction
RBF
:   Radial Basis Function
RF
:   Random Forest
RPART
:   Recursive Partitioning and Regression Trees
SOV
:   Source of variance
SVD
:   Singular value decomposition
SVM
:   Supportive Vector Machine
TSS
:   Transcription start site
uHSCs
:   Umbilical cord blood hematopoietic stem cells
UMAP
:   Uniform Manifold Approximation and Projection
VLCAD
:   Very long-chain acyl-CoA dehydrogenase
WGCNA
:   Weighted Gene Co-expression Network Analysis

*   Received July 27, 2024.
*   Revision received July 27, 2024.
*   Accepted July 27, 2024.


*   © 2024, Posted by Cold Spring Harbor Laboratory

This pre-print is available under a Creative Commons License (Attribution-NonCommercial-NoDerivs 4.0 International), CC BY-NC-ND 4.0, as described at [http://creativecommons.org/licenses/by-nc-nd/4.0/](http://creativecommons.org/licenses/by-nc-nd/4.0/)

## References

1.  1.Leddy, M.A., Power, M.L., and Schulkin, J. (2008). The impact of maternal obesity on maternal and fetal health. Rev. Obstet. Gynecol. 1, 170–178.
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=19173021&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F27%2F2024.07.27.24310936.atom) 

2.  2.Hjalgrim, L.L., Westergaard, T., Rostgaard, K., Schmiegelow, K., Melbye, M., Hjalgrim, H., and Engels, E.A. (2003). Birth weight as a risk factor for childhood leukemia: a meta-analysis of 18 epidemiologic studies. Am. J. Epidemiol. 158, 724–735.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/aje/kwg210&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=14561661&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F27%2F2024.07.27.24310936.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000185897700002&link_type=ISI) 

3.  3.Harder, T., Plagemann, A., and Harder, A. (2010). Birth weight and risk of neuroblastoma: a meta-analysis. Int. J. Epidemiol. 39, 746–756.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/ije/dyq040&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=20236985&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F27%2F2024.07.27.24310936.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000278438500013&link_type=ISI) 

4.  4.Cnattingius, S., Lundberg, F., Sandin, S., Grönberg, H., and Iliadou, A. (2009). Birth characteristics and risk of prostate cancer: the contribution of genetic factors. Cancer Epidemiol. Biomarkers Prev. 18, 2422–2426.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoiY2VicCI7czo1OiJyZXNpZCI7czo5OiIxOC85LzI0MjIiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyNC8wNy8yNy8yMDI0LjA3LjI3LjI0MzEwOTM2LmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 

5.  5.Eriksson, M., Wedel, H., Wallander, M.-A., Krakau, I., Hugosson, J., Carlsson, S., and Svärdsudd, K. (2007). The impact of birth weight on prostate cancer incidence and mortality in a population-based study of men born in 1913 and followed up from 50 to 85 years of age. Prostate 67, 1247–1254.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/pros.20428&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=17570499&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F27%2F2024.07.27.24310936.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000248537500012&link_type=ISI) 

6.  6.Michos, A., Xue, F., and Michels, K.B. (2007). Birth weight and the risk of testicular cancer: a meta-analysis. Int. J. Cancer 121, 1123–1131.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/ijc.22771&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=17471572&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F27%2F2024.07.27.24310936.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000248242500026&link_type=ISI) 

7.  7.Silva, I. dos S., De Stavola, B., McCormack, V., and Collaborative Group on Pre-Natal Risk Factors and Subsequent Risk of Breast Cancer (2008). Birth size and breast cancer risk: re-analysis of individual participant data from 32 studies. PLoS Med. 5, e193.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pmed.0050193&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=18828667&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F27%2F2024.07.27.24310936.atom) 

8.  8.8. Van Cleave, J., Gortmaker, S.L., and Perrin, J.M. (2010). Dynamics of obesity and chronic health conditions among children and youth. JAMA 303, 623–630.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1001/jama.2010.104&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=20159870&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F27%2F2024.07.27.24310936.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000274518600020&link_type=ISI) 

9.  9.Godfrey, K.M., Reynolds, R.M., Prescott, S.L., Nyirenda, M., Jaddoe, V.W., Eriksson, J.G., and Broekman, B.F. (2017). Influence of maternal obesity on the long-term health of offspring. The lancet. Diabetes & endocrinology 5. doi:10.1016/S2213-8587(16)30107-3.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S2213-8587(16)30107-3&link_type=DOI) 

10. 10.Barker, D.J. (2007). The origins of the developmental origins theory. J. Intern. Med. 261. doi:10.1111/j.1365-2796.2007.01809.x.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1111/j.1365-2796.2007.01809.x&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=17444880&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F27%2F2024.07.27.24310936.atom) 

11. 11.Barker, D.J. (2000). In utero programming of cardiovascular disease. Theriogenology 53. doi:10.1016/s0093-691x(99)00258-7.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/s0093-691x(99)00258-7&link_type=DOI) 

12. 12.Qiu, L., Low, H.P., Chang, C.-I., Strohsnitter, W.C., Anderson, M., Edmiston, K., Adami, H.-O., Ekbom, A., Hall, P., Lagiou, P., et al. (2012). Novel measurements of mammary stem cells in human umbilical cord blood as prospective predictors of breast cancer susceptibility in later life. Ann. Oncol. 23, 245–250.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/annonc/mdr153&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=21515664&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F27%2F2024.07.27.24310936.atom) 

13. 13.Strohsnitter, W.C., Savarese, T.M., Low, H.P., Chelmow, D.P., Lagiou, P., Lambe, M., Edmiston, K., Liu, Q., Baik, I., Noller, K.L., et al. (2008). Correlation of umbilical cord blood haematopoietic stem and progenitor cell levels with birth weight: implications for a prenatal influence on cancer risk. Br. J. Cancer 98, 660–663.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/sj.bjc.6604183&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=18256588&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F27%2F2024.07.27.24310936.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000252933500022&link_type=ISI) 

14. 14.Apgar, V. (1953). A Proposal for a New Method of Evaluation of the Newborn Infant. Anesthesia & Analgesia 32, 260.
    
    

15. 15.Morris, T.J., Butcher, L.M., Feber, A., Teschendorff, A.E., Chakravarthy, A.R., Wojdacz, T.K., and Beck, S. (2014). ChAMP: 450k Chip Analysis Methylation Pipeline. Bioinformatics 30, 428–430.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/bioinformatics/btt684&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24336642&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F27%2F2024.07.27.24310936.atom) 

16. 16.Aryee, M.J., Jaffe, A.E., Corrada-Bravo, H., Ladd-Acosta, C., Feinberg, A.P., Hansen, K.D., and Irizarry, R.A. (2014). Minfi: a flexible and comprehensive Bioconductor package for the analysis of Infinium DNA methylation microarrays. Bioinformatics 30, 1363–1369.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/bioinformatics/btu049&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24478339&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F27%2F2024.07.27.24310936.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000336530000004&link_type=ISI) 

17. 17.Fortin, J.-P., Triche, T.J., Jr., and Hansen, K.D. (2017). Preprocessing, normalization and integration of the Illumina HumanMethylationEPIC array with minfi. Bioinformatics 33, 558–560.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/bioinformatics/btw691&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F27%2F2024.07.27.24310936.atom) 

18. 18.Zhou, W., Laird, P.W., and Shen, H. (2017). Comprehensive characterization, annotation and innovative use of Infinium DNA methylation BeadChip probes. Nucleic Acids Res. 45, e22.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/nar/gkw967&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F27%2F2024.07.27.24310936.atom) 

19. 19.Teschendorff, A.E., Marabita, F., Lechner, M., Bartlett, T., Tegner, J., Gomez-Cabrero, D., and Beck, S. (2013). A beta-mixture quantile normalization method for correcting probe design bias in Illumina Infinium 450 k DNA methylation data. Bioinformatics 29, 189–196.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/bioinformatics/bts680&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23175756&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F27%2F2024.07.27.24310936.atom) 

20. 20.Du, P., Kibbe, W.A., and Lin, S.M. (2007). nuID: a universal naming scheme of oligonucleotides for illumina, affymetrix, and other microarrays. Biol. Direct 2, 16.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/1745-6150-2-16&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=17540033&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F27%2F2024.07.27.24310936.atom) 

21. 21.Lin, S.M., Du, P., Huber, W., and Kibbe, W.A. (2008). Model-based variance-stabilizing transformation for Illumina microarray data. Nucleic Acids Res. 36, e11.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/nar/gkm1075&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=18178591&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F27%2F2024.07.27.24310936.atom) 

22. 22.Du, P., Kibbe, W.A., and Lin, S.M. (2008). lumi: a pipeline for processing Illumina microarray. Bioinformatics 24, 1547–1548.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/bioinformatics/btn224&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=18467348&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F27%2F2024.07.27.24310936.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000257169700015&link_type=ISI) 

23. 23.Du, P., Zhang, X., Huang, C.-C., Jafari, N., Kibbe, W.A., Hou, L., and Lin, S.M. (2010). Comparison of Beta-value and M-value methods for quantifying methylation levels by microarray analysis. BMC Bioinformatics 11, 587.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/1471-2105-11-587&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=21118553&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F27%2F2024.07.27.24310936.atom) 

24. 24.He, B., Liu, Y., Maurya, M.R., Benny, P., Lassiter, C., Li, H., Subramaniam, S., and Garmire, L.X. (2021). The maternal blood lipidome is indicative of the pathogenesis of severe preeclampsia. J. Lipid Res. 62. doi:10.1016/j.jlr.2021.100118.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.jlr.2021.100118&link_type=DOI) 

25. 25.Chen, Y., He, B., Liu, Y., Aung, M.T., Rosario-Pabón, Z., Vélez-Vega, C.M., Alshawabkeh, A., Cordero, J.F., Meeker, J.D., and Garmire, L.X. (2022). Maternal plasma lipids are involved in the pathogenesis of preterm birth. Gigascience 11. doi:10.1093/gigascience/giac004.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/gigascience/giac004&link_type=DOI) 

26. 26.Ritchie, M.E., Phipson, B., Wu, D., Hu, Y., Law, C.W., Shi, W., and Smyth, G.K. (2015). limma powers differential expression analyses for RNA-sequencing and microarray studies. Nucleic Acids Res. 43, e47.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/nar/gkv007&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25605792&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F27%2F2024.07.27.24310936.atom) 

27. 27.Jaffe, A.E., Murakami, P., Lee, H., Leek, J.T., Fallin, M.D., Feinberg, A.P., and Irizarry, R.A. (2012). Bump hunting to identify differentially methylated regions in epigenetic epidemiology studies. Int. J. Epidemiol. 41, 200–209.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/ije/dyr238&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22422453&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F27%2F2024.07.27.24310936.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000302026800023&link_type=ISI) 

28. 28.Novakovic, B., Yuen, R.K., Gordon, L., Penaherrera, M.S., Sharkey, A., Moffett, A., Craig, J.M., Robinson, W.P., and Saffery, R. (2011). Evidence for widespread changes in promoter methylation profile in human placenta in response to increasing gestational age and environmental/stochastic factors. BMC Genomics 12. doi:10.1186/1471-2164-12-529.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/1471-2164-12-529&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22032438&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F27%2F2024.07.27.24310936.atom) 

29. 29.GEO Accession viewer [https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE36829](https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE36829).
    
    

30. 30.Chu, T., Bunce, K., Shaw, P., Shridhar, V., Althouse, A., Hubel, C., and Peters, D. (2014). Comprehensive analysis of preeclampsia-associated DNA methylation in the placenta. PLoS One 9. doi:10.1371/journal.pone.0107318.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pone.0107318&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25247495&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F27%2F2024.07.27.24310936.atom) 

31. 31.Blair, J.D., Yuen, R.K., Lim, B.K., McFadden, D.E., von Dadelszen, P., and Robinson, W.P. (2013). Widespread DNA hypomethylation at gene enhancer regions in placentas associated with early-onset pre-eclampsia. Mol. Hum. Reprod. 19. doi:10.1093/molehr/gat044.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/molehr/gat044&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23770704&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F27%2F2024.07.27.24310936.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000324828900008&link_type=ISI) 

32. 32.Hanna, C.W., Peñaherrera, M.S., Saadeh, H., Andrews, S., McFadden, D.E., Kelsey, G., and Robinson, W.P. (2016). Pervasive polymorphic imprinted methylation in the human placenta. Genome Res. 26. doi:10.1101/gr.196139.115.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NjoiZ2Vub21lIjtzOjU6InJlc2lkIjtzOjg6IjI2LzYvNzU2IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjQvMDcvMjcvMjAyNC4wNy4yNy4yNDMxMDkzNi5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

33. 33.Blair, J.D., Langlois, S., McFadden, D.E., and Robinson, W.P. (2014). Overlapping DNA methylation profile between placentas with trisomy 16 and early-onset preeclampsia. Placenta 35. doi:10.1016/j.placenta.2014.01.001.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.placenta.2014.01.001&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24462402&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F27%2F2024.07.27.24310936.atom) 

34. 34.34. Price, E.M., Peñaherrera, M.S., Portales-Casamar, E., Pavlidis, P., Van Allen, M.I., McFadden, D.E., and Robinson, W.P. (2016). Profiling placental and fetal DNA methylation in human neural tube defects. Epigenetics Chromatin 9. doi:10.1186/s13072-016-0054-8.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/s13072-016-0054-8&link_type=DOI) 

35. 35.Leavey, K., Wilson, S.L., Bainbridge, S.A., Robinson, W.P., and Cox, B.J. (2018). Epigenetic regulation of placental gene expression in transcriptional subtypes of preeclampsia. Clin. Epigenetics 10. doi:10.1186/s13148-018-0463-6.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/s13148-018-0463-6&link_type=DOI) 

36. 36.Wilson, S.L., Leavey, K., Cox, B.J., and Robinson, W.P. (2018). Mining DNA methylation alterations towards a classification of placental pathologies. Hum. Mol. Genet. 27. doi:10.1093/hmg/ddx391.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/hmg/ddx391&link_type=DOI) 

37. 37.Hansen, K.D. IlluminaHumanMethylation450kanno. ilmn12. hg19: annotation for Illumina’s 450k methylation arrays. R package version 0.6. 0.
    
    

38. 38.Maksimovic, J., Gordon, L., and Oshlack, A. (2012). SWAN: Subset-quantile within array normalization for illumina infinium HumanMethylation450 BeadChips. Genome Biol. 13, R44.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/gb-2012-13-6-r44&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22703947&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F27%2F2024.07.27.24310936.atom) 

39. 39.Phipson, B., and Oshlack, A. (2014). DiffVar: a new method for detecting differential variability with application to methylation in cancer and aging. Genome Biol. 15, 465.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/s13059-014-0465-4&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25245051&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F27%2F2024.07.27.24310936.atom) 

40. 40.Maksimovic, J., Gagnon-Bartsch, J.A., Speed, T.P., and Oshlack, A. (2015). Removing unwanted variation in a differential methylation analysis of Illumina HumanMethylation450 array data. Nucleic Acids Res. 43, e106.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/nar/gkv526&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25990733&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F27%2F2024.07.27.24310936.atom) 

41. 41.Phipson, B., Maksimovic, J., and Oshlack, A. (2016). missMethyl: an R package for analyzing data from Illumina’s HumanMethylation450 platform. Bioinformatics 32, 286–288.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/bioinformatics/btv560&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26424855&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F27%2F2024.07.27.24310936.atom) 

42. 42.Kanehisa, M., and Goto, S. (2000). KEGG: kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 28, 27–30.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/nar/28.1.27&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=10592173&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F27%2F2024.07.27.24310936.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000084896300007&link_type=ISI) 

43. 43.Kanehisa, M. (2019). Toward understanding the origin and evolution of cellular organisms. Protein Sci. 28, 1947–1951.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/pro.3715&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=31441146&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F27%2F2024.07.27.24310936.atom) 

44. 44.Kanehisa, M., Furumichi, M., Sato, Y., Kawashima, M., and Ishiguro-Watanabe, M. (2022). KEGG for taxonomy-based analysis of pathways and genomes. Nucleic Acids Res. doi:10.1093/nar/gkac963.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/nar/gkac963&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=36300620&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F27%2F2024.07.27.24310936.atom) 

45. 45.Langfelder, P., and Horvath, S. (2008). WGCNA: an R package for weighted correlation network analysis. BMC Bioinformatics 9, 559.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/1471-2105-9-559&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=19114008&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F27%2F2024.07.27.24310936.atom) 

46. 46.Langfelder, P., and Horvath, S. (2012). Fast R Functions for Robust Correlations and Hierarchical Clustering. J. Stat. Softw. 46. doi:10.18637/jss.v046.i11.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.18637/jss.v046.i11&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23050260&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F27%2F2024.07.27.24310936.atom) 

47. 47.Szklarczyk, D., Gable, A.L., Nastou, K.C., Lyon, D., Kirsch, R., Pyysalo, S., Doncheva, N.T., Legeay, M., Fang, T., Bork, P., et al. (2021). The STRING database in 2021: customizable protein-protein networks, and functional characterization of user-uploaded gene/measurement sets. Nucleic Acids Res. 49, D605–D612.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/nar/gkaa1074&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=33237311&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F27%2F2024.07.27.24310936.atom) 

48. 48.Li, H., Tong, P., Gallegos, J., Dimmer, E., Cai, G., Molldrem, J.J., and Liang, S. (2015). PAND: A Distribution to Identify Functional Linkage from Networks with Preferential Attachment Property. PLoS One 10, e0127968.
    
    

49. 49.Shannon, P., Markiel, A., Ozier, O., Baliga, N.S., Wang, J.T., Ramage, D., Amin, N., Schwikowski, B., and Ideker, T. (2003). Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. 13, 2498–2504.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NjoiZ2Vub21lIjtzOjU6InJlc2lkIjtzOjEwOiIxMy8xMS8yNDk4IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjQvMDcvMjcvMjAyNC4wNy4yNy4yNDMxMDkzNi5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

50. 50.He, B., and Garmire, L.X. (2021). ASGARD: A Single-cell Guided pipeline to Aid Repurposing of Drugs. ArXiv.
    
    

51. 51.MacArthur, B.D., and Lemischka, I.R. (2013). Statistical mechanics of pluripotency. Cell 154, 484–489.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.cell.2013.07.024&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23911316&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F27%2F2024.07.27.24310936.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000322629900004&link_type=ISI) 

52. 52.Martínez, O., and Reyes-Valdés, M.H. (2008). Defining diversity, specialization, and gene specificity in transcriptomes through information theory. Proc. Natl. Acad. Sci. U. S. A. 105, 9709–9714.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoicG5hcyI7czo1OiJyZXNpZCI7czoxMToiMTA1LzI4Lzk3MDkiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyNC8wNy8yNy8yMDI0LjA3LjI3LjI0MzEwOTM2LmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 

53. 53.Kannan, S., Farid, M., Lin, B.L., Miyamoto, M., and Kwon, C. (2021). Transcriptomic entropy benchmarks stem cell-derived cardiomyocyte maturation against endogenous tissue at single cell level. PLoS Comput. Biol. 17, e1009305.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pcbi.1009305&link_type=DOI) 

54. 54.Bushnell, B. (2014). BBMap: A Fast, Accurate, Splice-Aware Aligner.
    
    

55. 55.Andrews, S. FastQC: a quality control tool for high throughput sequence data. Available online. Retrieved May.
    
    

56. 56.Dobin, A., Davis, C.A., Schlesinger, F., Drenkow, J., Zaleski, C., Jha, S., Batut, P., Chaisson, M., and Gingeras, T.R. (2013). STAR: ultrafast universal RNA-seq aligner. Bioinformatics 29, 15–21.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/bioinformatics/bts635&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23104886&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F27%2F2024.07.27.24310936.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000312654600003&link_type=ISI) 

57. 57.Liao, Y., Smyth, G.K., and Shi, W. (2014). featureCounts: an efficient general purpose program for assigning sequence reads to genomic features. Bioinformatics 30, 923–930.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/bioinformatics/btt656&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24227677&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F27%2F2024.07.27.24310936.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000334078300005&link_type=ISI) 

58. 58.Liao, Y., Smyth, G.K., and Shi, W. (2013). The Subread aligner: fast, accurate and scalable read mapping by seed-and-vote. Nucleic Acids Res. 41, e108.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/nar/gkt214&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23558742&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F27%2F2024.07.27.24310936.atom) 

59. 59.Love, M.I., Huber, W., and Anders, S. (2014). Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 15, 550.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/s13059-014-0550-8&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25516281&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F27%2F2024.07.27.24310936.atom) 

60. 60.Chen, J., Xu, H., Aronow, B.J., and Jegga, A.G. (2007). Improved human disease candidate gene prioritization using mouse phenotype. BMC Bioinformatics 8, 392.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/1471-2105-8-392&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=17939863&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F27%2F2024.07.27.24310936.atom) 

61. 61.Chen, J., Aronow, B.J., and Jegga, A.G. (2009). Disease candidate gene identification and prioritization using protein interaction networks. BMC Bioinformatics 10, 73.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/1471-2105-10-73&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=19245720&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F27%2F2024.07.27.24310936.atom) 

62. 62.Chen, J., Bardes, E.E., Aronow, B.J., and Jegga, A.G. (2009). ToppGene Suite for gene list enrichment analysis and candidate gene prioritization. Nucleic Acids Res. 37, W305–W311.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/nar/gkp427&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=19465376&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F27%2F2024.07.27.24310936.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000267889100054&link_type=ISI) 

63. 63.Schlueter, R.J., Al-Akwaa, F.M., Benny, P.A., Gurary, A., Xie, G., Jia, W., Chun, S.J., Chern, I., and Garmire, L.X. (2020). Prepregnant Obesity of Mothers in a Multiethnic Cohort Is Associated with Cord Blood Metabolomic Changes in Offspring. J. Proteome Res. 19, 1361–1374.
    
    

64. 64.Johnson, W.E., Li, C., and Rabinovic, A. (2006). Adjusting batch effects in microarray expression data using empirical Bayes methods. Biostatistics 8, 118–127.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/biostatistics/kxj037&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=16632515&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F27%2F2024.07.27.24310936.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000242715400008&link_type=ISI) 

65. 65.Singh, A., Shannon, C.P., Gautier, B., Rohart, F., Vacher, M., Tebbutt, S.J., and Cao, K.-A.L. (2019). DIABLO: an integrative approach for identifying key molecular drivers from multi-omics assays. Bioinformatics 35, 3055.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/bioinformatics/bty1054&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F27%2F2024.07.27.24310936.atom) 

66. 66.Luo, W., and Brouwer, C. (2013). Pathview: an R/Bioconductor package for pathway-based data integration and visualization. Bioinformatics 29, 1830–1831.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/bioinformatics/btt285&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23740750&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F27%2F2024.07.27.24310936.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000321747800025&link_type=ISI) 

67. 67.Heard, E., and Martienssen, R.A. (2014). Transgenerational epigenetic inheritance: myths and mechanisms. Cell 157, 95–109.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.cell.2014.02.045&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24679529&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F27%2F2024.07.27.24310936.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000333968900010&link_type=ISI) 

68. 68.King, S.E., and Skinner, M.K. (2020). Epigenetic Transgenerational Inheritance of Obesity Susceptibility. Trends Endocrinol. Metab. 31, 478–494.
    
    

69. 69.Morisako, A.K., Tauali‘i, M., Ambrose, A.J.H., and Withy, K. (2017). Beyond the Ability to Pay: The Health Status of Native Hawaiians and Other Pacific Islanders in Relationship to Health Insurance. Hawaii J. Med. Public Health 76, 36.
    
    

70. 70.Heslehurst, N., Vieira, R., Akhter, Z., Bailey, H., Slack, E., Ngongalah, L., Pemu, A., and Rankin, J. (2019). The association between maternal body mass index and child obesity: A systematic review and meta-analysis. PLoS Med. 16, e1002817.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pmed.1002817&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=31185012&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F27%2F2024.07.27.24310936.atom) 

71. 71.Liu, W., Yang, X., Mao, Z., Du, Y., Lassiter, C., AlAkwaa, F.M., Benny, P.A., and Garmire, L.X. Severe preeclampsia is not associated with significant DNA methylation changes but cell proportion changes in the cord blood - caution on the importance of confounding adjustment. medRxiv. doi:10.1101/2023.08.31.23294898.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NzoibWVkcnhpdiI7czo1OiJyZXNpZCI7czoyMToiMjAyMy4wOC4zMS4yMzI5NDg5OHYxIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjQvMDcvMjcvMjAyNC4wNy4yNy4yNDMxMDkzNi5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

72. 72.Abu-Farha, M., Tiss, A., Abubaker, J., Khadir, A., Al-Ghimlas, F., Al-Khairi, I., Baturcam, E., Cherian, P., Elkum, N., Hammad, M., et al. (2013). Proteomics Analysis of Human Obesity Reveals the Epigenetic Factor HDAC4 as a Potential Target for Obesity. PLoS One 8. doi:10.1371/journal.pone.0075342.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pone.0075342&link_type=DOI) 

73. 73.Rönn, T., Volkov, P., Gillberg, L., Kokosar, M., Perfilyev, A., Jacobsen, A.L., Jørgensen, S.W., Brøns, C., Jansson, P.A., Eriksson, K.F., et al. (2015). Impact of age, BMI and HbA1c levels on the genome-wide DNA methylation and mRNA expression patterns in human adipose tissue and identification of epigenetic biomarkers in blood. Hum. Mol. Genet. 24. doi:10.1093/hmg/ddv124.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/hmg/ddv124&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25861810&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F27%2F2024.07.27.24310936.atom) 

74. 74.Lim, W.-J., Kim, K.H., Kim, J.-Y., Jeong, S., and Kim, N. (2019). Identification of DNA-Methylated CpG Islands Associated With Gene Silencing in the Adult Body Tissues of the Ogye Chicken Using RNA-Seq and Reduced Representation Bisulfite Sequencing. Front. Genet. 10, 346.
    
    

75. 75.Ching, T., Song, M.-A., Tiirikainen, M., Molnar, J., Berry, M., Towner, D., and Garmire, L.X. (2014). Genome-wide hypermethylation coupled with promoter hypomethylation in the chorioamniotic membranes of early onset pre-eclampsia. Mol. Hum. Reprod. 20, 885–904.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/molehr/gau046&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24944161&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F27%2F2024.07.27.24310936.atom) 

76. 76.Meng, D., Frank, A.R., and Jewell, J.L. (2018). mTOR signaling in stem and progenitor cells. Development 145. doi:10.1242/dev.152595.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NzoiZGV2ZWxvcCI7czo1OiJyZXNpZCI7czoxNToiMTQ1LzEvZGV2MTUyNTk1IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjQvMDcvMjcvMjAyNC4wNy4yNy4yNDMxMDkzNi5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

77. 77.Guo, M., Bao, E.L., Wagner, M., Whitsett, J.A., and Xu, Y. (2017). SLICE: determining cell differentiation and lineage based on single cell entropy. Nucleic Acids Res. 45, e54.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/nar/gkw1278&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=27998929&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F27%2F2024.07.27.24310936.atom) 

78. 78.Dai, D., Li, L., Huebner, A., Zeng, H., Guevara, E., Claypool, D.J., Liu, A., and Chen, J. (2012). Planar cell polarity effector gene Intu regulates cell fate-specific differentiation of keratinocytes through the primary cilia. Cell Death Differ. 20, 130–138.
    
    

79. 79.Liu, J., and Fuchs, S.Y. (2006). Cross-talk between APC/C and CBP/p300. Cancer Biol. Ther. 5. doi:10.4161/cbt.5.7.3103.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.4161/cbt.5.7.3103&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=16861917&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F27%2F2024.07.27.24310936.atom) 

80. 80.Corbi, S.C.T., Bastos, A.S., Nepomuceno, R., Cirelli, T., Dos Santos, R.A., Takahashi, C.S., Rocha, C.S., Orrico, S.R.P., Maurer-Morelli, C.V., and Scarel-Caminaga, R.M. (2017). Expression Profile of Genes Potentially Associated with Adequate Glycemic Control in Patients with Type 2 Diabetes Mellitus. Journal of diabetes research 2017. doi:10.1155/2017/2180819.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1155/2017/2180819&link_type=DOI) 

81. 81.Behl, T., Kaur, I., Sehgal, A., Singh, S., Zengin, G., Negrut, N., Nistor-Cseppento, D.C., Pavel, F.M., Aron, R.A.C., and Bungau, S. (2021). Exploring the Genetic Conception of Obesity via the Dual Role of FoxO. Int. J. Mol. Sci. 22. doi:10.3390/ijms22063179.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.3390/ijms22063179&link_type=DOI) 

82. 82.Kadakia, R., Zheng, Y., Zhang, Z., Zhang, W., Hou, L., and Josefson, J.L. (2017). Maternal pre-pregnancy BMI downregulates neonatal cord blood LEP methylation. Pediatr. Obes. 12 *Suppl 1*, 57–64.
    
    

83. 83.Yu, Z., Han, S., Zhu, J., Sun, X., Ji, C., and Guo, X. (2013). Pre-pregnancy body mass index in relation to infant birth weight and offspring overweight/obesity: a systematic review and meta-analysis. PLoS One 8, e61627.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pone.0061627&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23613888&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F27%2F2024.07.27.24310936.atom) 

84. 84.Sureshchandra, S., Marshall, N.E., and Messaoudi, I. (2019). Impact of pregravid obesity on maternal and fetal immunity: Fertile grounds for reprogramming. J. Leukoc. Biol. 106, 1035–1050.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/jlb.3ri0619-181r&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=31483523&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F27%2F2024.07.27.24310936.atom) 

85. 85.Levy, O. (2005). Innate immunity of the human newborn: distinct cytokine responses to LPS and other Toll-like receptor agonists. J. Endotoxin Res. 11, 113–116.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1177/09680519050110020701&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=15949138&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F27%2F2024.07.27.24310936.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000229439100007&link_type=ISI) 

86. 86.Ong, P.S., Wang, L.Z., Dai, X., Tseng, S.H., Loo, S.J., and Sethi, G. (2016). Judicious Toggling of mTOR Activity to Combat Insulin Resistance and Cancer: Current Evidence and Perspectives. Front. Pharmacol. 7. doi:10.3389/fphar.2016.00395.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.3389/fphar.2016.00395&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=27826244&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F27%2F2024.07.27.24310936.atom) 

87. 87.Mack, R., Zhang, L., Breslin, P., and Zhang, J. (2021). The fetal-to-adult hematopoietic stem cell transition and its role in childhood hematopoietic malignancies. Stem cell reviews and reports 17, 2059.
    
    

88. 88.Kwok, A.C., and Wong, J.T. (2005). Lipid biosynthesis and its coordination with cell cycle progression. Plant Cell Physiol. 46. doi:10.1093/pcp/pci213.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/pcp/pci213&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=16239308&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F27%2F2024.07.27.24310936.atom) 

89. 89.Denizli, M., Capitano, M.L., and Kua, K.L. (2022). Maternal obesity and the impact of associated early-life inflammation on long-term health of offspring. Front. Cell. Infect. Microbiol. 12, 940937.
    
    

90. 90.Meir, A.Y., Huang, W., Cao, T., Hong, X., Wang, G., Pearson, C., Adams, W.G., Wang, X., and Liang, L. (2023). Umbilical cord DNA methylation is associated with body mass index trajectories from birth to adolescence. EBioMedicine 91, 104550.
    
    

91. 91.Oesterreich, S., Allredl, D.C., Mohsin, S.K., Zhang, Q., Wong, H., Lee, A.V., Osborne, C.K., and O’Connell, P. (2001). High rates of loss of heterozygosity on chromosome 19p13 in human breast cancer. Br. J. Cancer 84, 493–498.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1054/bjoc.2000.1606&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=11207044&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F27%2F2024.07.27.24310936.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000167294600010&link_type=ISI) 

92. 92.Hong, E.A., Gautrey, H.L., Elliott, D.J., and Tyson-Capper, A.J. (2012). SAFB1- and SAFB2-mediated transcriptional repression: relevance to cancer. Biochem. Soc. Trans. 40, 826–830.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NzoicHBiaW9zdCI7czo1OiJyZXNpZCI7czo4OiI0MC80LzgyNiI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDI0LzA3LzI3LzIwMjQuMDcuMjcuMjQzMTA5MzYuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 

93. 93.Hammerich-Hille, S., Kaipparettu, B.A., Tsimelzon, A., Creighton, C.J., Jiang, S., Polo, J.M., Melnick, A., Meyer, R., and Oesterreich, S. (2010). SAFB1 mediates repression of immune regulators and apoptotic genes in breast cancer cells. J. Biol. Chem. 285, 3608–3616.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MzoiamJjIjtzOjU6InJlc2lkIjtzOjEwOiIyODUvNi8zNjA4IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjQvMDcvMjcvMjAyNC4wNy4yNy4yNDMxMDkzNi5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=)

 [1]: /embed/graphic-1.gif