Person:

Carey, Vincent

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.

Recent breakthroughs in next-generation sequencing technologies allow cost-effective methods for measuring a growing list of cellular properties, including DNA sequence and structural variation. Next-generation sequencing has the potential to revolutionize complex trait genetics by directly measuring common and rare genetic variants within a genome-wide context. Because for a given gene both rare and common causal variants can coexist and have independent effects on a trait, strategies that model the effects of both common and rare variants could enhance the power of identifying disease-associated genes. To date, little work has been done on integrating signals from common and rare variants into powerful statistics for finding disease genes in genome-wide association studies. In this analysis of the Genetic Analysis Workshop 17 data, we evaluate various strategies for association of rare, common, or a combination of both rare and common variants on quantitative phenotypes in unrelated individuals. We show that the analysis of common variants only using classical approaches can achieve higher power to detect causal genes than recently proposed rare variant methods and that strategies that combine association signals derived independently in rare and common variants can slightly increase the power compared to strategies that focus on the effect of either the rare variants or the common variants.

Background: Smoking while pregnant is associated with a myriad of negative health outcomes in the child. Some of the detrimental effects may be due to epigenetic modifications, although few studies have investigated this hypothesis in detail. Objectives: To characterize site-specific epigenetic modifications conferred by prenatal smoking exposure within asthmatic children. Methods: Using Illumina HumanMethylation27 microarrays, we estimated the degree of methylation at 27,578 distinct DNA sequences located primarily in gene promoters using whole blood DNA samples from the Childhood Asthma Management Program (CAMP) subset of Asthma BRIDGE childhood asthmatics (n = 527) ages 5–12 with prenatal smoking exposure data available. Using beta-regression, we screened loci for differential methylation related to prenatal smoke exposure, adjusting for gender, age and clinical site, and accounting for multiple comparisons by FDR. Results: Of 27,578 loci evaluated, 22,131 (80%) passed quality control assessment and were analyzed. Sixty-five children (12%) had a history of prenatal smoke exposure. At an FDR of 0.05, we identified 19 CpG loci significantly associated with prenatal smoke, of which two replicated in two independent populations. Exposure was associated with a 2% increase in mean CpG methylation in FRMD4A (p = 0.01) and Cllorf52 (p = 0.001) compared to no exposure. Four additional genes, XPNPEP1, PPEF2, SMPD3 and CRYGN, were nominally associated in at least one replication group. Conclusions: These data suggest that prenatal exposure to tobacco smoke is associated with reproducible epigenetic changes that persist well into childhood. However, the biological significance of these altered loci remains unknown.

There is intense interest in the role of vitamin D in the development of asthma and allergies. However, studies differ on whether a higher vitamin D intake or status in pregnancy or at birth is protective against asthma and allergies. To address this uncertainty, the Vitamin D Antenatal Asthma Reduction Trial (VDAART) was developed. VDAART is a randomized, double-blind, placebo-controlled trial of vitamin D supplementation in pregnant women to determine whether prenatal supplementation can prevent the development of asthma and allergies in the women’s offspring. A secondary aim is to determine whether vitamin D supplementation can prevent the development of pregnancy complications, such as preeclampsia, preterm birth, and gestational diabetes. Women were randomized to the treatment arm of 4,000 IU/day of vitamin D3 plus a daily multivitamin that contained 400 IU of vitamin D3 or the placebo arm of placebo plus a multivitamin that contained 400 IU daily of vitamin D3. Women who were between the gestational ages of 10–18 weeks were randomized from three clinical centers across the United States – Boston Medical Center, Washington University in St. Louis, and Kaiser Permanente Southern California Region (San Diego, CA). Supplementation took place throughout pregnancy. Monthly monitoring of urinary calcium to creatinine ratio was performed in addition to medical record review for adverse events. Offspring are being evaluated quarterly through questionnaires and yearly during in-person visits until the 3rd birthday of the child. Ancillary studies will investigate neonatal T-regulatory cell function, maternal vaginal flora, and maternal and child intestinal flora.

