Compressed Sparse FM-Index: Fast Sequence Alignment Using Large K-Steps

The FM-index is a data structure used in genomics for exact search of input sequences over large reference genomes. Algorithms based on the FM-index show an irregular memory access pattern, resulting in a memory bound problem. We analyze a recent implementation of the FM-index and highlight existing throughput-memory trade-offs, showing that memory requirements limit implementation of large k-steps. We propose COFI, a COmpressed FM-Index for large K-steps. COFI enables a 15-step FM-index using less than 16 GB for a human genome reference of 3 giga base pairs. An algorithm based on this new layout is evaluated on both a Knights Landing (KNL) and an Skylake-based system (SKX). We achieve average speed-ups of 1.46× and 1.39×, respectively, with respect to an state-of-the-art FM-index implementation that is already well optimized.

Mutant Strains of Escherichia coli and Methicillin-Resistant Staphylococcus aureus Obtained by Laboratory Selection To Survive on Metallic Copper Surfaces

Artificial laboratory evolution was used to produce mutant strains of Escherichia coli and methicillin-resistant Staphylococcus aureus (MRSA) able to survive on antimicrobial metallic copper surfaces. These mutants were 12- and 60-fold less susceptible to the copper-mediated contact killing process than their respective parent strains. Growth levels of the mutant and its parent in complex growth medium were similar. Tolerance to copper ions of the mutants was unchanged. The mutant phenotype remained stable over about 250 generations under nonstress conditions. The mutants and their respective parental strains accumulated copper released from the metallic surfaces to similar extents. Nevertheless, only the parental strains succumbed to copper stress when challenged on metallic copper surfaces, suffering complete destruction of the cell structure. Whole-genome sequencing and global transcriptome analysis were used to decipher the genetic alterations in the mutant strains; however, these results did not explain the copper-tolerance phenotypes on the systemic level. Instead, the mutants shared features with those of stressed bacterial subpopulations entering the early or "shallow" persister state. In contrast to the canonical persister state, however, the ability to survive on solid copper surfaces was adopted by the majority of the mutant strain population. This indicated that application of solid copper surfaces in hospitals and elsewhere has to be accompanied by strict cleaning regimens to keep the copper surfaces active and prevent evolution of tolerant mutant strains.IMPORTANCE Microbes are rapidly killed on solid copper surfaces by contact killing. Copper surfaces thus have an important role to play in preventing the spread of nosocomial infections. Bacteria adapt to challenging natural and clinical environments through evolutionary processes, for instance, by acquisition of beneficial spontaneous mutations. We wish to address the question of whether mutants can be selected that have evolved to survive contact killing on solid copper surfaces. We isolated such mutants from Escherichia coli and methicillin-resistant Staphylococcus aureus (MRSA) by artificial laboratory evolution. The ability to survive on solid copper surfaces was a stable phenotype of the mutant population and not restricted to a small subpopulation. As a consequence, standard operation procedures with strict hygienic measures are extremely important to prevent the emergence and spread of copper-surface-tolerant persister-like bacterial strains if copper surfaces are to be sustainably used to limit the spread of pathogenic bacteria, e.g., to curb nosocomial infections.

Accelerating Sequence Alignments Based on FM-Index Using the Intel KNL Processor

FM-index is a compact data structure suitable for fast matches of short reads to large reference genomes. The matching algorithm using this index exhibits irregular memory access patterns that cause frequent cache misses, resulting in a memory bound problem. This paper analyzes different FM-index versions presented in the literature, focusing on those computing aspects related to the data access. As a result of the analysis, we propose a new organization of FM-index that minimizes the demand for memory bandwidth, allowing a great improvement of performance on processors with high-bandwidth memory, such as the second-generation Intel Xeon Phi (Knights Landing, or KNL), integrating ultra high-bandwidth stacked memory technology. As the roofline model shows, our implementation reaches 95 percent of the peak random access bandwidth limit when executed on the KNL and almost all of the available bandwidth when executed on other Intel Xeon architectures with conventional DDR memory. In addition, the obtained throughput in KNL is much higher than the results reported for GPUs in the literature.

SRPRISM (Single Read Paired Read Indel Substitution Minimizer): an efficient aligner for assemblies with explicit guarantees

BACKGROUND

Alignment of sequence reads generated by next-generation sequencing is an integral part of most pipelines analyzing next-generation sequencing data. A number of tools designed to quickly align a large volume of sequences are already available. However, most existing tools lack explicit guarantees about their output. They also do not support searching genome assemblies, such as the human genome assembly GRCh38, that include primary and alternate sequences and placement information for alternate sequences to primary sequences in the assembly.

FINDINGS

This paper describes SRPRISM (Single Read Paired Read Indel Substitution Minimizer), an alignment tool for aligning reads without splices. SRPRISM has features not available in most tools, such as (i) support for searching genome assemblies with alternate sequences, (ii) partial alignment of reads with a specified region of reads to be included in the alignment, (iii) choice of ranking schemes for alignments, and (iv) explicit criteria for search sensitivity. We compare the performance of SRPRISM to GEM, Kart, STAR, BWA-MEM, Bowtie2, Hobbes, and Yara using benchmark sets for paired and single reads of lengths 100 and 250 bp generated using DWGSIM. SRPRISM found the best results for most benchmark sets with error rate of up to ∼2.5% and GEM performed best for higher error rates. SRPRISM was also more sensitive than other tools even when sensitivity was reduced to improve run time performance.

CONCLUSIONS

We present SRPRISM as a flexible read mapping tool that provides explicit guarantees on results.

Drug-resistance profiling and transmission dynamics of multidrug-resistant Mycobacterium tuberculosis in Saudi Arabia revealed by whole genome sequencing

Background

In Saudi Arabia, cross-border transmission of multidrug-resistant (MDR) Mycobacterium tuberculosis complex (MTBC) strains might be particularly fostered by high immigration rates. Herein, we aimed to elucidate the transmission dynamics of MDR-MTBC strains and reveal a detailed prediction of all resistance-conferring mutations for the first- and second-line drugs.

Methods

We investigated all MDR-MTBC strains collected between 2015 and 2017 from provincial mycobacteria referral laboratories and compared demographic and clinical parameters to a cohort of non-MDR-TB patients using a whole genome sequencing approach. Clusters were defined based on a maximum strain-to-strain genetic distance of five single-nucleotide polymorphisms (SNPs) as surrogate marker for recent transmission, and then investigated molecular drug-resistance markers (37 genes).

Results

Forty-eight (67.6%) MDR-MTBC strains were grouped in 14 different clusters, ranging in size from two to six strains; 22.5% (16/71) of all MDR-MTBC isolates were predicted to be fully resistant to all five first-line drugs, and five strains (7.0%) exhibited fluoroquinolone resistance. Moreover, we revealed the presence of 12 compensatory mutations as well as 26 non-synonymous SNPs in the rpoC gene and non-hotspot region in rpoB, respectively.

Conclusion

Optimized TB molecular surveillance, diagnosis, and patient management are urgently needed to contain MDR-MTBC transmission and prevent the development of additional drug resistance.

Redox-Sensing Under Hypochlorite Stress and Infection Conditions by the Rrf2-Family Repressor HypR in Staphylococcus aureus

AIMS

Staphylococcus aureus is a major human pathogen and has to cope with reactive oxygen and chlorine species (ROS, RCS) during infections, which requires efficient protection mechanisms to avoid destruction. Here, we have investigated the changes in the RNA-seq transcriptome by the strong oxidant sodium hypochlorite (NaOCl) in S. aureus USA300 to identify novel redox-sensing mechanisms that provide protection under infection conditions.

RESULTS

NaOCl stress caused an oxidative stress response in S. aureus as indicated by the induction of the PerR, QsrR, HrcA, and SigmaB regulons in the RNA-seq transcriptome. The hypR-merA (USA300HOU_0588-87) operon was most strongly upregulated under NaOCl stress, which encodes for the Rrf2-family regulator HypR and the pyridine nucleotide disulfide reductase MerA. We have characterized HypR as a novel redox-sensitive repressor that controls MerA expression and directly senses and responds to NaOCl and diamide stress via a thiol-based mechanism in S. aureus. Mutational analysis identified Cys33 and the conserved Cys99 as essential for NaOCl sensing, while Cys99 is also important for repressor activity of HypR in vivo. The redox-sensing mechanism of HypR involves Cys33-Cys99 intersubunit disulfide formation by NaOCl stress both in vitro and in vivo. Moreover, the HypR-controlled flavin disulfide reductase MerA was shown to protect S. aureus against NaOCl stress and increased survival in J774A.1 macrophage infection assays. Conclusion and Innovation: Here, we identified a new member of the widespread Rrf2 family as redox sensor of NaOCl stress in S. aureus that uses a thiol/disulfide switch to regulate defense mechanisms against the oxidative burst under infections in S. aureus. Antioxid. Redox Signal. 29, 615-636.

