Control of Artifactual Variation in Reported Intersample Relatedness during Clinical Use of a Mycobacterium tuberculosis Sequencing Pipeline

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.

Pangenome databases improve host removal and mycobacteria classification from clinical metagenomic data

BACKGROUND

Culture-free real-time sequencing of clinical metagenomic samples promises both rapid pathogen detection and antimicrobial resistance profiling. However, this approach introduces the risk of patient DNA leakage. To mitigate this risk, we need near-comprehensive removal of human DNA sequences at the point of sequencing, typically involving the use of resource-constrained devices. Existing benchmarks have largely focused on the use of standardized databases and largely ignored the computational requirements of depletion pipelines as well as the impact of human genome diversity.

RESULTS

We benchmarked host removal pipelines on simulated and artificial real Illumina and Nanopore metagenomic samples. We found that construction of a custom kraken database containing diverse human genomes results in the best balance of accuracy and computational resource usage. In addition, we benchmarked pipelines using kraken and minimap2 for taxonomic classification of Mycobacterium reads using standard and custom databases. With a database representative of the Mycobacterium genus, both tools obtained improved specificity and sensitivity, compared to the standard databases for classification of Mycobacterium tuberculosis. Computational efficiency of these custom databases was superior to most standard approaches, allowing them to be executed on a laptop device.

CONCLUSIONS

Customized pangenome databases provide the best balance of accuracy and computational efficiency when compared to standard databases for the task of human read removal and M. tuberculosis read classification from metagenomic samples. Such databases allow for execution on a laptop, without sacrificing accuracy, an especially important consideration in low-resource settings. We make all customized databases and pipelines freely available.

Mycobacterium tuberculosis transmission in Birmingham, UK, 2009-19: An observational study

Background

Over 10-years of whole-genome sequencing (WGS) of Mycobacterium tuberculosis in Birmingham presents an opportunity to explore epidemiological trends and risk factors for transmission in new detail.

Methods

Between 1st January 2009 and 15th June 2019, we obtained the first WGS isolate from every patient resident in a postcode district covered by Birmingham's centralised tuberculosis service. Data on patients' sex, country of birth, social risk-factors, anatomical locus of disease, and strain lineage were collected. Poisson harmonic regression was used to assess seasonal variation in case load and a mixed-effects multivariable Cox proportionate hazards model was used to assess risk factors for a future case arising in clusters defined by a 5 single nucleotide polymorphism (SNP) threshold, and by 12 SNPs in a sensitivity analysis.

Findings

511/1653 (31%) patients were genomically clustered with another. A seasonal variation in diagnoses was observed, peaking in spring, but only among clustered cases. Risk-factors for a future clustered case included UK-birth (aHR=2·03 (95%CI 1·35-3·04), p < 0·001), infectious (pulmonary/laryngeal/miliary) tuberculosis (aHR=3·08 (95%CI 1·98-4·78), p < 0·001), and M. tuberculosis lineage 3 (aHR=1·91 (95%CI 1·03-3·56), p = 0·041) and 4 (aHR=2·27 (95%CI 1·21-4·26), p = 0·011), vs. lineage 1. Similar results pertained to 12 SNP clusters, for which social risk-factors were also significant (aHR 1·72 (95%CI 1·02-2·93), p = 0·044). There was marked heterogeneity in transmission patterns between postcode districts.

Interpretation

There is seasonal variation in the diagnosis of genomically clustered, but not non-clustered, cases. Risk factors for clustering include UK-birth, infectious forms of tuberculosis, and infection with lineage 3 or 4.

Funding

Wellcome Trust, MRC, UKHSA.

M. tuberculosis microvariation is common and is associated with transmission: analysis of three years prospective universal sequencing in England

BACKGROUND

The prevalence, association with disease status, and public health impact of infection with mixtures of M. tuberculosis strains is unclear, in part due to limitations of existing methods for detecting mixed infections.

METHODS

We developed an algorithm to identify mixtures of M. tuberculosis strains using next generation sequencing data, assessing performance using simulated sequences. We identified mixed M. tuberculosis strains when there was at least one mixed nucleotide position, and where both the mixture's components were present in similar isolates from other individuals, compatible with transmission of the component strains. We determined risk factors for mixed infection among isolations of M. tuberculosis in England using logistic regression. We used survival analyses to assess the association between mixed infection and putative transmission.

