Phosphorus-Use-Efficiency Gene Identification in Fabaceae and RSL2 Expansion in Lupinus albus Is Associated with Low-Phosphorus Adaptation

Phosphorus is critical for plant growth but often becomes less accessible due to its precipitation with cations in soil. Fabaceae, a diverse plant family, exhibits robust adaptability and includes species like Lupinus albus, known for its efficient phosphorus utilization via cluster roots. Here, we systematically identified phosphorus-utilization-efficiency (PUE) gene families across 35 Fabaceae species, highlighting significant gene amplification in PUE pathways in Fabaceae. Different PUE pathways exhibited variable amplification, evolution, and retention patterns among various Fabaceae crops. Additionally, the number of homologous genes of the root hair development gene RSL2 in L. albus was far more than that in other Fabaceae species. Multiple copies of the RSL2 gene were amplified and retained in L. albus after whole genome triplication. The gene structure and motifs specifically retained in L. albus were different from homologous genes in other plants. Combining transcriptome analysis under low-phosphorus treatment, it was found that most of the homologous genes of RSL2 in L. albus showed high expression in the cluster roots, suggesting that the RSL2 gene family plays an important role in the adaptation process of L. albus to low-phosphorus environments and the formation of cluster roots.

The high-quality genome of Grona styracifolia uncovers the genomic mechanism of high levels of schaftoside, a promising drug candidate for treatment of COVID-19

Recent study has evidenced that traditional Chinese medicinal (TCM) plant-derived schaftoside shows promise as a potential drug candidate for COVID-19 treatment. However, the biosynthetic pathway of schaftoside in TCM plants remains unknown. In this study, the genome of the TCM herb Grona styracifolia (Osbeck) H.Ohashi & K.Ohashi (GSO), which is rich in schaftoside, was sequenced, and a high-quality assembly of GSO genome was obtained. Our findings revealed that GSO did not undergo recent whole genome duplication (WGD) but shared an ancestral papilionoid polyploidy event, leading to the gene expansion of chalcone synthase (CHS) and isoflavone 2'-hydroxylase (HIDH). Furthermore, GSO-specific tandem gene duplication resulted in the gene expansion of C-glucosyltransferase (CGT). Integrative analysis of the metabolome and transcriptome identified 13 CGTs and eight HIDHs involved in the biosynthetic pathway of schaftoside. Functional studies indicated that CGTs and HIDHs identified here are bona fide responsible for the biosynthesis of schaftoside in GSO, as confirmed through hairy root transgenic system and in vitro enzyme activity assay. Taken together, the ancestral papilionoid polyploidy event expanding CHSs and HIDHs, along with the GSO-specific tandem duplication of CGT, contributes, partially if not completely, to the robust biosynthesis of schaftoside in GSO. These findings provide insights into the genomic mechanisms underlying the abundant biosynthesis of schaftoside in GSO, highlighting the potential of GSO as a source of bioactive compounds for pharmaceutical development.

OROSOMUCOID PROTEIN 1 regulation of sphingolipid synthesis is required for nodulation in Aeschynomene evenia

Legumes establish symbiotic interactions with nitrogen-fixing rhizobia that are accommodated in root-derived organs known as nodules. Rhizobial recognition triggers a plant symbiotic signaling pathway that activates 2 coordinated processes: infection and nodule organogenesis. How these processes are orchestrated in legume species utilizing intercellular infection and lateral root base nodulation remains elusive. Here, we show that Aeschynomene evenia OROSOMUCOID PROTEIN 1 (AeORM1), a key regulator of sphingolipid biosynthesis, is required for nodule formation. Using A. evenia orm1 mutants, we demonstrate that alterations in AeORM1 function trigger numerous early aborted nodules, defense-like reactions, and shorter lateral roots. Accordingly, AeORM1 is expressed during lateral root initiation and elongation, including at lateral root bases where nodule primordium form in the presence of symbiotic bradyrhizobia. Sphingolipidomics revealed that mutations in AeORM1 lead to sphingolipid overaccumulation in roots relative to the wild type, particularly for very long-chain fatty acid-containing ceramides. Taken together, our findings reveal that AeORM1-regulated sphingolipid homeostasis is essential for rhizobial infection and nodule organogenesis, as well as for lateral root development in A. evenia.

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.

A reference-grade genome assembly for Astragalus mongholicus and insights into the biosynthesis and high accumulation of triterpenoids and flavonoids in its roots

Astragalus membranaceus var. mongholicus (AMM), a member of the Leguminosae, is one of the most important medicinal plants worldwide. The dried roots of AMM have a wide range of pharmacological effects and are a traditional Chinese medicine. Here, we report the first chromosome-level reference genome of AMM, comprising nine pseudochromosomes with a total size of 1.47 Gb and 27 868 protein-encoding genes. Comparative genomic analysis reveals that AMM has not experienced an independent whole-genome duplication (WGD) event after the WGD event shared by the Papilionoideae species. Analysis of long terminal repeat retrotransposons suggests a recent burst of these elements at approximately 0.13 million years ago, which may explain the large size of the AMM genome. Multiple gene families involved in the biosynthesis of triterpenoids and flavonoids were expanded, and our data indicate that tandem duplication has been the main driver for expansion of these families. Among the expanded families, the phenylalanine ammonia-lyase gene family was primarily expressed in the roots of AMM, suggesting their roles in the biosynthesis of phenylpropanoid compounds. The functional versatility of 2,3-oxidosqualene cyclase genes in cluster III may play a critical role in the diversification of triterpenoids in AMM. Our findings provide novel insights into triterpenoid and flavonoid biosynthesis and can facilitate future research on the genetics and medical applications of AMM.

Genes Associated with Biological Nitrogen Fixation Efficiency Identified Using RNA Sequencing in Red Clover (Trifolium pratense L.)

