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

Genome-wide discovery of di-nucleotide SSR markers based on whole genome re-sequencing data of Cicer arietinum L. and Cicer reticulatum Ladiz

Simple sequence repeats (SSRs) are valuable genetic markers due to their co-dominant inheritance, multi-allelic and reproducible nature. They have been largely used for exploiting genetic architecture of plant germplasms, phylogenetic analysis, and mapping studies. Among the SSRs, di-nucleotide repeats are the most frequent of the simple repeats distributed throughout the plant genomes. In present study, we aimed to discover and develop di-nucleotide SSR markers by using the whole genome re-sequencing (WGRS) data from Cicer arietinum L. and C. reticulatum Ladiz. A total of 35,329 InDels were obtained in C. arietinum, whereas 44,331 InDels in C. reticulatum. 3387 InDels with 2 bp length were detected in C. arietinum, there were 4704 in C. reticulatum. Among 8091 InDels, 58 di-nucleotide regions that were polymorphic between two species were selected and used for validation. We tested primers for evaluation of genetic diversity in 30 chickpea genotypes including C. arietinum, C. reticulatum, C. echinospermum P.H. Davis, C. anatolicum Alef., C. canariense A. Santos & G.P. Lewis, C. microphyllum Benth., C. multijugum Maesen, C. oxyodon Boiss. & Hohen. and C. songaricum Steph ex DC. A total of 244 alleles were obtained for 58 SSR markers giving an average of 2.36 alleles per locus. The observed heterozygosity was 0.08 while the expected heterozygosity was 0.345. Polymorphism information content was found to be 0.73 across all loci. Phylogenetic tree and principal coordinate analysis clearly divided the accessions into four groups. The SSR markers were also evaluated in 30 genotypes of a RIL population obtained from an interspecific cross between C. arietinum and C. reticulatum. Chi-square (χ²) test revealed an expected 1:1 segregation ratio in the population. These results demonstrated the success of SSR identification and marker development for chickpea with the use of WGRS data. The newly developed 58 SSR markers are expected to be useful for chickpea breeders.

The Landscape of Genome-Wide and Gender-Specific Microsatellites in Indo-Pacific Humpback Dolphin and Potential Applications in Cetacean Resource Investigation

Microsatellites are one of the important genome characterizations that can be a valuable resource for variety identification, genetic diversity, phylogenetic analysis, as well as comparative and conservation genomics research. Here, we developed comprehensive microsatellites through genome-wide mining for the threatened cetacean Indo-Pacific humpback dolphin (Sousa chinensis). We found 87,757 microsatellites with 2–6 bp nucleotide motifs, showing that about 32.5 microsatellites per megabase comprises microsatellites sequences. Approximately 97.8% of the markers developed in this study were consistent with the published identified markers. About 75.3% microsatellites were with dinucleotide motifs, followed by tetranucleotide motifs (17.4%), sharing the same composition pattern as other cetaceans. The microsatellites were not evenly distributed in the S. chinensis genome, mainly in non-coding regions, with only about 0.5% of the markers located in coding regions. The microsatellite-containing genes were mainly functionally enriched in the methylation process, probably demonstrating the potential impacts of microsatellites on biological functions. Polymorphic microsatellites were developed between different genders of S. chinensis, which was expected to lay the foundation for genetic diversity investigation in cetaceans. The specific markers for a male Indo-Pacific humpback dolphin will provide comprehensive and representative male candidate markers for sex identification, providing a potential biomolecular tool for further analysis of population structure and social behavior of wild populations, population trend evaluation, and species conservation management.

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.

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.

Genetic variation of a global germplasm collection of chickpea (Cicer arietinum L.) including Italian accessions at risk of genetic erosion

Chickpea (Cicer arietinum L.) is one of the most important legumes worldwide. We addressed this study to the genetic characterization of a germplasm collection from main chickpea growing countries. Several Italian traditional landraces at risk of genetic erosion were included in the analysis. Twenty-two simple sequence repeat (SSR) markers, widely used to explore genetic variation in plants, were selected and yielded 218 different alleles. Structure analysis and hierarchical clustering indicated that a model with three distinct subpopulations best fits the data. The composition of two subpopulations, named K1 and K2, broadly reflects the commercial classification of chickpea in the two types desi and kabuli, respectively. The third subpopulation (K3) is composed by both desi and kabuli genotypes. Italian accessions group both in K2 and K3. Interestingly, this study highlights genetic distance between desi genotypes cultivated in Asia and Ethiopia, which respectively represent the chickpea primary and the secondary centres of diversity. Moreover, European desi are closer to the Ethiopian gene pool. Overall, this study will be of importance for chickpea conservation genetics and breeding, which is limited by the poor characterization of germplasm collection.

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.

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.

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.