A comparison of parallel pyrosequencing and sanger clone-based sequencing and its impact on the characterization of the genetic diversity of HIV-1

Genetic Diversity of Grapevine Virus A in Three Australian Vineyards Using Amplicon High Throughput Sequencing (Amplicon-HTS)

Shiraz disease (SD) is one of the most destructive viral diseases of grapevines in Australia and is known to cause significant economic loss to local growers. Grapevine virus A (GVA) was reported to be the key pathogen associated with this disease. This study aimed to better understand the diversity of GVA variants both within and between individual SD and grapevine leafroll disease (LRD) affected grapevines located at vineyards in South Australia. Amplicon high throughput sequencing (Amplicon-HTS) combined with median-joining networks (MJNs) was used to analyze the variability in specific gene regions of GVA variants. Several GVAII variant groups contain samples from both vineyards studied, suggesting that these GVAII variants were from a common origin. Variant groups analyzed by MJNs using the overall data set denote that there may be a possible relationship between variant groups of GVA and the geographical location of the grapevines.

A sensitive and inexpensive high-resolution melting-based testing algorithm for diagnosis of transient abnormal myelopoiesis and myeloid leukemia of Down syndrome

Patients with Down syndrome (DS) are commonly affected by a pre-leukemic disorder known as transient abnormal myelopoiesis (TAM). This condition usually undergoes spontaneous remission within the first 2 months after birth; however, in children under 5, 20%-30% of cases evolve to myeloid leukemia of Down syndrome (ML-DS). TAM and ML-DS are caused by co-operation between trisomy 21 and acquired mutations in the GATA1 gene. Currently, only next-generation sequencing (NGS)-based methodologies are sufficiently sensitive for diagnosis in samples with small GATA1 mutant clones (≤10% blasts). Alternatively, this study presents research on a new, fast, sensitive, and inexpensive high-resolution melting (HRM)-based diagnostic approach that allows the detection of most cases of GATA1 mutations, including silent TAM. The algorithm first uses flow cytometry for blast count, followed by HRM and Sanger sequencing to search for mutations on exons 2 and 3 of GATA1. We analyzed 138 samples of DS patients: 110 of asymptomatic neonates, 10 suspected of having TAM, and 18 suspected of having ML-DS. Our algorithm enabled the identification of 33 mutant samples, among them five cases of silent TAM (5/110) and seven cases of ML-DS (7/18) with blast count ≤10%, in which GATA1 alterations were easily detected by HRM. Depending on the type of genetic variation and its location, our methodology reached sensitivity similar to that obtained by NGS (0.3%) at a considerably reduced time and cost, thus making it accessible worldwide.

CD8+ T-Cell Exhaustion Phenotype in Chronic Hepatitis C Virus Infection Is Associated With Epitope Sequence Variation

Background and Aims

During chronic hepatitis C virus (HCV) infection, CD8⁺ T-cells become functionally exhausted, undergoing progressive phenotypic changes, i.e., overexpression of “inhibitory” molecules such as PD-1 (programmed cell death protein 1) and/or Tim-3 (T-cell immunoglobulin and mucin domain-containing molecule-3). The extreme intrahost genetic diversity of HCV is a major mechanism of immune system evasion, facilitating epitope escape. The aim of the present study was to determine whether T-cell exhaustion phenotype in chronic HCV infection is related to the sequence repertoire of NS3 viral immunodominant epitopes.

Methods

The study population was ninety prospective patients with chronic HCV genotype 1b infection. Populations of peripheral blood CD8⁺ T-cells expressing PD-1/Tim-3 were assessed by multiparametric flow cytometry, including HCV-specific T-cells after magnetic-based enrichment using MHC-pentamer. Autologous epitope sequences were inferred from next-generation sequencing. The correction of sequencing errors and genetic variants reconstruction was performed using Quasirecomb.

Results

There was an interplay between the analyzed epitopes sequences and exhaustion phenotype of CD8⁺ T-cells. A predominance of NS3₁₄₀₆ epitope sequence, representing neither prototype KLSGLGLNAV nor cross-reactive variants (KLSSLGLNAV, KLSGLGINAV or KLSALGLNAV), was associated with higher percentage of HCV-specific CD8⁺PD-1⁺Tim-3⁺ T-cells, P=0.0102. Variability (at least two variants) of NS3₁₄₀₆ epitope sequence was associated with increased frequencies of global CD8⁺PD-1⁺Tim-3⁺ T-cells (P=0.0197) and lower frequencies of CD8⁺PD-1⁻Tim-3⁻ T-cells (P=0.0079). In contrast, infection with NS3₁₀₇₃ dominant variant epitope (other than prototype CVNGVCWTV) was associated with lower frequency of global CD8⁺PD-1⁺Tim-3⁺ T-cells (P=0.0054).

Conclusions