We investigate 71 single nucleotide polymorphisms (SNPs) identified in meta-analytic studies of genome-wide association studies (GWAS) of breast cancer, the majority of which are located in intergenic or intronic regions. To explore regulatory impacts of these variants we conducted expression quantitative loci (eQTL) analyses on tissue samples from 376 invasive postmenopausal breast cancer cases in the Nurses’ Health Study (NHS) diagnosed from 1990–2004. Expression analysis was conducted on all formalin-fixed paraffin-embedded (FFPE) tissue samples (and on 264 adjacent normal samples) using the Affymetrix Human Transcriptome Array. Significance and ranking of associations between tumor receptor status and expression variation was preserved between NHS FFPE and TCGA fresh-frozen sample sets (Spearman r = 0.85, p<10^-10 for 17 of the 21 Oncotype DX recurrence signature genes). At an FDR threshold of 10%, we identified 27 trans-eQTLs associated with expression variation in 217 distinct genes. SNP-gene associations can be explored using an open-source interactive browser distributed in a Bioconductor package. Using a new a procedure for testing hypotheses relating SNP content to expression patterns in gene sets, defined as molecular function pathways, we find that loci on 6q14 and 6q25 affect various gene sets and molecular pathways (FDR < 10%). Although the ultimate biological interpretation of the GWAS-identified variants remains to be uncovered, this study validates the utility of expression analysis of this FFPE expression set for more detailed integrative analyses.

Human-associated microbial communities have a crucial role in determining our health and well-being and this has led to the continuing development of microbiome-based therapies such as faecal microbiota transplantation. These microbial communities are very complex, dynamic and highly personalized ecosystems exhibiting a high degree of inter-individual variability in both species assemblages and abundance profiles9. It is not known whether the underlying ecological dynamics of these communities, which can be parameterized by growth rates, and intra- and inter-species interactions in population dynamics models, are largely host-independent (that is, universal) or host-specific. If the inter-individual variability reflects host-specific dynamics due to differences in host lifestyle, physiology or genetics, then generic microbiome manipulations may have unintended consequences, rendering them ineffective or even detrimental. Alternatively, microbial ecosystems of different subjects may exhibit universal dynamics, with the inter-individual variability mainly originating from differences in the sets of colonizing species Here we develop a new computational method to characterize human microbial dynamics. By applying this method to cross-sectional data from two large-scale metagenomic studies—the Human Microbiome Project and the Student Microbiome Project—we show that gut and mouth microbiomes display pronounced universal dynamics, whereas communities associated with certain skin sites are probably shaped by differences in the host environment. Notably, the universality of gut microbial dynamics is not observed in subjects with recurrent Clostridium difficile infection but is observed in the same set of subjects after faecal microbiota transplantation. These results fundamentally improve our understanding of the processes that shape human microbial ecosystems, and pave the way to designing general microbiome-based therapies.

Background: Little is known about the role of most asthma susceptibility genes during human lung development. Genetic determinants for normal lung development are not only important early in life, but also for later lung function.

Objective: To investigate the role of expression patterns of well-defined asthma susceptibility genes during human and murine lung development. We hypothesized that genes influencing normal airways development would be over-represented by genes associated with asthma.

Methods: Asthma genes were first identified via comprehensive search of the current literature. Next, we analyzed their expression patterns in the developing human lung during the pseudoglandular (gestational age, 7-16 weeks) and canalicular (17-26 weeks) stages of development, and in the complete developing lung time series of 3 mouse strains: A/J, SW, C57BL6.

Results: In total, 96 genes with association to asthma in at least two human populations were identified in the literature. Overall, there was no significant over-representation of the asthma genes among genes differentially expressed during lung development, although trends were seen in the human (Odds ratio, OR 1.22, confidence interval, CI 0.90-1.62) and C57BL6 mouse (OR 1.41, CI 0.92-2.11) data. However, differential expression of some asthma genes was consistent in both developing human and murine lung, e.g. NOD1, EDN1, CCL5, RORA and HLA-G. Among the asthma genes identified in genome wide association studies, ROBO1, RORA, HLA-DQB1, IL2RB and PDE10A were differentially expressed during human lung development.

Conclusions: Our data provide insight about the role of asthma susceptibility genes during lung development and suggest common mechanisms underlying lung morphogenesis and pathogenesis of respiratory diseases.

