Genome-wide association and large-scale follow up identifies 16 new loci influencing lung function

Heritability and genome-wide association study of diffusing capacity of the lung

Although several genome-wide association studies (GWAS) have investigated the genetics of pulmonary ventilatory function, little is known about the genetic factors that influence gas exchange. The aim of the study was to investigate the heritability of, and genetic variants associated with the diffusing capacity of the lung.GWAS was performed on diffusing capacity of the lung measured by carbon monoxide uptake (D_LCO) and per alveolar volume (V_A) using the single-breath technique, in 8372 individuals from two population-based cohort studies, the Rotterdam Study and the Framingham Heart Study. Heritability was estimated in related (n=6246) and unrelated (n=3286) individuals.Heritability of D_LCO and D_LCO/V_A ranged between 23% and 28% in unrelated individuals and between 45% and 49% in related individuals. Meta-analysis identified a genetic variant in ADGRG6 that is significantly associated with D_LCO/V_A Gene expression analysis of ADGRG6 in human lung tissue revealed a decreased expression in patients with chronic obstructive pulmonary disease (COPD) and subjects with decreased D_LCO/V_AD_LCO and D_LCO/V_A are heritable traits, with a considerable proportion of variance explained by genetics. A functional variant in ADGRG6 gene region was significantly associated with D_LCO/V_A Pulmonary ADGRG6 expression was decreased in patients with COPD.

Influence of genetic variants on childhood lung function - The Generation R Study

Genetic variants associated with adult lung function could already exert the effects on childhood lung function. We aimed to examine the associations of adult lung function-related genetic variants with childhood lung function and asthma, and whether these associations were modified by atopic predisposition, tobacco smoke exposure, or early growth characteristics. In a population-based prospective cohort study among 3347 children, we selected 7 and 20 single nucleotide polymorphisms (SNPs) associated with adult forced expiratory volume in 1 second (FEV₁ ) and FEV₁ /forced vital capacity (FEV₁ /FVC), respectively. Weighted genetic risk scores (GRSs) for FEV₁ and FEV₁ /FVC were constructed. At age 10, FEV₁ , FVC, FEV₁ /FVC, forced expiratory flow between 25% and 75% (FEF_25-75 ), and forced expiratory flow at 75% (FEF₇₅ ) of FVC were measured, and information on asthma was obtained by parental-reported questionnaires. The FEV₁ -GRS was associated with lower childhood FEV₁ , FEV₁ /FVC, and FEF₇₅ (Z-score (95% CI): -0.03 (-0.05, -0.01), -0.03 (-0.05, -0.01), and -0.04 (-0.05, -0.01), respectively, per additional risk allele). The FEV₁ /FVC-GRS was associated with lower childhood FEV₁ /FVC and FEF₇₅ (Z-score (95% CI): -0.04 (-0.05, -0.03) and -0.03 (-0.05, -0.02), respectively, per additional risk allele). Effect estimates of FEV₁ -GRS with FEF_25-75 , FEV₁ , FEF₇₅ , and FVC, and of FEV₁ /FVC-GRS with FEV₁ /FVC and FEF_25-75 were stronger among children exposed to non-atopic mothers, smoking during pregnancy or in childhood, or those born with a lower birthweight, respectively (P-values for interaction < .05). Genetic risk scores were not associated with asthma. Adult lung function-related genetic variants were associated with childhood lung function. Maternal atopy, smoking during pregnancy or in childhood, and birthweight modified the observed effects.

Retinoic acid signaling is essential for airway smooth muscle homeostasis

Airway smooth muscle (ASM) is a dynamic and complex tissue involved in regulation of bronchomotor tone, but the molecular events essential for the maintenance of ASM homeostasis are not well understood. Observational and genome-wide association studies in humans have linked airway function to the nutritional status of vitamin A and its bioactive metabolite retinoic acid (RA). Here, we provide evidence that ongoing RA signaling is critical for the regulation of adult ASM phenotype. By using dietary, pharmacologic, and genetic models in mice and humans, we show that (a) RA signaling is active in adult ASM in the normal lung, (b) RA-deficient ASM cells are hypertrophic, hypercontractile, profibrotic, but not hyperproliferative, (c) TGF-β signaling, known to cause ASM hypertrophy and airway fibrosis in human obstructive lung diseases, is hyperactivated in RA-deficient ASM, (d) pharmacologic and genetic inhibition of the TGF-β activity in ASM prevents the development of the aberrant phenotype induced by RA deficiency, and (e) the consequences of transient RA deficiency in ASM are long-lasting. These results indicate that RA signaling actively maintains adult ASM homeostasis, and disruption of RA signaling leads to aberrant ASM phenotypes similar to those seen in human chronic airway diseases such as asthma.

Genome-wide association study of lung function and clinical implication in heavy smokers

BACKGROUND

The aim of this study is to identify genetic loci associated with post-bronchodilator FEV1/FVC and FEV1, and develop a multi-gene predictive model for lung function in COPD.

METHODS

Genome-wide association study (GWAS) of post-bronchodilator FEV1/FVC and FEV1 was performed in 1645 non-Hispanic White European descent smokers.

RESULTS

A functional rare variant in SERPINA1 (rs28929474: Glu342Lys) was significantly associated with post-bronchodilator FEV1/FVC (p = 1.2 × 10- 8) and FEV1 (p = 2.1 × 10- 9). In addition, this variant was associated with COPD (OR = 2.3; p = 7.8 × 10- 4) and severity (OR = 4.1; p = 0.0036). Heterozygous subjects (CT genotype) had significantly lower lung function and higher percentage of COPD and more severe COPD than subjects with the CC genotype. 8.6% of the variance of post-bronchodilator FEV1/FVC can be explained by SNPs in 10 genes with age, sex, and pack-years of cigarette smoking (P <  2.2 × 10- 16).

CONCLUSIONS

This study is the first to show genome-wide significant association of rs28929474 in SERPINA1 with lung function. Of clinical importance, heterozygotes of rs28929474 (4.7% of subjects) have significantly reduced pulmonary function, demonstrating a major impact in smokers. The multi-gene model is significantly associated with CT-based emphysema and clinical outcome measures of severity. Combining genetic information with demographic and environmental factors will further increase the predictive power for assessing reduced lung function and COPD severity.

Multiethnic meta-analysis identifies ancestry-specific and cross-ancestry loci for pulmonary function

Nearly 100 loci have been identified for pulmonary function, almost exclusively in studies of European ancestry populations. We extend previous research by meta-analyzing genome-wide association studies of 1000 Genomes imputed variants in relation to pulmonary function in a multiethnic population of 90,715 individuals of European (N = 60,552), African (N = 8429), Asian (N = 9959), and Hispanic/Latino (N = 11,775) ethnicities. We identify over 50 additional loci at genome-wide significance in ancestry-specific or multiethnic meta-analyses. Using recent fine-mapping methods incorporating functional annotation, gene expression, and differences in linkage disequilibrium between ethnicities, we further shed light on potential causal variants and genes at known and newly identified loci. Several of the novel genes encode proteins with predicted or established drug targets, including KCNK2 and CDK12. Our study highlights the utility of multiethnic and integrative genomics approaches to extend existing knowledge of the genetics of lung function and clinical relevance of implicated loci.

A Genome-Wide Association Study in Hispanics/Latinos Identifies Novel Signals for Lung Function. The Hispanic Community Health Study/Study of Latinos

RATIONALE

Lung function and chronic obstructive pulmonary disease (COPD) are heritable traits. Genome-wide association studies (GWAS) have identified numerous pulmonary function and COPD loci, primarily in cohorts of European ancestry.

OBJECTIVES

Perform a GWAS of COPD phenotypes in Hispanic/Latino populations to identify loci not previously detected in European populations.

METHODS

GWAS of lung function and COPD in Hispanic/Latino participants from a population-based cohort. We performed replication studies of novel loci in independent studies.

