msABC: a modification of Hudson’s ms to facilitate multi‐locus ABC analysis

Geogenomic Predictors of Genetree Heterogeneity Explain Phylogeographic and Introgression History: A Case Study in an Amazonian Bird (Thamnophilus aethiops)

Can knowledge about genome architecture inform biogeographic and phylogenetic inference? Selection, drift, recombination, and gene flow interact to produce a genomic landscape of divergence wherein patterns of differentiation and genealogy vary nonrandomly across the genomes of diverging populations. For instance, genealogical patterns that arise due to gene flow should be more likely to occur on smaller chromosomes, which experience high recombination, whereas those tracking histories of geographic isolation (reduced gene flow caused by a barrier) and divergence should be more likely to occur on larger and sex chromosomes. In Amazonia, populations of many bird species diverge and introgress across rivers, resulting in reticulated genomic signals. Herein, we used reduced representation genomic data to disentangle the evolutionary history of 4 populations of an Amazonian antbird, Thamnophilus aethiops, whose biogeographic history was associated with the dynamic evolution of the Madeira River Basin. Specifically, we evaluate whether a large river capture event ca. 200 Ka, gave rise to reticulated genealogies in the genome by making spatially explicit predictions about isolation and gene flow based on knowledge about genomic processes. We first estimated chromosome-level phylogenies and recovered 2 primary topologies across the genome. The first topology (T1) was most consistent with predictions about population divergence and was recovered for the Z-chromosome. The second (T2), was consistent with predictions about gene flow upon secondary contact. To evaluate support for these topologies, we trained a convolutional neural network to classify our data into alternative diversification models and estimate demographic parameters. The best-fit model was concordant with T1 and included gene flow between non-sister taxa. Finally, we modeled levels of divergence and introgression as functions of chromosome length and found that smaller chromosomes experienced higher gene flow. Given that (1) genetrees supporting T2 were more likely to occur on smaller chromosomes and (2) we found lower levels of introgression on larger chromosomes (and especially the Z-chromosome), we argue that T1 represents the history of population divergence across rivers and T2 the history of secondary contact due to barrier loss. Our results suggest that a significant portion of genomic heterogeneity arises due to extrinsic biogeographic processes such as river capture interacting with intrinsic processes associated with genome architecture. Future phylogeographic studies would benefit from accounting for genomic processes, as different parts of the genome reveal contrasting, albeit complementary histories, all of which are relevant for disentangling the intricate geogenomic mechanisms of biotic diversification. [Amazonia; biogeography; demographic modeling; gene flow; gene tree; genome architecture; geogenomics; introgression; linked selection; neural network; phylogenomic; phylogeography; reproductive isolation; speciation; species tree.].

Integrative methods reveal multiple drivers of diversification in rice paddy snakes

Divergence dating analyses in systematics provide a framework to develop and test biogeographic hypotheses regarding speciation. However, as molecular datasets grow from multilocus to genomic, sample sizes decrease due to computational burdens, and the testing of fine-scale biogeographic hypotheses becomes difficult. In this study, we use coalescent demographic models to investigate the diversification of poorly known rice paddy snakes from Southeast Asia (Homalopsidae: Hypsiscopus), which have conflicting dates of origin based on previous studies. We use coalescent modeling to test the hypothesis that Hypsiscopus diversified 2.5 mya during the Khorat Plateau uplift in Thailand. Additionally, we use ecological niche analyses to identify potential differences in the niche space of the two most widely distributed species in the past and present. Our results suggest Hypsiscopus diversified ~ 2.4 mya, supporting that the Khorat Plateau may have initiated the diversification of rice paddy snakes. We also find significant niche differentiation and shifts between species of Hypsiscopus, indicating that environmental differences may have sustained differentiation of this genus after the Khorat Plateau uplift. Our study expands on the diversification history of snakes in Southeast Asia, and highlights how results from smaller multilocus datasets can be useful in developing and testing biogeographic hypotheses alongside genomic datasets.

Genomic Architecture Predicts Tree Topology, Population Structuring, and Demographic History in Amazonian Birds

Geographic barriers are frequently invoked to explain genetic structuring across the landscape. However, inferences on the spatial and temporal origins of population variation have been largely limited to evolutionary neutral models, ignoring the potential role of natural selection and intrinsic genomic processes known as genomic architecture in producing heterogeneity in differentiation across the genome. To test how variation in genomic characteristics (e.g. recombination rate) impacts our ability to reconstruct general patterns of differentiation between species that cooccur across geographic barriers, we sequenced the whole genomes of multiple bird populations that are distributed across rivers in southeastern Amazonia. We found that phylogenetic relationships within species and demographic parameters varied across the genome in predictable ways. Genetic diversity was positively associated with recombination rate and negatively associated with species tree support. Gene flow was less pervasive in genomic regions of low recombination, making these windows more likely to retain patterns of population structuring that matched the species tree. We further found that approximately a third of the genome showed evidence of selective sweeps and linked selection, skewing genome-wide estimates of effective population sizes and gene flow between populations toward lower values. In sum, we showed that the effects of intrinsic genomic characteristics and selection can be disentangled from neutral processes to elucidate spatial patterns of population differentiation.

Inferring the Demographic History and Inheritance Mode of Tetraploid Species Using ABC

Genomic patterns of diversity and divergence are impacted by certain life history traits, reproductive systems, and demographic history. The latter is characterized by fluctuations in population sizes over time, as well as by temporal patterns of introgression. For a given organism, identifying a demographic history that deviates from the standard neutral model allows a better understanding of its evolution but also helps to reduce the risk of false positives when screening for molecular targets of natural selection. Tetraploid organisms and beyond have demographic histories that are complicated by the mode of polyploidization, the mode of inheritance, and different scenarios of gene flow between sub-genomes and diploid parental species. Here we provide guidelines for experimenters wishing to address these issues through a flexible statistical framework: approximate Bayesian computation (ABC). The emphasis is on the general philosophy of the approach to encourage future users to exploit the enormous flexibility of ABC beyond the limitations imposed by generalist data analysis pipelines.

