Bioinformatics tools and databases for whole genome sequence analysis of Mycobacterium tuberculosis

Analysis of Drug-Resistance Characteristics and Genetic Diversity of Multidrug-Resistant Tuberculosis Based on Whole-Genome Sequencing on the Hainan Island, China

Purpose

Given the high burden of Tuberculosis (TB) in China, the prevalence of multidrug-resistant tuberculosis (MDR-TB) is significant. Whole-genome sequencing (WGS) of Mycobacterium tuberculosis (MTB) enables the identification of lineages, drug-resistant mutations, and transmission patterns, offering valuable insights for TB control, clinical diagnosis, and treatment.

Methods

We collected 202 MDR-MTB strains from 3519 suspected pulmonary TB patients treated at The Second Affiliated Hospital of Hainan Medical University between July 2019 and June 2021. Proportional drug-susceptibility testing was performed using 8 common anti-tuberculosis drugs. Subsequently, the genotypic drug resistance and genetic characteristics were analyzed by the WGS.

Results

Lineages are identified by TB-profiler revealed 202 MDR-MTB strains, showcasing three predominant lineages, with lineage 2 being the most prevalent. Close genomic relatedness analysis and evidence of MTB transmission led to the formation of 15 clusters comprising 42 isolates, resulting in a clustering rate of 20.8%. Novelty, lineage 2.1 (non-Beijing) accounted for 27.2% of the MDR-MTB strains, which is rare in China and Neighboring countries. Regarding first-line anti-TB drugs, genes associated with rifampicin resistance, primarily the rpoB gene, were detected in 200 strains (99.0%). Genes conferring resistance to isoniazid, ethambutol, and streptomycin were identified in 191 (94.5%), 125 (61.9%), and 100 (49.5%) strains, respectively. Among the second-line drugs, 97 (48.0%) strains exhibited genes encoding resistance to fluoroquinolones. Comparing the results to phenotypic drug susceptibility-based testing, the sensitivity of WGS for detecting resistance to each of the six drugs (rifampicin, isoniazid, ethambutol, ofloxacin, kanamycin, capreomycin) was 90% or higher. With the exception of ethambutol, the specificity of WGS prediction for the remaining drugs exceeded 88%.

Conclusion

Our study provides crucial insights into genetic mutation types, genetic diversity, and transmission of MDR-MTB on Hainan Island, serving as a significant reference for MDR-MTB surveillance and clinical decision-making.

Role of Delta/Notch-like epidermal growth factor-related receptor in gastric cancer patients and cells and its clinical significance

Gastric cancer (GC) is a common digestive system malignancy. The aim of this study was to explore the role of Delta/Notch-like epidermal growth factor-related receptor (DNER) in GC patients and cells. Gene expression omnibus data base public databases were used to analyze the DNER expression in GC patient. A total of 30 cases of GC and adjacent tissue samples were retrospectively obtained to analyze the DNER expression. MTT assay was conducted to measure the cell viability. The apoptosis rate of GC cells was determined by flow cytometry. The migration and invasion were detected by transwell assay. Real-time polymerase chain reaction and western blot were performed to measure the DNER expression. Bioinformatics tools exhibited that DNER expression is significantly upregulated in the GC, which was also found in GC tissues and cells. The high levels of DNER were closely related the tumor size, sex and lymph node metastasis. Additionally, the survival rate of patients with high DNER expression is decreased. Furthermore, knockdown of DNER inhibits the proliferation, migration and invasion, and induces the apoptosis rate of the GC cells. DNER was upregulated in GC and knockdown of DNER inhibits the growth and metastasis of DNER. DNER may be a potential prognostic biomarker and therapeutic target of GC patients.

Structure-Aware Mycobacterium tuberculosis Functional Annotation Uncloaks Resistance, Metabolic, and Virulence Genes

Accurate and timely functional genome annotation is essential for translating basic pathogen research into clinically impactful advances. Here, through literature curation and structure-function inference, we systematically update the functional genome annotation of Mycobacterium tuberculosis virulent type strain H37Rv. First, we systematically curated annotations for 589 genes from 662 publications, including 282 gene products absent from leading databases. Second, we modeled 1,711 underannotated proteins and developed a semiautomated pipeline that captured shared function between 400 protein models and structural matches of known function on Protein Data Bank, including drug efflux proteins, metabolic enzymes, and virulence factors. In aggregate, these structure- and literature-derived annotations update 940/1,725 underannotated H37Rv genes and generate hundreds of functional hypotheses. Retrospectively applying the annotation to a recent whole-genome transposon mutant screen provided missing function for 48% (13/27) of underannotated genes altering antibiotic efficacy and 33% (23/69) required for persistence during mouse tuberculosis (TB) infection. Prospective application of the protein models enabled us to functionally interpret novel laboratory generated pyrazinamide (PZA)-resistant mutants of unknown function, which implicated the emerging coenzyme A depletion model of PZA action in the mutants’ PZA resistance. Our findings demonstrate the functional insight gained by integrating structural modeling and systematic literature curation, even for widely studied microorganisms. Functional annotations and protein structure models are available at https://tuberculosis.sdsu.edu/H37Rv in human- and machine-readable formats.

