Post-polyploid diploidization and diversification through dysploid changes

Repeated migration, interbreeding and bottlenecking shaped the phylogeography of the selfing grass Brachypodium stacei

Brachypodium stacei is the most ancestral lineage in the genus Brachypodium, a model system for grass functional genomics. B. stacei shows striking and sometimes contradictory biological and evolutionary features, including a high selfing rate yet extensive admixture, an ancient Miocene origin yet with recent evolutionary radiation, and adaptation to different dry climate conditions in its narrow distribution range. Therefore, it constitutes an ideal system to study these life history traits. We studied the phylogeography of 17 native circum-Mediterranean B. stacei populations (39 individuals) using genome-wide RADseq SNP data and complete plastome sequences. Nuclear SNP data revealed the existence of six distinct genetic clusters, low levels of intra-population genetic diversity and high selfing rates, albeit with signatures of admixture. Coalescence-based dating analysis detected a recent split between crown lineages in the Late Quaternary. Plastome sequences showed incongruent evolutionary relationships with those recovered by the nuclear data, suggesting interbreeding and chloroplast capture events between genetically distant populations. Demographic and population dispersal coalescent models identified an ancestral origin of B. stacei in the western-central Mediterranean islands, followed by an early colonization of the Canary Islands and two independent colonization events of the eastern Mediterranean region through long-distance dispersal and bottleneck events as the most likely evolutionary history. Climate niche data identified three arid niches of B. stacei in the southern Mediterranean region. Our findings indicate that the phylogeography of B. stacei populations was shaped by recent radiations, frequent extinctions, long-distance dispersal events, occasional interbreeding, and adaptation to local climates.

Biased Retention of Environment-Responsive Genes Following Genome Fractionation

The molecular underpinnings and consequences of cycles of whole-genome duplication (WGD) and subsequent gene loss through subgenome fractionation remain largely elusive. Endogenous drivers, such as transposable elements (TEs), have been postulated to shape genome-wide dominance and biased fractionation, leading to a conserved least-fractionated (LF) subgenome and a degenerated most-fractionated (MF) subgenome. In contrast, the role of exogenous factors, such as those induced by environmental stresses, has been overlooked. In this study, a chromosome-scale assembly of the alpine buckler mustard (Biscutella laevigata; Brassicaceae) that underwent a WGD event about 11 million years ago is coupled with transcriptional responses to heat, cold, drought, and herbivory to assess how gene expression is associated with differential gene retention across the MF and LF subgenomes. Counteracting the impact of TEs in reducing the expression and retention of nearby genes across the MF subgenome, dosage balance is highlighted as a main endogenous promoter of the retention of duplicated gene products under purifying selection. Consistent with the "turn a hobby into a job" model, about one-third of environment-responsive duplicates exhibit novel expression patterns, with one copy typically remaining conditionally expressed, whereas the other copy has evolved constitutive expression, highlighting exogenous factors as a major driver of gene retention. Showing uneven patterns of fractionation, with regions remaining unbiased, but with others showing high bias and significant enrichment in environment-responsive genes, this mesopolyploid genome presents evolutionary signatures consistent with an interplay of endogenous and exogenous factors having driven gene content following WGD-fractionation cycles.

Chromosome-level assembly of Lindenbergia philippensis and comparative genomic analyses shed light on genome evolution in Lamiales

Lamiales, comprising over 23,755 species across 24 families, stands as a highly diverse and prolific plant group, playing a significant role in the cultivation of horticultural, ornamental, and medicinal plant varieties. Whole-genome duplication (WGD) and its subsequent post-polyploid diploidization (PPD) process represent the most drastic type of karyotype evolution, injecting significant potential for promoting the diversity of this lineage. However, polyploidization histories, as well as genome and subgenome fractionation following WGD events in Lamiales species, are still not well investigated. In this study, we constructed a chromosome-level genome assembly of Lindenbergia philippensis (Orobanchaceae) and conducted comparative genomic analyses with 14 other Lamiales species. L. philippensis is positioned closest to the parasitic lineage within Orobanchaceae and has a conserved karyotype. Through a combination of Ks analysis and syntenic depth analysis, we reconstructed and validated polyploidization histories of Lamiales species. Our results indicated that Primulina huaijiensis underwent three rounds of diploidization events following the γ-WGT event, rather than two rounds as reported. Besides, we reconfirmed that most Lamiales species shared a common diploidization event (L-WGD). Subsequently, we constructed the Lamiales Ancestral Karyotype (LAK), comprising 11 proto-chromosomes, and elucidated its evolutionary trajectory, highlighting the highly flexible reshuffling of the Lamiales paleogenome. We identified biased fractionation of subgenomes following the L-WGD event across eight species, and highlighted the positive impacts of non-WGD genes on gene family expansion. This study provides novel genomic resources and insights into polyploidy and karyotype remodeling of Lamiales species, essential for advancing our understanding of species diversification and genome evolution.

The link between ancient whole-genome duplications and cold adaptations in the Caryophyllaceae

PREMISE

The Caryophyllaceae (the carnation family) have undergone multiple transitions into colder climates and convergence on cushion plant adaptation, indicating that they may provide a natural system for cold adaptation research. Previous research has suggested that putative ancient whole-genome duplications (WGDs) are correlated with niche shifts into colder climates across the Caryophyllales. Here, we explored the genomic changes potentially involved in one of these discovered shifts in the Caryophyllaceae.

METHODS

We constructed a data set combining 26 newly generated transcriptomes with 45 published transcriptomes, including 11 cushion plant species across seven genera. With this data set, we inferred a dated phylogeny for the Caryophyllaceae and mapped ancient WGDs and gene duplications onto the phylogeny. We also examined functional groups enriched for gene duplications related to the climatic shift.

RESULTS

The ASTRAL topology was mostly congruent with the current consensus of relationships within the family. We inferred 15 putative ancient WGDs in the family, including eight that have not been previously published. The oldest ancient WGD (ca. 64.4-56.7 million years ago), WGD1, was found to be associated with a shift into colder climates by previous research. Gene regions associated with ubiquitination were overrepresented in gene duplications retained after WGD1 and those convergently retained by cushion plants in Colobanthus and Eremogone, along with other functional annotations.

CONCLUSIONS

Gene family expansions induced by ancient WGDs may have contributed to the shifts to cold climatic niches in the Caryophyllaceae. Transcriptomic data are crucial resources that help unravel heterogeneity in deep-time evolutionary patterns in plants.

Investigating historical drivers of latitudinal gradients in polyploid plant biogeography: A multiclade perspective

PREMISE

The proportion of polyploid plants in a community increases with latitude, and different hypotheses have been proposed about which factors drive this pattern. Here, we aimed to understand the historical causes of the latitudinal polyploidy gradient using a combination of ancestral state reconstruction methods. Specifically, we assessed whether (1) polyploidization enables movement to higher latitudes (i.e., polyploidization precedes occurrences in higher latitudes) or (2) higher latitudes facilitate polyploidization (i.e., occurrence in higher latitudes precedes polyploidization).

METHODS

We reconstructed the ploidy states and ancestral niches of 1032 angiosperm species at four paleoclimatic time slices ranging from 3.3 million years ago to the present, comprising taxa from four well-represented clades: Onagraceae, Primulaceae, Solanum (Solanaceae), and Pooideae (Poaceae). We used ancestral niche reconstruction models alongside a customized discrete character evolution model to allow reconstruction of states at specific time slices. Patterns of latitudinal movement were reconstructed and compared in relation to inferred ploidy shifts.

RESULTS

No single hypothesis applied equally well across all analyzed clades. While significant differences in median latitudinal occurrence were detected in the largest clade, Poaceae, no significant differences were detected in latitudinal movement in any clade.

CONCLUSIONS

Our preliminary study is the first to attempt to connect ploidy changes to continuous latitudinal movement, but we cannot favor one hypothesis over another. Given that patterns seem to be clade-specific, more clades must be analyzed in future studies for generalities to be drawn.

Polyploids of Brassicaceae: Genomic Insights and Assembly Strategies

