A model-based method for identifying species hybrids using multilocus genetic data

Characterization of the cryptic interspecific hybrid Lemna×mediterranea by an integrated approach provides new insights into duckweed diversity

Lemnaceae taxonomy is challenged by the particular morphology of these tiny free-floating angiosperms. Although molecular taxonomy has helped clarify the phylogenetic history of this family, some inconsistency with morphological data leads to frequent misclassifications in the genus Lemna. Recently, the finding that Lemna japonica is an interspecific hybrid between Lemna minor and Lemna turionifera provided a clear explanation for one such taxonomic question. Here we demonstrated that L. minor is also capable of hybridizing with Lemna gibba, generating a cryptic but widespread taxon in the Mediterranean area. The nothotaxon Lemna ×mediterranea is described and compared with clones of the putative parental species L. minor and L. gibba. Genetic analysis by nuclear and plastid markers, as well as genome size measurement, revealed that two different cytotypes, diploid and triploid, originated by at least two independent hybridization events. Despite high overall similarity, morphometrical, physiological, and biochemical analyses showed an intermediate position of L. ×mediterranea between its parental species in most qualitative and quantitative characters, and also separation of the two hybrid cytotypes by some criteria. These data provide evidence that hybridization and polyploidization, driving forces of terrestrial plant evolution, contribute to duckweed genetic diversity and may have shaped the phylogenetic history of these mainly asexual, aquatic plants.

Ecological divergence despite common mating sites: Genotypes and symbiotypes shed light on cryptic diversity in the black bean aphid species complex

Different host plants represent ecologically dissimilar environments for phytophagous insects. The resulting divergent selection can promote the evolution of specialized host races, provided that gene flow is reduced between populations feeding on different plants. In black bean aphids belonging to the Aphis fabae complex, several morphologically cryptic taxa have been described based on their distinct host plant preferences. However, host choice and mate choice are largely decoupled in these insects: they are host-alternating and migrate between specific summer host plants and shared winter hosts, with mating occurring on the shared hosts. This provides a yearly opportunity for gene flow among aphids using different summer hosts, and raises the question if and to what extent the ecologically defined taxa are reproductively isolated. Here, we analyzed a geographically and temporally structured dataset of microsatellite genotypes from A. fabae that were mostly collected from their main winter host Euonymus europaeus, and additionally from another winter host and fourteen summer hosts. The data reveals multiple, strongly differentiated genetic clusters, which differ in their association with different summer and winter hosts. The clusters also differ in the frequency of infection with two heritable, facultative endosymbionts, separately hinting at reproductive isolation and divergent ecological selection. Furthermore, we found evidence for occasional hybridization among genetic clusters, with putative hybrids collected more frequently in spring than in autumn. This suggests that similar to host races in other phytophagous insects, both prezygotic and postzygotic barriers including selection against hybrids maintain genetic differentiation among A. fabae taxa, despite a common mating habitat.

Widespread White Matter Abnormalities in Concussed Athletes Detected by 7T Diffusion Magnetic Resonance Imaging

Sports-related concussions may cause white matter injuries and persistent post-concussive symptoms (PPCS). We hypothesized that athletes with PPCS would have neurocognitive impairments and white matter abnormalities that could be revealed by advanced neuroimaging using ultra-high field strength diffusion tensor (DTI) and diffusion kurtosis (DKI) imaging metrics and cerebrospinal fluid (CSF) biomarkers. A cohort of athletes with PPCS severity limiting the ability to work/study and participate in sport school and/or social activities for ≥6 months completed 7T magnetic resonance imaging (MRI) (morphological T1-weighed volumetry, DTI and DKI), extensive neuropsychological testing, symptom rating, and CSF biomarker sampling. Twenty-two athletes with PPCS and 22 controls were included. Concussed athletes performed below norms and significantly lower than controls on all but one of the psychometric neuropsychology tests. Supratentorial white and gray matter, as well as hippocampal volumes did not differ between concussed athletes and controls. However, of the 72 examined white matter tracts, 16% of DTI and 35% of DKI metrics (in total 28%) were significantly different between concussed athletes and controls. DKI fractional anisotropy and axial kurtosis were increased, and DKI radial diffusivity and radial kurtosis decreased in concussed athletes when compared with controls. CSF neurofilament light (NfL; an axonal injury marker), although not glial fibrillary acidic protein, correlated with several diffusion metrics. In this first 7T DTI and DKI study investigating PPCS, widespread microstructural alterations were observed in the white matter, correlating with CSF markers of axonal injury. More white matter changes were observed using DKI than using DTI. These white matter alterations may indicate persistent pathophysiological processes following concussion in sport.

Hybrids of two destructive subterranean termites established in the field, revealing a potential for gene flow between species

Hybridization between invasive pest species may lead to significant genetic and economic impacts that require close monitoring. The two most invasive and destructive termite species worldwide, Coptotermes formosanus Shiraki and Coptotermes gestroi (Wasmann), have the potential for hybridization in the field. A three-year field survey conducted during the dispersal flight season of Coptotermes in Taiwan identified alates with atypical morphology, which were confirmed as hybrids of the two Coptotermes species using microsatellite and mitochondrial analyses. Out of 27,601 alates collected over three years, 4.4% were confirmed as hybrid alates, and some advanced hybrids (>F1 generations) were identified. The hybrid alates had a dispersal flight season that overlapped with the two parental species 13 out of 15 times. Most of the hybrid alates were females, implying that mating opportunities beyond F1 may primarily be possible through female hybrids. However, the incipient colony growth results from all potential mating combinations suggest that only backcross colonies with hybrid males could sometimes lead to brood development. The observed asymmetrical viability and fertility of hybrid alates may critically reduce the probability of advanced-hybrid colonies being established in the field.

Spatial genetic structure in European marbled newts revealed with target enrichment by sequence capture

European marbled newts come in two species that have abutting ranges. The northern species, Triturus marmoratus, is found in France and the northern part of the Iberian Peninsula, whereas the southern species, T. pygmaeus, is found in the southwestern corner of the Iberian Peninsula. We study the intraspecific genetic differentiation of the group because morphological data show geographical variation and because the Iberian Peninsula is a recognized center of speciation and intraspecific genetic diversity for all kinds of organisms, amphibians included. We use target enrichment by sequence capture to generate c. 7 k nuclear DNA markers. We observe limited genetic exchange between the species, which confirms their distinctiveness. Both species show substantial genetic structuring that is only in part mirrored by morphological variation. Genetically differentiated groups are found in the south (T. marmoratus) and west (T. pygmaeus) of the species ranges. Our observations highlight the position of the Iberian Peninsula as a hotspot for genetic differentiation.

Yosemite toad (Anaxyrus canorus) transcriptome reveals interplay between speciation genes and adaptive introgression

Genomes are heterogeneous during the early stages of speciation, with small 'islands' of DNA appearing to reflect strong adaptive differences, surrounded by vast seas of relative homogeneity. As species diverge, secondary contact zones between them can act as an interface and selectively filter through advantageous alleles of hybrid origin. Such introgression is another important adaptive process, one that allows beneficial mosaics of recombinant DNA ('rivers') to flow from one species into another. Although genomic islands of divergence appear to be associated with reproductive isolation, and genomic rivers form by adaptive introgression, it is unknown whether islands and rivers tend to be the same or different loci. We examined three replicate secondary contact zones for the Yosemite toad (Anaxyrus canorus) using two genomic data sets and a morphometric data set to answer the questions: (1) How predictably different are islands and rivers, both in terms of genomic location and gene function? (2) Are the adaptive genetic trait loci underlying tadpole growth and development reliably islands, rivers or neither? We found that island and river loci have significant overlap within a contact zone, suggesting that some loci are first islands, and later are predictably converted into rivers. However, gene ontology enrichment analysis showed strong overlap in gene function unique to all island loci, suggesting predictability in overall gene pathways for islands. Genome-wide association study outliers for tadpole development included LPIN3, a lipid metabolism gene potentially involved in climate change adaptation, that is island-like for all three contact zones, but also appears to be introgressing (as a river) across one zone. Taken together, our results suggest that adaptive divergence and introgression may be more complementary forces than currently appreciated.

Rolling with the punches: Organism-environment interactions shape spatial pattern of adaptive differentiation in the widespread mantis shrimp Oratosquilla oratoria

