Blueprint for a high-performance biomaterial: full-length spider dragline silk genes

Spider dragline (major ampullate) silk outperforms virtually all other natural and manmade materials in terms of tensile strength and toughness. For this reason, the mass-production of artificial spider silks through transgenic technologies has been a major goal of biomimetics research. Although all known arthropod silk proteins are extremely large (>200 kiloDaltons), recombinant spider silks have been designed from short and incomplete cDNAs, the only available sequences. Here we describe the first full-length spider silk gene sequences and their flanking regions. These genes encode the MaSp1 and MaSp2 proteins that compose the black widow's high-performance dragline silk. Each gene includes a single enormous exon (>9000 base pairs) that translates into a highly repetitive polypeptide. Patterns of variation among sequence repeats at the amino acid and nucleotide levels indicate that the interaction of selection, intergenic recombination, and intragenic recombination governs the evolution of these highly unusual, modular proteins. Phylogenetic footprinting revealed putative regulatory elements in non-coding flanking sequences. Conservation of both upstream and downstream flanking sequences was especially striking between the two paralogous black widow major ampullate silk genes. Because these genes are co-expressed within the same silk gland, there may have been selection for similarity in regulatory regions. Our new data provide complete templates for synthesis of recombinant silk proteins that significantly improve the degree to which artificial silks mimic natural spider dragline fibers.     

Long-term balancing selection at the west nile virus resistance gene, Oas1b, maintains transspecific polymorphisms in the house mouse

Oligoadenylate synthetases (OASs) are interferon-inducible enzymes that participate in the first line of defense against a wide range of viral infection. Recent studies have determined that Oas1b, a member of the OAS gene family in the house mouse (Mus musculus), provides specific protection against flavivirus infection (e.g., West Nile virus, dengue fever virus, and yellow fever virus). We characterized the nucleotide sequence variation in coding and noncoding regions of the Oas1b gene for a large number of wild-derived strains of M. musculus and related species. Our sequence analyses determined that this gene is one of the most polymorphic genes ever described in any mammal. The level of variation in noncoding regions of Oas1b is an order of magnitude higher than the level reported for other regions of the mouse genome and is significantly different from the level of intraspecific variation expected under neutrality. Furthermore, a phylogenetic analysis of intronic sequences demonstrated that Oas1b alleles are ancient and that their divergence predates several speciation events, resulting in transspecific polymorphisms. The amino acid sequence of Oas1b is also extremely variable, with 1 out of 7 amino acid positions being polymorphic within M. musculus. Oas1b alleles are comparatively more divergent at synonymous positions than most autosomal genes and the ratio of nonsynonymous to synonymous substitution is remarkably high, suggesting that positive selection has been acting on Oas1b. The ancestry of Oas1b polymorphisms and the high level of amino acid polymorphisms strongly suggest that the allelic variation at Oas1b has been maintained in mouse populations by long-term balancing selection.     

Recombination and selection at Brassica self-incompatibility loci

In Brassica species, self-incompatibility is controlled genetically by haplotypes involving two known genes, SLG and SRK, and possibly an as yet unknown gene controlling pollen incompatibility types. Alleles at the incompatibility loci are maintained by frequency-dependent selection, and diversity at SLG and SRK appears to be very ancient, with high diversity at silent and replacement sites, particularly in certain "hypervariable" portions of the genes. It is important to test whether recombination occurs in these genes before inferences about function of different parts of the genes can be made from patterns of diversity within their sequences. In addition, it has been suggested that, to maintain the relationship between alleles within a given S-haplotype, recombination is suppressed in the S-locus region. The high diversity makes many population genetic measures of recombination inapplicable. We have analyzed linkage disequilibrium within the SLG gene of two Brassica species, using published coding sequences. The results suggest that intragenic recombination has occurred in the evolutionary history of these alleles. This is supported by patterns of synonymous nucleotide diversity within both the SLG and SRK genes, and between domains of the SRK gene. Finally, clusters of linkage disequilibrium within the SLG gene suggest that hypervariable regions are under balancing selection, and are not merely regions of relaxed selective constraint.     

Neutral evolution of the sex-determining gene transformer in Drosophila

