Proteogenomic Investigation of Strain Variation in Clinical Mycobacterium tuberculosis Isolates

Mycobacterium tuberculosis complex molecular networks and their regulation: Implications of strain heterogeneity on epigenetic diversity and transcriptome regulation

Tuberculosis has been a public health crisis since the 1900, which has caused the highest mortalities due to a single bacterial infection worldwide, that was recently further complicated by the Coronavirus disease 2019 pandemic. The causative agent of Tuberculosis, Mycobacterium tuberculosis, belongs to a genetically well-characterized family of strains known as the Mycobacterium tuberculosis complex, which has complicated progress made towards eradicating Tuberculosis due to pathogen-specific phenotypic differences in the members of this complex. Mycobacterium tuberculosis complex strains are genetically diverse human- and animal-adapted pathogens belonging to 7 lineages (Indo-Oceanic, East-Asian, East-African Indian, Euro-American, M. africanum West Africa 1, M. africanum West Africa 2 and Ethopia), respectively and the recently identified Lineage 8 and M. africanum Lineage 9. Genomic studies have revealed that Mycobacterium tuberculosis complex members are ∼99 % similar, however, due to selective pressure and adaptation to human host, they are prone to mutations that have resulted in development of drug resistance and phenotypic heterogeneity that impact strain virulence. Furthermore, members of the Mycobacterium tuberculosis complex have preferred geographic locations and possess unique phenotypic characteristics that is linked to their pathogenicity. Due to the recent advances in development next generation sequencing platforms, several studies have revealed epigenetic changes in genomic regions combined with "unique" gene regulatory mechanisms through non-coding RNAs that are responsible for strain-specific behaviour on in vitro and in vivo infection models. The current review provides up to date epigenetic patterns, gene regulation through non-coding RNAs, together with implications of these mechanisms in down-stream proteome and metabolome, which may be responsible for "unique" responses to infection by members of the Mycobacterium tuberculosis complex. Understanding lineage-specific molecular mechanisms during infection may provide novel drug targets and disease control measures towards World Health organization END-TB strategy.

Synthesis and Structure-Activity Relationship of 2,6-Disubstituted Thiosemicarbazone Derivatives of Pyridine as Potential Antituberculosis Agents

In this study, six new 2,6-disubstituted thiosemicarbazone derivatives of pyridine were synthesized (4−9), and their tuberculostatic activity was evaluated. All of them showed two- to eightfold higher activity (minimum inhibitory concentration (MIC) 0.5−4 µg/mL) against the resistant strain compared with the reference drug. Compounds 5 and 7, which contained the most basic substituents—pyrrolidine and piperidine—in their structure, strongly inhibited the growth of the standard strain (MIC 2 µg/mL). Furthermore, the same derivatives exhibited activity comparable to that of the reference drugs against some types of Gram-positive bacteria (MIC 0.49 µg/mL) and showed no cytotoxicity (IC50 > 50 µg/mL) in HaCaT cells. The zwitterionic structure of each compound was determined using X-ray crystallography. Absorption, distribution, metabolism, and excretion analyses showed that all compounds are good drug candidates. Thus, compounds 5 and 7 were identified as leading structures for further research on antituberculosis drugs with extended effects.

Types and functions of heterogeneity in mycobacteria

The remarkable ability of Mycobacterium tuberculosis to survive attacks from the host immune response and drug treatment is due to the resilience of a few bacilli rather than a result of survival of the entire population. Maintenance of mycobacterial subpopulations with distinct phenotypic characteristics is key for survival in the face of dynamic and variable stressors encountered during infection. Mycobacterial populations develop a wide range of phenotypes through an innate asymmetric growth pattern and adaptation to fluctuating microenvironments during infection that point to heterogeneity being a vital survival strategy. In this Review, we describe different types of mycobacterial heterogeneity and discuss how heterogeneity is generated and regulated in response to environmental cues. We discuss how this heterogeneity may have a key role in recording memory of their environment at both the single-cell level and the population level to give mycobacterial populations plasticity to withstand complex stressors.

Deep N-terminomics of Mycobacterium tuberculosis H37Rv extensively correct annotated encoding genes

