Establishment of the Lotus japonicus Gene Expression Atlas (LjGEA) and its use to explore legume seed maturation

Single-cell transcriptome atlases of soybean root and mature nodule reveal new regulatory programs that control the nodulation process

The soybean root system is complex. In addition to being composed of various cell types, the soybean root system includes the primary root, the lateral roots, and the nodule, an organ in which mutualistic symbiosis with N-fixing rhizobia occurs. A mature soybean root nodule is characterized by a central infection zone where atmospheric nitrogen is fixed and assimilated by the symbiont, resulting from the close cooperation between the plant cell and the bacteria. To date, the transcriptome of individual cells isolated from developing soybean nodules has been established, but the transcriptomic signatures of cells from the mature soybean nodule have not yet been characterized. Using single-nucleus RNA-seq and Molecular Cartography technologies, we precisely characterized the transcriptomic signature of soybean root and mature nodule cell types and revealed the co-existence of different sub-populations of B. diazoefficiens-infected cells in the mature soybean nodule, including those actively involved in nitrogen fixation and those engaged in senescence. Mining of the single-cell-resolution nodule transcriptome atlas and the associated gene co-expression network confirmed the role of known nodulation-related genes and identified new genes that control the nodulation process. For instance, we functionally characterized the role of GmFWL3, a plasma membrane microdomain-associated protein that controls rhizobial infection. Our study reveals the unique cellular complexity of the mature soybean nodule and helps redefine the concept of cell types when considering the infection zone of the soybean nodule.

An integrated transcriptome mapping the regulatory network of coding and long non-coding RNAs provides a genomics resource in chickpea

Large-scale transcriptome analysis can provide a systems-level understanding of biological processes. To accelerate functional genomic studies in chickpea, we perform a comprehensive transcriptome analysis to generate full-length transcriptome and expression atlas of protein-coding genes (PCGs) and long non-coding RNAs (lncRNAs) from 32 different tissues/organs via deep sequencing. The high-depth RNA-seq dataset reveal expression dynamics and tissue-specificity along with associated biological functions of PCGs and lncRNAs during development. The coexpression network analysis reveal modules associated with a particular tissue or a set of related tissues. The components of transcriptional regulatory networks (TRNs), including transcription factors, their cognate cis-regulatory motifs, and target PCGs/lncRNAs that determine developmental programs of different tissues/organs, are identified. Several candidate tissue-specific and abiotic stress-responsive transcripts associated with quantitative trait loci that determine important agronomic traits are also identified. These results provide an important resource to advance functional/translational genomic and genetic studies during chickpea development and environmental conditions.

A full-length transcriptome and expression atlas of protein-coding genes and long non-coding RNAs is generated in chickpea. Components of transcriptional regulatory networks and candidate tissue-specific transcripts associated with quantitative trait loci are identified.

Exploring Chickpea Germplasm Diversity for Broadening the Genetic Base Utilizing Genomic Resourses

Legume crops provide significant nutrition to humans as a source of protein, omega-3 fatty acids as well as specific macro and micronutrients. Additionally, legumes improve the cropping environment by replenishing the soil nitrogen content. Chickpeas are the second most significant staple legume food crop worldwide behind dry bean which contains 17%–24% protein, 41%–51% carbohydrate, and other important essential minerals, vitamins, dietary fiber, folate, β-carotene, anti-oxidants, micronutrients (phosphorus, calcium, magnesium, iron, and zinc) as well as linoleic and oleic unsaturated fatty acids. Despite these advantages, legumes are far behind cereals in terms of genetic improvement mainly due to far less effort, the bottlenecks of the narrow genetic base, and several biotic and abiotic factors in the scenario of changing climatic conditions. Measures are now called for beyond conventional breeding practices to strategically broadening of narrow genetic base utilizing chickpea wild relatives and improvement of cultivars through advanced breeding approaches with a focus on high yield productivity, biotic and abiotic stresses including climate resilience, and enhanced nutritional values. Desirable donors having such multiple traits have been identified using core and mini core collections from the cultivated gene pool and wild relatives of Chickpea. Several methods have been developed to address cross-species fertilization obstacles and to aid in inter-specific hybridization and introgression of the target gene sequences from wild Cicer species. Additionally, recent advances in “Omics” sciences along with high-throughput and precise phenotyping tools have made it easier to identify genes that regulate traits of interest. Next-generation sequencing technologies, whole-genome sequencing, transcriptomics, and differential genes expression profiling along with a plethora of novel techniques like single nucleotide polymorphism exploiting high-density genotyping by sequencing assays, simple sequence repeat markers, diversity array technology platform, and whole-genome re-sequencing technique led to the identification and development of QTLs and high-density trait mapping of the global chickpea germplasm. These altogether have helped in broadening the narrow genetic base of chickpeas.

A novel SCARECROW-LIKE3 transcription factor LjGRAS36 in Lotus japonicus regulates the development of arbuscular mycorrhizal symbiosis

The symbiosis with arbuscular mycorrhizal (AM) fungi improves plants' nutrient uptake. During this process, transcription factors have been highlighted to play crucial roles. Members of the GRAS transcription factor gene family have been reported involved in AM symbiosis, but little is known about SCARECROW-LIKE3 (SCL3) genes belonging to this family in Lotus japonicus. In this study, 67 LjGRAS genes were identified from the L. japonicus genome, seven of which were clustered in the SCL3 group. Three of the seven LjGRAS genes expression levels were upregulated by AM fungal inoculation, and our biochemical results showed that the expression of LjGRAS36 was specifically induced by AM colonization. Functional loss of LjGRAS36 in mutant ljgras36 plants exhibited a significantly reduced mycorrhizal colonization rate and arbuscular size. Transcriptome analysis showed a deficiency of LjGRAS36 led to the dysregulation of the gibberellic acid signal pathway associated with AM symbiosis. Together, this study provides important insights for understanding the important potential function of SCL3 genes in regulating AM symbiotic development.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s12298-022-01161-z.

Processing of NODULE INCEPTION controls the transition to nitrogen fixation in root nodules

[Figure: see text].

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.

Multi-species relationships in legume roots: From pairwise legume-symbiont interactions to the plant - microbiome - soil continuum

Mutualistic relationships of legume plants with, either bacteria (like rhizobia) or fungi (like arbuscular mycorrhizal fungi), have been investigated intensively, usually as bi-partite interactions. However, diverse symbiotic interactions take place simultaneously or sequentially under field conditions. Their collective, but not additive, contribution to plant growth and performance remains hard to predict, and appears to be furthermore affected by crop species and genotype, non-symbiotic microbial interactions and environmental variables. The challenge is: (i) to unravel the complex overlapping mechanisms that operate between the microbial symbionts as well as between them, their hosts and the rhizosphere (ii) to understand the dynamics of the respective mechanisms in evolutionary and ecological terms. The target for agriculture, food security and the environment, is to use this insight as a solid basis for developing new integrated technologies, practices and strategies for the efficient use of beneficial microbes in legumes and other plants. We review recent advances in our understanding of the symbiotic interactions in legumes roots brought about with the aid of molecular and bioinformatics tools. We go through single symbiont-host interactions, proceed to tripartite symbiont-host interactions, appraise interactions of symbiotic and associative microbiomes with plants in the root-rhizoplane-soil continuum of habitats and end up by examining attempts to validate community ecology principles in the legume-microbe-soil biosystem.

The functional characterization of LjNRT2.4 indicates a novel, positive role of nitrate for an efficient nodule N2 -fixation activity