The amino acid sequence of the transformer (tra) gene exhibits an extremely rapid rate of evolution among Drosophila species, although the gene performs a critical step in sex determination. These changes in amino acid sequence are the result of either natural selection or neutral evolution. To differentiate between selective and neutral causes of this evolutionary change, analyses of both intraspecific and interspecific patterns of molecular evolution of tra gene sequences are presented. Sequences of 31 tra alleles were obtained from Drosophila americana. Many replacement and silent nucleotide variants are present among the alleles; however, the distribution of this sequence variation is consistent with neutral evolution. Sequence evolution was also examined among six species representative of the genus Drosophila. For most lineages and most regions of the gene, both silent and replacement substitutions have accumulated in a constant, clock-like manner. In exon 3 of D. virilis and D. americana we find evidence for an elevated rate of nonsynonymous substitution, but no statistical support for a greater rate of nonsynonymous relative to synonymous substitutions. Both levels of analysis of the tra sequence suggest that, although the gene is evolving at a rapid pace, these changes are neutral in function.     

HLA DNA sequence variation among human populations: molecular signatures of demographic and selective events

Molecular differences between HLA alleles vary up to 57 nucleotides within the peptide binding coding region of human Major Histocompatibility Complex (MHC) genes, but it is still unclear whether this variation results from a stochastic process or from selective constraints related to functional differences among HLA molecules. Although HLA alleles are generally treated as equidistant molecular units in population genetic studies, DNA sequence diversity among populations is also crucial to interpret the observed HLA polymorphism. In this study, we used a large dataset of 2,062 DNA sequences defined for the different HLA alleles to analyze nucleotide diversity of seven HLA genes in 23,500 individuals of about 200 populations spread worldwide. We first analyzed the HLA molecular structure and diversity of these populations in relation to geographic variation and we further investigated possible departures from selective neutrality through Tajima's tests and mismatch distributions. All results were compared to those obtained by classical approaches applied to HLA allele frequencies.Our study shows that the global patterns of HLA nucleotide diversity among populations are significantly correlated to geography, although in some specific cases the molecular information reveals unexpected genetic relationships. At all loci except HLA-DPB1, populations have accumulated a high proportion of very divergent alleles, suggesting an advantage of heterozygotes expressing molecularly distant HLA molecules (asymmetric overdominant selection model). However, both different intensities of selection and unequal levels of gene conversion may explain the heterogeneous mismatch distributions observed among the loci. Also, distinctive patterns of sequence divergence observed at the HLA-DPB1 locus suggest current neutrality but old selective pressures on this gene. We conclude that HLA DNA sequences advantageously complement HLA allele frequencies as a source of data used to explore the genetic history of human populations, and that their analysis allows a more thorough investigation of human MHC molecular evolution.     

Analysis of the sequence variations in the Mhc DRB1-like gene of the endangered Humboldt penguin (Spheniscus humboldti)

The Major Histocompatibility Complex (Mhc) genomic region of many vertebrates is known to contain at least one highly polymorphic class II gene that is homologous in sequence to one or other of the human Mhc DRB1 class II genes. The diversity of the avian Mhc class II gene sequences have been extensively studied in chickens, quails, and some songbirds, but have been largely ignored in the oceanic birds, including the flightless penguins. We have previously reported that several penguin species have a high degree of polymorphism on exon 2 of the Mhc class II DRB1-like gene. In this study, we present for the first time the complete nucleotide sequences of exon 2, intron 2, and exon 3 of the DRB1-like gene of 20 Humboldt penguins, a species that is presently vulnerable to the dangers of extinction. The Humboldt DRB1-like nucleotide and amino acid sequences reveal at least eight unique alleles. Phylogenetic analysis of all the available avian DRB-like sequences showed that, of five penguin species and nine other bird species, the sequences of the Humboldt penguins grouped most closely to the Little penguin and the mallard, respectively. The present analysis confirms that the sequence variations of the Mhc class II gene, DRB1, are useful for discriminating among individuals within the same penguin population as well those within different penguin population groups and 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 AB088371–AB088374, AB089199, AB154393–AB154399, and AB162144.     

Genome-wide comparative analysis of simple sequence coding repeats among 25 insect species

