Rapid and Accurate Identification of SARS-CoV-2 Variants Using Real Time PCR Assays

Background

Genomic surveillance efforts for SARS-CoV-2 are needed to understand the epidemiology of the COVID-19 pandemic. Viral variants may impact routine diagnostic testing, increase viral transmissibility, cause differences in disease severity, have decreased susceptibility to therapeutics, and/or confer the ability to evade host immunity. While viral whole-genome sequencing (WGS) has played a leading role in surveillance programs, many laboratories lack the expertise and resources for performing WGS. This study describes the performance of multiplexed real-time reverse transcription-PCR (RT-PCR) assays for identification of SARS-CoV-2 variants.

Methods

SARS-CoV-2 specimens were tested for spike-gene variants using a combination of allele-specific primer and allele-specific detection technology (PlexPrime^® and PlexZyme^®). Targeted detection of spike gene mutations by RT-PCR was compared to variant detection in positive specimens by WGS, including the recently emerged SARS-CoV-2 Omicron variant.

Results

A total of 398 SAR-CoV-2 RT-PCR positive and 39 negative specimens previously tested by WGS were re-tested by RT-PCR genotyping. PCR detection of spike gene mutations N501Y, E484K, and S982A correlated 100% with WGS for the 29 lineages represented, including Alpha (B.1.1.7), Beta (B.1.351), and Gamma (P.1). Incorporating the P681R spike gene mutation also allowed screening for the SARS-CoV-2 Delta variant (B.1.617.2 and AY sublineages). Further sampling of 664 specimens that were screened by WGS between June and August 2021 and then re-tested by RT-PCR showed strong agreement for Delta variant positivity: 34.5% for WGS vs 32.9% for RT-PCR in June; 100% vs 97.8% in August. In a blinded panel of 16 Omicron and 16 Delta specimens, results of RT-PCR were 100% concordant with WGS results.

Conclusions

These data demonstrate that multiplexed real-time RT-PCR genotyping has strong agreement with WGS and may provide additional SARS-CoV-2 variant screening capabilities when WGS is unavailable or cost-prohibitive. RT-PCR genotyping assays may also supplement existing sequencing efforts while providing rapid results at or near the time of diagnosis to help guide patient management.

Development and validation of a high throughput SARS-CoV-2 whole genome sequencing workflow in a clinical laboratory

Monitoring new mutations in SARS-CoV-2 provides crucial information for identifying diagnostic and therapeutic targets and important insights to achieve a more effective COVID-19 control strategy. Next generation sequencing (NGS) technologies have been widely used for whole genome sequencing (WGS) of SARS-CoV-2. While various NGS methods have been reported, one chief limitation has been the complexity of the workflow, limiting the scalability. Here, we overcome this limitation by designing a laboratory workflow optimized for high-throughput studies. The workflow utilizes modified ARTIC network v3 primers for SARS-CoV-2 whole genome amplification. NGS libraries were prepared by a 2-step PCR method, similar to a previously reported tailed PCR method, with further optimizations to improve amplicon balance, to minimize amplicon dropout for viral genomes harboring primer-binding site mutation(s), and to integrate robotic liquid handlers. Validation studies demonstrated that the optimized workflow can process up to 2688 samples in a single sequencing run without compromising sensitivity and accuracy and with fewer amplicon dropout events compared to the standard ARTIC protocol. We additionally report results for over 65,000 SARS-CoV-2 whole genome sequences from clinical specimens collected in the United States between January and September of 2021, as part of an ongoing national genomics surveillance effort.

Performance of Unobserved Self-Collected Nasal Swabs for Detection of SARS-CoV-2 by RT-PCR Utilizing a Remote Specimen Collection Strategy

Background

The use of a remote specimen collection strategy employing a kit designed for unobserved self-collection for severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) reverse transcription polymerase chain reaction (RT-PCR) can decrease the use of personal protective equipment (PPE) and exposure risk. To assess the impact of unobserved specimen self-collection on test performance, we examined results from a SARS-CoV-2 qualitative RT-PCR test for self-collected specimens from participants in a return-to-work screening program and assessed the impact of a pooled testing strategy in this cohort.

Methods

Self-collected anterior nasal swabs from employee return-to-work programs were tested using the Quest Diagnostics Emergency Use Authorization SARS-CoV-2 RT-PCR. The cycle threshold (Ct) values for the N1 and N3 N-gene targets and a human RNase P (RP) gene control target were tabulated. For comparison, we utilized Ct values from a cohort of health care provider–collected specimens from patients with and without coronavirus disease 2019 symptoms.

Results

Among 47923 participants, 1.8% were positive. RP failed to amplify for 13/115435 (0.011%) specimens. The median (interquartile range) Cts were 32.7 (25.0–35.7) for N1 and 31.3 (23.8–34.2) for N3. Median Ct values in the self-collected cohort were significantly higher than those of symptomatic but not asymptomatic patients. Based on Ct values, pooled testing with 4 specimens would have yielded inconclusive results in 67/1268 (5.2%) specimens but only a single false-negative result.

