A Nonsynonymous SNP Catalog of Mycobacterium tuberculosis Virulence Genes and Its Use for Detecting New Potentially Virulent Sublineages

An overview of next generation sequencing strategies and genomics tools used for tuberculosis research

Tuberculosis (TB) is a grave public health concern and is considered the foremost contributor to human mortality resulting from infectious disease. Due to the stringent clonality and extremely restricted genomic diversity, conventional methods prove inefficient for in-depth exploration of minor genomic variations and the evolutionary dynamics operating in Mycobacterium tuberculosis (M.tb) populations. Until now, the majority of reviews have primarily focused on delineating the application of whole-genome sequencing (WGS) in predicting antibiotic resistant genes, surveillance of drug resistance strains, and M.tb lineage classifications. Despite the growing use of next generation sequencing (NGS) and WGS analysis in TB research, there are limited studies that provide a comprehensive summary of there role in studying macroevolution, minor genetic variations, assessing mixed TB infections, and tracking transmission networks at an individual level. This highlights the need for systematic effort to fully explore the potential of WGS and its associated tools in advancing our understanding of TB epidemiology and disease transmission. We delve into the recent bioinformatics pipelines and NGS strategies that leverage various genetic features and simultaneous exploration of host-pathogen protein expression profile to decipher the genetic heterogeneity and host-pathogen interaction dynamics of the M.tb infections. This review highlights the potential benefits and limitations of NGS and bioinformatics tools and discusses their role in TB detection and epidemiology. Overall, this review could be a valuable resource for researchers and clinicians interested in NGS-based approaches in TB research.

Large-scale Pan Genomic Analysis of Mycobacterium tuberculosis Reveals Key Insights Into Molecular Evolutionary Rate of Specific Processes and Functions

Mycobacterium tuberculosis (Mtb) is the causative agent of tuberculosis (TB), an infectious disease that is a major killer worldwide. Due to selection pressure caused by the use of antibacterial drugs, Mtb is characterised by mutational events that have given rise to multi drug resistant (MDR) and extensively drug resistant (XDR) phenotypes. The rate at which mutations occur is an important factor in the study of molecular evolution, and it helps understand gene evolution. Within the same species, different protein-coding genes evolve at different rates. To estimate the rates of molecular evolution of protein-coding genes, a commonly used parameter is the ratio dN/dS, where dN is the rate of non-synonymous substitutions and dS is the rate of synonymous substitutions. Here, we determined the estimated rates of molecular evolution of select biological processes and molecular functions across 264 strains of Mtb. We also investigated the molecular evolutionary rates of core genes of Mtb by computing the dN/dS values, and estimated the pan genome of the 264 strains of Mtb. Our results show that the cellular amino acid metabolic process and the kinase activity function evolve at a significantly higher rate, while the carbohydrate metabolic process evolves at a significantly lower rate for M. tuberculosis. These high rates of evolution correlate well with Mtb physiology and pathogenicity. We further propose that the core genome of M. tuberculosis likely experiences varying rates of molecular evolution which may drive an interplay between core genome and accessory genome during M. tuberculosis evolution.

Whole Genome Sequencing Highlights the Pathogenic Profile in Nocardia Keratitis

Purpose

Nocardia keratitis is a serious and sight-threatening condition. This study aims to reveal the virulence and antimicrobial resistance gene profile of Nocardia strains using whole genome sequencing.

Methods

Whole-genome sequencing was performed on 23 cornea-derived Nocardia strains. Together with genomic data from the respiratory tract and the environment, 141 genomes were then utilized for phylogenetic and pan-genome analyses, followed by virulence and antibiotic resistance analysis. The correlations between virulence genes and pathogenicity were experimentally validated, including the characteristics of Nocardia colonies and clinical and histopathological evaluations of Nocardia keratitis mice models.

Results

Whole-genome sequencing of 141 Nocardia strains revealed a mean of 220 virulence genes contributed to bacterial pathogenesis. The mce gene family analysis led to the categorization of strains from the cornea into groups A, B, and C. The colonies of group C had the largest diameter, height, and fastest growth rate. The size of corneal ulcers and the clinical scores showed a significant increase in mouse models induced by group C. The relative expression levels of pro-inflammatory cytokines (CD4, IFN-γ, IL-6Rα, and TNF-α) in the lesion area exhibited an increasing trend from group A to group C. Antibiotic resistance genes (ARGs) spanned nine distinct drug classes, four resistance mechanisms, and seven primary antimicrobial resistance gene families.

