Mammalian population genetics: why not Y?

Hybridization in an isolated population of blesbok and red hartebeest

Hybridization in antelope species has been widely reported in South African national parks and provincial reserves as well as on private land due to anthropogenic effects. In a closed management setting, hybridization may occur due to the crossbreeding of closely related species with unequal sex ratios, resulting in either sterile or fertile offspring. In this study, we used molecular techniques to evaluate the risk of anthropogenic hybridization between blesbok (Damaliscus pygargus phillipsi) and red hartebeest (Alcelaphus buselaphus caama) in an isolated group that purposely included the two species with unequal sex ratios (one red hartebeest male and 19 male and female blesbok). Genetic analysis based on microsatellites confirmed the presence of seven hybrid individuals. Mitochondrial analysis verified that hybridization occurred between blesbok females and the red hartebeest male. STRUCTURE and NEWHYBRIDS classified the hybrids as F1. It is suspected that the hybrid individuals were sterile as the males had undeveloped testes and only F1 hybrids were detected. Thus, the risk of hybridization between these two species may be limited in the wild. In captive settings, genetic monitoring should be included in management plans for blesbok and red hartebeest to ensure that the long-term consequences of wasted reproductive effort are limited.

Genomic evidence for two phylogenetic species and long-term population bottlenecks in red pandas

The red panda (Ailurus fulgens), an endangered Himalaya-endemic mammal, has been classified as two subspecies or even two species - the Himalayan red panda (A. fulgens) and the Chinese red panda (Ailurus styani) - based on differences in morphology and biogeography. However, this classification has remained controversial largely due to lack of genetic evidence, directly impairing scientific conservation management. Data from 65 whole genomes, 49 Y-chromosomes, and 49 mitochondrial genomes provide the first comprehensive genetic evidence for species divergence in red pandas, demonstrating substantial inter-species genetic divergence for all three markers and correcting species-distribution boundaries. Combined with morphological evidence, these data thus clearly define two phylogenetic species in red pandas. We also demonstrate different demographic trajectories in the two species: A. styani has experienced two population bottlenecks and one large population expansion over time, whereas A. fulgens has experienced three bottlenecks and one very small expansion, resulting in very low genetic diversity, high linkage disequilibrium, and high genetic load.

Accounting for kin sampling reveals genetic connectivity in Tasmanian and New Zealand school sharks, Galeorhinus galeus

Fishing represents a major problem for conservation of chondrichthyans, with a quarter of all species being overexploited. School sharks, Galeorhinus galeus, are targeted by commercial fisheries in Australia and New Zealand. The Australian stock has been depleted to below 20% of its virgin biomass, and the species is recorded as Conservation Dependent within Australia. Individuals are known to move between both countries, but it is disputed whether the stocks are reproductively linked. Accurate and unbiased determination of stock and population connectivity is crucial to inform effective management. In this study, we assess the genetic composition and population connectivity between Australian and New Zealand school sharks using genome-wide SNPs, while accounting for non-random kin sampling. Between 2009 and 2013, 88 neonate and juvenile individuals from Tasmanian and New Zealand nurseries were collected and genotyped. Neutral loci were analyzed to detect fine-scale signals of reproductive connectivity. Seven full-sibling groups were identified and removed for unbiased analysis. Based on 6,587 neutral SNPs, pairwise genetic differentiation from Tasmanian and New Zealand neonates was non-significant (F ST = 0.0003, CI95 = [-0.0002, 0.0009], p = 0.1163; D est = 0.0006 ± 0.0002). This pattern was supported by clustering results. In conclusion, we show a significant effect of non-random sampling of kin and identify fine-scale reproductive connectivity between Australian and New Zealand school sharks.

OPEN RESEARCH BADGES

This article has earned an Open Data Badge for making publicly available the digitally-shareable data necessary to reproduce the reported results. The data is available at https://doi.org/10.5061/dryad.pd8612j.

Paternal phylogeographic structure of the brown bear (Ursus arctos) in northeastern Asia and the effect of male-mediated gene flow to insular populations

BACKGROUND

Sex-biased dispersal is widespread among mammals, including the brown bear (Ursus arctos). Previous phylogeographic studies of the brown bear based on maternally inherited mitochondrial DNA have shown intraspecific genetic structuring around the northern hemisphere. The brown bears on Hokkaido Island, northern Japan, comprise three distinct maternal lineages that presumably immigrated to the island from the continent in three different periods. Here, we investigate the paternal genetic structure across northeastern Asia and assess the connectivity among and within intraspecific populations in terms of male-mediated gene flow.

RESULTS