Our results indicate that PD-1/Tim-3 receptor expression is largely determined by viral epitope sequence and is evident for both HCV-specific and global CD8⁺ T-cells, pointing to the importance of evaluating autologous viral epitope sequences in the investigation of CD8⁺ T-cell exhaustion in HCV infection.

Culture-independent studies on bacterial dysbiosis in oral and oropharyngeal squamous cell carcinoma: A systematic review

Imbalance within the resident bacterial community (dysbiosis), rather than the presence and activity of a single organism, has been proposed to be associated with, and to influence, the development and progression of various diseases; however, the existence and significance of dysbiosis in oral/oropharyngeal cancer is yet to be clearly established. A systematic search (conducted on 25/01/2018 and updated on 25/05/2018) was performed on three databases (Pubmed, Web of Science & Scopus) to identify studies employing culture-independent methods which investigated the bacterial community in oral/oropharyngeal cancer patients compared to control subjects. Of the 1546 texts screened, only fifteen publications met the pre-determined selection criteria. Data extracted from 731 cases and 809 controls overall, could not identify consistent enrichment of any particular taxon in oral/oropharyngeal cancers, although common taxa could be identified between studies. Six studies reported the enrichment of Fusobacteria in cancer at different taxonomic levels whereas four studies reported an increase in Parvimonas. Changes in microbial diversity remained inconclusive, with four studies showing a higher diversity in controls, three studies showing a higher diversity in tumors and three additional studies showing no difference between tumors and controls. Even though most studies identified a component of dysbiosis in oral/oropharyngeal cancer, methodological and analytical variations prevented a standardized summary, which highlights the necessity for studies of superior quality and magnitude employing standardized methodology and reporting. Indeed an holistic metagenomic approach is likely to be more meaningful, as is understanding of the overall metabolome, rather than a mere enumeration of the organisms present.

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.

Added Value of Next-Generation Sequencing for Multilocus Sequence Typing Analysis of a Pneumocystis jirovecii Pneumonia Outbreak1

Pneumocystis jirovecii is a major threat for immunocompromised patients, and clusters of pneumocystis pneumonia (PCP) have been increasingly described in transplant units during the past decade. Exploring an outbreak transmission network requires complementary spatiotemporal and strain-typing approaches. We analyzed a PCP outbreak and demonstrated the added value of next-generation sequencing (NGS) for the multilocus sequence typing (MLST) study of P. jirovecii strains. Thirty-two PCP patients were included. Among the 12 solid organ transplant patients, 5 shared a major and unique genotype that was also found as a minor strain in a sixth patient. A transmission map analysis strengthened the suspicion of nosocomial acquisition of this strain for the 6 patients. NGS-MLST enables accurate determination of subpopulation, which allowed excluding other patients from the transmission network. NGS-MLST genotyping approach was essential to deciphering this outbreak. This innovative approach brings new insights for future epidemiologic studies on this uncultivable opportunistic fungus.

Antibiotic Resistome: Improving Detection and Quantification Accuracy for Comparative Metagenomics

The unprecedented rise of life-threatening antibiotic resistance (AR), combined with the unparalleled advances in DNA sequencing of genomes and metagenomes, has pushed the need for in silico detection of the resistance potential of clinical and environmental metagenomic samples through the quantification of AR genes (i.e., genes conferring antibiotic resistance). Therefore, determining an optimal methodology to quantitatively and accurately assess AR genes in a given environment is pivotal. Here, we optimized and improved existing AR detection methodologies from metagenomic datasets to properly consider AR-generating mutations in antibiotic target genes. Through comparative metagenomic analysis of previously published AR gene abundance in three publicly available metagenomes, we illustrate how mutation-generated resistance genes are either falsely assigned or neglected, which alters the detection and quantitation of the antibiotic resistome. In addition, we inspected factors influencing the outcome of AR gene quantification using metagenome simulation experiments, and identified that genome size, AR gene length, total number of metagenomics reads and selected sequencing platforms had pronounced effects on the level of detected AR. In conclusion, our proposed improvements in the current methodologies for accurate AR detection and resistome assessment show reliable results when tested on real and simulated metagenomic datasets.

Multiplexed highly-accurate DNA sequencing of closely-related HIV-1 variants using continuous long reads from single molecule, real-time sequencing

Single Molecule, Real-Time (SMRT) Sequencing (Pacific Biosciences, Menlo Park, CA, USA) provides the longest continuous DNA sequencing reads currently available. However, the relatively high error rate in the raw read data requires novel analysis methods to deconvolute sequences derived from complex samples. Here, we present a workflow of novel computer algorithms able to reconstruct viral variant genomes present in mixtures with an accuracy of >QV50. This approach relies exclusively on Continuous Long Reads (CLR), which are the raw reads generated during SMRT Sequencing. We successfully implement this workflow for simultaneous sequencing of mixtures containing up to forty different >9 kb HIV-1 full genomes. This was achieved using a single SMRT Cell for each mixture and desktop computing power. This novel approach opens the possibility of solving complex sequencing tasks that currently lack a solution.

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.