Harmonized Genome Wide Typing of Tubercle Bacilli Using a Web-Based Gene-By-Gene Nomenclature System

Background

Global tuberculosis (TB) control is challenged by uncontrolled transmission of Mycobacterium tuberculosis complex (Mtbc) strains, esp. of multidrug (MDR) or extensively resistant (XDR) variants. Precise analysis of transmission networks is the basis to trace outbreak M/XDR clones and improve TB control. However, classical genotyping tools lack discriminatory power due to the high similarity of strains of particular successful lineages, e.g. Beijing or outbreak strains. This can be overcome by whole genome sequencing (WGS) approaches, but these are not yet standardized to facilitate larger investigations encompassing different laboratories or outbreak tracing across borders.

Methods

We established and improved a whole genome gene-by-gene multi locus sequence typing approach encompassing a stable set of core genome genes (cgMLST) and linked it to a web-based nomenclature server (cgMLST.org) facilitating assignment and storage of allele numbers.

Findings

We evaluated and refined a previously suggested cgMLST schema by using a reference strain set (n = 251) reflecting the global diversity of the Mtbc. A set of 2891 genes showed excellent performance with at least 97% of the genes reliably identified in strains of all Mtbc lineages and in discriminating outbreak strains. cgMLST allele numbers were automatically retrieved from and stored at cgMLST.org.

Interpretation

The refined cgMLST schema provides high resolution genome-based typing of clinical strains of all Mtbc lineages. Combined with a web-based nomenclature server, it facilitates rapid, high-resolution, and harmonized tracing of clinical Mtbc strains needed for prospective local and global surveillance.

GRIM-Filter: Fast seed location filtering in DNA read mapping using processing-in-memory technologies

Background

Seed location filtering is critical in DNA read mapping, a process where billions of DNA fragments (reads) sampled from a donor are mapped onto a reference genome to identify genomic variants of the donor. State-of-the-art read mappers 1) quickly generate possible mapping locations for seeds (i.e., smaller segments) within each read, 2) extract reference sequences at each of the mapping locations, and 3) check similarity between each read and its associated reference sequences with a computationally-expensive algorithm (i.e., sequence alignment) to determine the origin of the read. A seed location filter comes into play before alignment, discarding seed locations that alignment would deem a poor match. The ideal seed location filter would discard all poor match locations prior to alignment such that there is no wasted computation on unnecessary alignments.

Results

We propose a novel seed location filtering algorithm, GRIM-Filter, optimized to exploit 3D-stacked memory systems that integrate computation within a logic layer stacked under memory layers, to perform processing-in-memory (PIM). GRIM-Filter quickly filters seed locations by 1) introducing a new representation of coarse-grained segments of the reference genome, and 2) using massively-parallel in-memory operations to identify read presence within each coarse-grained segment. Our evaluations show that for a sequence alignment error tolerance of 0.05, GRIM-Filter 1) reduces the false negative rate of filtering by 5.59x–6.41x, and 2) provides an end-to-end read mapper speedup of 1.81x–3.65x, compared to a state-of-the-art read mapper employing the best previous seed location filtering algorithm.

Conclusion

GRIM-Filter exploits 3D-stacked memory, which enables the efficient use of processing-in-memory, to overcome the memory bandwidth bottleneck in seed location filtering. We show that GRIM-Filter significantly improves the performance of a state-of-the-art read mapper. GRIM-Filter is a universal seed location filter that can be applied to any read mapper. We hope that our results provide inspiration for new works to design other bioinformatics algorithms that take advantage of emerging technologies and new processing paradigms, such as processing-in-memory using 3D-stacked memory devices.

New Mycobacterium tuberculosis Complex Sublineage, Brazzaville, Congo

Tuberculosis is a leading cause of illness and death in Congo. No data are available about the population structure and transmission dynamics of the Mycobacterium tuberculosis complex strains prevalent in this central Africa country. On the basis of single-nucleotide polymorphisms detected by whole-genome sequencing, we phylogenetically characterized 74 MTBC isolates from Brazzaville, the capital of Congo. The diversity of the study population was high; most strains belonged to the Euro-American lineage, which split into Latin American Mediterranean, Uganda I, Uganda II, Haarlem, X type, and a new dominant sublineage named Congo type (n = 26). Thirty strains were grouped in 5 clusters (each within 12 single-nucleotide polymorphisms), from which 23 belonged to the Congo type. High cluster rates and low genomic diversity indicate recent emergence and transmission of the Congo type, a new Euro-American sublineage of MTBC.

Graphics processing units in bioinformatics, computational biology and systems biology

Several studies in Bioinformatics, Computational Biology and Systems Biology rely on the definition of physico-chemical or mathematical models of biological systems at different scales and levels of complexity, ranging from the interaction of atoms in single molecules up to genome-wide interaction networks. Traditional computational methods and software tools developed in these research fields share a common trait: they can be computationally demanding on Central Processing Units (CPUs), therefore limiting their applicability in many circumstances. To overcome this issue, general-purpose Graphics Processing Units (GPUs) are gaining an increasing attention by the scientific community, as they can considerably reduce the running time required by standard CPU-based software, and allow more intensive investigations of biological systems. In this review, we present a collection of GPU tools recently developed to perform computational analyses in life science disciplines, emphasizing the advantages and the drawbacks in the use of these parallel architectures. The complete list of GPU-powered tools here reviewed is available at http://bit.ly/gputools.

Physiological roles of sigma factor SigD in Corynebacterium glutamicum

BACKGROUND

Sigma factors are one of the components of RNA polymerase holoenzymes, and an essential factor of transcription initiation in bacteria. Corynebacterium glutamicum possesses seven genes coding for sigma factors, most of which have been studied to some detail; however, the role of SigD in transcriptional regulation in C. glutamicum has been mostly unknown.

RESULTS

In this work, pleiotropic effects of sigD overexpression at the level of phenotype, transcripts, proteins and metabolites were investigated. Overexpression of sigD decreased the growth rate of C. glutamicum cultures, and induced several physiological effects such as reduced culture foaming, turbid supernatant and cell aggregation. Upon overexpression of sigD, the level of Cmt1 (corynomycolyl transferase) in the supernatant was notably enhanced, and carbohydrate-containing compounds were excreted to the supernatant. The real-time PCR analysis revealed that sigD overexpression increased the expression of genes related to corynomycolic acid synthesis (fadD2, pks), genes encoding corynomycolyl transferases (cop1, cmt1, cmt2, cmt3), L, D-transpeptidase (lppS), a subunit of the major cell wall channel (porH), and the envelope lipid regulation factor (elrF). Furthermore, overexpression of sigD resulted in trehalose dicorynomycolate accumulation in the cell envelope.

CONCLUSIONS

This study demonstrated that SigD regulates the synthesis of corynomycolate and related compounds, and expanded the knowledge of regulatory functions of sigma factors in C. glutamicum.

Global outbreak of severe Mycobacterium chimaera disease after cardiac surgery: a molecular epidemiological study

BACKGROUND

Since 2013, over 100 cases of Mycobacterium chimaera prosthetic valve endocarditis and disseminated disease were notified in Europe and the USA, linked to contaminated heater-cooler units (HCUs) used during cardiac surgery. We did a molecular epidemiological investigation to establish the source of these patients' disease.

METHODS

We included 24 M chimaera isolates from 21 cardiac surgery-related patients in Switzerland, Germany, the Netherlands, and the UK, 218 M chimaera isolates from various types of HCUs in hospitals, from LivaNova (formerly Sorin; London, UK) and Maquet (Rastatt, Germany) brand HCU production sites, and unrelated environmental sources and patients, as well as eight Mycobacterium intracellulare isolates. Isolates were analysed by next-generation whole-genome sequencing using Illumina and Pacific Biosciences technologies, and compared with published M chimaera genomes.

FINDINGS

