Detection of hyper-conserved regions in hepatitis B virus X gene potentially useful for gene therapy

AIM To detect hyper-conserved regions in the hepatitis B virus(HBV) X gene(HBX) 5' region that could be candidates for gene therapy.METHODS The study included 27 chronic hepatitis B treatmentnaive patients in various clinical stages(from chronic infection to cirrhosis and hepatocellular carcinoma, both HBeA g-negative and HBeA g-positive), and infected with HBV genotypes A-F and H. In a serum sample from each patient with viremia > 3.5 log IU/m L, the HBX 5' end region [nucleotide(nt) 1255-1611] was PCRamplified and submitted to next-generation sequencing(NGS). We assessed genotype variants by phylogenetic analysis, and evaluated conservation of this region by calculating the information content of each nucleotide position in a multiple alignment of all unique sequences(haplotypes) obtained by NGS. Conservation at the HBx protein amino acid(aa) level was also analyzed.RESULTS NGS yielded 1333069 sequences from the 27 samples, with a median of 4578 sequences/sample(2487-9279, IQR 2817). In 14/27 patients(51.8%), phylogenetic analysis of viral nucleotide haplotypes showed a complex mixture of genotypic variants. Analysis of the information content in the haplotype multiple alignments detected 2 hyper-conserved nucleotide regions, one in the HBX upstream non-coding region(nt 1255-1286) and the other in the 5' end coding region(nt 1519-1603). This last region coded for a conserved amino acid region(aa 63-76) that partially overlaps a Kunitz-like domain.CONCLUSION Two hyper-conserved regions detected in the HBX 5' end may be of value for targeted gene therapy, regardless of the patients' clinical stage or HBV genotype.

Characterization of hepatitis B virus X gene quasispecies complexity in mono-infection and hepatitis delta virus superinfection

Hepatitis delta virus(HDV) seems to strongly suppress hepatitis B virus(HBV)replication, although little is known about the mechanism of this interaction. Both these viruses show a dynamic distribution of mutants, resulting in viral quasispecies. Next-generation sequencing is a viable approach for analyzing the composition of these mutant spectra. As the regulatory hepatitis B X protein(HBx) is essential for HBV replication, determination of HBV X gene(HBX)quasispecies complexity in HBV/HDV infection compared to HBV monoinfection may provide information on the interactions between these two viruses.AIM To compare HBV quasispecies complexity in the HBX 5' region between chronic hepatitis delta(CHD) and chronic HBV mono-infected patients.METHODS Twenty-four untreated patients were included: 7/24(29.2%) with HBeAgnegative chronic HBV infection(CI, previously termed inactive carriers), 8/24(33.3%) with HBeAg-negative chronic hepatitis B(CHB) and 9/24(37.5%) with CHD. A serum sample from each patient was first tested for HBV DNA levels.The HBX 5' region [nucleotides(nt) 1255-1611] was then PCR-amplified for subsequent next-generation sequencing(MiSeq, Illumina, United States). HBV quasispecies complexity in the region analyzed was evaluated using incidencebased indices(number of haplotypes and number of mutations), abundancebased indices(Hill numbers of order 1 and 2), and functional indices(mutation frequency and nucleotide diversity). We also evaluated the pattern of nucleotide changes to investigate which of them could be the cause of the quasispecies complexity.RESULTS CHB patients showed higher median HBV-DNA levels [5.4 logIU/mL,interquartile range(IQR) 3.5-7.9] than CHD(3.4 logIU/mL, IQR 3-7.6)(P = n.s.)or CI(3.2 logIU/mL, IQR 2.3-3.5)(P < 0.01) patients. The incidence and abundance indices indicated that HBV quasispecies complexity was significantly greater in CI than CHB. A similar trend was observed in CHD patients, although only Hill numbers of order 2 showed statistically significant differences(CHB2.81, IQR 1.11-4.57 vs CHD 8.87, 6.56-11.18, P = 0.038). There were no significant differences in the functional indices, but CI and CHD patients also showed a trend towards greater complexity than CHB. No differences were found for any HBV quasispecies complexity indices between CHD and CI patients. G-to-A and C-to-T nucleotide changes, characteristic of APOBEC3 G, were higher in CHD and CI than in CHB in genotype A haplotypes, but not in genotype D. The proportion of nt G-to-A vs A-to-G changes and C-to-T vs T-to-C changes in genotype A and D haplotypes in CHD patients showed no significant differences. In CHB and CI the results of these comparisons were dependent on HBV genotype.CONCLUSION The lower-replication CHD and CI groups show a trend to higher quasispecies complexity than the higher-replication CHB group. The mechanisms associated with this greater complexity require elucidation.

