A complete high-quality MinION nanopore assembly of an extensively drug-resistant Mycobacterium tuberculosis Beijing lineage strain identifies novel variation in repetitive PE/PPE gene regions

Interred mechanisms of resistance and host immune evasion revealed through network-connectivity analysis of M. tuberculosis complex graph pangenome

Mycobacterium tuberculosis complex successfully adapts to environmental pressures through mechanisms of rapid adaptation which remain poorly understood despite knowledge gained through decades of research. In this study, we used 110 reference-quality, complete de novo assembled, long-read sequenced clinical genomes to study patterns of structural adaptation through a graph-based pangenome analysis, elucidating rarely studied mechanisms that enable enhanced clinical phenotypes offering a novel perspective to the species' adaptation. Across isolates, we identified a pangenome of 4,325 genes (3,767 core and 558 accessory), revealing 290 novel genes, and a substantially more complete account of difficult-to-sequence esx/pe/pgrs/ppe genes. Seventy-four percent of core genes were deemed non-essential in vitro, 38% of which support the pathogen's survival in vivo, suggesting a need to broaden current perspectives on essentiality. Through information-theoretic analysis, we reveal the ppe genes that contribute most to the species' diversity-several with known consequences for antigenic variation and immune evasion. Construction of a graph pangenome revealed topological variations that implicate genes known to modulate host immunity (Rv0071-73, Rv2817c, cas2), defense against phages/viruses (cas2, csm6, and Rv2817c-2821c), and others associated with host tissue colonization. Here, the prominent trehalose transport pathway stands out for its involvement in caseous granuloma catabolism and the development of post-primary disease. We show paralogous duplications of genes implicated in bedaquiline (mmpL5 in all L1 isolates) and ethambutol (embC-A) resistance, with a paralogous duplication of its regulator (embR) in 96 isolates. We provide hypotheses for novel mechanisms of immune evasion and antibiotic resistance through gene dosing that can escape detection by molecular diagnostics.IMPORTANCEM. tuberculosis complex (MTBC) has killed over a billion people in the past 200 years alone and continues to kill nearly 1.5 million annually. The pathogen has a versatile ability to diversify under immune and drug pressure and survive, even becoming antibiotic persistent or resistant in the face of harsh chemotherapy. For proper diagnosis and design of an appropriate treatment regimen, a full understanding of this diversification and its clinical consequences is desperately needed. A mechanism of diversification that is rarely studied systematically is MTBC's ability to structurally change its genome. In this article, we have de novo assembled 110 clinical genomes (the largest de novo assembled set to date) and performed a pangenomic analysis. Our pangenome provides structural variation-based hypotheses for novel mechanisms of immune evasion and antibiotic resistance through gene dosing that can compromise molecular diagnostics and lead to further emergence of antibiotic resistance.

High-Accuracy Long-Read Sequencing of Mycobacterium tuberculosis PSNK363 Isolated From the Democratic People's Republic of Korea

Long-read sequencing is a valuable technique for high-precision genome analysis. Despite the widespread use of the Mycobacterium tuberculosis H37Rv genome sequence as a reference for genetic variation analysis, its suitability for comparing clinical strains is limited. Therefore, we constructed the first known whole genome of a clinical M. tuberculosis strain, PSNK363, isolated from the Democratic People's Republic of Korea, using high-quality high-fidelity (HiFi) read sequencing and compared its genetic variations to those of H37Rv. PSNK363 was cultured to obtain genomic DNA, which was subjected to de novo whole-genome assembly using PacBio Sequel II with long-read HiFi sequencing. The sequences were compared to the reference genome H37Rv. HiFi long-read sequencing of M. tuberculosis PSNK363, with an accuracy of 99.99%, revealed a single circular chromosome of 4,422,110 bp, which is 10,578 bp longer than the H37Rv chromosome. The assembly had an average G + C content of 65.6%, 4079 protein-coding sequences, 53 tRNA genes, and 3 rRNA genes. Most genes (72.7%) were assigned as putative functions, whereas the remaining 27.3% were annotated as hypothetical. Comparison with H37Rv revealed a large inversion in the PSNK363 genome, which contains most of the deletion and insertion variants. M. tuberculosis PSNK363 had a longer genome sequence, more protein-coding genes, and a larger inversion region than H37Rv. High-accuracy whole-genome sequencing of PSNK363 holds the potential for enriching virulence databases and identifying informative loci for drug resistance analysis in M. tuberculosis isolates in the Democratic People's Republic of Korea.

Comparisons of genome assembly tools for characterization of Mycobacterium tuberculosis genomes using hybrid sequencing technologies

Background

Next-generation sequencing of Mycobacterium tuberculosis, the infectious agent causing tuberculosis, is improving the understanding of genomic diversity of circulating lineages and strain-types, and informing knowledge of drug resistance mutations. An increasingly popular approach to characterizing M. tuberculosis genomes (size: 4.4 Mbp) and variants (e.g., single nucleotide polymorphisms (SNPs)) involves the de novo assembly of sequence data.

Methods

We compared the performance of genome assembly tools (Unicycler, RagOut, and RagTag) on sequence data from nine drug resistant M. tuberculosis isolates (multi-drug (MDR) n = 1; pre-extensively-drug (pre-XDR) n = 8) generated using Illumina HiSeq, Oxford Nanopore Technology (ONT) PromethION, and PacBio platforms.

Results

