Prokaryotic transcriptomics: a new view on regulation, physiology and pathogenicity

Transcriptomes of Frankia sp. strain CcI3 in growth transitions

Background

Frankia sp. strains are actinobacteria that form N₂-fixing root nodules on angiosperms. Several reference genome sequences are available enabling transcriptome studies in Frankia sp. Genomes from Frankia sp. strains differ markedly in size, a consequence proposed to be associated with a high number of indigenous transposases, more than 200 of which are found in Frankia sp. strain CcI3 used in this study. Because Frankia exhibits a high degree of cell heterogeneity as a consequence of its mycelial growth pattern, its transcriptome is likely to be quite sensitive to culture age. This study focuses on the behavior of the Frankia sp. strain CcI3 transcriptome as a function of nitrogen source and culture age.

Results

To study global transcription in Frankia sp. CcI3 grown under different conditions, complete transcriptomes were determined using high throughput RNA deep sequencing. Samples varied by time (five days vs. three days) and by culture conditions (NH₄⁺ added vs. N₂ fixing). Assembly of millions of reads revealed more diversity of gene expression between five-day and three-day old cultures than between three day old cultures differing in nitrogen sources. Heat map analysis organized genes into groups that were expressed or repressed under the various conditions compared to median expression values. Twenty-one SNPs common to all three transcriptome samples were detected indicating culture heterogeneity in this slow-growing organism. Significantly higher expression of transposase ORFs was found in the five-day and N₂-fixing cultures, suggesting that N starvation and culture aging provide conditions for on-going genome modification. Transposases have previously been proposed to participate in the creating the large number of gene duplication or deletion in host strains. Subsequent RT-qPCR experiments confirmed predicted elevated transposase expression levels indicated by the mRNA-seq data.

Conclusions

The overall pattern of gene expression in aging cultures of CcI3 suggests significant cell heterogeneity even during normal growth on ammonia. The detection of abundant transcription of nif (nitrogen fixation) genes likely reflects the presence of anaerobic, N-depleted microsites in the growing mycelium of the culture, and the presence of significantly elevated transposase transcription during starvation indicates the continuing evolution of the Frankia sp. strain CcI3 genome, even in culture, especially under stressed conditions. These studies also sound a cautionary note when comparing the transcriptomes of Frankia grown in root nodules, where cell heterogeneity would be expected to be quite high.

cis-antisense RNA, another level of gene regulation in bacteria

A substantial amount of antisense transcription is a hallmark of gene expression in eukaryotes. However, antisense transcription was first demonstrated in bacteria almost 50 years ago. The transcriptomes of bacteria as different as Helicobacter pylori, Bacillus subtilis, Escherichia coli, Synechocystis sp. strain PCC6803, Mycoplasma pneumoniae, Sinorhizobium meliloti, Geobacter sulfurreducens, Vibrio cholerae, Chlamydia trachomatis, Pseudomonas syringae, and Staphylococcus aureus have now been reported to contain antisense RNA (asRNA) transcripts for a high percentage of genes. Bacterial asRNAs share functional similarities with trans-acting regulatory RNAs, but in addition, they use their own distinct mechanisms. Among their confirmed functional roles are transcription termination, codegradation, control of translation, transcriptional interference, and enhanced stability of their respective target transcripts. Here, we review recent publications indicating that asRNAs occur as frequently in simple unicellular bacteria as they do in higher organisms, and we provide a comprehensive overview of the experimentally confirmed characteristics of asRNA actions and intimately linked quantitative aspects. Emerging functional data suggest that asRNAs in bacteria mediate a plethora of effects and are involved in far more processes than were previously anticipated. Thus, the functional impact of asRNAs should be considered when developing new strategies against pathogenic bacteria and when optimizing bacterial strains for biotechnology.

Deciphering the regulatory codes in bacterial genomes

Interactions between cis-regulatory elements and trans-acting factors are fundamental for cellular functions such as transcription. With the revolution in microarrays and sequencing technologies, genome-wide binding locations of trans-acting factors are being determined in large numbers. The richness of the genome-scale information has revealed that the nature of the bacterial transcriptome and regulome are considerably more complex than previously expected. In addition, the emerging view of the bacterial transcriptome is revising the concept of the operon organization of the genome. This review describes current advances in the genome-scale analysis of the interaction between cis-regulatory elements and trans-acting factors in microorganisms.

Bacterial transcriptomics: what is beyond the RNA horiz-ome?

Over the past 3 years, bacterial transcriptomics has undergone a massive revolution. Increased sequencing capacity and novel tools have made it possible to explore the bacterial transcriptome to an unprecedented depth, which has revealed that the transcriptome is more complex and dynamic than expected. Alternative transcripts within operons challenge the classic operon definition, and many small RNAs involved in the regulation of transcription, translation and pathogenesis have been discovered. Furthermore, mRNAs may localize to specific areas in the cell, and the spatial organization and dynamics of the chromosome have been shown to be important for transcription. Epigenetic modifications of DNA also affect transcription, and RNA processing affects translation. Therefore, transcription in bacteria resembles that in eukaryotes in terms of complexity more closely than was previously thought. Here we will discuss the contribution of 'omics' approaches to these discoveries as well as the possible impact that they are expected to have in the future.

Evidence-based annotation of transcripts and proteins in the sulfate-reducing bacterium Desulfovibrio vulgaris Hildenborough

