Evolution and post-transcriptional regulation insights of m6A writers, erasers, and readers in plant epitranscriptome.

As a dynamic and reversible post-transcriptional marker, N6-methyladenosine (m6A) plays an important role in the regulation of biological functions, which are mediated by m6A pathway components including writers (MT-A70, FIP37, VIR and HAKAI family), erasers (ALKBH family) and readers (YTH family). There is an urgent need for a comprehensive analysis of m6A pathway components across species at evolutionary levels. In this study, we identified 4062 m6A pathway components from 154 plant species including green algae, utilizing large-scale phylogenetic to explore their origin and evolution. We discovered that the copy number of writers was conserved among different plant lineages, with notable expansions in the ALKBH and YTH families. Synteny network analysis revealed conserved genomic contexts and lineage-specific transpositions. Furthermore, we used Direct RNA Sequencing (DRS) to reveal the Poly(A) length (PAL) and m6A ratio profiles in six angiosperms species, with a particular focus on the m6A pathway components. The ECT1/2-Poeaece4 sub-branches (YTH family) with unique genomic contexts exhibited significantly higher expression level than genes of other ECT1/2 poeaece sub-branches (ECT1/2-Poeaece1-3), accompanied by lower m6A modification and PAL. Besides, conserved m6A sites distributed in CDS and 3'UTR were detected in the ECT1/2-Poaceae4, and the dual-luciferase assay further demonstrated that these conserved m6A sites in the 3'UTR negatively regulated the expression of Firefly luciferase (LUC) gene. Finally, we developed transcription factor regulatory networks for m6A pathway components, using yeast one-hybrid assay demonstrated that PheBPC1 could interact with the PheECT1/2-5 promoter. Overall, this study presents a comprehensive evolutionary and functional analysis of m6A pathway components and their modifications in plants, providing a valuable resource for future functional analysis in this field.

Ancient Duplication and Lineage-Specific Transposition Determine Evolutionary Trajectory of ERF Subfamily across Angiosperms.

AP2/ERF transcription factor family plays an important role in plant development and stress responses. Previous studies have shed light on the evolutionary trajectory of the AP2 and DREB subfamilies. However, knowledge about the evolutionary history of the ERF subfamily in angiosperms still remains limited. In this study, we performed a comprehensive analysis of the ERF subfamily from 107 representative angiosperm species by combining phylogenomic and synteny network approaches. We observed that the expansion of the ERF subfamily was driven not only by whole-genome duplication (WGD) but also by tandem duplication (TD) and transposition duplication events. We also found multiple transposition events in Poaceae, Brassicaceae, Poales, Brassicales, and Commelinids. These events may have had notable impacts on copy number variation and subsequent functional divergence of the ERF subfamily. Moreover, we observed a number of ancient tandem duplications occurred in the ERF subfamily across angiosperms, e.g., in Subgroup IX, IXb originated from ancient tandem duplication events within IXa. These findings together provide novel insights into the evolution of this important transcription factor family.

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

BACKGROUND: In the evolutionary study of gene families, exploring the duplication mechanisms of gene families helps researchers understand their evolutionary history. The tubby-like protein (TLP) family is essential for growth and development in plants and animals. Much research has been done on its function; however, limited information is available with regard to the evolution of the TLP gene family. Herein, we systematically investigated the evolution of TLP genes in seven representative Poaceae lineages. RESULTS: Our research showed that the evolution of TLP genes was influenced not only by whole-genome duplication (WGD) and dispersed duplication (DSD) but also by transposed duplication (TRD), which has been neglected in previous research. For TLP family size, we found an evolutionary pattern of progressive shrinking in the grass family. Furthermore, the evolution of the TLP gene family was at least affected by evolutionary driving forces such as duplication, purifying selection, and base mutations. CONCLUSIONS: This study presents the first comprehensive evolutionary analysis of the TLP gene family in grasses. We demonstrated that the TLP gene family is also influenced by a transposed duplication mechanism. Several new insights into the evolution of the TLP gene family are presented. This work provides a good reference for studying gene evolution and the origin of duplication.

New insights into the phylogeny of the TMBIM superfamily across the tree of life: Comparative genomics and synteny networks reveal independent evolution of the BI and LFG families in plants.

