Dissecting the Root Nodule Transcriptome of Chickpea (Cicer arietinum L.)

The developmental dynamics in cool season legumes with focus on chickpea

Chickpea is one of the most widely consumed grain legume world-wide. Advances in next-generation sequencing and genomics tools have led to genetic dissection and identification of potential candidate genes regulating agronomic traits in chickpea. However, the developmental particularities and its potential in reforming the yield and nutritional value remain largely unexplored. Studies in crops such as rice, maize, tomato and pea have highlighted the contribution of key regulator of developmental events in yield related traits. A comprehensive knowledge on the development aspects of a crop can pave way for new vistas to explore. Pea and Medicago are the close relatives of genus Cicer and the basic developmental events in these legumes are similar. However, there are some distinct developmental features in chickpea which hold potential for future crop improvement endeavours. The global chickpea germplasm encompasses wide range of diversities in terms of morphology at both vegetative and reproductive stages. There is an immediate need for understanding the genetic and molecular basis of this diversity and utilizing them for the yield contributing trait improvement. The review discusses some of the key developmental events which have potential in yield enhancement and the lessons which can be learnt from model legumes in this regard.

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.

Defining the mutation sites in chickpea nodulation mutants PM233 and PM405

BACKGROUND

Like most legumes, chickpeas form specialized organs called root nodules. These nodules allow for a symbiotic relationship with rhizobium bacteria. The rhizobia provide fixed atmospheric nitrogen to the plant in a usable form. It is of both basic and practical interest to understand the host plant genetics of legume root nodulation. Chickpea lines PM233 and PM405, which harbor the mutationally identified nodulation genes rn1 and rn4, respectively, both display nodulation-deficient phenotypes. Previous investigators identified the rn1 mutation with the chickpea homolog of Medicago truncatula nodulation gene NSP2, but were unable to define the mutant rn1 allele. We used Illumina and Nanopore sequencing reads to attempt to identify and characterize candidate mutation sites responsible for the PM233 and PM405 phenotypes.

RESULTS

We aligned Illumina reads to the available desi chickpea reference genome, and did a de novo contig assembly of Nanopore reads. In mutant PM233, the Nanopore contigs allowed us to identify the breakpoints of a ~ 35 kb deleted region containing the CaNSP2 gene, the Medicago truncatula homolog of which is involved in nodulation. In mutant PM405, we performed variant calling in read alignments and identified 10 candidate mutations. Genotyping of a segregating progeny population narrowed that pool down to a single candidate gene which displayed homology to M. truncatula nodulation gene NIN.

CONCLUSIONS

We have characterized the nodulation mutation sites in chickpea mutants PM233 and PM405. In mutant PM233, the rn1 mutation was shown to be due to deletion of the entire CaNSP2 nodulation gene, while in mutant PM405 the rn4 mutation was due to a single base deletion resulting in a frameshift mutation between the predicted RWP-RK and PB1 domains of the NIN nodulation gene. Critical to characterization of the rn1 allele was the generation of Nanopore contigs for mutant PM233 and its wild type parent ICC 640, without which the deletional boundaries could not be defined. Our results suggest that efforts of prior investigators were hampered by genomic misassemblies in the CaNSP2 region of both the desi and kabuli reference genomes.

Evolutionary and functional analysis of two-component system in chickpea reveals CaRR13, a TypeB RR, as positive regulator of symbiosis

