Computational genomics-proteomics and Phylogeny analysis of twenty one mycobacterial genomes (Tuberculosis & non Tuberculosis strains)

Discovery of Species-Specific Proteotypic Peptides To Establish a Spectral Library Platform for Identification of Nontuberculosis Mycobacteria from Mass Spectrometry-Based Proteomics

Nontuberculous mycobacteria are opportunistic bacteria pulmonary and extra-pulmonary infections in humans that closely resemble Mycobacterium tuberculosis. Although genome sequencing strategies helped determine NTMs, a common assay for the detection of coinfection by multiple NTMs with M. tuberculosis in the primary attempt of diagnosis is still elusive. Such a lack of efficiency leads to delayed therapy, an inappropriate choice of drugs, drug resistance, disease complications, morbidity, and mortality. Although a high-resolution LC-MS/MS-based multiprotein panel assay can be developed due to its specificity and sensitivity, it needs a library of species-specific peptides as a platform. Toward this, we performed an analysis of proteomes of 9 NTM species with more than 20 million peptide spectrum matches gathered from 26 proteome data sets. Our metaproteomic analyses determined 48,172 species-specific proteotypic peptides across 9 NTMs. Notably, M. smegmatis (26,008), M. abscessus (12,442), M. vaccae (6487), M. fortuitum (1623), M. avium subsp. paratuberculosis (844), M. avium subsp. hominissuis (580), and M. marinum (112) displayed >100 species-specific proteotypic peptides. Finally, these peptides and corresponding spectra have been compiled into a spectral library, FASTA, and JSON formats for future reference and validation in clinical cohorts by the biomedical community for further translation.

Function-based classification of hazardous biological sequences: Demonstration of a new paradigm for biohazard assessments

Bioengineering applies analytical and engineering principles to identify functional biological building blocks for biotechnology applications. While these building blocks are leveraged to improve the human condition, the lack of simplistic, machine-readable definition of biohazards at the function level is creating a gap for biosafety practices. More specifically, traditional safety practices focus on the biohazards of known pathogens at the organism-level and may not accurately consider novel biodesigns with engineered functionalities at the genetic component-level. This gap is motivating the need for a paradigm shift from organism-centric procedures to function-centric biohazard identification and classification practices. To address this challenge, we present a novel methodology for classifying biohazards at the individual sequence level, which we then compiled to distinguish the biohazardous property of pathogenicity at the whole genome level. Our methodology is rooted in compilation of hazardous functions, defined as a set of sequences and associated metadata that describe coarse-level functions associated with pathogens (e.g., adherence, immune subversion). We demonstrate that the resulting database can be used to develop hazardous “fingerprints” based on the functional metadata categories. We verified that these hazardous functions are found at higher levels in pathogens compared to non-pathogens, and hierarchical clustering of the fingerprints can distinguish between these two groups. The methodology presented here defines the hazardous functions associated with bioengineering functional building blocks at the sequence level, which provide a foundational framework for classifying biological hazards at the organism level, thus leading to the improvement and standardization of current biosecurity and biosafety practices.

Mycobacterium abscessus: insights from a bioinformatic perspective

Mycobacterium abscessus is a nontuberculous mycobacterium, associated with broncho-pulmonary infections in individuals suffering from cystic fibrosis, bronchiectasis, and pulmonary diseases. The risk factors for transmission include biofilms, contaminated water resources, fomites, and infected individuals. M. abscessus is extensively resistant to antibiotics. To date, there is no vaccine and combination antibiotic therapy is followed. However, drug toxicities, low cure rates, and high cost of treatment make it imperfect. Over the last 20 years, bioinformatic studies on M. abscessus have advanced our understanding of the pathogen. This review integrates knowledge from the analysis of genomes, microbiomes, genomic variations, phylogeny, proteome, transcriptome, secretome, antibiotic resistance, and vaccine design to further our understanding. The utility of genome-based studies in comprehending disease progression, surveillance, tracing transmission routes, and epidemiological outbreaks on a global scale has been highlighted. Furthermore, this review underlined the importance of using computational methodologies for pinpointing factors responsible for pathogen survival and resistance. We reiterate the significance of interdisciplinary research to fight M. abscessus. In a nutshell, the outcome of computational studies can go a long way in creating novel therapeutic avenues to control M. abscessus mediated pulmonary infections.

