Genomic Sequencing from Sputum for Tuberculosis Disease Diagnosis, Lineage Determination, and Drug Susceptibility Prediction.

Universal access to drug susceptibility testing for newly diagnosed tuberculosis patients is recommended. Access to culture-based diagnostics remains limited, and targeted molecular assays are vulnerable to emerging resistance mutations. Improved protocols for direct-from-sputum Mycobacterium tuberculosis sequencing would accelerate access to comprehensive drug susceptibility testing and molecular typing. We assessed a thermo-protection buffer-based direct-from-sample M. tuberculosis whole-genome sequencing protocol. We prospectively analyzed 60 acid-fast bacilli smear-positive clinical sputum samples in India and Madagascar. A diversity of semiquantitative smear positivity-level samples were included. Sequencing was performed using Illumina and MinION (monoplex and multiplex) technologies. We measured the impact of bacterial inoculum and sequencing platforms on genomic read depth, drug susceptibility prediction performance, and typing accuracy. M. tuberculosis was identified by direct sputum sequencing in 45/51 samples using Illumina, 34/38 were identified using MinION-monoplex sequencing, and 20/24 were identified using MinION-multiplex sequencing. The fraction of M. tuberculosis reads from MinION sequencing was lower than from Illumina, but monoplexing grade 3+ samples on MinION produced higher read depth than Illumina (P < 0.05) and MinION multiplexing (P < 0.01). No significant differences in sensitivity and specificity of drug susceptibility predictions were seen across sequencing modalities or within each technology when stratified by smear grade. Illumina sequencing from sputum accurately identified 1/8 (rifampin) and 6/12 (isoniazid) resistant samples, compared to 2/3 (rifampin) and 3/6 (isoniazid) accurately identified with Nanopore monoplex. Lineage agreement levels between direct and culture-based sequencing were 85% (MinION-monoplex), 88% (Illumina), and 100% (MinION-multiplex). M. tuberculosis direct-from-sample whole-genome sequencing remains challenging. Improved and affordable sample treatment protocols are needed prior to clinical deployment.

Impact of COVID-19 vaccination on transmission risk of breakthrough infections: Lessons from adapted N95 mask sampling for emerging variants and interventions.

This study used an adapted N95 mask sampling to understand the effect of COVID-19 vaccination in the context of circulating variants on infected individuals to emit the virus into the air, a key risk factor of transmission. Mask, swab, and blood samples were collected from 92 COVID-19 patients vaccinated (Covishield/COVAXIN-partial/fully) or unvaccinated between July and September 2021 during the Delta-dominated period in Mumbai. Mask/swab samples were analyzed by reverse transcription polymerase chain reaction for viral RNA. Blood was evaluated for SARS-CoV-2 anti-spike and nucleocapsid antibody responses. At <48 h of diagnosis, 93% of the patients emitted detectable viral RNA, with 40% emitting >1000 copies in 30 min (high emitters). About 8% continued to be high emitters even after 8 days of symptom onset. No significant difference was observed in emission patterns between partial, full, and unvaccinated patients. However, when vaccinated patients were stratified based on spike protein neutralization and nucleocapsid immunoglobulin G (IgG), the group with moderate/high neutralization showed a significantly lower proportion of high emitters and viral RNA copies than the group with no/low neutralization, which further reduced in the group having antinucleocapsid IgG. In conclusion, mask sampling showed that Delta infections were associated with greater virus emission in patients, which was significantly reduced only in vaccinated patients with moderate/high SARS-CoV-2 neutralization, especially with evidence of past infection. The study demonstrated that mask sampling could be useful for understanding the transmission risk of emerging variants, screening vaccine/booster candidates, and guiding control interventions.

Genomics of postvaccination SARS-CoV-2 infections during the Delta dominated second wave of COVID-19 pandemic, from Mumbai Metropolitan Region (MMR), India.

