Investigation of the evolutionary history of the lyssaviruses

Dear Editor, We report the results of evolutionary history estimation of the lyssaviruses based on an analysis of the Glycoprotein (G) sequences gene using the BEAST software package.The most recent common ancestor (TMRCA) of all the lyssavirus strains was estimated to be approximately 5030 years (95％ HPD 3988-6069 years),and there was a significant spread of the rabies virus throughout the world range in the last 200 years,consistent with significant time points in development and migration of human civilizations.We speculate that increased and expansion of human migration during this time period may have promote the increase in lyssavirus diversity.In addition,evidence of host switching in lyssavirus history from the Chiroptera to the Carnivora orders was also identified.     

Complete genome sequence of ikoma lyssavirus

Lyssaviruses (family Rhabdoviridae) constitute one of the most important groups of viral zoonoses globally. All lyssaviruses cause the disease rabies, an acute progressive encephalitis for which, once symptoms occur, there is no effective cure. Currently available vaccines are highly protective against the predominantly circulating lyssavirus species. Using next-generation sequencing technologies, we have obtained the whole-genome sequence for a novel lyssavirus, Ikoma lyssavirus (IKOV), isolated from an African civet in Tanzania displaying clinical signs of rabies. Genetically, this virus is the most divergent within the genus Lyssavirus. Characterization of the genome will help to improve our understanding of lyssavirus diversity and enable investigation into vaccine-induced immunity and protection.     

Evidence of two Lyssavirus phylogroups with distinct pathogenicity and immunogenicity

The genetic diversity of representative members of the Lyssavirus genus (rabies and rabies-related viruses) was evaluated using the gene encoding the transmembrane glycoprotein involved in the virus-host interaction, immunogenicity, and pathogenicity. Phylogenetic analysis distinguished seven genotypes, which could be divided into two major phylogroups having the highest bootstrap values. Phylogroup I comprises the worldwide genotype 1 (classic Rabies virus), the European bat lyssavirus (EBL) genotypes 5 (EBL1) and 6 (EBL2), the African genotype 4 (Duvenhage virus), and the Australian bat lyssavirus genotype 7. Phylogroup II comprises the divergent African genotypes 2 (Lagos bat virus) and 3 (Mokola virus). We studied immunogenic and pathogenic properties to investigate the biological significance of this phylogenetic grouping. Viruses from phylogroup I (Rabies virus and EBL1) were found to be pathogenic for mice when injected by the intracerebral or the intramuscular route, whereas viruses from phylogroup II (Mokola and Lagos bat viruses) were only pathogenic by the intracerebral route. We showed that the glycoprotein R333 residue essential for virulence was naturally replaced by a D333 in the phylogroup II viruses, likely resulting in their attenuated pathogenicity. Moreover, cross-neutralization distinguished the same phylogroups. Within each phylogroup, the amino acid sequence of the glycoprotein ectodomain was at least 74% identical, and antiglycoprotein virus-neutralizing antibodies displayed cross-neutralization. Between phylogroups, the identity was less than 64.5% and the cross-neutralization was absent, explaining why the classical rabies vaccines (phylogroup I) cannot protect against lyssaviruses from phylogroup II. Our tree-axial analysis divided lyssaviruses into two phylogroups that more closely reflect their biological characteristics than previous serotypes and genotypes.     

Enhanced Passive Bat Rabies Surveillance in Indigenous Bat Species from Germany - A Retrospective Study

