Preferential gene retention increases the robustness of cold regulation in Brassicaceae and other plants after polyploidization

New insights into the cold tolerance of upland switchgrass by integrating a haplotype-resolved genome and multi-omics analysis

BACKGROUND

Switchgrass (Panicum virgatum L.) is a bioenergy and forage crop. Upland switchgrass exhibits superior cold tolerance compared to the lowland ecotype, but the underlying molecular mechanisms remain unclear.

RESULTS

Here, we present a high-quality haplotype-resolved genome of the upland ecotype "Jingji31." We then conduct multi-omics analysis to explore the mechanism underlying its cold tolerance. By comparative transcriptome analysis of the upland and lowland ecotypes, we identify many genes with ecotype-specific differential expression, particularly members of the cold-responsive (COR) gene family, under cold stress. Notably, AFB1, ATL80, HOS10, and STRS2 gene families show opposite expression changes between the two ecotypes. Based on the haplotype-resolved genome of "Jingji31," we detect more cold-induced allele-specific expression genes in the upland ecotype than in the lowland ecotype, and these genes are significantly enriched in the COR gene family. By genome-wide association study, we detect an association signal related to the overwintering rate, which overlaps with a selective sweep region containing a cytochrome P450 gene highly expressed under cold stress. Heterologous overexpression of this gene in rice alleviates leaf chlorosis and wilting under cold stress. We also verify that expression of this gene is suppressed by a structural variation in the promoter region.

CONCLUSIONS

Based on the high-quality haplotype-resolved genome and multi-omics analysis of upland switchgrass, we characterize candidate genes responsible for cold tolerance. This study advances our understanding of plant cold tolerance, which provides crop breeding for improved cold tolerance.

Systematic identification of R2R3-MYB S6 subfamily genes in Brassicaceae and its role in anthocyanin biosynthesis in Brassica crops

The Brassicaceae family includes Arabidopsis thaliana, various vegetables and oil crops. The R2R3-MYB genes of the S6 subfamily are crucial for regulating anthocyanin biosynthesis, however, their systematic identification in Brassicaceae plants is still incomplete. Here, we systematically identified homologous genes of R2R3-MYB transcription factors from the S6 subfamily across 31 Brassicaceae species. A total of 92 homologous genes were identified, with species representation ranging from 0 to 10 genes per species. Phylogenetic analysis classified these homologous genes into six distinct groups. Notably, approximately 70% of the homologous genes were found within the G6 group, indicating a high degree of evolutionary conservation. Furthermore, a phylogenetic analysis was conducted on 35 homologous genes obtained from six species within the U's triangle Brassica plants. The findings provided evidence of significant conservation among orthologous genes across species and demonstrated strong collinearity on subgenomic chromosomes, with notable tandem duplications observed on chromosomes A7 and C6. Subsequently, we predicted the cis-acting elements of these 35 homologous genes, and analyzed their structures, conserved motifs, and characteristic conserved domains, confirming the significant similarities between orthologous genes. Additionally, we employed white and purple flower rapeseed specimens to conduct qRT-PCR validation of the key genes and transcriptional regulators associated with the anthocyanin synthesis pathway. The results revealed significant differential expression of BnaPAP2.A7.b in purple flowers, alongside the differential expression of BnaPAP2.C6.d. Ultimately, based on previous research and the findings of this study, we propose a transcriptional regulatory framework to govern anthocyanin accumulation in distinct tissues or organs of B. napus. Our findings offer a novel perspective on the functional diversification of R2R3-MYB transcription factors within the S6 subfamily homologous genes, while also shedding light on the regulatory network governing anthocyanin biosynthesis in Brassicaceae species.

Mechanisms underlining Kelp (Saccharina japonica) adaptation to relative high seawater temperature

Saccharina japonica has been cultivated in China for almost a century. From Dalian to Fujian, the lowest and the highest seawater temperatures in the period of cultivation increased by 14℃ and 8℃, respectively. Its adaptation to elevated seawater temperature is an example of securing the natural habitats of a species. To decipher the mechanisms underlining S. japonica adaptation to relative high seawater temperature, we assembled ~ 516.3 Mb female gametophyte genome and ~ 540.3 Mb of the male, respectively. The gametophytes isolated from southern China kelp cultivars acclimated to the relative high seawater temperature by transforming amino acids, glycosylating protein, maintaining osmotic pressure, intensifying the innate immune system, and exhausting energy and reduction power through the PEP-pyruvate-oxaloacetate node and the iodine cycle. They adapted to the relative high seawater temperature by transforming amino acids, changing sugar metabolism and intensifying innate immune system. The sex of S. japonica was determined by HMG-sex, and around this male gametophyte determiner the stress tolerant genes become linked to or associated with.

Convergent evolution in angiosperms adapted to cold climates

Convergent and parallel evolution occur more frequently than previously thought. Here, we focus on the evolutionary adaptations of angiosperms at sub-zero temperatures. We begin by introducing the history of research on convergent and parallel evolution, defining all independent similarities as convergent evolution. Our analysis reveals that frost zones (periodic or constant), which cover 49.1% of Earth's land surface, host 137 angiosperm families, with over 90% of their species thriving in these regions. In this context, we revisit the global biogeography and evolutionary trajectories of plant traits, such as herbaceous form and deciduous leaves, that are thought to be evasion strategies for frost adaptation. At the physiological and molecular levels, many angiosperms have independently evolved cold acclimation mechanisms through multiple pathways in addition to the well-characterized C-repeat binding factor/dehydration-responsive element binding protein 1 (CBF/DREB1) regulatory pathway. These convergent adaptations have occurred across various molecular levels, including amino acid substitutions and changes in gene duplication and expression within the same or similar functional pathways; however, identical amino acid changes are rare. Our results also highlight the prevalence of polyploidy in frost zones and the occurrence of paleopolyploidization events during global cooling. These patterns suggest repeated evolution in cold climates. Finally, we discuss plant domestication and predict climate zone shifts due to global warming and their effects on plant migration and in situ adaptation. Overall, the integration of ecological and molecular perspectives is essential for understanding and forecasting plant responses to climate change.

Neopolyploidy-induced changes in giant duckweed (Spirodela polyrhiza) alter herbivore preference and performance and plant population performance

PREMISE

Polyploidy is a widespread mutational process in angiosperms that may alter population performance of not only plants but also their interacting species. Yet, knowledge of whether polyploidy affects plant-herbivore dynamics is scarce. Here, we tested whether aphid herbivores exhibit preference for diploid or neopolyploid plants, whether polyploidy impacts plant and herbivore performance, and whether these interactions depend on the plant genetic background.

METHODS

Using independently synthesized neotetraploid strains paired with their diploid progenitors of greater duckweed (Spirodela polyrhiza), we evaluated the effect of neopolyploidy on duckweed's interaction with the water-lily aphid (Rhopalosiphum nymphaeae). Using paired-choice experiments, we evaluated feeding preference of the herbivore. We then evaluated the consequences of polyploidy on aphid and plant performance by measuring population growth over multiple generations.

RESULTS

Aphids preferred neopolyploids when plants were provided at equal abundances but not at equal surface areas, suggesting the role of plant population surface area in driving this preference. Additionally, neopolyploidy increased aphid population performance, but this result was dependent on the plant's genetic lineage. Lastly, the impact of herbivory on neopolyploid vs. diploid duckweed varied greatly with genetic lineage, where neopolyploids appeared to be variably tolerant compared to diploids, sometimes mirroring the effect on herbivore performance.