The Transmembrane BAX Inhibitor Motif containing (TMBIM) superfamily, divided into BAX Inhibitor (BI) and Lifeguard (LFG) families, comprises a group of cytoprotective cell death regulators conserved in prokaryotes and eukaryotes. However, no research has focused on the evolution of this superfamily in plants. We identified 685 TMBIM proteins in 171 organisms from Archaea, Bacteria, and Eukarya, and provided a phylogenetic overview of the whole TMBIM superfamily. Then, we used orthology and synteny network analyses to further investigate the evolution and expansion of the BI and LFG families in 48 plants from diverse taxa. Plant BI family forms a single monophyletic group; however, monocot BI sequences transposed to another genomic context during evolution. Plant LFG family, which expanded trough whole genome and tandem duplications, is subdivided in LFG I, LFG IIA, and LFG IIB major phylogenetic groups, and retains synteny in angiosperms. Moreover, two orthologous groups (OGs) are shared between bryophytes and seed plants. Other several lineage-specific OGs are present in plants. This work clarifies the phylogenetic classification of the TMBIM superfamily across the three domains of life. Furthermore, it sheds new light on the evolution of the BI and LFG families in plants providing a benchmark for future research.

Exploring the evolution of CHS gene family in plants.

Chalcone synthase (CHS) is a key enzyme that catalyzes the first committed step of flavonoid biosynthetic pathway. It plays a vital role not only in maintaining plant growth and development, but also in regulating plant response to environmental hazards. However, the systematic phylogenomic analysis of CHS gene family in a wide range of plant species has not been reported yet. To fill this knowledge gap, a large-scale investigation of CHS genes was performed in 178 plant species covering green algae to dicotyledons. A total of 2,011 CHS and 293 CHS-like genes were identified and phylogenetically divided into four groups, respectively. Gene distribution patterns across the plant kingdom revealed the origin of CHS can be traced back to before the rise of algae. The gene length varied largely in different species, while the exon structure was relatively conserved. Selection pressure analysis also indicated the conserved features of CHS genes on evolutionary time scales. Moreover, our synteny analysis pinpointed that, besides genome-wide duplication and tandem duplication, lineage specific transposition events also occurred in the evolutionary trajectory of CHS gene family. This work provides novel insights into the evolution of CHS gene family and may facilitate further research to better understand the regulatory mechanism of traits relating to flavonoid biosynthesis in diverse plants.

Genomic and Transcriptomic Insights into the Evolution of C4 Photosynthesis in Grasses.

C4 photosynthesis has independently evolved over 62 times within 19 angiosperm families. The recurrent evolution of C4 photosynthesis appears to contradict the complex anatomical and biochemical modifications required for the transition from C3 to C4 photosynthesis. In this study, we conducted an integrated analysis of genomics and transcriptomics to elucidate the molecular underpinnings of convergent C4 evolution in the grass family. Our genome-wide exploration of C4-related gene families suggests that the expansion of these gene families may have played an important role in facilitating C4 evolution in the grass family. A phylogenomic synteny network analysis uncovered the emergence of C4 genes in various C4 grass lineages from a common ancestral gene pool. Moreover, through a comparison between non-C4 and C4  PEPCs, we pinpointed 14 amino acid sites exhibiting parallel adaptations. These adaptations, occurring post the BEP-PACMAD divergence, shed light on why all C4 origins in grasses are confined to the PACMAD clade. Furthermore, our study revealed that the ancestor of Chloridoideae grasses possessed a more favorable molecular preadaptation for C4 functions compared to the ancestor of Panicoideae grasses. This molecular preadaptation potentially explains why C4 photosynthesis evolved earlier in Chloridoideae than in Panicoideae and why the C3-to-C4 transition occurred once in Chloridoideae but multiple times in Panicoideae. Additionally, we found that C4 genes share similar cis-elements across independent C4 lineages. Notably, NAD-ME subtype grasses may have retained the ancestral regulatory machinery of the C4  NADP-ME gene, while NADP-ME subtype grasses might have undergone unique cis-element modifications.

Unraveling the evolutionary dynamics of the TPS gene family in land plants.