In Germany, rabies in bats is a notifiable zoonotic disease, which is caused by European bat lyssaviruses type 1 and 2 (EBLV-1 and 2), and the recently discovered new lyssavirus species Bokeloh bat lyssavirus (BBLV). As the understanding of bat rabies in insectivorous bat species is limited, in addition to routine bat rabies diagnosis, an enhanced passive surveillance study, i.e. the retrospective investigation of dead bats that had not been tested for rabies, was initiated in 1998 to study the distribution, abundance and epidemiology of lyssavirus infections in bats from Germany. A total number of 5478 individuals representing 21 bat species within two families were included in this study. The Noctule bat (Nyctalus noctula) and the Common pipistrelle (Pipistrellus pipistrellus) represented the most specimens submitted. Of all investigated bats, 1.17% tested positive for lyssaviruses using the fluorescent antibody test (FAT). The vast majority of positive cases was identified as EBLV-1, predominately associated with the Serotine bat (Eptesicus serotinus). However, rabies cases in other species, i.e. Nathusius'' pipistrelle bat (Pipistrellus nathusii), P. pipistrellus and Brown long-eared bat (Plecotus auritus) were also characterized as EBLV-1. In contrast, EBLV-2 was isolated from three Daubenton''s bats (Myotis daubentonii). These three cases contribute significantly to the understanding of EBLV-2 infections in Germany as only one case had been reported prior to this study. This enhanced passive surveillance indicated that besides known reservoir species, further bat species are affected by lyssavirus infections. Given the increasing diversity of lyssaviruses and bats as reservoir host species worldwide, lyssavirus positive specimens, i.e. both bat and virus need to be confirmed by molecular techniques.     

A simple sandwich ELISA (WELYSSA) for the detection of lyssavirus nucleocapsid in rabies suspected specimens using mouse monoclonal antibodies

Monoclonal antibody (MAb)-based capture enzyme-linked immunosorbent assays (ELISA) were developed for the diagnosis of rabies-suspect specimens. A combination of four mouse monoclonal antibodies directed against the rabies virus nucleocapsid was selected and used for the detection. The test was optimized and standardized so that maximum concordance could be maintained with the standard procedures of rabies diagnosis recommended by the WHO expert committee. Using prototype viruses from the different genotypes of lyssavirus and from various geographic origins and phylogenetic lineages, this paper presents a reliable, rapid and transferable diagnostic method, named WELYSSA that readily permits the detection of lyssaviruses belonging to the 7 genotypes of lyssavirus circulating in Europe, Africa, Asia and Oceania. The threshold of detection of lyssavirus nucleocapsids is low (0.8 ng/ml). With a panel of 1030 specimens received for rabies diagnostic testing, this test was found to be highly specific (0.999) and sensitive (0.970) when compared to other recommended rabies diagnostic methods.     

Isolation of Irkut Virus from a Murina leucogaster Bat in China

Background and objectives

Bats are recognized as a major reservoir of lyssaviruses; however, no bat lyssavirus has been isolated in Asia except for Aravan and Khujand virus in Central Asia. All Chinese lyssavirus isolates in previous reports have been of species rabies virus, mainly from dogs. Following at least two recent bat-associated human rabies-like cases in northeast China, we have initiated a study of the prevalence of lyssaviruses in bats in Jilin province and their public health implications. A bat lyssavirus has been isolated and its pathogenicity in mice and genomic alignment have been determined.

Results

We report the first isolation of a bat lyssavirus in China, from the brain of a northeastern bat, Murina leucogaster. Its nucleoprotein gene shared 92.4%/98.9% (nucleotide) and 92.2%/98.8% (amino acid) identity with the two known Irkut virus isolates from Russia, and was designated IRKV-THChina12. Following intracranial and intramuscular injection, IRKV-THChina12 produced rabies-like symptoms in adult mice with a short inoculation period and high mortality. Nucleotide sequence analysis showed that IRKV-THChina12 has the same genomic organization as other lyssaviruses and its isolation provides an independent origin for the species IRKV.

Conclusions

We have identified the existence of a bat lyssavirus in a common Chinese bat species. Its high pathogenicity in adult mice suggests that public warnings and medical education regarding bat bites in China should be increased, and that surveillance be extended to provide a better understanding of Irkut virus ecology and its significance for public health.     

