Climate-mediated adaptation after mainland colonization of an ancestrally subtropical island lizard, Anolis carolinensis.

Climate-mediated evolution plays an integral role in species migration and range expansion. Gaining a clearer understanding of how climate affects demographic history and adaptation provides fundamental insight into the generation of intra- and interspecific diversity. In this study, we used the natural colonization of the green anole (Anolis carolinensis) from the island of Cuba to mainland North America to investigate the role of evolution at the niche, phenotypic and genetic levels after long-term establishment in a novel environment. The North American green anole occupies a broader range of thermal habitats than its Cuban sister species. We documented niche expansion in the mainland green anole, mediated primarily through adaptation to winter temperatures. Common garden experiments strongly suggest a genetic component to differences in thermal performance found between populations in different temperature regimes. Analysis of geographic variation in population structure based on 53 486 single nucleotide variants from RAD loci revealed increased genetic isolation between populations in different vs. similar thermal environments. Selection scans for environment-allele correlations reveal 19 genomic loci of known function that may have played a role in the physiological adaptation of A. carolinensis to temperate environments on the mainland.

Sex-biased expression of sex-differentiating genes FOXL2 and FGF9 in American alligators, alligator Mississippiensis.

Across amniotes, sex-determining mechanisms exhibit great variation, yet the genes that govern sexual differentiation are largely conserved. Studies of evolution of sex-determining and sex-differentiating genes require an exhaustive characterization of functions of those genes such as FOXL2 and FGF9. FOXL2 is associated with ovarian development, and FGF9 is known to play a role in testicular organogenesis in mammals and other amniotes. As a step toward characterization of the evolutionary history of sexual development, we measured expression of FOXL2 and FGF9 across 3 developmental stages and 8 juvenile tissue types in male and female American alligators, Alligator mississippiensis. We report surprisingly high expression of FOXL2 before the stage of embryonic development when sex is determined in response to temperature, and sustained and variable expression of FGF9 in juvenile male, but not female tissue types. Novel characterization of gene expression in reptiles with temperature-dependent sex determination such as American alligators may inform the evolution of sex-determining and sex-differentiating gene networks, as they suggest alternative functions from which the genes may have been exapted. Future functional profiling of sex-differentiating genes should similarly follow other genes and other species to enable a broad comparison across sex-determining mechanisms.

Reptiles and mammals have differentially retained long conserved noncoding sequences from the amniote ancestor.

Many noncoding regions of genomes appear to be essential to genome function. Conservation of large numbers of noncoding sequences has been reported repeatedly among mammals but not thus far among birds and reptiles. By searching genomes of chicken (Gallus gallus), zebra finch (Taeniopygia guttata), and green anole (Anolis carolinensis), we quantified the conservation among birds and reptiles and across amniotes of long, conserved noncoding sequences (LCNS), which we define as sequences ≥500 bp in length and exhibiting ≥95% similarity between species. We found 4,294 LCNS shared between chicken and zebra finch and 574 LCNS shared by the two birds and Anolis. The percent of genomes comprised by LCNS in the two birds (0.0024%) is notably higher than the percent in mammals (<0.0003% to <0.001%), differences that we show may be explained in part by differences in genome-wide substitution rates. We reconstruct a large number of LCNS for the amniote ancestor (ca. 8,630) and hypothesize differential loss and substantial turnover of these sites in descendent lineages. By contrast, we estimated a small role for recruitment of LCNS via acquisition of novel functions over time. Across amniotes, LCNS are significantly enriched with transcription factor binding sites for many developmental genes, and 2.9% of LCNS shared between the two birds show evidence of expression in brain expressed sequence tag databases. These results show that the rate of retention of LCNS from the amniote ancestor differs between mammals and Reptilia (including birds) and that this may reflect differing roles and constraints in gene regulation.

Variability in sex-determining mechanisms influences genome complexity in reptilia.

