The evolutionary significance of polyploidy

Core cellular and tissue-specific mechanisms enable desiccation tolerance in Craterostigma

Resurrection plants can survive prolonged life without water (anhydrobiosis) in regions with seasonal drying. This desiccation tolerance requires the coordination of numerous cellular processes across space and time, and individual plant tissues face unique constraints related to their function. Here, we analyzed the complex, octoploid genome of the model resurrection plant Craterostigma (C. plantagineum), and surveyed spatial and temporal expression dynamics to identify genetic elements underlying desiccation tolerance. Homeologous genes within the Craterostigma genome have divergent expression profiles, suggesting the subgenomes contribute differently to desiccation tolerance traits. The Craterostigma genome contains almost 200 tandemly duplicated early light-induced proteins, a hallmark trait of desiccation tolerance, with massive upregulation under water deficit. We identified a core network of desiccation-responsive genes across all tissues, but observed almost entirely unique expression dynamics in each tissue during recovery. Roots and leaves have differential responses related to light and photoprotection, autophagy and nutrient transport, reflecting their divergent functions. Our findings highlight a universal set of likely ancestral desiccation tolerance mechanisms to protect cellular macromolecules under anhydrobiosis, with secondary adaptations related to tissue function.

Fewer chromosomes, more co-occurring species within plant lineages: A likely effect of local survival and colonization

PREMISE

Plant lineages differ markedly in species richness globally, regionally, and locally. Differences in whole-genome characteristics (WGCs) such as monoploid chromosome number, genome size, and ploidy level may explain differences in global species richness through speciation or global extinction. However, it is unknown whether WGCs drive species richness within lineages also in a recent, postglacial regional flora or in local plant communities through local extinction or colonization and regional species turnover.

METHODS

We tested for relationships between WGCs and richness of angiosperm families across the Netherlands/Germany/Czechia as a region, and within 193,449 local vegetation plots.

RESULTS

Families that are species-rich across the region have lower ploidy levels and small monoploid chromosomes numbers or both (interaction terms), but the relationships disappear after accounting for continental and local richness of families. Families that are species-rich within occupied localities have small numbers of polyploidy and monoploid chromosome numbers or both, independent of their own regional richness and the local richness of all other locally co-occurring species in the plots. Relationships between WGCs and family species-richness persisted after accounting for niche characteristics and life histories.

CONCLUSIONS

Families that have few chromosomes, either monoploid or holoploid, succeed in maintaining many species in local communities and across a continent and, as indirect consequence of both, across a region. We suggest evolutionary mechanisms to explain how small chromosome numbers and ploidy levels might decrease rates of local extinction and increase rates of colonization. The genome of a macroevolutionary lineage may ultimately control whether its species can ecologically coexist.

Evolution of sex in crops: recurrent scrap and rebuild

Sexuality is the main strategy for maintaining genetic diversity within a species. In flowering plants (angiosperms), sexuality is derived from ancestral hermaphroditism and multiple sexualities can be expressed in an individual. The mechanisms conferring chromosomal sex determination in plants (or dioecy) have been studied for over a century by both biologists and agricultural scientists, given the importance of this field for crop cultivation and breeding. Despite extensive research, the sex determining gene(s) in plants had not been identified until recently. In this review, we dissect plant sex evolution and determining systems, with a focus on crop species. We introduced classic studies with theoretical, genetic, and cytogenic approaches, as well as more recent research using advanced molecular and genomic techniques. Plants have undergone very frequent transitions into, and out of, dioecy. Although only a few sex determinants have been identified in plants, an integrative viewpoint on their evolutionary trends suggests that recurrent neofunctionalization events are potentially common, in a "scrap and (re)build" cycle. We also discuss the potential association between crop domestication and transitions in sexual systems. We focus on the contribution of duplication events, which are particularly frequent in plant taxa, as a trigger for the creation of new sexual systems.

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.

The alternative transcription and expression characterization of Dmc1 in autotetraploid Carassius auratus

Established autotetraploids often have a highly stable meiosis with high fertility compared with neo-autotetraploids. The autotetraploid Carassius auratus (4n = 200, RRRR) (4nRR), which stemmed from whole-genome duplication of Carassius auratus red var. (2n = 100, RR) (RCC), produces diploid gametes with an adopted diploid-like chromosome pairing in meiosis and maintains the formation of autotetraploid lineages. In this study, we focused on Dmc1, a meiosis-specific recombinase during the prophase of meiosis I, and elaborated on the genetic variation, alternative transcription, expression characterization, and epigenetic modification of Dmc1 in RCC and 4nRR. Two original Dmc1 from RCC were identified in 4nRR, and two duplicated Dmc1 differences in genetic composition were observed in 4nRR. Furthermore, we only noticed that one original and one duplicated Dmc1 were expressed in RCC and 4nRR, respectively. However, both possessed identical gene expression profiles, differential expression of sexual dimorphism, and hypomethylation levels. These results indicated that the specific expression of duplicated Dmc1 may be involve in the progression of meiosis of the diploid-like chromosome pairing in autotetraploid Carassius auratus. Herein, the findings significantly increase knowledge of meiosis of autopolyploid fish and provide meaningful insights into genetic breeding in polyploidy fish.

Patterns of Chromosomal Variation, Homoeologous Exchange, and Their Relationship with Genomic Features in Early Generations of a Synthetic Rice Segmental Allotetraploid

Polyploidization is a driving force in plant evolution. Chromosomal variation often occurs at early generations following polyploid formation due to meiotic pairing irregularity that may compromise segregation fidelity and cause homoeologous exchange (HE). The trends of chromosomal variation and especially factors affecting HE remain to be fully deciphered. Here, by whole-genome resequencing, we performed nuanced analyses of patterns of chromosomal number variation and explored genomic features that affect HE in two early generations of a synthetic rice segmental allotetraploid. We found a wide occurrence of whole-chromosome aneuploidy and, to a lesser extent, also large segment gains/losses in both generations (S2 and S4) of the tetraploids. However, while the number of chromosome gains was similar between S2 and S4, that of losses in S4 was lower than in S2. HEs were abundant across all chromosomes in both generations and showed variable correlations with different genomic features at chromosomal and/or local scales. Contents of genes and transposable elements (TEs) were positively and negatively correlated with HE frequencies, respectively. By dissecting TEs into different classes, retrotransposons were found to be negatively correlated with HE frequency to a stronger extent than DNA transposons, whereas miniature terminal inverted elements (MITEs) showed a strong positive correlation. Local HE frequencies in the tetraploids and homologous recombination (HR) rates in diploids within 1 Mb sliding windows were significantly correlated with each other and showed similar overall distribution profiles. Nonetheless, non-concordant trends between HE and HR rates were found at distal regions in some chromosomes. At local scale, both shared and polymorphic retrotransposons between parents were negatively correlated with HE frequency; in contrast, both shared and polymorphic MITEs showed positive correlations with HE frequency. Our results shed new light on the patterns of chromosomal number variation and reveal genomic features influencing HE frequency in early generations following plant polyploidization.

SCL14 Inhibits the Functions of the NAC043-MYB61 Signaling Cascade to Reduce the Lignin Content in Autotetraploid Populus hopeiensis

