Lost in parameter space: a road map for stacks

A spatial approach to jointly estimate Wright's neighborhood size and long-term effective population size

Spatially continuous patterns of genetic differentiation, which are common in nature, are often poorly described by existing population genetic theory or methods that assume either panmixia or discrete, clearly definable populations. There is therefore a need for statistical approaches in population genetics that can accommodate continuous geographic structure, and that ideally use georeferenced individuals as the unit of analysis, rather than populations or subpopulations. In addition, researchers are often interested in describing the diversity of a population distributed continuously in space; this diversity is intimately linked to both the dispersal potential and the population density of the organism. A statistical model that leverages information from patterns of isolation by distance to jointly infer parameters that control local demography (such as Wright's neighborhood size), and the long-term effective size (Ne) of a population would be useful. Here, we introduce such a model that uses individual-level pairwise genetic and geographic distances to infer Wright's neighborhood size and long-term Ne. We demonstrate the utility of our model by applying it to complex, forward-time demographic simulations as well as an empirical dataset of the two-form bumblebee (Bombus bifarius). The model performed well on simulated data relative to alternative approaches and produced reasonable empirical results given the natural history of bumblebees. The resulting inferences provide important insights into the population genetic dynamics of spatially structured populations.

Deciphering Pteronia's evolution in the Cape Floristic Region: A comprehensive study disputes polyploid deficiency and affirms diploid radiation

The Greater Cape Floristic Region (GCFR) is renowned for its exceptional biodiversity, accommodating over 11 000 plant species, notable degree of endemism, and substantial diversification within limited plant lineages, a phenomenon ascribed to historical radiation events. While both abiotic and biotic factors contribute to this diversification, comprehensive genomic alterations, recognized as pivotal in the diversification of angiosperms, are perceived as uncommon. This investigation focuses on the genus Pteronia, a prominent representative of the Asteraceae family in the GCFR. Employing NGS-based HybSeq and RADSeq methodologies, flow cytometry, karyology, and ecological modeling, we scrutinize the intricacies of its polyploid evolution. Phylogenetic reconstructions using 951 low-copy nuclear genes confirm Pteronia as a well-supported, distinct clade within the tribe Astereae. The ingroup displays a structure indicative of rapid radiation likely antedating polyploid establishment, with the two main groups demarcated by their presence or absence in the fynbos biome. Genome size analysis encompasses 1293 individuals across 347 populations, elucidating significant variation ranging from 6.1 to 34.2 pg (2C-value). Pteronia demonstrates substantially large genome sizes within Astereae and phanerophytes. Polyploidy is identified in 31% of the studied species, with four discerned ploidy levels (2x, 4x, 6x, 8x). Cytotypes exhibit marked distinctions in environmental traits, influencing their distribution across biomes and augmenting their niche differentiation. These revelations challenge the presumed scarcity of polyploidy in the Cape flora, underscoring the imperative need for detailed population studies. The intricate evolutionary history of Pteronia, characterized by recent polyploidy and genome size variation, contributes substantially to the comprehension of diversification patterns within the GCFR biodiversity hotspot.

Group vocal composition and decision-making during roost finding in Spix's disk-winged bats

Theoretical work suggests that having many informed individuals within social groups can promote efficient resource location. However, it may also give rise to group fragmentation if members fail to reach consensus on their direction of movement. In this study, we investigate whether the number of informed individuals, exemplified by bats emitting calls from different roosts, influences group cohesion in Spix's disk-winged bats (Thyroptera tricolor). Additionally, we explore the role of signal reliability, quantified through signalling rates, in group consensus on where to roost. These bats use contact calls to announce the location of a roost site and recruit conspecifics. The groups they form exhibit high levels of cohesion and consist of both vocal and non-vocal bats, with vocal behaviour being consistent over time. Our findings revealed that an increase in the number of roosts broadcasting calls is strongly associated with the likelihood of groups fragmenting among multiple roosts. Additionally, we found that a majority of group members enter the roost with higher calling rates. This phenomenon can mitigate the risk of group fragmentation, as bats emitting more calls may contribute to greater group consensus on roosting locations, thereby reducing the likelihood of individuals separating and enhancing overall group cohesion. Our results highlight the potential costs of having too many information producers for group coordination, despite their established role in finding critical resources. This article is part of the theme issue 'The power of sound: unravelling how acoustic communication shapes group dynamics'.

Calling to the collective: contact calling rates within groups of disc-winged bats do not vary by kinship or association

Many group-living animals coordinate social behaviours using contact calls, which can be produced for all group members or targeted at specific individuals. In the disc-winged bat, Thyroptera tricolor, group members use 'inquiry' and 'response' calls to coordinate daily movements into new roosts (furled leaves). Rates of both calls show consistent among-individual variation, but causes of within-individual variation remain unknown. Here, we tested whether disc-winged bats produce more contact calls towards group members with higher kinship or association. In 446 experimental trials, we recorded 139 random within-group pairs of one flying bat (producing inquiry calls for roost searching) and one roosting bat (producing response calls for roost advertising). Using generalized linear mixed-effect models (GLMM), we assessed how response and inquiry calling rates varied by sender, receiver, genetic kinship and co-roosting association rate. Calling rates varied consistently across senders but not by receiver. Response calling was influenced by inquiry calling rates, but neither calling rate was higher when the interacting pair had higher kinship or association. Rather than dyadic calling rates indicating within-group relationships, our findings are consistent with the hypothesis that bats produce contact calls to maintain contact with any or all individuals within a group while collectively searching for a new roost site. This article is part of the theme issue 'The power of sound: unravelling how acoustic communication shapes group dynamics'.

Genetic and morphological evidence support the specific status of the endemic Ericaandevalensis (Ericales, Ericaceae)

Assessing the taxonomic status of closely related taxa is crucial in plant systematics and can have important implications for conservation and human plant use. Ericaandevalensis Cabezudo & Rivera is a metallophyte endemic species from highly metal-polluted soils of SW Iberian Peninsula, an area with a mining history going back more than 5,000 years. Ericaandevalensis is closely related to Ericamackayana Bab., a northern Iberian species also present in western Ireland. The status of E.andevalensis as a species or subspecies subordinated to E.mackayana is subject to debate. Here, we assessed the genetic and phenotypic relationship between both species, including the population structure of E.andevalensis. We used high throughput sequencing to determine genome-wide Single Nucleotide Polymorphisms (SNPs), and morphometric analyses from 35 reproductive and vegetative traits. The morphological analysis showed at least eight characters that can discriminate the two species, from which ovary hairiness and the size of leaf glandular hairs were the most informative. Genetic analyses showed that each species formed a monophyletic cluster with full support, separated by an interspecific genetic distance >4-fold higher than intra-specific distance. Population genetic analyses of E.andevalensis shows that populations are highly structured, with the Portuguese one as the most isolated and less variable. These results support the recognition of E.andevalensis as a distinct species with a highly constrained ecological requirements and a narrow geographic distribution, but with a limited gene flow between populations. We discuss the implications of these outcomes in conservation policies and potential uses of E.andevalensis such as decontamination of polluted soils.

Montane diversification as a mechanism of speciation in neotropical butterflies

The mountains in the Atlantic Forest domain are environments that harbor a high biodiversity, including species adapted to colder climates that were probably influenced by the climatic variations of the Pleistocene. To understand the phylogeographic pattern and assess the taxonomic boundaries between two sister montane species, a genomic study of the butterflies Actinote mantiqueira and A. alalia (Nymphalidae: Acraeini) was conducted. Analyses based on partial sequences of the mitochondrial gene COI (barcode region) failed to recover any phylogenetic or genetic structure discriminating the two species or sampling localities. However, single nucleotide polymorphisms gathered using Genotyping-by-Sequencing provided a strong isolation pattern in all analyses (genetic distance, phylogenetic hypothesis, clustering analyses, and F ST statistics) which is consistent with morphology, separating all individuals of A. alalia from all populations of A. mantiqueira. The three sampled mountain ranges where A. mantiqueira populations occur-Serra do Mar, Serra da Mantiqueira, and Poços de Caldas Plateau-were identified as three isolated clusters. Paleoclimate simulations indicate that both species' distributions changed according to climatic oscillations in the Pleistocene period, with the two species potentially occurring in areas of lower altitude during glacial periods when compared to the interglacial periods (as the present). Besides, a potential path between their distribution through the Serra do Mar Mountain range was inferred. Therefore, the Pleistocene climatic fluctuation had a significant impact on the speciation process between A. alalia and A. mantiqueira, which was brought on by isolation at different mountain summits during interglacial periods, as shown by the modeled historical distribution and the observed genetic structure.

Different reference genomes determine different results: Comparing SNP calling in RAD-seq of Engelhardia roxburghiana using different reference genomes

Advances in next-generation sequencing (NGS) have significantly reduced the cost and improved the efficiency of obtaining single nucleotide polymorphism (SNP) markers, particularly through restriction site-associated DNA sequencing (RAD-seq). Meanwhile, the progression in whole genome sequencing has led to the utilization of an increasing number of reference genomes in SNP calling processes. This study utilized RAD-seq data from 242 individuals of Engelhardia roxburghiana, a tropical tree of the walnut family (Juglandaceae), with SNP calling conducted using the STACKS pipeline. We aimed to compare both reference-based approaches, namely, employing a closely related species as the reference genome versus the species itself as the reference genome, to evaluate their respective merits and limitations. Our findings indicate a substantial discrepancy in the number of obtained SNPs between using a closely related species as opposed to the species itself as reference genomes, the former yielded approximately an order of magnitude fewer SNPs compared to the latter. While the missing rate of individuals and sites of the final SNPs obtained in the two scenarios showed no significant difference. The results showed that using the reference genome of the species itself tends to be prioritized in RAD-seq studies. However, if this is unavailable, considering closely related genomes is feasible due to their wide applicability and low missing rate as alternatives. This study contributes to enrich the understanding of the impact of SNP acquisition when utilizing different reference genomes.

Next-generation data filtering in the genomics era

Genomic data are ubiquitous across disciplines, from agriculture to biodiversity, ecology, evolution and human health. However, these datasets often contain noise or errors and are missing information that can affect the accuracy and reliability of subsequent computational analyses and conclusions. A key step in genomic data analysis is filtering - removing sequencing bases, reads, genetic variants and/or individuals from a dataset - to improve data quality for downstream analyses. Researchers are confronted with a multitude of choices when filtering genomic data; they must choose which filters to apply and select appropriate thresholds. To help usher in the next generation of genomic data filtering, we review and suggest best practices to improve the implementation, reproducibility and reporting standards for filter types and thresholds commonly applied to genomic datasets. We focus mainly on filters for minor allele frequency, missing data per individual or per locus, linkage disequilibrium and Hardy-Weinberg deviations. Using simulated and empirical datasets, we illustrate the large effects of different filtering thresholds on common population genetics statistics, such as Tajima's D value, population differentiation (FST), nucleotide diversity (π) and effective population size (Ne).

Genomic insights into the conservation status of the Idle Crayfish Austropotamobius bihariensis Pârvulescu, 2019: low genetic diversity in the endemic crayfish species of the Apuseni Mountains

BACKGROUND

Biodiversity in freshwater ecosystems is declining due to an increased anthropogenic footprint. Freshwater crayfish are keystone species in freshwater ecosystems and play a crucial role in shaping the structure and function of their habitats. The Idle Crayfish Austropotamobius bihariensis is a native European species with a narrow distribution range, endemic to the Apuseni Mountains (Romania). Although its area is small, the populations are anthropogenically fragmented. In this context, the assessment of its conservation status is timely.

RESULTS

Using a reduced representation sequencing approach, we identified 4875 genomic SNPs from individuals belonging to 13 populations across the species distribution range. Subsequent population genomic analyses highlighted low heterozygosity levels, low number of private alleles and small effective population size. Our structuring analyses revealed that the genomic similarity of the populations is conserved within the river basins.

CONCLUSION

Genomic SNPs represented excellent tools to gain insights into intraspecific genomic diversity and population structure of the Idle Crayfish. Our study highlighted that the analysed populations are at risk due to their limited genetic diversity, which makes them extremely vulnerable to environmental alterations. Thus, our results emphasize the need for conservation measures and can be used as a baseline to establish species management programs.

Limited evidence for genetic differentiation or adaptation in two amphibian species across replicated rural-urban gradients

Urbanization leads to complex environmental changes and poses multiple challenges to organisms. Amphibians are highly susceptible to the effects of urbanization, with land use conversion, habitat destruction, and degradation ranked as the most significant threats. Consequently, amphibians are declining in urban areas, in both population numbers and abundance, however, the effect of urbanization on population genetic parameters remains unclear. Here, we studied the genomic response to urbanization in two widespread European species, the common toad Bufo bufo (26 localities, 480 individuals), and the smooth newt Lissotriton vulgaris (30 localities, 516 individuals) in three geographic regions: southern and northern Poland and southern Norway. We assessed genome-wide SNP variation using RADseq (ca. 42 and 552 thousand SNPs in toads and newts, respectively) and adaptively relevant major histocompatibility complex (MHC) class I and II genes. The results linked most of the genetic differentiation in both marker types to regional (latitudinal) effects, which also correspond to historical biogeography. Further, we did not find any association between genetic differentiation and level of urbanization at local scales for either species. However, urban smooth newts, but not toads, have lower levels of within-population genome-wide diversity, suggesting higher susceptibility to the negative effects of urbanization. A decreasing level of genetic diversity linked to increasing urbanization was also found for MHC II in smooth newts, while the relationship between MHC class I diversity and urbanization differed between geographic regions. We did not find any effects of urbanization on MHC diversity in the toad populations. Although two genetic environment association analyses of genome-wide data, LFMM and BayPass, revealed numerous (219 in B. bufo and 7040 in L. vulgaris) SNPs statistically associated with urbanization, we found a marked lack of repeatability between geographic regions, suggesting a complex and multifaceted response to natural selection elicited by life in the city.

First Evidence of Thalassochory in the Ficus Genus: Seed Dispersal Using the Kuroshio Oceanic Current

AIM

Plants distributed between southern Taiwan and the north of the Philippines are spread among numerous small islands in an area crossed by the powerful Kuroshio current. Oceanic currents can be effective seed-dispersal agents for coastal plant species. Moreover, the Luzon Strait is an area prone to tropical cyclones. The aim of this study is to look at the dispersal capability of an endangered coastal plant species, the Mearns fig (Ficus pedunculosa var. mearnsii), using both experimental and population genetics methods.

LOCATION

Southern Taiwan, the Philippines, and the islands between Luzon and Taiwan Island.

METHODS

This study combined two types of analysis, i.e., buoyancy experiments on syconia and double digest restriction-associated DNA sequencing (ddRAD), to analyze the population genetics of the Mearns fig.

RESULTS

We first discovered that mature Mearns fig syconia could float in seawater. They have a mean float duration of 10 days to a maximum of 21 days. Germination rates varied significantly between Mearns fig seeds that had undergone different durations of flotation treatment. Population genetic analysis shows a high degree of inbreeding among various Mearns fig populations. Moreover, no isolation by distance was found between the populations and individuals.

MAIN CONCLUSIONS

From our analysis of the genetic structure of the Mearns fig populations, we can clearly highlight the effect of the Kuroshio oceanic current on the seed dispersal of this fig tree. Comprehensive analysis has shown that Mearns fig seeds are still viable before the mature syconium sinks into the seawater, and so they could use the Kuroshio Current to float to the current population locations in Taiwan.

Diploid and tetraploid cytotypes of the flagship Cape species Dicerothamnus rhinocerotis (Asteraceae): variation in distribution, ecological niche, morphology and genetics

BACKGROUND AND AIMS

The Greater Cape Floristic Region is one of the world's biodiversity hotspots and is considered poor in polyploids. To test this assumption, ploidy variation was investigated in a widespread Cape shrub, Dicerothamnus rhinocerotis (renosterbos, Asteraceae). The aim was to elucidate the cytotype distribution and population composition across the species range, and to assess differences in morphology, environmental niches and genetics.

METHODS

Ploidy level and genome size were determined via flow cytometry and cytotype assignment was confirmed by chromosome counting. Restriction site-associated DNA sequencing (RADseq) analyses were used to infer genetic relationships. Cytotype climatic and environmental niches were compared using a range of environmental layers and a soil model, while morphological differences were examined using multivariate methods.

KEY RESULTS

The survey of 171 populations and 2370 individuals showed that the species comprises diploid and tetraploid cytotypes, no intermediates and only 16.8 % of mixed populations. Mean 2C values were 1.80-2.06 pg for diploids and 3.48-3.80 pg for tetraploids, with very similar monoploid genome sizes. Intra-cytotype variation showed a significant positive correlation with altitude and longitude in both cytotypes and with latitude in diploids. Although niches of both cytotypes were highly equivalent and similar, their optima and breadth were shifted due to differences mainly in isothermality and available water capacity. Morphometric analyses showed significant differences in the leaves and corolla traits, the number of florets per capitulum, and cypsela dimensions between the two cytotypes. Genetic analyses revealed four groups, three of them including both cytotypes.

CONCLUSIONS

Dicerothamnus rhinocerotis includes two distinct cytotypes that are genetically similar. While tetraploids arise several times independently within different genetic groups, morphological and ecological differences are evident between cytotypes. Our results open up new avenues for questions regarding the importance of ploidy in the megadiverse Cape flora, and exemplify the need for population-based studies focused on ploidy variation.

Unravelling some factors affecting sexual reproduction in rock-specialist shrub: Insight from an endemic Daphne arbuscula (Thymelaeaceae)

The role of endemic species in global biodiversity is pivotal, and understanding their biology and ecology is imperative for their fitness and long-term survival, particularly in the face of ongoing climatic oscillations. Our primary goal was to investigate the sexual reproduction level of the endangered Western Carpathian endemic Daphne arbuscula (Thymelaeaceae), which inhabits extreme rocky habitats, and to comprehend the influence of specific factors on its reproductive success. We conducted the research across four populations, varying in size and environmental conditions. Over two years, we monitored flower and fruit production, analyzed genetic variability within and among populations, and studied pollination mechanisms. Daphne arbuscula proved to be strictly self-incompatible, with significant variations in flower and fruit production among populations and seasons. The average fruit production percentage consistently remained below 50% across populations, indicating challenges in sexual reproduction. Cold and harsh weather during the reproductive phase had a substantial negative impact on sexual reproduction efficacy, leading to decreased fruit production. Nevertheless, several individuals in sheltered microhabitats displayed significantly higher fruit production, ranging from 60% to 83%, emphasizing the critical role of microhabitat heterogeneity in sustaining sexual reproduction in this species. We found no pronounced differences in genetic diversity within or among populations, suggesting that genetic factors may not critically influence the reproductive success of this endemic species. The implications of our findings might be of paramount importance for the long-term survival of D. arbuscula and offer valuable insights for the development of effective conservation strategies for this species.

Capturing and managing genetic diversity in ex situ collections of threatened tropical trees: A case study in Karomia gigas

Premise

Although ex situ collections of threatened plants are most useful when they contain maximal genetic variation, the conservation and maintenance of genetic diversity in collections are often poorly known. We present a case study using population genomic analyses of an ex situ collection of Karomia gigas, a critically endangered tropical tree from Tanzania. Only ~43 individuals are known in two wild populations, and ex situ collections containing 34 individuals were established in two sites from wild-collected seed. The study aimed to understand how much diversity is represented in the collection, analyze the parentage of ex situ individuals, and identify efficient strategies to capture and maintain genetic diversity.

Methods

We genotyped all known individuals using a 2b-RADseq approach, compared genetic diversity in wild populations and ex situ collections, and conducted parentage analysis of the collections.

Results

Wild populations were found to have greater levels of genetic diversity than ex situ populations as measured by number of private alleles, number of polymorphic sites, observed and expected heterozygosity, nucleotide diversity, and allelic richness. In addition, only 32.6% of wild individuals are represented ex situ and many individuals were found to be the product of selfing by a single wild individual.

Discussion

Population genomic analyses provided important insights into the conservation of genetic diversity in K. gigas, identifying gaps and inefficiencies, but also highlighting strategies to conserve genetic diversity ex situ. Genomic analyses provide essential information to ensure that collections effectively conserve genetic diversity in threatened tropical trees.

Effects of landscape, resource use, and body size on genetic structure in bee populations

Quantifying genetic structure and levels of genetic variation are fundamentally important to predicting the ability of populations to persist in human-altered landscapes and adapt to future environmental changes. Genetic structure reflects the dispersal of individuals over generations, which can be mediated by species-level traits or environmental factors. Dispersal distances are commonly positively associated with body size and negatively associated with the amount of degraded habitat between sites, motivating the investigation of these potential drivers of dispersal concomitantly. We quantified genetic structure and genetic variability within populations of seven bee species from the genus Euglossa across fragmented landscapes. We genotyped bees at SNP loci and tested the following predictions: (1) deforested areas restrict gene flow; (2) larger species have lower genetic structure; (3) species with greater resource specialization have higher genetic structure; and (4) sites surrounded by more intact habitat have higher genetic diversity. Contrasting with previous work on bees, we found no associations between body size and genetic structure. Genetic structure was higher for species with greater resource specialization, and the amount of intact habitat between or surrounding sites was positively associated with parameters reflecting gene flow and genetic diversity. These results challenge the dominant paradigm that individuals of larger species disperse farther, and they suggest that landscape and resource requirements are important factors mediating dispersal.

Estimates of heterozygosity from single nucleotide polymorphism markers are context-dependent and often wrong

Genetic diversity is frequently described using heterozygosity, particularly in a conservation context. Often, it is estimated using single nucleotide polymorphisms (SNPs); however, it has been shown that heterozygosity values calculated from SNPs can be biased by both study design and filtering parameters. Though solutions have been proposed to address these issues, our own work has found them to be inadequate in some circumstances. Here, we aimed to improve the reliability and comparability of heterozygosity estimates, specifically by investigating how sample size and missing data thresholds influenced the calculation of autosomal heterozygosity (heterozygosity calculated from across the genome, i.e. fixed and variable sites). We also explored how the standard practice of tri- and tetra-allelic site exclusion could bias heterozygosity estimates and influence eventual conclusions relating to genetic diversity. Across three distinct taxa (a frog, Litoria rubella; a tree, Eucalyptus microcarpa; and a grasshopper, Keyacris scurra), we found heterozygosity estimates to be meaningfully affected by sample size and missing data thresholds, partly due to the exclusion of tri- and tetra-allelic sites. These biases were inconsistent both between species and populations, with more diverse populations tending to have their estimates more severely affected, thus having potential to dramatically alter interpretations of genetic diversity. We propose a modified framework for calculating heterozygosity that reduces bias and improves the utility of heterozygosity as a measure of genetic diversity, whilst also highlighting the need for existing population genetic pipelines to be adjusted such that tri- and tetra-allelic sites be included in calculations.

Phylogenetic insights into Central European Chorthippus and Pseudochorthippus (Orthoptera: Acrididae) species using ddRADseq data

The evolution of several orthopteran groups, especially within the grasshopper family Acrididae, remains poorly understood. This is particularly true for the subfamily Gomphocerinae, which comprises cryptic sympatric and syntopic species. Previous mitochondrial studies have highlighted major discrepancies between taxonomic and phylogenetic hypotheses, thereby emphasizing the necessity of genome-wide approaches. In this study, we employ double-digest restriction site-associated DNA sequencing (ddRADseq) to reconstruct the evolution of Central European Chorthippus and Pseudochorthippus species, especially C.smardai, P.tatrae and the C.biguttulus group. Our phylogenomic analyses recovered deep discordance with mitochondrial DNA barcoding, emphasizing its unreliability in Gomphocerinae grasshoppers. Specifically, our data robustly distinguished the C.biguttulus group and confirmed the distinctiveness of C.eisentrauti, also shedding light on its presence in the Berchtesgaden Alps. Moreover, our results support the reclassification of C.smardai to the genus Pseudochorthippus and of P.tatrae to the genus Chorthippus. Our study demonstrates the efficiency of high-throughput genomic methods such as RADseq without prior optimization to elucidate the complex evolution of grasshopper radiations with direct taxonomic implications. While RADseq has predominantly been utilized for population genomics and within-genus phylogenomics, its application extends to resolve relationships between deeply-diverged clades representative of distinct genera.

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.

Genetic structure, UV-vision, wing coloration and size coincide with colour polymorphism in Fabriciana adippe butterflies

Colour polymorphisms have long served as model systems in evolutionary studies and continue to inform about processes involved in the origin and dynamics of biodiversity. Modern sequencing tools allow for evaluating whether phenotypic differences between morphs reflect genetic differentiation rather than developmental plasticity, and for investigating whether polymorphisms represent intermediate stages of diversification towards speciation. We investigated phenotypic and genetic differentiation between two colour morphs of the butterfly Fabriciana adippe using a combination of ddRAD-sequencing and comparisons of body size, colour patterns and optical properties of bright wing spots. The silvery-spotted adippe form had larger and darker wings and reflected UV light, while the yellow cleodoxa form displayed more green scales and reflected very little UV, showcasing that they constitute distinct and alternative integrated phenotypes. Genomic analyses revealed genetic structuring according to source population, and to colour morph, suggesting that the phenotypic differentiation reflects evolutionary modifications. We report 17 outlier loci associated with colour morph, including ultraviolet-sensitive visual pigment (UVRh1), which is associated with intraspecific communication and mate choice in butterflies. Together with the demonstration that the wings of the adippe (but essentially not the cleodoxa) morph reflect UV light, that UV reflectance is higher in females than males and that morphs differ in wing size, this suggests that these colour morphs might represent genetically integrated phenotypes, possibly adapted to different microhabitats. We propose that non-random mating might contribute to the differentiation and maintenance of the polymorphism.

Ex situ conservation of two rare oak species using microsatellite and SNP markers

Plant collections held by botanic gardens and arboreta are key components of ex situ conservation. Maintaining genetic diversity in such collections allows them to be used as resources for supplementing wild populations. However, most recommended minimum sample sizes for sufficient ex situ genetic diversity are based on microsatellite markers, and it remains unknown whether these sample sizes remain valid in light of more recently developed next-generation sequencing (NGS) approaches. To address this knowledge gap, we examine how ex situ conservation status and sampling recommendations differ when derived from microsatellites and single nucleotide polymorphisms (SNPs) in garden and wild samples of two threatened oak species. For Quercus acerifolia, SNPs show lower ex situ representation of wild allelic diversity and slightly lower minimum sample size estimates than microsatellites, while results for each marker are largely similar for Q. boyntonii. The application of missing data filters tends to lead to higher ex situ representation, while the impact of different SNP calling approaches is dependent on the species being analyzed. Measures of population differentiation within species are broadly similar between markers, but larger numbers of SNP loci allow for greater resolution of population structure and clearer assignment of ex situ individuals to wild source populations. Our results offer guidance for future ex situ conservation assessments utilizing SNP data, such as the application of missing data filters and the usage of a reference genome, and illustrate that both microsatellites and SNPs remain viable options for botanic gardens and arboreta seeking to ensure the genetic diversity of their collections.

Ultimate drivers of forced extra-pair copulations in birds lacking a penis: jackdaws as a case-study

Forced copulation is common, presumably because it can increase male reproductive success. Forced extra-pair copulation (FEPC) occurs in birds, even though most species lack a penis and are widely thought to require female cooperation for fertilization. How FEPC persists, despite a presumed lack of siring success and likely non-negligible costs to the male, is unknown. Using the jackdaw (Corvus monedula) as a case study, we use SNPs to quantify the extra-pair paternity rate through FEPC and evaluate explanations for the persistence of FEPC in species without a penis. We then collate evidence for FEPC across penis-lacking birds. Combining genetic and behavioural analyses, our study suggests that the most likely explanations for the maintenance of FEPC in jackdaws are that it provides a selective advantage to males or it is a relic. Our literature review shows that across birds lacking a penis, FEPC is taxonomically widespread, and yet, little is known about its evolution. A broader implementation of the approach used here, combining both genetic and behavioural data, may shed light on why this widespread sexual behaviour persists. Additional work is necessary to understand whether a penis is needed for paternity through forced copulation and to quantify the costs of FEPC.

Population genomic analyses reveal that salinity and geographic isolation drive diversification in a free-living protist

Protists make up the vast diversity of eukaryotic life and play a critical role in biogeochemical cycling and in food webs. Because of their small size, cryptic life cycles, and large population sizes, our understanding of speciation in these organisms is very limited. We performed population genomic analyses on 153 strains isolated from eight populations of the recently radiated dinoflagellate genus Apocalathium, to explore the drivers and mechanisms of speciation processes. Species of this genus inhabit both freshwater and saline habitats, lakes and seas, and are found in cold temperate environments across the world. RAD sequencing analyses revealed that the populations were overall highly differentiated, but morphological similarity was not congruent with genetic similarity. While geographic isolation was to some extent coupled to genetic distance, this pattern was not consistent. Instead, we found evidence that the environment, specifically salinity, is a major factor in driving ecological speciation in Apocalathium. While saline populations were unique in loci coupled to genes involved in osmoregulation, freshwater populations appear to lack these. Our study highlights that adaptation to freshwater through loss of osmoregulatory genes may be an important speciation mechanism in free-living aquatic protists.

De Novo Genome Assembly of the Whitespot Parrotfish (Scarus forsteni): A Valuable Scaridae Genomic Resource

Scarus forsteni, a whitespot parrotfish from the Scaridae family, is a herbivorous fish inhabiting coral reef ecosystems. The deterioration of coral reefs has highly affected the habitats of the parrotfish. The decline in genetic diversity of parrotfish emphasizes the critical importance of conserving their genetic variability to ensure the resilience and sustainability of marine ecosystems for future generations. In this study, a genome of S. forsteni was assembled de novo through using Illumina and Nanopore sequencing. The 1.71-Gb genome of S. forsteni, was assembled into 544 contigs (assembly level: contig). It exhibited an N50 length of 17.97 Mb and a GC content percentage of 39.32%. Our BUSCO analysis revealed that the complete protein of the S. forsteni genome had 98.10% integrity. Combined with structure annotation data, 34,140 (74.81%) genes were functionally annotated out of 45,638 predicted protein-coding genes. Upon comparing the genome size and TE content of teleost fishes, a roughly linear relationship was observed between these two parameters. However, TE content is not a decisive factor in determining the genome size of S. forsteni. Population history analysis results indicate that S. forsteni experienced two major population expansions, both of which occurred before the last interglacial period. In addition, through a comparative genomic analysis of the evolutionary relationship of other species, it was found that S. forsteni had the closest relationship with Cheilinus undulatus, another member of the Labridae family. Our expansion and contraction analysis of the gene family showed that the expansion genes were mainly associated with immune diseases, organismal systems, and cellular processes. At the same time, cell transcription and translation, sex hormone regulation, and other related pathways were also more prominent in the positive selection genes. The genomic sequence of S. forsteni offers valuable resources for future investigations on the conservation, evolution, and behavior of fish species.

Phylogeography and phenotypic wing shape variation in a damselfly across populations in Europe

BACKGROUND

Describing geographical variation in morphology of organisms in combination with data on genetic differentiation and biogeography can provide important information on how natural selection shapes such variation. Here we study genetic structure using ddRAD seq and wing shape variation using geometric morphometrics in 14 populations of the damselfly Lestes sponsa along its latitudinal range in Europe.

RESULTS

The genetic analysis showed a significant, yet relatively weak population structure with high genetic heterozygosity and low inbreeding coefficients, indicating that neutral processes contributed very little to the observed wing shape differences. The genetic analysis also showed that some regions of the genome (about 10%) are putatively shaped by selection. The phylogenetic analysis showed that the Spanish and French populations were the ancestral ones with northern Swedish and Finnish populations being the most derived ones. We found that wing shape differed significantly among populations and showed a significant quadratic (but weak) relationship with latitude. This latitudinal relationship was largely attributed to allometric effects of wing size, but non-allometric variation also explained a portion of this relationship. However, wing shape showed no phylogenetic signal suggesting that lineage-specific variation did not contribute to the variation along the latitudinal gradient. In contrast, wing size, which is correlated with body size in L. sponsa, had a strong negative correlation with latitude.

CONCLUSION

Our results suggest a relatively weak population structure among the sampled populations across Europe, but a clear differentiation between south and north populations. The observed geographic phenotypic variation in wing shape may have been affected by different local selection pressures or environmental effects.

Pleistocene diversification of unifoliolate-leaved Lupinus (Leguminosae: Papilionoideae) in Florida

The importance and prevalence of recent ice-age and post-glacial speciation and species diversification during the Pleistocene across many organismal groups and physiographic settings are well established. However, the extent to which Pleistocene diversification can be attributed to climatic oscillations and their effects on distribution ranges and population structure remains debatable. In this study, we use morphologic, geographic and genetic (RADseq) data to document Pleistocene speciation and intra-specific diversification of the unifoliolate-leaved clade of Florida Lupinus, a small group of species largely restricted to inland and coastal sand ridges across the Florida peninsula and panhandle. Phylogenetic and demographic analyses alongside morphological and geographic evidence suggest that recent speciation and intra-specific divergence within this clade were driven by a combination of non-adaptive allopatric divergence caused by edaphic niche conservatism and opportunities presented by the emergence of new post-glacial sand ridge habitats. These results highlight the central importance of even modest geographic isolation and short periods of allopatric divergence following range expansion in the emergence of new taxa and add to the growing evidence that Pleistocene climatic oscillations may contribute to rapid diversification in a myriad of physiographic settings. Furthermore, our results shed new light on long-standing taxonomic debate surrounding the number of species in the Florida unifoliate Lupinus clade providing support for recognition of five species and a set of intra-specific variants. The important conservation implications for the narrowly restricted, highly endangered species Lupinus aridorum, which we show to be genetically distinct from its sister species Lupinus westianus, are discussed.

A lack of genetic diversity and minimal adaptive evolutionary divergence in introduced Mysis shrimp after 50 years

The successes of introduced populations in novel habitats often provide powerful examples of evolution and adaptation. In the 1950s, opossum shrimp (Mysis diluviana) individuals from Clearwater Lake in Minnesota, USA were transported and introduced to Twin Lakes in Colorado, USA by fisheries managers to supplement food sources for trout. Mysis were subsequently introduced from Twin Lakes into numerous lakes throughout Colorado. Because managers kept detailed records of the timing of the introductions, we had the opportunity to test for evolutionary divergence within a known time interval. Here, we used reduced representation genomic data to investigate patterns of genetic diversity, test for genetic divergence between populations, and for evidence of adaptive evolution within the introduced populations in Colorado. We found very low levels of genetic diversity across all populations, with evidence for some genetic divergence between the Minnesota source population and the introduced populations in Colorado. There was little differentiation among the Colorado populations, consistent with the known provenance of a single founding population, with the exception of the population from Gross Reservoir, Colorado. Demographic modeling suggests that at least one undocumented introduction from an unknown source population hybridized with the population in Gross Reservoir. Despite the overall low genetic diversity we observed, F ST outlier and environmental association analyses identified multiple loci exhibiting signatures of selection and adaptive variation related to elevation and lake depth. The success of introduced species is thought to be limited by genetic variation, but our results imply that populations with limited genetic variation can become established in a wide range of novel environments. From an applied perspective, the observed patterns of divergence between populations suggest that genetic analysis can be a useful forensic tool to determine likely sources of invasive species.

Cytotaxonomic characterization and estimation of migration patterns of onchocerciasis vectors (Simulium damnosum sensu lato) in northwestern Ethiopia based on RADSeq data

BACKGROUND

While much progress has been made in the control and elimination of onchocerciasis across Africa, the extent to which vector migration might confound progress towards elimination or result in re-establishment of endemism in areas where transmission has been eliminated remains unclear. In Northern Ethiopia, Metema and Metekel-two foci located near the Sudan border-exhibit continuing transmission. While progress towards elimination has been faster in Metema, there remains a problematic hotspot of transmission. Whether migration from Metekel contributes to this is currently unknown.

METHODOLOGY/PRINCIPLE FINDINGS

To assess the role of vector migration from Metekel into Metema, we present a population genomics study of 151 adult female vectors using 47,638 RADseq markers and mtDNA CoI sequencing. From additional cytotaxonomy data we identified a new cytoform in Metema, closely related to S. damnosum s.str, here called the Gondar form. RADseq data strongly indicate the existence of two distinctly differentiated clusters within S. damnosum s.l.: one genotypic cluster found only in Metema, and the second found predominantly in Metekel. Because blackflies from both clusters were found in sympatry (in all four collection sites in Metema), but hybrid genotypes were not detected, there may be reproductive barriers preventing interbreeding. The dominant genotype in Metema was not found in Metekel while the dominant genotype in Metekel was found in Metema, indicating that (at the time of sampling) migration is primarily unidirectional, with flies moving from Metekel to Metema. There was strong differentiation between clusters but little genetic differentiation within clusters, suggesting migration and gene flow of flies within the same genetic cluster are sufficient to prevent genetic divergence between sites.

CONCLUSIONS/SIGNIFICANCE

Our results confirm that Metekel and Metema represent different transmission foci, but also indicate a northward movement of vectors between foci that may have epidemiological importance, although its significance requires further study.

Genomic differentiation and niche divergence in the Hetaerina americana (Odonata) cryptic species complex

The evolution of reproductive barriers, that is, the speciation process, implies the limitation of gene flow between populations. Different patterns of genomic differentiation throughout the speciation continuum may provide insights into the causal evolutionary forces of species divergence. In this study, we analysed a cryptic species complex of the genus Hetaerina (Odonata). This complex includes H. americana and H. calverti; however, in H. americana two highly differentiated genetic groups have been previously detected, which, we hypothesize, may correspond to different species with low morphological variation. We obtained single nucleotide polymorphism (SNP) data for 90 individuals belonging to the different taxa in the complex and carried out differentiation tests to identify genetic isolation. The results from STRUCTURE and discriminant analysis of principal components (DAPC), based on almost 5000 SNPs, confirmed the presence of three highly differentiated taxa. Also, we found FST values above 0.5 in pairwise comparisons, which indicates a considerable degree of genetic isolation among the suggested species. We also found low climatic niche overlap among all taxa, suggesting that each group occurs at specific conditions of temperature, precipitation and elevation. We propose that H. americana comprises two cryptic species, which may be reproductively isolated by ecological barriers related to niche divergence, since the morphological variation is minimal and, therefore, mechanical barriers are probably less effective compared to other related species such as H. calverti.

Species-specific effects of production practices on genetic diversity in plant reintroduction programs

Plant production practices can influence the genetic diversity of cultivated plant materials and, ultimately, their potential to adapt to a reintroduction site. A common step in the plant production process is the application of seed pretreatment to alleviate physiological seed dormancy and successfully germinate seeds. In production settings, the seeds that germinate more rapidly may be favored in order to fill plant quotas. In this study, we investigated how the application of cold-moist stratification treatments with different durations can lead to differences in the genetic diversity of the propagated plant materials. Specifically, we exposed seeds of three Viola species to two different cold stratification durations, and then we analyzed the genetic diversity of the resulting subpopulations through double-digestion restriction site-associated sequencing (ddRADseq). Our results show that, in two out of three species, utilizing a short stratification period will decrease the genetic diversity of neutral and expressed loci, likely due to the imposition of a genetic bottleneck and artificial selection. We conclude that, in some species, the use of minimal stratification practices in production may jeopardize the adaptive potential and long-term persistence of reintroduced populations and suggest that practitioners carefully consider the evolutionary implications of their production protocols. We highlight the need to consider the germination ecology of target species when selecting the length of dormancy-breaking pretreatments.

Additive and non-additive genetic variance in juvenile Sitka spruce (Picea sitchensis Bong. Carr)

Many quantitative genetic models assume that all genetic variation is additive because of a lack of data with sufficient structure and quality to determine the relative contribution of additive and non-additive variation. Here the fractions of additive (fa) and non-additive (fd) genetic variation were estimated in Sitka spruce for height, bud burst and pilodyn penetration depth. Approximately 1500 offspring were produced in each of three sib families and clonally replicated across three geographically diverse sites. Genotypes from 1525 offspring from all three families were obtained by RADseq, followed by imputation using 1630 loci segregating in all families and mapped using the newly developed linkage map of Sitka spruce. The analyses employed a new approach for estimating fa and fd, which combined all available genotypic and phenotypic data with spatial modelling for each trait and site. The consensus estimate for fa increased with age for height from 0.58 at 2 years to 0.75 at 11 years, with only small overlap in 95% support intervals (I95). The estimated fa for bud burst was 0.83 (I95=[0.78, 0.90]) and 0.84 (I95=[0.77, 0.92]) for pilodyn depth. Overall, there was no evidence of family heterogeneity for height or bud burst, or site heterogeneity for pilodyn depth, and no evidence of inbreeding depression associated with genomic homozygosity, expected if dominance variance was the major component of non-additive variance. The results offer no support for the development of sublines for crossing within the species. The models give new opportunities to assess more accurately the scale of non-additive variation.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11295-023-01627-5.

Genetic depletion does not prevent rapid evolution in island-introduced lizards

Experimental introductions of species have provided some of the most tractable examples of rapid phenotypic changes, which may reflect plasticity, the impact of stochastic processes, or the action of natural selection. Yet to date, very few studies have investigated the neutral and potentially adaptive genetic impacts of experimental introductions. We dissect the role of these processes in shaping the population differentiation of wall lizards in three Croatian islands (Sušac, Pod Kopište, and Pod Mrčaru), including the islet of Pod Mrčaru, where experimentally introduced lizards underwent rapid (~30 generations) phenotypic changes associated with a shift from an insectivorous to a plant-based diet. Using a genomic approach (~82,000 ddRAD loci), we confirmed a founder effect during introduction and very low neutral genetic differentiation between the introduced population and its source. However, genetic depletion did not prevent rapid population growth, as the introduced lizards exhibited population genetic signals of expansion and are known to have reached a high density. Our genome-scan analysis identified just a handful of loci showing large allelic shifts between ecologically divergent populations. This low overall signal of selection suggests that the extreme phenotypic differences observed among populations are determined by a small number of large-effect loci and/or that phenotypic plasticity plays a major role in phenotypic changes. Nonetheless, functional annotation of the outlier loci revealed some candidate genes relevant to diet-induced adaptation, in agreement with the hypothesis of directional selection. Our study provides important insights on the evolutionary potential of bottlenecked populations in response to new selective pressures on short ecological timescales.

First Insights into Population Structure and Genetic Diversity Versus Host Specificity in Trypanorhynch Tapeworms Using Multiplexed Shotgun Genotyping

Theory predicts relaxed host specificity and high host vagility should contribute to reduced genetic structure in parasites while strict host specificity and low host vagility should increase genetic structure. Though these predictions are intuitive, they have never been explicitly tested in a population genomic framework. Trypanorhynch tapeworms, which parasitize sharks and rays (elasmobranchs) as definitive hosts, are the only order of elasmobranch tapeworms that exhibit considerable variability in their definitive host specificity. This allows for unique combinations of host use and geographic range, making trypanorhynchs ideal candidates for studying how these traits influence population-level structure and genetic diversity. Multiplexed shotgun genotyping (MSG) data sets were generated to characterize component population structure and infrapopulation diversity for a representative of each trypanorhynch suborder: the ray-hosted Rhinoptericola megacantha (Trypanobatoida) and the shark-hosted Callitetrarhynchus gracilis (Trypanoselachoida). Adults of R. megacantha are more host-specific and less broadly distributed than adults of C. gracilis, allowing correlation between these factors and genetic structure. Replicate tapeworm specimens were sequenced from the same host individual, from multiple conspecific hosts within and across geographic regions, and from multiple definitive host species. For R. megacantha, population structure coincided with geography rather than host species. For C. gracilis, limited population structure was found, suggesting a potential link between degree of host specificity and structure. Conspecific trypanorhynchs from the same host individual were found to be as, or more, genetically divergent from one another as from conspecifics from different host individuals. For both species, high levels of homozygosity and positive FIS values were documented.

Contrasting whole-genome and reduced representation sequencing for population demographic and adaptive inference: an alpine mammal case study

Genomes capture the adaptive and demographic history of a species, but the choice of sequencing strategy and sample size can impact such inferences. We compared whole genome and reduced representation sequencing approaches to study the population demographic and adaptive signals of the North American mountain goat (Oreamnos americanus). We applied the restriction site-associated DNA sequencing (RADseq) approach to 254 individuals and whole genome resequencing (WGS) approach to 35 individuals across the species range at mid-level coverage (9X) and to 5 individuals at high coverage (30X). We used ANGSD to estimate the genotype likelihoods and estimated the effective population size (Ne), population structure, and explicitly modelled the demographic history with δaδi and MSMC2. The data sets were overall concordant in supporting a glacial induced vicariance and extremely low Ne in mountain goats. We evaluated a set of climatic variables and geographic location as predictors of genetic diversity using redundancy analysis. A moderate proportion of total variance (36% for WGS and 21% for RADseq data sets) was explained by geography and climate variables; both data sets support a large impact of drift and some degree of local adaptation. The empirical similarities of WGS and RADseq presented herein reassuringly suggest that both approaches will recover large demographic and adaptive signals in a population; however, WGS offers several advantages over RADseq, such as inferring adaptive processes and calculating runs-of-homozygosity estimates. Considering the predicted climate-induced changes in alpine environments and the genetically depauperate mountain goat, the long-term adaptive capabilities of this enigmatic species are questionable.

The first linkage map for Australo-Papuan Treefrogs (family: Pelodryadidae) reveals the sex-determination system of the Green-eyed Treefrog (Litoria serrata)

Amphibians represent a useful taxon to study the evolution of sex determination because of their highly variable sex-determination systems. However, the sex-determination system for many amphibian families remains unknown, in part because of a lack of genomic resources. Here, using an F1 family of Green-eyed Treefrogs (Litoria serrata), we produce the first genetic linkage map for any Australo-Papuan Treefrogs (family: Pelodryadidae). The resulting linkage map contains 8662 SNPs across 13 linkage groups. Using an independent set of sexed adults, we identify a small region in linkage group 6 matching an XY sex-determination system. These results suggest Litoria serrata possesses a male heterogametic system, with a candidate sex-determination locus on linkage group 6. Furthermore, this linkage map represents the first genomic resource for Australo-Papuan Treefrogs, an ecologically diverse family of over 220 species.

No evidence of inbreeding depression despite a historical severe bottleneck in the endangered Bermuda petrel (Pterodroma cahow)

The Bermuda petrel Pterodroma cahow is an island endemic seabird that belongs to the Procellariiformes, one of the most endangered orders of birds. Historical records suggest a significant population size decline following human settlement in Bermuda, bringing the species to near extinction. Since the 1950s, the population has been recovering aided by the implementation of an ongoing conservation plan. However, it still faces several threats, and negative genetic effects resulting from that drastic decline are to be expected, including inbreeding and genetic drift. We studied genetic diversity and levels of inbreeding, and their effects on individual fitness and mating choice. We also tested for a genetic signature of the recent demographic bottleneck. For this, we analyzed variation in thousands of nuclear single-nucleotide polymorphisms derived from double digest restriction site-associated DNA sequencing and 1 mitochondrial gene (cytochrome oxidase I). The results revealed that the Bermuda petrel suffered a recent genetic bottleneck and shows low mitochondrial diversity compared with other petrel species. Conversely, nuclear diversity was similar to that of other endangered petrels. Inbreeding levels were not high overall, although some individuals were highly inbred. However, we found no evidence that individual inbreeding or relatedness between mates affected hatching success, or that mate choice is influenced by kinship in this very small population.

History of the terrestrial isopod genus Ligidium in Japan based on phylogeographic analysis

BACKGROUND

Phylogeographical approaches explain the genetic diversity of local organisms in the context of their geological and geographic environments. Thus, genetic diversity can be a proxy for geological history. Here we propose a genus of woodland isopod, Ligidium, as a marker of geological history in relation to orogeny and the Quaternary glacial cycle.

RESULTS

Mitochondrial analysis of 721 individuals from 97 sites across Japan revealed phylogenetic divergence between the northeastern and southwestern Japan arcs. It also showed repeated population expansions in northeastern Japan in response to Quaternary glacial and interglacial cycles. Genome-wide analysis of 83 selected individuals revealed multiple genetic nuclear clusters. The genomic groupings were consistent with the local geographic distribution, indicating that the Ligidium phylogeny reflects its regional history.

CONCLUSION

Ligidium DNA sequence analysis can provide insight into the geological, geographical, and paleoenvironmental history of the studied region.

A Comparative Phylogeography of Three Marine Species with Different PLD Modes Reveals Two Genetic Breaks across the Southern Caribbean Sea

The comparative phylogeography of marine species with contrasting dispersal potential across the southern Caribbean Sea was evaluated by the presence of two putative barriers: the Magdalena River plume (MRP) and the combination of the absence of a rocky bottom and the almost permanent upwelling in the La Guajira Peninsula (ARB + PUG). Three species with varying biological and ecological characteristics (i.e., dispersal potentials) that inhabit shallow rocky bottoms were selected: Cittarium pica (PLD < 6 days), Acanthemblemaria rivasi (PLD < 22 days), and Nerita tessellata (PLD > 60 days). We generated a set of SNPs for the three species using the ddRad-seq technique. Samples of each species were collected in five locations from Capurganá to La Guajira. For the first time, evidence of a phylogeographic break caused by the MRP is provided, mainly for A. rivasi (AMOVA: ΦCT = 0.420). The ARB + PUG barrier causes another break for A. rivasi (ΦCT = 0.406) and C. pica (ΦCT = 0.224). Three populations (K = 3) were identified for A. rivasi and C. pica, while N. tessellata presented one population (K = 1). The Mantel correlogram indicated that A. rivasi and C. pica fit the hierarchical population model, and only the A. rivasi and C. pica comparisons showed phylogeographic congruence. Our results demonstrate how the biological traits of these three species and the biogeographic barriers have influenced their phylogeographic structure.

Exceptional population genomic homogeneity in the black brittle star Ophiocomina nigra (Ophiuroidea, Echinodermata) along the Atlantic-Mediterranean coast

The Atlantic-Mediterranean marine transition is characterised by strong oceanographic barriers and steep environmental gradients that generally result in connectivity breaks between populations from both basins and may lead to local adaptation. Here, we performed a population genomic study of the black brittle star, Ophiocomina nigra, covering most of its distribution range along the Atlantic-Mediterranean region. Interestingly, O. nigra is extremely variable in its coloration, with individuals ranging from black to yellow-orange, and different colour morphs inhabiting different depths and habitats. In this work, we used a fragment of the mitochondrial COI gene and 2,374 genome-wide ddRADseq-derived SNPs to explore: (a) whether the different colour morphs of O. nigra represent different evolutionary units; (b) the disruptive effects of major oceanographic fronts on its population structure; and (c) genomic signals of local adaptation to divergent environments. Our results revealed exceptional population homogeneity, barely affected by oceanographic fronts, with no signals of local adaptation nor genetic differentiation between colour morphs. This remarkable panmixia likely results from a long pelagic larval duration, a large effective population size and recent demographic expansions. Our study unveils an extraordinary phenotypic plasticity in O. nigra, opening further research questions on the ecological and molecular mechanisms underpinning coloration in Ophiuroidea.

Mapping the root systems of individual trees in a natural community using genotyping-by-sequencing

The architecture of root systems is an important driver of plant fitness, competition and ecosystem processes. However, the methodological difficulty of mapping roots hampers the study of these processes. Existing approaches to match individual plants to belowground samples are low throughput and species specific. Here, we developed a scalable sequencing-based method to map the root systems of individual trees across multiple species. We successfully applied it to a tropical dry forest community in the Brazilian Caatinga containing 14 species. We sequenced all 42 individual shrubs and trees in a 14 × 14 m plot using double-digest restriction site-associated sequencing (ddRADseq). We identified species-specific markers and individual-specific haplotypes from the data. We matched these markers to the ddRADseq data from 100 mixed root samples from across the centre (10 × 10 m) of the plot at four different depths using a newly developed R package. We identified individual root samples for all species and all but one individual. There was a strong significant correlation between belowground and aboveground size measurements, and we also detected significant species-level root-depth preference for two species. The method is more scalable and less labour intensive than the current techniques and is broadly applicable to ecology, forestry and agricultural biology.

Speciation of rock-dwelling snail species: disjunct ranges and mosaic patterns reveal the importance of long-distance dispersal in Chilostoma (Cingulifera) in the European Southern Alps

To better understand the origin of the high diversity and endemism in the Southern Alps of Europe, we investigated the phylogeny and population structure of the rock-dwelling snail group Chilostoma (Cingulifera) in the Southern Alps. We generated genomic ddRAD data and mitochondrial sequences of 104 Cingulifera specimens from 28 populations and 14 other Ariantinae. Until recently, about 30 Cingulifera taxa were classified as subspecies of a single polytypic species. The phylogenetic and population genetic analyses of the ddRAD data and mitochondrial sequences revealed that Cingulifera in the Southern Alps is differentiated into three species. Each of the three Chilostoma (Cingulifera) species occupies disjunct sub-areas, which are separated by areas occupied by other Chilostoma taxa. Neighbouring populations of different species show little or no admixture. Tests indicating that the genetic differentiation of the three Cingulifera taxa cannot be explained by isolation by distance confirmed their species status. The disjunct range patterns demonstrate the importance of stochastic events such as passive long-distance dispersal for the evolution of population structure and speciation in these snails, and of priority effects and ecological competition as important factors influencing species distributions.

Inter- and intra-island speciation and their morphological and ecological correlates in Aeonium (Crassulaceae), a species-rich Macaronesian radiation

BACKGROUND AND AIMS

The most species-rich and ecologically diverse plant radiation on the Canary Islands is the Aeonium alliance (Crassulaceae). In island radiations like this, speciation can take place either within islands or following dispersal between islands. Aiming at quantifying intra- and inter-island speciation events in the evolution of Aeonium, and exploring their consequences, we hypothesized that (1) intra-island diversification resulted in stronger ecological divergence of sister lineages, and that (2) taxa on islands with a longer history of habitation by Aeonium show stronger ecological differentiation and produce fewer natural hybrids.

METHODS

We studied the biogeographical and ecological setting of diversification processes in Aeonium with a fully sampled and dated phylogeny inferred using a ddRADseq approach. Ancestral areas and biogeographical events were reconstructed in BioGeoBEARS. Eleven morphological characters and three habitat characteristics were taken into account to quantify the morphological and ecological divergence between sister lineages. A co-occurrence matrix of all Aeonium taxa is presented to assess the spatial separation of taxa on each island.

KEY RESULTS

We found intra- and inter-island diversification events in almost equal numbers. In lineages that diversified within single islands, morphological and ecological divergence was more pronounced than in lineages derived from inter-island diversification, but only the difference in morphological divergence was significant. Those islands with the longest history of habitation by Aeonium had the lowest percentages of co-occurring and hybridizing taxon pairs compared with islands where Aeonium arrived later.

CONCLUSIONS

Our findings illustrate the importance of both inter- and intra-island speciation, the latter of which is potentially sympatric speciation. Speciation on the same island entailed significantly higher levels of morphological divergence compared with inter-island speciation, but ecological divergence was not significantly different. Longer periods of shared island habitation resulted in the evolution of a higher degree of spatial separation and stronger reproductive barriers.

Using genomics, morphometrics, and environmental niche modeling to test the validity of a narrow-range endemic snail, Pateranantahala (Gastropoda, Polygyridae)

Terrestrial gastropods are among the most imperiled groups of organisms on Earth. Many species have a complex taxonomic history, often including poorly defined subspecies, most of which have not been the focus of modern systematics research. Genomic tools, geometric morphometrics, and environmental niche modeling were used to assess the taxonomic status of Pateraclarkiinantahala (Clench & Banks, 1932), a subspecies of high conservation concern with a restricted range of approximately 3.3 km² in North Carolina, USA. A genome-scale dataset was generated that included individuals with morphologies matching P.c.nantahala, P.c.clarkii, and one individual with an intermediate form between P.c.nantahala and P.c.clarkii that was initially hypothesized as a potential hybrid. Mitochondrial phylogenetics, nuclear species tree inference, and phylogenetic networks were used to assess relationships and gene flow. Differences in shell shape via geometric morphometrics and whether the environmental niches of the two subspecies were significantly different were also examined. Molecular analyses indicated an absence of gene flow among lineages of P.clarkii sensu lato. Analyses rejected our hypothesis that the intermediate shelled form represented a hybrid, but instead indicated that it was a distinct lineage. Environmental niche models indicated significant differences in environmental niche between P.c.clarkii and P.c.nantahala, and geometric morphometrics indicated that P.c.nantahala had a significantly different shell shape. Given multiple lines of evidence, species-level recognition of P.nantahala is warranted.

Genome-environment associations along elevation gradients in two snowbed species of the North-Eastern Calcareous Alps

BACKGROUND

Anthropogenic climate change leads to increasing temperatures and altered precipitation and snowmelt patterns, especially in alpine ecosystems. To understand species' responses to climate change, assessment of genetic structure and diversity is crucial as the basis for the evaluation of migration patterns, genetic adaptation potential as well as the identification of adaptive alleles.

RESULTS

We studied genetic structure, diversity and genome-environment associations of two snowbed species endemic to the Eastern Alps with a large elevational range, Achillea clusiana Tausch and Campanula pulla L. Genotyping-by-sequencing was employed to assemble loci de novo, call variants and perform population genetic analyses. Populations of either species were distinguishable by mountain, and to some extent by elevation. We found evidence for gene flow between elevations. Results of genome-environment associations suggested similar selective pressures acting on both species, emanating mainly from precipitation and exposition rather than temperature.

CONCLUSIONS

Given their genetic structure and amount of gene flow among populations the two study species are suitable to serve as a model for genetic monitoring of climate change adaptation along an elevation gradient. Consequences of climate change will predominantly manifest via changes in precipitation and, thus, duration of snow cover in the snowbeds and indirectly via shrub encroachment accompanied by increasing shading of snowbeds at lower range margins. Assembling genomes of the study species and studying larger sample sizes and time series will be necessary to functionally characterize and validate the herein identified genomic loci putatively involved in adaptive processes.

Genomic phylogeography illuminates deep cyto-nuclear discordances in midwife toads (Alytes)

The advent of genomic methods allows us to revisit the evolutionary history of organismal groups for which robust phylogenies are still lacking, particularly in species complexes that frequently hybridize. In this study, we conduct RAD-sequencing (RAD-seq) analyses of midwife toads (genus Alytes), an iconic group of western Mediterranean amphibians famous for their parental care behavior, but equally infamous for the difficulties to reconstruct their evolutionary history. Through admixture and phylogenetic analyses of thousands of loci, we provide the most comprehensive phylogeographic framework for the A. obstetricans complex to date, as well as the first fully resolved phylogeny for the entire genus. As part of this effort, we carefully explore the influence of different sampling schemes and data filtering thresholds on tree reconstruction, showing that several, slightly different, yet robust topologies may be retrieved with small datasets obtained by stringent SNP calling parameters, especially when admixed individuals are included. In contrast, analyses of incomplete but larger datasets converged on the same phylogeny, irrespective of the reconstruction method used or the proportion of missing data. The Alytes tree features three Miocene-diverged clades corresponding to the proposed subgenera Ammoryctis (A. cisternasii), Baleaphryne (A. maurus, A. dickhilleni and A. muletensis), and Alytes (A. obstetricans complex). The latter consists of six evolutionary lineages, grouped into three clades of Pliocene origin, and currently delimited as two species: (1) A. almogavarii almogavarii and A. a. inigoi; (2) A. obstetricans obstetricans and A. o. pertinax; (3) A. o. boscai and an undescribed taxon (A. o. cf. boscai). These results contradict the mitochondrial tree, due to past mitochondrial captures in A. a. almogavarii (central Pyrenees) and A. o. boscai (central Iberia) by A. obstetricans ancestors during the Pleistocene. Patterns of admixture between subspecies appear far more extensive than previously assumed from microsatellites, causing nomenclatural uncertainties, and even underlying the reticulate evolution of one taxon (A. o. pertinax). All Ammoryctis and Baleaphryne species form shallow clades, so their taxonomy should remain stable. Amid the prevalence of cyto-nuclear discordance among terrestrial vertebrates and the usual lack of resolution of conventional nuclear markers, our study advocates for phylogeography based on next-generation sequencing, but also encourages properly exploring parameter space and sampling schemes when building and analyzing genomic datasets.

The Origin and Invasion Pathway of Brown Rats Rattus norvegicus on Dok-Do Island Revealed by Genome-Wide Markers from 3-RADseq Approach

Biological invasions are known to cause local extinctions on islands. Dok-do, a small, remote volcanic island in the East Sea of Korea in the western Pacific, has recently been invaded by rats, posing ecological problems. To infer their origin and invasion pathway, we collected rats from Dok-do and from the potential introduction source locations, Ulleung-do in the Pacific Ocean, and four east coastal ports. First, we identified that the brown rat (Rattus norvegicus) was the only rat species occurring at collecting sites based on the key morphological characteristics. To determine the population-level genetic diversity pattern, we applied the 3-RADseq approach. After a series of filtrations (minor allele frequency < 0.05, Hardy–Weinberg equilibrium p < 1 × 10−7), 4042 SNPs were retained for the final dataset from the 25,439 SNPs initially isolated. The spatial structure and genetic diversity pattern of brown rats suggested that the rat population on Dok-do was likely introduced from Ulleung-do. Our work provides practical information that will assist in the management of invasive brown rats in vulnerable island ecosystems.

Morphology, nuclear SNPs and mate selection reveal that COI barcoding overestimates species diversity in a Mediterranean freshwater amphipod by an order of magnitude

DNA sequence information has revealed many morphologically cryptic species worldwide. For animals, DNA-based assessments of species diversity usually rely on the mitochondrial cytochrome c oxidase subunit I (COI) gene. However, a growing amount of evidence indicate that mitochondrial markers alone can lead to misleading species diversity estimates due to mito-nuclear discordance. Therefore, reports of putative species based solely on mitochondrial DNA should be verified by other methods, especially in cases where COI sequences are identical for different morphospecies or where divergence within the same morphospecies is high. Freshwater amphipods are particularly interesting in this context because numerous putative cryptic species have been reported. Here, we investigated the species status of the numerous mitochondrial molecular operational taxonomic units (MOTUs) found within Echinogammarus sicilianus. We used an integrative approach combining DNA barcoding with mate selection observations, detailed morphometrics and genome-wide double digest restriction site-associated DNA sequencing (ddRAD-seq). Within a relatively small sampling area, we detected twelve COI MOTUs (divergence = 1.8-20.3%), co-occurring in syntopy at two-thirds of the investigated sites. We found that pair formation was random and there was extensive nuclear gene flow among the ten MOTUs co-occurring within the same river stretch. The four most common MOTUs were also indistinguishable with respect to functional morphology. Therefore, the evidence best fits the hypothesis of a single, yet genetically diverse, species within the main river system. The only two MOTUs sampled outside the focal area were genetically distinct at the nuclear level and may represent distinct species. Our study reveals that COI-based species delimitation can significantly overestimate species diversity, highlighting the importance of integrative taxonomy for species validation, especially in hyperdiverse complexes with syntopically occurring mitochondrial MOTUs.

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.

Genetic relationships among laboratory lines of the egg parasitoid Trissolcus japonicus from native and adventive populations

Candidate biological control agents of invasive insect pests are increasingly being found in new geographic regions as a result of unintentional introductions. However, testing the degree of genetic differentiation among adventive and native-range populations of these agents is rarely done. We used reduced-representation sequencing of genomic DNA to investigate the relationships among laboratory lines of Trissolcus japonicus (Ashmead) (Hymenoptera, Scelionidae), an egg parasitoid and biological control agent of the brown marmorated stink bug, Halyomorpha halys (Stål) (Hemiptera, Pentatomidae). We compared sequences from multiple adventive populations in North America (Canada, USA) and Europe (Switzerland) with populations sourced from part of its native range in China. We found considerably more genetic variation among lines sourced from adventive populations than among those within native populations. In the Pacific Northwest of North America (British Columbia, Canada and Washington State, USA), we found preliminary evidence of three distinct genetic clusters, two of which were highly dissimilar from all other lines we genotyped. In contrast, we found that other adventive lines with close geographic proximity (two from Ontario, Canada, three from Switzerland) had limited genetic variation. These findings provide a basis for testing biological differences among lines that will inform their use as biological control agents, and provide evidence to support a hypothesis of several independent introductions of T. japonicus in western North America from different source areas.

Large effective population size masks population genetic structure in Hirondellea amphipods within the deepest marine ecosystem, the Mariana Trench

The examination of genetic structure in the deep-ocean hadal zone has focused on divergence between tectonic trenches to understand how environment and geography may drive species divergence and promote endemism. There has been little attempt to examine localized genetic structure within trenches, partly because of logistical challenges associated with sampling at an appropriate scale, and the large effective population sizes of species that can be sampled adequately may mask underlying genetic structure. Here we examine genetic structure in the superabundant amphipod Hirondellea gigas in the Mariana Trench at depths of 8126-10,545 m. RAD sequencing was used to identify 3182 loci containing 43,408 single nucleotide polymorphisms (SNPs) across individuals after stringent pruning of loci to prevent paralogous multicopy genomic regions being erroneously merged. Principal components analysis of SNP genotypes resolved no genetic structure between sampling locations, consistent with a signature of panmixia. However, discriminant analysis of principal components identified divergence between all sites driven by 301 outlier SNPs in 169 loci and significantly associated with latitude and depth. Functional annotation of loci identified differences between singleton loci used in analysis and paralogous loci pruned from the data set and also between outlier and nonoutlier loci, all consistent with hypotheses explaining the role of transposable elements driving genome dynamics. This study challenges the traditional perspective that highly abundant amphipods within a trench form a single panmictic population. We discuss the findings in relation to eco-evolutionary and ontogenetic processes operating in the deep sea, and highlight key challenges associated with population genetic analysis in nonmodel systems with inherent large effective population sizes and genomes.

High Prevalence of Clonal Reproduction and Low Genetic Diversity in Scutellaria floridana, a Federally Threatened Florida-Endemic Mint

The threatened mint Florida skullcap (Scutellaria floridana) is endemic to four counties in the Florida panhandle. Because development and habitat modification extirpated several historical occurrences, only 19 remain to date. To inform conservation management and delisting decisions, a comprehensive investigation of the genetic diversity and relatedness, population structure, and clonal diversity was conducted using SNP data generated by ddRAD. Compared with other Lamiaceae, we detected low genetic diversity (H_E = 0.125-0.145), low to moderate evidence of inbreeding (F_IS = -0.02-0.555), and moderate divergence (F_ST = 0.05-0.15). We identified eight populations with most of the genetic diversity, which should be protected in situ, and four populations with low genetic diversity and high clonality. Clonal reproduction in our circular plots and in 92% of the sites examined was substantial, with average clonal richness of 0.07 and 0.59, respectively. Scutellaria floridana appears to have experienced a continued decline in the number of extant populations since its listing under the Endangered Species Act; still, the combination of sexual and asexual reproduction may be advantageous for maintaining the viability of extant populations. However, the species will likely require ongoing monitoring, management, and increased public awareness to ensure its survival and effectively conserve its genetic diversity.