MEASUREMENTS AND MAIN RESULTS

Among 11,822 Hispanic/Latino participants, we identified eight novel signals; three replicated in independent populations of European Ancestry. A novel locus for FEV1 in ZSWIM7 (rs4791658; P = 4.99 × 10-9) replicated. A rare variant (minor allele frequency = 0.002) in HAL (rs145174011) was associated with FEV1/FVC (P = 9.59 × 10-9) in a region previously identified for COPD-related phenotypes; it remained significant in conditional analyses but did not replicate. Admixture mapping identified a novel region, with a variant in AGMO (rs41331850), associated with Amerindian ancestry and FEV1, which replicated. A novel locus for FEV1 identified among ever smokers (rs291231; P = 1.92 × 10-8) approached statistical significance for replication in admixed populations of African ancestry, and a novel SNP for COPD in PDZD2 (rs7709630; P = 1.56 × 10-8) regionally replicated. In addition, loci previously identified for lung function in European samples were associated in Hispanic/Latino participants in the Hispanic Community Health Study/Study of Latinos at the genome-wide significance level.

CONCLUSIONS

We identified novel signals for lung function and COPD in a Hispanic/Latino cohort. Including admixed populations when performing genetic studies may identify variants contributing to genetic etiologies of COPD.

Inhibition of RAGE Attenuates Cigarette Smoke-Induced Lung Epithelial Cell Damage via RAGE-Mediated Nrf2/DAMP Signaling

The oxidative stress and cellular apoptosis by environmental factor including cigarette smoke induces alveolar airway remodeling leading to chronic obstructive pulmonary disease (COPD). Recently, the receptor for advanced glycan end products (RAGE) which is highly expressed in alveolar epithelium is emerging as a biomarker for COPD susceptibility or progression. However, it still remains unknown how RAGE plays a role in cigarette smoke extract (CSE)-exposed human alveolar type II epithelial cell line. Therefore, we determined the efficacy of RAGE-specific antagonist FPS-ZM1 in response to CSE-induced lung epithelial cells. CSE induced the elevated generation of RONS and release of pro-inflammatory cytokines, and impaired the cellular antioxidant defense system. Further, CSE induced the alteration of RAGE distribution via the activation of redox-sensitive DAMP (Damage-associated molecular patterns) signaling through Nrf2 in cells. Although pre-treatment with SB202190 (p38 inhibitor) or SP600125 (JNK inhibitor) failed to recover the alteration of RAGE distribution, treatment of FPS-ZM1 significantly exhibited anti-inflammatory and anti-oxidative/nitrosative effects, also inhibited the activation of redox-sensitive DAMP signaling through Nrf2 (nuclear factor erythroid 2-related factor 2) migration in the presence of CSE. Taken together, our data demonstrate that RAGE and Nrf2 play a pivotal role in maintenance of alveolar epithelial integrity.

Genome-wide association study identifies loci and candidate genes for internal organ weights in Simmental beef cattle

Cattle internal organs as accessible raw materials have a long history of being widely used in beef processing, feed and pharmaceutical industry. These traits not only are of economic interest to breeders, but they are intrinsically linked to many valuable traits, such as growth, health, and productivity. Using the Illumina Bovine HD 770K SNP array, we performed a genome-wide association study for heart weight, liver weight, spleen weight, lung weight, and kidney weight in 1,217 Simmental cattle. In our research, 38 significant single nucleotide polymorphisms (SNPs) ( P < 1.49 × 10−6) were identified for five internal organ weight traits. These SNPs are within or near 13 genes, and some of them have been reported previously, including NDUFAF4, LCORL, BT.94996, SLIT2, FAM184B, LAP3, BBS12, MECOM, CD300LF, HSD17B3, TLR4, MXI1, and MB21D2. In addition, we detected four haplotype blocks on BTA6 containing 18 significant SNPs associated with spleen weight. Our results offer worthy insights into understanding the genetic mechanisms of internal organs' development, with potential application in breeding programs of Simmental beef cattle.

A candidate gene identification strategy utilizing mouse to human big-data mining: "3R-tenet" in COPD genetic research

Background

Early life impairments leading to lower lung function by adulthood are considered as risk factors for chronic obstructive pulmonary disease (COPD). Recently, we compared the lung transcriptomic profile between two mouse strains with extreme total lung capacities to identify plausible pulmonary function determining genes using microarray analysis (GSE80078). Advancement of high-throughput techniques like deep sequencing (eg. RNA-seq) and microarray have resulted in an explosion of genomic data in the online public repositories which however remains under-exploited. Strategic curation of publicly available genomic data with a mouse-human translational approach can effectively implement “3R- Tenet” by reducing screening experiments with animals and performing mechanistic studies using physiologically relevant in vitro model systems. Therefore, we sought to analyze the association of functional variations within human orthologs of mouse lung function candidate genes in a publicly available COPD lung RNA-seq data-set.

Methods

Association of missense single nucleotide polymorphisms, insertions, deletions, and splice junction variants were analyzed for susceptibility to COPD using RNA-seq data of a Korean population (GSE57148). Expression of the associated genes were studied using the Gene Paint (mouse embryo) and Human Protein Atlas (normal adult human lung) databases. The genes were also assessed for replication of the associations and expression in COPD−/mouse cigarette smoke exposed lung tissues using other datasets.

Results

Significant association (p <  0.05) of variations in 20 genes to higher COPD susceptibility have been detected within the investigated cohort. Association of HJURP, MCRS1 and TLR8 are novel in relation to COPD. The associated ADAM19 and KIT loci have been reported earlier. The remaining 15 genes have also been previously associated to COPD. Differential transcript expression levels of the associated genes in COPD- and/ or mouse emphysematous lung tissues have been detected.

Conclusion

Our findings suggest strategic mouse-human datamining approaches can identify novel COPD candidate genes using existing datasets in the online repositories. The candidates can be further evaluated for mechanistic role through in vitro studies using appropriate primary cells/cell lines. Functional studies can be limited to transgenic animal models of only well supported candidate genes. This approach will lead to a significant reduction of animal experimentation in respiratory research.

Electronic supplementary material

The online version of this article (10.1186/s12931-018-0795-y) contains supplementary material, which is available to authorized users.

Integration of multi-omics datasets enables molecular classification of COPD

Chronic obstructive pulmonary disease (COPD) is an umbrella diagnosis caused by a multitude of underlying mechanisms, and molecular sub-phenotyping is needed to develop molecular diagnostic/prognostic tools and efficacious treatments.The objective of these studies was to investigate whether multi-omics integration improves the accuracy of molecular classification of COPD in small cohorts.Nine omics data blocks (comprising mRNA, micro RNA, proteomes and metabolomes) collected from several anatomical locations from 52 female subjects were integrated by similarity network fusion (SNF). Multi-omics integration significantly improved the accuracy of group classification of COPD patients from healthy never-smokers and from smokers with normal spirometry, reducing required group sizes from n=30 to n=6 at 95% power. Seven different combinations of four to seven omics platforms achieved >95% accuracy.For the first time, a quantitative relationship between multi-omics data integration and accuracy of data-driven classification power has been demonstrated across nine omics data blocks. Integrating five to seven omics data blocks enabled 100% correct classification of COPD diagnosis with groups as small as n=6 individuals, despite strong confounding effects of current smoking. These results can serve as guidelines for the design of future systems-based multi-omics investigations, with indications that integrating five to six data blocks from several molecular levels and anatomical locations suffices to facilitate unsupervised molecular classification in small cohorts.

Understanding the role of the chromosome 15q25.1 in COPD through epigenetics and transcriptomics