The open pan-genome architecture and virulence landscape of Mycobacterium bovis

Animal tuberculosis (TB) is an emergent disease caused by Mycobacterium bovis, one of the animal-adapted ecotypes of the Mycobacterium tuberculosis complex (MTC). In this work, whole-genome comparative analyses of 70 M. bovis were performed to gain insights into the pan-genome architecture. The comparison across M. bovis predicted genome composition enabled clustering into the core- and accessory-genome components, with 2736 CDS for the former, while the accessory moiety included 3897 CDS, of which 2656 are restricted to one/two genomes only. These analyses predicted an open pan-genome architecture, with an average of 32 CDS added by each genome and show the diversification of discrete M. bovis subpopulations supported by both core- and accessory-genome components. The functional annotation of the pan-genome classified each CDS into one or several COG (Clusters of Orthologous Groups) categories, revealing ‘transcription’ (total average CDSs, n=258), ‘lipid metabolism and transport’ (n=242), ‘energy production and conversion’ (n=214) and ‘unknown function’ (n=876) as the most represented. The closer analysis of polymorphisms in virulence-related genes in a restrict group of M. bovis from a multi-host system enabled the identification of clade-monomorphic non-synonymous SNPs, illustrating clade-specific virulence landscapes and correlating with disease severity. This first comparative pan-genome study of a diverse collection of M. bovis encompassing all clonal complexes indicates a high percentage of accessory genes and denotes an open, dynamic non-conservative pan-genome structure, with high evolutionary potential, defying the canons of MTC biology. Furthermore, it shows that M. bovis can shape its virulence repertoire, either by acquisition and loss of genes or by SNP-based diversification, likely towards host immune evasion, adaptation and persistence.

Pan and Core Genome Analysis of 183 Mycobacterium tuberculosis Strains Revealed a High Inter-Species Diversity among the Human Adapted Strains

Tuberculosis (TB) is an airborne communicable disease with high morbidity and mortality rates, especially in developing countries. The causal agents of TB belong to the complex Mycobacterium tuberculosis (MTBc), which is composed of different human and animal TB associated species. Some animal associated species have zoonotic potential and add to the burden of TB management. The BCG ("Bacillus Calmette-Guérin") vaccine is widely used for the prevention against TB, but its use is limited in immunocompromised patients and animals due to the adverse effects and disseminated life-threatening complications. In this study, we aimed to carry out a comparative genome analysis between the human adapted species including BCG vaccine strains to identify and pinpoint the conserved genes related to the virulence across all the species, which could add a new value for vaccine development. For this purpose, the sequences of 183 Mycobacterium tuberculosis (MTB) strains were retrieved from the freely available WGS dataset at NCBI. The species included: 168 sensu stricto MTB species with other human MTB complex associated strains: M. tuberculosis var. africanum (3), M. tuberculosis var. bovis (2 draft genomes) and 10 BCG species, which enabled the analysis of core genome which contains the conserved genes and some virulence factor determinants. Further, a phylogenetic tree was constructed including the genomes of human (183); animals MTB adapted strains (6) and the environmental Mycobacterium strain "M. canettii". Our results showed that the core genome consists of 1166 conserved genes among these species, which represents a small portion of the pangenome (7036 genes). The remaining genes in the pangenome (5870) are accessory genes, adding a high inter-species diversity. Further, the core genome includes several virulence-associated genes and this could explain the rare infectiousness potential of some attenuated vaccine strains in some patients. This study reveals that low number of conserved genes in human adapted MTBc species and high inter-species diversity of the pan-genome could be considered for vaccine candidate development.

A Comparative Genomics Approach for Shortlisting Broad-Spectrum Drug Targets in Nontuberculous Mycobacteria

