ipyrad: Interactive assembly and analysis of RADseq datasets

Population structure and interspecific hybridisation of two invasive blowflies (Diptera: Calliphoridae) following replicated incursions into New Zealand

Rates of biological invasion are increasing globally, with associated negative effects on native biodiversity and ecosystem services. Among other genetic processes, hybridisation can facilitate invasion by producing new combinations of genetic variation that increase adaptive potential and associated population fitness. Yet the role of hybridisation (and resulting gene flow) in biological invasion in invertebrate species is under-studied. Calliphora hilli and Calliphora stygia are blowflies proposed to have invaded New Zealand separately from Australia between 1779 and 1841, and are now widespread throughout the country. Here, we analysed genome-wide single nucleotide polymorphisms (SNPs), generating genotyping-by-sequencing data for 154 individuals collected from 24 populations across New Zealand and Australia to assess the extent of gene flow and hybridisation occurring within and between these blowflies and to better understand their overall population structure. We found that New Zealand populations of both species had weak genetic structure, suggesting high gene flow and an absence of dispersal limitations across the country. We also found evidence that interspecific hybridisation is occurring in the wild between C. hilli and C. stygia in both the native and invasive ranges, and that intraspecific admixture is occurring among populations at appreciable rates. Collectively, these findings provide new insights into the population structure of these two invasive invertebrates and highlight the potential importance of hybridisation and gene flow in biological invasion.

Go west: Population genomics reveals unexpected population fluctuations and little gene flow in Western hemisphere populations of the predatory lady beetle, Hippodamia convergens

Hippodamia convergens-the convergent lady beetle, has been used extensively in augmentative biological control of aphids, thrips, and whiteflies across its native range in North America, and was introduced into South America in the 1950s. Overwintering H. convergens populations from its native western range in the United States are commercially collected and released across its current range in the eastern USA, with little knowledge of the effectiveness of its augmentative biological control. Here we use a novel ddRADseq-based SNP/haplotype discovery approach to estimate its range-wide population diversity, differentiation, and recent evolutionary history. Our results indicate (1) significant population differentiation among eastern USA, western USA, and South American populations of H. convergens, with (2) little to no detectable recent admixture between them, despite repeated population augmentation, and (3) continued recent population size expansion across its range. These results contradict previous findings using microsatellite markers. In light of these new findings, the implications for the effectiveness of augmentative biological control using H. convergens are discussed. Additionally, because quantifying the non-target effects of augmentative biological control is a difficult problem in migratory beetles, our results could serve as a cornerstone in improving and predicting the efficacy of future releases of H. convergens across its range.

Exploring coral speciation: Multiple sympatric Stylophora pistillata taxa along a divergence continuum on the Great Barrier Reef

Understanding how biodiversity originates and is maintained are fundamental challenge in evolutionary biology. Speciation is a continuous process and progression along this continuum depends on the interplay between evolutionary forces driving divergence and forces promoting genetic homogenisation. Coral reefs are broadly connected yet highly heterogeneous ecosystems, and divergence with gene flow at small spatial scales might therefore be common. Genomic studies are increasingly revealing the existence of closely related and sympatric taxa within taxonomic coral species, but the extent to which these taxa might still be exchanging genes and sharing environmental niches is unclear. In this study, we sampled extensively across diverse habitats at multiple reefs of the Great Barrier Reef (GBR) and comprehensively examined genome-wide diversity and divergence histories within and among taxa of the Stylophora pistillata species complex. S. pistillata is one of the most abundant and well-studied coral species, yet we discovered five distinct taxa, with wide geographic ranges and extensive sympatry. Demographic modelling showed that speciation events have occurred with gene flow and that taxa are at different stages along a divergence continuum. We found significant correlations between genetic divergence and specific environmental variables, suggesting that niche partitioning may have played a role in speciation and that S. pistillata taxa might be differentially adapted to different environments. Conservation actions rely on estimates of species richness, population sizes and species ranges, which are biased if divergent taxa are lumped together. As coral reefs are rapidly degrading due to climate change, our study highlights the importance of recognising evolutionarily distinct and differentially adapted coral taxa to improve conservation and restoration efforts aiming at protecting coral genetic diversity.

Discerning structure versus speciation in phylogeographic analysis of Seepage Salamanders (Desmognathus aeneus) using demography, environment, geography, and phenotype

Numerous mechanisms can drive speciation, including isolation by adaptation, distance, and environment. These forces can promote genetic and phenotypic differentiation of local populations, the formation of phylogeographic lineages, and ultimately, completed speciation. However, conceptually similar mechanisms may also result in stabilizing rather than diversifying selection, leading to lineage integration and the long-term persistence of population structure within genetically cohesive species. Processes that drive the formation and maintenance of geographic genetic diversity while facilitating high rates of migration and limiting phenotypic differentiation may thereby result in population genetic structure that is not accompanied by reproductive isolation. We suggest that this framework can be applied more broadly to address the classic dilemma of "structure" versus "species" when evaluating phylogeographic diversity, unifying population genetics, species delimitation, and the underlying study of speciation. We demonstrate one such instance in the Seepage Salamander (Desmognathus aeneus) from the southeastern United States. Recent studies estimated up to 6.3% mitochondrial divergence and four phylogenomic lineages with broad admixture across geographic hybrid zones, which could potentially represent distinct species supported by our species-delimitation analyses. However, while limited dispersal promotes substantial isolation by distance, microhabitat specificity appears to yield stabilizing selection on a single, uniform, ecologically mediated phenotype. As a result, climatic cycles promote recurrent contact between lineages and repeated instances of high migration through time. Subsequent hybridization is apparently not counteracted by adaptive differentiation limiting introgression, leaving a single unified species with deeply divergent phylogeographic lineages that nonetheless do not appear to represent incipient species.

The Artefactual Branch Effect and Phylogenetic Conflict: Species Delimitation with Gene Flow in Mangrove Pit Vipers (Trimeresurus purpureomaculatus-erythrurus Complex)

Mangrove pit vipers of the Trimeresurus purpureomaculatus-erythrurus complex are the only species of viper known to naturally inhabit mangroves. Despite serving integral ecological functions in mangrove ecosystems, the evolutionary history, distribution, and species boundaries of mangrove pit vipers remain poorly understood, partly due to overlapping distributions, confusing phenotypic variations, and the lack of focused studies. Here, we present the first genomic study on mangrove pit vipers and introduce a robust hypothesis-driven species delimitation framework that considers gene flow and phylogenetic uncertainty in conjunction with a novel application of a new class of speciation-based delimitation model implemented through the program Delineate. Our results showed that gene flow produced phylogenetic conflict in our focal species and substantiates the artefactual branch effect where highly admixed populations appear as divergent nonmonophyletic lineages arranged in a stepwise manner at the basal position of clades. Despite the confounding effects of gene flow, we were able to obtain unequivocal support for the recognition of a new species based on the intersection and congruence of multiple lines of evidence. This study demonstrates that an integrative hypothesis-driven approach predicated on the consideration of multiple plausible evolutionary histories, population structure/differentiation, gene flow, and the implementation of a speciation-based delimitation model can effectively delimit species in the presence of gene flow and phylogenetic conflict.

Conservation Genetics of the Endangered Lompoc Yerba Santa (Eriodictyon capitatum Eastw., Namaceae), including Phylogenomic Insights into the Evolution of Eriodictyon

Eriodictyon capitatum (Namaceae) is a narrowly distributed shrub endemic to western Santa Barbara County, where it is known from only 10 extant California Natural Diversity Database element occurrences (EOs). Owing to low numbers of plants in nature, a limited overall extent, and multiple current threats, E. capitatum is listed as Endangered under the Federal Endangered Species Act and as Rare under the California Native Plant Protection Act. In the present study, high-throughput DNA sequence data were analyzed to investigate genetic diversity within and among all accessible EOs; to determine the extent of genetic isolation among EOs; to examine clonality within EOs; and to examine the taxonomic circumscriptions of E. capitatum, E. altissimum, E. angustifolium, and E. californicum through phylogenomic analysis. Population genetic analyses of E. capitatum reveal a pattern of strong genetic differentiation by location/EO. The clonality assessment shows that certain small EOs may support relatively few multilocus genotypes. The phylogenomic analyses strongly support the present-day taxonomic circumscriptions of both E. altissimum and E. capitatum, showing them to be reciprocally monophyletic and sister with strong support. Taken together, these results paint a picture of an evolutionarily and morphologically distinct species known from relatively few, genetically isolated stations.

Reference-Based Restriction-Site-Associated DNA Sequencing Data Are Useful for Species Delineation in a Recently Diverged Asexually Reproducing Species Complex (Parmeliaceae, Ascomycota)

