Nuclear phylotranscriptomics and phylogenomics support numerous polyploidization events and hypotheses for the evolution of rhizobial nitrogen-fixing symbiosis in Fabaceae

Nuclear phylogenomics of Asteraceae with increased sampling provides new insights into convergent morphological and molecular evolution

Convergent morphological evolution is widespread in flowering plants, and understanding this phenomenon relies on well-resolved phylogenies. Nuclear phylogenetic reconstruction using transcriptome datasets has been successful in various angiosperm groups, but it is limited to taxa with available fresh materials. Asteraceae, which are one of the two largest angiosperm families and are important for both ecosystems and human livelihood, show multiple examples of convergent evolution. Nuclear Asteraceae phylogenies have resolved relationships among most subfamilies and many tribes, but many phylogenetic and evolutionary questions regarding subtribes and genera remain, owing to limited sampling. Here, we increased the sampling for Asteraceae phylogenetic reconstruction using transcriptomes and genome-skimming datasets and produced nuclear phylogenetic trees with 706 species representing two-thirds of recognized subtribes. Ancestral character reconstruction supports multiple convergent evolutionary events in Asteraceae, with gains and losses of bilateral floral symmetry correlated with diversification of some subfamilies and smaller groups, respectively. Presence of the calyx-related pappus may have been especially important for the success of some subtribes and genera. Molecular evolutionary analyses support the likely contribution of duplications of MADS-box and TCP floral regulatory genes to innovations in floral morphology, including capitulum inflorescences and bilaterally symmetric flowers, potentially promoting the diversification of Asteraceae. Subsequent divergences and reductions in CYC2 gene expression are related to the gain and loss of zygomorphic flowers. This phylogenomic work with greater taxon sampling through inclusion of genome-skimming datasets reveals the feasibility of expanded evolutionary analyses using DNA samples for understanding convergent evolution.

High-integrity Pueraria montana var. lobata genome and population analysis revealed the genetic diversity of Pueraria genus

Pueraria montana var. lobata (P. lobata) is a traditional medicinal plant belonging to the Pueraria genus of Fabaceae family. Pueraria montana var. thomsonii (P. thomsonii) and Pueraria montana var. montana (P. montana) are its related species. However, evolutionary history of the Pueraria genus is still largely unknown. Here, a high-integrity, chromosome-level genome of P. lobata and an improved genome of P. thomsonii were reported. It found evidence for an ancient whole-genome triplication and a recent whole-genome duplication shared with Fabaceae in three Pueraria species. Population genomics of 121 Pueraria accessions demonstrated that P. lobata populations had substantially higher genetic diversity, and P. thomsonii was probably derived from P. lobata by domestication as a subspecies. Selection sweep analysis identified candidate genes in P. thomsonii populations associated with the synthesis of auxin and gibberellin, which potentially play a role in the expansion and starch accumulation of tubers in P. thomsonii. Overall, the findings provide new insights into the evolutionary and domestication history of the Pueraria genome and offer a valuable genomic resource for the genetic improvement of these species.

Chromosome-level genome assembly and characterization of the Calophaca sinica genome

Calophaca sinica is a rare plant endemic to northern China which belongs to the Fabaceae family and possesses rich nutritional value. To support the preservation of the genetic resources of this plant, we have successfully generated a high-quality genome of C. sinica (1.06 Gb). Notably, transposable elements (TEs) constituted ~73% of the genome, with long terminal repeat retrotransposons (LTR-RTs) dominating this group of elements (~54% of the genome). The average intron length of the C. sinica genome was noticeably longer than what has been observed for closely related species. The expansion of LTR-RTs and elongated introns emerged had the largest influence on the enlarged genome size of C. sinica in comparison to other Fabaceae species. The proliferation of TEs could be explained by certain modes of gene duplication, namely, whole genome duplication (WGD) and dispersed duplication (DSD). Gene family expansion, which was found to enhance genes associated with metabolism, genetic maintenance, and environmental stress resistance, was a result of transposed duplicated genes (TRD) and WGD. The presented genomic analysis sheds light on the genetic architecture of C. sinica, as well as provides a starting point for future evolutionary biology, ecology, and functional genomics studies centred around C. sinica and closely related species.

Phylogenomic profiles of whole-genome duplications in Poaceae and landscape of differential duplicate retention and losses among major Poaceae lineages

Poaceae members shared a whole-genome duplication called rho. However, little is known about the evolutionary pattern of the rho-derived duplicates among Poaceae lineages and implications in adaptive evolution. Here we present phylogenomic/phylotranscriptomic analyses of 363 grasses covering all 12 subfamilies and report nine previously unknown whole-genome duplications. Furthermore, duplications from a single whole-genome duplication were mapped to multiple nodes on the species phylogeny; a whole-genome duplication was likely shared by woody bamboos with possible gene flow from herbaceous bamboos; and recent paralogues of a tetraploid Oryza are implicated in tolerance of seawater submergence. Moreover, rho duplicates showing differential retention among subfamilies include those with functions in environmental adaptations or morphogenesis, including ACOT for aquatic environments (Oryzoideae), CK2β for cold responses (Pooideae), SPIRAL1 for rapid cell elongation (Bambusoideae), and PAI1 for drought/cold responses (Panicoideae). This study presents a Poaceae whole-genome duplication profile with evidence for multiple evolutionary mechanisms that contribute to gene retention and losses.

Haplotype-resolved genome of Mimosa bimucronata revealed insights into leaf movement and nitrogen fixation

BACKGROUND

Mimosa bimucronata originates from tropical America and exhibits distinctive leaf movement characterized by a relative slow speed. Additionally, this species possesses the ability to fix nitrogen. Despite these intriguing traits, comprehensive studies have been hindered by the lack of genomic resources for M. bimucronata.

RESULTS

To unravel the intricacies of leaf movement and nitrogen fixation, we successfully assembled a high-quality, haplotype-resolved, reference genome at the chromosome level, spanning 648 Mb and anchored in 13 pseudochromosomes. A total of 32,146 protein-coding genes were annotated. In particular, haplotype A was annotated with 31,035 protein-coding genes, and haplotype B with 31,440 protein-coding genes. Structural variations (SVs) and allele specific expression (ASE) analyses uncovered the potential role of structural variants in leaf movement and nitrogen fixation in M. bimucronata. Two whole-genome duplication (WGD) events were detected, that occurred ~ 2.9 and ~ 73.5 million years ago. Transcriptome and co-expression network analyses revealed the involvement of aquaporins (AQPs) and Ca2+-related ion channel genes in leaf movement. Moreover, we also identified nodulation-related genes and analyzed the structure and evolution of the key gene NIN in the process of symbiotic nitrogen fixation (SNF).

CONCLUSION

The detailed comparative genomic and transcriptomic analyses provided insights into the mechanisms governing leaf movement and nitrogen fixation in M. bimucronata. This research yielded genomic resources and provided an important reference for functional genomic studies of M. bimucronata and other legume species.

Out of Place: Phylogenomics resolves the placement of Eurasian taxa and sheds light on origin of Thermopsideae in North America

The North American Thermopsideae (Fabaceae: Papilionoideae), a monophyletic group comprising the North American endemic genus Baptisia, and the paraphyletic Eurasian-North American disjunct Thermopsis, is nested within the tribe Sophoreae. Previous phylogenetic studies have identified two East Asian taxa within the North American Thermopsideae, suggesting two independent dispersal events between North America-East Asia. More recent studies have also placed a third taxon, Vuralia turcica, an endemic species from Turkey, among the North American Thermopsideae. The presence of three geographically distant Eurasian taxa within a relatively young clade of North American origin is unprecedented among papilionoid legumes, and the biogeographic implications of this observation are not clear. To investigate this matter, 1540 low-copy nuclear genes and complete plastomes were obtained from 36 taxa across the core genistoids, including 26 newly sequenced taxa. Nuclear and plastome based maximum likelihood (ML) and ASTRAL analyses were conducted based on varying degrees of taxon coverage and read mapping consensus threshold values. Additional analyses were performed to estimate divergence times and to reconstruct biogeographic history. The results strongly support a previously undetected Old World clade, presently composed of V. turcica and T. chinensis, which diverged from the ancestor of the North American lineage during the mid to late Miocene. A single and recent North America-East Asia dispersal involving T. lupinoides is reported. Furthermore, the traditional inclusion of the genus Ammopiptanthus among Thermopsideae is not supported, and the monotypic generic status of Vuralia is called into question. A relatively high degree of cytonuclear discordance is reported within each sub-clade of the North American Thermopsideae. This finding is likely attributable to the high degree of interspecific hybridization reported within these groups and raises the need for more rigorous genome-scale testing to better delimit species within each of the reticulating subclades. Subjects: Biodiversity, Biogeography, Evolutionary Studies, Genetics, Plant Science.

Haplotype-resolved genome assembly of Phanera championii reveals molecular mechanisms of flavonoid synthesis and adaptive evolution