Conclusions

Unobserved self-collection of nasal swabs provides adequate sampling for SARS-CoV-2 RT-PCR testing. These findings alleviate concerns of increased false negatives in this context. Specimen pooling could be used for this population, as the likelihood of false-negative results is very low when using a sensitive, dual-target methodology.

Pooling of Upper Respiratory Specimens Using a SARS-CoV-2 Real-time RT-PCR Assay Authorized for Emergency Use in Low-Prevalence Populations for High-Throughput Testing

Background

Nucleic acid amplification testing is a critical tool for addressing the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) pandemic. Specimen pooling can increase throughput and conserve testing resources but requires validation to ensure that reduced sensitivity does not increase the false-negative rate. We evaluated the performance of a real-time reverse transcription polymerase chain reaction (RT-PCR) test authorized by the US Food and Drug Administration (FDA) for emergency use for pooled testing of upper respiratory specimens.

Methods

Positive specimens were selected from 3 prevalence groups, 1%–3%, >3%–6%, and >6%–10%. Positive percent agreement (PPA) was assessed by pooling single-positive specimens with 3 negative specimens; performance was assessed using Passing-Bablok regression. Additionally, we assessed the distributions of RT-PCR cycle threshold (Ct) values for 3091 positive specimens.

Results

PPA was 100% for the 101 pooled specimens. There was a linear relationship between Ct values for pooled and single-tested specimens (r = 0.96–0.99; slope ≈ 1). The mean pooled Ct shifts at 40 cycles were 2.38 and 1.90, respectively, for the N1 and N3 targets. The median Cts for 3091 positive specimens were 25.9 (N1) and 24.7 (N3). The percentage of positive specimens with Cts between 40 and the shifted Ct was 1.42% (N1) and 0.0% (N3).

Conclusions

Pooled and individual testing of specimens positive for SARS-CoV-2 demonstrated 100% agreement, which demonstrates the viability of pooled specimens for SARS-COV-2 testing using a dual-target RT-PCR system. Pooled specimen testing can help increase testing capacity for SARS-CoV-2 with a low risk of false-negative results.

Trends in HIV-1 Drug Resistance Mutations from a U.S. Reference Laboratory from 2006 to 2017

Trends in resistance to antiretroviral drugs for HIV-1 may inform clinical support and drug development. We evaluated drug resistance mutation (DRM) trends for nucleoside reverse transcriptase inhibitor (NRTI), non-nucleoside reverse transcriptase inhibitor (NNRTI), protease inhibitor (PI), and integrase strand transfer inhibitor (INSTI) in a large U.S. reference laboratory database. DRMs with a Stanford HIV Drug Resistance Database mutation score ≥10 from deidentified subtype B NRTI/NNRTI/PI specimens (2006-2017; >10,000/year) and INSTI specimens (2010-2017; >1,000/year) were evaluated. Sequences with NRTI, NNRTI, or PI single- or multiclass DRMs declined from 48.9% to 39.3%. High-level dual- and triple-class resistance declined from 43.3% (2006) to 17.1% (2017), while sequences with only single-class DRMs increased from 40.0% to 52.9%. The prevalence of DRMs associated with earlier treatment regimens declined, while prevalence of some DRMs associated with newer regimens increased. M184V/I decreased from 48.3% to 29.4%. K103N/S/T declined from 42.5% in 2012 to 36.4% in 2017. Rilpivirine and etravirine DRMs E138A/Q/R and E138K increased from 4.9% and 0.4% to 9.7% and 1.7%, respectively. Sequences with ≥1 darunavir DRM declined from 18.1% to 4.8% by 2017. INSTI DRM Q148H/K/R declined from 39.3% (2010) to 13.8% (2017). Prevalence of elvitegravir-associated DRMs T66A/I/K, E92Q, S147G, and the dolutegravir-associated DRM R263K increased. For a subset of patients with serial testing, 50% (2,646/5,290) of those who initially had no reportable DRM subsequently developed ≥1 DRM for NRTI/NNRTI/PI and 49.7% (159/320) for INSTI. These trends may inform the need for baseline genotypic resistance testing. The detection of treatment-emergent DRMs in serially tested patients confirms the value of genotypic testing following virologic failure.

Detailed Transmission Network Analysis of a Large Opiate-Driven Outbreak of HIV Infection in the United States

