Transcriptome sequencing of wild chickpea as a rich resource for marker development

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.

Comparative Transcriptome Analyses of Different Rheum officinale Tissues Reveal Differentially Expressed Genes Associated with Anthraquinone, Catechin, and Gallic Acid Biosynthesis

Rheum officinale Baill. is an important traditional Chinese medicinal herb, its dried roots and rhizomes being widely utilized to cure diverse diseases. However, previous studies mainly focused on the active compounds and their pharmacological effects, and the molecular mechanism underlying the biosynthesis of these ingredients in R. officinale is still elusive. Here, we performed comparative transcriptome analyses to elucidate the differentially expressed genes (DEGs) in the root, stem, and leaf of R. officinale. A total of 236,031 unigenes with N50 of 769 bp was generated, 136,329 (57.76%) of which were annotated. A total of 5884 DEGs was identified after the comparative analyses of different tissues; 175 and 126 key enzyme genes with tissue-specific expression were found in the anthraquinone, catechin/gallic acid biosynthetic pathway, respectively, and some of these key enzyme genes were verified by qRT-PCR. The phylogeny of the PKS III family in Polygonaceae indicated that probably only PL_741 PKSIII1, PL_11549 PKSIII5, and PL_101745 PKSIII6 encoded PKSIII in the polyketide pathway. These results will shed light on the molecular basis of the tissue-specific accumulation and regulation of secondary metabolites in R. officinale, and lay a foundation for the future genetic diversity, molecular assisted breeding, and germplasm resource improvement of this essential medicinal plant.

Harnessing the hidden allelic diversity of wild Cicer to accelerate genomics-assisted chickpea crop improvement

Chickpea, commonly called Bengal gram or Garbanzo bean, faces a productivity crisis around the globe due to numerous biotic and abiotic stresses. The eroded genetic base of the cultivated Cicer gene pool is becoming a significant bottleneck in developing stress-resilient chickpea cultivars. In this scenario, the crop wild relatives (CWR) of chickpea, with the useful genomic wealth of their wild adaptation, give a ray of hope to improve the genetic background of the cultivated Cicer gene pool. To extrapolate these unearthed genomic diversities of wild, we require a thorough understanding of the pre-historic domestication episodes that are changing their shape with the expansion of the available scientific evidence. Keeping aforesaid in view, the current review article provides a glimpsed overview on several efforts done so far to reveal the mysterious origin and evolution of the Cicer gene pool, along with the constraints in their utilization for chickpea crop improvement. It encapsulates various stress-resilient CWR of chickpea and their use in several pre-breeding programs to develop numerous breeding populations for crop genetic enhancement. Further, this review will recapitulate the significant contributions of structural, functional and comparative genomics, pan-genomics and diverse genomics-assisted breeding strategy in dissecting the untapped trait-specific allelic/gene diversity and domestication pattern behind the CWR of chickpea, along with their potential and promises. We expect the newly explored genetic variations may be used in the breeding programs for re-wilding the cultigens' genomic background to open a new avenue for genetic gain and crop improvement capacity of chickpea.

Genomics Associated Interventions for Heat Stress Tolerance in Cool Season Adapted Grain Legumes

Cool season grain legumes occupy an important place among the agricultural crops and essentially provide multiple benefits including food supply, nutrition security, soil fertility improvement and revenue for farmers all over the world. However, owing to climate change, the average temperature is steadily rising, which negatively affects crop performance and limits their yield. Terminal heat stress that mainly occurred during grain development phases severely harms grain quality and weight in legumes adapted to the cool season, such as lentils, faba beans, chickpeas, field peas, etc. Although, traditional breeding approaches with advanced screening procedures have been employed to identify heat tolerant legume cultivars. Unfortunately, traditional breeding pipelines alone are no longer enough to meet global demands. Genomics-assisted interventions including new-generation sequencing technologies and genotyping platforms have facilitated the development of high-resolution molecular maps, QTL/gene discovery and marker-assisted introgression, thereby improving the efficiency in legumes breeding to develop stress-resilient varieties. Based on the current scenario, we attempted to review the intervention of genomics to decipher different components of tolerance to heat stress and future possibilities of using newly developed genomics-based interventions in cool season adapted grain legumes.

Whole genome sequence analysis of rice genotypes with contrasting response to salinity stress

Salinity is a major abiotic constraint for rice farming. Abundant natural variability exists in rice germplasm for salt tolerance traits. Since few studies focused on the genome level variation in rice genotypes with contrasting response to salt stress, genomic resequencing in diverse genetic materials is needed to elucidate the molecular basis of salt tolerance mechanisms. The whole genome sequences of two salt tolerant (Pokkali and Nona Bokra) and three salt sensitive (Bengal, Cocodrie, and IR64) rice genotypes were analyzed. A total of 413 million reads were generated with a mean genome coverage of 93% and mean sequencing depth of 18X. Analysis of the DNA polymorphisms revealed that 2347 nonsynonymous SNPs and 51 frameshift mutations could differentiate the salt tolerant from the salt sensitive genotypes. The integration of genome-wide polymorphism information with the QTL mapping and expression profiling data led to identification of 396 differentially expressed genes with large effect variants in the coding regions. These genes were involved in multiple salt tolerance mechanisms, such as ion transport, oxidative stress tolerance, signal transduction, and transcriptional regulation. The genome-wide DNA polymorphisms and the promising candidate genes identified in this study represent a valuable resource for molecular breeding of salt tolerant rice varieties.

Updates on Genomic Resources in Chickpea for Crop Improvement

In recent years, rapid advancement has been done in generation of genomic resources for the important legume crop chickpea. Here, we provide an update on important advancements made on availability of genomic resources for this crop. The availability of reference genome and transcriptome sequences, and resequencing of several accessions have enabled the discovery of gene space and molecular markers in chickpea. These resources have helped in elucidating evolutionary relationship and identification of quantitative trait loci for important agronomic traits. Gene expression in different tissues/organs during development and under abiotic/biotic stresses has been interrogated. In addition, single-base resolution DNA methylation patterns in different organs have been analyzed to understand gene regulation. Overall, we provide a consolidated overview of available genomic resources of chickpea that may help in fulfilling the promises for improvement of this important crop.

The genetic diversity and population structure of Sophora alopecuroides (Faboideae) as determined by microsatellite markers developed from transcriptome

Sophora alopecuroides (Faboideae) is an endemic species, mainly distributed in northwest China. However, the limited molecular markers range for this species hinders breeding and genetic studies. A total of 20,324 simple sequence repeat (SSR) markers were identified from 118,197 assembled transcripts and 18 highly polymorphic SSR markers were used to explore the genetic diversity and population structure of S. alopecuroides from 23 different geographical populations. A relatively low genetic diversity was found in S. alopecuroides based on mean values of the number of effective alleles (Ne = 1.81), expected heterozygosity (He = 0.39) and observed heterozygosity (Ho = 0.55). The results of AMOVA indicated higher levels of variation within populations than between populations. Bayesian-based cluster analysis, principal coordinates analysis and Neighbor-Joining phylogeny analysis roughly divided all genotypes into four major groups with some admixtures. Meanwhile, geographic barriers would have restricted gene flow between the northern and southern regions (separated by Tianshan Mountains), wherein the two relatively ancestral and independent clusters of S. alopecuroides occur. History trade and migration along the Silk Road would together have promoted the spread of S. alopecuroides from the western to the eastern regions of the northwest plateau in China, resulting in the current genetic diversity and population structure. The transcriptomic SSR markers provide a valuable resource for understanding the genetic diversity and population structure of S. alopecuroides, and will assist effective conservation management.