In this review, we describe the history of amniote sex determination as a classic example of Darwinian evolution. We suggest that evolutionary changes in sex determination provide a foundation for understanding important aspects of chromosome and genome organization that otherwise appear haphazard in their origins and contents. Species with genotypic sex determination often possess heteromorphic sex chromosomes, whereas species with environmental sex determination lack them. Through a series of mutations followed by selection at key genes, sex-determining mechanisms have turned over many times throughout the amniote lineage. As a consequence, amniote genomes have undergone gains or losses of sex chromosomes. We review the genomic and ecological contexts in which either temperature-dependent or genotypic sex determination has evolved. Once genotypic sex determination emerges in a lineage, viviparity and heteromorphic sex chromosomes become more likely to evolve. For example, in extinct marine reptiles, genotypic sex determination apparently led to viviparity, which in turn facilitated their pelagic radiation. Sex chromosomes comprise genome regions that differ from autosomes in recombination rate, mutation rate, levels of polymorphism, and the presence of sex-determining and sexually antagonistic genes. In short, many aspects of amniote genome complexity, life history, and adaptive radiation appear contingent on evolutionary changes in sex-determining mechanisms.

Prevalence and evolutionary origins of autoimmune susceptibility alleles in natural mouse populations.

The evolutionary origin of genetic diversity in the SLAM/CD2 gene cluster, implicated in autoimmune lupus susceptibility in mice, was investigated by sequence analysis of exons from six members of the cluster in 48 wild mouse samples derived from the global mouse population. A total of 80 coding region SNPs were identified among the six genes analyzed, indicating that this gene cluster is highly polymorphic in natural mouse populations. Phylogenetic analyses of these allelic sequences revealed clustering of alleles derived from multiple Mus species and subspecies, indicating alleles at several SLAM/CD2 loci were present in ancestral Mus populations prior to speciation and have persisted as polymorphisms for more than 1 million years. Analyses of nonsynonymous/synonymous ratios using likelihood codon substitution models identified several segments in Cd229, Cd48 and Cd84 that were impacted by positive diversifying selective pressures. These findings support the interpretation that selection favoring the generation and retention of functional polymorphisms has played a role in the evolutionary origin of genetic polymorphisms that are predisposing to autoimmunity.

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 gamma-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.

Multiple origins of XY female mice (genus Akodon): phylogenetic and chromosomal evidence.

Despite the diversity in sex determination across organisms, theory predicts that the evolution of XY females is rare in mammals due to fitness consequences associated with infertility or the loss of YY zygotes. We investigated this hypothesis from a phylogenetic perspective by examining the inter- and intraspecific distribution of Y chromosomes in males and females (XY females) in South American field mice (Akodon). We found that XY females occurred at appreciable frequencies (10-66%) in at least eight Akodon species, raising the possibility that this system of sex determination has arisen multiple times independently. To determine the number of origins of XY females in Akodon, we constructed a molecular phylogeny of 16 species of Akodon based on mitochondrial DNA control region sequences. Both parsimony and maximum-likelihood reconstruction of ancestral states suggest that multiple steps (gains or losses of XY females) best explain the evolution of XY females, but do not clearly differentiate between single and multiple origins. We then directly compared functional and non-functional Y chromosomes in six species by Southern blot analysis. We found that male and female Y chromosome restriction fragment length polymorphism patterns were identical within species, but always differed between species, providing evidence that XY females arose at least six times within the Akodon lineage. To our knowledge, this pattern in Akodon is the first documentation of a novel sex-determining system arising multiple times within a tight clade of mammals. In addition, this system provides a clear test of the accuracy of phylogenetic methods to reconstruct ancestral states.

A 39-kb sequence around a blackbird Mhc class II gene: ghost of selection past and songbird genome architecture.

