Jody Phelan

Mycobacterium tuberculosis whole genome sequencing and protein structure modelling provides insights into anti-tuberculosis drug resistance

BackgroundCombating the spread of drug resistant tuberculosis is a global health priority. Whole genome association studies are being applied to identify genetic determinants of resistance to anti-tuberculosis drugs. Protein structure and interaction modelling are used to understand the functional effects of putative mutations and provide insight into the molecular mechanisms leading to resistance.MethodsTo investigate the potential utility of these approaches, we analysed the genomes of 144 Mycobacterium tuberculosis clinical isolates from The Special Programme for Research and Training in Tropical Diseases (TDR) collection sourced from 20 countries in four continents. A genome-wide approach was applied to 127 isolates to identify polymorphisms associated with minimum inhibitory concentrations for first-line anti-tuberculosis drugs. In addition, the effect of identified candidate mutations on protein stability and interactions was assessed quantitatively with well-established computational methods.ResultsThe analysis revealed that mutations in the genes rpoB (rifampicin), katG (isoniazid), inhA-promoter (isoniazid), rpsL (streptomycin) and embB (ethambutol) were responsible for the majority of resistance observed. A subset of the mutations identified in rpoB and katG were predicted to affect protein stability. Further, a strong direct correlation was observed between the minimum inhibitory concentration values and the distance of the mutated residues in the three-dimensional structures of rpoB and katG to their respective drugs binding sites.ConclusionsUsing the TDR resource, we demonstrate the usefulness of whole genome association and convergent evolution approaches to detect known and potentially novel mutations associated with drug resistance. Further, protein structural modelling could provide a means of predicting the impact of polymorphisms on drug efficacy in the absence of phenotypic data. These approaches could ultimately lead to novel resistance mutations to improve the design of tuberculosis control measures, such as diagnostics, and inform patient management.

Outcomes, infectiousness, and transmission dynamics of patients with extensively drug-resistant tuberculosis and home-discharged patients with programmatically incurable tuberculosis: a prospective cohort study

BACKGROUND The emergence of programmatically incurable tuberculosis threatens to destabilise control efforts. The aim of this study was to collect prospective patient-level data to inform treatment and containment strategies. METHODS In a prospective cohort study, 273 South African patients with extensively drug-resistant tuberculosis, or resistance beyond extensively drug-resistant tuberculosis, were followed up over a period of 6 years. Transmission dynamics, infectiousness, and drug susceptibility were analysed in a subset of patients from the Western Cape using whole-genome sequencing (WGS; n=149), a cough aerosol sampling system (CASS; n=26), and phenotypic testing for 18 drugs (n=179). FINDINGS Between Oct 1, 2008, and Oct 31, 2012, we enrolled and followed up 273 patients for a median of 20·3 months (IQR 9·6-27·8). 203 (74%) had programmatically incurable tuberculosis and unfavourable outcomes (treatment failure, relapse, default, or death despite treatment with a regimen based on capreomycin, aminosalicylic acid, or both). 172 (63%) patients were discharged home, of whom 104 (60%) had an unfavourable outcome. 54 (31%) home-discharged patients had failed treatment, with a median time to death after discharge of 9·9 months (IQR 4·2-17·4). 35 (20%) home-discharged cases were smear-positive at discharge. Using CASS, six (23%) of 26 home-discharged cases with data available expectorated infectious culture-positive cough aerosols in the respirable range (<5 μm), and most reported inter-person contact with suboptimal protective mask usage. WGS identified 17 (19%) of the 90 patients (with available sequence data) that were discharged home before the diagnosis of 20 downstream cases of extensively drug-resistant tuberculosis with almost identical sequencing profiles suggestive of community-based transmission (five or fewer single nucleotide polymorphisms different and with identical resistance-encoding mutations for 14 drugs). 11 (55%) of these downstream cases had HIV co-infection and ten (50%) had died by the end of the study. 22 (56%) of 39 isolates in patients discharged home after treatment failure were resistant to eight or more drugs. However, five (16%) of 31 isolates were susceptible to rifabutin and more than 90% were likely to be sensitive to linezolid, bedaquiline, and delamanid. INTERPRETATION More than half of the patients with programmatically incurable tuberculosis were discharged into the community where they remained for an average of 16 months, were at risk of expectorating infectious cough aerosols, and posed a threat of transmission of extensively drug-resistant tuberculosis. Urgent action, including appropriate containment strategies, is needed to address this situation. Access to delamanid, bedaquiline, linezolid, and rifabutin, when appropriate, must be accelerated along with comprehensive drug susceptibility testing. FUNDING UK Medical Research Council, South African Medical Research Council, South African National Research Foundation, European & Developing Countries Clinical Trials Partnership, Oppenheimer Foundation, Newton Fund, Biotechnology and Biological Sciences Research Council, King Abdullah University of Science & Technology.