Commonly studied in the context of legume-rhizobia symbiosis, biological nitrogen fixation (BNF) is a key component of the nitrogen cycle in nature. Despite its potential in plant breeding and many years of research, information is still lacking as to the regulation of hundreds of genes connected with plant-bacteria interaction, nodulation, and nitrogen fixation. Here, we compared root nodule transcriptomes of red clover (Trifolium pratense L.) genotypes with contrasting nitrogen fixation efficiency, and we found 491 differentially expressed genes (DEGs) between plants with high and low BNF efficiency. The annotation of genes expressed in nodules revealed more than 800 genes not yet experimentally confirmed. Among genes mediating nodule development, four nod-ule-specific cysteine-rich (NCR) peptides were confirmed in the nodule transcriptome. Gene duplication analyses revealed that genes originating from tandem and dispersed duplication are significantly over-represented among DEGs. Weighted correlation network analysis (WGCNA) organized expression profiles of the transcripts into 16 modules linked to the analyzed traits, such as nitrogen fixation efficiency or sample-specific modules. Overall, the results obtained broaden our knowledge about transcriptomic landscapes of red clover's root nodules and shift the phenotypic description of BNF efficiency on the level of gene expression in situ.

Genome wide association study to detect genetic regions related to isoflavone content in a mutant soybean population derived from radiation breeding

Isoflavones are major secondary metabolites that are exclusively produced by legumes, including soybean. Soy isoflavones play important roles in human health as well as in the plant defense system. The isoflavone content is influenced by minor-effect quantitative trait loci, which interact with polygenetic and environmental factors. It has been difficult to clarify the regulation of isoflavone biosynthesis because of its complex heritability and the influence of external factors. Here, using a genotype-by-sequencing-based genome-wide association mapping study, 189 mutant soybean genotypes (the mutant diversity pool, MDP) were genotyped on the basis of 25,646 high-quality single nucleotide polymorphisms (SNPs) with minor allele frequency of >0.01 except for missing data. All the accessions were phenotyped by determining the contents of 12 isoflavones in the soybean seeds in two consecutive years (2020 and 2021). Then, quantitative trait nucleotides (QTNs) related to isoflavone contents were identified and validated using multi-locus GWAS models. A total of 112 and 46 QTNs related to isoflavone contents were detected by multiple MLM-based models in 2020 and 2021, respectively. Of these, 12 and 5 QTNs were related to more than two types of isoflavones in 2020 and 2021, respectively. Forty-four QTNs were detected within the 441-Kb physical interval surrounding Gm05:38940662. Of them, four QTNs (Gm05:38936166, Gm05:38936167, Gm05:38940662, and Gm05:38940717) were located at Glyma.05g206900 and Glyma.05g207000, which encode glutathione S-transferase THETA 1 (GmGSTT1), as determined from previous quantitative trait loci annotations and the literature. We detected substantial differences in the transcript levels of GmGSTT1 and two other core genes (IFS1 and IFS2) in the isoflavone biosynthetic pathway between the original cultivar and its mutant. The results of this study provide new information about the factors affecting isoflavone contents in soybean seeds and will be useful for breeding soybean lines with high and stable concentrations of isoflavones.

The chromosome-level genome assembly of Astragalus sinicus and comparative genomic analyses provide new resources and insights for understanding legume-rhizobial interactions

The legume species Astragalus sinicus (Chinese milk vetch [CMV]) has been widely cultivated for centuries in southern China as one of the most important green manures/cover crops for improving rice productivity and preventing soil degeneration. In this study, we generated the first chromosome-scale reference genome of CMV by combining PacBio and Illumina sequencing with high-throughput chromatin conformation capture (Hi-C) technology. The CMV genome was 595.52 Mb in length, with a contig N50 size of 1.50 Mb. Long terminal repeats (LTRs) had been amplified and contributed to genome size expansion in CMV. CMV has undergone two whole-genome duplication (WGD) events, and the genes retained after the WGD shared by Papilionoideae species shaped the rhizobial symbiosis and the hormonal regulation of nodulation. The chalcone synthase (CHS) gene family was expanded and was expressed primarily in the roots of CMV. Intriguingly, we found that resistance genes were more highly expressed in roots than in nodules of legume species, suggesting that their expression may be increased to bolster plant immunity in roots to cope with pathogen infection in legumes. Our work sheds light on the genetic basis of nodulation and symbiosis in CMV and provides a benchmark for accelerating genetic research and molecular breeding in the future.

Differential light-dependent regulation of soybean nodulation by papilionoid-specific HY5 homologs

Legumes have evolved photosynthesis and symbiotic nitrogen fixation for the acquisition of energy and nitrogen nutrients. During the transition from heterotrophic to autotrophic growth, blue light primarily triggers photosynthesis and low soil nitrogen induces symbiotic nodulation. Whether and how darkness and blue light influence root symbiotic nodulation during this transition is unknown. Here, we show that short-term darkness promotes nodulation and that blue light inhibits nodulation through two soybean TGACG-motif-binding factors (STF1 and STF2), which are Papilionoideae-specific transcription factors and divergent orthologs of Arabidopsis ELONGATED HYPOCOTYL 5 (HY5). STF1 and STF2 negatively regulate soybean nodulation by repressing the transcription of nodule inception a (GmNINa), which is a central regulator of nodulation, in response to darkness and blue light. STF1 and STF2 are not capable of moving from the shoots to roots, and they act both locally and systemically to mediate darkness- and blue-light-regulated nodulation. We further show that cryptochromes GmCRY1s are required for nodulation in the dark and partially contribute to the blue light inhibition of nodulation. In addition, root GmCRY1s mediate blue-light-induced transcription of STF1 and STF2, and intriguingly, GmCRY1b can interact with STF1 and STF2 to stabilize the protein stability of STF1 and STF2. Our results establish that the blue light receptor GmCRY1s-STF1/2 module plays a pivotal role in integrating darkness/blue light and nodulation signals. Furthermore, our findings reveal a molecular basis by which photosensory pathways modulate nodulation and autotrophic growth through an intricate interplay facilitating seedling establishment in response to low nitrogen and light signals.