Many members of nontuberculous mycobacteria (NTM) are opportunistic pathogens causing several infections in animals. The incidence of NTM infections and emergence of drug-resistant NTM strains are rising worldwide, emphasizing the need to develop novel anti-NTM drugs. The present study is aimed to identify broad-spectrum drug targets in NTM using a comparative genomics approach. The study identified 537 core proteins in NTM of which 45 were pathogen specific and essential for the survival of pathogens. Furthermore, druggability analysis indicated that 15 were druggable among those 45 proteins. These 15 proteins, which were core proteins, pathogen-specific, essential, and druggable, were considered as potential broad-spectrum candidates. Based on their locations in cytoplasm and membrane, targets were classified as drug and vaccine targets. The identified 15 targets were different enzymes, carrier proteins, transcriptional regulator, two-component system protein, ribosomal, and binding proteins. The identified targets could further be utilized by researchers to design inhibitors for the discovery of antimicrobial agents.

Comparative Genomics of the Rhodococcus Genus Shows Wide Distribution of Biodegradation Traits

The genus Rhodococcus exhibits great potential for bioremediation applications due to its huge metabolic diversity, including biotransformation of aromatic and aliphatic compounds. Comparative genomic studies of this genus are limited to a small number of genomes, while the high number of sequenced strains to date could provide more information about the Rhodococcus diversity. Phylogenomic analysis of 327 Rhodococcus genomes and clustering of intergenomic distances identified 42 phylogenomic groups and 83 species-level clusters. Rarefaction models show that these numbers are likely to increase as new Rhodococcus strains are sequenced. The Rhodococcus genus possesses a small “hard” core genome consisting of 381 orthologous groups (OGs), while a “soft” core genome of 1253 OGs is reached with 99.16% of the genomes. Models of sequentially randomly added genomes show that a small number of genomes are enough to explain most of the shared diversity of the Rhodococcus strains, while the “open” pangenome and strain-specific genome evidence that the diversity of the genus will increase, as new genomes still add more OGs to the whole genomic set. Most rhodococci possess genes involved in the degradation of aliphatic and aromatic compounds, while short-chain alkane degradation is restricted to a certain number of groups, among which a specific particulate methane monooxygenase (pMMO) is only found in Rhodococcus sp. WAY2. The analysis of Rieske 2Fe-2S dioxygenases among rhodococci genomes revealed that most of these enzymes remain uncharacterized.

Comparative in-silico proteomic analysis discerns potential granuloma proteins of Yersinia pseudotuberculosis

Yersinia pseudotuberculosis is one of the three pathogenic species of the genus Yersinia. Most studies regarding pathogenesis of Y. pseudotuberculosis are based on the proteins related to Type III secretion system, which is a well-known primary virulence factor in pathogenic Gram-negative bacteria, including Y. pseudotuberculosis. Information related to the factors involved in Y. pseudotuberculosis granuloma formation is scarce. In the present study we have used a computational approach to identify proteins that might be potentially involved in formation of Y. pseudotuberculosis granuloma. A comparative proteome analysis and conserved orthologous protein identification was performed between two different genera of bacteria - Mycobacterium and Yersinia, their only common pathogenic trait being ability to form necrotizing granuloma. Comprehensive analysis of orthologous proteins was performed in proteomes of seven bacterial species. This included M. tuberculosis, M. bovis and M. avium paratuberculosis - the known granuloma forming Mycobacterium species, Y. pestis and Y. frederiksenii - the non-granuloma forming Yersinia species and, Y. enterocolitica - that forms micro-granuloma and, Y. pseudotuberculosis - a prominent granuloma forming Yersinia species. In silico proteome analysis indicated that seven proteins (UniProt id A0A0U1QT64, A0A0U1QTE0, A0A0U1QWK3, A0A0U1R1R0, A0A0U1R1Z2, A0A0U1R2S7, A7FMD4) might play some role in Y. pseudotuberculosis granuloma. Validation of the probable involvement of the seven proposed Y. pseudotuberculosis granuloma proteins was done using transcriptome data analysis and, by mapping on a composite protein-protein interaction map of experimentally proved M. tuberculosis granuloma proteins (RD1 locus proteins, ESAT-6 secretion system proteins and intra-macrophage secreted proteins). Though, additional experiments involving knocking out of each of these seven proteins are required to confirm their role in Y. pseudotuberculosis granuloma our study can serve as a basis for further studies on Y. pseudotuberculosis granuloma.