Transcript profiling of chickpea pod wall revealed the expression of floral homeotic gene AGAMOUS-like X2 (CaAGLX2)

The differentially expressed genes in the chickpea pod wall have been identified for the first time using a forward suppression subtractive hybridization (SSH) library. In all, 226 clones of SSH library were sequenced and analyzed. A total of 179 high-quality expressed sequence tags (ESTs) were generated and based on the CAP3 assembly of these ESTs, 126 genes (97 singletons and 29 contigs) were computationally annotated. The mapping of 88.26% ESTs by gene ontology (GO) annotation distributed them into 751 GO terms of three categories, cellular location, molecular function, and biological process. The KEGG pathway analysis revealed 45 ESTs are involved in 49 different biological pathways. Also, 67 ESTs encodes four different classes of enzymes such as oxidoreductases (29), transferase (20), hydrolases (16) and isomerase (2). Six genes were selected and subjected to qPCR analysis, of these, two genes (FHG Floral homeotic AGAMOUS-like isoform X2, MADS1 MADS-box transcription factor) showed significant up-regulation in the pod wall compared to leaves. Surprisingly, one of the MADS1 box gene, FHG (CaAGLX2), responsible for flower development expressed in the pod wall. Therefore, understanding its specific role in the pod wall could be interesting. Thus, the transcript dynamics of the chickpea pod wall revealed differentially expressed genes in the pod wall, which may be participating in the metabolic build-up of both pod wall and seeds.

Transcriptome sequencing, de novo assembly, characterisation of wild accession of blackgram (Vigna mungo var. silvestris) as a rich resource for development of molecular markers and validation of SNPs by high resolution melting (HRM) analysis

BACKGROUND

Blackgram [Vigna mungo (L.) Hepper], is an important legume crop of Asia with limited genomic resources. We report a comprehensive set of genic simple sequence repeat (SSR) and single nucleotide polymorphism (SNPs) markers using Illumina MiSeq sequencing of transcriptome and its application in genetic variation analysis and mapping.

RESULTS

Transcriptome sequencing of immature seeds of wild blackgram, V. mungo var. silvestris by Illumina MiSeq technology generated 1.9 × 10⁷ reads, which were assembled into 40,178 transcripts (TCS) with an average length of 446 bp covering 2.97 GB of the genome. A total of 38,753 CDS (Coding sequences) were predicted from 40,178 TCS and 28,984 CDS were annotated through BLASTX and mapped to GO and KEGG database resulting in 140 unique pathways. The tri-nucleotides were most abundant (39.9%) followed by di-nucleotide (30.2%). About 60.3 and 37.6% of SSR motifs were present in the coding sequences (CDS) and untranslated regions (UTRs) respectively. Among SNPs, the most abundant substitution type were transitions (Ts) (61%) followed by transversions (Tv) type (39%), with a Ts/Tv ratio of 1.58. A total of 2306 DEGs were identified by RNA Seq between wild and cultivar and validation was done by quantitative reverse transcription polymerase chain reaction. In this study, we genotyped SNPs with a validation rate of 78.87% by High Resolution Melting (HRM) Assay.

CONCLUSION

In the present study, 1621genic-SSR and 1844 SNP markers were developed from immature seed transcriptome sequence of blackgram and 31 genic-SSR markers were used to study genetic variations among different blackgram accessions. Above developed markers contribute towards enriching available genomic resources for blackgram and aid in breeding programmes.

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.

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.

Genome-wide discovery of DNA polymorphisms among chickpea cultivars with contrasting seed size/weight and their functional relevance

Seed size/weight is a major agronomic trait which determine crop productivity in legumes. To understand the genetic basis of seed size determination, we sought to identify DNA polymorphisms between two small (Himchana 1 and Pusa 362) and two large-seeded (JGK 3 and PG 0515) chickpea cultivars via whole genome resequencing. We identified a total of 75535 single nucleotide polymorphisms (SNPs), 6486 insertions and deletions (InDels), 1938 multi-nucleotide polymorphisms (MNPs) and 5025 complex variants between the two small and two large-seeded chickpea cultivars. Our analysis revealed 814, 244 and 72 seed-specific genes harboring DNA polymorphisms in promoter or non-synonymous and large-effect DNA polymorphisms, respectively. Gene ontology analysis revealed enrichment of cell growth and division related terms in these genes. Among them, at least 22 genes associated with quantitative trait loci, and those involved in cell growth and division and encoding transcription factors harbored promoter and/or large-effect/non-synonymous DNA polymorphisms. These also showed higher expression at late-embryogenesis and/or mid-maturation stages of seed development in the large-seeded cultivar, suggesting their role in seed size/weight determination in chickpea. Altogether, this study provided a valuable resource for large-scale genotyping applications and a few putative candidate genes that might play crucial role in governing seed size/weight in chickpea.

Genome-wide discovery of DNA polymorphisms by whole genome sequencing differentiates weedy and cultivated rice

Analyzing the genome level DNA polymorphisms between weedy and cultivated rice is crucial to elucidate the molecular basis of weedy and agronomic traits, which in turn can enhance our ability to control weedy rice and its utilization for rice improvement. Here, we presented the genome-wide genetic variations between a weedy rice accession PSRR-1 and two cultivated rice accessions, Bengal and Nona Bokra, belonging to japonica and indica subspecies, respectively. The total number of SNPs and InDels in PSRR/Bengal was similar to that of Nona Bokra/Bengal, but was three times greater than that of PSRR/Nona Bokra. There were 11546 large-effect SNPs/InDels affecting 5673 genes, which most likely differentiated weedy rice from cultivated rice. These large effect DNA polymorphisms were mostly resulted in stop codon gain and least by start codon loss. Analysis of the molecular functions and biological processes of weedy rice specific SNPs/InDels indicated that most of these genes were involved in protein modification/phosphorylation, protein kinase activity, and protein/nucleotide binding. By integrating previous QTL mapping results with the DNA polymorphisms data, the candidate genes for seed dormancy and seed shattering were narrowed down. The genomic resource generated in this study will facilitate discovery of functional variants for weedy and agronomic traits.

Genome-wide generation and genotyping of informative SNPs to scan molecular signatures for seed yield in chickpea

We discovered 2150 desi and 2199 kabuli accessions-derived SNPs by cultivar-wise individual assembling of sequence-reads generated through genotyping-by-sequencing of 92 chickpea accessions. Subsequent large-scale validation and genotyping of these SNPs discovered 619 desi accessions-derived (DAD) SNPs, 531 kabuli accessions-derived (KAD) SNPs, 884 multiple accessions-derived (MAD) SNPs and 1083 two accessions (desi ICC 4958 and kabuli CDC Frontier)-derived (TAD) SNPs that were mapped on eight chromosomes. These informative SNPs were annotated in coding/non-coding regulatory sequence components of genes. The MAD-SNPs were efficient to detect high intra-specific polymorphic potential and wide natural allelic diversity level including high-resolution admixed-population genetic structure and precise phylogenetic relationship among 291 desi and kabuli accessions. This signifies their effectiveness in introgression breeding and varietal improvement studies targeting useful agronomic traits of chickpea. Six trait-associated genes with SNPs including quantitative trait nucleotides (QTNs) in combination explained 27.5% phenotypic variation for seed yield per plant (SYP). A pentatricopeptide repeat (PPR) gene with a synonymous-coding SNP/QTN significantly associated with SYP trait was found most-promising in chickpea. The essential information delineated can be of immense utility in genomics-assisted breeding applications to develop high-yielding chickpea cultivars.

Development of Microsatellite Markers for the Nipa Palm Hispid Beetle, Octodonta nipae (Maulik)

