Disease-Associated Mutations That Alter the RNA Structural Ensemble

Genome-wide association studies (GWAS) often identify disease-associated mutations in intergenic and non-coding regions of the genome. Given the high percentage of the human genome that is transcribed, we postulate that for some observed associations the disease phenotype is caused by a structural rearrangement in a regulatory region of the RNA transcript. To identify such mutations, we have performed a genome-wide analysis of all known disease-associated Single Nucleotide Polymorphisms (SNPs) from the Human Gene Mutation Database (HGMD) that map to the untranslated regions (UTRs) of a gene. Rather than using minimum free energy approaches (e.g. mFold), we use a partition function calculation that takes into consideration the ensemble of possible RNA conformations for a given sequence. We identified in the human genome disease-associated SNPs that significantly alter the global conformation of the UTR to which they map. For six disease-states (Hyperferritinemia Cataract Syndrome, β-Thalassemia, Cartilage-Hair Hypoplasia, Retinoblastoma, Chronic Obstructive Pulmonary Disease (COPD), and Hypertension), we identified multiple SNPs in UTRs that alter the mRNA structural ensemble of the associated genes. Using a Boltzmann sampling procedure for sub-optimal RNA structures, we are able to characterize and visualize the nature of the conformational changes induced by the disease-associated mutations in the structural ensemble. We observe in several cases (specifically the 5′ UTRs of FTL and RB1) SNP–induced conformational changes analogous to those observed in bacterial regulatory Riboswitches when specific ligands bind. We propose that the UTR and SNP combinations we identify constitute a “RiboSNitch,” that is a regulatory RNA in which a specific SNP has a structural consequence that results in a disease phenotype. Our SNPfold algorithm can help identify RiboSNitches by leveraging GWAS data and an analysis of the mRNA structural ensemble

Genetic Association and Risk Scores in a COPD Meta-Analysis of 16,707 Subjects

The heritability of chronic obstructive pulmonary disease (COPD) cannot be fully explained by recognized genetic risk factors identified as achieving genome-wide significance. In addition, the combined contribution of genetic variation to COPD risk has not been fully explored. We sought to determine 1) whether studies of variants from previous studies of COPD or lung function in a larger sample could identify additional associated variants, particularly for severe COPD, and 2) the impact of genetic risk scores on COPD. We genotyped 3,346 single nucleotide polymorphisms (SNP) in 2,588 cases (1,803 severe COPD) and 1,782 controls from four cohorts, and performed association testing with COPD, combining these results with existing genotyping data from 6,633 cases (3,497 severe COPD) and 5,704 controls. Additionally, we developed genetic risk scores from SNPs associated with lung function and COPD and tested their discriminatory power for COPD-related measures. We identified significant associations between SNPs near PPIC (p=1.28x10-8) and PPP4R4/SERPINA1 (p=1.01x10-8) and severe COPD; the latter association may be driven by recognized variants in SERPINA1. Genetic risk scores based on SNPs previously associated with COPD and lung function had a modest ability to discriminate COPD (AUC ~0.6), and accounted for a mean 0.9-1.9% lower FEV1 percent-predicted for each additional risk allele. In a large genetic association analysis, we identified associations with severe COPD near PPIC and SERPINA1. A risk score based on combining genetic variants had modest but significant effects on risk of COPD and lung function

Nicotinic acetylcholine receptor variants associated with susceptibility to chronic obstructive pulmonary disease: a meta-analysis

<p>Abstract</p> <p>Background</p> <p>Only 10-15% of smokers develop chronic obstructive pulmonary disease (COPD) which indicates genetic susceptibility to the disease. Recent studies suggested an association between COPD and polymorphisms in <it>CHRNA </it>coding subunits of nicotinic acetylcholine receptor. Herein, we performed a meta-analysis to clarify the impact of <it>CHRNA </it>variants on COPD.</p> <p>Methods</p> <p>We searched Web of Knowledge and Medline from 1990 through June 2011 for COPD gene studies reporting variants on <it>CHRNA</it>. Pooled odds ratios (ORs) were calculated using the major allele or genotype as reference group.</p> <p>Results</p> <p>Among seven reported variants in <it>CHRNA</it>, rs1051730 was finally analyzed with sufficient studies. Totally 3460 COPD and 11437 controls from 7 individual studies were pooled-analyzed. A-allele of rs1051730 was associated with an increased risk of COPD regardless of smoking exposure (pooled OR = 1.26, 95% CI 1.18-1.34, p < 10^-5). At the genotypic level, the ORs gradually increased per A-allele (OR = 1.27 and 1.50 for GA and AA respectively, p < 10^-5). Besides, AA genotype exhibited an association with reduced FEV1% predicted (mean difference 3.51%, 95%CI 0.87-6.16%, p = 0.009) and increased risk of emphysema (OR 1.93, 95%CI 1.29-2.90, p = 0.001).</p> <p>Conclusions</p> <p>Our findings suggest that rs1051730 in <it>CHRNA </it>is a susceptibility variant for COPD, in terms of both airway obstruction and parenchyma destruction.</p

Polymorphism of SERPINE2 gene is associated with pulmonary emphysema in consecutive autopsy cases

