Functional divergence of duplicated genes formed by polyploidy during Arabidopsis evolution

Sugar transporters PpSWEET9a and PpSWEET14 synergistically mediate peach sucrose allocation from source leaves to fruit

Sugar content is a critical indicator of fruit quality and is mainly controlled by sugar transporters. Sugars will eventually be exported transporters (SWEET) proteins play an indispensable role in sugar allocation between and within plant organs. Sucrose is the major sugar in many fruits and the predominant form of sugar translocated in peach (Prunus persica). However, the role of the multiple peach SWEET genes in sucrose allocation to fruit remains elusive. In this study, a total of 19 SWEET candidates have been identified in the peach genome, and two Clade III SWEET genes, PpSWEET9a and PpSWEET14, are found to be highly expressed in mature source leaves and branches. Complementation assays, transgene manipulations, and protein interaction studies reveal that PpSWEET9a and PpSWEET14 serve as sucrose efflux proteins and form a heterooligomer that synergistically directs sucrose allocation from source leaves to fruits. Our findings provide insights into the effect of SWEETs on sugar accumulation in peach fruit and identify genetic candidates for improving fruit quality.

Phosphorus-Use-Efficiency Gene Identification in Fabaceae and RSL2 Expansion in Lupinus albus Is Associated with Low-Phosphorus Adaptation

Phosphorus is critical for plant growth but often becomes less accessible due to its precipitation with cations in soil. Fabaceae, a diverse plant family, exhibits robust adaptability and includes species like Lupinus albus, known for its efficient phosphorus utilization via cluster roots. Here, we systematically identified phosphorus-utilization-efficiency (PUE) gene families across 35 Fabaceae species, highlighting significant gene amplification in PUE pathways in Fabaceae. Different PUE pathways exhibited variable amplification, evolution, and retention patterns among various Fabaceae crops. Additionally, the number of homologous genes of the root hair development gene RSL2 in L. albus was far more than that in other Fabaceae species. Multiple copies of the RSL2 gene were amplified and retained in L. albus after whole genome triplication. The gene structure and motifs specifically retained in L. albus were different from homologous genes in other plants. Combining transcriptome analysis under low-phosphorus treatment, it was found that most of the homologous genes of RSL2 in L. albus showed high expression in the cluster roots, suggesting that the RSL2 gene family plays an important role in the adaptation process of L. albus to low-phosphorus environments and the formation of cluster roots.

Duplication and sub-functionalization of flavonoid biosynthesis genes plays important role in Leguminosae root nodule symbiosis evolution

Gene innovation plays an essential role in trait evolution. Rhizobial symbioses, the most important N2-fixing agent in agricultural systems that exists mainly in Leguminosae, is one of the most attractive evolution events. However, the gene innovations underlying Leguminosae root nodule symbiosis (RNS) remain largely unknown. Here, we investigated the gene gain event in Leguminosae RNS evolution through comprehensive phylogenomic analyses. We revealed that Leguminosae-gain genes were acquired by gene duplication and underwent a strong purifying selection. Kyoto Encyclopedia of Genes and Genomes analyses showed that the innovated genes were enriched in flavonoid biosynthesis pathways, particular downstream of chalcone synthase (CHS). Among them, Leguminosae-gain type Ⅱ chalcone isomerase (CHI) could be further divided into CHI1A and CHI1B clades, which resulted from the products of tandem duplication. Furthermore, the duplicated CHI genes exhibited exon-intron structural divergences evolved through exon/intron gain/loss and insertion/deletion. Knocking down CHI1B significantly reduced nodulation in Glycine max (soybean) and Medicago truncatula; whereas, knocking down its duplication gene CHI1A had no effect on nodulation. Therefore, Leguminosae-gain type Ⅱ CHI participated in RNS and the duplicated CHI1A and CHI1B genes exhibited RNS functional divergence. This study provides functional insights into Leguminosae-gain genetic innovation and sub-functionalization after gene duplication that contribute to the evolution and adaptation of RNS in Leguminosae.

Biased Retention of Environment-Responsive Genes Following Genome Fractionation

The molecular underpinnings and consequences of cycles of whole-genome duplication (WGD) and subsequent gene loss through subgenome fractionation remain largely elusive. Endogenous drivers, such as transposable elements (TEs), have been postulated to shape genome-wide dominance and biased fractionation, leading to a conserved least-fractionated (LF) subgenome and a degenerated most-fractionated (MF) subgenome. In contrast, the role of exogenous factors, such as those induced by environmental stresses, has been overlooked. In this study, a chromosome-scale assembly of the alpine buckler mustard (Biscutella laevigata; Brassicaceae) that underwent a WGD event about 11 million years ago is coupled with transcriptional responses to heat, cold, drought, and herbivory to assess how gene expression is associated with differential gene retention across the MF and LF subgenomes. Counteracting the impact of TEs in reducing the expression and retention of nearby genes across the MF subgenome, dosage balance is highlighted as a main endogenous promoter of the retention of duplicated gene products under purifying selection. Consistent with the "turn a hobby into a job" model, about one-third of environment-responsive duplicates exhibit novel expression patterns, with one copy typically remaining conditionally expressed, whereas the other copy has evolved constitutive expression, highlighting exogenous factors as a major driver of gene retention. Showing uneven patterns of fractionation, with regions remaining unbiased, but with others showing high bias and significant enrichment in environment-responsive genes, this mesopolyploid genome presents evolutionary signatures consistent with an interplay of endogenous and exogenous factors having driven gene content following WGD-fractionation cycles.

Halophytes and heavy metals: A multi-omics approach to understand the role of gene and genome duplication in the abiotic stress tolerance of Cakile maritima

PREMISE

The origin of diversity is a fundamental biological question. Gene duplications are one mechanism that provides raw material for the emergence of novel traits, but evolutionary outcomes depend on which genes are retained and how they become functionalized. Yet, following different duplication types (polyploidy and tandem duplication), the events driving gene retention and functionalization remain poorly understood. Here we used Cakile maritima, a species that is tolerant to salt and heavy metals and shares an ancient whole-genome triplication with closely related salt-sensitive mustard crops (Brassica), as a model to explore the evolution of abiotic stress tolerance following polyploidy.

METHODS

Using a combination of ionomics, free amino acid profiling, and comparative genomics, we characterize aspects of salt stress response in C. maritima and identify retained duplicate genes that have likely enabled adaptation to salt and mild levels of cadmium.

RESULTS

Cakile maritima is tolerant to both cadmium and salt treatments through uptake of cadmium in the roots. Proline constitutes greater than 30% of the free amino acid pool in C. maritima and likely contributes to abiotic stress tolerance. We find duplicated gene families are enriched in metabolic and transport processes and identify key transport genes that may be involved in C. maritima abiotic stress tolerance.

CONCLUSIONS

These findings identify pathways and genes that could be used to enhance plant resilience and provide a putative understanding of the roles of duplication types and retention on the evolution of abiotic stress response.

The Characterization of a Novel PrMADS11 Transcription Factor from Pinus radiata Induced Early in Bent Pine Stem

A novel MADS-box transcription factor from Pinus radiata D. Don was characterized. PrMADS11 encodes a protein of 165 amino acids for a MADS-box transcription factor belonging to group II, related to the MIKC protein structure. PrMADS11 was differentially expressed in the stems of pine trees in response to 45° inclination at early times (1 h). Arabidopsis thaliana was stably transformed with a 35S::PrMADS11 construct in an effort to identify the putative targets of PrMADS11. A massive transcriptome analysis revealed 947 differentially expressed genes: 498 genes were up-regulated, and 449 genes were down-regulated due to the over-expression of PrMADS11. The gene ontology analysis highlighted a cell wall remodeling function among the differentially expressed genes, suggesting the active participation of cell wall modification required during the response to vertical stem loss. In addition, the phenylpropanoid pathway was also indicated as a PrMADS11 target, displaying a marked increment in the expression of the genes driven to the biosynthesis of monolignols. The EMSA assays confirmed that PrMADS11 interacts with CArG-box sequences. This TF modulates the gene expression of several molecular pathways, including other TFs, as well as the genes involved in cell wall remodeling. The increment in the lignin content and the genes involved in cell wall dynamics could be an indication of the key role of PrMADS11 in the response to trunk inclination.