Cryptic species are common in lichen-forming fungi and have been reported from different genera in the most speciose family, Parmeliaceae. Herein, we address species delimitation in a group of mainly asexually reproducing Parmelina species. The morphologically distinct P. pastillifera was previously found nested within a morphologically circumscribed P. tiliacea based on several loci. However, these studies demonstrated a relatively high genetic diversity within P. tiliacea sensu lato. Here, we revisit the species delimitation in the group by analyzing single-nucleotide polymorphisms (SNPs) through genome-wide assessment using Restriction-Site-Associated sequencing and population genomic methods. Our data support previous studies and provide further insight into the phylogenetic relationships of the four clades found within the complex. Based on the evidence suggesting a lack of gene flow among the clades, we recognize the four clades as distinct species, P. pastillifera and P. tiliacea sensu stricto, and two new species, P. clandestina sp. nov. and P. mediterranea sp. nov.

Riverscape community genomics: A comparative analytical approach to identify common drivers of spatial structure

Genetic differentiation among local groups of individuals, that is, genetic β-diversity, is a key component of population persistence related to connectivity and isolation. However, most genetic investigations of natural populations focus on a single species, overlooking opportunities for multispecies conservation plans to benefit entire communities in an ecosystem. We present an approach to evaluate genetic β-diversity within and among many species and demonstrate how this riverscape community genomics approach can be applied to identify common drivers of genetic structure. Our study evaluated genetic β-diversity in 31 co-distributed native stream fishes sampled from 75 sites across the White River Basin (Ozarks, USA) using SNP genotyping (ddRAD). Despite variance among species in the degree of genetic divergence, general spatial patterns were identified corresponding to river network architecture. Most species (N = 24) were partitioned into discrete subpopulations (K = 2-7). We used partial redundancy analysis to compare species-specific genetic β-diversity across four models of genetic structure: Isolation by distance (IBD), isolation by barrier (IBB), isolation by stream hierarchy (IBH), and isolation by environment (IBE). A significant proportion of intraspecific genetic variation was explained by IBH (x̄ = 62%), with the remaining models generally redundant. We found evidence for consistent spatial modularity in that gene flow is higher within rather than between hierarchical units (i.e., catchments, watersheds, basins), supporting the generalization of the stream hierarchy model. We discuss our conclusions regarding conservation and management and identify the 8-digit hydrologic unit (HUC) as the most relevant spatial scale for managing genetic diversity across riverine networks.

Reduced representation approaches produce similar results to whole genome sequencing for some common phylogeographic analyses

When designing phylogeographic investigations researchers can choose to collect many different types of molecular markers, including mitochondrial genes or genomes, SNPs from reduced representation protocols, large sequence capture data sets, and even whole genomes. Given that the statistical power and accuracy of various analyses are expected to differ depending on both the type of marker and the amount of data collected, an exploration of the variance across methodological results as a function of marker type should provide valuable information to researchers. Here we collect mitochondrial Cytochrome b sequences, whole mitochondrial genomes, single nucleotide polymorphisms (SNP)s isolated using a genotype by sequencing (GBS) protocol, sequences from ultraconserved elements, and low-coverage nuclear genomes from the North American water vole (Microtus richardsoni). We estimate genetic distances, population genetic structure, and historical demography using data from each of these datasets and compare the results across markers. As anticipated, the results exhibit differences across marker types, particularly in terms of the resolution offered by different analyses. A cost-benefit analysis indicates that SNPs collected using a GBS protocol are the most cost-effective molecular marker, with inferences that mirror those collected from the whole genome data at a fraction of the cost per sample.

Divergent selection and climate adaptation fuel genomic differentiation between sister species of Sphagnum (peat moss)

BACKGROUND AND AIMS

New plant species can evolve through the reinforcement of reproductive isolation via local adaptation along habitat gradients. Peat mosses (Sphagnaceae) are an emerging model system for the study of evolutionary genomics and have well-documented niche differentiation among species. Recent molecular studies have demonstrated that the globally distributed species Sphagnum magellanicum is a complex of morphologically cryptic lineages that are phylogenetically and ecologically distinct. Here, we describe the architecture of genomic differentiation between two sister species in this complex known from eastern North America: the northern S. diabolicum and the largely southern S. magniae.

METHODS

We sampled plant populations from across a latitudinal gradient in eastern North America and performed whole genome and restriction-site associated DNA sequencing. These sequencing data were then analyzed computationally.

KEY RESULTS

Using sliding-window population genetic analyses we find that differentiation is concentrated within 'islands' of the genome spanning up to 400 kb that are characterized by elevated genetic divergence, suppressed recombination, reduced nucleotide diversity and increased rates of non-synonymous substitution. Sequence variants that are significantly associated with genetic structure and bioclimatic variables occur within genes that have functional enrichment for biological processes including abiotic stress response, photoperiodism and hormone-mediated signalling. Demographic modelling demonstrates that these two species diverged no more than 225 000 generations ago with secondary contact occurring where their ranges overlap.

CONCLUSIONS

We suggest that this heterogeneity of genomic differentiation is a result of linked selection and reflects the role of local adaptation to contrasting climatic zones in driving speciation. This research provides insight into the process of speciation in a group of ecologically important plants and strengthens our predictive understanding of how plant populations will respond as Earth's climate rapidly changes.

Scrutinising an inscrutable bark-nesting ant: Exploring cryptic diversity in the Rhopalomastix javana (Hymenoptera: Formicidae) complex using DNA barcodes, genome-wide MIG-seq and geometric morphometrics

Overlooking cryptic species diversity has grave implications on assessments of climate change impacts on biodiversity, ecosystems and organismal populations. Discriminating between cryptic species has long been challenging even for seasoned taxonomists, as interspecies morphological differences are often indiscernible by visual observation. Multi-disciplinary methods involving genetic analyses in conjunction with quantitative morphological data, should therefore be used to investigate boundaries between cryptic species. We adopted an integrated approach combining analyses of mitochondrial COI barcodes, a genome-wide dataset obtained via multiplexed inter-simple sequence repeats (ISSRs) genotyping by sequencing (MIG-seq), and geometric morphometrics to investigate species divergences in the inscrutable Rhopalomastix javana species complex. Objective clustering of COI suggested five putative molecular species units divergent from each other by thresholds within 4.2-10.6% uncorrected pairwise distance. Phylogenetic analyses based on concatenated MIG-seq data also recovered and strongly supported the monophyly of five major lineages in agreement with COI clusters. Co-ancestry analyses based on MIG-seq data using fineRADstructure resolved variable patterns of admixture linked to geography, and potential genetic drift within some putative species. Geometric morphometric analyses of specimen images further detected statistically significant differences in at least one of three anatomical aspects (Head, Meso, Profile) between all pairs of putative species. Head shape (full-face view) was determined to be the most informative character for species diagnosis, with relatively high classification accuracy. Thin-plate spline deformation grids highlighted areas of high variation between species in each shape for deeper taxonomic scrutiny. The presence of species from multiple distinct lineages existing in near-sympatry firmly demonstrates that R. javana comprises more than one closely-related species, but exact species boundaries are difficult to ascertain. Differences in elevation and its associated abiotic effects on ant adaptations and reproductive phenology may contribute to restricting gene flow and maintaining species boundaries between sympatric populations of the R. javana complex. We further assess the advantages and limitations of geometric morphometrics as a taxonomic tool. Despite its drawbacks, our combined approach has helped draw important insights on cryptic diversity in R. javana, and also identified gaps of knowledge that await address. Results from this study will inform and prime future in-depth taxonomic investigation on the R. javana complex, including formal descriptions and establishment of the five putative species.

Differences in the genomic diversity, structure, and inbreeding patterns in wild and managed populations of Agave potatorum Zucc. used in the production of Tobalá mezcal in Southern Mexico

Agave potatorum Zucc. locally known as Tobalá, is an important species for mezcal production. It is a perennial species that takes 10 to 15 years to reach reproductive age. Because of high demand of Tobalá mezcal and the slow maturation of the plants, its wild populations have been under intense anthropogenic pressure. The main objective of this study was to estimate the genome-wide diversity in A. potatorum and determine if the type of management has had any effect on its diversity, inbreeding and structure. We analyzed 174 individuals (105 wild, 42 cultivated and 27 from nurseries) from 34 sites with a reduced representation genomic method (ddRADseq), using 14,875 SNPs. The diversity measured as expected heterozygosity was higher in the nursery and wild plants than in cultivated samples. We did not find private alleles in the cultivated and nursery plants, which indicates that the individuals under management recently derived from wild populations, which was supported by higher gene flow estimated from wild populations to the managed plants. We found low but positive levels of inbreeding (FIS = 0.082), probably related to isolation of the populations. We detected low genetic differentiation among populations (FST = 0.0796), with positive and significant isolation by distance. The population genetic structure in the species seems to be related to elevation and ecology, with higher gene flow among populations in less fragmented areas. We detected an outlier locus related to the recognition of pollen, which is also relevant to self-incompatibility protein (SI). Due to seed harvest and long generation time, the loss of diversity in A. potatorum has been gradual and artificial selection and incipient management have not yet caused drastic differences between cultivated and wild plants. Also, we described an agroecological alternative to the uncontrolled extraction of wild individuals.