Retrospective Analysis of Archived Pyrazinamide Resistant Mycobacterium tuberculosis Complex Isolates from Uganda-Evidence of Interspecies Transmission

The contribution of Mycobacterium bovis to the proportion of tuberculosis cases in humans is unknown. A retrospective study was undertaken on archived Mycobacterium tuberculosis complex (MTBC) isolates from a reference laboratory in Uganda to identify the prevalence of human M. bovis infection. A total of 5676 isolates maintained in this repository were queried and 136 isolates were identified as pyrazinamide resistant, a hallmark phenotype of M. bovis. Of these, 1.5% (n = 2) isolates were confirmed as M. bovis by using regions of difference PCR analysis. The overall size of whole genome sequences (WGSs) of these two M. bovis isolates were ~4.272 Mb (M. bovis Bz_31150 isolated from a captive chimpanzee) and 4.17 Mb (M. bovis B2_7505 from a human patient), respectively. Alignment of these genomes against 15 MTBC genome sequences revealed 7248 single nucleotide polumorphisms (SNPs). Theses SNPs were used for phylogenetic analysis that indicated a strong relationship between M. bovis and the chimpanzee isolate (Bz_31150) while the other M. bovis genome from the human patient (B2_7505) analyzed did not cluster with any M. bovis or M. tuberculosis strains. WGS analysis also revealed multidrug resistance genotypes; these genomes revealed pncA mutations at positions H57D in Bz_31150 and B2_7505. Phenotypically, B2_7505 was an extensively drug-resistant strain and this was confirmed by the presence of mutations in the major resistance-associated proteins for all anti-tuberculosis (TB) drugs, including isoniazid (KatG (S315T) and InhA (S94A)), fluoroquinolones (S95T), streptomycin (rrs (R309C)), and rifampin (D435Y, a rare but disputed mutation in rpoB). The presence of these mutations exclusively in the human M. bovis isolate suggested that these occurred after transmission from cattle. Genome analysis in this study identified M. bovis in humans and great apes, suggesting possible transmission from domesticated ruminants in the area due to a dynamic and changing interface, which has created opportunity for exposure and transmission.

Cell-wall synthesis and ribosome maturation are co-regulated by an RNA switch in Mycobacterium tuberculosis

The success of Mycobacterium tuberculosis relies on the ability to switch between active growth and non-replicating persistence, associated with latent TB infection. Resuscitation promoting factors (Rpfs) are essential for the transition between these states. Rpf expression is tightly regulated as these enzymes are able to degrade the cell wall, and hence potentially lethal to the bacterium itself. We have identified a regulatory element in the 5' untranslated region (UTR) of rpfB. We demonstrate that this element is a transcriptionally regulated RNA switch/riboswitch candidate, which appears to be restricted to pathogenic mycobacteria, suggesting a role in virulence. We have used translation start site mapping to re-annotate the RpfB start codon and identified and validated a ribosome binding site that is likely to be targeted by an rpfB antisense RNA. Finally, we show that rpfB is co-transcribed with ksgA and ispE downstream. ksgA encodes a universally conserved methyltransferase involved in ribosome maturation and ispE encodes an essential kinase involved in cell wall synthesis. This arrangement implies co-regulation of resuscitation, cell wall synthesis and ribosome maturation via the RNA switch.

Decoding the similarities and differences among mycobacterial species

Mycobacteriaceae comprises pathogenic species such as Mycobacterium tuberculosis, M. leprae and M. abscessus, as well as non-pathogenic species, for example, M. smegmatis and M. thermoresistibile. Genome comparison and annotation studies provide insights into genome evolutionary relatedness, identify unique and pathogenicity-related genes in each species, and explore new targets that could be used for developing new diagnostics and therapeutics. Here, we present a comparative analysis of ten-mycobacterial genomes with the objective of identifying similarities and differences between pathogenic and non-pathogenic species. We identified 1080 core orthologous clusters that were enriched in proteins involved in amino acid and purine/pyrimidine biosynthetic pathways, DNA-related processes (replication, transcription, recombination and repair), RNA-methylation and modification, and cell-wall polysaccharide biosynthetic pathways. For their pathogenicity and survival in the host cell, pathogenic species have gained specific sets of genes involved in repair and protection of their genomic DNA. M. leprae is of special interest owing to its smallest genome (1600 genes and ~1300 psuedogenes), yet poor genome annotation. More than 75% of the pseudogenes were found to have a functional ortholog in the other mycobacterial genomes and belong to protein families such as transferases, oxidoreductases and hydrolases.