We analyzed paternally inherited Y-chromosomal DNA sequence data and Y-linked microsatellite data of 124 brown bears from Hokkaido, the southern Kuril Islands (Kunashiri and Etorofu), Sakhalin, and continental Eurasia (Kamchatka Peninsula, Ural Mountains, European Russia, and Tibet). The Hokkaido brown bear population is paternally differentiated from, and lacked recent genetic connectivity with, the continental Eurasian and North American populations. We detected weak spatial genetic structuring of the paternal lineages on Hokkaido, which may have arisen through male-mediated gene flow among natal populations. In addition, our results suggest that the different dispersal patterns between male and female brown bears, combined with the founder effect and subsequent genetic drift, contributed to the makeup of the Etorofu Island population, in which the maternal and paternal lineages show different origins.

CONCLUSIONS

Brown bears on Hokkaido and the adjacent southern Kuril Islands experienced different maternal and paternal evolutionary histories. Our results indicate that sex-biased dispersal has played a significant role in the evolutionary history of the brown bear in continental populations and in peripheral insular populations, such as on Hokkaido, the southern Kuril Islands, and Sakhalin.

The impact of mating systems and dispersal on fine-scale genetic structure at maternally, paternally and biparentally inherited markers

For decades, studies have focused on how dispersal and mating systems influence genetic structure across populations or social groups. However, we still lack a thorough understanding of how these processes and their interaction shape spatial genetic patterns over a finer scale (tens-hundreds of metres). Using uniparentally inherited markers may help answer these questions, yet their potential has not been fully explored. Here, we use individual-level simulations to investigate the effects of dispersal and mating system on fine-scale genetic structure at autosomal, mitochondrial and Y chromosome markers. Using genetic spatial autocorrelation analysis, we found that dispersal was the major driver of fine-scale genetic structure across maternally, paternally and biparentally inherited markers. However, when dispersal was restricted (mean distance = 100 m), variation in mating behaviour created strong differences in the comparative level of structure detected at maternally and paternally inherited markers. Promiscuity reduced spatial genetic structure at Y chromosome loci (relative to monogamy), whereas structure increased under polygyny. In contrast, mitochondrial and autosomal markers were robust to differences in the specific mating system, although genetic structure increased across all markers when reproductive success was skewed towards fewer individuals. Comparing males and females at Y chromosome vs. mitochondrial markers, respectively, revealed that some mating systems can generate similar patterns to those expected under sex-biased dispersal. This demonstrates the need for caution when inferring ecological and behavioural processes from genetic results. Comparing patterns between the sexes, across a range of marker types, may help us tease apart the processes shaping fine-scale genetic structure.

Inner and inter population structure construction of Chinese Jiangsu Han population based on Y23 STR system

In this study, we analyzed the genetic polymorphisms of 23 Y-STR loci from PowerPlex^® Y23 system in 916 unrelated healthy male individuals from Chinese Jiangsu Han, and observed 912 different haplotypes including 908 unique haplotypes and 4 duplicate haplotypes. The haplotype diversity reached 0.99999 and the discrimination capacity and match probability were 0.9956 and 0.0011, respectively. The gene diversity values ranged from 0.3942 at DYS438 to 0.9607 at DYS385a/b. Population differentiation within 10 Jiangsu Han subpopulations were evaluated by R_ST values and visualized in Neighbor-Joining trees and Multi-Dimensional Scaling plots as well as population relationships between the Jiangsu Han population and other 18 Eastern Asian populations. Such results indicated that the 23 Y-STR loci were highly polymorphic in Jiangsu Han population and played crucial roles in forensic application as well as population genetics. For the first time, we reported the genetic diversity of male lineages in Jiangsu Han population at a high-resolution level of 23 Y-STR set and consequently contributed to familial searching, offender tracking, and anthropology analysis of Jiangsu Han population.

New insights into the genetic composition and phylogenetic relationship of wolves and dogs in the Iberian Peninsula

This study investigates the gene pool of Portuguese autochthonous dog breeds and their wild counterpart, the Iberian wolf subspecies (Canis lupus signatus), using standard molecular markers. A combination of paternal and maternal molecular markers was used to investigate the genetic composition, genetic differentiation and genetic relationship of native Portuguese dogs and the Iberian wolf. A total of 196 unrelated dogs, including breed and village dogs from Portugal, and other dogs from Spain and North Africa, and 56 Iberian wolves (wild and captive) were analyzed for nuclear markers, namely Y chromosome SNPs, Y chromosome STR loci, autosomal STR loci, and a mitochondrial fragment of the control region I. Our data reveal new variants for the molecular markers and confirm significant genetic differentiation between Iberian wolf and native domestic dogs from Portugal. Based on our sampling, no signs of recent introgression between the two subspecies were detected. Y chromosome data do not reveal genetic differentiation among the analyzed dog breeds, suggesting they share the same patrilineal origin. Moreover, the genetic distinctiveness of the Iberian wolf from other wolf populations is further confirmed with the description of new mtDNA variants for this endemism. Our research also discloses new molecular markers for wolf and dog subspecies assignment, which might become particularly relevant in the case of forensic or noninvasive genetic studies. The Iberian wolf represents a relic of the once widespread wolf population in Europe and our study reveals that it is a reservoir of unique genetic diversity of the grey wolf, Canis lupus. These results stress the need for conservation plans that will guarantee the sustainability of this threatened top predator in Iberia.

