Abstract
Genomic studies of Mycobacterium tuberculosis complex (MTBC) might shed light on the dynamics of its transmission, especially in high-burden settings, where recent outbreaks are embedded in the complex natural history of the disease. We applied Whole-genome sequencing (WGS) to characterize the local population of MTBC, unravel potential transmission links and evaluate associations with host and pathogen factors.
Methods A one-year prospective study was conducted in Mozambique, a high HIV/TB burden country. WGS was applied to 295 positive cultures. We combined phylogenetic, geographical and clustering analysis, and investigated associations between risk factors of transmission.
Findings A significant high proportion of strains were in recent transmission (45.5%). We fully characterized MTBC isolates by using phylogenetic approaches and dating evaluation. We found two likely endemic clades, comprised of 67 strains, belonging to L1.2, dating from the late XIX century and associated with recent spread among PLHIV.
Interpretation Our results unveil the population structure of MTBC in our setting. The clustering analysis revealed an unexpected pattern of spread and high rates of progression, suggesting the failure of control measures. The long-term presence of local strains in Mozambique, which were responsible for large transmission among HIV/TB coinfected patients, hint at possible coevolution with sympatric host populations and challenge the role of HIV in TB transmission.
Funding Ministry of Enterprise and Knowledge (Government of Catalonia & European Social Fund, AGAUR fellowship); European Research Council (ERC) European Union’s Horizon 2020.
Introduction
Tuberculosis (TB) remains the deadliest and one of the most prevalent infectious diseases worldwide.1 More than twenty years after the introduction of molecular tools in TB, it is now undoubted that unraveling the transmission dynamics of local epidemics is essential to tackle the ongoing spread.2 A complete understanding of who transmits, where, how and why, is essential for designing effective control interventions.3
The development of Whole-genome Sequencing (WGS) techniques is transforming the classical vision of the infection and decoding links unreachable by basic epidemiology or traditional genotyping.4 However, WGS is still scarcely applied in high-burden settings, where recent outbreaks are embedded in the complex natural evolution of the disease.5
On one hand, high-resolution genomic data allows for accurate cluster analysis to investigate outbreaks and decrypt transmission profiles. By using the SNPs cut-off approach, we can delineate transmission clusters with a resolution not seen before.4 On the other hand, the study of the origin and genetic diversity of local strains also helps to survey the evolutionary history of the Mycobacterium tuberculosis complex (MTBC) population and the potential variations between lineages.5,6
Phylogeographical methods have identified different lineages and sub-lineages of MTBC with differential distribution worldwide.7 This striking difference in geographic distribution has led to hypotheses as to why some are more widespread than others. Whereas Lineage (L) 4 is the most ubiquitous globally, L1 and L3 have been described as endemic to several high-burden regions in Asia and Eastern Africa.8 Although part of those differences can be explained by historical contingency9, host-pathogen coevolution is assumed to play a decisive role.10 In San Francisco, In San Francisco, there seemed to be a preference of transmission between lineages and hosts with the same geographic origin.11 Furthermore, some works suggest that this link between human populations and lineages from the same geographic origin breaks down with HIV status, suggesting a role for coevolution at immune system level.11 However, identifying the causes why a genotype succeeds locally remains a challenge and there is limited information on the interaction with other factors, such as HIV infection.
All the insights we are gaining from WGS data should be transferred to high-burden settings and used to define effective public health interventions. Mozambique is a country with one of the highest rates of TB and HIV/TB coinfection1, but information on MTBC population structure is scarce in the region.12,13 This is, to our knowledge, the first comprehensive population-based study in Mozambique on the application of WGS to characterize local strains, unravel likely transmission links, and evaluate possible associations with host and pathogen factors.
Methods
Study design and study population
This was a one-year prospective surveillance-based study (TOSSE study) implemented from August 2013 to 2014. It was conducted in the District of Manhiça, Maputo province, a semi-rural area in Southern Mozambique, 80 kilometers north of the capital, with an estimated population of 186,241 inhabitants at the time of the study (2013-2014)14, living in an area of 2,373 km.15 This is high TB/HIV burden area with a long history of high TB transmission.16,17
Presumptive TB adults, without history of previous TB treatment, presenting with TB compatible symptoms18 (no time criteria if HIV-positive), who attended any of the health units belonging to Manhiça District Hospital’s catchment area, were consecutively enrolled.19 Only participants with confirmed MTBC in culture were included in the analysis.
We consecutively enrolled participants who met the following criteria: presumptive TB-infected adults, without history of previous TB treatment, presenting TB-compatible symptoms (no time criteria in the case of people living with HIV), attending any of the health units belonging to Manhiça District Hospital’s catchment area.
Diagnostic procedures
Participants provided two sputum samples at time of diagnosis. Extra-pulmonary specimens were collected at hospital level. Diagnostic tests were performed at the Centro de Investigação em Saúde de Manhiça (CISM) - Biosafety level 3 (BSL3) laboratory, which is subject to external quality control and is ISO certified.
Ziehl-Neelsen (ZN)
Smear microscopy was done by Ziehl-Neelsen staining. Results were reported as negative or on a scale of positive grades according to international standards.20
Xpert MTB/RIF (Xpert)
Raw samples were tested according to the manufacturer instructions.21 Invalid results were excluded. The semiquantitative results for Xpert fell under the following categories: very low, low, medium, or high.
Solid and liquid culture
Samples with positive Xpert results were cultured. Remnant raw samples were decontaminated by Kubica method23 and resuspended. Afterwards, 500 microliters were inoculated into Mycobacteria Growth Indicator Tubes (MGIT) liquid medium and incubated in the Bactec MGIT 960 mycobacterial detection instrument (Becton Dickinson Microbiology System, BD, USA). Additionally, 200 microlitres were cultured in BD Lowenstein Jensen solid medium. After 42 days (for liquid culture) or 8 weeks (for solid culture) without growth, samples were classified as negative. In case of positive results, MTBC was confirmed using ZN staining and BD TB Identification test (Becton Dickinson Microbiology System, USA). Isolates were stored at -80ºC.
Sequencing library construction and bioinformatics pipeline
Culture isolates were shipped to the BSL3 laboratory of the Instituto de Biomedicina de Valencia (Spain). After inactivation, samples were used to prepare WGS libraries. Genomic libraries were constructed with the Nextera XT Sample Preparation kit (Illumina Inc., San Diego, CA) according to the manufacturer’s protocol24, with 12 cycles for indexing PCR. WGS was carried out in the MiSeq platform (2×300 cycles paired-end run; Illumina).
Sequence analysis was performed following a validated, previously-described bioinformatics pipeline.25 Briefly, FASTQ files were preprocessed with fastp26 in order to trim poor quality bases and potential sequencing errors. Later, in order to reduce likely contaminant reads, we classified and filtered those that did not belong to the MTBC using Kraken.27 Samples with less than 90% of MTBC reads were discarded for posterior analyses. After this filtering, reads were mapped against the MTBC most probable ancestral genome28 using the BWA-mem algorithm.29 Later, we discarded reads with ambiguous mapping based on the BWA MAPQ score (keep those with MAPQ=60) as well as potential duplicate reads by using picard tools. Samples with genomic coverage <90% were discarded. The variant calling was performed by a combination of SAMtools and VarScan.30
In order to avoid mapping errors and false SNPs, we kept variants that (i) were supported by at least 20 reads, (ii) were found in a frequency of at least 0.9, (iii) were not found inside detected indel areas, or (iv) were found in areas of high accumulation of variants (more than three variants in a 10-bp defined window). Variants were annotated using SnpEff. 31 Variants present in PE/PPE genes, phages, or repeated sequences were not considered. With these high-quality variants detected, we generated the alignment. Samples with at least two phylogenetic variants at frequency >10% were classified as mixed infection cases.
Phylogenetics and geographical analysis
Maximum-likelihood phylogeny was constructed with IQ-TREE32 under the General Time Reversible (GTR) model of evolution, with a bootstrap of 1000 replicates and using the -fconst option for accounting for invariant sites. Known drug-resistant positions were not considered for generating the phylogeny as they are highly homoplastic. Later, we aimed to define the population structure of closely related strains by using the fast hierarchical Bayesian Analysis of Population Structure (BAPS) algorithm implemented in R library fastbaps. BAPS groups were defined under the second level of clustering hierarchy.33
Geographic origin of the clades from Mozambique (MZ)
In order to unravel the potential geographic origin of the BAPS groups, we used the RASP program.34 This is based on both, Bayesian and parsimony approaches, and aims to estimate the ancestral geographic origin. It requires a phylogeny and the geographic origin of the tips as input. We reconstructed a phylogeny combining MZ isolates and 8,263 genomes representative of the MTBC global genetic diversity. We marked MZ strains as having their geographical origin in Mozambique, and strains from other datasets as having non-Mozambique origin. We coupled the Statistical-Dispersal Vicariance Analysis (S-DIVA) with a Bayesian Binary MCMC. For each analysis, we run 5 MCMC chains with 500000 cycles.
RASP output estimates the likely origin of ancestral nodes. We manually defined likely endemic clades as follow: i) they consisted in 1 or more BAPS, ii) the clade contained more than 10 samples from Mozambique and constituted >60% of the total clade, iii) the origin of the most recent common ancestor (MRCA) obtained with RASP was Mozambique (>80% of probability). Non-endemic clades were defined as those in which i) the clade contained more than 10 samples from Mozambique but constituted <60% of the total clade, and ii) the MRCA was not Mozambican. Unknown clades were those which did not meet any of the previous criteria.
Estimation of times
Dating analysis was performed by applying the Bayesian inference implemented in BEAST2 v2.6.535 A multifasta file was generated for the three major lineages (L1, L2 and L4). We used them as partitions in beauti. A GTR substitution model was defined with gamma count = 4 and empirical nucleotide frequencies. A strict clock model was selected; since we do not have tip or node calibration, we fixed different clock rates for each lineage (L1: 1.57 × 10−7; for L2: 4.1 × 10−8; L4: 3.79 × 10−8), according to Menardo et al, 2019.36 We evaluated exponential and constant coalescent population models, and selected the former, since we found evidence against the constant model (95% HPD of the rate growth not including 0). We chose an exponential prior for the population size (mean = 1, offset = 0), and for the exponential growth rate prior, we used the standard Laplace distribution (mu = 0.001, scale = 30.7). We corrected the xml files to specify the number of invariant sites as indicated here: https://groups.google.com/g/beast-users/c/QfBHMOqImFE. We ran two runs with 2 × 107 steps-long chains, sampling every 2000 steps and removing the initial 10% as burn-in. We evaluated that the mixing and the estimated sample size (ESS) of the posterior and of all parameters were larger than 200, with Tracer 1.7.1.37 We used logcombiner to combine the tree files from the independent runs and ggtree38 R package to annotate and plot the trees.
Lastly, by combining results from molecular dating and the likely ancestral geographical origin from RASP, we aimed to assess for how long previously defined clades have been circulating within the country. We calculated the median and 95%HPD (Highest Posterior Density interval) for MRCAs with origin in Mozambique for all clades, and we compared the average of the median years of introduction in the country among endemic and foreign clades.
Identification of recent transmission events
Transmission was evaluated by clustering analysis from pairwise distance. We employed the strict relatedness cut-off of 5 SNPs for describing recent transmission39 and we evaluated up to 10 thresholds for broad cluster definition. Clusters obtained from Mozambique (estimated incidence rate (IR): 551/100,000 in 2014) were compared to published transmission data (from population-based studies) from settings with different TB burden (Malawi, IR 2009 243/100,000; Valencia, Spain, IR 2015: 9.13/100,000). From available datasets, we chose annual data from 200940 and 2015 (manuscript under revision) respectively, based on the largest datasets close to 2013/4 and to have comparable timespan.
Statistical analysis
Statistical analyses were performed using the R statistical language (R version 3.5.2, The R Foundation for Statistical Computing Platform). For the phylogenetic analysis, branch lengths were extracted using the geiger41package. The fastbaps42 package was applied for the BAPS algorithm. For transmission analysis, transmission events were identified by using ape and adegenet packages. Tidyverse, tigerstats, naniar, epiR, purr, broom, ggplot and parameters,were used to visualize, describe and analyze epidemiological data.
We aimed to explore the association of being in transmission (dependent variable) with a range of collected risk factors (Table 1). The non-parametric Fisher’s exact test was used to identify differences in the distribution of independent variables. We hypothesized that transmission patterns would differ depending on the origin of clades (endemicity). Saturated logistic regression models were used to estimate mutually adjusted odds ratio (OR). Covariates with p-value<0.2, passed for the adjusted analysis. Age, sex and HIV status were chosen as a priori related risk factors irrespective of univariate associations. A backward strategy was used to finalize the model. Afterwards, considering the large prevalence of HIV-positive patients in our dataset, and that this unbalanced distribution might confound our findings, we stratified the analysis by HIV status (Table 2).
Ethics
The CISM’s Internal Scientific and Internal Bioethics Committees and the National Bioethics Committee (Maputo, Mozambique) waived ethical approval for this work (Reference: 199/CNBS/13). Written informed consent was provided by all study participants
Results
1. Overall population structure of MTBC strains in Manhiça district
From those 580 patients who started treatment during the study period in the district, 302 were microbiologically confirmed by culture, although 7 of those strains were not available for WGS due to poor quality of cultures. After exclusions (Figure 1), 275 strains were included in the analysis on population structure.
WGS: Whole genome sequencing
Overall, when reconstructing the phylogenetic tree (Figure 2), most cases (49.5%, 136/275) were classified as L4. The L4.3.4 (LAM) sub-lineage was the most common within L4 (23.3%, 64/275) although, at this level the most prevalent sublineage was L1.2 (25.4%, 70/275). The Beijing L2 were classified as L2.2.1, representing 12.7% (35/275) of the total population (Supplementary S1. Table 1). Additionally, eight samples were determined as mixed infections.
Lineages are represented in branch colours (pink for L1, blue for L2, purple for L3 and red for L4). Twenty-five BAPS groups are denoted by different colours in the outer ring of the phylogeny.
2. SNPs-threshold approach to evaluate transmission events (Figure 4)
The clustering analysis revealed that 57.8% of strains (159/275), were in transmission, applying a cut-off of 5 SNPs. To evaluate the ‘transmission profile’ we compared genetic distances up to 10 SNPs with previously published data from Malawi (2009) and Valencia (2015) (See methods section). The distribution of pairwise distances revealed that, in all three regions, a large proportion of cases were in transmission, although patterns differed from one population to another, towering in Mozambique (64.7%) compared to Malawi (41.2%), and Valencia (29.8%).
While TB transmission seemed continuous for all three sites, the burden of very recent spread was significantly higher in Mozambique: 45.5% (125/275) of strains shared identical genotype (0 SNPs genetic distance), whereas only 9.2% (11/119) and 15.9% (41/246) were identified for Malawi and Valencia, respectively.
3. Fine-scale population structure of MTBC
By applying the BAPS algorithm, we recognized 25 groups of strains that shared high genomic similarity at level 2 of clustering hierarchy: 9 from L1, 6 L2, 2 L3, 8 L4 (1 strain did not assemble) (Figure 2). These groups were mapped in a global phylogeny (Figure 3) and, through the use of RASP, we inferred the probable geographic origin of the most recent common ancestor (MRCA) for each group (see methods).
Footnotes: E: Endemic; For: Foreign (non-endemic); MRCA: Most recent common ancestror.
BAPS 1, 2, 3, 12, 13 and 14 were enriched in Mozambique isolates. BAPS 1 to 3, and 12 to 14 formed two largest groups sharing a MRCA with origin in the country (probability of 80%). Therefore, we merged those BAPS into two clades, and considered them ‘likely endemic’ or ‘local’ in our study setting (E1 and E2, 67/275, 24.4%) (Supplementary S2). Eight clades were defined as probably ‘non-endemic’ candidates (164 /275, 59.6%), and 11 did not accomplish defined criteria to be classified (see Methods). All strains categorized as endemic belonged to L1.2.
We then inferred the age of probable endemic and non-endemic clades (Supplementary S3). We found that, for those considered likely endemic, MRCAs located in Mozambique dated further back in time than in non-endemic ones. The mean of the median time when those local clades began to circulate in the country was 1894 [IQR 1887.7; 1900.6], whereas for non-endemic, it was estimated to be 2009 [IQR 1876.9; 2014.0] (Figure 4, Supplementary S4. Figure 4.1 & Figure 4.2).
A) Frequency of samples found to each pairwise distance. B) Cumulative frequency (%) of isolates in transmission from 0 to 10 SNPs pairwise distance. Colours: the blue line represents data from Mozambique(2013-2014); orange from Malawi (2009) and purple from Valencia (2015).
Pie charts represent the probability that the geographical origin of the ancestor was MZ. The 95% Highest Posterior Density (HPD) is drawn as orange intervals. The posterior indicates the probability distribution over the parameter state space. The red colour (posterior close to one) indicates the maximum probability.
Footnotes: E: Endemic; For: Foreign (non-endemic); 95% Highest Posterior Density
4. Impact of collected risk factor and endemicity on transmission
From those patients with WGS data, 223 strains had laboratory and epidemiological data available (Figure 1). These were used to investigate potential covariates associated with clustering, defined using a restrictive cut-off of pairwise distance at 5 SNPs to evaluate the burden of very recent transmission (see methods).
The median age of patients was 35 years [IQR 28.2;45.0], 35.4% (79/223) were female and the majority (72.2 %, 161/223) were people living with HIV (PLHIV). The median CD4+ count was 145.5 [49.9; 321.8] cells /mm3. The most common symptom was cough presented in 93.7% of cases (209/223). A sputum smear for acid-fast bacilli was positive in 65.0% (145/223) and 52.1% of participants had a pathological chest X-ray (CXR) (87/167 available). In this subset of data, 52.9% of isolates (118/223), were included in a cluster. Univariate exploration of other sociodemographic and microbiological variables is displayed in Table 1. The distribution of strains in recent transmission did not vary by sex, age, HIV status, symptoms, type of TB, different areas of the district or sublineage. Although none of the resistant strains were clustered, we were unable to explore this association further due to the limited number of isolates.
Based on the hypothesis that local genotypes might be linked to transmission, we explored the association between ‘endemicity’ and clustering. Overall, 24.2% (54/223) of strains were classified as ‘likely endemic’ (hereinafter endemic) and 28.0 % of them were in cluster (33/223). The results of the best fitted logistic regression model included are shown in Table 2. The odds of clustering did not differ by clades, whether endemic or not. Subsequent stratification by HIV status revealed that PLHIV who were infected with endemic strains had more than two times higher odds of being in recent transmission, regardless sex and age (2.11, 95%CI [1.04;4.45], p-value=0.04), and that this association disappeared among HIV-negative individuals.
Lastly, in light of those results, we investigated whether the degree of immunosuppression, measured through CD4 levels, might shape the distribution of endemic clades. We found that although the proportion of endemic strains increased as decreased the CD4 counts, this raise did not result meaningful (p-value=0.134) (Supplementary S5).
Discussion
High resolution WGS data on MTBC population structure and transmission patterns is scarce in the country13,43, and similar to the situation for most high-burden countries. In order to provide a comprehensive insight on the molecular epidemiology of MTBC in Southern Mozambique, we have conducted the first population-based study using WGS data. The clustering analysis unmasked a high recent transmission profile. We also found two large likely-endemic clades that were associated with recent transmission among PLHIV. These clades were dated back to the end of the XIX century.
When pairwise distances were compared to one-year data from Malawi and Spain, profiles revealed that the highest proportion of strains involved in transmission was in Mozambique (64.7%)), and was also higher than that found in other high-burden settings (Liberia (39%, cut-off of 12 SNPs)44 and reported in Malawi for 15 years (31%).40 Although high rates of ongoing transmission were expected3,45, almost half of the strains were connected by 0 SNPs. Assuming a molecular clock of 0.3-0.5 SNPs/year, this would translate into very recent transmission events46 and fast rates of progression to active tuberculosis. Importantly, this proportion might be even underestimated as one-year population analyses are not enough to reveal the extent of recent transmission.
Interrupting ongoing spread is a critical public health intervention to bring down the current TB burden. Presented data reveals the suboptimal control of the epidemic in our setting and reinforce the need to explore the homegrown scenario. Most of participants were severely ill at the time of diagnosis, which might support previous works describing low awareness of disease47,48 and/or the existence of barriers in accessing medical care in our setting.49 Besides, programmatic active-case finding activities have been described as inconsistent in Mozambique, failing to reach hidden cases and breaking the transmission net.50 In our cohort of microbiologically-confirmed cases (highly transmitters), 20% of participants had contact with a TB case, although they sought health care because they were symptomatic (presumptive), proving the lack of active screening activities to identify patients in early stages.
The globally spread L4 was the most prevalent in our region, in agreement with previous works10, although at sublineage level, L1.2 was the most frequent. This finding is aligned with the global distribution of L1, concentrated around the Indian Ocean and that has been supposed to spread from South-Asia as a result of old migration pathways.51 To define clades specific to Mozambique,i.e. ‘probably endemic’, we applied a detailed, fine-scale population structure approach that allowed us to define two major endemic clades, both belonging to the L1.2 sub-lineage. In accordance to the hypothesis that they are endemic to our study setting, dating estimates confirmed they would have been circulating in Mozambique for longer times than non-endemic ones.
The long-term presence of those local strains in Mozambique suggests possible coevolution with the sympatric host population, as seen elsewhere.52 Thus, we tested the likelihood that endemic clades, as well as other collected covariates, could be related to a higher probability of recent spread using 5 SNPs as a threshold.39 Although through non-stratified models we did not find any remarkable result, posterior stratification by HIV status revealed that those classified as probably endemic presented higher odds of being responsible for recent transmission among HIV positive patients.
The extent to which TB/ HIV coinfection influences the structure of the MTBC population and the efficacy of transmission is currently under debate.53,54 On one hand, studies evaluating MTBC infectiousness among PLHIV are widely heterogeneous.55,54 Historically it has been argued that HIV decreases transmission, due to EPTB, lower cavitary lesions, or lower sputum bacillary load, among others.53,56 Nevertheless, recent works state that HIV-seropositive patients are as infectious as seronegative when they present cavitary disease or have a positive bacilloscopy.54,57 Furthermore, among PLHIV, most cases of active tuberculosis are the result of new infections.58 Therefore, for our setting, the large proportion of strains in recent spread may be due to the fact that half of PLHIV participants were severely immunosuppressed (highly susceptible to develop disease), microbiologically confirmed, and just 5% of them had extrapulmonary disease.
On the other hand, studies on coevolution have also demonstrated increased transmission between sympatric linage-host associations11, although, differently to our results, this was disrupted by HIV coinfections. Yet, the interdependence of both epidemics in high burden settings needs to embrace socioeconomical determinants that may shape the circulation and spread of MTBC and confound results. 56 Further studies on how HIV may influence sympatric associations in high burden areas are needed. Alternatively, our local population, regardless HIV status, may represent a reservoir for endemic strains with lower transmissibility, with some data suggesting that this is the case for the L1 specialist genotype.59 Limitations to understand this interplay includes the fact that the ancestry of individuals involved in this study is not known and that transmission estimates are based on a one-year sampling window. Although this may have limited our result, given the fast progression of the disease in PLHIV and the particularity of our dataset (45.5% of strains were 0 SNPs pairwise distance), one-year data may indeed be give us a hint of what is happening in our region.
Lastly, we cannot forget that host-related factors must be considered when interpreting genomic data. Information on how host and pathogen genotypes interplay is still scarce and few genotype-genotype interaction studies are available.7 Since the adaptive immune response is an essential mechanism for host recognition and control of MTBC60, we hypothesize that the interaction of HIV, MTBC and human populations is more complex in countries where the two epidemics collide with downstream impact on TB transmission. We highlight the need for population-based genome-to-genome association studies, including the MTBC, human and HIV genotype combined with sociodemographic data that could potentially confound results.
Conclusion
Overall, our results unveil the population structure of MTBC in a high TB and HIV setting. We found an unexpected pattern of spread, with a substantial amount of strains of recent transmission, suggesting the uncontrolled TB spread and high rates of progression to active disease. We also identified endemic groups that were estimated to be circulating in the country for more than one century and that were responsible for recent transmission among PLHIV.
Data Availability
Raw sequencing data will be available via online repository at the European Nucleotide Archive (ENA) under study accession number PRJEB27421 and SAMEA12029268. De–identified dataset containing patient–level data produced in the present study are available upon reasonable request to the authors
Contributors
All authors contributed as investigators to this study and to drafting and revising the final manuscript. All of them had full access to all the data in the study and had final responsibility for the decision to submit for publication.
Declaration of interests
IC received consultancy fees from Foundation for innovative new diagnostics. The authors have no other competing interests to declare.
Data availability
Raw sequencing data will be available via online repository at the European Nucleotide Archive (ENA) under study accession number PRJEB27421and SAMEA12029268. A limited de-identified dataset containing patient-level data will also be made available on publication.
Supplementary material
Supplementary S1. Table 1. Distribution of lineages and sublineages overall and stratify by clustering and likely endemic clades
Supplementary S2. Table 2. Summary of BAPs groups
Supplementary S3. Table 3. Median time to the Mozambican Most Recent Common C (MRCA) for each clade and calculation of mean’s median times, by BAPS type (endemic (E) or foreign (For))
Supplementary S4. Figure 4.1 BEAST tree for strains belonging to lineage 3; Figure 4.2. BEAST tree lineage 4.
Supplementary S5. Figure S5. Distribution of endemic/non-endemic strains by HIV status, stratified by CD4 counts.
Acknowledgments
This project has received funding from the European Research Council (ERC) under the European Union’s Horizon 2020 research and innovation programs 101001038 (TB-RECONNECT), PID2019-104477RB-I00 from Ministerio de Economía y Competitividad (Spanish Government) (to I.C.).
We acknowledge support from the Spanish Ministry of Science, Innovation and Universities through the “Centro de Excelencia Severo Ochoa 2019-2023” Program (CEX2018-000806-S), and support from the Generalitat de Catalunya through the CERCA Program.
B.S receives a pre-doctoral fellowship from the Secretariat of Universities and Research, Ministry of Enterprise and Knowledge of the Government of Catalonia and co-funded by European Social Fund (AGAUR).
References
