RNA sequencing identifies novel non-coding RNA and exon-specific effects associated with cigarette smoking

Overlap between COPD genetic association results and transcriptional quantitative trait loci

Rationale

Genome-wide association studies (GWAS) have identified multiple genetic loci associated with chronic obstructive pulmonary disease (COPD). When integrated with GWAS results, expression quantitative trait locus (eQTL) studies can provide insight into biological mechanisms involved in disease by identifying single nucleotide polymorphisms (SNPs) that contribute to whole gene expression. However, there are multiple genetically driven regulatory and isoform-specific effects which cannot be detected in traditional eQTL analyses. Here, we identify SNPs that are associated with alternative splicing (sQTL) in addition to eQTLs to identify novel functions for COPD associated genetic variants.

Methods

We performed RNA sequencing on whole blood from 3743 subjects in the COPDGene Study. RNA sequencing data from lung tissue of 1241 subjects from the Lung Tissue Research Consortium (LTRC), and whole genome sequencing data on all subjects. Associations between all SNPs within 1000 kb of a gene (cis-) and splice and gene expression quantifications were tested using tensorQTL. In COPDGene a total of 11,869,333 SNPs were tested for association with 58,318 splice clusters, and 8,792,206 SNPs were tested for association with 70,094 splice clusters in LTRC. We assessed colocalization with COPD-associated SNPs from a published GWAS[1].

Results

After adjustment for multiple statistical testing, we identified 28,110 splice-sites corresponding to 3,889 unique genes that were significantly associated with genotype in COPDGene whole blood, and 58,258 splice-sites corresponding to 10,307 unique genes associated with genotype in LTRC lung tissue. We found 7,576 sQTL splice-sites corresponding to 2,110 sQTL genes were shared between whole blood and lung, while 20,534 sQTL splice-sites in 3,518 genes were unique to blood and 50,682 splice-sites in 9,677 genes were unique to lung. To determine what proportion of COPD-associated SNPs were associated with transcriptional splicing, we performed colocalization analysis between COPD GWAS and sQTL data, and found that 38 genomic windows, corresponding to 38 COPD GWAS loci had evidence of colocalization between QTLs and COPD. The top five colocalizations between COPD and lung sQTLs include NPNT , FBXO38 , HHIP , NTN4 and BTC .

Conclusions

A total of 38 COPD GWAS loci contain evidence of sQTLs, suggesting that analysis of sQTLs in whole blood and lung tissue can provide novel insights into disease mechanisms.

Multi-center Korean cohort study based on RNA-sequencing data targeting COPD patients

In 2023, WHO ranked chronic obstructive pulmonary disease (COPD) as the third leading cause of death, with 3.23 million fatalities in 2019. The intricate nature of the disease, which is influenced by genetics, environment, and lifestyle, is evident. The effect of air pollution and changes in atmospheric substances because of global warming highlight the need for this research. These environmental shifts are associated with the emergence of various respiratory infections such as COVID-19. RNA sequencing is pivotal in airway diseases, including COPD, as it enables comprehensive transcriptome analysis, biomarker discovery, and uncovers novel pathways. It facilitates personalized medicine by tracking dynamic changes in gene expression in response to various triggers. However, the limited research on East Asian populations may overlook the unique nuances of COPD development and progression. Bridging this gap and using peripheral blood samples for systemic analysis are crucial for comprehensive and globally applicable COPD diagnosis and treatment.

Deep learning on graphs for multi-omics classification of COPD

Network approaches have successfully been used to help reveal complex mechanisms of diseases including Chronic Obstructive Pulmonary Disease (COPD). However despite recent advances, we remain limited in our ability to incorporate protein-protein interaction (PPI) network information with omics data for disease prediction. New deep learning methods including convolution Graph Neural Network (ConvGNN) has shown great potential for disease classification using transcriptomics data and known PPI networks from existing databases. In this study, we first reconstructed the COPD-associated PPI network through the AhGlasso (Augmented High-Dimensional Graphical Lasso Method) algorithm based on one independent transcriptomics dataset including COPD cases and controls. Then we extended the existing ConvGNN methods to successfully integrate COPD-associated PPI, proteomics, and transcriptomics data and developed a prediction model for COPD classification. This approach improves accuracy over several conventional classification methods and neural networks that do not incorporate network information. We also demonstrated that the updated COPD-associated network developed using AhGlasso further improves prediction accuracy. Although deep neural networks often achieve superior statistical power in classification compared to other methods, it can be very difficult to explain how the model, especially graph neural network(s), makes decisions on the given features and identifies the features that contribute the most to prediction generally and individually. To better explain how the spectral-based Graph Neural Network model(s) works, we applied one unified explainable machine learning method, SHapley Additive exPlanations (SHAP), and identified CXCL11, IL-2, CD48, KIR3DL2, TLR2, BMP10 and several other relevant COPD genes in subnetworks of the ConvGNN model for COPD prediction. Finally, Gene Ontology (GO) enrichment analysis identified glycosaminoglycan, heparin signaling, and carbohydrate derivative signaling pathways significantly enriched in the top important gene/proteins for COPD classifications.

Identifying chronic obstructive pulmonary disease from integrative omics and clustering in lung tissue

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.

Distinct COPD subtypes in former smokers revealed by gene network perturbation analysis

BACKGROUND

Chronic obstructive pulmonary disease (COPD) varies significantly in symptomatic and physiologic presentation. Identifying disease subtypes from molecular data, collected from easily accessible blood samples, can help stratify patients and guide disease management and treatment.

METHODS

Blood gene expression measured by RNA-sequencing in the COPDGene Study was analyzed using a network perturbation analysis method. Each COPD sample was compared against a learned reference gene network to determine the part that is deregulated. Gene deregulation values were used to cluster the disease samples.

RESULTS

The discovery set included 617 former smokers from COPDGene. Four distinct gene network subtypes are identified with significant differences in symptoms, exercise capacity and mortality. These clusters do not necessarily correspond with the levels of lung function impairment and are independently validated in two external cohorts: 769 former smokers from COPDGene and 431 former smokers in the Multi-Ethnic Study of Atherosclerosis (MESA). Additionally, we identify several genes that are significantly deregulated across these subtypes, including DSP and GSTM1, which have been previously associated with COPD through genome-wide association study (GWAS).

CONCLUSIONS

The identified subtypes differ in mortality and in their clinical and functional characteristics, underlining the need for multi-dimensional assessment potentially supplemented by selected markers of gene expression. The subtypes were consistent across cohorts and could be used for new patient stratification and disease prognosis.

Hepatitis C and HIV detection by blood RNA-sequencing in cohort of smokers

Detection of viruses by RNA and DNA sequencing has improved the understanding of the human virome. We sought to identify blood viral signatures through secondary use of RNA-sequencing (RNA-seq) data in a large study cohort. The ability to reveal undiagnosed infections with public health implications among study subjects with available sequencing data could enable epidemiologic surveys and may lead to diagnosis and therapeutic interventions, leveraging existing research data in a clinical context. We detected viral RNA in peripheral blood RNA-seq data from a COPD-enriched population of current and former smokers. Correlation between viral detection and both reported infections and relevant disease outcomes was evaluated. We identified Hepatitis C virus RNA in 228 subjects and HIV RNA in 30 subjects. Overall, we observed 31 viral species, including Epstein-Barr virus and Cytomegalovirus. We observed an enrichment of Hepatitis C and HIV infections among subjects reporting liver disease and HIV infections, respectively. Higher interferon expression scores were observed in the subjects with Hepatitis C and HIV infections. Through secondary use of RNA-seq from a cohort of current and former smokers, we detected peripheral blood viral signatures. We identified HIV and Hepatitis C virus (HCV), highlighting potential public health implications for the approach described this study. We observed correlations with reported infections, chronic infection outcomes and the host transcriptomic response, providing evidence to support the validity of the approach.

Mechanism of Action of Decitabine in the Treatment of Acute Myeloid Leukemia by Regulating LINC00599

Objective