Evolutionary divergence of subgenomes in common carp provides insights into speciation and allopolyploid success

Hybridization and polyploidization have made great contributions to speciation, heterosis, and agricultural production within plants, but there is still limited understanding and utilization in animals. Subgenome structure and expression reorganization and cooperation post hybridization and polyploidization are essential for speciation and allopolyploid success. However, the mechanisms have not yet been comprehensively assessed in animals. Here, we produced a high-fidelity reference genome sequence for common carp, a typical allotetraploid fish species cultured worldwide. This genome enabled in-depth analysis of the evolution of subgenome architecture and expression responses. Most genes were expressed with subgenome biases, with a trend of transition from the expression of subgenome A during the early stages to that of subgenome B during the late stages of embryonic development. While subgenome A evolved more rapidly, subgenome B contributed to a greater level of expression during development and under stressful conditions. Stable dominant patterns for homoeologous gene pairs both during development and under thermal stress suggest a potential fixed heterosis in the allotetraploid genome. Preferentially expressing either copy of a homoeologous gene at higher levels to confer development and response to stress indicates the dominant effect of heterosis. The plasticity of subgenomes and their shifting of dominant expression during early development, and in response to stressful conditions, provide novel insights into the molecular basis of the successful speciation, evolution, and heterosis of the allotetraploid common carp.

Genome-Wide Analysis of MYB Gene Family in Chrysanthemum ×morifolium Provides Insights into Flower Color Regulation

MYBs constitute the second largest transcription factor (TF) superfamily in flowering plants with substantial structural and functional diversity, which have been brought into focus because they affect flower colors by regulating anthocyanin biosynthesis. Up to now, the genomic data of several Chrysanthemum species have been released, which provides us with abundant genomic resources for revealing the evolution of the MYB gene family in Chrysanthemum species. In the present study, comparative analyses of the MYB gene family in six representative species, including C. lavandulifolium, C. seticuspe, C. ×morifolium, Helianthus annuus, Lactuca sativa, and Arabidopsis thaliana, were performed. A total of 1104 MYBs, which were classified into four subfamilies and 35 lineages, were identified in the three Chrysanthemum species (C. lavandulifolium, C. seticuspe, and C. ×morifolium). We found that whole-genome duplication and tandem duplication are the main duplication mechanisms that drove the occurrence of duplicates in CmMYBs (particularly in the R2R3-MYB subfamily) during the evolution of the cultivated chrysanthemums. Sequence structure and selective pressure analyses of the MYB gene family revealed that some of R2R3-MYBs were subjected to positive selection, which are mostly located on the distal telomere segments of the chromosomes and contain motifs 7 and 8. In addition, the gene expression analysis of CmMYBs in different organs and at various capitulum developmental stages of C. ×morifolium indicated that CmMYBS2, CmMYB96, and CmMYB109 might be the negative regulators for anthocyanin biosynthesis. Our results provide the phylogenetic context for research on the genetic and functional evolution of the MYB gene family in Chrysanthemum species and deepen our understanding of the regulatory mechanism of MYB TFs on the flower color of C. ×morifolium.

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.

OsGADD45a1: a multifaceted regulator of rice architecture, grain yield, and blast resistance

KEY MESSAGE

The homolog gene of the Growth Arrest and DNA Damage-inducible 45 (GADD45) in rice functions in the regulation of plant architecture, grain yield, and blast resistance. The Growth Arrest and DNA Damage-inducible 45 (GADD45) family proteins, well-established stress sensors and tumor suppressors in mammals, serve as pivotal regulators of genotoxic stress responses and tumorigenesis. In contrast, the homolog and role of GADD45 in plants have remained unclear. Herein, using forward genetics, we identified an activation tagging mutant AC13 exhibited dwarf characteristics resulting from the loss-of-function of the rice GADD45α homolog, denoted as OsGADD45a1. osgadd45a1 mutants displayed reduced plant height, shortened panicle length, and decreased grain yield compared to the wild-type Kitaake. Conversely, no obvious differences in plant height, panicle length, or grain yield were observed between wild-type and OsGADD45a1 overexpression plants. OsGADD45a1 displayed relatively high expression in germinated seeds and panicles, with localization in both the nucleus and cytoplasm. RNA-sequencing analysis suggested a potential role for OsGADD45a1 in the regulation of photosynthesis, and binding partner identification indicates OsGADD45a1 interacts with OsRML1 to regulate rice growth. Intriguingly, our study unveiled a novel role for OsGADD45a1 in rice blast resistance, as osgadd45a1 mutant showed enhanced resistance to Magnaporthe oryzae, and the expression of OsGADD45a1 was diminished upon blast fungus treatment. The involvement of OsGADD45a1 in rice blast fungus resistance presents a groundbreaking finding. In summary, our results shed light on the multifaceted role of OsGADD45a1 in rice, encompassing biotic stress response and the modulation of several agricultural traits, including plant height, panicle length, and grain yield.

Genome-wide identification of heat shock protein gene family and their responses to pathogen challenge in Trachinotus ovatus

Heat Shock Proteins (HSPs) are a widely distributed family of proteins produced in response to heat and other stresses. To develop a deeper understanding of the mechanisms governing expression of HSPs in the bony fish Trachinotus ovatus, we carried out a whole genome analysis and identified 43 HSP genes. Based on their phylogenetic relationships with Danio rerio, Seriola dumerili, and Seriola lalandi, they were divided into four subfamilies: HSP20, HSP60, HSP70, and HSP90. We performed an analysis of the predicted physicochemical properties and subcellular localization of proteins encoded by these genes. The chromosomal localization results showed that the HSP genes are distributed across 20 chromosomes of T. ovatus.These genes were found to be expressed in different tissues, and they showed differential expression in the immune response against Streptococcus agalactiae. However, there was no significant differential expression in the different skin tissue locations of T. ovatus after infection by Cryptocaryon irritans Brown. This study provides basic information for further research on the evolution and structure and function of HSPs in teleosts.

Characterization of the pyruvate kinase gene family in soybean and identification of a putative salt responsive gene GmPK21

BACKGROUND

As a key regulatory enzyme in the glycolysis pathway, pyruvate kinase (PK) plays crucial roles in multiple physiological processes during plant growth and is also involved in the abiotic stress response. However, little information is known about PKs in soybean.

RESULTS

In this study, we identified 27 PK family genes against the genome of soybean cultivar Zhonghuang13. They were classified into 2 subfamilies including PKc and PKp. 22 segmental duplicated gene pairs and 1 tandem duplicated gene pair were identified and all of them experienced a strong purifying selective pressure during evolution. Furthermore, the abiotic stresses (especially salt stress) and hormone responsive cis-elements were present in the promoters of GmPK genes, suggesting their potential roles in abiotic stress tolerance. By performing the qRT-PCR, 6 GmPK genes that continuously respond to both NaCl and ABA were identified. Subsequently, GmPK21, which represented the most significant change under NaCl treatment was chosen for further study. Its encoded protein GmPK21 was localized in the cytoplasm and plasma membrane. The transgenic Arabidopsis overexpressing GmPK21 exhibited weakened salinity tolerance.

CONCLUSIONS

This study provides genomic information of soybean PK genes and a molecular basis for mining salt tolerance function of PKs in the future.

GRASSY TILLERS1 (GT1) and SIX-ROWED SPIKE1 (VRS1) homologs share conserved roles in growth repression

