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.

Molecular characterization of oleosin genes in Cyperus esculentus, a Cyperaceae plant producing oil in underground tubers.

KEY MESSAGE: CeOLE genes exhibit a tuber-predominant expression pattern and their mRNA/protein abundances are positively correlated with oil accumulation during tuber development. Overexpression could significantly increase the oil content of tobacco leaves. Oleosins (OLEs) are abundant structural proteins of lipid droplets (LDs) that function in LD formation and stabilization in seeds of oil crops. However, little information is available on their roles in vegetative tissues. In this study, we present the first genome-wide characterization of the oleosin family in tigernut (Cyperus esculentus L., Cyperaceae), a rare example accumulating high amounts of oil in underground tubers. Six members identified represent three previously defined clades (i.e. U, SL and SH) or six out of seven orthogroups (i.e. U, SL1, SL2, and SH1-3) proposed in this study. Comparative genomics analysis reveals that lineage-specific expansion of Clades SL and SH was contributed by whole-genome duplication and dispersed duplication, respectively. Moreover, presence of SL2 and SH3 in Juncus effuses implies their appearance sometime before Cyperaceae-Juncaceae divergence, whereas SH2 appears to be Cyperaceae specific. Expression analysis showed that CeOLE genes exhibit a tuber-predominant expression pattern and transcript levels are considerably more abundant than homologs in the close relative Cyperus rotundus. Moreover, CeOLE mRNA and protein abundances were shown to positively correlate with oil accumulation during tuber development. Additionally, two dominant isoforms (i.e. CeOLE2 and -5) were shown to locate in LDs as well as the endoplasmic reticulum of tobacco (Nicotiana benthamiana) leaves, and are more likely to function in homo and heteromultimers. Furthermore, overexpression of CeOLE2 and -5 in tobacco leaves could significantly increase the oil content, supporting their roles in oil accumulation. These findings provide insights into lineage-specific family evolution and putative roles of CeOLE genes in oil accumulation of vegetative tissues, which facilitate further genetic improvement for tigernut.

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.

Complex patterns of ploidy in a holocentric plant clade (Schoenus, Cyperaceae) in the Cape biodiversity hotspot.

BACKGROUND AND AIMS: It is unclear how widespread polyploidy is throughout the largest holocentric plant family - the Cyperaceae. Because of the prevalence of chromosomal fusions and fissions, which affect chromosome number but not genome size, it can be impossible to distinguish if individual plants are polyploids in holocentric lineages based on chromosome count data alone. Furthermore, it is unclear how differences in genome size and ploidy levels relate to environmental correlates within holocentric lineages, such as the Cyperaceae. METHODS: We focus our analyses on tribe Schoeneae, and more specifically the southern African clade of Schoenus. We examine broad-scale patterns of genome size evolution in tribe Schoeneae and focus more intensely on determining the prevalence of polyploidy across the southern African Schoenus by inferring ploidy level with the program ChromEvol, as well as interpreting chromosome number and genome size data. We further investigate whether there are relationships between genome size/ploidy level and environmental variables across the nutrient-poor and summer-arid Cape biodiversity hotspot. KEY RESULTS: Our results show a large increase in genome size, but not chromosome number, within Schoenus compared to other species in tribe Schoeneae. Across Schoenus, there is a positive relationship between chromosome number and genome size, and our results suggest that polyploidy is a relatively common process throughout the southern African Schoenus. At the regional scale of the Cape, we show that polyploids are more often associated with drier locations that have more variation in precipitation between dry and wet months, but these results are sensitive to the classification of ploidy level. CONCLUSIONS: Polyploidy is relatively common in the southern African Schoenus, where a positive relationship is observed between chromosome number and genome size. Thus, there may be a high incidence of polyploidy in holocentric plants, whose cell division properties differ from monocentrics.

RAD sequencing resolves the phylogeny, taxonomy and biogeography of Trichophoreae despite a recent rapid radiation (Cyperaceae).

