Search

MotivationRecombinant protein production is a widely used technique in the biotechnology and biomedical industries, yet only a quarter of target proteins are soluble and can therefore be purified.ResultsWe have discovered that global structural flexibility, which can be modeled by normalised B-factors, accurately predicts the solubility of 12,216 recombinant proteins expressed in Escherichia coli. We have optimised B-factors, and derived a new set of values for solubility scoring that further improves prediction accuracy. We call this new predictor the ‘Solubility-Weighted Index’ (SWI). Importantly, SWI outperforms many existing protein solubility prediction tools. Furthermore, we have developed ‘SoDoPE’ (Soluble Domain for Protein Expression), a web interface that allows users to choose a protein region of interest for predicting and maximising both protein expression and solubility.AvailabilityThe SoDoPE web server and source code are freely available at https://tisigner.com/sodope and https://github.com/Gardner-BinfLab/TISIGNER-ReactJS, respectively. The code and data for reproducing our analysis can be found at https://github.com/Gardner-BinfLab/SoDoPE_paper2020.

CRISPR–Cas systems provide bacteria with adaptive immunity against bacteriophages1. However, DNA modification2,3, the production of anti-CRISPR proteins4,5 and potentially other strategies enable phages to evade CRISPR–Cas. Here, we discovered a Serratia jumbo phage that evades type I CRISPR–Cas systems, but is sensitive to type III immunity. Jumbo phage infection resulted in a nucleus-like structure enclosed by a proteinaceous phage shell—a phenomenon only reported recently for distantly related Pseudomonas phages6,7. All three native CRISPR–Cas complexes in Serratia—type I-E, I-F and III-A—were spatially excluded from the phage nucleus and phage DNA was not targeted. However, the type III-A system still arrested jumbo phage infection by targeting phage RNA in the cytoplasm in a process requiring Cas7, Cas10 and an accessory nuclease. Type III, but not type I, systems frequently targeted nucleus-forming jumbo phages that were identified in global viral sequence datasets. The ability to recognize jumbo phage RNA and elicit immunity probably contributes to the presence of both RNA- and DNA-targeting CRISPR–Cas systems in many bacteria1,8. Together, our results support the model that jumbo phage nucleus-like compartments serve as a barrier to DNA-targeting, but not RNA-targeting, defences, and that this phenomenon is widespread among jumbo phages. This paper identifies a Serratia jumbo phage that, on infection, leads to the formation of a nucleus-like structure that protects phage DNA from CRISPR–Cas defence systems. However, the phage is still susceptible to CRISPR–Cas RNA targeting in the cytoplasm.

In computational biology and other sciences, researchers are frequently faced with a choice between several computational methods for performing data analyses. Benchmarking studies aim to rigorously compare the performance of different methods using well-characterized benchmark datasets, to determine the strengths of each method or to provide recommendations regarding suitable choices of methods for an analysis. However, benchmarking studies must be carefully designed and implemented to provide accurate, unbiased, and informative results. Here, we summarize key practical guidelines and recommendations for performing high-quality benchmarking analyses, based on our experiences in computational biology., Comment: Minor updates

