Tuberculosis: Fighting an old disease with next-generation sequencing

  1. Anzaan Dippenaar
  2. Robin M Warren  Is a corresponding author
  1. Stellenbosch University, South Africa

Mycobacterium tuberculosis is the causative agent of tuberculosis, a disease that is a major threat to human health worldwide. It is estimated that approximately 9 million people were diagnosed with tuberculosis during 2013, and that 1.5 million died from the disease. The global tuberculosis epidemic is being driven by co-infection with HIV and by the emergence and spread of drug-resistant strains of M. tuberculosis. The World Health Organization has reported that 3.5% of new tuberculosis patients, and 20.5% of patients who had been treated before, had multidrug-resistant forms of the disease in 2013. The diagnosis and treatment of drug-resistant tuberculosis is clearly a major global health challenge (World Health Organization, 2014).

The enormity of the tuberculosis epidemic has created a desperate need to develop methods to monitor the dynamics of the disease. In the early 1990s, the discovery of repetitive elements in the genome of M. tuberculosis laid the foundation for the development of the science of molecular epidemiology (van Embden et al., 1993). These methods have shown that in situations where tuberculosis is common, epidemics are driven by the transmission of the bacteria between individuals. However, in low incidence settings, epidemics are driven by the ‘reactivation’ of bacteria that have been lying dormant in individuals since an earlier infection. It is also known that recurrent disease—when the symptoms reappear after a patient has apparently been cured—can occur through a second infection event, and that drug resistance is spread by transmission (Mathema et al., 2006).

Traditional methods to identify strains of M. tuberculosis rely on the analysis of small windows of the genome, and it has been assumed that the DNA sequences in these windows are variable enough to allow researchers to separate strains of M. tuberculosis that are evolutionarily close or distant. However, the true complexity of disease dynamics cannot be resolved by tracking strains using a small section of the genome. The development of next-generation sequencing platforms has made it possible to view the complete genetic information of the bacteria, which should improve the accuracy of efforts to monitor strains of M. tuberculosis as they move through space and time (Roetzer et al., 2013). Rapid whole genome sequencing promises to be the ultimate tool for epidemiological investigations, diagnosis, and for testing whether strains of bacteria are susceptible to particular drugs.

Now, in eLife, Judith Glynn of the London School of Hygiene and Tropical Medicine (LSHTM) and co-workers—with Guerra-Assunção as first author—report how a long-term large-scale whole genome sequencing strategy has been used to decipher the tuberculosis epidemic in a high prevalence setting with multiple sources of infection (Guerra-Assunção et al., 2015). They analysed the whole genome sequences of 1687 M. tuberculosis samples (isolates) collected from patients in the Karonga District of Malawi over a period of 15 years. This represents 72% of the total number of confirmed tuberculosis cases during that time. The various strains of M. tuberculosis can be grouped into seven ‘lineages’ that each contain bacteria descending from a common ancestor. Guerra-Assunção et al. found that the epidemic was largely driven by members of one lineage, which implies that either this lineage arrived in the area earlier than the others, or that the members of this lineage were more successful.

The genome of M. tuberculosis consists of ∼4.4 million bases and is generally believed to be relatively stable (Jagielski et al., 2014). To identify isolates that were directly related in a transmission network (i.e., recently transmitted from one patient to the next), Guerra-Assunção et al. used a cut-off point of up to ten differences in single nucleotide polymorphisms between the genomes of the isolates. Next, they developed a clustering formula to group together directly related isolates. Using this formula in combination with network-analysis (where isolates are linked according to genome sequence similarity), they found that strains from certain lineages were more likely to be transmitted between patients than others. This suggests that there are differences in the abilities of bacteria in the different lineages to cause disease. In this high-incidence setting, 66% of identified cases clustered together, of which 38% of the patients had evidence of recent infection, implying ongoing transmission of the bacteria. This indicates that reactivation of previous infection was the primary driving force behind this epidemic.

Glynn, Guerra-Assunção and co-workers—who are based at the LSHTM, the Karonga Prevention Study in Malawi and the Wellcome Trust Sanger Institute—also showed that the proportion of tuberculosis cases due to reactivation increased over the duration of the 15 year study, as demonstrated by a marked decrease in transmission between 1999–2001 (45%) and 2008–2010 (30%). Guerra-Assunção et al. suggest that this decrease is due to the implementation of antiretroviral therapy and isoniazid preventative therapy in Karonga. However, this is counter-intuitive because both treatments should protect against reactivation, thereby raising an important question as to how reactivation may work in this context. Significantly, this study shows that the tuberculosis control program in Karonga has reduced transmission of the bacteria. It also demonstrates that whole genome sequencing can provide new insights into tuberculosis epidemics, which could be used to advise and fine tune control programs.