In January 2015, an outbreak of undiagnosed human immunodeficiency virus (HIV) infections among persons who inject drugs (PWID) was recognized in rural Indiana. By September 2016, 205 persons in this community of approximately 4400 had received a diagnosis of HIV infection. We report results of new approaches to analyzing epidemiologic and laboratory data to understand transmission during this outbreak. HIV genetic distances were calculated using the polymerase region. Networks were generated using data about reported high-risk contacts, viral genetic similarity, and their most parsimonious combinations. Sample collection dates and recency assay results were used to infer dates of infection. Epidemiologic and laboratory data each generated large and dense networks. Integration of these data revealed subgroups with epidemiologic and genetic commonalities, one of which appeared to contain the earliest infections. Predicted infection dates suggest that transmission began in 2011, underwent explosive growth in mid-2014, and slowed after the declaration of a public health emergency. Results from this phylodynamic analysis suggest that the majority of infections had likely already occurred when the investigation began and that early transmission may have been associated with sexual activity and injection drug use. Early and sustained efforts are needed to detect infections and prevent or interrupt rapid transmission within networks of uninfected PWID.

Antiretroviral Drug Resistance Among Children and Youth in the United States With Perinatal HIV

Among 234 US youths with perinatal human immunodeficiency virus, 75% had antiretroviral resistance, substantially higher than that of the reference laboratory overall (36%-44%). Resistance to newer antiretrovirals and to all antiretrovirals in a class was uncommon. The only factor independently associated with future resistance was a higher peak viral load.

A genotypic HIV-1 proviral DNA coreceptor tropism assay: characterization in viremic subjects

BACKGROUND

HIV-1 coreceptor tropism testing is used to evaluate eligibility for CCR5 antagonist therapy. However, HIV-1 RNA-based tests are not suitable for virologically suppressed patients, therefore the use of proviral DNA tropism testing has been investigated. We describe a novel proviral DNA-based genotypic tropism assay and compare its performance to that of a sensitive HIV-1 RNA-based genotypic test.

METHODS

Tropism was determined using HIV-1 plasma RNA and proviral DNA from 42 paired samples from patients with plasma viral loads ≥1000 HIV-1 RNA copies/mL. Proviral DNA sample types included whole blood, separated peripheral blood mononuclear cells resuspended in phosphate-buffered saline and peripheral blood mononuclear cells resuspended in spun plasma. The HIV-1 envelope V3 region was PCR-amplified, sequenced in triplicate, and analyzed for tropism with the geno2pheno algorithm using a 10% false-positive rate (FPR).

RESULTS

Amplicons were obtained from proviral DNA and plasma RNA in 41/42 samples. Tropism predictions were highly concordant (93%-98%) between proviral DNA and plasma RNA, regardless of the proviral DNA isolation method. Non-R5 proviral DNA results were obtained for 100% of patients with detectable non-R5 plasma HIV-1 RNA results. Geno2pheno FPRs for proviral DNA and plasma RNA were highly correlated (Spearman rho = 0.86).

CONCLUSIONS

Our findings demonstrate that proviral DNA tropism determinations from whole blood or peripheral blood mononuclear cells were highly concordant with plasma HIV-1 RNA tropism determinations. This assay may be useful for screening virologically suppressed patients for CCR5-antagonist eligibility and for research purposes.

Comparison of genotypic and phenotypic HIV type 1 tropism assay: results from the screening samples of Cenicriviroc Study 202, a randomized phase II trial in treatment-naive subjects

Cenicriviroc is a once-daily oral CCR5/CCR2 antagonist in development for treatment of HIV infection. CVC Study 202 (652-2-202; NCT01338883) excluded treatment-naive subjects demonstrated to harbor non-R5 (CXCR4-tropic or dual-mixed) tropic HIV-1 by either genotypic or phenotypic tropism testing. Here we compare the results of genotypic and phenotypic tropism testing in Study 202. A total of 304 subjects screened had paired genotypic and phenotypic results. Genotypic tropism testing (GTT) incorporated triplicate population sequencing using the geno2pheno algorithm and the PSSM algorithm, followed by ultradeep sequencing (UDS) for samples with R5 results. All samples were further evaluated with a phenotypic test, the enhanced-sensitivity Trofile assay (ESTA). Concordance between GTT and ESTA was 80% and increased to 84% when only geno2pheno was used for triplicate population sequencing. GTT (geno2pheno) classified 18% of the samples as non-R5 compared to 16% by ESTA. Only one-third of samples with non-R5 results by either test were classified as non-R5 by both tests. Median CD4((+)) cell counts were lower in patients with concordant non-R5 results by UDS and ESTA than in subjects with an R5 result by either assay (p=0.0004). UDS detected non-R5 virus in an additional 27/304 subjects (median 15% non-R5, interquartile range: 3.7-62%) with R5 results by ESTA. In conclusion, the geno2pheno algorithm improves concordance of GTT with a clinically validated phenotypic tropism assay as does the use of UDS. These findings provide support for recent guidelines indicating that genotypic tropism testing may be considered as an alternative to phenotypic testing.

Hepatitis B virus reverse transcriptase sequence variant database for sequence analysis and mutation discovery

