Whole-genome microsynteny-based phylogeny of angiosperms

Oryza genome evolution through a tetraploid lens

Oryza is a remarkable genus comprising 27 species and 11 genome types, with ~3.4-fold genome size variation, that possesses a virtually untapped reservoir of genes that can be used for crop improvement and neodomestication. Here we present 11 chromosome-level assemblies (nine tetraploid, two diploid) in the context of ~15 million years of evolution and show that the core Oryza (sub)genome is only ~200 Mb and largely syntenic, whereas the remaining nuclear fractions (~80-600 Mb) are intermingled, plastic and rapidly evolving. For the halophyte Oryza coarctata, we found that despite detection of gene fractionation in the subgenomes, homoeologous genes were expressed at higher levels in one subgenome over the other in a mosaic form, demonstrating subgenome equivalence. The integration of these 11 new reference genomes with previously published genome datasets provides a nearly complete view of the consequences of evolution for genome diversification across the genus.

Strain tracking in complex microbiomes using synteny analysis reveals per-species modes of evolution

Microbial species diversify into strains through single-nucleotide mutations and structural changes, such as recombination, insertions and deletions. Most strain-comparison methods quantify differences in single-nucleotide polymorphisms (SNPs) and are insensitive to structural changes. However, recombination is an important driver of phenotypic diversification in many species, including human pathogens. We introduce SynTracker, a tool that compares microbial strains using genome synteny-the order of sequence blocks in homologous genomic regions-in pairs of metagenomic assemblies or genomes. Genome synteny is a rich source of genomic information untapped by current strain-comparison tools. SynTracker has low sensitivity to SNPs, has no database requirement and is robust to sequencing errors. It outperforms existing tools when tracking strains in metagenomic data and is particularly suited for phages, plasmids and other low-data contexts. Applied to single-species datasets and human gut metagenomes, SynTracker, combined with an SNP-based tool, detects strains enriched in either point mutations or structural changes, providing insights into microbial evolution in situ.

Haplotype-resolved and chromosome-level reference genome assembly of Diospyros deyangensis provides insights into the evolution and juvenile growth of persimmon

The Diospyros genus , which includes both wild and cultivated species such as Diospyros lotus and Diospyros kaki, represents a diverse genetic pool with significant agricultural value. In this study, we present a high-quality, haplotype-resolved, chromosome-level genome assembly for Diospyros deyangensis (hereinafter referred to as 'Deyangshi'), an autotetraploid wild species notable for its short juvenile phase, by integrating high-fidelity single-molecule, nanopore sequencing, and high-throughput chromosome conformation capture techniques. The assembled genome size is ~3.01 Gb, anchored onto 60 pseudochromosomes. Comparative genomic analysis revealed that the D. deyangensis genome underwent an additional whole-genome duplication (WGD) event following the eudicots shared ancient hexaploidy event. Resequencing and clustering on 63 samples representing 11 geographically diverse Diospyros accessions revealed significant genetic differentiation between D. deyangensis and D. kaki, as well as between D. kaki and other Diospyros species using population genomic analyses, suggesting that D. kaki followed an independent evolutionary pathway. Additionally, we identified DdELF4 (EARLY FLOWERING 4) from the 'Deyangshi' backcross population using bulked segregant RNA sequencing (BSR-seq) with 50 early-flowering and 50 non-early-flowering individuals. Overexpression of DdELF4 in Arabidopsis resulted in delayed flowering and downregulation of FT gene expression, indicating its role as a flowering repressor. This high-quality genome assembly of 'Deyangshi' provides an essential genomic resource for the Diospyros genus, particularly for breeding programs focused on developing early-flowering persimmon varieties.

Genomic investigation of plant secondary metabolism: insights from synteny network analysis of oxidosqualene cyclase flanking genes

The clustered distribution of genes involved in metabolic pathways within the plant genome has garnered significant attention from researchers. By comparing and analyzing changes in the flanking regions of metabolic genes across a diverse array of species, we can enhance our understanding of the formation and distribution of biosynthetic gene clusters (BGCs). In this study, we have designed a workflow that uncovers and assesses conserved positional relationships between genes in various species by using synteny neighborhood networks (SNN). This workflow is then applied to the analysis of flanking genes associated with oxidosqualene cyclases (OSCs). The method allows for the recognition and comparison of homologous blocks with unique flanking genes accompanying different subfamilies of OSCs. The examination of the flanking genes of OSCs in 122 plant species revealed multiple genes with conserved positional relationships with OSCs in angiosperms. Specifically, the earliest adjacency of OSC genes and CYP716 genes first appeared in basal eudicots, and the nonrandom occurrence of CYP716 genes in the flanking region of OSC persists across different lineages of eudicots. Our study showed the substitution of genes in the flanking region of the OSC varies across different plant lineages, and our approach facilitates the investigation of flanking gene rearrangements in the formation of OSC-related BGCs.

Genome-wide analyses provide insights into genetic variation, phylo- and co-phylogenetic relationships, and biogeography of the entomopathogenic nematode genus Heterorhabditis

Multigene, genus-wide phylogenetic studies have uncovered the limited taxonomic resolution power of commonly used gene markers, particularly of rRNA genes, to discriminate closely related species of the nematode genus Heterorhabditis. In addition, conflicting tree topologies are often obtained using the different gene markers, which limits our understanding of the phylo- and co-phylogenetic relationships and biogeography of the entomopathogenic nematode genus Heterorhabditis. Here we carried out phylogenomic reconstructions using whole nuclear and mitochondrial genomes, and whole ribosomal operon sequences, as well as multiple phylogenetic reconstructions using various single nuclear and mitochondrial genes. Using the inferred phylogenies, we then investigated co-phylogenetic relationships between Heterorhabditis and their Photorhabdus bacterial symbionts and biogeographical patterns. Robust, well-resolved, and highly congruent phylogenetic relationships were reconstructed using both whole nuclear and mitochondrial genomes. Similarly, whole ribosomal operon sequences proved valuable for phylogenomic reconstructions, though they have limited value to discriminate closely related species. In addition, two mitochondrial genes, the cytochrome c oxidase subunit I (cox-1) and the NADH dehydrogenase subunit 4 (nad-4), and two housekeeping genes, the fanconi-associated nuclease 1 (fan-1) and the serine/threonine-protein phosphatase 4 regulatory subunit 1 (ppfr-1), provided the most robust phylogenetic reconstructions compared to other individual genes. According to our findings, whole nuclear and/or mitochondrial genomes are strongly recommended for reconstructing phylogenetic relationships of the genus Heterorhabditis. If whole nuclear and/or mitochondrial genomes are unavailable, a combination of nuclear and mitochondrial genes can be used as an alternative. Under these circumstances, sequences of multiple conspecific isolates in a genus-wide phylogenetic context should be analyzed to avoid artefactual species over-splitting driven by the high intraspecific sequence divergence of mitochondrial genes and to avoid artefactual species lumping driven by the low interspecific sequence divergence of some nuclear genes. On the other hand, we observed that the genera Heterorhabditis and Photorhabdus exhibit diverse biogeographic patterns, ranging from cosmopolitan species to potentially endemic species, and show high phylogenetic congruence, although host switches have also occurred. Our study contributes to a better understanding of the biodiversity and phylo- and co-phylogenetic relationships of an important group of biological control agents and advances our efforts to develop more tools that are compatible with sustainable and eco-friendly agricultural practices.

