INTROGRESSIVE HYBRIDIZATION AND NONCONCORDANT EVOLUTIONARY HISTORY OF MATERNAL AND PATERNAL LINEAGES IN NORTH AMERICAN DEER

Antler tine homologies and cervid systematics: A review of past and present controversies with special emphasis on Elaphurus davidianus

Antlers are the most conspicuous trait of cervids and have been used in the past to establish a classification of their fossil and living representatives. Since the availability of molecular data, morphological characters have generally become less important for phylogenetic reconstructions. In recent years, however, the appreciation of morphological characters has increased, and they are now more frequently used in addition to molecular data to reconstruct the evolutionary history of cervids. A persistent challenge when using antler traits in deer systematics is finding a consensus on the homology of structures. Here, we review early and recent attempts to homologize antler structures and objections to this approach, compare and evaluate recent advances on antler homologies, and critically discuss these different views in order to offer a basis for further scientific exchange on the topic. We further present some developmental aspects of antler branching patterns and discuss their potential for reconstructing cervid systematics. The use of heterogeneous data for reconstructing phylogenies has resulted in partly conflicting hypotheses on the systematic position of certain cervid species, on which we also elaborate here. We address current discussions on the use of different molecular markers in cervid systematics and the question whether antler morphology and molecular data can provide a consistent picture on the evolutionary history of cervids. In this context, special attention is given to the antler morphology and the systematic position of the enigmatic Pere David's deer (Elaphurus davidianus).

Mitochondrial capture and subsequent genetic divergence generates a novel haplogroup: evidence from ancient and ongoing hybridization in mule and white-tailed deer

Odocoileus virginianus (white-tailed deer) and O. hemionus (mule deer) are sympatric across much of North America. Molecular evidence suggests that up to 24% of individuals in some populations are a product of hybrid ancestry. Several studies have alluded to ancient and recent introgression between Odocoileus spp.; however, no divergence dates were proposed. Herein, phylogenetic analyses of DNA sequences obtained from the mitochondrial Cytochrome b gene in 690 individuals identified three clades corresponding to black-tailed deer, white-tailed deer, or a unique combination of both white-tailed deer and mule deer. White-tailed deer and mule deer diverged from a common ancestor of approximately 3.13 mya followed by an ancient hybridization event of approximately 1.32 mya, in which the white-tailed deer mitochondrial genome was “captured” by mule deer. This hybridization event produced a novel haplogroup for white-tailed deer and mule deer located west of the Appalachian Mountains and east of the Cascade Range, south to Veracruz, Mexico, and north to the Yukon Territory, Canada. The ancestral mule deer-like mitochondrial genome appears to be restricted to black-tailed deer distributed along the western portion of the Cascade and Sierra Nevada Ranges of the United States and Canada, whereas the ancestral white-tailed deer-like mitochondrial genome is restricted to the eastern United States and portions of Latin America and Caribbean regions. The “captured mitochondrial genome” has continued on an independent evolutionary trajectory and represents a unique and broadly distributed haplogroup that is 7.25% and 2.84% different from the ancestral mule deer and ancestral white-tailed deer haplogroups, respectively.

Population genomics of free-ranging Great Plains white-tailed and mule deer reflects a long history of interspecific hybridization

Hybridization is a natural process at species-range boundaries that may variably promote the speciation process or break down species barriers but minimally will influence management outcomes of distinct populations. White-tailed deer (Odocoileus virginianus) and mule deer (Odocoileus hemionus) have broad and overlapping distributions in North America and a recognized capacity for interspecific hybridization. In response to contemporary environmental change to any of one or multiple still-unknown factors, mule deer range is contracting westward accompanied by a westward expansion of white-tailed deer, leading to increasing interactions, opportunities for gene flow, and associated conservation implications. To quantify genetic diversity, phylogenomic structure, and dynamics of hybridization in sympatric populations of white-tailed and mule deer, we used mitochondrial cytochrome b data coupled with SNP loci discovered with double-digest restriction site-associated DNA sequencing. We recovered 25,018 SNPs across 92 deer samples from both species, collected from two regions of western Kansas. Eight individuals with unambiguous external morphology representing both species were of hybrid origin (8.7%), and represented the product of multi-generational backcrossing. Mitochondrial data showed both ancient and recent directional discordance with morphological species assignments, reflecting a legacy of mule deer males mating with white-tailed deer females. Mule deer had lower genetic diversity than white-tailed deer, and both mitochondrial and nuclear data suggest contemporary mule deer effective population decline. Landscape genetic analyses show relative isolation between the two study regions for white-tailed deer, but greater connectivity among mule deer, with predominant movement from north to south. Collectively, our results suggest a long history of gene flow between these species in the Great Plains and hint at evolutionary processes that purge incompatible functional genomic elements as a result of hybridization. Surviving hybrids evidently may be reproductive, but with unknown consequences for the future integrity of these species, population trajectories, or relative susceptibility to emerging pathogens.

Technical Limitations Associated With Molecular Barcoding of Arthropod Bloodmeals Taken From North American Deer Species

