Mutation Detection with Next-Generation Resequencing through a Mediator Genome

The affordability of next generation sequencing (NGS) is transforming the field of mutation analysis in bacteria. The genetic basis for phenotype alteration can be identified directly by sequencing the entire genome of the mutant and comparing it to the wild-type (WT) genome, thus identifying acquired mutations. A major limitation for this approach is the need for an a-priori sequenced reference genome for the WT organism, as the short reads of most current NGS approaches usually prohibit de-novo genome assembly. To overcome this limitation we propose a general framework that utilizes the genome of relative organisms as mediators for comparing WT and mutant bacteria. Under this framework, both mutant and WT genomes are sequenced with NGS, and the short sequencing reads are mapped to the mediator genome. Variations between the mutant and the mediator that recur in the WT are ignored, thus pinpointing the differences between the mutant and the WT. To validate this approach we sequenced the genome of Bdellovibrio bacteriovorus 109J, an obligatory bacterial predator, and its prey-independent mutant, and compared both to the mediator species Bdellovibrio bacteriovorus HD100. Although the mutant and the mediator sequences differed in more than 28,000 nucleotide positions, our approach enabled pinpointing the single causative mutation. Experimental validation in 53 additional mutants further established the implicated gene. Our approach extends the applicability of NGS-based mutant analyses beyond the domain of available reference genomes

Nucleotide and phylogenetic analyses of the Chlamydia trachomatis ompA gene indicates it is a hotspot for mutation

<p>Abstract</p> <p>Background</p> <p>Serovars of the human pathogen <it>Chlamydia trachomatis </it>occupy one of three specific tissue niches. Genomic analyses indicate that the serovars have a phylogeny congruent with their pathobiology and have an average substitution rate of less than one nucleotide per kilobase. In contrast, the gene that determines serovar specificity, <it>ompA</it>, has a phylogenetic association that is not congruent with tissue tropism and has a degree of nucleotide variability much higher than other genomic loci. The <it>ompA </it>gene encodes the major surface-exposed antigenic determinant, and the observed nucleotide diversity at the <it>ompA </it>locus is thought to be due to recombination and host immune selection pressure. The possible contribution of a localized increase in mutation rate, however, has not been investigated.</p> <p>Results</p> <p>Nucleotide diversity and phylogenetic relationships of the five constant and four variable domains of the <it>ompA </it>gene, as well as several loci surrounding <it>ompA</it>, were examined for each serovar. The loci flanking the <it>ompA </it>gene demonstrated that nucleotide diversity increased monotonically as <it>ompA </it>is approached and that their gene trees are not congruent with either <it>ompA </it>or tissue tropism. The variable domains of the <it>ompA </it>gene had a very high level of non-synonymous change, which is expected as these regions encode the surface-exposed epitopes and are under positive selection. However, the synonymous changes are clustered in the variable regions compared to the constant domains; if hitchhiking were to account for the increase in synonymous changes, these substitutions should be more evenly distributed across the gene. Recombination also cannot entirely account for this increase as the phylogenetic relationships of the constant and variable domains are congruent with each other.</p> <p>Conclusions</p> <p>The high number of synonymous substitutions observed within the variable domains of <it>ompA </it>appears to be due to an increased mutation rate within this region of the genome, whereas the increase in nucleotide substitution rate and the lack of phylogenetic congruence in the regions flanking <it>ompA </it>are characteristic motifs of gene conversion. Together, the increased mutation rate in the <it>ompA </it>gene, in conjunction with gene conversion and positive selection, results in a high degree of variability that promotes host immune evasion.</p

Analysis of pmpD Expression and PmpD Post-Translational Processing during the Life Cycle of Chlamydia trachomatis Serovars A, D, and L2

