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

Genetic loci associated with chronic obstructive pulmonary disease overlap with loci for lung function and pulmonary fibrosis

Chronic obstructive pulmonary disease (COPD) is a leading cause of mortality worldwide. We performed a genetic association study in 15,256 cases and 47,936 controls, with replication of select top results (P < 5 × 10-6) in 9,498 cases and 9,748 controls. In the combined meta-analysis, we identified 22 loci associated at genome-wide significance, including 13 new associations with COPD. Nine of these 13 loci have been associated with lung function in general population samples, while 4 (EEFSEC, DSP, MTCL1, and SFTPD) are new. We noted two loci shared with pulmonary fibrosis (FAM13A and DSP) but that had opposite risk alleles for COPD. None of our loci overlapped with genome-wide associations for asthma, although one locus has been implicated in joint susceptibility to asthma and obesity. We also identified genetic correlation between COPD and asthma. Our findings highlight new loci associated with COPD, demonstrate the importance of specific loci associated with lung function to COPD, and identify potential regions of genetic overlap between COPD and other respiratory diseases.

Vagal innervation is required for pulmonary function phenotype in Htr4-/- mice

Human genome-wide association studies have identified over 50 loci associated with pulmonary function and related phenotypes, yet follow-up studies to determine causal genes or variants are rare. Single nucleotide polymorphisms in serotonin receptor 4 (HTR4) are associated with human pulmonary function in genome-wide association studies and follow-up animal work has demonstrated that Htr4 is causally associated with pulmonary function in mice, although the precise mechanisms were not identified. We sought to elucidate the role of neural innervation and pulmonary architecture in the lung phenotype of Htr4^-/- animals. We report here that the Htr4^-/- phenotype in mouse is dependent on vagal innervation to the lung. Both ex vivo tracheal ring reactivity and in vivo flexiVent pulmonary functional analyses demonstrate that vagotomy abrogates the Htr4^-/- airway hyperresponsiveness phenotype. Hyperpolarized ³He gas magnetic resonance imaging and stereological assessment of wild-type and Htr4^-/- mice reveal no observable differences in lung volume, inflation characteristics, or pulmonary microarchitecture. Finally, control of breathing experiments reveal substantive differences in baseline breathing characteristics between mice with/without functional HTR4 in breathing frequency, relaxation time, flow rate, minute volume, time of inspiration and expiration and breathing pauses. These results suggest that HTR4's role in pulmonary function likely relates to neural innervation and control of breathing.

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.

Combining omics data to identify genes associated with allergic rhinitis

Allergic rhinitis is a common chronic disorder characterized by immunoglobulin E-mediated inflammation. To identify new genes associated with this trait, we performed genome- and epigenome-wide association studies and linked marginally significant CpGs located in genes or its promoter and SNPs located 1 Mb from the CpGs, by identifying cis methylation quantitative trait loci (mQTL). This approach relies on functional cellular aspects rather than stringent statistical correction. We were able to identify one gene with significant cis-mQTL for allergic rhinitis, caudal-type homeobox 1 (CDX1). We also identified 11 genes with marginally significant cis-mQTLs (p < 0.05) including one with both allergic rhinitis with or without asthma (RNF39). Moreover, most SNPs identified were not located closest to the gene they were linked to through cis-mQTLs counting the one linked to CDX1 located in a gene previously associated with asthma and atopic dermatitis. By combining omics data, we were able to identify new genes associated with allergic rhinitis and better assess the genes linked to associated SNPs.

Association of nsv823469 copy number loss with decreased risk of chronic obstructive pulmonary disease and pulmonary function in Chinese

It is highly possible that copy number variations (CNVs) in susceptible regions have effects on chronic obstructive pulmonary disease (COPD) development, while long noncoding RNA (lncRNAs) have been shown to cause COPD. We hypothesized that the common CNV, named nsv823469 located on 6p22.1, and covering lncRNAs (major histocompatibility complex, class I, A (HLA-A) and HLA complex group 4B (HCG4B)) has an effect on COPD risk. This association was assessed through a two-stage case-control study, and was further confirmed with COPD and pulmonary function-based family analyses, respectively. The copy number loss (0-copy/1-copy) of nsv823469 significantly decreased risk of COPD compared with normal (2-copy) (OR = 0.77, 95% CI = 0.69–0.85). The loss allele, inducing copy number loss of nsv823469, has a tendency to transmit to offspring or siblings (P = 0.010) and is associated with forced expiratory volume in 1 second (FEV1) (P = 0.030). Furthermore, the copy number loss of nsv823469 in normal pulmonary tissue decreases the expression levels of HCG4B (r = 0.315, P = 0.031) and HLA-A (r = 0.296, P = 0.044). Our data demonstrates that nsv823469 plays a role in COPD and pulmonary function inheritance by potentially altering expression of HCG4B.

Whole genome prediction and heritability of childhood asthma phenotypes

Introduction

While whole genome prediction (WGP) methods have recently demonstrated successes in the prediction of complex genetic diseases, they have not yet been applied to asthma and related phenotypes. Longitudinal patterns of lung function differ between asthmatics, but these phenotypes have not been assessed for heritability or predictive ability. Herein, we assess the heritability and genetic predictability of asthma‐related phenotypes.

Methods

We applied several WGP methods to a well‐phenotyped cohort of 832 children with mild‐to‐moderate asthma from CAMP. We assessed narrow‐sense heritability and predictability for airway hyperresponsiveness, serum immunoglobulin E, blood eosinophil count, pre‐ and post‐bronchodilator forced expiratory volume in 1 sec (FEV₁), bronchodilator response, steroid responsiveness, and longitudinal patterns of lung function (normal growth, reduced growth, early decline, and their combinations). Prediction accuracy was evaluated using a training/testing set split of the cohort.

Results

We found that longitudinal lung function phenotypes demonstrated significant narrow‐sense heritability (reduced growth, 95%; normal growth with early decline, 55%). These same phenotypes also showed significant polygenic prediction (areas under the curve [AUCs] 56% to 62%). Including additional demographic covariates in the models increased prediction 4–8%, with reduced growth increasing from 62% to 66% AUC. We found that prediction with a genomic relatedness matrix was improved by filtering available SNPs based on chromatin evidence, and this result extended across cohorts.

Conclusions

Longitudinal reduced lung function growth displayed extremely high heritability. All phenotypes with significant heritability showed significant polygenic prediction. Using SNP‐prioritization increased prediction across cohorts. WGP methods show promise in predicting asthma‐related heritable traits.

Targeted high-throughput sequencing of candidate genes for chronic obstructive pulmonary disease

Background

Reduced lung function in patients with chronic obstructive pulmonary disease (COPD) is likely due to both environmental and genetic factors. We report here a targeted high-throughput DNA sequencing approach to identify new and previously known genetic variants in a set of candidate genes for COPD.

Methods

Exons in 22 genes implicated in lung development as well as 61 genes and 10 genomic regions previously associated with COPD were sequenced using individual DNA samples from 68 cases with moderate or severe COPD and 66 controls matched for age, gender and smoking. Cases and controls were selected from the Obstructive Lung Disease in Northern Sweden (OLIN) studies.

Results

