Interpreting 16S metagenomic data without clustering to achieve sub-OTU resolution

The standard approach to analyzing 16S tag sequence data, which relies on clustering reads by sequence similarity into Operational Taxonomic Units (OTUs), underexploits the accuracy of modern sequencing technology. We present a clustering-free approach to multi-sample Illumina datasets that can identify independent bacterial subpopulations regardless of the similarity of their 16S tag sequences. Using published data from a longitudinal time-series study of human tongue microbiota, we are able to resolve within standard 97% similarity OTUs up to 20 distinct subpopulations, all ecologically distinct but with 16S tags differing by as little as 1 nucleotide (99.2% similarity). A comparative analysis of oral communities of two cohabiting individuals reveals that most such subpopulations are shared between the two communities at 100% sequence identity, and that dynamical similarity between subpopulations in one host is strongly predictive of dynamical similarity between the same subpopulations in the other host. Our method can also be applied to samples collected in cross-sectional studies and can be used with the 454 sequencing platform. We discuss how the sub-OTU resolution of our approach can provide new insight into factors shaping community assembly.Comment: Updated to match the published version. 12 pages, 5 figures + supplement. Significantly revised for clarity, references added, results not change

Who Watches the Watchmen? An Appraisal of Benchmarks for Multiple Sequence Alignment

Multiple sequence alignment (MSA) is a fundamental and ubiquitous technique in bioinformatics used to infer related residues among biological sequences. Thus alignment accuracy is crucial to a vast range of analyses, often in ways difficult to assess in those analyses. To compare the performance of different aligners and help detect systematic errors in alignments, a number of benchmarking strategies have been pursued. Here we present an overview of the main strategies--based on simulation, consistency, protein structure, and phylogeny--and discuss their different advantages and associated risks. We outline a set of desirable characteristics for effective benchmarking, and evaluate each strategy in light of them. We conclude that there is currently no universally applicable means of benchmarking MSA, and that developers and users of alignment tools should base their choice of benchmark depending on the context of application--with a keen awareness of the assumptions underlying each benchmarking strategy.Comment: Revie

Effect of Dietary Components on Larval Life History Characteristics in the Medfly (Ceratitis capitata: Diptera, Tephritidae)

Background: The ability to respond to heterogenous nutritional resources is an important factor in the adaptive radiation of insects such as the highly polyphagous Medfly. Here we examined the breadth of the Medfly’s capacity to respond to different developmental conditions, by experimentally altering diet components as a proxy for host quality and novelty. Methodology/Principal Findings: We tested responses of larval life history to diets containing protein and carbohydrate components found in and outside the natural host range of this species. A 40% reduction in the quantity of protein caused a significant increase in egg to adult mortality by 26.5%±6% in comparison to the standard baseline diet. Proteins and carbohydrates had differential effects on larval versus pupal development and survival. Addition of a novel protein source, casein (i.e. milk protein), to the diet increased larval mortality by 19.4%±3% and also lengthened the duration of larval development by 1.93±0.5 days in comparison to the standard diet. Alteration of dietary carbohydrate, by replacing the baseline starch with simple sugars, increased mortality specifically within the pupal stage (by 28.2%±8% and 26.2%±9% for glucose and maltose diets, respectively). Development in the presence of the novel carbohydrate lactose (milk sugar) was successful, though on this diet there was a decrease of 29.8±1.6 µg in mean pupal weight in comparison to pupae reared on the baseline diet. Conclusions: The results confirm that laboratory reared Medfly retain the ability to survive development through a wide range of fluctuations in the nutritional environment. We highlight new facets of the responses of different stages of holometabolous life histories to key dietary components. The results are relevant to colonisation scenarios and key to the biology of this highly invasive species

Evolution of Th2 responses : Characterization of IL-4/13 in sea bass (Dicentrarchus labrax L.) and studies of expression and biological activity

Acknowledgements This research was funded by the European Commission under the 7th Framework Programme for Research and Technological Development (FP7) of the European Union (Grant Agreement 311993 TARGETFISH). T.W. received funding from the MASTS pooling initiative (The Marine Alliance for Science and Technology for Scotland). MASTS is funded by the Scottish Funding Council (grant reference number HR09011) and contributing institutions.Peer reviewedPublisher PD