BACKGROUND: The polymorphic membrane protein D (PmpD) in Chlamydia is structurally similar to autotransporter proteins described in other bacteria and may be involved in cellular and humoral protective immunity against Chlamydia. The mechanism of PmpD post-translational processing and the role of its protein products in the pathogenesis of chlamydial infection have not been very well elucidated to date. METHODOLOGY/PRINCIPAL FINDINGS: Here we examined the expression and post-translational processing of the protein product of the pmpD gene during the life cycle of C. trachomatis serovars A, D, and L2. Each of these three serovars targets different human organs and tissues and encodes a different pmpD gene nucleotide sequence. Our quantitative real-time reverse transcription polymerase chain reaction results demonstrate that the pmpD gene is up-regulated at 12-24 hours after infection regardless of the Chlamydia serovar. This up-regulation is coincidental with the period of exponential growth and replication of reticulate bodies (RB) of Chlamydia and indicates a probable similarity in function of pmpD in serovars A, D, and L2 of Chlamydia. Using mass spectrometry analysis, we identified the protein products of post-translational processing of PmpD of C. trachomatis serovar L2 and propose a double pathway model for PmpD processing, with one cleavage site between the passenger and autotransporter domains and the other site in the middle of the passenger domain. Notably, when Chlamydia infected culture cells were subjected to low (28 degrees C) temperature, PmpD post-translational processing and secretion was found to be uninhibited in the resulting persistent infection. In addition, confocal microscopy of cells infected with Chlamydia confirms our earlier hypothesis that PmpD is secreted outside Chlamydia and its secretion increases with growth of the chlamydial inclusion. CONCLUSION/SIGNIFICANCE: The results of this current study involving multiple Chlamydia serovars support the general consensus that the pmpD gene is maximally expressed at mid infection and provide new information about PmpD as an autotransporter protein which is post-translationally processed and secreted outside Chlamydia during normal and low temperature induced persistent chlamydial infection

Rapid Detection of Chlamydia trachomatis and Typing of the Lymphogranuloma venereum associated L-Serovars by TaqMan PCR

<p>Abstract</p> <p>Background</p> <p>Infection due to <it>Chlamydia trachomatis </it>is the most common sexually transmitted bacterial disease of global health significance, and especially the L-serovars causing lymphogranuloma venereum are increasingly being found in Europe in men who have sex with men.</p> <p>Results</p> <p>The design and evaluation of a rapid, multiplex, real-time PCR targeting the major outer membrane protein (<it>omp-1</it>) -gene and a L-serovar-specific region of the polymorphic protein H (<it>pmp-H</it>) -gene for the detection of <it>Chlamydia trachomatis </it>is reported here. The PCR takes place as a single reaction with an internal control. For L1-, L2- and L3-serovar differentiation a second set of real-time PCRs was evaluated based on the amplification of serovar-specific <it>omp-1</it>-regions. The detection limit of each real-time PCR, multiplexed or not, was 50 genome copies per reaction with an efficiency ranging from 90,5–95,2%.</p> <p>In a retrospective analysis of 50 ocular, rectal and urogenital specimens formerly tested to be positive for <it>C. trachomatis </it>we identified six L2-serovars in rectal specimens of HIV-positive men, one in a double-infection with L3, and one L2 in a urethral specimen of an HIV-negative male.</p> <p>Conclusion</p> <p>This unique real-time PCR is specific and convenient for the rapid routine-diagnostic detection of lymphogranuloma venereum-associated L-serovars and enables the subsequent differentiation of L1, L2 and L3 for epidemiologic studies.</p

A colorimetric strategy based on dynamic chemistry for direct detection of Trypanosomatid species

