Recombination between sequences of hepatitis B virus from different genotypes

A comparison of 25 hepatitis B virus (HBV) isolates for which complete genome sequences are available revealed two that occupied different positions in phylogenetic trees reconstructed from different open reading frames. Further analysis indicated that this incongruence was the result of recombination between viruses of different genomic and antigenic types. Both putative recombinants originated from geographic regions where multiple genotypes are known to cocirculate. A search of the sequence databases showed evidence of similar intergenotypic recombinants. These observations indicate that recombination between divergent strains may represent an important source of genetic variation in HBV.     

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     

Estimation of evolutionary distances from protein spatial structures

New equations are derived to estimate the number of amino acid substitutions per site between two homologous proteins from the root mean square (RMS) deviation between two spatial structures and from the fraction of identical residues between two sequences. The equations are based on evolutionary models, analyzing predominantly structural changes and not sequence changes. Evolution of spatial structure is treated as a diffusion in an elastic force field. Diffusion accounts for structural changes caused by amino acid substitutions, and elastic force reflects selection, which preserves protein fold. Obtained equations are supported by analysis of protein spatial structures. Received: 21 September 1995 / Accepted: 19 May 1997     

Determination of the evolutionary relationships in Rattus sensu lato (Rodentia : Muridae) using L1 (LINE-1) amplification events

We determined ∼215 bp of DNA sequence from the 3′-untranslated region (UTR) of 240 cloned L1 (LINE-1) elements isolated from 22 species of Rattus sensu lato and Rattus sensu stricto murine rodents. The sequences were sorted into different L1 subfamilies, and oligonucleotides cognate to them were hybridized to genomic DNA of various taxa. From the distribution of the L1 subfamilies in the various species, we inferred the partial phylogeny of Rattus sensu lato. The four Maxomys species comprise a well-defined clade separate from a monophyletic cluster that contains the two Leopoldamys and four Niviventer species. The Niviventer/Leopoldamys clade, in turn, shares a node with the clade that contains Berylmys, Sundamys, Bandicota, and Rattus sensu stricto. The evolutionary relationships that we deduced agree with and significantly extend the phylogeny of Rattus sensu lato established by other molecular criteria. Furthermore, the L1 amplification events scored here produced a unique phylogenetic tree, that is, in no case did a character (a given L1 amplification event) appear on more than one branch. The lack of homoplasy found in this study supports the robustness of L1 amplification events as phylogenetic markers for the study of mammalian evolution. Received: 8 November 1996 / Accepted: 11 April 1997     

免疫预防阻断乙型肝炎病毒母婴传播效果的观察

阻断乙型肝炎病毒(HBV)母婴传播是控制乙型肝炎的重大问题。为探讨免疫预防对阻断HBV母婴传播的效果及影响因素,对667例HBV表面抗原(HBsAg)阳性孕妇及其婴儿进行研究。这些孕妇按HBVe抗原(HBeAg)和HBVDNA检测结果,分为HBeAg阳性组及阴性组、DNA阳性组及阴性组;按是否于孕晚期注射乙型肝炎免疫球蛋白(HBIG),分为注射组及未注射组。婴儿于出生24h内均肌内注射HBIG100IU,并按0、1、6方案注射10μg重组酵母HBV疫苗;8～12月龄后随访婴儿,并进行HBV标志物(HBV-M)检测。667个婴儿中,20例感染HBV,免疫阻断失败率为3.0%。孕妇HBeAg阳性组免疫阻断失败率为8.7%,阴性组为0.2%,两组差异显著(P<0.001);两组婴儿对疫苗免疫应答率分别为83.0%和83.1%,无显著差异(P=0.988)。孕妇DNA阳性组免疫阻断失败率为8.1%,HBVDNA均≥6log10copies/ml。孕期注射与未注射HBIG组婴儿免疫阻断失败率分别为3.7%和2.7%,无显著差异(P=0.479);两组婴儿对疫苗免疫应答率分别为84.4%和82.4%,无显著差异(P=0.519)。孕妇HBeAg阳性注射HBIG组与未注射组的免疫阻断失败率分别为8.4%和8.9%,无显著差异(P=0.892)。孕妇HBeAg阴性注射与未注射HBIG组的免疫阻断失败率分别为0.0%和0.3%,也无显著差异(P=0.538)。11例免疫阻断失败的婴儿中,10例出生时血清HBsAg已为阳性;8～12个月后随访,HBsAg仍持续阳性,提示为宫内感染。本研究证实,孕期注射HBIG未能提高婴儿对HBV疫苗加HBIG的免疫阻断效果。宫内感染可能是疫苗加HBIG免疫阻断失败的主要原因。采用降低孕妇血清HBVDNA的措施,如对孕妇进行抗HBV治疗,也许能降低HBV宫内感染率。     

