Analysis of 454 sequencing error rate, error sources, and artifact recombination for detection of Low-frequency drug resistance mutations in HIV-1 DNA

Quantitative Prediction of Human Immunodeficiency Virus Drug Resistance

Drug resistance of pathogens, including viruses, is one of the reasons for decreased efficacy of therapy. Considering the impact of HIV type 1 (HIV-1) on the development of progressive immune dysfunction and the rapid development of drug resistance, the analysis of HIV-1 resistance is of high significance. Currently, a substantial amount of data has been accumulated on HIV-1 drug resistance that can be used to build both qualitative and quantitative models of HIV-1 drug resistance. Quantitative models of drug resistance can enrich the information about the efficacy of a particular drug in the scheme of antiretroviral therapy. In our study, we investigated the possibility of developing models for quantitative prediction of HIV-1 resistance to eight protease inhibitors based on the analysis of amino acid sequences of HIV-1 protease for 900 virus variants. We developed random forest regression (RFR), support vector regression (SVR), and self-consistent regression (SCR) models using binary vectors containing values from 0 or 1, depending on the presence of a specific peptide fragment in each amino acid sequence as independent variables, while fold ratio, reflecting the level of resistance, was the predicted variable. The SVR and SCR models showed the highest predictive performances. The models built demonstrate reasonable performances for eight out of nine (R2 varied from 0.828 to 0.909) protease inhibitors, while R2 for predicting tipranavir fold ratio was lower (R2 was 0.642). We believe that the developed approach can be applied to evaluate drug resistance of molecular targets of other viruses where appropriate experimental data are available.

Dual UMIs and Dual Barcodes With Minimal PCR Amplification Removes Artifacts and Acquires Accurate Antibody Repertoire

Antibody repertoire sequencing (Rep-seq) has been widely used to reveal repertoire dynamics and to interrogate antibodies of interest at single nucleotide-level resolution. However, polymerase chain reaction (PCR) amplification introduces extensive artifacts including chimeras and nucleotide errors, leading to false discovery of antibodies and incorrect assessment of somatic hypermutations (SHMs) which subsequently mislead downstream investigations. Here, a novel approach named DUMPArts, which improves the accuracy of antibody repertoires by labeling each sample with dual barcodes and each molecule with dual unique molecular identifiers (UMIs) via minimal PCR amplification to remove artifacts, is developed. Tested by ultra-deep Rep-seq data, DUMPArts removed inter-sample chimeras, which cause artifactual shared clones and constitute approximately 15% of reads in the library, as well as intra-sample chimeras with erroneous SHMs and constituting approximately 20% of the reads, and corrected base errors and amplification biases by consensus building. The removal of these artifacts will provide an accurate assessment of antibody repertoires and benefit related studies, especially mAb discovery and antibody-guided vaccine design.

Chemistry and Bioinformatics Considerations in Using Next-Generation Sequencing Technologies to Inferring HIV Proviral DNA Genome-Intactness

HIV persists via integration of the viral DNA into the human genome. The HIV DNA pool within an infected individual is a complex population that comprises both intact and defective viral genomes, each with a distinct integration site, in addition to a unique repertoire of viral quasi-species. Obtaining an accurate profile of the viral DNA pool is critical to understanding viral persistence and resolving interhost differences. Recent advances in next-generation deep sequencing (NGS) technologies have enabled the development of two sequencing assays to capture viral near-full- genome sequences at single molecule resolution (FLIP-seq) or to co-capture full-length viral genome sequences in conjunction with its associated viral integration site (MIP-seq). This commentary aims to provide an overview on both FLIP-seq and MIP-seq, discuss their strengths and limitations, and outline specific chemistry and bioinformatics concerns when using these assays to study HIV persistence.

Accurate reconstruction of bacterial pan- and core genomes with PEPPAN

Bacterial genomes can contain traces of a complex evolutionary history, including extensive homologous recombination, gene loss, gene duplications, and horizontal gene transfer. To reconstruct the phylogenetic and population history of a set of multiple bacteria, it is necessary to examine their pangenome, the composite of all the genes in the set. Here we introduce PEPPAN, a novel pipeline that can reliably construct pangenomes from thousands of genetically diverse bacterial genomes that represent the diversity of an entire genus. PEPPAN outperforms existing pangenome methods by providing consistent gene and pseudogene annotations extended by similarity-based gene predictions, and identifying and excluding paralogs by combining tree- and synteny-based approaches. The PEPPAN package additionally includes PEPPAN_parser, which implements additional downstream analyses, including the calculation of trees based on accessory gene content or allelic differences between core genes. To test the accuracy of PEPPAN, we implemented SimPan, a novel pipeline for simulating the evolution of bacterial pangenomes. We compared the accuracy and speed of PEPPAN with four state-of-the-art pangenome pipelines using both empirical and simulated data sets. PEPPAN was more accurate and more specific than any of the other pipelines and was almost as fast as any of them. As a case study, we used PEPPAN to construct a pangenome of approximately 40,000 genes from 3052 representative genomes spanning at least 80 species of Streptococcus The resulting gene and allelic trees provide an unprecedented overview of the genomic diversity of the entire Streptococcus genus.

Short Communication: HIV-DRLink: A Tool for Reporting Linked HIV-1 Drug Resistance Mutations in Large Single-Genome Data Sets Using the Stanford HIV Database

The prevalence of HIV-1 drug resistance is increasing worldwide and monitoring its emergence is important for the successful management of populations receiving combination antiretroviral therapy. It is likely that pre-existing drug resistance mutations linked on the same viral genomes are predictive of treatment failure. Because of the large number of sequences generated by ultrasensitive single-genome sequencing (uSGS) and other similar next-generation sequencing methods, it is difficult to assess each sequence individually for linked drug resistance mutations. Several software/programs exist to report the frequencies of individual mutations in large data sets, but they provide no information on linkage of resistance mutations. In this study, we report the HIV-DRLink program, a research tool that provides resistance mutation frequencies as well as their genetic linkage by parsing and summarizing the Sierra output from the Stanford HIV Database. The HIV-DRLink program should only be used on data sets generated by methods that eliminate artifacts due to polymerase chain reaction recombination, for example, standard single-genome sequencing or uSGS. HIV-DRLink is exclusively a research tool and is not intended to inform clinical decisions.

Comparison of sequencing methods and data processing pipelines for whole genome sequencing and minority single nucleotide variant (mSNV) analysis during an influenza A/H5N8 outbreak

As high-throughput sequencing technologies are becoming more widely adopted for analysing pathogens in disease outbreaks there needs to be assurance that the different sequencing technologies and approaches to data analysis will yield reliable and comparable results. Conversely, understanding where agreement cannot be achieved provides insight into the limitations of these approaches and also allows efforts to be focused on areas of the process that need improvement. This manuscript describes the next-generation sequencing of three closely related viruses, each analysed using different sequencing strategies, sequencing instruments and data processing pipelines. In order to determine the comparability of consensus sequences and minority (sub-consensus) single nucleotide variant (mSNV) identification, the biological samples, the sequence data from 3 sequencing platforms and the *.bam quality-trimmed alignment files of raw data of 3 influenza A/H5N8 viruses were shared. This analysis demonstrated that variation in the final result could be attributed to all stages in the process, but the most critical were the well-known homopolymer errors introduced by 454 sequencing, and the alignment processes in the different data processing pipelines which affected the consistency of mSNV detection. However, homopolymer errors aside, there was generally a good agreement between consensus sequences that were obtained for all combinations of sequencing platforms and data processing pipelines. Nevertheless, minority variant analysis will need a different level of careful standardization and awareness about the possible limitations, as shown in this study.

High-throughput sequencing (HTS) for the analysis of viral populations

The development of High-Throughput Sequencing (HTS) technologies is having a major impact on the genomic analysis of viral populations. Current HTS platforms can capture nucleic acid variation across millions of genes for both selected amplicons and full viral genomes. HTS has already facilitated the discovery of new viruses, hinted new taxonomic classifications and provided a deeper and broader understanding of their diversity, population and genetic structure. Hence, HTS has already replaced standard Sanger sequencing in basic and applied research fields, but the next step is its implementation as a routine technology for the analysis of viruses in clinical settings. The most likely application of this implementation will be the analysis of viral genomics, because the huge population sizes, high mutation rates and very fast replacement of viral populations have demonstrated the limited information obtained with Sanger technology. In this review, we describe new technologies and provide guidelines for the high-throughput sequencing and genetic and evolutionary analyses of viral populations and metaviromes, including software applications. With the development of new HTS technologies, new and refurbished molecular and bioinformatic tools are also constantly being developed to process and integrate HTS data. These allow assembling viral genomes and inferring viral population diversity and dynamics. Finally, we also present several applications of these approaches to the analysis of viral clinical samples including transmission clusters and outbreak characterization.

