Indel and Carryforward Correction (ICC): a new analysis approach for processing 454 pyrosequencing data

Empirical evidence for concerted evolution in the 18S rDNA region of the planktonic diatom genus Chaetoceros

Concerted evolution is a process of homogenisation of repetitive sequences within a genome through unequal crossing over and gene conversion. This homogenisation is never fully achieved because mutations always create new variants. Classically, concerted evolution has been detected as "noise" in electropherograms and these variants have been characterised through cloning and sequencing of subsamples of amplified products. However, this approach limits the number of detectable variants and provides no information about the abundance of each variant. In this study, we investigated concerted evolution by using environmental time-series metabarcoding data, single strain high-throughput sequencing (HTS) and a collection of Sanger reference barcode sequences. We used six species of the marine planktonic diatom genus Chaetoceros as study system. Abundance plots obtained from environmental metabarcoding and single strain HTS showed the presence of a haplotype far more abundant than all the others (the "dominant" haplotype) and identical to the reference sequences of that species obtained with Sanger sequencing. This distribution fitted best with Zipf's law among the rank abundance/ dominance models tested. Furthermore, in each strain 99% of reads showed a similarity of 99% with the dominant haplotype, confirming the efficiency of the homogenisation mechanism of concerted evolution. We also demonstrated that minor haplotypes found in the environmental samples are not only technical artefacts, but mostly intragenomic variation generated by incomplete homogenisation. Finally, we showed that concerted evolution can be visualised inferring phylogenetic networks from environmental data. In conclusion, our study provides an important contribution to the understanding of concerted evolution and to the interpretation of DNA barcoding and metabarcoding data based on multigene family markers.

An Inexpensive and Accurate Reverse Transcription-PCR-Melting Temperature Analysis Assay for Real-Time Influenza Virus B Lineage Discrimination

In this work, we describe a SYBR-Green one-step reverse transcription-PCR protocol coupled with a melting temperature analysis (RT-PCR-T_m ), which allows the discrimination of influenza B lineages Yamagata and Victoria. The assay is performed using a regular real-time thermocycler and is based on differences in melting temperature (T_m ) of a 131-bp amplicon, obtained from a conserved region of hemagglutinin gene. A total of 410 samples collected during the 2004, 2008, and 2010-2017 influenza seasons in Brazil were tested, and the lineages were correctly characterized using their melting profiles. The temperature range is significantly different between both lineages throughout the time (Mann-Whitney test; P < 0.0001, confidence interval = 95%), and the T_m is not affected by viral load (Spearman correlation test; r = 0.287, P = 2.245 × 10^-9). The simplicity and cost-effectiveness of this protocol make it an option for influenza B lineage surveillance worldwide.

High-throughput sequencing of murine immunoglobulin heavy chain repertoires using single side unique molecular identifiers on an Ion Torrent PGM

With the advent of high-throughput sequencing (HTS), profiling immunoglobulin (IG) repertoires has become an essential part of immunological research. Advances in sequencing technology enable the IonTorrent Personal Genome Machine (PGM) to cover the full-length of IG mRNA transcripts. Nucleotide insertions and deletions (indels) are the dominant errors of the PGM sequencing platform and can critically influence IG repertoire assessments. Here, we present a PGM-tailored IG repertoire sequencing approach combining error correction through unique molecular identifier (UID) barcoding and indel detection through ImMunoGeneTics (IMGT), the most commonly used sequence alignment database for IG sequences. Using artificially falsified sequences for benchmarking, we found that IMGT's underlying algorithms efficiently detect 98% of the introduced indels. Undetected indels are either located at the end of the sequences or produce masked frameshifts with an insertion and deletion in close proximity. The complementary determining regions 3 (CDR3s) are returned correct for up to 3 insertions or 3 deletions through conservative culling. We further show, that our PGM-tailored unique molecular identifiers result in highly accurate HTS data if combined with the presented processing strategy. In this regard, considering sequences with at least two copies from datasets with UID families of minimum 3 reads result in correct sequences with over 99% confidence. Finally, we show that the protocol can readily be used to generate homogenous datasets for bulk sequencing of murine bone marrow samples. Taken together, this approach will help to establish benchtop-scale sequencing of IG heavy chain transcripts in the field of IG repertoire research.