Deep sequencing of protease inhibitor resistant HIV patient isolates reveals patterns of correlated mutations in Gag and protease

While the role of drug resistance mutations in HIV protease has been studied comprehensively, mutations in its substrate, Gag, have not been extensively cataloged. Using deep sequencing, we analyzed a unique collection of longitudinal viral samples from 93 patients who have been treated with therapies containing protease inhibitors (PIs). Due to the high sequence coverage within each sample, the frequencies of mutations at individual positions were calculated with high precision. We used this information to characterize the variability in the Gag polyprotein and its effects on PI-therapy outcomes. To examine covariation of mutations between two different sites using deep sequencing data, we developed an approach to estimate the tight bounds on the two-site bivariate probabilities in each viral sample, and the mutual information between pairs of positions based on all the bounds. Utilizing the new methodology we found that mutations in the matrix and p6 proteins contribute to continued therapy failure and have a major role in the network of strongly correlated mutations in the Gag polyprotein, as well as between Gag and protease. Although covariation is not direct evidence of structural propensities, we found the strongest correlations between residues on capsid and matrix of the same Gag protein were often due to structural proximity. This suggests that some of the strongest inter-protein Gag correlations are the result of structural proximity. Moreover, the strong covariation between residues in matrix and capsid at the N-terminus with p1 and p6 at the C-terminus is consistent with residue-residue contacts between these proteins at some point in the viral life cycle.

Understanding the structure of HIV proteins and the function of drug-resistant mutations of these proteins is critical for the development of effective HIV treatments. Selected gag mutations have been shown to provide compensatory functions for protease resistance mutations and may directly contribute to the development of drug resistance. To determine associations between protease inhibitor mutations and gag, we utilized deep sequencing of HIV gag and protease from a collection of viral isolates from patients treated with highly active retroviral protease inhibitors. Deep sequencing allows for accurate measurement of mutation frequencies at each position, allowing estimation, using a novel method we developed, of the covariation between any two residues on gag. Using this information, we characterize the variation within gag and protease and identify the most strongly correlated pairs of inter- and intra-protein residues. Our results suggest that matrix and p1/p6 mutations form the core of a network of strongly correlated gag mutations and contribute to recurrent treatment failure. Extracting gag residue covariation information from the deep sequencing of patient viral samples may provide insight into structural aspects of the Gag polyprotein as well new areas for small molecule targeting to disrupt Gag function.

A comparison of 454 sequencing and clonal sequencing for the characterization of hepatitis C virus NS3 variants

We compared 454 amplicon sequencing with clonal sequencing for the characterization of intra-host hepatitis C virus (HCV) NS3 variants. Clonal and 454 sequences were obtained from 12 patients enrolled in a clinical phase I study for telaprevir, an NS3-4a protease inhibitor. Thirty-nine datasets were used to compare the consensus sequence, average pairwise distance, normalized Shannon entropy, phylogenetic tree topology and the number and frequency of variants derived from both sequencing techniques. In general, a good concordance was observed between both techniques for the majority of datasets. Discordant results were observed for 5 out of 39 clonal and 454 datasets, which could be attributed to primer-related selective amplification used for clonal sequencing. Both 454 and clonal datasets consisted of a few major variants and a large number of low-frequency variants. Telaprevir resistance-associated variants were observed in low frequencies and were detected more often by 454. We conclude that performance of 454 and clonal sequencing is comparable for the characterization of intra-host virus populations. Not surprisingly, 454 is superior for the detection of low frequency resistance-associated variants. However, despite the greater coverage, 454 failed to detect some low frequency variants detected by clonal sequencing.

The Quantification of Representative Sequences pipeline for amplicon sequencing: case study on within-population ITS1 sequence variation in a microparasite infecting Daphnia

Next generation sequencing (NGS) platforms are replacing traditional molecular biology protocols like cloning and Sanger sequencing. However, accuracy of NGS platforms has rarely been measured when quantifying relative frequencies of genotypes or taxa within populations. Here we developed a new bioinformatic pipeline (QRS) that pools similar sequence variants and estimates their frequencies in NGS data sets from populations or communities. We tested whether the estimated frequency of representative sequences, generated by 454 amplicon sequencing, differs significantly from that obtained by Sanger sequencing of cloned PCR products. This was performed by analysing sequence variation of the highly variable first internal transcribed spacer (ITS1) of the ichthyosporean Caullerya mesnili, a microparasite of cladocerans of the genus Daphnia. This analysis also serves as a case example of the usage of this pipeline to study within-population variation. Additionally, a public Illumina data set was used to validate the pipeline on community-level data. Overall, there was a good correspondence in absolute frequencies of C. mesnili ITS1 sequences obtained from Sanger and 454 platforms. Furthermore, analyses of molecular variance (amova) revealed that population structure of C. mesnili differs across lakes and years independently of the sequencing platform. Our results support not only the usefulness of amplicon sequencing data for studies of within-population structure but also the successful application of the QRS pipeline on Illumina-generated data. The QRS pipeline is freely available together with its documentation under GNU Public Licence version 3 at http://code.google.com/p/quantification-representative-sequences.

