Epidemiology of carbapenem-resistant organisms in Alameda County, California, 2019-2021.

Objective: Carbapenem-resistant organisms (CROs) are an urgent health threat. Since 2017, Alameda County Health Public Health Department (ACPHD) mandates reporting of carbapenem-resistant Enterobacterales (CRE) and encourages voluntary reporting of non-CRE CROs including carbapenem-resistant Acinetobacter baumannii (CRAB) and carbapenem-resistant Pseudomonas aeruginosa (CRPA). Surveillance data from ACPHD were analyzed to describe the epidemiology of CROs and target public health interventions. Methods: Healthcare facilities in Alameda County reported CRO cases and submitted isolates to ACPHD to characterize carbapenemase genes; deaths were identified via the California Electronic Death Registration System. CRO cases with isolates resistant to one or more carbapenems were analyzed from surveillance data from July 2019 to June 2021. Results: Four hundred and forty-two cases of CROs were reported to Alameda County from 408 patients. The county case rate for CROs was 29 cases per 100,000 population, and cases significantly increased over the 2-year period. CRPA was most commonly reported (157 cases, 36%), and cases of CRAB increased 1.83-fold. One-hundred eighty-six (42%) cases were identified among residents of long-term care facilities; 152 (37%) patients had died by January 2022. One hundred and seven (24%) cases produced carbapenemases. Conclusions: The high burden of CROs in Alameda County highlights the need for continued partnership on reporting, testing, and infection prevention to limit the spread of resistant organisms. A large proportion of cases were identified in vulnerable long-term care residents, and CRAB was an emerging CRO among this population. Screening for CROs and surveillance at the local level are important to understand epidemiology and implement public health interventions.

TheiaEuk: a species-agnostic bioinformatics workflow for fungal genomic characterization.

Introduction: The clinical incidence of antimicrobial-resistant fungal infections has dramatically increased in recent years. Certain fungal pathogens colonize various body cavities, leading to life-threatening bloodstream infections. However, the identification and characterization of fungal isolates in laboratories remain a significant diagnostic challenge in medicine and public health. Whole-genome sequencing provides an unbiased and uniform identification pipeline for fungal pathogens but most bioinformatic analysis pipelines focus on prokaryotic species. To this end, TheiaEuk_Illumina_PE_PHB (TheiaEuk) was designed to focus on genomic analysis specialized to fungal pathogens. Methods: TheiaEuk was designed using containerized components and written in the workflow description language (WDL) to facilitate deployment on the cloud-based open bioinformatics platform Terra. This species-agnostic workflow enables the analysis of fungal genomes without requiring coding, thereby reducing the entry barrier for laboratory scientists. To demonstrate the usefulness of this pipeline, an ongoing outbreak of C. auris in southern Nevada was investigated. We performed whole-genome sequence analysis of 752 new C. auris isolates from this outbreak. Furthermore, TheiaEuk was utilized to observe the accumulation of mutations in the FKS1 gene over the course of the outbreak, highlighting the utility of TheiaEuk as a monitor of emerging public health threats when combined with whole-genome sequencing surveillance of fungal pathogens. Results: A primary result of this work is a curated fungal database containing 5,667 unique genomes representing 245 species. TheiaEuk also incorporates taxon-specific submodules for specific species, including clade-typing for Candida auris (C. auris). In addition, for several fungal species, it performs dynamic reference genome selection and variant calling, reporting mutations found in genes currently associated with antifungal resistance (FKS1, ERG11, FUR1). Using genome assemblies from the ATCC Mycology collection, the taxonomic identification module used by TheiaEuk correctly assigned genomes to the species level in 126/135 (93.3%) instances and to the genus level in 131/135 (97%) of instances, and provided zero false calls. Application of TheiaEuk to actual specimens obtained in the course of work at a local public health laboratory resulted in 13/15 (86.7%) correct calls at the species level, with 2/15 called at the genus level. It made zero incorrect calls. TheiaEuk accurately assessed clade type of Candida auris in 297/302 (98.3%) of instances. Discussion: TheiaEuk demonstrated effectiveness in identifying fungal species from whole genome sequence. It further showed accuracy in both clade-typing of C. auris and in the identification of mutations known to associate with drug resistance in that organism.