Members of the Mycobacteriaceae family, which are known to adapt to different environmental niches, comprise bacterial species with varied genome sizes. They are unique in their cell-wall composition, which is remarkably thick and lipid-rich as compared to other bacteria. We performed a comparative analysis at the proteome level for ten mycobacterial species that differ in their pathogenicity, genome size and environmental niches. A total of 1080 orthologous clusters with representation from all ten species were obtained, and these were further examined for their domain annotations, domain architecture similarities and enriched GO terms. These core orthologous clusters are enriched in various biosynthetic pathways. The proteins that are specific to each of the ten species were also investigated for their GO functions. The M. leprae genome has a large number of pseudogenes and we searched for their functional orthologs in other mycobacterial species in order to understand the functions that are lost from the M. leprae genome. The proteins present exclusively in M. leprae genome were studied in more detail, in order to predict putative drug targets and diagnostic markers. These findings, which have implications in understanding evolution of mycobacterial genomes, identify species-specific proteins that have potential for use in developing new diagnostic tools and therapeutics.

Identification of Novel Seroreactive Antigens in Johne's Disease Cattle by Using the Mycobacterium tuberculosis Protein Array

Johne's disease, a chronic gastrointestinal inflammatory disease caused by Mycobacterium avium subspecies paratuberculosis, is endemic in dairy cattle and other ruminants worldwide and remains a challenge to diagnose using traditional serological methods. Given the close phylogenetic relationship between M. avium subsp. paratuberculosis and the human pathogen Mycobacterium tuberculosis, here, we applied a whole-proteome M. tuberculosis protein array to identify seroreactive and diagnostic M. avium subsp. paratuberculosis antigens. A genome-scale pairwise analysis of amino acid identity levels between orthologous proteins in M. avium subsp. paratuberculosis and M. tuberculosis showed an average of 62% identity, with more than half the orthologous proteins sharing >75% identity. Analysis of the M. tuberculosis protein array probed with sera from M. avium subsp. paratuberculosis-infected cattle showed antibody binding to 729 M. tuberculosis proteins, with 58% of them having ≥70% identity to M. avium subsp. paratuberculosis orthologs. The results showed that only 4 of the top 40 seroreactive M. tuberculosis antigens were orthologs of previously reported M. avium subsp. paratuberculosis antigens, revealing the existence of a large number of previously unrecognized candidate diagnostic antigens. Enzyme-linked immunosorbent assay (ELISA) testing of 20 M. avium subsp. paratuberculosis recombinant proteins, representing reactive and nonreactive M. tuberculosis orthologs, further confirmed that the M. tuberculosis array has utility as a screening tool for identifying candidate antigens for Johne's disease diagnostics. Additional ELISA testing of field serum samples collected from dairy herds around the United States revealed that MAP2942c had the strongest seroreactivity with Johne's disease-positive samples. Collectively, our studies have considerably expanded the number of candidate M. avium subsp. paratuberculosis proteins with potential utility in the next generation of rationally designed Johne's disease diagnostic assays.

Nanodiagnostics for tuberculosis detection

INTRODUCTION

Tuberculosis (TB) is a leading killer worldwide. End TB strategy aims at ending the TB epidemic by 2030. Early, accurate, and affordable diagnosis represents a cornerstone to achieve this goal. Innovative strategies for TB diagnostics have been introduced. However, the ideal assay is yet unavailable and conventional methods remain necessary for diagnosis. Unique properties of nanoparticles (NPs) have allowed their utilization in TB detection via targeting disease biomarkers. Area covered: Until now, around thirty-five TB NP-based assays have been partially or fully characterized. Accuracy, low-cost, and short time-to-result represent the common properties of proposed platforms. TB nanodiagnostics now encompass almost all clinical aspects of the disease including active TB, non-tuberculous mycobacteria, rifampicin resistant TB, TB/HIV co-infection, latent TB, and extra-pulmonary TB. This review summarizes state-of-the-art knowledge of TB nanodiagnostics for the last 10 years. Special consideration is given for fabrication concepts, detection strategies, and clinical performance using various clinical specimens. The potential of TB nanodiagnostics to fulfill the need for ideal MTB testing is assessed. Expert commentary: TB nanodiagnostics show promise to be ideal detection tools that can meet the rigorous demands to end the TB epidemic by 2030.

