Search

Abstract 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.

Abstract Background Lung adenocarcinoma (LUAD) has been observed to have significant sex differences in incidence, prognosis, and response to therapy. However, the molecular mechanisms responsible for these disparities have not been investigated extensively. Methods Sample-specific gene regulatory network methods were used to analyze RNA sequencing data from non-cancerous human lung samples from The Genotype Tissue Expression Project (GTEx) and lung adenocarcinoma primary tumor samples from The Cancer Genome Atlas (TCGA); results were validated on independent data. Results We found that genes associated with key biological pathways including cell proliferation, immune response and drug metabolism are differentially regulated between males and females in both healthy lung tissue and tumor, and that these regulatory differences are further perturbed by tobacco smoking. We also discovered significant sex bias in transcription factor targeting patterns of clinically actionable oncogenes and tumor suppressor genes, including AKT2 and KRAS. Using differentially regulated genes between healthy and tumor samples in conjunction with a drug repurposing tool, we identified several small-molecule drugs that might have sex-biased efficacy as cancer therapeutics and further validated this observation using an independent cell line database. Conclusions These findings underscore the importance of including sex as a biological variable and considering gene regulatory processes in developing strategies for disease prevention and management. Graphical Abstract

ABSTRACTModelling the regulatory mechanisms that determine cell fate, response to external perturbation, and disease state depends on measuring many factors, a task made more difficult by the plasticity of the epigenome. Scanning the genome for the sequence patterns defined by Position Weight Matrices (PWM) can be used to estimate transcription factor (TF) binding locations. However, this approach does not incorporate information regarding the epigenetic context necessary for TF binding. CpG methylation is an epigenetic mark influenced by environmental factors that is commonly assayed in human cohort studies. We developed a framework to score inferred TF binding locations using methylation data. We intersected motif locations identified using PWMs with methylation information captured in both whole-genome bisulfite sequencing and Illumina EPIC array data for six cell lines, scored motif locations based on these data, and compared with experimental data characterizing TF binding (ChIP-seq). We found that for most TFs, binding prediction improves using methylation-based scoring compared to standard PWM-scores. We also illustrate that our approach can be generalized to infer TF binding when methylation information is only proximally available, i.e. measured for nearby CpGs that do not directly overlap with a motif location. Overall, our approach provides a framework for inferring context-specific TF binding using methylation data. Importantly, the availability of DNA methylation data in existing patient populations provides an opportunity to use our approach to understand the impact of methylation on gene regulatory processes in the context of human disease.

BackgroundBeyond exposure to cigarette smoking and aging, the factors that influence lung function decline to incident chronic obstructive pulmonary disease (COPD) remain unclear. Advancements have been made in categorizing COPD into emphysema and airway predominant disease subtypes; however, predicting which healthy individuals will progress to COPD is difficult because they can exhibit profoundly different disease trajectories despite similar initial risk factors. This study aimed to identify clinical, genetic, and radiological features that are directly linked-and subsequently predict-abnormal lung function.Methods and findingsWe employed graph modeling on 2,643 COPDGene participants (aged 45 to 80 years, 51.25% female, 35.1% African Americans; enrollment 11/2007-4/2011) with smoking history but normal spirometry at study enrollment to identify variables that are directly linked to future lung function abnormalities. We developed logistic regression and random forest predictive models for distinguishing individuals who maintain lung function from those who decline. Of the 131 variables analyzed, 6 were identified as informative to future lung function abnormalities, namely forced expiratory flow in the middle range (FEF25-75%), average lung wall thickness in a 10 mm radius (Pi10), severe emphysema, age, sex, and height. We investigated whether these features predict individuals leaving GOLD 0 status (normal spirometry according to Global Initiative for Obstructive Lung Disease (GOLD) criteria). Linear models, trained with these features, were quite predictive (area under receiver operator characteristic curve or AUROC = 0.75). Random forest predictors performed similarly to logistic regression (AUROC = 0.7), indicating that no significant nonlinear effects were present. The results were externally validated on 150 participants from Specialized Center for Clinically Oriented Research (SCCOR) cohort (aged 45 to 80 years, 52.7% female, 4.7% African Americans; enrollment: 7/2007-12/2012) (AUROC = 0.89). The main limitation of longitudinal studies with 5- and 10-year follow-up is the introduction of mortality bias that disproportionately affects the more severe cases. However, our study focused on spirometrically normal individuals, who have a lower mortality rate. Another limitation is the use of strict criteria to define spirometrically normal individuals, which was unavoidable when studying factors associated with changes in normalized forced expiratory volume in 1 s (FEV1%predicted) or the ratio of FEV1/FVC (forced vital capacity).ConclusionsThis study took an agnostic approach to identify which baseline measurements differentiate and predict the early stages of lung function decline in individuals with previous smoking history. Our analysis suggests that emphysema affects obstruction onset, while airway predominant pathology may play a more important role in future FEV1 (%predicted) decline without obstruction, and FEF25-75% may affect both.

Chronic Obstructive Pulmonary Disease (COPD) is a heterogeneous, chronic inflammatory process of the lungs and, like other complex diseases, is caused by both genetic and environmental factors. Detailed understanding of the molecular mechanisms of complex diseases requires the study of the interplay among different biomolecular layers, and thus the integration of different omics data types. In this study, we investigated COPD-associated molecular mechanisms through a correlation-based network integration of lung tissue RNA-seq and DNA methylation data of COPD cases (n = 446) and controls (n = 346) derived from the Lung Tissue Research Consortium. First, we performed a SWIM-network based analysis to build separate correlation networks for RNA-seq and DNA methylation data for our case-control study population. Then, we developed a method to integrate the results into a coupled network of differentially expressed and differentially methylated genes to investigate their relationships across both molecular layers. The functional enrichment analysis of the nodes of the coupled network revealed a strikingly significant enrichment in Immune System components, both innate and adaptive, as well as immune-system component communication (interleukin and cytokine-cytokine signaling). Our analysis allowed us to reveal novel putative COPD-associated genes and to analyze their relationships, both at the transcriptomics and epigenomics levels, thus contributing to an improved understanding of COPD pathogenesis.

Abstract The consequences of returning infectious pathogen test results identified incidentally in research studies have not been well-studied. Concerns include identification of an important health issue for individuals, accuracy of research test results, public health impact, potential emotional distress for participants, and need for IRB permissions. Blood RNA-sequencing analysis for non-human RNA in 3984 participants from the COPDGene study identified 228 participants with evidence suggestive for hepatitis C virus (HCV) infection. We hypothesized that incidentally discovered HCV results could be effectively returned to COPDGene participants with attention to the identified concerns. In conjunction with a COPDGene Participant Advisory Panel, we developed and obtained IRB approval for a process of returning HCV research results and an HCV Follow-Up Study questionnaire to capture information about previous HCV diagnosis and treatment information and participant reactions to return of HCV results. During phone calls following the initial HCV notification letter, 84 of 124 participants who could be contacted (67.7%) volunteered that they had been previously diagnosed with HCV infection. Thirty-one of these 124 COPDGene participants were enrolled in the HCV Follow-Up Study. Five of the 31 HCV Follow-Up Study participants did not report a previous diagnosis of HCV. For four of these participants, subsequent clinical HCV testing confirmed HCV infection. Thus, 30/31 Follow-Up Study participants had confirmed HCV diagnoses, supporting the accuracy of the HCV research test results. However, the limited number of participants in the Follow-Up Study precludes an accurate assessment of the false-positive and false-negative rates of the research RNA sequencing evidence for HCV. Most HCV Follow-Up Study participants (29/31) were supportive of returning HCV research results, and most participants found the process for returning HCV results to be informative and not upsetting. Newly diagnosed participants were more likely to be pleased to learn about a potentially curable infection (p = 0.027) and showed a trend toward being more frightened by the potential health risks of HCV (p = 0.11). We conclude that HCV results identified incidentally during transcriptomic research studies can be successfully returned to research study participants with a carefully designed process.