Our investigation found that Unicycler-based assemblies had significantly higher genome completeness (~98.7%; p values = 0.01) compared to other assembler tools (RagOut = 98.6%, and RagTag = 98.6%). The genome assembly sizes (bp) across isolates and sequencers based on RagOut was significantly longer (p values < 0.001) (4,418,574 ± 8,824 bp) than Unicycler and RagTag assemblies (Unicycler = 4,377,642 ± 55,257 bp, and RagTag = 4,380,711 ± 51,164 bp). RagOut-based assemblies had the fewest contigs (~32) and the longest genome size (4,418,574 bp; vs. H37Rv reference size 4,411,532 bp) and therefore were chosen for downstream analysis. Pan-genome analysis of Illumina and PacBio hybrid assemblies revealed the greatest number of detected genes (4,639 genes; H37Rv reference contains 3,976 genes), while Illumina and ONT hybrid assemblies produced the highest number of SNPs. The number of genes from hybrid assemblies with ONT and PacBio long-reads (mean: 4,620 genes) was greater than short-read assembly alone (4,478 genes). All nine RagOut hybrid genome assemblies detected known mutations in genes associated with MDR-TB and pre-XDR-TB.

Conclusions

Unicycler software performed the best in terms of achieving contiguous genomes, whereas RagOut improved the quality of Unicycler's genome assemblies by providing a longer genome size. Overall, our approach has demonstrated that short-read and long-read hybrid assembly can provide a more complete genome assembly than short-read assembly alone by detecting pan-genomes and more genes, including IS6110, and SNPs.

Description of Bacterial RNA Transcripts Detected in Mycobacterium tuberculosis - Infected Cells from Peripheral Human Granulomas using Single Cell RNA Sequencing

Mycobacterium tuberculosis (Mtb) remains a global human health threat and a significant cause of human morbidity and mortality. We document here the capture of Mtb transcripts in libraries designed to amplify eukaryotic mRNA. These reads are often considered spurious or nuisance and are rarely investigated. Because of early literature suggesting the possible presence of polyadenylated transcripts in Mtb RNA, we included the H37Rv Mtb reference genome when assembling scRNA seq libraries from fine needle aspirate samples from patients presenting at the TB clinic, Port Moresby General Hospital, Papua New Guinea. We used 10X Genomics single-cell RNA sequencing transcriptomics pipeline, which initiates mRNA amplification with poly-T primers on ~30-micron beads designed to capture, in this case, human mRNA associated with individual cells in the clinical samples. Utilizing the 10X Genomics Cell Ranger tool to align sequencing reads, we consistently detected bacterial small and large ribosomal subunit RNA sequences (rrs and rrl, respectively) and other bacterial gene transcripts in the cell culture and patient samples. We interpret Mtb reads associated with the host cell's unique molecular identifier (UMI) and transcriptome to indicate infection of that individual host cell. The Mtb transcripts detected showed frequent sequence variation from the reference genome, with greater than 90% of the rrs or rrl reads from many clinical samples having at least 1 sequence difference compared to the H37Rv reference genome. The data presented includes only bacterial sequences from patients with TB infections that were confirmed by the hospital pathology lab using acid-fast microscopy and/or GeneXpert analysis. The repeated, non-random nature of the sequence variations detected in Mtb rrs and rrl transcripts from multiple patients, suggests that, even though this appears to be a stochastic process, there is possibly some selective pressure that limits the types and locations of sequence variation allowed. The variation does not appear to be entirely artefactual, and it is hypothesized that it could represent an additional mechanism of adaptation to enhance bacterial fitness against host defenses or chemotherapy.

A practical comparison of the next-generation sequencing platform and assemblers using yeast genome

Assembling fragmented whole-genomic information from the sequencing data is an inevitable process for further genome-wide research. However, it is intricate to select the appropriate assembly pipeline for unknown species because of the species-specific genomic properties. Therefore, our study focused on relatively more static proclivities of sequencing platforms and assembly algorithms than the fickle genome sequences. A total of 212 draft and polished de novo assemblies were constructed under the different sequencing platforms and assembly algorithms with the repetitive yeast genome. Our comprehensive data indicated that sequencing reads from Oxford Nanopore with R7.3 flow cells generated more continuous assemblies than those derived from the PacBio Sequel, although the homopolymer-based assembly errors and chimeric contigs exist. In addition, the comparison between two second-generation sequencing platforms showed that Illumina NovaSeq 6000 provides more accurate and continuous assembly in the second-generation-sequencing-first pipeline, but MGI DNBSEQ-T7 provides a cheap and accurate read in the polishing process. Furthermore, our insight into the relationship among the computational time, read length, and coverage depth provided clues to the optimal pipelines of yeast assembly.

The PE-PPE Family of Mycobacterium tuberculosis: Proteins in Disguise

Mycobacterium tuberculosis has thrived in parallel with humans for millennia, and despite our efforts, M. tuberculosis continues to plague us, currently infecting a third of the world's population. The success of M. tuberculosis has recently been attributed, in part, to the PE-PPE family; a unique collection of 168 proteins fundamentally involved in the pathogenesis of M. tuberculosis. The PE-PPE family proteins have been at the forefront of intense research efforts since their discovery in 1998 and whilst our knowledge and understanding has significantly advanced over the last two decades, many important questions remain to be elucidated. This review consolidates and examines the vast body of existing literature regarding the PE-PPE family proteins, with respect to the latest developments in elucidating their evolution, structure, subcellular localisation, function, and immunogenicity. This review also highlights significant inconsistencies and contradictions within the field. Additionally, possible explanations for these knowledge gaps are explored. Lastly, this review poses many important questions, which need to be addressed to complete our understanding of the PE-PPE family, as well as highlighting the challenges associated with studying this enigmatic family of proteins. Further research into the PE-PPE family, together with technological advancements in genomics and proteomics, will undoubtedly improve our understanding of the pathogenesis of M. tuberculosis, as well as identify key targets/candidates for the development of novel drugs, diagnostics, and vaccines.