Crop engineering and de novo domestication using gene editing are new frontiers in agriculture. However, outside of well-studied crops and model systems, prioritizing engineering targets remains challenging. Evolution can guide us, revealing genes with deeply conserved roles that have repeatedly been selected in the evolution of plant form. Homologs of the transcription factor genes GRASSY TILLERS1 (GT1) and SIX-ROWED SPIKE1 (VRS1) have repeatedly been targets of selection in domestication and evolution, where they repress growth in many developmental contexts. This suggests a conserved role for these genes in regulating growth repression. To test this, we determined the roles of GT1 and VRS1 homologs in maize (Zea mays) and the distantly related grass brachypodium (Brachypodium distachyon) using gene editing and mutant analysis. In maize, gt1; vrs1-like1 (vrl1) mutants have derepressed growth of floral organs. In addition, gt1; vrl1 mutants bore more ears and more branches, indicating broad roles in growth repression. In brachypodium, Bdgt1; Bdvrl1 mutants have more branches, spikelets, and flowers than wild-type plants, indicating conserved roles for GT1 and VRS1 homologs in growth suppression over ca. 59 My of grass evolution. Importantly, many of these traits influence crop productivity. Notably, maize GT1 can suppress growth in arabidopsis (Arabidopsis thaliana) floral organs, despite ca. 160 My of evolution separating the grasses and arabidopsis. Thus, GT1 and VRS1 maintain their potency as growth regulators across vast timescales and in distinct developmental contexts. This work highlights the power of evolution to inform gene editing in crop improvement.

Characteristics and expression of lncRNA and transposable elements in Drosophila aneuploidy

Aneuploidy can globally affect the expression of the whole genome, which is detrimental to organisms. Dosage-sensitive regulators usually have multiple intermolecular interactions, and changes in their stoichiometry are responsible for the dysregulation of the regulatory network. Currently, studies on noncoding genes in aneuploidy are relatively rare. We studied the characteristics and expression profiles of long noncoding RNAs (lncRNAs) and transposable elements (TEs) in aneuploid Drosophila. It is found that lncRNAs and TEs are affected by genomic imbalance and appear to be more sensitive to an inverse dosage effect than mRNAs. Several dosage-sensitive lncRNAs and TEs were detected for their expression patterns during embryogenesis, and their biological functions in the ovary and testes were investigated using tissue-specific RNAi. This study advances our understanding of the noncoding sequences in imbalanced genomes and provides a novel perspective for the study of aneuploidy-related human diseases such as cancer.

Divergent Retention of Sucrose Metabolism Genes after Whole Genome Triplication in the Tomato (Solanum lycopersicum)

Sucrose, the primary carbon transport mode and vital carbohydrate for higher plants, significantly impacts plant growth, development, yield, and quality formation. Its metabolism involves three key steps: synthesis, transport, and degradation. Two genome triplication events have occurred in Solanaceae, which have resulted in massive gene loss. In this study, a total of 48 and 65 genes from seven sucrose metabolism gene families in Vitis vinifera and Solanum lycopersicum were identified, respectively. The number of members comprising the different gene families varied widely. And there were significant variations in the pattern of gene duplication and loss in the tomato following two WGD events. Tandem duplication is a major factor in the expansion of the SWEET and Acid INV gene families. All the genes are irregularly distributed on the chromosomes, with the majority of the genes showing collinearity with the grape, particularly the CIN family. And the seven gene families were subjected to a purifying selection. The expression patterns of the different gene families exhibited notable variations. This study presents basic information about the sucrose metabolism genes in the tomato and grape, and paves the way for further investigations into the impact of SCT events on the phylogeny, gene retention duplication, and function of sucrose metabolism gene families in the tomato or Solanaceae, and the adaptive evolution of the tomato.

Genome resources for the elite bread wheat cultivar Aikang 58 and mining of elite homeologous haplotypes for accelerating wheat improvement

Despite recent progress in crop genomics studies, the genomic changes brought about by modern breeding selection are still poorly understood, thus hampering genomics-assisted breeding, especially in polyploid crops with compound genomes such as common wheat (Triticum aestivum). In this work, we constructed genome resources for the modern elite common wheat variety Aikang 58 (AK58). Comparative genomics between AK58 and the landrace cultivar Chinese Spring (CS) shed light on genomic changes that occurred through recent varietal improvement. We also explored subgenome diploidization and divergence in common wheat and developed a homoeologous locus-based genome-wide association study (HGWAS) approach, which was more effective than single homoeolog-based GWAS in unraveling agronomic trait-associated loci. A total of 123 major HGWAS loci were detected using a genetic population derived from AK58 and CS. Elite homoeologous haplotypes (HHs), formed by combinations of subgenomic homoeologs of the associated loci, were found in both parents and progeny, and many could substantially improve wheat yield and related traits. We built a website where users can download genome assembly sequence and annotation data for AK58, perform blast analysis, and run JBrowse. Our work enriches genome resources for wheat, provides new insights into genomic changes during modern wheat improvement, and suggests that efficient mining of elite HHs can make a substantial contribution to genomics-assisted breeding in common wheat and other polyploid crops.

Functional Diversification and the Plant Secondary Cell Wall

Much evidence exists suggesting the presence of genetic functional diversification in plants, though literature associated with the role of functional diversification in the evolution of the plant secondary cell wall (SCW) has sparsely been compiled and reviewed in a recent context. This review aims to elucidate, through the examination of gene phylogenies associated with its biosynthesis and maintenance, the role of functional diversification in shaping the critical, dynamic, and characteristic organelle, the secondary cell wall. It will be asserted that gene families resulting from gene duplication and subsequent functional divergence are present and are heavily involved in SCW biosynthesis and maintenance. Furthermore, diversification will be presented as a significant driver behind the evolution of the many functional characteristics of the SCW. The structure and function of the plant cell wall and its constituents will first be explored, followed by a discussion on the phenomenon of gene duplication and the resulting genetic functional divergence that can emerge. Finally, the major constituents of the SCW and their individual relationships with duplication and divergence will be reviewed to the extent of current knowledge on the subject.

The Soybean Expression Atlas v2: A comprehensive database of over 5000 RNA-seq samples

Soybean is a crucial crop worldwide, used as a source of food, feed, and industrial products due to its high protein and oil content. Previously, the rapid accumulation of soybean RNA-seq data in public databases and the computational challenges of processing raw RNA-seq data motivated us to develop the Soybean Expression Atlas, a gene expression database of over a thousand RNA-seq samples. Over the past few years, our database has allowed researchers to explore the expression profiles of important gene families, discover genes associated with agronomic traits, and understand the transcriptional dynamics of cellular processes. Here, we present the Soybean Expression Atlas v2, an updated version of our database with a fourfold increase in the number of samples, featuring transcript- and gene-level transcript abundance matrices for 5481 publicly available RNA-seq samples. New features in our database include the availability of transcript-level abundance estimates and equivalence classes to explore differential transcript usage, abundance estimates in bias-corrected counts to increase the accuracy of differential gene expression analyses, a new web interface with improved data visualization and user experience, and a reproducible and scalable pipeline available as an R package. The Soybean Expression Atlas v2 is available at https://soyatlas.venanciogroup.uenf.br/, and it will accelerate soybean research, empowering researchers with high-quality and easily accessible gene expression data.

Dynamics of cis-regulatory sequences and transcriptional divergence of duplicated genes in soybean

Transcriptional divergence of duplicated genes after whole genome duplication (WGD) has been described in many plant lineages and is often associated with subgenome dominance, a genome-wide mechanism. However, it is unknown what underlies the transcriptional divergence of duplicated genes in polyploid species that lack subgenome dominance. Soybean is a paleotetraploid with a WGD that occurred 5 to 13 Mya. Approximately 50% of the duplicated genes retained from this WGD exhibit transcriptional divergence. We developed accessible chromatin region (ACR) datasets from leaf, flower, and seed tissues using MNase-hypersensitivity sequencing. We validated enhancer function of several ACRs associated with known genes using CRISPR/Cas9-mediated genome editing. The ACR datasets were used to examine and correlate the transcriptional patterns of 17,111 pairs of duplicated genes in different tissues. We demonstrate that ACR dynamics are correlated with divergence of both expression level and tissue specificity of individual gene pairs. Gain or loss of flanking ACRs and mutation of cis-regulatory elements (CREs) within the ACRs can change the balance of the expression level and/or tissue specificity of the duplicated genes. Analysis of DNA sequences associated with ACRs revealed that the extensive sequence rearrangement after the WGD reshaped the CRE landscape, which appears to play a key role in the transcriptional divergence of duplicated genes in soybean. This may represent a general mechanism for transcriptional divergence of duplicated genes in polyploids that lack subgenome dominance.

