Application of RNA-seq to reveal the transcript profile in bacteria

De novo assembly and annotation of the Amblyomma hebraeum tick midgut transcriptome response to Ehrlichia ruminantium infection

The South African bont tick Amblyomma hebraeum is a hematophagous vector for the heartwater disease pathogen Ehrlichia ruminantium in southern Africa. During feeding, the tick's enterocytes express proteins that perform vital functions in blood digestion, including proteins that may be involved in E. ruminantium acquisition, colonization or immunity. To delineate the molecular mechanism of midgut response to E. ruminantium infection, we performed comparative analyses of midgut transcriptomes of E. ruminantium infected engorged A. hebraeum nymphs, and infected adult male and female ticks with their corresponding matched uninfected controls, before and during feeding. A total of 102,036 unigenes were annotated in public databases and their expression levels analyzed for engorged nymphs as well as unfed and partly-fed adult ticks. There were 2,025 differentially expressed genes (DEGs) in midguts, of which 1,225 unigenes were up-regulated and 800 unigenes were down-regulated in the midguts of infected ticks. Annotation of DEGs revealed an increase in metabolic and cellular processes among E. ruminantium infected ticks. Notably, among the infected ticks, there was up-regulation in the expression of genes involved in tick immunity, histone proteins and oxidative stress responses. We also observed up-regulation of glycoproteins that E. ruminantium could potentially use as docking sites for host cell entry. Insights uncovered in this study offer a platform for further investigations into the molecular interaction between E. ruminantium and A. hebraeum.

Transcript Profiling of Nitroxoline-Treated Biofilms Shows Rapid Up-regulation of Iron Acquisition Gene Clusters

Bacterial biofilms are surface-attached communities of slow- or non-replicating cells embedded within a protective matrix of biomolecules. Unlike free-floating planktonic bacteria, biofilms are innately tolerant to conventional antibiotics and are prevalent in recurring and chronic infections. Nitroxoline, a broad-spectrum biofilm-eradicating agent, was used to probe biofilm viability. Transcript profiling (RNA-seq) showed that 452 of 2594 genes (17.4%) in methicillin-resistant Staphylococcus aureus (MRSA) biofilms were differentially expressed after a 2 h treatment of nitroxoline. WoPPER analysis and time-course validation (RT-qPCR) revealed that gene clusters involved in iron acquisition (sbn, isd, MW2101, MW0695, fhu, and feo) were rapidly up-regulated following nitroxoline treatment, which is indicative of iron starvation in MRSA biofilms. In addition, genes related to oligopeptide transporters and riboflavin biosynthesis were found to be up-regulated, while genes related to carotenoid biosynthesis and nitrate assimilation were down-regulated. RT-qPCR experiments revealed that iron uptake transcripts were also up-regulated in established Staphylococcus epidermidis and Acinetobacter baumannii biofilms following nitroxoline treatment. Overall, we show RNA-seq to be an ideal platform to define cellular pathways critical for biofilm survival, in addition to demonstrating the need these bacterial communities have for iron.

Current State and Challenges of the Global Outcomes of Dental Caries Research in the Meta-Omics Era

Despite significant healthcare advances in the 21st century, the exact etiology of dental caries remains unsolved. The past two decades have witnessed a tremendous growth in our understanding of dental caries amid the advent of revolutionary omics technologies. Accordingly, a consensus has been reached that dental caries is a community-scale metabolic disorder, and its etiology is beyond a single causative organism. This conclusion was based on a variety of microbiome studies following the flow of information along the central dogma of biology from genomic data to the end products of metabolism. These studies were facilitated by the unprecedented growth of the next- generation sequencing tools and omics techniques, such as metagenomics and metatranscriptomics, to estimate the community composition of oral microbiome and its functional potential. Furthermore, the rapidly evolving proteomics and metabolomics platforms, including nuclear magnetic resonance spectroscopy and/or mass spectrometry coupled with chromatography, have enabled precise quantification of the translational outcomes. Although the majority supports ‘conserved functional changes’ as indicators of dysbiosis, it remains unclear how caries dynamics impact the microbiota functions and vice versa, over the course of disease onset and progression. What compounds the situation is the host-microbiota crosstalk. Genome-wide association studies have been undertaken to elucidate the interaction of host genetic variation with the microbiome. However, these studies are challenged by the complex interaction of host genetics and environmental factors. All these complementary approaches need to be orchestrated to capture the key players in this multifactorial disease. Herein, we critically review the milestones in caries research focusing on the state-of-art singular and integrative omics studies, supplemented with a bibliographic network analysis to address the oral microbiome, the host factors, and their interactions. Additionally, we highlight gaps in the dental literature and shed light on critical future research questions and study designs that could unravel the complexities of dental caries, the most globally widespread disease.

RNA-Seq-Based Transcriptomic Analysis of Saccharopolyspora spinosa Revealed the Critical Function of PEP Phosphonomutase in the Replenishment Pathway

Spinosyns, the secondary metabolites produced by Saccharopolyspora spinosa, are the active ingredients in a family of novel biological insecticides. Although the complete genome sequence of S. spinosa has been published, the transcriptome of S. spinosa remains poorly characterized. In this study, high-throughput RNA sequencing (RNA-seq) technology was applied to dissect the transcriptome of S. spinosa. Through transcriptomic analysis of different periods of S. spinosa growth, we found large numbers of differentially expressed genes and classified them according to their different functions. Based on the RNA-seq data, the CRISPR-Cas9 method was used to knock out the PEP phosphonomutase gene (orf 06952-4171). The yield of spinosyns A and D in S. spinosa-ΔPEP was 178.91 mg/L and 42.72 mg/L, which was 2.14-fold and 1.76-fold higher than that in the wild type (83.51 and 24.34 mg/L), respectively. The analysis of the mutant strains also verified the validity of the transcriptome data. The deletion of the PEP phosphonomutase gene leads to an increase in pyruvate content and affects the biosynthesis of spinosad. The replenishment of phosphoenol pyruvate in S. spinosa provides the substrate for the production of spinosad. We envision that these transcriptomic analysis results will contribute to the further study of secondary metabolites in actinomycetes.

Subinhibitory concentrations of resveratrol reduce alpha-hemolysin production in Staphylococcus aureus isolates by downregulating saeRS

Resveratrol is a natural phytoalexin. In recent studies, it has been shown to have beneficial effects on cardiovascular disease and cancer and has been deemed to have effective antiviral and immunomodulatory activities. Methicillin-resistant Staphylococcus aureus is a multidrug-resistant pathogen associated with skin and soft tissue infections. Alpha-hemolysin is known to play a key role in the symptoms caused by S. aureus, and the saeRS two-component system has been shown to be a major regulatory system of S. aureus virulence. The present study was designed to determine the effect of subinhibitory concentrations of resveratrol on the production of alpha-hemolysin in S. aureus. The effect of resveratrol on the transcription of S. aureus was studied by transcriptome sequencing. A total of 760 genes with >2-fold changes in expression were selected, including 479 upregulated genes and 281 downregulated genes. On the basis of transcriptome sequencing, the expression of alpha-hemolysin in the S. aureus strains of the resveratrol-treated group was downregulated. Our results showed that resveratrol weakly inhibited the growth of S. aureus strains, and subinhibitory concentration of resveratrol decreased the expression of hla and inhibited the regulation of saeRS. Hemolysis testing confirmed that resveratrol had an inhibitory effect on the hemolysis of rabbit erythrocytes infected with S. aureus strains in a dose-dependent manner. Resveratrol also decreased the hemolytic capacity by reducing the production of alpha-hemolysin. We found that resveratrol could decrease the expression of hla and reduce the secretion of alpha-hemolysin by downregulating saeRS. These findings have provided more evidence of the potential of resveratrol as a drug for resisting S. aureus infections.