Chronic obstructive pulmonary disease (COPD) is a major health burden in adults and cigarette smoking is considered the most important environmental risk factor of COPD. Chromosome 15q25.1 locus is associated with both COPD and smoking. Our study aims at understanding the mechanism underlying the association of chromosome 15q25.1 with COPD through epigenetic and transcriptional variation in a population-based setting. To assess if COPD-associated variants in 15q25.1 are methylation quantitative trait loci, epigenome-wide association analysis of four genetic variants, previously associated with COPD (P < 5 × 10-8) in the 15q25.1 locus (rs12914385:C>T-CHRNA3, rs8034191:T>C-HYKK, rs13180:C>T-IREB2 and rs8042238:C>T-IREB2), was performed in the Rotterdam study (n = 1489). All four variants were significantly associated (P < 1.4 × 10-6) with blood DNA methylation of IREB2, CHRNA3 and PSMA4, of which two, including IREB2 and PSMA4, were also differentially methylated in COPD cases and controls (P < 0.04). Further additive and multiplicative effects of smoking were evaluated and no significant effect was observed. To evaluate if these four genetic variants are expression quantitative trait loci, transcriptome-wide association analysis was performed in 1087 lung samples. All four variants were also significantly associated with differential expression of the IREB2 3'UTR in lung tissues (P < 5.4 × 10-95). We conclude that regulatory mechanisms affecting the expression of IREB2 gene, such as DNA methylation, may explain the association between genetic variants in chromosome 15q25.1 and COPD, largely independent of smoking.

Prenatal fine particulate exposure associated with reduced childhood lung function and nasal epithelia GSTP1 hypermethylation: Sex-specific effects

Background

In utero exposure to particulate matter with an aerodynamic diameter of less than 2.5 μm (PM_2.5) has been linked to child lung function. Overlapping evidence suggests that child sex and exposure timing may modify effects and associations may be mediated through glutathione S-transferase P1 (GSTP1) methylation.

Methods

We prospectively examined associations among prenatal PM_2.5 exposure and child lung function and GSTP1 methylation in an urban pregnancy cohort study. We employed a validated satellite-based spatiotemporally resolved prediction model to estimate daily prenatal PM_2.5 exposure over gestation. We used Baysian distributed lag interaction models (BDLIMs) to identify sensitive windows for prenatal PM_2.5 exposure on child lung function and nasal epithelia GSTP1 methylation at age 7 years, and to examine effect modification by child sex.

Results

BDLIMs identified a sensitive window for prenatal PM_2.5 exposure at 35–40 weeks gestation [cumulative effect estimate (CEE) = − 0.10, 95%CI = − 0.19 to − 0.01, per μg/m³ increase in PM_2.5] and at 36–40 weeks (CEE = − 0.12, 95%CI = − 0.20 to − 0.01) on FEV₁ and FVC, respectively, in boys. BDLIMs also identified a sensitive window of exposure at 37–40 weeks gestation between higher prenatal PM_2.5 exposure and increased GSTP1 percent methylation. The association between higher GSTP1 percent methylation and decreased FEV1 was borderline significant in the sample as a whole (β = − 0.37, SE = 0.20, p = 0.06) and in boys in stratified analyses (β = − 0.56, SE = 0.29, p = 0.05).

Conclusions

Prenatal PM_2.5 exposure in late pregnancy was associated with impaired early childhood lung function and hypermethylation of GSTPI in DNA isolated from nasal epithelial cells. There was a trend towards higher GSTP1 percent methylation being associated with reduced FEV1. All findings were most evident among boys.

Of pleiotropy and trajectories: Does the TGF-β pathway link childhood asthma and chronic obstructive pulmonary disease?

The study of developmental trajectories is where epigenetics truly shines. The "epi" in epigenetics captures the fact that although epigenetic processes also preside over the maintenance and termination of gene expression, the unfolding and remodeling of chromatin architecture are especially critical to prepare genes for regulated transcription. These properties imply being on a path, a trajectory to events that will occur later thanks to epigenetic programming. Thus epigenetics is about timed and timely events. In this article we discuss epigenetic and genetic evidence from several independent studies of asthma, chronic obstructive pulmonary disease, and lung function, which converge to highlight a potential role of the TGF-β gene pathway in these processes. These results raise the possibility that at least in a subset of subjects, these conditions might be functionally connected in ways that need to be further defined but that likely reflect the uniquely pleiotropic nature of TGF-β pathway genes, particularly their ability to control both lung development and immune responses essential for regulation and inflammation. Further characterization of this pathway in longitudinally phenotyped populations might unmask novel trajectories to lung disease that begin in utero and unfold into old age.

Antidepressant drug-specific prediction of depression treatment outcomes from genetic and clinical variables

Individuals with depression differ substantially in their response to treatment with antidepressants. Specific predictors explain only a small proportion of these differences. To meaningfully predict who will respond to which antidepressant, it may be necessary to combine multiple biomarkers and clinical variables. Using statistical learning on common genetic variants and clinical information in a training sample of 280 individuals randomly allocated to 12-week treatment with antidepressants escitalopram or nortriptyline, we derived models to predict remission with each antidepressant drug. We tested the reproducibility of each prediction in a validation set of 150 participants not used in model derivation. An elastic net logistic model based on eleven genetic and six clinical variables predicted remission with escitalopram in the validation dataset with area under the curve 0.77 (95%CI; 0.66-0.88; p = 0.004), explaining approximately 30% of variance in who achieves remission. A model derived from 20 genetic variables predicted remission with nortriptyline in the validation dataset with an area under the curve 0.77 (95%CI; 0.65-0.90; p < 0.001), explaining approximately 36% of variance in who achieves remission. The predictive models were antidepressant drug-specific. Validated drug-specific predictions suggest that a relatively small number of genetic and clinical variables can help select treatment between escitalopram and nortriptyline.

Association of Genetic Instrumental Variables for Lung Function on Coronary Artery Disease Risk

Background:

Lung function, assessed by forced expiratory volume in 1 second (FEV 1 ) and forced vital capacity (FVC), is inversely associated with coronary artery disease (CAD), but these associations could be because of confounding or reversed causality. We conducted a 2-sample Mendelian randomization study, using publicly available data from relevant genome-wide association studies, to examine the role of FEV 1 or FVC on CAD.

Methods:

We used the most recent genome-wide association studies on lung function to extract genetic instruments related to FEV 1 and FVC (n=92 749). Data on the association between genetic instruments and CAD were obtained from Coronary Artery Disease Genome wide Replication and Meta-analysis plus The Coronary Artery Disease Genetics 1000 Genomes-based genome-wide association studies (60 801 CAD cases and 123 504 controls). We used inverse-variance weighting with a multiplicative random effect to estimate the genetic instrumented association of FEV 1 and FVC on CAD. Sensitivity analyses included weighted median and MR-Egger methods.

Results:

Each SD greater FEV 1 was associated with a lower risk of CAD (odds ratio, 0.78 per SD; 95% confidence interval, 0.62–0.98) with a similar magnitude for FVC on CAD risk (odds ratio, 0.82 per SD; 95% confidence interval, 0.64–1.06). Estimates for FEV 1 were similar when using MR-Egger method (odds ratio, 0.80 per SD; 95% confidence interval, 0.33–1.94) although the magnitude was smaller for weighted median method (odds ratio, 0.93 per SD; 95% confidence interval, 0.75–1.17). Estimates for FVC in the sensitivity analyses were attenuated (median) or changed direction (MR-Egger).

Conclusions:

Our study suggested an inverse relation between FEV 1 and CAD, but for FVC, evidence is less clear.

Microfibril-associated glycoproteins MAGP-1 and MAGP-2 in disease

Microfibril-associated glycoproteins 1 and 2 (MAGP-1, MAGP-2) are protein components of extracellular matrix microfibrils. These proteins interact with fibrillin, the core component of microfibrils, and impart unique biological properties that influence microfibril function in vertebrates. MAGPs bind active forms of TGFβ and BMPs and are capable of modulating Notch signaling. Mutations in MAGP-1 or MAGP-2 have been linked to thoracic aneurysms and metabolic disease in humans. MAGP-2 has also been shown to be an important biomarker in several human cancers. Mice lacking MAGP-1 or MAGP-2 have defects in multiple organ systems, which reflects the widespread distribution of microfibrils in vertebrate tissues. This review summarizes our current understanding of the function of the MAGPs and their relationship to human disease.

