Evaluating mechanisms of diversification in a Guineo‐Congolian tropical forest frog using demographic model selection

Ecological speciation with gene flow followed initial large-scale geographic speciation in the enigmatic afroalpine giant senecios (Dendrosenecio)

Mountains have highly heterogeneous environments that generate ample opportunities for lineage differentiation through ecological adaptation, geographic isolation and secondary contact. The geographic and ecological isolation of the afroalpine vegetation fragments on the East African mountain tops makes them an excellent system to study speciation. The initial diversification within the afroalpine endemic genus Dendrosenecio was shown to occur via allopatric divergence among four isolated mountain groups, but the potential role of ecological speciation within these groups and the role of gene flow in speciation remained uncertain. Here we extend the sampling of Dendrosenecio and use phylogenomics to assess the importance of gene flow in the diversification of the genus. Then, population genomics, demographic modelling and habitat differentiation analyses are used to study ecological speciation in two sister species occurring on Mount Kenya. We found that two sympatric sister species on Mt Kenya occupy distinct microhabitats, and our analyses support that they originated in situ via ecological speciation with gene flow. In addition, we obtained signals of admixture history between mountain groups. Taken together, these results suggest that geographic isolation shaped main lineages, while ecologically mediated speciation occurred within a single mountain.

Novel island species elucidate a species complex of Neotropical crocodiles

The evolutionary history of Neotropical crocodiles has remained elusive. They inhabit a broad geographic range with populations spanning from coastal, inland, and insular locations. Using a selection of natural insular, coastal, and one inland population of C. acutus, coastal C. moreletii, and the single surviving population of C. rhombifer, we discovered a remarkable genetic diversity for the group. Moreover, geometric morphometric results of skull shapes shows that these crocodylus species span a morphological cline. We recovered a high genetic differentiation between C. moreletii, C. rhombifer, and five clusters of C. acutus. The genetic and geographic differences among the C. acutus clusters were used to suggest these may be a species complex. Several ecological, morphological and genetics traits are identified in the well-studied populations from Banco Chinchorro and Cozumel islands off the Mexican Yucatan Peninsula to support discrete species designations for these populations. This work suggests the presence of rapid, recent evolution of several cryptic Crocodylus species throughout the Neotropics.

Plio-Pleistocene Climatic Fluctuations and Divergence With Gene Flow Drive Continent-Wide Diversification in an African Bird

Diversification mechanisms in Sub-Saharan Africa have long attracted research interest, with varying support for either allopatric or parapatric models of speciation. However, studies have seldom been performed across the entire continent, a scale which could elucidate the relative importance of allopatric and parapatric models of divergence. To shed light on continental-scale patterns of African biogeography and diversification, we investigated the historical demography of a bird with a continent-wide distribution in Sub-Saharan Africa, the Yellow-Rumped Tinkerbird, Pogoniulus bilineatus. We sampled populations from across the continent and, using genomic data, assessed genetic diversity, structure, and differentiation, reconstructed the phylogeny, and performed alternative demographic model selection between neighbouring clade pairs. We uncovered substantial genetic structure and differentiation patterns which corroborated the phylogenetic topology. Structure was chiefly influenced by the arid corridor, a postulated biogeographical barrier in Sub-Saharan Africa. Moreover, peak genetic diversities coincided with postulated refugial areas while demographic reconstructions between genetic lineages supported allopatric models consistent with the Pleistocene Forest Refuge hypothesis. However, within lineages, divergence with gene flow was supported. Continent-wide patterns of diversification involve an integration of both allopatric and parapatric mechanisms, with a role for both periods of divergence in isolation and across ecological gradients. Furthermore, our study emphasises the importance of the arid corridor as a primary biogeographical feature across which diversification occurs, yet one that has hitherto received scant attention regarding its importance in avian diversification in Sub-Saharan Africa.

Demographic History, Genetic Load, and the Efficacy of Selection in the Globally Invasive Mosquito Aedes aegypti

Aedes aegypti is the main vector species of yellow fever, dengue, Zika, and chikungunya. The species is originally from Africa but has experienced a spectacular expansion in its geographic range to a large swath of the world, and the demographic effects of which have remained largely understudied. In this report, we examine whole-genome sequences from six countries in Africa, North America, and South America to investigate the demographic history of the spread of A. aegypti into the Americas and its impact on genomic diversity and deleterious genetic load. In the Americas, we observe patterns of strong population structure consistent with relatively low (but probably nonzero) levels of gene flow but occasional long-range dispersal and/or recolonization events. We also find evidence that the colonization of the Americas has resulted in introduction bottlenecks. However, while each sampling location shows evidence of a past population contraction and subsequent recovery, our results suggest that the bottlenecks in America have led to a reduction in genetic diversity of only ∼35% relative to African populations, and the American samples have retained high levels of genetic diversity (expected heterozygosity of ∼0.02 at synonymous sites). We additionally find that American populations of aegypti have experienced only a minor reduction in the efficacy of selection, with evidence for both an accumulation of deleterious alleles and some purging of strongly deleterious alleles. These results exemplify how an invasive species can expand its range with remarkable genetic resilience in the face of strong eradication pressure.

Pangenomics Links Boll Weevil Divergence With Ancient Mesoamerican Cotton Cultivation

The boll weevil, Anthonomus grandis grandis Boheman, and thurberia weevil, Anthonomus grandis thurberiae Pierce, together comprise a species complex that ranges throughout Mexico, the southwestern regions of the United States and parts of South America. The boll weevil is a historically damaging and contemporaneously threatening pest to commercial upland cotton, Gossypium hirsutum L. (Malvales: Malvaceae), whereas the thurberia weevil is regarded as an innocuous non-pest subspecies that is mostly found on non-cultivated Thurber's or Arizona cotton, Gossypium thurberi L., throughout its native range in western Mexico and the southwestern United States. Recent independent analyses, using mitochondrial and whole-genome markers, have suggested the independent evolution of these lineages is more attributable to geographic isolation than biotic factors. We suggest a combination of drivers after employing comparative genomic, population genetic and pangenome methodologies to identify large and small polymorphisms. By leveraging genetic differences, we determined 39,310 diagnostic loci between the subspecies, find genes under selection, and model the subspecies' shared and unique evolutionary history. Interestingly, structural variations capture a large proportion of genes at the population level and demographic reconstruction suggests a split between approximately 3,320-16,300 before present (YBP), which coincides with cotton cultivation in Mesoamerica, approximately 3,000-5,000 YBP. Observed polymorphisms are enriched for reproductive, regulatory, and metabolic genes, which may be attributed to the subspecies split and coevolution with cultivated cotton. Our results demonstrate the utility of a holistic, comparative framework utilising small and large polymorphisms to reconstruct demography and identify genetic novelty via pangenomics.

Population genomics of the gametophyte-only fern Vittaria appalachiana provides insights into clonal plant evolution

How asexually reproducing organisms maintain genetic diversity and adaptive potential is a long-standing question in evolutionary biology. Asexual lineages have historically been thought of as evolutionary dead ends, yet some exhibit remarkable persistence through time. The gametophyte-only fern Vittaria appalachiana is a clonal eukaryote, the focus of extensive study due to its peculiar habit and life history, and is an excellent system to explore the consequences of asexuality. Using reduced representation sequencing and life cycle simulations, we assess theoretical expectations for genomic consequences of long-term asexual reproduction and test hypotheses about its origin and demographic history. We show that V. appalachiana colonies are not patches of single genotypes but are mosaics of genetic diversity, and the accumulation of mutations in the absence of recombination plays an important role in driving this diversity. We identify increased genomic variation, excess heterozygosity, decreased population differentiation, and increased effective population size, all of which are consistent with the expectations for prolonged clonality. Our analyses support the hypothesis that the loss of sexual reproduction in V. appalachiana occurred during the Last Glacial Maximum. Our results from empirical and simulation-based analyses illuminate how an asexual eukaryote generates, retains, and partitions genomic diversity.

A highland-adaptation variant near MCUR1 reduces its transcription and attenuates erythrogenesis in Tibetans

