Genome structure of the legume, Lotus japonicus

Chromosome-level reference genome and resequencing of 322 accessions reveal evolution, genomic imprint and key agronomic traits in adzuki bean

Adzuki bean (Vigna angularis) is an important legume crop cultivated in over 30 countries worldwide. We developed a high-quality chromosome-level reference genome of adzuki bean cultivar Jingnong6 by combining PacBio Sequel long-read sequencing with short-read and Hi-C technologies. The assembled genome covers 97.8% of the adzuki bean genome with a contig N50 of approximately 16 Mb and a total of 32 738 protein-coding genes. We also generated a comprehensive genome variation map of adzuki bean by whole-genome resequencing (WGRS) of 322 diverse adzuki beans accessions including both wild and cultivated. Furthermore, we have conducted comparative genomics and a genome-wide association study (GWAS) on key agricultural traits to investigate the evolution and domestication. GWAS identified several candidate genes, including VaCycA3;1, VaHB15, VaANR1 and VaBm, that exhibited significant associations with domestication traits. Furthermore, we conducted functional analyses on the roles of VaANR1 and VaBm in regulating seed coat colour. We provided evidence for the highest genetic diversity of wild adzuki (Vigna angularis var. nipponensis) in China with the presence of the most original wild adzuki bean, and the occurrence of domestication process facilitating transition from wild to cultigen. The present study elucidates the genetic basis of adzuki bean domestication traits and provides crucial genomic resources to support future breeding efforts in adzuki bean.

Innovations in functional genomics and molecular breeding of pea: exploring advances and opportunities

The garden pea (Pisum sativum L.) is a significant cool-season legume, serving as crucial food sources, animal feed, and industrial raw materials. The advancement of functional genomics over the past two decades has provided substantial theoretical foundations and progress to pea breeding. Notably, the release of the pea reference genome has enhanced our understanding of plant architecture, symbiotic nitrogen fixation (SNF), flowering time, floral organ development, seed development, and stress resistance. However, a considerable gap remains between pea functional genomics and molecular breeding. This review summarizes the current advancements in pea functional genomics and breeding while highlighting the future challenges in pea molecular breeding.

A rare non-canonical splice site in Trema orientalis SYMRK does not affect its dual symbiotic functioning in endomycorrhiza and rhizobium nodulation

BACKGROUND

Nitrogen-fixing nodules occur in ten related taxonomic lineages interspersed with lineages of non-nodulating plant species. Nodules result from an endosymbiosis between plants and diazotrophic bacteria; rhizobia in the case of legumes and Parasponia and Frankia in the case of actinorhizal species. Nodulating plants share a conserved set of symbiosis genes, whereas related non-nodulating sister species show pseudogenization of several key nodulation-specific genes. Signalling and cellular mechanisms critical for nodulation have been co-opted from the more ancient plant-fungal arbuscular endomycorrhizal symbiosis. Studies in legumes and actinorhizal plants uncovered a key component in symbiotic signalling, the LRR-type SYMBIOSIS RECEPTOR KINASE (SYMRK). SYMRK is essential for nodulation and arbuscular endomycorrhizal symbiosis. To our surprise, however, despite its arbuscular endomycorrhizal symbiosis capacities, we observed a seemingly critical mutation in a donor splice site in the SYMRK gene of Trema orientalis, the non-nodulating sister species of Parasponia. This led us to investigate the symbiotic functioning of SYMRK in the Trema-Parasponia lineage and to address the question of to what extent a single nucleotide polymorphism in a donor splice site affects the symbiotic functioning of SYMRK.

RESULTS

We show that SYMRK is essential for nodulation and endomycorrhization in Parasponia andersonii. Subsequently, it is revealed that the 5'-intron donor splice site of SYMRK intron 12 is variable and, in most dicotyledon species, doesn't contain the canonical dinucleotide 'GT' signature but the much less common motif 'GC'. Strikingly, in T. orientalis, this motif is converted into a rare non-canonical 5'-intron donor splice site 'GA'. This SYMRK allele, however, is fully functional and spreads in the T. orientalis population of Malaysian Borneo. A further investigation into the occurrence of the non-canonical GA-AG splice sites confirmed that these are extremely rare.

CONCLUSION

SYMRK functioning is highly conserved in legumes, actinorhizal plants, and Parasponia. The gene possesses a non-common 5'-intron GC donor splice site in intron 12, which is converted into a GA in T. orientalis accessions of Malaysian Borneo. The discovery of this functional GA-AG splice site in SYMRK highlights a gap in our understanding of splice donor sites.

Chromosome-scale assembly and analysis of Melilotus officinalis genome for SSR development and nodulation genes analysis

Melilotus officinalis is an important legume crop with forage and Chinese medicinal value. The unknown genome of M. officinalis restricted the domestication and utilization of the species and its germplasm resource diversity. A chromosome-scale assembly of the M. officinalis genome was assembled and analysed. The 976.27 Mb of genome was divided into eight chromosomes covering 99.16% of the whole genome. A total of 50022 genes were predicted in the genome. M. officinalis and Melilotus albus shared a common ancestor 0.5-5.65 million years ago (MYA). A genome-wide doubling event occurred 68.93 MYA according to the synonymous nucleotide-substitution values. A total of 552102 tandem repeats were predicted, and 46004 SSR primers of TRs with 10 or more base pairs were developed and designed. The elucidation of the M. officinalis genome provides a compelling model system for studying the genetic, evolutionary and biosynthesis of this legume.

The conservation of gene models can support genome annotation

Many genome annotations include false-positive gene models, leading to errors in phylogenetic and comparative studies. Here, we propose a method to support gene model prediction based on evolutionary conservation and use it to identify potentially erroneous annotations. Using this method, we developed a set of 15,345 representative gene models from 12 legume assemblies that can be used to support genome annotations for other legumes.

A Mesorhizobium japonicum quorum sensing circuit that involves three linked genes and an unusual acyl-homoserine lactone signal

Members of the genus Mesorhizobium, which are core components of the rhizosphere and specific symbionts of legume plants, possess genes for acyl-homoserine lactone (AHL) quorum sensing (QS). Here we show Mesorhizobium japonicum MAFF 303099 (formerly M. loti) synthesizes and responds to N-[(2E, 4E)-2,4-dodecadienoyl] homoserine lactone (2E, 4E-C12:2-HSL). We show that the 2E, 4E-C12:2-HSL QS circuit involves one of four luxR-luxI-type genes found in the sequenced genome of MAFF 303099. We refer to this circuit, which appears to be conserved among Mesorhizobium species, as R1-I1. We show that two other Mesorhizobium strains also produce 2E, 4E-C12:2-HSL. The 2E, 4E-C12:2-HSL is unique among known AHLs in its arrangement of two trans double bonds. The R1 response to 2E, 4E-C12:2-HSL is extremely selective in comparison with other LuxR homologs, and the trans double bonds appear critical for R1 signal recognition. Most well-studied LuxI-like proteins use S-adenosylmethionine and an acyl-acyl carrier protein as substrates for synthesis of AHLs. Others that form a subgroup of LuxI-type proteins use acyl-coenzyme A substrates rather than acyl-acyl carrier proteins. I1 clusters with the acyl-coenzyme A-type AHL synthases. We show that a gene linked to the I1 AHL synthase is involved in the production of the QS signal. The discovery of the unique I1 product enforces the view that further study of acyl-coenzyme A-dependent LuxI homologs will expand our knowledge of AHL diversity. The involvement of an additional enzyme in AHL generation leads us to consider this system a three-component QS circuit. IMPORTANCE We report a Mesorhizobium japonicum quorum sensing (QS) system involving a novel acyl-homoserine lactone (AHL) signal. This system is known to be involved in root nodule symbiosis with host plants. The chemistry of the newly described QS signal indicated that there may be a dedicated cellular enzyme involved in its synthesis in addition to the types known for production of other AHLs. Indeed, we report that an additional gene is required for synthesis of the unique signal, and we propose that this is a three-component QS circuit as opposed to the canonical two-component AHL QS circuits. The signaling system is exquisitely selective. The selectivity may be important when this species resides in the complex microbial communities around host plants and may make this system useful in various synthetic biology applications of QS circuits.

Resequencing of global Lotus corniculatus accessions reveals population distribution and genetic loci, associated with cyanogenic glycosides accumulation and growth traits

BACKGROUND

