Whole genome deep sequencing analysis of viral quasispecies diversity and evolution in HBeAg seroconverters.

BACKGROUND & AIMS: We aimed to investigate how viral quasispecies of the HBV whole genome evolves and diversifies in response to HBeAg seroconversion and viral control utilising next-generation sequencing (NGS). METHODS: Fifty HBeAg-positive chronic hepatitis B patients, including 18 treatment-naïve and 32 interferon (IFN)-treated individuals, were recruited. Serial HBV whole genomes in serum were analysed by NGS to determine sequence characteristics and viral quasispecies. RESULTS: HBV quasispecies diversity, measured by nucleotide diversity, was negatively correlated with viral load and hepatitis activity. Spontaneous HBeAg seroconverters exhibited significantly greater viral quasispecies diversity than treatment-naïve non-seroconverters from >1 year before seroconversion (0.0112 vs. 0.0060, p <0.01) to >1 year after seroconversion (0.0103 vs. 0.0068, p <0.01). IFN-induced HBeAg seroconverters tended to have higher viral genetic diversity than non-seroconverters along with treatment. Particularly, the IFN responders, defined as IFN-induced HBeAg seroconversion with low viraemia, exhibited significantly greater genetic diversity of whole HBV genome at 6 months post-IFN treatment than IFN non-responders (0.0148 vs. 0.0106, p = 0.048). Moreover, spontaneous HBeAg seroconverters and IFN responders exhibited significantly higher evolutionary rates and more intra-host single-nucleotide variants. Interestingly, in spontaneous HBeAg seroconverters and IFN responders, there were distinct evolutionary patterns in the HBV genome. CONCLUSIONS: Higher HBV quasispecies diversity is associated with spontaneous HBeAg seroconversion and IFN-induced HBeAg seroconversion with low viraemia, conferring a favourable clinical outcome. LAY SUMMARY: HBeAg seroconversion is a landmark in the natural history of chronic HBV infection. Using next-generation sequencing, we found that the nucleotide diversity of HBV was negatively correlated with viral load and hepatitis activity. Patients undergoing HBeAg seroconversion had more diverse HBV genomes and a faster viral evolution rate. Our findings suggest HBeAg seroconversion is driven by host selection pressure, likely immune selection pressure.

Combined viral quasispecies diversity and hepatitis B core-related antigen predict off-nucleos(t)ide analog durability in HBeAg-negative patients.

BACKGROUND: Viral quasispecies dynamics between pre- and post-nucleos(t)ide analog (NA) therapy remains unclear. AIM: This study aimed to investigate the HBV quasispecies evolution and its relationship with durability of off-therapy responses in HBeAg-negative chronic hepatitis B (CHB) patients who stopped NA therapy. METHODS: Fifty-four HBeAg-negative CHB patients who stopped NAs, including 19 virological controllers (VC) who maintained serum HBV DNA < 2000 IU/mL beyond 1-year off-therapy, and 35 virological relapsers (VR) experiencing virological relapse within 1-year off-therapy were recruited. Viral quasispecies was analyzed by deep sequencing. Hepatitis B core-related antigen (HBcrAg) and HBsAg were also measured. RESULTS: VC had significantly higher baseline viral quasispecies diversity of the precore/core gene, measured by nucleotide diversity, than VR. Low baseline viral nucleotide diversity (< 0.01) and high HBcrAg (≧ 2.0 KU/mL), but not HBsAg, at end of treatment (EOT) were significantly associated with higher risk of 1-year virological relapse (hazard ratio [HR] 6.09 and 3.31, respectively). Combination of low baseline viral nucleotide diversity and high HBcrAg at EOT could identify patients at high risk (HR 15.82). Further analysis of the evolution of HBV whole genome showed that HBV nucleotide diversity negatively correlated with serum HBV DNA levels. Notably, the viral quasispecies diversity between pre- and post-NA treatment remained relatively unchanged. CONCLUSION: Higher baseline HBV quasispecies diversity associates with more durable off-therapy viral suppression in HBeAg-negative CHB patients. Combination of baseline viral nucleotide diversity and HBcrAg at EOT can identify patients at high risk for virological relapse after stopping NAs.

Functional analysis of mutations in a severe congenital neutropenia syndrome caused by glucose-6-phosphatase-ß deficiency.

