Fifty-million-year-old polymorphism at an immunoglobulin variable region gene locus in the rabbit evolutionary lineage

Molecular time estimates for the Lagomorpha diversification

Despite their importance as members of the Glires group, lagomorph diversification processes have seldom been studied using molecular data. Notably, only a few phylogenetic studies have included most of the examined lagomorph lineages. Previous studies that included a larger sample of taxa and markers used nonconservative tests to support the branches of their proposed phylogeny. The objective of this study was to test the monophyly of families and genera of lagomorphs and to evaluate the group diversification process. To that end, this work expanded the sampling of markers and taxa in addition to implementing the bootstrap, a more rigorous statistical test to measure branch support; hence, a more robust phylogeny was recovered. Our supermatrix included five mitochondrial genes and 14 nuclear genes for eighty-eight taxa, including three rodent outgroups. Our maximum likelihood tree showed that all tested genera and both families, Leporidae and Ochotonidae, were recovered as monophyletic. In the Ochotona genus, the subgenera Conothoa and Pika, but not Ochotona, were recovered as monophyletic. Six calibration points based on fossils were used to construct a time tree. A calibration test was performed (via jackknife) by removing one calibration at a time and estimating divergence times for each set. The diversification of the main groups of lagomorphs indicated that the origin of the order's crown group was dated from the beginning of the Palaeogene. Our diversification time estimates for Lagomorpha were compared with those for the largest mammalian order, i.e., rodent lineages in Muroidea. According to our time-resolved phylogenetic tree, the leporids underwent major radiation by evolving a completely new morphospace-larger bodies and an efficient locomotor system-that enabled them to cover wide foraging areas and outrun predators more easily than rodents and pikas.

Trans-specific polymorphism and the convergent evolution of supertypes in major histocompatibility complex class II genes in darters (Etheostoma)

Major Histocompatibility Complex (MHC) genes are one of the most polymorphic gene groups known in vertebrates. MHC genes also exhibit allelic variants that are shared among taxa, referred to as trans-specific polymorphism (TSP). The role that selection plays in maintaining such high diversity within species, as well as TSP, is an ongoing discussion in biology. In this study, we used deep-sequencing techniques to characterize MHC class IIb gene diversity in three sympatric species of darters. We found at least 5 copies of the MHC gene in darters, with 126 genetic variants encoding 122 unique amino acid sequences. We identified four supertypes based on the binding properties of proteins encoded by the sequences. Although each species had a unique pool of variants, many variants were shared between species pairs and across all three species. Phylogenetic analysis showed that the variants did not group together monophyletically based on species identity or on supertype. An expanded phylogenetic analysis showed that some darter alleles grouped together with alleles from other percid fishes. Our findings show that TSP occurs in darters, which suggests that balancing selection is acting at the genotype level. Supertypes, however, are most likely evolving convergently, as evidenced by the fact that alleles do not form monophyletic groups based on supertype. Our research demonstrates that selection may be acting differently on MHC genes at the genotype and supertype levels, selecting for the maintenance of high genotypic diversity while driving the convergent evolution of similar MHC phenotypes across different species.

An overview of the lagomorph immune system and its genetic diversity

Our knowledge of the lagomorph immune system remains largely based upon studies of the European rabbit (Oryctolagus cuniculus), a major model for studies of immunology. Two important and devastating viral diseases, rabbit hemorrhagic disease and myxomatosis, are affecting European rabbit populations. In this context, we discuss the genetic diversity of the European rabbit immune system and extend to available information about other lagomorphs. Regarding innate immunity, we review the most recent advances in identifying interleukins, chemokines and chemokine receptors, Toll-like receptors, antiviral proteins (RIG-I and Trim5), and the genes encoding fucosyltransferases that are utilized by rabbit hemorrhagic disease virus as a portal for invading host respiratory and gut epithelial cells. Evolutionary studies showed that several genes of innate immunity are evolving by strong natural selection. Studies of the leporid CCR5 gene revealed a very dramatic change unique in mammals at the second extracellular loop of CCR5 resulting from a gene conversion event with the paralogous CCR2. For the adaptive immune system, we review genetic diversity at the loci encoding antibody variable and constant regions, the major histocompatibility complex (RLA) and T cells. Studies of IGHV and IGKC genes expressed in leporids are two of the few examples of trans-species polymorphism observed outside of the major histocompatibility complex. In addition, we review some endogenous viruses of lagomorph genomes, the importance of the European rabbit as a model for human disease studies, and the anticipated role of next-generation sequencing in extending knowledge of lagomorph immune systems and their evolution.

Identification of Allorecognition Loci in Neurospora crassa by Genomics and Evolutionary Approaches

