Polygenic risk scores identify heterogeneity in asthma and chronic obstructive pulmonary disease.

BACKGROUND: Asthma and chronic obstructive pulmonary disease (COPD) have distinct and overlapping genetic and clinical features. OBJECTIVE: We sought to test the hypothesis that polygenic risk scores (PRSs) for asthma (PRSAsthma) and spirometry (FEV1 and FEV1/forced vital capacity; PRSspiro) would demonstrate differential associations with asthma, COPD, and asthma-COPD overlap (ACO). METHODS: We developed and tested 2 asthma PRSs and applied the higher performing PRSAsthma and a previously published PRSspiro to research (Genetic Epidemiology of COPD study and Childhood Asthma Management Program, with spirometry) and electronic health record-based (Mass General Brigham Biobank and Genetic Epidemiology Research on Adult Health and Aging [GERA]) studies. We assessed the association of PRSs with COPD and asthma using modified random-effects and binary-effects meta-analyses, and ACO and asthma exacerbations in specific cohorts. Models were adjusted for confounders and genetic ancestry. RESULTS: In meta-analyses of 102,477 participants, the PRSAsthma (odds ratio [OR] per SD, 1.16 [95% CI, 1.14-1.19]) and PRSspiro (OR per SD, 1.19 [95% CI, 1.17-1.22]) both predicted asthma, whereas the PRSspiro predicted COPD (OR per SD, 1.25 [95% CI, 1.21-1.30]). However, results differed by cohort. The PRSspiro was not associated with COPD in GERA and Mass General Brigham Biobank. In the Genetic Epidemiology of COPD study, the PRSAsthma (OR per SD: Whites, 1.3; African Americans, 1.2) and PRSspiro (OR per SD: Whites, 2.2; African Americans, 1.6) were both associated with ACO. In GERA, the PRSAsthma was associated with asthma exacerbations (OR, 1.18) in Whites; the PRSspiro was associated with asthma exacerbations in White, LatinX, and East Asian participants. CONCLUSIONS: PRSs for asthma and spirometry are both associated with ACO and asthma exacerbations. Genetic prediction performance differs in research versus electronic health record-based cohorts.

Intrauterine Smoke Exposure, microRNA Expression during Human Lung Development, and Childhood Asthma.

Intrauterine smoke (IUS) exposure during early childhood has been associated with a number of negative health consequences, including reduced lung function and asthma susceptibility. The biological mechanisms underlying these associations have not been established. MicroRNAs regulate the expression of numerous genes involved in lung development. Thus, investigation of the impact of IUS on miRNA expression during human lung development may elucidate the impact of IUS on post-natal respiratory outcomes. We sought to investigate the effect of IUS exposure on miRNA expression during early lung development. We hypothesized that miRNA-mRNA networks are dysregulated by IUS during human lung development and that these miRNAs may be associated with future risk of asthma and allergy. Human fetal lung samples from a prenatal tissue retrieval program were tested for differential miRNA expression with IUS exposure (measured using placental cotinine concentration). RNA was extracted and miRNA-sequencing was performed. We performed differential expression using IUS exposure, with covariate adjustment. We also considered the above model with an additional sex-by-IUS interaction term, allowing IUS effects to differ by male and female samples. Using paired gene expression profiles, we created sex-stratified miRNA-mRNA correlation networks predictive of IUS using DIABLO. We additionally evaluated whether miRNAs were associated with asthma and allergy outcomes in a cohort of childhood asthma. We profiled pseudoglandular lung miRNA in n = 298 samples, 139 (47%) of which had evidence of IUS exposure. Of 515 miRNAs, 25 were significantly associated with intrauterine smoke exposure (q-value < 0.10). The IUS associated miRNAs were correlated with well-known asthma genes (e.g., ORM1-Like Protein 3, ORDML3) and enriched in disease-relevant pathways (oxidative stress). Eleven IUS-miRNAs were also correlated with clinical measures (e.g., Immunoglobulin E andlungfunction) in children with asthma, further supporting their likely disease relevance. Lastly, we found substantial differences in IUS effects by sex, finding 95 significant IUS-miRNAs in male samples, but only four miRNAs in female samples. The miRNA-mRNA correlation networks were predictive of IUS (AUC = 0.78 in males and 0.86 in females) and suggested that IUS-miRNAs are involved in regulation of disease-relevant genes (e.g., A disintegrin and metalloproteinase domain 19 (ADAM19), LBH regulator of WNT signaling (LBH)) and sex hormone signaling (Coactivator associated methyltransferase 1(CARM1)). Our study demonstrated differential expression of miRNAs by IUS during early prenatal human lung development, which may be modified by sex. Based on their gene targets and correlation to clinical asthma and atopy outcomes, these IUS-miRNAs may be relevant for subsequent allergy and asthma risk. Our study provides insight into the impact of IUS in human fetal lung transcriptional networks and on the developmental origins of asthma and allergic disorders.

