Trans-species polymorphism of the Mhc class II DRB-like gene in banded penguins (genus Spheniscus)

The Major Histocompatibility Complex (Mhc) class II DRB locus of vertebrates is highly polymorphic and some alleles may be shared between closely related species as a result of balancing selection in association with resistance to parasites. In this study, we developed a new set of PCR primers to amplify, clone, and sequence overlapping portions of the Mhc class II DRB-like gene from the 5′UTR end to intron 3, including exons 1, 2, and 3 and introns 1 and 2 in four species (20 Humboldt, six African, five Magellanic, and three Galapagos penguins) of penguin from the genus Spheniscus (Sphe). Analysis of gene sequence variation by the neighbor-joining method of 21 Sphe sequences and 20 previously published sequences from four other penguin species revealed overlapping clades within the Sphe species, but species-specific clades for the other penguin species. The overlap of the DRB-like gene sequence variants between the four Sphe species suggests that, despite their allopatric distribution, the Sphe species are closely related and that some shared DRB1 alleles may have undergone a trans-species inheritance because of balancing selection and/or recent rapid speciation. The new primers and PCR assays that we have developed for the identification of the DRB1 DNA and protein sequence variations appear to be useful for the characterization of the molecular evolution of the gene in closely related Penguin species and might be helpful for the assessment of the genetic health and the management of the conservation and captivity of these endangered species. The nucleotide sequence and amino acid sequence data reported in this paper have been submitted to the DDBJ database and have been assigned the accession numbers AB301478, AB301944–AB301950, AB302087–AB302090, AB302190–AB302192, AB302843, AB302844, and AB303942–AB303945.     

Low MHC variation in the endangered Galápagos penguin (<Emphasis Type="Italic">Spheniscus mendiculus</Emphasis>)

The major histocompatibility complex (MHC) is one of the most polymorphic regions of the genome, likely due to balancing selection acting to maintain alleles over time. Lack of MHC variability has been attributed to factors such as genetic drift in small populations and relaxed selection pressure. The Galápagos penguin (Spheniscus mendiculus), endemic to the Galápagos Islands, is the only penguin that occurs on the equator. It relies upon cold, nutrient-rich upwellings and experiences severe population declines when ocean temperatures rise during El Niño events. These bottlenecks, occurring in an already small population, have likely resulted in reduced genetic diversity in this species. In this study, we used MHC class II exon 2 sequence data from a DRB1-like gene to characterize the amount of genetic variation at the MHC in 30 Galápagos penguins, as well as one Magellanic penguin (S. magellanicus) and two king penguins (Aptenodytes patagonicus), and compared it to that in five other penguin species for which published data exist. We found that the Galápagos penguin had the lowest MHC diversity (as measured by number of polymorphic sites and average divergence among alleles) of the eight penguin species studied. A phylogenetic analysis showed that Galápagos penguin MHC sequences are most closely related to Humboldt penguin (Spheniscus humboldti) sequences, its putative sister species based on other loci. An excess of non-synonymous mutations and a pattern of trans-specific evolution in the neighbor-joining tree suggest that selection is acting on the penguin MHC.     

Mhc class I genes of the cichlid fish Oreochromis niloticus

In terms of number of species, perciform (perch-like) fishes are one of the most diversified groups of modern vertebrates. Within this group, the family Cichlidae is best known for its spectacular adaptive radiation in the great lakes of East Africa. The molecular tool kit used in the study of this radiation includes the major histocompatibility complex (Mhc) genes. To refine this tool, information about the organization of the Mhc regions is badly needed. In this study, the first step was taken toward providing such information for the Mhc class one regions of Oreochromis niloticus, a representative species of the tilapiine branch of the Cichlidae, for which good bacterial artificial chromosome library is available. Screening of the library with class I gene probes led to the identification and isolation of 31 class-I-positive clones. Sequencing of one of these clones and partial characterization of the remaining clones for the presence of class I exons resulted in the construction of two contigs representing the class I region of this species as well as identification of seven additional class-I-positive singleton clones. The O. niloticus genome was shown to contain at least 28 class I genes or gene fragments. The shorter of the two contigs was approximately 330 kb long and contained eight class I genes/gene fragments; the longer contig encompassed 1,200 kb of sequence and contained minimally 17 class I genes/gene fragments; three additional class I genes were found to be borne by a clone that might be part of the shorter contig. Electronic supplementary material Supplementary material is available in the online version of this article at and is accessible for authorized users. This work had been carried out in part at the Max-Planck-Institut für Biologie, Abteilung Immungenetik, Tübingen, Germany (A.S., R.D., N.T., S.S., and J.K.). The sequences reported in this paper have been deposited in the GenBank database (accession nos. AB270803–AB270897).     