The nipa palm hispid beetle, Octodonta nipae (Maulik) (Coleoptera: Chrysomelidae), is an important invasive pest on palm plants in southern China. Based on existing transcriptome data, polymorphism simple sequence repeat (SSR) loci were identified. In total, 1274 SSR loci were identified from 49919 unigenes. The majority of them contained mononucleotide, dinucleotide, and trinucleotide motifs (43.56%, 26.14%, and 28.18%), in which A/T (41.21%) and AT/TA (15.86%) were the most abundant motifs. 104 pairs of the SSR primers produced amplification bands of expected sizes in O. nipae, 80 pairs of SSR primers were tested randomly for polymorphism, 9 loci of them were validated to be polymorphic markers, and the number of alleles ranged from 2 to 3, with an average of 2.56 per locus. The population of Zhangzhou and Fuzhou was analyzed by the 9 loci (On1-On9). These SSR transcriptome data can provide invaluable resource for SSR development, population genetics research, invasion and expansion mechanism, paternity testing, and other research on O. nipae and its related species.

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.

Transcriptome analysis of near-isogenic line provides novel insights into genes associated with panicle traits regulation in rice

Panicle traits in rice impact yield and quality. The OsGRF4 gene encodes a growth-regulating factor controlling panicle traits, and was recently cloned. Gene expression profiling analysis can be used to study the molecular mechanisms underlying OsGRF4 regulation. Use of near-isogenic lines (NILs) reduces genetic background noise in omics studies. We compared transcriptome profiling of 7 cm long young panicles of NIL-Osgrf4 and NIL-OsGRF4 using RNAs sequence analyses. Eighty differentially expressed genes (DEGs) were identified. Our target gene OsGRF4 was up-regulated in NIL-OsGRF4 plants, which is consistent with a previous qPCR analysis. Hierarchical cluster analysis showed OsGRF4 is tightly clustered with the up-regulated DEG LOC_Os02g47320. Gene Ontology (GO) and KEGG analysis suggested that DEGs were primarily involved in somatic embryogenesis and chitinase activity. Two up-regulated DEGs, LOC_Os04g41680 and LOC_Os04g41620, were significantly enriched in the top 8 GO terms, and were over_represented in term of seed development, and may play key roles in grain shape regulation. The transcription factor Osmyb1 also exhibited differential expression between NILs, and may be is an important regulator of panicle traits. By searching reported functions of DEGs and by co-localization with previous identified quantitative trait loci (QTL), we determined that the pleiotropic gene OsGRF4 may also be involve in abiotic stress resistance. This study provides new candidates genes for further understanding the molecular mechanisms underlying rice panicle trait regulation.

Transcriptomic resources for the medicinal legume Mucuna pruriens: de novo transcriptome assembly, annotation, identification and validation of EST-SSR markers

Background

The medicinal legume Mucuna pruriens (L.) DC. has attracted attention worldwide as a source of the anti-Parkinson’s drug L-Dopa. It is also a popular green manure cover crop that offers many agronomic benefits including high protein content, nitrogen fixation and soil nutrients. The plant currently lacks genomic resources and there is limited knowledge on gene expression, metabolic pathways, and genetics of secondary metabolite production. Here, we present transcriptomic resources for M. pruriens, including a de novo transcriptome assembly and annotation, as well as differential transcript expression analyses between root, leaf, and pod tissues. We also develop microsatellite markers and analyze genetic diversity and population structure within a set of Indian germplasm accessions.

Results

One-hundred ninety-one million two hundred thirty-three thousand two hundred forty-two bp cleaned reads were assembled into 67,561 transcripts with mean length of 626 bp and N50 of 987 bp. Assembled sequences were annotated using BLASTX against public databases with over 80% of transcripts annotated. We identified 7,493 simple sequence repeat (SSR) motifs, including 787 polymorphic repeats between the parents of a mapping population. 134 SSRs from expressed sequenced tags (ESTs) were screened against 23 M. pruriens accessions from India, with 52 EST-SSRs retained after quality control. Population structure analysis using a Bayesian framework implemented in fastSTRUCTURE showed nearly similar groupings as with distance-based (neighbor-joining) and principal component analyses, with most of the accessions clustering per geographical origins. Pair-wise comparison of transcript expression in leaves, roots and pods identified 4,387 differentially expressed transcripts with the highest number occurring between roots and leaves. Differentially expressed transcripts were enriched with transcription factors and transcripts annotated as belonging to secondary metabolite pathways.

Conclusions

The M. pruriens transcriptomic resources generated in this study provide foundational resources for gene discovery and development of molecular markers. Polymorphic SSRs identified can be used for genetic diversity, marker-trait analyses, and development of functional markers for crop improvement. The results of differential expression studies can be used to investigate genes involved in L-Dopa synthesis and other key metabolic pathways in M. pruriens.

Electronic supplementary material

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

Draft genome sequence of Cicer reticulatum L., the wild progenitor of chickpea provides a resource for agronomic trait improvement

Cicer reticulatum L. is the wild progenitor of the fourth most important legume crop chickpea (C. arietinum L.). We assembled short-read sequences into 416 Mb draft genome of C. reticulatum and anchored 78% (327 Mb) of this assembly to eight linkage groups. Genome annotation predicted 25,680 protein-coding genes covering more than 90% of predicted gene space. The genome assembly shared a substantial synteny and conservation of gene orders with the genome of the model legume Medicago truncatula. Resistance gene homologs of wild and domesticated chickpeas showed high sequence homology and conserved synteny. Comparison of gene sequences and nucleotide diversity using 66 wild and domesticated chickpea accessions suggested that the desi type chickpea was genetically closer to the wild species than the kabuli type. Comparative analyses predicted gene flow between the wild and the cultivated species during domestication. Molecular diversity and population genetic structure determination using 15,096 genome-wide single nucleotide polymorphisms revealed an admixed domestication pattern among cultivated (desi and kabuli) and wild chickpea accessions belonging to three population groups reflecting significant influence of parentage or geographical origin for their cultivar-specific population classification. The assembly and the polymorphic sequence resources presented here would facilitate the study of chickpea domestication and targeted use of wild Cicer germplasms for agronomic trait improvement in chickpea.

Development of molecular map and identification of QTLs linked to Fusarium wilt resistance in chickpea

A number of genetic maps for Fusarium wilt resistance in chickpea have been reported in earlier studies, however QTLs identified for Fusarium wilt resistance were unstable. Hence, the present study aims to map novel molecular markers and to identify QTLs for Fusarium wilt resistance in chickpea. An intraspecific linkage map of chickpea (Cicer arietinum L.) was constructed using F10-F11 recombinant inbred lines (RILs) derived from a cross between K850 and WR315 segregating for H2 locus. A set of 31 polymorphic simple sequence repeat (SSR) markers obtained by screening 300 SSRs and were used for genotyping. The linkage map had four linkage groups and coverage of 690 cM with a marker density of 5.72 cM. The RILs were screened for their wilt reaction across two seasons in wilt sick plot at International Crop Research Institute for Semi-Arid Tropics (ICRISAT), Hyderabad, India. Five major quantitative trait loci (QTLs) were detected in both seasons for late wilting (60 days after sowing). A stable QTL (GSSR 18-TC14801) for wilt resistance was identified in both the seasons, and the QTL explained a variance of 69.80 and 60.80% in 2007 and 2008 rabi respectively.

Genome-wide development and deployment of informative intron-spanning and intron-length polymorphism markers for genomics-assisted breeding applications in chickpea

