Person:

Tantisira, Kelan

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.

Bronchodilator response (BDR) is an important asthma phenotype that measures reversibility of airway obstruction by comparing lung function (i.e. FEV1) 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 FEV1 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.

Background: Poor maternal vitamin D intake is a risk factor for subsequent childhood asthma, suggesting that in utero changes related to vitamin D responsive genes might play a crucial role in later disease susceptibility. We hypothesized that vitamin D pathway genes are developmentally active in the fetal lung and that these developmental genes would be associated with asthma susceptibility and regulation in asthma. Methods: Vitamin D pathway genes were derived from PubMed and Gene Ontology surveys. Principal component analysis was used to identify characteristic lung development genes. Results: Vitamin D regulated genes were markedly over-represented in normal human (odds ratio OR 2.15, 95% confidence interval CI: 1.69-2.74) and mouse (OR 2.68, 95% CI: 2.12-3.39) developing lung transcriptomes. 38 vitamin D pathway genes were in both developing lung transcriptomes with >63% of genes more highly expressed in the later than earlier stages of development. In immortalized B-cells derived from 95 asthmatics and their unaffected siblings, 12 of the 38 (31.6%) vitamin D pathway lung development genes were significantly differentially expressed (OR 3.00, 95% CI: 1.43-6.21), whereas 11 (29%) genes were significantly differentially expressed in 43 control versus vitamin D treated immortalized B-cells from Childhood Asthma Management Program subjects (OR 2.62, 95% CI: 1.22-5.50). 4 genes, LAMP3, PIP5K1B, SCARB2 and TXNIP were identified in both groups; each displays significant biologic plausibility for a role in asthma. Conclusions: Our findings demonstrate a significant association between early lung development and asthma–related phenotypes for vitamin D pathway genes, supporting a genomic mechanistic basis for the epidemiologic observations relating maternal vitamin D intake and childhood asthma susceptibility.

Background: MicroRNAs are key transcriptional and network regulators previously associated with asthma susceptibility. However, their role in relation to asthma severity has not been delineated. Objective: We hypothesized that circulating microRNAs could serve as biomarkers of changes in lung function in asthma patients. Methods: We isolated microRNAs from serum samples obtained at randomization for 160 participants of the Childhood Asthma Management Program. Using a TaqMan microRNA array containing 754 microRNA primers, we tested for the presence of known asthma microRNAs, and assessed the association of the individual microRNAs with lung function as measured by FEV1/FVC, FEV1% and FVC%. We further tested the subset of FEV1/FVC microRNAs for sex-specific and lung developmental associations. Results: Of the 108 well-detected circulating microRNAs, 74 (68.5%) had previously been linked to asthma susceptibility. We found 22 (20.3%), 4 (3.7%) and 8 (7.4%) microRNAs to be associated with FEV1/FVC, FEV1% and FVC%, respectively. 8 (of 22) FEV1/FVC, 3 (of 4) FEV1% and 1 (of 8) FVC% microRNAs had functionally validated target genes that have been linked via genome wide association studies to asthma and FEV1 change. Among the 22 FEV1/FVC microRNAs, 9 (40.9%) remain associated with FEV1/FVC in boys alone in a sex-stratified analysis (compared with 3 FEV1/FVC microRNAs in girls alone), 7 (31.8%) were associated with fetal lung development, and 3 (13.6%) in both. Ontology analyses revealed enrichment for pathways integral to asthma, including PPAR signaling, G-protein coupled signaling, actin and myosin binding, and respiratory system development. Conclusions: Circulating microRNAs reflect asthma biology and are associated with lung function differences in asthmatics. They may represent biomarkers of asthma severity.

Asthma is a chronic inflammatory respiratory disease that affects over 300 million people worldwide. Glucocorticoids are a mainstay therapy for asthma because they exert anti-inflammatory effects in multiple lung tissues, including the airway smooth muscle (ASM). However, the mechanism by which glucocorticoids suppress inflammation in ASM remains poorly understood. Using RNA-Seq, a high-throughput sequencing method, we characterized transcriptomic changes in four primary human ASM cell lines that were treated with dexamethasone—a potent synthetic glucocorticoid (1 µM for 18 hours). Based on a Benjamini-Hochberg corrected p-value <0.05, we identified 316 differentially expressed genes, including both well known (DUSP1, KLF15, PER1, TSC22D3) and less investigated (C7, CCDC69, CRISPLD2) glucocorticoid-responsive genes. CRISPLD2, which encodes a secreted protein previously implicated in lung development and endotoxin regulation, was found to have SNPs that were moderately associated with inhaled corticosteroid resistance and bronchodilator response among asthma patients in two previously conducted genome-wide association studies. Quantitative RT-PCR and Western blotting showed that dexamethasone treatment significantly increased CRISPLD2 mRNA and protein expression in ASM cells. CRISPLD2 expression was also induced by the inflammatory cytokine IL1β, and small interfering RNA-mediated knockdown of CRISPLD2 further increased IL1β-induced expression of IL6 and IL8. Our findings offer a comprehensive view of the effect of a glucocorticoid on the ASM transcriptome and identify CRISPLD2 as an asthma pharmacogenetics candidate gene that regulates anti-inflammatory effects of glucocorticoids in the ASM.

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.

Abstract 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.

Introduction: Circulating microRNAs (miRNA) are promising biomarkers for human diseases. Our study hypothesizes that circulating miRNA would reveal candidate biomarkers related to airway hyperresponsiveness (AHR) and provide biologic insights into asthma epigenetic influences. Methods: Serum samples obtained at randomization for 160 children in the Childhood Asthma Management Program were profiled using a TaqMan miRNA array set. The association of the isolated miRNA with methacholine PC20 was assessed. Network and pathway analyses were performed. Functional validation of two significant miRNAs was performed in human airway smooth muscle cells (HASMs). Results: Of 155 well-detected circulating miRNAs, eight were significantly associated with PC20 with the strongest association with miR-296-5p. Pathway analysis revealed miR-16-5p as a network hub, and involvement of multiple miRNAs interacting with genes in the FoxO and Hippo signaling pathways by KEGG analysis. Functional validation of two miRNA in HASM showed effects on cell growth and diameter. Conclusion: Reduced circulatory miRNA expression at baseline is associated with an increase in PC20. These miRNA provide biologic insights into, and may serve as biomarkers of, asthma severity. miR-16-5p and -30d-5p regulate airway smooth muscle phenotypes critically involved in asthma pathogenesis, supporting a mechanistic link to these findings. Functional ASM phenotypes may be directly relevant to AHR.