Accurate species-level identification of the source of arthropod bloodmeals is important for deciphering blood feeding patterns of field-collected specimens. Cytochrome c oxidase I (COI) mitochondrial gene sequencing has been used for this purpose; however, species resolution can be difficult to obtain from certain vertebrate genera, including Odocoileus. Sanger sequencing of mitochondrial genes was employed to identify the bloodmeal source of wild-caught mosquitoes trapped in Greeley, Colorado. Initial sequencing of the COI gene of mitochondrial DNA in bloodmeals was inadequate for species-level resolution of bloodmeals from deer in the genus Odocoileus, with current databases returning low fidelity matches to multiple genera. The use of the hypervariable D loop of the control region provided species-level identification of white-tailed deer (Order: Artiodactyla, Family: Cervidae, Odocoileus virginianus); however, taxonomic identification was successful only to genus for mule (O. hemionus hemionus) and black-tailed deer (O. hemionus columbianus). We advocate the use of multiple loci for bloodmeal analysis and the buildout of available databases to include multiple mitochondrial reference genes for reliable host species identification.

Sperm quality, aggressiveness and generation turnover may facilitate unidirectional Y chromosome introgression across the European house mouse hybrid zone

The widespread and locally massive introgression of Y chromosomes of the eastern house mouse (Mus musculus musculus) into the range of the western subspecies (M. m. domesticus) in Central Europe calls for an explanation of its underlying mechanisms. Given the paternal inheritance pattern, obvious candidates for traits mediating the introgression are characters associated with sperm quantity and quality. We can also expect traits such as size, aggression or the length of generation cycles to facilitate the spread. We have created two consomic strains carrying the non-recombining region of the Y chromosome of the opposite subspecies, allowing us to study introgression in both directions, something impossible in nature due to the unidirectionality of introgression. We analyzed several traits potentially related to male fitness. Transmission of the domesticus Y onto the musculus background had negative effects on all studied traits. Likewise, domesticus males possessing the musculus Y had, on average, smaller body and testes and lower sperm count than the parental strain. However, the same consomic males tended to produce less- dissociated sperm heads, to win more dyadic encounters, and to have shorter generation cycles than pure domesticus males. These data suggest that the domesticus Y is disadvantageous on the musculus background, while introgression in the opposite direction can confer a recognizable, though not always significant, selective advantage. Our results are thus congruent with the unidirectional musculus → domesticus Y chromosome introgression in Central Europe. In addition to some previous studies, they show this to be a multifaceted phenomenon demanding a multidisciplinary approach.

Evidence for adaptive introgression of exons across a hybrid swarm in deer

Background

Secondary contact between closely related lineages can result in a variety of outcomes, including hybridization, depending upon the strength of reproductive barriers. By examining the extent to which different parts of the genome introgress, it is possible to infer the strength of selection and gain insight into the evolutionary trajectory of lineages. Following secondary contact approximately 8000 years ago in the Pacific Northwest, mule deer (Odocoileus hemionus hemionus) and black-tailed deer (O. h. columbianus) formed a hybrid swarm along the Cascade mountain range despite substantial differences in body size (up to two times) and habitat preference. In this study, we examined genetic population structure, extent of introgression, and selection pressures in freely interbreeding populations of mule deer and black-tailed deer using mitochondrial DNA sequences, 9 microsatellite loci, and 95 SNPs from protein-coding genes.

Results

We observed bi-directional hybridization and classified approximately one third of the 172 individuals as hybrids, almost all of which were beyond the F1 generation. High genetic differentiation between black-tailed deer and mule deer at protein-coding genes suggests that there is positive divergent selection, though selection on these loci is relatively weak. Contrary to predictions, there was not greater selection on protein-coding genes thought to be associated with immune function and mate choice. Geographic cline analyses were consistent across genetic markers, suggesting long-term stability (over hundreds of generations), and indicated that the center of the hybrid swarm is 20-30 km to the east of the Cascades ridgeline, where there is a steep ecological transition from wet, forested habitat to dry, scrub habitat.

Conclusions

Our data are consistent with a genetic boundary between mule deer and black-tailed deer that is porous but maintained by many loci under weak selection having a substantial cumulative effect. The absence of clear reproductive barriers and the consistent centering of geographic clines at a sharp ecotone suggests that ecology is a driver of hybrid swarm dynamics. Adaptive introgression in this study (and others) promotes gene flow and provides valuable insight into selection strength on specific genes and the evolutionary trajectory of hybridizing taxa.

Electronic supplementary material

The online version of this article (10.1186/s12862-019-1497-x) contains supplementary material, which is available to authorized users.

A continent-scale study of the social structure and phylogeography of the bent-wing bat, Miniopterus schreibersii (Mammalia: Chiroptera), using new microsatellite data

Miniopterus schreibersii is a cave-dwelling bat species with a wide distribution in the western Palearctic spanning southern and central Europe, North Africa, Anatolia, the Middle East, and the Caucasus. We investigated the social structure and its effects on the genetic makeup of this species, using 10 nuclear microsatellite markers and a partial fragment of the mitochondrial cytochrome b gene. Samples were examined from the species' entire circum-Mediterranean range. Local structuring that was previously detected among populations of M. schreibersii using mitochondrial markers was not observed for microsatellite markers, indicating male-biased dispersal for the species. Some support was found for postglacial expansions in Europe, with Anatolia potentially acting as the primary refugium during the Last Glacial Maximum (LGM). However, support for this hypothesis is not as strong as that previously detected using mitochondrial DNA markers. This is likely due to the diminishing effect of male-mediated dispersal, replenishing the nuclear diversity faster than the mitochondrial diversity in regions that are relatively far from the glacial refugia.

Widespread hybridization in the introduced hog deer population of Victoria, Australia, and its implications for conservation