Lung tissue shows divergent gene expression between chronic obstructive pulmonary disease and idiopathic pulmonary fibrosis.

BACKGROUND: Chronic obstructive pulmonary disease (COPD) and idiopathic pulmonary fibrosis (IPF) are characterized by shared exposures and clinical features, but distinct genetic and pathologic features exist. These features have not been well-studied using large-scale gene expression datasets. We hypothesized that there are divergent gene, pathway, and cellular signatures between COPD and IPF. METHODS: We performed RNA-sequencing on lung tissues from individuals with IPF (n = 231) and COPD (n = 377) compared to control (n = 267), defined as individuals with normal spirometry. We grouped the overlapping differential expression gene sets based on direction of expression and examined the resultant sets for genes of interest, pathway enrichment, and cell composition. Using gene set variation analysis, we validated the overlap group gene sets in independent COPD and IPF data sets. RESULTS: We found 5010 genes differentially expressed between COPD and control, and 11,454 genes differentially expressed between IPF and control (1% false discovery rate). 3846 genes overlapped between IPF and COPD. Several pathways were enriched for genes upregulated in COPD and downregulated in IPF; however, no pathways were enriched for genes downregulated in COPD and upregulated in IPF. There were many myeloid cell genes with increased expression in COPD but decreased in IPF. We found that the genes upregulated in COPD but downregulated in IPF were associated with lower lung function in the independent validation cohorts. CONCLUSIONS: We identified a divergent gene expression signature between COPD and IPF, with increased expression in COPD and decreased in IPF. This signature is associated with worse lung function in both COPD and IPF.

A unifying framework for rare variant association testing in family-based designs, including higher criticism approaches, SKATs, and burden tests.

MOTIVATION: Analysis of rare variants in family-based studies remains a challenge. Transmission-based approaches provide robustness against population stratification, but the evaluation of the significance of test statistics based on asymptotic theory can be imprecise. Also, power will depend heavily on the choice of the test statistic and on the underlying genetic architecture of the locus, which will be generally unknown. RESULTS: In our proposed framework, we utilize the FBAT haplotype algorithm to obtain the conditional offspring genotype distribution under the null hypothesis given the sufficient statistic. Based on this conditional offspring genotype distribution, the significance of virtually any association test statistic can be evaluated based on simulations or exact computations, without the need for asymptotic approximations. Besides standard linear burden-type statistics, this enables our approach to also evaluate other test statistics such as variance components statistics, higher criticism approaches, and maximum-single-variant-statistics, where asymptotic theory might be involved or does not provide accurate approximations for rare variant data. Based on these P-values, combined test statistics such as the aggregated Cauchy association test (ACAT) can also be utilized. In simulation studies, we show that our framework outperforms existing approaches for family-based studies in several scenarios. We also applied our methodology to a TOPMed whole-genome sequencing dataset with 897 asthmatic trios from Costa Rica. AVAILABILITY AND IMPLEMENTATION: FBAT software is available at https://sites.google.com/view/fbatwebpage. Simulation code is available at https://github.com/julianhecker/FBAT_rare_variant_test_simulations. Whole-genome sequencing data for 'NHLBI TOPMed: The Genetic Epidemiology of Asthma in Costa Rica' is available at https://www.ncbi.nlm.nih.gov/projects/gap/cgi-bin/study.cgi?study_id=phs000988.v4.p1. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