PacBio Amplicon Sequencing Method To Measure Pilin Antigenic Variation Frequencies of Neisseria gonorrhoeae

Gene diversification is a common mechanism pathogens use to alter surface structures to aid in immune avoidance. Neisseria gonorrhoeae uses a gene conversion-based diversification system to alter the primary sequence of the gene encoding the major subunit of the pilus, pilE Antigenic variation occurs when one of the nonexpressed 19 silent copies donates part of its DNA sequence to pilE We have developed a method using Pacific Biosciences (PacBio) amplicon sequencing and custom software to determine pilin antigenic variation frequencies. The program analyzes 37 variable regions across the strain FA1090 1-81-S2 pilE gene and can be modified to determine sequence variation from other starting pilE sequences or other diversity generation systems. Using this method, we measured pilin antigenic variation frequencies for various derivatives of strain FA1090 and showed we can also analyze pilin antigenic variation frequencies during macrophage infection.IMPORTANCE Diversity generation systems are used by many unicellular organism to provide subpopulations of cell with different properties that are available when needed. We have developed a method using the PacBio DNA sequencing technology and a custom computer program to analyze the pilin antigenic variation system of the organism that is the sole cause of the sexually transmitted infection, gonorrhea.

Proliferative memory SAMHD1low CD4+ T cells harbour high levels of HIV-1 with compartmentalized viral populations

We previously reported the presence of memory CD4⁺ T cells that express low levels of SAMHD1 (SAMHD1^low) in peripheral blood and lymph nodes from both HIV-1 infected and uninfected individuals. These cells are enriched in Th17 and Tfh subsets, two populations known to be preferentially targeted by HIV-1. Here we investigated whether SAMHD1^low CD4⁺ T-cells harbour replication-competent virus and compartimentalized HIV-1 genomes. We sorted memory CD4⁺CD45RO⁺SAMHD1^low, CD4⁺CD45RO⁺SAMHD1⁺ and naive CD4⁺CD45RO^-SAMHD1⁺ cells from HIV-1-infected patients on anti-retroviral therapy (c-ART) and performed HIV-1 DNA quantification, ultra-deep-sequencing of partial env (C2/V3) sequences and phenotypic characterization of the cells. We show that SAMHD1^low cells include novel Th17 CCR6⁺ subsets that lack CXCR3 and CCR4 (CCR6⁺DN). There is a decrease of the % of Th17 in SAMHD1^low compartment in infected compared to uninfected individuals (41% vs 55%, p<0.05), whereas the % of CCR6⁺DN increases (7.95% vs 3.8%, p<0.05). Moreover, in HIV-1 infected patients, memory SAMHD1^low cells harbour high levels of HIV-1 DNA compared to memory SAMHD1⁺ cells (4.5 vs 3.8 log/10⁶cells, respectively, p<0.001), while naïve SAMHD1⁺ showed significantly lower levels (3.1 log/10⁶cells, p<0.0001). Importantly, we show that SAMHD1^low cells contain p24-producing cells. Moreover, phylogenetic analyses revealed well-segregated HIV-1 DNA populations with compartmentalization between SAMHD1^low and SAMHD1⁺ memory cells, and limited viral exchange. As expected, the % of Ki67⁺ cells was significantly higher in SAMHD1^low compared to SAMHD1⁺ cells. There was positive association between levels of HIV-1 DNA and Ki67⁺ in memory SAMHD1^low cells, but not in memory and naïve SAMHD1⁺ CD4⁺ T-cells. Altogether, these data suggest that proliferative memory SAMHD1^low cells contribute to viral persistence.

In our previous results we reported that memory CD4⁺ T cells expressing low levels of SAMHD1 (SAMHD1^low) are present in peripheral blood and lymph nodes from HIV-1 infected and uninfected individuals. These cells were enriched in Th17 and Tfh, two populations targeted by HIV-1. Here we used purified memory CD4⁺CD45RO⁺SAMHD1^low, CD4⁺CD45RO⁺SAMHD1⁺ and naive CD4⁺CD45RO^-SAMHD1⁺ cells from HIV-1-infected and treated patients to perform cell-associated HIV-1 DNA quantification, p24-producing cells detection, ultra-deep-sequencing of partial env (C2/V3) HIV-1 DNA and further phenotypic characterization. Our results demonstrate that (i) Th17 and CCR6⁺DN-expressing transcriptional signature of early Th17, two major populations that are susceptible to HIV-1 infection, are present in SAMHD1^low cells, and while the former decreased significantly in c-ART HIV-1 infected compared to uninfected individuals, the latter significantly increased; (ii) memory SAMHD1^low cells from c-ART patients carry high levels of HIV-1 DNA compared to SAMHD1⁺ cells, and these levels positively and significantly correlated with Ki67 expression; (iii) memory SAMHD1^low cells from patients harbour p24-producing cells; (iv) phylogenetic analyses revealed well-segregated HIV-1 DNA populations with significant compartmentalization between SAMHD1^low and SAMHD1⁺ cells and limited viral exchange. Our data demonstrate that memory SAMHD1^low cells contribute to HIV-1 persistence.

PCR Amplification Strategies Towards Full-length HIV-1 Genome Sequencing

The advent of next-generation sequencing has enabled greater resolution of viral diversity and improved feasibility of full viral genome sequencing allowing routine HIV-1 full genome sequencing in both research and diagnostic settings. Regardless of the sequencing platform selected, successful PCR amplification of the HIV-1 genome is essential for sequencing template preparation. As such, full HIV-1 genome amplification is a crucial step in dictating the successful and reliable sequencing downstream. Here we reviewed existing PCR protocols leading to HIV-1 full genome sequencing. In addition to the discussion on basic considerations on relevant PCR design, the advantages as well as the pitfalls of the published protocols were reviewed.

Ultra-deep amplicon sequencing indicates absence of low-grade mosaicism with normal cells in patients with type-1 NF1 deletions

Different types of large NF1 deletion are distinguishable by breakpoint location and potentially also by the frequency of mosaicism with normal cells lacking the deletion. However, low-grade mosaicism with fewer than 10% normal cells has not yet been excluded for all NF1 deletion types since it is impossible to assess by the standard techniques used to identify such deletions, including MLPA and array analysis. Here, we used ultra-deep amplicon sequencing to investigate the presence of normal cells in the blood of 20 patients with type-1 NF1 deletions lacking mosaicism according to MLPA. The ultra-deep sequencing entailed the screening of 96 amplicons for heterozygous SNVs located within the NF1 deletion region. DNA samples from three previously identified patients with type-2 NF1 deletions and low-grade mosaicism with normal cells as determined by FISH or microsatellite marker analysis were used to validate our methodology. In these type-2 NF1 deletion samples, proportions of 5.3%, 6.6% and 15.0% normal cells, respectively, were detected by ultra-deep amplicon sequencing. However, using this highly sensitive method, none of the 20 patients with type-1 NF1 deletions included in our analysis exhibited low-grade mosaicism with normal cells in blood, thereby supporting the view that the vast majority of type-1 deletions are germline deletions.

Quasispecies dynamics and clinical significance of hepatitis C virus (HCV) antiviral resistance

Hepatitis C virus (HCV) follows quasispecies dynamics in infected hosts and this influences its biology, how the virus diversifies into several genotypes and many subtypes, and how viral populations respond to antiviral therapies. Despite current antiviral combinations being able to cure a great percentage of HCV-infected patients, the presence of resistance-associated substitutions (RASs) diminishes the success of antiviral therapies, which is a main concern in the re-treatment of patients treated with direct-acting antiviral agents. Current methodologies such as ultra deep sequencing are ideal tools to obtain a detailed representation of the mutant spectrum composition circulating in infected patients. Such knowledge should allow optimisation of rescue treatments. A new mechanism of antiviral resistance not based on the selection of RASs but on high viral fitness is discussed.

Pipeline for specific subtype amplification and drug resistance detection in hepatitis C virus

BACKGROUND