Leishmaniasis and Chagas disease are endemic in many countries, and re-emerging in the developed countries. A rapid and accurate diagnosis is important for early treatment for reducing the duration of infection as well as for preventing further potential health complications. In this work, we have developed a novel colorimetric molecular assay that integrates nucleic acid analysis by dynamic chemistry (ChemNAT) with reverse dot-blot hybridization in an array format for a rapid and easy discrimination of Leishmania major and Trypanosoma cruzi. The assay consists of a singleplex PCR step that amplifies a highly homologous DNA sequence which encodes for the RNA component of the large ribosome subunit. The amplicons of the two different parasites differ between them by single nucleotide variations, known as “Single Nucleotide Fingerprint” (SNF) markers. The SNF markers can be easily identified by naked eye using a novel micro Spin-Tube device "Spin-Tube", as each of them creates a specific spot pattern. Moreover, the direct use of ribosomal RNA without requiring the PCR pre-amplification step is also feasible, further increasing the simplicity of the assay. The molecular assay delivers sensitivity capable of identifying up to 8.7 copies per μL with single mismatch specificity. The Spin-Tube thus represents an innovative solution providing benefits in terms of time, cost, and simplicity, all of which are crucial for the diagnosis of infectious disease in developing countries.This research work has received funding from Junta de Andalucía, Consejería de Economía e Innovación (project number 2012-BIO1778), the Spanish Ministerio de Economía y Competitividad (Grants CTQ2012-34778, BIO2016-80519-R, FPI Grant BES-2013- 063020). This research was partially supported by the 7th European Community Framework Program (FP7-PEOPLE-2012-CIG-Project Number 322276)

The mitochondrial genome sequence of the ciliate Paramecium caudatum reveals a shift in nucleotide composition and codon usage within the genus Paramecium

<p>Abstract</p> <p>Background</p> <p>Despite the fact that the organization of the ciliate mitochondrial genome is exceptional, only few ciliate mitochondrial genomes have been sequenced until today. All ciliate mitochondrial genomes are linear. They are 40 kb to 47 kb long and contain some 50 tightly packed genes without introns. Earlier studies documented that the mitochondrial guanine + cytosine contents are very different between <it>Paramecium tetraurelia </it>and all studied <it>Tetrahymena </it>species. This raises the question of whether the high mitochondrial G+C content observed in <it>P. tetraurelia </it>is a characteristic property of <it>Paramecium </it>mtDNA, or whether it is an exception of the ciliate mitochondrial genomes known so far. To test this question, we determined the mitochondrial genome sequence of <it>Paramecium caudatum </it>and compared the gene content and sequence properties to the closely related <it>P. tetraurelia</it>.</p> <p>Results</p> <p>The guanine + cytosine content of the <it>P. caudatum </it>mitochondrial genome was significantly lower than that of <it>P. tetraurelia </it>(22.4% vs. 41.2%). This difference in the mitochondrial nucleotide composition was accompanied by significantly different codon usage patterns in both species, i.e. within <it>P. caudatum </it>clearly A/T ending codons dominated, whereas for <it>P. tetraurelia </it>the synonymous codons were more balanced with a higher number of G/C ending codons. Further analyses indicated that the nucleotide composition of most members of the genus <it>Paramecium </it>resembles that of <it>P. caudatum </it>and that the shift observed in <it>P. tetraurelia </it>is restricted to the <it>P. aurelia </it>species complex.</p> <p>Conclusions</p> <p>Surprisingly, the codon usage bias in the <it>P. caudatum </it>mitochondrial genome, exemplified by the effective number of codons, is more similar to the distantly related <it>T. pyriformis </it>and other single-celled eukaryotes such as <it>Chlamydomonas</it>, than to the closely related <it>P. tetraurelia</it>. These differences in base composition and codon usage bias were, however, not reflected in the amino acid composition. Most probably, the observed picture is best explained by a hitherto unknown (neutral or adaptive) mechanism that increased the guanine + cytosine content in <it>P. tetraurelia </it>mtDNA on the one hand, and strong purifying selection on the ancestral amino acid composition on the other hand. These contradicting forces are counterbalanced by a considerably altered codon usage pattern.</p

Genome-level homology and phylogeny of Shewanella (Gammaproteobacteria: lteromonadales: Shewanellaceae)