De novo mutations across 1,465 diverse genomes reveal mutational insights and reductions in the Amish founder population.

De novo mutations (DNMs), or mutations that appear in an individual despite not being seen in their parents, are an important source of genetic variation whose impact is relevant to studies of human evolution, genetics, and disease. Utilizing high-coverage whole-genome sequencing data as part of the Trans-Omics for Precision Medicine (TOPMed) Program, we called 93,325 single-nucleotide DNMs across 1,465 trios from an array of diverse human populations, and used them to directly estimate and analyze DNM counts, rates, and spectra. We find a significant positive correlation between local recombination rate and local DNM rate, and that DNM rate explains a substantial portion (8.98 to 34.92%, depending on the model) of the genome-wide variation in population-level genetic variation from 41K unrelated TOPMed samples. Genome-wide heterozygosity does correlate with DNM rate, but only explains <1% of variation. While we are underpowered to see small differences, we do not find significant differences in DNM rate between individuals of European, African, and Latino ancestry, nor across ancestrally distinct segments within admixed individuals. However, we did find significantly fewer DNMs in Amish individuals, even when compared with other Europeans, and even after accounting for parental age and sequencing center. Specifically, we found significant reductions in the number of C→A and T→C mutations in the Amish, which seem to underpin their overall reduction in DNMs. Finally, we calculated near-zero estimates of narrow sense heritability (h2), which suggest that variation in DNM rate is significantly shaped by nonadditive genetic effects and the environment.

Alpha-1 Antitrypsin PiMZ Genotype Is Associated with Chronic Obstructive Pulmonary Disease in Two Racial Groups.

RATIONALE: Alpha-1 antitrypsin deficiency, caused primarily by homozygosity for the Z allele of the SERPINA1 gene, is a well-established genetic cause of chronic obstructive pulmonary disease (COPD). Whether the heterozygous PiMZ genotype for alpha-1 antitrypsin confers increased risk for COPD has been debated. OBJECTIVES: We analyzed 8,271 subjects in the Genetic Epidemiology of COPD (COPDGene) Study, hypothesizing that PiMZ would independently associate with COPD and COPD-related phenotypes. METHODS: The COPDGene Study comprises a multiethnic, cross-sectional, observational cohort of non-Hispanic white and African American current and former smokers with at least 10 pack-years of smoking who were enrolled for detailed clinical and genetic studies of COPD and COPD-related traits. We performed multivariate logistic regression analysis for moderate to severe COPD and assessed Pi genotype with other relevant covariates in models stratified by race. We analyzed quantitative characteristics on the basis of volumetric computed tomography with generalized linear models controlling for genotype, scanner type, and similar covariates. RESULTS: White PiMZ COPDGene subjects had significantly lower lung function, FEV1 percent predicted (68 ± 28 vs. 75 ± 27; P = 0.0005), and FEV1/FVC ratio (0.59 ± 0.18 vs. 0.63 ± 0.17; P = 0.0008), as well as more radiographic emphysema (P = 0.001), than subjects without alpha-1 antitrypsin Z risk alleles. Similarly, African American PiMZ subjects had lower lung function, FEV1 percent predicted (65 ± 33 vs. 84 ± 25; P = 0.009) and FEV1/FVC (0.61 ± 0.21 vs. 0.71 ± 0.15; P = 0.03). CONCLUSIONS: In the COPDGene Study, we demonstrate that PiMZ heterozygous individuals who smoke are at increased risk for COPD and obstructive lung function impairment compared with Z-allele noncarriers, regardless of race. Although severe alpha-1 antitrypsin deficiency is uncommon in African Americans, our study adds further support for initial targeted detection of all subjects with COPD for alpha-1 antitrypsin deficiency, including African Americans. Clinical trial registered with www.clinicaltrials.gov (NCT00608784).

Factors influencing the infant gut microbiome at age 3-6 months: Findings from the ethnically diverse Vitamin D Antenatal Asthma Reduction Trial (VDAART).

