A genotypic assay for the amplification and sequencing of integrase from diverse HIV-1 group M subtypes

A genotypic assay for the amplification and sequencing of integrase from diverse HIV-1 group M subtypes

Journal of Virological Methods 153 (2008) 176–181 Contents lists available at ScienceDirect Journal of Virological Methods journal homepage: www.els...

156KB Sizes 0 Downloads 18 Views

Journal of Virological Methods 153 (2008) 176–181

Contents lists available at ScienceDirect

Journal of Virological Methods journal homepage: www.elsevier.com/locate/jviromet

A genotypic assay for the amplification and sequencing of integrase from diverse HIV-1 group M subtypes Kristel Van Laethem a,b,∗ , Yoeri Schrooten a,b , Kris Covens a , Nathalie Dekeersmaeker a , Paul De Munter c , Eric Van Wijngaerden c , Marc Van Ranst a,b , Anne-Mieke Vandamme a a b c

Rega Institute for Medical Research, Katholieke Universiteit Leuven, Leuven, Belgium Aids Reference Laboratory, University Hospitals Leuven, Leuven, Belgium Internal Medicine, University Hospitals Leuven, Leuven, Belgium

a b s t r a c t Article history: Received 6 March 2008 Received in revised form 14 July 2008 Accepted 17 July 2008 Available online 2 September 2008 Keywords: PCR Sequencing HIV-1 subtypes Resistance

Recently, the Food and Drug Administration (FDA) of the USA approved the first integrase inhibitor for inclusion in treatment regimens of HIV-1 patients failing their current regimens with multi-drug resistant strains. However, treatment failure has been observed during integrase inhibitor-containing therapy. Several mutational pathways have been described with signature mutations at integrase positions 66, 92, 148 and 155. Therefore, a genotypic assay for the amplification and sequencing of HIV-1 integrase was developed. The assay displayed a detection limit of 10 HIV-1 IIIB RNA copies/ml plasma. As the HIV-1 pandemic is characterised by a large genetic diversity, the new assay was evaluated on a panel of 74 genetically divergent samples belonging to the following genetic forms A, B, C, D, F, G, J, CRF01-AE, CRF02-AG, CRFF03-AB, CRF12-BF and CRF13-cpx. Their viral load ranged from 178 until >500,000 RNA copies/ml. The amplification and sequencing was successful for 70 samples (a success rate of 95%). The four failures were most probably due to low viral load or poor quality of RNA and not to subtype issues. Some of the sequences obtained from integrase inhibitor-naïve patients displayed polymorphisms at integrase positions associated with resistance: 74IV, 138D, 151I, 157Q and 163AE. The relevance of these polymorphisms in the absence of the signature mutations remains unclear. © 2008 Elsevier B.V. All rights reserved.

1. Introduction Antiretroviral drug combinations for the treatment of patients infected with HIV-1 generally comprise three drugs that target mostly the viral enzymes reverse transcriptase and protease (DHHS Guidelines, 2007). Recently, the Food and Drug Administration (FDA) approved the first integrase (IN) inhibitor (INI) for use in combination with other antiretroviral drugs for the treatment of therapy-experienced HIV-1 patients failing their current regimens with multi-drug resistant strains. This decision was based upon the efficacy by which raltegravir in combination with an optimized background regimen could improve the viral load and CD4+ count at week 24 in phase III clinical trials (Cooper et al., 2007; Steigbigel et al., 2007; Grinsztejn et al., 2007). However, treatment failure, whether caused by lack of adherence, insufficient

∗ Corresponding author at: Rega Institute for Medical Research and University Hospitals Leuven, Microbiology and Immunology, Clinical and Epidemiological Virology, AIDS Reference Laboratory, Minderbroedersstraat 10, 3000 Leuven, Belgium. Tel.: +32 16 332177; fax: +32 16 332131. E-mail address: [email protected] (K. Van Laethem). 0166-0934/$ – see front matter © 2008 Elsevier B.V. All rights reserved. doi:10.1016/j.jviromet.2008.07.008

potency of the regimen or pre-existing antiviral drug resistance, has been observed during raltegravir-containing therapy. Several mutational pathways have been described, mainly within subtype B strains, with signature mutations at integrase positions 148 and 155 (Hazuda et al., 2007). The majority of HIV-1 infections worldwide, and also in Belgium, are caused by non-B subtypes (Snoeck et al., 2004; Vercauteren et al., 2008). Therefore, an IN genotypic resistance assay performing well for all HIV-1 subtypes was developed. 2. Materials and methods 2.1. Samples Dilution series of HIV-1 (IIIB ) virus (subtype B) were used (kindly provided by R.C. Gallo when at the National Institutes of Health, Bethesda, MD) to fine-tune conditions for RNA extraction, cDNA synthesis and PCR. In addition to this laboratory strain, 74 plasma samples were included from patients who had attended the University Hospitals in Leuven and for whom a genotypic resistance analysis towards protease and reverse transcriptase inhibitors had been performed as part of their routine management. The sam-

