Prenatal black carbon exposure and DNA methylation in umbilical cord blood.

BACKGROUND/OBJECTIVES: Prenatal exposure to ambient air pollution is associated with adverse cardiometabolic outcomes in childhood. We previously observed that prenatal black carbon (BC) was inversely associated with adiponectin, a hormone secreted by adipocytes, in early childhood. Changes to DNA methylation have been proposed as a potential mediator linking in utero exposures to lasting health impacts. METHODS: Among 532 mother-child pairs enrolled in the Colorado-based Healthy Start study, we performed an epigenome-wide association study of the relationship between prenatal exposure to a component of air pollution, BC, and DNA methylation in cord blood. Average pregnancy ambient BC was estimated at the mother's residence using a spatiotemporal prediction model. DNA methylation was measured using the Illumina 450K array. We used multiple linear regression to estimate associations between prenatal ambient BC and 429,246 cysteine-phosphate-guanine sites (CpGs), adjusting for potential confounders. We identified differentially methylated regions (DMRs) using DMRff and ENmix-combp. In a subset of participants (n = 243), we investigated DNA methylation as a potential mediator of the association between prenatal ambient BC and lower adiponectin in childhood. RESULTS: We identified 44 CpGs associated with average prenatal ambient BC after correcting for multiple testing. Several genes annotated to the top CpGs had reported functions in the immune system. There were 24 DMRs identified by both DMRff and ENmix-combp. One CpG (cg01123250), located on chromosome 2 and annotated to the UNC80 gene, was found to mediate approximately 20% of the effect of prenatal BC on childhood adiponectin, though the confidence interval was wide (95% CI: 3, 84). CONCLUSIONS: Prenatal BC was associated with DNA methylation in cord blood at several sites and regions in the genome. DNA methylation may partially mediate associations between prenatal BC and childhood cardiometabolic outcomes.

Proteomic networks and related genetic variants associated with smoking and chronic obstructive pulmonary disease.

BACKGROUND: Studies have identified individual blood biomarkers associated with chronic obstructive pulmonary disease (COPD) and related phenotypes. However, complex diseases such as COPD typically involve changes in multiple molecules with interconnections that may not be captured when considering single molecular features. METHODS: Leveraging proteomic data from 3,173 COPDGene Non-Hispanic White (NHW) and African American (AA) participants, we applied sparse multiple canonical correlation network analysis (SmCCNet) to 4,776 proteins assayed on the SomaScan v4.0 platform to derive sparse networks of proteins associated with current vs. former smoking status, airflow obstruction, and emphysema quantitated from high-resolution computed tomography scans. We then used NetSHy, a dimension reduction technique leveraging network topology, to produce summary scores of each proteomic network, referred to as NetSHy scores. We next performed a genome-wide association study (GWAS) to identify variants associated with the NetSHy scores, or network quantitative trait loci (nQTLs). Finally, we evaluated the replicability of the networks in an independent cohort, SPIROMICS. RESULTS: We identified networks of 13 to 104 proteins for each phenotype and exposure in NHW and AA, and the derived NetSHy scores significantly associated with the variable of interests. Networks included known (sRAGE, ALPP, MIP1) and novel molecules (CA10, CPB1, HIS3, PXDN) and interactions involved in COPD pathogenesis. We observed 7 nQTL loci associated with NetSHy scores, 4 of which remained after conditional analysis. Networks for smoking status and emphysema, but not airflow obstruction, demonstrated a high degree of replicability across race groups and cohorts. CONCLUSIONS: In this work, we apply state-of-the-art molecular network generation and summarization approaches to proteomic data from COPDGene participants to uncover protein networks associated with COPD phenotypes. We further identify genetic associations with networks. This work discovers protein networks containing known and novel proteins and protein interactions associated with clinically relevant COPD phenotypes across race groups and cohorts.

Proteomic Networks and Related Genetic Variants Associated with Smoking and Chronic Obstructive Pulmonary Disease.