Atmospheric nitrogen (N₂₎ -fixing nodules are formed on the roots of legume plants as result of the symbiotic interaction with rhizobia. Nodule functioning requires high amounts of carbon and energy, and therefore legumes have developed finely tuned mechanisms to cope with changing external environmental conditions, including nutrient availability and flooding. The investigation of the role of nitrate as regulator of the symbiotic N₂ fixation has been limited to the inhibitory effects exerted by high external concentrations on nodule formation, development and functioning. We describe a nitrate-dependent route acting at low external concentrations that become crucial in hydroponic conditions to ensure an efficient nodule functionality. Combined genetic, biochemical and molecular studies are used to unravel the novel function of the LjNRT2.4 gene. Two independent null mutants are affected by the nitrate content of nodules, consistent with LjNRT2.4 temporal and spatial profiles of expression. The reduced nodular nitrate content is associated to a strong reduction of nitrogenase activity and a severe N-starvation phenotype observed under hydroponic conditions. We also report the effects of the mutations on the nodular nitric oxide (NO) production and content. We discuss the involvement of LjNRT2.4 in a nitrate-NO respiratory chain taking place in the N₂ -fixing nodules.

Identification of Alternative Mitochondrial Electron Transport Pathway Components in Chickpea Indicates a Differential Response to Salinity Stress between Cultivars

All plants contain an alternative electron transport pathway (AP) in their mitochondria, consisting of the alternative oxidase (AOX) and type 2 NAD(P)H dehydrogenase (ND) families, that are thought to play a role in controlling oxidative stress responses at the cellular level. These alternative electron transport components have been extensively studied in plants like Arabidopsis and stress inducible isoforms identified, but we know very little about them in the important crop plant chickpea. Here we identify AP components in chickpea (Cicer arietinum) and explore their response to stress at the transcript level. Based on sequence similarity with the functionally characterized proteins of Arabidopsis thaliana, five putative internal (matrix)-facing NAD(P)H dehydrogenases (CaNDA1-4 and CaNDC1) and four putative external (inter-membrane space)-facing NAD(P)H dehydrogenases (CaNDB1-4) were identified in chickpea. The corresponding activities were demonstrated for the first time in purified mitochondria of chickpea leaves and roots. Oxidation of matrix NADH generated from malate or glycine in the presence of the Complex I inhibitor rotenone was high compared to other plant species, as was oxidation of exogenous NAD(P)H. In leaf mitochondria, external NADH oxidation was stimulated by exogenous calcium and external NADPH oxidation was essentially calcium dependent. However, in roots these activities were low and largely calcium independent. A salinity experiment with six chickpea cultivars was used to identify salt-responsive alternative oxidase and NAD(P)H dehydrogenase gene transcripts in leaves from a three-point time series. An analysis of the Na:K ratio and Na content separated these cultivars into high and low Na accumulators. In the high Na accumulators, there was a significant up-regulation of CaAOX1, CaNDB2, CaNDB4, CaNDA3 and CaNDC1 in leaf tissue under long term stress, suggesting the formation of a stress-modified form of the mitochondrial electron transport chain (mETC) in leaves of these cultivars. In particular, stress-induced expression of the CaNDB2 gene showed a striking positive correlation with that of CaAOX1 across all genotypes and time points. The coordinated salinity-induced up-regulation of CaAOX1 and CaNDB2 suggests that the mitochondrial alternative pathway of respiration is an important facet of the stress response in chickpea, in high Na accumulators in particular, despite high capacities for both of these activities in leaf mitochondria of non-stressed chickpeas.

A High-Quality Genome Sequence of Model Legume Lotus japonicus (MG-20) Provides Insights into the Evolution of Root Nodule Symbiosis

Lotus japonicus is an important model legume for studying symbiotic nitrogen fixation as well as plant development. A genomic sequence of L. japonicus (MG20) has been available for more than ten years. However, the low quality of the genome limits its application in functional genomic studies. Therefore, it is necessary to assemble high-quality chromosome sequences of L. japonicus using new sequencing technology to facilitate the study of functional genomics. In this report, we used the third-generation sequencing combined with the Illumina HiSeq platform to sequence the genome of L. japonicus (MG20). We obtained 544 Mb of genomic sequence using third-generation assembly. Based on sequence analysis, 357 Mb of repeats, 28,251 genes, 626 tRNAs, 1409 rRNAs, and 1233 pseudogenes were predicted in the genome. A total of 27,991 genes were annotated into databases. Compared to the previously published data, the new genome database contains complete L. japonicus sequences in the proper order and orientation with a contig N50 2.81Mb and an excellent genome coverage, which provides more accurate genome information and more precise assembly for functional genomic study.

Extreme genetic signatures of local adaptation during Lotus japonicus colonization of Japan

Colonization of new habitats is expected to require genetic adaptations to overcome environmental challenges. Here, we use full genome re-sequencing and extensive common garden experiments to investigate demographic and selective processes associated with colonization of Japan by Lotus japonicus over the past ~20,000 years. Based on patterns of genomic variation, we infer the details of the colonization process where L. japonicus gradually spread from subtropical conditions to much colder climates in northern Japan. We identify genomic regions with extreme genetic differentiation between northern and southern subpopulations and perform population structure-corrected association mapping of phenotypic traits measured in a common garden. Comparing the results of these analyses, we find that signatures of extreme subpopulation differentiation overlap strongly with phenotype association signals for overwintering and flowering time traits. Our results provide evidence that these traits were direct targets of selection during colonization and point to associated candidate genes.

Local adaptation contributes to plant colonization across extreme environmental gradients. Here, the authors reconstruct the colonization history of Lotus japonicus in Japan and identify extreme genetic signatures of local adaptation to a cold climate using genome resequencing and common garden experiments.

Transcriptome analysis of the differential effect of the NADPH oxidase gene RbohB in Phaseolus vulgaris roots following Rhizobium tropici and Rhizophagus irregularis inoculation

BACKGROUND

Reactive oxygen species (ROS) are generated by NADPH oxidases known as respiratory burst oxidase homologs (RBOHs) in plants. ROS regulate various cellular processes, including the mutualistic interactions between legumes and nitrogen-fixing bacteria or arbuscular mycorrhizal (AM) fungi. Rboh is a multigene family comprising nine members (RbohA-I) in common bean (Phaseolus vulgaris). The RNA interference-mediated silencing of RbohB (PvRbohB-RNAi) in this species diminished its ROS production and greatly impaired nodulation. By contrast, the PvRbohB-RNAi transgenic roots showed early hyphal root colonization with enlarged fungal hypopodia; therefore, we proposed that PvRbohB positively regulates rhizobial infection (Rhizobium tropici) and inhibits AM colonization by Rhizophagus irregularis in P. vulgaris.

RESULTS

To corroborate this hypothesis, an RNA-Seq transcriptomic analysis was performed to identify the differentially expressed genes in the PvRbohB-RNAi roots inoculated with Rhizobium tropici or Rhizophagus irregularis. We found that, in the early stages, root nodule symbioses generated larger changes of the transcriptome than did AM symbioses in P. vulgaris. Genes related to ROS homeostasis and cell wall flexibility were markedly upregulated in the early stages of rhizobial colonization, but not during AM colonization. Compared with AM colonization, the rhizobia induced the expression of a greater number of genes encoding enzymes involved in the metabolism of auxins, cytokinins, and ethylene, which were typically repressed in the PvRbohB-RNAi roots.

CONCLUSIONS

Our research provides substantial insights into the genetic interaction networks in the early stages of rhizobia and AM symbioses with P. vulgaris, as well as the differential roles that RbohB plays in processes related to ROS scavenging, cell wall remodeling, and phytohormone homeostasis during nodulation and mycorrhization in this legume.

Insight Into the Prospects for the Improvement of Seed Starch in Legume-A Review

In addition to proteins and/or oils, mature seeds of most legume crops contain important carbohydrate components, including starches and sugars. Starch is also an essential nutritional component of human and animal diets and has various food and non-food industrial applications. Starch is a primary insoluble polymeric carbohydrate produced by higher plants and consists of amylose and amylopectin as a major fraction. Legume seeds are an affordable source of not only protein but also the starch, which has an advantage of being resistant starch compared with cereal, root, and tuber starch. For these reasons, legume seeds form a good source of resistant starch-rich healthy food with a high protein content and can be utilized in various food applications. The genetics and molecular details of starch and other carbohydrate components are well studied in cereal crops but have received little attention in legumes. In order to improve legume starch content, quality, and quantity, it is necessary to understand the genetic and molecular factors regulating carbohydrate metabolism in legume crops. In this review, we assessed the current literature reporting the genetic and molecular basis of legume carbohydrate components, primarily focused on seed starch content. We provided an overview of starch biosynthesis in the heterotrophic organs, the chemical composition of major consumable legumes, the factors influencing starch digestibility, and advances in the genetic, transcriptomic, and metabolomic studies in important legume crops. Further, we discussed breeding and biotechnological approaches for the improvement of the starch composition in major legume crops. The information reviewed in this study will be helpful in facilitating the food and non-food applications of legume starch and provide economic benefits to farmers and industries.

Data Integration in Poplar: 'Omics Layers and Integration Strategies

Populus trichocarpa is an important biofuel feedstock that has been the target of extensive research and is emerging as a model organism for plants, especially woody perennials. This research has generated several large ‘omics datasets. However, only few studies in Populus have attempted to integrate various data types. This review will summarize various ‘omics data layers, focusing on their application in Populus species. Subsequently, network and signal processing techniques for the integration and analysis of these data types will be discussed, with particular reference to examples in Populus.

Gene Expression Maps in Plants: Current State and Prospects

For many years, progress in the identification of gene functions has been based on classical genetic approaches. However, considerable recent omics developments have brought to the fore indirect but high-resolution methods of gene function identification such as transcriptomics, proteomics, and metabolomics. A transcriptome map is a powerful source of functional information and the result of the genome-wide expression analysis of a broad sampling of tissues and/or organs from different developmental stages and/or environmental conditions. In plant science, the application of transcriptome maps extends from the inference of gene regulatory networks to evolutionary studies. However, only some of these data have been integrated into databases, thus enabling analyses to be conducted without raw data; without this integration, extensive data preprocessing is required, which limits data usability. In this review, we summarize the state of plant transcriptome maps, analyze the problems associated with the combined analysis of large-scale data from various studies, and outline possible solutions to these problems.

A member of the ALOG gene family has a novel role in regulating nodulation in Lotus japonicus

Legumes can control the number of symbiotic nodules that form on their roots, thus balancing nitrogen assimilation and energy consumption. Two major pathways participate in nodulation: the Nod factor (NF) signaling pathway which involves recognition of rhizobial bacteria by root cells and promotion of nodulation, and the autoregulation of nodulation (AON) pathway which involves long-distance negative feedback between roots and shoots. Although a handful of genes have a clear role in the maintenance of nodule number, additional unknown factors may also be involved in this process. Here, we identify a novel function for a Lotus japonicus ALOG (Arabidopsis LSH1 and Oryza G1) family member, LjALOG1, involved in positively regulating nodulation. LjALOG1 expression increased substantially after inoculation with rhizobia, with high levels of expression in whole nodule primordia and in the base of developing nodules. The ljalog1 mutants, which have an insertion of the LORE1 retroelement in LjALOG1, had significantly fewer nodules compared with wild type, along with increased expression of LjCLE-RS1 (L. japonicus CLE Root Signal 1), which encodes a nodulation suppressor in the AON pathway. In summary, our findings identified a novel factor that participates in controlling nodulation, possibly by suppressing the AON pathway.

Lotus SHAGGY-like kinase 1 is required to suppress nodulation in Lotus japonicus

Glycogen synthase kinase/SHAGGY-like kinases (SKs) are a highly conserved family of signaling proteins that participate in many developmental, cell-differentiation, and metabolic signaling pathways in plants and animals. Here, we investigate the involvement of SKs in legume nodulation, a process requiring the integration of multiple signaling pathways. We describe a group of SKs in the model legume Lotus japonicus (LSKs), two of which respond to inoculation with the symbiotic nitrogen-fixing bacterium Mesorhizobium loti. RNAi knock-down plants and an insertion mutant for one of these genes, LSK1, display increased nodulation. Ηairy-root lines overexpressing LSK1 form only marginally fewer mature nodules compared with controls. The expression levels of genes involved in the autoregulation of nodulation (AON) mechanism are affected in LSK1 knock-down plants at low nitrate levels, both at early and late stages of nodulation. At higher levels of nitrate, these same plants show the opposite expression pattern of AON-related genes and lose the hypernodulation phenotype. Our findings reveal an additional role for the versatile SK gene family in integrating the signaling pathways governing legume nodulation, and pave the way for further study of their functions in legumes.

Transcriptomic Analysis With the Progress of Symbiosis in 'Crack-Entry' Legume Arachis hypogaea Highlights Its Contrast With 'Infection Thread' Adapted Legumes

In root-nodule symbiosis, rhizobial invasion and nodule organogenesis is host controlled. In most legumes, rhizobia enter through infection threads and nodule primordium in the cortex is induced from a distance. But in dalbergoid legumes like Arachis hypogaea, rhizobia directly invade cortical cells through epidermal cracks to generate the primordia. Herein, we report the transcriptional dynamics with the progress of symbiosis in A. hypogaea at 1 day postinfection (dpi) (invasion), 4 dpi (nodule primordia), 8 dpi (spread of infection in nodule-like structure), 12 dpi (immature nodules containing rod-shaped rhizobia), and 21 dpi (mature nodules with spherical symbiosomes). Expression of putative ortholog of symbiotic genes in 'crack entry' legume A. hypogaea was compared with infection thread-adapted model legumes. The contrasting features were i) higher expression of receptors like LYR3 and EPR3 as compared with canonical Nod factor receptors, ii) late induction of transcription factors like NIN and NSP2 and constitutive high expression of ERF1, EIN2, bHLH476, and iii) induction of divergent pathogenesis-responsive PR-1 genes. Additionally, symbiotic orthologs of SymCRK, ROP6, RR9, SEN1, and DNF2 were not detectable and microsynteny analysis indicated the absence of a RPG homolog in diploid parental genomes of A. hypogaea. The implications are discussed and a molecular framework that guides crack-entry symbiosis in A. hypogaea is proposed.