Whole-genome duplication often results in a reduction in the lignin content in autopolyploid plants compared with their diploid counterparts. However, the regulatory mechanism underlying variation in the lignin content in autopolyploid plants remains unclear. Here, we characterize the molecular regulatory mechanism underlying variation in the lignin content after the doubling of homologous chromosomes in Populus hopeiensis. The results showed that the lignin content of autotetraploid stems was significantly lower than that of its isogenic diploid progenitor throughout development. Thirty-six differentially expressed genes involved in lignin biosynthesis were identified and characterized by RNA sequencing analysis. The expression of lignin monomer synthase genes, such as PAL, COMT, HCT, and POD, was significantly down-regulated in tetraploids compared with diploids. Moreover, 32 transcription factors, including MYB61, NAC043, and SCL14, were found to be involved in the regulatory network of lignin biosynthesis through weighted gene co-expression network analysis. We inferred that SCL14, a key repressor encoding the DELLA protein GAI in the gibberellin (GA) signaling pathway, might inhibit the NAC043-MYB61 signaling functions cascade in lignin biosynthesis, which results in a reduction in the lignin content. Our findings reveal a conserved mechanism in which GA regulates lignin synthesis after whole-genome duplication; these results have implications for manipulating lignin production.

Genome-Wide Identification and Expression Analysis of Cysteine-Rich Polycomb-like Protein (CPP) Gene Family in Tomato

The cysteine-rich polycomb-like protein (CPP) gene family is a class of transcription factors containing conserved cysteine-rich CRC structural domains that is involved in the regulation of plant growth and stress tolerance to adversity. Relative to other gene families, the CPP gene family has not received sufficient attention. In this study, six SlCPPs were identified for the first time using the most recent genome-wide identification data of tomato. Subsequently, a phylogenetic analysis classified SlCPPs into four subfamilies. The analysis of cis-acting elements in the promoter indicates that SlCPPs are involved in plant growth and development and also stress response. We present for the first time the prediction of the tertiary structure of these SlCPPs proteins using the AlphaFold2 artificial intelligence system developed by the DeepMind team. Transcriptome data analysis showed that SlCPPs were differentially expressed in different tissues. Gene expression profiling showed that all SlCPPs except SlCPP5 were up-regulated under drought stress; SlCPP2, SlCPP3 and SlCPP4 were up-regulated under cold stress; SlCPP2 and SlCPP5 were up-regulated under salt stress; all SlCPPs were up-regulated under inoculation with Cladosporium fulvum; and SlCPP1, SlCPP3, and SlCPP4 were up-regulated under inoculation with Stemphylium lycopersici. We performed a virus-induced gene silencing experiment on SlCPP3, and the results indicated that SlCPP3 was involved in the response to drought stress. Finally, we predicted the interaction network of the key gene SlCPP3, and there was an interaction relationship between SlCPP3 and 10 genes, such as RBR1 and MSI1. The positive outcome showed that SlCPPs responded to environmental stress. This study provides a theoretical and empirical basis for the response mechanisms of tomato in abiotic stresses.

Evolutionary conservation of sequence motifs at sites of protein modification

Gene duplications are common in biology and are likely to be an important source of functional diversification and specialization. The yeast Saccharomyces cerevisiae underwent a whole genome duplication event early in evolution, and a substantial number of duplicated genes have been retained. We identified more than 3,500 instances where only one of two paralogous proteins undergoes post-translational modification despite having retained the same amino acid residue in both. We also developed a web-based search algorithm (CoSMoS.c.) that scores conservation of amino acid sequences based on 1011 wild and domesticated yeast isolates and used it to compare differentially-modified pairs of paralogous proteins. We found that the most common modifications - phosphorylation, ubiquitylation and acylation but not N-glycosylation - occur in regions of high sequence conservation. Such conservation is evident even for ubiquitylation and succinylation, where there is no established 'consensus site' for modification. Differences in phosphorylation were not associated with predicted secondary structure or solvent accessibility, but did mirror known differences in kinase-substrate interactions. Thus, differences in post-translational modification likely result from differences in adjoining amino acids and their interactions with modifying enzymes. By integrating data from large scale proteomics and genomics analysis, in a system with such substantial genetic diversity, we obtained a more comprehensive understanding of the functional basis for genetic redundancies that have persisted for 100 million years.

Tracing 100 million years of grass genome evolutionary plasticity

Grasses derive from a family of monocotyledonous plants that includes crops of major economic importance such as wheat, rice, sorghum and barley, sharing a common ancestor some 100 million years ago. The genomic attributes of plant adaptation remain obscure and the consequences of recurrent whole genome duplications (WGD) or polyploidization events, a major force in plant evolution, remain largely speculative. We conducted a comparative analysis of omics data from ten grass species to unveil structural (inversions, fusions, fissions, duplications, substitutions) and regulatory (expression and methylation) basis of genome plasticity, as possible attributes of plant long lasting evolution and adaptation. The present study demonstrates that diverged polyploid lineages sharing a common WGD event often present the same patterns of structural changes and evolutionary dynamics, but these patterns are difficult to generalize across independent WGD events as a result of non-WGD factors such as selection and domestication of crops. Polyploidy is unequivocally linked to the evolutionary success of grasses during the past 100 million years, although it remains difficult to attribute this success to particular genomic consequences of polyploidization, suggesting that polyploids harness the potential of genome duplication, at least partially, in lineage-specific ways. Overall, the present study clearly demonstrates that post-polyploidization reprogramming is more complex than traditionally reported in investigating single species and calls for a critical and comprehensive comparison across independently polyploidized lineages.

Evolutionary and Genomic Diversity of True Polyploidy in Tetrapods

True polyploid organisms have more than two chromosome sets in their somatic and germline cells. Polyploidy is a major evolutionary force and has played a significant role in the early genomic evolution of plants, different invertebrate taxa, chordates, and teleosts. However, the contribution of polyploidy to the generation of new genomic, ecological, and species diversity in tetrapods has traditionally been underestimated. Indeed, polyploidy represents an important pathway of genomic evolution, occurring in most higher-taxa tetrapods and displaying a variety of different forms, genomic configurations, and biological implications. Herein, we report and discuss the available information on the different origins and evolutionary and ecological significance of true polyploidy in tetrapods. Among the main tetrapod lineages, modern amphibians have an unparalleled diversity of polyploids and, until recently, they were considered to be the only vertebrates with closely related diploid and polyploid bisexual species or populations. In reptiles, polyploidy was thought to be restricted to squamates and associated with parthenogenesis. In birds and mammals, true polyploidy has generally been considered absent (non-tolerated). These views are being changed due to an accumulation of new data, and the impact as well as the different evolutionary and ecological implications of polyploidy in tetrapods, deserve a broader evaluation.

The Torreya grandis genome illuminates the origin and evolution of gymnosperm-specific sciadonic acid biosynthesis

Torreya plants produce dry fruits with assorted functions. Here, we report the 19-Gb chromosome-level genome assembly of T. grandis. The genome is shaped by ancient whole-genome duplications and recurrent LTR retrotransposon bursts. Comparative genomic analyses reveal key genes involved in reproductive organ development, cell wall biosynthesis and seed storage. Two genes encoding a C₁₈ Δ⁹-elongase and a C₂₀ Δ⁵-desaturase are identified to be responsible for sciadonic acid biosynthesis and both are present in diverse plant lineages except angiosperms. We demonstrate that the histidine-rich boxes of the Δ⁵-desaturase are crucial for its catalytic activity. Methylome analysis reveals that methylation valleys of the T. grandis seed genome harbor genes associated with important seed activities, including cell wall and lipid biosynthesis. Moreover, seed development is accompanied by DNA methylation changes that possibly fuel energy production. This study provides important genomic resources and elucidates the evolutionary mechanism of sciadonic acid biosynthesis in land plants.