In total, 37 genetic variants showed association with COPD (p < 0.05, uncorrected). Several variants previously discovered to be associated with COPD from genetic genome-wide analysis studies were replicated using our sample. Two high-risk variants were followed-up for functional characterization in a large eQTL mapping study of 1,111 human lung specimens. The C allele of a synonymous variant, rs8040868, predicting a p.(S45=) in the gene for cholinergic receptor nicotinic alpha 3 (CHRNA3) was associated with COPD (p = 8.8 x 10⁻³). This association remained (p = 0.003 and OR = 1.4, 95 % CI 1.1-1.7) when analysing all available cases and controls in OLIN (n = 1,534). The rs8040868 variant is in linkage disequilibrium with rs16969968 previously associated with COPD and altered expression of the CHRNA5 gene. A follow-up analysis for detection of expression quantitative trait loci revealed that rs8040868-C was found to be significantly associated with a decreased expression of the nearby gene cholinergic receptor, nicotinic, alpha 5 (CHRNA5) in lung tissue.

Conclusion

Our data replicate previous result suggesting CHRNA5 as a candidate gene for COPD and rs8040868 as a risk variant for the development of COPD in the Swedish population.

Electronic supplementary material

The online version of this article (doi:10.1186/s12890-016-0309-y) contains supplementary material, which is available to authorized users.

Differential DNA methylation marks and gene comethylation of COPD in African-Americans with COPD exacerbations

BACKGROUND

Chronic obstructive pulmonary disease (COPD) is the third-leading cause of death worldwide. Identifying COPD-associated DNA methylation marks in African-Americans may contribute to our understanding of racial disparities in COPD susceptibility. We determined differentially methylated genes and co-methylation network modules associated with COPD in African-Americans recruited during exacerbations of COPD and smoking controls from the Pennsylvania Study of Chronic Obstructive Pulmonary Exacerbations (PA-SCOPE) cohort.

METHODS

We assessed DNA methylation from whole blood samples in 362 African-American smokers in the PA-SCOPE cohort using the Illumina Infinium HumanMethylation27 BeadChip Array. Final analysis included 19302 CpG probes annotated to the nearest gene transcript after quality control. We tested methylation associations with COPD case-control status using mixed linear models. Weighted gene comethylation networks were constructed using weighted gene coexpression network analysis (WGCNA) and network modules were analyzed for association with COPD.

RESULTS

There were five differentially methylated CpG probes significantly associated with COPD among African-Americans at an FDR less than 5 %, and seven additional probes that approached significance at an FDR less than 10 %. The top ranked gene association was MAML1, which has been shown to affect NOTCH-dependent angiogenesis in murine lung. Network modeling yielded the "yellow" and "blue" comethylation modules which were significantly associated with COPD (p-value 4 × 10-10 and 4 × 10-9, respectively). The yellow module was enriched for gene sets related to inflammatory pathways known to be relevant to COPD. The blue module contained the top ranked genes in the concurrent differential methylation analysis (FXYD1/LGI4, gene significance p-value 1.2 × 10-26; MAML1, p-value 2.0 × 10-26; CD72, p-value 2.1 × 10-25; and LPO, p-value 7.2 × 10-25), and was significantly associated with lung development processes in Gene Ontology gene-set enrichment analysis.

CONCLUSION

We identified 12 differentially methylated CpG sites associated with COPD that mapped to biologically plausible genes. Network module comethylation patterns have identified candidate genes that may be contributing to racial differences in COPD susceptibility and severity. COPD-associated comethylation modules contained genes previously associated with lung disease and inflammation and recapitulated known COPD-associated genes. The genes implicated by differential methylation and WGCNA analysis may provide mechanistic targets contributing to COPD susceptibility, exacerbations, and outcomes among African-Americans.

TRIAL REGISTRATION

Trial Registration: NCT00774176 , Registry: ClinicalTrials.gov, URL: www.clinicaltrials.gov , Date of Enrollment of First Participant: June 2004, Date Registered: 04 January 2008 (retrospectively registered).

Atlas of human diseases influenced by genetic variants with extreme allele frequency differences

Genetic variants with extreme allele frequency differences (EAFD) may underlie some human health disparities across populations. To identify EAFD loci, we systematically analyzed and characterized 81 million genomic variants from 2504 unrelated individuals of 26 world populations (phase III of the 1000 Genomes Project). Our analyses revealed a total of 434 genes, 15 pathways, and 18 diseases and traits influenced by EAFD variants from five continental populations. They included known EAFD genes, such as LCT (lactose tolerance), SLC24A5 (skin pigmentation), and EDAR (hair morphology). We found many novel EAFD genes, including TBC1D2B (autophagy mediator), TRIM40 (gastrointestinal inflammatory regulator), KRT71, KRT75, KRT83, and KRTAP10-1 (hair and epithelial keratin synthesis), PIK3R3 (insulin receptor interaction), DARS (neurological disorders), and NACA2 (skin inflammatory response). Our results also showed four complex diseases significantly associated with EAFD loci, including asthma (adjusted enrichment P = 4 × 10^-8), type I diabetes (P = 6 × 10^-9), alcohol consumption (P = 0.0002), and attention deficit/hyperactivity disorder (P = 0.003). This study provides a comprehensive atlas of genes, pathways, and human diseases significantly influenced by EAFD variants.

DNA methylation profiling in human lung tissue identifies genes associated with COPD

Chronic obstructive pulmonary disease (COPD) is a smoking-related disease characterized by genetic and phenotypic heterogeneity. Although association studies have identified multiple genomic regions with replicated associations to COPD, genetic variation only partially explains the susceptibility to lung disease, and suggests the relevance of epigenetic investigations. We performed genome-wide DNA methylation profiling in homogenized lung tissue samples from 46 control subjects with normal lung function and 114 subjects with COPD, all former smokers. The differentially methylated loci were integrated with previous genome-wide association study results. The top 535 differentially methylated sites, filtered for a minimum mean methylation difference of 5% between cases and controls, were enriched for CpG shelves and shores. Pathway analysis revealed enrichment for transcription factors. The top differentially methylated sites from the intersection with previous GWAS were in CHRM1, GLT1D1, and C10orf11; sorted by GWAS P-value, the top sites included FRMD4A, THSD4, and C10orf11. Epigenetic association studies complement genetic association studies to identify genes potentially involved in COPD pathogenesis. Enrichment for genes implicated in asthma and lung function and for transcription factors suggests the potential pathogenic relevance of genes identified through differential methylation and the intersection with a broader range of GWAS associations.

Genome-wide association study on the FEV1/FVC ratio in never-smokers identifies HHIP and FAM13A

BACKGROUND

Although a striking proportion (25% to 45%) of patients with chronic obstructive pulmonary disease are never-smokers, most genetic susceptibility studies have not focused on this group exclusively.

OBJECTIVE

The aim of this study was to identify common genetic variants associated with FEV₁ and its ratio to forced vital capacity (FVC) in never-smokers.

METHODS

Genome-wide association studies were performed in 5070 never-smokers of the identification cohort LifeLines, and results (P < 10^-5) were verified by using a meta-analysis of the Vlagtwedde-Vlaardingen study and the Rotterdam Study I-III (total n = 1966). Furthermore, we aimed to assess the effects of the replicated variants in more detail by performing genetic risk score, expression quantitative trait loci, and variant*ever-smoking interaction analyses.

RESULTS