Terpenes and terpenoids are key natural compounds for plant defense, development, and composition of plant oil. The synthesis and accumulation of a myriad of volatile terpenoid compounds in these plants may dramatically alter the quality and flavor of the oils, which provide great commercial utilization value for oil-producing plants. Terpene synthases (TPSs) are important enzymes responsible for terpenic diversity. Investigating the differentiation of the TPS gene family could provide valuable theoretical support for the genetic improvement of oil-producing plants. While the origin and function of TPS genes have been extensively studied, the exact origin of the initial gene fusion event - it occurred in plants or microbes - remains uncertain. Furthermore, a comprehensive exploration of the TPS gene differentiation is still pending. Here, phylogenetic analysis revealed that the fusion of the TPS gene likely occurred in the ancestor of land plants, following the acquisition of individual C- and N- terminal domains. Potential mutual transfer of TPS genes was observed among microbes and plants. Gene synteny analysis disclosed a differential divergence pattern between TPS-c and TPS-e/f subfamilies involved in primary metabolism and those (TPS-a/b/d/g/h subfamilies) crucial for secondary metabolites. Biosynthetic gene clusters (BGCs) analysis suggested a correlation between lineage divergence and potential natural selection in structuring terpene diversities. This study provides fresh perspectives on the origin and evolution of the TPS gene family.

Applying Synteny Networks (SynNet) to Study Genomic Arrangements of Protein-Coding Genes in Plants.

In comparative genomics, the study of synteny can be a powerful method for exploring genome rearrangements, inferring genomic ancestry, defining orthology relationships, determining gene and genome duplications, and inferring gene positional conservation patterns across taxa. In this chapter, we present a step-by-step protocol for microsynteny network (SynNet) analysis, as an alternative to traditional methods of synteny comparison, where nodes in the network represent protein-coding genes and edges represent the pairwise syntenic relationships. The SynNet pipeline consists of six main steps: (1) pairwise genome comparisons between all the genomes being analyzed, (2) detection of inter- and intrasynteny blocks, (3) generation of an entire synteny database (i.e., edgelist), (4) network clustering, (5) phylogenomic profiling of the gene family of interest, and (6) evolutionary inference. The SynNet approach facilitates the rapid analysis and visualization of synteny relationships (from specific genes, specific gene families up to all genes) across a large number of genomes.

A Systematic Phylogenomic Classification of the Multidrug and Toxic Compound Extrusion Transporter Gene Family in Plants.

Multidrug and toxic compound extrusion (MATE) transporters comprise a multigene family that mediates multiple functions in plants through the efflux of diverse substrates including organic molecules, specialized metabolites, hormones, and xenobiotics. MATE classification based on genome-wide studies remains ambiguous, likely due to a lack of large-scale phylogenomic studies and/or reference sequence datasets. To resolve this, we established a phylogeny of the plant MATE gene family using a comprehensive kingdom-wide phylogenomic analysis of 74 diverse plant species. We identified more than 4,000 MATEs, which were classified into 14 subgroups based on a systematic bioinformatics pipeline using USEARCH, blast+ and synteny network tools. Our classification was performed using a four-step process, whereby MATEs sharing ≥ 60% protein sequence identity with a ≤ 1E-05 threshold at different sequence lengths (either full-length, ≥ 60% length, or ≥ 150 amino acids) or retaining in the similar synteny blocks were assigned to the same subgroup. In this way, we assigned subgroups to 95.8% of the identified MATEs, which we substantiated using synteny network clustering analysis. The subgroups were clustered under four major phylogenetic groups and named according to their clockwise appearance within each group. We then generated a reference sequence dataset, the usefulness of which was demonstrated in the classification of MATEs in additional species not included in the original analysis. Approximately 74% of the plant MATEs exhibited synteny relationships with angiosperm-wide or lineage-, order/family-, and species-specific conservation. Most subgroups evolved independently, and their distinct evolutionary trends were likely associated with the development of functional novelties or the maintenance of conserved functions. Together with the systematic classification and synteny network profiling analyses, we identified all the major evolutionary events experienced by the MATE gene family in plants. We believe that our findings and the reference dataset provide a valuable resource to guide future functional studies aiming to explore the key roles of MATEs in different aspects of plant physiology. Our classification framework can also be readily extendable to other (super) families.