CONCLUSIONS

By experimentally testing the impacts of polyploidy on trophic species interactions, we showed that polyploidization can impact the preference and performance of herbivores on their plant hosts. These results have significant implications for the establishment and persistence of plants and herbivores in the face of plant polyploidy.

ldi-miR396-LdPMaT1 enhances reactive oxygen species scavenging capacity and promotes drought tolerance in Lilium distichum Nakai autotetraploids

Drought is a key environmental stress that inhibits plant growth, development, yield and quality. Whole-genome replication is an effective method for breeding drought resistant cultivars. Here, we evaluated the tolerance of Lilium distichum Nakai diploids (2n = 2× = 24) and artificially induced autotetraploids (2n = 4× = 48) to drought simulated by polyethylene glycol (PEG) stress. Autotetraploids showed stronger drought tolerance than diploids, and high-throughput sequencing during PEG stress identified five differentially expressed miRNAs. Transcriptome analysis revealed significantly different reactive oxygen species (ROS)-scavenger expression levels between diploids and autotetraploids, which increased the drought tolerance of autotetraploids. Specifically, we identified ldi-miR396b and its only target gene (LdPMaT1) for further study based on its expression level and ROS-scavenging ability in response to drought stress (DS). Autotetraploids showed higher expression of LdPMaT1 and significantly downregulated expression of ldi-miR396b under DS compared with diploids. Through a short tandem target mimic (STTM) in transgenic lilies, functional studies revealed that miR396b silencing promotes LdPMaT1 expression and the DS response. Under PEG stress, STTM393 transgenic lines showed improved drought resistance mediated by lowered MDA content but exhibited high antioxidant enzyme activity, consistent with the autotetraploid results. Collectively, these findings suggest that ldi-miR396b-LdPMaT1 potentially enhances ROS-scavenging ability, which contributes to improved stress adaptation in autotetraploid lilies.

De novo assembling a high-quality genome sequence of Amur grape (Vitis amurensis Rupr.) gives insight into Vitis divergence and sex determination

To date, there has been no high-quality sequence for genomes of the East Asian grape species, hindering biological and breeding efforts to improve grape cultivars. This study presents ~522 Mb of the Vitis amurensis (Va) genome sequence containing 27 635 coding genes. Phylogenetic analysis indicated that Vitis riparia (Vr) may have first split from the other two species, Va and Vitis vinifera (Vv). Divergent numbers of duplicated genes reserved among grapes suggests that the core eudicot-common hexaploidy (ECH) and the subsequent genome instability still play a non-negligible role in species divergence and biological innovation. Prominent accumulation of sequence variants might have improved cold resistance in Va, resulting in a more robust network of regulatory cold resistance genes, explaining why it is extremely cold-tolerant compared with Vv and Vr. In contrast, Va has preserved many fewer nucleotide binding site (NBS) disease resistance genes than the other grapes. Notably, multi-omics analysis identified one trans-cinnamate 4-monooxygenase gene positively correlated to the resveratrol accumulated during Va berry development. A selective sweep analysis revealed a hypothetical Va sex-determination region (SDR). Besides, a PPR-containing protein-coding gene in the hypothetical SDR may be related to sex determination in Va. The content and arrangement order of genes in the putative SDR of female Va were similar to those of female Vv. However, the putative SDR of female Va has lost one flavin-containing monooxygenase (FMO) gene and contains one extra protein-coding gene uncharacterized so far. These findings will improve the understanding of Vitis biology and contribute to the improvement of grape breeding.

Integrated anatomical structure, physiological, and transcriptomic analyses to identify differential cold tolerance responses of Ziziphus jujuba mill. 'Yueguang' and its autotetraploid 'Hongguang'

Cold stress is a limiting stress factor that limits plant distribution and development; however, polyploid plants have specific characteristics such as higher resistance to abiotic stress, especially cold stress, that allow them to overcome this challenge. The cultivated cultivar Ziziphus jujuba Mill. 'Yueguang' (YG) and its autotetraploid counterpart 'Hongguang' (HG) exhibit differential cold tolerance. However, the underlying molecular mechanism and methods to enhance their cold tolerance remain unknown. Anatomical structure and physiological analysis indicated YG had a higher wood bark ratio, and xylem ratio under cold treatment compared to HG. However, the half-lethal temperature (LT50), cortex ratio, and malondialdehyde (MDA) content were significantly decreased in YG than HG, which indicated YG was cold tolerant than HG. Transcriptome analysis showed that 2084, 1725, 2888, and 2934 differentially expressed genes (DEGs) were identified in HC vs YC, H20 vs Y20, Y20 vs YC, and H20 vs HC treatment, respectively. Meanwhile, KEGG enrichment analysis of DEGs showed that several metabolic pathways, primarily plant hormone signal transduction and the MAPK signaling pathway, were involved in the differential regulation of cold tolerance between YG and HG. Furthermore, exogenous abscisic acid (ABA) and brassinolide (BR) treatments could improve their cold tolerance through increased SOD and POD activities, decreased relative electrical conductivity, and MDA content. All of these findings suggested that plant hormone signal transduction, particularly ABA and BR, might have an important role in the regulation of differential cold tolerance between YG and HG, laying the foundation for further improving cold tolerance in jujube and examining the molecular mechanisms underlying differences in cold tolerance among different ploidy cultivars.

Convergent and/or parallel evolution of RNA-binding proteins in angiosperms after polyploidization

Increasing studies suggest that the biased retention of stress-related transcription factors (TFs) after whole-genome duplications (WGDs) could rewire gene transcriptional networks, facilitating plant adaptation to challenging environments. However, the role of posttranscriptional factors (e.g. RNA-binding proteins, RBPs) following WGDs has been largely ignored. Uncovering thousands of RBPs in 21 representative angiosperm species, we integrate genomic, transcriptomic, regulatomic, and paleotemperature datasets to unravel their evolutionary trajectories and roles in adapting to challenging environments. We reveal functional enrichments of RBP genes in stress responses and identify their convergent retention across diverse angiosperms from independent WGDs, coinciding with global cooling periods. Numerous RBP duplicates derived from WGDs are then identified as cold-induced. A significant overlap of 29 orthogroups between WGD-derived and cold-induced RBP genes across diverse angiosperms highlights a correlation between WGD and cold stress. Notably, we unveil an orthogroup (Glycine-rich RNA-binding Proteins 7/8, GRP7/8) and relevant TF duplicates (CCA1/LHY, RVE4/8, CBF2/4, etc.), co-retained in different angiosperms post-WGDs. Finally, we illustrate their roles in rewiring circadian and cold-regulatory networks at both transcriptional and posttranscriptional levels during global cooling. Altogether, we underline the adaptive evolution of RBPs in angiosperms after WGDs during global cooling, improving our understanding of plants surviving periods of environmental turmoil.

Coping with alpine habitats: genomic insights into the adaptation strategies of Triplostegia glandulifera (Caprifoliaceae)