Late Pleistocene landscape changes and habitat specialization as promoters of population genomic divergence in Amazonian floodplain birds

Although vicariant processes are expected to leave similar genomic signatures among codistributed taxa, ecological traits such as habitat and stratum can influence genetic divergence within species. Here, we combined landscape history and habitat specialization to understand the historical and ecological factors responsible for current levels of genetic divergence in three species of birds specialized in seasonally flooded habitats in muddy rivers and which are widespread in the Amazon basin but have isolated populations in the Rio Branco. Populations of the white-bellied spinetail (Mazaria propinqua), lesser wagtail-tyrant (Stigmatura napensis) and bicolored conebill (Conirostrum bicolor) are currently isolated in the Rio Branco by the black-waters of the lower Rio Negro, offering a unique opportunity to test the effect of river colour as a barrier to gene flow. We used ultraconserved elements (UCEs) to test alternative hypotheses of population history in a comparative phylogeographical approach by modelling genetic structure, demographic history and testing for shared divergence time among codistributed taxa. Our analyses revealed that (i) all three populations from the Rio Branco floodplains are genetically distinct from other populations along the Amazon River floodplains; (ii) these divergences are the result of at least two distinct events, consistent with species habitat specialization; and (iii) the most likely model of population evolution includes lower population connectivity during the Late Pleistocene transition (~250,000 years ago), with gene flow being completely disrupted after the Last Glacial Maximum (~21,000 years ago). Our findings highlight how landscape evolution modulates population connectivity in habitat specialist species and how organisms can have different responses to the same historical processes of environmental change, depending on their habitat affinity.

Coalescent simulations indicate that the São Francisco River is a biogeographic barrier for six vertebrates in a seasonally dry South American forest

The riverine barrier hypothesis has been extensively explored in Neotropical rainforests, while its importance in drier regions such as the Caatinga, a seasonally dry tropical forest in northeastern Brazil, has only recently received more attention. The Caatinga is bisected by the São Francisco River (SFR), which has long been suggested to be an important biogeographic feature in the region. However, recent studies have found mixed support for the role of the SFR as a hard barrier, most of them relying on the presence or absence of genetic breaks congruent with its course. Here, we used published multilocus and next-generation data from six vertebrate species to test the SFR’s strength as a barrier. Using model-based approaches (approximate Bayesian computation and supervised machine learning), we tested demographic models incorporating full, intermediate, and zero migration across the SFR, estimating divergence times and migration rates for each species. We found support for the SFR’s role as a barrier, allowing gene flow for some species. Estimated divergence times varied among species but are limited to the late Pleistocene, coherent with one of several proposed paleocourse changes in the river’s geological history. Contrary to the mixed results of previous studies, our study supports the SFR as an important phylogeographic barrier across different taxonomic groups, driving diversification in the Caatinga.

Efficient ancestry and mutation simulation with msprime 1.0

Stochastic simulation is a key tool in population genetics, since the models involved are often analytically intractable and simulation is usually the only way of obtaining ground-truth data to evaluate inferences. Because of this, a large number of specialized simulation programs have been developed, each filling a particular niche, but with largely overlapping functionality and a substantial duplication of effort. Here, we introduce msprime version 1.0, which efficiently implements ancestry and mutation simulations based on the succinct tree sequence data structure and the tskit library. We summarize msprime’s many features, and show that its performance is excellent, often many times faster and more memory efficient than specialized alternatives. These high-performance features have been thoroughly tested and validated, and built using a collaborative, open source development model, which reduces duplication of effort and promotes software quality via community engagement.

Microevolutionary dynamics show tropical valleys are deeper for montane birds of the Atlantic Forest

Tropical mountains hold more biodiversity than their temperate counterparts, and this disparity is often associated with the latitudinal climatic gradient. However, distinguishing the impact of latitude versus the background effects of species history and traits is challenging due to the evolutionary distance between tropical and temperate assemblages. Here, we test whether microevolutionary processes are linked to environmental variation across a sharp latitudinal transition in 21 montane birds of the southern Atlantic Forest in Brazil. We find that effective dispersal within populations in the tropical mountains is lower and genomic differentiation is better predicted by the current environmental complexity of the region than within the subtropical populations. The concordant response of multiple co-occurring populations is consistent with spatial climatic variability as a major process driving population differentiation. Our results provide evidence for how a narrow latitudinal gradient can shape microevolutionary processes and contribute to broader scale biodiversity patterns.

There are many hypotheses for why the tropics are more biodiverse than higher latitudes. Phylogenomic analyses of 21 montane birds finds that tropical birds disperse less and have more genetically structured populations than their counterparts at higher latitudes, possibly due to a larger elevational climate gradient in the tropics

Rugged relief and climate promote isolation and divergence between two neotropical cold-associated birds

The role of historical factors in establishing patterns of diversity in tropical mountains is of interest to understand the buildup of megadiverse biotas. In these regions, the historical processes of range fragmentation and contraction followed by dispersal are thought to be mediated by the interplay between rugged relief (complex topography) and climate fluctuations and likely explain most of the dynamics of diversification in plants and animals. Although empirical studies addressing the interaction between climate and topography have provided invaluable insights into population divergence and speciation patterns in tropical montane organisms, a more detailed and robust test of such processes in an explicit spatio-temporal framework is still lacking. Consequently, our ability to gain insights into historical range shifts over time and the genomic footprint left by them is limited. Here, we used niche modeling and subgenomic population-level datasets to explore the evolution of two species of warbling finches (genus Microspingus) disjunctly distributed across the Montane Atlantic Forest, a Neotropical region with complex geological and environmental histories. Population structure inferences suggest a scenario of three genetically differentiated populations, which are congruent with both geography and phenotypic variation. Demographic simulations support asynchronous isolation of these populations as recently as ∼40,000 years ago, relatively stable population sizes over recent time, and past gene flow subsequent to divergence. Throughout the last 800,000 years, niche models predicted extensive expansion into lowland areas with increasing overlap of species distributions during glacial periods, with prominent retractions and isolation into higher altitudes during interglacials, which are in line with signs of introgression of currently isolated populations. These results support a dual role of cyclical climatic changes: population divergence and persistence in mountain tops during warm periods followed by periods of expansion and admixture in lower elevations during cold periods. Our results underscore the role of the interplay between landscape and climate as an important mechanism in the evolution of the Neotropical montane biota.