Host-Shift Speciation Proceeded with Gene Flow in Algae Covering Shells

AbstractHost shifts represent the advancement of a novel niche and often lead to speciation in symbionts. However, its mechanisms are not well understood. Here, we focused on the alga Pseudocladophora conchopheria growing on the shells of intertidal snails. Previous surveys have shown that the alga has host specificity-only attaching to the shell of Lunella correensis-but we discovered that the alga attaches to the shells of multiple sympatric snails. A genome-wide single-nucleotide polymorphism analysis (MIG-seq) was performed to determine whether host-associated speciation occurred in the algae. As a result, there was no gene flow or limited gene flow among the algae from different hosts, and some algae were genetically differentiated among hosts. In addition, the demographic estimate revealed that speciation with gene flow occurred between the algae from different hosts. Therefore, these results support the idea that host-shift speciation gradually proceeded with gene flow in the algae, providing insight into the early evolution of host shifts.

Lack of host specialization despite selective host use in brood parasitic cuckoo catfish

Host-parasite dynamics involve coevolutionary arms races, which may lead to host specialization and ensuing diversification. Our general understanding of the evolution of host specialization in brood parasites is compromised by a restricted focus on bird and insect lineages. The cuckoo catfish (Synodontis multipunctatus) is an obligate parasite of parental care of mouthbrooding cichlids in Lake Tanganyika. Given the ecological and taxonomic diversity of mouthbrooding cichlids in the lake, we hypothesized the existence of sympatric host-specific lineages in the cuckoo catfish. In a sample of 779 broods from 20 cichlid species, we found four species parasitized by cuckoo catfish (with prevalence of parasitism of 2%-18%). All parasitized cichlids were from the tribe Tropheini, maternal mouthbrooders that spawn over a substrate (rather than in open water). Phylogenetic analysis based on genomic (ddRAD sequencing) and mitochondrial (Dloop) data from cuckoo catfish embryos showed an absence of host-specific lineages. This was corroborated by analyses of genetic structure and co-ancestry matrix. Within host species, parasitism was not associated with any individual characteristic we recorded (parent size, water depth), but was costly as parasitized parents carried smaller clutches of their own offspring. We conclude that the cuckoo catfish is an intermediate generalist and discuss costs, benefits and constraints of host specialization in this species and brood parasites in general.

Speciation in Coastal Basins Driven by Staggered Headwater Captures: Dispersal of a Species Complex, Leporinus bahiensis, as Revealed by Genome-wide SNP Data

Past sea level changes and geological instability along watershed boundaries have largely influenced fish distribution across coastal basins, either by dispersal via palaeodrainages now submerged or by headwater captures, respectively. Accordingly, the South American Atlantic coast encompasses several small and isolated drainages that share a similar species composition, representing a suitable model to infer historical processes. Leporinus bahiensis is a freshwater fish species widespread along adjacent coastal basins over narrow continental shelf with no evidence of palaeodrainage connections at low sea level periods. Therefore, this study aimed to reconstruct its evolutionary history to infer the role of headwater captures in the dispersal process. To accomplish this, we employed molecular-level phylogenetic and population structure analyses based on Sanger sequences (5 genes) and genome-wide SNP data. Phylogenetic trees based on Sanger data were inconclusive, but SNPs data did support the monophyletic status of L. bahiensis. Both COI and SNP data revealed structured populations according to each hydrographic basin. Species delimitation analyses revealed from 3 (COI) to 5 (multilocus approach) MOTUs, corresponding to the sampled basins. An intricate biogeographic scenario was inferred and supported by Approximate Bayesian Computation (ABC) analysis. Specifically, a staggered pattern was revealed and characterized by sequential headwater captures from basins adjacent to upland drainages into small coastal basins at different periods. These headwater captures resulted in dispersal throughout contiguous coastal basins, followed by deep genetic divergence among lineages. To decipher such recent divergences, as herein represented by L. bahiensis populations, we used genome-wide SNPs data. Indeed, the combined use of genome-wide SNPs data and ABC method allowed us to reconstruct the evolutionary history and speciation of L. bahiensis. This framework might be useful in disentangling the diversification process in other neotropical fishes subject to a reticulate geological history.

Phylogenomic and morphological data reveal hidden patterns of diversity in the national tree of Brazil, Paubrasilia echinata

PREMISE

Paubrasilia echinata (common names, pau brasil, brazilwood) is the national tree of Brazil and an endangered species endemic to the Brazilian Atlantic Forest. Over its wide distribution of 2000 km, its leaflets morphology exhibits extensive plasticity. Three morphotypes are commonly identified based on leaf size, but it is unclear if they represent distinct taxa or a single polymorphic species. This study aims to clarify the taxonomic position of the three morphotypes to inform conservation decisions.

METHODS

A morphometric study of leaf characters of herbarium specimens was coupled with genetic analyses using genotype-by-sequencing data. We used maximum-likelihood and coalescent methods to evaluate the phylogenetic and population structure of the species. We compared these with a morphological dendrogram built from hierarchical clustering.

RESULTS

Two of the three morphotypes formed separately evolving lineages, the third morphotype formed two geographically separate lineages, and northern trees with intermediate leaf morphology formed a separate fifth lineage. Leaflet size varied by over 35-fold, and although morphological clustering generally matched the genetic patterns, there were some overlaps, highlighting the cryptic diversity within this group.

CONCLUSIONS

Our genetic and morphological results provide some evidence that cultivated trees from different states in Brazil seem to have a limited genetic origin and do not reflect the broader genetic and geographical diversity of the species. As a result, more care is likely needed to preserve the overall genomic diversity of this endangered and iconic species.

Interspecific gene flow obscures phylogenetic relationships in an important insect pest species complex

As genomic data proliferates, the prevalence of post-speciation gene flow is making species boundaries and relationships increasingly ambiguous. Although current approaches inferring fully bifurcating phylogenies based on concatenated datasets provide simple and robust answers to many species relationships, they may be inaccurate because the models ignore inter-specific gene flow and incomplete lineage sorting. To examine the potential error resulting from ignoring gene flow, we generated both a RAD-seq and a 500 protein-coding loci highly multiplexed amplicon (HiMAP) dataset for a monophyletic group of 12 species defined as the Bactrocera dorsalis sensu lato clade. With some of the world's worst agricultural pests, the taxonomy of the B. dorsalis s.l. clade is important for trade and quarantines. However, taxonomic confusion confounds resolution due to intra- and interspecific phenotypic variation and convergence, mitochondrial introgression across half of the species, and viable hybrids. We compared the topological convergence of our datasets using concatenated phylogenetic and various multispecies coalescent approaches, some of which account for gene flow. All analyses agreed on species delimitation, but there was incongruence between species relationships. Under concatenation, both datasets suggest identical species relationships with mostly high statistical support. However, multispecies coalescent and multispecies network approaches suggest markedly different hypotheses and detected significant gene flow. We suggest that the network approaches are likely more accurate because gene flow violates the assumptions of the concatenated phylogenetic analyses, but the data-reductive requirements of network approaches resulted in reduced statistical support and could not unambiguously resolve gene flow directions. Our study highlights the importance of testing for gene flow, particularly with phylogenomic datasets, even when concatenated approaches receive high statistical support.

Population genomics reveals how 5 ka of human occupancy led the Lima leaf-toed gecko (Phyllodactylus sentosus) to the brink of extinction

Small species with high home fidelity, high ecological specialization or low vagility are particularly prone to suffer from habitat modification and fragmentation. The Lima leaf-toed gecko (Phyllodactylus sentosus) is a critically endangered Peruvian species that shelters mostly in pre-Incan archeological areas called huacas, where the original environmental conditions are maintained. We used genotyping by sequencing to understand the population genomic history of P. sentosus. We found low genetic diversity (He 0.0406-0.134 and nucleotide diversity 0.0812-0.145) and deviations of the observed heterozygosity relative to the expected heterozygosity in some populations (Fis - 0.0202 to 0.0187). In all analyses, a clear population structuring was observed that cannot be explained by isolation by distance alone. Also, low levels of historical gene flow were observed between most populations, which decreased as shown in contemporary migration rate analysis. Demographic inference suggests these populations experienced bottleneck events during the last 5 ka. These results indicate that habitat modification since pre-Incan civilizations severely affected these populations, which currently face even more drastic urbanization threats. Finally, our predictions show that this species could become extinct in a decade without further intervention, which calls for urgent conservation actions being undertaken.

Parallel evolution, atavism, and extensive introgression explain the radiation of Epimedium sect. Diphyllon (Berberidaceae) in southern East Asia