Conclusions

Whole genome sequencing highlights the pathogenic role of mce gene family in Nocardia keratitis. Its distribution pattern may contribute to the distinct characteristics of the growth of Nocardia colonies and the clinical severity of the mice models.

PCR Development for Analysis of Some Type II Toxin-Antitoxin Systems, relJK, mazEF3, and vapBC3 Genes, in Mycobacterium tuberculosis and Mycobacterium bovis

Toxin-Antitoxin (TA) systems are some small genetic modules in bacteria that play significant roles in resistance and tolerance development to antibiotics. Whole genome sequencing (WGS) is an effective method to analyze TA systems in pathogenic Mycobacteria. However, this study aimed to use a simple and inexpensive PCR-Sequencing approach to investigate the type II TA system. Using data from the WGS of Mycobacterium tuberculosis (M. tuberculosis) strain H37Rv and Mycobacterium bovis (M. bovis) strain BCG, primers specific to the relJK, mazEF3, and vapBC3 gene families were designed by Primer3 software. Following that, a total of 90 isolates were examined using the newly developed PCR assay, consisting of 64 M. tuberculosis and 26 M. bovis isolates, encompassing both 45 rifampin-sensitive and 45 rifampin-resistant strains. Finally, 28 isolates (including 14 rifampin-resistant isolates) were sent for sequencing, and their sequences were aligned and compared to the mentioned reference sequences. The amplicons size of mazEF3, relJK, and vapBC3 genes were 825, 875, and 934 bp, respectively. Furthermore, all tested isolates showed the specific amplicons for these TA families. To evaluate the specificity of the primers, PCR was performed on S. aureus and E.coli isolates. None of the examined samples had the desired amplicons. Therefore, the primers had acceptable specificity. The results indicated that the developed PCR-Sequencing approach can be used to effectively investigate certain types of TA systems. Considering high costs of WGS and difficulty in interpreting its results, such a simple and inexpensive method is beneficial in the evaluation of TA systems in Mycobacteria.

Polyketide synthases mutation in tuberculosis transmission revealed by whole genomic sequence, China, 2011-2019

Introduction: Tuberculosis (TB) is an infectious disease caused by a bacterium called Mycobacterium tuberculosis (Mtb). Previous studies have primarily focused on the transmissibility of multidrug-resistant (MDR) or extensively drug-resistant (XDR) Mtb. However, variations in virulence across Mtb lineages may also account for differences in transmissibility. In Mtb, polyketide synthase (PKS) genes encode large multifunctional proteins which have been shown to be major mycobacterial virulence factors. Therefore, this study aimed to identify the role of PKS mutations in TB transmission and assess its risk and characteristics. Methods: Whole genome sequences (WGSs) data from 3,204 Mtb isolates was collected from 2011 to 2019 in China. Whole genome single nucleotide polymorphism (SNP) profiles were used for phylogenetic tree analysis. Putative transmission clusters (≤10 SNPs) were identified. To identify the role of PKS mutations in TB transmission, we compared SNPs in the PKS gene region between "clustered isolates" and "non-clustered isolates" in different lineages. Results: Cluster-associated mutations in ppsA, pks12, and pks13 were identified among different lineage isolates. They were statistically significant among clustered strains, indicating that they may enhance the transmissibility of Mtb. Conclusion: Overall, this study provides new insights into the function of PKS and its localization in M. tuberculosis. The study found that ppsA, pks12, and pks13 may contribute to disease progression and higher transmission of certain strains. We also discussed the prospective use of mutant ppsA, pks12, and pks13 genes as drug targets.

High-Throughput Variant Detection Using a Color-Mixing Strategy