Genome-wide analysis of multi- and extensively drug-resistant Mycobacterium tuberculosis

To characterize the genetic determinants of resistance to antituberculosis drugs, we performed a genome-wide association study (GWAS) of 6,465 Mycobacterium tuberculosis clinical isolates from more than 30 countries. A GWAS approach within a mixed-regression framework was followed by a phylogenetics-based test for independent mutations. In addition to mutations in established and recently described resistance-associated genes, novel mutations were discovered for resistance to cycloserine, ethionamide and para-aminosalicylic acid. The capacity to detect mutations associated with resistance to ethionamide, pyrazinamide, capreomycin, cycloserine and para-aminosalicylic acid was enhanced by inclusion of insertions and deletions. Odds ratios for mutations within candidate genes were found to reflect levels of resistance. New epistatic relationships between candidate drug-resistance-associated genes were identified. Findings also suggest the involvement of efflux pumps (drrA and Rv2688c) in the emergence of resistance. This study will inform the design of new diagnostic tests and expedite the investigation of resistance and compensatory epistatic mechanisms.A GWAS of multi- and extensively drug-resistant tuberculosis using 6,465 Mycobacterium tuberculosis clinical isolates from more than 30 countries identifies novel mutations associated with resistance. The capacity to detect resistance in particular to ethionamide, pyrazinamide, capreomycin, cycloserine and paraaminosalicylic acid was enhanced by inclusion of insertions and deletions.

The variability and reproducibility of whole genome sequencing technology for detecting resistance to anti-tuberculous drugs

BackgroundThe emergence of resistance to anti-tuberculosis drugs is a serious and growing threat to public health. Next-generation sequencing is rapidly gaining traction as a diagnostic tool for investigating drug resistance in Mycobacterium tuberculosis to aid treatment decisions. However, there are few little data regarding the precision of such sequencing for assigning resistance profiles.MethodsWe investigated two sequencing platforms (Illumina MiSeq, Ion Torrent PGM™) and two rapid analytic pipelines (TBProfiler, Mykrobe predictor) using a well characterised reference strain (H37Rv) and clinical isolates from patients with tuberculosis resistant to up to 13 drugs. Results were compared to phenotypic drug susceptibility testing. To assess analytical robustness individual DNA samples were subjected to repeated sequencing.ResultsThe MiSeq and Ion PGM systems accurately predicted drug-resistance profiles and there was high reproducibility between biological and technical sample replicates. Estimated variant error rates were low (MiSeq 1 per 77 kbp, Ion PGM 1 per 41 kbp) and genomic coverage high (MiSeq 51-fold, Ion PGM 53-fold). MiSeq provided superior coverage in GC-rich regions, which translated into incremental detection of putative genotypic drug-specific resistance, including for resistance to para-aminosalicylic acid and pyrazinamide. The TBProfiler bioinformatics pipeline was concordant with reported phenotypic susceptibility for all drugs tested except pyrazinamide and para-aminosalicylic acid, with an overall concordance of 95.3%. When using the Mykrobe predictor concordance with phenotypic testing was 73.6%.ConclusionsWe have demonstrated high comparative reproducibility of two sequencing platforms, and high predictive ability of the TBProfiler mutation library and analytical pipeline, when profiling resistance to first- and second-line anti-tuberculosis drugs. However, platform-specific variability in coverage of some genome regions may have implications for predicting resistance to specific drugs. These findings may have implications for future clinical practice and thus deserve further scrutiny, set within larger studies and using updated mutation libraries.

Primary macrophages and J774 cells respond differently to infection with Mycobacterium tuberculosis.

