RADseq provides evidence for parallel ecotypic divergence in the autotetraploid Cochlearia officinalis in Northern Norway

Genetic diversity and population structure of Piper nigrum (black pepper) accessions based on next-generation SNP markers

Despite the economic importance of Piper nigrum (black pepper), a highly valued crop worldwide, development and utilization of genomic resources have remained limited, with diversity assessments often relying on only a few samples or DNA markers. Here we employed restriction-site associated DNA sequencing to analyze 175 P. nigrum accessions from eight main black pepper growing regions in Sri Lanka. The sequencing effort resulted in 1,976 million raw reads, averaging 11.3 million reads per accession, revealing 150,356 high-quality single nucleotide polymorphisms (SNPs) distributed across 26 chromosomes. Population structure analysis revealed two subpopulations (K = 2): a dominant group consisting of 152 accessions sourced from both home gardens and large-scale cultivations, and a smaller group comprising 23 accessions exclusively from native collections in home gardens. This clustering was further supported by principal component analysis, with the first two principal components explaining 35.2 and 12.1% of the total variation. Genetic diversity analysis indicated substantial gene flow (Nm = 342.21) and a low fixation index (FST = 0.00073) between the two subpopulations, with no clear genetic differentiation among accessions from different agro-climatic regions. These findings demonstrate that most current black pepper genotypes grown in Sri Lanka share a common genetic background, emphasizing the necessity to broaden the genetic base to enhance resilience to biotic and abiotic stresses. This study represents the first attempt at analyzing black pepper genetic diversity using high-resolution SNP markers, laying the foundation for future genome-wide association studies for SNP-based gene discovery and breeding.

Impact of whole-genome duplications on structural variant evolution in Cochlearia

Polyploidy, the result of whole-genome duplication (WGD), is a major driver of eukaryote evolution. Yet WGDs are hugely disruptive mutations, and we still lack a clear understanding of their fitness consequences. Here, we study whether WGDs result in greater diversity of genomic structural variants (SVs) and how they influence evolutionary dynamics in a plant genus, Cochlearia (Brassicaceae). By using long-read sequencing and a graph-based pangenome, we find both negative and positive interactions between WGDs and SVs. Masking of recessive mutations due to WGDs leads to a progressive accumulation of deleterious SVs across four ploidal levels (from diploids to octoploids), likely reducing the adaptive potential of polyploid populations. However, we also discover putative benefits arising from SV accumulation, as more ploidy-specific SVs harbor signals of local adaptation in polyploids than in diploids. Together, our results suggest that SVs play diverse and contrasting roles in the evolutionary trajectories of young polyploids.

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

BACKGROUND AND AIMS

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

METHODS

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

KEY RESULTS

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

CONCLUSIONS

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

Homoplastic versus xenoplastic evolution: exploring the emergence of key intrinsic and extrinsic traits in the montane genus Soldanella (Primulaceae)

Specific ecological conditions in the high mountain environment exert a selective pressure that often leads to convergent trait evolution. Reticulations induced by incomplete lineage sorting and introgression can lead to discordant trait patterns among gene and species trees (hemiplasy/xenoplasy), providing a false illusion that the traits under study are homoplastic. Using phylogenetic species networks, we explored the effect of gene exchange on trait evolution in Soldanella, a genus profoundly influenced by historical introgression. At least three features evolved independently multiple times: the single-flowered dwarf phenotype, dysploid cytotype, and ecological generalism. The present analyses also indicated that the recurring occurrence of stoloniferous growth might have been prompted by an introgression event between an ancestral lineage and a still extant species, although its emergence via convergent evolution cannot be completely ruled out. Phylogenetic regression suggested that the independent evolution of larger genomes in snowbells is most likely a result of the interplay between hybridization events of dysploid and euploid taxa and hostile environments at the range margins of the genus. The emergence of key intrinsic and extrinsic traits in snowbells has been significantly impacted not only by convergent evolution but also by historical and recent introgression events.

Towards the Investigation of the Adaptive Divergence in a Species of Exceptional Ecological Plasticity: Chromosome-Scale Genome Assembly of Chouardia litardierei (Hyacinthaceae)