Molecular population genetics of Sex-lethal (Sxl) in the Drosophila melanogaster species group: a locus that genetically interacts with Wolbachia pipientis in Drosophila melanogaster

Sex-lethal (Sxl) is the sex determination switch in Drosophila, and also plays a critical role in germ-line stem cell daughter differentiation in Drosophila melanogaster. Three female-sterile alleles at Sxl in D. melanogaster were previously shown to genetically interact to varying degrees with the maternally inherited endosymbiont Wolbachia pipientis. Given this genetic interaction and W. pipientis’ ability to manipulate reproduction in Drosophila, we carried out a careful study of both the population genetics (within four Drosophila species) and molecular evolutionary analysis (across 20 Drosophila species) of Sxl. Consistent with earlier studies, we find that selective constraint has played a prominent role in Sxl’s molecular evolution within Drosophila, but we also observe patterns that suggest both episodic bursts of protein evolution and recent positive selection at Sxl. The episodic nature of Sxl’s protein evolution is discussed in light of its genetic interaction with W. pipientis.

An overview of current population genomics methods for the analysis of whole-genome resequencing data in eukaryotes

Characterizing the population history of a species and identifying loci underlying local adaptation is crucial in functional ecology, evolutionary biology, conservation and agronomy. The constant improvement of high-throughput sequencing techniques has facilitated the production of whole genome data in a wide range of species. Population genomics now provides tools to better integrate selection into a historical framework, and take into account selection when reconstructing demographic history. However, this improvement has come with a profusion of analytical tools that can confuse and discourage users. Such confusion limits the amount of information effectively retrieved from complex genomic data sets, and impairs the diffusion of the most recent analytical tools into fields such as conservation biology. It may also lead to redundancy among methods. To address these isssues, we propose an overview of more than 100 state-of-the-art methods that can deal with whole genome data. We summarize the strategies they use to infer demographic history and selection, and discuss some of their limitations. A website listing these methods is available at www.methodspopgen.com.

Phylogeographic model selection using convolutional neural networks

The discipline of phylogeography has evolved rapidly in terms of the analytical toolkit used to analyse large genomic data sets. Despite substantial advances, analytical tools that could potentially address the challenges posed by increased model complexity have not been fully explored. For example, deep learning techniques are underutilized for phylogeographic model selection. In non-model organisms, the lack of information about their ecology and evolution can lead to uncertainty about which demographic models are appropriate. Here, we assess the utility of convolutional neural networks (CNNs) for assessing demographic models in South American lizards in the genus Norops. Three demographic scenarios (constant, expansion, and bottleneck) were considered for each of four inferred population-level lineages, and we found that the overall model accuracy was higher than 98% for all lineages. We then evaluated a set of 26 models that accounted for evolutionary relationships, gene flow, and changes in effective population size among the four lineages, identifying a single model with an estimated overall accuracy of 87% when using CNNs. The inferred demography of the lizard system suggests that gene flow between non-sister populations and changes in effective population sizes through time, probably in response to Pleistocene climatic oscillations, have shaped genetic diversity in this system. Approximate Bayesian computation (ABC) was applied to provide a comparison to the performance of CNNs. ABC was unable to identify a single model among the larger set of 26 models in the subsequent analysis. Our results demonstrate that CNNs can be easily and usefully incorporated into the phylogeographer's toolkit.

Population Genomics Insights into the First Wave of COVID-19

Full-genome-sequence computational analyses of the SARS-coronavirus (CoV)-2 genomes allow us to understand the evolutionary events and adaptability mechanisms. We used population genetics analyses on human SARS-CoV-2 genomes available on 2 April 2020 to infer the mutation rate and plausible recombination events between the Betacoronavirus genomes in nonhuman hosts that may have contributed to the evolution of SARS-CoV-2. Furthermore, we localized the targets of recent and strong, positive selection during the first pandemic wave. The genomic regions that appear to be under positive selection are largely co-localized with regions in which recombination from nonhuman hosts took place. Our results suggest that the pangolin coronavirus genome may have contributed to the SARS-CoV-2 genome by recombination with the bat coronavirus genome. However, we find evidence for additional recombination events that involve coronavirus genomes from other hosts, i.e., hedgehogs and sparrows. We further infer that recombination may have recently occurred within human hosts. Finally, we estimate the parameters of a demographic scenario involving an exponential growth of the size of the SARS-CoV-2 populations that have infected European, Asian, and Northern American cohorts, and we demonstrate that a rapid exponential growth in population size from the first wave can support the observed polymorphism patterns in SARS-CoV-2 genomes.

Phylogenomic Data Reveal Widespread Introgression Across the Range of an Alpine and Arctic Specialist

