Global transcriptome analysis of developing chickpea (Cicer arietinum L.) seeds

MicroRNA164e suppresses NAC100 transcription factor-mediated synthesis of seed storage proteins in chickpea

Development of protein-enriched chickpea varieties necessitates an understanding of specific genes and key regulatory circuits that govern the synthesis of seed storage proteins (SSPs). Here, we demonstrated the novel involvement of Ca-miR164e-CaNAC100 in regulating SSP synthesis in chickpea. Ca-miRNA164e was significantly decreased during seed maturation, especially in high-protein accessions. The miRNA was found to directly target the transactivation conferring C-terminal region of a nuclear-localized transcription factor, CaNAC100 as revealed using RNA ligase-mediated-rapid amplification of cDNA ends and target mimic assays. The functional role of CaNAC100 was demonstrated through seed-specific overexpression (NACOE) resulting in significantly augmented seed protein content (SPC) consequential to increased SSP transcription. Further, NACOE lines displayed conspicuously enhanced seed weight but reduced numbers and yield. Conversely, a downregulation of CaNAC100 and SSP transcripts was evident in seed-specific overexpression lines of Ca-miR164e that culminated in significantly lowered SPC. CaNAC100 was additionally demonstrated to transactivate the SSP-encoding genes by directly binding to their promoters as demonstrated using electrophoretic mobility shift and dual-luciferase reporter assays. Taken together, our study for the first time established a distinct role of CaNAC100 in positively influencing SSP synthesis and its critical regulation by CamiR164e, thereby serving as an understanding that can be utilized for developing SPC-rich chickpea varieties.

De novo transcriptomic analysis of Doum Palm (Hyphaene compressa) revealed an insight into its potential drought tolerance

BACKGROUND

Doum palms (Hyphaene compressa) perform a crucial starring role in the lives of Kenya's arid and semi-arid people for empowerment and sustenance. Despite the crop's potential for economic gain, there is a lack of genetic resources and detailed information about its domestication at the molecular level. Given the doum palm's vast potential as a widely distributed plant in semi-arid and arid climates and a source of many applications, coupled with the current changing climate scenario, it is essential to understand the molecular processes that provide drought resistance to this plant.

RESULTS

Assembly of the first transcriptome of doum palms subjected to water stress generated about 39.97 Gb of RNA-Seq data. The assembled transcriptome revealed 193,167 unigenes with an average length of 1655 bp, with 128,708 (66.63%) successfully annotated in seven public databases. Unigenes exhibited significant differentially expressed genes (DEGs) in well-watered and stressed-treated plants, with 45071 and 42457 accounting for up-regulated and down-regulated DEGs, respectively. GO term, KEGG, and KOG analysis showed that DEGs were functionally enriched cellular processes, metabolic processes, cellular and catalytic activity, metabolism, genetic information processing, signal transduction mechanisms, and posttranslational modification pathways. Transcription factors (TF), such as the MYB, WRKY, NAC family, FAR1, B3, bHLH, and bZIP, were the prominent TF families identified as doum palm DEGs encoding drought stress tolerance.

CONCLUSIONS

This study provides a complete understanding of DEGs involved in drought stress at the transcriptome level in doum palms. This research is, therefore, the foundation for the characterization of potential genes, leading to a clear understanding of its drought stress responses and providing resources for improved genetic modification.

Transcriptome analysis of ovules offers early developmental clues after fertilization in Cicer arietinum L

Chickpea (Cicer arietinum L.) seeds are valued for their nutritional scores and limited information on the molecular mechanisms of chickpea fertilization and seed development is available. In the current work, comparative transcriptome analysis was performed on two different stages of chickpea ovules (pre- and post-fertilization) to identify key regulatory transcripts. Two-staged transcriptome sequencing was generated and over 208 million reads were mapped to quantify transcript abundance during fertilization events. Mapping to the reference genome showed that the majority (92.88%) of high-quality Illumina reads were aligned to the chickpea genome. Reference-guided genome and transcriptome assembly yielded a total of 28,783 genes. Of these, 3399 genes were differentially expressed after the fertilization event. These involve upregulated genes including a protease-like secreted in CO(2) response (LOC101500970), amino acid permease 4-like (LOC101506539), and downregulated genes MYB-related protein 305-like (LOC101493897), receptor like protein 29 (LOC101491695). WGCNA analysis and pairwise comparison of datasets, successfully constructed four co-expression modules. Transcription factor families including bHLH, MYB, MYB-related, C2H2 zinc finger, ERF, WRKY and NAC transcription factor were also found to be activated after fertilization. Activation of these genes and transcription factors results in the accumulation of carbohydrates and proteins by enhancing their trafficking and biosynthesis. Total 17 differentially expressed genes, were randomly selected for qRT-PCR for validation of transcriptome analysis and showed statistically significant correlations with the transcriptome data. Our findings provide insights into the regulatory mechanisms underlying changes in fertilized chickpea ovules. This work may come closer to a comprehensive understanding of the mechanisms that initiate developmental events in chickpea seeds after fertilization.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-023-03599-8.

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.

Comparative transcriptomic and metabolite profiling reveals genotype-specific responses to Fe starvation in chickpea

Iron deficiency is a major nutritional stress that severely impacts crop productivity worldwide. However, molecular intricacies and subsequent physiological and metabolic changes in response to Fe starvation, especially in leguminous crops like chickpea, remain elusive. In the present study, we investigated physiological, transcriptional, and metabolic reprogramming in two chickpea genotypes (H6013 and L4958) with contrasting seed iron concentrations upon Fe deficiency. Our findings revealed that iron starvation affected growth and physiological parameters of both chickpea genotypes. Comparative transcriptome analysis led to the identification of differentially expressed genes between the genotypes related to strategy I uptake, metal ions transporters, reactive oxygen species-associated genes, transcription factors, and protein kinases that could mitigate Fe deficiency. Our gene correlation network discovered several putative candidate genes like CIPK25, CKX3, WRKY50, NAC29, MYB4, and PAP18, which could facilitate the investigation of the molecular rationale underlying Fe tolerance in chickpea. Furthermore, the metabolite analysis also illustrated the differential accumulation of organic acids, amino acids and other metabolites associated with Fe mobilization in chickpea genotypes. Overall, our study demonstrated the comparative transcriptional dynamics upon Fe starvation. The outcomes of the current endeavor will enable the development of Fe deficiency tolerant chickpea cultivars.