Acute myeloid leukemia (AML) is a heterogeneous malignancy with a low long-term survival rate. The aim of this study was to investigate the effects of decitabine (DAC) treatment cell proliferation and apoptosis in AML and role of the expression of LINC00599 and, consequently, miR-135a-5p.

Materials and Methods

Human promyelocytic leukemia cells (HL-60) and human acute lymphatic leukemia (CCRF-CEM) cells were treated with various concentrations of DAC. Cell proliferation in each group was detected using the cell counting kit 8. For each group, apoptosis and reactive oxygen species (ROS) levels were detected using flow cytometry. Reverse transcription polymerase chain reaction (RT-PCR) was performed to examine the expression of lncRNA LINC00599. The expression of apoptosis-related proteins was analyzed using western blotting. The regulatory relationship between miR-135a-5p and LINC00599 was verified by constructing miR-135a-5p mimics, miR-135a-5p inhibit, wild type LINC00599 3'-untranslated region (UTR), and mutant LINC00599 3'-UTR. Ki-67 expression in the tumor tissues of nude mice was detected using immunofluorescent assays.

Results

Both DAC and LINC00599 Inhibit groups were able to significantly reduce the proliferation of HL60 and CCRF-CEM cells, increase apoptosis, upregulate the expression of Bad, cleaved caspase-3, and miR-135a-5p, downregulate the expression of Bcl-2, and elevate ROS levels in cells, with these effects being more pronounced after combined treatment with DAC and LINC00599 Inhibit. In comparison to mimic NC, the miR-135a-5p mimic group significantly decreased the relative fluorescence activity ratio of LINC00599 3'-UTR wild-type CCRF-CEM cells. The LINC00599 Inhibit and miR-135a-5p mimic groups exhibited substantially reduced proliferation of HL60 and CCRF-CEM cells, increased apoptosis, upregulated Bad, cleaved caspase-3, and miR-135a-5p expression, along with downregulated Bcl-2 and LINC00599 expression and increased ROS levels in cells; these effects were more pronounced after LINC00599 Inhibit was combined with miR-135a-5p mimics. In vivo experiments revealed that both DAC and LINC00599 Inhibit were able to considerably reduce the long diameter, short meridian, volume, and mass of tumors, increase miR-135a-5p expression, and decrease LINC00599 and ki-67 expression in tumor tissues of nude mice. This effect was more pronounced when the DAC and LINC00599 Inhibit were used in combination.

Conclusion

DAC regulates the expression of miR-135a-5p by regulating the expression of LINC00599, which in turn affects cell proliferation, apoptosis, and tumor proliferation. Our findings provide a theoretical basis for improving the clinical outcome of AML.

Multi-omics insights into the biological mechanisms underlying statistical gene-by-lifestyle interactions with smoking and alcohol consumption

Though both genetic and lifestyle factors are known to influence cardiometabolic outcomes, less attention has been given to whether lifestyle exposures can alter the association between a genetic variant and these outcomes. The Cohorts for Heart and Aging Research in Genomic Epidemiology (CHARGE) Consortium's Gene-Lifestyle Interactions Working Group has recently published investigations of genome-wide gene-environment interactions in large multi-ancestry meta-analyses with a focus on cigarette smoking and alcohol consumption as lifestyle factors and blood pressure and serum lipids as outcomes. Further description of the biological mechanisms underlying these statistical interactions would represent a significant advance in our understanding of gene-environment interactions, yet accessing and harmonizing individual-level genetic and 'omics data is challenging. Here, we demonstrate the coordinated use of summary-level data for gene-lifestyle interaction associations on up to 600,000 individuals, differential methylation data, and gene expression data for the characterization and prioritization of loci for future follow-up analyses. Using this approach, we identify 48 genes for which there are multiple sources of functional support for the identified gene-lifestyle interaction. We also identified five genes for which differential expression was observed by the same lifestyle factor for which a gene-lifestyle interaction was found. For instance, in gene-lifestyle interaction analysis, the T allele of rs6490056 (ALDH2) was associated with higher systolic blood pressure, and a larger effect was observed in smokers compared to non-smokers. In gene expression studies, this allele is associated with decreased expression of ALDH2, which is part of a major oxidative pathway. Other results show increased expression of ALDH2 among smokers. Oxidative stress is known to contribute to worsening blood pressure. Together these data support the hypothesis that rs6490056 reduces expression of ALDH2, which raises oxidative stress, leading to an increase in blood pressure, with a stronger effect among smokers, in whom the burden of oxidative stress is greater. Other genes for which the aggregation of data types suggest a potential mechanism include: GCNT4×current smoking (HDL), PTPRZ1×ever-smoking (HDL), SYN2×current smoking (pulse pressure), and TMEM116×ever-smoking (mean arterial pressure). This work demonstrates the utility of careful curation of summary-level data from a variety of sources to prioritize gene-lifestyle interaction loci for follow-up analyses.

sJIVE: Supervised Joint and Individual Variation Explained

Analyzing multi-source data, which are multiple views of data on the same subjects, has become increasingly common in molecular biomedical research. Recent methods have sought to uncover underlying structure and relationships within and/or between the data sources, and other methods have sought to build a predictive model for an outcome using all sources. However, existing methods that do both are presently limited because they either (1) only consider data structure shared by all datasets while ignoring structures unique to each source, or (2) they extract underlying structures first without consideration to the outcome. The proposed method, supervised joint and individual variation explained (sJIVE), can simultaneously (1) identify shared (joint) and source-specific (individual) underlying structure and (2) build a linear prediction model for an outcome using these structures. These two components are weighted to compromise between explaining variation in the multi-source data and in the outcome. Simulations show sJIVE to outperform existing methods when large amounts of noise are present in the multi-source data. An application to data from the COPDGene study explores gene expression and proteomic patterns associated with lung function.

Impact of the Exposome on the Epigenome in Inflammatory Bowel Disease Patients and Animal Models

Inflammatory bowel diseases (IBD) are chronic inflammatory disorders of the gastrointestinal tract that encompass two main phenotypes, namely Crohn's disease and ulcerative colitis. These conditions occur in genetically predisposed individuals in response to environmental factors. Epigenetics, acting by DNA methylation, post-translational histones modifications or by non-coding RNAs, could explain how the exposome (or all environmental influences over the life course, from conception to death) could influence the gene expression to contribute to intestinal inflammation. We performed a scoping search using Medline to identify all the elements of the exposome that may play a role in intestinal inflammation through epigenetic modifications, as well as the underlying mechanisms. The environmental factors epigenetically influencing the occurrence of intestinal inflammation are the maternal lifestyle (mainly diet, the occurrence of infection during pregnancy and smoking); breastfeeding; microbiota; diet (including a low-fiber diet, high-fat diet and deficiency in micronutrients); smoking habits, vitamin D and drugs (e.g., IBD treatments, antibiotics and probiotics). Influenced by both microbiota and diet, short-chain fatty acids are gut microbiota-derived metabolites resulting from the anaerobic fermentation of non-digestible dietary fibers, playing an epigenetically mediated role in the integrity of the epithelial barrier and in the defense against invading microorganisms. Although the impact of some environmental factors has been identified, the exposome-induced epimutations in IBD remain a largely underexplored field. How these environmental exposures induce epigenetic modifications (in terms of duration, frequency and the timing at which they occur) and how other environmental factors associated with IBD modulate epigenetics deserve to be further investigated.

Immune-Related LncRNAs as Prognostic Factors for Pediatric Rhabdoid Tumor of the Kidney

Background

Immune-related long noncoding RNAs (IrlncRNAs) are recognized as important prognostic factors in a variety of cancers, but thus far, their prognostic value in pediatric rhabdoid tumor of the kidney (pRTK) has not been reported. Here, we clarified the associations between IrlncRNAs and overall survival (OS) of pRTK patients and constructed a model to predict their prognosis.

Methods