Challenges with using primer IDs to improve accuracy of next generation sequencing

Next generation sequencing technologies, like ultra-deep pyrosequencing (UDPS), allows detailed investigation of complex populations, like RNA viruses, but its utility is limited by errors introduced during sample preparation and sequencing. By tagging each individual cDNA molecule with barcodes, referred to as Primer IDs, before PCR and sequencing these errors could theoretically be removed. Here we evaluated the Primer ID methodology on 257,846 UDPS reads generated from a HIV-1 SG3Δenv plasmid clone and plasma samples from three HIV-infected patients. The Primer ID consisted of 11 randomized nucleotides, 4,194,304 combinations, in the primer for cDNA synthesis that introduced a unique sequence tag into each cDNA molecule. Consensus template sequences were constructed for reads with Primer IDs that were observed three or more times. Despite high numbers of input template molecules, the number of consensus template sequences was low. With 10,000 input molecules for the clone as few as 97 consensus template sequences were obtained due to highly skewed frequency of resampling. Furthermore, the number of sequenced templates was overestimated due to PCR errors in the Primer IDs. Finally, some consensus template sequences were erroneous due to hotspots for UDPS errors. The Primer ID methodology has the potential to provide highly accurate deep sequencing. However, it is important to be aware that there are remaining challenges with the methodology. In particular it is important to find ways to obtain a more even frequency of resampling of template molecules as well as to identify and remove artefactual consensus template sequences that have been generated by PCR errors in the Primer IDs.

Next-Generation Sequencing of 5' Untranslated Region of Hepatitis C Virus in Search of Minor Viral Variant in a Patient Who Revealed New Genotype While on Antiviral Treatment

The role of mixed infections with different hepatitis C virus (HCV) genotypes in viral persistence, treatment effects, and tissue tropism is unclear. Next-generation sequencing (NGS), which is suitable for analysis of large, genetically diverse populations offers unparalleled advantages for the study of mixed infections. The aim of the study was to determine, using two different deep sequencing strategies (pyrosequencing - 454 Life Sciences/Roche and reversible terminator sequencing-by-synthesis by Illumina), the origin of a novel HCV genotype transiently detectable during antiviral therapy (pre-existing minor population vs. de novo superinfection). Secondly, we compared 5' untranslated region (5'-UTR) variants obtained by the two NGS approaches. 5' UTR amplification products from 9 samples collected from genotype 1b infected patient before, during, and after treatment (4 serum and 5 peripheral blood mononuclear cell - PBMC - samples) were subjected to the next-generation sequencing. The sequencing revealed the presence of two (454/Roche) and one (Illumina) genotype 4 variants in PBMC at Week 16. None of these variants were present either in the preceding or following samples as revealed by both platforms. 454/Roche sequencing detected 24 different 5'-UTR variants: 8 were present in serum and PBMC, 4 only in serum and 12 only in PBMC. Illumina sequencing detected 11 different 5'-UTR variants: 5 in serum and PBMC, 4 only in serum and 2 only in PBMC. Six variants were identical for both sequencing platforms. The difference in variants number was primarily due to variability in two 5'-UTR homopolymeric regions. In conclusion, longitudinal analysis of HCV variants, employing two independent deep sequencing methods, suggests that the transient presence of a different genotype strain in PBMC was a result of superinfection and not a selection of pre-existing minor variant.

Detection of HIV-1 matrix protein p17 quasispecies variants in plasma of chronic HIV-1-infected patients by ultra-deep pyrosequencing

BACKGROUND

The HIV-1 matrix protein p17 (p17MA) is a pleiotropic protein that plays a key role in the HIV-1 life cycle. It has been long believed to have a highly conserved primary amino acid sequence and a well-preserved structural integrity to avoid severe fitness consequences. However, recent data revealed that the carboxy (COOH)-terminus of p17MA possesses high levels of predicted intrinsic disorder, which would subtend to at least partially unfolded status of this region. This finding pointed to the need of investigating p17MA heterogeneity.

METHODS

The degree of intrapatient variations in the p17MA primary sequence was assessed on plasma viral RNA by using ultra-deep pyrosequencing.

RESULTS

Data obtained support a complex nature of p17MA quasispecies, with variants present at variable frequency virtually in all patients. Clusters of mutations were scattered along the entire sequence of the viral protein, but they were more frequently detected within the COOH-terminal region of p17MA. Moreover, deletions and insertions also occurred in a restricted area of the COOH-terminal region.