Transmission of HIV-1 drug resistance mutations within partner-pairs: A cross-sectional study of a primary HIV infection cohort

BACKGROUND

Transmission of human immunodeficiency virus type 1 (HIV-1) drug resistance mutations, particularly that of minority drug-resistant variants, remains poorly understood. Population-based studies suggest that drug-resistant HIV-1 is less transmissible than drug-susceptible viruses. We compared HIV-1 drug-resistant genotypes among partner-pairs in order to assess the likelihood of transmission of drug resistance mutations and investigate the role of minority variants in HIV transmission.

METHODS AND FINDINGS

From 1992-2010, 340 persons with primary HIV-1 infection and their partners were enrolled into observational research studies at the University of Washington Primary Infection Clinic (UWPIC). Out of 50 partner-pairs enrolled, 36 (72%) transmission relationships were confirmed by phylogenetic distance analysis of HIV-1 envelope (env) sequences, and 31 partner-pairs enrolled after 1995 met criteria for this study. Drug resistance mutations in the region of the HIV-1 polymerase gene (pol) that encodes protease and reverse transcriptase were assessed by 454-pyrosequencing. In 25 partner-pairs where the transmission direction could be determined, 12 (48%) transmitters had 1-4 drug resistance mutations (23 total) detected in their HIV-1 populations at a median frequency of 6.0% (IQR 1.5%-98.7%, range 1.0%-99.6%). Of 10 major mutations detected in five transmitters at a frequency >95%, 100% (95% CI 69.2%-100%) were detected in recipients. All of these transmitters were antiretroviral (ARV)-naïve at the time of specimen collection. Fourteen mutations (eight major mutations and six accessory mutations) were detected in nine transmitters at low frequencies (1.0%-11.8%); four of these transmitters had previously received ARV therapy. Two (14% [95% CI 1.8%-42.8%]) G73S accessory mutations were detected in both transmitter and recipient. This number is not significantly different from the number expected based on the observed frequencies of drug-resistant viruses in transmitting partners. Limitations of this study include the small sample size and uncertainties in determining the timing of virus transmission and mutation history.

CONCLUSIONS

Drug-resistant majority variants appeared to be commonly transmitted by ARV-naïve participants in our analysis and may contribute significantly to transmitted drug resistance on a population level. When present at low frequency, no major mutation was observed to be shared between partner-pairs; identification of accessory mutations shared within a pair could be due to transmission, laboratory artifact, or apolipoprotein B mRNA-editing enzyme, catalytic polypeptides (APOBECs), and warrants further study.

Analytical parameters and validation of homopolymer detection in a pyrosequencing-based next generation sequencing system

Background

Current technologies in next-generation sequencing are offering high throughput reads at low costs, but still suffer from various sequencing errors. Although pyro- and ion semiconductor sequencing both have the advantage of delivering long and high quality reads, problems might occur when sequencing homopolymer-containing regions, since the repeating identical bases are going to incorporate during the same synthesis cycle, which leads to uncertainty in base calling. The aim of this study was to evaluate the analytical performance of a pyrosequencing-based next-generation sequencing system in detecting homopolymer sequences using homopolymer-preintegrated plasmid constructs and human DNA samples originating from patients with cystic fibrosis.

Results

In the plasmid system average correct genotyping was 95.8% in 4-mers, 87.4% in 5-mers and 72.1% in 6-mers. Despite the experienced low genotyping accuracy in 5- and 6-mers, it was possible to generate amplicons with more than a 90% adequate detection rate in every homopolymer tract. When homopolymers in the CFTR gene were sequenced average accuracy was 89.3%, but varied in a wide range (52.2 – 99.1%). In all but one case, an optimal amplicon-sequencing primer combination could be identified. In that single case (7A tract in exon 14 (c.2046_2052)), none of the tested primer sets produced the required analytical performance.

Conclusions