The role of plasmids in carbapenem resistant E. coli in Alameda County, California.

BACKGROUND: Antimicrobial resistant infections continue to be a leading global public health crisis. Mobile genetic elements, such as plasmids, have been shown to play a major role in the dissemination of antimicrobial resistance (AMR) genes. Despite its ongoing threat to human health, surveillance of AMR in the United States is often limited to phenotypic resistance. Genomic analyses are important to better understand the underlying resistance mechanisms, assess risk, and implement appropriate prevention strategies. This study aimed to investigate the extent of plasmid mediated antimicrobial resistance that can be inferred from short read sequences of carbapenem resistant E. coli (CR-Ec) in Alameda County, California. E. coli isolates from healthcare locations in Alameda County were sequenced using an Illumina MiSeq and assembled with Unicycler. Genomes were categorized according to predefined multilocus sequence typing (MLST) and core genome multilocus sequence typing (cgMLST) schemes. Resistance genes were identified and corresponding contigs were predicted to be plasmid-borne or chromosome-borne using two bioinformatic tools (MOB-suite and mlplasmids). RESULTS: Among 82 of CR-Ec identified between 2017 and 2019, twenty-five sequence types (STs) were detected. ST131 was the most prominent (n = 17) followed closely by ST405 (n = 12). blaCTX-M were the most common ESBL genes and just over half (18/30) of these genes were predicted to be plasmid-borne by both MOB-suite and mlplasmids. Three genetically related groups of E. coli isolates were identified with cgMLST. One of the groups contained an isolate with a chromosome-borne blaCTX-M-15 gene and an isolate with a plasmid-borne blaCTX-M-15 gene. CONCLUSIONS: This study provides insights into the dominant clonal groups driving carbapenem resistant E. coli infections in Alameda County, CA, USA clinical sites and highlights the relevance of whole-genome sequencing in routine local genomic surveillance. The finding of multi-drug resistant plasmids harboring high-risk resistance genes is of concern as it indicates a risk of dissemination to previously susceptible clonal groups, potentially complicating clinical and public health intervention.

GAMBIT (Genomic Approximation Method for Bacterial Identification and Tracking): A methodology to rapidly leverage whole genome sequencing of bacterial isolates for clinical identification.

Whole genome sequencing (WGS) of clinical bacterial isolates has the potential to transform the fields of diagnostics and public health. To realize this potential, bioinformatic software that reports identification results needs to be developed that meets the quality standards of a diagnostic test. We developed GAMBIT (Genomic Approximation Method for Bacterial Identification and Tracking) using k-mer based strategies for identification of bacteria based on WGS reads. GAMBIT incorporates this algorithm with a highly curated searchable database of 48,224 genomes. Herein, we describe validation of the scoring methodology, parameter robustness, establishment of confidence thresholds and the curation of the reference database. We assessed GAMBIT by way of validation studies when it was deployed as a laboratory-developed test in two public health laboratories. This method greatly reduces or eliminates false identifications which are often detrimental in a clinical setting.

Clonal Lineages and Virulence Factors of Carbapenem Resistant E. coli in Alameda County, California, 2017-2019.

