Person:

Demeo, Dawn

Objectives: To investigate the association between methylation of transposable elements Alu and long-interspersed nuclear elements (LINE-1) and lung function. Design: Cohort study. Setting: Outpatient Veterans Administration facilities in greater Boston, Massachusetts, USA. Participants: Individuals from the Veterans Administration Normative Aging Study, a longitudinal study of aging in men, evaluated between 1999 and 2007. The majority (97%) were white. Primary and secondary outcome measures: Primary predictor was methylation, assessed using PCR-pyrosequencing after bisulphite treatment. Primary outcome was lung function as assessed by spirometry, performed according to American Thoracic Society/European Respiratory Society guidelines at the same visit as the blood draws. Results: In multivariable models adjusted for age, height, body mass index (BMI), pack-years of smoking, current smoking and race, Alu hypomethylation was associated with lower forced expiratory volume in 1 s (FEV1) (β=28 ml per 1% change in Alu methylation, p=0.017) and showed a trend towards association with a lower forced vital capacity (FVC) (β=27 ml, p=0.06) and lower FEV1/FVC (β=0.3%, p=0.058). In multivariable models adjusted for age, height, BMI, pack-years of smoking, current smoking, per cent lymphocytes, race and baseline lung function, LINE-1 hypomethylation was associated with more rapid decline of FEV1 (β=6.9 ml/year per 1% change in LINE-1 methylation, p=0.005) and of FVC (β=9.6 ml/year, p=0.002). Conclusions: In multiple regression analysis, Alu hypomethylation was associated with lower lung function, and LINE-1 hypomethylation was associated with more rapid lung function decline in a cohort of older and primarily white men from North America. Future studies should aim to replicate these findings and determine if Alu or LINE-1 hypomethylation may be due to specific and modifiable environmental exposures.

Asthma is a common chronic respiratory disease characterized by airway hyperresponsiveness (AHR). The genetics of asthma have been widely studied in mouse and human, and homologous genomic regions have been associated with mouse AHR and human asthma-related phenotypes. Our goal was to identify asthma-related genes by integrating AHR associations in mouse with human genome-wide association study (GWAS) data. We used Efficient Mixed Model Association (EMMA) analysis to conduct a GWAS of baseline AHR measures from males and females of 31 mouse strains. Genes near or containing SNPs with EMMA p-values <0.001 were selected for further study in human GWAS. The results of the previously reported EVE consortium asthma GWAS meta-analysis consisting of 12,958 diverse North American subjects from 9 study centers were used to select a subset of homologous genes with evidence of association with asthma in humans. Following validation attempts in three human asthma GWAS (i.e., Sepracor/LOCCS/LODO/Illumina, GABRIEL, DAG) and two human AHR GWAS (i.e., SHARP, DAG), the Kv channel interacting protein 4 (KCNIP4) gene was identified as nominally associated with both asthma and AHR at a gene- and SNP-level. In EVE, the smallest KCNIP4 association was at rs6833065 (P-value 2.9e-04), while the strongest associations for Sepracor/LOCCS/LODO/Illumina, GABRIEL, DAG were 1.5e-03, 1.0e-03, 3.1e-03 at rs7664617, rs4697177, rs4696975, respectively. At a SNP level, the strongest association across all asthma GWAS was at rs4697177 (P-value 1.1e-04). The smallest P-values for association with AHR were 2.3e-03 at rs11947661 in SHARP and 2.1e-03 at rs402802 in DAG. Functional studies are required to validate the potential involvement of KCNIP4 in modulating asthma susceptibility and/or AHR. Our results suggest that a useful approach to identify genes associated with human asthma is to leverage mouse AHR association data.

