Curare and GenExVis: a versatile toolkit for analyzing and visualizing RNA-Seq data.

Even though high-throughput transcriptome sequencing is routinely performed in many laboratories, computational analysis of such data remains a cumbersome process often executed manually, hence error-prone and lacking reproducibility. For corresponding data processing, we introduce Curare, an easy-to-use yet versatile workflow builder for analyzing high-throughput RNA-Seq data focusing on differential gene expression experiments. Data analysis with Curare is customizable and subdivided into preprocessing, quality control, mapping, and downstream analysis stages, providing multiple options for each step while ensuring the reproducibility of the workflow. For a fast and straightforward exploration and visualization of differential gene expression results, we provide the gene expression visualizer software GenExVis. GenExVis can create various charts and tables from simple gene expression tables and DESeq2 results without the requirement to upload data or install software packages. In combination, Curare and GenExVis provide a comprehensive software environment that supports the entire data analysis process, from the initial handling of raw RNA-Seq data to the final DGE analyses and result visualizations, thereby significantly easing data processing and subsequent interpretation.

Expansion and re-classification of the extracytoplasmic function (ECF) σ factor family.

Extracytoplasmic function σ factors (ECFs) represent one of the major bacterial signal transduction mechanisms in terms of abundance, diversity and importance, particularly in mediating stress responses. Here, we performed a comprehensive phylogenetic analysis of this protein family by scrutinizing all proteins in the NCBI database. As a result, we identified an average of ∼10 ECFs per bacterial genome and 157 phylogenetic ECF groups that feature a conserved genetic neighborhood and a similar regulation mechanism. Our analysis expands previous classification efforts ∼50-fold, enriches many original ECF groups with previously unclassified proteins and identifies 22 entirely new ECF groups. The ECF groups are hierarchically related to each other and are further composed of subgroups with closely related sequences. This two-tiered classification allows for the accurate prediction of common promoter motifs and the inference of putative regulatory mechanisms across subgroups composing an ECF group. This comprehensive, high-resolution description of the phylogenetic distribution of the ECF family, together with the massive expansion of classified ECF sequences and an openly accessible data repository called 'ECF Hub' (https://www.computational.bio.uni-giessen.de/ecfhub), will serve as a powerful hypothesis-generator to guide future research in the field.

ReadXplorer 2-detailed read mapping analysis and visualization from one single source.

MOTIVATION: The vast amount of already available and currently generated read mapping data requires comprehensive visualization, and should benefit from bioinformatics tools offering a wide spectrum of analysis functionality from just one source. Appropriate handling of multiple mapped reads during mapping analyses remains an issue that demands improvement. RESULTS: The capabilities of the read mapping analysis and visualization tool ReadXplorer were vastly enhanced. Here, we present an even finer granulated read mapping classification, improving the level of detail for analyses and visualizations. The spectrum of automatic analysis functions has been broadened to include genome rearrangement detection as well as correlation analysis between two mapping data sets. Existing functions were refined and enhanced, namely the computation of differentially expressed genes, the read count and normalization analysis and the transcription start site detection. Additionally, ReadXplorer 2 features a highly improved support for large eukaryotic data sets and a command line version, enabling its integration into workflows. Finally, the new version is now able to display any kind of tabular results from other bioinformatics tools. AVAILABILITY AND IMPLEMENTATION: http://www.readxplorer.org CONTACT: readxplorer@computational.bio.uni-giessen.deSupplementary information: Supplementary data are available at Bioinformatics online.

Microbial diversity in different compartments of an aquaponics system.

Aquaponics is a solution for sustainable production of fish and plants in a single semi-closed system, where nutrient-rich water from the aquaculture provides nutrients for plant growth. We examined the microbial communities within an experimental aquaponics system. Whereas the fish feces contained a separate community dominated by bacteria of the genus Cetobacterium, the samples from plant roots, biofilter, and periphyton were more similar to each other, while the communities were more diverse. Detailed examination of the data gave the first indications to functional groups of organisms in the different compartments of the aquaponic system. As other nitrifiers other than members of the genus Nitrospira were only present at low numbers, it was anticipated that Nitrospirae may perform the nitrification process in the biofilm.

The Sinorhizobium fredii HH103 Genome: A Comparative Analysis With S. fredii Strains Differing in Their Symbiotic Behavior With Soybean.

Sinorhizobium fredii HH103 is a fast-growing rhizobial strain infecting a broad range of legumes including both American and Asiatic soybeans. In this work, we present the sequencing and annotation of the HH103 genome (7.25 Mb), consisting of one chromosome and six plasmids and representing the structurally most complex sinorhizobial genome sequenced so far. Comparative genomic analyses of S. fredii HH103 with strains USDA257 and NGR234 showed that the core genome of these three strains contains 4,212 genes (61.7% of the HH103 genes). Synteny plot analysis revealed that the much larger chromosome of USDA257 (6.48 Mb) is colinear to the HH103 (4.3 Mb) and NGR324 chromosomes (3.9 Mb). An additional region of the USDA257 chromosome of about 2 Mb displays similarity to plasmid pSfHH103e. Remarkable differences exist between HH103 and NGR234 concerning nod genes, flavonoid effect on surface polysaccharide production, and quorum-sensing systems. Furthermore a number of protein secretion systems have been found. Two genes coding for putative type III-secreted effectors not previously described in S. fredii, nopI and gunA, have been located on the HH103 genome. These differences could be important to understand the different symbiotic behavior of S. fredii strains HH103, USDA257, and NGR234 with soybean.

Carbohydrate-active enzymes identified by metagenomic analysis of deep-sea sediment bacteria.

Subseafloor sediment samples derived from a sediment core of 60 m length were used to enrich psychrophilic aerobic bacteria on cellulose, xylan, chitin, and starch. A variety of species belonging to Alpha- and Gammaproteobacteria and to Flavobacteria were isolated from sediment depths between 12 and 42 mbsf. Metagenomic DNA purified from the pooled enrichments was sequenced and analyzed for phylogenetic composition and presence of genes encoding carbohydrate-active enzymes. More than 200 open reading frames coding for glycoside hydrolases were identified, and more than 60 of them relevant for enzymatic degradation of lignocellulose. Four genes encoding ß-glucosidases with less than 52% identities to characterized enzymes were chosen for recombinant expression in Escherichia coli. In addition one endomannanase, two endoxylanases, and three ß-xylosidases were produced recombinantly. All genes could be actively expressed. Functional analysis revealed discrepancies and additional variability for the recombinant enzymes as compared to the sequence-based predictions.

Whole genome sequencing versus traditional genotyping for investigation of a Mycobacterium tuberculosis outbreak: a longitudinal molecular epidemiological study.

BACKGROUND: Understanding Mycobacterium tuberculosis (Mtb) transmission is essential to guide efficient tuberculosis control strategies. Traditional strain typing lacks sufficient discriminatory power to resolve large outbreaks. Here, we tested the potential of using next generation genome sequencing for identification of outbreak-related transmission chains. METHODS AND FINDINGS: During long-term (1997 to 2010) prospective population-based molecular epidemiological surveillance comprising a total of 2,301 patients, we identified a large outbreak caused by an Mtb strain of the Haarlem lineage. The main performance outcome measure of whole genome sequencing (WGS) analyses was the degree of correlation of the WGS analyses with contact tracing data and the spatio-temporal distribution of the outbreak cases. WGS analyses of the 86 isolates revealed 85 single nucleotide polymorphisms (SNPs), subdividing the outbreak into seven genome clusters (two to 24 isolates each), plus 36 unique SNP profiles. WGS results showed that the first outbreak isolates detected in 1997 were falsely clustered by classical genotyping. In 1998, one clone (termed "Hamburg clone") started expanding, apparently independently from differences in the social environment of early cases. Genome-based clustering patterns were in better accordance with contact tracing data and the geographical distribution of the cases than clustering patterns based on classical genotyping. A maximum of three SNPs were identified in eight confirmed human-to-human transmission chains, involving 31 patients. We estimated the Mtb genome evolutionary rate at 0.4 mutations per genome per year. This rate suggests that Mtb grows in its natural host with a doubling time of approximately 22 h (400 generations per year). Based on the genome variation discovered, emergence of the Hamburg clone was dated back to a period between 1993 and 1997, hence shortly before the discovery of the outbreak through epidemiological surveillance. CONCLUSIONS: Our findings suggest that WGS is superior to conventional genotyping for Mtb pathogen tracing and investigating micro-epidemics. WGS provides a measure of Mtb genome evolution over time in its natural host context.

Draft whole genome sequence of the cyanide-degrading bacterium Pseudomonas pseudoalcaligenes CECT5344.

Pseudomonas pseudoalcaligenes CECT5344 is a Gram-negative bacterium able to tolerate cyanide and to use it as the sole nitrogen source. We report here the first draft of the whole genome sequence of a P. pseudoalcaligenes strain that assimilates cyanide. Three aspects are specially emphasized in this manuscript. First, some generalities of the genome are shown and discussed in the context of other Pseudomonadaceae genomes, including genome size, G + C content, core genome and singletons among other features. Second, the genome is analysed in the context of cyanide metabolism, describing genes probably involved in cyanide assimilation, like those encoding nitrilases, and genes related to cyanide resistance, like the cio genes encoding the cyanide insensitive oxidases. Finally, the presence of genes probably involved in other processes with a great biotechnological potential like production of bioplastics and biodegradation of pollutants also is discussed.

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.

Changes in root bacterial communities associated to two different development stages of canola (Brassica napus L. var oleifera) evaluated through next-generation sequencing technology.

Crop production may benefit from plant growth-promoting bacteria. The knowledge on bacterial communities is indispensable in agricultural systems that intend to apply beneficial bacteria to improve plant health and production of crops such as canola. In this work, the diversity of root bacterial communities associated to two different developmental phases of canola (Brassica napus L.) plants was assessed through the application of new generation sequencing technology. Total bacterial DNA was extracted from root samples from two different growth states of canola (rosette and flowering). It could be shown how bacterial communities inside the roots changed with the growing stage of the canola plants. There were differences in the abundance of the genera, family, and even the phyla identified for each sample. While in both root samples Proteobacteria was the most common phylum, at the rosette stage, the most common bacteria belonged to the family Pseudomonadaceae and the genus Pseudomonas, and in the flowering stage, the Xanthomonadaceae family and the genus Xanthomonas dominated the community. This implies in a switch in the predominant bacteria in the different developmental stages of the plant, suggesting that the plant itself interferes with the associated microbial community.

Draft genome sequence of Wickerhamomyces ciferrii NRRL Y-1031 F-60-10.

Wickerhamomyces ciferrii is a microorganism characterized by the production and secretion of large amounts of acetylated sphingoid bases, in particular tetraacetyl phytosphingosine. Here, we present the 15.90-Mbp draft genome sequence of W. ciferrii NRRL Y-1031 F-60-10 generated by pyrosequencing and de novo assembly. The draft genome sequence comprising 364 contigs in 150 scaffolds was annotated and covered 6,702 protein-coding sequences. This information will contribute to the metabolic engineering of this yeast to improve the yield and spectrum of acetylated sphingoid bases in biotechnological production.

Genome sequence of the filamentous bacterium Fibrisoma limi BUZ 3T.

Fibrisoma limi strain BUZ 3(T), a Gram-negative bacterium, was isolated from coastal mud from the North Sea (Fedderwardersiel, Germany) and characterized using a polyphasic approach in 2011. The genome consists of a chromosome of about 7.5 Mb and three plasmids.

Genome sequence of the soybean symbiont Sinorhizobium fredii HH103.

Sinorhizobium fredii HH103 is a fast-growing rhizobial strain that is able to nodulate legumes that develop determinate nodules, e.g., soybean, and legumes that form nodules of the indeterminate type. Here we present the genome of HH103, which consists of one chromosome and five plasmids with a total size of 7.22 Mb.

Complete genome sequence of the hydrogenotrophic, methanogenic archaeon Methanoculleus bourgensis strain MS2(T), Isolated from a sewage sludge digester.

Methanoculleus bourgensis, of the order Methanomicrobiales, is a dominant methanogenic archaeon in many biogas-producing reactor systems fed with renewable primary products. It is capable of synthesizing methane via the hydrogenotrophic pathway utilizing hydrogen and carbon dioxide or formate as the substrates. Here we report the complete and finished genome sequence of M. bourgensis strain MS2(T), isolated from a sewage sludge digester.

Pangenomic study of Corynebacterium diphtheriae that provides insights into the genomic diversity of pathogenic isolates from cases of classical diphtheria, endocarditis, and pneumonia.

Corynebacterium diphtheriae is one of the most prominent human pathogens and the causative agent of the communicable disease diphtheria. The genomes of 12 strains isolated from patients with classical diphtheria, endocarditis, and pneumonia were completely sequenced and annotated. Including the genome of C. diphtheriae NCTC 13129, we herewith present a comprehensive comparative analysis of 13 strains and the first characterization of the pangenome of the species C. diphtheriae. Comparative genomics showed extensive synteny and revealed a core genome consisting of 1,632 conserved genes. The pangenome currently comprises 4,786 protein-coding regions and increases at an average of 65 unique genes per newly sequenced strain. Analysis of prophages carrying the diphtheria toxin gene tox revealed that the toxoid vaccine producer C. diphtheriae Park-Williams no. 8 has been lysogenized by two copies of the ω(tox)(+) phage, whereas C. diphtheriae 31A harbors a hitherto-unknown tox(+) corynephage. DNA binding sites of the tox-controlling regulator DtxR were detected by genome-wide motif searches. Comparative content analysis showed that the DtxR regulons exhibit marked differences due to gene gain, gene loss, partial gene deletion, and DtxR binding site depletion. Most predicted pathogenicity islands of C. diphtheriae revealed characteristics of horizontal gene transfer. The majority of these islands encode subunits of adhesive pili, which can play important roles in adhesion of C. diphtheriae to different host tissues. All sequenced isolates contain at least two pilus gene clusters. It appears that variation in the distributed genome is a common strategy of C. diphtheriae to establish differences in host-pathogen interactions.

Complete genome sequence, lifestyle, and multi-drug resistance of the human pathogen Corynebacterium resistens DSM 45100 isolated from blood samples of a leukemia patient.

BACKGROUND: Corynebacterium resistens was initially recovered from human infections and recognized as a new coryneform species that is highly resistant to antimicrobial agents. Bacteremia associated with this organism in immunocompromised patients was rapidly fatal as standard minocycline therapies failed. C. resistens DSM 45100 was isolated from a blood culture of samples taken from a patient with acute myelocytic leukemia. The complete genome sequence of C. resistens DSM 45100 was determined by pyrosequencing to identify genes contributing to multi-drug resistance, virulence, and the lipophilic lifestyle of this newly described human pathogen. RESULTS: The genome of C. resistens DSM 45100 consists of a circular chromosome of 2,601,311 bp in size and the 28,312-bp plasmid pJA144188. Metabolic analysis showed that the genome of C. resistens DSM 45100 lacks genes for typical sugar uptake systems, anaplerotic functions, and a fatty acid synthase, explaining the strict lipophilic lifestyle of this species. The genome encodes a broad spectrum of enzymes ensuring the availability of exogenous fatty acids for growth, including predicted virulence factors that probably contribute to fatty acid metabolism by damaging host tissue. C. resistens DSM 45100 is able to use external L-histidine as a combined carbon and nitrogen source, presumably as a result of adaptation to the hitherto unknown habitat on the human skin. Plasmid pJA144188 harbors several genes contributing to antibiotic resistance of C. resistens DSM 45100, including a tetracycline resistance region of the Tet W type known from Lactobacillus reuteri and Streptococcus suis. The tet(W) gene of pJA144188 was cloned in Corynebacterium glutamicum and was shown to confer high levels of resistance to tetracycline, doxycycline, and minocycline in vitro. CONCLUSIONS: The detected gene repertoire of C. resistens DSM 45100 provides insights into the lipophilic lifestyle and virulence functions of this newly recognized pathogen. Plasmid pJA144188 revealed a modular architecture of gene regions that contribute to the multi-drug resistance of C. resistens DSM 45100. The tet(W) gene encoding a ribosomal protection protein is reported here for the first time in corynebacteria. Cloning of the tet(W) gene mediated resistance to second generation tetracyclines in C. glutamicum, indicating that it might be responsible for the failure of minocycline therapies in patients with C. resistens bacteremia.

Exact and complete short-read alignment to microbial genomes using Graphics Processing Unit programming.

MOTIVATION: The introduction of next-generation sequencing techniques and especially the high-throughput systems Solexa (Illumina Inc.) and SOLiD (ABI) made the mapping of short reads to reference sequences a standard application in modern bioinformatics. Short-read alignment is needed for reference based re-sequencing of complete genomes as well as for gene expression analysis based on transcriptome sequencing. Several approaches were developed during the last years allowing for a fast alignment of short sequences to a given template. Methods available to date use heuristic techniques to gain a speedup of the alignments, thereby missing possible alignment positions. Furthermore, most approaches return only one best hit for every query sequence, thus losing the potentially valuable information of alternative alignment positions with identical scores. RESULTS: We developed SARUMAN (Semiglobal Alignment of short Reads Using CUDA and NeedleMAN-Wunsch), a mapping approach that returns all possible alignment positions of a read in a reference sequence under a given error threshold, together with one optimal alignment for each of these positions. Alignments are computed in parallel on graphics hardware, facilitating an considerable speedup of this normally time-consuming step. Combining our filter algorithm with CUDA-accelerated alignments, we were able to align reads to microbial genomes in time comparable or even faster than all published approaches, while still providing an exact, complete and optimal result. At the same time, SARUMAN runs on every standard Linux PC with a CUDA-compatible graphics accelerator. AVAILABILITY: http://www.cebitec.uni-bielefeld.de/brf/saruman/saruman.html.

Genome sequence of Wickerhamomyces anomalus DSM 6766 reveals genetic basis of biotechnologically important antimicrobial activities.

The ascomycetous yeast Wickerhamomyces anomalus (formerly Pichia anomala and Hansenula anomala) exhibits antimicrobial activities and flavoring features that are responsible for its frequent association with food, beverage and feed products. However, limited information on the genetic background of this yeast and its multiple capabilities are currently available. Here, we present the draft genome sequence of the neotype strain W. anomalus DSM 6766. On the basis of pyrosequencing, a de novo assembly of this strain resulted in a draft genome sequence with a total size of 25.47 Mbp. An automatic annotation using RAPYD generated 11 512 protein-coding sequences. This annotation provided the basis to analyse metabolic capabilities, phylogenetic relationships, as well as biotechnologically important features and yielded novel candidate genes of W. anomalus DSM 6766 coding for proteins participating in antimicrobial activities.

Genomics of the proteorhodopsin-containing marine flavobacterium Dokdonia sp. strain MED134.

Proteorhodopsin phototrophy is expected to have considerable impact on the ecology and biogeochemical roles of marine bacteria. However, the genetic features contributing to the success of proteorhodopsin-containing bacteria remain largely unknown. We investigated the genome of Dokdonia sp. strain MED134 (Bacteroidetes) for features potentially explaining its ability to grow better in light than darkness. MED134 has a relatively high number of peptidases, suggesting that amino acids are the main carbon and nitrogen sources. In addition, MED134 shares with other environmental genomes a reduction in gene copies at the expense of important ones, like membrane transporters, which might be compensated by the presence of the proteorhodopsin gene. The genome analyses suggest Dokdonia sp. MED134 is able to respond to light at least partly due to the presence of a strong flavobacterial consensus promoter sequence for the proteorhodopsin gene. Moreover, Dokdonia sp. MED134 has a complete set of anaplerotic enzymes likely to play a role in the adaptation of the carbon anabolism to the different sources of energy it can use, including light or various organic matter compounds. In addition to promoting growth, proteorhodopsin phototrophy could provide energy for the degradation of complex or recalcitrant organic matter, survival during periods of low nutrients, or uptake of amino acids and peptides at low concentrations. Our analysis suggests that the ability to harness light potentially makes MED134 less dependent on the amount and quality of organic matter or other nutrients. The genomic features reported here may well be among the keys to a successful photoheterotrophic lifestyle.

Complete genome sequence of the fire blight pathogen Erwinia pyrifoliae DSM 12163T and comparative genomic insights into plant pathogenicity.

BACKGROUND: Erwinia pyrifoliae is a newly described necrotrophic pathogen, which causes fire blight on Asian (Nashi) pear and is geographically restricted to Eastern Asia. Relatively little is known about its genetics compared to the closely related main fire blight pathogen E. amylovora. RESULTS: The genome of the type strain of E. pyrifoliae strain DSM 12163T, was sequenced using both 454 and Solexa pyrosequencing and annotated. The genome contains a circular chromosome of 4.026 Mb and four small plasmids. Based on their respective role in virulence in E. amylovora or related organisms, we identified several putative virulence factors, including type III and type VI secretion systems and their effectors, flagellar genes, sorbitol metabolism, iron uptake determinants, and quorum-sensing components. A deletion in the rpoS gene covering the most conserved region of the protein was identified which may contribute to the difference in virulence/host-range compared to E. amylovora. Comparative genomics with the pome fruit epiphyte Erwinia tasmaniensis Et1/99 showed that both species are overall highly similar, although specific differences were identified, for example the presence of some phage gene-containing regions and a high number of putative genomic islands containing transposases in the E. pyrifoliae DSM 12163T genome. CONCLUSIONS: The E. pyrifoliae genome is an important addition to the published genome of E. tasmaniensis and the unfinished genome of E. amylovora providing a foundation for re-sequencing additional strains that may shed light on the evolution of the host-range and virulence/pathogenicity of this important group of plant-associated bacteria.