Molecular Breeding and Drought Tolerance in Chickpea

Cicer arietinum L. is the third greatest widely planted imperative pulse crop worldwide, and it belongs to the Leguminosae family. Drought is the utmost common abiotic factor on plants, distressing their water status and limiting their growth and development. Chickpea genotypes have the natural ability to fight drought stress using certain strategies viz., escape, avoidance and tolerance. Assorted breeding methods, including hybridization, mutation, and marker-aided breeding, genome sequencing along with omics approaches, could be used to improve the chickpea germplasm lines(s) against drought stress. Root features, for instance depth and root biomass, have been recognized as the greatest beneficial morphological factors for managing terminal drought tolerance in the chickpea. Marker-aided selection, for example, is a genomics-assisted breeding (GAB) strategy that can considerably increase crop breeding accuracy and competence. These breeding technologies, notably marker-assisted breeding, omics, and plant physiology knowledge, underlined the importance of chickpea breeding and can be used in future crop improvement programmes to generate drought-tolerant cultivars(s).

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.

Transcriptional control of Arabidopsis seed development

The entire process of embryo development is under the tight control of various transcription factors. Together with other proteins, they act in a combinatorial manner and control distinct events during embryo development. Seed development is a complex process that proceeds through sequences of events regulated by the interplay of various genes, prominent among them being the transcription factors (TFs). The members of WOX, HD-ZIP III, ARF, and CUC families have a preferential role in embryonic patterning. While WOX TFs are required for initiating body axis, HD-ZIP III TFs and CUCs establish bilateral symmetry and SAM. And ARF5 performs a major role during embryonic root, ground tissue, and vasculature development. TFs such as LEC1, ABI3, FUS3, and LEC2 (LAFL) are considered the master regulators of seed maturation. Furthermore, several new TFs involved in seed storage reserves and dormancy have been identified in the last few years. Their association with those master regulators has been established in the model plant Arabidopsis. Also, using chromatin immunoprecipitation (ChIP) assay coupled with transcriptomics, genome-wide target genes of these master regulators have recently been proposed. Many seed-specific genes, including those encoding oleosins and albumins, have appeared as the direct target of LAFL. Also, several other TFs act downstream of LAFL TFs and perform their function during maturation. In this review, the function of different TFs in different phases of early embryogenesis and maturation is discussed in detail, including information about their genetic and molecular interactors and target genes. Such knowledge can further be leveraged to understand and manipulate the regulatory mechanisms involved in seed development. In addition, the genomics approaches and their utilization to identify TFs aiming to study embryo development are discussed.

Transcriptome Analysis Reveals Key Pathways and Candidate Genes Controlling Seed Development and Size in Ricebean (Vigna umbellata)

Ricebean (Vigna umbellata) is a lesser known pulse with well-recognized potential. Recently, it has emerged as a legume with endowed nutritional potential because of high concentration of quality protein and other vital nutrients in its seeds. However, the genes and pathways involved in regulating seed development and size are not understood in this crop. In our study, we analyzed the transcriptome of two genotypes with contrasting grain size (IC426787: large seeded and IC552985: small seeded) at two different time points, namely, 5 and 10 days post-anthesis (DPA). The bold seeded genotype across the time points (B5_B10) revealed 6,928 differentially expressed genes (DEGs), whereas the small seeded genotype across the time point (S5_S10) contributed to 14,544 DEGs. We have also identified several candidate genes for seed development-related traits like seed size and 100-seed weight. On the basis of similarity search and domain analysis, some candidate genes (PHO1, cytokinin dehydrogenase, A-type cytokinin, and ARR response negative regulator) related to 100-seed weight and seed size showed downregulation in the small seeded genotype. The MapMan and KEGG analysis confirmed that auxin and cytokinin pathways varied in both the contrasting genotypes and can therefore be the regulators of the seed size and other seed development-related traits in ricebeans. A total of 51 genes encoding SCF TIR1/AFB , Aux/IAA, ARFs, E3 ubiquitin transferase enzyme, and 26S proteasome showing distinct expression dynamics in bold and small genotypes were also identified. We have also validated randomly selected SSR markers in eight accessions of the Vigna species (V. umbellata: 6; Vigna radiata: 1; and Vigna mungo: 1). Cross-species transferability pattern of ricebean-derived SSR markers was higher in V. radiata (73.08%) than V. mungo (50%). To the best of our knowledge, this is the first transcriptomic study conducted in this crop to understand the molecular basis of any trait. It would provide us a comprehensive understanding of the complex transcriptome dynamics during the seed development and gene regulatory mechanism of the seed size determination in ricebeans.

Characterization of Seed Proteome Profile of Wild and Cultivated Chickpeas of India

BACKGROUND

Chickpea is a widely grown legume in India, Australia, Canada, and Mediterranean regions. Seeds of chickpea are good source of protein for both human and animals. Wild relatives of chickpea (Cicer arietinum) are the potential gene pool for crop improvement; however, very little information is available on the seed proteome of these wild chickpeas.

OBJECTIVE

We aimed to analyze the seed proteome profiles of three wild relatives of chickpea, Cicer bijugum, Cicer judaicum and Cicer microphyllum along with two cultivated varieties JG11 and DCP 92/3.

METHODS

Total seed proteins were extracted using various extraction buffers for 2-D gel electrophoresis. Protein separated in a 2-D gels were subjected to image analyses, differentially expressed proteins were extracted from the gels and identified by the MALDI TOF/TOF. Seed protease inhibitors were analysed biochemically.

RESULTS

