Gene and genome duplications: the impact of dosage-sensitivity on the fate of nuclear genes

Genome-Wide Identification, Expression, and Protein Analysis of CKX and IPT Gene Families in Radish (Raphanus sativus L.) Reveal Their Involvement in Clubroot Resistance

Cytokinins (CKs) are a group of phytohormones that are involved in plant growth, development, and disease resistance. The isopentenyl transferase (IPT) and cytokinin oxidase/dehydrogenase (CKX) families comprise key enzymes controlling CK biosynthesis and degradation. However, an integrated analysis of these two gene families in radish has not yet been explored. In this study, 13 RsIPT and 12 RsCKX genes were identified and characterized, most of which had four copies in Brassica napus and two copies in radish and other diploid Brassica species. Promoter analysis indicated that the genes contained at least one phytohormone or defense and stress responsiveness cis-acting element. RsIPTs and RsCKXs were expanded through segmental duplication. Moreover, strong purifying selection drove the evolution of the two gene families. The expression of the RsIPT and RsCKX genes distinctly showed diversity in different tissues and developmental stages of the root. Expression profiling showed that RsCKX1-1/1-2/1-3 was significantly upregulated in club-resistant materials during primary infection, suggesting their vital function in clubroot resistance. The interaction network of CKX proteins with similar 3D structures also reflected the important role of RsCKX genes in disease resistance. This study provides a foundation for further functional study on the IPT and CKX genes for clubroot resistance improvement in Raphanus.

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.

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

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

Comprehensive identification of maize ZmE2F transcription factors and the positive role of ZmE2F6 in response to drought stress

BACKGROUND

The early 2 factor (E2F) family is characterized as a kind of transcription factor that plays an important role in cell division, DNA damage repair, and cell size regulation. However, its stress response has not been well revealed.

RESULTS

In this study, ZmE2F members were comprehensively identified in the maize genome, and 21 ZmE2F genes were identified, including eight E2F subclade members, seven DEL subfamily genes, and six DP genes. All ZmE2F proteins possessed the DNA-binding domain (DBD) characterized by conserved motif 1 with the RRIYD sequence. The ZmE2F genes were unevenly distributed on eight maize chromosomes, showed diversity in gene structure, expanded by gene duplication, and contained abundant stress-responsive elements in their promoter regions. Subsequently, the ZmE2F6 gene was cloned and functionally verified in drought response. The results showed that the ZmE2F6 protein interacted with ZmPP2C26, localized in the nucleus, and responded to drought treatment. The overexpression of ZmE2F6 enhanced drought tolerance in transgenic Arabidopsis with longer root length, higher survival rate, and biomass by upregulating stress-related gene transcription.

CONCLUSIONS

This study provides novel insights into a greater understanding and functional study of the E2F family in the stress response.

Adaptation in Unstable Environments and Global Gene Losses: Small but Stable Gene Networks by the May-Wigner Theory

Although gene loss is common in evolution, it remains unclear whether it is an adaptive process. In a survey of seven major mangrove clades that are woody plants in the intertidal zones of daily environmental perturbations, we noticed that they generally evolved reduced gene numbers. We then focused on the largest clade of Rhizophoreae and observed the continual gene set reduction in each of the eight species. A great majority of gene losses are concentrated on environmental interaction processes, presumably to cope with the constant fluctuations in the tidal environments. Genes of the general processes for woody plants are largely retained. In particular, fewer gene losses are found in physiological traits such as viviparous seeds, high salinity, and high tannin content. Given the broad and continual genome reductions, we propose the May-Wigner theory (MWT) of system stability as a possible mechanism. In MWT, the most effective solution for buffering continual perturbations is to reduce the size of the system (or to weaken the total genic interactions). Mangroves are unique as immovable inhabitants of the compound environments in the land-sea interface, where environmental gradients (such as salinity) fluctuate constantly, often drastically. Extending MWT to gene regulatory network (GRN), computer simulations and transcriptome analyses support the stabilizing effects of smaller gene sets in mangroves vis-à-vis inland plants. In summary, we show the adaptive significance of gene losses in mangrove plants, including the specific role of promoting phenotype innovation and a general role in stabilizing GRN in unstable environments as predicted by MWT.

Genome-wide analysis of cellulose synthase gene superfamily in Tectona grandis L.f

This study aimed to explore Cellulose synthase gene superfamily of teak, and its evolutionary relationship with homologous genes of other woody species. The incidence of evolutionary events like gene duplication and gene loss, influence of the selection pressure, and consequent adaptive functional divergence of the duplicated TgCes gene were assessed alongside it's role in wood coloration. This study identified 39 full-length non-redundant proteins belonging to CesA and Csl gene families. TgCesA and TgCsl proteins with Cellulose synthase domain repeats indicated tandem gene duplication and probable genetic variability, enabling local adaptation. Further, multi-domain protein (MYB-like DNA-binding domain and CesA domain) with maximum introns was also identified indicating gene fusion and formation of complex protein with novel functions. Phylogenetic analysis grouped the genes into seven subfamilies (CesA, CslA, CslC, CslD, CslE, CslG, and CslM) with each undergoing gene duplication and loss along their evolutionary history. Post-species gene duplications and probable neofunctionalization were identified in TgCesA and TgCsl gene families. Each subfamily was found to be under strong purifying selection with a few or no sites under positive selection. Functional divergence analysis further revealed site-specific selective constraints in CesA and Csl genes of the teak Cellulose synthase gene family. Furthermore, protein-protein interaction network analysis identified co-expression of Cellulose synthase gene with flavonoid 3',5'-hydroxylase (F3'5'H, CYP75A), involved in the biosynthesis of xylem anthocyanin compounds, probably responsible for wood coloration. This study thus offers a foundation for future research in wood formation and wood property traits specific to teak and its provenances.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-024-03927-6.

The self-compatibility is acquired after polyploidization: a case study of Brassica napus self-incompatible trilinear hybrid breeding system

Self-incompatibility plays a vital role in angiosperms, by preventing inbreeding depression and maintaining genetic diversity within populations. Following polyploidization, many angiosperm species transition from self-incompatibility to self-compatibility. Here, we investigated the S-locus in Brassicaceae and identified distinct origins for the sRNA loci, SMI and SMI2 (SCR Methylation Inducer 1 and 2), within the S-locus. The SMI loci were found to be widespread in Cruciferae, whereas the SMI2 loci were exclusive to Brassica species. Additionally, we discovered four major S-haplotypes (BnS-1, BnS-6, BnS-7, and BnS-1300) in rapeseed. Overexpression of BnSMI-1 in self-incompatible Brassica napus ('S-70S1300S6 ') resulted in a significant increase in DNA methylation in the promoter regions of BnSCR-6 and BnSCR-1300, leading to self-compatibility. Conversely, by overexpressing a point mutation of BnSmi-1 in the 'S-70S1300S6 ' line, we observed lower levels of DNA methylation in BnSCR-6 and BnSCR-1300 promoters. Furthermore, the overexpression of BnSMI2-1300 in the 'SI-326S7S6 ' line inhibited the expression of BnSCR-7 through transcriptional repression of the Smi2 sRNA from the BnS-1300 haplotype. Our study demonstrates that the self-compatibility of rapeseed is determined by S-locus sRNA-mediated silencing of SCR after polyploidization, which helps to further breed self-incompatible or self-compatible rapeseed lines, thereby facilitating the utilization of heterosis.

The genetic and molecular basis of haploinsufficiency in flowering plants

In diploid organisms, haploinsufficiency can be defined as the requirement for more than one fully functional copy of a gene. In contrast to most genes, whose loss-of-function alleles are recessive, loss-of-function alleles of haploinsufficient genes are dominant. However, forward and reverse genetic screens are biased toward obtaining recessive, loss-of-function mutations, and therefore, dominant mutations of all types are underrepresented in mutant collections. Despite this underrepresentation, haploinsufficient loci have intriguing implications for studies of genome evolution, gene dosage, stability of protein complexes, genetic redundancy, and gene expression. Here we review examples of haploinsufficiency in flowering plants and describe the underlying molecular mechanisms and evolutionary forces driving haploinsufficiency. Finally, we discuss the masking of haploinsufficiency by genetic redundancy, a widespread phenomenon among angiosperms.