We accessed RNA sequencing data and corresponding clinical data of pRTK from the Therapeutically Applicable Research to Generate Effective Treatments (TARGET) database. An expression profile of immune-related genes (Irgenes) and lncRNAs of pRTK was extracted from the RNA sequencing data. IrlncRNAs were defined by co-expression analysis of lncRNAs and Irgenes. The limma R package was used to identify differential expression IrlncRNAs. Univariate and multivariate Cox regression analyses were conducted to build a prognostic IrlncRNAs model. The performance of this prognostic model was validated by multimethods, like ROC curve analysis.

Results

A total of 1097 IrlncRNAs were defined. Univariate Cox regression analysis identified 7 IrlncRNAs (AC004791.2, AP003068.23, RP11-54O7.14, RP11-680F8.1, TBC1D3P1-DHX40P1, TUNAR, and XXbac-BPG308K3.5) and were significantly associated with OS. Multivariate regression analysis constructed the best prognostic model based on the expression of AC004791.2, AP003068.23, RP11-54O7.14, TBC1D3P1-DHX40P1, and TUNAR. According to the prognostic model, a risk score of each patient was calculated, and patients were divided into high-risk and low-risk groups accordingly. The survival time of low-risk patients was significantly better than high-risk patients (p < 0.001). Univariate (hazard ratio 1.098, 95% confidence interval 1.048-1.149, p value <0.001) and multivariate (hazard ratio 1.095, 95% confidence interval 1.043-1.150, p value <0.001) analyses confirmed that the prognostic model was reliable and independent in prediction of OS. Time-dependent ROC analysis showed that 1-year survival AUC of prognostic model, stage, age, and sex was 0.824, 0.673, 0.531, and 0.495, respectively, which suggested that the prognostic model was the best predictor of survival in pRTK patients.

Conclusions

The prognostic model based on 5 IrlncRNAs was robust and could better predict the survival of pRTK than other clinical factors. Additionally, the mechanism of regulation and action of prognosis-associated lncRNAs could provide new avenues for basic research to explore the mechanism of tumor initiation and development in order to prevent and treat pRTK.

Ablation of long noncoding RNA MALAT1 activates antioxidant pathway and alleviates sepsis in mice

The metastasis-associated lung adenocarcinoma transcript1 (MALAT1) is a long noncoding RNA (lncRNA) and is known for its role in cancer development and prognosis. In this study, we report that MALAT1 plays an important role in regulating acute inflammatory responses in sepsis. In patient samples, MALAT1 expression was positively correlated with severity of sepsis. In cultured macrophages, LPS treatment significantly induced MALAT1 expression, while genetic ablation of MALAT1 greatly reduced proinflammatory cytokine levels. Furthermore, MALAT1-ablated mice had significantly increased survival rates in cecal ligation and puncture (CLP)-induced sepsis and LPS-induced endotoxemia. One novel and salient feature of MALAT1-ablated mice is greatly reduced ROS level in macrophages and other cell types and increased glutathione/oxidized glutathione (GSH/GSSG) ratio in macrophages, suggesting an increased antioxidant capacity. We showed a mechanism for MALAT1 ablation leading to enhanced antioxidant capacity is through activation of methionine cycle by epitranscriptomical regulation of methionine adenosyltransferase 2A (MAT2A). MAT2A 3′UTR can be methylated by METTL16 which was known to directly bind to MALAT1. MALAT1 ablation was found to reduce methylation in MAT2A hairpin1 and increase MAT2A protein levels. Our results suggest a MALAT1-METTL16-MAT2A interactive axis which may be targeted for treatments of sepsis.

Blood gene expression risk profiles and interstitial lung abnormalities: COPDGene and ECLIPSE cohort studies

BACKGROUND

Interstitial lung abnormalities (ILA) are radiologic findings that may progress to idiopathic pulmonary fibrosis (IPF). Blood gene expression profiles can predict IPF mortality, but whether these same genes associate with ILA and ILA outcomes is unknown. This study evaluated if a previously described blood gene expression profile associated with IPF mortality is associated with ILA and all-cause mortality.

METHODS

In COPDGene and ECLIPSE study participants with visual scoring of ILA and gene expression data, we evaluated the association of a previously described IPF mortality score with ILA and mortality. We also trained a new ILA score, derived using genes from the IPF score, in a subset of COPDGene. We tested the association with ILA and mortality on the remainder of COPDGene and ECLIPSE.

RESULTS

In 1469 COPDGene (training n = 734; testing n = 735) and 571 ECLIPSE participants, the IPF score was not associated with ILA or mortality. However, an ILA score derived from IPF score genes was associated with ILA (meta-analysis of test datasets OR 1.4 [95% CI: 1.2-1.6]) and mortality (HR 1.25 [95% CI: 1.12-1.41]). Six of the 11 genes in the ILA score had discordant directions of effects compared to the IPF score. The ILA score partially mediated the effects of age on mortality (11.8% proportion mediated).

CONCLUSIONS

An ILA gene expression score, derived from IPF mortality-associated genes, identified genes with concordant and discordant effects on IPF mortality and ILA. These results suggest shared, and unique biologic processes, amongst those with ILA, IPF, aging, and death.

Learning processes in hierarchical pairs regulate entire gene expression in cells

Expression of numerous genes is precisely controlled in a cell in various contexts. While genetic and epigenetic mechanisms contribute to this regulation, how each mechanism cooperates to ensure the proper expression patterns of the whole gene remains unclear. Here, I theoretically show that the repetition of simple biological processes makes cells functional with the appropriate expression patterns of all genes if the inappropriateness of current expression ratios is roughly fed back to the epigenetic states. A learning pair model is developed, in which two factors autonomously approach the target ratio by repeating two stochastic processes; competitive amplification with a small addition term and decay depending on the difference between the current and target ratios. Furthermore, thousands of factors are self-regulated in a hierarchical-pair architecture, in which the activation degrees competitively amplify, while transducing the activation signal, and decay at four different probabilities. Changes in whole-gene expression during human early embryogenesis and hematopoiesis are reproduced in simulation using this epigenetic learning process in a single genetically-determined hierarchical-pair architecture of gene regulatory cascades. On the background of this learning process, I propose the law of biological inertia, which means that a living cell basically maintains the expression pattern while renewing its contents.

An interferon-inducible signature of airway disease from blood gene expression profiling

BACKGROUND

The molecular basis of airway remodelling in chronic obstructive pulmonary disease (COPD) remains poorly understood. We identified gene expression signatures associated with chest computed tomography (CT) scan airway measures to understand molecular pathways associated with airway disease.

METHODS

In 2396 subjects in the COPDGene Study, we examined the relationship between quantitative CT airway phenotypes and blood transcriptomes to identify airway disease-specific genes and to define an airway wall thickness (AWT) gene set score. Multivariable regression analyses were performed to identify associations of the AWT score with clinical phenotypes, bronchial gene expression and genetic variants.

RESULTS

Type 1 interferon (IFN)-induced genes were consistently associated with AWT, square root wall area of a hypothetical airway with 10 mm internal perimeter (Pi10) and wall area percentage, with the strongest enrichment in AWT. A score derived from 18 genes whose expression was associated with AWT was associated with COPD-related phenotypes including reduced lung function (forced expiratory volume in 1 s percentage predicted β= -3.4; p<0.05) and increased exacerbations (incidence rate ratio 1.7; p<0.05). The AWT score was reproducibly associated with AWT in bronchial samples from 23 subjects (β=3.22; p<0.05). The blood AWT score was associated with genetic variant rs876039, an expression quantitative trait locus for IKZF1, a gene that regulates IFN signalling and is associated with inflammatory diseases.

CONCLUSIONS

A gene expression signature with IFN-stimulated genes from peripheral blood and bronchial brushings is associated with CT AWT, lung function and exacerbations. Shared genes and genetic associations suggest viral responses and/or autoimmune dysregulation as potential underlying mechanisms of airway disease in COPD.

Associations of Monocyte Count and Other Immune Cell Types with Interstitial Lung Abnormalities

