Phylogenomic Synteny Network Analysis of MADS-Box Transcription Factor Genes Reveals Lineage-Specific Transpositions, Ancient Tandem Duplications, and Deep Positional Conservation

Vole genomics links determinate and indeterminate growth of teeth

Continuously growing teeth are an important innovation in mammalian evolution, yet genetic regulation of continuous growth by stem cells remains incompletely understood. Dental stem cells responsible for tooth crown growth are lost at the onset of tooth root formation. Genetic signaling that initiates this loss is difficult to study with the ever-growing incisor and rooted molars of mice, the most common mammalian dental model species, because signals for root formation overlap with signals that pattern tooth size and shape (i.e., cusp patterns). Different species of voles (Cricetidae, Rodentia, Glires) have evolved rooted and unrooted molars that have similar size and shape, providing alternative models for studying roots. We assembled a de novo genome of Myodes glareolus, a vole with high-crowned, rooted molars, and performed genomic and transcriptomic analyses in a broad phylogenetic context of Glires (rodents and lagomorphs) to assess differential selection and evolution in tooth forming genes. We identified 15 dental genes with changing synteny relationships and six dental genes undergoing positive selection across Glires, two of which were undergoing positive selection in species with unrooted molars, Dspp and Aqp1. Decreased expression of both genes in prairie voles with unrooted molars compared to bank voles supports the presence of positive selection and may underlie differences in root formation. Bulk transcriptomics analyses of embryonic molar development in bank voles also demonstrated conserved patterns of dental gene expression compared to mice, with species-specific variation likely related to developmental timing and morphological differences between mouse and vole molars. Our results support ongoing evolution of dental genes across Glires, revealing the complex evolutionary background of convergent evolution for ever-growing molars.

Transposable elements in Rosaceae: insights into genome evolution, expression dynamics, and syntenic gene regulation

Transposable elements (TEs) exert significant influence on plant genomic structure and gene expression. Here, we explored TE-related aspects across 14 Rosaceae genomes, investigating genomic distribution, transposition activity, expression patterns, and nearby differentially expressed genes (DEGs). Analyses unveiled distinct long terminal repeat retrotransposon (LTR-RT) evolutionary patterns, reflecting varied genome size changes among nine species over the past million years. In the past 2.5 million years, Rubus idaeus showed a transposition rate twice as fast as Fragaria vesca, while Pyrus bretschneideri displayed significantly faster transposition compared with Crataegus pinnatifida. Genes adjacent to recent TE insertions were linked to adversity resistance, while those near previous insertions were functionally enriched in morphogenesis, enzyme activity, and metabolic processes. Expression analysis revealed diverse responses of LTR-RTs to internal or external conditions. Furthermore, we identified 3695 pairs of syntenic DEGs proximal to TEs in Malus domestica cv. 'Gala' and M. domestica (GDDH13), suggesting TE insertions may contribute to varietal trait differences in these apple varieties. Our study across representative Rosaceae species underscores the pivotal role of TEs in plant genome evolution within this diverse family. It elucidates how these elements regulate syntenic DEGs on a genome-wide scale, offering insights into Rosaceae-specific genomic evolution.

Evolutionary origin and gradual accumulation with plant evolution of the LACS family

BACKGROUND

LACS (long-chain acyl-CoA synthetase) genes are widespread in organisms and have multiple functions in plants, especially in lipid metabolism. However, the origin and evolutionary dynamics of the LACS gene family remain largely unknown.

RESULTS

Here, we identified 1785 LACS genes in the genomes of 166 diverse plant species and identified the clades (I, II, III, IV, V, VI) of six clades for the LACS gene family of green plants through phylogenetic analysis. Based on the evolutionary history of plant lineages, we found differences in the origins of different clades, with Clade IV originating from chlorophytes and representing the origin of LACS genes in green plants. The structural characteristics of different clades indicate that clade IV is relatively independent, while the relationships between clades (I, II, III) and clades (V, VI) are closer. Dispersed duplication (DSD) and transposed duplication (TRD) are the main forces driving the evolution of plant LACS genes. Network clustering analysis further grouped all LACS genes into six main clusters, with genes within each cluster showing significant co-linearity. Ka/Ks results suggest that LACS family genes underwent purifying selection during evolution. We analyzed the phylogenetic relationships and characteristics of six clades of the LACS gene family to explain the origin, evolutionary history, and phylogenetic relationships of different clades and proposed a hypothetical evolutionary model for the LACS family of genes in plants.

CONCLUSIONS

Our research provides genome-wide insights into the evolutionary history of the LACS gene family in green plants. These insights lay an important foundation for comprehensive functional characterization in 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.

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.

Diversified molecular adaptations of inorganic nitrogen assimilation and signaling machineries in plants

Plants evolved sophisticated machineries to monitor levels of external nitrogen supply, respond to nitrogen demand from different tissues and integrate this information for coordinating its assimilation. Although roles of inorganic nitrogen in orchestrating developments have been studied in model plants and crops, systematic understanding of the origin and evolution of its assimilation and signaling machineries remains largely unknown. We expanded taxon samplings of algae and early-diverging land plants, covering all main lineages of Archaeplastida, and reconstructed the evolutionary history of core components involved in inorganic nitrogen assimilation and signaling. Most components associated with inorganic nitrogen assimilation were derived from the ancestral Archaeplastida. Improvements of assimilation machineries by gene duplications and horizontal gene transfers were evident during plant terrestrialization. Clusterization of genes encoding nitrate assimilation proteins might be an adaptive strategy for algae to cope with changeable nitrate availability in different habitats. Green plants evolved complex nitrate signaling machinery that was stepwise improved by domains shuffling and regulation co-option. Our study highlights innovations in inorganic nitrogen assimilation and signaling machineries, ranging from molecular modifications of proteins to genomic rearrangements, which shaped developmental and metabolic adaptations of plants to changeable nutrient availability in environments.

Origin and early divergence of tandem duplicated sorbitol transporter genes in Rosaceae: insights from evolutionary analysis of the SOT gene family in angiosperms

Sorbitol is a critical photosynthate and storage substance in the Rosaceae family. Sorbitol transporters (SOTs) play a vital role in facilitating sorbitol allocation from source to sink organs and sugar accumulation in sink organs. While prior research has addressed gene duplications within the SOT gene family in Rosaceae, the precise origin and evolutionary dynamics of these duplications remain unclear, largely due to the complicated interplay of whole genome duplications and tandem duplications. Here, we investigated the synteny relationships among all identified Polyol/Monosaccharide Transporter (PLT) genes in 61 angiosperm genomes and SOT genes in representative genomes within the Rosaceae family. By integrating phylogenetic analyses, we elucidated the lineage-specific expansion and syntenic conservation of PLTs and SOTs across diverse plant lineages. We found that Rosaceae SOTs, as PLT family members, originated from a pair of tandemly duplicated PLT genes within Class III-A. Furthermore, our investigation highlights the role of lineage-specific and synergistic duplications in Amygdaloideae in contributing to the expansion of SOTs in Rosaceae plants. Collectively, our findings provide insights into the genomic origins, duplication events, and subsequent divergence of SOT gene family members. Such insights lay a crucial foundation for comprehensive functional characterizations in future studies.