Insights into the Taxonomically Challenging Hexaploid Alpine Shrub Willows of Salix Sections Phylicifoliae and Nigricantes (Salicaceae)

The complex genomic composition of allopolyploid plants leads to morphologically diverse species. The traditional taxonomical treatment of the medium-sized, hexaploid shrub willows distributed in the Alps is difficult based on their variable morphological characters. In this study, RAD sequencing data, infrared-spectroscopy, and morphometric data are used to analyze the phylogenetic relationships of the hexaploid species of the sections Nigricantes and Phylicifoliae in a phylogenetic framework of 45 Eurasian Salix species. Both sections comprise local endemics as well as widespread species. Based on the molecular data, the described morphological species appeared as monophyletic lineages (except for S. phylicifolia s.str. and S. bicolor, which are intermingled). Both sections Phylicifoliae and Nigricantes are polyphyletic. Infrared-spectroscopy mostly confirmed the differentiation of hexaploid alpine species. The morphometric data confirmed the molecular results and supported the inclusion of S. bicolor into S. phylicifolia s.l., whereas the alpine endemic S. hegetschweileri is distinct and closely related to species of the section Nigricantes. The genomic structure and co-ancestry analyses of the hexaploid species revealed a geographical pattern for widespread S. myrsinifolia, separating the Scandinavian from the alpine populations. The newly described S. kaptarae is tetraploid and is grouped within S. cinerea. Our data reveal that both sections Phylicifoliae and Nigricantes need to be redefined.

A clue to the evolutionary history of modern East Asian flora: insights from phylogeography and diterpenoid alkaloid distribution pattern of the Spiraea japonica complex

Each subkingdom of East Asian flora (EAF) has a unique evolutionary history, but which has rarely been described based on phylogeographic studies of EAF species. The Spiraea japonica L. complex, which is widespread in East Asia (EA), has received considerable attention because of the presence of diterpenoid alkaloids (DAs). It provides a proxy for understanding the genetic diversity and DA distribution patterns of species under various environmental conditions associated with the geological background in EA. In the present study, the plastome and chloroplast/nuclear DNA of 71 populations belonging to the S. japonica complex and its congeners were sequenced, combined with DA identification, environmental analyses, and ecological niche modelling, to investigate their phylogenetic relationships, genetic and DAs distribution patterns, biogeography, and demographic dynamics. An "ampliative" S. japonica complex was put forward, comprising all species of Sect. Calospira Ser. Japonicae, of which three evolutionary units carrying their respective unique types of DAs were identified and associated with the regionalization of EAF (referring to the Hengduan Mountains, central China, and east China). Moreover, a transition belt in central China with its biogeographic significance was revealed by genetic and DA distribution patterns from the perspective of ecological adaptation. The origin and onset differentiation of the "ampliative" S. japonica complex was estimated in the early Miocene (22.01/19.44 Ma). The formation of Japanese populations (6.75 Ma) was facilitated by the land bridge, which subsequently had a fairly stable demographic history. The populations in east China have undergone a founder effect after the Last Glacial Maximum, which may have been promoted by the expansion potential of polyploidization. Overall, the in-situ origin and diversification of the "ampliative" S. japonica complex since the early Miocene is a vertical section of the formation and development of modern EAF and was shaped by the geological history of each subkingdom.

Ancestral self-compatibility facilitates the establishment of allopolyploids in Brassicaceae

Self-incompatibility systems based on self-recognition evolved in hermaphroditic plants to maintain genetic variation of offspring and mitigate inbreeding depression. Despite these benefits in diploid plants, for polyploids who often face a scarcity of mating partners, self-incompatibility can thwart reproduction. In contrast, self-compatibility provides an immediate advantage: a route to reproductive viability. Thus, diploid selfing lineages may facilitate the formation of new allopolyploid species. Here, we describe the mechanism of establishment of at least four allopolyploid species in Brassicaceae (Arabidopsis suecica, Arabidopsis kamchatica, Capsella bursa-pastoris, and Brassica napus), in a manner dependent on the prior loss of the self-incompatibility mechanism in one of the ancestors. In each case, the degraded S-locus from one parental lineage was dominant over the functional S-locus of the outcrossing parental lineage. Such dominant loss-of-function mutations promote an immediate transition to selfing in allopolyploids and may facilitate their establishment.

Development of homozygous tetraploid potato and whole genome doubling-induced the enrichment of H3K27ac and potentially enhanced resistance to cold-induced sweetening in tubers

Polyploid plants typically display advantages on some agronomically important traits over their diploid counterparts. Extensive studies have shown genetic, transcriptomic, and epigenetic dynamics upon polyploidization in multiple plant species. However, few studies have unveiled those alternations imposed only by ploidy level, without any interference from heterozygosity. Cultivated potato is highly heterozygous. Thus, in this study, we developed two homozygous autotetraploid lines and one homozygous diploid line in parallel from a homozygous diploid potato. We confirmed their ploidy levels using chloroplast counting and karyotyping. Oligo-FISH and genome re-sequencing validated that these potato lines are nearly homozygous. We investigated variations in phenotypes, transcription, and histone modifications between two ploidies. Both autotetraploid lines produced larger but fewer tubers than the diploid line. Interestingly, each autotetraploid line displayed ploidy-related differential expression for various genes. We also discovered a genome-wide enrichment of H3K27ac in genic regions upon whole-genome doubling (WGD). However, such enrichment was not associated with the differential gene expression between two ploidies. The tetraploid lines may exhibit better resistance to cold-induced sweetening (CIS) than the diploid line in tubers, potentially regulated through the expression of CIS-related key genes, which seems to be associated with the levels of H3K4me3 in cold-stored tubers. These findings will help to understand the impacts of autotetraploidization on dynamics of phenotypes, transcription, and histone modifications, as well as on CIS-related genes in response to cold storage.

Molecular mechanisms of adaptive evolution in wild animals and plants

Wild animals and plants have developed a variety of adaptive traits driven by adaptive evolution, an important strategy for species survival and persistence. Uncovering the molecular mechanisms of adaptive evolution is the key to understanding species diversification, phenotypic convergence, and inter-species interaction. As the genome sequences of more and more non-model organisms are becoming available, the focus of studies on molecular mechanisms of adaptive evolution has shifted from the candidate gene method to genetic mapping based on genome-wide scanning. In this study, we reviewed the latest research advances in wild animals and plants, focusing on adaptive traits, convergent evolution, and coevolution. Firstly, we focused on the adaptive evolution of morphological, behavioral, and physiological traits. Secondly, we reviewed the phenotypic convergences of life history traits and responding to environmental pressures, and the underlying molecular convergence mechanisms. Thirdly, we summarized the advances of coevolution, including the four main types: mutualism, parasitism, predation and competition. Overall, these latest advances greatly increase our understanding of the underlying molecular mechanisms for diverse adaptive traits and species interaction, demonstrating that the development of evolutionary biology has been greatly accelerated by multi-omics technologies. Finally, we highlighted the emerging trends and future prospects around the above three aspects of adaptive evolution.

Artificial Diploid Escherichia coli by a CRISPR Chromosome-Doubling Technique

