Fine-scale empirical data on niche divergence and homeolog expression patterns in an allopolyploid and its diploid progenitor species

Phylotranscriptomics provides a treasure trove of flood-tolerance mechanisms in the Cardamineae tribe

Flooding events are highly detrimental to most terrestrial plant species. However, there is an impressive diversity of plant species that thrive in flood-prone regions and represent a treasure trove of unexplored flood-resilience mechanisms. Here we surveyed a panel of four species from the Cardamineae tribe representing a broad tolerance range. This included the flood-tolerant Cardamine pratensis, Rorippa sylvestris and Rorippa palustris and the flood-sensitive species Cardamine hirsuta. All four species displayed a quiescent strategy, evidenced by the repression of shoot growth underwater. Comparative transcriptomics analyses between the four species and the sensitive model species Arabidopsis thaliana were facilitated via de novo transcriptome assembly and identification of 16 902 universal orthogroups at a high resolution. Our results suggest that tolerance likely evolved separately in the Cardamine and Rorippa species. While the Rorippa response was marked by a strong downregulation of cell-cycle genes, Cardamine minimized overall transcriptional regulation. However, a weak starvation response was a universal trait of tolerant species, potentially achieved in multiple ways. It could result from a strong decline in cell-cycle activity, but is also intertwined with autophagy, senescence, day-time photosynthesis and night-time fermentation capacity. Our data set provides a rich source to study adaptational mechanisms of flooding tolerance.

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

BACKGROUND AND AIMS

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

METHODS

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

KEY RESULTS

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

CONCLUSIONS

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

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.

Novel Tri-Segmented Rhabdoviruses: A Data Mining Expedition Unveils the Cryptic Diversity of Cytorhabdoviruses

Cytorhabdoviruses (genus Cytorhabdovirus, family Rhabdoviridae) are plant-infecting viruses with enveloped, bacilliform virions. Established members of the genus Cytorhabdovirus have unsegmented single-stranded negative-sense RNA genomes (ca. 10-16 kb) which encode four to ten proteins. Here, by exploring large publicly available metatranscriptomics datasets, we report the identification and genomic characterization of 93 novel viruses with genetic and evolutionary cues of cytorhabdoviruses. Strikingly, five unprecedented viruses with tri-segmented genomes were also identified. This finding represents the first tri-segmented viruses in the family Rhabdoviridae, and they should be classified in a novel genus within this family for which we suggest the name "Trirhavirus". Interestingly, the nucleocapsid and polymerase were the only typical rhabdoviral proteins encoded by those tri-segmented viruses, whereas in three of them, a protein similar to the emaravirus (family Fimoviridae) silencing suppressor was found, while the other predicted proteins had no matches in any sequence databases. Genetic distance and evolutionary insights suggest that all these novel viruses may represent members of novel species. Phylogenetic analyses, of both novel and previously classified plant rhabdoviruses, provide compelling support for the division of the genus Cytorhabdovirus into three distinct genera. This proposed reclassification not only enhances our understanding of the evolutionary dynamics within this group of plant rhabdoviruses but also illuminates the remarkable genomic diversity they encompass. This study not only represents a significant expansion of the genomics of cytorhabdoviruses that will enable future research on the evolutionary peculiarity of this genus but also shows the plasticity in the rhabdovirus genome organization with the discovery of tri-segmented members with a unique evolutionary trajectory.

The Draft Genome of the "Golden Tide" Seaweed, Sargassum horneri: Characterization and Comparative Analysis

Sargassum horneri, a prevalent species of brown algae found along the coast of the northwest Pacific Ocean, holds significant importance as a valuable source of bioactive compounds. However, its rapid growth can lead to the formation of a destructive "golden tide", causing severe damage to the local economy and coastal ecosystems. In this study, we carried out de novo whole-genome sequencing of S. horneri using next-generation sequencing to unravel the genetic information of this alga. By utilizing a reference-guided de novo assembly pipeline with a closely related species, we successfully established a final assembled genome with a total length of 385 Mb. Repetitive sequences made up approximately 30.6% of this genome. Among the identified putative genes, around 87.03% showed homology with entries in the NCBI non-redundant protein database, with Ectocarpus siliculosus being the most closely related species for approximately one-third of these genes. One gene encoding an alkaline phosphatase family protein was found to exhibit positive selection, which could give a clue for the formation of S. horneri golden tides. Additionally, we characterized putative genes involved in fucoidan biosynthesis metabolism, a significant pathway in S. horneri. This study represents the first genome-wide characterization of a S. horneri species, providing crucial insights for future investigations, such as ecological genomic analyses.

A low-coverage 3' RNA-seq to detect homeolog expression in polyploid wheat