Background: Studies have identified individual blood biomarkers associated with chronic obstructive pulmonary disease (COPD) and related phenotypes. However, complex diseases such as COPD typically involve changes in multiple molecules with interconnections that may not be captured when considering single molecular features. Methods: Leveraging proteomic data from 3,173 COPDGene Non-Hispanic White (NHW) and African American (AA) participants, we applied sparse multiple canonical correlation network analysis (SmCCNet) to 4,776 proteins assayed on the SomaScan v4.0 platform to derive sparse networks of proteins associated with current vs. former smoking status, airflow obstruction, and emphysema quantitated from high-resolution computed tomography scans. We then used NetSHy, a dimension reduction technique leveraging network topology, to produce summary scores of each proteomic network, referred to as NetSHy scores. We next performed genome-wide association study (GWAS) to identify variants associated with the NetSHy scores, or network quantitative trait loci (nQTLs). Finally, we evaluated the replicability of the networks in an independent cohort, SPIROMICS. Results: We identified networks of 13 to 104 proteins for each phenotype and exposure in NHW and AA, and the derived NetSHy scores significantly associated with the variable of interests. Networks included known (sRAGE, ALPP, MIP1) and novel molecules (CA10, CPB1, HIS3, PXDN) and interactions involved in COPD pathogenesis. We observed 7 nQTL loci associated with NetSHy scores, 4 of which remained after conditional analysis. Networks for smoking status and emphysema, but not airflow obstruction, demonstrated a high degree of replicability across race groups and cohorts. Conclusions: In this work, we apply state-of-the-art molecular network generation and summarization approaches to proteomic data from COPDGene participants to uncover protein networks associated with COPD phenotypes. We further identify genetic associations with networks. This work discovers protein networks containing known and novel proteins and protein interactions associated with clinically relevant COPD phenotypes across race groups and cohorts.

TreeKernel: interpretable kernel machine tests for interactions between -omics and clinical predictors with applications to metabolomics and COPD phenotypes.

BACKGROUND: In this paper, we are interested in interactions between a high-dimensional -omics dataset and clinical covariates. The goal is to evaluate the relationship between a phenotype of interest and a high-dimensional omics pathway, where the effect of the omics data depends on subjects' clinical covariates (age, sex, smoking status, etc.). For instance, metabolic pathways can vary greatly between sexes which may also change the relationship between certain metabolic pathways and a clinical phenotype of interest. We propose partitioning the clinical covariate space and performing a kernel association test within those partitions. To illustrate this idea, we focus on hierarchical partitions of the clinical covariate space and kernel tests on metabolic pathways. RESULTS: We see that our proposed method outperforms competing methods in most simulation scenarios. It can identify different relationships among clinical groups with higher power in most scenarios while maintaining a proper Type I error rate. The simulation studies also show a robustness to the grouping structure within the clinical space. We also apply the method to the COPDGene study and find several clinically meaningful interactions between metabolic pathways, the clinical space, and lung function. CONCLUSION: TreeKernel provides a simple and interpretable process for testing for relationships between high-dimensional omics data and clinical outcomes in the presence of interactions within clinical cohorts. The method is broadly applicable to many studies.

Adipocyte hypertrophy in mesenchymal stem cells from infants of mothers with obesity.

OBJECTIVE: Fat content of adipocytes derived from infant umbilical cord mesenchymal stem cells (MSCs) predicts adiposity in children through 4 to 6 years of age. This study tested the hypothesis that MSCs from infants born to mothers with obesity (Ob-MSCs) exhibit adipocyte hypertrophy and perturbations in genes regulating adipogenesis compared with MSCs from infants of mothers with normal weight (NW-MSCs). METHODS: Adipogenesis was induced in MSCs embedded in three-dimensional hydrogel structures, and cell size and number were measured by three-dimensional imaging. Proliferation and protein markers of proliferation and adipogenesis in undifferentiated and adipocyte differentiating cells were measured. RNA sequencing was performed to determine pathways linked to adipogenesis phenotype. RESULTS: In undifferentiated MSCs, greater zinc finger protein (Zfp)423 protein content was observed in Ob- versus NW-MSCs. Adipocytes from Ob-MSCs were larger but fewer than adipocytes from NW-MSCs. RNA sequencing analysis showed that Zfp423 protein correlated with mRNA expression of genes enriched for cell cycle, MSC lineage specification, inflammation, and metabolism pathways. MSC proliferation was not different before differentiation but declined faster in Ob-MSCs upon adipogenic induction. CONCLUSIONS: Ob-MSCs have an intrinsic propensity for adipocyte hypertrophy and reduced hyperplasia during adipogenesis in vitro, perhaps linked to greater Zfp423 content and changes in cell cycle pathway gene expression.

Novel Metabolic Subtypes in Pregnant Women and Risk of Early Childhood Obesity in Offspring.