Gene duplication and evidence for balancing selection acting on MHC class II DAA gene of the half-smooth tongue sole (Cynoglossus semilaevis)

Allelic polymorphism and evolution mechanism of major histocompatibility complex (MHC) genes has been investigated in many mammals, however, much less is known in teleost. In order to investigate the mechanisms creating and maintaining variability at the MHC class II DAA locus, we examined the polymorphism, gene duplication and balancing selection of MHC class II DAA gene of the half-smooth tongue sole (Cynoglossus semilaevis). We described 33 alleles in the C. semilaevis, recombination and gene duplication seems to play more important roles in the origin of new alleles. The rate of non-synonymous substitutions (d(N)) occurred at a significantly higher frequency than that of synonymous substitutions (d(S)) in peptide-binding region (PBR) and non-PBR, suggesting balancing selection for maintaining polymorphisms at the MHC II DAA locus. Many positive selection sites were found to act very intensively on antigen-binding sites. Our founding suggests a snapshot in an evolutionary process of MHC-DAA gene evolution of the C. semilaevis.     

Identification of a Role for Nucleolin in Rabies Virus Infection

Rabies virus replicates in the cytoplasm of host cells, but rabies virus phosphoprotein (P-protein) undergoes active nucleocytoplasmic trafficking. Here we show that the largely nuclear P-protein isoform P3 can localize to nucleoli and forms specific interactions with nucleolin. Importantly, depletion of nucleolin expression inhibits viral protein expression and infectious virus production by infected cells. This provides the first evidence that lyssaviruses interact with nucleolin and that nucleolin is important to lyssavirus infection.     

Experimental Lagos bat virus infection in straw-colored fruit bats: A suitable model for bat rabies in a natural reservoir species

Rabies is a fatal neurologic disease caused by lyssavirus infection. Bats are important natural reservoir hosts of various lyssaviruses that can be transmitted to people. The epidemiology and pathogenesis of rabies in bats are poorly understood, making it difficult to prevent zoonotic transmission. To further our understanding of lyssavirus pathogenesis in a natural bat host, an experimental model using straw-colored fruit bats (Eidolon helvum) and Lagos bat virus, an endemic lyssavirus in this species, was developed. To determine the lowest viral dose resulting in 100% productive infection, bats in five groups (four bats per group) were inoculated intramuscularly with one of five doses, ranging from 10^0.1 to 10^4.1 median tissue culture infectious dose (TCID₅₀). More bats died due to the development of rabies after the middle dose (10^2.1 TCID₅₀, 4/4 bats) than after lower (10^1.1, 2/4; 10^1.1, 2/4) or higher (10^3.1, 2/4; 10^4.1, 2/4) doses of virus. In the two highest dose groups, 4/8 bats developed rabies. Of those bats that remained healthy 3/4 bats seroconverted, suggesting that high antigen loads can trigger a strong immune response that abrogates a productive infection. In contrast, in the two lowest dose groups, 3/8 bats developed rabies, 1/8 remained healthy and seroconverted and 4/8 bats remained healthy and did not seroconvert, suggesting these doses are too low to reliably induce infection. The main lesion in all clinically affected bats was meningoencephalitis associated with lyssavirus-positive neurons. Lyssavirus antigen was detected in tongue epithelium (5/11 infected bats) rather than in salivary gland epithelium (0/11), suggesting viral excretion via the tongue. Thus, intramuscular inoculation of 10^2.1 TCID₅₀ of Lagos bat virus into straw-colored fruit bats is a suitable model for lyssavirus associated bat rabies in a natural reservoir host, and can help with the investigation of lyssavirus infection dynamics in bats.     

Raccoons (Procyon lotor) in Germany as potential reservoir species for Lyssaviruses