The prevalence of carbapenem-resistant Enterobacterales (CRE) has been increasing since the year 2000 and is considered a serious public health threat according to the Centers for Disease Control and Prevention. Limited studies have genotyped Carbapenem-resistant Escherichia coli using whole genome sequencing to characterize the most common lineages and resistance and virulence genes. The aim of this study was to characterize sequence data from carbapenem-resistant E. coli isolates (n = 82) collected longitudinally by the Alameda County Public Health Laboratory (ACPHL) between 2017 and 2019. E. coli genomes were screened for antibiotic resistance genes (ARGs) and extraintestinal pathogenic E. coli virulence factor genes (VFGs). The carbapenem-resistant E. coli lineages were diverse, with 24 distinct sequence types (STs) represented, including clinically important STs: ST131, ST69, ST95, and ST73. All Ambler classes of Carbapenemases were present, with NDM-5 being most the frequently detected. Nearly all isolates (90%) contained genes encoding resistance to third-generation cephalosporins; blaCTX-M genes were most common. The number of virulence genes present within pandemic STs was significantly higher than the number in non-pandemic lineages (p = 0.035). Virulence genes fimA (92%), trat (71%), kpsM (54%), and iutA (46%) were the most prevalent within the isolates. Considering the public health risk associated with CRE, these data enhance our understanding of the diversity of clinically important E. coli that are circulating in Alameda County, California.

Imipenem-Relebactam Susceptibility and Genotypic Characteristics of Carbapenem-Resistant Enterobacterales Identified during Population-Based Surveillance.

Laboratories submit all carbapenem-resistant Enterobacter, Escherichia coli, and Klebsiella species to the Alameda County Public Health Department (ACPHD). ACPHD evaluated 75 isolates submitted during 9 months for susceptibility to imipenem-relebactam (I-R) and, using whole-genome sequencing, identified ß-lactamase genes. Of 60 (80%) isolates susceptible to I-R, 8 (13%) had detectable carbapenemase genes, including 4 KPC, two NDM, and two OXA-48-like; we described the relationship between the presence of ß-lactamase resistance genes and susceptibility to I-R.

Comparison of an In Vitro Diagnostic Next-Generation Sequencing Assay with Sanger Sequencing for HIV-1 Genotypic Resistance Testing.

The ability of next-generation sequencing (NGS) technologies to detect low frequency HIV-1 drug resistance mutations (DRMs) not detected by dideoxynucleotide Sanger sequencing has potential advantages for improved patient outcomes. We compared the performance of an in vitro diagnostic (IVD) NGS assay, the Sentosa SQ HIV genotyping assay for HIV-1 genotypic resistance testing, with Sanger sequencing on 138 protease/reverse transcriptase (RT) and 39 integrase sequences. The NGS assay used a 5% threshold for reporting low-frequency variants. The level of complete plus partial nucleotide sequence concordance between Sanger sequencing and NGS was 99.9%. Among the 138 protease/RT sequences, a mean of 6.4 DRMs was identified by both Sanger and NGS, a mean of 0.5 DRM was detected by NGS alone, and a mean of 0.1 DRM was detected by Sanger sequencing alone. Among the 39 integrase sequences, a mean of 1.6 DRMs was detected by both Sanger sequencing and NGS and a mean of 0.15 DRM was detected by NGS alone. Compared with Sanger sequencing, NGS estimated higher levels of resistance to one or more antiretroviral drugs for 18.2% of protease/RT sequences and 5.1% of integrase sequences. There was little evidence for technical artifacts in the NGS sequences, but the G-to-A hypermutation was detected in three samples. In conclusion, the IVD NGS assay evaluated in this study was highly concordant with Sanger sequencing. At the 5% threshold for reporting minority variants, NGS appeared to attain a modestly increased sensitivity for detecting low-frequency DRMs without compromising sequence accuracy.

Diff-seq: A high throughput sequencing-based mismatch detection assay for DNA variant enrichment and discovery.

Much of the within species genetic variation is in the form of single nucleotide polymorphisms (SNPs), typically detected by whole genome sequencing (WGS) or microarray-based technologies. However, WGS produces mostly uninformative reads that perfectly match the reference, while microarrays require genome-specific reagents. We have developed Diff-seq, a sequencing-based mismatch detection assay for SNP discovery without the requirement for specialized nucleic-acid reagents. Diff-seq leverages the Surveyor endonuclease to cleave mismatched DNA molecules that are generated after cross-annealing of a complex pool of DNA fragments. Sequencing libraries enriched for Surveyor-cleaved molecules result in increased coverage at the variant sites. Diff-seq detected all mismatches present in an initial test substrate, with specific enrichment dependent on the identity and context of the variation. Application to viral sequences resulted in increased observation of variant alleles in a biologically relevant context. Diff-Seq has the potential to increase the sensitivity and efficiency of high-throughput sequencing in the detection of variation.