Importance: The in utero metabolic milieu is associated with offspring adiposity. Standard definitions of maternal obesity (according to prepregnancy body mass index [BMI]) and gestational diabetes (GDM) may not be adequate to capture subtle yet important differences in the intrauterine environment that could be involved in programming. Objectives: To identify maternal metabolic subgroups during pregnancy and to examine associations of subgroup classification with adiposity traits in their children. Design, Setting, and Participants: This cohort study included mother-offspring pairs in the Healthy Start prebirth cohort (enrollment: 2010-2014) recruited from University of Colorado Hospital obstetrics clinics in Aurora, Colorado. Follow-up of women and children is ongoing. Data were analyzed from March to December 2022. Exposures: Metabolic subtypes of pregnant women ascertained by applying k-means clustering on 7 biomarkers and 2 biomarker indices measured at approximately 17 gestational weeks: glucose, insulin, Homeostatic Model Assessment for Insulin Resistance, total cholesterol, high-density lipoprotein cholesterol (HDL-C), triglycerides, free fatty acids (FFA), HDL-C:triglycerides ratio, and tumor necrosis factor α. Main Outcomes and Measures: Offspring birthweight z score and neonatal fat mass percentage (FM%). In childhood at approximately 5 years of age, offspring BMI percentile, FM%, BMI in the 95th percentile or higher, and FM% in the 95th percentile or higher. Results: A total of 1325 pregnant women (mean [SD] age, 27.8 [6.2 years]; 322 [24.3%] Hispanic, 207 non-Hispanic Black [15.6%], and 713 [53.8%] non-Hispanic White), and 727 offspring with anthropometric data measured in childhood (mean [SD] age 4.81 [0.72] years, 48% female) were included. We identified the following 5 maternal metabolic subgroups: reference (438 participants), high HDL-C (355 participants), dyslipidemic-high triglycerides (182 participants), dyslipidemic-high FFA (234 participants), and insulin resistant (IR)-hyperglycemic (116 participants). Compared with the reference subgroup, women in the IR-hyperglycemic and dyslipidemic-high FFA subgroups had offspring with 4.27% (95% CI, 1.94-6.59) and 1.96% (95% CI, 0.45-3.47) greater FM% during childhood, respectively. There was a higher risk of high FM% among offspring of the IR-hyperglycemic (relative risk, 8.7; 95% CI, 2.7-27.8) and dyslipidemic-high FFA (relative risk, 3.4; 95% CI, 1.0-11.3) subgroups; this risk was of greater magnitude compared with prepregnancy obesity alone, GDM alone, or both conditions. Conclusions and Relevance: In this cohort study, an unsupervised clustering approach revealed distinct metabolic subgroups of pregnant women. These subgroups exhibited differences in risk of offspring adiposity in early childhood. Such approaches have the potential to refine understanding of the in utero metabolic milieu, with utility for capturing variation in sociocultural, anthropometric, and biochemical risk factors for offspring adiposity.

Fat content in infant mesenchymal stem cells prospectively associates with childhood adiposity and fasting glucose.

OBJECTIVE: In human studies, new model systems are needed for improved mechanistic investigation of developmental predisposition for metabolic disease but also to serve as benchmarks in early life prevention or intervention efforts. In this regard, human infant umbilical cord-derived mesenchymal stem cells (MSCs) are an emerging tool. However, long-term clinical relevance to in vivo markers of metabolic disease is unknown. METHODS: In a cohort of 124 mother/child dyads, this study tested the hypothesis that triglyceride content (TG) of infant MSCs undergoing adipogenesis in vitro (MSC-TG) is associated with in vivo adiposity (percent fat mass) from birth to early childhood and with fasting glucose and insulin in early childhood. RESULTS: MSC-TG was positively associated with in vivo child adiposity at birth, age 4 to 6 months, and age 4 to 6 years. MSC-TG was associated with fasting glucose, but not insulin, at 4 to 6 years. Importantly, MSC-TG explained an additional 13% variance in child adiposity at 4 to 6 years, after accounting for other established birth predictors (weight and percent fat mass at birth) and other established covariates related to child adiposity (e.g., breastfeeding exposure, physical activity). CONCLUSIONS: This work demonstrates the strength of the MSC model for predicting offspring metabolic phenotype into childhood, even when considering the important contribution of other early life risk factors.

Ambient air pollution during pregnancy and DNA methylation in umbilical cord blood, with potential mediation of associations with infant adiposity: The Healthy Start study.

