Recombination in the genesis and evolution of hepatitis B virus genotypes

Hepatitis B virus (HBV) infection is widely distributed in both human and ape populations throughout the world and is a major cause of human morbidity and mortality. HBV variants are currently classified into the human genotypes A to H and species-associated chimpanzee and gibbon/orangutan groups. To examine the role of recombination in the evolution of HBV, large-scale data retrieval and automated phylogenetic analysis (TreeOrder scanning) were carried out on all available published complete genome sequences of HBV. We detected a total of 24 phylogenetically independent potential recombinants (different genotype combinations or distinct breakpoints), eight of which were previously undescribed. Instances of intergenotype recombination were observed in all human and ape HBV variants, including evidence for a novel gibbon/genotype C recombinant among HBV variants from Vietnam. By recording sequence positions in trees generated from sequential fragments across the genome, violations of phylogeny between trees also provided evidence for frequent intragenotype recombination between members of genotypes A, D, F/H, and gibbon variants but not in B, C, or the Asian B/C recombinant group. In many cases, favored positions for both inter- and intragenotype recombination matched positions of phylogenetic reorganization between the human and ape genotypes, such as the end of the surface gene and the core gene, where sequence relationships between genotypes changed in the TreeOrder scan. These findings provide evidence for the occurrence of past, extensive recombination events in the evolutionary history of the currently classified genotypes of HBV and potentially in changes in its global epidemiology and associations with human disease.     

A revised evolutionary history of hepatitis B virus (HBV)

Previous studies of the evolutionary history of hepatitis B virus (HBV) have been compromised by intergenotype recombination and complex patterns of nucleotide substitution, perhaps caused by differential selection pressures. We examined the phylogenetic distribution of recombination events among human HBV genotypes and found that genotypes A plus D, and genotypes B plus C, had distinct patterns of recombination suggesting differing epidemiological relationships among them. By analyzing the nonoverlapping regions of the viral genome we found strong bootstrap support for some intergenotypic groupings, with evidence of a division between human genotypes A–E from the viruses sampled from apes and human genotype F. However, the earliest events in the divergence of HBV remain uncertain. These uncertainties could not be explained by differential selection pressures, as the ratio of nonsynonymous-to-synonymous substitutions (d _N/d _S) did not vary extensively among lineages and there is no strong evidence for positive selection across the whole tree. Finally, we provide a new estimate of the mean substitution rate in HBV, 4.2 × 10⁻⁵, which suggests that divergence of HBV in humans and apes has occurred only in the last 6000 years.     

Reconstructing the Complex Evolutionary History of Hepatitis B Virus

A detailed analysis of the evolutionary history of hepatitis B virus (HBV) was undertaken using 39 mammalian hepadnaviruses for which complete genome sequences were available, including representatives of all six human genotypes, as well as a large sample of small S gene sequences. Phylogenetic trees of these data were ambiguous, supporting no single place of origin for HBV, and depended heavily on the underlying model of DNA substitution. In some instances genotype F, predominant in the Americas, was the first to diverge, suggesting that the virus arose in the New World. In other trees, however, sequences from genotype B, prevalent in East Asia, were the most divergent. An attempt was also made to determine the rate of nucleotide substitution in the C open reading frame and then to date the origin of HBV. However, no relationship between time and number of substitutions was found in two independent data sets, indicating that a reliable molecular clock does not exist for these data. Both the pattern and the rate of nucleotide substitution are therefore complex phenomena in HBV and hinder any attempt to reconstruct the past spread of this virus. Received: 5 December 1998 / Accepted: 23 February 1999     

Hepatitis B virus of genotype B with or without recombination with genotype C over the precore region plus the core gene

The entire nucleotide sequences of 70 hepatitis B virus (HBV) isolates of genotype B (HBV/B), including 38 newly determined and 32 retrieved from the international DNA database (DDBJ/EMBL/GenBank), were compared phylogenetically. Two subgroups of HBV/B were identified based on sequence divergence in the precore region plus the core gene, one with the recombination with genotype C and the other without it. The analysis over the entire genome of HBV/B by the SimPlot program located the recombination with genotype C in the precore region plus the core gene spanning nucleotide positions from 1740 to 1838 to 2443 to 2485. Within this genomic area, HBV/B strains with the recombination had higher nucleotide and amino acid homology to genotype C than those without the recombination (96.9 versus 91.1% in nucleotides and 97.0 versus 92.9% in amino acids). There were 29 HBV/B strains without the recombination, and they were all recovered from carriers in Japan. The remaining 41 HBV/B isolates having the recombination with genotype C were from carriers in China (12 strains), Hong Kong (3 strains), Indonesia (4 strains), Japan (3 strains), Taiwan (4 strains), Thailand (3 strains), and Vietnam (12 strains). Due to the frequency of the distribution of HBV/B without the recombination (29 of 32 isolates, or 91%) and the fact that it was exclusive to Japan, it was provisionally classified into the Bj (j standing for Japan) subgroup, and HBV/B with the recombination was classified into the Ba (a for Asia) subgroup. Virological differences between HBV/Bj and HBV/Ba may be reflected in the severity of clinical disease in the patients infected with HBV of genotype B, which seems to be under strong geographic influences in Asia.     