Phanera championii is a medicinal liana plant that has successfully adapted to hostile karst habitats. Despite extensive research on its medicinal components and pharmacological effects, the molecular mechanisms underlying the biosynthesis of critical flavonoids and its adaptation to karst habitats remain elusive. In this study, we performed high-coverage PacBio and Hi-C sequencing of P. championii, which revealed its high heterozygosity and phased the genome into two haplotypes: Hap1 (384.60 Mb) and Hap2 (383.70 Mb), encompassing a total of 58 612 annotated genes. Comparative genomes analysis revealed that P. championii experienced two whole-genome duplications (WGDs), with approximately 59.59% of genes originating from WGD events, thereby providing a valuable genetic resource for P. championii. Moreover, we identified a total of 112 genes that were strongly positively selected. Additionally, about 81.60 Mb of structural variations between the two haplotypes. The allele-specific expression patterns suggested that the dominant effect of P. championii was the elimination of deleterious mutations and the promotion of beneficial mutations to enhance fitness. Moreover, our transcriptome and metabolome analysis revealed alleles in different tissues or different haplotypes collectively regulate the synthesis of flavonoid metabolites. In summary, our comprehensive study highlights the significance of genomic and morphological adaptation in the successful adaptation of P. championii to karst habitats. The high-quality phased genomes obtained in this study serve as invaluable genomic resources for various applications, including germplasm conservation, breeding, evolutionary studies, and elucidation of pathways governing key biological traits of P. championii.

Chromosome-level genome of the transformable northern wattle, Acacia crassicarpa

The genus Acacia is a large group of woody legumes containing an enormous amount of morphological diversity in leaf shape. This diversity is at least in part the result of an innovation in leaf development where many Acacia species are capable of developing leaves of both bifacial and unifacial morphologies. While not unique in the plant kingdom, unifaciality is most commonly associated with monocots, and its developmental genetic mechanisms have yet to be explored beyond this group. In this study, we identify an accession of Acacia crassicarpa with high regeneration rates and isolate a clone for genome sequencing. We generate a chromosome-level assembly of this readily transformable clone, and using comparative analyses, confirm a whole-genome duplication unique to Caesalpinoid legumes. This resource will be important for future work examining genome evolution in legumes and the unique developmental genetic mechanisms underlying unifacial morphogenesis in Acacia.

Advances in pharmacology, biosynthesis, and metabolic engineering of Scutellaria-specialized metabolites

Scutellaria Linn., which belongs to the family Lamiaceae, is a commonly used medicinal plant for heat clearing and detoxification. In particular, the roots of S. baicalensis and the entire herb of S. barbata have been widely used in traditional medicine for thousands of years. The main active components of Scutellaria, including: baicalein, wogonin, norwogonin, scutellarein, and their glycosides have potential or existing drug usage. However, the wild resources of Scutellaria plants have been overexploited, and degenerated germplasm resources cannot fulfill the requirements of chemical extraction and clinical usage. Metabolic engineering and green production via microorganisms provide alternative strategies for greater efficiency in the production of natural products. Here, we review the progress of: pharmacological investigations, multi-omics, biosynthetic pathways, and metabolic engineering of various Scutellaria species and their active compounds. In addition, based on multi-omics data, we systematically analyze the phylogenetic relationships of Scutellaria and predict candidate transcription factors related to the regulation of active flavonoids. Finally, we propose the prospects of directed evolution of core enzymes and genome-assisted breeding to alleviate the shortage of plant resources of Scutellaria. This review provides important insights into the sustainable utilization and development of Scutellaria resources.

Nuclear phylogenomics of angiosperms and insights into their relationships and evolution

Angiosperms (flowering plants) are by far the most diverse land plant group with over 300,000 species. The sudden appearance of diverse angiosperms in the fossil record was referred to by Darwin as the "abominable mystery," hence contributing to the heightened interest in angiosperm evolution. Angiosperms display wide ranges of morphological, physiological, and ecological characters, some of which have probably influenced their species richness. The evolutionary analyses of these characteristics help to address questions of angiosperm diversification and require well resolved phylogeny. Following the great successes of phylogenetic analyses using plastid sequences, dozens to thousands of nuclear genes from next-generation sequencing have been used in angiosperm phylogenomic analyses, providing well resolved phylogenies and new insights into the evolution of angiosperms. In this review we focus on recent nuclear phylogenomic analyses of large angiosperm clades, orders, families, and subdivisions of some families and provide a summarized Nuclear Phylogenetic Tree of Angiosperm Families. The newly established nuclear phylogenetic relationships are highlighted and compared with previous phylogenetic results. The sequenced genomes of Amborella, Nymphaea, Chloranthus, Ceratophyllum, and species of monocots, Magnoliids, and basal eudicots, have facilitated the phylogenomics of relationships among five major angiosperms clades. All but one of the 64 angiosperm orders were included in nuclear phylogenomics with well resolved relationships except the placements of several orders. Most families have been included with robust and highly supported placements, especially for relationships within several large and important orders and families. Additionally, we examine the divergence time estimation and biogeographic analyses of angiosperm on the basis of the nuclear phylogenomic frameworks and discuss the differences compared with previous analyses. Furthermore, we discuss the implications of nuclear phylogenomic analyses on ancestral reconstruction of morphological, physiological, and ecological characters of angiosperm groups, limitations of current nuclear phylogenomic studies, and the taxa that require future attention.

Whole Genome Duplication Events Likely Contributed to the Aquatic Adaptive Evolution of Parkerioideae

As the only aquatic lineage of Pteridaceae, Parkerioideae is distinct from many xeric-adapted species of the family and consists of the freshwater Ceratopteris species and the only mangrove ferns from the genus Acrostichum. Previous studies have shown that whole genome duplication (WGD) has occurred in Parkerioideae at least once and may have played a role in their adaptive evolution; however, more in-depth research regarding this is still required. In this study, comparative and evolutionary transcriptomics analyses were carried out to identify WGDs and explore their roles in the environmental adaptation of Parkerioideae. Three putative WGD events were identified within Parkerioideae, two of which were specific to Ceratopteris and Acrostichum, respectively. The functional enrichment analysis indicated that the lineage-specific WGD events have played a role in the adaptation of Parkerioideae to the low oxygen concentrations of aquatic habitats, as well as different aquatic environments of Ceratopteris and Acrostichum, such as the adaptation of Ceratopteris to reduced light levels and the adaptation of Acrostichum to high salinity. Positive selection analysis further provided evidence that the putative WGD events may have facilitated the adaptation of Parkerioideae to changes in habitat. Moreover, the gene family analysis indicated that the plasma membrane H+-ATPase (AHA), vacuolar H+-ATPase (VHA), and suppressor of K+ transport growth defect 1 (SKD1) may have been involved in the high salinity adaptation of Acrostichum. Our study provides new insights into the evolution and adaptations of Parkerioideae in different aquatic environments.

Genome-Wide Identification of the Paulownia fortunei Aux/IAA Gene Family and Its Response to Witches' Broom Caused by Phytoplasma

The typical symptom of Paulownia witches' broom (PaWB), caused by phytoplasma infection, is excessive branching, which is mainly triggered by auxin metabolism disorder. Aux/IAA is the early auxin-responsive gene that participates in regulating plant morphogenesis such as apical dominance, stem elongation, lateral branch development, and lateral root formation. However, no studies have investigated the response of the Aux/IAA gene family to phytoplasma infection in Paulownia fortunei. In this study, a total of 62 Aux/IAA genes were found in the genome. Phylogenetic analysis showed that PfAux/IAA genes could be divided into eight subgroups, which were formed by tandem duplication and fragment replication. Most of them had a simple gene structure, and several members lacked one or two conserved domains. By combining the expression of PfAux/IAA genes under phytoplasma stress and SA-treated phytoplasma-infected seedlings, we found that PfAux/IAA13/33/45 may play a vital role in the occurrence of PaWB. Functional analysis based on homologous relationships showed a strong correlation between PfAux/IAA45 and branching. Protein-protein interaction prediction showed that PfARF might be the binding partner of PfAux/IAA, and the yeast two-hybrid assay and bimolecular fluorescent complementary assay confirmed the interaction of PfAux/IAA45 and PfARF13. This study provides a theoretical basis for further understanding the function of the PfAux/IAA gene family and exploring the regulatory mechanism of branching symptoms caused by PaWB.

Characteristics of chloroplast and mitochondrial genomes and intracellular gene transfer in the Korean endemic shrub, Sophora koreensis Nakai (Fabaceae)

Sophora koreensis Nakai, an endemic species distributed only in the Korean Peninsula, is of great geographical, economic, and taxonomic importance. Although its complete chloroplast (cp) genome sequence has been reported, its mitochondrial (mt) genome sequence has not yet been studied. Therefore, in this study, we aimed to investigate its mt genome sequence and compare it with those reported for other Fabaceae species. Total genomic DNA was extracted from fresh S. koreensis leaves collected from natural habitats in Gangwon-do Province, South Korea. This was followed by polymerase chain reaction (PCR) amplification of cpDNA insertions in the mt genome and the detection of microsatellites and dispersed repeats in the cp and mt genomes. Finally, the cp and mt genomes of S. koreensis were compared with those reported for other Fabaceae species. The cp sequence of S. koreensis showed identical gene orders and contents as those previously reported. Only six substitutions and one deletion were detected with 99 % homology. Conversely, the complete mt genome sequence, which was 517,845 bp in length and encoded 61 genes, including 43 protein-coding, 15 transfer RNAs, and 3 ribosomal RNA genes, was considerably different from that of S. japonica in terms of gene order and composition. Further, the mt genome of S. koreensis included ca. 7 and 3 kb insertions, representing an intracellular gene transfer (IGT) event, and the regions with these insertions were determined to be originally present in the cp genome. This IGT event was also confirmed via PCR amplification. IGT events can be induced via biological gene expression control or the use of repetitive sequences, and they provide important insights into the evolutionary lineage of S. koreensis. However, further studies are needed to clarify the gene transfer mechanisms between the two organelles.