The present study was initiated to understand the proportion of predominant variants of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) in postvaccination infections during the Delta dominated second wave of coronavirus disease 2019 (COVID-19) in the Mumbai Metropolitan Region (MMR) in India and to understand any mutations selected in the postvaccination infections or showing association with any patient demographics. Samples were collected (n = 166) from severe/moderate/mild COVID-19 patients who were either vaccinated (COVISHIELD/COVAXIN-partial/fully vaccinated) or unvaccinated, from a city hospital and from home isolation patients in MMR. A total of 150 viral genomes were sequenced by Oxford Nanopore sequencing and the data of 136 viral genomes were analyzed for clade/lineage and for identifying mutations. The sequences belonged to three clades (21A, 21I, and 21J) and their lineage was identified as either Delta (B.1.617.2) or Delta+ (B.1.617.2 + K417N) or sub-lineages of Delta variant (AY.120/AY.38/AY.99). A total of 620 mutations were identified of which 10 mutations showed an increase in trend with time (May-October 2021). Associations of six mutations (two in spike, three in orf1a, and one in nucleocapsid) were shown with milder forms of the disease and one mutation (in orf1a) with partial vaccination status. The results indicate a trend toward reduction in disease severity as the wave progressed.

Deregulation of Genes Associated with Alternate Drug Resistance Mechanisms in Mycobacterium tuberculosis.

Alternate mechanisms of drug resistance involving intrinsic defense pathways play an important role in development of drug resistance. Deregulation of drug efflux, cellular metabolism, and DNA repair have been indicated to have effect on drug tolerance and persistence. Here we chose eight genes from these pathways to investigate their association with development of multidrug resistance (MDR). We generated mono drug resistant and MDR strains of rifampicin and isoniazid and examined the differential expression of genes belonging to efflux, DNA repair and cell wall lipid synthesis pathways. Rv1687c, recB, ppsD and embC genes showed significant (P <0.05) upregulation in mono-resistant (both rifampicin and isoniazid) as well as MDR strains. mmr showed significant upregulation with rifampicin resistance while Rv1457c showed significant upregulation only with mono-resistant strains. Highest expression change was observed with Rv1687c and ppsD. The study identified potential key genes that are significantly associated with development of drug resistance in vitro. These genes may help identify clinical strains predisposed to acquiring drug resistance in patients during the course of treatment or help in management of MDR forms of tuberculosis.

Basis of science policies for infectious disease challenges in India.

With increase in drug resistance and related challenges in infectious diseases globally, it has become more important for a country like India with a high population to develop strategies to deal with such challenges. The literature review was conducted using the Google search engine to explore the contemporary science policies in India (since 2012) which are designed for tackling infectious diseases challenges in the country. This review article gives an insight into the strengths and limitations of the basis of some of the contemporary science policies in India that are drafted and implemented to combat the challenges of infectious diseases. The new national plans for controlling infectious disease challenges with bold strategies for their rapid and effective implementation hold promise in India.

Same-Day Diagnostic and Surveillance Data for Tuberculosis via Whole-Genome Sequencing of Direct Respiratory Samples.

Routine full characterization of Mycobacterium tuberculosis is culture based, taking many weeks. Whole-genome sequencing (WGS) can generate antibiotic susceptibility profiles to inform treatment, augmented with strain information for global surveillance; such data could be transformative if provided at or near the point of care. We demonstrate a low-cost method of DNA extraction directly from patient samples for M. tuberculosis WGS. We initially evaluated the method by using the Illumina MiSeq sequencer (40 smear-positive respiratory samples obtained after routine clinical testing and 27 matched liquid cultures). M. tuberculosis was identified in all 39 samples from which DNA was successfully extracted. Sufficient data for antibiotic susceptibility prediction were obtained from 24 (62%) samples; all results were concordant with reference laboratory phenotypes. Phylogenetic placement was concordant between direct and cultured samples. With Illumina MiSeq/MiniSeq, the workflow from patient sample to results can be completed in 44/16 h at a reagent cost of £96/£198 per sample. We then employed a nonspecific PCR-based library preparation method for sequencing on an Oxford Nanopore Technologies MinION sequencer. We applied this to cultured Mycobacterium bovis strain BCG DNA and to combined culture-negative sputum DNA and BCG DNA. For flow cell version R9.4, the estimated turnaround time from patient to identification of BCG, detection of pyrazinamide resistance, and phylogenetic placement was 7.5 h, with full susceptibility results 5 h later. Antibiotic susceptibility predictions were fully concordant. A critical advantage of MinION is the ability to continue sequencing until sufficient coverage is obtained, providing a potential solution to the problem of variable amounts of M. tuberculosis DNA in direct samples.