A Rare HBV Subgenotype D4 with Unique Genomic Signatures Identified in North-Eastern India –An Emerging Clinical Challenge?

Background/Aims

HBV has been classified into ten genotypes (A–J) and multiple subgenotypes, some of which strongly influence disease outcome and their distribution also correlate with human migration. HBV infection is highly prevalent in India and its diverse population provides an excellent opportunity to study the distinctiveness of HBV, its evolution and disease biology in variegated ethnic groups. The North-East India, having international frontiers on three sides, is one of the most ethnically and linguistically diverse region of the country. Given the paucity of information on molecular epidemiology of HBV in this region, the study aimed to carry out an in-depth genetic characterization of HBV prevailing in North-East state of Tripura.

Methods

From sera of chronically HBV infected patients biochemical/serological tests, HBV DNA quantification, PCR-amplification, sequencing of PreS/S or full-length HBV genomes were done. HBV genotype/subgenotype determination and sequence variability were assessed by MEGA5-software. The evolutionary divergence times of different HBV subgenotypes were estimated by DNAMLK/PHYLIP program while jpHMM method was used to detect any recombination event in HBV genomes.

Results

HBV genotypes D (89.5%), C (6.6%) and A (3.9%) were detected among chronic carriers. While all HBV/A and HBV/C isolates belonged to subgenotype-A1 and C1 respectively, five subgenotypes of HBV/D (D1–D5) were identified including the first detection of rare D4. These non-recombinant Indian D4 (IndD4) formed a distinct phylogenetic clade, had 2.7% nucleotide divergence and recent evolutionary radiation than other global D4. Ten unique amino acids and 9 novel nucleotide substitutions were identified as IndD4 signatures. All IndD4 carried T¹²⁰ and R¹²⁹ in ORF-S that may cause immune/vaccine/diagnostic escape and N¹²⁸ in ORF-P, implicated as compensatory Lamivudine resistance mutation.

Conclusions

IndD4 has potential to undermine vaccination programs or anti-viral therapy and its introduction to North-East India is believed to be linked with the settlement of ancient Tibeto-Burman migrants from East-Asia.     

Cytoplasmic and genomic effects on meiotic pairing in brassica hybrids and allotetraploids from pair crosses of three cultivated diploids

Interspecific hybridization and allopolyploidization contribute to the origin of many important crops. Synthetic Brassica is a widely used model for the study of genetic recombination and "fixed heterosis" in allopolyploids. To investigate the effects of the cytoplasm and genome combinations on meiotic recombination, we produced digenomic diploid and triploid hybrids and trigenomic triploid hybrids from the reciprocal crosses of three Brassica diploids (B. rapa, AA; B. nigra, BB; B. oleracea, CC). The chromosomes in the resultant hybrids were doubled to obtain three allotetraploids (B. juncea, AA.BB; B. napus, AA.CC; B. carinata, BB.CC). Intra- and intergenomic chromosome pairings in these hybrids were quantified using genomic in situ hybridization and BAC-FISH. The level of intra- and intergenomic pairings varied significantly, depending on the genome combinations and the cytoplasmic background and/or their interaction. The extent of intragenomic pairing was less than that of intergenomic pairing within each genome. The extent of pairing variations within the B genome was less than that within the A and C genomes, each of which had a similar extent of pairing. Synthetic allotetraploids exhibited nondiploidized meiotic behavior, and their chromosomal instabilities were correlated with the relationship of the genomes and cytoplasmic background. Our results highlight the specific roles of the cytoplasm and genome to the chromosomal behaviors of hybrids and allopolyploids.     

Comparative Analysis of CpG Islands among HBV Genotypes

