Integrating cryptic diversity into coral evolution, symbiosis and conservation

Understanding how diversity evolves and is maintained is critical to predicting the future trajectories of ecosystems under climate change; however, our understanding of these processes is limited in marine systems. Corals, which engineer reef ecosystems, are critically threatened by climate change, and global efforts are underway to conserve and restore populations as attempts to mitigate ocean warming continue. Recently, sequencing efforts have uncovered widespread undescribed coral diversity, including 'cryptic lineages'-genetically distinct but morphologically similar coral taxa. Such cryptic lineages have been identified in at least 24 coral genera spanning the anthozoan phylogeny and across ocean basins. These cryptic lineages co-occur in many reef systems, but their distributions often differ among habitats. Research suggests that cryptic lineages are ecologically specialized and several examples demonstrate differences in thermal tolerance, highlighting the critical implications of this diversity for predicting coral responses to future warming. Here, we draw attention to recent discoveries, discuss how cryptic diversity affects the study of coral adaptation and acclimation to future environments, explore how it shapes symbiotic partnerships, and highlight challenges and opportunities for conservation and restoration efforts.

The genomic diversity of the Eliurus genus in northern Madagascar with a putative new species

Madagascar exhibits extraordinarily high level of species richness and endemism, while being severely threatened by habitat loss and fragmentation (HL&F). In front of these threats to biodiversity, conservation effort can be directed, for instance, in the documentation of species that are still unknown to science, or in investigating how species respond to HL&F. The tufted-tail rats genus (Eliurus spp.) is the most speciose genus of endemic rodents in Madagascar, with 13 described species, which occupy two major habitat types: dry or humid forests. The large species diversity and association to specific habitat types make the Eliurus genus a suitable model for investigating species adaptation to new environments, as well as response to HL&F (dry vs humid). In the present study, we investigated Eliurus spp. genomic diversity across northern Madagascar, a region covered by both dry and humid fragmented forests. From the mitochondrial DNA (mtDNA) and nuclear genomic (RAD-seq) data of 124 Eliurus individuals sampled in poorly studied forests of northern Madagascar, we identified an undescribed Eliurus taxon (Eliurus sp. nova). We tested the hypothesis of a new Eliurus species using several approaches: i) DNA barcoding; ii) phylogenetic inferences; iii) species delimitation tests based on the Multi-Species Coalescent (MSC) model, iv) genealogical divergence index (gdi); v) an ad-hoc test of isolation-by-distance within versus between sister-taxa, vi) comparisons of %GC content patterns and vii) morphological analyses. All analyses support the recognition of the undescribed lineage as a putative distinct species. In addition, we show that Eliurus myoxinus, a species known from the dry forests of western Madagascar, is, surprisingly, found mostly in humid forests in northern Madagascar. In conclusion, we discuss the implications of such findings in the context of Eliurus species evolution and diversification, and use the distribution of northern Eliurus species as a proxy for reconstructing past changes in forest cover and vegetation type in northern Madagascar.

Abundant incongruence in a clade endemic to a biodiversity hotspot: Phylogenetics of the scrub mint clade (Lamiaceae)

The Scrub Mint clade(Lamiaceae) provides a unique system for investigating the evolutionary processes driving diversification in the North American Coastal Plain from both a systematic and biogeographic context. The clade comprisesDicerandra, Conradina, Piloblephis, Stachydeoma, and four species of the broadly defined genus Clinopodium(Mentheae; Lamiaceae), almost all of which are endemic to the North American Eastern Coastal Plain. Most species of this clade are threatened or endangered and restricted to sandhill or a mosaic of scrub habitats. We analyzed relationships in this clade to understand the evolution of the group and identify evolutionary mechanisms acting on the clade, with important implications for conservation. We used a target-capture method to sequence and analyze 238 nuclear loci across all species of scrub mints, reconstructed the phylogeny, and calculated gene tree concordance, gene tree estimation error, and reticulation indices for every node in the tree using ML methods. Phylogenetic networks were used to determine reticulation events. Our nuclear phylogenetic estimates were consistent with previous results, while greatly increasing the robustness of taxon sampling. The phylogeny resolved the full relationship between Dicerandra and Conradina and the less-studied members of the clade (Piloblephis, Stachydeoma, Clinopodium spp.). We found hotspots of gene tree discordance and reticulation throughout the tree, especially in perennial Dicerandra. Several instances of reticulation events were uncovered between annual and perennial Dicerandra, and within the Conradina + allies clade. Incomplete lineage sorting also likely contributed to phylogenetic discordance. These results clarify phylogenetic relationships in the clade and provide insight on important evolutionary drivers in the clade, such as hybridization. General relationships in the group were confirmed, while the large amount of gene tree discordance is likely due to reticulation across the phylogeny.

Genetically Encoded Lizard Color Divergence for Camouflage and Thermoregulation

Local adaptation is critical in speciation and evolution, yet comprehensive studies on proximate and ultimate causes of local adaptation are generally scarce. Here, we integrated field ecological experiments, genome sequencing, and genetic verification to demonstrate both driving forces and molecular mechanisms governing local adaptation of body coloration in a lizard from the Qinghai-Tibet Plateau. We found dark lizards from the cold meadow population had lower spectrum reflectance but higher melanin contents than light counterparts from the warm dune population. Additionally, the colorations of both dark and light lizards facilitated the camouflage and thermoregulation in their respective microhabitat simultaneously. More importantly, by genome resequencing analysis, we detected a novel mutation in Tyrp1 that underpinned this color adaptation. The allele frequencies at the site of SNP 459# in the gene of Tyrp1 are 22.22% G/C and 77.78% C/C in dark lizards and 100% G/G in light lizards. Model-predicted structure and catalytic activity showed that this mutation increased structure flexibility and catalytic activity in enzyme TYRP1, and thereby facilitated the generation of eumelanin in dark lizards. The function of the mutation in Tyrp1 was further verified by more melanin contents and darker coloration detected in the zebrafish injected with the genotype of Tyrp1 from dark lizards. Therefore, our study demonstrates that a novel mutation of a major melanin-generating gene underpins skin color variation co-selected by camouflage and thermoregulation in a lizard. The resulting strong selection may reinforce adaptive genetic divergence and enable the persistence of adjacent populations with distinct body coloration.

Genetics and the Evolution of Prezygotic Isolation

The significance of prezygotic isolation for speciation has been recognized at least since the Modern Synthesis. However, fundamental questions remain. For example, how are genetic associations between traits that contribute to prezygotic isolation maintained? What is the source of genetic variation underlying the evolution of these traits? And how do prezygotic barriers affect patterns of gene flow? We address these questions by reviewing genetic features shared across plants and animals that influence prezygotic isolation. Emerging technologies increasingly enable the identification and functional characterization of the genes involved, allowing us to test established theoretical expectations. Embedding these genes in their developmental context will allow further predictions about what constrains the evolution of prezygotic isolation. Ongoing improvements in statistical and computational tools will reveal how pre- and postzygotic isolation may differ in how they influence gene flow across the genome. Finally, we highlight opportunities for progress by combining theory with appropriate data.

Hybridization in late stages of speciation: Strong but incomplete genome-wide reproductive isolation and 'large Z-effect' in a moving hybrid zone

