RNA-Seq reveals differential gene expression in Staphylococcus aureus with single-nucleotide resolution

Hybrid assembly using long reads resolves repeats and completes the genome sequence of a laboratory strain of Staphylococcus aureus subsp. Aureus RN4220

Staphylococcus aureus RN4220 has been extensively used by staphylococcal researchers as an intermediate strain for genetic manipulation due to its ability to accept foreign DNA. Despite its wide use in laboratories, its complete genome is not available. In this study, we used a hybrid genome assembly approach using minION long reads and Illumina short reads to sequence the complete genome of S. aureus RN4220. The comparative analysis of the annotated complete genome showed the presence of 39 genes fragmented in the previous assembly, many of which were located near the repeat regions. Using RNA-Seq reads, we showed that a higher number of reads could be mapped to the complete genome than the draft genome and the gene expression profile obtained using the complete genome also differs from that obtained from the draft genome. Furthermore, by comparative transcriptomic analysis, we showed the correlation between expression levels of staphyloxanthin biosynthetic genes and the production of yellow pigment. This study highlighted the importance of long reads in completing microbial genomes, especially those possessing repetitive elements.

•

S. aureus RN4220 is used as an intermediate strain for genetic manipulation.

•

We completed its genome and found 39 fragmented genes in previous genome assembly.

•

RNA-Seq analysis improved mapping of the reads with the use of complete genome.

•

Expression of staphyloxanthin biosynthetic genes was correlated with its production.

Leveraging big data bioinformatics approaches to extract knowledge from Staphylococcus aureus public omics data

Staphylococcus aureus is a notorious pathogen posing challenges in the medical industry due to drug resistance and biofilm formation. The horizon of knowledge on S. aureus pathogenesis has expanded with the advancement of data-driven bioinformatics techniques. Mining information from sequenced genomes and their expression data is an economic approach that alleviates wastage of resources and redundancy in experiments. The current review covers how big data bioinformatics has been used in the analysis of S. aureus from publicly available -omics data to uncover mechanisms of infection and inhibition. Particularly, advances in the past two decades in biomarker discovery, host responses, phenotype identification, consolidation of information, and drug development are discussed highlighting the challenges and shortcomings. Overall, the review summarizes the diverse aspects of scrupulous re-analysis of S. aureus proteomic and transcriptomic expression datasets retrieved from public repositories in terms of the efforts taken, benefits offered, and follow-up actions. The detailed review thus serves as a reference and aid for (i) Computational biologists by briefing the approaches utilized for bacterial omics re-analysis concerning S. aureus and (ii) Experimental biologists by elucidating the potential of bioinformatics in biological research to generate reliable postulates in a prompt and economical manner.

Analysis of the structures of confirmed and questionable CRISPR loci in 325 Staphylococcus genomes

The CRISPR-Cas (clustered regular interspaced short palindromic repeats and CRISPR-associated proteins) system is a newly discovered immune defense system in the genome of prokaryotes, which can resist the invasion of foreign genetic elements, such as plasmids or phage. In this study, 154 strains of Staphylococcus published in the CRISPRDatabase and 171 strains included in NCBI were downloaded, the confirmed and questionable CRISPR loci of which were analyzed by bioinformatics methods, including their distribution, characteristics of the structure (including the direct repeats, spacers and cas genes), and the relationship between the presence of CRISPR and the mecA gene. Meanwhile, a comprehensive analysis of orphan CRISPR arrays was performed on this basis. A total of 196 confirmed and 1757 questionable CRISPR loci were found in 325 Staphylococcus genomes. Only 25 strains contained cas genes, which were classified into III-A (48.1%) and II-C (51.9%). The difference between the presence of the cas gene and the carrying rate of mecA was statistically significant, and they were negatively correlated. A total of 137 confirmed and 1755 questionable CRISPR loci were assumed to be false-CRISPR. The present study also analyzed the questionable CRISPR array for the first time while analyzing the confirmed CRISPR array in the Staphylococcal genome and screened the false-CRISPR elements in the orphan CRISPR array.

Transcriptome analysis of Corynebacterium pseudotuberculosis biovar Equi in two conditions of the environmental stress

Corynebacterium pseudotuberculosis has been widely studied in an effort to understand its biological evolution. Transcriptomics has revealed possible candidates for virulence and pathogenicity factors of strain 1002 (biovar Ovis). Because C. pseudotuberculosis is classified into two biovars, Ovis and Equi, it was interesting to assess the transcriptional profile of biovar Equi strain 258, the causative agent of ulcerative lymphangitis. The genome of this strain was re-sequenced; the reassembly was completed using optical mapping technology, and the sequence was subsequently re-annotated. Two growth conditions that occur during the host infection process were simulated for the transcriptome: the osmotic and acid medium. Genes that may be associated with the microorganism's resilience under unfavorable conditions were identified through RNAseq, including genes present in pathogenicity islands. The RT-qPCR was performed to confirm the results in biological triplicate for each condition for some genes. The results extend our knowledge of the factors associated with the spread and persistence of C. pseudotuberculosis during the infection process and suggest possible avenues for studies related to the development of vaccines, diagnosis, and therapies that might help minimize damage to agribusinesses.