The critical role of cytokinin in early nodulation in legumes is well known. In our study, exogenous cytokinin application to roots of the important crop legume, chickpea (Cicer arietinum L.), led to the formation of pseudo-nodules even in the absence of rhizobia. Hence, a genome-wide analysis of the cytokinin signalling, two-component system (TCS) genes, was conducted in chickpea, Medicago and Cajanus cajan. The integrated phylogenetic, evolutionary and expression analysis of the TCS genes was carried out, which revealed that histidine kinases (HKs) were highly conserved, whereas there was diversification leading to neofunctionalization at the level of response regulators (RRs) especially the TypeB RRs. Further, the functional role of the CaHKs in nodulation was established by complementation of the sln1Δ mutant of yeast and cre1 mutants of (Medicago) which led to restoration of the nodule-deficient phenotype. Additionally, the highest expressing TypeB RR of chickpea, CaRR13, was functionally characterized. Its localization in the nucleus and its Y1H assay-based interaction with the promoter of the early nodulation gene CaNSP2 indicated its role as a transcription factor regulating early nodulation. Overexpression, RNAi lines and complementation of cre1 mutants with CaRR13 revealed its critical involvement as an important signalling molecule regulating early events of nodule organogenesis in chickpea.

Evolutionary and expression dynamics of LRR-RLKs and functional establishment of KLAVIER homolog in shoot mediated regulation of AON in chickpea symbiosis

Chickpea shoot exogenously treated with cytokinin showed stunted phenotype of root, shoot and significantly reduced nodule numbers. Genome-wide identification of LRR-RLKs in chickpea and Medicago resulted in 200 and 371 genes respectively. Gene duplication analysis revealed that LRR-RLKs family expanded through segmental duplications in chickpea and tandem duplications in Medicago. Expression profiling of LRR-RLKs revealed their involvement in cytokinin signaling and plant organ development. Overexpression of KLAVIER ortholog of chickpea, Ca_LRR-RLK147, in roots revealed its localization in the membrane but showed no effect on root nodulation despite increased cle peptide levels. Two findings (i) drastic effect on nodule number by exogenous cytokinin treatment to only shoot and restoration to normal nodulation by treatment to both root and shoot tissue and (ii) no effect on nodule number by overexpression of Ca_LRR-RLK147 establishes the fact that despite presence of cle peptides in root, the function of Ca_LRR-RLK147 was shoot mediated during AON.

Fulvic acid increases forage legume growth inducing preferential up-regulation of nodulation and signalling-related genes

The use of potential biostimulants is of broad interest in plant science for improving yields. The application of a humic derivative called fulvic acid (FA) may improve forage crop production. FA is an uncharacterized mixture of chemicals and, although it has been reported to increase growth parameters in many species including legumes, its mode of action remains unclear. Previous studies of the action of FA have lacked appropriate controls, and few have included field trials. Here we report yield increases due to FA application in three European Medicago sativa cultivars, in studies which include the appropriate nutritional controls which hitherto have not been used. No significant growth stimulation was seen after FA treatment in grass species in this study at the treatment rate tested. Direct application to bacteria increased Rhizobium growth and, in M. sativa trials, root nodulation was stimulated. RNA transcriptional analysis of FA-treated plants revealed up-regulation of many important early nodulation signalling genes after only 3 d. Experiments in plate, glasshouse, and field environments showed yield increases, providing substantial evidence for the use of FA to benefit M. sativa forage production.

Gene Expression in Nitrogen-Fixing Symbiotic Nodule Cells in Medicago truncatula and Other Nodulating Plants

Root nodules formed by plants of the nitrogen-fixing clade (NFC) are symbiotic organs that function in the maintenance and metabolic integration of large populations of nitrogen-fixing bacteria. These organs feature unique characteristics and processes, including their tissue organization, the presence of specific infection structures called infection threads, endocytotic uptake of bacteria, symbiotic cells carrying thousands of intracellular bacteria without signs of immune responses, and the integration of symbiont and host metabolism. The early stages of nodulation are governed by a few well-defined functions, which together constitute the common symbiosis-signaling pathway (CSSP). The CSSP activates a set of transcription factors (TFs) that orchestrate nodule organogenesis and infection. The later stages of nodule development require the activation of hundreds to thousands of genes, mostly expressed in symbiotic cells. Many of these genes are only active in symbiotic cells, reflecting the unique nature of nodules as plant structures. Although how the nodule-specific transcriptome is activated and connected to early CSSP-signaling is poorly understood, candidate TFs have been identified using transcriptomic approaches, and the importance of epigenetic and chromatin-based regulation has been demonstrated. We discuss how gene regulation analyses have advanced our understanding of nodule organogenesis, the functioning of symbiotic cells, and the evolution of symbiosis in the NFC.