Experiments that are planned using accurate prediction algorithms will mitigate failures in recombinant protein production. We have developed TISIGNER (https://tisigner.com) with the aim of addressing technical challenges to recombinant protein production. We offer three web services, TIsigner (Translation Initiation coding region designer), SoDoPE (Soluble Domain for Protein Expression) and Razor, which are specialised in synonymous optimisation of recombinant protein expression, solubility and signal peptide analysis, respectively. Importantly, TIsigner, SoDoPE and Razor are linked, which allows users to switch between the tools when optimising genes of interest., Graphical Abstract Graphical AbstractOptimising recombinant protein production using TISIGNER.com

MotivationSignal peptides are responsible for protein transport and secretion and are ubiquitous to all forms of life. The annotation of signal peptides is important for understanding protein translocation and toxin secretion and evolution.ResultsHere we explore the features of these signal sequences from eukaryotic proteins. Strikingly, we find that the signal peptides from secretory toxins have common features across kingdoms, supporting the idea of horizontal gene transfer or convergence of toxin genes across kingdoms. We leverage these features to build Razor, a simple yet powerful tool specialised in identifying signal peptides from toxins using the first 23 N-terminal residues. We demonstrate the usability of Razor by analysing all the sequences reviewed by UniProt. Indeed, Razor is able to identify toxins using their N-terminal sequences only. Strikingly, we also discover that many defensive proteins across kingdoms harbour a toxin-like signal peptide; some of these defensive proteins have emerged through convergent evolution, e.g. defensin and defensin-like protein families, and phospholipase families. In sum, Razor uses an approach independent of homology search to identify novel and known toxin classes across species using N-terminal residues.Availability and implementationRazor is available as a web application (https://tisigner.com/razor) and a command-line tool (https://github.com/Gardner-BinfLab/Razor).

Proteins and non-coding RNAs are functional products of the genome that carry out the bulk of crucial cellular processes. With recent technological advances, researchers can sequence genomes in the thousands as well as probe for specific genomic activities of multiple species and conditions. These studies have identified thousands of potential proteins, RNAs and associated activities, however there are conflicting conclusions on the functional implications depending upon the burden of evidence researchers use, leading to diverse interpretations of which regions of the genome are “functional”. Here we investigate the association between gene functionality and genomic features, by comparing established functional protein-coding and non-coding genes to non-genic regions of the genome. We find that the strongest and most consistent association between functional genes and any genomic feature is evolutionary conservation and transcriptional activity. Other strongly associated features include sequence alignment statistics, such as maximum between-site covariation. We have also identified some concerns with 1,000 Genomes Project and Genome Aggregation Database SNP densities, as short non-coding RNAs tend to have greater than expected SNP densities. Our results demonstrate the importance of evolutionary conservation and transcription for sequence functionality, which should both be taken into consideration when differentiating between functional sequences and noise.

Planning experiments using accurate prediction algorithms could mitigate failures in recombinant protein production. We have developed TISIGNER.com with the aim of addressing the technical challenges in recombinant protein production. We offer two web services, TIsigner (Translation Initiation coding region designer) and SoDoPE (Soluble Domain for Protein Expression), which are specialised in prediction/optimisation of recombinant protein expression and solubility, respectively. Importantly, TIsigner and SoDoPE are linked, which allows users to switch between the tools when optimising their genes of interest.

Variants within the non-coding genome are frequently associated with phenotypes in genome-wide association studies. These non-coding regions may be involved in the regulation of gene expression, encode functional non-coding RNAs, or influence splicing and other cellular functions. We have curated a list of characterized non-coding human genome variants based on the published evidence that indicates phenotypic consequences of the variation. In order to minimize annotation errors, two curators have independently verified the supporting evidence for pathogenicity of each non-coding variant in the published literature. The database consists of 721 non-coding variants linked to the published literature describing the evidence of functional consequences. We have also sampled 7228 covariate-matched benign controls, that have a population frequency of over 5%, from the single nucleotide polymorphism database (dbSNP151) database. These were sampled controlling for potential confounding factors such as linkage with pathogenic variants, annotation type (untranslated region, intron, intergenic, etc.) and variant type (substitution or indel). The dataset presented here represents a curated repository, with a potential use for the training or evaluation of algorithms used in the prediction of non-coding variant functionality.Database URL: https://github.com/Gardner-BinfLab/ncVarDB.

Ngatiwai Trust Board. et al., The tuatara (Sphenodon punctatus)—the only living member of the reptilian order Rhynchocephalia (Sphenodontia), once widespread across Gondwana1,2—is an iconic species that is endemic to New Zealand2,3. A key link to the now-extinct stem reptiles (from which dinosaurs, modern reptiles, birds and mammals evolved), the tuatara provides key insights into the ancestral amniotes2,4. Here we analyse the genome of the tuatara, which—at approximately 5 Gb—is among the largest of the vertebrate genomes yet assembled. Our analyses of this genome, along with comparisons with other vertebrate genomes, reinforce the uniqueness of the tuatara. Phylogenetic analyses indicate that the tuatara lineage diverged from that of snakes and lizards around 250 million years ago. This lineage also shows moderate rates of molecular evolution, with instances of punctuated evolution. Our genome sequence analysis identifies expansions of proteins, non-protein-coding RNA families and repeat elements, the latter of which show an amalgam of reptilian and mammalian features. The sequencing of the tuatara genome provides a valuable resource for deep comparative analyses of tetrapods, as well as for tuatara biology and conservation. Our study also provides important insights into both the technical challenges and the cultural obligations that are associated with genome sequencing., N.J.G. acknowledges the support of Ngatiwai iwi, Allan Wilson Centre, University of Otago, New Zealand Department of Conservation, New Zealand Genomics and Illumina. J.I.A. was supported by CONICYT National Doctoral Scholarship No. 21130515. M.W. was supported by NIH grant R35 GM124827. Ensembl annotation was supported by the Wellcome Trust (WT108749/Z/15/Z) and the European Molecular Biology Laboratory.

DNA-based techniques are increasingly used for measuring the biodiversity (species presence, identity, abundance and community composition) of terrestrial and aquatic ecosystems. While there are numerous reviews of molecular methods and bioinformatic steps, there has been little consideration of the methods used to collect samples upon which these later steps are based. This represents a critical knowledge gap, as methodologically sound field sampling is the foundation for subsequent analyses. We reviewed field sampling methods used for metabarcoding studies of both terrestrial and freshwater ecosystem biodiversity over a nearly three-year period (n = 75). We found that 95% (n = 71) of these studies used subjective sampling methods and inappropriate field methods and/or failed to provide critical methodological information. It would be possible for researchers to replicate only 5% of the metabarcoding studies in our sample, a poorer level of reproducibility than for ecological studies in general. Our findings suggest greater attention to field sampling methods, and reporting is necessary in eDNA-based studies of biodiversity to ensure robust outcomes and future reproducibility. Methods must be fully and accurately reported, and protocols developed that minimize subjectivity. Standardization of sampling protocols would be one way to help to improve reproducibility and have additional benefits in allowing compilation and comparison of data from across studies.

Recombinant protein production is a key process in generating proteins of interest in the pharmaceutical industry and biomedical research. However, about 50% of recombinant proteins fail to be expressed in a variety of host cells. Here we show that the accessibility of translation initiation sites modelled using the mRNA base-unpairing across the Boltzmann’s ensemble significantly outperforms alternative features. This approach accurately predicts the successes or failures of expression experiments, which utilised Escherichia coli cells to express 11,430 recombinant proteins from over 189 diverse species. On this basis, we develop TIsigner that uses simulated annealing to modify up to the first nine codons of mRNAs with synonymous substitutions. We show that accessibility captures the key propensity beyond the target region (initiation sites in this case), as a modest number of synonymous changes is sufficient to tune the recombinant protein expression levels. We build a stochastic simulation model and show that higher accessibility leads to higher protein production and slower cell growth, supporting the idea of protein cost, where cell growth is constrained by protein circuits during overexpression., Author summary Recombinant proteins are widely used as therapeutics, such as vaccines, monoclonal antibodies, hormones and enzymes. However, the success rate of recombinant protein production is about 50%. To address this problem, we propose optimising the unpairing propensities of nucleotides around translation initiation sites using a thermodynamic quantity called mRNA accessibility. Our study shows that this method is generalisable across prokaryotic and eukaryotic expression hosts. Importantly, we validated this method using laboratory experiments and computational modelling. Furthermore, we propose a low cost technique to tune protein expression by engineering minimal changes to genes of interest through our web application (https://tisigner.com/tisigner).

The tuatara (Sphenodon punctatus), the only living member of the archaic reptilian order Rhynchocephalia (Sphenodontia) once widespread across Gondwana, is an iconic and enigmatic terrestrial vertebrate endemic to New Zealand. A key link to the now extinct stem reptiles from which dinosaurs, modern reptiles, birds and mammals evolved, the tuatara provides exclusive insights into the ancestral amniotes. The tuatara genome, at ∼5 Gbp, is among the largest vertebrate genomes assembled. Analysis of this genome and comparisons to other vertebrates reinforces the uniqueness of the tuatara. Phylogenetic analyses indicate tuatara diverged from the snakes and lizards ∼250 MYA. This lineage also shows moderate rates of molecular evolution, with instances of punctuated evolution. Genome sequence analysis identifies expansions of protein, non-protein-coding RNA families, and repeat elements, the latter of which show an extraordinary amalgam of reptilian and mammalian features. Sequencing of this genome provides a valuable resource for deep comparative analyses of tetrapods, as well as for tuatara biology and conservation. It also provides important insights into both the technical challenges and the cultural obligations associated with genome sequencing.

SUMMARYCRISPR-Cas systems provide bacteria with adaptive immunity against bacteriophages1. However, DNA modification2,3, the production of anti-CRISPR proteins4,5 and potentially other strategies enable phages to evade CRISPR-Cas. Here we discovered a Serratia jumbophage that evaded type I CRISPR-Cas systems, but was sensitive to type III immunity. Jumbophage infection resulted in a nucleus-like structure enclosed by a proteinaceous phage shell – a phenomenon only reported recently for distantly related Pseudomonas phages6,7. All three native CRISPR-Cas complexes in Serratia – type I-E, I-F and III-A – were spatially excluded from the phage nucleus and phage DNA was not targeted. However, the type III-A system still arrested jumbophage infection by targeting phage RNA in the cytoplasm in a process requiring Cas7, Cas10 and an accessory nuclease. Type III, but not type I, systems frequently targeted nucleus-forming jumbophages that were identified in global viral sequence datasets. These findings explain why many bacteria harbour both RNA- and DNA-targeting CRISPR-Cas systems1,8. Together, our results indicate that jumbophage nucleus-like compartments serve as a barrier to DNA-targeting, but not RNA-targeting defences, and that this phenomenon is widespread amongst jumbophages.

Recombinant protein production is a widely used technique, yet half of these experiments fail at the expression phase. Failures are expected for 9difficult-to-express9 proteins, but for others, codon bias and mRNA folding have been proposed explanations for variable protein abundance. We question how practical these features are for solving protein expression failures. We discover that the energetics of RNA structure ensembles, which models the 9accessibility9 of translation initiation sites, show the strongest correlation with protein abundance across species. Importantly, accessibility outperforms other features when predicting the outcomes of 11,430 recombinant protein production experiments in Escherichia coli. Furthermore, we show that protein level is tunable by synonymous codon changes of the first few codons that alter accessibility.

Understanding how new genes originate and integrate into cellular networks is key to understanding evolution. Bacteria present unique opportunities for both the natural history and experimental study of gene origins, due to their large effective population sizes, rapid generation times, and ease of genetic manipulation. Bacterial small non-coding RNAs (sRNAs), in particular, many of which operate through a simple antisense regulatory logic, may serve as tractable models for exploring processes of gene origin and adaptation. Understanding how and on what timescales these regulatory molecules arise has important implications for understanding the evolution of bacterial regulatory networks, in particular, for the design of comparative studies of sRNA function. Here, we introduce relevant concepts from evolutionary biology and review recent work that has begun to shed light on the timescales and processes through which non-functional transcriptional noise is co-opted to provide regulatory functions. We explore possible scenarios for sRNA origin, focusing on the co-option, or exaptation, of existing genomic structures which may provide protected spaces for sRNA evolution.

Campylobacter jejuni is the most common cause of bacterial diarrheal disease in the world. Clinical outcomes of infection can range from asymptomatic infection to life-threatening extraintestinal infections. This variability in outcomes for infected patients has raised questions as to whether genetic differences between C. jejuni isolates contribute to their likelihood of causing severe disease. In this study, we compare the genomes of ten C. jejuni isolates that were implicated in extraintestinal infections with reference gastrointestinal isolates, in order to identify unusual patterns of sequence variation associated with infection outcome. We identified a collection of genes that display a higher burden of uncommon mutations in invasive isolates compared with gastrointestinal close relatives, including some that have been previously linked to virulence and invasiveness in C. jejuni. Among the top genes identified were mreB and pgp1, which are both involved in determining cell shape. Electron microscopy confirmed morphological differences in isolates carrying unusual sequence variants of these genes, indicating a possible relationship between extraintestinal infection and changes in cell morphology.

Metagenomic and meta-barcode DNA sequencing has rapidly become a widely-used technique for investigating a range of questions, particularly related to health and environmental monitoring. There has also been a proliferation of bioinformatic tools for analysing metagenomic and amplicon datasets, which makes selecting adequate tools a significant challenge. A number of benchmark studies have been undertaken; however, these can present conflicting results. In order to address this issue we have applied a robustZ-score ranking procedure and a network meta-analysis method to identify software tools that are consistently accurate for mapping DNA sequences to taxonomic hierarchies. Based upon these results we have identified some tools and computational strategies that produce robust predictions.

Motivation The aim of this study is to assess the performance of RNA–RNA interaction prediction tools for all domains of life. Results Minimum free energy (MFE) and alignment methods constitute most of the current RNA interaction prediction algorithms. The MFE tools that include accessibility (i.e. RNAup, IntaRNA and RNAplex) to the final predicted binding energy have better true positive rates (TPRs) with a high positive predictive values (PPVs) in all datasets than other methods. They can also differentiate almost half of the native interactions from background. The algorithms that include effects of internal binding energies to their model and alignment methods seem to have high TPR but relatively low associated PPV compared to accessibility based methods. Availability and Implementation We shared our wrapper scripts and datasets at Github (github.com/UCanCompBio/RNA_Interactions_Benchmark). All parameters are documented for personal use. Supplementary information Supplementary data are available at Bioinformatics online.

The tuatara (Sphenodon punctatus)-the only living member of the reptilian order Rhynchocephalia (Sphenodontia), once widespread across Gondwana

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

The reactions of functional molecules like proteins and RNAs to mutation affect both host cell viability and biomolecular evolution. These molecules are considered robust if function is maintained despite mutations. Proteins and RNAs have different structural and functional characteristics that affect their robustness, and to date, comparisons between them have been theoretical. In this work, we test the relative mutational robustness of RNA and protein pairs using three approaches: evolutionary, structural, and functional. We compare the nucleotide diversities of functional RNAs with those of matched proteins. Across different levels of conservation, we found the nucleotide-level variations between the biomolecules largely overlapped, with proteins generally supporting more variation than matched RNAs. We then directly tested the robustness of the protein and RNA pairs with in vitro and in silico mutagenesis of their respective genes. The in silico experiments showed that proteins and RNAs reacted similarly to point mutations and insertions or deletions, yet proteins are slightly more robust on average than RNAs. In vitro, mutated fluorescent RNAs retained greater levels of function than the proteins. Overall this suggests that proteins and RNAs have remarkably similar degrees of robustness, with the average protein having moderately higher robustness than RNA as a group.Significance StatementThe ability of proteins and non-coding RNAs to maintain function despite mutations in their respective genes is known as mutational robustness. Robustness impacts how molecules maintain and change phenotypes, which has a bearing on the evolution and the origin of life as well as influencing modern biotechnology. Both protein and RNA have mechanisms that allow them to absorb DNA-level changes. Proteins have a redundant genetic code and non-coding RNAs can maintain structure and function through flexible base-pairing possibilities. The few theoretical treatments comparing protein and RNA robustness differ in their conclusions. In this experimental comparison of protein and RNA, we find that they have remarkably similar degrees of overall genetic robustness.

A core question in evolutionary biology is whether convergent phenotypic evolution is driven by convergent molecular changes in proteins or regulatory regions. We combined phylogenomic, developmental, and epigenomic analysis of 11 new genomes of paleognathous birds, including an extinct moa, to show that convergent evolution of regulatory regions, more so than protein-coding genes, is prevalent among developmental pathways associated with independent losses of flight. A Bayesian analysis of 284,001 conserved noncoding elements, 60,665 of which are corroborated as enhancers by open chromatin states during development, identified 2355 independent accelerations along lineages of flightless paleognaths, with functional consequences for driving gene expression in the developing forelimb. Our results suggest that the genomic landscape associated with morphological convergence in ratites has a substantial shared regulatory component.

The relative roles of regulatory and protein evolution in the origin and loss of convergent phenotypic traits is a core question in evolutionary biology. Here we combine phylogenomic, epigenomic and developmental data to show that convergent evolution of regulatory regions, but not protein-coding genes, is associated with flightlessness in palaeognathous birds, a classic example of a convergent phenotype. Eleven new genomes, including a draft genome from an extinct moa, resolve palaeognath phylogeny and show that the incidence of independent, convergent accelerations among 284,000 conserved non-exonic elements is significantly more frequent in ratites than other bird lineages. Ratite-specific acceleration of conserved regions and measures of open chromatin across eight tissues in the developing chick identify candidate regulatory regions that may have modified or lost function in ratites. Enhancer activity assays conducted in the early developing chicken forelimb confirm that volant versions of a conserved element in the first intron of the TEAD1 gene display conserved enhancer activity, whereas an accelerated flightless version fails to drive reporter gene expression. Our results show that convergent molecular changes associated with loss of flight are largely regulatory in nature.

It is now over 10 years since the first avian genome [International Chicken Genome Sequencing Consortium, 2004] and the first complete avian karyotype [Masabanda et al., 2004] were both published; however, until 2014, avian cytogenetics has focused heavily on descriptive studies [e.g. Griffin et al., 2007, 2008; Skinner et al., 2009; Volker et al., 2010] with less attention to its functional relevance. Last year, however, saw 2 landmark efforts in the chromosomal studies of birds: a special issue of Chromosome Research in April and the announcement of recently completed sequences of multiple new avian genomes in Science and the BMC journals (taking the total number sequenced to over 50) in December. Studying the chromosomes of birds is, perhaps for the first time, telling us more about avian biology, function and evolution than it ever has... Conclusions. The most recent advances in avian cytogenetics have culminated in great promise not only for the study of bird karyotypes, but also for providing insight into the mechanisms of chromosome evolution in general. New avenues for investigation include gene regulation; for instance, it will become necessary to map accurately the physical location of polyadenylation and transcription start sites, important reference points that define promoters and post-transcriptional regulation. It will also become possible to sequence full-length transcripts, to allow accurate identification of alternate splicing events and their controlling elements. The ENCODE (Encyclopedia of DNA Elements) project has helped to define functional elements of the human genome, including those aforementioned as well as other chromatin signals, e.g. active chromatin, enhancers, insulators, methylation domains, etc. An effort of agENCODE is underway to include agriculturally important birds such as chicken, turkey, duck, quail, and perhaps ostrich. The study of cytogenetics will be essential here in helping to define higher-order structures in nuclear organization that show regulatory interactions within and between chromosomes. Finally, reconstruction of evolutionary events allows us to study genome organization and function not only in extant but, by extrapolation, in extinct species also. Reconstruction of avian-reptilian ancestral karyotypes will allow us to define chromosomal rearrangements in long-dead species that have captured the public imagination. Here be dragons!

Neisseria meningitidis (meningococcus) can cause meningococcal disease, a rapidly progressing and often fatal disease that can occur in previously healthy children. Meningococci are found in healthy carriers, where they reside in the nasopharynx as commensals. While carriage is relatively common, invasive disease, associated with hypervirulent strains, is a comparatively rare event. The basis of increased virulence in some strains is not well understood. New Zealand suffered a protracted meningococcal disease epidemic, from 1991 to 2008. During this time, a household carriage study was carried out in Auckland: household contacts of index meningococcal disease patients were swabbed for isolation of carriage strains. In many households, healthy carriers harbored strains identical, as determined by laboratory typing, to the ones infecting the associated patient. We carried out more-detailed analyses of carriage and disease isolates from a select number of households. We found that isolates, although indistinguishable by laboratory typing methods and likely closely related, had many differences. We identified multiple genome variants and transcriptional differences between isolates. These studies enabled the identification of two new phase-variable genes. We also found that several carriage strains had lost their type IV pili and that this loss correlated with reduced tumor necrosis factor alpha (TNF-α) expression when cultured with epithelial cells. While nonpiliated meningococcal isolates have been previously found in carriage strains, this is the first evidence of an association between type IV pili from meningococci and a proinflammatory epithelial response. We also identified potentially important metabolic differences between carriage and disease isolates, including the sulfate assimilation pathway. IMPORTANCE Neisseria meningitidis causes meningococcal disease but is frequently carried in the throats of healthy individuals; the factors that determine whether invasive disease develops are not completely understood. We carried out detailed studies of isolates, collected from patients and their household contacts, to identify differences between commensal throat isolates and those that caused invasive disease. Though isolates were identical by laboratory typing methods, we uncovered many differences in their genomes, in gene expression, and in their interactions with host cells. In particular, we found that several carriage isolates had lost their type IV pili, a surprising finding since pili are often described as essential for colonization. However, loss of type IV pili correlated with reduced secretion of a proinflammatory cytokine, TNF-α, when meningococci were cocultured with human bronchial epithelial cells; hence, the loss of pili could provide an advantage to meningococci, by resulting in a dampened localized host immune response.

Emerging pathogens are a major threat to public health, however understanding how pathogens adapt to new niches remains a challenge. New methods are urgently required to provide functional insights into pathogens from the massive genomic data sets now being generated from routine pathogen surveillance for epidemiological purposes. Here, we measure the burden of atypical mutations in protein coding genes across independently evolved Salmonella enterica lineages, and use these as input to train a random forest classifier to identify strains associated with extraintestinal disease. Members of the species fall along a continuum, from pathovars which cause gastrointestinal infection and low mortality, associated with a broad host-range, to those that cause invasive infection and high mortality, associated with a narrowed host range. Our random forest classifier learned to perfectly discriminate long-established gastrointestinal and invasive serovars of Salmonella. Additionally, it was able to discriminate recently emerged Salmonella Enteritidis and Typhimurium lineages associated with invasive disease in immunocompromised populations in sub-Saharan Africa, and within-host adaptation to invasive infection. We dissect the architecture of the model to identify the genes that were most informative of phenotype, revealing a common theme of degradation of metabolic pathways in extraintestinal lineages. This approach accurately identifies patterns of gene degradation and diversifying selection specific to invasive serovars that have been captured by more labour-intensive investigations, but can be readily scaled to larger analyses., Author summary Researchers are now collecting a wealth of genomic data from bacterial pathogens, and this will continue to grow with the introduction of routine sequencing for disease surveillance. However, our ability to use this data to predict how changes in genome sequence lead to differences in disease is limited. Here, we have used machine learning to detect an enrichment in functionally significant mutations in genes associated with a shift in pathogenic niche. This approach captures convergence in functional outcomes that does not necessarily result in a convergence in sequence, facilitating the inclusion of rare variants of large effect in an analysis, and allowing for complex interactions between genes. We apply this approach to Salmonella, showing that we can detect changes associated with disease phenotype in emerging lineages associated with the HIV epidemic. This approach should be applicable to other bacterial species with lineages independently adapting to similar niches. We provide open-source implementations of both the predictive model, and the workflow used to build it.

Background The New Zealand collembolan genus Holacanthella contains the largest species of springtails (Collembola) in the world. Using Illumina technology we have sequenced and assembled a draft genome and transcriptome from Holacanthella duospinosa (Salmon). We have used this annotated assembly to investigate the genetic basis of a range of traits critical to the evolution of the Hexapoda, the phylogenetic position of H. duospinosa and potential horizontal gene transfer events. Results Our genome assembly was ~375 Mbp in size with a scaffold N50 of ~230 Kbp and sequencing coverage of ~180×. DNA elements, LTRs and simple repeats and LINEs formed the largest components and SINEs were very rare. Phylogenomics (370,877 amino acids) placed H. duospinosa within the Neanuridae. We recovered orthologs of the conserved sex determination genes thought to play a role in sex determination. Analysis of CpG content suggested the absence of DNA methylation, and consistent with this we were unable to detect orthologs of the DNA methyltransferase enzymes. The small subunit rRNA gene contained a possible retrotransposon. The Hox gene complex was broken over two scaffolds. For chemosensory ability, at least 15 and 18 ionotropic glutamate and gustatory receptors were identified, respectively. However, we were unable to identify any odorant receptors or their obligate co-receptor Orco. Twenty-three chitinase-like genes were identified from the assembly. Members of this multigene family may play roles in the digestion of fungal cell walls, a common food source for these saproxylic organisms. We also detected 59 and 96 genes that blasted to bacteria and fungi, respectively, but were located on scaffolds that otherwise contained arthropod genes. Conclusions The genome of H. duospinosa contains some unusual features including a Hox complex broken over two scaffolds, in a different manner to other arthropod species, a lack of odorant receptor genes and an apparent lack of environmentally responsive DNA methylation, unlike many other arthropods. Our detection of candidate horizontal gene transfer candidates confirms that this phenomenon is occurring across Collembola. These findings allow us to narrow down the regions of the arthropod phylogeny where key innovations have occurred that have facilitated the evolutionary success of Hexapoda. Electronic supplementary material The online version of this article (10.1186/s12864-017-4197-1) contains supplementary material, which is available to authorized users.

Legionella is a highly diverse genus of intracellular bacterial pathogens that cause Legionnaire’s disease (LD), an often severe form of pneumonia. Two L. micdadei sp. clinical isolates, obtained from patients hospitalized with LD from geographically distinct areas, were sequenced using PacBio SMRT cell technology, identifying incomplete phage regions, which may impact virulence.

RNA performs a diverse array of important functions across all cellular life. These functions include important roles in translation, building translational machinery and maturing messenger RNA. More recent discoveries include the miRNAs and bacterial sRNAs that regulate gene expression, the thermosensors, riboswitches and other cis-regulatory elements that help prokaryotes sense their environment and eukaryotic piRNAs that suppress transposition. However, there can be a long period between the initial discovery of a RNA and determining its function. We present a bioinformatic approach to characterize RNA motifs, which are critical components of many RNA structure-function relationships. These motifs can, in some instances, provide researchers with functional hypotheses for uncharacterized RNAs. Moreover, we introduce a new profile-based database of RNA motifs--RMfam--and illustrate some applications for investigating the evolution and functional characterization of RNA. All the data and scripts associated with this work are available from: https://github.com/ppgardne/RMfam.

The Rfam database (available at http://rfam.xfam.org) is a collection of non-coding RNA families represented by manually curated sequence alignments, consensus secondary structures and annotation gathered from corresponding Wikipedia, taxonomy and ontology resources. In this article, we detail updates and improvements to the Rfam data and website for the Rfam 12.0 release. We describe the upgrade of our search pipeline to use Infernal 1.1 and demonstrate its improved homology detection ability by comparison with the previous version. The new pipeline is easier for users to apply to their own data sets, and we illustrate its ability to annotate RNAs in genomic and metagenomic data sets of various sizes. Rfam has been expanded to include 260 new families, including the well-studied large subunit ribosomal RNA family, and for the first time includes information on short sequence- and structure-based RNA motifs present within families.

BackgroundComputational biology provides widely used and powerful software tools for testing and making inferences about biological data. In the face of rapidly increasing volumes of data, heuristic methods that trade software speed for accuracy may be employed. We are have studied these trade-offs using the results of a large number of independent software benchmarks, and evaluated whether external factors are indicative of accurate software.MethodWe have extracted accuracy and speed ranks from independent benchmarks of different bioinformatic software tools, and evaluated whether the speed, author reputation, journal impact, recency and developer efforts are indicative of accuracy.ResultsWe found that software speed, author reputation, journal impact, number of citations and age are all unreliable predictors of software accuracy. This is unfortunate because citations, author and journal reputation are frequently cited reasons for selecting software tools. However, GitHub-derived records and high version numbers show that the accurate bioinformatic software tools are generally the product of many improvements over time, often from multiple developers.DiscussionWe also find that the field of bioinformatics has a large excess of slow and inaccurate software tools, and this is consistent across many sub-disciplines. Meanwhile, there are few tools that are middle-of-road in terms of accuracy and speed trade-offs. We hypothesise that a form of publication-bias influences the publication and development of bioinformatic software. In other words, software that is intermediate in terms of both speed and accuracy may be difficult to publish - possibly due to author, editor and reviewer practices. This leaves an unfortunate hole in the literature as the ideal tools may fall into this gap. For example, high accuracy tools are not always useful if years of CPU time are required, while high speed is not useful if the results are also inaccurate.

Toxin-antitoxin (TA) systems are gene modules that appear to be widely horizontally mobile. It has been proposed that type I TA systems, with an antisense RNA-antitoxin, are less mobile than other TAs but no direct comparisons have been made. We searched for type I, II and III toxin families on chromosomes, plasmids and phages across bacterial phyla.The distribution of type I TA systems were more narrow than most type II and III system families, though this was less true of more recently discovered families.We discuss how the function and phenotypes of type I TA systems as well as biases in our databases and discovery of these modules may account for differences in their distribution.

Emerging high-throughput technologies have led to a deluge of putative non-coding RNA (ncRNA) sequences identified in a wide variety of organisms. Systematic characterization of these transcripts will be a tremendous challenge. Homology detection is critical to making maximal use of functional information gathered about ncRNAs: identifying homologous sequence allows us to transfer information gathered in one organism to another quickly and with a high degree of confidence. ncRNA presents a challenge for homology detection, as the primary sequence is often poorly conserved and de novo secondary structure prediction and search remains difficult. This protocol introduces methods developed by the Rfam database for identifying "families" of homologous ncRNAs starting from single "seed" sequences using manually curated sequence alignments to build powerful statistical models of sequence and structure conservation known as covariance models (CMs), implemented in the Infernal software package. We provide a step-by-step iterative protocol for identifying ncRNA homologs, then constructing an alignment and corresponding CM. We also work through an example for the bacterial small RNA MicA, discovering a previously unreported family of divergent MicA homologs in genus Xenorhabdus in the process., Submitted as a chapter for "Protocols in Bioinformatics". 41 pages, 7 figures

Metagenome studies are becoming increasingly widespread, yielding important insights into microbial communities covering diverse environments from terrestrial and aquatic ecosystems to human skin and gut. With the advent of high-throughput sequencing platforms, the use of large scale shotgun sequencing approaches is now commonplace. However, a thorough independent benchmark comparing state-of-the-art metagenome analysis tools is lacking. Here, we present a benchmark where the most widely used tools are tested on complex, realistic data sets. Our results clearly show that the most widely used tools are not necessarily the most accurate, that the most accurate tool is not necessarily the most time consuming and that there is a high degree of variability between available tools. These findings are important as the conclusions of any metagenomics study are affected by errors in the predicted community composition and functional capacity. Data sets and results are freely available from http://www.ucbioinformatics.org/metabenchmark.html

Automated RNA alignment algorithms often fail to recapture the essential conserved sites that are critical for function. To assist in the refinement of these algorithms, we manually curated a set of 148 alignments with a total of 9600 unique sequences, in which each alignment was backed by at least one crystal or NMR structure. These alignments included both naturally and artificially selected molecules. We used principles of isostericity to improve the alignments from an average of 83%–94% isosteric base pairs. We expect that this alignment collection will assist in a wide range of benchmarking efforts and provide new insight into evolutionary principles governing change in RNA structural motifs. The improved alignments have been contributed to the Rfam database.