BACKGROUND: The gut microbiome in infancy influences immune system maturation, and may have an important impact on allergic disease risk. OBJECTIVE: We sought to determine how prenatal and early life factors impact the gut microbiome in a relatively large, ethnically diverse study population of infants at age 3 to 6 months, who were enrolled in Vitamin D Antenatal Asthma Reduction Trial, a clinical trial of vitamin D supplementation in pregnancy to prevent asthma and allergies in offspring. METHODS: We performed 16S rRNA gene sequencing on 333 infants' stool samples. Microbial diversity was computed using the Shannon index. Factor analysis applied to the top 25 most abundant taxa revealed 4 underlying bacterial coabundance groups; the first dominated by Firmicutes (Lachnospiraceae/Clostridiales), the second by Proteobacteria (Klebsiella/Enterobacter), the third by Bacteriodetes, and the fourth by Veillonella. Scores for coabundance groups were used as outcomes in regression models, with prenatal/birth and demographic characteristics as independent predictors. Multivariate analysis, using all microbial community members, was also conducted. RESULTS: White race/ethnicity was associated with lower diversity but higher Bacteroidetes coabundance scores. C-section birth was associated with higher diversity, but decreased Bacteroidetes coabundance scores. Firmicutes scores were higher for infants born by C-section. Breast-fed infants had lower proportions of Clostridiales. Cord blood vitamin D was linked to increased Lachnobacterium, but decreased Lactococcus. CONCLUSIONS: The findings presented here suggest that race, mode of delivery, breast-feeding, and cord blood vitamin D levels are associated with infant gut microbiome composition, with possible long-term implications for immune system modulation and asthma/allergic disease incidence.

Expression Quantitative Trait Loci Information Improves Predictive Modeling of Disease Relevance of Non-Coding Genetic Variation.

Disease-associated loci identified through genome-wide association studies (GWAS) frequently localize to non-coding sequence. We and others have demonstrated strong enrichment of such single nucleotide polymorphisms (SNPs) for expression quantitative trait loci (eQTLs), supporting an important role for regulatory genetic variation in complex disease pathogenesis. Herein we describe our initial efforts to develop a predictive model of disease-associated variants leveraging eQTL information. We first catalogued cis-acting eQTLs (SNPs within 100 kb of target gene transcripts) by meta-analyzing four studies of three blood-derived tissues (n = 586). At a false discovery rate < 5%, we mapped eQTLs for 6,535 genes; these were enriched for disease-associated genes (P < 10(-04)), particularly those related to immune diseases and metabolic traits. Based on eQTL information and other variant annotations (distance from target gene transcript, minor allele frequency, and chromatin state), we created multivariate logistic regression models to predict SNP membership in reported GWAS. The complete model revealed independent contributions of specific annotations as strong predictors, including evidence for an eQTL (odds ratio (OR) = 1.2-2.0, P < 10(-11)) and the chromatin states of active promoters, different classes of strong or weak enhancers, or transcriptionally active regions (OR = 1.5-2.3, P < 10(-11)). This complete prediction model including eQTL association information ultimately allowed for better discrimination of SNPs with higher probabilities of GWAS membership (6.3-10.0%, compared to 3.5% for a random SNP) than the other two models excluding eQTL information. This eQTL-based prediction model of disease relevance can help systematically prioritize non-coding GWAS SNPs for further functional characterization.

CMTR1 is associated with increased asthma exacerbations in patients taking inhaled corticosteroids.

