Gene conversions within the skeletal myosin multigene family.

Comparisons of the nucleotide sequences of the light meromyosin (LMM) region of developmentally regulated fast chicken myosin heavy chain (MHC) isoforms indicates that chicken MHC isoforms are more similar to each other than to MHC isoforms in other species. The sequence data provide evidence that gene conversion events have occurred recently among the isoforms. An embryonic (Cemb1) isoform and neonatal isoform have the most extensive regions of sequence identity. Similar gene conversion events are present in the rat alpha- and beta-cardiac MHCs, but were not obvious in the LMM of developmentally regulated fast human MHC isoforms. The data suggest that gene conversion events can play a significant role in the evolution of the MHC multigene families and that concerted evolution of the chicken multigene family occurred after the divergence of mammals and avians.     

Gene duplication and gene conversion in class II MHC genes of New Zealand robins (Petroicidae)

In contrast to mammals, the evolution of MHC genes in birds appears to be characterized by high rates of gene duplication and concerted evolution. To further our understanding of the evolution of passerine MHC genes, we have isolated class II B sequences from two species of New Zealand robins, the South Island robin (Petroica australis australis), and the endangered Chatham Island black robin (Petroica traversi). Using an RT-PCR based approach we isolated four transcribed class II B MHC sequences from the black robin, and eight sequences from the South Island robin. RFLP analysis indicated that all class II B loci were contained within a single linkage group. Analysis of 3-untranslated region sequences enabled putative orthologous loci to be identified in the two species, and indicated that multiple rounds of gene duplication have occurred within the MHC of New Zealand robins. The orthologous relationships are not retained within the coding region of the gene, instead the sequences group within species. A number of putative gene conversion events were identified across the length of our sequences that may account for this. Exon 2 sequences are highly diverse and appear to have diverged under balancing selection. It is also possible that gene conversion involving short stretches of sequence within exon 2 adds to this diversity. Our study is the first report of putative orthologous MHC loci in passerines, and provides further evidence for the importance of gene duplication and gene conversion in the evolution of the passerine MHC.Nucleotide sequence data reported in this paper are available in the GenBank database under the accession numbers AY258333–AY258335, AY428561–AY428570, and AY530534–AY530535     

Gene Conversion in the Evolution of Both the H-2 and Qa Class I Genes of the Murine Major Histocompatibility Complex

In order to better understand the role of gene conversion in the evolution of the class I gene family of the major histocompatibility complex (MHC), we have used a computer algorithm to detect clustered sequence similarities among 24 class I DNA sequences from the H-2, Qa, and Tla regions of the murine MHC. Thirty-four statistically significant clusters were detected; individual analysis of the clusters suggested at least 25 past gene conversion or recombination events. These clusters are comparable in size to the conversions observed in the spontaneously occurring H-2K(bm) and H-2K(km2) mutations, and are distributed throughout all exons of the class I gene. Thus, gene conversion does not appear to be restricted to the regions of the class I gene encoding their antigen-presentation function. Moreover, both the highly polymorphic H-2 loci and the relatively monomorphic Qa and Tla loci appear to have participated as donors and recipients in conversion events. If gene conversion is not limited to the highly polymorphic loci of the MHC, then another factor, presumably natural selection, must be responsible for maintaining the observed differences in level of variation.     

Evolution of MHC class IIB in the genome of wild and ornamental guppies, Poecilia reticulata

This is the first study to quantify genomic sequence variation of the major histocompatibility complex (MHC) in wild and ornamental guppies, Poecilia reticulata. We sequenced 196-219 bp of exon 2 MHC class IIB (DAB) in 56 wild Trinidadian guppies and 14 ornamental strain guppies. Each of two natural populations possessed high allelic richness (15-16 alleles), whereas only three or fewer DAB alleles were amplified from ornamental guppies. The disparity in allelic richness between wild and ornamental fish cannot be fully explained by fixation of alleles by inbreeding, nor by the presence of non-amplified sequences (ie null alleles). Rather, we suggest that the same allele is fixed at duplicated MHC DAB loci owing to gene conversion. Alternatively, the number of loci in the ornamental strains has contracted during >100 generations in captivity, a hypothesis consistent with the accordion model of MHC evolution. We furthermore analysed the substitution patterns by making pairwise comparisons of sequence variation at the putative peptide binding region (PBR). The rate of non-synonymous substitutions (dN) only marginally exceeded synonymous substitutions (dS) in PBR codons. Highly diverged sequences showed no evidence for diversifying selection, possibly because synonymous substitutions have accumulated since their divergence. Also, the substitution pattern of similar alleles did not show evidence for diversifying selection, plausibly because advantageous non-synonymous substitutions have not yet accumulated. Intermediately diverged sequences showed the highest relative rate of non-synonymous substitutions, with dN/dS>14 in some pairwise comparisons. Consequently, a curvilinear relationship was observed between the dN/dS ratio and the level of sequence divergence.     