Although allopolyploid species are common among natural and crop species, it is not easy to distinguish duplicated genes, known as homeologs, during their genomic analysis. Yet, cost-efficient RNA sequencing (RNA-seq) is to be developed for large-scale transcriptomic studies such as time-series analysis and genome-wide association studies in allopolyploids. In this study, we employed a 3' RNA-seq utilizing 3' untranslated regions (UTRs) containing frequent mutations among homeologous genes, compared to coding sequence. Among the 3' RNA-seq protocols, we examined a low-cost method Lasy-Seq using an allohexaploid bread wheat, Triticum aestivum. HISAT2 showed the best performance for 3' RNA-seq with the least mapping errors and quick computational time. The number of detected homeologs was further improved by extending 1 kb of the 3' UTR annotation. Differentially expressed genes in response to mild cold treatment detected by the 3' RNA-seq were verified with high-coverage conventional RNA-seq, although the latter detected more differentially expressed genes. Finally, downsampling showed that even a 2 million sequencing depth can still detect more than half of expressed homeologs identifiable by the conventional 32 million reads. These data demonstrate that this low-cost 3' RNA-seq facilitates large-scale transcriptomic studies of allohexaploid wheat and indicate the potential application to other allopolyploid species.

Ecophysiological Differentiation among Two Resurrection Ferns and Their Allopolyploid Derivative

Theoretically, the coexistence of diploids and related polyploids is constrained by reproductive and competitive mechanisms. Although niche differentiation can explain the commonly observed co-occurrence of cytotypes, the underlying ecophysiological differentiation among cytotypes has hardly been studied. We compared the leaf functional traits of the allotetraploid resurrection fern Oeosporangium tinaei (HHPP) and its diploid parents, O. hispanicum (HH) and O. pteridioides (PP), coexisting in the same location. Our experimental results showed that all three species can recover physiological status after severe leaf dehydration, which confirms their 'resurrection' ability. However, compared with PP, HH had much higher investment per unit area of light-capturing surface, lower carbon assimilation rate per unit mass for the same midday water potential, higher non-enzymatic antioxidant capacity, higher carbon content, and lower contents of nitrogen, phosphorus, and other macronutrients. These traits allow HH to live in microhabitats with less availability of water and nutrients (rock crevices) and to have a greater capacity for resurrection. The higher assimilation capacity and lower antioxidant capacity of PP explain its more humid and nutrient-rich microhabitats (shallow soils). HHPP traits were mostly intermediate between those of HH and PP, and they allow the allotetraploid to occupy the free niche space left by the diploids.

Neutral processes underlying the macro eco-evolutionary dynamics of mixed-ploidy systems

Polyploidy, i.e. the occurrence of multiple sets of chromosomes, is regarded as an important phenomenon in plant ecology and evolution, with all flowering plants likely having a polyploid ancestry. Owing to genome shock, minority cytotype exclusion and reduced fertility, polyploids emerging in diploid populations are expected to face significant challenges to successful establishment. Their establishment and persistence are often explained by possible fitness or niche differences that would relieve the competitive pressure with diploid progenitors. Experimental evidence for such advantages is, however, not unambiguous, and considerable niche overlap exists among most polyploid species and their diploid counterparts. Here, we develop a neutral spatially explicit eco-evolutionary model to understand whether neutral processes can explain the eco-evolutionary patterns of polyploids. We present a general mechanism for polyploid establishment by showing that sexually reproducing organisms assemble in space in an iterative manner, reducing frequency-dependent mating disadvantages and overcoming potential reduced fertility issues. Moreover, we construct a mechanistic theoretical framework that allows us to understand the long-term evolution of mixed-ploidy populations and show that our model is remarkably consistent with recent phylogenomic estimates of species extinctions in the Brassicaceae family.

Phenotypic variation of a new synthetic allotetraploid Arabidopsis kamchatica enhanced in natural environment

The phenotypic variation of vegetative organs and reproductive organs of newly synthesized and natural Arabidopsis kamchatica genotypes was investigated in both a controlled environment and a natural environment in an experimental garden. When we compared the variation of their leaf shape as a vegetative organ, the synthetic A. kamchatica individuals grown in the garden showed larger variation compared with the individuals incubated in a growth chamber, suggesting enhanced phenotypic variation in a natural fluctuating environment. In contrast, the natural A. kamchatica genotypes did not show significant change in variation by growth condition. The phenotypic variation of floral organs by growth condition was much smaller in both synthetic and natural A. kamchatica genotypes, and the difference in variation width between the growth chamber and the garden was not significant in each genotype as well as among genotypes. The higher phenotypic variation in synthetic leaf may imply flexible transcriptomic regulation of a newly synthesized polyploid compared with a natural polyploid.

Allopolyploid origin and niche expansion of Rhodiola integrifolia (Crassulaceae)