Synthetic biology has been represented by the creation of artificial life forms at the genomic scale. In this work, a CRISPR-based chromosome-doubling technique is designed to first construct an artificial diploid Escherichia coli cell. The stable single-cell diploid E. coli is isolated by both maximal dilution plating and flow cytometry, and confirmed with quantitative PCR, fluorescent in situ hybridization, and third-generation genome sequencing. The diploid E. coli has a greatly reduced growth rate and elongated cells at 4-5 µm. It is robust against radiation, and the survival rate after exposure to UV increased 40-fold relative to WT. As a novel life form, the artificial diploid E. coli is an ideal substrate for research fundamental questions in life science concerning polyploidy. And this technique may be applied to other bacteria.

Learning to tango with four (or more): the molecular basis of adaptation to polyploid meiosis

Polyploidy, which arises from genome duplication, has occurred throughout the history of eukaryotes, though it is especially common in plants. The resulting increased size, heterozygosity, and complexity of the genome can be an evolutionary opportunity, facilitating diversification, adaptation and the evolution of functional novelty. On the other hand, when they first arise, polyploids face a number of challenges, one of the biggest being the meiotic pairing, recombination and segregation of the suddenly more than two copies of each chromosome, which can limit their fertility. Both for developing polyploidy as a crop improvement tool (which holds great promise due to the high and lasting multi-stress resilience of polyploids), as well as for our basic understanding of meiosis and plant evolution, we need to know both the specific nature of the challenges polyploids face, as well as how they can be overcome in evolution. In recent years there has been a dramatic uptick in our understanding of the molecular basis of polyploid adaptations to meiotic challenges, and that is the focus of this review.

The evolution and immunomodulatory role of Zc3h12 proteins in zebrafish (Danio rerio)

Zc3h12 family is an important RNA-binding protein family regulating mRNA of inflammatory cytokines in mammals. However, there are few studies on their post-transcriptional level regulation of inflammatory cytokines in fish. Here, we investigated the evolution of zebrafish Zc3h12 family and explored their immunomodulatory role. Phylogenetic and syntenic analysis indicated the number of zc3h12 family members had increased ranging from a single member in invertebrates to a single copy of four members in mammals. As the most evolutionarily diverse group of vertebrates, the number of zc3h12 family members was more complex and diverse in the teleost, each member experienced different fates and followed different rules in multiple rounds of whole-genome duplication events. Thereinto, zebrafish contained three zc3h12 genes, among which zc3h12aa and zc3h12ab were duplicated from the same gene. Zebrafish Zc3h12 family could recognize the 3'-UTR regions of inflammatory cytokines through binding to the specific RNA secondary structure and negatively regulate their expression. Deletion of either Zc3h12 domains or mutation of the key amino acid in RNAase domain attenuated their modulatory effect, suggesting both domain and RNAase activity are important to the immunomodulatory role. These results elucidated the evolution of Zc3h12 family and uncovered Zc3h12-mediated post-transcriptional regulation of cytokines in zebrafish.

Translational research from models to crops: comparative genomics for plant breeding

The concept of translational research, which originated in the medical field in the 1980s, consists in improving the efficient transfer of research results obtained in a species (which can be considered as a model or pivot) to all the species for which these results are of interest for its improvement in Agriculture. In this context, comparative genomics is an important tool for translational research, effectively identifying genes controlling common functions between species. Editing and phenotyping tools must thus allow the functional validation of the gene conserved within the species for which the knowledge has been extrapolated, that is to say transferred, and the identification of the best alleles and associated genotypes for exploitation in current breeding programs.

Genetic Diversity of Oxytropis Species from the Center of the Genus Origin: Insight from Molecular Studies

The genus Oxytropis (Fabaceae) was formed from the ancient species of Astragalus presumably approximately 5.6 Ma ago in Southern Siberia. Our study summarized data on the genetic diversity of 69 populations of 31 Oxytropis species in the center of origin of the genus based on the sequencing of plastid genome markers. Most of the populations (82.6%) are characterized by high gene diversity (0.600–1.000), which indicates a relatively stable state. Phylogenetic relationships between most Oxytropis species remain unresolved. Three genetic complexes and four phyletic lineages have been identified. Some species form weakly differentiated complexes, which is probably caused by their relatively recent divergence and the demography processes, as well as interspecific hybridization and polyploidy characteristic of Oxytropis species.

Genome-wide identification and expression analysis of the cotton patatin-related phospholipase A genes and response to stress tolerance

Patatin-related phospholipase A genes were involved in the response of Gossypium hirsutum to drought and salt tolerance. pPLA (patatin-related phospholipase A) is a key enzyme that catalyzes the initial step of lipid hydrolysis, which is involved in biological processes, such as drought, salt stress, and freezing injury. However, a comprehensive analysis of the pPLA gene family in cotton, especially the role of pPLA in the response to drought and salt tolerance, has not been reported so far. A total of 33 pPLA genes were identified in this study using a genome-wide search approach, and phylogenetic analysis classified these genes into three groups. These genes are unevenly distributed on the 26 chromosomes of cotton, and most of them contain a few introns. The results of the collinear analysis showed that G. hirsutum contained 1-5 copies of each pPLA gene found in G. arboreum and G. raimondii. The subcellular localization analysis of Gh_D08G061200 showed that the protein was localized in the nucleus. In addition, analysis of published upland cotton transcriptome data revealed that six GhPLA genes were expressed in various tissues and organs. Two genes (Gh_A04G142100.1 and Gh_D04G181000.1) were highly expressed in all tissues under normal conditions, showing the expression characteristics of housekeeping genes. Under different drought and salt tolerance stresses, we detected four genes with different expression levels. This study helps to clarify the role of pPLA in the response to drought and salt tolerance.

Autopolyploids of Arabidopsis thaliana are more phenotypically plastic than their diploid progenitors

BACKGROUND AND AIMS

Polyploids are often hypothesized to have increased phenotypic plasticity compared with their diploid progenitors, but recent work suggests that the relationship between whole-genome duplication (WGD) and plasticity is not so straightforward. Impacts of WGD on plasticity are moderated by other evolutionary processes in nature, which has impeded generalizations regarding the effects of WGD alone. We assessed shifts in phenotypic plasticity and mean trait values accompanying WGD, as well as the adaptive consequences of these shifts.

METHODS

To isolate WGD effects, we compared two diploid lineages of Arabidopsis thaliana wiht corresponding autotetraploids grown across different salt and nutrient conditions in a growth chamber.

KEY RESULTS

For the few cases in which diploids and polyploids differed in plasticity, polyploids were more plastic, consistent with hypotheses that WGD increases plasticity. Under stress, increased plasticity was often adaptive (associated with higher total seed mass), but in other cases plasticity was unrelated to fitness. Mean trait values and plasticity were equally likely to be affected by WGD, but the adaptive consequences of these shifts were often context dependent or lineage specific. For example, polyploids had extended life spans, a shift that was adaptive in one polyploid lineage under amenable conditions but was maladaptive in the other lineage under stress.

CONCLUSIONS

Our work shows that increased phenotypic plasticity can result from WGD alone, independent of other evolutionary processes. We find that the effects of WGD can differ depending on the genotype of the progenitor and the environmental context. Though our experiment was limited to two genotypes of a single species, these findings support the idea that WGD can indeed increase plasticity.

Genotypic variation rather than ploidy level determines functional trait expression in a foundation tree species in the presence and absence of environmental stress

BACKGROUND AND AIMS

