Population-based sequencing of Mycobacterium tuberculosis reveals how current population dynamics are shaped by past epidemics

  1. Irving Cancino-Muñoz
  2. Mariana G López  Is a corresponding author
  3. Manuela Torres-Puente
  4. Luis M Villamayor
  5. Rafael Borrás
  6. María Borrás-Máñez
  7. Montserrat Bosque
  8. Juan J Camarena
  9. Caroline Colijn
  10. Ester Colomer-Roig
  11. Javier Colomina
  12. Isabel Escribano
  13. Oscar Esparcia-Rodríguez
  14. Francisco García-García
  15. Ana Gil-Brusola
  16. Concepción Gimeno
  17. Adelina Gimeno-Gascón
  18. Bárbara Gomila-Sard
  19. Damiana Gónzales-Granda
  20. Nieves Gonzalo-Jiménez
  21. María Remedios Guna-Serrano
  22. José Luis López-Hontangas
  23. Coral Martín-González
  24. Rosario Moreno-Muñoz
  25. David Navarro
  26. María Navarro
  27. Nieves Orta
  28. Elvira Pérez
  29. Josep Prat
  30. Juan Carlos Rodríguez
  31. Ma Montserrat Ruiz-García
  32. Hermelinda Vanaclocha
  33. Valencia Region Tuberculosis Working Group
  34. Iñaki Comas  Is a corresponding author
  1. Tuberculosis Genomics Unit, Instituto de Biomedicina de Valencia (IBV-CSIC), Spain
  2. Unidad Mixta "Infección y Salud Pública" (FISABIO-CSISP), Spain
  3. Microbiology Service, Hospital Clínico Universitario, Spain
  4. Microbiology and Parasitology Service, Hospital Universitario de La Ribera, Spain
  5. Microbiology Service, Hospital Arnau de Vilanova, Spain
  6. Microbiology Service, Hospital Universitario Dr Peset, Spain
  7. Department of Mathematics, Faculty of Science, Simon Fraser University, Canada
  8. Microbiology Laboratory, Hospital Virgen de los Lirios, Spain
  9. Microbiology Service, Hospital de Denia, Spain
  10. Computational Genomics Department, Centro de Investigación Príncipe Felipe, Spain
  11. Microbiology Service, Hospital Universitari i Politècnic La Fe, Spain
  12. Microbiology Service, Hospital General Universitario de Valencia, Spain
  13. Microbiology Service, Hospital General Universitario de Alicante, Spain
  14. Microbiology Service, Hospital General Universitario de Castellón, Spain
  15. Microbiology Service, Hospital Lluís Alcanyis, Spain
  16. Microbiology Service, Hospital General Universitario de Elche, Spain
  17. Microbiology Service, Hospital Universitario de San Juan de Alicante, Spain
  18. Microbiology Service, Hospital de la Vega Baixa, Spain
  19. Subdirección General de Epidemiología y Vigilancia de la Salud y Sanidad Ambiental de Valencia (DGSP), Spain
  20. Microbiology Service, Hospital de Sagunto, Spain
  21. CIBER of Epidemiology and Public Health (CIBERESP), Spain
10 figures, 2 tables and 8 additional files

Figures

Genomic characterization of the study region.

(A) Phylogeny of 775 tuberculosis (TB) isolates collected during the years 2014 and 2016. Each ring represents genomic clusters detected by different single nucleotide polymorphism (SNP) thresholds …

Figure 1—source data 1

Genomic cluster types; Spanish: clusters including only Spanish-born cases; foreign: clusters including only foreign-born cases; mix: clusters including Spanish and foreign-born cases.

Cluster ID; number of Spanish. Foreign and unknown origin cases and total cluster size are indicated.

https://cdn.elifesciences.org/articles/76605/elife-76605-fig1-data1-v1.xlsx
Comparison between epidemiological and genomic clustering.

(A) Clustered samples using different pairwise distance thresholds, bars denote the number of cases within clusters for each single nucleotide polymorphism (SNP) threshold. Gray dashed line …

Figure 3 with 6 supplements
Historical transmission dynamics analysis.

(A) Distribution of local-born cases clustered by different pairwise distance SNP thresholds. Cases are expressed as the percentage of the plotted samples. Pie charts represent the proportion of …

Figure 3—source data 1

Bayesian dating results for all clusters (CL) from Oxfordshire.

Data include information of cluster date (in years AD); cluster distances (min, max, mean, and median); most recent common ancestor (MRCA) time in years (AD) with corresponding 95% highest probability density (HPD) intervals and cluster time span or duration.

