Evolution reflects a balance between innovation and constraint, often repurposing existing components in new contexts. Convergent evolution exemplifies this interplay, with similar traits evolving independently in different species, yet the genomic mechanisms enabling this repeatability remain poorly understood. Here, by analyzing 10 chromosome-scale genome assemblies, including seven newly generated, we found that the S-locus supergene (a cluster of tightly linked genes controlling a floral dimorphism called distyly) arose independently multiple times within the primrose family, Darwin's iconic system for studying distyly. In each case, the same gene was independently duplicated and co-opted. However, the resulting genomic architectures differed, ranging from hemizygous (present on one chromosome copy) to heterozygous (on both copies), challenging the prevailing view that hemizygosity is intrinsic to S-loci and suggesting alternative evolutionary routes to distyly supergene formation. By uncovering multiple mechanisms for supergene origins, our work shows how convergent evolution can produce similar phenotypes by reusing the same genetic building blocks while exploring distinct genomic configurations.

Abstract More than 200 years ago, Alexander von Humboldt described a tree of the genus Clusia for its ability to perform crassulacean acid metabolism (CAM). This drought-adaptive metabolism allows plants to maintain photosynthesis under water limitation by temporally separating CO₂ uptake and fixation. The diversity of CAM physiotypes has fueled a debate about evolutionary constraints and the feasibility of engineering CAM into C₃ crops. The genus Clusia displays an exceptional diversity of photosynthetic physiotypes, yet genome sequences and genomic mechanisms generating this diversity remain unresolved. Here, we sequence and compare the genomes of three Clusia species spanning weak, inducible, and strong CAM. We show that polyploidization followed by transposon-mediated genic diploidization could have shaped CAM-related gene families, particularly those controlling phosphoenol-pyruvate recycling via phosphorolytic leaf starch metabolism. Our results indicate that whole-genome duplication coupled to diploidization might have driven diversification of CAM physiotypes in Clusia, providing a genomic framework for understanding CAM diversity and evolution.

The enormous wealth of data on the incidence of polyploidy and its significance for genome evolution in angiosperms has greatly advanced our understanding of plant evolution. This has enabled the identification of correlations between polyploidization and life history trait evolution as well as unraveling the complex interplay between genome multiplication and the biology of individual lineages giving rise to intricate and unique patterns of genome evolution despite shared evolutionary processes. Advances in sequencing and analytical methods result in the quickly increasing representation of other land plant lineages, such as bryophytes or ferns. These data continue to reveal that the levels of genome evolution complexity known from angiosperms are to be found across all land plants, albeit to lineage-specific extents. Greatly expanded taxonomic sampling across land plants will not only deepen our understanding of the interplay between polyploidy and the evolution of traits and life histories, e.g. in the context of the conquest of land, but also will allow an evolutionarily broader assessment of the role of polyploidy in shaping genome dynamics.

Piper is the largest genus of Piperaceae and second largest in the angiosperms, but the biogeography and timing of its global radiation has not yet been fully investigated. Here, we carried out phylogenetic analyses with a molecular clock to provide a framework to examine divergence times and global biogeography and dispersal patterns of the genus, with a focus on Palaeotropical groups. Target enrichment produced 269 genes for 122 samples of Piperaceae. For Piper, we included species representing all geographic regions, most subgenera/sections, and all clades previously identified in phylogenetic studies. We used maximum likelihood (RAxML, IQ-Tree) and coalescent (ASTRAL-III) analyses to produce phylogenetic trees for Piperaceae. We optimized branch lengths for the coalescent species tree and provided divergence time estimates using treePL calibrated with a Cretaceous fossil of P. margaritae Martínez-A. and a maximum divergence age of 90.8 Myr, derived from a previous estimate of Piperaceae crown age. BioGeoBEARS was used for the biogeographic analyses. Target enrichment produced a generally well-supported maximum likelihood tree in which there was good support for (Verhuellia ((Manekia + Zippelia) + (Piper + Peperomia))). Intercontinental dispersals of Piper from South America started during the Middle Eocene climatic optimum (c. 40 Ma), first to southern Africa and the Pacific islands, followed by New Guinea, Southeast Asia, and the Philippines. In Asia, Piper dispersed and diversified across the Wallace and Lydekker’s lines multiple times. Despite there being only four endemic species, mainland Africa was colonized twice, and it appears the genus never diversified there.

