Polymorphism and gene organization of water buffalo MHC-DQB genes show homology to the BoLA DQB region.

In cattle (Bos taurus), there is evidence of more than 50 alleles of BoLA-DQB (bovine lymphocyte antigen DQB) that are distributed across at least five DQB loci, making this region one of the most complex in the BoLA gene family. In this study, DQB alleles were analysed for the water buffalo (Bubalus bubalis), another economically important bovine species. Twelve alleles for Bubu-DQB (Bubalis bubalis DQB) were determined by nucleotide sequence analysis. A phylogenetic analysis revealed numerous trans-species polymorphisms, with alleles from water buffalo assigned to at least three different loci (BoLA-DQB1, BoLA-DQB3 and BoLA-DQB4) that are also found in cattle. These presumptive loci were analysed for patterns of synonymous (d(S)) and non-synonymous (d(N)) substitution. Like BoLA-DQB1, Bubu-DQB1 was observed to be under strong positive selection for polymorphism. We conclude that water buffalo and cattle share the current arrangement of their DQB region because of their common ancestry.

Assessing substitution patterns, rates and homoplasy at HVRI of Steller sea lions, Eumetopias jubatus.

Despite the widely recognized incidence of homoplasy characterizing this region, the hypervariable region I (HVRI) of the mitochondrial control region is one of the most frequently used genetic markers for population genetic and phylogeographic studies. We present an evolutionary analysis of HVRI and cytochrome b sequences from a range-wide survey of 1031 Steller sea lions, Eumetopias jubatus, to quantify homoplasy and substitution rate at HVRI. Variation in HVRI was distributed across 41 variable sites in the 238-bp segment examined. All variants at HVR1 were found to be transitions. However, our analyses suggest that a minimum of 101 changes have actually occurred within HVRI with as many as 18 substitutions occurring at a single site. By including this hidden variation into our analyses, several instances of apparent long-range dispersal were resolved to be homoplasies and 8.5-12% of observed HVRI haplotypes were found to have geographic distributions descriptive of convergent molecular evolution rather than identity by descent. We estimate the rate of substitution at HVRI in Steller sea lions to be approximately 24 times that of cytochrome b with an absolute rate of HVRI substitution estimated at 27.45% per million years. These findings have direct implications regarding the utility of HVRI data to generate a variety of evolutionary genetic hypotheses.

Glacial refugia and the phylogeography of Steller's sea lion (Eumatopias jubatus) in the North Pacific.

Mitochondrial DNA sequence data were used to examine the phylogeographic history of Steller's sea lions (Eumetopias jubatus) in relation to the presence of Plio-Pleistocene insular refugia. Cytochrome b and control region sequences from 336 Steller's sea lions reveal phylogenetic lineages associated with continental refugia south of the ice sheets in North America and Eurasia. Phylogenetic analysis suggests the genetic structure of E. jubatus is the result of Pleistocene glacial geology, which caused the elimination and subsequent reappearance of suitable rookery habitat during glacial and interglacial periods. The cyclic nature of geological change produced a series of independent population expansions, contractions and isolations that had analogous results on Steller's sea lions and other marine and terrestrial species. Our data show evidence of four glacial refugia in which populations of Steller's sea lions diverged. These events occurred from approximately 60,000 to 180,000 years BP and thus preceded the last glacial maximum.

Gradual and quantum genome size shifts in the hystricognath rodents.

We assessed genome size variation by flow cytometry within and among 31 species of nine families of African and South American hystricognath rodents. Interspecific variation was extensive and genome size was relatively high among the South American radiation whereas only moderate variation and smaller estimates of genome size were observed in the African counterparts. The largest genome size, indicating tetraploidy was recorded in the South American octodontid, Tympanoctomys barrerae (16.8 pg DNA). This quantum shift in DNA content represents a novel mechanism of genome evolution in mammals. As expected in polyploid organisms, varying nucleotypic effects were observed in the dimensions of the sperm cells and lymphocytes of T. barrerae. The role of control mechanisms that influence cell dimensions in polyploid organisms is discussed.

Polymorphisms in MHC-DRA and -DRB alleles of water buffalo (Bubalus bubalis) reveal different features from cattle DR alleles.