Phylogenomic and syntenic data demonstrate complex evolutionary processes in early radiation of the rosids

Some of the most vexing problems of deep level relationship that remain in angiosperms involve the superrosids. The superrosid clade contains a quarter of all angiosperm species, with 18 orders in three subclades (Vitales, Saxifragales and core rosids) exhibiting remarkable morphological and ecological diversity. To help resolve deep-level relationships, we constructed a high-quality chromosome-level genome assembly for Tiarella polyphylla (Saxifragaceae) thus providing broader genomic representation of Saxifragales. Whole genome microsynteny analysis of superrosids showed that Saxifragales shared more synteny clusters with core rosids than Vitales, further supporting Saxifragales as more closely related with core rosids. To resolve the ordinal phylogeny of superrosids, we screened 122 single copy nuclear genes from genomes of 36 species, representing all 18 superrosid orders. Vitales were recovered as sister to all other superrosids (Saxifragales + core rosids). Our data suggest dramatic differences in relationships compared to earlier studies within core rosids. Fabids should be restricted to the nitrogen-fixing clade, while Picramniales, the Celastrales-Malpighiales (CM) clade, Huerteales, Oxalidales, Sapindales, Malvales and Brassicales formed an "expanded" malvid clade. The Celastrales-Oxalidales-Malpighiales (COM) clade (sensu APG IV) was not monophyletic. Crossosomatales, Geraniales, Myrtales and Zygophyllales did not belong to either of our well-supported malvids or fabids. There is strong discordance between nuclear and plastid phylogenetic hypotheses for superrosid relationships; we show that this is best explained by a combination of incomplete lineage sorting and ancient reticulation.

Genome-wide identification of DUF506 gene family in Oryzasativa and expression profiling under abiotic stresses

The domain of unknown function 560 (DUF560), also known as the PDDEXK_6 family, is a ubiquitous plant protein that has been confirmed to play critical roles in Arabidopsis root development as well as ABA and abiotic responses. However, genome-wide identification and expression pattern analysis in rice (Oryza sativa) still need to be improved. Based on the phylogenetic relationship, 10 OsDUF506 genes were identified and classified into four subfamilies. Segmental duplication was essential to the expansion of OsDUF506s, which were subjected to purifying selective pressure. Except for OsDUF50609 and OsDUF50610, the OsDUF506s shared colinear gene pairs with five monocot species, showing that they were conserved in evolution. Furthermore, the conserved domains, gene structures, SNPs distribution, and targeting miRNAs were systematically investigated. Massive cis-regulatory elements were discovered in promoter regions, implying that OsDUF506s may be important in hormone regulation and abiotic stress response. Therefore, we analyzed plant hormone-induced transcriptome data and performed qRT-PCR on eight OsDUF506s under drought, cold, and phosphorus-deficient stresses. The results revealed that most OsDUF506s respond to ABA and JA treatment, as well as drought and cold conditions. In conclusion, our findings provided insights into the evolution and function of OsDUF506s, which could benefit crop breeding in the future.

Dosage-sensitivity shapes how genes transcriptionally respond to allopolyploidy and homoeologous exchange in resynthesized Brassica napus

The gene balance hypothesis proposes that selection acts on the dosage (i.e. copy number) of genes within dosage-sensitive portions of networks, pathways, and protein complexes to maintain balanced stoichiometry of interacting proteins, because perturbations to stoichiometric balance can result in reduced fitness. This selection has been called dosage balance selection. Dosage balance selection is also hypothesized to constrain expression responses to dosage changes, making dosage-sensitive genes (those encoding members of interacting proteins) experience more similar expression changes. In allopolyploids, where whole-genome duplication involves hybridization of diverged lineages, organisms often experience homoeologous exchanges that recombine, duplicate, and delete homoeologous regions of the genome and alter the expression of homoeologous gene pairs. Although the gene balance hypothesis makes predictions about the expression response to homoeologous exchanges, they have not been empirically tested. We used genomic and transcriptomic data from 6 resynthesized, isogenic Brassica napus lines over 10 generations to identify homoeologous exchanges, analyzed expression responses, and tested for patterns of genomic imbalance. Groups of dosage-sensitive genes had less variable expression responses to homoeologous exchanges than dosage-insensitive genes, a sign that their relative dosage is constrained. This difference was absent for homoeologous pairs whose expression was biased toward the B. napus A subgenome. Finally, the expression response to homoeologous exchanges was more variable than the response to whole-genome duplication, suggesting homoeologous exchanges create genomic imbalance. These findings expand our knowledge of the impact of dosage balance selection on genome evolution and potentially connect patterns in polyploid genomes over time, from homoeolog expression bias to duplicate gene retention.

Analysis of gene duplication within the Arabidopsis NUCLEAR FACTOR Y, subunit B (NF-YB) protein family reveals domains under both purifying and diversifying selection

Gene duplication is an evolutionary mechanism that provides new genetic material. Since gene duplication is a major driver for molecular evolution, examining the fate of duplicated genes is an area of active research. The fate of duplicated genes can include loss, subfunctionalization, and neofunctionalization. In this manuscript, we chose to experimentally study the fate of duplicated genes using the Arabidopsis NUCLEAR FACTOR Y (NF-Y) transcription factor family. NF-Y transcription factors are heterotrimeric complexes, composed of NF-YA, NF-YB, and NF-YC. NF-YA subunits are responsible for nucleotide-specific binding to a CCAAT cis-regulatory element. NF-YB and NF-YC subunits make less specific, but essential complex-stabilizing contacts with the DNA flanking the core CCAAT pentamer. While ubiquitous in eukaryotes, each NF-Y family has expanded by duplication in the plant lineage. For example, the model plant Arabidopsis contains 10 each of the NF-Y subunits. Here we examine the fate of duplicated NF-YB proteins in Arabidopsis, which are composed of central histone fold domains (HFD) and less conserved flanking regions (N- and C-termini). Specifically, the principal question we wished to address in this manuscript was to what extent can the 10 Arabidopsis NF-YB paralogs functionally substitute the genes NF-YB2 and NF-YB3 in the promotion of photoperiodic flowering? Our results demonstrate that the conserved histone fold domains (HFD) may be under pressure for purifying (negative) selection, while the non-conserved N- and C-termini may be under pressure for diversifying (positive) selection, which explained each paralog's ability to substitute. In conclusion, our data demonstrate that the N- and C-termini may have allowed the duplicated genes to undergo functional diversification, allowing the retention of the duplicated genes.

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.

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.

Intra-Varietal Diversity and Its Contribution to Wheat Evolution, Domestication, and Improvement in Wheat

Crop genetic diversity is essential for adaptation and productivity in agriculture. A previous study revealed that poor allele diversity in wheat commercial cultivars is a major barrier to its further improvement. Homologs within a variety, including paralogs and orthologs in polyploid, account for a large part of the total genes of a species. Homolog diversity, intra-varietal diversity (IVD), and their functions have not been elucidated. Common wheat, an important food crop, is a hexaploid species with three subgenomes. This study analyzed the sequence, expression, and functional diversity of homologous genes in common wheat based on high-quality reference genomes of two representative varieties, a modern commercial variety Aikang 58 (AK58) and a landrace Chinese Spring (CS). A total of 85,908 homologous genes, accounting for 71.9% of all wheat genes, including inparalogs (IPs), outparalogs (OPs), and single-copy orthologs (SORs), were identified, suggesting that homologs are an important part of the wheat genome. The levels of sequence, expression, and functional variation in OPs and SORs were higher than that of IPs, which indicates that polyploids have more homologous diversity than diploids. Expansion genes, a specific type of OPs, made a great contribution to crop evolution and adaptation and endowed crop with special characteristics. Almost all agronomically important genes were from OPs and SORs, demonstrating their essential functions for polyploid evolution, domestication, and improvement. Our results suggest that IVD analysis is a novel approach for evaluating intra-genomic variations, and exploitation of IVD might be a new road for plant breeding, especially for polyploid crops, such as wheat.

