Sunando Roy

Comparison of Next-Generation Sequencing Technologies for Comprehensive Assessment of Full-Length Hepatitis C Viral Genomes.

ABSTRACT Affordable next-generation sequencing (NGS) technologies for hepatitis C virus (HCV) may potentially identify both viral genotype and resistance genetic motifs in the era of directly acting antiviral (DAA) therapies. This study compared the ability of high-throughput NGS methods to generate full-length, deep, HCV sequence data sets and evaluated their utility for diagnostics and clinical assessment. NGS methods using (i) unselected HCV RNA (metagenomics), (ii) preenrichment of HCV RNA by probe capture, and (iii) HCV preamplification by PCR implemented in four United Kingdom centers were compared. Metrics of sequence coverage and depth, quasispecies diversity, and detection of DAA resistance-associated variants (RAVs), mixed HCV genotypes, and other coinfections were compared using a panel of samples with different viral loads, genotypes, and mixed HCV genotypes/subtypes [geno(sub)types]. Each NGS method generated near-complete genome sequences from more than 90% of samples. Enrichment methods and PCR preamplification generated greater sequence depth and were more effective for samples with low viral loads. All NGS methodologies accurately identified mixed HCV genotype infections. Consensus sequences generated by different NGS methods were generally concordant, and majority RAVs were consistently detected. However, methods differed in their ability to detect minor populations of RAVs. Metagenomic methods identified human pegivirus coinfections. NGS provided a rapid, inexpensive method for generating whole HCV genomes to define infecting genotypes, RAVs, comprehensive viral strain analysis, and quasispecies diversity. Enrichment methods are particularly suited for high-throughput analysis while providing the genotype and information on potential DAA resistance.

Norovirus Whole-Genome Sequencing by SureSelect Target Enrichment: a Robust and Sensitive Method

ABSTRACT Norovirus full-genome sequencing is challenging due to sequence heterogeneity among genomes. Previous methods have relied on PCR amplification, which is problematic due to primer design, and transcriptome sequencing (RNA-Seq), which nonspecifically sequences all RNA, including host and bacterial RNA, in stool specimens. Target enrichment uses a panel of custom-designed 120-mer RNA baits that are complementary to all publicly available norovirus sequences, with multiple baits targeting each position of the genome, which overcomes the challenge of primer design. Norovirus genomes are enriched from stool RNA extracts to minimize the sequencing of nontarget RNA. SureSelect target enrichment and Illumina sequencing were used to sequence full genomes from 507 norovirus-positive stool samples with reverse transcription–real-time PCR cycle threshold (CT ) values of 10 to 43. Sequencing on an Illumina MiSeq system in batches of 48 generated, on average, 81% on-target reads per sample and 100% genome coverage with >12,000-fold read depth. Samples included genotypes GI.1, GI.2, GI.3, GI.6, GI.7, GII.1, GII.2, GII.3, GII.4, GII.5, GII.6, GII.7, GII.13, GII.14, and GII.17. When outliers were accounted for, we generated >80% genome coverage for all positive samples, regardless of CT values. A total of 164 samples were tested in parallel with conventional PCR genotyping of the capsid shell domain; 164/164 samples were successfully sequenced, compared to 158/164 samples that were amplified by PCR. Four of the samples that failed capsid PCR analysis had low titers, which suggests that target enrichment is more sensitive than gel-based PCR. Two samples failed PCR due to primer mismatches; target enrichment uses multiple baits targeting each position, thus accommodating sequence heterogeneity among norovirus genomes.

The emerging GII.P16-GII.4 Sydney 2012 norovirus lineage is circulating worldwide, arose by late-2014 and contains polymerase changes that may increase virus transmission

Noroviruses are a leading cause of human gastroenteritis worldwide. The norovirus genotype GII.4 is the most prevalent genotype in the human population and has caused six pandemics since 1995. A novel norovirus lineage containing the GII.P16 polymerase and pandemic GII.4 Sydney 2012 capsid was recently detected in Asia and Germany. We demonstrate that this lineage is also circulating within the UK and USA and has been circulating since October 2014 or earlier. While the lineage does not contain unique substitutions in the capsid, it does contain polymerase substitutions close to positions known to influence polymerase function and virus transmission. These polymerase substitutions are shared with a GII.P16-GII.2 virus that dominated outbreaks in Germany in Winter 2016. We suggest that the substitutions in the polymerase may have resulted in a more transmissible virus and the combination of this polymerase and the pandemic GII.4 capsid may result in a highly transmissible virus. Further surveillance efforts will be required to determine whether the GII.P16-GII.4 Sydney 2012 lineage increases in frequency over the coming months.