DNA methylation is being increasingly recognized to play a role in regulation of hepatitis B virus (HBV) gene expression. The aim of this study was to compare the CpG island distribution among different HBV genotypes. We analyzed 176 full-length HBV genomic sequences obtained from the GenBank database, belonging to genotypes A through J, to identify the CpG islands in the HBV genomes. Our results showed that while 79 out of 176 sequences contained three conventional CpG islands (I–III) as previously described, 83 HBV sequences harbored only two of the three known islands. Novel CpG islands were identified in the remaining 14 HBV isolates and named as CpG island IV, V, and VI. Among the eight known HBV genotypes and two putative genotypes, while HBV genomes containing three CpG islands were predominant in genotypes A, B, D, E, and I; genotypes C, F, G, and H tended to contain only two CpG islands (II and III). In conclusion, the CpG islands, which are potential targets for DNA methylation mediated by the host functions, differ among HBV genotypes, and these genotype-specific differences in CpG island distribution could provide new insights into the understanding of epigenetic regulation of HBV gene expression and hepatitis B disease outcome.     

Genotype-specific genomic markers associated with primary hepatomas, based on complete genomic sequencing of hepatitis B virus

We aimed to identify genomic markers in hepatitis B virus (HBV) that are associated with hepatocellular carcinoma (HCC) development by comparing the complete genomic sequences of HBVs among patients with HCC and those without. One hundred patients with HBV-related HCC and 100 age-matched HBV-infected non-HCC patients (controls) were studied. HBV DNA from serum was directly sequenced to study the whole viral genome. Data mining and rule learning were employed to develop diagnostic algorithms. An independent cohort of 132 cases (43 HCC and 89 non-HCC) was used to validate the accuracy of these algorithms. Among the 100 cases of HCC, 37 had genotype B (all subgenotype Ba) and 63 had genotype C (16 subgenotype Ce and 47 subgenotype Cs) HBV infection. In the control group, 51 had genotype B and 49 had genotype C (10 subgenotype Ce and 39 subgenotype Cs) HBV infection. Genomic algorithms associated with HCC were derived based on genotype/subgenotype-specific mutations. In genotype B HBV, mutations C1165T, A1762T and G1764A, T2712C/A/G, and A/T2525C were associated with HCC. HCC-related mutations T31C, T53C, and A1499G were associated with HBV subgenotype Ce, and mutations G1613A, G1899A, T2170C/G, and T2441C were associated with HBV subgenotype Cs. Amino acid changes caused by these mutations were found in the X, envelope, and precore/core regions in association with HBV genotype B, Ce, and Cs, respectively. In conclusion, infections with different genotypes of HBV (B, Ce, and Cs) carry different genomic markers for HCC at different parts of the HBV genome. Different HBV genotypes may have different virologic mechanisms of hepatocarcinogenesis.     

Adaptive evolution of UGT2B17 copy-number variation

The human UGT2B17 gene varies in copy number from zero to two per individual and also differs in mean number between populations from Africa, Europe, and East Asia. We show that such a high degree of geographical variation is unusual and investigate its evolutionary history. This required first reinterpreting the reference sequence in this region of the genome, which is misassembled from the two different alleles separated by an artifactual gap. A corrected assembly identifies the polymorphism as a 117 kb deletion arising by nonallelic homologous recombination between ~4.9 kb segmental duplications and allows the deletion breakpoint to be identified. We resequenced ~12 kb of DNA spanning the breakpoint in 91 humans from three HapMap and one extended HapMap populations and one chimpanzee. Diversity was unusually high and the time to the most recent common ancestor was estimated at ~2.4 or ~3.0 million years by two different methods, with evidence of balancing selection in Europe. In contrast, diversity was low in East Asia where a single haplotype predominated, suggesting positive selection for the deletion in this part of the world.     

Specific mutations of basal core promoter are associated with chronic liver disease in hepatitis B virus subgenotype D1 prevalent in Turkey