EVALUATING PATTERNS OF CONVERGENT EVOLUTION AND TRANS‐SPECIES POLYMORPHISM AT MHC IMMUNOGENES IN TWO SYMPATRIC STICKLEBACK SPECIES

The immunologically important major histocompatibility complex (MHC) harbors some of the most polymorphic genes in vertebrates. These genes presumably evolve under parasite‐mediated selection and frequently show inconsistent allelic genealogies, where some alleles are more similar between species than within species. This phenomenon is thought to arise either from convergent evolution under parallel selection or from the preservation of ancient allelic lineages beyond speciation events (trans‐species polymorphism, TSP). Here, we examine natural populations of two sympatric stickleback species (Gasterosteus aculeatus and Pungitius pungitius) to investigate the contribution of these two mechanisms to the evolution of inconsistent allelic genealogies at the MHC. Overlapping parasite taxa between the two host species in three different habitats suggest contemporary parallel selection on the MHC genes. Accordingly, we detected a lack of species‐specific phylogenetic clustering in the immunologically relevant antigen‐binding residues of the MHC IIB genes which contrasted with the rest of the coding and noncoding sequence. However, clustering was not habitat‐specific and a codon‐usage analysis revealed patterns of similarity by descent. In this light, common descent via TSP, in combination with intraspecies gene conversion, rather than convergent evolution is the more strongly supported scenario for the inconsistent genealogy at the MHC.     

Independent evolution of functional MHC class II DRB genes in New World bat species

Genes of the major histocompatibility complex (MHC) play a pivotal role in the vertebrate immune system and are attractive markers for functional, fitness-related, genetic variation. Although bats (Chiroptera) represent the second largest mammalian order and are prone to various emerging infectious diseases, little is known about MHC evolution in bats. In the present study, we examined expressed MHC class II DRB sequences (exons 1 to 4) of New World bat species, Saccopteryx bilineata, Carollia perspicillata, Noctilio albiventris and Noctilio leporinus (only exon 2). We found a wide range of copy number variation of DRB loci with one locus detected in the genus Noctilio and up to ten functional loci observed in S. bilineata. Sequence variation between alleles of the same taxa was high with evidence for positive selection. We found statistical support for recombination or gene conversion events among sequences within the same but not between bat species. Phylogenetic relationships among DRB alleles provided strong evidence for independent evolution of the functional MHC class II DRB genes in the three investigated species, either by recent gene duplication, or homogenization of duplicated loci by frequent gene conversion events. Phylogenetic analysis of all available chiropteran DRB exon 2 sequences confirmed their monophyletic origin within families, but revealed a possible trans-species mode of evolution pattern in congeneric bat species, e.g. within the genera Noctilio and Myotis. This is the first study investigating phylogenetic relationships of MHC genes within bats and therefore contributes to a better understanding of MHC evolution in one of the most dominant mammalian order.     

The evolution of MHC diversity: evidence of intralocus gene conversion and recombination in a single-locus system