Genome-Wide Identification of PEBP Gene Family in Solanum lycopersicum

The PEBP gene family is crucial for the growth and development of plants, the transition between vegetative and reproductive growth, the response to light, the production of florigen, and the reaction to several abiotic stressors. The PEBP gene family has been found in numerous species, but the SLPEBP gene family has not yet received a thorough bioinformatics investigation, and the members of this gene family are currently unknown. In this study, bioinformatics was used to identify 12 members of the SLPEBP gene family in tomato and localize them on the chromosomes. The physicochemical characteristics of the proteins encoded by members of the SLPEBP gene family were also examined, along with their intraspecific collinearity, gene structure, conserved motifs, and cis-acting elements. In parallel, a phylogenetic tree was built and the collinear relationships of the PEBP gene family among tomato, potato, pepper, and Arabidopsis were examined. The expression of 12 genes in different tissues and organs of tomato was analyzed using transcriptomic data. It was also hypothesized that SLPEBP3, SLPEBP5, SLPEBP6, SLPEBP8, SLPEBP9, and SLPEBP10 might be related to tomato flowering and that SLPEBP2, SLPEBP3, SLPEBP7, and SLPEBP11 might be related to ovary development based on the tissue-specific expression analysis of SLPEBP gene family members at five different stages during flower bud formation to fruit set. This article's goal is to offer suggestions and research directions for further study of tomato PEBP gene family members.

A Study of GUS Expression in Arabidopsis as a Tool for the Evaluation of Gene Evolution, Function and the Role of Expression Derived from Gene Duplication

Gene duplication played a fundamental role in eukaryote evolution and different copies of a given gene can be present in extant species, often with expressions and functions differentiated during evolution. We assume that, when such differentiation occurs in a gene copy, this may be indicated by its maintenance in all the derived species. To verify this hypothesis, we compared the histological expression domains of the three β-glucuronidase genes (AtGUS) present in Arabidopsis thaliana with the GUS evolutionary tree in angiosperms. We found that AtGUS gene expression overlaps in the shoot apex, the floral bud and the root hairs. In the root apex, AtGUS3 expression differs completely from AtGUS1 and AtGUS2, whose transcripts are present in the root cap meristem and columella, in the staminal cell niche, in the epidermis and in the proximal cortex. Conversely, AtGUS3 transcripts are limited to the old border-like cells of calyptra and those found along the protodermal cell line. The GUS evolutionary tree reveals that the two main clusters (named GUS1 and GUS3) originate from a duplication event predating angiosperm radiation. AtGUS3 belongs to the GUS3 cluster, while AtGUS1 and AtGUS2, which originate from a duplication event that occurred in an ancestor of the Brassicaceae family, are found together in the GUS1 cluster. There is another, previously undescribed cluster, called GUS4, originating from a very ancient duplication event. While the copy of GUS4 has been lost in many species, copies of GUS3 and GUS1 have been conserved in all species examined.

Predicting gene expression divergence between single-copy orthologs in two species

Predicting gene expression divergence is integral to understanding the emergence of new biological functions and associated traits. Whereas several sophisticated methods have been developed for this task, their applications are either limited to duplicate genes or require expression data from more than two species. Thus, here we present PiXi, the first machine learning framework for predicting gene expression divergence between single-copy orthologs in two species. PiXi models gene expression evolution as an Ornstein-Uhlenbeck process, and overlays this model with multi-layer neural network, random forest, and support vector machine architectures for making predictions. It outputs the predicted class "conserved" or "diverged" for each pair of orthologs, as well as their predicted expression optima in the two species. We show that PiXi has high power and accuracy in predicting gene expression divergence between single-copy orthologs, as well as high accuracy and precision in estimating their expression optima in the two species, across a wide range of evolutionary scenarios, with the globally best performance achieved by a multi-layer neural network. Moreover, application of our best performing PiXi predictor to empirical gene expression data from single-copy orthologs residing at different loci in two species of Drosophila reveals that approximately 23% underwent expression divergence after positional relocation. Further analysis shows that several of these "diverged" genes are involved in the electron transport chain of the mitochondrial membrane, suggesting that new chromatin environments may impact energy production in Drosophila. Thus, by providing a toolkit for predicting gene expression divergence between single-copy orthologs in two species, PiXi can shed light on the origins of novel phenotypes across diverse biological processes and study systems.

OAF is a DAF-like gene that controls ovule development in plants

We previously found that the RING-type E3 ligase DEFECTIVE IN ANTHER DEHISCENCE1- (DAD1-) Activating Factor (DAF) controls anther dehiscence by activating the jasmonate biosynthetic pathway in Arabidopsis. Here, we show that in Arabidopsis, the DAF ancestor was duplicated into three genes (DAF, Ovule Activating Factor (OAF), DAFL2), which evolved divergent partial functions from their ancestor through subfunctionalization. In this case, DAF-DAD1-JA signaling regulates anther dehiscence, whereas OAF controls ovule development by negatively regulating cinnamyl alcohol dehydrogenase 9 (CAD9) activity and being negatively regulated by miR847 itself in Arabidopsis. Downregulation of OAF or upregulation of CAD9 and miR847 caused similar abortion of ovule formation due to precocious ovule lignification in transgenic Arabidopsis. Interestingly, only one DAF-like gene, PaOAF, exists in the monocot orchids, which has likely evolved through nonfunctionalization and maintains a conserved function as Arabidopsis OAF in regulating ovule development since defective ovules were observed in the virus-induced gene silencing (VIGS) PaOAF Phalaenopsis orchids. The absence of the DAF ortholog and its function in orchids is likely due to the evolution of stamens to a unique pollinium structure that lacks the feature of anther dehiscence. These findings expand the current knowledge underlying the multifunctional evolution and diverse functionalization of duplicate gene pairs within/among plants.

Expansion and adaptive evolution of the WRKY transcription factor family in Avicennia mangrove trees

Mangroves are adapted to intertidal zones, which present extreme environmental conditions. WRKYs are among the most prominent transcription factors (TFs) in higher plants and act through various interconnected networks to regulate responses to multiple abiotic stressors. Here, based on omic data, we investigated the landscape and evolutionary patterns of WRKYs in the main mangrove genus Avicennia. We found that both the number and the proportion of TFs and WRKYs in Avicennia species exceeded their inland relatives, indicating a significant expansion of WRKYs in Avicennia. We identified 109 WRKY genes in the representative species Avicennia marina. Comparative genomic analysis showed that two recent whole-genome duplication (WGD) events played a critical role in the expansion of WRKYs, and 88% of Avicennia marina WRKYs (AmWRKYs) have been retained following these WGDs. Applying comparative transcriptomics on roots under experimental salt gradients, we inferred that there is high divergence in the expression of WGD-retained AmWRKYs. Moreover, we found that the expression of 16 AmWRKYs was stable between freshwater and moderately saline water but increased when the trees were exposed to high salinity. In particular, 14 duplicates were retained following the two recent WGD events, indicating potential neo- and sub-functionalization. We also found that WRKYs could interact with other upregulated genes involved in signalling pathways and natural antioxidant biosynthesis to enhance salt tolerance, contributing to the adaptation to intertidal zones. Our omic data of the WRKY family in A. marina broadens the understanding of how a TF family relates to the adaptive evolution of mangroves.