To identify genomic regions subject to positive selection that might contain genes involved in high-altitude adaptation (HAA), we performed a genome-wide scan by whole-genome sequencing of Tibetan highlanders and Han lowlanders. We revealed a collection of candidate genes located in 30 genomic loci under positive selection. Among them, MCUR1 at 6p23 was a novel pronounced candidate. By single-cell RNA sequencing and comprehensive functional studies, we demonstrated that MCUR1 depletion leads to impairment of erythropoiesis under hypoxia and normoxia. Mechanistically, MCUR1 knockdown reduced mitochondrial Ca2+ uptake and then concomitantly increased cytosolic Ca2+ levels, which thereby reduced erythropoiesis via the CAMKK2-AMPK-mTOR axis. Further, we revealed rs61644582 at 6p23 as an expression quantitative trait locus for MCUR1 and a functional variant that confers an allele-specific transcriptional regulation of MCUR1. Overall, MCUR1-mediated mitochondrial Ca2+ homeostasis is highlighted as a novel regulator of erythropoiesis, deepening our understanding of the genetic mechanism of HAA.

Inferring the Demographic History of Aye-Ayes (Daubentonia madagascariensis) from High-Quality, Whole-Genome, Population-Level Data

The nocturnal aye-aye, Daubentonia madagascariensis, is one of the most elusive lemurs on the island of Madagascar. The timing of its activity and arboreal lifestyle has generally made it difficult to obtain accurate assessments of population size using traditional census methods. Therefore, alternative estimates provided by population genetic inference are essential for yielding much needed information for conservation measures and for enabling ecological and evolutionary studies of this species. Here, we utilize genomic data from 17 individuals-including 5 newly sequenced, high-coverage genomes-to estimate this history. Essential to this estimation are recently published annotations of the aye-aye genome which allow for variation at putatively neutral genomic regions to be included in the estimation procedures, and regions subject to selective constraints, or in linkage to such sites, to be excluded owing to the biasing effects of selection on demographic inference. By comparing a variety of demographic estimation tools to develop a well-supported model of population history, we find strong support for two demes, separating northern Madagascar from the rest of the island. Additionally, we find that the aye-aye has experienced two severe reductions in population size. The first occurred rapidly, ∼3,000 to 5,000 years ago, and likely corresponded with the arrival of humans to Madagascar. The second occurred over the past few decades and is likely related to substantial habitat loss, suggesting that the species is still undergoing population decline and remains at great risk for extinction.

Genomic Architecture Underlying the Striking Colour Variation in the Presence of Gene Flow for the Guinan Toad-Headed Lizard

How divergence occurs between closely related organisms in the absence of geographic barriers to gene flow stands as one of the long-standing questions in evolutionary biology. Previous studies suggested that the interplay between selection, gene flow and recombination strongly affected the process of divergence with gene flow. However, the extent to which these forces interact to drive divergence remains largely ambiguous. Guinan toad-headed lizards (Phrynocephalus guinanensis) in the Mugetan Desert exhibit striking colour differences from lizards outside the desert and provide an excellent model to address this question. Through extensive sampling and whole genome sequencing, we obtained genotypes for 191 samples from 14 populations inside and outside the desert. Despite the colour differences, continuous and asymmetric gene flow was detected across the desert border. More importantly, 273 highly diverged regions (HDRs) were identified between them, accounting only for 0.47% of the genome but widely distributed across 20 (out of the total 24) chromosomes. Strong signatures of selection were identified in HDRs, and local recombination rates were repressed. Furthermore, five HDRs exhibited significantly higher divergence, which contained key genes associated with crucial functions in animal coloration, including pteridine and melanocyte pigmentation. Genes related to retinal cells and steroid hormones were identified in other HDRs, which might have also contributed to the formation of colour variation in the presence of gene flow. This study provided novel insights into the understanding of the evolutionary mechanisms of genetic divergence in the presence of gene flow.

Testing the genomic overlap between intraspecific mating traits and interspecific mating barriers

Differences in interspecific mating traits, such as male sexual signals and female preferences, often evolve quickly as initial barriers to gene flow between nascent lineages, and they may also strengthen such barriers during secondary contact via reinforcement. However, it is an open question whether loci contributing to intraspecific variation in sexual traits are co-opted during the formation and strengthening of mating barriers between species. To test this, we used a population genomics approach in natural populations of Australian cricket sister species that overlap in a contact zone: Teleogryllus oceanicus and Teleogryllus commodus. First, we identified loci associated with intraspecific variation in T. oceanicus mating signals: advertisement song and cuticular hydrocarbon (CHC) pheromones. We then separately identified candidate interspecific barrier loci between the species. Genes showing elevated allelic divergence between species were enriched for neurological functions, indicating potential behavioral rewiring. Only two CHC-associated genes overlapped with these interspecific candidate barrier loci, and intraspecific CHC loci showed signatures of being under strong selective constraints between species. In contrast, 10 intraspecific song-associated genes showed high genetic differentiation between T. commodus and T. oceanicus, and 2 had signals of high genomic divergence. The overall lack of shared loci in intra vs. interspecific comparisons of mating trait and candidate barrier loci is consistent with limited co-option of the genetic architecture of interspecific mating signals during the establishment and maintenance of reproductive isolation.

Understanding Species Boundaries that Arise from Complex Histories: Gene Flow Across the Speciation Continuum in the Spotted Whiptail Lizards

-Gene flow between diverging lineages challenges the resolution of species boundaries and the understanding of evolutionary history in recent radiations. Here, we integrate phylogenetic and coalescent tools to resolve reticulate patterns of diversification and use a perspective focused on evolutionary mechanisms to distinguish interspecific and intraspecific taxonomic variation. We use this approach to resolve the systematics for one of the most intensively studied but difficult to understand groups of reptiles: the spotted whiptail lizards of the genus Aspidoscelis (A. gularis complex). Whiptails contain the largest number of unisexual species known within any vertebrate group and the spotted whiptail complex has played a key role in the generation of this diversity through hybrid speciation. Understanding lineage boundaries and the evolutionary history of divergence and reticulation within this group is therefore key to understanding the generation of unisexual diversity in whiptails. Despite this importance, long-standing confusion about their systematics has impeded understanding of which gonochoristic species have contributed to the formation of unisexual lineages. Using reduced representation genomic data, we resolve patterns of divergence and gene flow within the spotted whiptails and clarify patterns of hybrid speciation. We find evidence that biogeographically structured ecological and environmental variation has been important in morphological and genetic diversification, as well as the maintenance of species boundaries in this system. Our study elucidates how gene flow among lineages and the continuous nature of speciation can bias the practice of species delimitation and lead taxonomists operating under different frameworks to different conclusions (here we propose that a 2 species arrangement best reflects our current understanding). In doing so, this study provides conceptual and methodological insights into approaches to resolving diversification patterns and species boundaries in rapid radiations with complex histories, as well as long-standing taxonomic challenges in the field of systematic biology.

Genome-Wide Data Uncover Cryptic Diversity With Multiple Reticulation Events in the Balkan-Anatolian Cardamine (Brassicaceae) Species Complex

Plant species diversity may be considerably underestimated, especially in evolutionarily complex genera and in diversity hotspots that have enabled long-term species persistence and diversification, such as the Balkan Peninsula. Here, we address the topic of underexplored plant diversity and underlying evolutionary and biogeographic processes by investigating the hygrophytic mountain species complex of Cardamine acris s.l. distributed in the Balkans (three subspecies within C. acris) and northwestern Anatolia (C. anatolica). We performed a series of phylogenetic and phylogeographic analyses based on restriction-site associated DNA sequencing (RADseq) and target enrichment (Hyb-Seq) data in combination with habitat suitability modelling. We found C. anatolica as a clade nested within the Balkan C. acris, probably resulting from a founder event, and uncovered three allopatric cryptic lineages within C. acris subsp. acris, allowing us to recognise a total of six entities in this complex. We observed the deepest genetic split within C. acris subsp. acris in the western Balkans, which was at odds with taxonomy and showed no distribution gap. We inferred vicariance as the most likely process for population divergence in the Balkans, accompanied by gene flow between the recognised entities, which was consistent with the modelled habitat suitability dynamics. Furthermore, we discovered several polyploid populations in C. acris, representing both pure intra- and inter-lineage hybrid polyploids, but detected only minor traces of hybridization with related congeners. Overall, our results illustrate that diverse evolutionary processes may influence the history of mountain plant species in the Balkan Peninsula, including vicariance, reticulation, polyploidization and cryptic diversification.

Demographic History and Adaptive Evolution of Indo-Pacific Bottlenose Dolphins (Tursiops aduncus) in Western Australia