East Asia is the richest region of plant biodiversity in the northern temperate zone, and its radiation provides key insights for understanding rapid speciation, including evolutionary patterns and processes. However, it is challenging to investigate the recent evolutionary radiation among plants because of the lack of genetic divergence, phenotypic convergence, and interspecific gene flow. Epimedium sect. Diphyllon is a rarely studied plant lineage endemic to East Asia, especially highly diversified in its southern part. In this study, we report a robust phylogenomic analysis based on genotyping-by-sequencing data of this lineage. The results revealed a clear biogeographic pattern for Epimedium sect. Diphyllon with recognition into two major clades corresponding to the Sino-Himalayan and Sino-Japanese subkingdoms of East Asian Flora and rapid diversification of the extant species dated to the Pleistocene. Evolutionary radiation of Epimedium sect. Diphyllon is characterized by recent and predominant parallel evolution and atavism between the two subkingdom regions, with extensive reticulating hybridization within each region during the course of diversification in southern East Asia. A parallel-atavism-introgression hypothesis is referred to in explaining the radiation of plant diversity in southern East Asia, which represents a potential model for the rapid diversification of plants under global climate cooling in the late Tertiary. Our study advances our understanding of the evolutionary processes of plant radiation in East Asia as well as in other biodiversity hotspot regions.

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.

The holocentric chromosome microevolution: From phylogeographic patterns to genomic associations with environmental gradients

Geographic isolation and chromosome evolution are two of the major drivers of diversification in eukaryotes in general, and specifically, in plants. On one hand, range shifts induced by Pleistocene glacial oscillations deeply shaped the evolutionary trajectories of species in the Northern Hemisphere. On the other hand, karyotype variability within species or species complexes may have adaptive potential as different karyotypes may represent different recombination rates and linkage groups that may be associated with locally adapted genes or supergenes. Organisms with holocentric chromosomes are ideal to study the link between local adaptation and chromosome evolution, due to their high cytogenetic variability, especially when it seems to be related to environmental variation. Here, we integrate the study of the phylogeography, chromosomal evolution and ecological requirements of a plant species complex distributed in the Western Euro-Mediterranean region (Carex gr. laevigata, Cyperaceae). We aim to clarify the relative influence of these factors on population differentiation and ultimately on speciation. We obtained a well-resolved RADseq phylogeny that sheds light on the phylogeographic patterns of molecular and chromosome number variation, which are compatible with south-to-north postglacial migration. In addition, landscape genomics analyses identified candidate loci for local adaptation, and also strong significant associations between the karyotype and the environment. We conclude that karyotype distribution in C. gr. laevigata has been constrained by both range shift dynamics and local adaptation. Our study demonstrates that chromosome evolution may be responsible, at least partially, for microevolutionary patterns of population differentiation and adaptation in Carex.

Genome-wide association study to identify novel loci and genes for Fusarium root rot resistance in sweet potato using genotyping-by-sequencing

Fusarium root rot, caused by Fusarium solani, is a major post-harvest disease in sweet potatoes (Ipomoea batatas (L.) Lam.). An effective strategy for controlling this disease is the development of resistant varieties. In this study, a genome-wide association study (GWAS) was conducted on 96 sweet potato genotypes to identify novel candidate loci and dissect the genetic basis of Fusarium root rot resistance. Genotyping was performed using genotyping-by-sequencing (GBS), and 44,255 SNPs were identified after filtering. The genotypes (n = 96) were evaluated through resistance tests in 2021 and 2022, separately and combined. The GWAS identified two significant SNP markers (LG3_22903756 and LG4_2449919) on chromosomes 3 and 4 associated with Fusarium root rot resistance, respectively. Lesion length showed significant differences between homozygous A and G alleles of LG3_22903756, which can potentially be used to develop molecular markers for selecting accessions resistant to Fusarium root rot. Expression analysis of 11 putative genes flanking the significant SNPs revealed the alteration in the expression of nine genes, indicating their possible involvement in Fusarium root rot resistance. The results of this study will aid in the marker-assisted selection and functional analysis of candidate genes for Fusarium root rot resistance in sweet potatoes.

Clonality, local population structure and gametophyte sex ratios in cryptic species of the Sphagnum magellanicum complex

BACKGROUND AND AIMS

Sphagnum (peatmoss) comprises a moss (Bryophyta) clade with ~300-500 species. The genus has unparalleled ecological importance because Sphagnum-dominated peatlands store almost a third of the terrestrial carbon pool and peatmosses engineer the formation and microtopography of peatlands. Genomic resources for Sphagnum are being actively expanded, but many aspects of their biology are still poorly known. Among these are the degree to which Sphagnum species reproduce asexually, and the relative frequencies of male and female gametophytes in these haploid-dominant plants. We assess clonality and gametophyte sex ratios and test hypotheses about the local-scale distribution of clones and sexes in four North American species of the S. magellanicum complex. These four species are difficult to distinguish morphologically and are very closely related. We also assess microbial communities associated with Sphagnum host plant clones and sexes at two sites.

METHODS

Four hundred and five samples of the four species, representing 57 populations, were subjected to restriction site-associated DNA sequencing (RADseq). Analyses of population structure and clonality based on the molecular data utilized both phylogenetic and phenetic approaches. Multi-locus genotypes (genets) were identified using the RADseq data. Sexes of sampled ramets were determined using a molecular approach that utilized coverage of loci on the sex chromosomes after the method was validated using a sample of plants that expressed sex phenotypically. Sex ratios were estimated for each species, and populations within species. Difference in fitness between genets was estimated as the numbers of ramets each genet comprised. Degrees of clonality [numbers of genets/numbers of ramets (samples)] within species, among sites, and between gametophyte sexes were estimated. Sex ratios were estimated for each species, and populations within species. Sphagnum-associated microbial communities were assessed at two sites in relation to Sphagnum clonality and sex.

KEY RESULTS

All four species appear to engage in a mixture of sexual and asexual (clonal) reproduction. A single ramet represents most genets but two to eight ramets were dsumbers ansd text etected for some genets. Only one genet is represented by ramets in multiple populations; all other genets are restricted to a single population. Within populations ramets of individual genets are spatially clustered, suggesting limited dispersal even within peatlands. Sex ratios are male-biased in S. diabolicum but female-biased in the other three species, although significantly so only in S. divinum. Neither species nor males/females differ in levels of clonal propagation. At St Regis Lake (NY) and Franklin Bog (VT), microbial community composition is strongly differentiated between the sites, but differences between species, genets and sexes were not detected. Within S. divinum, however, female gametophytes harboured two to three times the number of microbial taxa as males.

CONCLUSIONS

These four Sphagnum species all exhibit similar reproductive patterns that result from a mixture of sexual and asexual reproduction. The spatial patterns of clonally replicated ramets of genets suggest that these species fall between the so-called phalanx patterns, where genets abut one another but do not extensively mix because of limited ramet fragmentation, and the guerrilla patterns, where extensive genet fragmentation and dispersal result in greater mixing of different genets. Although sex ratios in bryophytes are most often female-biased, both male and female biases occur in this complex of closely related species. The association of far greater microbial diversity for female gametophytes in S. divinum, which has a female-biased sex ratio, suggests additional research to determine if levels of microbial diversity are consistently correlated with differing patterns of sex ratio biases.

Urbanization drives adaptive evolution in a Neotropical bird

Urbanization has dramatic impacts on natural habitats and such changes may potentially drive local adaptation of urban populations. Behavioral change has been specifically shown to facilitate the fast adaptation of birds to changing environments, but few studies have investigated the genetic mechanisms of this process. Such investigations could provide insights into questions about both evolutionary theory and management of urban populations. In this study, we investigated whether local adaptation has occurred in urban populations of a Neotropical bird species, Coereba flaveola, specifically addressing whether observed behavioral adaptations are correlated to genetic signatures of natural selection. To answer this question, we sampled 24 individuals in urban and rural environments, and searched for selected loci through a genome-scan approach based on RADseq genomic data, generated and assembled using a reference genome for the species. We recovered 46 loci as putative selection outliers, and 30 of them were identified as associated with biological processes possibly related to urban adaptation, such as the regulation of energetic metabolism, regulation of genetic expression, and changes in the immunological system. Moreover, genes involved in the development of the nervous system showed signatures of selection, suggesting a link between behavioral and genetic adaptations. Our findings, in conjunction with similar results in previous studies, support the idea that cities provide a similar selective pressure on urban populations and that behavioral plasticity may be enhanced through genetic changes in urban populations.

Evidence of introgression, ecological divergence and adaptation in Asterias sea stars