Phylogenetic analysis based on whole-genome sequencing of 250 isolates revealed two major M chimaera groups. Cardiac surgery-related patient isolates were all classified into group 1, in which all, except one, formed a distinct subgroup. This subgroup also comprised isolates from 11 cardiac surgery-related patients reported from the USA, most isolates from LivaNova HCUs, and one from their production site. Isolates from other HCUs and unrelated patients were more widely distributed in the phylogenetic tree.

INTERPRETATION

HCU contamination with M chimaera at the LivaNova factory seems a likely source for cardiothoracic surgery-related severe M chimaera infections diagnosed in Switzerland, Germany, the Netherlands, the UK, the USA, and Australia. Protective measures and heightened clinician awareness are essential to guarantee patient safety.

FUNDING

Partly funded by the EU Horizon 2020 programme, its FP7 programme, the German Center for Infection Research (DZIF), the Swiss National Science Foundation, the Swiss Federal Office of Public Health, and National Institute of Health Research Oxford Health Protection Research Units on Healthcare Associated Infection and Antimicrobial Resistance.

A joint cross-border investigation of a cluster of multidrug-resistant tuberculosis in Austria, Romania and Germany in 2014 using classic, genotyping and whole genome sequencing methods: lessons learnt

Molecular surveillance of multidrug-resistant tuberculosis (MDR-TB) using 24-loci MIRU-VNTR in the European Union suggests the occurrence of international transmission. In early 2014, Austria detected a molecular MDR-TB cluster of five isolates. Links to Romania and Germany prompted the three countries to investigate possible cross-border MDR-TB transmission jointly. We searched genotyping databases, genotyped additional isolates from Romania, used whole genome sequencing (WGS) to infer putative transmission links, and investigated pairwise epidemiological links and patient mobility. Ten isolates from 10 patients shared the same 24-loci MIRU-VNTR pattern. Within this cluster, WGS defined two subgroups of four patients each. The first comprised an MDR-TB patient from Romania who had sought medical care in Austria and two patients from Austria. The second comprised patients, two of them epidemiologically linked, who lived in three different countries but had the same city of provenance in Romania. Our findings strongly suggested that the two cases in Austrian citizens resulted from a newly introduced MDR-TB strain, followed by domestic transmission. For the other cases, transmission probably occurred in the same city of provenance. To prevent further MDR-TB transmission, we need to ensure universal access to early and adequate therapy and collaborate closely in tuberculosis care beyond administrative borders.

Monitoring global protein thiol-oxidation and protein S-mycothiolation in Mycobacterium smegmatis under hypochlorite stress

Mycothiol (MSH) is the major low molecular weight (LMW) thiol in Actinomycetes. Here, we used shotgun proteomics, OxICAT and RNA-seq transcriptomics to analyse protein S-mycothiolation, reversible thiol-oxidations and their impact on gene expression in Mycobacterium smegmatis under hypochlorite stress. In total, 58 S-mycothiolated proteins were identified under NaOCl stress that are involved in energy metabolism, fatty acid and mycolic acid biosynthesis, protein translation, redox regulation and detoxification. Protein S-mycothiolation was accompanied by MSH depletion in the thiol-metabolome. Quantification of the redox state of 1098 Cys residues using OxICAT revealed that 381 Cys residues (33.6%) showed >10% increased oxidations under NaOCl stress, which overlapped with 40 S-mycothiolated Cys-peptides. The absence of MSH resulted in a higher basal oxidation level of 338 Cys residues (41.1%). The RseA and RshA anti-sigma factors and the Zur and NrdR repressors were identified as NaOCl-sensitive proteins and their oxidation resulted in an up-regulation of the SigH, SigE, Zur and NrdR regulons in the RNA-seq transcriptome. In conclusion, we show here that NaOCl stress causes widespread thiol-oxidation including protein S-mycothiolation resulting in induction of antioxidant defense mechanisms in M. smegmatis. Our results further reveal that MSH is important to maintain the reduced state of protein thiols.

Genome-wide determination of transcription start sites reveals new insights into promoter structures in the actinomycete Corynebacterium glutamicum

The genome-wide identification of transcription start sites, enabled by high-throughput sequencing of a cDNA library enriched for native 5' transcript ends, is ideally suited for the analysis of promoters. Here, the transcriptome of Corynebacterium glutamicum, a non-pathogenic soil bacterium from the actinomycetes branch that is used in industry for the production of amino acids, was analysed by transcriptome sequencing of the 5'-ends of native transcripts. Total RNA samples were harvested from the exponential phase of growth, therefore the study mainly addressed promoters recognized by the main house-keeping sigma factor σ^A. The identification of 2454 transcription start sites (TSS) allowed the detailed analysis of most promoters recognized by σ^A and furthermore enabled us to form different promoter groups according to their location relative to protein-coding regions. These groups included leaderless transcripts (546 promoters), short-leadered (<500 bases) transcripts (917), and long-leadered (>500 bases) transcripts (173) as well as intragenic (557) and antisense transcripts (261). All promoters and the individual groups were searched for information, e.g. conserved residues and promoter motifs, and general design features as well as group-specific preferences were identified. A purine was found highly favored as TSS, whereas the -1 position was dominated by pyrimidines. The spacer between TSS and -10 region were consistently 6-7 bases and the -10 promoter motif was generally visible, whereas a recognizable -35 region was only occurring in a smaller fraction of promoters (7.5%) and enriched for leadered and antisense transcripts but depleted for leaderless transcripts. Promoters showing an extended -10 region were especially frequent in case of non-canonical -10 motifs (45.5%). Two bases downstream of the -10 core region, a G was conserved, exceeding 40% abundance in most groups. This fraction reached 74.6% for a group of putative σ^B-dependent promoters, thus giving a hint to a specific property of these promoters. In addition, the high number of promoters analysed allowed finding of subtle signals only showing up significantly with this large set. This included the observation of a periodically changing A+T-content with maxima spaced by a full turn of the DNA helix. This periodic structure includes the A+T-rich UP-element of bacterial promoters known before but was found to extend up to -100, indicating hitherto unknown constraints influencing promoter architecture and possibly also promoter function.

Short Read Mapping: An Algorithmic Tour

Ultra-high-throughput next-generation sequencing (NGS) technology allows us to determine the sequence of nucleotides of many millions of DNA molecules in parallel. Accompanied by a dramatic reduction in cost since its introduction in 2004, NGS technology has provided a new way of addressing a wide range of biological and biomedical questions, from the study of human genetic disease to the analysis of gene expression, protein-DNA interactions, and patterns of DNA methylation. The data generated by NGS instruments comprise huge numbers of very short DNA sequences, or 'reads', that carry little information by themselves. These reads therefore have to be pieced together by well-engineered algorithms to reconstruct biologically meaningful measurments, such as the level of expression of a gene. To solve this complex, high-dimensional puzzle, reads must be mapped back to a reference genome to determine their origin Due to sequencing errors and to genuine differences between the reference genome and the individual being sequenced, this mapping process must be tolerant of mismatches, insertions, and deletions. Although optimal alignment algorithms to solve this problem have long been available, the practical requirements of aligning hundreds of millions of short reads to the 3 billion base pair long human genome have stimulated the development of new, more efficient methods, which today are used routinely throughout the world for the analysis of NGS data.

Three-base periodicity of sites of sequence variation in Pseudomonas aeruginosa and Staphylococcus aureus core genomes

The three-base periodicity property is characteristic of protein-coding sequences. Here, we report on three-base periodicity of sequence variation in the core genome of bacteria. Single nucleotide polymorphism (SNP) syntenies were extracted from pairwise genome alignments of 41 Staphylococcus aureus or 20 Pseudomonas aeruginosa strains. The length of fragment pairs with identical nucleotides at all SNP positions showed a length-dependent overrepresentation of multiples of three nucleotides at corresponding codon positions of the AT-rich S. aureus and the GC-rich P. aeruginosa. Three-base SNP periodicity seems to be a characteristic feature of the tightly arranged bacterial core genome.

SNP synteny analysis of Staphylococcus aureus and Pseudomonas aeruginosa population genomics