Many diseases are related to multiple genetic alterations along a single gene. Probing for highly multiple (>10) variants in a single quantitative PCR tube is not possible because of a limited number of fluorescence channels and one variant per channel, so many more tubes are needed. Herein, a novel color-mixing strategy that uses fluorescence combinations as digital color codes to probe multiple variants simultaneously was experimentally validated. The color-mixing strategy relies on a simple intratube assay that can probe for 15 variants as part of an intertube assay that can probe for an exponentially increased number of variants. This strategy is achieved by using multiplex double-stranded toehold probes modified with fluorophores and quenchers; the probes are designed to be quenched or luminous after binding to wild-type or variant templates. The color-mixing strategy was used to probe for 21 pathogenic variants in thalassemia and to distinguish between heterozygous and homozygous variants in six tubes, with a specificity of 99% and a sensitivity of 94%. To support tuberculosis diagnosis, the same strategy was applied to simultaneously probe in Mycobacterium tuberculosis for rifampicin-resistance mutations occurring within one 81-bp region and one 48-bp region in the rpoB gene, plus five isoniazid-resistance mutations in the inhA and katG genes.

Novel Screening System of Virulent Strains for the Establishment of a Mycobacterium avium Complex Lung Disease Mouse Model Using Whole-Genome Sequencing

The establishment of animal models reflecting human Mycobacterium avium complex (MAC) lung disease (LD) pathology has the potential to expand our understanding of the disease pathophysiology. However, inducing sustained infection in immunocompetent mice is difficult since MAC generally shows less virulence and higher genetic variability than M. tuberculosis. To overcome this hurdle, we developed a screening system for identifying virulent MAC strains using whole-genome sequencing (WGS). We obtained nine clinical strains from Mycobacterium avium complex lung disease (MAC-LD) patients and divided them into two groups to make the mixed strain inocula for infection. Intranasal infection with the strain mixture of both groups in BALB/c mice resulted in progressive infection and extensive granuloma formation in the lungs, suggesting the existence of highly pathogenic strains in each group. We hypothesized that the change in the abundance of strain-specific single-nucleotide variants (SNVs) reflects the change in bacterial number of each strain in infected lungs. Based on this hypothesis, we quantified individual strain-specific SNVs in bacterial DNA from infected lungs. Specific SNVs for four strains were detected, suggesting the pathogenicity of these four strains. Consistent with these results, individual infection with these four strains induced a high lung bacterial burden, forming extensive peribronchial granuloma, while the other strains showed a decreased lung bacterial burden. The current method combining mixed infection and WGS accurately identified virulent strains that induced sustained infection in mice. This method will contribute to the establishment of mouse models that reflect human MAC-LD and lead to antimycobacterial drug testing. IMPORTANCE To promote research on Mycobacterium avium complex (MAC) pathogenicity, animal models reflecting human progressive MAC lung disease (MAC-LD) are needed. Because there is high genetic and virulence diversity among clinical MAC strains, choosing a suitable strain is an important process for developing a mouse model. In this study, we developed a screening system for virulent strains in mice by combining mixed infection and whole-genome sequencing analysis. This approach is designed on the hypothesis that in vivo virulence of MAC strains can be examined simultaneously by comparing changes in the abundance of strain-specific single-nucleotide variants in the mouse lungs after infection with mixed strains. The identified strains were shown to induce high bacterial burdens and cause extensive peribronchial granuloma resembling the pulmonary pathology of human MAC-LD. The current method will help researchers develop mouse models that reflect human MAC-LD and will lead to further investigation of MAC pathogenicity.

A Murine Model of Mycobacterium kansasii Infection Reproducing Necrotic Lung Pathology Reveals Considerable Heterogeneity in Virulence of Clinical Isolates

Among non-tuberculous mycobacteria, Mycobacterium kansasii is one of the most pathogenic, able to cause pulmonary disease indistinguishable from tuberculosis in immunocompetent susceptible adults. The lack of animal models that reproduce human-like lung disease, associated with the necrotic lung pathology, impairs studies of M. kansasii virulence and pathogenicity. In this study, we examined the ability of the C57BL/6 mice, intratracheally infected with highly virulent M. kansasii strains, to produce a chronic infection and necrotic lung pathology. As a first approach, we evaluated ten M. kansasii strains isolated from Brazilian patients with pulmonary disease and the reference strain M. kansasii ATCC 12478 for virulence-associated features in macrophages infected in vitro; five of these strains differing in virulence were selected for in vivo analysis. Highly virulent isolates induced progressive lung disease in mice, forming large encapsulated caseous granulomas in later stages (120–150 days post-infection), while the low-virulent strain was cleared from the lungs by day 40. Two strains demonstrated increased virulence, causing premature death in the infected animals. These data demonstrate that C57BL/6 mice are an excellent candidate to investigate the virulence of M. kansasii isolates. We observed considerable heterogeneity in the virulence profile of these strains, in which the presence of highly virulent strains allowed us to establish a clinically relevant animal model. Comparing public genomic data between Brazilian isolates and isolates from other geographic regions worldwide demonstrated that at least some of the highly pathogenic strains isolated in Brazil display remarkable genomic similarities with the ATCC strain 12478 isolated in the United States 70 years ago (less than 100 SNPs of difference), as well as with some recent European clinical isolates. These data suggest that few pathogenic clones have been widely spread within M. kansasii population around the world.

