Martin Morgan

Orchestrating high-throughput genomic analysis with Bioconductor

Bioconductor is an open-source, open-development software project for the analysis and comprehension of high-throughput data in genomics and molecular biology. The project aims to enable interdisciplinary research, collaboration and rapid development of scientific software. Based on the statistical programming language R, Bioconductor comprises 934 interoperable packages contributed by a large, diverse community of scientists. Packages cover a range of bioinformatic and statistical applications. They undergo formal initial review and continuous automated testing. We present an overview for prospective users and contributors.

Software for Computing and Annotating Genomic Ranges

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.

Genome-wide MyoD Binding in Skeletal Muscle Cells: A Potential for Broad Cellular Reprogramming

Recent studies have demonstrated that MyoD initiates a feed-forward regulation of skeletal muscle gene expression, predicting that MyoD binds directly to many genes expressed during differentiation. We have used chromatin immunoprecipitation and high-throughput sequencing to identify genome-wide binding of MyoD in several skeletal muscle cell types. As anticipated, MyoD preferentially binds to a VCASCTG sequence that resembles the in vitro-selected site for a MyoD:E-protein heterodimer, and MyoD binding increases during differentiation at many of the regulatory regions of genes expressed in skeletal muscle. Unanticipated findings were that MyoD was constitutively bound to thousands of additional sites in both myoblasts and myotubes, and that the genome-wide binding of MyoD was associated with regional histone acetylation. Therefore, in addition to regulating muscle gene expression, MyoD binds genome wide and has the ability to broadly alter the epigenome in myoblasts and myotubes.

Bacterial Communities in Women with Bacterial Vaginosis: High Resolution Phylogenetic Analyses Reveal Relationships of Microbiota to Clinical Criteria

Background Bacterial vaginosis (BV) is a common condition that is associated with numerous adverse health outcomes and is characterized by poorly understood changes in the vaginal microbiota. We sought to describe the composition and diversity of the vaginal bacterial biota in women with BV using deep sequencing of the 16S rRNA gene coupled with species-level taxonomic identification. We investigated the associations between the presence of individual bacterial species and clinical diagnostic characteristics of BV. Methodology/Principal Findings Broad-range 16S rRNA gene PCR and pyrosequencing were performed on vaginal swabs from 220 women with and without BV. BV was assessed by Amsel’s clinical criteria and confirmed by Gram stain. Taxonomic classification was performed using phylogenetic placement tools that assigned 99% of query sequence reads to the species level. Women with BV had heterogeneous vaginal bacterial communities that were usually not dominated by a single taxon. In the absence of BV, vaginal bacterial communities were dominated by either Lactobacillus crispatus or Lactobacillus iners. Leptotrichia amnionii and Eggerthella sp. were the only two BV-associated bacteria (BVABs) significantly associated with each of the four Amsel’s criteria. Co-occurrence analysis revealed the presence of several sub-groups of BVABs suggesting metabolic co-dependencies. Greater abundance of several BVABs was observed in Black women without BV. Conclusions/Significance The human vaginal bacterial biota is heterogeneous and marked by greater species richness and diversity in women with BV; no species is universally present. Different bacterial species have different associations with the four clinical criteria, which may account for discrepancies often observed between Amsel and Nugent (Gram stain) diagnostic criteria. Several BVABs exhibited race-dependent prevalence when analyzed in separate groups by BV status which may contribute to increased incidence of BV in Black women. Tools developed in this project can be used to study microbial ecology in diverse settings at high resolution.

ShortRead: a bioconductor package for input, quality assessment and exploration of high-throughput sequence data

Summary: ShortRead is a package for input, quality assessment, manipulation and output of high-throughput sequencing data. ShortRead is provided in the R and Bioconductor environments, allowing ready access to additional facilities for advanced statistical analysis, data transformation, visualization and integration with diverse genomic resources. Availability and Implementation: This package is implemented in R and available at the Bioconductor web site; the package contains a ‘vignette’ outlining typical work flows. Contact: mtmorgan@fhcrc.org

Recent genetic selection in the ancestral admixture of Puerto Ricans.

Recent studies have used dense markers to examine the human genome in ancestrally homogeneous populations for hallmarks of selection. No genomewide studies have focused on recently admixed groups--populations that have experienced admixing among continentally divided ancestral populations within the past 200-500 years. New World admixed populations are unique in that they represent the sudden confluence of geographically diverged genomes with novel environmental challenges. Here, we present a novel approach for studying selection by examining the genomewide distribution of ancestry in the genetically admixed Puerto Ricans. We find strong statistical evidence of recent selection in three chromosomal regions, including the human leukocyte antigen region on chromosome 6p, chromosome 8q, and chromosome 11q. Two of these regions harbor genes for olfactory receptors. Interestingly, all three regions exhibit deficiencies in the European-ancestry proportion.

