Re-annotation of the Saccharopolyspora erythraea genome using a systems biology approach

Proteins à la carte: riboproteogenomic exploration of bacterial N-terminal proteoform expression

In recent years, it has become evident that the true complexity of bacterial proteomes remains underestimated. Gene annotation tools are known to propagate biases and overlook certain classes of truly expressed proteins, particularly proteoforms-protein isoforms arising from a single gene. Recent (re-)annotation efforts heavily rely on ribosome profiling by providing a direct readout of translation to fully describe bacterial proteomes. In this study, we employ a robust riboproteogenomic pipeline to conduct a systematic census of expressed N-terminal proteoform pairs, representing two isoforms encoded by a single gene raised by annotated and alternative translation initiation, in Salmonella. Intriguingly, conditional-dependent changes in relative utilization of annotated and alternative translation initiation sites (TIS) were observed in several cases. This suggests that TIS selection is subject to regulatory control, adding yet another layer of complexity to our understanding of bacterial proteomes.

IMPORTANCE

With the emerging theme of genes within genes comprising the existence of alternative open reading frames (ORFs) generated by translation initiation at in-frame start codons, mechanisms that control the relative utilization of annotated and alternative TIS need to be unraveled and our molecular understanding of resulting proteoforms broadened. Utilizing complementary ribosome profiling strategies to map ORF boundaries, we uncovered dual-encoding ORFs generated by in-frame TIS usage in Salmonella. Besides demonstrating that alternative TIS usage may generate proteoforms with different characteristics, such as differential localization and specialized function, quantitative aspects of conditional retapamulin-assisted ribosome profiling (Ribo-RET) translation initiation maps offer unprecedented insights into the relative utilization of annotated and alternative TIS, enabling the exploration of gene regulatory mechanisms that control TIS usage and, consequently, the translation of N-terminal proteoform pairs.

A Practical Guide to Small Protein Discovery and Characterization Using Mass Spectrometry

Small proteins of up to ∼50 amino acids are an abundant class of biomolecules across all domains of life. Yet due to the challenges inherent in their size, they are often missed in genome annotations, and are difficult to identify and characterize using standard experimental approaches. Consequently, we still know few small proteins even in well-studied prokaryotic model organisms. Mass spectrometry (MS) has great potential for the discovery, validation, and functional characterization of small proteins. However, standard MS approaches are poorly suited to the identification of both known and novel small proteins due to limitations at each step of a typical proteomics workflow, i.e., sample preparation, protease digestion, liquid chromatography, MS data acquisition, and data analysis. Here, we outline the major MS-based workflows and bioinformatic pipelines used for small protein discovery and validation. Special emphasis is placed on highlighting the adjustments required to improve detection and data quality for small proteins. We discuss both the unbiased detection of small proteins and the targeted analysis of small proteins of interest. Finally, we provide guidelines to prioritize novel small proteins, and an outlook on methods with particular potential to further improve comprehensive discovery and characterization of small proteins.

Joint engineering of SACE_Lrp and its target MarR enhances the biosynthesis and export of erythromycin in Saccharopolyspora erythraea