At the population level, genetic diversity is a key determinant of a tree species' capacity to cope with stress. However, little is known about the relative importance of the different components of genetic diversity for tree stress responses. We compared how two sources of genetic diversity, genotype and cytotype (i.e. differences in ploidy levels), influence growth, phytochemical and physiological traits of Populus tremuloides in the presence and absence of environmental stress.

METHODS

In a series of field studies, we first assessed variation in traits across diploid and triploid aspen genotypes from Utah and Wisconsin under non-stressed conditions. In two follow-up experiments, we exposed diploid and triploid aspen genotypes from Wisconsin to individual and interactive drought stress and defoliation treatments and quantified trait variations under stress.

KEY RESULTS

We found that (1) tree growth and associated traits did not differ significantly between ploidy levels under non-stressed conditions. Instead, variation in tree growth and most other traits was driven by genotypic and population differences. (2) Genotypic differences were critical for explaining variation of most functional traits and their responses to stress. (3) Ploidy level played a subtle role in shaping traits and trait stress responses, as its influence was typically obscured by genotypic differences. (4) As an exception to the third conclusion, we showed that triploid trees expressed 17 % higher foliar defence (tremulacin) levels, 11 % higher photosynthesis levels and 23 % higher rubisco activity under well-watered conditions. Moreover, triploid trees displayed greater drought resilience than diploids as they produced 35 % more new tissue than diploids when recovering from drought stress.

CONCLUSION

Although ploidy level can strongly influence the ecology of tree species, those effects may be relatively small in contrast to the effects of genotypic variation in highly diverse species.

New cup out of old coffee: contribution of parental gene expression legacy to phenotypic novelty in coffee beans of the allopolyploid Coffea arabica L

BACKGROUND AND AIMS

Allopolyploidization is a widespread phenomenon known to generate novel phenotypes by merging evolutionarily distinct parental genomes and regulatory networks in a single nucleus. The objective of this study was to investigate the transcriptional regulation associated with phenotypic novelty in coffee beans of the allotetraploid Coffea arabica.

METHODS

A genome-wide comparative transcriptomic analysis was performed in C. arabica and its two diploid progenitors, C. canephora and C. eugenioides. Gene expression patterns and homeologue expression were studied on seeds at five different maturation stages. The involvement of homeologue expression bias (HEB) in specific traits was addressed both by functional enrichment analyses and by the study of gene expression in the caffeine and chlorogenic acid biosynthesis pathways.

KEY RESULTS

Expression-level dominance in C. arabica seed was observed for most of the genes differentially expressed between the species. Approximately a third of the genes analysed showed HEB. This proportion increased during seed maturation but the biases remained equally distributed between the sub-genomes. The relative expression levels of homeologues remained relatively constant during maturation and were correlated with those estimated in leaves of C. arabica and interspecific hybrids between C. canephora and C. eugenioides. Functional enrichment analyses performed on genes exhibiting HEB enabled the identification of processes potentially associated with physiological traits. The expression profiles of the genes involved in caffeine biosynthesis mirror the differences observed in the caffeine content of mature seeds of C. arabica and its parental species.

CONCLUSIONS

Neither of the two sub-genomes is globally preferentially expressed in C. arabica seeds, and homeologues appear to be co-regulated by shared trans-regulatory mechanisms. The observed HEBs are thought to be a legacy of gene expression differences inherited from diploid progenitor species. Pre-existing functional divergences between parental species appear to play an important role in controlling the phenotype of C. arabica seeds.

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.

Polyploidy: its consequences and enabling role in plant diversification and evolution

BACKGROUND

Most, if not all, green plant (Virdiplantae) species including angiosperms and ferns are polyploids themselves or have ancient polyploid or whole genome duplication signatures in their genomes. Polyploids are not only restricted to our major crop species such as wheat, maize, potato and the brassicas, but also occur frequently in wild species and natural habitats. Polyploidy has thus been viewed as a major driver in evolution, and its influence on genome and chromosome evolution has been at the centre of many investigations. Mechanistic models of the newly structured genomes are being developed that incorporate aspects of sequence evolution or turnover (low-copy genes and regulatory sequences, as well as repetitive DNAs), modification of gene functions, the re-establishment of control of genes with multiple copies, and often meiotic chromosome pairing, recombination and restoration of fertility.

SCOPE

World-wide interest in how green plants have evolved under different conditions - whether in small, isolated populations, or globally - suggests that gaining further insight into the contribution of polyploidy to plant speciation and adaptation to environmental changes is greatly needed. Forward-looking research and modelling, based on cytogenetics, expression studies, and genomics or genome sequencing analyses, discussed in this Special Issue of the Annals of Botany, consider how new polyploids behave and the pathways available for genome evolution. They address fundamental questions about the advantages and disadvantages of polyploidy, the consequences for evolution and speciation, and applied questions regarding the spread of polyploids in the environment and challenges in breeding and exploitation of wild relatives through introgression or resynthesis of polyploids.

CONCLUSION

Chromosome number, genome size, repetitive DNA sequences, genes and regulatory sequences and their expression evolve following polyploidy - generating diversity and possible novel traits and enabling species diversification. There is the potential for ever more polyploids in natural, managed and disturbed environments under changing climates and new stresses.

Inference of Ancient Polyploidy from Genomic Data

Whole-genome sequence data have revealed that numerous eukaryotic organisms derive from distant polyploid ancestors, even when these same organisms are genetically and karyotypically diploid. Such ancient whole-genome duplications (WGDs) have been important for long-term genome evolution and are often speculatively associated with important evolutionary events such as key innovations, adaptive radiations, or survival after mass extinctions. Clearly, reliable methods for unveiling ancient WGDs are key toward furthering understanding of the long-term evolutionary significance of polyploidy. In this chapter, we describe a set of basic established comparative genomics approaches for the inference of ancient WGDs from genomic data based on empirical age distributions and collinearity analyses, explain the principles on which they are based, and illustrate a basic workflow using the software "wgd," geared toward a typical exploratory analysis of a newly obtained genome sequence.

Studying Whole-Genome Duplication Using Experimental Evolution of Chlamydomonas

In this chapter, we present the use of Chlamydomonas reinhardtii in experiments designed to study the evolutionary impacts of whole genome duplication. We shortly introduce the algal species and depict why it is an excellent model for experimental evolution. Subsequently, we discuss the most relevant steps and methods in the design of a ploidy-related Chlamydomonas experiment. These steps include strain selection, ploidy determination, different methods of making diplo- and polyploid Chlamydomonas cells, replication, culturing conditions, preservation, and the ways to quantify phenotypic and genotypic change.

Inference of Ancient Polyploidy Using Transcriptome Data

Polyploidizations, or whole-genome duplications (WGDs), in plants have increased biological complexity, facilitated evolutionary innovation, and likely enabled adaptation under harsh conditions. Besides genomic data, transcriptome data have been widely employed to detect WGDs, due to their efficient accessibility to the gene space of a species. Age distributions based on synonymous substitutions (so-called K_S age distributions) for paralogs assembled from transcriptome data have identified numerous WGDs in plants, paving the way for further studies on the importance of WGDs for the evolution of seed and flowering plants. However, it is still unclear how transcriptome-based age distributions compare to those based on genomic data. In this chapter, we implemented three different de novo transcriptome assembly pipelines with two popular assemblers, i.e., Trinity and SOAPdenovo-Trans. We selected six plant species with published genomes and transcriptomes to evaluate how assembled transcripts from different pipelines perform when using K_S distributions to detect previously documented WGDs in the six species. Further, using genes predicted in each genome as references, we evaluated the effects of missing genes, gene family clustering, and de novo assembled transcripts on the transcriptome-based K_S distributions. Our results show that, although the transcriptome-based K_S distributions differ from the genome-based ones with respect to their shapes and scales, they are still reasonably reliable for unveiling WGDs, except in species where most duplicates originated from a recent WGD. We also discuss how to overcome some possible pitfalls when using transcriptome data to identify WGDs.