<p>Abstract</p> <p>Background</p> <p>The <it>SERPINA1</it>, <it>SERPINA3</it>, and <it>SERPINE2 </it>genes, which encode antiproteases, have been proposed to be susceptible genes for of chronic obstructive pulmonary disease (COPD) and related phenotypes. Whether they are associated with emphysema is not known.</p> <p>Methods</p> <p>Twelve previously reported single nucleotide polymorphisms (SNPs) in <it>SERPINA1 </it>(rs8004738, rs17751769, rs709932, rs11832, rs1303, rs28929474, and rs17580), <it>SERPINA3 </it>(rs4934, rs17473, and rs1800463), and <it>SERPINE2 </it>(rs840088 and rs975278) were genotyped in samples obtained from 1,335 consecutive autopsies of elderly Japanese people. The association between these SNPs and the severity of emphysema, as assessed using macroscopic scores, was determined.</p> <p>Results</p> <p>Emphysema of more than moderate degree was detected in 189 subjects (14.1%) and showed a significant gender difference (males, 20.5% and females, 7.0%; p < 0.0001). Among the 12 examined SNPs, only rs975278 in the <it>SERPINE2 </it>gene was positively associated with emphysema. Unlike the major alleles, homozygous minor alleles of rs975278 were associated with emphysema (odds ratio (OR) = 1.54; 95% confidence interval (CI) = 1.02-2.30; p = 0.037) and the association was very prominent in smokers (OR = 2.02; 95% CI = 1.29-3.15; p = 0.002).</p> <p>Conclusions</p> <p><it>SERPINE2 </it>may be a risk factor for the development of emphysema and its association with emphysema may be stronger in smokers.</p

Body mass index change in gastrointestinal cancer and chronic obstructive pulmonary disease is associated with Dedicator of Cytokinesis 1

BACKGROUND: There have been a number of candidate gene association studies of cancer cachexia-related traits, but no genome-wide association study (GWAS) has been published to date. Cachexia presents in patients with a number of complex traits, including both cancer and COPD. The objective of the current investigation was to search for a shared genetic aetiology for change in body mass index (ΔBMI) among cancer and COPD by using GWAS data in the Framingham Heart Study. METHODS: A linear mixed effects model accounting for age, sex, and change in smoking status was used to calculate ΔBMI in participants over 40 years of age with three consecutive BMI time points (n = 4162). Four GWAS of ΔBMI using generalized estimating equations were performed among 1085 participants with a cancer diagnosis, 204 with gastrointestinal (GI) cancer, 112 with lung cancer, and 237 with COPD to test for association with 418 365 single-nucleotide polymorphisms (SNPs). RESULTS: Two SNPs reached a level of genome-wide significance (P < 5 × 10−8) with ΔBMI: (i) rs41526344 within the CNTN4 gene, among COPD cases (β = 0.13, P = 4.3 × 10−8); and (ii) rs4751240 in the gene Dedicator of Cytokinesis 1 (DOCK1) among GI cancer cases (β = 0.10, P = 1.9 × 10−8). The DOCK1 SNP association replicated in the ΔBMI GWAS among COPD cases (βmeta-analyis = 0.10, Pmeta-analyis = 9.3 × 10−10). The DOCK1 gene codes for the dedicator of cytokinesis 1 protein, which has a role in myoblast fusion. CONCLUSIONS: In sum, one statistically significant common variant in the DOCK1 gene was associated with ΔBMI in GI cancer and COPD cases providing support for at least partially shared aetiology of ΔBMI in complex diseases

Using regulatory variants to detect gene-gene interactions identifies networks of genes linked to cell immortalization

The extent to which the impact of regulatory genetic variants may depend on other factors, such as the expression levels of upstream transcription factors, remains poorly understood. Here we report a framework in which regulatory variants are first aggregated into sets, and using these as estimates of the total cis-genetic effects on a gene we model their non-additive interactions with the expression of other genes in the genome. Using 1220 lymphoblastoid cell lines across platforms and independent datasets we identify 74 genes where the impact of their regulatory variant-set is linked to the expression levels of networks of distal genes. We show that these networks are predominantly associated with tumourigenesis pathways, through which immortalised cells are able to rapidly proliferate. We consequently present an approach to define gene interaction networks underlying important cellular pathways such as cell immortalisation

Targeted sequencing of lung function loci in chronic obstructive pulmonary disease cases and controls