https://cdn.elifesciences.org/articles/76605/elife-76605-fig3-data1-v1.xlsx
Figure 3—source data 2

Number of transmission links (TLs) traced back to 150 years before 2016 for Oxfordshire dataset.

Median time of all the local TLs pointed in Figure 3—figure supplement 1 were classified in different time windows. The total number of links within each time window is indicated even if they were links within a genomic cluster (total_CL) or not (total_noCL). Sampling period for Oxfordshire is 2007–2012.

https://cdn.elifesciences.org/articles/76605/elife-76605-fig3-data2-v1.xlsx
Figure 3—source data 3

Bayesian dating results for all clusters (CL) from Valencia region.

Data include information of cluster date (in years AD); cluster distances (min, max, mean, and median); most recent common ancestor (MRCA) time in years (AD) with corresponding 95% highest probability density (HPD) intervals and cluster time span or duration.

https://cdn.elifesciences.org/articles/76605/elife-76605-fig3-data3-v1.xlsx
Figure 3—source data 4

Number of transmission links (TLs) traced back to 150 years before 2016 for Oxfordshire dataset.

Median time of all the local TLs pointed in Figure 3—figure supplement 1 were classified in different time windows. The total number of links within each time window is indicated even if they were links within a genomic cluster (total_CL) or not (total_noCL). Sampling period for Oxfordshire is 2007–2012.

https://cdn.elifesciences.org/articles/76605/elife-76605-fig3-data4-v1.xlsx
Figure 3—source data 5

Number of transmission links (TLs) traced back to 150 years before 2016 for Malawi dataset.

Analysis performed for Malawi. Median time of all the TLs pointed in Figure 3—figure supplement 2 were classified in different time windows. The total number of links within each time window is indicated even if they were links within a genomic cluster (total_CL) or not (total_noCL). Color code for links outside clusters (no_CL) are the same as in Figure 3—figure supplement 2. Sampling period for Malawi is 2008–2010.

https://cdn.elifesciences.org/articles/76605/elife-76605-fig3-data5-v1.xlsx
Figure 3—source data 6

Number of transmission links (TLs) traced up to 150 years before 2016 for Valencia dataset.

Median time of all the TLs pointed in Figure 3—figure supplements 36 were classified in different time windows. The total number of links within each time window is indicated even if they were links within a genomic cluster (total_CL) or not (total_noCL). Color code for links outside clusters (no_CL) are the same as in Figure 3—figure supplements 36. Sampling period for Valencia is 20014–2016. The percentage of Spanish cases within each TL was evaluated in a global phylogeny as a proxi of the origin of the node.

https://cdn.elifesciences.org/articles/76605/elife-76605-fig3-data6-v1.xlsx
Figure 3—figure supplement 1
Bayesian time tree for Oxfordshire dataset.

Sample names in bold represent local born cases. Local clusters (CLs) are highlighted in orange and mixed CLs in green (those including foreign cases). The CL numbers are crossreferenced in Figure …

Figure 3—figure supplement 2
Bayesian time tree for Malawi dataset.

Sample names in bold represent local born cases. Local clusters (CLs) are highlighted in orange and mixed CLs in green (those including foreign cases). The CL numbers are crossreferenced in Figure …

Figure 3—figure supplement 3
Bayesian time tree for Valencia region dataset.

Sample names in bold represent local born cases. Local clusters (CL) are highlighted in orange and mixed CL in green (those including foreign cases). The CL numbers are crossreferenced in Figure …

Figure 3—figure supplement 4
Bayesian time tree for Velncia Region dataset.

Sample names in bold represent local born cases. Local clusters are highlighted in orange and mixed clusters in green (those including foreign cases). CL numbers are crossreferencied in Figure …

Figure 3—figure supplement 5
Bayesian time tree for Valencia Region dataset.

Sample names in bold represent local born cases. Local clusters are highlighted in orange and mixed clusters in green (those including foreign cases). CL numbers are crossreferencied in Figure …

Figure 3—figure supplement 6
Bayesian time tree for Valencia Region dataset.

Sample names in bold represent local born cases. Local clusters are highlighted in orange and mixed clusters in green (those including foreign cases). CL numbers are crossreferencied in Figure …

Hypothetical time trees indicating transmission links (TLs).

(A) (Left) The complete phylogeny, including all bacterial isolates and displaying multiple transmission events over time (located at nodes for simplification). This scenario allows the …

Appendix 1—figure 1
Workflow for sample selection.