A Novel Approach to Phylogeny Reconstruction from Protein Sequences

The reliable reconstruction of tree topology from a set of homologous sequences is one of the main goals in the study of molecular evolution. If consistent estimators of distances from a multiple sequence alignment are known, the distance method is attractive because the tree reconstruction is consistent. To obtain a distance estimate d, the observed proportion of differences p (p-distance) is usually ``corrected' for multiple and back substitutions by means of a functional relationship d=f(p). In this paper the conditions under which this correction of p-distances will not alter the selection of the tree topology are specified. When these conditions are not fulfilled the selection of the tree topology may depend on the correction function applied. A novel method which includes estimates of distances not only between sequence pairs, but between triplets, quadruplets, etc., is proposed to strengthen the proper selection of correction function and tree topology. A ``super' tree that includes all tree topologies as special cases is introduced. Received: 17 February 1998 / Accepted: 20 July 1998     

Evidence of Gene Conversion Events Between Paralogous Sequences Produced by Tetraploidization in Salmoninae Fish

We investigated the occurrence of gene conversions between paralogous sequences of Salmoninae derived from ancestral tetraploidization and their effect on the evolutionary history of DNA sequences. A microsatellite with long flanking regions (750 bp) including both coding and noncoding sequences was analyzed. Microsatellite size polymorphism was used to detect the alleles of both paralogous counterparts and infer linkage arrangement between loci. DNA sequencing of seven Salmoninae species revealed that paralogous sequences were highly differentiated within species, especially for noncoding regions. Ten gene conversion events between paralogous sequences were inferred. While these events appears to have homogenized regions of otherwise highly differential paralogous sequences, they amplified the differentiation among orthologous sequences. Their effects were larger on coding than on noncoding regions. As a consequence, noncoding sequences grouped by orthologous lineages in phylogenetic trees, whereas coding regions grouped by taxa. Based upon these results, we present a model showing how gene conversion events may also result in the PCR amplification of nonorthologous sequences in different taxa, with obvious complications for phylogenetic inferences, comparative mapping, and population genetic studies. Received: 11 October 2000 / Accepted: 18 September 2001     

Ribosomal External and Internal Transcribed Spacers: Combined Use in the Phylogenetic Analysis of Medicago (Leguminosae)

We have estimated the potential phylogenetic utility of the ribosomal external transcribed spacer (ETS) from the nuclear ribosomal region. The ETS was sequenced from 13 annual Medicago (Fabaceae) species upstream a highly conserved motive which was found among many different organisms. In the genus Medicago, the ETS was found to evolve 1.5 times faster than the internal transcribed spacer and to be 1.5 times more informative. Reduced ribosomal maturation process constraints on ETS are proposed to explain the different evolutionary rates between the two spacers. Maximal phylogenetic resolution and support was obtained when the two spacers were analyzed together. No incongruence between the two spacers was found and ETS appears to be a valuable source of information for solidifying ITS plant phylogeny. The phylogeny obtained in Medicago suggests that none of the three subsections included in the study is monophyletic. Received: 15 April 1997 / Accepted: 29 July 1997     

Characterization of Mitochondrial Small-Subunit Ribosomal RNAs from Holoparasitic Plants