The discovery and large-scale genotyping of informative gene-based markers is essential for rapid delineation of genes/QTLs governing stress tolerance and yield component traits in order to drive genetic enhancement in chickpea. A genome-wide 119169 and 110491 ISM (intron-spanning markers) from 23129 desi and 20386 kabuli protein-coding genes and 7454 in silico InDel (insertion-deletion) (1-45-bp)-based ILP (intron-length polymorphism) markers from 3283 genes were developed that were structurally and functionally annotated on eight chromosomes and unanchored scaffolds of chickpea. A much higher amplification efficiency (83%) and intra-specific polymorphic potential (86%) detected by these markers than that of other sequence-based genetic markers among desi and kabuli chickpea accessions was apparent even by a cost-effective agarose gel-based assay. The genome-wide physically mapped 1718 ILP markers assayed a wider level of functional genetic diversity (19-81%) and well-defined phylogenetics among domesticated chickpea accessions. The gene-derived 1424 ILP markers were anchored on a high-density (inter-marker distance: 0.65cM) desi intra-specific genetic linkage map/functional transcript map (ICC 4958×ICC 2263) of chickpea. This reference genetic map identified six major genomic regions harbouring six robust QTLs mapped on five chromosomes, which explained 11-23% seed weight trait variation (7.6-10.5 LOD) in chickpea. The integration of high-resolution QTL mapping with differential expression profiling detected six including one potential serine carboxypeptidase gene with ILP markers (linked tightly to the major seed weight QTLs) exhibiting seed-specific expression as well as pronounced up-regulation especially in seeds of high (ICC 4958) as compared to low (ICC 2263) seed weight mapping parental accessions. The marker information generated in the present study was made publicly accessible through a user-friendly web-resource, "Chickpea ISM-ILP Marker Database". The designing of multiple ISM and ILP markers (2-5 markers/gene) from an individual gene (transcription factor) with numerous aforementioned desirable genetic attributes can widen the user-preference to select suitable primer combination for simultaneous large-scale assaying of functional allelic variation, natural allelic diversity, molecular mapping and expression profiling of genes among chickpea accessions. This will essentially accelerate the identification of functionally relevant molecular tags regulating vital agronomic traits for genomics-assisted crop improvement by optimal resource expenses in chickpea.

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.

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.

Transcriptome Characterization and Functional Marker Development in Sorghum Sudanense

Sudangrass, Sorghum sudanense, is an important forage in warm regions. But little is known about its genome. In this study, the transcriptomes of sudangrass S722 and sorghum Tx623B were sequenced by Illumina sequencing. More than 4Gb bases were sequenced for each library. For Tx623B and S722, 88.79% and 83.88% reads, respectively were matched to the Sorghum bicolor genome. A total of 2,397 differentially expressed genes (DEGs) were detected by RNA-Seq between the two libraries, including 849 up-regulated genes and 1,548 down-regulated genes. These DEGs could be divided into three groups by annotation analysis. A total of 44,495 single nucleotide polymorphisms (SNPs) were discovered by aligning S722 reads to the sorghum reference genome. Of these SNPs, 61.37% were transition, and this value did not differ much between different chromosomes. In addition, 16,928 insertion and deletion (indel) loci were identified between the two genomes. A total of 5,344 indel markers were designed, 15 of which were selected to construct the genetic map derived from the cross of Tx623A and Sa. It was indicated that the indel markers were useful and versatile between sorghum and sudangrass. Comparison of synonymous base substitutions (Ks) and non-synonymous base substitutions (Ka) between the two libraries showed that 95% orthologous pairs exhibited Ka/Ks<1.0, indicating that these genes were influenced by purifying selection. The results from this study provide important information for molecular genetic research and a rich resource for marker development in sudangrass and other Sorghum species.

EcoTILLING-Based Association Mapping Efficiently Delineates Functionally Relevant Natural Allelic Variants of Candidate Genes Governing Agronomic Traits in Chickpea

The large-scale mining and high-throughput genotyping of novel gene-based allelic variants in natural mapping population are essential for association mapping to identify functionally relevant molecular tags governing useful agronomic traits in chickpea. The present study employs an alternative time-saving, non-laborious and economical pool-based EcoTILLING approach coupled with agarose gel detection assay to discover 1133 novel SNP allelic variants from diverse coding and regulatory sequence components of 1133 transcription factor (TF) genes by genotyping in 192 diverse desi and kabuli chickpea accessions constituting a seed weight association panel. Integrating these SNP genotyping data with seed weight field phenotypic information of 192 structured association panel identified eight SNP alleles in the eight TF genes regulating seed weight of chickpea. The associated individual and combination of all SNPs explained 10-15 and 31% phenotypic variation for seed weight, respectively. The EcoTILLING-based large-scale allele mining and genotyping strategy implemented for association mapping is found much effective for a diploid genome crop species like chickpea with narrow genetic base and low genetic polymorphism. This optimized approach thus can be deployed for various genomics-assisted breeding applications with optimal expense of resources in domesticated chickpea. The seed weight-associated natural allelic variants and candidate TF genes delineated have potential to accelerate marker-assisted genetic improvement of chickpea.

Genomics of crop wild relatives: expanding the gene pool for crop improvement

Plant breeders require access to new genetic diversity to satisfy the demands of a growing human population for more food that can be produced in a variable or changing climate and to deliver the high-quality food with nutritional and health benefits demanded by consumers. The close relatives of domesticated plants, crop wild relatives (CWRs), represent a practical gene pool for use by plant breeders. Genomics of CWR generates data that support the use of CWR to expand the genetic diversity of crop plants. Advances in DNA sequencing technology are enabling the efficient sequencing of CWR and their increased use in crop improvement. As the sequencing of genomes of major crop species is completed, attention has shifted to analysis of the wider gene pool of major crops including CWR. A combination of de novo sequencing and resequencing is required to efficiently explore useful genetic variation in CWR. Analysis of the nuclear genome, transcriptome and maternal (chloroplast and mitochondrial) genome of CWR is facilitating their use in crop improvement. Genome analysis results in discovery of useful alleles in CWR and identification of regions of the genome in which diversity has been lost in domestication bottlenecks. Targeting of high priority CWR for sequencing will maximize the contribution of genome sequencing of CWR. Coordination of global efforts to apply genomics has the potential to accelerate access to and conservation of the biodiversity essential to the sustainability of agriculture and food production.

Genome-wide generation and use of informative intron-spanning and intron-length polymorphism markers for high-throughput genetic analysis in rice

We developed genome-wide 84634 ISM (intron-spanning marker) and 16510 InDel-fragment length polymorphism-based ILP (intron-length polymorphism) markers from genes physically mapped on 12 rice chromosomes. These genic markers revealed much higher amplification-efficiency (80%) and polymorphic-potential (66%) among rice accessions even by a cost-effective agarose gel-based assay. A wider level of functional molecular diversity (17-79%) and well-defined precise admixed genetic structure was assayed by 3052 genome-wide markers in a structured population of indica, japonica, aromatic and wild rice. Six major grain weight QTLs (11.9-21.6% phenotypic variation explained) were mapped on five rice chromosomes of a high-density (inter-marker distance: 0.98 cM) genetic linkage map (IR 64 x Sonasal) anchored with 2785 known/candidate gene-derived ISM and ILP markers. The designing of multiple ISM and ILP markers (2 to 4 markers/gene) in an individual gene will broaden the user-preference to select suitable primer combination for efficient assaying of functional allelic variation/diversity and realistic estimation of differential gene expression profiles among rice accessions. The genomic information generated in our study is made publicly accessible through a user-friendly web-resource, "Oryza ISM-ILP marker" database. The known/candidate gene-derived ISM and ILP markers can be enormously deployed to identify functionally relevant trait-associated molecular tags by optimal-resource expenses, leading towards genomics-assisted crop improvement in rice.

Transcriptome Sequencing and Development of Genic SSR Markers of an Endangered Chinese Endemic Genus Dipteronia Oliver (Aceraceae)