Identification of the Important Genes of Bradyrhizobium diazoefficiens 113-2 Involved in Soybean Nodule Development and Senescence

Legume nodule development and senescence directly affect nitrogen fixation efficiency and involve a programmed series of molecular events. These molecular events are carried out synchronously by legumes and rhizobia. The characteristics and molecular mechanisms of nitrogen fixation at soybean important developmental stages play critical roles in soybean cultivation and fertilizer application. Although the gene expression of soybean were analyzed in nodules at five important soybean developmental stages, information on the expression of rhizobial genes in these nodule samples is limited. In the present study, we investigated the expression of Bradyrhizobium diazoefficiens 113-2 genes in the nodule samples from five developmental stages of soybean (Branching stage, flowering stage, fruiting stage, pod stage and harvest stage). Similar gene expression patterns of B. diazoefficiens 113-2 were existed during optimal symbiotic functioning, while different expression patterns were found among early nodule development, nitrogen fixation progress and nodule senescence. Besides, we identified 164 important different expression genes (DEGs) associated with nodule development and senescence. These DEGs included those encoding nod, nif, fix proteins and T3SS secretion system-related proteins, as well as proteins involved in nitrogen metabolism, ABC transporters and two-component system pathways. Gene Ontology, KEGG pathway and homology analysis of the identified DEGs revealed that most of these DEGs are uncharacterized genes associated with nodule development and senescence, and they are not core genes among the rhizobia genomes. Our results provide new clues for the understanding of the genetic determinants of soil rhizobia in nodule development and senescence, and supply theoretical basis for the creation of high efficiency soybean cultivation technology.

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

Fabaceae are the third largest angiosperm family, with 765 genera and ∼19 500 species. They are important both economically and ecologically, and global Fabaceae crops are intensively studied in part for their nitrogen-fixing ability. However, resolution of the intrasubfamilial Fabaceae phylogeny and divergence times has remained elusive, precluding a reconstruction of the evolutionary history of symbiotic nitrogen fixation in Fabaceae. Here, we report a highly resolved phylogeny using >1500 nuclear genes from newly sequenced transcriptomes and genomes of 391 species, along with other datasets, for a total of 463 legumes spanning all 6 subfamilies and 333 of 765 genera. The subfamilies are maximally supported as monophyletic. The clade comprising subfamilies Cercidoideae and Detarioideae is sister to the remaining legumes, and Duparquetioideae and Dialioideae are successive sisters to the clade of Papilionoideae and Caesalpinioideae. Molecular clock estimation revealed an early radiation of subfamilies near the K/Pg boundary, marked by mass extinction, and subsequent divergence of most tribe-level clades within ∼15 million years. Phylogenomic analyses of thousands of gene families support 28 proposed putative whole-genome duplication/whole-genome triplication events across Fabaceae, including those at the ancestors of Fabaceae and five of the subfamilies, and further analyses supported the Fabaceae ancestral polyploidy. The evolution of rhizobial nitrogen-fixing nodulation in Fabaceae was probed by ancestral character reconstruction and phylogenetic analyses of related gene families and the results support the hypotheses of one or two switch(es) to rhizobial nodulation followed by multiple losses. Collectively, these results provide a foundation for further morphological and functional evolutionary analyses across Fabaceae.

Genome-wide survey of soybean papain-like cysteine proteases and their expression analysis in root nodule symbiosis

BACKGROUND

Plant papain-like cysteine proteases (PLCPs) are a large class of proteolytic enzymes and play important roles in root nodule symbiosis (RNS), while the whole-genome studies of PLCP family genes in legume are quite limited, and the roles of Glycine max PLCPs (GmPLCPs) in nodulation, nodule development and senescence are not fully understood.

RESULTS

In the present study, we identified 97 GmPLCPs and performed a genome-wide survey to explore the expansion of soybean PLCP family genes and their relationships to RNS. Nineteen paralogous pairs of genomic segments, consisting of 77 GmPLCPs, formed by whole-genome duplication (WGD) events were identified, showing a high degree of complexity in duplication. Phylogenetic analysis among different species showed that the lineage differentiation of GmPLCPs occurred after family expansion, and large tandem repeat segment were specifically in soybean. The expression patterns of GmPLCPs in symbiosis-related tissues and nodules identified RNS-related GmPLCPs and provided insights into their putative symbiotic functions in soybean. The symbiotic function analyses showed that a RNS-related GmPLCP gene (Glyma.04G190700) really participate in nodulation and nodule development.

CONCLUSIONS

Our findings improved our understanding of the functional diversity of legume PLCP family genes, and provided insights into the putative roles of the legume PLCPs in nodulation, nodule development and senescence.

Genome-Wide Identification, Characterization, and Regulation of RWP-RK Gene Family in the Nitrogen-Fixing Clade

