Chromosome size matters: genome evolution in the cyperid clade

The smallest angiosperm genomes may be the price for effective traps of bladderworts

BACKGROUND

Species of the carnivorous family Lentibulariaceae exhibit the smallest genomes in flowering plants. We explored the hypothesis that their minute genomes result from the unique mitochondrial cytochrome c oxidase (COX) mutation. The mutation may boost mitochondrial efficiency, which is especially useful for suction-bladder traps of Utricularia, but also increase DNA-damaging reactive oxygen species, leading to genome shrinkage through deletion-biased DNA repair. We aimed to explore the impact of this mutation on genome size, providing insights into genetic mutation roles in plant genome evolution under environmental pressures.

METHODS

We compiled and measured genome and mean chromosome sizes for 127 and 67 species, respectively, representing all three genera (Genlisea, Pinguicula and Utricularia) of Lentibulariaceae. We also isolated and analysed COX sequences to detect the mutation. Through phylogenetic regressions and Ornstein-Uhlenbeck models of trait evolution, we assessed the impact of the COX mutation on the genome and chromosome sizes across the family.

RESULTS

Our findings reveal significant correlations between the COX mutation and smaller genome and chromosome sizes. Specifically, species carrying the ancestral COX sequence exhibited larger genomes and chromosomes than those with the novel mutation. This evidence supports the notion that the COX mutation contributes to genome downsizing, with statistical analyses confirming a directional evolution towards smaller genomes in species harbouring these mutations.

CONCLUSIONS

Our study confirms that the COX mutation in Lentibulariaceae is associated with genome downsizing, probably driven by increased reactive oxygen species production and subsequent DNA damage requiring deletion-biased repair mechanisms. While boosting mitochondrial energy output, this genetic mutation compromises genome integrity and may potentially affect recombination rates, illustrating a complex trade-off between evolutionary advantages and disadvantages. Our results highlight the intricate processes by which genetic mutations and environmental pressures shape genome size evolution in carnivorous plants.

Centromere drive may propel the evolution of chromosome and genome size in plants

BACKGROUND

Genome size is influenced by natural selection and genetic drift acting on variations from polyploidy and repetitive DNA sequences. We hypothesized that centromere drive, where centromeres compete for inclusion in the functional gamete during meiosis, may also affect genome and chromosome size. This competition occurs in asymmetric meiosis, where only one of the four meiotic products becomes a gamete. If centromere drive influences chromosome size evolution, it may also impact post-polyploid diploidization, where a polyploid genome is restructured to function more like a diploid through chromosomal rearrangements, including fusions. We tested if plant lineages with asymmetric meiosis exhibit faster chromosome size evolution compared to those with only symmetric meiosis, which lack centromere drive as all four meiotic products become gametes. We also examined if positive selection on centromeric histone H3 (CENH3), a protein that can suppress centromere drive, is more frequent in these asymmetric lineages.

METHODS

We analysed plant groups with different meiotic modes: asymmetric in gymnosperms and angiosperms, and symmetric in bryophytes, lycophytes and ferns. We selected species based on available CENH3 gene sequences and chromosome size data. Using Ornstein-Uhlenbeck evolutionary models and phylogenetic regressions, we assessed the rates of chromosome size evolution and the frequency of positive selection on CENH3 in these clades.

RESULTS

Our analyses showed that clades with asymmetric meiosis have a higher frequency of positive selection on CENH3 and increased rates of chromosome size evolution compared to symmetric clades.

CONCLUSIONS

Our findings support the hypothesis that centromere drive accelerates chromosome and genome size evolution, potentially also influencing the process of post-polyploid diploidization. We propose a model which in a single framework helps explain the stability of chromosome size in symmetric lineages (bryophytes, lycophytes and ferns) and its variability in asymmetric lineages (gymnosperms and angiosperms), providing a foundation for future research in plant genome evolution.

Macroevolutionary inference of complex modes of chromosomal speciation in a cosmopolitan plant lineage

The effects of single chromosome number change-dysploidy - mediating diversification remain poorly understood. Dysploidy modifies recombination rates, linkage, or reproductive isolation, especially for one-fifth of all eukaryote lineages with holocentric chromosomes. Dysploidy effects on diversification have not been estimated because modeling chromosome numbers linked to diversification with heterogeneity along phylogenies is quantitatively challenging. We propose a new state-dependent diversification model of chromosome evolution that links diversification rates to dysploidy rates considering heterogeneity and differentiates between anagenetic and cladogenetic changes. We apply this model to Carex (Cyperaceae), a cosmopolitan flowering plant clade with holocentric chromosomes. We recover two distinct modes of chromosomal evolution and speciation in Carex. In one diversification mode, dysploidy occurs frequently and drives faster diversification rates. In the other mode, dysploidy is rare, and diversification is driven by hidden, unmeasured factors. When we use a model that excludes hidden states, we mistakenly infer a strong, uniformly positive effect of dysploidy on diversification, showing that standard models may lead to confident but incorrect conclusions about diversification. This study demonstrates that dysploidy can have a significant role in speciation in a large plant clade despite the presence of other unmeasured factors that simultaneously affect diversification.