IMPORTANCEMycobacterium tuberculosis, the primary causative agent of tuberculosis, kills more humans than any other infectious bacterium. Yet 40% of its genome is functionally uncharacterized, leaving much about the genetic basis of its resistance to antibiotics, capacity to withstand host immunity, and basic metabolism yet undiscovered. Irregular literature curation for functional annotation contributes to this gap. We systematically curated functions from literature and structural similarity for over half of poorly characterized genes, expanding the functionally annotated Mycobacterium tuberculosis proteome. Applying this updated annotation to recent in vivo functional screens added functional information to dozens of clinically pertinent proteins described as having unknown function. Integrating the annotations with a prospective functional screen identified new mutants resistant to a first-line TB drug, supporting an emerging hypothesis for its mode of action. These improvements in functional interpretation of clinically informative studies underscore the translational value of this functional knowledge. Structure-derived annotations identify hundreds of high-confidence candidates for mechanisms of antibiotic resistance, virulence factors, and basic metabolism and other functions key in clinical and basic tuberculosis research. More broadly, they provide a systematic framework for improving prokaryotic reference annotations.

Trends in biological data integration for the selection of enzymes and transcription factors related to cellulose and hemicellulose degradation in fungi

Fungi are key players in biotechnological applications. Although several studies focusing on fungal diversity and genetics have been performed, many details of fungal biology remain unknown, including how cellulolytic enzymes are modulated within these organisms to allow changes in main plant cell wall compounds, cellulose and hemicellulose, and subsequent biomass conversion. With the advent and consolidation of DNA/RNA sequencing technology, different types of information can be generated at the genomic, structural and functional levels, including the gene expression profiles and regulatory mechanisms of these organisms, during degradation-induced conditions. This increase in data generation made rapid computational development necessary to deal with the large amounts of data generated. In this context, the origination of bioinformatics, a hybrid science integrating biological data with various techniques for information storage, distribution and analysis, was a fundamental step toward the current state-of-the-art in the postgenomic era. The possibility of integrating biological big data has facilitated exciting discoveries, including identifying novel mechanisms and more efficient enzymes, increasing yields, reducing costs and expanding opportunities in the bioprocess field. In this review, we summarize the current status and trends of the integration of different types of biological data through bioinformatics approaches for biological data analysis and enzyme selection.

Whole-genome analysis of drug-resistant Mycobacterium tuberculosis reveals novel mutations associated with fluoroquinolone resistance

Multidrug-resistant and extensively drug-resistant tuberculosis (M/XDR-TB) remains a global public-health challenge. Known mutations in quinolone resistance-determination regions cannot fully explain phenotypic fluoroquinolone (FQ) resistance in Mycobacterium tuberculosis (Mtb). The aim of this study was to look for novel mutations in Mtb associated with resistance to FQ drugs using whole-genome sequencing analysis. Whole-genome sequences of 659 Mtb strains, including 214 with phenotypic FQ resistance and 445 pan-susceptible isolates, were explored for mutations associated with FQ resistance overall and with resistance to individual FQ drugs (ofloxacin, levofloxacin, moxifloxacin and gatifloxacin). Three novel genes (recC, Rv2005c and PPE59) associated with FQ resistance were identified (P < 0.00001 based on screening analysis and absence of relevant mutations in a pan-susceptible validation set of 360 strains). Nine novel single nucleotide polymorphisms (SNPs), including in gyrB (G5383A and G6773A), gyrA (G7892A), recC (G725900C and G726857T/C), Rv2005c (C2251373G, G2251420C and C2251725T) and PPE59 (C3847269T), were used for diagnostic performance analysis. Enhancing the known SNP set with five of these novel SNPs, including gyrA [G7892A (Leu247Leu)], recC [G725900C (Leu893Leu) and G726857T/C (Arg484Arg)], Rv2005c [G2251420C (Pro205Arg)] and PPE59 [C3847269T (Asn35Asn)] increased the sensitivity of detection of FQ-resistant Mtb from 83.2% (178/214) to 86.9% (186/214) while maintaining 100% specificity (360/360). No specific mutation associated with resistance to only a single drug (ofloxacin, levofloxacin, moxifloxacin or gatifloxacin) was found. In conclusion, this study reports possible additional mutations associated with FQ resistance in Mtb.