RWP-RK is a plant-specific family of transcription factors, involved in nitrate response, gametogenesis, and nodulation. However, genome-wide characterization, phylogeny, and the regulation of RWP-RK genes in the nodulating and non-nodulating plant species of nitrogen-fixing clade (NFC) are widely unknown. Therefore, we identified a total of 292 RWP-RKs, including 278 RWP-RKs from 25 NFC species and 14 RWP-RKs from the outgroup, Arabidopsis thaliana. We classified the 292 RWP-RKs in two subfamilies: the NIN-like proteins (NLPs) and the RWP-RK domain proteins (RKDs). The transcriptome and phylogenetic analysis of RWP-RKs suggested that, compared to RKD genes, the NLP genes were just upregulated in nitrate response and nodulation. Moreover, nodule-specific NLP genes of some nodulating NFC species may have a common ancestor (OG0002084) with AtNLP genes in A. thaliana. Further, co-expression networks of A.thaliana under N-starvation and N-supplementation conditions revealed that there is a higher correlation between expression of AtNLP genes and symbiotic genes during N-starvation. In P. vulgaris, we confirmed that N-starvation stimulated nodulation by regulating expression of PvNLP2, closely related to AtNLP6 and AtNLP7 with another common origin (OG0004041). Taken together, we concluded that different origins of the NLP genes involved in both N-starvation response and specific expression of nodulation would contribute to the evolution of nodulation in NFC plant species. Our results shed light on the phylogenetic relationships of NLP genes and their differential regulation in nitrate response of A. thaliana and nodulation of NFC.

Polyploid plants obtain greater fitness benefits from a nutrient acquisition mutualism

Polyploidy is a key driver of ecological and evolutionary processes in plants, yet little is known about its effects on biotic interactions. This gap in knowledge is especially profound for nutrient acquisition mutualisms, despite the fact that they regulate global nutrient cycles and structure ecosystems. Generalism in mutualistic interactions depends on the range of potential partners (niche breadth), the benefits obtained and ability to maintain benefits across a variety of partners (fitness plasticity). Here, we determine how each of these is influenced by polyploidy in the legume-rhizobium mutualism. We inoculated a broad geographic sample of natural diploid and autotetraploid alfalfa (Medicago sativa) lineages with a diverse panel of Sinorhizobium bacterial symbionts. To analyze the extent and mechanism of generalism, we measured host growth benefits and functional traits. Autotetraploid plants obtained greater fitness enhancement from mutualistic interactions and were better able to maintain this across diverse rhizobial partners (i.e. low plasticity in fitness) relative to diploids. These benefits were not attributed to increases in niche breadth, but instead reflect increased rewards from investment in the mutualism. Polyploid plants displayed greater generalization in bacterial mutualisms relative to diploids, illustrating another axis of advantage for polyploids over diploids.

Autopolyploidy alters nodule-level interactions in the legume-rhizobium mutualism

PREMISE

Polyploidy is a major genetic driver of ecological and evolutionary processes in plants, yet its effects on plant interactions with mutualistic microbes remain unresolved. The legume-rhizobium symbiosis regulates global nutrient cycles and plays a role in the diversification of legume species. In this mutualism, rhizobia bacteria fix nitrogen in exchange for carbon provided by legume hosts. This exchange occurs inside root nodules, which house bacterial cells and represent the interface of legume-rhizobium interactions. Although polyploidy may directly impact the legume-rhizobium mutualism, no studies have explored how it alters the internal structure of nodules.

METHODS

We created synthetic autotetraploids using Medicago sativa subsp. caerulea. Neotetraploid plants and their diploid progenitors were singly inoculated with two strains of rhizobia, Sinorhizobium meliloti and S. medicae. Confocal microscopy was used to quantify internal traits of nodules produced by diploid and neotetraploid plants.

RESULTS

Autotetraploid plants produced larger nodules with larger nitrogen fixation zones than diploids for both strains of rhizobia, although the significance of these differences was limited by power. Neotetraploid M. sativa subsp. caerulea plants also produced symbiosomes that were significantly larger, nearly twice the size, than those present in diploids.

CONCLUSIONS

This study sheds light on how polyploidy directly affects a plant-bacterium mutualism and uncovers novel mechanisms. Changes in plant-microbe interactions that directly result from polyploidy likely contribute to the increased ability of polyploid legumes to establish in diverse environments.

A global coexpression network of soybean genes gives insights into the evolution of nodulation in nonlegumes and legumes

A coexpression network is a powerful tool for revealing genes' relationship with many biological processes. Mass transcriptomic and genomic data from different plant species provide the foundation for understanding the evolution of nodulation across the Viridiplantae at a systematic level. We used weighted coexpression network analysis (WGCNA) to mine a nodule-related module (NRM) in Glycine max. Comparative genomic analysis of 78 green plant species revealed that NRM genes are recruited from different evolutionary nodes along with gene duplication events. A set of core coexpressed genes within legumes may play vital roles in regulating nodule environments essential for nitrogen fixation, including oxygen concentrations, sulfur transport, and iron homeostasis (such as GmCHY). The regulation of these genes occurred mainly at the transcription level, although some of them, such as sulfate transporters, may also undergo positive selection at protein level. We revealed that ancient orthologs and duplication events before the origin of legumes were preadapted for symbiosis. Conserved coregulated genes found within legumes paved the way for nodule formation and nitrogen fixation. These findings provide significant insights into the evolution of nodulation and indicate promising candidates for identifying other key components of legume nodulation and nitrogen fixation.

Bradyrhizobium diazoefficiens USDA 110- Glycine max Interactome Provides Candidate Proteins Associated with Symbiosis