Gene conversion, the unidirectional exchange of genetic material between homologous sequences, is thought to strongly influence patterns of genetic diversity. The high diversity of major histocompatibility complex (MHC) genes in many species is thought to reflect a long history of gene conversion events both within and among loci. Theoretical work suggests that intra- and interlocus gene conversion leave characteristic signatures of nucleotide diversity, but empirical studies of MHC variation have rarely been able to analyze the effects of conversion events in isolation, due to the presence of multiple gene copies in most species. The potbellied seahorse (Hippocampus abdominalis), a species with a single copy of the MH class II beta-chain gene (MHIIb), provides an ideal system in which to explore predictions on the effects of intralocus gene conversion on patterns of genetic diversity. The genetic diversity of the MHIIb peptide binding region (PBR) is high in the seahorse, similar to other vertebrate species. In contrast, the remainder of the gene shows a total absence of synonymous variation and low levels of intronic sequence diversity, concentrated in 3 short repetitive regions and 1-12 SNPs per intron. The distribution of substitutions across the gene results in a patchwork pattern of shared polymorphism between otherwise divergent sequences. The pattern of nucleotide diversity observed in the seahorse MHIIb gene is congruent with theoretical expectations for intralocus gene conversion, indicating that this evolutionary mechanism has played an important role in MHC gene evolution, contributing to both the high diversity in the PBR and the low diversity outside this region. Neutral variation at this locus may be further reduced due to biases in nucleotide composition and functional constraints.     

Gene conversion rapidly generates major histocompatibility complex diversity in recently founded bird populations

Population bottlenecks can restrict variation at functional genes, reducing the ability of populations to adapt to new and changing environments. Understanding how populations generate adaptive genetic variation following bottlenecks is therefore central to evolutionary biology. Genes of the major histocompatibility complex (MHC) are ideal models for studying adaptive genetic variation due to their central role in pathogen recognition. While de novo MHC sequence variation is generated by point mutation, gene conversion can generate new haplotypes by transferring sections of DNA within and across duplicated MHC loci. However, the extent to which gene conversion generates new MHC haplotypes in wild populations is poorly understood. We developed a 454 sequencing protocol to screen MHC class I exon 3 variation across all 13 island populations of Berthelot's pipit (Anthus berthelotii). We reveal that just 11-15 MHC haplotypes were retained when the Berthelot's pipit dispersed across its island range in the North Atlantic ca. 75,000 years ago. Since then, at least 26 new haplotypes have been generated in situ across populations. We show that most of these haplotypes were generated by gene conversion across divergent lineages, and that the rate of gene conversion exceeded that of point mutation by an order of magnitude. Gene conversion resulted in significantly more changes at nucleotide sites directly involved with pathogen recognition, indicating selection for functional variants. We suggest that the creation of new variants by gene conversion is the predominant mechanism generating MHC variation in genetically depauperate populations, thus allowing them to respond to pathogenic challenges.     

Complete nucleotide and encoded amino acid sequence of a mammalian myosin heavy chain gene. Evidence against intron-dependent evolution of the rod

The complete nucleotide sequence and exon/intron structure of the rat embryonic skeletal muscle myosin heavy chain (MHC) gene has been determined. This gene comprises 24 X 10(3) bases of DNA and is split into 41 exons. The exons encode a 6035 nucleotide (nt) long mRNA consisting of 90 nt of 5' untranslated, 5820 nt of protein coding and 125 nt of 3' untranslated sequence. The rat embryonic MHC polypeptide is encoded by exons 3 to 41 and contains 1939 amino acid residues with a calculated Mr of 223,900. Its amino acid sequence displays the structural features typical for all sarcomeric MHCs, i.e. an amino-terminal "globular" head region and a carboxy-terminal alpha-helical rod portion that shows the characteristics of a coiled coil with a superimposed 28-residue repeat pattern interrupted at only four positions by "skip" residues. The complex structure of the rat embryonic MHC gene and the conservation of intron locations in this and other MHC genes are indicative of a highly split ancestral sarcomeric MHC gene. Introns in the rat embryonic gene interrupt the coding sequence at the boundaries separating the proteolytic subfragments of the head, but not at the head/rod junction or between the 28-residue repeats present within the rod. Therefore, there is little evidence for exon shuffling and intron-dependent evolution by gene duplication as a mechanism for the generation of the ancestral MHC gene. Rather, intron insertion into a previously non-split ancestral MHC rod gene consisting of multiple tandemly arranged 28-residue-encoding repeats, or convergent evolution of an originally non-repetitive ancestral MHC rod gene must account for the observed structure of the rod-encoding portion of present-day MHC genes.     

