Comprehensive identification of pathogenic microbes and antimicrobial resistance genes in food products using nanopore sequencing-based metagenomics.

Foodborne pathogens, particularly antimicrobial-resistant (AMR) bacteria, remain a significant threat to global health. Given the limitations of conventional culture-based approaches, which are limited in scope and time-consuming, metagenomic sequencing of food products emerges as a promising solution. This method provides a fast and comprehensive way to detect the presence of pathogenic microbes and antimicrobial resistance genes (ARGs). Notably, nanopore long-read sequencing provides more accurate bacterial taxonomic classification in comparison to short-read sequencing. Here, we revealed the impact of food types and attributes (origin, retail place, and food processing methods) on microbial communities and the AMR profile using nanopore metagenomic sequencing. We analyzed a total of 260 food products, including raw meat, sashimi, and ready-to-eat (RTE) vegetables. Clostridium botulinum, Acinetobacter baumannii, and Vibrio parahaemolyticus were identified as the top three foodborne pathogens in raw meat and sashimi. Importantly, even with low pathogen abundance, higher percentages of samples containing carbapenem and cephalosporin resistance genes were identified in chicken and RTE vegetables, respectively. In parallel, our results demonstrated that fresh, peeled, and minced foods exhibited higher levels of pathogenic bacteria. In conclusion, this comprehensive study offers invaluable data that can contribute to food safety assessments and serve as a basis for quality indicators.

Long-Read Sequencing with Hierarchical Clustering for Antiretroviral Resistance Profiling of Mixed Human Immunodeficiency Virus Quasispecies.

BACKGROUND: HIV infections often develop drug resistance mutations (DRMs), which can increase the risk of virological failure. However, it has been difficult to determine if minor mutations occur in the same genome or in different virions using Sanger sequencing and short-read sequencing methods. Oxford Nanopore Technologies (ONT) sequencing may improve antiretroviral resistance profiling by allowing for long-read clustering. METHODS: A new ONT sequencing-based method for profiling DRMs in HIV quasispecies was developed and validated. The method used hierarchical clustering of long amplicons that cover regions associated with different types of antiretroviral drugs. A gradient series of an HIV plasmid and 2 plasma samples was prepared to validate the clustering performance. The ONT results were compared to those obtained with Sanger sequencing and Illumina sequencing in 77 HIV-positive plasma samples to evaluate the diagnostic performance. RESULTS: In the validation study, the abundance of detected quasispecies was concordant with the predicted result with the R2 of > 0.99. During the diagnostic evaluation, 59/77 samples were successfully sequenced for DRMs. Among 18 failed samples, 17 were below the limit of detection of 303.9 copies/µL. Based on the receiver operating characteristic analysis, the ONT workflow achieved an F1 score of 0.96 with a cutoff of 0.4 variant allele frequency. Four cases were found to have quasispecies with DRMs, in which 2 harbored quasispecies with more than one class of DRMs. Treatment modifications were recommended for these cases. CONCLUSIONS: Long-read sequencing coupled with hierarchical clustering could differentiate the quasispecies resistance profiles in HIV-infected samples, providing a clearer picture for medical care.

Probable Animal-to-Human Transmission of Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) Delta Variant AY.127 Causing a Pet Shop-Related Coronavirus Disease 2019 (COVID-19) Outbreak in Hong Kong.

BACKGROUND: Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) can infect human and other mammals, including hamsters. Syrian (Mesocricetus auratus) and dwarf (Phodopus sp.) hamsters are susceptible to SARS-CoV-2 infection in the laboratory setting. However, pet shop-related Coronavirus Disease 2019 (COVID-19) outbreaks have not been reported. METHODS: We conducted an investigation of a pet shop-related COVID-19 outbreak due to Delta variant AY.127 involving at least 3 patients in Hong Kong. We tested samples collected from the patients, environment, and hamsters linked to this outbreak and performed whole genome sequencing analysis of the reverse transcription polymerase chain reaction (RT-PCR)-positive samples. RESULTS: The patients included a pet shop keeper (Patient 1), a female customer of the pet shop (Patient 2), and the husband of Patient 2 (Patient 3). Investigation showed that 17.2% (5/29) and 25.5% (13/51) environmental specimens collected from the pet shop and its related warehouse, respectively, tested positive for SARS-CoV-2 RNA by RT-PCR. Among euthanized hamsters randomly collected from the storehouse, 3% (3/100) tested positive for SARS-CoV-2 RNA by RT-PCR and seropositive for anti-SARS-CoV-2 antibody by enzyme immunoassay. Whole genome analysis showed that although all genomes from the outbreak belonged to the Delta variant AY.127, there were at least 3 nucleotide differences among the genomes from different patients and the hamster cages. Genomic analysis suggests that multiple strains have emerged within the hamster population, and these different strains have likely transmitted to human either via direct contact or via the environment. CONCLUSIONS: Our study demonstrated probable hamster-to-human transmission of SARS-CoV-2. As pet trading is common around the world, this can represent a route of international spread of this pandemic virus.