We used high-resolution tiling microarrays and 5' RNA sequencing to identify transcripts in Desulfovibrio vulgaris Hildenborough, a model sulfate-reducing bacterium. We identified the first nucleotide position for 1,124 transcripts, including 54 proteins with leaderless transcripts and another 72 genes for which a major transcript initiates within the upstream protein-coding gene, which confounds measurements of the upstream gene's expression. Sequence analysis of these promoters showed that D. vulgaris prefers -10 and -35 boxes different from those preferred by Escherichia coli. A total of 549 transcripts ended at intrinsic (rho-independent) terminators, but most of the other transcripts seemed to have variable ends. We found low-level antisense expression of most genes, and the 5' ends of these transcripts mapped to promoter-like sequences. Because antisense expression was reduced for highly expressed genes, we suspect that elongation of nonspecific antisense transcripts is suppressed by transcription of the sense strand. Finally, we combined the transcript results with comparative analysis and proteomics data to make 505 revisions to the original annotation of 3,531 proteins: we removed 255 (7.5%) proteins, changed 123 (3.6%) start codons, and added 127 (3.7%) proteins that had been missed. Tiling data had higher coverage than shotgun proteomics and hence led to most of the corrections, but many errors probably remain. Our data are available at http://genomics.lbl.gov/supplemental/DvHtranscripts2011/.

Identification and Role of Regulatory Non-Coding RNAs in Listeria monocytogenes

Bacterial regulatory non-coding RNAs control numerous mRNA targets that direct a plethora of biological processes, such as the adaption to environmental changes, growth and virulence. Recently developed high-throughput techniques, such as genomic tiling arrays and RNA-Seq have allowed investigating prokaryotic cis- and trans-acting regulatory RNAs, including sRNAs, asRNAs, untranslated regions (UTR) and riboswitches. As a result, we obtained a more comprehensive view on the complexity and plasticity of the prokaryotic genome biology. Listeria monocytogenes was utilized as a model system for intracellular pathogenic bacteria in several studies, which revealed the presence of about 180 regulatory RNAs in the listerial genome. A regulatory role of non-coding RNAs in survival, virulence and adaptation mechanisms of L. monocytogenes was confirmed in subsequent experiments, thus, providing insight into a multifaceted modulatory function of RNA/mRNA interference. In this review, we discuss the identification of regulatory RNAs by high-throughput techniques and in their functional role in L. monocytogenes.

Characterization of 3' termini of human cytomegalovirus UL138-UL145 transcripts in a clinical strain

The functions of some proteins encoded by human cytomegalovirus (HCMV) UL/b' genes have been studied; however, systematic analysis of the transcripts for this region is still insufficient. The results of both rapid amplification of cDNA ends (RACE) and cDNA library screening in this study proved that 3' termini of all transcripts in the UL138-UL145 region were located approximately 20 bp downstream from each potential poly (A) signal, which were at the positions of nucleotides 7184, 9954 and 12848 in the UL/b' sequence of the H strain, respectively. Thus, there were at least two large families of polycistronic transcripts in this gene region. The first family of 3'-coterminal transcripts contained UL139, UL140 and UL141 genes, and the second one consisted of UL142, UL143, UL144 and UL145 genes. The 3'-coterminal characterization further confirmed that multiple uses of polyadenylation signals were commonly used by HCMV to utilize genetic information.

Next-generation sequencing and its applications in molecular diagnostics

DNA sequencing is a powerful approach for decoding a number of human diseases, including cancers. The advent of next-generation sequencing (NGS) technologies has reduced sequencing cost by orders of magnitude and significantly increased the throughput, making whole-genome sequencing a possible way for obtaining global genomic information about patients on whom clinical actions may be taken. However, the benefits offered by NGS technologies come with a number of challenges that must be adequately addressed before they can be transformed from research tools to routine clinical practices. This article provides an overview of four commonly used NGS technologies from Roche Applied Science//454 Life Sciences, Illumina, Life Technologies and Helicos Biosciences. The challenges in the analysis of NGS data and their potential applications in clinical diagnosis are also discussed.

Advances in bacterial transcriptome and transposon insertion-site profiling using second-generation sequencing

The arrival of second-generation sequencing has revolutionized the study of bacteria within a short period. The sequence information generated from these platforms has helped in our understanding of bacterial development, adaptation and diversity and how bacteria cause disease. Furthermore, these technologies have quickly been adapted for high-throughput studies that were previously performed using DNA cloning or microarray-based applications. This has facilitated a more comprehensive study of bacterial transcriptomes through RNA sequencing (RNA-Seq) and the systematic determination of gene function by 'transposon monitoring'. In this review, we provide an outline of these powerful tools and the in silico analyses used in their application, and also highlight the biological questions being addressed in these approaches.

Isolation of high-quality total RNA from rumen anaerobic bacteria and fungi, and subsequent detection of glycoside hydrolases

The rumen is one of the most powerful fibrolytic fermentation systems known. Gene expression analyses, such as reverse transcription PCR (RT-PCR), microarrays, and metatranscriptomics, are techniques that could significantly expand our understanding of this ecosystem. The ability to isolate and stabilize representative RNA samples is critical to obtaining reliable results with these procedures. In this study, we successfully isolated high-quality total RNA from the solid phase of ruminal contents by using an improved RNA extraction method. This method is based on liquid nitrogen grinding of whole ruminal solids without microbial detachment and acid guanidinium - phenol - chloroform extraction combined with column purification. Yields of total RNA were as high as 150 µg per g of fresh ruminal content. The typical large subunit/small subunit rRNA ratio ranged from 1.8 to 2.0 with an RNA integrity number (Agilent Technologies) greater than 8.5. By eliminating the detachment step, the resulting RNA was more representative of the complete ecosystem. Our improved method removed a major barrier limiting analysis of rumen microbial function from a gene expression perspective. The polyA-tailed eukaryotic mRNAs obtained have successfully been applied to next-generation sequencing, and metatranscriptomic analysis of the solid fraction of rumen contents revealed abundant sequences related to rumen fungi.

Parallel evolution of transcriptome architecture during genome reorganization

