Single-cell landscape of functionally cured chronic hepatitis B patients reveals activation of innate and altered CD4-CTL-driven adaptive immunity.

BACKGROUND & AIMS: Hepatitis B surface antigen (HBsAg) loss or functional cure (FC) is considered the optimal therapeutic outcome for patients with chronic hepatitis B (CHB). However, the immune-pathological biomarkers and underlying mechanisms of FC remain unclear. In this study we comprehensively interrogate disease-associated cell states identified within intrahepatic tissue and matched PBMCs (peripheral blood mononuclear cells) from patients with CHB or after FC, at the resolution of single cells, to provide novel insights into putative mechanisms underlying FC. METHODS: We combined single-cell transcriptomics (single-cell RNA sequencing) with multiparametric flow cytometry-based immune phenotyping, and multiplexed immunofluorescence to elucidate the immunopathological cell states associated with CHB vs. FC. RESULTS: We found that the intrahepatic environment in CHB and FC displays specific cell identities and molecular signatures that are distinct from those found in matched PBMCs. FC is associated with the emergence of an altered adaptive immune response marked by CD4 cytotoxic T lymphocytes, and an activated innate response represented by liver-resident natural killer cells, specific Kupffer cell subtypes and marginated neutrophils. Surprisingly, we found MHC class II-expressing hepatocytes in patients achieving FC, as well as low but persistent levels of covalently closed circular DNA and pregenomic RNA, which may play an important role in FC. CONCLUSIONS: Our study provides conceptually novel insights into the immuno-pathological control of HBV cure, and opens exciting new avenues for clinical management, biomarker discovery and therapeutic development. We believe that the discoveries from this study, as it relates to the activation of an innate and altered immune response that may facilitate sustained, low-grade inflammation, may have broader implications in the resolution of chronic viral hepatitis. IMPACT AND IMPLICATIONS: This study dissects the immuno-pathological cell states associated with functionally cured chronic hepatitis B (defined by the loss of HBV surface antigen or HBsAg). We identified the sustained presence of very low viral load, accessory antigen-presenting hepatocytes, adaptive-memory-like natural killer cells, and the emergence of helper CD4 T cells with cytotoxic or effector-like signatures associated with functional cure, suggesting previously unsuspected alterations in the adaptive immune response, as well as a key role for the innate immune response in achieving or maintaining functional cure. Overall, the insights generated from this study may provide new avenues for the development of alternative therapies as well as patient surveillance for better clinical management of chronic hepatitis B.

Patterns of mutation within an emerging endemic lineage of HEV-3a.

The hepatitis E virus can cause chronic infections in immuno-suppressed patients, and cases have been on the rise globally. Viral mutations during such infections are difficult to characterize. We deep-sequenced viral populations from 15 immunocompromised patients with chronic HEV to identify the viral lineage and describe viral mutational hotspots within and across patients. A total of 21 viral RNA samples were collected between 2012 and 2017 from a single hospital in Singapore. Sequences covering a total of 3894 bp of the HEV genome were obtained. Phylogenetic analyses identified all sequences as belonging to the HEV-3a sub-clade and clearly indicate a unique local lineage. Deep sequencing reveals variable viral population complexity during infections. Comparisons of viral samples from the same patients spaced 2-19 months apart revealed rapid nucleotide replacements in the dominant viral sequence in both ribavirin treated and treatment-naive patients. Mutational hotspots were identified within ORF3 and the PCP/HVR domain of ORF1.

Single-virion sequencing of lamivudine-treated HBV populations reveal population evolution dynamics and demographic history.