Although the legume-rhizobium symbiosis is a most-important biological process, there is a limited knowledge about the protein interaction network between host and symbiont. Using interolog- and domain-based approaches, we constructed an interspecies protein interactome containing 5115 protein-protein interactions between 2291 Glycine max and 290 Bradyrhizobium diazoefficiens USDA 110 proteins. The interactome was further validated by the expression pattern analysis in nodules, gene ontology term semantic similarity, co-expression analysis, and luciferase complementation image assay. In the G. max-B. diazoefficiens interactome, bacterial proteins are mainly ion channel and transporters of carbohydrates and cations, while G. max proteins are mainly involved in the processes of metabolism, signal transduction, and transport. We also identified the top 10 highly interacting proteins (hubs) for each species. Kyoto Encyclopedia of Genes and Genomes pathway analysis for each hub showed that a pair of 14-3-3 proteins (SGF14g and SGF14k) and 5 heat shock proteins in G. max are possibly involved in symbiosis, and 10 hubs in B. diazoefficiens may be important symbiotic effectors. Subnetwork analysis showed that 18 symbiosis-related soluble N-ethylmaleimide sensitive factor attachment protein receptor proteins may play roles in regulating bacterial ion channels, and SGF14g and SGF14k possibly regulate the rhizobium dicarboxylate transport protein DctA. The predicted interactome provide a valuable basis for understanding the molecular mechanism of nodulation in soybean.

The direct effects of plant polyploidy on the legume-rhizobia mutualism

BACKGROUND

Polyploidy is known to significantly alter plant genomes, phenotypes and interactions with the abiotic environment, yet the impacts of polyploidy on plant-biotic interactions are less well known. A particularly important plant-biotic interaction is the legume-rhizobia mutualism, in which rhizobia fix atmospheric nitrogen in exchange for carbon provided by legume hosts. This mutualism regulates nutrient cycles in natural ecosystems and provides nitrogen to agricultural environments. Despite the ecological, evolutionary and agricultural importance of plant polyploidy and the legume-rhizobia mutualism, it is not yet fully understood whether plant polyploidy directly alters mutualism traits or the consequences on plant growth.

SCOPE

The aim was to propose a conceptual framework to understand how polyploidy might directly enhance the quantity and quality of rhizobial symbionts hosted by legume plants, resulting in increased host access to fixed nitrogen (N). Mechanistic hypotheses have been devised to examine how polyploidy can directly alter traits that impact the quantity (e.g. nodule number, nodule size, terminal bacteroid differentiation) and quality of symbionts (e.g. nodule environment, partner choice, host sanctions). To evaluate these hypotheses, an exhaustive review of studies testing the effects of plant polyploidy on the mutualism was conducted. In doing so, overall trends were synthesized, highlighting the limited understanding of the mechanisms that underlie variation in results achieved thus far, revealing striking gaps in knowledge and uncovering areas ripe for future research.

CONCLUSIONS

Plant polyploidy can immediately alter nodule size, N fixation rate and the identity of rhizobial symbionts hosted by polyploid legumes, but many of the mechanistic hypotheses proposed here, such as bacteroid number and enhancements of the nodule environment, remain unexplored. Although current evidence supports a role of plant polyploidy in enhancing key aspects of the legume-rhizobia mutualism, the underlying mechanisms and effects on host benefit from the mutualism remain unresolved.

Non-Additive Transcriptomic Responses to Inoculation with Rhizobia in a Young Allopolyploid Compared with Its Diploid Progenitors

Root nodule symbioses (nodulation) and whole genome duplication (WGD, polyploidy) are both important phenomena in the legume family (Leguminosae). Recently, it has been proposed that polyploidy may have played a critical role in the origin or refinement of nodulation. However, while nodulation and polyploidy have been studied independently, there have been no direct studies of mechanisms affecting the interactions between these phenomena in symbiotic, nodule-forming species. Here, we examined the transcriptome-level responses to inoculation in the young allopolyploid Glycine dolichocarpa (T2) and its diploid progenitor species to identify underlying processes leading to the enhanced nodulation responses previously identified in T2. We assessed the differential expression of genes and, using weighted gene co-expression network analysis (WGCNA), identified modules associated with nodulation and compared their expression between species. These transcriptomic analyses revealed patterns of non-additive expression in T2, with evidence of transcriptional responses to inoculation that were distinct from one or both progenitors. These differential responses elucidate mechanisms underlying the nodulation-related differences observed between T2 and the diploid progenitors. Our results indicate that T2 has reduced stress-related transcription, coupled with enhanced transcription of modules and genes implicated in hormonal signaling, both of which are important for nodulation.

Evolutionary History of Subtilases in Land Plants and Their Involvement in Symbiotic Interactions

Subtilases, a family of proteases involved in a variety of developmental processes in land plants, are also involved in both mutualistic symbiosis and host-pathogen interactions in different angiosperm lineages. We examined the evolutionary history of subtilase genes across land plants through a phylogenetic analysis integrating amino acid sequence data from full genomes, transcriptomes, and characterized subtilases of 341 species of diverse green algae and land plants along with subtilases from 12 species of other eukaryotes, archaea, and bacteria. Our analysis reconstructs the subtilase gene phylogeny and identifies 11 new gene lineages, six of which have no previously characterized members. Two large, previously unnamed, subtilase gene lineages that diverged before the origin of angiosperms accounted for the majority of subtilases shown to be associated with symbiotic interactions. These lineages expanded through both whole-genome and tandem duplication, with differential neofunctionalization and subfunctionalization creating paralogs associated with different symbioses, including nodulation with nitrogen-fixing bacteria, arbuscular mycorrhizae, and pathogenesis in different plant clades. This study demonstrates for the first time that a key gene family involved in plant-microbe interactions proliferated in size and functional diversity before the explosive radiation of angiosperms.

Function and evolution of a Lotus japonicus AP2/ERF family transcription factor that is required for development of infection threads