Comparative genomics reveals an SNP potentially leading to phenotypic diversity of Salmonella enterica serovar Enteritidis

An SNP is a spontaneous genetic change having a potential to modify the functions of the original genes and to lead to phenotypic diversity of bacteria in nature. In this study, a phylogenetic analysis of Salmonella enterica serovar Enteritidis, a major food-borne pathogen, showed that eight strains of S. Enteritidis isolated in South Korea, including FORC_075 and FORC_078, have almost identical genome sequences. Interestingly, however, the abilities of FORC_075 to form biofilms and red, dry and rough (RDAR) colonies were significantly impaired, resulting in phenotypic differences among the eight strains. Comparative genomic analyses revealed that one of the non-synonymous SNPs unique to FORC_075 has occurred in envZ, which encodes a sensor kinase of the EnvZ/OmpR two-component system. The SNP in envZ leads to an amino acid change from Pro248 (CCG) in other strains including FORC_078 to Leu248 (CTG) in FORC_075. Allelic exchange of envZ between FORC_075 and FORC_078 identified that the SNP in envZ is responsible for the impaired biofilm- and RDAR colony-forming abilities of S. Enteritidis. Biochemical analyses demonstrated that the SNP in envZ significantly increases the phosphorylated status of OmpR in S. Enteritidis and alters the expression of the OmpR regulon. Phenotypic analyses further identified that the SNP in envZ decreases motility of S. Enteritidis but increases its adhesion and invasion to both human epithelial cells and murine macrophage cells. In addition to an enhancement of infectivity to the host cells, survival under acid stress was also elevated by the SNP in envZ. Together, these results suggest that the natural occurrence of the SNP in envZ could contribute to phenotypic diversity of S. Enteritidis, possibly improving its fitness and pathogenesis.

Exact mapping of Illumina blind spots in the Mycobacterium tuberculosis genome reveals platform-wide and workflow-specific biases

Whole-genome sequencing (WGS) is fundamental to Mycobacterium tuberculosis basic research and many clinical applications. Coverage across Illumina-sequenced M. tuberculosis genomes is known to vary with sequence context, but this bias is poorly characterized. Here, through a novel application of phylogenomics that distinguishes genuine coverage bias from deletions, we discern Illumina ‘blind spots’ in the M. tuberculosis reference genome for seven sequencing workflows. We find blind spots to be widespread, affecting 529 genes, and provide their exact coordinates, enabling salvage of unaffected regions. Fifty-seven pe/ppe genes (the primary families assumed to exhibit Illumina bias) lack blind spots entirely, while the remaining pe/ppe genes account for 55.1 % of blind spots. Surprisingly, we find coverage bias persists in homopolymers as short as 6 bp, shorter tracts than previously reported. While G+C-rich regions challenge all Illumina sequencing workflows, a modified Nextera library preparation that amplifies DNA with a high-fidelity polymerase markedly attenuates coverage bias in G+C-rich and homopolymeric sequences, expanding the ‘Illumina-sequenceable’ genome. Through these findings, and by defining workflow-specific exclusion criteria, we spotlight effective strategies for handling bias in M. tuberculosis Illumina WGS. This empirical analysis framework may be used to systematically evaluate coverage bias in other species using existing sequencing data.

The ATP-Binding Cassette (ABC) Transport Systems in Mycobacterium tuberculosis: Structure, Function, and Possible Targets for Therapeutics

Simple Summary