Multivariate Statistical Analyses Demonstrate Unique Host Immune Responses to Single and Dual Lentiviral Infection

Background Feline immunodeficiency virus (FIV) and human immunodeficiency virus (HIV) are recently identified lentiviruses that cause progressive immune decline and ultimately death in infected cats and humans. It is of great interest to understand how to prevent immune system collapse caused by these lentiviruses. We recently described that disease caused by a virulent FIV strain in cats can be attenuated if animals are first infected with a feline immunodeficiency virus derived from a wild cougar. The detailed temporal tracking of cat immunological parameters in response to two viral infections resulted in high-dimensional datasets containing variables that exhibit strong co-variation. Initial analyses of these complex data using univariate statistical techniques did not account for interactions among immunological response variables and therefore potentially obscured significant effects between infection state and immunological parameters. Methodology and Principal Findings Here, we apply a suite of multivariate statistical tools, including Principal Component Analysis, MANOVA and Linear Discriminant Analysis, to temporal immunological data resulting from FIV superinfection in domestic cats. We investigated the co-variation among immunological responses, the differences in immune parameters among four groups of five cats each (uninfected, single and dual infected animals), and the “immune profiles” that discriminate among them over the first four weeks following superinfection. Dual infected cats mount an immune response by 24 days post superinfection that is characterized by elevated levels of CD8 and CD25 cells and increased expression of IL4 and IFNγ, and FAS. This profile discriminates dual infected cats from cats infected with FIV alone, which show high IL-10 and lower numbers of CD8 and CD25 cells. Conclusions Multivariate statistical analyses demonstrate both the dynamic nature of the immune response to FIV single and dual infection and the development of a unique immunological profile in dual infected cats, which are protected from immune decline.

On the effective depth of viral sequence data.

Abstract Genome sequence data are of great value in describing evolutionary processes in viral populations. However, in such studies, the extent to which data accurately describes the viral population is a matter of importance. Multiple factors may influence the accuracy of a dataset, including the quantity and nature of the sample collected, and the subsequent steps in viral processing. To investigate this phenomenon, we sequenced replica datasets spanning a range of viruses, and in which the point at which samples were split was different in each case, from a dataset in which independent samples were collected from a single patient to another in which all processing steps up to sequencing were applied to a single sample before splitting the sample and sequencing each replicate. We conclude that neither a high read depth nor a high template number in a sample guarantee the precision of a dataset. Measures of consistency calculated from within a single biological sample may also be insufficient; distortion of the composition of a population by the experimental procedure or genuine within-host diversity between samples may each affect the results. Where it is possible, data from replicate samples should be collected to validate the consistency of short-read sequence data.

Use of Whole-Genome Sequencing of Adenovirus in Immunocompromised Pediatric Patients to Identify Nosocomial Transmission and Mixed-Genotype Infection

Background Adenoviruses are significant pathogens for the immunocompromised, arising from primary infection or reinfection. Serotyping is insufficient to support nosocomial transmission investigations. We investigate whether whole-genome sequencing (WGS) provides clinically relevant information on transmission among patients in a pediatric tertiary hospital. Methods We developed a target-enriched adenovirus WGS technique for clinical samples and retrospectively sequenced 107 adenovirus-positive residual diagnostic samples, including viremias (>5 × 104 copies/mL), from 37 patients collected January 2011-March 2016. Whole-genome sequencing was used to determine genotype and for phylogenetic analysis. Results Adenovirus sequences were recovered from 105 of 107 samples. Full genome sequences were recovered from all 20 nonspecies C samples and from 36 of 85 species C viruses, with partial genome sequences recovered from the rest. Whole-genome phylogenetic analysis suggested linkage of 3 genotype A31 cases and uncovered an unsuspected epidemiological link to an A31 infection first detected on the same ward 4 years earlier. In 9 samples from 1 patient who died, we identified a mixed genotype adenovirus infection. Conclusions Adenovirus WGS from clinical samples is possible and useful for genotyping and molecular epidemiology. Whole-genome sequencing identified likely nosocomial transmission with greater resolution than conventional genotyping and distinguished between adenovirus disease due to single or multiple genotypes.