<p>Abstract</p> <p>Background</p> <p>The explosion in availability of whole genome data provides the opportunity to build phylogenetic hypotheses based on these data as well as the ability to learn more about the genomes themselves. The biological history of genes and genomes can be investigated based on the taxomonic history provided by the phylogeny. A phylogenetic hypothesis based on complete genome data is presented for the genus <it>Shewanella </it>(Gammaproteobacteria: Alteromonadales: Shewanellaceae). Nineteen taxa from <it>Shewanella </it>(16 species and 3 additional strains of one species) as well as three outgroup species representing the genera <it>Aeromonas </it>(Gammaproteobacteria: Aeromonadales: Aeromonadaceae), <it>Alteromonas </it>(Gammaproteobacteria: Alteromonadales: Alteromonadaceae) and <it>Colwellia </it>(Gammaproteobacteria: Alteromonadales: Colwelliaceae) are included for a total of 22 taxa.</p> <p>Results</p> <p>Putatively homologous regions were found across unannotated genomes and tested with a phylogenetic analysis. Two genome-wide data-sets are considered, one including only those genomic regions for which all taxa are represented, which included 3,361,015 aligned nucleotide base-pairs (bp) and a second that additionally includes those regions present in only subsets of taxa, which totaled 12,456,624 aligned bp. Alignment columns in these large data-sets were then randomly sampled to create smaller data-sets. After the phylogenetic hypothesis was generated, genome annotations were projected onto the DNA sequence alignment to compare the historical hypothesis generated by the phylogeny with the functional hypothesis posited by annotation.</p> <p>Conclusions</p> <p>Individual phylogenetic analyses of the 243 locally co-linear genome regions all failed to recover the genome topology, but the smaller data-sets that were random samplings of the large concatenated alignments all produced the genome topology. It is shown that there is not a single orthologous copy of 16S rRNA across the taxon sampling included in this study and that the relationships among the multiple copies are consistent with 16S rRNA undergoing concerted evolution. Unannotated whole genome data can provide excellent raw material for generating hypotheses of historical homology, which can be tested with phylogenetic analysis and compared with hypotheses of gene function.</p

Plasmids and Rickettsial Evolution: Insight from Rickettsia felis

BACKGROUND: The genome sequence of Rickettsia felis revealed a number of rickettsial genetic anomalies that likely contribute not only to a large genome size relative to other rickettsiae, but also to phenotypic oddities that have confounded the categorization of R. felis as either typhus group (TG) or spotted fever group (SFG) rickettsiae. Most intriguing was the first report from rickettsiae of a conjugative plasmid (pRF) that contains 68 putative open reading frames, several of which are predicted to encode proteins with high similarity to conjugative machinery in other plasmid-containing bacteria. METHODOLOGY/PRINCIPAL FINDINGS: Using phylogeny estimation, we determined the mode of inheritance of pRF genes relative to conserved rickettsial chromosomal genes. Phylogenies of chromosomal genes were in agreement with other published rickettsial trees. However, phylogenies including pRF genes yielded different topologies and suggest a close relationship between pRF and ancestral group (AG) rickettsiae, including the recently completed genome of R. bellii str. RML369-C. This relatedness is further supported by the distribution of pRF genes across other rickettsiae, as 10 pRF genes (or inactive derivatives) also occur in AG (but not SFG) rickettsiae, with five of these genes characteristic of typical plasmids. Detailed characterization of pRF genes resulted in two novel findings: the identification of oriV and replication termination regions, and the likelihood that a second proposed plasmid, pRFδ, is an artifact of the original genome assembly. CONCLUSION/SIGNIFICANCE: Altogether, we propose a new rickettsial classification scheme with the addition of a fourth lineage, transitional group (TRG) rickettsiae, that is unique from TG and SFG rickettsiae and harbors genes from possible exchanges with AG rickettsiae via conjugation. We offer insight into the evolution of a plastic plasmid system in rickettsiae, including the role plasmids may have played in the acquirement of virulence traits in pathogenic strains, and the likely origin of plasmids within the rickettsial tree