Genomic attributes of airway commensal bacteria and mucosa.

Microbial communities at the airway mucosal barrier are conserved and highly ordered, in likelihood reflecting co-evolution with human host factors. Freed of selection to digest nutrients, the airway microbiome underpins cognate management of mucosal immunity and pathogen resistance. We show here the initial results of systematic culture and whole-genome sequencing of the thoracic airway bacteria, identifying 52 novel species amongst 126 organisms that constitute 75% of commensals typically present in heathy individuals. Clinically relevant genes encode antimicrobial synthesis, adhesion and biofilm formation, immune modulation, iron utilisation, nitrous oxide (NO) metabolism and sphingolipid signalling. Using whole-genome content we identify dysbiotic features that may influence asthma and chronic obstructive pulmonary disease. We match isolate gene content to transcripts and metabolites expressed late in airway epithelial differentiation, identifying pathways to sustain host interactions with microbiota. Our results provide a systematic basis for decrypting interactions between commensals, pathogens, and mucosa in lung diseases of global significance.

Airway microbial communities, smoking and asthma in a general population sample.

BACKGROUND: Normal airway microbial communities play a central role in respiratory health but are poorly characterized. Cigarette smoking is the dominant global environmental influence on lung function, and asthma has become the most prevalent chronic respiratory disease worldwide. Both conditions have major microbial components that are incompletely defined. METHODS: We investigated airway bacterial communities in a general population sample of 529 Australian adults. Posterior oropharyngeal swabs were analyzed by sequencing of the 16S rRNA gene. The microbiota were characterized according to their prevalence, abundance and network memberships. FINDINGS: The microbiota were similar across the general population, and were strongly organized into co-abundance networks. Smoking was associated with diversity loss, negative effects on abundant taxa, profound alterations to network structure and expansion of Streptococcus spp. By contrast, the asthmatic microbiota were selectively affected by an increase in Neisseria spp. and by reduced numbers of low abundance but prevalent organisms. INTERPRETATION: Our study shows that the healthy airway microbiota in this population were contained within a highly structured ecosystem, suggesting balanced relationships between the microbiome and human host factors. The marked abnormalities in smokers may contribute to chronic obstructive pulmonary disease (COPD) and lung cancer. The narrow spectrum of abnormalities in asthmatics encourages investigation of damaging and protective effects of specific bacteria. FUNDING: The study was funded by the Asmarley Trust and a Wellcome Joint Senior Investigator Award to WOCC and MFM (WT096964MA and WT097117MA). The Busselton Healthy Ageing Study is supported by the Government of Western Australia (Office of Science, Department of Health) the City of Busselton, and private donations.

Y disruption, autosomal hypomethylation and poor male lung cancer survival.

Lung cancer is the most frequent cause of cancer death worldwide. It affects more men than women, and men generally have worse survival outcomes. We compared gene co-expression networks in affected and unaffected lung tissue from 126 consecutive patients with Stage IA-IV lung cancer undergoing surgery with curative intent. We observed marked degradation of a sex-associated transcription network in tumour tissue. This disturbance, detected in 27.7% of male tumours in the discovery dataset and 27.3% of male tumours in a further 123-sample replication dataset, was coincident with partial losses of the Y chromosome and extensive autosomal DNA hypomethylation. Central to this network was the epigenetic modifier and regulator of sexually dimorphic gene expression, KDM5D. After accounting for prognostic and epidemiological covariates including stage and histology, male patients with tumour KDM5D deficiency showed a significantly increased risk of death (Hazard Ratio [HR] 3.80, 95% CI 1.40-10.3, P = 0.009). KDM5D deficiency was confirmed as a negative prognostic indicator in a further 1100 male lung tumours (HR 1.67, 95% CI 1.4-2.0, P = 1.2 × 10-10). Our findings identify tumour deficiency of KDM5D as a prognostic marker and credible mechanism underlying sex disparity in lung cancer.

ORMDL3 regulates poly I:C induced inflammatory responses in airway epithelial cells.