Evolution of gag and gp41 in Patients Receiving Ritonavir-Boosted Protease Inhibitors.

Several groups have proposed that genotypic determinants in gag and the gp41 cytoplasmic domain (gp41-CD) reduce protease inhibitor (PI) susceptibility without PI-resistance mutations in protease. However, no gag and gp41-CD mutations definitively responsible for reduced PI susceptibility have been identified in individuals with virological failure (VF) while receiving a boosted PI (PI/r)-containing regimen. To identify gag and gp41 mutations under selective PI pressure, we sequenced gag and/or gp41 in 61 individuals with VF on a PI/r (n = 40) or NNRTI (n = 20) containing regimen. We quantified nonsynonymous and synonymous changes in both genes and identified sites exhibiting signal for directional or diversifying selection. We also used published gag and gp41 polymorphism data to highlight mutations displaying a high selection index, defined as changing from a conserved to an uncommon amino acid. Many amino acid mutations developed in gag and in gp41-CD in both the PI- and NNRTI-treated groups. However, in neither gene, were there discernable differences between the two groups in overall numbers of mutations, mutations displaying evidence of diversifying or directional selection, or mutations with a high selection index. If gag and/or gp41 encode PI-resistance mutations, they may not be confined to consistent mutations at a few sites.

Prevalence of Drug-Resistant Minority Variants in Untreated HIV-1-Infected Individuals With and Those Without Transmitted Drug Resistance Detected by Sanger Sequencing.

Minority variant human immunodeficiency virus type 1 (HIV-1) nonnucleoside reverse transcriptase inhibitor (NNRTI) resistance mutations are associated with an increased risk of virological failure during treatment with NNRTI-containing regimens. To determine whether individuals to whom variants with isolated NNRTI-associated drug resistance were transmitted are at increased risk of virological failure during treatment with a non-NNRTI-containing regimen, we identified minority variant resistance mutations in 33 individuals with isolated NNRTI-associated transmitted drug resistance and 49 matched controls. We found similar proportions of overall and nucleoside reverse transcriptase inhibitor-associated minority variant resistance mutations in both groups, suggesting that isolated NNRTI-associated transmitted drug resistance may not be a risk factor for virological failure during treatment with a non-NNRTI-containing regimen.

Mutational Correlates of Virological Failure in Individuals Receiving a WHO-Recommended Tenofovir-Containing First-Line Regimen: An International Collaboration.

Tenofovir disoproxil fumarate (TDF) genotypic resistance defined by K65R/N and/or K70E/Q/G occurs in 20% to 60% of individuals with virological failure (VF) on a WHO-recommended TDF-containing first-line regimen. However, the full spectrum of reverse transcriptase (RT) mutations selected in individuals with VF on such a regimen is not known. To identify TDF regimen-associated mutations (TRAMs), we compared the proportion of each RT mutation in 2873 individuals with VF on a WHO-recommended first-line TDF-containing regimen to its proportion in a cohort of 50,803 antiretroviral-naïve individuals. To identify TRAMs specifically associated with TDF-selection pressure, we compared the proportion of each TRAM to its proportion in a cohort of 5805 individuals with VF on a first-line thymidine analog-containing regimen. We identified 83 TRAMs including 33 NRTI-associated, 40 NNRTI-associated, and 10 uncommon mutations of uncertain provenance. Of the 33 NRTI-associated TRAMs, 12 - A62V, K65R/N, S68G/N/D, K70E/Q/T, L74I, V75L, and Y115F - were more common among individuals receiving a first-line TDF-containing compared to a first-line thymidine analog-containing regimen. These 12 TDF-selected TRAMs will be important for monitoring TDF-associated transmitted drug-resistance and for determining the extent of reduced TDF susceptibility in individuals with VF on a TDF-containing regimen.