Lotus corniculatus is a widely distributed perennial legume whose great adaptability to different environments and resistance to barrenness make it an excellent forage and ecological restoration plant. However, its molecular genetics and genomic relationships among populations are yet to be uncovered.

RESULT

Here we report on a genomic variation map from worldwide 272 L. corniculatus accessions by genome resequencing. Our analysis suggests that L. corniculatus accessions have high genetic diversity and could be further divided into three subgroups, with the genetic diversity centers were located in Transcaucasia. Several candidate genes and SNP site associated with CNglcs content and growth traits were identified by genome-wide associated study (GWAS). A non-synonymous in LjMTR was responsible for the decreased expression of CNglcs synthesis genes and LjZCD was verified to positively regulate CNglcs synthesis gene CYP79D3. The LjZCB and an SNP in LjZCA promoter were confirmed to be involved in plant growth.

CONCLUSION

This study provided a large number of genomic resources and described genetic relationship and population structure among different accessions. Moreover, we attempt to provide insights into the molecular studies and breeding of CNglcs and growth traits in L. corniculatus.

Development of a knowledge graph framework to ease and empower translational approaches in plant research: a use-case on grain legumes

While the continuing decline in genotyping and sequencing costs has largely benefited plant research, some key species for meeting the challenges of agriculture remain mostly understudied. As a result, heterogeneous datasets for different traits are available for a significant number of these species. As gene structures and functions are to some extent conserved through evolution, comparative genomics can be used to transfer available knowledge from one species to another. However, such a translational research approach is complex due to the multiplicity of data sources and the non-harmonized description of the data. Here, we provide two pipelines, referred to as structural and functional pipelines, to create a framework for a NoSQL graph-database (Neo4j) to integrate and query heterogeneous data from multiple species. We call this framework Orthology-driven knowledge base framework for translational research (Ortho_KB). The structural pipeline builds bridges across species based on orthology. The functional pipeline integrates biological information, including QTL, and RNA-sequencing datasets, and uses the backbone from the structural pipeline to connect orthologs in the database. Queries can be written using the Neo4j Cypher language and can, for instance, lead to identify genes controlling a common trait across species. To explore the possibilities offered by such a framework, we populated Ortho_KB to obtain OrthoLegKB, an instance dedicated to legumes. The proposed model was evaluated by studying the conservation of a flowering-promoting gene. Through a series of queries, we have demonstrated that our knowledge graph base provides an intuitive and powerful platform to support research and development programmes.

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.

DNA methylation signatures of duplicate gene evolution in angiosperms

Gene duplication is a source of evolutionary novelty. DNA methylation may play a role in the evolution of duplicate genes (paralogs) through its association with gene expression. While this relationship has been examined to varying extents in a few individual species, the generalizability of these results at either a broad phylogenetic scale with species of differing duplication histories or across a population remains unknown. We applied a comparative epigenomics approach to 43 angiosperm species across the phylogeny and a population of 928 Arabidopsis (Arabidopsis thaliana) accessions, examining the association of DNA methylation with paralog evolution. Genic DNA methylation was differentially associated with duplication type, the age of duplication, sequence evolution, and gene expression. Whole genome duplicates were typically enriched for CG-only gene-body methylated or unmethylated genes, while single-gene duplications were typically enriched for non-CG methylated or unmethylated genes. Non-CG methylation, in particular, was characteristic of more recent single-gene duplicates. Core angiosperm gene families differentiated into those which preferentially retain paralogs and 'duplication-resistant' families, which convergently reverted to singletons following duplication. Duplication-resistant families that still have paralogous copies were, uncharacteristically for core angiosperm genes, enriched for non-CG methylation. Non-CG methylated paralogs had higher rates of sequence evolution, higher frequency of presence-absence variation, and more limited expression. This suggests that silencing by non-CG methylation may be important to maintaining dosage following duplication and be a precursor to fractionation. Our results indicate that genic methylation marks differing evolutionary trajectories and fates between paralogous genes and have a role in maintaining dosage following duplication.

Legume-wide comparative analysis of pod shatter locus PDH1 reveals phaseoloid specificity, high cowpea expression, and stress responsive genomic context

Pod dehiscence is a major source of yield loss in legumes, which is exacerbated by aridity. Disruptive mutations in "Pod indehiscent 1" (PDH1), a pod sclerenchyma-specific lignin biosynthesis gene, has been linked to significant reductions in dehiscence in several legume species. We compared syntenic PDH1 regions across 12 legumes and two outgroups to uncover key historical evolutionary trends at this important locus. Our results clarified the extent to which PDH1 orthologs are present in legumes, showing the typical genomic context surrounding PDH1 has only arisen relatively recently in certain phaseoloid species (Vigna, Phaseolus, Glycine). The notable absence of PDH1 in Cajanus cajan may be a major contributor to its indehiscent phenotype compared with other phaseoloids. In addition, we identified a novel PDH1 ortholog in Vigna angularis and detected remarkable increases in PDH1 transcript abundance during Vigna unguiculata pod development. Investigation of the shared genomic context of PDH1 revealed it lies in a hotspot of transcription factors and signaling gene families that respond to abscisic acid and drought stress, which we hypothesize may be an additional factor influencing expression of PDH1 under specific environmental conditions. Our findings provide key insights into the evolutionary history of PDH1 and lay the foundation for optimizing the pod dehiscence role of PDH1 in major and understudied legume species.

Medicago truncatula resources to study legume biology and symbiotic nitrogen fixation

Medicago truncatula is a chosen model for legumes towards deciphering fundamental legume biology, especially symbiotic nitrogen fixation. Current genomic resources for M. truncatula include a completed whole genome sequence information for R108 and Jemalong A17 accessions along with the sparse draft genome sequences for other 226 M. truncatula accessions. These genomic resources are complemented by the availability of mutant resources such as retrotransposon (Tnt1) insertion mutants in R108 and fast neutron bombardment (FNB) mutants in A17. In addition, several M. truncatula databases such as small secreted peptides (SSPs) database, transporter protein database, gene expression atlas, proteomic atlas, and metabolite atlas are available to the research community. This review describes these resources and provide information regarding how to access these resources.

An improved reference genome for Trifolium subterraneum L. provides insight into molecular diversity and intra-specific phylogeny

Subterranean clover (Trifolium subterraneum L., Ts) is a geocarpic, self-fertile annual forage legume with a compact diploid genome (n = x = 8, 544 Mb/1C). Its resilience and climate adaptivity have made it an economically important species in Mediterranean and temperate zones. Using the cultivar Daliak, we generated higher resolution sequence data, created a new genome assembly TSUd_3.0, and conducted molecular diversity analysis for copy number variant (CNV) and single-nucleotide polymorphism (SNP) among 36 cultivars. TSUd_3.0 substantively improves prior genome assemblies with new Hi-C and long-read sequence data, covering 531 Mb, containing 41,979 annotated genes and generating a 94.4% BUSCO score. Comparative genomic analysis among select members of the tribe Trifolieae indicated TSUd 3.0 corrects six assembly-error inversion/duplications and confirmed phylogenetic relationships. Its synteny with T. pratense, T. repens, Medicago truncatula and Lotus japonicus genomes were assessed, with the more distantly related T. repens and M. truncatula showing higher levels of co-linearity with Ts than between Ts and its close relative T. pratense. Resequencing of 36 cultivars discovered 7,789,537 SNPs subsequently used for genomic diversity assessment and sequence-based clustering. Heterozygosity estimates ranged from 1% to 21% within the 36 cultivars and may be influenced by admixture. Phylogenetic analysis supported subspecific genetic structure, although it indicates four or five groups, rather than the three recognized subspecies. Furthermore, there were incidences where cultivars characterized as belonging to a particular subspecies clustered with another subspecies when using genomic data. These outcomes suggest that further investigation of Ts sub-specific classification using molecular and morpho-physiological data is needed to clarify these relationships. This upgraded reference genome, complemented with comprehensive sequence diversity analysis of 36 cultivars, provides a platform for future gene functional analysis of key traits, and genome-based breeding strategies for climate adaptation and agronomic performance. Pangenome analysis, more in-depth intra-specific phylogenomic analysis using the Ts core collection, and functional genetic and genomic studies are needed to further augment knowledge of Trifolium genomes.

Genetic Diversity, Population Structure and Subset Development in a Sesbania sesban Collection