A comprehensive analysis of the B3 superfamily identifies tissue-specific and stress-responsive genes in chickpea (Cicer arietinum L.)

The aim of this study was to provide a comprehensive analysis of the plant-specific B3 domain-containing transcription factors (TFs) in chickpea. Scanning of the chickpea genome resulted in the identification of 51 B3 domain-containing TFs that were located on seven out of eight chickpea chromosomes. Based on the presence of additional domains other than the B3 domain, the candidates were classified into four subfamilies, i.e., ARF (24), REM (19), LAV (6) and RAV (2). Phylogenetic analysis classified them into four groups in which members of the same group had similar intron-exon organization and motif composition. Genome duplication analysis of the candidate B3 genes revealed an event of segmental duplication that was instrumental in the expansion of the B3 gene family. Ka/Ks analysis showed that the B3 gene family was under purifying selection. Further, chickpea B3 genes showed maximum orthology with Medicago followed by soybean and Arabidopsis. Promoter analyses of the B3 genes led to the identification of several tissue-specific and stress-responsive cis-regulatory elements. Expression profiling of the candidate B3 genes using publicly available RNA-seq data of several chickpea tissues indicated their putative role in plant development and abiotic stress response. These findings were further validated by real-time expression analysis. Overall, this study provides a comprehensive analysis of the B3 domain-containing proteins in chickpea that would aid in devising strategies for crop manipulation in chickpea.

Genome-wide identification and evolutionary analysis of TGA transcription factors in soybean

The gain of function in genes and gene families is a continuous process and is a key factor in understanding gene and genome evolution in plants. TGACG-Binding (TGA) transcription factors (TFs) have long been known for their essential roles in plant defence in Arabidopsis, but their roles in legume symbiosis are yet to be explored. Here, we identified a total of 25 TGA (named GmTGA1-GmTGA25) genes in soybean. Through phylogenetic analysis, we discovered a clade of GmTGA proteins that appear to be legume-specific. Among them, two GmTGAs were unique by possessing the autophagy sequence in their proteins, while the third one was an orphan gene in soybean. GmTGAs were structurally different from AtTGAs, and their expression patterns also differed with the dominant expression of AtTGAs and GmTGAs in aerial and underground parts, respectively. Moreover, twenty-five GmTGAs showed a strong correlation among the gene expression in roots, nodules, and root hairs. The qRT-PCR analysis results revealed that among 15 tested GmTGAs, six were induced and four were suppressed by rhizobia inoculation, while 11 of these GmTGAs were induced by high nitrate. Our findings suggested the important roles of GmTGAs in symbiotic nodulation and in response to nitrogen availability in soybean.

A Toolbox for Nodule Development Studies in Chickpea: A Hairy-Root Transformation Protocol and an Efficient Laboratory Strain of Mesorhizobium sp

A Mesorhizobium sp. produces root nodules in chickpea. Chickpea and model legume Medicago truncatula are members of the inverted repeat-lacking clade (IRLC). The rhizobia, after internalization into the plant cell, are called bacteroids. Nodule-specific cysteine-rich peptides in IRLC legumes guide bacteroids to a terminally differentiated swollen (TDS) form. Bacteroids in chickpea are less TDS than those in Medicago spp. Nodule development in chickpea indicates recent evolutionary diversification and merits further study. A hairy-root transformation protocol and an efficient laboratory strain are prerequisites for performing any genetic study on nodulation. We have standardized a protocol for composite plant generation in chickpea with a transformation frequency above 50%, as shown by fluorescent markers. This protocol also works well in different ecotypes of chickpea. Localization of subcellular markers in these transformed roots is similar to the localization observed in transformed Medicago roots. When checked inside transformed nodules, peroxisomes were concentrated along the periphery of the nodules, while endoplasmic reticulum and Golgi bodies surrounded the symbiosomes. Different Mesorhizobium strains were evaluated for their ability to initiate nodule development and efficiency of nitrogen fixation. Inoculation with different strains resulted in different shapes of TDS bacteroids with variable nitrogen fixation. Our study provides a toolbox to study nodule development in the crop legume chickpea.

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.

