The MAGMA pipeline for comprehensive genomic analyses of clinical Mycobacterium tuberculosis samples

Phenotypic and genotypic resistance to bedaquiline in patients with multi-drug-resistant tuberculosis-experiences from Armenia

Risk factors for baseline bedaquiline (BDQ) resistance, amplification during treatment, and correlations with treatment outcomes are not fully understood. This cohort included Armenian patients with multidrug-resistant TB predominantly fluoroquinolone-resistant enrolled between 2013 and 2015 in a BDQ compassionate use program. BDQ resistance at baseline and during treatment was assessed using MGIT (pDSTMGIT), minimal inhibitory concentration in 7H11 (MIC7H11), and whole-genome sequencing. Risk factors, such as treatment effectiveness or stage of the disease, were analyzed for association with baseline BDQ resistance, acquired BDQ resistance, and treatment outcome. Among 39 patients, baseline BDQ resistance was 6% (2/33) by pDSTMGIT and 7% (2/29) by MIC7H11. All four baseline isolates with an Rv0678 mutation were phenotypically resistant. During treatment, 48% of the patients acquired BDQ resistance by pDSTMGIT, and 52% acquired mutations at various frequencies (97% in Rv0678). None of the factors significantly contributed to baseline or acquired BDQ resistance. Unfavorable treatment outcome (41%) was more frequent in the presence of acquired Rv0678 mutations [odds ratio (OR) 132, 95% confidence interval (CI) 7.43, 2375], phenotypic BDQ resistance (OR 176, 95% CI 6.48, 2423), or MIC increase above or below the critical concentration (both OR 84.3, 95% CI 2.93, 2423) during treatment. For these highly treatment-experienced patients, low baseline prevalence but high incidence of acquired BDQ resistance was observed. Acquisition of mutations in BDQ candidate resistance genes, regardless of their frequency, or increased MICs during treatment, even below the critical concentration, should be seen as a warning sign of resistance amplification and increased risk of unfavorable treatment outcome.

Comprehensive genomic surveillance reveals transmission profiles of extensively drug-resistant tuberculosis cases in Pará, Brazil

Bedaquiline, an antimicrobial used to treat drug-resistant tuberculosis (DR-TB), was introduced in Brazil in October 2021. Monitoring the emergence and transmission of DR-TB is crucial for implementing public health to control the spread of DR strains of Mycobacterium tuberculosis. To measure its impact on the multi-drug treatment scheme in the state of Pará, we aimed to conduct genomic surveillance of DR-TB after bedaquiline was introduced in Brazil. Individuals treated for DR-TB between October 2021 and December 2022, in the reference hospital to treat DR-TB cases from the state of Pará, were included in the study. Clinical and bacteriological information was obtained from the National Laboratory Management Environment and the Special TB Treatment Information System. Genomic DNA was extracted from bacterial cultures performed at the Pará Central Laboratory (LACEN-PA). Whole-genome sequencing (WGS) was obtained using Illumina Nextera-XT and NextSeq 550 and genomes were analyzed using the MAGMA and TB-Profiler pipelines interpreted according to the World Health Organization (WHO) mutations catalog 2nd edition. Geoprocessing was performed based on the patient's residences. Cutoffs of 5-12 single nucleotide polymorphisms (SNPs) were used for transmission analysis. From the 103 patients reported as DR-TB, viable cultures were obtained from 67. Forty isolates were selected randomly for WGS. Among these, a mixed infection of M. tuberculosis L1 and L4 and a co-infection of M. tuberculosis and Mycobacterium chelonae were observed. The genotypic drug susceptibility profile of TB stains (39/40) was as follows: sensitive (1/2, 5%), rifampicin mono-resistant (RR) (4/10%), isoniazid mono-resistant (1/2%), multidrug-resistant (MDR) (21/52%), extensively drug-resistant (XDR) (3/7%), pre-XDR (8/20%), and other (1/2%). Among the 38 isolates of M. tuberculosis strains without mixed infection, using a cutoff of 12 SNPs and suggestive of recent TB transmission, 14 (37%) were grouped into five clusters (C1-C5) and included RR (C5), MDR (C3, C4, C5), pre-XDR, and XDR (C2) strains. We recommend greater attention from the regional public health authorities to detect and track resistance to new drugs, especially in areas with pre-XDR and XDR cases. This is the first report on the detection and transmission of XDR-TB in Pará, Brazil, after the recent re-definition of XDR-TB by the WHO in 2021.