Maternal dominance contributes to subgenome differentiation in allopolyploid fishes

Teleost fishes, which are the largest and most diverse group of living vertebrates, have a rich history of ancient and recent polyploidy. Previous studies of allotetraploid common carp and goldfish (cyprinids) reported a dominant subgenome, which is more expressed and exhibits biased gene retention. However, the underlying mechanisms contributing to observed 'subgenome dominance' remains poorly understood. Here we report high-quality genomes of twenty-one cyprinids to investigate the origin and subsequent subgenome evolution patterns following three independent allopolyploidy events. We identify the closest extant relatives of the diploid progenitor species, investigate genetic and epigenetic differences among subgenomes, and conclude that observed subgenome dominance patterns are likely due to a combination of maternal dominance and transposable element densities in each polyploid. These findings provide an important foundation to understanding subgenome dominance patterns observed in teleost fishes, and ultimately the role of polyploidy in contributing to evolutionary innovations.

Genome-Wide Identification of BrCAX Genes and Functional Analysis of BrCAX1 Involved in Ca2+ Transport and Ca2+ Deficiency-Induced Tip-Burn in Chinese Cabbage (Brassica rapa L. ssp. pekinensis)

Calcium (Ca2+) plays essential roles in plant growth and development. Ca2+ deficiency causes a physiological disorder of tip-burn in Brassiceae crops and is involved in the regulation of cellular Ca2+ homeostasis. Although the functions of Ca2+/H+ exchanger antiporters (CAXs) in mediating transmembrane transport of Ca2+ have been extensively characterized in multiple plant species, the potential roles of BrCAX genes remain unclear in Chinese cabbage. In this study, eight genes of the BrCAX family were genome-widely identified in Chinese cabbage. These BrCAX proteins contained conserved Na_Ca_ex domain and belonged to five members of the CAX family. Molecular evolutionary analysis and sequence alignment revealed the evolutionary conservation of BrCAX family genes. Expression profiling demonstrated that eight BrCAX genes exhibited differential expression in different tissues and under heat stress. Furthermore, Ca2+ deficiency treatment induced the typical symptoms of tip-burn in Chinese cabbage seedlings and a significant decrease in total Ca2+ content in both roots and leaves. The expression changes in BrCAX genes were related to the response to Ca2+ deficiency-induced tip-burn of Chinese cabbage. Specially, BrCAX1-1 and BrCAX1-2 genes were highly expressed gene members of the BrCAX family in the leaves and were significantly differentially expressed under Ca2+ deficiency stress. Moreover, overexpression of BrCAX1-1 and BrCAX1-2 genes in yeast and Chinese cabbage cotyledons exhibited a higher Ca2+ tolerance, indicating the Ca2+ transport capacity of BrCAX1-1 and BrCAX1-2. In addition, suppression expression of BrCAX1-1 and BrCAX1-2 genes reduced cytosolic Ca2+ levels in the root tips of Chinese cabbage. These results provide references for functional studies of BrCAX genes and to investigate the regulatory mechanisms underlying Ca2+ deficiency disorder in Brassiceae vegetables.

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.

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.

Functional and evolutionary study of MLO gene family in the regulation of Sclerotinia stem rot resistance in Brassica napus L

BACKGROUND

Oilseed rape (Brassica napus L.) is known as one of the most important oilseed crops cultivated around the world. However, its production continuously faces a huge challenge of Sclerotinia stem rot (SSR), a destructive disease caused by the fungus Sclerotinia sclerotiorum, resulting in huge yield loss annually. The SSR resistance in B. napus is quantitative and controlled by a set of minor genes. Identification of these genes and pyramiding them into a variety are a major strategy for SSR resistance breeding in B. napus.

RESULTS

Here, we performed a genome-wide association study (GWAS) using a natural population of B. napus consisting of 222 accessions to identify BnaA08g25340D (BnMLO2_2) as a candidate gene that regulates the SSR resistance. BnMLO2_2 was a member of seven homolog genes of Arabidopsis Mildew Locus O 2 (MLO2) and the significantly SNPs were mainly distributed in the promoter of BnMLO2_2, suggesting a role of BnMLO2_2 expression level in the regulation of SSR resistance. We expressed BnMLO2_2 in Arabidopsis and the transgenic plants displayed an enhanced SSR resistance. Transcriptome profiling of different tissues of B. napus revealed that BnMLO2_2 had the most expression level in leaf and silique tissues among all the 7 BnMLO2 members and also expressed higher in the SSR resistant accession than in the susceptible accession. In Arabidopsis, mlo2 plants displayed reduced resistance to SSR, whereas overexpression of MLO2 conferred plants an enhanced SSR resistance. Moreover, a higher expression level of MLO2 showed a stronger SSR resistance in the transgenic plants. The regulation of MLO2 in SSR resistance may be associated with the cell death. Collinearity and phylogenetic analysis revealed a large expansion of MLO family in Brassica crops.

CONCLUSION

Our study revealed an important role of BnMLO2 in the regulation of SSR resistance and provided a new gene candidate for future improvement of SSR resistance in B. napus and also new insights into understanding of MLO family evolution in Brassica crops.

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.

Dosage balance acts as a time-dependent selective barrier to subfunctionalization

BACKGROUND

Gene duplication is an important process for genome expansion, sometimes allowing for new gene functions to develop. Duplicate genes can be retained through multiple processes, either for intermediate periods of time through processes such as dosage balance, or over extended periods of time through processes such as subfunctionalization and neofunctionalization.

RESULTS

Here, we built upon an existing subfunctionalization Markov model by incorporating dosage balance to describe the interplay between subfunctionalization and dosage balance to explore selective pressures on duplicate copies. Our model incorporates dosage balance using a biophysical framework that penalizes the fitness of genetic states with stoichiometrically imbalanced proteins. These imbalanced states cause increased concentrations of exposed hydrophobic surface areas, which cause deleterious mis-interactions. We draw comparison between our Subfunctionalization + Dosage-Balance Model (Sub + Dos) and the previous Subfunctionalization-Only (Sub-Only) Model. This comparison includes how the retention probabilities change over time, dependent upon the effective population size and the selective cost associated with spurious interaction of dosage-imbalanced partners. We show comparison between Sub-Only and Sub + Dos models for both whole-genome duplication and small-scale duplication events.

CONCLUSION

These comparisons show that following whole-genome duplication, dosage balance serves as a time-dependent selective barrier to the subfunctionalization process, by causing an overall delay but ultimately leading to a larger portion of the genome retained through subfunctionalization. This higher percentage of the genome that is ultimately retained is caused by the alternative competing process, nonfunctionalization, being selectively blocked to a greater extent. In small-scale duplication, the reverse pattern is seen, where dosage balance drives faster rates of subfunctionalization, but ultimately leads to a smaller portion of the genome retained as duplicates. This faster rate of subfunctionalization is because the dosage balance of interacting gene products is negatively affected immediately after duplication and the loss of a duplicate restores the stoichiometric balance. Our findings provide support that the subfunctionalization of genes that are susceptible to dosage balance effects, such as proteins involved in complexes, is not a purely neutral process. With stronger selection against stoichiometrically imbalanced gene partners, the rates of subfunctionalization and nonfunctionalization slow; however, this ultimately leads to a greater proportion of subfunctionalized gene pairs.

Dynamics of Gene Loss following Ancient Whole-Genome Duplication in the Cryptic Paramecium Complex