Comparative genome analyses of mycobacteria give better insights into their evolution

Mycobacteria a genus of Actinobacteria are widespread in nature ranging from soil-dwelling saprophytes to human and animal pathogens. The rate of growth has been a classifying factor for the Mycobacterium spp., dividing them into the rapid growers and the slow growers. Here we have performed a comparative genome study of mycobacterial species in order to get better understanding of their evolution, particularly to understand the distinction between the rapid and slow growers. Our study shows that the slow growers had generally gained and lost more genes compared to the rapid growers. The slow growers might haved eventually lost genes (LivFGMH operon, shaACDEFG genes and MspA porin) that could contribute to the slow growth rate of the slow growers. The genes gained and lost in mycobacteria had eventually helped these bacteria to adapt to different environments and have led to the evolution of the present day rapid and slow growers. Our results also show high number of Mycobacterium abscessus specific genes (811 genes) and some of them are associated with the known bacterial quorum sensing genes that might be important for Mycobacterium abscessus to adapt and survive in variety of unfavorable environments. Mycobacterium abscessus also does not contains genes involved in the bacterial defense system and together with the quorum sensing genes may have contributed to the high gene gain rate of Mycobacterium abscessus.

Comparative pan genome analysis of oral Prevotella species implicated in periodontitis

Prevotella is part of the oral bacterial community implicated in periodontitis. Pan genome analyses of eight oral Prevotella species, P. dentalis, P. enoeca, P. fusca, P. melaninogenica, P. denticola, P. intermedia 17, P. intermedia 17-2 and P. sp. oral taxon 299 are presented in this study. Analysis of the Prevotella pan genome revealed features such as secretion systems, resistance to oxidative stress and clustered regularly interspaced short palindromic repeat (CRISPR)-Cas systems that enable the bacteria to adapt to the oral environment. We identified the presence of type VI secretion system (T6SS) in P. fusca and P. intermedia strains. For some VgrG and Hcp proteins which were not part of the core T6SS loci, we used gene neighborhood analysis and identified putative effector proteins and putative polyimmunity loci in P. fusca and polymorphic toxin systems in P. intermedia strains. Earlier studies have identified the presence of Por secretion system (PorSS) in P. gingivalis, P. melaninogenica and P. intermedia. We noted the presence of their homologs in six other oral Prevotella studied here. We suggest that in Prevotella, PorSS is used to secrete cysteine proteases such as interpain and C-terminal domain containing proteins with a "Por_secre_tail" domain. We identified subtype I-B CRISPR-Cas system in P. enoeca. Putative CRISPR-Cas system subtypes for 37 oral Prevotella and 30 non-oral Prevotella species were also predicted. Further, we performed a BLASTp search of the Prevotella proteins which are also conserved in the red-complex pathogens, against the human proteome to identify potential broad-spectrum drug targets. In summary, the use of a pan genome approach enabled identification of secretion systems and defense mechanisms in Prevotella that confer adaptation to the oral cavity.

Phylogenetic Analysis of Prevalent Tuberculosis and Non-Tuberculosis Mycobacteria in Isfahan, Iran, Based on a 360 bp Sequence of the rpoB Gene

Background

Taxonomic and phylogenetic studies of Mycobacterium species have been based around the 16sRNA gene for many years. However, due to the high strain similarity between species in the Mycobacterium genus (94.3% - 100%), defining a valid phylogenetic tree is difficult; consequently, its use in estimating the boundaries between species is limited. The sequence of the rpoB gene makes it an appropriate gene for phylogenetic analysis, especially in bacteria with limited variation.

Objectives

In the present study, a 360bp sequence of rpoB was used for precise classification of Mycobacterium strains isolated in Isfahan, Iran.

