Search

Abstract Background Skeletal muscle dysfunction is a common extrapulmonary manifestation of chronic obstructive pulmonary disease (COPD). Alterations in skeletal muscle myosin heavy chain expression, with reduced type I and increased type II myosin heavy chain expression, are associated with COPD severity when studied in largely male cohorts. The objectives of this study were (1) to define an abnormal myofibre proportion phenotype in both males and females with COPD and (2) to identify transcripts and transcriptional networks associated with abnormal myofibre proportion in COPD. Methods Forty‐six participants with COPD were assessed for body composition, strength, endurance and pulmonary function. Skeletal muscle biopsies from the vastus lateralis were assayed for fibre‐type distribution and cross‐sectional area via immunofluorescence microscopy and RNA‐sequenced to generate transcriptome‐wide gene expression data. Sex‐stratified k‐means clustering of type I and IIx/IIax fibre proportions was used to define abnormal myofibre proportion in participants with COPD and contrasted with previously defined criteria. Single transcripts and weighted co‐expression network analysis modules were tested for correlation with the abnormal myofibre proportion phenotype. Results Abnormal myofibre proportion was defined in males with COPD (n = 29) as 22% type IIx/IIax fibres and in females with COPD (n = 17) as 12% type IIx/IIax fibres. Half of the participants with COPD were classified as having an abnormal myofibre proportion. Participants with COPD and an abnormal myofibre proportion had lower median handgrip strength (26.1 vs. 34.0 kg, P = 0.022), 6‐min walk distance (300 vs. 353 m, P = 0.039) and forced expiratory volume in 1 s‐to‐forced vital capacity ratio (0.42 vs. 0.48, P = 0.041) compared with participants with COPD and normal myofibre proportions. Twenty‐nine transcripts were associated with abnormal myofibre proportions in participants with COPD, with the upregulated NEB, TPM1 and TPM2 genes having the largest fold differences. Co‐expression network analysis revealed that two transcript modules were significantly positively associated with the presence of abnormal myofibre proportions. One of these co‐expression modules contained genes classically associated with muscle atrophy, as well as transcripts associated with both type I and type II myofibres, and was enriched for genetic loci associated with bone mineral density. Conclusions Our findings indicate that there are significant transcriptional alterations associated with abnormal myofibre proportions in participants with COPD. Transcripts canonically associated with both type I and type IIa fibres were enriched in a co‐expression network associated with abnormal myofibre proportion, suggesting altered transcriptional regulation across multiple fibre types.

Background: Coronary artery calcium (CAC) is a marker of subclinical atherosclerosis and is a complex heritable trait with both genetic and environmental risk factors, including sex and smoking. Methods: We performed genome-wide association (GWA) analyses for CAC among all participants and stratified by sex in the COPDGene study (n = 6144 participants of European ancestry and n = 2589 participants of African ancestry) with replication in the Diabetes Heart Study (DHS). We adjusted for age, sex, current smoking status, BMI, diabetes, self-reported high blood pressure, self-reported high cholesterol, and genetic ancestry (as summarized by principal components computed within each racial group). For the significant signals from the GWA analyses, we examined the single nucleotide polymorphism (SNP) by sex interactions, stratified by smoking status (current vs. former), and tested for a SNP by smoking status interaction on CAC. Results: We identified genome-wide significant associations for CAC in the chromosome 9p21 region [CDKN2B-AS1] among all COPDGene participants (p = 7.1 × 10−14) and among males (p = 1.0 × 10−9), but the signal was not genome-wide significant among females (p = 6.4 × 10−6). For the sex stratified GWA analyses among females, the chromosome 6p24 region [PHACTR1] had a genome-wide significant association (p = 4.4 × 10−8) with CAC, but this signal was not genome-wide significant among all COPDGene participants (p = 1.7 × 10−7) or males (p = 0.03). There was a significant interaction for the SNP rs9349379 in PHACTR1 with sex (p = 0.02), but the interaction was not significant for the SNP rs10757272 in CDKN2B-AS1 with sex (p = 0.21). In addition, PHACTR1 had a stronger association with CAC among current smokers (p = 6.2 × 10−7) than former smokers (p = 7.5 × 10−3) and the SNP by smoking status interaction was marginally significant (p = 0.03). CDKN2B-AS1 had a strong association with CAC among both former (p = 7.7 × 10−8) and current smokers (p = 1.7 × 10−7) and the SNP by smoking status interaction was not significant (p = 0.40). Conclusions: Among current and former smokers of European ancestry in the COPDGene study, we identified a genome-wide significant association in the chromosome 6p24 region [PHACTR1] with CAC among females, but not among males. This region had a significant SNP by sex and SNP by smoking interaction on CAC.