Mycobacterium tuberculosis is a bacterium of great medical importance because it causes tuberculosis, a disease that affects millions of people worldwide. Two important features are related to this bacterium: its ability to infect and survive inside the host, minimizing the immune response, and the burden of clinical isolates that are highly resistant to antibiotics treatment. These two phenomena are directly affected by cell envelope proteins, such as proteins from the ATP-Binding Cassette (ABC transporters) superfamily. In this review, we have compiled information on all the M. tuberculosis ABC transporters described so far, both from a functional and structural point of view, and show their relevance for the bacillus and the potential targets for studies aiming to control the microorganism and structural features.

Abstract

Mycobacterium tuberculosis is the etiological agent of tuberculosis (TB), a disease that affects millions of people in the world and that is associated with several human diseases. The bacillus is highly adapted to infect and survive inside the host, mainly because of its cellular envelope plasticity, which can be modulated to adapt to an unfriendly host environment; to manipulate the host immune response; and to resist therapeutic treatment, increasing in this way the drug resistance of TB. The superfamily of ATP-Binding Cassette (ABC) transporters are integral membrane proteins that include both importers and exporters. Both types share a similar structural organization, yet only importers have a periplasmic substrate-binding domain, which is essential for substrate uptake and transport. ABC transporter-type importers play an important role in the bacillus physiology through the transport of several substrates that will interfere with nutrition, pathogenesis, and virulence. Equally relevant, exporters have been involved in cell detoxification, nutrient recycling, and antibiotics and drug efflux, largely affecting the survival and development of multiple drug-resistant strains. Here, we review known ABC transporters from M. tuberculosis, with particular focus on the diversity of their structural features and relevance in infection and drug resistance.

CAPRIB: a user-friendly tool to study amino acid changes and selection for the exploration of intra-genus evolution

BACKGROUND

The evolution of bacteria is shaped by different mechanisms such as mutation, gene deletion, duplication, or insertion of foreign DNA among others. These genetic changes can accumulate in the descendants as a result of natural selection. Using phylogeny and genome comparisons, evolutionary paths can be somehow retraced, with recent events being much easier to detect than older ones. For this reason, multiple tools are available to study the evolutionary events within genomes of single species, such as gene composition alterations, or subtler mutations such as SNPs. However, these tools are generally designed to compare similar genomes and require advanced skills in bioinformatics. We present CAPRIB, a unique tool developed in Java that allows to determine the amino acid changes, at the genus level, that correlate with phenotypic differences between two groups of organisms.

RESULTS

CAPRIB has a user-friendly graphical interface and uses databases in SQL, making it easy to compare several genomes without the need for programming or thorough knowledge in bioinformatics. This intuitive software narrows down a list of amino acid changes that are concomitant with a given phenotypic divergence at the genus scale. Each permutation found by our software is associated with two already described statistical values that indicate its potential impact on the protein's function, helping the user decide which promising candidates to further investigate. We show that CAPRIB is able to detect already known mutations and uncovers many more, and that this tool can be used to question molecular phylogeny. Finally, we exemplify the utility of CAPRIB by pinpointing amino acid changes that coincided with the emergence of slow-growing mycobacteria from their fast-growing counterparts. The software is freely available at https://github.com/BactSymEvol/Caprib .

CONCLUSIONS

CAPRIB is a new bioinformatics software aiming to make genus-scale comparisons accessible to all. With its intuitive graphical interface, this tool identifies key amino acid changes concomitant with a phenotypic divergence. By comparing fast and slow-growing mycobacteria, we shed light on evolutionary hotspots, such as the cytokinin pathway, that are interesting candidates for further experimentations.

The whole genome sequence data analyses of a Mycobacterium tuberculosis strain SBH321 isolated in Sabah, Malaysia, belongs to Ural family of Lineage 4

In 2019, 10 million new cases of tuberculosis have been reported worldwide. Our data reports genetic analyses of a Mycobacterium tuberculosis strain SBH321 isolated from a 31-year-old female with pulmonary tuberculosis. The genomic DNA of the strain was extracted from pure culture and subjected to sequencing using Illumina platform. M. tuberculosis strain SBH321 consists of 4,374,895 bp with G+C content of 65.59%. The comparative analysis by SNP-based phylogenetic analysis using maximum-likelihood method showed that our strain belonging to sublineage of the Ural family of Europe-America-Africa lineage (Lineage 4) and clustered with M. tuberculosis strain OFXR-4 from Taiwan. The whole genome sequence is deposited at DDBJ/ENA/GenBank under the accession WCJH00000000 (SRR10230353).