Whole-genome duplications (WGDs) have shaped the gene repertoire of many eukaryotic lineages. The redundancy created by WGDs typically results in a phase of massive gene loss. However, some WGD-derived paralogs are maintained over long evolutionary periods, and the relative contributions of different selective pressures to their maintenance are still debated. Previous studies have revealed a history of three successive WGDs in the lineage of the ciliate Paramecium tetraurelia and two of its sister species from the Paramecium aurelia complex. Here, we report the genome sequence and analysis of 10 additional P. aurelia species and 1 additional out group, revealing aspects of post-WGD evolution in 13 species sharing a common ancestral WGD. Contrary to the morphological radiation of vertebrates that putatively followed two WGD events, members of the cryptic P. aurelia complex have remained morphologically indistinguishable after hundreds of millions of years. Biases in gene retention compatible with dosage constraints appear to play a major role opposing post-WGD gene loss across all 13 species. In addition, post-WGD gene loss has been slower in Paramecium than in other species having experienced genome duplication, suggesting that the selective pressures against post-WGD gene loss are especially strong in Paramecium. A near complete lack of recent single-gene duplications in Paramecium provides additional evidence for strong selective pressures against gene dosage changes. This exceptional data set of 13 species sharing an ancestral WGD and 2 closely related out group species will be a useful resource for future studies on Paramecium as a major model organism in the evolutionary cell biology.

Evolution and expression analysis of the caffeoyl-CoA 3-O-methyltransferase (CCoAOMT) gene family in jute (Corchorus L.)

BACKGROUND

Jute is considered one of the most important crops for fiber production and multipurpose usages. Caffeoyl-CoA 3-O-methyltransferase (CCoAOMT) is a crucial enzyme involved in lignin biosynthesis in plants. The potential functions of CCoAOMT in lignin biosynthesis of jute have been reported in several studies. However, little is known about the evolution of the CCoAOMT gene family, and either their expression level at different developing stages in different jute cultivars, as well as under abiotic stresses including salt and drought stress.

RESULTS

In the present study, 66 CCoAOMT genes from 12 species including 12 and eight CCoAOMTs in Corchorus olitorius and C. capsularis were identified. Phylogenetic analysis revealed that CCoAOMTs could be divided into six groups, and gene expansion was observed in C. olitorius. Furthermore, gene expression analysis of developing jute fibers was conducted at different developmental stages (15, 30, 45, 60, and 90 days after sowing [DAS]) in six varieties (Jute-179 [J179], Lubinyuanguo [LB], and Qiongyueqing [QY] for C. capsularis; Funong No.5 [F5], Kuanyechangguo [KY], and Cvlv [CL] for C. olitorius). The results showed that CCoAOMT1 and CCoAOMT2 were the dominant genes in the CCoAOMT family. Of these two dominant CCoAOMTs, CCoAOMT2 showed a constitutive expression level during the entire growth stages, while CCoAOMT1 exhibited differential expression patterns. These two genes showed higher expression levels in C. olitorius than in C. capsularis. The correlation between lignin content and CCoAOMT gene expression levels indicated that this gene family influences the lignin content of jute. Using real-time quantitative reverse transcription PCR (qRT-PCR), a substantial up-regulation of CCoAOMTs was detected in stem tissues of jute 24 h after drought treatment, with an up to 17-fold increase in expression compared to that of untreated plants.

CONCLUSIONS

This study provides a basis for comprehensive genomic studies of the entire CCoAOMT gene family in C. capsularis and C. olitorius. Comparative genomics analysis among the CCoAOMT gene families of 12 species revealed the close evolutionary relationship among Corchorus, Theobroma cacao and Gossypium raimondii. This study also shows that CCoAOMTs are not only involved in lignin biosynthesis, but also are associated with the abiotic stress response in jute, and suggests the potential use of these lignin-related genes to genetically improve the fiber quality of jute.

Molecular mechanisms of adaptive evolution in wild animals and plants

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

Chromosomal-level genome assembly of Melastoma candidum provides insights into trichome evolution

Melastoma, consisting of ~100 species diversified in tropical Asia and Oceania in the past 1-2 million years, represents an excellent example of rapid speciation in flowering plants. Trichomes on hypanthia, twigs and leaves vary markedly among species of this genus and are the most important diagnostic traits for species identification. These traits also play critical roles in contributing to differential adaptation of these species to their own habitats. Here we sequenced the genome of M. candidum, a common, erect-growing species from southern China, with the aim to provide genomic insights into trichome evolution in this genus. We generated a high-quality, chromosome-level genome assembly of M. candidum, with the genome size of 256.2 Mb and protein-coding gene number of 40,938. The gene families specific to, and significantly expanded in Melastoma are enriched for GO terms related to trichome initiation and differentiation. We provide evidence that Melastoma and its sister genus Osbeckia have undergone two whole genome duplications (WGDs) after the triplication event (γ) shared by all core eudicots. Preferential retention of trichome development-related transcription factor genes such as C2H2, bHLH, HD-ZIP, WRKY, and MYB after both WGDs might provide raw materials for trichome evolution and thus contribute to rapid species diversification in Melastoma. Our study provides candidate transcription factor genes related to trichome evolution in Melastoma, which can be used to evolutionary and functional studies of trichome diversification among species of this genus.

The genome sequence of Hirschfeldia incana, a new Brassicaceae model to improve photosynthetic light-use efficiency

Photosynthesis is a key process in sustaining plant and human life. Improving the photosynthetic capacity of agricultural crops is an attractive means to increase their yields. While the core mechanisms of photosynthesis are highly conserved in C₃ plants, these mechanisms are very flexible, allowing considerable diversity in photosynthetic properties. Among this diversity is the maintenance of high photosynthetic light-use efficiency at high irradiance as identified in a small number of exceptional C₃ species. Hirschfeldia incana, a member of the Brassicaceae family, is such an exceptional species, and because it is easy to grow, it is an excellent model for studying the genetic and physiological basis of this trait. Here, we present a reference genome of H. incana and confirm its high photosynthetic light-use efficiency. While H. incana has the highest photosynthetic rates found so far in the Brassicaceae, the light-saturated assimilation rates of closely related Brassica rapa and Brassica nigra are also high. The H. incana genome has extensively diversified from that of B. rapa and B. nigra through large chromosomal rearrangements, species-specific transposon activity, and differential retention of duplicated genes. Duplicated genes in H. incana, B. rapa, and B. nigra that are involved in photosynthesis and/or photoprotection show a positive correlation between copy number and gene expression, providing leads into the mechanisms underlying the high photosynthetic efficiency of these species. Our work demonstrates that the H. incana genome serves as a valuable resource for studying the evolution of high photosynthetic light-use efficiency and enhancing photosynthetic rates in crop species.

A Chromosome-Scale Genome Assembly of Mitragyna speciosa (Kratom) and the Assessment of Its Genetic Diversity in Thailand

Mitragyna speciosa (Kratom) is a tropical narcotic plant native to Southeast Asia with unique pharmacological properties. Here, we report the first chromosome-scale assembly of the M. speciosa genome. We employed PacBio sequencing to obtain a preliminary assembly, which was subsequently scaffolded using the chromatin contact mapping technique (Hi-C) into 22 pseudomolecules. The final assembly was 692 Mb with a scaffold N50 of 26 Mb. We annotated a total of 39,708 protein-coding genes, and our gene predictions recovered 98.4% of the highly conserved orthologs based on the BUSCO analysis. The phylogenetic analysis revealed that M. speciosa diverged from the last common ancestors of Coffea arabica and Coffea canephora approximately 47.6 million years ago. Our analysis of the sequence divergence at fourfold-degenerate sites from orthologous gene pairs provided evidence supporting a genome-wide duplication in M. speciosa, agreeing with the report that members of the genus Mitragyna are tetraploid. The STRUCTURE and principal component analyses demonstrated that the 85 M. speciosa accessions included in this study were an admixture of two subpopulations. The availability of our high-quality chromosome-level genome assembly and the transcriptomic resources will be useful for future studies on the alkaloid biosynthesis pathway, as well as comparative phylogenetic studies in Mitragyna and related species.

Genome-wide characterization of ovate family protein gene family associated with number of seeds per silique in Brassica napus