Evaluation of a Bayesian inference network for ligand-based virtual screening

Background Bayesian inference networks enable the computation of the probability that an event will occur. They have been used previously to rank textual documents in order of decreasing relevance to a user-defined query. Here, we modify the approach to enable a Bayesian inference network to be used for chemical similarity searching, where a database is ranked in order of decreasing probability of bioactivity. Results Bayesian inference networks were implemented using two different types of network and four different types of belief function. Experiments with the MDDR and WOMBAT databases show that a Bayesian inference network can be used to provide effective ligand-based screening, especially when the active molecules being sought have a high degree of structural homogeneity; in such cases, the network substantially out-performs a conventional, Tanimoto-based similarity searching system. However, the effectiveness of the network is much less when structurally heterogeneous sets of actives are being sought. Conclusion A Bayesian inference network provides an interesting alternative to existing tools for ligand-based virtual screening

The genome and transcriptome of Trichormus sp NMC-1: insights into adaptation to extreme environments on the Qinghai-Tibet Plateau

The Qinghai-Tibet Plateau (QTP) has the highest biodiversity for an extreme environment worldwide, and provides an ideal natural laboratory to study adaptive evolution. In this study, we generated a draft genome sequence of cyanobacteria Trichormus sp. NMC-1 in the QTP and performed whole transcriptome sequencing under low temperature to investigate the genetic mechanism by which T. sp. NMC-1 adapted to the specific environment. Its genome sequence was 5.9 Mb with a G+C content of 39.2% and encompassed a total of 5362 CDS. A phylogenomic tree indicated that this strain belongs to the Trichormus and Anabaena cluster. Genome comparison between T. sp. NMC-1 and six relatives showed that functionally unknown genes occupied a much higher proportion (28.12%) of the T. sp. NMC-1 genome. In addition, functions of specific, significant positively selected, expanded orthogroups, and differentially expressed genes involved in signal transduction, cell wall/membrane biogenesis, secondary metabolite biosynthesis, and energy production and conversion were analyzed to elucidate specific adaptation traits. Further analyses showed that the CheY-like genes, extracellular polysaccharide and mycosporine-like amino acids might play major roles in adaptation to harsh environments. Our findings indicate that sophisticated genetic mechanisms are involved in cyanobacterial adaptation to the extreme environment of the QTP

Evolutionary Toggling of Vpx/Vpr Specificity Results in Divergent Recognition of the Restriction Factor SAMHD1

SAMHD1 is a host restriction factor that blocks the ability of lentiviruses such as HIV-1 to undergo reverse transcription in myeloid cells and resting T-cells. This restriction is alleviated by expression of the lentiviral accessory proteins Vpx and Vpr (Vpx/Vpr), which target SAMHD1 for proteasome-mediated degradation. However, the precise determinants within SAMHD1 for recognition by Vpx/Vpr remain unclear. Here we show that evolution of Vpx/Vpr in primate lentiviruses has caused the interface between SAMHD1 and Vpx/Vpr to alter during primate lentiviral evolution. Using multiple HIV-2 and SIV Vpx proteins, we show that Vpx from the HIV-2 and SIVmac lineage, but not Vpx from the SIVmnd2 and SIVrcm lineage, require the C-terminus of SAMHD1 for interaction, ubiquitylation, and degradation. On the other hand, the N-terminus of SAMHD1 governs interactions with Vpx from SIVmnd2 and SIVrcm, but has little effect on Vpx from HIV-2 and SIVmac. Furthermore, we show here that this difference in SAMHD1 recognition is evolutionarily dynamic, with the importance of the N- and C-terminus for interaction of SAMHD1 with Vpx and Vpr toggling during lentiviral evolution. We present a model to explain how the head-to-tail conformation of SAMHD1 proteins favors toggling of the interaction sites by Vpx/Vpr during this virus-host arms race. Such drastic functional divergence within a lentiviral protein highlights a novel plasticity in the evolutionary dynamics of viral antagonists for restriction factors during lentiviral adaptation to its hosts. © 2013 Fregoso et al

Conserved presence of G-quadruplex forming sequences in the Long Terminal Repeat Promoter of Lentiviruses