Bedaquiline and clofazimine resistance in Mycobacterium tuberculosis: an in-vitro and in-silico data analysis.

BACKGROUND: Bedaquiline is a core drug for the treatment of multidrug-resistant tuberculosis; however, the understanding of resistance mechanisms is poor, which is hampering rapid molecular diagnostics. Some bedaquiline-resistant mutants are also cross-resistant to clofazimine. To decipher bedaquiline and clofazimine resistance determinants, we combined experimental evolution, protein modelling, genome sequencing, and phenotypic data. METHODS: For this in-vitro and in-silico data analysis, we used a novel in-vitro evolutionary model using subinhibitory drug concentrations to select bedaquiline-resistant and clofazimine-resistant mutants. We determined bedaquiline and clofazimine minimum inhibitory concentrations and did Illumina and PacBio sequencing to characterise selected mutants and establish a mutation catalogue. This catalogue also includes phenotypic and genotypic data of a global collection of more than 14 000 clinical Mycobacterium tuberculosis complex isolates, and publicly available data. We investigated variants implicated in bedaquiline resistance by protein modelling and dynamic simulations. FINDINGS: We discerned 265 genomic variants implicated in bedaquiline resistance, with 250 (94%) variants affecting the transcriptional repressor (Rv0678) of the MmpS5-MmpL5 efflux system. We identified 40 new variants in vitro, and a new bedaquiline resistance mechanism caused by a large-scale genomic rearrangement. Additionally, we identified in vitro 15 (7%) of 208 mutations found in clinical bedaquiline-resistant isolates. From our in-vitro work, we detected 14 (16%) of 88 mutations so far identified as being associated with clofazimine resistance and also seen in clinically resistant strains, and catalogued 35 new mutations. Structural modelling of Rv0678 showed four major mechanisms of bedaquiline resistance: impaired DNA binding, reduction in protein stability, disruption of protein dimerisation, and alteration in affinity for its fatty acid ligand. INTERPRETATION: Our findings advance the understanding of drug resistance mechanisms in M tuberculosis complex strains. We have established an extended mutation catalogue, comprising variants implicated in resistance and susceptibility to bedaquiline and clofazimine. Our data emphasise that genotypic testing can delineate clinical isolates with borderline phenotypes, which is essential for the design of effective treatments. FUNDING: Leibniz ScienceCampus Evolutionary Medicine of the Lung, Deutsche Forschungsgemeinschaft, Research Training Group 2501 TransEvo, Rhodes Trust, Stanford University Medical Scientist Training Program, National Institute for Health and Care Research Oxford Biomedical Research Centre, Oxford University Hospitals NHS Foundation Trust, Bill & Melinda Gates Foundation, Wellcome Trust, and Marie Sklodowska-Curie Actions.

High fluoroquinolone resistance proportions among multidrug-resistant tuberculosis driven by dominant L2 Mycobacterium tuberculosis clones in the Mumbai Metropolitan Region.