Sesbania sesban (L.) Merr. is a multipurpose legume tree grown primarily for fodder and forage in the tropical and subtropical world. In this study, the Sesbania sesban collection maintained in the International Livestock Research Institute (ILRI) forage Genebank was studied using genome-wide markers generated on the DArTseq platform. Genotyping produced 84,673 and 60,626 SNP and SilicoDArT markers with a mean polymorphic information content of 0.153 and 0.123, respectively. From the generated markers, 7587 and 15,031 highly informative SNP and SilicoDArT markers, respectively, were filtered and used for genetic diversity analysis and subset development. Analysis of molecular variance (AMOVA) revealed higher variability 'within' (52.73% for SNP markers and 67.36% for SilicoDArT markers) than 'between' accessions. Hierarchical cluster analysis showed the presence of four main clusters in the collection. Mantel correlation analysis showed a lack of relationship between genetic variation of the germplasm and their geographical origin. A representative subset of 34 accessions containing germplasm from diverse origins and agro-ecologies was developed using SNP markers. The genetic diversity information generated in this study could be used for marker-assisted screening for stress tolerance, gap analysis and identification and acquisition of new distinct genotype(s) to broaden the genetic basis of the collection for future improvement programs to develop high-yielding, stress-tolerant varieties for enhancing food and environmental security in crop-livestock-based production systems.

Characterization of zinc uptake and translocation visualized with positron-emitting 65Zn tracer and analysis of transport-related gene expression in two Lotus japonicus accessions

BACKGROUND AND AIMS

Zinc (Zn) is an essential element for humans and plants. However, Zn deficiency is widespread and 25 % of the world's population is at risk of Zn deficiency. To overcome the deficiency of Zn intake, crops with high Zn content are required. However, most crop-producing areas have Zn-deficient soils, therefore crops with excellent Zn uptake/transport characteristics (i.e. high Zn efficiency) are needed. Our objective was to identify the crucial factors responsible for high Zn efficiency in the legume Lotus japonicus.

METHODS

We evaluated Zn efficiency by static and real-time visualization of radioactive Zn (65Zn) uptake/transport in two L. japonicus accessions, MG-20 and B-129, that differ in Zn efficiency. The combination of visualization methods verified the dynamics of Zn accumulation and transport within the plant. We compared gene expression under a normal Zn concentration (control) and Zn deficiency to evaluate genetic factors that may determine the differential Zn efficiency of the accessions.

KEY RESULTS

The accession B-129 accumulated almost twice the amount of Zn as MG-20. In the static 65Zn images, 65Zn accumulated in meristematic tissues, such as root tips and the shoot apex, in both accessions. The positron-emitting tracer imaging system (PETIS), which follows the transport process in real time, revealed that 65Zn transport to the shoot was more rapid in B-129 than in MG-20. Many genes associated with Zn uptake and transport were more highly expressed in B-129 than in MG-20 under the control condition. These gene expression patterns under Zn deficiency differed from those under the control Zn condition.

CONCLUSIONS

PETIS confirmed that the real-time transport of 65Zn to the shoot was faster in B-129 than in MG-20. The high Zn efficiency of B-129 may be due to the elevated expression of a suite of Zn uptake- and transport-related genes.

Mendel: From genes to genome

Two hundred years after the birth of Gregor Mendel, it is an appropriate time to reflect on recent developments in the discipline of genetics, particularly advances relating to the prescient friar's model species, the garden pea (Pisum sativum L.). Mendel's study of seven characteristics established the laws of segregation and independent assortment. The genes underlying four of Mendel's loci (A, LE, I, and R) have been characterized at the molecular level for over a decade. However, the three remaining genes, influencing pod color (GP), pod form (V/P), and the position of flowers (FA/FAS), have remained elusive for a variety of reasons, including a lack of detail regarding the loci with which Mendel worked. Here, we discuss potential candidate genes for these characteristics, in light of recent advances in the genetic resources for pea. These advances, including the pea genome sequence and reverse-genetics techniques, have revitalized pea as an excellent model species for physiological-genetic studies. We also discuss the issues that have been raised with Mendel's results, such as the recent controversy regarding the discrete nature of the characters that Mendel chose and the perceived overly-good fit of his segregations to his hypotheses. We also consider the relevance of these controversies to his lasting contribution. Finally, we discuss the use of Mendel's classical results to teach and enthuse future generations of geneticists, not only regarding the core principles of the discipline, but also its history and the role of hypothesis testing.

A pan-genome and chromosome-length reference genome of narrow-leafed lupin (Lupinus angustifolius) reveals genomic diversity and insights into key industry and biological traits

Narrow-leafed lupin (NLL; Lupinus angustifolius) is a key rotational crop for sustainable farming systems, whose grain is high in protein content. It is a gluten-free, non-genetically modified, alternative protein source to soybean (Glycine max) and as such has gained interest as a human food ingredient. Here, we present a chromosome-length reference genome for the species and a pan-genome assembly comprising 55 NLL lines, including Australian and European cultivars, breeding lines and wild accessions. We present the core and variable genes for the species and report on the absence of essential mycorrhizal associated genes. The genome and pan-genomes of NLL and its close relative white lupin (Lupinus albus) are compared. Furthermore, we provide additional evidence supporting LaRAP2-7 as the key alkaloid regulatory gene for NLL and demonstrate the NLL genome is underrepresented in classical NLR disease resistance genes compared to other sequenced legume species. The NLL genomic resources generated here coupled with previously generated RNA sequencing datasets provide new opportunities to fast-track lupin crop improvement.

Genome sequencing and analysis uncover the regulatory elements involved in the development and oil biosynthesis of Pongamia pinnata (L.) - A potential biodiesel feedstock

Due to rapid industrialization, the consumption of petro-products has increased, while fossil fuel resources have been gradually depleted. There has been a resurgence of interest in plant-derived biofuels as a sustainable alternative to fossil fuels for the purpose of reducing greenhouse gas emissions. Pongamia pinnata L., which is also known as Millettia pinnata is an oil-yielding, leguminous tree with a large and complex genome. Despite its multiple industrial applications, this orphan tree species has inconsistent yields and a limited understanding of its functional genomics. We assessed physiological and morphological characteristics of five high-yielding pongamia accessions and deduced important yield descriptors. Furthermore, we sequenced the genome of this potential biofuel feedstock using Illumina HiSeq, NextSeq, and MiSeq platforms to generate paired-end reads. Around 173 million processed reads amounting to 65.2 Gb were assembled into a 685 Mb genome, with a gap rate of 0.02%. The sequenced scaffolds were used to identify 30,000 gene models, 406,385 Simple-Sequence-Repeat (SSR) markers, and 43.6% of repetitive sequences. We further analyzed the structural information of genes belonging to certain key metabolic pathways, including lipid metabolism, photosynthesis, circadian rhythms, plant-pathogen interactions, and karanjin biosynthesis, all of which are commercially significant for pongamia. A total of 2,219 scaffolds corresponding to 29 transcription factor families provided valuable information about gene regulation in pongamia. Similarity studies and phylogenetic analysis revealed a monophyletic group of Fabaceae members wherein pongamia out-grouped from Glycine max and Cajanus cajan, revealing its unique ability to synthesize oil for biodiesel. This study is the first step toward completing the genome sequence of this imminent biofuel tree species. Further attempts at re-sequencing with different read chemistry will certainly improve the genetic resources at the chromosome level and accelerate the molecular breeding programs.

Assembly of high-quality genomes of the locoweed Oxytropis ochrocephala and its endophyte Alternaria oxytropis provides new evidence for their symbiotic relationship and swainsonine biosynthesis