Deciphering Plant NLR Genomic Evolution: Synteny-Informed Classification Unveils Insights into TNL Gene Loss

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.

Phylogenomic and synteny analysis of BAHD and SCP/SCPL gene families reveal their evolutionary histories in plant specialized metabolism

Plant chemical diversity is largely owing to a number of enzymes which catalyse reactions involved in the assembly, and in the subsequent chemical modifications, of the core structures of major classes of plant specialized metabolites. One such reaction is acylation. With this in mind, to study the deep evolutionary history of BAHD and the serine-carboxypeptidase-like (SCPL) acyltransferase genes, we assembled phylogenomic synteny networks based on a large-scale inference analysis of orthologues across whole-genome sequences of 126 species spanning Stramenopiles and Archaeplastida, including Arabidopsis thaliana, tomato (Solanum lycopersicum) and maize (Zea mays). As such, this study combined the study of genomic location with changes in gene sequences. Our analyses revealed that serine-carboxypeptidase (SCP)/serine-carboxypeptidase-like (SCPL) genes had a deeper evolutionary origin than BAHD genes, which expanded massively on the transition to land and with the development of the vascular system. The two gene families additionally display quite distinct patterns of copy number variation across phylogenies as well as differences in cross-phylogenetic syntenic network components. In unlocking the above observations, our analyses demonstrate the possibilities afforded by modern phylogenomic (syntenic) networks, but also highlight their current limitations, as demonstrated by the inability of phylogenetic methods to separate authentic SCPL acyltransferases from standard SCP peptide hydrolases.This article is part of the theme issue 'The evolution of plant metabolism'.

Evolutionary dynamics and functional characterization of proximal duplicated sorbitol-6-phosphate dehydrogenase genes in Rosaceae

Sorbitol is a critical photosynthate and storage substance in the Rosaceae family. Sorbitol 6-phosphate dehydrogenase (S6PDH) functions as the pivotal rate-limiting enzyme in sorbitol synthesis. The origin and functional diversification of S6PDH in Rosaceae remain unclear, largely due to the complicated interplay of gene duplications. Here, we investigated the synteny relationships among all identified S6PDH genes in representative genomes within the Rosaceae family. By integrating phylogenetic analyses, we elucidated the lineage-specific expansion and syntenic conservation of S6PDH across diverse Rosaceae plant lineages. We found that S6PDH can be traced back to a pair of proximal duplicated genes of the common ancestor of the Rosaceae, and the further amplification of S6PDH in the Maleae primarily relies on WGD events in their lineages. In Rosaceae species, multiple copies of the S6PDH gene are preliminarily divided into two main clades (Clade 1 and Clade 2) based on sequence similarity. These clades have evolved to acquire different functional directions. In Clade 1, lineage-specific transposition events in the Amygdaloideae have led to changes in gene expression patterns and promoted lineage evolution. This is mainly characterized by a decrease in enzymatic activity and transcriptional expression in the leaves, but also includes specific functional diversification, such as sustained post-harvest fruit expression and enhanced expression under biotic stress in certain tissues. In contrast, S6PDH in the Rosoideae and Dryadoideae has not undergone additional duplications beyond early proximal duplication. The loss of exons and variations in exon length might the key factor leading to reduced enzymatic activity in the Clade 2 proximal gene pairs. Collectively, our findings illuminate the dynamic nature of S6PDH evolution and reveal the intricate interplay between duplication, transposition, and functional diversification. This work not only contributes valuable insights into the genetic mechanisms underlying sorbitol metabolism but also establish a crucial foundation for future investigations aimed at comprehensively characterizing the variations of sorbitol metabolism across different subfamilies within the Rosaceae family.

A molecular phylogeny of the Petaluridae (Odonata: Anisoptera): A 160-Million-Year-Old story of drift and extinction

Petaluridae (Odonata: Anisoptera) is a relict dragonfly family, having diverged from its sister family in the Jurassic, of eleven species that are notable among odonates (dragonflies and damselflies) for their exclusive use of fen and bog habitats, their burrowing behavior as nymphs, large body size as adults, and extended lifespans. To date, several nodes within this family remain unresolved, limiting the study of the evolution of this peculiar family. Using an anchored hybrid enrichment dataset of over 900 loci we reconstructed the species tree of Petaluridae. To estimate the temporal origin of the genera within this family, we used a set of well-vetted fossils and a relaxed molecular clock model in a divergence time estimation analysis. We estimate that Petaluridae originated in the early Cretaceous and confirm the existence of monophyletic Gondwanan and Laurasian clades within the family. Our relaxed molecular clock analysis estimated that these clades diverged from their MRCA approximately 160 mya. Extant lineages within this family were identified to have persisted from 6 (Uropetala) to 120 million years (Phenes). Our biogeographical analyses focusing on a set of key regions suggest that divergence within Petaluridae is largely correlated with continental drift, the exposure of land bridges, and the development of mountain ranges. Our results support the hypothesis that species within Petaluridae have persisted for tens of millions of years, with little fossil evidence to suggest widespread extinction in the family, despite optimal conditions for the fossilization of nymphs. Petaluridae appear to be a rare example of habitat specialists that have persisted for tens of millions of years.

Genome-Wide Transcriptional Response of Avocado to Fusarium sp. Infection

The avocado crop is relevant for its economic importance and because of its unique evolutionary history. However, there is a lack of information regarding the molecular processes during the defense response against fungal pathogens. Therefore, using a genome-wide approach in this work, we investigated the transcriptional response of the Mexican horticultural race of avocado (Persea americana var. drymifolia), including miRNAs profile and their possible targets. For that, we established an avocado-Fusarium hydroponic pathosystem and studied the response for 21 days. To guarantee robustness in the analysis, first, we improved the avocado genome assembly available for this variety, resulting in 822.49 Mbp in length with 36,200 gene models. Then, using an RNA-seq approach, we identified 13,778 genes differentially expressed in response to the Fusarium infection. According to their expression profile across time, these genes can be clustered into six groups, each associated with specific biological processes. Regarding non-coding RNAs, 8 of the 57 mature miRNAs identified in the avocado genome are responsive to infection caused by Fusarium, and the analysis revealed a total of 569 target genes whose transcript could be post-transcriptionally regulated. This study represents the first research in avocados to comprehensively explore the role of miRNAs in orchestrating defense responses against Fusarium spp. Also, this work provides valuable data about the genes involved in the intricate response of the avocado during fungal infection.

CLOCI: unveiling cryptic fungal gene clusters with generalized detection