Our results show that pyrosequencing is the most reliable in case of 4-mers and as homopolymer length gradually increases, accuracy deteriorates. With careful primer selection, the NGS system was able to correctly genotype all but one of the homopolymers in the CFTR gene. In conclusion, we configured a plasmid test system that can be used to assess genotyping accuracy of NGS devices and developed an accurate NGS assay for the molecular diagnosis of CF using self-designed primers for amplification and sequencing.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4544-x) contains supplementary material, which is available to authorized users.

Short Communication: Analysis of Minor Populations of Human Immunodeficiency Virus by Primer Identification and Insertion-Deletion and Carry Forward Correction Pipelines

Accurate analysis of minor populations of drug-resistant HIV requires analysis of a sufficient number of viral templates. We assessed the effect of experimental conditions on the analysis of HIV pol 454 pyrosequences generated from plasma using (1) the "Insertion-deletion (indel) and Carry Forward Correction" (ICC) pipeline, which clusters sequence reads using a nonsubstitution approach and can correct for indels and carry forward errors, and (2) the "Primer Identification (ID)" method, which facilitates construction of a consensus sequence to correct for sequencing errors and allelic skewing. The Primer ID and ICC methods produced similar estimates of viral diversity, but differed in the number of sequence variants generated. Sequence preparation for ICC was comparably simple, but was limited by an inability to assess the number of templates analyzed and allelic skewing. The more costly Primer ID method corrected for allelic skewing and provided the number of viral templates analyzed, which revealed that amplifiable HIV templates varied across specimens and did not correlate with clinical viral load. This latter observation highlights the value of the Primer ID method, which by determining the number of templates amplified, enables more accurate assessment of minority species in the virus population, which may be relevant to prescribing effective antiretroviral therapy.

Comparison of Major and Minor Viral SNPs Identified through Single Template Sequencing and Pyrosequencing in Acute HIV-1 Infection

Massively parallel sequencing (MPS) technologies, such as 454-pyrosequencing, allow for the identification of variants in sequence populations at lower levels than consensus sequencing and most single-template Sanger sequencing experiments. We sought to determine if the greater depth of population sampling attainable using MPS technology would allow detection of minor variants in HIV founder virus populations very early in infection in instances where Sanger sequencing detects only a single variant. We compared single nucleotide polymorphisms (SNPs) during acute HIV-1 infection from 32 subjects using both single template Sanger and 454-pyrosequencing. Pyrosequences from a median of 2400 viral templates per subject and encompassing 40% of the HIV-1 genome, were compared to a median of five individually amplified near full-length viral genomes sequenced using Sanger technology. There was no difference in the consensus nucleotide sequences over the 3.6kb compared in 84% of the subjects infected with single founders and 33% of subjects infected with multiple founder variants: among the subjects with disagreements, mismatches were found in less than 1% of the sites evaluated (of a total of nearly 117,000 sites across all subjects). The majority of the SNPs observed only in pyrosequences were present at less than 2% of the subject’s viral sequence population. These results demonstrate the utility of the Sanger approach for study of early HIV infection and provide guidance regarding the design, utility and limitations of population sequencing from variable template sources, and emphasize parameters for improving the interpretation of massively parallel sequencing data to address important questions regarding target sequence evolution.

Resistance detected by pyrosequencing following zidovudine monotherapy for prevention of HIV-1 mother-to-child-transmission

To prevent mother-to-child-transmission of HIV-1, the 2010 WHO guidelines recommended prenatal zidovudine (ZDV) monotherapy (option A). To determine if ZDV monotherapy selects for HIV resistance in antiretroviral-naive women during pregnancy, specimens from 50 individuals were examined using pyrosequencing. ZDV-resistance mutations were detected at delivery in seven women (14%, 95% confidence interval 6.6-26.5%). These data raise the question whether women administered ZDV monotherapy for prevention of mother-to-child-transmission could have higher risk of virologic failure when later started on combination antiretroviral therapy, as has been demonstrated following single-dose nevirapine prophylaxis.

Fitness-Balanced Escape Determines Resolution of Dynamic Founder Virus Escape Processes in HIV-1 Infection