A comprehensive update to the Mycobacterium tuberculosis H37Rv reference genome

H37Rv is the most widely used Mycobacterium tuberculosis strain, and its genome is globally used as the M. tuberculosis reference sequence. Here, we present Bact-Builder, a pipeline that uses consensus building to generate complete and accurate bacterial genome sequences and apply it to three independently cultured and sequenced H37Rv aliquots of a single laboratory stock. Two of the 4,417,942 base-pair long H37Rv assemblies are 100% identical, with the third differing by a single nucleotide. Compared to the existing H37Rv reference, the new sequence contains ~6.4 kb additional base pairs, encoding ten new regions that include insertions in PE/PPE genes and new paralogs of esxN and esxJ, which are differentially expressed compared to the reference genes. New sequencing and de novo assemblies with Bact-Builder confirm that all 10 regions, plus small additional polymorphisms, are also present in the commonly used H37Rv strains NR123, TMC102, and H37Rv1998. Thus, Bact-Builder shows promise as an improved method to perform accurate and reproducible de novo assemblies of bacterial genomes, and our work provides important updates to the primary M. tuberculosis reference genome.

Nanopore Electrochemistry for Pathogen Detection

Pathogen infections have seriously threatened human health, and there is an urgent demand for rapid and efficient pathogen identification to provide instructions in clinical diagnosis and therapeutic intervention. Recently, nanopore technology, a rapidly maturing technology which delivers ultrasensitive sensing and high throughput in real-time and at low cost, has achieved success in pathogen detection. Furthermore, the remarkable development of nanopore sequencing, for example, the MinION sequencer from Oxford Nanopore Technologies (ONT) as a competitive sequencing technology, has facilitated the rapid analysis of disease-related microbiomes at the whole-genome level and on a large scale. Here, we highlighted recent advances in nanopore approaches for pathogen detection at the single-molecule level. We also overviewed the applications of nanopore sequencing in pathogenic bacteria identification and diagnosis. In the end, we discussed the challenges and future developments of nanopore technology as promising tools for the management of infections, which may be helpful to aid understanding as well as decision-making.

Accuracy of an amplicon-sequencing nanopore approach to identify variants in tuberculosis drug-resistance-associated genes

A rapid and accurate diagnostic assay represents an important means to detect Mycobacterium tuberculosis, identify drug-resistant strains and ensure treatment success. Currently employed techniques to diagnose drug-resistant tuberculosis include slow phenotypic tests or more rapid molecular assays that evaluate a limited range of drugs. Whole-genome-sequencing-based approaches can detect known drug-resistance-conferring mutations and novel variations; however, the dependence on growing samples in culture, and the associated delays in achieving results, represents a significant limitation. As an alternative, targeted sequencing strategies can be directly performed on clinical samples at high throughput. This study proposes a targeted sequencing assay to rapidly detect drug-resistant strains of M. tuberculosis using the Nanopore MinION sequencing platform. We designed a single-tube assay that targets nine genes associated with drug resistance to seven drugs and two phylogenetic-determining regions to determine strain lineage and tested it in nine clinical isolates and six sputa. The study's main aim is to calibrate MinNION variant calling to detect drug-resistance-associated mutations with different frequencies to match the accuracy of Illumina (the current gold-standard sequencing technology) from both culture and sputum samples. After calibrating Nanopore MinION variant calling, we demonstrated 100% agreement between Illumina WGS and our MinION set up to detect known drug resistance and phylogenetic variants in our dataset. Importantly, other variants in the amplicons are also detected, decreasing the recall. We identify minority variants and insertions/deletions as crucial bioinformatics challenges to fully reproduce Illumina WGS results.

Nucleotide amplification and sequencing of the GC-rich region between matrix and fusion protein genes of peste des petits ruminants virus

Peste des petits ruminants virus (PPRV) causes a highly devastating disease of sheep and goats, that threatens the conservation of small wild ruminants. The development of PPRV vaccines, diagnostics and therapeutics, greatly depends on in-depth genomic data. Yet, high guanine-cytosine (GC) content between matrix (M) and fusion (F) genes of PPRV poses difficulty for both primer design and nucleotide amplification. In turn, this has led into absence or low nucleotide sequence coverage in this region. This poses a risk of missing important part of the genome that could help to infer viral evolution. Here, an overlapping long-read primer-based amplification strategy was developed to amplify the GC-rich fragments between M-F gene junction using nexus gradient polymerase chain reaction (PCR). The resulting amplicons were sequenced by dideoxynucleotide cycle sequencing and compared with other PPRV nucleotide sequences available at GenBank. Our findings indicate clear PCR amplification products with expected size of the GC-rich fragments on agarose gel electrophoresis. The sequencing results of these fragments indicate 99.5 % nucleotide identity with PPRV strain KY628761. An extremely difficult PCR target of 67.4 % GC contents was successfully amplified and sequenced using this long-read primer approach. The long-read primer set may be used in tiling multiplex PCR for complete genome sequencing of PPRV.

Evaluation of whole-genome sequence data analysis approaches for short- and long-read sequencing of Mycobacterium tuberculosis