The Clinical Utility of Two High-Throughput 16S rRNA Gene Sequencing Workflows for Taxonomic Assignment of Unidentifiable Bacterial Pathogens in Matrix-Assisted Laser Desorption Ionization-Time of Flight Mass Spectrometry.

Bacterial pathogens that cannot be identified using matrix-assisted laser desorption ionization-time of flight mass spectrometry (MALDI-TOF MS) are occasionally encountered in clinical laboratories. The 16S rRNA gene is often used for sequence-based analysis to identify these bacterial species. Nevertheless, traditional Sanger sequencing is laborious, time-consuming, and low throughput. Here, we compared two commercially available 16S rRNA gene sequencing tests that are based on Illumina and Nanopore sequencing technologies, respectively, in their ability to identify the species of 172 clinical isolates that failed to be identified by MALDI-TOF MS. Sequencing data were analyzed by the respective built-in programs (MiSeq Reporter software of Illumina and Epi2me of Nanopore) and BLAST+ (v2.11.0). Their agreement with Sanger sequencing on species-level identification was determined. Discrepancies were resolved by whole-genome sequencing. The diagnostic accuracy of each workflow was determined using the composite sequencing result as the reference standard. Despite the high base-calling accuracy of Illumina sequencing, we demonstrated that the Nanopore workflow had a higher taxonomic resolution at the species level. Using built-in analysis algorithms, the concordance of Sanger 16S with the Illumina and Nanopore workflows was 33.14% and 87.79%, respectively. The agreement was 65.70% and 83.14%, respectively, when BLAST+ was used for analysis. Compared with the reference standard, the diagnostic accuracy of Nanopore 16S was 96.36%, which was identical to that of Sanger 16S and better than that of Illumina 16S (69.07%). The turnaround time of the Illumina workflow and the Nanopore workflow was 78 h and 8.25 h, respectively. The per-sample cost of the Illumina and Nanopore workflows was US$28.5 and US$17.7, respectively.

Explosive outbreak of SARS-CoV-2 Omicron variant is associated with vertical transmission in high-rise residential buildings in Hong Kong.

The phenomenon of vertical transmission of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) in high-rise residential buildings (HRRBs) is unique in our densely populated cosmopolitan city. The compulsory testing of a whole building under the scheme of restriction-testing declaration (RTD) during the fourth wave (non-Omicron variant) and fifth wave (mostly Omicron variant) of COVID-19 outbreak in Hong Kong allowed us to study the prevalence of this phenomenon, which may represent a form of airborne transmission. From 23 January 2021 to 24 March 2022, 25,450 (5.8%) of 436,397 residents from 223 (63.0%) of 354 HRRBs under RTD were test-positive for SARS-CoV-2. Using the clustering of cases among vertically aligned flats with shared drainage stack and lightwell as a surrogate marker of vertical transmission, the number of vertically aligned flats with positive COVID-19 cases was significantly higher in the fifth wave compared with the fourth wave (14.2%, 6471/45,531 vs 0.24%, 3/1272; p < 0.001; or 2212 vs 1 per-million-flats; p < 0.001). Excluding 22,801 residents from 38 HRRBs who were tested negative outside the 12-week periods selected in fourth and fifth waves, the positive rate among residents was significantly higher among residents during the fifth wave than the fourth wave (6.5%, 25,434/389,700 vs 0.07%, 16/23,896; p < 0.001). Within the flats with COVID-19 cases, the proportion of vertically aligned flats was also significantly higher in the fifth wave than in the fourth wave (95.6%, 6471/6766 vs 30.0%, 3/10, p < 0.001). The proportion of HRRBs with COVID-19 cases was significantly higher during the corresponding 12-week period chosen for comparison (78.2%, 219/280 vs 11.1%, 4/36; p < 0.001). Whole-genome phylogenetic analysis of 332 viral genomes showed that Omicron BA.2 was the predominant strain, supporting the high transmissibility of BA.2 by airborne excreta-aerosol route in HRRBs of Hong Kong.