Hybrid zones are important windows into the evolutionary dynamics of populations, revealing how processes like introgression and adaptation structure population genomic variation. Importantly, they are useful for understanding speciation and how species respond to their environments. Here, we investigate two closely related sea star species, Asterias rubens and A. forbesi, distributed along rocky European and North American coastlines of the North Atlantic, and use genome-wide molecular markers to infer the distribution of genomic variation within and between species in this group. Using genomic data and environmental niche modelling, we document hybridization occurring between northern New England and the southern Canadian Maritimes. We investigate the factors that maintain this hybrid zone, as well as the environmental variables that putatively drive selection within and between species. We find that the two species differ in their environmental niche breadth; Asterias forbesi displays a relatively narrow environmental niche while conversely, A. rubens has a wider niche breadth. Species distribution models accurately predict hybrids to occur within environmental niche overlap, thereby suggesting environmental selection plays an important role in the maintenance of the hybrid zone. Our results imply that the distribution of genomic variation in North Atlantic sea stars is influenced by the environment, which will be crucial to consider as the climate changes.

Reduced-representation libraries in insect genetics

Genotyping-by-sequencing of reduced-representation libraries has ushered in an era where genome-wide data can be gotten for any species. Here, I review research on this topic during the last two years, report meta-analysis of the results, and discuss analysis methods and issues. Scanning the literature from 2021 to 2022 identified 21 papers, the majority of which were on population differences, including local adaptation and migration, but several papers were on genetic maps and their use in assembly scaffolding or analysis of quantitative trait loci, on the origin of incursions of pest insects, or on infection rates of a pathogen in a disease vector. The research reviewed includes 33 species from 25 families and 11 orders. Meta-analysis showed that less than 16%, and most often, less than 1% of the genome was implicated in local adaptation and that the number of adaptive loci correlated with genetic divergence among populations.

Loss of species and genetic diversity during colonization: Insights from acanthocephalan parasites in northern European seals

Studies on host-parasite systems that have experienced distributional shifts, range fragmentation, and population declines in the past can provide information regarding how parasite community richness and genetic diversity will change as a result of anthropogenic environmental changes in the future. Here, we studied how sequential postglacial colonization, shifts in habitat, and reduced host population sizes have influenced species richness and genetic diversity of Corynosoma (Acanthocephala: Polymorphidae) parasites in northern European marine, brackish, and freshwater seal populations. We collected Corynosoma population samples from Arctic, Baltic, Ladoga, and Saimaa ringed seal subspecies and Baltic gray seals, and then applied COI barcoding and triple-enzyme restriction-site associated DNA (3RAD) sequencing to delimit species, clarify their distributions and community structures, and elucidate patterns of intraspecific gene flow and genetic diversity. Our results showed that Corynosoma species diversity reflected host colonization histories and population sizes, with four species being present in the Arctic, three in the Baltic Sea, two in Lake Ladoga, and only one in Lake Saimaa. We found statistically significant population-genetic differentiation within all three Corynosoma species that occur in more than one seal (sub)species. Genetic diversity tended to be high in Corynosoma populations originating from Arctic ringed seals and low in the landlocked populations. Our results indicate that acanthocephalan communities in landlocked seal populations are impoverished with respect to both species and intraspecific genetic diversity. Interestingly, the loss of genetic diversity within Corynosoma species seems to have been less drastic than in their seal hosts, possibly due to their large local effective population sizes resulting from high infection intensities and effective intra-host population mixing. Our study highlights the utility of genomic methods in investigations of community composition and genetic diversity of understudied parasites.

Rearranging and completing the puzzle: Phylogenomic analysis of bark beetles Dendroctonus reveals new hypotheses about genus diversification

Studies carried out on bark beetles within Dendroctonus have been extensive and revealed diverse information in different areas of their natural history, taxonomy, evolution, and interactions, among others. Despite these efforts, phylogenetic hypotheses have remained obscured mainly due to limited information analyzed (taxonomic, gene sampling, or both) in studies focused on obtaining evolutionary hypotheses for this genus. With the aim of filling these gaps in the evolutionary history for Dendroctonus, we analyzed ∼1800 loci mapped to a reference genome obtained for 20 of the 21 species recognized to date, minimizing the impact of missing information and improving the assumption of orthology in a phylogenomic framework. We obtained congruent phylogenetic topologies from two phylogenomic inference strategies: loci concatenation (ML framework) and a multispecies coalescent model (MSC) through the analysis of site pattern frequencies (SNPs). Dendroctonus is composed of two major clades (A and B), each containing five and four subclades, respectively. According to our divergence dating analysis, the MRCA for Dendroctonus dates back to the early Eocene, while the MRCA for each major clade diverged in the mid-Eocene. Interestingly, most of the speciation events of extant species occurred during the Miocene, which could be correlated with the diversification of pine trees (Pinus). The MRCA for Dendroctonus inhabited large regions of North America, with all ancestors and descendants of clade A having diversified within this region. The Mexican Transition Zone is important in the diversification processes for the majority of clade A species. For clade B, we identified two important colonization events to the Old World from America: the first in the early Oligocene from the Arctic to Asia (via Beringia), and the second during the Miocene from the Arctic-Western-Alleghany region to Europe and Siberia (also via Beringia). Our genomic analyses also supported the existence of hidden structured lineages within the frontalis complex, and also that D. beckeri represent a lineage independent from D. valens, as previously suggested. The information presented here updates the knowledge concerning the diversification of a genus with remarkable ecological and economic importance.

Intra-island patterns of cryptic diversity within an oceanic island: Insights from the differentiation of Schizidium Verhoeff, 1901 (Oniscidea, Armadillidiidae) within Cyprus, with descriptions of two new species

Cyprus is a Mediterranean island of volcanic origin isolated for at least 5.3 Myr from surrounding continental areas. The present study focuses on the diversification of the isopod genus Schizidium within the island, including also specimens from surrounding continental areas. The genus Schizidium sensu lato is probably non monophyletic, comprising 26 species distributed from Greece to Iran. Up to date the only representative of the genus reported from Cyprus was Schizidium fissum. Aiming to investigate the patterns of genetic diversity within the focal island, to evaluate the morphology-based taxonomy of the species in the genus, and at the same time to explore phylogenetic relationships with mainland populations, we applied genome-wide ddRADseq as well as Sanger sequencing targeting three mitochondrial (16S, COI and 12S) and the nuclear NaK loci. Results of phylogenetic analyses support the existence of two distinct epigean Schizidium clades with well-defined geographic boundaries that conform to the known paleogeography of Cyprus, plus one endogean clade with restricted distribution within the island. Genetic data and morphology corroborate the assignment of this latter endogean clade to a new species, Schizidium myrrae n. sp. The two epigean clades are also considered as distinct species, one corresponding to the known S. fissum at the eastern part of the island (Pentadaktylos massif) and the other to the newly described S. christosi n. sp. distributed along the western part of the island (Troodos massif). Even though detailed examination of many specimens could not retrieve any morphological differences among representatives of these two clades, clado-chronological analysis indicates a long isolation between them, estimated at ∼ 9 Mya, as well as the sharing of a common ancestor with S. tiberianum from Israel at ∼ 15 Mya. Hence, we can consider these epigean Schizidium species as one more case of cryptic diversity on Cyprus, exhibiting similar patterns with the recently described case in the genus Armadillo.

Genomic analysis of Tibetan ground tits identifies molecular adaptations associated with cooperative breeding

Cooperative breeding is a sophisticated altruistic social behavior that helps social animals to adapt to harsh environments. The Tibetan ground tit, Pseudopodoces humilis, is a high-altitude bird endemic to the Tibetan plateau. Recently, it has become an exciting system for studying the evolution of facultative cooperative breeding. To test for molecular adaptations associated with cooperative breeding, we resequenced the whole genome of ground tits from 6 wild populations that display remarkable variation in the frequency of cooperative breeding. Population structure analyses showed that the 6 populations were divided into 4 lineages, which is congruent with the major geographical distribution of the sampling sites. Using genome-wide selective sweep analysis, we identified putative positively selected genes (PSGs) in groups of tits that displayed high and low cooperative breeding rates. The total number of PSGs varied from 146 to 722 in high cooperative breeding rate populations, and from 272 to 752 in low cooperative breeding rate populations. Functional enrichment analysis of these PSGs identified several significantly enriched ontologies related to oxytocin signaling, estrogen signaling, and insulin secretion. PSGs involved in these functional ontologies suggest that molecular adaptations in hormonal regulation may have played important roles in shaping the evolution of cooperative breeding in the ground tit. Taken together, our study provides candidate genes and functional ontologies involved in molecular adaptations associated with cooperative breeding in Tibetan ground tits, and calls for a better understanding of the genetic roles in the evolution of cooperative breeding.

Phylogenetic discordance and integrative species delimitation in the Mammillaria haageana species complex (Cactaceae)