Whole-genome sequencing (WGS) of Mycobacterium tuberculosis (MTB) isolates can be used to get an accurate diagnosis, to guide clinical decision making, to control tuberculosis (TB) and for outbreak investigations. We evaluated the performance of long-read (LR) and/or short-read (SR) sequencing for anti-TB drug-resistance prediction using the TBProfiler and Mykrobe tools, the fraction of genome recovery, assembly accuracies and the robustness of two typing approaches based on core-genome SNP (cgSNP) typing and core-genome multi-locus sequence typing (cgMLST). Most of the discrepancies between phenotypic drug-susceptibility testing (DST) and drug-resistance prediction were observed for the first-line drugs rifampicin, isoniazid, pyrazinamide and ethambutol, mainly with LR sequence data. Resistance prediction to second-line drugs made by both TBProfiler and Mykrobe tools with SR- and LR-sequence data were in complete agreement with phenotypic DST except for one isolate. The SR assemblies were more accurate than the LR assemblies, having significantly (P <0.05) fewer indels and mismatches per 100 kbp. However, the hybrid and LR assemblies had slightly higher genome fractions. For LR assemblies, Canu followed by Racon, and Medaka polishing was the most accurate approach. The cgSNP approach, based on either reads or assemblies, was more robust than the cgMLST approach, especially for LR sequence data. In conclusion, anti-TB drug-resistance prediction, particularly with only LR sequence data, remains challenging, especially for first-line drugs. In addition, SR assemblies appear more accurate than LR ones, and reproducible phylogeny can be achieved using cgSNP approaches.

Whole Genome Sequencing in the Management of Non-Tuberculous Mycobacterial Infections

Infections caused by non-tuberculous mycobacteria (NTM) have been a public health problem in recent decades and contribute significantly to the clinical and economic burden globally. The diagnosis of infections is difficult and time-consuming and, in addition, the conventional diagnostics tests do not have sufficient discrimination power in species identification due to cross-reactions and not fully specific probes. However, technological advances have been made and the whole genome sequencing (WGS) method has been shown to be an essential part of routine diagnostics in clinical mycobacteriology laboratories. The use of this technology has contributed to the characterization of new species of mycobacteria, as well as the identification of gene mutations encoding resistance and virulence factors. Sequencing data also allowed to track global outbreaks of nosocomial NTM infections caused by M. abscessus complex and M. chimaera. To highlight the utility of WGS, we summarize recent scientific studies on WGS as a tool suitable for the management of NTM-induced infections in clinical practice.

The within-host evolution of antimicrobial resistance in Mycobacterium tuberculosis

Tuberculosis (TB) has been responsible for the greatest number of human deaths due to an infectious disease in general, and due to antimicrobial resistance (AMR) in particular. The etiological agents of human TB are a closely-related group of human-adapted bacteria that belong to the Mycobacterium tuberculosis complex (MTBC). Understanding how MTBC populations evolve within-host may allow for improved TB treatment and control strategies. In this review, we highlight recent works that have shed light on how AMR evolves in MTBC populations within individual patients. We discuss the role of heteroresistance in AMR evolution, and review the bacterial, patient and environmental factors that likely modulate the magnitude of heteroresistance within-host. We further highlight recent works on the dynamics of MTBC genetic diversity within-host, and discuss how spatial substructures in patients' lungs, spatiotemporal heterogeneity in antimicrobial concentrations and phenotypic drug tolerance likely modulates the dynamics of MTBC genetic diversity in patients during treatment. We note the general characteristics that are shared between how the MTBC and other bacterial pathogens evolve in humans, and highlight the characteristics unique to the MTBC.

Genomic epidemiology of tuberculosis in eastern Malaysia: insights for strengthening public health responses

Tuberculosis is a leading public health priority in eastern Malaysia. Knowledge of the genomic epidemiology of tuberculosis can help tailor public health interventions. Our aims were to determine tuberculosis genomic epidemiology and characterize resistance mutations in the ethnically diverse city of Kota Kinabalu, Sabah, located at the nexus of Malaysia, Indonesia, Philippines and Brunei. We used an archive of prospectively collected Mycobacterium tuberculosis samples paired with epidemiological data. We collected sputum and demographic data from consecutive consenting outpatients with pulmonary tuberculosis at the largest tuberculosis clinic from 2012 to 2014, and selected samples from tuberculosis inpatients from the tertiary referral centre during 2012-2014 and 2016-2017. Two hundred and eight M. tuberculosis sequences were available for analysis, representing 8 % of cases notified during the study periods. Whole-genome phylogenetic analysis demonstrated that most strains were lineage 1 (195/208, 93.8 %), with the remainder being lineages 2 (8/208, 3.8 %) or 4 (5/208, 2.4 %). Lineages or sub-lineages were not associated with patient ethnicity. The lineage 1 strains were diverse, with sub-lineage 1.2.1 being dominant (192, 98 %). Lineage 1.2.1.3 isolates were geographically most widely distributed. The greatest diversity occurred in a border town sub-district. The time to the most recent common ancestor for the three major lineage 1.2.1 clades was estimated to be the year 1966 (95 % HPD 1948-1976). An association was found between failure of culture conversion by week 8 of treatment and infection with lineage 2 (4/6, 67 %) compared with lineage 1 strains (4/83, 5 %) (P<0.001), supporting evidence of greater virulence of lineage 2 strains. Eleven potential transmission clusters (SNP difference ≤12) were identified; at least five included people living in different sub-districts. Some linked cases spanned the whole 4-year study period. One cluster involved a multidrug-resistant tuberculosis strain matching a drug-susceptible strain from 3 years earlier. Drug resistance mutations were uncommon, but revealed one phenotype-genotype mismatch in a genotypically multidrug-resistant isolate, and rare nonsense mutations within the katG gene in two isolates. Consistent with the regionally mobile population, M. tuberculosis strains in Kota Kinabalu were diverse, although several lineage 1 strains dominated and were locally well established. Transmission clusters - uncommonly identified, likely attributable to incomplete sampling - showed clustering occurring across the community, not confined to households or sub-districts. The findings indicate that public health priorities should include active case finding and early institution of tuberculosis management in mobile populations, while there is a need to upscale effective contact investigation beyond households to include other contacts within social networks.