Nosocomial Outbreak of Coronavirus Disease 2019 by Possible Airborne Transmission Leading to a Superspreading Event.

BACKGROUND: Nosocomial outbreaks with superspreading of coronavirus disease 2019 due to a possible airborne transmission have not been reported. METHODS: Epidemiological analysis, environmental samplings, and whole-genome sequencing (WGS) were performed for a hospital outbreak. RESULTS: A superspreading event that involved 12 patients and 9 healthcare workers (HCWs) occurred within 9 days in 3 of 6 cubicles at an old-fashioned general ward with no air exhaust built within the cubicles. The environmental contamination by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) RNA was significantly higher in air grilles (>2 m from patients' heads and not within reach) than on high-touch clinical surfaces (36.4%, 8 of 22 vs 3.4%, 1 of 29, P = .003). Six (66.7%) of 9 contaminated air exhaust grilles were located outside patient cubicles. The clinical attack rate of patients was significantly higher than of HCWs (15.4%, 12 of 78 exposed patients vs 4.6%, 9 of 195 exposed HCWs, P = .005). Moreover, the clinical attack rate of ward-based HCWs was significantly higher than of nonward-based HCWs (8.1%, 7 of 68 vs 1.8%, 2 of 109, P = .045). The episodes (mean ±â€…standard deviation) of patient-care duty assignment in the cubicles was significantly higher among infected ward-based HCWs than among noninfected ward-based HCWs (6.0 ±â€…2.4 vs 3.0 ±â€…2.9, P = .012) during the outbreak period. The outbreak strains belong to SARS-CoV-2 lineage B.1.36.27 (GISAID clade GH) with the unique S-T470N mutation on WGS. CONCLUSIONS: This nosocomial point source superspreading event due to possible airborne transmission demonstrates the need for stringent SARS-CoV-2 screening at admission to healthcare facilities and better architectural design of ventilation systems to prevent such outbreaks. Portable high-efficiency particulate filters were installed in each cubicle to improve ventilation before resumption of clinical service.

Territorywide Study of Early Coronavirus Disease Outbreak, Hong Kong, China.

Initial cases of coronavirus disease in Hong Kong were imported from mainland China. A dramatic increase in case numbers was seen in February 2020. Most case-patients had no recent travel history, suggesting the presence of transmission chains in the local community. We collected demographic, clinical, and epidemiologic data from 50 patients, who accounted for 53.8% of total reported case-patients as of February 28, 2020. We performed whole-genome sequencing to determine phylogenetic relationship and transmission dynamics of severe acute respiratory syndrome coronavirus 2 infections. By using phylogenetic analysis, we attributed the community outbreak to 2 lineages; 1 harbored a common mutation, Orf3a-G251V, and accounted for 88.0% of the cases in our study. The estimated time to the most recent common ancestor of local coronavirus disease outbreak was December 24, 2019, with an evolutionary rate of 3.04 × 10-3 substitutions/site/year. The reproduction number was 1.84, indicating ongoing community spread.

Nanopore Sequencing Reveals Novel Targets for Detection and Surveillance of Human and Avian Influenza A Viruses.

Accurate detection of influenza A virus (IAV) is crucial for patient management, infection control, and epidemiological surveillance. The World Health Organization and the Centers for Disease Control and Prevention have recommended using the M gene as the diagnostic gene target for reverse-transcription-PCR (RT-PCR). However, M gene RT-PCR has reduced sensitivity for recent IAV due to novel gene mutations. Here, we sought to identify novel diagnostic targets for the molecular detection of IAV using long-read third-generation sequencing. Direct nanopore sequencing from 18 nasopharyngeal specimens and one saliva specimen showed that the 5' and 3' ends of the PB2 gene and the entire NS gene were highly abundant. Primers selected for PB2 and NS genes were well matched with seasonal or avian IAV gene sequences. Our novel PB2 and NS gene real-time RT-PCR assays showed limits of detection similar to or lower than that of M gene RT-PCR and achieved 100% sensitivity and specificity in the detection of A(H1N1), A(H3N2), and A(H7N9) in nasopharyngeal and saliva specimens. For 10 patients with IAV detected by M gene RT-PCR conversion in sequentially collected specimens, NS and/or PB2 gene RT-PCR was positive in 2 (20%) of the initial specimens that were missed by M gene RT-PCR. In conclusion, we have shown that PB2 or NS gene RT-PCRs are suitable alternatives to the recommended M gene RT-PCR for diagnosis of IAV. Long-read nanopore sequencing facilitates the identification of novel diagnostic targets.

Targeted-Sequencing Workflows for Comprehensive Drug Resistance Profiling of Mycobacterium tuberculosis Cultures Using Two Commercial Sequencing Platforms: Comparison of Analytical and Diagnostic Performance, Turnaround Time, and Cost.

BACKGROUND: The emergence of Mycobacterium tuberculosis with complex drug resistance profiles necessitates a rapid and comprehensive drug susceptibility test for guidance of patient treatment. We developed two targeted-sequencing workflows based on Illumina MiSeq and Nanopore MinION for the prediction of drug resistance in M. tuberculosis toward 12 antibiotics. METHODS: A total of 163 M. tuberculosis isolates collected from Hong Kong and Ethiopia were subjected to a multiplex PCR for simultaneous amplification of 19 drug resistance-associated genetic regions. The amplicons were then barcoded and sequenced in parallel on MiSeq and MinION in respective batch sizes of 24 and 12 samples. A web-based bioinformatics pipeline, BacterioChek-TB, was developed to translate the raw datasets into clinician-friendly reports. RESULTS: Both platforms successfully sequenced all samples with mean read depths of 1,127× and 1,649×, respectively. The variant calling by MiSeq and MinION could achieve 100% agreement if variants with an allele frequency of <40% reported by MinION were excluded. Both workflows achieved a mean clinical sensitivity of 94.8% and clinical specificity of 98.0% when compared with phenotypic drug susceptibility test (pDST). Turnaround times for the MiSeq and MinION workflows were 38 and 15 h, facilitating the delivery of treatment guidance at least 17-18 days earlier than pDST, respectively. The higher cost per sample on the MinION platform ($71.56) versus the MiSeq platform ($67.83) was attributed to differences in batching capabilities. CONCLUSION: Our study demonstrates the interchangeability of MiSeq and MinION platforms for generation of accurate and actionable results for the treatment of tuberculosis.

Direct Detection of Pyrazinamide Resistance in Mycobacterium tuberculosis by Use of pncA PCR Sequencing.

An in-house-developed pncA sequencing assay for analysis of pyrazinamide (PZA) resistance was evaluated using 162 archived Mycobacterium tuberculosis complex (MTBC) isolates with phenotypic PZA susceptibility profiles that were well defined by analysis of Bactec MGIT 960 PZA kit and PZase activity data. Preliminary results showed 100% concordance between pncA sequencing and phenotypic PZA drug susceptibility test (DST) results among archived isolates. Also, 637 respiratory specimens were prospectively collected, and 158 were reported as MTBC positive by the Abbott Realtime MTB assay (96.3% sensitivity [95% confidence interval {CI}: 92.2% to 98.7%]; 100% specificity [95% CI: 99.2% to 100.0%]). Genotypic and phenotypic PZA resistance profiles of these 158 MTBC-positive specimens were analyzed by pncA sequencing and Bactec MGIT 960 PZA kit, respectively. For analysis of PZA resistance, pncA sequencing detected pncA mutations in 5/5 (100%) phenotypic PZA-resistant respiratory specimens within 4 working days. No pncA mutations were detected among PZA-susceptible specimens. Combining archived isolates with prospective specimens, 27 were identified as phenotypic PZA resistant with pncA mutation. Among these 27 samples, 6/27 (22.2%) phenotypic PZA-resistant strains carried novel pncA mutations without rpsA and panD mutations. These included 5 with mutations (a deletion [Del] at 383T [Del383T], Del 380 to 390, insertion of A [A Ins] at position 127, A Ins at position 407, and G Ins at position 508) in pncA structural genes and 1 with a mutation (T-12C) at the pncA promoter region. All six of these strains had no or reduced PZase activities, indicating that the novel mutations might confer PZA resistance. Additionally, 25/27 phenotypic PZA-resistant strains were confirmed multidrug-resistant tuberculosis (MDR-TB) strains. As PZA is commonly used in MDR-TB treatment regimens, direct pncA sequencing will rapidly detect PZA resistance and facilitate judicious use of PZA in treating PZA-susceptible MDR-TB.

GAAP: Genome-organization-framework-Assisted Assembly Pipeline for prokaryotic genomes.

BACKGROUND: Next-generation sequencing (NGS) technologies have greatly promoted the genomic study of prokaryotes. However, highly fragmented assemblies due to short reads from NGS are still a limiting factor in gaining insights into the genome biology. Reference-assisted tools are promising in genome assembly, but tend to result in false assembly when the assigned reference has extensive rearrangements. RESULTS: Herein, we present GAAP, a genome assembly pipeline for scaffolding based on core-gene-defined Genome Organizational Framework (cGOF) described in our previous study. Instead of assigning references, we use the multiple-reference-derived cGOFs as indexes to assist in order and orientation of the scaffolds and build a skeleton structure, and then use read pairs to extend scaffolds, called local scaffolding, and distinguish between true and chimeric adjacencies in the scaffolds. In our performance tests using both empirical and simulated data of 15 genomes in six species with diverse genome size, complexity, and all three categories of cGOFs, GAAP outcompetes or achieves comparable results when compared to three other reference-assisted programs, AlignGraph, Ragout and MeDuSa. CONCLUSIONS: GAAP uses both cGOF and pair-end reads to create assemblies in genomic scale, and performs better than the currently available reference-assisted assembly tools as it recovers more assemblies and makes fewer false locations, especially for species with extensive rearranged genomes. Our method is a promising solution for reconstruction of genome sequence from short reads of NGS.

An upstream truncation of the furA-katG operon confers high-level isoniazid resistance in a Mycobacterium tuberculosis clinical isolate with no known resistance-associated mutations.

Although the major causes of isoniazid (INH) resistance in Mycobacterium tuberculosis are confined to structural mutations in katG and promoter mutations in the mabA-inhA operon, a significant proportion of INH-resistant strains have unknown resistance mechanisms. Recently, we identified a high-level INH-resistant M. tuberculosis clinical isolate, GB005, with no known resistance-associated mutations. A comprehensive study was performed to investigate the molecular basis of drug resistance in this strain. Although no mutations were found throughout the katG and furA-katG intergenic region, the katG expression and the catalase activity were greatly diminished compared to those in H37Rv (P < 0.01). Northern blotting revealed that the katG transcript from the isolate was smaller than that of H37Rv. Sequencing analysis of furA and upstream genes discovered a 7.2-kb truncation extended from the 96th base preceding the initiation codon of katG. Complementation of the M. tuberculosis Δ(furA-katG) strain with katG and different portions of the truncated region identified a 134-bp upstream fragment of furA that was essential for full catalase activity and INH susceptibility in M. tuberculosis. The promoter activity of this fragment was also shown to be stronger than that of the furA-katG intergenic region (P < 0.01). Collectively, these findings demonstrate that deletion of the 134-bp furA upstream fragment is responsible for the reduction in katG expression, resulting in INH resistance in GB005. To our knowledge, this is the first report showing that deletion of the upstream region preceding the furA-katG operon causes high-level INH resistance in a clinical isolate of M. tuberculosis.

Community Engagement in Vaccination Promotion: Systematic Review and Meta-Analysis.

BACKGROUND: Community engagement plays a vital role in global immunization strategies, offering the potential to overcome vaccination hesitancy and enhance vaccination confidence. Although there is significant backing for community engagement in health promotion, the evidence supporting its effectiveness in vaccination promotion is fragmented and of uncertain quality. OBJECTIVE: This review aims to systematically examine the effectiveness of different contents and extent of community engagement for promoting vaccination rates. METHODS: This study was performed in accordance with the PRISMA (Preferred Reporting Items for Systematic Reviews and Meta-Analyses) guidelines. A comprehensive and exhaustive literature search was performed in 4 English databases (PubMed, Embase, Web of Science, and Cochrane Library) and 2 Chinese databases (CNKI and Wan Fang) to identify all possible articles. Original research articles applying an experimental study design that investigated the effectiveness of community engagement in vaccination promotion were eligible for inclusion. Two reviewers independently performed the literature search, study selection, quality assessment, and data extraction. Discrepancies were resolved through discussion, with the arbitration of a third reviewer where necessary. RESULTS: A total of 20 articles out of 11,404 records from 2006 to 2021 were retrieved. The studies used various designs: 12 applied single-group pre-post study designs, 5 were cluster randomized controlled trials (RCTs), and 3 were non-RCTs. These studies targeted multiple vaccines, with 8 focusing on children's immunization, 8 on human papillomavirus vaccine, 3 on hepatitis B virus vaccine, and 1 on COVID-19 vaccine. The meta-analysis revealed significant increases in vaccination rates both in pre-post comparison (rate difference [RD] 0.34, 95% CI 0.21-0.47, I2=99.9%, P<.001) and between-group comparison (RD 0.18, 95% CI 0.07-0.29, I2=98.4%, P<.001). The meta-analysis revealed that participant recruitment had the largest effect size (RD 0.51, 95% CI 0.36-0.67, I2=99.9%, P<.001), followed by intervention development (RD 0.36, 95% CI 0.23-0.50, I2=100.0%, P<.001), intervention implementation (RD 0.35, 95% CI 0.22-0.47, I2=99.8%, P<.001), and data collection (RD 0.34, 95% CI 0.19-0.50, I2=99.8%, P<.001). The meta-analysis indicated that high community engagement extent yielded the largest effect size (RD 0.49, 95% CI 0.17-0.82, I2=100.0%, P<.001), followed by moderate community engagement extent (RD 0.45, 95% CI 0.33-0.58, I2=99.6%, P<.001) and low community engagement extent (RD 0.15, 95% CI 0.05-0.25, I2=99.2%, P<.001). The meta-analysis revealed that "health service support" demonstrated the largest effect sizes (RD 0.45, 95% CI 0.25-0.65, I2=99.9%, P<.001), followed by "health education and discussion" (RD 0.39, 95% CI 0.20-0.58, I2=99.7%, P<.001), "follow-up and reminder" (RD 0.33, 95% CI 0.23-0.42, I2=99.3%, P<.001), and "social marketing campaigns and community mobilization" (RD 0.24, 95% CI 0.06-0.41, I2=99.9%, P<.001). CONCLUSIONS: The results of this meta-analysis supported the effectiveness of community engagement in vaccination promotion with variations in terms of engagement contents and extent. Community engagement required a "fit-for-purpose" approach rather than a "one-size-fits-all" approach to maximize the effectiveness of vaccine promotion. TRIAL REGISTRATION: PROSPERO CRD42022339081; https://www.crd.york.ac.uk/prospero/display_record.php?RecordID=339081.

An AI-empowered indoor digital contact tracing system for COVID-19 outbreaks in residential care homes.

An AI-empowered indoor digital contact-tracing system was developed using a centralized architecture and advanced low-energy Bluetooth technologies for indoor positioning, with careful preservation of privacy and data security. We analyzed the contact pattern data from two RCHs and investigated a COVID-19 outbreak in one study site. To evaluate the effectiveness of the system in containing outbreaks with minimal contacts under quarantine, a simulation study was conducted to compare the impact of different quarantine strategies on outbreak containment within RCHs. The significant difference in contact hours between weekdays and weekends was observed for some pairs of RCH residents and staff during the two-week data collection period. No significant difference between secondary cases and uninfected contacts was observed in a COVID-19 outbreak in terms of their demographics and contact patterns. Simulation results based on the collected contact data indicated that a threshold of accumulative contact hours one or two days prior to diagnosis of the index case could dramatically increase the efficiency of outbreak containment within RCHs by targeted isolation of the close contacts. This study demonstrated the feasibility and efficiency of employing an AI-empowered system in indoor digital contact tracing of outbreaks in RCHs in the post-pandemic era.

Genetic mechanisms of co-emergence of INH-resistant Mycobacterium tuberculosis strains during the standard course of antituberculosis therapy.

The incidence of isoniazid (INH) resistant Mycobacterium tuberculosis is increasing globally. This study aimed to identify the molecular mechanisms behind the development of INH resistance in M. tuberculosis strains collected from the same patients during the standard course of treatment. Three M. tuberculosis strains were collected from a patient before and during antituberculosis (anti-TB) therapy. The strains were characterized using phenotypic drug susceptibility tests, Mycobacterial Interspersed Repeated Unit-Variable-Number Tandem Repeats (MIRU-VNTR), and whole-genome sequencing (WGS) to identify mutations associated with INH resistance. To validate the role of the novel mutations in INH resistance, the mutated katG genes were electroporated into a KatG-deleted M. tuberculosis strain (GA03). Three-dimensional structures of mutated KatG were modeled to predict their impact on INH binding. The pre-treatment strain was susceptible to INH. However, two INH-resistant strains were isolated from the patient after anti-TB therapy. MIRU-VNTR and WGS revealed that the three strains were clonally identical. A missense mutation (P232L) and a nonsense mutation (Q461Stop) were identified in the katG of the two post-treatment strains, respectively. Transformation experiments showed that katG of the pre-treatment strain restored INH susceptibility in GA03, whereas the mutated katG genes from the post-treatment strains rendered negative catalase activity and INH resistance. The protein model indicated that P232L reduced INH-KatG binding affinity while Q461Stop truncated gene transcription. Our results showed that the two katG mutations, P232L and Q461Stop, accounted for the co-emergence of INH-resistant clones during anti-TB therapy. The inclusion of these mutations in the design of molecular assays could increase the diagnostic performance.IMPORTANCEThe evolution of drug-resistant strains of Mycobacterium tuberculosis within the lung lesions of a patient has a detrimental impact on treatment outcomes. This is particularly concerning for isoniazid (INH), which is the most potent first-line antimycobacterial drug. However, the precise genetic factors responsible for drug resistance in patients have not been fully elucidated, with approximately 15% of INH-resistant strains harboring unknown genetic factors. This raises concerns about the emergence of drug-resistant clones within patients, further contributing to the global epidemic of resistance. In this study, we revealed the presence of two novel katG mutations, which emerged independently due to the stress exerted by antituberculosis (anti-TB) treatment on a parental strain. Importantly, we experimentally demonstrated the functional significance of both mutations in conferring resistance to INH. Overall, this research sheds light on the genetic mechanisms underlying the evolution of INH resistance within patients and provides valuable insights for improving diagnostic performance by targeting specific mutations.

Transmission Patterns of Co-Circulation of Omicron Sub-Lineages in Hong Kong SAR, China, a City with Rigorous Social Distancing Measures, in 2022.

OBJECTIVE: This study aimed to characterize the changing landscape of circulating SARS-CoV-2 lineages in the local community of Hong Kong throughout 2022. We examined how adjustments to quarantine arrangements influenced the transmission pattern of Omicron variants in a city with relatively rigorous social distancing measures at that time. METHODS: In 2022, a total of 4684 local SARS-CoV-2 genomes were sequenced using the Oxford Nanopore GridION sequencer. SARS-CoV-2 consensus genomes were generated by MAFFT, and the maximum likelihood phylogeny of these genomes was determined using IQ-TREE. The dynamic changes in lineages were depicted in a time tree created by Nextstrain. Statistical analysis was conducted to assess the correlation between changes in the number of lineages and adjustments to quarantine arrangements. RESULTS: By the end of 2022, a total of 83 SARS-CoV-2 lineages were identified in the community. The increase in the number of new lineages was significantly associated with the relaxation of quarantine arrangements (One-way ANOVA, F(5, 47) = 18.233, p < 0.001)). Over time, Omicron BA.5 sub-lineages replaced BA.2.2 and became the predominant Omicron variants in Hong Kong. The influx of new lineages reshaped the dynamics of Omicron variants in the community without fluctuating the death rate and hospitalization rate (One-way ANOVA, F(5, 47) = 2.037, p = 0.091). CONCLUSION: This study revealed that even with an extended mandatory quarantine period for incoming travelers, it may not be feasible to completely prevent the introduction and subsequent community spread of highly contagious Omicron variants. Ongoing molecular surveillance of COVID-19 remains essential to monitor the emergence of new recombinant variants.