BACKGROUND: Viral populations are complex, dynamic, and fast evolving. The evolution of groups of closely related viruses in a competitive environment is termed quasispecies. To fully understand the role that quasispecies play in viral evolution, characterizing the trajectories of viral genotypes in an evolving population is the key. In particular, long-range haplotype information for thousands of individual viruses is critical; yet generating this information is non-trivial. Popular deep sequencing methods generate relatively short reads that do not preserve linkage information, while third generation sequencing methods have higher error rates that make detection of low frequency mutations a bioinformatics challenge. Here we applied BAsE-Seq, an Illumina-based single-virion sequencing technology, to eight samples from four chronic hepatitis B (CHB) patients - once before antiviral treatment and once after viral rebound due to resistance. RESULTS: With single-virion sequencing, we obtained 248-8796 single-virion sequences per sample, which allowed us to find evidence for both hard and soft selective sweeps. We were able to reconstruct population demographic history that was independently verified by clinically collected data. We further verified four of the samples independently through PacBio SMRT and Illumina Pooled deep sequencing. CONCLUSIONS: Overall, we showed that single-virion sequencing yields insight into viral evolution and population dynamics in an efficient and high throughput manner. We believe that single-virion sequencing is widely applicable to the study of viral evolution in the context of drug resistance and host adaptation, allows differentiation between soft or hard selective sweeps, and may be useful in the reconstruction of intra-host viral population demographic history.

Chronic Infection With Camelid Hepatitis E Virus in a Liver Transplant Recipient Who Regularly Consumes Camel Meat and Milk.

There have been increasing reports of food-borne zoonotic transmission of hepatitis E virus (HEV) genotype 3, which causes chronic infections in immunosuppressed patients. We performed phylogenetic analyses of the HEV sequence (partial and full-length) from 1 patient from the Middle East who underwent liver transplantation, and compared it with other orthohepevirus A sequences. We found the patient to be infected by camelid HEV. This patient regularly consumed camel meat and milk, therefore camelid HEV, which is genotype 7, might infect human beings. Our finding links consumption of camel-derived food products to post-transplantation hepatitis E, which, if detected at early stages, can be cured with antiviral therapy and reduced administration of immunosuppressive agents.

LoFreq: a sequence-quality aware, ultra-sensitive variant caller for uncovering cell-population heterogeneity from high-throughput sequencing datasets.

The study of cell-population heterogeneity in a range of biological systems, from viruses to bacterial isolates to tumor samples, has been transformed by recent advances in sequencing throughput. While the high-coverage afforded can be used, in principle, to identify very rare variants in a population, existing ad hoc approaches frequently fail to distinguish true variants from sequencing errors. We report a method (LoFreq) that models sequencing run-specific error rates to accurately call variants occurring in <0.05% of a population. Using simulated and real datasets (viral, bacterial and human), we show that LoFreq has near-perfect specificity, with significantly improved sensitivity compared with existing methods and can efficiently analyze deep Illumina sequencing datasets without resorting to approximations or heuristics. We also present experimental validation for LoFreq on two different platforms (Fluidigm and Sequenom) and its application to call rare somatic variants from exome sequencing datasets for gastric cancer. Source code and executables for LoFreq are freely available at http://sourceforge.net/projects/lofreq/.

A randomized, double-blind placebo controlled trial of balapiravir, a polymerase inhibitor, in adult dengue patients.

BACKGROUND: Dengue is the most common arboviral infection of humans. There are currently no specific treatments for dengue. Balapiravir is a prodrug of a nucleoside analogue (called R1479) and an inhibitor of hepatitis C virus replication in vivo. METHODS: We conducted in vitro experiments to determine the potency of balapiravir against dengue viruses and then an exploratory, dose-escalating, randomized placebo-controlled trial in adult male patients with dengue with <48 hours of fever. RESULTS: The clinical and laboratory adverse event profile in patients receiving balapiravir at doses of 1500 mg (n = 10) or 3000 mg (n = 22) orally for 5 days was similar to that of patients receiving placebo (n = 32), indicating balapiravir was well tolerated. However, twice daily assessment of viremia and daily assessment of NS1 antigenemia indicated balapiravir did not measurably alter the kinetics of these virological markers, nor did it reduce the fever clearance time. The kinetics of plasma cytokine concentrations and the whole blood transcriptional profile were also not attenuated by balapiravir treatment. CONCLUSIONS: Although this trial, the first of its kind in dengue, does not support balapiravir as a candidate drug, it does establish a framework for antiviral treatment trials in dengue and provides the field with a clinically evaluated benchmark molecule. CLINICAL TRIALS REGISTRATION: NCT01096576.

Large-scale evolutionary surveillance of the 2009 H1N1 influenza A virus using resequencing arrays.