We have standardized the 2-D gel electrophoresis method and separated seed proteins using the modified method. We identified a large number (400) of protein proteins which were differentially expressed in cultivated and wild type species of chickpea. A comparative analysis between C. bijugum and JG 11 revealed the presence of 9 over-expressed and 22 under-expressed proteins, while the comparison between C. bijugum with DCP 92/3 showed 8 over-expressed and 18 under-- expressed proteins. Similarly, comparative analysis between C. microphyllum with DCP 92/3 showed 8 over-expressed proteins along with 22 under-expressed proteins, while the comparative study of C. microphyllum with JG11 displayed 9 over-expressed and 24 under-expressed proteins. We also compared C. judaicum with DCP 92/3 which revealed 15 overexpressed and 11 under-expressed proteins. On the other hand, the comparative analysis of C. judaicum with JG11 showed 10 over-expressed proteins, while the numbers of under-expressed proteins were 14. Among the differentially expressed protein proteins, 19 proteins were analyzed by the MS/MS, and peptides were identified using the MASCOT search engine. In the wild relatives the differentially expressed proteins are phosphatidylinositol 4-phosphate 5- kinase, β-1-6 galactosyltransferase, RNA helicase, phenyl alanine ammonia lyase 2, flavone 3'-0-methyl transferase, Argonaute 2, Myb related protein, Tubulin beta-2 chain and others. The most important one was legumin having α- amylase inhibition activity which was up regulated in C. bijugum. We also studied the activity of protease inhibitor (trypsin and α- amylase inhibitors) in these seed lines which showed differential activity of protease inhibitors. The highest trypsin and α- amylase inhibition was observed in C. judaicum and C. bijugum, respectively.

CONCLUSION

The differentially expressed proteins of wild relatives of chickpea appeared to be involved in various metabolic pathways. The study provides us information about the differences in the seed proteome of these wild species and cultivated varieties for the first time.

Cicer turcicum: A New Cicer Species and Its Potential to Improve Chickpea

Genetic resources of the genus Cicer L. are not only limited when compared to other important food legumes and major cereal crops but also, they include several endemic species with endangered status based on the criteria of the International Union for Conservation of Nature. The chief threats to endemic and endangered Cicer species are over-grazing and habitat change in their natural environments driven by climate changes. During a collection mission in east and south-east Anatolia (Turkey), a new Cicer species was discovered, proposed here as C. turcicum Toker, Berger & Gokturk. Here, we describe the morphological characteristics, images, and ecology of the species, and present preliminary evidence of its potential utility for chickpea improvement. C. turcicum is an annual species, endemic to southeast Anatolia and to date has only been located in a single population distant from any other known annual Cicer species. It belongs to section Cicer M. Pop. of the subgenus Pseudononis M. Pop. of the genus Cicer L. (Fabaceae) and on the basis of internal transcribed spacer (ITS) sequence similarity appears to be a sister species of C. reticulatum Ladiz. and C. echinospermum P.H. Davis, both of which are inter-fertile with domestic chickpea (C. arietinum L.). With the addition of C. turcicum, the genus Cicer now comprises 10 annual and 36 perennial species. As a preliminary evaluation of its potential for chickpea improvement two accessions of C. turcicum were field screened for reproductive heat tolerance and seeds were tested for bruchid resistance alongside a representative group of wild and domestic annual Cicer species. C. turcicum expressed the highest heat tolerance and similar bruchid resistance as C. judaicum Boiss. and C. pinnatifidum Juab. & Spach, neither of which are in the primary genepool of domestic chickpea. Given that C. arietinum and C. reticulatum returned the lowest and the second lowest tolerance and resistance scores, C. turcicum may hold much potential for chickpea improvement if its close relatedness supports interspecific hybridization with the cultigen. Crossing experiments are currently underway to explore this question.

Identification and molecular characterization of miRNAs and their target genes associated with seed development through small RNA sequencing in chickpea

Multiple studies have attempted to dissect the molecular mechanism underlying seed development in chickpea (Cicer arietinum L.). These studies highlight the need to focus on the role of miRNAs in regulating storage protein accumulation in seeds. Therefore, a total of 8,856,691 short-read sequences were generated from a small RNA library of developing chickpea seeds and were analyzed using miRDeep-P to identify 74 known and 26 novel miRNA sequences. Known miRNAs were classified into 22 miRNA families with miRNA156 family being most abundant. Of the 26 putative novel miRNAs identified, only 22 could be experimentally validated using stem loop end point PCR. Differential expression analyses led to the identification of known as well as novel miRNAs that could regulate various stages of chickpea seed development. In silico target prediction revealed several important target genes and transcription factors like SPL, mediator of RNA Polymerase II transcription subunit 12, aspartic proteinase and NACs, which were further validated by real-time PCR analysis. A comparative expression analysis in chickpea genotypes with contrasting seed protein content revealed one known (Car-miR156h) and two novel miRNA (Car-novmiR7 and Car-novmiR23) candidates to be highly expressed in the LPC (low protein content) chickpea genotypes, targets of which are known to regulate seed storage protein accumulation. Therefore, this study provides a useful resource in the form of miRNA and their targets which can be further utilized to understand and manipulate various regulatory mechanisms involved in seed development with the overall aim of improving yield and nutrition attributes in chickpea.

Transcriptomic Insights into Mechanisms of Early Seed Maturation in the Garden Pea (Pisum sativum L.)