Understanding the genetic and molecular bases of the ability to distinguish self from nonself (allorecognition) and mechanisms underlying evolution of allorecognition systems is an important endeavor for understanding cases where it becomes dysfunctional, such as in autoimmune disorders. In filamentous fungi, allorecognition can result in vegetative or heterokaryon incompatibility, which is a type of programmed cell death that occurs following fusion of genetically different cells. Allorecognition is genetically controlled by het loci, with coexpression of any combination of incompatible alleles triggering vegetative incompatibility. Herein, we identified, characterized, and inferred the evolutionary history of candidate het loci in the filamentous fungus Neurospora crassa. As characterized het loci encode proteins carrying an HET domain, we annotated HET domain genes in 25 isolates from a natural population along with the N. crassa reference genome using resequencing data. Because allorecognition systems can be affected by frequency-dependent selection favoring rare alleles (i.e., balancing selection), we mined resequencing data for HET domain loci whose alleles displayed elevated levels of variability, excess of intermediate frequency alleles, and deep gene genealogies. From these analyses, 34 HET domain loci were identified as likely to be under balancing selection. Using transformation, incompatibility assays and genetic analyses, we determined that one of these candidates functioned as a het locus (het-e). The het-e locus has three divergent allelic groups that showed signatures of positive selection, intra- and intergroup recombination, and trans-species polymorphism. Our findings represent a compelling case of balancing selection functioning on multiple alleles across multiple loci potentially involved in allorecognition.

Trans-Species Polymorphism in Immune Genes: General Pattern or MHC-Restricted Phenomenon?

Immunity exhibits extraordinarily high levels of variation. Evolution of the immune system in response to host-pathogen interactions in particular ecological contexts appears to be frequently associated with diversifying selection increasing the genetic variability. Many studies have documented that immunologically relevant polymorphism observed today may be tens of millions years old and may predate the emergence of present species. This pattern can be explained by the concept of trans-species polymorphism (TSP) predicting the maintenance and sharing of favourable functionally important alleles of immune-related genes between species due to ongoing balancing selection. Despite the generality of this concept explaining the long-lasting adaptive variation inherited from ancestors, current research in TSP has vastly focused only on major histocompatibility complex (MHC). In this review we summarise the evidence available on TSP in human and animal immune genes to reveal that TSP is not a MHC-specific evolutionary pattern. Further research should clearly pay more attention to the investigation of TSP in innate immune genes and especially pattern recognition receptors which are promising candidates for this type of evolution. More effort should also be made to distinguish TSP from convergent evolution and adaptive introgression. Identification of balanced TSP variants may represent an accurate approach in evolutionary medicine to recognise disease-resistance alleles.

Evolution of dimorphisms of the proteasome subunit beta type 8 gene (PSMB8) in basal ray-finned fish

The proteasome subunit beta type 8 (PSMB8) gene encodes a catalytic subunit of immunoproteasome that plays a central role in the processing of antigenic peptides presented by major histocompatibility complex class I molecules. The A- and F-type alleles defined by the 31st amino acid residue determining cleaving specificity have been identified from ray-finned fish, amphibia, and reptiles. These two types show extremely long-term trans-species polymorphism in Polypteriformes, Cypriniformes, and Salmoniformes, suggesting the presence of very ancient lineages termed A and F. To elucidate the evolution of the PSMB8 dimorphism in basal ray-finned fish, we analyzed Pantodon buchholzi (Osteoglossiformes), seven species of Anguilliformes, and Hypomesus nipponensis (Osmeriformes). Both A and F lineage sequences were identified from P. buchholzi and H. nipponensis, confirming that these two lineages have been conserved by basal ray-finned fish. However, both the A- and F-type alleles found in Anguilliformes species belonged to the F lineage irrespective of their types. This apparently suggests that the A lineage was lost in the common ancestor of Anguilliformes, and recovery of the A type within the F lineage occurred in Anguilliformes. The apparent loss of the F lineage and recovery of the F type within the A lineage have already been reported from tetrapods and higher teleosts. However, this is the first report on the reverse situation and reveals the dynamic evolution of the PSMB8 dimorphism.

A novel rabbit monoclonal antibody platform to dissect the diverse repertoire of antibody epitopes for HIV-1 Env immunogen design

The majority of available monoclonal antibodies (MAbs) in the current HIV vaccine field are generated from HIV-1-infected people. In contrast, preclinical immunogenicity studies have mainly focused on polyclonal antibody responses in experimental animals. Although rabbits have been widely used for antibody studies, there has been no report of using rabbit MAbs to dissect the specificity of antibody responses for AIDS vaccine development. Here we report on the production of a panel of 12 MAbs from a New Zealand White (NZW) rabbit that was immunized with an HIV-1 JR-FL gp120 DNA prime and protein boost vaccination regimen. These rabbit MAbs recognized a diverse repertoire of envelope (Env) epitopes ranging from the highly immunogenic V3 region to several previously underappreciated epitopes in the C1, C4, and C5 regions. Nine MAbs showed cross-reactivity to gp120s of clades other than clade B. Increased somatic mutation and extended CDR3 were observed with Ig genes of several molecularly cloned rabbit MAbs. Phylogenic tree analysis showed that the heavy chains of MAbs recognizing the same region on gp120 tend to segregate into an independent subtree. At least three rabbit MAbs showed neutralizing activities with various degrees of breadth and potency. The establishment of this rabbit MAb platform will significantly enhance our ability to test optimal designs of Env immunogens to gain a better understanding of the structural specificity and evolution process of Env-specific antibody responses elicited by candidate AIDS vaccines.