Genomic structure and expression of alternative oxidase genes in legumes

Mitochondria isolated from chickpea (Cicer arietinum) possess substantial alternative oxidase (AOX) activity, even in non-stressed plants, and one or two AOX protein bands were detected immunologically, depending on the organ. Four different AOX isoforms were identified in the chickpea genome: CaAOX1 and CaAOX2A, B and D. CaAOX2A was the most highly expressed form and was strongly expressed in photosynthetic tissues, whereas CaAOX2D was found in all organs examined. These results are very similar to those of previous studies with soybean and siratro. Searches of available databases showed that this pattern of AOX genes and their expression was common to at least 16 different legume species. The evolution of the legume AOX gene family is discussed, as is the in vivo impact of an inherently high AOX capacity in legumes on growth and responses to environmental stresses.

Analysis of genes encoding seed storage proteins (SSPs) in chickpea (Cicer arietinum L.) reveals co-expressing transcription factors and a seed-specific promoter

Improvement of the quality and quantity of chickpea seed protein can be greatly facilitated by an understanding of the genic organization and the genetic architecture of the genes encoding seed storage proteins (SSPs). The aim of this study was to provide a comprehensive analysis of the chickpea SSP genes, putative co-expressing transcription factors (TFs), and to identify a seed-specific SSP gene promoter. A genome-wide identification of SSP genes in chickpea led to the identification of 21 non-redundant SSP encoding genes located on 6 chromosomes. Phylogenetic analysis grouped SSP genes into 3 subgroups where members within the same clade demonstrated similar motif composition and intron-exon organization. Tandem duplications were identified to be the major contributors to the expansion of the SSP gene family in chickpea. Co-expression analysis revealed 14 TFs having expression profiles similar to the SSP genes that included members of important TF families that are known to regulate seed development. Expression analysis of SSP genes and TFs revealed significantly higher expression in late stages of seed development as well as in high seed protein content (HPC) genotypes. In silico analysis of the promoter regions of the SSP encoding genes revealed several seed-specific cis-regulatory elements such as RY repeats, ACGT motifs, CAANTG, and GCN4. A candidate promoter was analyzed for seed specificity by generating stable transgenics in Arabidopsis. Overall, this study provides a useful resource to explore the regulatory networks involved in SSP synthesis and/or accumulation for utilization in developing nutritionally improved chickpea genotypes.

Legume Cytosolic and Plastid Acetyl-Coenzyme-A Carboxylase Genes Differ by Evolutionary Patterns and Selection Pressure Schemes Acting before and after Whole-Genome Duplications