Rationale: Higher blood monocyte counts are associated with worse survival in adults with clinically diagnosed pulmonary fibrosis. Their association with the development and progression of interstitial lung abnormalities (ILA) in humans is unknown. Objectives: We evaluated the associations of blood monocyte count, and other immune cell types, with ILA, high-attenuation areas, and FVC in four independent cohorts. Methods: We included participants with measured monocyte counts and computed tomographic (CT) imaging enrolled in MESA (Multi-Ethnic Study of Atherosclerosis, n = 484), AGES-Reykjavik (Age/Gene Environment Susceptibility Study, n = 3,547), COPDGene (Genetic Epidemiology of COPD, n = 2,719), and the ECLIPSE (Evaluation of COPD Longitudinally to Identify Predictive Surrogate End-points, n = 646). Measurements and Main Results: After adjustment for covariates, a 1-SD increment in blood monocyte count was associated with ILA in MESA (odds ratio [OR], 1.3; 95% confidence interval [CI], 1.0-1.8), AGES-Reykjavik (OR, 1.2; 95% CI, 1.1-1.3), COPDGene (OR, 1.3; 95% CI, 1.2-1.4), and ECLIPSE (OR, 1.2; 95% CI, 1.0-1.4). A higher monocyte count was associated with ILA progression over 5 years in AGES-Reykjavik (OR, 1.2; 95% CI, 1.0-1.3). Compared with participants without ILA, there was a higher percentage of activated monocytes among those with ILA in MESA. Higher monocyte count was associated with greater high-attenuation areas in MESA and lower FVC in MESA and COPDGene. Associations of other immune cell types were less consistent. Conclusions: Higher blood monocyte counts were associated with the presence and progression of interstitial lung abnormalities and lower FVC.

Blood RNA sequencing shows overlapping gene expression across COPD phenotype domains

RATIONALE

COPD can be assessed using multidimensional grading systems with components from three domains: pulmonary function tests, symptoms and systemic features. Clinically, measures may be used interchangeably, though it is not known if they share similar pathobiology.

OBJECTIVE

To use RNA sequencing (RNA-seq) to determine if there is an overlap in the underlying biological mechanisms and consequences driving different components of the multidimensional grading systems.

METHODS

Whole blood was collected for RNA-seq from current and former smokers in the Genetic Epidemiology of COPD study. We tested the overlap in gene expression and biological pathways associated with case-control status and quantitative COPD phenotypes within and between the three domains.

RESULTS

In 2647 subjects, there were 3030 genes differentially expressed in any of the three domains or case-control status. There were five genes that overlapped between the three domains and case-control status, including G protein-coupled receptor 15(GPR15), sestrin 1 (SESN1) and interferon-induced guanylate-binding protein 1 (GBP1), which were associated with longitudinal decline in FEV1. The overlap between the three domains was enriched for pathways related to cellular components.

CONCLUSIONS

We identified gene sets and pathways that overlap between 12 COPD-related phenotypes and case-control status. There were no pathways represented in the overlap between the three domains and case-control status, but we identified multiple genes that demonstrated a consistent pattern of expression across several of the phenotypes. Patterns of gene expression correlation were generally similar to the correlation of clinical phenotypes in the PFT and symptom domains but not the systemic features.

Development of a Blood-based Transcriptional Risk Score for Chronic Obstructive Pulmonary Disease

Rationale: The ability of peripheral blood biomarkers to assess chronic obstructive pulmonary disease (COPD) risk and progression is unknown. Genetics and gene expression may capture important aspects of COPD-related biology that predict disease activity. Objectives: Develop a transcriptional risk score (TRS) for COPD and assess the contribution of the TRS and a polygenic risk score (PRS) for disease susceptibility and progression. Methods: We randomly split 2,569 COPDGene (Genetic Epidemiology of COPD) participants with whole-blood RNA sequencing into training (n = 1,945) and testing (n = 624) samples and used 468 ECLIPSE (Evaluation of COPD Longitudinally to Identify Predictive Surrogate End-points) COPD cases with microarray data for replication. We developed a TRS using penalized regression (least absolute shrinkage and selection operator) to model FEV1/FVC and studied the predictive value of TRS for COPD (Global Initiative for Chronic Obstructive Lung Disease 2-4), prospective FEV1 change (ml/yr), and additional COPD-related traits. We adjusted for potential confounders, including age and smoking. We evaluated the predictive performance of the TRS in the context of a previously derived PRS and clinical factors. Measurements and Main Results: The TRS included 147 transcripts and was associated with COPD (odds ratio, 3.3; 95% confidence interval [CI], 2.4-4.5; P < 0.001), FEV1 change (β, -17 ml/yr; 95% CI, -28 to -6.6; P = 0.002), and other COPD-related traits. In ECLIPSE cases, we replicated the association with FEV1 change (β, -8.2; 95% CI, -15 to -1; P = 0.025) and the majority of other COPD-related traits. Models including PRS, TRS, and clinical factors were more predictive of COPD (area under the receiver operator characteristic curve, 0.84) and annualized FEV1 change compared with models with one risk score or clinical factors alone. Conclusions: Blood transcriptomics can improve prediction of COPD and lung function decline when added to a PRS and clinical risk factors.

DNA Methylation Abnormalities and Altered Whole Transcriptome Profiles after Switching from Combustible Tobacco Smoking to Heated Tobacco Products

BACKGROUND

The use of heated tobacco products (HTP) has increased exponentially in Japan since 2016; however, their effects on health remain a major concern.

METHODS

Tsuruoka Metabolome Cohort Study participants (n = 11,002) were grouped on the basis of their smoking habits as never smokers (NS), past smokers (PS), combustible tobacco smokers (CS), and HTP users for <2 years. Peripheral blood mononuclear cells were collected from 52 participants per group matched to HTP users using propensity scores, and DNA and RNA were purified from the samples. DNA methylation (DNAm) analysis of the 17 smoking-associated DNAm biomarker genes (such as AHRR, F2RL3, LRRN3, and GPR15), as well as whole transcriptome analysis, was performed.

RESULTS

Ten of the 17 genes were significantly hypomethylated in CS and HTP users compared with NS, among which AHRR, F2RL3, and RARA showed intermediate characteristics between CS and NS; nonetheless, AHRR expression was significantly higher in CS than in the other three groups. Conversely, LRRN3 and GPR15 were more hypomethylated in HTP users than in NS, and GPR15 expression was markedly upregulated in all the groups when compared with that in NS.

CONCLUSIONS

HTP users (switched from CS <2 years) display abnormal DNAm and transcriptome profiles, albeit to a lesser extent than the CS. However, because the molecular genetic effects of long-term HTP use are still unknown, long-term molecular epidemiologic studies are needed.

IMPACT

This study provides new insights into the molecular genetic effects on DNAm and transcriptome profiles in HTP users who switched from CS.

LINC00599 influences smoke-related chronic obstructive pulmonary disease and regulates CSE-induced epithelial cell apoptosis and inflammation by targeting miR-212-5p/BASP1 axis

LINC00599 has been reported to be upregulated in response to cigarette smoking. However, the effect and underlying mechanism of LINC00599 in chronic obstructive pulmonary disease (COPD) are still under exploration. In this study, LINC00599 was upregulated in the COPD patients and was of clinical value to distinguish COPD patients. COPD cell models were established using 16HBE cells under cigarette smoke extract (CSE) treatment. LINC00599 levels were elevated in a dose and time-dependent way in response to CSE stimulation. The effect of LINC00599 on CSE-induced 16HBE cells was explored. The results showed that LINC00599 deficiency reversed the CSE-induced inhibition on cell viability and proliferation, and rescued the CSE-induced enhancement on cell 16HBE cell apoptosis and inflammation response. Moreover, LINC00599 bound with miR-212-5p to upregulate the BASP1 (brain abundant membrane attached signal protein 1) expression. MiR-212-5p was expressed at a low level in the tissue samples of COPD patients, and its levels were upregulated in LINC00599 silenced cells. BASP1 was targeted by miR-212-5p and its upregulation was identified in the tissue samples of COPD patients and cell models. BASP1 levels were downregulated after miR-212-5p overexpression or LINC00599 silencing. Moreover, the rescue assays demonstrated that BASP1 overexpression reversed the effect of silenced LINC00599 on 16HBE cells after CSE treatment, which indicated that LINC00599 promoted the COPD development by regulating BASP1 expression. In conclusion, LINC00599 facilitated CSE-induced cell apoptosis and inflammation response, while inhibiting the cell viability and proliferation in COPD progression via modulating miR-212-5p/BASP1 axis.