Sequencing of modern Lepus VDJ genes shows that the usage of VHn genes has been retained in both Oryctolagus and Lepus that diverged 12 million years ago

Among mammals, the European rabbit (Oryctolagus cuniculus) has a unique mechanism of generating the primary antibody repertoire. Despite having over 200 VH genes, the VH1 gene, the most D-proximal VH gene, is used in 80-90 % of VDJ rearrangements, while the remaining 10-20 % is encoded by the VHn genes that map at least 100 Kb upstream of VH1. The maintenance of the VHn genes usage in low frequency in VDJ rearrangements has been suggested to represent a relic of an ancestral immunologic response to pathogens. To address this question, we sequenced VDJ genes for another leporid, genus Lepus, which separated from European rabbit 12 million years ago. Approximately 25 VDJ gene sequences were obtained for each one of three Lepus europaeus individuals. We found that Lepus also uses the VHn genes in 5-10 % of its VDJ rearrangements. Our results show that the VHn genes are a conserved ancestral polymorphism that has been maintained in the leporids genome and is being used for the generation of VDJ rearrangements by both modern Lepus and Oryctolagus.

Where are the wormy mice? A reexamination of hybrid parasitism in the European house mouse hybrid zone

Wormy mice in a hybrid zone have been interpreted as evidence of low hybrid fitness, such that parasites contribute to species separation. However, because of its natural heterogeneity, observations of parasite load must be numerous with good field area coverage. We sampled 689 mice from 107 localities across the Bavaria-Bohemia region of the European house mouse hybrid zone and calculated their hybrid indices using 1401 diagnostic single nucleotide polymorphisms (SNPs). We tested whether hybrids have greater or lesser diversity and load of parasite helminths than additive expectations, performing load analyses on the four most common taxa. We found hybrids have significantly reduced diversity and load of each of the commonest helminths; rarer helminths further support reduced load. Although within-locality comparisons have little power, randomization tests show the repeated pattern is unlikely to be due to local parasite heterogeneity, and simulations show a patch of low parasite diversity is unlikely to fall by chance just so in the field area, such that it produces the observed effects. Our data therefore contradict the idea that helminths reduce hybrid fitness through increased load. We discuss a vicariant Red Queen model that implies immune genes tracking parasites will escape Dobzhansky-Muller incompatibilities, generating hybrid variants untargeted by parasites.

Long-lived dichotomous lineages of the proteasome subunit beta type 8 (PSMB8) gene surviving more than 500 million years as alleles or paralogs

On an evolutionary time scale, polymorphic alleles are believed to have a short life, persisting at most tens of millions of years even under long-term balancing selection. Here, we report highly diverged trans-species dimorphism of the proteasome subunit beta type 8 (PSMB8) gene, which encodes a catalytic subunit of the immunoproteasome responsible for the generation of peptides presented by major histocompatibility complex (MHC) class I molecules, in lower teleosts including Cypriniformes (zebrafish and loach) and Salmoniformes (trout and salmon), whose last common ancestor dates to 300 Ma. Moreover, phylogenetic analyses indicated that these dimorphic alleles share lineages with two shark paralogous genes, suggesting that these two lineages have been maintained for more than 500 My either as alleles or as paralogs, and that conversion between alleles and paralogs has occurred at least once during vertebrate evolution. Two lineages termed PSMB8A and PSMB8F show an A(31)F substitution that would probably affect their cleaving specificity, and whereas the PSMB8A lineage has been retained by all analyzed jawed vertebrates, the PSMB8F lineage has been lost by most jawed vertebrates except for cartilaginous fish and basal teleosts. However, a possible functional equivalent of the PSMB8F lineage has been revived as alleles within the PSMB8A lineage at least twice during vertebrate evolution in the amphibian Xenopus and teleostean Oryzias species. Dynamic evolution of the PSMB8 polymorphism through long-term persistence, loss, and regaining of dimorphism and conversion between alleles and paralogs implies the presence of strong selective pressure for functional polymorphism of this gene.

Molecular bases of genetic diversity and evolution of the immunoglobulin heavy chain variable region (IGHV) gene locus in leporids