Tackling maize fusariosis: in search of Fusarium graminearum biosuppressors

This review presents biocontrol agents employed to alleviate the deleterious effect of the pathogen Fusarium graminearum on maize. The control of this mycotoxigenic phytopathogen remains elusive despite the elaborate research conducted on its detection, identification, and molecular fingerprinting. This could be attributed to the fact that in vitro and greenhouse biocontrol studies on F. graminearum have exceeded the number of field studies done. Furthermore, along with the variances seen among these F. graminearum suppressing biocontrol strains, it is also clear that the majority of research done to tackle F. graminearum outbreaks was on wheat and barley cultivars. Most fusariosis management related to maize targeted other members of Fusarium such as Fusarium verticillioides, with biocontrol strains from the genera Bacillus and Pseudomonas being used frequently in the experiments. We highlight relevant current techniques needed to identify an effective biofungicide for maize fusariosis and recommend alternative approaches to reduce the scarcity of data for indigenous maize field trials.

The oral bacterial microbiome of occlusal surfaces in children and its association with diet and caries

Dental caries is the most prevalent disease in humans globally. Efforts to control it have been invigorated by an increasing knowledge of the oral microbiome composition. This study aimed to evaluate the bacterial diversity in occlusal biofilms and its relationship with clinical surface diagnosis and dietary habits. Anamneses were recorded from thirteen 12-year-old children. Biofilm samples collected from occlusal surfaces of 46 permanent second molars were analyzed by 16S rRNA amplicon sequencing combined with the BLASTN-based search algorithm for species identification. The overall mean decayed, missing and filled surfaces modified index [DMFSm Index, including active white spot lesions (AWSL)] value was 8.77±7.47. Biofilm communities were highly polymicrobial collectively, representing 10 bacterial phyla, 25 classes, 29 orders, 58 families, 107 genera, 723 species. Streptococcus sp_Oral_Taxon_065, Corynebacterium matruchotii, Actinomyces viscosus, Actinomyces sp_Oral_Taxon_175, Actinomyces sp_Oral_Taxon_178, Actinomyces sp_Oral_Taxon_877, Prevotella nigrescens, Dialister micraerophilus, Eubacterium_XI G 1 infirmum were more abundant among surfaces with AWSL, and Streptococcus gordonii, Streptococcus sp._Oral_Taxon_058, Enterobacter sp._str._638 Streptococcus australis, Yersinia mollaretii, Enterobacter cloacae, Streptococcus sp._Oral_Taxon_71, Streptococcus sp._Oral_Taxon_F11, Centipeda sp._Oral_Taxon_D18 were more abundant among sound surfaces. Streptococcus mutans was detected on all surfaces in all patients, while Streptococcus sobrinus was detected only in three patients (mean relative abundances 7.1% and 0.6%, respectively). Neither species differentiated healthy from diseased sites. Diets of nine of the subjects were scored as high in fermentable carbohydrates (≧2X/day between meals). A direct association between relative abundances of bacteria and carbohydrate consumption was observed among 18 species. High consumption of fermentable carbohydrates and sound surfaces were associated with a reduction in bacterial diversity. PCoA plots displayed differences in bacterial community profiles between sound and diseased surfaces. Our study showed that, in addition to mutans streptococci, other species may be associated with the initiation of dental caries on occlusal surfaces, and that biofilm diversity of tooth surfaces is influenced by carbohydrate consumption and a surface's health status.

New Insights into the Anti-pathogenic Potential of Lactococcus garvieae against Staphylococcus aureus Based on RNA Sequencing Profiling

The bio-preservation potential of Lactococcus garvieae lies in its capacity to inhibit the growth of staphylococci, especially Staphylococcus aureus, in dairy products and in vitro. In vitro, inhibition is modulated by the level of aeration, owing to hydrogen peroxide (H₂O₂) production by L. garvieae under aeration. The S. aureus response to this inhibition has already been studied. However, the molecular mechanisms of L. garvieae underlying the antagonism against S. aureus have never been explored. This study provides evidence of the presence of another extracellular inhibition effector in vitro. This effector was neither a protein, nor a lipid, nor a polysaccharide, nor related to an L-threonine deficiency. To better understand the H₂O₂-related inhibition mechanism at the transcriptome level and to identify other mechanisms potentially involved, we used RNA sequencing to determine the transcriptome response of L. garvieae to different aeration levels and to the presence or absence of S. aureus. The L. garvieae transcriptome differed radically between different aeration levels mainly in biological processes related to fundamental functions and nutritional adaptation. The transcriptomic response of L. garvieae to aeration level differed according to the presence or absence of S. aureus. The higher concentration of H₂O₂ with high aeration was not associated with a higher expression of L. garvieae H₂O₂-synthesis genes (pox, sodA, and spxA1) but rather with a repression of L. garvieae H₂O₂-degradation genes (trxB1, ahpC, ahpF, and gpx). We showed that L. garvieae displayed an original, previously undiscovered, H₂O₂ production regulation mechanism among bacteria. In addition to the key factor H₂O₂, the involvement of another extracellular effector in the antagonism against S. aureus was shown. Future studies should explore the relation between H₂O₂-metabolism, H₂O₂-producing LAB and the pathogen they inhibit. The nature of the other extracellular effector should also be determined.

Ribosomal mutations promote the evolution of antibiotic resistance in a multidrug environment

Antibiotic resistance arising via chromosomal mutations is typically specific to a particular antibiotic or class of antibiotics. We have identified mutations in genes encoding ribosomal components in Mycobacterium smegmatis that confer resistance to several structurally and mechanistically unrelated classes of antibiotics and enhance survival following heat shock and membrane stress. These mutations affect ribosome assembly and cause large-scale transcriptomic and proteomic changes, including the downregulation of the catalase KatG, an activating enzyme required for isoniazid sensitivity, and upregulation of WhiB7, a transcription factor involved in innate antibiotic resistance. Importantly, while these ribosomal mutations have a fitness cost in antibiotic-free medium, in a multidrug environment they promote the evolution of high-level, target-based resistance. Further, suppressor mutations can then be easily acquired to restore wild-type growth. Thus, ribosomal mutations can serve as stepping-stones in an evolutionary path leading to the emergence of high-level, multidrug resistance.

DOI:http://dx.doi.org/10.7554/eLife.20420.001