Sequence analysis of a polymorphic Mhc class II gene in Pacific salmon

Polymorphism of the nucleotide sequences encoding 149 amino acids of linked major histocompatibility complex (Mhc) class II 131 and 132 peptides, and of the intervening intron (548–773 base pairs), was examined within and among seven Pacific salmon (Oncorhynchus) species. Levels of nucleotide diversity were higher for theB1 sequence than forB2 or the intron in comparisons both within and between species. For the codons of the peptide binding region of the BI sequence, the level of nonsynonymous nucleotide substitution (d_N) exceeded the level of synonymous substitution (d_S) by a factor of ten for within-species comparisons, and by a factor of four for between-species comparisons. The excess of d_N indicates that balancing selection maintains diversity at this salmonidMhc class II locus, as is common forMhc loci in other vertebrates. Levels of nucleotide diversity for both the exon and intron sequences were greater among than within species, and there were numerous species-specific nucleotides present in both the coding and noncoding regions. Thus, neighbor-joining analysis of both the intron and exon regions provided phylogenies in which the sequences clustered strongly by species. There was little evidence of shared ancestral (trans-species) polymorphism in the exon phylogeny, and the intron phylogeny depicted standard relationships among the Pacific salmon species. The lack of shared allelicB1 lineages in these closely related species may result from severe bottlenecks that occurred during speciation or during the ice ages that glaciated the rim of the north Pacific Ocean approximately every 100 000 years in the Pleistocene.The nucleotide sequence data reported in this paper have been submitted to the GenBank nucleotide sequence database and have been assigned the accession numbers U34692-U34720     

Molecular cloning of major histocompatibility complex class II B gene cDNA from the Bengalese finch Lonchura striata

The only avian major histocompatibility complex (Mhc) genes thus far identified are from species of the relatively small order of Galliformes, while by far the largest order of Passeriformes (songbirds), containing some 60% of extant bird species, has not been studied at all in this regard. The Galliformes emerged more than 55 million years (my) ago, the Passeriformes some 25 my ago. Because of the potential for the use of Mhc genes as markers in the study of songbird populations, an attempt was made to clone class II B genes of a passeriform species, the Bangalese finch Lonchura striata acuticauda. Using a set of primers designed on the basis of known sequences, a probe corresponding to part of exon II was obtained by the polymerase chain reaction. The probe was then used to screen a Bengalese finch cDNA library and to isolate and sequence two nearly full-length clones. The sequences reveal the presence of one presumably functional class II B locus in this bird species.     

Convergent evolution of major histocompatibility complex molecules in humans and New World monkeys

  In both Old World and New World monkeys Mhc-DRB sequences have been found which resemble human DRB1*03 and DRB3 genes in their second exon. The resemblance is shared sequence motifs and clustering of the genes or the encoded proteins in phylogenetic trees. This similarity could be due to common ancestry, convergence at the molecular level, or chance. To test which of these three explanations applies, we sequenced segments of New World monkey and macaque genes which encompass the entire second exon and large parts of both flanking introns. The test strongly supports the monophyly of New World monkey DRB intron sequences. The phylogenies of introns 1 and 2 from DRB1*03-like and DRB3-like genes are congruent, but both are incongruent with the exon 2-based phylogeny. The matching of intron 1- and intron 2-based phylogenies with each other suggests that reciprocal recombination has not played a major role in exon 2 evolution. Statistical comparisons of exon 2 from different DRB1*03 and DRB3 lineages indicate that it was neither gene conversion (descent), nor chance, but molecular convergence that has shaped their characteristic motifs. The demonstration of convergence in anthropoid Mhc-DRB genes has implications for the classification, age, and mechanism of generation of DRB allelic lineages. Received: 30 August 1999 / Revised: 19 October 1999     

