Gene retention, fractionation and subgenome differences in polyploid plants

A CRISPR/Cas9-induced male-sterile line facilitating easy hybrid production in polyploid rapeseed (Brassica napus)

Rapeseed is a globally significant oilseed crop cultivated to meet the increasing demand for vegetable oil. In order to enhance yield and sustainability, breeders have adopted the development of rapeseed hybrids as a common strategy. However, current hybrid production systems in rapeseed have various limitations, necessitating the development of a simpler and more efficient approach. In this study, we propose a novel method involving the targeted disruption of Defective in Anther Dehiscence1 of Brassica napus (BnDAD1), an essential gene in the jasmonic acid biosynthesis pathway, using CRISPR/Cas9 technology, to create male-sterile lines. BnDAD1 was found to be dominantly expressed in the stamen of rapeseed flower buds. Disrupting BnDAD1 led to decreased levels of α-linolenic acid and jasmonate in the double mutants, resulting in defects in anther dehiscence and pollen maturation. By crossing the double mutant male-sterile lines with male-fertile lines, a two-line system was demonstrated, enabling the production of F 1 seeds. The male-sterile trait of the bndad1 double mutant lines was maintainable by applying exogenous methyl jasmonate and subsequently self-pollinating the flowers. This breakthrough holds promising potential for harnessing heterosis in rapeseed and offers a simpler and more efficient method for producing hybrid seeds.

Complex petal spot formation in the Beetle Daisy (Gorteria diffusa) relies on spot-specific accumulation of malonylated anthocyanin regulated by paralogous GdMYBSG6 transcription factors

Gorteria diffusa has elaborate petal spots that attract pollinators through sexual deception, but how G. diffusa controls spot development is largely unknown. Here, we investigate how pigmentation is regulated during spot formation. We determined the anthocyanin composition of G. diffusa petals and combined gene expression analysis with protein interaction assays to characterise R2R3-MYBs that likely regulate pigment production in G. diffusa petal spots. We found that cyanidin 3-glucoside pigments G. diffusa ray floret petals. Unlike other petal regions, spots contain a high proportion of malonylated anthocyanin. We identified three subgroup 6 R2R3-MYB transcription factors (GdMYBSG6-1,2,3) that likely activate the production of spot pigmentation. These genes are upregulated in developing spots and induce ectopic anthocyanin production upon heterologous expression in tobacco. Interaction assays suggest that these transcription factors regulate genes encoding three anthocyanin synthesis enzymes. We demonstrate that the elaboration of complex spots in G. diffusa begins with the accumulation of malonylated pigments at the base of ray floret petals, positively regulated by three paralogous R2R3-MYB transcription factors. Our results indicate that the functional diversification of these GdMYBSG6s involved changes in the spatial control of their transcription, and modification of the duration of GdMYBSG6 gene expression contributes towards floral variation within the species.

Allopolyploid subgenome identification and implications for evolutionary analysis

Whole-genome duplications (WGDs) are widespread genomic events in eukaryotes that are hypothesized to contribute to the evolutionary success of many lineages, including flowering plants, Saccharomyces yeast, and vertebrates. WGDs generally can be classified into autopolyploids (ploidy increase descended from one species) or allopolyploids (ploidy increase descended from multiple species). Assignment of allopolyploid progenitor species (called subgenomes in the polyploid) is important to understanding the biology and evolution of polyploids, including the asymmetric subgenome evolution following hybridization (biased fractionation). Here, I review the different methodologies used to identify the ancestors of allopolyploid subgenomes, discuss the advantages and disadvantages of these methods, and outline the implications of how these methods affect the subsequent evolutionary analysis of these genomes.

CDC48 in plants and its emerging function in plant immunity

Protein homeostasis, namely the balance between protein synthesis and degradation, must be finely controlled to ensure cell survival, notably through the ubiquitin-proteasome system (UPS). In all species, including plants, homeostasis is disrupted by biotic and abiotic stresses. A key player in the maintenance of protein balance, the protein CDC48, shows emerging functions in plants, particularly in response to biotic stress. In this review on CDC48 in plants, we detail its highly conserved structure, describe a gene expansion that is only present in Viridiplantae, discuss its various functions and regulations, and finally highlight its recruitment, still not clear, during the plant immune response.

Genome evolution of the ancient hexaploid Platanus × acerifolia (London planetree)

Whole-genome duplication (WGD; i.e., polyploidy) and chromosomal rearrangement (i.e., genome shuffling) significantly influence genome structure and organization. Many polyploids show extensive genome shuffling relative to their pre-WGD ancestors. No reference genome is currently available for Platanaceae (Proteales), one of the sister groups to the core eudicots. Moreover, Platanus × acerifolia (London planetree; Platanaceae) is a widely used street tree. Given the pivotal phylogenetic position of Platanus and its 2-y flowering transition, understanding its flowering-time regulatory mechanism has significant evolutionary implications; however, the impact of Platanus genome evolution on flowering-time genes remains unknown. Here, we assembled a high-quality, chromosome-level reference genome for P. × acerifolia using a phylogeny-based subgenome phasing method. Comparative genomic analyses revealed that P. × acerifolia (2n = 42) is an ancient hexaploid with three subgenomes resulting from two sequential WGD events; Platanus does not seem to share any WGD with other Proteales or with core eudicots. Each P. × acerifolia subgenome is highly similar in structure and content to the reconstructed pre-WGD ancestral eudicot genome without chromosomal rearrangements. The P. × acerifolia genome exhibits karyotypic stasis and gene sub-/neo-functionalization and lacks subgenome dominance. The copy number of flowering-time genes in P. × acerifolia has undergone an expansion compared to other noncore eudicots, mainly via the WGD events. Sub-/neo-functionalization of duplicated genes provided the genetic basis underlying the unique flowering-time regulation in P. × acerifolia. The P. × acerifolia reference genome will greatly expand understanding of the evolution of genome organization, genetic diversity, and flowering-time regulation in angiosperms.

Characteristics of duplicated gene expression and DNA methylation regulation in different tissues of allopolyploid Brassica napus

Plant polyploidization increases the complexity of epigenomes and transcriptional regulation, resulting in genome evolution and enhanced adaptability. However, few studies have been conducted on the relationship between gene expression and epigenetic modification in different plant tissues after allopolyploidization. In this study, we studied gene expression and DNA methylation modification patterns in four tissues (stems, leaves, flowers and siliques) of Brassica napusand its diploid progenitors. On this basis, the alternative splicing patterns and cis-trans regulation patterns of four tissues in B. napus and its diploid progenitors were also analyzed. It can be seen that the number of alternative splicing occurs in the B. napus is higher than that in the diploid progenitors, and the IR type increases the most during allopolyploidy. In addition, we studied the fate changes of duplicated genes after allopolyploidization in B. napus. We found that the fate of most duplicated genes is conserved, but the number of neofunctionalization and specialization is also large. The genetic fate of B. napus was classified according to five replication types (WGD, PD, DSD, TD, TRD). This study also analyzed generational transmission analysis of expression and DNA methylation patterns. Our study provides a reference for the fate differentiation of duplicated genes during allopolyploidization.

Telomere-to-telomere reference genome for Panax ginseng highlights the evolution of saponin biosynthesis