To gain an understanding of the evolution and genomic context of avian major histocompatibility complex (Mhc) genes, we sequenced a 38.8-kb Mhc-bearing cosmid insert from a red-winged blackbird (Agelaius phoeniceus). The DNA sequence, the longest yet retrieved from a bird other than a chicken, provides a detailed view of the process of gene duplication, divergence, and degeneration ("birth and death") in the avian Mhc, as well as a glimpse into major noncoding features of a songbird genome. The peptide-binding region (PBR) of the single Mhc class II B gene in this region, Agph-DAB2, is almost devoid of polymorphism, and a still-segregating single-base-pair deletion and other features suggest that it is nonfunctional. Agph-DAB2 is estimated to have diverged about 40 MYA from a previously characterized and highly polymorphic blackbird Mhc gene, Aph-DAB1, and is therefore younger than most mammalian Mhc paralogs and arose relatively late in avian evolution. Despite its nonfunctionality, Agph-DAB2 shows very high levels of nonsynonymous divergence from Agph-DAB1 and from reconstructed ancestral sequences in antigen-binding PBR codons-a strong indication of a period of adaptive divergence preceding loss of function. We also found that the region sequenced contains very few other unambiguous genes, a partial Mhc- class II gene fragment, and a paucity of simple-sequence and other repeats. Thus, this sequence exhibits some of the genomic streamlining expected for avian as compared with mammalian genomes, but is not as densely packed with functional genes as is the chicken Mhc.

Dynamics and phylogenetic implications of MtDNA control region sequences in New World Jays (Aves: Corvidae).

To study the evolution of mtDNA and the intergeneric relationships of New World Jays (Aves: Corvidae), we sequenced the entire mitochondrial DNA control region (CR) from 21 species representing all genera of New World jays, an Old World jay, crows, and a magpie. Using maximum likelihood methods, we found that both the transition/transversion ratio (kappa) and among site rate variation (alpha) were higher in flanking domains I and II than in the conserved central domain and that the frequency of indels was highest in domain II. Estimates of kappa and alpha were much more influenced by the density of taxon sampling than by alternative optimal tree topologies. We implemented a successive approximation method incorporating these parameters into phylogenetic analysis. In addition we compared our study in detail to a previous study using cytochrome b and morphology to examine the effect of taxon sampling, evolutionary rates of genes, and combined data on tree resolution. We found that the particular weighting scheme used had no effect on tree topology and little effect on tree robustness. Taxon sampling had a significant effect on tree robustness but little effect on the topology of the best tree. The CR data set differed nonsignificantly from the tree derived from the cytochrome b/morphological data set primarily in the placement of the genus Gymnorhinus, which is near the base of the CR tree. However, contrary to conventional taxonomy, the CR data set suggested that blue and black jays (Cyanocorax sensu lato) might be paraphyletic and that the brown jay Psilorhinus (=Cyanocorax) morio is the sister group to magpie jays (Calocitta), a phylogenetic hypothesis that is likely as parsimonious with regard to nonmolecular characters as monophyly of Cyanocorax. The CR tree also suggests that the common ancestor of NWJs was likely a cooperative breeder. Consistent with recent systematic theory, our data suggest that DNA sequences with high substitution rates such as the CR may nonetheless be useful in reconstructing relatively deep phylogenetic nodes in avian groups.

MHC class II pseudogene and genomic signature of a 32-kb cosmid in the house finch (Carpodacus mexicanus).

Large-scale sequencing studies in vertebrates have thus far focused primarily on the genomes of a few model organisms. Birds are of interest to genomics because of their much smaller and highly streamlined genomes compared to mammals. However, large-scale genetic work has been confined almost exclusively to the chicken; we know little about general aspects of genomes in nongame birds. This study examines the organization of a genomic region containing an Mhc class II B gene in a representative of another important lineage of the avian tree, the songbirds (Passeriformes). We used a shotgun sequencing approach to determine the sequence of a 32-kb cosmid insert containing a strongly hybridizing Mhc fragment from house finches (Carpodacus mexicanus). There were a total of three genes found on the cosmid clone, about the gene density expected for the mammalian Mhc: a class II Mhc beta-chain gene (Came-DAB1), a serine-threonine kinase, and a zinc finger motif. Frameshift mutations in both the second and third exons of Came-DAB1 and the unalignability of the gene after the third exon suggest that it is a nonfunctional pseudogene. In addition, the identifiable introns of Came-DAB1 are more than twice as large as those of chickens. Nucleotide diversity in the peptide-binding region of Came-DAB1 (Pi = 0.03) was much lower than polymorphic chicken and other functional Mhc genes but higher than the expected diversity for a neutral locus in birds, perhaps because of hitchhiking on a selected Mhc locus close by. The serine-threonine kinase gene is likely functional, whereas the zinc finger motif is likely nonfunctional. A paucity of long simple-sequence repeats and retroelements is consistent with emerging rules of chicken genomics, and a pictorial analysis of the "genomic signature" of this sequence, the first of its kind for birds, bears strong similarity to mammalian signatures, suggesting common higher-order structures in these homeothermic genomes. The house finch sequence is among a very few of its kind from nonmodel vertebrates and provides insight into the evolution of the avian Mhc and of avian genomes generally.