Evolutionary relationship ofHLA-DRB genes inferred from intron sequences

The major histocompatibility complex (Mhc) consists of class I and class II genes. In the humanMhc (HLA) class II genes, nineDRB loci have been identified. To elucidate the origin of these duplicated loci and allelic divergences at the most polymorphicDRBI locus, introns 4 and 5 as well as the 3′ untranslated region (altogether approximately 1,000 base pairs) of sevenHLA-DRB loci, threeHLA-DRBI alleles, and nine nonhuman primateDRB genes were examined. It is shown that there were two major diversification events inHLA-DRB genes, each involving gene duplications and allelic divergences. Approximately 50 million years (my) ago,DRBI ^*04 and an ancestor of theDRB1 ^*03 cluster (DRBI ^*03, DRBI^*15, andDRB3) diverged from each other andDRB5, DRB7, DRB8, and an ancestor of theDRB2 cluster (DRB2, DRB4, andDRB6) arose by gene duplication. Later, about 25 my ago,DRBI ^*15 diverged fromDRBI*03, andDRB3 was duplicated fromDRBI ^*03. Then, some 20 my ago, the lineage leading to theDRB2 cluster produced two new loci,DRB4 andDRB6. TheDRBI ^*03 andDRBI ^*04 allelic lineages are extraordinarily old and have persisted longer than some duplicated genes. The orthologous relationships ofDRB genes between human and Old World monkeys are apparent, but those between Catarrhini and New World monkeys are equivocal because of a rather rapid expansion and contraction of primateDRB genes by duplication and deletion. Correspondence to: Y. Satta     

Molecular evolutionary characterization of an Activator (Ac)-like transposable element sequence from pearl millet (Pennisetum glaucum) (Poaceae)

We present data on the evolution of the Ac/Ds family of transposable elements in select grasses (Poaceae). A defective Ac-like element was cloned from a DNA library of the grass Pennisetum glaucum (pearl millet) and its entire 4531 bp sequence has been determined. When the pearl millet Ac-like sequence is aligned with the maize Ac sequence, it is found that there is approximately 70% DNA similarity in the central region spanning most of maize Ac exon II and all of exon III. In addition, there are two smaller regions of similarity at the Ac terminii. Besides these three major structural similarities, Pennisetum Ac has two large regions, one 5 and one 3, that show little similarity to Zea Ac. Furthermore, most of the sequences corresponding to intron II in maize Ac are absent in pearl millet Ac. Kimura's evolutionary distance between the central region of maize and pearl millet Ac sequences is estimated to be 0.429±0.020 nucleotide substitutions per site. This value is not significantly different from the average number of synonymous substitutions for coding regions of the Adh1 gene between maize and pearl millet, which is 0.395±0.051 nucleotide substitutions per site. If we assume Ac and Adh1 divergence times are equivalent between maize and pearl millet, then the above calculations suggest Ac-like sequences have probably not been strongly constrained by natural selection. Conserved DNA and amino acid sequence motifs are also examined. The level of DNA sequence divergence between maize and pearl millet Ac sequences, the estimated date when maize and pearl millet diverged (25–40 million years ago), coupled with their reproductive isolation/lack of current genetic exchange, all support the theory that Ac-like sequences have not been recently introduced into pearl millet from maize. Instead, Ac-like sequences were probably present in the progenitor of maize and pearl millet and have thus existed in the grasses for at least 25 million years.     

Molecular study of Mhc-DRB in wild chacma baboons reveals high variability and evidence for trans-species inheritance