Gene clusters are genomic loci that contain multiple genes that are functionally and genetically linked. Gene clusters collectively encode diverse functions, including small molecule biosynthesis, nutrient assimilation, metabolite degradation, and production of proteins essential for growth and development. Identifying gene clusters is a powerful tool for small molecule discovery and provides insight into the ecology and evolution of organisms. Current detection algorithms focus on canonical 'core' biosynthetic functions many gene clusters encode, while overlooking uncommon or unknown cluster classes. These overlooked clusters are a potential source of novel natural products and comprise an untold portion of overall gene cluster repertoires. Unbiased, function-agnostic detection algorithms therefore provide an opportunity to reveal novel classes of gene clusters and more precisely define genome organization. We present CLOCI (Co-occurrence Locus and Orthologous Cluster Identifier), an algorithm that identifies gene clusters using multiple proxies of selection for coordinated gene evolution. Our approach generalizes gene cluster detection and gene cluster family circumscription, improves detection of multiple known functional classes, and unveils non-canonical gene clusters. CLOCI is suitable for genome-enabled small molecule mining, and presents an easily tunable approach for delineating gene cluster families and homologous loci.

Annelid Comparative Genomics and the Evolution of Massive Lineage-Specific Genome Rearrangement in Bilaterians

The organization of genomes into chromosomes is critical for processes such as genetic recombination, environmental adaptation, and speciation. All animals with bilateral symmetry inherited a genome structure from their last common ancestor that has been highly conserved in some taxa but seemingly unconstrained in others. However, the evolutionary forces driving these differences and the processes by which they emerge have remained largely uncharacterized. Here, we analyze genome organization across the phylum Annelida using 23 chromosome-level annelid genomes. We find that while many annelid lineages have maintained the conserved bilaterian genome structure, the Clitellata, a group containing leeches and earthworms, possesses completely scrambled genomes. We develop a rearrangement index to quantify the extent of genome structure evolution and show that, compared to the last common ancestor of bilaterians, leeches and earthworms have among the most highly rearranged genomes of any currently sampled species. We further show that bilaterian genomes can be classified into two distinct categories-high and low rearrangement-largely influenced by the presence or absence, respectively, of chromosome fission events. Our findings demonstrate that animal genome structure can be highly variable within a phylum and reveal that genome rearrangement can occur both in a gradual, stepwise fashion, or rapid, all-encompassing changes over short evolutionary timescales.

Topological change of soil microbiota networks for forest resilience under global warming

Forest management by thinning can mitigate the detrimental impact of increasing drought caused by global warming. Growing evidence shows that the soil microbiota can coordinate the dynamic relationship between forest functions and drought intensity, but how they function as a cohesive whole remains elusive. We outline a statistical topology model to chart the roadmap of how each microbe acts and interacts with every other microbe to shape the dynamic changes of microbial communities under forest management. To demonstrate its utility, we analyze a soil microbiota data collected from a two-way longitudinal factorial experiment involving three stand densities and three levels of rainfall over a growing season in artificial plantations of a forest tree - larix (Larix kaempferi). We reconstruct the most sophisticated soil microbiota networks that code maximally informative microbial interactions and trace their dynamic trajectories across time, space, and environmental signals. By integrating GLMY homology theory, we dissect the topological architecture of these so-called omnidirectional networks and identify key microbial interaction pathways that play a pivotal role in mediating the structure and function of soil microbial communities. The statistical topological model described provides a systems tool for studying how microbial community assembly alters its structure, function and evolution under climate change.

The Spiral Model of Evolution: Stable Life Forms of Organisms and Unstable Life Forms of Cancers

If one must prioritize among the vast array of contributing factors to cancer evolution, environmental-stress-mediated chromosome instability (CIN) should easily surpass individual gene mutations. CIN leads to the emergence of genomically unstable life forms, enabling them to grow dominantly within the stable life form of the host. In contrast, stochastic gene mutations play a role in aiding the growth of the cancer population, with their importance depending on the initial emergence of the new system. Furthermore, many specific gene mutations among the many available can perform this function, decreasing the clinical value of any specific gene mutation. Since these unstable life forms can respond to treatment differently than stable ones, cancer often escapes from drug treatment by forming new systems, which leads to problems during the treatment for patients. To understand how diverse factors impact CIN-mediated macroevolution and genome integrity-ensured microevolution, the concept of two-phased cancer evolution is used to reconcile some major characteristics of cancer, such as bioenergetic, unicellular, and multicellular evolution. Specifically, the spiral of life function model is proposed, which integrates major historical evolutionary innovations and conservation with information management. Unlike normal organismal evolution in the microevolutionary phase, where a given species occupies a specific location within the spiral, cancer populations are highly heterogenous at multiple levels, including epigenetic levels. Individual cells occupy different levels and positions within the spiral, leading to supersystems of mixed cellular populations that exhibit both macro and microevolution. This analysis, utilizing karyotype to define the genetic networks of the cellular system and CIN to determine the instability of the system, as well as considering gene mutation and epigenetics as modifiers of the system for information amplification and usage, explores the high evolutionary potential of cancer. It provides a new, unified understanding of cancer as a supersystem, encouraging efforts to leverage the dynamics of CIN to develop improved treatment options. Moreover, it offers a historically contingent model for organismal evolution that reconciles the roles of both evolutionary innovation and conservation through macroevolution and microevolution, respectively.

Global identification of LIM genes in response to different heat stress regimes in Lactuca sativa

BACKGROUND

LIM (Lineage-11 (LIN-11), Insulin-1 (ISL-1), and Mechanotransduction-3 (MEC-3)) genes belong to a family that hold ubiquitous properties contributing to organ, seed, and pollen development as well as developmental and cellular responses to biotic and abiotic stresses. Lettuce (Lactuca sativa) is a highly consumed vegetable crop susceptible heat stress. High temperatures limit lettuce's overall yield, quality and marketability. Lettuce LIM genes have not been identified and their role in response to high temperatures is not known. Aiming to identify potential new targets for thermoresilience, we searched for LIM genes in lettuce and compared them with orthologous of several dicotyledons and monocotyledons plant species.

RESULTS

We identified fourteen lettuce LIM genes distributed into eight different subgroups using a genome-wide analysis strategy. Three belonging to DAR (DA means "large" in Chinese) class I, two DAR class II, one in the WLIM1, two in the WLIM2, one in the PLIM1, two in the PLIM2 class, one ßLIM and two δLIMs. No DAR-like were identified in any of the species analyzed including lettuce. Interestingly, unlike other gene families in lettuce which underwent large genome tandem duplications, LIM genes did not increase in number compared to other plant species. The response to heat stress induced a dynamic transcriptional response on LsLIM genes. All heat stress regimes, including night stress, day stress and day and night stress were largely responsible for changes in LIM transcriptional expression.

CONCLUSIONS

Our global analysis at the genome level provides a detailed identification of LIM genes in lettuce and other dicotyledonous and monocotyledonous plant species. Gene structure, physical and chemical properties as well as chromosomal location and Cis-regulatory element analysis together with our gene expression analysis under different temperature regimes identified LsWLIM1, LsWLIM2b, LsDAR3 and LsDAR5 as candidate genes that could be used by breeding programs aiming to produce lettuce varieties able to withstand high temperatures.

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

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

Doubling down on polyploid discoveries: Global advances in genomics and ecological impacts of polyploidy