K. Van Laethem et al. / Journal of Virological Methods 153 (2008) 176–181 Table 1 Characteristics of the patient samples tested with the IN assay

Table 1 (Continued ) Genetic forma

Genetic forma

Sample

A

AR07-2743 AR02-689 AR04-1158 AR04-1145 AR02-852 AR04-320 AR05-073 P00-235 AR02-025

479 1,521 6,908 13,204 17,390 18,226 56,855 77,055 94,885

AR07-1935 AR07-2127 AR02-781 AR06-1245 AR07-1965 AR00-175 AR03-127 AR01-751 AR01-273 P00-470 P00-200 AR02-092 H97-3457 AR02-884 AR02-471

398 447 526 575 1,479 3,009 4,172 4,273 12,492 18,644 41,032 44,411 61,722 254,839 >500,000

1(+) 1(+) 1(−)/2(+) 1(+) 1(+) 1(+) 1(+) 1(+) 1(+) 1(+) 1(+) 1(+) 1(+) 1(+) 1(+)

AR07-955 AR05-124 AR07-050 AR07-1554 AR05-390 AR05-1284 AR04-354 AR05-756 AR06-505 AR01-760 AR01-491

562 569 1,096 1,380 2,649 4,778 7,324 11,378 23,988 50,118 138,739

1(+) 1(−)/2(−) 1(−)/2(+) 1(+) 1(+) 1(+) 1(+) 1(+) 1(+) 1(+) 1(+)

CRF01

AR02-766 AR02-969 AR04-1031 AR03-1252 AR03-1146 AR05-781 AR05-673 P98-015

2,017 3,629 9,765 10,554 11,400 59,544 141,551 145,465

1(+) 1(+) 1(+) 1(+) 1(−) 1(+) 1(+) 1(+)

CRF02

AR06-602 AR07-1637 AR06-309 AR05-569 AR05-1056 AR03-228 AR05-747 AR03-1196 AR05-1089

269 2,042 2,750 4,883 5,446 7,124 14,523 45,293 318,615

1(+) 1(−)/2(+) 1(−) 1(+) 1(+) 1(−)/2(+) 1(+) 1(+) 1(+)

B

C

Viral loadb

177

Successc 1(+) 1(+) 1(+) 1(+) 1(+) 1(+) 1(+) 1(+) 1(+)

CRF03

AR03-506

2,000

1(−)

CRF12

P99-538

CRF13

AR03-754 P00-530 AR06-225

80,577 767 29,979 102,329

1(+) 1(−)/2(+) 1(+) 1(+)

D

AR02-071 AR01-109 AR04-342 AR01-699 AR06-2632 AR04-1171

7,943 12,589 14,921 23,000 25,704 100,670

1(+) 1(+) 1(+) 1(+) 1(+) 1(+)

F

AR06-659 AR04-538 AR05-661

178 9,422 13,714

G

AR06-1056 AR01-066 AR06-1186

1,318 2,268 2,818

1(−)/2(+) 1(+) 1(+) 1(+) 1(+) 1(+)

J

Sample

Viral loadb

Successc

AR06-009 AR06-2205

25,349 27,542

1(+) 1(+)

AR00-287 AR03-1144 AR01-485

2,219 7,939 419,945

1(+) 1(+) 1(+)

a Genetic forms as obtained by submitting the PR-RT-IN sequences to the Rega HIV-1 subtyping tool (de Oliveira et al., 2005), except for samples that failed, for which only the PR-RT sequence was used. b RNA copies/ml. c 1(+) and 1(−)/2(+): successful amplification and sequencing in the first, respectively, second round. 1(−)/2(−): unsuccessful amplification and sequencing in both rounds. 1(−): unsuccessful amplification and sequencing in the first round and for which no samples were available to perform a second round.