Population genetic insights into the social organization of Guinea baboons (Papio papio): Evidence for female-biased dispersal

Sex differences in philopatry and dispersal have important consequences on the genetic structure of populations, social groups, and social relationships within groups. Among mammals, male dispersal and female philopatry are most common and closely related taxa typically exhibit similar dispersal patterns. However, among four well-studied species of baboons, only hamadryas baboons exhibit female dispersal, thus differing from their congenerics, which show female philopatry and close-knit female social relationships. Until recently, knowledge of the Guinea baboon social system and dispersal pattern remained sparse. Previous observations suggested that the high degree of tolerance observed among male Guinea baboons could be due to kinship. This led us to hypothesize that this species exhibits male philopatry and female dispersal, conforming to the hamadryas pattern. We genotyped 165 individuals from five localities in the Niokolo-Koba National Park, Senegal, at 14 autosomal microsatellite loci and sequenced a fragment of the mitochondrial hypervariable region I (HVRI) of 55 individuals. We found evidence for higher population structuring in males than in females, as expected if males are the more philopatric sex. A comparison of relatedness between male-male and female-female dyads within and among communities did not yield conclusive results. HVRI diversity within communities was high and did not differ between the sexes, also suggesting female gene flow. Our study is the first comprehensive analysis of the genetic population structure in Guinea baboons and provides evidence for female-biased dispersal in this species. In conjunction with their multilevel social organization, this finding parallels the observations for human hunter-gatherers and strengthens baboons as an intriguing model to elucidate the processes that shaped the highly cooperative societies of Homo.

Y-chromosomal variation of local goat breeds of Turkey close to the domestication centre

Genetic variations in chromosome Y are enabling researchers to identify paternal lineages, which are informative for introgressions and migrations. In this study, the male-specific region markers, sex-determining region-Y (SRY), amelogenin (AMELY) and zinc finger (ZFY) were analysed in seven Turkish native goat breeds, Angora, Kilis, Hair, Honamlı, Norduz, Gürcü and Abaza. A SNP in the ZFY gene defined a new haplotype Y2C. All domestic haplogroups originate from Capra aegagrus, while the finding of Y1A, Y1B, Y2A and Y2C in 32, 4, 126 and 2 Turkish domestic goats, respectively, appears to indicate a predomestic origin of the major haplotypes. The occurrence of four haplotypes in the Hair goat and, in contrast, a frequency of 96% of Y1A in the Kilis breed illustrate that Y-chromosomal variants have a more breed-dependent distribution than mitochondrial or autosomal DNA. This probably reflects male founder effects, but a role in adaptation cannot be excluded.

Bears in a forest of gene trees: phylogenetic inference is complicated by incomplete lineage sorting and gene flow

Ursine bears are a mammalian subfamily that comprises six morphologically and ecologically distinct extant species. Previous phylogenetic analyses of concatenated nuclear genes could not resolve all relationships among bears, and appeared to conflict with the mitochondrial phylogeny. Evolutionary processes such as incomplete lineage sorting and introgression can cause gene tree discordance and complicate phylogenetic inferences, but are not accounted for in phylogenetic analyses of concatenated data. We generated a high-resolution data set of autosomal introns from several individuals per species and of Y-chromosomal markers. Incorporating intraspecific variability in coalescence-based phylogenetic and gene flow estimation approaches, we traced the genealogical history of individual alleles. Considerable heterogeneity among nuclear loci and discordance between nuclear and mitochondrial phylogenies were found. A species tree with divergence time estimates indicated that ursine bears diversified within less than 2 My. Consistent with a complex branching order within a clade of Asian bear species, we identified unidirectional gene flow from Asian black into sloth bears. Moreover, gene flow detected from brown into American black bears can explain the conflicting placement of the American black bear in mitochondrial and nuclear phylogenies. These results highlight that both incomplete lineage sorting and introgression are prominent evolutionary forces even on time scales up to several million years. Complex evolutionary patterns are not adequately captured by strictly bifurcating models, and can only be fully understood when analyzing multiple independently inherited loci in a coalescence framework. Phylogenetic incongruence among gene trees hence needs to be recognized as a biologically meaningful signal.

Brown and polar bear Y chromosomes reveal extensive male-biased gene flow within brother lineages