Trichophoreae is a nearly cosmopolitan Cyperaceae tribe that contains ~17 species displaying striking variation in size, inflorescence complexity, and perianth morphology. Although morphologically distinct, the status of its three genera (Cypringlea, Oreobolopsis and Trichophorum) is controversial because recent phylogenetic studies have suggested they might not be reciprocally monophyletic. However, previous analyses have shown conflicting topologies and consistently poor support due to an initial rapid diversification of the tribe. We analysed restriction-site associated DNA sequencing (RADseq) data from nearly all species of the clade, combined with five Sanger-based markers (matK, ndhF, rps16, ETS-1f, ITS) sampled extensively within species. This approach allowed us to resolve deep and shallow relationships within Trichophoreae for the first time, despite an anomaly zone spanning several successive short branches that produced considerable gene tree incongruence. Analyses reveal a primary phylogenetic split of the tribe into two clades roughly corresponding to an East Asian-North American disjunction that dates back to the mid-Miocene, with both clades comprised of a mixture of reduced unispicate and larger taxa with highly compound inflorescences. Morphological characters traditionally used in the circumscription of Trichophoreae genera are shown to be homoplasious. Several of these characters correlate best with climatic conditions, with the most reduced species occurring in open habitats at high latitudes and altitudes. Close relatives with highly compound inflorescences are found in temperate or subtropical forest understories. Cypringlea and Oreobolopsis are deeply nested within Trichophorum, and we merge all three genera into a more broadly circumscribed Trichophorum. We also show that Scirpus filipes is another previously unrecognized East Asian species of Trichophorum with highly compound inflorescences.

Inferring hypothesis-based transitions in clade-specific models of chromosome number evolution in sedges (Cyperaceae).

Large-scale changes in chromosome number have been associated with diversification rate shifts in many lineages of plants. For instance, several ancient rounds of polyploidization events have been inferred to promote genomic differentiation and/or isolation and, consequently, angiosperm diversification. Dysploidy, although less studied, has been suggested to also play an important role in angiosperm diversification. In this article, we aim to elucidate the role of chromosomal rearrangements on lineage diversification by analyzing a new comprehensive sedge (Cyperaceae) phylogenetic tree. Our null hypothesis is that the mode and tempo of chromosome evolution are to be homogeneous across the complete phylogeny. In order to discern patterns of diversification shifts and chromosome number changes within the family tree, we tested clade-specific chromosome evolution models for several subtrees according to previously reported increments of diversification rates. Results show that a complex, heterogeneous model composed of different clade-specific chromosome evolution transitions are significantly supported against the null hypothesis of a model with no chromosome number model transition events along the phylogeny. This could suggest a link between diversification and changes in chromosome number evolution although other possibilities are not discarded. Our methodological approach may allow identifying different patterns of chromosome evolution, as found for Cyperaceae, for other lineages at different evolutionary levels.

Are holocentrics doomed to change? Limited chromosome number variation in Rhynchospora Vahl (Cyperaceae).

Karyotype evolution in species with non-localised centromeres (holocentric chromosomes) is usually very dynamic and associated with recurrent fission and fusion (also termed agmatoploidy/symploidy) events. In Rhynchospora (Cyperaceae), one of the most species-rich sedge genera, all analysed species have holocentric chromosomes and their numbers range from 2n = 4 to 2n = 84. Agmatoploidy/symploidy and polyploidy were suggested as the main processes in the reshuffling of Rhynchospora karyotypes, although testing different scenarios of chromosome number evolution in a phylogenetic framework has not been attempted until now. Here, we used maximum likelihood and model-based analyses, in combination with genome size estimation and ribosomal DNA distribution, to understand chromosome evolution in Rhynchospora. Overall, chromosome number variation showed a significant phylogenetic signal and the majority of the lineages maintained a karyotype of 2n = 10 (~48% of the species), the most likely candidate for the ancestral number of the genus. Higher and lower chromosome numbers were restricted to specific clades, whilst polyploidy and/or fusion/fission events were present in specific branches. Variation in genome size and ribosomal DNA site number showed no correlation with ploidy level or chromosome number. Although different mechanisms of karyotype evolution (polyploidy, fusion and fission) seem to be acting in distinct lineages, the degree of chromosome variation and the main mechanisms involved are comparable to those found in some monocentric genera and lower than expected for a holocentric genus.

Radiation and repeated transoceanic dispersal of Schoeneae (Cyperaceae) through the southern hemisphere.