Materials and Methods

From February to October 2013, 57 clinical and environmental isolates were collected, subcultured, and identified by phenotypic methods. After DNA extraction, a 360bp fragment was PCR-amplified and sequenced. The phylogenetic tree was constructed based on consensus sequence data, using MEGA5 software.

Results

Slow and fast-growing groups of the Mycobacterium strains were clearly differentiated based on the constructed tree of 56 common Mycobacterium isolates. Each species with a unique title in the tree was identified; in total, 13 nods with a bootstrap value of over 50% were supported. Among the slow-growing group was Mycobacterium kansasii, with M. tuberculosis in a cluster with a bootstrap value of 98% and M. gordonae in another cluster with a bootstrap value of 90%. In the fast-growing group, one cluster with a bootstrap value of 89% was defined, including all fast-growing members present in this study.

Conclusions

The results suggest that only the application of the rpoB gene sequence is sufficient for taxonomic categorization and definition of a new Mycobacterium species, due to its high resolution power and proper variation in its sequence (85% - 100%); the resulting tree has high validity.

Validation of qPCR Methods for the Detection of Mycobacterium in New World Animal Reservoirs

Zoonotic pathogens that cause leprosy (Mycobacterium leprae) and tuberculosis (Mycobacterium tuberculosis complex, MTBC) continue to impact modern human populations. Therefore, methods able to survey mycobacterial infection in potential animal hosts are necessary for proper evaluation of human exposure threats. Here we tested for mycobacterial-specific single- and multi-copy loci using qPCR. In a trial study in which armadillos were artificially infected with M. leprae, these techniques were specific and sensitive to pathogen detection, while more traditional ELISAs were only specific. These assays were then employed in a case study to detect M. leprae as well as MTBC in wild marmosets. All marmosets were negative for M. leprae DNA, but 14 were positive for the mycobacterial rpoB gene assay. Targeted capture and sequencing of rpoB and other MTBC genes validated the presence of mycobacterial DNA in these samples and revealed that qPCR is useful for identifying mycobacterial-infected animal hosts.

Comparative whole-genome analysis of clinical isolates reveals characteristic architecture of Mycobacterium tuberculosis pangenome

The tubercle complex consists of closely related mycobacterium species which appear to be variants of a single species. Comparative genome analysis of different strains could provide useful clues and insights into the genetic diversity of the species. We integrated genome assemblies of 96 strains from Mycobacterium tuberculosis complex (MTBC), which included 8 Indian clinical isolates sequenced and assembled in this study, to understand its pangenome architecture. We predicted genes for all the 96 strains and clustered their respective CDSs into homologous gene clusters (HGCs) to reveal a hard-core, soft-core and accessory genome component of MTBC. The hard-core (HGCs shared amongst 100% of the strains) was comprised of 2,066 gene clusters whereas the soft-core (HGCs shared amongst at least 95% of the strains) comprised of 3,374 gene clusters. The change in the core and accessory genome components when observed as a function of their size revealed that MTBC has an open pangenome. We identified 74 HGCs that were absent from reference strains H37Rv and H37Ra but were present in most of clinical isolates. We report PCR validation on 9 candidate genes depicting 7 genes completely absent from H37Rv and H37Ra whereas 2 genes shared partial homology with them accounting to probable insertion and deletion events. The pangenome approach is a promising tool for studying strain specific genetic differences occurring within species. We also suggest that since selecting appropriate target genes for typing purposes requires the expected target gene be present in all isolates being typed, therefore estimating the core-component of the species becomes a subject of prime importance.

Comparative functional pan-genome analyses to build connections between genomic dynamics and phenotypic evolution in polycyclic aromatic hydrocarbon metabolism in the genus Mycobacterium

Background

The bacterial genus Mycobacterium is of great interest in the medical and biotechnological fields. Despite a flood of genome sequencing and functional genomics data, significant gaps in knowledge between genome and phenome seriously hinder efforts toward the treatment of mycobacterial diseases and practical biotechnological applications. In this study, we propose the use of systematic, comparative functional pan-genomic analysis to build connections between genomic dynamics and phenotypic evolution in polycyclic aromatic hydrocarbon (PAH) metabolism in the genus Mycobacterium.