Mitochondrial small-subunit (19S) rDNA sequences were obtained from 10 angiosperms to further characterize sequence divergence levels and structural variation in this molecule. These sequences were derived from seven holoparasitic (nonphotosynthetic) angiosperms as well as three photosynthetic plants. 19S rRNA is composed of a conservative core region (ca. 1450 nucleotides) as well as two variable regions (V1 and V7). In pairwise comparisons of photosynthetic angiosperms to Glycine, the core 19S rDNA sequences differed by less than 1.4%, thus supporting the observation that variation in mitochondrial rDNA is 3–4 times lower than seen in protein coding and rDNA genes of other subcellular organelles. Sequences representing four distinct lineages of nonasterid holoparasites showed significantly increased numbers of substitutions in their core 19S rDNA sequences (2.3–7.6%), thus paralleling previous findings that showed accelerated rates in nuclear (18S) and plastid (16S) rDNA from the same plants. Relative rate tests confirmed the accelerated nucleotide substitution rates in the holoparasites whereas rates in nonparasitic plants were not significantly increased. Among comparisons of both parasitic and nonparasitic plants, transversions outnumbered transitions, in many cases more than two to one. The core 19S rRNA is conserved in sequence and structure among all nonparasitic angiosperms whereas 19S rRNA from members of holoparasitic Balanophoraceae have unique extensions to the V5 and V6 variable domains. Substitution and insertion/deletion mutations characterized the V1 and V7 regions of the nonasterid holoparasites. The V7 sequence of one holoparasite (Scybalium) contained repeat motifs. The cause of substitution rate increases in the holoparasites does not appear to be a result of RNA editing, hence the underlying molecular mechanism remains to be fully documented. Received: 18 May 1997 / Accepted: 11 July 1997     

Phylogenetic analysis of the first complete hepatitis E virus (HEV) genome from Africa

Hepatitis E virus (HEV) is globally distributed, transmitted enterically and between humans and animals. Phylogenetic analysis has identified five distinct HEV genotypes. The first full-length sequence of an African strain (Chad) is presented and compared to 31 complete HEV genomes available, including the fulminant hepatitis strain from India, swine strains and a strain from Morocco. The two African strains are more closely related to genotype 1 than to any other genotypes and together they possibly form a sub-genotype or sixth genotype. The first evidence for recombination between divergent HEV strains is presented.     

Antigenic Drift of Viruses Within a Host: A Finite Site Model with Demographic Stochasticity

We theoretically study the antigenic drift of viruses within an infected host, as observed in human immunodeficiency virus (HIV) and equine infectious anemia virus (EIAV) infections, assuming that a finite number of antigen-determining sites at the viral envelop gene are responsible for the specific immune response. The pattern of antigen evolution becomes more complex than that predicted from the previous one-dimensional antigen space models. If the viral growth rate is sufficiently large, the demographic stochasticity for the fate of a new antigen mutant can be neglected. The high dimensionality in the way a virus escapes the immune defense in genotype space could then causes a rapid increase in the antigenic diversity and the total viral density, until finally the whole antigen genotypes are used up. The viral population is then driven to extinction in a host by the enhanced immune response to all genotypes. In contrast, if the viral growth rate is moderate or small so that only a small fraction of new antigen mutants can survive during the initial endangered period of random extinction, the viral antigenic diversity and the total density remain bounded, thereby enabling them to persist for a prolonged period by shifting the dominant antigen types. The phylogenetic pattern of antigen divergence is well characterized by the mean number of surviving antigen mutants from an antigen genotype. The substitution rate at antigen-determining sites increases as the efficiency of host immune response increases. Received: 11 January 2000 / Accepted: 25 May 2000     

Vaccine- and hepatitis B immune globulin-induced escape mutations of hepatitis B virus surface antigen

Hepatitis B virus surface antigen (HBsAg) vaccination has been shown to be effective in preventing hepatitis B virus (HBV) infection. The protection is based on the induction of anti-HBs antibodies against a major cluster of antigenic epitopes of HBsAg, defined as the 'a' determinant region of small HBsAg. Prophylaxis of recurrent HBV infection in patients who have undergone liver transplantation for hepatitis B-related end-stage liver disease is achieved by the administration of hepatitis B immune globulins (HBIg) derived from HBsAg-vaccinated subjects. The anti-HBs-mediated immune pressure on HBV, however, seems to go along with the emergence and/or selection of immune escape HBV mutants that enable viral persistence in spite of adequate antibody titers. These HBsAg escape mutants harbor single or double point mutations that may significantly alter the immunological characteristics of HBsAg. Most escape mutations that influence HBsAg recognition by anti-HBs antibodies are located in the second 'a' determinant loop. Notably, HBsAg with an arginine replacement for glycine at amino acid 145 is considered the quintessential immune escape mutant because it has been isolated consistently in clinical samples of HBIg-treated individuals and vaccinated infants of chronically infected mothers. Direct binding studies with monoclonal antibodies demonstrated a more dramatic impact of this mutation on anti-HBs antibody recognition, compared with other point mutations in this antigenic domain. The clinical and epidemiological significance of these emerging HBsAg mutants will be a matter of research for years to come, especially as data available so far document that these mutants are viable and infectious strains. Strategies for vaccination programs and posttransplantation prophylaxis of recurrent hepatitis need to be developed that may prevent immune escape mutant HBV from spreading and to prevent these strains from becoming dominant during the next decennia.     