BACKGROUND: Oroscomucoid 3 (ORMDL3) has been linked to susceptibility of childhood asthma and respiratory viral infection. Polyinosinic-polycytidylic acid (poly I:C) is a synthetic analog of viral double-stranded RNA, a toll-like receptor 3 (TLR3) ligand and mimic of viral infection. METHODS: To investigate the functional role of ORMDL3 in the poly I:C-induced inflammatory response in airway epithelial cells, ORMDL3 knockdown and over-expression models were established in human A549 epithelial cells and primary normal human bronchial epithelial (NHBE) cells. The cells were stimulated with poly I:C or the Th17 cytokine IL-17A. IL-6 and IL-8 levels in supernatants,  mRNA levels of genes in the TLR3 pathway and inflammatory response from cell pellets were measured. ORMDL3 knockdown models in A549 and BEAS-2B epithelial cells were then infected with live human rhinovirus (HRV16) followed by IL-6 and IL-8 measurement. RESULTS: ORMDL3 knockdown and over-expression had little influence on the transcript levels of TLR3 in airway epithelial cells. Time course studies showed that ORMDL3-deficient A549 and NHBE cells had an attenuated IL-6 and IL-8 response to poly I:C stimulation. A549 and NHBE cells over-expressing ORMDL3 released relatively more IL-6 and IL-8 following poly I:C stimulation. IL-17A exhibited a similar inflammatory response in ORMDL3 knockdown and over-expressing cells, but co-stimulation of poly I:C and IL-17A did not significantly enhance the IL-6 and IL-8 response. Transcript abundance of IFNB following poly I:C stimulation was not significantly altered by ORMDL3 knockdown or over-expression. Dampening of the IL-6 response by ORMDL3 knockdown was confirmed in HRV16 infected BEAS-2B and A549 cells. CONCLUSIONS: ORMDL3 regulates the viral inflammatory response in airway epithelial cells via mechanisms independent of the TLR3 pathway.

Metabolomic, transcriptomic and genetic integrative analysis reveals important roles of adenosine diphosphate in haemostasis and platelet activation in non-small-cell lung cancer.

Lung cancer is the leading cause of cancer-related deaths in the world. The most prevalent subtype, accounting for 85% of cases, is non-small-cell lung cancer (NSCLC). Lung squamous cell carcinoma (LUSC) and lung adenocarcinoma (LUAD) are the most common subtypes. Despite recent advances in treatment, the low 5-year survival rate of NSCLC patients (approximately 13%) reflects the lack of early diagnostic biomarkers and incomplete understanding of the underlying disease mechanisms. We hypothesized that integration of metabolomic, transcriptomic and genetic profiles of tumours and matched normal tissues could help to identify important factors and potential therapeutic targets that contribute to tumorigenesis. We integrated omics profiles in tumours and matched adjacent normal tissues of patients with LUSC (N = 20) and LUAD (N = 17) using multiple system biology approaches. We confirmed the presence of previously described metabolic pathways in NSCLC, particularly those mediating the Warburg effect. In addition, through our combined omics analyses we found that metabolites and genes that contribute to haemostasis, angiogenesis, platelet activation and cell proliferation were predominant in both subtypes of NSCLC. The important roles of adenosine diphosphate in promoting cancer metastasis through platelet activation and angiogenesis suggest this metabolite could be a potential therapeutic target.

The ORMDL3 Asthma Gene Regulates ICAM1 and Has Multiple Effects on Cellular Inflammation.