Glucose-6-phosphatase-ß (G6Pase-ß or G6PC3) deficiency is characterized by neutropenia and dysfunction in both neutrophils and macrophages. G6Pase-ß is an enzyme embedded in the endoplasmic reticulum membrane that catalyzes the hydrolysis of glucose-6-phosphate (G6P) to glucose and phosphate. To date, 33 separate G6PC3 mutations have been identified in G6Pase-ß-deficient patients but only the p.R253H and p.G260R missense mutations have been characterized functionally for pathogenicity. Here we functionally characterize 16 of the 19 known missense mutations using a sensitive assay, based on a recombinant adenoviral vector-mediated expression system, to demonstrate pathogenicity. Fourteen missense mutations completely abolish G6Pase-ß enzymatic activity while the p.S139I and p.R189Q mutations retain 49% and 45%, respectively of wild type G6Pase-ß activity. A database of residual enzymatic activity retained by the G6Pase-ß mutations will serve as a reference for evaluating genotype-phenotype relationships.

The CRISPR/Cas9 System Facilitates Clearance of the Intrahepatic HBV Templates In Vivo.

Persistence of hepatitis B virus (HBV) covalently closed circular DNA (cccDNA) under current antiviral therapy is a major barrier to eradication of chronic hepatitis B (CHB). Curing CHB will require novel strategies for specific disruption of cccDNA. The clustered regularly interspaced short palindromic repeats (CRISPR)/Cas9 system is a newly developed tool for site-specific cleavage of DNA targets directed by a synthetic guide RNA (gRNA) base-paired to the target DNA sequence. To examine whether this system can cleave HBV genomes, we designed eight gRNAs against HBV of genotype A. With the HBV-specific gRNAs, the CRISPR/Cas9 system significantly reduced the production of HBV core and surface proteins in Huh-7 cells transfected with an HBV-expression vector. Among eight screened gRNAs, two effective ones were identified. Interestingly, one gRNA targeting the conserved HBV sequence acted against different genotypes. Using a hydrodynamics-HBV persistence mouse model, we further demonstrated that this system could cleave the intrahepatic HBV genome-containing plasmid and facilitate its clearance in vivo, resulting in reduction of serum surface antigen levels. These data suggest that the CRISPR/Cas9 system could disrupt the HBV-expressing templates both in vitro and in vivo, indicating its potential in eradicating persistent HBV infection.

Prevention of hepatocellular adenoma and correction of metabolic abnormalities in murine glycogen storage disease type Ia by gene therapy.

UNLABELLED: Glycogen storage disease type Ia (GSD-Ia), which is characterized by impaired glucose homeostasis and chronic risk of hepatocellular adenoma (HCA), is caused by deficiencies in the endoplasmic reticulum (ER)-associated glucose-6-phosphatase-α (G6Pase-α or G6PC) that hydrolyzes glucose-6-phosphate (G6P) to glucose. G6Pase-α activity depends on the G6P transporter (G6PT) that translocates G6P from the cytoplasm into the ER lumen. The functional coupling of G6Pase-α and G6PT maintains interprandial glucose homeostasis. We have shown previously that gene therapy mediated by AAV-GPE, an adeno-associated virus (AAV) vector expressing G6Pase-α directed by the human G6PC promoter/enhancer (GPE), completely normalizes hepatic G6Pase-α deficiency in GSD-Ia (G6pc(-/-) ) mice for at least 24 weeks. However, a recent study showed that within 78 weeks of gene deletion, all mice lacking G6Pase-α in the liver develop HCA. We now show that gene therapy mediated by AAV-GPE maintains efficacy for at least 70-90 weeks for mice expressing more than 3% of wild-type hepatic G6Pase-α activity. The treated mice displayed normal hepatic fat storage, had normal blood metabolite and glucose tolerance profiles, had reduced fasting blood insulin levels, maintained normoglycemia over a 24-hour fast, and had no evidence of hepatic abnormalities. After a 24-hour fast, hepatic G6PT messenger RNA levels in G6pc(-/-) mice receiving gene therapy were markedly increased. Because G6PT transport is the rate-limiting step in microsomal G6P metabolism, this may explain why the treated G6pc(-/-) mice could sustain prolonged fasts. The low fasting blood insulin levels and lack of hepatic steatosis may explain the absence of HCA. CONCLUSION: These results confirm that AAV-GPE-mediated gene transfer corrects hepatic G6Pase-α deficiency in murine GSD-Ia and prevents chronic HCA formation.

SLC37A1 and SLC37A2 are phosphate-linked, glucose-6-phosphate antiporters.

