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.

Skeletal muscle DNA methylation: Effects of exercise and HIV.

Aging, human immunodeficiency virus (HIV) infection, and antiretroviral therapy modify the epigenetic profile and function of cells and tissues, including skeletal muscle (SkM). In some cells, accelerated epigenetic aging begins very soon after the initial HIV infection, potentially setting the stage for the early onset of frailty. Exercise imparts epigenetic modifications in SkM that may underpin some health benefits, including delayed frailty, in people living with HIV (PWH). In this first report of exercise-related changes in SkM DNA methylation among PWH, we investigated the impact of 24 weeks of aerobic and resistance exercise training on SkM (vastus lateralis) DNA methylation profiles and epigenetic age acceleration (EAA) in older, virally suppressed PWH (n = 12) and uninfected controls (n = 18), and associations of EAA with physical function at baseline. We identified 983 differentially methylated positions (DMPs) in PWH and controls at baseline and 237 DMPs after training. The influence of HIV serostatus on SkM methylation was more pronounced than that of exercise training. There was little overlap in the genes associated with the probes most significantly differentiated by exercise training within each group. Baseline EAA (mean ± SD) was similar between PWH (-0.4 ± 2.5 years) and controls (0.2 ± 2.6 years), and the exercise effect was not significant (p = 0.79). EAA and physical function at baseline were not significantly correlated (all p ≥ 0.10). This preliminary investigation suggests HIV-specific epigenetic adaptations in SkM with exercise training but confirmation in a larger study that includes transcriptomic analysis is warranted.

TestDimorph: An R package for analysis of interpopulation sexual dimorphism differences using summary statistics.

OBJECTIVE: The degree of sexual dimorphism in certain traits between males and females differ from one sample to another. Although trait differences by sex are often reported in bioanthropological research, few studies test for statistical significance or make raw data available. TestDimorph is the first R package dedicated to testing and comparing the degree of sexual dimorphism among different samples by leveraging summary statistics. MATERIALS AND METHODS: We provide two approaches of analysis of inter-sample differences in degree of sexual dimorphism: univariate and multivariate for two or more samples. The methods follow upon publications primarily from the AJBA. Within-sex size variability between samples is compared using one-way ANOVA followed by control for multiple pairwise comparisons. In addition, we compute the overlapping area between the density functions of two normal distributions from the mixture intersection index or the non-overlapping area using the dissimilarity index as well as Hedges' g with inferential support using the 95% confidence interval. Finally, we use a multivariate analysis of differences in patterning of sexual dimorphism between samples. RESULTS: We demonstrate various results from applying TestDimorph functions to data supplied with the package. DISCUSSION: The package has many features including functionality for working with summary statistics, simulating data from summary statistics, and the extraction of summary statistics from raw data, so that the entire analysis can be performed through the package.

The Microbiome, Epigenome, and Diet in Adults with Obesity during Behavioral Weight Loss.

Obesity has been linked to the gut microbiome, epigenome, and diet, yet these factors have not been studied together during obesity treatment. Our objective was to evaluate associations among gut microbiota (MB), DNA methylation (DNAme), and diet prior to and during a behavioral weight loss intervention. Adults (n = 47, age 40.9 ± 9.7 years, body mass index (BMI) 33.5 ± 4.5 kg/m2, 77% female) with data collected at baseline (BL) and 3 months (3 m) were included. Fecal MB was assessed via 16S sequencing and whole blood DNAme via the Infinium EPIC array. Food group and nutrient intakes and Healthy Eating Index (HEI) scores were calculated from 7-day diet records. Linear models were used to test for the effect of taxa relative abundance on DNAme and diet cross-sectionally at each time point, adjusting for confounders and a false discovery rate of 5%. Mean weight loss was 6.2 ± 3.9% at 3 m. At BL, one MB taxon, Ruminiclostridium, was associated with DNAme of the genes COL20A1 (r = 0.651, p = 0.029), COL18A1 (r = 0.578, p = 0.044), and NT5E (r = 0.365, p = 0.043). At 3 m, there were 14 unique MB:DNAme associations, such as Akkermansia with DNAme of GUSB (r = -0.585, p = 0.003), CRYL1 (r = -0.419, p = 0.007), C9 (r = -0.439, p = 0.019), and GMDS (r = -0.559, p = 0.046). Among taxa associated with DNAme, no significant relationships were seen with dietary intakes of relevant nutrients, food groups, or HEI scores. Our findings indicate that microbes linked to mucin degradation, short-chain fatty acid production, and body weight are associated with DNAme of phenotypically relevant genes. These relationships offer an initial understanding of the possible routes by which alterations in gut MB may influence metabolism during weight loss.

