Person:

Litonjua, Augusto A.

Background: The mechanisms for the association between birth by cesarean section and atopy and asthma are largely unknown. Objective: To examine whether cesarean section results in neonatal secretion of cytokines that are associated with increased risk of atopy and/or asthma in childhood. To examine whether the association between mode of delivery and neonatal immune responses is explained by exposure to the maternal gut flora (a marker of the vaginal flora). Methods: CBMCs were isolated from 37 neonates at delivery, and secretion of IL-13, IFN-γ, and IL-10 (at baseline and after stimulation with antigens [dust mite and cat dander allergens, phytohemagglutinin, and lipopolysaccharide]) was quantified by ELISA. Total and specific microbes were quantified in maternal stool. The relation between mode of delivery and cord blood cytokines was examined by linear regression. The relation between maternal stool microbes and cord blood cytokines was examined by Spearman's correlation coefficients. Results: Cesarean section was associated with increased levels of IL-13 and IFN-γ. In multivariate analyses, cesarean section was associated with an increment of 79.4 pg/ml in secretion of IL-13 by CBMCs after stimulation with dust mite allergen (P < 0.001). Among children born by vaginal delivery, gram-positive anaerobes and total anaerobes in maternal stool were positively correlated with levels of IL-10, and gram-negative aerobic bacteria in maternal stool were negatively correlated with levels of IL-13 and IFN-γ. Conclusion: Cesarean section is associated with increased levels of IL-13 and IFN-γ, perhaps because of lack of labor and/or reduced exposure to specific microbes (e.g., gram-positive anaerobes) at birth.

Background: Identifying genetic determinants for lung function is important in providing insight into the pathophysiology of asthma. Signal transducer and activator of transcription 3 is a transcription factor latent in the cytoplasm; the gene (STAT3) is activated by a wide range of cytokines, and may play a role in lung development and asthma pathogenesis. Methods: We genotyped six single nucleotide polymorphisms (SNPs) in the STAT3 gene in a cohort of 401 Caucasian adult asthmatics. The associations between each SNP and forced expiratory volume in 1 second (FEV({1})), as a percent of predicted, at the baseline exam were tested using multiple linear regression models. Longitudinal analyses involving repeated measures of FEV({1}) were conducted with mixed linear models. Haplotype analyses were conducted using imputed haplotypes. We completed a second association study by genotyping the same six polymorphisms in a cohort of 652 Caucasian children with asthma. Results: We found that three polymorphisms were significantly associated with baseline FEV({1}): homozygotes for the minor alleles of each polymorphism had lower FEV({1}) than homozygotes for the major alleles. Moreover, these associations persisted when we performed an analysis on repeated measures of FEV({1}) over 8 weeks. A haplotypic analysis based on the six polymorphisms indicated that two haplotypes were associated with baseline FEV({1}). Among the childhood asthmatics, one polymorphism was associated with both baseline FEV({1}) and the repeated measures of FEV({1}) over 4 years. Conclusion: Our results indicate that genetic variants in STAT3, independent of asthma treatment, are determinants of FEV(_{1}) in both adults and children with asthma, and suggest that STAT3 may participate in inflammatory pathways that have an impact on level of lung function.

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: Sequence variants in genes functioning in folate-mediated one-carbon metabolism are hypothesized to lead to changes in levels of homocysteine and DNA methylation, which, in turn, are associated with risk of cardiovascular disease. Methods: 330 SNPs in 52 genes were studied in relation to plasma homocysteine and global genomic DNA methylation. SNPs were selected based on functional effects and gene coverage, and assays were completed on the Illumina Goldengate platform. Age-, smoking-, and nutrient-adjusted genotype--phenotype associations were estimated in regression models. Results: Using a nominal P (\leq) 0.005 threshold for statistical significance, 20 SNPs were associated with plasma homocysteine, 8 with Alu methylation, and 1 with LINE-1 methylation. Using a more stringent false discovery rate threshold, SNPs in FTCD, SLC19A1, and SLC19A3 genes remained associated with plasma homocysteine. Gene by vitamin B-6 interactions were identified for both Alu and LINE-1 methylation, and epistatic interactions with the MTHFR rs1801133 SNP were identified for the plasma homocysteine phenotype. Pleiotropy involving the MTHFD1L and SARDH genes for both plasma homocysteine and Alu methylation phenotypes was identified. Conclusions: No single gene was associated with all three phenotypes, and the set of the most statistically significant SNPs predictive of homocysteine or Alu or LINE-1 methylation was unique to each phenotype. Genetic variation in folate-mediated one-carbon metabolism, other than the well-known effects of the MTHFR c.665C>T (known as c.677 C>T, rs1801133, p.Ala222Val), is predictive of cardiovascular disease biomarkers.