Genomic sequence diversity of a bacterial species mainly results from the frequency distribution of single nucleotide polymorphisms (SNPs). Here we report on an SNP matrix-based binary algorithm to determine the intraclonal or interclonal genomic diversity by the number of shared sequential SNPs, the so-called SNP synteny or haplotype. All SNP positions and the frequency and length distribution of haplotypes are determined from pairwise alignment of completely sequenced genomes. This metric is invariant regarding the reference genome chosen. Information is obtained about the size of haplotypes, genomic gradients of recombination frequency, relatedness of strains and population composition of a taxon or clonal populations. The approach is illustrated with whole genome data sets of Staphylococcus aureus and Pseudomonas aeruginosa strains.

Tracing Mycobacterium tuberculosis transmission by whole genome sequencing in a high incidence setting: a retrospective population-based study in East Greenland

In East Greenland, a dramatic increase of tuberculosis (TB) incidence has been observed in recent years. Classical genotyping suggests a genetically similar Mycobacterium tuberculosis (Mtb) strain population as cause, however, precise transmission patterns are unclear. We performed whole genome sequencing (WGS) of Mtb isolates from 98% of culture-positive TB cases through 21 years (n = 182) which revealed four genomic clusters of the Euro-American lineage (mainly sub-lineage 4.8 (n = 134)). The time to the most recent common ancestor of lineage 4.8 strains was found to be 100 years. This sub-lineage further diversified in the 1970s, and massively expanded in the 1990s, a period of lowered TB awareness in Greenland. Despite the low genetic strain diversity, WGS data revealed several recent short-term transmission events in line with the increasing incidence in the region. Thus, the isolated setting and the uniformity of circulating Mtb strains indicated that the majority of East Greenlandic TB cases originated from one or few strains introduced within the last century. Thereby, the study shows the consequences of even short interruptions in TB control efforts in previously TB high incidence areas and demonstrates the potential role of WGS in detecting ongoing micro epidemics, thus guiding public health efforts in the future.

Deciphering the Transcriptional Response Mediated by the Redox-Sensing System HbpS-SenS-SenR from Streptomycetes

The secreted protein HbpS, the membrane-embedded sensor kinase SenS and the cytoplasmic response regulator SenR from streptomycetes have been shown to form a novel type of signaling pathway. Based on structural biology as well as different biochemical and biophysical approaches, redox stress-based post-translational modifications in the three proteins were shown to modulate the activity of this signaling pathway. In this study, we show that the homologous system, named here HbpSc-SenSc-SenRc, from the model species Streptomyces coelicolor A3(2) provides this bacterium with an efficient defense mechanism under conditions of oxidative stress. Comparative analyses of the transcriptomes of the Streptomyces coelicolor A3(2) wild-type and the generated hbpSc-senSc-senRc mutant under native and oxidative-stressing conditions allowed to identify differentially expressed genes, whose products may enhance the anti-oxidative defense of the bacterium. Amongst others, the results show an up-regulated transcription of genes for biosynthesis of cysteine and vitamin B₁₂, transport of methionine and vitamin B₁₂, and DNA synthesis and repair. Simultaneously, transcription of genes for degradation of an anti-oxidant compound is down-regulated in a HbpSc-SenSc-SenRc-dependent manner. It appears that HbpSc-SenSc-SenRc controls the non-enzymatic response of Streptomyces coelicolor A3(2) to counteract the hazardous effects of oxidative stress. Binding of the response regulator SenRc to regulatory regions of some of the studied genes indicates that the regulation is direct. The results additionally suggest that HbpSc-SenSc-SenRc may act in concert with other regulatory modules such as a transcriptional regulator, a two-component system and the Streptomyces B₁₂ riboswitch. The transcriptomics data, together with our previous in vitro results, enable a profound characterization of the HbpS-SenS-SenR system from streptomycetes. Since homologues to HbpS-SenS-SenR are widespread in different actinobacteria with ecological and medical relevance, the data presented here will serve as a basis to elucidate the biological role of these homologues.

Emended description of Mycobacterium abscessus, Mycobacterium abscessus subsp. abscessus and Mycobacteriumabscessus subsp. bolletii and designation of Mycobacteriumabscessus subsp. massiliense comb. nov

The taxonomic position of members of the Mycobacterium abscessus complex has been the subject of intensive investigation and, in some aspects confusion, in recent years as a result of varying approaches to genetic data interpretation. Currently, the former species Mycobacterium massiliense and Mycobacterium bolletii are grouped together as Mycobacterium abscessus subsp. bolletii. They differ greatly, however, as the former M. bolletii has a functional erm(41) gene that confers inducible resistance to macrolides, the primary therapeutic antimicrobials for M. abscessus, while in the former M. massiliense the erm(41) gene is non-functional. Furthermore, previous whole genome studies of the M. abscessus group support the separation of M. bolletii and M. massiliense. To shed further light on the population structure of Mycobacterium abscessus, 43 strains and three genomes retrieved from GenBank were subjected to pairwise comparisons using three computational approaches: verage ucleotide dentity, enome to enome istance and single nucleotide polymorphism analysis. The three methods produced overlapping results, each demonstrating three clusters of strains corresponding to the same number of taxonomic entities. The distances were insufficient to warrant distinction at the species level, but met the criteria for differentiation at the subspecies level. Based on prior erm(41)-related phenotypic data and current genomic data, we conclude that the species M. abscessus encompasses, in adjunct to the presently recognized subspecies M. abscessus subsp. abscessus and M. abscessus subsp. bolletii, a third subspecies for which we suggest the name M. abscessus subsp. massiliense comb. nov. (type strain CCUG 48898T=CIP 108297T=DSM 45103T=KCTC 19086T).

Comparative transcriptome analysis of the biocontrol strain Bacillus amyloliquefaciens FZB42 as response to biofilm formation analyzed by RNA sequencing

The strain Bacillus amyloliquefaciens FZB42 is a plant growth promoting rhizobacterium (PGPR) and biocontrol agent known to keep infections of lettuce (Lactuca sativa) by the phytopathogen Rhizoctonia solani down. Several mechanisms, including the production of secondary metabolites possessing antimicrobial properties and induction of the host plant's systemic resistance (ISR), were proposed to explain the biocontrol effect of the strain. B. amyloliquefaciens FZB42 is able to form plaques (biofilm-like structures) on plant roots and this feature was discussed to be associated with its biocontrol properties. For this reason, formation of B. amyloliquefaciens biofilms was studied at the transcriptional level using high-throughput sequencing of whole transcriptome cDNA libraries from cells grown under biofilm-forming conditions vs. planktonic growth. Comparison of the transcriptional profiles of B. amyloliquefaciens FZB42 under these growth conditions revealed a common set of highly transcribed genes mostly associated with basic cellular functions. The lci gene, encoding an antimicrobial peptide (AMP), was among the most highly transcribed genes of cells under both growth conditions suggesting that AMP production may contribute to biocontrol. In contrast, gene clusters coding for synthesis of secondary metabolites with antimicrobial properties were only moderately transcribed and not induced in biofilm-forming cells. Differential gene expression revealed that 331 genes were significantly up-regulated and 230 genes were down-regulated in the transcriptome of B. amyloliquefaciens FZB42 under biofilm-forming conditions in comparison to planktonic cells. Among the most highly up-regulated genes, the yvqHI operon, coding for products involved in nisin (class I bacteriocin) resistance, was identified. In addition, an operon whose products play a role in fructosamine metabolism was enhanced in its transcription. Moreover, genes involved in the production of the extracellular biofilm matrix including exopolysaccharide genes (eps) and the yqxM-tasA-sipW operon encoding amyloid fiber synthesis were up-regulated in the B. amyloliquefaciens FZB42 biofilm. On the other hand, highly down-regulated genes in biofilms are associated with synthesis, assembly and regulation of the flagellar apparatus, the degradation of aromatic compounds and the export of copper. The obtained transcriptional profile for B. amyloliquefaciens biofilm cells uncovered genes involved in its development and enabled the assessment that synthesis of secondary metabolites among other factors may contribute to the biocontrol properties of the strain.

Genome wide transcription start sites analysis of Xanthomonas campestris pv. campestris B100 with insights into the gum gene cluster directing the biosynthesis of the exopolysaccharide xanthan