Genome-wide methylation and expression analyses reveal the epigenetic landscape of immune-related diseases for tobacco smoking

BACKGROUND

Smoking is a major causal risk factor for lung cancer, chronic obstructive pulmonary disease (COPD), cardiovascular disease (CVD), and is the main preventable cause of deaths in the world. The components of cigarette smoke are involved in immune and inflammatory processes, which may increase the prevalence of cigarette smoke-related diseases. However, the underlying molecular mechanisms linking smoking and diseases have not been well explored. This study was aimed to depict a global map of DNA methylation and gene expression changes induced by tobacco smoking and to explore the molecular mechanisms between smoking and human diseases through whole-genome bisulfite sequencing (WGBS) and RNA-sequencing (RNA-seq).

RESULTS

We performed WGBS on 72 samples (36 smokers and 36 nonsmokers) and RNA-seq on 75 samples (38 smokers and 37 nonsmokers), and cytokine immunoassay on plasma from 22 males (9 smokers and 13 nonsmokers) who were recruited from the city of Jincheng in China. By comparing the data of the two groups, we discovered a genome-wide methylation landscape of differentially methylated regions (DMRs) associated with smoking. Functional enrichment analyses revealed that both smoking-related hyper-DMR genes (DMGs) and hypo-DMGs were related to synapse-related pathways, whereas the hypo-DMGs were specifically related to cancer and addiction. The differentially expressed genes (DEGs) revealed by RNA-seq analysis were significantly enriched in the "immunosuppression" pathway. Correlation analysis of DMRs with their corresponding gene expression showed that genes affected by tobacco smoking were mostly related to immune system diseases. Finally, by comparing cytokine concentrations between smokers and nonsmokers, we found that vascular endothelial growth factor (VEGF) was significantly upregulated in smokers.

CONCLUSIONS

In sum, we found that smoking-induced DMRs have different distribution patterns in hypermethylated and hypomethylated areas between smokers and nonsmokers. We further identified and verified smoking-related DMGs and DEGs through multi-omics integration analysis of DNA methylome and transcriptome data. These findings provide us a comprehensive genomic map of the molecular changes induced by smoking which would enhance our understanding of the harms of smoking and its relationship with diseases.

Improved prediction of smoking status via isoform-aware RNA-seq deep learning models

Most predictive models based on gene expression data do not leverage information related to gene splicing, despite the fact that splicing is a fundamental feature of eukaryotic gene expression. Cigarette smoking is an important environmental risk factor for many diseases, and it has profound effects on gene expression. Using smoking status as a prediction target, we developed deep neural network predictive models using gene, exon, and isoform level quantifications from RNA sequencing data in 2,557 subjects in the COPDGene Study. We observed that models using exon and isoform quantifications clearly outperformed gene-level models when using data from 5 genes from a previously published prediction model. Whereas the test set performance of the previously published model was 0.82 in the original publication, our exon-based models including an exon-to-isoform mapping layer achieved a test set AUC (area under the receiver operating characteristic) of 0.88, which improved to an AUC of 0.94 using exon quantifications from a larger set of genes. Isoform variability is an important source of latent information in RNA-seq data that can be used to improve clinical prediction models.

Predictive models based on gene expression are already a part of medical decision making for selected situations such as early breast cancer treatment. Most of these models are based on measures that do not capture critical aspects of gene splicing, but with RNA sequencing it is possible to capture some of these aspects of alternative splicing and use them to improve clinical predictions. Building on previous models to predict cigarette smoking status, we show that measures of alternative splicing significantly improve the accuracy of these predictive models.

Peripheral blood microbial signatures in current and former smokers

The human microbiome has a role in the development of multiple diseases. Individual microbiome profiles are highly personalized, though many species are shared. Understanding the relationship between the human microbiome and disease may inform future individualized treatments. We hypothesize the blood microbiome signature may be a surrogate for some lung microbial characteristics. We sought associations between the blood microbiome signature and lung-relevant host factors. Based on reads not mapped to the human genome, we detected microbial nucleic acids through secondary use of peripheral blood RNA-sequencing from 2,590 current and former smokers with and without chronic obstructive pulmonary disease (COPD) from the COPDGene study. We used the Genome Analysis Toolkit (GATK) microbial pipeline PathSeq to infer microbial profiles. We tested associations between the inferred profiles and lung disease relevant phenotypes and examined links to host gene expression pathways. We replicated our analyses using a second independent set of blood RNA-seq data from 1,065 COPDGene study subjects and performed a meta-analysis across the two studies. The four phyla with highest abundance across all subjects were Proteobacteria, Actinobacteria, Firmicutes and Bacteroidetes. In our meta-analysis, we observed associations (q-value < 0.05) between Acinetobacter, Serratia, Streptococcus and Bacillus inferred abundances and Modified Medical Research Council (mMRC) dyspnea score. Current smoking status was associated (q < 0.05) with Acinetobacter, Serratia and Cutibacterium abundance. All 12 taxa investigated were associated with at least one white blood cell distribution variable. Abundance for nine of the 12 taxa was associated with sex, and seven of the 12 taxa were associated with race. Host-microbiome interaction analysis revealed clustering of genera associated with mMRC dyspnea score and smoking status, through shared links to several host pathways. This study is the first to identify a bacterial microbiome signature in the peripheral blood of current and former smokers. Understanding the relationships between systemic microbial signatures and lung-related phenotypes may inform novel interventions and aid understanding of the systemic effects of smoking.

Analysis of the polarization states of the alveolar macrophages in chronic obstructive pulmonary disease samples based on miRNA-mRNA network signatures

Background

Multiple gene expression studies have been performed to investigate the biomarkers of chronic obstructive pulmonary disease (COPD). However, few studies have related COPD to macrophage cells.

Methods

The gene expression levels of clinical samples of COPD smokers (COPD; n=6), healthy smokers (Smoke; n=11), and never smokers (Never; n=4) were downloaded from the Gene Expression Omnibus (GEO) repository of GSE124180. The expression levels of messenger RNAs (mRNAs) and microRNAs (miRNAs) in macrophage cells of M0 (n=7), M1 (n=7), and M2 (n=7) were downloaded from the GEO repository of GSE46903 and GSE51307. Differentially expressed (DE) mRNAs (DEmRNAs) were identified by edgeR and GEO2R, with an adjusted P value <0.05 and |log₂fold change (FC)| ≥1 chosen as the cut-off threshold. The potential target genes of miRNA were identified using miRanda (v3.3a) and TargetScan (v6.0) with default settings. Gene Ontology (GO) and Reactome pathway analyses were performed.

Results

The composition of macrophages was quite different between COPD, Never, and Smoke samples. The proportion of M1 cells was lower than that of M0 and M2 cells in Smokers and COPD samples. Most of the genes specifically up-regulated in M1 are related to inflammation/immunity. The expression levels of miR-30a-5p, miR-200c-3p, miR-20b-5p, miR-199b-5p, and miR-301b-3p in M1 macrophages were all lower than that of M0. Their expression levels in M2 macrophages compared with M1 varied, with higher expression in miR-30a-5p, miR-20b-5p, and lower expression in miR-200c-3p, and miR-301b-3p. The mRNAs of the fms related receptor tyrosine kinase 1 (FLT1), cardiotrophin like cytokine factor 1 (CLCF1), phosphodiesterase 4D (PDE4D), coagulation factor III, and tissue factor (F3) were dysregulated in COPD and macrophage cells.

Conclusions

The present study mined the miRNA-mRNA signature which might play an essential role in COPD and macrophage polarization.

LINC01414/LINC00824 genetic polymorphisms in association with the susceptibility of chronic obstructive pulmonary disease

Objective