Drug resistance resulting from reverse transcriptase (RT) mutations is one of the main obstacles to successful hepatitis B virus (HBV) therapy. Indeed, HBV treatment guidelines recommend HBV genotypic resistance testing for patients receiving nucleos(t)ide RT inhibitors (N(t)RTIs) who develop virological failure. N(t)RTI-resistance mutations at 10 RT positions have been well characterized in phenotypic studies, however, data are lacking on the relative frequency of these mutations in N(t)RTI-treated and untreated individuals. There are also few published data on the extent of amino acid variation at most of the 344 positions of HBV RT and the extent to which this variation is influenced by N(t)RTI treatment. We retrieved 23,871 HBV RT sequences from GenBank and reviewed the published reports of these sequences to ascertain the number of individuals from whom the sequences were obtained, the N(t)RTI treatments of these individuals, and the year and region of virus sampling. We then used these data to populate a relational database we named HBVrtDB. As of July 2010, HBVrtDB contained 6811 sequences from 3869 individuals reported in 281 references. Among these 3869 individuals, 73% were N(t)RTI-naïve and 27% received one or more N(t)RTIs. Among the 10 well-characterized N(t)RTI-resistance mutations, L80I/V, V173L, L180M, A181T, T184S, S202G and M204I/V were significantly associated with treatment with lamivudine, an l-nucleoside analog, and A181S/T/V and N236T were significantly associated with treatment with adefovir, an acyclic nucleoside phosphonate. A similar analysis of ten additional less well-characterized resistance mutations demonstrated a significant association with N(t)RTI treatment for four of the mutations: L82M, S85A, A200V, and Q215S. We also created an interactive program, HBVseq, to enable users to identify mutations in submitted sequences and retrieve the prevalence of these mutations in HBVrtDB according to genotype and N(t)RTI treatment. HBVrtDB and HBVseq are available at http://hivdb.stanford.edu/HBV/releaseNotes/.

Enzymatic and structural analysis of the I47A mutation contributing to the reduced susceptibility to HIV protease inhibitor lopinavir

Lopinavir (LPV) is a second-generation HIV protease inhibitor (PI) designed to overcome resistance development in patients undergoing long-term antiviral therapy. The mutation of isoleucine at position 47 of the HIV protease (PR) to alanine is associated with a high level of resistance to LPV. In this study, we show that recombinant PR containing a single I47A substitution has the inhibition constant (K(i) ) value for lopinavir by two orders of magnitude higher than for the wild-type PR. The addition of the I47A substitution to the background of a multiply mutated PR species from an AIDS patient showed a three-order-of-magnitude increase in K(i) in vitro relative to the patient PR without the I47A mutation. The crystal structure of I47A PR in complex with LPV showed the loss of van der Waals interactions in the S2/S2' subsites. This is caused by the loss of three side-chain methyl groups due to the I47A substitution and by structural changes in the A47 main chain that lead to structural changes in the flap antiparallel beta-strand. Furthermore, we analyzed possible interaction of the I47A mutation with secondary mutations V32I and I54V. We show that both mutations in combination with I47A synergistically increase the relative resistance to LPV in vitro. The crystal structure of the I47A/I54V PR double mutant in complex with LPV shows that the I54V mutation leads to a compaction of the flap, and molecular modeling suggests that the introduction of the I54V mutation indirectly affects the strain of the bound inhibitor in the PR binding cleft.

Identification of a novel resistance (E40F) and compensatory (K43E) substitution in HIV-1 reverse transcriptase

Background

HIV-1 nucleoside reverse transcriptase inhibitors (NRTIs) have been used in the clinic for over twenty years. Interestingly, the complete resistance pattern to this class has not been fully elucidated. Novel mutations in RT appearing during treatment failure are still being identified. To unravel the role of two of these newly identified changes, E40F and K43E, we investigated their effect on viral drug susceptibility and replicative capacity.

Results

A large database (Quest Diagnostics database) was analysed to determine the associations of the E40F and K43E changes with known resistance mutations. Both amino acid changes are strongly associated with the well known NRTI-resistance mutations M41L, L210W and T215Y. In addition, a strong positive association between these changes themselves was observed. A panel of recombinant viruses was generated by site-directed mutagenesis and phenotypically analysed. To determine the effect on replication capacity, competition and in vitro evolution experiments were performed. Introduction of E40F results in an increase in Zidovudine resistance ranging from nine to fourteen fold depending on the RT background and at the same time confers a decrease in viral replication capacity. The K43E change does not decrease the susceptibility to Zidovudine but increases viral replication capacity, when combined with E40F, demonstrating a compensatory role for this codon change.

Conclusion

In conclusion, we have identified a novel resistance (E40F) and compensatory (K43E) change in HIV-1 RT. Further research is indicated to analyse the clinical importance of these changes.

Molecular basis of antagonism between K70E and K65R tenofovir-associated mutations in HIV-1 reverse transcriptase