ples were selected based on the subtype that was determined by submitting the pol sequences to the Rega HIV-1 subtyping tool (de Oliveira et al., 2005) (Table 1). The plasma viral load was determined using the VERSANT HIV-1 RNA 3.0 Assay (bDNA, Bayer HealthCare, Brussel, Belgium) and Abbott RealTime HIV-1 (Abbott Molecular, Louvain-La-Neuve, Belgium). 2.2. Primer development and synthesis The PCR and sequencing primers were designed using the polvif gene alignment of several HIV-1 group M strains (Leitner et al., 2005). The primer sequences are listed in Table 2. Primers were developed and analysed using the Oligo software (Medprobe, Oslo, Norway). The primers were synthesised by Invitrogen (Merelbeke, Belgium). 2.3. Extraction For HIV-1 viral RNA detection, 1 ml plasma samples or 1 ml virus IIIB supernatant were ultracentrifuged at 37,100 × g for 1 h to pellet the virus and were extracted subsequently using the extraction procedure of the Viroseq HIV-1 Genotyping System (Abbott Molecular, Louvain-La-Neuve, Belgium). 2.4. cDNA synthesis and amplification of the IN region Ten microliters out of 50 ␮l RNA extract was reverse transcribed and amplified in a one-step RT-PCR using the SuperScriptTM III One-Step RT-PCR System with Platinum® Taq High Fidelity (Invitrogen, Merelbeke, Belgium) and the following conditions: 1× Reaction Mix, 2.25 mM MgSO4 , 0.2 ␮M primer KVL068, 0.2 ␮M primer KVL069, 10 U Protector RNase Inhibitor (Roche Diagnostics, Vilvoorde, Belgium) and 1 ␮l SuperScriptTM III RT/Platinum® Taq High Fidelity Enzyme Mix. Cycling conditions were: a cDNAsynthesis step of 30 min at 55 ◦ C, a denaturation step of 2 min at 94 ◦ C, 40 cycles of 15 s at 94 ◦ C, 30 s at 53 ◦ C, 2 min 30 s at 68 ◦ C and a final extension step of 5 min at 68 ◦ C. A 2128-bp nucleotide fragment, encompassing the integrase region, was herewith obtained. A 1254-bp nucleotide fragment of integrase was amplified by adding 1 ␮l outer PCR product into in an inner PCR using the Expand High Fidelity PCR System and the following conditions: 1× Expand HF Buffer, 2 mM MgCl2 , 200 ␮M dNTPs, 0.4 ␮M primer KVL070, 0.4 ␮M primer KVL084 and 2.625 U Expand High Fidelity PCR System enzyme mix. Cycling conditions were: a denaturation step of 2 min at 94 ◦ C, 10 cycles of 15 s at 94 ◦ C, 30 s at 53 ◦ C, 1 min 30 s at 72 ◦ C, 30 cycles of 15 s at 94 ◦ C, 30 s at 53 ◦ C, 1 min 30 s at 72 ◦ C + cycle elongation of 5 s for each cycle and a final extension step of 7 min at 72 ◦ C. Amplification products were separated on a 1% agarose gel and visualized by ethidium bromide staining.

178

K. Van Laethem et al. / Journal of Virological Methods 153 (2008) 176–181

Table 2 Primers for the amplification and sequencing of the HIV-1 integrase region Primers

Sequence 5 –3

Positiona

Description

KVL068 KVL069 KVL070 KVL084 KVL076 KVL082 KVL083

AGGAGCAGAAACTTWCTATGTAGATGG TTCTTCCTGCCATAGGARATGCCTAAG TTCRGGATYAGAAGTAAAYATAGTAACAG TCCTGTATGCARACCCCAATATG GCACAYAAAGGRATTGGAGGAAATGAAC GGVATTCCCTACAATCCCCAAAG GAATACTGCCATTTGTACTGCTG

3854–3880 5955–5981 4013–4042 5243–5266 4161–4188 4647–4669 4750–4772

Sense outer primer for amplification Antisense outer primer for amplification Sense inner primer for amplification Antisense inner primer for amplification and sequencing Sense primer for sequencing Sense primer for sequencing Antisense primer for sequencing

a

Positions according to pNL4.3 (AF324493).