The Suf Iron-Sulfur Cluster Biosynthetic System Is Essential in Staphylococcus aureus, and Decreased Suf Function Results in Global Metabolic Defects and Reduced Survival in Human Neutrophils

Staphylococcus aureus remains a causative agent for morbidity and mortality worldwide. This is in part a result of antimicrobial resistance, highlighting the need to uncover novel antibiotic targets and to discover new therapeutic agents. In the present study, we explored the possibility that iron-sulfur (Fe-S) cluster synthesis is a viable antimicrobial target. RNA interference studies established that Suf (sulfur mobilization)-dependent Fe-S cluster synthesis is essential in S. aureus We found that sufCDSUB were cotranscribed and that suf transcription was positively influenced by sigma factor B. We characterized an S. aureus strain that contained a transposon inserted in the intergenic space between sufC and sufD (sufD*), resulting in decreased transcription of sufSUB Consistent with the transcriptional data, the sufD* strain had multiple phenotypes associated with impaired Fe-S protein maturation. They included decreased activities of Fe-S cluster-dependent enzymes, decreased growth in media lacking metabolites that require Fe-S proteins for synthesis, and decreased flux through the tricarboxylic acid (TCA) cycle. Decreased Fe-S cluster synthesis resulted in sensitivity to reactive oxygen and reactive nitrogen species, as well as increased DNA damage and impaired DNA repair. The sufD* strain also exhibited perturbed intracellular nonchelated Fe pools. Importantly, the sufD* strain did not exhibit altered exoprotein production or altered biofilm formation, but it was attenuated for survival upon challenge by human polymorphonuclear leukocytes. The results presented are consistent with the hypothesis that Fe-S cluster synthesis is a viable target for antimicrobial development.

Using Quantitative Spectrometry to Understand the Influence of Genetics and Nutritional Perturbations On the Virulence Potential of Staphylococcus aureus

Staphylococcus aureus (Sa) is the leading cause of a variety of bacterial infections ranging from superficial skin infections to invasive and life threatening diseases such as septic bacteremia, necrotizing pneumonia, and endocarditis. The success of Sa as a human pathogen is contributed to its ability to adapt to different environments by changing expression, production, or secretion of virulence factors. Although Sa immune evasion is well-studied, the regulation of virulence factors under different nutrient and growth conditions is still not well understood. Here, we used label-free quantitative mass spectrometry to quantify and compare the Sa exoproteins (i.e. exoproteomes) of master regulator mutants or established reference strains. Different environmental conditions were addressed by growing the bacteria in rich or minimal media at different phases of growth. We observed clear differences in the composition of the exoproteomes depending on the genetic background or growth conditions. The relative abundance of cytotoxins determined in our study correlated well with differences in cytotoxicity measured by lysis of human neutrophils. Our findings demonstrate that label-free quantitative mass spectrometry is a versatile tool for predicting the virulence of bacterial strains and highlights the importance of the experimental design for in vitro studies. Furthermore, the results indicate that label-free proteomics can be used to cluster isolates into groups with similar virulence properties, highlighting the power of label-free quantitative mass spectrometry to distinguish Sa strains.

Not all predicted CRISPR-Cas systems are equal: isolated cas genes and classes of CRISPR like elements

Background

The CRISPR–Cas systems in prokaryotes are RNA-guided immune systems that target and deactivate foreign nucleic acids. A typical CRISPR–Cas system consists of a CRISPR array of repeat and spacer units, and a locus of cas genes. The CRISPR and the cas locus are often located next to each other in the genomes. However, there is no quantitative estimate of the co-location. In addition, ad-hoc studies have shown that some non-CRISPR genomic elements contain repeat-spacer-like structures and are mistaken as CRISPRs.

Results

Using available genome sequences, we observed that a significant number of genomes have isolated cas loci and/or CRISPRs. We found that 11%, 22% and 28% of the type I, II and III cas loci are isolated (without CRISPRs in the same genomes at all or with CRISPRs distant in the genomes), respectively. We identified a large number of genomic elements that superficially reassemble CRISPRs but don’t contain diverse spacers and have no companion cas genes. We called these elements false-CRISPRs and further classified them into groups, including tandem repeats and Staphylococcus aureus repeat (STAR)-like elements.

Conclusion

This is the first systematic study to collect and characterize false-CRISPR elements. We demonstrated that false-CRISPRs could be used to reduce the false annotation of CRISPRs, therefore showing them to be useful for improving the annotation of CRISPR–Cas systems.

Electronic supplementary material

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

Genomic and transcriptomic comparison between Staphylococcus aureus strains associated with high and low within herd prevalence of intra-mammary infection

Background