Evaluation of the Aptima HIV-1 Quant Dx Assay Using Plasma and Dried Blood Spots.

HIV-1 RNA quantitation in plasma, or virus load testing, is the primary method by which the response to antiretroviral therapy is monitored. Here we describe evaluation of the Aptima HIV-1 Quant Dx assay (Aptima) performed on the automated Panther system. The clinical performance of Aptima was compared to that of the Cobas AmpliPrep/Cobas TaqMan HIV-1 Test v2.0 (CAP/CTM) using 162 EDTA plasma samples collected from patients undergoing HIV-1 monitoring. Overall agreement was 84.0% (136/162), with a kappa statistic of 0.723 (standard error, 0.047; 95% confidence interval [CI], 0.630 to 0.815), indicating substantial agreement. Using the 86 clinical samples quantifiable by both methods, Passing-Bablok regression revealed a regression line of Y = (1.069 × X) - 0.346 (95% CI of the slope [1.003 to 1.139] and intercept [-0.666 to -0.074]), and Bland-Altman analysis demonstrated a mean difference (Aptima-CAP/CTM) of -0.075 log10 copies/ml (95% limits of agreement of -0.624 to 0.475), consistent with negative bias. Comparison of Aptima results for paired dried blood spot (DBS) and plasma specimens archived from participants in the Peninsula AIDS Research Cohort Study (PARC) demonstrated an overall agreement of 94.7% (90/95) when 1,000 copies/ml was used as the threshold. In conclusion, the Aptima HIV-1 Quant Dx assay provides a suitable alternative for HIV-1 monitoring in plasma and DBS.

HIV-1 Protease, Reverse Transcriptase, and Integrase Variation.

UNLABELLED: HIV-1 protease (PR), reverse transcriptase (RT), and integrase (IN) variability presents a challenge to laboratories performing genotypic resistance testing. This challenge will grow with increased sequencing of samples enriched for proviral DNA such as dried blood spots and increased use of next-generation sequencing (NGS) to detect low-abundance HIV-1 variants. We analyzed PR and RT sequences from >100,000 individuals and IN sequences from >10,000 individuals to characterize variation at each amino acid position, identify mutations indicating APOBEC-mediated G-to-A editing, and identify mutations resulting from selective drug pressure. Forty-seven percent of PR, 37% of RT, and 34% of IN positions had one or more amino acid variants with a prevalence of ≥1%. Seventy percent of PR, 60% of RT, and 60% of IN positions had one or more variants with a prevalence of ≥0.1%. Overall 201 PR, 636 RT, and 346 IN variants had a prevalence of ≥0.1%. The median intersubtype prevalence ratios were 2.9-, 2.1-, and 1.9-fold for these PR, RT, and IN variants, respectively. Only 5.0% of PR, 3.7% of RT, and 2.0% of IN variants had a median intersubtype prevalence ratio of ≥10-fold. Variants at lower prevalences were more likely to differ biochemically and to be part of an electrophoretic mixture compared to high-prevalence variants. There were 209 mutations indicative of APOBEC-mediated G-to-A editing and 326 mutations nonpolymorphic treatment selected. Identification of viruses with a high number of APOBEC-associated mutations will facilitate the quality control of dried blood spot sequencing. Identifying sequences with a high proportion of rare mutations will facilitate the quality control of NGS. IMPORTANCE: Most antiretroviral drugs target three HIV-1 proteins: PR, RT, and IN. These proteins are highly variable: many different amino acids can be present at the same position in viruses from different individuals. Some of the amino acid variants cause drug resistance and occur mainly in individuals receiving antiretroviral drugs. Some variants result from a human cellular defense mechanism called APOBEC-mediated hypermutation. Many variants result from naturally occurring mutation. Some variants may represent technical artifacts. We studied PR and RT sequences from >100,000 individuals and IN sequences from >10,000 individuals to quantify variation at each amino acid position in these three HIV-1 proteins. We performed analyses to determine which amino acid variants resulted from antiretroviral drug selection pressure, APOBEC-mediated editing, and naturally occurring variation. Our results provide information essential to clinical, research, and public health laboratories performing genotypic resistance testing by sequencing HIV-1 PR, RT, and IN.