Identification of microbial community in the urban environment: The concordance between conventional culture and nanopore 16S rRNA sequencing.

Introduction: Microbes in the built environment have been implicated as a source of infectious diseases. Bacterial culture is the standard method for assessing the risk of exposure to pathogens in urban environments, but this method only accounts for <1% of the diversity of bacteria. Recently, full-length 16S rRNA gene analysis using nanopore sequencing has been applied for microbial evaluations, resulting in a rise in the development of long-read taxonomic tools for species-level classification. Regarding their comparative performance, there is, however, a lack of information. Methods: Here, we aim to analyze the concordance of the microbial community in the urban environment inferred by multiple taxonomic classifiers, including ARGpore2, Emu, Kraken2/Bracken and NanoCLUST, using our 16S-nanopore dataset generated by MegaBLAST, as well as assess their abilities to identify culturable species based on the conventional culture results. Results: According to our results, NanoCLUST was preferred for 16S microbial profiling because it had a high concordance of dominant species and a similar microbial profile to MegaBLAST, whereas Kraken2/Bracken, which had similar clustering results as NanoCLUST, was also desirable. Second, for culturable species identification, Emu with the highest accuracy (81.2%) and F1 score (29%) for the detection of culturable species was suggested. Discussion: In addition to generating datasets in complex communities for future benchmarking studies, our comprehensive evaluation of the taxonomic classifiers offers recommendations for ongoing microbial community research, particularly for complex communities using nanopore 16S rRNA sequencing.

