- Research article
- Open Access
The origin and spread of CRF85_BC, driven by heterosexual transmission among older people in Sichuan, China
BMC Infectious Diseases volume 20, Article number: 772 (2020)
CRF_BC recombinants, including CRF07_BC and CRF08_BC, were considered the predominant subtypes in China. Since the discovery of HIV-1 circulating recombinant form CRF 85_BC in Southwest China in 2016, this BC recombinant forms had been reported in different regions of China. However, the history and magnitude of CRF85_BC transmission were still to be investigated.
We conducted the most recent molecular epidemiology of HIV-1 among newly reported HIV-1 infected patients in Sichuan in 2019 by sequencing and phylogenetic analysis of 1291 pol sequences. Then, we used maximum likelihood approach and the Bayesian Markov chain Monte Carlo (MCMC) sampling of pol sequences to reconstruct the phylogeographic and demographic dynamics of the CRF85_BC.
HIV-1 CRF85_BC (68/1291, 5.27%) became the fourth most prevalent strain revealing a significant increase in local population. CRF85_BC were only found in heterosexually infected individuals and the majority of CRF85_BC (95.45%) were circulating among the people living with HIV aged 50 years and over (PLHIV50+), suggesting a unique prevalent pattern. The founder lineages of CRF85_BC were likely to have first emerged in Yunnan, a province of Southwest China bordering Sichuan, in the early 2000s. It then spread exponentially to various places (including Guangxi, Sichuan, et al) and became endemic around 2008.6 (2006.7–2010.2) in Sichuan.
Taken together, our findings on HIV-1 subtype CRF85_BC infections provided new insights into the spread of this virus and extended the understanding of the HIV epidemic in China.
The human immunodeficiency virus type 1 (HIV-1) epidemic remains a major public health burden in China. By the end of 2018, it is estimated that there are approximately 1.25 [1.10–1.40] million persons living with HIV/AIDS in China. During 2018, 81 [60–105] thousand new HIV infections occurred .
In early 1990s, HIV-1 Thai subtypes B′ and C were co-circulated among injection drug users (IDUs)  in Yunnan province in southwestern China bordering the ‘Golden Triangle’, leading to active inter-subtype recombination which generated CRF07_BC, CRF08_BC, and a large number of unique recombinant forms [3, 4]. In the initial period, drug trafficking activities were considered to be responsible for the rapid spread of CRF07_BC and CRF08_BC across China . From 1978, as commercial sex activity increased across the country , the proportion of newly reported HIV infections/acquired immune deficiency syndrome (AIDS) patients through heterosexual transmission increased dramatically. Nowadays, CRF07_BC and CRF08_BC have become the most epidemic strains in China [7, 8] and drug-driven epidemic has been shifting to sexually transmission-driven epidemic . There were other B/C CRFs reported in the last few years, such as CRF57_BC (identified in Baoshan, Western Yunnan) , CRF62_BC and CRF64_BC (identified in Ruili/Dehong, Western Yunnan) . However, there was scarce evidence that such B/C recombinant forms had spread to other areas. In 2016, we reported a novel B/C recombinant form, CRF85_BC, among heterosexually infected patients in southern Sichuan (Yibin city) . CRF85_BC was then reported in Anhui , Yunnan  and Northern Sichuan  (Guangyuan City). Due to lack of B or C epidemic history in the local high-risk population, the origin and the spread of CRF85_BC was still to be investigated.
Sichuan province is in the interior of Southwest China, adjacent to Yunnan and Guizhou in the South and Tibet in the West. Sichuan has a population of approximately 85 million people, and is a developing area in China. Since the initial HIV epidemic was fueled by IDU in 1990s, the HIV epidemic in Sichuan had become one of the worst in China. By 2018, there were 132, 680 HIV infections reported in Sichuan province, which ranked the first among all provinces/autonomous regions in China . Several studies have focused on the HIV epidemic in IDUs [17, 18] and CRF07_BC [19, 20] in Liangshan prefecture in Sichuan, and reported that the infected IDUs in this region could serve as a source of transmission to other regions of China . However, limited attention has been paid to the recent shift of the HIV-related risk behaviors from IDU to sexual contact, and also to other circulating HIV strains in this area.
Genetic and temporal dynamic analyses have been widely used to reconstruct the history of the HIV epidemic and have provided important information to aid in the development of strategies for the prevention of HIV-1 transmission. Here, we reported an investigation of the temporal and spatial dynamics of HIV CRF85_BC to give more clues about transmission of HIV strains in China.
A molecular epidemiological surveillance of Sichuan conducted among all newly HIV-1 diagnosed in April 2019. A total of 1573 antiretroviral-naive patient samples were collected at the first clinical visit. EDTA blood was collected with informed consent and 1291 HIV-1 pol sequences covering 1, 060 base pairs (HXB2: 2, 253–3, 312) were obtained as described previously . Demographic data (age, gender & self-identified ethnicity) and relevant information associated with HIV-1 transmission were abstracted from the national surveillance database. This study was reviewed and approved by the ethics committees of Sichuan Center for Disease Prevention and Control.
HIV-1 subtyping and pairwise distance analysis
HIV subtypes were assessed by phylogenetic analysis. The HIV-1 pol sequences were aligned with reference sequences of various subtypes and circulating recombinant forms (CRFs) from Los Alamos HIV sequence database (LANL, http://www.hiv.lanl.gov/). Multiple alignments were made automatically using the Mega version 7.0  with minor manual adjustments. FastTree 2.3 was used to estimate an approximately-maximum likelihood phylogenetic tree for pol sequences using the GTR + G + I nucleotide substitution model . The phylogenetic tree’s reliability was determined with local support values based on the Shimodaira-Hasegawa (SH) test  and presented using FigTree v1.4.3 (http://tree.bio.ed.ac.uk/software/figtree/). For sequences of unknown HIV-1 subtype, the jumping profile hidden Markov model (jpHMM) (http://jphmm.gobics.de/submission_hiv) was applied to screen for recombination breakpoints.
Pairwise distances of CRF85_BC sequences collected from two cross-sectional surveys (in 2014 and 2019) were computed using the Tamura-Nei model in the MEGA v7.0 software  .
CRF85_BC sequences dataset
CRF85_BC HIV-1 pol sequences covering 1, 060 base pairs (HXB2: 2, 253–3, 312) were derived from multiple ways, form (i) molecular epidemiological surveillance of Sichuan in 2019(N = 68), form (ii) molecular epidemiological surveillance of Sichuan in 2014  (N = 14), from (iii) Drug resistance surveillance database in Sichuan CDC between 2012 and 2019 (N = 58), (iv) from other Molecular epidemiological investigation conducted in Sichuan  (N = 9), and (v) all other CRF85_BC (Yunnan: n = 2; CQ: n = 1) and the sequences (Yunnan: n = 7; Anhui: n = 2; Guangxi: n = 2; HN: n = 1) that are most closely related to each of the Sichuan CRF85_BC sequences with known sampling dates by BLAST search in the Los Alamos HIV database (LANL, http://www.hiv.lanl.gov/). All the sequences were confirmed as CRF85_BC by phylogenetic analysis mentioned above.
To avoid biases due to convergent evolution introduced by drug resistance selection rather than genealogical similarity, known surveillance drug resistance mutations (SDRM) have been excluded from the sequences prior to any analyses as described previously [26, 27].
Bayesian MCMC evolutionary analysis
Estimation of evolutionary rate and the time of the most recent common ancestor (tMRCA) for CRF85_BC were performed by BEAST v1.8 . A Bayesian skygrid model with Markov chain Monte Carlo (MCMC) inference under the relaxed lognormal molecular clock was selected as a reliable mode for this analysis. The MCMC chains were run 200 million times with 20,000 generations logged in, and the first 10% generations were discarded as burn-in. Bayesian MCMC output was analyzed using TRACER v1.5 , and all parameters were estimated from an ESS > 200. The trees were summarized in a target tree using the Tree Annotator program and scanned using the Fig. Tree 1.4.3.
Chi-square test was used to compare the distributions of the population and demographic information among various HIV-1 strains. P values less than 0.05 were considered as statistical significance. All statistical analyses were performed using SPSS v.20.0 software (IBM Company, New York, USA).
Molecular epidemiological surveillance of Sichuan in 2019
Of the 1291 available sequences in the molecular epidemiological surveillance of Sichuan in 2019, CRF07_BC (676/1291, 52.36%), CRF01_AE (341/1291, 26.41%), CRF08_BC (114/1291, 8.83%), and CRF85_BC (68/1291, 5.27%) were the 4 most epidemic strains (Fig. 1a). CRF85_BC dominated the HIV-1 epidemics in Yibin (58/146), a southeast city in Sichuan Province, and was also detected in Luzhou (4/174), Leshan (3/120), Chengdu (2/312) and Mianyang (1/70) (Fig. 1b). Compare with the data of Sichuan molecular epidemiological surveillance in 2014 , the proportion of CRF85_BC in Sichuan increased from 3.66% (14/383) to 5.27%, and it became more widely distributed. To further understand whether CRF85_BC was involved in ongoing transmission, we examined the sequence diversity by estimating pairwise distances between sequences in 2014(n = 14) and sequences in 2019(n = 68). We detected a much more diversity in 2019 (mean distance 0.023, SD 0.002) compared with 2014 (mean distance 0.010, SD 0.001), supporting an ongoing evolution of CRF85_BC overtime.
A heterosexual transmission of CRF85_BC in Sichuan Province
The epidemiological and clinical features of the patients harboring CRF85_BC did not change significantly between 2014 and 2019 (Table 1). CRF85_BC were only found in heterosexually infected individuals, including commercial sexual transmission (64.7% in 2019 and 42.9% in 2014), non-marital and non-commercial sexual transmission (30.9% in 2019 and 21.4% in 2014) and spousal transmission (4.4% in 2019 and 35.7% in 2014). The majority of CRF85_BC, in 2014(78.6%) and 2019 (95.6%), was circulating among the people living with HIV aged 50 years and over (PLHIV50+), significantly higher than other HIV CRFs (P = 0.000, χ2 test). Patients with low education levels took a relatively larger proportion.
Estimated timeline and migration pathway of CRF85_BC
To elucidate the sources and dissemination routes of endemic CRF85_BC geographically, we constructed maximum clade credibility (MCC) tree based on 149 Sichuan sequences and 15 sequences from other provinces in China using Bayesian geographical evolution method. The corresponding median coefficient of rate variation was 0.33(95% highest posterior density (HPD): 0.18–0.49), supporting the selection of a relaxed molecular clock model. Both Bayesian skyline plot and Bayesian skygrid analysis showed that HIV-1 CRF85_BC experienced 2 fast growth phases during 2011–2014 and 2017–2019 (Fig. 2a) with a mean growth rate of 2.5*10− 3 per year [95% HPD: 2.0*10− 3 - 3.0*10− 3]. The tMRCA was inferred to be 2005.4 (95% HPD: 2002.3–2007.9), implying the origin time of the CRF85_BC (Table 2).
The MCC tree showed that CRF85_BC from Yunnan were located at the root of the tree (posterior probability (PP = 1). Sequentially, virus spread form Yunnan to Sichuan (PP = 0.9911), Guangxi (PP = 0.9978), and other regions. The first place of the viral introduction into Sichuan was probably Yibin in 2008.6 (2006.7–2010.2). An explosive outbreak was then found in the heterosexual elderly people in this city. Several small clusters further suggested that the virus became epidemic in many cities of Sichuan (Luzhou, Guangyuan, et al.) and had rapidly spread to other provinces (such as Liaoning, Anhui and Chongqing) (Fig. 2b).
Notably, most Sichuan sequences originated from Yibin, except for 4 sequences. A Guangyuan and Mianyang sequence directly formed a sub-cluster with a sequence from Yunnan (PP = 0.9965) while a Luzhou strain grouped significantly (PP = 0.9977) with the sequences from Guangxi (Fig. 2b).
Consistent with the differences in viral diversity among geographical regions in China , over 66.5% circulating HIV strains in Sichuan were defined as CRF_BC including CRF07_BC, CRF08_BC and CRF85_BC et al. Hence, CRF_BC infections became a more serious public health issue, which necessitated a deeper understanding of the CRF_BC epidemic patterns in China. CRF85_BC, a newly identified CRF_BC, were the fourth most prevalent HIV strains (68/1291, 5.27%) in Sichuan and was dominating the HIV-1 epidemic in a southwestern city, i e. Yibin city.
In the present study, we only found CRF85_BC prevalent in heterosexually infected individuals. In molecular epidemiological surveillance of Sichuan in 2019, 64.7% of CRF85_BC were transmission via commercial sex. Other early studies about CRF07_BC and CRF08_BC suggested that CRF_BC had spread from IDU to the general population by sexual transmission . Due to lack of the early CRF85_BC sequence, we could not rule out the possibility that this kind of CRF_BC originated from other high-risk groups such as IDUs or former plasma donors (FPDs) . However, commercial sex transmission was at least an important driving factor of CRF85_BC becoming endemic in some local areas. Notably, the majority of cases infected by CRF85_BC were PLHIV50+ cases, most of whom were low-educated. PLHIV50+ have been significantly increasing in the past two decades globally . The proportion of PLHIV50+ around the world increased substantially from 8% in 2000 to 16% in 2016 and is expected to reach an estimated 21% by 2020 . Poor awareness about safer sex practices and underestimating their risk of contracting HIV made them extremely vulnerable to HIV.
Analysis of 148 PR/RT sequences from CRF85_BC infected individuals sampled between 2011 and 2019 revealed that this strain had not just become endemic in Sichuan but also spread to some other provinces in China. Until now, at least 6 provinces have been affected and had a most recent common ancestor (tMRCA) is estimated at around 2005. Yunnan is considered to be the origin area. Previous studies in Yunnan between 2009 and 2012 revealed several CRF_BC HIV-1 strains circulating among local sexually active youth, which is coincident with the gateway position of Yunnan in China’s HIV epidemic history . By Bayesian skyline plot and Bayesian skygrid analysis, the exponential increase in the effective number of CRF85_BC indicated that the virus experienced 2 fast growth phases during 2011–2014 and 2017–2019.
CRF85_BC dispersion showed a strong geographic compartmentalization. After its origin in Yunnan, CRF85_BC almost at the same time spread to neighboring provinces such as Sichuan around 2008.6 (2006.7–2010.2), Guangxi around 2012.3 (2011.3–2013.2). The first place where the virus was introduced into Sichuan was probably Yibin city, which was located in the junction of Yunnan, Guizhou and Sichuan Provinces. In recent years, the local epidemic situation had deteriorated rapidly . The explosive epidemic of CRF85_BC was found in the heterosexual elderly people in this city, and 39.7% (58/146) HIV cases locally were infected by CRF85_BC. Phylogeographical reconstruction further suggested that CRF85_BC spread from Yibin to other cities of Sichuan and even other provinces. This virus was identified in 9 cities in Sichuan, 3 of which probably had an endemic, such Luzhou, Leshan and Guangyuan. Leshan and Luzhou were both bordering with Yibin. Guangyuan is a northeastern city in Sichuan. According to another research, CRF85_BC infected people with Yibin household registration travelled to and resided in Guangyuan, further adding fuel to local endemic.
CRF85_BC entered Sichuan probably through multiple introduction events. In phylogenetic trees, some single strains of Sichuan also formed sub-clusters with Yunnan or Guangxi sequences. However, not all imported cases of Sichuan caused local epidemics.
Our study is the first report to clarify the phylogeographic and demographic dynamics of CRF 85_BC, a newly-identified and increasingly prevalent HIV-1 strain, in Sichuan, China. This study also highlights that heterosexual behavior among older people played an important role in the contribution of CRF 85_BC transmission. The origin of CRF 85_BC was traced back to Yunnan around 2005, spread exponentially to various places (including Guangxi, Sichuan, et al), and later on became endemic around 2008.5 (2007.0–2010.5) in Sichuan. Our findings provide novel insights into the dissemination of a novel HIV-1 virus and extend the understanding of the HIV epidemic in China.
Availability of data and materials
The CRF85_BC sequences determined in this study are available in GenBank under Accession Numbers MT522696 - MT522844. Other datasets used and/or analyzed during the current study available from the corresponding author on reasonable request.
- PLHIV50 +:
The people living with HIV aged 50 years and over
Circulating recombinant forms
Injecting drug users
Los Alamos National Laboratory
Bayesian Markov chain Monte Carlo
Surveillance drug resistance mutations (SDRM)
The most recent common ancestor
Highest posterior density
Lyu P, Chen FF. National HIV/AIDS epidemic estimation and interpretation in China. Zhonghua liu xing bing xue za zhi = Zhonghua liuxingbingxue zazhi. 2019;40(10):1191–6.
Feng Y, Takebe Y, Wei H, He X, Hsi JH, Li Z, Xing H, Ruan Y, Yang Y, Li F, et al. Geographic origin and evolutionary history of China's two predominant HIV-1 circulating recombinant forms, CRF07_BC and CRF08_BC. Sci Rep. 2016;6:19279.
Yang R, Xia X, Kusagawa S, Zhang C, Ben K, Takebe Y. On-going generation of multiple forms of HIV-1 intersubtype recombinants in the Yunnan Province of China. AIDS. 2002;16(10):1401–7.
Nikolopoulos GK, Kostaki EG, Paraskevis D. Overview of HIV molecular epidemiology among people who inject drugs in Europe and Asia. Infect Genet Evol. 2016;46:256–68.
Xiao Y, Kristensen S, Sun J, Lu L, Vermund SH. Expansion of HIV/AIDS in China: lessons from Yunnan Province. Soc Sci Med. 2007;64(3):665–75.
Zhang L, Chow EP, Jing J, Zhuang X, Li X, He M, Sun H, Gorgens M, Wilson D, Wang L, et al. HIV prevalence in China: integration of surveillance data and a systematic review. Lancet Infect Dis. 2013;13(11):955–63.
Xiao P, Li J, Fu G, Zhou Y, Huan X, Yang H. Geographic distribution and temporal trends of HIV-1 subtypes through heterosexual transmission in China: a systematic review and meta-analysis. Int J Environ Res Public Health. 2017;14(7):830.
He X, Xing H, Ruan Y, Hong K, Cheng C, Hu Y, Xin R, Wei J, Feng Y, Hsi JH, et al. A comprehensive mapping of HIV-1 genotypes in various risk groups and regions across China based on a nationwide molecular epidemiologic survey. PLoS One. 2012;7(10):e47289.
Chen H, Wu X, Chen L, Lu H, Tang Z, Shen Z, Pan SW, Ruan Y, Shao Y. Rapidly spreading human immunodeficiency virus epidemic among older males and associated factors: a large-scale prospective cohort study in rural Southwest China. Sex Transm Dis. 2019;46(4):234–9.
Wei H, Liu Y, Feng Y, Hsi J, Xing H, He X, Liao L, Yutaka T, Li J, Shao Y. Genome sequence of a novel HIV-1 circulating recombinant form (CRF57_BC) identified from Yunnan, China. AIDS Res Hum Retrovir. 2014;30(4):384–8.
Hsi J, Wei H, Xing H, Feng Y, He X, Liao L, Jia M, Wang N, Ning C, Shao Y. Genome sequence of a novel HIV-1 circulating recombinant form (CRF64_BC) identified from Yunnan, China. AIDS Res Hum Retrovir. 2014;30(4):389–93.
Su L, Wei D, Yang H, Zeng Y, Hu Y, Yuan D, Feng L, Ruan Y, Qin G, Liang S. Identification of a novel HIV-1 circulating recombinant form (CRF85_BC) in Sichuan, China. AIDS Res Hum Retrovir. 2016;32(9):895–9.
Wu J, Zhang Y, Shen Y, Wang X, Xing H, Yang X, Ding X, Hu B, Li H, Han J, et al. Phylogenetic analysis highlights the role of older people in the transmission of HIV-1 in Fuyang, Anhui Province, China. BMC Infect Dis. 2019;19(1):562.
Chen M, Ma Y, Chen H, Dai J, Luo H, Yang C, Dong L, Jin X, Yang M, Yang L, et al. Spatial clusters of HIV-1 genotypes in a recently infected population in Yunnan, China. BMC Infect Dis. 2019;19(1):669.
Su L, Liang S, Hou X, Zhong P, Wei D, Fu Y, Ye L, Xiong L, Zeng Y, Hu Y, et al. Impact of worker emigration on HIV epidemics in labour export areas: a molecular epidemiology investigation in Guangyuan, China. Sci Rep. 2018;8(1):16046.
NCSTD N. National AIDS venereal disease in December 2017; 2018. p. 02.
Yang Y, Latkin CA, Luan R, Yang C. A cross-sectional study of the feasibility of pharmacy-delivered harm reduction services among people who inject drugs in Xichang, China. BMC Public Health. 2015;15:885.
Yang Y, Latkin C, Luan R, Yang C. Reality and feasibility for pharmacy-delivered services for people who inject drugs in Xichang, China: comparisons between pharmacy staff and people who inject drugs. Int J Drug Policy. 2016;27:113–20.
Li L, Wei D, Hsu WL, Li T, Gui T, Wood C, Liu Y, Li H, Bao Z, Liu S, et al. CRF07_BC strain dominates the HIV-1 epidemic in injection drug users in Liangshan prefecture of Sichuan, China. AIDS Res Hum Retrovir. 2015;31(5):479–87.
Meng Z, Xin R, Abubakar YF, Sun J, Wu H, Lu J, Ni Y, Zhang X, Xu J. Five new CRF07_BC near full-length sequences isolated from Sichuan, China. AIDS Res Hum Retrovir. 2013;29(1):191–7.
Kumar S, Stecher G, Tamura K. MEGA7: molecular evolutionary genetics analysis version 7.0 for bigger datasets. Mol Biol Evol. 2016;33(7):1870–4.
Price MN, Dehal PS, Arkin AP. FastTree 2 – approximately maximum-likelihood trees for large alignments. PLoS One. 2010;5(3):e9490.
Liu K, Linder CR, Warnow T. RAxML and FastTree: comparing two methods for large-scale maximum likelihood phylogeny estimation. PLoS One. 2011;6(11):e27731.
Ling S, Hong Y, Ying H, Li Y, Dan Y, Wei D, Zhou G, Jia Y, Liu H, Qin G. Molecular epidemiological study on HIV-1 in Sichuan province in 2014. Chin J Aids Std. 2016;22(12):946–9.
Ling S, Li Y, Yang L, Jiafeng Z, Shundong W, Shuibin L, Xueqin H, Dan Y, Hong Y, Dongbing W. HIV-1 molecular epidemiological survey of people over 50 years old in Sichuan. Chinese J Aids Std. 2018;24(8):788–92.
Bennett DE, Camacho RJ, Otelea D, Kuritzkes DR, Fleury H, Kiuchi M, Heneine W, Kantor R, Jordan MR, Schapiro JM, et al. Drug resistance mutations for surveillance of transmitted HIV-1 drug-resistance: 2009 update. PLoS One. 2009;4(3):e4724.
Paraskevis D, Beloukas A, Stasinos K, Pantazis N, de Mendoza C, Bannert N, Meyer L, Zangerle R, Gill J, Prins M, et al. HIV-1 molecular transmission clusters in nine European countries and Canada: association with demographic and clinical factors. BMC Med. 2019;17(1):4.
Drummond AJ, Rambaut A. BEAST: Bayesian evolutionary analysis by sampling trees. BMC Evol Biol. 2007;7:214.
Rambaut A SM, Xie D, Drummond AJ.: Tracer v1.6. . http://beast.bio.ed.ac.uk/Tracer 2014.
Lu L, Jia M, Ma Y, Yang L, Chen Z, Ho DD, Jiang Y, Zhang L. The changing face of HIV in China. Nature. 2008;455(7213):609–11.
Mahy M, Autenrieth CS, Stanecki K, Wynd S. Increasing trends in HIV prevalence among people aged 50 years and older: evidence from estimates and survey data. AIDS. 2014;28(Suppl 4):S453–9.
Autenrieth CS, Beck EJ, Stelzle D, Mallouris C, Mahy M, Ghys P. Global and regional trends of people living with HIV aged 50 and over: estimates and projections for 2000-2020. PLoS One. 2018;13(11):e0207005.
Wei H, Xing H, Hsi JH, Jia M, Feng Y, Duan S, He C, Yao S, Ruan Y, He X, et al. The sexually driven epidemic in youths in China's southwestern border region was caused by dynamic emerging multiple recombinant HIV-1 strains. Sci Rep. 2015;5:11323.
Chunlei L. EpidemiologicalCharacteristics ofHIV/AIDS in Yibin City, 2004-2013. J Occup Health Damage. 2015;5:275–77.
The authors would like to thank staffs in Division of Virology and Immunology, National Center for AIDS/STD Control and Prevention (NCAIDS), China CDC for technical assistance. The opinions expressed herein show the collective views of the coauthors and do not necessarily represent the official position of the Sichuan Provincial Center for Disease Control and Prevention, Center for AIDS/STD Control and Prevention.
This study was supported by the Sichuan Science and Technology Project (2017JY0302). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Ethics approval and consent to participate
This study was reviewed and approved by the ethics committees of Sichuan Center for Disease Prevention and Control. All participants signed written informed consents prior to sample donations.
Consent for publication
The authors declare no competing financial interests exist.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Su, L., Feng, Y., Liang, S. et al. The origin and spread of CRF85_BC, driven by heterosexual transmission among older people in Sichuan, China. BMC Infect Dis 20, 772 (2020). https://0-doi-org.brum.beds.ac.uk/10.1186/s12879-020-05488-4
- Heterosexually transmission