Dataset Information

Analysis of unusual and signature APOBEC-mutations in HIV-1 pol next-generation sequences.

ABSTRACT:

Introduction

At low mutation-detection thresholds, next generation sequencing (NGS) for HIV-1 genotypic resistance testing is susceptible to artifactual detection of mutations arising from PCR error and APOBEC-mediated G-to-A hypermutation.

Methods

We analyzed published HIV-1 pol Illumina NGS data to characterize the distribution of mutations at eight NGS mutation detection thresholds: 20%, 10%, 5%, 2%, 1%, 0.5%, 0.2%, and 0.1%. At each threshold, we determined proportions of amino acid mutations that were unusual (defined as having a prevalence <0.01% in HIV-1 group M sequences) or signature APOBEC mutations.

Results

Eight studies, containing 855 samples, in the NCBI Sequence Read Archive were analyzed. As detection thresholds were lowered, there was a progressive increase in the proportion of positions with usual and unusual mutations and in the proportion of all mutations that were unusual. The median proportion of positions with an unusual mutation increased gradually from 0% at the 20% threshold to 0.3% at the 1% threshold and then exponentially to 1.3% (0.5% threshold), 6.9% (0.2% threshold), and 23.2% (0.1% threshold). In two of three studies with available plasma HIV-1 RNA levels, the proportion of positions with unusual mutations was negatively associated with virus levels. Although the complete set of signature APOBEC mutations was much smaller than that of unusual mutations, the former outnumbered the latter in one-sixth of samples at the 0.5%, 1%, and 2% thresholds.

Conclusions

The marked increase in the proportion of positions with unusual mutations at thresholds below 1% and in samples with lower virus loads suggests that, at low thresholds, many unusual mutations are artifactual, reflecting PCR error or G-to-A hypermutation. Profiling the numbers of unusual and signature APOBEC pol mutations at different NGS mutation detection thresholds may be useful to avoid selecting a threshold that is too low and poses an unacceptable risk of identifying artifactual mutations.

SUBMITTER: Tzou PL

PROVIDER: S-EPMC7043932 | biostudies-literature |

REPOSITORIES: biostudies-literature

ACCESS DATA

Similar Datasets

Project description:Patients with antiretroviral therapy interruption have a high risk of virological failure when re-initiating antiretroviral therapy (ART), especially those with HIV drug resistance. Next-generation sequencing may provide close scrutiny on their minority drug resistance variant. A cross-sectional study was conducted in patients with ART interruption in five regions in China in 2016. Through Sanger and next-generation sequencing in parallel, HIV drug resistance was genotyped on their plasma samples. Rates of HIV drug resistance were compared by the McNemar tests. In total, 174 patients were included in this study, with a median 12 (interquartile range (IQR), 6-24) months of ART interruption. Most (86.2%) of them had received efavirenz (EFV)/nevirapine (NVP)-based first-line therapy for a median 16 (IQR, 7-26) months before ART interruption. Sixty-one (35.1%) patients had CRF07_BC HIV-1 strains, 58 (33.3%) CRF08_BC and 35 (20.1%) CRF01_AE. Thirty-four (19.5%) of the 174 patients were detected to harbor HIV drug-resistant variants on Sanger sequencing. Thirty-six (20.7%), 37 (21.3%), 42 (24.1%), 79 (45.4%) and 139 (79.9) patients were identified to have HIV drug resistance by next-generation sequencing at 20% (v.s. Sanger, p = 0.317), 10% (v.s. Sanger, p = 0.180), 5% (v.s. Sanger, p = 0.011), 2% (v.s. Sanger, p < 0.001) and 1% (v.s. Sanger, p < 0.001) of detection thresholds, respectively. K65R was the most common minority mutation, of 95.1% (58/61) and 93.1% (54/58) in CRF07_BC and CRF08_BC, respectively, when compared with 5.7% (2/35) in CRF01_AE (p < 0.001). In 49 patients that followed-up a median 10 months later, HIV drug resistance mutations at >20% frequency such as K103N, M184VI and P225H still existed, but with decreased frequencies. The prevalence of HIV drug resistance in ART interruption was higher than 15% in the survey. Next-generation sequencing was able to detect more minority drug resistance variants than Sanger. There was a sharp increase in minority drug resistance variants when the detection threshold was below 5%.

Project description:Merkel cell carcinoma is a rare but highly aggressive cutaneous neuroendocrine carcinoma. Cytokeratin 20 (CK20) is expressed in ~95% of Merkel cell carcinomas and is useful for distinction from morphologically similar entities including metastatic small-cell lung carcinoma. Lack of CK20 expression may make diagnosis of Merkel cell carcinoma more challenging, and has unknown biological significance. Approximately 80% of CK20-positive Merkel cell carcinomas are associated with the oncogenic Merkel cell polyomavirus. Merkel cell carcinomas lacking Merkel cell polyomavirus display distinct genetic changes from Merkel cell polyomavirus-positive Merkel cell carcinoma, including RB1 inactivating mutations. Unlike CK20-positive Merkel cell carcinoma, the majority of CK20-negative Merkel cell carcinomas are Merkel cell polyomavirus-negative, suggesting CK20-negative Merkel cell carcinomas predominantly arise through virus-independent pathway(s) and may harbor additional genetic differences from conventional Merkel cell carcinoma. Hence, we analyzed 15 CK20-negative Merkel cell carcinoma tumors (10 Merkel cell polyomavirus-negative, four Merkel cell polyomavirus-positive, and one undetermined) using the Ion Ampliseq Comprehensive Cancer Panel, which assesses copy number alterations and mutations in 409 cancer-relevant genes. Twelve tumors displayed prioritized high-level chromosomal gains or losses (average 1.9 per tumor). Non-synonymous high-confidence somatic mutations were detected in 14 tumors (average 11.9 per tumor). Assessing all somatic coding mutations, an ultraviolet-signature mutational profile was present, and more prevalent in Merkel cell polyomavirus-negative tumors. Recurrent deleterious tumor suppressor mutations affected TP53 (9/15, 60%), RB1 (3/15, 20%), and BAP1 (2/15, 13%). Oncogenic activating mutations included PIK3CA (3/15, 20%), AKT1 (1/15, 7%) and EZH2 (1/15, 7%). In conclusion, CK20-negative Merkel cell carcinoma display overlapping genetic changes with CK20-positive Merkel cell carcinoma, including RB1 mutations restricted to Merkel cell polyomavirus-negative tumors. However, some CK20-negative Merkel cell carcinomas harbor mutations not previously described in Merkel cell carcinoma. Hence, CK20-negative Merkel cell carcinomas harbor diverse oncogenic drivers which may represent therapeutic targets in individual tumors.