Raccoons can be found almost everywhere in Germany since their first successful introduction in 1934. Although the animal is a well-known reservoir species for rabies in the USA, during the last European fox rabies epizootic, only a few rabid raccoons were reported from Germany. In recent years, the raccoon population density has increased tremendously, especially in (semi) urban settings. Presently, Germany is free of terrestrial wildlife rabies. To assess the potential risk that the raccoon population in Germany could act as a reservoir species upon reemergence of rabies, the susceptibility of the local raccoon population was investigated. Wild-caught animals were inoculated with the most likely lyssavirus variants to infect the local population. It was shown that the raccoons were fully susceptible for a dog and raccoon rabies virus isolate. Also, five of six raccoons inoculated with a fox rabies virus isolate showed clinical signs. However, none of the raccoons infected with European Bat Lyssavirus type 1 succumbed to rabies; meanwhile, all these raccoons seroconverted. It is concluded that the highest risk for the raccoon population in Germany to become infected with lyssaviruses is through the importation of rabies infected dogs.     

Genomic diversity and evolution of the lyssaviruses

Lyssaviruses are RNA viruses with single-strand, negative-sense genomes responsible for rabies-like diseases in mammals. To date, genomic and evolutionary studies have most often utilized partial genome sequences, particularly of the nucleoprotein and glycoprotein genes, with little consideration of genome-scale evolution. Herein, we report the first genomic and evolutionary analysis using complete genome sequences of all recognised lyssavirus genotypes, including 14 new complete genomes of field isolates from 6 genotypes and one genotype that is completely sequenced for the first time. In doing so we significantly increase the extent of genome sequence data available for these important viruses. Our analysis of these genome sequence data reveals that all lyssaviruses have the same genomic organization. A phylogenetic analysis reveals strong geographical structuring, with the greatest genetic diversity in Africa, and an independent origin for the two known genotypes that infect European bats. We also suggest that multiple genotypes may exist within the diversity of viruses currently classified as 'Lagos Bat'. In sum, we show that rigorous phylogenetic techniques based on full length genome sequence provide the best discriminatory power for genotype classification within the lyssaviruses.     

Quantifying Antigenic Relationships among the Lyssaviruses