The garden pea (Pisum sativum L.) is a legume crop of immense economic value. Extensive breeding has led to the emergence of numerous pea varieties, of which some are distinguished by accelerated development in various stages of ontogenesis. One such trait is rapid seed maturation, which, despite novel insights into the genetic control of seed development in legumes, remains poorly studied. This article presents an attempt to dissect mechanisms of early maturation in the pea line Sprint-2 by means of whole transcriptome RNA sequencing in two developmental stages. By using a de novo assembly approach, we have obtained a reference transcriptome of 25,756 non-redundant entries expressed in pea seeds at either 10 or 20 days after pollination. Differential expression in Sprint-2 seeds has affected 13,056 transcripts. A comparison of the two pea lines with a common maturation rate demonstrates that while at 10 days after pollination, Sprint-2 seeds show development retardation linked to intensive photosynthesis, morphogenesis, and cell division, and those at 20 days show a rapid onset of desiccation marked by the cessation of translation and cell anabolism and accumulation of dehydration-protective and -storage moieties. Further inspection of certain transcript functional categories, including the chromatin constituent, transcription regulation, protein turnover, and hormonal regulation, has revealed transcriptomic trends unique to specific stages and cultivars. Among other remarkable features, Sprint-2 demonstrated an enhanced expression of transposable element-associated open reading frames and an altered expression of major maturation regulators and DNA methyltransferase genes. To the best of our knowledge, this is the first comparative transcriptomic study in which the issue of the seed maturation rate is addressed.

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.

Comprehensive Stress-Based De Novo Transcriptome Assembly and Annotation of Guar (Cyamopsis tetragonoloba (L.) Taub.): An Important Industrial and Forage Crop

The forage crop Guar (Cyamopsis tetragonoloba (L.) Taub.) has the ability to endure heat, drought, and mild salinity. A complete image on its genic architecture will promote our understanding about gene expression networks and different tolerance mechanisms at the molecular level. Therefore, whole mRNA sequence approach on the Guar plant was conducted to provide a snapshot of the mRNA information in the cell under salinity, heat, and drought stresses to be integrated with previous transcriptomic studies. RNA-Seq technology was employed to perform a 2 × 100 paired-end sequencing using an Illumina HiSeq 2500 platform for the transcriptome of leaves of C. tetragonoloba under normal, heat, drought, and salinity conditions. Trinity was used to achieve a de novo assembly followed by gene annotation, functional classification, metabolic pathway analysis, and identification of SSR markers. A total of 218.2 million paired-end raw reads (~44 Gbp) were generated. Of those, 193.5M paired-end reads of high quality were used to reconstruct a total of 161,058 transcripts (~266 Mbp) with N50 of 2552 bp and 61,508 putative genes. There were 6463 proteins having >90% full-length coverage against the Swiss-Prot database and 94% complete orthologs against Embryophyta. Approximately, 62.87% of transcripts were blasted, 50.46% mapped, and 43.50% annotated. A total of 4715 InterProScan families, 3441 domains, 74 repeats, and 490 sites were detected. Biological processes, molecular functions, and cellular components comprised 64.12%, 25.42%, and 10.4%, respectively. The transcriptome was associated with 985 enzymes and 156 KEGG pathways. A total of 27,066 SSRs were gained with an average frequency of one SSR/9.825 kb in the assembled transcripts. This resulting data will be helpful for the advanced analysis of Guar to multi-stress tolerance.

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.

Resistance to Plant-Parasitic Nematodes in Chickpea: Current Status and Future Perspectives

Plant-parasitic nematodes constrain chickpea (Cicer arietinum) production, with annual yield losses estimated to be 14% of total global production. Nematode species causing significant economic damage in chickpea include root-knot nematodes (Meloidogyne artiella, M. incognita, and M. javanica), cyst nematode (Heterodera ciceri), and root-lesion nematode (Pratylenchus thornei). Reduced functionality of roots from nematode infestation leads to water stress and nutrient deficiency, which in turn lead to poor plant growth and reduced yield. Integration of resistant crops with appropriate agronomic practices is recognized as the safest and most practical, economic and effective control strategy for plant-parasitic nematodes. However, breeding for resistance to plant-parasitic nematodes has numerous challenges that originate from the narrow genetic diversity of the C. arietinum cultigen. While levels of resistance to M. artiella, H. ciceri, and P. thornei have been identified in wild Cicer species that are superior to resistance levels in the C. arietinum cultigen, barriers to interspecific hybridization restrict the use of these crop wild relatives, as sources of nematode resistance. Wild Cicer species of the primary genepool, C. reticulatum and C. echinospermum, are the only species that have been used to introgress resistance genes into the C. arietinum cultigen. The availability of genomic resources, including genome sequence and re-sequence information, the chickpea reference set and mini-core collections, and new wild Cicer collections, provide unprecedented opportunities for chickpea improvement. This review surveys progress in the identification of novel genetic sources of nematode resistance in international germplasm collections and recommends genome-assisted breeding strategies to accelerate introgression of nematode resistance into elite chickpea cultivars.

Central role of the LEAFY COTYLEDON1 transcription factor in seed development

Seed development is a complex period of the flowering plant life cycle. After fertilization, the three main regions of the seed, embryo, endosperm and seed coat, undergo a series of developmental processes that result in the production of a mature seed that is developmentally arrested, desiccated, and metabolically quiescent. These processes are highly coordinated, both temporally and spatially, to ensure the proper growth and development of the seed. The transcription factor, LEAFY COTYLEDON1 (LEC1), is a central regulator that controls several aspects of embryo and endosperm development, including embryo morphogenesis, photosynthesis, and storage reserve accumulation. Thus, LEC1 regulates distinct sets of genes at different stages of seed development. Despite its critical importance for seed development, an understanding of the mechanisms underlying LEC1's multifunctionality is only beginning to be obtained. Recent studies describe the roles of specific transcription factors and the hormones, gibberellic acid and abscisic acid, in controlling the activity and transcriptional specificity of LEC1 across seed development. Moreover, studies indicate that LEC1 acts as a pioneer transcription factor to promote epigenetic reprogramming during embryogenesis. In this review, we discuss the mechanisms that enable LEC1 to serve as a central regulator of seed development.

Overexpression of Chickpea Defensin Gene Confers Tolerance to Water-Deficit Stress in Arabidopsis thaliana