In organisms reproducing sexually, speciation occurs when increasing divergence results in pre- or post-zygotic reproductive isolation between lineages. Studies focusing on reproductive isolation origin in early stages of speciation are common and many rely on genomic scans to infer introgression providing limited information on the genomic architecture of reproductive isolation long-term maintenance. This study analyses a natural hybrid zone between two species in a late stage of speciation. We used ddRADseq genotyping in the contact between Podarcis bocagei and P. carbonelli to examine admixture extent, analyse hybrid zone stability and assess genome-wide variation in selection against introgression. We confirmed strong but incomplete reproductive isolation in a bimodal hybrid zone. New findings revealed population genetic structure within P. carbonelli in the contact zone; geographical and genomic clines analysis suggested strong selection against gene flow, but a relatively small proportion of the loci can introgress, mostly within the narrow contact zone. However, geographical clines revealed that a few introgressed loci show signs of potential positive selection, particularly into P. bocagei. Geographical clines also detected a signal of hybrid zone movement towards P. bocagei distribution. Genomic cline analysis revealed heterogeneous patterns of introgression among loci within the syntopy zone, but the majority maintain a strong association with the genomic background of origin. However, incongruences between both cline approaches were found, potentially driven by confounding effects on genomic clines. Last, an important role of the Z chromosome in reproductive isolation is suggested. Importantly, overall patterns of restricted introgression seem to result from numerous strong intrinsic barriers across the genome.

Characteristics of plastid genomes in the genus Ceratostigma inhabiting arid habitats in China and their phylogenomic implications

BACKGROUND

Ceratostigma, a genus in the Plumbaginaceae, is an ecologically dominant group of shrubs, subshrub and herb mainly distributed in Qinghai-Tibet Plateau and North China. Ceratostigma has been the focal group in several studies, owing to their importance in economic and ecological value and unique breeding styles. Despite this, the genome information is limited and interspecific relationships within the genus Cerotastigma remains unexplored. Here we sequenced, assembled and characterized the 14 plastomes of five species, and conducted phylogenetic analyses of Cerotastigma using plastomes and nuclear ribosomal DNA (nrDNA) data.

RESULTS

Fourteen Cerotastigma plastomes possess typical quadripartite structures with lengths from 164,076 to 168,355 bp that consist of a large single copy, a small single copy and a pair of inverted repeats, and contain 127-128 genes, including 82-83 protein coding genes, 37 transfer RNAs and eight ribosomal RNAs. All plastomes are highly conservative and similar in gene order, simple sequence repeats (SSRs), long repeat repeats and codon usage patterns, but some structural variations in the border of single copy and inverted repeats. Mutation hotspots in coding (Pi values > 0.01: matK, ycf3, rps11, rps3, rpl22 and ndhF) and non-coding regions (Pi values > 0.02: trnH-psbA, rps16-trnQ, ndhF-rpl32 and rpl32-trnL) were identified among plastid genomes that could be served as potential molecular markers for species delimitation and genetic variation studies in Cerotastigma. Gene selective pressure analysis showed that most protein-coding genes have been under purifying selection except two genes. Phylogenetic analyses based on whole plastomes and nrDNA strongly support that the five species formed a monophyletic clade. Moreover, interspecific delimitation was well resolved except C. minus, individuals of which clustered into two main clades corresponding to their geographic distributions. The topology inferred from the nrDNA dataset was not congruent with the tree derived from the analyses of the plastid dataset.

CONCLUSION

These findings represent the first important step in elucidating plastome evolution in this widespread distribution genus Cerotastigma in the Qinghai-Tibet Plateau. The detailed information could provide a valuable resource for understanding the molecular dynamics and phylogenetic relationship in the family Plumbaginaceae. Lineage genetic divergence within C. minus was perhaps promoted by geographic barriers in the Himalaya and Hengduan Mountains region, but introgression or hybridization could not be completely excluded.

REFMAKER: make your own reference to target nuclear loci in low coverage genome skimming libraries. Phylogenomic application in Sapotaceae

Genome skimming approach is widely used in plant systematics to infer phylogenies mostly from organelle genomes. However, organelles represent only 10% of the produced libraries, and the low coverage associated with these libraries (< 3X) prevents the capture of nuclear sequences, which are not always available in non-model organisms or limited to the ribosomal regions. We developed REFMAKER, a user-friendly pipeline, to create specific sets of nuclear loci that can next be extracted directly from the genome skimming libraries. For this, a catalogue is built from the meta-assembly of each library contigs and cleaned by selecting the nuclear regions and removing duplicates from clustering steps. Libraries are next mapped onto this catalogue and consensus sequences are generated to produce a ready-to-use phylogenetic matrix following different filtering parameters aiming at removing putative errors and paralogous sequences. REFMAKER allowed us to infer a well resolved phylogeny in Capurodendron (Sapotaceae) on 67 nuclear loci from low-coverage libraries (<1X). The resulting phylogeny is concomitant with one previously inferred on 638 nuclear genes from target enrichment libraries. While it remains preliminary because of this low sequencing depth, REFMAKER therefore opens perspectives in phylogenomics by allowing nuclear phylogeny reconstructions with genome skimming datasets.

Completing the TRB family: newly characterized members show ancient evolutionary origins and distinct localization, yet similar interactions

Telomere repeat binding proteins (TRBs) belong to a family of proteins possessing a Myb-like domain which binds to telomeric repeats. Three members of this family (TRB1, TRB2, TRB3) from Arabidopsis thaliana have already been described as associated with terminal telomeric repeats (telomeres) or short interstitial telomeric repeats in gene promoters (telo-boxes). They are also known to interact with several protein complexes: telomerase, Polycomb repressive complex 2 (PRC2) E(z) subunits and the PEAT complex (PWOs-EPCRs-ARIDs-TRBs). Here we characterize two novel members of the TRB family (TRB4 and TRB5). Our wide phylogenetic analyses have shown that TRB proteins evolved in the plant kingdom after the transition to a terrestrial habitat in Streptophyta, and consequently TRBs diversified in seed plants. TRB4-5 share common TRB motifs while differing in several others and seem to have an earlier phylogenetic origin than TRB1-3. Their common Myb-like domains bind long arrays of telomeric repeats in vitro, and we have determined the minimal recognition motif of all TRBs as one telo-box. Our data indicate that despite the distinct localization patterns of TRB1-3 and TRB4-5 in situ, all members of TRB family mutually interact and also bind to telomerase/PRC2/PEAT complexes. Additionally, we have detected novel interactions between TRB4-5 and EMF2 and VRN2, which are Su(z)12 subunits of PRC2.

Origin and diversification of a Himalayan orchid genus Pleione

Pleione is an orchid endemically distributed in high mountain areas across the Hengduan Mountains (HDM), Himalayas, Southeast Asia and South of China. The unique flower shapes, rich colors and immense medicinal importance of Pleione are valuable ornamental and economic resources. However, the phylogenetic relationships and evolutionary history of the genus have not yet been comprehensively resolved. Here, the evolutionary history of Pleione was investigated using single-copy gene single nucleotide polymorphisms and chloroplast genome datasets. The data revealed that Pleione could be divided into five clades. Discordance in topology between the two phylogenetic trees and network and D-statistic analyses indicated the occurrence of reticulate evolution in the genus. The evolution could be attributed to introgression and incomplete lineage sorting. Ancestral area reconstruction suggested that Pleione was originated from the HDM. Uplifting of the HDM drove rapid diversification by creating conditions favoring rapid speciation. This coincided with two periods of consolidation of the Asian monsoon climate, which caused the first rapid diversification of Pleione from 8.87 to 7.83 Mya, and a second rapid diversification started at around 4.05 Mya to Pleistocene. The interaction between Pleione and climate changes, especially the monsoons, led to the current distribution pattern and shaped the dormancy characteristic of the different clades. In addition to revealing the evolutionary relationship of Pleione with orogeny and climate changes, the findings of this study provide insights into the speciation and diversification mechanisms of plants in the East Asian flora.

Genome architecture and diverged selection shaping pattern of genomic differentiation in wild barley