The rise of antibiotic resistant bacteria is challenging clinicians, and some infections are now resistant to almost all of the drugs that are currently available. Some types of bacteria – such as mycobacteria, which include the bacteria that cause tuberculosis and leprosy – can only acquire antibiotic resistance from mutations that alter their existing genes. The process by which bacteria develop resistance to multiple drugs is generally viewed as a stepwise accumulation of different mutations. However, the role of individual mutations that increase a bacterium’s resistance to multiple antibiotics has not been fully explored.

Gomez, Kaufmann-Malaga et al. exposed bacteria from the species Mycobacterium smegmatis, a cousin of the bacterium that causes tuberculosis, to a mixture of relatively low concentrations of different antibiotics that should kill the bacteria relatively slowly. Hundreds of small bacteria cultures were grown in parallel, and only a fraction of them developed antibiotic-resistant members. Gomez, Kaufmann-Malaga et al. identified mutations in these bacteria that unexpectedly gave the bacteria resistance to several unrelated classes of antibiotics.

Individual mutants carried single mutations in different components of the ribosome, a complex molecular machine that helps to build proteins inside cells. As well as increasing their resistance to antibiotics, these mutations also reduced the growth rate of the bacteria. This meant that when the bacteria were grown in an antibiotic-free environment they survived less well than non-mutant bacteria. However, the mutations gave the bacteria an advantage in environments that contained many different antibiotics, as they could more easily develop mutations that made them more resistant to other drugs. Thus, the mutant bacteria can serve as stepping-stones toward the development of high-level resistance to multiple drugs.

Further work will now explore whether this phenomenon occurs in a range of other bacterial species, including the bacteria that cause tuberculosis. While new antibiotics are desperately needed, a better understanding of how bacteria evolve the ability to resist the effects of antibiotics will help us to preserve the usefulness of existing and future drugs.

DOI:http://dx.doi.org/10.7554/eLife.20420.002

RETRACTED ARTICLE: Comparative analysis of the Corynebacterium glutamicum transcriptome in response to changes in dissolved oxygen levels

The dissolved oxygen (DO) level of a culture of Corynebacterium glutamicum (C. glutamicum) in a bioreactor has a significant impact on the cellular redox potential and the distribution of energy and metabolites. In this study, to gain a deeper understanding of the effects of DO on the metabolism of C. glutamicum, we sought to systematically explore the influence of different DO concentrations on genetic regulation and metabolism through transcriptomic analysis. The results revealed that after 20 h of fermentation, oxygen limitation enhanced the glucose metabolism, pyruvate metabolism and carbon overflow, and restricted NAD+ availability. A high oxygen supply enhanced the TCA cycle and reduced glyoxylate metabolism. Several key genes involved in response of C. glutamicum to different oxygen concentrations were examined, which provided suggestions for target site modifications in developing optimized oxygen supply strategies. These data provided new insights into the relationship between oxygen supply and metabolism of C. glutamicum.

Transcriptome Analysis of the Intracellular Facultative Pathogen Piscirickettsia salmonis: Expression of Putative Groups of Genes Associated with Virulence and Iron Metabolism

The intracellular facultative bacteria Piscirickettsia salmonis is one of the most important pathogens of the Chilean aquaculture. However, there is a lack of information regarding the whole genomic transcriptional response according to different extracellular environments. We used next generation sequencing (NGS) of RNA (RNA-seq) to study the whole transcriptome of an isolate of P. salmonis (FAVET-INBIOGEN) using a cell line culture and a modified cell-free liquid medium, with or without iron supplementation. This was done in order to obtain information about the factors there are involved in virulence and iron acquisition. First, the isolate was grown in the Sf21 cell line; then, the bacteria were cultured into a cell-free liquid medium supplemented or not with iron. We identified in the transcriptome, genes associated with type IV secretion systems, genes related to flagellar structure assembly, several proteases and sigma factors, and genes related to the development of drug resistance. Additionally, we identified for the first time several iron-metabolism associated genes including at least two iron uptake pathways (ferrous iron and ferric iron uptake) that are actually expressed in the different conditions analyzed. We further describe putative genes that are related with the use and storage of iron in the bacteria, which have not been previously described. Several sets of genes related to virulence were expressed in both the cell line and cell-free culture media (for example those related to flagellar structure; such as basal body, MS-ring, C-ring, proximal and distal rod, and filament), which may play roles in other basic processes rather than been restricted to virulence.

Replicates, Read Numbers, and Other Important Experimental Design Considerations for Microbial RNA-seq Identified Using Bacillus thuringiensis Datasets

RNA-seq is being used increasingly for gene expression studies and it is revolutionizing the fields of genomics and transcriptomics. However, the field of RNA-seq analysis is still evolving. Therefore, we specifically designed this study to contain large numbers of reads and four biological replicates per condition so we could alter these parameters and assess their impact on differential expression results. Bacillus thuringiensis strains ATCC10792 and CT43 were grown in two Luria broth medium lots on four dates and transcriptomics data were generated using one lane of sequence output from an Illumina HiSeq2000 instrument for each of the 32 samples, which were then analyzed using DESeq2. Genome coverages across samples ranged from 87 to 465X with medium lots and culture dates identified as major variation sources. Significantly differentially expressed genes (5% FDR, two-fold change) were detected for cultures grown using different medium lots and between different dates. The highly differentially expressed iron acquisition and metabolism genes, were a likely consequence of differing amounts of iron in the two media lots. Indeed, in this study RNA-seq was a tool for predictive biology since we hypothesized and confirmed the two LB medium lots had different iron contents (~two-fold difference). This study shows that the noise in data can be controlled and minimized with appropriate experimental design and by having the appropriate number of replicates and reads for the system being studied. We outline parameters for an efficient and cost effective microbial transcriptomics study.

RNA-seq reveals the critical role of CspA in regulating Brucella melitensis metabolism and virulence

Brucella melitensis is a facultative intracellular bacterium that replicates within macrophages. The ability of Brucella to survive and multiply in the hostile environment of host macrophages is essential for its virulence. The cold shock protein CspA plays an important role in the virulence of B. melitensis. To analyze the genes regulated by CspA, the whole transcriptomes of B. melitensis NIΔcspA and its parental wild-type strain, B. melitensis NI, were sequenced and analyzed using the Solexa/Illumina sequencing platform. A total of 446 differentially expressed genes were identified, including 324 up-regulated and 122 down-regulated genes. Numerous genes identified are involved in amino acid, fatty acid, nitrogen, and energy metabolism. Interestingly, all genes involved in the type IV secretion system and LuxR-type regulatory protein VjbR were significantly down-regulated in NIΔcspA. In addition, an effector translocation assay confirmed that the function of T4SS in NIΔcspA is influenced by deletion of the cspA gene. These results revealed the differential phenomena associated with virulence and metabolism in NIΔcspA and NI, providing important information for understanding detailed CspA-regulated interaction networks and Brucella pathogenesis.

Interactions between closely related bacterial strains are revealed by deep transcriptome sequencing