Investigating spatial pattern of adaptive variation and its underlying processes can inform the adaptive potential distributed within species ranges, which is increasingly important in the context of a changing climate. A correct interpretation of adaptive variation pattern requires that population history and the ensuing population genetic structure are taken into account. Here we carried out such a study by integrating population genomic analyses, demographic model testing and species distribution modeling to investigate patterns and causes of adaptive differentiation in a widespread mantis shrimp, Oratosquilla oratoria, along a replicated, broad-scale temperature gradient in the northwestern Pacific (NWP). Our results supported a strong hierarchical ecogeographic structure dominated by habitat-linked divergence among O. oratoria populations accompanied with introgressive hybridization. A combined FST outlier and environmental correlation analyses revealed remarkable temperature-associated clines in allele frequency across paired North-South populations on Chinese and Japanese coasts, and identified a suite of loci associated with temperature adaptation. Further demographic model testing revealed the observed clinal variation derived partly from Pleistocene divergence followed by recent secondary contact. More importantly, the likelihood of hybridization is predicted to increase as climate change progresses, which would break barriers to gene flow and enable the spread of adaptive genetic variation. These results support that not only is temperature-driven adaptive differentiation occurs in O. oratoria but that such pattern is likely attributed to ancient adaptive variation, sustained by contemporary ocean conditions and a semi-permeable barrier to gene flow maintained by selection. They moreover provide genomic insights into the distribution of adaptive potential across O. oratoria' s species range. This work can serve as a case study to characterize adaptive diversity of marine species in the NWP by integrating environmental and genetic data at temporal and spatial scales in a population genomic framework, which would improve management and conservation actions under climate change.

Sympatry in a nightingale contact zone has no effect on host-specific blood parasite prevalence and lineage diversity

Parasites are a key driving force behind many ecological and evolutionary processes. Prevalence and diversity of parasites, as well as their effects on hosts, are not uniform across host species. As such, the potential parasite spillover between species can significantly influence outcomes of interspecific interactions. We screened two species of Luscinia nightingales for haemosporidian blood parasites (Plasmodium, Leucocytozoon and Haemoproteus) along an approximately 3000 km transect in Europe, incorporating areas of host distant allopatry, close allopatry and sympatry. We found significant differences in infection rates between the two host species, with common nightingales having much lower parasite prevalence than thrush nightingales (36.7% versus 83.8%). This disparity was mostly driven by Haemoproteus prevalence, which was significantly higher in thrush nightingales while common nightingales had a small, but significantly higher, Plasmodium prevalence. Furthermore, we found no effect of proximity to the contact zone on infection rate in either host species. Despite having lower infection prevalence, common nightingales were infected with a significantly higher diversity of parasite lineages than thrush nightingales, and lineage assemblages differed considerably between the two species, even in sympatry. This pattern was mostly driven by the large diversity of comparatively rare lineages, while the most abundant lineages were shared between the two host species. This suggests that, despite the close evolutionary relationships between the two nightingales, there are significant differences in parasite prevalence and diversity, regardless of the distance from the contact zone. This suggests that spillover of haemosporidian blood parasites is unlikely to contribute towards interspecific interactions in this system.

Glacial vicariance and secondary contact shape demographic histories in a freshwater mussel species complex

Characterizing the mechanisms influencing the distribution of genetic variation in aquatic species can be difficult due to the dynamic nature of hydrological landscapes. In North America's Central Highlands, a complex history of glacial dynamics, long-term isolation, and secondary contact have shaped genetic variation in aquatic species. Although the effects of glacial history have been demonstrated in many taxa, responses are often lineage- or species-specific and driven by organismal ecology. In this study, we reconstruct the evolutionary history of a freshwater mussel species complex using a suite of mitochondrial and nuclear loci to resolve taxonomic and demographic uncertainties. Our findings do not support Pleurobema rubrum as a valid species, which is proposed for listing as threatened under the U.S. Endangered Species Act. We synonymize P. rubrum under Pleurobema sintoxia-a common and widespread species found throughout the Mississippi River Basin. Further investigation of patterns of genetic variation in P. sintoxia identified a complex demographic history, including ancestral vicariance and secondary contact, within the Eastern Highlands. We hypothesize these patterns were shaped by ancestral vicariance driven by the formation of Lake Green and subsequent secondary contact after the last glacial maximum. Our inference aligns with demographic histories observed in other aquatic taxa in the region and mirrors patterns of genetic variation of a freshwater fish species (Erimystax dissimilis) confirmed to serve as a parasitic larval host for P. sintoxia. Our findings directly link species ecology to observed patterns of genetic variation and may have significant implications for future conservation and recovery actions of freshwater mussels.

simRestore: A decision-making tool for adaptive management of the native genetic status of wild populations

Anthropogenic hybridization, or higher and non-natural rates of gene flow directly and indirectly induced by human activities, is considered a significant threat to biodiversity. The primary concern for conservation is the potential for genomic extinction and loss of adaptiveness for native species due to the extensive introgression of non-native genes. To alleviate or reverse trends for such scenarios requires the direct integration of genomic data within a model framework for effective management. Towards this end, we developed the simRestore R program as a decision-making tool that integrates ecological and genomic information to simulate ancestry outcomes from optimized conservation strategies. In short, the program optimizes supplementation and removal strategies across generations until a set native genetic threshold is reached within the studied population. Importantly, in addition to helping with initial decision-making, simulations can be updated with the outcomes of ongoing efforts, allowing for the adaptive management of populations. After demonstrating functionality, we apply and optimize among actionable management strategies for the endangered Hawaiian duck for which the current primary threat is genetic extinction through ongoing anthropogenic hybridization with feral mallards. Simulations demonstrate that supplemental and removal efforts can be strategically tailored to move the genetic ancestry of Hawaii's hybrid populations towards Hawaiian duck without the need to completely start over. Further, we discuss ecological parameter sensitivity, including which factors are most important to ensure genetic outcomes (i.e. number of offspring). Finally, to facilitate use, the program is also available online as a Shiny Web application.

Patterns of admixture and introgression in a mosaic Mytilus galloprovincialis and Mytilus edulis hybrid zone in SW England

The study of hybrid zones offers important insights into speciation. Earlier studies on hybrid populations of the marine mussel species Mytilus edulis and Mytilus galloprovincialis in SW England provided evidence of admixture but were constrained by the limited number of molecular markers available. We use 57 ancestry-informative SNPs, most of which have been mapped genetically, to provide evidence of distinctive differences between admixed populations in SW England and asymmetrical introgression from M. edulis to M. galloprovincialis. We combine the genetic study with analysis of phenotypic traits of potential ecological and adaptive significance. We demonstrate that hybrid individuals have brown mantle edges unlike the white or purple in the parental species, suggesting allelic or non-allelic genomic interactions. We report differences in gonad development stage between the species consistent with a prezygotic barrier between the species. By incorporating results from publications dating back to 1980, we confirm the long-term stability of the hybrid zone despite higher viability of M. galloprovincialis. This stability coincides with a dramatic change in temperature of UK coastal waters and suggests that these hybrid populations might be resisting the effects of global warming. However, a single SNP locus associated with the Notch transmembrane signalling protein shows a markedly different pattern of variation to the others and might be associated with adaptation of M. galloprovincialis to colder northern temperatures.

Discovery of ex situ individuals of Andrias sligoi, an extremely endangered species and one of the largest amphibians worldwide

The South China giant salamander, Andrias sligoi, is one of the largest extant amphibian species worldwide. It was recently distinguished from another Chinese species, the Chinese giant salamander, Andrias davidianus, which is considered Critically Endangered according to the International Union for Conservation of Nature (IUCN) Red List. It appears too late to save this extremely rare and large amphibian in situ. Another extant species of the same genus, Andrias japonicus, inhabits Japan. However, the introduction of Chinese giant salamanders into some areas of Japan has resulted in hybridization between the Japanese and Chinese species. During our genetic screening of giant salamanders in Japan, we unexpectedly discovered four individuals of the South China giant salamander: two were adult males in captivity, and one had recently died. The last individual was a preserved specimen. In this study, we report these extremely rare individuals of A. sligoi in Japan and discuss the taxonomic and conservational implications of these introduced individuals.

