Deep Sequencing of the HIV-1 Polymerase Gene for Characterisation of Cytotoxic T-Lymphocyte Epitopes During Early and Chronic Disease Stages

Background Despite multiple attempts, there is still no effective HIV-1 vaccine available. The HIV-1 polymerase (pol) gene is highly conserved and encodes cytotoxic T-lymphocyte (CTL) epitopes. The aim of the study was to characterise HIV-1 Pol CTL epitopes in mostly sample pairs obtained during early and chronic stages of infection. Methods Illumina deep sequencing was conducted for all samples while Sanger sequencing was only performed on baseline samples. Codons under immune selection pressure were assessed by computing nonsynonymous to synonymous mutation ratios using MEGA. Minority CTL epitope variants occurring at \(\ge\)5% were detected using low-frequency variant tool in CLC Genomics. Los Alamos HIV database was used for mapping mutations to known HIV-1 CTL epitopes. Results Fifty-two participants were enrolled in the study. Their median age was 28 years (interquartile range: 24–32 years) and majority of participants (92.3%) were female. Illumina minority variant analysis identified a significantly higher number of CTL epitopes (n = 65) compared to epitopes (n = 8) identified through Sanger sequencing. Most of the identified epitopes mapped to reverse transcriptase (RT) and integrase (IN) regardless of sequencing method. There was a significantly higher proportion of minority variant epitopes in RT (n = 39, 60.0%) compared to IN (n = 17, 26.2%) and PR (n = 9, 13.8%), p = 0.002 and < 0.0001, respectively. However, no significant difference was observed between the proportion of minority variant epitopes in IN versus PR, p = 0.06. Some epitopes were detected in either early or chronic HIV-1 infection whereas others were detected in both stages. Different distribution patterns of minority variant epitopes were observed in sample pairs; with some increasing or decreasing over time, while others remained constant. Some of the identified epitopes have not been previously reported for HIV-1 subtype C. There were also variants that could not be mapped to reported CTL epitopes in the Los Alamos HIV database. Conclusion Deep sequencing revealed many Pol CTL epitopes, including some not previously reported for HIV-1 subtype C. The findings of this study support the inclusion of RT and IN epitopes in HIV-1 vaccine candidates as these proteins harbour many CTL epitopes.

Deep sequencing of the HIV-1 polymerase gene for characterisation of cytotoxic T-lymphocyte epitopes during early and chronic disease stages

Background Despite multiple attempts, there is still no effective HIV-1 vaccine available. The HIV-1 polymerase (pol) gene is highly conserved and encodes cytotoxic T-lymphocyte (CTL) epitopes. In this study, deep sequencing was employed for characterisation of HIV-1 Pol CTL epitopes in mostly paired samples obtained during early and chronic stages of infection. Deep sequencing data was then compared to Sanger sequencing data only in samples obtained at baseline. Results Fifty-two participants were enrolled in the study. Their median age was 28 years (interquartile range: 24–32 years) and the majority of participants (92.3%) were female. Illumina minority variant analysis identified a significantly higher number of CTL epitopes (n = 65) compared to epitopes (n = 8) identified through Sanger sequencing. Most of the identified epitopes mapped to reverse transcriptase (RT) and integrase (IN) regardless of the method of sequencing. There was a significantly higher proportion of minority variant epitopes in RT (n = 39, 60.0%) compared to IN (n = 17, 26.2%) and PR (n = 9, 13.8%), p = 0.002 and < 0.0001, respectively. However, no significant difference was observed between the proportion of minority variant epitopes in IN versus PR, p = 0.06. Some epitopes were detected in either early or chronic HIV-1 infection whereas others were detected in both stages. Different distribution patterns of minority variant epitopes were observed in sample pairs; with some increasing or decreasing over time, while others remained constant. Some of the identified epitopes have not been previously reported for HIV-1 subtype C. There were also variants that could not be mapped to reported CTL epitopes in the Los Alamos HIV database. Conclusion Deep sequencing revealed many Pol CTL epitopes, including some not previously reported for HIV-1 subtype C. The findings of this study support the inclusion of RT and IN epitopes in HIV-1 vaccine candidates as these proteins harbour many CTL epitopes. Variants that were not mapped within CTL epitopes could represent new epitopes.

A Low Fidelity Virus Shows Increased Recombination during the Removal of an Alphavirus Reporter Gene

Reporter genes for RNA viruses are well-known to be unstable due to putative RNA recombination events that excise inserted nucleic acids. RNA recombination has been demonstrated to be co-regulated with replication fidelity in alphaviruses, but it is unknown how recombination events at the minority variant level act, which is important for vaccine and trans-gene delivery design. Therefore, we sought to characterize the removal of a reporter gene by a low-fidelity alphavirus mutant over multiple replication cycles. To examine this, GFP was inserted into TC-83, a live-attenuated vaccine for the alphavirus Venezuelan equine encephalitis virus, as well as a low-fidelity variant of TC-83, and passaged until fluorescence was no longer observed. Short-read RNA sequencing using ClickSeq was performed to determine which regions of the viral genome underwent recombination and how this changed over multiple replication cycles. A rapid removal of the GFP gene was observed, where minority variants in the virus population accumulated small deletions that increased in size over the course of passaging. Eventually, these small deletions merged to fully remove the GFP gene. The removal was significantly enhanced during the passaging of low-fidelity TC-83, suggesting that increased levels of recombination are a defining characteristic of this mutant.

Evaluating the performance of tools used to call minority variants from whole genome short-read data

