Deciphering the growth behaviour of Mycobacterium africanum.

BACKGROUND: Human tuberculosis (TB) in West Africa is not only caused by M. tuberculosis but also by bacteria of the two lineages of M. africanum. For instance, in The Gambia, 40% of TB is due to infections with M. africanum West African 2. This bacterial lineage is associated with HIV infection, reduced ESAT-6 immunogenicity and slower progression to active disease. Although these characteristics suggest an attenuated phenotype of M. africanum, no underlying mechanism has been described. From the first descriptions of M. africanum in the literature in 1969, the time to a positive culture of M. africanum on solid medium was known to be longer than the time to a positive culture of M. tuberculosis. However, the delayed growth of M. africanum, which may correlate with the less virulent phenotype in the human host, has not previously been studied in detail. METHODOLOGY/PRINCIPAL FINDINGS: We compared the growth rates of M. tuberculosis and M. africanum isolates from The Gambia in two liquid culture systems. M. africanum grows significantly slower than M. tuberculosis, not only when grown directly from sputa, but also in growth experiments under defined laboratory conditions. We also sequenced four M. africanum isolates and compared their whole genomes with the published M. tuberculosis H37Rv genome. M. africanum strains have several non-synonymous SNPs or frameshift mutations in genes that were previously associated with growth-attenuation. M. africanum strains also have a higher mutation frequency in genes crucial for transport of sulphur, ions and lipids/fatty acids across the cell membrane into the bacterial cell. Surprisingly, 5 of 7 operons, recently described as essential for intracellular survival of H37Rv in the host macrophage, showed at least one non-synonymously mutated gene in M. africanum. CONCLUSIONS/SIGNIFICANCE: The altered growth behaviour of M. africanum might indicate a different survival strategy within host cells

Gene expression responses to anti-tuberculous drugs in a whole blood model.

BACKGROUND: There is a need for better tools to evaluate new or repurposed TB drugs. The whole blood bactericidal activity (WBA) assay has been advocated for this purpose. We investigated whether transcriptional responses in the WBA assay resemble TB responses in vivo, and whether the approach might additionally reveal mechanisms of action. RESULTS: 1422 of 1798 (79%) of differentially expressed genes in WBA incubated with the standard combination of rifampicin, isoniazid, pyrazinamide and ethambutol were also expressed in sputum (P < 0.0001) obtained from patients taking the same combination of drugs; these comprised well-established treatment-response genes. Gene expression profiles in WBA incubated with the standard drugs individually, or with moxifloxacin or faropenem (with amoxicillin and clavulanic acid) clustered by individual drug exposure. Distinct pathways were detected for individual drugs, although only with isoniazid did these relate to known mechanisms of drug action. CONCLUSIONS: Substantial agreement between whole blood cultures and sputum and the ability to differentiate individual drugs suggest that transcriptomics may add value to the whole blood assay for evaluating new TB drugs

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-genom

EthA/R-Independent Killing of Mycobacterium tuberculosis by Ethionamide.

Ethionamide (ETH) is part of the drug arsenal available to treat multi-drug resistant tuberculosis. The current paradigm of this pro-drug activation involves the mycobacterial enzyme EthA and the transcriptional repressor, EthR. However, several lines of evidence suggest the involvement of additional players. The ethA/R locus was deleted in Mycobacterium bovis BCG and three Mycobacterium tuberculosis (MTB) strains. While complete resistance to ETH was observed with BCG ethA/R KO, drug susceptibility and dose-dependent killing were retained in the ethA/R KO MTB mutants, suggesting the existence of an alternative pathway of ETH bio-activation in MTB. We further demonstrated that this alternative pathway is EthR-independent, whereby re-introduction of ethR in ethA/R KO MTB did not lead to increased resistance to ETH. Consistently, ethA KO MTB (with intact ethR expression) displayed similar ETH susceptibility profile as their ethA/R KO counterparts. To identify the alternative ETH bio-activator, spontaneous ETH-resistant mutants were obtained from ethA/R KO MTB and whole genome sequencing identified single nucleotide polymorphisms in mshA, involved in mycothiol biosynthesis and previously linked to ETH resistance. Deletion of mshA in ethA/R KO MTB led to complete ETH resistance, supporting that the role of MshA in ETH killing is EthA/R-independent. Furthermore mshA single KO MTB displayed levels of ETH resistance similar or greater than those obtained with ethA/R KO strains, supporting that mshA is as critical as ethA/R for ETH killing efficacy

Whole genome sequencing of drug resistant Mycobacterium tuberculosis isolates from a high burden tuberculosis region of North West Pakistan

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

Single-virion sequencing of lamivudine-treated HBV populations reveal population evolution dynamics and demographic history.