Locoweeds are perennial forbs poisonous to livestock and cause extreme losses to animal husbandry. Locoweed toxicity is attributed to the symbiotic endophytes in Alternaria sect. Undifilum, which produce a mycotoxin swainsonine (SW). We performed a de novo whole genome sequencing of the most common locoweed in China, Oxytropis ochrocephala (2n = 16), and assembled a high-quality, chromosome-level reference genome. Its genome size is 958.83 Mb with 930.94 Mb (97.09 %) anchored and oriented onto 8 chromosomes, and 31,700 protein-coding genes were annotated. Phylogenetic and collinearity analysis showed it is closely related to Medicago truncatula with a pair of large interchromosomal rearrangements, and both species underwent a whole-genome duplication event. We also derived the genome of A. oxytropis at 74.48 Mb with a contig N50 of 8.87 Mb and 10,657 protein-coding genes, and refined the genes of SW biosynthesis. Multiple Alternaria species containing the swnK gene were grouped into a single clade, but in other genera, swnK's homologues are diverse. Resequencing of 41 A. oxytropis strains revealed one SNP in the SWN cluster causing changes in SW concentration. Comparing the transcriptomes of symbiotic and non-symbiotic interactions identified differentially expressed genes (DEG) linked to defense and secondary metabolism in the host. Within the endophyte DEGs were linked to cell wall degradation, fatty acids and nitrogen metabolism. Symbiosis induced the up-regulation of most of the SW biosynthetic genes. These two genomes and relevant sequencing data should provide valuable genetic resources for the study of the evolution, interaction, and SW biosynthesis in the symbiont.

The Genome of the Mimosoid Legume Prosopis cineraria, a Desert Tree

The mimosoid legumes are a clade of ~40 genera in the Caesalpinioideae subfamily of the Fabaceae that grow in tropical and subtropical regions. Unlike the better studied Papilionoideae, there are few genomic resources within this legume group. The tree Prosopis cineraria is native to the Near East and Indian subcontinent, where it thrives in very hot desert environments. To develop a tool to better understand desert plant adaptation mechanisms, we sequenced the P. cineraria genome to near-chromosomal assembly, with a total sequence length of ~691 Mb. We predicted 77,579 gene models (76,554 CDS, 361 rRNAs and 664 tRNAs) from the assembled genome, among them 55,325 (~72%) protein-coding genes that were functionally annotated. This genome was found to consist of over 58% repeat sequences, primarily long terminal repeats (LTR-)-retrotransposons. We find an expansion of terpenoid metabolism genes in P. cineraria and its relative Prosopis alba, but not in other legumes. We also observed an amplification of NBS-LRR disease-resistance genes correlated with LTR-associated retrotransposition, and identified 410 retrogenes with an active burst of chimeric retrogene creation that approximately occurred at the same time of divergence of P. cineraria from a common lineage with P. alba~23 Mya. These retrogenes include many biotic defense responses and abiotic stress stimulus responses, as well as the early Nodulin 93 gene. Nodulin 93 gene amplification is consistent with an adaptive response of the species to the low nitrogen in arid desert soil. Consistent with these results, our differentially expressed genes show a tissue specific expression of isoprenoid pathways in shoots, but not in roots, as well as important genes involved in abiotic salt stress in both tissues. Overall, the genome sequence of P. cineraria enriches our understanding of the genomic mechanisms of its disease resistance and abiotic stress tolerance. Thus, it is a very important step in crop and legume improvement.

Gene-Editing Technologies and Applications in Legumes: Progress, Evolution, and Future Prospects

Legumes are rich in protein and phytochemicals and have provided a healthy diet for human beings for thousands of years. In recognition of the important role they play in human nutrition and agricultural production, the researchers have made great efforts to gain new genetic traits in legumes such as yield, stress tolerance, and nutritional quality. In recent years, the significant increase in genomic resources for legume plants has prepared the groundwork for applying cutting-edge breeding technologies, such as transgenic technologies, genome editing, and genomic selection for crop improvement. In addition to the different genome editing technologies including the CRISPR/Cas9-based genome editing system, this review article discusses the recent advances in plant-specific gene-editing methods, as well as problems and potential benefits associated with the improvement of legume crops with important agronomic properties. The genome editing technologies have been effectively used in different legume plants including model legumes like alfalfa and lotus, as well as crops like soybean, cowpea, and chickpea. We also discussed gene-editing methods used in legumes and the improvements of agronomic traits in model and recalcitrant legumes. Despite the immense opportunities genome editing can offer to the breeding of legumes, governmental regulatory restrictions present a major concern. In this context, the comparison of the regulatory framework of genome editing strategies in the European Union and the United States of America was also discussed. Gene-editing technologies have opened up new possibilities for the improvement of significant agronomic traits in legume breeding.

Low Temperature Stress Tolerance: An Insight Into the Omics Approaches for Legume Crops

The change in climatic conditions is the major cause for decline in crop production worldwide. Decreasing crop productivity will further lead to increase in global hunger rate. Climate change results in environmental stress which has negative impact on plant-like deficiencies in growth, crop yield, permanent damage, or death if the plant remains in the stress conditions for prolonged period. Cold stress is one of the main abiotic stresses which have already affected the global crop production. Cold stress adversely affects the plants leading to necrosis, chlorosis, and growth retardation. Various physiological, biochemical, and molecular responses under cold stress have revealed that the cold resistance is more complex than perceived which involves multiple pathways. Like other crops, legumes are also affected by cold stress and therefore, an effective technique to mitigate cold-mediated damage is critical for long-term legume production. Earlier, crop improvement for any stress was challenging for scientific community as conventional breeding approaches like inter-specific or inter-generic hybridization had limited success in crop improvement. The availability of genome sequence, transcriptome, and proteome data provides in-depth sight into different complex mechanisms under cold stress. Identification of QTLs, genes, and proteins responsible for cold stress tolerance will help in improving or developing stress-tolerant legume crop. Cold stress can alter gene expression which further leads to increases in stress protecting metabolites to cope up the plant against the temperature fluctuations. Moreover, genetic engineering can help in development of new cold stress-tolerant varieties of legume crop. This paper provides a general insight into the "omics" approaches for cold stress in legume crops.

Spatiotemporal changes in gibberellin content are required for soybean nodulation

The plant hormone gibberellin (GA) is required at different stages of legume nodule development, with its spatiotemporal distribution tightly regulated. Transcriptomic and bioinformatic analyses established that several key GA biosynthesis and catabolism enzyme encoding genes are critical to soybean (Glycine max) nodule formation. We examined the expression of several GA oxidase genes and used a Förster resonance energy transfer-based GA biosensor to determine the bioactive GA content of roots inoculated with DsRed-labelled Bradyrhizobium diazoefficiens. We manipulated the level of GA by genetically disrupting the expression of GA oxidase genes. Moreover, exogenous treatment of soybean roots with GA3 induced the expression of key nodulation genes and altered infection thread and nodule phenotypes. GmGA20ox1a, GmGA3ox1a, and GmGA2ox1a are upregulated in soybean roots inoculated with compatible B. diazoefficiens. GmGA20ox1a expression is predominately localized to the transient meristem of soybean nodules and coincides with the spatiotemporal distribution of bioactive GA occurring throughout nodule organogenesis. GmGA2ox1a exhibits a nodule vasculature-specific expression pattern, whereas GmGA3ox1a can be detected throughout the nodule and root. Disruptions in the level of GA resulted in aberrant rhizobia infection and reduced nodule numbers. Collectively, our results establish a central role for GAs in root hair infection by symbiotic rhizobia and in nodule organogenesis.

Genome and systems biology of Melilotus albus provides insights into coumarins biosynthesis

Melilotus species are used as green manure and rotation crops worldwide and contain abundant pharmacologically active coumarins. However, there is a paucity of information on its genome and coumarin production and function. Here, we reported a chromosome-scale assembly of Melilotus albus genome with 1.04 Gb in eight chromosomes, containing 71.42% repetitive elements. Long terminal repeat retrotransposon bursts coincided with declining of population sizes during the Quaternary glaciation. Resequencing of 94 accessions enabled insights into genetic diversity, population structure, and introgression. Melilotus officinalis had relatively larger genetic diversity than that of M. albus. The introgression existed between M. officinalis group and M. albus group, and gene flows was from M. albus to M. officinalis. Selection sweep analysis identified candidate genes associated with flower colour and coumarin biosynthesis. Combining genomics, BSA, transcriptomics, metabolomics, and biochemistry, we identified a β-glucosidase (BGLU) gene cluster contributing to coumarin biosynthesis. MaBGLU1 function was verified by overexpression in M. albus, heterologous expression in Escherichia coli, and substrate feeding, revealing its role in scopoletin (coumarin derivative) production and showing that nonsynonymous variation drives BGLU enzyme activity divergence in Melilotus. Our work will accelerate the understanding of biologically active coumarins and their biosynthetic pathways, and contribute to genomics-enabled Melilotus breeding.

The Diesel Tree Sindora glabra Genome Provides Insights Into the Evolution of Oleoresin Biosynthesis