RATIONALE: Polymorphisms on chromosome 17q21 confer the major genetic susceptibility to childhood-onset asthma. Risk alleles positively correlate with ORMDL3 (orosomucoid-like 3) expression. The locus influences disease severity and the frequency of human rhinovirus (HRV)-initiated exacerbations. ORMDL3 is known to regulate sphingolipid synthesis by binding serine palmitoyltransferase, but its role in inflammation is incompletely understood. OBJECTIVES: To investigate the role of ORMDL3 in cellular inflammation. METHODS: We modeled a time series of IL1B-induced inflammation in A549 cells, using cytokine production as outputs and testing effects of ORMDL3 siRNA knockdown, ORMDL3 overexpression, and the serine palmitoyltransferase inhibitor myriocin. We replicated selected findings in normal human bronchial epithelial cells. Cytokine and metabolite levels were analyzed by analysis of variance. Transcript abundances were analyzed by group means parameterization, controlling the false discovery rate below 0.05. MEASUREMENTS AND MAIN RESULTS: Silencing ORMDL3 led to steroid-independent reduction of IL6 and IL8 release and reduced endoplasmic reticulum stress after IL1B stimulation. Overexpression and myriocin conversely augmented cytokine release. Knockdown reduced expression of genes regulating host-pathogen interactions, stress responses, and ubiquitination: in particular, ORMDL3 knockdown strongly reduced expression of the HRV receptor ICAM1. Silencing led to changes in levels of transcripts and metabolites integral to glycolysis. Increased levels of ceramides and the immune mediator sphingosine-1-phosphate were also observed. CONCLUSIONS: The results show ORMDL3 has pleiotropic effects during cellular inflammation, consistent with its substantial genetic influence on childhood asthma. Actions on ICAM1 provide a mechanism for the locus to confer susceptibility to HRV-induced asthma.

The Genetics and Genomics of Asthma.

Asthma is a common, clinically heterogeneous disease with strong evidence of heritability. Progress in defining the genetic underpinnings of asthma, however, has been slow and hampered by issues of inconsistency. Recent advances in the tools available for analysis-assaying transcription, sequence variation, and epigenetic marks on a genome-wide scale-have substantially altered this landscape. Applications of such approaches are consistent with heterogeneity at the level of causation and specify patterns of commonality with a wide range of alternative disease traits. Looking beyond the individual as the unit of study, advances in technology have also fostered comprehensive analysis of the human microbiome and its varied roles in health and disease. In this article, we consider the implications of these technological advances for our current understanding of the genetics and genomics of asthma.

COPD is accompanied by co-ordinated transcriptional perturbation in the quadriceps affecting the mitochondria and extracellular matrix.

Skeletal muscle dysfunction is a frequent extra-pulmonary manifestation of Chronic Obstructive Pulmonary Disease (COPD) with implications for both quality of life and survival. The underlying biology nevertheless remains poorly understood. We measured global gene transcription in the quadriceps using Affymetrix HuGene1.1ST arrays in an unselected cohort of 79 stable COPD patients in secondary care and 16 healthy age- and gender-matched controls. We detected 1,826 transcripts showing COPD-related variation. Eighteen exhibited ≥2fold changes (SLC22A3, FAM184B, CDKN1A, FST, LINC01405, MUSK, PANX1, ANKRD1, C12orf75, MYH1, POSTN, FRZB, TNC, ACTC1, LINC00310, MYH3, MYBPH and AREG). Thirty-one transcripts possessed previous reported evidence of involvement in COPD through genome-wide association, including FAM13A. Network analysis revealed a substructure comprising 6 modules of co-expressed genes. We identified modules with mitochondrial and extracellular matrix features, of which IDH2, a central component of the mitochondrial antioxidant pathway, and ABI3BP, a proposed switch between proliferation and differentiation, represent hubs respectively. COPD is accompanied by coordinated patterns of transcription in the quadriceps involving the mitochondria and extracellular matrix and including genes previously implicated in primary disease processes.

Host-Microbial Interactions in Idiopathic Pulmonary Fibrosis.