How plants find a way to thrive in alpine habitats remains largely unknown. Here we present a chromosome-level genome assembly for an alpine medicinal herb, Triplostegia glandulifera (Caprifoliaceae), and 13 transcriptomes from other species of Dipsacales. We detected a whole-genome duplication event in T. glandulifera that occurred prior to the diversification of Dipsacales. Preferential gene retention after whole-genome duplication was found to contribute to increasing cold-related genes in T. glandulifera. A series of genes putatively associated with alpine adaptation (e.g. CBFs, ERF-VIIs, and RAD51C) exhibited higher expression levels in T. glandulifera than in its low-elevation relative, Lonicera japonica. Comparative genomic analysis among five pairs of high- vs low-elevation species, including a comparison of T. glandulifera and L. japonica, indicated that the gene families related to disease resistance experienced a significantly convergent contraction in alpine plants compared with their lowland relatives. The reduction in gene repertory size was largely concentrated in clades of genes for pathogen recognition (e.g. CNLs, prRLPs, and XII RLKs), while the clades for signal transduction and development remained nearly unchanged. This finding reflects an energy-saving strategy for survival in hostile alpine areas, where there is a tradeoff with less challenge from pathogens and limited resources for growth. We also identified candidate genes for alpine adaptation (e.g. RAD1, DMC1, and MSH3) that were under convergent positive selection or that exhibited a convergent acceleration in evolutionary rate in the investigated alpine plants. Overall, our study provides novel insights into the high-elevation adaptation strategies of this and other alpine plants.

Chromosome-level genome of putative autohexaploid Actinidia deliciosa provides insights into polyploidisation and evolution

Actinidia ('Mihoutao' in Chinese) includes species with complex ploidy, among which diploid Actinidia chinensis and hexaploid Actinidia deliciosa are economically and nutritionally important fruit crops. Actinidia deliciosa has been proposed to be an autohexaploid (2n = 174) with diploid A. chinensis (2n = 58) as the putative parent. A CCS-based assembly anchored to a high-resolution linkage map provided a chromosome-resolved genome for hexaploid A. deliciosa yielded a 3.91-Gb assembly of 174 pseudochromosomes comprising 29 homologous groups with 6 members each, which contain 39 854 genes with an average of 4.57 alleles per gene. Here we provide evidence that much of the hexaploid genome matches diploid A. chinensis; 95.5% of homologous gene pairs exhibited >90% similarity. However, intragenome and intergenome comparisons of synteny indicate chromosomal changes. Our data, therefore, indicate that if A. deliciosa is an autoploid, chromosomal rearrangement occurred following autohexaploidy. A highly diversified pattern of gene expression and a history of rapid population expansion after polyploidisation likely facilitated the adaptation and niche differentiation of A. deliciosa in nature. The allele-defined hexaploid genome of A. deliciosa provides new genomic resources to accelerate crop improvement and to understand polyploid genome evolution.

Physio-Biochemical Insights into the Cold Resistance Variations among Nectarine (Prunus persica (L.) Batsch var. nectarina) Cultivars

Cold stress occurs in late winter and early spring threatens greatly the nectarine industry. In this study, the semi-lethal low temperature (LT50) and thirteen cold resistance related parameters of five nectarine cultivars, including 'Nonglehong little princess' (LP), 'Luyou No. 5' (LY), 'Nonglehong No. 6' (NL), 'Zhongyou No. 20' (ZY) and 'Qiuhongzhu' (QH), were determined. Based on these parameters, they were categorized into high-(HR, including NL and LP), moderate-(MR, including QH) and low-cold resistant (LR, including ZY and LY) groups. The relative water (RW), proline (PRO), soluble sucrose (SS) and soluble protein (SP) contents, and superoxide dismutase (SOD) and peroxidase (POD) activities of HR cultivars were higher while their relative electronic conductivity (RE), malondialdehyde (MDA) and gibberellin acid (GA3) contents and catalase (CAT) activity were lower than other cultivars during natural overwintering. Redundancy analysis revealed that the lowest temperature in a day (LT) and LT50 significantly explains 69.8% and 10.9% of these physiological variables, respectively. Moreover, GA3 and indoleacetic acid (IAA) contents and CAT activity were positively correlated, while PRO, SS, ABA and RW contents were negatively correlated with both LT and LT50. Our study will be helpful in understanding the cold resistance variations of nectarine germplasm resources.

The genome of Stephania japonica provides insights into the biosynthesis of cepharanthine

Stephania japonica is an early-diverging eudicotyledon plant with high levels of cepharanthine, proven to be effective in curing coronavirus infections. Here, we report a high-quality S. japonica genome. The genome size is 688.52 Mb, and 97.37% sequences anchor to 11 chromosomes. The genome comprises 67.46% repetitive sequences and 21,036 genes. It is closely related to two Ranunculaceae species, which diverged from their common ancestor 55.90-71.02 million years ago (Mya) with a whole-genome duplication 85.59-96.75 Mya. We further reconstruct ancestral karyotype of Ranunculales. Several cepharanthine biosynthesis genes are identified and verified by western blot. Two genes (Sja03G0243 and Sja03G0241) exhibit catalytic activity as shown by liquid chromatography-mass spectrometry. Then, cepharanthine biosynthesis genes, transcription factors, and CYP450 family genes are used to construct a comprehensive network. Finally, we construct an early-diverging eudicotyledonous genome resources (EEGR) database. As the first genome of the Menispermaceae family to be released, this study provides rich resources for genomic studies.

Adaptive evolution of chloroplast division mechanisms during plant terrestrialization

Despite extensive research, the origin and evolution of the chloroplast division machinery remain unclear. Here, we employ recently sequenced genomes and transcriptomes of Archaeplastida clades to identify the core components of chloroplast division and reconstruct their evolutionary histories, respectively. Our findings show that complete division ring structures emerged in Charophytes. We find that Glaucophytes experienced strong selection pressure, generating diverse variants adapted to the changing terrestrial environments. By integrating the functions of chloroplast division genes (CDGs) annotated in a workflow developed using large-scale multi-omics data, we further show that dispersed duplications acquire more species-specific functions under stronger selection pressures. Notably, PARC6, a dispersed duplicate CDG, regulates leaf color and plant growth in Solanum lycopersicum, demonstrating neofunctionalization. Our findings provide an integrated perspective on the functional evolution of chloroplast division machinery and highlight the potential of dispersed duplicate genes as the primary source of adaptive evolution of chloroplast division.

Functional Differentiation of the Duplicated Gene BrrCIPK9 in Turnip (Brassica rapa var. rapa)

Gene duplication is a key biological process in the evolutionary history of plants and an important driving force for the diversification of genomic and genetic systems. Interactions between the calcium sensor calcineurin B-like protein (CBL) and its target, CBL-interacting protein kinase (CIPK), play important roles in the plant's response to various environmental stresses. As a food crop with important economic and research value, turnip (Brassica rapa var. rapa) has been well adapted to the environment of the Tibetan Plateau and become a traditional crop in the region. The BrrCIPK9 gene in turnip has not been characterized. In this study, two duplicated genes, BrrCIPK9.1 and BrrCIPK9.2, were screened from the turnip genome. Based on the phylogenetic analysis, BrrCIPK9.1 and BrrCIPK9.2 were found located in different sub-branches on the phylogenetic tree. Real-time fluorescence quantitative PCR analyses revealed their differential expression levels between the leaves and roots and in response to various stress treatments. The differences in their interactions with BrrCBLs were also revealed by yeast two-hybrid analyses. The results indicate that BrrCIPK9.1 and BrrCIPK9.2 have undergone Asparagine-alanine-phenylalanine (NAF) site divergence during turnip evolution, which has resulted in functional differences between them. Furthermore, BrrCIPK9.1 responded to high-pH (pH 8.5) stress, while BrrCIPK9.2 retained its ancestral function (low K+), thus providing further evidence of their functional divergence. These functional divergence genes facilitate turnip's good adaptation to the extreme environment of the Tibetan Plateau. In summary, the results of this study reveal the characteristics of the duplicated BrrCIPK9 genes and provide a basis for further functional studies of BrrCBLs-BrrCIPKs in turnip.