Evaluation of Mycobacterium tuberculosis enrichment in metagenomic samples using ONT adaptive sequencing and amplicon sequencing for identification and variant calling.

Sensitive detection of Mycobacterium tuberculosis (TB) in small percentages in metagenomic samples is essential for microbial classification and drug resistance prediction. However, traditional methods, such as bacterial culture and microscopy, are time-consuming and sometimes have limited TB detection sensitivity. Oxford nanopore technologies (ONT) MinION sequencing allows rapid and simple sample preparation for sequencing. Its recently developed adaptive sequencing selects reads from targets while allowing real-time base-calling to achieve sequence enrichment or depletion during sequencing. Another common enrichment method is PCR amplification of the target TB genes. In this study, we compared both methods using ONT MinION sequencing for TB detection and variant calling in metagenomic samples using both simulation runs and those with synthetic and patient samples. We found that both methods effectively enrich TB reads from a high percentage of human (95%) and other microbial DNA. Adaptive sequencing with readfish and UNCALLDE achieved a 3.9-fold and 2.2-fold enrichment compared to the control run. We provide a simple automatic analysis framework to support the detection of TB for clinical use, openly available at https://github.com/HKU-BAL/ONT-TB-NF . Depending on the patient's medical condition and sample type, we recommend users evaluate and optimize their workflow for different clinical specimens to improve the detection limit.

