Species Delimitation and Conservation in Taxonomically Challenging Lineages: The Case of Two Clades of Capurodendron (Sapotaceae) in Madagascar.

Capurodendron is the largest endemic genus of plants from Madagascar, with around 76% of its species threatened by deforestation and illegal logging. However, some species are not well circumscribed and many of them remain undescribed, impeding a confident evaluation of their conservation status. Here we focus on taxa delimitation and conservation of two species complexes within Capurodendron: the Arid and Western complexes, each containing undescribed morphologies as well as intermediate specimens alongside well-delimited taxa. To solve these taxonomic issues, we studied 381 specimens morphologically and selected 85 of them to obtain intergenic, intronic, and exonic protein-coding sequences of 794 nuclear genes and 227 microsatellite loci. These data were used to test species limits and putative hybrid patterns using different approaches such as phylogenies, PCA, structure analyses, heterozygosity level, FST, and ABBA-BABA tests. The potential distributions were furthermore estimated for each inferred species. The results show that the Capurodendron Western Complex contains three well-delimited species, C. oblongifolium, C. perrieri, and C. pervillei, the first two hybridizing sporadically with the last and producing morphologies similar to, but genetically distinct from C. pervillei. The Arid Complex shows a more intricate situation, as it contains three species morphologically well-delimited but genetically intermixed. Capurodendron mikeorum nom. prov. is shown to be an undescribed species with a restricted distribution, while C. androyense and C. mandrarense have wider and mostly sympatric distributions. Each of the latter two species contains two major genetic pools, one showing interspecific admixture in areas where both taxa coexist, and the other being less admixed and comprising allopatric populations having fewer contacts with the other species. Only two specimens out of 172 showed clear genetic and morphological signals of recent hybridization, while all the others were morphologically well-delimited, independent of their degree of genetic admixture. Hybridization between Capurodendron androyense and C. microphyllum, the sister species of the Arid Complex, was additionally detected in areas where both species coexist, producing intermediate morphologies. Among the two complexes, species are well-defined morphologically with the exception of seven specimens (1.8%) displaying intermediate patterns and genetic signals compatible with a F1 hybridization. A provisional conservation assessment for each species is provided.

Dietary thiols accelerate aging of C. elegans.

Glutathione (GSH) is the most abundant cellular antioxidant. As reactive oxygen species (ROS) are widely believed to promote aging and age-related diseases, and antioxidants can neutralize ROS, it follows that GSH and its precursor, N-acetyl cysteine (NAC), are among the most popular dietary supplements. However, the long- term effects of GSH or NAC on healthy animals have not been thoroughly investigated. We employed C. elegans to demonstrate that chronic administration of GSH or NAC to young or aged animals perturbs global gene expression, inhibits skn-1-mediated transcription, and accelerates aging. In contrast, limiting the consumption of dietary thiols, including those naturally derived from the microbiota, extended lifespan. Pharmacological GSH restriction activates the unfolded protein response and increases proteotoxic stress resistance in worms and human cells. It is thus advantageous for healthy individuals to avoid excessive dietary antioxidants and, instead, rely on intrinsic GSH biosynthesis, which is fine-tuned to match the cellular redox status and to promote homeostatic ROS signaling.

New genetic markers for Sapotaceae phylogenomics: More than 600 nuclear genes applicable from family to population levels.

Some tropical plant families, such as the Sapotaceae, have a complex taxonomy, which can be resolved using Next Generation Sequencing (NGS). For most groups however, methodological protocols are still missing. Here we identified 531 monocopy genes and 227 Short Tandem Repeats (STR) markers and tested them on Sapotaceae using target capture and NGS. The probes were designed using two genome skimming samples from Capurodendron delphinense and Bemangidia lowryi, both from the Tseboneae tribe, as well as the published Manilkara zapota transcriptome from the Sapotoideae tribe. We combined our probes with 261 additional ones previously published and designed for the entire angiosperm group. On a total of 792 low-copy genes, 638 showed no signs of paralogy and were used to build a phylogeny of the family with 231 individuals from all main lineages. A highly supported topology was obtained at high taxonomic ranks but also at the species level. This phylogeny revealed the existence of more than 20 putative new species. Single nucleotide polymorphisms (SNPs) extracted from the 638 genes were able to distinguish lineages within a species complex and to highlight geographical structuration. STR were recovered efficiently for the species used as reference (C. delphinense) but the recovery rate decreased dramatically with the phylogenetic distance to the focal species. Altogether, the new loci will help reaching a sound taxonomic understanding of the family Sapotaceae for which many circumscriptions and relationships are still debated, at the species, genus and tribe levels.