Ginseng (Panax ginseng) is a representative of Chinese traditional medicine, also used worldwide, while the triterpene saponin ginsenoside is the most important effective compound within it. Ginseng is an allotetraploid, with complex genetic background, making the study of its metabolic evolution challenging. In this study, we assembled a telomere-to-telomere ginseng reference genome, constructed of 3.45 Gb with 24 chromosomes and 77 266 protein-coding genes. Additionally, the reference genome was divided into two subgenomes, designated as subgenome A and B. Subgenome A contains a larger number of genes, whereas subgenome B has a general expression advantage, suggesting that ginseng subgenomes experienced asymmetric gene loss with biased gene expression. The two subgenomes separated approximately 6.07 million years ago, and subgenome B shows the closest relation to Panax vietnamensis var. fuscidiscus. Comparative genomics revealed an expansion of gene families associated with ginsenoside biosynthesis in both ginseng subgenomes. Furthermore, both tandem duplications and proximal duplications play crucial roles in ginsenoside biosynthesis. We also screened functional genes identified in previous research and found that some of these genes located in colinear regions between subgenomes have divergence functions, revealing an unbalanced evolution in both subgenomes and the saponin biosynthesis pathway in ginseng. Our work provides important resources for future genetic studies and breeding programs of ginseng, as well as the biosynthesis of ginsenosides.

ATAD3 Proteins: Unique Mitochondrial Proteins Essential for Life in Diverse Eukaryotic Lineages

ATPase family AAA domain-containing 3 (ATAD3) proteins are unique mitochondrial proteins that arose deep in the eukaryotic lineage but that are surprisingly absent in Fungi and Amoebozoa. These ∼600-amino acid proteins are anchored in the inner mitochondrial membrane and are essential in metazoans and Arabidopsis thaliana. ATAD3s comprise a C-terminal ATPases Associated with a variety of cellular Activities (AAA+) matrix domain and an ATAD3_N domain, which is located primarily in the inner membrane space but potentially extends to the cytosol to interact with the ER. Sequence and structural alignments indicate that ATAD3 proteins are most similar to classic chaperone unfoldases in the AAA+ family, suggesting that they operate in mitochondrial protein quality control. A. thaliana has four ATAD3 genes in two distinct clades that appear first in the seed plants, and both clades are essential for viability. The four genes are generally coordinately expressed, and transcripts are highest in growing apices and imbibed seeds. Plants with disrupted ATAD3 have reduced growth, aberrant mitochondrial morphology, diffuse nucleoids and reduced oxidative phosphorylation complex I. These and other pleiotropic phenotypes are also observed in ATAD3 mutants in metazoans. Here, we discuss the distribution of ATAD3 proteins as they have evolved in the plant kingdom, their unique structure, what we know about their function in plants and the challenges in determining their essential roles in mitochondria.

The era of panomics-driven gene discovery in plants

Panomics is an approach to integrate multiple 'omics' datasets, generated using different individuals or natural variations. Considering their diverse phenotypic spectrum, the phenome is inherently associated with panomics-based science, which is further combined with genomics, transcriptomics, metabolomics, and other omics techniques, either independently or collectively. Panomics has been accelerated through recent technological advancements in the field of genomics that enable the detection of population-wide structural variations (SVs) and hence offer unprecedented insights into the genetic variations contributing to important agronomic traits. The present review provides the recent trends of panomics-driven gene discovery toward various traits related to plant development, stress tolerance, accumulation of specialized metabolites, and domestication/dedomestication. In addition, the success stories are highlighted in the broader context of enhancing crop productivity.

Phylogenomic profiles of whole-genome duplications in Poaceae and landscape of differential duplicate retention and losses among major Poaceae lineages

Poaceae members shared a whole-genome duplication called rho. However, little is known about the evolutionary pattern of the rho-derived duplicates among Poaceae lineages and implications in adaptive evolution. Here we present phylogenomic/phylotranscriptomic analyses of 363 grasses covering all 12 subfamilies and report nine previously unknown whole-genome duplications. Furthermore, duplications from a single whole-genome duplication were mapped to multiple nodes on the species phylogeny; a whole-genome duplication was likely shared by woody bamboos with possible gene flow from herbaceous bamboos; and recent paralogues of a tetraploid Oryza are implicated in tolerance of seawater submergence. Moreover, rho duplicates showing differential retention among subfamilies include those with functions in environmental adaptations or morphogenesis, including ACOT for aquatic environments (Oryzoideae), CK2β for cold responses (Pooideae), SPIRAL1 for rapid cell elongation (Bambusoideae), and PAI1 for drought/cold responses (Panicoideae). This study presents a Poaceae whole-genome duplication profile with evidence for multiple evolutionary mechanisms that contribute to gene retention and losses.

The evolution and functional divergence of FT-related genes in controlling flowering time in Brassica rapa ssp. rapa

KEY MESSAGE

The BrrFT paralogues exhibit distinct expression patterns and play different roles in regulating flowering time, and BrrFT4 competes with BrrFT1 and BrrFT2 to interact with BrrFD proteins. Flowering time is an important agricultural trait for Brassica crops, and early bolting strongly affects the yield and quality of Brassica rapa ssp. rapa. Flowering Locus T paralogues play an important role in regulating flowering time. In this study, we identified FT-related genes in turnip by phylogenetic classification, and four BrrFT homoeologs that shared with high identities with BraFT genes were isolated. The different gene structures, promoter binding sites, and expression patterns observed indicated that these genes may play different roles in flowering time regulation. Further genetic and biochemical experiments showed that as for FT-like paralogues, BrrFT2 acted as the key floral inducer, and BrrFT1 seems to act as a mild 'florigen' protein. However, BrrFT4 acts as a floral repressor and antagonistically regulates flowering time by competing with BrrFT1 and BrrFT2 to bind BrrFD proteins. BrrFT3 may have experienced loss of function via base shift mutation. Our results revealed the potential roles of FT-related genes in flowering time regulation in turnip.

Cytonuclear Interactions and Subgenome Dominance Shape the Evolution of Organelle-Targeted Genes in the Brassica Triangle of U

The interaction and coevolution between nuclear and cytoplasmic genomes are one of the fundamental hallmarks of eukaryotic genome evolution and, 2 billion yr later, are still major contributors to the formation of new species. Although many studies have investigated the role of cytonuclear interactions following allopolyploidization, the relative magnitude of the effect of subgenome dominance versus cytonuclear interaction on genome evolution remains unclear. The Brassica triangle of U features 3 diploid species that together have formed 3 separate allotetraploid species on similar evolutionary timescales, providing an ideal system for understanding the contribution of the cytoplasmic donor to hybrid polyploid. Here, we investigated the evolutionary pattern of organelle-targeted genes in Brassica carinata (BBCC) and 2 varieties of Brassica juncea (AABB) at the whole-genome level, with particular focus on cytonuclear enzyme complexes. We found partial evidence that plastid-targeted genes experience selection to match plastid genomes, but no obvious corresponding signal in mitochondria-targeted genes from these 2 separately formed allopolyploids. Interestingly, selection acting on plastid genomes always reduced the retention rate of plastid-targeted genes encoded by the B subgenome, regardless of whether the Brassica nigra (BB) subgenome was contributed by the paternal or maternal progenitor. More broadly, this study illustrates the distinct selective pressures experienced by plastid- and mitochondria-targeted genes, despite a shared pattern of inheritance and natural history. Our study also highlights an important role for subgenome dominance in allopolyploid genome evolution, even in genes whose function depends on separately inherited molecules.

Hybridization and gene expression: Beyond differentially expressed genes

Gene expression has a key role in reproductive isolation, and studies of hybrid gene expression have identified mechanisms causing hybrid sterility. Here, we review the evidence for altered gene expression following hybridization and outline the mechanisms shown to contribute to altered gene expression in hybrids. Transgressive gene expression, transcending that of both parental species, is pervasive in early generation sterile hybrids, but also frequently observed in viable, fertile hybrids. We highlight studies showing that hybridization can result in transgressive gene expression, also in established hybrid lineages or species. Such extreme patterns of gene expression in stabilized hybrid taxa suggest that altered hybrid gene expression may result in hybridization-derived evolutionary novelty. We also conclude that while patterns of misexpression in hybrids are well documented, the understanding of the mechanisms causing misexpression is lagging. We argue that jointly assessing differences in cell composition and cell-specific changes in gene expression in hybrids, in addition to assessing changes in chromatin and methylation, will significantly advance our understanding of the basis of altered gene expression. Moreover, uncovering to what extent evolution of gene expression results in altered expression for individual genes, or entire networks of genes, will advance our understanding of how selection moulds gene expression. Finally, we argue that jointly studying the dual roles of altered hybrid gene expression, serving both as a mechanism for reproductive isolation and as a substrate for hybrid ecological adaptation, will lead to significant advances in our understanding of the evolution of gene expression.