One of the central goals of evolutionary biology is to understand the genomic basis of adaptive divergence. Different aspects of evolutionary processes should be studied through genome-wide approaches, therefore maximizing the investigated genomic space. However, in-depth genome-scale analyses often are restricted to a model or economically important species and their closely related wild congeners with available reference genomes. Here, we present the high-quality chromosome-level genome assembly of Chouardia litardierei, a plant species with exceptional ecological plasticity. By combining PacBio and Hi-C sequencing technologies, we generated a 3.7 Gbp genome with a scaffold N50 size of 210 Mbp. Over 80% of the genome comprised repetitive elements, among which the LTR retrotransposons prevailed. Approximately 86% of the 27,257 predicted genes were functionally annotated using public databases. For the comparative analysis of different ecotypes' genomes, the whole-genome sequencing of two individuals, each from a distinct ecotype, was performed. The detected above-average SNP density within coding regions suggests increased adaptive divergence-related mutation rates, therefore confirming the assumed divergence processes within the group. The constructed genome presents an invaluable resource for future research activities oriented toward the investigation of the genetics underlying the adaptive divergence that is likely unfolding among the studied species' ecotypes.

Pervasive Introgression During Rapid Diversification of the European Mountain Genus Soldanella (L.) (Primulaceae)

Hybridization is a key mechanism involved in lineage diversification and speciation, especially in ecosystems that experienced repeated environmental oscillations. Recently radiated plant groups, which have evolved in mountain ecosystems impacted by historical climate change provide an excellent model system for studying the impact of gene flow on speciation. We combined organellar (whole-plastome) and nuclear genomic data (RAD-seq) with a cytogenetic approach (rDNA FISH) to investigate the effects of hybridization and introgression on evolution and speciation in the genus Soldanella (snowbells, Primulaceae). Pervasive introgression has already occurred among ancestral lineages of snowbells and has persisted throughout the entire evolutionary history of the genus, regardless of the ecology, cytotype, or distribution range size of the affected species. The highest extent of introgression has been detected in the Carpathian species, which is also reflected in their extensive karyotype variation. Introgression occurred even between species with dysploid and euploid cytotypes, which were considered to be reproductively isolated. The magnitude of introgression detected in snowbells is unprecedented in other mountain genera of the European Alpine System investigated hitherto. Our study stresses the prominent evolutionary role of hybridization in facilitating speciation and diversification on the one hand, but also enriching previously isolated genetic pools. [chloroplast capture; diversification; dysploidy; European Alpine system; introgression; nuclear-cytoplasmic discordance; ribosomal DNA.].

Genetic structure and geneflow of Malus across the Korean Peninsula using genotyping-by-sequencing

This study was to understand the genetic structure and diversity of the Korean Malus species. We used genotyping-by-sequencing (GBS) technology to analyze samples of 112 individuals belonging to 18 populations of wild Malus spp. Using GBS, we identified thousands of single nucleotide polymorphisms in the species analyzed. M. baccata and M. toringo, two dominant mainland species of the Korean Peninsula, were distinguishable based on their genetic structure. However, M. toringo collected from Jeju Island exhibited a different genetic profile than that from the mainland. We identified M. cf. micromalus as a hybrid resulting from the Jeju Island M. toringo (pollen donor) and the mainland M. baccata, (pollen recipient). Putative M. mandshurica distributed on the Korean Peninsula showed a high structural and genetic similarity with M. baccata, indicating that it might be an ecotype. Overall, this study contributes to the understanding of the population history and genetic structure of Malus in the Korean Peninsula.

Little hope for the polyploid endemic Pyrenean Larkspur ( Delphinium montanum ): Evidences from population genomics and Ecological Niche Modeling

Species endemic to restricted geographical ranges represent a particular conservation issue, be it for their heritage interest. In a context of global change, this is particularly the case for plants which belong to high‐mountain ecosystems and, because of their ecological requirements, are doomed to survive or disappear on their “sky islands”. The Pyrenean Larkspur (Delphinium montanum, Ranunculaceae) is endemic to the Eastern part of the Pyrenees (France and Spain). It is now only observable at a dozen of localities and some populations show signs of decline, such as a recurrent lack of flowering. Implementing population genomics approach (e.g., RAD‐seq like) is particularly useful to understand genomic patterns of diversity and differentiation in order to provide recommendations in term of conservation. However, it remains challenging for species such as D. montanum that are autotetraploid with a large genome size (1C‐value >10 pg) as most methods currently available were developed for diploid species. A Bayesian framework able to call genotypes with uncertainty allowed us to assess genetic diversity and population structure in this system. Our results show evidence for inbreeding (mean G_IS = 0.361) within all the populations and substantial population structure (mean G_ST = 0.403) at the metapopulation level. In addition to a lack of connectivity between populations, spatial projections of Ecological Niche Modeling (ENM) analyses under different climatic scenarios predict a dramatic decrease of suitable habitat for D. montanum in the future. Based on these results, we discuss the relevance and feasibility of different conservation measures.