The role of hepatitis B virus (HBV) genetics in the clinical manifestations of infection is being increasingly recognized. Genotype D is one of eight currently recognized major HBV genotypes. The virus is ubiquitous worldwide, but shows different features in different regions. One hundred and ninety‐eight patients with chronic HBV infection were enrolled in this study, 38 of whom had been diagnosed with cirrhosis of the liver and/or hepatocellular carcinoma. HBV DNA was isolated from the patients' blood samples and the entire genome and/or the basal core promoter/core promoter region sequenced. Phylogenetic analysis of the complete genomes revealed that subgenotype D1 is the most prevalent subgenotype in Turkey, but there was no definite phylogenetic grouping according to geography for isolates from different regions within Turkey, or for isolates in Turkey relative to other parts of the world. Turkish isolates tended to be genetically similar to European and central Asian isolates. Overall, HBV‐infection in Turkey appears to be characterized by early HBeAg seroconversion, a high incidence of the A1896 core promoter mutation and a small viral load. Genotype D characteristic mutations A1757 and T1764/G1766 were found in the BCP region. T1773 was associated with T1764/G1766 and a larger viral load. In conclusion, infection with HBV genotype D in Turkey has a similar clinical outcome to that of Europe and central Asia. Genotypic mutations in genotype D may be linked with disease prognosis in Turkey, but further studies with higher sample numbers and balanced clinical groups are needed to confirm this.     

Non-hybrid offspring from matings between hemiclonal hybrid waterfrogs suggest occasional recombination between clonal genomes

The hemiclonal waterfrog Rana esculenta , a hybrid between R. ridibunda and R. lessonae , eliminates the lessonae genome from the germline and clonally transmits the ridibunda genome (hybridogenesis). Such genomes are prone to accumulate deleterious mutations, which may explain why offspring from matings between hybrids are typically inviable. Here I present field data from a population for which experimental crossings showed that some R. esculenta pairs produce viable R. ridibunda offspring. I demonstrate: (1) that R. ridibunda metamorphs are also produced and survive under natural conditions; (2) that their genotypes are consistent with combinations of clonal ridibunda genomes found in hybrids; and (3) that all R. ridibunda are female. These females possibly recombine the clonal genomes they inherited and, upon mating with syntopic R. lessonae , produce new hemiclones with novel combinations of alleles. Hence, occasional recombination between otherwise clonal ridibunda genomes seems plausible and may provide an escape from the evolutionary dead end they were proposed to be trapped in.     

Evaluation of intra-host variants of the entire hepatitis B virus genome

Genetic analysis of hepatitis B virus (HBV) frequently involves study of intra-host variants, identification of which is commonly achieved using short regions of the HBV genome. However, the use of short sequences significantly limits evaluation of genetic relatedness among HBV strains. Although analysis of HBV complete genomes using genetic cloning has been developed, its application is highly labor intensive and practiced only infrequently. We describe here a novel approach to whole genome (WG) HBV quasispecies analysis based on end-point, limiting-dilution real-time PCR (EPLD-PCR) for amplification of single HBV genome variants, and their subsequent sequencing. EPLD-PCR was used to analyze WG quasispecies from serum samples of patients (n = 38) infected with HBV genotypes A, B, C, D, E and G. Phylogenetic analysis of the EPLD-isolated HBV-WG quasispecies showed the presence of mixed genotypes, recombinant variants and sub-populations of the virus. A critical observation was that HBV-WG consensus sequences obtained by direct sequencing of PCR fragments without EPLD are genetically close, but not always identical to the major HBV variants in the intra-host population, thus indicating that consensus sequences should be judiciously used in genetic analysis. Sequence-based studies of HBV WG quasispecies should afford a more accurate assessment of HBV evolution in various clinical and epidemiological settings.     

Subtype ayw variant of hepatitis B virus. DNA primary structure analysis

The entire genome of human hepatitis B virus (HBV) occurring in Latvia was sequenced. This sequence, which is 3182 nucleotides long, was compared with the other previously published HBV genomes and was shown to share maximum homology with HBV subtype ayw DNA. The coordinates of 4 main open reading frames as well as hairpin structures are very well conserved in the two genomes. The distribution of nucleotide substitutions among different HBV genomes suggest that the open reading frames P and X can fulfil a coding function. On the basis of primary structure comparison for hepadnaviral DNAs several evolutionary conclusions can be drawn.     

Evolution of Buchlo? dactyloides based on cloning and sequencing of matK, rbcL, and cob genes from plastid and mitochondrial genomes.

Buffalograss (Buchlo? dactyloides (Nutt.) Englem), a C4 turfgrass species, is native to the Great Plains region of North America. The evolutionary implications of buffalograss are unclear. Sequencing of rbcL and matK genes from plastid and the cob gene from mitochondrial genomes was examined to elucidate buffalo grass evolution. This study is the first to report sequencing of these genes from organelle genomes in the genus Buchlo?. Comparisons of sequence data from the mitochondrial and plastid genome revealed that all genotypes contained the same cytoplasmic origin. There were some rearrangements detected in mitochondrial genome. The buffalograss genome appears to have evolved through the rearrangements of convergent subgenomic domains. Combined analyses of plastid genes suggest that the evolutionary process in Buchlo? accessions studied was monophyletic rather than polyphyletic. However, since plastid and mitochondrial genomes are generally uniparentally inherited, the evolutionary history of these genomes may not reflect the evolutionary history of the organism, especially in a species in which out-crossing is common. The sequence information obtained from this study can be used as a genome-specific marker for investigation of the buffalograss polyploidy complex and testing of the mode of plastid and mitochondrial transmission in genus Buchlo?.     