Chronic obstructive pulmonary disease (COPD) is the third leading cause of death worldwide; smoking is the main risk factor for COPD, but genetic factors are also relevant contributors. Genome-wide association studies (GWAS) of the lung function measures used in the diagnosis of COPD have identified a number of loci, however association signals are often broad and collectively these loci only explain a small proportion of the heritability. In order to examine the association with COPD risk of genetic variants down to low allele frequencies, to aid fine-mapping of association signals and to explain more of the missing heritability, we undertook a targeted sequencing study in 300 COPD cases and 300 smoking controls for 26 loci previously reported to be associated with lung function. We used a pooled sequencing approach, with 12 pools of 25 individuals each, enabling high depth (30x) coverage per sample to be achieved. This pooled design maximised sample size and therefore power, but led to challenges during variant-calling since sequencing error rates and minor allele frequencies for rare variants can be very similar. For this reason we employed a rigorous quality control pipeline for variant detection which included the use of 3 independent calling algorithms. In order to avoid false positive associations we also developed tests to detect variants with potential batch effects and removed them before undertaking association testing. We tested for the effects of single variants and the combined effect of rare variants within a locus. We followed up the top signals with data available (only 67% of collapsing methods signals) in 4,249 COPD cases and 11,916 smoking controls from UK Biobank. We provide suggestive evidence for the combined effect of rare variants on COPD risk in TNXB and in sliding windows within MECOM and upstream of HHIP. These findings can lead to an improved understanding of the molecular pathways involved in the development of COPD

Cluster analysis in severe emphysema subjects using phenotype and genotype data: an exploratory investigation

Background: Numerous studies have demonstrated associations between genetic markers and COPD, but results have been inconsistent. One reason may be heterogeneity in disease definition. Unsupervised learning approaches may assist in understanding disease heterogeneity. Methods: We selected 31 phenotypic variables and 12 SNPs from five candidate genes in 308 subjects in the National Emphysema Treatment Trial (NETT) Genetics Ancillary Study cohort. We used factor analysis to select a subset of phenotypic variables, and then used cluster analysis to identify subtypes of severe emphysema. We examined the phenotypic and genotypic characteristics of each cluster. Results: We identified six factors accounting for 75% of the shared variability among our initial phenotypic variables. We selected four phenotypic variables from these factors for cluster analysis: 1) post-bronchodilator FEV1 percent predicted, 2) percent bronchodilator responsiveness, and quantitative CT measurements of 3) apical emphysema and 4) airway wall thickness. K-means cluster analysis revealed four clusters, though separation between clusters was modest: 1) emphysema predominant, 2) bronchodilator responsive, with higher FEV1; 3) discordant, with a lower FEV1 despite less severe emphysema and lower airway wall thickness, and 4) airway predominant. Of the genotypes examined, membership in cluster 1 (emphysema-predominant) was associated with TGFB1 SNP rs1800470. Conclusions: Cluster analysis may identify meaningful disease subtypes and/or groups of related phenotypic variables even in a highly selected group of severe emphysema subjects, and may be useful for genetic association studies

Comparison of Epithelial Differentiation and Immune Regulatory Properties of Mesenchymal Stromal Cells Derived from Human Lung and Bone Marrow

Mesenchymal stromal cells (MSCs) reside in many organs including lung, as shown by their isolation from fetal lung tissues, bronchial stromal compartment, bronchial-alveolar lavage and transplanted lung tissues. It is still controversial whether lung MSCs can undergo mesenchymal-to-epithelial-transition (MET) and possess immune regulatory properties. To this aim, we isolated, expanded and characterized MSCs from normal adult human lung (lung-hMSCs) and compared with human bone marrow-derived MSCs (BM-hMSCs). Our results show that lung-MSCs reside at the perivascular level and do not significantly differ from BM-hMSCs in terms of immunophenotype, stemness gene profile, mesodermal differentiation potential and modulation of T, B and NK cells. However, lung-hMSCs express higher basal level of the stemness-related marker nestin and show, following in vitro treatment with retinoic acid, higher epithelial cell polarization, which is anyway partial when compared to a control epithelial bronchial cell line. Although these results question the real capability of acquiring epithelial functions by MSCs and the feasibility of MSC-based therapeutic approaches to regenerate damaged lung tissues, the characterization of this lung-hMSC population may be useful to study the involvement of stromal cell compartment in lung diseases in which MET plays a role, such as in chronic obstructive pulmonary disease and idiopathic pulmonary fibrosis

Lung function associated gene Integrator Complex subunit 12 regulates protein synthesis pathways

Background: Genetic studies of human lung function and Chronic Obstructive Pulmonary Disease have identified a highly significant and reproducible signal on 4q24. It remains unclear which of the two candidate genes within this locus may regulate lung function: GSTCD, a gene with unknown function, and/or INTS12, a member of the Integrator Complex which is currently thought to mediate 3'end processing of small nuclear RNAs.Results: We found that, in lung tissue, 4q24 polymorphisms associated with lung function correlate with INTS12 but not neighbouring GSTCD expression. In contrast to the previous reports in other species, we only observed a minor alteration of snRNA processing following INTS12 depletion. RNAseq analysis of knockdown cells instead revealed dysregulation of a core subset of genes relevant to airway biology and a robust downregulation of protein synthesis pathways. Consistent with this, protein translation was decreased in INTS12 knockdown cells. In addition, ChIPseq experiments demonstrated INTS12 binding throughout the genome, which was enriched in transcriptionally active regions. Finally, we defined the INTS12 regulome which includes genes belonging to the protein synthesis pathways.Conclusion: INTS12 has functions beyond the canonical snRNA processing. We show that it regulates translation by regulating the expression of genes belonging to protein synthesis pathways. This study provides a detailed analysis of INTS12 activities on a genome-wide scale and contributes to the biology behind the genetic association for lung function at 4q24.</p