Understanding how gene flow affects population divergence and speciation remains challenging. Differentiating one evolutionary process from another can be difficult because multiple processes can produce similar patterns, and more than one process can occur simultaneously. Although simple population models produce predictable results, how these processes balance in taxa with patchy distributions and complicated natural histories is less certain. These types of populations might be highly connected through migration (gene flow), but can experience stronger effects of genetic drift and inbreeding, or localized selection. Although different signals can be difficult to separate, the application of high-throughput sequence data can provide the resolution necessary to distinguish many of these processes. We present whole-genome sequence data for an avian species group with an alpine and arctic tundra distribution to examine the role that different population genetic processes have played in their evolutionary history. Rosy-finches inhabit high elevation mountaintop sky islands and high-latitude island and continental tundra. They exhibit extensive plumage variation coupled with low levels of genetic variation. Additionally, the number of species within the complex is debated, making them excellent for studying the forces involved in the process of diversification, as well as an important species group in which to investigate species boundaries. Total genomic variation suggests a broadly continuous pattern of allele frequency changes across the mainland taxa of this group in North America. However, phylogenomic analyses recover multiple distinct, well supported, groups that coincide with previously described morphological variation and current species-level taxonomy. Tests of introgression using D-statistics and approximate Bayesian computation reveal significant levels of introgression between multiple North American taxa. These results provide insight into the balance between divergent and homogenizing population genetic processes and highlight remaining challenges in interpreting conflict between different types of analytical approaches with whole-genome sequence data. [ABBA-BABA; approximate Bayesian computation; gene flow; phylogenomics; speciation; whole-genome sequencing.].

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.

SGID: a comprehensive and interactive database of the silkworm

Although the domestic silkworm (Bombyx mori) is an important model and economic animal, there is a lack of comprehensive database for this organism. Here, we developed the silkworm genome informatics database (SGID). It aims to bring together all silkworm-related biological data and provide an interactive platform for gene inquiry and analysis. The function annotation in SGID is thorough and covers 98% of the silkworm genes. The annotation details include function description, Gene Ontology, Kyoto Encyclopedia of Genes and Genomes pathway, subcellular location, transmembrane topology, protein secondary/tertiary structure, homologous group and transcription factor. SGID provides genome-scale visualization of population genetics test results based on high-depth resequencing data of 158 silkworm samples. It also provides interactive analysis tools of transcriptomic and epigenomic data from 79 NCBI BioProjects. SGID will be extremely useful to silkworm research in the future.

Climatic dynamics and topography control genetic variation in Atlantic Forest montane birds

Montane organisms responded to Quaternary climate change by tracking suitable habitat along elevational gradients. However, it is unclear whether these past climatic dynamics generated predictable patterns of genetic diversity in co-occurring montane taxa. To test if the genetic variation is associated with historical changes in the elevational distribution of montane habitats, we integrated paleoclimatic data and a model selection approach for testing the demographic history of five co-distributed bird species occurring in the southern Atlantic Forest sky islands. We found that changes in historical population sizes and current genetic diversity are attributable to habitat dynamics among time periods and the current elevational distribution of populations. Taxa with populations restricted to the more climatically dynamic southern mountain block (SMB) had, on average, a six-fold demographic expansion, whereas the populations from the northern mountain block (NMB) remained constant. In the current configuration of the southern Atlantic Forest montane habitats, populations in the SMB have more widespread elevational distributions, occur at lower elevations, and harbor higher levels of genetic diversity than NMB populations. Despite the apparent coupling of demographic and climatic oscillations, our data rejected simultaneous population structuring due to historical habitat fragmentation. Demographic modeling indicated that the species had different modes of differentiation, and varied in the timing of divergence and the degree of gene flow across mountain blocks. Our results suggest that the heterogeneous distribution of genetic variation in birds of the Atlantic Forest sky islands is associated with the interplay between topography and climate of distinct mountains, leading to predictable patterns of genetic diversity.

Independent domestication events in the blue-cheese fungus Penicillium roqueforti

Domestication provides an excellent framework for studying adaptive divergence. Using population genomics and phenotypic assays, we reconstructed the domestication history of the blue cheese mould Penicillium roqueforti. We showed that this fungus was domesticated twice independently. The population used in Roquefort originated from an old domestication event associated with weak bottlenecks and exhibited traits beneficial for pre‐industrial cheese production (slower growth in cheese and greater spore production on bread, the traditional multiplication medium). The other cheese population originated more recently from the selection of a single clonal lineage, was associated with all types of blue cheese worldwide except Roquefort, and displayed phenotypes more suited for industrial cheese production (high lipolytic activity, efficient cheese cavity colonization ability and salt tolerance). We detected genomic regions affected by recent positive selection and putative horizontal gene transfers. This study sheds light on the processes of rapid adaptation and raises questions about genetic resource conservation.

see also the Perspective by Brigida Gallone, Jan Steensels and Kevin J. Verstrepen.

Demographic Histories and Genome-Wide Patterns of Divergence in Incipient Species of Shorebirds

Understanding how incipient species are maintained with gene flow is a fundamental question in evolutionary biology. Whole genome sequencing of multiple individuals holds great potential to illustrate patterns of genomic differentiation as well as the associated evolutionary histories. Kentish (Charadrius alexandrinus) and the white-faced (C. dealbatus) plovers, which differ in their phenotype, ecology and behavior, are two incipient species and parapatrically distributed in East Asia. Previous studies show evidence of genetic diversification with gene flow between the two plovers. Under this scenario, it is of great importance to explore the patterns of divergence at the genomic level and to determine whether specific regions are involved in reproductive isolation and local adaptation. Here we present the first population genomic analysis of the two incipient species based on the de novo Kentish plover reference genome and resequenced populations. We show that the two plover lineages are distinct in both nuclear and mitochondrial genomes. Using model-based coalescence analysis, we found that population sizes of Kentish plover increased whereas white-faced plovers declined during the Last Glaciation Period. Moreover, the two plovers diverged allopatrically, with gene flow occurring after secondary contact. This has resulted in low levels of genome-wide differentiation, although we found evidence of a few highly differentiated genomic regions in both the autosomes and the Z-chromosome. This study illustrates that incipient shorebird species with gene flow after secondary contact can exhibit discrete divergence at specific genomic regions and provides basis to further exploration on the genetic basis of relevant phenotypic traits.

Demographic inference in barn swallows using whole-genome data shows signal for bottleneck and subspecies differentiation during the Holocene