BACKGROUND: Multidrug-resistant (MDR) Mycobacterium tuberculosis complex (MTBC) strains are a serious health problem in India, also contributing to one-fourth of the global MDR tuberculosis (TB) burden. About 36% of the MDR MTBC strains are reported fluoroquinolone (FQ) resistant leading to high pre-extensively drug-resistant (pre-XDR) and XDR-TB (further resistance against bedaquiline and/or linezolid) rates. Still, factors driving the MDR/pre-XDR epidemic in India are not well defined. METHODS: In a retrospective study, we analyzed 1852 consecutive MTBC strains obtained from patients from a tertiary care hospital laboratory in Mumbai by whole genome sequencing (WGS). Univariate and multivariate statistics was used to investigate factors associated with pre-XDR. Core genome multi locus sequence typing, time scaled haplotypic density (THD) method and homoplasy analysis were used to analyze epidemiological success, and positive selection in different strain groups, respectively. RESULTS: In total, 1016 MTBC strains were MDR, out of which 703 (69.2%) were pre-XDR and 45 (4.4%) were XDR. Cluster rates were high among MDR (57.8%) and pre-XDR/XDR (79%) strains with three dominant L2 (Beijing) strain clusters (Cl 1-3) representing half of the pre-XDR and 40% of the XDR-TB cases. L2 strains were associated with pre-XDR/XDR-TB (P < 0.001) and, particularly Cl 1-3 strains, had high first-line and FQ resistance rates (81.6-90.6%). Epidemic success analysis using THD showed that L2 strains outperformed L1, L3, and L4 strains in short- and long-term time scales. More importantly, L2 MDR and MDR + strains had higher THD success indices than their not-MDR counterparts. Overall, compensatory mutation rates were highest in L2 strains and positive selection was detected in genes of L2 strains associated with drug tolerance (prpB and ppsA) and virulence (Rv2828c). Compensatory mutations in L2 strains were associated with a threefold increase of THD indices, suggesting improved transmissibility. CONCLUSIONS: Our data indicate a drastic increase of FQ resistance, as well as emerging bedaquiline resistance which endangers the success of newly endorsed MDR-TB treatment regimens. Rapid changes in treatment and control strategies are required to contain transmission of highly successful pre-XDR L2 strains in the Mumbai Metropolitan region but presumably also India-wide.

Genome Sequences of Five SARS-CoV-2 Variants from Mumbai, India, Obtained by Nanopore Sequencing.

We report here the genome sequences of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) variants from five coronavirus disease 2019 (COVID-19) patients in Mumbai, India. Viral genomic RNA was isolated from nasopharyngeal swabs and/or respiratory particles from the masks of the patients. Genomic variant analysis determined 8 to 22 mutations, and the variants belong to lineages previously associated with Indian variants.

Crystallization and preliminary X-ray crystallographic analysis of PhoK, an extracellular alkaline phosphatase from Sphingomonas sp. BSAR-1.

Alkaline phosphatases (APs) are widely distributed from microbes to humans and are involved in several important biological processes such as phosphate nutrition, signal transduction and pathogenesis. Alkaline phosphatases are also useful in various industrial applications and in recombinant DNA technology. A new AP enzyme from Sphingomonas sp. strain BSAR-1, termed PhoK, has been shown to be useful in uranium bioprecipitation. PhoK was expressed, purified and crystallized. The crystals belonged to space group P4(3)2(1)2 or P4(1)2(1)2, with unit-cell parameters a = b = 87.37, c = 168.16 A, and contained one enzyme molecule in the asymmetric unit. Native diffraction data have been collected to 1.95 A resolution at the ESRF.

Cloning and overexpression of alkaline phosphatase PhoK from Sphingomonas sp. strain BSAR-1 for bioprecipitation of uranium from alkaline solutions.

Cells of Sphingomonas sp. strain BSAR-1 constitutively expressed an alkaline phosphatase, which was also secreted in the extracellular medium. A null mutant lacking this alkaline phosphatase activity was isolated by Tn5 random mutagenesis. The corresponding gene, designated phoK, was cloned and overexpressed in Escherichia coli strain BL21(DE3). The resultant E. coli strain EK4 overexpressed cellular activity 55 times higher and secreted extracellular PhoK activity 13 times higher than did BSAR-1. The recombinant strain very rapidly precipitated >90% of input uranium in less than 2 h from alkaline solutions (pH, 9 +/- 0.2) containing 0.5 to 5 mM of uranyl carbonate, compared to BSAR-1, which precipitated uranium in >7 h. In both strains BSAR-1 and EK4, precipitated uranium remained cell bound. The EK4 cells exhibited a much higher loading capacity of 3.8 g U/g dry weight in <2 h compared to only 1.5 g U/g dry weight in >7 h in BSAR-1. The data demonstrate the potential utility of genetically engineering PhoK for the bioprecipitation of uranium from alkaline solutions.

Antibiotic Resistance Characterization of Environmental E. coli Isolated from River Mula-Mutha, Pune District, India.