Chronic obstructive pulmonary disease (COPD) is a complicated multi-factor, multi-gene disease. Here, we aimed to assess the association of genetic polymorphisms in LINC01414/ LINC00824 and interactions with COPD susceptibility.

Methods

Three single nucleotide polymorphisms (SNPs) in LINC01414/LINC00824 was genotyped by Agena MassARRAY platform among 315 COPD patients and 314 controls. Logistic analysis adjusted by age and gender were applied to estimate the genetic contribution of selected SNPs to COPD susceptibility.

Results

LINC01414 rs699467 (OR = 0.73, 95% CI 0.56–0.94, p = 0.015) and LINC00824 rs7815944 (OR = 0.56, 95% CI 0.31–0.99, p = 0.046) might be protective factors for COPD occurrence, while LINC01414 rs298207 (OR = 2.88, 95% CI 1.31–6.31, p = 0.008) risk-allele was related to the increased risk of COPD in the whole population. Rs7815944 was associated with the reduced risk of COPD in the subjects aged > 70 years (OR = 0.29, p = 0.005). Rs6994670 (OR = 0.57, p = 0.007) contribute to a reduced COPD risk, while rs298207 (OR = 7.94, p = 0.009) was related to a higher susceptibility to COPD at age ≤ 70 years. Rs298207 (OR = 2.54, p = 0.043) and rs7815944 (OR = 0.43, p = 0.028) variants was associated COPD risk among males. Rs7815944 (OR = 0.16, p = 0.031) was related to the reduced susceptibility of COPD in former smokers. Moreover, the association between rs298207 genotype and COPD patients with dyspnea was found (OR = 0.50, p = 0.016), and rs7815944 was related to COPD patients with wheezing (OR = 0.22, p = 0.008).

Conclusion

Our finding provided further insights into LINC01414/LINC00824 polymorphisms at risk of COPD occurrence and accumulated evidence for the genetic susceptibility of COPD.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12890-021-01579-3.

Cigarette smoke-induced alterations in blood: A review of research on DNA methylation and gene expression

Worldwide, smoking remains a threat to public health, causing preventable diseases and premature mortality. Cigarette smoke is a powerful inducer of DNA methylation and gene expression alterations, which have been associated with negative health consequences. Here, we review the current knowledge on smoking-related changes in DNA methylation and gene expression in human blood samples. We identified 30 studies focused on the association between active smoking, DNA methylation modifications, and gene expression alterations. Overall, we identified 1,758 genes with differentially methylated sites (DMS) and differentially expressed genes (DEG) between smokers and nonsmokers, of which 261 were detected in multiple studies (≥4). The most frequently (≥10 studies) reported genes were AHRR, GPR15, GFI1, and RARA. Functional enrichment analysis of the 261 genes identified the aryl hydrocarbon receptor repressor and T cell pathways (T helpers 1 and 2) as influenced by smoking status. These results highlight specific genes for future mechanistic and translational research that may be associated with cigarette smoke exposure and smoking-related diseases. (PsycInfo Database Record (c) 2021 APA, all rights reserved).

Dysregulation of microRNA Modulatory Network in Abdominal Aortic Aneurysm

Abdominal artery aneurysm (AAA) refers to abdominal aortic dilatation of 3 cm or greater. AAA is frequently underdiagnosed due to often asymptomatic character of the disease, leading to elevated mortality due to aneurysm rupture. MiRNA constitute a pool of small RNAs controlling gene expression and is involved in many pathologic conditions in human. Targeted panel detecting altered expression of miRNA and genes involved in AAA would improve early diagnosis of this disease. In the presented study, we selected and analyzed miRNA and gene expression signatures in AAA patients. Next, generation sequencing was applied to obtain miRNA and gene-wide expression profiles from peripheral blood mononuclear cells in individuals with AAA and healthy controls. Differential expression analysis was performed using DESeq2 and uninformative variable elimination by partial least squares (UVE-PLS) methods. A total of 31 miRNAs and 51 genes were selected as the most promising biomarkers of AAA. Receiver operating characteristics (ROC) analysis showed good diagnostic ability of proposed biomarkers. Genes regulated by selected miRNAs were determined in silico and associated with functional terms closely related to cardiovascular and neurological diseases. Proposed biomarkers may be used for new diagnostic and therapeutic approaches in management of AAA. The findings will also contribute to the pool of knowledge about miRNA-dependent regulatory mechanisms involved in pathology of that disease.

Review article: impact of cigarette smoking on intestinal inflammation-direct and indirect mechanisms

BACKGROUND

The inflammatory bowel diseases, Crohn's disease and ulcerative colitis are related multifactorial diseases. Their pathogenesis is influenced by each individual's immune system, the environmental factors within exposome and genetic predisposition. Smoking habit is the single best-established environmental factor that influences disease phenotype, behaviour and response to therapy.

AIM

To assess current epidemiological, experimental and clinical evidence that may explain how smoking impacts on the pathogenesis of inflammatory bowel disease.

METHODS

A Medline search for 'cigarette smoking', in combination with terms including 'passive', 'second-hand', 'intestinal inflammation', 'Crohn's disease', 'ulcerative colitis', 'colitis'; 'intestinal epithelium', 'immune system', 'intestinal microbiota', 'tight junctions', 'mucus', 'goblet cells', 'Paneth cells', 'autophagy'; 'epigenetics', 'genes', 'DNA methylation', 'histones', 'short noncoding/long noncoding RNAs'; 'carbon monoxide/CO' and 'nitric oxide/NO' was performed.

RESULTS

Studies found evidence of direct and indirect effects of smoking on various parameters, including oxidative damage, impairment of intestinal barrier and immune cell function, epigenetic and microbiota composition changes, that contribute to the pathogenesis of inflammatory bowel disease.

CONCLUSIONS

Cigarette smoking promotes intestinal inflammation by affecting the function and interactions among intestinal epithelium, immune system and microbiota/microbiome.

Exosomal Long Non-Coding RNAs in Lung Diseases

Within the non-coding genome landscape, long non-coding RNAs (lncRNAs) and their secretion within exosomes are a window that could further explain the regulation, the sustaining, and the spread of lung diseases. We present here a compilation of the current knowledge on lncRNAs commonly found in Chronic Obstructive Pulmonary Disease (COPD), asthma, Idiopathic Pulmonary Fibrosis (IPF), or lung cancers. We built interaction networks describing the mechanisms of action for COPD, asthma, and IPF, as well as private networks for H19, MALAT1, MEG3, FENDRR, CDKN2B-AS1, TUG1, HOTAIR, and GAS5 lncRNAs in lung cancers. We identified five signaling pathways targeted by these eight lncRNAs over the lung diseases mentioned above. These lncRNAs were involved in ten treatment resistances in lung cancers, with HOTAIR being itself described in seven resistances. Besides, five of them were previously described as promising biomarkers for the diagnosis and prognosis of asthma, COPD, and lung cancers. Additionally, we describe the exosomal-based studies on H19, MALAT1, HOTAIR, GAS5, UCA1, lnc-MMP2-2, GAPLINC, TBILA, AGAP2-AS1, and SOX2-OT. This review concludes on the need for additional studies describing the lncRNA mechanisms of action and confirming their potential as biomarkers, as well as their involvement in resistance to treatment, especially in non-cancerous lung diseases.

Heme metabolism genes Downregulated in COPD Cachexia

INTRODUCTION

Cachexia contributes to increased mortality and reduced quality of life in Chronic Obstructive Pulmonary Disease (COPD) and may be associated with underlying gene expression changes. Our goal was to identify differential gene expression signatures associated with COPD cachexia in current and former smokers.

METHODS