SNP genotyping revealed a hybrid zone between Mytilus chilensis and M. platensis in southern South America (the Strait of Magellan, Isla Grande de Tierra del Fuego and the Falkland Islands)

Hybrid zones among mussel species have been extensively studied in the northern hemisphere. In South America, it has only recently become possible to study the natural hybrid zones, due to the clarification of the taxonomy of native mussels of the Mytilus genus. Analysing 54 SNP markers, we show the genetic species composition and admixture in the hybrid zone between M. chilensis and M. platensis in the southern end of South America. Bayesian, non-Bayesian clustering and re-assignment algorithms showed that the natural hybrid zone between M. chilensis and M. platensis in the Strait of Magellan, Isla Grande de Tierra del Fuego and the Falkland Islands shows clinal architecture. The hybrid zone can be divided into three different areas: the first one is on the Atlantic coast where only pure M. platensis and hybrid were found. In the second one, inside the Strait of Magellan, pure individuals of both species and mussels with variable degrees of hybridisation coexist. In the last area at the Strait in front of Punta Arenas City, fjords on the Isla Grande de Tierra del Fuego, and at the Beagle Channel, only M. chilensis and a low number of hybrids were found. According to the proportion of hybrids, bays with protected conditions away from strong currents would give better conditions for hybridisation. We do not find evidence of any other mussel species such as M. edulis, M. galloprovincialis, M. planulatus or M. trossulus in the zone.

Molecular and morphological evidence of hybridization between two dimorphic sympatric species of Fuchsia (Onagraceae)

Hybridization is commonly reported in angiosperms, generally based on morphology, and in few cases confirmed by molecular markers. Fuchsia has a long tradition of ornamental cultivars with different hybrids produced by artificial crosses, so natural hybridization between sympatric Fuchsia species could be common. Natural hybridization between F. microphylla and F. thymifolia was tested using six newly developed microsatellites for F. microphylla in addition to other molecular markers with codominant and maternal inheritance. Geometric morphometrics of leaves and floral structures were also used to identify putative hybrids. Hybrids showed a different degree of genetic admixture between both parental species. Chloroplast DNA (cpDNA) sequences indicated that hybridization occurs in both directions, in fact, some of the hybrids showed new haplotypes for cpDNA and ITS (internal transcriber spacer of nuclear ribosomal RNA genes) sequences. The morphology of hybrid individuals varied between the two parental species, but they could be better identified by their leaves and floral tubes. Our study is the first to confirm the hybridization in natural populations of Fuchsia species and suggests that hybridization has probably occurred repeatedly throughout the entire distribution of the species. Phylogeographic analysis of both species will be essential to understanding the impact of hybridization throughout their complete distribution.

Assessment of genetic diversity, population structure and wolf-dog hybridisation in the Eastern Romanian Carpathian wolf population

The Carpathian Mountains have been constantly inhabited by grey wolves and present one of the largest distribution areas in Europe, comprising between 2300 and 2700 individuals in Romania. To date, however, relatively little is known about the Romanian wolf population. We aimed to provide a first assessment of genetic diversity, population structure and wolf-dog hybridisation based on 444 mostly non-invasively collected samples in the Eastern Romanian Carpathians. Pack reconstruction and analysis of population genetic parameters were performed with mitochondrial DNA control-region sequencing and microsatellite genotyping. We found relatively high levels of genetic diversity, which is similar to values found in previous studies on Carpathian wolves from Poland and Slovakia, as well as to the long-lasting Dinaric-Balkan wolf population. We found no significant population structure in our study region, suggesting effective dispersal and admixture. Analysis of wolf-dog hybridisation using a Single Nucleotide Polymorphism panel optimised for hybrid detection revealed low rates of admixture between wolves and domestic dogs. Our results provide evidence for the existence of a genetically viable wolf population in the Romanian Carpathians. The genetic data obtained in this study may serve as valuable baseline information for the elaboration of monitoring standards and management plans for wolves in Romania.

Phylogeography, hybrid zones and contemporary species boundaries in the south-eastern Australian smooth frogs (Anura: Myobatrachidae: Geocrinia)

Paleo-climatic fluctuations have driven episodic changes in species distributions, providing opportunities for populations to diverge in isolation and hybridise following secondary contact. Studies of phylogeographic diversity and patterns of gene flow across hybrid zones can provide insight into contemporary species boundaries and help to inform taxonomic and conservation inferences. Here we explore geographic diversity within the acoustically divergent yet morphologically conserved south-eastern Australian smooth frog complex and assess gene flow across a narrow hybrid zone using mitochondrial nucleotide sequences and nuclear genome-wide single nucleotide polymorphisms. Our analyses reveal the presence of an evolutionarily distinct taxon restricted to the Otway Plains and Ranges, Victoria, which forms a narrow (9-30 km wide), spatiotemporally stable (>50 years) hybrid zone with Geocrinia laevis, which we describe herein as a new species.

Genomic signatures of ecological divergence between savanna and forest populations of a Neotropical tree

BACKGROUND AND AIMS

In eastern Neotropical South America, the Cerrado, a large savanna vegetation, and the Atlantic Forest harbour high biodiversity levels, and their habitats are rather different from each other. The biomes have intrinsic evolutionary relationships, with high lineage exchange that can be attributed, in part, to a large contact zone between them. The genomic study of ecotypes, i.e. populations adapted to divergent habitats, can be a model to study the genomic signatures of ecological divergence. Here, we investigated two ecotypes of the tree Plathymenia reticulata, one from the Cerrado and the other from the Atlantic Forest, which have a hybrid zone in the ecotonal zone of Atlantic Forest-Cerrado.

METHODS

The ecotypes were sampled in the two biomes and their ecotone. The evolutionary history of the divergence of the species was analysed with double-digest restriction site-associated DNA sequencing. The genetic structure and the genotypic composition of the hybrid zone were determined. Genotype-association analyses were performed, and the loci under putative selection and their functions were investigated.

KEY RESULTS

High divergence between the two ecotypes was found, and only early-generation hybrids were found in the hybrid zone, suggesting a partial reproductive barrier. Ancient introgression between the Cerrado and Atlantic Forest was not detected. The soil and climate were associated with genetic divergence in Plathymenia ecotypes and outlier loci were found to be associated with the stress response, with stomatal and root development and with reproduction.

CONCLUSIONS

The high genomic, ecological and morphophysiological divergence between ecotypes, coupled with partial reproductive isolation, indicate that the ecotypes represent two species and should be managed as different evolutionary lineages. We advise that the forest species should be re-evaluated and restated as vulnerable. Our results provide insights into the genomic mechanisms underlying the diversification of species across savanna and forest habitats and the evolutionary forces acting in the species diversification in the Neotropics.

Morphological differentiation in giant salamanders, Andrias japonicus, A. davidianus, and their hybrids (Urodela, Cryptobranchidae), and its taxonomic implications

For a long time, it has been debated whether the two giant salamanders, Andrias japonicus from Japan and A. davidianus from China, are conspecific or heterospecific. Morphological information about their diagnostic characteristics has been limited, without considering sexual dimorphism and/or body size variation. Recently, A. davidianus, which was introduced into Japan sometime in the past, has been found to hybridize with A. japonicus in situ. Taxonomic identification of individuals involved in this unusual breeding is made based on mitochondrial and nuclear DNA analyses. This identification method is time-consuming and costly. Thus, developing easier methods of identification, such as utilizing external morphological characteristics, is urgently needed. In this study, we verify previous descriptions showing that A. davidianus has a longer relative tail length than A. japonicus, and the tubercles on the lower jaw and throat were present in both sexes of A. davidianus. In addition, many head characteristics were found to be relatively larger in A. davidianus than in A. japonicus, which were new distinguishing characters. These morphological differences help support the idea that these are heterospecific lineages. In hybrids, relative values of head width and tail length were larger than those of A. japonicus, and the tubercles on their lower jaw and throat were present as in A. davidianus, suggesting that the hybrids and A. davidianus are distinguishable from A. japonicus.

Disentangling the causes of high polymorphism sharing in sympatric Petunia species from subtropical highland grasslands: insights from nuclear diversity