Despite the high sustained virological response rates achieved with current directly-acting antiviral agents (DAAs) against hepatitis C virus (HCV), around 5-10% of treated patients do not respond to current antiviral therapies, and basal resistance to DAAs is increasingly detected among treatment-naïve infected individuals. Identification of amino acid substitutions (including those in minority variants) associated with treatment failure requires analytical designs that take into account the high diversification of HCV in more than 86 subtypes according to the ICTV website (June 2017).

METHODS

The methodology has involved five sequential steps: (i) to design 280 oligonucleotide primers (some including a maximum of three degenerate positions), and of which 120 were tested to amplify NS3, NS5A-, and NS5B-coding regions in a subtype-specific manner, (ii) to define a reference sequence for each subtype, (iii) to perform experimental controls to define a cut-off value for detection of minority amino acids, (iv) to establish bioinformatics' tools to quantify amino acid replacements, and (v) to validate the procedure with patient samples.

RESULTS

A robust ultra-deep sequencing procedure to analyze HCV circulating in serum samples from patients infected with virus that belongs to the ten most prevalent subtypes worldwide: 1a, 1b, 2a, 2b, 2c, 2j, 3a, 4d, 4e, 4f has been developed. Oligonucleotide primers are subtype-specific. A cut-off value of 1% mutant frequency has been established for individual mutations and haplotypes.

CONCLUSION

The methodological pipeline described here is adequate to characterize in-depth mutant spectra of HCV populations, and it provides a tool to understand HCV diversification and treatment failures. The pipeline can be periodically extended in the event of HCV diversification into new genotypes or subtypes, and provides a framework applicable to other RNA viral pathogens, with potential to couple detection of drug-resistant mutations with treatment planning.

Evaluating the accuracy and sensitivity of detecting minority HIV-1 populations by Illumina next-generation sequencing

The accuracy and sensitivity of deep sequencing were assessed using viral standards (pNL4-3 and pLAI.2) of both DNA and RNA. The sequencing accuracy did not depend on the type of nucleic acid, but critically depended on the number of reads and threshold of sensitivity to minor viral populations. With coverage of more than 236 reads, the accuracy of viral RNA sequencing was equal to or exceeded 99.9%, with a sensitivity threshold to minor nucleotides of 20%. When the sensitivity threshold was below 1%, reduced accuracy dynamics were clearly visible even when the coverage was massive (more than 9.000 reads). It was found that the floating sensitivity threshold allowed the sequencing accuracy to be maintained at an acceptable level in cases of low coverage (less than 1.500-2.000) of reads. These results indicate the quality that can be expected with a specific number of reads and sensitivity threshold. Deep sequencing is a very powerful tool that can significantly improve the value of study results, but despite its superior performance, it should be used with caution regarding its sensitivity to minor populations below 1%.

Analysis of drug resistance mutations in whole blood DNA from HIV-1 infected patients by single genome and ultradeep sequencing analysis

In HIV-1 infected patients blood CD4⁺ T lymphocytes could be a valuable target to analyse drug resistance mutations (DRM) selected over the course of the infection. However, detection of viral resistance mutations in cellular DNA by standard genotype resistance techniques (SGRT) is suboptimal. Whole blood DNA (wbDNA) from 12 HIV-1 infected patients on ART was studied by Single Genome Sequencing (SGS) and 8 of them also by Ultradeep pyrosequencing (UDP). Results were compared with contemporary and historical DRM detected in plasma by SGRT during follow up. All the contemporary DRM detected in plasma from the viremic patients were detected by SGS and UDP (20 from 7 patients and 4 from 5 patients respectively). Out of the 67 historical DRM detected in plasma and no longer present at the time of testing, 38 (57%) were detected by SGS in 12 patients and 27 out of 46 (59%) by UDP in 8 patients. Additional DRM never reported in plasma by SGRT were detected by SGS (12 from 8 patients) and UDP (10 from 8 patients). Analysis of wbDNA from HIV-1 infected patients by SGS and UDP provides proof of concept of the value of blood DNA to investigate current and archived DRM in HIV-1 infected patients on ART.

MERIT reveals the impact of genomic context on sequencing error rate in ultra-deep applications

BACKGROUND

Rapid progress in high-throughput sequencing (HTS) and the development of novel library preparation methods have improved the sensitivity of detecting mutations in heterogeneous samples, specifically in high-depth (> 500×) clinical applications. However, HTS methods are bounded by their technical and theoretical limitations and sequencing errors cannot be completely eliminated. Comprehensive quantification of the background noise can highlight both the efficiency and the limitations of any HTS methodology, and help differentiate true mutations at low abundance from artifacts.

RESULTS

We introduce MERIT (Mutation Error Rate Inference Toolkit), designed for in-depth quantification of erroneous substitutions and small insertions and deletions. MERIT incorporates an all-inclusive variant caller and considers genomic context, including the nucleotides immediately at 5 'and 3 ', thereby establishing error rates for 96 possible substitutions as well as four single-base and 16 double-base indels. We applied MERIT to ultra-deep sequencing data (1,300,000 ×) obtained from the amplification of multiple clinically relevant loci, and showed a significant relationship between error rates and genomic contexts. In addition to observing significant difference between transversion and transition rates, we identified variations of more than 100-fold within each error type at high sequencing depths. For instance, T >G transversions in trinucleotide GTCs occurred 133.5 ± 65.9 more often than those in ATAs. Similarly, C >T transitions in GCGs were observed at 73.8 ± 10.5 higher rate than those in TCTs. We also devised an in silico approach to determine the optimal sequencing depth, where errors occur at rates similar to those of expected true mutations. Our analyses showed that increasing sequencing depth might improve sensitivity for detecting some mutations based on their genomic context. For example, T >G rate of error in GTCs did not change when sequenced beyond 10,000 ×; in contrast, T >G rate in TTAs consistently improved even at above 500,000 ×.

CONCLUSIONS

Our results demonstrate significant variation in nucleotide misincorporation rates, and suggest that genomic context should be considered for comprehensive profiling of specimen-specific and sequencing artifacts in high-depth assays. This data provide strong evidence against assigning a single allele frequency threshold to call mutations, for it can result in substantial false positive as well as false negative variants, with important clinical consequences.

An innovative diagnostic technology for the codon mutation C580Y in kelch13 of Plasmodium falciparum with MinION nanopore sequencer

BACKGROUND

The recent spread of artemisinin (ART)-resistant Plasmodium falciparum represents an emerging global threat to public health. In Southeast Asia, the C580Y mutation of kelch13 (k13) is the dominant mutation of ART-resistant P. falciparum. Therefore, a simple method for the detection of C580Y mutation is urgently needed to enable widespread routine surveillance in the field. The aim of this study is to develop a new diagnostic procedure for the C580Y mutation using loop-mediated isothermal amplification (LAMP) combined with the MinION nanopore sequencer.

RESULTS

A LAMP assay for the k13 gene of P. falciparum to detect the C580Y mutation was successfully developed. The detection limit of this procedure was 10 copies of the reference plasmid harboring the k13 gene within 60 min. Thereafter, amplicon sequencing of the LAMP products using the MinION nanopore sequencer was performed to clarify the nucleotide sequences of the gene. The C580Y mutation was identified based on the sequence data collected from MinION reads 30 min after the start of sequencing. Further, clinical evaluation of the LAMP assay in 34 human blood samples collected from patients with P. falciparum malaria in Indonesia revealed a positive detection rate of 100%. All LAMP amplicons of up to 12 specimens were simultaneously sequenced using MinION. The results of sequencing were consistent with those of the conventional PCR and Sanger sequencing protocol. All procedures from DNA extraction to variant calling were completed within 3 h. The C580Y mutation was not found among these 34 P. falciparum isolates in Indonesia.

CONCLUSIONS

An innovative method combining LAMP and MinION will enable simple, rapid, and high-sensitivity detection of the C580Y mutation of P. falciparum, even in resource-limited situations in developing countries.

A Survey of Virus Recombination Uncovers Canonical Features of Artificial Chimeras Generated During Deep Sequencing Library Preparation

