Viral quasispecies diversity and evolution : a Bioinformatics molecular approach

El grup de hepatitis virals del Vall d'Hebron Institut de Recerca (VHIR) de Barcelona, en els darrers 10 anys, ha estat desenvolupant solucions metodològiques experimentals i computacionals per a l'estudi de poblacions complexes de virus (quasispecies) mitjançant l'aplicació de les tècniques de seqüenciació de nova generació (NGS). Aquest llibre consisteix en una selecció de treballs empírics sobre les quasispecies virals. En oferir aquest format obert de publicació, l'objectiu és que, per una banda, pugui ser una eina útil per a tots aquells investigadors interessats en aquest camp i, d'altra banda, divulgar aquesta àrea de coneixement a tota la comunitat científica que, sense ser necessàriament experta, vulgui conèixer amb més detall l'evolució i diversitat dels virus. En els tres primers treballs s'aprofundeix en l'ús, interpretació i utilitat d'índexs de biodiversitat, alguns específics per a poblacions genètiques i d'altres importats del camp de l'ecologia. La segona part posa de manifest algunes limitacions en aquests índexs de diversitat i aborda el desenvolupament d'eines integradores que proporcionen una interpretació més directa en termes biològics i clínics. Les seccions prèvies als sis treballs esmentats, situen al lector en el context en què es realitzen els desenvolupaments i expliquen la necessitat i utilitat. El llibre es tanca amb una secció que recull les observacions i conclusions generals dels treballs, i amb una altra que reflexiona sobre les limitacions que comporta l'estudi de sistemes complexos i dinàmics com les quasispecies virals

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

Altres ajuts: Cofinanced by the European Regional Development Fund (ERDF); and the Gilead Fellowship Program, No. GLD14-00296.To determine the variability/conservation of the domain of hepatitis B virus (HBV) preS1 region that interacts with sodium-taurocholate cotransporting polypeptide (hereafter, NTCP-interacting domain) and the prevalence of the rs2296651 polymorphism (S267F, NTCP variant) in a Spanish population. 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 NTCP-interacting domain, (aa 2-48 in viral genotype D) and a highly conserved preS1 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. The HBV preS1 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 NTCP-interacting 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 GenBank 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. 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

Study of Quasispecies Complexity and Liver Damage Progression after Liver Transplantation in Hepatitis C Virus Infected Patients

Fibrosis; Hepatitis C virus; Viral loadFibrosis; Virus de la hepatitis C; Carga viralFibrosi; Virus de l'hepatitis C; Càrrega viralCirrhosis derived from chronic hepatitis C virus (HCV) infection is still a common indication for liver transplantation (LT). Reinfection of the engrafted liver is universal in patients with detectable viral RNA at the time of transplant and causes fast progression of cirrhosis (within 5 years) in around one-third of these patients. To prevent damage to the liver graft, effective direct-acting antiviral (DAA) therapy is required as soon as possible. However, because of post-LT clinical instability, it is difficult to determine the optimal time to start DAAs with a low risk of complications. Evaluate changes in quasispecies complexity following LT and seek a predictive index of fast liver damage progression to determine the timing of DAA initiation. HCV genomes isolated from pre-LT and 15-day post-LT serum samples of ten patients, who underwent orthotopic LT, were quantified and sequenced using a next-generation sequencing platform. Sequence alignments, phylogenetic trees, quasispecies complexity measures, biostatistics analyses, adjusted R2 values, and analysis of variance (ANOVA) were carried out. Three different patterns of reinfection were observed (viral bottlenecking, conserved pre-LT population, and mixed populations), suggesting that bottlenecking or homogenization of the viral population is not a generalized effect after liver graft reinfection. None of the quasispecies complexity measures predicted the future degree of liver damage. Higher and more uniform viral load (VL) values were observed in all pre-LT samples, but values were more dispersed in post-LT samples. However, VL increased significantly from the pre-LT to 15-day post-LT samples in patients with advanced fibrosis at 1-year post-LT, suggesting that a VL increase on day 15 may be a predictor of fast liver fibrosis progression. HCV kinetics after LT differ between patients and are not fibrosis-dependent. Higher VL at day 15 post-LT versus pre-LT samples may predict fast liver fibrosis progression.This study was supported by grants from Instituto de Salud Carlos III cofinanced by the European Regional Development Fund (ERDF) with grant numbers PI19/00533, PI19/00301, Clinical Trial Gov. Identifier: NCT01707849, and from Centro para el Desarrollo Tecnológico Industrial-CDTI of the Spanish Ministry of Economics and Competitiveness (MINECO) grant number, IDI-20200297. C.P. is supported by the Miguel Servet program of Instituto de Salud Carlos III, grant CP14/00121, cofinanced by the ERDF. Astellas Pharma Inc and Novartis Pharma also provided funding for the study, but these companies had had no role in the study design, data collection or analysis, decision to publish, or preparation of the manuscript