Accounting for historical demographic features is vital for many types of evolutionary inferences, including the estimation of divergence times between closely related populations. In barn swallow, Hirundo rustica, inferring historical population sizes and subspecies divergence times can shed light on the recent co-evolution of this species with humans. Pairwise sequentially Markovian coalescent uncovered population growth beginning on the order of one million years ago-which may reflect the radiation of the broader Hirundo genus-and a more recent population decline. Additionally, we used approximate Bayesian computation to evaluate hypotheses about recent timescale barn swallow demography, including population growth due to human commensalism, and a potential founder event associated with the onset of nesting on human structures. We found signal for a bottleneck event approximately 7,700 years ago, near the time that humans began building substantial structures, although there was considerable uncertainty associated with this estimate. Subspecies differentiation and subsequent growth occurred after the bottleneck in the best-supported model, an order of magnitude more recently than previous estimates in this system. We also compared results obtained from whole-genome sequencing versus reduced representation sequencing, finding many similar results despite substantial allelic dropout in the reduced representation data, which may have affected estimates of some parameters. This study presents the first genetic evidence of a potential barn swallow founder effect and subspecies divergence coinciding with the Holocene, which is an important step in analysing the biogeographical history of a well-known human commensal species.

The genomic landscape at a late stage of stickleback speciation: High genomic divergence interspersed by small localized regions of introgression

Speciation is a continuous process and analysis of species pairs at different stages of divergence provides insight into how it unfolds. Previous genomic studies on young species pairs have revealed peaks of divergence and heterogeneous genomic differentiation. Yet less known is how localised peaks of differentiation progress to genome-wide divergence during the later stages of speciation in the presence of persistent gene flow. Spanning the speciation continuum, stickleback species pairs are ideal for investigating how genomic divergence builds up during speciation. However, attention has largely focused on young postglacial species pairs, with little knowledge of the genomic signatures of divergence and introgression in older stickleback systems. The Japanese stickleback species pair, composed of the Pacific Ocean three-spined stickleback (Gasterosteus aculeatus) and the Japan Sea stickleback (G. nipponicus), which co-occur in the Japanese islands, is at a late stage of speciation. Divergence likely started well before the end of the last glacial period and crosses between Japan Sea females and Pacific Ocean males result in hybrid male sterility. Here we use coalescent analyses and Approximate Bayesian Computation to show that the two species split approximately 0.68-1 million years ago but that they have continued to exchange genes at a low rate throughout divergence. Population genomic data revealed that, despite gene flow, a high level of genomic differentiation is maintained across the majority of the genome. However, we identified multiple, small regions of introgression, occurring mainly in areas of low recombination rate. Our results demonstrate that a high level of genome-wide divergence can establish in the face of persistent introgression and that gene flow can be localized to small genomic regions at the later stages of speciation with gene flow.

Speciation history of a species complex of Primulina eburnea (Gesneriaceae) from limestone karsts of southern China, a biodiversity hot spot

Limestone karsts in southern China are characterized by high edaphic and topographic heterogeneity and host high levels of species richness and endemism. However, the evolutionary mechanisms for generating such biodiversity remain poorly understood. Here, we performed species delimitation, population genetic analyses, simulations of gene flow and analyses of floral morphological traits to infer the geographic history of speciation in a species complex of Primulina eburnea from limestone karsts of southern China. Using Bayesian species delimitation, we determined that there are seven distinct species that correspond well to the putative morphological species. Species tree reconstruction, Structure and Neighbor‐Net analyses all recovered four lineages in agreement with currently species geographic boundaries. High levels of genetic differentiation were observed both within and among species. Isolation–migration coalescent analysis provides evidence for significant but low gene flow among species. Approximate Bayesian computation (ABC) analysis supports a scenario of historical gene flow rather than recent contemporary gene flow for most species divergences. Finally, we found no evidence of divergent selection contributing to population differentiation of a suite of flower traits. These results support the prevalence of allopatric speciation and highlight the role of geographic isolation in the diversification process. At small geographic scales, limited hybridization occurred in the past between proximate populations but did not eliminate species boundaries. We conclude that limited gene flow might have been the predominant evolutionary force in promoting population differentiation and speciation.

Population demographic history of a temperate shrub, Rhododendron weyrichii (Ericaceae), on continental islands of Japan and South Korea

Continental islands provide opportunities for testing the effects of isolation and migration on genetic variation in plant populations. In characteristic of continental islands is that the geographic connections between these islands, which are currently distinguished by seaways, have experienced fluctuations caused by sea‐level changes due to climate oscillations during the Quaternary. Plant populations on the islands have migrated between these islands via the exposed seafloors or been isolated. Here, we examined the demographic history of a temperate shrub, Rhododendron weyrichii, which is distributed in the southwestern parts of the Japanese archipelago and on an island of South Korea, using statistical phylogeographic approaches based on the DNA sequences of two chloroplast and eight nuclear loci in samples analyzed from 18 populations on eight continental islands, and palaeodistribution modeling. Time estimates for four island populations indicate that the durations of vicariance history are different between these populations, and these events have continued since the last glacial or may have predated the last glacial. The constancy or expansion of population sizes on the Japanese islands, and in contrast a bottleneck in population size on the Korean island Jeju, suggests that these islands may have provided different conditions for sustaining populations. The result of palaeodistribution modeling indicates that the longitudinal range of the species as a whole has not changed greatly since the last glacial maximum. These results indicate that exposed seafloors during the glacial period formed both effective and ineffective migration corridors. These findings may shed light on the effects of seafloor exposure on the migration of plants distributed across continental islands.

Sequencing of the genus Arabidopsis identifies a complex history of nonbifurcating speciation and abundant trans-specific polymorphism

The notion of species as reproductively isolated units related through a bifurcating tree implies that gene trees should generally agree with the species tree and that sister taxa should not share polymorphisms unless they diverged recently and should be equally closely related to outgroups. It is now possible to evaluate this model systematically. We sequenced multiple individuals from 27 described taxa representing the entire Arabidopsis genus. Cluster analysis identified seven groups, corresponding to described species that capture the structure of the genus. However, at the level of gene trees, only the separation of Arabidopsis thaliana from the remaining species was universally supported, and, overall, the amount of shared polymorphism demonstrated that reproductive isolation was considerably more recent than the estimated divergence times. We uncovered multiple cases of past gene flow that contradict a bifurcating species tree. Finally, we showed that the pattern of divergence differs between gene ontologies, suggesting a role for selection.

