Bioconductor: open software development for computational biology and bioinformatics

The Bioconductor project is an initiative for the collaborative creation of extensible software for computational biology and bioinformatics. The goals of the project include: fostering collaborative development and widespread use of innovative software, reducing barriers to entry into interdisciplinary scientific research, and promoting the achievement of remote reproducibility of research results. We describe details of our aims and methods, identify current challenges, compare Bioconductor to other open bioinformatics projects, and provide working examples.

Enterotypes of the human gut microbiome

Our knowledge of species and functional composition of the human gut microbiome is rapidly increasing, but it is still based on very few cohorts and little is known about variation across the world. By combining 22 newly sequenced faecal metagenomes of individuals from four countries with previously published data sets, here we identify three robust clusters (referred to as enterotypes hereafter) that are not nation or continent specific. We also confirmed the enterotypes in two published, larger cohorts, indicating that intestinal microbiota variation is generally stratified, not continuous. This indicates further the existence of a limited number of well-balanced host-microbial symbiotic states that might respond differently to diet and drug intake. The enterotypes are mostly driven by species composition, but abundant molecular functions are not necessarily provided by abundant species, highlighting the importance of a functional analysis to understand microbial communities. Although individual host properties such as body mass index, age, or gender cannot explain the observed enterotypes, data-driven marker genes or functional modules can be identified for each of these host properties. For example, twelve genes significantly correlate with age and three functional modules with the body mass index, hinting at a diagnostic potential of microbial markers.

affy--analysis of Affymetrix GeneChip data at the probe level

MOTIVATION

The processing of the Affymetrix GeneChip data has been a recent focus for data analysts. Alternatives to the original procedure have been proposed and some of these new methods are widely used.

RESULTS

The affy package is an R package of functions and classes for the analysis of oligonucleotide arrays manufactured by Affymetrix. The package is currently in its second release, affy provides the user with extreme flexibility when carrying out an analysis and make it possible to access and manipulate probe intensity data. In this paper, we present the main classes and functions in the package and demonstrate how they can be used to process probe-level data. We also demonstrate the importance of probe-level analysis when using the Affymetrix GeneChip platform.

Comparative genomics of Listeria species

Listeria monocytogenes is a food-borne pathogen with a high mortality rate that has also emerged as a paradigm for intracellular parasitism. We present and compare the genome sequences of L. monocytogenes (2,944,528 base pairs) and a nonpathogenic species, L. innocua (3,011,209 base pairs). We found a large number of predicted genes encoding surface and secreted proteins, transporters, and transcriptional regulators, consistent with the ability of both species to adapt to diverse environments. The presence of 270 L. monocytogenes and 149 L. innocua strain-specific genes (clustered in 100 and 63 islets, respectively) suggests that virulence in Listeria results from multiple gene acquisition and deletion events.

A new non-linear normalization method for reducing variability in DNA microarray experiments

BACKGROUND

Microarray data are subject to multiple sources of variation, of which biological sources are of interest whereas most others are only confounding. Recent work has identified systematic sources of variation that are intensity-dependent and non-linear in nature. Systematic sources of variation are not limited to the differing properties of the cyanine dyes Cy(5) and Cy(3) as observed in cDNA arrays, but are the general case for both oligonucleotide microarray (Affymetrix GeneChips) and cDNA microarray data. Current normalization techniques are most often linear and therefore not capable of fully correcting for these effects.

RESULTS

We present here a simple and robust non-linear method for normalization using array signal distribution analysis and cubic splines. These methods compared favorably to normalization using robust local-linear regression (lowess). The application of these methods to oligonucleotide arrays reduced the relative error between replicates by 5-10% compared with a standard global normalization method. Application to cDNA arrays showed improvements over the standard method and over Cy(3)-Cy(5) normalization based on dye-swap replication. In addition, a set of known differentially regulated genes was ranked higher by the t-test. In either cDNA or Affymetrix technology, signal-dependent bias was more than ten times greater than the observed print-tip or spatial effects.

