Inferring the geographic mode of speciation by contrasting autosomal and sex-linked genetic diversity

Cryptic species in a colorful genus: Integrative taxonomy of the bush robins (Aves, Muscicapidae, Tarsiger) suggests two overlooked species

Several cryptic avian species have been validated by recent integrative taxonomic efforts in the Sino-Himalayan mountains, indicating that avian diversity in this global biodiversity hotspot may be underestimated. In the present study, we investigated species limits in the genus Tarsiger, the bush robins, a group of montane forest specialists with high species richness in the Sino-Himalayan region. Based on comprehensive sampling of all 11 subspecies of the six currently recognized species, we applied an integrative taxonomic approach by combining multilocus, acoustic, plumage and morphometric analyses. Our results reveal that the isolated north-central Chinese populations of Tarsiger cyanurus, described as the subspecies albocoeruleus but usually considered invalid, is distinctive in genetics and vocalisation, but only marginally differentiated in morphology. We also found the Taiwan endemic T. indicus formosanus to be distinctive in genetics, song and morphology from T. i. indicus and T. i. yunnanensis of the Sino-Himalayan mountains. Moreover, Bayesian species delimitation using BPP suggests that both albocoeruleus and formosanus merit full species status. We propose their treatment as 'Qilian Bluetail' T. albocoeruleus and 'Taiwan Bush Robin' T. formosanus, respectively.

The rates of introgression and barriers to genetic exchange between hybridizing species: sex chromosomes vs autosomes

Interspecific crossing experiments have shown that sex chromosomes play a major role in reproductive isolation between many pairs of species. However, their ability to act as reproductive barriers, which hamper interspecific genetic exchange, has rarely been evaluated quantitatively compared to Autosomes. This genome-wide limitation of gene flow is essential for understanding the complete separation of species, and thus speciation. Here, we develop a mainland-island model of secondary contact between hybridizing species of an XY (or ZW) sexual system. We obtain theoretical predictions for the frequency of introgressed alleles, and the strength of the barrier to neutral gene flow for the two types of chromosomes carrying multiple interspecific barrier loci. Theoretical predictions are obtained for scenarios where introgressed alleles are rare. We show that the same analytical expressions apply for sex chromosomes and autosomes, but with different sex-averaged effective parameters. The specific features of sex chromosomes (hemizygosity and absence of recombination in the heterogametic sex) lead to reduced levels of introgression on the X (or Z) compared to autosomes. This effect can be enhanced by certain types of sex-biased forces, but it remains overall small (except when alleles causing incompatibilities are recessive). We discuss these predictions in the light of empirical data comprising model-based tests of introgression and cline surveys in various biological systems.

Strict allopatric speciation of sky island Pyrrhula erythaca species complex

Increasing evidence of post-divergence gene flow between taxa is shifting our understanding on the mode of speciation. A fundamental question arises concerning the circumstances under which strict allopatric speciation still holds true. Sky island populations might undergo reduced gene flow by niche conservatism to highland habitats and follow divergence in an allopatric manner. In this study, we tested this hypothesis in the sky island Grey-headed Bullfinch (Pyrrhula erythaca) species complex via statistical analyses of both genetic and ecological data. Results of coalescent-based analysis of multiple nuclear loci suggested that P. e. owstoni likely colonized Taiwan island during the severe mid-Pleistocene glacial climate followed by strictly allopatric divergence from P. e. erythaca distributed in Himalayas-Hengduan mountains and central North China. Results of ecological niche modeling suggested that their speciation may be attributed to the niche conservatism of these birds and the lack of a suitable ecological corridor during subsequent milder glacial episodes. In addition, we delimited the traditionally defined P. erythaca into two full species, P. erythaca in the Asian mainland and P. owstoni on the island of Taiwan, based on both genetic and behavioural evidences. These results suggest that ecology can have a dynamic role in allowing highland populations to expand their ranges and isolated by habitat barriers to diversify in a strictly allopatric manner.