Q148N, a Novel Integrase Inhibitor Resistance Mutation Associated with Low-Level Reduction in Elvitegravir Susceptibility.

The integrase strand transfer inhibitor (INSTI)-resistance mutations Q148H/K/R are arguably the most important INSTI-resistance mutations as they represented the first step to high-level dolutegravir cross-resistance. We describe an individual with transmitted four-class drug resistance whose virus sequence had the previously uncharacterized mutation Q148N. Infectious molecular HIV-1 clones containing Q148N alone and in combination with G140S demonstrated ∼2.4-4.5 reduced elvitegravir susceptibility depending on the virus's genetic context but retained susceptibility to raltegravir and dolutegravir. This level of reduced elvitegravir susceptibility is lower than that observed with Q148H/K/R and in fact the infected individual responded to an initial treatment regimen containing tenofovir/emtricitabine/elvitegravir/cobicistat. Q148N was associated with a higher replication capacity than Q148H, suggesting that this mutation may be more fit in the absence of selective INSTI therapy.

Geographic and temporal trends in the molecular epidemiology and genetic mechanisms of transmitted HIV-1 drug resistance: an individual-patient- and sequence-level meta-analysis.

BACKGROUND: Regional and subtype-specific mutational patterns of HIV-1 transmitted drug resistance (TDR) are essential for informing first-line antiretroviral (ARV) therapy guidelines and designing diagnostic assays for use in regions where standard genotypic resistance testing is not affordable. We sought to understand the molecular epidemiology of TDR and to identify the HIV-1 drug-resistance mutations responsible for TDR in different regions and virus subtypes. METHODS AND FINDINGS: We reviewed all GenBank submissions of HIV-1 reverse transcriptase sequences with or without protease and identified 287 studies published between March 1, 2000, and December 31, 2013, with more than 25 recently or chronically infected ARV-naïve individuals. These studies comprised 50,870 individuals from 111 countries. Each set of study sequences was analyzed for phylogenetic clustering and the presence of 93 surveillance drug-resistance mutations (SDRMs). The median overall TDR prevalence in sub-Saharan Africa (SSA), south/southeast Asia (SSEA), upper-income Asian countries, Latin America/Caribbean, Europe, and North America was 2.8%, 2.9%, 5.6%, 7.6%, 9.4%, and 11.5%, respectively. In SSA, there was a yearly 1.09-fold (95% CI: 1.05-1.14) increase in odds of TDR since national ARV scale-up attributable to an increase in non-nucleoside reverse transcriptase inhibitor (NNRTI) resistance. The odds of NNRTI-associated TDR also increased in Latin America/Caribbean (odds ratio [OR] = 1.16; 95% CI: 1.06-1.25), North America (OR = 1.19; 95% CI: 1.12-1.26), Europe (OR = 1.07; 95% CI: 1.01-1.13), and upper-income Asian countries (OR = 1.33; 95% CI: 1.12-1.55). In SSEA, there was no significant change in the odds of TDR since national ARV scale-up (OR = 0.97; 95% CI: 0.92-1.02). An analysis limited to sequences with mixtures at less than 0.5% of their nucleotide positions­a proxy for recent infection­yielded trends comparable to those obtained using the complete dataset. Four NNRTI SDRMs­K101E, K103N, Y181C, and G190A­accounted for >80% of NNRTI-associated TDR in all regions and subtypes. Sixteen nucleoside reverse transcriptase inhibitor (NRTI) SDRMs accounted for >69% of NRTI-associated TDR in all regions and subtypes. In SSA and SSEA, 89% of NNRTI SDRMs were associated with high-level resistance to nevirapine or efavirenz, whereas only 27% of NRTI SDRMs were associated with high-level resistance to zidovudine, lamivudine, tenofovir, or abacavir. Of 763 viruses with TDR in SSA and SSEA, 725 (95%) were genetically dissimilar; 38 (5%) formed 19 sequence pairs. Inherent limitations of this study are that some cohorts may not represent the broader regional population and that studies were heterogeneous with respect to duration of infection prior to sampling. CONCLUSIONS: Most TDR strains in SSA and SSEA arose independently, suggesting that ARV regimens with a high genetic barrier to resistance combined with improved patient adherence may mitigate TDR increases by reducing the generation of new ARV-resistant strains. A small number of NNRTI-resistance mutations were responsible for most cases of high-level resistance, suggesting that inexpensive point-mutation assays to detect these mutations may be useful for pre-therapy screening in regions with high levels of TDR. In the context of a public health approach to ARV therapy, a reliable point-of-care genotypic resistance test could identify which patients should receive standard first-line therapy and which should receive a protease-inhibitor-containing regimen.