G-quadruplexes (G4s) are secondary structures of nucleic acids that epigenetically regulate cellular processes. In the human immunodeficiency lentivirus 1 (HIV-1), dynamic G4s are located in the unique viral LTR promoter. Folding of HIV-1 LTR G4s inhibits viral transcription; stabilization by G4 ligands intensifies this effect. Cellular proteins modulate viral transcription by inducing/unfolding LTR G4s. We here expanded our investigation on the presence of LTR G4s to all lentiviruses. G4s in the 5'-LTR U3 region were completely conserved in primate lentiviruses. A G4 was also present in a cattle-infecting lentivirus. All other non-primate lentiviruses displayed hints of less stable G4s. In primate lentiviruses, the possibility to fold into G4s was highly conserved among strains. LTR G4 sequences were very similar among phylogenetically related primate viruses, while they increasingly differed in viruses that diverged early from a common ancestor. A strong correlation between primate lentivirus LTR G4s and Sp1/NF\u3baB binding sites was found. All LTR G4s folded: their complexity was assessed by polymerase stop assay. Our data support a role of the lentiviruses 5'-LTR G4 region as control centre of viral transcription, where folding/unfolding of G4s and multiple recruitment of factors based on both sequence and structure may take place

Development of the preterm gut microbiome in twins at risk of necrotising enterocolitis and sepsis

The preterm gut microbiome is a complex dynamic community influenced by genetic and environmental factors and is implicated in the pathogenesis of necrotising enterocolitis (NEC) and sepsis. We aimed to explore the longitudinal development of the gut microbiome in preterm twins to determine how shared environmental and genetic factors may influence temporal changes and compared this to the expressed breast milk (EBM) microbiome. Stool samples (n = 173) from 27 infants (12 twin pairs and 1 triplet set) and EBM (n = 18) from 4 mothers were collected longitudinally. All samples underwent PCR-DGGE (denaturing gradient gel electrophoresis) analysis and a selected subset underwent 454 pyrosequencing. Stool and EBM shared a core microbiome dominated by Enterobacteriaceae, Enterococcaceae, and Staphylococcaceae. The gut microbiome showed greater similarity between siblings compared to unrelated individuals. Pyrosequencing revealed a reduction in diversity and increasing dominance of Escherichia sp. preceding NEC that was not observed in the healthy twin. Antibiotic treatment had a substantial effect on the gut microbiome, reducing Escherichia sp. and increasing other Enterobacteriaceae. This study demonstrates related preterm twins share similar gut microbiome development, even within the complex environment of neonatal intensive care. This is likely a result of shared genetic and immunomodulatory factors as well as exposure to the same maternal microbiome during birth, skin contact and exposure to EBM. Environmental factors including antibiotic exposure and feeding are additional significant determinants of community structure, regardless of host genetics

Limited congruence exhibited across microbial, meiofaunal and macrofaunal benthic assemblages in a heterogeneous coastal environment

One of the most common approaches for investigating the ecology of spatially complex environments is to examine a single biotic assemblage present, such as macroinvertebrates. Underlying this approach are assumptions that sampled and unsampled taxa respond similarly to environmental gradients and exhibit congruence across different sites. These assumptions were tested for five benthic groups of various sizes (archaea, bacteria, microbial eukaryotes/protists, meiofauna and macrofauna) in Plymouth Sound, a harbour with many different pollution sources. Sediments varied in granulometry, hydrocarbon and trace metal concentrations. Following variable reduction, canonical correspondence analysis did not identify any associations between sediment characteristics and assemblage composition of archaea or macrofauna. In contrast, variation in bacteria was associated with granulometry, trace metal variations and bioturbation (e.g. community bioturbation potential). Protists varied with granulometry, hydrocarbon and trace metal predictors. Meiofaunal variation was associated with hydrocarbon and bioturbation predictors. Taxon turnover between sites varied with only three out of 10 group pairs showing congruence (meiofauna-protists, meiofauna-macrofauna and protists-macrofauna). While our results support using eukaryotic taxa as proxies for others, the lack of congruence suggests caution should be applied to inferring wider indicator or functional interpretations from studies of a single biotic assemblage