FINDINGS

6,560 isolations of TB were successfully sequenced in England 2016-2018. Of 3,691 (56%) specimens for which similar sequences had been isolated from at least two other individuals, 341 (9.2%) were mixed. Mixed infection was more common in lineages other than Lineage 4. Among the 1,823 individuals with pulmonary infection with Lineage 4 M. tuberculosis, mixed infection was associated with significantly increased risk of subsequent isolation of closely related organisms from a different individual (HR 1.43, 95% CI 1.05,1.94), indicative of transmission.

INTERPRETATION

Mixtures of transmissible strains occur in at least 5% of tuberculosis infections in England; when present in pulmonary disease, such mixtures are associated with an increased risk of tuberculosis transmission.

FUNDING

Public Health England; NIHR Health Protection Research Units; European Union.

DNA Thermo-Protection Facilitates Whole-Genome Sequencing of Mycobacteria Direct from Clinical Samples

Mycobacterium tuberculosis is the leading cause of death from bacterial infection. Improved rapid diagnosis and antimicrobial resistance determination, such as by whole-genome sequencing, are required. Our aim was to develop a simple, low-cost method of preparing DNA for sequencing direct from M. tuberculosis-positive clinical samples (without culture). Simultaneous sputum liquefaction, bacteria heat inactivation (99°C/30 min), and enrichment for mycobacteria DNA were achieved using an equal volume of thermo-protection buffer (4 M KCl, 0.

Mycobacterium tuberculosis is the leading cause of death from bacterial infection. Improved rapid diagnosis and antimicrobial resistance determination, such as by whole-genome sequencing, are required. Our aim was to develop a simple, low-cost method of preparing DNA for sequencing direct from M. tuberculosis-positive clinical samples (without culture). Simultaneous sputum liquefaction, bacteria heat inactivation (99°C/30 min), and enrichment for mycobacteria DNA were achieved using an equal volume of thermo-protection buffer (4 M KCl, 0.05 M HEPES buffer, pH 7.5, 0.1% dithiothreitol [DTT]). The buffer emulated intracellular conditions found in hyperthermophiles, thus protecting DNA from rapid thermodegradation, which renders it a poor template for sequencing. Initial validation experiments employed mycobacteria DNA, either extracted or intracellular. Next, mock clinical samples (infection-negative human sputum spiked with 0 to 10⁵Mycobacterium bovis BCG cells/ml) underwent liquefaction in thermo-protection buffer and heat inactivation. DNA was extracted and sequenced. Human DNA degraded faster than mycobacteria DNA, resulting in target enrichment. Four replicate experiments achieved M. tuberculosis detection at 10¹ BCG cells/ml, with 31 to 59 M. tuberculosis complex reads. Maximal genome coverage (>97% at 5× depth) occurred at 10⁴ BCG cells/ml; >91% coverage (1× depth) occurred at 10³ BCG cells/ml. Final validation employed M. tuberculosis-positive clinical samples (n = 20), revealing that initial sample volumes of ≥1 ml typically yielded higher mean depths of M. tuberculosis genome coverage, with an overall range of 0.55 to 81.02. A mean depth of 3 gave >96% 1-fold tuberculosis (TB) genome coverage (in 15/20 clinical samples). A mean depth of 15 achieved >99% 5-fold genome coverage (in 9/20 clinical samples). In summary, direct-from-sample sequencing of M. tuberculosis genomes was facilitated by a low-cost thermo-protection buffer.

High precision Neisseria gonorrhoeae variant and antimicrobial resistance calling from metagenomic Nanopore sequencing

The rise of antimicrobial-resistant Neisseria gonorrhoeae is a significant public health concern. Against this background, rapid culture-independent diagnostics may allow targeted treatment and prevent onward transmission. We have previously shown metagenomic sequencing of urine samples from men with urethral gonorrhea can recover near-complete N. gonorrhoeae genomes. However, disentangling the N. gonorrhoeae genome from metagenomic samples and robustly identifying antimicrobial resistance determinants from error-prone Nanopore sequencing is a substantial bioinformatics challenge. Here, we show an N. gonorrhoeae diagnostic workflow for analysis of metagenomic sequencing data obtained from clinical samples using R9.4.1 Nanopore sequencing. We compared results from simulated and clinical infections with data from known reference strains and Illumina sequencing of isolates cultured from the same patients. We evaluated three Nanopore variant callers and developed a random forest classifier to filter called SNPs. Clair was the most suitable variant caller after SNP filtering. A minimum depth of 20× reads was required to confidently identify resistant determinants over the entire genome. Our findings show that metagenomic Nanopore sequencing can provide reliable diagnostic information in N. gonorrhoeae infection.

Role and value of whole genome sequencing in studying tuberculosis transmission

BACKGROUND

Tuberculosis (TB) remains a serious public health threat worldwide. Theoretically ultimate resolution of whole genome sequencing (WGS) for Mycobacterium tuberculosis complex (MTBC) strain classification makes this technology very attractive for epidemiological investigations.

OBJECTIVES

To summarize the evidence available in peer-reviewed publications on the role and place of WGS in detection of TB transmission.

SOURCES

A total of 69 peer-reviewed publications identified in Pubmed database.

CONTENT

Evidence from >30 publications suggests that a cut-off value of fewer than six single nucleotide polymorphisms between strains efficiently excludes cases that are not the result of recent transmission and could be used for the identification of drug-sensitive isolates involved in direct human-to-human TB transmission. Sensitivity of WGS to identify epidemiologically linked isolates is high, reaching 100% in eight studies with specificity (17%-95%) highly dependent on the settings. Drug resistance and specific phylogenetic lineages may be associated with accelerated mutation rates affecting genetic distances. WGS can be potentially used to distinguish between true relapses and re-infections but in high-incidence low-diversity settings this would require consideration of epidemiological links and minority alleles. Data from four studies looking into within-host diversity highlight a need for developing criteria for acceptance or rejection of WGS relatedness results depending on the proportion of minority alleles.

IMPLICATIONS

WGS will potentially allow for more targeted public health actions preventing unnecessary investigations of false clusters. Consensus on standardization of raw data quality control processing criteria, analytical pipelines and reporting language is yet to be reached.

A Quantitative Evaluation of MIRU-VNTR Typing Against Whole-Genome Sequencing for Identifying Mycobacterium tuberculosis Transmission: A Prospective Observational Cohort Study

Background

Mycobacterial Interspersed Repetitive Unit-Variable Number Tandem Repeat (MIRU-VNTR) typing is widely used in high-income countries to determine Mycobacterium tuberculosis relatedness. Whole-genome sequencing (WGS) is known to deliver greater specificity, but no quantitative prospective comparison has yet been undertaken.

Methods

We studied isolates from the English Midlands, sampled consecutively between 1 January 2012 and 31 December 2015. In addition to routinely performed MIRU-VNTR typing, DNA was extracted from liquid cultures and sequenced using Illumina technology. Demographic and epidemiological data for the relevant patients were extracted from the Enhanced Tuberculosis Surveillance system run by Public Health England. Closely related samples, defined using a threshold of five single nucleotide variants (SNVs), were compared to samples with identical MIRU-VNTR profiles, to samples from individuals with shared epidemiological risk factors, and to those with both characteristics.

Findings

1999 patients were identified for whom at least one M. tuberculosis isolate had been MIRU-VNTR typed and sequenced. Comparing epidemiological risk factors with close genetic relatedness, only co-residence had a positive predictive value of over 5%. Excluding co-resident individuals, 18.6% of patients with identical MIRU-VNTR profiles were within 5 SNVs. Where patients also shared social risk factors and ethnic group, this rose to 48%. Only 8% of MIRU-VNTR linked pairs in lineage 1 were within 5 SNV, compared to 31% in lineage 4.

Interpretation

In the setting studied, this molecular epidemiological study shows MIRU-VNTR typing and epidemiological risk factors are poorly predictive of close genomic relatedness, assessed by SNV. MIRU-VNTR performance varies markedly by lineage.

Funding

Public Health England, Health Innovation Challenge Fund, NIHR Health Protection Research Unit Oxford, NIHR Oxford Biomedical Research Centre.

•

M. tuberculosis relatedness was measured by MIRU-VNTR typing and single nucleotide variants (SNV).

•

In consecutive isolates in England, only 20% of samples with identical MIRU-VNTR profiles were within 5 SNVs.

•

MIRU-VNTR is more predictive of close relatedness in lineage 4 (Euro-American) lineages than in other lineages.