Seventy-five individuals of Bubalus bubalis belonging to four different breeds, three of river buffalo and one of swamp buffalo, were studied for polymorphism in MHC DRB (Bubu-DRB) and DRA (Bubu-DRA) loci. Eight alleles of Bubu-DRB were found, and all alleles in the swamp type were shared with the three river breeds. All alleles sampled from the breed of European origin (Mediterranean) were present in breeds sampled in Brazil, thus variability of this locus may have been preserved to a great extent in the more recently founded Brazilian population. Bubu-DRB alleles contained higher proportions of synonymous vs. non-synonymous substitutions in the non-peptide-binding sites (PBS) region, in contrast to the pattern of variation found in BoLA-DRB3, the orthologous locus in cattle. This indicated that either the first domain exon (exon 2) of Bubu-DRB has not undergone as much recombination and/or gene conversion as in cattle alleles, or Bubu-DRB may be more ancient than BoLA-DRB3 alleles. Phylogenetic analysis of DRB alleles from Bubalus, Syncerus c. caffer, the Cape buffalo, and domestic cattle demonstrated transspecies polymorphism. Water buffalo contained two alleles of DRA that differed from each other in two amino acid positions, including one in the PBS (alpha22) that was also shared with Anoa depressicornis, the anoa. Discovery of variation in DRA was surprising as the first domain of DRA is a highly conserved polypeptide in mammals in general and especially in ruminants, where no other substitution in PBS was seen.

Historical population size change of bowhead whales inferred from DNA sequence polymorphism data.

Nucleotide sequence data from the mitochondrial control region were used from a phylogenetic context to investigate the long-term history of a population of bowhead whales (Balaena mysticetus). In addition, the coalescence time of these sequences was used to estimate the age of the inferred patterns of population size change. The results indicate that mitochondrial genetic polymorphism was not affected by a recent bottleneck that occurred near the turn of the 20th century, thereby preserving the signature of historical population size change in the mitochondrial genome. Further analysis showed that this population underwent an expansion initiated in the Middle to Late Pleistocene. As such, early Holocene changes in Arctic sea ice distribution appear to have had little influence on patterns of genetic variability in this population.

Molecular phylogeny and divergence time estimates for major rodent groups: evidence from multiple genes.

The order Rodentia contains half of all extant mammal species, and from an evolutionary standpoint, there are persistent controversies surrounding the monophyly of the order, divergence dates for major lineages, and relationships among families. Exons of growth hormone receptor (GHR) and breast cancer susceptibility (BRCA1) genes were sequenced for a wide diversity of rodents and other mammals and combined with sequences of the mitochondrial 12S rRNA gene and previously published sequences of von Willebrand factor (vWF). Rodents exhibit rates of amino acid replacement twice those observed for nonrodents, and this rapid rate of evolution influences estimates of divergence dates. Based on GHR sequences, monophyly is supported, with the estimated divergence between hystricognaths and most sciurognaths dating to about 75 MYA. Most estimated dates of divergence are consistent with the fossil record, including a date of 23 MYA for Mus-Rattus divergence. These dates are considerably later than those derived from some other molecular studies. Among combined and separate analyses of the various gene sequences, moderate to strong support was found for several clades. GHR appears to have greater resolving power than do 12S or vWF. Despite its complete unresponsiveness to growth hormone, Cavia (and other hystricognaths) exhibits a conservative rate of change in the intracellular domain of GHR.

Whence the red panda?

The evolutionary history of the red panda (Ailurus fulgens) plays a pivotal role in the higher-level phylogeny of the "bear-like" arctoid carnivoran mammals. Characters from morphology and molecules have provided inconsistent evidence for placement of the red panda. Whereas it certainly is an arctoid, there has been major controversy about whether it should be placed with the bears (ursids), ursids plus pinnipeds (seals, sea lions, walrus), raccoons (procyonids), musteloids (raccoons plus weasels, skunks, otters, and badgers [mustelids]), or as a monotypic lineage of uncertain phylogenetic affinities. Nucleotide sequence data from three mitochondrial genes and one nuclear intron were analyzed, with more complete taxonomic sampling of relevant taxa (arctoids) than previously available in analyses of primary molecular data, to clarify the phylogenetic relationships of the red panda to other arctoid carnivorans. This study provides detailed phylogenetic analyses (both parsimony and maximum-likelihood) of primary character data for arctoid carnivorans, including bootstrap and decay indices for all arctoid nodes, and three statistical tests of alternative phylogenetic hypotheses for the placement of the red panda. Combined phylogenetic analyses reject the hypotheses that the red panda is most closely related to the bears (ursids) or to the raccoons (procyonids). Rather, evidence from nucleotide sequences strongly support placement of the red panda within a broad Musteloidea (sensu lato) clade, including three major lineages (the red panda, the skunks [mephitids], and a clearly monophyletic clade of procyonids plus mustelids [sensu stricto, excluding skunks]). Within the Musteloidea, interrelationships of the three major lineages are unclear and probably are best considered an unresolved trichotomy. These data provide compelling evidence for the relationships of the red panda and demonstrate that small taxonomic sample sizes can result in misleading or possibly erroneous (based on prior modeling, as well as conflict between the results of our analyses of less and more complete data sets) conclusions about phylogenetic relationships and taxonomy.