Background: High-throughput whole genome sequencing facilitates investigation of minority virus sub-populations from virus positive samples. Minority variants are useful in understanding within and between host diversity, population dynamics and can potentially assist in elucidating person-person transmission pathways. Several minority variant callers have been developed to describe low frequency sub-populations from whole genome sequence data. These callers differ based on bioinformatics and statistical methods used to discriminate sequencing errors from low-frequency variants. Methods: We evaluated the diagnostic performance and concordance between published minority variant callers used in identifying minority variants from whole-genome sequence data from virus samples. We used the ART-Illumina read simulation tool to generate three artificial short-read datasets of varying coverage and error profiles from an RSV reference genome. The datasets were spiked with nucleotide variants at predetermined positions and frequencies. Variants were called using FreeBayes, LoFreq, Vardict, and VarScan2. The variant callers’ agreement in identifying known variants was quantified using two measures; concordance accuracy and the inter-caller concordance. Results: The variant callers reported differences in identifying minority variants from the datasets. Concordance accuracy and inter-caller concordance were positively correlated with sample coverage. FreeBayes identified the majority of variants although it was characterised by variable sensitivity and precision in addition to a high false positive rate relative to the other minority variant callers and which varied with sample coverage. LoFreq was the most conservative caller. Conclusions: We conducted a performance and concordance evaluation of four minority variant calling tools used to identify and quantify low frequency variants. Inconsistency in the quality of sequenced samples impacts on sensitivity and accuracy of minority variant callers. Our study suggests that combining at least three tools when identifying minority variants is useful in filtering errors when calling low frequency variants.

Evaluating the performance of tools used to call minority variants from whole genome short-read data

Background: High-throughput whole genome sequencing facilitates investigation of minority sub-populations from virus positive samples. Minority variants are useful in understanding within and between host diversity, population dynamics and can potentially help to elucidate person-person transmission chains. Several minority variant callers have been developed to describe the minority variants sub-populations from whole genome sequence data. However, they differ on bioinformatics and statistical approaches used to discriminate sequencing errors from low-frequency variants. Methods: We evaluated the diagnostic performance and concordance between published minority variant callers used in identifying minority variants from whole-genome sequence data. The ART-Illumina read simulation tool was used to generate three artificial short-read datasets of varying coverage and error profiles from an RSV reference genome. The datasets were spiked with nucleotide variants at predetermined positions and frequencies. Variants were called using FreeBayes, LoFreq, Vardict, and VarScan2. The variant callers’ agreement in identifying known variants was quantified using two measures; concordance accuracy and the inter-caller concordance. Results: The variant callers reported differences in identifying minority variants from the datasets. Concordance accuracy and inter-caller concordance were positively correlated with sample coverage. FreeBayes identified majority of the variants although it was characterised by variable sensitivity and precision in addition to a high false positive rate relative to the other minority variant callers and which varied with sample coverage. LoFreq was the most conservative caller. Conclusions: We conducted a performance and concordance evaluation of four minority variant calling tools used to identify and quantify low frequency variants. Inconsistency in the quality of sequenced samples impact on sensitivity and accuracy of minority variant callers. Our study suggests that combining at least three tools when identifying minority variants is useful in filtering errors when calling low frequency variants.

The unstoppable glottal: tracking rapid change in an iconic British variable

This article presents a sociolinguistic investigation of a rapidly expanding innovation in the UK, glottal replacement, in a variety spoken in northeast Scotland. Quantitative analysis of the form shows a dramatic change in apparent time: from a minority variant in the older generation to a full 90 per cent use in the younger generation. Further analysis of the constraints on use provide a detailed snapshot of how this variant moves through social and linguistic space. Males use higher rates of the non-standard form in the older generations but this constraint is neutralised in the younger generations as the form increases. Styleshifting according to interlocutor also neutralises through time. While these results across the social constraints are in line with previous analyses, the linguistic constraints differ in this variety. In contrast to most other varieties, intervocalic contexts such asbottleshow high rates of glottal replacement. Moreover, word-internal foot-initial contexts (e.g.sometimes) also frequently allow the non-standard variant, despite this being rare in other dialects. Although glottal replacement is largely considered to be a ‘torchbearer’ of geographical diffusion, this in-depth analysis suggests that different varieties may have different pathways of change in the rapid transition from [t] to [ʔ] throughout the UK.

Viral population analysis and minority-variant detection using short read next-generation sequencing

RNA viruses within infected individuals exist as a population of evolutionary-related variants. Owing to evolutionary change affecting the constitution of this population, the frequency and/or occurrence of individual viral variants can show marked or subtle fluctuations. Since the development of massively parallel sequencing platforms, such viral populations can now be investigated to unprecedented resolution. A critical problem with such analyses is the presence of sequencing-related errors that obscure the identification of true biological variants present at low frequency. Here, we report the development and assessment of the Quality Assessment of Short Read (QUASR) Pipeline ( http://sourceforge.net/projects/quasr ) specific for virus genome short read analysis that minimizes sequencing errors from multiple deep-sequencing platforms, and enables post-mapping analysis of the minority variants within the viral population. QUASR significantly reduces the error-related noise in deep-sequencing datasets, resulting in increased mapping accuracy and reduction of erroneous mutations. Using QUASR, we have determined influenza virus genome dynamics in sequential samples from an in vitro evolution of 2009 pandemic H1N1 (A/H1N1/09) influenza from samples sequenced on both the Roche 454 GSFLX and Illumina GAIIx platforms. Importantly, concordance between the 454 and Illumina sequencing allowed unambiguous minority-variant detection and accurate determination of virus population turnover in vitro .