Large scale proteomic studies create novel privacy considerations.

Privacy protection is a core principle of genomic but not proteomic research. We identified independent single nucleotide polymorphism (SNP) quantitative trait loci (pQTL) from COPDGene and Jackson Heart Study (JHS), calculated continuous protein level genotype probabilities, and then applied a naïve Bayesian approach to link SomaScan 1.3K proteomes to genomes for 2812 independent subjects from COPDGene, JHS, SubPopulations and InteRmediate Outcome Measures In COPD Study (SPIROMICS) and Multi-Ethnic Study of Atherosclerosis (MESA). We correctly linked 90-95% of proteomes to their correct genome and for 95-99% we identify the 1% most likely links. The linking accuracy in subjects with African ancestry was lower (~ 60%) unless training included diverse subjects. With larger profiling (SomaScan 5K) in the Atherosclerosis Risk Communities (ARIC) correct identification was > 99% even in mixed ancestry populations. We also linked proteomes-to-proteomes and used the proteome only to determine features such as sex, ancestry, and first-degree relatives. When serial proteomes are available, the linking algorithm can be used to identify and correct mislabeled samples. This work also demonstrates the importance of including diverse populations in omics research and that large proteomic datasets (> 1000 proteins) can be accurately linked to a specific genome through pQTL knowledge and should not be considered unidentifiable.

Multi-Omic Signatures of Sarcoidosis and Progression in Bronchoalveolar Lavage Cells.

Introduction: Sarcoidosis is a heterogeneous, granulomatous disease that can prove difficult to diagnose, with no accurate biomarkers of disease progression. Therefore, we profiled and integrated the DNA methylome, mRNAs, and microRNAs to identify molecular changes associated with sarcoidosis and disease progression that might illuminate underlying mechanisms of disease and potential genomic biomarkers. Methods: Bronchoalveolar lavage cells from 64 sarcoidosis subjects and 16 healthy controls were used. DNA methylation was profiled on Illumina HumanMethylationEPIC arrays, mRNA by RNA-sequencing, and miRNAs by small RNA-sequencing. Linear models were fit to test for effect of diagnosis and phenotype, adjusting for age, sex, and smoking. We built a supervised multi-omics model using a subset of features from each dataset. Results: We identified 46,812 CpGs, 1,842 mRNAs, and 5 miRNAs associated with sarcoidosis versus controls and 1 mRNA, SEPP1 - a protein that supplies selenium to cells, associated with disease progression. Our integrated model emphasized the prominence of the PI3K/AKT1 pathway in sarcoidosis, which is important in T cell and mTOR function. Novel immune related genes and miRNAs including LYST, RGS14, SLFN12L, and hsa-miR-199b-5p, distinguished sarcoidosis from controls. Our integrated model also demonstrated differential expression/methylation of IL20RB, ABCC11, SFSWAP, AGBL4, miR-146a-3p, and miR-378b between non-progressive and progressive sarcoidosis. Conclusions: Leveraging the DNA methylome, transcriptome, and miRNA-sequencing in sarcoidosis BAL cells, we detected widespread molecular changes associated with disease, many which are involved in immune response. These molecules may serve as diagnostic/prognostic biomarkers and/or drug targets, although future testing will be required for confirmation.

Gene expression associations with body mass index in the Multi-Ethnic Study of Atherosclerosis.