We identified associations between the FEV₁/FVC ratio and 5 common genetic variants in the identification cohort, and 2 of these associations were replicated. The 2 variants annotated to the genes hedgehog interacting protein (HHIP) and family with sequence similarity 13 member A (FAM13A) were shown to have an additive effect on FEV₁/FVC levels in the genetic risk score analysis; were associated with gene expression of HHIP and FAM13A in lung tissue, respectively; and were genome-wide significant in a meta-analysis including both identification and 4 verification cohorts (P < 2.19 × 10^-7). Finally, we did not identify significant interactions between the variants and ever smoking. Results of the FEV₁ identification analysis were not replicated.

CONCLUSION

The genes HHIP and FAM13A confer a risk for airway obstruction in general that is not driven exclusively by cigarette smoking, which is the main risk factor for chronic obstructive pulmonary disease.

Machine learning, statistical learning and the future of biological research in psychiatry

Psychiatric research has entered the age of 'Big Data'. Datasets now routinely involve thousands of heterogeneous variables, including clinical, neuroimaging, genomic, proteomic, transcriptomic and other 'omic' measures. The analysis of these datasets is challenging, especially when the number of measurements exceeds the number of individuals, and may be further complicated by missing data for some subjects and variables that are highly correlated. Statistical learning-based models are a natural extension of classical statistical approaches but provide more effective methods to analyse very large datasets. In addition, the predictive capability of such models promises to be useful in developing decision support systems. That is, methods that can be introduced to clinical settings and guide, for example, diagnosis classification or personalized treatment. In this review, we aim to outline the potential benefits of statistical learning methods in clinical research. We first introduce the concept of Big Data in different environments. We then describe how modern statistical learning models can be used in practice on Big Datasets to extract relevant information. Finally, we discuss the strengths of using statistical learning in psychiatric studies, from both research and practical clinical points of view.

Genome-wide association study of copy number variation with lung function identifies a novel signal of association near BANP for forced vital capacity

BACKGROUND

Genome-wide association studies of Single Nucleotide Polymorphisms (SNPs) have identified 55 SNPs associated with lung function. However, little is known about the effect of copy number variants (CNVs) on lung function, although CNVs represent a significant proportion of human genetic polymorphism. To assess the effect of CNVs on lung function quantitative traits, we measured copy number at 2788 previously characterised, common copy number variable regions in 6 independent cohorts (n = 24,237) using intensity data from SNP genotyping experiments. We developed a pipeline for genome-wide association analysis and meta-analysis of CNV genotypes measured across multiple studies using SNP genotype array intensity data from different platform technologies. We then undertook cohort-level genome-wide association studies of CNV with lung function in a subset of 4 cohorts (n < =12,403) with lung function measurements and meta-analysed the results. Follow-up was undertaken for CNVs which were well tagged by SNPs, in up to 146,871 individuals.

RESULTS

We generated robust copy number calls for 1962 out of 2788 (70 %) known CNV regions genome-wide, with 1103 measured with compatible class frequencies in at least 2 cohorts. We report a novel CNV association (discovery P = 0.0007) with Forced Vital Capacity (FVC) downstream of BANP on chromosome 16 that shows evidence of replication by a tag SNP in two independent studies (replication P = 0.004). In addition, we provide suggestive evidence (discovery P = 0.0002) for a role of complex copy number variation at a previously reported lung function locus, containing the rootletin gene CROCC, that is not tagged by SNPs.

CONCLUSIONS

We demonstrate how common CNV regions can be reliably and consistently called across cohorts, using an existing calling algorithm and rigorous quality control steps, using SNP genotyping array intensity data. Although many common biallelic CNV regions were well-tagged by common SNPs, we also identified associations with untagged mulitallelic CNV regions thereby illustrating the potential of our approach to identify some of the missing heritability of complex traits.

Hhip haploinsufficiency sensitizes mice to age-related emphysema

Genetic variants in Hedgehog interacting protein (HHIP) have consistently been associated with the susceptibility to develop chronic obstructive pulmonary disease and pulmonary function levels, including the forced expiratory volume in 1 s (FEV1), in general population samples by genome-wide association studies. However, in vivo evidence connecting Hhip to age-related FEV1 decline and emphysema development is lacking. Herein, using Hhip heterozygous mice (Hhip(+/-)), we observed increased lung compliance and spontaneous emphysema in Hhip(+/-) mice starting at 10 mo of age. This increase was preceded by increases in oxidative stress levels in the lungs of Hhip(+/-) vs. Hhip(+/+) mice. To our knowledge, these results provide the first line of evidence that HHIP is involved in maintaining normal lung function and alveolar structures. Interestingly, antioxidant N-acetyl cysteine treatment in mice starting at age of 5 mo improved lung function and prevented emphysema development in Hhip(+/-) mice, suggesting that N-acetyl cysteine treatment limits the progression of age-related emphysema in Hhip(+/-) mice. Therefore, reduced lung function and age-related spontaneous emphysema development in Hhip(+/-) mice may be caused by increased oxidative stress levels in murine lungs as a result of haploinsufficiency of Hhip.

Risk of Lung Disease in PI MZ Heterozygotes. Current Status and Future Research Directions

The potential for increased chronic obstructive pulmonary disease (COPD) risk among PI MZ subjects was initially recognized decades ago. However, despite many studies of this topic, it has remained controversial whether such increased risk exists. Several recent studies in large populations strongly support increased risk for COPD among PI MZ subjects. This increased PI MZ risk will need to be understood in the context of other identified COPD genetic determinants and investigations of COPD phenotypic heterogeneity.

Genes associated with polymorphic variants predicting lung function are differentially expressed during human lung development

BACKGROUND

Recent meta-analyses of genome-wide association studies have identified single nucleotide polymorphisms (SNPs) within/near 54 genes associated with lung function measures. Current understanding of the contribution of these genes to human lung development is limited. We set out to further define i) the expression profile of these genes during human lung development using a unique set of resources to examine both mRNA and protein expression and ii) the link between key polymorphisms and genes using expression quantitative trait loci (eQTL) approaches.

METHODS

The mRNA expression profile of lung function associated genes across pseudoglandular and canalicular stages of lung development were determined using expression array data of 38 human fetal lungs. eQTLs were investigated for selected genes using blood cell and lung tissue data. Immunohistochemistry of the top 5 candidates was performed in a panel of 24 fetal lung samples.

RESULTS

Twenty-nine lung function associated genes were differentially expressed during lung development at the mRNA level. The greatest magnitude of effect was observed for 5 genes; TMEM163, FAM13A and HHIP which had increasing expression and CDC123 and PTCH1 with decreased expression across developmental stages. Focussed eQTL analyses investigating SNPs in these five loci identified several cis-eQTL's. Protein expression of TMEM163 increased and CDC123 decreased with fetal lung age in agreement with mRNA data. Protein expression in FAM13A, HHIP and PTCH1 remained relatively constant throughout lung development.

CONCLUSIONS

We have identified that > 50 % of lung function associated genes show evidence of differential expression during lung development and we show that in particular TMEM163 and CDC123 are differentially expressed at both the mRNA and protein levels. Our data provides a systematic evaluation of lung function associated genes in this context and offers some insight into the potential role of several of these genes in contributing to human lung development.

Exome Array Analysis Identifies a Common Variant in IL27 Associated with Chronic Obstructive Pulmonary Disease

RATIONALE