Species complexes consist of very close phylogenetic relatives, where morphological similarities make it difficult to distinguish between them using traditional taxonomic methods. Here, we focused on the long-standing challenge of species delimitation in the Mammillaria haageana complex, a group that presents great morphological diversity that makes its taxonomy a puzzle. Our work integrates genomic, morphological, and ecological data to establish the taxonomic limits in the M. haageana complex, and we also studied the evolutionary relationships with the remainder of the M. ser. Supertextae species. Our genetic analyses, as well as morphological and ecological evidence, led us to propose that the M. haageana complex is made up of six distinct entities (M. acultzingensis, M. conspicua, M. haageana, M. lanigera, M. meissneri, and M. san-angelensis), mainly as a result of ecological speciation. A recent taxonomic proposal considered these taxa as a single species; therefore, we propose their recognition at the species level. Our results also show a high level of incomplete lineage sorting rather than reticulation, which is especially likely in recently diverged species such as those comprising M. ser. Supertextae. The species hypotheses proposed here may be useful in future extinction risk assessments and conservation strategies.

Optimizing ddRAD sequencing for population genomic studies with ddgRADer

Double-digest Restriction-site Associated DNA sequencing (ddRADseq) is widely used to generate genomic data for non-model organisms in evolutionary and ecological studies. Along with affordable paired-end sequencing, this method makes population genomic analyses more accessible. However, multiple factors should be considered when designing a ddRADseq experiment, which can be challenging for new users. The generated data often suffer from substantial read overlaps and adaptor contamination, severely reducing sequencing efficiency and affecting data quality. Here, we analyse diverse datasets from the literature and carry out controlled experiments to understand the effects of enzyme choice and size selection on sequencing efficiency. The empirical data reveal that size selection is imprecise and has limited efficacy. In certain scenarios, a substantial proportion of short fragments pass below the lower size-selection cut-off resulting in low sequencing efficiency. However, enzyme choice can considerably mitigate inadvertent inclusion of these shorter fragments. A simple model based on these experiments is implemented to predict the number of genomic fragments generated after digestion and size selection, number of SNPs genotyped, number of samples that can be multiplexed and the expected sequencing efficiency. We developed ddgRADer - http://ddgrader.haifa.ac.il/ - a user-friendly webtool and incorporated these calculations to aid in ddRADseq experimental design while optimizing sequencing efficiency. This tool can also be used for single enzyme protocols such as Genotyping-by-Sequencing. Given user-defined study goals, ddgRADer recommends enzyme pairs and allows users to compare and choose enzymes and size-selection criteria. ddgRADer improves the accessibility and ease of designing ddRADseq experiments and increases the probability of success of the first population genomic study conducted in labs with no prior experience in genomics.

A few essential genetic loci distinguish Penstemon species with flowers adapted to pollination by bees or hummingbirds

In the formation of species, adaptation by natural selection generates distinct combinations of traits that function well together. The maintenance of adaptive trait combinations in the face of gene flow depends on the strength and nature of selection acting on the underlying genetic loci. Floral pollination syndromes exemplify the evolution of trait combinations adaptive for particular pollinators. The North American wildflower genus Penstemon displays remarkable floral syndrome convergence, with at least 20 separate lineages that have evolved from ancestral bee pollination syndrome (wide blue-purple flowers that present a landing platform for bees and small amounts of nectar) to hummingbird pollination syndrome (bright red narrowly tubular flowers offering copious nectar). Related taxa that differ in floral syndrome offer an attractive opportunity to examine the genomic basis of complex trait divergence. In this study, we characterized genomic divergence among 229 individuals from a Penstemon species complex that includes both bee and hummingbird floral syndromes. Field plants are easily classified into species based on phenotypic differences and hybrids displaying intermediate floral syndromes are rare. Despite unambiguous phenotypic differences, genome-wide differentiation between species is minimal. Hummingbird-adapted populations are more genetically similar to nearby bee-adapted populations than to geographically distant hummingbird-adapted populations, in terms of genome-wide dXY. However, a small number of genetic loci are strongly differentiated between species. These approximately 20 "species-diagnostic loci," which appear to have nearly fixed differences between pollination syndromes, are sprinkled throughout the genome in high recombination regions. Several map closely to previously established floral trait quantitative trait loci (QTLs). The striking difference between the diagnostic loci and the genome as whole suggests strong selection to maintain distinct combinations of traits, but with sufficient gene flow to homogenize the genomic background. A surprisingly small number of alleles confer phenotypic differences that form the basis of species identity in this species complex.

The East Asian wild apples, Malus baccata (L.) Borkh and Malus hupehensis (Pamp.) Rehder., are additional contributors to the genomes of cultivated European and Chinese varieties

The domestication process in long-lived plant perennials differs dramatically from that of annuals, with a huge amount of genetic exchange between crop and wild populations. Though apple is a major fruit crop grown worldwide, the contribution of wild apple species to the genetic makeup of the cultivated apple genome remains a topic of intense study. We used population genomics approaches to investigate the contributions of several wild apple species to European and Chinese rootstock and dessert genomes, with a focus on the extent of wild-crop gene flow. Population genetic structure inferences revealed that the East Asian wild apples, Malus baccata (L.) Borkh and M. hupehensis (Pamp.), form a single panmictic group, and that the European dessert and rootstock apples form a specific gene pool whereas the Chinese dessert and rootstock apples were a mixture of three wild gene pools, suggesting different evolutionary histories of European and Chinese apple varieties. Coalescent-based inferences and gene flow estimates indicated that M. baccata - M. hupehensis contributed to the genome of both European and Chinese cultivated apples through wild-to-crop introgressions, and not as an initial contributor as previously supposed. We also confirmed the contribution through wild-to-crop introgressions of Malus sylvestris Mill. to the cultivated apple genome. Apple tree domestication is therefore one example in woody perennials that involved gene flow from several wild species from multiple geographical areas. This study provides an example of a complex protracted process of domestication in long-lived plant perennials, and is a starting point for apple breeding programmes.

Parsing variance by marker type: Testing biogeographic hypotheses and differential contribution of historical processes to population structure in a desert lizard

A fundamental goal of population genetic studies is to identify historical biogeographic patterns and understand the processes that generate them. However, localized demographic events can skew population genetic inference. Assessing populations with multiple types of genetic markers, each with unique mutation rates and responses to changes in population size, can help to identify potentially confounding population-specific demographic processes. Here, we compared population structure and connectivity inferred from microsatellites and restriction site-associated DNA loci among 17 populations of an arid-specialist lizard, the desert night lizard, Xantusia vigilis, in central California to test among historical processes structuring population genetic diversity. We found that both marker types yielded generally concordant insights into population genetic structure including a major phylogenetic break maintained between two populations separated by less than 10 km, suggesting that either marker type could be used to understand generalized demographic patterns across the region for management purposes. However, we also found that the effects of demography on marker discordance could be used to elucidate population histories and distinguish among competing biogeographic hypotheses. Our results suggest that comparisons of within-population diversity across marker types provide powerful opportunities for leveraging marker discordance, particularly for understanding the creation and maintenance of contact zones among clades.

Warming appears as the main risk of non-adaptedness for western Mediterranean relict fir forests under expected climate change scenarios

Circum-Mediterranean firs are considered among the most drought-sensitive species to climate change. Understanding the genetic basis of trees' adaptive capacity and intra-specific variability to drought avoidance is mandatory to define conservation measures, thus potentially preventing their extinction. We focus here on Abies pinsapo and Abies marocana, both relict tree species, endemic from south Spain and north Morocco, respectively. A total of 607 samples were collected from eight nuclei: six from Spanish fir and two from Moroccan fir. A genotyping by sequencing technique called double digestion restriction site-associated DNA sequencing (ddRAD-seq) was performed to obtain a genetic matrix based on single-nucleotide polymorphisms (SNPs). This matrix was utilized to study the genetic structure of A. pinsapo populations and to carry out selection signature studies. In order to understand how Spanish fir and Moroccan fir cope with climate change, genotype-environment associations (GEAs) were identified. Further, the vulnerability of these species to climate variations was estimated by the risk of non-adaptedness (RONA). The filtering of the de novo assembly of A. pinsapo provided 3,982 SNPs from 504 out of 509 trees sequenced. Principal component analysis (PCA) genetically separated Grazalema from the rest of the Spanish populations. However, FST values showed significant differences among the sampling points. We found 51 loci potentially under selection. Homolog sequences were found for some proteins related to abiotic stress response, such as dehydration-responsive element binding transcription factor, regulation of abscisic acid signaling, and methylation pathway. A total of 15 associations with 11 different loci were observed in the GEA studies, with the maximum temperature of the warmest month being the variable with the highest number of associated loci. This temperature sensitivity was also supported by the risk of non-adaptedness, which yielded a higher risk for both A. pinsapo and A. marocana under the high emission scenario (Representative Concentration Pathway (RCP) 8.5). This study sheds light on the response to climate change of these two endemic species.

Dynamic evolution of size and colour in the highly specialized Zodarion ant-eating spiders