To understand the interplay between host cytotoxic T-lymphocyte (CTL) responses and the mechanisms by which HIV-1 evades them, we studied viral evolutionary patterns associated with host CTL responses in six linked transmission pairs. HIV-1 sequences corresponding to full-length p17 and p24 gag were generated by 454 pyrosequencing for all pairs near the time of transmission, and seroconverting partners were followed for a median of 847 days postinfection. T-cell responses were screened by gamma interferon/interleukin-2 (IFN-γ/IL-2) FluoroSpot using autologous peptide sets reflecting any Gag variant present in at least 5% of sequence reads in the individual's viral population. While we found little evidence for the occurrence of CTL reversions, CTL escape processes were found to be highly dynamic, with multiple epitope variants emerging simultaneously. We found a correlation between epitope entropy and the number of epitope variants per response (r = 0.43; P = 0.05). In cases in which multiple escape mutations developed within a targeted epitope, a variant with no fitness cost became fixed in the viral population. When multiple mutations within an epitope achieved fitness-balanced escape, these escape mutants were each maintained in the viral population. Additional mutations found to confer escape but undetected in viral populations incurred high fitness costs, suggesting that functional constraints limit the available sites tolerable to escape mutations. These results further our understanding of the impact of CTL escape and reversion from the founder virus in HIV infection and contribute to the identification of immunogenic Gag regions most vulnerable to a targeted T-cell attack.

IMPORTANCE

Rapid diversification of the viral population is a hallmark of HIV-1 infection, and understanding the selective forces driving the emergence of viral variants can provide critical insight into the interplay between host immune responses and viral evolution. We used deep sequencing to comprehensively follow viral evolution over time in six linked HIV transmission pairs. We then mapped T-cell responses to explore if mutations arose due to adaption to the host and found that escape processes were often highly dynamic, with multiple mutations arising within targeted epitopes. When we explored the impact of these mutations on replicative capacity, we found that dynamic escape processes only resolve with the selection of mutations that conferred escape with no fitness cost to the virus. These results provide further understanding of the complicated viral-host interactions that occur during early HIV-1 infection and may help inform the design of future vaccine immunogens.

Oligonucleotide ligation assay detects HIV drug resistance associated with virologic failure among antiretroviral-naive adults in Kenya

BACKGROUND

Transmitted drug resistance (TDR) is increasing in some areas of Africa. Detection of TDR may predict virologic failure of first-line nonnucleoside reverse transcriptase inhibitor (NNRTI)-based antiretroviral therapy (ART). We evaluated the utility of a relatively inexpensive oligonucleotide ligation assay (OLA) to detect clinically relevant TDR at the time of ART initiation.

METHODS

Pre-ART plasmas from ART-naive Kenyans initiating an NNRTI-based fixed-dose combination ART in a randomized adherence trial conducted in 2006 were retrospectively analyzed by OLA for mutations conferring resistance to NNRTI (K103N, Y181C, and G190A) and lamivudine (M184V). Post-ART plasmas were analyzed for virologic failure (≥1000 copies/mL) at 6-month intervals over 18-month follow-up. Pre-ART plasmas of those with virologic failure were evaluated for drug resistance by consensus and 454-pyrosequencing.

RESULTS

Among 386 participants, TDR was detected by OLA in 3.89% (95% confidence interval: 2.19 to 6.33) and was associated with a 10-fold higher rate of virologic failure (hazard ratio: 10.39; 95% confidence interval: 3.23 to 32.41; P < 0.001) compared with those without TDR. OLA detected 24 TDR mutations (K103N: n = 13; Y181C: n = 5; G190A: n = 3; M184V: n = 3) in 15 subjects (NNRTI: n = 15; 3TC: n = 3). Among 51 participants who developed virologic failure, consensus sequencing did not detect additional TDR mutations conferring high-level resistance, and pyrosequencing only detected additional mutations at frequencies <2%. Mutant frequencies <2% at ART initiation were significantly less likely to be found at the time of virologic failure compared with frequencies ≥2% (22% vs. 63%; P < 0.001).

CONCLUSIONS

Detection of TDR by a point mutation assay may prevent the use of suboptimal ART.

Toward a more accurate view of human B-cell repertoire by next-generation sequencing, unbiased repertoire capture and single-molecule barcoding