Chronic obstructive pulmonary disease (COPD) susceptibility is in part related to genetic variants. Most genetic studies have been focused on genome-wide common variants without a specific focus on coding variants, but common and rare coding variants may also affect COPD susceptibility.

OBJECTIVES

To identify coding variants associated with COPD.

METHODS

We tested nonsynonymous, splice, and stop variants derived from the Illumina HumanExome array for association with COPD in five study populations enriched for COPD. We evaluated single variants with a minor allele frequency greater than 0.5% using logistic regression. Results were combined using a fixed effects meta-analysis. We replicated novel single-variant associations in three additional COPD cohorts.

MEASUREMENTS AND MAIN RESULTS

We included 6,004 control subjects and 6,161 COPD cases across five cohorts for analysis. Our top result was rs16969968 (P = 1.7 × 10(-14)) in CHRNA5, a locus previously associated with COPD susceptibility and nicotine dependence. Additional top results were found in AGER, MMP3, and SERPINA1. A nonsynonymous variant, rs181206, in IL27 (P = 4.7 × 10(-6)) was just below the level of exome-wide significance but attained exome-wide significance (P = 5.7 × 10(-8)) when combined with results from other cohorts. Gene expression datasets revealed an association of rs181206 and the surrounding locus with expression of multiple genes; several were differentially expressed in COPD lung tissue, including TUFM.

CONCLUSIONS

In an exome array analysis of COPD, we identified nonsynonymous variants at previously described loci and a novel exome-wide significant variant in IL27. This variant is at a locus previously described in genome-wide associations with diabetes, inflammatory bowel disease, and obesity and appears to affect genes potentially related to COPD pathogenesis.

Meta-analysis of 375,000 individuals identifies 38 susceptibility loci for migraine

Migraine is a debilitating neurological disorder affecting around one in seven people worldwide, but its molecular mechanisms remain poorly understood. There is some debate about whether migraine is a disease of vascular dysfunction or a result of neuronal dysfunction with secondary vascular changes. Genome-wide association (GWA) studies have thus far identified 13 independent loci associated with migraine. To identify new susceptibility loci, we carried out a genetic study of migraine on 59,674 affected subjects and 316,078 controls from 22 GWA studies. We identified 44 independent single-nucleotide polymorphisms (SNPs) significantly associated with migraine risk (P < 5 × 10(-8)) that mapped to 38 distinct genomic loci, including 28 loci not previously reported and a locus that to our knowledge is the first to be identified on chromosome X. In subsequent computational analyses, the identified loci showed enrichment for genes expressed in vascular and smooth muscle tissues, consistent with a predominant theory of migraine that highlights vascular etiologies.

Post-GWAS Prioritization Through Data Integration Provides Novel Insights on Chronic Obstructive Pulmonary Disease

Rich collections of genomic and epigenomic annotations, availabilities of large population cohorts for genome-wide association studies (GWAS), and advancements in data integration techniques provide the unprecedented opportunity to accelerate discoveries in complex disease studies through integrative analyses. In this paper, we apply a variety of approaches to integrate GWAS summary statistics of chronic obstructive pulmonary disease (COPD) with functional annotations to illustrate how data integration could help researchers understand complex human diseases. We show that incorporating functional annotations can better prioritize GWAS signals at both the global and the local levels. Signal prioritization on severe COPD GWAS reveals multiple potential risk loci that are linked with pulmonary functions. Enrichment analysis provides novel insights on the pathogenesis of COPD and hints the existence of genetic contributions to muscle dysfuncion and chronic lung inflammation, two symptoms that are often co-morbid with COPD. Our results suggest that rich signals for COPD genetics are still buried under the Bonferroni-corrected genome-wide significance threshold. Many more biological findings are expected to emerge as more samples are recruited for COPD studies.

PR Domain-containing Protein 7 (PRDM7) Is a Histone 3 Lysine 4 Trimethyltransferase

PR domain-containing protein 7 (PRDM7) is a primate-specific histone methyltransferase that is the result of a recent gene duplication of PRDM9. The two proteins are highly homologous, especially in the catalytic PR/SET domain, where they differ by only three amino acid residues. Here we report that PRDM7 is an efficient methyltransferase that selectively catalyzes the trimethylation of H3 lysine 4 (H3K4) both in vitro and in cells. Through selective mutagenesis we have dissected the functional roles of each of the three divergent residues between the PR domains of PRDM7 and PRDM9. These studies indicate that after a single serine to tyrosine mutation at residue 357 (S357Y), PRDM7 regains the substrate specificities and catalytic activities similar to its evolutionary predecessor, including the ability to efficiently methylate H3K36.

T cell factor-4 functions as a co-activator to promote NF-κB-dependent MMP-15 expression in lung carcinoma cells

Both TCF-4 and MMP-15 are closely linked to the development of lung cancer, while the regulatory role of TCF-4 in MMP-15 expression is still obscure. Here we found that expression of TCF-4 and MMP-15 was increased in lung cancer cells or tissues versus the normal ones. With gain-or loss-of -function studies, we demonstrated that TCF-4 positively regulated MMP-15 expression in mRNA and protein levels. With reporter gene assay, we found that TCF-4 regulated MMP-15 expression via a potential NF-κB binding element locating at -2833/-2824 in the mouse MMP-15 promoter. With ChIP and immunoblotting assays, we identified that TCF-4 functioned as a co-activator to potentiate the binding between p65 and MMP-15 promoter. Functionally, TCF-4 silence attenuated the migration activity of LLC cells, while additional overexpression of MMP-15 rescued this effect in cell scratch test and transwell migration assay. In xenograft model, TCF-4 silence-improved tumor lesions in lungs and survival time of LLC-tumor bearing mice were abolished by MMP-15 overexpression. In conclusion, we are the first to identify TCF-4 as a co-activator of NF-κB p65 to promote MMP-15 transcription and potentiate the migration activity of the lung cancer cells. Our findings shed light on the therapeutic strategies of this malignancy.

Genetic variants of CDH13 determine the susceptibility to chronic obstructive pulmonary disease in a Chinese population

AIM

Adiponectin has been implicated in the development of chronic obstructive pulmonary disease (COPD). The CDH13 gene encodes T-cadherin that is an adiponectin receptor, and genetic variants of CDH13 determine blood adiponectin levels. The aim of this study was to investigate the effects of CDH13 variants on COPD susceptibility in a Chinese population.

METHODS

Ten single-nucleotide polymorphisms (SNP) in CDH13 were screened using the SNaPshot method in 279 COPD patients and 367 control subjects. Association of genotypes or haplotypes constructed from these loci with COPD was analyzed in different genetic models.

RESULTS

Among the 10 SNPs tested, rs4783244 and rs12922394 exhibited significant differences in allele or genotype frequencies between COPD patients and control subjects, whereas 8 other SNPs did not. The minor allele T was associated with decreased risk of COPD in the recessive model at rs4783244 (OR=0.42, P=0.023) and in the dominant model at rs12922394 (OR=0.70, P=0.022). The genotype TT at either rs4783244 or rs12922394 was associated with a significantly low level of plasma adiponectin when compared to genotypes GG and CC (P<0.05). Haplotypes GC in block 1 (rs4783244-rs12922394) as well as GTAC and ATGT in block 3 (rs4783266-rs11640522-rs11646849-rs11860282) significantly increased the risk of COPD, whereas haplotypes TT in block 1, TG in block 2 (rs11646011- rs11640875) and ATGC in block 3 were protective against COPD.