Perspective: gene divergence, population divergence, and the variance in coalescence time in phylogeographic studies.

Molecular methods as applied to the biogeography of single species (phylogeography) or multiple codistributed species (comparative phylogeography) have been productively and extensively used to elucidate common historical features in the diversification of the Earth's biota. However, only recently have methods for estimating population divergence times or their confidence limits while taking into account the critical effects of genetic polymorphism in ancestral species become available, and earlier methods for doing so are underutilized. We review models that address the crucial distinction between the gene divergence, the parameter that is typically recovered in molecular phylogeographic studies, and the population divergence, which is in most cases the parameter of interest and will almost always postdate the gene divergence. Assuming that population sizes of ancestral species are distributed similarly to those of extant species, we show that phylogeographic studies in vertebrates suggest that divergence of alleles in ancestral species can comprise from less than 10% to over 50% of the total divergence between sister species, suggesting that the problem of ancestral polymorphism in dating population divergence can be substantial. The variance in the number of substitutions (among loci for a given species or among species for a given gene) resulting from the stochastic nature of DNA change is generally smaller than the variance due to substitutions along allelic lines whose coalescence times vary due to genetic drift in the ancestral population. Whereas the former variance can be reduced by further DNA sequencing at a single locus, the latter cannot. Contrary to phylogeographic intuition, dating population divergence times when allelic lines have achieved reciprocal monophyly is in some ways more challenging than when allelic lines have not achieved monophyly, because in the former case critical data on ancestral population size provided by residual ancestral polymorphism is lost. In the former case differences in coalescence time between species pairs can in principle be explained entirely by differences in ancestral population size without resorting to explanations involving differences in divergence time. Furthermore, the confidence limits on population divergence times are severely underestimated when those for number of substitutions per site in the DNA sequences examined are used as a proxy. This uncertainty highlights the importance of multilocus data in estimating population divergence times; multilocus data can in principle distinguish differences in coalescence time (T) resulting from differences in population divergence time and differences in T due to differences in ancestral population sizes and will reduce the confidence limits on the estimates. We analyze the contribution of ancestral population size (theta) to T and the effect of uncertainty in theta on estimates of population divergence (tau) for single loci under reciprocal monophyly using a simple Bayesian extension of Takahata and Satta's and Yang's recent coalescent methods. The confidence limits on tau decrease when the range over which ancestral population size theta is assumed to be distributed decreases and when tau increases; they generally exclude zero when tau/(4Ne) > 1. We also apply a maximum-likelihood method to several single and multilocus data sets. With multilocus data, the criterion for excluding tau = 0 is roughly that l tau/(4Ne) > 1, where l is the number of loci. Our analyses corroborate recent suggestions that increasing the number of loci is critical to decreasing the uncertainty in estimates of population divergence time.

Polymorphism across an exon-intron boundary in an avian Mhc class II B gene.