PHGD: An integrative and user-friendly database for plant hormone-related genes

Plant Hormone Gene Database (PHGD) database platform construction pipeline. First, we collected all reported hormone-related genes in the model plant Arabidopsis thaliana, and combined with the existing experimental background, mapped the hormone-gene interaction network to provide a blueprint. Next, we collected 469 high-quality plant genomes. Then, bioinformatics was used to identify hormone-related genes in these plants. Finally, these genetic data were programmed to be stored in a database and a platform website PHGD was built. PHGD was divided into eight modules, namely Home, Browse, Search, Resources, Download, Tools, Help, and Contact. We provided data resources and platform services to facilitate the study of plant hormones.

Large-scale analysis of trihelix transcription factors reveals their expansion and evolutionary footprint in plants

The trihelix transcription factor (TTF) gene family is an important class of transcription factors that play key roles in regulating developmental processes and responding to various stresses. To date, no comprehensive analysis of the TTF gene family in large-scale species has been performed. A cross-genome exploration of its origin, copy number variation, and expression pattern in plants is also unavailable. Here, we identified and characterized the TTF gene family in 110 species representing typical plant phylogenetic taxa. Interestingly, we found that the number of TTF genes was significantly expanded in Chara braunii compared to other species. Based on the available plant genomic datasets, our comparative analysis suggested that the TTF gene family likely originated from the GT-1-1 group and then expanded to form other groups through duplication or deletion of some domains. We found evidence that whole-genome duplication/triplication contributed most to the expansion of the TTF gene family in dicots, monocots and basal angiosperms. In contrast, dispersed and proximal duplications contributed to the expansion of the TTF gene family in algae and bryophyta. The expression patterns of TTF genes and their upstream and downstream genes in different treatments showed a functional divergence of TTF-related genes. Furthermore, we constructed the interaction network between TTF genes and the corresponding upstream and downstream genes, providing a blueprint for their regulatory pathways. This study provided a cross-genome comparative analysis of TTF genes in 110 species, which contributed to understanding their copy number expansion and evolutionary footprint in plants.

Whole-genome and dispersed duplication, including transposed duplication, jointly advance the evolution of TLP genes in seven representative Poaceae lineages

BACKGROUND

In the evolutionary study of gene families, exploring the duplication mechanisms of gene families helps researchers understand their evolutionary history. The tubby-like protein (TLP) family is essential for growth and development in plants and animals. Much research has been done on its function; however, limited information is available with regard to the evolution of the TLP gene family. Herein, we systematically investigated the evolution of TLP genes in seven representative Poaceae lineages.

RESULTS

Our research showed that the evolution of TLP genes was influenced not only by whole-genome duplication (WGD) and dispersed duplication (DSD) but also by transposed duplication (TRD), which has been neglected in previous research. For TLP family size, we found an evolutionary pattern of progressive shrinking in the grass family. Furthermore, the evolution of the TLP gene family was at least affected by evolutionary driving forces such as duplication, purifying selection, and base mutations.

CONCLUSIONS

This study presents the first comprehensive evolutionary analysis of the TLP gene family in grasses. We demonstrated that the TLP gene family is also influenced by a transposed duplication mechanism. Several new insights into the evolution of the TLP gene family are presented. This work provides a good reference for studying gene evolution and the origin of duplication.

CHH methylation of genes associated with fatty acid and jasmonate biosynthesis contributes to cold tolerance in autotetraploids of Poncirus trifoliata

Polyploids have elevated stress tolerance, but the underlying mechanisms remain largely elusive. In this study, we showed that naturally occurring tetraploid plants of trifoliate orange (Poncirus trifoliata (L.) Raf.) exhibited enhanced cold tolerance relative to their diploid progenitors. Transcriptome analysis revealed that whole-genome duplication was associated with higher expression levels of a range of well-characterized cold stress-responsive genes. Global DNA methylation profiling demonstrated that the tetraploids underwent more extensive DNA demethylation in comparison with the diploids under cold stress. CHH methylation in the promoters was associated with up-regulation of related genes, whereas CG, CHG, and CHH methylation in the 3'-regions was relevant to gene down-regulation. Of note, genes involved in unsaturated fatty acids (UFAs) and jasmonate (JA) biosynthesis in the tetraploids displayed different CHH methylation in the gene flanking regions and were prominently up-regulated, consistent with greater accumulation of UFAs and JA when exposed to the cold stress. Collectively, our findings explored the difference in cold stress response between diploids and tetraploids at both transcriptional and epigenetic levels, and gained new insight into the molecular mechanisms underlying enhanced cold tolerance of the tetraploid. These results contribute to uncovering a novel regulatory role of DNA methylation in better cold tolerance of polyploids.

Convergent evolution of AP2/ERF III and IX subfamilies through recurrent polyploidization and tandem duplication during eudicot adaptation to paleoenvironmental changes

Whole-genome duplication (WGD or polyploidization) has been suggested as a genetic contributor to angiosperm adaptation to environmental changes. However, many eudicot lineages did not undergo recent WGD (R-WGD) around and/or after the Cretaceous-Paleogene (K-Pg) boundary, times of severe environmental changes; how those plants survived has been largely ignored. Here, we collected 22 plants from major branches of the eudicot phylogeny and classified them into two groups according to the occurrence or absence of R-WGD: 12 R-WGD-containing plants (R-WGD-Y) and 10 R-WGD-lacking plants (R-WGD-N). Subsequently, we identified 496 gene-rich families in R-WGD-Y and revealed that members of the AP2/ERF transcription factor family were convergently over-retained after R-WGDs and showed exceptional cold stimulation. The evolutionary trajectories of the AP2/ERF family were then compared between R-WGD-Y and R-WGD-N to reveal convergent expansions of the AP2/ERF III and IX subfamilies through recurrent independent WGDs and tandem duplications (TDs) after the radiation of the plants. The expansions showed coincident enrichments in- times around and/or after the K-Pg boundary, when global cooling was a major environmental stressor. Consequently, convergent expansions and co-retentions of AP2/ERF III C-repeat binding factor (CBF) duplicates and their regulons in different eudicot lineages contributed to the rewiring of cold-specific regulatory networks. Moreover, promoter analysis of cold-responsive AP2/ERF genes revealed an underlying cis-regulatory code (G-box: CACGTG). We propose a seesaw model of WGDs and TDs in the convergent expansion of AP2/ERF III and IX genes that has contributed to eudicot adaptation during paleoenvironmental changes, and we suggest that TD may be a reciprocal/alternative mechanism for genetic innovation in plants that lack WGD.

Gene Recruitments and Dismissals in the Argonaut Genome Provide Insights into Pelagic Lifestyle Adaptation and Shell-like Eggcase Reacquisition