Ovate family proteins (OFPs) were firstly identified in tomato as proteins controlling the pear shape of the fruit. Subsequent studies have successively proved that OFPs are a class of negative regulators of plant development, and are involved in the regulation of complex traits in different plants. However, there has been no report about the functions of OFPs in rapeseed growth to date. Here, we identified the OFPs in rapeseed at the genomic level. As a result, a total of 67 members were obtained. We then analyzed the evolution from Arabidopsis thaliana to Brassica napus, illustrated their phylogenetic and syntenic relationships, and compared the gene structure and conserved domains between different copies. We also analyzed their expression patterns in rapeseed, and found significant differences in the expression of different members and in different tissues. Additionally, we performed a GWAS for the number of seeds per silique (NSPS) in a rapeseed population consisting of 204 natural accessions, and identified a new gene BnOFP13_2 significantly associated with NSPS, which was identified as a novel function of OFPs. Haplotype analysis revealed that the accessions with haplotype 3 had a higher NSPS than other accessions, suggesting that BnOFP13_2 is associated with NSPS. Transcript profiling during the five stages of silique development demonstrated that BnOFP13_2 negatively regulates NSPS. These findings provide evidence for functional diversity of OFP gene family and important implications for oilseed rape breeding.

A temporal gradient of cytonuclear coordination of chaperonins and chaperones during RuBisCo biogenesis in allopolyploid plants

Ribulose-1,5-bisphosphate carboxylase/oxygenase (RuBisCo), consisting of subunits encoded by nuclear and cytoplasmic genes, is a model for cytonuclear evolution in plant allopolyploids. To date, coordinated cytonuclear evolutionary responses of auxiliary cofactors involved in RuBisCo biogenesis remain unexplored. This study characterized and compared genomic and transcriptional cytonuclear coevolutionary responses of chaperonin/chaperones in RuBisCo folding and assembly processes across different allopolyploids. We discovered significant cytonuclear evolutionary responses in folding cofactors, with diminishing or attenuated responses later during assembly. Our results have general significance for understanding the unrecognized cytonuclear evolution of chaperonin/chaperone genes, structural and functional features of intermediate complexes, and the functioning stage of the Raf2 cofactor. Generally, the results reveal a hitherto unexplored dimension of allopolyploidy in plants.

Ribulose-1,5-bisphosphate carboxylase/oxygenase (RuBisCo) has long been studied from many perspectives. As a multisubunit (large subunits [LSUs] and small subunits[SSUs]) protein encoded by genes residing in the chloroplast (rbcL) and nuclear (rbcS) genomes, RuBisCo also is a model for cytonuclear coevolution following allopolyploid speciation in plants. Here, we studied the genomic and transcriptional cytonuclear coordination of auxiliary chaperonin and chaperones that facilitate RuBisCo biogenesis across multiple natural and artificially synthesized plant allopolyploids. We found similar genomic and transcriptional cytonuclear responses, including respective paternal-to-maternal conversions and maternal homeologous biased expression, in chaperonin/chaperon-assisted folding and assembly of RuBisCo in different allopolyploids. One observation is about the temporally attenuated genomic and transcriptional cytonuclear evolutionary responses during early folding and later assembly process of RuBisCo biogenesis, which were established by long-term evolution and immediate onset of allopolyploidy, respectively. Our study not only points to the potential widespread and hitherto unrecognized features of cytonuclear evolution but also bears implications for the structural interaction interface between LSU and Cpn60 chaperonin and the functioning stage of the Raf2 chaperone.

Polycomb-mediated repression of paternal chromosomes maintains haploid dosage in diploid embryos of Marchantia

Complex mechanisms regulate gene dosage throughout eukaryotic life cycles. Mechanisms controlling gene dosage have been extensively studied in animals, however it is unknown how generalizable these mechanisms are to diverse eukaryotes. Here, we use the haploid plant Marchantia polymorpha to assess gene dosage control in its short-lived diploid embryo. We show that throughout embryogenesis, paternal chromosomes are repressed resulting in functional haploidy. The paternal genome is targeted for genomic imprinting by the Polycomb mark H3K27me3 starting at fertilization, rendering the maternal genome in control of embryogenesis. Maintaining haploid gene dosage by this new form of imprinting is essential for embryonic development. Our findings illustrate how haploid-dominant species can regulate gene dosage through paternal chromosome inactivation and initiates the exploration of the link between life cycle history and gene dosage in a broader range of organisms.

Analysis of flavonol regulator evolution in the Brassicaceae reveals MYB12, MYB111 and MYB21 duplications and MYB11 and MYB24 gene loss

BACKGROUND

Flavonols are the largest subgroup of flavonoids, possessing multiple functions in plants including protection against ultraviolet radiation, antimicrobial activities, and flower pigmentation together with anthocyanins. They are of agronomical and economical importance because the major off-taste component in rapeseed protein isolates is a flavonol derivative, which limits rapeseed protein use for human consumption. Flavonol production in Arabidopsis thaliana is mainly regulated by the subgroup 7 (SG7) R2R3-MYB transcription factors MYB11, MYB12, and MYB111. Recently, the SG19 MYBs MYB21, MYB24, and MYB57 were shown to regulate flavonol accumulation in pollen and stamens. The members of each subgroup are closely related, showing gene redundancy and tissue-specific expression in A. thaliana. However, the evolution of these flavonol regulators inside the Brassicaceae, especially inside the Brassiceae, which include the rapeseed crop species, is not fully understood.

RESULTS

We studied the SG7 and SG19 MYBs in 44 species, including 31 species of the Brassicaceae, by phylogenetic analyses followed by synteny and gene expression analyses. Thereby we identified a deep MYB12 and MYB111 duplication inside the Brassicaceae, which likely occurred before the divergence of Brassiceae and Thelypodieae. These duplications of SG7 members were followed by the loss of MYB11 after the divergence of Eruca vesicaria from the remaining Brassiceae species. Similarly, MYB21 experienced duplication before the emergence of the Brassiceae tribe, where the gene loss of MYB24 is also proposed to have happened. The members of each subgroup revealed frequent overlapping spatio-temporal expression patterns in the Brassiceae member B. napus, which are assumed to compensate for the loss of MYB11 and MYB24 in the analysed tissues.

CONCLUSIONS

We identified a duplication of MYB12, MYB111, and MYB21 inside the Brassicaceae and MYB11 and MYB24 gene loss inside the tribe Brassiceae. We propose that polyploidization events have shaped the evolution of the flavonol regulators in the Brassicaceae, especially in the Brassiceae.

Convergent evolution of polyploid genomes from across the eukaryotic tree of life

By modeling the homoeologous gene losses that occurred in 50 genomes deriving from ten distinct polyploidy events, we show that the evolutionary forces acting on polyploids are remarkably similar, regardless of whether they occur in flowering plants, ciliates, fishes, or yeasts. We show that many of the events show a relative rate of duplicate gene loss before the first postpolyploidy speciation that is significantly higher than in later phases of their evolution. The relatively weak selective constraint experienced by the single-copy genes these losses produced leads us to suggest that most of the purely selectively neutral duplicate gene losses occur in the immediate postpolyploid period. Nearly all of the events show strong evidence of biases in the duplicate losses, consistent with them being allopolyploidies, with 2 distinct progenitors contributing to the modern species. We also find ongoing and extensive reciprocal gene losses (alternative losses of duplicated ancestral genes) between these genomes. With the exception of a handful of closely related taxa, all of these polyploid organisms are separated from each other by tens to thousands of reciprocal gene losses. As a result, it is very unlikely that viable diploid hybrid species could form between these taxa, since matings between such hybrids would tend to produce offspring lacking essential genes. It is, therefore, possible that the relatively high frequency of recurrent polyploidies in some lineages may be due to the ability of new polyploidies to bypass reciprocal gene loss barriers.

Characterization of silicon transporter gene family in Saccharum and functional analysis of the ShLsi6 gene in biotic stress