Sharing Programming Resources Between Bio* Projects.

Open-source software encourages computer programmers to reuse software components written by others. In evolutionary bioinformatics, open-source software comes in a broad range of programming languages, including C/C++, Perl, Python, Ruby, Java, and R. To avoid writing the same functionality multiple times for different languages, it is possible to share components by bridging computer languages and Bio* projects, such as BioPerl, Biopython, BioRuby, BioJava, and R/Bioconductor.In this chapter, we compare the three principal approaches for sharing software between different programming languages: by remote procedure call (RPC), by sharing a local "call stack," and by calling program to programs. RPC provides a language-independent protocol over a network interface; examples are SOAP and Rserve. The local call stack provides a between-language mapping, not over the network interface but directly in computer memory; examples are R bindings, RPy, and languages sharing the Java virtual machine stack. This functionality provides strategies for sharing of software between Bio* projects, which can be exploited more often.Here, we present cross-language examples for sequence translation and measure throughput of the different options. We compare calling into R through native R, RSOAP, Rserve, and RPy interfaces, with the performance of native BioPerl, Biopython, BioJava, and BioRuby implementations and with call stack bindings to BioJava and the European Molecular Biology Open Software Suite (EMBOSS).In general, call stack approaches outperform native Bio* implementations, and these, in turn, outperform "RPC"-based approaches. To test and compare strategies, we provide a downloadable Docker container with all examples, tools, and libraries included.

Chromosome copy number variation in telomerized human bone marrow stromal cells; insights for monitoring safe ex-vivo expansion of adult stem cells.

Adult human bone marrow stromal cells (hBMSC) cultured for cell therapy require evaluation of potency and stability for safe use. Chromosomal aberrations upsetting genomic integrity in such cells have been contrastingly described as "Limited" or "Significant". Previously reported stepwise acquisition of a spontaneous neoplastic phenotype during three-year continuous culture of telomerized cells (hBMSC-TERT20) didn't alter a diploid karyotype measured by spectral karyotype analysis (SKY). Such screening may not adequately monitor abnormal and potentially tumorigenic hBMSC in clinical scenarios. We here used array comparative genomic hybridization (aCGH) to more stringently compare non-tumorigenic parental hBMSC-TERT strains with their tumorigenic subcloned populations. Confirmation of a known chromosome 9p21 microdeletion at locus CDKN2A/B, showed it also impinged upon the adjacent MTAP gene. Compared to reference diploid human fibroblast genomic DNA, the non-tumorigenic hBMSC-TERT4 cells had a copy number variation (CNV) in at least 14 independent loci. The pre-tumorigenic hBMSC-TERT20 cell strain had further CNV including 1q44 gain enhancing SMYD3 expression and 11q13.1 loss downregulating MUS81 expression. Bioinformatic analysis of gene products reflecting 11p15.5 CNV gain in tumorigenic hBMSC-TERT20 cells highlighted networks implicated in tumorigenic progression involving cell cycle control and mis-match repair. We provide novel biomarkers for prospective risk assessment of expanded stem cell cultures.

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.

Glycogen controls Caenorhabditis elegans lifespan and resistance to oxidative stress.

A high-sugar diet has been associated with reduced lifespan in organisms ranging from worms to mammals. However, the mechanisms underlying the harmful effects of glucose are poorly understood. Here we establish a causative relationship between endogenous glucose storage in the form of glycogen, resistance to oxidative stress and organismal aging in Caenorhabditis elegans. We find that glycogen accumulated on high dietary glucose limits C. elegans longevity. Glucose released from glycogen and used for NADPH/glutathione reduction renders nematodes and human hepatocytes more resistant against oxidative stress. Exposure to low levels of oxidants or genetic inhibition of glycogen synthase depletes glycogen stores and extends the lifespan of animals fed a high glucose diet in an AMPK-dependent manner. Moreover, glycogen interferes with low insulin signalling and accelerates aging of long-lived daf-2 worms fed a high glucose diet. Considering its extensive evolutionary conservation, our results suggest that glycogen metabolism might also have a role in mammalian aging.

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 species1. Studies of microbial co-occurrence in the gut have indicated community substructures that may reflect functional and metabolic interactions between cohabiting species2,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 H2 and CO2 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.

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

[This corrects the article on p. 362 in vol. 5, PMID: 25520736.].

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.

Identification and assembly of genomes and genetic elements in complex metagenomic samples without using reference genomes.