The paper nautilus or greater argonaut, Argonauta argo, is a species of octopods which is characterized by its pelagic lifestyle and by the presence of a protective spiral-shaped shell-like eggcase in females. To reveal the genomic background of how the species adapted to the pelagic lifestyle and acquired its shell-like eggcase, we sequenced the draft genome of the species. The genome size was 1.1 Gb, which is the smallest among the cephalopods known to date, with the top 215 scaffolds (average length 5,064,479 bp) covering 81% (1.09 Gb) of the total assembly. A total of 26,433 protein-coding genes were predicted from 16,802 assembled scaffolds. From these, we identified nearly intact HOX, Parahox, Wnt clusters, and some gene clusters that could probably be related to the pelagic lifestyle, such as reflectin, tyrosinase, and opsin. The gene models also revealed several homologous genes related to calcified shell formation in Conchiferan mollusks, such as Pif-like, SOD, and TRX. Interestingly, comparative genomics analysis revealed that the homologous genes for such genes were also found in the genome of the shell-less octopus, as well as Nautilus, which has a true outer shell. Therefore, the draft genome sequence of Arg. argo presented here has helped us to gain further insights into the genetic background of the dynamic recruitment and dismissal of genes to form an important, converging extended phenotypic structure such as the shell and the shell-like eggcase. Additionally, it allows us to explore the evolution of from benthic to pelagic lifestyles in cephalopods and octopods.

Two independent allohexaploidizations and genomic fractionation in Solanales

Solanales, an order of flowering plants, contains the most economically important vegetables among all plant orders. To date, many Solanales genomes have been sequenced. However, the evolutionary processes of polyploidization events in Solanales and the impact of polyploidy on species diversity remain poorly understood. We compared two representative Solanales genomes (Solanum lycopersicum L. and Ipomoea triloba L.) and the Vitis vinifera L. genome and confirmed two independent polyploidization events. Solanaceae common hexaploidization (SCH) and Convolvulaceae common hexaploidization (CCH) occurred ∼43-49 and ∼40-46 million years ago (Mya), respectively. Moreover, we identified homologous genes related to polyploidization and speciation and constructed multiple genomic alignments with V. vinifera genome, providing a genomic homology framework for future Solanales research. Notably, the three polyploidization-produced subgenomes in both S. lycopersicum and I. triloba showed significant genomic fractionation bias, suggesting the allohexaploid nature of the SCH and CCH events. However, we found that the higher genomic fractionation bias of polyploidization-produced subgenomes in Solanaceae was likely responsible for their more abundant species diversity than that in Convolvulaceae. Furthermore, through genomic fractionation and chromosomal structural variation comparisons, we revealed the allohexaploid natures of SCH and CCH, both of which were formed by two-step duplications. In addition, we found that the second step of two paleohexaploidization events promoted the expansion and diversity of β-amylase (BMY) genes in Solanales. These current efforts provide a solid foundation for future genomic and functional exploration of Solanales.

FLOWERING LOCUS T indel variants confer vernalization-independent and photoperiod-insensitive flowering of yellow lupin (Lupinus luteus L.)

Ongoing climate change has considerably reduced the seasonal window for crop vernalization, concurrently expanding cultivation area into northern latitudes with long-day photoperiod. To address these changes, cool season legume breeders need to understand molecular control of vernalization and photoperiod. A key floral transition gene integrating signals from these pathways is the Flowering locus T (FT). Here, a recently domesticated grain legume, yellow lupin (Lupinus luteus L.), was explored for potential involvement of FT homologues in abolition of vernalization and photoperiod requirements. Two FTa (LlutFTa1a and LlutFTa1b) and FTc (LlutFTc1 and LlutFTc2) homologues were identified and sequenced for two contrasting parents of a reference recombinant inbred line (RIL) population, an early-flowering cultivar Wodjil and a late-flowering wild-type P28213. Large deletions were detected in the 5' promoter regions of three FT homologues. Quantitative trait loci were identified for flowering time and vernalization response in the RIL population and in a diverse panel of wild and domesticated accessions. A 2227 bp deletion found in the LlutFTc1 promoter was linked with early phenology and vernalization independence, whereas LlutFTa1a and LlutFTc2 indels with photoperiod responsiveness. Comparative mapping highlighted convergence of FTc1 indel evolution in two Old World lupin species, addressing both artificial selection during domestication and natural adaptation to short season environmental conditions. We concluded that rapid flowering in yellow lupin is associated with the de-repression of the LlutFTc1 homologue from the juvenile phase, putatively due to the elimination of all binding sites in the promoter region for the AGAMOUS-like 15 transcription factor.

Transcriptomics and Genomics Analysis Uncover the Differentially Expressed Chlorophyll and Carotenoid-Related Genes in Celery

Celery (Apium graveolens L.), a plant from Apiaceae, is one of the most important vegetables and is grown worldwide. Carotenoids can capture light energy and transfer it to chlorophyll, which plays a central role in photosynthesis. Here, by performing transcriptomics and genomics analysis, we identified and conducted a comprehensive analysis of chlorophyll and carotenoid-related genes in celery and six representative species. Significantly, different contents and gene expression patterns were found among three celery varieties. In total, 237 and 290 chlorophyll and carotenoid-related genes were identified in seven species. No notable gene expansion of chlorophyll biosynthesis was detected in examined species. However, the gene encoding ζ-carotene desaturase (ZDS) enzyme in carotenoid was expanded in celery. Comparative genomics and RNA-seq analyses revealed 16 and 5 key genes, respectively, regulating chlorophyll and carotenoid. An intriguing finding is that chlorophyll and carotenoid-related genes were coordinately regulated by transcriptional factors, which could be distinctively classified into positive- and negative-regulation groups. Six CONSTANS (CO)-like transcription factors co-regulated chlorophyll and carotenoid-related genes were identified in celery. In conclusion, this study provides new insights into the regulation of chlorophyll and carotenoid by transcription factors.

Diversity, phylogeny, and adaptation of bryophytes: insights from genomic and transcriptomic data

Bryophytes including mosses, liverworts, and hornworts are among the earliest land plants, and occupy a crucial phylogenetic position to aid in the understanding of plant terrestrialization. Despite their small size and simple structure, bryophytes are the second largest group of extant land plants. They live ubiquitously in various habitats and are highly diversified, with adaptive strategies to modern ecosystems on Earth. More and more genomes and transcriptomes have been assembled to address fundamental questions in plant biology. Here, we review recent advances in bryophytes associated with diversity, phylogeny, and ecological adaptation. Phylogenomic studies have provided increasing supports for the monophyly of bryophytes, with hornworts sister to the Setaphyta clade including liverworts and mosses. Further comparative genomic analyses revealed that multiple whole-genome duplications might have contributed to the species richness and morphological diversity in mosses. We highlight that the biological changes through gene gain or neofunctionalization that primarily evolved in bryophytes have facilitated the adaptation to early land environments; among the strategies to adapt to modern ecosystems in bryophytes, desiccation tolerance is the most remarkable. More genomic information for bryophytes would shed light on key mechanisms for the ecological success of these 'dwarfs' in the plant kingdom.

Pervasive genome duplications across the plant tree of life and their links to major evolutionary innovations and transitions

Whole-genome duplication (WGD) has occurred repeatedly during plant evolution and diversification, providing genetic layers for evolving new functions and phenotypes. Advances in long-read sequencing technologies have enabled sequencing and assembly of over 1000 plant genomes spanning nearly 800 species, in which a large set of ancient WGDs has been uncovered. Here, we review the recently reported WGDs that occurred in major plant lineages and key evolutionary positions, and highlight their contributions to morphological innovation and adaptive evolution. Current gaps and challenges in integrating enormous volumes of sequenced plant genomes, accurately inferring WGDs, and developing web-based analysis tools are emphasized. Looking to the future, ambitious genome sequencing projects and global efforts may substantially recapitulate the plant tree of life based on broader sampling of phylogenetic diversity, reveal much of the timetable of ancient WGDs, and address the biological significance of WGDs in plant adaptation and radiation.