Demographic processes can substantially affect a species' response to changing ecological conditions, necessitating the combined consideration of genetic responses to environmental variables and neutral genetic variation. Using a seascape genomics approach combined with population demographic modelling, we explored the interplay of demographic and environmental factors that shaped the current population structure in Indo-Pacific bottlenose dolphins (Tursiops aduncus) along the Western Australian coastline. We combined large-scale environmental data gathered via remote sensing with RADseq genomic data from 133 individuals at 19 sampling sites. Using population genetic and outlier detection analyses, we identified three distinct genetic clusters, coinciding with tropical, subtropical and temperate provincial bioregions. In contrast to previous studies, our demographic models indicated that populations occupying the paleo-shoreline split into two demographically independent lineages before the last glacial maximum (LGM). A subsequent split after the LGM 12-15 kya gave rise to the Shark Bay population, thereby creating the three currently observed clusters. Although multi-locus heterozygosity declined from north to south, dolphins from the southernmost cluster inhabiting temperate waters had higher heterozygosity in potentially adaptive loci compared to dolphins from subtropical and tropical waters. These findings suggest ongoing adaptation to cold-temperate waters in the southernmost cluster, possibly linked to distinct selective pressures between the different bioregions. Our study demonstrated that in the marine realm, without apparent physical boundaries, only a combined approach can fully elucidate the intricate environmental and genetic interactions shaping the evolutionary trajectory of marine mammals.

Population structure and demographic history of two highly-trafficked species of pangolin in the Congo Basin

African pangolins are hunted for their meat and for use in local traditional medicine, as well as for their scales, which are trafficked internationally, especially to growing Asian markets. Pangolin's population genetic structure can be used to trace the geographic origins of trafficked scales, but substantial sampling gaps across pangolins' ranges hinder these efforts. In this study, we documented population structure and dynamics in the two species of African pangolin, the white-bellied pangolin (Phataginus tricuspis) and the giant pangolin (Smutsia gigantea) in the underexplored Republic of Congo. Using the mitochondrial control region and two nuclear markers (beta-fibrinogen and titin), we identified high genetic diversity in both species. We document a distinct mitochondrial lineage of the white-bellied pangolin, which was most likely shaped by river barriers together with dynamics of forest refugia related to the climatic shifts during the Pleistocene. We detected population growth in the white-bellied pangolin coinciding with a dry period during the Pleistocene, suggesting some ability for this typically forest-dwelling species to persist under diverse environmental conditions. Using landscape genetics, we found all but one of the pangolins we sampled at bush meat markets originated locally. A single individual appeared to have been imported to Congo from Cameroon. These findings significantly contribute to our understanding of pangolin population biology and local trade dynamics. In addition, our data from a previously unstudied part of pangolins' ranges will help us to better understand international wildlife trafficking patterns and to target conservation and protection strategies for these highly vulnerable species.

Population genomic evidence of a putative 'far-west' African cryptic taxon in the Anopheles gambiae complex

The two main Afrotropical malaria vectors - Anopheles coluzzii and An. gambiae - are genetically distinct and reproductively isolated across West Africa. However, populations at the western extreme of their range are assigned as "intermediate" between the two species by whole genome sequence (WGS) data, and as hybrid forms by conventional molecular diagnostics. By exploiting WGS data from 1190 specimens collected across west Africa via the Anopheles gambiae 1000 Genomes network, we identified a putative taxon in the far-west (provisionally named Bissau molecular form), which did not arise by admixture but rather may have originated at the same time as the split between An. coluzzii and An. gambiae. Intriguingly, this taxon lacks insecticide resistance mechanisms commonly observed in the two main species. These findings lead to a change of perspective on malaria vector species in the far-west region with potential for epidemiological implications, and a new challenge for genetic-based mosquito control approaches.

Genomic divergence and mutation load in the Begonia masoniana complex from limestone karsts

Understanding genome-wide diversity, inbreeding, and the burden of accumulated deleterious mutations in small and isolated populations is essential for predicting and enhancing population persistence and resilience. However, these effects are rarely studied in limestone karst plants. Here, we re-sequenced the nuclear genomes of 62 individuals of the Begonia masoniana complex (B. liuyanii, B. longgangensis, B. masoniana and B. variegata) and investigated genomic divergence and genetic load for these four species. Our analyses revealed four distinct clusters corresponding to each species within the complex. Notably, there was only limited admixture between B. liuyanii and B. longgangensis occurring in overlapping geographic regions. All species experienced historical bottlenecks during the Pleistocene, which were likely caused by glacial climate fluctuations. We detected an asymmetric historical gene flow between group pairs within this timeframe, highlighting a distinctive pattern of interspecific divergence attributable to karst geographic isolation. We found that isolated populations of B. masoniana have limited gene flow, the smallest recent population size, the highest inbreeding coefficients, and the greatest accumulation of recessive deleterious mutations. These findings underscore the urgency to prioritize conservation efforts for these isolated population. This study is among the first to disentangle the genetic differentiation and specific demographic history of karst Begonia plants at the whole-genome level, shedding light on the potential risks associated with the accumulation of deleterious mutations over generations of inbreeding. Moreover, our findings may facilitate conservation planning by providing critical baseline genetic data and a better understanding of the historical events that have shaped current population structure of rare and endangered karst plants.

Divergent lineages in a young species: The case of datilillo (Yucca valida), a broadly distributed plant from the Baja California Peninsula

PREMISE

Globally, barriers triggered by climatic changes have caused habitat fragmentation and population allopatric divergence. Across North America, oscillations during the Quaternary have played important roles in the distribution of wildlife. Notably, diverse plant species from the Baja California Peninsula in western North America, isolated during the Pleistocene glacial-interglacial cycles, exhibit strong genetic structure and highly concordant divergent lineages across their ranges. A representative plant genus of the peninsula is Yucca, with Y. valida having the widest range. Although a dominant species, it has an extensive distribution discontinuity between 26° N and 27° N, suggesting restricted gene flow. Moreover, historical distribution models indicate the absence of an area with suitable conditions for the species during the Last Interglacial, making it an interesting model for studying genetic divergence.

METHODS

We assembled 4411 SNPs from 147 plants of Y. valida throughout its range to examine its phylogeography to identify the number of genetic lineages, quantify their genetic differentiation, reconstruct their demographic history and estimate the age of the species.

RESULTS

Three allopatric lineages were identified based on the SNPs. Our analyses support that genetic drift is the driver of genetic differentiation among these lineages. We estimated an age of less than 1 million years for the common ancestor of Y. valida and its sister species.

CONCLUSIONS

Habitat fragmentation caused by climatic changes, low dispersal, and an extensive geographical range gap acted as cumulative mechanisms leading to allopatric divergence in Y. valida.

Population genomics of Fusarium graminearum isolates from the Americas

Fusarium head blight (FHB) is a major disease of wheat and barley worldwide and is caused by different species in the genus Fusarium, Fusarium graminearum being the most important. We conducted population genomics analyses using SNPs obtained through genotyping by sequencing of over 500 isolates of F. graminearum from the US Upper Midwest, New York, Louisiana, and Uruguay. PCA and STRUCTURE analyses group our isolates into four previously described populations: NA1, NA2, Southern Louisiana (SLA) and Gulf Coast (GC). Some isolates were not assigned to populations because of mixed ancestry. Population structure was associated with toxin genotype and geographic origin. The NA1, NA2, and SLA populations are differentiated (FST 0.385 - 0.551) but the presence of admixed isolates indicates that the populations are not reproductively isolated. Patterns of linkage disequilibrium (LD) decay suggest frequent recombination within populations. Fusarium graminearum populations from the US have great evolutionary potential given the high recombination rate and a large proportion of admixed isolates. The NA1, NA2, and Southern Louisiana (SLA) populations separated from their common ancestral population roughly at the same time in the past and are evolving with moderate levels of subsequent gene flow between them. Genome-wide selection scans in all three populations revealed outlier regions with the strongest signatures of recent positive natural selection. These outlier regions include many genes with unknown function and some genes with known roles in plant-microbe interaction, fungicide/drug resistance, cellular transport and genes that are related to cellular organelles. Only a very small proportion of outlier regions are shared as outliers among the three populations, suggesting unique host-pathogen interactions and environmental adaptation.

Modeling biases from low-pass genome sequencing to enable accurate population genetic inferences