The evolutionary landscape of prokaryotic chromosome/plasmid balance

The balance between chromosomal and plasmid DNAs determines the genomic plasticity of prokaryotes. Natural selections, acting on the level of organisms or plasmids, shape the abundances of plasmid DNAs in prokaryotic genomes. Despite the importance of plasmids in health and engineering, there have been rare systematic attempts to quantitatively model and predict the determinants underlying the strength of different selection forces. Here, we develop a metabolic flux model that describes the intracellular resource competition between chromosomal and plasmid-encoded reactions. By coarse graining, this model predicts a landscape of natural selections on chromosome/plasmid balance, which is featured by the tradeoff between phenotypic and non-phenotypic selection pressures. This landscape is further validated by the observed pattern of plasmid distributions in the vast collection of prokaryotic genomes retrieved from the NCBI database. Our results establish a universal paradigm to understand the prokaryotic chromosome/plasmid interplay and provide insights into the evolutionary origin of plasmid diversity.

Genome size variation in Cape schoenoid sedges (Schoeneae) and its ecophysiological consequences

PREMISE

Increases in genome size in plants-often associated with larger, low-density stomata and greater water-use efficiency (WUE)-could affect plant ecophysiological and hydraulic function. Variation in plant genome size is often due to polyploidy, having occurred repeatedly in the austral sedge genus Schoenus in the Cape Floristic Region (CFR), while species in the other major schoenoid genus in the region, Tetraria, have smaller genomes. Comparing these genera is useful as they co-occur at the landscape level, under broadly similar bioclimatic conditions. We hypothesized that CFR Schoenus have greater WUE, with lower maximum stomatal conductance (gwmax) imposed by larger, less-dense stomata.

METHODS

We investigated relationships between genome size and stomatal parameters in a phylogenetic context, reconstructing a phylogeny of CFR-occurring Schoeneae (Cyperaceae). Species' stomatal and functional traits were measured from field-collected and herbarium specimens. Carbon stable isotopes were used as an index of WUE. Genome size was derived from flow-cytometric measurements of leafy shoots.

RESULTS

Evolutionary regressions demonstrated that stomatal size and density covary with genome size, positively and negatively, respectively, with genome size explaining 72-75% of the variation in stomatal size. Larger-genomed species had lower gwmax and C:N ratios, particularly in culms.

CONCLUSIONS

We interpret differences in vegetative physiology between the genera as evidence of more-conservative strategies in CFR Schoenus compared to the more-acquisitive Tetraria. Because Schoenus have smaller, reduced leaves, they likely rely more on culm photosynthesis than Tetraria. Across the CFR Schoeneae, ecophysiology correlates with genome size, but confounding sources of trait variation limit inferences about causal relationships between traits.

Chromosome-Scale Genome Assembly for Soft-Stem Bulrush (Schoenoplectus tabernaemontani) Confirms a Clade-Specific Whole-Genome Duplication in Cyperaceae

Schoenoplectus tabernaemontani (C. C. Gmelin) Palla is a typical macrophyte in diverse wetland ecosystems. This species holds great potential in decontamination applications and carbon sequestration. Previous studies have shown that this species may have experienced recent polyploidization. This would make S. tabernaemontani a unique model to study the processes and consequences of whole-genome duplications in the context of the well-documented holocentric chromosomes and dysploidy events in Cyperaceae. However, the inference was not completely solid because it lacked homology information that is essential to ascertain polyploidy. We present here the first chromosome-level genome assembly for S. tabernaemontani. By combining Oxford Nanopore Technologies (ONT) long reads and Illumina short reads, plus chromatin conformation via the Hi-C method, we assembled a genome spanning 507.96 Mb, with 99.43% of Hi-C data accurately mapped to the assembly. The assembly contig N50 value was 3.62 Mb. The overall BUSCO score was 94.40%. About 68.94% of the genome was comprised of repetitive elements. A total of 36,994 protein-coding genes were predicted and annotated. Long terminal repeat retrotransposons accounted for ∼26.99% of the genome, surpassing the content observed in most sequenced Cyperid genomes. Our well-supported haploid assembly comprised 21 pseudochromosomes, each harboring putative holocentric centromeres. Our findings corroborated a karyotype of 2n = 2X = 42. We also confirmed a recent whole-genome duplication occurring after the divergence between Schoenoplecteae and Bolboschoeneae. Our genome assembly expands the scope of sequenced genomes within the Cyperaceae family, encompassing the fifth genus. It also provides research resources on Cyperid evolution and wetland conservation.