In the current study, ceftazidime- and ciprofloxacin-resistant—or dual drug-resistant (DDR)—E. coli were isolated from river Mula-Mutha, which flows through rural Pune district and Pune city. The DDR E. coli were further examined for antibiotic resistance to six additional antibiotics. The study also included detection of genes responsible for ceftazidime and ciprofloxacin resistance and vectors for horizontal gene transfer. Twenty-eight percent of the identified DDR E. coli were resistant to more than six antibiotics, with 12% being resistant to all eight antibiotics tested. Quinolone resistance was determined through the detection of qnrA, qnrB, qnrS and oqxA genes, whereas cephalosporin resistance was confirmed through detection of TEM, CTX-M-15, CTX-M-27 and SHV genes. Out of 219 DDR E. coli, 8.2% were qnrS positive and 0.4% were qnrB positive. Percentage of isolates positive for the TEM, CTX-M-15 and CTX-M-27 genes were 32%, 46% and 0.9%, respectively. None of the DDR E. coli tested carried the qnrA, SHV and oqxA genes. Percentage of DDR E. coli carrying Class 1 and 2 integrons (mobile genetic elements) were 47% and 8%, respectively. The results showed that antibiotic resistance genes (ARGs) and integrons were present in the E. coli isolated from the river at points adjoining and downstream of Pune city.

Cost-effectiveness of GeneXpert and LED-FM for diagnosis of pulmonary tuberculosis: A systematic review.

BACKGROUND: Early and accurate diagnosis of tuberculosis is a priority for TB programs globally to initiate treatment early and improve treatment outcomes. Currently, Ziehl-Neelsen (ZN) stain-based microscopy, GeneXpert and Light Emitting Diode-Fluorescence Microscopy (LED-FM) are used for diagnosing pulmonary drug sensitive tuberculosis. Published evidence synthesising the cost-effectiveness of these diagnostic tools is scarce. METHODOLOGY: PubMed, EMBASE and Cost-effectiveness analysis registry were searched for studies that reported on the cost-effectiveness of GeneXpert and LED-FM, compared to ZN microscopy for diagnosing pulmonary TB. Risk of bias was assessed independently by four authors using the Consensus Health Economic Criteria (CHEC) extended checklist. The data variables included the study settings, population, type of intervention, type of comparator, year of study, duration of study, type of study design, costs for the test and the comparator and effectiveness indicators. Incremental cost-effectiveness ratio (ICER) was used for assessing the relative cost-effectiveness in this review. RESULTS: Of the 496 studies identified by the search, thirteen studies were included after removing duplicates and studies that did not fulfil inclusion criteria. Four studies compared LED-FM with ZN and nine studies compared GeneXpert with ZN. Three studies used patient cohorts and eight were modelling studies with hypothetical cohorts used to evaluate cost-effectiveness. All these studies were conducted from a health system perspective, with four studies utilising cost utility analysis. There were considerable variations in costing parameters and effectiveness indicators that precluded meta-analysis. The key findings from the included studies suggest that LED-FM and GeneXpert may be cost effective for pulmonary TB diagnosis from a health system perspective. CONCLUSION: Our review identifies a consistent trend of the cost effectiveness of LED-FM and GeneXpert for pulmonary TB diagnosis in different countries with diverse context of socio-economic condition, HIV burden and geographical distribution. However, all the studies used different parameters to estimate the impact of these tools and this underscores the need for improving the methodological issues related to the conduct and reporting of cost-effectiveness studies.

Whole genome sequencing of clinical strains of Mycobacterium tuberculosis from Mumbai, India: A potential tool for determining drug-resistance and strain lineage.

Amplification of drug resistance in Mycobacterium tuberculosis (M.tb) and its transmission are significant barriers in controlling tuberculosis (TB) globally. Diagnostic inaccuracies and delays impede appropriate drug administration, which exacerbates primary and secondary drug resistance. Increasing affordability of whole genome sequencing (WGS) and exhaustive cataloguing of drug resistance mutations is poised to revolutionise TB diagnostics and facilitate personalized drug therapy. However, application of WGS for diagnostics in high endemic areas is yet to be demonstrated. We report WGS of 74 clinical TB isolates from Mumbai, India, characterising genotypic drug resistance to first- and second-line anti-TB drugs. A concordance analysis between phenotypic and genotypic drug susceptibility of a subset of 29 isolates and the sensitivity of resistance prediction to the 4 drugs was calculated, viz. isoniazid-100%, rifampicin-100%, ethambutol-100% and streptomycin-85%. The whole genome based phylogeny showed almost equal proportion of East Asian (27/74) and Central Asian (25/74) strains. Interestingly we also found a clonal group of 9 isolates, of which 7 patients were found to be from the same geographical location and accessed the same health post. This provides the first evidence of epidemiological linkage for tracking TB transmission in India, an approach which has the potential to significantly improve chances of End-TB goals. Finally, the use of Mykrobe Predictor, as a standalone drug resistance and strain typing tool, requiring just few minutes to analyse raw WGS data into tabulated results, implies the rapid clinical applicability of WGS based TB diagnosis.

