Comparative genome analysis of Wolbachia strain wAu

BACKGROUND: Wolbachia intracellular bacteria can manipulate the reproduction of their arthropod hosts, including inducing sterility between populations known as cytoplasmic incompatibility (CI). Certain strains have been identified that are unable to induce or rescue CI, including wAu from Drosophila. Genome sequencing and comparison with CI-inducing related strain wMel was undertaken in order to better understand the molecular basis of the phenotype. RESULTS: Although the genomes were broadly similar, several rearrangements were identified, particularly in the prophage regions. Many orthologous genes contained single nucleotide polymorphisms (SNPs) between the two strains, but a subset containing major differences that would likely cause inactivation in wAu were identified, including the absence of the wMel ortholog of a gene recently identified as a CI candidate in a proteomic study. The comparative analyses also focused on a family of transcriptional regulator genes implicated in CI in previous work, and revealed numerous differences between the strains, including those that would have major effects on predicted function. CONCLUSIONS: The study provides support for existing candidates and novel genes that may be involved in CI, and provides a basis for further functional studies to examine the molecular basis of the phenotype

BamView: visualizing and interpretation of next-generation sequencing read alignments.

So-called next-generation sequencing (NGS) has provided the ability to sequence on a massive scale at low cost, enabling biologists to perform powerful experiments and gain insight into biological processes. BamView has been developed to visualize and analyse sequence reads from NGS platforms, which have been aligned to a reference sequence. It is a desktop application for browsing the aligned or mapped reads [Ruffalo, M, LaFramboise, T, Koyutürk, M. Comparative analysis of algorithms for next-generation sequencing read alignment. Bioinformatics 2011;27:2790-6] at different levels of magnification, from nucleotide level, where the base qualities can be seen, to genome or chromosome level where overall coverage is shown. To enable in-depth investigation of NGS data, various views are provided that can be configured to highlight interesting aspects of the data. Multiple read alignment files can be overlaid to compare results from different experiments, and filters can be applied to facilitate the interpretation of the aligned reads. As well as being a standalone application it can be used as an integrated part of the Artemis genome browser, BamView allows the user to study NGS data in the context of the sequence and annotation of the reference genome. Single nucleotide polymorphism (SNP) density and candidate SNP sites can be highlighted and investigated, and read-pair information can be used to discover large structural insertions and deletions. The application will also calculate simple analyses of the read mapping, including reporting the read counts and reads per kilobase per million mapped reads (RPKM) for genes selected by the user

Circlator: automated circularization of genome assemblies using long sequencing reads

The assembly of DNA sequence data is undergoing a renaissance thanks to emerging technologies capable of producing reads tens of kilobases long. Assembling complete bacterial and small eukaryotic genomes is now possible, but the final step of circularizing sequences remains unsolved. Here we present Circlator, the first tool to automate assembly circularization and produce accurate linear representations of circular sequences. Using Pacific Biosciences and Oxford Nanopore data, Circlator correctly circularized 26 of 27 circularizable sequences, comprising 11 chromosomes and 12 plasmids from bacteria, the apicoplast and mitochondrion of Plasmodium falciparum and a human mitochondrion. Circlator is available at http://sanger-pathogens.github.io/circlator/

Mucin glycosylation and sulphation in airway epithelial cells is not influenced by cystic fibrosis transmembrane conductance regulator expression

Abnormalities in mucus properties and clearance make a major contribution to the pathology of cystic fibrosis (CF). Our aim was to test the hypothesis that the defects in CF mucus are a direct result of mutations in the CF transmembrane conductance regulator (CFTR) protein. We evaluated a single mucin molecule MUC1F/5ACTR that carries tandem repeat sequence from MUC5AC, a major secreted airway mucin, in a MUC1 mucin vector. To establish whether the presence of mutant or normal CFTR directly influences the O-glycosylation and sulphation of mucins in airway epithelial cells, we used the CFT1-LC3 (DeltaF508 CFTR mutant) and CFT1-LCFSN (wild-type CFTR corrected) human airway epithelial cell lines. MUC1F/5ACTR mucin was immunoprecipitated, centricon purified, and O-glycosylation was evaluated by Matrix-assisted laser desorption ionization and electrospray tandem mass spectrometry to determine the composition of different carbohydrate structures. Mass spectrometry data showed the same O-glycans in both CFTR mutant and wild-type CFTR corrected cells. Metabolic labeling assays were performed to evaluate gross glycosylation and sulphation of the mucins and showed no significant difference in mucin synthesized in six independent clones of these cell lines. Our results show that the absence of functional CFTR protein causes neither an abnormality in mucin O-glycosylation nor an increase in mucin sulphation

PinR mediates the generation of reversible population diversity in Streptococcus zooepidemicus

Opportunistic pathogens must adapt to and survive in a wide range of complex ecosystems. Streptococcus zooepidemicus is an opportunistic pathogen of horses and many other animals, including humans. The assembly of different surface architecture phenotypes from one genotype is likely to be crucial to the successful exploitation of such an opportunistic lifestyle. Construction of a series of mutants revealed that a serine recombinase, PinR, inverts 114 bp of the promoter of SZO_08560, which is bordered by GTAGACTTTA and TAAAGTCTAC inverted repeats. Inversion acts as a switch, controlling the transcription of this sortase-processed protein, which may enhance the attachment of S. zooepidemicus to equine trachea. The genome of a recently sequenced strain of S. zooepidemicus, 2329 (Sz2329), was found to contain a disruptive internal inversion of 7 kb of the FimIV pilus locus, which is bordered by TAGAAA and TTTCTA inverted repeats. This strain lacks pinR and this inversion may have become irreversible following the loss of this recombinase. Active inversion of FimIV was detected in three strains of S. zooepidemicus, 1770 (Sz1770), B260863 (SzB260863) and H050840501 (SzH050840501), all of which encoded pinR. A deletion mutant of Sz1770 that lacked pinR was no longer capable of inverting its internal region of FimIV. The data highlight redundancy in the PinR sequence recognition motif around a short TAGA consensus and suggest that PinR can reversibly influence the wider surface architecture of S. zooepidemicus, providing this organism with a bet-hedging solution to survival in fluctuating environments