Historical demography of four gecko species specializing in boulder cave habitat: Implications in the evolutionary dead end hypothesis and conservation

Specialization in narrow ecological niches may not only help species to survive in competitive or unique environments but also contribute to their extermination over evolutionary time. Although the "evolutionary dead end" hypothesis has long been debated, empirical evidence from species with detailed information on niche specialization and evolutionary history remains rare. Here we use a group of four closely related Cnemaspis gecko species that depend highly on granite boulder caves in the Mekong Delta to investigate the potential impact of ecological specialization on their evolution and population dynamics. Isolated by unsuitable floodplain habitats, these boulder-dwelling geckos are among the most narrowly distributed Squamata in the world. We applied several coalescence-based approaches combined with the RAD-seq technique to estimate their divergence times, gene flow and demographic fluctuations during the speciation and population differentiation processes. Our results reveal long-term population shrinkage in the four geckos and limited gene flow during their divergence. The results suggest that the erosion and fragmentation of the granite boulder hills have greatly impacted population divergence and declines. The habitat specialization of these geckos has led to fine-scaled speciation in these granite rocky hills; in contrast, specialization might also have pushed these species toward the edge of extinction. Our study also emphasizes the conservation urgency of these vulnerable, cave-dependent geckos.

The first set of universal nuclear protein-coding loci markers for avian phylogenetic and population genetic studies

Multiple nuclear markers provide genetic polymorphism data for molecular systematics and population genetic studies. They are especially required for the coalescent-based analyses that can be used to accurately estimate species trees and infer population demographic histories. However, in avian evolutionary studies, these powerful coalescent-based methods are hindered by the lack of a sufficient number of markers. In this study, we designed PCR primers to amplify 136 nuclear protein-coding loci (NPCLs) by scanning the published Red Junglefowl (Gallus gallus) and Zebra Finch (Taeniopygia guttata) genomes. To test their utility, we amplified these loci in 41 bird species representing 23 Aves orders. The sixty-three best-performing NPCLs, based on high PCR success rates, were selected which had various mutation rates and were evenly distributed across 17 avian autosomal chromosomes and the Z chromosome. To test phylogenetic resolving power of these markers, we conducted a Neoavian phylogenies analysis using 63 concatenated NPCL markers derived from 48 whole genomes of birds. The resulting phylogenetic topology, to a large extent, is congruence with results resolved by previous whole genome data. To test the level of intraspecific polymorphism in these makers, we examined the genetic diversity in four populations of the Kentish Plover (Charadrius alexandrinus) at 17 of NPCL markers chosen at random. Our results showed that these NPCL markers exhibited a level of polymorphism comparable with mitochondrial loci. Therefore, this set of pan-avian nuclear protein-coding loci has great potential to facilitate studies in avian phylogenetics and population genetics.

Genetic Diversity on the Sex Chromosomes

Levels and patterns of genetic diversity can provide insights into a population’s history. In species with sex chromosomes, differences between genomic regions with unique inheritance patterns can be used to distinguish between different sets of possible demographic and selective events. This review introduces the differences in population history for sex chromosomes and autosomes, provides the expectations for genetic diversity across the genome under different evolutionary scenarios, and gives an introductory description for how deviations in these expectations are calculated and can be interpreted. Predominantly, diversity on the sex chromosomes has been used to explore and address three research areas: 1) Mating patterns and sex-biased variance in reproductive success, 2) signatures of selection, and 3) evidence for modes of speciation and introgression. After introducing the theory, this review catalogs recent studies of genetic diversity on the sex chromosomes across species within the major research areas that sex chromosomes are typically applied to, arguing that there are broad similarities not only between male-heterogametic (XX/XY) and female-heterogametic (ZZ/ZW) sex determination systems but also any mating system with reduced recombination in a sex-determining region. Further, general patterns of reduced diversity in nonrecombining regions are shared across plants and animals. There are unique patterns across populations with vastly different patterns of mating and speciation, but these do not tend to cluster by taxa or sex determination system.