Divergent selection of populations in contrasting environments leads to functional genomic divergence. However, the genomic architecture underlying heterogeneous genomic differentiation remains poorly understood. Here, we de novo assembled two high-quality wild barley (Hordeum spontaneum K. Koch) genomes and examined genomic differentiation and gene expression patterns under abiotic stress in two populations. These two populations had a shared ancestry and originated in close geographic proximity but experienced different selective pressures due to their contrasting micro-environments. We identified structural variants that may have played significant roles in affecting genes potentially associated with well-differentiated phenotypes such as flowering time and drought response between two wild barley genomes. Among them, a 29-bp insertion into the promoter region formed a cis-regulatory element in the HvWRKY45 gene, which may contribute to enhanced tolerance to drought. A single SNP mutation in the promoter region may influence HvCO5 expression and be putatively linked to local flowering time adaptation. We also revealed significant genomic differentiation between the two populations with ongoing gene flow. Our results indicate that SNPs and small SVs link to genetic differentiation at the gene level through local adaptation and are maintained through divergent selection. In contrast, large chromosome inversions may have shaped the heterogeneous pattern of genomic differentiation along the chromosomes by suppressing chromosome recombination and gene flow. Our research offers novel insights into the genomic basis underlying local adaptation and provides valuable resources for the genetic improvement of cultivated barley.

Climate influences the genetic structure and niche differentiation among populations of the olive field mouse Abrothrix olivacea (Cricetidae: Abrotrichini)

Even when environmental variation over time and space is commonly considered as an important driver of population divergence, few evaluations of intraspecific genetic variation explicitly assess whether observed structure has been caused by or is correlated with landscape heterogeneity. Several phylogeographic studies have characterized the mitochondrial diversity of Abrothrix olivacea, but none has incorporated landscape genetics analyses and ecological niche modeling, leaving a gap in the understanding of the species evolutionary history. Here, these aspects were addressed based on 186 single nucleotide polymorphisms, extracted from sequences of 801 bp of Cytb gene, gathered from 416 individuals collected at 103 localities in Argentina and Chile. Employing multivariate statistical analyses (gPCA, Mantel and Partial Mantel Tests, Procrustes Analysis, and RDA), associations between genetic differences and geographic and climatic distances were evaluated. Presence data was employed to estimate the potential geographic distribution of this species during historical and contemporary climatic scenarios, and to address differences among the climatic niches of their main mitochondrial lineages. The significant influence of landscape features in structuring mitochondrial variability was evidenced at different spatial scales, as well as the role of past climatic dynamics in driving geographic range shifts, mostly associated to Quaternary glaciations. Overall, these results suggest that throughout geographic range gene flow is unevenly influenced by climatic dissimilarity and the geographic distancing, and that studied lineages do not exhibit distributional signals of climatic niche conservatism. Additionally, genetic differentiation occurred by more complex evolutionary processes than mere disruption of gene flow or drift.

Complex population history affects admixture analyses in nine-spined sticklebacks

Introgressive hybridization is an important process in evolution but challenging to identify, undermining the efforts to understand its role and significance. On the contrary, many analytical methods assume direct descent from a single common ancestor, and admixture among populations can violate their assumptions and lead to seriously biased results. A detailed analysis of 888 whole-genome sequences of nine-spined sticklebacks (Pungitius pungitius) revealed a complex pattern of population ancestry involving multiple waves of gene flow and introgression across northern Europe. The two recognized lineages were found to have drastically different histories, and their secondary contact zone was wider than anticipated, displaying a smooth gradient of foreign ancestry with some curious deviations from the expected pattern. Interestingly, the freshwater isolates provided peeks into the past and helped to understand the intermediate states of evolutionary processes. Our analyses and findings paint a detailed picture of the complex colonization history of northern Europe and provide backdrop against which introgression and its role in evolution can be investigated. However, they also expose the challenges in analyses of admixed populations and demonstrate how hidden admixture and colonization history misleads the estimation of admixture proportions and population split times.

Variable hybridization between two Lake Tanganyikan cichlid species in recent secondary contact

Closely related taxa frequently exist in sympatry before the evolution of robust reproductive barriers, which can lead to substantial gene flow. Post-divergence gene flow can promote several disparate trajectories of divergence ranging from the erosion of distinctiveness and eventual collapse of the taxa to the strengthening of reproductive isolation. Among many relevant factors, understanding the demographic history of divergence (e.g. divergence time, extent of historical gene flow) can be particularly informative when examining contemporary gene flow between closely related taxa because this history can influence gene flow's prevalence and consequences. Here, we used genotyping-by-sequencing data to investigate speciation and contemporary hybridization in two closely related and sympatrically distributed Lake Tanganyikan cichlid species in the genus Petrochromis. Demographic modeling supported a speciation scenario involving divergence in isolation followed by secondary contact with bidirectional gene flow. Further investigation of this recent gene flow found evidence of ongoing hybridization between the species that varied in extent between different co-occurring populations. Relationships between abundance and the degree of admixture across populations suggest that the availability of conspecific mates may influence patterns of hybridization. These results, together with the observation that sets of recently diverged cichlid taxa are generally geographically separated in the lake, suggest that ongoing speciation in Lake Tanganyikan cichlids relies on initial spatial isolation. Additionally, the spatially heterogeneous patterns of admixture between the Petrochromis species illustrates the complexities of hybridization when species are in recent secondary contact.

Maintenance of local adaptation despite gene flow in a coastal songbird

Adaptation to local environments is common in widespread species and the basis of ecological speciation. The song sparrow (Melospiza melodia) is a widespread, polytypic passerine that occurs in shrubland habitats throughout North America. We examined the population structure of two parapatric subspecies that inhabit different environments: the Atlantic song sparrow (M. m. atlantica), a coastal specialist, and the eastern song sparrow (M. m. melodia), a shrubland generalist. These populations lacked clear mitochondrial population structure, yet coastal birds formed a distinct nuclear genetic cluster. We found weak overall genomic differentiation between these subspecies, suggesting either recent divergence, extensive gene flow, or a combination thereof. There was a steep genetic cline at the transition to coastal habitats, consistent with isolation by environment, not isolation by distance. A phenotype under divergent selection, bill size, varied with the amount of coastal ancestry in transitional areas, but larger bill size was maintained in coastal habitats regardless of ancestry, further supporting a role for selection in the maintenance of these subspecies. Demographic modeling suggested a divergence history of limited gene flow followed by secondary contact, which has emerged as a common theme in adaptive divergence across taxa.

Does the Pachytene Checkpoint, a Feature of Meiosis, Filter Out Mistakes in Double-Strand DNA Break Repair and as a side-Effect Strongly Promote Adaptive Speciation?

This essay aims to explain two biological puzzles: why eukaryotic transcription units are composed of short segments of coding DNA interspersed with long stretches of non-coding (intron) DNA, and the near ubiquity of sexual reproduction. As is well known, alternative splicing of its coding sequences enables one transcription unit to produce multiple variants of each encoded protein. Additionally, padding transcription units with non-coding DNA (often many thousands of base pairs long) provides a readily evolvable way to set how soon in a cell cycle the various mRNAs will begin being expressed and the total amount of mRNA that each transcription unit can make during a cell cycle. This regulation complements control via the transcriptional promoter and facilitates the creation of complex eukaryotic cell types, tissues, and organisms. However, it also makes eukaryotes exceedingly vulnerable to double-strand DNA breaks, which end-joining break repair pathways can repair incorrectly. Transcription units cover such a large fraction of the genome that any mis-repair producing a reorganized chromosome has a high probability of destroying a gene. During meiosis, the synaptonemal complex aligns homologous chromosome pairs and the pachytene checkpoint detects, selectively arrests, and in many organisms actively destroys gamete-producing cells with chromosomes that cannot adequately synapse; this creates a filter favoring transmission to the next generation of chromosomes that retain the parental organization, while selectively culling those with interrupted transcription units. This same meiotic checkpoint, reacting to accidental chromosomal reorganizations inflicted by error-prone break repair, can, as a side effect, provide a mechanism for the formation of new species in sympatry. It has been a long-standing puzzle how something as seemingly maladaptive as hybrid sterility between such new species can arise. I suggest that this paradox is resolved by understanding the adaptive importance of the pachytene checkpoint, as outlined above.