Legume-rhizobium symbiosis is achieved by two major events evolutionarily acquired: root hair infection and organogenesis. Infection thread (IT) development is a distinct element for rhizobial infection. Through ITs, rhizobia are efficiently transported from infection foci on root hairs to dividing meristematic cortical cells. To unveil this process, we performed genetic screening using Lotus japonicus MG-20 and isolated symbiotic mutant lines affecting nodulation, root hair morphology, and IT development. Map-based cloning identified an AP2/ERF transcription factor gene orthologous to Medicago truncatula ERN1. LjERN1 was activated in response to rhizobial infection and depended on CYCLOPS and NSP2. Legumes conserve an ERN1 homolog, ERN2, that functions redundantly with ERN1 in M. truncatula. Phylogenetic analysis showed that the lineages of ERN1 and ERN2 genes originated from a gene duplication event in the common ancestor of legume plants. However, genomic analysis suggested the lack of ERN2 gene in the L. japonicus genome, consistent with Ljern1 mutants exhibited a root hair phenotype that is observed in ern1/ern2 double mutants in M. truncatula. Molecular evolutionary analysis suggested that the nonsynonymous/synonymous rate ratios of legume ERN1 genes was almost identical to that of non-legume plants, whereas the ERN2 genes experienced a relaxed selective constraint.

Comparative transcriptome analyses of flower development in four species of Achimenes (Gesneriaceae)

BACKGROUND

Flowers have an amazingly diverse display of colors and shapes, and these characteristics often vary significantly among closely related species. The evolution of diverse floral form can be thought of as an adaptive response to pollination and reproduction, but it can also be seen through the lens of morphological and developmental constraints. To explore these interactions, we use RNA-seq across species and development to investigate gene expression and sequence evolution as they relate to the evolution of the diverse flowers in a group of Neotropical plants native to Mexico-magic flowers (Achimenes, Gesneriaceae).

RESULTS

The assembled transcriptomes contain between 29,000 and 42,000 genes expressed during development. We combine sequence orthology and coexpression clustering with analyses of protein evolution to identify candidate genes for roles in floral form evolution. Over 25% of transcripts captured were distinctive to Achimenes and overrepresented by genes involved in transcription factor activity. Using a model-based clustering approach we find dynamic, temporal patterns of gene expression among species. Selection tests provide evidence of positive selection in several genes with roles in pigment production, flowering time, and morphology. Combining these approaches to explore genes related to flower color and flower shape, we find distinct patterns that correspond to transitions of floral form among Achimenes species.

CONCLUSIONS

The floral transcriptomes developed from four species of Achimenes provide insight into the mechanisms involved in the evolution of diverse floral form among closely related species with different pollinators. We identified several candidate genes that will serve as an important and useful resource for future research. High conservation of sequence structure, patterns of gene coexpression, and detection of positive selection acting on few genes suggests that large phenotypic differences in floral form may be caused by genetic differences in a small set of genes. Our characterized floral transcriptomes provided here should facilitate further analyses into the genomics of flower development and the mechanisms underlying the evolution of diverse flowers in Achimenes and other Neotropical Gesneriaceae.

Enhanced rhizobial symbiotic capacity in an allopolyploid species of Glycine (Leguminosae)

PREMISE OF THE STUDY

Previous studies have shown that polyploidy can alter biotic interactions, and it has been suggested that these effects may contribute to the increased ability for colonization of new habitats shown by many allopolyploids. Little is known, however, about the effects of allopolyploidy, which combines hybridity and genome doubling, on symbiotic interactions with rhizobial bacteria.

METHODS

We examined interactions of the allopolyploid Glycine dolichocarpa (designated T2) with novel rhizobial partners, such as might occur in a context of colonization, and compared these with the responses of its diploid progenitors, G. tomentella (D3) and G. syndetika (D4). We assessed root hair response, nodule formation, nodule mass, nodule number, and plant biomass.

KEY RESULTS

The allopolyploid (T2) showed a greater root hair deformation response when exposed to rhizobia, compared with either diploid. T2 had a greater probability of forming nodules with NGR234 compared with diploid D4, and greater total nodule mass per nodulated plant compared with diploid D3. T2 also had greater plant biomass responses to nitrogen and when exposed to NGR234.

CONCLUSIONS

The allopolyploid is characterized by transgressive responses to rhizobia for some variables, while also combining certain parental diploid responses such that its capacity for interactions with rhizobia appears to be greater than for either diploid progenitor. This overall enhanced nodulation capacity and the ability to make greater gains from exposure to both rhizobia and additional nitrogen indicate a greater potential of the allopolyploid to benefit from these factors both generally and in a context of colonization.

Genome-wide characterization of the β-1,3-glucanase gene family in Gossypium by comparative analysis

The β-1,3-glucanase gene family is involved in a wide range of plant developmental processes as well as pathogen defense mechanisms. Comprehensive analyses of β-1,3-glucanase genes (GLUs) have not been reported in cotton. Here, we identified 67, 68, 130 and 158 GLUs in four sequenced cotton species, G. raimondii (D5), G. arboreum (A2), G. hirsutum acc. TM-1 (AD1), and G. barbadense acc. 3-79 (AD2), respectively. Cotton GLUs can be classified into the eight subfamilies (A-H), and their protein domain architecture and intron/exon structure are relatively conserved within each subfamily. Sixty-seven GLUs in G. raimondii were anchored onto 13 chromosomes, with 27 genes involved in segmental duplications, and 13 in tandem duplications. Expression patterns showed highly developmental and spatial regulation of GLUs in TM-1. In particular, the expression of individual member of GLUs in subfamily E was limited to roots, leaves, floral organs or fibers. Members of subfamily E also showed more protein evolution and subgenome expression bias compared with members of other subfamilies. We clarified that GLU42 and GLU43 in subfamily E were preferentially expressed in root and leaf tissues and significantly upregulated after Verticillium dahliae inoculation. Silencing of GLU42 and GLU43 significantly increased the susceptibility of cotton to V. dahliae.

The impacts of polyploidy, geographic and ecological isolations on the diversification of Panax (Araliaceae)

BACKGROUND