Spatiotemporal Distribution of Tuberculosis in the Oromia Region of Ethiopia: A Hotspot Analysis.

Tuberculosis (TB) is a major public health concern in low- and middle-income countries including Ethiopia. This study aimed to assess the spatiotemporal distribution of TB and identify TB risk factors in Ethiopia's Oromia region. Descriptive and spatiotemporal analyses were conducted. Bayesian spatiotemporal modeling was used to identify covariates that accounted for variability in TB and its spatiotemporal distribution. A total of 206,278 new pulmonary TB cases were reported in the Oromia region between 2018 and 2022, with the lowest annual TB case notification (96.93 per 100,000 population) reported in 2020 (i.e., during the COVID-19 pandemic) and the highest TB case notification (106.19 per 100,000 population) reported in 2019. Substantial spatiotemporal variations in the distribution of notified TB case notifications were observed at zonal and district levels with most of the hotspot areas detected in the northern and southern parts of the region. The spatiotemporal distribution of notified TB incidence was positively associated with different ecological variables including temperature (ß = 0.142; 95% credible interval (CrI): 0.070, 0.215), wind speed (ß = -0.140; 95% CrI: -0.212, -0.068), health service coverage (ß = 0.426; 95% CrI: 0.347, 0.505), and population density (ß = 0.491; 95% CrI: 0.390, 0.594). The findings of this study indicated that preventive measures considering socio-demographic and health system factors can be targeted to high-risk areas for effective control of TB in the Oromia region. Further studies are needed to develop effective strategies for reducing the burden of TB in hotspot areas.