Acetyl-coenzyme A carboxylase (ACCase, E.C.6.4.1.2) catalyzes acetyl-coenzyme A carboxylation to malonyl coenzyme A. Plants possess two distinct ACCases differing by cellular compartment and function. Plastid ACCase contributes to de novo fatty acid synthesis, whereas cytosolic enzyme to the synthesis of very long chain fatty acids, phytoalexins, flavonoids, and anthocyanins. The narrow leafed lupin (Lupinus angustifolius L.) represents legumes, a plant family which evolved by whole-genome duplications (WGDs). The study aimed on the contribution of these WGDs to the multiplication of ACCase genes and their further evolutionary patterns. The molecular approach involved bacterial artificial chromosome (BAC) library screening, fluorescent in situ hybridization, linkage mapping, and BAC sequencing. In silico analysis encompassed sequence annotation, comparative mapping, selection pressure calculation, phylogenetic inference, and gene expression profiling. Among sequenced legumes, the highest number of ACCase genes was identified in lupin and soybean. The most abundant plastid ACCase subunit genes were accB. ACCase genes in legumes evolved by WGDs, evidenced by shared synteny and Bayesian phylogenetic inference. Transcriptional activity of almost all copies was confirmed. Gene duplicates were conserved by strong purifying selection, however, positive selection occurred in Arachis (accB2) and Lupinus (accC) lineages, putatively predating the WGD event(s). Early duplicated accA and accB genes underwent transcriptional sub-functionalization.

The RNA-Seq-based high resolution gene expression atlas of chickpea (Cicer arietinum L.) reveals dynamic spatio-temporal changes associated with growth and development

Chickpea is one of the world's largest cultivated food legumes and is an excellent source of high-quality protein to the human diet. Plant growth and development are controlled by programmed expression of a suite of genes at the given time, stage, and tissue. Understanding how the underlying genome sequence translates into specific plant phenotypes at key developmental stages, information on gene expression patterns is crucial. Here, we present a comprehensive Cicer arietinum Gene Expression Atlas (CaGEA) across different plant developmental stages and organs covering the entire life cycle of chickpea. One of the widely used drought tolerant cultivars, ICC 4958 has been used to generate RNA-Seq data from 27 samples at 5 major developmental stages of the plant. A total of 816 million raw reads were generated and of these, 794 million filtered reads after quality control (QC) were subjected to downstream analysis. A total of 15,947 unique number of differentially expressed genes across different pairwise tissue combinations were identified. Significant differences in gene expression patterns contributing in the process of flowering, nodulation, and seed and root development were inferred in this study. Furthermore, differentially expressed candidate genes from "QTL-hotspot" region associated with drought stress response in chickpea were validated.

Systemic control of legume susceptibility to rhizobial infection by a mobile microRNA

Nitrogen-fixing root nodules on legumes result from two developmental processes, bacterial infection and nodule organogenesis. To balance symbiosis and plant growth, legume hosts restrict nodule numbers through an inducible autoregulatory process. Here, we present a mechanism where repression of a negative regulator ensures symbiotic susceptibility of uninfected roots of the host Lotus japonicus We show that microRNA miR2111 undergoes shoot-to-root translocation to control rhizobial infection through posttranscriptional regulation of the symbiosis suppressor TOO MUCH LOVE in roots. miR2111 maintains a susceptible default status in uninfected hosts and functions as an activator of symbiosis downstream of LOTUS HISTIDINE KINASE1-mediated cytokinin perception in roots and HYPERNODULATION ABERRANT ROOT FORMATION1, a shoot factor in autoregulation. The miR2111-TML node ensures activation of feedback regulation to balance infection and nodulation events.

Temporal transcriptome profiling of developing seeds reveals a concerted gene regulation in relation to oil accumulation in Pongamia (Millettia pinnata)

BACKGROUND

Pongamia (Millettia pinnata syn. Pongamia pinnata), an oilseed legume species, is emerging as potential feedstock for sustainable biodiesel production. Breeding Pongamia for favorable traits in commercial application will rely on a comprehensive understanding of molecular mechanism regulating oil accumulation during its seed development. To date, only limited genomic or transcript sequences are available for Pongamia, while a temporal transcriptome profiling of developing seeds is still lacking in this species.

RESULTS

In this work, we conducted a time-series analysis of morphological and physiological characters, oil contents and compositions, as well as global gene expression profiles in developing Pongamia seeds. Firstly, three major developmental phases were characterized based on the combined evidences from embryonic shape, seed weight, seed moisture content, and seed color. Then, the gene expression levels at these three phases were quantified by RNA-Seq analyses with three biological replicates from each phase. Nearly 94% of unigenes were expressed at all three phases, whereas only less than 2% of unigenes were exclusively expressed at one of these phases. A total of 8881 differentially expressed genes (DEGs) were identified between phases. Furthermore, the qRT-PCR analyses for 10 DEGs involved in lipid metabolism demonstrated a good reliability of our RNA-Seq data in temporal gene expression profiling. We observed a dramatic increase in seed oil content from the embryogenesis phase to the early seed-filling phase, followed by a steady and moderate increase towards the maximum at the desiccation phase. We proposed that a highly active expression of most genes related to fatty acid (FA) and triacylglycerol (TAG) biosynthesis at the embryogenesis phase might trigger both the substantial oil accumulation and the membrane lipid synthesis for rapid cell proliferation at this phase, while a concerted reactivation of TAG synthesis-related genes at the desiccation phase might further promote storage lipid synthesis to achieve the maximum content of seed oils.

CONCLUSIONS

This study not only built a bridge between gene expression profiles and oil accumulation in developing seeds, but also laid a foundation for future attempts on genetic engineering of Pongamia varieties to acquire higher oil yield or improved oil properties for biofuel applications.

Lotus japonicus NOOT-BOP-COCH-LIKE1 is essential for nodule, nectary, leaf and flower development

The NOOT-BOP-COCH-LIKE (NBCL) genes are orthologs of Arabidopsis thaliana BLADE-ON-PETIOLE1/2. The NBCLs are developmental regulators essential for plant shaping, mainly through the regulation of organ boundaries, the promotion of lateral organ differentiation and the acquisition of organ identity. In addition to their roles in leaf, stipule and flower development, NBCLs are required for maintaining the identity of indeterminate nitrogen-fixing nodules with persistent meristems in legumes. In legumes forming determinate nodules, without persistent meristem, the roles of NBCL genes are not known. We thus investigated the role of Lotus japonicus NOOT-BOP-COCH-LIKE1 (LjNBCL1) in determinate nodule identity and studied its functions in aerial organ development using LORE1 insertional mutants and RNA interference-mediated silencing approaches. In Lotus, LjNBCL1 is involved in leaf patterning and participates in the regulation of axillary outgrowth. Wild-type Lotus leaves are composed of five leaflets and possess a pair of nectaries at the leaf axil. Legumes such as pea and Medicago have a pair of stipules, rather than nectaries, at the base of their leaves. In Ljnbcl1, nectary development is abolished, demonstrating that nectaries and stipules share a common evolutionary origin. In addition, ectopic roots arising from nodule vascular meristems and reorganization of the nodule vascular bundle vessels were observed on Ljnbcl1 nodules. This demonstrates that NBCL functions are conserved in both indeterminate and determinate nodules through the maintenance of nodule vascular bundle identity. In contrast to its role in floral patterning described in other plants, LjNBCL1 appears essential for the development of both secondary inflorescence meristem and floral meristem.

Diel pattern of circadian clock and storage protein gene expression in leaves and during seed filling in cowpea (Vigna unguiculata)

BACKGROUND

Cowpea (Vigna unguiculata) is an important source of protein supply for animal and human nutrition. The major storage globulins VICILIN and LEGUMIN (LEG) are synthesized from several genes including LEGA, LEGB, LEGJ and CVC (CONVICILIN). The current hypothesis is that the plant circadian core clock genes are conserved in a wide array of species and that primary metabolism is to a large extent controlled by the plant circadian clock. Our aim was to investigate a possible link between gene expression of storage proteins and the circadian clock.

RESULTS