In Australia, many species have been introduced that have since undergone drastic declines in their native range. One species of note is the hog deer (Axis porcinus) which was introduced in the 1860s to Victoria, Australia, and has since become endangered in its native range throughout South-East Asia. There is increased interest in using non-native populations as a source for genetic rescue; however, considerations need to be made of the genetic suitability of the non-native population. Three mitochondrial markers and two nuclear markers were sequenced to assess the genetic variation of the Victorian population of hog deer, which identified that the Victorian population has hybrid origins with the closely related chital (Axis axis), a species that is no longer present in the wild in Victoria. In addition, the mitochondrial D-loop region within the Victorian hog deer is monomorphic, demonstrating that mitochondrial genetic diversity is very low within this population. This study is the first to report of long-term persistence of hog deer and chital hybrids in a wild setting, and the continual survival of this population suggests that hybrids of these two species are fertile. Despite the newly discovered hybrid status in Victorian hog deer, this population may still be beneficial for future translocations within the native range. However, more in-depth analysis of genetic diversity within the Victorian hog deer population and investigation of hybridization rates within the native range are necessary before translocations are attempted.

Development of a Novel Mule Deer Genomic Assembly and Species-Diagnostic SNP Panel for Assessing Introgression in Mule Deer, White-Tailed Deer, and Their Interspecific Hybrids

Mule deer (Odocoileus hemionus) are endemic to a wide variety of habitats in western North America, many of which are shared in sympatry with their closely related sister-species white-tailed deer (Odocoileus virginianus), whom they hybridize with in wild populations. Although mule deer meet many ideal conditions for a molecular ecological research species, such as high abundance, ecological importance, and broad dispersal and gene flow, conservation genetic studies have been limited by a relative lack of existing genomic resources and inherent difficulties caused by introgression with white-tailed deer. Many molecular tools currently available for the study of cervids were designed using reference assemblies of divergent model species, specifically cattle (Bos taurus). Bovidae and Cervidae diverged approximately 28 million years ago, therefore, we sought to ameliorate the available resources by contributing the first mule deer whole genome sequence draft assembly with an average genome-wide read depth of 25X, using the white-tailed genome assembly (Ovir.te_1.0) as a reference. Comparing the two assemblies, we identified ∼33 million single nucleotide polymorphisms (SNPs) and insertion/deletion variants. We then verified fixed SNP differences between the two species and developed a 40-loci SNP assay capable of identifying pure mule deer, white-tailed deer, and interspecific hybrids. Assignment capacity of the panel, which was tested on simulated datasets, is reliable up to and including the third backcross hybrid generation. Identification of post-F1 hybrids will be necessary for hybrid zone population studies going forward, and the new mule deer assembly will be a valuable resource for genetic and comparative genomics studies.

Size-assortative choice and mate availability influences hybridization between red wolves (Canis rufus) and coyotes (Canis latrans)

Anthropogenic hybridization of historically isolated taxa has become a primary conservation challenge for many imperiled species. Indeed, hybridization between red wolves (Canis rufus) and coyotes (Canis latrans) poses a significant challenge to red wolf recovery. We considered seven hypotheses to assess factors influencing hybridization between red wolves and coyotes via pair‐bonding between the two species. Because long‐term monogamy and defense of all‐purpose territories are core characteristics of both species, mate choice has long‐term consequences. Therefore, red wolves may choose similar‐sized mates to acquire partners that behave similarly to themselves in the use of space and diet. We observed multiple factors influencing breeding pair formation by red wolves and found that most wolves paired with similar‐sized conspecifics and wolves that formed congeneric pairs with nonwolves (coyotes and hybrids) were mostly female wolves, the smaller of the two sexes. Additionally, we observed that lower red wolf abundance relative to nonwolves and the absence of helpers increased the probability that wolves consorted with nonwolves. However, successful pairings between red wolves and nonwolves were associated with wolves that maintained small home ranges. Behaviors associated with territoriality are energetically demanding and behaviors (e.g., aggressive interactions, foraging, and space use) involved in maintaining territories are influenced by body size. Consequently, we propose the hypothesis that size disparities between consorting red wolves and coyotes influence positive assortative mating and may represent a reproductive barrier between the two species. We offer that it may be possible to maintain wild populations of red wolves in the presence of coyotes if management strategies increase red wolf abundance on the landscape by mitigating key threats, such as human‐caused mortality and hybridization with coyotes. Increasing red wolf abundance would likely restore selection pressures that increase mean body and home‐range sizes of red wolves and decrease hybridization rates via reduced occurrence of congeneric pairs.

The genomic ancestry, landscape genetics and invasion history of introduced mice in New Zealand

The house mouse (Mus musculus) provides a fascinating system for studying both the genomic basis of reproductive isolation, and the patterns of human-mediated dispersal. New Zealand has a complex history of mouse invasions, and the living descendants of these invaders have genetic ancestry from all three subspecies, although most are primarily descended from M. m. domesticus. We used the GigaMUGA genotyping array (approximately 135 000 loci) to describe the genomic ancestry of 161 mice, sampled from 34 locations from across New Zealand (and one Australian city-Sydney). Of these, two populations, one in the south of the South Island, and one on Chatham Island, showed complete mitochondrial lineage capture, featuring two different lineages of M. m. castaneus mitochondrial DNA but with only M. m. domesticus nuclear ancestry detectable. Mice in the northern and southern parts of the North Island had small traces (approx. 2-3%) of M. m. castaneus nuclear ancestry, and mice in the upper South Island had approximately 7-8% M. m. musculus nuclear ancestry including some Y-chromosomal ancestry-though no detectable M. m. musculus mitochondrial ancestry. This is the most thorough genomic study of introduced populations of house mice yet conducted, and will have relevance to studies of the isolation mechanisms separating subspecies of mice.