A Novel Complex A/C/G Intergenotypic Recombinant of Hepatitis B Virus Isolated in Southern China

Hepatitis B virus (HBV) genotypes and subgenotypes may vary in geographical distribution and virological features. Previous investigations, including ours, showed that HBV genotypes B and C were respectively predominant in South and North China, while genotypes A and D were infrequently detected and genotype G was not found. In this study, a novel A/C/G intergenotype was identified in patients with chronic HBV infection in Guilin, a city in southern China. Initial phylogenetic analysis based on the S gene suggested the HBV recombinant to be genotype G. However, extended genotyping based on the entire HBV genome indicated it to be an A/C/G intergenotype with a closer relation to genotype C. Breakpoint analysis using the SIMPLOT program revealed that the recombinant had a recombination with a arrangement of genotypes A, G, A and C fragments. Compared with the HBV recombinants harboring one or two genotype G fragments found in Asian countries, this Guilin recombinant was highly similar to the Vietnam (98–99%) and Long An recombinants (96–99%), but had a relatively low similarity to the Thailand one (89%). Unlike those with the typical genotype G of HBV, the patients with the Guilin recombinant were seropositive for HBeAg. Moreover, a relatively high HBV DNA viral load (>2×10⁶ IU/ml) was detected in the patients, and the analysis of viral replication capacity showed that the Guilin recombinant strains had a competent replication capacity similar to genotypes B and C strains. These findings can aid in not only the clarification of the phylogenetic origin of the HBV recombinants with the genotype G fragment found in Asian countries, but also the understanding of the virological properties of these complicated HBV recombinants.     

The variability of the hepatitis B virus genome: statistical analysis and biological implications

A statistical analysis of the nucleotide sequence variability in 14 published hepatitis B virus (HBV) genomes was carried out using parametric and nonparametric methods. A parametric statistical model revealed that the different regions of the genome differed significantly in their variability. The conclusion was supported by a nonparametric kernel-density model of the HBV genome. Genes S, C, and P, region X, the precore region, and the pre-S2/pre-S1 regions were ranked in order of increasing variability. In many instances, conserved regions of the genome identified with sequences of known function in HBV biology. However, other characterized regions (such as pre-S) showed much variability despite the involvement of their encoded peptides in specific functions. Point mutations that may result in the formation of stop codons and amino acid changes may affect the clinical picture of HBV infection and may be reflected in atypical serological patterns.      

Novel hepatitis B virus genotype a subtyping assay that distinguishes subtype Aa from Ae and its application in epidemiological studies

The eight genotypes of hepatitis B virus (HBV) have different geographical distributions, virological characteristics, and clinical manifestations. A unique subtype of HBV genotype A (HBV/A) was reported in sub-Saharan Africa, raising the possibility that patients infected with this subtype (HBV/Aa ["a" for African and Asian]) may have different clinical outcomes than other HBV/A isolates (HBV/Ae ["e" for European]). Comparison between 30 HBV/Aa and 30 HBV/Ae isolates indicated that almost all HBV/Ae isolates had G at nucleotide (nt) 1809 and C at nt 1812, whereas HBV/Aa isolates had T1809/T1812. Taking advantage of these two single nucleotide polymorphisms (SNPs), a novel subtype-specific PCR assay in the X/precore/core region was developed. This assay was combined with a restriction fragment length polymorphism assay using BglII in a different region (nt 1984 to 1989), which has a SNP distinguishing HBV/Aa from HBV/Ae, resulting in 100% specificity for the combined assay. Application of the subtyping assay using sera from 109 paid donors in the United States indicated significantly different distributions of HBV/A subtypes among races; African-Americans, Caucasians, and Hispanics had HBV/Ae, whereas Asians had mainly HBV/Aa, suggesting that the HBV/Aa isolates may have been imported by recent immigration from Asia. In conclusion, the specificity and sensitivity of the combined subtyping assay were confirmed, and its usefulness was demonstrated in a practical context.