We present a detailed genome-scale comparative analysis of simple sequence repeats within protein coding regions among 25 insect genomes. The repetitive sequences in the coding regions primarily represented single codon repeats and codon pair repeats. The CAG triplet is highly repetitive in the coding regions of insect genomes. It is frequently paired with the synonymous codon CAA to code for polyglutamine repeats. The codon pairs that are least repetitive code for polyalanine repeats. The frequency of hexanucleotide and dinucleotide motifs of codon pair repeats is significantly (p<0.001) different in the Drosophila species compared to the non-Drosophila species. However, the frequency of synonymous and non-synonymous codon pair repeats varies in a correlated manner (r(2)=0.79) among all the species. Results further show that perfect and imperfect repeats have significant association with the trinucleotide and hexanucleotide coding repeats in most of these insects. However, only select species show significant association between the numbers of perfect/imperfect hexamers and repeat coding for single amino acid/amino acid pair runs. Our data further suggests that genes containing simple sequence coding repeats may be under negative selection as they tend to be poorly conserved across species. The sequences of coding repeats of orthologous genes vary according to the known phylogeny among the species. In conclusion, the study shows that simple sequence coding repeats are important features of genome diversity among insects.     

Genomic structure of the luciferase gene and phylogenetic analysis in the Hotaria-group fireflies

The luminescent fireflies have species specific flash patterns, being recognized as sexual communication. The luciferase gene is the sole enzyme responsible for bioluminescence. We describe here the complete nucleotide sequence and the exon-intron structure of the luciferase gene of the Hotaria-group fireflies, H. unmunsana, H. papariensis and H. tsushimana. The luciferase gene of the Hotaria-group firefly including the known H. parvula spans 1950 bp and consisted of six introns and seven exons coding for 548 amino acid residues, suggesting highly conserved structure among the Hotaria-group fireflies. Although only one luciferase gene was cloned from H. papariensis, each of the two sequences of the gene was found in H. unmunsana (U1 and Uc) and H. tsushimana (T1 and T2). The amino acid sequence divergence among H. unmunsana, H. papariensis, and H. tsushimana only ranged from zero to three amino acid residues, but H. parvula differed by 10-11 amino acid residues from the other Hotaria-group fireflies, suggesting a divergent relationship of this species. Phylogenetic analysis using the deduced amino acid sequences of the luciferase gene resulted in a monophyletic group in the Hotaria excluding H. parvula, suggesting a close relationship among H. unmunsana, H. papariensis and H. tsushimana. Additionally, we also analyzed the mitochondrial cytochrome oxidase I (COI) gene of the Hotaria-group fireflies. The deduced amino acid sequence of the COI gene of H. unmunsana was identical to that of H. papariensis and H. tsushimana, but different by three positions from H. parvula. In terms of nucleotide sequences of the COI gene, intraspecific sequence divergence was sometimes larger than interspecies level, and phylogenetic analysis placed the three species into monophyletic groups unresolved among them, but excluded H. parvula. In conclusion, our results suggest that H. unmunsana, H. papariensis and H. tsushimana are very closely related or might be an identical species, at least based on the luciferase and COI genes.     

Molecular evolution of two linked genes, Est-6 and Sod, in Drosophila melanogaster.

We have obtained 15 sequences of Est-6 from a natural population of Drosophila melanogaster to test whether linkage disequilibrium exists between Est-6 and the closely linked Sod, and whether natural selection may be involved. An early experiment with allozymes had shown linkage disequilibrium between these two loci, while none was detected between other gene pairs. The Sod sequences for the same 15 haplotypes were obtained previously. The two genes exhibit similar levels of nucleotide polymorphism, but the patterns are different. In Est-6, there are nine amino acid replacement polymorphisms, one of which accounts for the S-F allozyme polymorphism. In Sod, there is only one replacement polymorphism, which corresponds to the S-F allozyme polymorphism. The transversion/transition ratio is more than five times larger in Sod than in Est-6. At the nucleotide level, the S and F alleles of Est-6 make up two allele families that are quite different from each other, while there is relatively little variation within each of them. There are also two families of alleles in Sod, one consisting of a subset of F alleles, and the other consisting of another subset of F alleles, designed F(A), plus all the S alleles. The Sod F(A) and S alleles are completely or nearly identical in nucleotide sequence, except for the replacement mutation that accounts for the allozyme difference. The two allele families have independent evolutionary histories in the two genes. There are traces of statistically significant linkage disequilibrium between the two genes that, we suggest, may have arisen as a consequence of selection favoring one particular sequence at each locus.     

Mutational bias affects protein evolution in flowering plants

