Prevalence of resistance to integrase strand-transfer inhibitors (INSTIs) among untreated HIV-1 infected patients in Morocco

Objective The integrase strand-transfer inhibitors (INSTIs) are an important class in the arsenal of antiretroviral drugs designed to block the integration of HIV-1 cDNA into the host DNA through the inhibition of DNA strand transfer. In this study for the first time in Morocco, the complete HIV-1 integrase gene was analysed from newly diagnosed patients to evaluate the prevalence of natural polymorphisms and INSTIs resistance-associated mutations in the integrase gene. Results The 864pb IN coding region was successfully sequenced from plasma sample for 77 among 80 antiretroviral naïve patients. The sequences were interpreted for drug resistance according to the Stanford algorithm. Sixty samples were HIV-1 subtype B (78%), fourteen CRF02_AG (18%), two subtype C and one subtype A. Overall 81 of 288 (28%) amino acid IN positions presented at least one polymorphism each. We found 18 (36.73%), 42 (25.76%) and 21 (27.27%) of polymorphic residues assigned to the N-Terminal Domain, Catalytic Core Domaine and the C-Terminal Domain positions respectively. Primary INSTIs resistance mutation were absent, however secondary mutations L74IM, T97A were detected in four samples (5.2%). These results demonstrate that untreated HIV-1 infected Moroccans will be susceptible to INSTIs. Electronic supplementary material The online version of this article (10.1186/s13104-018-3492-5) contains supplementary material, which is available to authorized users.


Introduction
Since the introduction of combination therapy (highly active antiretroviral therapy, HAART) with protease inhibitors (PIs), nucleoside reverse transcriptase inhibitors (NRTIs) and non-nucleoside reverse transcriptase inhibitors (NNRTIs) in Morocco in 1998, the mortality and morbidity of HIV/AIDS patients has reduced [1]. These drugs suppress viral replication and reduce HIV viral RNA loads in the plasma of patients, thus helping to maintain the immune system, but they do not prevent Therefore, they have become an essential component of HAART used in many countries. Morocco has installed national programs following World Health Organization ARV guidelines and newer ARV salvage regimens including third-line drugs such as INSTIs will been introduced in the upcoming years. The aim of this study is to analyze, for the first time in Morocco, integrase (IN) sequence variability among ARV treatment naïve patients to determine the frequency of resistance mutations and the prevalence of natural polymorphisms of the IN gene and in order to estimate INSTIs efficacy prior to their introduction into the country.

Study population and samples
Plasma samples were collected for genotypic assay of the IN gene region from eighty seven HIV-infected, antiretroviral therapy-naïve patients originating from different geographic parts of the country enrolling at the dermatology department of Mohammed V Military Teaching Hospital in Rabat between the years 2009 and 2015. Demographic, clinical and laboratory data were collected for all patients. The quantitative HIV-RNA tests were performed using Cobas TaqMan HIV-1 Test, version 1.0 (Roche Diagnostics Systems, Germany, P/N: 03542998 190). CD4 cell enumerations were performed using the FacsCount instrumentation. (FacsCount, Becton-Dickinson, P/N: 339010). After viral load testing,eighty samples were used for sequencing assay. The detection of IN mutations by sequencing was unsuccessful for three samples.
The study was approved by the Ethical Committee of Biological Research, Faculty of Medicine and Pharmacy-Rabat, and was conducted with respect to legal aspects. Written informed consent was obtained from all participants before any data analysis procedure.