The Brassicaceae family is distinguished by its inclusion of high-value crops such as cabbage, broccoli, mustard, and wasabi, all noted for their glucosinolates. In this family, many polyploidy species are distributed and shaped by numerous whole-genome duplications, independent genome doublings, and hybridization events. The evolutionary trajectory of the family is marked by enhanced diversification and lineage splitting after paleo- and meso-polyploidization, with discernible remnants of whole-genome duplications within their genomes. The recent neopolyploidization events notably increased the proportion of polyploid species within the family. Although sequencing efforts for the Brassicaceae genome have been robust, accurately distinguishing sub-genomes remains a significant challenge, frequently complicating the assembly process. Assembly strategies include comparative analyses with ancestral species and examining k-mers, long terminal repeat retrotransposons, and pollen sequencing. This review comprehensively explores the unique genomic characteristics of the Brassicaceae family, with a particular emphasis on polyploidization events and the latest strategies for sequencing and assembly. This review will significantly improve our understanding of polyploidy in the Brassicaceae family and assist in future genome assembly methods.

Striking variation in chromosome structure within Musa acuminata subspecies, diploid cultivars, and F1 diploid hybrids

The majority of cultivated bananas originated from inter- and intra(sub)specific crosses between two wild diploid species, Musa acuminata and Musa balbisiana. Hybridization and polyploidization events during the evolution of bananas led to the formation of clonally propagated cultivars characterized by a high level of genome heterozygosity and reduced fertility. The combination of low fertility in edible clones and differences in the chromosome structure among M. acuminata subspecies greatly hampers the breeding of improved banana cultivars. Using comparative oligo-painting, we investigated large chromosomal rearrangements in a set of wild M. acuminata subspecies and cultivars that originated from natural and human-made crosses. Additionally, we analyzed the chromosome structure of F1 progeny that resulted from crosses between Mchare bananas and the wild M. acuminata 'Calcutta 4' genotype. Analysis of chromosome structure within M. acuminata revealed the presence of a large number of chromosomal rearrangements showing a correlation with banana speciation. Chromosome painting of F1 hybrids was complemented by Illumina resequencing to identify the contribution of parental subgenomes to the diploid hybrid clones. The balanced presence of both parental genomes was revealed in all F1 hybrids, with the exception of one clone, which contained only Mchare-specific SNPs and thus most probably originated from an unreduced diploid gamete of Mchare.

Uniparental silencing of 5S rRNA genes in plant allopolyploids - insights from Cardamine (Brassicaceae)

While the phenomenon of uniparental silencing of 35S rDNA in interspecific hybrids and allopolyploids is well documented, there is a notable absence of information regarding whether such silencing extends to the 5S RNA component of ribosomes. To address this gap in knowledge, we analyzed the 5S and 35S rDNA expression in Cardamine (Brassicaceae) allopolyploids, namely C. × insueta (2n = 3x = 24, genome composition RRA), C. flexuosa (2n = 4x = 32, AAHH), and C. scutata (2n = 4x = 32, PPAA) which share a common diploid ancestor (AA). We employed high-throughput sequencing of transcriptomes and genomes and phylogenetic analyses of 5S rRNA variants. The genomic organization of rDNA was further scrutinized through clustering and fluorescence in situ hybridization. In the C. × insueta allotriploid, we observed uniparental dominant expression of 5S and 35S rDNA loci. In the C. flexuosa and C. scutata allotetraploids, the expression pattern differed, with the 35S rDNA being expressed from the A subgenome, whereas the 5S rDNA was expressed from the partner subgenome. Both C. flexuosa and C. scutata but not C. × insueta showed copy and locus number changes. We conclude that in stabilized allopolyploids, transcription of ribosomal RNA components occurs from different subgenomes. This phenomenon appears to result in the formation of chimeric ribosomes comprising rRNA molecules derived from distinct parental origins. We speculate that the interplay of epigenetic silencing and rDNA rearrangements introduces an additional layer of variation in multimolecule ribosomal complexes, potentially contributing to the evolutionary success of allopolyploids.

First investigation into the genetic control of meiosis in sugarcane

The sugarcane (Saccharum spp.) genome is one of the most complex of all. Modern varieties are highly polyploid and aneuploid as a result of hybridization between Saccharum officinarum and S. spontaneum. Little research has been done on meiotic control in polyploid species, with the exception of the wheat Ph1 locus harboring the ZIP4 gene (TaZIP4-B2) which promotes pairing between homologous chromosomes while suppressing crossover between homeologs. In sugarcane, despite its interspecific origin, bivalent association is favored, and multivalents, if any, are resolved at the end of prophase I. Thus, our aim herein was to investigate the purported genetic control of meiosis in the parental species and in sugarcane itself. We investigated the ZIP4 gene and immunolocalized meiotic proteins, namely synaptonemal complex proteins Zyp1 and Asy1. The sugarcane ZIP4 gene is located on chromosome 2 and expressed more abundantly in flowers, a similar profile to that found for TaZIP4-B2. ZIP4 expression is higher in S. spontaneum a neoautopolyploid, with lower expression in S. officinarum, a stable octoploid species. The sugarcane Zip4 protein contains a TPR domain, essential for scaffolding. Its 3D structure was also predicted, and it was found to be very similar to that of TaZIP4-B2, reflecting their functional relatedness. Immunolocalization of the Asy1 and Zyp1 proteins revealed that S. officinarum completes synapsis. However, in S. spontaneum and SP80-3280 (a modern variety), no nuclei with complete synapsis were observed. Importantly, our results have implications for sugarcane cytogenetics, genetic mapping, and genomics.

Phylogenomics and morphology of Celmisiinae (Asteraceae: Astereae): Taxonomic and evolutionary implications

The tribe Astereae (Asteraceae) includes 36 subtribes and 252 genera, and is distributed worldwide in temperate and tropical regions. One of the subtribes, Celmisiinae Saldivia, has been recently circumscribed to include six genera and ca. 160 species, and is restricted to eastern Australia, New Zealand, and New Guinea. The species show an impressive range of growth habit, from small herbs and ericoid subshrubs to medium-sized trees. They live in a wide range of habitats and are often dominant in subalpine and alpine vegetation. Despite the well-supported circumscription of Celmisiinae, uncertainties have remained about their internal relationships and classification at genus and species levels. This study exploited recent advances in high-throughput sequencing to build a robust multi-gene phylogeny for the subtribe Celmisiinae. The target enrichment Angiosperms353 bait set and the hybpiper-nf and paragone-nf pipelines were used to retrieve, infer, and assemble orthologous loci from 75 taxa representing all the main putative clades within the subtribe. Because of the diploidised ploidy level in Celmisiinae, as well as missing data in the assemblies, uncertainty remains surrounding the inference of orthology detection. However, based on a variety of gene-family sets, coalescent and concatenation-based phylogenetic reconstructions recovered similar topologies. Paralogy and missing data in the gene-families caused some problems, but the estimated phylogenies were well-supported and well-resolved. The phylogenomic evidence supported Celmisiinae and three main clades: the Pleurophyllum clade (Pleurophyllum, Macrolearia and Damnamenia), mostly in the New Zealand Subantarctic Islands, Celmisia of mainland New Zealand and Australia, and Shawia (including 'Olearia pro parte' and Pachystegia) of New Zealand, Australia and New Guinea. The results presented here add to the accumulating support for the Angiosperms353 bait set as an efficient method for documenting plant diversity.

Regulation of gene-edited plants in Europe: from the valley of tears into the shining sun?

Some 20 years ago, the EU introduced complex regulatory rules for the growth of transgenic crops, which resulted in a de facto ban to grow these plants in fields within most European countries. With the rise of novel genome editing technologies, it has become possible to improve crops genetically in a directed way without the need for incorporation of foreign genes. Unfortunately, in 2018, the European Court of Justice ruled that such gene-edited plants are to be regulated like transgenic plants. Since then, European scientists and breeders have challenged this decision and requested a revision of this outdated law. Finally, after 5 years, the European Commission has now published a proposal on how, in the future, to regulate crops produced by new breeding technologies. The proposal tries to find a balance between the different interest groups in Europe. On one side, genetically modified plants, which cannot be discerned from their natural counterparts, will exclusively be used for food and feed and are-besides a registration step-not to be regulated at all. On the other side, plants expressing herbicide resistance are to be excluded from this regulation, a concession to the strong environmental associations and NGOs in Europe. Moreover, edited crops are to be excluded from organic farming to protect the business interests of the strong organic sector in Europe. Nevertheless, if this law passes European parliament and council, unchanged, it will present a big step forward toward establishing a more sustainable European agricultural system. Thus, it might soon be possible to develop and grow crops that are more adapted to global warming and whose cultivation will require lower amounts of pesticides. However, there is still a long way to go until the law is passed. Too often, the storm of arguments raised by the opponents, based on irrational fears of mutations and a naive understanding of nature, has fallen on fruitful ground in Europe.

