On the Analysis of Genome-Wide Association Studies in Family-Based Designs: A Universal, Robust Analysis Approach and an Application to Four Genome-Wide Association Studies

For genome-wide association studies in family-based designs, we propose a new, universally applicable approach. The new test statistic exploits all available information about the association, while, by virtue of its design, it maintains the same robustness against population admixture as traditional family-based approaches that are based exclusively on the within-family information. The approach is suitable for the analysis of almost any trait type, e.g. binary, continuous, time-to-onset, multivariate, etc., and combinations of those. We use simulation studies to verify all theoretically derived properties of the approach, estimate its power, and compare it with other standard approaches. We illustrate the practical implications of the new analysis method by an application to a lung-function phenotype, forced expiratory volume in one second (FEV1) in 4 genome-wide association studies

CRISPR interference interrogation of COPD GWAS genes reveals the functional significance of desmoplakin in iPSC-derived alveolar epithelial cells

Genome-wide association studies (GWAS) have identified dozens of loci associated with chronic obstructive pulmonary disease (COPD) susceptibility; however, the function of associated genes in the cell type(s) affected in disease remains poorly understood, partly due to a lack of cell models that recapitulate human alveolar biology. Here, we apply CRISPR interference to interrogate the function of nine genes implicated in COPD by GWAS in induced pluripotent stem cell–derived type 2 alveolar epithelial cells (iAT2s). We find that multiple genes implicated by GWAS affect iAT2 function, including differentiation potential, maturation, and/or proliferation. Detailed characterization of the GWAS gene DSP demonstrates that it regulates iAT2 cell-cell junctions, proliferation, mitochondrial function, and response to cigarette smoke–induced injury. Our approach thus elucidates the biological function, as well as disease-relevant consequences of dysfunction, of genes implicated in COPD by GWAS in type 2 alveolar epithelial cells.This work was supported by a CJ Martin Early Career Fellowship from the Australian National Health and Medical Research Council awarded to R.B.W.; NIH grant F30HL147426 awarded to K.M.A.; NIH grants U01TR001810, R01DK101501, and R01DK117940 awarded to A.A.W.; NIH grants R01HL135142, R01HL137927, and R01HL147148 awarded to M.H.C.; and NIH grants R01HL127200 and R01HL148667 awarded to X.Z

GAWMerge expands GWAS sample size and diversity by combining array-based genotyping and whole-genome sequencing

Genome-wide association studies (GWAS) have made impactful discoveries for complex diseases, often by amassing very large sample sizes. Yet, GWAS of many diseases remain underpowered, especially for non-European ancestries. One cost-effective approach to increase sample size is to combine existing cohorts, which may have limited sample size or be case-only, with public controls, but this approach is limited by the need for a large overlap in variants across genotyping arrays and the scarcity of non-European controls. We developed and validated a protocol, Genotyping Array-WGS Merge (GAWMerge), for combining genotypes from arrays and whole-genome sequencing, ensuring complete variant overlap, and allowing for diverse samples like Trans-Omics for Precision Medicine to be used. Our protocol involves phasing, imputation, and filtering. We illustrated its ability to control technology driven artifacts and type-I error, as well as recover known disease-associated signals across technologies, independent datasets, and ancestries in smoking-related cohorts. GAWMerge enables genetic studies to leverage existing cohorts to validly increase sample size and enhance discovery for understudied traits and ancestries

A Genome-Wide Association Study in Chronic Obstructive Pulmonary Disease (COPD): Identification of Two Major Susceptibility Loci

There is considerable variability in the susceptibility of smokers to develop chronic obstructive pulmonary disease (COPD). The only known genetic risk factor is severe deficiency of α1-antitrypsin, which is present in 1–2% of individuals with COPD. We conducted a genome-wide association study (GWAS) in a homogenous case-control cohort from Bergen, Norway (823 COPD cases and 810 smoking controls) and evaluated the top 100 single nucleotide polymorphisms (SNPs) in the family-based International COPD Genetics Network (ICGN; 1891 Caucasian individuals from 606 pedigrees) study. The polymorphisms that showed replication were further evaluated in 389 subjects from the US National Emphysema Treatment Trial (NETT) and 472 controls from the Normative Aging Study (NAS) and then in a fourth cohort of 949 individuals from 127 extended pedigrees from the Boston Early-Onset COPD population. Logistic regression models with adjustments of covariates were used to analyze the case-control populations. Family-based association analyses were conducted for a diagnosis of COPD and lung function in the family populations. Two SNPs at the α-nicotinic acetylcholine receptor (CHRNA 3/5) locus were identified in the genome-wide association study. They showed unambiguous replication in the ICGN family-based analysis and in the NETT case-control analysis with combined p-values of 1.48×10−10, (rs8034191) and 5.74×10−10 (rs1051730). Furthermore, these SNPs were significantly associated with lung function in both the ICGN and Boston Early-Onset COPD populations. The C allele of the rs8034191 SNP was estimated to have a population attributable risk for COPD of 12.2%. The association of hedgehog interacting protein (HHIP) locus on chromosome 4 was also consistently replicated, but did not reach genome-wide significance levels. Genome-wide significant association of the HHIP locus with lung function was identified in the Framingham Heart study (Wilk et al., companion article in this issue of PLoS Genetics; doi:10.1371/journal.pgen.1000429). The CHRNA 3/5 and the HHIP loci make a significant contribution to the risk of COPD. CHRNA3/5 is the same locus that has been implicated in the risk of lung cancer