Comparative genomics, metagenomics, and single-cell technologies have shown that populations of microbial species encompass assemblages of closely related strains. This raises the question of whether individual bacterial lineages respond to the presence of their close relatives by modifying their gene expression or, instead, whether assemblages simply act as the arithmetic addition of their individual components. Here, we took advantage of transcriptome sequencing to address this question. For this, we analyzed the transcriptomes of two closely related strains of the extremely halophilic bacterium Salinibacter ruber grown axenically and in coculture. These organisms dominate bacterial assemblages in hypersaline environments worldwide. The strains used here cooccurred in the natural environment and are 100% identical in their 16S rRNA genes, and each strain harbors an accessory genome representing 10% of its complete genome. Overall, transcriptomic patterns from pure cultures were very similar for both strains. Expression was detected along practically the whole genome albeit with some genes at low levels. A subset of genes was very highly expressed in both strains, including genes coding for the light-driven proton pump xanthorhodopsin, genes involved in the stress response, and genes coding for transcriptional regulators. Expression differences between pure cultures affected mainly genes involved in environmental sensing. When the strains were grown in coculture, there was a modest but significant change in their individual transcription patterns compared to those in pure culture. Each strain sensed the presence of the other and responded in a specific manner, which points to fine intraspecific transcriptomic modulation.

RNA-seq reveals the critical role of OtpR in regulating Brucella melitensis metabolism and virulence under acidic stress

The response regulator OtpR is critical for the growth, morphology and virulence of Brucella melitensis. Compared to its wild type strain 16 M, B. melitensis 16 MΔotpR mutant has decreased tolerance to acid stress. To analyze the genes regulated by OtpR under acid stress, we performed RNA-seq whole transcriptome analysis of 16 MΔotpR and 16 M. In total, 501 differentially expressed genes were identified, including 390 down-regulated and 111 up-regulated genes. Among these genes, 209 were associated with bacterial metabolism, including 54 genes involving carbohydrate metabolism, 13 genes associated with nitrogen metabolism, and seven genes associated with iron metabolism. The 16 MΔotpR also decreased capacity to utilize different carbon sources and to tolerate iron limitation in culture experiments. Notably, OtpR regulated many Brucella virulence factors essential for B. melitensis intracellular survival. For instance, the virB operon encoding type IV secretion system was significantly down-regulated, and 36 known transcriptional regulators (e.g., vjbR and blxR) were differentially expressed in 16 MΔotpR. Selected RNA-seq results were experimentally confirmed by RT-PCR and RT-qPCR. Overall, these results deciphered differential phenomena associated with virulence, environmental stresses and cell morphology in 16 MΔotpR and 16 M, which provided important information for understanding the detailed OtpR-regulated interaction networks and Brucella pathogenesis.

Microevolution of Virulence-Related Genes in Helicobacter pylori Familial Infection

Helicobacter pylori, a bacterial pathogen that can infect human stomach causing gastritis, ulcers and cancer, is known to have a high degree of genome/epigenome diversity as the result of mutation and recombination. The bacteria often infect in childhood and persist for the life of the host. One of the reasons of the rapid evolution of H. pylori is that it changes its genome drastically for adaptation to a new host. To investigate microevolution and adaptation of the H. pylori genome, we undertook whole genome sequencing of the same or very similar sequence type in multi-locus sequence typing (MLST) with seven genes in members of the same family consisting of parents and children in Japan. Detection of nucleotide substitutions revealed likely transmission pathways involving children. Nonsynonymous (amino acid changing) mutations were found in virulence-related genes (cag genes, vacA, hcpDX, tnfα, ggt, htrA and the collagenase gene), outer membrane protein (OMP) genes and other cell surface-related protein genes, signal transduction genes and restriction-modification genes. We reconstructed various pathways by which H. pylori can adapt to a new human host, and our results raised the possibility that the mutational changes in virulence-related genes have a role in adaptation to a child host. Changes in restriction-modification genes might remodel the methylome and transcriptome to help adaptation. This study has provided insights into H. pylori transmission and virulence and has implications for basic research as well as clinical practice.

Global Gene Expression of Kosteletzkya virginica Seedlings Responding to Salt Stress

Soil salinization is becoming a serious threat to crop yield all over the world. Nowadays, acquainting the specific molecular mechanisms underlying various abiotic stresses especially to salt stress should be of great importance. While the development of the high-throughout sequencing technology promoted the progress powerfully. The intricate perception, transduction and regulation mechanisms underlying salt stress are being illustrated more and more clearly. As a perennial halophytic plant, Kosteletzkya virginica is able to help us to understand the mechanisms more directly and effectively. We carried out the whole transcriptome analysis on young seedlings with or without salt treatment through high-throughout sequencing technology. The results revealed that the numbers of different expressed transcripts between control and different treatments are 4145 and 9134, respectively. The ORF prediction suggested that there were 94308 ORF out of the 103489 (91.10%) total transcripts. We also carried out further differential expression analysis through gene ontology (GO) classification, cluster of orthologous groups (COG) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis. In a word, our transcriptome study on Kosteletzkya virginica would provide direct and effective references for researches on molecular mechanisms of salt-tolerance, extending our view of salt tolerance in plant further. Above all, the related report in this paper is the first about Kosteletzkya virginica.

Transcriptomic analysis of three Veillonella spp. present in carious dentine and in the saliva of caries-free individuals

Veillonella spp. are predominant bacteria found in all oral biofilms. In this study, a metatranscriptomic approach was used to investigate the gene expression levels of three oral Veillonella spp. (V. parvula, V. dispar and V. atypica) in whole stimulated saliva from caries-free volunteers and in carious lesions (n = 11 for each group). In the lesions the greatest proportion of reads were assigned to V. parvula and genes with the highest level of expression in carious samples were those coding for membrane transport systems. All three Veillonella spp. increased expression of genes involved in the catabolism of lactate and succinate, notably the alpha- and beta-subunits of L(+)-tartrate dehydratase (EC 4.2.1.32). There was also significantly increased expression of histidine biosynthesis pathway in V. parvula, suggesting higher intra-cellular levels of histidine that could provide intra-cellular buffering capacity and, therefore, assist survival in the acidic environment. Various other systems such as potassium uptake systems were also up regulated that may aid in the survival and proliferation of V. parvula in carious lesions.

Analysis of strand-specific RNA-seq data using machine learning reveals the structures of transcription units in Clostridium thermocellum

Identification of transcription units (TUs) encoded in a bacterial genome is essential to elucidation of transcriptional regulation of the organism. To gain a detailed understanding of the dynamically composed TU structures, we have used four strand-specific RNA-seq (ssRNA-seq) datasets collected under two experimental conditions to derive the genomic TU organization of Clostridium thermocellum using a machine-learning approach. Our method accurately predicted the genomic boundaries of individual TUs based on two sets of parameters measuring the RNA-seq expression patterns across the genome: expression-level continuity and variance. A total of 2590 distinct TUs are predicted based on the four RNA-seq datasets. Among the predicted TUs, 44% have multiple genes. We assessed our prediction method on an independent set of RNA-seq data with longer reads. The evaluation confirmed the high quality of the predicted TUs. Functional enrichment analyses on a selected subset of the predicted TUs revealed interesting biology. To demonstrate the generality of the prediction method, we have also applied the method to RNA-seq data collected on Escherichia coli and achieved high prediction accuracies. The TU prediction program named SeqTU is publicly available at https://code.google.com/p/seqtu/. We expect that the predicted TUs can serve as the baseline information for studying transcriptional and post-transcriptional regulation in C. thermocellum and other bacteria.