Genomic and phenotypic differentiation of Arabidopsis thaliana along altitudinal gradients in the North Italian Alps

Altitudinal gradients in mountain regions are short-range clines of different environmental parameters such as temperature or radiation. We investigated genomic and phenotypic signatures of adaptation to such gradients in five Arabidopsis thaliana populations from the North Italian Alps that originated from 580 to 2350 m altitude by resequencing pools of 19-29 individuals from each population. The sample includes two pairs of low- and high-altitude populations from two different valleys. High-altitude populations showed a lower nucleotide diversity and negative Tajima's D values and were more closely related to each other than to low-altitude populations from the same valley. Despite their close geographic proximity, demographic analysis revealed that low- and high-altitude populations split between 260 000 and 15 000 years before present. Single nucleotide polymorphisms whose allele frequencies were highly differentiated between low- and high-altitude populations identified genomic regions of up to 50 kb length where patterns of genetic diversity are consistent with signatures of local selective sweeps. These regions harbour multiple genes involved in stress response. Variation among populations in two putative adaptive phenotypic traits, frost tolerance and response to light/UV stress was not correlated with altitude. Taken together, the spatial distribution of genetic diversity reflects a potentially adaptive differentiation between low- and high-altitude populations, whereas the phenotypic differentiation in the two traits investigated does not. It may resemble an interaction between adaptation to the local microhabitat and demographic history influenced by historical glaciation cycles, recent seed dispersal and genetic drift in local populations.

CoMuS: simulating coalescent histories and polymorphic data from multiple species

The simultaneous analysis of intra- and interspecies variation is challenging mainly because our knowledge about patterns of polymorphisms where both intra- and interspecies samples coexist is limited. In this study, we present CoMuS (Coalescent of Multiple Species), a multispecies coalescent software that can simulate intra- and interspecies polymorphisms. CoMuS supports a variety of speciation models and demographic scenarios related to the history of each species. In CoMuS, speciation can be accompanied by either instant or gradual isolation between sister species. Sampling may also occur in the past, and thus, we can study simultaneously extinct and extant species. Our software supports both the infinite- and the finite-site model, with substitution rate heterogeneity among sites and a user-defined proportion of invariable sites. We demonstrate the usage of CoMuS in various applications: species delimitation, software testing, model selection and parameter inference involving present-day and ancestral samples, comparison between gradual and instantaneous isolation models, estimation of speciation time between human and chimpanzee using both intra- and interspecies variation. We expect that CoMuS will be particularly useful for studies where species have been separated recently from their common ancestor and phenomena such as incomplete lineage sorting or introgression still occur.

Contributions of multiple refugia during the last glacial period to current mainland populations of Korean pine (Pinus koraiensis)

The northern microrefugia that existed during the Last Glacial Maximum (LGM) are a key factor in the demographic history of species. Pinus koraiensis has a unique distribution in northeast Asia. The Changbai Mountains and the Korean peninsula (CM/KP) are usually considered to be the LGM refugia for P. koraiensis. However, the Xiaoxingan Range (XR), at the northern part of this species' distribution, is another possible refugium. We used chloroplast sequencing and ten nuclear single-copy gene loci to calculate the genetic diversity pattern of P. koraiensis. The probabilities of a single LGM refugium and of multiple LGM refugia were calculated based on approximate Bayesian computation. The effect of the latitudinal gradient on genetic diversity was not significant. However, unique alleles occurred at low frequencies in CM/KP and XR. A conservative estimate of the coalescence time between CM/KP and XR is 0.4 million years ago, a time prior to the LGM. Gene flow between CM/KP and XR was estimated to be more than one in per generation, an amount that may be sufficient to limit genetic divergence between the regions. Our study strongly supports the hypothesis that XR was another LGM refugium in addition to CM/KP.

Speciation with gene flow in whiptail lizards from a Neotropical xeric biome

Two main hypotheses have been proposed to explain the diversification of the Caatinga biota. The riverine barrier hypothesis (RBH) claims that the São Francisco River (SFR) is a major biogeographic barrier to gene flow. The Pleistocene climatic fluctuation hypothesis (PCH) states that gene flow, geographic genetic structure and demographic signatures on endemic Caatinga taxa were influenced by Quaternary climate fluctuation cycles. Herein, we analyse genetic diversity and structure, phylogeographic history, and diversification of a widespread Caatinga lizard (Cnemidophorus ocellifer) based on large geographical sampling for multiple loci to test the predictions derived from the RBH and PCH. We inferred two well-delimited lineages (Northeast and Southwest) that have diverged along the Cerrado-Caatinga border during the Mid-Late Miocene (6-14 Ma) despite the presence of gene flow. We reject both major hypotheses proposed to explain diversification in the Caatinga. Surprisingly, our results revealed a striking complex diversification pattern where the Northeast lineage originated as a founder effect from a few individuals located along the edge of the Southwest lineage that eventually expanded throughout the Caatinga. The Southwest lineage is more diverse, older and associated with the Cerrado-Caatinga boundaries. Finally, we suggest that C. ocellifer from the Caatinga is composed of two distinct species. Our data support speciation in the presence of gene flow and highlight the role of environmental gradients in the diversification process.

Climate niche differentiation between two passerines despite ongoing gene flow