In April 2009, a new influenza A (H1N1 2009) virus emerged that rapidly spread around the world. While current variants of this virus have caused widespread disease, particularly in vulnerable groups, there remains the possibility that future variants may cause increased virulence, drug resistance or vaccine escape. Early detection of these virus variants may offer the chance for increased containment and potentially prevention of the virus spread. We have developed and field-tested a resequencing kit that is capable of interrogating all eight segments of the 2009 influenza A(H1N1) virus genome and its variants, with added focus on critical regions such as drug-binding sites, structural components and mutation hotspots. The accompanying base-calling software (EvolSTAR) introduces novel methods that utilize neighbourhood hybridization intensity profiles and substitution bias of probes on the microarray for mutation confirmation and recovery of ambiguous base queries. Our results demonstrate that EvolSTAR is highly accurate and has a much improved call rate. The high throughput and short turn-around time from sample to sequence and analysis results (30 h for 24 samples) makes this kit an efficient large-scale evolutionary biosurveillance tool.

Dengue virus surveillance for early warning, Singapore.

In Singapore, after a major outbreak of dengue in 2005, another outbreak occurred in 2007. Laboratory-based surveillance detected a switch from dengue virus serotype 1 (DENV-1) to DENV-2. Phylogenetic analysis showed a clade replacement within DENV-2 cosmopolitan genotype, which accompanied the predominant serotype switch, and cocirculation of multiple genotypes of DENV-3.

Genomic epidemiology of a dengue virus epidemic in urban Singapore.

Dengue is one of the most important emerging diseases of humans, with no preventative vaccines or antiviral cures available at present. Although one-third of the world's population live at risk of infection, little is known about the pattern and dynamics of dengue virus (DENV) within outbreak situations. By exploiting genomic data from an intensively studied major outbreak, we are able to describe the molecular epidemiology of DENV at a uniquely fine-scaled temporal and spatial resolution. Two DENV serotypes (DENV-1 and DENV-3), and multiple component genotypes, spread concurrently and with similar epidemiological and evolutionary profiles during the initial outbreak phase of a major dengue epidemic that took place in Singapore during 2005. Although DENV-1 and DENV-3 differed in viremia and clinical outcome, there was no evidence for adaptive evolution before, during, or after the outbreak, indicating that ecological or immunological rather than virological factors were the key determinants of epidemic dynamics.

Multifactorial heterogeneity of virus-specific T cells and association with the progression of human chronic hepatitis B infection.

Associations between chronic antigen stimulation, T cell dysfunction, and the expression of various inhibitory receptors are well characterized in several mouse and human systems. During chronic hepatitis B virus (HBV) infection (CHB), T cell responses are blunted with low frequencies of virus-specific T cells observed, making these parameters difficult to study. Here, using mass cytometry and a highly multiplexed combinatorial peptide-major histocompatibility complex (pMHC) tetramer strategy that allows for the detection of rare antigen-specific T cells, we simultaneously probed 484 unique HLA-A*1101-restricted epitopes spanning the entire HBV genome on T cells from patients at various stages of CHB. Numerous HBV-specific T cell populations were detected, validated, and profiled. T cells specific for two epitopes (HBVpol387 and HBVcore169) displayed differing and complex heterogeneities that were associated with the disease progression, and the expression of inhibitory receptors on these cells was not linearly related with their extent of T cell dysfunction. For HBVcore169-specific CD8+ T cells, we found cellular markers associated with long-term memory, polyfunctionality, and the presence of several previously unidentified public TCR clones that correlated with viral control. Using high-dimensional trajectory analysis of these cellular phenotypes, a pseudo-time metric was constructed that fit with the status of viral infection in corresponding patients. This was validated in a longitudinal cohort of patients undergoing antiviral therapy. Our study uncovers complex relationships of inhibitory receptors between the profiles of antigen-specific T cells and the status of CHB with implications for new strategies of therapeutic intervention.

Dengue virus not detected in human semen.

Zika and Ebola viruses can persist in semen and pose a risk for sexual transmission. To determine if dengue virus, another flavivirus, is similarly detectable in semen, we performed dengue PCR on semen in five patients with acute dengue virus infection. All five tested negative, suggesting that dengue does not pose a risk for sexual transmission.