Xanthomonas campestris pv. campestris (Xcc) is the major producer of the exopolysaccharide xanthan, the commercially most important natural polysaccharide of microbial origin. The current work provides deeper insights into the yet uncharacterized transcriptomic features of the xanthan producing strain Xcc-B100. Towards this goal, RNA sequencing of a library based on the selective enrichment of the 5' ends of native transcripts was performed. This approach resulted in the genome wide identification of 3067 transcription start sites (TSSs) that were further classified based on their genomic positions. Among them, 1545 mapped upstream of an actively transcribed CDS and 1363 were classified as novel TSSs representing antisense, internal, and TSSs belonging to previously unidentified genomic features. Analyzing the transcriptional strength of primary and antisense TSSs revealed that in some instances antisense transcription seemed to be initiated at a higher level than its sense counterpart. Mapping the exact positions of TSSs aided in the identification of promoter consensus motifs, ribosomal binding sites, and enhanced the genome annotation of 159 in silico predicted translational start (TLS) sites. The global view on length distribution of the 5' untranslated regions (5'-UTRs) deduced from the data pointed to the occurrence of leaderless transcripts and transcripts with unusually long 5'-UTRs, in addition to identifying seven putative riboswitch elements for Xcc-B100. Concerning the biosynthesis of xanthan, we focused on the transcriptional organization of the gum gene cluster. Under the conditions tested, we present evidence for a complex transcription pattern of the gum genes with multiple TSSs and an obvious considerable role of antisense transcription. The gene gumB, encoding an outer membrane xanthan exporter, is presented here as an example for genes that possessed a strong antisense TSS.

Phenotypic Heterogeneity Affects Stenotrophomonas maltophilia K279a Colony Morphotypes and β-Lactamase Expression

Phenotypic heterogeneity at the cellular level in response to various stresses, e.g., antibiotic treatment has been reported for a number of bacteria. In a clonal population, cell-to-cell variation may result in phenotypic heterogeneity that is a mechanism to survive changing environments including antibiotic therapy. Stenotrophomonas maltophilia has been frequently isolated from cystic fibrosis patients, can cause numerous infections in other organs and tissues, and is difficult to treat due to antibiotic resistances. S. maltophilia K279a produces the L1 and L2 β-lactamases in response to β-lactam treatment. Here we report that the patient isolate S. maltophilia K279a diverges into cellular subpopulations with distinct but reversible morphotypes of small and big colonies when challenged with ampicillin. This observation is consistent with the formation of elongated chains of bacteria during exponential growth phase and the occurrence of mainly rod-shaped cells in liquid media. RNA-seq analysis of small versus big colonies revealed differential regulation of at least seven genes among the colony morphotypes. Among those, bla L1 and bla L2 were transcriptionally the most strongly upregulated genes. Promoter fusions of bla L1 and bla L2 genes indicated that expression of both genes is also subject to high levels of phenotypic heterogeneous expression on a single cell level. Additionally, the comE homolog was found to be differentially expressed in homogenously versus heterogeneously bla L2 expressing cells as identified by RNA-seq analysis. Overexpression of comE in S. maltophilia K279a reduced the level of cells that were in a bla L2-ON mode to 1% or lower. Taken together, our data provide strong evidence that S. maltophilia K279a populations develop phenotypic heterogeneity in an ampicillin challenged model. This cellular variability is triggered by regulation networks including bla L1, bla L2, and comE.

Alignment of Short Reads: A Crucial Step for Application of Next-Generation Sequencing Data in Precision Medicine

Precision medicine or personalized medicine has been proposed as a modernized and promising medical strategy. Genetic variants of patients are the key information for implementation of precision medicine. Next-generation sequencing (NGS) is an emerging technology for deciphering genetic variants. Alignment of raw reads to a reference genome is one of the key steps in NGS data analysis. Many algorithms have been developed for alignment of short read sequences since 2008. Users have to make a decision on which alignment algorithm to use in their studies. Selection of the right alignment algorithm determines not only the alignment algorithm but also the set of suitable parameters to be used by the algorithm. Understanding these algorithms helps in selecting the appropriate alignment algorithm for different applications in precision medicine. Here, we review current available algorithms and their major strategies such as seed-and-extend and q-gram filter. We also discuss the challenges in current alignment algorithms, including alignment in multiple repeated regions, long reads alignment and alignment facilitated with known genetic variants.

Identification of two mutations increasing the methanol tolerance of Corynebacterium glutamicum

Background

Methanol is present in most ecosystems and may also occur in industrial applications, e.g. as an impurity of carbon sources such as technical glycerol. Methanol often inhibits growth of bacteria, thus, methanol tolerance may limit fermentative production processes.

Results

The methanol tolerance of the amino acid producing soil bacterium Corynebacterium glutamicum was improved by experimental evolution in the presence of methanol. The resulting strain Tol1 exhibited significantly increased growth rates in the presence of up to 1 M methanol. However, neither transcriptional changes nor increased enzyme activities of the linear methanol oxidation pathway were observed, which was in accordance with the finding that tolerance to the downstream metabolites formaldehyde and formate was not improved. Genome sequence analysis of strain Tol1 revealed two point mutations potentially relevant to enhanced methanol tolerance: one leading to the amino acid exchange A165T of O-acetylhomoserine sulfhydrolase MetY and the other leading to shortened CoA transferase Cat (Q342*). Introduction of either mutation into the genome of C. glutamicum wild type increased methanol tolerance and introduction of both mutations into C. glutamicum was sufficient to achieve methanol tolerance almost indistinguishable from that of strain Tol1.

Conclusion

The methanol tolerance of C. glutamicum can be increased by two point mutations leading to amino acid exchange of O-acetylhomoserine sulfhydrolase MetY and shortened CoA transferase Cat. Introduction of these mutations into producer strains may be helpful when using carbon sources containing methanol as component or impurity.

Electronic supplementary material

The online version of this article (doi:10.1186/s12866-015-0558-6) contains supplementary material, which is available to authorized users.

Boosting the FM-Index on the GPU: Effective Techniques to Mitigate Random Memory Access

The recent advent of high-throughput sequencing machines producing big amounts of short reads has boosted the interest in efficient string searching techniques. As of today, many mainstream sequence alignment software tools rely on a special data structure, called the FM-index, which allows for fast exact searches in large genomic references. However, such searches translate into a pseudo-random memory access pattern, thus making memory access the limiting factor of all computation-efficient implementations, both on CPUs and GPUs. Here, we show that several strategies can be put in place to remove the memory bottleneck on the GPU: more compact indexes can be implemented by having more threads work cooperatively on larger memory blocks, and a k-step FM-index can be used to further reduce the number of memory accesses. The combination of those and other optimisations yields an implementation that is able to process about two Gbases of queries per second on our test platform, being about 8 × faster than a comparable multi-core CPU version, and about 3 × to 5 × faster than the FM-index implementation on the GPU provided by the recently announced Nvidia NVBIO bioinformatics library.

Direct DNA Extraction from Mycobacterium tuberculosis Frozen Stocks as a Reculture-Independent Approach to Whole-Genome Sequencing

Culturing before DNA extraction represents a major time-consuming step in whole-genome sequencing of slow-growing bacteria, such as Mycobacterium tuberculosis. We report a workflow to extract DNA from frozen isolates without reculturing. Prepared libraries and sequence data were comparable with results from recultured aliquots of the same stocks.

Engineering L-arabinose metabolism in triacylglycerol-producing Rhodococcus opacus for lignocellulosic fuel production

Advanced biofuels from lignocellulosic biomass have been considered as a potential solution for the issues of energy sustainability and environmental protection. Triacylglycerols (TAGs) are potential precursors for the production of lipid-based liquid biofuels. Rhodococcus opacus PD630 can accumulate large amounts of TAGs when grown under physiological conditions of high carbon and low nitrogen. However, R. opacus PD630 does not utilize the sugar L-arabinose present in lignocellulosic hydrolysates. Here, we report the engineering of R. opacus to produce TAGs on L-arabinose. We constructed a plasmid (pASC8057) harboring araB, araD and araA genes derived from a Streptomyces bacterium, and introduced the genes into R. opacus PD630. One of the engineered strains, MITAE-348, was capable of growing on high concentrations (up to 100 g/L) of L-arabinose. MITAE-348 was grown in a defined medium containing 16 g/L L-arabinose or a mixture of 8 g/L L-arabinose and 8 g/L D-glucose. In a stationary phase occurring 3 days post-inoculation, the strain was able to completely utilize the sugar, and yielded 2.0 g/L for L-arabinose and 2.2 g/L for L-arabinose/D-glucose of TAGs, corresponding to 39.7% or 42.0%, respectively, of the cell dry weight.

Detoxifying Escherichia coli for endotoxin-free production of recombinant proteins

BACKGROUND

Lipopolysaccharide (LPS), also referred to as endotoxin, is the major constituent of the outer leaflet of the outer membrane of virtually all Gram-negative bacteria. The lipid A moiety, which anchors the LPS molecule to the outer membrane, acts as a potent agonist for Toll-like receptor 4/myeloid differentiation factor 2-mediated pro-inflammatory activity in mammals and, thus, represents the endotoxic principle of LPS. Recombinant proteins, commonly manufactured in Escherichia coli, are generally contaminated with endotoxin. Removal of bacterial endotoxin from recombinant therapeutic proteins is a challenging and expensive process that has been necessary to ensure the safety of the final product.

RESULTS

As an alternative strategy for common endotoxin removal methods, we have developed a series of E. coli strains that are able to grow and express recombinant proteins with the endotoxin precursor lipid IVA as the only LPS-related molecule in their outer membranes. Lipid IVA does not trigger an endotoxic response in humans typical of bacterial LPS chemotypes. Hence the engineered cells themselves, and the purified proteins expressed within these cells display extremely low endotoxin levels.

CONCLUSIONS

This paper describes the preparation and characterization of endotoxin-free E. coli strains, and demonstrates the direct production of recombinant proteins with negligible endotoxin contamination.

Evolutionary history and global spread of the Mycobacterium tuberculosis Beijing lineage

Mycobacterium tuberculosis strains of the Beijing lineage are globally distributed and are associated with the massive spread of multidrug-resistant (MDR) tuberculosis in Eurasia. Here we reconstructed the biogeographical structure and evolutionary history of this lineage by genetic analysis of 4,987 isolates from 99 countries and whole-genome sequencing of 110 representative isolates. We show that this lineage initially originated in the Far East, from where it radiated worldwide in several waves. We detected successive increases in population size for this pathogen over the last 200 years, practically coinciding with the Industrial Revolution, the First World War and HIV epidemics. Two MDR clones of this lineage started to spread throughout central Asia and Russia concomitantly with the collapse of the public health system in the former Soviet Union. Mutations identified in genes putatively under positive selection and associated with virulence might have favored the expansion of the most successful branches of the lineage.

Transcriptome analysis of thermophilic methylotrophic Bacillus methanolicus MGA3 using RNA-sequencing provides detailed insights into its previously uncharted transcriptional landscape

Background

Bacillus methanolicus MGA3 is a thermophilic, facultative ribulose monophosphate (RuMP) cycle methylotroph. Together with its ability to produce high yields of amino acids, the relevance of this microorganism as a promising candidate for biotechnological applications is evident. The B. methanolicus MGA3 genome consists of a 3,337,035 nucleotides (nt) circular chromosome, the 19,174 nt plasmid pBM19 and the 68,999 nt plasmid pBM69. 3,218 protein-coding regions were annotated on the chromosome, 22 on pBM19 and 82 on pBM69. In the present study, the RNA-seq approach was used to comprehensively investigate the transcriptome of B. methanolicus MGA3 in order to improve the genome annotation, identify novel transcripts, analyze conserved sequence motifs involved in gene expression and reveal operon structures. For this aim, two different cDNA library preparation methods were applied: one which allows characterization of the whole transcriptome and another which includes enrichment of primary transcript 5′-ends.

Results

Analysis of the primary transcriptome data enabled the detection of 2,167 putative transcription start sites (TSSs) which were categorized into 1,642 TSSs located in the upstream region (5′-UTR) of known protein-coding genes and 525 TSSs of novel antisense, intragenic, or intergenic transcripts. Firstly, 14 wrongly annotated translation start sites (TLSs) were corrected based on primary transcriptome data. Further investigation of the identified 5′-UTRs resulted in the detailed characterization of their length distribution and the detection of 75 hitherto unknown cis-regulatory RNA elements. Moreover, the exact TSSs positions were utilized to define conserved sequence motifs for translation start sites, ribosome binding sites and promoters in B. methanolicus MGA3. Based on the whole transcriptome data set, novel transcripts, operon structures and mRNA abundances were determined. The analysis of the operon structures revealed that almost half of the genes are transcribed monocistronically (940), whereas 1,164 genes are organized in 381 operons. Several of the genes related to methylotrophy had highly abundant transcripts.

Conclusion

The extensive insights into the transcriptional landscape of B. methanolicus MGA3, gained in this study, represent a valuable foundation for further comparative quantitative transcriptome analyses and possibly also for the development of molecular biology tools which at present are very limited for this organism.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1239-4) contains supplementary material, which is available to authorized users.