Integrated transcriptomic and proteomic analysis of the global response of Synechococcus to high light stress

Sufficient light is essential for the growth and physiological functions of photosynthetic organisms, but prolonged exposure to high light (HL) stress can cause cellular damage and ultimately result in the death of these organisms. Synechococcus sp. PCC 7002 (hereafter Synechococcus 7002) is a unicellular cyanobacterium with exceptional tolerance to HL intensities. However, the molecular mechanisms involved in HL response by Synechococcus 7002 are not well understood. Here, an integrated RNA sequencing transcriptomic and quantitative proteomic analysis was performed to investigate the cellular response to HL in Synechococcus 7002. A total of 526 transcripts and 233 proteins were identified to be differentially regulated under HL stress. Data analysis revealed major changes in mRNAs and proteins involved in the photosynthesis pathways, resistance to light-induced damage, DNA replication and repair, and energy metabolism. A set of differentially expressed mRNAs and proteins were validated by quantitative RT-PCR and Western blot, respectively. Twelve genes differentially regulated under HL stress were selected for knockout generation and growth analysis of these mutants led to the identification of key genes involved in the response of HL in Synechococcus 7002. Taken altogether, this study established a model for global response mechanisms to HL in Synechococcus 7002 and may be valuable for further studies addressing HL resistance in photosynthetic organisms.

Choice of next-generation sequencing pipelines

The next-generation sequencing (NGS) technologies are revolutionary tools which have made possible achieving remarkable advances in genetics since the beginning of the twenty-first century. Thanks to the possibility to produce large amount of sequence data, these tools are going to completely substitute other high-throughput technologies. Moreover, the large applications of NGS protocols are increasing the genetic decoding of biological systems through studies of genome anatomy and gene mapping, coupled to the transcriptome pictures. The application of NGS pipelines such as (1) de-novo genomic sequencing by mate-paired and whole-genome shotgun strategies; (2) specific gene sequencing on large bacterial communities; and (3) RNA-seq methods including whole transcriptome sequencing and Serial Analysis of Gene Expression (Sage-analysis) are fundamental in the genome-wide fields like metagenomics. Recently, the availability of these advanced protocols has allowed to overcome the usual sequencing technical issues related to the mapping specificity over standard shotgun library sequencing, the detection of large structural genomes variations and bridging sequencing gaps, as well as more precise gene annotation. In this chapter we will discuss how to manage a successful NGS pipeline from the planning of sequencing projects through the choice of the platforms up to the data analysis management.

Omics profiles used to evaluate the gene expression of Exiguobacterium antarcticum B7 during cold adaptation

BACKGROUND

Exiguobacterium antarcticum strain B7 is a Gram-positive psychrotrophic bacterial species isolated in Antarctica. Although this bacteria has been poorly studied, its genome has already been sequenced. Therefore, it is an appropriate model for the study of thermal adaptation. In the present study, we analyzed the transcriptomes and proteomes of E. antarcticum B7 grown at 0°C and 37°C by SOLiD RNA-Seq, Ion Torrent RNA-Seq and two-dimensional difference gel electrophoresis tandem mass spectrometry (2D-DIGE-MS/MS).

RESULTS

We found expression of 2,058 transcripts in all replicates from both platforms and differential expression of 564 genes (absolute log2FC≥1, P-value<0.001) comparing the two temperatures by RNA-Seq. A total of 73 spots were differentially expressed between the two temperatures on 2D-DIGE, 25 of which were identified by MS/MS. Some proteins exhibited patterns of dispersion in the gel that are characteristic of post-translational modifications.

CONCLUSIONS

Our findings suggest that the two sequencing platforms yielded similar results and that different omic approaches may be used to improve the understanding of gene expression. To adapt to low temperatures, E. antarcticum B7 expresses four of the six cold-shock proteins present in its genome. The cold-shock proteins were the most abundant in the bacterial proteome at 0°C. Some of the differentially expressed genes are required to preserve transcription and translation, while others encode proteins that contribute to the maintenance of the intracellular environment and appropriate protein folding. The results denote the complexity intrinsic to the adaptation of psychrotrophic organisms to cold environments and are based on two omic approaches. They also unveil the lifestyle of a bacterial species isolated in Antarctica.

RNA sequencing provides evidence for functional variability between naturally co-existing Alteromonas macleodii lineages

BACKGROUND

Alteromonas macleodii is a ubiquitous gammaproteobacterium shown to play a biogeochemical role in marine environments. Two A. macleodii strains (AltDE and AltDE1) isolated from the same sample (i.e., the same place at the same time) show considerable genomic differences. In this study, we investigate the transcriptional response of these two strains to varying growth conditions in order to investigate differences in their ability to adapt to varying environmental parameters.

RESULTS

RNA sequencing revealed transcriptional changes between all growth conditions examined (e.g., temperature and medium) as well as differences between the two A. macleodii strains within a given condition. The main inter-strain differences were more marked in the adaptation to grow on minimal medium with glucose and, even more so, under starvation. These differences suggested that AltDE1 may have an advantage over AltDE when glucose is the major carbon source, and co-culture experiments confirmed this advantage. Additional differences were observed between the two strains in the expression of ncRNAs and phage-related genes, as well as motility.

CONCLUSIONS

This study shows that the genomic diversity observed in closely related strains of A. macleodii from a single environment result in different transcriptional responses to changing environmental parameters. This data provides additional support for the idea that greater diversity at the strain level of a microbial community could enhance the community's ability to adapt to environmental shifts.

Global transcriptome analysis of hexavalent chromium stress responses in Staphylococcus aureus LZ-01

Staphylococcus aureus strain LZ-01, isolated from the Lanzhou reaches of the Yellow River, is capable of reducing Cr(VI) to Cr(III) aerobically. We employed transcriptome sequencing analysis to identify genes involved in Cr(VI) stress responses in S. aureus LZ-01. Our results showed that 512 of the 2,370 predicted genes displayed up-regulation (>2-fold), and 49 genes were down-regulated (<50 %) after Cr(VI) exposure. Among up-regulated genes, 128 genes were annotated to encode proteins involved in cellular processes; 68 were categorized to transport and binding proteins; 26 were involved in DNA repair and 32 were associated with regulatory functions. To further elucidate the Cr(VI) resistance and reduction mechanism, we carried out physiological tests and quantitative PCR analysis. Both RNA-seq and qRT-PCR data showed genes encoding a thioredoxin reductase and main subunits of cytochrome c oxidase complex were up-regulated upon Cr(VI) treatment. Either cadmium or NaN3 treatment could inhibit Cr(VI) reduction which indicates that thioredoxin and cytochrome are involved in Cr(VI) reduction strain LZ-01. 29 ABC-type metal/multidrug transporters and efflux pumps were up-regulated, suggesting that they are involved in Cr(VI) resistance by pumping chromium ions out of cells. The up-regulation of 26 DNA repair genes demonstrate that Cr(VI) is toxic to DNA and those DNA protection proteins need to be responded for Cr(VI) stress. Based on these results, the mechanism of strain LZ-01 resists and reduces Cr(VI) is revealed.

