Skip to main content
medRxiv
  • Home
  • About
  • Submit
  • ALERTS / RSS
Advanced Search

Genomic evidence for prisons as amplifiers of community tuberculosis epidemics

View ORCID ProfileKatharine S. Walter, Paulo César Pereira dos Santos, Thais Oliveira Gonçalves, Brunna Oliveira da Silva, Andrea da Silva Santos, Alessandra de Cássia Leite, Alessandra Moura da Silva, Flora Martinez Figueira Moreira, Roberto Dias de Oliveira, Everton Ferreira Lemos, Eunice Cunha, Yiran E. Liu, Albert Ko, Caroline Colijn, Ted Cohen, Barun Mathema, Julio Croda, View ORCID ProfileJason R. Andrews
doi: https://doi.org/10.1101/2021.06.22.21259360
Katharine S. Walter
1Division of Infectious Diseases and Geographic Medicine, Stanford University School of Medicine; Stanford, United States
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Katharine S. Walter
  • For correspondence: kwalter@stanford.edu
Paulo César Pereira dos Santos
2Health Sciences Research Laboratory, Federal University of Grande Dourados; Dourados, Brazil
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Thais Oliveira Gonçalves
3Laboratory of Bacteriology, Central Laboratory of Mato Grosso do Sul; Campo Grande, Brazil
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Brunna Oliveira da Silva
2Health Sciences Research Laboratory, Federal University of Grande Dourados; Dourados, Brazil
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Andrea da Silva Santos
2Health Sciences Research Laboratory, Federal University of Grande Dourados; Dourados, Brazil
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Alessandra de Cássia Leite
2Health Sciences Research Laboratory, Federal University of Grande Dourados; Dourados, Brazil
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Alessandra Moura da Silva
4School of Medicine, Federal University of Mato Grosso do Sul, School of Medicine; Campo Grande, Brazil
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Flora Martinez Figueira Moreira
2Health Sciences Research Laboratory, Federal University of Grande Dourados; Dourados, Brazil
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Roberto Dias de Oliveira
5Nursing Course, State University of Mato Grosso do Sul; Dourados, Brazil
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Everton Ferreira Lemos
4School of Medicine, Federal University of Mato Grosso do Sul, School of Medicine; Campo Grande, Brazil
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Eunice Cunha
3Laboratory of Bacteriology, Central Laboratory of Mato Grosso do Sul; Campo Grande, Brazil
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Yiran E. Liu
1Division of Infectious Diseases and Geographic Medicine, Stanford University School of Medicine; Stanford, United States
6Cancer Biology Graduate Program, Stanford University School of Medicine; Stanford, United States
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Albert Ko
7Department of Epidemiology of Microbial Diseases, Yale School of Public Health; New Haven, United States
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Caroline Colijn
8Department of Mathematics, Simon Fraser University; Burnaby, Canada
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Ted Cohen
7Department of Epidemiology of Microbial Diseases, Yale School of Public Health; New Haven, United States
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Barun Mathema
9Department of Epidemiology, Columbia University Mailman School of Public Health; New York, United States
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Julio Croda
4School of Medicine, Federal University of Mato Grosso do Sul, School of Medicine; Campo Grande, Brazil
7Department of Epidemiology of Microbial Diseases, Yale School of Public Health; New Haven, United States
10Mato Grosso do Sul Office, Oswaldo Cruz Foundation; Campo Grande, Brazil
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Jason R. Andrews
1Division of Infectious Diseases and Geographic Medicine, Stanford University School of Medicine; Stanford, United States
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Jason R. Andrews
  • Abstract
  • Full Text
  • Info/History
  • Metrics
  • Supplementary material
  • Data/Code
  • Preview PDF
Loading

Abstract

Prisons are high-incidence settings for tuberculosis around the world, yet the contribution of spillover from prisons in driving community epidemics has not been quantified. We whole genome sequenced 1,152 M. tuberculosis isolates from participants diagnosed with tuberculosis within prisons and in the community in Mato Grosso do Sul, Brazil from 2014 to 2019. By integrating timed phylogenies and detailed location data, we reconstructed probabilistic transmission histories. M. tuberculosis sequences from incarcerated and non-incarcerated people were closely phylogenetically related. We found that 57% of recent community-wide tuberculosis cases were attributable to transmission from individuals with an incarceration history, 2.6% of the population. Further, we find genomic evidence that the prison system disseminates M. tuberculosis genotypes through frequent transfers across the state. This population-wide genomic transmission reconstruction framework can be applied to identify key environments amplifying infectious disease transmission to prioritize public health interventions.

Teaser Prison spillover plays an outsized role in community tuberculosis transmission.

Introduction

Tuberculosis remains one of the leading causes of death by an infectious disease and led to 1.4 million deaths in 2019.(1) Locating where M. tuberculosis transmission occurs is critical for public health interventions, yet in high-incidence settings, where most tuberculosis is transmitted remains unknown.(1–6) Prisons, hospitals, or mines are known to be high-risk environments(7–9) and may have epidemiological impacts that transcend their walls, serving as institutional amplifiers of infection.(10, 11) Identifying institutional amplifiers of tuberculosis is of particular concern in Central and South America, where the incarcerated population has increased by over 200% since 2000 and 11% percent of tuberculosis cases now occur among incarcerated people, who comprise <1% of the total population.(12)

The full impact of increasing incarceration rates on the tuberculosis epidemic in the Americas and elsewhere has been difficult to quantify. Many people are infected in prison, but only diagnosed later.(11, 13, 14) Further, spillover may occur when incarcerated people, prison staff, or visitors are infected within high-incidence prison environments and transmit M. tuberculosis onwards in communities outside prison. Previous studies have found that the excess risk of tuberculosis within prisons extends to surrounding communities(15, 16) and that M. tuberculosis genotypes responsible for jail and prison outbreaks are also found in the surrounding communities.(13, 15, 17, 18)

Pathogen genomes can be powerfully harnessed to infer high-resolution transmission histories, including who-infected-whom.(19) However, studies have not yet combined M. tuberculosis genomic and epidemiologic data to reconstruct transmission linkages between potential institutional amplifiers and the community, nor have they estimated the fraction of community transmission attributable to such spillover. To quantify the role of prisons in population-wide tuberculosis transmission, we conducted prospective genomic surveillance for tuberculosis in the Brazilian state of Mato Grosso do Sul and used full M. tuberculosis genomes, combined with detailed meta-data on incarceration history, to reconstruct tuberculosis transmission chains spanning prisons and the community.

Results

Tuberculosis and incarceration trends in Mato Grosso do Sul1

Brazil’s incarcerated population grew by 107% (from 361,402 to 748,009) from 2005 to 2019, an increase closely paralleled in the Central West state of Mato Grosso do Sul, where the incarcerated population has grown by 115% (7,891 to 16,976) over the same period (Fig. 1a-c). In 2019, prisons in the state were at 203% occupancy. While the state’s notification rate of new and retreatment tuberculosis in the general population (38 per 100,000) is similar to Brazil’s national tuberculosis notification rate (40 per 100,000) in 2019,(20) the notification rate within prisons was 43 times as high (1,666 per 100,000), again similar to Brazil’s national notification rate within prisons (1,596 per 100,000) in 2019 (Fig. 1d,e).

Figure 1.
  • Download figure
  • Open in new tab
Figure 1. Tuberculosis is increasingly concentrated within prisons in Mato Grosso do Sul, Brazil.