Implementing precision medicine for complex diseases such as chronic obstructive lung disease (COPD) will require extensive use of biomarkers and an in-depth understanding of how genetic, epigenetic, and environmental variations contribute to phenotypic diversity and disease progression. A meta-analysis from two large cohorts of current and former smokers with and without COPD [SPIROMICS (N = 750); COPDGene (N = 590)] was used to identify single nucleotide polymorphisms (SNPs) associated with measurement of 88 blood proteins (protein quantitative trait loci; pQTLs). PQTLs consistently replicated between the two cohorts. Features of pQTLs were compared to previously reported expression QTLs (eQTLs). Inference of causal relations of pQTL genotypes, biomarker measurements, and four clinical COPD phenotypes (airflow obstruction, emphysema, exacerbation history, and chronic bronchitis) were explored using conditional independence tests. We identified 527 highly significant (p < 8 X 10−10) pQTLs in 38 (43%) of blood proteins tested. Most pQTL SNPs were novel with low overlap to eQTL SNPs. The pQTL SNPs explained >10% of measured variation in 13 protein biomarkers, with a single SNP (rs7041; p = 10−392) explaining 71%-75% of the measured variation in vitamin D binding protein (gene = GC). Some of these pQTLs [e.g., pQTLs for VDBP, sRAGE (gene = AGER), surfactant protein D (gene = SFTPD), and TNFRSF10C] have been previously associated with COPD phenotypes. Most pQTLs were local (cis), but distant (trans) pQTL SNPs in the ABO blood group locus were the top pQTL SNPs for five proteins. The inclusion of pQTL SNPs improved the clinical predictive value for the established association of sRAGE and emphysema, and the explanation of variance (R2) for emphysema improved from 0.3 to 0.4 when the pQTL SNP was included in the model along with clinical covariates. Causal modeling provided insight into specific pQTL-disease relationships for airflow obstruction and emphysema. In conclusion, given the frequency of highly significant local pQTLs, the large amount of variance potentially explained by pQTL, and the differences observed between pQTLs and eQTLs SNPs, we recommend that protein biomarker-disease association studies take into account the potential effect of common local SNPs and that pQTLs be integrated along with eQTLs to uncover disease mechanisms. Large-scale blood biomarker studies would also benefit from close attention to the ABO blood group.

Genome Wide Association Studies (GWAS) and expression quantitative trait locus (eQTL) analyses have identified genetic associations with a wide range of human phenotypes. However, many of these variants have weak effects and understanding their combined effect remains a challenge. One hypothesis is that multiple SNPs interact in complex networks to influence functional processes that ultimately lead to complex phenotypes, including disease states. Here we present CONDOR, a method that represents both cis- and trans-acting SNPs and the genes with which they are associated as a bipartite graph and then uses the modular structure of that graph to place SNPs into a functional context. In applying CONDOR to eQTLs in chronic obstructive pulmonary disease (COPD), we found the global network “hub” SNPs were devoid of disease associations through GWAS. However, the network was organized into 52 communities of SNPs and genes, many of which were enriched for genes in specific functional classes. We identified local hubs within each community (“core SNPs”) and these were enriched for GWAS SNPs for COPD and many other diseases. These results speak to our intuition: rather than single SNPs influencing single genes, we see groups of SNPs associated with the expression of families of functionally related genes and that disease SNPs are associated with the perturbation of those functions. These methods are not limited in their application to COPD and can be used in the analysis of a wide variety of disease processes and other phenotypic traits.

Alterations in microRNA (miRNA) expression may contribute to COPD pathogenesis. In COPD, lung fibroblast repair functions are altered in multiple ways, including extracellular mediator release. Our prior study revealed miR-503 expression is decreased in COPD lung fibroblasts, although the exact role played by miR-503 is undetermined. The current study examined a role of miR-503 in cytokine, growth factor and fibronectin production by lung fibroblasts from patients with and without COPD. Primary adult lung fibroblasts were isolated from patients with or without COPD. MiR-503 expression and interleukin (IL)-6, -8, PGE2, HGF, KGF, VEGF and fibronectin release were examined with or without inflammatory cytokines, IL-1β and tumor necrosis factor (TNF)-α. MiR-503 expression was decreased in COPD lung fibroblasts. The expression of miR-503 was positively correlated with %FVC, %FEV1, and %DLco as well as IL-6, -8, PGE2, HGF, KGF, and VEGF in the absence or presence of IL-1ß/TNF-α. In addition, IL-8 and VEGF release from COPD lung fibroblasts were increased compared to those from control. Exogenous miR-503 inhibited VEGF release from primary adult and fetal lung fibroblasts but not IL-8 release. As expected, COPD fibroblasts proliferated more slowly than control fibroblasts. MiR-503 did not affect proliferation of either control or COPD lung fibroblasts. MiR-503 inhibition of VEGF protein production and mRNA was mediated by direct binding to the 3’ untranslated region of VEGF mRNA. Endogenous miR-503 was differently regulated by exogenous stimulants associated with COPD pathogenesis, including IL-1ß/TNF-α, TGF-ß1 and PGE2. Endogenous miR-503 inhibition augmented VEGF release by IL-1ß/TNF-α and TGF-ß1 but not by PGE2, demonstrating selectivity of miR-503 regulation of VEGF. In conclusions, reduced miR-503 augments VEGF release from lung fibroblasts from patients with COPD. Since VEGF contributes to disturbed vasculature in COPD, altered miR-503 production might play a role in modulating fibroblast-mediated vascular homeostasis in COPD.