Assembly of genes into operons is generally viewed as an important process during the continual adaptation of microbes to changing environmental challenges. However, the genome reorganization events that drive this process are also the roots of instability for existing operons. We have determined that there exists a statistically significant trend that correlates the proportion of genes encoded in operons in archaea to their phylogenetic lineage. We have further characterized how microbes deal with operon instability by mapping and comparing transcriptome architectures of four phylogenetically diverse extremophiles that span the range of operon stabilities observed across archaeal lineages: a photoheterotrophic halophile (Halobacterium salinarum NRC-1), a hydrogenotrophic methanogen (Methanococcus maripaludis S2), an acidophilic and aerobic thermophile (Sulfolobus solfataricus P2), and an anaerobic hyperthermophile (Pyrococcus furiosus DSM 3638). We demonstrate how the evolution of transcriptional elements (promoters and terminators) generates new operons, restores the coordinated regulation of translocated, inverted, and newly acquired genes, and introduces completely novel regulation for even some of the most conserved operonic genes such as those encoding subunits of the ribosome. The inverse correlation (r=-0.92) between the proportion of operons with such internally located transcriptional elements and the fraction of conserved operons in each of the four archaea reveals an unprecedented view into varying stages of operon evolution. Importantly, our integrated analysis has revealed that organisms adapted to higher growth temperatures have lower tolerance for genome reorganization events that disrupt operon structures.

Directional RNA deep sequencing sheds new light on the transcriptional response of Anabaena sp. strain PCC 7120 to combined-nitrogen deprivation

Background

Cyanobacteria are potential sources of renewable chemicals and biofuels and serve as model organisms for bacterial photosynthesis, nitrogen fixation, and responses to environmental changes. Anabaena (Nostoc) sp. strain PCC 7120 (hereafter Anabaena) is a multicellular filamentous cyanobacterium that can "fix" atmospheric nitrogen into ammonia when grown in the absence of a source of combined nitrogen. Because the nitrogenase enzyme is oxygen sensitive, Anabaena forms specialized cells called heterocysts that create a microoxic environment for nitrogen fixation. We have employed directional RNA-seq to map the Anabaena transcriptome during vegetative cell growth and in response to combined-nitrogen deprivation, which induces filaments to undergo heterocyst development. Our data provide an unprecedented view of transcriptional changes in Anabaena filaments during the induction of heterocyst development and transition to diazotrophic growth.

Results

Using the Illumina short read platform and a directional RNA-seq protocol, we obtained deep sequencing data for RNA extracted from filaments at 0, 6, 12, and 21 hours after the removal of combined nitrogen. The RNA-seq data provided information on transcript abundance and boundaries for the entire transcriptome. From these data, we detected novel antisense transcripts within the UTRs (untranslated regions) and coding regions of key genes involved in heterocyst development, suggesting that antisense RNAs may be important regulators of the nitrogen response. In addition, many 5' UTRs were longer than anticipated, sometimes extending into upstream open reading frames (ORFs), and operons often showed complex structure and regulation. Finally, many genes that had not been previously identified as being involved in heterocyst development showed regulation, providing new candidates for future studies in this model organism.

Conclusions

Directional RNA-seq data were obtained that provide comprehensive mapping of transcript boundaries and abundance for all transcribed RNAs in Anabaena filaments during the response to nitrogen deprivation. We have identified genes and noncoding RNAs that are transcriptionally regulated during heterocyst development. These data provide detailed information on the Anabaena transcriptome as filaments undergo heterocyst development and begin nitrogen fixation.

Transcript amplification from single bacterium for transcriptome analysis

Total transcript amplification (TTA) from single eukaryotic cells for transcriptome analysis is established, but TTA from a single prokaryotic cell presents additional challenges with much less starting material, the lack of poly(A)-tails, and the fact that the messages can be polycistronic. Here, we describe a novel method for single-bacterium TTA using a model organism, Burkholderia thailandensis, exposed to a subinhibitory concentration of the antibacterial agent, glyphosate. Utilizing a B. thailandensis microarray to assess the TTA method showed low fold-change bias (less than twofold difference and Pearson correlation coefficient R ≈ 0.87-0.89) and drop-outs (4%-6% of 2842 detectable genes), compared with data obtained from the larger-scale nonamplified RNA samples. Further analysis of the microarray data suggests that B. thailandensis, when exposed to the aromatic amino acid biosynthesis inhibitor glyphosate, induces (or represses) genes to possibly recuperate and balance the intracellular amino acid pool. We validated our single-cell microarray data at the multi-cell and single-cell levels with lacZ and gfp reporter-gene fusions, respectively. Sanger sequencing of 192 clones generated from the TTA product of a single cell, with and without enrichment by elimination of rRNA and tRNA, detected only B. thailandensis sequences with no contamination. These data indicate that RNA-seq of TTA from a single cell is possible using this novel method.

Small RNAs as regulators of primary and secondary metabolism in Pseudomonas species

Small RNAs (sRNAs) exert important functions in pseudomonads. Classical sRNAs comprise the 4.5S, 6S, 10Sa and 10Sb RNAs, which are known in enteric bacteria as part of the signal recognition particle, a regulatory component of RNA polymerase, transfer-messenger RNA (tmRNA) and the RNA component of RNase P, respectively. Their homologues in pseudomonads are presumed to have analogous functions. Other sRNAs of pseudomonads generally have little or no sequence similarity with sRNAs of enteric bacteria. Numerous sRNAs repress or activate the translation of target mRNAs by a base-pairing mechanism. Examples of this group in Pseudomonas aeruginosa are the iron-repressible PrrF1 and PrrF2 sRNAs, which repress the translation of genes encoding iron-containing proteins, and PhrS, an anaerobically inducible sRNA, which activates the expression of PqsR, a regulator of the Pseudomonas quinolone signal. Other sRNAs sequester RNA-binding proteins that act as translational repressors. Examples of this group in P. aeruginosa include RsmY and RsmZ, which are central regulatory elements in the GacS/GacA signal transduction pathway, and CrcZ, which is a key regulator in the CbrA/CbrB signal transduction pathway. These pathways largely control the extracellular activities (including virulence traits) and the selection of the energetically most favourable carbon sources, respectively, in pseudomonads.