Synthetic tunable amplifying buffer circuit in E. coli.

While predictable design of a genetic circuit's output is a major goal of synthetic biology, it remains a significant challenge because DNA binding sites in the cell affect the concentration of available transcription factors (TF). To mitigate this problem, we propose to use a TF that results from the (reversible) phosphorylation of protein substrate as a circuit's output. We demonstrate that by comparatively increasing the amounts of substrate and phosphatase, the TF concentration becomes robust to the presence of DNA binding sites and can be kept at a desired value. The circuit's input/output gain can, in turn, be tuned by changing the relative amounts of the substrate and phosphatase, realizing an amplifying buffer circuit with tunable gain. In our experiments in E. coli, we employ phospho-NRI as the output TF, phosphorylated by the NRII kinase, and dephosphorylated by the NRII phosphatase. Amplifying buffer circuits such as ours could be used to insulate a circuit's output from the context, bringing synthetic biology one step closer to modular design.

Retroactivity controls the temporal dynamics of gene transcription.

Just like in many engineering systems, impedance-like effects, called retroactivity, arise at the interconnection of biomolecular circuits, leading to unexpected changes in a circuit's behavior. In this paper, we provide a combined experimental and theoretical study to characterize the effects of retroactivity on the temporal dynamics of a gene transcription module in vivo. The response of the module to an inducer was measured both in isolation and when the module was connected to downstream clients. The connected module, when compared to the isolated module, responded selectively to the introduction of the inducer versus its withdrawal. Specifically, a "sign-sensitive delay" appeared, in which the connected module displayed a time delay in the response to induction and anticipation in the response to de-induction. The extent of these effects can be made larger by increasing the amounts of downstream clients and/or their binding affinity to the output protein of the module. Our experimental results and mathematical formulas make it possible to predict the extent of the change in the dynamic behavior of a module after interconnection. They can be employed to both recover the predictive power of a modular approach to understand systems or as an additional design tool to shape the temporal behavior of gene transcription.

X-ray structure reveals a new class and provides insight into evolution of alkaline phosphatases.

The alkaline phosphatase (AP) is a bi-metalloenzyme of potential applications in biotechnology and bioremediation, in which phosphate monoesters are nonspecifically hydrolysed under alkaline conditions to yield inorganic phosphate. The hydrolysis occurs through an enzyme intermediate in which the catalytic residue is phosphorylated. The reaction, which also requires a third metal ion, is proposed to proceed through a mechanism of in-line displacement involving a trigonal bipyramidal transition state. Stabilizing the transition state by bidentate hydrogen bonding has been suggested to be the reason for conservation of an arginine residue in the active site. We report here the first crystal structure of alkaline phosphatase purified from the bacterium Sphingomonas. sp. Strain BSAR-1 (SPAP). The crystal structure reveals many differences from other APs: 1) the catalytic residue is a threonine instead of serine, 2) there is no third metal ion binding pocket, and 3) the arginine residue forming bidentate hydrogen bonding is deleted in SPAP. A lysine and an aspargine residue, recruited together for the first time into the active site, bind the substrate phosphoryl group in a manner not observed before in any other AP. These and other structural features suggest that SPAP represents a new class of APs. Because of its direct contact with the substrate phosphoryl group, the lysine residue is proposed to play a significant role in catalysis. The structure is consistent with a mechanism of in-line displacement via a trigonal bipyramidal transition state. The structure provides important insights into evolutionary relationships between members of AP superfamily.