A chromosome-scale genome sequence of Aeonium(Aeonium arboreum 'Velour') provides novel insights into the evolution of anthocyanin synthesis

Anthocyanin glycoside is a water-soluble flavonoid compound that colors plants and aids in stress resistance. The mechanism driving the evolution of the anthocyanin synthesis pathway in plants remains unclear. Aeonium plants are highly regarded as model organisms for studying adaptive evolution. These plants can be categorized into various types, each distinguished by the content and distribution of anthocyanins in their leaves. The categories include red leaves, green leaves, black leaves, yellow leaves, and a classification known as the 'spot brocade series. In this study, we successfully assembled and annotated the genome of cultivar 'Aeonium arboreum 'Velour'' at chromosomal level. The genome size is 1,334.85 Mb containing 18 chromosomes in a single set, with a contig N50 of 23.47 Mb and a Scaffold N50 of 25.07 Mb. Through homology prediction, de novo prediction, and transcriptome prediction, we identified 166,228 coding genes, 161,656 of which were successfully annotated in the database. Comparative genomic analysis revealed that Aeonium arboreum 'Velour' underwent an independent genome-wide replication event after differentiating from Sedum album, Kalanchoe laxiflora, and Kalanchoe fedtschenkoi. It also shared a genome-wide replication event with Sedum album and Kalanchoe laxiflora. Aeonium arboreum 'Velour' exhibits a higher number of multi-copy gene families compared to other species. A total of 5,129 gene families unique to Aeonium arboreum 'Velour' were identified, primarily enriched in various metabolic pathways, including monoterpenoid biosynthesis, sesquiterpene and triterpene biosynthesis, cyanamide acid metabolism, flavonoid and flavonol biosynthesis, phosphonate and phosphinate metabolism, fatty acid degradation, biosynthesis of unsaturated fatty acid, ether lipid metabolism, tyrosine metabolism, and isoflavone biosynthesis according to the KEGG pathway analysis. Aeonium arboreum 'Velour' and Sedum album diversion dates back to approximately 43.11 million years ago during the Paleogene period, marked by the expansion of 2,807 gene families. In contrast, the divergence from Kalanchoe laxiflora and Kalanchoe fedtschenkoi began around 57.28 million years ago, with 219 gene families expanding. GO analysis highlighted that most of the expansion or contraction gene families were predominantly enriched in flower organs, leaf organ development, anthocyanin metabolism regulation, and light energy absorption and utilization. Remarkably, anthocyanin metabolism regulation is enriched to 80 expanded genes, including 36 bHLH transcription factors, possibly functioning as photosensitive pigment interaction factors (PIFs). We speculate that flavonoids play a pivotal role in the adaptation of Aeonium arboreum 'Velour' to environmental stress. Moreover, the evolution of the anthocyanin synthesis pathway is potentially driven by the plant's capability to absorb and utilize light energy, especially in high CO2 and high-temperature settings characteristic of the early Paleogene.

Unifying duplication episode clustering and gene-species mapping inference

We present a novel problem, called MetaEC, which aims to infer gene-species assignments in a collection of partially leaf-labeled gene trees labels by minimizing the size of duplication episode clustering (EC). This problem is particularly relevant in metagenomics, where incomplete data often poses a challenge in the accurate reconstruction of gene histories. To solve MetaEC, we propose a polynomial time dynamic programming (DP) formulation that verifies the existence of a set of duplication episodes from a predefined set of episode candidates. In addition, we design a method to infer distributions of gene-species mappings. We then demonstrate how to use DP to design an algorithm that solves MetaEC. Although the algorithm is exponential in the worst case, we introduce a heuristic modification of the algorithm that provides a solution with the knowledge that it is exact. To evaluate our method, we perform two computational experiments on simulated and empirical data containing whole genome duplication events, showing that our algorithm is able to accurately infer the corresponding events.

Phenotyping, genetics, and "-omics" approaches to unravel and introgress enhanced resistance against apple scab (Venturia inaequalis) in apple cultivars (Malus × domestica)

Apple scab disease, caused by the fungus Venturia inaequalis, endangers commercial apple production globally. It is predominantly managed by frequent fungicide sprays that can harm the environment and promote the development of fungicide-resistant strains. Cultivation of scab-resistant cultivars harboring diverse qualitative Rvi resistance loci and quantitative trait loci associated with scab resistance could reduce the chemical footprint. A comprehensive understanding of the host-pathogen interaction is, however, needed to efficiently breed cultivars with enhanced resistance against a variety of pathogenic strains. Breeding efforts should not only encompass pyramiding of Rvi loci and their corresponding resistance alleles that directly or indirectly recognize pathogen effectors, but should also integrate genes that contribute to effective downstream defense mechanisms. This review provides an overview of the phenotypic and genetic aspects of apple scab resistance, and currently known corresponding defense mechanisms. Implementation of recent "-omics" approaches has provided insights into the complex network of physiological, molecular, and signaling processes that occur before and upon scab infection, thereby revealing the importance of both constitutive and induced defense mechanisms. Based on the current knowledge, we outline advances toward more efficient introgression of enhanced scab resistance into novel apple cultivars by conventional breeding or genetic modification techniques. However, additional studies integrating different "-omics" approaches combined with functional studies will be necessary to unravel effective defense mechanisms as well as key regulatory genes underpinning scab resistance in apple. This crucial information will set the stage for successful knowledge-based breeding for enhanced scab resistance.

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.

Asymmetric and parallel subgenome selection co-shape common carp domestication

BACKGROUND

The common carp (Cyprinus carpio) might best represent the domesticated allopolyploid animals. Although subgenome divergence which is well-known to be a key to allopolyploid domestication has been comprehensively characterized in common carps, the link between genetic architecture underlying agronomic traits and subgenome divergence is unknown in the selective breeding of common carps globally.

RESULTS

We utilized a comprehensive SNP dataset in 13 representative common carp strains worldwide to detect genome-wide genetic variations associated with scale reduction, vibrant skin color, and high growth rate in common carp domestication. We identified numerous novel candidate genes underlie the three agronomically most desirable traits in domesticated common carps, providing potential molecular targets for future genetic improvement in the selective breeding of common carps. We found that independently selective breeding of the same agronomic trait (e.g., fast growing) in common carp domestication could result from completely different genetic variations, indicating the potential advantage of allopolyploid in domestication. We observed that candidate genes associated with scale reduction, vibrant skin color, and/or high growth rate are repeatedly enriched in the immune system, suggesting that domestication of common carps was often accompanied by the disease resistance improvement.

CONCLUSIONS

In common carp domestication, asymmetric subgenome selection is prevalent, while parallel subgenome selection occurs in selective breeding of common carps. This observation is not due to asymmetric gene retention/loss between subgenomes but might be better explained by reduced pleiotropy through transposable element-mediated expression divergence between ohnologs. Our results demonstrate that domestication benefits from polyploidy not only in plants but also in animals.

A genome assembly for Orinus kokonorica provides insights into the origin, adaptive evolution and further diversification of two closely related grass genera