Deep sequencing-based identification of small non-coding RNAs in Streptomyces coelicolor

Streptomyces coelicolor is considered the model organism among Gram positive, GC rich bacteria. Its genome has been sequenced but little is known about the occurrence and distribution of small non-coding RNAs in this biotechnologically relevant organism. Using deep sequencing we analyzed the transcriptome at the end of exponential growth, which corresponds to the onset of secondary metabolism. We mapped 193 transcriptional start sites of mRNA genes and identified putative new and alternative open reading frames. We identified 63 non-coding RNAs including 29 cis encoded antisense RNAs, and confirmed expression for 11, most of them being growth-phase dependent. A comparison between the sequencing results and bioinformatic sRNA predictions using Dynalign and RNAz revealed only a small overlap between the different approaches.

Probing bacterial pathogenesis with genetics, genomics, and chemical biology: past, present, and future approaches

Classical genetic approaches for studying bacterial pathogenesis have provided a solid foundation for our current understanding of microbial physiology and the interactions between pathogen and host. During the past decade however, advances in several arenas have expanded the ways in which the biology of pathogens can be studied. This review discussed the impact of these advances on bacterial genetics, including the application of genomics and chemical biology to the study of pathogenesis.

An integrated approach for finding overlooked genes in Shigella

BACKGROUND

The completion of numerous genome sequences introduced an era of whole-genome study. However, many genes are missed during genome annotation, including small RNAs (sRNAs) and small open reading frames (sORFs). In order to improve genome annotation, we aimed to identify novel sRNAs and sORFs in Shigella, the principal etiologic agents of bacillary dysentery.

METHODOLOGY/PRINCIPAL FINDINGS

We identified 64 sRNAs in Shigella, which were experimentally validated in other bacteria based on sequence conservation. We employed computer-based and tiling array-based methods to search for sRNAs, followed by RT-PCR and northern blots, to identify nine sRNAs in Shigella flexneri strain 301 (Sf301) and 256 regions containing possible sRNA genes. We found 29 candidate sORFs using bioinformatic prediction, array hybridization and RT-PCR verification. We experimentally validated 557 (57.9%) DOOR operon predictions in the chromosomes of Sf301 and 46 (76.7%) in virulence plasmid.We found 40 additional co-expressed gene pairs that were not predicted by DOOR.

CONCLUSIONS/SIGNIFICANCE

We provide an updated and comprehensive annotation of the Shigella genome. Our study increased the expected numbers of sORFs and sRNAs, which will impact on future functional genomics and proteomics studies. Our method can be used for large scale reannotation of sRNAs and sORFs in any microbe with a known genome sequence.

The challenge of relating gene expression to the virulence of Salmonella enterica serovar Typhimurium

The first decade of transcriptomic studies of Salmonella enterica serovar Typhimurium focused upon gene expression in vitro, and during the infection of mammalian cells. The published regulons and stimulons show that the three Type Three Secretion Systems of S. Typhimurium respond to a diverse range of environmental conditions, and are controlled by a hierarchy of regulatory proteins. The integration of in vitro generated transcriptomic data with global gene expression of S. Typhimurium during infection is beginning to yield valuable information. The coordinated regulation of Salmonella gene expression is a key process for survival, adaptation and virulence capacities of the pathogen.

Metatranscriptomic approach to analyze the functional human gut microbiota

The human gut is the natural habitat for a large and dynamic bacterial community that has a great relevance for health. Metagenomics is increasing our knowledge of gene content as well as of functional and genetic variability in this microbiome. However, little is known about the active bacteria and their function(s) in the gastrointestinal tract. We performed a metatranscriptomic study on ten healthy volunteers to elucidate the active members of the gut microbiome and their functionality under conditions of health. First, the microbial cDNAs obtained from each sample were sequenced using 454 technology. The analysis of 16S transcripts showed the phylogenetic structure of the active microbial community. Lachnospiraceae, Ruminococcaceae, Bacteroidaceae, Prevotellaceae, and Rickenellaceae were the predominant families detected in the active microbiota. The characterization of mRNAs revealed a uniform functional pattern in healthy individuals. The main functional roles of the gut microbiota were carbohydrate metabolism, energy production and synthesis of cellular components. In contrast, housekeeping activities such as amino acid and lipid metabolism were underrepresented in the metatranscriptome. Our results provide new insights into the functionality of the complex gut microbiota in healthy individuals. In this RNA-based survey, we also detected small RNAs, which are important regulatory elements in prokaryotic physiology and pathogenicity.

Sequencing illustrates the transcriptional response of Legionella pneumophila during infection and identifies seventy novel small non-coding RNAs