Nuclear and mtDNA phylogenetic analyses clarify the evolutionary history of two species of native Hawaiian bats and the taxonomy of Lasiurini (Mammalia: Chiroptera)

Previous studies on genetics of hoary bats produced differing conclusions on the timing of their colonization of the Hawaiian Islands and whether or not North American (Aeorestes cinereus) and Hawaiian (A. semotus) hoary bats are distinct species. One study, using mtDNA COI and nuclear Rag2 and CMA1, concluded that hoary bats colonized the Hawaiian Islands no more than 10,000 years ago based on indications of population expansion at that time using Extended Bayesian Skyline Plots. The other study, using 3 mtDNA and 1 Y-chromosome locus, concluded that the Hawaiian Islands were colonized about 1 million years ago. To address the marked inconsistencies between those studies, we examined DNA sequences from 4 mitochondrial and 2 nuclear loci in lasiurine bats to investigate the timing of colonization of the Hawaiian Islands by hoary bats, test the hypothesis that Hawaiian and North American hoary bats belong to different species, and further investigate the generic level taxonomy within the tribe. Phylogenetic analysis and dating of the nodes of mtDNA haplotypes and of nuclear CMA1 alleles show that A. semotus invaded the Hawaiian Islands approximately 1.35 Ma and that multiple arrivals of A. cinereus occurred much more recently. Extended Bayesian Skyline plots show population expansion at about 20,000 years ago in the Hawaiian Islands, which we conclude does not represent the timing of colonization of the Hawaiian Islands given the high degree of genetic differentiation among A. cinereus and A. semotus (4.2% divergence at mtDNA Cytb) and the high degree of genetic diversity within A. semotus. Rather, population expansion 20,000 years ago could have resulted from colonization of additional islands, expansion after a bottleneck, or other factors. New genetic data also support the recognition of A. semotus and A. cinereus as distinct species, a finding consistent with previous morphological and behavioral studies. The phylogenetic analysis of CMA1 alleles shows the presence of 2 clades that are primarily associated with A. semotus mtDNA haplotypes, and are unique to the Hawaiian Islands. There is evidence for low levels of hybridization between A. semotus and A. cinereus on the Hawaiian Islands, but it is not extensive (<15% of individuals are of hybrid origin), and clearly each species is able to maintain its own genetic distinctiveness. Both mtDNA and nuclear DNA sequences show deep divergence between the 3 groups (genera) of lasiurine bats that correspond to the previously recognized morphological differences between them. We show that the Tribe Lasiurini contains the genera Aeorestes (hoary bats), Lasiurus (red bats), and Dasypterus (yellow bats).

A gene-tree test of the traditional taxonomy of American deer: the importance of voucher specimens, geographic data, and dense sampling

The taxonomy of American deer has been established almost entirely on the basis of morphological data and without the use of explicit phylogenetic methods; hence, phylogenetic analyses including data for all of the currently recognized species, even if based on a single gene, might improve current understanding of their taxonomy. We tested the monophyly of the morphology-defined genera and species of New World deer (Odocoileini) with phylogenetic analyses of mitochondrial DNA sequences. This is the first such test conducted using extensive geographic and taxonomic sampling. Our results do not support the monophyly of Mazama, Odocoileus, Pudu, M. americana, M. nemorivaga, Od. hemionus, and Od. virginianus. Mazama contains species that belong to other genera. We found a novel sister-taxon relationship between "Mazama" pandora and a clade formed by Od. hemionus columbianus and Od. h. sitkensis, and transfer pandora to Odocoileus. The clade formed by Od. h. columbianus and Od. h. sitkensis may represent a valid species, whereas the remaining subspecies of Od. hemionus appear closer to Od. virginianus. Pudu (Pudu) puda was not found sister to Pudu (Pudella) mephistophiles. If confirmed, this result will prompt the recognition of the monotypic Pudella as a distinct genus. We provide evidence for the existence of an undescribed species now confused with Mazama americana, and identify other instances of cryptic, taxonomically unrecognized species-level diversity among populations here regarded as Mazama temama, "Mazama" nemorivaga, and Hippocamelus antisensis. Noteworthy records that substantially extend the known distributions of M. temama and "M." gouazoubira are provided, and we unveil a surprising ambiguity regarding the distribution of "M." nemorivaga, as it is described in the literature. The study of deer of the tribe Odocoileini has been hampered by the paucity of information regarding voucher specimens and the provenance of sequences deposited in GenBank. We pinpoint priorities for future systematic research on the tribe Odocoileini.

Human-modified habitats change patterns of population genetic structure and group relatedness in Peter's tent-roosting bats

Although coloniality is widespread among mammals, it is still not clear what factors influence composition of social groups. As animals need to adapt to multiple habitat and environmental conditions throughout their range, variation in group composition should be influenced by adaptive adjustment to different ecological factors. Relevant to anthropogenic disturbance, increased habitat modification by humans can alter species' presence, density, and population structure. Therefore, it is important to understand the consequences of changes to landscape composition, in particular how habitat modification affects social structure of group-forming organisms. Here, we combine information on roosting associations with genetic structure of Peter's tent-roosting bats, Uroderma bilobatum to address how different habitat characteristics at different scales affect structure of social groups. By dividing analyses by age and sex, we determined that genetic structure was greater for adult females than adult males or offspring. Habitat variables explained 80% of the variation in group relatedness (mainly influenced by female relatedness) with roost characteristics contributing the most explained variation. This suggests that females using roosts of specific characteristics exhibit higher relatedness and seem to be philopatric. These females mate with more males than do more labile female groups. Results describe ecological and microevolutionary processes, which affect relatedness and social structure; findings are highly relevant to species distributions in both natural and human-modified environments.