CONCLUSIONS

On the whole, our data show that the intrapatient level of sequence diversity in the p17MA is much higher than predicted before. Our results pave the way for further studies aimed at unraveling possible correlations between the presence of distinct p17MA variants and peculiar clinical evolutions of HIV-1 disease. The presence of p17MA quasispecies diversity may offer new tools to improve our understanding of the viral adaptation during the natural history of HIV-1 infection.

Minority drug-resistant HIV-1 variants in treatment naïve East-African and Caucasian patients detected by allele-specific real-time PCR

Objective

To assess the presence of two major non-nucleoside reverse transcriptase inhibitors (NNRTI) drug resistance mutations (DRMs), Y181C and K103N, in minor viral quasispecies of treatment naïve HIV-1 infected East-African and Swedish patients by allele-specific polymerase chain reaction (AS-PCR).

Methods

Treatment naïve adults (n = 191) with three epidemiological backgrounds were included: 92 Ethiopians living in Ethiopia; 55 East-Africans who had migrated to Sweden; and 44 Caucasians living in Sweden. The pol gene was analysed by standard population sequencing and by AS-PCR for the detection of Y181C and K103N.

Results

The Y181C was detected in the minority quasispecies of six Ethiopians (6.5%), in two Caucasians (4.5%), and in one East-African (1.8%). The K103N was detected in one East- African (1.8%), by both methods. The proportion of mutants ranged from 0.25% to 17.5%. Additional DRMs were found in all three treatment naïve patient groups by population sequencing.

Conclusions

Major NNRTI mutations can be found by AS-PCR in minor quasispecies of treatment naïve HIV-1 infected Ethiopians living in Ethiopia, in East-African and Caucasian patients living in Sweden in whom population sequencing reveal wild-type virus only. Surveys with standard sequencing are likely to underestimate transmitted drug resistance and the presence of resistant minor quasispecies in treatment naïve patients should be topic for future large scale studies.

Microbial ecology in anaerobic digestion at agitated and non-agitated conditions

To investigate the distribution and dynamics of microbial community in anaerobic digestion at agitated and non-agitated condition, 454 pyrosequencing of 16s rRNA was conducted. It revealed the distinct community compositions between the two digesters and their progressive shifting over time. Methanogens and syntrophic bacteria were found much less abundant in the agitated digester, which was mainly attributed to the presence of bacterial genera Acetanaerobacterium and Ruminococcus with relatively high abundance. The characterization of the microbial community corroborated the digestion performance affected at the agitated condition, where lower methane yield and delayed methane production rate were observed. This was further verified by the accumulation of propionic acid in the agitated digester.

The utility of PacBio circular consensus sequencing for characterizing complex gene families in non-model organisms

BACKGROUND

Molecular characterization of highly diverse gene families can be time consuming, expensive, and difficult, especially when considering the potential for relatively large numbers of paralogs and/or pseudogenes. Here we investigate the utility of Pacific Biosciences single molecule real-time (SMRT) circular consensus sequencing (CCS) as an alternative to traditional cloning and Sanger sequencing PCR amplicons for gene family characterization. We target vomeronasal gene receptors, one of the most diverse gene families in mammals, with the goal of better understanding intra-specific V1R diversity of the gray mouse lemur (Microcebus murinus). Our study compares intragenomic variation for two V1R subfamilies found in the mouse lemur. Specifically, we compare gene copy variation within and between two individuals of M. murinus as characterized by different methods for nucleotide sequencing. By including the same individual animal from which the M. murinus draft genome was derived, we are able to cross-validate gene copy estimates from Sanger sequencing versus CCS methods.

RESULTS

We generated 34,088 high quality circular consensus sequences of two diverse V1R subfamilies (here referred to as V1RI and V1RIX) from two individuals of Microcebus murinus. Using a minimum threshold of 7× coverage, we recovered approximately 90% of V1RI sequences previously identified in the draft M. murinus genome (59% being identical at all nucleotide positions). When low coverage sequences were considered (i.e. < 7× coverage) 100% of V1RI sequences identified in the draft genome were recovered. At least 13 putatively novel V1R loci were also identified using CCS technology.

CONCLUSIONS

Recent upgrades to the Pacific Biosciences RS instrument have improved the CCS technology and offer an alternative to traditional sequencing approaches. Our results suggest that the Microcebus murinus V1R repertoire has been underestimated in the draft genome. In addition to providing an improved understanding of V1R diversity in the mouse lemur, this study demonstrates the utility of CCS technology for characterizing complex regions of the genome. We anticipate that long-read sequencing technologies such as PacBio SMRT will allow for the assembly of multigene family clusters and serve to more accurately characterize patterns of gene copy variation in large gene families, thus revealing novel micro-evolutionary patterns within non-model organisms.