Low-pass genome sequencing is cost-effective and enables analysis of large cohorts. However, it introduces biases by reducing heterozygous genotypes and low-frequency alleles, impacting subsequent analyses such as demographic history inference. We developed a probabilistic model of low-pass biases from the Genome Analysis Toolkit (GATK) multi-sample calling pipeline, and we implemented it in the population genomic inference software dadi. We evaluated the model using simulated low-pass datasets and found that it alleviated low-pass biases in inferred demographic parameters. We further validated the model by downsampling 1000 Genomes Project data, demonstrating its effectiveness on real data. Our model is widely applicable and substantially improves model-based inferences from low-pass population genomic data.

Complex patterns of genetic population structure in the mouthbrooding marine catfish, Bagre marinus, in the Gulf of Mexico and U.S. Atlantic

Patterns of genetic variation reflect interactions among microevolutionary forces that vary in strength with changing demography. Here, patterns of variation within and among samples of the mouthbrooding gafftopsail catfish (Bagre marinus, Family Ariidae) captured in the U.S. Atlantic and throughout the Gulf of Mexico were analyzed using genomics to generate neutral and non-neutral SNP data sets. Because genomic resources are lacking for ariids, linkage disequilibrium network analysis was used to examine patterns of putatively adaptive variation. Finally, historical demographic parameters were estimated from site frequency spectra. The results show four differentiated groups, corresponding to the (1) U.S. Atlantic, and the (2) northeastern, (3) northwestern, and (4) southern Gulf of Mexico. The non-neutral data presented two contrasting signals of structure, one due to increases in diversity moving west to east and north to south, and another to increased heterozygosity in the Atlantic. Demographic analysis suggested that recently reduced long-term effective population size in the Atlantic is likely an important driver of patterns of genetic variation and is consistent with a known reduction in population size potentially due to an epizootic. Overall, patterns of genetic variation resemble that of other fishes that use the same estuarine habitats as nurseries, regardless of the presence/absence of a larval phase, supporting the idea that adult/juvenile behavior and habitat are important predictors of contemporary patterns of genetic structure.

Computationally Efficient Demographic History Inference from Allele Frequencies with Supervised Machine Learning

Inferring past demographic history of natural populations from genomic data is of central concern in many studies across research fields. Previously, our group had developed dadi, a widely used demographic history inference method based on the allele frequency spectrum (AFS) and maximum composite-likelihood optimization. However, dadi's optimization procedure can be computationally expensive. Here, we present donni (demography optimization via neural network inference), a new inference method based on dadi that is more efficient while maintaining comparable inference accuracy. For each dadi-supported demographic model, donni simulates the expected AFS for a range of model parameters then trains a set of Mean Variance Estimation neural networks using the simulated AFS. Trained networks can then be used to instantaneously infer the model parameters from future genomic data summarized by an AFS. We demonstrate that for many demographic models, donni can infer some parameters, such as population size changes, very well and other parameters, such as migration rates and times of demographic events, fairly well. Importantly, donni provides both parameter and confidence interval estimates from input AFS with accuracy comparable to parameters inferred by dadi's likelihood optimization while bypassing its long and computationally intensive evaluation process. donni's performance demonstrates that supervised machine learning algorithms may be a promising avenue for developing more sustainable and computationally efficient demographic history inference methods.

Survival in nunatak and peripheral glacial refugia of three alpine plant species is partly predicted by altitudinal segregation

Mountain biota survived the Quaternary cold stages most probably in peripheral refugia and/or ice-free peaks within ice-sheets (nunataks). While survival in peripheral refugia has been broadly demonstrated, evidence for nunatak refugia is still scarce. We generated RADseq data from three mountain plant species occurring at different elevations in the southeastern European Alps to investigate the role of different glacial refugia during the Last Glacial Maximum (LGM). We tested the following hypotheses. (i) The deep Piave Valley forms the deepest genetic split in the species distributed across it, delimiting two peripheral refugia. (ii) The montane to alpine species Campanula morettiana and Primula tyrolensis survived the LGM in peripheral refugia, while high-alpine to subnival Saxifraga facchinii likely survived in several nunatak refugia. (iii) The lower elevation species suffered a strong population decline during the LGM. By contrast, the higher elevation species shows long-term stability of population sizes due to survival on permanently ice-free peaks and small population sizes at present. We found peripheral refugia on both sides of the Piave Valley, which acted as a major genetic barrier. Demographic modelling confirmed nunatak survival not only for S. facchinii but also for montane to alpine C. morettiana. Altitudinal segregation influenced the species' demographic fluctuations, with the lower elevation species showing a significant population increase at the end of the LGM, and the higher elevation species either showing decrease towards the present or stable population sizes with a short bottleneck. Our results highlight the role of nunatak survival and species ecology in the demographic history of mountain species.

Population genomics of flat-tailed horned lizards (Phrynosoma mcallii) informs conservation and management across a fragmented Colorado Desert landscape

Phrynosoma mcallii (flat-tailed horned lizards) is a species of conservation concern in the Colorado Desert of the United States and Mexico. We analysed ddRADseq data from 45 lizards to estimate population structure, infer phylogeny, identify migration barriers, map genetic diversity hotspots, and model demography. We identified the Colorado River as the main geographic feature contributing to population structure, with the populations west of this barrier further subdivided by the Salton Sea. Phylogenetic analysis confirms that northwestern populations are nested within southeastern populations. The best-fit demographic model indicates Pleistocene divergence across the Colorado River, with significant bidirectional gene flow, and a severe Holocene population bottleneck. These patterns suggest that management strategies should focus on maintaining genetic diversity on both sides of the Colorado River and the Salton Sea. We recommend additional lands in the United States and Mexico that should be considered for similar conservation goals as those in the Rangewide Management Strategy. We also recommend periodic rangewide genomic sampling to monitor ongoing attrition of diversity, hybridization, and changing structure due to habitat fragmentation, climate change, and other long-term impacts.

Genomic evidence for evolutionary history and local adaptation of two endemic apricots: Prunus hongpingensis and P. zhengheensis

Apricot, belonging to the Armeniaca section of Rosaceae, is one of the economically important crop fruits that has been extensively cultivated. The natural wild apricots offer valuable genetic resources for crop improvement. However, some of them are endemic, with small populations, and are even at risk of extinction. In this study we unveil chromosome-level genome assemblies for two southern China endemic apricots, Prunus hongpingensis (PHP) and P. zhengheensis (PZH). We also characterize their evolutionary history and the genomic basis of their local adaptation using whole-genome resequencing data. Our findings reveal that PHP and PZH are closely related to Prunus armeniaca and form a distinct lineage. Both species experienced a decline in effective population size following the Last Glacial Maximum (LGM), which likely contributed to their current small population sizes. Despite the observed decrease in genetic diversity and heterozygosity, we do not observe an increased accumulation of deleterious mutations in these two endemic apricots. This is likely due to the combined effects of a low inbreeding coefficient and strong purifying selection. Furthermore, we identify a set of genes that have undergone positive selection and are associated with local environmental adaptation in PHP and PZH, respectively. These candidate genes can serve as valuable genetic resources for targeted breeding and improvement of cultivated apricots. Overall, our study not only enriches our comprehension of the evolutionary history of apricot species but also offers crucial insights for the conservation and future breeding of other endemic species amidst rapid climate changes.

The opposed forces of differentiation and admixture across glacial cycles in the butterfly Aglais urticae

Glacial cycles lead to periodic population interbreeding and isolation in warm-adapted species, which impact genetic structure and evolution. However, the effects of these processes on highly mobile and more cold-tolerant species are not well understood. This study aims to shed light on the phylogeographic history of Aglais urticae, a butterfly species with considerable dispersal ability, and a wide Palearctic distribution reaching the Arctic. Through the analysis of genomic data, four main genetic lineages are identified: European, Sierra Nevada, Sicily/Calabria/Peloponnese, and Eastern. The results indicate that the Sardo-Corsican endemic taxon ichnusa is a distinct species. The split between the relict lineages in southern Europe and the main European lineage is estimated to have happened 400-450 thousand years ago, with admixture observed during the Quaternary glacial cycles, and still ongoing, albeit to a much smaller extent. These results suggest that these lineages may be better treated as subspecific parapatric taxa. Ecological niche modelling supported the existence of both Mediterranean and extra-Mediterranean refugia during the glacial periods, with the main one located on the Atlantic coast. Nevertheless, gene flow between populations was possible, indicating that both differentiation and admixture have acted continuously across glacial cycles in this cold-tolerant butterfly, generally balancing each other but producing differentiated lineages in the southern peninsulas. We conclude that the population dynamics and the processes shaping the population genetic structure of cold-adapted species during the Quaternary ice ages may be different than those classically accepted for warm-adapted species.