A uniquely prevalent nonnucleoside reverse transcriptase inhibitor resistance mutation in Russian subtype A HIV-1 viruses.

BACKGROUND: The subtype A variant in the Former Soviet Union (A(FSU)) causes most of Russia's HIV-1 infections. However, the spectrum of drug-resistance mutations (DRMs) in antiretroviral experienced patients with this variant has not been studied. METHODS: Between 2010 and 2013, genotypic resistance testing was performed on plasma samples from 366 antiretroviral-experienced patients in Siberia. RESULTS: Three-hundred patients (82%) had subtype A(FSU) and 55 (15%) had CRF02_AG viruses. The pattern of DRMs was consistent with patient antiretroviral history with one exception. G190S was the most common nonnucleoside reverse transcriptase inhibitor (NNRTI) resistance mutation, occurring in 55 (33%) subtype A(FSU) viruses from 167 NNRTI-experienced patients compared with none of 37 CRF02_AG viruses from NNRTI-experienced patients (P < 0.001). The next most common subtype A(FSU) NNRTI-resistance mutation, K103N, occurred in 25 (15%) viruses. Wild-type glycine (G) at position 190 is encoded by GGC in more than 99% of published A(FSU) strains. By contrast, G190 is encoded by GGA or GGG in 97% of other subtypes and in subtype A strains outside of the FSU. Therefore, G190S results from a single G→A transition: G (GGC) → S (AGC) almost exclusively in subtype A(FSU) viruses. CONCLUSION: The predisposition of subtype A(FSU) to G190S is concerning because G→A is the most common HIV-1 mutation and because G190S causes higher levels of nevirapine and efavirenz resistance than K103N. This study exemplifies the need for characterizing the genetic mechanisms of resistance in diverse populations and warrants studies to verify that NRTI/NNRTI regimens are as efficacious in treating subtype A(FSU) as viruses belonging to other subtypes.

Non-nucleoside reverse transcriptase inhibitor (NNRTI) cross-resistance: implications for preclinical evaluation of novel NNRTIs and clinical genotypic resistance testing.