A phylogeny of Antirrhinum reveals parallel evolution of alpine morphology

Parallel evolution of similar morphologies in closely related lineages provides insight into the repeatability and predictability of evolution. In the genus Antirrhinum (snapdragons), as in other plants, a suite of morphological characters are associated with adaptation to alpine environments. We tested for parallel trait evolution in Antirrhinum by investigating phylogenetic relationships using restriction-site associated DNA (RAD) sequencing. We then associated phenotypic information to our phylogeny to reconstruct the patterns of morphological evolution and related this to evidence for hybridisation between emergent lineages. Phylogenetic analyses showed that the alpine character syndrome is present in multiple groups, suggesting that Antirrhinum has repeatedly colonised alpine habitats. Dispersal to novel environments happened in the presence of intraspecific and interspecific gene flow. We found support for a model of parallel evolution in Antirrhinum. Hybridisation in natural populations, and a complex genetic architecture underlying the alpine morphology syndrome, support an important role of natural selection in maintaining species divergence in the face of gene flow.

Evolutionary footprints of a cold relic in a rapidly warming world

With accelerating global warming, understanding the evolutionary dynamics of plant adaptation to environmental change is increasingly urgent. Here, we reveal the enigmatic history of the genus Cochlearia (Brassicaceae), a Pleistocene relic that originated from a drought-adapted Mediterranean sister genus during the Miocene. Cochlearia rapidly diversified and adapted to circum-Arctic regions and other cold-characterized habitat types during the Pleistocene. This sudden change in ecological preferences was accompanied by a highly complex, reticulate polyploid evolution, which was apparently triggered by the impact of repeated Pleistocene glaciation cycles. Our results illustrate that two early diversified Arctic-alpine diploid gene pools contributed differently to the evolution of this young polyploid genus now captured in a cold-adapted niche. Metabolomics revealed central carbon metabolism responses to cold in diverse species and ecotypes, likely due to continuous connections to cold habitats that may have facilitated widespread adaptation to alpine and subalpine habitats, and which we speculate were coopted from existing drought adaptations. Given the growing scientific interest in the adaptive evolution of temperature-related traits, our results provide much-needed taxonomic and phylogenomic resolution of a model system as well as first insights into the origins of its adaptation to cold.

Species boundaries in the messy middle-A genome-scale validation of species delimitation in a recently diverged lineage of coastal fog desert lichen fungi

Species delimitation among closely related species is challenging because traditional phenotype-based approaches, for example, using morphology, ecological, or chemical characteristics, may not coincide with natural groupings. With the advent of high-throughput sequencing, it has become increasingly cost-effective to acquire genome-scale data which can resolve previously ambiguous species boundaries. As the availability of genome-scale data has increased, numerous species delimitation analyses, such as BPP and SNAPP+Bayes factor delimitation (BFD*), have been developed to delimit species boundaries. However, even empirical molecular species delimitation approaches can be biased by confounding evolutionary factors, for example, hybridization/introgression and incomplete lineage sorting, and computational limitations. Here, we investigate species boundaries and the potential for micro-endemism in a lineage of lichen-forming fungi, Niebla Rundel & Bowler, in the family Ramalinaceae by analyzing single-locus and genome-scale data consisting of (a) single-locus species delimitation analysis using ASAP, (b) maximum likelihood-based phylogenetic tree inference, (c) genome-scale species delimitation models, e.g., BPP and SNAPP+BFD, and (d) species validation using the genealogical divergence index (gdi). We specifically use these methods to cross-validate results between genome-scale and single-locus datasets, differently sampled subsets of genomic data and to control for population-level genetic divergence. Our species delimitation models tend to support more speciose groupings that were inconsistent with traditional taxonomy, supporting a hypothesis of micro-endemism, which may include morphologically cryptic species. However, the models did not converge on robust, consistent species delimitations. While the results of our analysis are somewhat ambiguous in terms of species boundaries, they provide a valuable perspective on how to use these empirical species delimitation methods in a nonmodel system. This study thus highlights the challenges inherent in delimiting species, particularly in groups such as Niebla, with complex, relatively recent phylogeographic histories.

Spatial and Ecological Drivers of Genetic Structure in Greek Populations of Alkanna tinctoria (Boraginaceae), a Polyploid Medicinal Herb