Globally, asthma is a chronic inflammatory respiratory disease affecting over 300 million people. Some asthma patients remain poorly controlled by conventional therapies and experience more life-threatening exacerbations. Vitamin D, as an adjunct therapy, may improve disease control in severe asthma patients since vitamin D enhances glucocorticoid responsiveness and mitigates airway smooth muscle (ASM) hyperplasia. We sought to characterize differences in transcriptome responsiveness to vitamin D between fatal asthma- and non-asthma-derived ASM by using RNA-Seq to measure ASM transcript expression in five donors with fatal asthma and ten non-asthma-derived donors at baseline and with vitamin D treatment. Based on a Benjamini-Hochberg corrected p-value <0.05, 838 genes were differentially expressed in fatal asthma vs. non-asthma-derived ASM at baseline, and vitamin D treatment compared to baseline conditions induced differential expression of 711 and 867 genes in fatal asthma- and non-asthma-derived ASM, respectively. Functional gene categories that were highly represented in all groups included extracellular matrix, and responses to steroid hormone stimuli and wounding. Genes differentially expressed by vitamin D also included cytokine and chemokine activity categories. Follow-up qPCR and individual analyte ELISA experiments were conducted for four cytokines (i.e. CCL2, CCL13, CXCL12, IL8) to measure TNFα-induced changes by asthma status and vitamin D treatment. Vitamin D inhibited TNFα-induced IL8 protein secretion levels to a comparable degree in fatal asthma- and non-asthma-derived ASM even though IL8 had significantly higher baseline levels in fatal asthma-derived ASM. Our findings identify vitamin D-specific gene targets and provide transcriptomic data to explore differences in the ASM of fatal asthma- and non-asthma-derived donors.

Background: The response to treatment for asthma is characterized by wide interindividual variability, with a significant number of patients who have no response. We hypothesized that a genomewide association study would reveal novel pharmacogenetic determinants of the response to inhaled glucocorticoids. Methods: We analyzed a small number of statistically powerful variants selected on the basis of a family-based screening algorithm from among 534,290 single-nucleotide polymorphisms (SNPs) to determine changes in lung function in response to inhaled glucocorticoids. A significant, replicated association was found, and we characterized its functional effects. Results: We identified a significant pharmacogenetic association at SNP rs37972, replicated in four independent populations totaling 935 persons (P=0.0007), which maps to the glucocorticoid-induced transcript 1 gene (GLCCI1) and is in complete linkage disequilibrium (i.e., perfectly correlated) with rs37973. Both rs37972 and rs37973 are associated with decrements in GLCCI1 expression. In isolated cell systems, the rs37973 variant is associated with significantly decreased luciferase reporter activity. Pooled data from treatment trials indicate reduced lung function in response to inhaled glucocorticoids in subjects with the variant allele (P=0.0007 for pooled data). Overall, the mean (±SE) increase in forced expiratory volume in 1 second in the treated subjects who were homozygous for the mutant rs37973 allele was only about one third of that seen in similarly treated subjects who were homozygous for the wild-type allele (3.2±1.6% vs. 9.4±1.1%), and their risk of a poor response was significantly higher (odds ratio, 2.36; 95% confidence interval, 1.27 to 4.41), with genotype accounting for about 6.6% of overall inhaled glucocorticoid response variability. Conclusions: A functional GLCCI1 variant is associated with substantial decrements in the response to inhaled glucocorticoids in patients with asthma.

Background: A modest number of prospective studies of the composition of the intestinal microbiota and eczema in early life have yielded conflicting results. Objective: To examine the relationship between the bacterial diversity of the gut and the development of eczema in early life by methods other than stool culture. Methods: Fecal samples were collected from 21 infants at 1 and 4 months of life. Nine infants were diagnosed with eczema by the age of 6 months (cases) and 12 infants were not (controls). After conducting denaturating gradient gel electrophoresis (DGGE) of stool samples, we compared the microbial diversity of cases and controls using the number of electrophoretic bands and the Shannon index of diversity (H') as indicators. Results: Control subjects had significantly greater fecal microbial diversity than children with eczema at ages 1 (mean H' for controls = 0.75 vs. 0.53 for cases, P = 0.01) and 4 months (mean H' for controls = 0.92 vs. 0.59 for cases, P = 0.02). The increase in diversity from 1 to 4 months of age was significant in controls (P = 0.04) but not in children who developed eczema by 6 months of age (P = 0.32). Conclusion: Our findings suggest that reduced microbial diversity is associated with the development of eczema in early life.

Background— It has recently been shown that vitamin D deficiency can increase asthma development and severity and that variations in vitamin D receptor genes are associated with asthma susceptibility. Objective— We sought to find genetic factors that might interact with vitamin D levels to affect the risk of asthma exacerbation. Methods: We conducted a genome-wide study of gene–vitamin D interaction on asthma exacerbations using population-based and family-based approaches on 403 subjects and trios from the Childhood Asthma Management Program. Twenty-three polymorphisms with significant interactions were studied in a replication analysis in 584 children from a Costa Rican cohort. Results: We identified 3 common variants in the class I MHC– restricted T cell–associated molecule gene (CRTAM) that were associated with an increased rate of asthma exacerbations based on the presence of a low circulating vitamin D level. These results were replicated in a second independent population (unadjusted combined interaction, P =.00028–. 00097; combined odds ratio, 3.28–5.38). One variant, rs2272094, is a nonsynonymous coding polymorphism of CRTAM. Functional studies on cell lines confirmed the interaction of vitamin D and rs2272094 on CRTAM expression. CRTAM is highly expressed in activated human CD8 + and natural killer T cells, both of which have been implicated in asthmatic patients. Conclusion— The findings highlight an important gene-environment interaction that elucidates the role of vitamin D and CD8 + and natural killer T cells in asthma exacerbation in a genome-wide gene-environment interaction study that has been replicated in an independent population. The results suggest the potential importance of maintaining adequate vitamin D levels in subsets of high-risk asthmatic patients.