Inspecting the Ribozyme Region of Hepatitis Delta Virus Genotype 1: Conservation and Variability

Gene silencing; Quasispecies; RibozymeSilenciament gènic; Quasiespècie; RibozimaSilenciamiento de genes; Cuasiespecies; RibozimaThe hepatitis delta virus (HDV) genome has an autocatalytic region called the ribozyme, which is essential for viral replication. The aim of this study was to use next-generation sequencing (NGS) to analyze the ribozyme quasispecies (QS) in order to study its evolution and identify highly conserved regions potentially suitable for a gene-silencing strategy. HDV RNA was extracted from 2 longitudinal samples of chronic HDV patients and the ribozyme (nucleotide, nt 688–771) was analyzed using NGS. QS conservation, variability and genetic distance were analyzed. Mutations were identified by aligning sequences with their specific genotype consensus. The main relevant mutations were tested in vitro. The ribozyme was conserved overall, with a hyper-conserved region between nt 715–745. No difference in QS was observed over time. The most variable region was between nt 739–769. Thirteen mutations were observed, with three showing a higher frequency: T23C, T69C and C64 deletion. This last strongly reduced HDV replication by more than 1 log in vitro. HDV Ribozyme QS was generally highly conserved and was maintained during follow-up. The most conserved portion may be a valuable target for a gene-silencing strategy. The presence of the C64 deletion may strongly impair viral replication, as it is a potential mechanism of viral persistence.This research was funded by Institute of Health Carlos III, grant number PI20/01692 and co-financed by the European Regional Development Fund (ERDF), and by the European Regional Development Fund (ERDF)- Ministry of Economy, Industry and Competitiveness, grantRTI2018-101936-B-I00

Ultra Deep Sequencing of Circulating Cell-Free DNA as a Potential Tool for Hepatocellular Carcinoma Management

Biomarkers; cfDNA; Liquid biopsyBiomarcadores; cfDNA; Biopsia líquidaBiomarcadors; cfDNA; Biòpsia líquidaBackground: Cell-free DNA (cfDNA) concentrations have been described to be inversely correlated with prognosis in cancer. Mutations in HCC-associated driver genes in cfDNA have been reported, but their relation with patient’s outcome has not been described. Our aim was to elucidate whether mutations found in cfDNA could be representative from those present in HCC tissue, providing the rationale to use the cfDNA to monitor HCC. Methods: Tumoral tissue, paired nontumor adjacent tissue and blood samples were collected from 30 HCC patients undergoing curative therapies. Deep sequencing targeting HCC driver genes was performed. Results: Patients with more than 2 ng/µL of cfDNA at diagnosis had higher mortality (mean OS 24.6 vs. 31.87 months, p = 0.01) (AUC = 0.782). Subjects who died during follow-up, had a significantly higher number of mutated genes (p = 0.015) and number of mutations (p = 0.015) on cfDNA. Number of mutated genes (p = 0.001), detected mutations (p = 0.001) in cfDNA and ratio (number of mutations/cfDNA) (p = 0.003) were significantly associated with recurrence. However, patients with a ratio (number of mutations/cfDNA) above 6 (long-rank p = 0.0003) presented a higher risk of recurrence than those with a ratio under 6. Detection of more than four mutations in cfDNA correlated with higher risk of death (long-rank p = 0.042). Conclusions: In summary, cfDNA and detection of prevalent HCC mutations could have prognostic implications in early-stage HCC patientsThis research was funded by Instituto de Salud Carlos III (ISCIII) (PI18/00961) and (PI21/00714) to B.M. PI19/00301 by Instituto de Salud Carlos III (ISCIII) and Centro para el Desarrollo Tecnológico Industrial (CDTI) from the Spanish Ministry of Economy and Business, grant number IDI-20200297 to J.Q., E.V.-A. is a recipient of a predoctoral fellowship from Instituto de Salud Carlos III (ISCIII) (FI18/00027)

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

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. To detect, by next-generation sequencing, HBC hyper-conserved regions that could potentially be prognostic factors and targets for new therapies. 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 MiSeq 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. 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 (P79Q) 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). The differentially conserved HBC and HBV core protein regions and the P79Q substitution could be involved in disease progression. The hyper-conserved regions detected could be targets for future therapeutic and diagnostic strategies