The rabbit has long been a model for studies of the immune system. Work using rabbits contributed both to the battle against infectious diseases such as rabies and syphilis, and to our knowledge, of antibodies' structure, function, and regulated expression. With the description of rabbit Ig allotypes, the discovery of different gene segments encoding immunoglobulins became possible. This challenged the "one gene-one protein" dogma. The observation that rabbit allotypic specificities of the variable regions were present on IgM and IgG molecules also led to the hypothesis of Ig class switching. Rabbit allotypes contributed to the documentation of phenomena such as allelic exclusion and imbalance in production of allelic gene products. During the last 30 years, the rabbit Ig allotypes revealed a number of unique features, setting them apart from mice, humans, and other mammals. Here, we review the most relevant findings concerning the rabbit IGHV. Among these are the preferential usage of one VH gene in VDJ rearrangements, the existence of trans-species polymorphism in the IGHV locus revealed by serology and confirmed by sequencing IGHV genes in Lepus, the unusually large genetic distances between allelic lineages and the fact that the antibody repertoire is diversified in this species only after birth. The whole genome sequence of a rabbit, plus re-sequencing of additional strains and related genera, will allow further evolutionary investigations of antibody variation. Continued research will help define the roles that genetic, allelic, and population diversity at antibody loci may play in host-parasite interactions.

Factors important in evolutionary shaping of immunoglobulin gene loci

Background

The extraordinary diversity characterizing the antibody repertoire is generated by both evolution and lymphocyte development. Much of this diversity is due to the existence of immunoglobulin (Ig) variable region gene segment libraries, which were diversified during evolution and, in higher vertebrates, are used in generating the combinatorial diversity of antibody genes. The aim of the present study was to address the following questions: What evolutionary parameters affect the size and structure of gene libraries? Are the number of genes in libraries of contemporary species, and the corresponding gene locus structure, a random result of evolutionary history, or have these properties been optimized with respect to individual or population fitness? If a larger number of genes or different genome structures do not increase the fitness, then the current structure is probably optimized.

Results

We used a simulation of variable region gene library evolution. We measured the effect of different parameters on gene library size and diversity, and the corresponding fitness. We found compensating relationships between parameters, which optimized Ig library size and diversity.

Conclusions

We conclude that contemporary species' Ig libraries have been optimized by evolution in terms of Ig sequence lengths, the number and diversity of Ig genes, and antibody-antigen affinities.

Transspecies dimorphic allelic lineages of the proteasome subunit beta-type 8 gene (PSMB8) in the teleost genus Oryzias

The proteasome subunit β-type 8 (PSMB8) gene in the jawed vertebrate MHC genomic region encodes a catalytic subunit of the immunoproteasome involved in the generation of peptides to be presented by the MHC class I molecules. A teleost, the medaka (Oryzias latipes), has highly diverged dimorphic allelic lineages of the PSMB8 gene with only about 80% amino acid identity, termed "PSMB8d" and "PSMB8N," which have been retained by most wild populations analyzed. To elucidate the evolutionary origin of these two allelic lineages, seven species of the genus Oryzias were analyzed for their PSMB8 allelic sequences using a large number of individuals from wild populations. All the PSMB8 alleles of these species were classified into one of these two allelic lineages based on their nucleotide sequences of exons and introns, indicating that the Oryzias PSMB8 gene has a truly dichotomous allelic lineage. Retention of both allelic lineages was confirmed except for one species. The PSMB8d lineage showed a higher frequency than the PSMB8N lineage in all seven species. The two allelic lineages showed curious substitutions at the 31st and 53rd residues of the mature peptide, probably involved in formation of the S1 pocket, suggesting that these allelic lineages show a functional difference in cleavage specificity. These results indicate that the PSMB8 dimorphism was established before speciation within the genus Oryzias and has been maintained for more than 30-60 million years under a strict and asymmetric balancing selection through several speciation events.

Evolutionary history of an MHC gene in two leporid species: characterisation of Mhc-DQA in the European brown hare and comparison with the European rabbit

We surveyed the genetic diversity of the expressed major histocompatibility complex class II DQA locus in natural populations of European brown hares, Lepus europaeus, from Austria and Belgium (267 individuals in total). Based on cDNA sequences, we designed hare-specific primers to amplify the highly variable second exon of the DQA gene. Using cloning-sequencing methodology and capillary electrophoresis single-strand conformation polymorphism, we found ten alleles of the DQA exon 2 locus across these two European regions, of which eight are described for the first time. To search for signals of selection and recombination in the evolution of the DQA gene within the leporids, we augmented our sample with orthologous DQA alleles from the European rabbit, Oryctolagus cuniculus, in order to carry out a species level, species pairwise comparison. We found evidence of recombination in the history of the DQA sequences in leporids with some recombinant alleles bridging the species divide. In both species, selection on peptide binding site codons can be detected, though stronger for the rabbit. This result suggests that there may be a differential selection pressure in the deeper evolutionary history of these two species due to differences in several demographic and ecological traits likely subjecting them to differential selection by parasites. Finally, evolutionary relationships show a widespread and statistically significant intermingling of alleles from the two species. The many macroparasites shared between hares and rabbits may explain this pattern of trans-species polymorphism.