The images were processed on a videoimager (ImageMaster VDS, Amersham Biosciences, Roosendaal, The Netherlands). 2.5. Sequencing and data analysis PCR products for population sequencing were purified with Microspin S-400 (Amersham Biosciences, Roosendaal, The Netherlands). Sequencing was performed using the ABI PRISM BigDye Terminator v3.1 Ready Reaction Cycle Sequencing Kit and the primers KVL084, KVL076, KVL082 and KVL083 (Table 2). The sequencing reactions were run on an ABI3100 Genetic Analyzer (Applera, Nieuwerkerk a/d Ijssel, The Netherlands). The sequences were analysed using Sequence Analysis version 3.7 and SeqScape version 2.0 (Applera, Nieuwerkerk a/d Ijssel, The Netherlands). The genetic form of the newly obtained IN fragments was obtained by submitting the sequence to the Rega HIV-1 subtyping tool (de Oliveira et al., 2005). The variability within the IN was defined as the number of observed amino acids divided by the frequency of the most prevalent amino acid at a certain position. The variability can range theoretically between 1 (1/1; completely conserved) and 400 (20/0.05). The accession numbers of the sequences obtained in this manuscript are EU882742–EU882811. 3. Results 3.1. Performance of the assay Experimental conditions were optimized starting from dilution series in PBS of HIV-1 IIIB, ranging from 1,000,000 until 10 RNA copies/ml (final conditions as described in Section 2). The optimized nested PCR procedure was specific, generating only a single amplification band, and it was very sensitive (10 RNA copies/ml). Subsequently, it was validated on a set of 74 clinical samples from INI-naïve patients, belonging to a broad range of HIV-1 group M subtypes (A, B, C, D, F, G, J, CRF01, CRF02, CRF03, CRF12, CRF13) and displaying a viral load range between 178 and >500,000 RNA copies/ml (Table 1). At first, the amplification was successful for 64 samples (64/74; 86% success rate). The negative samples were of subtype B, C, F, CRF01, CRF02, CRF03 and CRF13, and half had a viral load below 1000 RNA copies/ml. As the IN amplification was performed on RNA extracts, obtained previously for routine resistance testing purposes and stored for many years at −80 ◦ C, 7 of the 10 negative samples were re-extracted. Plasma was not available for the three remaining negative samples. In this run, only one sample remained unamplifiable (AR05-124, 569 RNA copies/ml), giving a total success rate of 95% (70/74). Full-length sequences were obtained for all 70 amplified samples (70/74; 95%). Ninety-four percent of the IN sequences were covered fully in both directions (66 samples sequenced in sense and antisense direction). The four other IN sequences were partly unidirectional. When amplification was successful, primers KVL083 and KVL084 were always able to generate a good quality sequence, whereas primers KVL076 and KVL082 failed in, respectively, three

and one instance. Two failures were due to a shift in the fluorograms after 160, respectively, 220 covered nucleotides (KVL082 on AR07-2743 and KVL076 on AR06-1186). The primer KVL076 was unsuccessful for two samples that displayed a deletion of nine nucleotides at the primer-binding site (nucleotide string GAGGAAATG) (samples AR07-050 and AR07-1554). 3.2. Variability of the HIV-1 IN region The analysis of the 70 full-length IN sequences revealed that 110 positions of IN were polymorphic (110/288; 38%) (Table 3). The variability per site ranged from 1 until 12.7. The most variable positions were 136 (12.7 variability), 125 (9.8), 112 (9.5), 124 (8.8) and 119 (8.4). None of these high polymorphic positions are associated with INI resistance (Hazuda et al., 2007; McColl et al., 2007). Of the secondary resistance mutations against raltegravir and elvitegravir, only 151I and 157Q were observed in, respectively, four and five of our INI-naïve patients. In addition, other polymorphic changes at positions associated with resistance were observed, i.e. 74IV, 138D and 163AE. 4. Discussion The approval of raltegravir provides clinicians with the means of designing sufficiently potent regimens to attain long-term success for HIV-1 patients failing their current therapy with multi-drug resistant strains (Cooper et al., 2007; Steigbigel et al., 2007; Grinsztejn et al., 2007). Nevertheless, development of drug resistance against integrase inhibitors is expected, especially in patients in whom these new drugs cannot be supported by a potent backbone regimen due to broad cross-resistance in existing drug classes (Hazuda et al., 2007). The aim of this study was to develop a genotypic drug resistance test applicable to all HIV-1 group M subtypes. Therefore, amplification and sequencing primers were chosen at the most conserved sites surrounding and within the IN region. Within the primers, ambiguity positions were also included to take into account the genetic diversity at the primer binding sites. This approach resulted in a successful amplification of 70 genetically divergent samples belonging to the genetic forms A, B, C, D, F, G, J, CRF01-AE, CRF02-AG, CRF12-BF and CRF13-cpx. Their viral load ranged from 178 to >500,000 RNA copies/ml. However, for 4 samples (1 subtype C, 1 CRF01, 1 CRF02 and 1 CRF03) no amplification product could be obtained. These negative results could be due to the genetic differences between target and primer sequence, low viral load or the quality of the RNA extracts. Indeed, in the first attempt, when using old RNA extracts, only 64 samples could be amplified, while when using freshly extracted RNA, all but one of the failures could be resolved, indicating that the quality of the RNA is probably one of the most important determinants of success. For three of the remaining failures, a sample was not available to carry out a fresh extraction, and the fourth repeated failing sample had a low viral load (569 RNA copies/ml). The most likely rea-

Table 3 Natural polymorphisms within the HIV-1 integrase region Genetic form

Position in NTD of IN K7

Variability Genetic form

E10 D2

E11 D2 ; ED1 D4 ; ED1 D1 D2

E13

K14

S17

R20

A21

N1 C3 ; SN1 N6

K1 K1

T3 ; AT1 T1

D3

R5 R1 ; KR1 R2 ; KR1 R7 R5 ; KR1

KR1 D1 D1 D4

Q1

R1 Q1 3.2

S24

V2 ; AV1

H2 G1 ; N1 D1 ; G1 ; N1

D25 E1 ; DE1 E5