The Fabaceae in Northeastern Mexico (Subfamily Caesalpinioideae, Mimosoideae Clade, Tribes Mimoseae, Acacieae, and Ingeae)

A synoptic compendium of the legumes of the Mimosoideae clade in northeastern Mexico is presented for the first time, including changes in their botanical nomenclature and retypification of genera. Furthermore, based on new information recently published, the taxonomic limits of several new genera segregated from Acacia (Acaciella, Mariosousa, Senegalia, and Vachellia) and Prosopis (Neltuma and Strombocarpa) are clarified and included. Based on field work, collection of botanical samples over the past 30 years, and reviewing botanical materials in national and international herbaria, we have completed the diversity of legumes of the Mimosoideae clade of northeastern Mexico. Three tribes (Acacieae, Ingeae, and Mimosaeae), 22 genera, 92 species, and 19 infraspecific categories were recorded. Only the genus Painteria is endemic to Mexico. Eighty-eight species are native to Mexico, and four are exotic: Acacia salicina, Neptunia prostrata, Neltuma chilensis and Albizia lebbeck. Twenty-eight species are endemic to Mexico, nine species are endemic to northeastern Mexico, and four species are endemic to only one state in Mexico. The 22 registered genera represent 44% and 65% of the generic flora of the Mimosoideae clade for Mexico and the planet, respectively, while the 92 species registered represent 3% and 18% of the species of the clade Mimosoideae for the planet and Mexico, respectively. According to the new nomenclature of legumes, the number of genera in the Mimosoideae clade in northern Mexico has increased from 19 to 24.

The complete mitochondrial genome of Bauhinia variegata (Leguminosae)

The mitogenome of Bauhinia variegate was assembled and characterized in this study. The mitogenome size was 437,271 bp, and its GC content was 45.5%. 36 protein-coding genes, 17 tRNAs and 3 rRNAs were annotated in the mitogenome. A total of 12 MTPTs, ranging from 71 bp to 3562 bp, were identified in the mitogenome and covered 1.46% (6373 bp) of the mitogenome. Phylogenetic analysis of 15 species of Leguminosae based on 23 core protein-coding genes showed that B. variegata was sister to Tylosema esculentum, another member from the subfamily Cercidoideae. The mitogenome of B. variegata provides a valuable genetic resource for further phylogenetic studies of this family.

Identification and Expression Analysis of the Nucleotidyl Transferase Protein (NTP) Family in Soybean (Glycine max) under Various Abiotic Stresses

Nucleotidyl transferases (NTPs) are common transferases in eukaryotes and play a crucial role in nucleotide modifications at the 3' end of RNA. In plants, NTPs can regulate RNA stability by influencing 3' end modifications, which in turn affect plant growth, development, stress responses, and disease resistance. Although the functions of NTP family members have been extensively studied in Arabidopsis, rice, and maize, there is limited knowledge about NTP genes in soybeans. In this study, we identified 16 members of the NTP family in soybeans, including two subfamilies (G1 and G2) with distinct secondary structures, conserved motifs, and domain distributions at the protein level. Evolutionary analysis of genes in the NTP family across multiple species and gene collinearity analysis revealed a relatively conserved evolutionary pattern. Analysis of the tertiary structure of the proteins showed that NTPs have three conserved aspartic acids that bind together to form a possible active site. Tissue-specific expression analysis indicated that some NTP genes exhibit tissue-specific expression, likely due to their specific functions. Stress expression analysis showed significant differences in the expression levels of NTP genes under high salt, drought, and cold stress. Additionally, RNA-seq analysis of soybean plants subjected to salt and drought stress further confirmed the association of soybean NTP genes with abiotic stress responses. Subcellular localization experiments revealed that GmNTP2 and GmNTP14, which likely have similar functions to HESO1 and URT1, are located in the nucleus. These research findings provide a foundation for further investigations into the functions of NTP family genes in soybeans.

Comparative phylogenomics and phylotranscriptomics provide insights into the genetic complexity of nitrogen-fixing root-nodule symbiosis

Plant root-nodule symbiosis (RNS) with mutualistic nitrogen-fixing bacteria is restricted to a single clade of angiosperms, the Nitrogen-Fixing Nodulation Clade (NFNC), and is best understood in the legume family. Nodulating species share many commonalities, explained either by divergence from a common ancestor over 100 million years ago or by convergence following independent origins over that same time period. Regardless, comparative analyses of diverse nodulation syndromes can provide insights into constraints on nodulation-what must be acquired or cannot be lost for a functional symbiosis-and the latitude for variation in the symbiosis. However, much remains to be learned about nodulation, especially outside of legumes. Here, we employed a large-scale phylogenomic analysis across 88 species, complemented by 151 RNA-seq libraries, to elucidate the evolution of RNS. Our phylogenomic analyses further emphasize the uniqueness of the transcription factor NIN as a master regulator of nodulation and identify key mutations that affect its function across the NFNC. Comparative transcriptomic assessment revealed nodule-specific upregulated genes across diverse nodulating plants, while also identifying nodule-specific and nitrogen-response genes. Approximately 70% of symbiosis-related genes are highly conserved in the four representative species, whereas defense-related and host-range restriction genes tend to be lineage specific. Our study also identified over 900 000 conserved non-coding elements (CNEs), over 300 000 of which are unique to sampled NFNC species. NFNC-specific CNEs are enriched with the active H3K9ac mark and are correlated with accessible chromatin regions, thus representing a pool of candidate regulatory elements for genes involved in RNS. Collectively, our results provide novel insights into the evolution of nodulation and lay a foundation for engineering of RNS traits in agriculturally important crops.

Phylotranscriptomic analyses reveal deep gene tree discordance in Camellia (Theaceae)

Gene tree discordance is a significant legacy of biological evolution. Multiple factors can result in incongruence among genes, such as introgression, incomplete lineage sorting (ILS), gene duplication or loss. Resolving the background of gene tree discordance is a critical way to uncover the process of species diversification. Camellia, the largest genus in Theaceae, has controversial taxonomy and systematics due in part to a complex evolutionary history. We used 60 transcriptomes of 55 species, which represented 15 sections of Camellia to investigate its phylogeny and the possible causes of gene tree discordance. We conducted gene tree discordance analysis based on 1,617 orthologous low-copy nuclear genes, primarily using coalescent species trees and polytomy tests to distinguish hard and soft conflict. A selective pressure analysis was also performed to assess the impact of selection on phylogenetic topology reconstruction. Our results detected different levels of gene tree discordance in the backbone of Camellia, and recovered rapid diversification as one of the possible causes of gene tree discordance. Furthermore, we confirmed that none of the currently proposed sections of Camellia was monophyletic. Comparisons among datasets partitioned under different selective pressure regimes showed that integrating all orthologous genes provided the best phylogenetic resolution of the species tree of Camellia. The findings of this study reveal rapid diversification as a major source of gene tree discordance in Camellia and will facilitate future investigation of reticulate relationships at the species level in this important plant genus.

Phylogenomic Analyses Reveal an Allopolyploid Origin of Core Didymocarpinae (Gesneriaceae) Followed by Rapid Radiation

Allopolyploid plants have long been regarded as possessing genetic advantages under certain circumstances due to the combined effects of their hybrid origins and duplicated genomes. However, the evolutionary consequences of allopolyploidy in lineage diversification remain to be fully understood. Here, we investigate the evolutionary consequences of allopolyploidy using 138 transcriptomic sequences of Gesneriaceae, including 124 newly sequenced, focusing particularly on the largest subtribe Didymocarpinae. We estimated the phylogeny of Gesneriaceae using concatenated and coalescent-based methods based on five different nuclear matrices and 27 plastid genes, focusing on relationships among major clades. To better understand the evolutionary affinities in this family, we applied a range of approaches to characterize the extent and cause of phylogenetic incongruence. We found that extensive conflicts between nuclear and chloroplast genomes and among nuclear genes were caused by both incomplete lineage sorting (ILS) and reticulation, and we found evidence of widespread ancient hybridization and introgression. Using the most highly supported phylogenomic framework, we revealed multiple bursts of gene duplication throughout the evolutionary history of Gesneriaceae. By incorporating molecular dating and analyses of diversification dynamics, our study shows that an ancient allopolyploidization event occurred around the Oligocene-Miocene boundary, which may have driven the rapid radiation of core Didymocarpinae.

Parallel selection of loss-of-function alleles of Pdh1 orthologous genes in warm-season legumes for pod indehiscence and plasticity is related to precipitation