All lyssaviruses cause fatal encephalitis in mammals. There is sufficient antigenic variation within the genus to cause variable vaccine efficacy, but this variation is difficult to characterize quantitatively: sequence analysis cannot yet provide detailed antigenic information, and antigenic neutralization data have been refractory to high-resolution robust interpretation. Here, we address these issues by using state-of-the-art antigenic analyses to generate a high-resolution antigenic map of a global panel of 25 lyssaviruses. We compared the calculated antigenic distances with viral glycoprotein ectodomain sequence data. Although 67% of antigenic variation was predictable from the glycoprotein amino acid sequence, there are in some cases substantial differences between genetic and antigenic distances, thus highlighting the risk of inferring antigenic relationships solely from sequence data at this time. These differences included epidemiologically important antigenic differences between vaccine strains and wild-type rabies viruses. Further, we quantitatively assessed the antigenic relationships measured by using rabbit, mouse, and human sera, validating the use of nonhuman experimental animals as a model for determining antigenic variation in humans. The use of passive immune globulin is a crucial component of rabies postexposure prophylaxis, and here we also show that it is possible to predict the reactivity of immune globulin against divergent lyssaviruses.Rabies remains a globally important zoonosis, despite being one of the oldest recognized infectious diseases (27, 55). The majority of rabies in terrestrial animals and humans is caused by classical rabies virus (RABV), a lyssavirus in the family Rhabdoviridae. Since the 1950s, many related lyssaviruses which are capable of causing clinical rabies have been identified. The majority of those viruses have been isolated from bats (Chiroptera), including four divergent viruses, which were isolated in separate geographic locations throughout Eurasia in the past 18 years (2, 29, 31). The Chiroptera, therefore, represent a global reservoir for lyssaviruses, creating the potential for “spillover” infection to terrestrial mammals, including humans. Occasionally transmission between members of a new host species will occur, with potential for a subsequent adaptation in that species (35). Phylogenetic evidence suggests that one or more host-switching events from bats into terrestrial mammals were originally responsible for the ongoing global epidemic of terrestrial RABV (6).Pre- or postexposure prophylaxis, using vaccination and passive immune globulin administration according to World Health Organization (WHO) guidelines, is currently the only effective way to prevent rabies after infection with a lyssavirus (1). The efficacy of both active and passive immunization is likely to be affected by antigenic differences between viruses. The lyssavirus trimeric glycoprotein is the primary surface antigen, the major target for neutralizing antibodies (8), and is involved in cell binding and entry (34, 36, 53). Antigenic sites on the glycoprotein have been described using monoclonal antibody escape mutants (8, 16, 47, 51). These studies have elucidated two major sites (sites II and III) and multiple minor sites. Although estimates of antigenic differences can be made using information regarding these known antigenic sites, protein structure, and amino acid properties, predictions of the relative importance of those sites and specific mutations within those sites cannot be quantitatively tested without a method to reliably measure antigenic effect.The use of serological cross neutralization data to measure antigenic difference is limited by the reliability of the serological test and, more importantly, by paradoxes, or irregularities in the data. These irregularities include higher heterologous than homologous titers and individual variations between sera raised against the same antigen (22, 52). Hence, serological data are considered to have low resolution, and they are often used only qualitatively. Despite these difficulties, studies have attempted to further quantify antigenic differences among lyssaviruses. Badrane et al. (5) showed a linear correlation between the glycoprotein amino acid identity and four cross neutralization titers. Other studies have demonstrated variable serological cross-reactivity between European bat lyssaviruses (EBLV) and RABVs (10, 11) and suggested that antigenic relationships between EBLV-1 and EBLV-2 may not be fully reflected in the genetic relationships (41). Recent investigations into the efficacy of biologics against the Eurasian lyssaviruses showed an array of relatedness between lyssavirus species, with, for example, a murine anti-Aravan virus (anti-ARAV) serum neutralizing Khujand virus (KHUV) and ARAV equally but an anti-KHUV serum being less effective at neutralizing ARAV than KHUV (22). Until recently, however, there were no established tools for the quantitative analysis of antigenic data.Here we resolve the issue of quantitative interpretation of antigenic data using antigenic cartography. Antigenic cartography is a theory and associated computational method that resolves the paradoxes in the interpretation of antigenic data and makes possible high-resolution quantitative analyses and visualizations of binding assay data (15, 20, 25, 44, 49, 52).Integrating antigenic data with direct sequencing data, here we quantify the antigenic and genetic variation among a global panel of lyssaviruses, including representatives from all lyssavirus species. Furthermore, we address two key issues in the development of antilyssavirus biologics: the appropriateness of animal models and the development of efficacious alternatives to human rabies immune globulins (HRIGs).     

Bat rabies surveillance in France: first report of unusual mortality among serotine bats

Background

Rabies is a fatal viral encephalitic disease that is caused by lyssaviruses which can affect all mammals, including human and bats. In Europe, bat rabies cases are attributed to five different lyssavirus species, the majority of rabid bats being attributed to European bat 1 lyssavirus (EBLV-1), circulating mainly in serotine bats (Eptesicus serotinus). In France, rabies in bats is under surveillance since 1989, with 77 positive cases reported between 1989 and 2016.

Case presentation

In the frame of the bat rabies surveillance, an unusual mortality of serotine bats was reported in 2009 in a village in North-East France. Six juvenile bats from an E. serotinus maternity colony counting ~200 individuals were found to be infected with EBLV-1. The active surveillance of the colony by capture sessions of bats from July to September 2009 showed a high detection rate of neutralising EBLV-1 antibodies (≈ 50%) in the colony. Moreover, one out of 111 animals tested was found to shed viable virus in saliva, while lyssavirus RNA was detected by RT-PCR for five individuals.