Chimeric reads can be generated by in vitro recombination during the preparation of high-throughput sequencing libraries. Our attempt to detect biological recombination between the genomes of dengue virus (DENV; +ssRNA genome) and its mosquito host using the Illumina Nextera sequencing library preparation kit revealed that most, if not all, detected host-virus chimeras were artificial. Indeed, these chimeras were not more frequent than with control RNA from another species (a pillbug), which was never in contact with DENV RNA prior to the library preparation. The proportion of chimera types merely reflected those of the three species among sequencing reads. Chimeras were frequently characterized by the presence of 1-20 bp microhomology between recombining fragments. Within-species chimeras mostly involved fragments in opposite orientations and located less than 100 bp from each other in the parental genome. We found similar features in published datasets using two other viruses: Ebola virus (EBOV; -ssRNA genome) and a herpesvirus (dsDNA genome), both produced with the Illumina Nextera protocol. These canonical features suggest that artificial chimeras are generated by intra-molecular template switching of the DNA polymerase during the PCR step of the Nextera protocol. Finally, a published Illumina dataset using the Flock House virus (FHV; +ssRNA genome) generated with a protocol preventing artificial recombination revealed the presence of 1-10 bp microhomology motifs in FHV-FHV chimeras, but very few recombining fragments were in opposite orientations. Our analysis uncovered sequence features characterizing recombination breakpoints in short-read sequencing datasets, which can be helpful to evaluate the presence and extent of artificial recombination.

HIV evolution and diversity in ART-treated patients

Characterizing HIV genetic diversity and evolution during antiretroviral therapy (ART) provides insights into the mechanisms that maintain the viral reservoir during ART. This review describes common methods used to obtain and analyze intra-patient HIV sequence data, the accumulation of diversity prior to ART and how it is affected by suppressive ART, the debate on viral replication and evolution in the presence of ART, HIV compartmentalization across various tissues, and mechanisms for the emergence of drug resistance. It also describes how CD4+ T cells that were likely infected with latent proviruses prior to initiating treatment can proliferate before and during ART, providing a renewable source of infected cells despite therapy. Some expanded cell clones carry intact and replication-competent proviruses with a small fraction of the clonal siblings being transcriptionally active and a source for residual viremia on ART. Such cells may also be the source for viral rebound after interrupting ART. The identical viral sequences observed for many years in both the plasma and infected cells of patients on long-term ART are likely due to the proliferation of infected cells both prior to and during treatment. Studies on HIV diversity may reveal targets that can be exploited in efforts to eradicate or control the infection without ART.

Lower pre-ART intra-participant HIV-1 pol diversity may not be associated with virologic failure in adults

BACKGROUND

Identifying pre-ART factors associated with the emergence of HIV-1 drug resistance is critical for optimizing strategies to prevent virologic failure. A previous study reported that lower pre-ART HIV-1 pol diversity was associated with higher risk of virologic failure in HIV-1-infected children. To investigate this association in adults, we measured HIV-1 diversity with deep sequencing in pre-ART samples from adults with well-characterized virologic outcomes in a study (A5142) of initial ART conducted by the AIDS Clinical Trials Group (ACTG).

METHODS

We identified 22 cases in ACTG A5142 who experienced virologic failure with drug resistance mutations in RT and 44 matched controls who did not experience virologic failure. cDNA was synthesized from plasma HIV-1 RNA. Each cDNA molecule was tagged with a unique primer ID and RT codons 41-103 were amplified and deep sequenced. Sequences with the same tag were aligned and a consensus was generated to reduce PCR and sequencing errors. Diversity was calculated by measuring average pairwise distance (APD) of the consensus sequences. An exact conditional logistic regression model with percent APD as the risk factor estimated the odds ratio for VF and the corresponding 95% confidence interval.

RESULTS

Consensus single-genome sequences and diversity estimates of pol were obtained for pre-ART samples from 21 cases and 42 controls. The median (IQR) pre-ART percent APD was 0.71 (0.31-1.13) in cases and 0.58 (0.32-0.94) in controls. A possible trend was found for higher diversity being associated with greater risk of virologic failure in adults (OR = 2.2 per one percent APD increase, 95% CI = [0.8, 7.2]; p = 0.15).

CONCLUSIONS

This study in adults suggests there is a positive association between higher pre-ART pol diversity and the risk of virologic failure in adults rather than an inverse relationship reported in children.

Higher sequence diversity in the vaginal tract than in blood at early HIV-1 infection

In the majority of cases, human immunodeficiency virus type 1 (HIV-1) infection is transmitted through sexual intercourse. A single founder virus in the blood of the newly infected donor emerges from a genetic bottleneck, while in rarer instances multiple viruses are responsible for systemic infection. We sought to characterize the sequence diversity at early infection, between two distinct anatomical sites; the female reproductive tract vs. systemic compartment. We recruited 72 women from Uganda and Zimbabwe within seven months of HIV-1 infection. Using next generation deep sequencing, we analyzed the total genetic diversity within the C2-V3-C3 envelope region of HIV-1 isolated from the female genital tract at early infection and compared this to the diversity of HIV-1 in plasma. We then compared intra-patient viral diversity in matched cervical and blood samples with three or seven months post infection. Genetic analysis of the C2-V3-C3 region of HIV-1 env revealed that early HIV-1 isolates within blood displayed a more homogeneous genotype (mean 1.67 clones, range 1–5 clones) than clones in the female genital tract (mean 5.7 clones, range 3–10 clones) (p<0.0001). The higher env diversity observed within the genital tract compared to plasma was independent of HIV-1 subtype (A, C and D). Our analysis of early mucosal infections in women revealed high HIV-1 diversity in the vaginal tract but few transmitted clones in the blood. These novel in vivo finding suggest a possible mucosal sieve effect, leading to the establishment of a homogenous systemic infection.

During chronic HIV-1 infection, high viral diversity can be found in the blood and semen of donors. However, a single HIV-1 clone establishes productive infection in the recipient following heterosexual transmission. To investigate the genetic bottleneck occurring at the earliest stages of heterosexual HIV-1 transmission, we characterized the HIV-1 envelope sequence diversity at very early and early stages of infection in the female reproductive tract and matched plasma samples from a cohort of Ugandan and Zimbabwean women. A more diverse viral population was observed in the endocervical swab samples compared to plasma. Endocervical samples harbored a larger number of viral clones, while in the majority of plasma samples only a single clone was present early in infection. Interestingly, these observations were independent of HIV-1 subtype, hormonal contraceptive use or the number of sex acts and partners. Furthermore, in the cases of higher HIV-1 diversity in the blood during early infection, faster CD4 T cell decline were observed during chronic disease suggesting faster disease progression. Our findings provide novel in vivo evidence for the existence of an intra-patient genetic bottleneck restricting the HIV-1 from the vaginal tract to the blood during early heterosexual HIV-1 transmission.

Size, Composition, and Evolution of HIV DNA Populations during Early Antiretroviral Therapy and Intensification with Maraviroc

Residual viremia is common during antiretroviral therapy (ART) and could be caused by ongoing low-level virus replication or by release of viral particles from infected cells. ART intensification should impact ongoing viral propagation but not virion release. Eighteen acutely infected men were enrolled in a randomized controlled trial and monitored for a median of 107 weeks. Participants started ART with (n = 9) or without (n = 9) intensification with maraviroc (MVC) within 90 days of infection. Levels of HIV DNA and cell-free RNA were quantified by droplet digital PCR. Deep sequencing of C2-V3 env, gag, and pol (454 Roche) was performed on longitudinally collected plasma and peripheral blood mononuclear cell (PBMC) samples while on ART. Sequence data were analyzed for evidence of evolution by (i) molecular diversity analysis, (ii) nonparametric test for panmixia, and (iii) tip date randomization within a Bayesian framework. There was a longitudinal decay of HIV DNA after initiation of ART with no difference between MVC intensification groups (-0.08 ± 0.01 versus -0.09 ± 0.01 log₁₀ copies/week in MVC⁺ versus MVC^- groups; P = 0.62). All participants had low-level residual viremia (median, 2.8 RNA copies/ml). Across participants, medians of 56 (interquartile range [IQR], 36 to 74), 29 (IQR, 25 to 35), and 40 (IQR, 31 to 54) haplotypes were generated for env, gag, and pol regions, respectively. There was no clear evidence of viral evolution during ART and no difference in viral diversity or population structure from individuals with or without MVC intensification. Further efforts focusing on elucidating the mechanism(s) of viral persistence in various compartments using recent sequencing technologies are still needed, and potential low-level viral replication should always be considered in cure strategies.IMPORTANCE Residual viremia is common among HIV-infected people on ART. It remains controversial if this viremia is a consequence of propagating infection. We hypothesized that molecular evolution would be detectable during viral propagation and that therapy intensified with the entry inhibitor maraviroc would demonstrate less evolution. We performed a randomized double-blinded treatment trial with 18 acutely infected men (standard ART versus standard ART plus maraviroc). From longitudinally collected blood plasma and cells, levels of HIV DNA and cell-free HIV RNA were quantified by droplet digital PCR, and HIV DNA (env, gag, and pol coding regions) was deep sequenced (454 Roche). Investigating people who started ART during the earliest stages of their HIV infection, when viral diversity is low, provides an opportunity to detect evidence of viral evolution. Despite using a battery of analytical techniques, no clear and consistent evidence of viral propagation for over 90 weeks of observation could be discerned.