Ancient whole-genome duplication (WGD) or polyploidization is prevalent in plants and has played a crucial role in plant adaptation. However, the underlying genomic basis of ecological adaptation and subsequent diversification after WGD are still poorly understood in most plants. Here, we report a chromosome-scale genome assembly for the genus Orinus (Orinus kokonorica as representative) and preform comparative genomics with its closely related genus Cleistogenes (Cleistogenes songorica as representative), both belonging to a newly named subtribe Orininae of the grass subfamily Chloridoideae. The two genera may share one paleo-allotetraploidy event before 10 million years ago, and the two subgenomes of O. kokonorica display neither fractionation bias nor global homoeolog expression dominance. We find substantial genome rearrangements and extensive structural variations (SVs) between the two species. With comparative transcriptomics, we demonstrate that functional innovations of orthologous genes may have played an important role in promoting adaptive evolution and diversification of the two genera after polyploidization. In addition, copy number variations and extensive SVs between orthologs of flower and rhizome related genes may contribute to the morphological differences between the two genera. Our results provide new insights into the adaptive evolution and subsequent diversification of the two genera after polyploidization.

Evolutionary analysis of conserved non-coding elements subsequent to whole-genome duplication in opium poppy

Whole-genome duplication (WGD) leads to the duplication of both coding and non-coding sequences within an organism's genome, providing an abundant supply of genetic material that can drive evolution, ultimately contributing to plant adaptation and speciation. Although non-coding sequences contain numerous regulatory elements, they have been understudied compared to coding sequences. In order to address this gap, we explored the evolutionary patterns of regulatory sequences, coding sequences and transcriptomes using conserved non-coding elements (CNEs) as regulatory element proxies following the recent WGD event in opium poppy (Papaver somniferum). Our results showed similar evolutionary patterns in subgenomes of regulatory and coding sequences. Specifically, the biased or unbiased retention of coding sequences reflected the same pattern as retention levels in regulatory sequences. Further, the divergence of gene expression patterns mediated by regulatory element variations occurred at a more rapid pace than that of gene coding sequences. However, gene losses were purportedly dependent on relaxed selection pressure in coding sequences. Specifically, the rapid evolution of tissue-specific benzylisoquinoline alkaloid production in P. somniferum was associated with regulatory element changes. The origin of a novel stem-specific ACR, which utilized ancestral cis-elements as templates, is likely to be linked to the evolutionary trajectory behind the transition of the PSMT1-CYP719A21 cluster from high levels of expression solely in P. rhoeas root tissue to its elevated expression in P. somniferum stem tissue. Our findings demonstrate that rapid regulatory element evolution can contribute to the emergence of new phenotypes and provide valuable insights into the high evolvability of regulatory elements.

Subgenome dominance shapes novel gene evolution in the decaploid pitcher plant Nepenthes gracilis

Subgenome dominance after whole-genome duplication generates distinction in gene number and expression at the level of chromosome sets, but it remains unclear how this process may be involved in evolutionary novelty. Here we generated a chromosome-scale genome assembly of the Asian pitcher plant Nepenthes gracilis to analyse how its novel traits (dioecy and carnivorous pitcher leaves) are linked to genomic evolution. We found a decaploid karyotype and a clear indication of subgenome dominance. A male-linked and pericentromerically located region on the putative sex chromosome was identified in a recessive subgenome and was found to harbour three transcription factors involved in flower and pollen development, including a likely neofunctionalized LEAFY duplicate. Transcriptomic and syntenic analyses of carnivory-related genes suggested that the paleopolyploidization events seeded genes that subsequently formed tandem clusters in recessive subgenomes with specific expression in the digestive zone of the pitcher, where specialized cells digest prey and absorb derived nutrients. A genome-scale analysis suggested that subgenome dominance likely contributed to evolutionary innovation by permitting recessive subgenomes to diversify functions of novel tissue-specific duplicates. Our results provide insight into how polyploidy can give rise to novel traits in divergent and successful high-ploidy lineages.

Subgenomic Stability of Progenitor Genomes During Repeated Allotetraploid Origins of the Same Grass Brachypodium hybridum

Both homeologous exchanges and homeologous expression bias are generally found in most allopolyploid species. Whether homeologous exchanges and homeologous expression bias differ between repeated allopolyploid speciation events from the same progenitor species remains unknown. Here, we detected a third independent and recent allotetraploid origin for the model grass Brachypodium hybridum. Our homeologous exchange with replacement analyses indicated the absence of significant homeologous exchanges in any of the three types of wild allotetraploids, supporting the integrity of their progenitor subgenomes and the immediate creation of the amphidiploids. Further homeologous expression bias tests did not uncover significant subgenomic dominance in different tissues and conditions of the allotetraploids. This suggests a balanced expression of homeologs under similar or dissimilar ecological conditions in their natural habitats. We observed that the density of transposons around genes was not associated with the initial establishment of subgenome dominance; rather, this feature is inherited from the progenitor genome. We found that drought response genes were highly induced in the two subgenomes, likely contributing to the local adaptation of this species to arid habitats in the third allotetraploid event. These findings provide evidence for the consistency of subgenomic stability of parental genomes across multiple allopolyploidization events that led to the same species at different periods. Our study emphasizes the importance of selecting closely related progenitor species genomes to accurately assess homeologous exchange with replacement in allopolyploids, thereby avoiding the detection of false homeologous exchanges when using less related progenitor species genomes.

LHP1-mediated epigenetic buffering of subgenome diversity and defense responses confers genome plasticity and adaptability in allopolyploid wheat

Polyploidization is a major driver of genome diversification and environmental adaptation. However, the merger of different genomes may result in genomic conflicts, raising a major question regarding how genetic diversity is interpreted and regulated to enable environmental plasticity. By analyzing the genome-wide binding of 191 trans-factors in allopolyploid wheat, we identified like heterochromatin protein 1 (LHP1) as a master regulator of subgenome-diversified genes. Transcriptomic and epigenomic analyses of LHP1 mutants reveal its role in buffering the expression of subgenome-diversified defense genes by controlling H3K27me3 homeostasis. Stripe rust infection releases latent subgenomic variations by eliminating H3K27me3-related repression. The simultaneous inactivation of LHP1 homoeologs by CRISPR-Cas9 confers robust stripe rust resistance in wheat seedlings. The conditional repression of subgenome-diversified defenses ensures developmental plasticity to external changes, while also promoting neutral-to-non-neutral selection transitions and adaptive evolution. These findings establish an LHP1-mediated buffering system at the intersection of genotypes, environments, and phenotypes in polyploid wheat. Manipulating the epigenetic buffering capacity offers a tool to harness cryptic subgenomic variations for crop improvement.

Evolution and diversification of the O-methyltransferase (OMT) gene family in Solanaceae

O-methyltransferases (OMTs) are a group of enzymes involved in several fundamental biological processes in plants, including lignin biosynthesis, pigmentation, and aroma production. Despite the intensive investigation of the role of OMTs in plant secondary metabolism, the evolution and diversification of this gene family in Solanaceae remain poorly understood. Here, we conducted a genome-wide survey of OMT genes in six Solanaceae species, reconstructing gene phylogenetic trees, predicting the potential involvement in biological processes, and investigating the exon/intron structure and chromosomal location. We identified 57 caffeoyl-CoA OMTs (CCoAOMTs) and 196 caffeic acid OMTs (COMTs) in the studied species. We observed a conserved gene block on chromosome 2 that consisted of tandem duplicated copies of OMT genes. Our results suggest that the expansion of the OMT gene family in Solanaceae was driven by whole genome duplication, segmental duplication, and tandem duplication, with multiple genes being retained by neofunctionalization and subfunctionalization. This study represents an essential first step in unraveling the evolutionary history of OMTs in Solanaceae. Our findings deepen our understanding of the crucial role of OMTs in several biological processes and highlight their significance as potential biotechnological targets.

Physiological Control and Genetic Basis of Leaf Curvature and Heading in Brassica rapa L

BACKGROUND

Heading is an important agronomic feature for Chinese cabbage, cabbage, and lettuce. The heading leaves function as nutrition storage organs, which contribute to the high quality and economic worth of leafy heads. Leaf development is crucial during the heading stage, most genes previously predicted to be involved in the heading process are based on Arabidopsis leaf development studies.