(a) Map of Brazil with states colored by the 2019 incarceration rate with Mato Grosso do Sul outlined in black. (b) Map of Mato Grosso do Sul state with the cities with passive surveillance only and where active surveillance was additionally conducted within prisons. (c) The incarcerated population size in Mato Grosso do Sul grew by 142% from 2005 to 2020, increasing from 7,891 to 19,065 people. (d) Mato Grosso do Sul state’s annual new and retreatment tuberculosis notifications, colored by incarceration status and (e) the notification rate per 100,000 people for the incarcerated population and non-incarcerated populations from 2009 to 2019. The y-axis is log-scaled. The mean tuberculosis notification rate was more than 46 times greater among the incarcerated population compared to the non-incarcerated population. The green and yellow bars (d) and (e) indicate the period of active surveillance in prisons in Campo Grande and Dourados respectively.

Population-based genomic surveillance across Mato Grosso do Sul state

From 2014 to 2019, 3,491 new and retreatment cases of tuberculosis were reported in the Campo Grande and Dourados, the two largest cities in Mato Grosso do Sul, Brazil. Of these, 1,249 had positive cultures and we sequenced 787 M. tuberculosis isolates. We sequenced additional isolates from three other cities in the state and from earlier years in Dourados for a total of 1,090 isolates from unique tuberculosis notifications (Fig. 1b). Whole genome sequences (WGS) for 1,043 isolates met our coverage and quality criteria (Methods). We excluded 10% (108/1043) of isolates with mixed infection, resulting in 935 high-quality genomes from 935 unique tuberculosis episodes from 918 individuals.

Of the 935 isolates in our analyses, 50% (465/935) were incarcerated; 16% (150/935) were from patients who were formerly incarcerated; and 34% (320/935) of the study population did not have an incarceration history at the time of tuberculosis notification. Among those who did not have an incarceration history, 32 people reported contact with incarcerated individuals. Additional population characteristics are in Table S1. Isolates were largely from M. tuberculosis lineage 4 and predominantly fell into three sublineages: lineage 4.1, 4.3, and 4.4 (Fig. 2). Overall, we found a low prevalence of antibiotic resistance across isolates, with 93.7% of isolates susceptible to all antibiotics (Fig. 2).

Figure 2.
  • Download figure
  • Open in new tab
Figure 2. M. tuberculosis isolates from incarcerated and non-incarcerated people are closely related in Mato Grosso do Sul, Brazil.

A maximum likelihood phylogeny of 932 tuberculosis isolates from Lineage 4 inferred from a multiple sequence alignment of 19,753 SNPs. Branch lengths are in units of substitutions per site. From the inside, rings are colored by incarceration status, antimicrobial resistance prediction, and genomic cluster (Methods). “Other” resistance indicates resistance to at least one antibiotic other than isoniazid or rifampicin. Branches are colored by sublineage. The clade containing the largest predicted genomic cluster (Cluster 5), including 170 isolates, is highlighted in purple.

Phylogenetic structure of M. tuberculosis from prisons and the community

A maximum likelihood phylogeny (Fig. 2) demonstrates that M. tuberculosis isolates sampled from incarcerated and non-incarcerated people do not form distinct clades and are closely phylogenetically related. Terminal branch lengths, a proxy for the degree of recent transmission,(21) are significantly longer for isolates from people with no incarceration history at the time of tuberculosis notification (9.4 × 10−4 substitutions/site) than isolates from incarcerated individuals (2.2 × 10−4; two-sample t-test p-value = 1.14 × 10−7) or formerly incarcerated individuals (3.8 × 10−4; p = 1.6 × 10−4), suggesting that individuals with an incarceration history were more likely to be recently infected (Fig. S1).

Genetic clustering of M. tuberculosis from prisons and the community

To identify potential clusters of recent transmission, we applied a commonly used 12-SNP threshold,(22) including all isolates within the threshold distance of at least one other clustered isolate. Eighty-three percent (777/935) of the isolates fell into 84 genomic clusters (each including 2 to 170 isolates), providing evidence that notifications were largely attributable to recent, local transmission rather than travel-associated importation or re-activation of genetically distinct latent infections (Fig. 2). Isolates from incarcerated people were more frequently clustered (93.3%, 434/465), than those from formerly incarcerated (86.0%, 129/150; p <0.0001), or never incarcerated people (66.9%, 214/320; p <0.0001), again suggesting more recent transmission within prisons.

We predicted that if prison and community-associated epidemics were distinct, isolates from the community would be most closely related to and cluster with other isolates from the community and vice versa. Of the 46 clusters with three or more isolates, including 703 participants, 40 included isolates from community members who had no reported incarceration history. Of those potential transmission clusters, 85.0% (34/40) included isolates from people who were currently or previously incarcerated.

We found a similar pattern of clustering of isolates from community members and individuals with an incarceration history when using an alternative 5-SNP threshold for genomic clustering(23). With this threshold, of the 80 clusters with three or more isolates, 56 included isolates from people with no incarceration history and 78.6% (44/56) also included isolates from people with an incarceration history.

A M. tuberculosis clone spans prisons and the community across the state

The largest potential transmission cluster, including 170 isolates sampled from September 2010 to April 2019, was distributed across the state, including cases found across state’s two largest cities, Campo Grande and Dourados, as well as the smaller cities Corumbá (on the Bolivian border) and Ponta Porã (on the Paraguayan border) (Fig. 1b, Fig. 3a). The cluster had a most recent common ancestor in 1996 (95% HPD: 1989-2003; Fig. 3a), indicating that the drug-susceptible clone had circulated for approximately 23 years before the most recent samples were collected. 103 isolates were from people incarcerated at the time of diagnosis, 35 from people formerly incarcerated, and 32 from people with no incarceration history at the time of tuberculosis notification.

Figure 3.
  • Download figure
  • Open in new tab
Figure 3. A single, recently emerged M. tuberculosis clone spans prisons and the community across the major cities Mato Grosso do Sul state.

(a) A Bayesian time-calibrated phylogeny of the largest sampled genomic cluster of 170 isolates. The cluster emerged approximately 23 years before the most recently sampled isolate, with a most recent common ancestor in 1996 (95% HPD: 1989-2003). Tip point color indicates patient’s incarceration status at the time of tuberculosis notification. Annotation bar colors indicate city. Grey error bars indicate the 95% Bayesian highest posterior density intervals for node date. Clade posterior support values are shown on the middle of branches for clades with posterior support > 0.5. (b) A haplotype network of the single largest genomic cluster, including 170 isolates. Nodes represent unique haplotypes and are scaled to size. Points along branches indicate SNP distances between haplotypes. Node color indicates incarceration status at the time of diagnosis.

As observed in other clusters, isolates from people who are currently and formerly incarcerated were closely genetically related—and often, identical—to isolates from people who were not incarcerated at the time of notification, suggesting they are linked through recent transmission (Fig. 3b).

We hypothesized that the criminal justice system’s frequent transfer of people between prisons and jails could disseminate TB across the state. Brazil’s national criminal law defines different “regimes” or stages of incarceration(30), which may facilitate the potential for spillover and dissemination of infection. People are incarcerated under closed regimes, within prisons; semi-open regimes, in which people may work outside prisons and return at night; and open regimes, under which people serve sentences outside of prison, but are required to make periodic court appearances.(30) We analyzed movements in the state’s incarceration database and found that the average duration of incarceration was 230 days in a prison and 25 days in a jail. The incarceration database documents 7,982 mean yearly transfers between prisons (including closed and semi-open prisons) from 2015-2018 in a prison population of 17,221 in August 2018, including frequent transfers between cities across the state (Fig. 4a).

Figure 4.
  • Download figure
  • Open in new tab