Second generation sequencing has prompted a number of groups to re-interrogate the transcriptomes of several bacterial and archaeal species. One of the central findings has been the identification of complex networks of small non-coding RNAs that play central roles in transcriptional regulation in all growth conditions and for the pathogen's interaction with and survival within host cells. Legionella pneumophila is a Gram-negative facultative intracellular human pathogen with a distinct biphasic lifestyle. One of its primary environmental hosts in the free-living amoeba Acanthamoeba castellanii and its infection by L. pneumophila mimics that seen in human macrophages. Here we present analysis of strand specific sequencing of the transcriptional response of L. pneumophila during exponential and post-exponential broth growth and during the replicative and transmissive phase of infection inside A. castellanii. We extend previous microarray based studies as well as uncovering evidence of a complex regulatory architecture underpinned by numerous non-coding RNAs. Over seventy new non-coding RNAs could be identified; many of them appear to be strain specific and in configurations not previously reported. We discover a family of non-coding RNAs preferentially expressed during infection conditions and identify a second copy of 6S RNA in L. pneumophila. We show that the newly discovered putative 6S RNA as well as a number of other non-coding RNAs show evidence for antisense transcription. The nature and extent of the non-coding RNAs and their expression patterns suggests that these may well play central roles in the regulation of Legionella spp. specific traits and offer clues as to how L. pneumophila adapts to its intracellular niche. The expression profiles outlined in the study have been deposited into Genbank's Gene Expression Omnibus (GEO) database under the series accession GSE27232.

The intracellular sRNA transcriptome of Listeria monocytogenes during growth in macrophages

Small non-coding RNAs (sRNAs) are widespread effectors of post-transcriptional gene regulation in bacteria. Currently extensive information exists on the sRNAs of Listeria monocytogenes expressed during growth in extracellular environments. We used deep sequencing of cDNAs obtained from fractioned RNA (<500 nt) isolated from extracellularly growing bacteria and from L. monocytogenes infected macrophages to catalog the sRNA repertoire during intracellular bacterial growth. Here, we report on the discovery of 150 putative regulatory RNAs of which 71 have not been previously described. A total of 29 regulatory RNAs, including small non-coding antisense RNAs, are specifically expressed intracellularly. We validated highly expressed sRNAs by northern blotting and demonstrated by the construction and characterization of isogenic mutants of rli31, rli33-1 and rli50* for intracellular expressed sRNA candidates, that their expression is required for efficient growth of bacteria in macrophages. All three mutants were attenuated when assessed for growth in mouse and insect models of infection. Comparative genomic analysis revealed the presence of lineage specific sRNA candidates and the absence of sRNA loci in genomes of naturally occurring infection-attenuated bacteria, with additional loss in non-pathogenic listerial genomes. Our analyses reveal extensive sRNA expression as an important feature of bacterial regulation during intracellular growth.

A simple and efficient method to search for selected primary transcripts: non-coding and antisense RNAs in the human pathogen Enterococcus faecalis

Enterococcus faecalis is a commensal bacterium and a major opportunistic human pathogen. In this study, we combined in silico predictions with a novel 5′RACE-derivative method coined ‘5′tagRACE’, to perform the first search for non-coding RNAs (ncRNAs) encoded on the E. faecalis chromosome. We used the 5′tagRACE to simultaneously probe and characterize primary transcripts, and demonstrate here the simplicity, the reliability and the sensitivity of the method. The 5′tagRACE is complementary to tiling arrays or RNA-sequencing methods, and is also directly applicable to deep RNA sequencing and should significantly improve functional studies of bacterial RNA landscapes. From 45 selected loci of the E. faecalis chromosome, we discovered and mapped 29 novel ncRNAs, 10 putative novel mRNAs and 16 antisense transcriptional organizations. We describe in more detail the oxygen-dependent expression of one ncRNA located in an E. faecalis pathogenicity island, the existence of an ncRNA that is antisense to the ncRNA modulator of the RNA polymerase, SsrS and provide evidences for the functional interplay between two distinct toxin–antitoxin modules.

Application of high-throughput DNA sequencing in phytopathology

The new sequencing technologies are already making a big impact in academic research on medically important microbes and may soon revolutionize diagnostics, epidemiology, and infection control. Plant pathology also stands to gain from exploiting these opportunities. This manuscript reviews some applications of these high-throughput sequencing methods that are relevant to phytopathology, with emphasis on the associated computational and bioinformatics challenges and their solutions. Second-generation sequencing technologies have recently been exploited in genomics of both prokaryotic and eukaryotic plant pathogens. They are also proving to be useful in diagnostics, especially with respect to viruses.

Sequencing antibody repertoires: the next generation

Genomic studies have been revolutionized by the use of next generation sequencing (NGS) that delivers huge amounts of sequence information in a short span of time. The number of applications for NGS is rapidly expanding and significantly transforming many areas of life sciences. The field of antibody research and discovery is no exception. Several recent studies have harnessed the power of NGS for analyzing natural or synthetic immunoglobulin repertoires with unprecedented resolution and exploring alternative paths for antibody discovery. Thus, appreciating, and then exploiting, these advances is essential for staying at the edge of antibody innovation.

Metagenomic analyses: past and future trends

Metagenomics has revolutionized microbiology by paving the way for a cultivation-independent assessment and exploitation of microbial communities present in complex ecosystems. Metagenomics comprising construction and screening of metagenomic DNA libraries has proven to be a powerful tool to isolate new enzymes and drugs of industrial importance. So far, the majority of the metagenomically exploited habitats comprised temperate environments, such as soil and marine environments. Recently, metagenomes of extreme environments have also been used as sources of novel biocatalysts. The employment of next-generation sequencing techniques for metagenomics resulted in the generation of large sequence data sets derived from various environments, such as soil, the human body, and ocean water. Analyses of these data sets opened a window into the enormous taxonomic and functional diversity of environmental microbial communities. To assess the functional dynamics of microbial communities, metatranscriptomics and metaproteomics have been developed. The combination of DNA-based, mRNA-based, and protein-based analyses of microbial communities present in different environments is a way to elucidate the compositions, functions, and interactions of microbial communities and to link these to environmental processes.

Metagenomic analyses: past and future trends