Inhaled corticosteroids (ICS) are the most effective controller medications for asthma, and variability in ICS response is associated with genetic variation. Despite ICS treatment, some patients with poor asthma control experience severe asthma exacerbations, defined as a hospitalization or emergency room visit. We hypothesized that some individuals may be at increased risk of asthma exacerbations, despite ICS use, due to genetic factors. A GWAS of 237,726 common, independent markers was conducted in 806 Caucasian asthmatic patients from two population-based biobanks: BioVU, at Vanderbilt University Medical Center (VUMC) in Tennessee (369 patients), and Personalized Medicine Research Project (PMRP) at the Marshfield Clinic in Wisconsin (437 patients). Using a case-control study design, the association of each SNP locus with the outcome of asthma exacerbations (defined as asthma-related emergency department visits or hospitalizations concurrent with oral corticosteroid use), was evaluated for each population by logistic regression analysis, adjusting for age, gender and the first four principal components. A meta-analysis of the results was conducted. Validation of expression of selected candidate genes was determined by evaluating an independent microarray expression data set. Our study identified six novel SNPs associated with differential risk of asthma exacerbations (P < 10(-05)). The top GWAS result, rs2395672 in CMTR1, was associated with an increased risk of exacerbations in both populations (OR = 1.07, 95% CI 1.03-1.11; joint P = 2.3 × 10(-06)). Two SNPs (rs2395672 and rs279728) were associated with increased risk of exacerbations, while the remaining four SNPs (rs4271056, rs6467778, rs2691529, and rs9303988) were associated with decreased risk. Three SNPs (rs2395672, rs6467778, and rs2691529) were present in three genes: CMTR1, TRIM24 and MAGI2. The CMTR1 mRNA transcript was significantly differentially expressed in nasal lavage samples from asthmatics during acute exacerbations, suggesting potential involvement of this gene in the development of this phenotype. We show that genetic variability may contribute to asthma exacerbations in patients taking ICS. Furthermore, our studies implicate CMTR1 as a novel candidate gene with potential roles in the pathogenesis of asthma exacerbations.

Integrated genome-wide association, coexpression network, and expression single nucleotide polymorphism analysis identifies novel pathway in allergic rhinitis.

BACKGROUND: Allergic rhinitis is a common disease whose genetic basis is incompletely explained. We report an integrated genomic analysis of allergic rhinitis. METHODS: We performed genome wide association studies (GWAS) of allergic rhinitis in 5633 ethnically diverse North American subjects. Next, we profiled gene expression in disease-relevant tissue (peripheral blood CD4+ lymphocytes) collected from subjects who had been genotyped. We then integrated the GWAS and gene expression data using expression single nucleotide (eSNP), coexpression network, and pathway approaches to identify the biologic relevance of our GWAS. RESULTS: GWAS revealed ethnicity-specific findings, with 4 genome-wide significant loci among Latinos and 1 genome-wide significant locus in the GWAS meta-analysis across ethnic groups. To identify biologic context for these results, we constructed a coexpression network to define modules of genes with similar patterns of CD4+ gene expression (coexpression modules) that could serve as constructs of broader gene expression. 6 of the 22 GWAS loci with P-value ≤ 1x10-6 tagged one particular coexpression module (4.0-fold enrichment, P-value 0.0029), and this module also had the greatest enrichment (3.4-fold enrichment, P-value 2.6 × 10-24) for allergic rhinitis-associated eSNPs (genetic variants associated with both gene expression and allergic rhinitis). The integrated GWAS, coexpression network, and eSNP results therefore supported this coexpression module as an allergic rhinitis module. Pathway analysis revealed that the module was enriched for mitochondrial pathways (8.6-fold enrichment, P-value 4.5 × 10-72). CONCLUSIONS: Our results highlight mitochondrial pathways as a target for further investigation of allergic rhinitis mechanism and treatment. Our integrated approach can be applied to provide biologic context for GWAS of other diseases.

Risk loci for chronic obstructive pulmonary disease: a genome-wide association study and meta-analysis.