The K70E mutation in HIV-1 reverse transcriptase was observed in 10% of virologic non-responders of the abacavir/lamivudine/tenofovir arm of ESS30009, alone, or in mixtures with K65R by population sequencing. Clonal analysis of six ESS30009 K70E isolates failed to identify double mutants carrying K65R+K70E. Site-directed K70E mutants had a replication capacity of 97+/-29%, but only 2.4+/-0.9% for K65R+K70E and 0.01% for K65R+K70E+M184V mutants. K65R+K70E phenotypic fold changes for abacavir, lamivudine and tenofovir were comparable to reported values for K65R alone. In molecular dynamic simulations, the epsilon-amino group of K65 was positioned 2.7+/-0.1A from the gamma-phosphate of the dTTP ligand and stabilized the triphosphate. In the R65 mutant, this distance increased to 4.2+/-0.4A and the interaction energy with the ligand was less favorable, but the K70 epsilon-amino group was repositioned closer to the gamma-phosphate and had a more favorable interaction energy. In the double mutant, E70 could not stabilize the gamma-phosphate, resulting in a more severe defect. The net effect of the atomic-level changes in the double mutant may be to destabilize the pyrophosphate leaving group of the ligand, more severely affecting the catalytic rate of the polymerization reaction than the R65 single mutation.

The probable source of both the primary multidrug-resistant (MDR) HIV-1 strain found in a patient with rapid progression to AIDS and a second recombinant MDR strain found in a chronically HIV-1-infected patient

BACKGROUND

Rapid progression to AIDS after acute infection with a multidrug-resistant (MDR), dual-tropic strain of human immunodeficiency virus type 1 (HIV-1) was reported in a New York City man (hereafter referred to as "NYC") who has sex with men. The probable source of this HIV-1 (hereafter referred to as "CT01") and the development of a recombinant MDR HIV-1 in the source's partner (hereafter referred to as "CT02") are described.

METHODS

After identification of the epidemiological link of CT01 and CT02 to NYC, viral sequences and phenotypic analyses were compared. Confirmatory genotypic and phenotypic analyses, replicative capacity, and viral coreceptor use were assessed. Viral recombination was assessed using a sliding window technique and phylogenetic tree analysis.

RESULTS

NYC and CT01 were linked historically and epidemiologically and were genetically confirmed from CT01's samples acquired 2 days before and subsequent to the transmission event. Genotypic, recombination, and phylogenetic analyses suggest CT02 became superinfected by CT01 with subsequent production of a recombinant panresistant HIV-1.

CONCLUSION

The probable source of a dual-tropic, MDR HIV-1 that was associated with rapid progression to AIDS is illustrated, suggesting progression was not explained by the HIV-1 variant alone. A probable second finding of a chronically infected host becoming superinfected with MDR HIV-1 with subsequent formation of a panresistant recombinant HIV-1 is described. This case illustrates the public health implications of unsafe sex between serodiscordant and seroconcordant partners.

Increasing prevalence of HIV-1 protease inhibitor-associated mutations correlates with long-term non-suppressive protease inhibitor treatment

Treatment of human immunodeficiency virus type 1 with protease inhibitors (PIs) is associated with the emergence of resistance-associated mutations. Treatment-characterized datasets have been used to identify novel treatment-associated protease mutations. In this study, we utilized two large reference laboratory databases (>115,000 viral sequences) to identify non-established resistance-associated protease mutations. We found 20 non-established protease mutations occurring in 82% of viruses with a PI resistance score of 4-7, 62% of viruses with a resistance score of 1-3, and 35% of viruses with no predicted PI resistance. We correlated mutational prevalence to treatment duration in a treatment-characterized dataset of 2161 patients undergoing non-suppressive PI therapy. In the non-suppressed dataset, 24 mutations became more prevalent and three mutations became less prevalent after more than 48 months of non-suppressive PI-therapy. Longer durations of non-suppressive treatment correlated with higher PI resistance scores. Mutations at eight non-established positions that were more common in viruses with the longest duration of non-suppressive therapy were also more common in viruses with the highest PI resistance score. Covariation analysis of 3036 protease amino acid substitutions identified 75 positive and nine negative correlations between resistance associated positions. Our findings support the utility of reference laboratory datasets for surveillance of mutation prevalence and covariation.

Structural analysis of an HIV-1 protease I47A mutant resistant to the protease inhibitor lopinavir