Metagenomics has revolutionized microbiology by paving the way for a cultivation-independent assessment and exploitation of microbial communities present in complex ecosystems. Metagenomics comprising construction and screening of metagenomic DNA libraries has proven to be a powerful tool to isolate new enzymes and drugs of industrial importance. So far, the majority of the metagenomically exploited habitats comprised temperate environments, such as soil and marine environments. Recently, metagenomes of extreme environments have also been used as sources of novel biocatalysts. The employment of next-generation sequencing techniques for metagenomics resulted in the generation of large sequence data sets derived from various environments, such as soil, the human body, and ocean water. Analyses of these data sets opened a window into the enormous taxonomic and functional diversity of environmental microbial communities. To assess the functional dynamics of microbial communities, metatranscriptomics and metaproteomics have been developed. The combination of DNA-based, mRNA-based, and protein-based analyses of microbial communities present in different environments is a way to elucidate the compositions, functions, and interactions of microbial communities and to link these to environmental processes.

Structure and dynamics of the operon map of Buchnera aphidicola sp. strain APS

BACKGROUND

Gene expression regulation is still poorly documented in bacteria with highly reduced genomes. Understanding the evolution and mechanisms underlying the regulation of gene transcription in Buchnera aphidicola, the primary endosymbiont of aphids, is expected both to enhance our understanding of this nutritionally based association and to provide an intriguing case-study of the evolution of gene expression regulation in a reduced bacterial genome.

RESULTS

A Bayesian predictor was defined to infer the B. aphidicola transcription units, which were further validated using transcriptomic data and RT-PCR experiments. The characteristics of B. aphidicola predicted transcription units (TUs) were analyzed in order to evaluate the impact of operon map organization on the regulation of gene transcription.On average, B. aphidicola TUs contain more genes than those of E. coli. The global layout of B. aphidicola operon map was mainly shaped by the big reduction and the rearrangements events, which occurred at the early stage of the symbiosis. Our analysis suggests that this operon map may evolve further only by small reorganizations around the frontiers of B. aphidicola TUs, through promoter and/or terminator sequence modifications and/or by pseudogenization events. We also found that the need for specific transcription regulation exerts some pressure on gene conservation, but not on gene assembling in the operon map in Buchnera. Our analysis of the TUs spacing pointed out that a selection pressure is maintained on the length of the intergenic regions between divergent adjacent gene pairs.

CONCLUSIONS

B. aphidicola can seemingly only evolve towards a more polycistronic operon map. This implies that gene transcription regulation is probably subject to weak selection pressure in Buchnera conserving operons composed of genes with unrelated functions.

The dawning era of comprehensive transcriptome analysis in cellular microbiology

Bacteria rapidly change their transcriptional patterns during infection in order to adapt to the host environment. To investigate host–bacteria interactions, various strategies including the use of animal infection models, in vitro assay systems and microscopic observations have been used. However, these studies primarily focused on a few specific genes and molecules in bacteria. High-density tiling arrays and massively parallel sequencing analyses are rapidly improving our understanding of the complex host–bacterial interactions through identification and characterization of bacterial transcriptomes. Information resulting from these high-throughput techniques will continue to provide novel information on the complexity, plasticity, and regulation of bacterial transcriptomes as well as their adaptive responses relative to pathogenecity. Here we summarize recent studies using these new technologies and discuss the utility of transcriptome analysis.

ODB: a database for operon organizations, 2011 update

ODB (Operon DataBase) aims to collect data of all known and conserved operons in completely sequenced genomes. Three newly updated features of this database have been added as follows: (i) Data from included operons were updated. The genome-wide analysis of transcription and transcriptional units has become popular recently and ODB successfully integrates these high-throughput operon data, including genome-wide transcriptional units of five prokaryotes and two eukaryotes. The current version of our database contains information from about 10 000 known operons in more than 50 genomes, and more than 400 000 conserved operons obtained from more than 1000 bacterial genomes. (ii) ODB proposes the idea of reference operons as a new operon prediction tool. A reference operon, a set of possible orthologous genes that organize operons, is defined by clustering all known operons. A large number of known operons, including the recently added genome-wide analysis of operons, allowed us to define more reliable reference operons. (iii) ODB also provides new graphical interfaces. One is for comparative analyses of operon structures in multiple genomes. The other is for visualization of possible operons in multiple genomes obtained from the reference operons. The 2011 updated version of ODB is now available at http://operondb.jp/.

Regulatory RNA in bacterial pathogens

Bacteria constitute a large and diverse class of infectious agents, causing devastating diseases in humans, animals, and plants. Our understanding of gene expression control, which forms the basis for successful prevention and treatment strategies, has until recently neglected the many roles that regulatory RNAs might have in bacteria. In recent years, several such regulators have been found to facilitate host-microbe interactions and act as key switches between saprophytic and pathogenic lifestyles. This review covers the versatile regulatory RNA mechanisms employed by bacterial pathogens and highlights the dynamic interplay between riboregulation and virulence factor expression.

Metagenomic era for biocatalyst identification

Microbial enzymes have many known applications as biocatalysts. However, only a few of them are currently employed for biocatalysis even though an annotated collection of more than 190 billion bases is available in metagenome sequence databases from uncultured and highly diverse microbial populations. This review aims at providing conceptual and technical bases for the translation of metagenome data into both experimental and computational frameworks that facilitates a comprehensive analysis of the biocatalysts diversity space. We will also briefly present the status of the current capabilities that assess and predict catalytic potential of environmental sites and track its diversity and evolution in large-scale biocatalysis process resulting from studies applying metagenomics in association with gene fingerprinting, catabolic arrays and complementary '-omics'.

High-throughput sequencing and clinical microbiology: progress, opportunities and challenges

