Conserved genomic context provides critical information for comparative evolutionary analysis. With the increase in numbers of sequenced plant genomes, synteny analysis can provide new insights into gene family evolution. Here, we exploit a network analysis approach to organize and interpret massive pairwise syntenic relationships. Specifically, we analyzed synteny networks of the MADS-box transcription factor gene family using 51 completed plant genomes. In combination with phylogenetic profiling, several novel evolutionary patterns were inferred and visualized from synteny network clusters. We found lineage-specific clusters that derive from transposition events for the regulators of floral development (APETALA3 and PI) and flowering time (FLC) in the Brassicales and for the regulators of root development (AGL17) in Poales. We also identified two large gene clusters that jointly encompass many key phenotypic regulatory Type II MADS-box gene clades (SEP1, SQUA, TM8, SEP3, FLC, AGL6, and TM3). Gene clustering and gene trees support the idea that these genes are derived from an ancient tandem gene duplication that likely predates the radiation of the seed plants and then expanded by subsequent polyploidy events. We also identified angiosperm-wide conservation of synteny of several other less studied clades. Combined, these findings provide new hypotheses for the genomic origins, biological conservation, and divergence of MADS-box gene family members.

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.

A comprehensive analysis of relative gene order, or microsynteny, can provide valuable information for understanding the evolutionary history of genes and genomes, and ultimately traits and species, across broad phylogenetic groups and divergence times. We have used our network-based phylogenomic synteny analysis pipeline to first analyze the overall patterns and major differences between 87 mammalian and 107 angiosperm genomes. These two important groups have both evolved and radiated over the last ∼170 MYR. Secondly, we identified the genomic outliers or “rebel genes” within each clade. We theorize that rebel genes potentially have influenced trait and lineage evolution. Microsynteny networks use genes as nodes and syntenic relationships between genes as edges. Networks were decomposed into clusters using the Infomap algorithm, followed by phylogenomic copy-number profiling of each cluster. The differences in syntenic properties of all annotated gene families, including BUSCO genes, between the two clades are striking: most genes are single copy and syntenic across mammalian genomes, whereas most genes are multicopy and/or have lineage-specific distributions for angiosperms. We propose microsynteny scores as an alternative and complementary metric to BUSCO for assessing genome assemblies. We further found that the rebel genes are different between the two groups: lineage-specific gene transpositions are unusual in mammals, whereas single-copy highly syntenic genes are rare for flowering plants. We illustrate several examples of mammalian transpositions, such as brain-development genes in primates, and syntenic conservation across angiosperms, such as single-copy genes related to photosynthesis. Future experimental work can test if these are indeed rebels with a cause.

Freezing limits plant growth and crop productivity, and plant species in temperate zones have the capacity to develop freezing tolerance through complex modulation of gene expression affecting various aspects of metabolism and physiology. While many components of freezing tolerance have been identified in model species under controlled laboratory conditions, little is known about the mechanisms that impart freezing tolerance in natural populations of wild species. Here, we performed a quantitative trait locus (QTL) study of acclimated freezing tolerance in seedlings of Boechera stricta, a highly adapted relative of Arabidopsis (Arabidopsis thaliana) native to the Rocky Mountains. A single QTL was identified that contained the gene encoding ACYL-COENZYME A:DIACYLGLYCEROL ACYLTRANSFERASE1 (BstDGAT1), whose expression is highly cold responsive. The primary metabolic enzyme DGAT1 catalyzes the final step in assembly of triacylglycerol (TAG) by acyl transfer from acyl-CoA to diacylglycerol. Freezing tolerant plants showed higher DGAT1 expression during cold acclimation than more sensitive plants, and this resulted in increased accumulation of TAG in response to subsequent freezing. Levels of oligogalactolipids that are produced by SFR2 (SENSITIVE TO FREEZING2), an indispensable element of freezing tolerance in Arabidopsis, were also higher in freezing-tolerant plants. Furthermore, overexpression of AtDGAT1 led to increased freezing tolerance. We propose that DGAT1 confers freezing tolerance in plants by supporting SFR2-mediated remodeling of chloroplast membranes.

Nucleotide-binding leucine-rich repeat receptor (NLR) genes encode a pivotal class of plant immune receptors. However, their rampant duplication and loss have made inferring their genomic evolutionary trajectory difficult, exemplified by the loss of TNL family genes in monocots. In this study, we introduce a novel classification system for angiosperm NLR genes, grounded in network analysis of microsynteny information. This refined classification categorizes these genes into five classes: CNL_A, CNL_B, CNL_C, TNL, and RNL. Compared to the previous classification, we further subdivided CNLs into three subclasses. The credibility of this classification is supported by phylogenetic analysis and examination of protein domain structures. Importantly, this classification enabled a model to explain the extinction of TNL genes in monocots. Compelling microsynteny evidence underscores this revelation, indicating a clear synteny correspondence between the non-TNLs in monocots and the extinct TNL subclass. Our study provides crucial insights into the genomic origin and divergence of plant NLR subfamilies, unveiling the malleability-driven journey that has shaped the functionality and diversity of plant NLR genes.

A complex regulatory network governs fruit ripening, but natural variations and functional differentiation of fruit ripening genes remain largely unknown. Utilizing a genome-wide association study (GWAS), we identified the NAC family transcription factor MdNAC18.1, whose expression is closely associated with fruit ripening in apple (Malus × domestica Borkh.). MdNAC18.1 activated the transcription of genes related to fruit softening (Polygalacturonase, PG) and ethylene biosynthesis (1-aminocyclopropane-1-carboxylic acid synthase, ACS), thereby promoting fruit ripening of apple and tomato (Solanum lycopersicum). There were two single-nucleotide polymorphisms (SNP-1,545 and SNP-2,002) and a 58-bp insertion–deletion (InDel-58) in the promoter region of MdNAC18.1. Among these, InDel-58 serves as the main effector in activating the expression of MdNAC18.1 and driving fruit ripening. InDel-58 determines the binding affinity of the class D MADS-box protein AGAMOUS-LIKE 11 (MdAGL11), a negative regulator of fruit ripening. The InDel-58 deletion in the early-ripening genotype reduces the inhibitory effect of MdAGL11 on MdNAC18.1. Moreover, MdNAC18.1 and its homologous genes originated from a common ancestor across 61 angiosperms, with functional diversification attributed to tandem replications that occurred in basal angiosperms. In summary, our study revealed how a set of natural variations influence fruit ripening and explored the functional diversification of MdNAC18.1 during evolution.