Blood glucose homeostasis between meals depends upon production of glucose within the endoplasmic reticulum (ER) of the liver and kidney by hydrolysis of glucose-6-phosphate (G6P) into glucose and phosphate (P(i)). This reaction depends on coupling the G6P transporter (G6PT) with glucose-6-phosphatase-α (G6Pase-α). Only one G6PT, also known as SLC37A4, has been characterized, and it acts as a P(i)-linked G6P antiporter. The other three SLC37 family members, predicted to be sugar-phosphate:P(i) exchangers, have not been characterized functionally. Using reconstituted proteoliposomes, we examine the antiporter activity of the other SLC37 members along with their ability to couple with G6Pase-α. G6PT- and mock-proteoliposomes are used as positive and negative controls, respectively. We show that SLC37A1 and SLC37A2 are ER-associated, P(i)-linked antiporters, that can transport G6P. Unlike G6PT, neither is sensitive to chlorogenic acid, a competitive inhibitor of physiological ER G6P transport, and neither couples to G6Pase-α. We conclude that three of the four SLC37 family members are functional sugar-phosphate antiporters. However, only G6PT/SLC37A4 matches the characteristics of the physiological ER G6P transporter, suggesting the other SLC37 proteins have roles independent of blood glucose homeostasis.

The helical domains of the stem region of dengue virus envelope protein are involved in both virus assembly and entry.

The envelope (E) of dengue virus (DENV) is a determinant of tropism and virulence. At the C terminus of E protein, there is a stem region containing two amphipathic α-helical domains (EH1 and EH2) and a stretch of conserved sequences in between. The crystal structure of E protein at the postfusion state suggested the involvement of the stem during the fusion; however, the critical domains or residues involved remain unknown. Site-directed mutagenesis was carried out to replace each of the stem residues at the hydrophobic face with an alanine or proline in a DENV serotype 4 (DENV4) precursor membrane (prM)/E expression construct. Most of the 15 proline mutations at either EH1 or EH2 severely affected the assembly of virus-like particles (VLPs). Radioimmunoprecipitation and membrane flotation assays revealed that EH1 mutations primarily affect prM-E heterodimerization and EH2 mutations affect the membrane binding of the stem. Introducing four proline mutations at either EH1 or EH2 into a DENV2 replicon packaging system greatly affects assembly and entry. Moreover, introducing these mutations into a DENV2 infectious clone confirmed the impairment in assembly and infectivity. Sequencing analysis of adaptive mutations in passage 5 viruses revealed a change to a leucine or wild-type residue at the original site, suggesting the importance of maintaining the helical structure. Collectively, these findings suggest that the EH1 and EH2 domains are involved in both assembly and entry steps of the DENV replication cycle; this feature, together with the high degree of sequence conservation, suggests that the stem region is a potential target of antiviral strategies.

High levels of mecA DNA detected by a quantitative real-time PCR assay are associated with mortality in patients with methicillin-resistant Staphylococcus aureus bacteremia.

Persistent methicillin-resistant Staphylococcus aureus (MRSA) bacteremia is known to be a poor prognostic factor. While several PCR assays for the detection of MRSA in various clinical samples were recently reported, the possibility that a quantitative PCR assay could be used to quantify and monitor MRSA bacteremia has not been explored. In this study, we established a quantitative real-time PCR assay for the mecA gene using known copy numbers of a plasmid containing mecA DNA as a standard and the previously described mecA-specific primers and probe (P. Francois et al., J. Clin. Microbiol. 41:254-260, 2003). We employed this assay to examine 250 sequential whole-blood samples from 20 adult patients, including 13 survivors and 7 nonsurvivors, with culture-proven MRSA bacteremia at the intensive care units of National Taiwan University Hospital between 1 July 2006 and 31 January 2007. The levels of mecA DNA in the nonsurvivors were significantly higher than those in the survivors during the three periods of bacteremia examined (days 0 to 2, 3 to 5, and 6 to 8) (P = 0.003 by two-tailed Mann-Whitney U test). Moreover, the nonsurvivors had higher mecA DNA levels than the survivors after 3 days and 7 days of anti-MRSA therapy (medians for nonsurvivors and survivors at 3 days, 5.86 and 4.30 log copies/ml, respectively; medians for nonsurvivors and survivors at 7 days, 5.21 and 4.36 log copies/ml, respectively; P = 0.02 and P = 0.04, respectively, by two-tailed Mann-Whitney U test). Together, these findings suggest that the level of mecA DNA in blood could potentially be used to monitor MRSA bacteremia and evaluate responses to therapy.

Evolution of dengue virus type 2 during two consecutive outbreaks with an increase in severity in southern Taiwan in 2001-2002.