Integration of genomic and ecological methods inform management of an undescribed, yet highly exploited, sardine species

Assessing genetic diversity within species is key for conservation strategies in the context of human-induced biotic changes. This is important in marine systems, where many species remain undescribed while being overfished, and conflicts between resource-users and conservation agencies are common. Combining niche modelling with population genomics can contribute to resolving those conflicts by identifying management units and understanding how past climatic cycles resulted in current patterns of genetic diversity. We addressed these issues on an undescribed but already overexploited species of sardine of the genus Harengula. We find that the species distribution is determined by salinity and depth, with a continuous distribution along the Brazilian mainland and two disconnected oceanic archipelagos. Genomic data indicate that such biogeographic barriers are associated with two divergent intraspecific lineages. Changes in habitat availability during the last glacial cycle led to different demographic histories among stocks. One coastal population experienced a 3.6-fold expansion, whereas an island-associated population contracted 3-fold, relative to the size of the ancestral population. Our results indicate that the island population should be managed separately from the coastal population, and that a Marine Protected Area covering part of the island population distribution can support the viability of this lineage.

Landscape Heterogeneity Explains the Genetic Differentiation of a Forest Bird across the Sino-Himalayan Mountains

Mountains are the world's most important centers of biodiversity. The Sino-Himalayan Mountains are global biodiversity hotspot due to their extremely high species richness and endemicity. Ample research investigated the impact of the Qinghai-Tibet Plateau uplift and Quaternary glaciations in driving species diversification in plants and animals across the Sino-Himalayan Mountains. However, little is known about the role of landscape heterogeneity and other environmental features in driving diversification in this region. We utilized whole genomes and phenotypic data in combination with landscape genetic approaches to investigate population structure, demography, and genetic diversity in a forest songbird species native to the Sino-Himalayan Mountains, the red-billed leiothrix (Leiothrix lutea). We identified 5 phylogeographic clades, including 1 in the East of China, 1 in Yunnan, and 3 in Tibet, roughly consistent with differences in song and plumage coloration but incongruent with traditional subspecies boundaries. Isolation-by-resistance model best explained population differentiation within L. lutea, with extensive secondary contact after allopatric isolation leading to admixture among clades. Ecological niche modeling indicated relative stability in the extent of suitable distribution areas of the species across Quaternary glacial cycles. Our results underscore the importance of mountains in the diversification of this species, given that most of the distinct genetic clades are concentrated in a relatively small area in the Sino-Himalayan Mountain region, while a single shallow clade populates vast lower-lying areas to the east. This study highlights the crucial role of landscape heterogeneity in promoting differentiation and provides a deep genomic perspective on the mechanisms through which diversity hotspots form.

Genomic Analyses Capture the Human-Induced Demographic Collapse and Recovery in a Wide-Ranging Cervid

The glacial cycles of the Quaternary heavily impacted species through successions of population contractions and expansions. Similarly, populations have been intensely shaped by human pressures such as unregulated hunting and land use changes. White-tailed and mule deer survived in different refugia through the Last Glacial Maximum, and their populations were severely reduced after the European colonization. Here, we analyzed 73 resequenced deer genomes from across their North American range to understand the consequences of climatic and anthropogenic pressures on deer demographic and adaptive history. We found strong signals of climate-induced vicariance and demographic decline; notably, multiple sequentially Markovian coalescent recovers a severe decline in mainland white-tailed deer effective population size (Ne) at the end of the Last Glacial Maximum. We found robust evidence for colonial overharvest in the form of a recent and dramatic drop in Ne in all analyzed populations. Historical census size and restocking data show a clear parallel to historical Ne estimates, and temporal Ne/Nc ratio shows patterns of conservation concern for mule deer. Signatures of selection highlight genes related to temperature, including a cold receptor previously highlighted in woolly mammoth. We also detected immune genes that we surmise reflect the changing land use patterns in North America. Our study provides a detailed picture of anthropogenic and climatic-induced decline in deer diversity and clues to understanding the conservation concerns of mule deer and the successful demographic recovery of white-tailed deer.

Computationally efficient demographic history inference from allele frequencies with supervised machine learning

Inferring past demographic history of natural populations from genomic data is of central concern in many studies across research fields. Previously, our group had developed dadi, a widely used demographic history inference method based on the allele frequency spectrum (AFS) and maximum composite likelihood optimization. However, dadi's optimization procedure can be computationally expensive. Here, we developed donni (demography optimization via neural network inference), a new inference method based on dadi that is more efficient while maintaining comparable inference accuracy. For each dadi-supported demographic model, donni simulates the expected AFS for a range of model parameters then trains a set of Mean Variance Estimation neural networks using the simulated AFS. Trained networks can then be used to instantaneously infer the model parameters from future input data AFS. We demonstrated that for many demographic models, donni can infer some parameters, such as population size changes, very well and other parameters, such as migration rates and times of demographic events, fairly well. Importantly, donni provides both parameter and confidence interval estimates from input AFS with accuracy comparable to parameters inferred by dadi's likelihood optimization while bypassing its long and computationally intensive evaluation process. donni's performance demonstrates that supervised machine learning algorithms may be a promising avenue for developing more sustainable and computationally efficient demographic history inference methods.

Unidirectional trans-Atlantic gene flow and a mixed spawning area shape the genetic connectivity of Atlantic bluefin tuna

The commercially important Atlantic bluefin tuna (Thunnus thynnus), a large migratory fish, has experienced notable recovery aided by accurate resource assessment and effective fisheries management efforts. Traditionally, this species has been perceived as consisting of eastern and western populations, spawning respectively in the Mediterranean Sea and the Gulf of Mexico, with mixing occurring throughout the Atlantic. However, recent studies have challenged this assumption by revealing weak genetic differentiation and identifying a previously unknown spawning ground in the Slope Sea used by Atlantic bluefin tuna of uncertain origin. To further understand the current and past population structure and connectivity of Atlantic bluefin tuna, we have assembled a unique dataset including thousands of genome-wide single-nucleotide polymorphisms (SNPs) from 500 larvae, young of the year and spawning adult samples covering the three spawning grounds and including individuals of other Thunnus species. Our analyses support two weakly differentiated but demographically connected ancestral populations that interbreed in the Slope Sea. Moreover, we also identified signatures of introgression from albacore (Thunnus alalunga) into the Atlantic bluefin tuna genome, exhibiting varied frequencies across spawning areas, indicating strong gene flow from the Mediterranean Sea towards the Slope Sea. We hypothesize that the observed genetic differentiation may be attributed to increased gene flow caused by a recent intensification of westward migration by the eastern population, which could have implications for the genetic diversity and conservation of western populations. Future conservation efforts should consider these findings to address potential genetic homogenization in the species.

Signals of differential introgression in the genome of natural hybrids of Caribbean anoles

Hybridization facilitates recombination between divergent genetic lineages and can be shaped by both neutral and selective processes. Upon hybridization, loci with no net fitness effects introgress randomly from parental species into the genomes of hybrid individuals. Conversely, alleles from one parental species at some loci may provide a selective advantage to hybrids, resulting in patterns of introgression that do not conform to random expectations. We investigated genomic patterns of differential introgression in natural hybrids of two species of Caribbean anoles, Anolis pulchellus and A. krugi in Puerto Rico. Hybrids exhibit A. pulchellus phenotypes but possess A. krugi mitochondrial DNA, originated from multiple, independent hybridization events, and appear to have replaced pure A. pulchellus across a large area in western Puerto Rico. Combining genome-wide SNP datasets with bioinformatic methods to identify signals of differential introgression in hybrids, we demonstrate that the genomes of hybrids are dominated by pulchellus-derived alleles and show only 10%-20% A. krugi ancestry. The majority of A. krugi loci in hybrids exhibit a signal of non-random differential introgression and include loci linked to genes involved in development and immune function. Three of these genes (delta like canonical notch ligand 1, jagged1 and notch receptor 1) affect cell differentiation and growth and interact with mitochondrial function. Our results suggest that differential non-random introgression for a subset of loci may be driven by selection favouring the inheritance of compatible mitochondrial and nuclear-encoded genes in hybrids.