SN1

E1

N27

V1

P30

V31

H1

SN1

V32

I9 I1 I9 I7 I7

I4 ; L1

K34

V37

A38

KR1

I1 ; VI1 I2

D1

S39

D41

L45

N1 C1 ; SC1 C1 ; N1 N1 ; SN1 N2

N1

I1 V1 ; LV1 Q3

I1 R2

D1 2.1

T5

A23

R4 R3

2.8

2.1

3.7

N2 ; ST1 C1 ; N1 ; SN1 N3

K1

5.7

2.1

2.4

L74

V77

AT1

N2

G1

VI1

N1 2.1

6.1

2.3

3.1

S1

I1 I3

2.0

3.8

3.5

K103

A105

X1 T1

2.0

2.1

2.0

4.8

2.0

4.4

Position in CCD of IN M50

A9 B15 C10 CRF017 CRF028 CRF121 CRF133 D6 F3 G5 J3

I1 ; MI1 I4 ; T1 I2 ; T1 MI2 ; MIT2

Variability

5.2

S57 SG1

G59 GE1

I60

L63

V72

V1

I3 I6 ; VI2 I8 VI1 I8

M1 ; L1

I73

LV1 I1

I2

3.2

3.8

I1

2.0

2.0

3.1

IV1

A91

E96

F100

VI1

I1 I3

M5 ; MI1 ; ML1 M1 ; V1 IM1 L1

I1

IM1

Y5

AG1

D1

I1 I2

2.0

Y99

M1

L1 I3 I2 I2 I4 I1

I84

F2

3.5

2.0

4.9

2.0

2.0

2.1

2.1

L101 I1 I5 ; LI1 I9 I3 ; V1 I6 ; LI1 I1 I2 ; V1 I2 I2 I3 ; V1 I1 ; IV2

6.0

2.1

T122

T124

T125

V126

G134

I135

K136

E138

I141

I151

A7 ; TA1 ; GS1 A2 ; N3 ; TA2 ; DN1 N3 ; A2 ; S3 ; GS1 A7 A6 ; N2 A1 A3 A4 ; TA1 ; GS1 A2 A3 ; N2 A3

A7 ; S1 A1 ; MV1 A8 A7 A8

F1

N3 ; D2 ; S1

V4

Q2 ; R2 ; HQ1

ED1

V1

V9 V13 ; IV1 V10 V7 V8 V1 V1 V6 V3 V5 V3

Variability

2.4

8.8

9.8

A1 ; TA1 A1 ; S1 ; TA1 A1 ; V1 ; TA1 A5 A3

N2 ; X1 F1

2.1

ED1

IV1 D1 N5 D2

V3 V3

8.5

2.9

H1 ; Q2 T3 ; Q1 ; RT1 12.7

2.1

2.0

2.1

S119 P2 P4 ; T1 ; SG1 P1

GA1

R1 R3 ; KT1

3.8

8.4

V1

I1 I4 ; TI1

V4 V7

V113 I7 I15 I10 I6 ; VI1 I7 I1 I3 VI1 I2 ; L1 I4 ; L1 I1 ; VI2

R2

A9 B15 C10 CRF017 CRF028 CRF121 CRF133 D6 F3 G5 J3

N6 N5 ; DN1

T112 V7 ; A1 ; IV1 A1 ; TI1 V10 V7 V8

R1

Position in CCD of IN

TI1 I1 I1 I1

K111 T1 ; KR1 Q1

KR1

Genetic form

Q8 ; X1 Q4 ; R2 ; KQ1 T5 Q1 T2

G106 A1

2.0

K156 R1

2.1

E157

4.6

K160

Q1

E1

Q2

Q1

Q2

N2 ; KQ1

V1 V4 ; I1 ; M1 R1 ; IL1 V5 V2

9.5

SG1 P1 P1 ; X1 A1 ; R1 P1 P2

G163

V165

D167

E1

I1 I1

E7

I1

E6

K. Van Laethem et al. / Journal of Virological Methods 153 (2008) 176–181

A9 B15 C10 CRF017 CRF028 CRF121 CRF133 D6 F3 G5 J3

A1 E1 ; DE1 E2 I1 E1

2.0

2.2

4.3

2.9

E1 2.1

2.7

179

Table 3 (Continued ) 180

Genetic form

Position in CCD of IN H171

9

Q1

Variability

2.0

R R2 R3 ; KR1

G193

E198

I200

V201

V1

GE1 1

KR

3.2

1

X

M

KR1

2.1

A205

T206

L2 2.0

2.1

2.0

3.1

2.0

K219

I220

N222

3.2

D207

I208

T210

K211

1

I I4 I9 ; VI1 I7 I8 I1 I2 I6 I3 I5 I3

D1

1

IV1

2.3

I203

8

S S4 S2

1

