{'fraction_non_alphanumeric': 0.037528307990941444, 'fraction_numerical': 0.03105791006146878, 'mean_word_length': 4.919100047869795, 'pattern_counts': {'":': 0, '<': 1, '<?xml version=': 0, '>': 8, 'https://': 0, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 1, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
The prevalence of drug resistance among treatment-naïve HIV-1-infected individuals in China during pre- and post- 2004

Background The widespread use of antiretroviral therapies has led to considerable concerns about the prevalence of drug-resistant, as transmission of drug-resistant (TDR) strains poses a challenge for the control of the HIV-1 epidemic. Methods We conducted an epidemiological study enrolling treatment-naïve HIV-1-positive subjects at the Peking Union Medical College Hospital since 1991. Drug resistance was determined by submitting the sequences to the Stanford University Network HIV-1 database. Results Of 521 participants, 478 samples were amplified and sequenced successfully. HIV Transmitted drug resistance prevalence in China was determined to be 6.7 %. We did not find significant differences in the TDR rate by demographic characteristics. No significant time trend in the prevalence of overall TDR was observed (p > 0.05). Conclusions We identified an intermediate prevalence of transmitted drug resistance (TDR), exhibiting a stable time trend. These findings enhance our understanding of HIV-1 drug resistance prevalence and time trend, and provide some guidelines for the comprehensive public health strategy of TDR prevention.


Background
The "Four Frees and One Care" policy implemented in China in 2003, which provided free antiretroviral (ARV) drugs, has been credited for significantly reducing the rate of morbidity and mortality attributed to HIV infection in the country [1,2]. By August 31st 2013, HIV-1 patients living in 31 provinces and 2215 counties had received antiretroviral therapy (ART) and greatly benefited from this program. The accumulative total of treatment-experienced patients was 234,655, with 188,126 being treated in China.
Characterized by a high replication rate [3] and the error-prone reverse transcriptase [4,5], HIV-1 is extensively genetically diverse. Like those in developed countries, the emergence and transmission of drugresistant HIV-1 strains is of considerable concern with the widespread use of ART in developing countries. It has been reported that many treated patients developed resistance to a specific drug within a class or an entire class of ART drugs [6], which may lead to treatment failure and limitation of alteration of treatment regimen. Furthermore, transmitted drug resistance (TDR), which is the transmission of drug-resistant variants to individuals who have never undergone ART can occur. TDR can compromise the efficacy of combination ARV regimens in the management of HIV disease, becoming a major obstacle to ART. The prevalence of TDR varies worldwide ranging from 3.2 % to 24.1 % [2,[7][8][9][10][11][12][13][14][15][16][17][18][19][20][21][22]. In China, free therapy was provided to AIDS patients in 2004, which means that there might be transmission of drug-resistant variants. Some reports have demonstrated the increasing TDR rate [22,23]. So in this study we focused on the compare of the prevalent ratio of drug-resistant variants before and after 2004 in China to provide some guidelines for the TDR prevention.

Patient samples
Plasma specimens were collected from HIV-1 treatmentnaïve patients enrolled in the Peking Union Medical College Hospital. Participants had received no antiretroviral drugs at the time of sampling. The patients recruited in this study were AIDS stage for their CD4 cells being below 350 counts per microlitre and viral load exceeding 10 3 log. Plasma samples were shipped on dry ice and stored at -80°C until analysis. All patients provided a written informed consent and their basic personal and clinical information was acquired by face to face interview. A total of 521 patients were enrolled, and this study was approved by the Ethics Committee at the Institutional Review Board of the Academy of Military Medical Sciences.

Subtype identification
Subtype was determined by phylogenetic tree analysis of the PR-RT fragments. Alignments also included the reference sequences representing HIV-1 genetic circulating recombinant forms obtained from the HIV Database [25]. The phylogenetic tree was constructed using the Neighbor-Joining method by MEGA (Molecular Evolutionary Genetic Analysis Software, Version 5.4) and genetic distance reliability was evaluated using the bootstrap test (500 bootstrap replicates).

Statistical analysis
Fisher test analysis was performed using SAS version 9.2. All analyses were performed using a two-tailed test and p values less than 0.05 were defined as statistically significant.

Patient characteristics
The  Table 1 summarizes demographic and clinical characteristics of the study population. The average age was 35.49 ± 10.01 (18-65 years) and males dominated the study population (69.3 %). More than half (77.9 %) of patients were infected by sexual contact, of which 214 (41.1 %) was through heterosexual contact and 182 (34.9 %) through sex between men. The mean HIV viral RNA load was 4.54 log copies/ml.
Comparison of the resistant rate in two major risk groups revealed that there was no significant difference (p > 0.05) between sexually infected individuals (28/373, 7.5 %) and patients infected by blood transfusions (3/43, 7.0 %). In addition, we did not find a significant difference in resistance in patients infected through heterosexual contact versus homosexual contact (p > 0.05). The prevalence of resistance in male and female was 8.1 % and 2.7 %, respectively (p > 0.05). Shenzhen had the highest rate of drug resistance at 11.6 %, and there was no significant difference (p > 0.05) among these regions (Table 3). There were different prevalence of resistance among subtype B, CRF_01AE and CRF_07BC/08BC (p = 0.0037). The susceptibility of CRF01_AE was significantly higher than subtype B (p = 0.0030).
As a result of the introduction of free ART in 2004, we compared the resistance between the two periods. No significant difference (p > 0.05) on resistance was observed in the 1991-2004 (2/28, 7.1 %) period and 2005-2009 (30/450, 6.7 %) period. And there is also no difference in TDR prevalence between the two periods in each origins (p > 0.05). Additionally, the differences on resistance to NRTIs, NNRTIs and PIs between the two periods were not statistically significant (p > 0.05, Fig. 3).

Discussion
We report here the drug susceptibility of 521 treatmentnaïve HIV-1-infected individuals at the Peking Union Medical College Hospital from 1991 to 2009. Four hundred and seventy eight gag-pol gene sequences (478/ 521) were sequenced successfully and the results of subtype and genetic resistance were obtained.
Various HIV-1 strains have been identified in China, of which CRF07_BC, CRF01_AE, CRF08_BC, and B/ B'are the most common. Other reported strains include subtype G, subtype C, CRF55_01B and numerous unique recombinant forms (URFs) [26,27]. Various HIV-1 subtypes/CRFs have been disseminated among different high risk populations. Our findings found that the most common subtype was CRF01_AE among the patients who were infected by sexual contact, and subtype B was  found to circulate mainly in the blood-borne transmission population, which were consistent with previous Chinese studies [28,29]. In the present study, we found that the main route of transmission remained blood-born route for the subjects in Henan province, of which subtype B was dominant, as a result of HIV infection outbreak in the mid-1990s in paid blood donors in Henan.
Our studies determined that the TDR prevalence in China was 6.7 %, a level classified as "intermediate" according to the WHO thresholds, and which was slightly lower than that observed in some industrialized countries. In the present study, we found no significant differences in the presence of TDR by demographic characteristic, which was similar with a previous study [8]. Here we showed that there was no significant difference in the prevalence of drug-resistant mutations among the three antiretroviral classes (PI, 2.7 %; NRTI, 2.5 %; and NNRTI, 2.3 %). In addition, we noted that the prevalence of HIV-1 strains with high level of resistance mutations to NNRTIs were significantly higher than those to NRTIs and PIs (p < 0.0050). The higher prevalence may be related to mutations associated with decreased susceptibility to NNRTIs with a low genetic barrier, which are rapidly generated and emerge early in the selection process [30].
It is known that there is an approximate 8-10 % genetic diversity exists in pol gene among various subtypes; therefore drug resistance susceptibility of genotype virus may be different. Here we showed a strong association between the prevalence of TDR and subtypes. There was different prevalence of resistance among subtype B, CRF_01AE and CRF_07BC/08BC. Resistance was significantly higher in ARV-naïve patients who were infected with subtype B than in those infected with CRF01_AE viruses (p = 0.0030), which was similar with previously published studies [17-19, 31, 32]. Although subtypes may vary in mutational pathways, it is unclear whether subtype B was more prone to mutation [33]. As the predominant prevalence in China, HIV-1 subtype B viruses have been exposed to ARV drugs for a longer period than other subtypes [26,34].
This study was performed on ART-naïve individuals. It has been reported that resistance testing prior to the initiation of treatment in untreated populations with HIV infection is cost-effective and may be beneficial to the patient [35,36]. It is commonly considered that resistant HIV-1 variants in ART-naïve patients are transmitted from ART-experienced patients or other ART-naïve individuals with drug-resistant strains. There are some studies indicating that drug resistant mutations among treatment-naïve patients were associated with ART used prior to and early in the HAART era [16]. As the ART use has been extended in China in recent years, it becomes increasingly important to monitor HIV-1 genetic diversity and TDR. Hence, we also tested whether the free ART program, which was implemented in 2003, has influenced the prevalence of overall TDR among ARTnaïve patients in China. It should be noted that contrary to the reports from other countries that have implemented broad-access programs to ART [37,38], results of our analysis showed that there was no timedependent trend in the overall TDR, origins, or in any class of drug resistance. Although we did not observe evidence of an increasing prevalence of TDR over years in the sampling patients, the number of HIV-1-infected individuals in our study may be small. While we did not study a random sample of HIV-1infected individuals, this was a large, diverse population with a wide time span from a widespread referral network in China, and our study well reflects the time trend and demographic characteristics of TDR.

Conclusions
In summary, in this study, we fulfilled the comprehensive investigation on the circulation of TDR of three ARV classes in China. The overall prevalence of TDR remains intermediate in ART-naïve individuals, and has a stable TDR time trend. These findings enhance our understanding of HIV-1 drug resistance prevalence and time trend, and provide some guidelines for proposing efficacious and effective programs to prevent transmission of HIV-1drug resistant strains.

Acknowledgments
We thank the patients who participate in the study; and the physicians and nurses for excellent patient care.

Funding
This study was funded by The National Key S&T Special Projects on Major Infectious Disease grant (2012ZX10001-002).

Availability of data and materials
The datasets during and/or analysed during the current study available from the corresponding author on reasonable request.

Authors' contributions
HL and SC conceived this study, analyze data, and drafted the manuscript. YH collected samples. DZ, LL, YL, SL, and ZB carried out the sequence alignment, and edit manuscript. WZ and HS participated in the design of the study, and check all data. TL and JL conceived the study, and participated in its design and coordination. All authors read and approved the final manuscript.

Competing interests
The authors declare that they have no competing interests.

Consent for publication
Not applicable.
Ethics approval and consent to participate Academy of Military Medical Science ethics committee that approved the study.