Genetic analysis and drugs resistance
HIV RNA was extracted from plasma using High Pure Viral RNA Kit (Roche Diagnostics Systems, Germany, P/N: 11858882001) and the integrase coding region (867 bp) was amplified by one-step reverse transcriptase polymerase chain reaction (RT-PCR) using MyTaq One- Step RT-PCR kit (Bioline, London, UK, P/N:  and the primer set KVL068 -KVL069 [16]. The Nested-PCR assay was carried out using MyFi DNA Polymerase kit (Bioline, London, UK, P/N: BIO-21118) and the primer set KVL070 and KVL084 [16]. Sequencing reaction was performed using BigDye Terminator v3.1 Ready Reaction Cycle Sequencing Kit (P/N: 4337455) with an ABI PRISM 3130XL Genetic Analyzer (Applied Biosystems) using the POP-7 polymer (P/N: 4393708). Data were analyzed by sequencing Analysis Software version 5.3.1 (Applied Biosystems, P/N: 4360967).
IN Sequences were assembled and aligned using DNA Dragon Sequence Assembler version 1.6.0 (Sequentix-Digital DNA Processing, Germany) and Muscle method in MEGA 6 software [17], respectively. All sequences were submitted to GenBank and registered under accession numbers: KU609274-KU609350.
HIV-1 subtyping and screening of IN polymorphisms in comparison with the HxB2 HIV-1 clade B consensus sequence (GenBank accession number K03455.1) were done using geno2pheno subtyping tools [18]. Phylogenetic tree was constructed by using Maximum Likelihood method, and Bootstrap resampling was performed 1.000 times for all sequences with MEGA 6 software ( Fig. 1). INSTIs resistance mutations in IN sequences were interpreted using the Stanford HIVdb Program (Version September 23,2016) and screened for the presence of additional changes in 17 positions (V72I, T112I, S119PRTG, T124A, T125K, A128T, Q146K, M154I, K156N, V165I, V201I, I203M,T206S, S230N, D232N, V249I and C280Y) previously related to INIs resistance in vitro and frequently reported by different studies [19][20][21][22][23]. The variability at the D 64 D 116 E 152 and H 12 H 16 C 40 C 43 motifs, the residues that interact with the human lens epithelium derived growth factor (LEDGF/P75) and the degree of variability of the three functional areas of IN were also investigated.

Results
Out of eighty patients, sixty-two (77.5%) were male. Median age was 36 years old, where patients with an age between 25 and 44 years old represent 61%. Sixty-five (81.25%) patients were suspected to have acquired HIV infection through heterosexual contact, four (5%) perinatally, and the mode of infection was unknown for 11 patients. The median of CD4+ T cell count and viral load at the time of sequencing for the available values were 409 cells/mm 3 and 95,800 copies/ml respectively (Table 1).
In 3 samples the sequencing failed and subtypes were unavailable. The final dataset for the baseline IN resistance included 77 individuals. IN sequencing was consistently successful at HIV viral loads higher than 66 copies/ ml. Sixty (77.92%) were of HIV-1 subtype B, fourteen (18.18%) CRF02_AG, two (2.6%) subtype C and one (1.3%) subtype A. These data are also summarised in the IN phylogenetic tree shown in listed in the Stanford resistance algorithm were found in this study, while secondary mutations L74IM and T97A associated with drug resistance to RAL and EVG were observed only in four patients (5.2%) (Additional file 1). The L74I was observed in two strains, a subtype B (KU609297) and a subtype A (KU609282). The L74M was observed in one subtype B strain (KU609337) and T97A was observed in one CRF02_AG strain (KU609303). No mutation associated with DTG resistance was observed among all studied patients. Other mutations in IN positions not included in Stanford list but frequently reported  Not determined 16 20 in regards to INIs resistance in vitro V72I, T112I, S119PRT, T124A, K156N, V165I, V201I, I203M,T206S and S230N have been detected with different frequencies (Fig. 2). Finally there were no significant associations (p > 0.05) between resistance mutations with exposure category, viral load, Transmission route and CD4 T cell count.

Discussion
Three INSTIs, RAL, EVG and DTG have been approved for clinical use by the FDA and European Medicines Agency [6][7][8]. These compounds have proven to be highly efficient for the treatment of both ARV-naïve and -experienced individuals even with preexisting drug resistance or treatment complication [9][10][11][12][13][14][15], thus INIs has rapidly became an important class in the arsenal of ARV drugs. Eighty HIV-1 untreated patients were recruited in this study in order to examine the variability within the IN gene at positions associated with resistance to INIs. Despite the limited size of the study population, these results are in good agreement with the current situation of HIV in Morocco as reported in the national report on AIDS [1]. The high prevalence of subtype B among infected Moroccan population has been previously reported by other studies which suggested that it was due to the presence of a close relationship between Morocco and European countries [3,4]. Detailed phylogenetic analysis of IN sequences showed that all CRF02_ AG isolates were related to strains found in central Africa and Europe, which agreed with previous reports that suggested that increasing prevalence of CRF02_AG might be associated with increasing immigration from sub-Saharan Africa to Europe via Morocco [27,28].  E75 D59  R10  N 2  K 46 T6  V1 N4 E43  I 6  I18 I6  Q1  I2  T1  C9  N1  Q45  I 7  T1  S1  G44  M 1  H 1  I 1  L 4  V1  T2   51  60  70  80  90 100 I2  I1  M4  G1  T2  A1 G2  N1  Q43 T1  M 1   101  110  120  130  1 40 150  or ARV-experienced [29]. This digit indicates the relative conservation of the protein in untreated Moroccan's patients. The analysis showed conservation of the HHCC Zn + -binding motif, the catalytic triad DDE, and several important IN residues involved in the chemical bond and hydrophobic contact with LEDGF/P75; which an essential HIV integration cofactor linking IN to chromatin [24][25][26]. The conservation of these specific structural domains is strictly necessary for the correct performance of HIV-1 IN functions [30]. Importantly the major resistance mutations with reduced susceptibility to RAL, EVG and DTG were totally absent. The absence of such mutations in our study is consistent with the results of several studies in treatmentnaïve patients [19,20,[31][32][33][34][35] and with the fact that the transmission of the drug resistance is unlikely in populations previously unexposed to INIs treatment [36]. Only three secondary drug resistance mutations included in Stanford list were observed in 4 strains. These mutations have been previously been described as polymorphic, occurring in 1-2% of IN sequences, observed in subtypes A, B, C, D, CRF01_AE and CRF02_AG [29]. These mutations contribute to INI resistance only in the presence of primary INI resistance mutations [19,37]. Moreover, there were three amino acid substitutions of unknown significance at position 163 that were encountered in subtype B and CRF02-A/G strains: G163E, G163T and G163Q. In IN residues, it is usually considered non polymorphic in all subtypes except subtype F [38]. Regarding DTG R263K resistance mutation, no strain from our study exhibited this mutation, whereas the L101I and T124A mutations were found in 26 and 12 strains, respectively. These mutations were previously shown to be selected in vitro in the presence of DTG and have shown little impact on virological response to DTG [39,40]. Conversely, specific additional mutations in amino acid IN positions 72, 112, 119,156, 165, 201, 203, 206, and 230 occurred with different prevalence in subtype B and non-B HIV-1 variant were more common. These mutations have not been described to be associated with RAL or EGV resistance [20,21]. In the same way, it has been reported that in the absence of primary mutation, all these secondary mutations had little if any effect on drug susceptibility in vitro, thus suggesting rather a secondary role for viral fitness rescue and/or increasing resistance [19]. Furthermore many previous genotypic studies on HIV-1 IN in treatment-naïve patients living with various viral subtypes in different countries: England, Spain, South Africa, Sub-Saharan Countries, Thailand, Indonesia and Korea have showed that these differences are natural polymorphisms [32][33][34][41][42][43]. According to previous studies, and the fact that INIs have not yet been introduced in Morocco, all secondary and additional mutations identified in this study are also likely natural polymorphisms.
In conclusion, these results demonstrate that untreated HIV-1 infected Moroccans are likely to benefit from INSTI-based drug regimens, particularly given the rising issues related to drug resistance against reverse transcriptase inhibitors that are currently used in Morocco.