Extended Evaluation of Virological, Immunological and Pharmacokinetic Endpoints of CELADEN: A Randomized, Placebo-Controlled Trial of Celgosivir in Dengue Fever Patients.

UNLABELLED: CELADEN was a randomized placebo-controlled trial of 50 patients with confirmed dengue fever to evaluate the efficacy and safety of celgosivir (A study registered at ClinicalTrials.gov, number NCT01619969). Celgosivir was given as a 400 mg loading dose and 200 mg bid (twice a day) over 5 days. Replication competent virus was measured by plaque assay and compared to reverse transcription quantitative PCR (qPCR) of viral RNA. Pharmacokinetics (PK) correlations with viremia, immunological profiling, next generation sequence (NGS) analysis and hematological data were evaluated as exploratory endpoints here to identify possible signals of pharmacological activity. Viremia by plaque assay strongly correlated with qPCR during the first four days. Immunological profiling demonstrated a qualitative shift in T helper cell profile during the course of infection. NGS analysis did not reveal any prominent signature that could be associated with drug treatment; however the phylogenetic spread of patients' isolates underlines the importance of strain variability that may potentially confound interpretation of dengue drug trials conducted during different outbreaks and in different countries. Celgosivir rapidly converted to castanospermine (Cast) with mean peak and trough concentrations of 5727 ng/mL (30.2 µM) and 430 ng/mL (2.3 µM), respectively and cleared with a half-life of 2.5 (± 0.6) hr. Mean viral log reduction between day 2 and 4 (VLR2-4) was significantly greater in secondary dengue than primary dengue (p = 0.002). VLR2-4 did not correlate with drug AUC but showed a trend of greater response with increasing Cmin. PK modeling identified dosing regimens predicted to achieve 2.4 to 4.5 times higher Cmin. than in the CELADEN trial for only 13% to 33% increase in overall dose. A small, non-statistical trend towards better outcome on platelet nadir and difference between maximum and minimum hematocrit was observed in celgosivir-treated patients with secondary dengue infection. Optimization of the dosing regimen and patient stratification may enhance the ability of a clinical trial to demonstrate celgosivir activity in treating dengue fever based on hematological endpoints. A new clinical trial with a revised dosing regimen is slated to start in 2016 (NCT02569827). Furthermore celgosivir's potential value for treatment of other flaviruses such as Zika virus should be investigated urgently. TRIAL REGISTRATION: ClinicalTrials.gov NCT01619969.

Tracking Dengue Virus Intra-host Genetic Diversity during Human-to-Mosquito Transmission.

Dengue virus (DENV) infection of an individual human or mosquito host produces a dynamic population of closely-related sequences. This intra-host genetic diversity is thought to offer an advantage for arboviruses to adapt as they cycle between two very different host species, but it remains poorly characterized. To track changes in viral intra-host genetic diversity during horizontal transmission, we infected Aedes aegypti mosquitoes by allowing them to feed on DENV2-infected patients. We then performed whole-genome deep-sequencing of human- and matched mosquito-derived DENV samples on the Illumina platform and used a sensitive variant-caller to detect single nucleotide variants (SNVs) within each sample. >90% of SNVs were lost upon transition from human to mosquito, as well as from mosquito abdomen to salivary glands. Levels of viral diversity were maintained, however, by the regeneration of new SNVs at each stage of transmission. We further show that SNVs maintained across transmission stages were transmitted as a unit of two at maximum, suggesting the presence of numerous variant genomes carrying only one or two SNVs each. We also present evidence for differences in selection pressures between human and mosquito hosts, particularly on the structural and NS1 genes. This analysis provides insights into how population drops during transmission shape RNA virus genetic diversity, has direct implications for virus evolution, and illustrates the value of high-coverage, whole-genome next-generation sequencing for understanding viral intra-host genetic diversity.

Next-generation whole genome sequencing of dengue virus.

