Colocalization analysis of 3' UTR alternative polyadenylation quantitative trait loci reveals novel mechanisms underlying associations with lung function.

While many disease-associated single nucleotide polymorphisms (SNPs) are expression quantitative trait loci (eQTLs), a large proportion of genome-wide association study (GWAS) variants are of unknown function. Alternative polyadenylation (APA) plays an important role in posttranscriptional regulation by allowing genes to shorten or extend 3' untranslated regions (UTRs). We hypothesized that genetic variants that affect APA in lung tissue may lend insight into the function of respiratory associated GWAS loci. We generated alternative polyadenylation (apa) QTLs using RNA sequencing and whole genome sequencing on 1241 subjects from the Lung Tissue Research Consortium (LTRC) as part of the NHLBI TOPMed project. We identified 56 179 APA sites corresponding to 13 582 unique genes after filtering out APA sites with low usage. We found that a total of 8831 APA sites were associated with at least one SNP with q-value < 0.05. The genomic distribution of lead APA SNPs indicated that the majority are intronic variants (33%), followed by downstream gene variants (26%), 3' UTR variants (17%), and upstream gene variants (within 1 kb region upstream of transcriptional start site, 10%). APA sites in 193 genes colocalized with GWAS data for at least one phenotype. Genes containing the top APA sites associated with GWAS variants include membrane associated ring-CH-type finger 2 (MARCHF2), nectin cell adhesion molecule 2 (NECTIN2), and butyrophilin subfamily 3 member A2 (BTN3A2). Overall, these findings suggest that APA may be an important mechanism for genetic variants in lung function and chronic obstructive pulmonary disease (COPD).

Blood-based Transcriptomic and Proteomic Biomarkers of Emphysema.

Rationale: Emphysema is a chronic obstructive pulmonary disease phenotype with important prognostic implications. Identifying blood-based biomarkers of emphysema will facilitate early diagnosis and development of targeted therapies. Objectives: To discover blood omics biomarkers for chest computed tomography-quantified emphysema and develop predictive biomarker panels. Methods: Emphysema blood biomarker discovery was performed using differential gene expression, alternative splicing, and protein association analyses in a training sample of 2,370 COPDGene participants with available blood RNA sequencing, plasma proteomics, and clinical data. Internal validation was conducted in a COPDGene testing sample (n = 1,016), and external validation was done in the ECLIPSE study (n = 526). Because low body mass index (BMI) and emphysema often co-occur, we performed a mediation analysis to quantify the effect of BMI on gene and protein associations with emphysema. Elastic net models with bootstrapping were also developed in the training sample sequentially using clinical, blood cell proportions, RNA-sequencing, and proteomic biomarkers to predict quantitative emphysema. Model accuracy was assessed by the area under the receiver operating characteristic curves for subjects stratified into tertiles of emphysema severity. Measurements and Main Results: Totals of 3,829 genes, 942 isoforms, 260 exons, and 714 proteins were significantly associated with emphysema (false discovery rate, 5%) and yielded 11 biological pathways. Seventy-four percent of these genes and 62% of these proteins showed mediation by BMI. Our prediction models demonstrated reasonable predictive performance in both COPDGene and ECLIPSE. The highest-performing model used clinical, blood cell, and protein data (area under the receiver operating characteristic curve in COPDGene testing, 0.90; 95% confidence interval, 0.85-0.90). Conclusions: Blood transcriptome and proteome-wide analyses revealed key biological pathways of emphysema and enhanced the prediction of emphysema.

Lung Transcriptomics Links Emphysema to Barrier Dysfunction and Macrophage Subpopulations.