Panax L. is a medicinally important genus within family Araliaceae, where almost all species are of cultural significance for traditional Chinese medicine. Previous studies suggested two independent origins of the East Asia and North America disjunct distribution of this genus and multiple rounds of whole genome duplications (WGDs) might have occurred during the evolutionary process.

RESULTS

We employed multiple chloroplast and nuclear markers to investigate the evolution and diversification of Panax. Our phylogenetic analyses confirmed previous observations of the independent origins of disjunct distribution and both ancient and recent WGDs have occurred within Panax. The estimations of divergence time implied that the ancient WGD might have occurred before the establishment of Panax. Thereafter, at least two independent recent WGD events have occurred within Panax, one of which has led to the formation of three geographically isolated tetraploid species P. ginseng, P. japonicus and P. quinquefolius. Population genetic analyses showed that the diploid species P. notoginseng harbored significantly lower nucleotide diversity than those of the two tetraploid species P. ginseng and P. quinquefolius and the three species showed distinct nucleotide variation patterns at exon regions.

CONCLUSION

Our findings based on the phylogenetic and population genetic analyses, coupled with the species distribution patterns of Panax, suggested that the two rounds of WGD along with the geographic and ecological isolations might have together contributed to the evolution and diversification of this genus.

SCARN a Novel Class of SCAR Protein That Is Required for Root-Hair Infection during Legume Nodulation

Rhizobial infection of legume root hairs requires a rearrangement of the actin cytoskeleton to enable the establishment of plant-made infection structures called infection threads. In the SCAR/WAVE (Suppressor of cAMP receptor defect/WASP family verpolin homologous protein) actin regulatory complex, the conserved N-terminal domains of SCAR proteins interact with other components of the SCAR/WAVE complex. The conserved C-terminal domains of SCAR proteins bind to and activate the actin-related protein 2/3 (ARP2/3) complex, which can bind to actin filaments catalyzing new actin filament formation by nucleating actin branching. We have identified, SCARN (SCAR-Nodulation), a gene required for root hair infection of Lotus japonicus by Mesorhizobium loti. Although the SCARN protein is related to Arabidopsis thaliana SCAR2 and SCAR4, it belongs to a distinct legume-sub clade. We identified other SCARN-like proteins in legumes and phylogeny analyses suggested that SCARN may have arisen from a gene duplication and acquired specialized functions in root nodule symbiosis. Mutation of SCARN reduced formation of infection-threads and their extension into the root cortex and slightly reduced root-hair length. Surprisingly two of the scarn mutants showed constitutive branching of root hairs in uninoculated plants. However we observed no effect of scarn mutations on trichome development or on the early actin cytoskeletal accumulation that is normally seen in root hair tips shortly after M. loti inoculation, distinguishing them from other symbiosis mutations affecting actin nucleation. The C-terminal domain of SCARN binds to ARPC3 and ectopic expression of the N-terminal SCAR-homology domain (but not the full length protein) inhibited nodulation. In addition, we found that SCARN expression is enhanced by M. loti in epidermal cells and that this is directly regulated by the NODULE INCEPTION (NIN) transcription factor.

Characterization of Lotus japonicus mutants defective for nodule infection by rhizobia led to the identification of a gene we named SCARN. Two of the five alleles caused formation of branched root-hairs in uninoculated seedlings, suggesting SCARN plays a role in the microtubule and actin-regulated polar growth of root hairs. SCARN is one of three L. japonicus proteins containing the conserved N and C terminal domains predicted to be required for rearrangement of the actin cytoskeleton. SCARN expression is induced in response to rhizobial nodulation factors by the NIN (NODULE INCEPTION) transcription factor and appears to be adapted to promoting rhizobial infection, possibly arising from a gene duplication event. SCARN binds to ARPC3, one of the predicted components in the actin-related protein complex involved in the activation of actin nucleation.

Evolutionary signals of symbiotic persistence in the legume-rhizobia mutualism

Understanding the origins and evolutionary trajectories of symbiotic partnerships remains a major challenge. Why are some symbioses lost over evolutionary time whereas others become crucial for survival? Here, we use a quantitative trait reconstruction method to characterize different evolutionary stages in the ancient symbiosis between legumes (Fabaceae) and nitrogen-fixing bacteria, asking how labile is symbiosis across different host clades. We find that more than half of the 1,195 extant nodulating legumes analyzed have a high likelihood (>95%) of being in a state of high symbiotic persistence, meaning that they show a continued capacity to form the symbiosis over evolutionary time, even though the partnership has remained facultative and is not obligate. To explore patterns associated with the likelihood of loss and retention of the N2-fixing symbiosis, we tested for correlations between symbiotic persistence and legume distribution, climate, soil and trait data. We found a strong latitudinal effect and demonstrated that low mean annual temperatures are associated with high symbiotic persistence in legumes. Although no significant correlations between soil variables and symbiotic persistence were found, nitrogen and phosphorus leaf contents were positively correlated with legumes in a state of high symbiotic persistence. This pattern suggests that highly demanding nutrient lifestyles are associated with more stable partnerships, potentially because they "lock" the hosts into symbiotic dependency. Quantitative reconstruction methods are emerging as a powerful comparative tool to study broad patterns of symbiont loss and retention across diverse partnerships.

Identification of domestication-related loci associated with flowering time and seed size in soybean with the RAD-seq genotyping method

Flowering time and seed size are traits related to domestication. However, identification of domestication-related loci/genes of controlling the traits in soybean is rarely reported. In this study, we identified a total of 48 domestication-related loci based on RAD-seq genotyping of a natural population comprising 286 accessions. Among these, four on chromosome 12 and additional two on chromosomes 11 and 15 were associated with flowering time, and four on chromosomes 11 and 16 were associated with seed size. Of the five genes associated with flowering time and the three genes associated with seed size, three genes Glyma11g18720, Glyma11g15480 and Glyma15g35080 were homologous to Arabidopsis genes, additional five genes were found for the first time to be associated with these two traits. Glyma11g18720 and Glyma05g28130 were co-expressed with five genes homologous to flowering time genes in Arabidopsis, and Glyma11g15480 was co-expressed with 24 genes homologous to seed development genes in Arabidopsis. This study indicates that integration of population divergence analysis, genome-wide association study and expression analysis is an efficient approach to identify candidate domestication-related genes.