The MHC class II genes of many primate species were investigated extensively in recent years. However, while Mhc-DRB genes were studied in Old World monkeys such as rhesus macaques, the Mhc-DRB of baboons was only studied in a limited way. Because of their close anatomical and physiological relationship to humans, baboons are often used as models for reproduction and transplantation research. Baboons are also studied as a model species in behavioural ecology. Thus, identification of MHC genes would provide a foundation for studies of Mhc, biology and behaviour. Here, we describe the use of PCR, cloning, denaturing gradient gel electrophoresis (DGGE) and sequencing to identify Mhc-DRB sequences in wild chacma baboons (Papio ursinus). We amplified the highly variable second exon of baboon Mhc-DRB sequences using generic DRB primers. To validate and optimize the DGGE protocol, four DNA samples were initially studied using cloning and sequencing. Clones were screened using a novel RFLP approach to increase the number of clones identified for each individual. Results from cloning and sequencing were used to optimise DGGE conditions for Mhc-DRB genotyping of the remaining study subjects. Using these techniques, we identified 16 Paur-DRB sequences from 30 chacma baboons. On the basis of phylogenetic tree analyses, representatives of the Mhc-DRB1 and Mhc-DRB5 loci, and 13 different DRB lineages were identified. Evidence for trans-species inheritance of some Mhc-DRB sequences comes from high identity between the new Paur-DRB sequences and sequences from Papio cynocephalus, Macaca mulatta and possibly Galago moholi.Nucleotide sequence data reported are available in the GenBank/EMBL/DDBJ databases under the accession numbers DQ339722–DQ339737.     

Major histocompatibility complex class II A gene polymorphism in the striped bass

Adaptions of the polymerase chain reaction were used to isolate cDNA sequences encoding the Major histocompatibility complex(Mhc) class II A gene(s) of the striped bass (Morone saxatilis). Four complete Mhc class II A genes were cloned and sequenced from a specimen originating in the Roanoke River, North Carolina, and another three A genes from a specimen originating from the Santee-Cooper Reservoir, South Carolina, identifying a total of seven unique sequences. The sequence suggests the presence of at least two Mhc class II A loci. The extensive sequence variability observed between the seven different Mhc class II clones was concentrated in the 1 encoding domain. The encoded 2, transmembrane, and cytoplasmic regions of all seven striped genes correlated well with those of known vertebrate Mhc class II proteins. Overall, the striped bass sequences showed greatest similarity to the Mhc class II A genes of the zebrafish. Southern blot analysis demonstrated extensive polymorphism in the Mhc class II A genes in members of a Roanoke river-caught population of striped bass versus a lesser degree of polymorphism in an aquacultured Santee-Cooper population of striped bass.The nucleotide sequence data reported in this paper have been submitted to the GenBank nucleotide sequence database and have been assigned the accession numbers (Mosa-A-S5) L35062, (Mosa-A-S8) L35066, (Mosa-A-R7) L35067, and (Mosa-A-S7) L35072 L35066, (Mossa-A-R7) L35067, and (Mosa-A-S7) L35072     

Songbird Genomics: Analysis of 45 kb Upstream of a Polymorphic Mhc Class II Gene in Red-Winged Blackbirds (Agelaius phoeniceus)

Here we present the sequence of a 45 kb cosmid containing a previously characterized poly-morphic Mhc class II B gene (Agph-DAB1) from the red-winged blackbird (Agelaius phoeniceus). We compared it with a previously sequenced cosmid from this species, revealing two regions of 7.5 kb and 13.0 kb that averaged greater than 97% similarity to each another, indicating a very recent shared duplication. We found 12 retroelements, including two chicken repeat 1 (CR1) elements, constituting 6.4% of the sequence and indicating a lower frequency of retroelements than that found in mammalian genomic DNA. Agph-DAB3, a new class II B gene discovered in the cosmid, showed a low rate of polymorphism and may be functional. In addition, we found a Mhc class II B gene fragment and three genes likely to be functional (encoding activin receptor type II, a zinc finger, and a putative γ-filamin). Phylogenetic analysis of exon 2 alleles of all three known blackbird Mhc genes indicated strong clustering of alleles by locus, implying that large amounts of interlocus gene conversion have not occurred since these genes have been diverging. Despite this, interspecific comparisons indicate that all three blackbird Mhc genes diverged from one another less than 35 million years ago and are subject to concerted evolution in the long term. Comparison of blackbird and chicken Mhc promoter regions revealed songbird promoter elements for the first time. The high gene density of this cosmid confirms similar findings for the chicken Mhc, but the segment duplications and diversity of retroelements resembles mammalian sequences.     