Phylogenomic classification and synteny network analyses deciphered the evolutionary landscape of aldo-keto reductase (AKR) gene superfamily in the plant kingdom.

Aldo-keto reductase-domain (PF00248) containing proteins (AKRs) are NAD(P)(H)-dependent oxidoreductases of a multigene superfamily that mediate versatile functions in plants ranging from detoxification, metal chelation, potassium ion efflux to specialized metabolism. To uncover the complete repertoire of AKR gene superfamily in plants, a systematic kingdom-wide identification, phylogeny reconstruction, classification and synteny network clustering analyses were performed in this study using 74 diverse plant genomes. Plant AKRs were omnipresent, legitimately classified into 4 groups (based on phylogeny) and 14 subgroups (based on the ≥ 60% of protein sequence identity). Species composition of AKR subgroups highlights their distinct emergence during plant evolution. Loss of AKR subgroups among plants was apparent and that various lineage-, order/family- and species-specific losses were observed. The subgroups IA, IVB and IVF were flourished and diversified well during plant evolution, likely related to the complexity of plant's specialized metabolism and environmental adaptation. About 65% of AKRs were in genomic synteny regions across the plant kingdom and the AKRs relevant to important functions (e.g. vitamin B6 metabolism) were in profoundly conserved angiosperm-wide synteny communities. This study underscores the evolutionary landscape of plant AKRs and provides a comprehensive resource to facilitate the functional characterization of them.

Distinctive Features of PipX, a Unique Signaling Protein of Cyanobacteria.

PipX is a unique cyanobacterial protein identified by its ability to bind to PII and NtcA, two key regulators involved in the integration of signals of the nitrogen/carbon and energy status, with a tremendous impact on nitrogen assimilation and gene expression in cyanobacteria. PipX provides a mechanistic link between PII, the most widely distributed signaling protein, and NtcA, a global transcriptional regulator of cyanobacteria. PII, required for cell survival unless PipX is inactivated or down-regulated, functions by protein-protein interactions with transcriptional regulators, transporters, and enzymes. In addition, PipX appears to be involved in a wider signaling network, supported by the following observations: (i) PII-PipX complexes interact with PlmA, an as yet poorly characterized transcriptional regulator also restricted to cyanobacteria; (ii) the pipX gene is functionally connected with pipY, a gene encoding a universally conserved pyridoxal phosphate binding protein (PLPBP) involved in vitamin B6 and amino acid homeostasis, whose loss-of-function mutations cause B6-dependent epilepsy in humans, and (iii) pipX is part of a relatively robust, six-node synteny network that includes pipY and four additional genes that might also be functionally connected with pipX. In this overview, we propose that the study of the protein-protein interaction and synteny networks involving PipX would contribute to understanding the peculiarities and idiosyncrasy of signaling pathways that are conserved in cyanobacteria.

Dissecting the Genome-Wide Evolution and Function of R2R3-MYB Transcription Factor Family in Rosa chinensis.

Rosa chinensis, an important ancestor species of Rosa hybrida, the most popular ornamental plant species worldwide, produces flowers with diverse colors and fragrances. The R2R3-MYB transcription factor family controls a wide variety of plant-specific metabolic processes, especially phenylpropanoid metabolism. Despite their importance for the ornamental value of flowers, the evolution of R2R3-MYB genes in plants has not been comprehensively characterized. In this study, 121 predicted R2R3-MYB gene sequences were identified in the rose genome. Additionally, a phylogenomic synteny network (synnet) was applied for the R2R3-MYB gene families in 35 complete plant genomes. We also analyzed the R2R3-MYB genes regarding their genomic locations, Ka/Ks ratio, encoded conserved motifs, and spatiotemporal expression. Our results indicated that R2R3-MYBs have multiple synteny clusters. The RcMYB114a gene was included in the Rosaceae-specific Cluster 54, with independent evolutionary patterns. On the basis of these results and an analysis of RcMYB114a-overexpressing tobacco leaf samples, we predicted that RcMYB114a functions in the phenylpropanoid pathway. We clarified the relationship between R2R3-MYB gene evolution and function from a new perspective. Our study data may be relevant for elucidating the regulation of floral metabolism in roses at the transcript level.