An ancestry-based approach for detecting interactions

BACKGROUND

Epistasis and gene-environment interactions are known to contribute significantly to variation of complex phenotypes in model organisms. However, their identification in human association studies remains challenging for myriad reasons. In the case of epistatic interactions, the large number of potential interacting sets of genes presents computational, multiple hypothesis correction, and other statistical power issues. In the case of gene-environment interactions, the lack of consistently measured environmental covariates in most disease studies precludes searching for interactions and creates difficulties for replicating studies.

RESULTS

In this work, we develop a new statistical approach to address these issues that leverages genetic ancestry, defined as the proportion of ancestry derived from each ancestral population (e.g., the fraction of European/African ancestry in African Americans), in admixed populations. We applied our method to gene expression and methylation data from African American and Latino admixed individuals, respectively, identifying nine interactions that were significant at P<5×10-8. We show that two of the interactions in methylation data replicate, and the remaining six are significantly enriched for low P-values (P<1.8×10-6).

CONCLUSION

We show that genetic ancestry can be a useful proxy for unknown and unmeasured covariates in the search for interaction effects. These results have important implications for our understanding of the genetic architecture of complex traits.

Human airway branch variation and chronic obstructive pulmonary disease

Susceptibility to chronic obstructive pulmonary disease (COPD) beyond cigarette smoking is incompletely understood, although several genetic variants associated with COPD are known to regulate airway branch development. We demonstrate that in vivo central airway branch variants are present in 26.5% of the general population, are unchanged over 10 y, and exhibit strong familial aggregation. The most common airway branch variant is associated with COPD in two cohorts (n = 5,054), with greater central airway bifurcation density, and with emphysema throughout the lung. The second most common airway branch variant is associated with COPD among smokers, with narrower airway lumens in all lobes, and with genetic polymorphisms within the FGF10 gene. We conclude that central airway branch variation, readily detected by computed tomography, is a biomarker of widely altered lung structure with a genetic basis and represents a COPD susceptibility factor.

Meta-analysis of exome array data identifies six novel genetic loci for lung function

Background: Over 90 regions of the genome have been associated with lung function to date, many of which have also been implicated in chronic obstructive pulmonary disease. Methods: We carried out meta-analyses of exome array data and three lung function measures: forced expiratory volume in one second (FEV 1), forced vital capacity (FVC) and the ratio of FEV 1 to FVC (FEV 1/FVC). These analyses by the SpiroMeta and CHARGE consortia included 60,749 individuals of European ancestry from 23 studies, and 7,721 individuals of African Ancestry from 5 studies in the discovery stage, with follow-up in up to 111,556 independent individuals. Results: We identified significant (P<2·8x10 -7) associations with six SNPs: a nonsynonymous variant in RPAP1, which is predicted to be damaging, three intronic SNPs ( SEC24C, CASC17 and UQCC1) and two intergenic SNPs near to LY86 and FGF10. Expression quantitative trait loci analyses found evidence for regulation of gene expression at three signals and implicated several genes, including TYRO3 and PLAU. Conclusions: Further interrogation of these loci could provide greater understanding of the determinants of lung function and pulmonary disease.

Meta-analysis of exome array data identifies six novel genetic loci for lung function

Background: Over 90 regions of the genome have been associated with lung function to date, many of which have also been implicated in chronic obstructive pulmonary disease. Methods: We carried out meta-analyses of exome array data and three lung function measures: forced expiratory volume in one second (FEV ₁), forced vital capacity (FVC) and the ratio of FEV ₁ to FVC (FEV ₁/FVC). These analyses by the SpiroMeta and CHARGE consortia included 60,749 individuals of European ancestry from 23 studies, and 7,721 individuals of African Ancestry from 5 studies in the discovery stage, with follow-up in up to 111,556 independent individuals. Results: We identified significant (P<2·8x10 ^-7) associations with six SNPs: a nonsynonymous variant in RPAP1, which is predicted to be damaging, three intronic SNPs ( SEC24C, CASC17 and UQCC1) and two intergenic SNPs near to LY86 and FGF10. Expression quantitative trait loci analyses found evidence for regulation of gene expression at three signals and implicated several genes, including TYRO3 and PLAU. Conclusions: Further interrogation of these loci could provide greater understanding of the determinants of lung function and pulmonary disease.

Meta-analysis of exome array data identifies six novel genetic loci for lung function

Background: Over 90 regions of the genome have been associated with lung function to date, many of which have also been implicated in chronic obstructive pulmonary disease. Methods: We carried out meta-analyses of exome array data and three lung function measures: forced expiratory volume in one second (FEV 1), forced vital capacity (FVC) and the ratio of FEV 1 to FVC (FEV 1/FVC). These analyses by the SpiroMeta and CHARGE consortia included 60,749 individuals of European ancestry from 23 studies, and 7,721 individuals of African Ancestry from 5 studies in the discovery stage, with follow-up in up to 111,556 independent individuals. Results: We identified significant (P<2·8x10 -7) associations with six SNPs: a nonsynonymous variant in RPAP1, which is predicted to be damaging, three intronic SNPs ( SEC24C, CASC17 and UQCC1) and two intergenic SNPs near to LY86 and FGF10. Expression quantitative trait loci analyses found evidence for regulation of gene expression at three signals and implicated several genes, including TYRO3 and PLAU. Conclusions: Further interrogation of these loci could provide greater understanding of the determinants of lung function and pulmonary disease.

A decade of research on the 17q12-21 asthma locus: Piecing together the puzzle

Chromosome 17q12–21 remains the most highly replicated and significant asthma locus. Genotypes in the core region defined by the first genome-wide association study correlate with expression of 2 genes, ORM1-like 3 (ORMDL3) and gasdermin B (GSDMB), making these prime candidate asthma genes, although recent studies have implicated gasdermin A (GSDMA) distal to and post-GPI attachment to proteins 3 (PGAP3) proximal to the core region as independent loci. We review 10 years of studies on the 17q12–21 locus and suggest that genotype-specific risks for asthma at the proximal and distal loci are not specific to early-onset asthma and mediated by PGAP3, ORMDL3, and/or GSDMA expression. We propose that the weak and inconsistent associations of 17q single nucleotide polymorphisms with asthma in African Americans is due to the high frequency of some 17q alleles, the breakdown of linkage disequilibrium on African-derived chromosomes, and possibly different early-life asthma endotypes in these children. Finally, the inconsistent association between asthma and gene expression levels in blood or lung cells from older children and adults suggests that genotype effects may mediate asthma risk or protection during critical developmental windows and/or in response to relevant exposures in early life. Thus studies of young children and ethnically diverse populations are required to fully understand the relationship between genotype and asthma phenotype and the gene regulatory architecture at this locus. (J Allergy Clin Immunol 2018;142:749–64.)

Clinical approaches towards asthma and chronic obstructive pulmonary disease based on the heterogeneity of disease pathogenesis

Asthma and chronic obstructive pulmonary disease (COPD) are each heterogeneous disease classifications that include several clinical and pathophysiological phenotypes. This heterogeneity complicates characterization of each disease and, in some cases, hinders the selection of appropriate treatment. Therefore, in recent years, emphasis has been placed on improving our understanding of the various phenotypes of asthma and of COPD and identifying biomarkers for each phenotype. Likewise, the concept of the endotype has been gaining acceptance; an endotype is a disease subtype that is defined by unique or distinctive functional or pathophysiological mechanisms. Endotypes of asthma or COPD may be primarily characterized by increased susceptibility to type 2 inflammation, increased susceptibility to viral infections, bacterial colonization or impaired lung development. The 'Dutch hypothesis' is as follows: gene variants underlying particular endotypes interact with detrimental environmental stimuli (e.g. smoking, viral infection and air pollution) and contribute to the ultimate development of asthma, COPD or both. Novel approaches that involve multidimensional assessment should facilitate identification and management of the components that generate this heterogeneity. Ultimately, patients with chronic inflammatory lung diseases may be treated based on these endotypes as determined by the respective biomarkers that correspond to individual endotypes instead of on disease labels such as asthma, COPD or even asthma-COPD overlap syndrome (ACOS).