Systematic Review of Mutations Associated with Isoniazid Resistance Points to Continuing Evolution and Subsequent Evasion of Molecular Detection, and Potential for Emergence of Multidrug Resistance in Clinical Strains of Mycobacterium tuberculosis

Molecular testing is rapidly becoming an integral component of global tuberculosis (TB) control. Uncommon mechanisms of resistance escape detection by these platforms and undermine our ability to contain outbreaks. This article is a systematic review of published articles that reported isoniazid (INH) resistance-conferring mutations between September 2013 and December 2019. The genes katG, inhA, and fabG1, and the intergenic region oxyR'-ahpC were considered in this review. Fifty-two articles were included that described 9,306 clinical isolates (5,804 INH resistant [INH^r] and 3,502 INH susceptible [INH^s]) from 31 countries. The three most frequently mutated loci continue to be locus 315 of katG (katG315; n = 4,271), locus -15 of inhA (inhA-15; n = 787), and locus -8 of inhA (inhA-8; 106). However, the diagnostic value of inhA-8 is far lower than previously thought, as it only appears in 25 (0.4%) of the INH^r isolates lacking the first two mutations. I catalogued 45 new loci (29 katG, nine inhA, and seven ahpC) associated with INH resistance and identified 59 loci (common to this and previous reviews) as a reliable basis for molecular diagnostics. Including all observed mutations provides a cumulative sensitivity of 85.6%. In 14.4% of resistant isolates, no mechanism of resistance was detected, making them likely to escape molecular detection, and in the case of INH monoresistance, likely to convert to multidrug-resistant TB (MDR-TB). Integrating the information cataloged in this study into current diagnostic tools is essential for combating the emergence of MDR-TB, and its exclusion can lead to an unintended selection against common mechanisms and to diversifying evolution. Observation of many low-frequency resistance-conferring mutations points to an advantage of whole-genome sequencing (WGS) for diagnostics. Finally, I provide five recommendations for future diagnostic platforms.

The Applications of Nanopore Sequencing Technology in Pathogenic Microorganism Detection

Infectious diseases are major threats to human health and lead to a serious public health burden. The emergence of new pathogens and the mutation of known pathogens challenge our ability to diagnose and control infectious diseases. Nanopore sequencing technology exhibited versatile applications in pathogenic microorganism detection due to its flexible data throughput. This review article introduced the applications of nanopore sequencing in clinical microbiology and infectious diseases management, including the monitoring of emerging infectious diseases outbreak, identification of pathogen drug resistance, and disease-related microbial communities characterization.

Characterization of Fecal Microbiota with Clinical Specimen Using Long-Read and Short-Read Sequencing Platform

Accurate and rapid identification of microbiotic communities using 16S ribosomal (r)RNA sequencing is a critical task for expanding medical and clinical applications. Next-generation sequencing (NGS) is widely considered a practical approach for direct application to communities without the need for in vitro culturing. In this report, a comparative evaluation of short-read (Illumina) and long-read (Oxford Nanopore Technologies (ONT)) platforms toward 16S rRNA sequencing with the same batch of total genomic DNA extracted from fecal samples is presented. Different 16S gene regions were amplified, bar-coded, and sequenced using the Illumina MiSeq and ONT MinION sequencers and corresponding kits. Mapping of the sequenced amplicon using MinION to the entire 16S rRNA gene was analyzed with the cloud-based EPI2ME algorithm. V3–V4 reads generated using MiSeq were aligned by applying the CLC genomics workbench. More than 90% of sequenced reads generated using distinct sequencers were accurately classified at the genus or species level. The misclassification of sequenced reads at the species level between the two approaches was less substantial as expected. Taken together, the comparative results demonstrate that MinION sequencing platform coupled with the corresponding algorithm could function as a practicable strategy in classifying bacterial community to the species level.

Application of Targeted Next-Generation Sequencing Assay on a Portable Sequencing Platform for Culture-Free Detection of Drug-Resistant Tuberculosis from Clinical Samples

Targeted next-generation sequencing (tNGS) has emerged as a comprehensive alternative to existing methods for drug susceptibility testing (DST) of Mycobacterium tuberculosis from patient sputum samples for clinical diagnosis of drug-resistant tuberculosis (DR-TB). However, the complexity of sequencing platforms has limited their uptake in low-resource settings. The goal of this study was to evaluate the use of the tNGS-based DST solution Genoscreen Deeplex Myc-TB, for use on the compact, low-cost Oxford Nanopore Technologies MinION sequencer. One hundred four DNA samples extracted from smear-positive sputum sediments, previously sequenced using the Deeplex assay on an Illumina MiniSeq, were resequenced on MinION after applying a custom library preparation. MinION read quality, mapping statistics, and variant calling were computed using an in-house pipeline and compared to the reference MiniSeq data. The average percentage of MinION reads mapped to an H37RV reference genome was 90.8%, versus 99.5% on MiniSeq. The mean depths of coverage were 4,151× and 4,177× on MinION and MiniSeq, respectively, with heterogeneous distribution across targeted genes. Composite reference coverage breadth was >99% for both platforms. We observed full concordance between technologies in reporting the clinically relevant drug-resistant markers, including full gene deletions. In conclusion, we demonstrated that the workflow and sequencing data obtained from Deeplex on MinION are comparable to those for the MiniSeq, despite the higher raw error rates on MinION, with the added advantage of MinION's portability, versatility, and low capital costs. Targeted NGS on MinION is a promising DST solution for rapidly providing clinically relevant data to manage complex DR-TB cases.

Genetic Identification and Drug-Resistance Characterization of Mycobacterium tuberculosis Using a Portable Sequencing Device. A Pilot Study