Brown and polar bears have become prominent examples in phylogeography, but previous phylogeographic studies relied largely on maternally inherited mitochondrial DNA (mtDNA) or were geographically restricted. The male-specific Y chromosome, a natural counterpart to mtDNA, has remained underexplored. Although this paternally inherited chromosome is indispensable for comprehensive analyses of phylogeographic patterns, technical difficulties and low variability have hampered its application in most mammals. We developed 13 novel Y-chromosomal sequence and microsatellite markers from the polar bear genome and screened these in a broad geographic sample of 130 brown and polar bears. We also analyzed a 390-kb-long Y-chromosomal scaffold using sequencing data from published male ursine genomes. Y chromosome evidence support the emerging understanding that brown and polar bears started to diverge no later than the Middle Pleistocene. Contrary to mtDNA patterns, we found 1) brown and polar bears to be reciprocally monophyletic sister (or rather brother) lineages, without signals of introgression, 2) male-biased gene flow across continents and on phylogeographic time scales, and 3) male dispersal that links the Alaskan ABC islands population to mainland brown bears. Due to female philopatry, mtDNA provides a highly structured estimate of population differentiation, while male-biased gene flow is a homogenizing force for nuclear genetic variation. Our findings highlight the importance of analyzing both maternally and paternally inherited loci for a comprehensive view of phylogeographic history, and that mtDNA-based phylogeographic studies of many mammals should be reevaluated. Recent advances in sequencing technology render the analysis of Y-chromosomal variation feasible, even in nonmodel organisms.

Sex-linked and autosomal microsatellites provide new insights into island populations of the tammar wallaby

The emerging availability of microsatellite markers from mammalian sex chromosomes provides opportunities to investigate both male- and female-mediated gene flow in wild populations, identifying patterns not apparent from the analysis of autosomal markers alone. Tammar wallabies (Macropus eugenii), once spread over the southern mainland, have been isolated on several islands off the Western Australian and South Australian coastlines for between 10,000 and 13,000 years. Here, we combine analyses of autosomal, Y-linked and X-linked microsatellite loci to investigate genetic variation in populations of this species on two islands (Kangaroo Island, South Australia and Garden Island, Western Australia). All measures of diversity were higher for the larger Kangaroo Island population, in which genetic variation was lowest at Y-linked markers and highest at autosomal markers (θ=3.291, 1.208 and 0.627 for autosomal, X-linked and Y-linked data, respectively). Greater relatedness among females than males provides evidence for male-biased dispersal in this population, while sex-linked markers identified genetic lineages not apparent from autosomal data alone. Overall genetic diversity in the Garden Island population was low, especially on the Y chromosome where most males shared a common haplotype, and we observed high levels of inbreeding and relatedness among individuals. Our findings highlight the utility of this approach for management actions, such as the selection of animals for translocation or captive breeding, and the ecological insights that may be gained by combining analyses of microsatellite markers on sex chromosomes with those derived from autosomes.

Sex-specific clines support incipient speciation in a common European mammal

Hybrid zones provide excellent opportunities to study processes and mechanisms underlying reproductive isolation and speciation. Here we investigated sex-specific clines of molecular markers in hybrid zones of morphologically cryptic yet genetically highly-diverged evolutionary lineages of the European common vole (Microtus arvalis). We analyzed the position and width of four secondary contact zones along three independent transects in the region of the Alps using maternally (mitochondrial DNA) and paternally (Y-chromosome) inherited genetic markers. Given male-biased dispersal in the common vole, a selectively neutral secondary contact would show broader paternal marker clines than maternal ones. In a selective case, for example, involving a form of Haldane's rule, Y-chromosomal clines would not be expected to be broader than maternal markers because they are transmitted by the heterogametic sex and thus gene flow would be restricted. Consistent with the selective case, paternal clines were significantly narrower or at most equal in width to maternal clines in all contact zones. In addition, analyses using maximum likelihood cline-fitting detected a shift of paternal relative to maternal clines in three of four contact zones. These patterns suggest that processes at the contact zones in the common vole are not selectively neutral, and that partial reproductive isolation is already established between these evolutionary lineages. We conclude that hybrid zone movement, sexual selection and/or genetic incompatibilities are likely associated with an unusual unidirectional manifestation of Haldane's rule in this common European mammal.

Paternally inherited genetic markers reveal new insights into genetic structuring within Macropus fuliginosus and hybridisation with sympatric Macropus giganteus