Homoplastic single nucleotide polymorphisms contributed to phenotypic diversity in Mycobacterium tuberculosis

Homoplastic mutations are mutations independently occurring in different clades of an organism. The homoplastic changes may be a result of convergence evolution due to selective pressures. Reports on the analysis of homoplastic mutations in Mycobacterium tuberculosis have been limited. Here we characterized the distribution of homoplastic single nucleotide polymorphisms (SNPs) among genomes of 1,170 clinical M. tuberculosis isolates. They were present in all functional categories of genes, with pe/ppe gene family having the highest ratio of homoplastic SNPs compared to the total SNPs identified in the same functional category. Among the pe/ppe genes, the homoplastic SNPs were common in a relatively small number of homologous genes, including ppe18, the protein of which is a component of a promising candidate vaccine, M72/AS01E. The homoplastic SNPs in ppe18 were particularly common among M. tuberculosis Lineage 1 isolates, suggesting the need for caution in extrapolating the results of the vaccine trial to the population where L1 is endemic in Asia. As expected, homoplastic SNPs strongly associated with drug resistance. Most of these mutations are already well known. However, a number of novel mutations associated with streptomycin resistance were identified, which warrants further investigation. A SNP in the intergenic region upstream of Rv0079 (DATIN) was experimentally shown to increase transcriptional activity of the downstream gene, suggesting that intergenic homoplastic SNPs should have effects on the physiology of the bacterial cells. Our study highlights the potential of homoplastic mutations to produce phenotypic changes. Under selective pressure and during interaction with the host, homoplastic mutations may confer advantages to M. tuberculosis and deserve further characterization.

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.

Comparative genomics and association analysis identifies virulence genes of Cercospora sojina in soybean

BACKGROUND

Recently, a new strain of Cercospora sojina (Race15) has been identified, which has caused the breakdown of resistance in most soybean cultivars in China. Despite this serious yield reduction, little is known about why this strain is more virulent than others. Therefore, we sequenced the Race15 genome and compared it to the Race1 genome sequence, as its virulence is significantly lower. We then re-sequenced 30 isolates of C. sojina from different regions to identifying differential virulence genes using genome-wide association analysis (GWAS).

RESULTS

The 40.12-Mb Race15 genome encodes 12,607 predicated genes and contains large numbers of gene clusters that have annotations in 11 different common databases. Comparative genomics revealed that although these two genomes had a large number of homologous genes, their genome structures have evolved to introduce 245 specific genes. The most important 5 candidate virulence genes were located on Contig 3 and Contig 1 and were mainly related to the regulation of metabolic mechanisms and the biosynthesis of bioactive metabolites, thereby putatively affecting fungi self-toxicity and reducing host resistance. Our study provides insight into the genomic basis of C. sojina pathogenicity and its infection mechanism, enabling future studies of this disease.

CONCLUSIONS

Via GWAS, we identified five candidate genes using three different methods, and these candidate genes are speculated to be related to metabolic mechanisms and the biosynthesis of bioactive metabolites. Meanwhile, Race15 specific genes may be linked with high virulence. The genes highly prevalent in virulent isolates should also be proposed as candidates, even though they were not found in our SNP analysis. Future work should focus on using a larger sample size to confirm and refine candidate gene identifications and should study the functional roles of these candidates, in order to investigate their potential roles in C. sojina pathogenicity.

Genomic Analysis of Mycobacterium tuberculosis Isolates and Construction of a Beijing Lineage Reference Genome