Pod dehiscence facilitates seed dispersal in wild legumes but results in yield loss in cultivated legumes. The evolutionary genetics of the legume pod dehiscence trait remain largely elusive. We characterized the pod dehiscence of chromosome segment substitution lines of Glycine max crossed with Glycine soja and found that the gene underlying the predominant quantitative trait locus (QTL) of soybean pod-shattering trait was Pod dehiscence 1 (Pdh1). A few rare loss-of-function (LoF) Pdh1 alleles were identified in G. soja, while only an allele featuring a premature stop codon was selected for pod indehiscence in cultivated soybean and spread to low-precipitation regions with increased frequency. Moreover, correlated interactions among Pdh1's haplotype, gene expression, and environmental changes for the developmental plasticity of the pod dehiscence trait were revealed in G. max. We found that orthologous Pdh1 genes specifically originated in warm-season legumes and their LoF alleles were then parallel-selected during the domestication of legume crops. Our results provide insights into the convergent evolution of pod dehiscence in warm-season legumes, facilitate an understanding of the intricate interactions between genetic robustness and environmental adaptation for developmental plasticity, and guide the breeding of new legume varieties with pod indehiscence.

The complete plastid genome sequence of Lysidice brevicalyx (Fabaceae: Detarioideae), an arbor species endemic to China

The plastid genome of Lysidice brevicalyx was successfully assembled using Illumina sequencing reads for the first time. The complete plastid genome of L. brevicalyx is a circular structure of 159,084 bp with a GC content of 36.4%. It comprises a large single-copy (LSC) region of 87,783 bp, a small single-copy (SSC) region of 19,557 bp, and two inverted repeat regions (IRA and IRB) of 25,872 bp, each. The plastome of L. brevicalyx contains a total of 128 genes, including 83 protein-coding genes, 37 tRNAs, and 8 rRNAs. The phylogenetic analysis strongly supports the monophyly of Lysidice. This study provides the first complete plastid genome sequence of L. brevicalyx and contributes to our understanding of the molecular characteristics and evolutionary relationships of this plant species.

Interploidy Introgression Shaped Adaptation during the Origin and Domestication History of Brassica napus

Polyploidy is recurrent across the tree of life and known as an evolutionary driving force in plant diversification and crop domestication. How polyploid plants adapt to various habitats has been a fundamental question that remained largely unanswered. Brassica napus is a major crop cultivated worldwide, resulting from allopolyploidy between unknown accessions of diploid B. rapa and B. oleracea. Here, we used whole-genome resequencing data of accessions representing the majority of morphotypes and ecotypes from the species B. rapa, B. oleracea, and B. napus to investigate the role of polyploidy during domestication. To do so, we first reconstructed the phylogenetic history of B. napus, which supported the hypothesis that the emergence of B. napus derived from the hybridization of European turnip of B. rapa and wild B. oleracea. These analyses also showed that morphotypes of swede and Siberian kale (used as vegetable and fodder) were domesticated before rapeseed (oil crop). We next observed that frequent interploidy introgressions from sympatric diploids were prominent throughout the domestication history of B. napus. Introgressed genomic regions were shown to increase the overall genetic diversity and tend to be localized in regions of high recombination. We detected numerous candidate adaptive introgressed regions and found evidence that some of the genes in these regions contributed to phenotypic diversification and adaptation of different morphotypes. Overall, our results shed light on the origin and domestication of B. napus and demonstrate interploidy introgression as an important mechanism that fuels rapid diversification in polyploid species.

Chromosome-level genome assembly and population genomics of Robinia pseudoacacia reveal the genetic basis for its wide cultivation

Urban greening provides important ecosystem services and ideal places for urban recreation and is a serious consideration for municipal decision-makers. Among the tree species cultivated in urban green spaces, Robinia pseudoacacia stands out due to its attractive flowers, fragrances, high trunks, wide adaptability, and essential ecosystem services. However, the genomic basis and consequences of its wide-planting in urban green spaces remains unknown. Here, we report the chromosome-level genome assembly of R. pseudoacacia, revealing a genome size of 682.4 Mb and 33,187 protein-coding genes. More than 99.3% of the assembly is anchored to 11 chromosomes with an N50 of 59.9 Mb. Comparative genomic analyses among 17 species reveal that gene families related to traits favoured by urbanites, such as wood formation, biosynthesis, and drought tolerance, are notably expanded in R. pseudoacacia. Our population genomic analyses further recover 11 genes that are under recent selection. Ultimately, these genes play important roles in the biological processes related to flower development, water retention, and immunization. Altogether, our results reveal the evolutionary forces that shape R. pseudoacacia cultivated for urban greening. These findings also present a valuable foundation for the future development of agronomic traits and molecular breeding strategies for R. pseudoacacia.

Lineage-specific amplification and epigenetic regulation of LTR-retrotransposons contribute to the structure, evolution, and function of Fabaceae species

BACKGROUND

Long terminal repeat (LTR)-retrotransposons (LTR-RTs) are ubiquitous and make up the majority of nearly all sequenced plant genomes, whereas their pivotal roles in genome evolution, gene expression regulation as well as their epigenetic regulation are still not well understood, especially in a large number of closely related species.

RESULTS

Here, we analyzed the abundance and dynamic evolution of LTR-RTs in 54 species from an economically and agronomically important family, Fabaceae, and also selected two representative species for further analysis in expression of associated genes, transcriptional activity and DNA methylation patterns of LTR-RTs. Annotation results revealed highly varied proportions of LTR-RTs in these genomes (5.1%~68.4%) and their correlation with genome size was highly positive, and they were significantly contributed to the variance in genome size through species-specific unique amplifications. Almost all of the intact LTR-RTs were inserted into the genomes 4 Mya (million years ago), and more than 50% of them were inserted in the last 0.5 million years, suggesting that recent amplifications of LTR-RTs were an important force driving genome evolution. In addition, expression levels of genes with intronic, promoter, and downstream LTR-RT insertions of Glycine max and Vigna radiata, two agronomically important crops in Fabaceae, showed that the LTR-RTs located in promoter or downstream regions suppressed associated gene expression. However, the LTR-RTs within introns promoted gene expression or had no contribution to gene expression. Additionally, shorter and younger LTR-RTs maintained higher mobility and transpositional potential. Compared with the transcriptionally silent LTR-RTs, the active elements showed significantly lower DNA methylation levels in all three contexts. The distributions of transcriptionally active and silent LTR-RT methylation varied across different lineages due to the position of LTR-RTs located or potentially epigenetic regulation.

CONCLUSION

Lineage-specific amplification patterns were observed and higher methylation level may repress the activity of LTR-RTs, further influence evolution in Fabaceae species. This study offers valuable clues into the evolution, function, transcriptional activity and epigenetic regulation of LTR-RTs in Fabaceae genomes.

The draft genome reveals early duplication event in Pterocarpus santalinus: an endemic timber species

MAIN CONCLUSION

A 541 Mb draft genome of Pterocarpus santalinus is presented and evidence of whole-genome duplication in the Eocene period with expansion of drought responsive gene families is documented. Pterocarpus santalinus Linn. f., popularly known as Red Sanders, is a deciduous tree, endemic to southern parts of Eastern Ghats in India. The heartwood is highly valued in the international market due to its deep red colour, fragrant heartwood and wavy grained texture. In the present study, a high-quality draft genome of P. santalinus was assembled using short and long reads generated from Illumina and Oxford Nanopore Sequencing platforms, respectively. The haploid genome size was estimated at 541 Mb and the hybrid assembly showed 99.60% genome completeness. A total of 51,713 consensus gene set were predicted with 31,437 annotated genes. The age of the whole-genome duplication event in the species was dated at 30-39 mya with 95% confidence suggesting early genome duplication event during the Eocene period. Concurrently, phylogenomic assessment of seven Papilionoideae members including P. santalinus grouped the species based on the tribal classification and established divergence of the tribe Dalbergieae from tribe Trifolieae at ~ 54.20 mya. A significant expansion of water deprivation/drought responsive gene families documented in the study probably explains the occurrence of the species in dry rocky patches. Additionally, re-sequencing of six diverse genotypes predicted one variant every 27 bases. This report presents the first draft genome in the genus Pterocarpus and the unprecedented genomic information generated is expected to accelerate population divergence studies in the species in relation to its endemic nature, support trait-based breeding programme and aid in development of diagnostic tools for timber forensics.

Diversity and regulation of symbiotic nitrogen fixation in plants

Plants associate with nitrogen-fixing bacteria to secure nitrogen, which is generally the most limiting nutrient for plant growth. Endosymbiotic nitrogen-fixing associations are widespread among diverse plant lineages, ranging from microalgae to angiosperms, and are primarily one of three types: cyanobacterial, actinorhizal or rhizobial. The large overlap in the signaling pathways and infection components of arbuscular mycorrhizal, actinorhizal and rhizobial symbioses reflects their evolutionary relatedness. These beneficial associations are influenced by environmental factors and other microorganisms in the rhizosphere. In this review, we summarize the diversity of nitrogen-fixing symbioses, key signal transduction pathways and colonization mechanisms relevant to such interactions, and compare and contrast these interactions with arbuscular mycorrhizal associations from an evolutionary standpoint. Additionally, we highlight recent studies on environmental factors regulating nitrogen-fixing symbioses to provide insights into the adaptation of symbiotic plants to complex environments.