Evidence from intron 1 of the nuclear transthyretin (Prealbumin) gene for the phylogeny of African mole-rats (Bathyergidae).

A 900- to 1100-bp fragment encompassing intron 1 of the nuclear transthyretin (prealbumin) gene was examined in 12 taxa of Old World hystricognath rodents of the families Bathyergidae, Petromuridae, Thryonomyidae, and Hystricidae. Within the Bathyergidae, Heterocephalus glaber (naked mole-rat) was basal, and the other East African species, Heliophobius argenteocinereus (silvery mole-rat), was sister to a southern African clade containing Bathyergus, Cryptomys, and Georychus (dune, common, and cape mole-rats). These results are congruent with studies using mitochondrial 12S rRNA gene sequences. A combined analysis of transthyretin and 12S rRNA data resulted in a well-supported topology with better resolution than either gene analyzed separately. These data support the findings by M. W. Allard and R. L. Honeycutt (1992, Mol. Biol. Evol. 9: 27-40) and R. L. Honeycutt (1992, Am. Sci. 80: 43-53) that complex social systems evolved independently at least twice, in the common and naked mole-rats.

Evaluating a putative bottleneck in a population of bowhead whales from patterns of microsatellite diversity and genetic disequilibria.

A size-selected Balaena mysticetus genomic library was screened for clones containing simple sequence repeat, or microsatellite, loci. A total of 11 novel loci was identified. These loci were combined with a set of 9 published loci, for a total of 20 markers, and were scored across a sample of 108 bowhead whales from the Bering-Chukchi-Beaufort Seas population of bowhead whales. Genetic variability was measured in terms of polymorphism information content values and unbiased heterozygosity. From the latter, estimates of long-term effective population size were obtained. In addition, gametic phase disequilibrium among loci was investigated. Moderate to high levels of polymorphism were found overall, and the long-term effective size estimates were large relative to total population size. Tests of heterozygosity excess (Cornuet and Luikart 1996) and allele frequency distribution (Luikart et al. 1998) indicated that the possibility of a recent genetic bottleneck in the Bering-Chukchi-Beaufort Seas population of bowhead whales is highly unlikely. However, the fact that five loci displayed a statistically significant heterozygote deficiency remains to be explained.

Genetic evidence of Dobrava virus in Apodemus agrarius in Hungary.

Using nested polymerase chain reaction, we sequenced Dobrava virus (DOB) from the rodent Apodemus agrarius in Hungary. The samples we isolated group with DOB samples previously isolated from A. flavicollis. This grouping may indicate host switching.

Nucleotide polymorphism at the alcohol dehydrogenase locus of pocket gophers, genus Geomys.

Using the strictly neutral model as a null hypothesis, we tested for deviations from expected levels of nucleotide polymorphism at the alcohol dehydrogenase locus (Adh-1) within and among four species of pocket gophers (Geomys bursarius major, G. knoxjonesi, G. texensis llanensis, and G. attwateri). The complete protein-encoding region was examined, and 10 unique alleles, representing both electromorphic and cryptic alleles, were used to test hypotheses (e.g., the neutral model) concerning the maintenance of genetic variation. Nineteen variable sites were identified among the 10 alleles examined, including 9 segregating sites occurring in synonymous positions and 10 that were nonsynonymous. Several statistical methods, including those that test for within-species variation as well as those that examine variation within and among species, failed to reject the null hypothesis that variation (both within and between species of Geomys) at the Adh locus is consistent with the neutral theory. However, there was significant heterogeneity in the ratio of polymorphism to divergence across the gene, with polymorphisms clustered in the first half of the coding region and fixed differences clustered in the second half of the gene. Two alternative hypotheses are discussed as possible explanations for this heterogeneity: an old balanced polymorphism in the first half of the gene or a recent selective sweep in the second half of the gene.