Metabolic Signatures of Bacterial Vaginosis

ABSTRACT Bacterial vaginosis (BV) is characterized by shifts in the vaginal microbiota from Lactobacillus dominant to a microbiota with diverse anaerobic bacteria. Few studies have linked specific metabolites with bacteria found in the human vagina. Here, we report dramatic differences in metabolite compositions and concentrations associated with BV using a global metabolomics approach. We further validated important metabolites using samples from a second cohort of women and a different platform to measure metabolites. In the primary study, we compared metabolite profiles in cervicovaginal lavage fluid from 40 women with BV and 20 women without BV. Vaginal bacterial representation was determined using broad-range PCR with pyrosequencing and concentrations of bacteria by quantitative PCR. We detected 279 named biochemicals; levels of 62% of metabolites were significantly different in women with BV. Unsupervised clustering of metabolites separated women with and without BV. Women with BV have metabolite profiles marked by lower concentrations of amino acids and dipeptides, concomitant with higher levels of amino acid catabolites and polyamines. Higher levels of the signaling eicosanoid 12-hydroxyeicosatetraenoic acid (12-HETE), a biomarker for inflammation, were noted in BV. Lactobacillus crispatus and Lactobacillus jensenii exhibited similar metabolite correlation patterns, which were distinct from correlation patterns exhibited by BV-associated bacteria. Several metabolites were significantly associated with clinical signs and symptoms (Amsel criteria) used to diagnose BV, and no metabolite was associated with all four clinical criteria. BV has strong metabolic signatures across multiple metabolic pathways, and these signatures are associated with the presence and concentrations of particular bacteria. IMPORTANCE Bacterial vaginosis (BV) is a common but highly enigmatic condition that is associated with adverse outcomes for women and their neonates. Small molecule metabolites in the vagina may influence host physiology, affect microbial community composition, and impact risk of adverse health outcomes, but few studies have comprehensively studied the metabolomics profile of BV. Here, we used mass spectrometry to link specific metabolites with particular bacteria detected in the human vagina by PCR. BV was associated with strong metabolic signatures across multiple pathways affecting amino acid, carbohydrate, and lipid metabolism, highlighting the profound metabolic changes in BV. These signatures were associated with the presence and concentrations of particular vaginal bacteria, including some bacteria yet to be cultivated, thereby providing clues as to the microbial origin of many metabolites. Insights from this study provide opportunities for developing new diagnostic markers of BV and novel approaches for treatment or prevention of BV. Bacterial vaginosis (BV) is a common but highly enigmatic condition that is associated with adverse outcomes for women and their neonates. Small molecule metabolites in the vagina may influence host physiology, affect microbial community composition, and impact risk of adverse health outcomes, but few studies have comprehensively studied the metabolomics profile of BV. Here, we used mass spectrometry to link specific metabolites with particular bacteria detected in the human vagina by PCR. BV was associated with strong metabolic signatures across multiple pathways affecting amino acid, carbohydrate, and lipid metabolism, highlighting the profound metabolic changes in BV. These signatures were associated with the presence and concentrations of particular vaginal bacteria, including some bacteria yet to be cultivated, thereby providing clues as to the microbial origin of many metabolites. Insights from this study provide opportunities for developing new diagnostic markers of BV and novel approaches for treatment or prevention of BV.

More Than Meets the Eye: Associations of Vaginal Bacteria with Gram Stain Morphotypes Using Molecular Phylogenetic Analysis