BACKGROUND: The genetic risk factors for susceptibility to chronic obstructive pulmonary disease (COPD) are still largely unknown. Additional genetic variants are likely to be identified by genome-wide association studies in larger cohorts or specific subgroups. We sought to identify risk loci for moderate to severe and severe COPD with data from several cohort studies. METHODS: We combined genome-wide association analysis data from participants in the COPDGene study (non-Hispanic white and African-American ethnic origin) and the ECLIPSE, NETT/NAS, and Norway GenKOLS studies (self-described white ethnic origin). We did analyses comparing control individuals with individuals with moderate to severe COPD and with a subset of individuals with severe COPD. Single nucleotide polymorphisms yielding a p value of less than 5 × 10(-7) in the meta-analysis at loci not previously described were genotyped in individuals from the family-based ICGN study. We combined results in a joint meta-analysis (threshold for significance p<5 × 10(-8)). FINDINGS: Analysis of 6633 individuals with moderate to severe COPD and 5704 control individuals confirmed association at three known loci: CHRNA3 (p=6·38 × 10(-14)), FAM13A (p=1·12 × 10(-14)), and HHIP (p=1·57 × 10(-12)). We also showed significant evidence of association at a novel locus near RIN3 (p=5·25 × 10(-9)). In the overall meta-analysis (ie, including data from 2859 ICGN participants), the association with RIN3 remained significant (p=5·4 × 10(-9)). 3497 individuals were included in our analysis of severe COPD. The effect estimates for the loci near HHIP and CHRNA3 were significantly stronger in severe disease than in moderate to severe disease (p<0·01). We also identified associations at two additional loci: MMP12 (overall joint meta-analysis p=2·6 × 10(-9)) and TGFB2 (overall joint meta-analysis p=8·3 × 10(-9)). INTERPRETATION: We have confirmed associations with COPD at three known loci and identified three new genome-wide significant associations. Genetic variants other than in α-1 antitrypsin increase the risk of COPD. FUNDING: US National Heart, Lung, and Blood Institute; the Alpha-1 Foundation; the COPD Foundation through contributions from AstraZeneca, Boehringer Ingelheim, Novartis, and Sepracor; GlaxoSmithKline; Centers for Medicare and Medicaid Services; Agency for Healthcare Research and Quality; and US Department of Veterans Affairs.

ITGB5 and AGFG1 variants are associated with severity of airway responsiveness.

BACKGROUND: Airway hyperresponsiveness (AHR), a primary characteristic of asthma, involves increased airway smooth muscle contractility in response to certain exposures. We sought to determine whether common genetic variants were associated with AHR severity. METHODS: A genome-wide association study (GWAS) of AHR, quantified as the natural log of the dosage of methacholine causing a 20% drop in FEV1, was performed with 994 non-Hispanic white asthmatic subjects from three drug clinical trials: CAMP, CARE, and ACRN. Genotyping was performed on Affymetrix 6.0 arrays, and imputed data based on HapMap Phase 2, was used to measure the association of SNPs with AHR using a linear regression model. Replication of primary findings was attempted in 650 white subjects from DAG, and 3,354 white subjects from LHS. Evidence that the top SNPs were eQTL of their respective genes was sought using expression data available for 419 white CAMP subjects. RESULTS: The top primary GWAS associations were in rs848788 (P-value 7.2E-07) and rs6731443 (P-value 2.5E-06), located within the ITGB5 and AGFG1 genes, respectively. The AGFG1 result replicated at a nominally significant level in one independent population (LHS P-value 0.012), and the SNP had a nominally significant unadjusted P-value (0.0067) for being an eQTL of AGFG1. CONCLUSIONS: Based on current knowledge of ITGB5 and AGFG1, our results suggest that variants within these genes may be involved in modulating AHR. Future functional studies are required to confirm that our associations represent true biologically significant findings.

Integration of mouse and human genome-wide association data identifies KCNIP4 as an asthma gene.

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

Genome-wide association analysis in asthma subjects identifies SPATS2L as a novel bronchodilator response gene.