Twenty-three sequence haplotypes spanning the boundary of the second exon and intron of a red-winged blackbird Mhc class II B gene, Agph-DAB1, are presented. The polymorphism of the exon segment is distributed in two divergent allelic lineages which appear to be maintained by balancing selection. The silent nucleotide diversity of the exon (pi = 0.101) is more than five times that of the intron (pi = 0.018) and decays rapidly across the exon-intron boundary. Additionally, genealogical reconstruction indicates that divergence from a common ancestor in the exon sample is over four times that of the intron. The intron sequences reveal a pattern of polymorphism which is characteristic of directional selection, rather than a pattern expected from linkage to a balanced polymorphism. These results suggest that the evolutionary histories of these two adjacent regions have been disassociated by recombination or gene conversion. The estimated population recombination parameter between the exon and the intron is sufficiently high (4NeC = 8.545) to explain the homogenization of intron sequences. Compatibility analyses estimate that these events primarily occur from the exon-intron boundary to about 20-30 bases into the intron. Additionally, the observation that divergent exon alleles share identical intron sequence supports the conclusion of disassociation of exon and intron evolutionary histories by recombination.

Toward an evolutionary genomics of the avian Mhc.

We review recent developments in the ongoing study of the evolution of the Mhc gene family in birds, with emphasis on class II B genes and results from songbirds obtained in our laboratory. Southern blots suggest a surprising diversity in Mhc class II gene number among various songbird species (Passeriformes). We have sequenced approximately 30 kb contigs from Mhc bearing cosmid clones from two species, red-winged blackbirds (Agelaius phoeniceus) and house finches (Carpodacus mexicanus), whose demography, lifetime reproductive success, epizootics, parasitology and mate choice are among the best studied for natural populations of birds. Of three genes cloned from these species, only one appears strongly polymorphic, and one (from the house finch) is likely a pseudogene. All are similar in structure to those in chickens, albeit with introns intermediate in length between chickens and mammals. Phylogenetic analysis of available class II B peptide-binding region exons suggests that the overwhelming long-term force operating on avian genes sampled thus far has been post-speciation gene duplication and/or concerted evolution. These and other results suggest that the evolution of class II B genes in birds conforms to a mixture of several models of multigene family evolution proposed for the mammalian Mhc, incorporating ongoing homogenization, duplication and pseudogene formation. Large-scale sequencing studies in these and other species, though still in their infancy, will prove invaluable for studying the comparative structures of avian Mhcs, as well as patterns of selection, mutation and linkage disequilibrium at several scales.

Genomics and polymorphism of Agph-DAB1, an Mhc class II B gene in red-winged blackbirds (Agelaius phoeniceus).

To further our understanding of the evolution of avian Mhc genes at the genomic level, we screened a cosmid library made from a red-winged blackbird (Agelaius phoeniceus) with a blackbird cDNA probe and subcloned from one of the Mhc-containing cosmids a gene which we designate Agph-DAB1. The structure of the gene is similar to that found for chicken class II B genes, except that the introns are surprisingly large, ranging from 98 to over 600 bp, making this the longest avian class II B gene to date. Using primers targeted toward the introns flanking the peptide-binding region (PBR), we amplified the entirety of the second exon and determined nucleotide sequences of 41 PCR products from eight individual blackbirds. The 10 sequence types found, among which were two probable pseudogene sequences, exhibit the classic hallmarks for evolution of PBRs, namely, an excess of nonsynonymous over synonymous substitutions and evidence of gene conversion events in polymorphic subdomains. Despite these patterns and our use of intron primers, the distribution of sequences among individuals suggests that more than one locus was amplified in most individuals, and the bushlike tree of sequences provides little information as to locus-specific clusters. These results imply a complex history of gene conversion, recent duplication, or possibly, concerted evolution among multiple loci, although Agph-DAB1, the first genomic Mhc sequence from a bird other than chicken, provides important clues in the quest for locus-specific Mhc primers in birds.

Can weighting improve bushy trees? Models of cytochrome b evolution and the molecular systematics of pipits and wagtails (Aves: Motacillidae).