There are several aspects of biology in which the contribution of males and females is unequal. In these instances the examination of Y chromosome markers may be used to elucidate male-specific attributes. Here, male dispersal patterns and genetic structuring were examined using four Y-microsatellite loci in 186 male western grey kangaroos, Macropus fuliginosus, from throughout the species’ trans-continental distribution. In addition, 52 male grey kangaroos were examined to investigate hybridisation between M. fuliginosus and the eastern grey kangaroo, Macropus giganteus, in their region of sympatry in eastern Australia. Detected Y chromosome diversity was low, resulting from low effective male population size due to skewed sex ratios and a polygynous mating system. As expected, male dispersal was high across the range. However, the Lake Torrens–Flinders Ranges region appears to have significantly restricted male movement between eastern and central/western Australia. There was little evidence to suggest that other barriers (Nullarbor Plain and Swan River Valley) previously identified by nuclear and mitochondrial DNA marker studies restrict male movement. Hence, the admixture events previously identified may be associated with high male dispersal. Within the region of sympatry between M. fuliginosus and M. giganteus in eastern Australia, four M. giganteus individuals were found to possess M. fuliginosus Y-haplotypes. These results confirm the occurrence of hybridisation between male M. fuliginosus and female M. giganteus. Additionally, the introgression of M. fuliginosus Y-haplotypes into M. giganteus populations indicates that at least some male hybrids are fertile, despite evidence to the contrary from captive studies. This study has provided insights into the male contribution to population history, structure and hybridisation in M. fuliginosus, which were not predicted by comparisons between biparentally and maternally inherited markers. This highlights the importance of direct examination of the Y chromosome to provide novel insights into male-mediated processes, especially where the contribution of the sexes may differ.

Spatial framework of nine distinct local populations of the Japanese dormouse Glirulus japonicus based on matrilineal cytochrome b and patrilineal SRY gene sequences

We previously revealed the presence of six genetically distinct matrilineal populations of the Japanese dormouse Glirulus japonicus in the distribution range of Honshu, Shikoku, and Kyushu islands. In this study, we extended this analysis using mitochondrial cytochrome b gene sequences (n = 96) and Y-chromosome-specific SRY gene sequences (n = 22) from individuals collected from Honshu, Shikoku, Kyushu, and Oki Dogo I. The cytochrome b sequence data allowed us to define precise geographic ranges of the six previously known and three newly found distinct matrilineal lineages: northeastern Honshu (I), east-central Honshu (II), west-central Honshu and the Kii Peninsula (III), the western part of Honshu (IV), Shikoku (V), westernmost Honshu and Kyushu (VI), the northern part of central Honshu (VII), the southern part of central Honshu (VIII), and Oki Dogo I. (IX). Our inference of geographic borders suggests that regions of lower and higher altitudes in the mountain systems played important roles in driving the hosting and separation of lineages, respectively. Six matrilineal lineages (I, II, V, VI, VIII, and XI) were shown to possess their own SRY haplotypes, while lineages III and IV shared one haplotype. These data together with our previous observation of nuclear ribosomal RNA gene variation indicate advanced populational subdivision in this species. It is thus evident that each of the populations, including those living at high latitudes and in limited geographic spaces, have survived for several million years. A specific ability to tolerate cold may have permitted G. japonicus to preserve anciently diverged lineages in each locality.

Y-chromosome phylogeny in the evolutionary net of chamois (genus Rupicapra)

Background

The chamois, distributed over most of the medium to high altitude mountain ranges of southern Eurasia, provides an excellent model for exploring the effects of historical and evolutionary events on diversification. Populations have been grouped into two species, Rupicapra pyrenaica from southwestern Europe and R. rupicapra from eastern Europe. The study of matrilineal mitochondrial DNA (mtDNA) and biparentally inherited microsatellites showed that the two species are paraphyletic and indicated alternate events of population contraction and dispersal-hybridization in the diversification of chamois. Here we investigate the pattern of variation of the Y-chromosome to obtain information on the patrilineal phylogenetic position of the genus Rupicapra and on the male-specific dispersal of chamois across Europe.

Results

We analyzed the Y-chromosome of 87 males covering the distribution range of the Rupicapra genus. We sequenced a fragment of the SRY gene promoter and characterized the male specific microsatellites UMN2303 and SRYM18. The SRY promoter sequences of two samples of Barbary sheep (Ammotragus lervia) were also determined and compared with the sequences of Bovidae available in the GenBank. Phylogenetic analysis of the alignment showed the clustering of Rupicapra with Capra and the Ammotragus sequence obtained in this study, different from the previously reported sequence of Ammotragus which groups with Ovis. Within Rupicapra, the combined data define 10 Y-chromosome haplotypes forming two haplogroups, which concur with taxonomic classification, instead of the three clades formed for mtDNA and nuclear microsatellites. The variation shows a west-to-east geographical cline of ancestral to derived alleles.

Conclusions

The phylogeny of the SRY-promoter shows an association between Rupicapra and Capra. The position of Ammotragus needs a reinvestigation. The study of ancestral and derived characters in the Y-chromosome suggests that, contrary to the presumed Asian origin, the paternal lineage of chamois originated in the Mediterranean, most probably in the Iberian Peninsula, and dispersed eastwards through serial funding events during the glacial-interglacial cycles of the Quaternary. The diversity of Y-chromosomes in chamois is very low. The differences in patterns of variation among Y-chromosome, mtDNA and biparental microsatellites reflect the evolutionary characteristics of the different markers as well as the effects of sex-biased dispersal and species phylogeography.

Identifying Y-chromosomal diversity by long-template PCR