We have identified a rare HIV-1 protease (PR) mutation, I47A, associated with a high level of resistance to the protease inhibitor lopinavir (LPV) and with hypersusceptibility to the protease inhibitor saquinavir (SQV). The I47A mutation was found in 99 of 112,198 clinical specimens genotyped after LPV became available in late 2000, but in none of 24,426 clinical samples genotyped from 1998 to October 2000. Phenotypic data obtained for five I47A mutants showed unexpected resistance to LPV (86- to >110-fold) and hypersusceptibility to SQV (0.1- to 0.7-fold). Molecular modeling and energy calculations for these mutants using our structural phenotyping methodology showed an increase in the binding energy of LPV by 1.9-3.1 kcal/mol with respect to the wild type complex, corresponding to a 20- to >100-fold decrease in binding affinity, consistent with the observed high levels of LPV resistance. In the WT PR-LPV complex, the Ile 47 side chain is positioned close to the phenoxyacetyl moiety of LPV and its van der Waals interactions contribute significantly to the ligand binding. These interactions are lost for the smaller Ala 47 residue. Calculated binding energy changes for SQV ranged from -0.4 to -1.2 kcal/mol. In the mutant I47A PR-SQV complexes, the PR flaps are packed more tightly around SQV than in the WT complex, resulting in the formation of additional hydrogen bonds that increase binding affinity of SQV consistent with phenotypic hypersusceptibility. The emergence of mutations at PR residue 47 strongly correlates with increasing prescriptions of LPV (Spearman correlation r(s) = 0.96, P < .0001).

New two-amino acid insertion near codon 70 of the HIV type 1 protease gene

A two-amino acid insertion near codon 70 of the HIV-1 protease gene was found in an individual failing protease inhibitor therapy. Susceptibility of this strain to protease inhibitors was similar to that of non-insert-containing strains with comparable resistance mutations and one or more other major PI mutations.

Rare One and Two Amino Acid Inserts Adjacent to Codon 103 of the HIV-1 Reverse Transcriptase (Rt) Affect Susceptibility to Non-Nucleoside Rt Inhibitors

HIV-1 strains that possess a one or two amino acid insert between codons 102 and 103 of the reverse transcriptase (RT) gene were identified in three HIV-1-infected individuals. Each strain also had one or more known mutations associated with nucleoside RT inhibitors (NRTIs) and non-nucleoside RT inhibitors (NNRTIs). Recombinant viruses from these strains had reduced susceptibility to efavirenz and nevirapine, and homology modelling predicted a loss of binding contacts with efavirenz. Mutagenesis studies indicated that replication of insert-containing strains was dependent on RT gene mutations and polymorphisms that co-evolved with the insert. These results suggest that inserts in the NNRTI-binding pocket contribute to NNRTI resistance, but are tolerated only under specific genetic conditions.

Rare one and two amino acid inserts adjacent to codon 103 of the HIV-1 reverse transcriptase (RT) affect susceptibility to non-nucleoside RT inhibitors

HIV-1 strains that possess a one or two amino acid insert between codons 102 and 103 of the reverse transcriptase (RT) gene were identified in three HIV-1-infected individuals. Each strain also had one or more known mutations associated with nucleoside RT inhibitors (NRTIs) and non-nucleoside RT inhibitors (NNRTIs). Recombinant viruses from these strains had reduced susceptibility to efavirenz and nevirapine, and homology modelling predicted a loss of binding contacts with efavirenz. Mutagenesis studies indicated that replication of insert-containing strains was dependent on RT gene mutations and polymorphisms that co-evolved with the insert. These results suggest that inserts in the NNRTI-binding pocket contribute to NNRTI resistance, but are tolerated only under specific genetic conditions.

Structure-based phenotyping predicts HIV-1 protease inhibitor resistance

Mutations in HIV-1 drug targets lead to resistance and consequent therapeutic failure of antiretroviral drugs. Phenotypic resistance assays are time-consuming and costly, and genotypic rules-based interpretations may fail to predict the effects of multiple mutations. We have developed a computational procedure that rapidly evaluates changes in the binding energy of inhibitors to mutant HIV-1 PR variants. Models of WT complexes were produced from crystal structures. Mutant complexes were built by amino acid substitutions in the WT complexes with subsequent energy minimization of the ligand and PR binding site residues. Accuracy of the models was confirmed by comparison with available crystal structures and by prediction of known resistance-related mutations. PR variants from clinical isolates were modeled in complex with six FDA-approved PIs, and changes in the binding energy (DeltaE(bind)) of mutant versus WT complexes were correlated with the ratios of phenotypic 50% inhibitory concentration (IC(50)) values. The calculated DeltaE(bind) of five PIs showed significant correlations (R(2) = 0.7-0.8) with IC(50) ratios from the Virco Antivirogram assay, and the DeltaE(bind) of six PIs showed good correlation (R(2) = 0.76-0.85) with IC(50) ratios from the Virologic PhenoSense assay. DeltaE(bind) cutoffs corresponding to a four-fold increase in IC(50) were used to define the structure-based phenotype as susceptible, resistant, or equivocal. Blind predictions for 78 PR variants gave overall agreement of 92% (kappa = 0.756) and 86% (kappa = 0.666) with PhenoSense and Antivirogram phenotypes, respectively. The structural phenotyping predicted drug resistance of clinical HIV-1 PR variants with an accuracy approaching that of frequently used cell-based phenotypic assays.

The use of data mining to investigate a possible quality problem with ultrasensitive HIV viral load data at a large reference laboratory