Plant defensins are mainly known for their antifungal activity. However, limited information is available regarding their function in abiotic stresses. In this study, a defensin gene, Ca-AFP, from Cicer arietinum, commonly known as chickpea, was cloned and transformed in Arabidopsis thaliana for its functional characterization under simulated water-deficit conditions. Under simulated water-deficit conditions (mannitol and polyethylene glycol-6000 induced), the transgenic A. thaliana plants had higher accumulation of the Ca-AFP transcript compared to that under non-stress condition and showed higher germination rate, root length, and biomass than the wild-type (WT) plants. To get further insights into the role of Ca-AFP in conferring tolerance to water-deficit stress, we determined various physiological parameters and found significant reduction in the transpiration rate and stomatal conductance whereas the net photosynthesis and water use efficiency was increased in the transgenic plants compared to that in the WT plants under water deficit conditions. The transgenic plants showed enhanced superoxide dismutase, ascorbate peroxidase, and catalase activities, had higher proline, chlorophyll, and relative water content, and exhibited reduced ion leakage and malondialdehyde content under water-deficit conditions. Overall, our results indicate that overexpression of Ca-AFP could be an efficient approach for conferring tolerance to water-deficit stress in plants.

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.

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.

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.

RNA-Seq analysis revealed genes associated with drought stress response in kabuli chickpea (Cicer arietinum L.)

Drought is the most important constraint that effects chickpea production globally. RNA-Seq has great potential to dissect the molecular mechanisms of tolerance to environmental stresses. Transcriptome profiles in roots and shoots of two contrasting Iranian kabuli chickpea genotypes (Bivanij and Hashem) were investigated under water-limited conditions at early flowering stage using RNA-Seq approach. A total of 4,572 differentially expressed genes (DEGs) were identified. Of these, 261 and 169 drought stress responsive genes were identified in the shoots and the roots, respectively, and 17 genes were common in the shoots and the roots. Gene Ontology (GO) analysis revealed several sub-categories related to the stress, including response to stress, defense response and response to stimulus in the tolerant genotype Bivanij as compared to the sensitive genotype Hashem under drought stress. In addition, several Transcription factors (TFs) were identified in major metabolic pathways such as, ABA, proline and flavonoid biosynthesis. Furthermore, a number of the DEGs were observed in "QTL-hotspot" regions which were reported earlier in chickpea. Drought tolerance dissection in the genotypes revealed that the genes and the pathways involved in shoots of Bivanij were the most important factor to make a difference between the genotypes for drought tolerance. The identified TFs in the experiment, particularly those which were up-regulated in shoots of Bivanij during drought stress, were potential candidates for enhancing tolerance to drought.

Isolation of Physiologically Active and Intact Mitochondria from Chickpea

Chickpea is an important leguminous crop that belongs to Fabaceae family, highly valued for its nutritious seeds. Seeds contain reserve food for the developing embryos. Mitochondria are crucial organelle for generation of chemical energy in the form of ATP which is required for achieving metabolically active state; therefore, investigating mitochondrial function and respiration rate is crucial for exploring various metabolic and physio-biochemical changes that occur during seed germination. Here we describe a method for isolation of mitochondria from germinating seeds of two chickpea varieties, i.e., Desi and Kabuli. Structure of Mitotracker-stained isolated mitochondria was observed by confocal microscopy and respiration rate was measured using an oxygen microsensor.

Comparative transcriptomics uncovers differences in photoautotrophic versus photoheterotrophic modes of nutrition in relation to secondary metabolites biosynthesis in Swertia chirayita

Swertia chirayita is a high-value medicinal herb exhibiting antidiabetic, hepatoprotective, anticancer, antiediematogenic and antipyretic properties. Scarcity of its plant material has necessitated in vitro production of therapeutic metabolites; however, their yields were low compared to field grown plants. Possible reasons for this could be differences in physiological and biochemical processes between plants grown in photoautotrophic versus photoheterotrophic modes of nutrition. Comparative transcriptomes of S. chirayita were generated to decipher the crucial molecular components associated with the secondary metabolites biosynthesis. Illumina HiSeq sequencing yielded 57,460 and 43,702 transcripts for green house grown (SCFG) and tissue cultured (SCTC) plants, respectively. Biological role analysis (GO and COG assignments) revealed major differences in SCFG and SCTC transcriptomes. KEGG orthology mapped 351 and 341 transcripts onto secondary metabolites biosynthesis pathways for SCFG and SCTC transcriptomes, respectively. Nineteen out of 30 genes from primary metabolism showed higher in silico expression (FPKM) in SCFG versus SCTC, possibly indicating their involvement in regulating the central carbon pool. In silico data were validated by RT-qPCR using a set of 16 genes, wherein 10 genes showed similar expression pattern across both the methods. Comparative transcriptomes identified differentially expressed transcription factors and ABC-type transporters putatively associated with secondary metabolism in S. chirayita. Additionally, functional classification was performed using NCBI Biosystems database. This study identified the molecular components implicated in differential modes of nutrition (photoautotrophic vs. photoheterotrophic) in relation to secondary metabolites production in S. chirayita.

Transcriptomic profiling of developing fiber in levant cotton (Gossypium herbaceum L.)

Cotton (Gossypium spp.) is an imperative economic crop of the globe due to its natural textile fiber. Molecular mechanisms of fiber development have been greatly revealed in allotetraploid cotton but remained unexplored in Gossypium herbaceum. G. herbaceum can withstand the rigors of nature like drought and pests but produce coarse lint. This undesirable characteristic strongly needs the knowledge of fiber development at molecular basis. The present study reported the transcriptome sequence of the developing fiber of G. herbaceum on pyrosequencing and its analysis. About 1.38 million raw and 1.12 million quality trimmed reads were obtained followed by de novo assembly-generated 20,125 unigenes containing 14,882 coding sequences (CDs). BLASTx-based test of homology indicated that A1-derived transcripts shared a high similarity with Gossypium arboreum (A2). Functional annotation of the CDs using the UniProt categorized them into biological processes, cellular components, and molecular function, COG classification showed that a large number of CDs have significant homology in COG database (6215 CDs), and mapping of CDs with Kyoto Encyclopedia of Genes and Genomes (KEGG) database generated 200 pathways ultimately showing predominant engagement in the fiber development process. Transcription factors were predicted by comparison with Plant Transcription Factor Database, and their differential expression between stages exposed their important regulatory role in fiber development. Differential expression analysis based on reads per kilobase of transcript per million mapped reads (RPKM) value revealed activities of specific gene related to carbohydrate and lipid synthesis, carbon metabolism, energy metabolism, signal transduction, etc., at four stages of fiber development, and was validated by qPCR. Overall, this study will help as a valuable foundation for diploid cotton fiber improvement.