CONCLUSIONS

Intensity-dependent normalization is important for both high-density oligonucleotide array and cDNA array data. Both the regression and spline-based methods described here performed better than existing linear methods when assessed on the variability of replicate arrays. Dye-swap normalization was less effective at Cy(3)-Cy(5) normalization than either regression or spline-based methods alone.

Bacterial nitric oxide extends the lifespan of C. elegans

Nitric oxide (NO) is an important signaling molecule in multicellular organisms. Most animals produce NO from L-arginine via a family of dedicated enzymes known as NO synthases (NOSes). A rare exception is the roundworm Caenorhabditis elegans, which lacks its own NOS. However, in its natural environment, C. elegans feeds on Bacilli that possess functional NOS. Here, we demonstrate that bacterially derived NO enhances C. elegans longevity and stress resistance via a defined group of genes that function under the dual control of HSF-1 and DAF-16 transcription factors. Our work provides an example of interspecies signaling by a small molecule and illustrates the lifelong value of commensal bacteria to their host.

Alternative mapping of probes to genes for Affymetrix chips

BACKGROUND

Short oligonucleotide arrays have several probes measuring the expression level of each target transcript. Therefore the selection of probes is a key component for the quality of measurements. However, once probes have been selected and synthesized on an array, it is still possible to re-evaluate the results using an updated mapping of probes to genes, taking into account the latest biological knowledge available.

METHODS

We investigated how probes found on recent commercial microarrays for human genes (Affymetrix HG-U133A) were matching a recent curated collection of human transcripts: the NCBI RefSeq database. We also built mappings and used them in place of the original probe to genes associations provided by the manufacturer of the arrays.

RESULTS

In a large number of cases, 36%, the probes matching a reference sequence were consistent with the grouping of probes by the manufacturer of the chips. For the remaining cases there were discrepancies and we show how that can affect the analysis of data.

CONCLUSIONS

While the probes on Affymetrix arrays remain the same for several years, the biological knowledge concerning the genomic sequences evolves rapidly. Using up-to-date knowledge can apparently change the outcome of an analysis.

A DNA metabarcoding study of a primate dietary diversity and plasticity across its entire fragmented range

In tropical regions, most primary ecosystems have been replaced by mosaic landscapes in which species must cope with a large shift in the distribution of their habitat and associated food resources. Primates are particularly vulnerable to habitat modifications. Most species persist in small fragments surrounded by complex human-mediated matrices whose structure and connectivity may strongly influence their dispersal and feeding behavior. Behavioral plasticity appears to be a crucial parameter governing the ability of organisms to exploit the resources offered by new matrix habitats and thus to persist in fragmented habitats. In this study, we were interested in the dietary plasticity of the golden-crowned sifaka (Propithecus tattersalli), an endangered species of lemur, found only in the Daraina region in north-eastern Madagascar. We used a DNA-based approach combining the barcoding concept and Illumina next-generation sequencing to (i) describe the species diet across its entire range and (ii) evaluate the influence of landscape heterogeneity on diet diversity and composition. Faeces from 96 individuals were sampled across the entire species range and their contents were analyzed using the trnL metabarcoding approach. In parallel, we built a large DNA reference database based on a checklist of the plant species of the Daraina region. Our results suggest that golden-crowned sifakas exhibit remarkable dietary diversity with at least 130 plant species belonging to 80 genera and 49 different families. We highlighted an influence of both habitat type and openness on diet composition suggesting a high flexibility of foraging strategies. Moreover, we observed the presence of numerous cultivated and naturalized plants in the faeces of groups living in forest edge areas. Overall, our findings support our initial expectation that P. tattersalli is able to cope with the current level of alteration of the landscape and confirm our previous results on the distribution and the dispersal ability of this species.

Photoperiodic induction of synchronous flowering near the Equator