Nucleotide sequence evolution at the kappa-casein locus: evidence for positive selection within the family Bovidae.

Kappa-casein is a mammalian milk protein involved in a number of important physiological processes. In the gut, the ingested protein is split into an insoluble peptide (para kappa-casein) and a soluble hydrophilic glycopeptide (caseinomacropeptide). Caseinomacropeptide is responsible for increased efficiency of digestion, prevention of neonate hypersensitivity to ingested proteins, and inhibition of gastric pathogens. Variation within this peptide has significant effects associated with important traits such as milk production. The nucleotide sequences for regions of kappa-casein exon and intron four were determined for representatives of the artiodactyl family Bovidae. The pattern of nucleotide substitution in kappa-casein sequences for distantly related bovid taxa demonstrates that positive selection has accelerated their divergence at the amino acid sequence level. This selection has differentially influenced the molecular evolution of the two kappa-casein split peptides and is focused within a 34-codon region of caseinomacropeptide.

Multiple and ancient origins of the domestic dog.

Mitochondrial DNA control region sequences were analyzed from 162 wolves at 27 localities worldwide and from 140 domestic dogs representing 67 breeds. Sequences from both dogs and wolves showed considerable diversity and supported the hypothesis that wolves were the ancestors of dogs. Most dog sequences belonged to a divergent monophyletic clade sharing no sequences with wolves. The sequence divergence within this clade suggested that dogs originated more than 100,000 years before the present. Associations of dog haplotypes with other wolf lineages indicated episodes of admixture between wolves and dogs. Repeated genetic exchange between dog and wolf populations may have been an important source of variation for artificial selection.

Evolution of eutherian cytochrome c oxidase subunit II: heterogeneous rates of protein evolution and altered interaction with cytochrome c.

Cytochrome c oxidase subunit II (COII), encoded by the mitochondrial genome, exhibits one of the most heterogeneous rates of amino acid replacement among placental mammals. Moreover, it has been demonstrated that cytochrome c oxidase has undergone a structural change in higher primates which has altered its physical interaction with cytochrome c. We collected a large data set of COII sequences from several orders of mammals with emphasis on primates, rodents, and artiodactyls. Using phylogenetic hypotheses based on data independent of the COII gene, we demonstrated that an increased number of amino acid replacements are concentrated among higher primates. Incorporating approximate divergence dates derived from the fossil record, we find that most of the change occurred independently along the New World monkey lineage and in a rapid burst before apes and Old World monkeys diverged. There is some evidence that Old World monkeys have undergone a faster rate of nonsynonymous substitution than have apes. Rates of substitution at four-fold degenerate sites in primates are relatively homogeneous, indicating that the rate heterogeneity is restricted to nondegenerate sites. Excluding the rate acceleration mentioned above, primates, rodents, and artiodactyls have remarkably similar nonsynonymous replacement rates. A different pattern is observed for transversions at four-fold degenerate sites, for which rodents exhibit a higher rate of replacement than do primates and artiodactyls. Finally, we hypothesize specific amino acid replacements which may account for much of the structural difference in cytochrome c oxidase between higher primates and other mammals.

Mitochondrial gene sequences and the molecular systematics of the artiodactyl subfamily bovinae.

Nucleotide sequence evolution of the mitochondrial cytochrome c oxidase subunit II (COII) gene was used to examine the molecular phylogenetics and evolution of the Bovinae, a subfamily within the mammalian order Artiodactyla. The COII gene was sequenced in representatives of three bovine tribes (Bovini, Boselaphini, and Tragelaphini) and the outgroup taxon Capra (subfamily Caprinae). Although the phylogenetic analyses grouped Bison as sister to Bos, the genus Bison was paraphyletic, with the American bison being most closely related to species of Bos rather than to the European bison. COII data also supported a close relationship between African (Syncerus) and Asian (bubalus) buffaloes, the monophyly of the tribe Bovini, and a sister-group relationship between the tribes Bovini and Boselaphini. Analysis of nucleotide substitutions in the COII gene prompted a system of differential weighting of nucleotide substitutions for inferring phylogenetic relationships across the range of divergence times examined here (2-20 million years). Rates of evolution in the COII gene are examined and compared to evolutionary rates in mtDNA tRNA/rRNA genes and the D-loop among other artiodactyl taxa.