PREMISE OF THE STUDY: The broad austral distribution of Schoeneae is almost certainly a product of long-distance dispersal. Owing to the inadequacies of existing phylogenetic data and a lack of rigorous biogeographic analysis, relationships within the tribe remain poorly resolved and its pattern of radiation and dispersal uncertain. We employed an expanded sampling of taxa and markers and a rigorous analytic approach to address these limitations. We evaluated the roles of geography and ecology in stimulating the initial radiation of the group and its subsequent dispersal across the southern hemisphere. METHODS: A dated tree was reconstructed using reversible-jump Markov chain Monte Carlo (MCMC) with a polytomy prior and molecular dating, applied to data from two nuclear and three cpDNA regions. Ancestral areas and habitats were inferred using dispersal-extinction-cladogenesis models. KEY RESULTS: Schoeneae originated in Australia in the Paleocene. The existence of a "hard" polytomy at the base of the clade reflects the rapid divergence of six principal lineages ca. 50 Ma, within Australia. From this ancestral area, Schoeneae have traversed the austral oceans with remarkable frequency, a total of 29 distinct dispersal events being reported here. Dispersal rates between landmasses are not explicable in terms of the geographical distances separating them. Transoceanic dispersal generally involved habitat stasis. CONCLUSIONS: Although the role of dispersal in explaining global distribution patterns is now widely accepted, the apparent ease with which such dispersal may occur has perhaps been under-appreciated. In Schoeneae, transoceanic dispersal has been remarkably frequent, with ecological opportunity, rather than geography, being most important in dictating dispersal patterns.

Characterisation of ALS genes in the polyploid species Schoenoplectus mucronatus and implications for resistance management.

BACKGROUND: The polyploid weed Schoenoplectus mucronatus (L.) Palla has evolved target-site resistance to ALS-inhibiting herbicides in Italian rice crops. Molecular and genetic characterisation of the resistance mechanism is relevant to the evolution and management of herbicide resistance. The authors aimed (a) to study the organisation of the target-site loci in two field-selected S. mucronatus populations with different cross-resistance patterns, (b) to identify the mutations endowing resistance to ALS inhibitors and determine the role of these mutations by using transgenesis and (c) to analyse the implications for the management of the S. mucronatus populations. RESULTS: Two complete ALS genes (ALS1 and ALS2) having an intron and a third partial intronless ALS gene (ALS3) were identified. The presence of multiple ALS genes was confirmed by Southern blot analyses, and ALS loci were characterised by examining cytosine methylation. In S. mucronatus leaves, the transcripts of ALS1, ALS2 and ALS3 were detected. Two mutations endowing resistance (Pro(197) to His and Trp(574) to Leu) were found in both resistant populations, but at different frequencies. Tobacco plants transformed with the two resistant alleles indicated that the Pro(197)-to-His substitution conferred resistance to SU and TP herbicides, while the allele with the Trp(574)-to-Leu substitution conferred cross-resistance to SU, TP, IMI and PTB herbicides. CONCLUSION: Schoenoplectus mucronatus has multiple ALS genes characterised by methylated sites that can influence the expression profile. The two mutated alleles proved to be responsible for ALS resistance. At population level, the resistance pattern depends on the frequency of various resistant genotypes, and this influences the efficacy of various ALS-inhibiting herbicides.

A phylogeny of the schoenoid sedges (Cyperaceae: Schoeneae) based on plastid DNA sequences, with special reference to the genera found in Africa.

Despite its large size (about 700 species), the australy-centred sedge tribe Schoeneae has received little explicit phylogenetic study, especially using molecular data. As a result, generic relationships are poorly understood, and even the monophyly of the tribe is open to question. In this study, plastid DNA sequences (rbcL, trnL-trnF, and rps16) drawn from a broad array of Schoeneae are analysed using Bayesian and parsimony-based approaches to infer a framework phylogeny for the tribe. Both analytical methods broadly support the monophyly of Schoeneae, Bayesian methods doing so with good support. Within the schoenoid clade, there is strong support for a series of monophyletic generic groupings whose interrelationships are unclear. These lineages form a large polytomy at the base of Schoeneae that may be indicative of past radiation, probably following the fragmentation of Gondwana. Most of these lineages contain both African and non-African members, suggesting a history of intercontinental dispersal. The results of this study clearly identify the relationships of the African-endemic schoenoid genera and demonstrate that the African-Australasian genus Tetraria, like Costularia, is polyphyletic. This pattern is morphologically consistent and suggests that these genera require realignment.