The chosen few: Mycobacterium tuberculosis isolates for IMPAc-TB

The three programs that make up the Immune Mechanisms of Protection Against Mycobacterium tuberculosis Centers (IMPAc-TB) had to prioritize and select strains to be leveraged for this work. The CASCADE team based at Seattle Children's Research Institute are leveraging M.tb H37Rv, M.tb CDC1551, and M.tb SA161. The HI-IMPACT team based at Harvard T.H. Chan School of Public Health, Boston, have selected M.tb Erdman as well as a novel clinical isolate recently characterized during a longitudinal study in Peru. The PHOENIX team also based at Seattle Children's Research Institute have selected M.tb HN878 and M.tb Erdman as their isolates of choice. Here, we describe original source isolation, genomic references, key virulence characteristics, and relevant tools that make these isolates attractive for use. The global context for M.tb lineage 2 and 4 selection is reviewed including what is known about their relative abundance and acquisition of drug resistance. Host-pathogen interactions seem driven by genomic differences on each side, and these play an important role in pathogenesis and immunity. The few M.tb strains chosen for this work do not reflect the vast genomic diversity within this species. They do, however, provide specific virulence, pathology, and growth kinetics of interest to the consortium. The strains selected should not be considered as "representative" of the growing available array of M.tb isolates, but rather tools that are being used to address key outstanding questions in the field.

Exploring the potential of Oxford Nanopore Technologies sequencing for Mycobacterium tuberculosis sequencing: An assessment of R10 flowcells and V14 chemistry

Oxford Nanopore Technologies (ONT) sequencing is a promising technology. We assessed the performance of the new ONT R10 flowcells and V14 rapid sequencing chemistry for Mtb whole genome sequencing of Mycobacterium tuberculosis (Mtb) DNA extracted from clinical primary liquid cultures (CPLCs). Using the recommended protocols for MinION Mk1C, R10.4.1 MinION flowcells, and the ONT Rapid Sequencing Kit V14 on six CPLC samples, we obtained a pooled library yield of 10.9 ng/μl, generated 1.94 Gb of sequenced bases and 214k reads after 48h in a first sequencing run. Only half (49%) of all generated reads met the Phred Quality score threshold (>8). To assess if the low data output and sequence quality were due to impurities present in DNA extracted directly from CPLCs, we added a pre-library preparation bead-clean-up step and included purified DNA obtained from an Mtb subculture as a control sample in a second sequencing run. The library yield for DNA extracted from four CPLCs and one Mtb subculture (control) was similar (10.0 ng/μl), 2.38 Gb of bases and 822k reads were produced. The quality was slightly better with 66% of the produced reads having a Phred Quality >8. A third run of DNA from six CPLCs with bead clean-up pre-processing produced a low library yield (±1 Gb of bases, 166k reads) of low quality (51% of reads with a Phred Quality score >8). A median depth of coverage above 10× was only achieved for five of 17 (29%) sequenced libraries. Compared to Illumina WGS of the same samples, accurate lineage predictions and full drug resistance profiles from the generated ONT data could not be determined by TBProfiler. Further optimization of the V14 ONT rapid sequencing chemistry and library preparation protocol is needed for clinical Mtb WGS applications.

Droplet based whole genome amplification for sequencing minute amounts of purified Mycobacterium tuberculosis DNA