Background and Aims: Quantifying genetic variation is fundamental to understand a species' demographic trajectory and its ability to adapt to future changes. In comparison with diploids, however, genetic variation and factors fostering genetic divergence remain poorly studied in polyploids due to analytical challenges. Here, by employing a ploidy-aware framework, we investigated the genetic structure and its determinants in polyploid Alkanna tinctoria (Boraginaceae), an ancient medicinal herb that is the source of bioactive compounds known as alkannin and shikonin (A/S). From a practical perspective, such investigation can inform biodiversity management strategies. Methods: We collected 14 populations of A. tinctoria within its main distribution range in Greece and genotyped them using restriction site-associated DNA sequencing. In addition, we included two populations of A. sieberi. By using a ploidy-aware genotype calling based on likelihoods, we generated a dataset of 16,107 high-quality SNPs. Classical and model-based analysis was done to characterize the genetic structure within and between the sampled populations, complemented by genome size measurements and chromosomal counts. Finally, to reveal the drivers of genetic structure, we searched for associations between allele frequencies and spatial and climatic variables. Key Results: We found support for a marked regional structure in A. tinctoria along a latitudinal gradient in line with phytogeographic divisions. Several analyses identified interspecific admixture affecting both mainland and island populations. Modeling of spatial and climatic variables further demonstrated a larger contribution of neutral processes and a lesser albeit significant role of selection in shaping the observed genetic structure in A. tinctoria. Conclusion: Current findings provide evidence of strong genetic structure in A. tinctoria mainly driven by neutral processes. The revealed natural genomic variation in Greek Alkanna can be used to further predict variation in A/S production, whereas our bioinformatics approach should prove useful for the study of other non-model polyploid species.

Parallel ddRAD and Genome Skimming Analyses Reveal a Radiative and Reticulate Evolutionary History of the Temperate Bamboos

Rapid evolutionary radiations are among the most challenging phylogenetic problems, wherein different types of data (e.g., morphology and molecular) or genetic markers (e.g., nuclear and organelle) often yield inconsistent results. The tribe Arundinarieae, that is, the temperate bamboos, is a clade of tetraploid originated 22 Ma and subsequently radiated in East Asia. Previous studies of Arundinarieae have found conflicting relationships and/or low support. Here, we obtain nuclear markers from ddRAD data for 213 Arundinarieae taxa and parallel sampling of chloroplast genomes from genome skimming for 147 taxa. We first assess the feasibility of using ddRAD-seq data for phylogenetic estimates of paleopolyploid and rapidly radiated lineages, optimize clustering thresholds, and analysis workflow for orthology identification. Reference-based ddRAD data assembly approaches perform well and yield strongly supported relationships that are generally concordant with morphology-based taxonomy. We recover five major lineages, two of which are notable (the pachymorph and leptomorph lineages), in that they correspond with distinct rhizome morphologies. By contrast, the phylogeny from chloroplast genomes differed significantly. Based on multiple lines of evidence, the ddRAD tree is favored as the best species tree estimation for temperate bamboos. Using a time-calibrated ddRAD tree, we find that Arundinarieae diversified rapidly around the mid-Miocene corresponding with intensification of the East Asian monsoon and the evolution of key innovations including the leptomorph rhizomes. Our results provide a highly resolved phylogeny of Arundinarieae, shed new light on the radiation and reticulate evolutionary history of this tribe, and provide an empirical example for the study of recalcitrant plant radiations. [Arundinarieae; ddRAD; paleopolyploid; genome skimming; rapid diversification; incongruence.]

Phylogenomic Relationships and Evolution of Polyploid Salix Species Revealed by RAD Sequencing Data

Polyploidy is common in the genus Salix. However, little is known about the origin, parentage and genomic composition of polyploid species because of a lack of suitable molecular markers and analysis tools. We established a phylogenomic framework including species of all described sections of Eurasian shrub willows. We analyzed the genomic composition of seven polyploid willow species in comparison to putative diploid parental species to draw conclusions on their origin and the effects of backcrossing and post-origin evolution. We applied recently developed programs like SNAPP, HyDe, and SNiPloid to establish a bioinformatic pipeline for unravelling the complexity of polyploid genomes. RAD sequencing revealed 23,393 loci and 320,010 high quality SNPs for the analysis of relationships of 35 species of Eurasian shrub willows (Salix subg. Chamaetia/Vetrix). Polyploid willow species appear to be predominantly of allopolyploid origin. More ancient allopolyploidization events were observed for two hexaploid and one octoploid species, while our data suggested a more recent allopolyploid origin for the included tetraploids and identified putative parental taxa. SNiPloid analyses disentangled the different genomic signatures resulting from hybrid origin, backcrossing, and secondary post-origin evolution in the polyploid species. Our RAD sequencing data demonstrate that willow genomes are shaped by ancient and recent reticulate evolution, polyploidization, and post-origin divergence of species.