Multiancestry association study identifies new asthma risk loci that colocalize with immune-cell enhancer marks

We examined common variation in asthma risk by conducting a meta-analysis of worldwide asthma genome-wide association studies (23,948 cases, 118,538 controls) from ethnically-diverse populations. We identified five new asthma loci, uncovered two additional novel associations at two known asthma loci, established asthma associations at two loci implicated previously in comorbidity of asthma plus hay fever, and confirmed nine known loci. Investigation of pleiotropy showed large overlaps in genetic variants with autoimmune and inflammatory diseases. Enrichment of asthma risk loci in enhancer marks, especially in immune cells, suggests a major role of these loci in the regulation of immune-related mechanisms.

Lung function discordance in monozygotic twins and associated differences in blood DNA methylation

Background

Lung function is an important predictor of morbidity and mortality, with accelerated lung function decline reported to have immense consequences for the world's healthcare systems. The lung function decline across individual's lifetime is a consequence of age-related changes in lung anatomical structure and combination of various environmental factors; however, the exact molecular mechanisms contributing to this decline are not fully understood. DNA methylation is an epigenetic modification that changes across individual's lifetime, as well as allows for interplay between environmental and genetic factors. DNA methylation plays a crucial role in regulation of gene expression, with increasing evidence linking aberrant DNA methylation levels with a number of common human diseases. In this study, we investigated possible associations between genome-wide DNA methylation levels and lung function in 169 pairs of middle-aged monozygotic twins (86 male pairs: mean age (min-max) = 66 years (57-79); 83 female pairs: mean age (min-max) = 66 years (56-78)). The twins were collected from the Danish Twin Registry and were examined at baseline (1998-1999) and follow-up (2008-2011) visits. Using the twin design, we correlated intra-pair differences in cross-sectional and longitudinal lung function with intra-pair blood DNA methylation differences at follow-up by linear regression analyses adjusted for sex, age, BMI, smoking, and blood cell composition measured for each individual with the use of flow cytometry.

Results

We identified several differentially methylated CpG sites associated with forced expiratory volume the first second (FEV1) and forced vital capacity (FVC). Three probes identified for level of FVC were located in GLIPR1L2 gene (lowest p value = 7.14 × 10-8), involved in innate immunity and tumour-suppressor/pro-oncogenic mechanisms. Change in FEV1 during the 11-year follow-up period was associated with blood DNA methylation level in TRIM27 gene (p value = 1.55 × 10-6), a negative regulator of CD4 T cells, and also involved in cancer development. Several enriched pathways were identified, especially for FEV1, with one being "TGFBR" (Benjamini-Hochbergadjp value = 0.045), the receptor for TGFβ, a growth factor involved in normal lung tissue repair through pro-fibrotic effects.

Conclusions

Our findings suggest that epigenetic regulation of immunological- and cancer-related genes, as well as TGF-β-receptor-related genes, may be involved in the cross-sectional level and longitudinal change in lung function in middle-aged monozygotic twins.

What do polymorphisms tell us about the mechanisms of COPD?

COPD (chronic obstructive pulmonary disease) is characterized by irreversible lung airflow obstruction. Cigarette smoke is the major risk factor for COPD development. However, only a minority number of smokers develop COPD, and there are substantial variations in lung function among smokers, suggesting that genetic determinants in COPD susceptibility. During the past decade, genome-wide association studies and exome sequencing have been instrumental to identify the genetic determinants of complex traits, including COPD. Focused studies have revealed mechanisms by which genetic variants contribute to COPD and have led to novel insights in COPD pathogenesis. Through functional investigations of causal variants in COPD, from the proteinase-antiproteinase theory to emerging roles of developmental pathways (such as Hedgehog and Wnt pathways) in COPD, we have greatly expanded our understanding on this complex pulmonary disease. In this review, we critically review functional investigations on roles of genetic polymorphisms in COPD, and discuss future challenges and opportunities in discovering novel mechanisms of functional variants.

Age-related gene expression changes, and transcriptome wide association study of physical and cognitive aging traits, in the Lothian Birth Cohort 1936

Gene expression is influenced by both genetic variants and the environment. As individuals age, changes in gene expression may be associated with decline in physical and cognitive abilities. We measured transcriptome-wide expression levels in lymphoblastoid cell lines derived from members of the Lothian Birth Cohort 1936 at mean ages 70 and 76 years. Changes in gene expression levels were identified for 1,741 transcripts in 434 individuals. Gene Ontology enrichment analysis indicated an enrichment of biological processes involved in the immune system. Transcriptome-wide association analysis was performed for eleven cognitive, fitness, and biomedical aging-related traits at age 70 years (N=665 to 781) and with mortality. Transcripts for genes (F2RL3, EMILIN1 and CDC42BPA) previously identified as being differentially methylated or expressed in smoking or smoking-related cancers were overexpressed in smokers compared to non-smokers and the expression of transcripts for genes (HERPUD1, GAB2, FAM167A and GLS) previously associated with stress response, autoimmune disease and cancer were associated with telomere length. No associations between expression levels and other traits, or mortality were identified.

Genetic and Environmental Effects on Telomere Length and Lung Function: A Twin Study

Background

The purpose of the study was to estimate the heritability of leukocyte telomere length (LTL) and lung function and to examine whether LTL and lung function share genetic or environmental effects in common.

Methods

386 monozygotic and dizygotic Finnish twin sisters (age 68.4±3.4 years) were included. Relative LTL was determined from peripheral blood DNA by qPCR. Lung function measures of FEV1, FVC, FEV1/FVC, and PEF were derived from spirometry. Genetic modeling was performed with MPlus statistical software.

Results

Univariate analysis revealed that in LTL, 62% (95% confidence interval 50-72) of the variance was explained by additive genetic and 38% (28-50) by unique environmental factors. For FEV1, FVC, and PEF, the corresponding estimates were 65%-67% for additive genetic and 33%-35% for unique environmental factors. Across the sample, the phenotypic correlation between LTL and FEV1 was modest (r = .104, p = .041). Bivariate correlated factors model revealed that the genetic correlation between LTL and FEV1 was .18 (-0.19 to 0.64) and environmental correlation was -.10 (-0.84 to 0.55).

Conclusions

Both LTL and lung function variables are moderately to highly genetically determined. The associations between LTL and the lung function variables were weak. However, the positive genetic correlation point estimate gave minor suggestions that, in a larger sample, genetic factors in common might play a role in the phenotypic correlation between LTL and FEV1. Future studies with larger samples are needed to confirm these preliminary findings.

Genome-wide association study for genetic variants related with maximal voluntary ventilation reveals two novel genomic signals associated with lung function

Genome-wide association studies (GWAS) for spirometry parameters have been limited to forced vital capacity (FVC), forced expiratory volume in 1 second (FEV1), and their ratio. This study examined to identify genetic variants associated with maximal voluntary ventilation (MVV), an important spirometry parameter presenting inspiratory muscle strength.A total of 8842 Korean subjects participated in the Korean Association REsource Consortium were used to identify nucleotide variants associated with MVV and other spirometry parameters through a GWAS. Genetic associations were determined by employing a mixed model that can control background polygenic effects.The analysis revealed 3 nucleotide variants associated with MVV (P < 5 × 10). One (rs1496255) was also associated with FVC and FEV1. The other 2 variants were identified only for MVV and located in the genes of LOC102724340 (rs41434646) and FHIT (rs9833533). In particular, FHIT represses transcriptional activity of β-catenin, a critical protein for growth of skeletal muscle, and thus might have influenced the level of MVV.The current study revealed 2 novel nucleotide variants as genetic association signals for MVV. The association signals were suggested specific for neuromuscular diseases with a restrictive ventilatory impairment. Further studies are required to understand underlying mechanisms for their influence to restrictive lung diseases.