We identified cowpea orthologues of the core clock genes VunLHY, VunTOC1, VunGI and VunELF3, the protein storage genes VunLEG, VunLEGJ, and VunCVC as well as nine candidate reference genes used in RT-PCR. ELONGATION FACTOR 1-A (ELF1A) resulted the most suitable reference gene. The clock genes VunELF3, VunGI, VunTOC1 and VunLHY showed a rhythmic expression profile in leaves with a typical evening/night and morning/midday phased expression. The diel patterns were not completely robust and only VungGI and VungELF3 retained a rhythmic pattern under free running conditions of darkness. Under field conditions, rhythmicity and phasing apparently faded during early pod and seed development and was regained in ripening pods for VunTOC1 and VunLHY. Mature seeds showed a rhythmic expression of VunGI resembling leaf tissue under controlled growth chamber conditions. Comparing time windows during developmental stages we found that VunCVC and VunLEG were significantly down regulated during the night in mature pods as compared to intermediate ripe pods, while changes in seeds were non-significant due to high variance. The rhythmic expression under field conditions was lost under growth chamber conditions.

CONCLUSIONS

The core clock gene network is conserved in cowpea leaves showing a robust diel expression pattern except VunELF3 under growth chamber conditions. There appears to be a clock transcriptional reprogramming in pods and seeds compared to leaves. Storage protein deposition may be circadian regulated under field conditions but the strong environmental signals are not met under artificial growth conditions. Diel expression pattern in field conditions may result in better usage of energy for protein storage.

Genistein-Specific G6DT Gene for the Inducible Production of Wighteone in Lotus japonicus

Prenylated isoflavonoids have been found in several legume plants, and they possess various biological activities that play important roles in both plant defense and human health. However, it is still unknown whether prenylated isoflavonoids are present in the model legume plant Lotus japonicus. In the present study, we found that the prenylated isoflavonoid wighteone was produced in L. japonicus when leaf was supplemented with genistein. Furthermore, a novel prenyltransferase gene, LjG6DT, was identified, which shared high similarity with and was closely related to several known prenyltransferase genes involved in isoflavonoid biosynthesis. The recombinant LjG6DT protein expressed in yeast exhibited prenylation activity toward genistein as an exclusive substrate, which produced wighteone, a prenylated genistein at the C-6 position that occurs normally in legume plants. The LjG6DT-green fluorescent protein (GFP) fusion protein is targeted to plastids. The transcript level of LjG6DT is induced by glutathione, methyl jasmonate and salicylic acid, implying that LjG6DT is involved in stress response. Overexpression of LjG6DT in L. japonicus hairy roots led to increased accumulation of wighteone when genistein was supplied, indicating that LjG6DT is functional in vivo. Feeding assays with the upstream intermediate naringenin revealed that accumulation of wighteone in L. japonicus was dependent on genistein supplementation, and accumulation of wighteone is competed by genistein methylation. This study demonstrated that phytoalexin wighteone is inducibly produced in L. japonicus, and it provides new insight into the biosynthesis and accumulation of prenylated isoflavonoids in legume plants.

Genome-Wide Identification, Characterization, and Expression Profiling of the Legume BZR Transcription Factor Gene Family

The BRASSINAZOLE-RESISTANT (BZR) family of transcription factors (TFs) are positive regulators in the biosynthesis of brassinosteroids. The latter is a class of steroid hormones that affect a variety of developmental and physiological processes in plants. BZR TFs play essential roles in the regulation of plant growth and development, including multiple stress-resistance functions. However, the evolutionary history and individual expression patterns of the legume BZR genes has not been determined. In this study, we performed a genome-wide investigation of the BZR gene family in seven legume species. In total, 52 BZR genes were identified and characterized. By analyzing their phylogeny, we divided these BZR genes into five groups by comparison with orthologs/paralogs in Arabidopsis thaliana. The intron/exon structural patterns and conserved protein motifs of each gene were analyzed and showed high group-specificities. Legume BZR genes were unevenly distributed among their corresponding genomes. Genome and gene sequence comparisons revealed that gene expansion of the BZR TF family in legumes mainly resulted from segmental duplications and that this family has undergone purifying selection. Synteny analysis showed that BZR genes tended to localize within syntenic blocks conserved across legume genomes. The expression patterns of BZR genes among various legume vegetative tissues and in response to different abiotic stresses were analyzed using a combination of public transcriptome data and quantitative PCR. The patterns indicated that many BZR genes regulate legume organ development and differentiation, and significantly respond to drought and salt stresses. This study may provide valuable information for understanding the evolution of BZR gene structure and expression, and lays a foundation for future functional analysis of the legume BZR genes by species and by gene.

Distinct Lotus japonicus Transcriptomic Responses to a Spectrum of Bacteria Ranging From Symbiotic to Pathogenic

Lotus japonicus is a well-studied nodulating legume and a model organism for the investigation of plant-microbe interactions. The majority of legume transcriptome studies have focused on interactions with compatible symbionts, whereas responses to non-adapted rhizobia and pathogenic bacteria have not been well-characterized. In this study, we first characterized the transcriptomic response of L. japonicus to its compatible symbiont, Mesorhizobium loti R7A, through RNA-seq analysis of various plant tissues. Early symbiotic signaling was largely Nod factor-dependent and enhanced within root hairs, and we observed large-scale transcriptional reprogramming in nodule primordia and mature nitrogen-fixing nodules. We then characterized root transcriptional responses to a spectrum of L. japonicus interacting bacteria ranging from semi-compatible symbionts to pathogens. M. loti R7A and the semi-compatible strain Sinorhizobium fredii HH103 showed remarkably similar responses, allowing us to identify a small number of genes potentially involved in differentiating between fully and semi-compatible symbionts. The incompatible symbiont Bradyrhizobium elkanii USDA61 induced a more attenuated response, but the weakest response was observed for the foliar pathogen Pseudomonas syringae pv. tomato DC3000, where the affected genes also responded to other tested bacteria, pointing to a small set of common bacterial response genes. In contrast, the root pathogen Ralstonia solanacearum JS763 induced a pronounced and distinct transcriptomic pathogen response, which we compared to the results of the other treatments. This comparative analysis did not support the concept that an early defense-like response is generally evoked by compatible rhizobia during establishment of symbiosis.

The Medicago truncatula GRAS protein RAD1 supports arbuscular mycorrhiza symbiosis and Phytophthora palmivora susceptibility

The roots of most land plants are colonized by symbiotic arbuscular mycorrhiza (AM) fungi. To facilitate this symbiosis, plant genomes encode a set of genes required for microbial perception and accommodation. However, the extent to which infection by filamentous root pathogens also relies on some of these genes remains an open question. Here, we used genome-wide association mapping to identify genes contributing to colonization of Medicago truncatula roots by the pathogenic oomycete Phytophthora palmivora. Single-nucleotide polymorphism (SNP) markers most significantly associated with plant colonization response were identified upstream of RAD1, which encodes a GRAS transcription regulator first negatively implicated in root nodule symbiosis and recently identified as a positive regulator of AM symbiosis. RAD1 transcript levels are up-regulated both in response to AM fungus and, to a lower extent, in infected tissues by P. palmivora where its expression is restricted to root cortex cells proximal to pathogen hyphae. Reverse genetics showed that reduction of RAD1 transcript levels as well as a rad1 mutant are impaired in their full colonization by AM fungi as well as by P. palmivora. Thus, the importance of RAD1 extends beyond symbiotic interactions, suggesting a general involvement in M. truncatula microbe-induced root development and interactions with unrelated beneficial and detrimental filamentous microbes.

A single amino acid substitution in a chitinase of the legume Medicago truncatula is sufficient to gain Nod-factor hydrolase activity