Porcine Y-chromosome variation is consistent with the occurrence of paternal gene flow from non-Asian to Asian populations

Pigs (Sus scrofa) originated in Southeast Asia and expanded to Europe and North Africa approximately 1 MYA. Analyses of porcine Y-chromosome variation have shown the existence of two main haplogroups that are highly divergent, a result that is consistent with previous mitochondrial and autosomal data showing that the Asian and non-Asian pig populations remained geographically isolated until recently. Paradoxically, one of these Y-chromosome haplogroups is extensively shared by pigs and wild boars from Asia and Europe, an observation that is difficult to reconcile with a scenario of prolonged geographic isolation. To shed light on this issue, we genotyped 33 Y-linked SNPs and one indel in a worldwide sample of pigs and wild boars and sequenced a total of 9903 nucleotide sites from seven loci distributed along the Y-chromosome. Notably, the nucleotide diversity per site at the Y-linked loci (0.0015 in Asian pigs) displayed the same order of magnitude as that described for autosomal loci (~0.0023), a finding compatible with a process of sustained and intense isolation. We performed an approximate Bayesian computation analysis focused on the paternal diversity of wild boars and local pig breeds in which we compared three demographic models: two isolation models (I models) differing in the time of isolation and a model of isolation with recent unidirectional migration (IM model). Our results suggest that the most likely explanation for the extensive sharing of one Y-chromosome haplogroup between non-Asian and Asian populations is a recent and unidirectional (non-Asian > Asian) paternal migration event.

Sequence diversity patterns suggesting balancing selection in partially sex-linked genes of the plant Silene latifolia are not generated by demographic history or gene flow

DNA sequence diversity in genes in the partially sex-linked pseudoautosomal region (PAR) of the sex chromosomes of the plant Silene latifolia is higher than expected from within-species diversity of other genes. This could be the footprint of sexually antagonistic (SA) alleles that are maintained by balancing selection in a PAR gene (or genes) and affect polymorphism in linked genome regions. SA selection is predicted to occur during sex chromosome evolution, but it is important to test whether the unexpectedly high sequence polymorphism could be explained without it, purely by the combined effects of partial linkage with the sex-determining region and the population's demographic history, including possible introgression from Silene dioica. To test this, we applied approximate Bayesian computation-based model choice to autosomal sequence diversity data, to find the most plausible scenario for the recent history of S. latifolia and then to estimate the posterior density of the most relevant parameters. We then used these densities to simulate variation to be expected at PAR genes. We conclude that an excess of variants at high frequencies at PAR genes should arise in S. latifolia populations only for genes with strong associations with fully sex-linked genes, which requires closer linkage with the fully sex-linked region than that estimated for the PAR genes where apparent deviations from neutrality were observed. These results support the need to invoke selection to explain the S. latifolia PAR gene diversity, and encourage further work to test the possibility of balancing selection due to sexual antagonism.

Likelihood-Free Inference in High-Dimensional Models

Methods that bypass analytical evaluations of the likelihood function have become an indispensable tool for statistical inference in many fields of science. These so-called likelihood-free methods rely on accepting and rejecting simulations based on summary statistics, which limits them to low-dimensional models for which the value of the likelihood is large enough to result in manageable acceptance rates. To get around these issues, we introduce a novel, likelihood-free Markov chain Monte Carlo (MCMC) method combining two key innovations: updating only one parameter per iteration and accepting or rejecting this update based on subsets of statistics approximately sufficient for this parameter. This increases acceptance rates dramatically, rendering this approach suitable even for models of very high dimensionality. We further derive that for linear models, a one-dimensional combination of statistics per parameter is sufficient and can be found empirically with simulations. Finally, we demonstrate that our method readily scales to models of very high dimensionality, using toy models as well as by jointly inferring the effective population size, the distribution of fitness effects (DFE) of segregating mutations, and selection coefficients for each locus from data of a recent experiment on the evolution of drug resistance in influenza.

Incipient speciation with gene flow on a continental island: Species delimitation of the Hainan Hwamei (Leucodioptron canorum owstoni, Passeriformes, Aves)