A locus-dependent mixed inheritance in the segmental allohexaploid sweetpotato (Ipomoea batatas [L.] Lam)

Two interrelated aspects of the sweetpotato genome, its polyploid origin and inheritance type, remain uncertain. We recently proposed a segmental allohexaploid sweetpotato and thus sought to clarify its inheritance type by direct analyses of homoeolog segregations at selected single-copy loci. For such analyses, we developed a digital quantitative PCR genotyping method using one nondiscriminatory and three discriminatory probes for each selected locus to discriminate and quantify three homoeolog-differentiating variation types (homoeolog-types) in genomic DNA samples for genotype fitting and constructed a F2 population for segregation analyses. We confirmed inter-subgenomic distinctions of three identified homoeolog-types at each of five selected loci by their interspecific differentiations among 14 species in Ipomoea section batatas and genotyped the loci in 549 F2 lines, selected F1 progenies, and their founding parents. Segregation and genotype analyses revealed a locus-dependent mixed inheritance (disomic, polysomic, and intermediate types) of the homoeolog-types at 4 loci in the F2 population, displaying estimated disomic-inheritance frequencies of 0, 2.72%, 14.52%, and 36.92%, and probably in the F1 population too. There were also low-frequency non-hexaploid F1 and F2 genotypes that were probably derived from double-reduction recombination or partially unreduced gametes, and F2 genotypes of apparent aneuploids/dysploids with neopolyploid-like frequencies. Additional analyses of homoeolog-type genotypes at the 5 loci in 46 lines from various regions revealed locus-dependent selection biases, favoring genotypes having more of one homoeolog-type, i.e. more of di- or homogenized homoeolog-type composition, and one-direction ploidy trending among apparent aneuploids/dysploids. These inheritance features pointed to an evolving segmental allohexaploid sweetpotato impacted by selection biases.

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.

The genetic consequences of range expansion and its influence on diploidization in polyploids

Despite newly formed polyploids being subjected to myriad fitness consequences, the relative prevalence of polyploidy both contemporarily and in ancestral branches of the tree of life suggests alternative advantages that outweigh these consequences. One proposed advantage is that polyploids may more easily colonize novel habitats such as deglaciated areas. However, previous research conducted in diploids suggests that range expansion comes with a fitness cost as deleterious mutations may fix rapidly on the expansion front. Here, we interrogate the potential consequences of expansion in polyploids by conducting spatially explicit forward-in-time simulations to investigate how ploidy and inheritance patterns impact the relative ability of polyploids to expand their range. We show that under realistic dominance models, autopolyploids suffer greater fitness reductions than diploids as a result of range expansion due to the fixation of increased mutational load that is masked in the range core. Alternatively, the disomic inheritance of allopolyploids provides a shield to this fixation resulting in minimal fitness consequences. In light of this advantage provided by disomy, we investigate how range expansion may influence cytogenetic diploidization through the reversion to disomy in autotetraploids. We show that under a wide range of parameters investigated for two models of diploidization, disomy frequently evolves more rapidly on the expansion front than in the range core, and that this dynamic inheritance model has additional effects on fitness. Together our results point to a complex interaction between dominance, ploidy, inheritance, and recombination on fitness as a population spreads across a geographic range.

The global distribution of angiosperm genome size is shaped by climate

Angiosperms, which inhabit diverse environments across all continents, exhibit significant variation in genome sizes, making them an excellent model system for examining hypotheses about the global distribution of genome size. These include the previously proposed large genome constraint, mutational hazard, polyploidy-mediated, and climate-mediated hypotheses. We compiled the largest genome size dataset to date, encompassing 16 017 (> 5% of known) angiosperm species, and analyzed genome size distribution using a comprehensive geographic distribution dataset for all angiosperms. We observed that angiosperms with large range sizes generally had small genomes, supporting the large genome constraint hypothesis. Climate was shown to exert a strong influence on genome size distribution along the global latitudinal gradient, while the frequency of polyploidy and the type of growth form had negligible effects. In contrast to the unimodal patterns along the global latitudinal gradient shown by plant size traits and polyploid proportions, the increase in angiosperm genome size from the equator to 40-50°N/S is probably mediated by different (mostly climatic) mechanisms than the decrease in genome sizes observed from 40 to 50°N northward. Our analysis suggests that the global distribution of genome sizes in angiosperms is mainly shaped by climatically mediated purifying selection, genetic drift, relaxed selection, and environmental filtering.

Genome architecture and genetic diversity of allopolyploid okra (Abelmoschus esculentus)

The allopolyploid okra (Abelmoschus esculentus) unveiled telomeric repeats flanking distal gene-rich regions and short interstitial TTTAGGG telomeric repeats, possibly representing hallmarks of chromosomal speciation. Ribosomal RNA (rRNA) genes organize into 5S clusters, distinct from the 18S-5.8S-28S units, indicating an S-type rRNA gene arrangement. The assembly, in line with cytogenetic and cytometry observations, identifies 65 chromosomes and a 1.45 Gb genome size estimate in a haploid sibling. The lack of aberrant meiotic configurations implies limited to no recombination among sub-genomes. k-mer distribution analysis reveals 75% has a diploid nature and 15% heterozygosity. The configurations of Benchmarking Universal Single-Copy Ortholog (BUSCO), k-mer, and repeat clustering point to the presence of at least two sub-genomes one with 30 and the other with 35 chromosomes, indicating the allopolyploid nature of the okra genome. Over 130 000 putative genes, derived from mapped IsoSeq data and transcriptome data from public okra accessions, exhibit a low genetic diversity of one single nucleotide polymorphisms per 2.1 kbp. The genes are predominantly located at the distal chromosome ends, declining toward central scaffold domains. Long terminal repeat retrotransposons prevail in central domains, consistent with the observed pericentromeric heterochromatin and distal euchromatin. Disparities in paralogous gene counts suggest potential sub-genome differentiation implying possible sub-genome dominance. Amino acid query sequences of putative genes facilitated phenol biosynthesis pathway annotation. Comparison with manually curated reference KEGG pathways from related Malvaceae species reveals the genetic basis for putative enzyme coding genes that likely enable metabolic reactions involved in the biosynthesis of dietary and therapeutic compounds in okra.

Genomic evidence for rediploidization and adaptive evolution following the whole-genome triplication

Whole-genome duplication (WGD), or polyploidy, events are widespread and significant in the evolutionary history of angiosperms. However, empirical evidence for rediploidization, the major process where polyploids give rise to diploid descendants, is still lacking at the genomic level. Here we present chromosome-scale genomes of the mangrove tree Sonneratia alba and the related inland plant Lagerstroemia speciosa. Their common ancestor has experienced a whole-genome triplication (WGT) approximately 64 million years ago coinciding with a period of dramatic global climate change. Sonneratia, adapting mangrove habitats, experienced extensive chromosome rearrangements post-WGT. We observe the WGT retentions display sequence and expression divergence, suggesting potential neo- and sub-functionalization. Strong selection acting on three-copy retentions indicates adaptive value in response to new environments. To elucidate the role of ploidy changes in genome evolution, we improve a model of the polyploidization-rediploidization process based on genomic evidence, contributing to the understanding of adaptive evolution during climate change.

Subgenome evolutionary dynamics in allotetraploid ferns: insights from the gene expression patterns in the allotetraploid species Phegopteris decursivepinnata (Thelypteridacea, Polypodiales)