Sindora glabra is an economically important tree that produces abundant oleoresin in the trunk. Here, we present a high-quality chromosome-scale assembly of S. glabra genome by combining Illumina HiSeq, Pacific Biosciences sequencing, and Hi-C technologies. The size of S. glabra genome was 1.11 Gb, with a contig N50 of 1.27 Mb and 31,944 predicted genes. This is the first sequenced genome of the subfamily Caesalpinioideae. As a sister taxon to Papilionoideae, S. glabra underwent an ancient genome triplication shared by core eudicots and further whole-genome duplication shared by early-legume in the last 73.3 million years. S. glabra harbors specific genes and expanded genes largely involved in stress responses and biosynthesis of secondary metabolites. Moreover, 59 terpene backbone biosynthesis genes and 64 terpene synthase genes were identified, which together with co-expressed transcription factors could contribute to the diversity and specificity of terpene compounds and high terpene content in S. glabra stem. In addition, 63 disease resistance NBS-LRR genes were found to be unique in S. glabra genome and their expression levels were correlated with the accumulation of terpene profiles, suggesting potential defense function of terpenes in S. glabra. These together provide new resources for understanding genome evolution and oleoresin production.

Phenolic Acids Induce Nod Factor Production in <i>Lotus japonicus</i>–<i>Mesorhizobium</i> Symbiosis

In legume–rhizobia symbiosis, partner recognition and the initiation of symbiosis processes require the mutual exchange of chemical signals. Chemicals, generally (iso)flavonoids, in the root exudates of the host plant induce the expression of nod genes in rhizobia, and, thus, are called nod gene inducers. The expression of nod genes leads to the production of lipochitooligosaccharides (LCOs) called Nod factors. Natural nod gene inducer(s) in Lotus japonicus–Mesorhizobium symbiosis remain unknown. Therefore, we developed an LCO detection method based on ultra-high-performance liquid chromatography–tandem-quadrupole mass spectrometry (UPLC-TQMS) to identify these inducers and used it herein to screen 40 phenolic compounds and aldonic acids for their ability to induce LCOs in Mesorhizobium japonicum MAFF303099. We identified five phenolic acids with LCO-inducing activities, including p-coumaric, caffeic, and ferulic acids. The induced LCOs caused root hair deformation, and nodule numbers in L. japonicus inoculated with M. japonicum were increased by these phenolic acids. The three phenolic acids listed above induced the expression of the nodA, nodB, and ttsI genes in a strain harboring a multicopy plasmid encoding NodD1, but not that encoding NodD2. The presence of p-coumaric and ferulic acids in the root exudates of L. japonicus was confirmed by UPLC-TQMS, and the induction of ttsI::lacZ in the strain harboring the nodD1 plasmid was detected in the rhizosphere of L. japonicus. Based on these results, we propose that phenolic acids are a novel type of nod gene inducer in L. japonicus–Mesorhizobium symbiosis.

Pathogen-driven coevolution across the CBP60 plant immune regulator subfamilies confers resilience on the regulator module

Components of the plant immune signaling network need mechanisms that confer resilience against fast-evolving pathogen effectors that target them. Among eight Arabidopsis CaM-Binding Protein (CBP) 60 family members, AtCBP60g and AtSARD1 are partially functionally redundant, major positive immune regulators, and AtCBP60a is a negative immune regulator. We investigated possible resilience-conferring evolutionary mechanisms among the CBP60a, CBP60g and SARD1 immune regulatory subfamilies. Phylogenetic analysis was used to investigate the times of CBP60 subfamily neofunctionalization. Then, using the pairwise distance rank based on the newly developed analytical platform Protein Evolution Analysis in a Euclidean Space (PEAES), hypotheses of specific coevolutionary mechanisms that could confer resilience on the regulator module were tested. The immune regulator subfamilies diversified around the time of angiosperm divergence and have been evolving very quickly. We detected significant coevolutionary interactions across the immune regulator subfamilies in all of 12 diverse core eudicot species lineages tested. The coevolutionary interactions were consistent with the hypothesized coevolution mechanisms. Despite their unusually fast evolution, members across the CBP60 immune regulator subfamilies have influenced the evolution of each other long after their diversification in a way that could confer resilience on the immune regulator module against fast-evolving pathogen effectors.

A multispecies amplicon sequencing approach for genetic diversity assessments in grassland plant species

Grasslands are widespread and economically relevant ecosystems at the basis of sustainable roughage production. Plant genetic diversity (PGD; i.e., within‐species diversity) is related to many beneficial effects on the ecosystem functioning of grasslands. The monitoring of PGD in temperate grasslands is complicated by the multiplicity of species present and by a shortage of methods for large‐scale assessments. However, the continuous advancement of high‐throughput DNA sequencing approaches has improved the prospects of broad, multispecies PGD monitoring. Among them, amplicon sequencing stands out as a robust and cost‐effective method. Here, we report a set of 12 multispecies primer pairs that can be used for high‐throughput PGD assessments in multiple grassland plant species. The target loci were selected and tested in two phases: a “discovery phase” based on a sequence capture assay (611 nuclear loci assessed in 16 grassland plant species), which resulted in the selection of 11 loci; and a “validation phase”, in which the selected loci were targeted and sequenced using multispecies primers in test populations of Dactylis glomerata L., Lolium perenne L., Festuca pratensis Huds., Trifolium pratense L. and T. repens L. The multispecies amplicons had nucleotide diversities per species from 5.19 × 10⁻³ to 1.29 × 10⁻², which is in the range of flowering‐related genes but slightly lower than pathogen resistance genes. We conclude that the methodology, the DNA sequence resources, and the primer pairs reported in this study provide the basis for large‐scale, multispecies PGD monitoring in grassland plants.

Stress-regulated elements in Lotus spp., as a possible starting point to understand signalling networks and stress adaptation in legumes

Although legumes are of primary economic importance for human and livestock consumption, the information regarding signalling networks during plant stress response in this group is very scarce. Lotus japonicus is a major experimental model within the Leguminosae family, whereas L. corniculatus and L. tenuis are frequent components of natural and agricultural ecosystems worldwide. These species display differences in their perception and response to diverse stresses, even at the genotype level, whereby they have been used in many studies aimed at achieving a better understanding of the plant stress-response mechanisms. However, we are far from the identification of key components of their stress-response signalling network, a previous step for implementing transgenic and editing tools to develop legume stress-resilient genotypes, with higher crop yield and quality. In this review we scope a body of literature, highlighting what is currently known on the stress-regulated signalling elements so far reported in Lotus spp. Our work includes a comprehensive review of transcription factors chaperones, redox signals and proteins of unknown function. In addition, we revised strigolactones and genes regulating phytochelatins and hormone metabolism, due to their involvement as intermediates in several physiological signalling networks. This work was intended for a broad readership in the fields of physiology, metabolism, plant nutrition, genetics and signal transduction. Our results suggest that Lotus species provide a valuable information platform for the study of specific protein-protein (PPI) interactions, as a starting point to unravel signalling networks underlying plant acclimatation to bacterial and abiotic stressors in legumes. Furthermore, some Lotus species may be a source of genes whose regulation improves stress tolerance and growth when introduced ectopically in other plant species.

Deciphering Evolutionary Dynamics of WRKY I Genes in Rosaceae Species

WRKY transcription factors participate in various regulation processes at different developmental stages in higher plants. Here, 98 WRKY I genes were identified in seven Rosaceae species. The WRKY I genes are highly enriched in some subgroups and are selectively expanded in Chinese pear [Pyrus bretschneideri (P. bretschneideri)] and apple [Malus domestica (M. domestica)]. By searching for intra-species gene microsynteny, we found the majority of chromosomal segments for WRKY I-containing segments in both P. bretschneideri and M. domestica genomes, while paired segments were hardly identified in the other five genomes. Furthermore, we analyzed the environmental selection pressure of duplicated WRKY I gene pairs, which indicated that the strong purifying selection for WRKY domains may contribute to the stability of its structure and function. The expression patterns of duplication PbWRKY genes revealed that functional redundancy for some of these genes was derived from common ancestry and neo-functionalization or sub-functionalization for some of them. This study traces the evolution of WRKY I genes in Rosaceae genomes and lays the foundation for functional studies of these genes in the future. Our results also show that the rates of gene loss and gain in different Rosaceae genomes are far from equilibrium.