Highly differentiated genomic properties underpin the different cell walls of Poaceae and eudicots

Plant cell walls of Poaceae and eudicots differ substantially, both in the content and composition of their components. However, the genomic and genetic basis underlying these differences is not fully resolved. In this research, we analyzed multiple genomic properties of 150 cell wall gene families across 169 angiosperm genomes. The properties analyzed include gene presence/absence, copy number, synteny, occurrence of tandem gene clusters, and phylogenetic gene diversity. Results revealed a profound genomic differentiation of cell wall genes between Poaceae and eudicots, often associated with the cell wall diversity between these plant groups. For example, overall patterns of gene copy number variation and synteny were clearly divergent between Poaceae and eudicot species. Moreover, differential Poaceae-eudicot copy number and genomic contexts were observed for all the genes within the BEL1-like HOMEODOMAIN 6 regulatory pathway, which respectively induces and represses secondary cell wall synthesis in Poaceae and eudicots. Similarly, divergent synteny, copy number, and phylogenetic gene diversification were observed for the major biosynthetic genes of xyloglucans, mannans, and xylans, potentially contributing to the differences in content and types of hemicellulosic polysaccharides differences in Poaceae and eudicot cell walls. Additionally, the Poaceae-specific tandem clusters and/or higher copy number of PHENYLALANINE AMMONIA-LYASE, CAFFEIC ACID O-METHYLTRANSFERASE, or PEROXIDASE genes may underly the higher content and larger variety of phenylpropanoid compounds observed in Poaceae cell walls. All these patterns are discussed in detail in this study, along with their evolutionary and biological relevance for cell wall (genomic) diversification between Poaceae and eudicots.

Genome-wide identification of late embryogenesis abundant protein family and their key regulatory network in Pinus tabuliformis cold acclimation

Cold acclimation is a crucial biological process that enables conifers to overwinter safely. The late embryogenesis abundant (LEA) protein family plays a pivotal role in enhancing freezing tolerance during this process. Despite its importance, the identification, molecular functions and regulatory networks of the LEA protein family have not been extensively studied in conifers or gymnosperms. Pinus tabuliformis, a conifer with high ecological and economic values and with high-quality genome sequence, is an ideal candidate for such studies. Here, a total of 104 LEA genes were identified from P. tabuliformis, and we renamed them according to their subfamily group: PtLEA1-PtLEA92 (group LEA1-LEA6), PtSMP1-PtSMP6 (group seed maturation protein) and PtDHN1-PtDHN6 (group Dehydrin). While the sequence structure of P. tabuliformis  LEA genes are conserved, their physicochemical properties exhibit unique characteristics within different subfamily groupings. Notably, the abundance of low-temperature responsive elements in PtLEA genes was observed. Using annual rhythm and temperature gradient transcriptome data, PtLEA22 was identified as a key gene that responds to low-temperature induction while conforming to the annual cycle of cold acclimation. Overexpression of PtLEA22 enhanced Arabidopsis freezing tolerance. Furthermore, several transcription factors potentially co-expressed with PtLEA22 were validated using yeast one-hybrid and dual-luciferase assays, revealing that PtDREB1 could directly bind PtLEA22 promoter to positively regulate its expression. These findings reveal the genome-wide characterization of P. tabuliformis  LEA genes and their importance in the cold acclimation, while providing a theoretical basis for studying the molecular mechanisms of cold acclimation in conifers.

Dynamic gene regulatory networks improving spike fertility through regulation of floret primordia fate in wheat

The developmental process of spike is critical for spike fertility through affecting floret primordia fate in wheat; however, the genetic regulation of this dynamic and complex developmental process remains unclear. Here, we conducted a high temporal-resolution analysis of spike transcriptomes and monitored the number and morphology of floret primordia within spike. The development of all floret primordia in a spike was clearly separated into three distinct phases: differentiation, pre-dimorphism and dimorphism. Notably, we identified that floret primordia with meiosis ability at the pre-dimorphism phase usually develop into fertile floret primordia in the next dimorphism phase. Compared to control, increasing plant space treatment achieved the maximum increasement range (i.e., 50%) in number of fertile florets by accelerating spike development. The process of spike fertility improvement was directed by a continuous and dynamic regulatory network involved in transcription factor and genes interaction. This was based on the coordination of genes related to heat shock protein and jasmonic acid biosynthesis during differentiation phase, and genes related to lignin, anthocyanin and chlorophyll biosynthesis during dimorphism phase. The multi-dimensional association with high temporal-resolution approach reported here allows rapid identification of genetic resource for future breeding studies to realise the maximum spike fertility potential in more cereal crops.

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.

Genome sequencing provides insights into the evolution of gene families encoding plant cell wall-degrading enzymes in longhorned beetles

With more than 36,000 species, the longhorned beetles (family Cerambycidae) are a mega-diverse lineage of mostly xylophagous insects, all of which are represented by the sole sequenced genome of the Asian longhorned beetle (Anoplophora glabripennis; Lamiinae). Their successful radiation has been linked to their ability to degrade plant cell wall components using a range of so-called plant cell wall-degrading enzymes (PCWDEs). Our previous analysis of larval gut transcriptomes demonstrated that cerambycid beetles horizontally acquired genes encoding PCWDEs from various microbial donors; these genes evolved through multiple duplication events to form gene families. To gain further insights into the evolution of these gene families during the Cerambycidae radiation, we assembled draft genomes for four beetle species belonging to three subfamilies using long-read nanopore sequencing. All the PCWDE-encoding genes we annotated from the corresponding larval gut transcriptomes were present in these draft genomes. We confirmed that the newly discovered horizontally acquired glycoside hydrolase family 7 (GH7), subfamily 26 of GH43 (GH43_26), and GH53 (all of which are absent from the A. glabripennis genome) were indeed encoded by these beetles' genome. Most of the PCWDE-encoding genes of bacterial origin gained introns after their transfer into the beetle genome. Altogether, we show that draft genome assemblies generated from nanopore long-reads offer meaningful information to the study of the evolution of gene families in insects. We anticipate that our data will support studies aiming to better understand the biology of the Cerambycidae and other beetles in general.

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

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

Pangenome-based trajectories of intracellular gene transfers in Poaceae unveil high cumulation in Triticeae

Intracellular gene transfers (IGTs) between the nucleus and organelles, including plastids and mitochondria, constantly reshape the nuclear genome during evolution. Despite the substantial contribution of IGTs to genome variation, the dynamic trajectories of IGTs at the pangenomic level remain elusive. Here, we developed an approach, IGTminer, that maps the evolutionary trajectories of IGTs using collinearity and gene reannotation across multiple genome assemblies. We applied IGTminer to create a nuclear organellar gene (NOG) map across 67 genomes covering 15 Poaceae species, including important crops. The resulting NOGs were verified by experiments and sequencing data sets. Our analysis revealed that most NOGs were recently transferred and lineage specific and that Triticeae species tended to have more NOGs than other Poaceae species. Wheat (Triticum aestivum) had a higher retention rate of NOGs than maize (Zea mays) and rice (Oryza sativa), and the retained NOGs were likely involved in photosynthesis and translation pathways. Large numbers of NOG clusters were aggregated in hexaploid wheat during 2 rounds of polyploidization, contributing to the genetic diversity among modern wheat accessions. We implemented an interactive web server to facilitate the exploration of NOGs in Poaceae. In summary, this study provides resources and insights into the roles of IGTs in shaping interspecies and intraspecies genome variation and driving plant genome evolution.