Abstract Fibrosis is a leading cause of morbidity and mortality worldwide. Although fibrosis may involve different organ systems, transforming growth factor-β (TGFβ) has been established as a master regulator of fibrosis across organs. Pirfenidone and Nintedanib are the only currently-approved drugs to treat fibrosis, specifically idiopathic pulmonary fibrosis, but their mechanisms of action remain poorly understood. To identify novel drug targets and uncover potential mechanisms by which these drugs attenuate fibrosis, we performed an integrative ‘omics analysis of transcriptomic and proteomic responses to TGFβ1-stimulated lung fibroblasts. Significant findings were annotated as associated with pirfenidone and nintedanib treatment in silico via Coremine. Integrative ‘omics identified a co-expressed transcriptomic and proteomic module significantly correlated with TGFβ1 treatment that was enriched (FDR-p = 0.04) with genes associated with pirfenidone and nintedanib treatment. While a subset of genes in this module have been implicated in fibrogenesis, several novel TGFβ1 signaling targets were identified. Specifically, four genes (BASP1, HSD17B6, CDH11, and TNS1) have been associated with pirfenidone, while five genes (CLINT1, CADM1, MTDH, SYDE1, and MCTS1) have been associated with nintedanib, and MYDGF has been implicated with treatment using both drugs. Using the Clue Drug Repurposing Hub, succinic acid was highlighted as a metabolite regulated by the protein encoded by HSD17B6 . This study provides new insights into the anti-fibrotic actions of pirfenidone and nintedanib and identifies novel targets for future mechanistic studies.

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

Abstract Background The large airway epithelial barrier provides one of the first lines of defense against respiratory viruses, including SARS-CoV-2 that causes COVID-19. Substantial inter-individual variability in individual disease courses is hypothesized to be partially mediated by the differential regulation of the genes that interact with the SARS-CoV-2 virus or are involved in the subsequent host response. Here, we comprehensively investigated non-genetic and genetic factors influencing COVID-19-relevant bronchial epithelial gene expression. Methods We analyzed RNA-sequencing data from bronchial epithelial brushings obtained from uninfected individuals. We related ACE2 gene expression to host and environmental factors in the SPIROMICS cohort of smokers with and without chronic obstructive pulmonary disease (COPD) and replicated these associations in two asthma cohorts, SARP and MAST. To identify airway biology beyond ACE2 binding that may contribute to increased susceptibility, we used gene set enrichment analyses to determine if gene expression changes indicative of a suppressed airway immune response observed early in SARS-CoV-2 infection are also observed in association with host factors. To identify host genetic variants affecting COVID-19 susceptibility in SPIROMICS, we performed expression quantitative trait (eQTL) mapping and investigated the phenotypic associations of the eQTL variants. Results We found that ACE2 expression was higher in relation to active smoking, obesity, and hypertension that are known risk factors of COVID-19 severity, while an association with interferon-related inflammation was driven by the truncated, non-binding ACE2 isoform. We discovered that expression patterns of a suppressed airway immune response to early SARS-CoV-2 infection, compared to other viruses, are similar to patterns associated with obesity, hypertension, and cardiovascular disease, which may thus contribute to a COVID-19-susceptible airway environment. eQTL mapping identified regulatory variants for genes implicated in COVID-19, some of which had pheWAS evidence for their potential role in respiratory infections. Conclusions These data provide evidence that clinically relevant variation in the expression of COVID-19-related genes is associated with host factors, environmental exposures, and likely host genetic variation.