Fine-scale genetic correlates to condition and migration in a wild cervid

The relationship between genetic variation and phenotypic traits is fundamental to the study and management of natural populations. Such relationships often are investigated by assessing correlations between phenotypic traits and heterozygosity or genetic differentiation. Using an extensive data set compiled from free-ranging mule deer (Odocoileus hemionus), we combined genetic and ecological data to (i) examine correlations between genetic differentiation and migration timing, (ii) screen for mitochondrial haplotypes associated with migration timing, and (iii) test whether nuclear heterozygosity was associated with condition. Migration was related to genetic differentiation (more closely related individuals migrated closer in time) and mitochondrial haplogroup. Body fat was related to heterozygosity at two nuclear loci (with antagonistic patterns), one of which is situated near a known fat metabolism gene in mammals. Despite being focused on a widespread panmictic species, these findings revealed a link between genetic variation and important phenotypes at a fine scale. We hypothesize that these correlations are either the result of mixing refugial lineages or differential mitochondrial haplotypes influencing energetics. The maintenance of phenotypic diversity will be critical to enable the potential tracking of changing climatic conditions, and these correlates highlight the need to consider evolutionary mechanisms in management, even in widely distributed panmictic species.

Timeframes of speciation, reticulation, and hybridization in the bulldog bat explained through phylogenetic analyses of all genetic transmission elements

Phylogenetic comparisons of the different mammalian genetic transmission elements (mtDNA, X-, Y-, and autosomal DNA) is a powerful approach for understanding the process of speciation in nature. Through such comparisons the unique inheritance pathways of each genetic element and gender-biased processes can link genomic structure to the evolutionary process, especially among lineages which have recently diversified, in which genetic isolation may be incomplete. Bulldog bats of the genus Noctilio are an exemplar lineage, being a young clade, widely distributed, and exhibiting unique feeding ecologies. In addition, currently recognized species are paraphyletic with respect to the mtDNA gene tree and contain morphologically identifiable clades that exhibit mtDNA divergences as great as among many species. To test taxonomic hypotheses and understand the contribution of hybridization to the extant distribution of genetic diversity in Noctilio, we used phylogenetic, coalescent stochastic modeling, and divergence time estimates using sequence data from cytochrome-b, cytochrome c oxidase-I, zinc finger Y, and zinc finger X, as well as evolutionary reconstructions based on amplified fragment length polymorphisms (AFLPs) data. No evidence of ongoing hybridization between the two currently recognized species was identified. However, signatures of an ancient mtDNA capture were recovered in which an mtDNA lineage of one species was captured early in the noctilionid radiation. Among subspecific mtDNA clades, which were generally coincident with morphology and statistically definable as species, signatures of ongoing hybridization were observed in sex chromosome sequences and AFLP. Divergence dating of genetic elements corroborates the diversification of extant Noctilio beginning about 3 Ma, with ongoing hybridization between mitochondrial lineages separated by 2.5 myr. The timeframe of species' divergence within Noctilio supports the hypothesis that shifts in the dietary strategies of gleaning insects (N. albiventris) or fish (N. leporinus) are among the most rapid instances of dietary evolution observed in mammals. This study illustrates the complex evolutionary dynamics shaping gene pools in nature, how comparisons of genetic elements can serve for understanding species boundaries, and the complex considerations for accurate taxonomic assignment.

An Account of the Taxonomy of North American Wolves From Morphological and Genetic Analyses

The available scientific literature was reviewed to assess the taxonomic standing of North American wolves, including subspecies of the gray wolf, Canis lupus. The recent scientific proposal that the eastern wolf, C. l. lycaon, is not a subspecies of gray wolf, but a full species, Canis lycaon, is well-supported by both morphological and genetic data. This species' range extends westward to Minnesota, and it hybridizes with gray wolves where the two species are in contact in eastern Canada and the Upper Peninsula of Michigan, Wisconsin, and Minnesota. Genetic data support a close relationship between eastern wolf and red wolf Canis rufus, but do not support the proposal that they are the same species; it is more likely that they evolved independently from different lineages of a common ancestor with coyotes. The genetic distinctiveness of the Mexican wolf Canis lupus baileyi supports its recognition as a subspecies. The available genetic and morphometric data do not provide clear support for the recognition of the Arctic wolf Canis lupus arctos, but the available genetic data are almost entirely limited to one group of genetic markers (microsatellite DNA) and are not definitive on this question. Recognition of the northern timber wolf Canis lupus occidentalis and the plains wolf Canis lupus nubilus as subspecies is supported by morphological data and extensive studies of microsatellite DNA variation where both subspecies are in contact in Canada. The wolves of coastal areas in southeastern Alaska and British Columbia should be assigned to C. lupus nubilus. There is scientific support for the taxa recognized here, but delineation of exact geographic boundaries presents challenges. Rather than sharp boundaries between taxa, boundaries should generally be thought of as intergrade zones of variable width.

Identification of novel single nucleotide polymorphisms (SNPs) in deer (Odocoileus spp.) using the BovineSNP50 BeadChip