We analyzed whole-blood gene expression data from participants with COPD in a discovery cohort (COPDGene, N = 400) and assessed replication (ECLIPSE, N = 114). To approximate the consensus definition using available criteria, cachexia was defined as weight-loss > 5% in the past 12 months or low body mass index (BMI) (< 20 kg/m2) and 1/3 criteria: decreased muscle strength (six-minute walk distance < 350 m), anemia (hemoglobin < 12 g/dl), and low fat-free mass index (FFMI) (< 15 kg/m2 among women and < 17 kg/m2 among men) in COPDGene. In ECLIPSE, cachexia was defined as weight-loss > 5% in the past 12 months or low BMI and 3/5 criteria: decreased muscle strength, anorexia, abnormal biochemistry (anemia or high c-reactive protein (> 5 mg/l)), fatigue, and low FFMI. Differential gene expression was assessed between cachectic and non-cachectic subjects, adjusting for age, sex, white blood cell counts, and technical covariates. Gene set enrichment analysis was performed using MSigDB.

RESULTS

The prevalence of COPD cachexia was 13.7% in COPDGene and 7.9% in ECLIPSE. Fourteen genes were differentially downregulated in cachectic versus non-cachectic COPD patients in COPDGene (FDR < 0.05) and ECLIPSE (FDR < 0.05).

DISCUSSION

Several replicated genes regulating heme metabolism were downregulated among participants with COPD cachexia. Impaired heme biosynthesis may contribute to cachexia development through free-iron buildup and oxidative tissue damage.

Omics and the Search for Blood Biomarkers in Chronic Obstructive Pulmonary Disease. Insights from COPDGene

There is an unmet need for blood biomarkers in diagnosis and prognosis of chronic obstructive pulmonary disease (COPD). The search for these biomarkers has been revolutionized by high-throughput sequencing techniques and multiplex platforms that can measure thousands of gene transcripts, proteins, or metabolites. We review COPDGene (Genetic Epidemiology of COPD) project publications that include DNA methylation, transcriptomic, proteomic, and metabolomic blood biomarkers and discuss their impact on COPD. Key contributions from COPDGene include identification of DNA methylation effects from smoking and genetic variation, new transcriptomic signatures in the blood, identification of protein biomarkers associated with severity and progression (e.g., sRAGE [soluble receptor for advanced glycosylation end products], inflammatory cytokines IL-6 and IL-8), and identification of small molecules (ceramides and sphingomyelin) that may be pathogenic. COPDGene studies have revealed that some of the COPD genome-wide association study polymorphisms are strongly associated with blood biomarkers (e.g., rs2070600 in AGER is a pQTL [protein quantitative trait locus] for sRAGE), underscoring the importance of combining omics results. Investigators have developed molecular networks identifying lower CD4⁺ resting memory cells associated with COPD. Genes, proteins, and metabolite networks are particularly important because the explanatory value of any single molecule is small (1-10%) compared with panels of multiple markers. COPDGene has been a useful resource in the identification and validation of multiple biomarkers for COPD. These biomarkers, either combined in multiple biomarker panels or integrated with other omics data types, may lead to novel diagnostic and prognostic tests for COPD phenotypes and may be relevant for assessing novel therapies.

Integrative Genomics Analysis Identifies ACVR1B as a Candidate Causal Gene of Emphysema Distribution

Genome-wide association studies (GWAS) have identified multiple associations with emphysema apicobasal distribution (EABD), but the biological functions of these variants are unknown. To characterize the functions of EABD-associated variants, we integrated GWAS results with 1) expression quantitative trait loci (eQTL) from the Genotype Tissue Expression (GTEx) project and subjects in the COPDGene (Genetic Epidemiology of COPD) study and 2) cell type epigenomic marks from the Roadmap Epigenomics project. On the basis of these analyses, we selected a variant near ACVR1B (activin A receptor type 1B) for functional validation. SNPs from 168 loci with P values less than 5 × 10^-5 in the largest GWAS meta-analysis of EABD were analyzed. Eighty-four loci overlapped eQTL, with 12 of these loci showing greater than 80% likelihood of harboring a single, shared GWAS and eQTL causal variant. Seventeen cell types were enriched for overlap between EABD loci and Roadmap Epigenomics marks (permutation P < 0.05), with the strongest enrichment observed in CD4⁺, CD8⁺, and regulatory T cells. We selected a putative causal variant, rs7962469, associated with ACVR1B expression in lung tissue for additional functional investigation, and reporter assays confirmed allele-specific regulatory activity for this variant in human bronchial epithelial and Jurkat immune cell lines. ACVR1B expression levels exhibit a nominally significant association with emphysema distribution. EABD-associated loci are preferentially enriched in regulatory elements of multiple cell types, most notably T-cell subsets. Multiple EABD loci colocalize to regulatory elements that are active across multiple tissues and cell types, and functional analyses confirm the presence of an EABD-associated functional variant that regulates ACVR1B expression, indicating that transforming growth factor-β signaling plays a role in the EABD phenotype. Clinical trial registered with www.clinicaltrials.gov (NCT00608764).

Long Noncoding Transcriptome in Chronic Obstructive Pulmonary Disease

Chronic airway inflammation from recurring exposures to noxious environmental stimuli results in a progressive and irreversible airflow limitation and the lung parenchymal damage that characterizes chronic obstructive pulmonary disease (COPD). The large variability observed in the onset and progression of COPD is primarily driven by complex gene-environment interactions. The transcriptomic and epigenetic memory potential of lung epithelial and innate immune cells drive responses, such as mucus hyperreactivity and airway remodeling, that are tightly regulated by various molecular mechanisms, for which several candidate susceptibility genes have been described. However, the recently described noncoding RNA species, in particular the long noncoding RNAs, may also have an important role in modulating pulmonary responses to chronic inhalation of toxic substances and the development of COPD. This review outlines the features of long noncoding RNAs that have been implicated in regulating the airway inflammatory responses to cigarette smoke exposure and their possible association with COPD pathogenesis. As COPD continues to debilitate the increasingly aging population and contribute to higher morbidity and mortality rates worldwide, the search for better biomarkers and alternative therapeutic options is pivotal.

Identification of an emphysema-associated genetic variant near TGFB2 with regulatory effects in lung fibroblasts

Murine studies have linked TGF-β signaling to emphysema, and human genome-wide association studies (GWAS) studies of lung function and COPD have identified associated regions near genes in the TGF-β superfamily. However, the functional regulatory mechanisms at these loci have not been identified. We performed the largest GWAS of emphysema patterns to date, identifying 10 GWAS loci including an association peak spanning a 200 kb region downstream from TGFB2. Integrative analysis of publicly available eQTL, DNaseI, and chromatin conformation data identified a putative functional variant, rs1690789, that may regulate TGFB2 expression in human fibroblasts. Using chromatin conformation capture, we confirmed that the region containing rs1690789 contacts the TGFB2 promoter in fibroblasts, and CRISPR/Cas-9 targeted deletion of a ~ 100 bp region containing rs1690789 resulted in decreased TGFB2 expression in primary human lung fibroblasts. These data provide novel mechanistic evidence linking genetic variation affecting the TGF-β pathway to emphysema in humans.

Analysis of genetically driven alternative splicing identifies FBXO38 as a novel COPD susceptibility gene

While many disease-associated single nucleotide polymorphisms (SNPs) are associated with gene expression (expression quantitative trait loci, eQTLs), a large proportion of complex disease genome-wide association study (GWAS) variants are of unknown function. Some of these SNPs may contribute to disease by regulating gene splicing. Here, we investigate whether SNPs that are associated with alternative splicing (splice QTL or sQTL) can identify novel functions for existing GWAS variants or suggest new associated variants in chronic obstructive pulmonary disease (COPD). RNA sequencing was performed on whole blood from 376 subjects from the COPDGene Study. Using linear models, we identified 561,060 unique sQTL SNPs associated with 30,333 splice sites corresponding to 6,419 unique genes. Similarly, 708,928 unique eQTL SNPs involving 15,913 genes were detected at 10% FDR. While there is overlap between sQTLs and eQTLs, 55.3% of sQTLs are not eQTLs. Co-localization analysis revealed that 7 out of 21 loci associated with COPD (p<1x10-6) in a published GWAS have at least one shared causal variant between the GWAS and sQTL studies. Among the genes identified to have splice sites associated with top GWAS SNPs was FBXO38, in which a novel exon was discovered to be protective against COPD. Importantly, the sQTL in this locus was validated by qPCR in both blood and lung tissue, demonstrating that splice variants relevant to lung tissue can be identified in blood. Other identified genes included CDK11A and SULT1A2. Overall, these data indicate that analysis of alternative splicing can provide novel insights into disease mechanisms. In particular, we demonstrated that SNPs in a known COPD GWAS locus on chromosome 5q32 influence alternative splicing in the gene FBXO38.