Chromosome-length genome assemblies of six legume species provide insights into genome organization, evolution, and agronomic traits for crop improvement

INTRODUCTION

Legume crops are an important source of protein and oil for human health and in fixing atmospheric N₂ for soil enrichment. With an objective to accelerate much-needed genetic analyses and breeding applications, draft genome assemblies were generated in several legume crops; many of them are not high quality because they are mainly based on short reads. However, the superior quality of genome assembly is crucial for a detailed understanding of genomic architecture, genome evolution, and crop improvement.

OBJECTIVES

Present study was undertaken with an objective of developing improved chromosome-length genome assemblies in six different legumes followed by their systematic investigation to unravel different aspects of genome organization and legume evolution.

METHODS

We employed in situ Hi-C data to improve the existing draft genomes and performed different evolutionary and comparative analyses using improved genome assemblies.

RESULTS

We have developed chromosome-length genome assemblies in chickpea, pigeonpea, soybean, subterranean clover, and two wild progenitor species of cultivated groundnut (A. duranensis and A. ipaensis). A comprehensive comparative analysis of these genome assemblies offered improved insights into various evolutionary events that shaped the present-day legume species. We highlighted the expansion of gene families contributing to unique traits such as nodulation in legumes, gravitropism in groundnut, and oil biosynthesis in oilseed legume crops such as groundnut and soybean. As examples, we have demonstrated the utility of improved genome assemblies for enhancing the resolution of "QTL-hotspot" identification for drought tolerance in chickpea and marker-trait associations for agronomic traits in pigeonpea through genome-wide association study. Genomic resources developed in this study are publicly available through an online repository, 'Legumepedia'.

CONCLUSION

This study reports chromosome-length genome assemblies of six legume species and demonstrates the utility of these assemblies in crop improvement. The genomic resources developed here will have significant role in accelerating genetic improvement applications of legume crops.

Genome and transcriptome-based characterization of high energy carbon-ion beam irradiation induced delayed flower senescence mutant in Lotus japonicus

BACKGROUND

Flower longevity is closely related to pollen dispersal and reproductive success in all plants, as well as the commercial value of ornamental plants. Mutants that display variation in flower longevity are useful tools for understanding the mechanisms underlying this trait. Heavy-ion beam irradiation has great potential to improve flower shapes and colors; however, few studies are available on the mutation of flower senescence in leguminous plants.

RESULTS

A mutant (C416) exhibiting blossom duration eight times longer than that of the wild type (WT) was isolated in Lotus japonicus derived from carbon ion beam irradiation. Genetic assays supported that the delayed flower senescence of C416 was a dominant trait controlled by a single gene, which was located between 4,616,611 Mb and 5,331,876 Mb on chromosome III. By using a sorting strategy of multi-sample parallel genome sequencing, candidate genes were narrowed to the gene CUFF.40834, which exhibited high identity to ethylene receptor 1 in other model plants. A physiological assay demonstrated that C416 was insensitive to ethylene precursor. Furthermore, the dynamic changes of phytohormone regulatory network in petals at different developmental stages was compared by using RNA-seq. In brief, the ethylene, jasmonic acid (JA), and salicylic acid (SA) signaling pathways were negatively regulated in C416, whereas the brassinosteroid (BR) and cytokinin signaling pathways were positively regulated, and auxin exhibited dual effects on flower senescence in Lotus japonicus. The abscisic acid (ABA) signaling pathway is positively regulated in C416.

CONCLUSION

So far, C416 might be the first reported mutant carrying a mutation in an endogenous ethylene-related gene in Lotus japonicus, rather than through the introduction of exogenous genes by transgenic techniques. A schematic of the flower senescence of Lotus japonicus from the perspective of the phytohormone regulatory network was provided based on transcriptome profiling of petals at different developmental stages. This study is informative for elucidating the molecular mechanism of delayed flower senescence in C416, and lays a foundation for candidate flower senescence gene identification in Lotus japonicus. It also provides another perspective for the improvement of flower longevity in legume plants by heavy-ion beam.

The Dynamism of Transposon Methylation for Plant Development and Stress Adaptation

Plant development processes are regulated by epigenetic alterations that shape nuclear structure, gene expression, and phenotypic plasticity; these alterations can provide the plant with protection from environmental stresses. During plant growth and development, these processes play a significant role in regulating gene expression to remodel chromatin structure. These epigenetic alterations are mainly regulated by transposable elements (TEs) whose abundance in plant genomes results in their interaction with genomes. Thus, TEs are the main source of epigenetic changes and form a substantial part of the plant genome. Furthermore, TEs can be activated under stress conditions, and activated elements cause mutagenic effects and substantial genetic variability. This introduces novel gene functions and structural variation in the insertion sites and primarily contributes to epigenetic modifications. Altogether, these modifications indirectly or directly provide the ability to withstand environmental stresses. In recent years, many studies have shown that TE methylation plays a major role in the evolution of the plant genome through epigenetic process that regulate gene imprinting, thereby upholding genome stability. The induced genetic rearrangements and insertions of mobile genetic elements in regions of active euchromatin contribute to genome alteration, leading to genomic stress. These TE-mediated epigenetic modifications lead to phenotypic diversity, genetic variation, and environmental stress tolerance. Thus, TE methylation is essential for plant evolution and stress adaptation, and TEs hold a relevant military position in the plant genome. High-throughput techniques have greatly advanced the understanding of TE-mediated gene expression and its associations with genome methylation and suggest that controlled mobilization of TEs could be used for crop breeding. However, development application in this area has been limited, and an integrated view of TE function and subsequent processes is lacking. In this review, we explore the enormous diversity and likely functions of the TE repertoire in adaptive evolution and discuss some recent examples of how TEs impact gene expression in plant development and stress adaptation.

A chromosome-scale draft genome sequence of horsegram (Macrotyloma uniflorum)

Horsegram (Macrotyloma uniflorum [Lam.] Verdc.) is an underutilized warm-season diploid legume (2n = 20, 22). Because of its ability to grow under water-deficient and marginal soil conditions, horsegram is a preferred choice in the era of global climate change. In recognition of its potential as a crop species, we generated and analyzed a draft genome sequence for a horsegram variety, HPK-4. Ten chromosome-scale pseudomolecules were created by aligning Illumina scaffold sequences onto a linkage map. The total length of the ten pseudomolecules was 259.2 Mbp, covering 89% of the total length of the assembled sequences. A total of 36,105 genes were predicted on the assembled sequences. Diversity analysis of 89 horsegram accessions by dd-RAD-Seq identified 277 single nucleotide polymorphisms (SNPs), suggesting narrow genetic diversity among the horsegram accessions. This is the first attempt to generate a draft genome sequence of horsegram and will provide a reference for sequence-based analysis of horsegram germplasm.

Genome features of common vetch (Vicia sativa) in natural habitats

Wild plants are often tolerant to biotic and abiotic stresses in their natural environments, whereas domesticated plants such as crops frequently lack such resilience. This difference is thought to be due to the high levels of genome heterozygosity in wild plant populations and the low levels of heterozygosity in domesticated crop species. In this study, common vetch (Vicia sativa) was used as a model to examine this hypothesis. The common vetch genome (2n = 14) was estimated as 1.8 Gb in size. Genome sequencing produced a reference assembly that spanned 1.5 Gb, from which 31,146 genes were predicted. Using this sequence as a reference, 24,118 single nucleotide polymorphisms were discovered in 1243 plants from 12 natural common vetch populations in Japan. Common vetch genomes exhibited high heterozygosity at the population level, with lower levels of heterozygosity observed at specific genome regions. Such patterns of heterozygosity are thought to be essential for adaptation to different environments. The resources generated in this study will provide insights into de novo domestication of wild plants and agricultural enhancement.

Fab Advances in Fabaceae for Abiotic Stress Resilience: From 'Omics' to Artificial Intelligence