M

DE1

TI1

S8 S1

S1 AS1

2.5

2.0

2.1

S1 S2 S5 S3

E1

3.2

2.1

L3 IM1

TA1

L1 3.2

KR1 KR1 R1 KR1 Q1 ; R1

R1 3.1

3.3

Position in CTD of IN K215 2

Q216

I217

1

T218 1

R ;N N2 KN2

1

1

S ; TI S1 I3 ; TI1

L L1 V1

N1

R224

1

Y227

S230

2

K K1

X1 QR1 A1 A1 3.1

3.4

2.0

2.0

I3 ; TI3 ; IL1 I1 I1 M1 ; TI1 ; TM1 I1 I2 ; A1 ; S1 I1

Q1

F F1

I L12 ; I3 I10 I7 I8

N1

7.8

4.3

K2 ; NH1

KQ1

L1

I3 I5 ; T1 I1 ; T1 I5 I3

Q1 YF1 F2

K1

E1 Q1 3.2

3.3

2.0

V234 8

IV1

2.2

2.0

5.3

K240

D253

1

N254 1

R

K K1 Q1

G2 ; N1 G1

ND1

N3

E1

G1 N1 N3 ; G1

KR1 2.1

S255

2.0

4.3

3.7

Position in CTD of IN

9

Variability

G189

R L1 ; Y1

R1

255-ins A B15 C10 CRF017 CRF028 CRF121 CRF133 D6 F3 G5 J3

K188

1

D256

V259

2

E E3 E3

A265 1

I1

AV V1 V4 ; AV2

I268

R269

D270

1

1

L

K273

M275

D278

1

A282

1

H ; DN 6

V281

8

M

2

8

K ; RK

1

A

P VM1

E1 E1 E2

I1

G1

V2 AT1

E3 E4 E3 –

2.9

V2

K2

V1 IL2

H3

Q1

N1 DH1

3.7

3.2

1.7

3.4

2.0

2.0

G G1 G3 ; SG1 G7 G7 G1 G2 ; SG1 G3 G3 G3

V1 2.1

S283

2.3

2.1

2.0

3.7

R284

D286

E287

D288

1

G

N1 G2 ; RG1

DN1

RG1 G1

N1

Q1

2.2

2.1

2.0

2.0

Wild-type amino acids and positions according to pNL4.3 (AF324493). The number of sequences belonging to a certain genetic form (subtype or CRF) or displaying a certain amino acid at a particular position is displayed in superscript. A mixture is indicated by a 2-letter code listing both possible amino acids, when containing more than 2 amino acids it is annotated as X. Amino acid changes and positions associated with raltegravir or elvitegravir resistance are in bold or underlined, respectively. Variability is defined as the number of observed amino acids at a certain position divided by the frequency of the most prevalent amino acid. 255-ins: insertion between amino acids 255 and 256. IN: integrase; NTD: amino-terminal domain; CCD: catalytic core domain; CTD: carboxyl-terminal domain.

K. Van Laethem et al. / Journal of Virological Methods 153 (2008) 176–181

Genetic form

K186

D1

9

Variability

I182

L1

E212 A B15 C10 CRF017 CRF028 CRF121 CRF133 D6 F3 G5 J3

F181

1

A B15 C10 CRF017 CRF028 CRF121 CRF133 D6 F3 G5 J3

Genetic form

K173

K. Van Laethem et al. / Journal of Virological Methods 153 (2008) 176–181