Implementation of whole genome sequencing (WGS) for patient care is hindered by limited Mycobacterium tuberculosis (Mtb) in clinical specimens and slow Mtb growth. We evaluated droplet multiple displacement amplification (dMDA) for amplification of minute amounts of Mtb DNA to enable WGS as an alternative to other Mtb enrichment methods. Purified genomic Mtb-DNA (0.1, 0.5, 1, and 5 pg) was encapsulated and amplified using the Samplix Xdrop-instrument and sequenced alongside a control sample using standard Illumina protocols followed by MAGMA-analysis. The control and 5 pg input dMDA samples underwent nanopore sequencing followed by Nanoseq and TB-profiler analysis. dMDA generated 105-2400 ng DNA from the 0.1-5 pg input DNA, respectively. Followed by Illumina WGS, dMDA raised mean sequencing depth from 7 × for 0.1 pg input DNA to ≥ 60 × for 5 pg input and the control sample. Bioinformatic analysis revealed a high number of false positive and false negative variants when amplifying ≤ 0.5 pg input DNA. Nanopore sequencing of the 5 pg dMDA sample presented excellent coverage depth, breadth, and accurate strain characterization, albeit elevated false positive and false negative variants compared to Illumina-sequenced dMDA sample with identical Mtb DNA input. dMDA coupled with Illumina WGS for samples with ≥ 5 pg purified Mtb DNA, equating to approximately 1000 copies of the Mtb genome, offers precision for drug resistance, phylogeny, and transmission insights.

Systematic Analysis of the Relationship Between Elevated Zinc and Epilepsy

Previous studies have indicated a potential relationship between zinc and epilepsy. The aim of this study is to investigate the causal relationship between zinc, zinc-dependent carbonic anhydrase, and gray matter volume in brain regions enriched with zinc and epilepsy, as well as explore the possible mechanisms by which zinc contributes to epilepsy. First, this study assessed the risk causality between zinc, carbonic anhydrase, and gray matter volume alterations in zinc-enriched brain regions and various subtypes of epilepsy based on Two-sample Mendelian randomization analysis. And then, this study conducted GO/KEGG analysis based on colocalization analysis, MAGMA analysis, lasso regression, random forest model, and XGBoost model. The results of Mendelian randomization analyses showed a causal relationship between zinc, carbonic anhydrase-4, and generalized epilepsy (p = 0.044 , p = 0.010). Additionally, carbonic anhydrase-1 and gray matter volume of the caudate nucleus were found to be associated with epilepsy and focal epilepsy (p = 0.014, p = 0.003 and p = 0.022, p = 0.009). A colocalization relationship was found between epilepsy and focal epilepsy (PP.H4.abf = 97.7e - 2). Meanwhile, the MAGMA analysis indicated that SNPs associated with epilepsy and focal epilepsy were functionally localized to zinc-finger-protein-related genes (p < 1.0e - 5). The genes associated with focal epilepsy were found to have a molecular function of zinc ion binding (FDR = 2.3e - 6). After the onset of epilepsy, the function of the gene whose expression changed in the rats with focal epilepsy was enriched in the biological process of vascular response (FDR = 4.0e - 5). These results revealed mechanism of the increased risk of epilepsy caused by elevated zinc may be related to the increase of zinc ion-dependent carbonic anhydrase or the increase of the volume of zinc-rich caudate gray matter.

Characteristic SNPs defining the major multidrug-resistant Mycobacterium tuberculosis clusters identified by EuSeqMyTB to support routine surveillance, EU/EEA, 2017 to 2019

BackgroundThe EUSeqMyTB project, conducted in 2020, used whole genome sequencing (WGS) for surveillance of drug-resistant Mycobacterium tuberculosis in the European Union/European Economic Area (EU/EEA) and identified 56 internationally clustered multidrug-resistant (MDR) tuberculosis (TB) clones.AimWe aimed to define and establish a rapid and computationally simple screening method to identify probable members of the main cross-border MDR-TB clusters in WGS data to facilitate their identification and track their future spread.MethodsWe screened 34 of the larger cross-border clusters identified in the EuSeqMyTB pilot study (2017-19) for characteristic single nucleotide polymorphism (SNP) signatures that could identify and define members of each cluster. We also linked this analysis with published clusters identified in previous studies and identified more distant genetic relationships between some of the current clusters.ResultsA panel of 30 characteristic SNPs is presented that can be used as an initial (routine) screen for members of each cluster. For four of the clusters, no unique defining SNP could be identified; three of these are closely related (within approximately 20 SNPs) to one or more other clusters and likely represent a single established MDR-TB clade composed of multiple recent subclusters derived from the previously described ECDC0002 cluster.ConclusionThe identified SNP signatures can be integrated into routine pipelines and contribute to the more effective monitoring, rapid and widespread screening for TB. This SNP panel will also support accurate communication between laboratories about previously identified internationally transmitted MDR-TB genotypes.