An evolutionarily distinct ringed seal in the Ilulissat Icefjord

The Earth's polar regions are low rates of inter- and intraspecific diversification. An extreme mammalian example is the Arctic ringed seal (Pusa hispida hispida), which is assumed to be panmictic across its circumpolar Arctic range. Yet, local Inuit communities in Greenland and Canada recognize several regional variants; a finding supported by scientific studies of body size variation. It is however unclear whether this phenotypic variation reflects plasticity, morphs or distinct ecotypes. Here, we combine genomic, biologging and survey data, to document the existence of a unique ringed seal ecotype in the Ilulissat Icefjord (locally 'Kangia'), Greenland; a UNESCO World Heritage site, which is home to the most productive marine-terminating glacier in the Arctic. Genomic analyses reveal a divergence of Kangia ringed seals from other Arctic ringed seals about 240 kya, followed by secondary contact since the Last Glacial Maximum. Despite ongoing gene flow, multiple genomic regions appear under strong selection in Kangia ringed seals, including candidate genes associated with pelage coloration, growth and osmoregulation, potentially explaining the Kangia seal's phenotypic and behavioural uniqueness. The description of 'hidden' diversity and adaptations in yet another Arctic species merits a reassessment of the evolutionary processes that have shaped Arctic diversity and the traditional view of this region as an evolutionary freezer. Our study highlights the value of indigenous knowledge in guiding science and calls for efforts to identify distinct populations or ecotypes to understand how these might respond differently to environmental change.

Evolutionary dynamics of Euphorbia carniolica suggest a complex Plio-Pleistocene history of understorey species of deciduous forest in southeastern Europe

Deciduous forests form the dominant natural vegetation of Europe today, but were restricted to small refugia during Pleistocene cold stages, implying an evolutionary past shaped by recurrent range contractions and expansions. Cold-stage forest refugia were probably widespread in southern and central Europe, with the northwestern Balkan Peninsula being of particular importance. However, the actual number and location of deciduous forest refugia, as well as the connections between them, remain disputed. Here, we address the evolutionary dynamics of the deciduous forest understorey species Euphorbia carniolica as a proxy for past forest dynamics. To do so, we obtained genomic and morphometric data from populations representing the species' entire range, investigated phylogenetic position and intraspecific genetic variation, tested explicit demographic scenarios and applied species distribution models. Our data support two disjoint groups linked to separate refugia on the northwestern and central Balkan Peninsula. We find that genetic differentiation between groups started in the early Pleistocene via vicariance, suggesting a larger distribution in the past. Both refugia acted as sources for founder events to the southeastern Alps and the Carpathians; the latter were likely colonised before the last cold stage. In line with traditional views on the pre-Pleistocene origin of many southeastern European deciduous forest species, the origin of E. carniolica was dated to the late Pliocene. The fact that E. carniolica evolved at a time when a period of continuous forestation was ending in much of Eurasia provides an interesting biogeographical perspective on the past links between Eurasian deciduous forests and their biota.

Integrative taxonomy in a rapid speciation group associated with mating system transition: A case study in the Primula cicutariifolia complex

Accurate species delimitation is the key to biodiversity conservation and is fundamental to most branches of biology. However, species delimitation remains challenging in those evolutionary radiations associated with mating system transition from outcrossing to self-fertilization, which have frequently occurred in angiosperms and are usually accompanied by rapid speciation. Here, using the Primula cicutariifolia complex as a case, we integrated molecular, morphological and reproductive isolation evidence to test and verify whether its outcrossing (distylous) and selfing (homostylous) populations have developed into independent evolutionary lineages. Phylogenetic trees based on whole plastomes and SNPs of the nuclear genome both indicated that the distylous and homostylous populations grouped into two different clades. Multispecies coalescent, gene flow and genetic structure analyses all supported such two clades as two different genetic entities. In morphology, as expected changes in selfing syndrome, homostylous populations have significantly fewer umbel layers and smaller flower and leaf sizes compared to distylous populations, and the variation range of some floral traits, such as corolla diameter and umbel layers, show obvious discontinuity. Furthermore, hand-pollinated hybridization between the two clades produced almost no seeds, indicating that well post-pollination reproductive isolation has been established between them. Therefore, the distylous and homostylous populations in this studied complex are two independent evolutionary lineages, and thus these distylous populations should be treated as a distinct species, here named Primula qiandaoensis W. Zhang & J.W. Shao sp. nov.. Our empirical study of the P. cicutariifolia complex highlights the importance of applying multiple lines of evidence, in particular genomic data, to delimit species in pervasive evolutionary plant radiations associated with mating system transition.

Means, motive and opportunity for biological invasions: Genetic introgression in a fungal pathogen

Invasions by fungal plant pathogens pose a significant threat to the health of agricultural ecosystems. Despite limited standing genetic variation, many invasive fungal species can adapt and spread rapidly, resulting in significant losses to crop yields. Here, we report on the population genomics of Colletotrichum truncatum, a polyphagous pathogen that can infect more than 460 plant species, and an invasive pathogen of soybean in Brazil. We study the whole-genome sequences of 18 isolates representing 10 fields from two major regions of soybean production. We show that Brazilian C. truncatum is subdivided into three phylogenetically distinct lineages that exchange genetic variation through hybridization. Introgression affects 2%-30% of the nucleotides of genomes and varies widely between the lineages. We find that introgressed regions comprise secreted protein-encoding genes, suggesting possible co-evolutionary targets for selection in those regions. We highlight the inherent vulnerability of genetically uniform crops in the agro-ecological environment, particularly when faced with pathogens that can take full advantage of the opportunities offered by an increasingly globalized world. Finally, we discuss "the means, motive and opportunity" of fungal pathogens and how they can become invasive species of crops. We call for more population genomic studies because such analyses can help identify geographical areas and pathogens that pose a risk, thereby helping to inform control strategies to better protect crops in the future.

Evidence for ecotone speciation across an African rainforest-savanna gradient

Accelerating climate change and habitat loss make it imperative that plans to conserve biodiversity consider species' ability to adapt to changing environments. However, in biomes where biodiversity is highest, the evolutionary mechanisms responsible for generating adaptative variation and, ultimately, new species are frequently poorly understood. African rainforests represent one such biome, as decadal debates continue concerning the mechanisms generating African rainforest biodiversity. These debates hinge on the relative importance of geographic isolation versus divergent natural selection across environmental gradients. Hindering progress is a lack of robust tests of these competing hypotheses. Because African rainforests are severely at-risk due to climate change and other anthropogenic activities, addressing this long-standing debate is critical for making informed conservation decisions. We use demographic inference and allele frequency-environment relationships to investigate mechanisms of diversification in an African rainforest skink, Trachylepis affinis, a species inhabiting the gradient between rainforest and rainforest-savanna mosaic (ecotone). We provide compelling evidence of ecotone speciation, in which gene flow has all but ceased between rainforest and ecotone populations, at a level consistent with infrequent hybridization between sister species. Parallel patterns of genomic, morphological, and physiological divergence across this environmental gradient and pronounced allele frequency-environment correlation indicate speciation is mostly probably driven by ecological divergence, supporting a central role for divergent natural selection. Our results provide strong evidence for the importance of ecological gradients in African rainforest speciation and inform conservation strategies that preserve the processes that produce and maintain biodiversity.

Deep genome-wide divergences among species in White Cloud Mountain minnow Tanichthys albonubes (Cypriniformes: Tanichthyidae) complex: Conservation and species management implications

Identifying cryptic species is important for the assessments of biodiversity. Further, untangling mechanisms underlying the origins of cryptic species can facilitate our understanding of evolutionary processes. Advancements in genomic approaches for non-model systems have offered unprecedented opportunities to investigate these areas. The White Cloud Mountain minnow (Tanichthys albonubes) is a popular freshwater pet fish worldwide but its wild populations in China are critically endangered. Recent research based on a few molecular markers suggested that this species in fact comprised seven cryptic species, of which six were previously unknown. Here, we tested six of these cryptic species and quantified genomic interspecific divergences between species in the T. albonubes complex by analyzing genome-wide restriction site-associated DNA sequencing (RADseq) data generated from 189 individuals sampled from seven populations (including an outgroup congeneric species, T. micagemmae). We found that six cryptic species previously suggested were well supported by RADseq data. The genetic diversity of each species in the T. albonubes complex was low compared with T. micagemmae and the contemporary effective population sizes (Ne) of each cryptic species were small. Phylogenetic analysis showed seven clades with high support values confirmed with Neighbor-Net trees. The pairwise divergences between species in T. albonubes complex were deep and the highly differentiated loci were evenly distributed across the genome. We proposed that the divergence level of T. albonubes complex is at a late stage of cryptic speciation and lacking gene flow. Our findings provide new insights into cryptic speciation and have important implications for conservation and species management of T. albonubes complex.