Silicon, one of the most prevalent elements in the soil, is beneficial for plant growth and defense against different stresses. The silicon transporter gene (Lsi) plays an important role in the uptake and transport of silicon in higher plants. In this study, a total of 32 Lsi genes, including 20 SsLsi in sugarcane wild species Saccharum spontaneum, 5 ShLsi in Saccharum hybrid cultivar R570 and 7 SbLsi in sugarcane related species Sorghum bicolor, were identified and classified into three groups. Bioinformatics analysis showed that instability, hydrophobicity, localization of cell membranes and vacuoles were the main features of the Lsi proteins. Whole genome and segmental duplication contributed to the main expansion of Lsi gene family. Collinearity analysis of the Lsi genes showed that S. spontanum and R570 had a collinear relationship with monocotyledonous plants S. bicolor and Oryza sativa, but not with dicotyledonous plants Arabidopsis thaliana and Vitis vinifera. The replicated Lsi genes were mainly subjected to strong selection pressure for purification. The diverse cis-regulatory elements in the promoter of SsLsi, ShLsi and SbLsi genes suggested that they were widely involved in the response of plants to various stresses and the regulation of the growth and development. Transcriptome data and real time quantitative PCR analysis showed that the Lsi genes exhibited different expression profiles in sugarcane tissues and under Sporisorium scitamineum, drought and cold stresses. In addition, the cDNA and genomic DNA sequences of ShLsi6 that was homologous to SsLsi1b gene was cloned from Saccharum hybrid cultivar ROC22. Transient expression analysis showed that, compared with the control, Nicotiana benthamiana leaves which overexpressed the ShLsi6 gene showed a high sensitivity after inoculation with tobacco pathogens Ralstonia solanacearum and Fusarium solani var. coeruleum. This study provides important information for further functional analysis of Lsi genes and resistant breeding in sugarcane.

CHD chromatin remodeling protein diversification yields novel clades and domains absent in classic model organisms

Chromatin remodelers play a fundamental role in the assembly of chromatin, regulation of transcription, and DNA repair. Biochemical and functional characterizations of the CHD family of chromatin remodelers from a variety of model organisms have shown that these remodelers participate in a wide range of activities. However, because the evolutionary history of CHD homologs is unclear, it is difficult to predict which of these activities are broadly conserved and which have evolved more recently in individual eukaryotic lineages. Here, we performed a comprehensive phylogenetic analysis of 8,042 CHD homologs from 1,894 species to create a model for the evolution of this family across eukaryotes with a particular focus on the timing of duplications that gave rise to the diverse copies observed in plants, animals, and fungi. Our analysis confirms that the three major subfamilies of CHD remodelers originated in the eukaryotic last common ancestor, and subsequent losses occurred independently in different lineages. Improved taxon sampling identified several subfamilies of CHD remodelers in plants that were absent or highly divergent in the model plant Arabidopsis thaliana. Whereas the timing of CHD subfamily expansions in vertebrates corresponds to whole genome duplication events, the mechanisms underlying CHD diversification in land plants appear more complicated. Analysis of protein domains reveals that CHD remodeler diversification has been accompanied by distinct transitions in domain architecture, contributing to the functional differences observed between these remodelers. This study demonstrates the importance of proper taxon sampling when studying ancient evolutionary events to prevent misinterpretation of subsequent lineage-specific changes and provides an evolutionary framework for functional and comparative analysis of this critical chromatin remodeler family across eukaryotes.

Evolution of binding preferences among whole-genome duplicated transcription factors

Throughout evolution, new transcription factors (TFs) emerge by gene duplication, promoting growth and rewiring of transcriptional networks. How TF duplicates diverge was studied in a few cases only. To provide a genome-scale view, we considered the set of budding yeast TFs classified as whole-genome duplication (WGD)-retained paralogs (~35% of all specific TFs). Using high-resolution profiling, we find that ~60% of paralogs evolved differential binding preferences. We show that this divergence results primarily from variations outside the DNA-binding domains (DBDs), while DBD preferences remain largely conserved. Analysis of non-WGD orthologs revealed uneven splitting of ancestral preferences between duplicates, and the preferential acquiring of new targets by the least conserved paralog (biased neo/sub-functionalization). Interactions between paralogs were rare, and, when present, occurred through weak competition for DNA-binding or dependency between dimer-forming paralogs. We discuss the implications of our findings for the evolutionary design of transcriptional networks.

Plant cytoplasmic ribosomal proteins: an update on classification, nomenclature, evolution and resources