Mycobacterium tuberculosis (MTB) is a severe causing agent of tuberculosis (TB). Although H37Rv, the type strain of M. tuberculosis was sequenced in 1998, annotation errors of encoding genes have been frequently reported in hundreds of papers. This phenomenon is particularly severe at the 5' end of the genes. Here, we applied a TMPP [(N-Succinimidyloxycarbonylmethyl) tris (2,4,6-trimethoxyphenyl) phosphonium bromide] labeling combined with StageTip separating strategy on M. tuberculosis H37Rv to characterize the N-terminal start sites of its annotated encoding genes. Totally, 1047 proteins were identified with 2058 TMPP labeled N-terminal peptides from all the 2625 mass spectrometer (MS) sequenced proteins. Comparative genomics analysis allowed the re-annotation of 43 proteins' N-termini in H37Rv and 762 proteins in Mycobacteriaceae. All revised N-termini start sites were distributed in 5'-UTR of annotated genes due to over-annotation of previous N-terminal initiation codon, especially the ATG. In addition, we identified and verified a novel gene Rv1078A in +3 frame different from the annotated gene Rv1078 in +2 frame. Altogether, our findings contribute to the better understanding of N-terminal of H37Rv and other species from Mycobacteriaceae that can assist future studies on biological study.

Benzofuroxan Derivatives as Potent Agents against Multidrug-Resistant Mycobacterium tuberculosis

Tuberculosis (TB) is currently the leading cause of death related to infectious diseases worldwide, as reported by the World Health Organization. Moreover, the increasing number of multidrug-resistant tuberculosis (MDR-TB) cases has alarmed health agencies, warranting extensive efforts to discover novel drugs that are effective and also safe. In this study, 23 new compounds were synthesized and evaluated in vitro against the drug-resistant strains of M. tuberculosis. The compound 6-((3-fluoro-4-thiomorpholinophenyl)carbamoyl)benzo[c][1,2,5]oxadiazole 1-N-oxide (5 b) was particularly remarkable in this regard as it demonstrated MIC₉₀ values below 0.28 μM against all the MDR strains evaluated, thus suggesting that this compound might have a different mechanism of action. Benzofuroxans are an attractive new class of anti-TB agents, exemplified by compound 5 b, with excellent potency against the replicating and drug-resistant strains of M. tuberculosis.

Proteogenomic Approach for Mycobacterium tuberculosis Investigation

Recent advances in MS/MS technology have made it possible to use proteomic data to predict protein-coding sequences. This approach is called proteogenomics, and it allows to correctly translate start and stop sites and to reveal new open reading frames. Here, we focus on using proteogenomics to improve the annotation of Mycobacteriumtuberculosis strains. We also describe detail procedures of the extraction of proteins and their further preparation for LC-MS/MS analysis and outline the main steps of data analysis.

Identification of novel protein-coding sequences in Eucalyptus grandis plants by high-resolution mass spectrometry

Eucalyptus species are widely used in the forestry industry, and a significant increase in the number of sequences available in database repositories has been observed for these species. In proteomics, a protein is identified by correlating the theoretical fragmentation spectrum derived from genomic/transcriptomic data against the experimental fragmentation mass spectrum acquired from large-scale analysis of protein mixtures. Proteogenomics is an alternative approach that can identify novel proteins encoded by regions previously considered as non-coding. This study aimed to confidently identify and confirm the existence of previously unknown protein-coding sequences in the Eucalyptus grandis genome. To this end, we used a modified spectral correlation strategy and a dedicated de novo peptide sequencing pipeline. Upon the strategy used here, we confidently identified 41 novel peptide forms and six peptides containing at least one single amino acid substitution. The most representative genomic class of novel peptides was identified as originating from alternative reading frames. In contrast, no clear single amino acid substitution pattern was identified. Validation of the identifications was carried out using a parallel reaction monitoring approach that provided further mass spectrometry support for the existence of the novel peptide sequences. Data are available via ProteomeXchange with identifier PXD022110.

Molecular Basis of Pathogenesis of Coronaviruses: A Comparative Genomics Approach to Planetary Health to Prevent Zoonotic Outbreaks in the 21st Century

In the first quarter of the 21st century, we are already facing the third emergence of a coronavirus outbreak, the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) responsible for the coronavirus disease 2019 (COVID-19) pandemic. Comparative genomics can inform a deeper understanding of the pathogenesis of COVID-19. Previous strains of coronavirus, SARS-CoV, and Middle-East respiratory syndrome-coronavirus (MERS-CoV), have been known to cause acute lung injuries in humans. SARS-CoV-2 shares genetic similarity with SARS-CoV with some modification in the S protein leading to their enhanced binding affinity toward the angiotensin-converting enzyme 2 (ACE2) receptors of human lung cells. This expert review examines the features of all three coronaviruses through a conceptual lens of comparative genomics. In particular, the life cycle of SARS-CoV-2 that enables its survival within the host is highlighted. Susceptibility of humans to coronavirus outbreaks in the 21st century calls for comparisons of the transmission history, hosts, reservoirs, and fatality rates of these viruses so that evidence-based and effective planetary health interventions can be devised to prevent future zoonotic outbreaks. Comparative genomics offers new insights on putative and novel viral targets with an eye to both therapeutic innovation and prevention. We conclude the expert review by (1) articulating the lessons learned so far, whereas the research is still being actively sought after in the field, and (2) the challenges and prospects in deciphering the linkages among multiomics biological variability and COVID-19 pathogenesis.