Introgressive hybridisation between domestic pigs (Sus scrofa domesticus) and endemic Corsican wild boars (S. s. meridionalis): effects of human-mediated interventions

Owing to the intensified domestication process with artificial trait selection, introgressive hybridisation between domestic and wild species poses a management problem. Traditional free-range livestock husbandry, as practiced in Corsica and Sardinia, is known to facilitate hybridisation between wild boars and domestic pigs (Sus scrofa). Here, we assessed the genetic distinctness and genome-wide domestic pig ancestry levels of the Corsican wild boar subspecies S. s. meridionalis, with reference to its Sardinian conspecifics, employing a genome-wide single nucleotide polymorphism (SNP) assay and mitochondrial control region (mtCR) haplotypes. We also assessed the reliance of morphological criteria and the melanocortin-1 receptor (MC1R) coat colour gene to identify individuals with domestic introgression. While Corsican wild boars showed closest affinity to Sardinian and Italian wild boars compared to other European populations based on principal component analysis, the observation of previously undescribed mtCR haplotypes and high levels of nuclear divergence (Weir’s θ > 0.14) highlighted the genetic distinctness of Corsican S. s. meridionalis. Across three complementary analyses of mixed ancestry (i.e., STRUCTURE, PCADMIX, and ELAI), proportions of domestic pig ancestry were estimated at 9.5% in Corsican wild boars, which was significantly higher than in wild boars in Sardinia, where free-range pig keeping was banned in 2012. Comparison of morphologically pure- and hybrid-looking Corsican wild boars suggested a weak correlation between morphological criteria and genome-wide domestic pig ancestry. The study highlights the usefulness of molecular markers to assess the direct impacts of management practices on gene flow between domestic and wild species.

Divergence time estimation using ddRAD data and an isolation-with-migration model applied to water vole populations of Arvicola

Molecular dating methods of population splits are crucial in evolutionary biology, but they present important difficulties due to the complexity of the genealogical relationships of genes and past migrations between populations. Using the double digest restriction-site associated DNA (ddRAD) technique and an isolation-with-migration (IM) model, we studied the evolutionary history of water vole populations of the genus Arvicola, a group of complex evolution with fossorial and semi-aquatic ecotypes. To do this, we first estimated mutation rates of ddRAD loci using a phylogenetic approach. An IM model was then used to estimate split times and other relevant demographic parameters. A set of 300 ddRAD loci that included 85 calibrated loci resulted in good mixing and model convergence. The results showed that the two populations of A. scherman present in the Iberian Peninsula split 34 thousand years ago, during the last glaciation. In addition, the much greater divergence from its sister species, A. amphibius, may help to clarify the controversial taxonomy of the genus. We conclude that this approach, based on ddRAD data and an IM model, is highly useful for analyzing the origin of populations and species.

Faster-haplodiploid evolution under divergence-with-gene-flow: simulations and empirical data from pine-feeding hymenopterans

Although haplodiploidy is widespread in nature, the evolutionary consequences of this mode of reproduction are not well characterized. Here, we examine how genome-wide hemizygosity and a lack of recombination in haploid males affects genomic differentiation in populations that diverge via natural selection while experiencing gene flow. First, we simulated diploid and haplodiploid "genomes" (500-kb loci) evolving under an isolation-with-migration model with mutation, drift, selection, migration, and recombination; and examined differentiation at neutral sites both tightly and loosely linked to a divergently selected site. So long as there is divergent selection and migration, sex-limited hemizygosity and recombination cause elevated differentiation (i.e., produce a "faster-haplodiploid effect") in haplodiploid populations relative to otherwise equivalent diploid populations, for both recessive and codominant mutations. Second, we used genome-wide SNP data to model divergence history and describe patterns of genomic differentiation between sympatric populations of Neodiprion lecontei and N. pinetum, a pair of pine sawfly species (order: Hymenoptera; family: Diprionidae) that are specialized on different pine hosts. These analyses support a history of continuous gene exchange throughout divergence and reveal a pattern of heterogeneous genomic differentiation that is consistent with divergent selection on many unlinked loci. Third, using simulations of haplodiploid and diploid populations evolving according to the estimated divergence history of N. lecontei and N. pinetum, we found that divergent selection would lead to higher differentiation in haplodiploids. Based on these results, we hypothesize that haplodiploids undergo divergence-with-gene-flow and sympatric speciation more readily than diploids.

The first steps toward a global pandemic: Reconstructing the demographic history of parasite host switches in its native range

Host switching allows parasites to expand their niches. However, successful switching may require suites of adaptations and also may decrease performance on the old host. As a result, reductions in gene flow accompany many host switches, driving speciation. Because host switches tend to be rapid, it is difficult to study them in real-time, and their demographic parameters remain poorly understood. As a result, fundamental factors that control subsequent parasite evolution, such as the size of the switching population or the extent of immigration from the original host, remain largely unknown. To shed light on the host switching process, we explored how host switches occur in independent host shifts by two ectoparasitic honey bee mites (Varroa destructor and V. jacobsoni). Both switched to the western honey bee (Apis mellifera) after being brought into contact with their ancestral host (Apis cerana), ~70 and ~12 years ago, respectively. Varroa destructor subsequently caused worldwide collapses of honey bee populations. Using whole-genome sequencing on 63 mites collected in their native ranges from both the ancestral and novel hosts, we were able to reconstruct the known temporal dynamics of the switch. We further found multiple previously undiscovered mitochondrial lineages on the novel host, along with the genetic equivalent of tens of individuals that were involved in the initial host switch. Despite being greatly reduced, some gene flow remains between mites adapted to different hosts. Our findings suggest that while reproductive isolation may facilitate the fixation of traits beneficial for exploiting the new host, ongoing genetic exchange may allow genetic amelioration of inbreeding effects.

Population genetic structure of the carrot weevil ( Listronotus oregonensis ) in North America

Population genetic studies of insect pests enhance our ability to anticipate problems in agroecosystems, such as pest outbreaks, insecticide resistance, or expansions of the host range. This study focuses on geographic distance and host plant selection as potential determinants of genetic differentiation of the carrot weevil Listronotus oregonensis, a major pest of several apiaceous crops in North America. To undertake genetic studies on this species, we assembled the first complete genome sequence for L. oregonensis. Then, we used both haplotype discrimination with mitochondrial DNA (mtDNA) and a genotyping‐by‐sequencing (GBS) approach to characterize the genetic population structure. A total of 220 individuals were sampled from 17 localities in the provinces of Québec, Ontario, Nova Scotia (Canada), and the state of Ohio (USA). Our results showed significant genetic differences between distant populations across North America, indicating that geographic distance represents an important factor of differentiation for the carrot weevil. Furthermore, the GBS analysis revealed more different clusters than COI analysis between Québec and Nova Scotia populations, suggesting a recent differentiation in the latter province. In contrast, we found no clear evidence of population structure associated with the four cultivated apiaceous plants tested (carrot, parsley, celery, and celeriac) using populations from Québec. This first characterization of the genetic structure of the carrot weevil contributes to a better understanding of the gene flow of the species and helps to adapt local pest management measures to better control this agricultural pest.