Dipteronia Oliver (Aceraceae) is an endangered Chinese endemic genus consisting of two living species, Dipteronia sinensis and Dipteronia dyeriana. However, studies on the population genetics and evolutionary analyses of Dipteronia have been hindered by limited genomic resources and genetic markers. Here, the generation, de novo assembly and annotation of transcriptome datasets, and a large set of microsatellite or simple sequence repeat (SSR) markers derived from Dipteronia have been described. After Illumina pair-end sequencing, approximately 93.2 million reads were generated and assembled to yield a total of 99,358 unigenes. A majority of these unigenes (53%, 52,789) had at least one blast hit against the public protein databases. Further, 12,377 SSR loci were detected and 4179 primer pairs were designed for experimental validation. Of these 4179 primer pairs, 435 primer pairs were randomly selected to test polymorphism. Our results show that products from 132 primer pairs were polymorphic, in which 97 polymorphic SSR markers were further selected to analyze the genetic diversity of 10 natural populations of Dipteronia. The identification of SSR markers during our research will provide the much valuable data for population genetic analyses and evolutionary studies in Dipteronia.

Transcriptome analyses reveal genotype- and developmental stage-specific molecular responses to drought and salinity stresses in chickpea

Drought and salinity are the major factors that limit chickpea production worldwide. We performed whole transcriptome analyses of chickpea genotypes to investigate the molecular basis of drought and salinity stress response/adaptation. Phenotypic analyses confirmed the contrasting responses of the chickpea genotypes to drought or salinity stress. RNA-seq of the roots of drought and salinity related genotypes was carried out under control and stress conditions at vegetative and/or reproductive stages. Comparative analysis of the transcriptomes revealed divergent gene expression in the chickpea genotypes at different developmental stages. We identified a total of 4954 and 5545 genes exclusively regulated in drought-tolerant and salinity-tolerant genotypes, respectively. A significant fraction (~47%) of the transcription factor encoding genes showed differential expression under stress. The key enzymes involved in metabolic pathways, such as carbohydrate metabolism, photosynthesis, lipid metabolism, generation of precursor metabolites/energy, protein modification, redox homeostasis and cell wall component biogenesis, were affected by drought and/or salinity stresses. Interestingly, transcript isoforms showed expression specificity across the chickpea genotypes and/or developmental stages as illustrated by the AP2-EREBP family members. Our findings provide insights into the transcriptome dynamics and components of regulatory network associated with drought and salinity stress responses in chickpea.

An Efficient Strategy Combining SSR Markers- and Advanced QTL-seq-driven QTL Mapping Unravels Candidate Genes Regulating Grain Weight in Rice

Development and use of genome-wide informative simple sequence repeat (SSR) markers and novel integrated genomic strategies are vital to drive genomics-assisted breeding applications and for efficient dissection of quantitative trait loci (QTLs) underlying complex traits in rice. The present study developed 6244 genome-wide informative SSR markers exhibiting in silico fragment length polymorphism based on repeat-unit variations among genomic sequences of 11 indica, japonica, aus, and wild rice accessions. These markers were mapped on diverse coding and non-coding sequence components of known cloned/candidate genes annotated from 12 chromosomes and revealed a much higher amplification (97%) and polymorphic potential (88%) along with wider genetic/functional diversity level (16-74% with a mean 53%) especially among accessions belonging to indica cultivar group, suggesting their utility in large-scale genomics-assisted breeding applications in rice. A high-density 3791 SSR markers-anchored genetic linkage map (IR 64 × Sonasal) spanning 2060 cM total map-length with an average inter-marker distance of 0.54 cM was generated. This reference genetic map identified six major genomic regions harboring robust QTLs (31% combined phenotypic variation explained with a 5.7-8.7 LOD) governing grain weight on six rice chromosomes. One strong grain weight major QTL region (OsqGW5.1) was narrowed-down by integrating traditional QTL mapping with high-resolution QTL region-specific integrated SSR and single nucleotide polymorphism markers-based QTL-seq analysis and differential expression profiling. This led us to delineate two natural allelic variants in two known cis-regulatory elements (RAV1AAT and CARGCW8GAT) of glycosyl hydrolase and serine carboxypeptidase genes exhibiting pronounced seed-specific differential regulation in low (Sonasal) and high (IR 64) grain weight mapping parental accessions. Our genome-wide SSR marker resource (polymorphic within/between diverse cultivar groups) and integrated genomic strategy can efficiently scan functionally relevant potential molecular tags (markers, candidate genes and alleles) regulating complex agronomic traits (grain weight) and expedite marker-assisted genetic enhancement in rice.

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.

De novo Transcriptome Analysis and Molecular Marker Development of Two Hemarthria Species

Hemarthria R. Br. is an important genus of perennial forage grasses that is widely used in subtropical and tropical regions. Hemarthria grasses have made remarkable contributions to the development of animal husbandry and agro-ecosystem maintenance; however, there is currently a lack of comprehensive genomic data available for these species. In this study, we used Illumina high-throughput deep sequencing to characterize of two agriculturally important Hemarthria materials, H. compressa "Yaan" and H. altissima "1110." Sequencing runs that used each of four normalized RNA samples from the leaves or roots of the two materials yielded more than 24 million high-quality reads. After de novo assembly, 137,142 and 77,150 unigenes were obtained for "Yaan" and "1110," respectively. In addition, a total of 86,731 "Yaan" and 48,645 "1110" unigenes were successfully annotated. After consolidating the unigenes for both materials, 42,646 high-quality SNPs were identified in 10,880 unigenes and 10,888 SSRs were identified in 8330 unigenes. To validate the identified markers, high quality PCR primers were designed for both SNPs and SSRs. We randomly tested 16 of the SNP primers and 54 of the SSR primers and found that the majority of these primers successfully amplified the desired PCR product. In addition, high cross-species transferability (61.11-87.04%) of SSR markers was achieved for four other Poaceae species. The amount of RNA sequencing data that was generated for these two Hemarthria species greatly increases the amount of genomic information available for Hemarthria and the SSR and SNP markers identified in this study will facilitate further advancements in genetic and molecular studies of the Hemarthria genus.

Characterization of Global Transcriptome Using Illumina Paired-End Sequencing and Development of EST-SSR Markers in Two Species of Gynostemma (Cucurbitaceae)

Gynostemma pentaphyllum is an important medicinal herb of the Cucurbitaceae family, but limited genomic data have hindered genetic studies. In this study, transcriptomes of two closely-related Gynostemma species, Gynostemma cardiospermum and G. pentaphyllum, were sequenced using Illumina paired-end sequencing technology. A total of 71,607 nonredundant unigenes were assembled. Of these unigenes, 60.45% (43,288) were annotated based on sequence similarity search with known proteins. A total of 11,059 unigenes were identified in the Kyoto Encyclopedia of Genes and Genomes Pathway (KEGG) database. A total of 3891 simple sequence repeats (SSRs) were detected in 3526 nonredundant unigenes, 2596 primer pairs were designed and 360 of them were randomly selected for validation. Of these, 268 primer pairs yielded clear products among six G. pentaphyllum samples. Thirty polymorphic SSR markers were used to test polymorphism and transferability in Gynostemma. Finally, 15 SSR makers that amplified in all 12 Gynostemma species were used to assess genetic diversity. Our results generated a comprehensive sequence resource for Gynostemma research.

A genome-scale integrated approach aids in genetic dissection of complex flowering time trait in chickpea