Ancient trans-specific polymorphism at pheromone receptor genes in basidiomycetes

In the majority of sexual organisms, reproduction occurs almost exclusively through the combination of distinct and alternate forms, called sexes or mating types. In some fungi, there can be dozens to hundreds of alternate alleles that determine compatible mating types. Such extensive polymorphism is expected to be maintained by balancing selection, and in extreme cases may give rise to trans-specific polymorphism. Here, we analyzed sequences of two pheromone receptors in the Microbotryum fungal species complex (Basidiomycota), which has only two alternate mating types. Several lines of evidence strongly suggest that the pheromone receptors are two allelic sequences acting to determine the alternate A1 and A2 mating types required for mating in Microbotryum. Phylogenetic trees of pheromone receptors in the Microbotryum species complex indicated a trans-specific polymorphism: the Microbotryum sequences from a given mating type were all more similar to the pheromone receptors of distantly related classes of fungi than to the alternate pheromone receptor in the Microbotryum species. A phylogenetic tree built using other known pheromone receptors from basidiomycetes showed that trans-specific polymorphism is widespread. The pheromone receptor alleles from Microbotryum appeared as the oldest, being at least 370 million years old. This represents the oldest known trans-specific polymorphism known in any organism so far, which may be due to the existence of sex chromosomes, obligate sexuality, mitochondrial inheritance linked to the mating type, and a highly selfing mating system in Microbotryum.

Conservation of inter-protein binding sites in RUSH and RFBP, an ATP11B isoform

Isoforms of RUSH interact with a RING-finger binding protein (RFBP), which is a splice variant of the Type IV P-type ATPase, ATP11B. Splice arrays and RT-PCR showed that although most splice variants in RUSH and ATP11B are conserved in human and rabbit, the RFBP isoform is specific to rabbit. Interactions between the discontinuous PVITHC-HAKCPL sequence in the RING-domain of RUSH and the KVIRLIKIS sequence in the catalytic loop of RFBP were first identified with pull-down assays. Fine mapping involved probing CLIPS-constrained RING peptides with GST-tagged KVIRLIKIS. When the companion site in RFBP was fine mapped by replacement analysis with MBP-tagged RING, a four-fold increase in binding was noted for the KVIRLDKIS mutant. Direct comparison of splicing events in the RUSH and ATP11B genes between human and rabbit shows high structural stability in these protein interactions sites, which are 100% conserved in all mammalian orthologs.

Balancing selection at a frog antimicrobial peptide locus: fluctuating immune effector alleles?

Balancing selection is common on many defense genes, but it has rarely been reported for immune effector proteins such as antimicrobial peptides (AMPs). We describe genetic diversity at a brevinin-1 AMP locus in three species of leopard frogs (Rana pipiens, Rana blairi, and Rana palustris). Several highly divergent allelic lineages are segregating at this locus. That this unusual pattern results from balancing selection is demonstrated by multiple lines of evidence, including a ratio of nonsynonymous/synonymous polymorphism significantly higher than 1, the ZnS test, incongruence between the number of segregating sites and haplotype diversity, and significant Tajima's D values. Our data are more consistent with a model of fluctuating selection in which alleles change frequencies over time than with a model of stable balancing selection such as overdominance. Evidence for fluctuating selection includes skewed allele frequencies, low levels of synonymous variation, nonneutral values of Tajima's D within allelic lineages, an inverse relationship between the frequency of an allelic lineage and its degree of polymorphism, and divergent allele frequencies among populations. AMP loci could be important sites of adaptive genetic diversity, with consequences for host–pathogen coevolution and the ability of species to resist disease epidemics.

Genetic contributions to the autoantibody profile in a rabbit model of Systemic Lupus Erythematosus (SLE)

For the development of rabbit models of Systemic Lupus Erythematosus (SLE), immunoglobulin allotype-defined pedigreed rabbits from the National Institute of Allergy and Infectious Diseases rabbit resource more closely approximate human populations due to their non-inbred pedigreed structure. In an initial study from this laboratory, peptides (SM and GR) from the spliceosomal Smith (Sm) and the NMDA glutamate receptor NR2b, on branched polylysine backbones (BB) elicited antinuclear and anti-dsDNA autoantibodies typical of SLE, as well as seizures and nystagmus sometimes observed as neurological manifestations in SLE patients. This suggested the feasibility of further selective breeding to develop a more reproducible rabbit model for investigations of SLE. Here we report the results of GR-MAP-8 and control BB immunization on autoantibody responses in a group of 24 rabbits specifically bred and developed from parents and ancestors tested for autoantibody responses. The changes in hematological profile and blood chemistry in the experimental rabbits were evaluated along with autoantibody responses. Elevations of total white blood cell (WBC), monocyte, eosinophil and basophil counts that developed following immunizations were moderately influenced by litter and presence of the antibody heavy chain allotype VH1a1. Autoantibody development followed a sequential pattern with anti-nuclear antibodies (ANA) followed by anti-dsDNA and subsequently anti-Sm and anti-RNP similar to SLE patients. High autoantibody levels to one autoantigen were not always associated with antibody response to another. Female rabbits had higher prevalence and levels of autoantibodies similar to human SLE. Higher autoantibody levels of anti-dsDNA and -ANA were observed among some full sibs and the presence of high responder ancestors in the pedigree was associated the augmented responses. We observed significant association between highest antibody responses to GR-MAP-8 and highest anti-dsDNA levels. Naturally occurring autoantibodies were found in some pre-immune sera and some unique ANA fluorescent staining patterns within the experimental group were observed. Background immunofluorescence in pre-immune sera, distinct patterns of programmed autoantibody responses unique among individual rabbits may have been modulated by genetic constitution, gender and environmental factors including exposure to antigens. The high incidence and intensity of autoantibody responses among descendants of high responders suggest that there may be an additive mode of inheritance with high heritability. It is conceivable that further rigorous pedigree selection for autoantibody responses could lead to development of rabbit models with spontaneous occurrence of SLE like serology and disease phenotypes.