Advancements toward a systems level understanding of the human oral microbiome

Oral microbes represent one of the most well studied microbial communities owing to the fact that they are a fundamental part of human development influencing health and disease, an easily accessible human microbiome, a highly structured and remarkably resilient biofilm as well as a model of bacteria-bacteria and bacteria-host interactions. In the last 80 years since oral plaque was first characterized for its functionally stable physiological properties such as the highly repeatable rapid pH decrease upon carbohydrate addition and subsequent recovery phase, the fundamental approaches to study the oral microbiome have cycled back and forth between community level investigations and characterizing individual model isolates. Since that time, many individual species have been well characterized and the development of the early plaque community, which involves many cell-cell binding interactions, has been carefully described. With high throughput sequencing enabling the enormous diversity of the oral cavity to be realized, a number of new challenges to progress were revealed. The large number of uncultivated oral species, the high interpersonal variability of taxonomic carriage and the possibility of multiple pathways to dysbiosis pose as major hurdles to obtain a systems level understanding from the community to the gene level. It is now possible however to start connecting the insights gained from single species with community wide approaches. This review will discuss some of the recent insights into the oral microbiome at a fundamental level, existing knowledge gaps, as well as challenges that have surfaced and the approaches to address them.

Transcriptome analysis of Acidovorax avenae subsp. avenae cultivated in vivo and co-culture with Burkholderia seminalis

Response of bacterial pathogen to environmental bacteria and its host is critical for understanding of microbial adaption and pathogenesis. Here, we used RNA-Seq to comprehensively and quantitatively assess the transcriptional response of Acidovorax avenae subsp. avenae strain RS-1 cultivated in vitro, in vivo and in co-culture with rice rhizobacterium Burkholderia seminalis R456. Results revealed a slight response to other bacteria, but a strong response to host. In particular, a large number of virulence associated genes encoding Type I to VI secretion systems, 118 putative non-coding RNAs, and 7 genomic islands (GIs) were differentially expressed in vivo based on comparative genomic and transcriptomic analyses. Furthermore, the loss of virulence for knockout mutants of 11 differentially expressed T6SS genes emphasized the importance of these genes in bacterial pathogenicity. In addition, the reliability of expression data obtained by RNA-Seq was supported by quantitative real-time PCR of the 25 selected T6SS genes. Overall, this study highlighted the role of differentially expressed genes in elucidating bacterial pathogenesis based on combined analysis of RNA-Seq data and knockout of T6SS genes.

Integrated transcriptomic and proteomic analysis of the bile stress response in a centenarian-originated probiotic Bifidobacterium longum BBMN68

Bifidobacteria are natural inhabitants of the human gastrointestinal tract and well known for their health-promoting effects. Tolerance to bile stress is crucial for bifidobacteria to survive in the colon and to exert their beneficial actions. In this work, RNA-Seq transcriptomic analysis complemented with proteomic analysis was used to investigate the cellular response to bile in Bifidobacterium longum BBMN68. The transcript levels of 236 genes were significantly changed (≥ threefold, p < 0.001) and 44 proteins were differentially abundant (≥1.6-fold, p < 0.01) in B. longum BBMN68 when exposed to 0.75 g l(-1) ox-bile. The hemolysin-like protein and bile efflux systems were significantly over produced, which might prevent bile adsorption and exclude bile, respectively. The cell membrane composition was modified probably by an increase of cyclopropane fatty acid and a decrease of transmembrane proteins, resulting in a cell membrane more impermeable to bile salts. Our hypothesis was later confirmed by surface hydrophobicity assay. The transcription of genes related to xylose utilization and bifid shunt were up-regulated, which increased the production of ATP and reducing equivalents to cope with bile-induced damages in a xylan-rich colon environment. Bile salts signal the B. longum BBMN68 to gut entrance and enhance the expression of esterase and sortase associated with adhesion and colonization in intestinal tract, which was supported by a fivefold increased adhesion ability to HT-29 cells by BBMN68 upon bile exposure. Notably, bacterial one-hybrid and EMSA assay revealed that the two-component system senX3-regX3 controlled the expression of pstS in bifidobacteria and the role of this target gene in bile resistance was further verified by heterologous expression in Lactococcus lactis. Taken altogether, this study established a model for global response mechanisms in B. longum to bile.

Transcriptomic analysis of Staphylococcus aureus using microarray and advanced next-generation RNA-seq technologies

The transcriptome has shown tremendous potential for the comprehensive investigation of gene expression profiles and transcriptional levels in comparative biology, the identification of regulatory mechanism of transcriptional regulators, and the evaluation of target gene for developing new chemotherapeutic agents, vaccine, and diagnostic methods. The traditional microarray and advanced next-generation RNA sequencing technologies (RNA-seq) provide powerful and effective tools for the determination of the transcriptome of bacterial cells. In this chapter, we provide a detailed protocol for scientists who want to investigate gene expression profiles by performing microarray and/or RNA-seq analysis, including different RNA purification methods, mRNA enrichment, decontamination, cDNA synthesis, fragmentation, biotin labeling for hybridization using Affymetrix Staphylococcus aureus chips, quantitative real-time reverse transcription PCR, and RNA-seq data analysis.

SearchDOGS bacteria, software that provides automated identification of potentially missed genes in annotated bacterial genomes

We report the development of SearchDOGS Bacteria, software to automatically detect missing genes in annotated bacterial genomes by combining BLAST searches with comparative genomics. Having successfully applied the approach to yeast genomes, we redeveloped SearchDOGS to function as a standalone, downloadable package, requiring only a set of GenBank annotation files as input. The software automatically generates a homology structure using reciprocal BLAST and a synteny-based method; this is followed by a scan of the entire genome of each species for unannotated genes. Results are provided in a HTML interface, providing coordinates, BLAST results, syntenic location, omega values (Ka/Ks, where Ks is the number of synonymous substitutions per synonymous site and Ka is the number of nonsynonymous substitutions per nonsynonymous site) for protein conservation estimates, and other information for each candidate gene. Using SearchDOGS Bacteria, we identified 155 gene candidates in the Shigella boydii sb227 genome, including 56 candidates of length < 60 codons. SearchDOGS Bacteria has two major advantages over currently available annotation software. First, it outperforms current methods in terms of sensitivity and is highly effective at identifying small or highly diverged genes. Second, as a freely downloadable package, it can be used with unpublished or confidential data.

Transcriptome analysis of Vibrio parahaemolyticus in type III secretion system 1 inducing conditions