Contrasting environmental drivers of genetic and phenotypic divergence in an Andean poison frog (Epipedobates anthonyi)

Phenotypic and genetic divergence are shaped by the homogenizing effects of gene flow and the differentiating processes of genetic drift and local adaptation. Herein, we examined the mechanisms that underlie phenotypic (size and color) and genetic divergence in 35 populations (535 individuals) of the poison frog Epipedobates anthonyi along four elevational gradients (0-1800 m asl) in the Ecuadorian Andes. We found phenotypic divergence in size and color despite relatively low genetic divergence at neutral microsatellite loci. Genetic and phenotypic divergence were both explained by landscape resistance between sites (isolation-by-resistance, IBR), likely due to a cold and dry mountain ridge between the northern and southern elevational transects that limits dispersal and separates two color morphs. Moreover, environmental differences among sites also explained genetic and phenotypic divergence, suggesting isolation-by-environment (IBE). When northern and southern transects were analyzed separately, genetic divergence was predicted either by distance (isolation-by-distance, IBD; northern) or environmental resistance between sites (IBR; southern). In contrast, phenotypic divergence was primarily explained by environmental differences among sites, supporting the IBE hypothesis. These results indicate that although distance and geographic barriers are important drivers of population divergence, environmental variation has a two-fold effect on population divergence. On the one hand, landscape resistance between sites reduces gene flow (IBR), while on the other hand, environmental differences among sites exert divergent selective pressures on phenotypic traits (IBE). Our work highlights the importance of studying both genetic and phenotypic divergence to better understand the processes of population divergence and speciation along ecological gradients.

OUP accepted manuscript

Introgressive hybridization plays a key role in adaptive evolution and species diversification in many groups of species. However, frequent hybridization and gene flow between species make estimation of the species phylogeny and key population parameters challenging. Here, we show that by accounting for phasing and using full-likelihood methods, introgression histories and population parameters can be estimated reliably from whole-genome sequence data. We employ the multispecies coalescent (MSC) model with and without gene flow to infer the species phylogeny and cross-species introgression events using genomic data from six members of the erato-sara clade of Heliconius butterflies. The methods naturally accommodate random fluctuations in genealogical history across the genome due to deep coalescence. To avoid heterozygote phasing errors in haploid sequences commonly produced by genome assembly methods, we process and compile unphased diploid sequence alignments and use analytical methods to average over uncertainties in heterozygote phase resolution. There is robust evidence for introgression across the genome, both among distantly related species deep in the phylogeny and between sister species in shallow parts of the tree. We obtain chromosome-specific estimates of key population parameters such as introgression directions, times and probabilities, as well as species divergence times and population sizes for modern and ancestral species. We confirm ancestral gene flow between the sara clade and an ancestral population of Heliconius telesiphe, a likely hybrid speciation origin for Heliconius hecalesia, and gene flow between the sister species Heliconius erato and Heliconius himera. Inferred introgression among ancestral species also explains the history of two chromosomal inversions deep in the phylogeny of the group. This study illustrates how a full-likelihood approach based on the MSC makes it possible to extract rich historical information of species divergence and gene flow from genomic data. [3s; bpp; gene flow; Heliconius; hybrid speciation; introgression; inversion; multispecies coalescent]

Genetic and morphological data demonstrate hybridization and backcrossing in a pair of salamanders at the far end of the speciation continuum

Deeply diverged yet hybridizing species provide a system to investigate the final stages of the speciation process. We study a hybridizing pair of salamander species-the morphologically and genetically drastically different newts Triturus cristatus and T. marmoratus-with a panel of 32 nuclear and mitochondrial genetic markers. Morphologically identified hybrids are mostly of the F1 generation and mothered by T. cristatus. The sex ratio of the F1 hybrid class is reciprocally skewed, with a preponderance of females in T. cristatus-mothered hybrids and males in T. marmoratus-mothered hybrids. This amounts to the Haldane effect operating in one direction of the cross. Deeper generation hybrids are occasionally produced, possibly including F1 hybrid × backcross hybrid offspring. Interspecific gene flow is low, yet skewed toward T. cristatus. This asymmetry may be caused by hybrid zone movement, with the superseding species being predisposed to introgression. The persisting gene flow between deeply differentiated species supports the notion that full genetic isolation may be selected against. Conversely, published morphological data suggest that introgressive hybridization is detrimental, with digital malformations occurring more frequently in the area of sympatry. Finally, to assist field identification, both within the area of natural range overlap and concerning anthropogenic introductions elsewhere, we document the phenotypical variation of two generations of hybrids compared with both parental species. We suggest that fluctuating range boundaries, ecological segregation, cytonuclear incompatibilities and hybrid breakdown through Bateson-Dobzhansky-Muller incompatibilities all contribute to species integrity, despite incomplete isolation during secondary contact.

Estimation of the strength of mate preference from mated pairs observed in the wild

A number of key processes in evolution are driven by individuals preferring mates with particular phenotypes. However, despite long‐standing interest, it is difficult to quantify the strength of mate preference from phenotypic observations in nature in a way that connects directly to key parameters in theoretical models. To bridge the gap between mathematical models and empirical data, we develop a novel maximum likelihood‐based method to estimate the strength and form of mate preference, where preference depends on traits expressed in both males and females. Using simulated data, we demonstrate that our method accurately infers model parameters, including the strength of mate preference and the optimal offset match between trait values in mated pairs when model assumptions are satisfied. Applying our method to two previous studies of assortative mating in marine gastropods and the European common frog, we support previous findings, but also give additional insight into the role of mate preference in each system. Our method can be generalized to a variety of plant and animal taxa that exhibit mating preferences to facilitate the testing of evolutionary hypotheses and link empirical data to theoretical models of assortative mating, sexual selection, and speciation.

Sea surface temperature, rather than land mass or geographic distance, may drive genetic differentiation in a species complex of highly dispersive seabirds

Seabirds, particularly Procellariiformes, are highly mobile organisms with a great capacity for long dispersal, though simultaneously showing high philopatry, two conflicting life-history traits that may lead to contrasted patterns of genetic population structure. Landmasses were suggested to explain differentiation patterns observed in seabirds, but philopatry, isolation by distance, segregation between breeding and nonbreeding zones, and oceanographic conditions (sea surface temperatures) may also contribute to differentiation patterns. To our knowledge, no study has simultaneously contrasted the multiple factors contributing to the diversification of seabird species, especially in the gray zone of speciation. We conducted a multilocus phylogeographic study on a widespread seabird species complex, the little shearwater complex, showing highly homogeneous morphology, which led to considerable taxonomic debate. We sequenced three mitochondrial and six nuclear markers on all extant populations from the Atlantic (lherminieri) and Indian Oceans (bailloni), that is, five nominal lineages from 13 populations, along with one population from the eastern Pacific Ocean (representing the dichrous lineage). We found sharp differentiation among populations separated by the African continent with both mitochondrial and nuclear markers, while only mitochondrial markers allowed characterizing the five nominal lineages. No differentiation could be detected within these five lineages, questioning the strong level of philopatry showed by these shearwaters. Finally, we propose that Atlantic populations likely originated from the Indian Ocean. Within the Atlantic, a stepping-stone process accounts for the current distribution. Based on our divergence time estimates, we suggest that the observed pattern of differentiation mostly resulted from historical and current variation in sea surface temperatures.

Genetic structure in Louisiana Iris species reveals patterns of recent and historical admixture

PREMISE