Comparing Y-chromosomal and mitochondrial haplotype variation is a promising approach to independently investigate paternal and maternal evolutionary histories in wild mammal populations. However, the difficulty of developing male-specific genetic markers, because of its distinctive genetic architecture and the general low level of polymorphisms observed on the Y chromosome, hampers usually an effective application of this approach. Here, we present a further method of the established Y chromosome conserved anchored tagged sequences strategy to develop Y-chromosomal markers by screening introns of male-specific region (MSY) genes for sequence polymorphisms. By applying long-template PCR using target species-specific primers, adequate sequence information of several kb in size can be obtained. We applied this method in the snow vole (Chionomys nivalis) and obtained 12.4 kb of male-specific sequence data for nine males representing four populations in the Swiss Alps. A total of 28 single nucleotide polymorphisms, four indels (> 1 bp) and one polymorphic microsatellite were identified in introns of the SMCY and DBY genes. Based on this information, we developed a Y-chromosomal genotyping assay and identified four different paternal lineages within one local snow vole population. The method we present is straightforward and as such will probably be suitable to detect adequate Y-chromosomal diversity in a wide range of mammalian species.

The quest for Y-chromosomal markers - methodological strategies for mammalian non-model organisms

Tracing maternal and paternal lineages independently to explore breeding systems and dispersal strategies in natural populations has been high on the wish-list of evolutionary biologists. As males are the heterogametic sex in mammals, such sex-specific patterns can be indirectly observed when Y chromosome polymorphism is combined with mitochondrial sequence information. Over the past decade, Y-chromosomal markers applied to human populations have revealed remarkable differences in the demographic history and behaviour between the sexes. However, with a few exceptions, genetic data tracing the paternal line are lacking in most other mammalian species. This deficit can be attributed to the difficulty of developing Y-specific genetic markers in non-model organisms and the general low levels of polymorphisms observed on the Y chromosome. Here, we present an overview of the currently employed strategies for developing paternal markers in mammals. Moreover, we review the practical feasibility and requirements of various methodological strategies and highlight their future prospects when combined with new molecular techniques such as next generation sequencing.

A highly polymorphic insertion in the Y-chromosome amelogenin gene can be used for evolutionary biology, population genetics and sexing in Cetacea and Artiodactyla

Background

The early radiation of the Cetartiodactyla is complex, and unambiguous molecular characters are needed to clarify the positions of hippotamuses, camels and pigs relative to the remaining taxa (Cetacea and Ruminantia). There is also a need for informative genealogic markers for Y-chromosome population genetics as well as a sexing method applicable to all species from this group. We therefore studied the sequence variation of a partial sequence of the evolutionary conserved amelogenin gene to assess its potential use in each of these fields.

Results and discussion

We report a large interstitial insertion in the Y amelogenin locus in most of the Cetartiodactyla lineages (cetaceans and ruminants). This sex-linked size polymorphism is the result of a 460–465 bp inserted element in intron 4 of the amelogenin gene of Ruminants and Cetaceans. Therefore, this polymorphism can easily be used in a sexing assay for these species.

When taking into account this shared character in addition to nucleotide sequence, gene genealogy follows sex-chromosome divergence in Cetartiodactyla whereas it is more congruent with zoological history when ignoring these characters. This could be related to a loss of homology between chromosomal copies given the old age of the insertion.

The 1 kbp Amel-Y amplified fragment is also characterized by high nucleotide diversity (64 polymorphic sites spanning over 1 kbp in seven haplotypes) which is greater than for other Y-chromosome sequence markers studied so far but less than the mitochondrial control region.

Conclusion

The gender-dependent polymorphism we have identified is relevant not only for phylogenic inference within the Cetartiodactyla but also for Y-chromosome based population genetics and gender determination in cetaceans and ruminants. One single protocol can therefore be used for studies in population and evolutionary genetics, reproductive biotechnologies, and forensic science.

Sex-biased evolutionary forces shape genomic patterns of human diversity

Comparisons of levels of variability on the autosomes and X chromosome can be used to test hypotheses about factors influencing patterns of genomic variation. While a tremendous amount of nucleotide sequence data from across the genome is now available for multiple human populations, there has been no systematic effort to examine relative levels of neutral polymorphism on the X chromosome versus autosomes. We analyzed ∼210 kb of DNA sequencing data representing 40 independent noncoding regions on the autosomes and X chromosome from each of 90 humans from six geographically diverse populations. We correct for differences in mutation rates between males and females by considering the ratio of within-human diversity to human-orangutan divergence. We find that relative levels of genetic variation are higher than expected on the X chromosome in all six human populations. We test a number of alternative hypotheses to explain the excess polymorphism on the X chromosome, including models of background selection, changes in population size, and sex-specific migration in a structured population. While each of these processes may have a small effect on the relative ratio of X-linked to autosomal diversity, our results point to a systematic difference between the sexes in the variance in reproductive success; namely, the widespread effects of polygyny in human populations. We conclude that factors leading to a lower male versus female effective population size must be considered as important demographic variables in efforts to construct models of human demographic history and for understanding the forces shaping patterns of human genomic variability.