Conclusion

This study demonstrated that the lyssavirus infection in the serotine maternity colony was followed by a high rate of bat rabies immunity after circulation of the virus in the colony. The ratio of seropositive bats is probably indicative of an efficient virus transmission coupled to a rapid circulation of EBLV-1 in the colony.

     

Development of a Mouse Monoclonal Antibody Cocktail for Post-exposure Rabies Prophylaxis in Humans

As the demand for rabies post-exposure prophylaxis (PEP) treatments has increased exponentially in recent years, the limited supply of human and equine rabies immunoglobulin (HRIG and ERIG) has failed to provide the required passive immune component in PEP in countries where canine rabies is endemic. Replacement of HRIG and ERIG with a potentially cheaper and efficacious alternative biological for treatment of rabies in humans, therefore, remains a high priority. In this study, we set out to assess a mouse monoclonal antibody (MoMAb) cocktail with the ultimate goal to develop a product at the lowest possible cost that can be used in developing countries as a replacement for RIG in PEP. Five MoMAbs, E559.9.14, 1112-1, 62-71-3, M727-5-1, and M777-16-3, were selected from available panels based on stringent criteria, such as biological activity, neutralizing potency, binding specificity, spectrum of neutralization of lyssaviruses, and history of each hybridoma. Four of these MoMAbs recognize epitopes in antigenic site II and one recognizes an epitope in antigenic site III on the rabies virus (RABV) glycoprotein, as determined by nucleotide sequence analysis of the glycoprotein gene of unique MoMAb neutralization-escape mutants. The MoMAbs were produced under Good Laboratory Practice (GLP) conditions. Unique combinations (cocktails) were prepared, using different concentrations of the MoMAbs that were capable of targeting non-overlapping epitopes of antigenic sites II and III. Blind in vitro efficacy studies showed the MoMab cocktails neutralized a broad spectrum of lyssaviruses except for lyssaviruses belonging to phylogroups II and III. In vivo, MoMAb cocktails resulted in protection as a component of PEP that was comparable to HRIG. In conclusion, all three novel combinations of MoMAbs were shown to have equal efficacy to HRIG and therefore could be considered a potentially less expensive alternative biological agent for use in PEP and prevention of rabies in humans.     

Detection and genetic characterization of rabies virus from human patients

Saliva and blood were collected from two patients who had not received post exposure prophylaxis in the cities of Wenzhou and Xinning respectively. Both patients were confirmed as positive for rabies by detection of rabies virus specific nucleoprotein antibodies in the sera by Western Blot. However, rabies virus specific RNA was only identified in the saliva collected from the patient in Wenzhou. Furthermore, the isolate Zhejiang Wz0 (H) was obtained by inoculating one-day-old suckling mice. Both nucleoprotein (N) and glycoprotein (G) genes from the isolate were amplified by RT-PCR and sequenced. Phylogenetic analysis indicated that the isolate belonged to classic rabies virus, and shared a higher homology with the street viruses from dogs in the main endemic areas in China and the street virus from dogs in Indonesia than with other known strains. Further comparison of the deduced amino acid sequences between the isolate and the vaccine strains used in China showed that the virus had a higher level of homology with the vaccine strain CTN than with the other vaccine strains (3aG, PV, PM and ERA). In particular, amino acid residues substitutions located in antigenic site III in the G protein, which could react with the neutralizing antibodies, were observed. These results suggested that the virus belonged to the classic rabies virus, and both N and G genes diverged from the current vaccine strains used in China at either the nucleotide or the amino acid level. Foundation items: The project was supported by the Grants (2003BA712A08-02 and 2004BA718B03) from the Chinese Ministry of Science and Technology.     

Codon-based tests of positive selection,branch lengths,and the evolution of mammalian immune system genes