Analyses of Genome Regulatory Evolution Following Whole-Genome Duplication Using the Phylogenetic EVE Model

Whole-genome duplications (WGDs) are important in shaping the evolution of complex genomes, including rewiring of genome regulation. To address key questions about how WGDs impact the evolution of genome regulation, we need to understand the relative importance of selection versus drift and temporal evolutionary dynamics. One promising class of statistical models that can help address such questions are phylogenetic Ornstein-Uhlenbeck (OU) models.Here we present a computational pipeline for the comparative phylogenetic analyses of genome regulation using an OU model. We have implemented this model in R and provide a step-by-step protocol for the use of this model, including example scripts and simulated test data. We provide the nonspecialist a brief overview of how this model works and how to perform tests for signatures of selection on genome regulation as well as power simulations to aid in experimental design and interpretation of results. We believe that these resources could help polyploidy research move forward in an era of rapidly increasing functional genomics data across the tree of life.

Efficiently Editing Multiple Duplicated Homeologs and Alleles for Recurrent Polyploids

Research on the evolutionary fate of duplicated genes in recurrent polyploids is scarce due to the difficulties in disentangling the different homeologs and alleles of duplicated genes. This chapter describes the detailed procedures to identify different homeologs and alleles of duplicated genes, to analyze their molecular characteristics, and to reveal their functional divergence by gene editing with CRISPR/Cas9 (clustered regularly interspaced short palindromic repeats/CRISPR-associated system 9). Using the gene editing approach, we efficiently constructed multiple knockout mutant lines with single or simultaneously disrupted different homeologs or alleles in a recurrent polyploid fish, demonstrating its usability for targeting and mutating multiple divergent homeologs and alleles in recurrent duplicated genomes.

Revisiting ancient polyploidy in leptosporangiate ferns

Ferns, and particularly homosporous ferns, have long been assumed to have experienced recurrent whole-genome duplication (WGD) events because of their substantially large genome sizes, surprisingly high chromosome numbers, and high degrees of polyploidy among many extant members. As the number of sequenced fern genomes is limited, recent studies have employed transcriptome data to find evidence for WGDs in ferns. However, they have reached conflicting results concerning the occurrence of ancient polyploidy, for instance, in the lineage of leptosporangiate ferns. Because identifying WGDs in a phylogenetic context is the foremost step in studying the contribution of ancient polyploidy to evolution, we here revisited earlier identified WGDs in leptosporangiate ferns, mainly the core leptosporangiate ferns, by building K_S -age distributions and applying substitution rate corrections and by conducting statistical gene tree-species tree reconciliation analyses. Our integrative analyses not only identified four ancient WGDs in the sampled core leptosporangiate ferns but also identified false positives and false negatives for WGDs that recent studies have reported earlier. In conclusion, we underscore the significance of substitution rate corrections and uncertainties in gene tree-species tree reconciliations in calling WGD events and advance an exemplar workflow to overcome such often-overlooked issues.

Biased mutations and gene losses underlying diploidization of the tetraploid broomcorn millet genome

Broomcorn millet (Panicum miliaceum L.) is one of the earliest domesticated crops, and is a valuable resource to secure food diversity and combat drought stresses under the global warming scenario. However, due to the absence of extant diploid progenitors, the polyploidy genome of broomcorn millet remains poorly understood. Here, we report the chromosome-scale genome assembly of broomcorn millet. We divided the broomcorn millet genome into two subgenomes using the genome sequence of Panicum hallii, a diploid relative of broomcorn millet. Our analyses revealed that the two subgenomes diverged at ~4.8 million years ago (Mya), while the allotetraploidization of broomcorn millet may have occurred about ~0.48 Mya, suggesting that broomcorn millet is a relatively recent allotetraploid. Comparative analyses showed that subgenome B was larger than subgenome A in size, which was caused by the biased accumulation of long terminal repeat retrotransposons in the progenitor of subgenome B before polyploidization. Notably, the accumulation of biased mutations in the transposable element-rich subgenome B led to more gene losses. Although no significant dominance of either subgenome was observed in the expression profiles of broomcorn millet, we found the minimally expressed genes in P. hallii tended to be lost during diploidization of broomcorn millet. These results suggest that broomcorn millet is at the early stage of diploidization and that mutations likely occurred more on genes that were marked with lower expression levels.

Autoploid origin and rapid diploidization of the tetraploid Thinopyrum elongatum revealed by genome differentiation and chromosome pairing in meiosis

Polyploidy is a common mode of evolution in flowering plants. Both the natural tetraploid Thinopyrum elongatum and the diploid one from the same population show a diploid-like pairing in meiosis. However, debate on the chromosome composition and origin of the tetraploid Th. elongatum is ongoing. In the present study, we obtained the induced tetraploid Th. elongatum and found that the induced and natural tetraploids are morphologically close, except for slower development and lower seed setting. Using probes developed from single chromosome microdissection and a Fosmid library, obvious differentiations were discovered between two chromosome sets (E₁ and E₂ ) of the natural tetraploid Th. elongatum but not the induced one. Interestingly, hybrid F₁ derived from the two different wheat-tetraploid Th. elongatum amphiploids 8802 and 8803 produced seeds well. More importantly, analysis of meiosis in F₂ individuals revealed that chromosomes from E₁ and E₂ could pair well on the durum wheat background with the presence of Ph1. No chromosome set differentiation on the FISH level was discovered from the S1 to S4 generations in the induced one. In metaphase of the meiosis first division in the natural tetraploid, more pairings were bivalents and fewer quadrivalents with ratio of 13.94 II + 0.03 IV (n = 31). Chromosome pairing configuration in the induced tetraploid is 13.05 II + 0.47 IV (n = 19), with the quadrivalent ratio being only slightly higher than the ratio in the natural tetraploid. Therefore, the natural tetraploid Th. elongatum is of autoploid origin and the induced tetraploid Th. elongatum evolutionarily underwent rapid diploidization in the low generation.

ploidyfrost: Reference-free estimation of ploidy level from whole genome sequencing data based on de Bruijn graphs

Polyploidy is ubiquitous and its consequences are complex and variable. A change of ploidy level generally influences genetic diversity and results in morphological, physiological and ecological differences between cells or organisms with different ploidy levels. To avoid cumbersome experiments and take advantage of the less biased information provided by the vast amounts of genome sequencing data, computational tools for ploidy estimation are urgently needed. Until now, although a few such tools have been developed, many aspects of this estimation, such as the requirement of a reference genome, the lack of informative results and objective inferences, and the influence of false positives from errors and repeats, need further improvement. We have developed ploidyfrost, a de Bruijn graph-based method, to estimate ploidy levels from whole genome sequencing data sets without a reference genome. ploidyfrost provides a visual representation of allele frequency distribution generated using the ggplot2 package as well as quantitative results using the Gaussian mixture model. In addition, it takes advantage of colouring information encoded in coloured de Bruijn graphs to analyse multiple samples simultaneously and to flexibly filter putative false positives. We evaluated the performance of ploidyfrost by analysing highly heterozygous or repetitive samples of Cyclocarya paliurus and a complex allooctoploid sample of Fragaria × ananassa. Moreover, we demonstrated that the accuracy of analysis results can be improved by constraining a threshold such as Cramér's V coefficient on variant features, which may significantly reduce the side effects of sequencing errors and annoying repeats on the graphical structure constructed.