Bronchodilator response (BDR) is an important asthma phenotype that measures reversibility of airway obstruction by comparing lung function (i.e. FEV(1)) before and after the administration of a short-acting ß(2)-agonist, the most common rescue medications used for the treatment of asthma. BDR also serves as a test of ß(2)-agonist efficacy. BDR is a complex trait that is partly under genetic control. A genome-wide association study (GWAS) of BDR, quantified as percent change in baseline FEV(1) after administration of a ß(2)-agonist, was performed with 1,644 non-Hispanic white asthmatic subjects from six drug clinical trials: CAMP, LOCCS, LODO, a medication trial conducted by Sepracor, CARE, and ACRN. Data for 469,884 single-nucleotide polymorphisms (SNPs) were used to measure the association of SNPs with BDR using a linear regression model, while adjusting for age, sex, and height. Replication of primary P-values was attempted in 501 white subjects from SARP and 550 white subjects from DAG. Experimental evidence supporting the top gene was obtained via siRNA knockdown and Western blotting analyses. The lowest overall combined P-value was 9.7E-07 for SNP rs295137, near the SPATS2L gene. Among subjects in the primary analysis, those with rs295137 TT genotype had a median BDR of 16.0 (IQR = [6.2, 32.4]), while those with CC or TC genotypes had a median BDR of 10.9 (IQR = [5.0, 22.2]). SPATS2L mRNA knockdown resulted in increased ß(2)-adrenergic receptor levels. Our results suggest that SPATS2L may be an important regulator of ß(2)-adrenergic receptor down-regulation and that there is promise in gaining a better understanding of the biological mechanisms of differential response to ß(2)-agonists through GWAS.

Genome-wide association study of the age of onset of childhood asthma.