RATIONALE: While many studies have examined gene expression in lung tissue, the gene regulatory processes underlying emphysema are still not well understood. Finding efficient non-imaging screening methods and disease-modifying therapies has been challenging, but knowledge of the transcriptomic features of emphysema may help in this effort. OBJECTIVES: Our goals were to identify emphysema-associated biological pathways through transcriptomic analysis of bulk lung tissue, to determine the lung cell types in which these emphysema-associated pathways are altered, and to detect unique and overlapping transcriptomic signatures in blood and lung samples. METHODS: Using RNA-sequencing data from 446 samples in the Lung Tissue Research Consortium (LTRC) and 3,606 blood samples from the COPDGene study, we examined the transcriptomic features of chest computed tomography-quantified emphysema. We also leveraged publicly available lung single-cell RNA-sequencing data to identify cell types showing COPD-associated differential expression of the emphysema pathways found in the bulk analyses. MEASUREMENTS AND MAIN RESULTS: In the bulk lung RNA-seq analysis, 1,087 differentially expressed genes and 34 dysregulated pathways were significantly associated with emphysema. We observed alternative splicing of several genes and increased activity in pluripotency and cell barrier function pathways. Lung tissue and blood samples shared differentially expressed genes and biological pathways. Multiple lung cell types displayed dysregulation of epithelial barrier function pathways, and distinct pathway activities were observed among various macrophage subpopulations. CONCLUSIONS: This study identified emphysema-related changes in gene expression and alternative splicing, cell-type specific dysregulated pathways, and instances of shared pathway dysregulation between blood and lung.

Bridging the splicing gap in human genetics with long-read RNA sequencing: finding the protein isoform drivers of disease.

Aberrant splicing underlies many human diseases, including cancer, cardiovascular diseases and neurological disorders. Genome-wide mapping of splicing quantitative trait loci (sQTLs) has shown that genetic regulation of alternative splicing is widespread. However, identification of the corresponding isoform or protein products associated with disease-associated sQTLs is challenging with short-read RNA-seq, which cannot precisely characterize full-length transcript isoforms. Furthermore, contemporary sQTL interpretation often relies on reference transcript annotations, which are incomplete. Solutions to these issues may be found through integration of newly emerging long-read sequencing technologies. Long-read sequencing offers the capability to sequence full-length mRNA transcripts and, in some cases, to link sQTLs to transcript isoforms containing disease-relevant protein alterations. Here, we provide an overview of sQTL mapping approaches, the use of long-read sequencing to characterize sQTL effects on isoforms, the linkage of RNA isoforms to protein-level functions and comment on future directions in the field. Based on recent progress, long-read RNA sequencing promises to be part of the human disease genetics toolkit to discover and treat protein isoforms causing rare and complex diseases.

Clonal Somatic Mutations in Chronic Lung Diseases Are Associated with Reduced Lung Function.

Rationale: Constantly exposed to the external environment and mutagens such as tobacco smoke, human lungs have one of the highest somatic mutation rates among all human organs. However, the relationship of these mutations to lung disease and function is not known. Objectives: To identify the prevalence and significance of clonal somatic mutations in chronic lung diseases. Methods: We analyzed the clonal somatic mutations from 1,251 samples of normal and diseased noncancerous lung tissue RNA sequencing with paired whole-genome sequencing from the Lung Tissue Research Consortium. We examined the associations of somatic mutations with lung function, disease status, and computationally deconvoluted cell types in two of the most common diseases represented in our dataset, chronic obstructive pulmonary disease (COPD; 29%) and idiopathic pulmonary fibrosis (IPF; 13%). Measurements and Main Results: Clonal somatic mutational burden was associated with reduced lung function in both COPD and IPF. We identified an increased prevalence of clonal somatic mutations in individuals with IPF compared with normal control subjects and individuals with COPD independent of age and smoking status. IPF clonal somatic mutations were enriched in disease-related and airway epithelial-expressed genes such as MUC5B in IPF. Patients who were MUC5B risk variant carriers had increased odds of developing somatic mutations of MUC5B that were explained by increased expression of MUC5B. Conclusions: Our identification of an increased prevalence of clonal somatic mutation in diseased lung that correlates with airway epithelial gene expression and disease severity highlights for the first time the role of somatic mutational processes in lung disease genetics.

Activation of CD8+ T Cells in Chronic Obstructive Pulmonary Disease Lung.