The clinical utility of Nanopore 16S rRNA gene sequencing for direct bacterial identification in normally sterile body fluids.

The prolonged incubation period of traditional culture methods leads to a delay in diagnosing invasive infections. Nanopore 16S rRNA gene sequencing (Nanopore 16S) offers a potential rapid diagnostic approach for directly identifying bacteria in infected body fluids. To evaluate the clinical utility of Nanopore 16S, we conducted a study involving the collection and sequencing of 128 monomicrobial samples, 65 polymicrobial samples, and 20 culture-negative body fluids. To minimize classification bias, taxonomic classification was performed using 3 analysis pipelines: Epi2me, Emu, and NanoCLUST. The result was compared to the culture references. The limit of detection of Nanopore 16S was also determined using simulated bacteremic blood samples. Among the three classifiers, Emu demonstrated the highest concordance with the culture results. It correctly identified the taxon of 125 (97.7%) of the 128 monomicrobial samples, compared to 109 (85.2%) for Epi2me and 102 (79.7%) for NanoCLUST. For the 230 cultured species in the 65 polymicrobial samples, Emu correctly identified 188 (81.7%) cultured species, compared to 174 (75.7%) for Epi2me and 125 (54.3%) for NanoCLUST. Through ROC analysis on the monomicrobial samples, we determined a threshold of relative abundance at 0.058 for distinguishing potential pathogens from background in Nanopore 16S. Applying this threshold resulted in the identification of 107 (83.6%), 117 (91.4%), and 114 (91.2%) correctly detected samples for Epi2me, Emu, and NanoCLUST, respectively, in the monomicrobial samples. Nanopore 16S coupled with Epi2me could provide preliminary results within 6 h. However, the ROC analysis of polymicrobial samples exhibited a random-like performance, making it difficult to establish a threshold. The overall limit of detection for Nanopore 16S was found to be about 90 CFU/ml.