Hemizygous Deletion on Chromosome 3p26.1 Is Associated with Heavy Smoking among African American Subjects in the COPDGene Study

Many well-powered genome-wide association studies have identified genetic determinants of self-reported smoking behaviors and measures of nicotine dependence, but most have not considered the role of structural variants, such as copy number variation (CNVs), influencing these phenotypes. Here, we included 2,889 African American and 6,187 non-Hispanic White subjects from the COPDGene cohort (http://www.copdgene.org) to carefully investigate the role of polymorphic CNVs across the genome on various measures of smoking behavior. We identified a CNV component (a hemizygous deletion) on chromosome 3p26.1 associated with two quantitative phenotypes related to smoking behavior among African Americans. This polymorphic hemizygous deletion is significantly associated with pack-years and cigarettes smoked per day among African American subjects in the COPDGene study. We sought evidence of replication in African Americans from the population based Atherosclerosis Risk in Communities (ARIC) study. While we observed similar CNV counts, the extent of exposure to cigarette smoking among ARIC subjects was quite different and the smaller sample size of heavy smokers in ARIC severely limited statistical power, so we were unable to replicate our findings from the COPDGene cohort. But meta-analyses of COPDGene and ARIC study subjects strengthened our association signal. However, a few linkage studies have reported suggestive linkage to the 3p26.1 region, and a few genome-wide association studies (GWAS) have reported markers in the gene (GRM7) nearest to this 3p26.1 area of polymorphic deletions are associated with measures of nicotine dependence among subjects of European ancestry

Folliculin mutations are not associated with severe COPD

<p>Abstract</p> <p>Background</p> <p>Rare loss-of-function folliculin (<it>FLCN</it>) mutations are the genetic cause of Birt-Hogg-Dubé syndrome, a monogenic disorder characterized by spontaneous pneumothorax, fibrofolliculomas, and kidney tumors. Loss-of-function folliculin mutations have also been described in pedigrees with familial spontaneous pneumothorax. Because the majority of patients with folliculin mutations have radiographic evidence of pulmonary cysts, folliculin has been hypothesized to contribute to the development of emphysema.</p> <p>To determine whether folliculin sequence variants are risk factors for severe COPD, we genotyped seven previously reported Birt-Hogg-Dubé or familial spontaneous pneumothorax associated folliculin mutations in 152 severe COPD probands participating in the Boston Early-Onset COPD Study. We performed bidirectional resequencing of all 14 folliculin exons in a subset of 41 probands and subsequently genotyped four identified variants in an independent sample of345 COPD subjects from the National Emphysema Treatment Trial (cases) and 420 male smokers with normal lung function from the Normative Aging Study (controls).</p> <p>Results</p> <p>None of the seven previously reported Birt-Hogg-Dubé or familial spontaneous pneumothorax mutations were observed in the 152 severe, early-onset COPD probands. Exon resequencing identified 31 variants, including two non-synonymous polymorphisms and two common non-coding polymorphisms. No significant association was observed for any of these four variants with presence of COPD or emphysema-related phenotypes.</p> <p>Conclusion</p> <p>Genetic variation in folliculin does not appear to be a major risk factor for severe COPD. These data suggest that familial spontaneous pneumothorax and COPD have distinct genetic causes, despite some overlap in radiographic characteristics.</p

Characteristics of stable chronic obstructive pulmonary disease patients in the pulmonology clinics of seven Asian cities

Background and objectives Chronic obstructive pulmonary disease (COPD) is responsible for significant morbidity and mortality worldwide. We evaluated the characteristics of stable COPD patients in the pulmonology clinics of seven Asian cities and also evaluated whether the exposure to biomass fuels and dusty jobs were related to respiratory symptoms, airflow limitation, and quality of life in the COPD patients. Methods: This cross-sectional observational study recruited 922 COPD patients from seven cities of Asia. The patients underwent spirometry and were administered questionnaires about their exposure to cigarette smoking, biomass fuels, and dusty jobs in addition to respiratory symptoms and health related quality of life. Results: Of the patients, there appeared to be variations from city to city in the history of exposure to biomass fuels and dusty jobs and also in respiratory symptoms of cough, phlegm, wheeze, and dyspnea. These symptoms were more frequent in those COPD patients with a history of exposure to biomass fuels than without and those with a history of exposure to dusty jobs than without (P < 0.01 for all comparisons). Airflow limitation was more severe in those COPD patients with a history of exposure to biomass fuels than without (52.2% predicted versus 55.9% of post-bronchodilator forced expiratory volume in 1 second [FEV1], P = 0.009); quality of life was poorer in those with exposure to biomass fuels than without (40.4 versus 36.2 of the St George’s Respiratory Questionnaire [SGRQ] total score, P = 0.001). Airflow limitation was more severe in those COPD patients with a history of exposure to dusty jobs than without (51.2% predicted versus 57.3% of post-bronchodilator FEV1, P < 0.001); quality of life was poorer in those with dusty jobs than without (41.0 versus 34.6 of SGRQ score, P = 0.006). Conclusion: In Asian cities, the characteristics of COPD patients vary and the history of exposure to biomass fuels or dusty jobs was related to frequency of symptoms, severe airflow limitation, and poor quality of life

A Genome-Wide Association Study of Pulmonary Function Measures in the Framingham Heart Study

The ratio of forced expiratory volume in one second to forced vital capacity (FEV1/FVC) is a measure used to diagnose airflow obstruction and is highly heritable. We performed a genome-wide association study in 7,691 Framingham Heart Study participants to identify single-nucleotide polymorphisms (SNPs) associated with the FEV1/FVC ratio, analyzed as a percent of the predicted value. Identified SNPs were examined in an independent set of 835 Family Heart Study participants enriched for airflow obstruction. Four SNPs in tight linkage disequilibrium on chromosome 4q31 were associated with the percent predicted FEV1/FVC ratio with p-values of genome-wide significance in the Framingham sample (best p-value = 3.6e-09). One of the four chromosome 4q31 SNPs (rs13147758; p-value 2.3e-08 in Framingham) was genotyped in the Family Heart Study and produced evidence of association with the same phenotype, percent predicted FEV1/FVC (p-value = 2.0e-04). The effect estimates for association in the Framingham and Family Heart studies were in the same direction, with the minor allele (G) associated with higher FEV1/FVC ratio levels. Results from the Family Heart Study demonstrated that the association extended to FEV1 and dichotomous airflow obstruction phenotypes, particularly among smokers. The SNP rs13147758 was associated with the percent predicted FEV1/FVC ratio in independent samples from the Framingham and Family Heart Studies producing a combined p-value of 8.3e-11, and this region of chromosome 4 around 145.68 megabases was associated with COPD in three additional populations reported in the accompanying manuscript. The associated SNPs do not lie within a gene transcript but are near the hedgehog-interacting protein (HHIP) gene and several expressed sequence tags cloned from fetal lung. Though it is unclear what gene or regulatory effect explains the association, the region warrants further investigation

On the Analysis of Genome-Wide Association Studies in Family-Based Designs: A Universal, Robust Analysis Approach and an Application to Four Genome-Wide Association Studies

For genome-wide association studies in family-based designs, we propose a new, universally applicable approach. The new test statistic exploits all available information about the association, while, by virtue of its design, it maintains the same robustness against population admixture as traditional family-based approaches that are based exclusively on the within-family information. The approach is suitable for the analysis of almost any trait type, e.g. binary, continuous, time-to-onset, multivariate, etc., and combinations of those. We use simulation studies to verify all theoretically derived properties of the approach, estimate its power, and compare it with other standard approaches. We illustrate the practical implications of the new analysis method by an application to a lung-function phenotype, forced expiratory volume in one second (FEV1) in 4 genome-wide association studies

Asthma-susceptibility variants identified using probands in case-control and family-based analyses

<p>Abstract</p> <p>Background</p> <p>Asthma is a chronic respiratory disease whose genetic basis has been explored for over two decades, most recently via genome-wide association studies. We sought to find asthma-susceptibility variants by using probands from a single population in both family-based and case-control association designs.</p> <p>Methods</p> <p>We used probands from the Childhood Asthma Management Program (CAMP) in two primary genome-wide association study designs: (1) probands were combined with publicly available population controls in a case-control design, and (2) probands and their parents were used in a family-based design. We followed a two-stage replication process utilizing three independent populations to validate our primary findings.</p> <p>Results</p> <p>We found that single nucleotide polymorphisms with similar case-control and family-based association results were more likely to replicate in the independent populations, than those with the smallest p-values in either the case-control or family-based design alone. The single nucleotide polymorphism that showed the strongest evidence for association to asthma was rs17572584, which replicated in 2/3 independent populations with an overall p-value among replication populations of 3.5E-05. This variant is near a gene that encodes an enzyme that has been implicated to act coordinately with modulators of Th2 cell differentiation and is expressed in human lung.</p> <p>Conclusions</p> <p>Our results suggest that using probands from family-based studies in case-control designs, and combining results of both family-based and case-control approaches, may be a way to augment our ability to find SNPs associated with asthma and other complex diseases.</p