Amino acid sequences from several thousand homologous gene pairs were compared for two plant genomes, Oryza sativa and Arabidopsis thaliana. The Arabidopsis genes all have similar G+C (guanine plus cytosine) contents, whereas their homologs in rice span a wide range of G+C levels. The results show that those rice genes that display increased divergence in their nucleotide composition (specifically, increased G+C content) showed a corresponding, predictable change in the amino acid compositions of the encoded proteins relative to their Arabidopsis homologs. This trend was not seen in a "control" set of rice genes that had nucleotide contents closer to their Arabidopsis homologs. In addition to showing an overall difference in the amino acid composition of the homologous proteins, we were also able to investigate the biased patterns of amino acid substitution since the divergence of these two species. We found that the amino acid exchange matrix was highly asymmetric when comparing the High G+C rice genes with their Arabidopsis homologs. Finally, we investigated the possible causes of this biased pattern of sequence evolution. Our results indicate that the biased pattern of protein evolution is the consequence, rather than the cause, of the corresponding changes in nucleotide content. In fact, there is an even more marked asymmetry in the patterns of substitution at synonymous nucleotide sites. Surprisingly, there is a very strong negative correlation between the level of nucleotide bias and the length of the coding sequences within the rice genome. This difference in gene length may provide important clues about the underlying mechanisms.     

Conservation of alternative splicing and genomic organization of the myosin alkali light-chain (Mlc1) gene among Drosophila species

The Mlc1 gene of Drosophila melanogaster encodes two MLC1 isoforms via developmentally regulated alternative pre-mRNA splicing. In larval muscle and tubular and abdominal muscles of adults, all of the six exons are included in the spliced mRNA, whereas, in the fibrillar indirect flight muscle of adult, exon 5 is excluded from the mRNA. We show that this tissue-specific pattern of alternative splicing of the Mlc1 pre-mRNA is conserved in D. simulans, D. pseudoobscura, and D. virilis. Isolation and sequencing of the Mlc1 genes from these three other Drosophila species have revealed that the overall organization of the genes is identical and that the genes have maintained a very high level of sequence identity within the coding region. Pairwise amino acid identities are 94%-99%, and there are no charge changes among the proteins. Total nucleotide divergence within the coding region of the four genes supports the accepted genealogy of these species, but the data indicate a significantly higher rate of amino acid replacement in the branch leading to D. pseudoobscura. A comparison of nucleotide substitutions in the coding portions of exon 5 and exon 6, which encode the alternative carboxyl termini of the two MLC1 isoforms, suggests that exon 5 is subject to greater evolutionary constraints than is exon 6. In addition to the coding sequences, there is significant sequence conservation within the 5' and 3' noncoding DNA and two of the introns, including one that flanks exon 5. These regions are candidates for cis- regulatory elements. Our results suggest that evolutionary constraints are acting on both the coding and noncoding sequences of the Mlc1 gene to maintain proper expression and function of the two MLC1 polypeptides.      

Characterization of cryptic flagellin genes in Shigella boydii and Shigella dysenteriae.

Flagellin (fliC) genes of 12 Shigella boydii and five Shigella dysenteriae strains were characterized. Though these strains are nonmotile, the cryptic fliCSB gene, cloned from S. boydii strain C3, is functional for expression of flagellin. It consists of 1,704 bp, and encodes 568 amino acid residues (57,918 Da). The fliCSD gene from S. dysenteriae strain 16 consists of 1,650 bp encoding 549 amino acid residues (57,591 Da) and contains an IS1 element inserted in its 3' end. The two genes are composed of the 5'-constant, central variable and 3'-constant sequences, like other known fliC genes. The two genes share high homology in nucleotide and amino acid sequences with each other and also with the Escherichia coli fliCE gene, indicating that both genes are closely related to the fliCE gene. Comparison of the central variable sequences of six different fliC genes showed that the fliCSB and fliCSD genes share low homology in amino acid sequence with the other fliC genes, suggesting that they encode antigenic determinants intrinsic to respective subgroups. However, Southern blotting using as probes the central variable sequences of several fliC genes showed that four of 12 S. boydii strains have a fliC gene similar to that of Shigella flexneri, and that among five fliC genes from S. dysenteriae strains, one is similar to that of S. flexneri, two are similar to that of S. boydii, and only one is unique to S. dysenteriae. Some of these variant alleles were verified by immunoblotting with flagellins produced from cloned fliC genes. The presence of variant fliC alleles in S. boydii and S. dysenteriae indicates that subdivision into subgroups does not reflect the ancestral flagella H antigenic relationships. These data will be useful in considering the evolutionary divergence of the Shigella spp..     