Analysis of hepatitis B virus preS1 variability and prevalence of the rs2296651 polymorphism in a Spanish population

AIM To determine the variability/conservation of the domain of hepatitis B virus(HBV) pre S1 region that interacts with sodium-taurocholate cotransporting polypeptide(hereafter, NTCP-interacting domain) and the prevalence of the rs2296651 polymorphism(S267 F, NTCP variant) in a Spanish population. METHODS Serum samples from 246 individuals were included and divided into 3 groups: patients with chronic HBV infection(CHB)(n = 41, 73% Caucasians), patients with resolved HBV infection(n = 100, 100% Caucasians) and an HBV-uninfected control group(n = 105, 100% Caucasians). Variability/conservation of the amino acid(aa) sequences of the NTCPinteracting domain,(aa 2-48 in viral genotype D) and a highly conserved pre S1 domain associated with virion morphogenesis(aa 92-103 in viral genotype D) were analyzed by next-generation sequencing and compared in 18 CHB patients with viremia > 4 log IU/mL. The rs2296651 polymorphism was determined in all individuals in all 3 groups using an in-house real-time PCR melting curve analysis.RESULTS The HBV pre S1 NTCP-interacting domain showed a high degree of conservation among the examined viral genomes especially between aa 9 and 21(in the genotype D consensus sequence). As compared with the virion morphogenesis domain, the NTCPinteracting domain had a smaller proportion of HBV genotype-unrelated changes comprising > 1% of the quasispecies(25.5% vs 31.8%), but a larger proportion of genotype-associated viral polymorphisms(34% vs 27.3%), according to consensus sequences from Gen Bank patterns of HBV genotypes A to H. Variation/conservation in both domains depended on viral genotype, with genotype C being the most highly conserved and genotype E the most variable(limited finding, only 2 genotype E included). Of note, proline residues were highly conserved in both domains, and serine residues showed changes only to threonine or tyrosine in the virion morphogenesis domain. The rs2296651 polymorphism was not detected in any participant.CONCLUSION In our CHB population, the NTCP-interacting domain was highly conserved, particularly the proline residues and essential amino acids related with the NTCP interaction, and the prevalence of rs2296651 was low/null.

Quasispecies dynamics in main core epitopes of hepatitis B virus by ultra-deep-pyrosequencing

AIM:To investigate the variability of the main immunodominant motifs of hepatitis B virus(HBV) core gene by ultra-deep-pyrosequencing(UDPS).METHODS:Four samples(2 genotype A and 2 genotype D) from 4 treatment-na ve patients were assessed for baseline variability.Two additional samples from one patient(patient 4,genotype D) were selected for analysis:one sample corresponded to a 36-mo treatment-free period from baseline and the other to the time of viral breakthrough after 18 mo of lamivudine treatment.The HBV region analyzed covered amino acids 40 to 95 of the core gene,and included the two main epitopic regions,Th50-69 and B74-84.UDPS was carried out in the Genome Sequencer FLX system(454 Life Sciences,Roche).After computer filtering of UDPS data based on a Poisson statistical model,122 813 sequences were analyzed.The most conserved position detected by UDPS was analyzed by site-directed mutagenesis and evaluated in cell culture.RESULTS:Positions with highest variability rates were mainly located in the main core epitopes,confirming their role as immune-stimulating regions.In addition,the distribution of variability showed a relationship with HBV genotype.Patient 1(genotype A) presented the lowest variability rates and patient 2(genotype A) had 3 codons with variability higher than 1%.Patient 3 and 4(both genotype D) presented 5 and 8 codons with variability higher than 1%,respectively.The median baseline frequencies showed that genotype A samples had higher variability in epitopic positions than in the other positions analyzed,approaching significance(P = 0.07,sample 1 and P = 0.05,sample 2).In contrast,there were no significant differences in variability between the epitopic and other positions in genotype D cases.Interestingly,patient 1 presented a completely mutated motif from amino acid 64 to 67(E 64 LMT 67),which is commonly recognized by T helper cells.Additionally,the variability observed in all 4 patients was particularly associated with the E 64 LMT 67 motif.Codons 78 and 79 were highly conserved in all samples,in keeping with their involvement in the interaction between the HBV virion capsid and the surface antigens(HBsAg).Of note,codon 76 was even more conserved than codons 78 and 79,suggesting a possible role in HBsAg interactions or even in hepatitis B e antigen conformation.Sequential analysis of samples from patient 4(genotype D) illustrated the dynamism of the HBV quasispecies,with strong selection of one minor baseline variant coinciding with a decrease in core variability during the treatment-free and lamivudinetreated period.The drop in variability seemed to result from a "steady state" situation of the HBV quasispecies after selection of the variant with greatest fitness.CONCLUSION:Host immune pressure seems to be the main cause of HBV core evolution.UDPS analysis is a useful technique for studying viral quasispecies.