RNA Sequencing and Coexpression Analysis Reveal Key Genes Involved in α-Linolenic Acid Biosynthesis in Perilla frutescens Seed

Perilla frutescen is used as traditional food and medicine in East Asia. Its seeds contain high levels of α-linolenic acid (ALA), which is important for health, but is scarce in our daily meals. Previous reports on RNA-seq of perilla seed had identified fatty acid (FA) and triacylglycerol (TAG) synthesis genes, but the underlying mechanism of ALA biosynthesis and its regulation still need to be further explored. So we conducted Illumina RNA-sequencing in seven temporal developmental stages of perilla seeds. Sequencing generated a total of 127 million clean reads, containing 15.88 Gb of valid data. The de novo assembly of sequence reads yielded 64,156 unigenes with an average length of 777 bp. A total of 39,760 unigenes were annotated and 11,693 unigenes were found to be differentially expressed in all samples. According to Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis, 486 unigenes were annotated in the “lipid metabolism” pathway. Of these, 150 unigenes were found to be involved in fatty acid (FA) biosynthesis and triacylglycerol (TAG) assembly in perilla seeds. A coexpression analysis showed that a total of 104 genes were highly coexpressed (r > 0.95). The coexpression network could be divided into two main subnetworks showing over expression in the medium or earlier and late phases, respectively. In order to identify the putative regulatory genes, a transcription factor (TF) analysis was performed. This led to the identification of 45 gene families, mainly including the AP2-EREBP, bHLH, MYB, and NAC families, etc. After coexpression analysis of TFs with highly expression of FAD2 and FAD3 genes, 162 TFs were found to be significantly associated with two FAD genes (r > 0.95). Those TFs were predicted to be the key regulatory factors in ALA biosynthesis in perilla seed. The qRT-PCR analysis also verified the relevance of expression pattern between two FAD genes and partial candidate TFs. Although it has been reported that some TFs are involved in seed development, more direct evidence is still needed to verify their function. However, these findings can provide clues to reveal the possible molecular mechanisms of ALA biosynthesis and its regulation in perilla seed.

Genome-wide analysis of the CCCH zinc finger family identifies tissue specific and stress responsive candidates in chickpea (Cicer arietinum L.)

The CCCH zinc finger is a group of proteins characterised by a typical motif consisting of three cysteine residues and one histidine residue. These proteins have been reported to play important roles in regulation of plant growth, developmental processes and environmental responses. In the present study, genome wide analysis of the CCCH zinc finger gene family was carried out in the available chickpea genome. Various bioinformatics tools were employed to predict 58 CCCH zinc finger genes in chickpea (designated CarC3H1-58), which were analysed for their physio-chemical properties. Phylogenetic analysis classified the proteins into 12 groups in which members of a particular group had similar structural organization. Further, the numbers as well as the types of CCCH motifs present in the CarC3H proteins were compared with those from Arabidopsis and Medicago truncatula. Synteny analysis revealed valuable information regarding the evolution of this gene family. Tandem and segmental duplication events were identified and their Ka/Ks values revealed that the CarC3H gene family in chickpea had undergone purifying selection. Digital, as well as real time qRT-PCR expression analysis was performed which helped in identification of several CarC3H members that expressed preferentially in specific chickpea tissues as well as during abiotic stresses (desiccation, cold, salinity). Moreover, molecular characterization of an important member CarC3H45 was carried out. This study provides comprehensive genomic information about the important CCCH zinc finger gene family in chickpea. The identified tissue specific and abiotic stress specific CCCH genes could be potential candidates for further characterization to delineate their functional roles in development and stress.

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.

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.

Deciphering Transcriptional Programming during Pod and Seed Development Using RNA-Seq in Pigeonpea (Cajanus cajan)

Seed development is an important event in plant life cycle that has interested humankind since ages, especially in crops of economic importance. Pigeonpea is an important grain legume of the semi-arid tropics, used mainly for its protein rich seeds. In order to understand the transcriptional programming during the pod and seed development, RNA-seq data was generated from embryo sac from the day of anthesis (0 DAA), seed and pod wall (5, 10, 20 and 30 DAA) of pigeonpea variety "Asha" (ICPL 87119) using Illumina HiSeq 2500. About 684 million sequencing reads have been generated from nine samples, which resulted in the identification of 27,441 expressed genes after sequence analysis. These genes have been studied for their differentially expression, co-expression, temporal and spatial gene expression. We have also used the RNA-seq data to identify important seed-specific transcription factors, biological processes and associated pathways during seed development process in pigeonpea. The comprehensive gene expression study from flowering to mature pod development in pigeonpea would be crucial in identifying candidate genes involved in seed traits directly or indirectly related to yield and quality. The dataset will serve as an important resource for gene discovery and deciphering the molecular mechanisms underlying various seed related traits.

Transcriptome landscape of perennial wild Cicer microphyllum uncovers functionally relevant molecular tags regulating agronomic traits in chickpea