son for failure is thus probably not related to divergence of the subtype. Full-length sequences covered in both directions were obtained for almost all sequences. Intriguingly, one of the sequencing primers failed in two epidemiologically linked sequences due to a deletion of 9 nucleotides at the primer-binding site. This would result into a yet unreported deletion of three amino acids at the 3 prime end of RNase H (amino acids GGN at position 543–545). In the study of Hazuda et al. (2007) in which they investigated drug resistance development in patients failing raltegravir, the majority of patients failed with either 155H or 148HKR within IN. Each of these signature mutations were associated often with secondary resistance mutations, i.e. 74M, 92Q, 97A, 143HR, 151I, 163KR or 232N for 155H and 74M, 138AK or 140AS for 148HKR. In a phase II clinical trial of elvitegravir, another integrase inhibitor, the most common mutations were 66AIK, 92Q, 138K, 148HKR and 155H, suggesting cross-resistance between the first generation integrase inhibitors (McColl et al., 2007). In vitro selection experiments with elvitegravir revealed two distinct resistance pathways, i.e. 66I with 95K, 138K, 146P and 147G; and 92Q with 51Y, 147Q and 157Q (Shimura et al., 2008). In general, the signature mutations lead to a decrease in susceptibility towards the respective integrase inhibitors, whereas secondary resistance mutations augment the resistance level (Hazuda et al., 2007; Malet et al., 2008; Shimura et al., 2008). Within this study of IN sequences from INI-naïve patients, the majority of the positions associated with IN resistance were conserved. Only positions 74, 138, 151, 157 and 163 displayed some low-level variability (range 2.1–3.5). With the exception of 151I in four and 157Q in five patients, they all showed other amino acid changes than those reported previously (Hazuda et al., 2007; McColl et al., 2007; Shimura et al., 2008). The relevance of these secondary mutations in absence of signature mutations remains unclear. However, the single mutants 151I and 157Q revealed a rather highly reduced susceptibility towards raltegravir and elvitegravir in one particular in vitro study (Shimura et al., 2008). Additionally, from studies investigating resistance development against protease and reverse transcriptase inhibitors, it is known that the genetic background and the presence of secondary mutations can influence drug susceptibility, resistance development and clinical response (Perno et al., 2001; Abecasis et al., 2005, 2006; Deforche et al., 2006, 2007). In conclusion, this integrase assay displayed a high amplification and sequencing success rate and was shown to be a useful tool for the genotypic characterisation of diverse HIV-1 group M subtypes. Acknowledgements This work was supported by the AIDS Reference Laboratory of Leuven that receives support from the Belgian Ministry of Social Affairs through a fund within the Health Insurance System, by the Fonds voor Wetenschappelijk Onderzoek Vlaanderen (grant G.0266.04) and by the Centers of Excellence of the Katholieke Universiteit Leuven (krediet EF/05/015). Kris Covens was funded by a PhD grant of the Institute for the Promotion of Innovation through Sciences and Technology in Flanders (IWT). References Abecasis, A., Deforche, K., Snoeck, J., Bachelor, L.T., McKenna, P., Carvalho, A.P., Gomes, P., Camacho, R., Vandamme, A.-M., 2005. Protease mutation M89I/V is linked to therapy failure in patients infected with the HIV-1 non-B subtypes C, F or G. AIDS 19, 1799–1806.

181