MTBC, Mycobacterium tuberculosis complex; WGS, whole genome sequencing.

Appendix 1—figure 2
Hypothetical time tree.

(A) Complete phylogeny, including all bacterial isolates. Dashed lines represent missing samples that could not be retrieved during the sampling period (gray dashed lines). Transmission events are …

Appendix 1—figure 3
Boxplots of ages of cases from spanish-born genomic clusters in the Valencia region vs. the inferred ancestor of the clusters.

Bars represent the percentage of cases in gClusters (by 12 SNPs) for each time period. Boxplots represent the age distribution of patients within the clusters. Differences between the age cases for …

Appendix 1—figure 4
Correlation between the height median estimated by two ascertainment bias correction approaches; ‘adjusting clock rate’ and ‘including invariant positions’.

Correlation was calculated for each dataset.

Author response image 1
Author response image 2

Tables

Table 1
Dating of local genomic clusters (gCluster).

Times of the oldest and youngest local gClusters obtained by a Bayesian analysis are presented, with values in years (AD) and 95% highest posterior density given in brackets. The number of gClusters …

DatasetSampling periodLocal samplesN local gClusterLocal clusteringMedian distance rangeOldest gClusterYoungest gCluster
Oxfordshire2006–201274627%0–71993 (1982–2003)2009 (2003–2012)
Malawi2008–20101064049.80%0–141979 (1968–1988)2009 (2004–2010)
Valencia region2014–20164566547.40%0–111985 (1972–1996)2015 (2012–2016)
Author response table 1
Total cases expected: 3 cases, cases observed: 5*based on (Chaoui et al., 2014; Coscolla et al., 2021; de Jong et al., 2010; Gehre et al., 2016; Yeboah-Manu et al., 2017).
Origen#total_patients%_reported*cases_expected
MAURITANIA1no data
SENEGAL7201.4
GUINEA EQ6no data
CAMEROON1120.12
MALI4200.8
NIGERIA2110.2
MOROCCO47no Africanum reported
IVORY1220.2
ALGERIA5no data

Additional files

Supplementary file 1

Samples included in this study and all clinical.

Laboratory and epidemiological data associated. WGS data includes all the information obtained by WGS; run accession number are included in the ENA projects PRJEB29604 and PRJEB38719; depth coverage. MTBC lineage. genomic and epidemiological clusters ID. Patient information including: age, gender, country of birth, year of arrival to Spain, and province of residence in the Valencia region. NA indicates not available epidemiological information.

https://cdn.elifesciences.org/articles/76605/elife-76605-supp1-v1.xlsx
Supplementary file 2

Beast input file (.xml) indicating ascertainment bias correction implemented.

https://cdn.elifesciences.org/articles/76605/elife-76605-supp2-v1.zip
Supplementary file 3

Demographic.

clinical and epidemiological data for all culture positive cases and sequenced samples reported between 2014 and 2016 in Valencia region.

https://cdn.elifesciences.org/articles/76605/elife-76605-supp3-v1.xlsx
Supplementary file 4

Statistical analysis of clustered vs. unique TB cases.

Age, sex, place of birth, sputum smear, and disease type were considered in the analysis for all TB cases and for Spanish-born only.

https://cdn.elifesciences.org/articles/76605/elife-76605-supp4-v1.xlsx
Supplementary file 5

Epidemiological cluster evaluation.

Clusters identified by epidemiological links were compared against genomic clustering. Similarities and discrepancies between both methods are described for all samples. Cluster size and pairwise distance among the samples included in the cluster are mentioned. Median time and 95% HPD (highest posterior density) of the most common ancestor (tMRCA) for the cluster obtained by both approaches are provided. For mean and median values. Only the clusters with tMRCA were considered.

https://cdn.elifesciences.org/articles/76605/elife-76605-supp5-v1.xlsx
Supplementary file 6

Performance of the WGS method for transmission.

Values were extracted from 724 clinical TB isolates and assumed epidemiological links from contact tracing as the gold standard. NPV, negative predictive value; PPV, positive predictive value.

https://cdn.elifesciences.org/articles/76605/elife-76605-supp6-v1.xlsx
Supplementary file 7

Patients mean age for those cases within genomic clusters which ancestors’ age are included in different time windows in Valencia region.

https://cdn.elifesciences.org/articles/76605/elife-76605-supp7-v1.xlsx
MDAR checklist
https://cdn.elifesciences.org/articles/76605/elife-76605-mdarchecklist1-v1.pdf

Download links