The Lrp and MarR families are two groups of transcriptional regulators widely distributed among prokaryotes. However, the hierarchical-regulatory relationship between the Lrp family and the MarR family remains unknown. Our previous study found that an Lrp (SACE_Lrp) from Saccharopolyspora erythraea indirectly repressed the biosynthesis of erythromycin. In this study, we characterized a novel MarR family protein (SACE_6745) from S. erythraea, which is controlled by SACE_Lrp and plays a direct regulatory role in erythromycin biosynthesis and export. SACE_Lrp directly regulated the expression of marR by specifically binding a precise site OM (5'-CTCCGGGAACCATT-3'). Gene disruption of marR increased the production of erythromycin by 45% in S. erythraea A226. We found that MarR has direct DNA-binding activity for the promoter regions of the erythromycin biosynthetic genes, as well as an ABC exporter SACE_2701-2702 which was genetically proved to be responsible for erythromycin efflux. Disruption of SACE_Lrp in industrial S. erythraea WB was an efficient strategy to enhance erythromycin production. Herein, we jointly engineered SACE_Lrp and its target MarR by deleting marR in WBΔSACE_Lrp, resulting in 20% increase in erythromycin yield in mutant WBΔLrpΔmarR compared to WBΔSACE_Lrp, and 39% to WB. Overall, our findings provide new insights into the hierarchical-regulatory relationship of Lrp and MarR proteins and new avenues for coordinating antibiotic biosynthesis and export by joint engineering regulators in actinomycetes. KEY POINTS: • The hierarchical-regulatory relationship between SACE_Lrp and MarR was identified. • MarR directly controlled the expression of erythromycin biosynthesis and export genes. • Joint engineering of SACE_Lrp-MarR regulatory element enhanced erythromycin production.

Bacterial riboproteogenomics: the era of N-terminal proteoform existence revealed

With the rapid increase in the number of sequenced prokaryotic genomes, relying on automated gene annotation became a necessity. Multiple lines of evidence, however, suggest that current bacterial genome annotations may contain inconsistencies and are incomplete, even for so-called well-annotated genomes. We here discuss underexplored sources of protein diversity and new methodologies for high-throughput genome reannotation. The expression of multiple molecular forms of proteins (proteoforms) from a single gene, particularly driven by alternative translation initiation, is gaining interest as a prominent contributor to bacterial protein diversity. In consequence, riboproteogenomic pipelines were proposed to comprehensively capture proteoform expression in prokaryotes by the complementary use of (positional) proteomics and the direct readout of translated genomic regions using ribosome profiling. To complement these discoveries, tailored strategies are required for the functional characterization of newly discovered bacterial proteoforms.

Characterization and engineering of the Lrp/AsnC family regulator SACE_5717 for erythromycin overproduction in Saccharopolyspora erythraea

In this work, we found that the Lrp/AsnC family protein SACE_5717 negatively regulated erythromycin biosynthesis in S. erythraea. Disruption of SACE_5717 led to a 27% improvement in the yield of erythromycin in S. erythraea A226. SACE_5717 directly repressed its own gene expression, as well as that of the adjacent gene SACE_5716 by binding to the target sequence 5'-GAACGTTCGCCGTCACGCC-3'. The predicted LysE superfamily protein SACE_5716 directly influenced the export of lysine, histidine, threonine and glycine in S. erythraea. Arginine, tyrosine and tryptophan were characterized as the effectors of SACE_5717 by weakening the binding affinity of SACE_5717. In the industrial S. erythraea WB strain, deletion of SACE_5717 (WBΔSACE_5717) increased erythromycin yield by 20%, and by 36% when SACE_5716 was overexpressed in WBΔSACE_5717 (WBΔSACE_5717/5716). In large-scale 5-L fermentation experiment, erythromycin yield in the engineered strain WBΔSACE_5717/5716 reached 4686 mg/L, a 41% enhancement over 3323 mg/L of the parent WB strain.

An integrative strategy to identify the entire protein coding potential of prokaryotic genomes by proteogenomics