Genetic polymorphism sharing between closely related and sympatric plant species could result from common ancestry, ancient or recent hybridization. Here we analyzed four Petunia species from the subtropical highland grasslands in southern South America based on nuclear diversity to disentangle the causes of high polymorphism sharing between them. We genotyped microsatellite loci, employed population genetic methods to estimate variability, species limits, and ancient and recent gene flow, and assigned individuals to genetic and taxonomic groups. Finally, we modeled evolutionary processes to determine the impact of Quaternary climate changes on species phylogenetic relationships. Our results indicated that genetic diversity was strongly influenced by expansion and habitat fragmentation during the Quaternary cycles. The extensive polymorphism sharing is mainly due to species' common ancestry, and we did not discard ancient hybridization. Nowadays, niche differentiation is the primary driver for maintaining genetic and morphological limits between the four analysed Petunia species and there is no recent gene flow between them.

Demography-driven and adaptive introgression in a hybrid zone of the Armeria syngameon

Syngameons represent networks of otherwise distinct species connected by limited gene exchange. Although most studies have focused on how species maintain their cohesiveness despite gene flow, there are additional relevant questions regarding the evolutionary dynamics of syngameons and their drivers, as well as the success of their members and the network as a whole. Using a ddRADseq approach, we analysed the genetic structure, genomic clines and demographic history of a coastal hybrid zone involving two species of the Armeria (Plumbaginaceae) syngameon in southern Spain. We inferred that a peripheral population of the sand dune-adapted A. pungens diverged from the rest of the conspecific populations and subsequently hybridized with a locally more abundant pinewood congener, A. macrophylla. Both species display extensive plastid DNA haplotype sharing. Genomic cline analysis identified bidirectional introgression, but more outlier loci with excess A. pungens than A. macrophylla ancestry, suggesting the possibility of selection for A. pungens alleles. This is consistent with the finding that the A. pungens phenotype is selected for in open habitats, and with the strong correlation found between ancestry and phenotype. Taken together, our analyses suggest an intriguing scenario in which bidirectional introgression may, on the one hand, help to avoid reduced levels of genetic diversity due to the small size and isolated location of the A. pungens range-edge population, thereby minimizing demographic risks of stochastic extinction. On the other hand, the data also suggest that introgression into A. macrophylla may allow individuals to grow in open, highly irradiated, deep sandy, salt-exposed habitats.

Revision of Litoria rothii (Anura: Pelodryadidae) from northern Australia

Litoria rothii is a widespread pelodryadid frog with a charismatic laughing advertisement call, distributed across the Australian Monsoon Tropics and southern New Guinea. Given its large distribution spanning well-known biogeographic barriers, variation in male advertisement calls and the prevalence of unresolved species complexes in the Australian frog fauna, we examine the genetic, morphological and acoustic diversity in the species from across its range. Our analyses reveal the presence of a previously unrecognised species in western parts of the range of L. rothii sensu lato, which we describe herein as a new species. Litoria ridibunda sp. nov. is distinguished from L. rothii on the basis of paraphyly of nuclear gene trees with L. everetti from Indonesia, colour patterns on the posterior thigh and male advertisement calls. Compared to L. rothii, the new species has a less contrasting pattern on the posterior thigh and a male advertisement call with a greater number of notes per call and a greater call duration. In particular, the magnitude of call differences between the species is highest where the ranges of the two species are in proximity in north-western Queensland. Our study further emphasises the undiagnosed diversity that remains in Australian frogs, even in relatively large, charismatic, frequently encountered species that often share human dwellings.

Estimating admixture pedigrees of recent hybrids without a contiguous reference genome

The genome of recently admixed individuals or hybrids has characteristic genetic patterns that can be used to learn about their recent admixture history. One of these are patterns of interancestry heterozygosity, which can be inferred from SNP data from either called genotypes or genotype likelihoods, without the need for information on genomic location. This makes them applicable to a wide range of data that are often used in evolutionary and conservation genomic studies, such as low-depth sequencing mapped to scaffolds and reduced representation sequencing. Here we implement maximum likelihood estimation of interancestry heterozygosity patterns using two complementary models. We furthermore develop apoh (Admixture Pedigrees of Hybrids), a software that uses estimates of paired ancestry proportions to detect recently admixed individuals or hybrids, and to suggest possible admixture pedigrees. It furthermore calculates several hybrid indices that make it easier to identify and rank possible admixture pedigrees that could give rise to the estimated patterns. We implemented apoh both as a command line tool and as a Graphical User Interface that allows the user to automatically and interactively explore, rank and visualize compatible recent admixture pedigrees, and calculate the different summary indices. We validate the performance of the method using admixed family trios from the 1000 Genomes Project. In addition, we show its applicability on identifying recent hybrids from RAD-seq data of Grant's gazelle (Nanger granti and Nanger petersii) and whole genome low-depth data of waterbuck (Kobus ellipsiprymnus) which shows complex admixture of up to four populations.

Seasonal differences in the timing of flight between the invasive winter moth and native Bruce spanworm promotes reproductive isolation

The European winter moth, Operophtera brumata L. (Lepidoptera: Geometridae), was accidentally introduced to North America on at least 4 separate occasions, where it has been hybridizing with the native Bruce spanworm, O. bruceata Hulst, at rates up to 10% per year. Both species are known to respond to the same sex pheromones and to produce viable offspring, but whether they differ in the seasonal timing of their mating flights is unknown. Therefore, we collected adult male moths weekly along 2 transects in the northeastern United States and genotyped individuals using polymorphic microsatellite markers as males of these 2 species cannot be differentiated morphologically. Along each transect, we then estimated the cumulative proportions (i.e., the number of individuals out of the total collected) of each species on each calendar day. Our results indicate that there are significant differences between the species regarding their seasonal timing of flight, and these allochronic differences likely are acting to promote reproductive isolation between these 2 species. Lastly, our results suggest that the later flight observed by winter moth compared to Bruce spanworm may be limiting its inland spread in the northeastern United States because of increased exposure to extreme winter events.

Hybridization boosters diversification in a Neotropical Bulbophyllum (Orchidaceae) group

Genetic data shows that cryptic hybrids are more common than previously thought and that hybridization and introgression are widespread processes. Regardless, studies on hybridization are scarce for the highly speciose Bulbophyllum. The genus presents more than 2,200 species and many examples of recent radiations, in which hybridization is expected to be frequent. Currently, only four natural Bulbophyllum hybrids are recognized, all of them recently described based on morphological evidence. Here we test whether genomic evidence supports the hybrid status of two Neotropical Bulbophyllum species, while also evaluating the impact of this phenomenon on the genomes of the putative parental species. We also assess if there is evidence of hybridization among B. involutum and B. exaltatum, sister species that diverged recently. We leverage the power of next-generation sequence data, associated with model-based analysis for three systems putatively constituted by two parental species and one hybrid. All taxa belong to the Neotropical B. sect. Didactyle clade. We found evidence of hybridization in all studied systems. Despite the occurrence of hybridization, there are no signs of backcrossing. Because of the high propensity of hybridization across many taxa, the common occurrence of hybridization during the evolutionary history of B. sect. Didactyle means it is time to account for and examine its evolutionary role in these orchids.

Detecting introgressive hybridization to maintain genetic integrity in endangered large waterbird: a case study in milky stork

Hybridization between milky stork (Mycteria cinerea) and painted stork (M. leucocephala) occurs frequently in captivity. Dusit Zoo is a captive breeding facility where storks with phenotypically ambiguous patterns have recently been observed, and their status remaining inconclusive. Here, we used a combination of phenotypic characters and genetic markers (cytochrome b and 14 microsatellite markers) to distinguish and identify hybrids from the two parental species (n = 114). Haplotype analysis revealed asymmetric mtDNA introgression from M. cinerea to M. leucocephala, with twelve morphologically classified M. leucocephala individuals carrying heterospecific mtDNA. Comprehensive biparental genetic assessments identified 33% of all three genetic clusters as admixed individuals, of which most were either F2 hybrids, backcrosses with M. leucocephala, or hybrids of unknown generation, implying weak premating isolation with the absence of intrinsic postzygotic isolation between parentals. Morphological analysis demonstrated that the absence or indistinctness of a black bar across the breast is the most noticeable trait to identify these hybrids. The endangered M. cinerea was found to have genomic contamination from M. leucocephala and vice versa, with at least 41 hybrid individuals being identified. These findings provide critical information for detecting hybrids and identifying suitable breeding stocks with genetic purity for future reintroduction and conservation management.