Molecular characterization of a human G20P[28] rotavirus a strain with multiple genes related to bat rotaviruses

Group A rotaviruses are the major cause of severe gastroenteritis in the young of mammals and birds. This report describes characterization of an unusual G20P[28] rotavirus strain detected in a 24month old child from Suriname. Genomic sequence analyses revealed that the genotype constellation of the Suriname strain RVA/Human-wt/SUR/2014735512/2013/G20P[28] was G20-P[28]-I13-R13-C13-M12-A23-N13-T15-E20-H15. Genes VP1, VP2, VP3, NSP1, NSP2, NSP3, NSP4 and NSP5 were recently assigned novel genotypes by the Rotavirus Classification Working Group (RCWG). Three of the 11 gene segments (VP7, VP4, VP6) were similar to cognate gene sequences of bat-like human rotavirus strain Ecu534 from Ecuador and the VP7, NSP3 and NSP5 gene segments of strain RVA/Human-wt/SUR/2014735512/2013/G20P[28] were found to be closely related to gene sequences of bat rotavirus strain 3081/BRA detected in Brazil. Although distantly related, the VP1 gene of the study strain and bat strain BatLi09 detected in Cameroon in 2014 are monophyletic. The NSP1 gene was found to be most closely related to human strain QUI-35-F5 from Brazil. These findings suggest that strain RVA/Human-wt/SUR/2014735512/2013/G20P[28] represents a zoonotic infection from a bat host.

Super-infections and relapses occur in chronic norovirus infections

BACKGROUND Norovirus causes chronic infections in immunocompromised patients with considerable associated morbidity. It is not known whether chronic infections involve super- or re-infections or relapses. OBJECTIVES To retrospectively investigate whether longitudinal sampling in chronically infected patients demonstrates persistent infection with the same virus, or super- or re-infection. STUDY DESIGN Norovirus full genomes were generated from 86 longitudinal samples from 25 paediatric patients. Consensus sequences were used for phylogenetic analysis and genotyping. RESULTS Super-infections occurred in 17% of chronically infected patients who were continuously PCR positive; including two with mixed norovirus infections. The median duration of infection was 107days longer in those with super-infections; however this was not statistically significant. A third of patients with interrupted norovirus shedding continued to be infected with the same virus despite up to 2 months of PCR negative stools, classified as a relapse. The majority (67%) of patients with interrupted shedding were re-infected with a different genotype. CONCLUSIONS Chronically infected patients who are continuously PCR positive are most likely to remain infected with the same virus; however super-infections do occur leading to mixed infection. Patients with interrupted shedding are likely to represent re-infection with a different genotype, however relapsing infections also occur. Our findings have implications for infection control as immunosuppressed patients remain susceptible to new norovirus infections despite current or recent infection and may continue to be infectious after norovirus is undetectable in stool. The relevance to children without co-morbidities remains to be determined.

Haplotype assignment of virus NGS data using co-variation of variant frequencies

Longitudinal deep sequencing of viruses can provide detailed information about intra-host evolutionary dynamics including how viruses interact with and transmit between hosts. Many analyses require haplotype reconstruction, identifying which variants are co-located on the same genomic element. Most current methods to perform this reconstruction are based on a high density of variants and cannot perform this reconstruction for slowly evolving viruses. We present a new approach, HaROLD (HAplotype Reconstruction Of Longitudinal Deep sequencing data), which performs this reconstruction based on identifying co-varying variant frequencies using a probabilistic framework. We test this method with synthetic data sets of mixed cytomegalovirus and norovirus genomes, demonstrating high accuracy when longitudinal samples are available.

Norovirus Transmission Dynamics in a Pediatric Hospital Using Full Genome Sequences

Norovirus genome sequencing identified 33% of patients whose sequences were linked phylogenetically to those of another patient in the study. An additional 24% of nosocomially infected patients had unlinked sequences, suggesting infection from unsampled sources. Genome sequencing identifies unsuspected nosocomial norovirus infection.