Phylogenomics and the flowering plant tree of life

The advances accelerated by next-generation sequencing and long-read sequencing technologies continue to provide an impetus for plant phylogenetic study. In the past decade, a large number of phylogenetic studies adopting hundreds to thousands of genes across a wealth of clades have emerged and ushered plant phylogenetics and evolution into a new era. In the meantime, a roadmap for researchers when making decisions across different approaches for their phylogenomic research design is imminent. This review focuses on the utility of genomic data (from organelle genomes, to both reduced representation sequencing and whole-genome sequencing) in phylogenetic and evolutionary investigations, describes the baseline methodology of experimental and analytical procedures, and summarizes recent progress in flowering plant phylogenomics at the ordinal, familial, tribal, and lower levels. We also discuss the challenges, such as the adverse impact on orthology inference and phylogenetic reconstruction raised from systematic errors, and underlying biological factors, such as whole-genome duplication, hybridization/introgression, and incomplete lineage sorting, together suggesting that a bifurcating tree may not be the best model for the tree of life. Finally, we discuss promising avenues for future plant phylogenomic studies.

Recent allopolyploidy alters Spartina microRNA expression in response to xenobiotic-induced stress

Environmental contamination by xenobiotics represents a major threat for natural ecosystems and public health. In response, xenobiotic detoxification is a fundamental trait of organisms for developmental plasticity and stress tolerance, but the underlying molecular mechanisms remain poorly understood in plants. To decipher this process, we explored the consequences of allopolyploidy on xenobiotic tolerance in the genus Spartina Schreb. Specifically, we focused on microRNAs (miRNAs) owing to their central function in the regulation of gene expression patterns, including responses to stress. Small RNA-Seq was conducted on the parents S. alterniflora and S. maritima, their F1 hybrid S. x townsendii and the allopolyploid S. anglica under phenanthrene-induced stress (phe), a model Polycyclic Aromatic Hydrocarbon (PAH) compound. Differentially expressed miRNAs in response to phe were specifically identified within species. In complement, the respective impacts of hybridization and genome doubling were detected, through changes in miRNA expression patterns between S. x townsendii, S. anglica and the parents. The results support the impact of allopolyploidy in miRNA-guided regulation of plant response to phe. In total, we identified 17 phe-responsive miRNAs in Spartina among up-regulated MIR156 and down-regulated MIR159. We also describe novel phe-responsive miRNAs as putative Spartina-specific gene expression regulators in response to stress. Functional validation using Arabidopsis (L.) Heynh. T-DNA lines inserted in homologous MIR genes was performed, and the divergence of phe-responsive miRNA regulatory networks between Arabidopsis and Spartina was discussed.

The immediate effects of polyploidization of Spirodela polyrhiza change in a strain-specific way along environmental gradients

The immediate effects of plant polyploidization are well characterized and it is generally accepted that these morphological, physiological, developmental, and phenological changes contribute to polyploid establishment. Studies on the environmental dependence of the immediate effects of whole-genome duplication (WGD) are, however, scarce but suggest that these immediate effects are altered by stressful conditions. As polyploid establishment seems to be associated with environmental disturbance, the relationship between ploidy-induced phenotypical changes and environmental conditions is highly relevant. Here, we use a common garden experiment on the greater duckweed Spirodela polyrhiza to test whether the immediate effects of WGD can facilitate the establishment of tetraploid duckweed along gradients of two environmental stressors. Because successful polyploid establishment often depends on recurrent polyploidization events, we include four genetically diverse strains and assess whether these immediate effects are strain-specific. We find evidence that WGD can indeed confer a fitness advantage under stressful conditions and that the environment affects ploidy-induced changes in fitness and trait reaction norms in a strain-specific way.

Deciphering evolutionary dynamics of WRKY genes in Arachis species

BACKGROUND

Cultivated peanut (Arachis hypogaea), a progeny of the cross between A. duranensis and A. ipaensis, is an important oil and protein crop from South America. To date, at least six Arachis genomes have been sequenced. WRKY transcription factors (TFs) play crucial roles in plant growth, development, and response to abiotic and biotic stresses. WRKY TFs have been identified in A. duranensis, A. ipaensis, and A. hypogaea cv. Tifrunner; however, variations in their number and evolutionary patterns across various Arachis spp. remain unclear.

RESULTS

WRKY TFs were identified and compared across different Arachis species, including A. duranensis, A. ipaensis, A. monticola, A. hypogaea cultivars (cv.) Fuhuasheng, A. hypogaea cv. Shitouqi, and A. hypogaea cv. Tifrunner. The results showed that the WRKY TFs underwent dynamic equilibrium between diploid and tetraploid peanut species, characterized by the loss of old WRKY TFs and retention of the new ones. Notably, cultivated peanuts inherited more conserved WRKY orthologs from wild tetraploid peanuts than their wild diploid donors. Analysis of the W-box elements and protein-protein interactions revealed that different domestication processes affected WRKY evolution across cultivated peanut varieties. WRKY TFs of A. hypogaea cv. Fuhuasheng and Shitouqi exhibited a similar domestication process, while those of cv. Tifrunner of the same species underwent a different domestication process based on protein-protein interaction analysis.

CONCLUSIONS

This study provides new insights into the evolution of WRKY TFs in Arachis spp.

Transposable elements are associated with genome-specific gene expression in bread wheat

INTRODUCTION

Recent studies in wheat emphasized the importance of TEs, which occupy ~85% of the wheat genome, as a major source of intraspecific genetic variation due to their recent activity and involvement in genomic rearrangements. The contribution of TEs to structural and functional variations in bread wheat genes is not fully understood.

METHODS

Here, publicly available RNA-Seq databases of bread wheat were integrated to identify TE insertions within gene bodies (exons\ introns) and assess the impact of TE insertions on gene expression variations of homoeologs gene groups. Overall, 70,818 homoeologs genes were analyzed: 55,170 genes appeared in each one of the three subgenomes (termed ABD), named triads; 12,640 genes appeared in two of the three subgenomes (in A and B only, termed AB; or in A and D only, termed AD; or in B and D only, termed BD);, named dyads; and 3,008 genes underwent duplication in one of the three subgenomes (two copies in: subgenome A, termed AABD; subgenome B, termed ABBD; or subgenome D, termed ABDD), named tetrads.

RESULTS

To this end, we found that ~36% of the 70,818 genes contained at least one TE insertion within the gene body, mostly in triads. Analysis of 14,258 triads revealed that the presence of TE insertion in at least one of the triad genes (7,439 triads) was associated with balanced expression (similar expression levels) between the homoeolog genes. TE insertions within the exon or in the untranslated regions (UTRs) of one or more of the homoeologs in a triad were significantly associated with homoeolog expression bias. Furthermore, we found a statistically significant correlation between the presence\absence of TEs insertions belonging to six TE superfamilies and 17 TE subfamilies and the suppression of a single homoeolog gene. A significant association was observed between the presence of TE insertions from specific superfamilies and the expression of genes that are associated with biotic and abiotic stress responses.

CONCLUSION