A community genomics approach to natural hybridization

Hybridization is a complicated, oft-misunderstood process. Once deemed unnatural and uncommon, hybridization is now recognized as ubiquitous among species. But hybridization rates within and among communities are poorly understood despite the relevance to ecology, evolution and conservation. To clarify, we examined hybridization across 75 freshwater fish communities within the Ozarks of the North American Interior Highlands (USA) by single nucleotide polymorphism (SNP) genotyping 33 species (N = 2865 individuals; double-digest restriction site-associated DNA sequencing (ddRAD)). We found evidence of hybridization (70 putative hybrids; 2.4% of individuals) among 18 species-pairs involving 73% (24/33) of study species, with the majority being concentrated within one family (Leuciscidae/minnows; 15 species; 66 hybrids). Interspecific genetic exchange-or introgression-was evident from 24 backcrossed individuals (10/18 species-pairs). Hybrids occurred within 42 of 75 communities (56%). Four selected environmental variables (species richness, protected area extent, precipitation (May and annually)) exhibited 73-78% accuracy in predicting hybrid occurrence via random forest classification. Our community-level assessment identified hybridization as spatially widespread and environmentally dependent (albeit predominantly within one diverse, omnipresent family). Our approach provides a more holistic survey of natural hybridization by testing a wide range of species-pairs, thus contrasting with more conventional evaluations.

Ocean warming favours a northern Argyrosomus species over its southern congener, whereas preliminary metabolic evidence suggests that hybridization may promote their adaptation

Anthropogenic-induced climate change is having profound impacts on aquatic ecosystems, and the resilience of fish populations will be determined by their response to these impacts. The northern Namibian coast is an ocean warming hotspot, with temperatures rising faster than the global average. The rapid warming in Namibia has had considerable impacts on marine fauna, such as the southern extension of the distribution of Argyrosomus coronus from southern Angola into northern Namibian waters, where it now overlaps and hybridizes with the closely related Namibian species, A. inodorus. Understanding how these species (and their hybrids) perform at current and future temperatures is vital to optimize adaptive management for Argyrosomus species. Intermittent flow-through respirometry was used to quantify standard and maximum metabolic rates for Argyrosomus individuals across a range of temperatures. The modelled aerobic scope (AS) of A. inodorus was notably higher at cooler temperatures (12, 15, 18 and 21°C) compared with that of A. coronus, whereas the AS was similar at 24°C. Although only five hybrids were detected and three modelled, their AS was in the upper bounds of the models at 15, 18 and 24°C. These findings suggest that the warming conditions in northern Namibia may increasingly favour A. coronus and promote the poleward movement of the leading edge of their southern distribution. In contrast, the poor aerobic performance of both species at cold temperatures (12°C) suggests that the cold water associated with the permanent Lüderitz Upwelling Cell in the south may constrain both species to central Namibia. This is most concerning for A. inodorus because it may be subjected to a considerable coastal squeeze.

Citizen science facilitates first ever genetic detection of wolf-dog hybridization in Indian savannahs

Human demographic expansion has confined wildlife to fragmented habitats, often in proximity to human-modified landscapes. Such interfaces facilitate increased interactions between feral or domesticated animals and wildlife, posing a high risk to wild species. This is especially relevant for free-ranging dogs (Canis lupus familiaris) and wild canids like gray wolves (Canis lupus) and golden jackals (Canis aureus). Wolf-dog hybridization may lead to a significant reduction of specific adaptations in wolves that could result in the decline of wolf populations. Detection and genetic discrimination of hybrids between dogs and wolves are challenging because of their complex demographic history and close ancestry. Citizen scientists identified two phenotypically different-looking individuals and subsequently collected non-invasive samples that were used by geneticists to test wolf-dog hybridization. Genomic data from shed hair samples of suspected hybrid individuals using double-digest restriction-site-associated DNA (ddRAD) sequencing resulted in 698 single nucleotide polymorphism (SNP) markers. We investigated the genetic origin of these two individuals analyzed with genetically known dogs, wolves, and other canid species including jackals and dholes (Cuon alpinus). Our results provide the first genetic evidence of one F2 hybrid and the other individual could be a complex hybrid between dogs and wolves. Our results re-iterate the power of next-generation sequencing (NGS) for non-invasive samples as an efficient tool for detecting hybrids. Our results suggest the need for more robust monitoring of wolf populations and highlight the tremendous potential for collaborative approaches between citizens and conservation scientists to detect and monitor threats to biodiversity.

From genomics to integrative species delimitation? The case study of the Indo-Pacific Pocillopora corals

With the advent of genomics, sequencing thousands of loci from hundreds of individuals now appears feasible at reasonable costs, allowing complex phylogenies to be resolved. This is particularly relevant for cnidarians, for which insufficient data is available due to the small number of currently available markers and obscures species boundaries. Difficulties in inferring gene trees and morphological incongruences further blur the study and conservation of these organisms. Yet, can genomics alone be used to delimit species? Here, focusing on the coral genus Pocillopora, whose colonies play key roles in Indo-Pacific reef ecosystems but have challenged taxonomists for decades, we explored and discussed the usefulness of multiple criteria (genetics, morphology, biogeography and symbiosis ecology) to delimit species of this genus. Phylogenetic inferences, clustering approaches and species delimitation methods based on genome-wide single-nucleotide polymorphisms (SNP) were first used to resolve Pocillopora phylogeny and propose genomic species hypotheses from 356 colonies sampled across the Indo-Pacific (western Indian Ocean, tropical southwestern Pacific and south-east Polynesia). These species hypotheses were then compared to other lines of evidence based on genetic, morphology, biogeography and symbiont associations. Out of 21 species hypotheses delimited by genomics, 13 were strongly supported by all approaches, while six could represent either undescribed species or nominal species that have been synonymised incorrectly. Altogether, our results support (1) the obsolescence of macromorphology (i.e., overall colony and branches shape) but the relevance of micromorphology (i.e., corallite structures) to refine Pocillopora species boundaries, (2) the relevance of the mtORF (coupled with other markers in some cases) as a diagnostic marker of most species, (3) the requirement of molecular identification when species identity of colonies is absolutely necessary to interpret results, as morphology can blur species identification in the field, and (4) the need for a taxonomic revision of the genus Pocillopora. These results give new insights into the usefulness of multiple criteria for resolving Pocillopora, and more widely, scleractinian species boundaries, and will ultimately contribute to the taxonomic revision of this genus and the conservation of its species.

Hybridization in the Fringed Orchids: An Analysis of Species Boundaries in the Face of Gene Flow

Natural hybridization between closely related species in sympatry is an evolutionary process that is common in orchids. Once seen as a threat to parent species, interspecific genetic change is increasingly viewed as a source of novel variation in some ecological contexts. Terrestrial fringed orchids in the genus Platanthera contain several clades with high genetic compatibility among species and many putative hybrids. We used biallelic SNPs generated with 3RAD sequencing to study the hybrid complex formed from the parent species P. blephariglottis, P. ciliaris, and P. cristata with high resolution. The genetic structure and phylogenetic relationship of the hybrid complex revealed site-dependent gene flow between species. We documented extensive hybridization and cryptic hybrids in sympatric sites. Interspecific genetic exchange is particularly common between P. blephariglottis and P. ciliaris, with cryptic hybrids among putative P. ciliaris samples being more common than parental assignments in sympatric sites. Hybridization across the triad species complex can reticulate lineages and introduce adaptive alleles. Conversely, it can reduce diversification rates and introduce maladaptive alleles. Investigation into whether anthropogenic forces are eroding species boundaries, particularly the permeable P. blephariglottis and P. ciliaris boundary, is appropriate for conservation efforts.

Hybridization despite elaborate courtship behavior and female choice in Neotropical tree frogs