Metabolic engineering of an ATP-neutral Embden-Meyerhof-Parnas pathway in Corynebacterium glutamicum: growth restoration by an adaptive point mutation in NADH dehydrogenase

Corynebacterium glutamicum uses the Embden-Meyerhof-Parnas pathway of glycolysis and gains 2 mol of ATP per mol of glucose by substrate-level phosphorylation (SLP). To engineer glycolysis without net ATP formation by SLP, endogenous phosphorylating NAD-dependent glyceraldehyde-3-phosphate dehydrogenase (GAPDH) was replaced by nonphosphorylating NADP-dependent glyceraldehyde-3-phosphate dehydrogenase (GapN) from Clostridium acetobutylicum, which irreversibly converts glyceraldehyde-3-phosphate (GAP) to 3-phosphoglycerate (3-PG) without generating ATP. As shown recently (S. Takeno, R. Murata, R. Kobayashi, S. Mitsuhashi, and M. Ikeda, Appl Environ Microbiol 76:7154-7160, 2010, http://dx.doi.org/10.1128/AEM.01464-10), this ATP-neutral, NADPH-generating glycolytic pathway did not allow for the growth of Corynebacterium glutamicum with glucose as the sole carbon source unless hitherto unknown suppressor mutations occurred; however, these mutations were not disclosed. In the present study, a suppressor mutation was identified, and it was shown that heterologous expression of udhA encoding soluble transhydrogenase from Escherichia coli partly restored growth, suggesting that growth was inhibited by NADPH accumulation. Moreover, genome sequence analysis of second-site suppressor mutants that were able to grow faster with glucose revealed a single point mutation in the gene of non-proton-pumping NADH:ubiquinone oxidoreductase (NDH-II) leading to the amino acid change D213G, which was shared by these suppressor mutants. Since related NDH-II enzymes accepting NADPH as the substrate possess asparagine or glutamine residues at this position, D213G, D213N, and D213Q variants of C. glutamicum NDH-II were constructed and were shown to oxidize NADPH in addition to NADH. Taking these findings together, ATP-neutral glycolysis by the replacement of endogenous NAD-dependent GAPDH with NADP-dependent GapN became possible via oxidation of NADPH formed in this pathway by mutant NADPH-accepting NDH-II(D213G) and thus by coupling to electron transport phosphorylation (ETP).

CLAST: CUDA implemented large-scale alignment search tool

BACKGROUND

Metagenomics is a powerful methodology to study microbial communities, but it is highly dependent on nucleotide sequence similarity searching against sequence databases. Metagenomic analyses with next-generation sequencing technologies produce enormous numbers of reads from microbial communities, and many reads are derived from microbes whose genomes have not yet been sequenced, limiting the usefulness of existing sequence similarity search tools. Therefore, there is a clear need for a sequence similarity search tool that can rapidly detect weak similarity in large datasets.

RESULTS

We developed a tool, which we named CLAST (CUDA implemented large-scale alignment search tool), that enables analyses of millions of reads and thousands of reference genome sequences, and runs on NVIDIA Fermi architecture graphics processing units. CLAST has four main advantages over existing alignment tools. First, CLAST was capable of identifying sequence similarities ~80.8 times faster than BLAST and 9.6 times faster than BLAT. Second, CLAST executes global alignment as the default (local alignment is also an option), enabling CLAST to assign reads to taxonomic and functional groups based on evolutionarily distant nucleotide sequences with high accuracy. Third, CLAST does not need a preprocessed sequence database like Burrows-Wheeler Transform-based tools, and this enables CLAST to incorporate large, frequently updated sequence databases. Fourth, CLAST requires <2 GB of main memory, making it possible to run CLAST on a standard desktop computer or server node.

CONCLUSIONS

CLAST achieved very high speed (similar to the Burrows-Wheeler Transform-based Bowtie 2 for long reads) and sensitivity (equal to BLAST, BLAT, and FR-HIT) without the need for extensive database preprocessing or a specialized computing platform. Our results demonstrate that CLAST has the potential to be one of the most powerful and realistic approaches to analyze the massive amount of sequence data from next-generation sequencing technologies.

Massively parallel read mapping on GPUs with the q-group index and PEANUT

We present the q-group index, a novel data structure for read mapping tailored towards graphics processing units (GPUs) with a small memory footprint and efficient parallel algorithms for querying and building. On top of the q-group index we introduce PEANUT, a highly parallel GPU-based read mapper. PEANUT provides the possibility to output both the best hits or all hits of a read. Our benchmarks show that PEANUT outperforms other state-of-the-art read mappers in terms of speed while maintaining or slightly increasing precision, recall and sensitivity.

Whole-genome-based Mycobacterium tuberculosis surveillance: a standardized, portable, and expandable approach

Whole-genome sequencing (WGS) allows for effective tracing of Mycobacterium tuberculosis complex (MTBC) (tuberculosis pathogens) transmission. However, it is difficult to standardize and, therefore, is not yet employed for interlaboratory prospective surveillance. To allow its widespread application, solutions for data standardization and storage in an easily expandable database are urgently needed. To address this question, we developed a core genome multilocus sequence typing (cgMLST) scheme for clinical MTBC isolates using the Ridom SeqSphere(+) software, which transfers the genome-wide single nucleotide polymorphism (SNP) diversity into an allele numbering system that is standardized, portable, and not computationally intensive. To test its performance, we performed WGS analysis of 26 isolates with identical IS6110 DNA fingerprints and spoligotyping patterns from a longitudinal outbreak in the federal state of Hamburg, Germany (notified between 2001 and 2010). The cgMLST approach (3,041 genes) discriminated the 26 strains with a resolution comparable to that of SNP-based WGS typing (one major cluster of 22 identical or closely related and four outlier isolates with at least 97 distinct SNPs or 63 allelic variants). Resulting tree topologies are highly congruent and grouped the isolates in both cases analogously. Our data show that SNP- and cgMLST-based WGS analyses facilitate high-resolution discrimination of longitudinal MTBC outbreaks. cgMLST allows for a meaningful epidemiological interpretation of the WGS genotyping data. It enables standardized WGS genotyping for epidemiological investigations, e.g., on the regional public health office level, and the creation of web-accessible databases for global TB surveillance with an integrated early warning system.

ReadXplorer--visualization and analysis of mapped sequences

MOTIVATION

Fast algorithms and well-arranged visualizations are required for the comprehensive analysis of the ever-growing size of genomic and transcriptomic next-generation sequencing data.

RESULTS

ReadXplorer is a software offering straightforward visualization and extensive analysis functions for genomic and transcriptomic DNA sequences mapped on a reference. A unique specialty of ReadXplorer is the quality classification of the read mappings. It is incorporated in all analysis functions and displayed in ReadXplorer's various synchronized data viewers for (i) the reference sequence, its base coverage as (ii) normalizable plot and (iii) histogram, (iv) read alignments and (v) read pairs. ReadXplorer's analysis capability covers RNA secondary structure prediction, single nucleotide polymorphism and deletion-insertion polymorphism detection, genomic feature and general coverage analysis. Especially for RNA-Seq data, it offers differential gene expression analysis, transcription start site and operon detection as well as RPKM value and read count calculations. Furthermore, ReadXplorer can combine or superimpose coverage of different datasets.

AVAILABILITY AND IMPLEMENTATION

ReadXplorer is available as open-source software at http://www.readxplorer.org along with a detailed manual.

Improving the genome annotation of the acarbose producer Actinoplanes sp. SE50/110 by sequencing enriched 5'-ends of primary transcripts

Actinoplanes sp. SE50/110 is the producer of the alpha-glucosidase inhibitor acarbose, which is an economically relevant and potent drug in the treatment of type-2 diabetes mellitus. In this study, we present the detection of transcription start sites on this genome by sequencing enriched 5'-ends of primary transcripts. Altogether, 1427 putative transcription start sites were initially identified. With help of the annotated genome sequence, 661 transcription start sites were found to belong to the leader region of protein-coding genes with the surprising result that roughly 20% of these genes rank among the class of leaderless transcripts. Next, conserved promoter motifs were identified for protein-coding genes with and without leader sequences. The mapped transcription start sites were finally used to improve the annotation of the Actinoplanes sp. SE50/110 genome sequence. Concerning protein-coding genes, 41 translation start sites were corrected and 9 novel protein-coding genes could be identified. In addition to this, 122 previously undetermined non-coding RNA (ncRNA) genes of Actinoplanes sp. SE50/110 were defined. Focusing on antisense transcription start sites located within coding genes or their leader sequences, it was discovered that 96 of those ncRNA genes belong to the class of antisense RNA (asRNA) genes. The remaining 26 ncRNA genes were found outside of known protein-coding genes. Four chosen examples of prominent ncRNA genes, namely the transfer messenger RNA gene ssrA, the ribonuclease P class A RNA gene rnpB, the cobalamin riboswitch RNA gene cobRS, and the selenocysteine-specific tRNA gene selC, are presented in more detail. This study demonstrates that sequencing of enriched 5'-ends of primary transcripts and the identification of transcription start sites are valuable tools for advanced genome annotation of Actinoplanes sp. SE50/110 and most probably also for other bacteria.

Comprehensive analysis of the Corynebacterium glutamicum transcriptome using an improved RNAseq technique

Background

The use of RNAseq to resolve the transcriptional organization of an organism was established in recent years and also showed the complexity and dynamics of bacterial transcriptomes. The aim of this study was to comprehensively investigate the transcriptome of the industrially relevant amino acid producer and model organism Corynebacterium glutamicum by RNAseq in order to improve its genome annotation and to describe important features for transcription and translation.

Results

RNAseq data sets were obtained by two methods, one that focuses on 5′-ends of primary transcripts and another that provides the overall transcriptome with an improved resolution of 3′-ends of transcripts. Subsequent data analysis led to the identification of more than 2,000 transcription start sites (TSSs), the definition of 5′-UTRs (untranslated regions) for annotated protein-coding genes, operon structures and many novel transcripts located between or in antisense orientation to protein-coding regions. Interestingly, a high number of mRNAs (33%) is transcribed as leaderless transcripts. From the data, consensus promoter and ribosome binding site (RBS) motifs were identified and it was shown that the majority of genes in C. glutamicum are transcribed monocistronically, but operons containing up to 16 genes are also present.

Conclusions

The comprehensive transcriptome map of C. glutamicum established in this study represents a major step forward towards a complete definition of genetic elements (e.g. promoter regions, gene starts and stops, 5′-UTRs, RBSs, transcript starts and ends) and provides the ideal basis for further analyses on transcriptional regulatory networks in this organism. The methods developed are easily applicable for other bacteria and have the potential to be used also for quantification of transcriptomes, replacing microarrays in the near future.

Whole genome sequencing reveals complex evolution patterns of multidrug-resistant Mycobacterium tuberculosis Beijing strains in patients

Multidrug-resistant (MDR) Mycobacterium tuberculosis complex (MTBC) strains represent a major threat for tuberculosis (TB) control. Treatment of MDR-TB patients is long and less effective, resulting in a significant number of treatment failures. The development of further resistances leads to extensively drug-resistant (XDR) variants. However, data on the individual reasons for treatment failure, e.g. an induced mutational burst, and on the evolution of bacteria in the patient are only sparsely available. To address this question, we investigated the intra-patient evolution of serial MTBC isolates obtained from three MDR-TB patients undergoing longitudinal treatment, finally leading to XDR-TB. Sequential isolates displayed identical IS6110 fingerprint patterns, suggesting the absence of exogenous re-infection. We utilized whole genome sequencing (WGS) to screen for variations in three isolates from Patient A and four isolates from Patient B and C, respectively. Acquired polymorphisms were subsequently validated in up to 15 serial isolates by Sanger sequencing. We determined eight (Patient A) and nine (Patient B) polymorphisms, which occurred in a stepwise manner during the course of the therapy and were linked to resistance or a potential compensatory mechanism. For both patients, our analysis revealed the long-term co-existence of clonal subpopulations that displayed different drug resistance allele combinations. Out of these, the most resistant clone was fixed in the population. In contrast, baseline and follow-up isolates of Patient C were distinguished each by eleven unique polymorphisms, indicating an exogenous re-infection with an XDR strain not detected by IS6110 RFLP typing. Our study demonstrates that intra-patient microevolution of MDR-MTBC strains under longitudinal treatment is more complex than previously anticipated. However, a mutator phenotype was not detected. The presence of different subpopulations might confound phenotypic and molecular drug resistance tests. Furthermore, high resolution WGS analysis is necessary to accurately detect exogenous re-infection as classical genotyping lacks discriminatory power in high incidence settings.

Discovering epistasis in large scale genetic association studies by exploiting graphics cards

Despite the enormous investments made in collecting DNA samples and generating germline variation data across thousands of individuals in modern genome-wide association studies (GWAS), progress has been frustratingly slow in explaining much of the heritability in common disease. Today's paradigm of testing independent hypotheses on each single nucleotide polymorphism (SNP) marker is unlikely to adequately reflect the complex biological processes in disease risk. Alternatively, modeling risk as an ensemble of SNPs that act in concert in a pathway, and/or interact non-additively on log risk for example, may be a more sensible way to approach gene mapping in modern studies. Implementing such analyzes genome-wide can quickly become intractable due to the fact that even modest size SNP panels on modern genotype arrays (500k markers) pose a combinatorial nightmare, require tens of billions of models to be tested for evidence of interaction. In this article, we provide an in-depth analysis of programs that have been developed to explicitly overcome these enormous computational barriers through the use of processors on graphics cards known as Graphics Processing Units (GPU). We include tutorials on GPU technology, which will convey why they are growing in appeal with today's numerical scientists. One obvious advantage is the impressive density of microprocessor cores that are available on only a single GPU. Whereas high end servers feature up to 24 Intel or AMD CPU cores, the latest GPU offerings from nVidia feature over 2600 cores. Each compute node may be outfitted with up to 4 GPU devices. Success on GPUs varies across problems. However, epistasis screens fare well due to the high degree of parallelism exposed in these problems. Papers that we review routinely report GPU speedups of over two orders of magnitude (>100x) over standard CPU implementations.

Comprehensive discovery and characterization of small RNAs in Corynebacterium glutamicum ATCC 13032

Background

Recent discoveries on bacterial transcriptomes gave evidence that small RNAs (sRNAs) have important regulatory roles in prokaryotic cells. Modern high-throughput sequencing approaches (RNA-Seq) enable the most detailed view on transcriptomes offering an unmatched comprehensiveness and single-base resolution. Whole transcriptome data obtained by RNA-Seq can be used to detect and characterize all transcript species, including small RNAs. Here, we describe an RNA-Seq approach for comprehensive detection and characterization of small RNAs from Corynebacterium glutamicum, an actinobacterium of high industrial relevance and model organism for medically important Corynebacterianeae, such as C. diphtheriae and Mycobacterium tuberculosis.

Results

In our RNA-Seq approach, total RNA from C. glutamicum ATCC 13032 was prepared from cultures grown in minimal medium at exponential growth or challenged by physical (heat shock, cold shock) or by chemical stresses (diamide, H₂O₂, NaCl) at this time point. Total RNA samples were pooled and sequencing libraries were prepared from the isolated small RNA fraction. High throughput short read sequencing and mapping yielded over 800 sRNA genes. By determining their 5′- and 3′-ends and inspection of their locations, these potential sRNA genes were classified into UTRs of mRNAs (316), cis-antisense sRNAs (543), and trans-encoded sRNAs (262). For 77 of trans-encoded sRNAs significant sequence and secondary structure conservation was found by a computational approach using a whole genome alignment with the closely related species C. efficiens YS-314 and C. diphtheriae NCTC 13129. Three selected trans-encoded sRNAs were characterized by Northern blot analysis and stress-specific transcript patterns were found.

Conclusions

The study showed comparable numbers of sRNAs known from genome-wide surveys in other bacteria. In detail, our results give deep insight into the comprehensive equipment of sRNAs in C. glutamicum and provide a sound basis for further studies concerning the functions of these sRNAs.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-14-714) contains supplementary material, which is available to authorized users.