OBJECTIVES: The introduction of two new non-nucleoside reverse transcriptase inhibitors (NNRTIs) in the past 5 years and the identification of novel NNRTI-associated mutations have made it necessary to reassess the extent of phenotypic NNRTI cross-resistance. METHODS: We analysed a dataset containing 1975, 1967, 519 and 187 genotype-phenotype correlations for nevirapine, efavirenz, etravirine and rilpivirine, respectively. We used linear regression to estimate the effects of RT mutations on susceptibility to each of these NNRTIs. RESULTS: Sixteen mutations at 10 positions were significantly associated with the greatest contribution to reduced phenotypic susceptibility (≥10-fold) to one or more NNRTIs, including: 14 mutations at six positions for nevirapine (K101P, K103N/S, V106A/M, Y181C/I/V, Y188C/L and G190A/E/Q/S); 10 mutations at six positions for efavirenz (L100I, K101P, K103N, V106M, Y188C/L and G190A/E/Q/S); 5 mutations at four positions for etravirine (K101P, Y181I/V, G190E and F227C); and 6 mutations at five positions for rilpivirine (L100I, K101P, Y181I/V, G190E and F227C). G190E, a mutation that causes high-level nevirapine and efavirenz resistance, also markedly reduced susceptibility to etravirine and rilpivirine. K101H, E138G, V179F and M230L mutations, associated with reduced susceptibility to etravirine and rilpivirine, were also associated with reduced susceptibility to nevirapine and/or efavirenz. CONCLUSIONS: The identification of novel cross-resistance patterns among approved NNRTIs illustrates the need for a systematic approach for testing novel NNRTIs against clinical virus isolates with major NNRTI-resistance mutations and for testing older NNRTIs against virus isolates with mutations identified during the evaluation of a novel NNRTI.

Prototypical Recombinant Multi-Protease-Inhibitor-Resistant Infectious Molecular Clones of Human Immunodeficiency Virus Type 1.

The many genetic manifestations of HIV-1 protease inhibitor (PI) resistance present challenges to research into the mechanisms of PI resistance and the assessment of new PIs. To address these challenges, we created a panel of recombinant multi-PI-resistant infectious molecular clones designed to represent the spectrum of clinically relevant multi-PI-resistant viruses. To assess the representativeness of this panel, we examined the sequences of the panel's viruses in the context of a correlation network of PI resistance amino acid substitutions in sequences from more than 10,000 patients. The panel of recombinant infectious molecular clones comprised 29 of 41 study-defined PI resistance amino acid substitutions and 23 of the 27 tightest amino acid substitution clusters. Based on their phenotypic properties, the clones were classified into four groups with increasing cross-resistance to the PIs most commonly used for salvage therapy: lopinavir (LPV), tipranavir (TPV), and darunavir (DRV). The panel of recombinant infectious molecular clones has been made available without restriction through the NIH AIDS Research and Reference Reagent Program. The public availability of the panel makes it possible to compare the inhibitory activities of different PIs with one another. The diversity of the panel and the high-level PI resistance of its clones suggest that investigational PIs active against the clones in this panel will retain antiviral activity against most if not all clinically relevant PI-resistant viruses.

Collinearity of protease mutations in HIV-1 samples with high-level protease inhibitor class resistance.

OBJECTIVES: To determine whether pan-protease inhibitor (PI)-resistant virus populations are composed predominantly of viruses with resistance to all PIs or of diverse virus populations with resistance to different subsets of PIs. METHODS: We performed deep sequencing of plasma virus samples from nine patients with high-level genotypic and/or phenotypic resistance to all licensed PIs. The nine virus samples had a median of 12 PI resistance mutations by direct PCR Sanger sequencing. RESULTS: For each of the nine virus samples, deep sequencing showed that each of the individual viruses within a sample contained nearly all of the mutations detected by Sanger sequencing. Indeed, a median of 94.9% of deep sequence reads had each of the PI resistance mutations present as a single chromatographic peak in the Sanger sequence. A median of 5.0% of reads had all but one of the Sanger mutations that were not part of an electrophoretic mixture. CONCLUSIONS: The collinearity of PI resistance mutations in the nine virus samples demonstrated that pan-PI-resistant viruses are able to replicate in vivo despite their highly mutated protease enzymes. We hypothesize that the marked collinearity of PI resistance mutations in pan-PI-resistant virus populations results from the unique requirements for multi-PI resistance and the extensive cross-resistance conferred by many of the accessory PI resistance mutations.