Supplementary Information

The online version contains supplementary material available at 10.1007/s42995-023-00177-y.

Genome-Wide Effects on Gene Expression Between Parental and Filial Generations of Trisomy 11 and 12 of Rice

Aneuploid refers to the gene dosage imbalance due to copy number alterations. Aneuploidy is generally harmful to the growth, development and reproduction of organisms according to the numerous research. However, it has rarely been reported on whether aneuploid have a relevant pattern of genome expression between the parental and its offspring generations. In this study, mRNA sequencing analysis was performed on rice (Oryza sativa L.) primary trisomes 11 and 12, same primary trisomes and normal individuals in their filial generation. We systematically summarized the changes in gene expression patterns that occur on cis genes and on trans genes between parental and filial generations. In T11 and T12, the ratio of cis-gene expression showed intermediate type in parents and dosage compensation in filial generations, which maybe due to more genes being downregulated. The trans genes were also affected by aneuploidy and manifested as cis-related. The strains with normal chromosomes in filial generations, there are still aneuploid-sensitive genes differentially expressed in their genomes, indicating that the effect of aneuploidy is far-reaching and could not be easily eliminated. Meanwhile, among these differentially expressed genes, genes with low-expression level were more likely to be upregulated, while genes with medium- and high-expression level were easy to be downregulated. For the different types of rice aneuploid, upregulated genes were mainly associated with genomic imbalance while downregulated genes were mainly influenced by the specific added chromosome. In conclusion, our results provide new insights into the genetic characterization and evolution of biological aneuploidy genomes.

Recent genome-wide replication promoted expansion and functional differentiation of the JAZs in soybeans

Jasmonate Zim-domain (JAZ) protein is an inhibitor of the jasmonate (JA) signal transduction pathway, and plays an important role in regulating plant growth, development, and defense. However, there have been few studies on its function under environmental stress in soybeans. In this study, a total of 275 JAZs protein-coding genes were identified in 29 soybean genomes. SoyC13 contained the least JAZ family members (26 JAZs), which was twice as high as AtJAZs. The genes are mainly generated by recent genome-wide replication (WGD), which replicated during the Late Cenozoic Ice Age. In addition, transcriptome analysis showed that the differences in gene expression patterns in the roots, stems, and leaves of the 29 cultivars at the V1 stage were not significant, but there was a significant difference among the three seed development stages. Finally, qRT-PCR results showed that GmJAZs responded the most strongly to heat stress, followed by drought and cold stress. This is consistent with the reason for their expansion and promoter analysis results. Therefore, we explored the significant role of conserved, duplicated, and neofunctionalized JAZs in the evolution of soybeans, which will contribute to the functional characterization of GmJAZ and the improvement of crops.

Enrichment of intrinsically disordered residues in ohnologs facilitates abiotic stress resilience in Brassica rapa

Arabidopsis thaliana and Brassica rapa are in the same evolutionary lineage, although the latter experienced an additional whole genome triplication event. Therefore, it would be intriguing to investigate the traits that gene duplication imposes to mediate plant stress tolerance. Here, we noticed that B. rapa abiotic stress resistance (ASR) genes which code at least one stress responsive domain have a significantly higher number of paralogs than A. thaliana. Analysing the disordered content of the ASR genes in both species, we found that intrinsically disordered residues (IDR) are specifically enriched in whole genome duplication (WGD) derived paralogs. Subsequently, domain similarity analysis between WGD pairs of both species has revealed that majority of WGD pairs in B. rapa did not share domains with each other. Furthermore, domain enrichment analysis has shown that B. rapa paralogs contain 36 distinct stress responsive enriched domains, significantly higher than A. thaliana paralogs. Next, we performed MSA to investigate the domain conservation between orthologs and ohnologs pairs, we found that 80.13% of B. rapa ohnologs acquire new domains, depicting the fact that ohnologs play a significant role in stress-related behaviours. The average IDR content of the ohnologs enriching new domains after gene duplication in B. rapa (0.19), is also significantly higher than A. thaliana (0.04). Interestingly, we also found that all of these attributes i.e., exhibiting higher number of WGD paralogs and enhancement of IDR in ASR genes of B. rapa compared to A. thaliana is exclusive for ASR genes only. No such significant differences were observed in randomly selected non-ASR genes between the two species. Together these results provide strong support for the hypothesis that augmentation of IDR content followed by a whole genome duplication event imposes the stress resistance potentiality in B. rapa. This research will shed light on the mechanism of how B. rapa is able to successfully adapt to stress over the evolutionary timescale.

Acquisition of new function through gene duplication in the metallocarboxypeptidase family

Gene duplication is a key first step in the process of expanding the functionality of a multigene family. In order to better understand the process of gene duplication and its role in the formation of new enzymes, we investigated recent duplication events in the M14 family of proteolytic enzymes. Within vertebrates, four of 23 M14 genes were frequently found in duplicate form. While AEBP1, CPXM1, and CPZ genes were duplicated once through a large-scale, likely whole-genome duplication event, the CPO gene underwent many duplication events within fish and Xenopus lineages. Bioinformatic analyses of enzyme specificity and conservation suggested a greater amount of neofunctionalization and purifying selection in CPO paralogs compared with other CPA/B enzymes. To examine the functional consequences of evolutionary changes on CPO paralogs, the four CPO paralogs from Xenopus tropicalis were expressed in Sf9 and HEK293T cells. Immunocytochemistry showed subcellular distribution of Xenopus CPO paralogs to be similar to that of human CPO. Upon activation with trypsin, the enzymes demonstrated differential activity against three substrates, suggesting an acquisition of new function following duplication and subsequent mutagenesis. Characteristics such as gene size and enzyme activation mechanisms are possible contributors to the evolutionary capacity of the CPO gene.

Comparative gene retention analysis in barley, wild emmer, and bread wheat pangenome lines reveals factors affecting gene retention following gene duplication

BACKGROUND

Gene duplication is a prevalent phenomenon and a major driving force underlying genome evolution. The process leading to the fixation of gene duplicates following duplication is critical to understand how genome evolves but remains fragmentally understood. Most previous studies on gene retention are based on gene duplicate analyses in single reference genome. No population-based comparative gene retention analysis has been performed to date.

RESULTS

Taking advantage of recently published genomic data in Triticeae, we dissected a divergent homogentisate phytyltransferase (HPT2) lineage caught in the middle stage of gene fixation following duplication. The presence/absence of HPT2 in barley (diploid), wild emmer (tetraploid), and bread wheat (hexaploid) pangenome lines appears to be associated with gene dosage constraint and environmental adaption. Based on these observations, we adopted a phylogeny-based orthology inference approach and performed comparative gene retention analyses across barley, wild emmer, and bread wheat. This led to the identification of 326 HPT2-pattern-like genes at whole genome scale, representing a pool of gene duplicates in the middle stage of gene fixation. Majority of these HPT2-pattern-like genes were identified as small-scale duplicates, such as dispersed, tandem, and proximal duplications. Natural selection analyses showed that HPT2-pattern-like genes have experienced relaxed selection pressure, which is generally accompanied with partial positive selection and transcriptional divergence. Functional enrichment analyses showed that HPT2-pattern-like genes are over-represented with molecular-binding and defense response functions, supporting the potential role of environmental adaption during gene retention. We also observed that gene duplicates from larger gene family are more likely to be lost, implying a gene dosage constraint effect. Further comparative gene retention analysis in barley and bread wheat pangenome lines revealed combined effects of species-specific selection and gene dosage constraint.

CONCLUSIONS

Comparative gene retention analyses at the population level support gene dosage constraint, environmental adaption, and species-specific selection as three factors that may affect gene retention following gene duplication. Our findings shed light on the evolutionary process leading to the retention of newly formed gene duplicates and will greatly improve our understanding on genome evolution via duplication.

MicroRNAs play regulatory roles in genomic balance