Macrophages play an essential role in the early immune response to Mycobacterium tuberculosis and are the cell type preferentially infected in vivo. Primary macrophages and macrophage-like cell lines are commonly used as infection models, although the physiological relevance of cell lines, particularly for host-pathogen interaction studies, is debatable. Here we use high-throughput RNA-sequencing to analyse transcriptome dynamics of two macrophage models in response to M. tuberculosis infection. Specifically, we study the early response of bone marrow-derived mouse macrophages and cell line J774 to infection with live and γ-irradiated (killed) M. tuberculosis. We show that infection with live bacilli specifically alters the expression of host genes such as Rsad2, Ifit1/2/3 and Rig-I, whose potential roles in resistance to M. tuberculosis infection have not yet been investigated. In addition, the response of primary macrophages is faster and more intense than that of J774 cells in terms of number of differentially expressed genes and magnitude of induction/repression. Our results point to potentially novel processes leading to immune containment early during M. tuberculosis infection, and support the idea that important differences exist between primary macrophages and cell lines, which should be taken into account when choosing a macrophage model to study host-pathogen interactions.

The draft genome of Mycobacterium aurum, a potential model organism for investigating drugs against Mycobacterium tuberculosis and Mycobacterium leprae

Mycobacterium aurum (M. aurum) is an environmental mycobacteria that has previously been used in studies of anti-mycobacterial drugs due to its fast growth rate and low pathogenicity. The M. aurum genome has been sequenced and assembled into 46 contigs, with a total length of 6.02Mb containing 5684 annotated protein-coding genes. A phylogenetic analysis using whole genome alignments positioned M. aurum close to Mycobacterium vaccae and Mycobacterium vanbaalenii, within a clade related to fast-growing mycobacteria. Large-scale genomic rearrangements were identified by comparing the M. aurum genome to those of Mycobacterium tuberculosis and Mycobacterium leprae. M. aurum orthologous genes implicated in resistance to anti-tuberculosis drugs in M. tuberculosis were observed. The sequence identity at the DNA level varied from 68.6% for pncA (pyrazinamide drug-related) to 96.2% for rrs (streptomycin, capreomycin). We observed two homologous genes encoding the catalase-peroxidase enzyme (katG) that is associated with resistance to isoniazid. Similarly, two embB homologues were identified in the M. aurum genome. In addition to describing for the first time the genome of M. aurum, this work provides a resource to aid the use of M. aurum in studies to develop improved drugs for the pathogenic mycobacteria M. tuberculosis and M. leprae.

Genome-wide diversity and differentiation in New World populations of the human malaria parasite Plasmodium vivax

Background The Americas were the last continent colonized by humans carrying malaria parasites. Plasmodium falciparum from the New World shows very little genetic diversity and greater linkage disequilibrium, compared with its African counterparts, and is clearly subdivided into local, highly divergent populations. However, limited available data have revealed extensive genetic diversity in American populations of another major human malaria parasite, P. vivax. Methods We used an improved sample preparation strategy and next-generation sequencing to characterize 9 high-quality P. vivax genome sequences from northwestern Brazil. These new data were compared with publicly available sequences from recently sampled clinical P. vivax isolates from Brazil (BRA, total n = 11 sequences), Peru (PER, n = 23), Colombia (COL, n = 31), and Mexico (MEX, n = 19). Principal findings/Conclusions We found that New World populations of P. vivax are as diverse (nucleotide diversity π between 5.2 × 10−4 and 6.2 × 10−4) as P. vivax populations from Southeast Asia, where malaria transmission is substantially more intense. They display several non-synonymous nucleotide substitutions (some of them previously undescribed) in genes known or suspected to be involved in antimalarial drug resistance, such as dhfr, dhps, mdr1, mrp1, and mrp-2, but not in the chloroquine resistance transporter ortholog (crt-o) gene. Moreover, P. vivax in the Americas is much less geographically substructured than local P. falciparum populations, with relatively little between-population genome-wide differentiation (pairwise FST values ranging between 0.025 and 0.092). Finally, P. vivax populations show a rapid decline in linkage disequilibrium with increasing distance between pairs of polymorphic sites, consistent with very frequent outcrossing. We hypothesize that the high diversity of present-day P. vivax lineages in the Americas originated from successive migratory waves and subsequent admixture between parasite lineages from geographically diverse sites. Further genome-wide analyses are required to test the demographic scenario suggested by our data.

Genome-wide mutagenesis and multi-drug resistance in American trypanosomes induced by the front-line drug benznidazole.