Suppression of HIV viral load to <50 copies/mL, the lower limit of detection for the ultrasensitive assays, has been shown to correlate with favorable clinical outcome. Patients periodically exhibited transient or sustained low-level viremia based on this test. In order to investigate a possible quality concern, we used our corporate data warehouse to examine the patterns in our data over time as well as across geographic regions.

Functional correlates of insertion mutations in the protease gene of human immunodeficiency virus type 1 isolates from patients

Twenty-four of over 24,000 patients genotyped over the past 3 years were found to have human immunodeficiency virus (HIV) isolates that possess an insert in the protease gene. In this report, we evaluated the spectrum of protease gene insertion mutations in patient isolates and analyzed the effect of these various insertion mutations on viral phenotypes. The inserts were composed of 1, 2, 5, or 6 amino acids that mapped at or between codons 35 and 38, 17 and 18, 21 and 25, or 95 and 96. Reduced susceptibility to protease inhibitors was found in isolates which possess previously reported drug resistance mutations. Fitness assays, including replication and competition experiments, showed that most of the isolates with inserts grew somewhat better than their counterparts with a deletion of the insert. These experiments demonstrate that, rarely, insertion mutations can develop in the HIV type 1 protease gene, are no more resistant than any other sequences which have similar associated resistance mutations, and can provide a borderline advantage in replication.

Targeted gene disruption of the Caenorhabditis elegans L-isoaspartyl protein repair methyltransferase impairs survival of dauer stage nematodes

The methylation of abnormal L-isoaspartyl residues by protein L-isoaspartate (D-aspartate) O-methyltransferase (EC 2.1.1.77) can lead to their conversion to L-aspartyl residues. For polypeptides damaged by spontaneous reactions that generate L-isoaspartyl residues, these steps represent a protein repair pathway that can limit the accumulation of potentially detrimental proteins in the aging process. We report here the construction and the characterization of an animal model deficient in this methyltransferase. We utilized Tc1-transposon-mediated mutagenesis in the nematode Caenorhabditis elegans to construct a homozygous excision mutant lacking exons 2-5 of the pcm-1 gene encoding this enzyme. Nematodes carrying this deletion exhibited no detectable L-isoaspartyl methyltransferase activity. These worms demonstrated normal morphology and behavior and adult mutant nematodes exhibited a normal lifespan. However, the survival of dauer-phase mutants was diminished by 3.5-fold relative to wild-type dauers after 50 days in the dauer phase. The fitness of the pcm-1 deletion nematodes was reduced by about 16% relative to that of wild-type nematodes as measured by the ability of these mutants to compete reproductively against a wild-type population. We found that the absence of the functional methyltransferase gene leads to a modest accumulation of altered protein substrates in aged dauer worms. However, in the viable fraction of these dauer worms, no differences were seen in the levels of altered substrate proteins in the parent and methyltransferase-deficient worms, suggesting that the enzyme in wild-type cells does not efficiently catalyze the repair of spontaneously damaged proteins.

Molecular phylogenetics of a protein repair methyltransferase

Protein-L-isoaspartyl (D-aspartyl) O-methyltransferase (E.C. 2.1.1.77) is a well-conserved and widely distributed protein repair enzyme that methylates isomerized or racemized aspartyl residues in age-damaged proteins. We exploited the availability of protein sequences from 10 diverse animal, plant and bacterial taxa to construct a phylogenetic tree and determine the rates of amino acid substitution for this enzyme. We used a likelihood ratio test to show that this enzyme fulfills the conditions for a molecular clock. We found that the rate of substitution is 0.39 amino acid substitutions per site per 10(9) years and remains relatively constant from bacteria to humans. We argue that this degree of sequence conservation may result from the functional constraints necessitated by the requirement to specifically recognize altered aspartyl but not normal aspartyl residues in proteins. Relative rate analysis of the Caenorhabditis elegans sequence suggests that the amino acid substitution rate in the nematode lineage may be higher than that in other lineages and that the divergence of nematodes may have been a more recent event than suggested by previous analysis.

Protein L-isoaspartyl methyltransferase from the nematode Caenorhabditis elegans: genomic structure and substrate specificity