When divergent lineages come into secondary contact, reproductive isolation may be incomplete, thus providing an opportunity to investigate how speciation is manifested in the genome. The Louisiana Irises (Iris, series Hexagonae) comprise a group of three or more ecologically and reproductively divergent lineages that can produce hybrids where they come into contact. We estimated standing genetic variation to understand the current distribution of population structure in the Louisiana Irises.

METHODS

We used genotyping-by-sequencing techniques to sample the genomes of Louisiana Iris species across their ranges. We sampled 20 populations (n = 632 individuals) across 11,249 loci and used Entropy and PCA models to assess population genetic data.

RESULTS

We discovered evidence for interspecific gene flow in parts of the range. Our analysis revealed patterns of population structure at odds with widely accepted nominal taxonomy. We discovered undescribed hybrid populations, designated as belonging to the I. brevicaulis lineage. Iris nelsonii shared significant ancestry with only one of the purported parent species, I. fulva, evidence inconsistent with a hybrid origin.

CONCLUSIONS

This study provides several key findings important to the investigation of standing genetic variation in the Louisiana Iris species complex. Compared to the other nominal species, I. brevicaulis contains a large amount of genetic diversity. In addition, we discovered a previously unknown hybrid zone between I. brevicaulis and I. hexagona along the Texas coast. Finally, our results do not support the long-standing hypothesis that I. nelsonii has mixed ancestry from three parental taxa.

Phylogenomics and biogeography of Catalpa (Bignoniaceae) reveal incomplete lineage sorting and three dispersal events

Catalpa Scop. (Bignoniaceae) is a small genus (8 spp.) of trees that is disjunctly distributed among eastern Asia, eastern United States, and the West Indies. Catalpa bears beautiful inflorescences and have been cultivated as important ornamental trees for landscaping, gardening, and timber. However, the phylogenetic relationships and biogeographic history of the genus have remained unresolved. In this study, we used a large genomic dataset that includes data from the chloroplast (plastomes), and nuclear genomes (ITS and 5,759 single-copy nuclear genes) to reconstruct phylogenetic relationship within Catalpa, test interspecific gene flow events within the genus, and infer its biogeographic history. Our phylogenetic results indicate that Catalpa is monophyletic containing two main clades, section Catalpa and section Macrocatalpa. Section Catalpa is further divided into three subclades. While most relationships are congruent between the chloroplast and nuclear datasets, the position of C. ovata differs, likely due to incomplete lineage sorting. Interspecific gene flow events include C. bungei s.s. with vectors of inheritance from C. duclouxii and C. fargesii, supporting a combination of these three species and recognizing a broadly circumscribed C. bungei s.l. Our biogeographic study suggests three main dispersal events, two of which occurred during the Oligocene. The first dispersal event occurred from southwestern North America and Mexico into the Greater Antilles giving rise to the ancestor of the section of Macrocatalpa. The second dispersal event also occurred from southwestern North America and Mexico, but led to central and northern North America, subsequently reaching China through the Bering land bridge, and also reaching Europe through the North Atlantic land bridge. The third dispersal event took place in the Miocene from China to North America and gave rise to a clade composed of C. bignonioides and C. speciosa. This study uses a phylogenomic approach and biogeographical methods to infer the evolutionary history of Catalpa, highlighting issues associated with gene tree discordance, and suggesting that incomplete lineage sorting likely played an important role in the evolutionary history of Catalpa.

Delimitation despite discordance: Evaluating the species limits of a confounding species complex in the face of mitonuclear discordance

The delimitation of species is an essential pursuit of biology, and proper taxonomies are crucial for the assessment and conservation management of organismal diversity. However, delimiting species can be hindered by a number of factors including highly conserved morphologies (e.g., cryptic species), differences in criteria of species concepts, lineages being in the early stages of the speciation or divergence process, and discordance between gene topologies (e.g., mitonuclear discordance). Here we use a taxonomically confounded species complex of toads in Central America that exhibits extensive mitonuclear discordance to test delimitation hypotheses. Our investigation integrates mitochondrial sequences, nuclear SNPs, morphology, and macroecological data to determine which taxonomy best explains the divergence and evolutionary relationships among these toads. We found that a three species taxonomy following the distributions of the nuclear SNP haplotypes offers the best explanation of the species in this complex based off of the integrated data types. Due to the taxonomic instability of this group, we also discuss conservation concerns in the face of improper taxonomic delimitation. Our study provides an empirical and integrative hypothesis testing framework to assess species delimitation hypotheses in the face of cryptic morphology and mitonuclear discordance and highlights the importance that a stable taxonomy has over conservation-related actions.

Genetic variation for adaptive traits is associated with polymorphic inversions in Littorina saxatilis

Chromosomal inversions have long been recognized for their role in local adaptation. By suppressing recombination in heterozygous individuals, they can maintain coadapted gene complexes and protect them from homogenizing effects of gene flow. However, to fully understand their importance for local adaptation we need to know their influence on phenotypes under divergent selection. For this, the marine snail Littorina saxatilis provides an ideal study system. Divergent ecotypes adapted to wave action and crab predation occur in close proximity on intertidal shores with gene flow between them. Here, we used F2 individuals obtained from crosses between the ecotypes to test for associations between genomic regions and traits distinguishing the Crab‐/Wave‐adapted ecotypes including size, shape, shell thickness, and behavior. We show that most of these traits are influenced by two previously detected inversion regions that are divergent between ecotypes. We thus gain a better understanding of one important underlying mechanism responsible for the rapid and repeated formation of ecotypes: divergent selection acting on inversions. We also found that some inversions contributed to more than one trait suggesting that they may contain several loci involved in adaptation, consistent with the hypothesis that suppression of recombination within inversions facilitates differentiation in the presence of gene flow.

A Revised Phylogeny of the Mentha spicata Clade Reveals Cryptic Species

The genus Mentha is taxonomically and phylogenetically challenging due to complex genomes, polyploidization and an extensive historical nomenclature, potentially hiding cryptic taxa. A straightforward interpretation of phylogenetic relationships within the section Mentha is further hindered by dominant but outdated concepts on historically identified hybrid taxa. Mentha spicata is traditionally considered to be of hybrid origin, but the evidence for this is weak. Here, we aim to understand the phylogenetic relationships within the section Mentha using large sample sizes and to revisit the hybrid status and identity of M. spicata. We show that two of three traditional species in the subsection Spicatae are polyphyletic, as is the subsection as a whole, while the real number of cryptic species was underestimated. Compared to previous studies we present a fundamentally different phylogeny, with a basal split between M. spicata s.s. and M. longifolia s.s. Cluster analyses of morphological and genotypic data demonstrate that there is a dissociation between morphologically and genotypically defined groups of samples. We did not find any evidence that M. spicata is of hybrid origin, and we conclude its taxonomic status should be revised. The combination of genetic and phenotypic information is essential when evaluating hyperdiverse taxonomic groups.

Adaptive Divergence under Gene Flow along an Environmental Gradient in Two Coexisting Stickleback Species

There is a general and solid theoretical framework to explain how the interplay between natural selection and gene flow affects local adaptation. Yet, to what extent coexisting closely related species evolve collectively or show distinctive evolutionary responses remains a fundamental question. To address this, we studied the population genetic structure and morphological differentiation of sympatric three-spined and nine-spined stickleback. We conducted genotyping-by-sequencing and morphological trait characterisation using 24 individuals of each species from four lowland brackish water (LBW), four lowland freshwater (LFW) and three upland freshwater (UFW) sites in Belgium and the Netherlands. This combination of sites allowed us to contrast populations from isolated but environmentally similar locations (LFW vs. UFW), isolated but environmentally heterogeneous locations (LBW vs. UFW), and well-connected but environmentally heterogenous locations (LBW vs. LFW). Overall, both species showed comparable levels of genetic diversity and neutral genetic differentiation. However, for all three spatial scales, signatures of morphological and genomic adaptive divergence were substantially stronger among populations of the three-spined stickleback than among populations of the nine-spined stickleback. Furthermore, most outlier SNPs in the two species were associated with local freshwater sites. The few outlier SNPs that were associated with the split between brackish water and freshwater populations were located on one linkage group in three-spined stickleback and two linkage groups in nine-spined stickleback. We conclude that while both species show congruent evolutionary and genomic patterns of divergent selection, both species differ in the magnitude of their response to selection regardless of the geographical and environmental context.