Abecasis, A.B., Deforche, K., Bachelor, L.T., McKenna, P., Carvalho, A.P., Gomes, P., Vandamme, A.-M., Camacho, R.J., 2006. Investigation of baseline susceptibility to protease inhibitors in HIV-1 subtypes C, F G and CRF02 AG. Antivir. Ther. 11, 581–589. Cooper, D., Gatell, J., Rocktroh, J., Katlama, C., Yeni, P., Lazzarin, A., Chen, J., Isaaces, R., Teppler, H., Nguyen, B., and for the BENCHMRK-1 Study Group, 2007. Results of BENCHMRK-1, a phase III study evaluating the efficacy and safety of MK0518, a novel HIV-1 integrase inhibitor, in patients with triple-class resistant virus. In: Proceedings of the 15th Conference on Retroviruses and Opportunistic Infections, Abstract 105aLB, Boston, US. Deforche, K., Silander, T., Camacho, R., Grossman, Z., Soares, M.A., Van Laethem, K., Kantor, R., Moreau, Y., Vandamme, A.-M., 2006. Analysis of HIV-1 pol sequences using Bayesian Networks: implications for drug resistance. Bioinformatics 22, 2975–2979. Deforche, K., Camacho, R., Grossman, Z., Silander, T., Soares, M.A., Moreau, Y., Shafer, R.W., Van Laethem, K., Carvalho, A.P., Wynhoven, B., Cane, P., Snoeck, J., Clarke, J., Sirivichayakul, S., Ariyoshi, K., Holguin, A., Rudich, H., Rodrigues, R., Bouzas, M.B., Cahn, P., Brigido, L.F., Soriano, V., Sugiura, W., Phanuphak, P., Morris, L., Weber, J., Pillay, D., Tanuri, A., Harrigan, P.R., Shapiro, J.M., Katzenstein, D.A., Kantor, R., Vandamme, A.-M., 2007. Bayesian network analysis of resistance pathways against HIV-1 protease inhibitors. Infect. Genet. Evol. 7, 382–390. de Oliveira, T., Deforche, K., Cassol, S., Salminen, M., Paraskevis, D., Seebregts, C., Snoeck, J., van Rensburg, E.J., Wensing, A.M.J., van de Vijver, D.A., Boucher, C.A., Camacho, R., Vandamme, A.-M., 2005. An automated genotyping system for analysis of HIV-1 and other microbial sequences. Bioinformatics 21, 3797–3800, http://jose.med.kuleuven.be/subtypetool/html/index.html. DHHS panel on antiretroviral guidelines for adults and adolescents—a working group of the office of AIDS research advisory council (OARAC), 2007. Guidelines for the use of antiretroviral agents in HIV-1-infected adults and adolescents (01.12.07, http://AIDSinfo.nih.gov). Grinsztejn, B., Nguyen, B.-Y., Katlama, C., Gatell, J.M., Lazzarin, A., Vittecoq, D., Gonzalez, C.J., Chen, J., Harvey, C.M., Isaacs, R.D., for the Protocol 005 Team, 2007. Safety and efficacy of the HIV-1 integrase inhibitor raltegravir (MK-0518) in treatmentexperienced patients with multidrug-resistant virus: a phase II randomised controlled trial. Lancet 369, 1261–1269. Hazuda, D.J., Miller, M.D., Nguyen, B.Y., Zhao, J., for the P005 Study Team, 2007. Resistance to the HIV-integrase inhibitor raltegravir: analysis of protocol 005, a phase II study in patients with triple-class resistant HIV-1 infection. Antivir. Ther. 12, S10. Leitner, T., Foley, B., Hahn, B., Marx, P., McCutchan, F., Mellors, J., Wolinsky, S., Korber, B. (Eds.), 2005. HIV Sequence Compendium. Published by Theoretical Biology and Biophysics Group, Los Alamos National Laboratory, NM, LA-UR 06–0680. Malet, I., Delelis, O., Valantin, M.-A., Montes, B., Soulie, C., Wirden, M., Tchertanov, L., Peytavin, G., Reynes, J., Mouscadet, J.-F., Katlama, C., Calvez, V., Marcelin, A.-G., 2008. Mutations associated with failure to raltegravir treatment affect integrase sensitivity to the inhibitor in vitro. Antimicrob. Agents Chemother. 52, 1351–1358. McColl, D.J., Fransen, S., Gupta, S., Parkin, N., Margot, N., Chuck, S., Cheng, S., Miller, M.D., 2007. Resistance and cross-resistance to first generation integrase inhibitors: insights from a phase II study of elvitegravir (GS-9137). Antivir. Ther. 12, S11. Perno, C.F., Cozzi-Lepri, A., Balotta, C., Forbici, F., Violin, M., Bertoli, A., Facchi, G., Pezzotti, P., Cadeo, G., Tositti, G., Pasquinucci, S., Pauluzzi, S., Scalzini, A., Salassa, B., Vincenti, A., Phillips, A.N., Dianzani, F., Appice, A., Angarano, G., Monno, L., Ippolito, G., Moroni, M., d’Arminio Monforte, A., the Italian Cohort Naive Antiretroviral (I.C.O.N.A.) Study Group, 2001. Secondary mutations in the protease region of human immunodeficiency virus and virologic failure in drugnaive patients treated with protease inhibitor-based therapy. J. Infect. Dis. 184, 983–991. Shimura, K., Kodama, E., Sakagami, Y., Matsuzaki, Y., Watanabe, W., Yamataka, K., Watanabe, Y., Ohata, Y., Doi, S., Sato, M., Kano, M., Ikeda, S., Matsuoka, M., 2008. Broad antiretroviral activity and resistance profile of the novel human immunodeficiency virus integrase inhibitor elvitegravir (JTK-303/GS-9137). J. Virol. 82, 764–774. Snoeck, J., Van Laethem, K., Hermans, P., Van Wijngaerden, E., Derdelinckx, I., Schrooten, Y., van de Vijver, D.A., De Wit, S., Clumeck, N., Vandamme, A.-M., 2004. Rising prevalence of HIV-1 non-B subtypes in Belgium: 1983–2001. J. Acquir. Immune Defic. Syndr. 35, 279–285. Steigbigel, R., Kumar, P., Eron, J., Schechter, M., Markowitz, M., Loufty, M., Zhao, J., Isaacs, R., Nguyen, B., Teppler, H., and for the BENCHMRK-2 Study Group, 2007. Results of BENCHMRK-2, a phase III study evaluating the efficacy and safety of MK-0518, a novel HIV-1 integrase inhibitor, in patients with triple-class resistant virus. In: Proceedings of the 15th Conference on Retroviruses and Opportunistic Infections, Abstract 105bLB, Boston, US. Vercauteren, J., Derdelinckx, I., Sasse, A., Bogaert, M., Ceunen, H., De Roo, A., De Wit, S., Deforche, K., Echahidi, F., Fransen, K., Goffard, J.C., Goubau, P., Goudeseune, E., Yombi, J.C., Lacor, P., Liesnard, C., Moutschen, M., Pierard, D., Rens, R., Schrooten, Y., Vaira, D., van den Heuvel, A., van der Gucht, B., Van Ranst, M., Van Wijngaerden, E., Vandercam, B., Vekemans, M., Verhofstede, C., Clumeck, N., Vandamme, A.-M., Van Laethem, K., 2008. Prevalence and epidemiology of HIV type 1 drug resistance among newly diagnosed therapy-naive patients in Belgium from 2003 to 2006. AIDS Res. Hum. Retroviruses 24, 355–362.