BACKGROUND: Prenatal exposure to ambient air pollution has been associated with adverse offspring health outcomes. Childhood health effects of prenatal exposures may be mediated through changes to DNA methylation detectable at birth. METHODS: Among 429 non-smoking women in a cohort study of mother-infant pairs in Colorado, USA, we estimated associations between prenatal exposure to ambient fine particulate matter (PM2.5) and ozone (O3), and epigenome-wide DNA methylation of umbilical cord blood cells at delivery (2010-2014). We calculated average PM2.5 and O3 in each trimester of pregnancy and the full pregnancy using inverse-distance-weighted interpolation. We fit linear regression models adjusted for potential confounders and cell proportions to estimate associations between air pollutants and methylation at each of 432,943 CpGs. Differentially methylated regions (DMRs) were identified using comb-p. Previously in this cohort, we reported positive associations between 3rd trimester O3 exposure and infant adiposity at 5 months of age. Here, we quantified the potential for mediation of that association by changes in DNA methylation in cord blood. RESULTS: We identified several DMRs for each pollutant and period of pregnancy. The greatest number of significant DMRs were associated with third trimester PM2.5 (21 DMRs). No single CpGs were associated with air pollutants at a false discovery rate <0.05. We found that up to 8% of the effect of 3rd trimester O3 on 5-month adiposity may be mediated by locus-specific methylation changes, but mediation estimates were not statistically significant. CONCLUSIONS: Differentially methylated regions in cord blood were identified in association with maternal exposure to PM2.5 and O3. Genes annotated to the significant sites played roles in cardiometabolic disease, immune function and inflammation, and neurologic disorders. We found limited evidence of mediation by DNA methylation of associations between third trimester O3 exposure and 5-month infant adiposity.

Metabolomic Profiles in Childhood and Adolescence Are Associated with Fetal Overnutrition.

Fetal overnutrition predisposes offspring to increased metabolic risk. The current study used metabolomics to assess sustained differences in serum metabolites across childhood and adolescence among youth exposed to three typologies of fetal overnutrition: maternal obesity only, gestational diabetes mellitus (GDM) only, and obesity + GDM. We included youth exposed in utero to obesity only (BMI ≥ 30; n = 66), GDM only (n = 56), obesity + GDM (n = 25), or unexposed (n = 297), with untargeted metabolomics measured at ages 10 and 16 years. We used linear mixed models to identify metabolites across both time-points associated with exposure to any overnutrition, using a false-discovery-rate correction (FDR) <0.20. These metabolites were included in a principal component analysis (PCA) to generate profiles and assess metabolite profile differences with respect to overnutrition typology (adjusted for prenatal smoking, offspring age, sex, and race/ethnicity). Fetal overnutrition was associated with 52 metabolites. PCA yielded four factors accounting for 17−27% of the variance, depending on age of measurement. We observed differences in three factor patterns with respect to overnutrition typology: sphingomyelin-mannose (8−13% variance), skeletal muscle metabolism (6−10% variance), and 3-carboxy-4-methyl-5-propyl-2-furanpropanoic acid (CMPF; 3−4% variance). The sphingomyelin-mannose factor score was higher among offspring exposed to obesity vs. GDM. Exposure to obesity + GDM (vs. GDM or obesity only) was associated with higher skeletal muscle metabolism and CMPF scores. Fetal overnutrition is associated with metabolic changes in the offspring, but differences between typologies of overnutrition account for a small amount of variation in the metabolome, suggesting there is likely greater pathophysiological overlap than difference.

Longitudinal analysis of SARS-CoV-2 spike and RNA-dependent RNA polymerase protein sequences reveals the emergence and geographic distribution of diverse mutations.

Amid the ongoing COVID-19 pandemic, it has become increasingly important to monitor the mutations that arise in the SARS-CoV-2 virus, to prepare public health strategies and guide the further development of vaccines and therapeutics. The spike (S) protein and the proteins comprising the RNA-Dependent RNA Polymerase (RdRP) are key vaccine and drug targets, respectively, making mutation surveillance of these proteins of great importance. Full protein sequences were downloaded from the GISAID database, aligned, and the variants identified. 437,006 unique viral genomes were analyzed. Polymorphisms in the protein sequence were investigated and examined longitudinally to identify sequence and strain variants appearing between January 5th, 2020 and January 16th, 2021. A structural analysis was also performed to investigate mutations in the receptor binding domain and the N-terminal domain of the spike protein. Within the spike protein, there were 766 unique mutations observed in the N-terminal domain and 360 in the receptor binding domain. Four residues that directly contact ACE2 were mutated in more than 100 sequences, including positions K417, Y453, S494, and N501. Within the furin cleavage site of the spike protein, a high degree of conservation was observed, but the P681H mutation was observed in 10.47% of sequences analyzed. Within the RNA dependent RNA polymerase complex proteins, 327 unique mutations were observed in Nsp8, 166 unique mutations were observed in Nsp7, and 1157 unique mutations were observed in Nsp12. Only 4 sequences analyzed contained mutations in the 9 residues that directly interact with the therapeutic Remdesivir, suggesting limited mutations in drug interacting residues. The identification of new variants emphasizes the need for further study on the effects of the mutations and the implications of increased prevalence, particularly for vaccine or therapeutic efficacy.