CONCLUSION

CDH13 genetic variants determine Chinese individuals' susceptibility to COPD and thus are efficient genetic biomarkers for early detection of COPD.

Exome-wide analysis of rare coding variation identifies novel associations with COPD and airflow limitation in MOCS3, IFIT3 and SERPINA12

Background

Several regions of the genome have shown to be associated with COPD in genome-wide association studies of common variants.

Objective

To determine rare and potentially functional single nucleotide polymorphisms (SNPs) associated with the risk of COPD and severity of airflow limitation.

Methods

3226 current or former smokers of European ancestry with lung function measures indicative of Global Initiative for Chronic Obstructive Lung Disease (GOLD) 2 COPD or worse were genotyped using an exome array. An analysis of risk of COPD was carried out using ever smoking controls (n=4784). Associations with %predicted FEV₁ were tested in cases. We followed-up signals of interest (p<10⁻⁵) in independent samples from a subset of the UK Biobank population and also undertook a more powerful discovery study by meta-analysing the exome array data and UK Biobank data for variants represented on both arrays.

Results

Among the associated variants were two in regions previously unreported for COPD; a low frequency non-synonymous SNP in MOCS3 (rs7269297, p_discovery=3.08×10⁻⁶, p_replication=0.019) and a rare SNP in IFIT3, which emerged in the meta-analysis (rs140549288, p_meta=8.56×10⁻⁶). In the meta-analysis of % predicted FEV₁ in cases, the strongest association was shown for a splice variant in a previously unreported region, SERPINA12 (rs140198372, p_meta=5.72×10⁻⁶). We also confirmed previously reported associations with COPD risk at MMP12, HHIP, GPR126 and CHRNA5. No associations in novel regions reached a stringent exome-wide significance threshold (p<3.7×10⁻⁷).

Conclusions

This study identified several associations with the risk of COPD and severity of airflow limitation, including novel regions MOCS3, IFIT3 and SERPINA12, which warrant further study.

The role of WWOX polymorphisms on COPD susceptibility and pulmonary function traits in Chinese: a case-control study and family-based analysis

Single nucleotide polymorphisms (SNPs) in the WW domain containing oxidoreductase (WWOX) gene were recently identified to be quantitative trait loci for lung function and thus likely to be susceptible biomarkers for COPD. However, the associations between WWOX SNPs and COPD risk are still unclear. Here, by conducting a two-center case-control study including 1511 COPD cases and 1677 controls and a family-based analysis comprising 95 nuclear pedigrees, we tested the associations between five SNPs that are rs10220974C >T, rs3764340C >G, rs12918952G >A, rs383362G >T, rs12828G >A of WWOX and COPD risk as well as the hereditary inclination of these loci among COPD families. We found that the SNP rs383362G >T was significantly associated with an increased risk of COPD in a T allele-number dependent-manner (OR = 1.30, 95%CI = 1.11 - 1.52). The T allele was more prone to over transmit to sick children and sibs than the G allele (Z = 2.900, P = 0.004). Moreover, the forced expiratory volume in one second/forced vital capacity (FEV1/FVC), FEV1/predicted-FEV1 and annual FEV1 also significantly decreased in the rs383362T carriers compared to the rs383362GG carriers. For other SNPs, no significant association was observed for COPD and pulmonary function. Taken together, our data demonstrated that the SNP rs383362G >T of WWOX plays a role in COPD inheritance.

Ndufc2 Gene Inhibition Is Associated With Mitochondrial Dysfunction and Increased Stroke Susceptibility in an Animal Model of Complex Human Disease

BACKGROUND

The genetic basis of stroke susceptibility remains to be elucidated. STR1 quantitative trait locus (STR1/QTL) was identified on rat chromosome 1 of stroke-prone spontaneously hypertensive rat (SHRSP) upon Japanese-style stroke-permissive diet (JD), and it contributes to 20% of the stroke phenotype variance.

METHODS AND RESULTS

Nine hundred eighty-six probe sets mapping on STR1 were selected from the Rat RAE230A array and screened through a microarray differential expression analysis in brains of SHRSP and stroke-resistant SHR (SHRSR) fed with either regular diet or JD. The gene encoding Ndufc2 (NADH dehydrogenase [ubiquinone] 1 subunit), mapping 8 Mb apart from STR1/QTL Lod score peak, was found significantly down-regulated under JD in SHRSP compared to SHRSR. Ndufc2 disruption altered complex I assembly and activity, reduced mitochondrial membrane potential and ATP levels, and increased reactive oxygen species production and inflammation both in vitro and in vivo. SHRSR carrying heterozygous Ndufc2 deletion showed renal abnormalities and stroke occurrence under JD, similarly to SHRSP. In humans, T allele variant at NDUFC2/rs11237379 was associated with significant reduction in gene expression and with increased occurrence of early-onset ischemic stroke by recessive mode of transmission (odds ratio [OR], 1.39; CI, 1.07-1.80; P=0.012). Subjects carrying TT/rs11237379 and A allele variant at NDUFC2/rs641836 had further increased risk of stroke (OR=1.56; CI, 1.14-2.13; P=0.006).

CONCLUSIONS

A significant reduction of Ndufc2 expression causes complex I dysfunction and contributes to stroke susceptibility in SHRSP. Moreover, our current evidence may suggest that Ndufc2 can contribute to an increased occurrence of early-onset ischemic stroke in humans.

Polymorphisms associated with expression of BPIFA1/BPIFB1 and lung disease severity in cystic fibrosis

BPI fold containing family A, member 1 (BPIFA1) and BPIFB1 are putative innate immune molecules expressed in the upper airways. Because of their hypothesized roles in airway defense, these molecules may contribute to lung disease severity in cystic fibrosis (CF). We interrogated BPIFA1/BPIFB1 single-nucleotide polymorphisms in data from an association study of CF modifier genes and found an association of the G allele of rs1078761 with increased lung disease severity (P = 2.71 × 10(-4)). We hypothesized that the G allele of rs1078761 is associated with decreased expression of BPIFA1 and/or BPIFB1. Genome-wide lung gene expression and genotyping data from 1,111 individuals with lung disease, including 51 patients with CF, were tested for associations between genotype and BPIFA1 and BPIFB1 gene expression levels. Findings were validated by quantitative PCR in a subset of 77 individuals. Western blotting was used to measure BPIFA1 and BPIFB1 protein levels in 93 lung and 101 saliva samples. The G allele of rs1078761 was significantly associated with decreased mRNA levels of BPIFA1 (P = 4.08 × 10(-15)) and BPIFB1 (P = 0.0314). These findings were confirmed with quantitative PCR and Western blotting. We conclude that the G allele of rs1078761 may be detrimental to lung function in CF owing to decreased levels of BPIFA1 and BPIFB1.

Association of Forced Vital Capacity with the Developmental Gene NCOR2

Background

Forced Vital Capacity (FVC) is an important predictor of all-cause mortality in the absence of chronic respiratory conditions. Epidemiological evidence highlights the role of early life factors on adult FVC, pointing to environmental exposures and genes affecting lung development as risk factors for low FVC later in life. Although highly heritable, a small number of genes have been found associated with FVC, and we aimed at identifying further genetic variants by focusing on lung development genes.