The symbiotic interaction between nitrogen-fixing rhizobia and legumes depends on lipo-chitooligosaccharidic Nod-factors (NFs). The NF hydrolase MtNFH1 of Medicago truncatula is a symbiotic enzyme that hydrolytically inactivates NFs with a C16 : 2 acyl chain produced by the microsymbiont Sinorhizobium meliloti 1021. MtNFH1 is related to class V chitinases (glycoside hydrolase family 18) but lacks chitinase activity. Here, we investigated the substrate specificity of MtNFH1-related proteins. MtCHIT5a and MtCHIT5b of M. truncatula as well as LjCHIT5 of Lotus japonicus showed chitinase activity, suggesting a role in plant defence. The enzymes failed to hydrolyse NFs from S. meliloti. NFs from Rhizobium leguminosarum with a C18 : 4 acyl moiety were neither hydrolysed by these chitinases nor by MtNFH1. Construction of chimeric proteins and further amino acid replacements in MtCHIT5b were performed to identify chitinase variants that gained the ability to hydrolyse NFs. A single serine-to-proline substitution was sufficient to convert MtCHIT5b into an NF-cleaving enzyme. MtNFH1 with the corresponding proline-to-serine substitution failed to hydrolyse NFs. These results are in agreement with a substrate-enzyme model that predicts NF cleavage when the C16 : 2 moiety is placed into a distinct fatty acid-binding cleft. Our findings support the view that MtNFH1 evolved from the ancestral MtCHIT5b by gene duplication and subsequent symbiosis-related neofunctionalization.

Legume isoflavone synthase genes have evolved by whole-genome and local duplications yielding transcriptionally active paralogs

Isoflavone synthase (IFS) is the key enzyme of isoflavonoid biosynthesis. IFS genes were identified in numerous species, although their evolutionary patterns have not yet been reconstructed. To address this issue, we performed structural and functional genomic analysis. Narrow leafed lupin, Lupinus angustifolius L., was used as a reference species for the genus, because it has the most developed molecular tools available. Nuclear genome BAC library clones carrying IFS homologs were localized by linkage mapping and fluorescence in situ hybridization in three chromosome pairs. Annotation of BAC, scaffold and transcriptome sequences confirmed the presence of three full-length IFS genes in the genome. Microsynteny analysis and Bayesian inference provided clear evidence that IFS genes in legumes have evolved by lineage-specific whole-genome and tandem duplications. Gene expression profiling and RNA-seq data mining showed that the vast majority of legume IFS copies have maintained their transcriptional activity. L. angustifolius IFS homologs exhibited organ-specific expression patterns similar to those observed in other Papilionoideae. Duplicated lupin IFS homologs retained non-negligible levels of substitutions in conserved motifs, putatively due to positive selection acting during early evolution of the genus, before the whole-genome duplication. Strong purifying selection preserved newly arisen IFS duplicates from further nonsynonymous changes.

CLE peptide-encoding gene families in Medicago truncatula and Lotus japonicus, compared with those of soybean, common bean and Arabidopsis

CLE peptide hormones are critical regulators of many cell proliferation and differentiation mechanisms in plants. These 12-13 amino acid glycosylated peptides play vital roles in a diverse range of plant tissues, including the shoot, root and vasculature. CLE peptides are also involved in controlling legume nodulation. Here, the entire family of CLE peptide-encoding genes was identified in Medicago truncatula (52) and Lotus japonicus (53), including pseudogenes and non-functional sequences that were identified. An array of bioinformatic techniques were used to compare and contrast these complete CLE peptide-encoding gene families with those of fellow legumes, Glycine max and Phaseolus vulgaris, in addition to the model plant Arabidopsis thaliana. This approach provided insight into the evolution of CLE peptide families and enabled us to establish putative M. truncatula and L. japonicus orthologues. This includes orthologues of nodulation-suppressing CLE peptides and AtCLE40 that controls the stem cell population of the root apical meristem. A transcriptional meta-analysis was also conducted to help elucidate the function of the CLE peptide family members. Collectively, our analyses considerably increased the number of annotated CLE peptides in the model legume species, M. truncatula and L. japonicus, and substantially enhanced the knowledgebase of this critical class of peptide hormones.

Evolutionary Analysis of the LAFL Genes Involved in the Land Plant Seed Maturation Program

Seeds are one of the most significant innovations in the land plant lineage, critical to the diversification and adaptation of plants to terrestrial environments. From perspective of seed evo-devo, the most crucial developmental stage in this innovation is seed maturation, which includes accumulation of storage reserves, acquisition of desiccation tolerance, and induction of dormancy. Based on previous studies of seed development in the model plant Arabidopsis thaliana, seed maturation is mainly controlled by the LAFL regulatory network, which includes LEAFY COTYLEDON1 (LEC1) and LEC1-LIKE (L1L) of the NF-YB gene family, and ABSCISIC ACID INSENSITIVE3 (ABI3), FUSCA3 (FUS3), and LEC2 (LEAFY COTYLEDON2) of the B3-AFL gene family. In the present study, molecular evolution of these LAFL genes was analyzed, using representative species from across the major plant lineages. Additionally, to elucidate the molecular mechanisms of the seed maturation program, co-expression pattern analyses of LAFL genes were conducted across vascular plants. The results show that the origin of AFL gene family dates back to a common ancestor of bryophytes and vascular plants, while LEC1-type genes are only found in vascular plants. LAFL genes of vascular plants likely specify their co-expression in two different developmental phrases, spore and seed maturation, respectively, and expression patterns vary slightly across the major vascular plants lineages. All the information presented in this study will provide insights into the origin and diversification of seed plants.

Gene expression atlas of pigeonpea and its application to gain insights into genes associated with pollen fertility implicated in seed formation

Pigeonpea (Cajanus cajan) is an important grain legume of the semi-arid tropics, mainly used for its protein rich seeds. To link the genome sequence information with agronomic traits resulting from specific developmental processes, a Cajanus cajan gene expression atlas (CcGEA) was developed using the Asha genotype. Thirty tissues/organs representing developmental stages from germination to senescence were used to generate 590.84 million paired-end RNA-Seq data. The CcGEA revealed a compendium of 28 793 genes with differential, specific, spatio-temporal and constitutive expression during various stages of development in different tissues. As an example to demonstrate the application of the CcGEA, a network of 28 flower-related genes analysed for cis-regulatory elements and splicing variants has been identified. In addition, expression analysis of these candidate genes in male sterile and male fertile genotypes suggested their critical role in normal pollen development leading to seed formation. Gene network analysis also identified two regulatory genes, a pollen-specific SF3 and a sucrose-proton symporter, that could have implications for improvement of agronomic traits such as seed production and yield. In conclusion, the CcGEA provides a valuable resource for pigeonpea to identify candidate genes involved in specific developmental processes and to understand the well-orchestrated growth and developmental process in this resilient crop.

The Ethylene Responsive Factor Required for Nodulation 1 (ERN1) Transcription Factor Is Required for Infection-Thread Formation in Lotus japonicus

Several hundred genes are transcriptionally regulated during infection-thread formation and development of nitrogen-fixing root nodules. We have characterized a set of Lotus japonicus mutants impaired in root-nodule formation and found that the causative gene, Ern1, encodes a protein with a characteristic APETALA2/Ethylene Responsive Factor (AP2/ERF) transcription-factor domain. Phenotypic characterization of four ern1 alleles shows that infection pockets are formed but root-hair infection threads are absent. Formation of root-nodule primordia is delayed and no normal transcellular infection threads are found in the infected nodules. Corroborating the role of ERN1 (ERF Required for Nodulation1) in nodule organogenesis, spontaneous nodulation induced by an autoactive CCaMK and cytokinin-induced nodule primordia were not observed in ern1 mutants. Expression of Ern1 is induced in the susceptible zone by Nod factor treatment or rhizobial inoculation. At the cellular level, the pErn1:GUS reporter is highly expressed in root epidermal cells of the susceptible zone and in the cortical cells that form nodule primordia. The genetic regulation of this cellular expression pattern was further investigated in symbiotic mutants. Nod factor induction of Ern1 in epidermal cells was found to depend on Nfr1, Cyclops, and Nsp2 but was independent of Nin and Nf-ya1. These results suggest that ERN1 functions as a transcriptional regulator involved in the formation of infection threads and development of nodule primordia and may coordinate these two processes.

Transcriptome profiles reveal gene regulation of peanut (Arachis hypogaea L.) nodulation