All flowering plants are now recognized as diploidized paleopolyploids (Jiao et al., 2011; One Thousand Plant Transcriptomes Initiative, 2019), and polyploid species comprise approximately 30% of contemporary plant species (Wood et al., 2009; Barker et al., 2016a). A major implication of these discoveries is that, to appreciate the evolution of plant diversity, we need to understand the fundamental biology of polyploids and diploidization. This need is broadly recognized by our community as there is a continued, growing interest in polyploidy as a research topic. Over the past 25 years, the sequencing and analysis of plant genomes has revolutionized our understanding of the importance of polyploid speciation to the evolution of land plants.

Species-tree topology impacts the inference of ancient whole-genome duplications across the angiosperm phylogeny

PREMISE

The history of angiosperms is marked by repeated rounds of ancient whole-genome duplications (WGDs). Here we used state-of-the-art methods to provide an up-to-date view of the distribution of WGDs in the history of angiosperms that considers both uncertainty introduced by different WGD inference methods and different underlying species-tree hypotheses.

METHODS

We used the distribution synonymous divergences (Ks) of paralogs and orthologs from transcriptomic and genomic data to infer and place WGDs across two hypothesized angiosperm phylogenies. We further tested these WGD hypotheses with syntenic inferences and Bayesian models of duplicate gene gain and loss.

RESULTS

The predicted number of WGDs in the history of angiosperms (~170) based on the current taxon sampling is largely similar across different inference methods, but varies in the precise placement of WGDs on the phylogeny. Ks-based methods often yield alternative hypothesized WGD placements due to variation in substitution rates among lineages. Phylogenetic models of duplicate gene gain and loss are more robust to topological variation. However, errors in species-tree inference can still produce spurious WGD hypotheses, regardless of method used.

CONCLUSIONS

Here we showed that different WGD inference methods largely agree on an average of 3.5 WGD in the history of individual angiosperm species. However, the precise placement of WGDs on the phylogeny is subject to the WGD inference method and tree topology. As researchers continue to test hypotheses regarding the impacts ancient WGDs have on angiosperm evolution, it is important to consider the uncertainty of the phylogeny as well as WGD inference methods.

Evolution of Whirly1 in the angiosperms: sequence, splicing, and expression in a clade of early transitional mycoheterotrophic orchids

The plastid-targeted transcription factor Whirly1 (WHY1) has been implicated in chloroplast biogenesis, plastid genome stability, and fungal defense response, which together represent characteristics of interest for the study of autotrophic losses across the angiosperms. While gene loss in the plastid and nuclear genomes has been well studied in mycoheterotrophic plants, the evolution of the molecular mechanisms impacting genome stability is completely unknown. Here, we characterize the evolution of WHY1 in four early transitional mycoheterotrophic orchid species in the genus Corallorhiza by synthesizing the results of phylogenetic, transcriptomic, and comparative genomic analyses with WHY1 genomic sequences sampled from 21 orders of angiosperms. We found an increased number of non-canonical WHY1 isoforms assembled from all but the greenest Corallorhiza species, including intron retention in some isoforms. Within Corallorhiza, phylotranscriptomic analyses revealed the presence of tissue-specific differential expression of WHY1 in only the most photosynthetically capable species and a coincident increase in the number of non-canonical WHY1 isoforms assembled from fully mycoheterotrophic species. Gene- and codon-level tests of WHY1 selective regimes did not infer significant signal of either relaxed selection or episodic diversifying selection in Corallorhiza but did so for relaxed selection in the late-stage full mycoheterotrophic orchids Epipogium aphyllum and Gastrodia elata. Additionally, nucleotide substitutions that most likely impact the function of WHY1, such as nonsense mutations, were only observed in late-stage mycoheterotrophs. We propose that our findings suggest that splicing and expression changes may precede the selective shifts we inferred for late-stage mycoheterotrophic species, which therefore does not support a primary role for WHY1 in the transition to mycoheterotrophy in the Orchidaceae. Taken together, this study provides the most comprehensive view of WHY1 evolution across the angiosperms to date.

Pan-evolutionary and regulatory genome architecture delineated by an integrated macro- and microsynteny approach