Standardized naming systems are essential to integrate and unify distinct research fields, and to link multi-species data within and across kingdoms. We conducted a comprehensive survey of cytoplasmic ribosomal proteins (CRPs) in the dicot model Arabidopsis thaliana and the monocot model rice, noting that the standardized naming system has not been widely adopted in the plant community. We generated a database linking the old classical names to their updated and compliant names. We also explored the sequences, molecular evolution, and structural and functional characteristics of all plant CRP families, emphasizing evolutionarily conserved and plant-specific features through cross-kingdom comparisons. Unlike fungal CRP paralogs that were mainly created by whole-genome duplication (WGD) or retroposition under a concerted evolution mode, plant CRP genes evolved primarily through both WGD and tandem duplications in a rapid birth-and-death process. We also provide a web-based resource (http://www.plantcrp.cn/) with the aim of sharing the latest knowledge on plant CRPs and facilitating the continued development of a standardized framework across the entire community.

Genome-wide identification and expression analysis of HAK/KUP/KT potassium transporter provides insights into genes involved in responding to potassium deficiency and salt stress in pepper (Capsicum annuum L.)

In plants, the HAK/KUP/KT family is the largest group of potassium transporters, and it plays an important role in mineral element absorption, plant growth, environmental stress adaptation, and symbiosis. Although these important genes have been investigated in many plant species, limited information is currently available on the HAK/KUP/KT genes for Pepper (Capsicum annuum L.). In the present study, a total of 20 CaHAK genes were identified from the pepper genome and the CaHAK genes were numbered 1 - 20 based on phylogenetic analysis. For the genes and their corresponding proteins, the physicochemical properties, phylogenetic relationship, chromosomal distribution, gene structure, conserved motifs, gene duplication events, and expression patterns were analyzed. Phylogenetic analysis divided CaHAK genes into four cluster (I-IV) based on their structural features and the topology of the phylogenetic tree. Purifying selection played a crucial role in the evolution of CaHAK genes, while whole-genome triplication contributed to the expansion of the CaHAK gene family. The expression patterns showed that CaHAK proteins exhibited functional divergence in terms of plant K⁺ uptake and salt stress response. In particular, transcript abundance of CaHAK3 and CaHAK7 was strongly and specifically up-regulated in pepper roots under low K⁺ or high salinity conditions, suggesting that these genes are candidates for high-affinity K⁺ uptake transporters and may play crucial roles in the maintenance of the Na⁺/K⁺ balance during salt stress in pepper. In summary, the results not only provided the important information on the characteristics and evolutionary relationships of CaHAKs, but also provided potential genes responding to potassium deficiency and salt stress.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s13205-022-03136-z.

Machine learning approaches to identify core and dispensable genes in pangenomes

A gene in a given taxonomic group is either present in every individual (core) or absent in at least a single individual (dispensable). Previous pangenomic studies have identified certain functional differences between core and dispensable genes. However, identifying if a gene belongs to the core or dispensable portion of the genome requires the construction of a pangenome, which involves sequencing the genomes of many individuals. Here we aim to leverage the previously characterized core and dispensable gene content for two grass species [Brachypodium distachyon (L.) P. Beauv. and Oryza sativa L.] to construct a machine learning model capable of accurately classifying genes as core or dispensable using only a single annotated reference genome. Such a model may mitigate the need for pangenome construction, an expensive hurdle especially in orphan crops, which often lack the adequate genomic resources.

All Ways Lead to Rome—Meiotic Stabilization Can Take Many Routes in Nascent Polyploid Plants

Newly formed polyploids often show extensive meiotic defects, resulting in aneuploid gametes, and thus reduced fertility. However, while many neopolyploids are meiotically unstable, polyploid lineages that survive in nature are generally stable and fertile; thus, those lineages that survive are those that are able to overcome these challenges. Several genes that promote polyploid stabilization are now known in plants, allowing speculation on the evolutionary origin of these meiotic adjustments. Here, I discuss results that show that meiotic stability can be achieved through the differentiation of certain alleles of certain genes between ploidies. These alleles, at least sometimes, seem to arise by novel mutation, while standing variation in either ancestral diploids or related polyploids, from which alleles can introgress, may also contribute. Growing evidence also suggests that the coevolution of multiple interacting genes has contributed to polyploid stabilization, and in allopolyploids, the return of duplicated genes to single copies (genome fractionation) may also play a role in meiotic stabilization. There is also some evidence that epigenetic regulation may be important, which can help explain why some polyploid lineages can partly stabilize quite rapidly.

Ploidy level enhances the photosynthetic capacity of a tetraploid variety of Acer buergerianum Miq

Background

Polyploidy plays an important role in plant breeding and has widespread effects on photosynthetic capacity. To determine the photosynthetic capacity of the tetraploid variety Acer buergerianum Miq. ‘Xingwang’, we compared the gas exchange parameters, chloroplast structure, chlorophyll contents, and chlorophyll fluorescence parameters between the tetraploid Acer buergerianum ‘Xingwang’ and the diploid ‘S4’. To evaluate the effects of genome duplication on the photosynthetic capacity of Acer buergerianum ‘Xingwang’, the transcriptomes of the autotetraploid ‘Xingwang’ and the diploid ‘S4’ of A. buergerianum were compared.

Methods

The ploidy of Acer buergerianum ‘Xingwang’ was identified by flow cytometry and the chromosome counting method. An LI-6800 portable photosynthesis system analyzer was used to assess the gas exchange parameters of the tetraploid variety ‘Xingwang’ and diploid variety ‘S4’ of A. buergerianum. We used a BioMate 3S ultraviolet-visible spectrophotometer and portable modulated fluorometer to measure the chlorophyll contents and chlorophyll fluorescence parameters, respectively, of ‘Xingwang’ and ‘S4’. Illumina high-throughput sequencing technology was used to identify the differences in the genes involved in the photosynthetic differences and determine their expression characteristics.

Results

The single-cell DNA content and chromosome number of the tetraploid ‘Xingwang’ were twice those found in the normal diploid ‘S4’. In terms of gas exchange parameters, the change in stomatal conductance, change in intercellular CO₂ concentration, transpiration rate and net photosynthetic rate of ‘Xingwang’ were higher than those of the diploid ‘S4’. The chlorophyll contents, the maximal photochemical efficiency of PSII and the potential photochemical efficiency of PSII in ‘Xingwang’ were higher than those of ‘S4’. The chloroplasts of ‘Xingwang’ contained thicker thylakoid lamellae. By the use of Illumina sequencing technology, a total of 51,807 unigenes were obtained; they had an average length of 1,487 nt, and the average N50 was 2,034 nt. The lengths of most of the unigenes obtained ranged from 200–300 bp, with an average value of 5,262, followed by those longer than 3,000 bp, with an average value of 4,791. The data revealed numerous differences in gene expression between the two transcriptomes. In total, 24,221 differentially expressed genes were screened, and the percentage of differentially expressed genes was as high as 46.75% (24,224/51,807), of which 10,474 genes were upregulated and 13,747 genes were downregulated. We analyzed the key genes in the photosynthesis pathway and the porphyrin and chlorophyll metabolism pathway; the upregulation of HemB may promote an increase in the chlorophyll contents of ‘Xingwang’, and the upregulation of related genes in PSII and PSI may enhance the light harvesting of ‘Xingwang’, increasing its light energy conversion efficiency.

MicroRNAs as Indicators into the Causes and Consequences of Whole-Genome Duplication Events

Whole-genome duplications (WGDs) have long been considered the causal mechanism underlying dramatic increases to morphological complexity due to the neo-functionalization of paralogs generated during these events. Nonetheless, an alternative hypothesis suggests that behind the retention of most paralogs is not neo-functionalization, but instead the degree of the inter-connectivity of the intended gene product, as well as the mode of the WGD itself. Here, we explore both the causes and consequences of WGD by examining the distribution, expression, and molecular evolution of microRNAs (miRNAs) in both gnathostome vertebrates as well as chelicerate arthropods. We find that although the number of miRNA paralogs tracks the number of WGDs experienced within the lineage, few of these paralogs experienced changes to the seed sequence, and thus are functionally equivalent relative to their mRNA targets. Nonetheless, in gnathostomes, although the retention of paralogs following the 1R autotetraploidization event is similar across the two subgenomes, the paralogs generated by the gnathostome 2R allotetraploidization event are retained in higher numbers on one subgenome relative to the second, with the miRNAs found on the preferred subgenome showing both higher expression of mature miRNA transcripts and slower molecular evolution of the precursor miRNA sequences. Importantly, WGDs do not result in the creation of miRNA novelty, nor do WGDs correlate to increases in complexity. Instead, it is the number of miRNA seed sequences in the genome itself that not only better correlate to instances in complexification, but also mechanistically explain why complexity increases when new miRNA families are established.

Genetic Divergence of Lineage-Specific Tandemly Duplicated Gene Clusters in Four Diploid Potato Genotypes

Potato (Solanum tuberosum L.) is the most important non-grain food crop. Tandem duplication significantly contributes to genome evolution. The objectives of this study were to (i) identify tandemly duplicated genes and compare their genomic distributions across potato genotypes, (ii) investigate the bias in functional specificities, (iii) explore the relationships among coding sequence, promoter and expression divergences associated with tandemly duplicated genes, (iv) examine the role of tandem duplication in generating and expanding lineage-specific gene families, (v) investigate the evolutionary forces affecting tandemly duplicated genes, and (vi) assess the similarities and differences with respect to above mentioned aspects between cultivated genotypes and their wild-relative. In this study, we used well-annotated and chromosome-scale de novo genome assemblies of multiple potato genotypes. Our results showed that tandemly duplicated genes are abundant and dispersed through the genome. We found that several functional specificities, such as disease resistance, stress-tolerance, and biosynthetic pathways of tandemly duplicated genes were differentially enriched across multiple potato genomes. Our results indicated the existence of a significant correlation among expression, promoter, and protein divergences in tandemly duplicated genes. We found about one fourth of tandemly duplicated gene clusters as lineage-specific among multiple potato genomes, and these tended to localize toward centromeres and revealed distinct selection signatures and expression patterns. Furthermore, our results showed that a majority of duplicated genes were retained through sub-functionalization followed by genetic redundancy, while only a small fraction of duplicated genes was retained though neo-functionalization. The lineage-specific expansion of gene families by tandem duplication coupled with functional bias might have significantly contributed to potato's genotypic diversity, and, thus, to adaption to environmental stimuli.

Evolutionary implications of the RNA N6-methyladenosine methylome in plants

Epigenetic modifications play important roles in genome evolution and innovation. However, most analyses have focused on the evolutionary role of DNA modifications, and little is understood about the influence of post-transcriptional RNA modifications on genome evolution. To explore the evolutionary significance of RNA modifications, we generated transcriptome-wide profiles of N6-methyladenosine (m6A), the most prevalent internal modification of mRNA, for 13 representative plant species spanning over half a billion years of evolution. These data reveal the evolutionary conservation and divergence of m6A methylomes in plants, uncover the preference of m6A modifications on ancient orthologous genes, and demonstrate less m6A divergence between orthologous gene pairs with earlier evolutionary origins. Further investigation revealed that the evolutionary divergence of m6A modifications is related to sequence variation between homologs from whole genome duplication and gene family expansion from local genome duplication. Unexpectedly, a significant negative correlation was found between the retention ratio of m6A modifications and the number of family members. Moreover, the divergence of m6A modifications is accompanied by variation in the expression level and translation efficiency of duplicated genes from whole and local genome duplication. Our work reveals new insights into evolutionary patterns of m6A methylomes in plant species and their implications, and provides a resource of plant m6A profiles for further studies of m6A regulation and function in an evolutionary context.

The genome of Shorea leprosula (Dipterocarpaceae) highlights the ecological relevance of drought in aseasonal tropical rainforests

Hyperdiverse tropical rainforests, such as the aseasonal forests in Southeast Asia, are supported by high annual rainfall. Its canopy is dominated by the species-rich tree family of Dipterocarpaceae (Asian dipterocarps), which has both ecological (e.g., supports flora and fauna) and economical (e.g., timber production) importance. Recent ecological studies suggested that rare irregular drought events may be an environmental stress and signal for the tropical trees. We assembled the genome of a widespread but near threatened dipterocarp, Shorea leprosula, and analyzed the transcriptome sequences of ten dipterocarp species representing seven genera. Comparative genomic and molecular dating analyses suggested a whole-genome duplication close to the Cretaceous-Paleogene extinction event followed by the diversification of major dipterocarp lineages (i.e. Dipterocarpoideae). Interestingly, the retained duplicated genes were enriched for genes upregulated by no-irrigation treatment. These findings provide molecular support for the relevance of drought for tropical trees despite the lack of an annual dry season.

Trajectories of Homoeolog-Specific Expression in Allotetraploid Tragopogon castellanus Populations of Independent Origins

Polyploidization can have a significant ecological and evolutionary impact by providing substantially more genetic material that may result in novel phenotypes upon which selection may act. While the effects of polyploidization are broadly reviewed across the plant tree of life, the reproducibility of these effects within naturally occurring, independently formed polyploids is poorly characterized. The flowering plant genus Tragopogon (Asteraceae) offers a rare glimpse into the intricacies of repeated allopolyploid formation with both nascent (< 90 years old) and more ancient (mesopolyploids) formations. Neo- and mesopolyploids in Tragopogon have formed repeatedly and have extant diploid progenitors that facilitate the comparison of genome evolution after polyploidization across a broad span of evolutionary time. Here, we examine four independently formed lineages of the mesopolyploid Tragopogon castellanus for homoeolog expression changes and fractionation after polyploidization. We show that expression changes are remarkably similar among these independently formed polyploid populations with large convergence among expressed loci, moderate convergence among loci lost, and stochastic silencing. We further compare and contrast these results for T. castellanus with two nascent Tragopogon allopolyploids. While homoeolog expression bias was balanced in both nascent polyploids and T. castellanus, the degree of additive expression was significantly different, with the mesopolyploid populations demonstrating more non-additive expression. We suggest that gene dosage and expression noise minimization may play a prominent role in regulating gene expression patterns immediately after allopolyploidization as well as deeper into time, and these patterns are conserved across independent polyploid lineages.

HAG1 and SWI3A/B control of male germ line development in P. patens suggests conservation of epigenetic reproductive control across land plants

KEY MESSAGE

Bryophytes as models to study the male germ line: loss-of-function mutants of epigenetic regulators HAG1 and SWI3a/b demonstrate conserved function in sexual reproduction. With the water-to-land transition, land plants evolved a peculiar haplodiplontic life cycle in which both the haploid gametophyte and the diploid sporophyte are multicellular. The switch between these phases was coined alternation of generations. Several key regulators that control the bauplan of either generation are already known. Analyses of such regulators in flowering plants are difficult due to the highly reduced gametophytic generation, and the fact that loss of function of such genes often is embryo lethal in homozygous plants. Here we set out to determine gene function and conservation via studies in bryophytes. Bryophytes are sister to vascular plants and hence allow evolutionary inferences. Moreover, embryo lethal mutants can be grown and vegetatively propagated due to the dominance of the bryophyte gametophytic generation. We determined candidates by selecting single copy orthologs that are involved in transcriptional control, and of which flowering plant mutants show defects during sexual reproduction, with a focus on the under-studied male germ line. We selected two orthologs, SWI3a/b and HAG1, and analyzed loss-of-function mutants in the moss P. patens. In both mutants, due to lack of fertile spermatozoids, fertilization and hence the switch to the diploid generation do not occur. Pphag1 additionally shows arrested male and impaired female gametangia development. We analyzed HAG1 in the dioecious liverwort M. polymorpha and found that in Mphag1 the development of gametangiophores is impaired. Taken together, we find that involvement of both regulators in sexual reproduction is conserved since the earliest divergence of land plants.

Genome-Wide Analysis of Coding and Non-coding RNA Reveals a Conserved miR164-NAC-mRNA Regulatory Pathway for Disease Defense in Populus

MicroRNAs (miRNAs) contribute to plant defense responses by increasing the overall genetic diversity; however, their origins and functional importance in plant defense remain unclear. Here, we employed Illumina sequencing technology to assess how miRNA and messenger RNA (mRNA) populations vary in the Chinese white poplar (Populus tomentosa) during a leaf black spot fungus (Marssonina brunnea) infection. We sampled RNAs from infective leaves at conidia germinated stage [12 h post-inoculation (hpi)], infective vesicles stage (24 hpi), and intercellular infective hyphae stage (48 hpi), three essential stages associated with plant colonization and biotrophic growth in M. brunnea fungi. In total, 8,938 conserved miRNA-target gene pairs and 3,901 Populus-specific miRNA-target gene pairs were detected. The result showed that Populus-specific miRNAs (66%) were more involved in the regulation of the disease resistance genes. By contrast, conserved miRNAs (>80%) target more whole-genome duplication (WGD)-derived transcription factors (TFs). Among the 1,023 WGD-derived TF pairs, 44.9% TF pairs had only one paralog being targeted by a miRNA that could be due to either gain or loss of a miRNA binding site after the WGD. A conserved hierarchical regulatory network combining promoter analyses and hierarchical clustering approach uncovered a miR164–NAM, ATAF, and CUC (NAC) transcription factor–mRNA regulatory module that has potential in Marssonina defense responses. Furthermore, analyses of the locations of miRNA precursor sequences reveal that pseudogenes and transposon contributed a certain proportion (∼30%) of the miRNA origin. Together, these observations provide evolutionary insights into the origin and potential roles of miRNAs in plant defense and functional innovation.

Insights of Phenolic Pathway in Fruits: Transcriptional and Metabolic Profiling in Apricot (Prunus armeniaca)

There is an increasing interest in polyphenols, plant secondary metabolites, in terms of fruit quality and diet, mainly due to their antioxidant effect. However, the identification of key gene enzymes and their roles in the phenylpropanoid pathway in temperate fruits species remains uncertain. Apricot (Prunus armeniaca) is a Mediterranean fruit with high diversity and fruit quality properties, being an excellent source of polyphenol compounds. For a better understanding of the phenolic pathway in these fruits, we selected a set of accessions with genetic-based differences in phenolic compounds accumulation. HPLC analysis of the main phenolic compounds and transcriptional analysis of the genes involved in key steps of the polyphenol network were carried out. Phenylalanine ammonia-lyase (PAL), dihydroflavonol-4-reductase (DFR) and flavonol synthase (FLS) were the key enzymes selected. Orthologous of the genes involved in transcription of these enzymes were identified in apricot: ParPAL1, ParPAL2, ParDFR, ParFLS1 and ParFLS2. Transcriptional data of the genes involved in those critical points and their relationships with the polyphenol compounds were analyzed. Higher expression of ParDFR and ParPAL2 has been associated with red-blushed accessions. Differences in expression between paralogues could be related to the presence of a BOXCOREDCPAL cis-acting element related to the genes involved in anthocyanin synthesis ParFLS2, ParDFR and ParPAL2.

Learning Retention Mechanisms and Evolutionary Parameters of Duplicate Genes from Their Expression Data

Learning about the roles that duplicate genes play in the origins of novel phenotypes requires an understanding of how their functions evolve. A previous method for achieving this goal, CDROM, employs gene expression distances as proxies for functional divergence and then classifies the evolutionary mechanisms retaining duplicate genes from comparisons of these distances in a decision tree framework. However, CDROM does not account for stochastic shifts in gene expression or leverage advances in contemporary statistical learning for performing classification, nor is it capable of predicting the parameters driving duplicate gene evolution. Thus, here we develop CLOUD, a multi-layer neural network built on a model of gene expression evolution that can both classify duplicate gene retention mechanisms and predict their underlying evolutionary parameters. We show that not only is the CLOUD classifier substantially more powerful and accurate than CDROM, but that it also yields accurate parameter predictions, enabling a better understanding of the specific forces driving the evolution and long-term retention of duplicate genes. Further, application of the CLOUD classifier and predictor to empirical data from Drosophila recapitulates many previous findings about gene duplication in this lineage, showing that new functions often emerge rapidly and asymmetrically in younger duplicate gene copies, and that functional divergence is driven by strong natural selection. Hence, CLOUD represents a major advancement in classifying retention mechanisms and predicting evolutionary parameters of duplicate genes, thereby highlighting the utility of incorporating sophisticated statistical learning techniques to address long-standing questions about evolution after gene duplication.

Evolution of Conserved Noncoding Sequences in Arabidopsis thaliana

Recent pangenome studies have revealed a large fraction of the gene content within a species exhibits presence-absence variation (PAV). However, coding regions alone provide an incomplete assessment of functional genomic sequence variation at the species level. Little to no attention has been paid to noncoding regulatory regions in pangenome studies, though these sequences directly modulate gene expression and phenotype. To uncover regulatory genetic variation, we generated chromosome-scale genome assemblies for thirty Arabidopsis thaliana accessions from multiple distinct habitats and characterized species level variation in Conserved Noncoding Sequences (CNS). Our analyses uncovered not only PAV and positional variation (PosV) but that diversity in CNS is nonrandom, with variants shared across different accessions. Using evolutionary analyses and chromatin accessibility data, we provide further evidence supporting roles for conserved and variable CNS in gene regulation. Additionally, our data suggests that transposable elements contribute to CNS variation. Characterizing species-level diversity in all functional genomic sequences may later uncover previously unknown mechanistic links between genotype and phenotype.

The Effects of Gene Duplication Modes on the Evolution of Regulatory Divergence in Wild and Cultivated Soybean

Regulatory changes include divergence in both cis-elements and trans-factors, which play roles in organismal evolution. Whole genome duplications (WGD) followed by diploidization are a recurrent feature in the evolutionary history of angiosperms. Prior studies have shown that duplicated genes have different evolutionary fates due to variable selection constraints and results in genomic compositions with hallmarks of paleopolyploidy. The recent sequential WGDs and post-WGD evolution in the common ancestor of cultivated soybean (Glycine max) and wild soybean (Glycine soja), together with other models of gene duplication, have resulted in a highly duplicated genome. In this study, we investigated the transcriptional changes in G. soja and G. max. We identified a sizable proportion of interspecific differentially expressed genes (DEGs) and found parental expression level dominance of G. max in their F1 hybrids. By classifying genes into different regulatory divergence types, we found the trans-regulatory changes played a predominant role in transcriptional divergence between wild and cultivated soybean. The same gene ontology (GO) and protein family (Pfam) terms were found to be over-represented in DEGs and genes of cis-only between JY47 and GS, suggesting the substantial contribution of cis-regulatory divergences to the evolution of wild and cultivated soybeans. By further dissecting genes into five different duplication modes, we found genes in different duplication modes tend to accumulate different types of regulatory differences. A relatively higher proportion of cis-only regulatory divergences was detected in singleton, dispersed, proximal, and tandem duplicates than WGD duplicates and genome-wide level, which is in line with the prediction of gene balance hypothesis for the differential fates of duplicated genes post-WGD. The numbers of cis-only and trans-only regulated genes were similar for singletons, whereas there were more genes of trans-only than cis-only in the rest duplication types, especially in WGD in which there were two times more trans-only genes than that in cis-only type. Tandem duplicates showed the highest proportion of trans-only genes probably due to some special features of this class. In summary, our results demonstrate that genes in different duplication modes have different fates in transcriptional evolution underpinned by cis- or trans-regulatory divergences in soybean and likely in other paleopolyploid higher organisms.

Identification and Characterization of PLATZ Transcription Factors in Wheat

The PLATZ (plant AT-rich protein and zinc-binding protein) transcription factor family is a class of plant-specific zinc-dependent DNA-binding proteins. PLATZ has essential roles in seed endosperm development, as well as promoting cell proliferation duration in the earlier stages of the crops. In the present study, 62 TaPLATZ genes were identified from the wheat genome, and they were unequally distributed on 15 chromosomes. According to the phylogenetic analysis, 62 TaPLATZ genes were classified into six groups, including two groups that were unique in wheat. Members in the same groups shared similar exon-intron structures. The polyploidization, together with genome duplication of wheat, plays a crucial role in the expansion of the TaPLATZs family. Transcriptome data indicated a distinct divergence expression pattern of TaPLATZ genes that could be clustered into four modules. The TaPLATZs in Module b possessed a seed-specific expression pattern and displayed obvious high expression in the earlier development stage of seeds. Subcellular localization data of TaPLATZs suggesting that they likely perform a function as a conventional transcription factor. This study provides insight into understanding the structure divergence, evolutionary features, expression profiles, and potential function of PLATZ in wheat.

Evolution of Homeologous Gene Expression in Polyploid Wheat

Polyploidization has played a prominent role in the evolutionary history of plants. Two recent and sequential allopolyploidization events have resulted in the formation of wheat species with different ploidies, and which provide a model to study the effects of polyploidization on the evolution of gene expression. In this study, we identified differentially expressed genes (DEGs) between four BBAA tetraploid wheats of three different ploidy backgrounds. DEGs were found to be unevenly distributed among functional categories and duplication modes. We observed more DEGs in the extracted tetraploid wheat (ETW) than in natural tetraploid wheats (TD and TTR13) as compared to a synthetic tetraploid (AT2). Furthermore, DEGs showed higher Ka/Ks ratios than those that did not show expression changes (non-DEGs) between genotypes, indicating DEGs and non-DEGs experienced different selection pressures. For A-B homeolog pairs with DEGs, most of them had only one differentially expressed copy, however, when both copies of a homeolog pair were DEGs, the A and B copies were more likely to be regulated to the same direction. Our results suggest that both cis- and inter-subgenome trans-regulatory changes are important drivers in the evolution of homeologous gene expression in polyploid wheat, with ploidy playing a significant role in the process.

Molecular regulatory mechanisms underlying the adaptability of polyploid plants

Polyploidization influences the genetic composition and gene expression of an organism. This multi-level genetic change allows the formation of new regulatory pathways leading to increased adaptability. Although both forms of polyploidization provide advantages, autopolyploids were long thought to have little impact on plant divergence compared to allopolyploids due to their formation through genome duplication only, rather than in combination with hybridization. Recent advances have begun to clarify the molecular regulatory mechanisms such as microRNAs, alternative splicing, RNA-binding proteins, histone modifications, chromatin remodelling, DNA methylation, and N⁶ -methyladenosine (m6A) RNA methylation underlying the evolutionary success of polyploids. Such research is expanding our understanding of the evolutionary adaptability of polyploids and the regulatory pathways that allow adaptive plasticity in a variety of plant species. Herein we review the roles of individual molecular regulatory mechanisms and their potential synergistic pathways underlying plant evolution and adaptation. Notably, increasing interest in m6A methylation has provided a new component in potential mechanistic coordination that is still predominantly unexplored. Future research should attempt to identify and functionally characterize the evolutionary impact of both individual and synergistic pathways in polyploid plant species.

Evidence for sub-functionalization of tandemly duplicated XPB nucleotide excision repair genes in Arabidopsis thaliana

Plants are continuously exposed to agents that can generate DNA lesions. Nucleotide Excision Repair (NER) is one of the repair pathways employed by plants to protect their genome, including from sunlight. The Xeroderma Pigmentosum type B (XPB) protein is a DNA helicase shown to be involved in NER and is also an essential subunitof the Transcription Factor IIH (TFIIH) complex. XPB was found to be a single copy gene in eukaryotes, but found as a tandem duplication in the plant Arabidopsis thaliana, AtXPB1 and AtXPB2. We aimed to investigate whether the XPB in tandem duplication was common within members of the Brassicaceae. We analyzed genomic DNA of species from different tribes of the family and the results indicate that the tandem duplication occurred in Camelineae tribe ancestor, of which A. thaliana belongs, at approximately 8 million years ago. Further experiments were devised to study possible functional roles for the A. thaliana AtXPB paralogs. A non-coincident expression profile of the paralogs was observed in various plant organs, developmental and cell cycle stages. AtXPB2 expression was observed in proliferating cells and clustered with the transcription of other components of the TFIIH such as p44, p52 and XPD/UVH6 along the cell cycle. AtXPB1 gene transcription, on the other hand, was enhanced specifically after UV-B irradiation in leaf trichomes. Altogether, our results reported herein suggest a functional specialization for the AtXPB paralogs: while the AtXPB2 paralog may have a role in cell proliferation and repair as XPB of other eukaryotes, the AtXPB1 paralog is most likely implicated in repair functions in highly specialized A. thaliana cells.