Search

Background: Malaria continues to be a major threat to global public health. Whole genome sequencing (WGS) of the underlying Plasmodium parasites has provided insights into the genomic epidemiology of malaria. Genome sequencing is rapidly gaining traction as a diagnostic and surveillance tool for clinical settings, where the profiling of co-infections, identification of imported malaria parasites, and detection of drug resistance are crucial for infection control and disease elimination. To support this informatically, we have developed the Malaria-Profiler tool, which rapidly (within minutes) predicts Plasmodium species, geographical source, and resistance to antimalarial drugs directly from WGS data. Results: The online and command line versions of Malaria-Profiler detect ~ 250 markers from genome sequences covering Plasmodium speciation, likely geographical source, and resistance to chloroquine, sulfadoxine-pyrimethamine (SP), and other anti-malarial drugs for P. falciparum, but also providing mutations for orthologous resistance genes in other species. The predictive performance of the mutation library was assessed using 9321 clinical isolates with WGS and geographical data, with most being single-species infections (P. falciparum 7152/7462, P. vivax 1502/1661, P. knowlesi 143/151, P. malariae 18/18, P. ovale ssp. 5/5), but co-infections were identified (456/9321; 4.8%). The accuracy of the predicted geographical profiles was high to both continental (96.1%) and regional levels (94.6%). For P. falciparum, markers were identified for resistance to chloroquine (49.2%; regional range: 24.5% to 100%), sulfadoxine (83.3%; 35.4– 90.5%), pyrimethamine (85.4%; 80.0–100%) and combined SP (77.4%). Markers associated with the partial resistance of artemisinin were found in WGS from isolates sourced from Southeast Asia (30.6%). Conclusions: Malaria-Profiler is a user-friendly tool that can rapidly and accurately predict the geographical regional source and anti-malarial drug resistance profiles across large numbers of samples with WGS data. The software is flexible with modifiable bioinformatic pipelines. For example, it is possible to select the sequencing platform, display specific variants, and customise the format of outputs. With the increasing application of next-generation sequencing platforms on Plasmodium DNA, Malaria-Profiler has the potential to be integrated into point-of-care and surveillance settings, thereby assisting malaria control. Malaria-Profiler is available online (bioinformatics.lshtm.ac.uk/malaria-profiler) and as standalone software (https://github.com/jodyphelan/malaria-profiler). [ABSTRACT FROM AUTHOR]

Flaviviruses are a diverse group of RNA viruses, which include the etiological agents of Zika, dengue and yellow fever that are transmitted by mosquitoes. Flaviviruses do not encode reverse transcriptase and cannot reverse transcribe into DNA, yet DNA sequences of flaviviruses are found both integrated in the chromosomes of Aedes aegypti mosquitoes and as extrachromosomal sequences. We have previously examined the Ae. aegypti reference genome to identify flavivirus integrations and analyzed conservation of these sequences among whole-genome data of 464 Ae. aegypti collected across 10 countries globally. Here, we extended this analysis by identifying flavivirus sequences in these samples independently of the Ae. aegypti reference assembly. Our aim was to identify the complete set of viral sequences, including those absent in the reference genome, and their geographical distribution. We compared the identified sequences using BLASTn and applied machine learning methods to identify clusters of similar sequences. Apart from clusters of sequences that correspond to the four viral integration events that we had previously described, we identified 19 smaller clusters. The only cluster with a strong geographic association consisted of Cell-fusing agent virus-like sequences specific to Thailand. The remaining clusters did not have a geographic association and mostly consisted of near identical short sequences without strong similarity to any known flaviviral genomes. The short read sequencing data did not permit us to determine whether identified sequences were extrachromosomal or integrated into Ae. aegypti chromosomes. Our results suggest that Liverpool strain and field Ae. aegypti mosquitoes have a similar variety of conserved flaviviral DNA, whose functional role should be investigated in follow-up studies. [ABSTRACT FROM AUTHOR]

Malaria has exhibited the strongest known selective pressure on the human genome in recent history and is the evolutionary driving force behind genetic conditions, such as sickle-cell disease, glucose-6-phosphatase deficiency, and some other erythrocyte defects. Genomic studies (e.g., The 1000 Genomes project) have provided an invaluable baseline for human genetics, but with an estimated two thousand ethno-linguistic groups thought to exist across the African continent, our understanding of the genetic differences between indigenous populations and their implications on disease is still limited. Low-cost sequencing-based approaches make it possible to target specific molecular markers and genes of interest, leading to potential insights into genetic diversity. Here we demonstrate the versatility of custom dual-indexing technology and Illumina next generation sequencing to generate a genetic profile of human polymorphisms associated with malaria pathology. For 100 individuals diagnosed with severe malaria in Northeast Tanzania, variants were successfully characterised on the haemoglobin subunit beta (HBB), glucose-6-phosphate dehydrogenase (G6PD), atypical chemokine receptor 1 (ACKR1) genes, and the intergenic Dantu genetic blood variant, then validated using pre-existing genotyping data. High sequencing coverage was observed across all amplicon targets in HBB, G6PD, ACKR1, and the Dantu blood group, with variants identified at frequencies previously observed within this region of Tanzania. Sequencing data exhibited high concordance rates to pre-existing genotyping data (> 99.5%). Our work demonstrates the potential utility of amplicon sequencing for applications in human genetics, including to personalise medicine and understand the genetic diversity of loci linked to important host phenotypes, such as malaria susceptibility. [ABSTRACT FROM AUTHOR]

Malaria control initiatives require rapid and reliable methods for the detection and monitoring of molecular markers associated with antimalarial drug resistance in Plasmodium falciparum parasites. Ngodhe island, Kenya, presents a unique malaria profile, with lower P. falciparum incidence rates than the surrounding region, and a high proportion of sub-microscopic and low-density infections. Here, using custom dual-indexing and Illumina next generation sequencing, we generate resistance profiles on seventy asymptomatic and low-density P. falciparum infections from a mass drug administration program implemented on Ngodhe island between 2015 and 2016. Our assay encompasses established molecular markers on the Pfcrt, Pfmdr1, Pfdhps, Pfdhfr, and Pfk13 genes. Resistance markers for sulfadoxine-pyrimethamine were identified at high frequencies, including a quintuple mutant haplotype (Pfdhfr/Pfdhps: N51I, C59R, S108N/A437G, K540E) identified in 62.2% of isolates. The Pfdhps K540E biomarker, used to inform decision making for intermittent preventative treatment in pregnancy, was identified in 79.2% of isolates. Several variants on Pfmdr1, associated with reduced susceptibility to quinolones and lumefantrine, were also identified (Y184F 47.1%; D1246Y 16.0%; N86 98%). Overall, we have presented a low-cost and extendable approach that can provide timely genetic profiles to inform clinical and surveillance activities, especially in settings with abundant low-density infections, seeking malaria elimination. [ABSTRACT FROM AUTHOR]

Motivation Tuberculosis (TB) is caused by members of the Mycobacterium tuberculosis complex (MTBC), which has a strain- or lineage-based clonal population structure. The evolution of drug-resistance in the MTBC poses a threat to successful treatment and eradication of TB. Machine learning approaches are being increasingly adopted to predict drug-resistance and characterize underlying mutations from whole genome sequences. However, such approaches may not generalize well in clinical practice due to confounding from the population structure of the MTBC. Results To investigate how population structure affects machine learning prediction, we compared three different approaches to reduce lineage dependency in random forest (RF) models, including stratification, feature selection, and feature weighted models. All RF models achieved moderate-high performance (area under the ROC curve range: 0.60–0.98). First-line drugs had higher performance than second-line drugs, but it varied depending on the lineages in the training dataset. Lineage-specific models generally had higher sensitivity than global models which may be underpinned by strain-specific drug-resistance mutations or sampling effects. The application of feature weights and feature selection approaches reduced lineage dependency in the model and had comparable performance to unweighted RF models. Availability and implementation https://github.com/NinaMercedes/RF_lineages. [ABSTRACT FROM AUTHOR]

Since its first detection in 2012 in Djibouti, Anopheles stephensi has invaded and established in the Horn of Africa, and more recently Nigeria. The expansion of this vector poses a significant threat to malaria control and elimination efforts. Integrated vector management is the primary strategy used to interrupt disease transmission; however, growing insecticide resistance is threatening to reverse gains in global malaria control. We present a next-generation amplicon-sequencing approach, for high-throughput monitoring of insecticide resistance genes (ace1, GSTe2, vgsc and rdl), species identification and characterization of genetic diversity (its2 and cox1) in An. stephensi. Ninety-five An. stephensi mosquitoes, collected in Ethiopia, were screened, identifying 104 SNPs, including the knock-down mutation L958F (L1014F in Musca domestica), and for the first time in this vector species, the A296S substitution (A301S in Drosophila melanogaster) in the rdl locus. Two other amino acid substitutions (ace1-N177D, GSTe2-V189L) were also identified but have not been previously implicated in insecticide resistance. Genetic diversity in the mitochondrial cox1 gene revealed shared haplotypes between Ethiopian An. stephensi with samples from Pakistan, Sudan, and Djibouti. Overall, we present a reliable, cost-effective strategy using amplicon-sequencing to monitor known insecticide resistance mutations, with the potential to identify new genetic variants, to assist in the high-throughput surveillance of insecticide resistance in An. stephensi populations. [ABSTRACT FROM AUTHOR]

Artificial intelligence (AI)-based approaches can now use electrocardiograms (ECGs) to provide expert-level performance in detecting heart abnormalities and diagnosing disease. Additionally, patient age predicted from ECGs by AI models has shown great potential as a biomarker for cardiovascular age, where recent work has found its deviation from chronological age ("delta age") to be associated with mortality and co-morbidities. However, despite being crucial for understanding underlying individual risk, the genetic underpinning of delta age is unknown. In this work we performed a genome-wide association study using UK Biobank data (n=34,432) and identified eight loci associated with delta age ( p ≤ 5 × 10 - 8 ), including genes linked to cardiovascular disease (CVD) (e.g. SCN5A) and (heart) muscle development (e.g. TTN). Our results indicate that the genetic basis of cardiovascular ageing is predominantly determined by genes directly involved with the cardiovascular system rather than those connected to more general mechanisms of ageing. Our insights inform the epidemiology of CVD, with implications for preventative and precision medicine. [ABSTRACT FROM AUTHOR]

Aedes mosquito vectors transmit many viruses of global health concern, including dengue, chikungunya and Zika. These vector-borne viral diseases have a limited number of treatment options, and vaccines vary in their effectiveness. Consequently, integrated vector management is a primary strategy for disease control. However, the increasing emergence and spread of insecticide resistance is threatening the efficacy of vector control methods. Identifying mutations associated with resistance in vector populations is important to monitor the occurrence and evolution of insecticide resistance and inform control strategies. Rapid and cost-effective genome sequencing approaches are urgently needed. Here we present an adaptable targeted amplicon approach for cost-effective implementation within next generation sequencing platforms. This approach can identify single nucleotide polymorphisms (SNPs) and small insertions and deletions (indels) in genes involved in insecticide resistance in Aedes aegypti mosquitoes. We designed and tested eleven amplicons, which included segments of the ace-1 (carbamate target), the Voltage-Gated Sodium Channel (vgsc; pyrethroids, DDT and organochlorines), and rdl (dieldrin) genes; thereby covering established knockdown resistance (kdr) mutations (e.g., S989P, I1011M/V, V1016G/I and F1534C), with the potential to identify novel ones. The amplicon assays were designed with internal barcodes, to facilitate multiplexing of large numbers of mosquitoes at low cost, and were sequenced using an Illumina platform. Our approach was evaluated on 152 Ae. aegypti mosquitoes collected in Cabo Verde, an archipelago with a history of arbovirus outbreaks. The amplicon sequence data revealed 146 SNPs, including four non-synonymous polymorphisms in the vgsc gene, one in ace-1 and the 296S rdl mutation previously associated with resistance to organochlorines. The 296S rdl mutation was identified in 98% of mosquitoes screened, consistent with the past use of an organochlorine compound (e.g., DDT). Overall, our work shows that targeted amplicon sequencing is a rapid, robust, and cost-effective tool that can be used to perform high throughput monitoring of insecticide resistance. Author summary: Many viruses, such as dengue and Zika, are transmitted by Aedes aegypti mosquitoes. These vector-borne diseases are a major public health problem in the tropics and sub-tropics worldwide. The primary strategy to reduce their burden is vector control, including through the application of insecticides. However, many mosquito populations have developed resistance to insecticides. Here, we present a rapid, robust, and cost-effective tool that allows the screening of genes associated with insecticide resistance across many Aedes aegypti mosquitoes, identifying known resistance and novel mutations. This assay will support vector control strategies by informing on the emergence and spread of insecticide resistance mutations across Aedes aegypti populations. [ABSTRACT FROM AUTHOR]

Mosquitoes of the genus Aedes are the main vectors of many viruses, e.g. dengue and Zika, which affect millions of people each year and for which there are limited treatment options. Understanding how Aedes mosquitoes tolerate high viral loads may lead to better disease control strategies. Elucidating endogenous viral elements (EVEs) within vector genomes may give exploitable biological insights. Previous studies have reported the presence of a large number of EVEs in Aedes genomes. Here we investigated if flavivirus EVEs are conserved across populations and different Aedes species by using ~ 500 whole genome sequence libraries from Aedes aegypti and Aedes albopictus, sourced from colonies and field mosquitoes across continents. We found that nearly all flavivirus EVEs in the Ae. aegypti reference genome originate from four separate putative viral integration events, and that they are highly conserved across geographically diverse samples. By contrast, flavivirus EVEs in the Ae. albopictus reference genome originate from up to nine distinct integration events and show low levels of conservation, even within samples from narrow geographical ranges. Our analysis suggests that flaviviruses integrated as long sequences and were subsequently fragmented and shuffled by transposable elements. Given that EVEs of Ae. aegypti and Ae. albopictus belong to different phylogenetic clades and have very differing levels of conservation, they may have different evolutionary origins and potentially different functional roles. [ABSTRACT FROM AUTHOR]

Tuberculosis (TB), caused by Mycobacterium tuberculosis, is endemic in Pakistan. Resistance to both firstline rifampicin and isoniazid drugs (multidrug-resistant TB; MDR-TB) is hampering disease control. Rifampicin resistance is attributed to rpoB gene mutations, but rpoA and rpoC loci may also be involved. To characterise underlying rifampicin resistance mutations in the TB endemic province of Khyber Pakhtunkhwa, we sequenced 51 M. tuberculosis isolates collected between 2016 and 2019; predominantly, MDR-TB (n = 44; 86.3%) and lineage 3 (n = 30, 58.8%) strains. We found that known mutations in rpoB (e.g. S405L), katG (e.g. S315T), or inhA promoter loci explain the MDR-TB. There were 24 unique mutations in rpoA, rpoB, and rpoC genes, including four previously unreported. Five instances of within-host resistance diversity were observed, where two were a mixture of MDR-TB strains containing mutations in rpoB, katG, and the inhA promoter region, as well as compensatory mutations in rpoC. Heteroresistance was observed in two isolates with a single lineage. Such complexity may reflect the high transmission nature of the Khyber Pakhtunkhwa setting. Our study reinforces the need to apply sequencing approaches to capture the full-extent of MDR-TB genetic diversity, to understand transmission, and to inform TB control activities in the highly endemic setting of Pakistan. [ABSTRACT FROM AUTHOR]

Tuberculosis (TB), caused by Mycobacterium tuberculosis, is endemic in Pakistan. Resistance to both firstline rifampicin and isoniazid drugs (multidrug-resistant TB; MDR-TB) is hampering disease control. Rifampicin resistance is attributed to rpoB gene mutations, but rpoA and rpoC loci may also be involved. To characterise underlying rifampicin resistance mutations in the TB endemic province of Khyber Pakhtunkhwa, we sequenced 51 M. tuberculosis isolates collected between 2016 and 2019; predominantly, MDR-TB (n = 44; 86.3%) and lineage 3 (n = 30, 58.8%) strains. We found that known mutations in rpoB (e.g. S405L), katG (e.g. S315T), or inhA promoter loci explain the MDR-TB. There were 24 unique mutations in rpoA, rpoB, and rpoC genes, including four previously unreported. Five instances of within-host resistance diversity were observed, where two were a mixture of MDR-TB strains containing mutations in rpoB, katG, and the inhA promoter region, as well as compensatory mutations in rpoC. Heteroresistance was observed in two isolates with a single lineage. Such complexity may reflect the high transmission nature of the Khyber Pakhtunkhwa setting. Our study reinforces the need to apply sequencing approaches to capture the full-extent of MDR-TB genetic diversity, to understand transmission, and to inform TB control activities in the highly endemic setting of Pakistan. [ABSTRACT FROM AUTHOR]

Mosquitoes of the genus Aedes are the main vectors of many viruses, e.g. dengue and Zika, which affect millions of people each year and for which there are limited treatment options. Understanding how Aedes mosquitoes tolerate high viral loads may lead to better disease control strategies. Elucidating endogenous viral elements (EVEs) within vector genomes may give exploitable biological insights. Previous studies have reported the presence of a large number of EVEs in Aedes genomes. Here we investigated if flavivirus EVEs are conserved across populations and different Aedes species by using ~ 500 whole genome sequence libraries from Aedes aegypti and Aedes albopictus, sourced from colonies and field mosquitoes across continents. We found that nearly all flavivirus EVEs in the Ae. aegypti reference genome originate from four separate putative viral integration events, and that they are highly conserved across geographically diverse samples. By contrast, flavivirus EVEs in the Ae. albopictus reference genome originate from up to nine distinct integration events and show low levels of conservation, even within samples from narrow geographical ranges. Our analysis suggests that flaviviruses integrated as long sequences and were subsequently fragmented and shuffled by transposable elements. Given that EVEs of Ae. aegypti and Ae. albopictus belong to different phylogenetic clades and have very differing levels of conservation, they may have different evolutionary origins and potentially different functional roles. [ABSTRACT FROM AUTHOR]

Background: Approximately 10 % of the Mycobacterium tuberculosis genome is made up of two families of genes that are poorly characterized due to their high GC content and highly repetitive nature. The PE and PPE families are typified by their highly conserved N-terminal domains that incorporate proline-glutamate (PE) and proline-proline-glutamate (PPE) signature motifs. They are hypothesised to be important virulence factors involved with host-pathogen interactions, but their high genetic variability and complexity of analysis means they are typically disregarded in genome studies. Results: To elucidate the structure of these genes, 518 genomes from a diverse international collection of clinical isolates were de novo assembled. A further 21 reference M. tuberculosis complex genomes and long read sequence data were used to validate the approach. SNP analysis revealed that variation in the majority of the 168 pe/ppe genes studied was consistent with lineage. Several recombination hotspots were identified, notably pe_pgrs3 and pe_pgrs17. Evidence of positive selection was revealed in 65 pe/ppe genes, including epitopes potentially binding to major histocompatibility complex molecules. Conclusions: This, the first comprehensive study of the pe and ppe genes, provides important insight into M. tuberculosis diversity and has significant implications for vaccine development. [ABSTRACT FROM AUTHOR]

Resistance to isoniazid (INH) and rifampicin (RIF) first-line drugs in Mycobacterium tuberculosis (Mtb), together called multi-drug resistance, threatens tuberculosis control. Resistance mutations in katG (for INH) and rpoB (RIF) genes often come with fitness costs. To overcome these costs, Mtb compensatory mutations have arisen in rpoC/rpoA (RIF) and ahpC (INH) loci. By leveraging the presence of known compensatory mutations, we aimed to detect novel resistance mutations occurring in INH and RIF target genes. Across ~ 32 k Mtb isolates with whole genome sequencing (WGS) data, there were 6262 (35.7%) with INH and 5435 (30.7%) with RIF phenotypic resistance. Known mutations in katG and rpoB explained ~ 99% of resistance. However, 188 (0.6%) isolates had ahpC compensatory mutations with no known resistance mutations in katG, leading to the identification of 31 putative resistance mutations in katG, each observed in at least 3 isolates. These putative katG mutations can co-occur with other INH variants (e.g., katG-Ser315Thr, fabG1 mutations). For RIF, there were no isolates with rpoC/rpoA compensatory mutations and unknown resistance mutations. Overall, using WGS data we identified putative resistance markers for INH that could be used for genotypic drug-resistance profiling. Establishing the complete repertoire of Mtb resistance mutations will assist the clinical management of tuberculosis. [ABSTRACT FROM AUTHOR]

Malaria, caused by Plasmodium parasites, is a major global health challenge. Whole genome sequencing (WGS) of Plasmodium falciparum and Plasmodium vivax genomes is providing insights into parasite genetic diversity, transmission patterns, and can inform decision making for clinical and surveillance purposes. Advances in sequencing technologies are helping to generate timely and big genomic datasets, with the prospect of applying Artificial Intelligence analytical techniques (e.g., machine learning) to support programmatic malaria control and elimination. Here, we assess the potential of applying deep learning convolutional neural network approaches to predict the geographic origin of infections (continents, countries, GPS locations) using WGS data of P. falciparum (n = 5957; 27 countries) and P. vivax (n = 659; 13 countries) isolates. Using identified high-quality genome-wide single nucleotide polymorphisms (SNPs) (P. falciparum: 750 k, P. vivax: 588 k), an analysis of population structure and ancestry revealed clustering at the country-level. When predicting locations for both species, classification (compared to regression) methods had the lowest distance errors, and > 90% accuracy at a country level. Our work demonstrates the utility of machine learning approaches for geo-classification of malaria parasites. With timelier WGS data generation across more malaria-affected regions, the performance of machine learning approaches for geo-classification will improve, thereby supporting disease control activities. [ABSTRACT FROM AUTHOR]

With >1 million associated deaths in 2020, human tuberculosis (TB) caused by the bacteria Mycobacterium tuberculosis remains one of the deadliest infectious diseases. A plethora of genomic tools and bioinformatics pipelines have become available in recent years to assist the whole genome sequencing of M. tuberculosis. The Oxford Nanopore Technologies (ONT) portable sequencer is a promising platform for cost-effective application in clinics, including personalizing treatment through detection of drug resistance-associated mutations, or in the field, to assist epidemiological and transmission investigations. In this study, we performed a comparison of 10 clinical isolates with DNA sequenced on both long-read ONT and (gold standard) short-read Illumina HiSeq platforms. Our analysis demonstrates the robustness of the ONT variant calling for single nucleotide polymorphisms, despite the high error rate. Moreover, because of improved coverage in repetitive regions where short sequencing reads fail to align accurately, ONT data analysis can incorporate additional regions of the genome usually excluded (e.g. pe / ppe genes). The resulting extra resolution can improve the characterization of transmission clusters and dynamics based on inferring closely related isolates. High concordance in variants in loci associated with drug resistance supports its use for the rapid detection of resistant mutations. Overall, ONT sequencing is a promising tool for TB genomic investigations, particularly to inform clinical and surveillance decision-making to reduce the disease burden. [ABSTRACT FROM AUTHOR]

Tuberculosis (TB), caused by Mycobacterium tuberculosis bacteria, is a leading infectious cause of mortality worldwide, including in Pakistan. Drug resistant M. tuberculosis is an emerging threat for TB control, making it important to detect the underlying genetic mutations, and thereby inform treatment decision making and prevent transmission. Whole genome sequencing has emerged as the new diagnostic to reliably predict drug resistance within a clinically relevant time frame, and its deployment will have the greatest impact on TB control in highly endemic regions. To evaluate the mutations leading to drug resistance and to assess for evidence of the transmission of resistant strains, 81 M. tuberculosis samples from Khyber Pakhtunkhwa province (North West Pakistan) were subjected to whole genome sequencing and standard drug susceptibility testing for eleven anti-TB drugs. We found the majority of M. tuberculosis isolates were the CAS/Delhi strain-type (lineage 3; n = 57; 70.4%) and multi-drug resistant (MDR; n = 62; 76.5%). The most frequent resistance mutations were observed in the katG and rpoB genes, conferring resistance to isoniazid and rifampicin respectively. Mutations were also observed in genes conferring resistance to other first and second-line drugs, including in pncA (pyrazinamide), embB (ethambutol), gyrA (fluoroquinolones), rrs (aminoglycosides), rpsL, rrs and giB (streptomycin) loci. Whilst the majority of mutations have been reported in global datasets, we describe unreported putative resistance markers in katG, ethA (ethionamide), gyrA and gyrB (fluoroquinolones), and pncA. Analysis of the mutations revealed that acquisition of rifampicin resistance often preceded isoniazid in our isolates. We also observed a high proportion (17.6%) of pre-MDR isolates with fluoroquinolone resistance markers, potentially due to unregulated anti-TB drug use. Our isolates were compared to previously sequenced strains from Pakistan in a combined phylogenetic tree analysis. The presence of lineage 2 was only observed in our isolates. Using a cut-off of less than ten genome-wide mutation differences between isolates, a transmission analysis revealed 18 M. tuberculosis isolates clustering within eight networks, thereby providing evidence of drug-resistant TB transmission in the Khyber Pakhtunkhwa province. Overall, we have demonstrated that drug-resistant TB isolates are circulating and transmitted in North West Pakistan. Further, we have shown the usefulness of whole genome sequencing as a diagnostic tool for characterizing M. tuberculosis isolates, which will assist future epidemiological studies and disease control activities in Pakistan. [ABSTRACT FROM AUTHOR]

Whole genome sequencing (WGS) of Mycobacterium tuberculosis offers valuable insights for tuberculosis (TB) control. High throughput platforms like Illumina and Oxford Nanopore Technology (ONT) are increasingly used globally, although ONT is known for higher error rates and is less established for genomic studies. Here we present a study comparing the sequencing outputs of both Illumina and ONT platforms, analysing DNA from 59 clinical isolates in highly endemic TB regions of Thailand. The resulting sequence data were used to profile the M. tuberculosis pairs for their lineage, drug resistance and presence in transmission chains, and were compared to publicly available WGS data from Thailand (n = 1456). Our results revealed isolates that are predominantly from lineages 1 and 2, with consistent drug resistance profiles, including six multidrug-resistant strains; however, analysis of ONT data showed longer phylogenetic branches, emphasising the technologies higher error rate. An analysis incorporating the larger dataset identified fifteen of our samples within six potential transmission clusters, including a significant clade of 41 multi-drug resistant isolates. ONT's extended sequences also revealed strain-specific structural variants in pe/ppe genes (e.g. ppe50), which are candidate loci for vaccine development. Despite some limitations, our results show that ONT sequencing is a promising approach for TB genomic research, supporting precision medicine and decision-making in areas with less developed infrastructure, which is crucial for tackling the disease's significant regional burden. [ABSTRACT FROM AUTHOR]

Analysis of genome sequencing data from >100,000 genomes of Mycobacterium tuberculosis complex using TB-Annotator software revealed a previously unknown lineage, proposed name L10, in central Africa. Phylogenetic reconstruction suggests L10 could represent a missing link in the evolutionary and geographic migration histories of M. africanum. [ABSTRACT FROM AUTHOR]

The Philippines has a high incidence of tuberculosis disease (TB), with an increasing prevalence of multidrug-resistant Mycobacterium tuberculosis (MDR-TB) strains making its control difficult. Although the M. tuberculosis "Manila" ancient lineage 1 strain-type is thought to be prevalent in the country, with evidence of export to others, little is known about the genetic diversity of circulating strains. By whole genome sequencing (WGS) 178 isolates from the Philippines National Drug Resistance Survey, we found the majority (143/178; 80.3%) belonged to the lineage 1 Manila clade, with the minority belonging to lineages 4 (European-American; n = 33) and 2 (East Asian; n = 2). A high proportion were found to be multidrug-resistant (34/178; 19.1%), established through highly concordant laboratory drug susceptibility testing and in silico prediction methods. Some MDR-TB isolates had near identical genomic variation, providing potential evidence of transmission. By placing the Philippine isolates within a phylogeny of global M. tuberculosis (n > 17,000), we established that they are genetically similar to those observed outside the country, including a clade of Manila-like strain-types in Thailand. An analysis of the phylogeny revealed a set of ~200 SNPs that are specific for the Manila strain-type, and a subset can be used within a molecular barcode. Sixty-eight mutations known to be associated with 10 anti-TB drug resistance were identified in the Philippine strains, and all have been observed in other populations. Whilst nine putative streptomycin resistance conferring markers in gid (8) and rrs (1) genes appear to be novel and with functional consequences. Overall, this study provides an important baseline characterisation of M. tuberculosis genetic diversity for the Philippines, and will fill a gap in global datasets and aid the development of a nation-wide database for epidemiological studies and clinical decision making. Further, by establishing a molecular barcode for detecting Manila strains it will assist with the design of diagnostic tools for disease control activities. [ABSTRACT FROM AUTHOR]

Background: Mycobacterium tuberculosis resistance to anti-tuberculosis drugs is a major threat to global public health. Whole genome sequencing (WGS) is rapidly gaining traction as a diagnostic tool for clinical tuberculosis settings. To support this informatically, previous work led to the development of the widely used TBProfiler webtool, which predicts resistance to 14 drugs from WGS data. However, for accurate and rapid high throughput of samples in clinical or epidemiological settings, there is a need for a stand-alone tool and the ability to analyse data across multiple WGS platforms, including Oxford Nanopore MinION. Results: We present a new command line version of the TBProfiler webserver, which includes hetero-resistance calling and will facilitate the batch processing of samples. The TBProfiler database has been expanded to incorporate 178 new markers across 16 anti-tuberculosis drugs. The predictive performance of the mutation library has been assessed using > 17,000 clinical isolates with WGS and laboratory-based drug susceptibility testing (DST) data. An integrated MinION analysis pipeline was assessed by performing WGS on 34 replicates across 3 multi-drug resistant isolates with known resistance mutations. TBProfiler accuracy varied by individual drug. Assuming DST as the gold standard, sensitivities for detecting multi-drug-resistant TB (MDR-TB) and extensively drug-resistant TB (XDR-TB) were 94% (95%CI 93–95%) and 83% (95%CI 79–87%) with specificities of 98% (95%CI 98–99%) and 96% (95%CI 95–97%) respectively. Using MinION data, only one resistance mutation was missed by TBProfiler, involving an insertion in the tlyA gene coding for capreomycin resistance. When compared to alternative platforms (e.g. Mykrobe predictor TB, the CRyPTIC library), TBProfiler demonstrated superior predictive performance across first- and second-line drugs. Conclusions: The new version of TBProfiler can rapidly and accurately predict anti-TB drug resistance profiles across large numbers of samples with WGS data. The computing architecture allows for the ability to modify the core bioinformatic pipelines and outputs, including the analysis of WGS data sourced from portable technologies. TBProfiler has the potential to be integrated into the point of care and WGS diagnostic environments, including in resource-poor settings. [ABSTRACT FROM AUTHOR]

Genomic studies of Mycobacterium tuberculosis bacteria have revealed loci associated with resistance to anti-tuberculosis drugs. However, the molecular consequences of polymorphism within these candidate loci remain poorly understood. To address this, we have used computational tools to quantify the effects of point mutations conferring resistance to three major anti-tuberculosis drugs, isoniazid (n = 189), rifampicin (n = 201) and D-cycloserine (n = 48), within their primary targets, katG, rpoB, and alr. Notably, mild biophysical effects brought about by high incidence mutations were considered more tolerable, while different structural effects brought about by haplotype combinations reflected differences in their functional importance. Additionally, highly destabilising mutations such as alr Y388, highlighted a functional importance of the wildtype residue. Our qualitative analysis enabled us to relate resistance mutations onto a theoretical landscape linking enthalpic changes with phenotype. Such insights will aid the development of new resistance-resistant drugs and, via an integration into predictive tools, in pathogen surveillance. [ABSTRACT FROM AUTHOR]

During COVID-19, diagnostic serological tools and vaccines have been developed. To inform control activities in a post-vaccine surveillance setting, we have developed an online "immuno-analytics" resource that combines epitope, sequence, protein and SARS-CoV-2 mutation analysis. SARS-CoV-2 spike and nucleocapsid proteins are both vaccine and serological diagnostic targets. Using the tool, the nucleocapsid protein appears to be a sub-optimal target for use in serological platforms. Spike D614G (and nsp12 L314P) mutations were most frequent (> 86%), whilst spike A222V/L18F have recently increased. Also, Orf3a proteins may be a suitable target for serology. The tool can accessed from: http://genomics.lshtm.ac.uk/immuno (online); https://github.com/dan-ward-bio/COVID-immunoanalytics (source code). [ABSTRACT FROM AUTHOR]

Vector control strategies have been successful in reducing the number of malaria cases and deaths globally, but the spread of insecticide resistance represents a significant threat to disease control. Insecticide resistance has been reported across Anopheles (An.) vector populations, including species within the An. funestus group. These mosquitoes are responsible for intense malaria transmission across sub-Saharan Africa, including in the Democratic Republic of the Congo (DRC), a country contributing > 12% of global malaria infections and mortality events. To support the continuous efficacy of vector control strategies, it is essential to monitor insecticide resistance using molecular surveillance tools. In this study, we developed an amplicon sequencing ("Amp-seq") approach targeting An. funestus, and using multiplex PCR, dual index barcoding, and next-generation sequencing for high throughput and low-cost applications. Using our Amp-seq approach, we screened 80 An. funestus field isolates from the DRC across a panel of nine genes with mutations linked to insecticide resistance (ace-1, CYP6P4, CYP6P9a, GSTe2, vgsc, and rdl) and mosquito speciation (cox-1, mtND5, and ITS2). Amongst the 18 non-synonymous mutations detected, was N485I, in the ace-1 gene associated with carbamate resistance. Overall, our panel represents an extendable and much-needed method for the molecular surveillance of insecticide resistance in An. funestus populations. [ABSTRACT FROM AUTHOR]

Background Carbapenem-resistant Klebsiella pneumoniae (CRKP) strains are of particular concern, especially strains with mobilizable carbapenemase genes such as bla KPC, bla NDM or bla OXA-48, given that carbapenems are usually the last line drugs in the β-lactam class and, resistance to this sub-class is associated with increased mortality and frequently co-occurs with resistance to other antimicrobial classes. Objectives To characterize the genomic diversity and international dissemination of CRKP strains from tertiary care hospitals in Lisbon, Portugal. Methods Twenty CRKP isolates obtained from different patients were subjected to WGS for species confirmation, typing, drug resistance gene detection and phylogenetic reconstruction. Two additional genomic datasets were included for comparative purposes: 26 isolates (ST13, ST17 and ST231) from our collection and 64 internationally available genomic assemblies (ST13). Results By imposing a 21 SNP cut-off on pairwise comparisons we identified two genomic clusters (GCs): ST13/GC1 (n  = 11), all bearing bla KPC-3, and ST17/GC2 (n  = 4) harbouring bla OXA-181 and bla CTX-M-15 genes. The inclusion of the additional datasets allowed the expansion of GC1/ST13/KPC-3 to 23 isolates, all exclusively from Portugal, France and the Netherlands. The phylogenetic tree reinforced the importance of the GC1/KPC-3-producing clones along with their rapid emergence and expansion across these countries. The data obtained suggest that the ST13 branch emerged over a decade ago and only more recently did it underpin a stronger pulse of transmission in the studied population. Conclusions This study identifies an emerging OXA-181/ST17-producing strain in Portugal and highlights the ongoing international dissemination of a KPC-3/ST13-producing clone from Portugal. [ABSTRACT FROM AUTHOR]

Background: Second-line drug resistance (SLD) among tuberculosis (TB) patients is a serious emerging challenge towards global control of the disease. We characterized SLD-resistance conferring-mutations among TB patients with rifampicin and/or isoniazid (RIF and/or INH) drug-resistance tested at the Uganda National TB Reference Laboratory (NTRL) between June 2017 and December 2019.Methods: This was a descriptive cross-sectional secondary data analysis of 20,508 M. tuberculosis isolates of new and previously treated patients' resistant to RIF and/or INH. DNA strips with valid results to characterise the SLD resistance using the commercial Line Probe Assay Genotype MTBDRsl Version 2.0 Assay (Hain Life Science, Nehren, Germany) were reviewed. Data were analysed with STATAv15 using cross-tabulation for frequency and proportions of known resistance-conferring mutations to injectable agents (IA) and fluoroquinolones (FQ).Results: Among the eligible participants, 12,993/20,508 (63.4%) were male and median (IQR) age 32 (24-43). A total of 576/20,508 (2.8%) of the M. tuberculosis isolates from participants had resistance to RIF and/or INH. These included; 102/576 (17.7%) single drug-resistant and 474/576 (82.3%) multidrug-resistant (MDR) strains. Only 102 patients had test results for FQ of whom 70/102 (68.6%) and 01/102 (0.98%) had resistance-conferring mutations in the gyrA locus and gyrB locus respectively. Among patients with FQ resistance, gyrAD94G 42.6% (30.0-55.9) and gyrA A90V 41.1% (28.6-54.3) mutations were most observed. Only one mutation, E540D was detected in the gyrB locus. A total of 26 patients had resistance-conferring mutations to IA in whom, 20/26 77.0% (56.4-91.0) had A1401G mutation in the rrs gene locus.Conclusions: Our study reveals a high proportion of mutations known to confer high-level fluoroquinolone drug-resistance among patients with rifampicin and/or isoniazid drug resistance. Utilizing routinely generated laboratory data from existing molecular diagnostic methods may aid real-time surveillance of emerging tuberculosis drug-resistance in resource-limited settings. [ABSTRACT FROM AUTHOR]

Multidrug-resistant tuberculosis (MDR TB), pre-extensively drug-resistant tuberculosis (pre-XDR TB), and extensively drug-resistant tuberculosis (XDR TB) complicate disease control. We analyzed whole-genome sequence data for 579 phenotypically drug-resistant M. tuberculosis isolates (28% of available MDR/pre-XDR and all culturable XDR TB isolates collected in Thailand during 2014-2017). Most isolates were from lineage 2 (n = 482; 83.2%). Cluster analysis revealed that 281/579 isolates (48.5%) formed 89 clusters, including 205 MDR TB, 46 pre-XDR TB, 19 XDR TB, and 11 poly-drug-resistant TB isolates based on genotypic drug resistance. Members of most clusters had the same subset of drug resistance-associated mutations, supporting potential primary resistance in MDR TB (n = 176/205; 85.9%), pre-XDR TB (n = 29/46; 63.0%), and XDR TB (n = 14/19; 73.7%). Thirteen major clades were significantly associated with geography (p<0.001). Clusters of clonal origin contribute greatly to the high prevalence of drug-resistant TB in Thailand. [ABSTRACT FROM AUTHOR]

Background: Tuberculosis (TB), particularly multi- and or extensive drug resistant TB, is still a global medical emergency. Whole genome sequencing (WGS) is a current alternative to the WHO-approved probe-based methods for TB diagnosis and detection of drug resistance, genetic diversity and transmission dynamics of Mycobacterium tuberculosis complex (MTBC). This study compared WGS and clinical data in participants with TB. Results: This cohort study performed WGS on 87 from MTBC DNA isolates, 57 (66%) and 30 (34%) patients with drug resistant and susceptible TB, respectively. Drug resistance was determined by Xpert® MTB/RIF assay and phenotypic culture-based drug-susceptibility-testing (DST). WGS and bioinformatics data that predict phenotypic resistance to anti-TB drugs were compared with participant's clinical outcomes. They were 47 female participants (54%) and the median age was 35 years (IQR): 29–44). Twenty (23%) and 26 (30%) of participants had TB/HIV co-infection BMI < 18 kg/m2 respectively. MDR-TB participants had MTBC with multiple mutant genes, compared to those with mono or polyresistant TB, and the majority belonged to lineage 3 Central Asian Strain (CAS). Also, MDR-TB was associated with delayed culture-conversion (median: IQR (83: 60–180 vs. 51:30–66) days). WGS had high concordance with both culture-based DST and Xpert® MTB/RIF assay in detecting drug resistance (kappa = 1.00). Conclusion: This study offers comparison of mutations detected by Xpert and WGS with phenotypic DST of M. tuberculosis isolates in Tanzania. The high concordance between the different methods and further insights provided by WGS such as PZA-DST, which is not routinely performed in most resource-limited-settings, provides an avenue for inclusion of WGS into diagnostic matrix of TB including drug-resistant TB. [ABSTRACT FROM AUTHOR]

Portugal is a low incidence country for tuberculosis (TB) disease. Now figuring among TB low incidence countries, it has since the 1990s reported multidrug resistant and extensively drug resistant (XDR) TB cases, driven predominantly by two strain-types: Lisboa3 and Q1. This study describes the largest characterization of the evolutionary trajectory of M/XDR-TB strains in Portugal, spanning a time-period of two decades. By combining whole-genome sequencing and phenotypic susceptibility data for 207 isolates, we report the geospatial patterns of drug resistant TB, particularly the dispersion of Lisboa3 and Q1 clades, which underly 64.2% and 94.0% of all MDR-TB and XDR-TB isolates, respectively. Genomic-based similarity and a phylogenetic analysis revealed multiple clusters (n = 16) reflecting ongoing and uncontrolled recent transmission of M/XDR-TB, predominantly associated with the Lisboa3 and Q1 clades. These clades are now thought to be evolving in a polycentric mode across multiple geographical districts. The inferred evolutionary history is compatible with MDR- and XDR-TB originating in Portugal in the 70's and 80's, respectively, but with subsequent multiple emergence events of MDR and XDR-TB particularly involving the Lisboa3 clade. A SNP barcode was defined for Lisboa3 and Q1 and comparison with a phylogeny of global strain-types (n = 28 385) revealed the presence of Lisboa3 and Q1 strains in Europe, South America and Africa. In summary, Portugal displays an unusual and unique epidemiological setting shaped by >40 years of uncontrolled circulation of two main phylogenetic clades, leading to a sympatric evolutionary trajectory towards XDR-TB with the potential for global reach. [ABSTRACT FROM AUTHOR]

Human tuberculosis disease (TB), caused by Mycobacterium tuberculosis (Mtb), is a complex disease, with a spectrum of outcomes. Genomic, transcriptomic and methylation studies have revealed differences between Mtb lineages, likely to impact on transmission, virulence and drug resistance. However, so far no studies have integrated sequence-based genomic, transcriptomic and methylation characterisation across a common set of samples, which is critical to understand how DNA sequence and methylation affect RNA expression and, ultimately, Mtb pathogenesis. Here we perform such an integrated analysis across 22 M. tuberculosis clinical isolates, representing ancient (lineage 1) and modern (lineages 2 and 4) strains. The results confirm the presence of lineage-specific differential gene expression, linked to specific SNP-based expression quantitative trait loci: with 10 eQTLs involving SNPs in promoter regions or transcriptional start sites; and 12 involving potential functional impairment of transcriptional regulators. Methylation status was also found to have a role in transcription, with evidence of differential expression in 50 genes across lineage 4 samples. Lack of methylation was associated with three novel variants in mamA, likely to cause loss of function of this enzyme. Overall, our work shows the relationship of DNA sequence and methylation to RNA expression, and differences between ancient and modern lineages. Further studies are needed to verify the functional consequences of the identified mechanisms of gene expression regulation. [ABSTRACT FROM AUTHOR]

Background: Mixed, polyclonal Mycobacterium tuberculosis infection occurs in natural populations. Developing an effective method for detecting such cases is important in measuring the success of treatment and reconstruction of transmission between patients. Using whole genome sequence (WGS) data, we assess two methods for detecting mixed infection: (i) a combination of the number of heterozygous sites and the proportion of heterozygous sites to total SNPs, and (ii) Bayesian model-based clustering of allele frequencies from sequencing reads at heterozygous sites. Results: In silico and in vitro artificially mixed and known pure M. tuberculosis samples were analysed to determine the specificity and sensitivity of each method. We found that both approaches were effective in distinguishing between pure strains and mixed infection where there was relatively high (> 10%) proportion of a minor strain in the mixture. A large dataset of clinical isolates (n = 1963) from the Karonga Prevention Study in Northern Malawi was tested to examine correlations with patient characteristics and outcomes with mixed infection. The frequency of mixed infection in the population was found to be around 10%, with an association with year of diagnosis, but no association with age, sex, HIV status or previous tuberculosis. Conclusions: Mixed Mycobacterium tuberculosis infection was identified in silico using whole genome sequence data. The methods presented here can be applied to population-wide analyses of tuberculosis to estimate the frequency of mixed infection, and to identify individual cases of mixed infections. These cases are important when considering the evolution and transmission of the disease, and in patient treatment. [ABSTRACT FROM AUTHOR]