New insights into Chlamydomonas reinhardtii hydrogen production processes by combined microarray/RNA-seq transcriptomics

Hydrogen production with Chlamydomonas reinhardtii induced by sulphur starvation is a multiphase process while the cell internal metabolism is completely remodelled. The first cellular response is characterized by induction of genes with regulatory functions, followed by a total remodelling of the metabolism to provide reduction equivalents for cellular processes. We were able to characterize all major processes that provide energy and reduction equivalents during hydrogen production. Furthermore, C. reinhardtii showed a strong transcript increase for gene models responsible for stress response and detoxification of oxygen radicals. Finally, we were able to determine potential bottlenecks and target genes for manipulation to increase hydrogen production or to prolong the hydrogen production phase. The investigation of transcriptomic changes during the time course of hydrogen production in C. reinhardtii with microarrays and RNA-seq revealed new insights into the regulation and remodelling of the cell internal metabolism. Both methods showed a good correlation. The microarray platform can be used as a reliable standard tool for routine gene expression analysis. RNA-seq additionally allowed a detailed time-dependent study of gene expression and determination of new genes involved in the hydrogen production process.

Complex regulation of the phosphoenolpyruvate carboxykinase gene pck and characterization of its GntR-type regulator IolR as a repressor of myo-inositol utilization genes in Corynebacterium glutamicum