BACKGROUND: Childhood asthma is a complex disease with known heritability and phenotypic diversity. Although an earlier onset has been associated with more severe disease, there has been no genome-wide association study of the age of onset of asthma in children. OBJECTIVE: We sought to identify genetic variants associated with earlier onset of childhood asthma. METHODS: We conducted the first genome-wide association study of the age of onset of childhood asthma among participants in the Childhood Asthma Management Program (CAMP) and used 3 independent cohorts from North America, Costa Rica, and Sweden for replication. RESULTS: Two single nucleotide polymorphisms (SNPs) were associated with earlier onset of asthma in the combined analysis of CAMP and the replication cohorts: rs9815663 (Fisher P= 2.31 × 10(-8)) and rs7927044 (P= 6.54 × 10(-9)). Of these 2 SNPs, rs9815663 was also significantly associated with earlier asthma onset in an analysis including only the replication cohorts. Ten SNPs in linkage disequilibrium with rs9815663 were also associated with earlier asthma onset (2.24 × 10(-7)

A genome-wide association study of COPD identifies a susceptibility locus on chromosome 19q13.

The genetic risk factors for chronic obstructive pulmonary disease (COPD) are still largely unknown. To date, genome-wide association studies (GWASs) of limited size have identified several novel risk loci for COPD at CHRNA3/CHRNA5/IREB2, HHIP and FAM13A; additional loci may be identified through larger studies. We performed a GWAS using a total of 3499 cases and 1922 control subjects from four cohorts: the Evaluation of COPD Longitudinally to Identify Predictive Surrogate Endpoints (ECLIPSE); the Normative Aging Study (NAS) and National Emphysema Treatment Trial (NETT); Bergen, Norway (GenKOLS); and the COPDGene study. Genotyping was performed on Illumina platforms with additional markers imputed using 1000 Genomes data; results were summarized using fixed-effect meta-analysis. We identified a new genome-wide significant locus on chromosome 19q13 (rs7937, OR = 0.74, P = 2.9 × 10(-9)). Genotyping this single nucleotide polymorphism (SNP) and another nearby SNP in linkage disequilibrium (rs2604894) in 2859 subjects from the family-based International COPD Genetics Network study (ICGN) demonstrated supportive evidence for association for COPD (P = 0.28 and 0.11 for rs7937 and rs2604894), pre-bronchodilator FEV(1) (P = 0.08 and 0.04) and severe (GOLD 3&4) COPD (P = 0.09 and 0.017). This region includes RAB4B, EGLN2, MIA and CYP2A6, and has previously been identified in association with cigarette smoking behavior.

Identification of FGF7 as a novel susceptibility locus for chronic obstructive pulmonary disease.

RATIONALE: Traditional genome-wide association studies (GWASs) of large cohorts of subjects with chronic obstructive pulmonary disease (COPD) have successfully identified novel candidate genes, but several other plausible loci do not meet strict criteria for genome-wide significance after correction for multiple testing. OBJECTIVES: The authors hypothesise that by applying unbiased weights derived from unique populations we can identify additional COPD susceptibility loci. Methods The authors performed a homozygosity haplotype analysis on a group of subjects with and without COPD to identify regions of conserved homozygosity haplotype (RCHHs). Weights were constructed based on the frequency of these RCHHs in case versus controls, and used to adjust the p values from a large collaborative GWAS of COPD. RESULTS: The authors identified 2318 RCHHs, of which 576 were significantly (p<0.05) over-represented in cases. After applying the weights constructed from these regions to a collaborative GWAS of COPD, the authors identified two single nucleotide polymorphisms (SNPs) in a novel gene (fibroblast growth factor-7 (FGF7)) that gained genome-wide significance by the false discovery rate method. In a follow-up analysis, both SNPs (rs12591300 and rs4480740) were significantly associated with COPD in an independent population (combined p values of 7.9E-7 and 2.8E-6, respectively). In another independent population, increased lung tissue FGF7 expression was associated with worse measures of lung function. CONCLUSION: Weights constructed from a homozygosity haplotype analysis of an isolated population successfully identify novel genetic associations from a GWAS on a separate population. This method can be used to identify promising candidate genes that fail to meet strict correction for multiple testing.

Association of SERPINE2 with asthma.

BACKGROUND: The "Dutch hypothesis" suggests that asthma and COPD have common genetic determinants. The serpin peptidase inhibitor, clade E (nexin, plasminogen activator inhibitor type 1), member 2 (SERPINE2) gene previously has been associated with COPD. We sought to determine whether SERPINE2 is associated with asthma and asthma-related phenotypes. METHODS: We measured the association of 39 SERPINE2 single-nucleotide polymorphisms (SNPs) with asthma-related phenotypes in 655 parent-child trios from the Childhood Asthma Management Program (CAMP), and we measured the association of 19 SERPINE2 SNPs with asthma in a case-control design of 359 CAMP probands and 846 population control subjects. We attempted to replicate primary asthma-related phenotype findings in one independent population and primary asthma affection status findings in two independent populations. We compared association results with CAMP proband expression quantitative trait loci. RESULTS: Nine of 39 SNPs had P < .05 for at least one phenotype in CAMP, and two of these replicated in an independent population of 426 people with childhood asthma. Six of 19 SNPs had P < .05 for association with asthma in CAMP/Illumina. None of these replicated in two independent populations. The expression quantitative trait loci revealed that five SNPs associated with asthma in CAMP/Illumina and one SNP associated with FEV(1) in CAMP are strongly correlated with SERPINE2 expression levels. Comparison of results to previous COPD studies identified five SNPs associated with both asthma- and COPD-related phenotypes. CONCLUSIONS: Our results weakly support SERPINE2 as a Dutch hypothesis candidate gene through nominally significant associations with asthma and related traits. Further study of SERPINE2 is necessary to verify its involvement in asthma and COPD.

The impact of self-identified race on epidemiologic studies of gene expression.

Although population differences in gene expression have been established, the impact on differential gene expression studies in large populations is not well understood. We describe the effect of self-reported race on a gene expression study of lung function in asthma. We generated gene expression profiles for 254 young adults (205 non-Hispanic whites and 49 African Americans) with asthma on whom concurrent total RNA derived from peripheral blood CD4(+) lymphocytes and lung function measurements were obtained. We identified four principal components that explained 62% of the variance in gene expression. The dominant principal component, which explained 29% of the total variance in gene expression, was strongly associated with self-identified race (P<10(-16)). The impact of these racial differences was observed when we performed differential gene expression analysis of lung function. Using multivariate linear models, we tested whether gene expression was associated with a quantitative measure of lung function: pre-bronchodilator forced expiratory volume in one second (FEV(1)). Though unadjusted linear models of FEV(1) identified several genes strongly correlated with lung function, these correlations were due to racial differences in the distribution of both FEV(1) and gene expression, and were no longer statistically significant following adjustment for self-identified race. These results suggest that self-identified race is a critical confounding covariate in epidemiologic studies of gene expression and that, similar to genetic studies, careful consideration of self-identified race in gene expression profiling studies is needed to avoid spurious association.