RNA viruses are notorious for their ability to quickly adapt to selective pressure from the host immune system and/or antivirals. This adaptability is likely due to the error-prone characteristics of their RNA-dependent, RNA polymerase [1, 2]. Dengue virus, a member of the Flaviviridae family of positive-strand RNA viruses, is also known to share these error-prone characteristics [3]. Utilizing high-throughput, massively parallel sequencing methodologies, or next-generation sequencing (NGS), we can now accurately quantify these populations of viruses and track the changes to these populations over the course of a single infection. The aim of this chapter is twofold: to describe the methodologies required for sample preparation prior to sequencing and to describe the bioinformatics analyses required for the resulting data.

BAsE-Seq: a method for obtaining long viral haplotypes from short sequence reads.

We present a method for obtaining long haplotypes, of over 3 kb in length, using a short-read sequencer, Barcode-directed Assembly for Extra-long Sequences (BAsE-Seq). BAsE-Seq relies on transposing a template-specific barcode onto random segments of the template molecule and assembling the barcoded short reads into complete haplotypes. We applied BAsE-Seq on mixed clones of hepatitis B virus and accurately identified haplotypes occurring at frequencies greater than or equal to 0.4%, with >99.9% specificity. Applying BAsE-Seq to a clinical sample, we obtained over 9,000 viral haplotypes, which provided an unprecedented view of hepatitis B virus population structure during chronic infection. BAsE-Seq is readily applicable for monitoring quasispecies evolution in viral diseases.

A method for full genome sequencing of all four serotypes of the dengue virus.

The availability of whole genome sequencing has contributed to many aspects of dengue research, and its use in dengue virus (DENV) surveillance for early epidemic warning has been proposed. Methods to sequence the genomes of individual dengue serotypes have been described previously, but no single method is known to be applicable for all four serotypes. This report describes a method for sequencing the entire genome of all four DENV serotypes. Using tagged oligonucleotide primers designed for the 3' end, viral RNA was reverse transcribed into a cDNA spanning the entire genome of each of the four serotypes (DENV-1 to -4). This was followed by amplification of the entire cDNA in five overlapping amplicons. A sequence tag was added to the sense primer annealing to the 5' UTR sequence and the antisense primer annealing to the 3' UTR sequence to ensure no terminal nucleotides were omitted during PCR. Sixty-one virus isolates were sequenced: 58 DENV-2, one DENV-1, one DENV-4 and one DENV-3 published previously. The method described could be applied readily for viral biology studies and incorporated into proactive dengue virologic surveillance.

Host gene expression profiling of dengue virus infection in cell lines and patients.

BACKGROUND: Despite the seriousness of dengue-related disease, with an estimated 50-100 million cases of dengue fever and 250,000-500,000 cases of dengue hemorrhagic fever/dengue shock syndrome each year, a clear understanding of dengue pathogenesis remains elusive. Because of the lack of a disease model in animals and the complex immune interaction in dengue infection, the study of host response and immunopathogenesis is difficult. The development of genomics technology, microarray and high throughput quantitative PCR have allowed researchers to study gene expression changes on a much broader scale. We therefore used this approach to investigate the host response in dengue virus-infected cell lines and in patients developing dengue fever. METHODOLOGY/PRINCIPAL FINDINGS: Using microarray and high throughput quantitative PCR method to monitor the host response to dengue viral replication in cell line infection models and in dengue patient blood samples, we identified differentially expressed genes along three major pathways; NF-kappaB initiated immune responses, type I interferon (IFN) and the ubiquitin proteasome pathway. Among the most highly upregulated genes were the chemokines IP-10 and I-TAC, both ligands of the CXCR3 receptor. Increased expression of IP-10 and I-TAC in the peripheral blood of ten patients at the early onset of fever was confirmed by ELISA. A highly upregulated gene in the IFN pathway, viperin, was overexpressed in A549 cells resulting in a significant reduction in viral replication. The upregulation of genes in the ubiquitin-proteasome pathway prompted the testing of proteasome inhibitors MG-132 and ALLN, both of which reduced viral replication. CONCLUSION/SIGNIFICANCE: Unbiased gene expression analysis has identified new host genes associated with dengue infection, which we have validated in functional studies. We showed that some parts of the host response can be used as potential biomarkers for the disease while others can be used to control dengue viral replication, thus representing viable targets for drug therapy.