Ecological specialists constitute relevant case studies for understanding the mechanisms, potential and limitations of evolution. The species-rich and strictly myrmecophagous spiders of the genus Zodarion show diversified defence mechanisms, including myrmecomorphy of different ant species and nocturnality. Through Hybridization Capture Using RAD Probes (hyRAD), a phylogenomic technique designed for sequencing poorly preserved specimens, we reconstructed a phylogeny of Zodarion using 52 (approx. a third of the nominal) species that cover its phylogenetic and distributional diversity. We then estimated the evolution of body size and colour, traits that have diversified noticeably and are linked to defence mechanisms, across the group. Our genomic matrix of 300 loci led to a well-supported phylogenetic hypothesis that uncovered two main clades inside Zodarion. Ancestral state estimation revealed the highly dynamic evolution of body size and colour across the group, with multiple transitions and convergences in both traits, which we propose is likely indicative of multiple transitions in ant specialization across the genus. Our study will allow the informed targeted selection of Zodarion taxa of special interest for research into the group's remarkable adaptations to ant specialization. It also exemplifies the utility of hyRAD for phylogenetic studies using museum material.

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.

Complex Patterns of Diversification in the Gray Zone of Speciation: Model-Based Approaches Applied to Patagonian Liolaemid Lizards (Squamata: Liolaemus kingii clade)

In this study we detangled the evolutionary history of the Patagonian lizard clade Liolaemus kingii, coupling dense geographic sampling and novel computational analytical approaches. We analyzed nuclear and mitochondrial data (restriction site-associated DNA sequencing and cytochrome b) to hypothesize and evaluate species limits, phylogenetic relationships, and demographic histories. We complemented these analyses with posterior predictive simulations to assess the fit of the genomic data to the multispecies coalescent model. We also employed a novel approach to time-calibrate a phylogenetic network. Our results show several instances of mito-nuclear discordance and consistent support for a reticulated history, supporting the view that the complex evolutionary history of the kingii clade is characterized by extensive gene flow and rapid diversification events. We discuss our findings in the contexts of the "gray zone" of speciation, phylogeographic patterns in the Patagonian region, and taxonomic outcomes. [Model adequacy; multispecies coalescent; multispecies network coalescent; phylogenomics; species delimitation.].

The importance of cryptic diversity in the conservation of wide-ranging species: The red-footed tortoise Chelonoidis carbonarius in Colombia

Wide-ranging species are seldom considered conservation priorities, yet they have the potential to harbour genetically deeply differentiated units across environments or ecological barriers, including some that warrant taxonomic recognition. Documenting such cryptic genetic diversity is especially important for wide-ranging species that are in decline, as they may comprise a set of even more endangered lineages or species with small distributions. However, studies of wide-ranging species, particularly when they cross political borders, are extremely challenging. One approach to overcoming these challenges is to conduct detailed local analyses in combination with less detailed, range-wide studies. We used this approach with the red-footed tortoise (Chelonoidis carbonarius), a threatened species likely to contain cryptic diversity given its vast range and the distinctive ecoregions that it inhabits. Previous single-gene molecular studies indicated the presence of at least five lineages, two of which occur in different ecoregions separated by the Andes within Colombia. We used a comprehensive genomic analysis to test the hypothesis of cryptic diversity within the single jurisdiction of Colombia. We used a combination of restriction-site-associated DNA sequencing and environmental niche modelling to provide three independent lines of evidence that support the presence of important cryptic diversity that may deserve taxonomic recognition: allopatric reproductive isolation, local adaptation and ecological divergence. We also provide a fine-scale genetic map with the distribution of conservation units in Colombia. As we complete ongoing range-wide analyses and make taxonomic adjustments, we recommend that the two lineages in Colombia be treated as separate units for conservation purposes.

Genomic analysis unveils reduced genetic variability but increased proportion of heterozygotic genotypes of the intensively managed mezcal agave, Agave angustifolia

PREMISE

The central Oaxaca Basin has a century-long history of agave cultivation and is hypothesized to be the region of origin of other cultivated crops. Widely cultivated for mezcal production, the perennial crop known as "espadín" is putatively derived from wild Agave angustifolia. Nevertheless, little is known about its genetic relationship to the wild A. angustifolia or how the decades-long clonal propagation has affected its genetics.

METHODS

Using restriction-site-associated DNA sequencing and over 8000 single-nucleotide polymorphisms, we studied aspects of the population genomics of wild and cultivated A. angustifolia in Puebla and Oaxaca, Mexico. We assessed patterns of genetic diversity, inbreeding, distribution of genetic variation, and differentiation among and within wild populations and plantations.

RESULTS

Genetic differentiation between wild and cultivated plants was strong, and both gene pools harbored multiple unique alleles. Nevertheless, we found several cultivated individuals with high genetic affinity with wild samples. Higher heterozygosity was observed in the cultivated individuals, while in total, they harbored considerably fewer alleles and presented higher linkage disequilibrium compared to the wild plants. Independently of geographic distance among sampled plantations, the genetic relatedness of the cultivated plants was high, suggesting a common origin and prevalent role of clonal propagation.

CONCLUSIONS

The considerable heterozygosity found in espadín is contained within a network of highly related individuals, displaying high linkage disequilibrium generated by decades of clonal propagation and possibly by the accumulation of somatic mutations. Wild A. angustifolia, on the other hand, represents a significant genetic diversity reservoir that should be carefully studied and conserved.

Speciation rates are unrelated to the formation of population structure in Malagasy gemsnakes

Speciation rates vary substantially across the tree of life. These rates should be linked to the rate at which population structure forms if a continuum between micro and macroevolutionary patterns exists. Previous studies examining the link between speciation rates and the degree of population formation in clades have been shown to be either correlated or uncorrelated depending on the group, but no study has yet examined the relationship between speciation rates and population structure in a young group that is constrained spatially to a single-island system. We examine this correlation in 109 gemsnakes (Pseudoxyrhophiidae) endemic to Madagascar and originating in the early Miocene, which helps control for extinction variation across time and space. We find no relationship between rates of speciation and the formation rates of population structure over space in 33 species of gemsnakes. Rates of speciation show low variation, yet population structure varies widely across species, indicating that speciation rates and population structure are disconnected. We suspect this is largely due to the persistence of some lineages not susceptible to extinction. Importantly, we discuss how delimiting populations versus species may contribute to problems understanding the continuum between shallow and deep evolutionary processes.

Climate and Soil Microsite Conditions Determine Local Adaptation in Declining Silver Fir Forests

Ongoing climatic change is threatening the survival of drought-sensitive tree species, such as silver fir (Abies alba). Drought-induced dieback had been previously explored in this conifer, although the role played by tree-level genetic diversity and its relationship with growth patterns and soil microsite conditions remained elusive. We used double digest restriction-site-associated DNA sequencing (ddRADseq) to describe different genetic characteristics of five silver fir forests in the Spanish Pyrenees, including declining and non-declining trees. Single nucleotide polymorphisms (SNPs) were used to investigate the relationships between genetics, dieback, intraspecific trait variation (functional dendrophenotypic traits and leaf traits), local bioclimatic conditions, and rhizosphere soil properties. While there were no noticeable genetic differences between declining and non-declining trees, genome-environment associations with selection signatures were abundant, suggesting a strong influence of climate, soil physicochemical properties, and soil microbial diversity on local adaptation. These results provide novel insights into how genetics and diverse environmental factors are interrelated and highlight the need to incorporate genetic data into silver fir forest dieback studies to gain a better understanding of local adaptation.

Demographic responses of hybridizing cinquefoils to changing climate in the Colorado Rocky Mountains

Hybridization between taxa generates new pools of genetic variation that can lead to different environmental responses and demographic trajectories over time than seen in parental lineages. The potential for hybrids to have novel environmental tolerances may be increasingly important in mountainous regions, which are rapidly warming and drying due to climate change. Demographic analysis makes it possible to quantify within- and among-species responses to variation in climate and to predict population growth rates as those conditions change. We estimated vital rates and population growth in 13 natural populations of two cinquefoil taxa (Potentilla hippiana and P. pulcherrima) and their hybrid across elevation gradients in the Southern Rockies. Using three consecutive years of environmental and demographic data, we compared the demographic responses of hybrid and parental taxa to environmental variation across space and time. All three taxa had lower predicted population growth rates under warm, dry conditions. However, the magnitude of these responses varied among taxa and populations. Hybrids had consistently lower predicted population growth rates than P. hippiana. In contrast, hybrid performance relative to P. pulcherrima varied with population and climate, with the hybrid maintaining relatively stable growth rates while populations of P. pulcherrima shrank under warm, dry conditions. Our findings demonstrate that hybrids in this system are neither intrinsically unfit nor universally more vigorous than parents, suggesting that the demographic consequences of hybridization are context-dependent. Our results also imply that shifts to warmer and drier conditions could have particularly negative repercussions for P. pulcherrima, which is currently the most abundant taxon in the study area, possibly as a legacy of more favorable historical climates. More broadly, the distributions of these long-lived taxa are lagging behind their demographic trajectories, such that the currently less common P. hippiana could become the most abundant of the Potentilla taxa as this region continues to warm and dry.