Tuberculosis (TB), an infectious disease caused by Mycobacterium tuberculosis, kills over 1 million people worldwide annually. Development of drug resistance (DR) in the pathogen is a major challenge for TB control. We conducted whole-genome analysis of seven Taiwan M. tuberculosis isolates: One drug susceptible (DS) and five DR Beijing lineage isolates and one DR Euro-American lineage isolate. Developing a new method for DR mutation identification and applying it to the next-generation sequencing (NGS) data from the 6 Beijing lineage isolates, we identified 13 known and 6 candidate DR mutations and provided experimental support for 4 of them. We assembled the genomes of one DS and two DR Beijing lineage isolates and the Euro-American lineage isolate using NGS data. Moreover, using both PacBio and NGS sequencing data, we obtained a high-quality assembly of an extensive DR Beijing lineage isolate. Comparative analysis of these five newly assembled genomes and two published complete genomes revealed a large number of genetic changes, including gene gains and losses, indels and translocations, suggesting rapid evolution of M. tuberculosis. We found the MazEF toxin-antitoxin system in all the seven isolates studied and several interesting mutations in MazEF proteins. Finally, we used the four assembled Beijing lineage genomes to construct a high-quality Beijing lineage reference genome that is DS and contains all the genes in the four genomes. It contains 212 genes not found in the standard reference H37Rv, which is Euro-American. It is therefore a better reference than H37Rv for the Beijing lineage, the predominant lineage in Asia.

Draft Genome Sequences of Mycobacterium tuberculosis Clinical Isolates from the Ural Region of Russia That Carry the pks15/1 Gene

Here, we report the draft genome sequences of 15 Mycobacterium tuberculosis isolates of the Beijing-B0/W-148 sublineage that carry a 7-bp insertion within the pks15 gene, which leads to the synthesis of Pks15/1 fusion protein. Pks15/1 is involved in phenolglycolipid synthesis and biofilm formation, thus potentially contributing to the B0/W-148 lineage's enhanced virulence and drug resistance.

Sequencing and Analysis of Three Mycobacterium tuberculosis Genomes of the B0/N-90 Sublineage

We report the draft genome sequences of three Mycobacterium tuberculosis isolates belonging to the B0/N-90 sublineage, EKB34, EKB53, and EKB79. The B0/N-90 sublineage belongs to the prevalent (in Russia) and highly virulent Beijing-B0/W148 sublineage. Isolates EKB34 and EKB79 were obtained from people with immune deficiency.

Whole genome sequencing data and analysis of a rifampicin-resistant Mycobacterium tuberculosis strain SBH162 from Sabah, Malaysia

A Mycobacterium tuberculosis strain SBH162 was isolated from a 49-year-old male with pulmonary tuberculosis. GeneXpert MDR/RIF identified the strain as rifampicin-resistant M. tuberculosis. The whole genome sequencing was performed using Illumina HiSeq 4000 system to further investigate and verify the mutation sites of the strain through genetic analyses namely variant calling using bioinformatics tools. The de novo assembly of genome generated 100 contigs with N50 of 156,381bp. The whole genome size was 4,343,911 bp with G + C content of 65.58% and consisted of 4,306 predicted genes. The mutation site, S450L, for rifampicin resistance was detected in the rpoB gene. Based on the phylogenetic analysis using the Maximum Likelihood method, the strain was identified as belonging to the Europe America Africa lineage (Lineage 4). The genome dataset has been deposited at DDBJ/ENA/GenBank under the accession number SMOE00000000.

Evaluation of a gene-by-gene approach for prospective whole-genome sequencing-based surveillance of multidrug resistant Mycobacterium tuberculosis

Whole-genome sequencing (WGS) offers unprecedented resolution for tracking Mycobacterium tuberculosis transmission and antibiotic-resistance spread. Still, the establishment of standardized WGS-based pipelines and the definition of epidemiological clusters based on genetic relatedness are under discussion. We aimed to implement a dynamic gene-by-gene approach, fully relying on freely available software, for prospective WGS-based tuberculosis surveillance, demonstrating its application for detecting transmission chains by retrospectively analysing all M/XDR strains isolated in 2013-2017 in Portugal. We observed a good correlation between genetic relatedness and epidemiological links, with strongly epilinked clusters displaying mean pairwise allele differences (AD) always below 0.3% (ratio of mean AD over the total number of shared loci between same-cluster strains). This data parallels the genetic distances acquired by the core-SNV analysis, while providing higher resolution and epidemiological concordance than MIRU-VNTR genotyping. The dynamic analysis of strain sub-sets (i.e., increasing the number of shared loci within each sub-set) also strengthens the confidence in detecting epilinked clusters. This gene-by-gene strategy also offers several practical benefits (e.g., reliance on freely-available software, scalability and low computational requirements) that further consolidated its suitability for a timely and robust prospective WGS-based laboratory surveillance of M/XDR-TB cases.