Background: Genome-wide association studies have yet to identify the majority of genetic variants involved in asthma. We hypothesized that expression quantitative trait locus (eQTL) mapping can identify novel asthma genes by enabling prioritization of putative functional variants for association testing. Objective: We evaluated 6,706 cis-acting expression-associated variants (eSNP) identified through a genome-wide eQTL survey of CD4+ lymphocytes for association with asthma. Methods: eSNP were tested for association with asthma in 359 asthma cases and 846 controls from the Childhood Asthma Management Program, with verification using family-based testing. Significant associations were tested for replication in 579 parent-child trios with asthma from Costa Rica. Further functional validation was performed by Formaldehyde Assisted Isolation of Regulatory Elements (FAIRE)-qPCR and Chromatin-Immunoprecipitation (ChIP)-PCR in lung derived epithelial cell lines (Beas-2B and A549) and Jurkat cells, a leukemia cell line derived from T lymphocytes. Results: Cis-acting eSNP demonstrated associations with asthma in both cohorts. We confirmed the previously-reported association of ORMDL3/GSDMB variants with asthma (combined p=2.9 × 108). Reproducible associations were also observed for eSNP in three additional genes: FADS2 (p=0.002), NAGA (p=0.0002), and F13A1 (p=0.0001). We subsequently demonstrated that FADS2 mRNA is increased in CD4+ lymphocytes in asthmatics, and that the associated eSNPs reside within DNA segments with histone modifications that denote open chromatin status and confer enhancer activity. Conclusions: Our results demonstrate the utility of eQTL mapping in the identification of novel asthma genes, and provide evidence for the importance of FADS2, NAGA, and F13A1 in the pathogenesis of asthma.

We describe Bioconductor infrastructure for representing and computing on annotated genomic ranges and integrating genomic data with the statistical computing features of R and its extensions. At the core of the infrastructure are three packages: IRanges, GenomicRanges, and GenomicFeatures. These packages provide scalable data structures for representing annotated ranges on the genome, with special support for transcript structures, read alignments and coverage vectors. Computational facilities include efficient algorithms for overlap and nearest neighbor detection, coverage calculation and other range operations. This infrastructure directly supports more than 80 other Bioconductor packages, including those for sequence analysis, differential expression analysis and visualization.

Introduction: Management of persistently non-adherent youth living with HIV (YLHIV) with virologic failure (VF) on combination antiretroviral therapy (cART) remains challenging. One strategy has been using 3TC/ FTC monotherapy (3TC/FTC), which in the presence of the M184V resistance mutation, does not suppress viral replication nor select for additional drug resistance mutations, and reduces viral fitness with limited side effects. P1094 compared the immunologic outcome of continuing failing cART vs. switching to 3TC/FTC as a “bridging strategy” to subsequent suppressive cART for non-adherent YLHIV with pre-existing M184V resistance. Materials & methods Participants with documented nonadherence, M184V mutation, CD4+ T cell count ≥100 cells/mm3 and VF (HIV-1 plasma RNA ≥400 copies/mL (2.6 log10 HIV-1 RNA) were enrolled and randomized to continue failing cART vs. switch to 3TC/FTC. The primary endpoint (time to ≥30% CD4+ T cell decline or development of CDC class C events) at 28-weeks were assessed by Kaplan-Meier (K-M) curves in an intent-to-treat analysis. Results: Thirty-three perinatally acquired YLHIV participants (16 continuing cART and 17 3TC/FTC) enrolled in the study. The median age, entry CD4+ T cell count, and viral load were 15 years (Inter-quartile range (IQR) 14–20), 472 cells/mm3 (IQR 384–651), and 4.0 log10HIV-1 RNA copies/ml (IQR 3.2–4.5), respectively. Five participants, all in the 3TC/FTC arm, reached the primary endpoint for absolute CD4+ T cell decline (p = 0.02, exact log-rank test comparing monotherapy to cART). The Kaplan-Meier estimate of probability of primary endpoint on 3TC/FTC at 28 weeks was 0.41 (standard error 0.14). There were no CDC class C events or deaths and no statistically significant difference in frequencies of adverse events between the arms. Conclusions: Non-adherent participants randomized to 3TC/FTC were more likely than those maintained on failing cART to experience a confirmed decline in CD4+ count of ≥30%. Although this study suffers from limitations of small sample size and premature discontinuation, the randomized comparison to continuing failing cART indicates that 3TC/FTC provides inferior protection from immunologic deterioration for non-adherent youth with M184V resistance. Better alternatives to 3TC/FTC such as ART with higher barriers to resistance and novel adherence and treatment strategies for nonadherent youth are urgently needed. Trial registration Clinical Trials.gov NCT01338025