Polyploidy after hybridization between species can lead to immediate post-zygotic isolation, causing saltatory origin of new species. Although the incidence of polyploidization in plants is high, it is thought that a new polyploid lineage can succeed only if it establishes a new ecological niche divergent from its progenitor lineages. We tested the hypothesis that Rhodiola integrifolia from North America is an allopolyploid produced by R. rhodantha and R. rosea and determined whether its survival can be explained by the niche divergence hypothesis. To this end, we sequenced two low-copy nuclear genes (ncpGS and rpb2) in a phylogenetic analysis of 42 Rhodiola species and tested for niche equivalency and similarity using Schoener's D as the index of niche overlap. Our phylogeny-based approach showed that R. integrifolia possesses alleles from both R. rhodantha and R. rosea. Dating analysis showed that the hybridization event that led to R. integrifolia occurred ca. 1.67 Mya and niche modeling analysis showed that at this time, both R. rosea and R. rhodantha may have been present in Beringia, providing the opportunity for the hybridization event. We also found that the niche of R. integrifolia differs from that of its progenitors in both niche breadth and optimum. Taken together, these results confirm the hybrid origin of R. integrifolia and support the niche divergence hypothesis for this tetraploid species. Our results underscore the fact that lineages with no current overlapping distribution could produce hybrid descendants in the past, when climate oscillations made their distributions overlap.

Niche divergence at the intraspecific level in an endemic rare peony (Paeonia rockii): A phylogenetic, climatic and environmental survey

Ecological factors have received increasing attention as drivers of speciation but also in the maintenance of postspeciation divergence. However, the relative significance of the responses of species to climate oscillations for driving niche divergence or conservatism in the evolution of many species that pass through diverse environments and limited geographical boundaries remains poorly understood. Paeonia rockii (one of the ancient species of Paeonia) comprising two subspecies called Paeonia rockii subsp. rockii and Paeonia rockii subsp. taibaishanica is an endemic, rare, and endangered medicinal plant in China. In this study, we integrated whole chloroplast genomes, and ecological factors to obtain insights into ecological speciation and species divergence in this endemic rare peony. RAxML analysis indicated that the topological trees recovered from three different data sets were identical, where P. rockii subsp. rockii and P. rockii subsp. taibaishanica clustered together, and molecular dating analyses suggested that the two subspecies diverged 0.83 million years ago. In addition, ecological niche modeling showed that the predicted suitable distribution areas for P. rockii subsp. rockii and P. rockii subsp. taibaishanica differed considerably, although the predicted core distribution areas were similar, where the population contracted in the last interglacial and expanded in the last glacial maximum. Under the emissions scenarios for the 2050s and 2070s, the suitable distribution areas were predicted to contract significantly, where the migration routes of the two subspecies tended to migrate toward high latitudes and elevations, thereby suggesting strong responses of the distributions of the two subspecies to climate change. These findings combined with the phylogeographic relationships provide comprehensive insights into niche variation and differentiation in this endemic rare peony, and they highlight the importance of geological and climatic changes for species divergence and changes in the population geographic patterns of rare and endangered medicinal plants in East Asia.

Robustness and the generalist niche of polyploid species: Genome shock or gradual evolution?

The prevalence of polyploidy in wild and crop species has stimulated debate over its evolutionary advantages and disadvantages. Previous studies have focused on changes occurring at the polyploidization events, including genome-wide changes termed "genome shock," as well as ancient polyploidy. Recent bioinformatics advances and empirical studies of Arabidopsis and wheat relatives are filling a research gap: the functional evolutionary study of polyploid species using RNA-seq, DNA polymorphism, and epigenomics. Polyploid species can become generalists in natura through environmental robustness by inheriting and merging parental stress responses. Their evolvability is enhanced by mutational robustness working on inherited standing variation. The identification of key genes responsible for gradual adaptive evolution will encourage synthetic biological approaches to transfer polyploid advantages to other species.

Novel Approaches for Species Concepts and Delimitation in Polyploids and Hybrids

Hybridization and polyploidization are important processes for plant evolution. However, classification of hybrid or polyploid species has been notoriously difficult because of the complexity of processes and different evolutionary scenarios that do not fit with classical species concepts. Polyploid complexes are formed via combinations of allopolyploidy, autopolyploidy and homoploid hybridization with persisting sexual reproduction, resulting in many discrete lineages that have been classified as species. Polyploid complexes with facultative apomixis result in complicated net-work like clusters, or rarely in agamospecies. Various case studies illustrate the problems that apply to traditional species concepts to hybrids and polyploids. Conceptual progress can be made if lineage formation is accepted as an inevitable consequence of meiotic sex, which is established already in the first eukaryotes as a DNA restoration tool. The turnaround of the viewpoint that sex forms species as lineages helps to overcome traditional thinking of species as "units". Lineage formation and self-sustainability is the prerequisite for speciation and can also be applied to hybrids and polyploids. Species delimitation is aided by the improved recognition of lineages via various novel -omics methods, by understanding meiosis functions, and by recognizing functional phenotypes by considering morphological-physiological-ecological adaptations.