Like many primate species, the mating system of humans is considered to be moderately polygynous (i.e., males exhibit a higher variance in reproductive success than females). As a consequence, males are expected to have a lower effective population size (N_e) than females, and the proportion of neutral genetic variation on the X chromosome (relative to the autosomes) should be higher than expected under the assumption of strict neutrality and an equal breeding sex ratio. We test for the effects of polygyny by measuring levels of neutral polymorphism at 40 independent loci on the X chromosome and autosomes in six human populations. To correct for mutation rate heterogeneity among loci, we divide our diversity estimates within human populations by divergence with orangutan at each locus. Consistent with expectations under a model of polygyny, we find elevated levels of X-linked versus autosomal diversity. While it is possible that multiple demographic processes may contribute to the observed patterns of genomic diversity (i.e., background selection, changes in population size, and sex-specific migration), we conclude that an historical excess of breeding females over the number of breeding males can by itself explain most of the observed increase in effective population size of the X chromosome.

Differentiation of tundra/taiga and boreal coniferous forest wolves: genetics, coat colour and association with migratory caribou

The grey wolf has one of the largest historic distributions of any terrestrial mammal and can disperse over great distances across imposing topographic barriers. As a result, geographical distance and physical obstacles to dispersal may not be consequential factors in the evolutionary divergence of wolf populations. However, recent studies suggest ecological features can constrain gene flow. We tested whether wolf-prey associations in uninterrupted tundra and forested regions of Canada explained differences in migratory behaviour, genetics, and coat colour of wolves. Satellite-telemetry data demonstrated that tundra wolves (n = 19) migrate annually with caribou (n = 19) from denning areas in the tundra to wintering areas south of the treeline. In contrast, nearby boreal coniferous forest wolves are territorial and associated year round with resident prey. Spatially explicit analysis of 14 autosomal microsatellite loci (n = 404 individuals) found two genetic clusters corresponding to tundra vs. boreal coniferous forest wolves. A sex bias in gene flow was inferred based on higher levels of mtDNA divergence (F(ST) = 0.282, 0.028 and 0.033; P < 0.0001 for mitochondrial, nuclear autosomal and Y-chromosome markers, respectively). Phenotypic differentiation was substantial as 93% of wolves from tundra populations exhibited light colouration whereas only 38% of boreal coniferous forest wolves did (chi(2) = 64.52, P < 0.0001). The sharp boundary representing this discontinuity was the southern limit of the caribou migration. These findings show that substantial genetic and phenotypic differentiation in highly mobile mammals can be caused by prey-habitat specialization rather than distance or topographic barriers. The presence of a distinct wolf ecotype in the tundra of North America highlights the need to preserve migratory populations.

Sex-biased dispersal in western lowland gorillas (Gorilla gorilla gorilla)

We explored two hypotheses related to potential differences between sexes in dispersal behaviour in western lowland gorillas (Gorilla gorilla gorilla). Direct observations suggest that immature females have more opportunities to move between breeding groups than immature males. The distribution of kin dyadic relationships within and between groups does not, however, support this hypothesis. At larger geographical scales, dispersal is likely to be easier for males than females because of the solitary phase most blackbacks experience before founding their own breeding group. However, previous work indicates that males settle preferentially close to male kin. By specifically tracing female and male lineages with mitochondrial and Y-chromosomal genetic markers, we found that male gorillas in the 6000 km2 area we surveyed form a single population whereas females are restricted to the individual sites we sampled and do not freely move around this area. These differences are more correctly described as differences in dispersal distances, rather than differences in dispersal rates between sexes (both sexes emigrate from their natal group in this species). Differences in resource competition and dispersal costs between female and male gorillas are compatible with the observed pattern, but more work is needed to understand if these ultimate causes are responsible for sex-biased dispersal distances in western lowland gorillas.

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.

Dispersal is not female biased in a resource-defence mating ungulate, the European roe deer

Dispersal is frequently more prevalent in one sex compared to the other. Greenwood proposed that patterns of sex-biased dispersal among birds and mammals are linked to their mating strategies. For species where males defend resources rather than females, he predicted female-biased dispersal, because males should remain at their birth site where they are familiar with the distribution of the resources that they must defend. Greenwood's hypothesis has been extensively supported among birds, where most species exhibit a resource-defence mating strategy. However, almost no equivalent information is available for mammals as males generally defend mates in this group. An exception is the European roe deer, a resource-defence mating ungulate. We thus tested Greenwood's hypothesis on this atypical mammalian model, looking for female-biased dispersal using sex-specific inter-individual genetic distances. We conclusively show that gene flow is not higher among females compared to males in the studied roe deer population, and hence that dispersal is not female-biased, suggesting that male mating strategy is not the primary selective force driving the evolution of dispersal in roe deer. We discuss the role of female mate choice and intra-sexual competition as possible alternative selective pressures involved.