Molecular epidemiology and putative origin of hepatitis C virus in random volunteers from Argentina

AIM:To study the subtype prevalence and the phylogenetic relatedness of hepatitis C virus(HCV)sequences obtained from the Argentine general population,a large cohort of individuals was analyzed.METHODS:Healthy Argentinian volunteers(n=6251)from 12 provinces representing all geographical regions of the country were studied.All parents or legal guardians of individuals younger than 18 years provided informed written consent for participation.The corresponding written permission from all municipal authorities was obtained from each city or town where subjects were to be included.HCV RNA reverse transcription-polymerase chain reaction products were sequenced and phylogenetically analyzed.The 5’untranslated region(5’UTR)was used for RNA detection and initial genotype classification.The NS5B polymerase region,encompassing nt 8262-8610,was used for subtyping.RESULTS:An unexpectedly low prevalence of HCV infection in the general population(0.32%)was observed.Our data contrasted with previous studies that reported rates ranging from 1.5%to 2.5%,mainly performed in selected populations of blood donors or vulnerable groups.The latter values are in keeping with the prevalence reported by the 2007 Argentinian HCV Consensus(approximately 2%).HCV subtypes weredistributed as follows:1a(25%),1b(25%),2c(25%),3a(5%),and 2j(5%).Two isolates ascribed either to genotype 1(5%)or to genotype 3(5%)by 5’UTR phylogenetic analysis could not be subtyped.Subtype 1a sequences comprised a highly homogeneous population and clustered with United States sequences.Genotype1b sequences represented a heterogeneous population,suggesting that this genotype might have been introduced from different sources.Most subtype 2c sequences clustered close to the 2c reported from Italy and Southern France.CONCLUSION:HCV has a low prevalence of 0.32%in the studied general population of Argentina.The pattern of HCV introduction and transmission in Argentina appears to be a consequence of multiple events and different for each subtype.

Cross-sectional evaluation of circulating hepatitis B virus RNA and DNA: Different quasispecies?

BACKGROUND Different forms of pregenomic and other hepatitis B virus(HBV)RNA have been detected in patients’sera.These circulating HBV-RNAs may be useful for monitoring covalently closed circular DNA activity,and predicting hepatitis B eantigen seroconversion or viral rebound after nucleos(t)ide analog cessation.Data on serum HBV-RNA quasispecies,however,is scarce.It is therefore important to develop methodologies to thoroughly analyze this quasispecies,ensuring the elimination of any residual HBV-DNA.Studying circulating HBV-RNA quasispecies may facilitate achieving functional cure of HBV infection.AIM To establish a next-generation sequencing(NGS)methodology for analyzing serum HBV-RNA and comparing it with DNA quasispecies.METHODS Thirteen untreated chronic hepatitis B patients,showing different HBV-genotypes and degrees of severity of liver disease were enrolled in the study and a serum sample with HBV-DNA>5 Log10 IU/mL and HBV-RNA>4 Log10 copies/mL was taken from each patient.HBV-RNA was treated with DNAse I to remove any residual DNA,and the region between nucleotides(nt)1255-1611 was amplified using a 3-nested polymerase chain reaction protocol,and analyzed with NGS.Variability/conservation and complexity was compared between HBV-DNA and RNA quasispecies.RESULTS No HBV-DNA contamination was detected in cDNA samples from HBV-RNA quasispecies.HBV quasispecies complexity showed heterogeneous behavior among patients.The Rare Haplotype Load at 1%was greater in DNA than in RNA quasispecies,with no statistically significant differences(P=0.1641).Regarding conservation,information content was equal in RNA and DNA quasispecies in most nt positions[218/357(61.06%)].In 102 of the remaining 139(73.38%),HBV-RNA showed slightly higher variability.Sliding window analysis identified 4 hyper-conserved sequence fragments in each quasispecies,3 of them coincided between the 2 quasispecies:nts 1258-1286,1545-1573 and 1575-1604.The 2 hyper-variable sequence fragments also coincided:nts 1311-1344 and 1461-1485.Sequences between nts 1519-1543 and 1559-1587 were only hyper-conserved in HBV-DNA and RNA,respectively.CONCLUSION Our methodology allowed analyzing HBV-RNA quasispecies complexity and conservation without interference from HBV-DNA.Thanks to this,we have been able to compare both quasispecies in the present study.