RATIONALE: Changes in the respiratory microbiome are associated with disease progression in idiopathic pulmonary fibrosis (IPF). The role of the host response to the respiratory microbiome remains unknown. OBJECTIVES: To explore the host-microbial interactions in IPF. METHODS: Sixty patients diagnosed with IPF were prospectively enrolled together with 20 matched control subjects. Subjects underwent bronchoalveolar lavage (BAL), and peripheral whole blood was collected into PAXgene tubes for all subjects at baseline. For subjects with IPF, additional samples were taken at 1, 3, and 6 months and (if alive) 1 year. Gene expression profiles were generated using Affymetrix Human Gene 1.1 ST arrays. MEASUREMENTS AND MAIN RESULTS: By network analysis of gene expression data, we identified two gene modules that strongly associated with a diagnosis of IPF, BAL bacterial burden (determined by 16S quantitative polymerase chain reaction), and specific microbial operational taxonomic units, as well as with lavage and peripheral blood neutrophilia. Genes within these modules that are involved in the host defense response include NLRC4, PGLYRP1, MMP9, and DEFA4. The modules also contain two genes encoding specific antimicrobial peptides (SLPI and CAMP). Many of these particular transcripts were associated with survival and showed longitudinal overexpression in subjects experiencing disease progression, further strengthening the relationship of the transcripts with disease. CONCLUSIONS: Integrated analysis of the host transcriptome and microbial signatures demonstrated an apparent host response to the presence of an altered or more abundant microbiome. These responses remained elevated in longitudinal follow-up, suggesting that the bacterial communities of the lower airways may act as persistent stimuli for repetitive alveolar injury in IPF.

A mechanistic target of rapamycin complex 1/2 (mTORC1)/V-Akt murine thymoma viral oncogene homolog 1 (AKT1)/cathepsin H axis controls filaggrin expression and processing in skin, a novel mechanism for skin barrier disruption in patients with atopic dermatitis.

BACKGROUND: Filaggrin, which is encoded by the filaggrin gene (FLG), is an important component of the skin's barrier to the external environment, and genetic defects in FLG strongly associate with atopic dermatitis (AD). However, not all patients with AD have FLG mutations. OBJECTIVE: We hypothesized that these patients might possess other defects in filaggrin expression and processing contributing to barrier disruption and AD, and therefore we present novel therapeutic targets for this disease. RESULTS: We describe the relationship between the mechanistic target of rapamycin complex 1/2 protein subunit regulatory associated protein of the MTOR complex 1 (RAPTOR), the serine/threonine kinase V-Akt murine thymoma viral oncogene homolog 1 (AKT1), and the protease cathepsin H (CTSH), for which we establish a role in filaggrin expression and processing. Increased RAPTOR levels correlated with decreased filaggrin expression in patients with AD. In keratinocyte cell cultures RAPTOR upregulation or AKT1 short hairpin RNA knockdown reduced expression of the protease CTSH. Skin of CTSH-deficient mice and CTSH short hairpin RNA knockdown keratinocytes showed reduced filaggrin processing, and the mouse had both impaired skin barrier function and a mild proinflammatory phenotype. CONCLUSION: Our findings highlight a novel and potentially treatable signaling axis controlling filaggrin expression and processing that is defective in patients with AD.

Global gene regulation during activation of immunoglobulin class switching in human B cells.

Immunoglobulin class switch recombination (CSR) to IgE is a tightly regulated process central to atopic disease. To profile the B-cell transcriptional responses underlying the activation of the germinal centre activities leading to the generation of IgE, naïve human B-cells were stimulated with IL-4 and anti-CD40. Gene expression and alternative splicing were profiled over 12 days using the Affymetrix Human Exon 1.0 ST Array. A total of 1,399 genes, forming 13 temporal profiles were differentially expressed. CCL22 and CCL17 were dramatically induced but followed a temporal trajectory distinct from classical mediators of isotype switching. AICDA, NFIL3, IRF4, XBP1 and BATF3 shared a profile with several genes involved in innate immunity, but with no recognised role in CSR. A transcription factor BHLHE40 was identified at the core of this profile. B-cell activation was also accompanied by variation in exon retention affecting >200 genes including CCL17. The data indicate a circadian component and central roles for the Th2 chemokines CCL22 and CCL17 in the activation of CSR.

An epigenome-wide association study of total serum immunoglobulin E concentration.