Gene duplication,allelic diversity,selection processes and adaptive value of MHC class II <Emphasis Type="Italic">DRB</Emphasis> genes of the bank vole, <Emphasis Type="Italic">Clethrionomys glareolus</Emphasis>

The generation and maintenance of allelic polymorphism in genes of the major histocompatibility complex (MHC) is a central issue in evolutionary genetics. Recently, the focus has changed from ex situ to in situ populations to understand the mechanisms that determine adaptive MHC polymorphism under natural selection. Birth-and-death evolution and gene conversion events are considered to generate sequence diversity in MHC genes, which subsequently is maintained by balancing selection through parasites. The ongoing arms race between the host and parasites leads to an adaptive selection pressure upon the MHC, evident in high rates of non-synonymous vs synonymous substitution rates. We characterised the MHC class II DRB exon 2 of free living bank voles, Clethrionomys glareolus by single-strand conformation polymorphism and direct sequencing. Unlike other arvicolid species, the DRB locus of the bank vole is at least quadruplicated. No evidence for gene conversion events in the Clgl-DRB sequences was observed. We found not only high allelic polymorphism with 26 alleles in 36 individuals but also high rates of silent polymorphism. Exceptional for MHC class II genes is a purifying selection pressure upon the majority of MHC-DRB sequences. Further, we analysed the association between certain DRB alleles and the parasite burden with gastrointestinal trichostrongyle nematodes Heligmosomum mixtum and Heligmosomoides glareoli and found significant quality differences between specific alleles with respect to infection intensity. Our findings suggest a snapshot in an evolutionary process of ongoing birth-and-death evolution. One allele cluster has lost its function and is already silenced, another is loosing its adaptive value in terms of gastrointestinal nematode resistance, while a third group of alleles indicates all signs of classical functional MHC alleles.     

The nucleotide sequence of the murine I-E beta b immune response gene: evidence for gene conversion events in class II genes of the major histocompatibility complex.

We have determined the DNA sequence of the murine I-E beta b immune response gene of the major histocompatibility complex (MHC) of the C57BL/10 mouse and compared it with the sequence of allelic I-E and non-allelic I-A genes from the d and k haplotypes. The polymorphic exon sequences which encode the first extracellular globular domain of the E beta domain show approximately 8% nucleotide substitutions between the E beta b and E beta d alleles compared with only approximately 2% substitutions for the intron sequences. This suggests that an active mechanism such as micro gene conversion events drive the accumulation of these mutations in the polymorphic exons. The fact that several of the nucleotide changes are clustered supports this hypothesis. The E beta b and E beta k genes show approximately 2-fold fewer nucleotide substitutions than the E beta d/E beta b pair. The A beta bm12, a mutant I-A beta b gene from the C57BL/6 mouse, has been shown to result from three nucleotide changes clustered in a short region of the beta 1 domain, which suggests that a micro gene conversion event caused this mutation. We show here that the E beta b gene is identical to the non-allelic A beta bm12 DNA sequence in the mutated region and suggest, therefore, that the E beta b gene was the donor sequence for this intergenic transfer of genetic information. Diversity in class II MHC genes appears therefore to be generated, at least in part, by the same mechanism proposed for class I genes: intergenic transfer of short DNA regions between non-allelic genes.     

Characterization and Evolution of MHC Class II B Genes in Ardeid Birds

Major histocompatibility complex (MHC) is a multi-gene family that is very suitable to investigate a wide range of open questions in evolutionary ecology. In this study, we characterized two expressed MHC class II B genes (DAB1 and DAB2) in the Grey Heron (Aves: Ardea cinerea). We further developed the primer pairs to amplify and sequence two MHC class II B loci in ten ardeid birds. Phylogenetic analysis revealed that different parts of the genes showed different evolutionary patterns. The exon 2 sequences tended to cluster two gene-specific lineages. In each lineage, exon 2 sequences from several species showed closer relationships than sequences within species, and two shared identical alleles were found between species (Egretta sacra and Nycticorax nycticorax; Egretta garzetta and Bubulcus ibis), supporting the hypothesis of trans-species polymorphism. In contrast, the species-specific intron 2 plus partial exon 3 tree suggested that DAB1 and DAB2 were subject to concerted evolution. GENECONV analyses showed the gene exchange played an important role in the ardeid MHC evolution.     