Abattoir survey of bovine tuberculosis in tanta, centre of the Nile delta, with in silico analysis of gene mutations and protein-protein interactions of the involved mycobacteria

Bovine tuberculosis is a transboundary disease of high economic and public health burden worldwide. In this study, post-mortem examination of 750 cattle and buffalo in Tanta abattoir, Centre of the Nile Delta, revealed visible TB in 4% of animals and a true prevalence of 6.85% (95% CI: 5.3%-8.9%). Mycobacterial culture, histopathology and RT-PCR targeting all members of M. tuberculosis complex were performed, upon which 85%, 80% and 100% of each tested lesions were confirmed as TB, respectively. Mpb70-targeting PCR was conducted on ten RT-PCR positive samples for sequencing and identified nine Mycobacterium (M.) bovis strains and, interestingly, one M. tuberculosis (Mtb) strain from a buffalo. Bioinformatics tools were used for prediction of mutations, nucleotide polymorphisms, lineages, drug resistance and protein-protein interactions (PPI) of the sequenced strains. The Mtb strain was resistant to rifampicin, isoniazid and streptomycin, and to the best of our knowledge, this is the first report of multidrug resistant (MDR)-Mtb originating from buffaloes. Seven M. bovis strains were resistant to ethambutol and ethionamide. Such resistances were associated with KatG, rpoB, rpsL, embB and ethA genes mutations. Other mutations and nucleotide polymorphisms were also predicted, some are reported for the first time and require experimental work for validation. PPI revealed more interactions than what would be expected for a random set of proteins of similar size and had dense interactions between nodes that are biologically connected, as a group. Two M. bovis strains belonged to BOV AFRI lineage (Spoligotypes BOV 1; BOV 2) and eight strains belonged to East-Asian (Beijing) lineage. In conclusion, visible TB was prevalent in the study area, RT-PCR is the best to confirm the disease, MDR-Mtb is associated with buffalo TB, and mycobacteria of different lineages carry many resistance genes to chemotherapeutic agents used in treatment of human TB constituting a major public health risk.

Whole-genome sequencing analysis of multidrug-resistant Mycobacterium tuberculosis from Java, Indonesia

Introduction. Multidrug-resistant tuberculosis (MDR-TB) is a major public health problem globally, including in Indonesia. Whole-genome sequencing (WGS) analysis has rarely been used for the study of TB and MDR-TB in Indonesia.Aim. We evaluated the use of WGS for drug-susceptibility testing (DST) and to investigate the population structure of drug-resistant Mycobacterium tuberculosis in Java, Indonesia.Methodology. Thirty suspected MDR-TB isolates were subjected to MGIT 960 system (MGIT)-based DST and to WGS. Phylogenetic analysis was done using the WGS data. Results obtained using MGIT-based DST and WGS-based DST were compared.Results. Agreement between WGS and MGIT was 93.33 % for rifampicin, 83.33 % for isoniazid and 76.67 % for streptomycin but only 63.33 % for ethambutol. Moderate WGS-MGIT agreement was found for second-line drugs including amikacin, kanamycin and fluoroquinolone (73.33-76.67 %). MDR-TB was more common in isolates of the East Asian Lineage (63.3%). No evidence of clonal transmission of DR-TB was found among members of the tested population.Conclusion. Our study demonstrated the applicability of WGS for DST and molecular epidemiology of DR-TB in Java, Indonesia. We found no transmission of DR-TB in Indonesia.

Current topics of molecular mycobacteriology

The 2nd St. Petersburg Symposium on Tuberculosis and Mycobacteria: Molecular Approach, was held in St. Petersburg, Russia on 5-6 December 2018. A special issue of Infection, Genetics and Evolution will publish articles based on the selected presentations. In this paper, I will discuss some of the hot topics of molecular mycobacteriology highlighted at this meeting that I had the pleasure to organize and honor to chair. The symposium addressed interrelated fundamental and applied issues of modern mycobacteriology such as molecular evolution and phylogenomics, host-microbe interactions and pathogenesis, coevolution of M. tuberculosis with humans, new genomic and postgenomic technologies. Molecular methods for TB diagnostics and drug resistance detection are supported by WHO and whole genome/next generation sequencing presents a comprehensive approach. At the same time, cost and implementation of new methods for direct analysis of clinical samples and/or in low-resource settings remain a great challenge. A due attention was also given to the medically important nontuberculous mycobacteria. Assessment of spectrum of the circulating mycobacterial species in the Russian Federation and the countries of the European Union was presented and the underlying reasons of the observed diversity were discussed. To conclude, the symposium became a multidisciplinary event that was useful to promote networking and exchange of knowledge and experience. The next (third) symposium was planned to be organized in 2021.