Using basic probability theory, we show that there is a substantial likelihood that even in the presence of strong purifying selection, there will be a number of codons in which the number of synonymous nucleotide substitutions per site (d (S)) exceeds the number of non-synonymous nucleotide substitutions per site (d (N)). In an empirical study, we examined the numbers of synonymous (b (S)) and non-synonymous substitutions (b (N)) along branches of the phylogenies of 69 single-copy orthologous genes from seven species of mammals. A pattern of b (N) > b (S) was most commonly seen in the shortest branches of the tree and was associated with a high coefficient of variation in both b (N) and b (S), suggesting that high stochastic error in b (N) and b (S) on short branches, rather than positive Darwinian selection, is the explanation of most cases where b (N) is greater than b (S) on a given branch. The branch-site method of Zhang et al. (Zhang, Nielsen, Yang, Mol Biol Evol, 22:2472-2479, 2005) identified 117 codons on 35 branches as "positively selected," but a majority of these codons lacked synonymous substitutions, while in the others, synonymous and non-synonymous differences per site occurred in approximately equal frequencies. Thus, it was impossible to rule out the hypothesis that chance variation in the pattern of mutation across sites, rather than positive selection, accounted for the observed pattern. Our results showed that b (N)/b (S) was consistently elevated in immune system genes, but neither the search for branches with b (N) > b (S) nor the branch-site method revealed this trend.     

Human immunodeficiency virus seroconversion and evolution of the hepatitis C virus quasispecies

When chronic hepatitis C virus (HCV) infections are complicated by acquisition of human immunodeficiency virus (HIV), liver disease appears to accelerate and serum levels of HCV RNA may rise. We hypothesized that HIV might affect the HCV quasispecies by decreasing both complexity (if HIV-induced immunosuppression lessens pressure for selecting HCV substitutions) and the ratio of nonsynonymous (d(N)) to synonymous (d(S)) substitutions, because d(N) may be lower (if there is less selective pressure). To test this hypothesis, we studied the evolution of HCV sequences in 10 persons with chronic HCV infection who seroconverted to HIV and, over the next 3 years, had slow or rapid progression of HIV-associated disease. From each subject, four serum specimens were selected with reference to HIV seroconversion: (i) more than 2 years prior, (ii) less than 2 years prior, (iii) less than 2 years after, and (iv) more than 2 years after. The HCV quasispecies in these specimens was characterized by generating clones containing 1 kb of cDNA that spanned the E1 gene and the E2 hypervariable region 1 (HVR1), followed by analysis of clonal frequencies (via electrophoretic migration) and nucleotide sequences. We examined 1,320 cDNA clones (33 per time point) and 287 sequences (median of 7 per time point). We observed a trend toward lower d(N)/d(S) after HIV seroconversion in 7 of 10 subjects and lower d(N)/d(S) in those with rapid HIV disease progression. However, the magnitude of these differences was small. These results are consistent with the hypothesis that HIV infection alters the HCV quasispecies, but the number of subjects and observation time may be too low to characterize the full effect.     

Allelic variation, balancing selection and positive selected sites detected from MHC class Iα gene of olive flounder

Major histocompatibility complex (MHC) genes play an important role in the immune response of vertebrates. Allelic polymorphism and pattern of evolution in MHC genes has been investigated in many mammals, however, much less is known in teleost. In the present study, we have investigated complete MHC Iα gene consists of 7 exons and 6 introns in Olive flounder (Paralichthys olivaceus). Genetic variation in the MHC class Iα gene was also tested in flounder. In 32 individuals, a total of 62 alleles were detected from exon 2 of MHC class Iα gene. The rate of non-synonymous substitutions (d ( N )) occurred at a significantly higher frequency than that of synonymous substitutions (d ( S )) in PBR and non-PBR, suggesting that balancing selection for maintaining polymorphisms at the MHC Iα locus. Many positive selection sites were found to act very intensively on antigen binding sites. Our founding suggests a snapshot in an evolutionary process of MHC Iα gene evolution of the P. olivaceus.