Different selective pressures shape the molecular evolution of color vision in chimpanzee and human populations

A population genetic analysis of the long-wavelength opsin (OPN1LW, "red") color vision gene in a global sample of 236 human nucleotide sequences had previously discovered nine amino acid replacement single nucleotide polymorphisms, which were found at high frequencies in both African and non-African populations and associated with an unusual haplotype diversity. Although this pattern of nucleotide diversity is consistent with balancing selection, it has been argued that a recombination "hot spot" or gene conversion within and between X-linked color vision genes alone may explain these patterns. The current analysis investigates a closely related primate with trichromatism to determine whether color vision gene amino acid polymorphism and signatures of adaptive evolution are characteristic of humans alone. Our population sample of 56 chimpanzee (Pan troglodytes) OPN1LW sequences shows three singleton amino acid polymorphisms and no unusual recombination or linkage disequilibrium patterns across the approximately 5.5-kb region analyzed. Our comparative population genetic approach shows that the patterns of OPN1LW variation in humans and chimpanzees are consistent with positive and purifying selection within the two lineages, respectively. Although the complex role of color vision has been greatly documented in primate evolution in general, it is surprising that trichromatism has followed very different selective trajectories even between humans and our closest relatives.     

Comparative analysis of the prion protein (PrP) gene in cetacean species

The partial PrP gene sequence and the deduced protein of eight cetacean species, seven of which have never been reported so far, have been determined in order to extend knowledge of sequence variability of the PrP genes in different species and to aid in speculation on cetacean susceptibility to prions. Both the nucleotide and the deduced amino acid sequences have been analysed in comparison with some of the known mammalian PrPs. Cetacean PrPs present typical features of eutherian PrPs. The PrP gene from the species of the family Delphinidae gave identical nucleic acid sequences, while differences in the PrP gene were found in Balaenopteridae and Ziphidae. The phylogenetic tree resulting from analysis of the cetacean PrP gene sequences, together with reported sequences of some ungulates, carnivores and primates, showed that the PrP gene phylogenesis mirrors the species phylogenesis. The PrP gene of cetaceans is very close to species where natural forms of TSEs are known. From an analysis of the sequences and the phylogenesis of the PrP gene, susceptibility to or occurrence of prion diseases in cetaceans can not be excluded.     

Molecular evolution of the nontandemly repeated genes of the histone 3 multigene family.

In some species, histone gene clusters consist of tandem arrays of each type of histone gene, whereas in other species the genes may be clustered but not arranged in tandem. In certain species, however, histone genes are found scattered across several different chromosomes. This study examines the evolution of histone 3 (H3) genes that are not arranged in large clusters of tandem repeats. Although H3 amino acid sequences are highly conserved both within and between species, we found that the nucleotide sequence divergence at synonymous sites is high, indicating that purifying selection is the major force for maintaining H3 amino acid sequence homogeneity over long-term evolution. In cases where synonymous-site divergence was low, recent gene duplication appeared to be a better explanation than gene conversion. These results, and other observations on gene inactivation, organization, and phylogeny, indicated that these H3 genes evolve according to a birth-and-death process under strong purifying selection. Thus, we found little evidence to support previous claims that all H3 proteins, regardless of their genome organization, undergo concerted evolution. Further analyses of the structure of H3 proteins revealed that the histones of higher eukaryotes might have evolved from a replication-independent-like H3 gene.     

Structural diversity of the classical H-2 genes: K, D, and L.

Twenty-three class I DNA sequences, representing alleles of the H-2K, D, and L loci, were analyzed to assess patterns of nucleotide and amino acid diversity. Comparisons of the allelic and nonallelic sequences revealed locus specificity in regions encoding the leader peptides and the carboxyl-terminal segments of the Ag presenting molecules. Analyses focusing on the sequences that determine the Ag binding domains revealed weak or insignificant allelic associations, a finding that is in sharp contrast to previously observed relationships among the homologous human sequences. The amino acid positions exhibiting high diversity in the encoded glycoproteins in both mice and humans are localized primarily to the Ag binding site. In the mouse, diverse amino acids were positioned similarly in the K and D/L glycoproteins, although in humans, the A and B glycoproteins exhibit distinctive differences in their locations within the Ag binding site. The absence of locus specificity among the sequences that determine the Ag binding domains of the mouse is consistent with the hypothesis that ectopic gene conversion leads to interlocus exchange of class I sequences. Comparable interlocus exchanges among human class I genes have not played a similar role in shaping human A and B sequences. The basis of this difference between mice and humans is not clear. The nature of amino acid substitutions distinguishing class I loci in mice and humans are comparable, and the role of natural selection in determining diversity appears to be similar in the two species.     