Single nucleotide polymorphisms (SNPs) are growing in popularity as a genetic marker for investigating evolutionary processes. A panel of SNPs is often developed by comparing large quantities of DNA sequence data across multiple individuals to identify polymorphic sites. For non-model species, this is particularly difficult, as performing the necessary large-scale genomic sequencing often exceeds the resources available for the project. In this study, we trial the Bovine SNP50 BeadChip developed in cattle (Bos taurus) for identifying polymorphic SNPs in cervids Odocoileus hemionus (mule deer and black-tailed deer) and O. virginianus (white-tailed deer) in the Pacific Northwest. We found that 38.7% of loci could be genotyped, of which 5% (n = 1068) were polymorphic. Of these 1068 polymorphic SNPs, a mixture of putatively neutral loci (n = 878) and loci under selection (n = 190) were identified with the F_ST-outlier method. A range of population genetic analyses were implemented using these SNPs and a panel of 10 microsatellite loci. The three types of deer could readily be distinguished with both the SNP and microsatellite datasets. This study demonstrates that commercially developed SNP chips are a viable means of SNP discovery for non-model organisms, even when used between very distantly related species (the Bovidae and Cervidae families diverged some 25.1−30.1 million years before present).

Tracking of wisent-bison-yak mitochondrial evolution

One of the most informative sources which allow the drawing of far-reaching conclusions about the origins and phylogenetics of many species, including domestic animals and humans, is mitochondrial DNA (mtDNA). One of the important research targets should include the identification of similarities between wild and domestic species. The analysis involved the nucleotide sequences of mtDNA of wisent, auroch, bison, yak, bovine reference sequence (BRS) T3, T3a, T3b, T1, T1a, T1’2’3, T2, T3, T4, T5, Q, Q1, P, R, I1, and I2 bovine haplotypes. The non-coding D-loop regions were excluded from the evolutionary analysis and 15,419-bp coding sequences were used in the final dataset. Trees constructed on the basis of whole mitochondrial genomes or on total mtDNA coding sequences alignment were generally in agreement with previous studies on the Bovini tribe. American bison shows stronger maternal relationships to yak than to wisent. It seems that the isolation and divergence of wisent took place early, almost 2 to 1.6 million years ago. This appears to be compatible with the paleontological date, indicating Late Pleistocene speciation of Bison bonasus. The yak/bison mitochondrial transfer model is in agreement with our mutation analysis and phylogenetic tree. The bison/yak mutations were collected in the bison mitochondrial genome before the transfer. After the transfer, the parallel accumulation of unique mutations took place. According to our assessment, the transfer took place at about 700 ky. The characteristic feature of the wisent and bison evolution is the maintenance of mtDNA variability, despite the fact that both species underwent population bottlenecks. Our studies did not reveal any impact of these phenomena populations in the analyzed mitochondrial genomes.

Species-wide phylogeography of North American mule deer (Odocoileus hemionus): cryptic glacial refugia and postglacial recolonization

Quaternary climatic oscillations greatly influenced the present-day population genetic structure of animals and plants. For species with high dispersal and reproductive potential, phylogeographic patterns resulting from historical processes can be cryptic, overshadowed by contemporary processes. Here we report a study of the phylogeography of Odocoileus hemionus, a large, vagile ungulate common throughout western North America. We examined sequence variation of mitochondrial DNA (control region and cytochrome b) within and among 70 natural populations across the entire range of the species. Among the 1766 individual animals surveyed, we recovered 496 haplotypes. Although fine-scale phylogenetic structure was weakly resolved using phylogenetic methods, network analysis clearly revealed the presence of 12 distinct haplogroups. The spatial distribution of haplogroups showed a strong genetic discontinuity between the two morphological types of O. hemionus, mule deer and black-tailed deer, east and west of the Cascade Mountains in the Pacific Northwest. Within the mule deer lineage, we identified several haplogroups that expanded before or during the Last Glacial Maximum, suggesting that mule deer persisted in multiple refugia south of the ice sheets. Patterns of genetic diversity within the black-tailed deer lineage suggest a single refugium along the Pacific Northwest coast, and refute the hypothesis that black-tailed deer persisted in one or more northern refugia. Our data suggest that black-tailed deer recolonized areas in accordance with the pattern of glacial retreat, with initial recolonization northward along a coastal route and secondary recolonization inland.

Homogeneous differentiation of hepatocyte-like cells from embryonic stem cells: applications for the treatment of liver failure

One of the major hurdles of cellular therapies for the treatment of liver failure is the low availability of functional human hepatocytes. While embryonic stem (ES) cells represent a potential cell source for therapy, current methods for differentiation result in mixed cell populations or low yields of the cells of interest. Here we describe a rapid, direct differentiation method that yields a homogeneous population of endoderm-like cells with 95% purity. Mouse ES cells cultured on top of collagen-sandwiched hepatocytes differentiated and proliferated into a uniform and homogeneous cell population of endoderm-like cells. The endoderm-like cell population was positive for Foxa2, Sox17, and AFP and could be further differentiated into hepatocyte-like cells, demonstrating hepatic morphology, functionality, and gene and protein expression. Incorporating the hepatocyte-like cells into a bioartificial liver device to treat fulminant hepatic failure improved animal survival, thereby underscoring the therapeutic potential of these cells.

Development of Y chromosome intraspecific polymorphic markers in the Felidae