Rapid diversification underlying the global dominance of a cosmopolitan phytoplankton

Marine phytoplankton play important roles in the global ecosystem, with a limited number of cosmopolitan keystone species driving their biomass. Recent studies have revealed that many of these phytoplankton are complexes composed of sibling species, but little is known about the evolutionary processes underlying their formation. Gephyrocapsa huxleyi, a widely distributed and abundant unicellular marine planktonic algae, produces calcified scales (coccoliths), thereby significantly affects global biogeochemical cycles via sequestration of inorganic carbon. This species is composed of morphotypes defined by differing degrees of coccolith calcification, the evolutionary ecology of which remains unclear. Here, we report an integrated morphological, ecological and genomic survey across globally distributed G. huxleyi strains to reconstruct evolutionary relationships between morphotypes in relation to their habitats. While G. huxleyi has been considered a single cosmopolitan species, our analyses demonstrate that it has evolved to comprise at least three distinct species, which led us to formally revise the taxonomy of the G. huxleyi complex. Moreover, the first speciation event occurred before the onset of the last interglacial period (~140 ka), while the second followed during this interglacial. Then, further rapid diversifications occurred during the most recent ice-sheet expansion of the last glacial period and established morphotypes as dominant populations across environmental clines. These results suggest that glacial-cycle dynamics contributed to the isolation of ocean basins and the segregations of oceans fronts as extrinsic drivers of micro-evolutionary radiations in extant marine phytoplankton.

Phylogenomic inference and demographic model selection suggest peripatric separation of the cryptic steppe ant species Plagiolepis pyrenaica stat. rev

The ant Plagiolepis taurica Santschi, 1920 (Hymenoptera, Formicidae) is a typical species of the Eurasian steppes, a large grassland dominated biome that stretches continuously from Central Asia to Eastern Europe and is represented by disjunct outposts also in Central and Western Europe. The extent of this biome has been influenced by the Pleistocene climate, and steppes expanded recurrently during cold stages and contracted in warm stages. Consequently, stenotopic steppe species such as P. taurica repeatedly went through periods of demographic expansion and severe isolation. Here, we explore the impact of these dynamics on the genetic diversification within P. taurica. Delimitation of P. taurica from other Plagiolepis species has been unclear since its initial description, which raised questions on both its classification and its spatiotemporal diversification early on. We re-evaluate species limits and explore underlying mechanisms driving speciation by using an integrative approach based on genomic and morphometric data. We found large intraspecific divergence within P. taurica and resolved geographically coherent western and eastern genetic groups, which likewise differed morphologically. A morphometric survey of type material showed that Plagiolepis from the western group were more similar to P. barbara pyrenaica Emery, 1921 than to P. taurica; we thus lift the former from synonymy and establish it as separate species, P. pyrenaica stat. rev. Explicit evolutionary model testing based on genomic data supported a peripatric speciation for the species pair, probably as a consequence of steppe contraction and isolation during the mid-Pleistocene. We speculate that this scenario could be exemplary for many stenotopic steppe species, given the emphasized dynamics of Eurasian steppes.

Climate-induced range shifts drive adaptive response via spatio-temporal sieving of alleles

Quaternary climate fluctuations drove many species to shift their geographic ranges, in turn shaping their genetic structures. Recently, it has been argued that adaptation may have accompanied species range shifts via the "sieving" of genotypes during colonisation and establishment. However, this has not been directly demonstrated, and knowledge remains limited on how different evolutionary forces, which are typically investigated separately, interacted to jointly mediate species responses to past climatic change. Here, through whole-genome re-sequencing of over 1200 individuals of the carnation Dianthus sylvestris coupled with integrated population genomic and gene-environment models, we reconstruct the past neutral and adaptive landscape of this species as it was shaped by the Quaternary glacial cycles. We show that adaptive responses emerged concomitantly with the post-glacial range shifts and expansions of this species in the last 20 thousand years. This was due to the heterogenous sieving of adaptive alleles across space and time, as populations expanded out of restrictive glacial refugia into the broader and more heterogeneous range of habitats available in the present-day inter-glacial. Our findings reveal a tightly-linked interplay of migration and adaptation under past climate-induced range shifts, which we show is key to understanding the spatial patterns of adaptive variation we see in species today.

Whole-genome resequencing reveals complex effects of geographical-palaeoclimatic interactions on diversification of moustache toads in East Asia

Geographical features and palaeoclimatic fluctuations are two classical evolutionary forces that shape genetic diversification within species. Fine-grained analysis of the mechanisms involved through population demographic processes, however, remains limited. Taking advantage of two recently published reference genomes, we resequenced the genomes and examined the evolutionary history of the moustache toads, a group endemic to East Asia where complex topography and fluctuating palaeoclimate are known to have had profound impacts on organisms. Moustache toads probably originated in southeast Yunnan, China, and diversified towards the northwestern of Yunnan, as well as central and eastern China. Further exploration based on three widespread species (Leptobrachium ailaonicum, L. boringii and L. liui) using demographic modelling and species distribution models revealed that mountains and river valleys in East Asia not only functioned as geographical barriers, but also provided dispersal corridors and facilitated continuous migration or post-glacial secondary contact among moustache toad populations. Furthermore, periodic oscillation of effective population sizes accompanying fluctuations of historical temperature and population contraction at the Last Glacial Maximum support the widespread impact of climatic changes of the Pleistocene on species diversification in East Asia. This impact was moderate for populations of L. ailaonicum and L. boringii in the southwestern mountains but severe for populations of L. liui in the eastern lowland regions of continental East Asia, which is supported by different degrees of change of their effective population sizes. Our findings reveal mechanisms underlying genetic diversification among moustache toads, and highlight the power of genomic data and demographic modelling for examining complex historical population-level processes and for understanding how geographical and palaeoclimatic factors interactively shape current intraspecific diversity.

Cryptic diversity and spatial genetic variation in the coral Acropora tenuis and its endosymbionts across the Great Barrier Reef

Genomic studies are uncovering extensive cryptic diversity within reef-building corals, suggesting that evolutionarily and ecologically relevant diversity is highly underestimated in the very organisms that structure coral reefs. Furthermore, endosymbiotic algae within coral host species can confer adaptive responses to environmental stress and may represent additional axes of coral genetic variation that are not constrained by taxonomic divergence of the cnidarian host. Here, we examine genetic variation in a common and widespread, reef-building coral, Acropora tenuis, and its associated endosymbiotic algae along the entire expanse of the Great Barrier Reef (GBR). We use SNPs derived from genome-wide sequencing to characterize the cnidarian coral host and organelles from zooxanthellate endosymbionts (genus Cladocopium). We discover three distinct and sympatric genetic clusters of coral hosts, whose distributions appear associated with latitude and inshore-offshore reef position. Demographic modelling suggests that the divergence history of the three distinct host taxa ranges from 0.5 to 1.5 million years ago, preceding the GBR's formation, and has been characterized by low-to-moderate ongoing inter-taxon gene flow, consistent with occasional hybridization and introgression typifying coral evolution. Despite this differentiation in the cnidarian host, A. tenuis taxa share a common symbiont pool, dominated by the genus Cladocopium (Clade C). Cladocopium plastid diversity is not strongly associated with host identity but varies with reef location relative to shore: inshore colonies contain lower symbiont diversity on average but have greater differences between colonies as compared with symbiont communities from offshore colonies. Spatial genetic patterns of symbiont communities could reflect local selective pressures maintaining coral holobiont differentiation across an inshore-offshore environmental gradient. The strong influence of environment (but not host identity) on symbiont community composition supports the notion that symbiont community composition responds to habitat and may assist in the adaptation of corals to future environmental change.

Testing assertions of widespread introgressive hybridization in a clade of neotropical toads with low mate selectivity (Rhinella granulosa species group)