Inference of gene flow in the process of speciation: Efficient maximum-likelihood implementation of a generalised isolation-with-migration model

The 'isolation with migration' (IM) model has been extensively used in the literature to detect gene flow during the process of speciation. In this model, an ancestral population split into two or more descendant populations which subsequently exchanged migrants at a constant rate until the present. Of course, the assumption of constant gene flow until the present is often over-simplistic in the context of speciation. In this paper, we consider a 'generalised IM' (GIM) model: a two-population IM model in which migration rates and population sizes are allowed to change at some point in the past. By developing a maximum-likelihood implementation of this model, we enable inference on both historical and contemporary rates of gene flow between two closely related populations or species. The GIM model encompasses both the standard two-population IM model and the 'isolation with initial migration' (IIM) model as special cases, as well as a model of secondary contact. We examine for simulated data how our method can be used, by means of likelihood ratio tests or AIC scores, to distinguish between the following scenarios of population divergence: (a) divergence in complete isolation; (b) divergence with a period of gene flow followed by isolation; (c) divergence with a period of isolation followed by secondary contact; (d) divergence with ongoing gene flow. Our method is based on the coalescent and is suitable for data sets consisting of the number of nucleotide differences between one pair of DNA sequences at each of a large number of independent loci. As our method relies on an explicit expression for the likelihood, it is computationally very fast.

Demographic consequences of dispersal-related trait shift in two recently diverged taxa of montane grasshoppers

Although the pervasiveness of intraspecific wing-size polymorphism and transitions to flightlessness have long captivated biologists, the demographic outcomes of shifts in dispersal ability are not yet well understood and have been seldom studied at early stages of diversification. Here, we use genomic data to infer the consequences of dispersal-related trait variation in the taxonomically controversial short-winged (Chorthippus corsicus corsicus) and long-winged (Chorthippus corsicus pascuorum) Corsican grasshoppers. Our analyses revealed lack of contemporary hybridization between sympatric long- and short-winged forms and phylogenomic reconstructions supported their taxonomic distinctiveness, rejecting the hypothesis of intraspecific wing polymorphism. Statistical evaluation of alternative models of speciation strongly supported a scenario of Pleistocene divergence (<1.5 Ma) with ancestral gene flow. According to neutral expectations from differences in dispersal capacity, historical effective migration rates from the long- to the short-winged taxon were threefold higher than in the opposite direction. Although populations of the two taxa present a marked genetic structure and have experienced parallel demographic histories, our coalescent-based analyses suggest that reduced dispersal has fueled diversification in the short-winged C. c. corsicus. Collectively, our study illustrates how dispersal reduction can speed up geographical diversification and increase the opportunity for allopatric speciation in topographically complex landscapes.

Extensive Genome-Wide Phylogenetic Discordance Is Due to Incomplete Lineage Sorting and Not Ongoing Introgression in a Rapidly Radiated Bryophyte Genus

The relative importance of introgression for diversification has long been a highly disputed topic in speciation research and remains an open question despite the great attention it has received over the past decade. Gene flow leaves traces in the genome similar to those created by incomplete lineage sorting (ILS), and identification and quantification of gene flow in the presence of ILS is challenging and requires knowledge about the true phylogenetic relationship among the species. We use whole nuclear, plastid, and organellar genomes from 12 species in the rapidly radiated, ecologically diverse, actively hybridizing genus of peatmoss (Sphagnum) to reconstruct the species phylogeny and quantify introgression using a suite of phylogenomic methods. We found extensive phylogenetic discordance among nuclear and organellar phylogenies, as well as across the nuclear genome and the nodes in the species tree, best explained by extensive ILS following the rapid radiation of the genus rather than by postspeciation introgression. Our analyses support the idea of ancient introgression among the ancestral lineages followed by ILS, whereas recent gene flow among the species is highly restricted despite widespread interspecific hybridization known in the group. Our results contribute to phylogenomic understanding of how speciation proceeds in rapidly radiated, actively hybridizing species groups, and demonstrate that employing a combination of diverse phylogenomic methods can facilitate untangling complex phylogenetic patterns created by ILS and introgression.

Speciation in Daphnia

The microcrustacean Daphnia is arguably one of the most studied zooplankton species, having a well understood ecology, life history, and a relatively well studied evolutionary history. Despite this wealth of knowledge, species boundaries within closely related species in this genus often remain elusive and the major evolutionary forces driving the diversity of daphniids remain controversial. This genus contains more than 80 species with multiple cryptic species complexes, with many closely related species able to hybridize. Here, we review speciation research in Daphnia within the framework of current speciation theory. We evaluate the role of geography, ecology, and biology in restricting gene flow and promoting diversification. Of the 253 speciation studies on Daphnia, the majority of studies examine geographic barriers (55%). While evidence shows that geographic barriers play a role in species divergence, ecological barriers are also probably prominent in Daphnia speciation. We assess the contribution of ecological and nonecological reproductive isolating barriers between closely related species of Daphnia and found that none of the reproductive isolating barriers are restricting gene flow completely. Research on reproductive isolating barriers has disproportionally focused on two species complexes, the Daphnia pulex and Daphnia longispina species complexes. Finally, we identify areas of research that remain relatively unexplored and discuss future research directions that build our understanding of speciation in daphniids.

Variable levels of introgression between the endangered Podarcis carbonelli and highly divergent congeneric species

Recent empirical studies have demonstrated that speciation with gene flow is more common than previously thought. From a conservation perspective, the potential negative effects of hybridization raise concerns on the genetic integrity of endangered species. However, introgressive hybridization has also been growingly recognized as a source of diversity and new advantageous alleles. Carbonell's wall lizard (Podarcis carbonelli) is an endangered species whose distribution overlaps with four other congeneric species. Our goal here was to determine whether P. carbonelli is completely reproductively isolated from its congeners and to evaluate the relevance of hybridization and interspecific gene flow for developing a conservation plan. We used restriction site associated DNA (RAD) sequencing to discover SNPs in samples from four contact zones between P. carbonelli and four other species. Principal component analysis, multilocus genotype assignment and interspecific heterozygosity suggest incomplete reproductive isolation and ongoing gene flow between species. However, hybridization dynamics vary across all pairs, suggesting complex interactions between multiple intrinsic and extrinsic barriers. Despite seemingly ubiquitous interspecific gene flow, we found evidence of strong reproductive isolation across most contact zones. Instead, indirect effects of hybridization like waste of reproductive effort in small isolated populations may be more problematic. Our results highlight the need to further evaluate the consequences of introgression for P. carbonelli, both on a geographic and genomic level and included in a comprehensive and urgently needed conservation plan. Besides, those findings will add important insights on the potential effects of hybridization and introgression for endangered species.

Flower evolution in the presence of heterospecific gene flow and its contribution to lineage divergence