The conundrum of species delimitation: a genomic perspective on a mitogenetically super-variable butterfly

The Palaearctic butterfly Melitaea didyma stands out as one of the most striking cases of intraspecific genetic differentiation detected in Lepidoptera: 11 partially sympatric mitochondrial lineages have been reported, displaying levels of divergence of up to 7.4%. To better understand the evolutionary processes underlying the diversity observed in mtDNA, we compared mtDNA and genome-wide SNP data using double-digest restriction site-associated DNA sequencing (ddRADseq) results from 93 specimens of M. didyma ranging from Morocco to eastern Kazakhstan. We found that, between ddRADseq and mtDNA results, there is a match only in populations that probably remained allopatric for long periods of time. Other mtDNA lineages may have resulted from introgression events and were probably affected by Wolbachia infection. The five main ddRADseq clades supported by STRUCTURE were parapatric or allopatric and showed high pairwise F_ST values, but some were also estimated to display various levels of gene flow. Melitaea didyma represents one of the first cases of deep mtDNA splits among European butterflies assessed by a genome-wide DNA analysis and reveals that the interpretation of patterns remains challenging even when a high amount of genomic data is available. These findings actualize the ongoing debate of species delimitation in allopatry, an issue probably of relevance to a significant proportion of global biodiversity.

Multiple auto- and allopolyploidisations marked the Pleistocene history of the widespread Eurasian steppe plant Astragalus onobrychis (Fabaceae)

The Eurasian steppes occupy a significant portion of the worldwide land surface and their biota have been affected by specific past range dynamics driven by ice ages-related climatic fluctuations. The dynamic alterations in conditions during the Pleistocene often triggered reticulate evolution and whole genome duplication events. Employing genomic, genetic and cytogenetic tools as well as morphometry we investigate the intricate evolution of Astragalus onobrychis, a widespread Eurasian steppe plant with diploid, tetraploid and octoploid cytotypes. To analyse the heteroploid RADseq dataset we employ both genotype-based and genotype-free methods that result in highly consistent results, and complement our inference with information from the plastid ycf1 region. We uncover a complex and reticulate evolutionary history, including at least one auto-tetraploidization event and two allo-octoploidization events; one of them involved also genetic contributions from other species, most likely A. goktschaicus. The present genetic structure points to the existence of four main clades within A. onobrychis, which only partly correspond to different ploidies. Time-calibrated diffusion models suggest that diversification within A. onobrychis was associated with ice age-related climatic fluctuations during the last million years. We finally argue for the usefulness of uniparentally inherited plastid markers, even in the genomic era, especially when investigating heteroploid systems.

Phylogenomics resolves evolutionary relationships and provides insights into floral evolution in the tribe Shoreeae (Dipterocarpaceae)

A supra-annual, community-level synchronous flowering prevails in several parts of the tropical forests of Southeast Asia and its evolution has been hypothesized to be linked to pollinator shifts. The aseasonal Southeast Asian lowland rainforests are dominated by Dipterocarpaceae, which exhibit great floral diversity, a range of pollination syndromes and include species with annual and supra-annual gregarious flowering. Phylogenetic relationships within this family are still unclear, especially in the tribe Shoreeae. Here, we develop a pipeline to maximize recovery of genome-wide SNPs from restriction-site associated DNA sequencing (RADseq) in non-model organisms across wide phylogenetic scales. We then infer phylogenomic relationships in the tribe Shoreeae using both traditional and coalescent analyses. The phylogenetic trees obtained with these methods are congruent to each other and highly resolved. They allow reconstructing the evolutionary patterns of floral traits (number of stamens, anther structure and anther/appendage size) in the group. Our inferences indicate that species with many stamens, but smaller, globose anthers and longer appendages and have evolved multiple times from species with fewer stamens, but larger, oblong anthers and shorter appendages. This could have happened in parallel to iterative shifts in pollinators across the uncovered phylogeny from larger, longer generation to smaller, shorter-generation insects that can quickly build up the necessary population sizes during mass flowering episodes.