Methods

Per-allele effects of 24,728 SNPs in 403 genes involved in lung development were tested in 7,749 adults from three studies (NFBC1966, ECRHS, EGEA). The most significant SNP for the top 25 genes was followed-up in 46,103 adults (CHARGE and SpiroMeta consortia) and 5,062 children (ALSPAC). Associations were considered replicated if the replication p-value survived Bonferroni correction (p<0.002; 0.05/25), with a nominal p-value considered as suggestive evidence. For SNPs with evidence of replication, effects on the expression levels of nearby genes in lung tissue were tested in 1,111 lung samples (Lung eQTL consortium), with further functional investigation performed using public epigenomic profiling data (ENCODE).

Results

NCOR2-rs12708369 showed strong replication in children (p = 0.0002), with replication unavailable in adults due to low imputation quality. This intronic variant is in a strong transcriptional enhancer element in lung fibroblasts, but its eQTL effects could not be tested due to low imputation quality in the eQTL dataset. SERPINE2-rs6754561 replicated at nominal level in both adults (p = 0.036) and children (p = 0.045), while WNT16-rs2707469 replicated at nominal level only in adults (p = 0.026). The eQTL analyses showed association of WNT16-rs2707469 with expression levels of the nearby gene CPED1. We found no statistically significant eQTL effects for SERPINE2-rs6754561.

Conclusions

We have identified a new gene, NCOR2, in the retinoic acid signalling pathway pointing to a role of vitamin A metabolism in the regulation of FVC. Our findings also support SERPINE2, a COPD gene with weak previous evidence of association with FVC, and suggest WNT16 as a further promising candidate.

Complex genetics of pulmonary diseases: lessons from genome-wide association studies and next-generation sequencing

The advent of high-throughput technologies has provided exceptional assistance for lung scientists to discover novel genetic variants underlying the development and progression of complex lung diseases. However, the discovered variants thus far do not explain much of the estimated heritability of complex lung diseases. Here, we review the literature of successfully used genome-wide association studies (GWASs) and identified the polymorphisms that reproducibly underpin the susceptibility to various noncancerous complex lung diseases or affect therapeutic responses. We also discuss the inherent limitations of GWAS approaches and how the use of next-generation sequencing technologies has furthered our understanding about the genetic determinants of these diseases. Next, we describe the contribution of the metagenomics to understand the interactions of the airways microbiome with lung diseases. We then highlight the urgent need for new integrative genomics-phenomics methods to more effectively interrogate and understand multiple downstream "omics" (eg, chromatin modification patterns). Finally, we address the scarcity of genetic studies addressing under-represented populations such as African Americans and Hispanics.

Role of Lung Function Genes in the Development of Asthma

Although our previous GWAS failed to identify SNPs associated with pulmonary function at the level of genomewide significance, it did show that the heritability for FEV1/FVC was 41.6% in a Japanese population, suggesting that the heritability of pulmonary function traits can be explained by the additive effects of multiple common SNPs. In addition, our previous study indicated that pulmonary function genes identified in previous GWASs in non-Japanese populations accounted for 4.3% to 12.0% of the entire estimated heritability of FEV1/FVC in a Japanese population. Therefore, given that many loci with individual weak effects may contribute to asthma risk, in this study, we created a quantitative score of genetic load based on 16 SNPs implicated in lower lung function in both Japanese and non-Japanese populations. This genetic risk score (GRS) for lower FEV1/FVC was consistently associated with the onset of asthma (P = 9.6 × 10(-4)) in 2 independent Japanese populations as well as with the onset of COPD (P = 0.042). Clustering of asthma patients based on GRS levels indicated that an increased GRS may be responsible for the development of a particular phenotype of asthma characterized by early onset, atopy, and severer airflow obstruction.

Sixteen new lung function signals identified through 1000 Genomes Project reference panel imputation

Lung function measures are used in the diagnosis of chronic obstructive pulmonary disease. In 38,199 European ancestry individuals, we studied genome-wide association of forced expiratory volume in 1 s (FEV1), forced vital capacity (FVC) and FEV1/FVC with 1000 Genomes Project (phase 1)-imputed genotypes and followed up top associations in 54,550 Europeans. We identify 14 novel loci (P<5 × 10(-8)) in or near ENSA, RNU5F-1, KCNS3, AK097794, ASTN2, LHX3, CCDC91, TBX3, TRIP11, RIN3, TEKT5, LTBP4, MN1 and AP1S2, and two novel signals at known loci NPNT and GPR126, providing a basis for new understanding of the genetic determinants of these traits and pulmonary diseases in which they are altered.

A genome-wide association study identifies risk loci for spirometric measures among smokers of European and African ancestry

BACKGROUND

Pulmonary function decline is a major contributor to morbidity and mortality among smokers. Post bronchodilator FEV1 and FEV1/FVC ratio are considered the standard assessment of airflow obstruction. We performed a genome-wide association study (GWAS) in 9919 current and former smokers in the COPDGene study (6659 non-Hispanic Whites [NHW] and 3260 African Americans [AA]) to identify associations with spirometric measures (post-bronchodilator FEV1 and FEV1/FVC). We also conducted meta-analysis of FEV1 and FEV1/FVC GWAS in the COPDGene, ECLIPSE, and GenKOLS cohorts (total n = 13,532).

RESULTS

Among NHW in the COPDGene cohort, both measures of pulmonary function were significantly associated with SNPs at the 15q25 locus [containing CHRNA3/5, AGPHD1, IREB2, CHRNB4] (lowest p-value = 2.17 × 10(-11)), and FEV1/FVC was associated with a genomic region on chromosome 4 [upstream of HHIP] (lowest p-value = 5.94 × 10(-10)); both regions have been previously associated with COPD. For the meta-analysis, in addition to confirming associations to the regions near CHRNA3/5 and HHIP, genome-wide significant associations were identified for FEV1 on chromosome 1 [TGFB2] (p-value = 8.99 × 10(-9)), 9 [DBH] (p-value = 9.69 × 10(-9)) and 19 [CYP2A6/7] (p-value = 3.49 × 10(-8)) and for FEV1/FVC on chromosome 1 [TGFB2] (p-value = 8.99 × 10(-9)), 4 [FAM13A] (p-value = 3.88 × 10(-12)), 11 [MMP3/12] (p-value = 3.29 × 10(-10)) and 14 [RIN3] (p-value = 5.64 × 10(-9)).

CONCLUSIONS

In a large genome-wide association study of lung function in smokers, we found genome-wide significant associations at several previously described loci with lung function or COPD. We additionally identified a novel genome-wide significant locus with FEV1 on chromosome 9 [DBH] in a meta-analysis of three study populations.

Chronic obstructive pulmonary disease