Distinct evolutionary trajectories of MHC class I and class II genes in Old World finches and buntings

Major histocompatibility complex (MHC) genes code for key proteins of the adaptive immune system, which present antigens from intra-cellular (MHC class I) and extra-cellular (MHC class II) pathogens. Because of their unprecedented diversity, MHC genes have long been an object of scientific interest, but due to methodological difficulties in genotyping of duplicated loci, our knowledge on the evolution of the MHC across different vertebrate lineages is still limited. Here, we compared the evolution of MHC class I and class II genes in three sister clades of common passerine birds, finches (Fringillinae and Carduelinae) and buntings (Emberizidae) using a uniform methodological (genotyping and data processing) approach and uniform sample sizes. Our analyses revealed contrasting evolutionary trajectories of the two MHC classes. We found a stronger signature of pervasive positive selection and higher allele diversity (allele numbers) at the MHC class I than class II. In contrast, MHC class II genes showed greater allele divergence (in terms of nucleotide diversity) and a much stronger recombination (gene conversion) signal. Gene copy numbers at both MHC class I and class II evolved via fluctuating selection and drift (Brownian Motion evolution), but the evolutionary rate was higher at class I. Our study constitutes one of few existing examples, where evolution of MHC class I and class II genes was directly compared using a multi-species approach. We recommend that re-focusing MHC research from single-species and single-class approaches towards multi-species analyses of both MHC classes can substantially increase our understanding MHC evolution in a broad phylogenetic context.Subject terms: Molecular evolution, Immunogenetics     

Non-neutral evolution of the major histocompatibility complex class II gene DRB1 in the sac-winged bat Saccopteryx bilineata

The immune genes of the major histocompatibility complex (MHC) are classical examples for high levels of genetic diversity and non-neutral evolution. This is particularly true for the regions containing the antigen-binding sites as, for instance, in the exon 2 of the MHC class II gene DRB. We surveyed, for the first time in the order Chiroptera, the genetic diversity within this exon in the sac-winged bat Saccopteryx bilineata. We detected 11 alleles among 85 bats, of which 79 were sampled in one population. Pairwise comparisons revealed that interallelic sequence differences ranged between 3 and 22%, although nucleotide substitutions were not evenly distributed along the exon sequence. This was most probably the result of intragenic recombination. High levels of sequence divergence and significantly more nonsynonymous than synonymous substitutions (d(N)/d(S)>1) suggest long-term balancing selection. Thus, the data are consistent with the hypothesis that recombination gives rise to new alleles at the DRB locus of the sac-winged bat, and these are maintained in the population through balancing selection. In this respect, the sac-winged bat closely resembles other mammalian species.     

Contribution of gene conversion in the evolution of the human β-like globin gene family

Gene conversion is referred to as one of two types of mechanisms known to act on gene families, mainly to maintain their sequence homogeneity or, in certain cases, to produce sequence diversity. The concept of gene conversion was established 20 years ago by researchers working with fungi. A few years later, gene conversion was also observed in the human genome, i.e. the γ-globin locus. The aim of this article is to emphasize the role of genetic recombination, particularly of gene conversion, in the evolution of the human β-like globin genes and further to summarize its contribution to the convergent evolution of the fetal globin genes. Finally, this article attempts to re-examine the origin and spread of specific mutations of the β-globin cluster, such as the sickle cell or β-thalassemia mutations, on the basis of repeated gene conversion events. Received: 13 February 1997 / Accepted: 15 May 1998     

Coevolution of PERB11 (MIC) and HLA Class I Genes with HERV-16 and Retroelements by Extended Genomic Duplication