Identifying miRNA-mRNA Networks Associated With COPD Phenotypes.

Chronic obstructive pulmonary disease (COPD) is characterized by expiratory airflow limitation and symptoms such as shortness of breath. Although many studies have demonstrated dysregulated microRNA (miRNA) and gene (mRNA) expression in the pathogenesis of COPD, how miRNAs and mRNAs systematically interact and contribute to COPD development is still not clear. To gain a deeper understanding of the gene regulatory network underlying COPD pathogenesis, we used Sparse Multiple Canonical Correlation Network (SmCCNet) to integrate whole blood miRNA and RNA-sequencing data from 404 participants in the COPDGene study to identify novel miRNA-mRNA networks associated with COPD-related phenotypes including lung function and emphysema. We hypothesized that phenotype-directed interpretable miRNA-mRNA networks from SmCCNet would assist in the discovery of novel biomarkers that traditional single biomarker discovery methods (such as differential expression) might fail to discover. Additionally, we investigated whether adjusting -omics and clinical phenotypes data for covariates prior to integration would increase the statistical power for network identification. Our study demonstrated that partial covariate adjustment for age, sex, race, and CT scanner model (in the quantitative emphysema networks) improved network identification when compared with no covariate adjustment. However, further adjustment for current smoking status and relative white blood cell (WBC) proportions sometimes weakened the power for identifying lung function and emphysema networks, a phenomenon which may be due to the correlation of smoking status and WBC counts with the COPD-related phenotypes. With partial covariate adjustment, we found six miRNA-mRNA networks associated with COPD-related phenotypes. One network consists of 2 miRNAs and 28 mRNAs which had a 0.33 correlation (p = 5.40E-12) to forced expiratory volume in 1 s (FEV1) percent predicted. We also found a network of 5 miRNAs and 81 mRNAs that had a 0.45 correlation (p = 8.80E-22) to percent emphysema. The miRNA-mRNA networks associated with COPD traits provide a systems view of COPD pathogenesis and complements biomarker identification with individual miRNA or mRNA expression data.

Soluble receptor for advanced glycation end products (sRAGE) as a biomarker of COPD.

BACKGROUND: Soluble receptor for advanced glycation end products (sRAGE) is a proposed emphysema and airflow obstruction biomarker; however, previous publications have shown inconsistent associations and only one study has investigate the association between sRAGE and emphysema. No cohorts have examined the association between sRAGE and progressive decline of lung function. There have also been no evaluation of assay compatibility, receiver operating characteristics, and little examination of the effect of genetic variability in non-white population. This manuscript addresses these deficiencies and introduces novel data from Pittsburgh COPD SCCOR and as well as novel work on airflow obstruction. A meta-analysis is used to quantify sRAGE associations with clinical phenotypes. METHODS: sRAGE was measured in four independent longitudinal cohorts on different analytic assays: COPDGene (n = 1443); SPIROMICS (n = 1623); ECLIPSE (n = 2349); Pittsburgh COPD SCCOR (n = 399). We constructed adjusted linear mixed models to determine associations of sRAGE with baseline and follow up forced expiratory volume at one second (FEV1) and emphysema by quantitative high-resolution CT lung density at the 15th percentile (adjusted for total lung capacity). RESULTS: Lower plasma or serum sRAGE values were associated with a COPD diagnosis (P < 0.001), reduced FEV1 (P < 0.001), and emphysema severity (P < 0.001). In an inverse-variance weighted meta-analysis, one SD lower log10-transformed sRAGE was associated with 105 ± 22 mL lower FEV1 and 4.14 ± 0.55 g/L lower adjusted lung density. After adjusting for covariates, lower sRAGE at baseline was associated with greater FEV1 decline and emphysema progression only in the ECLIPSE cohort. Non-Hispanic white subjects carrying the rs2070600 minor allele (A) and non-Hispanic African Americans carrying the rs2071288 minor allele (A) had lower sRAGE measurements compare to those with the major allele, but their emphysema-sRAGE regression slopes were similar. CONCLUSIONS: Lower blood sRAGE is associated with more severe airflow obstruction and emphysema, but associations with progression are inconsistent in the cohorts analyzed. In these cohorts, genotype influenced sRAGE measurements and strengthened variance modelling. Thus, genotype should be included in sRAGE evaluations.