Conservation and variability of hepatitis B core at different chronic hepatitis stages

BACKGROUND Since it is currently not possible to eradicate hepatitis B virus(HBV)infection with existing treatments,research continues to uncover new therapeutic strategies.HBV core protein,encoded by the HBV core gene(HBC),intervenes in both structural and functional processes,and is a key protein in the HBV life cycle.For this reason,both the protein and the gene could be valuable targets for new therapeutic and diagnostic strategies.Moreover,alterations in the protein sequence could serve as potential markers of disease progression.AIM To detect,by next-generation sequencing,HBC hyper-conserved regions that could potentially be prognostic factors and targets for new therapies.METHODS Thirty-eight of 45 patients with chronic HBV initially selected were included and grouped according to liver disease stage[chronic hepatitis B infection without liver damage(CHB,n=16),liver cirrhosis(LC,n=5),and hepatocellular carcinoma(HCC,n=17)].HBV DNA was extracted from patients’plasma.A region between nucleotide(nt)1863 and 2483,which includes HBC,was amplified and analyzed by next-generation sequencing(Illumina Mi Seq platform).Sequences were genotyped by distance-based discriminant analysis.General and intergroup nt and amino acid(aa)conservation was determined by sliding window analysis.The presence of nt insertion and deletions and/or aa substitutions in the different groups was determined by aligning the sequences with genotype-specific consensus sequences.RESULTS Three nt(nt 1900-1929,2249-2284,2364-2398)and 2 aa(aa 117-120,159-167)hyper-conserved regions were shared by all the clinical groups.All groups showed a similar pattern of conservation,except for five nt regions(nt 1946-1992,2060-2095,2145-2175,2230-2250,2270-2293)and one aa region(aa 140-160),where CHB and LC,respectively,were less conserved(P<0.05).Some group-specific conserved regions were also observed at both nt(2306-2334 in CHB and 1935-1976 and 2402-2435 in LC)and aa(between aa 98-103 in CHB and 28-30 and 51-54 in LC)levels.No differences in insertion and deletions frequencies were observed.An aa substitution(P79 Q)was observed in the HCC group with a median(interquartile range)frequency of 15.82(0-78.88)vs 0(0-0)in the other groups(P<0.05 vs CHB group).CONCLUSION The differentially conserved HBC and HBV core protein regions and the P79 Q substitution could be involved in disease progression.The hyper-conserved regions detected could be targets for future therapeutic and diagnostic strategies.

New real-time-PCR method to identify single point mutations in hepatitis C virus

AIM To develop a fast, low-cost diagnostic strategy to identify single point mutations in highly variable genomes such as hepatitis C virus(HCV).METHODS In patients with HCV infection, resistance-associated amino acid substitutions within the viral quasispecies prior to therapy can confer decreased susceptibility to direct-acting antiviral agents and lead to treatment failure and virological relapse. One such naturally occurring mutation is the Q80 K substitution in the HCV-NS3 protease gene, which confers resistance to PI inhibitors, particularly simeprevir. Low-cost, highly sensitive techniques enabling routine detection of these single point mutations would be useful to identify patients at a risk of treatment failure. Light Cycler methods, based on real-time PCR with sequencespecific probe hybridization, have been implemented in most diagnostic laboratories. However, this technique cannot identify single point mutations in highly variable genetic environments, such as the HCV genome. To circumvent this problem, we developed a new method to homogenize all nucleotides present in a region except the point mutation of interest. RESULTS Using nucleotide-specific probes Q, K, and R substitutions at position 80 were clearly identified at a sensitivity of 10%(mutations present at a frequency of at least 10% were detected). The technique was successfully applied to identify the Q80 K substitution in 240 HCV G1 serum samples, with performance comparable to that of direct Sanger sequencing, the current standard procedure for this purpose. The new method was then validated in a Catalonian population of 202 HCV G1-infected individuals. Q80 K was detected in 14.6% of G1 a patients and 0% of G1 b in our setting. CONCLUSION A fast, low-cost diagnostic strategy based on real-time PCR and fluorescence resonance energy transfer probe melting curve analysis has been successfully developed to identify single point mutations in highly variable genomes such as hepatitis C virus. This technique can be adapted to detect any single point mutation in highly variable genomes.