Figure 4. Prisons disseminate tuberculosis through the frequent movement of incarcerated people between prisons.

(a) Incarceration locates people within highly connected contact networks. Map of Mato Grosso do Sul state with orange points indicating the 21 prisons which most frequently transfer people between prisons. Edges are colored by the mean number of yearly transfers made between prisons from 2015-2018, including transfers in both directions. Edges were drawn between prison pairs with more than 10 transfers in the study period. Transfers between prisons within a city are not shown. (b) An M. tuberculosis clone (depicted in Figure 3) spreads throughout multiple prisons, jails, and the community. Each line represents an individual patient infected with an isolate in the largest genomic cluster, ordered by notification date. Line color indicates patient location, determined by matching patient names with the state incarceration database (Methods). Points indicate TB notification dates; point color indicates patient location at the time of TB notification. Individuals with notification dates after December 31, 2018 are not shown.

To investigate the role of the criminal justice system in disseminating infection, we tracked the spread of the largest sampled M. tuberculosis cluster (Fig. 3) across prisons in Mato Grosso do Sul (Fig. 4b). From early 2011 to 2014, sampled isolates from the cluster were concentrated in the maximum security prison in Dourados, after which it was identified within two prisons in Campo Grande and then exported to prisons and jails in other cities. The cluster was not contained within the state’s prison network; notified tuberculosis cases among community members outside of prison occurred over the duration of the clusters’ state-wide spread.

Transmission trees reveal who-infected-whom, where, and when transmission occurred

We further investigated genomic clusters by integrating timed phylogenies with epidemiological information to infer transmission trees, allowing us to reconstruct not only transmission linkages of who-infected-whom, but also locate when and where transmission occurred. Reconstructed transmission linkages revealed frequent transmission within prisons and spillover from prisons to the community. For example, we inferred likely transmission (89% posterior probability) between two individuals (Fig. 5a,b) and occurred between July 2014 and July 2015 with 94% posterior probability (Fig. 5c). Both the recipient and predicted infector were incarcerated in the same prison in Campo Grande at the time of the recipient’s tuberculosis notification and that transmission occurred within this prison with a posterior probability of 97%.

Figure 5.
  • Download figure
  • Open in new tab
Figure 5. Frequent transmission within prisons and spillover of infection from prisons perpetuates the tuberculosis epidemic in Mato Grosso do Sul, Brazil.

We integrated M. tuberculosis phylogenetic trees with epidemiologic data to infer probable transmission linkages among sampled patients as demonstrated by an example cluster (a-c). (a) Transmission histories can be represented by a phylogenetic trees with branches colored by host in which tips represent sampled isolates and transmission events are indicated by asterisks. (b) We summarized sampled transmission trees into who-acquired-infection-from-whom matrices in which rows represent individual recipients, columns represent individual infectors, and each cell is colored by pairwise transmission probability. (c) The posterior distribution of infection times (grey histogram) for individual d with bars below representing location information for individual d and individual c, the predicted infector. Bars are colored by location over time and circles indicate date of tuberculosis notification. (d) The population-wide WAIFW matrix describing the proportions of observed transmissions to each recipient population from each infector population, corrected for genomic sampling (Methods). Rows represent recipient populations, columns represent infector populations, and cells indicate the proportion of total infections to each recipient population from each infector population defined by their incarceration status at the time of transmission and corrected for sampling rate.

Spillover from prisons plays a disproportionate role in community transmission

To estimate the total contribution of different populations to community-wide tuberculosis transmission, we inferred transmission trees for all clustered isolates. We identified the location of both the infector and recipient over the distribution of likely infection times, allowing us to determine location-specific transmission probabilities. Finally, we summed transmission probabilities over all potential transmission pairs by the incarceration status of both the infector and recipient at the time of transmission, giving a population-wide WAIFW matrix (Fig. 5d). Because we conducted active surveillance inside prisons and not in the general population, we corrected pairwise transmission probabilities to account for higher genomic sampling proportion of incarcerated populations (Fig. S2, Methods).

The population-wide WAIFW matrix reveals that a small proportion of the state’s population, 0.84% (22,706 incarcerated people of a total population of 2,713,147 in 2017), were associated with 18% of observed transmissions to people with no incarceration history (Fig. 5d). When additionally including people who were previously incarcerated, 2.6% of the population (71,703 of 2,713,147) were associated with 39% of the transmission to people with no incarceration history (Fig. 5d), demonstrating the outsized impact of prisons on tuberculosis transmission in the state. In total, we found that 57% of all infections captured in our study (including incarcerated, formerly incarcerated and never incarcerated individuals), including 46% of infections in the community (among formerly incarcerated and never incarcerated individuals), were attributable to transmission from incarcerated or formerly incarcerated individuals, who comprise a small fraction of the state’s population.

Transmission inferences are robust to epidemiological priors

We conducted a sensitivity analysis to determine the effect of epidemiological priors (Table S2) on transmission inferences for the largest genomic cluster (170 isolates). Transmission probabilities inferred for pairs of individuals were closely correlated across generation time priors (Pearson’s r: 0.94 – 0.98, Fig. S3) and sampling time priors (Pearson’s r: 0.84 – 0.98, Fig. S4), although the total number of observed transmission events was strongly influenced by both priors (Figure S5). This suggests that, as expected, epidemiological priors impact the magnitude of inferred transmission probabilities rather than the identity of predicted transmission pairs. Further, alternative priors did not strongly influence the population-wide WAIWF matrix for the largest genomic cluster (Fig. S6). The fraction of transmission to those with no incarceration history attributable to spillover (from infectors who are currently or previously incarcerated) in the largest cluster ranges from 63.5% (medium sampling and generation times) to 72.3% (short generation time and long sampling time).

Discussion

Here, we find genomic evidence that spillover from a key, high-risk environment can drive tuberculosis transmission risk in the broader population. Controlling tuberculosis thus urgently requires reducing risk created by institutional amplifiers, which have outsized effects on the general population.(11) Our population-level transmission reconstruction framework can be applied to identify other high-risk environments or populations and quantify their contribution to community tuberculosis transmission, informing targeted public health interventions.(11)

The disproportionate role that prisons play in driving community tuberculosis transmission is not unique to Brazil. Incidence of tuberculosis in prisons is extremely high globally.(7, 8) The rapid growth in incarceration rates across Central and South America, in particular, puts an increasing population at heightened risk of tuberculosis—risk that we have found can extend into neighboring communities. Although the prevalence of antibiotic resistance is low in Central West, Brazil, prisons have also amplified drug-resistant tuberculosis in Russia(24) and elsewhere. Further, the role of prisons and other detention centers as disease reservoirs is not unique to tuberculosis. Previous studies have identified spillover of meningococcal disease(25) and SARS-CoV-2(26) from prisons to the community.

Our results emphasize that reducing the extraordinarily high transmission risk created by prison environments is an urgent priority. Public health programs need to work to reduce the extreme overcrowding within prisons; improve the unsanitary, often inhumane conditions of incarceration; and expand access to primary healthcare and nutrition. Tuberculosis control programs should expand routine active screenings during(11, 27) and following incarceration so that cases can be diagnosed early and linked to treatment. Individuals leaving prison are at heightened risk of tuberculosis and may be a particularly vulnerable population. Further, preliminary evidence suggests that isoniazid protective therapy is protective against tuberculosis infection and should be considered as a protective intervention in high-risk prison environments.(28)

Reducing this excess burden of disease will require work that extends beyond biomedical interventions. The most direct way to mitigate the excess tuberculosis risk created by prisons is for governments to identify alternatives to detention so that incarceration is used only as a last resort. Prison releases have already begun as a means of reducing infectious disease transmission risk:(29) in an attempt to reduce the risk of COVID-19 in prisons, Brazil released more than thirty thousand people early by mid-2020.(30)