Comprehensive phylogenetic analyses of Orchidaceae using nuclear genes and evolutionary insights into epiphytism

Orchidaceae (with >28,000 orchid species) are one of the two largest plant families, with economically and ecologically important species, and occupy global and diverse niches with primary distribution in rainforests. Among orchids, 70% grow on other plants as epiphytes; epiphytes contribute up to ~50% of the plant diversity in rainforests and provide food and shelter for diverse animals and microbes, thereby contributing to the health of these ecosystems. Orchids account for over two-thirds of vascular epiphytes and provide an excellent model for studying evolution of epiphytism. Extensive phylogenetic studies of Orchidaceae and subgroups have ;been crucial for understanding relationships among many orchid lineages, although some uncertainties remain. For example, in the largest subfamily Epidendroideae with nearly all epiphytic orchids, relationships among some tribes and many subtribes are still controversial, hampering evolutionary analyses of epiphytism. Here we obtained 1,450 low-copy nuclear genes from 610 orchid species, including 431 with newly generated transcriptomes, and used them for the reconstruction of robust Orchidaceae phylogenetic trees with highly supported placements of tribes and subtribes. We also provide generally well-supported phylogenetic placements of 131 genera and 437 species that were not sampled by previous plastid and nuclear phylogenomic studies. Molecular clock analyses estimated the Orchidaceae origin at ~132 million years ago (Ma) and divergences of most subtribes from 52 to 29 Ma. Character reconstruction supports at least 14 parallel origins of epiphytism; one such origin was placed at the most recent common ancestor of ~95% of epiphytic orchids and linked to modern rainforests. Ten occurrences of rapid increase in the diversification rate were detected within Epidendroideae near and after the K-Pg boundary, contributing to ~80% of the Orchidaceae diversity. This study provides a robust and the largest family-wide Orchidaceae nuclear phylogenetic tree thus far and new insights into the evolution of epiphytism in vascular plants.

Nuclear phylogeny and insights into whole-genome duplications and reproductive development of Solanaceae plants

Solanaceae, the nightshade family, have ∼2700 species, including the important crops potato and tomato, ornamentals, and medicinal plants. Several sequenced Solanaceae genomes show evidence for whole-genome duplication (WGD), providing an excellent opportunity to investigate WGD and its impacts. Here, we generated 93 transcriptomes/genomes and combined them with 87 public datasets, for a total of 180 Solanaceae species representing all four subfamilies and 14 of 15 tribes. Nearly 1700 nuclear genes from these transcriptomic/genomic datasets were used to reconstruct a highly resolved Solanaceae phylogenetic tree with six major clades. The Solanaceae tree supports four previously recognized subfamilies (Goetzeioideae, Cestroideae, Nicotianoideae, and Solanoideae) and the designation of three other subfamilies (Schizanthoideae, Schwenckioideae, and Petunioideae), with the placement of several previously unassigned genera. We placed a Solanaceae-specific whole-genome triplication (WGT1) at ∼81 million years ago (mya), before the divergence of Schizanthoideae from other Solanaceae subfamilies at ∼73 mya. In addition, we detected two gene duplication bursts (GDBs) supporting proposed WGD events and four other GDBs. An investigation of the evolutionary histories of homologs of carpel and fruit developmental genes in 14 gene (sub)families revealed that 21 gene clades have retained gene duplicates. These were likely generated by the Solanaceae WGT1 and may have promoted fleshy fruit development. This study presents a well-resolved Solanaceae phylogeny and a new perspective on retained gene duplicates and carpel/fruit development, providing an improved understanding of Solanaceae evolution.

A plastid phylogenomic framework for the palm family (Arecaceae)

BACKGROUND

Over the past decade, phylogenomics has greatly advanced our knowledge of angiosperm evolution. However, phylogenomic studies of large angiosperm families with complete species or genus-level sampling are still lacking. The palms, Arecaceae, are a large family with ca. 181 genera and 2600 species and are important components of tropical rainforests bearing great cultural and economic significance. Taxonomy and phylogeny of the family have been extensively investigated by a series of molecular phylogenetic studies in the last two decades. Nevertheless, some phylogenetic relationships within the family are not yet well-resolved, especially at the tribal and generic levels, with consequent impacts for downstream research.

RESULTS

Plastomes of 182 palm species representing 111 genera were newly sequenced. Combining these with previously published plastid DNA data, we were able to sample 98% of palm genera and conduct a plastid phylogenomic investigation of the family. Maximum likelihood analyses yielded a robustly supported phylogenetic hypothesis. Phylogenetic relationships among all five palm subfamilies and 28 tribes were well-resolved, and most inter-generic phylogenetic relationships were also resolved with strong support.

CONCLUSIONS

The inclusion of nearly complete generic-level sampling coupled with nearly complete plastid genomes strengthened our understanding of plastid-based relationships of the palms. This comprehensive plastid genome dataset complements a growing body of nuclear genomic data. Together, these datasets form a novel phylogenomic baseline for the palms and an increasingly robust framework for future comparative biological studies of this exceptionally important plant family.

Legumes from the Paleocene sediments of India and their ecological significance

During the early Paleogene, greenhouse gases created warm global climates. These warm climates redistributed the habitat of marine and terrestrial biota globally. Understanding the ecology of biotas under extremely warm climates is important to decipher their behavior in future climate warming. Here we report two new legume fossils (Leguminocarpum meghalayensis Bhatia, Srivastava et Mehrotra sp. nov., and Parvileguminophyllum damalgiriensis Bhatia, Srivastava et Mehrotra sp. nov.) from the late Paleocene sediments of Tura Formation of Meghalaya, northeast India. Globally, the Paleocene legume fossil records indicate that legumes most likely immigrated to India from Africa via the Ladakh-Kohistan Arc during the early Paleogene. Moreover, previously reconstructed climate data from the Tura Formation indicate that legumes were well adapted to a warm seasonal climate with monsoon rains.

Transcriptome-based variations effectively untangling the intraspecific relationships and selection signals in Xinyang Maojian tea population

As one of the world's top three popular non-alcoholic beverages, tea is economically and culturally valuable. Xinyang Maojian, this elegant green tea, is one of the top ten famous tea in China and has gained prominence for thousands of years. However, the cultivation history of Xinyang Maojian tea population and selection signals of differentiation from the other major variety Camellia sinensis var. assamica (CSA) remain unclear. We newly generated 94 Camellia sinensis (C. sinensis) transcriptomes including 59 samples in the Xinyang area and 35 samples collected from 13 other major tea planting provinces in China. Comparing the very low resolution of phylogeny inferred from 1785 low-copy nuclear genes with 94 C. sinensis samples, we successfully resolved the phylogeny of C. sinensis samples by 99,115 high-quality SNPs from the coding region. The sources of tea planted in the Xinyang area were extensive and complex. Specifically, Shihe District and Gushi County were the two earliest tea planting areas in Xinyang, reflecting a long history of tea planting. Furthermore, we identified numerous selection sweeps during the differentiation of CSA and CSS and these positive selection genes are involved in many aspects such as regulation of secondary metabolites synthesis, amino acid metabolism, photosynthesis, etc. Numerous specific selective sweeps of modern cultivars were annotated with functions in various different aspects, indicating the CSS and CSA populations possibly underwent independent specific domestication processes. Our study indicated that transcriptome-based SNP-calling is an efficient and cost-effective method in untangling intraspecific phylogenetic relationships. This study provides a significant understanding of the cultivation history of the famous Chinese tea Xinyang Maojian and unravels the genetic basis of physiological and ecological differences between the two major tea subspecies.

Phylotranscriptomic analyses reveal multiple whole-genome duplication events, the history of diversification and adaptations in the Araceae

BACKGROUND AND AIMS

The Araceae are one of the most diverse monocot families with numerous morphological and ecological novelties. Plastid and mitochondrial genes have been used to investigate the phylogeny and to interpret shifts in the pollination biology and biogeography of the Araceae. In contrast, the role of whole-genome duplication (WGD) in the evolution of eight subfamilies remains unclear.

METHODS

New transcriptomes or low-depth whole-genome sequences of 65 species were generated through Illumina sequencing. We reconstructed the phylogenetic relationships of Araceae using concatenated and species tree methods, and then estimated the age of major clades using TreePL. We inferred the WGD events by Ks and gene tree methods. We investigated the diversification patterns applying time-dependent and trait-dependent models. The expansions of gene families and functional enrichments were analysed using CAFE and InterProScan.

KEY RESULTS

Gymnostachydoideae was the earliest diverging lineage followed successively by Orontioideae, Lemnoideae and Lasioideae. In turn, they were followed by the clade of 'bisexual climbers' comprised of Pothoideae and Monsteroideae, which was resolved as the sister to the unisexual flowers clade of Zamioculcadoideae and Aroideae. A special WGD event ψ (psi) shared by the True-Araceae clade occurred in the Early Cretaceous. Net diversification rates first declined and then increased through time in the Araceae. The best diversification rate shift along the stem lineage of the True-Araceae clade was detected, and net diversification rates were enhanced following the ψ-WGD. Functional enrichment analyses revealed that some genes, such as those encoding heat shock proteins, glycosyl hydrolase and cytochrome P450, expanded within the True-Araceae clade.

CONCLUSIONS