A Major Mycobacterium tuberculosis outbreak caused by one specific genotype in a low-incidence country: Exploring gene profile virulence explanations

Denmark, a tuberculosis low burden country, still experiences significant active Mycobacterium tuberculosis (Mtb) transmission, especially with one specific genotype named Cluster 2/1112-15 (C2), the most prevalent lineage in Scandinavia. In addition to environmental factors, antibiotic resistance, and human genetics, there is increasing evidence that Mtb strain variation plays a role for the outcome of infection and disease. In this study, we explore the reasons for the success of the C2 genotype by analysing strain specific polymorphisms identified through whole genome sequencing of all C2 isolates identified in Denmark between 1992 and 2014 (n = 952), and the demographic distribution of C2. Of 234 non-synonymous (NS) monomorphic SNPs found in C2 in comparison with Mtb reference strain H37Rv, 23 were in genes previously reported to be involved in Mtb virulence. Of these 23 SNPs, three were specific for C2 including a NS mutation in a gene associated with hyper-virulence. We show that the genotype is readily transmitted to different ethnicities and is also found outside Denmark. Our data suggest that strain specific virulence factor variations are important for the success of the C2 genotype. These factors, likely in combination with poor TB control, seem to be the main drivers of C2 success.

Single nucleotide polymorphisms of Beijing lineage Mycobacterium tuberculosis toxin-antitoxin system genes: Their role in the changes of protein activity and evolution

The article investigates SNP in genes of toxin-antitoxin systems type II in Mycobacterium tuberculosis Beijing lineage strains and their possible role in the development and formation of new sublineages. We established the catalog of SNPs in 142 TA systems genes in 1349 sequenced genomes of the M. tuberculosis Beijing lineage. Based on the catalog, 15 new sublineages were identified as part of Beijing lineages by non-synonymous SNP in 21 genes of TA systems. We discovered three toxin genes with mutations specific for epidemiologically dangerous sublineages Beijing-modern (vapC37 A46G, vapC38 T143C) and Beijing-B0/W148 (vapC12 A95G). We proved the functional significance of these polymorphisms by cloning these genes wild-type and with marker mutations for the Beijing lineage vapC12 (A95G), vapC37 (A46G), vapC38 (T143C). In vitro study of their activities revealed effect of mutations on the RNase activity of toxin proteins. Mutations in vapC37 and vapC38 decreased toxin activity, and mutation in the vapC12 increased it. We cloned the toxin vapC37 gene of Mycobacterium smegmatis mc² 155 in both allelic variants: without mutation and with A46G mutation, specific for the Beijing-modern lineage. It was shown that this mutation leads to a loss of toxicity.

Mycobacterium tuberculosis and whole-genome sequencing: how close are we to unleashing its full potential?

BACKGROUND

Nearly two decades after completion of the genome sequence of Mycobacterium tuberculosis (MTB), and with the advent of next generation sequencing technologies (NGS), whole-genome sequencing (WGS) has been applied to a wide range of clinical scenarios. Starting in 2017, England is the first country in the world to pioneer its use on a national scale for the diagnosis of tuberculosis, detection of drug resistance, and typing of MTB.

AIMS

This narrative review critically analyses the current applications of WGS for MTB and explains how close we are to realizing its full potential as a diagnostic, epidemiologic, and research tool.

SOURCES

We searched for reports (both original articles and reviews) published in English up to 31 May 2017, with combinations of the following keywords: whole-genome sequencing, Mycobacterium, and tuberculosis. MEDLINE, Embase, and Scopus were used as search engines. We included articles that covered different aspects of whole-genome sequencing in relation to MTB.

CONTENT

This review focuses on three main themes: the role of WGS for the prediction of drug susceptibility, MTB outbreak investigation and genetic diversity, and research applications of NGS.

IMPLICATIONS

Many of the original expectations have been accomplished, and we believe that with its unprecedented sensitivity and power, WGS has the potential to address many unanswered questions in the near future. However, caution is still needed when interpreting WGS data as there are some important limitations to be aware of, from correct interpretation of drug susceptibilities to the bioinformatic support needed.