Accurate annotation of all protein-coding sequences (CDSs) is an essential prerequisite to fully exploit the rapidly growing repertoire of completely sequenced prokaryotic genomes. However, large discrepancies among the number of CDSs annotated by different resources, missed functional short open reading frames (sORFs), and overprediction of spurious ORFs represent serious limitations. Our strategy toward accurate and complete genome annotation consolidates CDSs from multiple reference annotation resources, ab initio gene prediction algorithms and in silico ORFs (a modified six-frame translation considering alternative start codons) in an integrated proteogenomics database (iPtgxDB) that covers the entire protein-coding potential of a prokaryotic genome. By extending the PeptideClassifier concept of unambiguous peptides for prokaryotes, close to 95% of the identifiable peptides imply one distinct protein, largely simplifying downstream analysis. Searching a comprehensive Bartonella henselae proteomics data set against such an iPtgxDB allowed us to unambiguously identify novel ORFs uniquely predicted by each resource, including lipoproteins, differentially expressed and membrane-localized proteins, novel start sites and wrongly annotated pseudogenes. Most novelties were confirmed by targeted, parallel reaction monitoring mass spectrometry, including unique ORFs and single amino acid variations (SAAVs) identified in a re-sequenced laboratory strain that are not present in its reference genome. We demonstrate the general applicability of our strategy for genomes with varying GC content and distinct taxonomic origin. We release iPtgxDBs for B. henselae, Bradyrhizobium diazoefficiens and Escherichia coli and the software to generate both proteogenomics search databases and integrated annotation files that can be viewed in a genome browser for any prokaryote.

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.

Inactivation of SACE_3446, a TetR family transcriptional regulator, stimulates erythromycin production in Saccharopolyspora erythraea

Erythromycin A is a widely used antibiotic produced by Saccharopolyspora erythraea; however, its biosynthetic cluster lacks a regulatory gene, limiting the yield enhancement via regulation engineering of S. erythraea. Herein, six TetR family transcriptional regulators (TFRs) belonging to three genomic context types were individually inactivated in S. erythraea A226, and one of them, SACE_3446, was proved to play a negative role in regulating erythromycin biosynthesis. EMSA and qRT-PCR analysis revealed that SACE_3446 covering intact N-terminal DNA binding domain specifically bound to the promoter regions of erythromycin biosynthetic gene eryAI, the resistant gene ermE and the adjacent gene SACE_3447 (encoding a long-chain fatty-acid CoA ligase), and repressed their transcription. Furthermore, we explored the interaction relationships of SACE_3446 and previously identified TFRs (SACE_3986 and SACE_7301) associated with erythromycin production. Given demonstrated relatively independent regulation mode of SACE_3446 and SACE_3986 in erythromycin biosynthesis, we individually and concomitantly inactivated them in an industrial S. erythraea WB. Compared with WB, the WBΔ3446 and WBΔ3446Δ3986 mutants respectively displayed 36% and 65% yield enhancement of erythromycin A, following significantly elevated transcription of eryAI and ermE. When cultured in a 5 L fermentor, erythromycin A of WBΔ3446 and WBΔ3446Δ3986 successively reached 4095 mg/L and 4670 mg/L with 23% and 41% production improvement relative to WB. The strategy reported here will be useful to improve antibiotics production in other industrial actinomycete.

Systems Biology Approaches to Understand Natural Products Biosynthesis

Actinomycetes populate soils and aquatic sediments that impose biotic and abiotic challenges for their survival. As a result, actinomycetes metabolism and genomes have evolved to produce an overwhelming diversity of specialized molecules. Polyketides, non-ribosomal peptides, post-translationally modified peptides, lactams, and terpenes are well-known bioactive natural products with enormous industrial potential. Accessing such biological diversity has proven difficult due to the complex regulation of cellular metabolism in actinomycetes and to the sparse knowledge of their physiology. The past decade, however, has seen the development of omics technologies that have significantly contributed to our better understanding of their biology. Key observations have contributed toward a shift in the exploitation of actinomycete’s biology, such as using their full genomic potential, activating entire pathways through key metabolic elicitors and pathway engineering to improve biosynthesis. Here, we review recent efforts devoted to achieving enhanced discovery, activation, and manipulation of natural product biosynthetic pathways in model actinomycetes using genome-scale biological datasets.

Proteogenomics of rare taxonomic phyla: A prospective treasure trove of protein coding genes