Allopolyploidization often leads to disruptive conflicts among more than two sets of subgenomes, leading to genomic modifications and changes in gene expression. Although the evolutionary trajectories of subgenomes in allopolyploids have been studied intensely in angiosperms, the dynamics of subgenome evolution remain poorly understood in ferns, despite the prevalence of allopolyploidization. In this study, we have focused on an allotetraploid fern-Phegopteris decursivepinnata-and its diploid parental species, P. koreana (K) and P. taiwaniana (T). Using RNA-seq analyses, we have compared the gene expression profiles for 9,540 genes among parental species, synthetic F1 hybrids, and natural allotetraploids. The changes in gene expression patterns were traced from the F1 hybrids to the natural allopolyploids. This study has revealed that the expression patterns observed in most genes in the F1 hybrids are largely conserved in the allopolyploids; however, there were substantial differences in certain genes between these groups. In the allopolyploids compared with the F1 hybrids, the number of genes showing a transgressive pattern in total expression levels was increased. There was a slight reduction in T-dominance and a slight increase in K-dominance, in terms of expression level dominance. Interestingly, there is no obvious bias toward the T- or K-subgenomes in the number and expression levels overall, showing the absence of subgenome dominance. These findings demonstrated the impacts of the substantial transcriptome change after hybridization and the moderate modification during allopolyploid establishment on gene expression in ferns and provided important insights into subgenome evolution in polyploid ferns.

Plant chromosome engineering - past, present and future

Spontaneous chromosomal rearrangements (CRs) play an essential role in speciation, genome evolution and crop domestication. To be able to use the potential of CRs for breeding, plant chromosome engineering was initiated by fragmenting chromosomes by X-ray irradiation. With the rise of the CRISPR/Cas system, it became possible to induce double-strand breaks (DSBs) in a highly efficient manner at will at any chromosomal position. This has enabled a completely new level of predesigned chromosome engineering. The genetic linkage between specific genes can be broken by inducing chromosomal translocations. Natural inversions, which suppress genetic exchange, can be reverted for breeding. In addition, various approaches for constructing minichromosomes by downsizing regular standard A or supernumerary B chromosomes, which could serve as future vectors in plant biotechnology, have been developed. Recently, a functional synthetic centromere could be constructed. Also, different ways of genome haploidization have been set up, some based on centromere manipulations. In the future, we expect to see even more complex rearrangements, which can be combined with previously developed engineering technologies such as recombinases. Chromosome engineering might help to redefine genetic linkage groups, change the number of chromosomes, stack beneficial genes on mini cargo chromosomes, or set up genetic isolation to avoid outcrossing.

Cytological analysis of the diploid-like inheritance of newly synthesized allotetraploid wheat

Polyploidization is a process which is related to species hybridization and whole genome duplication. It is widespread among angiosperm evolution and is essential for speciation and diversification. Allopolyploidization is mainly derived from interspecific hybridization and is believed to pose chromosome imbalances and genome instability caused by meiotic irregularity. However, the self-compatible allopolyploid in wild nature is cytogenetically and genetically stable. Whether this stabilization form was achieved in initial generation or a consequence of long term of evolution was largely unknown. Here, we synthesized a series of nascent allotetraploid wheat derived from three diploid genomes of A, S*, and D. The chromosome numbers of the majority of the progeny derived from these newly formed allotetraploid wheat plants were found to be relatively consistent, with each genome containing 14 chromosomes. In meiosis, bivalent was the majority of the chromosome configuration in metaphase I which supports the stable chromosome number inheritance in the nascent allotetraploid. These findings suggest that diploidization occurred in the newly formed synthetic allotetraploid wheat. However, we still detected aneuploids in a proportion of newly formed allotetraploid wheat, and meiosis of these materials present more irregular chromosome behavior than the euploid. We found that centromere pairing and centromere clustering in meiosis was affected in the aneuploids, which suggest that aneuploidy may trigger the irregular interactions of centromere in early meiosis which may take participate in promoting meiosis stabilization in newly formed allotetraploid wheat.

Chromosomal dynamics in Senna: comparative PLOP-FISH analysis of tandem repeats and flow cytometric nuclear genome size estimations

Introduction

Tandem repeats (TRs) occur abundantly in plant genomes. They play essential roles that affect genome organization and evolution by inducing or generating chromosomal rearrangements such as duplications, deletions, inversions, and translocations. These impact gene expression and chromosome structure and even contribute to the emergence of new species.

Method

We investigated the effects of TRs on speciation in Senna genus by performing a comparative analysis using fluorescence in situ hybridization (FISH) with S. tora-specific TR probes. We examined the chromosomal distribution of these TRs and compared the genome sizes of seven Senna species (estimated using flow cytometry) to better understand their evolutionary relationships.

Results

Two (StoTR03_159 and StoTR04_55) of the nine studied TRs were not detected in any of the seven Senna species, whereas the remaining seven were found in all or some species with patterns that were similar to or contrasted with those of S. tora. Of these studies species, only S. angulata showed significant genome rearrangements and dysploid karyotypes resembling those of S. tora. The genome sizes varied among these species and did not positively correlate with chromosome number. Notably, S. angulata had the fewest chromosomes (2n = 22) but a relatively large genome size.

Discussion

These findings reveal the dynamics of TRs and provide a cytogenetic depiction of chromosomal rearrangements during speciation in Senna. To further elucidate the dynamics of repeat sequences in Senna, future studies must include related species and extensive repeatomic studies, including those on transposable elements.

The Impact of Chromosomal Rearrangements in Speciation: From Micro- to Macroevolution

Chromosomal rearrangements (CRs) have been known since almost the beginning of genetics. While an important role for CRs in speciation has been suggested, evidence primarily stems from theoretical and empirical studies focusing on the microevolutionary level (i.e., on taxon pairs where speciation is often incomplete). Although the role of CRs in eukaryotic speciation at a macroevolutionary level has been supported by associations between species diversity and rates of evolution of CRs across phylogenies, these findings are limited to a restricted range of CRs and taxa. Now that more broadly applicable and precise CR detection approaches have become available, we address the challenges in filling some of the conceptual and empirical gaps between micro- and macroevolutionary studies on the role of CRs in speciation. We synthesize what is known about the macroevolutionary impact of CRs and suggest new research avenues to overcome the pitfalls of previous studies to gain a more comprehensive understanding of the evolutionary significance of CRs in speciation across the tree of life.

Global Brassicaceae phylogeny based on filtering of 1,000-gene dataset

The mustard family (Brassicaceae) is a scientifically and economically important family, containing the model plant Arabidopsis thaliana and numerous crop species that feed billions worldwide. Despite its relevance, most phylogenetic trees of the family are incompletely sampled and often contain poorly supported branches. Here, we present the most complete Brassicaceae genus-level family phylogenies to date (Brassicaceae Tree of Life or BrassiToL) based on nuclear (1,081 genes, 319 of the 349 genera; 57 of the 58 tribes) and plastome (60 genes, 265 genera; all tribes) data. We found cytonuclear discordance between the two, which is likely a result of rampant hybridization among closely and more distantly related lineages. To evaluate the impact of such hybridization on the nuclear phylogeny reconstruction, we performed five different gene sampling routines, which increasingly removed putatively paralog genes. Our cleaned subset of 297 genes revealed high support for the tribes, whereas support for the main lineages (supertribes) was moderate. Calibration based on the 20 most clock-like nuclear genes suggests a late Eocene to late Oligocene origin of the family. Finally, our results strongly support a recently published new family classification, dividing the family into two subfamilies (one with five supertribes), together representing 58 tribes. This includes five recently described or re-established tribes, including Arabidopsideae, a monogeneric tribe accommodating Arabidopsis without any close relatives. With a worldwide community of thousands of researchers working on Brassicaceae and its diverse members, our new genus-level family phylogeny will be an indispensable tool for studies on biodiversity and plant biology.

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.

The meso-octoploid Heliophila variabilis genome sheds a new light on the impact of polyploidization and diploidization on the diversity of the Cape flora

Although the South African Cape flora is one of the most remarkable biodiversity hotspots, its high diversity has not been associated with polyploidy. Here, we report the chromosome-scale genome assembly of an ephemeral cruciferous species Heliophila variabilis (~334 Mb, n = 11) adapted to South African semiarid biomes. Two pairs of differently fractionated subgenomes suggest an allo-octoploid origin of the genome at least 12 million years ago. The ancestral octoploid Heliophila genome (2n = 8x = ~60) has probably originated through hybridization between two allotetraploids (2n = 4x = ~30) formed by distant, intertribal, hybridization. Rediploidization of the ancestral genome was marked by extensive reorganization of parental subgenomes, genome downsizing, and speciation events in the genus Heliophila. We found evidence for loss-of-function changes in genes associated with leaf development and early flowering, and over-retention and sub/neofunctionalization of genes involved in pathogen response and chemical defense. The genomic resources of H. variabilis will help elucidate the role of polyploidization and genome diploidization in plant adaptation to hot arid environments and origin of the Cape flora. The sequenced H. variabilis represents the first chromosome-scale genome assembly of a meso-octoploid representative of the mustard family.