BACKGROUND/OBJECTIVES: Obesity, defined as excessive fat accumulation that represents a health risk, is increasing in adults and children, reaching global epidemic proportions. Body mass index (BMI) correlates with body fat and future health risk, yet differs in prediction by fat distribution, across populations and by age. Nonetheless, few genetic studies of BMI have been conducted in ancestrally diverse populations. Gene expression association with BMI was assessed in the Multi-Ethnic Study of Atherosclerosis (MESA) in four self-identified race and ethnicity (SIRE) groups to identify genes associated with obesity. SUBJECTS/METHODS: RNA-sequencing was performed on 1096 MESA participants (37.8% white, 24.3% Hispanic, 28.4% African American, and 9.5% Chinese American) and linear models were used to assess the association of expression from each gene for its effect on BMI, adjusting for age, sex, sequencing center, study site, five expression and four genetic principal components in each self-identified race group. Sample-size-weighted meta-analysis was performed to identify genes with BMI-associated expression across ancestry groups. RESULTS: Within individual SIRE groups, there were zero to three genes whose expression is significantly (p < 1.97 × 10-6) associated with BMI. Across all groups, 45 genes were identified by meta-analysis whose expression was significantly associated with BMI, explaining 29.7% of BMI variation. The 45 genes are expressed in a variety of tissues and cell types and are enriched for obesity-related processes including erythrocyte function, oxygen binding and transport, and JAK-STAT signaling. CONCLUSIONS: We have identified genes whose expression is significantly associated with obesity in a multi-ethnic cohort. We have identified novel genes associated with BMI as well as confirmed previously identified genes from earlier genetic analyses. These novel genes and their biological pathways represent new targets for understanding the biology of obesity as well as new therapeutic intervention to reduce obesity and improve global public health.

Multiomic Signatures of Chronic Beryllium Disease Bronchoalveolar Lavage Cells Relate to T-Cell Function and Innate Immunity.

Chronic beryllium disease (CBD) is a Th1 granulomatous lung disease preceded by sensitization to beryllium (BeS). We profiled the methylome, transcriptome, and selected proteins in the lung to identify molecular signatures and networks associated with BeS and CBD. BAL cell DNA and RNA were profiled using microarrays from CBD (n = 30), BeS (n = 30), and control subjects (n = 12). BAL fluid proteins were measured using Olink Immune Response Panel proteins from CBD (n = 22) and BeS (n = 22) subjects. Linear models identified features associated with CBD, adjusting for covariation and batch effects. Multiomic integration methods identified correlated features between datasets. We identified 1,546 differentially expressed genes in CBD versus control subjects and 204 in CBD versus BeS. Of the 101 shared transcripts, 24 have significant cis relationships between gene expression and DNA methylation, assessed using expression quantitative trait methylation analysis, including genes not previously identified in CBD. A multiomic model of top DNA methylation and gene expression features demonstrated that the first component separated CBD from other samples and the second component separated control subjects from remaining samples. The top features on component one were enriched for T-lymphocyte function, and the top features on component two were enriched for innate immune signaling. We identified six differentially abundant proteins in CBD versus BeS, with two (SIT1 and SH2D1A) selected as important RNA features in the multiomic model. Our integrated analysis of DNA methylation, gene expression, and proteins in the lung identified multiomic signatures of CBD that differentiated it from BeS and control subjects.

Colocalization of Gene Expression and DNA Methylation with Genetic Risk Variants Supports Functional Roles of MUC5B and DSP in Idiopathic Pulmonary Fibrosis.

Rationale: Common genetic variants have been associated with idiopathic pulmonary fibrosis (IPF). Objectives: To determine functional relevance of the 10 IPF-associated common genetic variants we previously identified. Methods: We performed expression quantitative trait loci (eQTL) and methylation quantitative trait loci (mQTL) mapping, followed by co-localization of eQTL and mQTL with genetic association signals and functional validation by luciferase reporter assays. Illumina multi-ethnic genotyping arrays, mRNA sequencing, and Illumina 850k methylation arrays were performed on lung tissue of participants with IPF (234 RNA and 345 DNA samples) and non-diseased controls (188 RNA and 202 DNA samples). Measurements and Main Results: Focusing on genetic variants within 10 IPF-associated genetic loci, we identified 27 eQTLs in controls and 24 eQTLs in cases (false-discovery-rate-adjusted P < 0.05). Among these signals, we identified associations of lead variants rs35705950 with expression of MUC5B and rs2076295 with expression of DSP in both cases and controls. mQTL analysis identified CpGs in gene bodies of MUC5B (cg17589883) and DSP (cg08964675) associated with the lead variants in these two loci. We also demonstrated strong co-localization of eQTL/mQTL and genetic signal in MUC5B (rs35705950) and DSP (rs2076295). Functional validation of the mQTL in MUC5B using luciferase reporter assays demonstrates that the CpG resides within a putative internal repressor element. Conclusions: We have established a relationship of the common IPF genetic risk variants rs35705950 and rs2076295 with respective changes in MUC5B and DSP expression and methylation. These results provide additional evidence that both MUC5B and DSP are involved in the etiology of IPF.

Epigenetics and sarcoidosis.