Results

Phylogenetic, phenotypic, and genomic information for 27 completely genome-sequenced mycobacteria was systematically integrated to reconstruct a mycobacterial phenotype network (MPN) with a pan-genomic concept at a network level. In the MPN, mycobacterial phenotypes show typical scale-free relationships. PAH degradation is an isolated phenotype with the lowest connection degree, consistent with phylogenetic and environmental isolation of PAH degraders. A series of functional pan-genomic analyses provide conserved and unique types of genomic evidence for strong epistatic and pleiotropic impacts on evolutionary trajectories of the PAH-degrading phenotype. Under strong natural selection, the detailed gene gain/loss patterns from horizontal gene transfer (HGT)/deletion events hypothesize a plausible evolutionary path, an epistasis-based birth and pleiotropy-dependent death, for PAH metabolism in the genus Mycobacterium. This study generated a practical mycobacterial compendium of phenotypic and genomic changes, focusing on the PAH-degrading phenotype, with a pan-genomic perspective of the evolutionary events and the environmental challenges.

Conclusions

Our findings suggest that when selection acts on PAH metabolism, only a small fraction of possible trajectories is likely to be observed, owing mainly to a combination of the ambiguous phenotypic effects of PAHs and the corresponding pleiotropy- and epistasis-dependent evolutionary adaptation. Evolutionary constraints on the selection of trajectories, like those seen in PAH-degrading phenotypes, are likely to apply to the evolution of other phenotypes in the genus Mycobacterium.

Electronic supplementary material

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

Comparative functional genomics and the bovine macrophage response to strains of the mycobacterium genus

Mycobacterial infections are major causes of morbidity and mortality in cattle and are also potential zoonotic agents with implications for human health. Despite the implementation of comprehensive animal surveillance programs, many mycobacterial diseases have remained recalcitrant to eradication in several industrialized countries. Two major mycobacterial pathogens of cattle are Mycobacterium bovis and Mycobacterium avium subspecies paratuberculosis (MAP), the causative agents of bovine tuberculosis (BTB) and Johne's disease (JD), respectively. BTB is a chronic, granulomatous disease of the respiratory tract that is spread via aerosol transmission, while JD is a chronic granulomatous disease of the intestines that is transmitted via the fecal-oral route. Although these diseases exhibit differential tissue tropism and distinct complex etiologies, both M. bovis and MAP infect, reside, and replicate in host macrophages - the key host innate immune cell that encounters mycobacterial pathogens after initial exposure and mediates the subsequent immune response. The persistence of M. bovis and MAP in macrophages relies on a diverse series of immunomodulatory mechanisms, including the inhibition of phagosome maturation and apoptosis, generation of cytokine-induced necrosis enabling dissemination of infection through the host, local pathology, and ultimately shedding of the pathogen. Here, we review the bovine macrophage response to infection with M. bovis and MAP. In particular, we describe how recent advances in functional genomics are shedding light on the host macrophage-pathogen interactions that underlie different mycobacterial diseases. To illustrate this, we present new analyses of previously published bovine macrophage transcriptomics data following in vitro infection with virulent M. bovis, the attenuated vaccine strain M. bovis BCG, and MAP, and discuss our findings with respect to the differing etiologies of BTB and JD.

A novel molecular typing method of Mycobacteria based on DNA barcoding visualization

Different subtypes of Mycobacterium tuberculosis (MTB) may induce diverse severe human infections, and some of their symptoms are similar to other pathogenes, e.g. Nontuberculosis mycobacteria (NTM). So determination of mycobacterium subtypes facilitates the effective control of MTB infection and proliferation. This study exploits a novel DNA barcoding visualization method for molecular typing of 17 mycobacteria genomes published in the NCBI prokaryotic genome database. Three mycobacterium genes (Rv0279c, Rv3508 and Rv3514) from the PE/PPE family of MT Band were detected to best represent the inter-strain pathogenetic variations. An accurate and fast MTB substrain typing method was proposed based on the combination of the aforementioned three biomarker genes and the 16S rRNA gene. The protocol of establishing a bacterial substrain typing system used in this study may also be applied to the other pathogenes.