Chagas disease is caused by the protozoan parasite Trypanosoma cruzi and affects 5–8 million people in Latin America. Although the nitroheterocyclic compound benznidazole has been the front-line drug for several decades, treatment failures are common. Benznidazole is a pro-drug and is bio-activated within the parasite by the mitochondrial nitroreductase TcNTR-1, leading to the generation of reactive metabolites that have trypanocidal activity. To better assess drug action and resistance, we sequenced the genomes of T. cruzi Y strain (35.5 Mb) and three benznidazole-resistant clones derived from a single drug-selected population. This revealed the genome-wide accumulation of mutations in the resistant parasites, in addition to variations in DNA copy-number. We observed mutations in DNA repair genes, linked with increased susceptibility to DNA alkylating and inter-strand cross-linking agents. Stop-codon-generating mutations in TcNTR-1 were associated with cross-resistance to other nitroheterocyclic drugs. Unexpectedly, the clones were also highly resistant to the ergosterol biosynthesis inhibitor posaconazole, a drug proposed for use against T. cruzi infections, in combination with benznidazole. Our findings therefore identify the highly mutagenic activity of benznidazole metabolites in T. cruzi, demonstrate that this can result in multi-drug resistance, and indicate that vigilance will be required if benznidazole is used in combination therapy.

Global analysis of Plasmodium falciparum histidine-rich protein-2 (pfhrp2) and pfhrp3 gene deletions using whole-genome sequencing data and meta-analysis

Many rapid diagnostic tests (RDT) used on suspected malaria cases are based on the detection of the protein encoded by the Plasmodium falciparum histidine-rich protein-2 (pfhrp2) gene, which shares a high sequence homology with pfhrp3 in the 3D7 reference genome. Parasite isolates showing pfhrp2 and pfhrp3 gene deletions have been emerging over the years, but a comprehensive genetic analysis of these variants is still lacking. With this purpose, genomic data from experimental P. falciparum genetic crosses between different laboratory lines (3D7, HB3, DD2, 7G8 and GB4) were first analysed (n = 98). The frequency of pfhrp2 deletions was consistent with a Mendelian prediction in HB3 × DD2 (56.7%; 95%CI = (39.5%-72.9%)). Moreover, the pfhrp2 and pfhrp3 deletions segregated independently of each other in the same genetic cross. Analysis of 3D7 × HB3 and 7G8 × GB4 estimated the probability of spontaneously generating a pfhrp2 deletion during sexual recombination to be up to 6.2%. Next, whole genome sequence data from 1970 P. falciparum isolates collected globally were analysed. Nine samples displayed depth of coverage consistent with pfhrp2 deletions (0.5%), but the corresponding split-read analysis could not confirm deletions in seven of these samples. Twenty-eight isolates had evidence of pfhrp3 deletions (1.4%), which are widespread in Southeast Asia. Finally, a meta-analysis of published data revealed a positive mean association between the frequencies of pfhrp2 and pfhrp3 deletions in Africa and South America. This result suggested a shared selective pressure acting on these genetic variants. In conclusion, evidence of genetic selection on both pfhrp2 and pfhrp3 deletions was presented, but experimental crosses do not provide evidence of a fitness cost of these variants. Further work is urgently needed to accurately determine the prevalence and the degree of association between these genetic variants, and the respective impact on diagnostic accuracy of many in-use RDT.

Methylation in Mycobacterium tuberculosis is lineage specific with associated mutations present globally.

DNA methylation is an epigenetic modification of the genome involved in regulating crucial cellular processes, including transcription and chromosome stability. Advances in PacBio sequencing technologies can be used to robustly reveal methylation sites. The methylome of the Mycobacterium tuberculosis complex is poorly understood but may be involved in virulence, hypoxic survival and the emergence of drug resistance. In the most extensive study to date, we characterise the methylome across the 4 major lineages of M. tuberculosis and 2 lineages of M. africanum, the leading causes of tuberculosis disease in humans. We reveal lineage-specific methylated motifs and strain-specific mutations that are abundant globally and likely to explain loss of function in the respective methyltransferases. Our work provides a set of sixteen new complete reference genomes for the Mycobacterium tuberculosis complex, including complete lineage 5 genomes. Insights into lineage-specific methylomes will further elucidate underlying biological mechanisms and other important phenotypes of the epi-genome.