Chronic obstructive pulmonary disease (COPD) is a common disease with high global morbidity and mortality. COPD is characterized by poorly reversible airway obstruction, which is confirmed by spirometry, and includes obstruction of the small airways (chronic obstructive bronchiolitis) and emphysema, which lead to air trapping and shortness of breath in response to physical exertion. The most common risk factor for the development of COPD is cigarette smoking, but other environmental factors, such as exposure to indoor air pollutants - especially in developing countries - might influence COPD risk. Not all smokers develop COPD and the reasons for disease susceptibility in these individuals have not been fully elucidated. Although the mechanisms underlying COPD remain poorly understood, the disease is associated with chronic inflammation that is usually corticosteroid resistant. In addition, COPD involves accelerated ageing of the lungs and an abnormal repair mechanism that might be driven by oxidative stress. Acute exacerbations, which are mainly triggered by viral or bacterial infections, are important as they are linked to a poor prognosis. The mainstay of the management of stable disease is the use of inhaled long-acting bronchodilators, whereas corticosteroids are beneficial primarily in patients who have coexisting features of asthma, such as eosinophilic inflammation and more reversibility of airway obstruction. Apart from smoking cessation, no treatments reduce disease progression. More research is needed to better understand disease mechanisms and to develop new treatments that reduce disease activity and progression.

Genome-wide patterns of selection in 230 ancient Eurasians

Ancient DNA makes it possible to directly witness natural selection by analyzing samples from populations before, during and after adaptation events. Here we report the first scan for selection using ancient DNA, capitalizing on the largest genome-wide dataset yet assembled: 230 West Eurasians dating to between 6500 and 1000 BCE, including 163 with newly reported data. The new samples include the first genome-wide data from the Anatolian Neolithic culture whose genetic material we extracted from the DNA-rich petrous bone and who we show were members of the population that was the source of Europe’s first farmers. We also report a complete transect of the steppe region in Samara between 5500 and 1200 BCE that allows us to recognize admixture from at least two external sources into steppe populations during this period. We detect selection at loci associated with diet, pigmentation and immunity, and two independent episodes of selection on height.

Shared Genetic Signals of Hypoxia Adaptation in Drosophila and in High-Altitude Human Populations

The ability to withstand low oxygen (hypoxia tolerance) is a polygenic and mechanistically conserved trait that has important implications for both human health and evolution. However, little is known about the diversity of genetic mechanisms involved in hypoxia adaptation in evolving populations. We used experimental evolution and whole-genome sequencing in Drosophila melanogaster to investigate the role of natural variation in adaptation to hypoxia. Using a generalized linear mixed model we identified significant allele frequency differences between three independently evolved hypoxia-tolerant populations and normoxic control populations for approximately 3,800 single nucleotide polymorphisms. Around 50% of these variants are clustered in 66 distinct genomic regions. These regions contain genes that are differentially expressed between hypoxia-tolerant and normoxic populations and several of the differentially expressed genes are associated with metabolic processes. Additional genes associated with respiratory and open tracheal system development also show evidence of directional selection. RNAi-mediated knockdown of several candidate genes’ expression significantly enhanced survival in severe hypoxia. Using genomewide single nucleotide polymorphism data from four high-altitude human populations—Sherpas, Tibetans, Ethiopians, and Andeans, we found that several human orthologs of the genes under selection in flies are also likely under positive selection in all four high-altitude human populations. Thus, our results indicate that selection for hypoxia tolerance can act on standing genetic variation in similar genes and pathways present in organisms diverged by hundreds of millions of years.

Systematic Integration of Brain eQTL and GWAS Identifies ZNF323 as a Novel Schizophrenia Risk Gene and Suggests Recent Positive Selection Based on Compensatory Advantage on Pulmonary Function

Genome-wide association studies have identified multiple risk variants and loci that show robust association with schizophrenia. Nevertheless, it remains unclear how these variants confer risk to schizophrenia. In addition, the driving force that maintains the schizophrenia risk variants in human gene pool is poorly understood. To investigate whether expression-associated genetic variants contribute to schizophrenia susceptibility, we systematically integrated brain expression quantitative trait loci and genome-wide association data of schizophrenia using Sherlock, a Bayesian statistical framework. Our analyses identified ZNF323 as a schizophrenia risk gene (P = 2.22×10(-6)). Subsequent analyses confirmed the association of the ZNF323 and its expression-associated single nucleotide polymorphism rs1150711 in independent samples (gene-expression: P = 1.40×10(-6); single-marker meta-analysis in the combined discovery and replication sample comprising 44123 individuals: P = 6.85×10(-10)). We found that the ZNF323 was significantly downregulated in hippocampus and frontal cortex of schizophrenia patients (P = .0038 and P = .0233, respectively). Evidence for pleiotropic effects was detected (association of rs1150711 with lung function and gene expression of ZNF323 in lung: P = 6.62×10(-5) and P = 9.00×10(-5), respectively) with the risk allele (T allele) for schizophrenia acting as protective allele for lung function. Subsequent population genetics analyses suggest that the risk allele (T) of rs1150711 might have undergone recent positive selection in human population. Our findings suggest that the ZNF323 is a schizophrenia susceptibility gene whose expression may influence schizophrenia risk. Our study also illustrates a possible mechanism for maintaining schizophrenia risk variants in the human gene pool.

Candidate genes for COPD: current evidence and research

COPD is a common complex disease characterized by progressive airflow limitation. Several genome-wide association studies (GWASs) have discovered genes that are associated with COPD. Recently, candidate genes for COPD identified by GWASs include CHRNA3/5 (cholinergic nicotine receptor alpha 3/5), IREB2 (iron regulatory binding protein 2), HHIP (hedgehog-interacting protein), FAM13A (family with sequence similarity 13, member A), and AGER (advanced glycosylation end product–specific receptor). Their association with COPD susceptibility has been replicated in multiple populations. Since these candidate genes have not been considered in COPD, their pathological roles are still largely unknown. Herein, we review some evidences that they can be effective drug targets or serve as biomarkers for diagnosis or subtyping. However, more study is required to understand the functional roles of these candidate genes. Future research is needed to characterize the effect of genetic variants, validate gene function in humans and model systems, and elucidate the genes’ transcriptional and posttranscriptional regulatory mechanisms.

Family-Based Association Study of Pulmonary Function in a Population in Northeast Asia

The spirometric measurement of pulmonary function by measuring the forced expiratory volume in one second (FEV1) is a heritable trait that reflects the physiological condition of the lung and airways. Genome-wide linkage and association studies have identified a number of genes and genetic loci associated with pulmonary function. However, limited numbers of studies have been reported for Asian populations. In this study, we aimed to investigate genetic evidence of pulmonary function in a population in northeast Asia. We conducted a family-based association test with 706 GENDISCAN study participants from 72 Mongolian families to determine candidate genetic determinants of pulmonary function. For the replication, we chose seven candidate single nucleotide polymorphisms (SNPs) from the 5 loci, and tested 1062 SNPs for association with FEV1 from 2,729 subjects of the Korea Healthy Twin study. We identified TMEM132C as a potential candidate gene at 12q24.3, which is a previously reported locus of asthma and spirometric indices. We also found two adjacent candidate genes (UNC93A and TTLL2) in the 6q27 region, which has been previously identified as a pulmonary function locus in the Framingham cohort study. Our findings suggest that novel candidate genes (TMEM132C, UNC93A and TTLL2) in two different regions are associated with pulmonary function in a population in northeast Asia.

Novel insights into the genetics of smoking behaviour, lung function, and chronic obstructive pulmonary disease (UK BiLEVE): a genetic association study in UK Biobank

BACKGROUND

Understanding the genetic basis of airflow obstruction and smoking behaviour is key to determining the pathophysiology of chronic obstructive pulmonary disease (COPD). We used UK Biobank data to study the genetic causes of smoking behaviour and lung health.