Staphylococcus aureus (Staph. aureus) is one of the major pathogens causing mastitis in dairy ruminants worldwide. The chronic nature of Staph. aureus infection enhances the contagiousness risk and diffusion in herds. In order to identify the factors involved in intra-mammary infection (IMI) and diffusion in dairy cows, we investigated the molecular characteristics of two groups of Staph. aureus strains belonging to ST8 and ST398, differing in clinical properties, through comparison of whole genome and whole transcriptome sequencing.

Results

The two groups of strains, one originated from high IMI prevalence herds and the other from low IMI prevalence herds, present a peculiar set of genes and polymorphisms related to phenotypic features, such as bacterial invasion of mammary epithelial cells and host adaptation. Transcriptomic analysis supports the high propensity of ST8 strain to chronicity of infection and to a higher potential cytotoxicity.

Conclusions

Our data are consistent with the invasiveness and host adaptation feature for the strains GTB/ST8 associated to high within-herd prevalence of mastitis. Variation in genes coding for surface exposed proteins and those associated to virulence and defence could constitute good targets for further research.

Electronic supplementary material

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

A Transcriptome Map of Actinobacillus pleuropneumoniae at Single-Nucleotide Resolution Using Deep RNA-Seq

Actinobacillus pleuropneumoniae is the pathogen of porcine contagious pleuropneumoniae, a highly contagious respiratory disease of swine. Although the genome of A. pleuropneumoniae was sequenced several years ago, limited information is available on the genome-wide transcriptional analysis to accurately annotate the gene structures and regulatory elements. High-throughput RNA sequencing (RNA-seq) has been applied to study the transcriptional landscape of bacteria, which can efficiently and accurately identify gene expression regions and unknown transcriptional units, especially small non-coding RNAs (sRNAs), UTRs and regulatory regions. The aim of this study is to comprehensively analyze the transcriptome of A. pleuropneumoniae by RNA-seq in order to improve the existing genome annotation and promote our understanding of A. pleuropneumoniae gene structures and RNA-based regulation. In this study, we utilized RNA-seq to construct a single nucleotide resolution transcriptome map of A. pleuropneumoniae. More than 3.8 million high-quality reads (average length ~90 bp) from a cDNA library were generated and aligned to the reference genome. We identified 32 open reading frames encoding novel proteins that were mis-annotated in the previous genome annotations. The start sites for 35 genes based on the current genome annotation were corrected. Furthermore, 51 sRNAs in the A. pleuropneumoniae genome were discovered, of which 40 sRNAs were never reported in previous studies. The transcriptome map also enabled visualization of 5'- and 3'-UTR regions, in which contained 11 sRNAs. In addition, 351 operons covering 1230 genes throughout the whole genome were identified. The RNA-Seq based transcriptome map validated annotated genes and corrected annotations of open reading frames in the genome, and led to the identification of many functional elements (e.g. regions encoding novel proteins, non-coding sRNAs and operon structures). The transcriptional units described in this study provide a foundation for future studies concerning the gene functions and the transcriptional regulatory architectures of this pathogen.

Understanding bacterial resistance to antimicrobial peptides: From the surface to deep inside

Resistant bacterial infections are a major health problem in many parts of the world. The major commercial antibiotic classes often fail to combat common bacteria. Although antimicrobial peptides are able to control bacterial infections by interfering with microbial metabolism and physiological processes in several ways, a large number of cases of resistance to antibiotic peptide classes have also been reported. To gain a better understanding of the resistance process various technologies have been applied. Here we discuss multiple strategies by which bacteria could develop enhanced antimicrobial peptide resistance, focusing on sub-cellular regions from the surface to deep inside, evaluating bacterial membranes, cell walls and cytoplasmic metabolism. Moreover, some high-throughput methods for antimicrobial resistance detection and discrimination are also examined. This article is part of a Special Issue entitled: Bacterial Resistance to Antimicrobial Peptides.

Staphylococcus aureus competence genes: mapping of the SigH, ComK1 and ComK2 regulons by transcriptome sequencing

Staphylococcus aureus is a major human pathogen. Hospital infections caused by methicillin-resistant strains (MRSA), which have acquired resistance to a broad spectrum of antibiotics through horizontal gene transfer (HGT), are of particular concern. In S. aureus, virulence and antibiotic resistance genes are often encoded on mobile genetic elements that are disseminated by HGT. Conjugation and phage transduction have long been known to mediate HGT in this species, but it is unclear whether natural genetic transformation contributes significantly to the process. Recently, it was reported that expression of the alternative sigma factor SigH induces the competent state in S. aureus. The transformation efficiency obtained, however, was extremely low, indicating that the optimal conditions for competence development had not been found. We therefore used transcriptome sequencing to determine whether the full set of genes known to be required for competence in other naturally transformable bacteria is part of the SigH regulon. Our results show that several essential competence genes are not controlled by SigH. This presumably explains the low transformation efficiency previously reported, and demonstrates that additional regulating mechanisms must be involved. We found that one such mechanism involves ComK1, a transcriptional activator that acts synergistically with SigH.