Vibrio parahaemolyticus is an emerging bacterial pathogen capable of causing inflammatory gastroenteritis, wound infections, and septicemia. As a food-borne illness, infection is most frequently associated with the consumption of raw or undercooked seafood, particularly shellfish. It is the primary cause of Vibrio-associated food-borne illness in the United States and the leading cause of food-borne illness in Japan. The larger of its two chromosomes harbors a set of genes encoding type III section system 1 (T3SS1), a virulence factor present in all V. parahaemolyticus strains that is similar to the Yersinia ysc T3SS. T3SS1 translocates effector proteins into eukaryotic cells where they induce changes to cellular physiology and modulate host-pathogen interactions. T3SS1 is also responsible for cytotoxicity toward several different cultured cell lines as well as mortality in a mouse model. Herein we used RNA-seq to obtain global transcriptome patterns of V. parahaemolyticus under conditions that either induce [growth in Dulbecco's Modified Eagle Medium (DMEM) media, in trans expression of transcriptional regulator exsA] or repress T3SS1 expression (growth in LB-S media, in trans exsD expression) and during infection of HeLa cells over time. Comparative transcriptomic analysis demonstrated notable differences in the expression patterns under inducing conditions and was also used to generate an expression profile of V. parahaemolyticus during infection of HeLa cells. In addition, we identified several new genes that are associated with T3SS1 expression and may warrant further study.

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.

Antimicrobial nodule-specific cysteine-rich peptides induce membrane depolarization-associated changes in the transcriptome of Sinorhizobium meliloti

Leguminous plants establish symbiosis with nitrogen-fixing alpha- and betaproteobacteria, collectively called rhizobia, which provide combined nitrogen to support plant growth. Members of the inverted repeat-lacking clade of legumes impose terminal differentiation on their endosymbiotic bacterium partners with the help of the nodule-specific cysteine-rich (NCR) peptide family composed of close to 600 members. Among the few tested NCR peptides, cationic ones had antirhizobial activity measured by reduction or elimination of the CFU and uptake of the membrane-impermeable dye propidium iodide. Here, the antimicrobial spectrum of two of these peptides, NCR247 and NCR335, was investigated, and their effect on the transcriptome of the natural target Sinorhizobium meliloti was characterized. Both peptides were able to kill quickly a wide range of Gram-negative and Gram-positive bacteria; however, their spectra were only partially overlapping, and differences were found also in their efficacy on given strains, indicating that the actions of NCR247 and NCR335 might be similar though not identical. Treatment of S. meliloti cultures with either peptide resulted in a quick downregulation of genes involved in basic cellular functions, such as transcription-translation and energy production, as well as upregulation of genes involved in stress and oxidative stress responses and membrane transport. Similar changes provoked mainly in Gram-positive bacteria by antimicrobial agents were coupled with the destruction of membrane potential, indicating that it might also be a common step in the bactericidal actions of NCR247 and NCR335.

Detection of phase-dependent transcriptomic changes and Rubisco-mediated CO2 fixation into poly (3-hydroxybutyrate) under heterotrophic condition in Ralstonia eutropha H16 based on RNA-seq and gene deletion analyses

Background

Ralstonia eutropha H16 is well known to produce polyhydroxyalkanoates (PHAs), which are potential bio-based biodegradable plastics, in an efficient manner as an energy storage material under unbalanced growth conditions. To obtain further knowledge of PHA biosynthesis, this study performed a quantitative transcriptome analysis based on deep sequencing of the complementary DNA generated from the RNA (RNA-seq) of R. eutropha H16.

Results

Total RNAs were extracted from R. eutropha cells in growth, PHA production, and stationary phases on fructose. rRNAs in the preparation were removed by repeated treatments with magnetic beads specific to bacterial rRNAs, and then the 36 bp sequences were determined using an Illumina high-throughput sequencer. The RNA-seq results indicated the induction of gene expression for transcription, translation, cell division, peptidoglycan biosynthesis, pilus and flagella assembly, energy conservation, and fatty acid biosynthesis in the growth phase; and the repression trends of genes involved in central metabolisms in the PHA production phase. Interestingly, the transcription of genes for Calvin-Benson-Bassham (CBB) cycle and several genes for β-oxidation were significantly induced in the PHA production phase even when the cells were grown on fructose. Moreover, incorporation of ¹³C was observed in poly(3-hydroxybutyrate) synthesized by R. eutropha H16 from fructose in the presence of NaH¹³CO₃, and further gene deletion analyses revealed that both of the two ribulose 1,5-bisphosphate carboxylase (Rubiscos) in CBB cycle were actually functional in CO₂ fixation under the heterotrophic condition.

Conclusions

The results revealed the phase-dependent transcriptomic changes and a CO₂ fixation capability under heterotrophic conditions by PHA-producing R. eutropha.

Current challenges in bacterial transcriptomics

Over the past decade or so, dramatic developments in our ability to experimentally determine the content and function of genomes have taken place. In particular, next-generation sequencing technologies are now inspiring a new understanding of bacterial transcriptomes on a global scale. In bacterial cells, whole-transcriptome studies have not received attention, owing to the general view that bacterial genomes are simple. However, several recent RNA sequencing results are revealing unexpected levels of complexity in bacterial transcriptomes, indicating that the transcribed regions of genomes are much larger and complex than previously anticipated. In particular, these data show a wide array of small RNAs, antisense RNAs, and alternative transcripts. Here, we review how current transcriptomics are now revolutionizing our understanding of the complexity and regulation of bacterial transcriptomes.

Transcriptome profiling analysis of Vibrio vulnificus during human infection

Vibrio vulnificus is a waterborne pathogen that was responsible for an outbreak of severe soft-tissue infections among fish farmers and fish consumers in Israel. Several factors have been shown to be associated with virulence. However, the transcriptome profile of the pathogen during human infection has not been determined yet. We compared the transcriptome profile, using RNA sequencing, of a human-pathogenic strain harvested directly from tissue of a patient suffering from severe soft-tissue infection with necrotizing fasciitis, with the same strain and three other environmental strains grown in vitro. The five sequenced libraries were aligned to the reference genomes of V. vulnificus strains CMCP6 and YJ016. Approximately 47.8 to 62.3 million paired-end raw reads were generated from the five runs. Nearly 84 % of the genome was covered by reads from at least one of the five runs, suggesting that nearly 16 % of the genome is not transcribed or is transcribed at low levels. We identified 123 genes that were differentially expressed during the acute phase of infection. Sixty-three genes were mapped to the large chromosome, 47 genes mapped to the small chromosome and 13 genes mapped to the YJ016 plasmid. The 123 genes fell into a variety of functional categories including transcription, signal transduction, cell motility, carbohydrate metabolism, intracellular trafficking and cell envelope biogenesis. Among the genes differentially expressed during human infection we identified genes encoding bacterial toxin (RtxA1) and genes involved in flagellar components, Flp-coding region, GGDEF family protein, iron acquisition system and sialic acid metabolism.

Imaging of the skeletal muscle metastases

Preclinical Research

Copyright 2011 Wiley-Liss, Inc., A Wiley CompanyThis article is being made freely available through PubMed Central as part of the COVID-19 public health emergency response. It can be used for unrestricted research re-use and analysis in any form or by any means with acknowledgement of the original source, for the duration of the public health emergency.

Omics technologies include genomics, transcriptomics, proteomics, metabolomics, and immunomics. These technologies have been used in vaccine research, which can be summarized using the term “vaccinomics.” These omics technologies combined with advanced bioinformatics analysis form the core of “systems vaccinology.” Omics technologies provide powerful methods in vaccine target identification. The genomics‐based reverse vaccinology starts with predicting vaccine protein candidates through in silico bioinformatics analysis of genome sequences. The VIOLIN Vaxign vaccine design program (http://www.violinet.org/vaxign) is the first web‐based vaccine target prediction software based on the reverse vaccinology strategy. Systematic transcriptomics and proteomics analyses facilitate rational vaccine target identification by detesting genome‐wide gene expression profiles. Immunomics is the study of the set of antigens recognized by host immune systems and has also been used for efficient vaccine target prediction. With the large amount of omics data available, it is necessary to integrate various vaccine data using ontologies, including the Gene Ontology (GO) and Vaccine Ontology (VO), for more efficient vaccine target prediction and assessment. All these omics technologies combined with advanced bioinformatics analysis methods for a systems biology‐based vaccine target prediction strategy. This article reviews the various omics technologies and how they can be used in vaccine target identification.

cDNA normalization by hydroxyapatite chromatography to enrich transcriptome diversity in RNA-seq applications

Second-generation sequencing (SGS) has become the preferred method for RNA transcriptome profiling of organisms and single cells. However, SGS analysis of transcriptome diversity (including protein-coding transcripts and regulatory non-coding RNAs) is inefficient unless the sample of interest is first depleted of nucleic acids derived from ribosomal RNA (rRNA), which typically account for up to 95% of total intracellular RNA content. Here we describe a novel microscale hydroxyapatite chromatography (HAC) normalization method to remove eukaryotic and prokaryotic high abundant rRNA species, thereby increasing sequence coverage depth and transcript diversity across non-rRNA populations. RNA-seq analysis of Escherichia coli K-12 and human intracellular total RNA showed that HAC-based normalization enriched for all non-ribosomal RNA species regardless of RNA transcript abundance or length when compared with untreated controls. Microcolumn HAC normalization generated rRNA-depleted cDNA libraries comparable to the well-established duplex specific nuclease (DSN) normalization and Ribo-Zero rRNA-depletion methods, thus establishing microscale HAC as an effective, cost saving, and non-destructive alternative normalization technique.

Emerging tools for synthetic genome design

Synthetic biology is an emerging discipline for designing and synthesizing predictable, measurable, controllable, and transformable biological systems. These newly designed biological systems have great potential for the development of cheaper drugs, green fuels, biodegradable plastics, and targeted cancer therapies over the coming years. Fortunately, our ability to quickly and accurately engineer biological systems that behave predictably has been dramatically expanded by significant advances in DNA-sequencing, DNA-synthesis, and DNA-editing technologies. Here, we review emerging technologies and methodologies in the field of building designed biological systems, and we discuss their future perspectives.

Next-generation sequencing technologies and their impact on microbial genomics

Next-generation sequencing technologies have had a dramatic impact in the field of genomic research through the provision of a low cost, high-throughput alternative to traditional capillary sequencers. These new sequencing methods have surpassed their original scope and now provide a range of utility-based applications, which allow for a more comprehensive analysis of the structure and content of microbial genomes than was previously possible. With the commercialization of a third generation of sequencing technologies imminent, we discuss the applications of current next-generation sequencing methods and explore their impact on and contribution to microbial genome research.

Multidisciplinary Research Agenda for Novel Antimicrobial Agents for Caries Prevention and Treatment

Antimicrobial methods to augment fluoride-mediated caries inhibition are necessary. Several methods are described here, but none was considered likely to be as effective as fluoride usage. None had been tested in effective models to demonstrate their ability to act either additively or synergistically with fluoride-containing toothpastes. Dental caries is a biofilm-mediated disease: The composition of the biofilm associated with caries initiation and progression is diverse. Caries-associated taxa – including mutans streptococci, lactobacilli, bifidobacteria, and yeasts – may be useful surrogate markers for in vivo investigations. In vitro testing should progress from single-species planktonic cells to multi-species biofilms prior to essential testing in randomized control trials (RCTs). Modern high-throughput sequencing techniques need to be applied to the study of bacterial acquisition from birth and of the composition of the biofilm associated with the formation of white-spot lesions. The determination of the functions of the biofilm and the phenotype of the bacterial components may be determined by RNA-seq techniques, since they must be conserved between caries lesions and will include the ability to produce acids and survive and proliferate in acidic conditions. The application of such methods will significantly improve our understanding of the etiology and progression of dental caries.

FunSys: Software for functional analysis of prokaryotic transcriptome and proteome

The vast amount of data produced by next-generation sequencing (NGS) has necessitated the development of computational tools to assist in understanding the myriad functions performed by the biological macromolecules involved in heredity. In this work, we developed the FunSys programme, a stand-alone tool with an user friendly interface that enables us to evaluate and correlate differential expression patterns from RNA sequencing and proteomics datasets. The FunSys generates charts and reports based on the results of the analysis of differential expression to aid the interpretation of the results.

AVAILABILITY

The database is available for free at https://sourceforge.net/projects/funsysufpa/

The promoter architectural landscape of the Salmonella PhoP regulon

The DNA-binding protein PhoP controls virulence and Mg²⁺ homeostasis in the Gram-negative pathogen Salmonella enterica serovar Typhimurium. PhoP regulates expression of a large number of genes that differ both in their ancestry and in the biochemical functions and physiological roles of the encoded products. This suggests that PhoP-regulated genes are differentially expressed. To understand how a bacterial activator might generate varied gene expression behaviour, we investigated the cis-acting promoter features (i.e. the number of PhoP binding sites, as well as their orientation and location with respect to the sites bound by RNA polymerase and the sequences that constitute the PhoP binding sites) in 23 PhoP-activated promoters. Our results show that natural PhoP-activated promoters utilize only a limited number of combinations of cis-acting features – or promoter architectures. We determine that PhoP activates transcription by different mechanisms, and that ancestral and horizontally acquired PhoP-activated genes have distinct promoter architectures.

Response to copper stress in Streptomyces lividans extends beyond genes under direct control of a copper-sensitive operon repressor protein (CsoR)

A copper-sensitive operon repressor protein (CsoR) has been identified in Streptomyces lividans (CsoR(Sl)) and found to regulate copper homeostasis with attomolar affinity for Cu(I). Solution studies reveal apo- and Cu(I)-CsoR(Sl) to be a tetramer assembly, and a 1.7-Å resolution crystal structure of apo-CsoR(Sl) reveals that a significant conformational change is necessary to enable Cu(I) binding. In silico prediction of the CsoR regulon was confirmed in vitro (EMSA) and in vivo (RNA-seq), which highlighted that next to the csoR gene itself, the regulon consists of two Cu(I) efflux systems involving a CopZ-like copper metallochaperone protein and a CopA P(1)-type ATPase. Although deletion of csoR has only minor effects on S. lividans development when grown under high copper concentrations, mutations of the Cu(I) ligands decrease tolerance to copper as a result of the Cu(I)-CsoR mutants failing to disengage from the DNA targets, thus inhibiting the derepression of the regulon. RNA-seq experiments carried out on samples incubated with exogenous copper and a ΔcsoR strain showed that the set of genes responding to copper stress is much wider than anticipated and largely extends beyond genes targeted by CsoR. This suggests more control levels are operating and directing other regulons in copper homeostasis beside the CsoR regulon.