Molecular characterization and evolutionary study of spider tubuliform (eggcase) silk protein

As a result of hundreds of millions of years of evolution, orb-web-weaving spiders have developed the use of seven different silks produced by different abdominal glands for various functions. Tubuliform silk (eggcase silk) is unique among these spider silks due to its high serine and very low glycine content. In addition, tubuliform silk is the only silk produced just during a short period of time, the reproductive season, in the spider's life. To understand the molecular characteristics of the proteins composing this silk, we constructed tubuliform-gland-specific cDNA libraries from three different spider families, Nephila clavipes, Argiope aurantia, and Araneus gemmoides. Sequencing of tubuliform silk cDNAs reveals the repetitive architecture of its coding sequence and novel amino acid motifs. The inferred protein, tubuliform spidroin 1 (TuSp1), contains highly homogenized repeats in all three spiders. Amino acid composition comparison of the predicted tubuliform silk protein sequence to tubuliform silk indicates that TuSp1 is the major component of tubuliform silk. Repeat unit alignment of TuSp1 among three spider species shows high sequence conservation among tubuliform silk protein orthologue groups. Sequence comparison among TuSp1 repetitive units within species suggests intragenic concerted evolution, presumably through gene conversion and unequal crossover events. Comparative analysis demonstrates that TuSp1 represents a new orthologue in the spider silk gene family.     

The nucleotide sequence of the mouse immunoglobulin epsilon gene: comparison with the human epsilon gene sequence

We have determined the nucleotide sequence of the immunoglobulin epsilon gene cloned from newborn mouse DNA. The epsilon gene sequence allows prediction of the amino acid sequence of the constant region of the epsilon chain and comparison of it with sequences of the human epsilon and other mouse immunoglobulin genes. The epsilon gene was shown to be under the weakest selection pressure at the protein level among the immunoglobulin genes although the divergence at the synonymous position is similar. Our results suggest that the epsilon gene may be dispensable, which is in accord with the fact that IgE has only obscure roles in the immune defense system but has an undesirable role as a mediator of hypersensitivity. The sequence data suggest that the human and murine epsilon genes were derived from different ancestors duplicated a long time ago. The amino acid sequence of the epsilon chain is more homologous to those of the gamma chains than the other mouse heavy chains. Two membrane exons, separated by an 80-base intron, were identified 1.7 kb 3' to the CH4 domain of the epsilon gene and shown to conserve a hydrophobic portion similar to those of other heavy chain genes. RNA blot hybridization showed that the epsilon membrane exons are transcribed into two species of mRNA in an IgE hybridoma.     

Diversity and duplication of DQB and DRB-like genes of the MHC in baleen whales (suborder: Mysticeti)

The molecular diversity and phylogenetic relationships of two class II genes of the baleen whale major histocompatibility complex were investigated and compared to toothed whales and out-groups. Amplification of the DQB exon 2 provided sequences showing high within-species and between-species nucleotide diversity and uninterrupted reading frames consistent with functional class II loci found in related mammals (e.g., ruminants). Cloning of amplified products indicated gene duplication in the humpback whale and triplication in the southern right whale, with average nucleotide diversity of 5.9 and 6.3%, respectively, for alleles of each species. Significantly higher nonsynonymous divergence at sites coding for peptide binding (32% for humpback and 40% for southern right) suggested that these loci were subject to positive (overdominant) selection. A population survey of humpback whales detected 23 alleles, differing by up to 21% of their inferred amino acid sequences. Amplification of the DRB exon 2 resulted in two groups of sequences. One was most similar to the DRB3 of the cow and present in all whales screened to date, including toothed whales. The second was most similar to the DRB2 of the cow and was found only in the bowhead and right whales. Both loci showed low diversity among species and apparent loss of function or altered function including interruption of reading frames. Finally, comparison of inferred protein sequence of the DRB3-like locus suggested convergence with the DQB, perhaps resulting from intergenic conversion or recombination.Electronic Supplementary Material Supplementary material is available for this article at