Transcriptional-Readthrough RNAs Reflect the Phenomenon of "A Gene Contains Gene(s)" or "Gene(s) within a Gene" in the Human Genome, and Thus Are Not Chimeric RNAs

Tens of thousands of chimeric RNAs, i.e., RNAs with sequences of two genes, have been identified in human cells. Most of them are formed by two neighboring genes on the same chromosome and are considered to be derived via transcriptional readthrough, but a true readthrough event still awaits more evidence and trans-splicing that joins two transcripts together remains as a possible mechanism. We regard those genomic loci that are transcriptionally read through as unannotated genes, because their transcriptional and posttranscriptional regulations are the same as those of already-annotated genes, including fusion genes formed due to genetic alterations. Therefore, readthrough RNAs and fusion-gene-derived RNAs are not chimeras. Only those two-gene RNAs formed at the RNA level, likely via trans-splicing, without corresponding genes as genomic parents, should be regarded as authentic chimeric RNAs. However, since in human cells, procedural and mechanistic details of trans-splicing have never been disclosed, we doubt the existence of trans-splicing. Therefore, there are probably no authentic chimeras in humans, after readthrough and fusion-gene derived RNAs are all put back into the group of ordinary RNAs. Therefore, it should be further determined whether in human cells all two-neighboring-gene RNAs are derived from transcriptional readthrough and whether trans-splicing truly exists.

Organellar genome assembly methods and comparative analysis of horticultural plants

Although organellar genomes (including chloroplast and mitochondrial genomes) are smaller than nuclear genomes in size and gene number, organellar genomes are very important for the investigation of plant evolution and molecular ecology mechanisms. Few studies have focused on the organellar genomes of horticultural plants. Approximately 1193 chloroplast genomes and 199 mitochondrial genomes of land plants are available in the National Center for Biotechnology Information (NCBI), of which only 39 are from horticultural plants. In this paper, we report an innovative and efficient method for high-quality horticultural organellar genome assembly from next-generation sequencing (NGS) data. Sequencing reads were first assembled by Newbler, Amos, and Minimus software with default parameters. The remaining gaps were then filled through BLASTN search and PCR. The complete DNA sequence was corrected based on Illumina sequencing data using BWA (Burrows-Wheeler Alignment tool) software. The advantage of this approach is that there is no need to isolate organellar DNA from total DNA during sample preparation. Using this procedure, the complete mitochondrial and chloroplast genomes of an ornamental plant, Salix suchowensis, and a fruit tree, Ziziphus jujuba, were identified. This study shows that horticultural plants have similar mitochondrial and chloroplast sequence organization to other seed plants. Most horticultural plants demonstrate a slight bias toward A+T rich features in the mitochondrial genome. In addition, a phylogenetic analysis of 39 horticultural plants based on 15 protein-coding genes showed that some mitochondrial genes are horizontally transferred from chloroplast DNA. Our study will provide an important reference for organellar genome assembly in other horticultural plants. Furthermore, phylogenetic analysis of the organellar genomes of horticultural plants could accurately clarify the unanticipated relationships among these plants.

Ultra-Deep Sequencing Analysis on HIV Drug-Resistance-Associated Mutations Among HIV-Infected Individuals: First Report from the Philippines

A sharp increase in the number of people living with HIV has been documented in the Philippines. In response, the government has instituted antiretroviral therapy (ART) nationwide through HIV treatment hubs. However, no data presently exist on the status of ART drug-resistance-associated mutations (DRMs). In this study, we aim at analyzing DRM profiles in the Philippines and at providing comprehensive data on DRMs to guide treatment decisions and prevent viral failures. We conducted a cross-sectional study in 119 volunteers who tested positive for HIV from more than 8,000 participants screened for HIV across the nation through the 2013 Integrated HIV Behavioral and Serologic Surveillance (IHBSS) program. Amplicons were generated from plasma RNA by using primers designed to analyze diverse HIV-1 isolates targeting the reverse transcriptase region and sequenced on a 454 ultra-deep sequencing (UDS) platform to assess DRMs. DRMs were defined by using the Stanford HIV drug resistance database, and we found only 2 from 110 evaluable individuals with major HIV variants (>20% prevalence) that were highly resistant to the non-nucleoside reverse transcriptase inhibitor (NNRTI: efavirenz and nevirapine). However, a larger fraction of individuals harbored minority drug-resistant HIV variants (0.5%-20% prevalence) and they were highly resistant to NNRTI nevirapine (89/110), rilpivirine (5/110), and efavirenz (49/110). This study is the first report on the presence of HIV drug resistance in the Philippines and demonstrates the utility of UDS in assisting the detection of HIV minor variants. Monitoring for ART-DRMs will assist in improving HIV management strategies in curtailing the evolving epidemic in the Philippines.

Deep sequencing in the management of hepatitis virus infections

The hepatitis viruses represent a major public health problem worldwide. Procedures for characterization of the genomic composition of their populations, accurate diagnosis, identification of multiple infections, and information on inhibitor-escape mutants for treatment decisions are needed. Deep sequencing methodologies are extremely useful for these viruses since they replicate as complex and dynamic quasispecies swarms whose complexity and mutant composition are biologically relevant traits. Population complexity is a major challenge for disease prevention and control, but also an opportunity to distinguish among related but phenotypically distinct variants that might anticipate disease progression and treatment outcome. Detailed characterization of mutant spectra should permit choosing better treatment options, given the increasing number of new antiviral inhibitors available. In the present review we briefly summarize our experience on the use of deep sequencing for the management of hepatitis virus infections, particularly for hepatitis B and C viruses, and outline some possible new applications of deep sequencing for these important human pathogens.

A high-throughput assay for quantitative measurement of PCR errors

The accuracy with which DNA polymerase can replicate a template DNA sequence is an extremely important property that can vary by an order of magnitude from one enzyme to another. The rate of nucleotide misincorporation is shaped by multiple factors, including PCR conditions and proofreading capabilities, and proper assessment of polymerase error rate is essential for a wide range of sensitive PCR-based assays. In this paper, we describe a method for studying polymerase errors with exceptional resolution, which combines unique molecular identifier tagging and high-throughput sequencing. Our protocol is less laborious than commonly-used methods, and is also scalable, robust and accurate. In a series of nine PCR assays, we have measured a range of polymerase accuracies that is in line with previous observations. However, we were also able to comprehensively describe individual errors introduced by each polymerase after either 20 PCR cycles or a linear amplification, revealing specific substitution preferences and the diversity of PCR error frequency profiles. We also demonstrate that the detected high-frequency PCR errors are highly recurrent and that the position in the template sequence and polymerase-specific substitution preferences are among the major factors influencing the observed PCR error rate.

MAGERI: Computational pipeline for molecular-barcoded targeted resequencing

Unique molecular identifiers (UMIs) show outstanding performance in targeted high-throughput resequencing, being the most promising approach for the accurate identification of rare variants in complex DNA samples. This approach has application in multiple areas, including cancer diagnostics, thus demanding dedicated software and algorithms. Here we introduce MAGERI, a computational pipeline that efficiently handles all caveats of UMI-based analysis to obtain high-fidelity mutation profiles and call ultra-rare variants. Using an extensive set of benchmark datasets including gold-standard biological samples with known variant frequencies, cell-free DNA from tumor patient blood samples and publicly available UMI-encoded datasets we demonstrate that our method is both robust and efficient in calling rare variants. The versatility of our software is supported by accurate results obtained for both tumor DNA and viral RNA samples in datasets prepared using three different UMI-based protocols.

Aligning coding sequences with frameshift extension penalties

BACKGROUND