The mechanisms of hybridization can be elucidated by analyzing genotypes as well as phenotypes that could act as premating barriers, as the reproductive interactions among heterospecifics can alter the evolutionary history of species. In frogs, hybrids typically occur among species that reproduce explosively (in dense aggregations) with few opportunities for mate selection but are rare in species with elaborate courtship behaviors that may prevent erroneous mating. Using 21 microsatellite markers, we examined hybridization in the prolonged-breeding tree frogs Bokermannohyla ibitiguara and B. sazimai sampled within a contact zone in the Brazilian savanna (72 tadpoles; 74 adults). We also compared acoustic and morphological data. We confirmed both parental species genetically; STRUCTURE results confirmed 14 hybrids, 11 of which were second-generation according to NEWHYBRIDS, all with intermediate values of genetic dissimilarities compared to the parentals. Morphological and acoustic analyses revealed that hybrids showed variable but not necessarily intermediate phenotypes. Moreover, 2 hybrids exhibited call types different from parentals. The reproduction of B. ibitiguara involves territorial and aggressive males, elaborate courtships with acoustic and tactile stimuli, choosy females, and opportunistic strategies. Our study uncovers a rare case of viable hybridization among closely related frogs with such a combination of complex courtship behaviors and mate choice. We discuss the likely directionality and mechanisms behind this phenomenon, and highlight the importance of investigating hybridization even in species that show elaborate reproduction and female choice to advance our understanding of animal diversification.

Homoploid hybrids are common but evolutionary dead ends, whereas polyploidy is not linked to hybridization in a group of Pyrenean saxifrages

Hybridization and polyploidy are major forces in plant evolution. Homoploid hybridization can generate new species via hybrid speciation, or modify extant evolutionary lineages through introgression. Polyploidy enables instantaneous reproductive isolation from the parental lineage(s) and is often coupled with evolutionary innovations, especially when linked to hybridization. While allopolyploidy is a well-known and common mechanism of plant speciation, the evolutionary role of autopolyploidy might have been underestimated. Here, we studied the saxifrages of Saxifraga subsection Saxifraga in the Pyrenees, which easily hybridise and include polyploid populations of uncertain origin, as a model to unravel evolutionary consequences and origin of hybridization and polyploidy. Additionally, we investigate the phylogenetic relationship between the two subspecies of the endemic S. pubescens to ascertain whether they should rather be treated as different species. For these purposes, we combined ploidy-informed restriction associated DNA analyses, plastid DNA sequences and morphological data on a comprehensive population sample of seven species. Our results unravel multiple homoploid hybridization events at the diploid level between different species pairs, but with limited evolutionary impact. The ploidy-informed analyses reveal that all tetraploid populations detected in the present study belong to the widespread alpine species S. moschata. Although of autopolyploid origin, they are to some extent morphologically differentiated and underwent a different evolutionary pathway than their diploid parent. However, the high plastid DNA diversity and the internal structure within eastern and western population groups suggest multiple origins of the polyploids. Finally, our phylogenetic analyses show that S. pubescens and S. iratiana are clearly not sister lineages, and should consequently be considered as independent species.

Introgressive hybridization in the west Pacific pen shells (genus Atrina): Restricted interspecies gene flow within the genome

A compelling interest in marine biology is to elucidate how species boundaries between sympatric free-spawning marine invertebrates such as bivalve molluscs are maintained in the face of potential hybridization. Hybrid zones provide the natural resources for us to study the underlying genetic mechanisms of reproductive isolation between hybridizing species. Against this backdrop, we examined the occurrence of introgressive hybridization (introgression) between two bivalves distributed in the western Pacific margin, Atrina japonica and Atrina lischkeana, based on single-nucleotide polymorphisms (SNPs) derived from restriction site-associated DNA sequencing. Using 1066 ancestry-informative SNP sites, we also investigated the extent of introgression within the genome to search for SNP sites with reduced interspecies gene flow. A series of our individual-level clustering analyses including the principal component analysis, Bayesian model-based clustering, and triangle plotting based on ancestry-heterozygosity relationships for an admixed population sample from the Seto Inland Sea (Japan) consistently suggested the presence of specimens with varying degrees of genomic admixture, thereby implying that the two species are not completely isolated. The Bayesian genomic cline analysis identified 10 SNP sites with reduced introgression, each of which was located within a genic region or an intergenic region physically close to a functional gene. No, or very few, heterozygotes were observed at these sites in the hybrid zone, suggesting that selection acts against heterozygotes. Accordingly, we raised the possibility that the SNP sites are within genomic regions that are incompatible between the two species. Our finding of restricted interspecies gene flow at certain genomic regions gives new insight into the maintenance of species boundaries in hybridizing broadcast-spawning molluscs.

Hybridization between closely related songbirds is related to human habitat disturbance

Human habitat disturbances can promote hybridization between closely related, but typically reproductively isolated, species. We explored whether human habitat disturbances are related to hybridization between two closely related songbirds, black-capped and mountain chickadees, using both genomic and citizen science data sets. First, we genotyped 409 individuals from across both species' ranges using reduced-representation genome sequencing and compared measures of genetic admixture to a composite measure of human landscape disturbance. Then, using eBird observations, we compared human landscape disturbance values for sites where phenotypically diagnosed hybrids were observed to locations where either parental species was observed to determine whether hybrid chickadees are reported in more disturbed areas. We found that hybridization between black-capped and mountain chickadees positively correlates with human habitat disturbances. From genomic data, we found that (1) hybrid index (HI) significantly increased with habitat disturbance, (2) more hybrids were sampled in disturbed habitats, (3) mean HIs were higher in disturbed habitats versus wild habitats, and (4) hybrids were detected in habitats with significantly higher disturbance values than parentals. Using eBird data, we found that both hybrid and black-capped chickadees were significantly more disturbance-associated than mountain chickadees. Surprisingly, we found that nearly every black-capped chickadee we sampled contained some proportion of hybrid ancestry, while we detected very few mountain chickadee backcrosses. Our results highlight that hybridization between black-capped and mountain chickadees is widespread, but initial hybridization is rare (few F1s were detected). We conclude that human habitat disturbances can erode pre-zygotic reproductive barriers between chickadees and that post-zygotic isolation is incomplete. Understanding what becomes of recently hybridizing species following large-scale habitat disturbances is a new, but pressing, consideration for successfully preserving genetic biodiversity in a rapidly changing world.

Appraising the Genetic Makeup of an Allochthonous Southern Pike Population: An Opportunity to Predict the Evolution of Introgressive Hybridization in Isolated Populations?

Biological invasions are a major threat to the conservation of biodiversity, as invasive species affect native biota through competition, predation, pathogen introduction, habitat alteration, and hybridisation. The present study focuses on a southern pike population, Esox cisalpinus (Teleostei: Esocidae), that has been introduced outside the species' native range. Using microsatellite markers, this study's objective was to gather baseline genetic information and assess the presence of hybrids between this species and E. lucius in the introduced population. The resulting estimates of genetic diversity and effective population size are comparable to those observed in the species' native range. Although different methods yield contrasting and uncertain evidence regarding introgressive hybridization, the presence of late-generation hybrids cannot be completely ruled out. Large numbers of breeders as well as multiple introductions of genetically divergent cohorts and introgressive hybridisation may explain the high genetic diversity of this recently introduced southern pike population. The present study issues a warning that the conservation of southern pike' introgressive hybridisation between northern and southern pike might be underestimated. The genetic information gathered herein may unravel the origin, number of introduction events, and evolutionary trajectory of the introduced population. This information may help us understand the evolution of introgressive hybridisation in the southern pike's native areas.

Evolution of a hybrid zone of two willow species (Salix L.) in the European Alps analyzed by RAD-seq and morphometrics

Natural hybridization of plants can result in many outcomes with several evolutionary consequences, such as hybrid speciation and introgression. Natural hybrid zones can arise in mountain systems as a result of fluctuating climate during the exchange of glacial and interglacial periods, where species retract and expand their territories, resulting in secondary contacts. Willows are a large genus of woody plants with an immense capability of interspecific crossing. In this study, the sympatric area of two diploid sister species, S. foetida and S. waldsteiniana in the eastern European Alps, was investigated to study the genomic structure of populations within and outside their contact zone and to analyze congruence of morphological phenotypes with genetic data. Eleven populations of the two species were sampled across the Alps and examined using phylogenetic network and population genetic structure analyses of RAD Seq data and morphometric analyses of leaves. The results showed that a homoploid hybrid zone between the two species was established within their sympatric area. Patterns of genetic admixture in homoploid hybrids indicated introgression with asymmetric backcrossing to not only one of the parental species but also one hybrid population forming a separate lineage. The lack of F1 hybrids indicated a long-term persistence of the hybrid populations. Insignificant isolation by distance suggests that gene flow can act over large geographical scales. Morphometric characteristics of hybrids supported the molecular data and clearly separated populations of the parental species, but showed intermediacy in the hybrid zone populations with a bias toward S. waldsteiniana. The homoploid hybrid zone might have been established via secondary contact hybridization, and its establishment was fostered by the low genetic divergence of parental species and a lack of strong intrinsic crossing barriers. Incomplete ecological separation and the ability of long-distance dispersal of willows could have contributed to the spatial expansion of the hybrid zone.