Comparisons of whole-genome sequencing and phenotypic drug susceptibility testing for Mycobacterium tuberculosis causing MDR-TB and XDR-TB in Thailand

Drug-resistant tuberculosis (TB) is a major public health problem. There is little information regarding the genotypic-phenotypic association of anti-TB drugs, especially for second-line drugs. This study compared phenotypic drug susceptibility testing (DST) with predictions based on whole-genome sequencing (WGS) data for 266 Mycobacterium tuberculosis isolates. Phenotypic DST used the standard proportional method. Clinical isolates of M. tuberculosis collected in Thailand between 1998 and 2013 comprised 51 drug-sensitive strains, six mono-resistant strains, two multiple-resistant strains, 88 multi-drug-resistant strains, 95 pre-extensively drug-resistant strains and 24 extensively drug-resistant strains. WGS analysis was performed using the computer programs PhyResSE and TB-Profiler. TB-Profiler had higher average concordance with phenotypic DST than PhyResSE for both first-line (91.96% vs. 91.4%) and second-line (79.67% vs. 78.20%) anti-TB drugs. The average sensitivity for all anti-TB drugs was also higher (83.13% vs. 72.08%) with slightly lower specificity (83.50% vs. 86.68%). Regardless of the program used, isoniazid, rifampicin and amikacin had the highest concordance with phenotypic DST (96.2%, 93.5% and 95.6%, respectively). Ethambutol, ethionamide and fluoroquinolones had the lowest concordance (87.34%, 81.44% and 73.85%, respectively). Concordance rates of ofloxacin (a second-generation fluoroquinolone), levofloxacin, moxifloxacin and gatifloxacin (third- and fourth-generation fluoroquinolones) were 91.79%, 76.62%, 72.64% and 57.35%, respectively. Discordance between phenotypic and WGS-based DSTs may be due, in part, to the choice of critical concentration and variable reproducibility of the phenotypic tests. It may also be due to limitations of the mutation databases (especially for the second-line drugs) and the analysis program used. Mutations related to fluoroquinolone resistance, especially the later generations, need to be identified.

Drug resistance mechanisms and drug susceptibility testing for tuberculosis

Tuberculosis (TB) caused by Mycobacterium tuberculosis (MTB) is the deadliest infectious disease and the associated global threat has worsened with the emergence of drug resistance, in particular multidrug-resistant TB (MDR-TB) and extensively drug-resistant TB (XDR-TB). Although the World Health Organization (WHO) End-TB Strategy advocates for universal access to antimicrobial susceptibility testing, this is not widely available and/or it is still underused. The majority of drug resistance in clinical MTB strains is attributed to chromosomal mutations. Resistance-related mutations could also exert certain fitness cost to the drug-resistant MTB strains and growth fitness could be restored by the presence of compensatory mutations. Understanding these underlying mechanisms could provide an important insight into TB pathogenesis and predict the future trend of MDR-TB global pandemic. This review covers the mechanisms of resistance in MTB and provides a comprehensive overview of current phenotypic and molecular approaches for drug susceptibility testing, with particular attention to the methods endorsed and recommended by the WHO.

Use of whole genome sequencing in surveillance of drug resistant tuberculosis

INTRODUCTION

The threat of resistance to anti-tuberculosis drugs is of global concern. Current efforts to monitor resistance rely on phenotypic testing where cultured bacteria are exposed to critical concentrations of the drugs. Capacity for such testing is low in TB endemic countries. Drug resistance is caused by mutations in the Mycobacterium tuberculosis genome and whole genome sequencing to detect these mutations offers an alternative means of assessing resistance. Areas covered: The challenges of assessing TB drug resistance are discussed. Progress in elucidating the M. tuberculosis resistome and evidence of the accuracy of next generation sequencing for detecting resistance is reviewed. Expert Commentary: There are considerable advantages to using next generation sequencing for TB drug resistance surveillance. Accuracy is high for detecting resistance to the major first-line drugs but is currently lower for the second-line drugs due to our incomplete knowledge regarding resistance causing mutations. With the advances in sequencing technology and the opportunity to replace phenotypic drug susceptibility testing with safer and more cost effective methods it would appear that the question is when to implement. Current bottlenecks are sample extraction to allow whole genome sequencing directly from sputum and the lack of bioinformatics expertise in some TB endemic countries.