BACKGROUND: Viral populations are complex, dynamic, and fast evolving. The evolution of groups of closely related viruses in a competitive environment is termed quasispecies. To fully understand the role that quasispecies play in viral evolution, characterizing the trajectories of viral genotypes in an evolving population is the key. In particular, long-range haplotype information for thousands of individual viruses is critical; yet generating this information is non-trivial. Popular deep sequencing methods generate relatively short reads that do not preserve linkage information, while third generation sequencing methods have higher error rates that make detection of low frequency mutations a bioinformatics challenge. Here we applied BAsE-Seq, an Illumina-based single-virion sequencing technology, to eight samples from four chronic hepatitis B (CHB) patients - once before antiviral treatment and once after viral rebound due to resistance. RESULTS: With single-virion sequencing, we obtained 248-8796 single-virion sequences per sample, which allowed us to find evidence for both hard and soft selective sweeps. We were able to reconstruct population demographic history that was independently verified by clinically collected data. We further verified four of the samples independently through PacBio SMRT and Illumina Pooled deep sequencing. CONCLUSIONS: Overall, we showed that single-virion sequencing yields insight into viral evolution and population dynamics in an efficient and high throughput manner. We believe that single-virion sequencing is widely applicable to the study of viral evolution in the context of drug resistance and host adaptation, allows differentiation between soft or hard selective sweeps, and may be useful in the reconstruction of intra-host viral population demographic history

A reference genome and methylome for the Plasmodium knowlesi A1-H.1 line.

Plasmodium knowlesi, a common parasite of macaques, is recognised as a significant cause of human malaria in Malaysia. The P. knowlesi A1H1 line has been adapted to continuous culture in human erythrocytes, successfully providing an in vitro model to study the parasite. We have assembled a reference genome for the PkA1-H.1 line using PacBio long read combined with Illumina short read sequence data. Compared with the H-strain reference, the new reference has improved genome coverage and a novel description of methylation sites. The PkA1-H.1 reference will enhance the capabilities of the in vitro model to improve the understanding of P. knowlesi infection in humans

Methylation analysis of Klebsiella pneumoniae from Portuguese hospitals

Klebsiella pneumoniae is an important nosocomial infectious agent with a high antimicrobial resistance (AMR) burden. The application of long read sequencing technologies is providing insights into bacterial chromosomal and putative extra-chromosomal genetic elements (PEGEs) associated with AMR, but also epigenetic DNA methylation, which is thought to play a role in cleavage of foreign DNA and expression regulation. Here, we apply the PacBio sequencing platform to eight Portuguese hospital isolates, including one carbapenemase producing isolate, to identify methylation motifs. The resulting assembled chromosomes were between 5.2 and 5.5Mbp in length, and twenty-six PEGEs were found. Four of our eight samples carry blaCTX-M-15, a dominant Extended Spectrum Beta Lactamase in Europe. We identified methylation motifs that control Restriction-Modification systems, including GATC of the DNA adenine methylase (Dam), which methylates N6-methyladenine (m6A) across all our K. pneumoniae assemblies. There was a consistent lack of methylation by Dam of the GATC motif downstream of two genes: fosA, a locus associated with low level fosfomycin resistance, and tnpB transposase on IncFIB(K) plasmids. Overall, we have constructed eight high quality reference genomes of K. pneumoniae, with insights into horizontal gene transfer and methylation m6A motifs

Global genetic diversity of var2csa in Plasmodium falciparum with implications for malaria in pregnancy and vaccine development

Malaria infection during pregnancy, caused by the sequestering of Plasmodium falciparum parasites in the placenta, leads to high infant mortality and maternal morbidity. The parasite-placenta adherence mechanism is mediated by the VAR2CSA protein, a target for natural occurring immunity. Currently, vaccine development is based on its ID1-DBL2Xb domain however little is known about the global genetic diversity of the encoding var2csa gene, which could influence vaccine efficacy. In a comprehensive analysis of the var2csa gene in >2,000 P. falciparum field isolates across 23 countries, we found that var2csa is duplicated in high prevalence (>25%), African and Oceanian populations harbour a much higher diversity than other regions, and that insertions/deletions are abundant leading to an underestimation of the diversity of the locus. Further, ID1-DBL2Xb haplotypes associated with adverse birth outcomes are present globally, and African-specific haplotypes exist, which should be incorporated into vaccine design

A blood RNA transcript signature for TB exposure in household contacts.

BACKGROUND: Current tools for diagnosing latent TB infection (LTBI) detect immunological memory of past exposure but are unable to determine whether exposure is recent. We sought to identify a whole-blood transcriptome signature of recent TB exposure. METHODS: We studied household contacts of TB patients; healthy volunteers without recent history of TB exposure; and patients with active TB. We performed whole-blood RNA sequencing (in all), an interferon gamma release assay (IGRA; in contacts and healthy controls) and PET/MRI lung scans (in contacts only). We evaluated differentially-expressed genes in household contacts (log2 fold change ≥1 versus healthy controls; false-discovery rate  0.19). CONCLUSIONS: Transcriptomics can detect TB exposure and, with further development, may be an approach of value for epidemiological research and targeting public health interventions