Y chromosome haplotyping based on microsatellites and single nucleotide polymorphisms (SNPs) has proved to be a powerful tool for population genetic studies of humans. However, the promise of the approach is hampered in the majority of nonhuman mammals by the lack of Y-specific polymorphic markers. We were able to identify new male-specific polymorphisms in the domestic cat Felis catus and 6 additional Felidae species with a combination of molecular genetic and cytogenetic approaches including 1) identifying domestic cat male-specific microsatellites from markers generated from a male cat microsatellite-enriched genomic library, a flow-sorted Y cosmid library, or a Y-specific cat bacteria artificial chromosome (BAC) clone, (2) constructing microsatellite-enriched libraries from flow-sorted Y chromosomes isolated directly from focal wildcat species, and (3) screening Y chromosome conserved anchored tagged sequences primers in Felidae species. Forty-one male-specific microsatellites were identified, but only 6 were single-copy loci, consistent with the repetitive nature of the Y chromosome. Nucleotide diversity (pi) of Y-linked intron sequences (2.1 kbp) was in the range of 0 (tiger) to 9.95 x 10(-4) (marbled cat), and the number of SNPs ranged from none in the tiger to 7 in the Asian leopard cat. The Y haplotyping system described here, consisting of 4 introns (SMCY3, SMCY7, UTY11, and DBY7) and 1 polymorphic microsatellite (SMCY-STR), represents the first available markers for tracking intraspecific male lineage polymorphisms in Felidae species and promises to provide significant insights to evolutionary and population genetic studies of the species.

A genomic perspective on the shortcomings of mitochondrial DNA for "barcoding" identification

Approximately 600-bp sequences of mitochondrial DNA (mtDNA) have been designated as "DNA barcodes" and have become one of the most contentious and animated issues in the application of genetic information to global biodiversity assessment and species identification. Advocates of DNA barcodes have received extensive attention and promotion in many popular and refereed scientific publications. However, we suggest that the utility of barcodes is suspect and vulnerable to technical challenges that are particularly pertinent to mtDNA. We review the natural history of mtDNA and discuss problems for barcoding which are particularly associated with mtDNA and inheritance, including reduced effective population size, maternal inheritance, recombination, inconsistent mutation rate, heteroplasmy, and compounding evolutionary processes. The aforementioned could significantly limit the application and utility of mtDNA barcoding efforts. Furthermore, global use of barcodes will require application and acceptance of a barcode-based species concept that has not been evaluated in the context of the extensive literature concerning species designation. Implementation of mtDNA barcodes in spite of technical and practical shortcomings we discuss may degrade the longstanding synthesis of genetic and organism-based research and will not advance studies ranging from genomic evolution to biodiversity assessment.

SPECIATION IN MAMMALS AND THE GENETIC SPECIES CONCEPT

We define a genetic species as a group of genetically compatible interbreeding natural populations that is genetically isolated from other such groups. This focus on genetic isolation rather than reproductive isolation distinguishes the Genetic Species Concept from the Biological Species Concept. Recognition of species that are genetically isolated (but not reproductively isolated) results in an enhanced understanding of biodiversity and the nature of speciation as well as speciation-based issues and evolution of mammals. We review criteria and methods for recognizing species of mammals and explore a theoretical scenario, the Bateson-Dobzhansky-Muller (BDM) model, for understanding and predicting genetic diversity and speciation in mammals. If the BDM model is operating in mammals, then genetically defined phylogroups would be predicted to occur within species defined by morphology, and phylogroups experiencing stabilizing selection will evolve genetic isolation without concomitant morphological diversification. Such species will be undetectable using classical skin and skull morphology (Morphological Species Concept). Using cytochrome-b data from sister species of mammals recognized by classical morphological studies, we estimated the number of phylogroups that exist within mammalian species and hypothesize that there will be >2,000 currently unrecognized species of mammals. Such an underestimation significantly affects conclusions on the nature of speciation in mammals, barriers associated with evolution of genetic isolation, estimates of biodiversity, design of conservation initiatives, zoonoses, and so on. A paradigm shift relative to this and other speciation-based issues will be needed. Data that will be effective in detecting these "morphologically cryptic genetic species" are genetic, especially DNA-sequence data. Application of the Genetic Species Concept uses genetic data from mitochondrial and nuclear genomes to identify species and species boundaries, the extent to which the integrity of the gene pool is protected, nature of hybridization (if present), and introgression. Genetic data are unique in understanding species because the use of genetic data 1) can quantify genetic divergence from different aspects of the genome (mitochondrial and nuclear genes, protein coding genes, regulatory genes, mobile DNA, microsatellites, chromosomal rearrangements, heterochromatin, etc.); 2) can provide divergence values that increase with time, providing an estimate of time since divergence; 3) can provide a population genetics perspective; 4) is less subject to convergence and parallelism relative to other sets of characters; 5) can identify monophyly, sister taxa, and presence or absence of introgression; and 6) can accurately identify hybrid individuals (kinship and source of hybrid individuals, F(1)s, backcrosses, direction of hybridization, and in concert with other data identify which hybrids are sterile or fertile). The proposed definition of the Genetic Species Concept is more compatible with a description of biodiversity of mammals than is "reproductively isolated species." Genetic profiles of mammalian species will result in a genetic description of species and mammalian diversity, and such studies are being accelerated by technological advances that reduce cost and increase speed and efficiency of generating genetic data. We propose that this genetic revolution remain museum- and voucher specimen-based and that new names are based on a holotype (including associated tissues) deposited in an accredited museum.

Hybridization in the recent past