Discordance between different genomic regions, often identified through multilocus sequencing of selected markers, presents particular difficulties in identifying historical processes which drive species diversity and boundaries. Mechanisms causing discordance, such as incomplete lineage sorting or introgression due to interspecific hybridization, are better identified based on population-level genomic datasets. In the toads of the Rhinella granulosa species group, patterns of mito-nuclear discordance and potential hybridization have been reported by several studies. However, these patterns were proposed based on few loci, such that alternative mechanisms behind gene-tree heterogeneity cannot be ruled out. Using genome-wide ddRADseq loci from a subset of species within this clade, we found only partial concordance between currently recognized species-level taxon boundaries and patterns of genetic structure. While most taxa within the R. granulosa group correspond to clades, genetic clustering analyses sometimes grouped distinct taxonomic units into a single cluster. Moreover, levels of admixture between inferred clusters were limited and restricted to a single taxon pair which is best explained by incomplete lineage sorting as opposed to introgressive hybridization, according to D-statistics results. These findings contradict previous assertions of widespread cryptic diversity and gene flow within the R. granulosa clade. Lastly, our analyses suggest that diversification events within the Rhinella granulosa group mostly dated back to the early Pliocene, being generally younger than species divergences in other closely related clades that present high levels of cross-species gene flow. This finding uniquely contradicts common assertions that this young clade of toads exhibits interspecific hybridization.

Diversification and historical demography of Rhampholeon spectrum in West-Central Africa

Pygmy Chameleons of the genus Rhampholeon represent a moderately diverse, geographically circumscribed radiation, with most species (18 out of 19 extant taxa) limited to East Africa. The one exception is Rhampholeon spectrum, a species restricted to West-Central African rainforests. We set out to characterize the geographic basis of genetic variation in this disjunctly distributed Rhampholeon species using a combination of multilocus Sanger data and genomic sequences to explore population structure and range-wide phylogeographic patterns. We also employed demographic analyses and niche modeling to distinguish between alternate explanations to contextualize the impact of past geological and climatic events on the present-day distribution of intraspecific genetic variation. Phylogenetic analyses suggest that R. spectrum is a complex of five geographically delimited populations grouped into two major clades (montane vs. lowland). We found pronounced population structure suggesting that divergence and, potentially, speciation began between the late Miocene and the Pleistocene. Sea level changes during the Pleistocene climatic oscillations resulted in allopatric divergence associated with dispersal over an ocean channel barrier and colonization of Bioko Island. Demographic inferences and range stability mapping each support diversification models with secondary contact due to population contraction in lowland and montane refugia during the interglacial period. Allopatric divergence, congruent with isolation caused by geologic uplift of the East African rift system, the "descent into the Icehouse," and aridification of sub-Saharan Africa during the Eocene-Oligocene are identified as the key events explaining the population divergence between R. spectrum and its closely related sister clade from the Eastern Arc Mountains. Our results unveil cryptic genetic diversity in R. spectrum, suggesting the possibility of a species complex distributed across the Lower Guinean Forest and the Island of Bioko. We highlight the major element of species diversification that modelled today's diversity and distributions in most West-Central African vertebrates.

Estimation of site frequency spectra from low-coverage sequencing data using stochastic EM reduces overfitting, runtime, and memory usage

The site frequency spectrum is an important summary statistic in population genetics used for inference on demographic history and selection. However, estimation of the site frequency spectrum from called genotypes introduces bias when working with low-coverage sequencing data. Methods exist for addressing this issue but sometimes suffer from 2 problems. First, they can have very high computational demands, to the point that it may not be possible to run estimation for genome-scale data. Second, existing methods are prone to overfitting, especially for multidimensional site frequency spectrum estimation. In this article, we present a stochastic expectation-maximization algorithm for inferring the site frequency spectrum from NGS data that address these challenges. We show that this algorithm greatly reduces runtime and enables estimation with constant, trivial RAM usage. Furthermore, the algorithm reduces overfitting and thereby improves downstream inference. An implementation is available at github.com/malthesr/winsfs.

Genome-wide SNPs in the spiny lobster Panulirus homarus reveal a hybrid origin for its subspecies

Background

Evolutionary divergence and speciation often occur at a slower rate in the marine realm due to the higher potential for long-distance reproductive interaction through larval dispersal. One common evolutionary pattern in the Indo-Pacific, is divergence of populations and species at the peripheries of widely-distributed organisms. However, the evolutionary and demographic histories of such divergence are yet to be well understood. Here we address these issues by coupling genome-wide SNP data with mitochondrial DNA sequences to test the patterns of genetic divergence and possible secondary contact among geographically distant populations of the highly valuable spiny lobster Panulirus homarus species complex, distributed widely through the Indo-Pacific, from South Africa to the Marquesas Islands.

Result

After stringent filtering, 2020 SNPs were used for population genetic and demographic analyses, revealing strong regional structure (FST = 0.148, P < 0001), superficially in accordance with previous analyses. However, detailed demographic analyses supported a much more complex evolutionary history of these populations, including a hybrid origin of a North-West Indian Ocean (NWIO) population, which has previously been discriminated morphologically, but not genetically. The best-supported demographic models suggested that the current genetic relationships among populations were due to a complex series of past divergences followed by asymmetric migration in more recent times.

Conclusion

Overall, this study suggests that alternating periods of marine divergence and gene flow have driven the current genetic patterns observed in this lobster and may help explain the observed wider patterns of marine species diversity in the Indo-Pacific.

Pre-domestication bottlenecks of the cultivated seaweed Gracilaria chilensis

Gracilaria chilensis is the main cultivated seaweed in Chile. The low genetic diversity observed in the Chilean populations has been associated with the over-exploitation of natural beds and/or the founder effect that occurred during post-glacial colonization from New Zealand. How these processes have affected its evolutionary trajectory before farming and incipient domestication is poorly understood. In this study, we used 2232 single nucleotide polymorphisms (SNPs) to assess how the species' evolutionary history in New Zealand (its region of origin), the founder effect linked to transoceanic dispersion and colonization of South America, and the recent over-exploitation of natural populations have influenced the genetic architecture of G. chilensis in Chile. The contrasting patterns of genetic diversity and structure observed between the two main islands in New Zealand attest to the important effects of Quaternary glacial cycles on G. chilensis. Approximate Bayesian Computation (ABC) analyses indicated that Chatham Island and South America were colonized independently near the end of the Last Glacial Maximum and emphasized the importance of coastal and oceanic currents during that period. Furthermore, ABC analyses inferred the existence of a recent and strong genetic bottleneck in Chile, matching the period of over-exploitation of the natural beds during the 1970s, followed by rapid demographic expansion linked to active clonal propagation used in farming. Recurrent genetic bottlenecks strongly eroded the genetic diversity of G. chilensis prior to its cultivation, raising important challenges for the management of genetic resources in this incipiently domesticated species.

Genomic divergence, local adaptation, and complex demographic history may inform management of a popular sportfish species complex

The Neosho Bass (Micropterus velox), a former subspecies of the keystone top-predator and globally popular Smallmouth Bass (M. dolomieu), is endemic and narrowly restricted to small, clear streams of the Arkansas River Basin in the Central Interior Highlands (CIH) ecoregion, USA. Previous studies have detected some morphological, genetic, and genomic differentiation between the Neosho and Smallmouth Basses; however, the extent of neutral and adaptive divergence and patterns of intraspecific diversity are poorly understood. Furthermore, lineage diversification has likely been impacted by gene flow in some Neosho populations, which may be due to a combination of natural biogeographic processes and anthropogenic introductions. We assessed: (1) lineage divergence, (2) local directional selection (adaptive divergence), and (3) demographic history among Smallmouth Bass populations in the CIH using population genomic analyses of 50,828 single-nucleotide polymorphisms (SNPs) obtained through ddRAD-seq. Neosho and Smallmouth Bass formed monophyletic clades with 100% bootstrap support. We identified two major lineages within each species. We discovered six Neosho Bass populations (two nonadmixed and four admixed) and three nonadmixed Smallmouth Bass populations. We detected 29 SNPs putatively under directional selection in the Neosho range, suggesting populations may be locally adapted. Two populations were admixed via recent asymmetric secondary contact, perhaps after anthropogenic introduction. Two other populations were likely admixed via combinations of ancient and recent processes. These species comprise independently evolving lineages, some having experienced historical and natural admixture. These results may be critical for management of Neosho Bass as a distinct species and may aid in the conservation of other species with complex biogeographic histories.