Niche evolution underpins the generation and maintenance of biological diversity, but niche conservatism, in which niches remain little changed over time in closely related taxa, and the role of ecology in niche evolution are continually debated. To test whether climate niches are conserved in two closely related passerines in East Asia - the vinous-throated (Paradoxornis webbianus) and ashy-throated (P. alphonsianus) parrotbills - we established their potential allopatric and sympatric regions using ecological niche models and compared differences in their climate niches using niche overlap indices in background tests and multivariate statistical analyses. We also used polymorphism data on 44 nuclear genes to infer their divergence demography. We found that these two parrotbills occupy different climate niches, in both their allopatric and potential sympatric regions. Because the potential sympatric region is the area predicted to be suitable for both parrotbills based on the ecological niche models, it can serve as a natural common garden. Therefore, their observed niche differences in this potential sympatry were not simply rendered by phenotypic plasticity and probably had a genetic basis. Our genetic analyses revealed that the two parrotbills are not evolutionarily independent for the most recent part of their divergence history. The two parrotbills diverged c. 856,000 years ago and have had substantial gene flow since a presumed secondary contact c. 290,000 years ago. This study provides an empirical case demonstrating that climate niches may not be homogenized in nascent species in spite of substantial, ongoing gene flow, which in turn suggests a role for ecology in promoting and maintaining diversification among incipient species.

Adaptive evolution of genes involved in the regulation of germline stem cells in Drosophila melanogaster and D. simulans

Population genetic and comparative analyses in diverse taxa have shown that numerous genes involved in reproduction are adaptively evolving. Two genes involved in germline stem cell regulation, bag of marbles (bam) and benign gonial cell neoplasm (bgcn), have been shown previously to experience recurrent, adaptive evolution in both Drosophila melanogaster and D. simulans. Here we report a population genetic survey on eight additional genes involved in germline stem cell regulation in D. melanogaster and D. simulans that reveals all eight of these genes reject a neutral model of evolution in at least one test and one species after correction for multiple testing using a false-discovery rate of 0.05. These genes play diverse roles in the regulation of germline stem cells, suggesting that positive selection in response to several evolutionary pressures may be acting to drive the adaptive evolution of these genes.

Combined next-generation sequencing and morphology reveal fine-scale speciation in Crocodile Skinks (Squamata: Scincidae: Tribolonotus)

Next-generation sequencing has vast potential to revolutionize the fields of phylogenetics and population genetics through its ability to collect genomic scale data sets of thousands of orthologous loci. Despite this potential, other types of data (e.g. morphology, ecology) remain important, particularly for studies endeavouring to delimit species. Here, we integrate next-generation sequencing with morphology to examine divergence between populations of Tribolonotus pseudoponceleti on the islands of Buka and Bougainville in the Solomon Archipelago. We used the Ion Torrent PGM to collect over 648 Mbp of sequence data for 12 samples, representing 1526 loci recovered from all samples, and 3342 were recovered from at least six samples. Genetic structure analyses strongly support the distinctiveness of these two populations, and Bayes factor delimitations decisively select speciation between Buka and Bougainville. Principal components and discriminant function analyses reveal concordant morphological divergence. Finally, demographic analyses via diffusion approximation and approximate Bayesian computation prefer a complex model of mid-Pleistocene divergence with migration, and a later decrease or cessation of migration and population size shift, suggesting a scenario in which migration was enabled by Pleistocene merging of these two islands, and limited when isolated by higher sea levels. Further analysis of four Sanger sequenced loci in IMa2 had limited power to distinguish among models including and excluding migration, but resulted in similar population size and divergence time estimates, although with much broader confidence intervals. This study represents a framework for how next-generation sequencing and morphological data can be combined and leveraged towards validating putative species and testing demographic scenarios for speciation.

Influence of gene flow on divergence dating - implications for the speciation history of Takydromus grass lizards

Dating the time of divergence and understanding speciation processes are central to the study of the evolutionary history of organisms but are notoriously difficult. The difficulty is largely rooted in variations in the ancestral population size or in the genealogy variation across loci. To depict the speciation processes and divergence histories of three monophyletic Takydromus species endemic to Taiwan, we sequenced 20 nuclear loci and combined with one mitochondrial locus published in GenBank. They were analysed by a multispecies coalescent approach within a Bayesian framework. Divergence dating based on the gene tree approach showed high variation among loci, and the divergence was estimated at an earlier date than when derived by the species-tree approach. To test whether variations in the ancestral population size accounted for the majority of this variation, we conducted computer inferences using isolation-with-migration (IM) and approximate Bayesian computation (ABC) frameworks. The results revealed that gene flow during the early stage of speciation was strongly favoured over the isolation model, and the initiation of the speciation process was far earlier than the dates estimated by gene- and species-based divergence dating. Due to their limited dispersal ability, it is suggested that geographical isolation may have played a major role in the divergence of these Takydromus species. Nevertheless, this study reveals a more complex situation and demonstrates that gene flow during the speciation process cannot be overlooked and may have a great impact on divergence dating. By using multilocus data and incorporating Bayesian coalescence approaches, we provide a more biologically realistic framework for delineating the divergence history of Takydromus.

ABC inference of multi-population divergence with admixture from unphased population genomic data

Rapidly developing sequencing technologies and declining costs have made it possible to collect genome-scale data from population-level samples in nonmodel systems. Inferential tools for historical demography given these data sets are, at present, underdeveloped. In particular, approximate Bayesian computation (ABC) has yet to be widely embraced by researchers generating these data. Here, we demonstrate the promise of ABC for analysis of the large data sets that are now attainable from nonmodel taxa through current genomic sequencing technologies. We develop and test an ABC framework for model selection and parameter estimation, given histories of three-population divergence with admixture. We then explore different sampling regimes to illustrate how sampling more loci, longer loci or more individuals affects the quality of model selection and parameter estimation in this ABC framework. Our results show that inferences improved substantially with increases in the number and/or length of sequenced loci, while less benefit was gained by sampling large numbers of individuals. Optimal sampling strategies given our inferential models included at least 2000 loci, each approximately 2 kb in length, sampled from five diploid individuals per population, although specific strategies are model and question dependent. We tested our ABC approach through simulation-based cross-validations and illustrate its application using previously analysed data from the oak gall wasp, Biorhiza pallida.

Rabbit genome analysis reveals a polygenic basis for phenotypic change during domestication