AIM OF REVIEW

Till date, there is no published review article that demonstrated a complete layout of all the identified regulators of leaf curvature and heading. In this review, we have summarized all the identified physiological and genetic regulators that are directly or indirectly involved in leaf curvature and heading in Brassica crops. By integrating all identified regulators that provide a coherent logic of leaf incurvature and heading, we proposed a molecular mechanism in Brassica crops with graphical illustrations. This review adds value to future breeding of distinct heading kinds of cabbage and Chinese cabbage by providing unique insights into leaf development.

KEY SCIENTIFIC CONCEPTS OF REVIEW

Leaf curvature and heading are established by synergistic interactions among genes, transcription factors, microRNAs, phytohormones, and environmental stimuli that regulate primary and secondary morphogenesis. Various genes have been identified using transformation and genome editing that are responsible for the formation of leaf curvature and heading in Brassica crops. A range of leaf morphologies have been observed in Brassica, which are established because of the mutated determinants that are responsible for cell division and leaf polarity.

Systematic thermal analysis of the Arabidopsis proteome: Thermal tolerance, organization, and evolution

Understanding the thermal stability of the plant proteome in the context of the native cellular environment would aid the design of crops with high thermal tolerance, but only limited such data are available. Here, we applied quantitative mass spectrometry to profile the thermal stability of the Arabidopsis proteome and identify thermo-sensitive and thermo-resilient protein networks in Arabidopsis, providing a basis for understanding heat-induced damage. We also show that the similarities of the protein-melting curves can be used as a proxy to evaluate system-wide protein-protein interactions in non-engineered plants and enable the identification of transient interactions exhibited by metabolons in the context of the cellular environment. Finally, we report a systematic comparison of the thermal stability of paralogs in Arabidopsis to aid the investigation and understanding of gene duplication and protein evolution. Taken together, our results could have broad implications for the fields of plant thermal tolerance, plant protein assemblies, and evolution.

The lack of negative association between TE load and subgenome dominance in synthesized Brassica allotetraploids

Polyploidization is important to the evolution of plants. Subgenome dominance is a distinct phenomenon associated with most allopolyploids. A gene on the dominant subgenome tends to express to higher RNA levels in all organs as compared to the expression of its syntenic paralogue (homoeolog). The mechanism that underlies the formation of subgenome dominance remains unknown, but there is evidence for the involvement of transposon/DNA methylation density differences nearby the genes of parents as being causal. The subgenome with lower density of transposon and methylation near genes is positively associated with subgenome dominance. Here, we generated eight generations of allotetraploid progenies from the merging of parental genomes Brassica rapa and Brassica oleracea. We found that transposon/methylation density differ near genes between the parental (rapa:oleracea) existed in the wide hybrid, persisted in the neotetraploids (the synthetic Brassica napus), but these neotetraploids expressed no expected subgenome dominance. This absence of B. rapa vs. B. oleracea subgenome dominance is particularly significant because, while there is no negative relationship between transposon/methylation level and subgenome dominance in the neotetraploids, the more ancient parental subgenomes for all Brassica did show differences in transposon/methylation densities near genes and did express, in the same samples of cells, biased gene expression diagnostic of subgenome dominance. We conclude that subgenome differences in methylated transposon near genes are not sufficient to initiate the biased gene expressions defining subgenome dominance. Our result was unexpected, and we suggest a "nuclear chimera" model to explain our data.

From genome size to trait evolution during angiosperm radiation

Angiosperm diversity arises from trait flexibility and repeated evolutionary radiations, but the role of genomic characters in these radiations remains unclear. In this opinion article, we discuss how genome size can influence angiosperm diversification via its intricate link with cell size, tissue packing, and physiological processes which, in turn, influence the macroevolution of functional traits. We propose that integrating genome size, functional traits, and phylogenetic data across a wide range of lineages allows us to test whether genome size decrease consistently leads to increased trait flexibility, while genome size increase constrains trait evolution. Combining theories from molecular biology, functional ecology and macroevolution, we provide a framework to better understand the role of genome size in trait evolution, evolutionary radiations, and the global distribution of angiosperms.

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

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

Genome shock in a synthetic allotetraploid wheat invokes subgenome-partitioned gene regulation, meiotic instability, and karyotype variation

It is becoming increasingly evident that interspecific hybridization at the homoploid level or coupled with whole-genome duplication (i.e. allopolyploidization) has played a major role in biological evolution. However, the direct impacts of hybridization and allopolyploidization on genome structure and function, phenotype, and fitness remains to be fully understood. Synthetic hybrids and allopolyploids are trackable experimental systems that can be used to address this issue. In this study, we resynthesized a pair of reciprocal F1 hybrids and corresponding reciprocal allotetraploids using the two diploid progenitor species of bread wheat (Triticum aestivum, BBAADD), namely T. urartu (AA) and Aegilops tauschii (DD). By comparing phenotypes related to growth, development, and fitness, and by analysing genome expression in both hybrids and allotetraploids in relation to the parents, we found that the types and trends of karyotype variation in the immediately formed allotetraploids were correlated with both instability of meiosis and chromosome- and subgenome-biased expression. We determined clear advantages of allotetraploids over diploid F1 hybrids in several morphological traits including fitness that mirrored the tissue- and developmental stage-dependent subgenome-partitioning of the allotetraploids. The allotetraploids were meiotically unstable primarily due to homoeologous pairing that varied dramatically among the chromosomes. Nonetheless, the manifestation of organismal karyotype variation and the occurrence of meiotic irregularity were not concordant, suggesting a role of functional constraints probably imposed by subgenome- and chromosome-biased gene expression. Our results provide new insights into the direct impacts and consequences of hybridization and allopolyploidization that are relevant to evolution and likely to be informative for future crop improvement approaches using synthetic polyploids.

DNA methylation role in subgenome expression dominance of Juglans regia and its wild relative J. mandshurica

Subgenome expression dominance plays a crucial role in the environmental adaptation of polyploids. However, the epigenetic molecular mechanism underlying this process has not been thoroughly investigated, particularly in perennial woody plants. Persian walnut (Juglans regia) and its wild relative, Manchurian walnut (Juglans mandshurica), are woody plants of great economic importance and are both paleopolyploids that have undergone whole-genome duplication events. In this study, we explored the characteristics of subgenome expression dominance in these 2 Juglans species and examined its epigenetic basis. We divided their genomes into dominant subgenome (DS) and submissive subgenome (SS) and found that the DS-specific genes might play critical roles in biotic stress response or pathogen defense. We comprehensively elucidated the characteristics of biased gene expression, asymmetric DNA methylation, transposable elements (TEs), and alternative splicing (AS) events of homoeologous gene pairs between subgenomes. The results showed that biased expression genes (BEGs) in 2 Juglans species were mainly related to external stimuli response, while non-BEGs were related to complexes that might be involved in signal transduction. DS genes had higher expression and more AS events while having less DNA methylation and TEs than homoeologous genes from the SS in the 2 Juglans species. Further studies showed that DNA methylation might contribute to the biased expression of gene pairs by modifying LTR/TIR/nonTIR TEs and improving the AS efficiency of corresponding precursor mRNAs in a particular context. Our study contributes to understanding the epigenetic basis of subgenome expression dominance and the environmental adaptation of perennial woody plants.

New insights on the evolution of nucleolar dominance in newly resynthesized hexaploid wheat Triticum zhukovskyi