Structural and Functional Analysis of the MADS-Box Genes Reveals Their Functions in Cold Stress Responses and Flower Development in Tea Plant (Camellia sinensis)

MADS-box genes comprise a large family of transcription factors that play crucial roles in all aspects of plant growth and development. However, no detailed information on the evolutionary relationship and functional characterization of MADS-box genes is currently available for some representative lineages, such as the Camellia plant. In this study, 136 MADS-box genes were detected from a reference genome of the tea plant (Camellia sinensis) by employing a 569 bp HMM (Hidden Markov Model) developed using nucleotide sequencing including 73 type I and 63 type II genes. An additional twenty-seven genes were identified, with five MIKC-type genes. Truncated and/or inaccurate gene models were manually verified and curated to improve their functional characterization. Subsequently, phylogenetic relationships, chromosome locations, conserved motifs, gene structures, and gene expression profiles were systematically investigated. Tea plant MIKC genes were divided into all 14 major eudicot subfamilies, and no gene was found in Mβ. The expansion of MADS-box genes in the tea plant was mainly contributed by WGD/fragment and tandem duplications. The expression profiles of tea plant MADS-box genes in different tissues and seasons were analyzed, revealing widespread evolutionary conservation and genetic redundancy. The expression profiles linked to cold stress treatments suggested the wide involvement of MADS-box genes from the tea plant in response to low temperatures. Moreover, a floral 'ABCE' model was proposed in the tea plant and proved to be both conserved and ancient. Our analyses offer a detailed overview of MADS-box genes in the tea plant, allowing us to hypothesize the potential functions of unknown genes and providing a foundation for further functional characterizations.

Phylogenomic analysis provides insights into MADS-box and TCP gene diversification and floral development of the Asteraceae, supported by de novo genome and transcriptome sequences from dandelion (Taraxacum officinale)

The Asteraceae is the largest angiosperm family with more than 25,000 species. Individual studies have shown that MADS-box and TCP transcription factors are regulators of the development and symmetry of flowers, contributing to their iconic flower-head (capitulum) and floret. However, a systematic study of MADS-box and TCP genes across the Asteraceae is lacking. We performed a comparative analysis of genome sequences of 33 angiosperm species including our de novo assembly of diploid sexual dandelion (Taraxacum officinale) and 11 other Asteraceae to investigate the lineage-specific evolution of MADS-box and TCP genes in the Asteraceae. We compared the phylogenomic results of MADS-box and TCP genes with their expression in T. officinale floral tissues at different developmental stages to demonstrate the regulation of genes with Asteraceae-specific attributes. Here, we show that MADS-box MIKC^c and TCP-CYCLOIDEA (CYC) genes have expanded in the Asteraceae. The phylogenomic analysis identified AGAMOUS-like (AG-like: SEEDSTICK [STK]-like), SEPALATA-like (SEP3-like), and TCP-PROLIFERATING CELL FACTOR (PCF)-like copies with lineage-specific genomic contexts in the Asteraceae, Cichorioideae, or dandelion. Different expression patterns of some of these gene copies suggest functional divergence. We also confirm the presence and revisit the evolutionary history of previously named “Asteraceae-Specific MADS-box genes (AS-MADS).” Specifically, we identify non-Asteraceae homologs, indicating a more ancient origin of this gene clade. Syntenic relationships support that AS-MADS is paralogous to FLOWERING LOCUS C (FLC) as demonstrated by the shared ancient duplication of FLC and SEP3.

A conserved gene regulatory network controls root epidermal cell patterning in superrosid species

In superrosid species, root epidermal cells differentiate into root hair cells and nonhair cells. In some superrosids, the root hair cells and nonhair cells are distributed randomly (Type I pattern), and in others, they are arranged in a position-dependent manner (Type III pattern). The model plant Arabidopsis (Arabidopsis thaliana) adopts the Type III pattern, and the gene regulatory network (GRN) that controls this pattern has been defined. However, it is unclear whether the Type III pattern in other species is controlled by a similar GRN as in Arabidopsis, and it is not known how the different patterns evolved. In this study, we analyzed superrosid species Rhodiola rosea, Boehmeria nivea, and Cucumis sativus for their root epidermal cell patterns. Combining phylogenetics, transcriptomics, and cross-species complementation, we analyzed homologs of the Arabidopsis patterning genes from these species. We identified R. rosea and B. nivea as Type III species and C. sativus as Type I species. We discovered substantial similarities in structure, expression, and function of Arabidopsis patterning gene homologs in R. rosea and B. nivea, and major changes in C. sativus. We propose that in superrosids, diverse Type III species inherited the patterning GRN from a common ancestor, whereas Type I species arose by mutations in multiple lineages.

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.

Synteny Identifies Reliable Orthologs for Phylogenomics and Comparative Genomics of the Brassicaceae

Large genomic data sets are becoming the new normal in phylogenetic research, but the identification of true orthologous genes and the exclusion of problematic paralogs is still challenging when applying commonly used sequencing methods such as target enrichment. Here, we compared conventional ortholog detection using OrthoFinder with ortholog detection through genomic synteny in a data set of 11 representative diploid Brassicaceae whole-genome sequences spanning the entire phylogenetic space. Then, we evaluated the resulting gene sets regarding gene number, functional annotation, and gene and species tree resolution. Finally, we used the syntenic gene sets for comparative genomics and ancestral genome analysis. The use of synteny resulted in considerably more orthologs and also allowed us to reliably identify paralogs. Surprisingly, we did not detect notable differences between species trees reconstructed from syntenic orthologs when compared with other gene sets, including the Angiosperms353 set and a Brassicaceae-specific target enrichment gene set. However, the synteny data set comprised a multitude of gene functions, strongly suggesting that this method of marker selection for phylogenomics is suitable for studies that value downstream gene function analysis, gene interaction, and network studies. Finally, we present the first ancestral genome reconstruction for the Core Brassicaceae which predating the Brassicaceae lineage diversification ∼25 million years ago.

syntenet: an R/Bioconductor package for the inference and analysis of synteny networks

SUMMARY

Interpreting and visualizing synteny relationships across several genomes is a challenging task. We previously proposed a network-based approach for better visualization and interpretation of large-scale microsynteny analyses. Here, we present syntenet, an R package to infer and analyze synteny networks from whole-genome protein sequence data. The package offers a simple and complete framework, including data preprocessing, synteny detection and network inference, network clustering and phylogenomic profiling, and microsynteny-based phylogeny inference. Graphical functions are also available to create publication-ready plots. Synteny networks inferred with syntenet can highlight taxon-specific gene clusters that likely contributed to the evolution of important traits, and microsynteny-based phylogenies can help resolve phylogenetic relationships under debate.

AVAILABILITY AND IMPLEMENTATION