Abstract 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) ( 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

Polygenic risk scores (PRSs) aggregate the effects of genetic variants across the genome and are used to predict risk of complex diseases, such as obesity. Current PRSs only include common variants (minor allele frequency (MAF) ≥1%), whereas the contribution of rare variants in PRSs to predict disease remains unknown. Here, we examine whether augmenting the standard common variant PRS (PRScommon) with a rare variant PRS (PRSrare) improves prediction of obesity. We used genome-wide genotyped and imputed data on 451,145 European-ancestry participants of the UK Biobank, as well as whole exome sequencing (WES) data on 184,385 participants. We performed single variant analyses (for both common and rare variants) and gene-based analyses (for rare variants) for association with BMI (kg/m2), obesity (BMI ≥ 30 kg/m2), and extreme obesity (BMI ≥ 40 kg/m2). We built PRSscommon and PRSsrare using a range of methods (Clumping+Thresholding [C+T], PRS-CS, lassosum, gene-burden test). We selected the best-performing PRSs and assessed their performance in 36,757 European-ancestry unrelated participants with whole genome sequencing (WGS) data from the Trans-Omics for Precision Medicine (TOPMed) program. The best-performing PRScommon explained 10.1% of variation in BMI, and 18.3% and 22.5% of the susceptibility to obesity and extreme obesity, respectively, whereas the best-performing PRSrare explained 1.49%, and 2.97% and 3.68%, respectively. The PRSrare was associated with an increased risk of obesity and extreme obesity (ORobesity = 1.37 per SDPRS, Pobesity = 1.7x10-85; ORextremeobesity = 1.55 per SDPRS, Pextremeobesity = 3.8x10-40), which was attenuated, after adjusting for PRScommon (ORobesity = 1.08 per SDPRS, Pobesity = 9.8x10-6; ORextremeobesity= 1.09 per SDPRS, Pextremeobesity = 0.02). When PRSrare and PRScommon are combined, the increase in explained variance attributed to PRSrare was small (incremental Nagelkerke R2 = 0.24% for obesity and 0.51% for extreme obesity). Consistently, combining PRSrare to PRScommon provided little improvement to the prediction of obesity (PRSrare AUC = 0.591; PRScommon AUC = 0.708; PRScombined AUC = 0.710). In summary, while rare variants show convincing association with BMI, obesity and extreme obesity, the PRSrare provides limited improvement over PRScommon in the prediction of obesity risk, based on these large populations.

Abstract Background Cachexia is associated with increased mortality risk among chronic obstructive pulmonary disease (COPD) patients. However, low body mass index (BMI) as opposed to cachexia is often used, particularly when calculating the BODE (BMI, Obstruction, Dyspnea and Exercise) index. For this reason, we examined mortality using a consensus definition and a weight-loss definition of cachexia among COPD cases and compared two new COPD severity indices with BODE. Methods In the current report, the consensus definition for cachexia incorporated weight-loss > 5% in 12-months or low BMI in addition to 3/5 of decreased muscle strength, fatigue, anorexia, low FFMI and inflammation. The weight-loss definition incorporated weight-loss > 5% or weight-loss > 2% (if low BMI) in 12-months. The low BMI component in BODE was replaced with the consensus definition to create the CODE (Consensus cachexia, Obstruction, Dyspnea and Exercise) index and the weight-loss definition to create the WODE (Weight loss, Obstruction, Dyspnea and Exercise) index. Mortality was assessed using Kaplan-Meier survival and Cox Regression. Performance of models was compared using C-statistics. Results Among 1483 COPD cases, the prevalences of cachexia by the consensus and weight-loss definitions were 4.7 and 10.4%, respectively. Cachectic patients had a greater than three-fold increased mortality by either the consensus or the weight-loss definition of cachexia independent of BMI and lung function. The CODE index predicted mortality slightly more accurately than the BODE and WODE indices. Conclusions Cachexia is associated with increased mortality among COPD patients. Monitoring cachexia using weight-loss criteria is relatively simple and predictive of mortality among COPD cases who may be missed if only low BMI is used.

Parkinson’s disease (PD) is the most common motor neurodegenerative disease, and neuromuscular function deficits associated with PD contribute to disability. Targeting these symptoms, our laboratory has previously evaluated 16-week high-intensity resistance exercise as rehabilitative training (RT) in individuals with PD. We reported significant improvements in muscle mass, neuromuscular function (strength, power, and motor unit activation), indices of neuromuscular junction integrity, total and motor scores on the unified Parkinson’s disease rating scale (UPDRS), and total and sub-scores on the 39-item PD Quality of Life Questionnaire (PDQ-39), supporting the use of RT to reverse symptoms. Our objective was to identify transcriptional networks that may contribute to RT-induced neuromuscular remodeling in PD. We generated transcriptome-wide skeletal muscle RNA-sequencing in 5 participants with PD [4M/1F, 67 ± 2 years, Hoehn and Yahr stages 2 (n = 3) and 3 (n = 2)] before and after 16-week high intensity RT to identify transcriptional networks that may in part underpin RT-induced neuromuscular remodeling in PD. Following RT, 304 genes were significantly upregulated, notably related to remodeling and nervous system/muscle development. Additionally, 402 genes, primarily negative regulators of muscle adaptation, were downregulated. We applied the recently developed Pathway-Level Information ExtractoR (PLIER) method to reveal coordinated gene programs (as latent variables, LVs) that differed in skeletal muscle among young (YA) and old (OA) healthy adults and PD (n = 12 per cohort) at baseline and in PD pre- vs. post-RT. Notably, one LV associated with angiogenesis, axon guidance, and muscle remodeling was significantly lower in PD than YA at baseline and was significantly increased by exercise. A different LV annotated to denervation, autophagy, and apoptosis was increased in both PD and OA relative to YA and was also reduced by 16-week RT in PD. Thus, this analysis identified two novel skeletal muscle transcriptional programs that are dysregulated by PD and aging, respectively. Notably, RT has a normalizing effect on both programs in individuals with PD. These results identify potential molecular transducers of the RT-induced improvements in neuromuscular remodeling and motor function that may aid in optimizing exercise rehabilitation strategies for individuals with PD.

Chromosomal deletions or reciprocal duplications of the 16p13.1 region have been implicated in a variety of neuropsychiatric disorders such as autism, schizophrenia, epilepsies, and attention-deficit hyperactivity disorder (ADHD). In this study, we investigated the association of recurrent genomic copy number variants (CNVs) with thoracic aortic aneurysms and dissections (TAAD). By using SNP arrays to screen and comparative genomic hybridization microarrays to validate, we identified 16p13.1 duplications in 8 out of 765 patients of European descent with adult-onset TAAD compared with 4 of 4,569 controls matched for ethnicity (P = 5.0 × 10⁻⁵, OR = 12.2). The findings were replicated in an independent cohort of 467 patients of European descent with TAAD (P = 0.005, OR = 14.7). Patients with 16p13.1 duplications were more likely to harbor a second rare CNV (P = 0.012) and to present with aortic dissections (P = 0.010) than patients without duplications. Duplications of 16p13.1 were identified in 2 of 130 patients with familial TAAD, but the duplications did not segregate with TAAD in the families. MYH11, a gene known to predispose to TAAD, lies in the duplicated region of 16p13.1, and increased MYH11 expression was found in aortic tissues from TAAD patients with 16p13.1 duplications compared with control aortas. These data suggest chromosome 16p13.1 duplications confer a risk for TAAD in addition to the established risk for neuropsychiatric disorders. It also indicates that recurrent CNVs may predispose to disorders involving more than one organ system, an observation critical to the understanding of the role of recurrent CNVs in human disease and a finding that may be common to other recurrent CNVs involving multiple genes.