A combinatorial approach of candidate gene-based association analysis and genome-wide association study (GWAS) integrated with QTL mapping, differential gene expression profiling and molecular haplotyping was deployed in the present study for quantitative dissection of complex flowering time trait in chickpea. Candidate gene-based association mapping in a flowering time association panel (92 diverse desi and kabuli accessions) was performed by employing the genotyping information of 5724 SNPs discovered from 82 known flowering chickpea gene orthologs of Arabidopsis and legumes as well as 832 gene-encoding transcripts that are differentially expressed during flower development in chickpea. GWAS using both genome-wide GBS- and candidate gene-based genotyping data of 30,129 SNPs in a structured population of 92 sequenced accessions (with 200-250 kb LD decay) detected eight maximum effect genomic SNP loci (genes) associated (34% combined PVE) with flowering time. Six flowering time-associated major genomic loci harbouring five robust QTLs mapped on a high-resolution intra-specific genetic linkage map were validated (11.6-27.3% PVE at 5.4-11.7 LOD) further by traditional QTL mapping. The flower-specific expression, including differential up- and down-regulation (>three folds) of eight flowering time-associated genes (including six genes validated by QTL mapping) especially in early flowering than late flowering contrasting chickpea accessions/mapping individuals during flower development was evident. The gene haplotype-based LD mapping discovered diverse novel natural allelic variants and haplotypes in eight genes with high trait association potential (41% combined PVE) for flowering time differentiation in cultivated and wild chickpea. Taken together, eight potential known/candidate flowering time-regulating genes [efl1 (early flowering 1), FLD (Flowering locus D), GI (GIGANTEA), Myb (Myeloblastosis), SFH3 (SEC14-like 3), bZIP (basic-leucine zipper), bHLH (basic helix-loop-helix) and SBP (SQUAMOSA promoter binding protein)], including novel markers, QTLs, alleles and haplotypes delineated by aforesaid genome-wide integrated approach have potential for marker-assisted genetic improvement and unravelling the domestication pattern of flowering time in chickpea.

Development of gene-based markers for use in construction of the chickpea (Cicer arietinum L.) genetic linkage map and identification of QTLs associated with seed weight and plant height

Seed weight and plant height are important agronomic traits and contribute to seed yield. The objective of this study was to identify QTLs underlying these traits using an intra-specific mapping population of chickpea. A F11 population of 177 recombinant inbred lines derived from a cross between SBD377 (100-seed weight--48 g and plant height--53 cm) and BGD112 (100-seed weight--15 g and plant height--65 cm) was used. A total of 367 novel EST-derived functional markers were developed which included 187 EST-SSRs, 130 potential intron polymorphisms (PIPs) and 50 expressed sequence tag polymorphisms (ESTPs). Along with these, 590 previously published markers including 385 EST-based markers and 205 genomic SSRs were utilized. Of the 957 markers tested for analysis of parental polymorphism between the two parents of the mapping population, 135 (14.64%) were found to be polymorphic. Of these, 131 polymorphic markers could be mapped to the 8 linkage groups. The linkage map had a total length of 1140.54 cM with an average marker density of 8.7 cM. The map was further used for QTL identification using composite interval mapping method (CIM). Two QTLs each for seed weight, qSW-1 and qSW-2 (explaining 11.54 and 19.24% of phenotypic variance, respectively) and plant height, qPH-1 and qPH-2 (explaining 13.98 and 12.17% of phenotypic variance, respectively) were detected. The novel set of genic markers, the intra-specific linkage map and the QTLs identified in the present study will serve as valuable genomic resources in improving the chickpea seed yield using marker-assisted selection (MAS) strategies.

CTDB: An Integrated Chickpea Transcriptome Database for Functional and Applied Genomics

Chickpea is an important grain legume used as a rich source of protein in human diet. The narrow genetic diversity and limited availability of genomic resources are the major constraints in implementing breeding strategies and biotechnological interventions for genetic enhancement of chickpea. We developed an integrated Chickpea Transcriptome Database (CTDB), which provides the comprehensive web interface for visualization and easy retrieval of transcriptome data in chickpea. The database features many tools for similarity search, functional annotation (putative function, PFAM domain and gene ontology) search and comparative gene expression analysis. The current release of CTDB (v2.0) hosts transcriptome datasets with high quality functional annotation from cultivated (desi and kabuli types) and wild chickpea. A catalog of transcription factor families and their expression profiles in chickpea are available in the database. The gene expression data have been integrated to study the expression profiles of chickpea transcripts in major tissues/organs and various stages of flower development. The utilities, such as similarity search, ortholog identification and comparative gene expression have also been implemented in the database to facilitate comparative genomic studies among different legumes and Arabidopsis. Furthermore, the CTDB represents a resource for the discovery of functional molecular markers (microsatellites and single nucleotide polymorphisms) between different chickpea types. We anticipate that integrated information content of this database will accelerate the functional and applied genomic research for improvement of chickpea. The CTDB web service is freely available at http://nipgr.res.in/ctdb.html.

High density linkage mapping of genomic and transcriptomic SNPs for synteny analysis and anchoring the genome sequence of chickpea

This study presents genome-wide discovery of SNPs through next generation sequencing of the genome of Cicer reticulatum. Mapping of the C. reticulatum sequenced reads onto the draft genome assembly of C. arietinum (desi chickpea) resulted in identification of 842,104 genomic SNPs which were utilized along with an additional 36,446 genic SNPs identified from transcriptome sequences of the aforementioned varieties. Two new chickpea Oligo Pool All (OPAs) each having 3,072 SNPs were designed and utilized for SNP genotyping of 129 Recombinant Inbred Lines (RILs). Using Illumina GoldenGate Technology genotyping data of 5,041 SNPs were generated and combined with the 1,673 marker data from previously published studies, to generate a high resolution linkage map. The map comprised of 6698 markers distributed on eight linkage groups spanning 1083.93 cM with an average inter-marker distance of 0.16 cM. Utility of the present map was demonstrated for improving the anchoring of the earlier reported draft genome sequence of desi chickpea by ~30% and that of kabuli chickpea by 18%. The genetic map reported in this study represents the most dense linkage map of chickpea , with the potential to facilitate efficient anchoring of the draft genome sequences of desi as well as kabuli chickpea varieties.

Development of genome-wide informative simple sequence repeat markers for large-scale genotyping applications in chickpea and development of web resource