Nucleolar dominance (ND) is a widespread epigenetic phenomenon in hybridizations where nucleolus transcription fails at the nucleolus organizer region (NOR). However, the dynamics of NORs during the formation of Triticum zhukovskyi (GGAu Au Am Am ), another evolutionary branch of allohexaploid wheat, remains poorly understood. Here, we elucidated genetic and epigenetic changes occurring at the NOR loci within the Am , G, and D subgenomes during allopolyploidization by synthesizing hexaploid wheat GGAu Au Am Am and GGAu Au DD. In T. zhukovskyi, Au genome NORs from T. timopheevii (GGAu Au ) were lost, while the second incoming NORs from T. monococcum (Am Am ) were retained. Analysis of the synthesized T. zhukovskyi revealed that rRNA genes from the Am genome were silenced in F1 hybrids (GAu Am ) and remained inactive after genome doubling and subsequent self-pollinations. We observed increased DNA methylation accompanying the inactivation of NORs in the Am genome and found that silencing of NORs in the S1 generation could be reversed by a cytidine methylase inhibitor. Our findings provide insights into the ND process during the evolutionary period of T. zhukovskyi and highlight that inactive rDNA units may serve as a 'first reserve' in the form of R-loops, contributing to the successful evolution of T. zhukovskyi.

Genome assembly of Polygala tenuifolia provides insights into its karyotype evolution and triterpenoid saponin biosynthesis

Polygala tenuifolia is a perennial medicinal plant that has been widely used in traditional Chinese medicine for treating mental diseases. However, the lack of genomic resources limits the insight into its evolutionary and biological characterization. In the present work, we reported the P. tenuifolia genome, the first genome assembly of the Polygalaceae family. We sequenced and assembled this genome by a combination of Illumnina, PacBio HiFi, and Hi-C mapping. The assembly includes 19 pseudochromosomes covering ~92.68% of the assembled genome (~769.62 Mb). There are 36 463 protein-coding genes annotated in this genome. Detailed comparative genome analysis revealed that P. tenuifolia experienced two rounds of whole genome duplication that occurred ~39-44 and ~18-20 million years ago, respectively. Accordingly, we systematically reconstructed ancestral chromosomes of P. tenuifolia and inferred its chromosome evolution trajectories from the common ancestor of core eudicots to the present species. Based on the transcriptomics data, enzyme genes and transcription factors involved in the synthesis of triterpenoid saponin in P. tenuifolia were identified. Further analysis demonstrated that whole-genome duplications and tandem duplications play critical roles in the expansion of P450 and UGT gene families, which contributed to the synthesis of triterpenoid saponins. The genome and transcriptome data will not only provide valuable resources for comparative and functional genomic researches on Polygalaceae, but also shed light on the synthesis of triterpenoid saponin.

Genome-wide identification, classification, and expression analysis of heat shock transcription factor family in switchgrass (Panicum virgatum L.)

Switchgrass is one of the most promising bioenergy crops and is generally cultivated in arid climates and poor soils. Heat shock transcription factors (Hsfs) are key regulators of plant responses to abiotic and biotic stressors. However, their role and mechanism of action in switchgrass have not been elucidated. Hence, this study aimed to identify the Hsf family in switchgrass and understand its functional role in heat stress signal transduction and heat tolerance by using bioinformatics and RT-PCR analysis. Forty-eight PvHsfs were identified and divided into three main classes based on their gene structure and phylogenetic relationships: HsfA, HsfB, and HsfC. The results of the bioinformatics analysis showed a DNA-binding domain (DBD) at the N-terminal in PvHsfs, and they were not evenly distributed on all chromosomes except for chromosomes 8 N and 8 K. Many cis-elements related to plant development, stress responses, and plant hormones were identified in the promoter sequence of each PvHsf. Segmental duplication is the primary force underlying Hsf family expansion in switchgrass. The results of the expression pattern of PvHsfs in response to heat stress showed that PvHsf03 and PvHsf25 might play critical roles in the early and late stages of switchgrass response to heat stress, respectively, and HsfB mainly showed a negative response to heat stress. Ectopic expression of PvHsf03 in Arabidopsis significantly increased the heat resistance of seedlings. Overall, our research lays a notable foundation for studying the regulatory network in response to deleterious environments and for further excavating tolerance genes in switchgrass.

Whole-genome Duplications and the Long-term Evolution of Gene Regulatory Networks in Angiosperms

Angiosperms have a complex history of whole-genome duplications (WGDs), with varying numbers and ages of WGD events across clades. These WGDs have greatly affected the composition of plant genomes due to the biased retention of genes belonging to certain functional categories following their duplication. In particular, regulatory genes and genes encoding proteins that act in multiprotein complexes have been retained in excess following WGD. Here, we inferred protein-protein interaction (PPI) networks and gene regulatory networks (GRNs) for seven well-characterized angiosperm species and explored the impact of both WGD and small-scale duplications (SSDs) in network topology by analyzing changes in frequency of network motifs. We found that PPI networks are enriched in WGD-derived genes associated with dosage-sensitive intricate systems, and strong selection pressures constrain the divergence of WGD-derived genes at the sequence and PPI levels. WGD-derived genes in network motifs are mostly associated with dosage-sensitive processes, such as regulation of transcription and cell cycle, translation, photosynthesis, and carbon metabolism, whereas SSD-derived genes in motifs are associated with response to biotic and abiotic stress. Recent polyploids have higher motif frequencies than ancient polyploids, whereas WGD-derived network motifs tend to be disrupted on the longer term. Our findings demonstrate that both WGD and SSD have contributed to the evolution of angiosperm GRNs, but in different ways, with WGD events likely having a more significant impact on the short-term evolution of polyploids.

Homoeologous evolution of the allotetraploid genome of Poa annua L

BACKGROUND

Poa annua (annual bluegrass) is an allotetraploid turfgrass, an agronomically significant weed, and one of the most widely dispersed plant species on earth. Here, we report the chromosome-scale genome assemblies of P. annua's diploid progenitors, P. infirma and P. supina, and use multi-omic analyses spanning all three species to better understand P. annua's evolutionary novelty.

RESULTS

We find that the diploids diverged from their common ancestor 5.5 - 6.3 million years ago and hybridized to form P. annua ≤ 50,000 years ago. The diploid genomes are similar in chromosome structure and most notably distinguished by the divergent evolutionary histories of their transposable elements, leading to a 1.7 × difference in genome size. In allotetraploid P. annua, we find biased movement of retrotransposons from the larger (A) subgenome to the smaller (B) subgenome. We show that P. annua's B subgenome is preferentially accumulating genes and that its genes are more highly expressed. Whole-genome resequencing of several additional P. annua accessions revealed large-scale chromosomal rearrangements characterized by extensive TE-downsizing and evidence to support the Genome Balance Hypothesis.

CONCLUSIONS

The divergent evolutions of the diploid progenitors played a central role in conferring onto P. annua its remarkable phenotypic plasticity. We find that plant genes (guided by selection and drift) and transposable elements (mostly guided by host immunity) each respond to polyploidy in unique ways and that P. annua uses whole-genome duplication to purge highly parasitized heterochromatic sequences. The findings and genomic resources presented here will enable the development of homoeolog-specific markers for accelerated weed science and turfgrass breeding.

Diploid and tetraploid genomes of Acorus and the evolution of monocots

Monocots are a major taxon within flowering plants, have unique morphological traits, and show an extraordinary diversity in lifestyle. To improve our understanding of monocot origin and evolution, we generate chromosome-level reference genomes of the diploid Acorus gramineus and the tetraploid Ac. calamus, the only two accepted species from the family Acoraceae, which form a sister lineage to all other monocots. Comparing the genomes of Ac. gramineus and Ac. calamus, we suggest that Ac. gramineus is not a potential diploid progenitor of Ac. calamus, and Ac. calamus is an allotetraploid with two subgenomes A, and B, presenting asymmetric evolution and B subgenome dominance. Both the diploid genome of Ac. gramineus and the subgenomes A and B of Ac. calamus show clear evidence of whole-genome duplication (WGD), but Acoraceae does not seem to share an older WGD that is shared by most other monocots. We reconstruct an ancestral monocot karyotype and gene toolkit, and discuss scenarios that explain the complex history of the Acorus genome. Our analyses show that the ancestors of monocots exhibit mosaic genomic features, likely important for that appeared in early monocot evolution, providing fundamental insights into the origin, evolution, and diversification of monocots.