Alcohol Use Disorders Are Associated With a Unique Impact on Airway Epithelial Cell Gene Expression.

BACKGROUND: Alcohol use disorders (AUDs) and cigarette smoking both increase risk for the development of community-acquired pneumonia (CAP), likely through adverse effects on proximal airway mucociliary clearance and pathogen recognition. Smoking-related alterations on airway gene expression are well described, but little is known about the impact of AUDs. We measured gene expression in human airway epithelial cells (AECs), hypothesizing that AUDs would be associated with novel differences in gene expression that could alter risk for CAP. METHODS: Bronchoscopy with airway brushings was performed in participants with AUDs and controls to obtain AECs. An AUD Identification Test was used to define AUD. RNA was extracted from AECs, and mRNA expression data were collected on an Agilent micro-array. Differential expression analyses were performed on the filtered and normalized data with correction for multiple testing. Enrichment analyses were performed using clusterProfiler. RESULTS: Expression data from 19 control and 18 AUD participants were evaluated. After adjustment for smoking, AUDs were associated with significant differential expression of 520 AEC genes, including genes for ribosomal proteins and genes involved in protein folding. Enrichment analyses indicated significant differential expression of 24 pathways in AUDs, including those implicated in protein targeting to membrane and viral gene expression. Smoking-associated AEC gene expression differences mirrored previous reports, but differed from those associated with AUDs. CONCLUSIONS: AUDs have a distinct impact on AEC gene expression that may influence proximal airway function independent of smoking. Alcohol-associated alterations may influence risk for CAP through modifying key mechanisms important in protecting proximal airway integrity.

DNA methylation near the INS gene is associated with INS genetic variation (rs689) and type 1 diabetes in the Diabetes Autoimmunity Study in the Young.

OBJECTIVE: Mechanisms underlying the role of non-human leukocyte antigen (HLA) genetic risk variants in type 1 diabetes (T1D) are poorly understood. We aimed to test the association between methylation and non-HLA genetic risk. METHODS: We conducted a methylation quantitative trait loci (mQTL) analysis in a nested case-control study from the Dietary Autoimmunity Study in the Young. Controls (n = 83) were frequency-matched to T1D cases (n = 83) based on age, race/ethnicity, and sample availability. We evaluated 13 non-HLA genetic markers known be associated with T1D. Genome-wide methylation profiling was performed on peripheral blood samples collected prior to T1D using the Illumina 450 K (discovery set) and infinium methylation EPIC beadchip (EPIC validation) platforms. Linear regression models, adjusting for age and sex, were used to test to each single nucleotide polymorphism (SNP) -probe combination. Logistic regression models were used to test the association between T1D and methylation levels among probes with a significant mQTL. A meta-analysis was used to combine odds ratios from the two platforms. RESULTS: We identified 10 SNP-methylation probe pairs (false discovery rate (FDR) adjusted P < .05 and validation P < .05). Probes were associated with the GSDMB, C1QTNF6, IL27, and INS genes. The cg03366382 (OR: 1.9, meta-P = .0495), cg21574853 (OR: 2.5, meta-P = .0232), and cg25336198 (odds ratio: 6.6, meta-P = .0081) probes were significantly associated with T1D. The three probes were located upstream from the INS transcription start site. CONCLUSIONS: We confirmed an association between DNA methylation and rs689 that has been identified in related studies. Measurements in our study preceded the onset of T1D suggesting methylation may have a role in the relationship between INS variation and T1D development.

A prospective study of associations between in utero exposure to gestational diabetes mellitus and metabolomic profiles during late childhood and adolescence.