Our results improve our understanding of aroid phylogeny using the large number of single-/low-copy nuclear genes. In contrast to the Proto-Araceae group and the lemnoid clade adaption to aquatic environments, our analyses of WGD, diversification and functional enrichment indicated that WGD may play a more important role in the evolution of adaptations to tropical, terrestrial environments in the True-Araceae clade. These insights provide us with new resources to interpret the evolution of the Araceae.

Phylogenomics and the flowering plant tree of life

The advances accelerated by next-generation sequencing and long-read sequencing technologies continue to provide an impetus for plant phylogenetic study. In the past decade, a large number of phylogenetic studies adopting hundreds to thousands of genes across a wealth of clades have emerged and ushered plant phylogenetics and evolution into a new era. In the meantime, a roadmap for researchers when making decisions across different approaches for their phylogenomic research design is imminent. This review focuses on the utility of genomic data (from organelle genomes, to both reduced representation sequencing and whole-genome sequencing) in phylogenetic and evolutionary investigations, describes the baseline methodology of experimental and analytical procedures, and summarizes recent progress in flowering plant phylogenomics at the ordinal, familial, tribal, and lower levels. We also discuss the challenges, such as the adverse impact on orthology inference and phylogenetic reconstruction raised from systematic errors, and underlying biological factors, such as whole-genome duplication, hybridization/introgression, and incomplete lineage sorting, together suggesting that a bifurcating tree may not be the best model for the tree of life. Finally, we discuss promising avenues for future plant phylogenomic studies.

The nearly complete assembly of the Cercis chinensis genome and Fabaceae phylogenomic studies provide insights into new gene evolution

Fabaceae is a large family of angiosperms with high biodiversity that contains a variety of economically important crops and model plants for the study of biological nitrogen fixation. Polyploidization events have been extensively studied in some Fabaceae plants, but the occurrence of new genes is still concealed, owing to a lack of genomic information on certain species of the basal clade of Fabaceae. Cercis chinensis (Cercidoideae) is one such species; it diverged earliest from Fabaceae and is essential for phylogenomic studies and new gene predictions in Fabaceae. To facilitate genomic studies on Fabaceae, we performed genome sequencing of C. chinensis and obtained a 352.84 Mb genome, which was further assembled into seven pseudochromosomes with 30 612 predicted protein-coding genes. Compared with other legume genomes, that of C. chinensis exhibits no lineage-specific polyploidization event. Further phylogenomic analyses of 22 legumes and 11 other angiosperms revealed that many gene families are lineage specific before and after the diversification of Fabaceae. Among them, dozens of genes are candidates for new genes that have evolved from intergenic regions and are thus regarded as de novo-originated genes. They differ significantly from established genes in coding sequence length, exon number, guanine-cytosine content, and expression patterns among tissues. Functional analysis revealed that many new genes are related to asparagine metabolism. This study represents an important advance in understanding the evolutionary pattern of new genes in legumes and provides a valuable resource for plant phylogenomic studies.

Gain or Loss? Evidence for Legume Predisposition to Symbiotic Interactions with Rhizobia via Loss of Pathogen-Resistance-Related Gene Families

Nodulation is a hallmark yet non-universal characteristic of legumes. It is unknown whether the mechanisms underlying nitrogen-fixing symbioses evolved within legumes and the broader nitrogen-fixing clade (NFC) repeatedly de novo or based on common ancestral pathways. Ten new transcriptomes representing members from the Cercidoideae and Caesalpinioideae subfamilies were supplemented with published omics data from 65 angiosperms, to investigate how gene content correlates with nodulation capacity within Fabaceae and the NFC. Orthogroup analysis categorized annotated genes into 64150 orthogroups, of which 19 were significantly differentially represented between nodulating versus non-nodulating NFC species and were most commonly absent in nodulating taxa. The distribution of six over-represented orthogroups within Viridiplantae representatives suggested that genomic evolution events causing gene family expansions, including whole-genome duplications (WGDs), were unlikely to have facilitated the development of stable symbioses within Fabaceae as a whole. Instead, an absence of representation of 13 orthogroups indicated that losses of genes involved in trichome development, defense and wounding responses were strongly associated with rhizobial symbiosis in legumes. This finding provides novel evidence of a lineage-specific predisposition for the evolution and/or stabilization of nodulation in Fabaceae, in which a loss of pathogen resistance genes may have allowed for stable mutualistic interactions with rhizobia.

Comparative cytogenomics reveals genome reshuffling and centromere repositioning in the legume tribe Phaseoleae

The tribe Phaseoleae includes several legume crops with assembled genomes. Comparative genomic studies have evidenced the preservation of large genomic blocks among legumes, although chromosome dynamics during Phaseoleae evolution has not been investigated. We conducted a comparative genomic analysis to define an informative genomic block (GB) system and to reconstruct the ancestral Phaseoleae karyotype (APK). We identified GBs based on the orthologous genes between Phaseolus vulgaris and Vigna unguiculata and searched for GBs in different genomes of the Phaseolinae (P. lunatus) and Glycininae (Amphicarpaea edgeworthii) subtribes and Spatholobus suberectus (sister to Phaseolinae and Glycininae), using Medicago truncatula as the outgroup. We also used oligo-FISH probes of two P. vulgaris chromosomes to paint the orthologous chromosomes of two non-sequenced Phaseolinae species. We inferred the APK as having n = 11 and 19 GBs (A to S), hypothesizing five chromosome fusions that reduced the ancestral legume karyotype to n = 11. We identified the rearrangements among the APK and the subtribes and species, with extensive centromere repositioning in Phaseolus. We also reconstructed the chromosome number reduction in S. suberectus. The development of the GB system and the proposed APK provide useful approaches for future comparative genomic analyses of legume species.

Tree2GD: a phylogenomic method to detect large-scale gene duplication events

MOTIVATION

Whole-genome duplication events have long been discovered throughout the evolution of eukaryotes, contributing to genome complexity and biodiversity and leaving traces in the descending organisms. Therefore, an accurate and rapid phylogenomic method is needed to identify the retained duplicated genes on various lineages across the target taxonomy.

RESULTS

Here, we present Tree2GD, an integrated method to identify large-scale gene duplication events by automatically perform multiple procedures, including sequence alignment, recognition of homolog, gene tree/species tree reconciliation, Ks distribution of gene duplicates and synteny analyses. Application of Tree2GD on 2 datasets, 12 metazoan genomes and 68 angiosperms, successfully identifies all reported whole-genome duplication events exhibited by these species, showing effectiveness and efficiency of Tree2GD on phylogenomic analyses of large-scale gene duplications.

AVAILABILITY AND IMPLEMENTATION

Tree2GD is written in Python and C++ and is available at https://github.com/Dee-chen/Tree2gd.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Interplay of Ecological Opportunities and Functional Traits Drives the Evolution and Diversification of Millettiod Legumes (Fabaceae)

Understanding the striking diversity of the angiosperms is a paramount issue in biology and of interest to biologists. The Millettiod legumes is one of the most hyper-diverse groups of the legume family, containing many economically important medicine, furniture and craft species. In the present study, we explore how the interplay of past climate change, ecological opportunities and functional traits' evolution may have triggered diversification of the Millettiod legumes. Using a comprehensive species-level phylogeny from three plastid markers, we estimate divergence times, infer habit shifts, test the phylogenetic and temporal diversification heterogeneity, and reconstruct ancestral biogeographical ranges. We found that three dramatic accumulations of the Millettiod legumes occurred during the Miocene. The rapid diversification of the Millettiod legumes in the Miocene was driven by ecological opportunities created by the emergence of new niches and range expansion. Additionally, habit shifts and the switch between biomes might have facilitated the rapid diversification of the Millettiod legumes. The Millettiod legumes provide an excellent case for supporting the idea that the interplay of functional traits, biomes, and climatic and geographic factors drives evolutionary success.

Plastid phylogenomic insights into relationships, divergence, and evolution of Apiales

Members of Apiales are monophyletic and radiated in the Late Cretaceous. Fruit morphologies are critical for Apiales evolution and negative selection and mutation pressure play important roles in environmental adaptation. Apiales include many foods, spices, medicinal, and ornamental plants, but the phylogenetic relationships, origin and divergence, and adaptive evolution remain poorly understood. Here, we reconstructed Apiales phylogeny based on 72 plastid genes from 280 species plastid genomes representing six of seven families of this order. Highly supported phylogenetic relationships were detected, which revealed that each family of Apiales is monophyletic and confirmed that Pennanticeae is a member of Apiales. Genera Centella and Dickinsia are members of Apiaceae, and the genus Hydrocotyle previously classified into Apiaceae is confirmed to belong to Araliaceae. Besides, coalescent phylogenetic analysis and gene trees cluster revealed ten genes that can be used for distinguishing species among families of Apiales. Molecular dating suggested that the Apiales originated during the mid-Cretaceous (109.51 Ma), with the families' radiation occurring in the Late Cretaceous. Apiaceae species exhibit higher differentiation compared to other families. Ancestral trait reconstruction suggested that fruit morphological evolution may be related to shifts in plant types (herbaceous or woody), which in turn is related to the distribution areas and species numbers. Codon bias and positive selection analyses suggest that negative selection and mutation pressure may play important roles in environmental adaptation of Apiales members. Our results improve the phylogenetic framework of Apiales and provide insights into the origin, divergence, and adaptive evolution of this order and its members.