High rates of structural rearrangements have shaped the chromosome evolution in dysploid Phaseolus beans

KEY MESSAGE

Karyotypes evolve through numerical and structural chromosome rearrangements. We show that Phaseolus leptostachyus, a wild bean, underwent a rapid genome reshuffling associated with the reduction from 11 to 10 chromosome pairs, but without whole genome duplication, the highest chromosome evolution rate known for plants. Plant karyotypes evolve through structural rearrangements often associated with polyploidy or dysploidy. The genus Phaseolus comprises ~ 90 species, five of them domesticated due to their nutritional relevance. Most of the species have 2n = 22 karyotypes and are highly syntenic, except for three dysploid karyotypes of species from the Leptostachyus group (2n = 20) that have accumulated several rearrangements. Here, we investigated the degrees of structural rearrangements among Leptostachyus and other Phaseolus groups by estimating their chromosomal evolution rates (CER). For this, we combined our oligo-FISH barcode system for beans and chromosome-specific painting probes for chromosomes 2 and 3, with rDNA and a centromeric probe to establish chromosome orthologies and identify structural rearrangements across nine Phaseolus species. We also integrated the detected rearrangements with a phylogenomic approach to estimate the CERs for each Phaseolus lineage. Our data allowed us to identify translocations, inversions, duplications and deletions, mostly in species belonging to the Leptostachyus group. Phaseolus leptostachyus showed the highest CER (12.31 rearrangements/My), a tenfold increase in contrast to the 2n = 22 species analysed. This is the highest rate known yet for plants, making it a model species for investigating the mechanisms behind rapid genome reshuffling in early species diversification.

Overcoming genetic paucity of Camelina sativa: possibilities for interspecific hybridization conditioned by the genus evolution pathway

Camelina or false flax (Camelina sativa) is an emerging oilseed crop and a feedstock for biofuel production. This species is believed to originate from Western Asian and Eastern European regions, where the center of diversity of the Camelina genus is located. Cultivated Camelina species arose via a series of polyploidization events, serving as bottlenecks narrowing genetic diversity of the species. The genetic paucity of C. sativa is foreseen as the most crucial limitation for successful breeding and improvement of this crop. A potential solution to this challenge could be gene introgression from Camelina wild species or from resynthesized allohexaploid C. sativa. However, both approaches would require a complete comprehension of the evolutionary trajectories that led to the C. sativa origin. Although there are some studies discussing the origin and evolution of Camelina hexaploid species, final conclusions have not been made yet. Here, we propose the most complete integrated evolutionary model for the Camelina genus based on the most recently described findings, which enables efficient improvement of C. sativa via the interspecific hybridization with its wild relatives. We also discuss issues of interspecific and intergeneric hybridization, aimed on improving C. sativa and overcoming the genetic paucity of this crop. The proposed comprehensive evolutionary model of Camelina species indicates that a newly described species Camelina neglecta has a key role in origin of tetra- and hexaploids, all of which have two C. neglecta-based subgenomes. Understanding of species evolution within the Camelina genus provides insights into further research on C. sativa improvements via gene introgression from wild species, and a potential resynthesis of this emerging oilseed crop.

The Abundant and Unique Transcripts and Alternative Splicing of the Artificially Autododecaploid London Plane (Platanus × acerifolia)

Transcription and alternative splicing (AS) are now appreciated in plants, but few studies have examined the effects of changing ploidy on transcription and AS. In this study, we showed that artificially autododecaploid plants of London plane (Platanus × acerifolia (Aiton) Willd) had few flowers relative to their hexaploid progenitors. Transcriptome analysis based on full-length Oxford Nanopore Technologies (ONTs) and next-generation sequencing (NGS) revealed that the increased ploidy level in P. × acerifolia led to more transcribed isoforms, accompanied by an increase in the number of isoforms per gene. The functional enrichment of genes indicated that novel genes transcribed specifically in the dodecaploids may have been highly correlated with the ability to maintain genome stability. The dodecaploids showed a higher number of genes with upregulated differentially expressed genes (DEGs) compared with the hexaploid counterpart. The genome duplication of P. × acerifolia resulted mainly in the DEGs involved in basic biological pathways. It was noted that there was a greater abundance of alternative splicing (AS) events and AS genes in the dodecaploids compared with the hexaploids in P. × acerifolia. In addition, a significant difference between the structure and expression of AS events between the hexaploids and dodecaploids of Platanus was found. Of note, some DEGs and differentially spliced genes (DSGs) related to floral transition and flower development were consistent with the few flower traits in the dodecaploids of P. × acerifolia. Collectively, our findings explored the difference in transcription and AS regulation between the hexaploids and dodecaploids of P. × acerifolia and gained new insight into the molecular mechanisms underlying the few-flower phenotype of P. × acerifolia. These results contribute to uncovering the regulatory role of transcription and AS in polyploids and breeding few-flower germplasms.

Genomics of plant speciation

Studies of plants have been instrumental for revealing how new species originate. For several decades, botanical research has complemented and, in some cases, challenged concepts on speciation developed via the study of other organisms while also revealing additional ways in which species can form. Now, the ability to sequence genomes at an unprecedented pace and scale has allowed biologists to settle decades-long debates and tackle other emerging challenges in speciation research. Here, we review these recent genome-enabled developments in plant speciation. We discuss complications related to identification of reproductive isolation (RI) loci using analyses of the landscape of genomic divergence and highlight the important role that structural variants have in speciation, as increasingly revealed by new sequencing technologies. Further, we review how genomics has advanced what we know of some routes to new species formation, like hybridization or whole-genome duplication, while casting doubt on others, like population bottlenecks and genetic drift. While genomics can fast-track identification of genes and mutations that confer RI, we emphasize that follow-up molecular and field experiments remain critical. Nonetheless, genomics has clarified the outsized role of ancient variants rather than new mutations, particularly early during speciation. We conclude by highlighting promising avenues of future study. These include expanding what we know so far about the role of epigenetic and structural changes during speciation, broadening the scope and taxonomic breadth of plant speciation genomics studies, and synthesizing information from extensive genomic data that have already been generated by the plant speciation community.

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.].

Genomic Insights into Adaptation to Karst Limestone and Incipient Speciation in East Asian Platycarya spp. (Juglandaceae)

When challenged by similar environmental conditions, phylogenetically distant taxa often independently evolve similar traits (convergent evolution). Meanwhile, adaptation to extreme habitats might lead to divergence between taxa that are otherwise closely related. These processes have long existed in the conceptual sphere, yet molecular evidence, especially for woody perennials, is scarce. The karst endemic Platycarya longipes and its only congeneric species, Platycarya strobilacea, which is widely distributed in the mountains in East Asia, provide an ideal model for examining the molecular basis of both convergent evolution and speciation. Using chromosome-level genome assemblies of both species, and whole-genome resequencing data from 207 individuals spanning their entire distribution range, we demonstrate that P. longipes and P. strobilacea form two species-specific clades, which diverged around 2.09 million years ago. We find an excess of genomic regions exhibiting extreme interspecific differentiation, potentially due to long-term selection in P. longipes, likely contributing to the incipient speciation of the genus Platycarya. Interestingly, our results unveil underlying karst adaptation in both copies of the calcium influx channel gene TPC1 in P. longipes. TPC1 has previously been identified as a selective target in certain karst-endemic herbs, indicating a convergent adaptation to high calcium stress among karst-endemic species. Our study reveals the genic convergence of TPC1 among karst endemics and the driving forces underneath the incipient speciation of the two Platycarya lineages.

Patterns, mechanisms, and consequences of homoeologous exchange in allopolyploid angiosperms: a genomic and epigenomic perspective