A Reexamination of the Phylogenetic Position of Callimico (Primates) Incorporating New Mitochondrial DNA Sequence Data

The New World monkeys are divided into two main groups, Callitrichidae and Cebidae. Callimico goeldii shares traits with both the Cebidae and the Callitrichidae. Recent morphological phyletic studies generally place Callimico as the most basal member of the Callitrichidae. In contrast, genetic studies (immunological, restriction fragment, and sequence data) have consistently placed Callimico somewhere within the Callitrichidae, not basal to this clade. A DNA sequence data set from the terminal 236 codons of the mitochondrial ND4 gene and the tRNA^His, tRNA^Ser, and tRNA^Leu genes was generated to clarify the position of Callimico. The sequences of 887 base pairs were analyzed by maximum-parsimony, neighbor-joining, and maximum-likelihood methods. The results of these various methods are generally congruent and place Callimico within the Callitrichidae between the marmosets (Callithrix and Cebuella) and the tamarins (Saguinus and Leontopithecus). Combined analyses of all suitable nuclear and mitochondrial gene sequences confirm the position of Callimico between the marmosets and the tamarins. As available molecular evidence indicates that Callimico is more closely related to the marmosets than to the tamarins, a reconsideration of the morphological evidence in light of the consensus tree from DNA sequence analyses is warranted. The marmosets and tamarins share four morphological characters (loss of the third molar, loss of the hypocone, reduced body size, reproductive twinning). Dwarfism may have evolved repeatedly among the Callitrichidae. It is well-known that the loss of a character can occur many times independently. The reproduction of marmosets and tamarins is extremely specialized and it is difficult to imagine that this complex and unique twinning system evolved separately in marmosets and tamarins. However, it is possible that a secondary reversal to single offspring took place in Callimico. Received: 20 March 1997 / Accepted: 17 December 1997     

Phylogenetic Position of the Subgenus Lordiphosa of the Genus Drosophila (Diptera: Drosophilidae) Inferred from Alcohol Dehydrogenase (Adh) Gene Sequences

We analyzed the phylogenetic relationship between the species of Lordiphosa and other Drosophilidae using alcohol dehydrogenase (Adh) gene sequences. The phylogenetic trees consistently show that the four species Drosophila kurokawai, D. collinella, D. stackelbergi, and D. clarofinis, which include three species groups of Lordiphosa, form a monophyletic clade. This clade is placed as a sister group to the willistoni and saltans groups of Sophophora. On the other hand, three species of Lordiphosa, D. tenuicauda, D. pseudotenuicauda, and D. acutissima, all of which belong to the tenuicauda group, are not shown to be related to the major Lordiphosa lineage. In the phylogenetic trees, these species are included into the clade comprised of Drosophila and Hirtodrosophila, although it remains uncertain whether the tenuicauda group is a monophyletic group or not. These results indicate that Lordiphosa is polyphyletic and that most of the members of the subgenus have a close relationship to the neotropical groups of Sophophora. The above conclusion is compatible with the hypothesis of Okada (Mushi [1963] 37:79–100) and Lastovka and Máca (Acta Ent Bohemoslov [1978] 75:404–420) that Lordiphosa is most closely related to Sophophora; in contrast, our results contradict the hypothesis of Grimaldi (Bull Am Mus Nat Hist [1990] 197:1–139) that Lordiphosa is a sister group to the genus Scaptomyza. Received: 12 May 1999 / Accepted: 14 April 2000     

Phylogenetic Analysis of Vertebrate Lactate Dehydrogenase (LDH) Multigene Families