Clinical management of tuberculosis (TB) in endemic areas is often challenged by a lack of resources including laboratories for Mycobacterium tuberculosis (Mtb) culture. Traditional phenotypic drug susceptibility testing for Mtb is costly and time consuming, while PCR-based methods are limited to selected target loci. We herein utilized a portable, USB-powered, long-read sequencing instrument (MinION), to investigate Mtb genomic DNA from clinical isolates to determine the presence of anti-TB drug-resistance conferring mutations. Data analysis platform EPI2ME and antibiotic-resistance analysis using the real time ARMA workflow, identified Mtb species as well as extensive resistance gene profiles. The approach was highly sensitive, being able to detect almost all described drug resistance conferring mutations based on previous whole genome sequencing analysis. Our findings are supportive of the practical use of this system as a suitable method for the detection of antimicrobial resistance genes, and effective in providing Mtb genomic information. Future improvements in the error rate through statistical analysis, drug resistance prediction algorithms and reference databases would make this a platform suited for the clinical setting. The small size, relatively inexpensive cost of the device, as well as its rapid and simple library preparation protocol and analysis, make it an attractive option for settings with limited laboratory infrastructure.

Comparison of two molecular diagnostic methods for identifying Beijing genotype of Mycobacterium tuberculosis

BACKGROUND AND OBJECTIVES

The Beijing family of Mycobacterium tuberculosis has been identified as a severe pathogen among this species and found in many clinical isolates during the last decade. Early identification of such genotype is important for better prevention and treatment of tuberculosis. The present study performed to compare the efficiency of Real-Time PCR and IS6110-Based Inverse PCR methods to identify the Beijing family.

MATERIALS AND METHODS

This study was carried out on 173 clinical isolates of Mycobacterium tuberculosis complex in Golestan Province, northern Iran. DNA extraction performed by boiling and determining the Beijing and non-Beijing strains carried out using Real-Time PCR and IS6110-Based Inverse PCR.

RESULTS

In both Real-Time PCR and IS6110-Based Inverse PCR method, 24 specimens (13.9%) of the Beijing family were identified and the result of the IS6110-Based Inverse PCR method showed that all the Beijing strains in this region belonged to the Ancient Beijing sub-lineage.

CONCLUSION

Although the efficacy of the two methods in the diagnosis of the Beijing family is similar, the IS6110-Based Inverse PCR is more applicable to the ability to detect new and old Beijing family.

Comparison of routine field epidemiology and whole genome sequencing to identify tuberculosis transmission in a remote setting

Yukon Territory (YT) is a remote region in northern Canada with ongoing spread of tuberculosis (TB). To explore the utility of whole genome sequencing (WGS) for TB surveillance and monitoring in a setting with detailed contact tracing and interview data, we used a mixed-methods approach. Our analysis included all culture-confirmed cases in YT (2005-2014) and incorporated data from 24-locus Mycobacterial Interspersed Repetitive Units-Variable Number of Tandem Repeats (MIRU-VNTR) genotyping, WGS and contact tracing. We compared field-based (contact investigation (CI) data + MIRU-VNTR) and genomic-based (WGS + MIRU-VNTR + basic case data) investigations to identify the most likely source of each person's TB and assessed the knowledge, attitudes and practices of programme personnel around genotyping and genomics using online, multiple-choice surveys (n = 4) and an in-person group interview (n = 5). Field- and genomics-based approaches agreed for 26 of 32 (81%) cases on likely location of TB acquisition. There was less agreement in the identification of specific source cases (13/22 or 59% of cases). Single-locus MIRU-VNTR variants and limited genetic diversity complicated the analysis. Qualitative data indicated that participants viewed genomic epidemiology as a useful tool to streamline investigations, particularly in differentiating latent TB reactivation from the recent transmission. Based on this, genomic data could be used to enhance CIs, focus resources, target interventions and aid in TB programme evaluation.

New insights into the mycobacterial PE and PPE proteins provide a framework for future research

The PE and PPE proteins of Mycobacterium tuberculosis have been studied with great interest since their discovery. Named after the conserved proline (P) and glutamic acid (E) residues in their N-terminal domains, these proteins are postulated to perform wide-ranging roles in virulence and immune modulation. However, technical challenges in studying these proteins and their encoding genes have hampered the elucidation of molecular mechanisms and leave many open questions regarding the biological functions mediated by these proteins. Here, I review the shared and unique characteristics of PE and PPE proteins from a molecular perspective linking this information to their functions in mycobacterial virulence. I discuss how the different subgroups (PE_PGRS, PPE-PPW, PPE-SVP and PPE-MPTR) are defined and why this classification of paramount importance to understand the PE and PPE proteins as individuals and or groups. The goal of this MicroReview is to summarize and structure the existing information on this gene family into a simplified framework of thinking about PE and PPE proteins and genes. Thereby, I hope to provide helpful starting points in studying these genes and proteins for researchers with different backgrounds. This has particular implications for the design and monitoring of novel vaccine candidates and in understanding the evolution of the M. tuberculosis complex.

Complete nontuberculous mycobacteria whole genomes using an optimized DNA extraction protocol for long-read sequencing

Background

Nontuberculous mycobacteria (NTM) are a major cause of pulmonary and systemic disease in at-risk populations. Gaps in knowledge about transmission patterns, evolution, and pathogenicity during infection have prompted a recent surge in genomic NTM research. Increased availability and affordability of whole genome sequencing (WGS) techniques provide new opportunities to sequence and construct complete bacterial genomes faster and at a lower cost. However, extracting large quantities of pure genomic DNA is particularly challenging with NTM due to its slow growth and recalcitrant cell wall. Here we report a DNA extraction protocol that is optimized for long-read WGS of NTM, yielding large quantities of highly pure DNA with no additional clean-up steps.

Results