Genetic variation at chemokine receptor CCR5 in leporids: alteration at the 2nd extracellular domain by gene conversion with CCR2 in Oryctolagus, but not in Sylvilagus and Lepus species

Whereas in its natural host (Sylvilagus sps.) the effects of myxoma virus infections are benign, in European rabbit (Oryctolagus cuniculus), it causes a highly infectious disease with very high mortality rate, known as myxomatosis. There is evidence that, as with HIV-1 virus in human, myxoma virus may use chemokine receptors such as CCR5 of the host target cell for entry and activation of pathways of immune avoidance. We have characterized and compared CCR5 genes of leporid species with different susceptibility levels to myxomatosis. The CCR5 protein of O. cuniculus differs markedly from all those known from other species. The most striking was the replacement of a specific peptide motif of the second extracellular loop (ECL2) by a motif, which in other species characterizes the CCR2 molecules. While absent in Sylvilagus and Lepus species, this CCR2 imposed CCR5-ECL2 alteration was observed in all genomes of 25 European rabbits, representing the subspecies O. cuniculus algirus and O. cuniculus cuniculus. Allelic variation at the rabbit CCR5 locus confirmed that the gene conversion predates the subspecies split (1-2 Ma).

Concerted and birth-and-death evolution of multigene families

Until around 1990, most multigene families were thought to be subject to concerted evolution, in which all member genes of a family evolve as a unit in concert. However, phylogenetic analysis of MHC and other immune system genes showed a quite different evolutionary pattern, and a new model called birth-and-death evolution was proposed. In this model, new genes are created by gene duplication and some duplicate genes stay in the genome for a long time, whereas others are inactivated or deleted from the genome. Later investigations have shown that most non-rRNA genes including highly conserved histone or ubiquitin genes are subject to this type of evolution. However, the controversy over the two models is still continuing because the distinction between the two models becomes difficult when sequence differences are small. Unlike concerted evolution, the model of birth-and-death evolution can give some insights into the origins of new genetic systems or new phenotypic characters.

The evolution of the immunoglobulin heavy chain variable region (IgV H ) in Leporids: an unusual case of transspecies polymorphism

In domestic rabbit (Oryctolagus cuniculus), three serological types have been distinguished at the variable domain of the antibody H chain, the so-called V(H) a allotypes a1, a2, and a3. They correspond to highly divergent allelic lineages of the V(H) 1 gene, which is the gene rabbit utilizes in more than 80% of VDJ rearrangements. The sharing of serological V(H) a markers between rabbit and snowshoe hare (Lepus americanus) has suggested that the large genetic distances between rabbit V(H) 1 alleles (9-14% nucleotide differences) can be explained by unusually long lineage persistence times (transspecies polymorphism). Because this interpretation of the serological data is uncertain, we have determined the nucleotide sequences of V(H) genes expressed in specimens of Lepus species. Two sequence groups were distinguished, one of which occurred only in hare specimen displaying serological motifs of the rabbit V(H) a-a2 allotype. Sequences of this group are part of a monophyletic cluster containing the V(H) 1 sequences of the rabbit a2 allotype. The fact that this "transspecies a2 cluster" did not include genes of other rabbit V(H) a allotypes (a1, a3, and a4) is incompatible with the existence of a common V(H) a ancestor gene within the species, and suggests that the divergence of the V(H) a lineages preceded the Lepus vs Oryctolagus split. The sequence data are furthermore compatible with the hypothesis that the V(H)a polymorphism can be two times older than the divergence time between the Lepus and Oryctolagus lineages, which was estimated at 16-24 million years.

Dramatic variation of the vomeronasal pheromone receptor gene repertoire among five orders of placental and marsupial mammals