1.AbstractThe subgenusTillandsia(Bromeliaceae) belongs to one of the fastest radiating clades in the plant kingdom and is characterised by the repeated evolution of Crassulacean Acid Metabolism (CAM). Despite its complex genetic basis, this water-conserving trait has evolved independently across many plant families and is regarded as a key innovation trait and driver of ecological diversification in Bromeliaceae. By producing high-quality genome assemblies of aTillandsiaspecies pair displaying divergent photosynthetic phenotypes, and combining genome-wide investigations of synteny, TE dynamics, sequence evolution, gene family evolution and temporal differential expression, we were able to pinpoint the genomic drivers of CAM evolution inTillandsia. Several large-scale rearrangements associated with karyotype changes between the two genomes and a highly dynamic TE landscape shaped the genomes ofTillandsia. However, our analyses show that rewiring of photosynthetic metabolism is mainly obtained through regulatory evolution rather than coding sequence evolution, as CAM-related genes are differentially expressed across a 24-hour cycle between the two species, but are no candidates of positive selection. Gene orthology analyses reveal that CAM-related gene families manifesting differential expression underwent accelerated gene family expansion in the constitutive CAM species, further supporting the view of gene family evolution as a driver of CAM evolution.

Effective protection of endangered species is often limited by taxonomic discrepancies across state borders. This is also the case of the Dactylorhiza maculata agg. in Central Europe, where one to three species and several infraspecific taxa are recognized in various countries. Based on an extensive analysis of morphological variation, ploidy levels, environmental traits and habitats of 64 populations in Central Europe and adjacent regions, we aimed to propose a unified taxonomic concept applicable throughout the study area. Multivariate analysis of morphological traits revealed continuous variation at the individual level and only minor differences between particular clusters of populations. Four DNA-ploidy levels were detected using flow cytometry. Diploids (2 n = 40) and tetraploids (2 n = 80) were the most abundant and usually formed single-cytotype populations whereas DNA-triploids and DNA-hexaploids occurred only sporadically as minority cytotypes. The inferred patterns of morphological and ploidy variation were not congruent with traditional taxonomic treatment regarding diploid D. fuchsii and tetraploid D. maculata as two species with several infraspecific taxa. Instead, all taxa analysed in the current study are best treated at the subspecies level within D. maculata s. lat. due to somewhat continuous morphological variation between morphotypes. A total of eight D. maculata subspecies may be recognized in Central Europe, of which one is newly described here as D. maculata subsp. arcana , subsp. nov. Some nomenclatural riddles have been resolved, and the threat status of the recognized taxa is discussed.

Patterns of genetic variation in crops are the result of multiple processes that have occurred during their domestication and improvement, and are influenced by their wild progenitors that often remain understudied. The locoto chile, Capsicum pubescens, is a crop grown mainly in mid-highlands of South-Central America. This species is not known from the wild and exists only as a cultigen. The evolutionary affinities and exact origin of C. pubescens have still not been elucidated, with hypotheses suggesting its genetic relatedness and origin to two wild putative ancestral Capsicum species from the Central Andes, C. eximium and C. cardenasii. In the current study, RAD-sequencing was applied to obtain genome-wide data for 48 individuals of C. pubescens and its wild allies representing different geographical areas. Bayesian, Maximum Likelihood and coalescent-based analytical approaches were used to reconstruct population genetic patterns and phylogenetic relationships of the studied species. The results revealed that C. pubescens forms a well-defined monotypic lineage closely related to wild C. cardenasii and C. eximium, and also to C. eshbaughii. The primary lineages associated with the diversification under domestication of C. pubescens were also identified. Although direct ancestor-descendant relationship could not be inferred within this group of taxa, hybridization events were detected between C. pubescens and both C. cardenasii and C. eximium. Therefore, although hybrid origin of C. pubescens could not be inferred, gene flow involving its wild siblings was shown to be an important factor contributing to its contemporary genetic diversity. The data allowed for the inference of the center of origin of C. pubescens in central-western Bolivia highlands and for better understanding of the dynamics of its gene pool. The results of this study are essential for germplasm conservation and breeding purposes, and provide excellent basis for further research of the locoto chile and its wild relatives.

Chromosome number and genome size changes via dysploidy and polyploidy accompany plant diversification and speciation. Such changes often impact also morphological characters. An excellent system to address the questions of how extensive and structured chromosomal changes within one species complex affect the phenotype is the monocot species complex of Barnardia japonica. This taxon contains two well established and distinct diploid cytotypes differing in base chromosome numbers (AA: x = 8, BB: x = 9) and their allopolyploid derivatives on several ploidy levels (from 3x to 6x). This extensive and structured genomic variation, however, is not mirrored by gross morphological differentiation. The current study aims to analyze the correlations between the changes of chromosome numbers and genome sizes with palynological and leaf micromorphological characters in diploids and selected allopolyploids of the B. japonica complex. The chromosome numbers varied from 2n = 16 and 18 (2n = 25 with the presence of supernumerary B chromosomes), and from 2n = 26 to 51 in polyploids on four different ploidy levels (3x, 4x, 5x, and 6x). Despite additive chromosome numbers compared to diploid parental cytotypes, all polyploid cytotypes have experienced genome downsizing. Analyses of leaf micromorphological characters did not reveal any diagnostic traits that could be specifically assigned to individual cytotypes. The variation of pollen grain sizes correlated positively with ploidy levels. This study clearly demonstrates that karyotype and genome size differentiation does not have to be correlated with morphological differentiation of cytotypes.