Because of their isolation, continental islands (e.g., Madagascar) are often thought of as ideal systems to study allopatric speciation. However, many such islands have been connected intermittently to their neighboring continent during recent periods of glaciation, which may cause frequent contact between the diverging populations on the island and continent. As a result, the speciation processes on continental islands may not meet the prerequisites for strictly allopatric speciation. We used multiple lines of evidence to re-evaluate the taxonomic status of the Hainan Hwamei (Leucodioptron canorum owstoni), which is endemic to Hainan, the largest continental island in the South China Sea. Our analysis of mitochondrial DNA and twelve nuclear loci suggests that the Hainan Hwamei can be regarded as an independent species (L. owstoni); the morphological traits of the Hainan Hwamei also showed significant divergence from those of their mainland sister taxon, the Chinese Hwamei (L. canorum). We also inferred the divergence history of the Hainan and Chinese Hwamei to see whether their divergence was consistent with a strictly allopatric model. Our results suggest that the two Hwameis split only 0.2 million years ago with limited asymmetrical post-divergence gene flow. This implies that the Hainan Hwamei is an incipient species and that speciation occurred through ecologically divergent selection and/or assortative mating rather than a strictly allopatric process.

Sex-linked genomic variation and its relationship to avian plumage dichromatism and sexual selection

BACKGROUND

Sexual dichromatism is the tendency for sexes to differ in color pattern and represents a striking form of within-species morphological variation. Conspicuous intersexual differences in avian plumage are generally thought to result from Darwinian sexual selection, to the extent that dichromatism is often treated as a surrogate for the intensity of sexual selection in phylogenetic comparative studies. Intense sexual selection is predicted to leave a footprint on genetic evolution by reducing the relative genetic diversity on sex chromosome to that on the autosomes.

RESULTS

In this study, we test the association between plumage dichromatism and sex-linked genetic diversity using eight species pairs with contrasting levels of dichromatism. We estimated Z-linked and autosomal genetic diversity for these non-model avian species using restriction-site associated (RAD) loci that covered ~3 % of the genome. We find that monochromatic birds consistently have reduced sex-linked genomic variation relative to phylogenetically-paired dichromatic species and this pattern is robust to mutational biases.

CONCLUSIONS

Our results are consistent with several interpretations. If present-day sexual selection is stronger in dichromatic birds, our results suggest that its impact on sex-linked genomic variation is offset by other processes that lead to proportionately lower Z-linked variation in monochromatic species. We discuss possible factors that may contribute to this discrepancy between phenotypes and genomic variation. Conversely, it is possible that present-day sexual selection -- as measured by the variance in male reproductive success -- is stronger in the set of monochromatic taxa we have examined, potentially reflecting the importance of song, behavior and other non-plumage associated traits as targets of sexual selection. This counterintuitive finding suggests that the relationship between genomic variation and sexual selection is complex and highlights the need for a more comprehensive survey of genomic variation in avian taxa that vary markedly in social and genetic mating systems.

Drosophila suzukii: the genetic footprint of a recent, worldwide invasion

Native to Asia, the soft-skinned fruit pest Drosophila suzukii has recently invaded the United States and Europe. The eastern United States represents the most recent expansion of their range, and presents an opportunity to test alternative models of colonization history. Here, we investigate the genetic population structure of this invasive fruit fly, with a focus on the eastern United States. We sequenced six X-linked gene fragments from 246 individuals collected from a total of 12 populations. We examine patterns of genetic diversity within and between populations and explore alternative colonization scenarios using approximate Bayesian computation. Our results indicate high levels of nucleotide diversity in this species and suggest that the recent invasions of Europe and the continental United States are independent demographic events. More broadly speaking, our results highlight the importance of integrating population structure into demographic models, particularly when attempting to reconstruct invasion histories. Finally, our simulation results illustrate the general challenge in reconstructing invasion histories using genetic data and suggest that genome-level data are often required to distinguish among alternative demographic scenarios.