Rationale: Despite the importance of inflammation in chronic obstructive pulmonary disease (COPD), the immune cell landscape in the lung tissue of patients with mild-moderate disease has not been well characterized at the single-cell and molecular level. Objectives: To define the immune cell landscape in lung tissue from patients with mild-moderate COPD at single-cell resolution. Methods: We performed single-cell transcriptomic, proteomic, and T-cell receptor repertoire analyses on lung tissue from patients with mild-moderate COPD (n = 5, Global Initiative for Chronic Obstructive Lung Disease I or II), emphysema without airflow obstruction (n = 5), end-stage COPD (n = 2), control (n = 6), or donors (n = 4). We validated in an independent patient cohort (N = 929) and integrated with the Hhip+/- murine model of COPD. Measurements and Main Results: Mild-moderate COPD lungs have increased abundance of two CD8+ T cell subpopulations: cytotoxic KLRG1+TIGIT+CX3CR1+ TEMRA (T effector memory CD45RA+) cells, and DNAM-1+CCR5+ T resident memory (TRM) cells. These CD8+ T cells interact with myeloid and alveolar type II cells via IFNG and have hyperexpanded T-cell receptor clonotypes. In an independent cohort, the CD8+KLRG1+ TEMRA cells are increased in mild-moderate COPD lung compared with control or end-stage COPD lung. Human CD8+KLRG1+ TEMRA cells are similar to CD8+ T cells driving inflammation in an aging-related murine model of COPD. Conclusions: CD8+ TEMRA cells are increased in mild-moderate COPD lung and may contribute to inflammation that precedes severe disease. Further study of these CD8+ T cells may have therapeutic implications for preventing severe COPD.

Integrating Genetics, Transcriptomics, and Proteomics in Lung Tissue to Investigate Chronic Obstructive Pulmonary Disease.

The integration of transcriptomic and proteomic data from lung tissue with chronic obstructive pulmonary disease (COPD)-associated genetic variants could provide insight into the biological mechanisms of COPD. Here, we assessed associations between lung transcriptomics and proteomics with COPD in 98 subjects from the Lung Tissue Research Consortium. Low correlations between transcriptomics and proteomics were generally observed, but higher correlations were found for COPD-associated proteins. We integrated COPD risk SNPs or SNPs near COPD-associated proteins with lung transcripts and proteins to identify regulatory cis-quantitative trait loci (QTLs). Significant expression QTLs (eQTLs) and protein QTLs (pQTLs) were found regulating multiple COPD-associated biomarkers. We investigated mediated associations from significant pQTLs through transcripts to protein levels of COPD-associated proteins. We also attempted to identify colocalized effects between COPD genome-wide association studies and eQTL and pQTL signals. Evidence was found for colocalization between COPD genome-wide association study signals and a pQTL for RHOB and an eQTL for DSP. We applied weighted gene co-expression network analysis to find consensus COPD-associated network modules. Two network modules generated by consensus weighted gene co-expression network analysis were associated with COPD with a false discovery rate lower than 0.05. One network module is related to the catenin complex, and the other module is related to plasma membrane components. In summary, multiple cis-acting determinants of transcripts and proteins associated with COPD were identified. Colocalization analysis, mediation analysis, and correlation-based network analysis of multiple omics data may identify key genes and proteins that work together to influence COPD pathogenesis.

Heterogeneity and Progression of Chronic Obstructive Pulmonary Disease: Emphysema-Predominant and Non-Emphysema-Predominant Disease.