AIMS/HYPOTHESIS: This study aimed to: (1) identify metabolite patterns during late childhood that differ with respect to exposure to maternal gestational diabetes mellitus (GDM); (2) examine the persistence of GDM/metabolite associations 5 years later, during adolescence; and (3) investigate the associations of metabolite patterns with adiposity and metabolic biomarkers from childhood through adolescence. METHODS: This study included 592 mother-child pairs with information on GDM exposure (n = 92 exposed), untargeted metabolomics data at age 6-14 years (T1) and at 12-19 years (T2), and information on adiposity and metabolic risk biomarkers at T1 and T2. We first consolidated 767 metabolites at T1 into factors (metabolite patterns) via principal component analysis (PCA) and used multivariable regression to identify factors that differed by GDM exposure, at α = 0.05. We then examined associations of GDM with individual metabolites within factors of interest at T1 and T2, and investigated associations of GDM-related factors at T1 with adiposity and metabolic risk throughout T1 and T2 using mixed-effects linear regression models. RESULTS: Of the six factors retained from PCA, GDM exposure was associated with greater odds of being in quartile (Q)4 (vs Q1-3) of 'Factor 4' at T1 after accounting for age, sex, race/ethnicity, maternal smoking habits during pregnancy, Tanner stage, physical activity and total energy intake, at α = 0.05 (OR 1.78 [95% CI 1.04, 3.04]; p = 0.04). This metabolite pattern comprised phosphatidylcholines, diacylglycerols and phosphatidylethanolamines. GDM was consistently associated with elevations in a subset of individual compounds within this pattern at T1 and T2. While this metabolite pattern was not related to the health outcomes in boys, it corresponded with greater adiposity and a worse metabolic profile among girls throughout the follow-up period. Each 1-unit increment in Factor 4 corresponded with 0.17 (0.08, 0.25) units higher BMI z score, 8.83 (5.07, 12.59) pmol/l higher fasting insulin, 0.28 (0.13, 0.43) units higher HOMA-IR, and 4.73 (2.15, 7.31) nmol/l higher leptin. CONCLUSIONS/INTERPRETATION: Exposure to maternal GDM was nominally associated with a metabolite pattern characterised by elevated serum phospholipids in late childhood and adolescence at α = 0.05. This metabolite pattern was associated with greater adiposity and metabolic risk among female offspring throughout the late childhood-to-adolescence transition. Future studies are warranted to confirm our findings.

Bronchoalveolar Lavage Fluid from COPD Patients Reveals More Compounds Associated with Disease than Matched Plasma.

Smoking causes chronic obstructive pulmonary disease (COPD). Though recent studies identified a COPD metabolomic signature in blood, no large studies examine the metabolome in bronchoalveolar lavage (BAL) fluid, a more direct representation of lung cell metabolism. We performed untargeted liquid chromatography-mass spectrometry (LC-MS) on BAL and matched plasma from 115 subjects from the SPIROMICS cohort. Regression was performed with COPD phenotypes as the outcome and metabolites as the predictor, adjusted for clinical covariates and false discovery rate. Weighted gene co-expression network analysis (WGCNA) grouped metabolites into modules which were then associated with phenotypes. K-means clustering grouped similar subjects. We detected 7939 and 10,561 compounds in BAL and paired plasma samples, respectively. FEV1/FVC (Forced Expiratory Volume in One Second/Forced Vital Capacity) ratio, emphysema, FEV1 % predicted, and COPD exacerbations associated with 1230, 792, eight, and one BAL compounds, respectively. Only two plasma compounds associated with a COPD phenotype (emphysema). Three BAL co-expression modules associated with FEV1/FVC and emphysema. K-means BAL metabolomic signature clustering identified two groups, one with more airway obstruction (34% of subjects, median FEV1/FVC 0.67), one with less (66% of subjects, median FEV1/FVC 0.77; p < 2 × 10-4). Associations between metabolites and COPD phenotypes are more robustly represented in BAL compared to plasma.

Unsupervised discovery of phenotype-specific multi-omics networks.

MOTIVATION: Complex diseases often involve a wide spectrum of phenotypic traits. Better understanding of the biological mechanisms relevant to each trait promotes understanding of the etiology of the disease and the potential for targeted and effective treatment plans. There have been many efforts towards omics data integration and network reconstruction, but limited work has examined the incorporation of relevant (quantitative) phenotypic traits. RESULTS: We propose a novel technique, sparse multiple canonical correlation network analysis (SmCCNet), for integrating multiple omics data types along with a quantitative phenotype of interest, and for constructing multi-omics networks that are specific to the phenotype. As a case study, we focus on miRNA-mRNA networks. Through simulations, we demonstrate that SmCCNet has better overall prediction performance compared to popular gene expression network construction and integration approaches under realistic settings. Applying SmCCNet to studies on chronic obstructive pulmonary disease (COPD) and breast cancer, we found enrichment of known relevant pathways (e.g. the Cadherin pathway for COPD and the interferon-gamma signaling pathway for breast cancer) as well as less known omics features that may be important to the diseases. Although those applications focus on miRNA-mRNA co-expression networks, SmCCNet is applicable to a variety of omics and other data types. It can also be easily generalized to incorporate multiple quantitative phenotype simultaneously. The versatility of SmCCNet suggests great potential of the approach in many areas. AVAILABILITY AND IMPLEMENTATION: The SmCCNet algorithm is written in R, and is freely available on the web at https://cran.r-project.org/web/packages/SmCCNet/index.html. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Metabolomics and transcriptomics pathway approach reveals outcome-specific perturbations in COPD.