B-cell repertoire analysis using next-generation sequencing has become a valuable tool for interrogating the genetic record of humoral response to infection. However, key obstacles such as low throughput, short read length, high error rate, and undetermined bias of multiplex PCR method have hindered broader application of this technology. In this study, we report several technical advances in antibody repertoire sequencing. We first demonstrated the ability to sequence antibody variable domains using the Ion Torrent PGM platform. As a test case, we analyzed the PGT121 class of antibodies from IAVI donor 17, an HIV-1-infected individual. We then obtained "unbiased" antibody repertoires by sequencing the 5'-RACE PCR products of B-cell transcripts from IAVI donor 17 and two HIV-1-uninfected individuals. We also quantified the bias of previously published gene-specific primers by comparing the repertoires generated by 5'-RACE PCR and multiplex PCR. We further developed a single-molecule barcoding strategy to reduce PCR-based amplification noise. Lastly, we evaluated several new PGM technologies in the context of antibody sequencing. We expect that, based upon long-read and high-fidelity next-generation sequencing technologies, the unbiased analysis will provide a more accurate view of the overall antibody repertoire while the barcoding strategy will facilitate high-resolution analysis of individual antibody families.

Validation of an oligonucleotide ligation assay for quantification of human immunodeficiency virus type 1 drug-resistant mutants by use of massively parallel sequencing

Global HIV treatment programs need sensitive and affordable tests to monitor HIV drug resistance. We compared mutant detection by the oligonucleotide ligation assay (OLA), an economical and simple test, to massively parallel sequencing. Nonnucleoside reverse transcriptase inhibitor (K103N, V106M, Y181C, and G190A) and lamivudine (M184V) resistance mutations were quantified in blood-derived plasma RNA and cell DNA specimens by OLA and 454 pyrosequencing. A median of 1,000 HIV DNA or RNA templates (range, 163 to 1,874 templates) from blood specimens collected in Mozambique (n = 60) and Kenya (n = 51) were analyzed at 4 codons in each sample (n = 441 codons assessed). Mutations were detected at 75 (17%) codons by OLA sensitive to 2.0%, at 71 codons (16%; P = 0.78) by pyrosequencing using a cutoff value of ≥ 2.0%, and at 125 codons (28%; P < 0.0001) by pyrosequencing sensitive to 0.1%. Discrepancies between the assays included 15 codons with mutant concentrations of ∼2%, one at 8.8% by pyrosequencing and not detected by OLA, and one at 69% by OLA and not detected by pyrosequencing. The latter two cases were associated with genetic polymorphisms in the regions critical for ligation of the OLA probes and pyrosequencing primers, respectively. Overall, mutant concentrations quantified by the two methods correlated well across the codons tested (R(2) > 0.8). Repeat pyrosequencing of 13 specimens showed reproducible detection of 5/24 mutations at <2% and 6/6 at ≥ 2%. In conclusion, the OLA and pyrosequencing performed similarly in the quantification of nonnucleoside reverse transcriptase inhibitor and lamivudine mutations present at >2% of the viral population in clinical specimens. While pyrosequencing was more sensitive, detection of mutants below 2% was not reproducible.

Improved detection of rare HIV-1 variants using 454 pyrosequencing

454 pyrosequencing, a massively parallel sequencing (MPS) technology, is often used to study HIV genetic variation. However, the substantial mismatch error rate of the PCR required to prepare HIV-containing samples for pyrosequencing has limited the detection of rare variants within viral populations to those present above ~1%. To improve detection of rare variants, we varied PCR enzymes and conditions to identify those that combined high sensitivity with a low error rate. Substitution errors were found to vary up to 3-fold between the different enzymes tested. The sensitivity of each enzyme, which impacts the number of templates amplified for pyrosequencing, was shown to vary, although not consistently across genes and different samples. We also describe an amplicon-based method to improve the consistency of read coverage over stretches of the HIV-1 genome. Twenty-two primers were designed to amplify 11 overlapping amplicons in the HIV-1 clade B gag-pol and env gp120 coding regions to encompass 4.7 kb of the viral genome per sample at sensitivities as low as 0.01-0.2%.