Despite the advantages of whole genome sequencing, it is important to acknowledge the complexity of the technology and data analysis. This questions how useful it could be in high-incidence settings where tens of thousands of cases are diagnosed annually. Furthermore, the current technology is restricted to clinical isolates that need to undergo a lengthy culturing and DNA extraction process, which prevents its use as a real-time monitoring tool. Additionally, whole genome sequencing is labor intensive and financially demanding, although costs have decreased significantly over the last decade. Regardless of these challenges, this technology has the potential to immediately revolutionise drug susceptibility testing by identifying the complete repertoire of mutations in target genes that confer drug resistance (Steiner et al., 2014). Application of this technology would decrease diagnostic delay, thereby reducing transmission, morbidity and mortality and, at the same time, improving treatment outcome.

References

    1. van Embden JD
    2. Cave MD
    3. Crawford JT
    4. Dale JW
    5. Eisenach KD
    6. Gicquel B
    7. Hermans P
    8. Martin C
    9. McAdam R
    10. Shinnick TM
    (1993)
    Strain identification of Mycobacterium tuberculosis by DNA fingerprinting: recommendations for a standardized methodology
    Journal of Clinical Microbiology 31:406–409.

Article and author information

Author details

  1. Anzaan Dippenaar

    DST/NRF Centre of Excellence for Biomedical Tuberculosis Research, SAMRC Centre for TB Research and the Division of Molecular Biology and Human Genetics, Stellenbosch University, Stellenbosch, South Africa
    Competing interests
    No competing interests declared.
  2. Robin M Warren

    DST/NRF Centre of Excellence for Biomedical Tuberculosis Research, SAMRC Centre for TB Research and the Division of Molecular Biology and Human Genetics, Stellenbosch University, Stellenbosch, South Africa
    For correspondence
    rw1@sun.ac.za
    Competing interests
    No competing interests declared.

Publication history

  1. Version of Record published: March 3, 2015 (version 1)

Copyright

© 2015, Dippenaar and Warren

This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.

Metrics

  • 1,905
    Page views
  • 249
    Downloads
  • 2
    Citations

Article citation count generated by polling the highest count across the following sources: Scopus, Crossref, PubMed Central.

Download links

A two-part list of links to download the article, or parts of the article, in various formats.

Downloads (link to download the article as PDF)

Open citations (links to open the citations from this article in various online reference manager services)

Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)

  1. Anzaan Dippenaar
  2. Robin M Warren
(2015)
Tuberculosis: Fighting an old disease with next-generation sequencing
eLife 4:e06782.
https://doi.org/10.7554/eLife.06782

Further reading

  1. Whole genome sequencing is providing new insights into the spread of different lineages of tuberculosis.

    1. Ecology
    2. Epidemiology and Global Health
    Aleksandra Kovacevic, David RM Smith ... Lulla Opatowski
    Research Article

    Non-pharmaceutical interventions implemented to block SARS-CoV-2 transmission in early 2020 led to global reductions in the incidence of invasive pneumococcal disease (IPD). By contrast, most European countries reported an increase in antibiotic resistance among invasive Streptococcus pneumoniae isolates from 2019 to 2020, while an increasing number of studies reported stable pneumococcal carriage prevalence over the same period. To disentangle the impacts of the COVID-19 pandemic on pneumococcal epidemiology in the community setting, we propose a mathematical model formalizing simultaneous transmission of SARS-CoV-2 and antibiotic-sensitive and -resistant strains of S. pneumoniae. To test hypotheses underlying these trends five mechanisms were built into the model and examined: (1) a population-wide reduction of antibiotic prescriptions in the community, (2) lockdown effect on pneumococcal transmission, (3) a reduced risk of developing an IPD due to the absence of common respiratory viruses, (4) community azithromycin use in COVID-19 infected individuals, (5) and a longer carriage duration of antibiotic-resistant pneumococcal strains. Among 31 possible pandemic scenarios involving mechanisms individually or in combination, model simulations surprisingly identified only two scenarios that reproduced the reported trends in the general population. They included factors (1), (3), and (4). These scenarios replicated a nearly 50% reduction in annual IPD, and an increase in antibiotic resistance from 20% to 22%, all while maintaining a relatively stable pneumococcal carriage. Exploring further, higher SARS-CoV-2 R0 values and synergistic within-host virus-bacteria interaction mechanisms could have additionally contributed to the observed antibiotic resistance increase. Our work demonstrates the utility of the mathematical modeling approach in unraveling the complex effects of the COVID-19 pandemic responses on AMR dynamics.