Our DNA extraction method was compared to 6 other methods with variations in timing of mechanical disruption and enzymatic digestion of the cell wall, quantity of matrix material, and reagents used in extraction and precipitation. We tested our optimized method on 38 clinical isolates from the M. avium and M. abscessus complexes, which yielded optimal quality and quantity measurements for Oxford Nanopore Technologies sequencing. We also present the efficient completion of circularized M. avium subspecies hominissuis genomes using our extraction technique and the long-read sequencing MinION platform, including the identification of a novel plasmid.

Conclusions

Our optimized extraction protocol and assembly pipeline was both sufficient and efficient for genome closure. We expect that our finely-tuned extraction method will prove to be a valuable tool in long-read sequencing and completion of mycobacterial genomes going forward. Utilization of comprehensive, long-read based approaches will advance the understanding evolution and pathogenicity of NTM infections.

Diversity of Mycobacterium tuberculosis in the Middle Fly District of Western Province, Papua New Guinea: microbead-based spoligotyping using DNA from Ziehl-Neelsen-stained microscopy preparations

Tuberculosis remains the world's leading cause of death from an infectious agent, and is a serious health problem in Papua New Guinea (PNG) with an estimated 36,000 new cases each year. This study describes the genetic diversity of Mycobacterium tuberculosis among tuberculosis patients in the Balimo/Bamu region in the Middle Fly District of Western Province in PNG, and investigates rifampicin resistance-associated mutations. Archived Ziehl-Neelsen-stained sputum smears were used to conduct microbead-based spoligotyping and assess genotypic resistance. Among the 162 samples included, 80 (49.4%) generated spoligotyping patterns (n = 23), belonging predominantly to the L2 Lineage (44%) and the L4 Lineage (30%). This is consistent with what has been found in other PNG regions geographically distant from Middle Fly District of Western Province, but is different from neighbouring South-East Asian countries. Rifampicin resistance was identified in 7.8% of the successfully sequenced samples, with all resistant samples belonging to the L2/Beijing Lineage. A high prevalence of mixed L2/L4 profiles was suggestive of polyclonal infection in the region, although this would need to be confirmed. The method described here could be a game-changer in resource-limited countries where large numbers of archived smear slides could be used for retrospective (and prospective) studies of M. tuberculosis genetic epidemiology.

Whole genome sequencing, analyses of drug resistance-conferring mutations, and correlation with transmission of Mycobacterium tuberculosis carrying katG-S315T in Hanoi, Vietnam

Drug-resistant tuberculosis (TB) is a serious global problem, and pathogen factors involved in the transmission of isoniazid (INH)-resistant TB have not been fully investigated. We performed whole genome sequencing of 332 clinical Mycobacterium tuberculosis (Mtb) isolates collected from patients newly diagnosed with smear-positive pulmonary TB in Hanoi, Vietnam. Using a bacterial genome-wide approach based on linear mixed models, we investigated the associations between 31-bp k-mers and clustered strains harboring katG-S315T, a major INH-resistance mutation in the present cohort and in the second panel previously published in South Africa. Five statistically significant genes, namely, PPE18/19, gid, emrB, Rv1588c, and pncA, were shared by the two panels. We further identified variants of the genes responsible for these k-mers, which are relevant to the spread of INH-resistant strains. Phylogenetic convergence test showed that variants relevant to PPE46/47-like chimeric genes were significantly associated with the same phenotype in Hanoi. The associations were further confirmed after adjustment for the confounders. These findings suggest that genomic variations of the pathogen facilitate the expansion of INH-resistance TB, at least in part, and our study provides a new insight into the mechanisms by which drug-resistant Mtb maintains fitness and spreads in Asia and Africa.

Completing Circular Bacterial Genomes With Assembly Complexity by Using a Sampling Strategy From a Single MinION Run With Barcoding

The Oxford Nanopore MinION is an affordable and portable DNA sequencer that can produce very long reads (tens of kilobase pairs), which enable de novo bacterial genome assembly. Although many algorithms and tools have been developed for base calling, read mapping, de novo assembly, and polishing, an automated pipeline is not available for one-stop analysis for circular bacterial genome reconstruction. In this paper, we present the pipeline CCBGpipe for completing circular bacterial genomes. Raw current signals are demultiplexed and base called to generate sequencing data. Sequencing reads are de novo assembled several times by using a sampling strategy to produce circular contigs that have a sequence in common between their start and end. The circular contigs are polished by using raw signals and sequencing reads; then, duplicated sequences are removed to form a linear representation of circular sequences. The circularized contigs are finally rearranged to start at the start position of dnaA/repA or a replication origin based on the GC skew. CCBGpipe implemented in Python is available at https://github.com/jade-nhri/CCBGpipe. Using sequencing data produced from a single MinION run, we obtained 48 circular sequences, comprising 12 chromosomes and 36 plasmids of 12 bacteria, including Acinetobacter nosocomialis, Acinetobacter pittii, and Staphylococcus aureus. With adequate quantities of sequencing reads (80×), CCBGpipe can provide a complete and automated assembly of circular bacterial genomes.

Applications of sequencing technology in clinical microbial infection

Infectious diseases are a type of disease caused by pathogenic microorganisms. Although the discovery of antibiotics changed the treatment of infectious diseases and reduced the mortality of bacterial infections, resistant bacterial strains have emerged. Anti‐infective therapy based on aetiological evidence is the gold standard for clinical treatment, but the time lag and low positive culture rate of traditional methods of pathogen diagnosis leads to relative difficulty in obtaining the evidence of pathogens. Compared with traditional methods of pathogenic diagnosis, next‐generation and third‐generation sequencing technologies have many advantages in the detection of pathogenic microorganisms. In this review, we mainly introduce recent progress in research on pathogenic diagnostic technology and the applications of sequencing technology in the diagnosis of pathogenic microorganisms. This review provides new insights into the application of sequencing technology in the clinical diagnosis of microorganisms.