Chronic obstructive pulmonary disease (COPD) comprises multiple phenotypes such as airflow obstruction, emphysema, and frequent episodes of acute worsening of respiratory symptoms, known as exacerbations. The goal of this pilot study was to test the usefulness of unbiased metabolomics and transcriptomics approaches to delineate biological pathways associated with COPD phenotypes and outcomes. Blood was collected from 149 current or former smokers with or without COPD and separated into peripheral blood mononuclear cells (PBMC) and plasma. PBMCs and plasma were analyzed using microarray and liquid chromatography mass spectrometry, respectively. Statistically significant transcripts and compounds were mapped to pathways using IMPaLA. Results showed that glycerophospholipid metabolism was associated with worse airflow obstruction and more COPD exacerbations. Sphingolipid metabolism was associated with worse lung function outcomes and exacerbation severity requiring hospitalizations. The strongest associations between a pathway and a certain COPD outcome were: fat digestion and absorption and T cell receptor signaling with lung function outcomes; antigen processing with exacerbation frequency; arginine and proline metabolism with exacerbation severity; and oxidative phosphorylation with emphysema. Overlaying transcriptomic and metabolomics datasets across pathways enabled outcome and phenotypic differences to be determined. Findings are relevant for identifying molecular targets for animal intervention studies and early intervention markers in human cohorts.

Condition-adaptive fused graphical lasso (CFGL): An adaptive procedure for inferring condition-specific gene co-expression network.

Co-expression network analysis provides useful information for studying gene regulation in biological processes. Examining condition-specific patterns of co-expression can provide insights into the underlying cellular processes activated in a particular condition. One challenge in this type of analysis is that the sample sizes in each condition are usually small, making the statistical inference of co-expression patterns highly underpowered. A joint network construction that borrows information from related structures across conditions has the potential to improve the power of the analysis. One possible approach to constructing the co-expression network is to use the Gaussian graphical model. Though several methods are available for joint estimation of multiple graphical models, they do not fully account for the heterogeneity between samples and between co-expression patterns introduced by condition specificity. Here we develop the condition-adaptive fused graphical lasso (CFGL), a data-driven approach to incorporate condition specificity in the estimation of co-expression networks. We show that this method improves the accuracy with which networks are learned. The application of this method on a rat multi-tissue dataset and The Cancer Genome Atlas (TCGA) breast cancer dataset provides interesting biological insights. In both analyses, we identify numerous modules enriched for Gene Ontology functions and observe that the modules that are upregulated in a particular condition are often involved in condition-specific activities. Interestingly, we observe that the genes strongly associated with survival time in the TCGA dataset are less likely to be network hubs, suggesting that genes associated with cancer progression are likely to govern specific functions or execute final biological functions in pathways, rather than regulating a large number of biological processes. Additionally, we observed that the tumor-specific hub genes tend to have few shared edges with normal tissue, revealing tumor-specific regulatory mechanism.

Maternal obesity alters fatty acid oxidation, AMPK activity, and associated DNA methylation in mesenchymal stem cells from human infants.

OBJECTIVE: Infants born to mothers with obesity have greater adiposity, ectopic fat storage, and are at increased risk for childhood obesity and metabolic disease compared with infants of normal weight mothers, though the cellular mechanisms mediating these effects are unclear. METHODS: We tested the hypothesis that human, umbilical cord-derived mesenchymal stem cells (MSCs) from infants born to obese (Ob-MSC) versus normal weight (NW-MSC) mothers demonstrate altered fatty acid metabolism consistent with adult obesity. In infant MSCs undergoing myogenesis in vitro, we measured cellular lipid metabolism and AMPK activity, AMPK activation in response to cellular nutrient stress, and MSC DNA methylation and mRNA content of genes related to oxidative metabolism. RESULTS: We found that Ob-MSCs exhibit greater lipid accumulation, lower fatty acid oxidation (FAO), and dysregulation of AMPK activity when undergoing myogenesis in vitro. Further experiments revealed a clear phenotype distinction within the Ob-MSC group where more severe MSC metabolic perturbation corresponded to greater neonatal adiposity and umbilical cord blood insulin levels. Targeted analysis of DNA methylation array revealed Ob-MSC hypermethylation in genes regulating FAO (PRKAG2, ACC2, CPT1A, SDHC) and corresponding lower mRNA content of these genes. Moreover, MSC methylation was positively correlated with infant adiposity. CONCLUSIONS: These data suggest that greater infant adiposity is associated with suppressed AMPK activity and reduced lipid oxidation in MSCs from infants born to mothers with obesity and may be an important, early marker of underlying obesity risk.