Pheromones are chemicals emitted and sensed by conspecifics to elicit social and sexual responses and are perceived in terrestrial vertebrates primarily by the vomeronasal organ (VNO). Pheromone receptors in the mammalian VNO are encoded by the V1R and V2R gene superfamilies. The V1R superfamily contains 187 and 102 putatively functional genes in the mouse and rat, respectively. To investigate whether this large repertoire size is typical among mammals with functional VNOs, we here describe the V1R repertoires of dog, cow, and opossum based on their draft genome sequences. The dog and cow have only 8 and 32 intact V1R genes, respectively. Thus, the intact V1R repertoire size varies by at least 23-fold among placental mammals with functional VNOs. To our knowledge, this size ratio represents the greatest among-species variation in gene family size of all mammalian gene families. Phylogenetic analysis of placental V1R genes suggests multiple losses of ancestral genes in carnivores and artiodactyls and gains of many new genes by gene duplication in rodents, manifesting massive gene births and deaths. We also identify 49 intact opossum V1R genes and discover independent expansions of the repertoire in placentals and marsupials. We further show a concordance between the V1R repertoire size and the complexity of VNO morphology, suggesting that the latter could indicate the sophistication of pheromone communications within species. In sum, our results demonstrate tremendous diversity and rapid evolution of mammalian V1R gene inventories and caution the generalization of VNO biology from rodents to all mammals.

A molecular supermatrix of the rabbits and hares (Leporidae) allows for the identification of five intercontinental exchanges during the Miocene

The hares and rabbits belonging to the family Leporidae have a nearly worldwide distribution and approximately 72% of the genera have geographically restricted distributions. Despite several attempts using morphological, cytogenetic, and mitochondrial DNA evidence, a robust phylogeny for the Leporidae remains elusive. To provide phylogenetic resolution within this group, a molecular supermatrix was constructed for 27 taxa representing all 11 leporid genera. Five nuclear (SPTBN1, PRKCI, THY, TG, and MGF) and two mitochondrial (cytochrome b and 12S rRNA) gene fragments were analyzed singly and in combination using parsimony, maximum likelihood, and Bayesian inference. The analysis of each gene fragment separately as well as the combined mtDNA data almost invariably failed to provide strong statistical support for intergeneric relationships. In contrast, the combined nuclear DNA topology based on 3601 characters greatly increased phylogenetic resolution among leporid genera, as was evidenced by the number of topologies in the 95% confidence interval and the number of significantly supported nodes. The final molecular supermatrix contained 5483 genetic characters and analysis thereof consistently recovered the same topology across a range of six arbitrarily chosen model specifications. Twelve unique insertion-deletions were scored and all could be mapped to the tree to provide additional support without introducing any homoplasy. Dispersal-vicariance analyses suggest that the most parsimonious solution explaining the current geographic distribution of the group involves an Asian or North American origin for the Leporids followed by at least nine dispersals and five vicariance events. Of these dispersals, at least three intercontinental exchanges occurred between North America and Asia via the Bering Strait and an additional three independent dispersals into Africa could be identified. A relaxed Bayesian molecular clock applied to the seven loci used in this study indicated that most of the intercontinental exchanges occurred between 14 and 9 million years ago and this period is broadly coincidental with the onset of major Antarctic expansions causing land bridges to be exposed.

The evolution of parasite recognition genes in the innate immune system: purifying selection on Drosophila melanogaster peptidoglycan recognition proteins

Genes involved in the recognition of parasites by the acquired immune system are often subject to intense selection pressures. In some cases, selection to recognize a diverse range of parasites has resulted in high levels of polymorphism, while elsewhere the protein sequence has changed rapidly under directional selection. We tested whether parasite recognition genes in the innate immune system show similar patterns of evolution. We sequenced seven peptidoglycan recognition protein genes (PGRPs) from 12 lines of Drosophila melanogaster and one line of D. simulans and used a variety of tests to determine whether the observed mutations were selectively neutral. We were unable to detect either balancing or directional selection. This suggests that the molecular cues used by insects to detect parasites are highly conserved and probably under strong functional constraints which prevent their evolving to evade the host immune response. Therefore, interactions between these genes are unlikely to be the focus of host-parasite coevolution, at least in Drosophila. We also found evidence of gene conversion occurring between two genes, PGRP-SC1A and PGRP-SC1B.

Allelic Variation at the VHa Locus in Natural Populations of Rabbit (Oryctolagus cuniculus, L.)