To investigate viral determinants and evolution linked to outbreak with increased severity, we examined dengue virus type 2 (DENV-2) sequences from plasma of 31 patients (14 dengue fever, 17 dengue hemorrhagic fever, DHF) continuously during the 2001 and 2002 outbreaks in southern Taiwan, in which both the total cases and proportion of DHF cases increased. Analysis of envelope (E) and full-genome sequences between viruses of the two outbreaks revealed 5 nucleotide changes in E, NS1, NS4A, and NS5 genes. None was identical to those reported in the DENV-2 outbreak in Cuba in 1997, suggesting viral determinants linked to severe outbreak are genotype dependent. Compared with previous reports of lineage turnover years apart, our findings that the 2002 viruses descended from a minor variant of the 2001 viruses in less than 6 months was novel, and may represent a mechanism of evolution of DENV from one outbreak to another.

Antibodies to envelope glycoprotein of dengue virus during the natural course of infection are predominantly cross-reactive and recognize epitopes containing highly conserved residues at the fusion loop of domain II.

The antibody response to the envelope (E) glycoprotein of dengue virus (DENV) is known to play a critical role in both protection from and enhancement of disease, especially after primary infection. However, the relative amounts of homologous and heterologous anti-E antibodies and their epitopes remain unclear. In this study, we examined the antibody responses to E protein as well as to precursor membrane (PrM), capsid, and nonstructural protein 1 (NS1) of four serotypes of DENV by Western blot analysis of DENV serotype 2-infected patients with different disease severity and immune status during an outbreak in southern Taiwan in 2002. Based on the early-convalescent-phase sera tested, the rates of antibody responses to PrM and NS1 proteins were significantly higher in patients with secondary infection than in those with primary infection. A blocking experiment and neutralization assay showed that more than 90% of anti-E antibodies after primary infection were cross-reactive and nonneutralizing against heterologous serotypes and that only a minor proportion were type specific, which may account for the type-specific neutralization activity. Moreover, the E-binding activity in sera of 10 patients with primary infection was greatly reduced by amino acid replacements of three fusion loop residues, tryptophan at position 101, leucine at position 107, and phenylalanine at position 108, but not by replacements of those outside the fusion loop of domain II, suggesting that the predominantly cross-reactive anti-E antibodies recognized epitopes involving the highly conserved residues at the fusion loop of domain II. These findings have implications for our understanding of the pathogenesis of dengue and for the future design of subunit vaccine against DENV as well.

Study of sequence variation of dengue type 3 virus in naturally infected mosquitoes and human hosts: implications for transmission and evolution.

Dengue virus is an arbovirus that replicates alternately in the mosquito vector and human host. We investigated sequences of dengue type 3 virus in naturally infected Aedes aegypti mosquitoes and in eight patients from the same outbreak and reported that the extent of sequence variation seen with the mosquitoes was generally lower than that seen with the patients (mean diversity, 0.21 versus 0.38% and 0.09 versus 0.23% for the envelope [E] and capsid [C] genes, respectively). This was further verified with five experimentally infected mosquitoes (mean diversity, 0.09 and 0.10% for the E and C genes, respectively). Examination of the quasispecies structures of the E sequences of the mosquitoes and of the patients revealed that the sequences of the major variants were the same, suggesting that the major variant was transmitted. These findings support our hypothesis that mosquitoes contribute to the evolutionary conservation of dengue virus by maintaining a more homogenous viral population and a dominant variant during transmission.

Concurrent infections by two dengue virus serotypes among dengue patients in Taiwan.

The co-circulation of multiple dengue virus serotypes in the same region has been reported in several countries in Southeast Asia as well as in Central and South America for decades. Although outbreaks involving more than one serotype of dengue virus have been reported in Taiwan since 1987, concurrent infection in the same individual by multiple serotypes of dengue virus have never been identified. Using a modified multiplex reverse transcription-polymerase chain reaction assay, we detected and determined the serotypes of 21 dengue patients during an outbreak in southern Taiwan in 2000. While either dengue type 2 or type 3 virus was found in most cases, 2 of the 21 cases were concurrent infections by both dengue type 2 and dengue type 3 viruses. This was further confirmed by sequence analysis of the amplified products. This study reports for the first time that concurrent infections occur in Taiwan, and suggests that future virologic surveillance of dengue outbreak in Taiwan should take this into consideration.

Dengue type 3 virus in plasma is a population of closely related genomes: quasispecies.

Using reverse transcription-PCR and clonal sequencing of the dengue virus envelope gene derived from the plasma samples of six patients, we reported for the first time that dengue virus circulates as a population of closely related genomes. The extent of sequence diversity varied among patients, with the mean pairwise proportions of difference ranging from 0.21 to 1.67%. Genome-defective viruses were found in 5.8% of the total number of clones analyzed. Our findings on the quasispecies nature of dengue virus and the defective virus in vivo have implications with regard to the pathogenesis of dengue virus.