Allopolyploids result from hybridization between different evolutionary lineages coupled with genome doubling. Homoeologous chromosomes (chromosomes with common shared ancestry) may undergo recombination immediately after allopolyploid formation and continue over successive generations. The outcome of this meiotic pairing behavior is dynamic and complex. Homoeologous exchanges (HEs) may lead to the formation of unbalanced gametes, reduced fertility, and selective disadvantage. By contrast, HEs could act as sources of novel evolutionary substrates, shifting the relative dosage of parental gene copies, generating novel phenotypic diversity, and helping the establishment of neo-allopolyploids. However, HE patterns vary among lineages, across generations, and even within individual genomes and chromosomes. The causes and consequences of this variation are not fully understood, though interest in this evolutionary phenomenon has increased in the last decade. Recent technological advances show promise in uncovering the mechanistic basis of HEs. Here, we describe recent observations of the common patterns among allopolyploid angiosperm lineages, underlying genomic and epigenomic features, and consequences of HEs. We identify critical research gaps and discuss future directions with far-reaching implications in understanding allopolyploid evolution and applying them to the development of important phenotypic traits of polyploid crops.

Independent rediploidization masks shared whole genome duplication in the sturgeon-paddlefish ancestor

Whole genome duplication (WGD) is a dramatic evolutionary event generating many new genes and which may play a role in survival through mass extinctions. Paddlefish and sturgeon are sister lineages that both show genomic evidence for ancient WGD. Until now this has been interpreted as two independent WGD events due to a preponderance of duplicate genes with independent histories. Here we show that although there is indeed a plurality of apparently independent gene duplications, these derive from a shared genome duplication event occurring well over 200 million years ago, likely close to the Permian-Triassic mass extinction period. This was followed by a prolonged process of reversion to stable diploid inheritance (rediploidization), that may have promoted survival during the Triassic-Jurassic mass extinction. We show that the sharing of this WGD is masked by the fact that paddlefish and sturgeon lineage divergence occurred before rediploidization had proceeded even half-way. Thus, for most genes the resolution to diploidy was lineage-specific. Because genes are only truly duplicated once diploid inheritance is established, the paddlefish and sturgeon genomes are thus a mosaic of shared and non-shared gene duplications resulting from a shared genome duplication event.

The evolution of the hypotetraploid Catolobus pendulus genome - the poorly known sister species of Capsella

The establishment of Arabidopsis as the most important plant model has also brought other crucifer species into the spotlight of comparative research. While the genus Capsella has become a prominent crucifer model system, its closest relative has been overlooked. The unispecific genus Catolobus is native to temperate Eurasian woodlands, from eastern Europe to the Russian Far East. Here, we analyzed chromosome number, genome structure, intraspecific genetic variation, and habitat suitability of Catolobus pendulus throughout its range. Unexpectedly, all analyzed populations were hypotetraploid (2n = 30, ~330 Mb). Comparative cytogenomic analysis revealed that the Catolobus genome arose by a whole-genome duplication in a diploid genome resembling Ancestral Crucifer Karyotype (ACK, n = 8). In contrast to the much younger Capsella allotetraploid genomes, the presumably autotetraploid Catolobus genome (2n = 32) arose early after the Catolobus/Capsella divergence. Since its origin, the tetraploid Catolobus genome has undergone chromosomal rediploidization, including a reduction in chromosome number from 2n = 32 to 2n = 30. Diploidization occurred through end-to-end chromosome fusion and other chromosomal rearrangements affecting a total of six of 16 ancestral chromosomes. The hypotetraploid Catolobus cytotype expanded toward its present range, accompanied by some longitudinal genetic differentiation. The sister relationship between Catolobus and Capsella allows comparative studies of tetraploid genomes of contrasting ages and different degrees of genome diploidization.

Genomic and Transcriptional Profiles of Kelch-like (klhl) Gene Family in Polyploid Carassius Complex

Genome duplication supplies raw genetic materials and has been thought to be essential for evolutionary innovation and ecological adaptation. Here, we select Kelch-like (klhl) genes to study the evolution of the duplicated genes in the polyploid Carassius complex, including amphidiploid C. auratus and amphitriploid C. gibelio. Phylogenetic, chromosomal location and read coverage analyses indicate that most of Carassius klhl genes exhibit a 2:1 relationship with zebrafish orthologs and confirm two rounds of polyploidy, an allotetraploidy followed by an autotriploidy, occurred during Carassius evolution. The lineage-specific expansion and biased retention/loss of klhl genes are also found in Carassius. Transcriptome analyses across eight adult tissues and seven embryogenesis stages reveal varied expression dominance and divergence between the two species. The expression of klhls in response to Carassius herpesvirus 2 infection shows different expression changes corresponding to distinct herpesvirus resistances in three C. gibelio gynogenetic clones. Finally, we find that most C. gibelio klhl genes possess three alleles except eight genes that have lost one or two alleles due to genome rearrangement. The allele expression bias is prosperous for Cgklhl genes and varies during embryogenesis owning to the sequential expression manner of the alleles. The current study provides global insights into the genomic and transcriptional evolution of duplicated genes in a given superfamily resulting from multiple rounds of polyploidization.

A non-homogeneous model of chromosome-number evolution to reveal shifts in the transition patterns across the phylogeny

Changes in chromosome numbers, including polyploidy and dysploidy events, play a key role in eukaryote evolution as they could expediate reproductive isolation and have the potential to foster phenotypic diversification. Deciphering the pattern of chromosome-number change within a phylogeny currently relies on probabilistic evolutionary models. All currently available models assume time homogeneity, such that the transition rates are identical throughout the phylogeny. Here, we develop heterogeneous models of chromosome-number evolution that allow multiple transition regimes to operate in distinct parts of the phylogeny. The partition of the phylogeny to distinct transition regimes may be specified by the researcher or, alternatively, identified using a sequential testing approach. Once the number and locations of shifts in the transition pattern are determined, a second search phase identifies regimes with similar transition dynamics, which could indicate on convergent evolution. Using simulations, we study the performance of the developed model to detect shifts in patterns of chromosome-number evolution and demonstrate its applicability by analyzing the evolution of chromosome numbers within the Cyperaceae plant family. The developed model extends the capabilities of probabilistic models of chromosome-number evolution and should be particularly helpful for the analyses of large phylogenies that include multiple distinct subclades.

An in-depth investigation of cryptic taxonomic diversity in the rare endemic mustard Draba maguirei

PREMISE

Previously published evidence suggests that Draba maguirei, a mustard endemic to a few localities in the Bear River, Wellsville, and Wasatch Mountains of northern Utah, may represent a cryptic species complex rather than a single species. Conservation concerns prompted an in-depth systematic study of this taxon and its putative relatives.

METHODS

Sampling most known populations of D. maguirei s.l. (D. maguirei var. maguirei and D. maguirei var. burkei), we integrate data from geography, ecology, morphology, cytogenetics and pollen, enzyme electrophoresis, and the phylogenetic analysis of nuclear internal transcribed spacer sequences to explore potential taxonomic diversity in the species complex.

RESULTS

Draba maguirei var. burkei is shown here to be a distinct species (D. burkei) most closely related to D. globosa, rather than to D. maguirei. Within D. maguirei s.s., the northern (high elevation) and southern (low elevation) population clusters are genetically isolated and morphologically distinguishable, leading to the recognition here of the southern taxon as D. maguirei subsp. stonei.

CONCLUSIONS

Our study reveals that plants traditionally assigned to D. maguirei comprise three genetically divergent lineages (D. burkei and two newly recognized subspecies of D. maguirei), each exhibiting a different chromosome number and occupying a discrete portion of the geographic range. Although previously overlooked and underappreciated taxonomically, the three taxa are morphologically recognizable based on the distribution and types of trichomes present on the leaves, stems, and fruit. Our clarification of the diversity and distribution of these taxa provides an improved framework for conservation efforts.

Disentangling Crocus Series Verni and Its Polyploids