Comparison of the role of HHIP SNPs in susceptibility to chronic obstructive pulmonary disease between Chinese Han and Mongolian populations

The contribution of risk alleles to chronic obstructive pulmonary disease (COPD) may vary between populations. This study aimed to investigate the role of single nucleotide polymorphisms (SNPs) in HHIP on COPD susceptibility in the Chinese Mongolian and Han populations from Inner Mongolia autonomous region. In this case-control study, five SNPs in HHIP gene were detected in 700 COPD patients (350 Chinese Han and 350 Chinese Mongolian) and 700 healthy controls (350 Chinese Han and 350 Chinese Mongolian). The genotype, the association with COPD, and the differences between Chinese Han and Mongolian populations were evaluated using the chi-squared (χ²) test, genetic models, and logistic analysis. The minor allele C in SNP rs10519717 was associated with COPD in Mongolian (Odds ration (OR)=1.401, 95% confidence interval (CI): 1.110-1.769, P=0.044); however, not in the Han. The CC genotype in SNP rs10519717 was a risk factor for COPD in Mongolian (OR=2.667, 95% CI: 1.479-4.809, P=0.044); however, the TC genotype in Han played the same role (OR=1.396, 95% CI: 1.018-1.915, P=0.044). The GG genotype in SNP rs13147758 was protective in the Han (OR=0.546, 95% CI: 0.332-0.897, P=0.017). The homozygote of the minor alleles in SNPs rs12504628, rs1828591, and rs13118928 had a protective role in the both of COPD populations. For the minor allele distribution, the differences between the Han and Mongolian were presented only in the case group for rs12504628 (P=0.003), rs13147758 (P=0.002), rs1828591 (P=0.001), and rs13118928 (P=0.002); for the genotypes, differences was presented at the frequency of the minor allele homozygote in rs13147758 (P=0.048), rs10519717 (P=0.027), rs1828591 (P=0.041), and rs13118928 (P=0.044) in Mongolian. Our findings suggested that HHIP rs10519717 might be associated with susceptibility of Mongolian COPD. For the other SNPs, the differences between the two populations were represented by minor allele distribution and frequency of the minor allele homozygotes.

Revisiting asthma therapeutics: focus on WNT signal transduction

Asthma is a complex disease of the airways that develops as a consequence of both genetic and environmental factors. This interaction has highlighted genes important in early life, particularly those that control lung development, such as the Wingless/Integrase-1 (WNT) signalling pathway. Although aberrant WNT signalling is involved with an array of human conditions, it has received little attention within the context of asthma. Yet it is highly relevant, driving events involved with inflammation, airway remodelling, and airway hyper-responsiveness (AHR). In this review, we revisit asthma therapeutics by examining whether WNT signalling is a valid therapeutic target for asthma.

The Rotterdam Study: 2018 update on objectives, design and main results

The Rotterdam Study is a prospective cohort study ongoing since 1990 in the city of Rotterdam in The Netherlands. The study targets cardiovascular, endocrine, hepatic, neurological, ophthalmic, psychiatric, dermatological, otolaryngological, locomotor, and respiratory diseases. As of 2008, 14,926 subjects aged 45 years or over comprise the Rotterdam Study cohort. Since 2016, the cohort is being expanded by persons aged 40 years and over. The findings of the Rotterdam Study have been presented in over 1500 research articles and reports (see www.erasmus-epidemiology.nl/rotterdamstudy ). This article gives the rationale of the study and its design. It also presents a summary of the major findings and an update of the objectives and methods.

Epidemiology in Germany-general development and personal experience

Did you ever hear about epidemiology in Germany? Starting from an epidemiological desert the discipline has grown remarkably, especially during the last 10-15 years: research institutes have been established, research funding has improved, multiple curriculae in Epidemiology and Public Health are offered. This increase has been quite steep, and now the epidemiological infrastructure is much better. Several medium-sized and even big population cohorts are ongoing, and the number and quality of publications from German epidemiologists has reached a respectable level. My own career in epidemiology started in the field of environmental health. After German reunification I concentrated for many years on environmental problems in East Germany and observed the health benefits after improvement of the situation. Later, I concentrated on population-based cohorts in newborns (GINI/LISA) and adults (KORA, German National Cohort), and on biobanking. This Essay describes the development in Germany after worldwar 2, illustrated by examples of research results and build-up of epidemiological infractructures worth mentioning.

Transcriptomic analysis comparing mouse strains with extreme total lung capacities identifies novel candidate genes for pulmonary function

BACKGROUND

Failure to attain peak lung function by early adulthood is a risk factor for chronic lung diseases. Previously, we reported that C3H/HeJ mice have about twice total lung capacity (TLC) compared to JF1/MsJ mice. We identified seven lung function quantitative trait loci (QTL: Lfnq1-Lfnq7) in backcross/intercross mice derived from these inbred strains. We further demonstrated, superoxide dismutase 3, extracellular (Sod3), Kit oncogene (Kit) and secreted phosphoprotein 1 (Spp1) located on these Lfnqs as lung function determinants. Emanating from the concept of early origin of lung disease, we sought to identify novel candidate genes for pulmonary function by investigating lung transcriptome in C3H/HeJ and JF1/MsJ mice at the completion of embryonic development, bulk alveolar formation and maturity.

METHODS

Design-based stereological analysis was performed to study lung structure in C3H/HeJ and JF1/MsJ mice. Microarray was used for lung transcriptomic analysis [embryonic day 18, postnatal days 28, 70]. Quantitative real time polymerase chain reaction (qRT-PCR), western blot and immunohistochemical analysis were used to confirm selected differences.

RESULTS

Stereological analysis revealed decreased alveolar number density, elastin to collagen ratio and increased mean alveolar volume in C3H/HeJ mice compared to JF1/MsJ. Gene ontology term "extracellular region" was enriched among the decreased JF1/MsJ transcripts. Candidate genes identified using the expression-QTL strategy include: ATP-binding cassette, sub-family G (WHITE), member 1 (Abcg1), formyl peptide receptor 1 (Fpr1), gamma-aminobutyric acid (GABA) B receptor, 1 (Gabbr1); histocompatibility 2 genes: class II antigen E beta (H2-Eb1), D region locus 1 (H2-D1), and Q region locus 4 (H2-Q4); leucine rich repeat containing 6 (testis) (Lrrc6), radial spoke head 1 homolog (Rsph1), and surfactant associated 2 (Sfta2). Noteworthy genes selected as candidates for their consistent expression include: Wnt inhibitor factor 1 (Wif1), follistatin (Fst), chitinase-like 1 (Chil1), and Chil3.

CONCLUSIONS

Comparison of late embryonic, adolescent and adult lung transcript profiles between mouse strains with extreme TLCs lead to the identification of candidate genes for pulmonary function that has not been reported earlier. Further mechanistic investigations are warranted to elucidate their mode of action in determining lung function.

Lessons learned from birth cohort studies conducted in diverse environments

Childhood asthma develops from a complex interaction among host and environmental factors in early life. Birth cohort studies have provided valuable insight into asthma risk factors and the natural history of wheezing and asthma through childhood and beyond. Early life aeroallergen sensitization and wheezing illnesses associated with virus and bacterial infections have been identified as pivotal risk factors for asthma inception. Recently, focus has turned toward protective factors that promote lung health in children. Studies in a variety of environments, including farms and urban communities, suggest that diverse exposures to microbes in early life lead to a lower risk of allergy and asthma in childhood. The mechanisms underlying how these exposures and the gut and airway microbiomes alter the host response to allergens and viruses are of interest and an area of ongoing study. Longitudinal follow up of birth cohorts in diverse environments worldwide will continue to provide critical knowledge about the factors that impact the natural history of asthma.