Immunoglobulin E (IgE) is a central mediator of allergic (atopic) inflammation. Therapies directed against IgE can alleviate hay fever and allergic asthma. Genetic association studies have not yet identified novel therapeutic targets or pathways underlying IgE regulation. We therefore surveyed epigenetic associations between serum IgE concentrations and methylation at loci concentrated in CpG islands genome wide in 95 nuclear pedigrees, using DNA from peripheral blood leukocytes. We validated positive results in additional families and in subjects from the general population. Here we show replicated associations--with a meta-analysis false discovery rate less than 10(-4)--between IgE and low methylation at 36 loci. Genes annotated to these loci encode known eosinophil products, and also implicate phospholipid inflammatory mediators, specific transcription factors and mitochondrial proteins. We confirmed that methylation at these loci differed significantly in isolated eosinophils from subjects with and without asthma and high IgE levels. The top three loci accounted for 13% of IgE variation in the primary subject panel, explaining the tenfold higher variance found compared with that derived from large single-nucleotide polymorphism genome-wide association studies. This study identifies novel therapeutic targets and biomarkers for patient stratification for allergic diseases.

The role of bacteria in the pathogenesis and progression of idiopathic pulmonary fibrosis.

RATIONALE: Idiopathic pulmonary fibrosis (IPF) is a progressive lung disease of unknown cause that leads to respiratory failure and death within 5 years of diagnosis. Overt respiratory infection and immunosuppression carry a high morbidity and mortality, and polymorphisms in genes related to epithelial integrity and host defense predispose to IPF. OBJECTIVES: To investigate the role of bacteria in the pathogenesis and progression of IPF. METHODS: We prospectively enrolled patients diagnosed with IPF according to international criteria together with healthy smokers, nonsmokers, and subjects with moderate chronic obstructive pulmonary disease as control subjects. Subjects underwent bronchoalveolar lavage (BAL), from which genomic DNA was isolated. The V3-V5 region of the bacterial 16S rRNA gene was amplified, allowing quantification of bacterial load and identification of communities by 16S rRNA quantitative polymerase chain reaction and pyrosequencing. MEASUREMENTS AND MAIN RESULTS: Sixty-five patients with IPF had double the burden of bacteria in BAL fluid compared with 44 control subjects. Baseline bacterial burden predicted the rate of decline in lung volume and risk of death and associated independently with the rs35705950 polymorphism of the MUC5B mucin gene, a proven host susceptibility factor for IPF. Sequencing yielded 912,883 high-quality reads from all subjects. We identified Haemophilus, Streptococcus, Neisseria, and Veillonella spp. to be more abundant in cases than control subjects. Regression analyses indicated that these specific operational taxonomic units as well as bacterial burden associated independently with IPF. CONCLUSIONS: IPF is characterized by an increased bacterial burden in BAL that predicts decline in lung function and death. Trials of antimicrobial therapy are needed to determine if microbial burden is pathogenic in the disease.

Predicting DNA methylation level across human tissues.

Differences in methylation across tissues are critical to cell differentiation and are key to understanding the role of epigenetics in complex diseases. In this investigation, we found that locus-specific methylation differences between tissues are highly consistent across individuals. We developed a novel statistical model to predict locus-specific methylation in target tissue based on methylation in surrogate tissue. The method was evaluated in publicly available data and in two studies using the latest IlluminaBeadChips: a childhood asthma study with methylation measured in both peripheral blood leukocytes (PBL) and lymphoblastoid cell lines; and a study of postoperative atrial fibrillation with methylation in PBL, atrium and artery. We found that our method can greatly improve accuracy of cross-tissue prediction at CpG sites that are variable in the target tissue [R(2) increases from 0.38 (original R(2) between tissues) to 0.89 for PBL-to-artery prediction; from 0.39 to 0.95 for PBL-to-atrium; and from 0.81 to 0.98 for lymphoblastoid cell line-to-PBL based on cross-validation, and confirmed using cross-study prediction]. An extended model with multiple CpGs further improved performance. Our results suggest that large-scale epidemiology studies using easy-to-access surrogate tissues (e.g. blood) could be recalibrated to improve understanding of epigenetics in hard-to-access tissues (e.g. atrium) and might enable non-invasive disease screening using epigenetic profiles.