Phylogenomics as an effective approach to untangle cross-species hybridization event: A case study in the family Nymphaeaceae

Hybridization is common and considered as an important evolutionary force to increase intraspecific genetic diversity. Detecting hybridization events is crucial for understanding the evolutionary history of species and further improving molecular breeding. The studies on identifying hybridization events through the phylogenomic approach are still limited. We proposed the conception and method of identifying allopolyploidy events by phylogenomics. The reconciliation and summary of nuclear multi-labeled gene family trees were adopted to untangle hybridization events from next-generation data in our novel phylogenomic approach. Given horticulturalists’ relatively clear cultivated crossbreeding history, the water lily family is a suitable case for examining recent allopolyploidy events. Here, we reconstructed and confirmed the well-resolved nuclear phylogeny for the Nymphaeales family in the context of geological time as a framework for identifying hybridization signals. We successfully identified two possible allopolyploidy events with the parental lineages for the hybrids in the family Nymphaeaceae based on summarization from multi-labeled gene family trees of Nymphaeales. The lineages where species Nymphaea colorata and Nymphaea caerulea are located may be the progenitors of horticultural cultivated species Nymphaea ‘midnight’ and Nymphaea ‘Woods blue goddess’. The proposed hybridization hypothesis is also supported by horticultural breeding records. Our methodology can be widely applied to identify hybridization events and theoretically facilitate the genome breeding design of hybrid plants.

A genome resource for Acacia, Australia's largest plant genus

Acacia (Leguminosae, Caesalpinioideae, mimosoid clade) is the largest and most widespread genus of plants in the Australian flora, occupying and dominating a diverse range of environments, with an equally diverse range of forms. For a genus of its size and importance, Acacia currently has surprisingly few genomic resources. Acacia pycnantha, the golden wattle, is a woody shrub or tree occurring in south-eastern Australia and is the country's floral emblem. To assemble a genome for A. pycnantha, we generated long-read sequences using Oxford Nanopore Technology, 10x Genomics Chromium linked reads, and short-read Illumina sequences, and produced an assembly spanning 814 Mb, with a scaffold N50 of 2.8 Mb, and 98.3% of complete Embryophyta BUSCOs. Genome annotation predicted 47,624 protein-coding genes, with 62.3% of the genome predicted to comprise transposable elements. Evolutionary analyses indicated a shared genome duplication event in the Caesalpinioideae, and conflict in the relationships between Cercis (subfamily Cercidoideae) and subfamilies Caesalpinioideae and Papilionoideae (pea-flowered legumes). Comparative genomics identified a suite of expanded and contracted gene families in A. pycnantha, and these were annotated with both GO terms and KEGG functional categories. One expanded gene family of particular interest is involved in flowering time and may be associated with the characteristic synchronous flowering of Acacia. This genome assembly and annotation will be a valuable resource for all studies involving Acacia, including the evolution, conservation, breeding, invasiveness, and physiology of the genus, and for comparative studies of legumes.

The innovation of the symbiosome has enhanced the evolutionary stability of nitrogen fixation in legumes

Nitrogen-fixing symbiosis is globally important in ecosystem functioning and agriculture, yet the evolutionary history of nodulation remains the focus of considerable debate. Recent evidence suggesting a single origin of nodulation followed by massive parallel evolutionary losses raises questions about why a few lineages in the N₂ -fixing clade retained nodulation and diversified as stable nodulators, while most did not. Within legumes, nodulation is restricted to the two most diverse subfamilies, Papilionoideae and Caesalpinioideae, which show stable retention of nodulation across their core clades. We characterize two nodule anatomy types across 128 species in 56 of the 152 genera of the legume subfamily Caesalpinioideae: fixation thread nodules (FTs), where nitrogen-fixing bacteroids are retained within the apoplast in modified infection threads, and symbiosomes, where rhizobia are symplastically internalized in the host cell cytoplasm within membrane-bound symbiosomes (SYMs). Using a robust phylogenomic tree based on 997 genes from 147 Caesalpinioideae genera, we show that losses of nodulation are more prevalent in lineages with FTs than those with SYMs. We propose that evolution of the symbiosome allows for a more intimate and enduring symbiosis through tighter compartmentalization of their rhizobial microsymbionts, resulting in greater evolutionary stability of nodulation across this species-rich pantropical legume clade.

Unraveling Origin, History, Genetics, and Strategies for Accelerated Domestication and Diversification of Food Legumes

Domestication is a dynamic and ongoing process of transforming wild species into cultivated species by selecting desirable agricultural plant features to meet human needs such as taste, yield, storage, and cultivation practices. Human plant domestication began in the Fertile Crescent around 12,000 years ago and spread throughout the world, including China, Mesoamerica, the Andes and Near Oceania, Sub-Saharan Africa, and eastern North America. Indus valley civilizations have played a great role in the domestication of grain legumes. Crops, such as pigeon pea, black gram, green gram, lablab bean, moth bean, and horse gram, originated in the Indian subcontinent, and Neolithic archaeological records indicate that these crops were first domesticated by early civilizations in the region. The domestication and evolution of wild ancestors into today’s elite cultivars are important contributors to global food supply and agricultural crop improvement. In addition, food legumes contribute to food security by protecting human health and minimize climate change impacts. During the domestication process, legume crop species have undergone a severe genetic diversity loss, and only a very narrow range of variability is retained in the cultivars. Further reduction in genetic diversity occurred during seed dispersal and movement across the continents. In general, only a few traits, such as shattering resistance, seed dormancy loss, stem growth behavior, flowering–maturity period, and yield traits, have prominence in the domestication process across the species. Thus, identification and knowledge of domestication responsive loci were often useful in accelerating new species’ domestication. The genes and metabolic pathways responsible for the significant alterations that occurred as an outcome of domestication might aid in the quick domestication of novel crops. Further, recent advances in “omics” sciences, gene-editing technologies, and functional analysis will accelerate the domestication and crop improvement of new crop species without losing much genetic diversity. In this review, we have discussed about the origin, center of diversity, and seed movement of major food legumes, which will be useful in the exploration and utilization of genetic diversity in crop improvement. Further, we have discussed about the major genes/QTLs associated with the domestication syndrome in pulse crops and the future strategies to improve the food legume crops.

Phylogenomics Resolves the Phylogeny of Theaceae by Using Low-Copy and Multi-Copy Nuclear Gene Makers and Uncovers a Fast Radiation Event Contributing to Tea Plants Diversity

Simple Summary

The Theaceae includes more than 300 species of great morphological diversity and has immense economic, cultural, and ornamental values. However, the evolutionary history of this family remains elusive. We integrated 91 genomes and transcriptome datasets of Theaceae and successfully resolved the phylogeny of Theaceae including relatives of cultivated tea plants from both extensive low-copy and multi-copy nuclear gene markers. Bayes-based molecular dating revealed that the ancestor of the tea family originated slightly earlier than the K-Pg boundary (Mass extinction events including the extinction of dinosaurs) with early diversification of three tribes associated with the Early Eocene Climatic Optimum. Further speciation analysis suggested a sole significant diversification shift rate in the common ancestor of Camellia associated with the Mid-Miocene Climatic Optimum. Collectively, polyploidy events, and key morphological innovation characters, such as pericarp with seed coat hardening, could possibly contribute to the Theaceae species diversity.

Abstract

Tea is one of the three most popular nonalcoholic beverages globally and has extremely high economic and cultural value. Currently, the classification, taxonomy, and evolutionary history of the tea family are largely elusive, including phylogeny, divergence, speciation, and diversity. For understanding the evolutionary history and dynamics of species diversity in Theaceae, a robust phylogenetic framework based on 1785 low-copy and 79,103 multi-copy nuclear genes from 91 tea plant genomes and transcriptome datasets had been reconstructed. Our results maximumly supported that the tribes Stewartieae and Gordonieae are successive sister groups to the tribe Theeae from both coalescent and super matrix ML tree analyses. Moreover, in the most evolved tribe, Theeae, the monophyletic genera Pyrenaria, Apterosperma, and Polyspora are the successive sister groups of Camellia. We also yield a well-resolved relationship of Camellia, which contains the vast majority of Theaceae species richness. Molecular dating suggests that Theaceae originated in the late L-Cretaceous, with subsequent early radiation under the Early Eocene Climatic Optimal (EECO) for the three tribes. A diversification rate shift was detected in the common ancestors of Camellia with subsequent acceleration in speciation rate under the climate optimum in the early Miocene. These results provide a phylogenetic framework and new insights into factors that likely have contributed to the survival of Theaceae, especially a successful radiation event of genus Camellia members to subtropic/tropic regions. These novel findings will facilitate the efficient conservation and utilization of germplasm resources for breeding cultivated tea and oil-tea. Collectively, these results provide a foundation for further morphological and functional evolutionary analyses across Theaceae.

Vicariance Between Cercis siliquastrum L. and Ceratonia siliqua L. Unveiled by the Physical-Chemical Properties of the Leaves' Epicuticular Waxes