Environment Changes Genetic Effects on Respiratory Conditions and Allergic Phenotypes

The prevalence of asthma and allergic diseases is disproportionately distributed among different populations, with an increasing trend observed in Western countries. Here we investigated how the environment affected genotype-phenotype association in a genetically homogeneous, but geographically separated population. We evaluated 18 single nucleotide polymorphisms (SNPs) corresponding to 8 genes (ADAM33, ALOX5, LT-α, LTC4S, NOS1, ORMDL3, TBXA2R and TNF-α), the lung function and five respiratory/allergic conditions (ever asthma, bronchitis, rhinitis, dermatitis and atopy) in two populations of Inuit residing either in the westernized environment of Denmark or in the rural area of Greenland. Our results showed that lung function was associated with genetic variants in ORMDL3, with polymorphisms having a significant interaction with place of residence. LT-α SNP rs909253 and rs1041981 were significantly associated with bronchitis risk. LT-α SNP rs2844484 was related to dermatitis susceptibility and was significantly influenced by the place of residence. The observed gene-phenotype relationships were exclusively present in one population and absent in the other population. We conclude that the genotype-phenotype associations relating to bronchitis and allergy susceptibility are dependent on the environment and that environmental factors/lifestyles modify genetic predisposition and change the genetic effects on diseases.

Systematic Evaluation of Pleiotropy Identifies 6 Further Loci Associated With Coronary Artery Disease

Background

Genome-wide association studies have so far identified 56 loci associated with risk of coronary artery disease (CAD). Many CAD loci show pleiotropy; that is, they are also associated with other diseases or traits.

Objectives

This study sought to systematically test if genetic variants identified for non-CAD diseases/traits also associate with CAD and to undertake a comprehensive analysis of the extent of pleiotropy of all CAD loci.

Methods

In discovery analyses involving 42,335 CAD cases and 78,240 control subjects we tested the association of 29,383 common (minor allele frequency >5%) single nucleotide polymorphisms available on the exome array, which included a substantial proportion of known or suspected single nucleotide polymorphisms associated with common diseases or traits as of 2011. Suggestive association signals were replicated in an additional 30,533 cases and 42,530 control subjects. To evaluate pleiotropy, we tested CAD loci for association with cardiovascular risk factors (lipid traits, blood pressure phenotypes, body mass index, diabetes, and smoking behavior), as well as with other diseases/traits through interrogation of currently available genome-wide association study catalogs.

Results

We identified 6 new loci associated with CAD at genome-wide significance: on 2q37 (KCNJ13-GIGYF2), 6p21 (C2), 11p15 (MRVI1-CTR9), 12q13 (LRP1), 12q24 (SCARB1), and 16q13 (CETP). Risk allele frequencies ranged from 0.15 to 0.86, and odds ratio per copy of the risk allele ranged from 1.04 to 1.09. Of 62 new and known CAD loci, 24 (38.7%) showed statistical association with a traditional cardiovascular risk factor, with some showing multiple associations, and 29 (47%) showed associations at p < 1 × 10⁻⁴ with a range of other diseases/traits.

Conclusions

We identified 6 loci associated with CAD at genome-wide significance. Several CAD loci show substantial pleiotropy, which may help us understand the mechanisms by which these loci affect CAD risk.

Epigenome-wide association study of chronic obstructive pulmonary disease and lung function in Koreans

AIM

To identify differentially methylated probes (DMPs) and regions (DMRs) in relation to chronic obstructive pulmonary disease (COPD) and lung function traits.

METHODS

We performed an epigenome-wide association study of COPD and spirometric parameters, including forced expiratory volume in 1 s (FEV1), forced vital capacity (FVC) and FEV1/FVC, in blood DNA using the Infinium HumanMethylation450 (n = 100, a Korean COPD cohort).

RESULTS

We found one significant DMP (cg03559389, DIP2C) and 104 significant DMRs after multiple-testing correction. Of these, 34 DMRs mapped to genes differential expressed with respect to the same trait. Five of the genes were associated with more than two traits: CTU2, USP36, ZNF516, KLK10 and CPT1B.

CONCLUSION

We identified novel differential methylation loci related to COPD and lung function in blood DNA in Koreans and confirmed previous findings in non-Asians. Epigenetic modification could contribute to the etiology of these phenotypes.