Founder events and subsequent genetic bottlenecks underlie karyotype evolution in the Ibero-North African endemic Carex helodes

BACKGROUND AND AIMS

Despite chromosomal evolution being one of the major drivers of diversification in plants, we do not yet have a clear view of how new chromosome rearrangements become fixed within populations, which is a crucial step forward for understanding chromosomal speciation.

METHODS

In this study, we test the role of genetic drift in the establishment of new chromosomal variants in the context of hybrid dysfunction models of chromosomal speciation. We genotyped 178 individuals from seven populations (plus 25 seeds from one population) across the geographical range of Carex helodes (Cyperaceae). We also characterized karyotype geographical patterns of the species across its distribution range. For one of the populations, we performed a detailed study of the fine-scale, local spatial distribution of its individuals and their genotypes and karyotypes.

KEY RESULTS

Synergistically, phylogeographical and karyotypic evidence revealed two main genetic groups: southwestern Iberian Peninsula vs. northwestern African populations; and within Europe our results suggest a west-to-east expansion with signals of genetic bottlenecks. Additionally, we inferred a pattern of descending dysploidy, plausibly as a result of a west-to-east process of post-glacial colonization in Europe.

CONCLUSIONS

Our results give experimental support to the role of geographical isolation, drift and inbreeding in the establishment of new karyotypes, which is key in the speciation models of hybrid dysfunction.

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.

Global analysis of Poales diversification - parallel evolution in space and time into open and closed habitats

Poales are one of the most species-rich, ecologically and economically important orders of plants and often characterise open habitats, enabled by unique suites of traits. We test six hypotheses regarding the evolution and assembly of Poales in open and closed habitats throughout the world, and examine whether diversification patterns demonstrate parallel evolution. We sampled 42% of Poales species and obtained taxonomic and biogeographic data from the World Checklist of Vascular Plants database, which was combined with open/closed habitat data scored by taxonomic experts. A dated supertree of Poales was constructed. We integrated spatial phylogenetics with regionalisation analyses, historical biogeography and ancestral state estimations. Diversification in Poales and assembly of open and closed habitats result from dynamic evolutionary processes that vary across lineages, time and space, most prominently in tropical and southern latitudes. Our results reveal parallel and recurrent patterns of habitat and trait transitions in the species-rich families Poaceae and Cyperaceae. Smaller families display unique and often divergent evolutionary trajectories. The Poales have achieved global dominance via parallel evolution in open habitats, with notable, spatially and phylogenetically restricted divergences into strictly closed habitats.

Chromosome-Scale Genome Assembly for Clubrush (Bolboschoenus planiculmis) Indicates a Karyotype with High Chromosome Number and Heterogeneous Centromere Distribution

Bolboschoenus planiculmis (F.Schmidt) T.V.Egorova is a typical wetland plant in the species-rich Cyperaceae family. This species contributes prominently to carbon dynamics and trophic integration in wetland ecosystems. Previous studies have reported that the chromosomes of B. planiculmis are holocentric; i.e. they have kinetic activity along their entire length and carry multiple centromeres. This feature was suggested to lead to a rapid genome evolution through chromosomal fissions and fusions and participate to the diversification and ecological success of the Bolboschoenus genus. However, the specific mechanism remains uncertain, partly due to the scarcity of genetic information on Bolboschoenus. We present here the first chromosome-level genome assembly for B. planiculmis. Through the integration of high-quality long-read and short-read data, together with chromatin conformation using Hi-C technology, the ultimate genome assembly was 238.01 Mb with a contig N50 value of 3.61 Mb. Repetitive elements constituted 37.04% of the genome, and 18,760 protein-coding genes were predicted. The low proportion of long terminal repeat retrotransposons (∼9.62%) was similar to that reported for other Cyperaceae species. The Ks (synonymous substitutions per synonymous site) distribution suggested no recent large-scale genome duplication in this genome. The haploid assembly contained a large number of 54 pseudochromosomes with a small mean size of 4.10 Mb, covering most of the karyotype. The results of centromere detection support that not all the chromosomes in B. planiculmis have multiple centromeres, indicating more efforts are needed to fully reveal the specific style of holocentricity in cyperids and its evolutionary significance.

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.