The PRIDE database and related tools and resources in 2019: improving support for quantification data

The PRoteomics IDEntifications (PRIDE) database (https://www.ebi.ac.uk/pride/) is the world’s largest data repository of mass spectrometry-based proteomics data, and is one of the founding members of the global ProteomeXchange (PX) consortium. In this manuscript, we summarize the developments in PRIDE resources and related tools since the previous update manuscript was published in Nucleic Acids Research in 2016. In the last 3 years, public data sharing through PRIDE (as part of PX) has definitely become the norm in the field. In parallel, data re-use of public proteomics data has increased enormously, with multiple applications. We first describe the new architecture of PRIDE Archive, the archival component of PRIDE. PRIDE Archive and the related data submission framework have been further developed to support the increase in submitted data volumes and additional data types. A new scalable and fault tolerant storage backend, Application Programming Interface and web interface have been implemented, as a part of an ongoing process. Additionally, we emphasize the improved support for quantitative proteomics data through the mzTab format. At last, we outline key statistics on the current data contents and volume of downloads, and how PRIDE data are starting to be disseminated to added-value resources including Ensembl, UniProt and Expression Atlas.

Interrogating Fractionation and Other Sources of Variability in Shotgun Proteomes Using Quality Metrics

The increasing amount of publicly available proteomics data creates opportunities for data scientists to investigate quality metrics in novel ways. QuaMeter IDFree is used to generate quality metrics from 665 RAW files and 97 WIFF files representing publicly available "shotgun" mass spectrometry datasets. These experiments are selected to represent Mycobacterium tuberculosis lysates, mouse MDSCs, and exosomes derived from human cell lines. Machine learning techniques are demonstrated to detect outliers within experiments and it is shown that quality metrics may be used to distinguish sources of variability among these experiments. In particular, the findings demonstrate that according to nested ANOVA performed on an SDS-PAGE shotgun principal component analysis, runs of fractions from the same gel regions cluster together rather than technical replicates, close temporal proximity, or even biological samples. This indicates that the individual fraction may have had a higher impact on the quality metrics than other factors. In addition, sample type, instrument type, mass analyzer, fragmentation technique, and digestion enzyme are identified as sources of variability. From a quality control perspective, the importance of study design and in particular, the run order, is illustrated in seeking ways to limit the impact of technical variability.

Identification of gene fusion events in Mycobacterium tuberculosis that encode chimeric proteins

Mycobacterium tuberculosis is a facultative intracellular pathogen responsible for causing tuberculosis. The harsh environment in which M. tuberculosis survives requires this pathogen to continuously adapt in order to maintain an evolutionary advantage. However, the apparent absence of horizontal gene transfer in M. tuberculosis imposes restrictions in the ways by which evolution can occur. Large-scale changes in the genome can be introduced through genome reduction, recombination events and structural variation. Here, we identify a functional chimeric protein in the ppe38-71 locus, the absence of which is known to have an impact on protein secretion and virulence. To examine whether this approach was used more often by this pathogen, we further develop software that detects potential gene fusion events from multigene deletions using whole genome sequencing data. With this software we could identify a number of other putative gene fusion events within the genomes of M. tuberculosis isolates. We were able to demonstrate the expression of one of these gene fusions at the protein level using mass spectrometry. Therefore, gene fusions may provide an additional means of evolution for M. tuberculosis in its natural environment whereby novel chimeric proteins and functions can arise.

Lineage-specific differences in lipid metabolism and its impact on clinical strains of Mycobacterium tuberculosis

Mycobacterium tuberculosis (M. tb) is the causative agent of TB and its incidences has been on the rise since 1993. Lipid metabolism is an imperative metabolic process, which grants M. tb the ability to utilize host-derived lipids as a secondary source of nutrition during infection. In addition to degrading host lipids, M. tb is proficient at using lipids, such as cholesterol, to facilitate its entry into macrophages. Mycolic acids, constituents of the mycobacterial cell wall, offer protection and aid in persistence of the bacterium. These are effectively synthesized using a complex fatty acid synthase system. Many pathogenesis studies have reported differences in lipid-metabolism of clinical strains of M. tb that belongs to diverse lineages of the Mycobacterium tuberculosis complex (MTBC). East-Asian and Euro-American lineages possess "unique" cell wall-associated lipids compared to the less transmissible Ethiopian lineage, which may offer these lineages a competitive advantage. Therefore, it is crucial to comprehend the complexities among the MTBC lineages with lipid metabolism and their impact on virulence, transmissibility and pathogenesis. Thus, this review provides an insight into lipid metabolism in various lineages of the MTBC and their impact on virulence and persistence during infection, as this may provide critical insight into developing novel therapeutics to combat TB.

Comprehensive Characterization of the Attenuated Double Auxotroph Mycobacterium tuberculosisΔleuDΔpanCD as an Alternative to H37Rv

Although currently available model organisms such as Mycobacterium smegmatis and Mycobacterium bovis Bacillus Calmette-Guérin (BCG) have significantly contributed to our understanding of tuberculosis (TB) biology, these models have limitations such as differences in genome size, growth rates and virulence. However, attenuated Mycobacterium tuberculosis strains may provide more representative, safer models to study M. tuberculosis biology. For example, the M. tuberculosis ΔleuDΔpanCD double auxotroph, has undergone rigorous in vitro and in vivo safety testing. Like other auxotrophic strains, this has subsequently been approved for use in biosafety level (BSL) 2 facilities. Auxotrophic strains have been assessed as models for drug-resistant M. tuberculosis and for studying latent TB. These offer the potential as safe and useful models, but it is important to understand how well these recapitulate salient features of non-attenuated M. tuberculosis. We therefore performed a comprehensive comparison of M. tuberculosis H37Rv and M. tuberculosisΔleuDΔpanCD. These strains demonstrated similar in vitro and intra-macrophage replication rates, similar responses to anti-TB agents and whole genome sequence conservation. Shotgun proteomics analysis suggested that M. tuberculosisΔleuDΔpanCD has a heightened stress response that leads to reduced bacterial replication during exposure to acid stress, which has been verified using a dual-fluorescent replication reporter assay. Importantly, infection of human peripheral blood mononuclear cells with the 2 strains elicited comparable cytokine production, demonstrating the suitability of M. tuberculosisΔleuDΔpanCD for immunological assays. We provide comprehensive evidence to support the judicious use of M. tuberculosisΔleuDΔpanCD as a safe and suitable model organism for M. tuberculosis research, without the need for a BSL3 facility.

Dual Analysis of Virus-Host Interactions: The Case of Ostreid herpesvirus 1 and the Cupped Oyster Crassostrea gigas

Dual analyses of the interactions between Ostreid herpesvirus 1 (OsHV-1) and the bivalve Crassostrea gigas during infection can unveil events critical to the onset and progression of this viral disease and can provide novel strategies for mitigating and preventing oyster mortality. Among the currently used “omics” technologies, dual transcriptomics (dual RNA-seq) coupled with the analysis of viral DNA in the host tissues has greatly advanced the knowledge of genes and pathways mostly contributing to host defense responses, expression profiles of annotated and unknown OsHV-1 open reading frames (ORFs), and viral genome variability. In addition to dual RNA-seq, proteomics and metabolomics analyses have the potential to add complementary information, needed to understand how a malacoherpesvirus can redirect and exploit the vital processes of its host. This review explores our current knowledge of “omics” technologies in the study of host-pathogen interactions and highlights relevant applications of these fields of expertise to the complex case of C gigas infections by OsHV-1, which currently threaten the mollusk production sector worldwide.

Proteomics in non-human primates: utilizing RNA-Seq data to improve protein identification by mass spectrometry in vervet monkeys

Background

Shotgun proteomics utilizes a database search strategy to compare detected mass spectra to a library of theoretical spectra derived from reference genome information. As such, the robustness of proteomics results is contingent upon the completeness and accuracy of the gene annotation in the reference genome. For animal models of disease where genomic annotation is incomplete, such as non-human primates, proteogenomic methods can improve the detection of proteins by incorporating transcriptional data from RNA-Seq to improve proteomics search databases used for peptide spectral matching. Customized search databases derived from RNA-Seq data are capable of identifying unannotated genetic and splice variants while simultaneously reducing the number of comparisons to only those transcripts actively expressed in the tissue.

Results

We collected RNA-Seq and proteomic data from 10 vervet monkey liver samples and used the RNA-Seq data to curate sample-specific search databases which were analyzed in the program Morpheus. We compared these results against those from a search database generated from the reference vervet genome. A total of 284 previously unannotated splice junctions were predicted by the RNA-Seq data, 92 of which were confirmed by peptide spectral matches. More than half (53/92) of these unannotated splice variants had orthologs in other non-human primates, suggesting that failure to match these peptides in the reference analyses likely arose from incomplete gene model information. The sample-specific databases also identified 101 unique peptides containing single amino acid substitutions which were missed by the reference database. Because the sample-specific searches were restricted to actively expressed transcripts, the search databases were smaller, more computationally efficient, and identified more peptides at the empirically derived 1 % false discovery rate.

Conclusion

Proteogenomic approaches are ideally suited to facilitate the discovery and annotation of proteins in less widely studies animal models such as non-human primates. We expect that these approaches will help to improve existing genome annotations of non-human primate species such as vervet.

Electronic supplementary material

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