Genetic Association and Risk Scores in a Chronic Obstructive Pulmonary Disease 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 (SNPs) in 2,588 cases (1,803 severe COPD) and 1,782 control subjects 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 control subjects. In addition, 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.28 × 10^-8) and PPP4R4/SERPINA1 (P = 1.01 × 10^-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 (area under the curve, ∼0.6), and accounted for a mean 0.9-1.9% lower forced expiratory volume in 1 second 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.

Genetic Control of Fatty Acid β-Oxidation in Chronic Obstructive Pulmonary Disease

Bioenergetics homeostasis is important for cells to sustain normal functions and defend against injury. The genetic controls of bioenergetics homeostasis, especially lipid metabolism, remain poorly understood in chronic obstructive pulmonary disease (COPD), the third leading cause of death in the world. Additionally, the biological function of most of the susceptibility genes identified from genome-wide association studies (GWASs) in COPD remains unclear. Here, we aimed to address (1) how fatty acid oxidation (FAO), specifically β-oxidation, a key lipid metabolism pathway that provides energy to cells, contributes to cigarette smoke (CS)-induced COPD; and (2) whether-and if so, how-FAM13A (family with sequence similarity 13 member A), a well-replicated COPD GWAS gene, modulates the FAO pathway. We demonstrated that CS induced expression of carnitine palmitoyltransferase 1A (CPT1A), a key mitochondrial enzyme for the FAO pathway, thereby enhancing FAO. Pharmacological inhibition of FAO by etomoxir blunted CS-induced reactive oxygen species accumulation and cell death in lung epithelial cells. FAM13A promoted FAO, possibly by interacting with and activating sirutin 1, and increasing expression of CPT1A. Furthermore, CS-induced cell death was reduced in lungs from Fam13a-/- mice. Our results suggest that FAM13A, the COPD GWAS gene, shapes the cellular metabolic response to CS exposure by promoting the FAO pathway, which may contribute to COPD development.

Evidence for large-scale gene-by-smoking interaction effects on pulmonary function

Background

Smoking is the strongest environmental risk factor for reduced pulmonary function. The genetic component of various pulmonary traits has also been demonstrated, and at least 26 loci have been reproducibly associated with either FEV 1 (forced expiratory volume in 1 second) or FEV 1 /FVC (FEV 1 /forced vital capacity). Although the main effects of smoking and genetic loci are well established, the question of potential gene-by-smoking interaction effect remains unanswered. The aim of the present study was to assess, using a genetic risk score approach, whether the effect of these 26 loci on pulmonary function is influenced by smoking.

Methods

We evaluated the interaction between smoking exposure, considered as either ever vs never or pack-years, and a 26-single nucleotide polymorphisms (SNPs) genetic risk score in relation to FEV 1 or FEV 1 /FVC in 50 047 participants of European ancestry from the Cohorts for Heart and Aging Research in Genomic Epidemiology (CHARGE) and SpiroMeta consortia.

Results

We identified an interaction ( βint  = -0.036, 95% confidence interval, -0.040 to -0.032, P  = 0.00057) between an unweighted 26 SNP genetic risk score and smoking status (ever/never) on the FEV 1 /FVC ratio. In interpreting this interaction, we showed that the genetic risk of falling below the FEV /FVC threshold used to diagnose chronic obstructive pulmonary disease is higher among ever smokers than among never smokers. A replication analysis in two independent datasets, although not statistically significant, showed a similar trend in the interaction effect.

Conclusions

This study highlights the benefit of using genetic risk scores for identifying interactions missed when studying individual SNPs and shows, for the first time, that persons with the highest genetic risk for low FEV 1 /FVC may be more susceptible to the deleterious effects of smoking.

Exome Sequencing Analysis in Severe, Early-Onset Chronic Obstructive Pulmonary Disease

RATIONALE

Genomic regions identified by genome-wide association studies explain only a small fraction of heritability for chronic obstructive pulmonary disease (COPD). Alpha-1 antitrypsin deficiency shows that rare coding variants of large effect also influence COPD susceptibility. We hypothesized that exome sequencing in families identified through a proband with severe, early-onset COPD would identify additional rare genetic determinants of large effect.

OBJECTIVES

To identify rare genetic determinants of severe COPD.

METHODS

We applied filtering approaches to identify potential causal variants for COPD in whole exomes from 347 subjects in 49 extended pedigrees from the Boston Early-Onset COPD Study. We assessed the power of this approach under different levels of genetic heterogeneity using simulations. We tested genes identified in these families using gene-based association tests in exomes of 204 cases with severe COPD and 195 resistant smokers from the COPDGene study. In addition, we examined previously described loci associated with COPD using these datasets.

MEASUREMENTS AND MAIN RESULTS

We identified 69 genes with predicted deleterious nonsynonymous, stop, or splice variants that segregated with severe COPD in at least two pedigrees. Four genes (DNAH8, ALCAM, RARS, and GBF1) also demonstrated an increase in rare nonsynonymous, stop, and/or splice mutations in cases compared with resistant smokers from the COPDGene study; however, these results were not statistically significant. We demonstrate the limitations of the power of this approach under genetic heterogeneity through simulation.

CONCLUSIONS

Rare deleterious coding variants may increase risk for COPD, but multiple genes likely contribute to COPD susceptibility.

Genetic variants affecting cross-sectional lung function in adults show little or no effect on longitudinal lung function decline

BACKGROUND

Genome-wide association studies have identified numerous genetic regions that influence cross-sectional lung function. Longitudinal decline in lung function also includes a heritable component but the genetic determinants have yet to be defined.

OBJECTIVES

We aimed to determine whether regions associated with cross-sectional lung function were also associated with longitudinal decline and to seek novel variants which influence decline.

METHODS

We analysed genome-wide data from 4167 individuals from the Busselton Health Study cohort, who had undergone spirometry (12 695 observations across eight time points). A mixed model was fitted and weighted risk scores were calculated for the joint effect of 26 known regions on baseline and longitudinal changes in FEV₁ and FEV₁/FVC. Potential additional regions of interest were identified and followed up in two independent cohorts.

RESULTS

The 26 regions previously associated with cross-sectional lung function jointly showed a strong effect on baseline lung function (p=4.44×10^-16 for FEV₁/FVC) but no effect on longitudinal decline (p=0.160 for FEV₁/FVC). This was replicated in an independent cohort. 39 additional regions of interest (48 variants) were identified; these associations were not replicated in two further cohorts.

CONCLUSIONS

Previously identified genetic variants jointly have a strong effect on cross-sectional lung function in adults but little or no effect on the rate of decline of lung function. It is possible that they influence COPD risk through lung development. Although no genetic variants have yet been associated with lung function decline at stringent genome-wide significance, longitudinal change in lung function is heritable suggesting that there is scope for future discoveries.

Missing heritability: is the gap closing? An analysis of 32 complex traits in the Lifelines Cohort Study

Despite the recent explosive rise in number of genetic markers for complex disease traits identified in genome-wide association studies, there is still a large gap between the known heritability of these traits and the part explained by these markers. To gauge whether this 'heritability gap' is closing, we first identified genome-wide significant SNPs from the literature and performed replication analyses for 32 highly relevant traits from five broad disease areas in 13 436 subjects of the Lifelines Cohort. Next, we calculated the variance explained by multi-SNP genetic risk scores (GRSs) for each trait, and compared it to their broad- and narrow-sense heritabilities captured by all common SNPs. The majority of all previously-associated SNPs (median=75%) were significantly associated with their respective traits. All GRSs were significant, with unweighted GRSs generally explaining less phenotypic variance than weighted GRSs, for which the explained variance was highest for height (15.5%) and varied between 0.02 and 6.7% for the other traits. Broad-sense common-SNP heritability estimates were significant for all traits, with the additive effect of common SNPs explaining 48.9% of the variance for height and between 5.6 and 39.2% for the other traits. Dominance effects were uniformly small (0-1.5%) and not significant. On average, the variance explained by the weighted GRSs accounted for only 10.7% of the common-SNP heritability of the 32 traits. These results indicate that GRSs may not yet be ready for accurate personalized prediction of complex disease traits limiting widespread adoption in clinical practice.

Genetic pleiotropy between age-related macular degeneration and 16 complex diseases and traits

BACKGROUND

Age-related macular degeneration (AMD) is a common condition of vision loss with disease development strongly influenced by environmental and genetic factors. Recently, 34 loci were associated with AMD at genome-wide significance. So far, little is known about a genetic overlap between AMD and other complex diseases or disease-relevant traits.

METHODS

For each of 60 complex diseases/traits with publicly available genome-wide significant association data, the lead genetic variant per independent locus was extracted and a genetic score was calculated for each disease/trait as the weighted sum of risk alleles. The association with AMD was estimated based on 16,144 AMD cases and 17,832 controls using logistic regression.

RESULTS

Of the respective disease/trait variance, the 60 genetic scores explained on average 4.8% (0.27-20.69%) and 16 of them were found to be significantly associated with AMD (Q-values < 0.01, p values from < 1.0 × 10^-16 to 1.9 × 10^-3). Notably, an increased risk for AMD was associated with reduced risk for cardiovascular diseases, increased risk for autoimmune diseases, higher HDL and lower LDL levels in serum, lower bone-mineral density as well as an increased risk for skin cancer. By restricting the analysis to 1824 variants initially used to compute the 60 genetic scores, we identified 28 novel AMD risk variants (Q-values < 0.01, p values from 1.1 × 10^-7 to 3.0 × 10^-4), known to be involved in cardiovascular disorders, lipid metabolism, autoimmune diseases, anthropomorphic traits, ocular disorders, and neurological diseases. The latter variants represent 20 novel AMD-associated, pleiotropic loci. Genes in the novel loci reinforce previous findings strongly implicating the complement system in AMD pathogenesis.

CONCLUSIONS

We demonstrate a substantial overlap of the genetics of several complex diseases/traits with AMD and provide statistically significant evidence for an additional 20 loci associated with AMD. This highlights the possibility that so far unrelated pathologies may have disease pathways in common.

Genome-wide association analyses for lung function and chronic obstructive pulmonary disease identify new loci and potential druggable targets

Chronic obstructive pulmonary disease (COPD) is characterized by reduced lung function and is the third leading cause of death globally. Through genome-wide association discovery in 48,943 individuals, selected from extremes of the lung function distribution in UK Biobank, and follow-up in 95,375 individuals, we increased the yield of independent signals for lung function from 54 to 97. A genetic risk score was associated with COPD susceptibility (odds ratio per 1 s.d. of the risk score (∼6 alleles) (95% confidence interval) = 1.24 (1.20-1.27), P = 5.05 × 10-49), and we observed a 3.7-fold difference in COPD risk between individuals in the highest and lowest genetic risk score deciles in UK Biobank. The 97 signals show enrichment in genes for development, elastic fibers and epigenetic regulation pathways. We highlight targets for drugs and compounds in development for COPD and asthma (genes in the inositol phosphate metabolism pathway and CHRM3) and describe targets for potential drug repositioning from other clinical indications.