Legumes are a better source of proteins and are richer in diverse micronutrients over the nutritional profile of widely consumed cereals. However, when exposed to a diverse range of abiotic stresses, their overall productivity and quality are hugely impacted. Our limited understanding of genetic determinants and novel variants associated with the abiotic stress response in food legume crops restricts its amelioration. Therefore, it is imperative to understand different molecular approaches in food legume crops that can be utilized in crop improvement programs to minimize the economic loss. 'Omics'-based molecular breeding provides better opportunities over conventional breeding for diversifying the natural germplasm together with improving yield and quality parameters. Due to molecular advancements, the technique is now equipped with novel 'omics' approaches such as ionomics, epigenomics, fluxomics, RNomics, glycomics, glycoproteomics, phosphoproteomics, lipidomics, regulomics, and secretomics. Pan-omics-which utilizes the molecular bases of the stress response to identify genes (genomics), mRNAs (transcriptomics), proteins (proteomics), and biomolecules (metabolomics) associated with stress regulation-has been widely used for abiotic stress amelioration in food legume crops. Integration of pan-omics with novel omics approaches will fast-track legume breeding programs. Moreover, artificial intelligence (AI)-based algorithms can be utilized for simulating crop yield under changing environments, which can help in predicting the genetic gain beforehand. Application of machine learning (ML) in quantitative trait loci (QTL) mining will further help in determining the genetic determinants of abiotic stress tolerance in pulses.

The Lotus japonicus AFB6 Gene Is Involved in the Auxin Dependent Root Developmental Program

Auxin is essential for root development, and its regulatory action is exerted at different steps from perception of the hormone up to transcriptional regulation of target genes. In legume plants there is an overlap between the developmental programs governing lateral root and N₂-fixing nodule organogenesis, the latter induced as the result of the symbiotic interaction with rhizobia. Here we report the characterization of a member of the L. japonicus TIR1/AFB auxin receptor family, LjAFB6. A preferential expression of the LjAFB6 gene in the aerial portion of L. japonicus plants was observed. Significant regulation of the expression was not observed during the symbiotic interaction with Mesorhizobium loti and the nodule organogenesis process. In roots, the LjAFB6 expression was induced in response to nitrate supply and was mainly localized in the meristematic regions of both primary and lateral roots. The phenotypic analyses conducted on two independent null mutants indicated a specialized role in the control of primary and lateral root elongation processes in response to auxin, whereas no involvement in the nodulation process was found. We also report the involvement of LjAFB6 in the hypocotyl elongation process and in the control of the expression profile of an auxin-responsive gene.

Effect of Three Types of Ion Beam Irradiation on Gerbera (Gerbera hybrida) In Vitro Shoots with Mutagenesis Efficiency

Gerbera in vitro shoots were irradiated using three types of ion beams with different line energy transfers (LETs) to investigate the effective LET and absorbed doses for mutagenesis. Furthermore, genomic mutation analyses were conducted on the obtained mutants. Survival rate analysis showed a lower lethal dose 50% (LD₅₀) with ion beams with higher LETs. Trait/morphological mutations exhibited changes in the color and shape of petals and male sterility. Irradiation conditions with the highest growth change and trait/morphological mutation rates in each ion were C irradiation at 10 Gy, Ar irradiation at 5 Gy, and Fe irradiation at 5 Gy, with a range of absorbed dose of around LD₅₀ to about 10 Gy lower. The highest trait/morphological mutation rate was 14.1% with Ar irradiation at 5 Gy, which was one of the criteria for ion beam irradiation of gerbera in vitro shoots. Furthermore, the genomic mutation in the flower color, petal shape, and male sterile mutants were confirmed by genotype analysis using Genotyping by Random Amplicon Sequencing-Direct technology. This is the first study to report the efficient production of gerbera mutants that could be analyzed. Our findings may lead to more efficient gerbera mutant production and analysis technology.

Genome-Wide Identification and Analysis of the Polycomb Group Family in Medicago truncatula

Polycomb group (PcG) proteins, which are important epigenetic regulators, play essential roles in the regulatory networks involved in plant growth, development, and environmental stress responses. Currently, as far as we know, no comprehensive and systematic study has been carried out on the PcG family in Medicago truncatula. In the present study, we identified 64 PcG genes with distinct gene structures from the M. truncatula genome. All of the PcG genes were distributed unevenly over eight chromosomes, of which 26 genes underwent gene duplication. The prediction of protein interaction network indicated that 34 M. truncatula PcG proteins exhibited protein-protein interactions, and MtMSI1;4 and MtVRN2 had the largest number of protein-protein interactions. Based on phylogenetic analysis, we divided 375 PcG proteins from 27 species into three groups and nine subgroups. Group I and Group III were composed of five components from the PRC1 complex, and Group II was composed of four components from the PRC2 complex. Additionally, we found that seven PcG proteins in M. truncatula were closely related to the corresponding proteins of Cicer arietinum. Syntenic analysis revealed that PcG proteins had evolved more conservatively in dicots than in monocots. M. truncatula had the most collinearity relationships with Glycine max (36 genes), while collinearity with three monocots was rare (eight genes). The analysis of various types of expression data suggested that PcG genes were involved in the regulation and response process of M. truncatula in multiple developmental stages, in different tissues, and for various environmental stimuli. Meanwhile, many differentially expressed genes (DEGs) were identified in the RNA-seq data, which had potential research value in further studies on gene function verification. These findings provide novel and detailed information on the M. truncatula PcG family, and in the future it would be helpful to carry out related research on the PcG family in other legumes.

Use of the rhizobial type III effector gene nopP to improve Agrobacterium rhizogenes-mediated transformation of Lotus japonicus

BACKGROUND

Protocols for Agrobacterium rhizogenes-mediated hairy root transformation of the model legume Lotus japonicus have been established previously. However, little efforts were made in the past to quantify and improve the transformation efficiency. Here, we asked whether effectors (nodulation outer proteins) of the nodule bacterium Sinorhizobium sp. NGR234 can promote hairy root transformation of L. japonicus. The co-expressed red fluorescent protein DsRed1 was used for visualization of transformed roots and for estimation of the transformation efficiency.

RESULTS

Strong induction of hairy root formation was observed when A. rhizogenes strain LBA9402 was used for L. japonicus transformation. Expression of the effector gene nopP in L. japonicus roots resulted in a significantly increased transformation efficiency while nopL, nopM, and nopT did not show such an effect. In nopP expressing plants, more than 65% of the formed hairy roots were transgenic as analyzed by red fluorescence emitted by co-transformed DsRed1. A nodulation experiment indicated that nopP expression did not obviously affect the symbiosis between L. japonicus and Mesorhizobium loti.

CONCLUSION

We have established a novel protocol for hairy root transformation of L. japonicus. The use of A. rhizogenes LBA9402 carrying a binary vector containing DsRed1 and nopP allowed efficient formation and identification of transgenic roots.

Insights into the evolution of symbiosis gene copy number and distribution from a chromosome-scale Lotus japonicus Gifu genome sequence

Lotus japonicus is a herbaceous perennial legume that has been used extensively as a genetically tractable model system for deciphering the molecular genetics of symbiotic nitrogen fixation. Our aim is to improve the L. japonicus reference genome sequence, which has so far been based on Sanger and Illumina sequencing reads from the L. japonicus accession MG-20 and contained a large fraction of unanchored contigs. Here, we use long PacBio reads from L. japonicus Gifu combined with Hi-C data and new high-density genetic maps to generate a high-quality chromosome-scale reference genome assembly for L. japonicus. The assembly comprises 554 megabases of which 549 were assigned to six pseudomolecules that appear complete with telomeric repeats at their extremes and large centromeric regions with low gene density. The new L. japonicus Gifu reference genome and associated expression data represent valuable resources for legume functional and comparative genomics. Here, we provide a first example by showing that the symbiotic islands recently described in Medicago truncatula do not appear to be conserved in L. japonicus.

Draft genome sequence of the pulse crop blackgram [Vigna mungo (L.) Hepper] reveals potential R-genes

Blackgram [Vigna mungo (L.) Hepper] (2n = 2x = 22), an important Asiatic legume crop, is a major source of dietary protein for the predominantly vegetarian population. Here we construct a draft genome sequence of blackgram, for the first time, by employing hybrid genome assembly with Illumina reads and third generation Oxford Nanopore sequencing technology. The final de novo whole genome of blackgram is ~ 475 Mb (82% of the genome) and has maximum scaffold length of 6.3 Mb with scaffold N50 of 1.42 Mb. Genome analysis identified 42,115 genes with mean coding sequence length of 1131 bp. Around 80.6% of predicted genes were annotated. Nearly half of the assembled sequence is composed of repetitive elements with retrotransposons as major (47.3% of genome) transposable elements, whereas, DNA transposons made up only 2.29% of the genome. A total of 166,014 SSRs, including 65,180 compound SSRs, were identified and primer pairs for 34,816 SSRs were designed. Out of the 33,959 proteins, 1659 proteins showed presence of R-gene related domains. KIN class was found in majority of the proteins (905) followed by RLK (239) and RLP (188). The genome sequence of blackgram will facilitate identification of agronomically important genes and accelerate the genetic improvement of blackgram.