We identified a protein L-isoaspartate (D-aspartate) O-methyltransferase (EC 2.1.1.77) in the nematode worm Caenorhabditis elegans. The methylation of abnormal L-isoaspartyl residues by this enzyme can lead to their conversion to L-aspartyl residues and represents a protein repair step for polypeptides damaged by spontaneous reactions during the aging process. We show that the levels of this enzyme increase 2-fold in C. elegans in the dauer larval form, a developmental stage where the organism can survive for extended periods of time. Utilizing degenerate oligonucleotide primers derived from conserved amino acid sequences of mammalian, plant, and bacterial L-isoaspartyl methyltransferases and PCR amplification, we made DNA probes that allowed us to obtain cDNA and genomic DNA clones encoding this enzyme in the nematode. The deduced amino acid sequence is 53% identical to the human enzyme and 29% identical to the Escherichia coli enzyme. Overexpression of the cDNA for the C. elegans enzyme in E. coli gave an active product with micromolar Km values for L-isoaspartyl-containing peptide substrates and for the methyl donor S-adenosyl-L-methionine. No methylation of D-aspartyl-containing peptides was detected under conditions where the human enzyme catalyzed this reaction, suggesting that the ability to methylate D-aspartyl residues in addition to L-isoaspartyl residues was a later evolutionary adaptation of this enzyme. The C. elegans gene for the methyltransferase, designated pcm-1, was mapped to a single site in a 31 kb region in the central portion of chromosome V. The gene is 3.2 kb in length and includes six introns. Although much smaller, its genomic organization is similar to that of the corresponding mouse gene, with identically positioned intron--exon splice junctions at five of seven sites. We propose that this gene plays an important role in facilitating the long term survival of this organisms.

Widespread occurrence of three sequence motifs in diverse S-adenosylmethionine-dependent methyltransferases suggests a common structure for these enzymes

Three regions of sequence similarity have been reported in several protein and small-molecule S-adenosylmethionine-dependent methyltransferases. Using multiple alignments, we have now identified these three regions in a much broader group of methyltransferases and have used these data to define a consensus for each region. Of the 84 non-DNA methyltransferase sequences in the GenBank, NBRF PIR, and Swissprot databases comprising 37 distinct enzymes, we have found 69 sequences possessing motif I. This motif is similar to a conserved region previously described in DNA adenine and cytosine methyltransferases. Motif II is found in 46 sequences, while motif III is found in 61 sequences. All three regions are found in 45 of these enzymes, and an additional 15 have motifs I and III. The motifs are always found in the same order on the polypeptide chain and are separated by comparable intervals. We suggest that these conserved regions contribute to the binding of the substrate S-adenosylmethionine and/or the product S-adenosylhomocysteine. These motifs can also be identified in certain nonmethyltransferases that utilize either S-adenosylmethionine or S-adenosylhomocysteine, including S-adenosylmethionine decarboxylase, S-adenosylmethionine synthetase, and S-adenosylhomocysteine hydrolase. In the latter two types of enzymes, motif I is similar to the conserved nucleotide binding motif of protein kinases and other nucleotide binding proteins. These motifs may be of use in predicting methyltransferases and related enzymes from the open reading frames generated by genomic sequencing projects.

Alternative splicing of the human isoaspartyl protein carboxyl methyltransferase RNA leads to the generation of a C-terminal -RDEL sequence in isozyme II

We have isolated two cDNA clones that correspond to the mRNAs for two isozymes of the human L-isoaspartyl/D-aspartyl protein carboxyl methyltransferase (EC 2.1.1.77). The DNA sequence of one of these encodes the amino acid sequence of the C-terminal half of the human erythrocyte isozyme I. The other cDNA clone includes the complete coding region of the more acidic isozyme II. With the exception of potential polymorphic sites at amino acid residues 119 and 205, the deduced amino acid sequences differ only at the C-terminus, where the -RWK sequence of isozyme I is replaced by a -RDEL sequence in isozyme II. The latter sequence is identical to a mammalian endoplasmic reticulum retention signal. With the previous evidence for only a single gene for the L-isoaspartyl/D-aspartyl methyltransferase in humans, and with evidence for consensus sites for alternative splicing in corresponding mouse genomic clones, we suggest that alternative splicing reactions can generate the major isozymes previously identified in human erythrocytes. The presence of alternative splicing leads us to predict the existence of a third isozyme with a -R C-terminus. The calculated isoelectric point of this third form is similar to that of a previously detected but uncharacterized minor methyltransferase activity.

Distinct C-terminal sequences of isozymes I and II of the human erythrocyte L-isoaspartyl/D-aspartyl protein methyltransferase

We have purified the more acidic major isozyme (II) of the human erythrocyte L-isoaspartyl/D-aspartyl methyltransferase and compared its structure to that of the previously sequenced isozyme I. These isozymes are both monomers of 25,000 molecular weight polypeptides and have similar enzymatic properties, but have isoelectric points that differ by one pH unit. Analysis of 16 tryptic peptides of isozyme II accounting for 89% of the sequence of isozyme I revealed no differences between these enzyme forms. However, analysis of a Staphylococcal V8 protease C-terminal fragment revealed that the last two residues of these proteins differed. The Trp-Lys-COOH terminus of isozyme I is replaced by a Asp-Asp-COOH terminus in isozyme II. Southern blot analysis of genomic DNA suggests that the human genome [corrected] may contain only a single gene encoding the enzyme. We propose that the distinct C-termini of isozymes I and II can arise from the generation of multiple mRNA's by alternative splicing.

Using redefinition and paradox with children in placement who provoke rejection

A model outlining an approach to helping placed children who provoke rejection analyzes the problem behavior from the perspective of the child and the family.