Floral scent divergence across an elevational hybrid zone with varying pollinators

Divergence in floral traits attractive to different pollinators can promote reproductive isolation in related species. When isolation is incomplete, hybridization may occur, which offers the opportunity to explore mechanisms underlying reproductive isolation. Recent work suggests that divergence in floral scent may frequently contribute to reproductive barriers, although such divergence has seldom been examined in species with generalized pollination. Here, we used two closely related Penstemon species, P. newberryi and P. davidsonii, and their natural hybrids from an elevational gradient with pollinator communities that are predicted to vary in their reliance on floral scent (i.e., primarily hummingbirds at low elevation vs. bees at high elevation). The species vary in a suite of floral traits, but scent is uncharacterized. To address whether scent varies along elevation and potentially contributes to reproductive isolation, we genetically characterized individuals collected at field and identified whether they were parental species or hybrids. We then characterized scent amount and composition. Although the parental species had similar total emissions, some scent characteristics (i.e., scent composition, aromatic emission) diverged between them and may contribute to their isolation. However, the species emitted similar compound sets which could explain hybridization in the contact area. Hybrids were similar to the parents for most scent traits, suggesting that their floral scent would not provide a strong barrier to backcrossing. Our study suggests floral scent may be a trait contributing to species boundaries even in plants with generalized pollination, and reinforces the idea that evolutionary pollinator transitions may involve changes in multiple floral traits.

Is there hybridization between diploid and tetraploid Euphrasia in a secondary contact zone?

PREMISE

Strong postzygotic reproductive isolating barriers are usually expected to limit the extent of natural hybridization between species with contrasting ploidy. However, genomic sequencing has revealed previously overlooked examples of natural cross-ploidy hybridization in some flowering plant genera, suggesting that the phenomenon may be more common than once thought. We investigated potential cross-ploidy hybridization in British eyebrights (Euphrasia, Orobanchaceae), a group from which 13 putative cross-ploidy hybrid combinations have been reported based on morphology.

METHODS

We analyzed a contact zone between diploid Euphrasia rostkoviana and tetraploid E. arctica in Wales. We sequenced part of the internal transcribed spacer (ITS) of nuclear ribosomal DNA and used genotyping by sequencing (GBS) to look for evidence of cross-ploidy hybridization and introgression.

RESULTS

Common variant sites in the ITS region were fixed between diploids and tetraploids, indicating a strong barrier to hybridization. Clustering analyses of 356 single-nucleotide polymorphisms (SNPs) generated using GBS clearly separated samples by ploidy and revealed strong genetic structure (F_ST  = 0.44). However, the F_ST distribution across all SNPs was bimodal, indicating potential differential selection on loci between diploids and tetraploids. Demographic inference suggested potential gene flow, limited to around one or fewer migrants per generation.

CONCLUSIONS

Our results suggest that recent cross-ploidy hybridization is rare or absent in a site of secondary contact in Euphrasia. While a strong ploidy barrier prevents hybridization over ecological timescales, such hybrids may form in stable populations over evolutionary timescales, potentially allowing cross-ploidy introgression to take place.

Reproductive isolation via divergent genital morphology due to cascade reinforcement in Ohomopterus ground beetles

Secondary contact between incipient species and selection against maladaptive hybridization can drive reinforcement between populations in contact and result in reproductive character displacement (RCD). Resultant divergence in mating traits within a species may generate downstream reproductive isolation between populations with displaced and non-displaced traits, referred to as the cascade reinforcement hypothesis. We examined this hypothesis using three allopatric populations of the ground beetle Carabus maiyasanus with a genital lock-and-key system. This species shows RCD in male and female genital morphologies in populations in contact with the sister species C. iwawakianus. In a reciprocal mating experiment using three allopatric populations with differences in male and female genital sizes, insemination failure increased as the difference in genital size increased. Based on the reproductive isolation index, insemination failure was the major postmating-prezygotic isolation barrier, at least in one population pair with comparable total isolation to those of other species pairs. By contrast, there was only incomplete premating isolation among populations. These results suggest that RCD in genital morphologies drives incipient allopatric speciation, supporting the cascade reinforcement hypothesis. These findings provide insight into the roles of interspecific interactions and subsequent trait diversification in speciation processes.

Hexabothriidae and Monocotylidae (Monogenoidea) from the gills of neonate hammerhead sharks (Sphyrnidae) Sphyrna gilberti, Sphyrna lewini and their hybrids from the western North Atlantic Ocean

Neonates of hammerhead sharks (Sphyrnidae), Sphyrna lewini (Griffith and Smith, 1834), the sympatric cryptic species, Sphyrna gilberti Quattro et al., 2013, and their hybrids were captured in the western North Atlantic, along the coast of South Carolina, USA, between 2018 and 2019 and examined for gill monogenoids. Parasites were identified and redescribed from the gills of 79 neonates, and DNA sequences from partial fragments of the nuclear 28S ribosomal RNA (rDNA) and cytochrome c oxidase I mitochondrial DNA (COI) genes were generated to confirm species identifications. Three species of monogenoids from Hexabothriidae Price, 1942 and Monocotylidae Taschenberg, 1879 were determined and redescribed. Two species of Hexabothriidae, Erpocotyle microstoma (Brooks, 1934) and Erpocotyle sphyrnae (MacCallum, 1931), infecting both species of Sphyrna and hybrids; and 1 species of Monocotylidae, Loimosina wilsoni Manter, 1944, infecting only S. lewini and hybrids. Loimosina wilsoni 28S rDNA sequences matched those of Loimosina sp. from the southern coast of Brazil. Based on limited morphological analysis, Loimosina parawilsoni is likely a junior synonym of L. wilsoni. This is the first taxonomic study of monogenoids infecting S. gilberti and hybrids of S. gilberti and S. lewini.

Real-time geographic settling of a hybrid zone between the invasive winter moth (Operophtera brumata L.) and the native Bruce spanworm (O. bruceata Hulst)

Hybridization plays an important and underappreciated role in shaping the evolutionary trajectories of species. Following the introduction of a non-native organism to a novel habitat, hybridization with a native congener may affect the probability of establishment of the introduced species. In most documented cases of hybridization between a native and a non-native species, a mosaic hybrid zone is formed, with hybridization occurring heterogeneously across the landscape. In contrast, most naturally occurring hybrid zones are clinal in structure. Here, we report on a long-term microsatellite data set that monitored hybridization between the invasive winter moth, Operophtera brumata (Lepidoptera: Geometridae), and the native Bruce spanworm, O. bruceata, over a 12-year period. Our results document one of the first examples of the real-time formation and geographic settling of a clinal hybrid zone. In addition, by comparing one transect in Massachusetts where extreme winter cold temperatures have been hypothesized to restrict the distribution of winter moth, and one in coastal Connecticut, where winter temperatures are moderated by Long Island Sound, we found that the location of the hybrid zone appeared to be independent of environmental variables and maintained under a tension model wherein the stability of the hybrid zone was constrained by population density, reduced hybrid fitness, and low dispersal rates. Documenting the formation of a contemporary clinal hybrid zone may provide important insights into the factors that shaped other well-established hybrid zones.

Contrasting levels of hybridization across the two contact zones between two hedgehog species revealed by genome-wide SNP data