Characterization of the horse (<Emphasis Type="Italic">Equus caballus</Emphasis>)<Emphasis Type="Italic"> IGHA</Emphasis> gene

Nucleotide sequences of the immunoglobulin constant heavy chain genes of the horse have been described for IGHM, IGHG and IGHE genes, but not for IGHA. Here, we provide the nucleotide sequence of the genomic IGHA gene of the horse (Equus caballus), including its secretion region and the transmembrane exon. The equine IGHA gene shows the typical structure of a mammalian IGHA gene, with only three exons, separated by two introns of similar size. The hinge exon is located at the 5 end of the CH2 exon and encodes a hinge region of 11 amino acids, which contains five proline residues. The coding nucleotide sequence of the secreted form of the equine IGHA gene shares around 72% identity with the human IGHA1 and IGHA2 genes, as well as the bovine, ovine, porcine and canine IGHA genes, without distinct preference for any of these species. The same species also cluster together in a phylogenetic tree of the IGHA coding regions of various mammals, whereas rodent, rabbit, marsupial and monotreme IGHA genes each build a separate cluster.The nucleotide sequences reported in this paper have been assigned the EMBL/GenBank accession numbers AY247966 and AY351982     

Evolution of Ac and Dsl elements in select grasses (Poaceae)

We present data on evolution of the Ac/Ds family of transposable elements in select grasses (Poaceae). An Ac-like element was cloned from a DNA library of the grass Pennisetum glaucum (pearl millet) and 2387 bp of it have been sequenced. When the pearl millet Ac-like sequence is aligned with the corresponding region of the maize Ac sequence, it is found that all sequences corresponding to intron II in maize Ac are absent in pearl millet Ac. Kimura's evolutionary distance between maize and pearl millet Ac sequences is estimated to be 0.429±0.020 nucleotide substitutions per site. This value is not significantly different from the average number of synonymous substitutions for coding regions of the Adh1 gene between maize and pearl millet, which is 0.395±0.051 nucleotide substitutions per site. If we can assume Ac and Adh1 divergence times are equivalent between maize and pearl millet, then the above calculations suggest Ac-like sequences have probably not been strongly constrained by natural selection. The level of DNA sequence divergence between maize and pearl millet Ac sequences, the estimated date when maize and pearl millet diverged (25–40 million years ago), coupled with their reproductive isolation/lack of current genetic exchange, all support the theory that Ac-like sequences have not been recently introduced into pearl millet from maize. Instead, Ac-like sequences were probably present in the progenitor of maize and pearl millet, and have thus existed in the grasses for at least 25 million years. Ac-like sequences may be widely distributed among the grasses. We also present the first 2 Dsl controlling element sequences from teosinte species: Zea luxurians and Zea perennis. A total of 10 Dsl elements had previously been sequenced from maize and a distant maize relative, Tripsacum. When a maximum likelihood network of genetic relationships is constructed for all 12 sequenced Dsl elements, the 2 teosinte Dsl elements are as distant from most maize Dsl elements and from each other, as the maize Dsl elements are from one another. Our new teosinte sequence data support the previous conclusion that Dsl elements have been accumulating mutations independently since maize and Tripsacum diverged. We present a scenario for the origin of Dsl elements.     

Evolutionary history of an MHC gene in two leporid species: characterisation of <Emphasis Type="Italic">Mhc-DQA</Emphasis> 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. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users. The nucleotide sequence data reported in this paper have been submitted to Genbank and have been assigned the accession numbers FJ225335–FJ225346.     