Our study has several limitations. First, while we were able to enroll a majority of M. tuberculosis isolates from incarcerated people with tuberculosis in the city of Dourados, overall sampling of the state’s ongoing tuberculosis epidemic was incomplete. To minimize sampling bias, we determined genomic sampling probabilities and adjusted transmission probabilities to account for over-sampling of prisons. Our probabilistic approach allows us to account for uncertainty in transmission linkages; however, we still cannot infer the identity or incarceration history of unsampled people who contribute to transmission. Terminal branch lengths may similarly be affected by sampling bias or pathogen population size and are thus an imperfect proxy of recency of transmission.

Second, transmission trees reflect information from timed M. tuberculosis phylogenies and epidemiological priors. It is difficult to distinguish genomic sampling proportion from generation time and sampling time distributions. Because we conducted transmission inference in a tuberculosis-endemic setting, rather than in a discrete outbreak, our rate of genomic sampling was lower than in some previous tuberculosis transmission studies. To reduce uncertainty in sampling proportion, we conducted transmission inference on recent subtrees (MRCA > 2012) for which we had sampled a greater proportion of isolates. To investigate the influence of epidemiological priors on transmission inferences, we conducted a sensitivity analysis (Table S2) and found that our conclusions were largely unchanged (Supplementary Text, Figs. S3-S6).

Third, we conducted whole genome sequencing of M. tuberculosis isolates cultured from sputum and generated single consensus genomes from each participant. Excluding within-host diversity could lead to some incorrect transmission inferences. For example, if we sampled an M. tuberculosis genome that was genetically distant from the genome that was transmitted, we could incorrectly exclude transmission.

Finally, we employed an existing bioinformatic approach to identify M. tuberculosis variants and excluded repetitive genomic regions as is common practice.(31, 32) Future studies that leverage within-host M. tuberculosis diversity and/or include the diverse PE/PPE genes could provide greater resolution and reduce the uncertainty of transmission inferences.(33)

In this study, we present genomic evidence that prisons act as tuberculosis reservoirs in Central West, Brazil. The dramatic expansion of incarceration in recent decades has put an increasing population at extremely high risk of tuberculosis; this risk extends to surrounding communities. Reducing the excess tuberculosis transmission risk within prisons and other detention centers is an urgent public health priority. Our findings additionally highlight the potential for pathogen genomic surveillance to identify other key environments or core populations with a disproportionate role in pathogen transmission and guide public health priorities.

Materials and Methods

Study design

We conducted population-based tuberculosis surveillance in the state of Mato Grosso do Sul in Central West Brazil from January 2014 through May 2019. Surveillance included active screening in three of the largest prisons in the state as well as ongoing passive surveillance focused on the two largest cities in the state, Campo Grande and Dourados, and three cities at the state’s border with Paraguay and Bolivia (Fig. 1a,b; Supplementary Methods). All participants provided written consent, and this study was conducted with the approval of the Research Ethics Committee from the Federal University of Grande Dourados, Federal University of Mato Grosso do Sul and National Research Ethics Committee (CONEP) (CAAE 37237814.4.0000.5160, 2676613.3.1001.5160, and 26620619.6.0000.0021) and Stanford University Institutional Review Board (IRB-40285).

Incarceration history

To investigate the incarceration history of tuberculosis patients more closely, we obtained permission from the Mato Grosso do Sul state prison administration agency, Agência Estadual de Administração do Sistema Penitenciário, to access a database of all prison entries, exits, and transfers within the criminal justice record system, Sistema Integrado de Gestão Operacional, from January 1, 2005 through December 31, 2018. Brazil’s national criminal law defines different “regimes” or stages of incarceration(34), which may facilitate the potential for spillover of infection. People are incarcerated under closed regimes, within prisons; semi-open regimes, in which people may work outside prisons and return at night; and open regimes, under which people serve sentences outside of prison, but are required to make periodic court appearances.(34)

Whole genome sequencing

We sequenced whole genomes from cultured isolates on an Illumina NextSeq (2 × 151-bp). Sequence data is available on the Sequence Read Archive (SRA), in BioProject PRJNA671770. We applied variant calling methods closely following those described in Menardo et al.(31) to be consistent with the methods used for molecular clock estimation (Supplementary Methods).

Phylogenetic and Bayesian evolutionary analysis

We fit maximum likelihood trees with RAxML-ng 1.0.1.(35) For each cluster including three or more isolates, we fit Bayesian trees with BEAST 2.6.2(36) with a strict clock, constant coalescent population size model using tuberculosis notification dates to calibrate tips (Supplementary Methods).

Transmission inference

M. tuberculosis phylogenetic trees represent patterns of evolutionary relatedness between the consensus bacterial genomes sampled from different individuals. Because most outbreaks are incompletely sampled (i.e. an M. tuberculosis sequence is not available for every case) and individuals may be infected with diverse populations of M. tuberculosis, phylogenies do not represent the underlying history of transmission.(19) We used TransPhylo(19) to infer transmission linkages—including unsampled hosts—that were consistent with the underlying timed phylogenies and transmission process (i.e. generation time and sampling intervals). To focus on densely sampled trees that would not be dominated by unsampled hosts, we sliced phylogenies at 2012 in order to generate subtrees with a most recent common ancestor of 2012 or later for transmission inference, including 122 subtrees each with 2 to 23 tips and comprising 56% of isolates (528/935).

By sampling from the posterior distribution of transmission trees, we estimated the posterior probability that each individual within a transmission tree was infected by each other individual in the tree. We summarized the posterior set of transmission trees in TransPhylo with who acquired infection from whom (WAIFW) matrices, W, where Wi,j is the posterior probability of transmission from individual j to individual i. Many transmission events are unobserved because of incomplete sampling; therefore, the sum of transmission probabilities to any individual is often less than one. We combined transmission trees with incarceration histories to quantify the proportion of observed transmission attributable to people with a history of incarceration. An individual’s location at the time of their tuberculosis notification may not reflect where transmission occurred. To understand where infector-recipient pairs were at the time transmission occurred, we integrated transmission trees with incarceration information (Fig. 4c). First, we used the transmission trees to infer the posterior distribution of infection times for each recipient i. We then partitioned the infector-recipient transmission probability Wi,j across the recipient’s infection time distribution, resulting in Embedded Image, the probability of transmission from individual j to individual i while j has incarceration status v and i has incarceration status u. Incarceration status was trichotomized as incarcerated, formerly incarcerated, and no incarceration history.

We then summed pairwise, location-specific transmission probabilities over all individual subtrees with respect to the incarceration status of both infector j and recipient i. We defined the population wide transmission matrix, Embedded Image, where Tu,v is the sum of the posterior probabilities of observed transmission events to individuals in population u from individuals in population v. We normalized Tu,v by dividing by the sum of posterior transmission probabilities to each recipient population u, Nu = ∑v Tu,v. We then defined the uncorrected proportion of the total posterior transmission probabilities to population u from population v, or the transmission fraction, as Embedded Image.

Genomic sampling varied across populations and specifically, we predicted that incarcerated people might have been sampled at a higher rate due to our active case finding and greater culture coverage in three prisons in the state. To correct for sampling bias, we determined the genomic sampling proportion rk, which we defined as the proportion of notified tuberculosis cases for which a high-quality genomic sequence was available, for population k at the time of tuberculosis notification. We calculated genomic sampling proportion as the number of sequenced genomes for each group divided by the number of tuberculosis notifications for each group from 2014-2019 in Campo Grande and Dourados, the two major cities in our prospective study (Fig. S2). We adjusted each pairwise transmission probability by the sampling rate of the infector, based on the incarceration status of the infector at the time of tuberculosis notification, for Embedded Image and calculated the sampling-adjusted transmission fractions as described above.