The RNA-sequencing followed by de-novo transcriptome assembly identified 11621 genes differentially xpressed in roots vs. shoots of a wild perennial Cicer microphyllum. Comparative analysis of transcriptomes between microphyllum and cultivated desi cv. ICC4958 detected 12772 including 3242 root- and 1639 shoot-specific microphyllum genes with 85% expression validation success rate. Transcriptional reprogramming of microphyllum root-specific genes implicates their possible role in regulating differential natural adaptive characteristics between wild and cultivated chickpea. The transcript-derived 5698 including 282 in-silico polymorphic SSR and 127038 SNP markers annotated at a genome-wide scale exhibited high amplification and polymorphic potential among cultivated (desi and kabuli) and wild accessions suggesting their utility in chickpea genomics-assisted breeding applications. The functional significance of markers was assessed based on their localization in non-synonymous coding and regulatory regions of microphyllum root-specific genes differentially expressed predominantly in ICC 4958 roots under drought stress. A high-density 490 genic SSR- and SNP markers-anchored genetic linkage map identified six major QTLs regulating drought tolerance-related traits, yield per plant and harvest-index in chickpea. The integration of high-resolution QTL mapping with comparative transcriptome profiling delineated five microphyllum root-specific genes with non-synonymous and regulatory SNPs governing drought-responsive yield traits. Multiple potential key regulators and functionally relevant molecular tags delineated can drive translational research and drought tolerance-mediated chickpea genetic enhancement.

De novo Transcriptome Assembly and Dynamic Spatial Gene Expression Analysis in Red Clover

Red clover ( L.) is a cool-season forage legume grown throughout the northeastern United States and is the most widely planted forage legume after alfalfa ( L.). Red clover provides high-value feed to the livestock because of high protein content and easy digestibility. To date, genomic resources for red clover are scarce. In the current study, a de novo transcriptome assembly of red clover was constructed representing different tissue types. The draft assembly consists of 37,565 contigs with N50 and average contig length of 1707 and 1262 bp, respectively. A comparative study with three other legume species displayed a high degree of sequence conservation between red clover and other legumes. The assembled transcriptome was annotated to allow identification of desirable genes. In particular, a genome-wide identification of red clover transcripts encoding putative transcription factors was performed. A comparative gene expression analysis between different tissue types was performed using the assembled transcriptome as the reference, which revealed dynamic gene expression patterns across different tissue types and also identified spatially dynamic gene coexpression clusters. Genes representing tissue-enriched clusters were subjected to gene ontology (GO) enrichment analysis to identify over-represented functional groups. Identification of these tissue-enriched gene coexpression clusters can help in future research focusing on developmental studies across tissues or in biotechnological improvement of red clover.

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

A hallmark trait of chickpea (Cicer arietinum L.), like other legumes, is the capability to convert atmospheric nitrogen (N2) into ammonia (NH3) in symbiotic association with Mesorhizobium ciceri. However, the complexity of molecular networks associated with the dynamics of nodule development in chickpea need to be analyzed in depth. Hence, in order to gain insights into the chickpea nodule development, the transcriptomes of nodules at early, middle and late stages of development were sequenced using the Roche 454 platform. This generated 490.84 Mb sequence data comprising 1,360,251 reads which were assembled into 83,405 unigenes. Transcripts were annotated using Gene Ontology (GO), Cluster of Orthologous Groups (COG) and Kyoto Encyclopedia of Genes and Genomes (KEGG) metabolic pathways analysis. Differential expression analysis revealed that a total of 3760 transcripts were differentially expressed in at least one of three stages, whereas 935, 117 and 2707 transcripts were found to be differentially expressed in the early, middle and late stages of nodule development respectively. MapMan analysis revealed enrichment of metabolic pathways such as transport, protein synthesis, signaling and carbohydrate metabolism during root nodulation. Transcription factors were predicted and analyzed for their differential expression during nodule development. Putative nodule specific transcripts were identified and enriched for GO categories using BiNGO which revealed many categories to be enriched during nodule development, including transcription regulators and transporters. Further, the assembled transcriptome was also used to mine for genic SSR markers. In conclusion, this study will help in enriching the transcriptomic resources implicated in understanding of root nodulation events in chickpea.

Identification of candidate genes and natural allelic variants for QTLs governing plant height in chickpea

In the present study, molecular mapping of high-resolution plant height QTLs was performed by integrating 3625 desi genome-derived GBS (genotyping-by-sequencing)-SNPs on an ultra-high resolution intra-specific chickpea genetic linkage map (dwarf/semi-dwarf desi cv. ICC12299 x tall kabuli cv. ICC8261). The identified six major genomic regions harboring six robust QTLs (11.5-21.3 PVE), associated with plant height, were mapped within <0.5 cM average marker intervals on six chromosomes. Five SNPs-containing genes tightly linked to the five plant height QTLs, were validated based upon their high potential for target trait association (12.9-20.8 PVE) in 65 desi and kabuli chickpea accessions. The vegetative tissue-specific expression, including higher differential up-regulation (>5-fold) of five genes especially in shoot, young leaf, shoot apical meristem of tall mapping parental accession (ICC8261) as compared to that of dwarf/semi-dwarf parent (ICC12299) was apparent. Overall, combining high-resolution QTL mapping with genetic association analysis and differential expression profiling, delineated natural allelic variants in five candidate genes (encoding cytochrome-c-biosynthesis protein, malic oxidoreductase, NADH dehydrogenase iron-sulfur protein, expressed protein and bZIP transcription factor) regulating plant height in chickpea. These molecular tags have potential to dissect complex plant height trait and accelerate marker-assisted genetic enhancement for developing cultivars with desirable plant height ideotypes in chickpea.

Genome-Wide Scans for Delineation of Candidate Genes Regulating Seed-Protein Content in Chickpea