The molecular mechanisms of symbiosis in cultivated peanut with a ‘crack entry’ infection process are largely understudied. In this study, we investigated the root transcriptional profiles of two pairs of non-nodulating (nod−) and nodulating (nod+) sister inbred peanut lines, E4/E5 and E7/E6, and their nod+ parents, F487A and PI262090 during rhizobial infection and nodule initiation by using RNA-seq technology. A total of 143, 101, 123, 215, 182, and 289 differentially expressed genes (DEGs) were identified in nod− E4, E7 and nod+ E5, E6, F487A, and PI262090 after inoculation with Bradyrhizobium sp. Different deficiencies at upstream of symbiotic signaling pathway were revealed in the two nod− genotypes. DEGs specific in nod+ genotypes included orthologs to some known symbiotic signaling pathway genes, such as NFR5, NSP2, NIN, ERN1, and many other novel and/or functionally unknown genes. Gene ontology (GO) enrichment analysis of nod+ specific DEGs revealed 54 significantly enriched GO terms, including oxidation-reduction process, metabolic process, and catalytic activity. Genes related with plant defense systems, hormone biosynthesis and response were particularly enriched. To our knowledge, this is the first report revealing symbiosis-related genes in a genome-wide manner in peanut representative of the ‘crack entry’ species.

A Combined Comparative Transcriptomic, Metabolomic, and Anatomical Analyses of Two Key Domestication Traits: Pod Dehiscence and Seed Dormancy in Pea (Pisum sp.)

The origin of the agriculture was one of the turning points in human history, and a central part of this was the evolution of new plant forms, domesticated crops. Seed dispersal and germination are two key traits which have been selected to facilitate cultivation and harvesting of crops. The objective of this study was to analyze anatomical structure of seed coat and pod, identify metabolic compounds associated with water-impermeable seed coat and differentially expressed genes involved in pea seed dormancy and pod dehiscence. Comparative anatomical, metabolomics, and transcriptomic analyses were carried out on wild dormant, dehiscent Pisum elatius (JI64, VIR320) and cultivated, indehiscent Pisum sativum non-dormant (JI92, Cameor) and recombinant inbred lines (RILs). Considerable differences were found in texture of testa surface, length of macrosclereids, and seed coat thickness. Histochemical and biochemical analyses indicated genotype related variation in composition and heterogeneity of seed coat cell walls within macrosclereids. Liquid chromatography-electrospray ionization/mass spectrometry and Laser desorption/ionization-mass spectrometry of separated seed coats revealed significantly higher contents of proanthocyanidins (dimer and trimer of gallocatechin), quercetin, and myricetin rhamnosides and hydroxylated fatty acids in dormant compared to non-dormant genotypes. Bulk Segregant Analysis coupled to high throughput RNA sequencing resulted in identification of 770 and 148 differentially expressed genes between dormant and non-dormant seeds or dehiscent and indehiscent pods, respectively. The expression of 14 selected dormancy-related genes was studied by qRT-PCR. Of these, expression pattern of four genes: porin (MACE-S082), peroxisomal membrane PEX14-like protein (MACE-S108), 4-coumarate CoA ligase (MACE-S131), and UDP-glucosyl transferase (MACE-S139) was in agreement in all four genotypes with Massive analysis of cDNA Ends (MACE) data. In case of pod dehiscence, the analysis of two candidate genes (SHATTERING and SHATTERPROOF) and three out of 20 MACE identified genes (MACE-P004, MACE-P013, MACE-P015) showed down-expression in dorsal and ventral pod suture of indehiscent genotypes. Moreover, MACE-P015, the homolog of peptidoglycan-binding domain or proline-rich extensin-like protein mapped correctly to predicted Dpo1 locus on PsLGIII. This integrated analysis of the seed coat in wild and cultivated pea provides new insight as well as raises new questions associated with domestication and seed dormancy and pod dehiscence.

Involvement of three putative glucosyltransferases from the UGT72 family in flavonol glucoside/rhamnoside biosynthesis in Lotus japonicus seeds

Flavonols are one of the largest groups of flavonoids that confer benefits for the health of plants and animals. Flavonol glycosides are the predominant flavonoids present in the model legume Lotus japonicus. The molecular mechanisms underlying the biosynthesis of flavonol glycosides as yet remain unknown in L. japonicus. In the present study, we identified a total of 188 UDP-glycosyltransferases (UGTs) in L. japonicus by genome-wide searching. Notably, 12 UGTs from the UGT72 family were distributed widely among L. japonicus chromosomes, expressed in all tissues, and showed different docking scores in an in silico bioinformatics docking analysis. Further enzymatic assays showed that five recombinant UGTs (UGT72AD1, UGT72AF1, UGT72AH1, UGT72V3, and UGT72Z2) exhibit activity toward flavonol, flavone, and isoflavone aglycones. In particular, UGT72AD1, UGT72AH1, and UGT72Z2 are flavonol-specific UGTs with different kinetic properties. In addition, the overexpression of UGT72AD1 and UGT72Z2 led to increased accumulation of flavonol rhamnosides in L. japonicus and Arabidopsis thaliana. Moreover, the increase of kaempferol 3-O-rhamnoside-7-O-rhamnoside in transgenic A. thaliana inhibited root growth as compared with the wild-type control. These results highlight the significance of the UGT72 family in flavonol glycosylation and the role of flavonol rhamnosides in plant growth.

Lotus Base: An integrated information portal for the model legume Lotus japonicus

Lotus japonicus is a well-characterized model legume widely used in the study of plant-microbe interactions. However, datasets from various Lotus studies are poorly integrated and lack interoperability. We recognize the need for a comprehensive repository that allows comprehensive and dynamic exploration of Lotus genomic and transcriptomic data. Equally important are user-friendly in-browser tools designed for data visualization and interpretation. Here, we present Lotus Base, which opens to the research community a large, established LORE1 insertion mutant population containing an excess of 120,000 lines, and serves the end-user tightly integrated data from Lotus, such as the reference genome, annotated proteins, and expression profiling data. We report the integration of expression data from the L. japonicus gene expression atlas project, and the development of tools to cluster and export such data, allowing users to construct, visualize, and annotate co-expression gene networks. Lotus Base takes advantage of modern advances in browser technology to deliver powerful data interpretation for biologists. Its modular construction and publicly available application programming interface enable developers to tap into the wealth of integrated Lotus data. Lotus Base is freely accessible at: https://lotus.au.dk.

Lotus japonicus NF-YA1 Plays an Essential Role During Nodule Differentiation and Targets Members of the SHI/STY Gene Family

Legume plants engage in intimate relationships with rhizobial bacteria to form nitrogen-fixing nodules, root-derived organs that accommodate the microsymbiont. Members of the Nuclear Factor Y (NF-Y) gene family, which have undergone significant expansion and functional diversification during plant evolution, are essential for this symbiotic liaison. Acting in a partially redundant manner, NF-Y proteins were shown, previously, to regulate bacterial infection, including selection of a superior rhizobial strain, and to mediate nodule structure formation. However, the exact mechanism by which these transcriptional factors exert their symbiotic functions has remained elusive. By carrying out detailed functional analyses of Lotus japonicus mutants, we demonstrate that LjNF-YA1 becomes indispensable downstream from the initial cortical cell divisions but prior to nodule differentiation, including cell enlargement and vascular bundle formation. Three affiliates of the SHORT INTERNODES/STYLISH transcription factor gene family, called STY1, STY2, and STY3, are demonstrated to be among likely direct targets of LjNF-YA1, and our results point to their involvement in nodule formation.

The relationship between thiamine and two symbioses: Root nodule symbiosis and arbuscular mycorrhiza

Lotus japonicus THIC is expressed in all organs, and the encoded protein catalyzes thiamine biosynthesis. Loss of function produces chlorosis, a typical thiamine-deficiency phenotype, and mortality. To investigate thiamine's role in symbiosis, we focused on THI1, a thiamine-biosynthesis gene expressed in roots, nodules, and seeds. The thi1 mutant had green leaves, but formed small nodules and immature seeds. These phenotypes were rescued by THI1 complementation and by exogenous thiamine. Thus, THI1 is required for nodule enlargement and seed maturation. On the other hand, colonization by arbuscular mycorrhiza (AM) fungus Rhizophagus irregularis was not affected in the thi1 mutant or by exogenous thiamine. However, spores of R. irregularis stored more thiamine than the source (host plants), despite lacking thiamine biosynthesis genes. Therefore, disturbance of the thiamine supply would affect progeny phenotypes such as spore formation and hyphal growth. Further investigation will be required to elucidate thiamine's effect on AM.