Mitochondrial introgression and mito-nuclear discordance obscured the closely related species boundaries in Cletus Stål from China (Heteroptera: Coreidae)

Accurate taxonomy and delimitation are of great importance for pest control strategies and management programs. Here, we focus on Cletus (Insecta: Hemiptera: Coreidae), which includes many crop pests. The species boundaries still conflict and only cytochrome c oxidase subunit I (COI) barcoding has been previously used for molecular studies. We generated new mitochondrial genome and nuclear genome-wide SNPs to explore the species boundaries of 46 Cletus samples from China using multiple species delimitation approaches. All results recovered a monophyly with high support, except for two closely related species in clade I - C. punctiger and C. graminis. Mitochondrial data demonstrated admixture in clade I, while genome-wide SNPs unambiguously identified two separate species, which were confirmed by morphological classification. Inconsistent nuclear and mitochondrial data indicated mito-nuclear discordance. Mitochondrial introgression is the most likely explanation, and more extensive sampling and more comprehensive data are needed to ascertain a pattern. Accurate species delimitation will shed light on species status; thus, an accurate taxonomy is of particular concern, as there is a pressing need to implement precise control of agricultural pests and to perform further research on diversification.

Molecular phylogeny and inflorescence evolution of Prunus (Rosaceae) based on RAD-seq and genome skimming analyses

Prunus is an economically important genus widely distributed in the temperate Northern Hemisphere. Previous studies on the genus using a variety of loci yielded conflicting phylogenetic hypotheses. Here, we generated nuclear reduced representation sequencing data and plastid genomes for 36 Prunus individuals and two outgroups. Both nuclear and plastome data recovered a well-resolved phylogeny. The species were divided into three main clades corresponding to their inflorescence types, - the racemose group, the solitary-flower group and the corymbose group - with the latter two sister to one another. Prunus was inferred to have diversified initially in the Late Cretaceous around 67.32 million years ago. The diversification of the three major clades began between the Paleocene and Miocene, suggesting that paleoclimatic events were an important driving force for Prunus diversification. Ancestral state reconstructions revealed that the most recent common ancestor of Prunus had racemose inflorescences, and the solitary-flower and corymb inflorescence types were derived by reduction of flower number and suppression of the rachis, respectively. We also tested the hybrid origin hypothesis of the racemose group proposed in previous studies. Prunus has undergone extensive hybridization events, although it is difficult to identify conclusively specific instances of hybridization when using SNP data, especially deep in the phylogeny. Our study provides well-resolved nuclear and plastid phylogenies of Prunus, reveals substantial cytonuclear discord at shallow scales, and sheds new light on inflorescence evolution in this economically important lineage.

Global phylogenomic assessment of Leptoseris and Agaricia reveals substantial undescribed diversity at mesophotic depths

BACKGROUND

Mesophotic coral communities are increasingly gaining attention for the unique biological diversity they host, exemplified by the numerous mesophotic fish species that continue to be discovered. In contrast, many of the photosynthetic scleractinian corals observed at mesophotic depths are assumed to be depth-generalists, with very few species characterised as mesophotic-specialists. This presumed lack of a specialised community remains largely untested, as phylogenetic studies on corals have rarely included mesophotic samples and have long suffered from resolution issues associated with traditional sequence markers.

RESULTS

Here, we used reduced-representation genome sequencing to conduct a phylogenomic assessment of the two dominant mesophotic genera of plating corals in the Indo-Pacific and Western Atlantic, respectively, Leptoseris and Agaricia. While these genome-wide phylogenies broadly corroborated the morphological taxonomy, they also exposed deep divergences within the two genera and undescribed diversity across the current taxonomic species. Five of the eight focal species consisted of at least two sympatric and genetically distinct lineages, which were consistently detected across different methods.

CONCLUSIONS

The repeated observation of genetically divergent lineages associated with mesophotic depths highlights that there may be many more mesophotic-specialist coral species than currently acknowledged and that an urgent assessment of this largely unstudied biological diversity is warranted.

Genomic data provide a robust phylogeny backbone for Rhodiola L. (Crassulaceae) and reveal extensive reticulate evolution during its rapid radiation

The Tibetan Plateau and adjacent mountain regions (TP; including the Tibetan Plateau, Himalaya, Hengduan Mountains and Mountains of Central Asia) harbor great biodiversity, some lineages on which may have undergone rapid radiations. However, only a few studies have investigated the evolutionary pattern of such diversification in depth using genomic data. In this study, we reconstructed a robust phylogeny backbone of Rhodiola, a lineage that may have undergone rapid radiation in the TP, using Genotyping-by-sequencing data, and conducted a series of gene flow and diversification analyses. The concatenation and coalescent-based methods yield similar tree topologies, and five well-supported clades were revealed. Potential gene flow and introgression events were detected, both between species from different major clades and closely related species, suggesting pervasive hybridization and introgression. An initial rapid and later slowdown of the diversification rate was revealed, indicating niche filling. Molecular dating and correlation analyses showed that the uplift of TP and global cooling in the mid-Miocene might have played an important role in promoting the rapid radiation of Rhodiola. Our work demonstrates that gene flow and introgression might be an important contributor to rapid radiation possibly by quickly reassembling old genetic variation into new combinations.

Species limits and diversification of the Dendropsophus rubicundulus subgroup (Anura, Hylidae) in Neotropical savannas

Understanding the processes that generate and maintain biodiversity at and below the species level is a central goal of evolutionary biology. Here we explore the spatial and temporal drivers of diversification of the treefrog subgroup Dendropsophus rubicundulus, a subgroup of the D. microcephalus species group, over periods of pronounced geological and climatic changes in the Neotropical savannas that they inhabit. This subgroup currently comprises 11 recognized species distributed across the Brazilian and Bolivian savannas, but the taxonomy has been in a state of flux, necessitating reexamination. Using newly generated single nucleotide polymorphism (SNP) data from restriction-site associated DNA sequencing (RADseq) and mitochondrial 16S sequence data for ∼150 specimens, we inferred phylogenetic relationships, tested species limits using a model-based approach, and estimated divergence times to gain insights into the geographic and climatic events that affected the diversification of this subgroup. Our results recognized at least nine species: D. anataliasiasi, D. araguaya, D. cerradensis, D. elianeae, D. jimi, D. rubicundulus, D. tritaeniatus, D. rozenmani, and D. sanborni. Although we did not collect SNP data for the latter two species, they are likely distinct based on mitochondrial data. In addition, we found genetic structure within the widespread species D. rubicundulus, which comprises three allopatric lineages connected by gene flow upon secondary contact. We also found evidence of population structure and perhaps undescribed diversity in D. elianeae, which warrants further study. The D. rubicundulus subgroup is estimated to have originated in the Late Miocene (∼5.45 million years ago), with diversification continuing through the Pliocene and Early Pleistocene, followed by the most recent divergence of D. rubicundulus lineages in the Middle Pleistocene. The epeirogenic uplift followed by erosion and denudation of the central Brazilian plateau throughout the Pliocene and Pleistocene, in combination with the increasing frequency and amplitude of climatic fluctuations during the Pleistocene, was important for generating and structuring diversity at or below the species level in the D. rubicundulus subgroup.

A Foundational Population Genetics Investigation of the Sexual Systems of Solanum (Solanaceae) in the Australian Monsoon Tropics Suggests Dioecious Taxa May Benefit from Increased Genetic Admixture via Obligate Outcrossing

Solanum section Leptostemonum is an ideal lineage to test the theoretical framework regarding proposed evolutionary benefits of outcrossing sexual systems in comparison to cosexuality. Theoretically, non-cosexual taxa should support more genetic diversity within populations, experience less inbreeding, and have less genetic structure due to a restricted ability to self-fertilize. However, many confounding factors present challenges for a confident inference that inherent differences in sexual systems influence observed genetic patterns among populations. This study provides a foundational baseline of the population genetics of several species of different sexual systems with the aim of generating hypotheses of any factor-including sexual system-that influences genetic patterns. Importantly, results indicate that dioecious S. asymmetriphyllum maintains less genetic structure and greater admixture among populations than cosexual S. raphiotes at the same three locations where they co-occur. This suggests that when certain conditions are met, the evolution of dioecy may have proceeded as a means to avoid genetic consequences of self-compatibility and may support hypotheses of benefits gained through differential resource allocation partitioned across sexes. Arguably, the most significant finding of this study is that all taxa are strongly inbred, possibly reflective of a shared response to recent climate shifts, such as the increased frequency and intensity of the region's fire regime.