Frameshift translation is an important phenomenon that contributes to the appearance of novel coding DNA sequences (CDS) and functions in gene evolution, by allowing alternative amino acid translations of gene coding regions. Frameshift translations can be identified by aligning two CDS, from a same gene or from homologous genes, while accounting for their codon structure. Two main classes of algorithms have been proposed to solve the problem of aligning CDS, either by amino acid sequence alignment back-translation, or by simultaneously accounting for the nucleotide and amino acid levels. The former does not allow to account for frameshift translations and up to now, the latter exclusively accounts for frameshift translation initiation, not considering the length of the translation disruption caused by a frameshift.

RESULTS

We introduce a new scoring scheme with an algorithm for the pairwise alignment of CDS accounting for frameshift translation initiation and length, while simultaneously considering nucleotide and amino acid sequences. The main specificity of the scoring scheme is the introduction of a penalty cost accounting for frameshift extension length to compute an adequate similarity score for a CDS alignment. The second specificity of the model is that the search space of the problem solved is the set of all feasible alignments between two CDS. Previous approaches have considered restricted search space or additional constraints on the decomposition of an alignment into length-3 sub-alignments. The algorithm described in this paper has the same asymptotic time complexity as the classical Needleman-Wunsch algorithm.

CONCLUSIONS

We compare the method to other CDS alignment methods based on an application to the comparison of pairs of CDS from homologous human, mouse and cow genes of ten mammalian gene families from the Ensembl-Compara database. The results show that our method is particularly robust to parameter changes as compared to existing methods. It also appears to be a good compromise, performing well both in the presence and absence of frameshift translations. An implementation of the method is available at https://github.com/UdeS-CoBIUS/FsePSA.

An Improved Barcoded Oligonucleotide Primers-based Next-generation Sequencing Approach for Direct Identification of Viral Pathogens in Clinical Specimens

OBJECTIVE

To provide a feasible and cost-effective next-generation sequencing (NGS) method for accurate identification of viral pathogens in clinical specimens, because enormous limitations impede the clinical use of common NGS, such as high cost, complicated procedures, tremendous data analysis, and high background noise in clinical samples.

METHODS

Viruses from cell culture materials or clinical specimens were identified following an improved NGS procedure: reduction of background noise by sample preprocessing, viral enrichment by barcoded oligonucleotide (random hexamer or non-ribosomal hexanucleotide) primer-based amplification, fragmentation-free library construction and sequencing of one-tube mixtures, as well as rapid data analysis using an in-house pipeline.

RESULTS

NGS data demonstrated that both barcoded primer sets were useful to simultaneously capture multiple viral pathogens in cell culture materials or clinical specimens and verified that hexanucleotide primers captured as many viral sequences as hexamers did. Moreover, direct testing of clinical specimens using this improved hexanucleotide primer-based NGS approach provided further detailed genotypes of enteroviruses causing hand, foot, and mouth disease (HFMD) and identified other potential viruses or differentiated misdiagnosis events.

CONCLUSION

The improved barcoded oligonucleotide primer-based NGS approach is simplified, time saving, cost effective, and appropriate for direct identification of viral pathogens in clinical practice.

Error rates, PCR recombination, and sampling depth in HIV-1 whole genome deep sequencing

Deep sequencing is a powerful and cost-effective tool to characterize the genetic diversity and evolution of virus populations. While modern sequencing instruments readily cover viral genomes many thousand fold and very rare variants can in principle be detected, sequencing errors, amplification biases, and other artifacts can limit sensitivity and complicate data interpretation. For this reason, the number of studies using whole genome deep sequencing to characterize viral quasi-species in clinical samples is still limited. We have previously undertaken a large scale whole genome deep sequencing study of HIV-1 populations. Here we discuss the challenges, error profiles, control experiments, and computational test we developed to quantify the accuracy of variant frequency estimation.

Ultrasensitive single-genome sequencing: accurate, targeted, next generation sequencing of HIV-1 RNA

BACKGROUND

Although next generation sequencing (NGS) offers the potential for studying virus populations in unprecedented depth, PCR error, amplification bias and recombination during library construction have limited its use to population sequencing and measurements of unlinked allele frequencies. Here we report a method, termed ultrasensitive Single-Genome Sequencing (uSGS), for NGS library construction and analysis that eliminates PCR errors and recombinants, and generates single-genome sequences of the same quality as the "gold-standard" of HIV-1 single-genome sequencing assay but with more than 100-fold greater depth.

RESULTS

Primer ID tagged cDNA was synthesized from mixtures of cloned BH10 wild-type and mutant HIV-1 transcripts containing ten drug resistance mutations. First, the resultant cDNA was divided and NGS libraries were generated in parallel using two methods: uSGS and a method applying long PCR primers to attach the NGS adaptors (LP-PCR-1). Second, cDNA was divided and NGS libraries were generated in parallel comparing 3 methods: uSGS and 2 methods adapted from more recent reports using variations of the long PCR primers to attach the adaptors (LP-PCR-2 and LP-PCR-3). Consistently, the uSGS method amplified a greater proportion of cDNAs, averaging 30% compared to 13% for LP-PCR-1, 21% for LP-PCR-2 and 14% for LP-PCR-3. Most importantly, when the uSGS sequences were binned according to their primer IDs, 94% of the bins did not contain PCR recombinant sequences versus only 55, 75 and 65% for LP-PCR-1, 2 and 3, respectively. Finally, when uSGS was applied to plasma samples from HIV-1 infected donors, both frequent and rare variants were detected in each sample and neighbor-joining trees revealed clusters of genomes driven by the linkage of these mutations, showing the lack of PCR recombinants in the datasets.

CONCLUSIONS

The uSGS assay can be used for accurate detection of rare variants and for identifying linkage of rare alleles associated with HIV-1 drug resistance. In addition, the method allows accurate in-depth analyses of the complex genetic relationships of viral populations in vivo.

Promises and pitfalls of Illumina sequencing for HIV resistance genotyping

Genetic sequencing ("genotyping") plays a critical role in the modern clinical management of HIV infection. This virus evolves rapidly within patients because of its error-prone reverse transcriptase and short generation time. Consequently, HIV variants with mutations that confer resistance to one or more antiretroviral drugs can emerge during sub-optimal treatment. There are now multiple HIV drug resistance interpretation algorithms that take the region of the HIV genome encoding the major drug targets as inputs; expert use of these algorithms can significantly improve to clinical outcomes in HIV treatment. Next-generation sequencing has the potential to revolutionize HIV resistance genotyping by lowering the threshold that rare but clinically significant HIV variants can be detected reproducibly, and by conferring improved cost-effectiveness in high-throughput scenarios. In this review, we discuss the relative merits and challenges of deploying the Illumina MiSeq instrument for clinical HIV genotyping.

An Efficient Microarray-Based Genotyping Platform for the Identification of Drug-Resistance Mutations in Majority and Minority Subpopulations of HIV-1 Quasispecies

The response of human immunodeficiency virus type 1 (HIV-1) quasispecies to antiretroviral therapy is influenced by the ensemble of mutants that composes the evolving population. Low-abundance subpopulations within HIV-1 quasispecies may determine the viral response to the administered drug combinations. However, routine sequencing assays available to clinical laboratories do not recognize HIV-1 minority variants representing less than 25% of the population. Although several alternative and more sensitive genotyping techniques have been developed, including next-generation sequencing (NGS) methods, they are usually very time consuming, expensive and require highly trained personnel, thus becoming unrealistic approaches in daily clinical practice. Here we describe the development and testing of a HIV-1 genotyping DNA microarray that detects and quantifies, in majority and minority viral subpopulations, relevant mutations and amino acid insertions in 42 codons of the pol gene associated with drug- and multidrug-resistance to protease (PR) and reverse transcriptase (RT) inhibitors. A customized bioinformatics protocol has been implemented to analyze the microarray hybridization data by including a new normalization procedure and a stepwise filtering algorithm, which resulted in the highly accurate (96.33%) detection of positive/negative signals. This microarray has been tested with 57 subtype B HIV-1 clinical samples extracted from multi-treated patients, showing an overall identification of 95.53% and 89.24% of the queried PR and RT codons, respectively, and enough sensitivity to detect minority subpopulations representing as low as 5–10% of the total quasispecies. The developed genotyping platform represents an efficient diagnostic and prognostic tool useful to personalize antiviral treatments in clinical practice.