High-throughput sequencing is sweeping through clinical microbiology, transforming our discipline in its wake. It is already providing an enhanced view of pathogen biology through rapid and inexpensive whole-genome sequencing and more sophisticated applications such as RNA-seq. It also promises to deliver high-resolution genomic epidemiology as the ultimate typing method for bacteria. However, the most revolutionary effect of this 'disruptive technology' is likely to be creation of a novel sequence-based, culture-independent diagnostic microbiology that incorporates microbial community profiling, metagenomics and single-cell genomics. We should prepare for the coming 'technological singularity' in sequencing, when this technology becomes so fast and so cheap that it threatens to out-compete existing diagnostic and typing methods in microbiology.

Deciphering the physiological blueprint of a bacterial cell: revelations of unanticipated complexity in transcriptome and proteome

During the last few months, several pioneer genome-wide transcriptomic, proteomic and metabolomic studies have revolutionised the understanding of bacterial biological processes, leading to a picture that resembles eukaryotic complexity. Technological advances such as next-generation high-throughput sequencing and high-density oligonucleotide microarrays have allowed the determination, in several bacteria, of the entire boundaries of all expressed transcripts. Consequently, novel RNA-mediated regulatory mechanisms have been discovered including multifunctional RNAs. Moreover, resolution of bacterial proteome organisation (interactome) and global protein localisation (localizome) have unveiled an unanticipated complexity that highlights the significance of protein multifunctionality and localisation in the cell. Also, analysis of a complete bacterial metabolic network has again revealed a high fraction of multifunctional enzymes and an unexpectedly high level of metabolic responses and adaptation. Altogether, these novel approaches have permitted the deciphering of the entire physiological landscape of one of the smallest bacteria, Mycoplasma pneumoniae. Here, we summarise and discuss recent findings aimed at defining the blueprint of any prokaryote.

Advantages and limitations of current network inference methods

Network inference, which is the reconstruction of biological networks from high-throughput data, can provide valuable information about the regulation of gene expression in cells. However, it is an underdetermined problem, as the number of interactions that can be inferred exceeds the number of independent measurements. Different state-of-the-art tools for network inference use specific assumptions and simplifications to deal with underdetermination, and these influence the inferences. The outcome of network inference therefore varies between tools and can be highly complementary. Here we categorize the available tools according to the strategies that they use to deal with the problem of underdetermination. Such categorization allows an insight into why a certain tool is more appropriate for the specific research question or data set at hand.

Archaeal/eukaryal RNase P: subunits, functions and RNA diversification

RNase P, a catalytic ribonucleoprotein (RNP), is best known for its role in precursor tRNA processing. Recent discoveries have revealed that eukaryal RNase P is also required for transcription and processing of select non-coding RNAs, thus enmeshing RNase P in an intricate network of machineries required for gene expression. Moreover, the RNase P RNA seems to have been subject to gene duplication, selection and divergence to generate two new catalytic RNPs, RNase MRP and MRP-TERT, which perform novel functions encompassing cell cycle control and stem cell biology. We present new evidence and perspectives on the functional diversification of the RNase P RNA to highlight it as a paradigm for the evolutionary plasticity that underlies the extant broad repertoire of catalytic and unexpected regulatory roles played by RNA-driven RNPs.

Unique cost dynamics elucidate the role of frameshifting errors in promoting translational robustness

There is now considerable evidence supporting the view that codon usage is frequently under selection for translational accuracy. There are, however, multiple forms of inaccuracy (missense, premature termination, and frameshifting errors) and pinpointing a particular error process behind apparently adaptive mRNA anatomy is rarely straightforward. Understanding differences in the fitness costs associated with different types of translational error can help us devise critical tests that can implicate one error process to the exclusion of others. To this end, we present a model that captures distinct features of frameshifting cost and apply this to 641 prokaryotic genomes. We demonstrate that, although it is commonly assumed that the ribosome encounters an off-frame stop codon soon after the frameshift and costs of mis-elongation are therefore limited, genomes with high GC content typically incur much larger per-error costs. We go on to derive the prediction, unique to frameshifting errors, that differences in translational robustness between the 5' and 3' ends of genes should be less pronounced in genomes with higher GC content. This prediction we show to be correct. Surprisingly, this does not mean that GC-rich organisms necessarily carry a greater fitness burden as a consequence of accidental frameshifting. Indeed, increased per-error costs are often more than counterbalanced by lower predicted error rates owing to more diverse anticodon repertoires in GC-rich genomes. We therefore propose that selection on tRNA repertoires may operate to reduce frameshifting errors.

Structural and operational complexity of the Geobacter sulfurreducens genome

Prokaryotic genomes can be annotated based on their structural, operational, and functional properties. These annotations provide the pivotal scaffold for understanding cellular functions on a genome-scale, such as metabolism and transcriptional regulation. Here, we describe a systems approach to simultaneously determine the structural and operational annotation of the Geobacter sulfurreducens genome. Integration of proteomics, transcriptomics, RNA polymerase, and sigma factor-binding information with deep-sequencing-based analysis of primary 5'-end transcripts allowed for a most precise annotation. The structural annotation is comprised of numerous previously undetected genes, noncoding RNAs, prevalent leaderless mRNA transcripts, and antisense transcripts. When compared with other prokaryotes, we found that the number of antisense transcripts reversely correlated with genome size. The operational annotation consists of 1453 operons, 22% of which have multiple transcription start sites that use different RNA polymerase holoenzymes. Several operons with multiple transcription start sites encoded genes with essential functions, giving insight into the regulatory complexity of the genome. The experimentally determined structural and operational annotations can be combined with functional annotation, yielding a new three-level annotation that greatly expands our understanding of prokaryotic genomes.

Understanding PRRSV infection in porcine lung based on genome-wide transcriptome response identified by deep sequencing