In tropical rainforests, 30-65% of tree species grow at densities of less than one individual per hectare. At these low population densities, successful cross-pollination relies on synchronous flowering. In rainforests with low climatic seasonality, photoperiodic control is the only reliable mechanism for inducing synchronous flowering. This poses a problem because there is no variation in day length at the Equator. Here we propose a new mechanism of photoperiodic timekeeping based on the perception of variation in sunrise or sunset time, which explains and predicts the annually repeated, staggered, synchronous and bimodal flowering of many tree species in Amazonian rainforests near the Equator.

Transcriptional interactions suggest niche segregation among microorganisms in the human gut

The human gastrointestinal (GI) tract is the habitat for hundreds of microbial species, of which many cannot be cultivated readily, presumably because of the dependencies between species¹. Studies of microbial co-occurrence in the gut have indicated community substructures that may reflect functional and metabolic interactions between cohabiting species^2,3. To move beyond species co-occurrence networks, we systematically identified transcriptional interactions between pairs of coexisting gut microbes using metagenomics and microarray-based metatranscriptomics data from 233 stool samples from Europeans. In 102 significantly interacting species pairs, the transcriptional changes led to a reduced expression of orthologous functions between the coexisting species. Specific species-species transcriptional interactions were enriched for functions important for H₂ and CO₂ homeostasis, butyrate biosynthesis, ATP-binding cassette (ABC) transporters, flagella assembly and bacterial chemotaxis, as well as for the metabolism of carbohydrates, amino acids and cofactors. The analysis gives the first insight into the microbial community-wide transcriptional interactions, and suggests that the regulation of gene expression plays an important role in species adaptation to coexistence and that niche segregation takes place at the transcriptional level.

Mechanism of biofilm-mediated stress resistance and lifespan extension in C. elegans

Bacteria naturally form communities of cells known as biofilms. However the physiological roles of biofilms produced by non-pathogenic microbiota remain largely unknown. To assess the impact of a biofilm on host physiology we explored the effect of several non-pathogenic biofilm-forming bacteria on Caenorhabditis elegans. We show that biofilm formation by Bacillus subtilis, Lactobacillus rhamnosus and Pseudomonas fluorescens induces C. elegans stress resistance. Biofilm also protects against pathogenic infection and prolongs lifespan. Total mRNA analysis identified a set of host genes that are upregulated in response to biofilm formation by B. subtilis. We further demonstrate that mtl-1 is responsible for the biofilm-mediated increase in oxidative stress resistance and lifespan extension. Induction of mtl-1 and hsp-70 promotes biofilm-mediated thermotolerance. ilys-2 activity accounts for biofilm-mediated resistance to Pseudomonas aeruginosa killing. These results reveal the importance of non-pathogenic biofilms for host physiology and provide a framework to study commensal biofilms in higher organisms.

TLR2 controls intestinal carcinogen detoxication by CYP1A1

Intestinal cytochrome P450 subclass 1A1 (CYP1A1) contributes to a metabolic “shield” protecting the host from ingested carcinogens such as polycyclic aromatic hydrocarbons (PAH). The expression of CYP1 (including CYP1A2 and CYP1B1) is considered to depend solely on a heterodimeric transcription factor consisting of the arylhydrocarbon receptor (AHR) and the AHR nuclear translocator (ARNT). So far, no interference has been noted between the regulation of CYP1 and the activation of Toll-like receptor 2 (TLR2), which modulates the inflammatory response to bacterial cell wall components in immune cells and enterocytes. Here we report that intestinal CYP1A1 is silenced in TLR2-deficient mice, even when under exposure to the carcinogenic PAH benzo[a]pyrene (BaP). In contrast, hepatic CYP1A1 was moderately induced in TLR2-deficient mice without restoring their ability to clear BaP from systemic circulation, as present in wild-type animals. After feeding of BaP for 21 days, only TLR2^−/− mice, but not their wild type littermates developed polyps in the colon. Gene expressions and protein concentrations of AHR and ARNT in the intestine did not differ between the genotypes. In conclusion, the presence of ligands for TLR2 of bacterial origin seems to be crucial for detoxication of luminal carcinogens by CYP1A1 in the intestine. This unprecedented finding indicates a complex interplay between the immune system of the host and intestinal bacteria with detoxication mechanisms. This highlights the relevance of intestinal microbiota when trying to unravel pathways present in mammals and opens new perspectives for research in human health.