Decreased phenotypic susceptibility to etravirine in patients with predicted genotypic sensitivity

Background

A sensitive, phenotypic reverse transcriptase (RT)-based drug susceptibility assay for the detection of etravirine (ETR) resistance in patient isolates was developed and compared with the results from direct sequencing and ultra-deep pyrosequencing (UDPS).

Methods

Samples were obtained from 15 patients with antiretroviral therapy (ART) failure and from five non-nucleoside reverse transcriptase inhibitor (NNRTI)-naïve patients of whom four were infected by an NNRTI-resistant strain (transmitted drug resistance, TDR). In five patients, two consecutive samples (a and b) were taken for follow up of the virological response. HIV-1 RT was purified and drug susceptibility (IC₅₀) to ETR was estimated. Direct sequencing was performed in all samples and UDPS in samples from nine patients.

Results

Increased IC₅₀ to ETR was found in samples from 13 patients where direct sequencing predicted resistance in only four. UDPS identified additional (N = 11) NNRTI resistance associated mutations (RAMs) in six of nine tested patients. During early failure, IC₅₀ increases were observed in three of six patients without any ETR-RAMs detected by direct sequencing. In further two patients, who stopped NNRTI before sampling, increased IC₅₀ values were found shortly after, despite absence of ETR-RAMs. In two patients who had stopped NNRTI for >1 year, a concordance between phenotype and genotypes was found. Two patients with TDR had increased IC₅₀ despite no ETR-RAMs were detected by direct sequencing. UDPS revealed additional ETR-RAMs in four patients with a discrepancy between phenotype and direct sequencing.

Conclusions

The RT-based phenotypic assay showed decreased ETR susceptibility in patients where direct sequencing predicted ETR-sensitive virus. This increased phenotypic sensitivity was to a large extent supported by UDPS and treatment history. Our method could be valuable for further studies on the phenotypic kinetics of NNRTI resistance. The clinical relevance remains to be studied in larger patient-populations.

Group M consensus Gag and Nef peptides are as efficient at detecting clade A1 and D cross-subtype T-cell functions as subtype-specific consensus peptides

OBJECTIVE

Evaluating HIV-1 specific T-cell response in African populations is sometimes compromised by extensive virus diversity and paucity of non-clade B reagents. We evaluated whether consensus group M (ConM) peptides could serve as comparable substitutes for detecting immune responses in clade A and clade D HIV-1 infection.

METHODS

Frequencies, breadths and polyfunctionality (≥ 3 functions: IFN-γ, IL-2, TNF-α and Perforin) of HIV-specific responses utilizing ConM, ConA and ConD Gag and Nef peptides was compared.

RESULTS

Median genetic distances of infecting gag sequences from consensus group M were (8.9%, IQR 8.2-9.7 and 9%, IQR 3.3-10) for consensus A and D, respectively. Of 24 subjects infected with A and D clade virus, Gag responses were detected in comparable proportions of subjects when using ConM peptides 22/24, ConA peptides 17/24, and ConD peptides 21/24; p=0.12. Nef responses were also detected at similar proportions of subjects when using ConM peptides 15/23, ConA peptides 19/23, and ConD peptides 16/23, p=0.39. Virus-specific CD4+ and CD8+ T-cell polyfunctionality were also detected in similar proportions of infected individuals when using different peptide sets.

CONCLUSIONS

These data support the use of consensus group M overlapping peptide sets as reagents for detecting HIV-specific responses in a clade A and D infected population, but underscore the limitations of utilizing these reagents when evaluating the breadth of virus-specific responses.

Next generation sequencing reveals the hidden diversity of zooplankton assemblages

Background

Zooplankton play an important role in our oceans, in biogeochemical cycling and providing a food source for commercially important fish larvae. However, difficulties in correctly identifying zooplankton hinder our understanding of their roles in marine ecosystem functioning, and can prevent detection of long term changes in their community structure. The advent of massively parallel next generation sequencing technology allows DNA sequence data to be recovered directly from whole community samples. Here we assess the ability of such sequencing to quantify richness and diversity of a mixed zooplankton assemblage from a productive time series site in the Western English Channel.

Methodology/Principle Findings

Plankton net hauls (200 µm) were taken at the Western Channel Observatory station L4 in September 2010 and January 2011. These samples were analysed by microscopy and metagenetic analysis of the 18S nuclear small subunit ribosomal RNA gene using the 454 pyrosequencing platform. Following quality control a total of 419,041 sequences were obtained for all samples. The sequences clustered into 205 operational taxonomic units using a 97% similarity cut-off. Allocation of taxonomy by comparison with the National Centre for Biotechnology Information database identified 135 OTUs to species level, 11 to genus level and 1 to order, <2.5% of sequences were classified as unknowns. By comparison a skilled microscopic analyst was able to routinely enumerate only 58 taxonomic groups.