Outgrowth of the bacterial airway microbiome after rhinovirus exacerbation of chronic obstructive pulmonary disease.

RATIONALE: Rhinovirus infection is followed by significantly increased frequencies of positive, potentially pathogenic sputum cultures in chronic obstructive pulmonary disease (COPD). However, it remains unclear whether these represent de novo infections or an increased load of organisms from the complex microbial communities (microbiome) in the lower airways. OBJECTIVES: To investigate the effect of rhinovirus infection on the airway bacterial microbiome. METHODS: Subjects with COPD (n = 14) and healthy control subjects with normal lung function (n = 17) were infected with rhinovirus. Induced sputum was collected at baseline before rhinovirus inoculation and again on Days 5, 15, and 42 after rhinovirus infection and DNA was extracted. The V3-V5 region of the bacterial 16S ribosomal RNA gene was amplified and pyrosequenced, resulting in 370,849 high-quality reads from 112 of the possible 124 time points. MEASUREMENTS AND MAIN RESULTS: At 15 days after rhinovirus infection, there was a sixfold increase in 16S copy number (P = 0.007) and a 16% rise in numbers of proteobacterial sequences, most notably in potentially pathogenic Haemophilus influenzae (P = 2.7 × 10(-20)), from a preexisting community. These changes occurred only in the sputum microbiome of subjects with COPD and were still evident 42 days after infection. This was in contrast to the temporal stability demonstrated in the microbiome of healthy smokers and nonsmokers. CONCLUSIONS: After rhinovirus infection, there is a rise in bacterial burden and a significant outgrowth of Haemophilus influenzae from the existing microbiota of subjects with COPD. This is not observed in healthy individuals. Our findings suggest that rhinovirus infection in COPD alters the respiratory microbiome and may precipitate secondary bacterial infections.

A genome-wide association study of atopic dermatitis identifies loci with overlapping effects on asthma and psoriasis.

Atopic dermatitis (AD) is the most common dermatological disease of childhood. Many children with AD have asthma and AD shares regions of genetic linkage with psoriasis, another chronic inflammatory skin disease. We present here a genome-wide association study (GWAS) of childhood-onset AD in 1563 European cases with known asthma status and 4054 European controls. Using Illumina genotyping followed by imputation, we generated 268 034 consensus genotypes and in excess of 2 million single nucleotide polymorphisms (SNPs) for analysis. Association signals were assessed for replication in a second panel of 2286 European cases and 3160 European controls. Four loci achieved genome-wide significance for AD and replicated consistently across all cohorts. These included the epidermal differentiation complex (EDC) on chromosome 1, the genomic region proximal to LRRC32 on chromosome 11, the RAD50/IL13 locus on chromosome 5 and the major histocompatibility complex (MHC) on chromosome 6; reflecting action of classical HLA alleles. We observed variation in the contribution towards co-morbid asthma for these regions of association. We further explored the genetic relationship between AD, asthma and psoriasis by examining previously identified susceptibility SNPs for these diseases. We found considerable overlap between AD and psoriasis together with variable coincidence between allergic rhinitis (AR) and asthma. Our results indicate that the pathogenesis of AD incorporates immune and epidermal barrier defects with combinations of specific and overlapping effects at individual loci.

Variants of DENND1B associated with asthma in children.