Patterns of diversification amongst tropical regions compared: a case study in Sapotaceae

Species diversity is unequally distributed across the globe, with the greatest concentration occurring in the tropics. Even within the tropics, there are significant differences in the numbers of taxa found in each continental region. Manilkara is a pantropical genus of trees in the Sapotaceae comprising c. 78 species. Its distribution allows for biogeographic investigation and testing of whether rates of diversification differ amongst tropical regions. The age and geographical origin of Manilkara are inferred to determine whether Gondwanan break-up, boreotropical migration or long distance dispersal have shaped its current disjunct distribution. Diversification rates through time are also analyzed to determine whether the timing and tempo of speciation on each continent coincides with geoclimatic events. Bayesian analyses of nuclear (ITS) and plastid (rpl32-trnL, rps16-trnK, and trnS-trnFM) sequences were used to reconstruct a species level phylogeny of Manilkara and related genera in the tribe Mimusopeae. Analyses of the nuclear data using a fossil-calibrated relaxed molecular clock indicate that Manilkara evolved 32-29 million years ago (Mya) in Africa. Lineages within the genus dispersed to the Neotropics 26-18 Mya and to Asia 28-15 Mya. Higher speciation rates are found in the Neotropical Manilkara clade than in either African or Asian clades. Dating of regional diversification correlates with known palaeoclimatic events. In South America, the divergence between Atlantic coastal forest and Amazonian clades coincides with the formation of drier Cerrado and Caatinga habitats between them. In Africa diversification coincides with Tertiary cycles of aridification and uplift of the east African plateaux. In Southeast Asia dispersal may have been limited by the relatively recent emergence of land in New Guinea and islands further east c. 10 Mya.

Genomic profiling of thousands of candidate polymorphisms predicts risk of relapse in 778 Danish and German childhood acute lymphoblastic leukemia patients

Childhood acute lymphoblastic leukemia survival approaches 90%. New strategies are needed to identify the 10-15% who evade cure. We applied targeted, sequencing-based genotyping of 25 000 to 34 000 preselected potentially clinically relevant single-nucleotide polymorphisms (SNPs) to identify host genome profiles associated with relapse risk in 352 patients from the Nordic ALL92/2000 protocols and 426 patients from the German Berlin-Frankfurt-Munster (BFM) ALL2000 protocol. Patients were enrolled between 1992 and 2008 (median follow-up: 7.6 years). Eleven cross-validated SNPs were significantly associated with risk of relapse across protocols. SNP and biologic pathway level analyses associated relapse risk with leukemia aggressiveness, glucocorticosteroid pharmacology/response and drug transport/metabolism pathways. Classification and regression tree analysis identified three distinct risk groups defined by end of induction residual leukemia, white blood cell count and variants in myeloperoxidase (MPO), estrogen receptor 1 (ESR1), lamin B1 (LMNB1) and matrix metalloproteinase-7 (MMP7) genes, ATP-binding cassette transporters and glucocorticosteroid transcription regulation pathways. Relapse rates ranged from 4% (95% confidence interval (CI): 1.6-6.3%) for the best group (72% of patients) to 76% (95% CI: 41-90%) for the worst group (5% of patients, P<0.001). Validation of these findings and similar approaches to identify SNPs associated with toxicities may allow future individualized relapse and toxicity risk-based treatments adaptation.

Mobilization of Philadelphia-negative peripheral blood progenitor cells with chemotherapy and rhuG-CSF in chronic myelogenous leukaemia patients with a poor response to interferon-alpha