Characterization and Coexpression Analysis of the TIFY Family Genes in Euryale ferox Related to Leaf Development

TIFYs are plant-specific transcription factors that contain the TIFY structural domain and play an important role in plant leaf growth and development. However, the role played by TIFY in E. ferox (Euryale ferox Salisb.) leaf development has not been investigated. In this study, 23 TIFY genes were identified in E. ferox. Phylogenetic analyses of the TIFY genes showed clustering into three groups (JAZ, ZIM, and PPD). The TIFY domain was shown to be conserved. JAZ was mainly expanded via wholegenome triplication (WGT) in E. ferox. Based on analyses of the TIFY genes in nine species, we found that JAZ has a closer relationship with PPD, in addition to appearing the most recently and expanding most rapidly, leading to the rapid expansion of TIFYs in Nymphaeaceae. In addition, their different evolution types were discovered. Different gene expressions showed the distinct and corresponsive expression patterns of the EfTIFYs in different stages of tissue and leaf development. Finally, The qPCR analysis revealed that the expression of EfTIFY7.2 and EfTIFY10.1 showed an upward trend and high expression throughout leaf development. Further co-expression analysis indicated that EfTIFY7.2 might be more important for the development of E. ferox leaves. This information will be valuable when exploring the molecular mechanisms of EfTIFYs in plants.

An Emerging Animal Model for Querying the Role of Whole Genome Duplication in Development, Evolution, and Disease

Whole genome duplication (WGD) or polyploidization can occur at the cellular, tissue, and organismal levels. At the cellular level, tetraploidization has been proposed as a driver of aneuploidy and genome instability and correlates strongly with cancer progression, metastasis, and the development of drug resistance. WGD is also a key developmental strategy for regulating cell size, metabolism, and cellular function. In specific tissues, WGD is involved in normal development (e.g., organogenesis), tissue homeostasis, wound healing, and regeneration. At the organismal level, WGD propels evolutionary processes such as adaptation, speciation, and crop domestication. An essential strategy to further our understanding of the mechanisms promoting WGD and its effects is to compare isogenic strains that differ only in their ploidy. Caenorhabditis elegans (C. elegans) is emerging as an animal model for these comparisons, in part because relatively stable and fertile tetraploid strains can be produced rapidly from nearly any diploid strain. Here, we review the use of Caenorhabditis polyploids as tools to understand important developmental processes (e.g., sex determination, dosage compensation, and allometric relationships) and cellular processes (e.g., cell cycle regulation and chromosome dynamics during meiosis). We also discuss how the unique characteristics of the C. elegans WGD model will enable significant advances in our understanding of the mechanisms of polyploidization and its role in development and disease.

Two-step model of paleohexaploidy, ancestral genome reshuffling and plasticity of heat shock response in Asteraceae

An ancient hexaploidization event in the most but not all Asteraceae plants, may have been responsible for shaping the genomes of many horticultural, ornamental, and medicinal plants that promoting the prosperity of the largest angiosperm family on the earth. However, the duplication process of this hexaploidy, as well as the genomic and phenotypic diversity of extant Asteraceae plants caused by paleogenome reorganization, are still poorly understood. We analyzed 11 genomes from 10 genera in Asteraceae, and redated the Asteraceae common hexaploidization (ACH) event ~70.7-78.6 million years ago (Mya) and the Asteroideae specific tetraploidization (AST) event ~41.6-46.2 Mya. Moreover, we identified the genomic homologies generated from the ACH, AST and speciation events, and constructed a multiple genome alignment framework for Asteraceae. Subsequently, we revealed biased fractionations between the paleopolyploidization produced subgenomes, suggesting the ACH and AST both are allopolyplodization events. Interestingly, the paleochromosome reshuffling traces provided clear evidence for the two-step duplications of ACH event in Asteraceae. Furthermore, we reconstructed ancestral Asteraceae karyotype (AAK) that has 9 paleochromosomes, and revealed a highly flexible reshuffling of Asteraceae paleogenome. Of specific significance, we explored the genetic diversity of Heat Shock Transcription Factors (Hsfs) associated with recursive whole-genome polyploidizations, gene duplications, and paleogenome reshuffling, and revealed that the expansion of Hsfs gene families enable heat shock plasticity during the genome evolution of Asteraceae. Our study provides insights on polyploidy and paleogenome remodeling for the successful establishment of Asteraceae, and is helpful for further communication and exploration of the diversification of plant families and phenotypes.

An efficient approach to synthesize sterile allopolyploids through the combined reproduction mode of ameiotic oogenesis and sperm-egg fusion in the polyploid Carassius complex

The association between polyploidy and reproduction transition, which is an intriguing issue in evolutionary genetics, can also be exploited as an approach for genetic improvement in agriculture. Recently, we generated novel amphitriploids (NA3n) by integrating the genomes of the gynogenetic Carassius gibelio and sexual C. auratus, and found gynogenesis was recovered in most NA3n females (NA3n♀I). Here, we discovered a unique reproduction mode, termed ameio-fusiongenesis, which combines the abilities of both ameiotic oogenesis and sperm-egg fusion, in a few NA3n females (NA3n♀II). These females inherited ameiotic oogenesis to produce unreduced eggs from gynogenetic C. gibelio and sperm-egg fusion from sexual C. auratus. Subsequently, we utilized this unique reproduction mode to generate a group of synthetic alloheptaploids by crossing NA3n♀II with Megalobrama amblycephala. They contained all chromosomes of maternal NA3n♀II and a chromosomal set of paternal M. amblycephala. Intergenomic chromosome translocations between NA3n♀II and M. amblycephala were also observed in a few somatic cells. Primary oocytes of the alloheptaploid underwent severe apoptosis owing to incomplete double-strand break repair at prophase I. Although spermatocytes displayed similar chromosome behavior at prophase I, they underwent apoptosis due to chromosome separation failure at metaphase I. Therefore, the alloheptaploid females and males were all sterile. Finally, we established a sustainable clone for the large-scale production of NA3n♀II and developed an efficient approach to synthesize diverse allopolyploids containing genomes of different cyprinid species. These findings not only broaden our understanding of reproduction transition but also offer a practical strategy for polyploidy breeding and heterosis fixing.

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

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

DNA methylation signatures of duplicate gene evolution in angiosperms

Gene duplication is a source of evolutionary novelty. DNA methylation may play a role in the evolution of duplicate genes (paralogs) through its association with gene expression. While this relationship has been examined to varying extents in a few individual species, the generalizability of these results at either a broad phylogenetic scale with species of differing duplication histories or across a population remains unknown. We applied a comparative epigenomics approach to 43 angiosperm species across the phylogeny and a population of 928 Arabidopsis (Arabidopsis thaliana) accessions, examining the association of DNA methylation with paralog evolution. Genic DNA methylation was differentially associated with duplication type, the age of duplication, sequence evolution, and gene expression. Whole genome duplicates were typically enriched for CG-only gene-body methylated or unmethylated genes, while single-gene duplications were typically enriched for non-CG methylated or unmethylated genes. Non-CG methylation, in particular, was characteristic of more recent single-gene duplicates. Core angiosperm gene families differentiated into those which preferentially retain paralogs and 'duplication-resistant' families, which convergently reverted to singletons following duplication. Duplication-resistant families that still have paralogous copies were, uncharacteristically for core angiosperm genes, enriched for non-CG methylation. Non-CG methylated paralogs had higher rates of sequence evolution, higher frequency of presence-absence variation, and more limited expression. This suggests that silencing by non-CG methylation may be important to maintaining dosage following duplication and be a precursor to fractionation. Our results indicate that genic methylation marks differing evolutionary trajectories and fates between paralogous genes and have a role in maintaining dosage following duplication.