Tangled up in two: a burst of genome duplications at the end of the Cretaceous and the consequences for plant evolution

Genome sequencing has demonstrated that besides frequent small-scale duplications, large-scale duplication events such as whole genome duplications (WGDs) are found on many branches of the evolutionary tree of life. Especially in the plant lineage, there is evidence for recurrent WGDs, and the ancestor of all angiosperms was in fact most likely a polyploid species. The number of WGDs found in sequenced plant genomes allows us to investigate questions about the roles of WGDs that were hitherto impossible to address. An intriguing observation is that many plant WGDs seem associated with periods of increased environmental stress and/or fluctuations, a trend that is evident for both present-day polyploids and palaeopolyploids formed around the Cretaceous–Palaeogene (K–Pg) extinction at 66 Ma. Here, we revisit the WGDs in plants that mark the K–Pg boundary, and discuss some specific examples of biological innovations and/or diversifications that may be linked to these WGDs. We review evidence for the processes that could have contributed to increased polyploid establishment at the K–Pg boundary, and discuss the implications on subsequent plant evolution in the Cenozoic.

Multiple polyploidy events in the early radiation of nodulating and nonnodulating legumes

Unresolved questions about evolution of the large and diverse legume family include the timing of polyploidy (whole-genome duplication; WGDs) relative to the origin of the major lineages within the Fabaceae and to the origin of symbiotic nitrogen fixation. Previous work has established that a WGD affects most lineages in the Papilionoideae and occurred sometime after the divergence of the papilionoid and mimosoid clades, but the exact timing has been unknown. The history of WGD has also not been established for legume lineages outside the Papilionoideae. We investigated the presence and timing of WGDs in the legumes by querying thousands of phylogenetic trees constructed from transcriptome and genome data from 20 diverse legumes and 17 outgroup species. The timing of duplications in the gene trees indicates that the papilionoid WGD occurred in the common ancestor of all papilionoids. The earliest diverging lineages of the Papilionoideae include both nodulating taxa, such as the genistoids (e.g., lupin), dalbergioids (e.g., peanut), phaseoloids (e.g., beans), and galegoids (=Hologalegina, e.g., clovers), and clades with nonnodulating taxa including Xanthocercis and Cladrastis (evaluated in this study). We also found evidence for several independent WGDs near the base of other major legume lineages, including the Mimosoideae-Cassiinae-Caesalpinieae (MCC), Detarieae, and Cercideae clades. Nodulation is found in the MCC and papilionoid clades, both of which experienced ancestral WGDs. However, there are numerous nonnodulating lineages in both clades, making it unclear whether the phylogenetic distribution of nodulation is due to independent gains or a single origin followed by multiple losses.

The wild side of a major crop: soybean's perennial cousins from Down Under

The accumulation of over 30 years of basic research on the biology, genetic variation, and evolution of the wild perennial relatives of soybean (Glycine max) provides a foundation to improve cultivated soybean. The cultivated soybean and its wild progenitor, G. soja, have a center of origin in eastern Asia and are the only two species in the annual subgenus Soja. Systematic and evolutionary studies of the ca. 30 perennial species of subgenus Glycine, native to Australia, have benefited from the availability of the G. max genomic sequence. The perennial species harbor many traits of interest to soybean breeders, among them resistance to major soybean pathogens such as cyst nematode and leaf rust. New species in the Australian subgenus continue to be described, due to the collection of new material and to insights gleaned through systematic studies of accessions in germplasm collections. Ongoing studies in perennial species focus on genomic regions that contain genes for key traits relevant to soybean breeding. These comparisons also include the homoeologous regions that are the result of polyploidy in the common ancestor of all Glycine species. Subgenus Glycine includes a complex of recently formed allopolyploids that are the focus of studies aimed at elucidating genomic, transcriptomic, physiological, taxonomic, morphological, developmental, and ecological processes related to polyploid evolution. Here we review what has been learned over the past 30 years and outline ongoing work on photosynthesis, nitrogen fixation, and floral biology, much of it drawing on new technologies and resources.

A single evolutionary innovation drives the deep evolution of symbiotic N2-fixation in angiosperms

Symbiotic associations occur in every habitat on earth, but we know very little about their evolutionary histories. Current models of trait evolution cannot adequately reconstruct the deep history of symbiotic innovation, because they assume homogenous evolutionary processes across millions of years. Here we use a recently developed, heterogeneous and quantitative phylogenetic framework to study the origin of the symbiosis between angiosperms and nitrogen-fixing (N₂) bacterial symbionts housed in nodules. We compile the largest database of global nodulating plant species and reconstruct the symbiosis’ evolution. We identify a single, cryptic evolutionary innovation driving symbiotic N₂-fixation evolution, followed by multiple gains and losses of the symbiosis, and the subsequent emergence of ‘stable fixers’ (clades extremely unlikely to lose the symbiosis). Originating over 100 MYA, this innovation suggests deep homology in symbiotic N₂-fixation. Identifying cryptic innovations on the tree of life is key to understanding the evolution of complex traits, including symbiotic partnerships.

Symbiotic associations are widespread, yet their evolutionary histories remain poorly understood. Here, Werner et. al. show a single evolutionary innovation driving symbiotic nitrogen fixation, followed by multiple gains and losses of the symbiosis and the emergence of groups with stable nitrogen fixers.