Structural variants identified by Oxford Nanopore PromethION sequencing of the human genome

We sequenced the genome of the Yoruban reference individual NA19240 on the long-read sequencing platform Oxford Nanopore PromethION for evaluation and benchmarking of recently published aligners and germline structural variant calling tools, as well as a comparison with the performance of structural variant calling from short-read sequencing data. The structural variant caller Sniffles after NGMLR or minimap2 alignment provides the most accurate results, but additional confidence or sensitivity can be obtained by a combination of multiple variant callers. Sensitive and fast results can be obtained by minimap2 for alignment and a combination of Sniffles and SVIM for variant identification. We describe a scalable workflow for identification, annotation, and characterization of tens of thousands of structural variants from long-read genome sequencing of an individual or population. By discussing the results of this well-characterized reference individual, we provide an approximation of what can be expected in future long-read sequencing studies aiming for structural variant identification.

Severe Clinical Outcomes of Tuberculosis in Kharkiv Region, Ukraine, Are Associated with Beijing Strains of Mycobacterium tuberculosis

Genotypic variation in Beijing lineages of Mycobacterium tuberculosis (MTB), the causative agent of tuberculosis (TB), has been associated with hyper virulence and the spread of extensively and multiple drug (X/MDR) resistant MTB strains in Eastern Europe, Central Asia, and East Asia. The clinical outcomes of 215 new cases of TB among the population of the Kharkiv region of Eastern Ukraine were analyzed to uncover factors associated with severe infection. Infecting MTB strains were profiled by 5 locus exact tandem repeats (ETRs) and 15 locus mycobacterial interspersed repetitive unit-variable number tandem repeat (MIRU-VNTR) genotyping. Among diverse MTB genotypes discovered in Ukraine, the Beijing genotype (MIRU-VNTR 42425) was significantly associated with risk factors for severe outcomes of disease in the study population, including TB/HIV co-infection and treatment failure. Strain replacement (superinfection) was observed in 10 patients, suggesting repeated exposure to novel MTB strains in hospital or community settings. Inclusion of MTB genotyping data may identify at-risk patients and improve treatment adherence to prevent X/MDR development for effective public health response against tuberculosis in Ukraine.

NanoPipe-a web server for nanopore MinION sequencing data analysis

BACKGROUND

The fast-moving progress of the third-generation long-read sequencing technologies will soon bring the biological and medical sciences to a new era of research. Altogether, the technique and experimental procedures are becoming more straightforward and available to biologists from diverse fields, even without any profound experience in DNA sequencing. Thus, the introduction of the MinION device by Oxford Nanopore Technologies promises to "bring sequencing technology to the masses" and also allows quick and operative analysis in field studies. However, the convenience of this sequencing technology dramatically contrasts with the available analysis tools, which may significantly reduce enthusiasm of a "regular" user. To really bring the sequencing technology to every biologist, we need a set of user-friendly tools that can perform a powerful analysis in an automatic manner.

FINDINGS

NanoPipe was developed in consideration of the specifics of the MinION sequencing technologies, providing accordingly adjusted alignment parameters. The range of the target species/sequences for the alignment is not limited, and the descriptive usage page of NanoPipe helps a user to succeed with NanoPipe analysis. The results contain alignment statistics, consensus sequence, polymorphisms data, and visualization of the alignment. Several test cases are used to demonstrate the efficiency of the tool.

CONCLUSIONS

Freely available NanoPipe software allows effortless and reliable analysis of MinION sequencing data for experienced bioinformaticians, as well for wet-lab biologists with minimum bioinformatics knowledge. Moreover, for the latter group, we describe the basic algorithm necessary for MinION sequencing analysis from the first to last step.

Evaluation of strategies for the assembly of diverse bacterial genomes using MinION long-read sequencing

Background

Short-read sequencing technologies have made microbial genome sequencing cheap and accessible. However, closing genomes is often costly and assembling short reads from genomes that are repetitive and/or have extreme %GC content remains challenging. Long-read, single-molecule sequencing technologies such as the Oxford Nanopore MinION have the potential to overcome these difficulties, although the best approach for harnessing their potential remains poorly evaluated.

Results

We sequenced nine bacterial genomes spanning a wide range of GC contents using Illumina MiSeq and Oxford Nanopore MinION sequencing technologies to determine the advantages of each approach, both individually and combined. Assemblies using only MiSeq reads were highly accurate but lacked contiguity, a deficiency that was partially overcome by adding MinION reads to these assemblies. Even more contiguous genome assemblies were generated by using MinION reads for initial assembly, but these assemblies were more error-prone and required further polishing. This was especially pronounced when Illumina libraries were biased, as was the case for our strains with both high and low GC content. Increased genome contiguity dramatically improved the annotation of insertion sequences and secondary metabolite biosynthetic gene clusters, likely because long-reads can disambiguate these highly repetitive but biologically important genomic regions.

Conclusions

Genome assembly using short-reads is challenged by repetitive sequences and extreme GC contents. Our results indicate that these difficulties can be largely overcome by using single-molecule, long-read sequencing technologies such as the Oxford Nanopore MinION. Using MinION reads for assembly followed by polishing with Illumina reads generated the most contiguous genomes with sufficient accuracy to enable the accurate annotation of important but difficult to sequence genomic features such as insertion sequences and secondary metabolite biosynthetic gene clusters. The combination of Oxford Nanopore and Illumina sequencing can therefore cost-effectively advance studies of microbial evolution and genome-driven drug discovery.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5381-7) contains supplementary material, which is available to authorized users.