Identification of potential genes/alleles governing complex seed-protein content (SPC) is essential in marker-assisted breeding for quality trait improvement of chickpea. Henceforth, the present study utilized an integrated genomics-assisted breeding strategy encompassing trait association analysis, selective genotyping in traditional bi-parental mapping population and differential expression profiling for the first-time to understand the complex genetic architecture of quantitative SPC trait in chickpea. For GWAS (genome-wide association study), high-throughput genotyping information of 16376 genome-based SNPs (single nucleotide polymorphism) discovered from a structured population of 336 sequenced desi and kabuli accessions [with 150-200 kb LD (linkage disequilibrium) decay] was utilized. This led to identification of seven most effective genomic loci (genes) associated [10-20% with 41% combined PVE (phenotypic variation explained)] with SPC trait in chickpea. Regardless of the diverse desi and kabuli genetic backgrounds, a comparable level of association potential of the identified seven genomic loci with SPC trait was observed. Five SPC-associated genes were validated successfully in parental accessions and homozygous individuals of an intra-specific desi RIL (recombinant inbred line) mapping population (ICC 12299 × ICC 4958) by selective genotyping. The seed-specific expression, including differential up-regulation (>four fold) of six SPC-associated genes particularly in accessions, parents and homozygous individuals of the aforementioned mapping population with a high level of contrasting SPC (21-22%) was evident. Collectively, the integrated genomic approach delineated diverse naturally occurring novel functional SNP allelic variants in six potential candidate genes regulating SPC trait in chickpea. Of these, a non-synonymous SNP allele-carrying zinc finger transcription factor gene exhibiting strong association with SPC trait was found to be the most promising in chickpea. The informative functionally relevant molecular tags scaled-down essentially have potential to accelerate marker-assisted genetic improvement by developing nutritionally rich chickpea cultivars with enhanced SPC.

High-density linkage map construction and mapping of seed trait QTLs in chickpea (Cicer arietinum L.) using Genotyping-by-Sequencing (GBS)

This study reports the use of Genotyping-by-Sequencing (GBS) for large-scale SNP discovery and simultaneous genotyping of recombinant inbred lines (RILs) of an intra-specific mapping population of chickpea contrasting for seed traits. A total of 119,672 raw SNPs were discovered, which after stringent filtering revealed 3,977 high quality SNPs of which 39.5% were present in genic regions. Comparative analysis using physically mapped marker loci revealed a higher degree of synteny with Medicago in comparison to soybean. The SNP genotyping data was utilized to construct one of the most saturated intra-specific genetic linkage maps of chickpea having 3,363 mapped positions including 3,228 SNPs on 8 linkage groups spanning 1006.98 cM at an average inter marker distance of 0.33 cM. The map was utilized to identify 20 quantitative trait loci (QTLs) associated with seed traits accounting for phenotypic variations ranging from 9.97% to 29.71%. Analysis of the genomic sequence corresponding to five robust QTLs led to the identification of 684 putative candidate genes whose expression profiling revealed that 101 genes exhibited seed specific expression. The integrated approach utilizing the identified QTLs along with the available genome and transcriptome could serve as a platform for candidate gene identification for molecular breeding of chickpea.

Comprehensive transcriptome analysis of Crocus sativus for discovery and expression of genes involved in apocarotenoid biosynthesis

Background

Crocus sativus stigmas form rich source of apocarotenoids like crocin, picrocrocin and saffranal which besides imparting color, flavour and aroma to saffron spice also have tremendous pharmacological properties. Inspite of their importance, the biosynthetic pathway of Crocus apocarotenoids is not fully elucidated. Moreover, the mechanism of their stigma specific accumulation remains unknown. Therefore, deep transcriptome sequencing of Crocus stigma and rest of the flower tissue was done to identify the genes and transcriptional regulators involved in the biosynthesis of these compounds.

Results

Transcriptome of stigma and rest of the flower tissue was sequenced using Illumina Genome Analyzer IIx platform which generated 64,604,402 flower and 51,350,714 stigma reads. Sequences were assembled de novo using trinity resulting in 64,438 transcripts which were classified into 32,204 unigenes comprising of 9853 clusters and 22,351 singletons. A comprehensive functional annotation and gene ontology (GO) analysis was carried out. 58.5 % of the transcripts showed similarity to sequences present in public databases while rest could be specific to Crocus. 5789 transcripts showed similarity to transcription factors representing 76 families out of which Myb family was most abundant. Many genes involved in carotenoid/apocarotenoid pathway were identified for the first time in this study which includes zeta-carotene isomerase and desaturase, carotenoid isomerase and lycopene epsilon-cyclase. GO analysis showed that the predominant classes in biological process category include metabolic process followed by cellular process and primary metabolic process. KEGG mapping analysis indicated that pathways involved in ribosome, carbon and starch and sucrose metabolism were highly represented. Differential expression analysis indicated that key carotenoid/apocarotenoid pathway genes including phytoene synthase, phytoene desaturase and carotenoid cleavage dioxygenase 2 are enriched in stigma thereby providing molecular proof for stigma to be the site of apocarotenoid biosynthesis.

Conclusions

This data would provide a rich source for understanding the carotenoid/apocarotenoid metabolism in Crocus. The database would also help in investigating many questions related to saffron biology including flower development.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1894-5) contains supplementary material, which is available to authorized users.

Comparative Transcriptomic Analyses of Vegetable and Grain Pea (Pisum sativum L.) Seed Development

Understanding the molecular mechanisms regulating pea seed developmental process is extremely important for pea breeding. In this study, we used high-throughput RNA-Seq and bioinformatics analyses to examine the changes in gene expression during seed development in vegetable pea and grain pea, and compare the gene expression profiles of these two pea types. RNA-Seq generated 18.7 G of raw data, which were then de novo assembled into 77,273 unigenes with a mean length of 930 bp. Our results illustrate that transcriptional control during pea seed development is a highly coordinated process. There were 459 and 801 genes differentially expressed at early and late seed maturation stages between vegetable pea and grain pea, respectively. Soluble sugar and starch metabolism related genes were significantly activated during the development of pea seeds coinciding with the onset of accumulation of sugar and starch in the seeds. A comparative analysis of genes involved in sugar and starch biosynthesis in vegetable pea (high seed soluble sugar and low starch) and grain pea (high seed starch and low soluble sugar) revealed that differential expression of related genes at late development stages results in a negative correlation between soluble sugar and starch biosynthetic flux in vegetable and grain pea seeds. RNA-Seq data was validated by using real-time quantitative RT-PCR analysis for 30 randomly selected genes. To our knowledge, this work represents the first report of seed development transcriptomics in pea. The obtained results provide a foundation to support future efforts to unravel the underlying mechanisms that control the developmental biology of pea seeds, and serve as a valuable resource for improving pea breeding.