The genetic changes underlying the initial steps of animal domestication are still poorly understood. We generated a high-quality reference genome for the rabbit and compared it to resequencing data from populations of wild and domestic rabbits. We identified more than 100 selective sweeps specific to domestic rabbits but only a relatively small number of fixed (or nearly fixed) single-nucleotide polymorphisms (SNPs) for derived alleles. SNPs with marked allele frequency differences between wild and domestic rabbits were enriched for conserved noncoding sites. Enrichment analyses suggest that genes affecting brain and neuronal development have often been targeted during domestication. We propose that because of a truly complex genetic background, tame behavior in rabbits and other domestic animals evolved by shifts in allele frequencies at many loci, rather than by critical changes at only a few domestication loci.

Association mapping, patterns of linkage disequilibrium and selection in the vicinity of the PHYTOCHROME C gene in pearl millet

Linkage analysis confirmed the association in the region of PHYC in pearl millet. The comparison of genes found in this region suggests that PHYC is the best candidate. Major efforts are currently underway to dissect the phenotype-genotype relationship in plants and animals using existing populations. This method exploits historical recombinations accumulated in these populations. However, linkage disequilibrium sometimes extends over a relatively long distance, particularly in genomic regions containing polymorphisms that have been targets for selection. In this case, many genes in the region could be statistically associated with the trait shaped by the selected polymorphism. Statistical analyses could help in identifying the best candidate genes into such a region where an association is found. In a previous study, we proposed that a fragment of the PHYTOCHROME C gene (PHYC) is associated with flowering time and morphological variations in pearl millet. In the present study, we first performed linkage analyses using three pearl millet F2 families to confirm the presence of a QTL in the vicinity of PHYC. We then analyzed a wider genomic region of ~100 kb around PHYC to pinpoint the gene that best explains the association with the trait in this region. A panel of 90 pearl millet inbred lines was used to assess the association. We used a Markov chain Monte Carlo approach to compare 75 markers distributed along this 100-kb region. We found the best candidate markers on the PHYC gene. Signatures of selection in this region were assessed in an independent data set and pointed to the same gene. These results foster confidence in the likely role of PHYC in phenotypic variation and encourage the development of functional studies.

Demographic divergence history of pied flycatcher and collared flycatcher inferred from whole-genome re-sequencing data

Profound knowledge of demographic history is a prerequisite for the understanding and inference of processes involved in the evolution of population differentiation and speciation. Together with new coalescent-based methods, the recent availability of genome-wide data enables investigation of differentiation and divergence processes at unprecedented depth. We combined two powerful approaches, full Approximate Bayesian Computation analysis (ABC) and pairwise sequentially Markovian coalescent modeling (PSMC), to reconstruct the demographic history of the split between two avian speciation model species, the pied flycatcher and collared flycatcher. Using whole-genome re-sequencing data from 20 individuals, we investigated 15 demographic models including different levels and patterns of gene flow, and changes in effective population size over time. ABC provided high support for recent (mode 0.3 my, range <0.7 my) species divergence, declines in effective population size of both species since their initial divergence, and unidirectional recent gene flow from pied flycatcher into collared flycatcher. The estimated divergence time and population size changes, supported by PSMC results, suggest that the ancestral species persisted through one of the glacial periods of middle Pleistocene and then split into two large populations that first increased in size before going through severe bottlenecks and expanding into their current ranges. Secondary contact appears to have been established after the last glacial maximum. The severity of the bottlenecks at the last glacial maximum is indicated by the discrepancy between current effective population sizes (20,000-80,000) and census sizes (5-50 million birds) of the two species. The recent divergence time challenges the supposition that avian speciation is a relatively slow process with extended times for intrinsic postzygotic reproductive barriers to evolve. Our study emphasizes the importance of using genome-wide data to unravel tangled demographic histories. Moreover, it constitutes one of the first examples of the inference of divergence history from genome-wide data in non-model species.

Demographic inference reveals African and European admixture in the North American Drosophila melanogaster population

Drosophila melanogaster spread from sub-Saharan Africa to the rest of the world colonizing new environments. Here, we modeled the joint demography of African (Zimbabwe), European (The Netherlands), and North American (North Carolina) populations using an approximate Bayesian computation (ABC) approach. By testing different models (including scenarios with continuous migration), we found that admixture between Africa and Europe most likely generated the North American population, with an estimated proportion of African ancestry of 15%. We also revisited the demography of the ancestral population (Africa) and found-in contrast to previous work-that a bottleneck fits the history of the population of Zimbabwe better than expansion. Finally, we compared the site-frequency spectrum of the ancestral population to analytical predictions under the estimated bottleneck model.

Population and sex differences in Drosophila melanogaster brain gene expression

Background

Changes in gene regulation are thought to be crucial for the adaptation of organisms to their environment. Transcriptome analyses can be used to identify candidate genes for ecological adaptation, but can be complicated by variation in gene expression between tissues, sexes, or individuals. Here we use high-throughput RNA sequencing of a single Drosophila melanogaster tissue to detect brain-specific differences in gene expression between the sexes and between two populations, one from the ancestral species range in sub-Saharan Africa and one from the recently colonized species range in Europe.

Results

Relatively few genes (<100) displayed sexually dimorphic expression in the brain, but there was an enrichment of sex-biased genes, especially male-biased genes, on the X chromosome. Over 340 genes differed in brain expression between flies from the African and European populations, with the inter-population divergence being highly correlated between males and females. The differentially expressed genes included those involved in stress response, olfaction, and detoxification. Expression differences were associated with transposable element insertions at two genes implicated in insecticide resistance (Cyp6g1 and CHKov1).

Conclusions

Analysis of the brain transcriptome revealed many genes differing in expression between populations that were not detected in previous studies using whole flies. There was little evidence for sex-specific regulatory adaptation in the brain, as most expression differences between populations were observed in both males and females. The enrichment of genes with sexually dimorphic expression on the X chromosome is consistent with dosage compensation mechanisms affecting sex-biased expression in somatic tissues.