Our data strongly indicate that TEs might play a prominent role in controlling gene expression in a genome-specific manner in bread wheat.

The small subunit of Rubisco and its potential as an engineering target

Rubisco catalyses the first rate-limiting step in CO2 fixation and is responsible for the vast majority of organic carbon present in the biosphere. The function and regulation of Rubisco remain an important research topic and a longstanding engineering target to enhance the efficiency of photosynthesis for agriculture and green biotechnology. The most abundant form of Rubisco (Form I) consists of eight large and eight small subunits, and is found in all plants, algae, cyanobacteria, and most phototrophic and chemolithoautotrophic proteobacteria. Although the active sites of Rubisco are located on the large subunits, expression of the small subunit regulates the size of the Rubisco pool in plants and can influence the overall catalytic efficiency of the Rubisco complex. The small subunit is now receiving increasing attention as a potential engineering target to improve the performance of Rubisco. Here we review our current understanding of the role of the small subunit and our growing capacity to explore its potential to modulate Rubisco catalysis using engineering biology approaches.

Genetic Redundancy in Rye Shows in a Variety of Ways

Fifty years ago Susumu Ohno formulated the famous C-value paradox, which states that there is no correlation between the physical sizes of the genome, i.e., the amount of DNA, and the complexity of the organism, and highlighted the problem of genome redundancy. DNA that does not have a positive effect on the fitness of organisms has been characterized as "junk or selfish DNA". The controversial concept of junk DNA remains viable. Rye is a convenient subject for yet another test of the correctness and scientific significance of this concept. The genome of cultivated rye, Secale cereale L., is considered one of the largest among species of the tribe Triticeae and thus it tops the average angiosperm genome and the genomes of its closest evolutionary neighbors, such as species of barley, Hordeum (by approximately 30-35%), and diploid wheat species, Triticum (approximately 25%). The review provides an analysis of the structural organization of various regions of rye chromosomes with a description of the molecular mechanisms contributing to their size increase during evolution and the classes of DNA sequences involved in these processes. The history of the development of the concept of eukaryotic genome redundancy is traced and the current state of this problem is discussed.

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.

Developing a CRISPR System in Nongenetic Model Polyploids

The genetic consequences following polyploidy (i.e., whole-genome duplication; WGD) vary greatly across organisms and through time since polyploidization. At the gene level in allopolyploids, changes include loss/retention of both parental gene copies, function/expression divergence between the two parental copies, and silencing of one parental copy. Functional studies of genes with different retention patterns contribute to a better understanding of the genetic factors underlying the success of polyploids. Most research on gene functions to date focuses on a few well-established genetic models or crops. However, many species that best exemplify the polyploidy process are nongenetic models; the lack of an efficient genome editing system hinders functional studies in these systems. In this chapter, we discuss the considerations of developing CRISPR, a robust and efficient genome editing system, in polyploid plants that are not genetic models. We use diploid and polyploid Tragopogon (Asteraceae) as examples of a well-studied evolutionary model system for which abundant genetic and genomic resources are lacking. Using this system, we provide our protocols for sgRNA design, plasmid construction, a useful protoplast transient assay, and a plant transformation method we developed for this system. We also provide suggestions for possible modifications to these protocols to help promote successful application to other non-models. With the rapid applications of CRISPR in plant sciences, the broad adaptation of CRISPR in studies of the evolutionary significance of WGD holds enormous potential. We hope our studies and methods developed for polyploid Tragopogon will provide a guideline for establishing a CRISPR system in other nongenetic model polyploids of evolutionary or other interest.

Sequencing and Assembly of Polyploid Genomes

Polyploidy has been observed throughout major eukaryotic clades and has played a vital role in the evolution of angiosperms. Recent polyploidizations often result in highly complex genome structures, posing challenges to genome assembly and phasing. Recent advances in sequencing technologies and genome assembly algorithms have enabled high-quality, near-complete chromosome-level assemblies of polyploid genomes. Advances in novel sequencing technologies include highly accurate single-molecule sequencing with HiFi reads, chromosome conformation capture with Hi-C technique, and linked reads sequencing. Additionally, new computational approaches have also significantly improved the precision and reliability of polyploid genome assembly and phasing, such as HiCanu, hifiasm, ALLHiC, and PolyGembler. Herein, we review recently published polyploid genomes and compare the various sequencing, assembly, and phasing approaches that are utilized in these genome studies. Finally, we anticipate that accurate and telomere-to-telomere chromosome-level assembly of polyploid genomes could ultimately become a routine procedure in the near future.

New insights into the genome and transmission of the microsporidian pathogen Nosema muscidifuracis

Introduction

Nosema is a diverse genus of unicellular microsporidian parasites of insects and other arthropods. Nosema muscidifuracis infects parasitoid wasp species of Muscidifurax zaraptor and M. raptor (Hymenoptera: Pteromalidae), causing ~50% reduction in longevity and ~90% reduction in fecundity.

Methods and Results

Here, we report the first assembly of the N. muscidifuracis genome (14,397,169 bp in 28 contigs) of high continuity (contig N50 544.3 Kb) and completeness (BUSCO score 97.0%). A total of 2,782 protein-coding genes were annotated, with 66.2% of the genes having two copies and 24.0% of genes having three copies. These duplicated genes are highly similar, with a sequence identity of 99.3%. The complex pattern suggests extensive gene duplications and rearrangements across the genome. We annotated 57 rDNA loci, which are highly GC-rich (37%) in a GC-poor genome (25% genome average). Nosema-specific qPCR primer sets were designed based on 18S rDNA annotation as a diagnostic tool to determine its titer in host samples. We discovered high Nosema titers in Nosema-cured M. raptor and M. zaraptor using heat treatment in 2017 and 2019, suggesting that the remedy did not completely eliminate the Nosema infection. Cytogenetic analyses revealed heavy infections of N. muscidifuracis within the ovaries of M. raptor and M. zaraptor, consistent with the titer determined by qPCR and suggesting a heritable component of infection and per ovum vertical transmission.

Discussion

The parasitoids-Nosema system is laboratory tractable and, therefore, can serve as a model to inform future genome manipulations of Nosema-host system for investigations of Nosemosis.

Where whole-genome duplication is most beneficial: Adaptation of mangroves to a wide salinity range between land and sea

Whole-genome duplication (WGD) is believed to increase the chance of adaptation to a new environment. This conjecture may apply particularly well to new environments that are not only different but also more variable than ancestral habitats. One such prominent environment is the interface between land and sea, which has been invaded by woody plants, collectively referred as mangroves, multiple times. Here, we use two distantly related mangrove species (Avicennia marina and Rhizophora apiculata) to explore the effects of WGD on the adaptive process. We found that a high proportion of duplicated genes retained after WGD have acquired derived differential expression in response to salt gradient treatment. The WGD duplicates differentially expressed in at least one copy usually (>90%) diverge from their paralogues' expression profiles. Furthermore, both species evolved in parallel to have one paralogue expressed at a high level in both fresh water and hypersaline conditions but at a lower level at medium salinity. The pattern contrasts with the conventional view of monotone increase/decrease as salinity increases. Differentially expressed copies have thus probably acquired a new role in salinity tolerance. Our results indicate that the WGD duplicates may have evolved to function collaboratively in coping with different salinity levels, rather than specializing in the intermediate salinity optimal for mangrove plants. In conclusion, WGD and the retained duplicates appear to be an effective solution for adaptation to new and unstable environments.