Castaldi, Peter J.

Most genetic variants associated with complex traits are hypothesized to regulate gene expression. To understand the genetics underlying gene expression variability, we characterized 14,324 RNA-sequencing samples from the Trans-Omics for Precision Medicine program and performed expression and splicing quantitative trait locus (e/sQTL) analyses in six tissues and cell types, including whole blood ( n = 6454) and lung ( n = 1291). We detected tens of thousands of secondary cis-e/sQTLs, showing that secondary cis-e/sQTL discovery remains unsaturated. We fine-mapped UK Biobank–derived genome-wide association study (GWAS) signals from 164 traits and identified e/sQTL colocalizations for 10,611 GWAS signals, including 7096 that colocalize with secondary e/sQTLs. Our results suggest that even larger e/sQTL analyses will uncover additional secondary e/sQTLs, further benefiting GWAS interpretation.

Abstract RATIONALE: Chronic Obstructive Pulmonary Disease (COPD) is a leading cause of mortality worldwide. We previously published the polygenic transcriptome risk score (PTRS) which uses the cumulative effect of predicted gene expression to construct genetic predictors for complex diseases, and we demonstrated the value of PTRS for COPD built on Genotype-Tissue Expression (GTEx) Lung tissue with significantly improved cross-ancestry portability. However, the postmortem collection and lack of lung disease status on GTEx samples likely affect results. We hypothesized that the performance of PTRS will be improved by using disease-relevant expression quantitative trait loci (eQTLs) to construct the score. Here, we aimed to improve performance of PTRS by leveraging disease-relevant eQTLs. METHODS: We constructed two transcriptome prediction models, Elastic Net (EN) and Prediction Using Models Informed by Chromatin conformation and Epigenomics (PUMICE), using TOPMed lung RNA-seq data from the Lung Tissue Research Consortium (LTRC), which included COPD patients undergoing lung surgery. We compared their performance with two published models built on RNA-seq from GTEx Lung tissue, GTEx-EN and GTEx-PUMICE. We first integrated each of four transcriptome models with multi-ancestry GWAS of FEV1/FVC ratio (Shrine et al. 2023) to generate transcriptome-wide association study (TWAS) results, which produced trait-associated genes. Pathway analysis was then conducted on significant TWAS genes. Finally, the PTRS was computed on TWAS genes that were included in the top 20 significantly enriched pathways. We tested performance of PTRS in race-stratified analysis for moderate-to-severe and severe COPD in COPDGene. RESULTS: Our models produced more significant TWAS genes (FDR<0.05) and had more overlaps with genes identified for FEV1/FVC ratio from Shrine et al. 2023 compared to two GTEx models (Fig. A). Additionally, our models had significantly higher AUC for predicting both moderate-to-severe and severe COPD in race-stratified analysis. The mean AUC from LTRC-EN for predicting severe COPD was 0.57 and 0.55 respectively for EUR and AFA, which was significantly higher than mean AUC of 0.54 and 0.51 from GTEx-EN (Delong p-value=8.15x10-10 for EUR and 1.42x10-3 for AFA) (Fig. B). Similarly, the PTRS derived from our models produced higher Odds Ratios (OR) per standard deviation. The mean OR from LTRC-PUMICE for moderate-to-severe COPD was 1.25 and 1.18 respectively for EUR and AFA, while it was 1.17 and 1.12 from GTEx-PUMICE (Fig. C). CONCLUSIONS: Our study demonstrates the value of leveraging disease-relevant RNA-seq to construct transcriptome prediction models and their improvement on identifying lung function genes and predicting COPD across ancestries.

Abstract Rationale Although many studies have examined gene expression in lung tissue, the gene regulatory processes underlying emphysema are still not well understood. Finding efficient nonimaging 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 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 chronic obstructive pulmonary disease–associated differential expression of the emphysema pathways found in the bulk analyses. Measurements and Main Results In the bulk lung RNA-sequencing 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.

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

Abstract Exonic variants present some of the strongest links between genotype and phenotype. However, these variants can have significant inter-individual pathogenicity differences, known as variable penetrance. In this study, we propose a model where genetically controlled mRNA splicing modulates the pathogenicity of exonic variants. By first cataloging exonic inclusion from RNA-sequencing data in GTEx V8, we find that pathogenic alleles are depleted on highly included exons. Using a large-scale phased whole genome sequencing data from the TOPMed consortium, we observe that this effect may be driven by common splice-regulatory genetic variants, and that natural selection acts on haplotype configurations that reduce the transcript inclusion of putatively pathogenic variants, especially when limiting to haploinsufficient genes. Finally, we test if this effect may be relevant for autism risk using families from the Simons Simplex Collection, but find that splicing of pathogenic alleles has a penetrance reducing effect here as well. Overall, our results indicate that common splice-regulatory variants may play a role in reducing the damaging effects of rare exonic variants.

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