Conclusions

Metagenetics reveals a previously hidden taxonomic richness, especially for Copepoda and hard-to-identify meroplankton such as Bivalvia, Gastropoda and Polychaeta. It also reveals rare species and parasites. We conclude that Next Generation Sequencing of 18S amplicons is a powerful tool for elucidating the true diversity and species richness of zooplankton communities. While this approach allows for broad diversity assessments of plankton it may become increasingly attractive in future if sequence reference libraries of accurately identified individuals are better populated.

Genetic analysis of West Nile virus isolates from an outbreak in Idaho, United States, 2006-2007

West Nile virus (WNV) appeared in the U.S. in 1999 and has since become endemic, with yearly summer epidemics causing tens of thousands of cases of serious disease over the past 14 years. Analysis of WNV strains isolated during the 2006–2007 epidemic seasons demonstrates that a new genetic variant had emerged coincidentally with an intense outbreak in Idaho during 2006. The isolates belonging to the new variant carry a 13 nt deletion, termed ID-Δ13, located at the variable region of the 3′UTR, and are genetically related. The analysis of deletions and insertions in the 3′UTR of two major lineages of WNV revealed the presence of conserved repeats and two indel motifs in the variable region of the 3′UTR. One human and two bird isolates from the Idaho 2006–2007 outbreaks were sequenced using Illumina technology and within-host variability was analyzed. Continued monitoring of new genetic variants is important for public health as WNV continues to evolve.

PCR-induced transitions are the major source of error in cleaned ultra-deep pyrosequencing data

Background

Ultra-deep pyrosequencing (UDPS) is used to identify rare sequence variants. The sequence depth is influenced by several factors including the error frequency of PCR and UDPS. This study investigated the characteristics and source of errors in raw and cleaned UDPS data.

Results

UDPS of a 167-nucleotide fragment of the HIV-1 SG3Δenv plasmid was performed on the Roche/454 platform. The plasmid was diluted to one copy, PCR amplified and subjected to bidirectional UDPS on three occasions. The dataset consisted of 47,693 UDPS reads. Raw UDPS data had an average error frequency of 0.30% per nucleotide site. Most errors were insertions and deletions in homopolymeric regions. We used a cleaning strategy that removed almost all indel errors, but had little effect on substitution errors, which reduced the error frequency to 0.056% per nucleotide. In cleaned data the error frequency was similar in homopolymeric and non-homopolymeric regions, but varied considerably across sites. These site-specific error frequencies were moderately, but still significantly, correlated between runs (r = 0.15–0.65) and between forward and reverse sequencing directions within runs (r = 0.33–0.65). Furthermore, transition errors were 48-times more common than transversion errors (0.052% vs. 0.001%; p<0.0001). Collectively the results indicate that a considerable proportion of the sequencing errors that remained after data cleaning were generated during the PCR that preceded UDPS.

Conclusions

A majority of the sequencing errors that remained after data cleaning were introduced by PCR prior to sequencing, which means that they will be independent of platform used for next-generation sequencing. The transition vs. transversion error bias in cleaned UDPS data will influence the detection limits of rare mutations and sequence variants.

Comparative genome analysis of ciprofloxacin-resistant Pseudomonas aeruginosa reveals genes within newly identified high variability regions associated with drug resistance development

The alarming rise of ciprofloxacin-resistant Pseudomonas aeruginosa has been reported in several clinical studies. Though the mutation of resistance genes and their role in drug resistance has been researched, the process by which the bacterium acquires high-level resistance is still not well understood. How does the genomic evolution of P. aeruginosa affect resistance development? Could the exposure of antibiotics to the bacteria enrich genomic variants that lead to the development of resistance, and if so, how are these variants distributed through the genome? To answer these questions, we performed 454 pyrosequencing and a whole genome analysis both before and after exposure to ciprofloxacin. The comparative sequence data revealed 93 unique resistance strain variation sites, which included a mutation in the DNA gyrase subunit A gene. We generated variation-distribution maps comparing the wild and resistant types, and isolated 19 candidates from three discrete resistance-associated high variability regions that had available transposon mutants, to perform a ciprofloxacin exposure assay. Of these region candidates with transposon disruptions, 79% (15/19) showed a reduction in the ability to gain high-level resistance, suggesting that genes within these high variability regions might enrich for certain functions associated with resistance development.

SIV genome-wide pyrosequencing provides a comprehensive and unbiased view of variation within and outside CD8 T lymphocyte epitopes