Classic genetics studies found that genomic imbalance caused by changing the dosage of part of the genome (aneuploidy) has more detrimental effects than altering the dosage of the whole genome (ploidy). Previous analysis revealed global modulation of gene expression triggered by aneuploidy across various species, including maize (Zea mays), Arabidopsis, yeast, mammals, etc. Plant microRNAs (miRNAs) are a class of 20- to 24-nt endogenous small noncoding RNAs that carry out post-transcriptional gene expression regulation. That miRNAs and their putative targets are preferentially retained as duplicates after whole-genome duplication, as are many transcription factors and signaling components, indicates miRNAs are likely to be dosage-sensitive and potentially involved in genomic balance networks. This review addresses the following questions regarding the role of miRNAs in genomic imbalance. (1) How do aneuploidy and polyploidy impact the expression of miRNAs? (2) Do miRNAs play a regulatory role in modulating the expression of their targets under genomic imbalance?

A genome wide association study to dissect the genetic architecture of agronomic traits in Andean lupin (Lupinus mutabilis)

Establishing Lupinus mutabilis as a protein and oil crop requires improved varieties adapted to EU climates. The genetic regulation of strategic breeding traits, including plant architecture, growing cycle length and yield, is unknown. This study aimed to identify associations between 16 669 single nucleotide polymorphisms (SNPs) and 9 agronomic traits on a panel of 223 L. mutabilis accessions, grown in four environments, by applying a genome wide association study (GWAS). Seven environment-specific QTLs linked to vegetative yield, plant height, pods number and flowering time, were identified as major effect QTLs, being able to capture 6 to 20% of the phenotypic variation observed in these traits. Furthermore, two QTLs across environments were identified for flowering time on chromosome 8. The genes FAF, GAMYB and LNK, regulating major pathways involved in flowering and growth habit, as well as GA30X1, BIM1, Dr1, HDA15, HAT3, interacting with these pathways in response to hormonal and environmental cues, were prosed as candidate genes. These results are pivotal to accelerate the development of L. mutabilis varieties adapted to European cropping conditions by using marker-assisted selection (MAS), as well as to provide a framework for further functional studies on plant development and phenology in this species.

The evolutionary patterns, expression profiles, and genetic diversity of expanded genes in barley

Gene duplication resulting from whole-genome duplication (WGD), small-scale duplication (SSD), or unequal hybridization plays an important role in the expansion of gene families. Gene family expansion can also mediate species formation and adaptive evolution. Barley (Hordeum vulgare) is the world's fourth largest cereal crop, and it contains valuable genetic resources due to its ability to tolerate various types of environmental stress. In this study, 27,438 orthogroups in the genomes of seven Poaceae were identified, and 214 of them were significantly expanded in barley. The evolutionary rates, gene properties, expression profiles, and nucleotide diversity between expanded and non-expanded genes were compared. Expanded genes evolved more rapidly and experienced lower negative selection. Expanded genes, including their exons and introns, were shorter, they had fewer exons, their GC content was lower, and their first exons were longer compared with non-expanded genes. Codon usage bias was also lower for expanded genes than for non-expanded genes; the expression levels of expanded genes were lower than those of non-expanded genes, and the expression of expanded genes showed higher tissue specificity than that of non-expanded genes. Several stress-response-related genes/gene families were identified, and these genes could be used to breed barley plants with greater resistance to environmental stress. Overall, our analysis revealed evolutionary, structural, and functional differences between expanded and non-expanded genes in barley. Additional studies are needed to clarify the functions of the candidate genes identified in our study and evaluate their utility for breeding barley plants with greater stress resistance.

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

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

Genome assembly of an Australian native grass species reveals a recent whole-genome duplication and biased gene retention of genes involved in stress response

BACKGROUND

The adaptive significance of polyploidy has been extensively debated, and chromosome-level genome assemblies of polyploids can provide insight into this. The Australian grass Bothriochloa decipiens belongs to the BCD clade, a group with a complex history of hybridization and polyploid. This is the first genome assembly and annotation of a species that belongs to this fascinating yet complex group.

FINDINGS

Using Illumina short reads, 10X Genomics linked reads, and Hi-C sequencing data, we assembled a highly contiguous genome of B. decipiens, with a total length of 1,218.22 Mb and scaffold N50 of 42.637 Mb. Comparative analysis revealed that the species experienced a relatively recent whole-genome duplication. We clustered the 20 major scaffolds, representing the 20 chromosomes, into the 2 subgenomes of the parental species using unique repeat signatures. We found evidence of biased fractionation and differences in the activity of transposable elements between the subgenomes prior to hybridization. Duplicates were enriched for genes involved in transcription and response to external stimuli, supporting a biased retention of duplicated genes following whole-genome duplication.

CONCLUSIONS

Our results support the hypotheses of a biased retention of duplicated genes following polyploidy and point to differences in repeat activity associated with subgenome dominance. B. decipiens is a widespread species with the ability to establish across many soil types, making it a prime candidate for climate change- resilient ecological restoration of Australian grasslands. This reference genome is a valuable resource for future population genomic research on Australian grasses.

Tubby-like proteins (TLPs) transcription factor in different regulatory mechanism in plants: a review

Tubby-like proteins (TLPs) transcription factors are found in single-celled to multi-cellular eukaryotes in the form of large multigene families. TLPs are identified through a specific signature of carboxyl terminal tubby domain, required for plasma membrane tethering and amino terminal F-box domain communicate as functional SCF-type E3 ligases. The comprehensive distribution of TLP gene family members in diverse species indicates some conserved functions of TLPs in multicellular organisms. Plant TLPs have higher gene members than animals and these members reported important role in multiple physiological and developmental processes and various environmental stress responses. Although the TLPs are suggested to be a putative transcription factors but their functional mechanism is not much clear. This review provides significant recent updates on TLP-mediated regulation with an insight into its functional roles, origin and evolution and also phytohormones related regulation to combat with various stresses and its involvement in adaptive stress response in crop plants.

Genome-Wide Analysis of SQUAMOSA-Promoter-Binding Protein-like Family in Flowering Pleioblastus pygmaeus

SQUAMOSA Promoter-Binding Protein-Like (SPL) family is well-known for playing an important role in plant growth and development, specifically in the reproductive process. Bamboo plants have special reproductive characteristics with a prolonged vegetative phase and uncertain flowering time. However, the underlying functions of SPL genes in reproductive growth are undisclosed in bamboo plants. In the study, a total of 28 SPLs were screened from an ornamental dwarf bamboo species, Pleioblastus pygmaeus. Phylogenetic analysis indicates that 183 SPLs from eight plant species can be classified into nine subfamilies, and the 28 PpSPLs are distributed among eight subfamilies. Homologous analysis shows that as many as 32 pairs of homologous genes were found between P. pygmaeus and rice, and 83 pairs were found between P. pygmaeus and Moso bamboo, whose Ka/Ks values are all <1. MiRNA target prediction reveals that 13 out of the 28 PpSPLs have recognition sites complementary to miRNA156. To screen the SPLs involved in the reproductive growth of bamboo plants, the mRNA abundance of the 28 PpSPLs was profiled in the different tissues of flowering P. pygmaeus and non-flowering plants by RNA-Seq. Moreover, the relative expression level of eight PpSPLs is significantly higher in flowering P. pygmaeus than that in non-flowering plants, which was also validated by RT-qPCR. Combined with phylogenetic analysis and homologous analysis, the eight significant, differentially expressed PpSPLs were identified to be associated with the reproductive process and flower organ development. Among them, there are four potential miRNA156-targeting PpSPLs involved in the flowering process. Of significant interest in the study is the identification of 28 SPLs and the exploration of four key flowering-related SPLs from P. pygmaeus, which provides a theoretic basis for revealing the underlying functions of SPLs in the reproductive growth of bamboo plants.

Hybridization order is not the driving factor behind biases in duplicate gene losses among the hexaploid Solanaceae