Fluorescence in situ hybridization allows for the mapping of various sequence types in the genomes and is thus widely used in structural, functional, and evolutionary studies. One particular type of in situ hybridization that specifically allows to map whole parental genomes in diploid and polyploid hybrids is genomic in situ hybridization (GISH). The efficiency of GISH, i.e., the specificity of hybridization of genomic DNA probes to the parental subgenomes in hybrids depends, among others, on the age of the polyploids and the similarity of the parental genomes, specifically their repetitive DNA fractions. Typically, high levels of overall repeat similarity between the parental genomes result in lower efficiency of GISH. Here, we present the formamide-free GISH (ff-GISH) protocol that can be applied to diploid and polyploid hybrids of both monocots and dicots. ff-GISH allows higher efficiency of the labeling of the putative parental genomes compared to the standard GISH protocol and allows discrimination of parental chromosome sets that share up to 80-90% repeat similarity. This modified method is nontoxic, is simple, and lends itself to modifications. It can also be used for standard FISH and mapping of individual sequence types in chromosomes/genomes.

Although Crepis was the first model plant group in which chromosomal changes were considered to play an important role in speciation, their chromosome structure and evolution have been barely investigated using molecular cytogenetic methods. The aim of the study was to provide a better understanding of the patterns and directions of Crepis chromosome evolution, using comparative analyses of rDNA loci number and localisation. The chromosome base number and chromosomal organisation of 5S and 35S rDNA loci were analysed in the phylogenetic background for 39 species of Crepis, which represent the evolutionary lineages of Crepis sensu stricto and Lagoseris, including Lapsana communis. The phylogenetic relationships among all the species were inferred from nrITS and newly obtained 5S rDNA NTS sequences. Despite high variations in rDNA loci chromosomal organisation, most species had a chromosome with both rDNA loci within the same (usually short) chromosomal arm. The comparative analyses revealed several independent rDNA loci number gains and loci repositioning that accompanied diversification and speciation in Crepis. Some of the changes in rDNA loci patterns were reconstructed for the same evolutionary lineages as descending dysploidy.

The Iris series Chinenses in Korea comprises four species (I. minutoaurea, I. odaesanensis, I. koreana, and I. rossii), and the group includes some endangered species, owing to their high ornamental, economic, and conservation values. Among them, the putative allotetraploid, Iris koreana (2n = 4x = 50), is hypothesized to have originated from the hybridization of the diploids I. minutoaurea (2n = 2x = 22) and I. odaesanensis (2n = 2x = 28) based on morphological characters, chromosome numbers, and genome size additivity. Despite extensive morphological and molecular phylogenetical studies on the genus Iris, little is known about Korean irises in terms of their complete chloroplast (cp) genomes and molecular cytogenetics that involve rDNA loci evolution based on fluorescence in situ hybridization (FISH). This study reports comparative analyses of the karyotypes of the three Iris species (I. koreana, I. odaesanensis, and I. minutoaurea), with an emphasis on the 5S and 35S rDNA loci number and localization using FISH together with the genome size and chromosome number. Moreover, the cp genomes of the same individuals were sequenced and assembled for comparative analysis. The rDNA loci numbers, which were localized consistently at the same position in all species, and the chromosome numbers and genome size values of tetraploid Iris koreana (four 5S and 35S loci; 2n = 50; 1C = 7.35 pg) were additively compared to its putative diploid progenitors, I. minutoaurea (two 5S and 35S loci; 2n = 22; 1C = 3.71 pg) and I. odaesanensis (two 5S and 35S loci; 2n = 28; 1C = 3.68 pg). The chloroplast genomes were 152,259–155,145 bp in length, and exhibited a conserved quadripartite structure. The Iris cp genomes were highly conserved and similar to other Iridaceae cp genomes. Nucleotide diversity analysis indicated that all three species had similar levels of genetic variation, but the cp genomes of I. koreana and I. minutoaurea were more similar to each other than to I. odaesanensis. Positive selection was inferred for psbK and ycf2 genes of the three Iris species. Phylogenetic analyses consistently recovered I. odaesanensis as a sister to a clade containing I. koreana and I. minutoaurea. Although the phylogenetic relationship, rDNA loci number, and localization, together with the genome size and chromosome number of the three species, allowed for the inference of I. minutoaurea as a putative maternal taxon and I. odaesanensis as a paternal taxon, further analyses involving species-specific molecular cytogenetic markers and genomic in situ hybridization are required to interpret the mechanisms involved in the origin of the chromosomal variation in Iris series Chinenses. This study contributes towards the genomic and chromosomal evolution of the genus Iris.