Spring crocuses, the eleven species within Crocus series Verni (Iridaceae), consist of di- and tetraploid cytotypes. Among them is a group of polyploids from southeastern Europe with yet-unclear taxonomic affiliation. Crocuses are generally characterized by complex dysploid chromosome number changes, preventing a clear correlation between these numbers and ploidy levels. To reconstruct the evolutionary history of series Verni and particularly its polyploid lineages associated with C. heuffelianus, we used an approach combining phylogenetic analyses of two chloroplast regions, 14 nuclear single-copy genes plus rDNA spacers, genome-wide genotyping-by-sequencing (GBS) data, and morphometry with ploidy estimations through genome size measurements, analysis of genomic heterozygosity frequencies and co-ancestry, and chromosome number counts. Chromosome numbers varied widely in diploids with 2n = 8, 10, 12, 14, 16, and 28 and tetraploid species or cytotypes with 2n = 16, 18, 20, and 22 chromosomes. Crocus longiflorus, the diploid with the highest chromosome number, possesses the smallest genome (2C = 3.21 pg), while the largest diploid genomes are in a range of 2C = 7-8 pg. Tetraploid genomes have 2C values between 10.88 pg and 12.84 pg. Heterozygosity distribution correlates strongly with genome size classes and allows discernment of di- and tetraploid cytotypes. Our phylogenetic analyses showed that polyploids in the C. heuffelianus group are allotetraploids derived from multiple and partly reciprocal crosses involving different genotypes of diploid C. heuffelianus (2n = 10) and C. vernus (2n = 8). Dysploid karyotype changes after polyploidization resulted in the tetraploid cytotypes with 20 and 22 chromosomes. The multi-data approach we used here for series Verni, combining evidence from nuclear and chloroplast phylogenies, genome sizes, chromosome numbers, and genomic heterozygosity for ploidy estimations, provides a way to disentangle the evolution of plant taxa with complex karyotype changes that can be used for the analysis of other groups within Crocus and beyond. Comparing these results with morphometric analysis results in characters that can discern the different taxa currently subsumed under C. heuffelianus.

Polyploidy on islands - concerted evolution and gene loss amid chromosomal stasis

BACKGROUND AND AIMS

Polyploidy is an important process that often generates genomic diversity within lineages, but it can also cause changes that result in loss of genomic material. Island lineages, while often polyploid, typically show chromosomal stasis but have not been investigated in detail regarding smaller-scale gene loss. Our aim was to investigate post-polyploidization genome dynamics in a chromosomally stable lineage of Malvaceae endemic to New Zealand.

METHODS

We determined chromosome numbers and used fluorescence in situ hybridization to localize 18S and 5S rDNA. Gene sequencing of 18S rDNA, the internal transcribed spacers (ITS) with intervening 5.8S rDNA, and a low-copy nuclear gene, GBSSI-1, was undertaken to determine if gene loss occurred in the New Zealand lineage following polyploidy.

KEY RESULTS

The chromosome number for all species investigated was 2n = 42, with the first published report for the monotypic Australian genus Asterotrichion. The five species investigated all had two 5S rDNA signals localized interstitially on the long arm of one of the largest chromosome pairs. All species, except Plagianthus regius, had two 18S rDNA signals localized proximally on the short arm of one of the smallest chromosome pairs. Plagianthus regius had two additional 18S rDNA signals on a separate chromosome, giving a total of four. Sequencing of nuclear ribosomal 18S rDNA and the ITS cistron indicated loss of historical ribosomal repeats. Phylogenetic analysis of a low-copy nuclear gene, GBSSI-1, indicated that some lineages maintained three copies of the locus, while others have lost one or two copies.

CONCLUSIONS

Although island endemic lineages show chromosomal stasis, with no additional changes in chromosome number, they may undergo smaller-scale processes of gene loss and concerted evolution ultimately leading to further genome restructuring and downsizing.

Phylotranscriptomic analyses reveal multiple whole-genome duplication events, the history of diversification and adaptations in the Araceae

BACKGROUND AND AIMS

The Araceae are one of the most diverse monocot families with numerous morphological and ecological novelties. Plastid and mitochondrial genes have been used to investigate the phylogeny and to interpret shifts in the pollination biology and biogeography of the Araceae. In contrast, the role of whole-genome duplication (WGD) in the evolution of eight subfamilies remains unclear.

METHODS

New transcriptomes or low-depth whole-genome sequences of 65 species were generated through Illumina sequencing. We reconstructed the phylogenetic relationships of Araceae using concatenated and species tree methods, and then estimated the age of major clades using TreePL. We inferred the WGD events by Ks and gene tree methods. We investigated the diversification patterns applying time-dependent and trait-dependent models. The expansions of gene families and functional enrichments were analysed using CAFE and InterProScan.

KEY RESULTS

Gymnostachydoideae was the earliest diverging lineage followed successively by Orontioideae, Lemnoideae and Lasioideae. In turn, they were followed by the clade of 'bisexual climbers' comprised of Pothoideae and Monsteroideae, which was resolved as the sister to the unisexual flowers clade of Zamioculcadoideae and Aroideae. A special WGD event ψ (psi) shared by the True-Araceae clade occurred in the Early Cretaceous. Net diversification rates first declined and then increased through time in the Araceae. The best diversification rate shift along the stem lineage of the True-Araceae clade was detected, and net diversification rates were enhanced following the ψ-WGD. Functional enrichment analyses revealed that some genes, such as those encoding heat shock proteins, glycosyl hydrolase and cytochrome P450, expanded within the True-Araceae clade.

CONCLUSIONS

Our results improve our understanding of aroid phylogeny using the large number of single-/low-copy nuclear genes. In contrast to the Proto-Araceae group and the lemnoid clade adaption to aquatic environments, our analyses of WGD, diversification and functional enrichment indicated that WGD may play a more important role in the evolution of adaptations to tropical, terrestrial environments in the True-Araceae clade. These insights provide us with new resources to interpret the evolution of the Araceae.

Chromoanagenesis in plants: triggers, mechanisms, and potential impact

Chromoanagenesis is a single catastrophic event that involves, in most cases, localized chromosomal shattering and reorganization, resulting in a dramatically restructured chromosome. First discovered in cancer cells, it has since been observed in various other systems, including plants. In this review, we discuss the origin, characteristics, and potential mechanisms underlying chromoanagenesis in plants. We report that multiple processes, including mutagenesis and genetic engineering, can trigger chromoanagenesis via a variety of mechanisms such as micronucleation, breakage-fusion-bridge (BFB) cycles, or chain-like translocations. The resulting rearranged chromosomes can be preserved during subsequent plant growth, and sometimes inherited to the next generation. Because of their high tolerance to genome restructuring, plants offer a unique system for investigating the evolutionary consequences and potential practical applications of chromoanagenesis.

Karyotype evolution of the Asterids insights from the first genome sequences of the family Cornaceae

Cornaceae is a core representative family in Cornales, the earliest branching lineage in the Asterids on the life tree of angiosperms. This family includes the only genus Cornus, a group of ~55 species. These species occur widely in Northern Hemisphere and have been used as resources for horticultural ornaments, medicinal and industrial manufacturing. However, no any genome sequences are available for this family. Here, we reported a chromosome-level genome for Cornus controversa. This was generated using high-fidelity plus Hi-C sequencing, and totally ~771.80 Mb assembled sequences and 39,886 protein-coding genes were obtained. We provided evidence for a whole-genome duplication event (WGD) unique to C. controversa. The evolutionary features of this genome indicated that the expanded and unique genes might have contributed to response to stress, stimulus and defense. By using chromosome-level syntenic blocks shared between eight living genomes, we found high degrees of genomic diversification from the ancestral core-eudicot genome to the present-day genomes, suggesting an important role of WGD in genomic plasticity that leads to speciation and diversification. These results provide foundational insights on the evolutionary history of Cornaceae, as well as on the Asterids diversification.

The genome sequence of Hirschfeldia incana, a new Brassicaceae model to improve photosynthetic light-use efficiency

Photosynthesis is a key process in sustaining plant and human life. Improving the photosynthetic capacity of agricultural crops is an attractive means to increase their yields. While the core mechanisms of photosynthesis are highly conserved in C₃ plants, these mechanisms are very flexible, allowing considerable diversity in photosynthetic properties. Among this diversity is the maintenance of high photosynthetic light-use efficiency at high irradiance as identified in a small number of exceptional C₃ species. Hirschfeldia incana, a member of the Brassicaceae family, is such an exceptional species, and because it is easy to grow, it is an excellent model for studying the genetic and physiological basis of this trait. Here, we present a reference genome of H. incana and confirm its high photosynthetic light-use efficiency. While H. incana has the highest photosynthetic rates found so far in the Brassicaceae, the light-saturated assimilation rates of closely related Brassica rapa and Brassica nigra are also high. The H. incana genome has extensively diversified from that of B. rapa and B. nigra through large chromosomal rearrangements, species-specific transposon activity, and differential retention of duplicated genes. Duplicated genes in H. incana, B. rapa, and B. nigra that are involved in photosynthesis and/or photoprotection show a positive correlation between copy number and gene expression, providing leads into the mechanisms underlying the high photosynthetic efficiency of these species. Our work demonstrates that the H. incana genome serves as a valuable resource for studying the evolution of high photosynthetic light-use efficiency and enhancing photosynthetic rates in crop species.