The purpose of this cooperative study was to evaluate the quantity and quality of Ph1-negative progenitor cells mobilized in the peripheral blood of patients with chronic myelogenous leukaemia soon after aplasia induced by chemotherapy. 32 patients ineligible for allografting who were cytogenetically refractory to interferon-alpha (IFN-alpha) were entered into this study. The chronic phase varied widely, with a median duration of 17 months (range 3-90 months). All patients were treated with intensive conventional chemotherapy regimens and recombinant human granulocyte colony-stimulating factor (rhuG-CSF, lenograstim). Peripheral blood progenitor cells (PBPC) were harvested by leukaphereses during early recovery from chemotherapy-induced aplasia. A total of 119 leukaphereses were performed. Median numbers of CD34+ cells and CFU-GM collected were 2.04 x 10(6)/kg and 2 9 x 10(4)/kg, respectively. There was a significant correlation between white cell count and number of CD34+ cells in the leukaphereses (P = 0.0001, r2 = 0.41, n = 104). A strict correlation between the number of CD34+ cells and CFU-GM in the leukapheretic product (P = 0.0001, r2 = 0.39, n = 110) was observed. 21% of evaluable patients (6/29) achieved a complete cytogenetic remission in the leukapheretic product and the other four patients achieved a major cytogenetic response for an overall response of 35% (10/22 patients). To date, 16 patients have been autografted and are alive. Five of them are Ph1-negative (three patients) or partially Ph1-negative (two patients). In conclusion, despite the high-risk characteristics of this study population, Ph1-negative PBPC were successfully mobilized in more than one-quarter of patients using a chemotherapy plus rhuG-CSF regimen. The importance of this achievement is increased by the current lack of other practical methods of rescuing Ph-negative cells in such patients.

Bioconductor: open software development for computational biology and bioinformatics

The Bioconductor project is an initiative for the collaborative creation of extensible software for computational biology and bioinformatics. The goals of the project include: fostering collaborative development and widespread use of innovative software, reducing barriers to entry into interdisciplinary scientific research, and promoting the achievement of remote reproducibility of research results. We describe details of our aims and methods, identify current challenges, compare Bioconductor to other open bioinformatics projects, and provide working examples.

Removal of the phage-shock protein PspB causes reduction of virulence in Salmonella enterica serovar Typhimurium independently of NRAMP1

The phage-shock protein (Psp) system is believed to manage membrane stress in all Enterobacteriaceae and has recently emerged as being important for virulence in several pathogenic species of this phylum. The core of the Psp system consists of the pspA–D operon and the distantly located pspG gene. In Salmonella enterica serovar Typhimurium (S. Typhimurium), it has recently been reported that PspA is essential for systemic infection of mice, but only in NRAMP1+ mice, signifying that attenuation is related to coping with divalent cation starvation in the intracellular environment. In the present study, we investigated the contribution of individual psp genes to virulence of S. Typhimurium. Interestingly, deletion of the whole pspA–D set of genes caused attenuation in both NRAMP1+ and NRAMP1− mice, indicating that one or more of the psp genes contribute to virulence independently of NRAMP1 expression in the host. Investigations of single gene mutants showed that knock out of pspB reduced virulence in both types of mice, while deletion of pspA only caused attenuation in NRAMP1+ mice, and deletion of pspD had a minor effect in NRAMP1− mice, while deletions of either pspC or pspG did not affect virulence. Experiments addressed at elucidating the role of PspB in virulence revealed that PspB is dispensable for uptake to and intracellular replication in cultured macrophages and resistance to complement-induced killing. Furthermore, the Psp system of S. Typhimurium was dispensable during pIV-induced secretin stress. In conclusion, our results demonstrate that removal of PspB reduces virulence in S. Typhimurium independently of host NRAMP1 expression, demonstrating that PspB has roles in intra-host survival distinct from the reported contributions of PspA.