In this paper we analyzed 49 lactate dehydrogenase (LDH) sequences, mostly from vertebrates. The amino acid sequence differences were found to be larger for a human–killifish pair than a human–lamprey pair. This indicates that some protein sequence convergence may occur and reduce the sequence differences in distantly related species. We also examined transitions and transversions separately for several species pairs and found that the transitions tend to be saturated in the distantly related species pair, while transversions are increasing. We conclude that transversions maintain a conservative rate through the evolutionary time. Kimura's two-parameter model for multiple-hit correction on transversions only was used to derive a distance measure and then construct a neighbor-joining (NJ) tree. Three findings were revealed from the NJ tree: (i) the branching order of the tree is consistent with the common branch pattern of major vertebrates; (ii) Ldh-A and Ldh-B genes were duplicated near the origin of vertebrates; and (iii) Ldh-C and Ldh-A in mammals were produced by an independent gene duplication in early mammalian history. Furthermore, a relative rate test showed that mammalian Ldh-C evolved more rapidly than mammalian Ldh-A. Under a two-rate model, this duplication event was calibrated to be approximately 247 million years ago (mya), dating back to the Triassic period. Other gene duplication events were also discovered in Xenopus, the first duplication occurring approximately 60–70 mya in both Ldh-A and Ldh-B, followed by another recent gene duplication event, approximately 20 mya, in Ldh-B. Received: 5 October 2001 / Accepted: 24 October 2001     

The Effect of Recombination on the Accuracy of Phylogeny Estimation

Phylogenetic studies based on DNA sequences typically ignore the potential occurrence of recombination, which may produce different alignment regions with different evolutionary histories. Traditional phylogenetic methods assume that a single history underlies the data. If recombination is present, can we expect the inferred phylogeny to represent any of the underlying evolutionary histories? We examined this question by applying traditional phylogenetic reconstruction methods to simulated recombinant sequence alignments. The effect of recombination on phylogeny estimation depended on the relatedness of the sequences involved in the recombinational event and on the extent of the different regions with different phylogenetic histories. Given the topologies examined here, when the recombinational event was ancient, or when recombination occurred between closely related taxa, one of the two phylogenies underlying the data was generally inferred. In this scenario, the evolutionary history corresponding to the majority of the positions in the alignment was generally recovered. Very different results were obtained when recombination occurred recently among divergent taxa. In this case, when the recombinational breakpoint divided the alignment in two regions of similar length, a phylogeny that was different from any of the true phylogenies underlying the data was inferred.     

The Complete Mitochondrial DNA Sequence of the Pig (Sus scrofa)

The complete mitochondrial genome sequence of the pig, Sus scrofa, was determined. The length of the sequence presented is 16,679 nucleotides. This figure is not absolute, however, due to pronounced heteroplasmy caused by variable numbers of the motif GTACACGTGC in the control region of different molecules. A phylogenetic study was performed on the concatenated amino acid and nucleotide sequences of 12 protein-coding genes of the mitochondrial genome. The analysis identified the pig (Suiformes) as a sister group of a cow/whale clade, making Artiodactyla paraphyletic. The split between pig and cow/whale was molecularly dated at 65 million years before present. Received: 2 December 1997 / Accepted: 20 February 1998     

乙型肝炎病毒X蛋白对细胞凋亡的影响存在剂量依赖性

乙型肝炎病毒(HBV)X蛋白(HBx)与HBV相关肝细胞肝癌(HCC)的发生和发展密切相关.深入研究HBx在HCC形成中的作用将为探索HBV致癌机制提供重要依据.HBx是多功能蛋白,其对细胞凋亡的影响至今仍存在分歧.许多研究表明,HBx既有促进细胞凋亡又有抑制细胞凋亡的功能,但原因不清楚.本研究中,将表达HBx的质粒短...     

Extreme Endemic Radiation of the Malagasy Vangas (Aves: Passeriformes)

Phylogenetic relationships of the family Vangidae and representatives of several other passeriform families were inferred from 882 base positions of mitochondrial DNA sequences of 12S and 16S rRNA genes. Results indicated the monophyly of the Vangidae, which includes the genus Tylas, hitherto often placed in the family Pycnonotidae. Our results also revealed the Malagasy endemic Newtonia, a genus never previously assigned to the Vangidae, to be a member of this family. These results suggest the occurrence of an extensive in situ radiation of this family within Madagascar, and that the extant high diversity of this family is not the result of multiple colonizations from outside. The extremely high morphological and ecological diversification of the family seems to have been enhanced through the use and ultimate occupancy of vacant niches in this island. Received: 8 September 2000 / Accepted: 13 February 2001