Transcriptional Reprogramming of Legume Genomes: Perspective and Challenges Associated With Single-Cell and Single Cell-Type Approaches During Nodule Development

Transcriptomic approaches revealed thousands of genes differentially or specifically expressed during nodulation, a biological process resulting from the symbiosis between leguminous plant roots and rhizobia, atmospheric nitrogen-fixing symbiotic bacteria. Ultimately, nodulation will lead to the development of a new root organ, the nodule. Through functional genomic studies, plant transcriptomes have been used by scientists to reveal plant genes potentially controlling nodulation. However, it is important to acknowledge that the physiology, transcriptomic programs, and biochemical properties of the plant cells involved in nodulation are continuously regulated. They also differ between the different cell-types composing the nodules. To generate a more accurate picture of the transcriptome, epigenome, proteome, and metabolome of the cells infected by rhizobia and cells composing the nodule, there is a need to implement plant single-cell and single cell-types strategies and methods. Accessing such information would allow a better understanding of the infection of plant cells by rhizobia and will help understanding the complex interactions existing between rhizobia and the plant cells. In this mini-review, we are reporting the current knowledge on legume nodulation gained by plant scientists at the level of single cell-types, and provide perspectives on single cell/single cell-type approaches when applied to legume nodulation.

Current advances in chickpea genomics: applications and future perspectives

Chickpea genomics promises to illuminate our understanding of genome organization, structural variations, evolutionary and domestication-related insights and fundamental biology of legume crops. Unprecedented advancements of next generation sequencing (NGS) technologies have enabled in decoding of multiple chickpea genome sequences and generating huge genomic resources in chickpea both at functional and structural level. This review is aimed to update the current progress of chickpea genomics ranging from high density linkage map development, genome-wide association studies (GWAS), functional genomics resources for various traits, emerging role of abiotic stress responsive coding and non-coding RNAs after the completion of draft chickpea genome sequences. Additionally, the current efforts of whole genome re-sequencing (WGRS) approach of global chickpea germplasm to capture the global genetic diversity existing in the historically released varieties across the world and increasing the resolution of the previously identified candidate gene(s) of breeding importance have been discussed. Thus, the outcomes of these genomics resources will assist in genomics-assisted selection and facilitate breeding of climate-resilient chickpea cultivars for sustainable agriculture.

Morphotype of bacteroids in different legumes correlates with the number and type of symbiotic NCR peptides

In legume nodules, rhizobia differentiate into nitrogen-fixing forms called bacteroids, which are enclosed by a plant membrane in an organelle-like structure called the symbiosome. In the Inverted Repeat-Lacking Clade (IRLC) of legumes, this differentiation is terminal due to irreversible loss of cell division ability and is associated with genome amplification and different morphologies of the bacteroids that can be swollen, elongated, spherical, and elongated-branched, depending on the host plant. In Medicago truncatula, this process is orchestrated by nodule-specific cysteine-rich peptides (NCRs) delivered into developing bacteroids. Here, we identified the predicted NCR proteins in 10 legumes representing different subclades of the IRLC with distinct bacteroid morphotypes. Analysis of their expression and predicted sequences establishes correlations between the composition of the NCR family and the morphotypes of bacteroids. Although NCRs have a single origin, their evolution has followed different routes in individual lineages, and enrichment and diversification of cationic peptides has resulted in the ability to impose major morphological changes on the endosymbionts. The wide range of effects provoked by NCRs such as cell enlargement, membrane alterations and permeabilization, and biofilm and vesicle formation is dependent on the amino acid composition and charge of the peptides. These effects are strongly influenced by the rhizobial surface polysaccharides that affect NCR-induced differentiation and survival of rhizobia in nodule cells.