Genome-wide characterization of ubiquitin-conjugating enzyme gene family explores its genetic effects on the oil content and yield of Brassica napus

Ubiquitin-conjugating enzyme (UBC) is a critical part of the ubiquitin-proteasome pathway and plays crucial roles in growth, development and abiotic stress response in plants. Although UBC genes have been detected in several plant species, characterization of this gene family at the whole-genome level has not been conducted in Brassica napus. In the present study, 200 putative BnUBCs were identified in B. napus, which were clustered into 18 subgroups based on phylogenetic analysis. BnUBCs within each subgroup showed relatively conserved gene architectures and motifs. Moreover, the gene expression patterns in various tissues as well as the identification of cis-acting regulatory elements in BnUBC promoters suggested further investigation of their potential functions in plant growth and development. Furthermore, three BnUBCs were predicted as candidate genes for regulating agronomic traits related to oil content and yield through association mapping. In conclusion, this study provided a wealth of information on the UBC family in B. napus and revealed their effects on oil content and yield, which will aid future functional research and genetic breeding of B. napus.

Genome-wide analysis of NBS-LRR genes revealed contribution of disease resistance from Saccharum spontaneum to modern sugarcane cultivar

INTRODUCTION

During plant evolution, nucleotide-binding sites (NBS) and leucine-rich repeat (LRR) genes have made significant contributions to plant disease resistance. With many high-quality plant genomes sequenced, identification and comprehensive analyses of NBS-LRR genes at whole genome level are of great importance to understand and utilize them.

METHODS

In this study, we identified the NBS-LRR genes of 23 representative species at whole genome level, and researches on NBS-LRR genes of four monocotyledonous grass species, Saccharum spontaneum, Saccharum officinarum, Sorghum bicolor and Miscanthus sinensis, were focused.

RESULTS AND DISCUSSION

We found that whole genome duplication, gene expansion, and allele loss could be factors affecting the number of NBS-LRR genes in the species, and whole genome duplication is likely to be the main cause of the number of NBS-LRR genes in sugarcane. Meanwhile, we also found a progressive trend of positive selection on NBS-LRR genes. These studies further elucidated the evolutionary pattern of NBS-LRR genes in plants. Transcriptome data from multiple sugarcane diseases revealed that more differentially expressed NBS-LRR genes were derived from S. spontaneum than from S. officinarum in modern sugarcane cultivars, and the proportion was significantly higher than the expected. This finding reveals that S. spontaneum has a greater contribution to disease resistance for modern sugarcane cultivars. In addition, we observed allelespecific expression of seven NBS-LRR genes under leaf scald, and 125 NBS-LRR genes responding to multiple diseases were identified. Finally, we built a plant NBS-LRR gene database to facilitate subsequent analysis and use of NBSLRR genes obtained here. In conclusion, this study complemented and completed the research of plant NBS-LRR genes, and discussed how NBS-LRR genes responding to sugarcane diseases, which provided a guide and genetic resources for further research and utilization of NBS-LRR genes.

Analyses of Cullin1 homologs reveal functional redundancy in S-RNase-based self-incompatibility and evolutionary relationships in eudicots

In Petunia (Solanaceae family), self-incompatibility (SI) is regulated by the polymorphic S-locus, which contains the pistil-specific S-RNase and multiple pollen-specific S-Locus F-box (SLF) genes. SLFs assemble into E3 ubiquitin ligase complexes known as Skp1-Cullin1-F-box complexes (SCFSLF). In pollen tubes, these complexes collectively mediate ubiquitination and degradation of all nonself S-RNases, but not self S-RNase, resulting in cross-compatible, but self-incompatible, pollination. Using Petunia inflata, we show that two pollen-expressed Cullin1 (CUL1) proteins, PiCUL1-P and PiCUL1-B, function redundantly in SI. This redundancy is lost in Petunia hybrida, not because of the inability of PhCUL1-B to interact with SSK1, but due to a reduction in the PhCUL1-B transcript level. This is possibly caused by the presence of a DNA transposon in the PhCUL1-B promoter region, which was inherited from Petunia axillaris, one of the parental species of Pe. hybrida. Phylogenetic and syntenic analyses of Cullin genes in various eudicots show that three Solanaceae-specific CUL1 genes share a common origin, with CUL1-P dedicated to S-RNase-related reproductive processes. However, CUL1-B is a dispersed duplicate of CUL1-P present only in Petunia, and not in the other species of the Solanaceae family examined. We suggest that the CUL1s involved (or potentially involved) in the SI response in eudicots share a common origin.

Trends in the evolution of intronless genes in Poaceae

Intronless genes (IGs), which are a feature of prokaryotes, are a fascinating group of genes that are also present in eukaryotes. In the current study, a comparison of Poaceae genomes revealed that the origin of IGs may have involved ancient intronic splicing, reverse transcription, and retrotranspositions. Additionally, IGs exhibit the typical features of rapid evolution, including recent duplications, variable copy numbers, low divergence between paralogs, and high non-synonymous to synonymous substitution ratios. By tracing IG families along the phylogenetic tree, we determined that the evolutionary dynamics of IGs differed among Poaceae subfamilies. IG families developed rapidly before the divergence of Pooideae and Oryzoideae and expanded slowly after the divergence. In contrast, they emerged gradually and consistently in the Chloridoideae and Panicoideae clades during evolution. Furthermore, IGs are expressed at low levels. Under relaxed selection pressure, retrotranspositions, intron loss, and gene duplications and conversions may promote the evolution of IGs. The comprehensive characterization of IGs is critical for in-depth studies on intron functions and evolution as well as for assessing the importance of introns in eukaryotes.

Advances in the Mining of Disease Resistance Genes from Aegilops tauschii and the Utilization in Wheat

Aegilops tauschii is one of the malignant weeds that affect wheat production and is also the wild species ancestor of the D genome of hexaploid wheat (Triticum aestivum, AABBDD). It contains many disease resistance genes that have been lost in the long-term evolution of wheat and is an important genetic resource for the mining and utilization of wheat disease resistance genes. In recent years, the genome sequence of Aegilops tauschii has been preliminarily completed, which has laid a good foundation for the further exploration of wheat disease resistance genes in Aegilops tauschii. There are many studies on disease resistance genes in Aegilops tauschii; in order to provide better help for the disease resistance breeding of wheat, this paper analyzes and reviews the relationship between Aegilops tauschii and wheat, the research progress of Aegilops tauschii, the discovery of disease resistance genes from Aegilops tauschii, and the application of disease resistance genes from Aegilops tauschii to modern wheat breeding, providing a reference for the further exploration and utilization of Aegilops tauschii in wheat disease resistance breeding.

Navigating the CoGe Online Software Suite for Polyploidy Research

The CoGe software suite at genomevolution.org hosts a number of tools that facilitate genomic research on plant and animal whole-genome multiplication-polyploidy. SynMap permits analysis and visualization of two-way syntenic dotplot alignments of genomes, includes many options and data/graphics download possibilities, and even permits three-genome synteny maps and interactive views. FractBias is a tool that operates within SynMap that permits calculation and graphic display of genome fragments (such as chromosomes) of one species mapped to another, displaying both blockwise homology depths and the extent of syntenic gene (syntelog) loss following polyploidy events. SynMap macrosynteny results can segue into the microsynteny tool GEvo, which provides genome-browser-like views of homologous genome blocks. CoGe FeatView allows call-up of given gene features already stored in the CoGe resource, and CoGeBlast permits searches for additional features that can be analyzed or downloaded further. Links from these tools can be fed into SynFind, which can find syntenic blocks surrounding a feature across multiple specified genomes while also simultaneously providing overall genome-wide syntenic depth calculations that can be interpreted to reflect polyploidy levels. Here, we describe basic use of these tools on the CoGe software suite.