The success of species depends on their ability to exploit ecological resources in order to optimize their reproduction. However, species are not usually found within single-species ecosystems but in complex communities. Because of their genetic relatedness, closely related lineages tend to cluster within the same ecosystem, rely on the same resources, and be phenotypically similar. In sympatry, they will therefore compete for the same resources and, in the case of flowering plants, exchange their genes through heterospecific pollen transfer. These interactions, nevertheless, pose significant challenges to species co-existence because they can lead to resource limitation and reproductive interference. In such cases, divergent selective pressures on floral traits will favour genotypes that isolate or desynchronize the reproduction of sympatric lineages. The resulting displacement of reproductive characters will, in turn, lead to pre-mating isolation and promote intraspecific divergence, thus initiating or reinforcing the speciation process. In this review, we discuss the current theoretical and empirical knowledge on the influence of heterospecific pollen transfer on flower evolution, highlighting its potential to uncover the ecological and genomic constraints shaping the speciation process.

Reconciling seascape genetics and fisheries science in three codistributed flatfishes

Uncertainty hampers innovative mixed-fisheries management by the scales at which connectivity dynamics are relevant to management objectives. The spatial scale of sustainable stock management is species-specific and depends on ecology, life history and population connectivity. One valuable approach to understand these spatial scales is to determine to what extent population genetic structure correlates with the oceanographic environment. Here, we compare the level of genetic connectivity in three codistributed and commercially exploited demersal flatfish species living in the North East Atlantic Ocean. Population genetic structure was analysed based on 14, 14 and 10 neutral DNA microsatellite markers for turbot, brill and sole, respectively. We then used redundancy analysis (RDA) to attribute the genetic variation to spatial (geographical location), temporal (sampling year) and oceanographic (water column characteristics) components. The genetic structure of turbot was composed of three clusters and correlated with variation in the depth of the pycnocline, in addition to spatial factors. The genetic structure of brill was homogenous, but correlated with average annual stratification and spatial factors. In sole, the genetic structure was composed of three clusters, but was only linked to a temporal factor. We explored whether the management of data poor commercial fisheries, such as in brill and turbot, might benefit from population-specific information. We conclude that the management of fish stocks has to consider species-specific genetic structures and may benefit from the documentation of the genetic seascape and life-history traits.

Recent diversification in the high Andes: unveiling the evolutionary history of the Ecuadorian hillstar,Oreotrochilus chimborazo(Apodiformes: Trochilidae)

Studying the genetic signatures of evolutionary diversification in young lineages is among the most promising approaches for unveiling the processes behind speciation. Here, we focus on Oreotrochilus chimborazo, a high Andean species of hummingbird that might have experienced rapid diversification in the recent past. To understand the evolution of this species, we generated a dataset of ten microsatellite markers and complementary data on morphometrics, plumage variation and ecological niches. We applied a series of population and coalescent-based analyses to understand the population structure and differentiation within the species, in addition to the signatures of current and historical gene flow, the location of potential contact zones and the relationships among lineages. We found that O. chimborazo comprises three genetic groups: one corresponding to subspecies O. c. chimborazo, from Chimborazo volcano and surroundings, and two corresponding to the northern and southern ranges of subspecies O. c. jamesonii, found from the extreme south of Colombia to southern Ecuador. We inferred modest levels of both contemporary and historical gene flow and proposed the location of a contact zone between lineages. Also, our coalescent-based analyses supported a rapid split among these three lineages during the mid-to-late Holocene. We discuss our results in the light of past and present potential distributions of the species, in addition to evolutionary trends seen in other Andean hummingbirds.

Natural selection maintains species despite frequent hybridization in the desert shrub Encelia

Natural selection is an important driver of genetic and phenotypic differentiation between species. For species in which potential gene flow is high but realized gene flow is low, adaptation via natural selection may be a particularly important force maintaining species. For a recent radiation of New World desert shrubs (Encelia: Asteraceae), we use fine-scale geographic sampling and population genomics to determine patterns of gene flow across two hybrid zones formed between two independent pairs of species with parapatric distributions. After finding evidence for extremely strong selection at both hybrid zones, we use a combination of field experiments, high-resolution imaging, and physiological measurements to determine the ecological basis for selection at one of the hybrid zones. Our results identify multiple ecological mechanisms of selection (drought, salinity, herbivory, and burial) that together are sufficient to maintain species boundaries despite high rates of hybridization. Given that multiple pairs of Encelia species hybridize at ecologically divergent parapatric boundaries, such mechanisms may maintain species boundaries throughout Encelia.

Molecular systematics and phylogeography of a widespread Neotropical avian lineage: evidence for cryptic speciation with protracted gene flow throughout the Late Quaternary

Here we use an integrative approach, including coalescent-based methods, isolation–migration and species distribution models, to infer population structure, divergence times and diversification in the two species of the genus Cymbilaimus (Aves, Thamnophilidae). Our results support a recent and rapid diversification with both incomplete lineage sorting and gene flow shaping the evolutionary history of Cymbilaimus. The spatio-temporal pattern of cladogenesis suggests that Cymbilaimus originated in the north/western portion of cis-Andean South America and then diversified into the Brazilian Shield and Central America after consolidation of the modern Amazonian drainage and the Andean range. This evolutionary scenario is explained by cycles of range expansion and dispersal, followed by isolation, and recurrent gene flow, during the last 1.2 Myr. Our results agree with those recently reported for other closely related suboscine lineages, whereby the window of introgression between closely related taxa remains open for up to a few million years after their original split. In Cymbilaimus, introgression was recurrent between C. lineatus and C. sanctaemariae, even after they acquired vocal and ecological differentiation, supporting the claim that at least in Neotropical suboscines, full reproductive compatibility may take millions of years to evolve and cannot be interpreted as synonymous with a lack of speciation.

The Perfect Storm: Gene Tree Estimation Error, Incomplete Lineage Sorting, and Ancient Gene Flow Explain the Most Recalcitrant Ancient Angiosperm Clade, Malpighiales

The genomic revolution offers renewed hope of resolving rapid radiations in the Tree of Life. The development of the multispecies coalescent model and improved gene tree estimation methods can better accommodate gene tree heterogeneity caused by incomplete lineage sorting (ILS) and gene tree estimation error stemming from the short internal branches. However, the relative influence of these factors in species tree inference is not well understood. Using anchored hybrid enrichment, we generated a data set including 423 single-copy loci from 64 taxa representing 39 families to infer the species tree of the flowering plant order Malpighiales. This order includes 9 of the top 10 most unstable nodes in angiosperms, which have been hypothesized to arise from the rapid radiation during the Cretaceous. Here, we show that coalescent-based methods do not resolve the backbone of Malpighiales and concatenation methods yield inconsistent estimations, providing evidence that gene tree heterogeneity is high in this clade. Despite high levels of ILS and gene tree estimation error, our simulations demonstrate that these two factors alone are insufficient to explain the lack of resolution in this order. To explore this further, we examined triplet frequencies among empirical gene trees and discovered some of them deviated significantly from those attributed to ILS and estimation error, suggesting gene flow as an additional and previously unappreciated phenomenon promoting gene tree variation in Malpighiales. Finally, we applied a novel method to quantify the relative contribution of these three primary sources of gene tree heterogeneity and demonstrated that ILS, gene tree estimation error, and gene flow contributed to 10.0$\%$, 34.8$\%$, and 21.4$\%$ of the variation, respectively. Together, our results suggest that a perfect storm of factors likely influence this lack of resolution, and further indicate that recalcitrant phylogenetic relationships like the backbone of Malpighiales may be better represented as phylogenetic networks. Thus, reducing such groups solely to existing models that adhere strictly to bifurcating trees greatly oversimplifies reality, and obscures our ability to more clearly discern the process of evolution. [Coalescent; concatenation; flanking region; hybrid enrichment, introgression; phylogenomics; rapid radiation, triplet frequency.].

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.].