Comparative cytogenomics reveals genome reshuffling and centromere repositioning in the legume tribe Phaseoleae

The tribe Phaseoleae includes several legume crops with assembled genomes. Comparative genomic studies have evidenced the preservation of large genomic blocks among legumes, although chromosome dynamics during Phaseoleae evolution has not been investigated. We conducted a comparative genomic analysis to define an informative genomic block (GB) system and to reconstruct the ancestral Phaseoleae karyotype (APK). We identified GBs based on the orthologous genes between Phaseolus vulgaris and Vigna unguiculata and searched for GBs in different genomes of the Phaseolinae (P. lunatus) and Glycininae (Amphicarpaea edgeworthii) subtribes and Spatholobus suberectus (sister to Phaseolinae and Glycininae), using Medicago truncatula as the outgroup. We also used oligo-FISH probes of two P. vulgaris chromosomes to paint the orthologous chromosomes of two non-sequenced Phaseolinae species. We inferred the APK as having n = 11 and 19 GBs (A to S), hypothesizing five chromosome fusions that reduced the ancestral legume karyotype to n = 11. We identified the rearrangements among the APK and the subtribes and species, with extensive centromere repositioning in Phaseolus. We also reconstructed the chromosome number reduction in S. suberectus. The development of the GB system and the proposed APK provide useful approaches for future comparative genomic analyses of legume species.

Potential allopolyploid origin of Ericales revealed with gene-tree reconciliation

Few incidents of ancient allopolyploidization (polyploidization by hybridization or merging diverged genomes) were previously revealed, although there is significant evidence for the accumulation of whole genome duplications (WGD) in plants. Here, we focused on Ericales, one of the largest and most diverse angiosperm orders with significant ornamental and economic value. Through integrating 24 high-quality whole genome data selected from ~ 200 Superasterids genomes/species and an algorithm of topology-based gene-tree reconciliation, we explored the evolutionary history of in Ericales with ancient complex. We unraveled the allopolyploid origin of Ericales and detected extensive lineage-specific gene loss following the polyploidization. Our study provided a new hypothesis regarding the origin of Ericales and revealed an instructive perspective of gene loss as a pervasive source of genetic variation and adaptive phenotypic diversity in Ericales.

An atlas of fish genome evolution reveals delayed rediploidization following the teleost whole-genome duplication

Teleost fishes are ancient tetraploids descended from an ancestral whole-genome duplication that may have contributed to the impressive diversification of this clade. Whole-genome duplications can occur via self-doubling (autopolyploidy) or via hybridization between different species (allopolyploidy). The mode of tetraploidization conditions evolutionary processes by which duplicated genomes return to diploid meiotic pairing, and subsequent genetic divergence of duplicated genes (cytological and genetic rediploidization). How teleosts became tetraploid remains unresolved, leaving a fundamental gap in the interpretation of their functional evolution. As a result of the whole-genome duplication, identifying orthologous and paralogous genomic regions across teleosts is challenging, hindering genome-wide investigations into their polyploid history. Here, we combine tailored gene phylogeny methodology together with a state-of-the-art ancestral karyotype reconstruction to establish the first high-resolution comparative atlas of paleopolyploid regions across 74 teleost genomes. We then leverage this atlas to investigate how rediploidization occurred in teleosts at the genome-wide level. We uncover that some duplicated regions maintained tetraploidy for more than 60 million years, with three chromosome pairs diverging genetically only after the separation of major teleost families. This evidence suggests that the teleost ancestor was an autopolyploid. Further, we find evidence for biased gene retention along several duplicated chromosomes, contradicting current paradigms that asymmetrical evolution is specific to allopolyploids. Altogether, our results offer novel insights into genome evolutionary dynamics following ancient polyploidizations in vertebrates.

The Chromosome Number and rDNA Loci Evolution in Onobrychis (Fabaceae)

The evolution of chromosome number and ribosomal DNA (rDNA) loci number and localisation were studied in Onobrychis Mill. Diploid and tetraploid species, as well as two basic chromosome numbers, x = 7 and x = 8, were observed among analysed taxa. The chromosomal distribution of rDNA loci was presented here for the first time using fluorescence in situ hybridisation (FISH) with 5S and 35S rDNA probes. Onobrychis species showed a high polymorphism in the number and localisation of rDNA loci among diploids, whereas the rDNA loci pattern was very similar in polyploids. Phylogenetic relationships among the species, inferred from nrITS sequences, were used as a framework to reconstruct the patterns of basic chromosome number and rDNA loci evolution. Analysis of the evolution of the basic chromosome numbers allowed the inference of x = 8 as the ancestral number and the descending dysploidy and polyploidisation as the major mechanisms of the chromosome number evolution. Analyses of chromosomal patterns of rRNA gene loci in a phylogenetic context resulted in the reconstruction of one locus of 5S rDNA and one locus of 35S rDNA in the interstitial chromosomal position as the ancestral state in this genus.

Comparative Analyses of Complete Chloroplast Genomes and Karyotypes of Allotetraploid Iris koreana and Its Putative Diploid Parental Species (Iris Series Chinenses, Iridaceae)

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.

Dynamic genome evolution in a model fern

The large size and complexity of most fern genomes have hampered efforts to elucidate fundamental aspects of fern biology and land plant evolution through genome-enabled research. Here we present a chromosomal genome assembly and associated methylome, transcriptome and metabolome analyses for the model fern species Ceratopteris richardii. The assembly reveals a history of remarkably dynamic genome evolution including rapid changes in genome content and structure following the most recent whole-genome duplication approximately 60 million years ago. These changes include massive gene loss, rampant tandem duplications and multiple horizontal gene transfers from bacteria, contributing to the diversification of defence-related gene families. The insertion of transposable elements into introns has led to the large size of the Ceratopteris genome and to exceptionally long genes relative to other plants. Gene family analyses indicate that genes directing seed development were co-opted from those controlling the development of fern sporangia, providing insights into seed plant evolution. Our findings and annotated genome assembly extend the utility of Ceratopteris as a model for investigating and teaching plant biology.

Analysis of flavonol regulator evolution in the Brassicaceae reveals MYB12, MYB111 and MYB21 duplications and MYB11 and MYB24 gene loss

BACKGROUND

Flavonols are the largest subgroup of flavonoids, possessing multiple functions in plants including protection against ultraviolet radiation, antimicrobial activities, and flower pigmentation together with anthocyanins. They are of agronomical and economical importance because the major off-taste component in rapeseed protein isolates is a flavonol derivative, which limits rapeseed protein use for human consumption. Flavonol production in Arabidopsis thaliana is mainly regulated by the subgroup 7 (SG7) R2R3-MYB transcription factors MYB11, MYB12, and MYB111. Recently, the SG19 MYBs MYB21, MYB24, and MYB57 were shown to regulate flavonol accumulation in pollen and stamens. The members of each subgroup are closely related, showing gene redundancy and tissue-specific expression in A. thaliana. However, the evolution of these flavonol regulators inside the Brassicaceae, especially inside the Brassiceae, which include the rapeseed crop species, is not fully understood.

RESULTS

We studied the SG7 and SG19 MYBs in 44 species, including 31 species of the Brassicaceae, by phylogenetic analyses followed by synteny and gene expression analyses. Thereby we identified a deep MYB12 and MYB111 duplication inside the Brassicaceae, which likely occurred before the divergence of Brassiceae and Thelypodieae. These duplications of SG7 members were followed by the loss of MYB11 after the divergence of Eruca vesicaria from the remaining Brassiceae species. Similarly, MYB21 experienced duplication before the emergence of the Brassiceae tribe, where the gene loss of MYB24 is also proposed to have happened. The members of each subgroup revealed frequent overlapping spatio-temporal expression patterns in the Brassiceae member B. napus, which are assumed to compensate for the loss of MYB11 and MYB24 in the analysed tissues.

CONCLUSIONS

We identified a duplication of MYB12, MYB111, and MYB21 inside the Brassicaceae and MYB11 and MYB24 gene loss inside the tribe Brassiceae. We propose that polyploidization events have shaped the evolution of the flavonol regulators in the Brassicaceae, especially in the Brassiceae.