Artemis: an integrated platform for visualization and analysis of high-throughput sequence-based experimental data.

MOTIVATION: High-throughput sequencing (HTS) technologies have made low-cost sequencing of large numbers of samples commonplace. An explosion in the type, not just number, of sequencing experiments has also taken place including genome re-sequencing, population-scale variation detection, whole transcriptome sequencing and genome-wide analysis of protein-bound nucleic acids. RESULTS: We present Artemis as a tool for integrated visualization and computational analysis of different types of HTS datasets in the context of a reference genome and its corresponding annotation. AVAILABILITY: Artemis is freely available (under a GPL licence) for download (for MacOSX, UNIX and Windows) at the Wellcome Trust Sanger Institute websites: http://www.sanger.ac.uk/resources/software/artemis/

Spatiotemporal co-existence of two Mycobacterium ulcerans clonal complexes in the Offin River Valley of Ghana

In recent years, comparative genome sequence analysis of African Mycobacterium ulcerans strains isolated from Buruli ulcer (BU) lesion specimen has revealed a very limited genetic diversity of closely related isolates and a striking association between genotype and geographical origin of the patients. Here, we compared whole genome sequences of five M. ulcerans strains isolated in 2004 or 2013 from BU lesions of four residents of the Offin river valley with 48 strains isolated between 2002 and 2005 from BU lesions of individuals residing in the Densu river valley of Ghana. While all M. ulcerans isolates from the Densu river valley belonged to the same clonal complex, members of two distinct clonal complexes were found in the Offin river valley over space and time. The Offin strains were closely related to genotypes from either the Densu region or from the Asante Akim North district of Ghana. These results point towards an occasional involvement of a mobile reservoir in the transmission of M. ulcerans, enabling the spread of bacteria across different regions

Bayesian inference of ancestral dates on bacterial phylogenetic trees

The sequencing and comparative analysis of a collection of bacterial genomes from a single species or lineage of interest can lead to key insights into its evolution, ecology or epidemiology. The tool of choice for such a study is often to build a phylogenetic tree, and more specifically when possible a dated phylogeny, in which the dates of all common ancestors are estimated. Here, we propose a new Bayesian methodology to construct dated phylogenies which is specifically designed for bacterial genomics. Unlike previous Bayesian methods aimed at building dated phylogenies, we consider that the phylogenetic relationships between the genomes have been previously evaluated using a standard phylogenetic method, which makes our methodology much faster and scalable. This two-step approach also allows us to directly exploit existing phylogenetic methods that detect bacterial recombination, and therefore to account for the effect of recombination in the construction of a dated phylogeny. We analysed many simulated datasets in order to benchmark the performance of our approach in a wide range of situations. Furthermore, we present applications to three different real datasets from recent bacterial genomic studies. Our methodology is implemented in a R package called BactDating which is freely available for download at https://github.com/xavierdidelot/BactDating

Heterogeneity in the Frequency and Characteristics of Homologous Recombination in Pneumococcal Evolution

The bacterium Streptococcus pneumoniae (pneumococcus) is one of the most important human bacterial pathogens, and a leading cause of morbidity and mortality worldwide. The pneumococcus is also known for undergoing extensive homologous recombination via transformation with exogenous DNA. It has been shown that recombination has a major impact on the evolution of the pathogen, including acquisition of antibiotic resistance and serotype-switching. Nevertheless, the mechanism and the rates of recombination in an epidemiological context remain poorly understood. Here, we proposed several mathematical models to describe the rate and size of recombination in the evolutionary history of two very distinct pneumococcal lineages, PMEN1 and CC180. We found that, in both lineages, the process of homologous recombination was best described by a heterogeneous model of recombination with single, short, frequent replacements, which we call micro-recombinations, and rarer, multi-fragment, saltational replacements, which we call macro-recombinations. Macro-recombination was associated with major phenotypic changes, including serotype-switching events, and thus was a major driver of the diversification of the pathogen. We critically evaluate biological and epidemiological processes that could give rise to the micro-recombination and macro-recombination processes

Multi-State RNA Design with Geometric Multi-Graph Neural Networks

Computational RNA design has broad applications across synthetic biology and therapeutic development. Fundamental to the diverse biological functions of RNA is its conformational flexibility, enabling single sequences to adopt a variety of distinct 3D states. Currently, computational biomolecule design tasks are often posed as inverse problems, where sequences are designed based on adopting a single desired structural conformation. In this work, we propose gRNAde, a geometric RNA design pipeline that operates on sets of 3D RNA backbone structures to explicitly account for and reflect RNA conformational diversity in its designs. We demonstrate the utility of gRNAde for improving native sequence recovery over single-state approaches on a new large-scale 3D RNA design dataset, especially for multi-state and structurally diverse RNAs. Our code is available at https://github.com/chaitjo/geometric-rna-desig