The fate of drought-related genes after polyploidization in Arachis hypogaea cv. Tifrunner

Drought stress affects plant growth and development. Cultivated peanut (Arachis hypogaea) was formed by a cross between A. duranensis and A. ipaensis. The drought tolerance of A. duranensis and A. ipaensis is reportedly stronger than that of cultivated peanut. However, there has been little study of drought tolerance genes in Arachis. In this study, we compared drought tolerance genes between A. hypogaea cv. Tifrunner and its diploid donors. We have observed that polyploidization does not generate more drought tolerance genes in A. hypogaea cv. Tifrunner but promotes the loss of many ancient drought tolerance genes. Although putative drought tolerance genes occurred on gene duplication events in A. hypogaea cv. Tifrunner, most copies lacked drought tolerance. These findings suggest that the loss of drought tolerance genes in A. hypogaea cv. Tifrunner could possibly result in weaker drought tolerance. In addition, we have observed that the three Arachis species stochastically lost putative drought tolerance genes. The evolution of drought tolerance genes could possibly have correlated with environmental changes. Our results enhance the current understanding of drought tolerance and polyploidy evolution in Arachis species.

Supplementary Information

The online version contains supplementary material available at 10.1007/s12298-022-01198-0.

Chromosome-level pepino genome provides insights into genome evolution and anthocyanin biosynthesis in Solanaceae

Pepino (Solanum muricatum, 2n = 2x = 24), a member of the Solanaceae family, is an important globally grown fruit. Herein, we report high-quality, chromosome-level pepino genomes. The 91.67% genome sequence is anchored to 12 chromosomes, with a total length of 1.20 Gb and scaffold N50 of 87.03 Mb. More than half the genome comprises repetitive sequences. In addition to the shared ancient whole-genome triplication (WGT) event in eudicots, an additional new WGT event was present in the pepino. Our findings suggest that pepinos experienced chromosome rearrangements, fusions, and gene loss after a WGT event. The large number of gene removals indicated the instability of Solanaceae genomes, providing opportunities for species divergence and natural selection. The paucity of disease-resistance genes (NBS) in pepino and eggplant has been explained by extensive loss and limited generation of genes after WGT events in Solanaceae. The outbreak of NBS genes was not synchronized in Solanaceae species, which occurred before the Solanaceae WGT event in pepino, tomato, and tobacco, whereas it was almost synchronized with WGT events in the other four Solanaceae species. Transcriptome and comparative genomic analyses revealed several key genes involved in anthocyanin biosynthesis. Although an extra WGT event occurred in Solanaceae, CHS genes related to anthocyanin biosynthesis in grapes were still significantly expanded compared with those in Solanaceae species. Proximal and tandem duplications contributed to the expansion of CHS genes. In conclusion, the pepino genome and annotation facilitate further research into important gene functions and comparative genomic analysis in Solanaceae.

Large-scale analyses of heat shock transcription factors and database construction based on whole-genome genes in horticultural and representative plants