Classically, vicariant phenomena have been essentially identified on the basis of biogeographical and ecological data. Here, we report unequivocal evidences that demonstrate that a physical-chemical characterization of the epicuticular waxes of the surface of plant leaves represents a very powerful strategy to get rich insight into vicariant events. We found vicariant similarity between Cercis siliquastrum L. (family Fabaceae, subfamily Cercidoideae) and Ceratonia siliqua L. (family Fabaceae, subfamily Caesalpinoideae). Both taxa converge in the Mediterranean basin (C. siliquastrum on the north and C. siliqua across the south), in similar habitats (sclerophyll communities of maquis) and climatic profiles. These species are the current representation of their subfamilies in the Mediterranean basin, where they overlap. Because of this biogeographic and ecological similarity, the environmental pattern of both taxa was found to be very significant. The physical-chemical analysis performed on the epicuticular waxes of C. siliquastrum and C. siliqua leaves provided relevant data that confirm the functional proximity between them. A striking resemblance was found in the epicuticular waxes of the abaxial surfaces of C. siliquastrum and C. siliqua leaves in terms of the dominant chemical compounds (1-triacontanol (C30) and 1-octacosanol (C28), respectively), morphology (intricate network of randomly organized nanometer-thick and micrometer-long plates), wettability (superhydrophobic character, with water contact angle values of 167.5 ± 0.5° and 162 ± 3°, respectively), and optical properties (in both species the light reflectance/absorptance of the abaxial surface is significantly higher/lower than that of the adaxial surface, but the overall trend in reflectance is qualitatively similar). These results enable us to include for the first time C. siliqua in the vicariant process exhibited by C. canadensis L., C. griffithii L., and C. siliquastrum.

Characterization of the complete chloroplast genome sequences of six Dalbergia species and its comparative analysis in the subfamily of Papilionoideae (Fabaceae)

Dalbergia spp. are numerous and widely distributed in pantropical areas in Asia, Africa and America, and most of the species have important economic and ecological value as precious timber. In this study, we determined and characterized six complete chloroplast genomes of Dalbergia species (Dalbergia obtusifolia, D. hupeana, D. mimosoides, D. sissoo, D. hancei, D. balansae), which displayed the typical quadripartite structure of angiosperms. The sizes of the genomes ranged from 155,698 bp (D. hancei) to 156,419 bp (D. obtusifolia). The complete chloroplast genomes of Dalbergia include 37 tRNA genes, eight rRNA genes and 84 protein-coding genes. We analysed the sequence diversity of Dalberigia chloroplast genomes coupled with previous reports. The results showed 12 noncoding regions (rps16-accD, trnR-UCU-trnG-UCC, ndhE-ndhG, trnG-UCC-psbZ, rps8-rpl14, trnP-UGG-psaJ, ndhH-rps15, trnQ-UUG-rps16, trnS-GCU-psbI, rps12-clpP, psbA-trnK-UUU, trnK-UUU-intron), and four coding regions (rps16, ycf1, rps15 and ndhF) showed many nucleotide variations that could be used as potential molecular markers. Based on a site-specific model, we analysed the selective pressure of chloroplast genes in Dalbergia species. Twenty-two genes with positively selected sites were detected, involving the photosynthetic system (ndhC, adhD, ndhF, petB, psaA, psaB, psbB, psbC, psbK and rbcL), self-replication category of genes (rpoA, rpoC2, rps3, rps12 and rps18) and others (accD, ccsA, cemA, clpP, matK, ycf1 and ycf2). Additionally, we identified potential RNA editing sites that were relatively conserved in the genus Dalbergia. Furthermore, the comparative analysis of cp genomes of Dalbergieae species indicated that the boundary of IRs/SSC was highly variable, which resulted in the size variation of cp genomes. Finally, phylogenetic analysis showed an inferred phylogenetic tree of Papilionoideae species with high bootstrap support and suggested that Amorpheae was the sister of the clade Dalbergieae. Moreover, three genera of the Pterocarpus clade showed a nested evolutionary relationship. These complete cp genomes provided valuable information for understanding the genetic variation and phylogenetic relationship of Dalbergia species with their relatives.

Plastid and mitochondrial phylogenomics reveal correlated substitution rate variation in Koenigia (Polygonoideae, Polygonaceae) and a reduced plastome for Koenigia delicatula including loss of all ndh genes

Koenigia, a genus proposed by Linnaeus, has a contentious taxonomic history. In particular, relationships among species and the circumscription of the genus relative to Aconogonon remain uncertain. To explore phylogenetic relationships of Koenigia with other members of tribe Persicarieae and to establish the timing of major evolutionary diversification events, genome skimming of organellar sequences was used to assemble plastomes and mitochondrial genes from 15 individuals representing 13 species. Most Persicarieae plastomes exhibit a conserved structure and content relative to other flowering plants. However, Koenigia delicatula has lost functional copies of all ndh genes and the intron from atpF. In addition, the rpl32 gene was relocated in the K. delicatula plastome, which likely occurred via overlapping inversions or differential expansion and contraction of the inverted repeat. The highly supported but conflicting relationships between plastome and mitochondrial trees and among gene trees complicates the circumscription of Koenigia, which could be caused by rapid diversification within a short period. Moreover, the plastome and mitochondrial trees revealed correlated variation in substitution rates among Persicarieae species, suggesting a shared underlying mechanism promoting evolutionary rate variation in both organellar genomes. The divergence of dwarf K. delicatula from other Koenigia species may be associated with the well-known Eocene Thermal Maximum 2 or Early Eocene Climatic Optimum event, while diversification of the core-Koenigia clade associates with the Mid-Miocene Climatic Optimum and the uplift of Qinghai-Tibetan Plateau and adjacent areas.

Diversity and phylogeny of the bacterial strains isolated from nodules of fenugreek (Trigonella foenum-graecum L.) in Iran

The diversity of fenugreek (Trigonella foenum-graecum L.) microsymbionts has been barely studied even though it is of great interest for being a spice and a medicinal plant. Here we analyzed fifty-nine bacterial strains isolated from fenugreek nodules originating from different geographic and climatic areas of Iran. Most of these strains exhibit phenotypic characteristics compatible with rhizobia and they nodulate fenugreek. Analysis of the recA and atpD genes shows that representative strains of ERIC-BOX-PCR groups cluster with the type strains of Ensifer meliloti and Ensifer kummerowiae as well as with strains capable of nodulating different Trigonella species found in other countries. The closeness of E. meliloti and E. kummerowiae suggests there is a need to revise the taxonomic status of the latter species. The nodC gene analysis shows that most Trigonella-nodulating strains belong to the symbiovar meliloti except those nodulating Trigonella arcuata in China, which belong to the symbiovar rigiduloides. This analysis shows that the type strains of E. kummerowiae, E. meliloti and E. medicae belonged to three well-defined groups within the symbiovar meliloti, with the Iranian strains belonging to the E. kummerowiae subgroup. The small group of strains unable to nodulate fenugreek isolated in this study belong to Enterobacter cloacae, reported for the first time as being a possible endophyte of fenugreek nodules.

Innovation and appropriation in mycorrhizal and rhizobial Symbioses

Most land plants benefit from endosymbiotic interactions with mycorrhizal fungi, including legumes and some nonlegumes that also interact with endosymbiotic nitrogen (N)-fixing bacteria to form nodules. In addition to these helpful interactions, plants are continuously exposed to would-be pathogenic microbes: discriminating between friends and foes is a major determinant of plant survival. Recent breakthroughs have revealed how some key signals from pathogens and symbionts are distinguished. Once this checkpoint has been passed and a compatible symbiont is recognized, the plant coordinates the sequential development of two types of specialized structures in the host. The first serves to mediate infection, and the second, which appears later, serves as sophisticated intracellular nutrient exchange interfaces. The overlap in both the signaling pathways and downstream infection components of these symbioses reflects their evolutionary relatedness and the common requirements of these two interactions. However, the different outputs of the symbioses, phosphate uptake versus N fixation, require fundamentally different components and physical environments and necessitated the recruitment of different master regulators, NODULE INCEPTION-LIKE PROTEINS, and PHOSPHATE STARVATION RESPONSES, for nodulation and mycorrhization, respectively.

A Reappraisal of Polyploidy Events in Grasses (Poaceae) in a Rapidly Changing World

Around 80% of megaflora species became extinct at the Cretaceous–Paleogene (K–Pg) boundary. Subsequent polyploidy events drove the survival of thousands of plant species and played a significant historical role in the development of the most successful modern cereal crops. However, current and rapid global temperature change poses an urgent threat to food crops worldwide, including the world’s big three cereals: rice, wheat, and maize, which are members of the grass family, Poaceae. Some minor cereals from the same family (such as teff) have grown in popularity in recent years, but there are important knowledge gaps regarding the similarities and differences between major and minor crops, including how polyploidy affects their biological processes under natural and (a)biotic stress conditions and thus the potential to harness polyploidization attributes for improving crop climate resilience. This review focuses on the impact of polyploidy events on the Poaceae family, which includes the world’s most important food sources, and discusses the past, present, and future of polyploidy research for major and minor crops. The increasing accessibility to genomes of grasses and their wild progenitors together with new tools and interdisciplinary research on polyploidy can support crop improvement for global food security in the face of climate change.