Among-site rate variation (alpha) and transition bias (kappa) have been shown, most often as independent parameters, to be important dynamics in DNA evolution. Accounting for these dynamics should result in better estimates of phylogenetic relationships. To test this idea, we simultaneously estimated overall (averaged over all codon positions) and codon-specific values of alpha and kappa, using maximum likelihood analyses of cytochrome b data from all genera of pipits and wagtails (Aves: Motacillidae), and six outgroup species, using initial trees generated with default values. Estimates of alpha and kappa were robust to initial tree topology and suggested substantial among-site rate variation even within codon classes; alpha was lowest (large among-site rate variation) at second-codon and highest (low among-site rate variation) at third-codon positions. When overall values were applied, there were shifts in tree topology and dramatic and statistically significant improvements in log-likelihood scores of trees compared with the scores from application of default values. Applying codon-specific values resulted in yet another highly significant increase in likelihood. However, although incorporating substitution dynamics into maximum likelihood, maximum parsimony, and neighbor-joining analyses resulted in increases in congruence among trees, there were only minor improvements in phylogenetic signal, and none of the successive approximations tree topologies were statistically distinguishable from one another by the data. We suggest that the bushlike nature of many higher-level phylogenies in birds makes estimating the dynamics of DNA evolution less sensitive to tree topology but also less susceptible to improvement via weighting.

Evolution and ecology of MHC molecules: from genomics to sexual selection.

In the past few years the DNA sequence database for molecules of the MHC (major histocompatibility complex) has expanded greatly, yielding a more complete picture of the long-term rates and patterns of evolution of the MHC in vertebrates. Sharing of MHC allelic lineages between long-diverged species (trans-species evolution) has been detected virtually wherever it is sought, but new analyses of linked neutral regions and the complexities of sequence convergence and microrecombination in the peptide binding region challenge traditional phylogenetic analyses. Methods for estimating the intensity of selection on MHC genes suggest that viability is important, but recent studies in natural populations of mammals give inconsistent results concerning mate choice. The complex and interacting roles of microrecombination, parasite-mediated selection and mating preferences for maintaining the extraordinary levels of MHC polymorphism observed are still difficult to evaluate.

Ancestral polymorphism of Mhc class II genes in mice: implications for balancing selection and the mammalian molecular clock.

To investigate the evolutionary dynamics at Mhc class II DR genes of mice (genus Mus), we sequenced the peptide binding regions (PBRs) of 41 DRB (= E beta) genes and eight DRA (= E alpha) genes from 15 strains representing eight species. As expected trees of these PBR sequences imply extensive maintenance of ancestral DRB alleles across species. We use a coalescent simulation model to show that the number of interspecific coalescent events (c) observed on these trees was higher than the number expected for neutral genealogies and similar sample sizes and is more consistent with balancing selection that with neutrality. Patterns of ancestral polymorphism in mouse DRB alleles were also used to examine the tempo of synonymous substitution in the PBR of mouse class II genes. Both absolute and relative rate tests on DRA and DRB genes imply increased substitution rates at two- and fourfold degenerate sites of mice and rats relative to primates, and decreased rates for the DRB genes of primates relative to ungulate and carnivore relatives. Thus rates of synonymous substitution at Mhc DR genes in mammals appear to be subject to generation time effects in ways similar to those found at other mammalian genes.

Re: Congruence and phylogenetic reanalysis of perching bird cytochrome b sequences.

In a study of the phylogenetic relationships among perching bird mtDNA sequences, Edwards et al. (1991; Proc. R. Soc. London Ser B. 243: 99-107) sequenced part of the mitochondrial cytochrome b gene from 13 perching birds (Passeriformes) and a woodpecker outgroup. However, recently the validity of part of the sequence of the hermit thrush (Catharus guttatus) in that study has been questioned. To determine the effect of inclusion of this apparently chimeric sequence on the conclusions of the original analysis, we reanalyzed these sequences using a different published hermit thrush sequence. In addition, we applied tests of congruence to examine the possibility that the aberrant phylogenetic behavior of chimeric mtDNA sequences might be detected. The reanalysis confirms the ability of slow evolving first and second codon positions of cytochrome b sequences to resolve deep branches in the avian tree. The fact that the new data set does not reject the DNA hybridization tree of these species probably indicates poor ability of the cytochrome b sequences to discriminate among alternative trees, rather than consistency among data sets. Statistical testing of trees based on individual amplified segments of mtDNA indicates that congruence tests may be one useful way of identifying chimeric mtDNA sequences when they have not been detected in the laboratory or during standard phylogenetic analysis.

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 beta 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 beta 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.