HIV-1 Protease, Reverse Transcriptase, and Integrase Variation

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.

Detecting primary drug-resistant mutations in Korean HIV patients using ultradeep pyrosequencing

HIV primary resistance, drug resistance in treatment-naïve patients, is an emerging public health issue. The prevalence of HIV primary resistance mutations down to the level of 1% minor variants was investigated using ultradeep pyrosequencing (UDPS) in HIV-positive Korean blood donors and in treatment naïve chronic patients for the comparison. The entire pol region was sequenced from 25 HIV-positive blood donors, and 18 treatment-naïve chronic HIV patients. UDPS was successful in 19 blood donors and 18 chronic patients. In total, 1,011,338 sequence reads were aligned, and 28,093 sequence reads were aligned on average per sample. The prevalence of HIV primary resistance mutations in the HIV-positive blood donors and chronic HIV patients were 63.2% and 44.4% according to UDPS, respectively. Protease inhibitor (PI) drugs demonstrated different patterns in HIV-positive blood donors and chronic HIV patients, whereas non-nucleoside reverse transcriptase inhibitors (NNRTI), nucleoside reverse transcriptase inhibitors (NRTI), and integrase inhibitor (INI) drugs showed similar patterns between the two groups. Higher level of primary resistance prevalence was observed mainly because UDPS method could detect mutations in minor variants with 1-10% frequency. The higher resistance prevalence was observed in HIV-positive blood donors than in chronic patients. Considering that treatments for HIV-infected patients were recently amended to start at an earlier stage, information about degree of drug resistance to each drug between the two groups would help to establish future policies, design additional clinical trials, assess HIV patient care in Korea.

New bioinformatic tool for quick identification of functionally relevant endogenous retroviral inserts in human genome

Endogenous retroviruses (ERVs) and LTR retrotransposons (LRs) occupy ∼8% of human genome. Deep sequencing technologies provide clues to understanding of functional relevance of individual ERVs/LRs by enabling direct identification of transcription factor binding sites (TFBS) and other landmarks of functional genomic elements. Here, we performed the genome-wide identification of human ERVs/LRs containing TFBS according to the ENCODE project. We created the first interactive ERV/LRs database that groups the individual inserts according to their familial nomenclature, number of mapped TFBS and divergence from their consensus sequence. Information on any particular element can be easily extracted by the user. We also created a genome browser tool, which enables quick mapping of any ERV/LR insert according to genomic coordinates, known human genes and TFBS. These tools can be used to easily explore functionally relevant individual ERV/LRs, and for studying their impact on the regulation of human genes. Overall, we identified ∼110,000 ERV/LR genomic elements having TFBS. We propose a hypothesis of “domestication” of ERV/LR TFBS by the genome milieu including subsequent stages of initial epigenetic repression, partial functional release, and further mutation-driven reshaping of TFBS in tight coevolution with the enclosing genomic loci.

Quantifying Next Generation Sequencing Sample Pre-Processing Bias in HIV-1 Complete Genome Sequencing

Genetic analyses play a central role in infectious disease research. Massively parallelized “mechanical cloning” and sequencing technologies were quickly adopted by HIV researchers in order to broaden the understanding of the clinical importance of minor drug-resistant variants. These efforts have, however, remained largely limited to small genomic regions. The growing need to monitor multiple genome regions for drug resistance testing, as well as the obvious benefit for studying evolutionary and epidemic processes makes complete genome sequencing an important goal in viral research. In addition, a major drawback for NGS applications to RNA viruses is the need for large quantities of input DNA. Here, we use a generic overlapping amplicon-based near full-genome amplification protocol to compare low-input enzymatic fragmentation (Nextera™) with conventional mechanical shearing for Roche 454 sequencing. We find that the fragmentation method has only a modest impact on the characterization of the population composition and that for reliable results, the variation introduced at all steps of the procedure—from nucleic acid extraction to sequencing—should be taken into account, a finding that is also relevant for NGS technologies that are now more commonly used. Furthermore, by applying our protocol to deep sequence a number of pre-therapy plasma and PBMC samples, we illustrate the potential benefits of a near complete genome sequencing approach in routine genotyping.

Evaluation of Automatic Analysis of Ultradeep Pyrosequencing Raw Data to Determine Percentages of HIV Resistance Mutations in Patients Followed-Up in Hospital

A major obstacle to using next generation sequencing (NGS) technology in clinical routine practice is reliable data analysis. Thousands of sequences need to be aligned and validated, to exclude sequencing artifacts and generate accurate results. We compared two analysis pipelines for Roche 454 ultradeep pyrosequencing (UDPS) raw data generated from HIV-1 clinical samples: a commercial and fully automated Web-based software NGS HIV-1 Module (SmartGene, Zug, Switzerland) vs. the Amplicon Variant Analyzer software (AVA, 454 Life Sciences; Roche). Results were also compared to those obtained with Sanger sequencing. HIV-1 reverse transcriptase and protease genes from 34 plasma samples were submitted to Sanger sequencing and GS Junior UDPS. Raw UDPS data (sff files) from all samples were analyzed with AVA 2.7 software plus manual review of the alignments and the fully automated SmartGene NGS HIV-1 Module prototype (SMG). Results obtained with both analysis pipelines showed good correlation (85.0%). Divergent results were mainly observed at homopolymer positions, such as K101, where the frame-aware alignment and error corrections of the automated approach were more efficient and more accurate, both in terms of detecting and quantifying drug resistance mutations. Our study shows that NGS data can easily be analyzed via a fully automated analysis pipeline, here the SmartGene NGS HIV-1 Module, thus minimizing the need for manual review of alignments by the user, otherwise essential to ensure accurate results. Such automated analysis pipelines may facilitate the adoption of NGS platforms in the routine clinical laboratory.

Rapid identification of pathogens using molecular techniques

Real-time PCR is the traditional face of nucleic acid detection in the diagnostic microbiology laboratory and is now generally regarded as robust enough to be widely adopted. Methods based on nucleic acid detection of this type are bringing increased accuracy to diagnosis in areas where culture is difficult and/or expensive, and these methods are often effective partners to other rapid molecular diagnostic tools such as matrix-assisted laser desorption ionisation-time of flight mass spectrometry (MALDI-TOF MS). This change in practice has particularly affected the recognition of viruses and fastidious or antibiotic-exposed bacteria, but has been also shown to be effective in the recognition of troublesome or specialised phenotypes such as antiviral resistance and transmissible antibiotic resistance in the Enterobacteriaceae. Quantitation and high-intensity sequencing (of multiple whole genomes) has brought new opportunities as well as new challenges to the microbiology community. Diagnostic microbiologists currently training might be expected to deal less with the culture-based techniques of the last half-century than with the high-volume data and complex analyses of the next.

Comparison of 454 Ultra-Deep Sequencing and Allele-Specific Real-Time PCR with Regard to the Detection of Emerging Drug-Resistant Minor HIV-1 Variants after Antiretroviral Prophylaxis for Vertical Transmission

BACKGROUND

Pregnant HIV-infected women were screened for the development of HIV-1 drug resistance after implementation of a triple-antiretroviral transmission prophylaxis as recommended by the WHO in 2006. The study offered the opportunity to compare amplicon-based 454 ultra-deep sequencing (UDS) and allele-specific real-time PCR (ASPCR) for the detection of drug-resistant minor variants in the HIV-1 reverse transcriptase (RT).

METHODS

Plasma samples from 34 Tanzanian women were previously analysed by ASPCR for key resistance mutations in the viral RT selected by AZT, 3TC, and NVP (K70R, K103N, Y181C, M184V, T215Y/F). In this study, the RT region of the same samples was investigated by amplicon-based UDS for resistance mutations using the 454 GS FLX System.

RESULTS

Drug-resistant HIV-variants were identified in 69% (20/29) of women by UDS and in 45% (13/29) by ASPCR. The absolute number of resistance mutations identified by UDS was twice that identified by ASPCR (45 vs 24). By UDS 14 of 24 ASPCR-detected resistance mutations were identified at the same position. The overall concordance between UDS and ASPCR was 61.0% (25/41). The proportions of variants quantified by UDS were approximately 2-3 times lower than by ASPCR. Amplicon generation from samples with viral loads below 20,000 copies/ml failed more frequently by UDS compared to ASPCR (limit of detection = 650 copies/ml), resulting in missing or insufficient sequence coverage.

CONCLUSIONS

