1
|
Oliveira RC, Souza JSMD, Alcântara LCJ, Guimarães ML, Brites C, Monteiro-Cunha JP. Molecular and Phylogenetic Analysis of HIV-1 Subtype C in Bahia State, Northeastern Brazil. AIDS Res Hum Retroviruses 2024; 40:37-41. [PMID: 37312563 DOI: 10.1089/aid.2023.0038] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/15/2023] Open
Abstract
HIV-1 subtype C is associated with more than half of infections in southern Brazil and has been increasing in other regions of the country. In a previous study carried out in northeastern Brazil, we found a prevalence of 4.1% of subtype C. This work investigates the origin of subtype C in the state of Bahia based on five new viral sequences. The phylogenetic analysis showed that subtype C viruses found in Bahia descend from the main lineage that circulates in other Brazilian regions.
Collapse
Affiliation(s)
- Rodrigo Cunha Oliveira
- Departamento de Bioquímica e Biofísica, Núcleo de Bioinformática, Universidade Federal da Bahia, Salvador, Brasil
| | | | | | | | - Carlos Brites
- Laboratório de Pesquisa em Infectologia, Complexo Hospitalar Prof. Edgard Santos, Universidade Federal da Bahia, Salvador, Brasil
| | - Joana Paixão Monteiro-Cunha
- Departamento de Bioquímica e Biofísica, Núcleo de Bioinformática, Universidade Federal da Bahia, Salvador, Brasil
| |
Collapse
|
2
|
de Almeida Rego FF, de Moraes L, Giovanetti M, Silva JAG, Torres FG, de Oliveira Silva M, da Purificação Pereira da Silva M, Van Weyenbergh J, Santos LA, Khouri R. Genomic Detection of the Emerging, Highly Pathogenic HIV-1 Subtype D in Bahia, Northeast Brazil. Viruses 2023; 15:1650. [PMID: 37631993 PMCID: PMC10458544 DOI: 10.3390/v15081650] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/31/2023] [Revised: 06/30/2023] [Accepted: 07/06/2023] [Indexed: 08/27/2023] Open
Abstract
(1) Background: The HIV subtype D is generally associated with a faster decline in CD4+ T cell counts, a higher viral load, and a faster progression to AIDS. However, it is still poorly characterized in Brazil. In this study, we used genomics and epidemiological data to investigate the transmission dynamics of HIV subtype D in the state of Bahia, Northeast Brazil. (2) Methods: To achieve this goal, we obtained four novel HIV-1 subtype D partial pol genome sequences using the Sanger method. To understand the emergence of this novel subtype in the state of Bahia, we used phylodynamic analysis on a dataset comprising 3704 pol genome sequences downloaded from the Los Alamos database. (3) Results: Our analysis revealed three branching patterns, indicating multiple introductions of the HIV-1 subtype D in Brazil from the late 1980s to the late 2000s and a single introduction event in the state of Bahia. Our data further suggest that these introductions most likely originated from European, Eastern African, Western African, and Southern African countries. (4) Conclusion: Understanding the distribution of HIV-1 viral strains and their temporal dynamics is crucial for monitoring the real-time evolution of circulating subtypes and recombinant forms, as well as for designing novel diagnostic and vaccination strategies. We advocate for a shift to active surveillance, to ensure adequate preparedness for future epidemics mediated by emerging viral strains.
Collapse
Affiliation(s)
- Filipe Ferreira de Almeida Rego
- Escola Bahiana de Medicina e Saúde Pública, Avenida Dom João VI, 275, Brotas, Salvador 40290-000, Bahia, Brazil; (F.F.d.A.R.); (L.A.S.)
| | - Laise de Moraes
- Instituto Gonçalo Moniz, Fundação Oswaldo Cruz, Rua Waldemar Falcão, 121, Candeal, Salvador 40296-710, Bahia, Brazil; (L.d.M.); (J.A.G.S.); (F.G.T.)
- Programa de Pós-Graduação em Ciências da Saúde, Faculdade de Medicina da Bahia, Universidade Federal da Bahia, Praça Ramos de Queirós, s/n, Largo do Terreiro de Jesus, Salvador 40026-010, Bahia, Brazil
| | - Marta Giovanetti
- Instituto Rene Rachou, Fundação Oswaldo Cruz, Avenida Augusto de Lima, 1715, Barro Preto, Belo Horizonte 30190-002, Minas Gerais, Brazil;
- Sciences and Technologies for Sustainable Development and One Health, University of Campus Bio-Medico, 21-00128 Rome, Italy
| | - José Adriano Góes Silva
- Instituto Gonçalo Moniz, Fundação Oswaldo Cruz, Rua Waldemar Falcão, 121, Candeal, Salvador 40296-710, Bahia, Brazil; (L.d.M.); (J.A.G.S.); (F.G.T.)
- Programa de Pós-Graduação em Ciências da Saúde, Faculdade de Medicina da Bahia, Universidade Federal da Bahia, Praça Ramos de Queirós, s/n, Largo do Terreiro de Jesus, Salvador 40026-010, Bahia, Brazil
- Centro Estadual Especializado em Diagnóstico, Assistência e Pesquisa, Secretaria de Saúde do Estado da Bahia, Salvador 40110-050, Bahia, Brazil; (M.d.O.S.); (M.d.P.P.d.S.)
| | - Felipe Guimarães Torres
- Instituto Gonçalo Moniz, Fundação Oswaldo Cruz, Rua Waldemar Falcão, 121, Candeal, Salvador 40296-710, Bahia, Brazil; (L.d.M.); (J.A.G.S.); (F.G.T.)
| | - Marcio de Oliveira Silva
- Centro Estadual Especializado em Diagnóstico, Assistência e Pesquisa, Secretaria de Saúde do Estado da Bahia, Salvador 40110-050, Bahia, Brazil; (M.d.O.S.); (M.d.P.P.d.S.)
| | - Maria da Purificação Pereira da Silva
- Centro Estadual Especializado em Diagnóstico, Assistência e Pesquisa, Secretaria de Saúde do Estado da Bahia, Salvador 40110-050, Bahia, Brazil; (M.d.O.S.); (M.d.P.P.d.S.)
| | - Johan Van Weyenbergh
- Rega Institute for Medical Research, Department of Immunology, Microbiology and Transplantation, KU Leuven, Herestraat 49-Box 1030, 3000 Leuven, Belgium;
| | - Luciane Amorim Santos
- Escola Bahiana de Medicina e Saúde Pública, Avenida Dom João VI, 275, Brotas, Salvador 40290-000, Bahia, Brazil; (F.F.d.A.R.); (L.A.S.)
- Instituto Gonçalo Moniz, Fundação Oswaldo Cruz, Rua Waldemar Falcão, 121, Candeal, Salvador 40296-710, Bahia, Brazil; (L.d.M.); (J.A.G.S.); (F.G.T.)
- Programa de Pós-Graduação em Ciências da Saúde, Faculdade de Medicina da Bahia, Universidade Federal da Bahia, Praça Ramos de Queirós, s/n, Largo do Terreiro de Jesus, Salvador 40026-010, Bahia, Brazil
| | - Ricardo Khouri
- Instituto Gonçalo Moniz, Fundação Oswaldo Cruz, Rua Waldemar Falcão, 121, Candeal, Salvador 40296-710, Bahia, Brazil; (L.d.M.); (J.A.G.S.); (F.G.T.)
- Programa de Pós-Graduação em Ciências da Saúde, Faculdade de Medicina da Bahia, Universidade Federal da Bahia, Praça Ramos de Queirós, s/n, Largo do Terreiro de Jesus, Salvador 40026-010, Bahia, Brazil
- Rega Institute for Medical Research, Department of Immunology, Microbiology and Transplantation, KU Leuven, Herestraat 49-Box 1030, 3000 Leuven, Belgium;
| |
Collapse
|
3
|
Avila-Rios S, García-Morales C, Reyes-Terán G, González-Rodríguez A, Matías-Florentino M, Mehta SR, Chaillon A. Phylodynamics of HIV in the Mexico City Metropolitan Region. J Virol 2022; 96:e0070822. [PMID: 35762759 PMCID: PMC9327710 DOI: 10.1128/jvi.00708-22] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/04/2022] [Accepted: 06/06/2022] [Indexed: 12/30/2022] Open
Abstract
Evolutionary analyses of viral sequences can provide insights into transmission dynamics, which in turn can optimize prevention interventions. Here, we characterized the dynamics of HIV transmission within the Mexico City metropolitan area. HIV pol sequences from persons recently diagnosed at the largest HIV clinic in Mexico City (between 2016 and 2021) were annotated with demographic/geographic metadata. A multistep phylogenetic approach was applied to identify putative transmission clades. A data set of publicly available sequences was used to assess international introductions. Clades were analyzed with a discrete phylogeographic model to evaluate the timing and intensity of HIV introductions and transmission dynamics among municipalities in the region. A total of 6,802 sequences across 96 municipalities (5,192 from Mexico City and 1,610 from the neighboring State of Mexico) were included (93.6% cisgender men, 5.0% cisgender women, and 1.3% transgender women); 3,971 of these sequences formed 1,206 clusters, involving 78 municipalities, including 89 clusters of ≥10 sequences. Discrete phylogeographic analysis revealed (i) 1,032 viral introductions into the region, over one-half of which were from the United States, and (ii) 354 migration events between municipalities with high support (adjusted Bayes factor of ≥3). The most frequent viral migrations occurred between northern municipalities within Mexico City, i.e., Cuauhtémoc to Iztapalapa (5.2% of events), Iztapalapa to Gustavo A. Madero (5.4%), and Gustavo A. Madero to Cuauhtémoc (6.5%). Our analysis illustrates the complexity of HIV transmission within the Mexico City metropolitan area but also identifies a spatially active transmission area involving a few municipalities in the north of the city, where targeted interventions could have a more pronounced effect on the entire regional epidemic. IMPORTANCE Phylogeographic investigation of the Mexico City HIV epidemic illustrates the complexity of HIV transmission in the region. An active transmission area involving a few municipalities in the north of the city, with transmission links throughout the region, is identified and could be a location where targeted interventions could have a more pronounced effect on the entire regional epidemic, compared with those dispersed in other manners.
Collapse
Affiliation(s)
- Santiago Avila-Rios
- Center for Research in Infectious Diseases, National Institute of Respiratory Diseases, Mexico City, Mexico
| | - Claudia García-Morales
- Center for Research in Infectious Diseases, National Institute of Respiratory Diseases, Mexico City, Mexico
| | - Gustavo Reyes-Terán
- Coordinating Commission of the National Institutes of Health and High Specialty Hospitals, Ministry of Health, Mexico City, Mexico
| | | | | | - Sanjay R. Mehta
- Division of Infectious Diseases and Global Public Health, University of California, San Diego, San Diego, California, USA
- Veterans Affairs Health System, San Diego, California, USA
| | - Antoine Chaillon
- Division of Infectious Diseases and Global Public Health, University of California, San Diego, San Diego, California, USA
| |
Collapse
|
4
|
Oliveira RC, Gräf T, Rego FFDA, Silva GPSA, Giovanetti M, Monteiro Cunha JP. Dynamic Dispersion of HIV-1 Subtype C Toward Brazilian Northeastern Region. AIDS Res Hum Retroviruses 2021; 37:913-921. [PMID: 34036794 DOI: 10.1089/aid.2020.0308] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/12/2022] Open
Abstract
The subtype C accounts for >50% of HIV type 1 (HIV-1) infections worldwide and it is currently the predominant viral form in South Brazil. Subtype C has been reported in all Brazilian regions; however, the phylogenetic relationship among strains circulating in those regions still remains unclear. This study aimed to investigate the origin and dynamic dispersion of HIV-1 subtype C toward Northeast Brazil. Our phylogenetic analysis suggests that most subtype C strains circulating in Brazil (99%) are descendant from the main lineage whose entrance in the country was previously described in the 1970s. According to the literature, additional introductions of subtype C were reported in the country through the Southeast region and in this study we identified another entry event that occurred most likely through the North region. Furthermore, our analysis suggests that the spread of subtype C to Brazilian Northeastern states occurred through multiple independent introductions of the main lineage that originated in South Brazil between mid-1980s and late 1990s. Despite the observation of eventual new HIV-1 subtype C introductions, our results highlight the predominance of a single lineage of this subtype in Brazil and the importance of South region in its dissemination throughout the country.
Collapse
Affiliation(s)
- Rodrigo Cunha Oliveira
- Núcleo de Bioinformática, Departamento de Bioquímica e Biofísica, Universidade Federal da Bahia, Salvador, Brazil
| | - Tiago Gräf
- Instituto Gonçalo Moniz, Fundação Oswaldo Cruz—FIOCRUZ, Salvador, Brazil
| | | | | | - Marta Giovanetti
- Instituto Oswaldo Cruz, Fundação Oswaldo Cruz—FIOCRUZ, Rio de Janeiro, Brazil
- Laboratório de Genética Celular e Molecular, Instituto de Ciências Biológicas, Universidade Federal de Minas Gerais, Belo Horizonte, Brazil
| | - Joana Paixão Monteiro Cunha
- Núcleo de Bioinformática, Departamento de Bioquímica e Biofísica, Universidade Federal da Bahia, Salvador, Brazil
| |
Collapse
|
5
|
An M, Zheng C, Li H, Chen L, Yang Z, Gan Y, Han X, Zhao J, Shang H. Independent epidemic patterns of HIV-1 CRF01_AE lineages driven by mobile population in Shenzhen, an immigrant city of China. Virus Evol 2021; 7:veab094. [PMID: 35299786 PMCID: PMC8923236 DOI: 10.1093/ve/veab094] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/05/2021] [Revised: 10/19/2021] [Accepted: 11/08/2021] [Indexed: 12/14/2022] Open
Abstract
Abstract
Shenzhen, a city with >12 million migrant population, may play a key role in the spread of human immunodeficiency virus (HIV)-1 in China. The transmission dynamics of CRF01_AE, a predominant subtype in Shenzhen, is a good model to characterize the impact of human mobility on HIV-1 epidemic locally and nationally. We used phylodynamic and phylogeographic methods to estimate the viral transmission dynamics and migration trajectory of variable lineages based on 1,423 CRF01_AE sequences in Shenzhen sampled between 2006 and 2015. Eleven lineages of CRF01_AE were detected in Shenzhen. Of those, four main lineages originated during the 1990s. Their basic viral reproduction number (R0) ranged 1.96–3.92. The effective viral reproduction number (Re) of two lineages prevalent among heterosexuals/people who inject drugs had reduced <1 at the end of sampling, and the main sources were the intra-provincial immigrants (72 per cent) for one and local residents of Shenzhen (91 per cent) for another. Within two lineages among men who have sex with men (MSM), Re had been above or close to 1 at the end of sampling, and the immigrants from Jiangxi/Shaanxi and Hubei as sources accounted for 93 per cent and 68 per cent of all viral migration events, respectively. Moreover, no obvious recipients were found throughout the viral migration history for any lineage. Our findings demonstrate that HIV epidemic is declining in Shenzhen, which coincided with the initiation of the interventions during the 2000s. However, the obvious differences of the epidemic patterns between lineages emphasize the importance of further targeting interventions and continued molecular tracing, focusing on high-risk transmission sources among MSM.
Collapse
Affiliation(s)
- Minghui An
- NHC Key Laboratory of AIDS Immunology (China Medical University), National Clinical Research Center for Laboratory Medicine, The First Affiliated Hospital of China Medical University, No 155, Nanjing North Street, Shenyang, Liaoning 110001, China
- Key Laboratory of AIDS Immunology, Chinese Academy of Medical Sciences, No 155, Nanjing North Street, Shenyang, Liaoning 110001, China
| | - Chenli Zheng
- Shenzhen Center for Disease Control and Prevention, No 8, Longyuan Road, Shenzhen, Guangdong 518055, China
| | - Hao Li
- Shenzhen Center for Disease Control and Prevention, No 8, Longyuan Road, Shenzhen, Guangdong 518055, China
| | - Lin Chen
- Shenzhen Center for Disease Control and Prevention, No 8, Longyuan Road, Shenzhen, Guangdong 518055, China
| | - Zhengrong Yang
- Shenzhen Center for Disease Control and Prevention, No 8, Longyuan Road, Shenzhen, Guangdong 518055, China
| | - Yongxia Gan
- Shenzhen Center for Disease Control and Prevention, No 8, Longyuan Road, Shenzhen, Guangdong 518055, China
| | - Xiaoxu Han
- NHC Key Laboratory of AIDS Immunology (China Medical University), National Clinical Research Center for Laboratory Medicine, The First Affiliated Hospital of China Medical University, No 155, Nanjing North Street, Shenyang, Liaoning 110001, China
- Key Laboratory of AIDS Immunology, Chinese Academy of Medical Sciences, No 155, Nanjing North Street, Shenyang, Liaoning 110001, China
| | - Jin Zhao
- Shenzhen Center for Disease Control and Prevention, No 8, Longyuan Road, Shenzhen, Guangdong 518055, China
| | - Hong Shang
- NHC Key Laboratory of AIDS Immunology (China Medical University), National Clinical Research Center for Laboratory Medicine, The First Affiliated Hospital of China Medical University, No 155, Nanjing North Street, Shenyang, Liaoning 110001, China
- Key Laboratory of AIDS Immunology, Chinese Academy of Medical Sciences, No 155, Nanjing North Street, Shenyang, Liaoning 110001, China
| |
Collapse
|
6
|
Fritsch HM, Tanuri A, Pinto AR, Gräf T. Transmission dynamics and molecular characterization of HIV-1 epidemic among therapeutic failure patients in Santa Catarina state, southern Brazil. INFECTION GENETICS AND EVOLUTION 2021; 92:104854. [PMID: 33839313 DOI: 10.1016/j.meegid.2021.104854] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/05/2021] [Revised: 03/28/2021] [Accepted: 04/05/2021] [Indexed: 10/21/2022]
Abstract
The HIV-1 epidemic in southern Brazil is mostly caused by subtype C, which contrasts the dominance of subtype B in the other regions of the country. Santa Catarina (SC), although the smallest state in the southern region, presents one of the highest incidences and mortality rates in Brazil due to AIDS. This work investigated the HIV-1 molecular diversity and phylogenetic transmission networks in SC state by analyzing a database of 3070 sequences of the national genotyping service. HIV-1C proved to be the most frequent subtype, with a significant increase in prevalence over time. HIV-1B was observed to be associated with highly educated men, suggesting a compartmentalization from other subtypes. Such observation was confirmed by the high frequency of HIV-1B circulating in MSM transmission networks. Identified transmission clusters were majority composed by individuals living up to 25 km away and interstate linkages were mainly between southern neighbor states. In general, individuals between 25 and 40 years old and sequences sampled after 2014 were more likely to be in transmission chains, in agreement with the universal treatment protocol launched in 2014. The present study brings new insights about HIV-1 transmission dynamics in southern Brazil.
Collapse
Affiliation(s)
- Hegger M Fritsch
- Laboratório de Imunologia Aplicada, Universidade Federal de Santa Catarina, Florianópolis, Santa Catarina, Brazil
| | - Amilcar Tanuri
- Laboratório de Virologia Molecular, Universidade Federal do Rio de Janeiro, Rio de Janeiro, Rio de Janeiro, Brazil
| | - Aguinaldo R Pinto
- Laboratório de Imunologia Aplicada, Universidade Federal de Santa Catarina, Florianópolis, Santa Catarina, Brazil.
| | - Tiago Gräf
- Instituto Gonçalo Moniz, Fundação Oswaldo Cruz, Salvador, Bahia, Brazil.
| |
Collapse
|
7
|
Chaillon A, Smith DM. Phylogenetic analyses of SARS-CoV-2 B.1.1.7 lineage suggest a single origin followed by multiple exportation events versus convergent evolution. Clin Infect Dis 2021; 73:2314-2317. [PMID: 33772259 PMCID: PMC8083653 DOI: 10.1093/cid/ciab265] [Citation(s) in RCA: 15] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/20/2021] [Indexed: 12/21/2022] Open
Abstract
The emergence of new variants of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) herald a new phase of the pandemic. This study used state-of-the-art phylodynamic methods to ascertain that the rapid rise of B.1.1.7 “Variant of Concern” most likely occurred by global dispersal rather than convergent evolution from multiple sources.
Collapse
Affiliation(s)
- A Chaillon
- Division of Infectious Diseases and Global Public Health, University of California San Diego, CA, USA
| | - D M Smith
- Division of Infectious Diseases and Global Public Health, University of California San Diego, CA, USA
| |
Collapse
|
8
|
Vrancken B, Zhao B, Li X, Han X, Liu H, Zhao J, Zhong P, Lin Y, Zai J, Liu M, Smith DM, Dellicour S, Chaillon A. Comparative Circulation Dynamics of the Five Main HIV Types in China. J Virol 2020; 94:e00683-20. [PMID: 32938762 PMCID: PMC7654276 DOI: 10.1128/jvi.00683-20] [Citation(s) in RCA: 17] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/14/2020] [Accepted: 09/02/2020] [Indexed: 01/17/2023] Open
Abstract
The HIV epidemic in China accounts for 3% of the global HIV incidence. We compared the patterns and determinants of interprovincial spread of the five most prevalent circulating types. HIV pol sequences sampled across China were used to identify relevant transmission networks of the five most relevant HIV-1 types (B and circulating recombinant forms [CRFs] CRF01_AE, CRF07_BC, CRF08_BC, and CRF55_01B) in China. From these, the dispersal history across provinces was inferred. A generalized linear model (GLM) was used to test the association between migration rates among provinces and several measures of human mobility. A total of 10,707 sequences were collected between 2004 and 2017 across 26 provinces, among which 1,962 are newly reported here. A mean of 18 (minimum and maximum, 1 and 54) independent transmission networks involving up to 17 provinces were identified. Discrete phylogeographic analysis largely recapitulates the documented spread of the HIV types, which in turn, mirrors within-China population migration flows to a large extent. In line with the different spatiotemporal spread dynamics, the identified drivers thereof were also heterogeneous but are consistent with a central role of human mobility. The comparative analysis of the dispersal dynamics of the five main HIV types circulating in China suggests a key role of large population centers and developed transportation infrastructures as hubs of HIV dispersal. This advocates for coordinated public health efforts in addition to local targeted interventions.IMPORTANCE While traditional epidemiological studies are of great interest in describing the dynamics of epidemics, they struggle to fully capture the geospatial dynamics and factors driving the dispersal of pathogens like HIV as they have difficulties capturing linkages between infections. To overcome this, we used a discrete phylogeographic approach coupled to a generalized linear model extension to characterize the dynamics and drivers of the across-province spread of the five main HIV types circulating in China. Our results indicate that large urbanized areas with dense populations and developed transportation infrastructures are facilitators of HIV dispersal throughout China and highlight the need to consider harmonized country-wide public policies to control local HIV epidemics.
Collapse
Affiliation(s)
- Bram Vrancken
- Department of Microbiology, Immunology and Transplantation, Rega Institute, Laboratory for Computational and Evolutionary Virology, KU Leuven, Leuven, Belgium
| | - Bin Zhao
- NHC Key Laboratory of AIDS Immunology (China Medical University), National Clinical Research Center for Laboratory Medicine, The First Affiliated Hospital of China Medical University, Shenyang, China
| | - Xingguang Li
- Department of Hospital Office, The First People's Hospital of Fangchenggang, Fangchenggang, China
| | - Xiaoxu Han
- NHC Key Laboratory of AIDS Immunology (China Medical University), National Clinical Research Center for Laboratory Medicine, The First Affiliated Hospital of China Medical University, Shenyang, China
| | - Haizhou Liu
- Centre for Emerging Infectious Diseases, The State Key Laboratory of Virology, Wuhan Institute of Virology, University of Chinese Academy of Sciences, Wuhan, China
| | - Jin Zhao
- Shenzhen Center for Disease Control and Prevention, Shenzhen, China
| | - Ping Zhong
- Department of AIDS and STD, Shanghai Municipal Center for Disease Control and Prevention; Shanghai Municipal Institutes for Preventive Medicine, Shanghai, China
| | - Yi Lin
- Department of AIDS and STD, Shanghai Municipal Center for Disease Control and Prevention; Shanghai Municipal Institutes for Preventive Medicine, Shanghai, China
| | - Junjie Zai
- Immunology innovation Team, School of Medicine, Ningbo University, Ningbo, Zhejiang China
| | - Mingchen Liu
- NHC Key Laboratory of AIDS Immunology (China Medical University), National Clinical Research Center for Laboratory Medicine, The First Affiliated Hospital of China Medical University, Shenyang, China
| | - Davey M Smith
- Division of Infectious Diseases and Global Public Health, University of California San Diego, California, USA
| | - Simon Dellicour
- Department of Microbiology, Immunology and Transplantation, Rega Institute, Laboratory for Computational and Evolutionary Virology, KU Leuven, Leuven, Belgium
- Spatial Epidemiology Lab (SpELL), Université Libre de Bruxelles, Brussels, Belgium
| | - Antoine Chaillon
- Division of Infectious Diseases and Global Public Health, University of California San Diego, California, USA
| |
Collapse
|
9
|
Vrancken B, Mehta SR, Ávila-Ríos S, García-Morales C, Tapia-Trejo D, Reyes-Terán G, Navarro-Álvarez S, Little SJ, Hoenigl M, Pines HA, Patterson T, Strathdee SA, Smith DM, Dellicour S, Chaillon A. Dynamics and Dispersal of Local HIV Epidemics Within San Diego and Across The San Diego-Tijuana Border. Clin Infect Dis 2020; 73:e2018-e2025. [PMID: 33079188 DOI: 10.1093/cid/ciaa1588] [Citation(s) in RCA: 16] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/15/2020] [Indexed: 01/25/2023] Open
Abstract
BACKGROUND Evolutionary analyses of well-annotated HIV sequence data can provide insights into viral transmission patterns and associated factors. Here, we explored the transmission dynamics of the HIV-1 subtype B epidemic across the San Diego (US) - Tijuana (Mexico) border region to identify factors that could help guide public health policy. METHODS HIV pol sequences were collected from people with HIV in San Diego County and from Tijuana between 1996-2018. A multistep phylogenetic approach was used to characterize the dynamics of spread. The contribution of geospatial factors and HIV risk group to the local dynamics were evaluated. RESULTS Phylogeographic analyses of the 2,034 sequences revealed an important contribution of local transmission in sustaining the epidemic, as well as a complex viral migration network across the region. Geospatial viral dispersal between San Diego communities occurred predominantly among men-who-have-sex with-men with central San Diego being the main source (34.9%) and recipient (39.5%) of migration events. HIV migration was more frequent from San Diego county towards Tijuana than vice versa. Migrations were best explained by driving time between locations. CONCLUSION The US-Mexico border may not be a major barrier to the spread of HIV, which may stimulate coordinated transnational intervention approaches. Whereas a focus on central San Diego has the potential to avert most spread, the substantial viral migration independent of central San Diego shows that county-wide efforts will be more effective. Combined, this work shows that epidemiological information gleaned from pathogen genomes can uncover mechanisms that underlie sustained spread and, in turn, can be a building block of public health decision making.
Collapse
Affiliation(s)
- Bram Vrancken
- Department of Microbiology, Immunology and Transplantation, Rega Institute, Laboratory for Computational and Evolutionary Virology, KU Leuven, Herestraat, Leuven, Belgium
| | - Sanjay R Mehta
- Division of Infectious Diseases and Global Public Health, University of California San Diego, CA
| | - Santiago Ávila-Ríos
- Centre for Research in Infectious Diseases, National Institute of Respiratory Diseases, Calzada de Tlalpan, Colonia Sección XVI, CP, Mexico City, Mexico
| | - Claudia García-Morales
- Centre for Research in Infectious Diseases, National Institute of Respiratory Diseases, Calzada de Tlalpan, Colonia Sección XVI, CP, Mexico City, Mexico
| | - Daniela Tapia-Trejo
- Centre for Research in Infectious Diseases, National Institute of Respiratory Diseases, Calzada de Tlalpan, Colonia Sección XVI, CP, Mexico City, Mexico
| | - Gustavo Reyes-Terán
- Coordinating Commission of the Mexican National Institutes of Health, Periférico Sur, Arenal Tepepan, Mexico City, Mexico
| | | | - Susan J Little
- Division of Infectious Diseases and Global Public Health, University of California San Diego, CA
| | - Martin Hoenigl
- Division of Infectious Diseases and Global Public Health, University of California San Diego, CA
| | - Heather A Pines
- Division of Infectious Diseases and Global Public Health, University of California San Diego, CA
| | - Thomas Patterson
- Division of Infectious Diseases and Global Public Health, University of California San Diego, CA
| | - Steffanie A Strathdee
- Division of Infectious Diseases and Global Public Health, University of California San Diego, CA
| | - Davey M Smith
- Division of Infectious Diseases and Global Public Health, University of California San Diego, CA
| | - Simon Dellicour
- Department of Microbiology, Immunology and Transplantation, Rega Institute, Laboratory for Computational and Evolutionary Virology, KU Leuven, Herestraat, Leuven, Belgium.,Spatial Epidemiology Lab (SpELL), Université Libre de Bruxelles, av. FD Roosevelt, Bruxelles, Belgium
| | - Antoine Chaillon
- Division of Infectious Diseases and Global Public Health, University of California San Diego, CA
| |
Collapse
|
10
|
Ishikawa SA, Zhukova A, Iwasaki W, Gascuel O. A Fast Likelihood Method to Reconstruct and Visualize Ancestral Scenarios. Mol Biol Evol 2019; 36:2069-2085. [PMID: 31127303 PMCID: PMC6735705 DOI: 10.1093/molbev/msz131] [Citation(s) in RCA: 117] [Impact Index Per Article: 23.4] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022] Open
Abstract
The reconstruction of ancestral scenarios is widely used to study the evolution of characters along phylogenetic trees. One commonly uses the marginal posterior probabilities of the character states, or the joint reconstruction of the most likely scenario. However, marginal reconstructions provide users with state probabilities, which are difficult to interpret and visualize, whereas joint reconstructions select a unique state for every tree node and thus do not reflect the uncertainty of inferences. We propose a simple and fast approach, which is in between these two extremes. We use decision-theory concepts (namely, the Brier score) to associate each node in the tree to a set of likely states. A unique state is predicted in tree regions with low uncertainty, whereas several states are predicted in uncertain regions, typically around the tree root. To visualize the results, we cluster the neighboring nodes associated with the same states and use graph visualization tools. The method is implemented in the PastML program and web server. The results on simulated data demonstrate the accuracy and robustness of the approach. PastML was applied to the phylogeography of Dengue serotype 2 (DENV2), and the evolution of drug resistances in a large HIV data set. These analyses took a few minutes and provided convincing results. PastML retrieved the main transmission routes of human DENV2 and showed the uncertainty of the human-sylvatic DENV2 geographic origin. With HIV, the results show that resistance mutations mostly emerge independently under treatment pressure, but resistance clusters are found, corresponding to transmissions among untreated patients.
Collapse
Affiliation(s)
- Sohta A Ishikawa
- Unité Bioinformatique Evolutive, Institut Pasteur, C3BI USR 3756 IP & CNRS, Paris, France
- Department of Biological Sciences, The University of Tokyo, Tokyo, Japan
- Evolutionary Genomics of RNA Viruses, Virology Department, Institut Pasteur, Paris, France
| | - Anna Zhukova
- Unité Bioinformatique Evolutive, Institut Pasteur, C3BI USR 3756 IP & CNRS, Paris, France
| | - Wataru Iwasaki
- Department of Biological Sciences, The University of Tokyo, Tokyo, Japan
| | - Olivier Gascuel
- Unité Bioinformatique Evolutive, Institut Pasteur, C3BI USR 3756 IP & CNRS, Paris, France
| |
Collapse
|
11
|
Bello G, Arantes I, Lacoste V, Ouka M, Boncy J, Césaire R, Liautaud B, Nacher M, Dos Santos G. Phylogeographic Analyses Reveal the Early Expansion and Frequent Bidirectional Cross-Border Transmissions of Non-pandemic HIV-1 Subtype B Strains in Hispaniola. Front Microbiol 2019; 10:1340. [PMID: 31333594 PMCID: PMC6622406 DOI: 10.3389/fmicb.2019.01340] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/27/2018] [Accepted: 05/29/2019] [Indexed: 11/13/2022] Open
Abstract
The human immunodeficiency virus-type 1 (HIV-1) subtype B has probably been circulating on the island of Hispaniola since the 1960s, but information about the early viral history on this Caribbean island is scarce. In this study, we reconstruct the dissemination dynamics of early divergent non-pandemic subtype B lineages (designated BCAR) on Hispaniola by analyzing a country-balanced dataset of HIV-1 BCARpol sequences from Haiti (n = 103) and the Dominican Republic (n = 123). Phylogenetic analyses supported that BCAR strains from Haiti and the Dominican Republic were highly intermixed between each other, although the null hypothesis of completely random mixing was rejected. Bayesian phylogeographic analyses placed the ancestral BCAR virus in Haiti and the Dominican Republic with the same posterior probability support. These analyses estimate frequent viral transmissions between Haiti and the Dominican Republic since the early 1970s onwards, and the presence of local BCAR transmission networks in both countries before first AIDS cases was officially recognized. Demographic reconstructions point that the BCAR epidemic in Hispaniola grew exponentially until the 1990s. These findings support that the HIV-1 epidemics in Haiti and the Dominican Republic have been connected by a recurrent bidirectional viral flux since the initial phase, which poses a great challenge in tracing the geographic origin of the BCAR epidemic within Hispaniola using only genetic data. These data also reinforce the notion that prevention programs have successfully reduced the rate of new HIV-1 transmissions in Hispaniola since the end of the 1990s.
Collapse
Affiliation(s)
- Gonzalo Bello
- Laboratório de AIDS e Imunologia Molecular, Instituto Oswaldo Cruz, FIOCRUZ, Rio de Janeiro, Brazil
| | - Ighor Arantes
- Laboratório de AIDS e Imunologia Molecular, Instituto Oswaldo Cruz, FIOCRUZ, Rio de Janeiro, Brazil
| | - Vincent Lacoste
- Laboratoire des Interactions Virus-Hôtes, Institut Pasteur de la Guyane, Cayenne, French Guiana
| | - Marlene Ouka
- Virology Laboratory, EA 4537, Martinique University Hospital, Fort de France, Martinique
| | - Jacques Boncy
- Laboratoire National de Santé Publique, Ministère de la Santé Publique et de la Population, Port-au-Prince, Haiti
| | - Raymond Césaire
- Virology Laboratory, EA 4537, Martinique University Hospital, Fort de France, Martinique
| | | | - Mathieu Nacher
- Coordination Régionale de la lutte contre le VIH (COREVIH) and Centre d'Investigation Clinique-CIC INSERM 1424, Centre Hospitalier de Cayenne "Andrée Rosemon", Cayenne, French Guiana
| | - Georges Dos Santos
- Virology Laboratory, EA 4537, Martinique University Hospital, Fort de France, Martinique
| |
Collapse
|
12
|
Alves BM, Siqueira JD, Prellwitz IM, Botelho OM, Da Hora VP, Sanabani S, Recordon-Pinson P, Fleury H, Soares EA, Soares MA. Estimating HIV-1 Genetic Diversity in Brazil Through Next-Generation Sequencing. Front Microbiol 2019; 10:749. [PMID: 31024510 PMCID: PMC6465556 DOI: 10.3389/fmicb.2019.00749] [Citation(s) in RCA: 16] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/15/2018] [Accepted: 03/25/2019] [Indexed: 12/26/2022] Open
Abstract
Approximately 36.7 million people were living with the human immunodeficiency virus (HIV) at the end of 2016 according to UNAIDS, representing a global prevalence rate of 0.8%. In Brazil, an HIV prevalence of 0.24% has been estimated, which represents approximately 830,000 individuals living with the virus. As a touristic and commercial hub in Latin America, Brazil harbors an elevated HIV genetic variability, further contributed by the selective pressure exerted by the host immune system and by antiretroviral treatment. Through the progress of the next-generation sequencing (NGS) techniques, it has been possible to expand the study of HIV genetic diversity, evolutionary, and epidemic processes, allowing the generation of HIV complete or near full-length genomes (NFLG) and improving the characterization of intra- and interhost diversity of viral populations. Greater sensitivity in the detection of viral recombinant forms represents one of the major improvements associated with this development. It is possible to identify unique or circulating recombinant forms using the near full-length viral genomes with increasing accuracy. It also permits the characterization of multiple viral infections within individual hosts. Previous Brazilian studies using NGS to analyze HIV diversity were able to identify several distinct unique and circulating recombinant forms and evidenced dual infections. These data unveiled unprecedented high rates of viral recombination and highlighted that novel recombinants are continually arising in the Brazilian epidemic. In the pooled analysis depicted in this report, HIV subtypes have been determined from HIV-positive patients in five states of Brazil with some of the highest HIV prevalence, three in the Southeast (Rio de Janeiro, São Paulo, and Minas Gerais), one in the Northeast (Pernambuco) and one in the South (Rio Grande do Sul). Combined data analysis showed a significant prevalence of recombinant forms (29%; 101/350), and a similar 26% when only NFLGs were considered. Moreover, the analysis was able to evidence the occurrence of multiple infections in some individuals. Our data highlight the great HIV genetic diversity found in Brazil and unveils a more accurate scenario of the HIV evolutionary dynamics in the region.
Collapse
Affiliation(s)
- Brunna M Alves
- Programa de Oncovirologia, Instituto Nacional de Câncer, Rio de Janeiro, Brazil
| | - Juliana D Siqueira
- Programa de Oncovirologia, Instituto Nacional de Câncer, Rio de Janeiro, Brazil
| | - Isabel M Prellwitz
- Programa de Oncovirologia, Instituto Nacional de Câncer, Rio de Janeiro, Brazil
| | - Ornella M Botelho
- Programa de Oncovirologia, Instituto Nacional de Câncer, Rio de Janeiro, Brazil
| | - Vanusa P Da Hora
- Laboratório de Biologia Molecular, Escola de Medicina, Universidade Federal do Rio Grande, Rio Grande do Sul, Brazil
| | - Sabri Sanabani
- LIM-3, Hospital das Clinicas FMUSP, Faculty of Medicine, University of São Paulo, São Paulo, Brazil
| | | | - Hervé Fleury
- CNRS MFP-UMR 5234, University Hospital of Bordeaux, University of Bordeaux, Bordeaux, France
| | - Esmeralda A Soares
- Programa de Oncovirologia, Instituto Nacional de Câncer, Rio de Janeiro, Brazil
| | - Marcelo A Soares
- Programa de Oncovirologia, Instituto Nacional de Câncer, Rio de Janeiro, Brazil.,Departamento de Genética, Universidade Federal do Rio de Janeiro, Rio de Janeiro, Brazil
| |
Collapse
|
13
|
Pérez AB, Vrancken B, Chueca N, Aguilera A, Reina G, García-del Toro M, Vera F, Von Wichman MA, Arenas JI, Téllez F, Pineda JA, Omar M, Bernal E, Rivero-Juárez A, Fernández-Fuertes E, de la Iglesia A, Pascasio JM, Lemey P, Garcia F, Cuypers L. Increasing importance of European lineages in seeding the hepatitis C virus subtype 1a epidemic in Spain. Euro Surveill 2019; 24:1800227. [PMID: 30862327 PMCID: PMC6402173 DOI: 10.2807/1560-7917.es.2019.24.9.1800227] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/29/2022] Open
Abstract
BackgroundReducing the burden of the hepatitis C virus (HCV) requires large-scale deployment of intervention programmes, which can be informed by the dynamic pattern of HCV spread. In Spain, ongoing transmission of HCV is mostly fuelled by people who inject drugs (PWID) infected with subtype 1a (HCV1a).AimOur aim was to map how infections spread within and between populations, which could help formulate more effective intervention programmes to halt the HCV1a epidemic in Spain.MethodsEpidemiological links between HCV1a viruses from a convenience sample of 283 patients in Spain, mostly PWID, collected between 2014 and 2016, and 1,317, 1,291 and 1,009 samples collected abroad between 1989 and 2016 were reconstructed using sequences covering the NS3, NS5A and NS5B genes. To efficiently do so, fast maximum likelihood-based tree estimation was coupled to a flexible Bayesian discrete phylogeographic inference method.ResultsThe transmission network structure of the Spanish HCV1a epidemic was shaped by continuous seeding of HCV1a into Spain, almost exclusively from North America and European countries. The latter became increasingly relevant and have dominated in recent times. Export from Spain to other countries in Europe was also strongly supported, although Spain was a net sink for European HCV1a lineages. Spatial reconstructions showed that the epidemic in Spain is diffuse, without large, dominant within-country networks.ConclusionTo boost the effectiveness of local intervention efforts, concerted supra-national strategies to control HCV1a transmission are needed, with a strong focus on the most important drivers of ongoing transmission, i.e. PWID and other high-risk populations.
Collapse
Affiliation(s)
- Ana Belen Pérez
- Department of Microbiology, Institute of Bio Sanitary Research (IBIS), AIDS Research Network, University Hospital of Granada, Granada, Spain,These authors contributed equally to the article
| | - Bram Vrancken
- These authors contributed equally to the article,KU Leuven, Department of Microbiology and Immunology, Rega Institute for Medical Research, Laboratory of Evolutionary and Computational Virology, Leuven, Belgium
| | - Natalia Chueca
- Department of Microbiology, Institute of Bio Sanitary Research (IBIS), AIDS Research Network, University Hospital of Granada, Granada, Spain
| | - Antonio Aguilera
- Department of Microbiology, University Hospital of Santiago, Santiago de Compostela, Spain
| | - Gabriel Reina
- Department of Microbiology, University Hospital of Navarra, Institute for Health Research (IdisNA), Pamplona, Spain
| | | | - Francisco Vera
- Unit of Infectious Diseases, Internal Medicine, General Hospital of Rosell, Cartagena, Murcia, Spain
| | | | - Juan Ignacio Arenas
- Unit of Infectious Diseases, Hospital Universitario de San Sebastian, San Sebastian, Spain
| | - Francisco Téllez
- Unit of Infectious Diseases and Microbiology, University Hospital of Puerto Real, Cádiz, Spain
| | - Juan A Pineda
- Unit of Infectious Diseases, University Hospital of Valme, Sevilla, Spain (J.A. Pineda)
| | | | - Enrique Bernal
- Unit of Infectious Diseases, General University Hospital, Murcia, Spain
| | - Antonio Rivero-Juárez
- Unit of Infectious Diseases, University Hospital Reina Sofía of Córdoba, Maimonides Institute of Biomedical Research of Córdoba, University of Córdoba, Córdoba, Spain
| | | | | | - Juan Manuel Pascasio
- Clinical Management Unit of Digestive Diseases, University Hospital of Virgen del Rocío, Sevilla, Spain
| | - Philippe Lemey
- KU Leuven, Department of Microbiology and Immunology, Rega Institute for Medical Research, Laboratory of Evolutionary and Computational Virology, Leuven, Belgium
| | - Féderico Garcia
- Department of Microbiology, Institute of Bio Sanitary Research (IBIS), AIDS Research Network, University Hospital of Granada, Granada, Spain,These authors contributed equally to the article
| | - Lize Cuypers
- These authors contributed equally to the article,KU Leuven, Department of Microbiology and Immunology, Rega Institute for Medical Research, Laboratory of Clinical and Epidemiological Virology, Leuven, Belgium
| |
Collapse
|
14
|
Oliveira MI, Romão de Souza Junior V, Fernanda de Lacerda Vidal C, Sérgio Ramos de Araújo P. Virologic suppression in response to antiretroviral therapy despite extensive resistance within HIV-1 reverse transcriptase after the first virologic failure. BMC Infect Dis 2018; 18:514. [PMID: 30314470 PMCID: PMC6186096 DOI: 10.1186/s12879-018-3400-6] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/28/2018] [Accepted: 09/20/2018] [Indexed: 11/20/2022] Open
Abstract
Background Incomplete virologic suppression results in mutations associated with resistance and is a major obstacle to disease control. We analyzed the genotypic profiles of HIV-1 patients at the time of the first virologic failure and the response to a salvage regimen after 48 weeks. Methods This work was a cross-sectional, retrospective, analytical study based on data collected from medical records and genotyping tests between 2006 and 2016. The sample consisted of data on individuals living with HIV (PLWH) from three major reference centers. Results A total of 184 patients were included in the data analysis. Viral subtype B was the most common (81.3%) as well as M184 V/I (85.3%) and K103 codon mutations (65.8%). Forty-eight weeks after switching to a salvage regimen, 67.3% of patients achieved an undetectable viral load. Discussion The number of mutations associated with nucleos(t)ide reverse transcriptase inhibitors (NRTI(t)s) did not affect virologic suppression (9.3% for zero NRTI(t)-associated mutations vs 48.6% for 1–2 NRTI(t)-associated mutations vs 42.1% for ≥3 NRTI(t)-associated mutations, p = 0.179). An ARV time (the beginning of the first ARV regimen up to genotyping) of > 36 months was a protective factor for detectable viral load (PR = 0.60, 95% CI = 0.39–0.92, p = 0.020) and a risk factor for developing ≥3 NRTI(t)-associated mutations (PR = 2.43, 95% CI 1.38–4.28, p = 0.002). Conclusions We found that extensive resistance to NRTI(t)s at the time of the first virologic failure did not impact virologic suppression at 48 weeks after switching to a second-line therapy based on NRTI(t)s plus protease inhibitors.
Collapse
Affiliation(s)
- Marta Iglis Oliveira
- Programa de Pós-graduação em Ciências da Saúde, Universidade Federal de Pernambuco, Av. Prof. Moraes Rego 1235, Recife, 50670-901, Brazil
| | - Valter Romão de Souza Junior
- Faculdade de Medicina do Recife, Universidade Federal de Pernambuco, Av. Prof. Moraes Rego, 1235, Recife, Pernambuco, Brazil.
| | | | - Paulo Sérgio Ramos de Araújo
- Programa de Pós-graduação em Ciências da Saúde, Universidade Federal de Pernambuco, Av. Prof. Moraes Rego 1235, Recife, 50670-901, Brazil.,Instituto Aggeu Magalhaes, FIOCRUZ, Av. Prof. Moraes Rego 1235, Recife, 50670-901, Brazil
| |
Collapse
|
15
|
Dellicour S, Vrancken B, Trovão NS, Fargette D, Lemey P. On the importance of negative controls in viral landscape phylogeography. Virus Evol 2018; 4:vey023. [PMID: 30151241 PMCID: PMC6101606 DOI: 10.1093/ve/vey023] [Citation(s) in RCA: 21] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Phylogeographic reconstructions are becoming an established procedure to evaluate the factors that could impact virus spread. While a discrete phylogeographic approach can be used to test predictors of transition rates among discrete locations, alternative continuous phylogeographic reconstructions can also be exploited to investigate the impact of underlying environmental layers on the dispersal velocity of a virus. The two approaches are complementary tools for studying pathogens' spread, but in both cases, care must be taken to avoid misinterpretations. Here, we analyse rice yellow mottle virus (RYMV) sequence data from West and East Africa to illustrate how both approaches can be used to study the impact of environmental factors on the virus’ dispersal frequency and velocity. While it was previously reported that host connectivity was a major determinant of RYMV spread, we show that this was a false positive result due to the lack of appropriate negative controls. We also discuss and compare the phylodynamic tools currently available for investigating the impact of environmental factors on virus spread.
Collapse
Affiliation(s)
- Simon Dellicour
- Laboratory for Clinical and Epidemiological Virology, Rega Institute, KU Leuven, Leuven, Belgium.,Spatial Epidemiology Lab (SpELL), Université Libre de Bruxelles, CP160/12 50, av. FD Roosevelt, 1050 Bruxelles, Belgium
| | - Bram Vrancken
- Laboratory for Clinical and Epidemiological Virology, Rega Institute, KU Leuven, Leuven, Belgium
| | - Nídia S Trovão
- Laboratory for Clinical and Epidemiological Virology, Rega Institute, KU Leuven, Leuven, Belgium
| | - Denis Fargette
- Institut de Recherche pour le Développement (IRD), UMR IPME (IRD, CIRAD, Université de Montpellier), BP 64051 34394 Montpellier cedex 5, France
| | - Philippe Lemey
- Laboratory for Clinical and Epidemiological Virology, Rega Institute, KU Leuven, Leuven, Belgium
| |
Collapse
|
16
|
Magee D, Scotch M. The effects of random taxa sampling schemes in Bayesian virus phylogeography. INFECTION GENETICS AND EVOLUTION 2018; 64:225-230. [PMID: 29991455 DOI: 10.1016/j.meegid.2018.07.003] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/02/2017] [Revised: 06/07/2018] [Accepted: 07/02/2018] [Indexed: 11/30/2022]
Abstract
Public health researchers are often tasked with accurately and quickly identifying the location and time when an epidemic originated from a representative sample of nucleotide sequences. In this paper, we investigate multiple approaches to subsampling the sequence set when employing a Bayesian phylogeographic generalized linear model. Our results indicate that near-categorical posterior MCC estimates on the root can be obtained with replicate runs using 25-50% of the sequence data, and that including 90% of sequences does not necessarily entail more accurate inferences. We present the first analysis of predictor signal suppression and show how the ability to detect the influence of predictor variables is limited when sample size predictors are included in the models.
Collapse
Affiliation(s)
- Daniel Magee
- Department of Biomedical Informatics, Arizona State University, 13212 E. Shea Blvd., Scottsdale 85259, AZ, USA; Biodesign Center for Environmental Health Engineering, Arizona State University, 1001 S. McAllister Ave, Tempe 85281, AZ, USA
| | - Matthew Scotch
- Department of Biomedical Informatics, Arizona State University, 13212 E. Shea Blvd., Scottsdale 85259, AZ, USA; Biodesign Center for Environmental Health Engineering, Arizona State University, 1001 S. McAllister Ave, Tempe 85281, AZ, USA.
| |
Collapse
|
17
|
Mir D, Gräf T, Esteves de Matos Almeida S, Pinto AR, Delatorre E, Bello G. Inferring population dynamics of HIV-1 subtype C epidemics in Eastern Africa and Southern Brazil applying different Bayesian phylodynamics approaches. Sci Rep 2018; 8:8778. [PMID: 29884822 PMCID: PMC5993807 DOI: 10.1038/s41598-018-26824-4] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/02/2018] [Accepted: 05/18/2018] [Indexed: 12/01/2022] Open
Abstract
The subtype C Eastern Africa clade (CEA), a particularly successful HIV-1 subtype C lineage, has seeded several sub-epidemics in Eastern African countries and Southern Brazil during the 1960s and 1970s. Here, we characterized the past population dynamics of the major CEA sub-epidemics in Eastern Africa and Brazil by using Bayesian phylodynamic approaches based on coalescent and birth-death models. All phylodynamic models support similar epidemic dynamics and exponential growth rates until roughly the mid-1980s for all the CEA sub-epidemics. Divergent growth patterns, however, were supported afterwards. The Bayesian skygrid coalescent model (BSKG) and the birth-death skyline model (BDSKY) supported longer exponential growth phases than the Bayesian skyline coalescent model (BSKL). The BDSKY model uncovers patterns of a recent decline for the CEA sub-epidemics in Burundi/Rwanda and Tanzania (Re < 1) and a recent growth for Southern Brazil (Re > 1); whereas coalescent models infer an epidemic stabilization. To the contrary, the BSKG model captured a decline of Ethiopian CEA sub-epidemic between the mid-1990s and mid-2000s that was not uncovered by the BDSKY model. These results underscore that the joint use of different phylodynamic approaches may yield complementary insights into the past HIV population dynamics.
Collapse
Affiliation(s)
- Daiana Mir
- Laboratório de AIDS e Imunologia Molecular, Instituto Oswaldo Cruz, Fiocruz, Rio de Janeiro, Brazil.
| | - Tiago Gräf
- Departamento de Genética, Instituto de Biologia, Universidade Federal do Rio de Janeiro, Rio de Janeiro, Brazil.,KwaZulu-Natal Research Innovation and Sequencing Platform (KRISP), College of Health Sciences, University of KwaZulu-Natal, Durban, South Africa
| | - Sabrina Esteves de Matos Almeida
- Centro de Desenvolvimento Científico e Tecnológico, Fundação Estadual de Produção e Pesquisa em Saúde, Porto Alegre, Brazil.,Programa de Pós-Graduação em Genética e Biologia Molecular, Universidade Federal do Rio Grande do Sul, Porto Alegre, Brazil.,Instituto de Ciências da Saúde, Universidade Feevale, Novo Hamburgo, Brazil
| | - Aguinaldo Roberto Pinto
- Laboratório de Imunologia Aplicada, Departamento de Microbiologia, Imunologia e Parasitologia, Universidade Federal de Santa Catarina, Florianópolis, Brazil
| | - Edson Delatorre
- Laboratório de AIDS e Imunologia Molecular, Instituto Oswaldo Cruz, Fiocruz, Rio de Janeiro, Brazil
| | - Gonzalo Bello
- Laboratório de AIDS e Imunologia Molecular, Instituto Oswaldo Cruz, Fiocruz, Rio de Janeiro, Brazil
| |
Collapse
|
18
|
The Effects of Sampling Location and Predictor Point Estimate Certainty on Posterior Support in Bayesian Phylogeographic Generalized Linear Models. Sci Rep 2018; 8:5905. [PMID: 29651124 PMCID: PMC5897398 DOI: 10.1038/s41598-018-24264-8] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/01/2017] [Accepted: 03/26/2018] [Indexed: 01/27/2023] Open
Abstract
The use of generalized linear models in Bayesian phylogeography has enabled researchers to simultaneously reconstruct the spatiotemporal history of a virus and quantify the contribution of predictor variables to that process. However, little is known about the sensitivity of this method to the choice of the discrete state partition. Here we investigate this question by analyzing a data set containing 299 sequences of the West Nile virus envelope gene sampled in the United States and fifteen predictors aggregated at four spatial levels. We demonstrate that although the topology of the viral phylogenies was consistent across analyses, support for the predictors depended on the level of aggregation. In particular, we found that the variance of the predictor support metrics was minimized at the most precise level for several predictors and maximized at more sparse levels of aggregation. These results suggest that caution should be taken when partitioning a region into discrete locations to ensure that interpretable, reproducible posterior estimates are obtained. These results also demonstrate why researchers should use the most precise discrete states possible to minimize the posterior variance in such estimates and reveal what truly drives the diffusion of viruses.
Collapse
|
19
|
Fritsch HM, Almeida SEM, Pinto AR, Gräf T. Spatiotemporal and demographic history of the HIV-1 circulating recombinant form CRF31_BC in Brazil. INFECTION GENETICS AND EVOLUTION 2018; 61:113-118. [PMID: 29601872 DOI: 10.1016/j.meegid.2018.03.024] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/10/2018] [Revised: 03/19/2018] [Accepted: 03/25/2018] [Indexed: 10/17/2022]
Abstract
CRF31_BC is an HIVs-1 recombinant form very prevalent in the southernmost capital city of Brazil, Porto Alegre. Recent studies have been describing a growing number of cases of infection by CRF31_BC in other Brazilian cities and countries, suggesting a process of expansion of this strain. Aiming to describe the city of origin, dispersion routes and demographic history of CRF31_BC, this study analyzed all HIV-1 CRF31_BC and Brazilian BC mosaic publicly available sequences. CRF31_BC classification was performed by bootscanning and tree reconstruction methods. Bayesian phylogeographic and phylodynamic model approaches were used to reconstruct the spatiotemporal and demographic history of 95 sequences identified as CRF31_BC-like. Porto Alegre was estimated to be the origin and center of the dispersion of the CRF31_BC for most of the analyzed locations. However, some viral transitions independent from Porto Alegre were observed in other cities from the Rio Grande do Sul state and also in other Brazilian states. The estimated CRF31_BC epidemic growth rate was similar to subtype C and B in Brazil. Our findings suggest that CRF31_BC, although mostly prevalent in south region, is circulating nation-wide with some localities presenting autochthonous transmissions.
Collapse
Affiliation(s)
- Hegger M Fritsch
- Universidade Federal do Rio Grande do Sul, Porto Alegre, RS, Brazil; Centro de Desenvolvimento Científico e Tecnológico, Secretaria de Saúde do Estado do Rio Grande do Sul, Porto Alegre, RS, Brazil
| | - Sabrina E M Almeida
- Centro de Desenvolvimento Científico e Tecnológico, Secretaria de Saúde do Estado do Rio Grande do Sul, Porto Alegre, RS, Brazil; Universidade Feevale, Novo Hamburgo, RS, Brazil
| | - Aguinaldo R Pinto
- Laboratório de Imunologia Aplicada, Departamento de Microbiologia, Imunologia e Parasitologia, Universidade Federal de Santa Catarina, Florianópolis, SC, Brazil
| | - Tiago Gräf
- KwaZulu-Natal Research Innovation and Sequencing Platform (KRISP), College of Health Sciences, University of KwaZulu-Natal, Durban, South Africa; Departamento de Genética, Instituto de Biologia, Universidade Federal do Rio de Janeiro, Rio de Janeiro, RJ, Brazil.
| |
Collapse
|
20
|
Delatorre E, Couto-Fernandez JC, Bello G. HIV-1 Genetic Diversity in Northeastern Brazil: High Prevalence of Non-B Subtypes. AIDS Res Hum Retroviruses 2017; 33:639-647. [PMID: 28325060 DOI: 10.1089/aid.2017.0045] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/16/2023] Open
Abstract
The Northeastern Brazilian region has experienced a constant increase in the number of newly reported AIDS cases over the last decade, but the genetic diversity of HIV-1 strains currently disseminated in this region remains poorly explored. HIV-1 pol sequences were obtained from 140 patients followed at outpatient clinics from four Northeastern Brazilian states (Alagoas, Bahia, Ceará, and Piauí) between 2014 and 2015. Subtype B was the most prevalent HIV-1 clade (72%) detected in the Northeastern region, followed by subtypes F1 (6%), C (5%), and D (1%). The remaining strains (16%) displayed a recombinant structure and were classified as follows: BF1 (11%), BC (4%), BCF1 (1%), and CRF02_AG like (1%). The 20 HIV-1 BF1 and BC recombinant sequences detected were distributed among 11 lineages classified as follows: CRF28/29_BF like (n = 5), CRF39_BF like (n = 1), URF_BF (n = 9), and URF_BC (n = 5). Non-B subtypes were detected in all Northeastern Brazilian states, but with variable prevalence, ranging from 16% in Ceará to 55% in Alagoas. Phylogenetic analyses support that subtype D and CRF02_AG strains detected in the Northeastern region resulted from the expansion of autochthonous transmission networks, rather than from exogenous introductions from other countries. These results reveal that HIV-1 epidemic spreading in the Northeastern Brazilian region comprised by multiple subtypes and recombinant strains and the molecular epidemiologic pattern in this Brazilian region is much more complex than originally estimated.
Collapse
Affiliation(s)
- Edson Delatorre
- Laboratório de AIDS & Imunologia Molecular, Instituto Oswaldo Cruz, FIOCRUZ, Rio de Janeiro, Brazil
| | | | - Gonzalo Bello
- Laboratório de AIDS & Imunologia Molecular, Instituto Oswaldo Cruz, FIOCRUZ, Rio de Janeiro, Brazil
| |
Collapse
|
21
|
The multi-faceted dynamics of HIV-1 transmission in Northern Alberta: A combined analysis of virus genetic and public health data. INFECTION GENETICS AND EVOLUTION 2017; 52:100-105. [PMID: 28427935 DOI: 10.1016/j.meegid.2017.04.005] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/18/2016] [Revised: 03/30/2017] [Accepted: 04/05/2017] [Indexed: 01/05/2023]
Abstract
Molecular epidemiology has become a key tool for tracking infectious disease epidemics. Here, the spread of the most prevalent HIV-1 subtypes in Northern Alberta, Canada, was characterized with a Bayesian phylogenetic approach using 1146 HIV-1 pol sequences collected between 2007 and 2013 for routine clinical management purposes. Available patient metadata were qualitatively interpreted and correlated with onwards transmission using Fisher exact tests and logistic regression. Most infections were from subtypes A (n=36), B (n=815) and C (n=211). Africa is the dominant origin location for subtypes A and C while the subtype B epidemic was seeded from the USA and Middle America and, from the early 1990s onwards, mostly by interprovincial spread. Subtypes A (77.8%) and C (74.0%) were usually heterosexually transmitted and circulate predominantly among Blacks (61.1% and 85% respectively). Subtype B was mostly found among Caucasians (48.6%) and First Nations (36.8%), and its modes of transmission were stratified by ethnic origin. Compared to subtypes A (5.6%) and C (3.8-10.0%), a larger portion of subtype B patients were found within putative provincial transmission networks (20.3-29.5%), and this almost doubled when focusing on nationwide transmission clusters (37.9-57.5%). No clear association between cluster membership and particular patient characteristics was found. This study reveals complex and multi-faceted transmission dynamics of the HIV-1 epidemic in this otherwise low HIV prevalence population in Northern Alberta, Canada. These findings can aid public health planning.
Collapse
|
22
|
Magee D, Suchard MA, Scotch M. Bayesian phylogeography of influenza A/H3N2 for the 2014-15 season in the United States using three frameworks of ancestral state reconstruction. PLoS Comput Biol 2017; 13:e1005389. [PMID: 28170397 PMCID: PMC5321473 DOI: 10.1371/journal.pcbi.1005389] [Citation(s) in RCA: 22] [Impact Index Per Article: 3.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/06/2016] [Revised: 02/22/2017] [Accepted: 01/27/2017] [Indexed: 12/17/2022] Open
Abstract
Ancestral state reconstructions in Bayesian phylogeography of virus pandemics have been improved by utilizing a Bayesian stochastic search variable selection (BSSVS) framework. Recently, this framework has been extended to model the transition rate matrix between discrete states as a generalized linear model (GLM) of genetic, geographic, demographic, and environmental predictors of interest to the virus and incorporating BSSVS to estimate the posterior inclusion probabilities of each predictor. Although the latter appears to enhance the biological validity of ancestral state reconstruction, there has yet to be a comparison of phylogenies created by the two methods. In this paper, we compare these two methods, while also using a primitive method without BSSVS, and highlight the differences in phylogenies created by each. We test six coalescent priors and six random sequence samples of H3N2 influenza during the 2014–15 flu season in the U.S. We show that the GLMs yield significantly greater root state posterior probabilities than the two alternative methods under five of the six priors, and significantly greater Kullback-Leibler divergence values than the two alternative methods under all priors. Furthermore, the GLMs strongly implicate temperature and precipitation as driving forces of this flu season and nearly unanimously identified a single root state, which exhibits the most tropical climate during a typical flu season in the U.S. The GLM, however, appears to be highly susceptible to sampling bias compared with the other methods, which casts doubt on whether its reconstructions should be favored over those created by alternate methods. We report that a BSSVS approach with a Poisson prior demonstrates less bias toward sample size under certain conditions than the GLMs or primitive models, and believe that the connection between reconstruction method and sampling bias warrants further investigation. For the better part of the last decade, epidemiological researchers have employed a Bayesian framework to reconstruct phylogenetic trees and determine the spatiotemporal relationships between clades of viruses. Recently, an extension of this framework has enabled direct assessment of how various demographic, geographic, genetic, and environmental variables play a role in these relationships, but there has yet to be a comparison between the former and the latter. Here, we aim to assess the differences between the two reconstruction techniques, as well as an additional primitive method, using the 2014–15 influenza season in the U.S. as a case study under a variety of population growth scenarios. We highlight how the new method demonstrates significant increases in commonly-reported trends in phylogenies and that the method identifies climate predictors that appear to be consistent with known trends in seasonal trends in influenza. However, we found that this method appears to be the most heavily influenced by the locations at which the viruses were obtained. Our work offers valuable insight for researchers wishing to study the evolutionary history of viruses and also may prove useful in determining the correct method to choose for a given application of virus phylogeography.
Collapse
Affiliation(s)
- Daniel Magee
- Department of Biomedical Informatics, Arizona State University, Tempe, Arizona, United States of America
- Biodesign Center for Environmental Security, Arizona State University, Tempe, Arizona, United States of America
| | - Marc A. Suchard
- Departments of Biomathematics and Human Genetics, David Geffen School of Medicine, University of California, Los Angeles, California, United States of America
- Department of Biostatistics, School of Public Health, University of California, Los Angeles, California, United States of America
| | - Matthew Scotch
- Department of Biomedical Informatics, Arizona State University, Tempe, Arizona, United States of America
- Biodesign Center for Environmental Security, Arizona State University, Tempe, Arizona, United States of America
- * E-mail:
| |
Collapse
|
23
|
Worobey M, Watts TD, McKay RA, Suchard MA, Granade T, Teuwen DE, Koblin BA, Heneine W, Lemey P, Jaffe HW. 1970s and 'Patient 0' HIV-1 genomes illuminate early HIV/AIDS history in North America. Nature 2016; 539:98-101. [PMID: 27783600 PMCID: PMC5257289 DOI: 10.1038/nature19827] [Citation(s) in RCA: 101] [Impact Index Per Article: 12.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/05/2016] [Accepted: 09/07/2016] [Indexed: 11/08/2022]
Abstract
The emergence of HIV-1 group M subtype B in North American men who have sex with men was a key turning point in the HIV/AIDS pandemic. Phylogenetic studies have suggested cryptic subtype B circulation in the United States (US) throughout the 1970s and an even older presence in the Caribbean. However, these temporal and geographical inferences, based upon partial HIV-1 genomes that postdate the recognition of AIDS in 1981, remain contentious and the earliest movements of the virus within the US are unknown. We serologically screened >2,000 1970s serum samples and developed a highly sensitive approach for recovering viral RNA from degraded archival samples. Here, we report eight coding-complete genomes from US serum samples from 1978-1979-eight of the nine oldest HIV-1 group M genomes to date. This early, full-genome 'snapshot' reveals that the US HIV-1 epidemic exhibited extensive genetic diversity in the 1970s but also provides strong evidence for its emergence from a pre-existing Caribbean epidemic. Bayesian phylogenetic analyses estimate the jump to the US at around 1970 and place the ancestral US virus in New York City with 0.99 posterior probability support, strongly suggesting this was the crucial hub of early US HIV/AIDS diversification. Logistic growth coalescent models reveal epidemic doubling times of 0.86 and 1.12 years for the US and Caribbean, respectively, suggesting rapid early expansion in each location. Comparisons with more recent data reveal many of these insights to be unattainable without archival, full-genome sequences. We also recovered the HIV-1 genome from the individual known as 'Patient 0' (ref. 5) and found neither biological nor historical evidence that he was the primary case in the US or for subtype B as a whole. We discuss the genesis and persistence of this belief in the light of these evolutionary insights.
Collapse
Affiliation(s)
- Michael Worobey
- Department of Ecology and Evolutionary Biology, University of Arizona, Tucson, Arizona, 85721, USA
| | - Thomas D. Watts
- Department of Ecology and Evolutionary Biology, University of Arizona, Tucson, Arizona, 85721, USA
| | - Richard A. McKay
- Department of History and Philosophy of Science, University of Cambridge, Cambridge, Cambridgeshire, CB2 3RH, UK
| | - Marc A. Suchard
- Departments of Biomathematics, Biostatistics and Human Genetics, David Geffen School of Medicine, University of California Los Angeles, Los Angeles, CA, 90095, USA
| | - Timothy Granade
- Centers for Disease Control and Prevention, Atlanta, Georgia 30333, USA
| | | | - Beryl A. Koblin
- Laboratory of Infectious Disease Prevention, The New York Blood Center, New York, New York 10065, USA
| | - Walid Heneine
- Centers for Disease Control and Prevention, Atlanta, Georgia 30333, USA
| | - Philippe Lemey
- Department of Microbiology and Immunology, Rega Institute, KU Leuven – University of Leuven, Minderbroedersstaat 10, 3000 Leuven, Belgium
| | - Harold W. Jaffe
- Centers for Disease Control and Prevention, Atlanta, Georgia 30333, USA
| |
Collapse
|
24
|
Comprehensive Characterization of HIV-1 Molecular Epidemiology and Demographic History in the Brazilian Region Most Heavily Affected by AIDS. J Virol 2016; 90:8160-8. [PMID: 27384663 DOI: 10.1128/jvi.00363-16] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/21/2016] [Accepted: 06/27/2016] [Indexed: 11/20/2022] Open
Abstract
UNLABELLED The high incidence of AIDS cases and the dominance of HIV-1 subtype C infections are two features that distinguish the HIV-1 epidemic in the two southernmost Brazilian states (Rio Grande do Sul [RS] and Santa Catarina [SC]) from the epidemic in other parts of the country. Nevertheless, previous studies on HIV molecular epidemiology were conducted mainly in capital cities, and a more comprehensive understanding of factors driving this unique epidemic in Brazil is necessary. Blood samples were collected from individuals in 13 municipalities in the Brazilian southern region. HIV-1 env and pol genes were submitted to phylogenetic analyses for assignment of subtype, and viral population phylodynamics were reconstructed by applying Skygrid and logistic coalescent models in a Bayesian analysis. A high prevalence of subtype C was observed in all sampled locations; however, an increased frequency of recombinant strains was found in RS, with evidence for new circulating forms (CRFs). In the SC state, subtype B and C epidemics were associated with distinct exposure groups. Although logistic models estimated similar growth rates for HIV-1 subtype C (HIV-1C) and HIV-1B, a Skygrid plot reveals that the former epidemic has been expanding for a longer time. Our results highlight a consistent expansion of HIV-1C in south Brazil, and we also discuss how heterosexual and men who have sex with men (MSM) transmission chains might have impacted the current prevalence of HIV-1 subtypes in this region. IMPORTANCE The AIDS epidemic in south Brazil is expanding rapidly, but the circumstances driving this condition are not well known. A high prevalence of HIV-1 subtype C was reported in the capital cities of this region, in contrast to the subtype B dominance in the rest of the country. This study sought to comparatively investigate the HIV-1 subtype B and C epidemics by sampling individuals from several cities in the two states with the highest AIDS incidences in Brazil. Our analyses showed distinct epidemic growth curves for the two epidemics, and we also found evidence suggesting that separate transmission chains may be impacting the viral phylodynamics and the emergence of new recombinant forms.
Collapse
|