Epigenetic modifications are emerging as important regulatory mechanisms of gene expression in lung disease, given that they are influenced by environmental exposures and genetic variants, and that they regulate immune and fibrotic processes. In this review, we introduce these concepts with a focus on the study of DNA methylation and histone modifications and discuss how they have been applied to lung disease, and how they can be applied to sarcoidosis. This information has implications for other exposure and immunologically mediated lung diseases, such as chronic beryllium disease, hypersensitivity pneumonitis, and asbestosis.

Molecular Signatures of Idiopathic Pulmonary Fibrosis.

Molecular patterns and pathways in idiopathic pulmonary fibrosis (IPF) have been extensively investigated, but few studies have assimilated multiomic platforms to provide an integrative understanding of molecular patterns that are relevant in IPF. Herein, we combine the coding and noncoding transcriptomes, DNA methylomes, and proteomes from IPF and healthy lung tissue to identify molecules and pathways associated with this disease. RNA sequencing, Illumina MethylationEPIC array, and liquid chromatography-mass spectrometry proteomic data were collected on lung tissue from 24 subjects with IPF and 14 control subjects. Significant differential features were identified by using linear models adjusting for age and sex, inflation, and bias when appropriate. Data Integration Analysis for Biomarker Discovery Using a Latent Component Method for Omics Studies was used for integrative multiomic analysis. We identified 4,643 differentially expressed transcripts aligning to 3,439 genes, 998 differentially abundant proteins, 2,500 differentially methylated regions, and 1,269 differentially expressed long noncoding RNAs (lncRNAs) that were significant after correcting for multiple tests (false discovery rate < 0.05). Unsupervised hierarchical clustering using 20 coding mRNA, protein, methylation, and lncRNA features with the highest loadings on the top latent variable from the four data sets demonstrates perfect separation of IPF and control lungs. Our analysis confirmed previously validated molecules and pathways known to be dysregulated in disease and implicated novel molecular features as potential drivers and modifiers of disease. For example, 4 proteins, 18 differentially methylated regions, and 10 lncRNAs were found to have strong correlations (|r| > 0.8) with MMP7 (matrix metalloproteinase 7). Therefore, by using a system biology approach, we have identified novel molecular relationships in IPF.

Multiomic Predictors of Short-Term Weight Loss and Clinical Outcomes During a Behavioral-Based Weight Loss Intervention.

OBJECTIVE: Identifying predictors of weight loss and clinical outcomes may increase understanding of individual variability in weight loss response. We hypothesized that baseline multiomic features, including DNA methylation (DNAme), metabolomics, and gut microbiome, would be predictive of short-term changes in body weight and other clinical outcomes within a comprehensive weight loss intervention. METHODS: Healthy adults with overweight or obesity (n = 62, age 18-55 years, BMI 27-45 kg/m2 , 75.8% female) participated in a 1-year behavioral weight loss intervention. To identify baseline omic predictors of changes in clinical outcomes at 3 and 6 months, whole-blood DNAme, plasma metabolites, and gut microbial genera were analyzed. RESULTS: A network of multiomic relationships informed predictive models for 10 clinical outcomes (body weight, waist circumference, fat mass, hemoglobin A1c , homeostatic model assessment of insulin resistance, total cholesterol, triglycerides, C-reactive protein, leptin, and ghrelin) that changed significantly (P < 0.05). For eight of these, adjusted R2 ranged from 0.34 to 0.78. Our models identified specific DNAme sites, gut microbes, and metabolites that were predictive of variability in weight loss, waist circumference, and circulating triglycerides and that are biologically relevant to obesity and metabolic pathways. CONCLUSIONS: These data support the feasibility of using baseline multiomic features to provide insight for precision nutrition-based weight loss interventions.

Host methylation predicts SARS-CoV-2 infection and clinical outcome.