Abstract Expression quantitative trait methylation (eQTM) analysis identifies DNA CpG sites at which methylation is associated with gene expression. The present study describes an eQTM resource of CpG-transcript pairs derived from whole blood DNA methylation and RNA sequencing gene expression data in 2115 Framingham Heart Study participants. We identified 70,047 significant cis CpG-transcript pairs at p

ABSTRACTBackground: Recent studies have identified thousands of associations between DNA methylation CpGs and complex diseases/traits, emphasizing the critical role of epigenetics in understanding disease aetiology and identifying biomarkers. However, association analyses based on methylation array data are susceptible to batch/slide effects, which can lead to inflated false positive rates or reduced statistical powerResults: We use multiple DNA methylation datasets based on the popular Illumina Infinium MethylationEPIC BeadChip array to describe consistent patterns and the joint distribution of slide effects across CpGs, confirming and extending previous results. The susceptible CpGs overlap with the Illumina Infinium HumanMethylation450 BeadChip array content.Conclusions: Our findings reveal systematic patterns in slide effects. The observations provide further insights into the characteristics of these effects and can improve existing adjustment approaches.

Summary: Background: Preeclampsia has been associated with maternal epigenetic changes, in particular DNA methylation changes in the placenta. It has been suggested that preeclampsia could also cause DNA methylation changes in the neonate. We examined DNA methylation in relation to gene expression in the cord blood of offspring born to mothers with preeclampsia. Methods: This study included 128 mother-child pairs who participated in the Vitamin D Antenatal Asthma Reduction Trial (VDAART), where assessment of preeclampsia served as secondary outcome. We performed an epigenome-wide association study of preeclampsia and cord blood DNA methylation (Illumina 450 K chip). We then examined gene expression of the same subjects for validation and replicated the gene signatures in independent DNA methylation datasets. Lastly, we applied functional enrichment and network analyses to identify biological pathways that could potentially be involved in preeclampsia. Findings: In the cord blood samples (n = 128), 263 CpGs were differentially methylated (FDR

Abstract Background Chronic obstructive pulmonary disease (COPD) is a highly morbid and heterogenous disease. While COPD is defined by spirometry, many COPD characteristics are seen in cigarette smokers with normal spirometry. The extent to which COPD and COPD heterogeneity is captured in omics of lung tissue is not known. Methods We clustered gene expression and methylation data in 78 lung tissue samples from former smokers with normal lung function or severe COPD. We applied two integrative omics clustering methods: (1) Similarity Network Fusion (SNF) and (2) Entropy-Based Consensus Clustering (ECC). Results SNF clusters were not significantly different by the percentage of COPD cases (48.8% vs. 68.6%, p = 0.13), though were different according to median forced expiratory volume in one second (FEV1) % predicted (82 vs. 31, p = 0.017). In contrast, the ECC clusters showed stronger evidence of separation by COPD case status (48.2% vs. 81.8%, p = 0.013) and similar stratification by median FEV1% predicted (82 vs. 30.5, p = 0.0059). ECC clusters using both gene expression and methylation were identical to the ECC clustering solution generated using methylation data alone. Both methods selected clusters with differentially expressed transcripts enriched for interleukin signaling and immunoregulatory interactions between lymphoid and non-lymphoid cells. Conclusions Unsupervised clustering analysis from integrated gene expression and methylation data in lung tissue resulted in clusters with modest concordance with COPD, though were enriched in pathways potentially contributing to COPD-related pathology and heterogeneity.

Abstract Inference and analysis of gene regulatory networks (GRNs) require software that integrates multi-omic data from various sources. The Network Zoo (netZoo; netzoo.github.io) is a collection of open-source methods to infer GRNs, conduct differential network analyses, estimate community structure, and explore the transitions between biological states. The netZoo builds on our ongoing development of network methods, harmonizing the implementations in various computing languages and between methods to allow better integration of these tools into analytical pipelines. We demonstrate the utility using multi-omic data from the Cancer Cell Line Encyclopedia. We will continue to expand the netZoo to incorporate additional methods.

Summary: Background: Identifying novel epigenetic signatures associated with serum immunoglobulin E (IgE) may improve our understanding of molecular mechanisms underlying asthma and IgE-mediated diseases. Methods: We performed an epigenome-wide association study using whole blood from Framingham Heart Study (FHS; n = 3,471, 46% females) participants and validated results using the Childhood Asthma Management Program (CAMP; n = 674, 39% females) and the Genetic Epidemiology of Asthma in Costa Rica Study (CRA; n = 787, 41% females). Using the closest gene to each IgE-associated CpG, we highlighted biologically plausible pathways underlying IgE regulation and analyzed the transcription patterns linked to IgE-associated CpGs (expression quantitative trait methylation loci; eQTMs). Using prior UK Biobank summary data from genome-wide association studies of asthma and allergy, we performed Mendelian randomization (MR) for causal inference testing using the IgE-associated CpGs from FHS with methylation quantitative trait loci (mQTLs) as instrumental variables. Findings: We identified 490 statistically significant differentially methylated CpGs associated with IgE in FHS, of which 193 (39.3%) replicated in CAMP and CRA (FDR < 0.05). Gene ontology analysis revealed enrichment in pathways related to transcription factor binding, asthma, and other immunological processes. eQTM analysis identified 124 cis-eQTMs for 106 expressed genes (FDR < 0.05). MR in combination with drug-target analysis revealed CTSB and USP20 as putatively causal regulators of IgE levels (Bonferroni adjusted P

Abstract Background Asthma is a heterogeneous disease with high morbidity. Advancement in high-throughput multi-omics approaches has enabled the collection of molecular assessments at different layers, providing a complementary perspective of complex diseases. Numerous computational methods have been developed for the omics-based patient classification or disease outcome prediction. Yet, a systematic benchmarking of those methods using various combinations of omics data for the prediction of asthma development is still lacking. Objective We aimed to investigate the computational methods in disease status prediction using multi-omics data. Method We systematically benchmarked 18 computational methods using all the 63 combinations of six omics data (GWAS, miRNA, mRNA, microbiome, metabolome, DNA methylation) collected in The Vitamin D Antenatal Asthma Reduction Trial (VDAART) cohort. We evaluated each method using standard performance metrics for each of the 63 omics combinations. Results Our results indicate that overall Logistic Regression, Multi-Layer Perceptron, and MOGONET display superior performance, and the combination of transcriptional, genomic and microbiome data achieves the best prediction. Moreover, we find that including the clinical data can further improve the prediction performance for some but not all the omics combinations. Conclusions Specific omics combinations can reach the optimal prediction of asthma development in children. And certain computational methods showed superior performance than other methods.

Abstract Background The association between genetic variants on the X chromosome to risk of COPD has not been fully explored. We hypothesize that the X chromosome harbors variants important in determining risk of COPD related phenotypes and may drive sex differences in COPD manifestations. Methods Using X chromosome data from three COPD-enriched cohorts of adult smokers, we performed X chromosome specific quality control, imputation, and testing for association with COPD case–control status, lung function, and quantitative emphysema. Analyses were performed among all subjects, then stratified by sex, and subsequently combined in meta-analyses. Results Among 10,193 subjects of non-Hispanic white or European ancestry, a variant near TMSB4X, rs5979771, reached genome-wide significance for association with lung function measured by FEV1/FVC ( $$\beta$$ β 0.020, SE 0.004, p 4.97 × 10–08), with suggestive evidence of association with FEV1 ( $$\beta$$ β 0.092, SE 0.018, p 3.40 × 10–07). Sex-stratified analyses revealed X chromosome variants that were differentially trending in one sex, with significantly different effect sizes or directions. Conclusions This investigation identified loci influencing lung function, COPD, and emphysema in a comprehensive genetic association meta-analysis of X chromosome genetic markers from multiple COPD-related datasets. Sex differences play an important role in the pathobiology of complex lung disease, including X chromosome variants that demonstrate differential effects by sex and variants that may be relevant through escape from X chromosome inactivation. Comprehensive interrogation of the X chromosome to better understand genetic control of COPD and lung function is important to further understanding of disease pathology. Trial registration Genetic Epidemiology of COPD Study (COPDGene) is registered at ClinicalTrials.gov, NCT00608764 (Active since January 28, 2008). Evaluation of COPD Longitudinally to Identify Predictive Surrogate Endpoints Study (ECLIPSE), GlaxoSmithKline study code SCO104960, is registered at ClinicalTrials.gov, NCT00292552 (Active since February 16, 2006). Genetics of COPD in Norway Study (GenKOLS) holds GlaxoSmithKline study code RES11080, Genetics of Chronic Obstructive Lung Disease.

Cigarette smoking is associated with epigenetic changes that may be reversible following smoking cessation. Whole blood DNA methylation was evaluated in Framingham Heart Study Offspring (n = 169) and Third Generation (n = 30) cohort participants at two study visits 6 years apart and in Atherosclerosis Risk in Communities (ARIC) study (n = 222) participants at two study visits 20 years apart. Changes in DNA methylation (delta β values) at 483,565 cytosine-phosphate-guanine (CpG) sites and differentially methylated regions (DMRs) were compared between participants who were current, former, or never smokers at both visits (current-current, former-former, never-never, respectively), versus those who quit in the interim (current-former). Interim quitters had more hypermethylation at four CpGs annotated to AHRR, one CpG annotated to F2RL3, and one intergenic CpG (cg21566642) compared with current-current smokers (FDR < 0.02 for all), and two significant DMRs were identified. While there were no significant differentially methylated CpGs in the comparison of interim quitters and former-former smokers, 106 DMRs overlapping with small nucleolar RNA were identified. As compared with all non-smokers, current-current smokers additionally had more hypermethylation at two CpG sites annotated to HIVEP3 and TMEM126A, respectively, and another intergenic CpG (cg14339116). Gene transcripts associated with smoking cessation were implicated in immune responses, cell homoeostasis, and apoptosis. Smoking cessation is associated with early reversion of blood DNA methylation changes at CpG sites annotated to AHRR and F2RL3 towards those of never smokers. Associated gene expression suggests a role of longitudinal smoking-related DNA methylation changes in immune response processes.

Abstract Background Chronic obstructive pulmonary disease (COPD) and idiopathic pulmonary fibrosis (IPF) are characterized by shared exposures and clinical features, but distinct genetic and pathologic features exist. These features have not been well-studied using large-scale gene expression datasets. We hypothesized that there are divergent gene, pathway, and cellular signatures between COPD and IPF. Methods We performed RNA-sequencing on lung tissues from individuals with IPF (n = 231) and COPD (n = 377) compared to control (n = 267), defined as individuals with normal spirometry. We grouped the overlapping differential expression gene sets based on direction of expression and examined the resultant sets for genes of interest, pathway enrichment, and cell composition. Using gene set variation analysis, we validated the overlap group gene sets in independent COPD and IPF data sets. Results We found 5010 genes differentially expressed between COPD and control, and 11,454 genes differentially expressed between IPF and control (1% false discovery rate). 3846 genes overlapped between IPF and COPD. Several pathways were enriched for genes upregulated in COPD and downregulated in IPF; however, no pathways were enriched for genes downregulated in COPD and upregulated in IPF. There were many myeloid cell genes with increased expression in COPD but decreased in IPF. We found that the genes upregulated in COPD but downregulated in IPF were associated with lower lung function in the independent validation cohorts. Conclusions We identified a divergent gene expression signature between COPD and IPF, with increased expression in COPD and decreased in IPF. This signature is associated with worse lung function in both COPD and IPF.

Abstract Background Chronic obstructive pulmonary disease (COPD) is a leading cause of death in adults that may have origins in early lung development. It is a complex disease, influenced by multiple factors including genetic variants and environmental factors. Maternal smoking during pregnancy may influence the risk for diseases during adulthood, potentially through epigenetic modifications including methylation. Methods In this work, we explore the fetal origins of COPD by utilizing lung DNA methylation marks associated with in utero smoke (IUS) exposure, and evaluate the network relationships between methylomic and transcriptomic signatures associated with adult lung tissue from former smokers with and without COPD. To identify potential pathobiological mechanisms that may link fetal lung, smoke exposure and adult lung disease, we study the interactions (physical and functional) of identified genes using protein–protein interaction networks. Results We build IUS-exposure and COPD modules, which identify connected subnetworks linking fetal lung smoke exposure to adult COPD. Studying the relationships and connectivity among the different modules for fetal smoke exposure and adult COPD, we identify enriched pathways, including the AGE-RAGE and focal adhesion pathways. Conclusions The modules identified in our analysis add new and potentially important insights to understanding the early life molecular perturbations related to the pathogenesis of COPD. We identify AGE-RAGE and focal adhesion as two biologically plausible pathways that may reveal lung developmental contributions to COPD. We were not only able to identify meaningful modules but were also able to study interconnections between smoke exposure and lung disease, augmenting our knowledge about the fetal origins of COPD.

Abstract Background Optimism is the general belief that good things will occur in the future; optimism is modifiable by cognitive behavioral therapy (CBT). Previous studies have associated higher optimism with improved health outcomes and lower all-cause mortality. Research question Investigate association between optimism and disease-related characteristics in chronic obstructive pulmonary disease (COPD). Study design and methods Current and former smokers with/without COPD and Preserved Ratio Impaired Spirometry (PRISm) from the 10-year follow-up visit for the Genetic Epidemiology of COPD (COPDGene) study were included. Optimism was assessed at the 10-year visit using the Life Orientation Test-Revised. Models of optimism as a predictor of lung function, COPD-associated phenotypes including exacerbations, and functional assessments, were adjusted for demographic confounders, smoking status, and comorbidities. Results Among 1967 subjects, higher optimism was significantly associated with older age, non-Hispanic white race, marital status, quitting smoking status, absence of COPD, and absence of depression. In multivariable analysis, higher optimism was independently associated with fewer prior exacerbations of COPD (coef = − 0.037, P

Abstract The biological processes that drive cellular function can be represented by a complex network of interactions between regulators (transcription factors) and their targets (genes). A cell’s epigenetic state plays an important role in mediating these interactions, primarily by influencing chromatin accessibility. However, how to effectively use epigenetic data when constructing a gene regulatory network remains an open question. Almost all existing network reconstruction approaches focus on estimating transcription factor to gene connections using transcriptomic data. In contrast, computational approaches for analyzing epigenetic data generally focus on improving transcription factor binding site predictions rather than deducing regulatory network relationships. We bridged this gap by developing SPIDER, a network reconstruction approach that incorporates epigenetic data into a message-passing framework to estimate gene regulatory networks. We validated SPIDER’s predictions using ChIP-seq data from ENCODE and found that SPIDER networks are both highly accurate and include cell-line-specific regulatory interactions. Notably, SPIDER can recover ChIP-seq verified transcription factor binding events in the regulatory regions of genes that do not have a corresponding sequence motif. The networks estimated by SPIDER have the potential to identify novel hypotheses that will allow us to better characterize cell-type and phenotype specific regulatory mechanisms.

Abstract Background Chronic obstructive pulmonary disease (COPD) is the third leading cause of death globally. COPD patients with cachexia or weight loss have increased risk of death independent of body mass index (BMI) and lung function. We tested the hypothesis genetic variation is associated with weight loss in COPD using a genome‐wide association study approach. Methods Participants with COPD (N = 4308) from three studies (COPDGene, ECLIPSE, and SPIROMICS) were analysed. Discovery analyses were performed in COPDGene with replication in SPIROMICS and ECLIPSE. In COPDGene, weight loss was defined as self‐reported unintentional weight loss > 5% in the past year or low BMI (BMI

Abstract Objective Illumina BeadChip arrays are commonly used to generate DNA methylation data for large epidemiological studies. Updates in technology over time create challenges for data harmonization within and between studies, many of which obtained data from the older 450K and newer EPIC platforms. The pre-processing pipeline for DNA methylation is not trivial, and influences the downstream analyses. Incorporating different platforms adds a new level of technical variability that has not yet been taken into account by recommended pipelines. Our study evaluated the performance of various tools on different versions of platform data harmonization at each step of pre-processing pipeline, including quality control (QC), normalization, batch effect adjustment, and genomic inflation. We illustrate our novel approach using 450K and EPIC data from the Diabetes Autoimmunity Study in the Young (DAISY) prospective cohort. Results We found normalization and probe filtering had the biggest effect on data harmonization. Employing a meta-analysis was an effective and easily executable method for accounting for platform variability. Correcting for genomic inflation also helped with harmonization. We present guidelines for studies seeking to harmonize data from the 450K and EPIC platforms, which includes the use of technical replicates for evaluating numerous pre-processing steps, and employing a meta-analysis.

The identification and understanding of gene-environment interactions can provide insights into the pathways and mechanisms underlying complex diseases. However, testing for gene-environment interaction remains a challenge since a.) statistical power is often limited and b.) modeling of environmental effects is nontrivial and such model misspecifications can lead to false positive interaction findings. To address the lack of statistical power, recent methods aim to identify interactions on an aggregated level using, for example, polygenic risk scores. While this strategy can increase the power to detect interactions, identifying contributing genes and pathways is difficult based on these relatively global results. Here, we propose RITSS (Robust Interaction Testing using Sample Splitting), a gene-environment interaction testing framework for quantitative traits that is based on sample splitting and robust test statistics. RITSS can incorporate sets of genetic variants and/or multiple environmental factors. Based on the user's choice of statistical/machine learning approaches, a screening step selects and combines potential interactions into scores with improved interpretability. In the testing step, the application of robust statistics minimizes the susceptibility to main effect misspecifications. Using extensive simulation studies, we demonstrate that RITSS controls the type 1 error rate in a wide range of scenarios, and we show how the screening strategy influences statistical power. In an application to lung function phenotypes and human height in the UK Biobank, RITSS identified highly significant interactions based on subcomponents of genetic risk scores. While the contributing single variant interaction signals are weak, our results indicate interaction patterns that result in strong aggregated effects, providing potential insights into underlying gene-environment interaction mechanisms.

Abstract Background Secondary polycythemia is associated with cigarette smoking and chronic obstructive pulmonary disease (COPD). However, the prevalence of polycythemia in COPD and the contributing risk factors for polycythemia in COPD have not been extensively studied. Methods We analyzed the presence of secondary polycythemia in current and former smokers with moderate to very severe COPD at the five-year follow-up visit in the observational COPDGene study. We used logistic regression to evaluate the association of polycythemia with age, sex, race, altitude, current smoking status, spirometry, diffusing capacity for carbon monoxide (DLCO), quantitative chest CT measurements (including emphysema, airway wall thickness, and pulmonary artery to aorta diameter ratio), resting hypoxemia, exercise-induced hypoxemia, and long-term oxygen therapy. Results In a total of 1928 COPDGene participants with moderate to very severe COPD, secondary polycythemia was found in 97 (9.2%) male and 31 (3.5%) female participants. In a multivariable logistic model, severe resting hypoxemia (OR 3.50, 95% CI 1.41–8.66), impaired DLCO (OR 1.28 for each 10-percent decrease in DLCO % predicted, CI 1.09–1.49), male sex (OR 3.60, CI 2.20–5.90), non-Hispanic white race (OR 3.33, CI 1.71–6.50), current smoking (OR 2.55, CI 1.49–4.38), and enrollment in the Denver clinical center (OR 4.42, CI 2.38–8.21) were associated with higher risk for polycythemia. In addition, continuous (OR 0.13, CI 0.05–0.35) and nocturnal (OR 0.46, CI 0.21–0.97) supplemental oxygen were associated with lower risk for polycythemia. Results were similar after excluding participants with anemia and participants enrolled at the Denver clinical center. Conclusions In a large cohort of individuals with moderate to very severe COPD, male sex, current smoking, enrollment at the Denver clinical center, impaired DLCO, and severe hypoxemia were associated with increased risk for secondary polycythemia. Continuous or nocturnal supplemental oxygen use were associated with decreased risk for polycythemia.

Cigarette smoking impacts DNA methylation, but the investigation of sex-specific features of lung tissue DNA methylation in smokers has been limited. Women appear more susceptible to cigarette smoke, and often develop more severe lung disease at an earlier age with less smoke exposure. We aimed to analyse whether there are sex differences in DNA methylation in lung tissue and whether these DNA methylation marks interact with smoking. We collected lung tissue samples from former smokers who underwent lung tissue resection. One hundred thirty samples from white subjects were included for this analysis. Regression models for sex as a predictor of methylation were adjusted for age, presence of COPD, smoking variables and technical batch variables revealed 710 associated sites. 294 sites demonstrated robust sex-specific methylation associations in foetal lung tissue. Pathway analysis identified 6 nominally significant pathways including the mitophagy pathway. Three CpG sites demonstrated a suggested interaction between sex and pack-years of smoking: GPR132, ANKRD44 and C19orf60. All of them were nominally significant in both male- and female-specific models, and the effect estimates were in opposite directions for male and female; GPR132 demonstrated significant association between DNA methylation and gene expression in lung tissue (P

Abstract 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

Abstract Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) has led to the global coronavirus disease 2019 (COVID-19) pandemic. SARS-CoV-2 enters cells via angiotensin-Converting Enzyme 2 (ACE2) receptors, highly expressed in nasal epithelium with parallel high infectivity.1,2 The nasal epigenome is in direct contact with the environment and could explain COVID-19 disparities by reflecting social and environmental influences on ACE2 regulation. We collected nasal swabs from anterior nares of 547 children, measured DNA methylation (DNAm), and tested differences at 15 ACE2 CpGs by sex, age, race/ethnicity and epigenetic age. ACE2 CpGs were differentially methylated by sex with 12 sites having lower DNAm (mean = 12.71%) and 3 sites greater DNAm (mean = 1.45%) among females relative to males. We observed differential DNAm at 5 CpGs for Hispanic females (mean absolute difference = 3.22%) and lower DNAm at 8 CpGs for Black males (mean absolute difference = 1.33%), relative to white participants. Longer DNAm telomere length was associated with greater ACE2 DNAm at 11 and 13 CpGs among males (mean absolute difference = 7.86%) and females (mean absolute difference = 8.21%), respectively. Nasal ACE2 DNAm differences could contribute to our understanding COVID-19 severity and disparities reflecting upstream environmental and social influences. Findings need to be confirmed among adults and patients with risk factors for COVID-19 severity.

Background: Sex-specific differences in fetal lung maturation have been well described; however, little is known about the sex-specific differences in microRNA (miRNA) expression during human fetal lung development. Interestingly, many adult chronic lung diseases also demonstrate sex-specific differences in prevalence. The developmental origins of health and disease hypothesis suggests that these sex-specific differences in fetal lung development may influence disease susceptibility later in life. In this study, we performed miRNA sequencing on human fetal lung tissue samples to investigate differential expression of miRNAs between males and females in the pseudoglandular stage of lung development. We hypothesized that differences in miRNA expression are present between sexes in early human lung development and may contribute to the sex-specific differences seen in pulmonary diseases later in life.Methods: RNA was isolated from human fetal lung tissue samples for miRNA sequencing. The count of each miRNA was modeled by sex using negative binomial regression models in DESeq2, adjusting for post-conception age, age2, smoke exposure, batch, and RUV factors. We tested for differential expression of miRNAs by sex, and for the presence of sex-by-age interactions to determine if miRNA expression levels by age were distinct between males and females.Results: miRNA expression profiles were generated on 298 samples (166 males and 132 females). Of the 809 miRNAs expressed in human fetal lung tissue during the pseudoglandular stage of lung development, we identified 93 autosomal miRNAs that were significantly differentially expressed by sex and 129 miRNAs with a sex-specific pattern of miRNA expression across the course of the pseudoglandular period.Conclusion: Our study demonstrates differential expression of numerous autosomal miRNAs between the male and female developing human lung. Additionally, the expression of some miRNAs are modified by age across the pseudoglandular stage in a sex-specific way. Some of these differences in miRNA expression may impact susceptibility to pulmonary disease later in life. Our results suggest that sex-specific miRNA expression during human lung development may be a potential mechanism to explain sex-specific differences in lung development and may impact subsequent disease susceptibility.

Evidence suggests there are roles for vitamin K in various chronic disease outcomes, but population-level diet and supplement recommendations are difficult to determine due to high levels of variability in measures of status and response to intake compared to other nutrients. In this preliminary investigation, a blood-based epigenome-wide association study (EWAS) comparing responders and non-responders to phylloquinone (vitamin K1) supplementation (NCT00183001) was undertaken in order to better understand the molecular underpinnings of this observed variability. Responders (n = 24) and non-responders (n = 24) were identified in a prior 3-year phylloquinone supplementation trial based on their changes in plasma phylloquinone concentrations. Differential DNA methylation was identified in multiple regions with previously unknown relationships to phylloquinone absorption and metabolism, such as at the TMEM263 locus. A hypothesis-driven analysis of lipid-related genes highlighted a site in the NPC1L1 gene, supplementing existing evidence for its role in phylloquinone absorption. Furthermore, an EWAS for baseline plasma phylloquinone concentrations revealed a strong correlation between the epigenomic signatures of phylloquinone baseline status and response to supplementation. This work can guide future epigenomic research on vitamin K and contributes to the development of more personalized dietary recommendations for vitamin K.

Abstract Background Preterm birth and shorter duration of pregnancy are associated with increased morbidity in neonatal and later life. As the epigenome is known to have an important role during fetal development, we investigated associations between gestational age and blood DNA methylation in children. Methods We performed meta-analysis of Illumina’s HumanMethylation450-array associations between gestational age and cord blood DNA methylation in 3648 newborns from 17 cohorts without common pregnancy complications, induced delivery or caesarean section. We also explored associations of gestational age with DNA methylation measured at 4–18 years in additional pediatric cohorts. Follow-up analyses of DNA methylation and gene expression correlations were performed in cord blood. DNA methylation profiles were also explored in tissues relevant for gestational age health effects: fetal brain and lung. Results We identified 8899 CpGs in cord blood that were associated with gestational age (range 27–42 weeks), at Bonferroni significance, P

Summary: Genetic studies on telomere length are important for understanding age-related diseases. Prior GWASs for leukocyte TL have been limited to European and Asian populations. Here, we report the first sequencing-based association study for TL across ancestrally diverse individuals (European, African, Asian, and Hispanic/Latino) from the NHLBI Trans-Omics for Precision Medicine (TOPMed) program. We used whole-genome sequencing (WGS) of whole blood for variant genotype calling and the bioinformatic estimation of telomere length in n = 109,122 individuals. We identified 59 sentinel variants (p < 5 × 10−9) in 36 loci associated with telomere length, including 20 newly associated loci (13 were replicated in external datasets). There was little evidence of effect size heterogeneity across populations. Fine-mapping at OBFC1 indicated that the independent signals colocalized with cell-type-specific eQTLs for OBFC1 (STN1). Using a multi-variant gene-based approach, we identified two genes newly implicated in telomere length, DCLRE1B (SNM1B) and PARN. In PheWAS, we demonstrated that our TL polygenic trait scores (PTSs) were associated with an increased risk of cancer-related phenotypes.

To increase power and minimize bias in statistical analyses, quantitative outcomes are often adjusted for precision and confounding variables using standard regression approaches. The outcome is modeled as a linear function of the precision variables and confounders; however, for many complex phenotypes, the assumptions of the linear regression models are not always met. As an alternative, we used neural networks for the modeling of complex phenotypes and covariate adjustments. We compared the prediction accuracy of the neural network models to that of classical approaches based on linear regression. Using data from the UK Biobank, COPDGene study, and Childhood Asthma Management Program (CAMP), we examined the features of neural networks in this context and compared them with traditional regression approaches for prediction of three outcomes: forced expiratory volume in one second (FEV1), age at smoking cessation, and log transformation of age at smoking cessation (due to age at smoking cessation being right-skewed). We used mean squared error to compare neural network and regression models, and found the models performed similarly unless the observed distribution of the phenotype was skewed, in which case the neural network had smaller mean squared error. Our results suggest neural network models have an advantage over standard regression approaches when the phenotypic distribution is skewed. However, when the distribution is not skewed, the approaches performed similarly. Our findings are relevant to studies that analyze phenotypes that are skewed by nature or where the phenotype of interest is skewed as a result of the ascertainment condition.

Rationale: PM2.5-induced adverse effects on respiratory health may be driven by epigenetic modifications in airway cells. The potential impact of exposure duration on epigenetic alterations in the airways is not yet known. Objectives: We aimed to study associations of fine particulate matter PM2.5 exposure with DNA methylation in nasal cells. Methods: We conducted nasal epigenome-wide association analyses within 503 children from Project Viva (mean age 12.9 y), and examined various exposure durations (1-day, 1-week, 1-month, 3-months and 1-year) prior to nasal sampling. We used residential addresses to estimate average daily PM2.5 at 1 km resolution. We collected nasal swabs from the anterior nares and measured DNA methylation (DNAm) using the Illumina MethylationEPIC BeadChip. We tested 719,075 high quality autosomal CpGs using CpG-by-CpG and regional DNAm analyses controlling for multiple comparisons, and adjusted for maternal education, household smokers, child sex, race/ethnicity, BMI z-score, age, season at sample collection and cell-type heterogeneity. We further corrected for bias and genomic inflation. We tested for replication in a cohort from the Netherlands (PIAMA). Results: In adjusted analyses, we found 362 CpGs associated with 1-year PM2.5 (FDR

Epigenetic differences in nasal epithelium have been proposed as a biomarker for lower airway disease and asthma. Here, in epigenome-wide association studies for asthma and other airway traits using nasal swabs, the authors identify differentially methylated CpGs that highlight genes involved in TH2 response.

Epigenetic mechanisms integrate both genetic variability and environmental exposures. However, comprehensive epigenome-wide analysis has not been performed across major childhood allergic phenotypes. We examined the association of epigenome-wide DNA methylation in mid-childhood peripheral blood (Illumina HumanMethyl450K) with mid-childhood atopic sensitization, environmental/inhalant and food allergen sensitization in 739 children in two birth cohorts (Project Viva–Boston, and the Generation R Study–Rotterdam). We performed covariate-adjusted epigenome-wide association meta-analysis and employed pathway and regional analyses of results. Seven-hundred and five methylation sites (505 genes) were significantly cross-sectionally associated with mid-childhood atopic sensitization, 1411 (905 genes) for environmental and 45 (36 genes) for food allergen sensitization (FDR

Birthweight has been found to associate with later-life health outcomes. Here the authors perform a meta-analysis of epigenome-wide association studies of 8,825 neonates from 24 birth cohorts in the Pregnancy And Childhood Epigenetics Consortium, identifying differentially methylated CpGs in neonatal blood that associate with birthweight.

Abstract Background Multiple gene expression studies have been performed separately in peripheral blood, lung, and airway tissues to study COPD. We performed RNA-sequencing gene expression profiling of large-airway epithelium, alveolar macrophage and peripheral blood samples from the same subset of COPD cases and controls from the COPDGene study who underwent bronchoscopy at a single center. Using statistical and gene set enrichment approaches, we sought to improve the understanding of COPD by studying gene sets and pathways across these tissues, beyond the individual genomic determinants. Methods We performed differential expression analysis using RNA-seq data obtained from 63 samples from 21 COPD cases and controls (includes four non-smokers) via the R package DESeq2. We tested associations between gene expression and variables related to lung function, smoking history, and CT scan measures of emphysema and airway disease. We examined the correlation of differential gene expression across the tissues and phenotypes, hypothesizing that this would reveal preserved and private gene expression signatures. We performed gene set enrichment analyses using curated databases and findings from prior COPD studies to provide biological and disease relevance. Results The known smoking-related genes CYP1B1 and AHRR were among the top differential expression results for smoking status in the large-airway epithelium data. We observed a significant overlap of genes primarily across large-airway and macrophage results for smoking and airway disease phenotypes. We did not observe specific genes differentially expressed in all three tissues for any of the phenotypes. However, we did observe hemostasis and immune signaling pathways in the overlaps across all three tissues for emphysema, and amyloid and telomere-related pathways for smoking. In peripheral blood, the emphysema results were enriched for B cell related genes previously identified in lung tissue studies. Conclusions Our integrative analyses across COPD-relevant tissues and prior studies revealed shared and tissue-specific disease biology. These replicated and novel findings in the airway and peripheral blood have highlighted candidate genes and pathways for COPD pathogenesis.

Abstract Background Maternal mood disorders and their treatment during pregnancy may have effects on the offspring epigenome. We aim to evaluate associations of maternal prenatal antidepressant use, anxiety, and depression with cord blood DNA methylation across the genome at birth and test for persistence of associations in early and mid-childhood blood DNA. Methods A discovery phase was conducted in Project Viva, a prospective pre-birth cohort study with external replication in an independent cohort, the Generation R Study. In Project Viva, pregnant women were recruited between 1999 and 2002 in Eastern Massachusetts, USA. In the Generation R Study, pregnant women were recruited between 2002 and 2006 in Rotterdam, the Netherlands. In Project Viva, 479 infants had data on maternal antidepressant use, anxiety, depression, and cord blood DNA methylation, 120 children had DNA methylation measured in early childhood (~ 3 years), and 460 in mid-childhood (~ 7 years). In the Generation R Study, 999 infants had data on maternal antidepressants and cord blood DNA methylation. The prenatal antidepressant prescription was obtained from medical records. At-mid pregnancy, symptoms of anxiety and depression were assessed with the Pregnancy-Related Anxiety Scale and the Edinburgh Postnatal Depression Scale in Project Viva and with the Brief Symptom Inventory in the Generation R Study. Genome-wide DNA methylation was measured using the Infinium HumanMethylation450 BeadChip in both cohorts. Results In Project Viva, 2.9% (14/479) pregnant women were prescribed antidepressants, 9.0% (40/445) experienced high pregnancy-related anxiety, and 8.2% (33/402) reported symptoms consistent with depression. Newborns exposed to antidepressants in pregnancy had 7.2% lower DNA methylation (95% CI, − 10.4, − 4.1; P = 1.03 × 10−8) at cg22159528 located in the gene body of ZNF575, and this association replicated in the Generation R Study (β = − 2.5%; 95% CI − 4.2, − 0.7; P = 0.006). In Project Viva, the association persisted in early (β = − 6.2%; 95% CI − 10.7, − 1.6) but not mid-childhood. We observed cohort-specific associations for maternal anxiety and depression in Project Viva that did not replicate. Conclusions The ZNF575 gene is involved in transcriptional regulation but specific functions are largely unknown. Given the widespread use of antidepressants in pregnancy, as well as the effects of exposure to anxiety and depression, implications of potential fetal epigenetic programming by these risk factors and their impacts on development merit further investigation.

Despite their known importance in clinical medicine, differences based on sex and gender are among the least studied factors affecting cancer susceptibility, progression, survival, and therapeutic response. In particular, the molecular mechanisms driving sex differences are poorly understood and so most approaches to precision medicine use mutational or other genomic data to assign therapy without considering how the sex of the individual might influence therapeutic efficacy. The mandate by the National Institutes of Health that research studies include sex as a biological variable has begun to expand our understanding on its importance. Sex differences in cancer may arise due to a combination of environmental, genetic, and epigenetic factors, as well as differences in gene regulation, and expression. Extensive sex differences occur genome-wide, and ultimately influence cancer biology and outcomes. In this review, we summarize the current state of knowledge about sex-specific genetic and genome-wide influences in cancer, describe how differences in response to environmental exposures and genetic and epigenetic alterations alter the trajectory of the disease, and provide insights into the importance of integrative analyses in understanding the interplay of sex and genomics in cancer. In particular, we will explore some of the emerging analytical approaches, such as the use of network methods, that are providing a deeper understanding of the drivers of differences based on sex and gender. Better understanding these complex factors and their interactions will improve cancer prevention, treatment, and outcomes for all individuals.

Background Epigenome‐wide association studies for cardiometabolic risk factors have discovered multiple loci associated with incident cardiovascular disease (CVD). However, few studies have sought to directly optimize a predictor of CVD risk. Furthermore, it is challenging to train multivariate models across multiple studies in the presence of study‐ or batch effects. Methods and Results Here, we analyzed existing DNA methylation data collected using the Illumina HumanMethylation450 microarray to create a predictor of CVD risk across 3 cohorts: Women's Health Initiative, Framingham Heart Study Offspring Cohort, and Lothian Birth Cohorts. We trained Cox proportional hazards‐based elastic net regressions for incident CVD separately in each cohort and used a recently introduced cross‐study learning approach to integrate these individual scores into an ensemble predictor. The methylation‐based risk score was associated with CVD time‐to‐event in a held‐out fraction of the Framingham data set (hazard ratio per SD=1.28, 95% CI, 1.10–1.50) and predicted myocardial infarction status in the independent REGICOR (Girona Heart Registry) data set (odds ratio per SD=2.14, 95% CI, 1.58–2.89). These associations remained after adjustment for traditional cardiovascular risk factors and were similar to those from elastic net models trained on a directly merged data set. Additionally, we investigated interactions between the methylation‐based risk score and both genetic and biochemical CVD risk, showing preliminary evidence of an enhanced performance in those with less traditional risk factor elevation. Conclusions This investigation provides proof‐of‐concept for a genome‐wide, CVD‐specific epigenomic risk score and suggests that DNA methylation data may enable the discovery of high‐risk individuals who would be missed by alternative risk metrics.

DNA methylation is an epigenetic regulator of gene transcription, which has been found to be both metastable and variable within human cohort studies. Currently, few studies have been done to identify metastable DNA methylation biomarkers associated with longitudinal lung function decline in humans. The identification of such biomarkers is important for screening vulnerable populations. We hypothesized that quantifiable blood-based DNA methylation alterations would serve as metastable biomarkers of lung function decline and aging, which may help to discover new pathways and/or mechanisms related to pulmonary pathogenesis. Using linear mixed models, we performed an Epigenome Wide Association Study (EWAS) between DNA methylation at CpG dinucleotides and longitudinal lung function (FVC, FEV1, FEF25–75%) decline and aging with initial discovery in the Normative Aging Study, and replication in the Cooperative Health Research in the Region of Augsburg cohort. We identified two metastable epigenetic loci associated with either poor lung function and aging, cg05575921 (AHRR gene), or lung function independently of aging, cg06126421 (IER3 gene). These loci may inform basic mechanisms associated with pulmonary function, pathogenesis, and aging. Human epigenomic variation, may help explain features of lung function decline and related pathophysiology not attributable to DNA sequence alone, such as accelerated pulmonary decline in smokers, former smokers, and perhaps non-smokers. Our EWAS across two cohorts, therefore, will likely have implications for the human population, not just the elderly.

Obesity is associated with higher cardio-metabolic risk even in childhood and adolescence; whether this association is mediated by epigenetic mechanisms remains unclear. We examined the extent to which mid-childhood body mass index (BMI) z-score (median age 7.7 years) was associated with cardio-metabolic risk score in early adolescence (median age 12.9 years) via mid-childhood DNA methylation among 265 children in the Project Viva. We measured DNA methylation in leukocytes using the Infinium Human Methylation450K BeadChip. We assessed mediation CpG-by-CpG using epigenome-wide association analyses, high-dimensional mediation analysis, and natural effect models. We observed mediation by mid-childhood DNA methylation at 6 CpGs for the association between mid-childhood BMI z-score and cardio-metabolic risk score in early adolescence in the high-dimensional mediation analysis (accounting for 10% of the total effect) and in the natural effect model (β = 0.04, P = 3.2e-2, accounting for 13% of the total effect). The natural direct effect of BMI z-score on cardio-metabolic risk score was still evident (β = 0.27, P = 1.1e-25). We also observed mediation by mid-childhood DNA methylation at 5 CpGs that was in the opposite direction from the total effect (natural effect model: β = −0.04, P = 2.0e-2). Mediation in different directions implies a complex role of DNA methylation in the association between BMI and cardio-metabolic risk and needs further investigation. Future studies with larger sample size and greater variability in cardio-metabolic risk will further help elucidate the role of DNA methylation for cardio-metabolic risk.

Abstract Background IgE-mediated sensitization may be epigenetically programmed in utero, but early childhood environment may further alter complex traits and disease phenotypes through epigenetic plasticity. However, the epigenomic footprint underpinning IgE-mediated type-I hypersensitivity has not been well-understood, especially under a longitudinal early-childhood life-course framework. Methods We used epigenome-wide DNA methylation (IlluminaHumanMethylation450 BeadChip) in cord blood and mid-childhood peripheral blood to investigate pre- and post-natal methylation marks associated with mid-childhood (age 6.7–10.2) total serum IgE levels in 217 mother-child pairs in Project Viva—a prospective longitudinal pre-birth cohort in eastern Massachusetts, USA. We identified methylation sites associated with IgE using covariate-adjusted robust linear regressions. Results Nineteen methylation marks in cord blood were associated with IgE in mid-childhood (FDR

Abstract Background Genome-wide association studies (GWAS) have identified single nucleotide polymorphisms (SNPs) significantly associated with chronic obstructive pulmonary disease (COPD). However, many genetic variants show suggestive evidence for association but do not meet the strict threshold for genome-wide significance. Integrative analysis of multiple omics datasets has the potential to identify novel genes involved in disease pathogenesis by leveraging these variants in a functional, regulatory context. Results We performed expression quantitative trait locus (eQTL) analysis using genome-wide SNP genotyping and gene expression profiling of lung tissue samples from 86 COPD cases and 31 controls, testing for SNPs associated with gene expression levels. These results were integrated with a prior COPD GWAS using an ensemble statistical and network methods approach to identify relevant genes and observe them in the context of overall genetic control of gene expression to highlight co-regulated genes and disease pathways. We identified 250,312 unique SNPs and 4997 genes in the cis(local)-eQTL analysis (5% false discovery rate). The top gene from the integrative analysis was MAPT, a gene recently identified in an independent GWAS of lung function. The genes HNRNPAB and PCBP2 with RNA binding activity and the gene ACVR1B were identified in network communities with validated disease relevance. Conclusions The integration of lung tissue gene expression with genome-wide SNP genotyping and subsequent intersection with prior GWAS and omics studies highlighted candidate genes within COPD loci and in communities harboring known COPD genes. This integration also identified novel disease genes in sub-threshold regions that would otherwise have been missed through GWAS.

Abstract Background Blood biomarkers are increasingly used to stratify high risk chronic obstructive pulmonary disease (COPD) patients; however, there are fewer studies that have investigated multiple biomarkers and replicated in multiple large well-characterized cohorts of susceptible current and former smokers. Methods We used two MSD multiplex panels to measure 9 cytokines and chemokines in 2123 subjects from COPDGene and 1117 subjects from SPIROMICS. These biomarkers included: interleukin (IL)-2, IL-6, IL-8, IL-10, tumor necrosis factor (TNF)-α, interferon (IFN)-γ, eotaxin/CCL-11, eotaxin-3/CCL-26, and thymus and activation-regulated chemokine (TARC)/CCL-17. Regression models adjusted for clinical covariates were used to determine which biomarkers were associated with the following COPD phenotypes: airflow obstruction (forced expiratory flow at 1 s (FEV1%) and FEV1/forced vital capacity (FEV1/FVC), chronic bronchitis, COPD exacerbations, and emphysema. Biomarker-genotype associations were assessed by genome-wide association of single nucleotide polymorphisms (SNPs). Results Eotaxin and IL-6 were strongly associated with airflow obstruction and accounted for 3–5% of the measurement variance on top of clinical variables. IL-6 was associated with progressive airflow obstruction over 5 years and both IL-6 and IL-8 were associated with progressive emphysema over 5 years. None of the biomarkers were consistently associated with chronic bronchitis or COPD exacerbations. We identified one novel SNP (rs9302690 SNP) that was associated with CCL17 plasma measurements. Conclusion When assessing smoking related pulmonary disease, biomarkers of inflammation such as IL-2, IL-6, IL-8, and eotaxin may add additional modest predictive value on top of clinical variables alone. Trial registration COPDGene (ClinicalTrials.gov Identifier: NCT02445183 ). Subpopulations and Intermediate Outcomes Measures in COPD Study (SPIROMICS) ( ClinicalTrials.gov Identifier: NCT 01969344 ).

Abstract Background Cell lines are an indispensable tool in biomedical research and often used as surrogates for tissues. Although there are recognized important cellular and transcriptomic differences between cell lines and tissues, a systematic overview of the differences between the regulatory processes of a cell line and those of its tissue of origin has not been conducted. The RNA-Seq data generated by the GTEx project is the first available data resource in which it is possible to perform a large-scale transcriptional and regulatory network analysis comparing cell lines with their tissues of origin. Results We compared 127 paired Epstein-Barr virus transformed lymphoblastoid cell lines (LCLs) and whole blood samples, and 244 paired primary fibroblast cell lines and skin samples. While gene expression analysis confirms that these cell lines carry the expression signatures of their primary tissues, albeit at reduced levels, network analysis indicates that expression changes are the cumulative result of many previously unreported alterations in transcription factor (TF) regulation. More specifically, cell cycle genes are over-expressed in cell lines compared to primary tissues, and this alteration in expression is a result of less repressive TF targeting. We confirmed these regulatory changes for four TFs, including SMAD5, using independent ChIP-seq data from ENCODE. Conclusions Our results provide novel insights into the regulatory mechanisms controlling the expression differences between cell lines and tissues. The strong changes in TF regulation that we observe suggest that network changes, in addition to transcriptional levels, should be considered when using cell lines as models for tissues.

Abstract Genetic variants annotated to the hedgehog interacting protein (HHIP) are robustly associated with chronic obstructive pulmonary disease (COPD). Hhip haploinsufficiency in mice leads to increased susceptibility towards the development of emphysema following exposure to chronic cigarette smoke (CS). To explore the molecular pathways which contribute to increased susceptibility, we performed metabolomic profiling using high performance liquid chromatography tandem mass spectroscopy (LC/MS-MS) on plasma, urine, and lung tissue of Hhip +/− heterozygotes and wild type (Hhip +/+) C57/BL6 mice exposed to either room-air or CS for six months. Univariate comparisons between groups were made with a combined fold change ≥2 and Student’s t-test p-value

Abstract Prenatal exposure to mercury, a known neurotoxic metal, is associated with lower cognitive performance during childhood. Disruption of fetal epigenetic programming could explain mercury’s neurodevelopmental effects. We screened for epigenome-wide methylation differences associated with maternal prenatal blood mercury levels in 321 cord blood DNA samples and examined the persistence of these alterations during early (n = 75; 2.9–4.9 years) and mid-childhood (n = 291; 6.7–10.5 years). Among males, prenatal mercury levels were associated with lower regional cord blood DNA methylation at the Paraoxonase 1 gene (PON1) that persisted in early childhood and was attenuated in mid-childhood blood. Cord blood methylation at the PON1 locus predicted lower cognitive test scores measured during early childhood. Methylation at the PON1 locus was associated with PON1 expression in an independent set of cord blood samples. The observed persistent epigenetic disruption of the PON1 gene may modulate mercury toxicity in humans and might serve as a biomarker of exposure and disease susceptibility.

Susceptibility and progression of lung disease, as well as response to treatment, often differ by sex, yet the metabolic mechanisms driving these sex-specific differences are still poorly understood. Women with chronic obstructive pulmonary disease (COPD) have less emphysema and more small airway disease on average than men, though these differences become less pronounced with more severe airflow limitation. While small studies of targeted metabolites have identified compounds differing by sex and COPD status, the sex-specific effect of COPD on systemic metabolism has yet to be interrogated. Significant sex differences were observed in 9 of the 11 modules identified in COPDGene. Sex-specific associations by COPD status and emphysema were observed in 3 modules for each phenotype. Sex stratified individual metabolite associations with COPD demonstrated male-specific associations in sphingomyelins and female-specific associations in acyl carnitines and phosphatidylethanolamines. There was high preservation of module assignments in SPIROMICS (SubPopulations and InteRmediate Outcome Measures In COPD Study) and similar female-specific shift in acyl carnitines. Several COPD associated metabolites differed by sex. Acyl carnitines and sphingomyelins demonstrate sex-specific abundances and may represent important metabolic signatures of sex differences in COPD. Accurately characterizing the sex-specific molecular differences in COPD is vital for personalized diagnostics and therapeutics.

Identifying differentially variable (DV) genomic probes is becoming a new approach to detect novel genomic risk factors for complex human diseases. The F test is the standard equal-variance test in statistics. For high-throughput genomic data, the probe-wise F test has been successfully used to detect biologically relevant DNA methylation marks that have different variances between two groups of subjects (e.g., cases versus controls). In addition to DNA methylation, microRNA (miRNA) is another important mechanism of epigenetics. However, to the best of our knowledge, no studies have identified DV miRNAs. In this article, we proposed a novel model-based clustering method to improve the power of the probe-wise F test to detect DV miRNAs. We imposed special structures on covariance matrices for each cluster of miRNAs based on the prior information about the relationship between variances in cases and controls and about the independence among them. Simulation studies showed that the proposed method seems promising in detecting DV probes. Based on two real datasets about human hepatocellular carcinoma (HCC), we identified 7 DV-only miRNAs (hsa-miR-1826, hsa-miR-191, hsa-miR-194-star, hsa-miR-222, hsa-miR-502-3p, hsa-miR-93, and hsa-miR-99b) using the proposed method, one (hsa-miR-1826) of which has not yet been reported to be related to HCC in the literature.