METHODS

We sampled individuals of European ancestry from UK Biobank, from the middle and extremes of the forced expiratory volume in 1 s (FEV1) distribution among heavy smokers (mean 35 pack-years) and never smokers. We developed a custom array for UK Biobank to provide optimum genome-wide coverage of common and low-frequency variants, dense coverage of genomic regions already implicated in lung health and disease, and to assay rare coding variants relevant to the UK population. We investigated whether there were shared genetic causes between different phenotypes defined by extremes of FEV1. We also looked for novel variants associated with extremes of FEV1 and smoking behaviour and assessed regions of the genome that had already shown evidence for a role in lung health and disease. We set genome-wide significance at p<5 × 10(-8).

FINDINGS

UK Biobank participants were recruited from March 15, 2006, to July 7, 2010. Sample selection for the UK BiLEVE study started on Nov 22, 2012, and was completed on Dec 20, 2012. We selected 50,008 unique samples: 10,002 individuals with low FEV1, 10,000 with average FEV1, and 5002 with high FEV1 from each of the heavy smoker and never smoker groups. We noted a substantial sharing of genetic causes of low FEV1 between heavy smokers and never smokers (p=2.29 × 10(-16)) and between individuals with and without doctor-diagnosed asthma (p=6.06 × 10(-11)). We discovered six novel genome-wide significant signals of association with extremes of FEV1, including signals at four novel loci (KANSL1, TSEN54, TET2, and RBM19/TBX5) and independent signals at two previously reported loci (NPNT and HLA-DQB1/HLA-DQA2). These variants also showed association with COPD, including in individuals with no history of smoking. The number of copies of a 150 kb region containing the 5' end of KANSL1, a gene that is important for epigenetic gene regulation, was associated with extremes of FEV1. We also discovered five new genome-wide significant signals for smoking behaviour, including a variant in NCAM1 (chromosome 11) and a variant on chromosome 2 (between TEX41 and PABPC1P2) that has a trans effect on expression of NCAM1 in brain tissue.

INTERPRETATION

By sampling from the extremes of the lung function distribution in UK Biobank, we identified novel genetic causes of lung function and smoking behaviour. These results provide new insight into the specific mechanisms underlying airflow obstruction, COPD, and tobacco addiction, and show substantial shared genetic architecture underlying airflow obstruction across individuals, irrespective of smoking behaviour and other airway disease.

FUNDING

Medical Research Council.

Molecular mechanisms underlying variations in lung function: a systems genetics analysis

BACKGROUND

Lung function measures reflect the physiological state of the lung, and are essential to the diagnosis of chronic obstructive pulmonary disease (COPD). The SpiroMeta-CHARGE consortium undertook the largest genome-wide association study (GWAS) so far (n=48,201) for forced expiratory volume in 1 s (FEV1) and the ratio of FEV1 to forced vital capacity (FEV1/FVC) in the general population. The lung expression quantitative trait loci (eQTLs) study mapped the genetic architecture of gene expression in lung tissue from 1111 individuals. We used a systems genetics approach to identify single nucleotide polymorphisms (SNPs) associated with lung function that act as eQTLs and change the level of expression of their target genes in lung tissue; termed eSNPs.

METHODS

The SpiroMeta-CHARGE GWAS results were integrated with lung eQTLs to map eSNPs and the genes and pathways underlying the associations in lung tissue. For comparison, a similar analysis was done in peripheral blood. The lung mRNA expression levels of the eSNP-regulated genes were tested for associations with lung function measures in 727 individuals. Additional analyses identified the pleiotropic effects of eSNPs from the published GWAS catalogue, and mapped enrichment in regulatory regions from the ENCODE project. Finally, the Connectivity Map database was used to identify potential therapeutics in silico that could reverse the COPD lung tissue gene signature.

FINDINGS

SNPs associated with lung function measures were more likely to be eQTLs and vice versa. The integration mapped the specific genes underlying the GWAS signals in lung tissue. The eSNP-regulated genes were enriched for developmental and inflammatory pathways; by comparison, SNPs associated with lung function that were eQTLs in blood, but not in lung, were only involved in inflammatory pathways. Lung function eSNPs were enriched for regulatory elements and were over-represented among genes showing differential expression during fetal lung development. An mRNA gene expression signature for COPD was identified in lung tissue and compared with the Connectivity Map. This in-silico drug repurposing approach suggested several compounds that reverse the COPD gene expression signature, including a nicotine receptor antagonist. These findings represent novel therapeutic pathways for COPD.

INTERPRETATION

The system genetics approach identified lung tissue genes driving the variation in lung function and susceptibility to COPD. The identification of these genes and the pathways in which they are enriched is essential to understand the pathophysiology of airway obstruction and to identify novel therapeutic targets and biomarkers for COPD, including drugs that reverse the COPD gene signature in silico.

FUNDING

The research reported in this article was not specifically funded by any agency. See Acknowledgments for a full list of funders of the lung eQTL study and the Spiro-Meta CHARGE GWAS.

The Effect of Ancestry and Genetic Variation on Lung Function Predictions: What Is "Normal" Lung Function in Diverse Human Populations?

Lung function measures are an invaluable screening test for respiratory health and have been associated with the morbidity and mortality related to different airway disease as well as all-cause mortality. Currently, reference values for spirometric measurements are obtained using equations derived from individual ethnic or racial groups. The rapid expansion of more racially and ethnically diverse populations will challenge current race-based lung function reference equations. Recent international general population studies and ancestry-based genetic studies have found that ancestry and genetic variation are determinants of lung function and have suggested a role for genetic ancestry or gene variants in future lung function reference equations. In this review, we discuss the potential limitations of current lung function reference equations in a global society which is becoming more ethnically, racially, and, thus, genetically diverse. We also focus on how an individual's ancestral background or genetic profile could provide the basis for more accurate, personalized predictions of an individual's baseline lung function.

Integrative pathway genomics of lung function and airflow obstruction

Chronic respiratory disorders are important contributors to the global burden of disease. Genome-wide association studies (GWASs) of lung function measures have identified several trait-associated loci, but explain only a modest portion of the phenotypic variability. We postulated that integrating pathway-based methods with GWASs of pulmonary function and airflow obstruction would identify a broader repertoire of genes and processes influencing these traits. We performed two independent GWASs of lung function and applied gene set enrichment analysis to one of the studies and validated the results using the second GWAS. We identified 131 significantly enriched gene sets associated with lung function and clustered them into larger biological modules involved in diverse processes including development, immunity, cell signaling, proliferation and arachidonic acid. We found that enrichment of gene sets was not driven by GWAS-significant variants or loci, but instead by those with less stringent association P-values. Next, we applied pathway enrichment analysis to a meta-analyzed GWAS of airflow obstruction. We identified several biologic modules that functionally overlapped with those associated with pulmonary function. However, differences were also noted, including enrichment of extracellular matrix (ECM) processes specifically in the airflow obstruction study. Network analysis of the ECM module implicated a candidate gene, matrix metalloproteinase 10 (MMP10), as a putative disease target. We used a knockout mouse model to functionally validate MMP10's role in influencing lung's susceptibility to cigarette smoke-induced emphysema. By integrating pathway analysis with population-based genomics, we unraveled biologic processes underlying pulmonary function traits and identified a candidate gene for obstructive lung disease.