The large interallelic distances between the three rabbit Ig V(H)a lineages, a1, a2 and a3, suggest that the persistence time of the V(H)a polymorphism could amount to 50 million years, which is much longer than that of MHC polymorphisms. Rabbit originated in the Iberian Peninsula where two subspecies coexist, one of which is confined to Southwestern Iberia (Oryctolagus cuniculus algirus). We studied the V(H) loci in the original species range to obtain a better understanding of the evolutionary history of this unusual polymorphism. Serological surveys revealed that sera from the subspecies algirus, when tested with V(H)a locus-specific alloantisera, showed either cross-reactivity ("a-positive" variants) or no reaction at all ("a-blank"). Using RT-PCR, we determined 120 sequences of rearranged V(H) genes expressed in seven algirus rabbits that were typed as either a-positive or a-blank. The data show that the V(H) genes transcribed in a-positive rabbits are closely related to the V(H)1 alleles of domestic rabbits. In contrast, a-blank rabbits were found to preferentially use V(H) genes that, although clearly related to the known V(H)a genes, define a new major allotypic lineage, designated a4. The a4 sequences have hallmark rabbit V(H)a residues together with a number of unprecedented amino acid changes in framework region 2 and 3. The net protein distances between the V(H)a4 and the V(H)a1, a2, and a3 lineages were 20, 29, and 21% respectively. We conclude that at least four distantly related lineages of the rabbit V(H)a locus exist, one of which seems to be endemic in the Iberian range.

Evolutionary dynamics of the T-cell receptor VB gene family as inferred from the human and mouse genomic sequences

The diversity of T-cell receptors is generated primarily by the variable-region gene families, each of which is composed of a large number of member genes. The entire genomic sequence of the variable region (VB) of the T- cell receptor beta chain from humans and mice has become available. To understand the evolutionary dynamics of the VB gene family, we conducted a phylogenetic analysis of all VB genes from humans and mice, as well as a detailed analysis of internal DNA duplications in the human genomic VB region. The phylogenetic tree obtained shows that human and mouse VB genes intermingle extensively rather than forming two separate clusters and that many gene duplications occurred both before and after the divergence between primates and rodents. Analyzing the genomic maps of transposable elements (e.g., LINEs and SINEs) and relic VB genes in the VB gene region, we present evidence that a 20-kb VB region duplicated tandemly four times in the human lineage during the last 32 Myr, and 6 out of the 15 VB genes in this region have become nonfunctional during this period. Our results show that the VB gene family is subject to evolution by a birth-and-death process rather than to concerted evolution.

Estimation of divergence times from multiprotein sequences for a few mammalian species and several distantly related organisms

When many protein sequences are available for estimating the time of divergence between two species, it is customary to estimate the time for each protein separately and then use the average for all proteins as the final estimate. However, it can be shown that this estimate generally has an upward bias, and that an unbiased estimate is obtained by using distances based on concatenated sequences. We have shown that two concatenation-based distances, i.e., average gamma distance weighted with sequence length (d(2)) and multiprotein gamma distance (d(3)), generally give more satisfactory results than other concatenation-based distances. Using these two distance measures for 104 protein sequences, we estimated the time of divergence between mice and rats to be approximately 33 million years ago. Similarly, the time of divergence between humans and rodents was estimated to be approximately 96 million years ago. We also investigated the dependency of time estimates on statistical methods and various assumptions made by using sequence data from eubacteria, protists, plants, fungi, and animals. Our best estimates of the times of divergence between eubacteria and eukaryotes, between protists and other eukaryotes, and between plants, fungi, and animals were 3, 1.7, and 1.3 billion years ago, respectively. However, estimates of ancient divergence times are subject to a substantial amount of error caused by uncertainty of the molecular clock, horizontal gene transfer, errors in sequence alignments, etc.

An ancient transpecific polymorphism shows extreme divergence in a multitrait cline in an intertidal snail (Nucella lapillus (L.))

Clines in intraspecific genetic variation are frequently associated with an environmental transition. Here, divergence among nucleotide sequences of two nuclear loci, cytosolic and mitochondrial malate dehydrogenase (cMDH and mMDH, respectively), is described, in a multitrait cline over a distance of ca. 3 km where shell phenotype, allozyme, mitochondrial DNA haplotype, and centric fusion (Robertsonian translocations) frequencies covary with temperature and humidity and change abruptly in a continuous population of the dog-whelk (Nucella lapillus), a common intertidal snail of the north temperate Atlantic. Protein electrophoresis has already shown two alleles of mMDH varying from fixation of one allele to near fixation of the other, whereas cMDH appears to be monomorphic. The results of this study show a striking disparity in nucleotide sequence divergence among alleles at the two loci, with extreme molecular differentiation in one of them. Four alleles of cMDH were found to have nucleotide and amino acid sequence divergences of 0.4% and 0.3%, respectively. In contrast, the two mMDH cDNA alleles differed by 23% and 20% at the nucleotide and amino acid levels, respectively. Analysis of a 91-bp partial nucleotide sequence of mMDH from Nucella freycineti, the closest relative of N. lapillus, revealed two similar alleles and indicated that the divergence in mMDH in N. lapillus represents an ancient transpecific polymorphism in these Nucella. Together with earlier studies on variation in N. lapillus, it is argued that the polymorphism in mMDH and the clines in N. lapillus represent the presence of two persistent coadapted gene complexes, multitrait coevolving genetic solutions to environmental variation, which may presently enable this snail to exploit a diverse environment successfully.