Dissecting whole-genome sequencing-based online tools for predicting resistance in Mycobacterium tuberculosis: can we use them for clinical decision guidance?

Whole-genome sequencing (WGS)-based bioinformatics platforms for the rapid prediction of resistance will soon be implemented in the Tuberculosis (TB) laboratory, but their accuracy assessment still needs to be strengthened. Here, we fully-sequenced a total of 54 multidrug-resistant (MDR) and five susceptible TB strains and performed, for the first time, a simultaneous evaluation of the major four free online platforms (TB Profiler, PhyResSE, Mykrobe Predictor and TGS-TB). Overall, the sensitivity of resistance prediction ranged from 84.3% using Mykrobe predictor to 95.2% using TB profiler, while specificity was higher and homogeneous among platforms. TB profiler revealed the best performance robustness (sensitivity, specificity, PPV and NPV above 95%), followed by TGS-TB (all parameters above 90%). We also observed a few discrepancies between phenotype and genotype, where, in some cases, it was possible to pin-point some "candidate" mutations (e.g., in the rpsL promoter region) highlighting the need for their confirmation through mutagenesis assays and potential review of the anti-TB genetic databases. The rampant development of the bioinformatics algorithms and the tremendously reduced time-frame until the clinician may decide for a definitive and most effective treatment will certainly trigger the technological transition where WGS-based bioinformatics platforms could replace phenotypic drug susceptibility testing for TB.

Molecular epidemiology of tuberculosis

The application of genotyping tools allowed us to discriminate between the Mycobacterium tuberculosis isolates obtained in the laboratory. The differentiation between single strains opened the door to molecular epidemiology studies, which had helped us to progress in our knowledge of how this pathogen is transmitted in the progressively more complex socio-epidemiological scenario. The genetic stability of this microorganism led to develop specific methodologies, which are thoroughly revised in this chapter. In addition to their application in epidemiology, we review, how they can offer a response to different diagnostic and clinical challenges. Finally, we focus on describing the novel genomic revolution we are experiencing in the analysis of tuberculosis, the methodology in which it is based and the novel possibilities it offers, including new routes of integrating both the molecular and genomic languages in innovative post-genomic proposals, better suited to our real-life context.

Prevalence and extent of heteroresistance by next generation sequencing of multidrug-resistant tuberculosis

Amplicon-based Next Generation Sequencing (NGS) is an emerging method for Mycobacterium tuberculosis drug susceptibility testing (DST) but has not been well described. We examined 158 clinical multidrug-resistant M. tuberculosis isolates via NGS of 11 resistance-associated gene regions covering 3519 nucleotides. Across these gene regions, complete resistance or heteroresistance (defined as 1%-99% mutation) was present in at least one isolate in 6.3% of loci. The number of isolates with heteroresistance was highest for gyrA codon 94, rpoB codons 526 and 531, and embB codons 306, 372 and 406 (range 11-26% of isolates exhibited heteroresistance). 57% of MDR strains had heteroresistance of one or more recognized resistance-associated mutation. Heteroresistant loci generally exhibited high or low degrees of mutation (>90% or <10%). The deep sensitivity of NGS for detecting low level pncA heteroresistance appeared to improve genotypic-phenotypic PZA susceptibility correlations over that of Sanger. NGS demonstrates that heteroresistance in TB in the regions of key genes is common and will need to be bioinformatically managed. The clinical significance of such heteroresistance is unclear, and further study of pncA should be pursued.

Removing the bottleneck in whole genome sequencing of Mycobacterium tuberculosis for rapid drug resistance analysis: a call to action

Whole genome sequencing (WGS) can provide a comprehensive analysis of Mycobacterium tuberculosis mutations that cause resistance to anti-tuberculosis drugs. With the deployment of bench-top sequencers and rapid analytical software, WGS is poised to become a useful tool to guide treatment. However, direct sequencing from clinical specimens to provide a full drug resistance profile remains a serious challenge. This article reviews current practices for extracting M. tuberculosis DNA and possible solutions for sampling sputum. Techniques under consideration include enzymatic digestion, physical disruption, chemical degradation, detergent solubilization, solvent extraction, ligand-coated magnetic beads, silica columns, and oligonucleotide pull-down baits. Selective amplification of genomic bacterial DNA in sputum prior to WGS may provide a solution, and differential lysis to reduce the levels of contaminating human DNA is also being explored. To remove this bottleneck and accelerate access to WGS for patients with suspected drug-resistant tuberculosis, it is suggested that a coordinated and collaborative approach be taken to more rapidly optimize, compare, and validate methodologies for sequencing from patient samples.