Both methods can provide useful information about drug-resistant minor HIV-1 variants. ASPCR has a higher sensitivity than UDS, but is restricted to single resistance mutations. In contrast, UDS is limited by its requirement for high viral loads to achieve sufficient sequence coverage, but the sequence information reveals the complete resistance patterns within the genomic region analysed. Improvements to the UDS limit of detection are in progress, and UDS could then facilitate monitoring of drug-resistant minor variants in the HIV-1 quasispecies.

Genome plasticity of triple-reassortant H1N1 influenza A virus during infection of vaccinated pigs

To gain insight into the evolution of influenza A viruses (IAVs) during infection of vaccinated pigs, we experimentally infected a 3-week-old naive pig with a triple-reassortant H1N1 IAV and placed the seeder pig in direct contact with a group of age-matched vaccinated pigs (n = 10). We indexed the genetic diversity and evolution of the virus at an intra-host level by deep sequencing the entire genome directly from nasal swabs collected at two separate samplings during infection. We obtained 13 IAV metagenomes from 13 samples, which included the virus inoculum and two samples from each of the six pigs that tested positive for IAV during the study. The infection produced a population of heterogeneous alleles (sequence variants) that was dynamic over time. Overall, 794 polymorphisms were identified amongst all samples, which yielded 327 alleles, 214 of which were unique sequences. A total of 43 distinct haemagglutinin proteins were translated, two of which were observed in multiple pigs, whereas the neuraminidase (NA) was conserved and only one dominant NA was found throughout the study. The genetic diversity of IAVs changed dynamically within and between pigs. However, most of the substitutions observed in the internal gene segments were synonymous. Our results demonstrated remarkable IAV diversity, and the complex, rapid and dynamic evolution of IAV during infection of vaccinated pigs that can only be appreciated with repeated sampling of individual animals and deep sequence analysis.

Artifact-Free Quantification and Sequencing of Rare Recombinant Viruses by Using Drop-Based Microfluidics

Recombination is an important driver in the evolution of viruses and thus is key to understanding viral epidemics and improving strategies to prevent future outbreaks. Characterization of rare recombinant subpopulations remains technically challenging because of artifacts such as artificial recombinants, known as chimeras, and amplification bias. To overcome this, we have developed a high-throughput microfluidic technique with a second verification step in order to amplify and sequence single recombinant viruses with high fidelity in picoliter drops. We obtained the first artifact-free estimate of in vitro recombination rate between murine norovirus strains MNV-1 and WU20 co-infecting a cell (P(rec) = 3.3 × 10(-4) ± 2 × 10(-5) ) for a 1205 nt region. Our approach represents a time- and cost-effective improvement over current methods, and can be adapted for genomic studies requiring artifact- and bias-free selective amplification, such as microbial pathogens, or rare cancer cells.

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.

Retrohoming of a Mobile Group II Intron in Human Cells Suggests How Eukaryotes Limit Group II Intron Proliferation

Mobile bacterial group II introns are evolutionary ancestors of spliceosomal introns and retroelements in eukaryotes. They consist of an autocatalytic intron RNA (a “ribozyme”) and an intron-encoded reverse transcriptase, which function together to promote intron integration into new DNA sites by a mechanism termed “retrohoming”. Although mobile group II introns splice and retrohome efficiently in bacteria, all examined thus far function inefficiently in eukaryotes, where their ribozyme activity is limited by low Mg²⁺ concentrations, and intron-containing transcripts are subject to nonsense-mediated decay (NMD) and translational repression. Here, by using RNA polymerase II to express a humanized group II intron reverse transcriptase and T7 RNA polymerase to express intron transcripts resistant to NMD, we find that simply supplementing culture medium with Mg²⁺ induces the Lactococcus lactis Ll.LtrB intron to retrohome into plasmid and chromosomal sites, the latter at frequencies up to ~0.1%, in viable HEK-293 cells. Surprisingly, under these conditions, the Ll.LtrB intron reverse transcriptase is required for retrohoming but not for RNA splicing as in bacteria. By using a genetic assay for in vivo selections combined with deep sequencing, we identified intron RNA mutations that enhance retrohoming in human cells, but <4-fold and not without added Mg²⁺. Further, the selected mutations lie outside the ribozyme catalytic core, which appears not readily modified to function efficiently at low Mg²⁺ concentrations. Our results reveal differences between group II intron retrohoming in human cells and bacteria and suggest constraints on critical nucleotide residues of the ribozyme core that limit how much group II intron retrohoming in eukaryotes can be enhanced. These findings have implications for group II intron use for gene targeting in eukaryotes and suggest how differences in intracellular Mg²⁺ concentrations between bacteria and eukarya may have impacted the evolution of introns and gene expression mechanisms.

Mobile group II introns are bacterial retrotransposons that are evolutionary ancestors of spliceosomal introns and retroelements in eukaryotes. They consist of an autocatalytic intron RNA (a ribozyme) and an intron-encoded reverse transcriptase, which together promote intron mobility to new DNA sites by a mechanism called retrohoming. Although found in bacteria, archaea and eukaryotic organelles, group II introns are absent from eukaryotic nuclear genomes, where host defenses impede their expression and lower intracellular Mg²⁺ concentrations limit their ribozyme activity. Here, we developed a mobile group II intron expression system that bypasses expression barriers and show that simply adding Mg²⁺ to culture medium enables group II intron retrohoming into plasmid and chromosomal target sites in human cells at appreciable frequencies. Genetic selections and deep sequencing identified intron RNA mutations that moderately enhance retrohoming in human cells, but not without added Mg²⁺. Thus, low Mg²⁺ concentrations in human cells are a natural barrier to efficient retrohoming that is not readily overcome by mutational variation and selection. Our results have implications for group II intron use for gene targeting in higher organisms and highlight the impact of different intracellular environments on intron evolution and gene expression mechanisms in bacteria and eukarya.

Similarities between Human Immunodeficiency Virus Type 1 and Hepatitis C Virus Genetic and Phenotypic Protease Quasispecies Diversity

Human immunodeficiency virus type 1 (HIV-1) and hepatitis C virus (HCV) are two highly variable RNA viruses that cause chronic infections in humans. Although HCV likely preceded the AIDS epidemic by some decades, the global spread of both viruses is a relatively recent event. Nevertheless, HCV global diversity is higher than that of HIV-1. To identify differences in mutant diversity, we compared the HIV-1 protease and HCV NS3 protease quasispecies. Three protease gene quasispecies samples per virus, isolated from a total of six infected patients, were genetically and phenotypically analyzed at high resolution (HIV-1, 308 individual clones; HCV, 299 clones). Single-nucleotide variant frequency did not differ between quasispecies from the two viruses (HIV-1, 2.4 × 10(-3) ± 0.4 × 10(-3); HCV, 2.1 × 10(-3) ± 0.5 × 10(-3)) (P = 0.1680). The proportion of synonymous substitutions to potential synonymous sites was similar (3.667 ± 0.6667 and 2.183 ± 0.9048, respectively) (P = 0.2573), and Shannon's entropy values did not differ between HIV-1 and HCV (0.84 ± 0.02 and 0.83 ± 0.12, respectively) (P = 0.9408). Of note, 65% (HIV-1) and 67% (HCV) of the analyzed enzymes displayed detectable protease activity, suggesting that both proteases have a similar mutational robustness. In both viruses, there was a rugged protease enzymatic activity landscape characterized by a sharp peak, representing the master sequence, surrounded by a collection of diverse variants present at lower frequencies. These results indicate that nucleotide quasispecies diversification during chronic infection is not responsible for the higher worldwide genetic diversity observed in HCV.

IMPORTANCE

HCV global diversity is higher than that of HIV-1. We asked whether HCV genetic diversification during infection is responsible for the higher worldwide genetic diversity observed in HCV. To this end, we analyzed and compared the genotype and enzymatic activities of HIV-1 and HCV protease quasispecies existing in infected individuals. Our results indicate that HIV-1 and HCV protease quasispecies have very similar genetic diversity and comparable rugged enzymatic activity landscapes. Therapy for HCV has expanded, with new therapeutic agents such as the direct-acting antivirals (DAAs). DAAs, which target HCV NS3 protease and other virus proteins, have improved cure rates. However, major questions remain to be elucidated regarding the virologic correlates of HCV eradication. The findings shown here may help our understanding of the different therapeutic responses observed during chronic HCV infection.