Role of the funding source

The funding source had no role in study design, data collection, data analysis, data interpretation, or writing of the report.

Statistical methods

We inferred transmission trees with TransPhylo(19) as described above. We compared phylogenetic tree terminal branch lengths between populations with two-sample t-tests and compared the rate of clustering between populations with two-sample proportion tests.

Data Availability

Sequence data is available on the Sequence Read Archive (SRA), in BioProject PRJNA671770.

Funding

National Institutes of Health grant R01AI130058 (JRA)

Brazil’s National Council for Scientific and Technological Development grant 404237/2012-6 (JC)

Author contributions

Conceptualization: JRA, JC, KSW, CC, TC, BM, AIK

Methodology: JRA, JC, KSW, CC, TC, BM, AIK

Investigation: PCPS, TOG, BOS, ASS, ACL, AMS, FMFM, RDO, EFL, EC, YL

Visualization: KSW, JRA

Funding acquisition: JRA, JC

Project administration: JRA, JC

Supervision: JRA, JC

Writing – original draft: KSW, JRA

Writing – review & editing: all authors

Competing interests

Authors declare that they have no competing interests.

Data and materials availability

Raw Illumina sequence data are available on the Sequence Read Archive under BioProject PRJNA671770.

Acknowledgments

We would like to acknowledge the contributions of AGEPEN, LACEN, and the Mato Grosso do Sul State Health Department. All authors had full access to all the data in the study and had final responsibility for the decision to submit