BACKGROUND: Since the onset of the SARS-CoV-2 pandemic, most clinical testing has focused on RT-PCR1. Host epigenome manipulation post coronavirus infection2-4 suggests that DNA methylation signatures may differentiate patients with SARS-CoV-2 infection from uninfected individuals, and help predict COVID-19 disease severity, even at initial presentation. METHODS: We customized Illumina's Infinium MethylationEPIC array to enhance immune response detection and profiled peripheral blood samples from 164 COVID-19 patients with longitudinal measurements of disease severity and 296 patient controls. RESULTS: Epigenome-wide association analysis revealed 13,033 genome-wide significant methylation sites for case-vs-control status. Genes and pathways involved in interferon signaling and viral response were significantly enriched among differentially methylated sites. We observe highly significant associations at genes previously reported in genetic association studies (e.g. IRF7, OAS1). Using machine learning techniques, models built using sparse regression yielded highly predictive findings: cross-validated best fit AUC was 93.6% for case-vs-control status, and 79.1%, 80.8%, and 84.4% for hospitalization, ICU admission, and progression to death, respectively. CONCLUSIONS: In summary, the strong COVID-19-specific epigenetic signature in peripheral blood driven by key immune-related pathways related to infection status, disease severity, and clinical deterioration provides insights useful for diagnosis and prognosis of patients with viral infections.

Host methylation predicts SARS-CoV-2 infection and clinical outcome.

BACKGROUND: Since the onset of the SARS-CoV-2 pandemic, most clinical testing has focused on RT-PCR1. Host epigenome manipulation post coronavirus infection2-4 suggests that DNA methylation signatures may differentiate patients with SARS-CoV-2 infection from uninfected individuals, and help predict COVID-19 disease severity, even at initial presentation. METHODS: We customized Illumina's Infinium MethylationEPIC array to enhance immune response detection and profiled peripheral blood samples from 164 COVID-19 patients with longitudinal measurements of disease severity and 296 patient controls. RESULTS: Epigenome-wide association analysis revealed 13,033 genome-wide significant methylation sites for case-vs-control status. Genes and pathways involved in interferon signaling and viral response were significantly enriched among differentially methylated sites. We observe highly significant associations at genes previously reported in genetic association studies (e.g. IRF7, OAS1). Using machine learning techniques, models built using sparse regression yielded highly predictive findings: cross-validated best fit AUC was 93.6% for case-vs-control status, and 79.1%, 80.8%, and 84.4% for hospitalization, ICU admission, and progression to death, respectively. CONCLUSIONS: In summary, the strong COVID-19-specific epigenetic signature in peripheral blood driven by key immune-related pathways related to infection status, disease severity, and clinical deterioration provides insights useful for diagnosis and prognosis of patients with viral infections.

Viral infections affect the body in many ways, including via changes to the epigenome, the sum of chemical modifications to an individual's collection of genes that affect gene activity. Here, we analyzed the epigenome in blood samples from people with and without COVID-19 to determine whether we could find changes consistent with SARS-CoV-2 infection. Using a combination of statistical and machine learning techniques, we identify markers of SARS-CoV-2 infection as well as of severity and progression of COVID-19 disease. These signals of disease progression were present from the initial blood draw when first walking into the hospital. Together, these approaches demonstrate the potential of measuring the epigenome for monitoring SARS-CoV-2 status and severity.

Chronic Hypersensitivity Pneumonitis, an Interstitial Lung Disease with Distinct Molecular Signatures.

Rationale: Chronic hypersensitivity pneumonitis (CHP) is caused by an immune response to antigen inhalation and is characterized by variable histopathological and clinical features. A subset of subjects with CHP have usual interstitial pneumonia and appear to be clinically similar to subjects with idiopathic pulmonary fibrosis (IPF).Objectives: To determine the common and unique molecular features of CHP and IPF.Methods: Transcriptome analysis of lung samples from CHP (n = 82), IPF (n = 103), and unaffected controls (n = 103) was conducted. Differential gene expression was determined adjusting for sex, race, age, and smoking history and using false discovery rate to control for multiple comparisons.Measurements and Main Results: When compared with controls, we identified 413 upregulated and 317 downregulated genes in CHP and 861 upregulated and 322 downregulated genes in IPF. Concordantly upregulated or downregulated genes in CHP and IPF were related to collagen catabolic processes and epithelial development, whereas genes specific to CHP (differentially expressed in CHP when compared with control and not differentially expressed in IPF) were related to chemokine-mediated signaling and immune responsiveness. Using weighted gene coexpression network analysis, we found that among subjects with CHP, genes involved in adaptive immunity or epithelial cell development were associated with improved or reduced lung function, respectively, and that MUC5B expression was associated with epithelial cell development. MUC5B expression was also associated with lung fibrosis and honeycombing.Conclusions: Gene expression analysis of CHP and IPF identified signatures common to CHP and IPF, as well as genes uniquely expressed in CHP. Select modules of gene expression are characterized by distinct clinical and pathological features of CHP.