Sustainable innovations in sequencing technologies have resulted in a torrent of microbial genome sequencing projects. However, the prokaryotic genomes sequenced so far are unequally distributed along their phylogenetic tree; few phyla contain the majority, the rest only a few representatives. Accurate genome annotation lags far behind genome sequencing. While automated computational prediction, aided by comparative genomics, remains a popular choice for genome annotation, substantial fraction of these annotations are erroneous. Proteogenomics utilizes protein level experimental observations to annotate protein coding genes on a genome wide scale. Benefits of proteogenomics include discovery and correction of gene annotations regardless of their phylogenetic conservation. This not only allows detection of common, conserved proteins but also the discovery of protein products of rare genes that may be horizontally transferred or taxonomy specific. Chances of encountering such genes are more in rare phyla that comprise a small number of complete genome sequences. We collated all bacterial and archaeal proteogenomic studies carried out to date and reviewed them in the context of genome sequencing projects. Here, we present a comprehensive list of microbial proteogenomic studies, their taxonomic distribution, and also urge for targeted proteogenomics of underexplored taxa to build an extensive reference of protein coding genes.

Proteogenomics in microbiology: taking the right turn at the junction of genomics and proteomics

High-accuracy and high-throughput proteomic methods have completely changed the way we can identify and characterize proteins. MS-based proteomics can now provide a unique supplement to genomic data and add a new level of information to the interpretation of genomic sequences. Proteomics-driven genome annotation has become especially relevant in microbiology where genomes are sequenced on a daily basis and limitations of an in silico driven annotation process are well recognized. In this review paper, we outline different strategies on how one can design a proteogenomic experiment, for example on genome-sequenced (synonymous proteogenomics) versus unsequenced organisms (ortho-proteogenomics) or with the aid of other "omic" data such as RNA-seq. We touch upon many challenges that are encountered during a typical proteogenomic study, mostly concerning bioinformatics methods and downstream data analysis, but also related to creation and use of sequence databases. A large list of proteogenomic case studies of different microorganisms is provided to illustrate the mapping of MS/MS-derived peptide spectra to genomic DNA sequences. These investigations have led to accurate determination of translational initiation sites, pointed out eventual read-throughs or programmed frameshifts, detected signal peptide processing or other protein maturation events, removed questionable annotation assignments, and provided evidence for predicted hypothetical proteins.

Control of chitin and N-acetylglucosamine utilization in Saccharopolyspora erythraea

Chitin degradation and subsequent N-acetylglucosamine (GlcNAc) catabolism is thought to be a common trait of a large majority of actinomycetes. Utilization of aminosugars had been poorly investigated outside the model strain Streptomyces coelicolor A3(2), and we examined here the genetic setting of the erythromycin producer Saccharopolyspora erythraea for GlcNAc and chitin utilization, as well as the transcriptional control thereof. Sacch. erythraea efficiently utilize GlcNAc most likely via the phosphotransferase system (PTS(GlcNAc)); however, this strain is not able to grow when chitin or N,N'-diacetylchitobiose [(GlcNAc)2] is the sole nutrient source, despite a predicted extensive chitinolytic system (chi genes). The inability of Sacch. erythraea to utilize chitin and (GlcNAc)2 is probably because of the loss of genes encoding the DasABC transporter for (GlcNAc)2 import, and genes for intracellular degradation of (GlcNAc)2 by β-N-acetylglucosaminidases. Transcription analyses revealed that in Sacch. erythraea all putative chi and GlcNAc utilization genes are repressed by DasR, whereas in Strep. coelicolor DasR displayed either activating or repressing functions whether it targets genes involved in the polymer degradation or genes for GlcNAc dimer and monomer utilization, respectively. A transcriptomic analysis further showed that GlcNAc not only activates the transcription of GlcNAc catabolism genes but also activates chi gene expression, as opposed to the previously reported GlcNAc-mediated catabolite repression in Strep. coelicolor. Finally, synteny exploration revealed an identical genetic background for chitin utilization in other rare actinomycetes, which suggests that screening procedures that used only the chitin-based protocol for selective isolation of antibiotic-producing actinomycetes could have missed the isolation of many industrially promising strains.