Porcine reproductive and respiratory syndrome (PRRS) has been one of the most economically important diseases affecting swine industry worldwide and causes great economic losses each year. PRRS virus (PRRSV) replicates mainly in porcine alveolar macrophages (PAMs) and dendritic cells (DCs) and develops persistent infections, antibody-dependent enhancement (ADE), interstitial pneumonia and immunosuppression. But the molecular mechanisms of PRRSV infection still are poorly understood. Here we report on the first genome-wide host transcriptional responses to classical North American type PRRSV (N-PRRSV) strain CH 1a infection using Solexa/Illumina's digital gene expression (DGE) system, a tag-based high-throughput transcriptome sequencing method, and analyse systematically the relationship between pulmonary gene expression profiles after N-PRRSV infection and infection pathology. Our results suggest that N-PRRSV appeared to utilize multiple strategies for its replication and spread in infected pigs, including subverting host innate immune response, inducing an anti-apoptotic and anti-inflammatory state as well as developing ADE. Upregulation expression of virus-induced pro-inflammatory cytokines, chemokines, adhesion molecules and inflammatory enzymes and inflammatory cells, antibodies, complement activation were likely to result in the development of inflammatory responses during N-PRRSV infection processes. N-PRRSV-induced immunosuppression might be mediated by apoptosis of infected cells, which caused depletion of immune cells and induced an anti-inflammatory cytokine response in which they were unable to eradicate the primary infection. Our systems analysis will benefit for better understanding the molecular pathogenesis of N-PRRSV infection, developing novel antiviral therapies and identifying genetic components for swine resistance/susceptibility to PRRS.

Identification of novel non-coding small RNAs from Streptococcus pneumoniae TIGR4 using high-resolution genome tiling arrays

Background

The identification of non-coding transcripts in human, mouse, and Escherichia coli has revealed their widespread occurrence and functional importance in both eukaryotic and prokaryotic life. In prokaryotes, studies have shown that non-coding transcripts participate in a broad range of cellular functions like gene regulation, stress and virulence. However, very little is known about non-coding transcripts in Streptococcus pneumoniae (pneumococcus), an obligate human respiratory pathogen responsible for significant worldwide morbidity and mortality. Tiling microarrays enable genome wide mRNA profiling as well as identification of novel transcripts at a high-resolution.

Results

Here, we describe a high-resolution transcription map of the S. pneumoniae clinical isolate TIGR4 using genomic tiling arrays. Our results indicate that approximately 66% of the genome is expressed under our experimental conditions. We identified a total of 50 non-coding small RNAs (sRNAs) from the intergenic regions, of which 36 had no predicted function. Half of the identified sRNA sequences were found to be unique to S. pneumoniae genome. We identified eight overrepresented sequence motifs among sRNA sequences that correspond to sRNAs in different functional categories. Tiling arrays also identified approximately 202 operon structures in the genome.

Conclusions

In summary, the pneumococcal operon structures and novel sRNAs identified in this study enhance our understanding of the complexity and extent of the pneumococcal 'expressed' genome. Furthermore, the results of this study open up new avenues of research for understanding the complex RNA regulatory network governing S. pneumoniae physiology and virulence.

Recent progress and new challenges in metagenomics for biotechnology

A brief historical perspective on metagenomics is given followed by a discussion of the rapid progress in this field largely defined by transition to the next generation sequencing technologies. Problems and challenges connected to this transition are also addressed. The review focuses on recent literature describing metagenomic approaches connecting sequence information to functionality that are especially relevant to biotechnological applications, including metagenomics of specialized or enriched microbial communities, metagenomics combined with specific labeling techniques, metatranscriptomics and metaproteomics.

Transcriptome analysis of Pseudomonas syringae identifies new genes, noncoding RNAs, and antisense activity

To fully understand how bacteria respond to their environment, it is essential to assess genome-wide transcriptional activity. New high-throughput sequencing technologies make it possible to query the transcriptome of an organism in an efficient unbiased manner. We applied a strand-specific method to sequence bacterial transcripts using Illumina's high-throughput sequencing technology. The resulting sequences were used to construct genome-wide transcriptional profiles. Novel bioinformatics analyses were developed and used in combination with proteomics data for the qualitative classification of transcriptional activity in defined regions. As expected, most transcriptional activity was consistent with predictions from the genome annotation. Importantly, we identified and confirmed transcriptional activity in areas of the genome inconsistent with the annotation and in unannotated regions. Further analyses revealed potential RpoN-dependent promoter sequences upstream of several noncoding RNAs (ncRNAs), suggesting a role for these ncRNAs in RpoN-dependent phenotypes. We were also able to validate a number of transcriptional start sites, many of which were consistent with predicted promoter motifs. Overall, our approach provides an efficient way to survey global transcriptional activity in bacteria and enables rapid discovery of specific areas in the genome that merit further investigation.

The primary transcriptome of the major human pathogen Helicobacter pylori

Genome sequencing of Helicobacter pylori has revealed the potential proteins and genetic diversity of this prevalent human pathogen, yet little is known about its transcriptional organization and noncoding RNA output. Massively parallel cDNA sequencing (RNA-seq) has been revolutionizing global transcriptomic analysis. Here, using a novel differential approach (dRNA-seq) selective for the 5' end of primary transcripts, we present a genome-wide map of H. pylori transcriptional start sites and operons. We discovered hundreds of transcriptional start sites within operons, and opposite to annotated genes, indicating that complexity of gene expression from the small H. pylori genome is increased by uncoupling of polycistrons and by genome-wide antisense transcription. We also discovered an unexpected number of approximately 60 small RNAs including the epsilon-subdivision counterpart of the regulatory 6S RNA and associated RNA products, and potential regulators of cis- and trans-encoded target messenger RNAs. Our approach establishes a paradigm for mapping and annotating the primary transcriptomes of many living species.