Bacterial vaginosis (BV) is a highly prevalent condition associated with adverse health outcomes. Gram stain analysis of vaginal fluid is the standard for confirming the diagnosis of BV, wherein abundances of key bacterial morphotypes are assessed. These Lactobacillus, Gardnerella, Bacteroides, and Mobiluncus morphotypes were originally linked to particular bacterial species through cultivation studies, but no studies have systematically investigated associations between uncultivated bacteria detected by molecular methods and Gram stain findings. In this study, 16S-rRNA PCR/pyrosequencing was used to examine associations between vaginal bacteria and bacterial morphotypes in 220 women with and without BV. Species-specific quantitative PCR (qPCR) and fluorescence in Situ hybridization (FISH) methods were used to document concentrations of two bacteria with curved rod morphologies: Mobiluncus and the fastidious BV-associated bacterium-1 (BVAB1). Rank abundance of vaginal bacteria in samples with evidence of curved gram-negative rods showed that BVAB1 was dominant (26.1%), while Mobiluncus was rare (0.2% of sequence reads). BVAB1 sequence reads were associated with Mobiluncus morphotypes (p<0.001). Among women with curved rods, mean concentration of BVAB1 DNA was 2 log units greater than Mobiluncus (p<0.001) using species-specific quantitative PCR. FISH analyses revealed that mean number of BVAB1 cells was 2 log units greater than Mobiluncus cells in women with highest Nugent score (p<0.001). Prevotella and Porphyromonas spp. were significantly associated with the “Bacteroides morphotype,” whereas Bacteroides species were rare. Gram-negative rods designated Mobiluncus morphotypes on Gram stain are more likely BVAB1. These findings provide a clearer picture of the bacteria associated with morphotypes on vaginal Gram stain.

MYC-driven tumorigenesis is inhibited by WRN syndrome gene deficiency

MYC-induced DNA damage is exacerbated in WRN-deficient cells, leading to replication stress and accelerated cellular senescence. To determine whether WRN deficiency impairs MYC-driven tumor development, we used both xenograft and autochthonous tumor models. Conditional silencing of WRN expression in c-MYC overexpressing non–small cell lung cancer xenografts impaired both tumor establishment and tumor growth. This inhibitory effect of WRN knockdown was accompanied by increased DNA damage, decreased proliferation, and tumor necrosis. In the Eμ-Myc mouse model of B-cell lymphoma, a germline mutation in the helicase domain of Wrn (WrnΔhel/Δhel) resulted in a significant delay in emergence of lethal lymphomas, extending tumor-free survival by more than 30%. Analysis of preneoplastic B cells from Eμ-Myc Wrn mutant mice revealed increased DNA damage, elevation of senescence markers, and decreased proliferation in comparison with cells from age-matched Eμ-Myc mice. Immunohistochemical and global gene expression analysis of overt Eμ-Myc WrnΔhel/Δhel lymphomas showed a marked increase in expression of the CDK inhibitor, p16Ink4a, as well as elevation of TAp63, a known mediator of senescence. Collectively, these studies show that in the context of Myc-associated tumorigenesis, loss of Wrn amplifies the DNA damage response, both in preneoplastic and neoplastic tissue, engaging activation of tumor suppressor pathways. This leads to inhibition of tumor growth and prolonged tumor-free survival. Targeting WRN or its enzymatic function could prove to be an effective strategy in the treatment of MYC-associated cancers. Mol Cancer Res; 10(4); 535–45. ©2012 AACR.

Gene expression changes in adipose tissue with diet- and/or exercise-induced weight loss

Adipose tissue plays a role in obesity-related cancers via increased production of inflammatory factors, steroid hormones, and altered adipokines. The impact of weight loss on adipose tissue gene expression may provide insights into pathways linking obesity with cancer risk. We conducted an ancillary study within a randomized trial of diet, exercise, or combined diet + exercise versus control among overweight/obese postmenopausal women. In 45 women, subcutaneous adipose tissue biopsies were conducted at baseline and after 6 months, and changes in adipose tissue gene expression were determined by microarray with an emphasis on prespecified candidate pathways as well as by unsupervised clustering of more than 37,000 transcripts (Illumina). Analyses were conducted first by randomization group and then by degree of weight change at 6-months in all women combined. At 6 months, diet, exercise, and diet + exercise participants lost a mean of 8.8, 2.5, and 7.9 kg (all P < 0.05 vs. no change in controls). There was no significant change in candidate gene expression by intervention group. In analysis by weight change category, greater weight loss was associated a decrease in 17β-hydroxysteroid dehydrogenase-1 (HSD17B1, Ptrend < 0.01) and leptin (LEP, Ptrend < 0.01) expression, and marginally significant increased expression of estrogen receptor-1 (ESR1, Ptrend = 0.08) and insulin-like growth factor–binding protein-3 (IGFBP3, Ptrend = 0.08). Unsupervised clustering revealed 83 transcripts with statistically significant changes. Multiple gene expression changes correlated with changes in associated serum biomarkers. Weight loss was associated with changes in adipose tissue gene expression after 6 months, particularly in two pathways postulated to link obesity and cancer, that is, steroid hormone metabolism and IGF signaling. Cancer Prev Res; 6(3); 217–31. ©2013 AACR.