Deep sequencing technology is revolutionizing our understanding of HIV/SIV evolution. It is known that acute SIV sequence variation within CD8 T lymphocyte (CD8-TL) epitopes is similar among MHC-identical animals, but we do not know whether this persists into the chronic phase. We now determine whether chronic viral variation in MHC-identical animals infected with clonal SIV is similar throughout the entire coding sequence when using a sensitive deep sequencing approach. We pyrosequenced the entire coding sequence of the SIV genome isolated from a unique cohort of four SIVmac239-infected, MHC-identical Mauritian cynomolgus macaques (MCM) 48 weeks after infection; one MCM in the cohort became an elite controller. Among the three non-controllers, we found that genome-wide sequences were similar between animals and we detected increased sequence complexity within 64% of CD8-TL epitopes when compared to Sanger sequencing methods. When we compared sequences between the MHC-matched controller and the three non-controllers, we found the viral population in the controller was less diverse and accumulated different variants than the viral populations in the non-controllers. Importantly, we found that initial PCR amplification of viral cDNA did not significantly affect the sequences detected, suggesting that data obtained by pyrosequencing PCR-amplified viral cDNA accurately represents the diversity of sequences replicating within an animal. This demonstrates that chronic sequence diversity across the entire SIV coding sequence is similar among MHC-identical animals with comparable viral loads when infected with the same clonal virus stock. Additionally, our approach to genome-wide SIV sequencing accurately reflects the diversity of sequences present in the replicating viral population. In sum, our study suggests that genome-wide pyrosequencing of immunodeficiency viruses captures a thorough and unbiased picture of sequence diversity, and may be a useful approach to employ when evaluating which sequences to include as part of a vaccine immunogen.

Comparison of a high-resolution melting assay to next-generation sequencing for analysis of HIV diversity

Next-generation sequencing (NGS) has recently been used for analysis of HIV diversity, but this method is labor-intensive, costly, and requires complex protocols for data analysis. We compared diversity measures obtained using NGS data to those obtained using a diversity assay based on high-resolution melting (HRM) of DNA duplexes. The HRM diversity assay provides a single numeric score that reflects the level of diversity in the region analyzed. HIV gag and env from individuals in Rakai, Uganda, were analyzed in a previous study using NGS (n = 220 samples from 110 individuals). Three sequence-based diversity measures were calculated from the NGS sequence data (percent diversity, percent complexity, and Shannon entropy). The amplicon pools used for NGS were analyzed with the HRM diversity assay. HRM scores were significantly associated with sequence-based measures of HIV diversity for both gag and env (P < 0.001 for all measures). The level of diversity measured by the HRM diversity assay and NGS increased over time in both regions analyzed (P < 0.001 for all measures except for percent complexity in gag), and similar amounts of diversification were observed with both methods (P < 0.001 for all measures except for percent complexity in gag). Diversity measures obtained using the HRM diversity assay were significantly associated with those from NGS, and similar increases in diversity over time were detected by both methods. The HRM diversity assay is faster and less expensive than NGS, facilitating rapid analysis of large studies of HIV diversity and evolution.

Interplay between HIV-1 and Host Genetic Variation: A Snapshot into Its Impact on AIDS and Therapy Response

As of February 2012, 50 circulating recombinant forms (CRFs) have been reported for HIV-1 while one CRF for HIV-2. Also according to HIV sequence compendium 2011, the HIV sequence database is replete with 414,398 sequences. The fact that there are CRFs, which are an amalgamation of sequences derived from six or more subtypes (CRF27_cpx (cpx refers to complex) is a mosaic with sequences from 6 different subtypes besides an unclassified fragment), serves as a testimony to the continual divergent evolution of the virus with its approximate 1% per year rate of evolution, and this phenomena per se poses tremendous challenge for vaccine development against HIV/AIDS, a devastating disease that has killed 1.8 million patients in 2010. Here, we explore the interaction between HIV-1 and host genetic variation in the context of HIV/AIDS and antiretroviral therapy response.

Ras oncogenes in oral cancer: the past 20 years

Oral squamous cell carcinoma (OSCC) of head and neck is associated with high morbidity and mortality in both Western and Asian countries. Several risk factors for the development of oral cancer are very well established, including tobacco chewing, betel quid, smoking, alcohol drinking and human papilloma virus (HPV) infection. Apart from these risk factors, many genetic factors such as oncogenes, tumor suppressor genes and regulatory genes are identified to involve in oral carcinogenesis with these risk factors dependent and independent manner. Ras is one of the most frequently genetically deregulated oncogene in oral cancer. In this review, we analyze the past 22years of literature on genetic alterations such as mutations and amplifications of the isoforms of the ras oncogene in oral cancer. Further, we addressed the isoform-specific role of the ras in oral carcinogenesis. We also discussed how targeting the Akt and MEK, downstream effectors of the PI3K/Akt and MAPK pathways, respectively, would probably pave the possible molecular therapeutic target for the ras driven tumorigenesis in oral cancer. Analysis of these ras isoforms may critically enlighten specific role of a particular ras isoform in oral carcinogenesis, enhance prognosis and pave the way for isoform-specific molecular targeted therapy in OSCC.