DNA affinity chromatography with the promoter region of the Corynebacterium glutamicum pck gene, encoding phosphoenolpyruvate carboxykinase, led to the isolation of four transcriptional regulators, i.e., RamA, GntR1, GntR2, and IolR. Determination of the phosphoenolpyruvate carboxykinase activity of the ΔramA, ΔgntR1 ΔgntR2, and ΔiolR deletion mutants indicated that RamA represses pck during growth on glucose about 2-fold, whereas GntR1, GntR2, and IolR activate pck expression about 2-fold irrespective of whether glucose or acetate served as the carbon source. The DNA binding sites of the four regulators in the pck promoter region were identified and their positions correlated with the predicted functions as repressor or activators. The iolR gene is located upstream and in a divergent orientation with respect to a iol gene cluster, encoding proteins involved in myo-inositol uptake and degradation. Comparative DNA microarray analysis of the ΔiolR mutant and the parental wild-type strain revealed strongly (>100-fold) elevated mRNA levels of the iol genes in the mutant, indicating that the primary function of IolR is the repression of the iol genes. IolR binding sites were identified in the promoter regions of iolC, iolT1, and iolR. IolR therefore is presumably subject to negative autoregulation. A consensus DNA binding motif (5'-KGWCHTRACA-3') which corresponds well to those of other GntR-type regulators of the HutC family was identified. Taken together, our results disclose a complex regulation of the pck gene in C. glutamicum and identify IolR as an efficient repressor of genes involved in myo-inositol catabolism of this organism.

Benchmarking short sequence mapping tools

BACKGROUND

The development of next-generation sequencing instruments has led to the generation of millions of short sequences in a single run. The process of aligning these reads to a reference genome is time consuming and demands the development of fast and accurate alignment tools. However, the current proposed tools make different compromises between the accuracy and the speed of mapping. Moreover, many important aspects are overlooked while comparing the performance of a newly developed tool to the state of the art. Therefore, there is a need for an objective evaluation method that covers all the aspects. In this work, we introduce a benchmarking suite to extensively analyze sequencing tools with respect to various aspects and provide an objective comparison.

RESULTS

We applied our benchmarking tests on 9 well known mapping tools, namely, Bowtie, Bowtie2, BWA, SOAP2, MAQ, RMAP, GSNAP, Novoalign, and mrsFAST (mrFAST) using synthetic data and real RNA-Seq data. MAQ and RMAP are based on building hash tables for the reads, whereas the remaining tools are based on indexing the reference genome. The benchmarking tests reveal the strengths and weaknesses of each tool. The results show that no single tool outperforms all others in all metrics. However, Bowtie maintained the best throughput for most of the tests while BWA performed better for longer read lengths. The benchmarking tests are not restricted to the mentioned tools and can be further applied to others.

CONCLUSION

The mapping process is still a hard problem that is affected by many factors. In this work, we provided a benchmarking suite that reveals and evaluates the different factors affecting the mapping process. Still, there is no tool that outperforms all of the others in all the tests. Therefore, the end user should clearly specify his needs in order to choose the tool that provides the best results.