The genome of a wild Medicago species provides insights into the tolerant mechanisms of legume forage to environmental stress

BACKGROUND

Medicago ruthenica, a wild and perennial legume forage widely distributed in semi-arid grasslands, is distinguished by its outstanding tolerance to environmental stress. It is a close relative of commonly cultivated forage of alfalfa (Medicago sativa). The high tolerance of M. ruthenica to environmental stress makes this species a valuable genetic resource for understanding and improving traits associated with tolerance to harsh environments.

RESULTS

We sequenced and assembled genome of M. ruthenica using an integrated approach, including PacBio, Illumina, 10×Genomics, and Hi-C. The assembled genome was 904.13 Mb with scaffold N50 of 99.39 Mb, and 50,162 protein-coding genes were annotated. Comparative genomics and transcriptomic analyses were used to elucidate mechanisms underlying its tolerance to environmental stress. The expanded FHY3/FAR1 family was identified to be involved in tolerance of M. ruthenica to drought stress. Many genes involved in tolerance to abiotic stress were retained in M. ruthenica compared to other cultivated Medicago species. Hundreds of candidate genes associated with drought tolerance were identified by analyzing variations in single nucleotide polymorphism using accessions of M. ruthenica with varying tolerance to drought. Transcriptomic data demonstrated the involvements of genes related to transcriptional regulation, stress response, and metabolic regulation in tolerance of M. ruthenica.

CONCLUSIONS

We present a high-quality genome assembly and identification of drought-related genes in the wild species of M. ruthenica, providing a valuable resource for genomic studies on perennial legume forages.

Phospholipases C and D and Their Role in Biotic and Abiotic Stresses

Plants, as sessile organisms, have adapted a fine sensing system to monitor environmental changes, therefore allowing the regulation of their responses. As the interaction between plants and environmental changes begins at the surface, these changes are detected by components in the plasma membrane, where a molecule receptor generates a lipid signaling cascade via enzymes, such as phospholipases (PLs). Phospholipids are the key structural components of plasma membranes and signaling cascades. They exist in a wide range of species and in different proportions, with conversion processes that involve hydrophilic enzymes, such as phospholipase-C (PLC), phospholipase-D (PLD), and phospholipase-A (PLA). Hence, it is suggested that PLC and PLD are highly conserved, compared to their homologous genes, and have formed clusters during their adaptive history. Additionally, they generate responses to different functions in accordance with their protein structure, which should be reflected in specific signal transduction responses to environmental stress conditions, including innate immune responses. This review summarizes the phospholipid systems associated with signaling pathways and the innate immune response.

BAC- and oligo-FISH mapping reveals chromosome evolution among Vigna angularis, V. unguiculata, and Phaseolus vulgaris

Cytogenomic resources have accelerated synteny and chromosome evolution studies in plant species, including legumes. Here, we established the first cytogenetic map of V. angularis (Va, subgenus Ceratotropis) and compared this new map with those of V. unguiculata (Vu, subgenus Vigna) and P. vulgaris (Pv) by BAC-FISH and oligopainting approaches. We mapped 19 Vu BACs and 35S rDNA probes to the 11 chromosome pairs of Va, Vu, and Pv. Vigna angularis shared a high degree of macrosynteny with Vu and Pv, with five conserved syntenic chromosomes. Additionally, we developed two oligo probes (Pv2 and Pv3) used to paint Vigna orthologous chromosomes. We confirmed two reciprocal translocations (chromosomes 2 and 3 and 1 and 8) that have occurred after the Vigna and Phaseolus divergence (~9.7 Mya). Besides, two inversions (2 and 4) and one translocation (1 and 5) have occurred after Vigna and Ceratotropis subgenera separation (~3.6 Mya). We also observed distinct oligopainting patterns for chromosomes 2 and 3 of Vigna species. Both Vigna species shared similar major rearrangements compared to Pv: one translocation (2 and 3) and one inversion (chromosome 3). The sequence synteny identified additional inversions and/or intrachromosomal translocations involving pericentromeric regions of both orthologous chromosomes. We propose chromosomes 2 and 3 as hotspots for chromosomal rearrangements and de novo centromere formation within and between Vigna and Phaseolus. Our BAC- and oligo-FISH mapping contributed to physically trace the chromosome evolution of Vigna and Phaseolus and its application in further studies of both genera.

The genome of Medicago polymorpha provides insights into its edibility and nutritional value as a vegetable and forage legume

Medicago polymorpha is a nutritious and palatable forage and vegetable plant that also fixes nitrogen. Here, we reveal the chromosome-scale genome sequence of M. polymorpha using an integrated approach including Illumina, PacBio and Hi-C technologies. We combined PacBio full-length RNA-seq, metabolomic analysis, structural anatomy analysis and related physiological indexes to elucidate the important agronomic traits of M. polymorpha for forage and vegetable usage. The assembled M. polymorpha genome consisted of 457.53 Mb with a long scaffold N50 of 57.72 Mb, and 92.92% (441.83 Mb) of the assembly was assigned to seven pseudochromosomes. Comparative genomic analysis revealed that expansion and contraction of the photosynthesis and lignin biosynthetic gene families, respectively, led to enhancement of nutritious compounds and reduced lignin biosynthesis in M. polymorpha. In addition, we found that several positively selected nitrogen metabolism-related genes were responsible for crude protein biosynthesis. Notably, the metabolomic results revealed that a large number of flavonoids, vitamins, alkaloids, and terpenoids were enriched in M. polymorpha. These results imply that the decreased lignin content but relatively high nutrient content of M. polymorpha enhance its edibility and nutritional value as a forage and vegetable. Our genomic data provide a genetic basis that will accelerate functional genomic and breeding research on M. polymorpha as well as other Medicago and legume plants.

Mutants of Lotus japonicus deficient in flavonoid biosynthesis

Spatiotemporal features of anthocyanin accumulation in a model legume Lotus japonicus (Regel) K.Larsen were elucidated to develop criteria for the genetic analysis of flavonoid biosynthesis. Artificial mutants and wild accessions, with lower anthocyanin accumulation in the stem than the standard wild type (B-129 'Gifu'), were obtained by ethyl methanesulfonate (EMS) mutagenesis and from a collection of wild-grown variants, respectively. The loci responsible for the green stem of the mutants were named as VIRIDICAULIS (VIC). Genetic and chemical analysis identified two loci, namely, VIC1 and VIC2, required for the production of both anthocyanins and proanthocyanidins (condensed tannins), and two loci, namely, VIC3 and VIC4, required for the steps specific to anthocyanin biosynthesis. A mutation in VIC5 significantly reduced the anthocyanin accumulation. These mutants will serve as a useful system for examining the effects of anthocyanins and proanthocyanidins on the interactions with herbivorous pests, pathogenic microorganisms and nitrogen-fixing symbiotic bacteria, Mesorhizobium loti.

Transcriptome analysis reveals major transcriptional changes during regrowth after mowing of red clover (Trifolium pratense)

BACKGROUND

Red clover (Trifolium pratense) is globally used as a fodder plant due its high nutritional value and soil improving qualities. In response to mowing, red clover exhibits specific morphological traits to compensate the loss of biomass. The morphological reaction is well described, but the underlying molecular mechanisms and its role for plants grown in the field are unclear.

RESULTS

Here, we characterize the global transcriptional response to mowing of red clover by comparing plants grown under greenhouse conditions with plants growing on agriculturally used fields. Unexpectedly, we found that biotic and abiotic stress related changes of plants grown in the field overlay their regrowth related transcriptional changes and characterized transcription related protein families involved in these processes. Further, we can show that gibberellins, among other phytohormones, also contribute to the developmental processes related to regrowth after biomass-loss.

CONCLUSIONS

Our findings show that massive biomass loss triggers less transcriptional changes in field grown plants than their struggle with biotic and abiotic stresses and that gibberellins also play a role in the developmental program related to regrowth after mowing in red clover. Our results provide first insights into the physiological and developmental processes of mowing on red clover and may serve as a base for red clover yield improvement.