Low LINC00599 expression is a poor prognostic factor in glioma

LINC00599 has been suggested to be involved in physiological and pathological processes including carcinogenesis. However, the clinical and prognostic significance of LINC00599 in glioma patients and the effect of LINC00599 on glioma cell migration and invasion remain unknown. In our results, we first observe the expression of LINC00599 in 31 types of human cancers including tumor tissues and corresponding normal tissues at The Cancer Genome Atlas (TCGA) database, and found that LINC00599 expression levels were only reduced in lower grade glioma (LGG) tissues and glioblastoma multiforme (GBM) tissues compared with normal brain tissues. Moreover, we confirmed levels of LINC00599 expression were decreased in glioma tissues and cell lines compared with matched adjacent normal tissues and normal human astrocytes (NHAs), respectively. Meanwhile, we found that glioma tissues with WHO III-IV grade exhibited lower levels of LINC00599 expression than glioma tissues with I-II grade. The survival analysis at TCGA data showed low LINC00599 expression was associated with poor disease-free survival and overall survival in glioma patients. In vitro study suggested up-regulation of LINC00599 depressed glioma cell migration and invasion through regulating epithelial–mesenchymal transition (EMT) process. In conclusion, LINC00599 acts as a tumor-suppressing long non-coding RNA (lncRNA) in glioma.

RNA-sequencing across three matched tissues reveals shared and tissue-specific gene expression and pathway signatures of COPD

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.

High-Throughput Sequencing in Respiratory, Critical Care, and Sleep Medicine Research. An Official American Thoracic Society Workshop Report

High-throughput, "next-generation" sequencing methods are now being broadly applied across all fields of biomedical research, including respiratory disease, critical care, and sleep medicine. Although there are numerous review articles and best practice guidelines related to sequencing methods and data analysis, there are fewer resources summarizing issues related to study design and interpretation, especially as applied to common, complex, nonmalignant diseases. To address these gaps, a single-day workshop was held at the American Thoracic Society meeting in May 2017, led by the American Thoracic Society Section on Genetics and Genomics. The aim of this workshop was to review the design, analysis, interpretation, and functional follow-up of high-throughput sequencing studies in respiratory, critical care, and sleep medicine research. This workshop brought together experts in multiple fields, including genetic epidemiology, biobanking, bioinformatics, and research ethics, along with physician-scientists with expertise in a range of relevant diseases. The workshop focused on application of DNA and RNA sequencing research in common chronic diseases and did not cover sequencing studies in lung cancer, monogenic diseases (e.g., cystic fibrosis), or microbiome sequencing. Participants reviewed and discussed study design, data analysis and presentation, interpretation, functional follow-up, and reporting of results. This report summarizes the main conclusions of the workshop, specifically addressing the application of these methods in respiratory, critical care, and sleep medicine research. This workshop report may serve as a resource for our research community as well as for journal editors and reviewers of sequencing-based manuscript submissions in our research field.

Genomics and response to long-term oxygen therapy in chronic obstructive pulmonary disease

Chronic obstructive pulmonary disease (COPD) is a leading cause of death worldwide, and long-term oxygen therapy has been shown to reduce mortality in COPD patients with severe hypoxemia. However, the Long-term Oxygen Treatment Trial (LOTT), a large randomized trial, found no benefit of oxygen therapy in COPD patients with moderate hypoxemia. We hypothesized that there may be differences in response to oxygen which depend on genotype or gene expression. In a genome-wide time-to-event analysis of the primary outcome of death or hospitalization in 331 subjects, 97 single nucleotide polymorphisms (SNPs) showed evidence of interaction with oxygen therapy at P < 1e-5, including 7 SNPs near arylsulfatase B (ARSB; P = 6e-6). In microarray expression profiling on 51 whole blood samples from 37 individuals, at screening and/or at 12-month follow-up, ARSB expression was associated with the primary outcome depending on oxygen treatment. The significant SNPs were conditional expression quantitative trait loci for ARSB expression. In a network analysis of genes affected by long-term oxygen, two observed clusters including 26 co-expressed genes were enriched in mitochondrial function. Using data from the observational COPDGene Study, we validated the expression of 25 of these 26 genes, plus ARSB. The effect of long-term oxygen therapy in COPD varied based on ARSB expression and genotype. ARSB has previously been shown to be associated with hypoxemia in human bronchial and colonic epithelial cells and in a mouse model. In peripheral blood, long-term oxygen treatment affected expression of mitochondrial-related genes, a biologically relevant pathway in COPD. SNPs and expression of ARSB are associated with response to long-term oxygen in COPD. The ARSB SNPs were expression quantitative trait loci depending on oxygen therapy. Genes differentially expressed by long-term oxygen were enriched in mitochondrial functions. This suggests a potential biomarker to personalize use of long-term oxygen in COPD.

Long non-coding RNA expression patterns in lung tissues of chronic cigarette smoke induced COPD mouse model

Long non-coding RNAs (lncRNAs) have critical regulatory roles in protein-coding gene expression. Aberrant expression profiles of lncRNAs have been observed in various human diseases. In this study, we investigated transcriptome profiles in lung tissues of chronic cigarette smoke (CS)-induced COPD mouse model. We found that 109 lncRNAs and 260 mRNAs were significantly differential expressed in lungs of chronic CS-induced COPD mouse model compared with control animals. GO and KEGG analyses indicated that differentially expressed lncRNAs associated protein-coding genes were mainly involved in protein processing of endoplasmic reticulum pathway, and taurine and hypotaurine metabolism pathway. The combination of high throughput data analysis and the results of qRT-PCR validation in lungs of chronic CS-induced COPD mouse model, 16HBE cells with CSE treatment and PBMC from patients with COPD revealed that NR_102714 and its associated protein-coding gene UCHL1 might be involved in the development of COPD both in mouse and human. In conclusion, our study demonstrated that aberrant expression profiles of lncRNAs and mRNAs existed in lungs of chronic CS-induced COPD mouse model. From animal models perspective, these results might provide further clues to investigate biological functions of lncRNAs and their potential target protein-coding genes in the pathogenesis of COPD.

Integrative genomics identifies new genes associated with severe COPD and emphysema

BACKGROUND

Genome-wide association studies have identified several genetic risk loci for severe chronic obstructive pulmonary disease (COPD) and emphysema. However, these studies do not fully explain disease heritability and in most cases, fail to implicate specific genes. Integrative methods that combine gene expression data with GWAS can provide more power in discovering disease-associated genes and give mechanistic insight into regulated genes.

METHODS

We applied a recently described method that imputes gene expression using reference transcriptome data to genome-wide association studies for two phenotypes (severe COPD and quantitative emphysema) and blood and lung tissue gene expression datasets. We further tested the potential causality of individual genes using multi-variant colocalization.

RESULTS

We identified seven genes significantly associated with severe COPD, and five genes significantly associated with quantitative emphysema in whole blood or lung. We validated results in independent transcriptome databases and confirmed colocalization signals for PSMA4, EGLN2, WNT3, DCBLD1, and LILRA3. Three of these genes were not located within previously reported GWAS loci for either phenotype. We also identified genetically driven pathways, including those related to immune regulation.

CONCLUSIONS

An integrative analysis of GWAS and gene expression identified novel associations with severe COPD and quantitative emphysema, and also suggested disease-associated genes in known COPD susceptibility loci.

TRIAL REGISTRATION

NCT00608764 , Registry: ClinicalTrials.gov, Date of Enrollment of First Participant: November 2007, Date Registered: January 28, 2008 (retrospectively registered); NCT00292552 , Registry: ClinicalTrials.gov, Date of Enrollment of First Participant: December 2005, Date Registered: February 14, 2006 (retrospectively registered).