Heat shock transcription factor (Hsf) plays a critical role in regulating heat resistance. Here, 2950 Hsf family genes were identified from 111 horticultural and representative plants. More Hsf genes were detected in higher plants than lower plants. Based on all Hsf genes, we constructed a phylogenetic tree, which indicated that Hsf genes of each branch evolved independently after species differentiation. Furthermore, we uncovered the evolutionary trajectories of Hsf genes by motif analysis. There were only 6 motifs (M1 to M6) in lower plants, and then 4 novel motifs (M7-M10) appeared in higher plants. However, the motifs of some Hsf genes were lost in higher plant, indicating that Hsf genes have undergone sequence variation during the evolution. The number of Hsf gene loss was more than duplication after whole-genome duplication in higher plants. The heat response network was constructed using 24 Hsf genes, 2421 downstream, and 222 upstream genes of Arabidopsis. Further enrichment analysis revealed that Hsf genes and other transcription factors interacted with each other to response heat resistance. The global expression maps were illustrated for Hsf genes under various abiotic, biotic stresses, and several developmental stages in Arabidopsis. The syntenic and phylogenetic analyses were conducted using Hsf genes of Arabidopsis and Pan-genome of 18 Brassica rapa accessions. We also performed the expression pattern analysis of Hsf and six Hsp family genes using expression values from different tissues and heat treatments in B. rapa. The interaction network between Hsf and Hsp gene families was constructed in B. rapa, and several core genes were detected in the network. Finally, we constructed a Hsf database (http://hsfdb.bio2db.com) for researchers to retrieve Hsf gene family information. Therefore, our study will provide rich resources for the evolution and functional study of Hsf genes.

Impact of polyploidy on plant tolerance to abiotic and biotic stresses

Polyploidy, defined as the coexistence of three or more complete sets of chromosomes in an organism's cells, is considered as a pivotal moving force in the evolutionary history of vascular plants and has played a major role in the domestication of several crops. In the last decades, improved cultivars of economically important species have been developed artificially by inducing autopolyploidy with chemical agents. Studies on diverse species have shown that the anatomical and physiological changes generated by either natural or artificial polyploidization can increase tolerance to abiotic and biotic stresses as well as disease resistance, which may positively impact on plant growth and net production. The aim of this work is to review the current literature regarding the link between plant ploidy level and tolerance to abiotic and biotic stressors, with an emphasis on the physiological and molecular mechanisms responsible for these effects, as well as their impact on the growth and development of both natural and artificially generated polyploids, during exposure to adverse environmental conditions. We focused on the analysis of those types of stressors in which more progress has been made in the knowledge of the putative morpho-physiological and/or molecular mechanisms involved, revealing both the factors in common, as well as those that need to be addressed in future research.

Plant protein-coding gene families: Their origin and evolution

Steady advances in genome sequencing methods have provided valuable insights into the evolutionary processes of several gene families in plants. At the core of plant biodiversity is an extensive genetic diversity with functional divergence and expansion of genes across gene families, representing unique phenomena. The evolution of gene families underpins the evolutionary history and development of plants and is the subject of this review. We discuss the implications of the molecular evolution of gene families in plants, as well as the potential contributions, challenges, and strategies associated with investigating phenotypic alterations to explain the origin of plants and their tolerance to environmental stresses.

Plant Breeding for Intercropping in Temperate Field Crop Systems: A Review

Monoculture cropping systems currently dominate temperate agroecosystems. However, intercropping can provide valuable benefits, including greater yield stability, increased total productivity, and resilience in the face of pest and disease outbreaks. Plant breeding efforts in temperate field crops are largely focused on monoculture production, but as intercropping becomes more widespread, there is a need for cultivars adapted to these cropping systems. Cultivar development for intercropping systems requires a systems approach, from the decision to breed for intercropping systems through the final stages of variety testing and release. Design of a breeding scheme should include information about species variation for performance in intercropping, presence of genotype × management interaction, observation of key traits conferring success in intercropping systems, and the specificity of intercropping performance. Together this information can help to identify an optimal selection scheme. Agronomic and ecological knowledge are critical in the design of selection schemes in cropping systems with greater complexity, and interaction with other researchers and key stakeholders inform breeding decisions throughout the process. This review explores the above considerations through three case studies: (1) forage mixtures, (2) perennial groundcover systems (PGC), and (3) soybean-pennycress intercropping. We provide an overview of each cropping system, identify relevant considerations for plant breeding efforts, describe previous breeding focused on the cropping system, examine the extent to which proposed theoretical approaches have been implemented in breeding programs, and identify areas for future development.

The reference genome and full-length transcriptome of pakchoi provide insights into cuticle formation and heat adaption

Brassica rapa includes various vegetables with high economic value. Among them, green petiole type pakchoi (B. rapa ssp. chinensis) is one of the major vegetables grown in southern China. Compared with other B. rapa varieties, green petiole type pakchoi shows a higher level of heat resistance, which is partially derived from the rich epicuticular wax. Here we sequence a high-quality genome of green petiole type pakchoi, which has been widely used as the parent in breeding. Our results reveal that long terminal repeat retrotransposon insertion plays critical roles in promoting the genome expansion and transcriptional diversity of pakchoi genes through preferential insertions, particularly in cuticle biosynthetic genes. After whole-genome triplication, over-retained pakchoi genes escape stringent selection pressure, and among them a set of cuticle-related genes are retained. Using bulked-segregant analysis of a heat-resistant pakchoi cultivar, we identify a frame-shift deletion across the third exon and the subsequent intron of BrcCER1 in candidate regions. Using Nanopore long-read sequencing, we analyze the full-length transcriptome of two pakchoi cultivars with opposite sensitivity to high temperature. We find that the heat-resistant pakchoi cultivar can mitigate heat-caused leaf damage by activating an unfolded protein response, as well as by inhibiting chloroplast development and energy metabolism, which are presumably mediated by both transcriptional regulation and splicing factors. Our study provides valuable resources for Brassica functional genomics and breeding research, and deepens our understanding of plant stress resistance.

Pan-genome of Raphanus highlights genetic variation and introgression among domesticated, wild, and weedy radishes

Post-polyploid diploidization associated with descending dysploidy and interspecific introgression drives plant genome evolution by unclear mechanisms. Raphanus is an economically and ecologically important Brassiceae genus and model system for studying post-polyploidization genome evolution and introgression. Here, we report the de novo sequence assemblies for 11 genomes covering most of the typical sub-species and varieties of domesticated, wild and weedy radishes from East Asia, South Asia, Europe, and America. Divergence among the species, sub-species, and South/East Asian types coincided with Quaternary glaciations. A genus-level pan-genome was constructed with family-based, locus-based, and graph-based methods, and whole-genome comparisons revealed genetic variations ranging from single-nucleotide polymorphisms (SNPs) to inversions and translocations of whole ancestral karyotype (AK) blocks. Extensive gene flow occurred between wild, weedy, and domesticated radishes. High frequencies of genome reshuffling, biased retention, and large-fragment translocation have shaped the genomic diversity. Most variety-specific gene-rich blocks showed large structural variations. Extensive translocation and tandem duplication of dispensable genes were revealed in two large rearrangement-rich islands. Disease resistance genes mostly resided on specific and dispensable loci. Variations causing the loss of function of enzymes modulating gibberellin deactivation were identified and could play an important role in phenotype divergence and adaptive evolution. This study provides new insights into the genomic evolution underlying post-polyploid diploidization and lays the foundation for genetic improvement of radish crops, biological control of weeds, and protection of wild species' germplasms.

Chloranthus genome provides insights into the early diversification of angiosperms

Chloranthales remain the last major mesangiosperm lineage without a nuclear genome assembly. We therefore assemble a high-quality chromosome-level genome of Chloranthus spicatus to resolve enigmatic evolutionary relationships, as well as explore patterns of genome evolution among the major lineages of mesangiosperms (eudicots, monocots, magnoliids, Chloranthales, and Ceratophyllales). We find that synteny is highly conserved between genomic regions of Amborella, Vitis, and Chloranthus. We identify an ancient single whole-genome duplication (WGD) (κ) prior to the divergence of extant Chloranthales. Phylogenetic inference shows Chloranthales as sister to magnoliids. Furthermore, our analyses indicate that ancient hybridization may account for the incongruent phylogenetic placement of Chloranthales + magnoliids relative to monocots and eudicots in nuclear and chloroplast trees. Long genes and long introns are found to be prevalent in both Chloranthales and magnoliids compared to other angiosperms. Overall, our findings provide an improved context for understanding mesangiosperm relationships and evolution and contribute a valuable genomic resource for future investigations.

The chromosome-scale genome of Magnolia officinalis provides insight into the evolutionary position of magnoliids

Magnolia officinalis, a representative tall aromatic tree of the Magnoliaceae family, is a medicinal plant that is widely used in diverse industries from medicine to cosmetics. We report a chromosome-scale draft genome of M. officinalis, in which ∼99.66% of the sequences were anchored onto 19 chromosomes with the scaffold N50 of 76.62 Mb. We found that a high proportion of repetitive sequences was a common feature of three Magnoliaceae with known genomic data. Magnoliids were a sister clade to eudicots-monocots, which provided more support for understanding the phylogenetic position among angiosperms. An ancient duplication event occurred in the genome of M. officinalis and was shared with Lauraceae. Based on RNA-seq analysis, we identified several key enzyme-coding gene families associated with the biosynthesis of lignans in the genome. The construction of the M. officinalis genome sequence will serve as a reference for further studies of Magnolia, as well as other Magnoliaceae.

Genome-Wide Comparative Analysis of Flowering-Time Genes; Insights on the Gene Family Expansion and Evolutionary Perspective

In polyploids, whole genome duplication (WGD) played a significant role in genome expansion, evolution and diversification. Many gene families are expanded following polyploidization, with the duplicated genes functionally diversified by neofunctionalization or subfunctionalization. These mechanisms may support adaptation and have likely contributed plant survival during evolution. Flowering time is an important trait in plants, which affects critical features, such as crop yields. The flowering-time gene family is one of the largest expanded gene families in plants, with its members playing various roles in plant development. Here, we performed genome-wide identification and comparative analysis of flowering-time genes in three palnt families i.e., Malvaceae, Brassicaceae, and Solanaceae, which indicate these genes were expanded following the event/s of polyploidization. Duplicated genes have been retained during evolution, although genome reorganization occurred in their flanking regions. Further investigation of sequence conservation and similarity network analyses provide evidence for functional diversification of duplicated genes during evolution. These functionally diversified genes play important roles in plant development and provide advantages to plants for adaptation and survival in response to environmental changes encountered during evolution. Collectively, we show that flowering-time genes were expanded following polyploidization and retained as large gene family by providing advantages from functional diversification during evolution.

Cytological and morphology characteristics of natural microsporogenesis within Camellia oleifera

Camellia oleifera is believed to exhibit a complex intraspecific polyploidy phenomenon. Abnormal microsporogenesis can promote the formation of unreduced gametes in plants and lead to sexual polyploidy, so it is hypothesized that improper meiosis probably results in the formation of natural polyploidy in Camellia oleifera. In this study, based on the cytological observation of meiosis in pollen mother cells (PMCs), we found natural 2n pollen for the first time in Camellia oleifera, which may lead to the formation of natural polyploids by sexual polyploidization. Additionally, abnormal cytological behaviour during meiosis, including univalent chromosomes, extraequatorial chromosomes, early segregation, laggard chromosomes, chromosome stickiness, asynchronous meiosis and deviant cytokinesis (monad, dyads, triads), was observed, which could be the cause of 2n pollen formation. Moreover, we confirmed a relationship among the length-width ratio of flower buds, stylet length and microsporogenesis. This result suggested that we can immediately determine the microsporogenesis stages by phenotypic characteristics, which may be applicable to breeding advanced germplasm in Camellia oleifera.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s12298-021-01002-5.

Comprehensive identification and characterization of simple sequence repeats based on the whole-genome sequences of 14 forest and fruit trees

Simple sequence repeats (SSRs) are popular and important molecular markers that exist widely in plants. Here, we conducted a comprehensive identification and comparative analysis of SSRs in 14 tree species. A total of 16, 298 SSRs were identified from 429, 449 genes, and primers were successfully designed for 99.44% of the identified SSRs. Our analysis indicated that tri-nucleotide SSRs were the most abundant, with an average of ~834 per species. Functional enrichment analysis by combining SSR-containing genes in all species, revealed 50 significantly enriched terms, with most belonging to transcription factor families associated with plant development and abiotic stresses such as Myeloblastosis_DNA-bind_4 (Myb_DNA-bind_4), APETALA2 (AP2), and Fantastic Four meristem regulator (FAF). Further functional enrichment analysis showed that 48 terms related to abiotic stress regulation and floral development were significantly enriched in ten species, whereas no significantly enriched terms were found in four species. Interestingly, the largest number of enriched terms was detected in Citrus sinensis (L.) Osbeck, accounting for 54.17% of all significantly enriched functional terms. Finally, we analyzed AP2 and trihelix gene families (Myb_DNA-bind_4) due to their significant enrichment in SSR-containing genes. The results indicated that whole-genome duplication (WGD) and whole genome triplication (WGT) might have played major roles in the expansion of the AP2 gene family but only slightly affected the expansion of the trihelix gene family during evolution. In conclusion, the identification and comprehensive characterization of SSR markers will greatly facilitate future comparative genomics and functional genomics studies.

Sucrose synthase gene family in Brassica juncea: genomic organization, evolutionary comparisons, and expression regulation

Sucrose synthase (SUS) plays an important role in sucrose metabolism and plant development. The SUS gene family has been identified in many plants, however, there is no definitive study of SUS gene in Brassica juncea. In this study, 14 SUS family genes were identified and comprehensively analyzed using bioinformatics tools. The analyzed parameters included their family member characteristics, chromosomal locations, gene structures and phylogenetic as well as transcript expression profiles. Phylogenetic analysis revealed that the 14 members could be allocated into three groups: SUS I, SUS II and SUS III. Comparisons of the exon/intron structure of the mustard SUS gene indicated that its structure is highly conserved. The conserved structure is attributed to purification selection during evolution. Expansion of the SUS gene family is associated with fragment and tandem duplications of the mustard SUS gene family. Collinearity analysis among species revealed that the SUS gene family could be lost or mutated to varying degrees after the genome was doubled, or when Brassica rapa and Brassica nigra hybridized to form Brassica juncea. The expression patterns of BjuSUSs vary among different stages of mustard stem swelling. Transcriptomics revealed that the BjuSUS01-04 expression levels were the most elevated. It has been hypothesized that they play an important role in sucrose metabolism during stem development. The expression levels of some BjuSUSs were significantly up-regulated when they were treated with plant hormones. However, when subjected to abiotic stress factors, their expression levels were suppressed. This study establishes SUS gene functions during mustard stem development and stress.

Molecular Evolution and Local Root Heterogeneous Expression of the Chenopodium quinoa ARF Genes Provide Insights into the Adaptive Domestication of Crops in Complex Environments

Auxin response factors (ARFs) influence plant growth and development via the coupling of basic biological processes. However, the evolution, expansion, and regulatory mechanisms of ARFs in the domesticated crop quinoa after artificial selection remain elusive. In this study, we systematically identified 30 Chenopodium quinoa ARFs (CqARFs). In this typical domesticated crop, ARFs divided into three subfamilies are subjected to strong purification selection and have a highly conserved evolutionary pattern. Polyploidy is the primary reason for the expansion of the ARF family after quinoa domestication. The expression patterns of CqARFs in different tissues have been differentiated, and CqARF2, 5, 9 and 10 from class A have the characteristics of local heterogeneous expression in different regions of roots, which may be the key factors for crops to respond in complex environments. Overall, we examined the evolution and expansion of ARFs in representative domesticated crops using the genome, transcriptome, and molecular biology and discovered a class A ARF-centered heterogeneous expression network that played an important role in auxin signaling and environmental responses. We provide new insights into how ARFs promote domesticated crop adaptation to artificial selection by polyploid expansion.

The ancient wave of polyploidization events in flowering plants and their facilitated adaptation to environmental stress

Flowering plants, or angiosperms, consist of more than 300,000 species, far more than any other land plant lineages. The accumulated evidence indicates that multiple ancient polyploidy events occurred around 100 to 120 million years ago during the Cretaceous and drove the early diversification of four major clades of angiosperms: gamma whole-genome triplication in the common ancestor of core eudicots, tau whole-genome duplication during the early diversification of monocots, lambda whole-genome duplication during the early diversification of magnoliids, and pi whole-genome duplication in the Nymphaeales lineage. These four polyploidy events have played essential roles in the adaptive evolution and diversification of major clades of flowering plants. Here, we specifically review the current understanding of this wave of ancient whole-genome duplications and their evolutionary significance. Notably, although these ancient whole-genome duplications occurred independently, they have contributed to the expansion of many stress-related genes (e.g., heat shock transcription factors and Arabidopsis response regulators)，and these genes could have been selected for by global environmental changes in the Cretaceous. Therefore, this ancient wave of paleopolyploidy events could have significantly contributed to the adaptation of angiosperms to environmental changes, and potentially promoted the wide diversification of flowering plants.

Identification, evolution and expression analyses of whole genome-wide TLP gene family in Brassica napus

BACKGROUND

Brassica is a very important genus of Brassicaceae, including many important oils, vegetables, forage crops, and ornamental horticultural plants. TLP family genes play important regulatory roles in the growth and development of plants. Therefore, this study used a bioinformatics approach to conduct the systematic comparative genomics analysis of TLP gene family in B. napus and other three important Brassicaceae crops.

RESULTS

Here, we identified a total of 29 TLP genes from B. napus genome, and they distributed on 16 chromosomes of B. napus. The evolutionary relationship showed that these genes could be divided into six groups from Group A to F. We found that the gene corresponding to Arabidopsis thaliana AT1G43640 was completely lost in B. rapa, B. oleracea and B. napus after whole genome triplication. The gene corresponding to AT1G25280 was retained in all the three species we analysed, belonging to 1:3:6 ratios. Our analyses suggested that there was a selective loss of some genes that might be redundant after genome duplication. This study proposed that the TLP genes in B. napus did not directly expansion compared with its diploid parents B. rapa, and B. oleracea. Instead, an indirect expansion of TLP gene family occurred in its two diploid parents. In addition, the study further utilized RNA-seq to detect the expression pattern of TLP genes between different tissues and two subgenomes.

CONCLUSIONS

This study systematically conducted the comparative analyses of TLP gene family in B. napus, discussed the loss and expansion of genes after genome duplication. It provided rich gene resources for exploring the molecular mechanism of TLP gene family. Meanwhile, it provided guidance and reference for the research of other gene families in B. napus.