The question we address in this article is how hybridization in the recent past can be detected in recently evolved species. Such species may not have evolved genetic incompatibilities and may hybridize with little or no fitness loss. Hybridization can be recognized by relatively small genetic differences between sympatric populations because sympatric populations have the opportunity to interbreed whereas allopatric populations do not. Using microsatellite DNA data from Darwin's finches in the Galapagos archipelago, we compare sympatric and allopatric genetic distances in pairs of Geospiza and Camarhynchus species. In agreement with the hybridization hypothesis, we found a statistically strong tendency for a species to be more similar genetically to a sympatric relative than to allopatric populations of that relative. Hybridization has been studied directly on two islands, but it is evidently more widespread in the archipelago. We argue that introgressive hybridization may have been a persistent feature of the adaptive radiation through most of its history, facilitating evolutionary diversification and occasionally affecting both the speed and direction of evolution.

mtDNA perspective of chromosomal diversification and hybridization in Peters' tent-making bat (Uroderma bilobatum: Phyllostomidae)

We compared sequence variation in the complete mitochondrial cytochrome-b gene with chromosomal and geographical variation for specimens of Peters' tent-making bat (Uroderma bilobatum). Three different chromosomal races have been described in this species: a 2n = 42 race from South America east of the Andes, a 2n = 44 from NW Central America and 2n = 38 from the rest of Central America and NW South America. The deepest nodes in the tree were found within the South American race (42 race), which is consistent with a longer history of this race. Average distance among races ranged from 2.5 to 2.9%, with the highest amount of intraracial variation found within the 2n = 42 race (1.7%), intermediate values within the 2n = 38 race (0.9%) and lowest within the 2n = 44 race (0.5%). Variation among chromosomal races accounted for over 55% of molecular variance, whereas variation among populations within races accounted for 6%. The 2n = 38 and 2n = 44 races hybridize in the coastal lowlands of Honduras, near the Gulf of Fonseca. Introgression between these two races is low (two introgressed individuals in 45 examined). Clinal variation across the hybrid zone for the cytochrome-b of U. bilobatum, is similar to clinal variation reported for chromosomes and isozymes of this species. Mismatch distribution analyses suggests that geographical isolation and karyological changes have interplayed in a synergistic fashion. Fixation of the alternative chromosomal rearrangements in geographical isolation and secondary contact is the most likely mechanism accounting for the hybrid zone between the 2n = 38 and 2n = 44 races. If a molecular clock is assumed, with rates ranging from 2.3 to 5.0% per million years, then isolation between these races occurred within the last million years, implying a relatively recent origin of the extant diversity in Uroderma bilobatum. None the less, the three chromosomal races probably represent three different biological species.

Phylogeographical structure, postglacial recolonization and barriers to gene flow in the distinctive Valais chromosome race of the common shrew (Sorex araneus)

Using one male-inherited and eight biparentally inherited microsatellite markers, we investigate the population genetic structure of the Valais chromosome race of the common shrew (Sorex araneus) in the Central Alps of Europe. Unexpectedly, the Y-chromosome microsatellite suggests nearly complete absence of male gene flow among populations from the St-Bernard and Simplon regions (Switzerland). Autosomal markers also show significant genetic structuring among these two geographical areas. Isolation by distance is significant and possible barriers to gene flow exist in the study area. Two different approaches are used to better understand the geographical patterns and the causes of this structuring. Using a principal component analysis for which testing procedure exists, and partial Mantel tests, we show that the St-Bernard pass does not represent a significant barrier to gene flow although it culminates at 2469 m, close to the highest altitudinal record for this species. Similar results are found for the Simplon pass, indicating that both passes represented potential postglacial recolonization routes into Switzerland from Italian refugia after the last Pleistocene glaciations. In contrast with the weak effect of these mountain passes, the Rhône valley lowlands significantly reduce gene flow in this species. Natural obstacles (the large Rhône river) and unsuitable habitats (dry slopes) are both present in the valley. Moreover, anthropogenic changes to landscape structures are likely to have strongly reduced available habitats for this shrew in the lowlands, thereby promoting genetic differentiation of populations found on opposite sides of the Rhône valley.

Haploid chromosomes in molecular ecology: lessons from the human Y

We review the potential use of haploid chromosomes in molecular ecology, using recent work on the human Y chromosome as a paradigm. Chromosomal sex-determination systems, and hence constitutively haploid chromosomes, which escape from recombination over much of their length, have evolved multiple times in the animal kingdom. In mammals, where males are the heterogametic sex, the patrilineal Y chromosome represents a paternal counterpart to mitochondrial DNA. Work on the human Y chromosome has shown it to contain the same range of polymorphic markers as the rest of the nuclear genome and these have rendered it the most informative haplotypic system in the human genome. Examples from research on the human Y chromosome are used to illustrate the common interests of anthropologists and ecologists in investigating the genetic impact of sex-specific behaviours and dispersals, as well as patterns of global diversity. We present some methodologies for extracting information from these uniquely informative yet under-utilized loci.

Characterization of a hybrid zone between two chromosomal races of the weta Hemideina thoracica following a geologically recent volcanic eruption

Two chromosomal races (2n=17 and 2n=15; XO) of the weta Hemideina thoracica meet at the centre of a volcanic region in North Island, New Zealand. Five independent polymorphic genetic markers showed broadly coinciding, steep frequency clines from north to south across this zone beside the flooded crater, Lake Taupo. Three unlinked nuclear gene markers provide estimates of zone width that are at least twice the width of the chromosomal and mitochondrial clines, with cline centres displaced at least 2.5 km. The different zone widths and centres suggest that this hybrid zone is a semipermeable barrier reducing the introgression of the chromosomal markers more than genic markers. We estimate that this species of weta must have a dispersal rate of at least 100 m per generation using the time since the last Taupo eruption (1850 years ago), which covered an area of about 20 000 km2 with pyroclastic flow.