While variation in emphysema severity between patients with chronic obstructive pulmonary disease (COPD) is well-recognized, clinically applicable definitions of the emphysema-predominant disease (EPD) and non-emphysema-predominant disease (NEPD) subtypes have not been established. To study the clinical relevance of the EPD and NEPD subtypes, we tested the association of these subtypes with prospective decline in forced expiratory volume in 1 second (FEV1) and mortality among 3,427 subjects with Global Initiative for Chronic Obstructive Lung Disease (GOLD) spirometric grade 2-4 COPD at baseline in the Genetic Epidemiology of COPD (COPDGene) Study, an ongoing national multicenter study that started in 2007. NEPD was defined as airflow obstruction with less than 5% computed tomography (CT) quantitative densitometric emphysema at -950 Hounsfield units, and EPD was defined as airflow obstruction with 10% or greater CT emphysema. Mixed-effects models for FEV1 demonstrated larger average annual FEV1 loss in EPD subjects than in NEPD subjects (-10.2 mL/year; P < 0.001), and subtype-specific associations with FEV1 decline were identified. Cox proportional hazards models showed higher risk of mortality among EPD patients versus NEPD patients (hazard ratio = 1.46, 95% confidence interval: 1.34, 1.60; P < 0.001). To determine whether the NEPD/EPD dichotomy is captured by previously described COPDGene subtypes, we used logistic regression and receiver operating characteristic (ROC) curve analysis to predict NEPD/EPD membership using these previous subtype definitions. The analysis generally showed excellent discrimination, with areas under the ROC curve greater than 0.9. The NEPD and EPD COPD subtypes capture important aspects of COPD heterogeneity and are associated with different rates of disease progression and mortality.

From pre-COPD to COPD: a Simple, Low cost and easy to IMplement (SLIM) risk calculator.

BACKGROUND: The lifetime risk of developing clinical COPD among smokers ranges from 13% to 22%. Identifying at-risk individuals who will develop overt disease in a reasonable timeframe may allow for early intervention. We hypothesised that readily available clinical and physiological variables could help identify ever-smokers at higher risk of developing chronic airflow limitation (CAL). METHODS: Among 2273 Lovelace Smokers' Cohort (LSC) participants, we included 677 (mean age 54 years) with normal spirometry at baseline and a minimum of three spirometries, each 1 year apart. Repeated spirometric measurements were used to determine incident CAL. Using logistic regression, demographics, anthropometrics, smoking history, modified Medical Research Council dyspnoea scale, St George's Respiratory Questionnaire, comorbidities and spirometry, we related variables obtained at baseline to incident CAL as defined by the Global Initiative for Chronic Obstructive Lung Disease and lower limit of normal criteria. The predictive model derived from the LSC was validated in subjects from the COPDGene study. RESULTS: Over 6.3 years, the incidence of CAL was 26 cases per 1000 person-years. The strongest independent predictors were forced expiratory volume in 1 s (FEV1)/forced vital capacity (FVC) <0.75, having smoked ≥30 pack-years, body mass index (BMI) ≤25 kg·m2 and symptoms of chronic bronchitis. Having all four predictors increased the risk of developing CAL over 6 years to 85% (area under the receiver operating characteristic curve (AUC ROC) 0.84, 95% CI 0.81-0.89). The prediction model showed similar results when applied to subjects in the COPDGene study with a follow-up period of 10 years (AUC ROC 0.77, 95% CI 0.72-0.81). CONCLUSION: In middle-aged ever-smokers, a simple predictive model with FEV1/FVC, smoking history, BMI and chronic bronchitis helps identify subjects at high risk of developing CAL.

Alpha-1 Antitrypsin MZ Heterozygosity Is an Endotype of Chronic Obstructive Pulmonary Disease.

Rationale: Multiple studies have demonstrated an increased risk of chronic obstructive pulmonary disease (COPD) in heterozygous carriers of the AAT (alpha-1 antitrypsin) Z allele. However, it is not known if MZ subjects with COPD are phenotypically different from noncarriers (MM genotype) with COPD. Objectives: To assess if MZ subjects with COPD have different clinical features compared with MM subjects with COPD. Methods: Genotypes of SERPINA1 were ascertained by using whole-genome sequencing data in three independent studies. We compared outcomes between MM subjects with COPD and MZ subjects with COPD in each study and combined the results in a meta-analysis. We performed longitudinal and survival analyses to compare outcomes in MM and MZ subjects with COPD over time. Measurements and Main Results: We included 290 MZ subjects with COPD and 6,184 MM subjects with COPD across the three studies. MZ subjects had a lower FEV1% predicted and greater quantitative emphysema on chest computed tomography scans compared with MM subjects. In a meta-analysis, the FEV1 was 3.9% lower (95% confidence interval [CI], -6.55% to -1.26%) and emphysema (the percentage of lung attenuation areas <-950 HU) was 4.14% greater (95% CI, 1.44% to 6.84%) in MZ subjects. We found one gene, PGF (placental growth factor), to be differentially expressed in lung tissue from one study between MZ subjects and MM subjects. Conclusions: Carriers of the AAT Z allele (those who were MZ heterozygous) with COPD had lower lung function and more emphysema than MM subjects with COPD. Taken with the subtle differences in gene expression between the two groups, our findings suggest that MZ subjects represent an endotype of COPD.

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.

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.

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.