The forthcoming massive genome data generated by the Earth BioGenome Project will open up a new era of comparative genomics, for which genome synteny analysis provides an important framework. Profiling genome synteny represents an essential step in elucidating genome architecture, regulatory blocks/elements and their evolutionary history. Here we describe PanSyn, ( https://github.com/yhw320/PanSyn ), the most comprehensive and up-to-date genome synteny pipeline, providing step-by-step instructions and application examples to demonstrate its usage. PanSyn inherits both basic and advanced functions from existing popular tools, offering a user-friendly, highly customized approach for genome macrosynteny analysis and integrated pan-evolutionary and regulatory analysis of genome architecture, which are not yet available in public synteny software or tools. The advantages of PanSyn include: (i) advanced microsynteny analysis by functional profiling of microsynteny genes and associated regulatory elements; (ii) comprehensive macrosynteny analysis, including the inference of karyotype evolution from ancestors to extant species; and (iii) functional integration of microsynteny and macrosynteny for pan-evolutionary profiling of genome architecture and regulatory blocks, as well as integration with external functional genomics datasets from three- or four-dimensional genome and ENCODE projects. PanSyn requires basic knowledge of the Linux environment and Perl programming language and the ability to access a computer cluster, especially for large-scale genomic comparisons. Our protocol can be easily implemented by a competent graduate student or postdoc and takes several days to weeks to execute for dozens to hundreds of genomes. PanSyn provides yet the most comprehensive and powerful tool for integrated evolutionary and functional genomics.

Detection of Ghost Introgression Requires Exploiting Topological and Branch Length Information

In recent years, the study of hybridization and introgression has made significant progress, with ghost introgression-the transfer of genetic material from extinct or unsampled lineages to extant species-emerging as a key area for research. Accurately identifying ghost introgression, however, presents a challenge. To address this issue, we focused on simple cases involving 3 species with a known phylogenetic tree. Using mathematical analyses and simulations, we evaluated the performance of popular phylogenetic methods, including HyDe and PhyloNet/MPL, and the full-likelihood method, Bayesian Phylogenetics and Phylogeography (BPP), in detecting ghost introgression. Our findings suggest that heuristic approaches relying on site-pattern counts or gene-tree topologies struggle to differentiate ghost introgression from introgression between sampled non-sister species, frequently leading to incorrect identification of donor and recipient species. The full-likelihood method BPP uses multilocus sequence alignments directly-hence taking into account both gene-tree topologies and branch lengths, by contrast, is capable of detecting ghost introgression in phylogenomic datasets. We analyzed a real-world phylogenomic dataset of 14 species of Jaltomata (Solanaceae) to showcase the potential of full-likelihood methods for accurate inference of introgression.

Common Methods for Phylogenetic Tree Construction and Their Implementation in R

A phylogenetic tree can reflect the evolutionary relationships between species or gene families, and they play a critical role in modern biological research. In this review, we summarize common methods for constructing phylogenetic trees, including distance methods, maximum parsimony, maximum likelihood, Bayesian inference, and tree-integration methods (supermatrix and supertree). Here we discuss the advantages, shortcomings, and applications of each method and offer relevant codes to construct phylogenetic trees from molecular data using packages and algorithms in R. This review aims to provide comprehensive guidance and reference for researchers seeking to construct phylogenetic trees while also promoting further development and innovation in this field. By offering a clear and concise overview of the different methods available, we hope to enable researchers to select the most appropriate approach for their specific research questions and datasets.

The promise and pitfalls of synteny in phylogenomics

Reconstructing the tree of life remains a central goal in biology. Early methods, which relied on small numbers of morphological or genetic characters, often yielded conflicting evolutionary histories, undermining confidence in the results. Investigations based on phylogenomics, which use hundreds to thousands of loci for phylogenetic inquiry, have provided a clearer picture of life's history, but certain branches remain problematic. To resolve difficult nodes on the tree of life, 2 recent studies tested the utility of synteny, the conserved collinearity of orthologous genetic loci in 2 or more organisms, for phylogenetics. Synteny exhibits compelling phylogenomic potential while also raising new challenges. This Essay identifies and discusses specific opportunities and challenges that bear on the value of synteny data and other rare genomic changes for phylogenomic studies. Synteny-based analyses of highly contiguous genome assemblies mark a new chapter in the phylogenomic era and the quest to reconstruct the tree of life.

Solving an enigma in the tree of life, at the origins of teleost fishes

Tracing the phylogenetic relationships between species is one of the fundamental objectives of evolutionary biology. Since Charles Darwin's seminal work in the 19th century, considerable progress has been made towards establishing a tree of life that summarises the evolutionary history of species. Nevertheless, substantial uncertainties still remain. Specifically, the relationships at the origins of teleost fishes have been the subject of extensive debate over the last 50 years. This question has major implications for various research fields: there are almost 30,000 species in the teleost group, which includes invaluable model organisms for biomedical, evolutionary and ecological studies. Here, we present the work in which we solved this enigma. We demonstrated that eels are more closely related to bony-tongued fishes than to the rest of teleost fishes. We achieved this by taking advantage of new genomic data and leveraging innovative phylogenetic markers. Notably, in addition to traditional molecular phylogeny methods based on the evolution of gene sequences, we also considered the evolution of gene order along the DNA molecule. We discuss the challenges and opportunities that these new markers represent for the field of molecular phylogeny, and in particular the possibilities they offer for re-examining other controversial branches in the tree of life.

Angiosperm-wide analysis of fruit and ovary evolution aided by a new nuclear phylogeny supports association of the same ovary type with both dry and fleshy fruits

Fruit functions in seed protection and dispersal and belongs to many dry and fleshy types, yet their evolutionary pattern remains unclear in part due to uncertainties in the phylogenetic relationships among several orders and families. Thus we used nuclear genes of 502 angiosperm species representing 231 families to reconstruct a well supported phylogeny, with resolved relationships for orders and families with previously uncertain placements. Using this phylogeny as a framework, molecular dating supports a Triassic origin of the crown angiosperms, followed by the emergence of most orders in the Jurassic and Cretaceous and their rise to ecological dominance during the Cretaceous Terrestrial Revolution. The robust phylogeny allowed an examination of the evolutionary pattern of fruit and ovary types, revealing a trend of parallel carpel fusions during early diversifications in eudicots, monocots, and magnoliids. Moreover, taxa in the same order or family with the same ovary type can develop either dry or fleshy fruits with strong correlations between specific types of dry and fleshy fruits; such associations of ovary, dry and fleshy fruits define several ovary-fruit "modules" each found in multiple families. One of the frequent modules has an ovary containing multiple ovules, capsules and berries, and another with an ovary having one or two ovules, achenes (or other single-seeded dry fruits) and drupes. This new perspective of relationships among fruit types highlights the closeness of specific dry and fleshy fruit types, such as capsule and berry, that develop from the same ovary type and belong to the same module relative to dry and fleshy fruits of other modules (such as achenes and drupes). Further analyses of gene families containing known genes for ovary and fruit development identified phylogenetic nodes with multiple gene duplications, supporting a possible role of whole-genome duplications, in combination with climate changes and animal behaviors, in angiosperm fruit and ovary diversification.

The chromosome-level genome of the submerged plant Cryptocoryne crispatula provides insights into the terrestrial-freshwater transition in Araceae

Plant terrestrialization (i.e. the transition to a terrestrial environment) is a significant evolutionary event that has been intensively studied. While certain plant lineages, particularly in angiosperms, have re-adapted to freshwater habitats after colonizing terrene, however, the molecular mechanism of the terrestrial-freshwater (T-F) transition remains limited. Here, the basal monocot Araceae was selected as the study object to explore the T-F transition adaptation mechanism by comparative genomic analysis. Our findings revealed that the substitution rates significantly increased in the lineage of freshwater Araceae, which may promote their adaptation to the freshwater habitat. Additionally, 20 gene sets across all four freshwater species displayed signs of positive selection contributing to tissue development and defense responses in freshwater plants. Comparative synteny analysis showed that genes specific to submerged plants were enriched in cellular respiration and photosynthesis. In contrast, floating plants were involved in regulating gene expression, suggesting that gene and genome duplications may provide the original material for plants to adapt to the freshwater environment. Our study provides valuable insights into the genomic aspects of the transition from terrestrial to aquatic environments in Araceae, laying the groundwork for future research in the angiosperm.

The Genome of Plasmodium gonderi: Insights into the Evolution of Human Malaria Parasites

Plasmodium species causing malaria in humans are not monophyletic, sharing common ancestors with nonhuman primate parasites. Plasmodium gonderi is one of the few known Plasmodium species infecting African old-world monkeys that are not found in apes. This study reports a de novo assembled P. gonderi genome with complete chromosomes. The P. gonderi genome shares codon usage, syntenic blocks, and other characteristics with the human parasites Plasmodium ovale s.l. and Plasmodium malariae, also of African origin, and the human parasite Plasmodium vivax and species found in nonhuman primates from Southeast Asia. Using phylogenetically aware methods, newly identified syntenic blocks were found enriched with conserved metabolic genes. Regions outside those blocks harbored genes encoding proteins involved in the vertebrate host-Plasmodium relationship undergoing faster evolution. Such genome architecture may have facilitated colonizing vertebrate hosts. Phylogenomic analyses estimated the common ancestor between P. vivax and an African ape parasite P. vivax-like, within the Asian nonhuman primates parasites clade. Time estimates incorporating P. gonderi placed the P. vivax and P. vivax-like common ancestor in the late Pleistocene, a time of active migration of hominids between Africa and Asia. Thus, phylogenomic and time-tree analyses are consistent with an Asian origin for P. vivax and an introduction of P. vivax-like into Africa. Unlike other studies, time estimates for the clade with Plasmodium falciparum, the most lethal human malaria parasite, coincide with their host species radiation, African hominids. Overall, the newly assembled genome presented here has the quality to support comparative genomic investigations in Plasmodium.

Gene Transfer-Based Phylogenetics: Analytical Expressions and Additivity via Birth-Death Theory

The genomic era has opened up vast opportunities in molecular systematics, one of which is deciphering the evolutionary history in fine detail. Under this mass of data, analyzing the point mutations of standard markers is often too crude and slow for fine-scale phylogenetics. Nevertheless, genome dynamics (GD) events provide alternative, often richer information. The synteny index (SI) between a pair of genomes combines gene order and gene content information, allowing the comparison of genomes of unequal gene content, together with order considerations of their common genes. Recently, genome dynamics has been modeled as a continuous-time Markov process, and gene distance in the genome as a birth-death-immigration process. Nevertheless, due to complexities arising in this setting, no precise and provably consistent estimators could be derived, resulting in heuristic solutions. Here, we extend this modeling approach by using techniques from birth-death theory to derive explicit expressions of the system's probabilistic dynamics in the form of rational functions of the model parameters. This, in turn, allows us to infer analytically accurate distances between organisms based on their SI. Subsequently, we establish additivity of this estimated evolutionary distance (a desirable property yielding phylogenetic consistency). Applying the new measure in simulation studies shows that it provides accurate results in realistic settings and even under model extensions such as gene gain/loss or over a tree structure. In the real-data realm, we applied the new formulation to unique data structure that we constructed-the ordered orthology DB-based on a new version of the EggNOG database, to construct a tree with more than 4.5K taxa. To the best of our knowledge, this is the largest gene-order-based tree constructed and it overcomes shortcomings found in previous approaches. Constructing a GD-based tree allows to confirm and contrast findings based on other phylogenetic approaches, as we show.

Genome-wide expansion and reorganization during grass evolution: from 30 Mb chromosomes in rice and Brachypodium to 550 Mb in Avena

BACKGROUND

The BOP (Bambusoideae, Oryzoideae, and Pooideae) clade of the Poaceae has a common ancestor, with similarities to the genomes of rice, Oryza sativa (2n = 24; genome size 389 Mb) and Brachypodium, Brachypodium distachyon (2n = 10; 271 Mb). We exploit chromosome-scale genome assemblies to show the nature of genomic expansion, structural variation, and chromosomal rearrangements from rice and Brachypodium, to diploids in the tribe Aveneae (e.g., Avena longiglumis, 2n = 2x = 14; 3,961 Mb assembled to 3,850 Mb in chromosomes).

RESULTS

Most of the Avena chromosome arms show relatively uniform expansion over the 10-fold to 15-fold genome-size increase. Apart from non-coding sequence diversification and accumulation around the centromeres, blocks of genes are not interspersed with blocks of repeats, even in subterminal regions. As in the tribe Triticeae, blocks of conserved synteny are seen between the analyzed species with chromosome fusion, fission, and nesting (insertion) events showing deep evolutionary conservation of chromosome structure during genomic expansion. Unexpectedly, the terminal gene-rich chromosomal segments (representing about 50 Mb) show translocations between chromosomes during speciation, with homogenization of genome-specific repetitive elements within the tribe Aveneae. Newly-formed intergenomic translocations of similar extent are found in the hexaploid A. sativa.

CONCLUSIONS

The study provides insight into evolutionary mechanisms and speciation in the BOP clade, which is valuable for measurement of biodiversity, development of a clade-wide pangenome, and exploitation of genomic diversity through breeding programs in Poaceae.

Spiralian genomics and the evolution of animal genome architecture

Recent developments in sequencing technologies have greatly improved our knowledge of phylogenetic relationships and genomic architectures throughout the tree of life. Spiralia, a diverse clade within Protostomia, is essential for understanding the evolutionary history of parasitism, gene conversion, nervous systems and animal body plans. In this review, we focus on the current hypotheses of spiralian phylogeny and investigate the impact of long-read sequencing on the quality of genome assemblies. We examine chromosome-level assemblies to highlight key genomic features that have driven spiralian evolution, including karyotype, synteny and the Hox gene organization. In addition, we show how chromosome rearrangement has influenced spiralian genomic structures. Although spiralian genomes have undergone substantial changes, they exhibit both conserved and lineage-specific features. We recommend increasing sequencing efforts and expanding functional genomics research to deepen insights into spiralian biology.

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.

Plant Tolerance to Drought Stress with Emphasis on Wheat

Environmental stresses, such as drought, have negative effects on crop yield. Drought is a stress whose impact tends to increase in some critical regions. However, the worldwide population is continuously increasing and climate change may affect its food supply in the upcoming years. Therefore, there is an ongoing effort to understand the molecular processes that may contribute to improving drought tolerance of strategic crops. These investigations should contribute to delivering drought-tolerant cultivars by selective breeding. For this reason, it is worthwhile to review regularly the literature concerning the molecular mechanisms and technologies that could facilitate gene pyramiding for drought tolerance. This review summarizes achievements obtained using QTL mapping, genomics, synteny, epigenetics, and transgenics for the selective breeding of drought-tolerant wheat cultivars. Synthetic apomixis combined with the msh1 mutation opens the way to induce and stabilize epigenomes in crops, which offers the potential of accelerating selective breeding for drought tolerance in arid and semi-arid regions.

Unique gene duplications and conserved microsynteny potentially associated with resistance to wood decay in the Lauraceae

Wood decay resistance (WDR) is marking the value of wood utilization. Many trees of the Lauraceae have exceptional WDR, as evidenced by their use in ancient royal palace buildings in China. However, the genetics of WDR remain elusive. Here, through comparative genomics, we revealed the unique characteristics related to the high WDR in Lauraceae trees. We present a 1.27-Gb chromosome-level assembly for Lindera megaphylla (Lauraceae). Comparative genomics integrating major groups of angiosperm revealed Lauraceae species have extensively shared gene microsynteny associated with the biosynthesis of specialized metabolites such as isoquinoline alkaloids, flavonoid, lignins and terpenoid, which play significant roles in WDR. In Lauraceae genomes, tandem and proximal duplications (TD/PD) significantly expanded the coding space of key enzymes of biosynthesis pathways related to WDR, which may enhance the decay resistance of wood by increasing the accumulation of these compounds. Among Lauraceae species, genes of WDR-related biosynthesis pathways showed remarkable expansion by TD/PD and conveyed unique and conserved motifs in their promoter and protein sequences, suggesting conserved gene collinearity, gene expansion and gene regulation supporting the high WDR. Our study thus reveals genomic profiles related to biochemical transitions among major plant groups and the genomic basis of WDR in the Lauraceae.

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.

Reconstruction of hundreds of reference ancestral genomes across the eukaryotic kingdom

Ancestral sequence reconstruction is a fundamental aspect of molecular evolution studies and can trace small-scale sequence modifications through the evolution of genomes and species. In contrast, fine-grained reconstructions of ancestral genome organizations are still in their infancy, limiting our ability to draw comprehensive views of genome and karyotype evolution. Here we reconstruct the detailed gene contents and organizations of 624 ancestral vertebrate, plant, fungi, metazoan and protist genomes, 183 of which are near-complete chromosomal gene order reconstructions. Reconstructed ancestral genomes are similar to their descendants in terms of gene content as expected and agree precisely with reference cytogenetic and in silico reconstructions when available. By comparing successive ancestral genomes along the phylogenetic tree, we estimate the intra- and interchromosomal rearrangement history of all major vertebrate clades at high resolution. This freely available resource introduces the possibility to follow evolutionary processes at genomic scales in chronological order, across multiple clades and without relying on a single extant species as reference.

The Saururus chinensis genome provides insights into the evolution of pollination strategies and herbaceousness in magnoliids

Saururus chinensis, an herbaceous magnoliid without perianth, represents a clade of early-diverging angiosperms that have gone through woodiness-herbaceousness transition and pollination obstacles: the characteristic white leaves underneath inflorescence during flowering time are considered a substitute for perianth to attract insect pollinators. Here, using the newly sequenced S. chinensis genome, we revisited the phylogenetic position of magnoliids within mesangiosperms, and recovered a sister relationship for magnoliids and Chloranthales. By considering differentially expressed genes, we identified candidate genes that are involved in the morphogenesis of the white leaves in S. chinensis. Among those genes, we verified - in a transgenic experiment with Arabidopsis - that increasing the expression of the "pseudo-etiolation in light" gene (ScPEL) can inhibit the biosynthesis of chlorophyll. ScPEL is thus likely responsible for the switches between green and white leaves, suggesting that changes in gene expression may underlie the evolution of pollination strategies. Despite being an herbaceous plant, S. chinensis still has vascular cambium and maintains the potential for secondary growth as a woody plant, because the necessary machinery, i.e., the entire gene set involved in lignin biosynthesis, is well preserved. However, similar expression levels of two key genes (CCR and CAD) between the stem and other tissues in the lignin biosynthesis pathway are possibly associated with the herbaceous nature of S. chinensis. In conclusion, the S. chinensis genome provides valuable insights into the adaptive evolution of pollination in Saururaceae and reveals a possible mechanism for the evolution of herbaceousness in magnoliids.

Phylogenetic studies of magnoliids: Advances and perspectives

Magnoliids are the largest flowering plant clades outside of the eudicots and monocots, which are distributed worldwide and have high economic, ornamental and ecological values. Eudicots, monocots and magnoliids are the three major clades of Mesangiospermae, and their phylogenetic relationship is one of the most interesting issues. In recent years, with the continuous accumulation of genomic information, the evolutionary status of magnoliids has become a hot spot in plant phylogenetic research. Although great efforts have been made to study the evolution of magnoliids using molecular data from several representative species such as nuclear genome, plastid genome, mitochondrial genome, and transcriptome, the results of current studies on the phylogenetic status of magnoliids are inconsistent. Here, we systematically describe the current understanding of the molecular research on magnoliid phylogeny and review the differences in the evolutionary state of magnoliids. Understanding the research approaches and limitations of magnoliid phylogeny can guide research strategies to further improve the study of the phylogenetic evolution of magnoliids.

Variation in plant Toll/Interleukin-1 receptor domain protein dependence on ENHANCED DISEASE SUSCEPTIBILITY 1

Toll/Interleukin-1 receptor (TIR) domains are integral to immune systems across all kingdoms. In plants, TIRs are present in nucleotide-binding leucine-rich repeat (NLR) immune receptors, NLR-like, and TIR-only proteins. Although TIR-NLR and TIR signaling in plants require the ENHANCED DISEASE SUSCEPTIBILITY 1 (EDS1) protein family, TIRs persist in species that have no EDS1 members. To assess whether particular TIR groups evolved with EDS1, we searched for TIR-EDS1 co-occurrence patterns. Using a large-scale phylogenetic analysis of TIR domains from 39 algal and land plant species, we identified 4 TIR families that are shared by several plant orders. One group occurred in TIR-NLRs of eudicots and another in TIR-NLRs across eudicots and magnoliids. Two further groups were more widespread. A conserved TIR-only group co-occurred with EDS1 and members of this group elicit EDS1-dependent cell death. In contrast, a maize (Zea mays) representative of TIR proteins with tetratricopeptide repeats was also present in species without EDS1 and induced EDS1-independent cell death. Our data provide a phylogeny-based plant TIR classification and identify TIRs that appear to have evolved with and are dependent on EDS1, while others have EDS1-independent activity.

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.

Phylogenomic resolution of order- and family-level monocot relationships using 602 single-copy nuclear genes and 1375 BUSCO genes

We assess relationships among 192 species in all 12 monocot orders and 72 of 77 families, using 602 conserved single-copy (CSC) genes and 1375 benchmarking single-copy ortholog (BUSCO) genes extracted from genomic and transcriptomic datasets. Phylogenomic inferences based on these data, using both coalescent-based and supermatrix analyses, are largely congruent with the most comprehensive plastome-based analysis, and nuclear-gene phylogenomic analyses with less comprehensive taxon sampling. The strongest discordance between the plastome and nuclear gene analyses is the monophyly of a clade comprising Asparagales and Liliales in our nuclear gene analyses, versus the placement of Asparagales and Liliales as successive sister clades to the commelinids in the plastome tree. Within orders, around six of 72 families shifted positions relative to the recent plastome analysis, but four of these involve poorly supported inferred relationships in the plastome-based tree. In Poales, the nuclear data place a clade comprising Ecdeiocoleaceae+Joinvilleaceae as sister to the grasses (Poaceae); Typhaceae, (rather than Bromeliaceae) are resolved as sister to all other Poales. In Commelinales, nuclear data place Philydraceae sister to all other families rather than to a clade comprising Haemodoraceae+Pontederiaceae as seen in the plastome tree. In Liliales, nuclear data place Liliaceae sister to Smilacaceae, and Melanthiaceae are placed sister to all other Liliales except Campynemataceae. Finally, in Alismatales, nuclear data strongly place Tofieldiaceae, rather than Araceae, as sister to all the other families, providing an alternative resolution of what has been the most problematic node to resolve using plastid data, outside of those involving achlorophyllous mycoheterotrophs. As seen in numerous prior studies, the placement of orders Acorales and Alismatales as successive sister lineages to all other extant monocots. Only 21.2% of BUSCO genes were demonstrably single-copy, yet phylogenomic inferences based on BUSCO and CSC genes did not differ, and overall functional annotations of the two sets were very similar. Our analyses also reveal significant gene tree-species tree discordance despite high support values, as expected given incomplete lineage sorting (ILS) related to rapid diversification. Our study advances understanding of monocot relationships and the robustness of phylogenetic inferences based on large numbers of nuclear single-copy genes that can be obtained from transcriptomes and genomes.

Potential allopolyploid origin of Ericales revealed with gene-tree reconciliation

Few incidents of ancient allopolyploidization (polyploidization by hybridization or merging diverged genomes) were previously revealed, although there is significant evidence for the accumulation of whole genome duplications (WGD) in plants. Here, we focused on Ericales, one of the largest and most diverse angiosperm orders with significant ornamental and economic value. Through integrating 24 high-quality whole genome data selected from ~ 200 Superasterids genomes/species and an algorithm of topology-based gene-tree reconciliation, we explored the evolutionary history of in Ericales with ancient complex. We unraveled the allopolyploid origin of Ericales and detected extensive lineage-specific gene loss following the polyploidization. Our study provided a new hypothesis regarding the origin of Ericales and revealed an instructive perspective of gene loss as a pervasive source of genetic variation and adaptive phenotypic diversity in Ericales.

Core circadian clock and light signaling genes brought into genetic linkage across the green lineage

The circadian clock is conserved at both the level of transcriptional networks as well as core genes in plants, ensuring that biological processes are phased to the correct time of day. In the model plant Arabidopsis (Arabidopsis thaliana), the core circadian SHAQKYF-type-MYB (sMYB) genes CIRCADIAN CLOCK ASSOCIATED 1 (CCA1) and REVEILLE (RVE4) show genetic linkage with PSEUDO-RESPONSE REGULATOR 9 (PRR9) and PRR7, respectively. Leveraging chromosome-resolved plant genomes and syntenic ortholog analysis enabled tracing this genetic linkage back to Amborella trichopoda, a sister lineage to the angiosperm, and identifying an additional evolutionarily conserved genetic linkage in light signaling genes. The LHY/CCA1-PRR5/9, RVE4/8-PRR3/7, and PIF3-PHYA genetic linkages emerged in the bryophyte lineage and progressively moved within several genes of each other across an array of angiosperm families representing distinct whole-genome duplication and fractionation events. Soybean (Glycine max) maintained all but two genetic linkages, and expression analysis revealed the PIF3-PHYA linkage overlapping with the E4 maturity group locus was the only pair to robustly cycle with an evening phase, in contrast to the sMYB-PRR morning and midday phase. While most monocots maintain the genetic linkages, they have been lost in the economically important grasses (Poaceae), such as maize (Zea mays), where the genes have been fractionated to separate chromosomes and presence/absence variation results in the segregation of PRR7 paralogs across heterotic groups. The environmental robustness model is put forward, suggesting that evolutionarily conserved genetic linkages ensure superior microhabitat pollinator synchrony, while wide-hybrids or unlinking the genes, as seen in the grasses, result in heterosis, adaptation, and colonization of new ecological niches.

A chromosome-level reference genome of Ensete glaucum gives insight into diversity and chromosomal and repetitive sequence evolution in the Musaceae

BACKGROUND

Ensete glaucum (2n = 2x = 18) is a giant herbaceous monocotyledonous plant in the small Musaceae family along with banana (Musa). A high-quality reference genome sequence assembly of E. glaucum is a resource for functional and evolutionary studies of Ensete, Musaceae, and the Zingiberales.

FINDINGS

Using Oxford Nanopore Technologies, chromosome conformation capture (Hi-C), Illumina and RNA survey sequence, supported by molecular cytogenetics, we report a high-quality 481.5 Mb genome assembly with 9 pseudo-chromosomes and 36,836 genes. A total of 55% of the genome is composed of repetitive sequences with predominantly LTR-retroelements (37%) and DNA transposons (7%). The single 5S ribosomal DNA locus had an exceptionally long monomer length of 1,056 bp, more than twice that of the monomers at multiple loci in Musa. A tandemly repeated satellite (1.1% of the genome, with no similar sequence in Musa) was present around all centromeres, together with a few copies of a long interspersed nuclear element (LINE) retroelement. The assembly enabled us to characterize in detail the chromosomal rearrangements occurring between E. glaucum and the x = 11 species of Musa. One E. glaucum chromosome has the same gene content as Musa acuminata, while others show multiple, complex, but clearly defined evolutionary rearrangements in the change between x= 9 and 11.

CONCLUSIONS

The advance towards a Musaceae pangenome including E. glaucum, tolerant of extreme environments, makes a complete set of gene alleles, copy number variation, and a reference for structural variation available for crop breeding and understanding environmental responses. The chromosome-scale genome assembly shows the nature of chromosomal fusion and translocation events during speciation, and features of rapid repetitive DNA change in terms of copy number, sequence, and genomic location, critical to understanding its role in diversity and evolution.

Unlocking plant metabolic diversity: A (pan)-genomic view

Plants produce a remarkable diversity of structurally and functionally diverse natural chemicals that serve as adaptive compounds throughout their life cycles. However, unlocking this metabolic diversity is significantly impeded by the size, complexity, and abundant repetitive elements of typical plant genomes. As genome sequencing becomes routine, we anticipate that links between metabolic diversity and genetic variation will be strengthened. In addition, an ever-increasing number of plant genomes have revealed that biosynthetic gene clusters are not only a hallmark of microbes and fungi; gene clusters for various classes of compounds have also been found in plants, and many are associated with important agronomic traits. We present recent examples of plant metabolic diversification that have been discovered through the exploration and exploitation of various genomic and pan-genomic data. We also draw attention to the fundamental genomic and pan-genomic basis of plant chemodiversity and discuss challenges and future perspectives for investigating metabolic diversity in the coming pan-genomics era.

Emergence of distinct syntenic density regimes is associated with early metazoan genomic transitions

Background

Animal genomes are strikingly conserved in terms of local gene order (microsynteny). While some of these microsyntenies have been shown to be coregulated or to form gene regulatory blocks, the diversity of their genomic and regulatory properties across the metazoan tree of life remains largely unknown.

Results

Our comparative analyses of 49 animal genomes reveal that the largest gains of synteny occurred in the last common ancestor of bilaterians and cnidarians and in that of bilaterians. Depending on their node of emergence, we further show that novel syntenic blocks are characterized by distinct functional compositions (Gene Ontology terms enrichment) and gene density properties, such as high, average and low gene density regimes. This is particularly pronounced among bilaterian novel microsyntenies, most of which fall into high gene density regime associated with higher gene coexpression levels. Conversely, a majority of vertebrate novel microsyntenies display a low gene density regime associated with lower gene coexpression levels.

Conclusions

Our study provides first evidence for evolutionary transitions between different modes of microsyntenic block regulation that coincide with key events of metazoan evolution. Moreover, the microsyntenic profiling strategy and interactive online application (Syntenic Density Browser, available at: http://synteny.csb.univie.ac.at/) we present here can be used to explore regulatory properties of microsyntenic blocks and predict their coexpression in a wide-range of animal genomes.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08304-2.

Why sequence all eukaryotes?

Life on Earth has evolved from initial simplicity to the astounding complexity we experience today. Bacteria and archaea have largely excelled in metabolic diversification, but eukaryotes additionally display abundant morphological innovation. How have these innovations come about and what constraints are there on the origins of novelty and the continuing maintenance of biodiversity on Earth? The history of life and the code for the working parts of cells and systems are written in the genome. The Earth BioGenome Project has proposed that the genomes of all extant, named eukaryotes-about 2 million species-should be sequenced to high quality to produce a digital library of life on Earth, beginning with strategic phylogenetic, ecological, and high-impact priorities. Here we discuss why we should sequence all eukaryotic species, not just a representative few scattered across the many branches of the tree of life. We suggest that many questions of evolutionary and ecological significance will only be addressable when whole-genome data representing divergences at all of the branchings in the tree of life or all species in natural ecosystems are available. We envisage that a genomic tree of life will foster understanding of the ongoing processes of speciation, adaptation, and organismal dependencies within entire ecosystems. These explorations will resolve long-standing problems in phylogenetics, evolution, ecology, conservation, agriculture, bioindustry, and medicine.