We model the post-hexaploidy evolution of four genomes from the Solanaceae, a group of flowering plants comprising tomatoes, potatoes and their relatives. The hexaploidy that these genomes descend from occurred through two sequential allopolyploidy events and was marked by the unequal losses of duplicated genes from the different progenitor subgenomes. In contrast with the hexaploid Brassiceae (broccoli and its relatives), where the subgenome with the most surviving genes arrived last in the hexaploidy, among the Solanaceae the most preserved subgenome descends from one of the original two tetraploid progenitors. In fact, the last-arriving subgenome in these plants actually has the fewest surviving genes in the modern genomes. We explore whether the distribution of repetitive elements (REs) in these genomes can explain the biases in gene losses, but while the signals we find are broadly consistent with a role for high RE density in driving gene losses, the REs turn over so quickly that little signal of the RE condition at the time of paleopolyploidy is extant in the modern genomes.

The MKK2a Gene Involved in the MAPK Signaling Cascades Enhances Populus Salt Tolerance

Mitogen-activated protein kinase (MAPK) cascades are highly conserved signal transduction modules, which transmit environmental signals in plant cells through stepwise phosphorylation and play indispensable roles in a wide range of physiological and biochemical processes. Here, we isolated and characterized a gene encoding MKK2 protein from poplar through the rapid amplification of cDNA ends (RACE). The full-length PeMKK2a gene was 1571 bp, including a 1068 bp open reading frame (ORF) encoding 355 amino acids, and the putative PeMKK2a protein belongs to the PKc_like (protein kinase domain) family (70–336 amino acids) in the PKc_MAPKK_plant subfamily and contains 62 sites of possible phosphorylation and two conserved domains, DLK and S/T-xxxxx-S/T. Detailed information about its gene structure, sequence similarities, subcellular localization, and transcript profiles under salt-stress conditions was revealed. Transgenic poplar lines overexpressing PeMKK2a exhibited higher activities of superoxide dismutase (SOD), catalase (CAT), and peroxidase (POD) than non-transgenic poplar under salt stress conditions. These results will provide insight into the roles of MAPK signaling cascades in poplar response to salt stress.

Evolution and function of ubiquitin-specific proteases (UBPs): Insight into seed development roles in tung tree (Vernicia fordii)

The tung oil produced by the tung tree (Vernicia fordii) provides resources for the manufacture of biodiesel. Ubiquitin-specific proteases (UBPs) are the largest group of deubiquitinases and play key roles in regulating development and stress responses. Here, 21 UBPs were identified in V. fordii, roughly one-half the number found in Manihot esculenta and Hevea brasiliensis. Most UBP duplications are produced from whole-genome duplication (WGD), and significant differences in gene retention existed among Euphorbiaceae. The great majority of UBP-containing blocks in V. fordii, V. montana, Ricinus communis, and Jatropha curcas exhibited extensive conservation with the duplicated regions of M. esculenta and H. brasiliensis. These blocks formed 14 orthologous groups, indicating they shared WGD with UBPs in M. esculenta and H. brasiliensis, but most of these UBPs copies were lost. The UBP orthologs contained significant functional divergence which explained the susceptibility of V. fordii to Fusarium wilt and the resistance of V. montana to Fusarium wilt. The expression patterns and experiments suggested that Vf03G1417 could affect the seed-related traits and positively regulate the seed oil accumulation. This study provided important insights into the evolution of UBPs in Euphorbiaceae and identified important candidate VfUBPs for marker-assisted breeding in V. fordii.

Identification and expression analysis of the lipid phosphate phosphatases gene family reveal their involvement in abiotic stress response in kiwifruit

Lipid phosphate phosphatases (LPPs) are a key enzyme in the production and degradation of phosphatidic acid (PA), which plays an important role in plant growth, development, stress resistance and plant hormone response. Thus far, little is known about the LPP family genes in kiwifruit (Actinidia spp.). According to this study, 7 members in the AcLPP family were identified from the whole genome of kiwifruit, the subcellular localization predictions were mainly on the plasma membrane. Chromosomal localization analysis showed that the AcLPP genes were unevenly distributed on 5 chromosomes, it was determined to have undergone strong purifying selection pressure. There were 5 duplicate gene pairs and all underwent segmental duplication events. The LPP genes of kiwifruit were conserved when compared with other plants, especially in terms of evolutionary relationships, conserved motifs, protein sequences, and gene structures. Cis-regulatory elements mainly included hormone response elements and abiotic response elements. Functional annotation of GO revealed that AcLPP genes were closely related to phosphatase/hydrolase activity, phosphorus metabolism and dephosphorylation. AcLPP genes family were predicted to be targets of miRNA. Transcript level analysis revealed that the AcLPP family played diverse functions in different tissues and during growth, development, and postharvest storage stages. qPCR analysis showed that the members of AcLPP gene family might be regulated by ETH, ABA, GA3, and IAA hormone signals. The family members were regulated by the stress of salt stress, osmotic stress, cold stress, and heat stress. These results would provide a basis and reference for studying the agricultural characteristics of kiwifruit and improving its stress resistance.

Engineering Ribosomes to Alleviate Abiotic Stress in Plants: A Perspective

As the centerpiece of the biomass production process, ribosome activity is highly coordinated with environmental cues. Findings revealing ribosome subgroups responsive to adverse conditions suggest this tight coordination may be grounded in the induction of variant ribosome compositions and the differential translation outcomes they might produce. In this perspective, we go through the literature linking ribosome heterogeneity to plants’ abiotic stress response. Once unraveled, this crosstalk may serve as the foundation of novel strategies to custom cultivars tolerant to challenging environments without the yield penalty.

Transcription factor retention through multiple polyploidization steps in wheat

Whole-genome duplication is widespread in plant evolutionary history and is followed by nonrandom gene loss to return to a diploid state. Across multiple angiosperm species, the retained genes tend to be dosage-sensitive regulatory genes such as transcription factors, yet data for younger polyploid species is sparse. Here, we analyzed the retention, expression, and genetic variation in transcription factors in the recent allohexaploid bread wheat (Triticum aestivum L.). By comparing diploid, tetraploid, and hexaploid wheat, we found that, following each of two hybridization and whole-genome duplication events, the proportion of transcription factors in the genome increased. Transcription factors were preferentially retained over other genes as homoeologous groups in tetraploid and hexaploid wheat. Across cultivars, transcription factor homoeologs contained fewer deleterious missense mutations than nontranscription factors, suggesting that transcription factors are maintained as three functional homoeologs in hexaploid wheat populations. Transcription factor homoeologs were more strongly coexpressed than nontranscription factors, indicating conservation of function between homoeologs. We found that the B3, MADS-M-type, and NAC transcription factor families were less likely to have three homoeologs present than other families, which was associated with low expression levels and high levels of tandem duplication. Together, our results show that transcription factors are preferentially retained in polyploid wheat genomes although there is variation between families. Knocking out one transcription factor homoeolog to alter gene dosage, using TILLING or CRISPR, could generate new phenotypes for wheat breeding.

The multiple fates of gene duplications: Deletion, hypofunctionalization, subfunctionalization, neofunctionalization, dosage balance constraints, and neutral variation

Gene duplications have long been recognized as a contributor to the evolution of genes with new functions. Multiple copies of genes can result from tandem duplication, from transposition to new chromosomes, or from whole-genome duplication (polyploidy). The most common fate is that one member of the pair is deleted to return the gene to the singleton state. Other paths involve the reduced expression of both copies (hypofunctionalization) that are held in duplicate to maintain sufficient quantity of function. The two copies can split functions (subfunctionalization) or can diverge to generate a new function (neofunctionalization). Retention of duplicates resulting from doubling of the whole genome occurs for genes involved with multicomponent interactions such as transcription factors and signal transduction components. In contrast, these classes of genes are underrepresented in small segmental duplications. This complementary pattern suggests that the balance of interactors affects the fate of the duplicate pair. We discuss the different mechanisms that maintain duplicated genes, which may change over time and intersect.