The recent availability of genomic sequence information for the class I region of the MHC has provided an opportunity to examine the genomic organization of HLA class I (HLAcI) and PERB11/MIC genes with a view to explaining their evolution from the perspective of extended genomic duplications rather than by simple gene duplications and/or gene conversion events. Analysis of genomic sequence from two regions of the MHC (the alpha- and beta-blocks) revealed that at least 6 PERB11 and 14 HLAcI genes, pseudogenes, and gene fragments are contained within extended duplicated segments. Each segment was searched for the presence of shared (paralogous) retroelements by RepeatMasker in order to use them as markers of evolution, genetic rearrangements, and evidence of segmental duplications. Shared Alu elements and other retroelements allowed the duplicated segments to be classified into five distinct groups (A to E) that could be further distilled down to an ancient preduplication segment containing a HLA and PERB11 gene, an endogenous retrovirus (HERV-16), and distinctive retroelements. The breakpoints within and between the different HLAcI segments were found mainly within the PERB11 and HLA genes, HERV-16, and other retroelements, suggesting that the latter have played a major role in duplication and indel events leading to the present organization of PERB11 and HLAcI genes. On the basis of the features contained within the segments, a coevolutionary model premised on tandem duplication of single and multipartite genomic segments is proposed. The model is used to explain the origins and genomic organization of retroelements, HERV-16, DNA transposons, PERB11, and HLAcI genes as distinct segmental combinations within the alpha- and beta-blocks of the human MHC. Received: 5 December 1998 / Accepted: 27 January 1999     

Retained orthologous relationships of the MHC Class I genes during euteleost evolution

Major histocompatibility complex (MHC) class I molecules play a pivotal role in immune defense system, presenting the antigen peptides to cytotoxic CD8+ T lymphocytes. Most vertebrates possess multiple MHC class I loci, but the analysis of their evolutionary relationships between distantly related species has difficulties because genetic events such as gene duplication, deletion, recombination, and/or conversion have occurred frequently in these genes. Human MHC class I genes have been conserved only within the primates for up to 46-66 My. Here, we performed comprehensive analysis of the MHC class I genes of the medaka fish, Oryzias latipes, and found that they could be classified into four groups of ancient origin. In phylogenetic analysis using these genes and the classical and nonclassical class I genes of other teleost fishes, three extracellular domains of the class I genes showed quite different evolutionary histories. The α1 domains generated four deeply diverged lineages corresponding to four medaka class I groups with high bootstrap values. These lineages were shared with salmonid and/or other acanthopterygian class I genes, unveiling the orthologous relationships between the classical MHC class I genes of medaka and salmonids, which diverged approximately 260 Ma. This suggested that the lineages must have diverged in the early days of the euteleost evolution and have been maintained for a long time in their genome. In contrast, the α3 domains clustered by species or fish groups, regardless of classical or nonclassical gene types, suggesting that this domain was homogenized in each species during prolonged evolution, possibly retaining the potential for CD8 binding even in the nonclassical genes. On the other hand, the α2 domains formed no apparent clusters with the α1 lineages or with species, suggesting that they were diversified partly by interlocus gene conversion, and that the α1 and α2 domains evolved separately. Such evolutionary mode is characteristic to the teleost MHC class I genes and might have contributed to the long-term conservation of the α1 domain.     

Characterization of a divergent non-classical MHC class I gene in sharks

Sharks are the most ancient group of vertebrates known to possess members of the major histocompatibility complex (MHC) gene family. For this reason, sharks provide a unique opportunity to gain insight into the evolution of the vertebrate immune system through comparative analysis. Two genes encoding proteins related to the MHC class I gene family were isolated from splenic cDNA derived from spiny dogfish shark ( Squalus acanthias). The genes have been designated MhcSqac-UAA*01 and MhcSqac-UAA*NC1. Comparative analysis demonstrates that the Sqac-UAA*01 protein sequence clusters with classical MHC class I of several shark species and has structural elements common to most classical MHC class I molecules. In contrast, Sqac-UAA*NC1 is highly divergent from all vertebrate classical MHC class I proteins, including the Sqac-UAA *01 sequence and those of other shark species. Although Sqac-UAA*NC1 is clearly related to the MHC class I gene family, no orthologous genes from other species were identified due to the high degree of sequence divergence. In fact, the Sqac NC1 protein sequence is the most divergent MHC class-I-like protein identified thus far in any shark species. This high degree of divergence is similar in magnitude to some of the MHC class-I-related genes found in mammals, such as MICA or CD1. These data support the existence of a class of highly divergent non-classical MHC class I genes in the most primitive vertebrates known to possess homologues of the MHC and other components of the adaptive immune system.