Lung tissue shows divergent gene expression between chronic obstructive pulmonary disease and idiopathic pulmonary fibrosis.

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

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.

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.

Cigarette smoking-associated isoform switching and 3' UTR lengthening via alternative polyadenylation.

Cigarette smoking induces a profound transcriptomic and systemic inflammatory response. Previous studies have focused on gene level differential expression of smoking, but the genome-wide effects of smoking on alternative isoform regulation have not yet been described. We conducted RNA sequencing in whole-blood samples of 454 current and 767 former smokers in the COPDGene Study, and we analyzed the effects of smoking on differential usage of isoforms and exons. At 10% FDR, we detected 3167 differentially expressed genes, 945 differentially used isoforms and 160 differentially used exons. Isoform switch analysis revealed widespread 3' UTR lengthening associated with cigarette smoking. The lengthening of these 3' UTRs was consistent with alternative usage of distal polyadenylation sites, and these extended 3' UTR regions were significantly enriched with functional sequence elements including microRNA and RNA-protein binding sites. These findings warrant further studies on alternative polyadenylation events as potential biomarkers and novel therapeutic targets for smoking-related diseases.

Statistical considerations for the analysis of massively parallel reporter assays data.

Noncoding DNA contains gene regulatory elements that alter gene expression, and the function of these elements can be modified by genetic variation. Massively parallel reporter assays (MPRA) enable high-throughput identification and characterization of functional genetic variants, but the statistical methods to identify allelic effects in MPRA data have not been fully developed. In this study, we demonstrate how the baseline allelic imbalance in MPRA libraries can produce biased results, and we propose a novel, nonparametric, adaptive testing method that is robust to this bias. We compare the performance of this method with other commonly used methods, and we demonstrate that our novel adaptive method controls Type I error in a wide range of scenarios while maintaining excellent power. We have implemented these tests along with routines for simulating MPRA data in the Analysis Toolset for MPRA (@MPRA), an R package for the design and analyses of MPRA experiments. It is publicly available at http://github.com/redaq/atMPRA.

Lung proteomic biomarkers associated with chronic obstructive pulmonary disease.

Identifying protein biomarkers for chronic obstructive pulmonary disease (COPD) has been challenging. Most previous studies have used individual proteins or preselected protein panels measured in blood samples. Mass spectrometry proteomic studies of lung tissue have been based on small sample sizes. We used mass spectrometry proteomic approaches to discover protein biomarkers from 150 lung tissue samples representing COPD cases and controls. Top COPD-associated proteins were identified based on multiple linear regression analysis with false discovery rate (FDR) < 0.05. Correlations between pairs of COPD-associated proteins were examined. Machine learning models were also evaluated to identify potential combinations of protein biomarkers related to COPD. We identified 4,407 proteins passing quality controls. Twenty-five proteins were significantly associated with COPD at FDR < 0.05, including interleukin 33, ferritin (light chain and heavy chain), and two proteins related to caveolae (CAV1 and CAVIN1). Multiple previously reported plasma protein biomarkers for COPD were not significantly associated with proteomic analysis of COPD in lung tissue, although RAGE was borderline significant. Eleven pairs of top significant proteins were highly correlated (r > 0.8), including several strongly correlated with RAGE (EHD2 and CAVIN1). Machine learning models using Random Forests with the top 5% of protein biomarkers demonstrated reasonable accuracy (0.707) and area under the curve (0.714) for COPD prediction. Mass spectrometry-based proteomic analysis of lung tissue is a promising approach for the identification of biomarkers for COPD.

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

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