Disentangling reasons for low Y chromosome variation in the greater white-toothed shrew (Crocidura russula)

Y chromosome variation is determined by several confounding factors including mutation rate, effective population size, demography, and selection. Disentangling these factors is essential to better understand the evolutionary properties of the Y chromosome. We analyzed genetic variation on the Y chromosome, X chromosome, and mtDNA of the greater white-toothed shrew, a species with low variance in male reproductive success and limited sex-biased dispersal, which enables us to control to some extent for life-history effects. We also compared ancestral (Moroccan) to derived (European) populations to investigate the role of demographic history in determining Y variation. Recent colonization of Europe by a small number of founders (combined with low mutation rates) is largely responsible for low diversity observed on the European Y and X chromosomes compared to mtDNA. After accounting for mutation rate, copy number, and demography, the Y chromosome still displays a deficit in variation relative to the X in both populations. This is possibly influenced by directional selection, but the slightly higher variance in male reproductive success is also likely to play a role, even though the difference is small compared to that in highly polygynous species. This study illustrates that demography and life-history effects should be scrutinized before inferring strong selective pressure as a reason for low diversity on the Y chromosome.

Low Y chromosome variation in Saudi-Arabian hamadryas baboons (Papio hamadryas hamadryas)

It is important to characterise the amount of variation on the mammalian Y chromosome in order to assess its potential for use in evolutionary studies. We report very low levels of polymorphism on the Y chromosome of Saudi-Arabian hamadryas baboons, Papio hamadryas hamadryas. We found no segregating sites on the Y, despite sequence analysis of 3 kb noncontiguous intron sequence in 16 males with divergent autosomal microsatellite genotypes, and a further analysis of 1.1 kb intron sequence in 97 males from four populations by SSCP. In addition, we tested seven human-derived Y-linked microsatellites in baboons. Only four of these loci were male-specific and only one was polymorphic in our 97 male sample set. Polymorphism on the Y chromosome of Arabian hamadryas appears to be low compared to other primate species for which data are available (eg humans, chimpanzees and bonobos). Low effective population size (Ne) of paternal genes due to polygyny and female-biased adult sex ratio is a potential reason for low Y chromosome variation in this species. However, low Ne for the Y should be counterbalanced to some extent by the species' atypical pattern of male philopatry and female-biased dispersal. Allelic richness averaged over seven loci was not significantly different between an African and an Arabian population, suggesting that loss of variation during the colonisation of Arabia does not explain low Y variation. Finally, in the absence of nucleotide polymorphism, it is unclear to what extent selection could be responsible for low Y variation in this species.

High levels of nucleotide diversity in the European rabbit (Oryctolagus cuniculus) SRY gene

We have sequenced 2,388 bp of the European rabbit sex determining region Y (SRY) gene. These data provide a 10-fold increase in the coverage of the Y chromosome in this species, including the entire open reading frame of the SRY, the polyadenylation signal, and two repetitive sequences in the 5' -region. A survey of 2021 bp of this gene in eight domestic breeds and four wild individuals revealed a total of nine single nucleotide polymorphisms and one indel, defining two deeply divergent lineages. The resulting estimation of nucleotide diversity (pi=1.34 x10(-3)) is very high when compared with other species, but no variability was detected among the domestic breeds. This study represents a first step in the characterization of the European rabbit Y chromosome and its variability. These sequences can be used in additional phylogeographical analyses of the European rabbit and other Leporid species, as well as in evolutionary studies of sex determination and the Y chromosome in wild species.

Analysis of sex-linked sequences supports a new mammal species in Europe

European mammals have been the focus of particularly detailed taxonomic studies by traditional morphological methods. However, DNA analyses have the potential to reveal additional, cryptic species. We describe two highly divergent evolutionary lineages within a small Eurasian mammal, the field vole (Microtus agrestis). We show that the two lineages can be detected not only with maternally (mitochondrial DNA), but also with paternally (Y chromosome) and biparentally (X chromosome) inherited DNA sequences. Reciprocal monophyly of all genealogies and their congruent geographical distributions is consistent with reproductive isolation. Our results suggest that the field vole should be reclassified as two separate species.

The expansion of conservation genetics

The 'crisis discipline' of conservation biology has voraciously incorporated many technologies to speed up and increase the accuracy of conservation decision-making. Genetic approaches to characterizing endangered species or areas that contain endangered species are prime examples of this. Technical advances in areas such as high-throughput sequencing, microsatellite analysis and non-invasive DNA sampling have led to a much-expanded role for genetics in conservation. Such expansion will allow for more precise conservation decisions to be made and, more importantly, will allow conservation genetics to contribute to area- and landscape-based decision-making processes.