Most current approaches for analyzing metagenomic data rely on comparisons to reference genomes, but the microbial diversity of many environments extends far beyond what is covered by reference databases. De novo segregation of complex metagenomic data into specific biological entities, such as particular bacterial strains or viruses, remains a largely unsolved problem. Here we present a method, based on binning co-abundant genes across a series of metagenomic samples, that enables comprehensive discovery of new microbial organisms, viruses and co-inherited genetic entities and aids assembly of microbial genomes without the need for reference sequences. We demonstrate the method on data from 396 human gut microbiome samples and identify 7,381 co-abundance gene groups (CAGs), including 741 metagenomic species (MGS). We use these to assemble 238 high-quality microbial genomes and identify affiliations between MGS and hundreds of viruses or genetic entities. Our method provides the means for comprehensive profiling of the diversity within complex metagenomic samples.

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.

Towards in vivo regulon kinetics: PurR activation by 5-phosphoribosyl-α-1-pyrophosphate during purine depletion in Lactococcus lactis.

Short-term adaptation to changing environments relies on regulatory elements translating shifting metabolite concentrations into a specifically optimized transcriptome. So far the focus of analyses has been divided between regulatory elements identified in vivo and kinetic studies of small molecules interacting with the regulatory elements in vitro. Here we describe how in vivo regulon kinetics can describe a regulon through the effects of the metabolite controlling it, exemplified by temporal purine exhaustion in Lactococcus lactis. We deduced a causal relation between the pathway precursor 5-phosphoribosyl-α-1-pyrophosphate (PRPP) and individual mRNA levels, whereby unambiguous and homogeneous relations could be obtained for PurR regulated genes, thus linking a specific regulon to a specific metabolite. As PurR activates gene expression upon binding of PRPP, the pur mRNA curves reflect the in vivo kinetics of PurR PRPP binding and activation. The method singled out the xpt-pbuX operon as kinetically distinct, which was found to be caused by a guanine riboswitch whose regulation was overlaying the PurR regulation. Importantly, genes could be clustered according to regulatory mechanism and long-term consequences could be distinguished from transient changes--many of which would not be seen in a long-term adaptation to a new environment. The strategy outlined here can be adapted to analyse the individual effects of members from larger metabolomes in virtually any organism, for elucidating regulatory networks in vivo.

The effect of network biology on drug toxicology.

INTRODUCTION: The high failure rate of drug candidates due to toxicity, during clinical trials, is a critical issue in drug discovery. Network biology has become a promising approach, in this regard, using the increasingly large amount of biological and chemical data available and combining it with bioinformatics. With this approach, the assessment of chemical safety can be done across multiple scales of complexity from molecular to cellular and system levels in human health. Network biology can be used at several levels of complexity. AREAS COVERED: This review describes the strengths and limitations of network biology. The authors specifically assess this approach across different biological scales when it is applied to toxicity. EXPERT OPINION: There has been much progress made with the amount of data that is generated by various omics technologies. With this large amount of useful data, network biology has the opportunity to contribute to a better understanding of a drug's safety profile. The authors believe that considering a drug action and protein's function in a global physiological environment may benefit our understanding of the impact some chemicals have on human health and toxicity. The next step for network biology will be to better integrate differential and quantitative data.

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.

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.

Low-bandwidth and non-compute intensive remote identification of microbes from raw sequencing reads.

Cheap DNA sequencing may soon become routine not only for human genomes but also for practically anything requiring the identification of living organisms from their DNA: tracking of infectious agents, control of food products, bioreactors, or environmental samples. We propose a novel general approach to the analysis of sequencing data where a reference genome does not have to be specified. Using a distributed architecture we are able to query a remote server for hints about what the reference might be, transferring a relatively small amount of data. Our system consists of a server with known reference DNA indexed, and a client with raw sequencing reads. The client sends a sample of unidentified reads, and in return receives a list of matching references. Sequences for the references can be retrieved and used for exhaustive computation on the reads, such as alignment. To demonstrate this approach we have implemented a web server, indexing tens of thousands of publicly available genomes and genomic regions from various organisms and returning lists of matching hits from query sequencing reads. We have also implemented two clients: one running in a web browser, and one as a python script. Both are able to handle a large number of sequencing reads and from portable devices (the browser-based running on a tablet), perform its task within seconds, and consume an amount of bandwidth compatible with mobile broadband networks. Such client-server approaches could develop in the future, allowing a fully automated processing of sequencing data and routine instant quality check of sequencing runs from desktop sequencers. A web access is available at http://tapir.cbs.dtu.dk. The source code for a python command-line client, a server, and supplementary data are available at http://bit.ly/1aURxkc.

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.