Generation and Characterisation of a Reference Transcriptome for Lentil (Lens culinaris Medik.)

RNA-Seq using second-generation sequencing technologies permits generation of a reference unigene set for a given species, in the absence of a well-annotated genome sequence, supporting functional genomics studies, gene characterisation and detailed expression analysis for specific morphophysiological or environmental stress response traits. A reference unigene set for lentil has been developed, consisting of 58,986 contigs and scaffolds with an N50 length of 1719 bp. Comparison to gene complements from related species, reference protein databases, previously published lentil transcriptomes and a draft genome sequence validated the current dataset in terms of degree of completeness and utility. A large proportion (98%) of unigenes were expressed in more than one tissue, at varying levels. Candidate genes associated with mechanisms of tolerance to both boron toxicity and time of flowering were identified, which can eventually be used for the development of gene-based markers. This study has provided a comprehensive, assembled and annotated reference gene set for lentil that can be used for multiple applications, permitting identification of genes for pathway-specific expression analysis, genetic modification approaches, development of resources for genotypic analysis, and assistance in the annotation of a future lentil genome sequence.

The Vigna unguiculata Gene Expression Atlas (VuGEA) from de novo assembly and quantification of RNA-seq data provides insights into seed maturation mechanisms

Legume research and cultivar development are important for sustainable food production, especially of high-protein seed. Thanks to the development of deep-sequencing technologies, crop species have been taken to the front line, even without completion of their genome sequences. Black-eyed pea (Vigna unguiculata) is a legume species widely grown in semi-arid regions, which has high potential to provide stable seed protein production in a broad range of environments, including drought conditions. The black-eyed pea reference genotype has been used to generate a gene expression atlas of the major plant tissues (i.e. leaf, root, stem, flower, pod and seed), with a developmental time series for pods and seeds. From these various organs, 27 cDNA libraries were generated and sequenced, resulting in more than one billion reads. Following filtering, these reads were de novo assembled into 36 529 transcript sequences that were annotated and quantified across the different tissues. A set of 24 866 unique transcript sequences, called Unigenes, was identified. All the information related to transcript identification, annotation and quantification were stored into a gene expression atlas webserver (http://vugea.noble.org), providing a user-friendly interface and necessary tools to analyse transcript expression in black-eyed pea organs and to compare data with other legume species. Using this gene expression atlas, we inferred details of molecular processes that are active during seed development, and identified key putative regulators of seed maturation. Additionally, we found evidence for conservation of regulatory mechanisms involving miRNA in plant tissues subjected to drought and seeds undergoing desiccation.

Global analysis of the developmental dynamics of Gossypium hirsutum based on strand-specific transcriptome

Cotton is an economically important crop that provides both natural fiber and by-products such as oil and protein. Its global gene expression could provide insight into the biological processes underlying growth and development, which involve suites of genes expressed with temporal and spatial control by regulatory networks. Generally, the goal for cotton breeding is improvement of the fiber; thus, most previous research has focused on identifying genes specific to the fiber. However, seeds may also play an important role in fiber development. In this study, we constructed and systematically analyzed 21 strand-specific RNA-Seq libraries for Gossypium hirsutum, covering different tissues, organs and development stages, from which approximately 970 million reads were generated to provide a global view of gene expression during cotton development. The organ (tissue)-specific gene expression patterns were investigated, providing further insight into the dynamic programming associated with developmental processes and a way to study the coordination of development between fiber cells and ovules. Series of transcription factors and seed-specific genes have been identified as candidate genes that could elucidate key mechanisms and regulatory networks in nutrient accumulation during ovule development and in fiber development. This study reports comprehensive transcriptome dynamics at various stages of cotton development and will serve as a valuable genome-wide transcriptome resource for initial gene discovery and functional characterization of genes in cotton.

Genome-Wide Identification, Evolutionary Analysis and Expression Profiles of LATERAL ORGAN BOUNDARIES DOMAIN Gene Family in Lotus japonicus and Medicago truncatula

The LATERAL ORGANBOUNDARIESDOMAIN (LBD) gene family has been well-studied in Arabidopsis and play crucial roles in the diverse growth and development processes including establishment and maintenance of boundary of developmental lateral organs. In this study we identified and characterized 38 LBD genes in Lotus japonicus (LjLBD) and 57 LBD genes in Medicago truncatula (MtLBD), both of which are model legume plants that have some specific development features absent in Arabidopsis. The phylogenetic relationships, their locations in the genome, genes structure and conserved motifs were examined. The results revealed that all LjLBD and MtLBD genes could be distinctly divided into two classes: Class I and II. The evolutionary analysis showed that Type I functional divergence with some significantly site-specific shifts may be the main force for the divergence between Class I and Class II. In addition, the expression patterns of LjLBD genes uncovered the diverse functions in plant development. Interestingly, we found that two LjLBD proteins that were highly expressed during compound leaf and pulvinus development, can interact via yeast two-hybrid assays. Taken together, our findings provide an evolutionary and genetic foundation in further understanding the molecular basis of LBD gene family in general, specifically in L. japonicus and M. truncatula.

An Expanded Maize Gene Expression Atlas based on RNA Sequencing and its Use to Explore Root Development

Comprehensive and systematic transcriptome profiling provides valuable insight into biological and developmental processes that occur throughout the life cycle of a plant. We have enhanced our previously published microarray-based gene atlas of maize ( L.) inbred B73 to now include 79 distinct replicated samples that have been interrogated using RNA sequencing (RNA-seq). The current version of the atlas includes 50 original array-based gene atlas samples, a time-course of 12 stalk and leaf samples postflowering, and an additional set of 17 samples from the maize seedling and adult root system. The entire dataset contains 4.6 billion mapped reads, with an average of 20.5 million mapped reads per biological replicate, allowing for detection of genes with lower transcript abundance. As the new root samples represent key additions to the previously examined tissues, we highlight insights into the root transcriptome, which is represented by 28,894 (73.2%) annotated genes in maize. Additionally, we observed remarkable expression differences across both the longitudinal (four zones) and radial gradients (cortical parenchyma and stele) of the primary root supported by fourfold differential expression of 9353 and 4728 genes, respectively. Among the latter were 1110 genes that encode transcription factors, some of which are orthologs of previously characterized transcription factors known to regulate root development in (L.) Heynh., while most are novel, and represent attractive targets for reverse genetics approaches to determine their roles in this important organ. This comprehensive transcriptome dataset is a powerful tool toward understanding maize development, physiology, and phenotypic diversity.

Genome and transcriptome analysis of the Mesoamerican common bean and the role of gene duplications in establishing tissue and temporal specialization of genes

BACKGROUND

Legumes are the third largest family of angiosperms and the second most important crop class. Legume genomes have been shaped by extensive large-scale gene duplications, including an approximately 58 million year old whole genome duplication shared by most crop legumes.

RESULTS

We report the genome and the transcription atlas of coding and non-coding genes of a Mesoamerican genotype of common bean (Phaseolus vulgaris L., BAT93). Using a comprehensive phylogenomics analysis, we assessed the past and recent evolution of common bean, and traced the diversification of patterns of gene expression following duplication. We find that successive rounds of gene duplications in legumes have shaped tissue and developmental expression, leading to increased levels of specialization in larger gene families. We also find that many long non-coding RNAs are preferentially expressed in germ-line-related tissues (pods and seeds), suggesting that they play a significant role in fruit development. Our results also suggest that most bean-specific gene family expansions, including resistance gene clusters, predate the split of the Mesoamerican and Andean gene pools.

CONCLUSIONS

The genome and transcriptome data herein generated for a Mesoamerican genotype represent a counterpart to the genomic resources already available for the Andean gene pool. Altogether, this information will allow the genetic dissection of the characters involved in the domestication and adaptation of the crop, and their further implementation in breeding strategies for this important crop.