The grass Deschampsia cespitosa is a variable taxon out of which many varieties, subspecies and endemic species have been separated. In this paper, the variation in genome size (GS) and ploidy of this grass including several of its subspecies and two related species in Eurasia was investigated by flow cytometric (FCM) measurements. GS and ploidy data were also related to specific environments and reproduction mode. Ploidy levels found by FCM were confirmed by chromosome counts of diploid (2n = 28) and tetraploid (2n = 52) samples. Seminiferous (seed bearing) D. cespitosa was mainly diploid (GS between 3.754 and 5.438 pg/1C). GS variation in diploids showed a geographic pattern with a significant difference (H = 41,441, P < 0.001) between European (median = 4.377 pg) and Asian (median = 4.881 pg) accessions. Genome size (1C) in tetraploids ranged from 7.9426 to 9.0399 pg. Tetraploid seminiferous D. cespitosa was found mostly in disturbed habitats in western and southern Europe, while tetraploids in Asia were registered in wet Arctic habitats. Genome size (1C between 8.3278 and 8.8603 pg) of the pseudoviviparous plants (spikelets produce plantlets asexually) of wet habitats in central and northern Europe indicated tetraploidy. A putative triploid (GS 6.6817 pg) was detected in Iceland. Summing up, we found a high variation in GS on the geographic scale with significant regional differences in diploid D. cespitosa. Among the tetraploids, the asexually reproducing plants were bound to specific habitats, while the seminiferous plants showed a habitat preference similar to the diploids.

Capsicum L. encompasses 43 American species, including the five domesticated worldwide consumed sweet and hot chiles. This study presents new, updated and age-calibrated phylogenetic hypothesis on the genus evolution incorporating nearly all currently accepted Capsicum species. A new model of spatial diversification of Capsicum is proposed based on analyses of several thousands of genome-wide RAD-seq derived SNPs. Maximum likelihood approaches were used to reconstruct phylogenies and to estimate dates of all major splits. Ancestral ranges were inferred and diversification events were modeled in a time frame using a Bayesian approach. Nine clades corresponding to genetically and (mostly) geographically well-defined lineages, which diversified starting around mid-upper Miocene, were recovered with strong support. The Northern and Central Andes were inferred to represent the most likely ancestral range of the genus Capsicum. A few early vicariant and dispersal events were estimated to have driven the geographic divergence of the main Capsicum clades. Each lineage was inferred to have diversified within a distinct region of South America and expanded geographically to different extent. Extant species diversification was inferred to have begun at the beginning of Pliocene and continued through the Pleistocene. The Central Andes, represented mainly by the territory of present-day Bolivia, were proposed to play a central role in the diversification of lineages comprising domesticated capsicums and their wild allies. The genome-wide approach allowed for high resolution and support of deep phylogenetic nodes providing novel insights into the affinities of major lineages and clades as well as on the geographic expansion of Capsicum. This study provides the first dated evolutionary history of the genus encompassing most of the chile species diversity.

To provide insights into the fate of transposable elements (TEs) across timescales in a post-polyploidization context, we comparatively investigate five sibling Dactylorhiza allotetraploids (Orchidaceae) formed independently and sequentially between 500 and 100K generations ago by unidirectional hybridization between diploids D. fuchsii and D. incarnata. Our results first reveal that the paternal D. incarnata genome shows a marked increased content of LTR retrotransposons compared to the maternal species, reflected in its larger genome size and consistent with a previously hypothesized bottleneck. With regard to the allopolyploids, in the youngest D. purpurella both genome size and TE composition appear to be largely additive with respect to parents, whereas for polyploids of intermediate ages we uncover rampant genome expansion on a magnitude of multiple entire genomes of some plants such as Arabidopsis. The oldest allopolyploids in the series are not larger than the intermediate ones. A putative tandem repeat, potentially derived from a non-autonomous miniature inverted-repeat TE (MITE) drives much of the genome dynamics in the allopolyploids. The highly dynamic MITE-like element is found in higher proportions in the maternal diploid, D. fuchsii, but is observed to increase in copy number in both subgenomes of the allopolyploids. Altogether, the fate of repeats appears strongly regulated and therefore predictable across multiple independent allopolyploidization events in this system. Apart from the MITE-like element, we consistently document a mild genomic shock following the allopolyploidizations investigated here, which may be linked to their relatively large genome sizes, possibly associated with strong selection against further genome expansions.