syntenet is available on Bioconductor (https://bioconductor.org/packages/syntenet), and the source code is available on a GitHub repository (https://github.com/almeidasilvaf/syntenet).

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Contrasting modes of macro and microsynteny evolution in a eukaryotic subphylum

Examination of the changes in order and arrangement of homologous genes is key for understanding the mechanisms of genome evolution in eukaryotes. Previous comparisons between eukaryotic genomes have revealed considerable conservation across species that diverged hundreds of millions of years ago (e.g., vertebrates,1,2,3 bilaterian animals,4,5 and filamentous fungi6). However, understanding how genome organization evolves within and between eukaryotic major lineages remains underexplored. We analyzed high-quality genomes of 120 representative budding yeast species (subphylum Saccharomycotina) spanning ∼400 million years of eukaryotic evolution to examine how their genome organization evolved and to compare it with the evolution of animal and plant genome organization.7 We found that the decay of both macrosynteny (the conservation of homologous chromosomes) and microsynteny (the conservation of local gene content and order) was strongly associated with evolutionary divergence across budding yeast major clades. However, although macrosynteny decayed very fast, within ∼100 million years, the microsynteny of many genes-especially genes in metabolic clusters (e.g., in the GAL gene cluster8)-was much more deeply conserved both within major clades and across the subphylum. We further found that when genomes with similar evolutionary divergence times were compared, budding yeasts had lower macrosynteny conservation than animals and filamentous fungi but higher conservation than angiosperms. In contrast, budding yeasts had levels of microsynteny conservation on par with mammals, whereas angiosperms exhibited very low conservation. Our results provide new insight into the tempo and mode of the evolution of gene and genome organization across an entire eukaryotic subphylum.

A Tea Plant (Camellia sinensis) FLOWERING LOCUS C-like Gene, CsFLC1, Is Correlated to Bud Dormancy and Triggers Early Flowering in Arabidopsis

Flowering and bud dormancy are crucial stages in the life cycle of perennial angiosperms in temperate climates. MADS-box family genes are involved in many plant growth and development processes. Here, we identified three MADS-box genes in tea plant belonging to the FLOWERING LOCUS C (CsFLC) family. We monitored CsFLC1 transcription throughout the year and found that CsFLC1 was expressed at a higher level during the winter bud dormancy and flowering phases. To clarify the function of CsFLC1, we developed transgenic Arabidopsis thaliana plants heterologously expressing 35S::CsFLC1. These lines bolted and bloomed earlier than the WT (Col-0), and the seed germination rate was inversely proportional to the increased CsFLC1 expression level. The RNA-seq of 35S::CsFLC1 transgenic Arabidopsis showed that many genes responding to ageing, flower development and leaf senescence were affected, and phytohormone-related pathways were especially enriched. According to the results of hormone content detection and RNA transcript level analysis, CsFLC1 controls flowering time possibly by regulating SOC1, AGL42, SEP3 and AP3 and hormone signaling, accumulation and metabolism. This is the first time a study has identified FLC-like genes and characterized CsFLC1 in tea plant. Our results suggest that CsFLC1 might play dual roles in flowering and winter bud dormancy and provide new insight into the molecular mechanisms of FLC in tea plants as well as other plant species.

SEP-like genes of Gossypium hirsutum promote flowering via targeting different loci in a concentration-dependent manner

SEP genes are famous for their function in the morphological novelty of bisexual flowers. Although the diverse functions of SEP genes were reported, only the regulatory mechanisms underlying floral organ development have been addressed. In this study, we identified SEP-like genes in Gossypium and found that SEP3 genes were duplicated in diploid cotton varieties. GhSEP4.1 and GhSEP4.2 were abundantly transcribed in the shoot apical meristem (SAM), but only GhSEP4.2 was expressed in the leaf vasculature. The expression pattern of GhSEPs in floral organs was conserved with that of homologs in Arabidopsis, except for GhSEP2 that was preponderantly expressed in ovules and fibers. The overexpression and silencing of each single GhSEP gene suggested their distinct role in promoting flowering via direct binding to GhAP1 and GhLFY genomic regions. The curly leaf and floral defects in overexpression lines with a higher expression of GhSEP genes revealed the concentration-dependent target gene regulation of GhSEP proteins. Moreover, GhSEP proteins were able to dimerize and interact with flowering time regulators. Together, our results suggest the dominant role of GhSEP4.2 in leaves to promote flowering via GhAP1-A04, and differently accumulated GhSEP proteins in the SAM alternately participate in forming the dynamic tetramer complexes to target at the different loci of GhAP1 and GhLFY to maintain reproductive growth. The regulatory roles of cotton SEP genes reveal their conserved and diversified functions.

The chromosome-level genome of female ginseng (Angelica sinensis) provides insights into molecular mechanisms and evolution of coumarin biosynthesis

Coumarins are natural products with important medicinal values, and include simple coumarins, furanocoumarins and pyranocoumarins. Female ginseng (Angelica sinensis) is a renowned herb with abundant coumarins, originated in China and known for the treatment of female ailments for thousands of years. The molecular basis of simple coumarin biosynthesis in A. sinensis and the evolutionary history of the genes involved in furanocoumarin biosynthesis are largely unknown. Here, we generated the first chromosome-scale genome of A. sinensis. It has a genome size of 2.37 Gb, which was generated by combining PacBio and Hi-C sequencing technologies. The genome was predicted to contain 43 202 protein-coding genes dispersed mainly on 11 pseudochromosomes. We not only provided evidence for whole-genome duplication (WGD) specifically occurring in the Apioideae subfamily, but also demonstrated the vital role of tandem duplication for phenylpropanoid biosynthesis in A. sinensis. Combined analyses of transcriptomic and metabolomic data revealed key genes and candidate transcription factors regulating simple coumarin biosynthesis. Furthermore, phylogenomic synteny network analyses suggested prenyltransferase genes involved in furanocoumarin biosynthesis evolved independently in the Moraceae, Fabaceae, Rutaceae and Apiaceae after ζ and ε WGD. Our work sheds light on coumarin biosynthesis, and provides a benchmark for accelerating genetic research and molecular breeding in A. sinensis.

Genome and transcriptome-wide study of carbamoyltransferase genes in major fleshy fruits: A multi-omics study of evolution and functional significance

The carbamoyltransferase or aspartate carbamoyltransferase (ATCase)/ornithine carbamoyltransferase (OTCase) is an evolutionary conserved protein family, which contains two genes, ATCase and OTCase. The ATCase catalyzes the committed step in the synthesis of UMP from which all pyrimidine molecules are synthesized. The second member, OTCase, catalytically regulates the conversion of ornithine to citrulline. This study traces the evolution of the carbomoyltransferase genes in the plant kingdom and their role during fruit ripening in fleshy fruits. These genes are highly conserved throughout the plant kingdom and, except for melon and watermelon, do not show gene expansion in major fleshy fruits. In this study, 393 carbamoyltransferase genes were identified in the plant kingdom, including 30 fleshy fruit representatives. Their detailed phylogeny, evolutionary patterns with their expression during the process of fruit ripening, was analyzed. The ATcase and OTcase genes were conserved throughout the plant kingdom and exhibited lineage-specific signatures. The expression analysis of the ATcase and OTcase genes during fruit development and ripening in climacteric and non-climacteric fruits showed their involvement in fruit ripening irrespective of the type of fruits. No direct role in relation to ethylene-dependent or -independent ripening was identified; however, the co-expression network suggests their involvement in the various ripening processes.

Tracing the Evolution of the SEPALLATA Subfamily across Angiosperms Associated with Neo- and Sub-Functionalization for Reproductive and Agronomically Relevant Traits

SEPALLATA transcription factors (SEP TFs) have been extensively studied in angiosperms as pivotal components of virtually all the MADS-box tetrameric complex master regulators of floral organ identities. However, there are published reports that suggest that some SEP members also regulate earlier reproductive events, such as inflorescence meristem determinacy and inflorescence architecture, with potential for application in breeding programs in crops. The SEP subfamily underwent a quite complex pattern of duplications during the radiation of the angiosperms. Taking advantage of the many whole genomic sequences now available, we present a revised and expanded SEP phylogeny and link it to the known functions of previously characterized genes. This snapshot supports the evidence that the major SEP3 clade is highly specialized for the specification of the three innermost floral whorls, while its sister LOFSEP clade is functionally more versatile and has been recruited for diverse roles, such as the regulation of extra-floral bract formation and inflorescence determinacy and shape. This larger pool of angiosperm SEP genes confirms previous evidence that their evolution was driven by whole-genome duplications rather than small-scale duplication events. Our work may help to identify those SEP lineages that are the best candidates for the improvement of inflorescence traits, even in far distantly related crops.

MbHY5-MbYSL7 mediates chlorophyll synthesis and iron transport under iron deficiency in Malus baccata

Iron (Fe) plays an important role in cellular respiration and catalytic reactions of metalloproteins in plants and animals. Plants maintain iron homeostasis through absorption, translocation, storage, and compartmentalization of iron via a cooperative regulative network. Here, we showed different physiological characteristics in the leaves and roots of Malus baccata under Fe sufficiency and Fe deficiency conditions and propose that MbHY5 (elongated hypocotyl 5), an important transcription factor for its function in photomorphogenesis, participated in Fe deficiency response in both the leaves and roots of M. baccata. The gene co-expression network showed that MbHY5 was involved in the regulation of chlorophyll synthesis and Fe transport pathway under Fe-limiting conditions. Specifically, we found that Fe deficiency induced the expression of MbYSL7 in root, which was positively regulated by MbHY5. Overexpressing or silencing MbYSL7 influenced the expression of MbHY5 in M. baccata.

The membrane associated NAC transcription factors ANAC060 and ANAC040 are functionally redundant in the inhibition of seed dormancy in Arabidopsis thaliana

The NAC family of transcription factors is involved in plant development and various biotic and abiotic stresses. The Arabidopsis thaliana ANAC genes ANAC060, ANAC040, and ANAC089 are highly homologous based on protein and nucleotide sequence similarity. These three genes are predicted to be membrane bound transcription factors (MTFs) containing a conserved NAC domain, but divergent C-terminal regions. The anac060 mutant shows increased dormancy when compared with the wild type. Mutations in ANAC040 lead to higher seed germination under salt stress, and a premature stop codon in ANAC089 Cvi allele results in seeds exhibiting insensitivity to high concentrations of fructose. Thus, these three homologous MTFs confer distinct functions, although all related to germination. To investigate whether the differences in function are caused by a differential spatial or temporal regulation, or by differences in the coding sequence (CDS), we performed swapping experiments in which the promoter and CDS of the three MTFs were exchanged. Seed dormancy and salt and fructose sensitivity analyses of transgenic swapping lines in mutant backgrounds showed that there is functional redundancy between ANAC060 and ANAC040, but not between ANAC060 and ANAC089.

The origin and evolution of RNase T2 family and gametophytic self-incompatibility system in plants

Ribonuclease (RNase) T2 genes are found widely in both eukaryotes and prokaryotes, and genes from this family have been revealed to have various functions in plants. In particular, S-RNase is known to be the female determinant in the S-RNase-based gametophytic self-incompatibility (GSI) system. However, the origin and evolution of the RNase T2 gene family and GSI system are not well understood. In this study, 785 RNase T2 genes were identified in 81 sequenced plant genomes representing broad-scale diversity and divided into three subgroups (Class I, II, and III) based on phylogenetic and synteny network analysis. Class I was found to be of ancient origin and to emerge in green algae, Class II was shown to originate with the appearance of angiosperms, while Class III was discovered to be eudicot-specific. Each of the three major classes could be further classified into several subclasses of which some subclasses were found to be lineage-specific. Furthermore, duplication, deletion, or inactivation of the S/S-like-locus was revealed to be linked to repeated loss and gain of self-incompatibility in different species from distantly related plant families with GSI. Finally, the origin and evolutionary history of S-locus in Rosaceae species was unraveled with independent loss and gain of S-RNase occurred in different subfamilies of Rosaceae. Our findings provide insights into the origin and evolution of the RNase T2 family and the GSI system in plants.

Genomic Architecture and Evolution of the Cellulose synthase Gene Superfamily as Revealed by Phylogenomic Analysis

The Cellulose synthase superfamily synthesizes cellulose and different hemicellulosic polysaccharides in plant cell walls. While much has been discovered about the evolution and function of these genes, their genomic architecture and relationship with gene (sub-)functionalization and evolution remains unclear. By using 242 genomes covering plant evolution from green algae to eudicots, we performed a large-scale analysis of synteny, phylogenetic, and functional data of the CesA superfamily. Results revealed considerable gene copy number variation across species and gene families, and also two patterns - singletons vs. tandem arrays - in chromosomic gene arrangement. Synteny analysis revealed exceptional conservation of gene architecture across species, but also lineage-specific patterns across gene (sub-)families. Synteny patterns correlated with gene sub-functionalization into primary and secondary CesAs and distinct CslD functional isoforms. Furthermore, a genomic context shift of a group of cotton secondary CesAs was associated with peculiar properties of cotton fiber synthesis. Finally, phylogenetics suggested that primary CesA sequences appeared before the secondary CesAs, while phylogenomic analyses unveiled the genomic trace of the CslD duplication that initiated the CslF family. Our results describe in detail the genomic architecture of the CesA superfamily in plants, highlighting its crucial relevance for gene diversification and sub-functionalization, and for understanding their evolution.

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 B₆ 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.

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.

Firming up your tomato: a natural promoter variation in a MADS-box gene is causing all-flesh tomatoes

This article comments on: Liu L, Zhang K, Bai JR, Lu J, Lu X, Hu J, Pan C, He S, Yuan J, Zhang Y, Zhang M, Guo Y, Wang X, Huang Z, Du Y, Cheng F, Li J. 2022. All-flesh fruit in tomato is controlled by reduced expression dosage of AFF through a structural variant mutation in the promoter. Journal of Experimental Botany 73, 123–138.

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.

New insights into the evolution of SPX gene family from algae to legumes; a focus on soybean

BACKGROUND

SPX-containing proteins have been known as key players in phosphate signaling and homeostasis. In Arabidopsis and rice, functions of some SPXs have been characterized, but little is known about their function in other plants, especially in the legumes.

RESULTS

We analyzed SPX gene family evolution in legumes and in a number of key species from algae to angiosperms. We found that SPX harboring proteins showed fluctuations in domain fusions from algae to the angiosperms with, finally, four classes appearing and being retained in the land plants. Despite these fluctuations, Lysine Surface Cluster (KSC), and the third residue of Phosphate Binding Sites (PBS) showed complete conservation in almost all of SPXs except few proteins in Selaginella moellendorffii and Papaver sumniferum, suggesting they might have different ligand preferences. In addition, we found that the WGD/segmentally or dispersed duplication types were the most frequent contributors to the SPX expansion, and that there is a positive correlation between the amount of WGD contribution to the SPX expansion in individual species and its number of EXS genes. We could also reveal that except SPX class genes, other classes lost the collinearity relationships among Arabidopsis and legume genomes. The sub- or neo-functionalization of the duplicated genes in the legumes makes it difficult to find the functional orthologous genes. Therefore, we used two different methods to identify functional orthologs in soybean and Medicago. High variance in the dynamic and spatial expression pattern of GmSPXs proved the new or sub-functionalization in the paralogs.

CONCLUSION

This comprehensive analysis revealed how SPX gene family evolved from algae to legumes and also discovered several new domains fused to SPX domain in algae. In addition, we hypothesized that there different phosphate sensing mechanisms might occur in S. moellendorffii and P. sumniferum. Finally, we predicted putative functional orthologs of AtSPXs in the legumes, especially, orthologs of AtPHO1, involved in long-distance Pi transportation. These findings help to understand evolution of phosphate signaling and might underpin development of new legume varieties with improved phosphate use efficiency.

Phylogenomic and Microsynteny Analysis Provides Evidence of Genome Arrangements of High-Affinity Nitrate Transporter Gene Families of Plants

Nitrate transporter 2 (NRT2) and NRT3 or nitrate-assimilation-related 2 (NAR2) proteins families form a two-component, high-affinity nitrate transport system, which is essential for the acquisition of nitrate from soils with low N availability. An extensive phylogenomic analysis across land plants for these families has not been performed. In this study, we performed a microsynteny and orthology analysis on the NRT2 and NRT3 genes families across 132 plants (Sensu lato) to decipher their evolutionary history. We identified significant differences in the number of sequences per taxonomic group and different genomic contexts within the NRT2 family that might have contributed to N acquisition by the plants. We hypothesized that the greater losses of NRT2 sequences correlate with specialized ecological adaptations, such as aquatic, epiphytic, and carnivory lifestyles. We also detected expansion on the NRT2 family in specific lineages that could be a source of key innovations for colonizing contrasting niches in N availability. Microsyntenic analysis on NRT3 family showed a deep conservation on land plants, suggesting a high evolutionary constraint to preserve their function. Our study provides novel information that could be used as guide for functional characterization of these gene families across plant lineages.

Reproductive developmental transcriptome analysis of Tripidium ravennae (Poaceae)

BACKGROUND

Tripidium ravennae is a cold-hardy, diploid species in the sugarcane complex (Poaceae subtribe Saccharinae) with considerable potential as a genetic resource for developing improved bioenergy and ornamental grasses. An improved understanding of the genetic regulation of reproductive processes (e.g., floral induction, inflorescence development, and seed development) will enable future applications of precision breeding and gene editing of floral and seed development. In particular, the ability to silence reproductive processes would allow for developing seedless forms of valuable but potentially invasive plants. The objective of this research was to characterize the gene expression environment of reproductive development in T. ravennae.

RESULTS

During the early phases of inflorescence development, multiple key canonical floral integrators and pathways were identified. Annotations of type II subfamily of MADS-box transcription factors, in particular, were over-represented in the GO enrichment analyses and tests for differential expression (FDR p-value < 0.05). The differential expression of floral integrators observed in the early phases of inflorescence development diminished prior to inflorescence determinacy regulation. Differential expression analysis did not identify many unique genes at mid-inflorescence development stages, though typical biological processes involved in plant growth and development expressed abundantly. The increase in inflorescence determinacy regulatory elements and putative homeotic floral development unigenes at mid-inflorescence development coincided with the expression of multiple meiosis annotations and multicellular organism developmental processes. Analysis of seed development identified multiple unigenes involved in oxidative-reductive processes.

CONCLUSION

Reproduction in grasses is a dynamic system involving the sequential coordination of complex gene regulatory networks and developmental processes. This research identified differentially expressed transcripts associated with floral induction, inflorescence development, and seed development in T. ravennae. These results provide insights into the molecular regulation of reproductive development and provide a foundation for future investigations and analyses, including genome annotation, functional genomics characterization, gene family evolutionary studies, comparative genomics, and precision breeding.

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.

Whole-genome microsynteny-based phylogeny of angiosperms

Plant genomes vary greatly in size, organization, and architecture. Such structural differences may be highly relevant for inference of genome evolution dynamics and phylogeny. Indeed, microsynteny-the conservation of local gene content and order-is recognized as a valuable source of phylogenetic information, but its use for the inference of large phylogenies has been limited. Here, by combining synteny network analysis, matrix representation, and maximum likelihood phylogenetic inference, we provide a way to reconstruct phylogenies based on microsynteny information. Both simulations and use of empirical data sets show our method to be accurate, consistent, and widely applicable. As an example, we focus on the analysis of a large-scale whole-genome data set for angiosperms, including more than 120 available high-quality genomes, representing more than 50 different plant families and 30 orders. Our 'microsynteny-based' tree is largely congruent with phylogenies proposed based on more traditional sequence alignment-based methods and current phylogenetic classifications but differs for some long-contested and controversial relationships. For instance, our synteny-based tree finds Vitales as early diverging eudicots, Saxifragales within superasterids, and magnoliids as sister to monocots. We discuss how synteny-based phylogenetic inference can complement traditional methods and could provide additional insights into some long-standing controversial phylogenetic relationships.

Plant design gets its details: Modulating plant architecture by phase transitions

Plants evolved different strategies to better adapt to the environmental conditions in which they live: the control of their body architecture and the timing of phase change are two important processes that can improve their fitness. As they age, plants undergo two major phase changes (juvenile to adult and adult to reproductive) that are a response to environmental and endogenous signals. These phase transitions are accompanied by alterations in plant morphology and also by changes in physiology and the behavior of gene regulatory networks. Six main pathways involving environmental and endogenous cues that crosstalk with each other have been described as responsible for the control of plant phase transitions: the photoperiod pathway, the autonomous pathway, the vernalization pathway, the temperature pathway, the GA pathway, and the age pathway. However, studies have revealed that sugar is also involved in phase change and the control of branching behavior. In this review, we discuss recent advances in plant biology concerning the genetic and molecular mechanisms that allow plants to regulate phase transitions in response to the environment. We also propose connections between phase transition and plant architecture control.

Kingdom-wide analysis of the evolution of the plant type III polyketide synthase superfamily

The emergence of type III polyketide synthases (PKSs) was a prerequisite for the conquest of land by the green lineage. Within the PKS superfamily, chalcone synthases (CHSs) provide the entry point reaction to the flavonoid pathway, while LESS ADHESIVE POLLEN 5 and 6 (LAP5/6) provide constituents of the outer exine pollen wall. To study the deep evolutionary history of this key family, we conducted phylogenomic synteny network and phylogenetic analyses of whole-genome data from 126 species spanning the green lineage including Arabidopsis thaliana, tomato (Solanum lycopersicum), and maize (Zea mays). This study thereby combined study of genomic location and context with changes in gene sequences. We found that the two major clades, CHS and LAP5/6 homologs, evolved early by a segmental duplication event prior to the divergence of Bryophytes and Tracheophytes. We propose that the macroevolution of the type III PKS superfamily is governed by whole-genome duplications and triplications. The combined phylogenetic and synteny analyses in this study provide insights into changes in the genomic location and context that are retained for a longer time scale with more recent functional divergence captured by gene sequence alterations.

All type II classic MADS-box genes in the lycophyte Selaginella moellendorffii are broadly yet discretely expressed in vegetative and reproductive tissues

The MADS-box genes constitute a large transcription factor family that appear to have evolved by duplication and diversification of function. Two types of MADS-box genes are distinguished throughout eukaryotes, types I and II. Type II classic MADS-box genes, also known as MIKC-type, are key developmental regulators in flowering plants and are particularly well-studied for their role in floral organ specification. However, very little is known about the role that these genes might play outside of the flowering plants. We investigated the evolution of type II classic MADS-box genes across land plants by performing a maximum likelihood analysis with a particular focus on lycophytes. Here, we present the expression patterns of all three type II classic MADS-box homologs throughout plant development in the lycophyte Selaginella moellendorffii: SmMADS1, SmMADS3, and SmMADS6. We used scanning electron microscopy and histological analyses to define stages of sporangia development in S. moellendorffii. We performed phylogenetic analyses of this gene lineage across land plants and found that lycophyte sequences appeared before the multiple duplication events that gave rise to the major MADS-box gene lineages in seed plants. Our expression analyses by in situ hybridization show that all type II classic MADS-box genes in S. moellendorffii have broad but distinct patterns of expression in vegetative and reproductive tissues, where SmMADS1 and SmMADS6 only differ during late sporangia development. The broad expression during S. moellendorffii development suggests that MADS-box genes have undergone neofunctionalization and subfunctionalization after duplication events in seed plants.

A gene expression atlas for kiwifruit (Actinidia chinensis) and network analysis of transcription factors

BACKGROUND

Transcriptomic studies combined with a well annotated genome have laid the foundations for new understanding of molecular processes. Tools which visualise gene expression patterns have further added to these resources. The manual annotation of the Actinidia chinensis (kiwifruit) genome has resulted in a high quality set of 33,044 genes. Here we investigate gene expression patterns in diverse tissues, visualised in an Electronic Fluorescent Pictograph (eFP) browser, to study the relationship of transcription factor (TF) expression using network analysis.

RESULTS

Sixty-one samples covering diverse tissues at different developmental time points were selected for RNA-seq analysis and an eFP browser was generated to visualise this dataset. 2839 TFs representing 57 different classes were identified and named. Network analysis of the TF expression patterns separated TFs into 14 different modules. Two modules consisting of 237 TFs were correlated with floral bud and flower development, a further two modules containing 160 TFs were associated with fruit development and maturation. A single module of 480 TFs was associated with ethylene-induced fruit ripening. Three "hub" genes correlated with flower and fruit development consisted of a HAF-like gene central to gynoecium development, an ERF and a DOF gene. Maturing and ripening hub genes included a KNOX gene that was associated with seed maturation, and a GRAS-like TF.

CONCLUSIONS

This study provides an insight into the complexity of the transcriptional control of flower and fruit development, as well as providing a new resource to the plant community. The Actinidia eFP browser is provided in an accessible format that allows researchers to download and work internally.

Genomic Characterization of WRKY Transcription Factors Related to Andrographolide Biosynthesis in Andrographis paniculata

Andrographolide, which is enriched in the leaves of Andrographis paniculata, has been known as “natural antibiotic” due to its pharmacological activities such as anti-inflammatory, antimicrobial and antioxidant effects. Several key enzymes in andrographolide biosynthetic pathway have been studied since the genome sequences were released, but its regulatory mechanism remains unknown. WRKY transcription factors proteins have been reported to regulate plant secondary metabolism, development as well as biotic and abiotic stresses. Here, WRKY transcription factors related to andrographolide biosynthesis were systematically identified, including sequences alignment, phylogenetic analysis, chromosomal distribution, gene structure, conserved motifs, synteny, alternative splicing event and Gene ontology (GO) annotation. A total of 58 WRKYs were identified in Chuanxinlian genome and phylogenetically classified into three groups. Moreover, nine WRKY genes underwent alternative splicing events. Furthermore, the combination of binding site prediction, gene-specific expression patterns, and phylogenetic analysis suggested that 7 WRKYs (ApWRKY01, ApWRKY08, ApWRKY12, ApWRKY14, ApWRKY19, ApWRKY20, and ApWRKY50) might regulate andrographolide biosynthesis. This study laid a foundation for understanding the regulatory mechanism of andrographolide biosynthesis and the improvement and breeding of Andrographis paniculata varieties.

Far-Red Light-Induced Azolla filiculoides Symbiosis Sexual Reproduction: Responsive Transcripts of Symbiont Nostoc azollae Encode Transporters Whilst Those of the Fern Relate to the Angiosperm Floral Transition

Water ferns of the genus Azolla and the filamentous cyanobacteria Nostoc azollae constitute a model symbiosis that enabled the colonization of the water surface with traits highly desirable for the development of more sustainable crops: their floating mats capture CO₂ and fix N₂ at high rates using light energy. Their mode of sexual reproduction is heterosporous. The regulation of the transition from the vegetative phase to the spore forming phase in ferns is largely unknown, yet a prerequisite for Azolla domestication, and of particular interest as ferns represent the sister lineage of seed plants. Sporocarps induced with far red light could be crossed so as to verify species attribution of strains from the Netherlands but not of the strain from the Anzali lagoon in Iran; the latter strain was assigned to a novel species cluster from South America. Red-dominated light suppresses the formation of dissemination stages in both gametophyte- and sporophyte-dominated lineages of plants, the response likely is a convergent ecological strategy to open fields. FR-responsive transcripts included those from MIKC^C homologues of CMADS1 and miR319-controlled GAMYB transcription factors in the fern, transporters in N. azollae, and ycf2 in chloroplasts. Loci of conserved microRNA (miRNA) in the fern lineage included miR172, yet FR only induced miR529 and miR535, and reduced miR319 and miR159. Phylogenomic analyses of MIKC^C TFs suggested that the control of flowering and flower organ specification may have originated from the diploid to haploid phase transition in the homosporous common ancestor of ferns and seed plants.

Characteristics of banana B genome MADS-box family demonstrate their roles in fruit development, ripening, and stress

MADS-box genes are critical regulators of growth and development in flowering plants. Sequencing of the Musa balbisiana (B) genome has provided a platform for the systematic analysis of the MADS-box gene family in the important banana ancestor Musa balbisiana. Seventy-seven MADS-box genes, including 18 type I and 59 type II, were strictly identified from the banana (Pisang Klutuk Wulung, PKW, 2n = 2x = 22) B genome. These genes have been preferentially placed on the banana B genome. Evolutionary analysis suggested that M. balbisiana MCM1-AGAMOUS-DEFICIENS-SRF (MbMADS) might be organized into the MIKC^c, MIKC*, Mα, Mβ, and Mγ groups according to the phylogeny. MIKC^c was then further categorized into 10 subfamilies according to conserved motif and gene structure analyses. The well-defined MADS-box genes highlight gene birth and death in banana. MbMADSes originated from the same ancestor as MaMADSes. Transcriptome analysis in cultivated banana (ABB) revealed that MbMADSes were conserved and differentially expressed in several organs, in various fruit developing and ripening stages, and in stress treatments, indicating the participation of these genes in fruit development, ripening, and stress responses. Of note, SEP/AGL2 and AG, as well as other several type II MADS-box genes, including the STMADS11 and TM3/SOC1 subfamilies, indicated elevated expression throughout banana fruit development, ripening, and stress treatments, indicating their new parts in controlling fruit development and ripening. According to the co-expression network analysis, MbMADS75 interacted with bZIP and seven other transcription factors to perform its function. This systematic analysis reveals fruit development, ripening, and stress candidate MbMADSes genes for additional functional studies in plants, improving our understanding of the transcriptional regulation of MbMADSes genes and providing a base for genetic modification of MADS-mediated fruit development, ripening, and stress.

The Transcriptome of Cunninghamia lanceolata male/female cone reveal the association between MIKC MADS-box genes and reproductive organs development

BACKGROUND

Cunninghamia lanceolata (Chinese fir), a member of the conifer family Cupressaceae, is one of the most popular cultivated trees for wood production in China. Continuous research is being performed to improve C. lanceolata breeding values. Given the high rate of seed abortion (one of the reasons being the failure of ovule and pollen development) in C. lanceolata, the proper formation of female/male cones could theoretically increase the number of offspring in future generations. MIKC MADS-box genes are well-known for their roles in the flower/cone development and comprise the typical/atypical floral development model for both angiosperms and gymnosperms.

RESULTS

We performed a transcriptomic analysis to find genes differentially expressed between female and male cones at a single, carefully determined developmental stage, focusing on the MIKC MADS-box genes. We finally obtained 47 unique MIKC MADS-box genes from C. lanceolata and divided these genes into separate branches. 27 out of the 47 MIKC MADS-box genes showed differential expression between female and male cones, and most of them were not expressed in leaves. Out of these 27 genes, most B-class genes (AP3/PI) were up-regulated in the male cone, while TM8 genes were up-regulated in the female cone. Then, with no obvious overall preference for AG (class C + D) genes in female/male cones, it seems likely that these genes are involved in the development of both cones. Finally, a small number of genes such as GGM7, SVP, AGL15, that were specifically expressed in female/male cones, making them candidate genes for sex-specific cone development.

CONCLUSIONS

Our study identified a number of MIKC MADS-box genes showing differential expression between female and male cones in C. lanceolata, illustrating a potential link of these genes with C. lanceolata cone development. On the basis of this, we postulated a possible cone development model for C. lanceolata. The gene expression library showing differential expression between female and male cones shown here, can be used to discover unknown regulatory networks related to sex-specific cone development in the future.

Fibrillarin evolution through the Tree of Life: Comparative genomics and microsynteny network analyses provide new insights into the evolutionary history of Fibrillarin

Fibrillarin (FIB), a methyltransferase essential for life in the vast majority of eukaryotes, is involved in methylation of rRNA required for proper ribosome assembly, as well as methylation of histone H2A of promoter regions of rRNA genes. RNA viral progression that affects both plants and animals requires FIB proteins. Despite the importance and high conservation of fibrillarins, there little is known about the evolutionary dynamics of this small gene family. We applied a phylogenomic microsynteny-network approach to elucidate the evolutionary history of FIB proteins across the Tree of Life. We identified 1063 non-redundant FIB sequences across 1049 completely sequenced genomes from Viruses, Bacteria, Archaea, and Eukarya. FIB is a highly conserved single-copy gene through Archaea and Eukarya lineages, except for plants, which have a gene family expansion due to paleopolyploidy and tandem duplications. We found a high conservation of the FIB genomic context during plant evolution. Surprisingly, FIB in mammals duplicated after the Eutheria split (e.g., ruminants, felines, primates) from therian mammals (e.g., marsupials) to form two main groups of sequences, the FIB and FIB-like groups. The FIB-like group transposed to another genomic context and remained syntenic in all the eutherian mammals. This transposition correlates with differences in the expression patterns of FIB-like proteins and with elevated Ks values potentially due to reduced evolutionary constraints of the duplicated copy. Our results point to a unique evolutionary event in mammals, between FIB and FIB-like genes, that led to non-redundant roles of the vital processes in which this protein is involved.

Phylogenomic analysis of the APETALA2 transcription factor subfamily across angiosperms reveals both deep conservation and lineage-specific patterns

The APETALA2 (AP2) subfamily of transcription factors are key regulators of angiosperm root, shoot, flower and embryo development. The broad diversity of anatomical and morphological structures is potentially associated with the genomic dynamics of the AP2 subfamily. However, a comprehensive phylogenomic analysis of the AP2 subfamily across angiosperms is lacking. We combined phylogenetic and synteny analysis of distinct AP2 subclades in the completed genomes of 107 angiosperm species. We identified major changes in copy number variation and genomic context within subclades across lineages, and discuss how these changes may have contributed to the evolution of lineage-specific traits. Multiple AP2 subclades show highly conserved patterns of copy number and synteny across angiosperms, while others are more dynamic and show distinct lineage-specific patterns. As examples of lineage-specific morphological divergence due to AP2 subclade dynamics, we hypothesize that loss of PLETHORA1/2 in monocots correlates with the absence of taproots, whereas independent lineage-specific changes of PLETHORA4/BABY BOOM and WRINKLED1 genes in Brassicaceae and monocots point towards regulatory divergence of embryogenesis between these lineages. Additionally, copy number expansion of TOE1 and TOE3/AP2 in asterids is implicated with differential regulation of flower development. Moreover, we show that the genomic context of AP2s is in general highly specialized per angiosperm lineage. To our knowledge, this study is the first to shed light on the evolutionary divergence of the AP2 subfamily subclades across major angiosperm lineages and emphasizes the need for lineage-specific characterization of developmental networks to understand trait variability further.

Systematic analysis of 1298 RNA-Seq samples and construction of a comprehensive soybean (Glycine max) expression atlas

Soybean (Glycine max [L.] Merr.) is a major crop in animal feed and human nutrition, mainly for its rich protein and oil contents. The remarkable rise in soybean transcriptome studies over the past 5 years generated an enormous amount of RNA-seq data, encompassing various tissues, developmental conditions and genotypes. In this study, we have collected data from 1298 publicly available soybean transcriptome samples, processed the raw sequencing reads and mapped them to the soybean reference genome in a systematic fashion. We found that 94% of the annotated genes (52 737/56 044) had detectable expression in at least one sample. Unsupervised clustering revealed three major groups, comprising samples from aerial, underground and seed/seed-related parts. We found 452 genes with uniform and constant expression levels, supporting their roles as housekeeping genes. On the other hand, 1349 genes showed heavily biased expression patterns towards particular tissues. A transcript-level analysis revealed that 95% (70 963 of 74 490) of the assembled transcripts have intron chains exactly matching those from known transcripts, whereas 3256 assembled transcripts represent potentially novel splicing isoforms. The dataset compiled here constitute a new resource for the community, which can be downloaded or accessed through a user-friendly web interface at http://venanciogroup.uenf.br/resources/. This comprehensive transcriptome atlas will likely accelerate research on soybean genetics and genomics.