Hybridization and introgression have played important roles in the history of various species, including lineage diversification and the evolution of adaptive traits. Hybridization can accelerate the development of reproductive isolation between diverging species, and thus valuable insight into the evolution of reproductive barrier formation may be gained by studying secondary contact zones. Hedgehogs of the genus Erinaceus, which are insectivores sensitive to changes in climate, are a pioneer model in Pleistocene phylogeography. The present study provides the first genome-wide SNP data regarding the Erinaceus hedgehogs species complex, offering a unique comparison of two secondary contact zones between Erinaceus europaeus and E. roumanicus. Results confirmed diversification of the genus during the Pleistocene period, and detected a new refugial lineage of E. roumanicus outside the Mediterranean basin, most likely in the Ponto-Caspian region. In the Central European zone, the level of hybridization was low, whereas in the Russian-Baltic zone, both species hybridise extensively. Asymmetrical gene flow from E. europaeus to E. roumanicus suggests that reproductive isolation varies according to the direction of the crosses in the hybrid zones. However, no loci with significantly different patterns of introgression were detected. Markedly different pre- and post-zygotic barriers, and thus diverse modes of species boundary maintenance in the two contact zones, likely exist. This pattern is probably a consequence of the different age and thus of the different stage of evolution of reproductive isolating mechanisms in each hybrid zone.

Unusual patterns of hybridization involving two alpine Salvia species: Absence of both F1 and backcrossed hybrids

Natural hybridization plays an important role in speciation; however, we still know little about the mechanisms underlying the early stages of hybrid speciation. Hybrid zones are commonly dominated by F₁s, or backcrosses, which impedes further speciation. In the present study, morphological traits and double digest restriction-site associated DNA sequencing (ddRAD-seq) data have been used to confirm natural hybridization between Salvia flava and S. castanea, the first case of identification of natural hybridization using combined phenotypic and molecular evidence in the East Asian clade of Salvia. We further examined several reproductive barriers in both pre-zygotic and post-zygotic reproductive stages to clarify the causes and consequences of the hybridization pattern. Our results revealed that reproductive isolation between the two species was strong despite the occurrence of hybridization. Interestingly, we found that most of the hybrids were likely to be F₂s. This is a very unusual pattern of hybridization, and has rarely been reported before. The prevalence of geitonogamy within these self-compatible hybrids due to short distance foraging by pollinators might explain the origin of this unusual pattern. F₂s can self-breed and develop further, therefore, we might be witnessing the early stages of hybrid speciation. Our study provides a new case for understanding the diversification of plants on the Qinghai-Tibet Plateau.

Strong bidirectional gene flow between fish lineages separated for over 100,000 years

Restoring levels of genetic diversity in small and declining populations is increasingly being considered in biodiversity conservation. Evidence-based genetic management requires assessment of risks and benefits of crossing populations. Because risks are challenging to assess experimentally, e.g. through multi-generational crosses, decision-support approaches utilize proxy risk factors such as time since separation of lineages. However, the paucity of empirical datasets on fitness consequences of longer separation times tends to favour crossing lineages with conservatively short separations, restricting wildlife managers’ options. Here, we assessed the genetic outcomes of interbreeding in the wild between lineages of a threatened Australian freshwater fish (Macquarie perch) separated by an estimated 119,000–385,000 years of evolution in distinct environments. Fish belonging to the Murray-Darling Basin (MDB) lineage escaped from Cataract Dam—into which they were translocated in ~ 1915—into the Cataract River, where they interbred with the local Hawkesbury-Nepean Basin (HNB) lineage. Analyses of reduced-representation genomic data revealed no evidence of genetic incompatibilities during interbreeding of the two lineages in the Cataract River: assignment to genotypic clusters indicated a spectrum of hybrid types including second generation hybrids and backcrosses to both parental lineages. Thus, no adverse effects were detected from genetic mixing of populations separated by > 100,000 years. We are not advocating purposely crossing the two lineages for management purposes under present cost–benefit considerations, because there are currently sufficient intra-lineage source populations to beneficially mix. Instead, this study presents a useful calibration point: two morphologically different lineages evolved in different habitats for 119,000–385,000 years can successfully interbreed.

Hybridization is strongly constrained by salinity during secondary contact between silverside fishes (Odontesthes, Atheriniformes)

This study investigates a contact zone between two silverside fish species (marine Odontesthes argentinensis and freshwater O. bonariensis) in the estuarine Mar Chiquita lagoon along the Atlantic coast in Argentina (MChL), in which intermediate morphs had been reported. It has been suggested that admixture and introgression occur in MChL between these two species, but direct genetic evidence is lacking. Leveraging samples collected over several years (n = 676), we document the spatial distribution of both species and intermediate morphs within this habitat and collect landmark-based morphometric and multilocus genetic data (9876 loci for n = 110 individuals) to test the hypothesis of hybridization. Our analysis unambiguously characterizes intermediate morphs as F1 or F2 hybrids. We show that the low frequency of hybrid individuals in MChL may be explained by uneven abundance of parental species, which in turn are strongly affected by water salinity, limiting the size of the contact zone. Although hybrids seem to be fertile, their fitness may be reduced by external and intrinsic factors that may limit their success and suggest that this is an unstable hybrid zone. Genetic distinctiveness of both parental species is strongly supported by genome-wide data, explaining a known pattern of mitonuclear discordance as a consequence of hybridization followed by mitochondrial introgression. A clear signature of population genetic structure was detected in O. argentinensis, distinguishing MChL residents from marine populations of this species, that also was supported by distinctive morphometric features among these groups. Previous hypotheses of speciation in these fishes are discussed in the light of the new findings.

Evolutionary processes in an undescribed eucalypt: implications for the translocation of a critically endangered species

BACKGROUND AND AIMS

Knowledge of the evolutionary processes responsible for the distribution of threatened and highly localized species is important for their conservation. Population genomics can provide insights into evolutionary processes to inform management practices, including the translocation of threatened plant species. In this study, we focus on a critically endangered eucalypt, Eucalyptus sp. Cattai, which is restricted to a 40-km2 area of Sydney, Australia, and is threatened by increased urbanization. Eucalyptus sp. Cattai has yet to be formally described in part due to its suspected hybrid origin. Here, we examined evolutionary processes and species boundaries in E. sp. Cattai to determine whether translocation was warranted.

METHODS

We used genome-wide scans to investigate the evolutionary relationships of E. sp. Cattai with related species, and to assess levels of genetic health and admixture. Morphological trait and genomic data were obtained from seedlings of E. sp. Cattai propagated in a common garden to assess their genetic provenance and hybrid status.

KEY RESULTS

All analyses revealed that E. sp. Cattai was strongly supported as a distinct species. Genetic diversity varied across populations, and clonality was unexpectedly high. Interspecific hybridization was detected, and was more prevalent in seedlings compared to in situ adult plants, indicating that post-zygotic barriers may restrict the establishment of hybrids.

CONCLUSIONS

Multiple evolutionary processes (e.g. hybridization and clonality) can operate within one rare and restricted species. Insights regarding evolutionary processes from our study were used to assist with the translocation of genetically 'pure' and healthy ex situ seedlings to nearby suitable habitat. Our findings demonstrate that it is vital to provide an understanding of evolutionary relationships and processes with an examination of population genomics in the design and implementation of an effective translocation strategy.

Population genomics and sexual signals support reproductive character displacement in Uperoleia (Anura: Myobatrachidae) in a contact zone

When closely related species come into contact via range expansion, both may experience reduced fitness as a result of the interaction. Selection is expected to favour traits that minimize costly interspecies reproductive interactions (such as mismating) via a phenomenon called reproductive character displacement (RCD). Research on RCD frequently assumes secondary contact between species, but the geographical history of species interactions is often unknown. Population genomic data permit tests of geographical hypotheses about species origins and secondary contact through range expansion. We used population genomic data from single nucleotide polymorphisms (SNPs), mitochondrial sequence data, advertisement call data and morphological data to investigate a species complex of toadlets (Uperoleia borealis, U. crassa, U. inundata) from northern Australia. Although the three species of frogs were morphologically indistinguishable in our analysis, we determined that U. crassa and U. inundata form a single species (synonymized here) based on an absence of genomic divergence. SNP data identified the phylogeographical origin of U. crassa as the Top End, with subsequent westward invasion into the range of U. borealis in the Kimberley. We identified six F1 hybrids, all of which had the U. borealis mitochondrial haplotype, suggesting unidirectional hybridization. Consistent with the RCD hypothesis, U. borealis and U. crassa sexual signals differ more in sympatry than in allopatry. Hybrid males have intermediate calls, which probably reduces attractiveness to females. Integrating population genomic data, mitochondrial sequencing, morphology and behavioural approaches provides an unusually detailed collection of evidence for reproductive character displacement following range expansion and secondary contact.