BACKGROUND: Asthma is a complex disease that has genetic and environmental causes. The genetic factors associated with susceptibility to asthma remain largely unknown. METHODS: We carried out a genomewide association study involving children with asthma. The sample included 793 North American children of European ancestry with persistent asthma who required daily inhaled glucocorticoid therapy and 1988 matched controls (the discovery set). We also tested for genomewide association in an independent cohort of 917 persons of European ancestry who had asthma and 1546 matched controls (the replication set). Finally, we tested for an association between 20 single-nucleotide polymorphisms (SNPs) at chromosome 1q31 and asthma in 1667 North American children of African ancestry who had asthma and 2045 ancestrally matched controls. RESULTS: In our meta-analysis of all samples from persons of European ancestry, we observed an association, with genomewide significance, between asthma and SNPs at the previously reported locus on 17q21 and an additional eight SNPs at a novel locus on 1q31. The SNP most strongly associated with asthma was rs2786098 (P=8.55x10(-9)). We observed replication of the association of asthma with SNP rs2786098 in the independent series of persons of European ancestry (combined P=9.3x10(-11)). The alternative allele of each of the eight SNPs on chromosome 1q31 was strongly associated with asthma in the children of African ancestry (P=1.6x10(-13) for the comparison across all samples). The 1q31 locus contains the 1q31 locus contains DENND1B, a gene expressed by natural killer cells and dendritic cells. DENND1B protein is predicted to interact with the tumor necrosis factor α receptor [corrected]. CONCLUSIONS: We have identified a locus containing DENND1B on chromosome 1q31.3 that is associated with susceptibility to asthma.

Genome-wide association study implicates chromosome 9q21.31 as a susceptibility locus for asthma in mexican children.

Many candidate genes have been studied for asthma, but replication has varied. Novel candidate genes have been identified for various complex diseases using genome-wide association studies (GWASs). We conducted a GWAS in 492 Mexican children with asthma, predominantly atopic by skin prick test, and their parents using the Illumina HumanHap 550 K BeadChip to identify novel genetic variation for childhood asthma. The 520,767 autosomal single nucleotide polymorphisms (SNPs) passing quality control were tested for association with childhood asthma using log-linear regression with a log-additive risk model. Eleven of the most significantly associated GWAS SNPs were tested for replication in an independent study of 177 Mexican case-parent trios with childhood-onset asthma and atopy using log-linear analysis. The chromosome 9q21.31 SNP rs2378383 (p = 7.10x10(-6) in the GWAS), located upstream of transducin-like enhancer of split 4 (TLE4), gave a p-value of 0.03 and the same direction and magnitude of association in the replication study (combined p = 6.79x10(-7)). Ancestry analysis on chromosome 9q supported an inverse association between the rs2378383 minor allele (G) and childhood asthma. This work identifies chromosome 9q21.31 as a novel susceptibility locus for childhood asthma in Mexicans. Further, analysis of genome-wide expression data in 51 human tissues from the Novartis Research Foundation showed that median GWAS significance levels for SNPs in genes expressed in the lung differed most significantly from genes not expressed in the lung when compared to 50 other tissues, supporting the biological plausibility of our overall GWAS findings and the multigenic etiology of childhood asthma.

Genome-wide association analysis identifies PDE4D as an asthma-susceptibility gene.

Asthma, a chronic airway disease with known heritability, affects more than 300 million people around the world. A genome-wide association (GWA) study of asthma with 359 cases from the Childhood Asthma Management Program (CAMP) and 846 genetically matched controls from the Illumina ICONdb public resource was performed. The strongest region of association seen was on chromosome 5q12 in PDE4D. The phosphodiesterase 4D, cAMP-specific (phosphodiesterase E3 dunce homolog, Drosophila) gene (PDE4D) is a regulator of airway smooth-muscle contractility, and PDE4 inhibitors have been developed as medications for asthma. Allelic p values for top SNPs in this region were 4.3 x 10(-07) for rs1588265 and 9.7 x 10(-07) for rs1544791. Replications were investigated in ten independent populations with different ethnicities, study designs, and definitions of asthma. In seven white and Hispanic replication populations, two PDE4D SNPs had significant results with p values less than 0.05, and five had results in the same direction as the original population but had p values greater than 0.05. Combined p values for 18,891 white and Hispanic individuals (4,342 cases) in our replication populations were 4.1 x 10(-04) for rs1588265 and 9.2 x 10(-04) for rs1544791. In three black replication populations, which had different linkage disequilibrium patterns than the other populations, original findings were not replicated. Further study of PDE4D variants might lead to improved understanding of the role of PDE4D in asthma pathophysiology and the efficacy of PDE4 inhibitor medications.