Development of informative polymorphic simple sequence repeat (SSR) markers at a genome-wide scale is essential for efficient large-scale genotyping applications. We identified genome-wide 1835 SSRs showing polymorphism between desi and kabuli chickpea. A total of 1470 polymorphic SSR markers from diverse coding and non-coding regions of the chickpea genome were developed. These physically mapped SSR markers exhibited robust amplification efficiency (73.9%) and high intra- and inter-specific polymorphic potential (63.5%), thereby suggesting their immense use in various genomics-assisted breeding applications. The SSR markers particularly derived from intergenic and intronic sequences revealed high polymorphic potential. Using the mapped SSR markers, a wider functional molecular diversity (16-94%, mean: 68%), and parentage- and cultivar-specific admixed domestication pattern and phylogenetic relationships in a structured population of desi and kabuli chickpea genotypes was evident. The intra-specific polymorphism (47.6%) and functional molecular diversity (65%) potential of polymorphic SSR markers developed in our study is much higher than that of previous documentations. Finally, we have developed a user-friendly web resource, Chickpea Microsatellite Database (CMsDB; http://www.nipgr.res.in/CMsDB.html), which provides public access to the data and results reported in this study. The developed informative SSR markers can serve as a resource for various genotyping applications, including genetic enhancement studies in chickpea.

An advanced draft genome assembly of a desi type chickpea (Cicer arietinum L.)

Chickpea (Cicer arietinum L.) is an important pulse legume crop. We previously reported a draft genome assembly of the desi chickpea cultivar ICC 4958. Here we report an advanced version of the ICC 4958 genome assembly (version 2.0) generated using additional sequence data and an improved genetic map. This resulted in 2.7-fold increase in the length of the pseudomolecules and substantial reduction of sequence gaps. The genome assembly covered more than 94% of the estimated gene space and predicted the presence of 30,257 protein-coding genes including 2230 and 133 genes encoding potential transcription factors (TF) and resistance gene homologs, respectively. Gene expression analysis identified several TF and chickpea-specific genes with tissue-specific expression and displayed functional diversification of the paralogous genes. Pairwise comparison of pseudomolecules in the desi (ICC 4958) and the earlier reported kabuli (CDC Frontier) chickpea assemblies showed an extensive local collinearity with incongruity in the placement of large sequence blocks along the linkage groups, apparently due to use of different genetic maps. Single nucleotide polymorphism (SNP)-based mining of intra-specific polymorphism identified more than four thousand SNPs differentiating a desi group and a kabuli group of chickpea genotypes.

Genome-wide high-throughput SNP discovery and genotyping for understanding natural (functional) allelic diversity and domestication patterns in wild chickpea

We identified 82489 high-quality genome-wide SNPs from 93 wild and cultivated Cicer accessions through integrated reference genome- and de novo-based GBS assays. High intra- and inter-specific polymorphic potential (66-85%) and broader natural allelic diversity (6-64%) detected by genome-wide SNPs among accessions signify their efficacy for monitoring introgression and transferring target trait-regulating genomic (gene) regions/allelic variants from wild to cultivated Cicer gene pools for genetic improvement. The population-specific assignment of wild Cicer accessions pertaining to the primary gene pool are more influenced by geographical origin/phenotypic characteristics than species/gene-pools of origination. The functional significance of allelic variants (non-synonymous and regulatory SNPs) scanned from transcription factors and stress-responsive genes in differentiating wild accessions (with potential known sources of yield-contributing and stress tolerance traits) from cultivated desi and kabuli accessions, fine-mapping/map-based cloning of QTLs and determination of LD patterns across wild and cultivated gene-pools are suitably elucidated. The correlation between phenotypic (agromorphological traits) and molecular diversity-based admixed domestication patterns within six structured populations of wild and cultivated accessions via genome-wide SNPs was apparent. This suggests utility of whole genome SNPs as a potential resource for identifying naturally selected trait-regulating genomic targets/functional allelic variants adaptive to diverse agroclimatic regions for genetic enhancement of cultivated gene-pools.

A genome-wide SNP scan accelerates trait-regulatory genomic loci identification in chickpea

We identified 44844 high-quality SNPs by sequencing 92 diverse chickpea accessions belonging to a seed and pod trait-specific association panel using reference genome- and de novo-based GBS (genotyping-by-sequencing) assays. A GWAS (genome-wide association study) in an association panel of 211, including the 92 sequenced accessions, identified 22 major genomic loci showing significant association (explaining 23-47% phenotypic variation) with pod and seed number/plant and 100-seed weight. Eighteen trait-regulatory major genomic loci underlying 13 robust QTLs were validated and mapped on an intra-specific genetic linkage map by QTL mapping. A combinatorial approach of GWAS, QTL mapping and gene haplotype-specific LD mapping and transcript profiling uncovered one superior haplotype and favourable natural allelic variants in the upstream regulatory region of a CesA-type cellulose synthase (Ca_Kabuli_CesA3) gene regulating high pod and seed number/plant (explaining 47% phenotypic variation) in chickpea. The up-regulation of this superior gene haplotype correlated with increased transcript expression of Ca_Kabuli_CesA3 gene in the pollen and pod of high pod/seed number accession, resulting in higher cellulose accumulation for normal pollen and pollen tube growth. A rapid combinatorial genome-wide SNP genotyping-based approach has potential to dissect complex quantitative agronomic traits and delineate trait-regulatory genomic loci (candidate genes) for genetic enhancement in crop plants, including chickpea.

Development and Integration of Genome-Wide Polymorphic Microsatellite Markers onto a Reference Linkage Map for Constructing a High-Density Genetic Map of Chickpea

The identification of informative in silico polymorphic genomic and genic microsatellite markers by comparing the genome and transcriptome sequences of crop genotypes is a rapid, cost-effective and non-laborious approach for large-scale marker validation and genotyping applications, including construction of high-density genetic maps. We designed 1494 markers, including 1016 genomic and 478 transcript-derived microsatellite markers showing in-silico fragment length polymorphism between two parental genotypes (Cicer arietinum ICC4958 and C. reticulatum PI489777) of an inter-specific reference mapping population. High amplification efficiency (87%), experimental validation success rate (81%) and polymorphic potential (55%) of these microsatellite markers suggest their effective use in various applications of chickpea genetics and breeding. Intra-specific polymorphic potential (48%) detected by microsatellite markers in 22 desi and kabuli chickpea genotypes was lower than inter-specific polymorphic potential (59%). An advanced, high-density, integrated and inter-specific chickpea genetic map (ICC4958 x PI489777) having 1697 map positions spanning 1061.16 cM with an average inter-marker distance of 0.625 cM was constructed by assigning 634 novel informative transcript-derived and genomic microsatellite markers on eight linkage groups (LGs) of our prior documented, 1063 marker-based genetic map. The constructed genome map identified 88, including four major (7–23 cM) longest high-resolution genomic regions on LGs 3, 5 and 8, where the maximum number of novel genomic and genic microsatellite markers were specifically clustered within 1 cM genetic distance. It was for the first time in chickpea that in silico FLP analysis at genome-wide level was carried out and such a large number of microsatellite markers were identified, experimentally validated and further used in genetic mapping. To best of our knowledge, in the presently constructed genetic map, we mapped highest number of new sequence-based robust microsatellite markers (634) which is an advancement over the previously documented (~300 markers) inter-specific genetic maps. This advanced high-density map will serve as a foundation for large-scale marker validation and genotyping applications, including identification and targeted mapping of trait-specific genes/QTLs (quantitative trait loci) with sub-optimal use of resources and labour in chickpea.

Ultra-high density intra-specific genetic linkage maps accelerate identification of functionally relevant molecular tags governing important agronomic traits in chickpea

We discovered 26785 and 16573 high-quality SNPs differentiating two parental genotypes of a RIL mapping population using reference desi and kabuli genome-based GBS assay. Of these, 3625 and 2177 SNPs have been integrated into eight desi and kabuli chromosomes, respectively in order to construct ultra-high density (0.20-0.37 cM) intra-specific chickpea genetic linkage maps. One of these constructed high-resolution genetic map has potential to identify 33 major genomic regions harbouring 35 robust QTLs (PVE: 17.9-39.7%) associated with three agronomic traits, which were mapped within <1 cM mean marker intervals on desi chromosomes. The extended LD (linkage disequilibrium) decay (~15 cM) in chromosomes of genetic maps have encouraged us to use a rapid integrated approach (comparative QTL mapping, QTL-region specific haplotype/LD-based trait association analysis, expression profiling and gene haplotype-based association mapping) rather than a traditional QTL map-based cloning method to narrow-down one major seed weight (SW) robust QTL region. It delineated favourable natural allelic variants and superior haplotype-containing one seed-specific candidate embryo defective gene regulating SW in chickpea. The ultra-high-resolution genetic maps, QTLs/genes and alleles/haplotypes-related genomic information generated and integrated strategy for rapid QTL/gene identification developed have potential to expedite genomics-assisted breeding applications in crop plants, including chickpea for their genetic enhancement.

Deploying QTL-seq for rapid delineation of a potential candidate gene underlying major trait-associated QTL in chickpea

A rapid high-resolution genome-wide strategy for molecular mapping of major QTL(s)/gene(s) regulating important agronomic traits is vital for in-depth dissection of complex quantitative traits and genetic enhancement in chickpea. The present study for the first time employed a NGS-based whole-genome QTL-seq strategy to identify one major genomic region harbouring a robust 100-seed weight QTL using an intra-specific 221 chickpea mapping population (desi cv. ICC 7184 × desi cv. ICC 15061). The QTL-seq-derived major SW QTL (CaqSW1.1) was further validated by single-nucleotide polymorphism (SNP) and simple sequence repeat (SSR) marker-based traditional QTL mapping (47.6% R² at higher LOD >19). This reflects the reliability and efficacy of QTL-seq as a strategy for rapid genome-wide scanning and fine mapping of major trait regulatory QTLs in chickpea. The use of QTL-seq and classical QTL mapping in combination narrowed down the 1.37 Mb (comprising 177 genes) major SW QTL (CaqSW1.1) region into a 35 kb genomic interval on desi chickpea chromosome 1 containing six genes. One coding SNP (G/A)-carrying constitutive photomorphogenic9 (COP9) signalosome complex subunit 8 (CSN8) gene of these exhibited seed-specific expression, including pronounced differential up-/down-regulation in low and high seed weight mapping parents and homozygous individuals during seed development. The coding SNP mined in this potential seed weight-governing candidate CSN8 gene was found to be present exclusively in all cultivated species/genotypes, but not in any wild species/genotypes of primary, secondary and tertiary gene pools. This indicates the effect of strong artificial and/or natural selection pressure on target SW locus during chickpea domestication. The proposed QTL-seq-driven integrated genome-wide strategy has potential to delineate major candidate gene(s) harbouring a robust trait regulatory QTL rapidly with optimal use of resources. This will further assist us to extrapolate the molecular mechanism underlying complex quantitative traits at a genome-wide scale leading to fast-paced marker-assisted genetic improvement in diverse crop plants, including chickpea.

Transcriptome sequencing of mung bean (Vigna radiate L.) genes and the identification of EST-SSR markers

Mung bean (Vigna radiate (L.) Wilczek) is an important traditional food legume crop, with high economic and nutritional value. It is widely grown in China and other Asian countries. Despite its importance, genomic information is currently unavailable for this crop plant species or some of its close relatives in the Vigna genus. In this study, more than 103 million high quality cDNA sequence reads were obtained from mung bean using Illumina paired-end sequencing technology. The processed reads were assembled into 48,693 unigenes with an average length of 874 bp. Of these unigenes, 25,820 (53.0%) and 23,235 (47.7%) showed significant similarity to proteins in the NCBI non-redundant protein and nucleotide sequence databases, respectively. Furthermore, 19,242 (39.5%) could be classified into gene ontology categories, 18,316 (37.6%) into Swiss-Prot categories and 10,918 (22.4%) into KOG database categories (E-value < 1.0E-5). A total of 6,585 (8.3%) were mapped onto 244 pathways using the Kyoto Encyclopedia of Genes and Genome (KEGG) pathway database. Among the unigenes, 10,053 sequences contained a unique simple sequence repeat (SSR), and 2,303 sequences contained more than one SSR together in the same expressed sequence tag (EST). A total of 13,134 EST-SSRs were identified as potential molecular markers, with mono-nucleotide A/T repeats being the most abundant motif class and G/C repeats being rare. In this SSR analysis, we found five main repeat motifs: AG/CT (30.8%), GAA/TTC (12.6%), AAAT/ATTT (6.8%), AAAAT/ATTTT (6.2%) and AAAAAT/ATTTTT (1.9%). A total of 200 SSR loci were randomly selected for validation by PCR amplification as EST-SSR markers. Of these, 66 marker primer pairs produced reproducible amplicons that were polymorphic among 31 mung bean accessions selected from diverse geographical locations. The large number of SSR-containing sequences found in this study will be valuable for the construction of a high-resolution genetic linkage maps, association or comparative mapping and genetic analyses of various Vigna species.

Employing genome-wide SNP discovery and genotyping strategy to extrapolate the natural allelic diversity and domestication patterns in chickpea

The genome-wide discovery and high-throughput genotyping of SNPs in chickpea natural germplasm lines is indispensable to extrapolate their natural allelic diversity, domestication, and linkage disequilibrium (LD) patterns leading to the genetic enhancement of this vital legume crop. We discovered 44,844 high-quality SNPs by sequencing of 93 diverse cultivated desi, kabuli, and wild chickpea accessions using reference genome- and de novo-based GBS (genotyping-by-sequencing) assays that were physically mapped across eight chromosomes of desi and kabuli. Of these, 22,542 SNPs were structurally annotated in different coding and non-coding sequence components of genes. Genes with 3296 non-synonymous and 269 regulatory SNPs could functionally differentiate accessions based on their contrasting agronomic traits. A high experimental validation success rate (92%) and reproducibility (100%) along with strong sensitivity (93-96%) and specificity (99%) of GBS-based SNPs was observed. This infers the robustness of GBS as a high-throughput assay for rapid large-scale mining and genotyping of genome-wide SNPs in chickpea with sub-optimal use of resources. With 23,798 genome-wide SNPs, a relatively high intra-specific polymorphic potential (49.5%) and broader molecular diversity (13-89%)/functional allelic diversity (18-77%) was apparent among 93 chickpea accessions, suggesting their tremendous applicability in rapid selection of desirable diverse accessions/inter-specific hybrids in chickpea crossbred varietal improvement program. The genome-wide SNPs revealed complex admixed domestication pattern, extensive LD estimates (0.54-0.68) and extended LD decay (400-500 kb) in a structured population inclusive of 93 accessions. These findings reflect the utility of our identified SNPs for subsequent genome-wide association study (GWAS) and selective sweep-based domestication trait dissection analysis to identify potential genomic loci (gene-associated targets) specifically regulating important complex quantitative agronomic traits in chickpea. The numerous informative genome-wide SNPs, natural allelic diversity-led domestication pattern, and LD-based information generated in our study have got multidimensional applicability with respect to chickpea genomics-assisted breeding.

A new set of ESTs from chickpea (Cicer arietinum L.) embryo reveals two novel F-box genes, CarF-box_PP2 and CarF-box_LysM, with potential roles in seed development

Considering the economic importance of chickpea (C. arietinum L.) seeds, it is important to understand the mechanisms underlying seed development for which a cDNA library was constructed from 6 day old chickpea embryos. A total of 8,186 ESTs were obtained from which 4,048 high quality ESTs were assembled into 1,480 unigenes that majorly encoded genes involved in various metabolic and regulatory pathways. Of these, 95 ESTs were found to be involved in ubiquitination related protein degradation pathways and 12 ESTs coded specifically for putative F-box proteins. Differential transcript accumulation of these putative F-box genes was observed in chickpea tissues as evidenced by quantitative real-time PCR. Further, to explore the role of F-box proteins in chickpea seed development, two F-box genes were selected for molecular characterization. These were named as CarF-box_PP2 and CarF-box_LysM depending on their C-terminal domains, PP2 and LysM, respectively. Their highly conserved structures led us to predict their target substrates. Subcellular localization experiment revealed that CarF-box_PP2 was localized in the cytoplasm and CarF-box_LysM was localized in the nucleus. We demonstrated their physical interactions with SKP1 protein, which validated that they function as F-box proteins in the formation of SCF complexes. Sequence analysis of their promoter regions revealed certain seed specific cis-acting elements that may be regulating their preferential transcript accumulation in the seed. Overall, the study helped in expanding the EST database of chickpea, which was further used to identify two novel F-box genes having a potential role in seed development.