Philip D. Butcher

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.

Nigeria has an emerging problem with multidrug-resistant tuberculosis (MDR-TB). Whole-genome sequencing was used to understand the epidemiology of tuberculosis and genetics of multi-drug resistance among patients from two tertiary referral centers in Southwest Nigeria. In line with previous molecular epidemiology studies, most isolates of Mycobacterium tuberculosis from this dataset belonged to the Cameroon clade within the Euro-American lineage. Phylogenetic analysis showed this clade was undergoing clonal expansion in this region, and suggests that it was involved in community transmission of sensitive and multidrug-resistant tuberculosis. Five patients enrolled for retreatment were infected with pre-extensively drug resistant (pre-XDR) due to fluoroquinolone resistance in isolates from the Cameroon clade. In all five cases resistance was conferred through a mutation in the gyrA gene. In some patients, genomic changes occurred in bacterial isolates during the course of treatment that potentially led to decreased drug susceptibility. We conclude that inter-patient transmission of resistant isolates, principally from the Cameroon clade, contributes to the spread of MDR-TB in this setting, underscoring the urgent need to curb the spread of multi-drug resistance in this region.

RIFAQUIN was a tuberculosis chemotherapy trial in southern Africa including regimens with high-dose rifapentine with moxifloxacin. Here, the application of whole-genome sequencing (WGS) is evaluated within RIFAQUIN for identifying new infections in treated patients as either relapses or reinfections. WGS is further compared with mycobacterial interspersed repetitive units-variable number tandem repeats (MIRU-VNTR) typing. This is the first report of WGS being used to evaluate new infections in a completed clinical trial for which all treatment and epidemiological data are available for analysis. DNA from 36 paired samples of Mycobacterium tuberculosis cultured from patients before and after treatment was typed using 24-loci MIRU-VNTR, in silico spoligotyping and WGS. Following WGS, the sequences were mapped against the reference strain H37Rv, the single-nucleotide polymorphism (SNP) differences between pairs were identified, and a phylogenetic reconstruction was performed. WGS indicated that 32 of the paired samples had a very low number of SNP differences (0–5; likely relapses). One pair had an intermediate number of SNP differences, and was likely the result of a mixed infection with a pre-treatment minor genotype that was highly related to the post-treatment genotype; this was reclassified as a relapse, in contrast to the MIRU-VNTR result. The remaining three pairs had very high SNP differences (>750; likely reinfections). WGS and MIRU-VNTR both similarly differentiated relapses and reinfections, but WGS provided significant extra information. The low proportion of reinfections seen suggests that in standard chemotherapy trials with up to 24 months of follow-up, typing the strains brings little benefit to an analysis of the trial outcome in terms of differentiating relapse and reinfection. However, there is a benefit to using WGS as compared to MIRU-VNTR in terms of the additional genotype information obtained, in particular for defining the presence of mixed infections and the potential to identify known and novel drug-resistance markers.

New treatment options are needed to maintain and improve therapy for tuberculosis, which caused the death of 1.5 million people in 2013 despite potential for an 86 % treatment success rate. A greater understanding of Mycobacterium tuberculosis (M.tb) bacilli that persist through drug therapy will aid drug development programs. Predictive biomarkers for treatment efficacy are also a research priority. Genome-wide transcriptional profiling was used to map the mRNA signatures of M.tb from the sputa of 15 patients before and 3, 7 and 14 days after the start of standard regimen drug treatment. The mRNA profiles of bacilli through the first 2 weeks of therapy reflected drug activity at 3 days with transcriptional signatures at days 7 and 14 consistent with reduced M.tb metabolic activity similar to the profile of pre-chemotherapy bacilli. These results suggest that a pre-existing drug-tolerant M.tb population dominates sputum before and after early drug treatment, and that the mRNA signature at day 3 marks the killing of a drug-sensitive sub-population of bacilli. Modelling patient indices of disease severity with bacterial gene expression patterns demonstrated that both microbiological and clinical parameters were reflected in the divergent M.tb responses and provided evidence that factors such as bacterial load and disease pathology influence the host-pathogen interplay and the phenotypic state of bacilli. Transcriptional signatures were also defined that predicted measures of early treatment success (rate of decline in bacterial load over 3 days, TB test positivity at 2 months, and bacterial load at 2 months). This study defines the transcriptional signature of M.tb bacilli that have been expectorated in sputum after two weeks of drug therapy, characterizing the phenotypic state of bacilli that persist through treatment. We demonstrate that variability in clinical manifestations of disease are detectable in bacterial sputa signatures, and that the changing M.tb mRNA profiles 0–2 weeks into chemotherapy predict the efficacy of treatment 6 weeks later. These observations advocate assaying dynamic bacterial phenotypes through drug therapy as biomarkers for treatment success.

Drug-resistant tuberculosis (TB) remains a major challenge to global health and to healthcare in the UK. In 2014, a total of 6,520 cases of TB were recorded in England, of which 1.4 % were multidrug-resistant TB (MDR-TB). Extensively drug-resistant TB (XDR-TB) occurs at a much lower rate, but the impact on the patient and hospital is severe. Current diagnostic methods such as drug susceptibility testing and targeted molecular tests are slow to return or examine only a limited number of target regions, respectively. Faster, more comprehensive diagnostics will enable earlier use of the most appropriate drug regimen, thus improving patient outcomes and reducing overall healthcare costs. Whole genome sequencing (WGS) has been shown to provide a rapid and comprehensive view of the genotype of the organism, and thus enable reliable prediction of the drug susceptibility phenotype within a clinically relevant timeframe. In addition, it provides the highest resolution when investigating transmission events in possible outbreak scenarios. However, robust software and database tools need to be developed for the full potential to be realized in this specialized area of medicine.

Increasing use of nucleic acid amplification tests (NAATs) as the primary means of diagnosing gonococcal infection has resulted in diminished availability of Neisseria gonorrhoeae antimicrobial susceptibility data. We conducted a prospective diagnostic assessment of a real-time PCR assay (NGSNP) enabling direct detection of gonococcal ciprofloxacin susceptibility from a range of clinical sample types. NGSNP, designed to discriminate an SNP associated with ciprofloxacin resistance within the N. gonorrhoeae genome, was validated using a characterized panel of geographically diverse isolates (n = 90) and evaluated to predict ciprofloxacin susceptibility directly on N. gonorrhoeae-positive NAAT lysates derived from genital (n = 174) and non-genital (n = 116) samples (n = 290), from 222 culture-confirmed clinical episodes of gonococcal infection. NGSNP correctly genotyped all phenotypically susceptible (n = 49) and resistant (n = 41) panel isolates. Ciprofloxacin-resistant N. gonorrhoeae was responsible for infection in 29.7% (n = 66) of clinical episodes evaluated. Compared with phenotypic susceptibility testing, NGSNP demonstrated sensitivity and specificity of 95.8% (95% CI 91.5%–98.3%) and 100% (95% CI 94.7%–100%), respectively, for detecting ciprofloxacin-susceptible N. gonorrhoeae, with a positive predictive value of 100% (95% CI 97.7%–100%). Applied to urogenital (n = 164), rectal (n = 40) and pharyngeal samples alone (n = 30), positive predictive values were 100% (95% CI 96.8%–100%), 100% (95% CI 87.2%–100%) and 100% (95% CI 82.4%–100%), respectively. Genotypic prediction of N. gonorrhoeae ciprofloxacin susceptibility directly from clinical samples was highly accurate and, in the absence of culture, will facilitate use of tailored therapy for gonococcal infection, sparing use of current empirical treatment regimens and enhancing acquisition of susceptibility data for surveillance.