Sophisticated viral quasispecies with a genotype-related pattern of mutations in the hepatitis B X gene of HBeAg-ve chronically infected patients

Haplotipos; Enfermedades infecciosas; MutaciónHaplotips; Enfermetats infeccioses; MutacióHaplotypes; Infectious diseases; MutationPatients with HBeAg-negative chronic infection (CI) have not been extensively studied because of low viremia. The HBx protein, encoded by HBX, has a key role in viral replication. Here, we analyzed the viral quasispecies at the 5′ end of HBX in CI patients and compared it with that of patients in other clinical stages. Fifty-eight HBeAg-negative patients were included: 16 CI, 19 chronic hepatitis B, 16 hepatocellular carcinoma and 6 liver cirrhosis. Quasispecies complexity and conservation were determined in the region between nucleotides 1255 and 1611. Amino acid changes detected were tested in vitro. CI patients showed higher complexity in terms of mutation frequency and nucleotide diversity and higher quasispecies conservation (p < 0.05). A genotype D-specific pattern of mutations (A12S/P33S/P46S/T36D-G) was identified in CI (median frequency, 81.7%), which determined a reduction in HBV DNA release of up to 1.5 log in vitro. CI patients showed a more complex and conserved viral quasispecies than the other groups. The genotype-specific pattern of mutations could partially explain the low viremia observed in these patients.This study was supported by grants from the Institute of Health Carlos III (grant PI15/00856 and PI18/01436) and co-financed by the European Regional Development Fund (ERDF)

Baseline hepatitis C virus resistance-associated substitutions present at frequencies lower than 15% may be clinically significant

Controversy is ongoing about whether a minority mutant present at frequencies below 15% may be clinically relevant and should be considered to guide treatment. Resistance-associated substitution (RAS) studies were performed in patients before and at failure of antiviral treatments using Next-generation hepatitis C virus (HCV) sequencing (NGS). We have found two patients with genotype 1a infection having RAS in 3.5%-7.1% of the viral population at baseline that were selected during ledipasvir + sofosbuvir treatment. Co-selection of RAS located in a region not directly affected by the antiviral treatment also occurred. This observation calls into question, the recommendations to guide RAS-based direct-acting antiviral (DAA) treatment only when RAS are present in >15% of the sequences generated. Our results suggests that RAS study should include all three HCV DAA target proteins and minority mutants should be considered as clinically relevant

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

To detect hyper-conserved regions in the hepatitis B virus (HBV) X gene (HBX) 5' region that could be candidates for gene therapy. The study included 27 chronic hepatitis B treatment-naive patients in various clinical stages (from chronic infection to cirrhosis and hepatocellular carcinoma, both HBeAg-negative and HBeAg-positive), and infected with HBV genotypes A-F and H. In a serum sample from each patient with viremia > 3.5 log IU/mL, the HBX 5' end region [nucleotide (nt) 1255-1611] was PCR-amplified 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. 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. 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

Quantifying In-Host Quasispecies Evolution

What takes decades, centuries or millennia to happen with a natural ecosystem, it takes only days, weeks or months with a replicating viral quasispecies in a host, especially when under treatment. Some methods to quantify the evolution of a quasispecies are introduced and discussed, along with simple simulated examples to help in the interpretation and understanding of the results. The proposed methods treat the molecules in a quasispecies as individuals of competing species in an ecosystem, where the haplotypes are the competing species, and the ecosystem is the quasispecies in a host, and the evolution of the system is quantified by monitoring changes in haplotype frequencies. The correlation between the proposed indices is also discussed, and the R code used to generate the simulations, the data and the plots is provided. The virtues of the proposed indices are finally shown on a clinical case