References

  1. 1.↵
    G. Churchyard, P. Kim, N. S. Shah, R. Rustomjee, N. Gandhi, B. Mathema, D. Dowdy, A. Kasmar, V. Cardenas, What We Know about Tuberculosis Transmission: An Overview, J. Infect. Dis. 216, S629–S635 (2017).
    OpenUrl
  2. 2.
    S. C. Auld, N. Sarita Shah, B. Mathema, T. S. Brown, N. Ismail, S. V. Omar, J. C. M. Brust, K. N. Nelson, S. Allana, A. Campbell, K. Mlisana, P. Moodley, N. R. Gandhi, Extensively drug-resistant tuberculosis in South Africa: Genomic evidence supporting transmission in communities, Eur. Respir. J. 52 (2018), doi:10.1183/13993003.00246-2018.
    OpenUrlAbstract/FREE Full Text
  3. 3.
    K. Middelkoop, B. Mathema, L. Myer, E. Shashkina, A. Whitelaw, G. Kaplan, B. Kreiswirth, R. Wood, L. G. Bekker, Transmission of tuberculosis in a south African community with a high prevalence of HIV infection, J. Infect. Dis. 211, 53–61 (2015).
    OpenUrlCrossRefPubMed
  4. 4.
    J. R. Andrews, C. Morrow, R. P. Walensky, R. Wood, Integrating social contact and environmental data in evaluating tuberculosis transmission in a South African township., J. Infect. Dis. 210, 597–603 (2014).
    OpenUrlCrossRefPubMed
  5. 5.
    T. A. Yates, P. Y. Khan, G. M. Knight, J. G. Taylor, T. D. McHugh, M. Lipman, R. G. White, T. Cohen, F. G. Cobelens, R. Wood, D. A. J. Moore, I. Abubakar, The transmission of Mycobacterium tuberculosis in high burden settings, Lancet Infect. Dis. 16, 227–238 (2016).
    OpenUrlCrossRef
  6. 6.↵
    J. R. Andrews, C. Morrow, R. Wood, Modeling the role of public transportation in sustaining tuberculosis transmission in South Africa, Am. J. Epidemiol. 177, 556–561 (2013).
    OpenUrlCrossRefPubMedWeb of Science
  7. 7.↵
    K. Dolan, A. L. Wirtz, B. Moazen, M. Ndeffo-mbah, A. Galvani, S. A. Kinner, R. Courtney, M. McKee, J. J. Amon, L. Maher, M. Hellard, C. Beyrer, F. L. Altice, Global burden of HIV, viral hepatitis, and tuberculosis in prisoners and detainees, Lancet 388, 1089–1102 (2016).
    OpenUrlCrossRefPubMed
  8. 8.↵
    O. Cords, L. Martinez, J. Warren, J. O’Marr, K. S. Walter, T. Cohen, J. Zheng, A. I. Ko, J. Croda, J. R. Andrews, Incidence and Prevalence of Tuberculosis in Incarcerated Populations: A Systematic Review and Meta-analysis, Lancet Public Heal..
  9. 9.↵
    J. Urrego, A. I. Ko, A. Da Silva Santos Carbone, D. S. G. Paião, R. V. E. Sgarbi, C. W. Yeckel, J. R. Andrews, J. Croda, The impact of ventilation and early diagnosis on tuberculosis transmission in Brazilian prisons, Am. J. Trop. Med. Hyg. 93, 739–746 (2015).
    OpenUrlAbstract/FREE Full Text
  10. 10.↵
    S. Basu, D. Stuckler, M. McKee, Addressing institutional amplifiers in the dynamics and control of tuberculosis epidemics, Am. J. Trop. Med. Hyg. 84, 30–37 (2011).
    OpenUrlAbstract/FREE Full Text
  11. 11.↵
    T. S. Mabud, M. de Lourdes Delgado Alves, A. I. Ko, S. Basu, K. S. Walter, T. Cohen, B. Mathema, C. Colijn, E. Lemos, J. Croda, J. R. Andrews, J. Z. Metcalfe, Ed. Evaluating strategies for control of tuberculosis in prisons and prevention of spillover into communities: An observational and modeling study from Brazil, PLoS Med. 16, e1002737 (2019).
    OpenUrlCrossRefPubMed
  12. 12.↵
    K. Walter, L. Martinez, Arakaki-Sanchez, V. Sequera, S. G. Estigarribia, T. Cohen, A. Ko, A. García-Basteiro, Z. Rueda, O. R. Lopez, M. Espinal, J. Croda, J. Andrews, The escalating tuberculosis crisis in Central and South American prisons. In press., Lancet (2020), doi:10.1016/S0140-6736(20)32578-2.
    OpenUrlCrossRef
  13. 13.↵
    F. P. C. Sacchi, R. M. Praça, M. B. Tatara, V. Simonsen, L. Ferrazoli, M. G. Croda, P. N. Suffys, A. I. Ko, J. R. Andrews, J. Croda, Prisons as reservoir for community transmission of tuberculosis, brazil, Emerg. Infect. Dis. 21, 452–455 (2015).
    OpenUrlCrossRef
  14. 14.↵
    R. Coker, M. McKee, R. Atun, B. Dimitrova, E. Dodonova, S. Kuznetsov, F. Drobniewski, Risk factors for pulmonary tuberculosis in Russia: Case-control study, Br. Med. J. 332, 85–87 (2006).
    OpenUrl
  15. 15.↵
    J. L. Warren, L. Grandjean, D. A. J. Moore, A. Lithgow, J. Coronel, P. Sheen, J. L. Zelner, J. R. Andrews, T. Cohen, Investigating spillover of multidrug-resistant tuberculosis from a prison: a spatial and molecular epidemiological analysis, BMC Med. 16, 122 (2018).
    OpenUrlCrossRef
  16. 16.↵
    H. E. Jenkins, V. Plesca, A. Ciobanu, V. Crudu, I. Galusca, V. Soltan, A. Serbulenco, M. Zignol, A. Dadu, M. Dara, T. Cohen, Assessing spatial heterogeneity of multidrug-resistant tuberculosis in a high-burden country, Eur. Respir. J. 42, 1291–1301 (2013).
    OpenUrlAbstract/FREE Full Text
  17. 17.↵
    M. N. Séraphin, X. Didelot, D. J. Nolan, J. R. May, M. S. R. Khan, E. R. Murray, M. Salemi, J. G. Morris, M. Lauzardo, Genomic investigation of a Mycobacterium tuberculosis outbreak involving prison and community cases in Florida, United States, Am. J. Trop. Med. Hyg. 99, 867– 874 (2018).
    OpenUrl
  18. 18.↵
    T. F. Jones, C. L. Woodley, F. F. Fountain, W. Schaffner, Increased incidence of the outbreak strain of Mycobacterium tuberculosis in the surrounding community after an outbreak in a jail. (Original Article), South. Med. J. 96, 155–158 (2003).
    OpenUrlPubMedWeb of Science
  19. 19.↵
    X. Didelot, C. Fraser, J. Gardy, C. Colijn, H. Malik, Genomic infectious disease epidemiology in partially sampled and ongoing outbreaks, Mol. Biol. Evol. 34, 997–1007 (2017).
    OpenUrl
  20. 20.↵
    World Health Organization, Global tuberculosis report 2020 (Geneva, 2020).
  21. 21.↵
    K. E. Holt, P. Mcadam, P. Vuong, K. Thai, N. Thuy, T. Thuong, D. Thi, M. Ha, N. N. Lan, N. H. Lan, N. Thi, Q. Nhu, H. T. Hai, V. Thi, N. Ha, G. Thwaites, D. J. Edwards, A. P. Nath, K. Pham, D. B. Ascher, J. Farrar, C. C. Khor, Y. Y. Teo, M. Inouye, Frequent transmission of the Mycobacterium tuberculosis Beijing lineage and positive selection for the EsxW Beijing variant in Vietnam, Nat. Genet. (2018), doi:10.1038/s41588-018-0117-9.
    OpenUrlCrossRef
  22. 22.↵
    T. M. Walker, T. A. Kohl, S. V. Omar, J. Hedge, C. Del Ojo Elias, P. Bradley, Z. Iqbal, S. Feuerriegel, K. E. Niehaus, D. J. Wilson, D. A. Clifton, G. Kapatai, C. L. C. Ip, R. Bowden, F. A. Drobniewski, C. Allix-Béguec, C. Gaudin, J. Parkhill, R. Diel, P. Supply, D. W. Crook, E. G. Smith, A. S. Walker, N. Ismail, S. Niemann, T. E. A. Peto, J. Davies, C. Crichton, M. Acharya, L. Madrid-Marquez, D. Eyre, D. Wyllie, T. Golubchik, M. Munang, Whole-genome sequencing for prediction of Mycobacterium tuberculosis drug susceptibility and resistance: A retrospective cohort study, Lancet Infect. Dis. 15, 1193–1202 (2015).
    OpenUrlCrossRefPubMed
  23. 23.↵
    T. M. Walker, C. L. C. Ip, R. H. Harrell, J. T. Evans, G. Kapatai, M. J. Dedicoat, D. W. Eyre, D. J. Wilson, P. M. Hawkey, D. W. Crook, J. Parkhill, D. Harris, A. S. Walker, R. Bowden, P. Monk, E. G. Smith, T. E. A. Peto, Whole-genome sequencing to delineate Mycobacterium tuberculosis outbreaks: A retrospective observational study, Lancet Infect. Dis. 13, 137–146 (2013).
    OpenUrlCrossRefPubMedWeb of Science
  24. 24.↵
    F. Drobniewski, Y. Balabanova, V. Nikolayevsky, M. Ruddy, S. Kuznetzov, S. Zakharova, A. Melentyev, I. Fedorin, Drug-resistant tuberculosis, clinical virulence, and the dominance of the Beijing strain family in Russia, J. Am. Med. Assoc. 293, 2726–2731 (2005).
    OpenUrlCrossRefPubMedWeb of Science
  25. 25.↵
    J. W. Tappero, R. Reporter, J. D. Wenger, B. A. Ward, M. W. Reeves, T. S. Missbach, B. D. Plikaytis, L. Mascola, A. Schuchat, Meningococcal Disease in Los Angeles County, California, and among Men in the County Jails, N. Engl. J. Med. 335, 833–841 (1996).
    OpenUrlCrossRefPubMedWeb of Science
  26. 26.↵
    E. Reinhart, D. L. Chen, Incarceration And Its Disseminations: COVID-19 Pandemic Lessons From Chicago’s Cook County Jail, Health Aff. 39, 1412–1418 (2020).
    OpenUrlPubMed
  27. 27.↵
    A. da S. Santos, R. D. de Oliveira, E. F. Lemos, F. Lima, T. Cohen, O. Cords, L. Martinez, C. Gonçalves, A. I. Ko, J. R. Andrews, J. Croda, Yield, Efficiency, and Costs of Mass Screening Algorithms for Tuberculosis in Brazilian Prisons, Clin. Infect. Dis. (2020), doi:10.1093/cid/ciaa135.
    OpenUrlCrossRef
  28. 28.↵
    R. D. de Oliveira, A. da Silva Santos, C. B. Reis, A. de Cássia Leite, F. P. C. Sacchi, R. C. P. de Araujo, P. C. P. dos Santos, V. C. Rolla, L. Martinez, J. Andrews, J. Croda, Primary prophylaxis to prevent tuberculosis infection in prison inmates: A randomized, double-blind, placebo-controlled trial, Am. J. Trop. Med. Hyg. 103, 1466–1472 (2020).
    OpenUrl
  29. 29.↵
    E. A. Wang, B. Western, D. M. Berwick, COVID-19, Decarceration, and the Role of Clinicians, Health Systems, and Payers: A Report From the National Academy of Sciences, Engineering, and Medicine., Jama 06511 (2020), doi:10.1001/jama.2020.22109.
    OpenUrlCrossRef
  30. 30.↵
    R. De Oliveira Andrade, Covid-19: Prisons exposed in Brazil’s crisis BMJ 370 (2020), doi:10.1136/bmj.m2884.
    OpenUrlFREE Full Text
  31. 31.↵
    F. Menardo, S. Duchêne, D. Brites, S. Gagneux, The molecular clock of mycobacterium tuberculosis, PLoS Pathog. 15 (2019), doi:10.1371/journal.ppat.1008067.
    OpenUrlCrossRef
  32. 32.↵
    C. J. Meehan, G. A. Goig, T. A. Kohl, L. Verboven, A. Dippenaar, M. Ezewudo, M. R. Farhat, J. L. Guthrie, K. Laukens, P. Miotto, B. Ofori-Anyinam, V. Dreyer, P. Supply, A. Suresh, C. Utpatel, D. van Soolingen, Y. Zhou, P. M. Ashton, D. Brites, A. M. Cabibbe, B. C. de Jong, M. de Vos, F. Menardo, S. Gagneux, Q. Gao, T. H. Heupink, Q. Liu, C. Loiseau, L. Rigouts, T. C. Rodwell, E. Tagliani, T. M. Walker, R. M. Warren, Y. Zhao, M. Zignol, M. Schito, J. Gardy, D. M. Cirillo, S. Niemann, I. Comas, A. Van Rie, Whole genome sequencing of Mycobacterium tuberculosis: current standards and open issues, Nat. Rev. Microbiol. 17, 533–545 (2019).
    OpenUrl
  33. 33.↵
    K. S. Walter, C. Colijn, T. Cohen, B. Mathema, Q. Liu, J. Bowers, D. M. Engelthaler, A. Narechania, D. Lemmer, J. Croda, J. R. Andrews, Genomic variant-identification methods may alter mycobacterium tuberculosis transmission inferences, Microb. Genomics 6, 1–16 (2020).
    OpenUrl
  34. 34.↵
    J. Mariner, J. Cavallaro, Behind Bars in Brazil (Human Rights Watch, New York, 1998).
  35. 35.↵
    A. M. Kozlov, D. Darriba, T. Flouri, B. Morel, A. Stamatakis, J. Wren, Ed. RAxML-NG: A fast, scalable, and user-friendly tool for maximum likelihood phylogenetic inference, Bioinformatics (2019), doi:10.1093/bioinformatics/btz305.
    OpenUrlCrossRefPubMed
  36. 36.↵
    R. Bouckaert, T. G. Vaughan, J. Barido-Sottani, S. Duchêne, M. Fourment, A. Gavryushkina, J. Heled, G. Jones, D. Kühnert, N. De Maio, M. Matschiner, F. K. Mendes, N. F. Müller, H. A. Ogilvie, L. Du Plessis, A. Popinga, A. Rambaut, D. Rasmussen, I. Siveroni, M. A. Suchard, C. H. Wu, D. Xie, C. Zhang, T. Stadler, A. J. Drummond, BEAST 2.5: An advanced software platform for Bayesian evolutionary analysis, PLoS Comput. Biol. 15, e1006650 (2019).
    OpenUrlCrossRefPubMed
  37. 37.
    A. da S. Santos, R. D. de Oliveira, E. F. Lemos, F. Lima, T. Cohen, O. Cords, L. Martinez, C. Gonçalves, A. Ko, J. R. Andrews, J. Croda, Yield, Efficiency and Costs of Mass Screening Algorithms for Tuberculosis in Brazilian Prisons., Clin. Infect. Dis. (2020), doi:10.1093/cid/ciaa135.
    OpenUrlCrossRef
  38. 38.
    D. Robinson, fuzzyjoin: Join Tables Together on Inexact Matching (2020) (available at https://github.com/dgrtwo/fuzzyjoin).
  39. 39.
    F. Krueger, Trim Galore (2019) (available at https://github.com/FelixKrueger/TrimGalore/blob/master/Docs/Trim_Galore_User_Guide.md).
  40. 40.
    M. Martin, Cutadapt Removes Adapter Sequences From High-Throughput Sequencing Reads, EMBnet.journal 17 (2011), doi:https://doi.org/10.14806/ej.17.1.200.
  41. 41.
    D. E. Wood, S. L. Salzberg, Kraken: Ultrafast metagenomic sequence classification using exact alignments, Genome Biol. 15, R46 (2014).
    OpenUrlCrossRefPubMed
  42. 42.
    H. Li, R. Durbin, Fast and accurate short read alignment with Burrows-Wheeler transform, Bioinformatics 25, 1754–1760 (2009).
    OpenUrlCrossRefPubMedWeb of Science
  43. 43.
    D. Brites, C. Loiseau, F. Menardo, S. Borrell, M. B. Boniotti, R. Warren, A. Dippenaar, S. D. C. Parsons, C. Beisel, M. A. Behr, J. A. Fyfe, M. Coscolla, S. Gagneux, A new phylogenetic framework for the animal-adapted mycobacterium tuberculosis complex, Front. Microbiol. 9, 2820 (2018).
    OpenUrlCrossRef
  44. 44.
    C. Anyansi, A. Keo, B. J. Walker, T. J. Straub, A. L. Manson, A. M. Earl, T. Abeel, QuantTB-A method to classify mixed Mycobacterium tuberculosis infections within whole genome sequencing data, BMC Genomics 21, 80 (2020).
    OpenUrl
  45. 45.
    A. J. Page, S. R. Harris, T. Seemann, B. Taylor, A. J. Delaney, J. A. Keane, J. Soares, SNP-sites: rapid efficient extraction of SNPs from multi-FASTA alignments, Microb. Genomics 2, 1–5 (2016).
    OpenUrl
  46. 46.
    P. Bradley, N. C. Gordon, T. M. Walker, L. Dunn, S. Heys, B. Huang, S. Earle, L. J. Pankhurst, L. Anson, M. De Cesare, P. Piazza, A. A. Votintseva, T. Golubchik, D. J. Wilson, D. H. Wyllie, R. Diel, S. Niemann, S. Feuerriegel, T. A. Kohl, N. Ismail, S. V Omar, E. G. Smith, D. Buck, G. McVean, A. S. Walker, T. E. A. Peto, D. W. Crook, Z. Iqbal, Rapid antibiotic-resistance predictions from genome sequence data for Staphylococcus aureus and Mycobacterium tuberculosis, Nat. Commun. 6 (2015), doi:10.1038/ncomms10063.
    OpenUrlCrossRefPubMed
  47. 47.
    C. Allix-Béguec, I. Arandjelovic, L. Bi, P. Beckert, M. Bonnet, P. Bradley, A. M. Cabibbe, I. Cancino-Muñoz, M. J. Caulfield, A. Chaiprasert, D. M. Cirillo, D. A. Clifton, I. Comas, D. W. Crook, M. R. De Filippo, H. de Neeling, R. Diel, F. A. Drobniewski, K. Faksri, M. R. Farhat, J. Fleming, P. Fowler, T. A. Fowler, Q. Gao, J. Gardy, D. Gascoyne-Binzi, A. L. Gibertoni-Cruz, A. Gil-Brusola, T. Golubchik, X. Gonzalo, L. Grandjean, G. He, J. L. Guthrie, S. Hoosdally, M. Hunt, Z. Iqbal, N. Ismail, J. Johnston, F. M. Khanzada, C. C. Khor, T. A. Kohl, C. Kong, S. Lipworth, Q. Liu, G. Maphalala, E. Martinez, V. Mathys, M. Merker, P. Miotto, N. Mistry, D. A. J. Moore, M. Murray, S. Niemann, R. T. H. Ong, T. E. A. Peto, J. E. Posey, T. Prammananan, A. Pym, C. Rodrigues, M. Rodrigues, T. Rodwell, G. M. Rossolini, E. S. Padilla, M. Schito, X. Shen, J. Shendure, V. Sintchenko, A. Sloutsky, E. G. Smith, M. Snyder, K. Soetaert, A. M. Starks, P. Supply, P. Suriyapol, S. Tahseen, P. Tang, Y. Y. Teo, T. N. T. Thuong, G. Thwaites, E. Tortoli, S. V. Omar, D. van Soolingen, A. S. Walker, T. M. Walker, M. Wilcox, D. J. Wilson, D. Wyllie, Y. Yang, H. Zhang, Y. Zhao, B. Zhu, Prediction of susceptibility to first-line tuberculosis drugs by DNA sequencing, N. Engl. J. Med. 379, 1403–1415 (2018).
    OpenUrlCrossRefPubMed
  48. 48.
    J. E. Phelan, D. M. O’Sullivan, D. Machado, J. Ramos, Y. E. A. Oppong, S. Campino, J. O’Grady, R. McNerney, M. L. Hibberd, M. Viveiros, J. F. Huggett, T. G. Clark, Integrating informatics tools and portable sequencing technology for rapid detection of resistance to anti-tuberculous drugs, Genome Med. 11, 41 (2019).
    OpenUrlCrossRef
  49. 49.
    E. Paradis, K. Schliep, R. Schwartz, Ed. ape 5.0: an environment for modern phylogenetics and evolutionary analyses in R, Bioinformatics 35, 526–528 (2019).
    OpenUrlCrossRefPubMed
  50. 50.
    A. D. Leaché, B. L. Banbury, J. Felsenstein, A. Nieto-Montes De Oca, A. Stamatakis, New Acquisition Bias Corrections for Inferring SNP Phylogenies, Syst. Biol 64, 1032–1047 (2015).
    OpenUrlCrossRefPubMed
  51. 51.
    A. Rambaut, A. J. Drummond, D. Xie, G. Baele, M. A. Suchard, Posterior summarization in Bayesian phylogenetics using Tracer 1.7, Syst. Biol. 67, 901–904 (2018).
    OpenUrlCrossRefPubMed
  52. 52.
    A. Rieux, C. E. Khatchikian, tipdatingbeast: an r package to assist the implementation of phylogenetic tip-dating tests using beast, Mol. Ecol. Resour. 17, 608–613 (2017).
    OpenUrlCrossRef
  53. 53.
    Y. I. Xu, I. Cancino-Muñoz, M. I. Torres-Puente, L. M. Villamayor, R. I. Borrá, M. Borrá s-Má ñez, M. Bosque, J. J. Camarena ID, E. Colomer-Roig, J. I. Colomina, I. Escribano, O. I. Esparcia-Rodríguez, A. I. Gil-Brusola, R. Moreno-Muñoz, D. NavarroID, M. Navarro, N. Orta, E. Pé rez, J. Prat, J. I. Carlos Rodríguez, M. Montserrat Ruiz-García, H. Vanaclocha, C. ColijnID, I. I. Comas, I. Cancino-Munoz, M. I. Torres-Puente, L. M. Villamayor, R. Borrás, M. Borrás-Máñez, M. Bosque, J. J. Camarena, E. Colomer-Roig, J. I. Colomina, I. Escribano, O. I. Esparcia-Rodríguez, A. I. Gil-Brusola, C. Gimeno, A. Gimeno-Gascón, B. Gomila-Sard, D. González-Granda, N. Gonzalo-Jiménez, M. R. Guna-Serrano, J. L. López-Hontangas, C. Martín-González, R. Moreno-Munoz, D. Navarro, M. Navarro, N. Orta, E. Pérez, J. Prat, J. C. Rodríguez, M. M. Ruiz-García, H. Vanaclocha, C. Colijn, I. I. Comas, High-resolution mapping of tuberculosis transmission: Whole genome sequencing and phylogenetic modelling of a cohort from Valencia Region, Spain, PLoS Med. 16 (2019), doi:10.1371/journal.pmed.1002961.
    OpenUrlCrossRef
  54. 54.
    D. Ayabina, J. O. Ronning, K. Alfsnes, N. Debech, O. B. Brynildsrud, T. Arnesen, G. Norheim, A.-T. T. Mengshoel, R. Rykkvin, U. R. Dahle, C. Colijn, V. Eldholm, Genome-based transmission modelling separates imported tuberculosis from recent transmission within an immigrant population, Microb. Genomics 4, 1–13 (2018).
    OpenUrl
Back to top
PreviousNext
Posted June 24, 2021.
Download PDF

Supplementary Material

Data/Code
Email

Thank you for your interest in spreading the word about medRxiv.

NOTE: Your email address is requested solely to identify you as the sender of this article.

Enter multiple addresses on separate lines or separate them with commas.
Genomic evidence for prisons as amplifiers of community tuberculosis epidemics
(Your Name) has forwarded a page to you from medRxiv
(Your Name) thought you would like to see this page from the medRxiv website.
CAPTCHA
This question is for testing whether or not you are a human visitor and to prevent automated spam submissions.
Share
Genomic evidence for prisons as amplifiers of community tuberculosis epidemics
Katharine S. Walter, Paulo César Pereira dos Santos, Thais Oliveira Gonçalves, Brunna Oliveira da Silva, Andrea da Silva Santos, Alessandra de Cássia Leite, Alessandra Moura da Silva, Flora Martinez Figueira Moreira, Roberto Dias de Oliveira, Everton Ferreira Lemos, Eunice Cunha, Yiran E. Liu, Albert Ko, Caroline Colijn, Ted Cohen, Barun Mathema, Julio Croda, Jason R. Andrews
medRxiv 2021.06.22.21259360; doi: https://doi.org/10.1101/2021.06.22.21259360
Digg logo Reddit logo Twitter logo Facebook logo Google logo LinkedIn logo Mendeley logo
Citation Tools
Genomic evidence for prisons as amplifiers of community tuberculosis epidemics
Katharine S. Walter, Paulo César Pereira dos Santos, Thais Oliveira Gonçalves, Brunna Oliveira da Silva, Andrea da Silva Santos, Alessandra de Cássia Leite, Alessandra Moura da Silva, Flora Martinez Figueira Moreira, Roberto Dias de Oliveira, Everton Ferreira Lemos, Eunice Cunha, Yiran E. Liu, Albert Ko, Caroline Colijn, Ted Cohen, Barun Mathema, Julio Croda, Jason R. Andrews
medRxiv 2021.06.22.21259360; doi: https://doi.org/10.1101/2021.06.22.21259360

Citation Manager Formats

  • BibTeX
  • Bookends
  • EasyBib
  • EndNote (tagged)
  • EndNote 8 (xml)
  • Medlars
  • Mendeley
  • Papers
  • RefWorks Tagged
  • Ref Manager
  • RIS
  • Zotero
  • Tweet Widget
  • Facebook Like
  • Google Plus One

Subject Area

  • Epidemiology
Subject Areas
All Articles
  • Addiction Medicine (215)
  • Allergy and Immunology (495)
  • Anesthesia (106)
  • Cardiovascular Medicine (1093)
  • Dentistry and Oral Medicine (195)
  • Dermatology (141)
  • Emergency Medicine (274)
  • Endocrinology (including Diabetes Mellitus and Metabolic Disease) (499)
  • Epidemiology (9757)
  • Forensic Medicine (5)
  • Gastroenterology (480)
  • Genetic and Genomic Medicine (2303)
  • Geriatric Medicine (222)
  • Health Economics (462)
  • Health Informatics (1553)
  • Health Policy (732)
  • Health Systems and Quality Improvement (602)
  • Hematology (236)
  • HIV/AIDS (501)
  • Infectious Diseases (except HIV/AIDS) (11631)
  • Intensive Care and Critical Care Medicine (616)
  • Medical Education (236)
  • Medical Ethics (67)
  • Nephrology (256)
  • Neurology (2139)
  • Nursing (134)
  • Nutrition (335)
  • Obstetrics and Gynecology (426)
  • Occupational and Environmental Health (517)
  • Oncology (1172)
  • Ophthalmology (363)
  • Orthopedics (128)
  • Otolaryngology (220)
  • Pain Medicine (145)
  • Palliative Medicine (50)
  • Pathology (309)
  • Pediatrics (694)
  • Pharmacology and Therapeutics (298)
  • Primary Care Research (265)
  • Psychiatry and Clinical Psychology (2172)
  • Public and Global Health (4645)
  • Radiology and Imaging (775)
  • Rehabilitation Medicine and Physical Therapy (455)
  • Respiratory Medicine (623)
  • Rheumatology (274)
  • Sexual and Reproductive Health (225)
  • Sports Medicine (208)
  • Surgery (250)
  • Toxicology (43)
  • Transplantation (120)
  • Urology (94)