Evolutionary relationships among the primate Mhc-DQA1 and DQA2 alleles

The variation of the Mhc-DQA1 and DQA2 loci of ten different primate species (hominoids and Old World monkeys) was studied in order to obtain an insight in the processes that generate polymorphism of major histocompatibility complex (Mhc) class II genes and to establish the evolutionary relationships of their alleles. To that end nucleotide sequences of 36 Mhc class II DQA1 and seven DQA2 second exons were determined and phylogenetic trees that illustrate their evolutionary relationships were constructed. We demonstrate the existence of four primate Mhc-DQA1 allele lineages, two of which probably existed before the separation of the ancestors of the hominoids and Old World monkeys (approximately 22–28 million years ago). Mhc-DQA2 sequences were found only in the hominoid species and showed little diversity. We found no evidence for a major contribution of recombinational events to the generation of allelic diversity of the primate Mhc-DQA1 locus. Instead, our data suggest that the primate Mhc-DQA1 and DQA2 loci are relatively stable entities that mutated primarily as a result of point mutations.The nucleotide sequence data reported in this paper have been submitted to the GenBank nucleotide sequence database and have been assigned the accession numbers M76186-M76229.     

Evolution of the transposable element Uhu in five species of Hawaiian Drosophila

The complete DNA sequence of three independent isolates of Uhu, a member of the Tc1-like class of transposable elements from D. heteroneura (Uhu-1, Uhu-3, and Uhu-4), has been determined. These isolates have between 95 and 96.4% nucleotide sequence identity indicating that Uhu is well conserved within this species. A comparison of the DNA sequences of Uhu and the D. melanogaster Hb1 transposable element shows that the nucleotide substitution rate for Uhu is comparable to the synonymous rate for the Adh gene in these species. Uhu has been identified in four other species of endemic Hawaiian Drosophila, D. silvestris, D. differens, D. planitibia and D. picticornis, and nine Uhu elements were isolated from genomic libraries of these four species. A 444 base pair region from within the coding region of the Uhu element, with well conserved ends, was amplified by the polymerase chain reaction and used for sequence comparison of elements from different species. The analysis of the sequence similarities between the elements within and between the species shows a grouping of the two pairs of most closely related species (D. heteroneura and D. silvestris, and D. differens and D. planitibia), but shows a much larger variation within the most recently diverged species (D. heteroneura and D. silvestris) than expected. There are extensive nucleotide substitutions and deletions in the Uhu elements from D. picticornis showing that they are degenerating and being lost in this species. These observations indicate that the Uhu element has been transmitted vertically and that transposition may have been activated at the time of formation of each species as it colonized the newly formed islands of the Hawaiian archipelago.     

Dynamics of Mhc evolution in birds and crocodilians: amplification of class II genes with degenerate primers

Genes of the major histocompatibility complex (Mhc) are the most polymorphic functional loci in mammalian populations, but little is known of Mhc variability in natural populations of nonmammalian vertebrates. To help extend such studies to birds and relatives, we present a pair of degenerate primers that amplify polymorphic segments of one chain (the β chain) of the class II genes from the major histocompatibility complex (Mhc) of archosaurs (birds + crocodilians). The primers target two conserved regions lying within portions of the antigen-binding site (ABS) encoded by the second exon and amplify multiple genes from both genomic DNA and cDNA. The pattern of nucleotide substitution in ABS codons of 51 sequences amplified and cloned from five species of passerine birds and an alligator (Alligator mississippiensis) indicates that archosaurian class II β genes are subject to selective forces similar to those operating in mammalian populations. Hybridization of a genomic clone generated by the primers revealed highly polymorphic bands in a sample of Florida scrub jays (Aphelocoma coerulescens coerulescens). Because the primers amplify only part of the ABS from multiple class II genes, they will be useful primarily for generating species specific clones, thereby providing a critical inroad to more detailed structural and evolutionary studies.