Search

Genetic diversity studies provide important details on target trait availability and its variability, for the success of breeding programs. In this study, GBS approach was used to reveal a new structuration of genetic diversity and population structure of pigeonpea in Benin. We used a total of 688 high-quality Single Nucleotide Polymorphism markers for a total of 44 pigeonpea genotypes. The distribution of SNP markers on the 11 chromosomes ranged from 14 on chromosome 5 to 133 on chromosome 2. The Polymorphism Information Content and gene diversity values were 0.30 and 0.34 respectively. The analysis of population structure revealed four clear subpopulations. The Weighted Neighbor Joining tree agreed with structure analyses by grouping the 44 genotypes into four clusters. The PCoA revealed that genotypes from subpopulations 1, 2 and 3 intermixed among themselves. The Analysis of Molecular Variance showed 7% of the total variation among genotypes while the rest of variation (93%) was within genotypes from subpopulations indicating a high gene exchange (Nm = 7.13) and low genetic differentiation (PhiPT = 0.07) between subpopulations. Subpopulation 2 presented the highest mean values of number of different alleles (Na = 1.57), number of loci with private alleles (Pa = 0.11) and the percentage of polymorphic loci (P = 57.12%). We discuss our findings and demonstrate how the genetic diversity and the population structure of this specie can be used through the Genome Wide Association Studies and Marker-Assisted Selection to enhance genetic gain in pigeonpea breeding programs in Benin., Competing Interests: The authors have declared that no competing interests exist.

Heat stress during reproductive stages has been leading to significant yield losses in chickpea (Cicer arietinum L.). With an aim of identifying the genomic regions or QTLs responsible for heat tolerance, 187 F8 recombinant inbred lines (RILs) derived from the cross GPF 2 (heat tolerant) × ILWC 292 (heat sensitive) were evaluated under late-sown irrigated (January-May) and timely-sown irrigated environments (November-April) at Ludhiana and Faridkot in Punjab, India for 13 heat tolerance related traits. The pooled ANOVA for both locations for the traits namely days to germination (DG), days to flowering initiation (DFI), days to 50% flowering (DFF), days to 100% flowering (DHF), plant height (PH), pods per plant (NPP), biomass (BIO), grain yield (YLD), 100-seed weight (HSW), harvest index (HI), membrane permeability index (MPI), relative leaf water content (RLWC) and pollen viability (PV)) showed a highly significant difference in RILs. The phenotyping data coupled with the genetic map comprising of 1365 ddRAD-Seq based SNP markers were used for identifying the QTLs for heat tolerance. Composite interval mapping provided a total of 28 and 23 QTLs, respectively at Ludhiana and Faridkot locations. Of these, 13 consensus QTLs for DG, DFI, DFF, DHF, PH, YLD, and MPI have been identified at both locations. Four QTL clusters containing QTLs for multiple traits were identified on the same genomic region at both locations. Stable QTLs for days to flowering can be one of the major factors for providing heat tolerance as early flowering has an advantage of more seed setting due to a comparatively longer reproductive period. Identified QTLs can be used in genomics-assisted breeding to develop heat stress-tolerant high yielding chickpea cultivars., Competing Interests: The authors have declared that no competing interests exist.

Unravelling the genetic architecture underlying yield components and agronomic traits is important for enhancing crop productivity. Here, a recombinant inbred line (RIL) population, developed from ICC 4958 and DCP 92-3 cross, was used for constructing linkage map and QTL mapping analysis. The RIL population was genotyped using a high-throughput Axiom®CicerSNP array, which enabled the development of a high-density genetic map consisting of 3,818 SNP markers and spanning a distance of 1064.14 cM. Analysis of phenotyping data for yield, yield components and agronomic traits measured across three years together with genetic mapping data led to the identification of 10 major-effect QTLs and six minor-effect QTLs explaining up to 59.70% phenotypic variance. The major-effect QTLs identified for 100-seed weight, and plant height possessed key genes, such as C3HC4 RING finger protein, pentatricopeptide repeat (PPR) protein, sugar transporter, leucine zipper protein and NADH dehydrogenase, amongst others. The gene ontology studies highlighted the role of these genes in regulating seed weight and plant height in crop plants. The identified genomic regions for yield, yield components, and agronomic traits, and the closely linked markers will help advance genetics research and breeding programs in chickpea., Competing Interests: The authors have declared that no competing interests exist.

Photosynthates generated after heading contributes to 60% - 80% of grain yield in rice. Delay in leaf senescence can contribute to a long grain-filling period and thereby increased yield. The objective of this study was to identify genomic region(s) responsible for delayed leaf senescence (DLS) and validate the role of underlying candidate genes in controlling target traits. 302 BC2F4 backcross-derived lines (BILs) developed from a cross between Swarna and Moroberekan were phenotyped for two seasons (DS2016 and WS2017) for chlorophyll content and yield parameters. KASPar-SNP assays based genotyping data with 193 SNPs of mapping population was used to identify the targeted genomic region(s). Significant positive correlation was observed between the two most important determinants of DLS traits viz., RDCF (reduced decline degree of chlorophyll content of flag leaf) and RDCS (reduced decline degree of chlorophyll content of second leaf) with plant height (PH), grain number per panicle (GPN), panicle length (PL), number of tiller (NT) and grain yield (GY). A total of 41 and 29 QTLs with phenotypic variance (PVE) ranging from 8.2 to 25.1% were detected for six DLS traits during DS2016 and WS2017, respectively. Out of these identified QTLs, 19 were considered as stable QTLs detected across seasons. 17 of the identified stable QTLs were found to be novel. In-silico analysis revealed five key genes regulating chlorophyll metabolism. Expression analysis of these genes confirmed their strong association with the senescence pattern in leaf tissue of parents as well as selected phenotypically extreme lines. The identified stable QTLs regulating DLS traits and validation of potential candidate genes provides insight into genetic basis of delayed senescence and is expected to contribute in enhancing grain yield through genomics-assisted breeding (GAB)., Competing Interests: The authors have declared that no competing interests exist.

Peanut (Arachis hypogaea L.) is an important nutrient-rich food legume and valued for its good quality cooking oil. The fatty acid content is the major determinant of the quality of the edible oil. The oils containing higher monounsaturated fatty acid are preferred for improved shelf life and potential health benefits. Therefore, a high oleic/linoleic fatty acid ratio is the target trait in an advanced breeding program. The two mutant alleles, ahFAD2A (on linkage group a09) and ahFAD2B (on linkage group b09) control fatty acid composition for higher oleic/linoleic ratio in peanut. In the present study, marker-assisted backcrossing was employed for the introgression of two FAD2 mutant alleles from SunOleic95R into the chromosome of ICGV06100, a high oil content peanut breeding line. In the marker-assisted backcrossing-introgression lines, a 97% increase in oleic acid, and a 92% reduction in linoleic acid content was observed in comparison to the recurrent parent. Besides, the oleic/linoleic ratio was increased to 25 with respect to the recurrent parent, which was only 1.2. The most significant outcome was the stable expression of oil-content, oleic acid, linoleic acid, and palmitic acid in the marker-assisted backcrossing-introgression lines over the locations. No significant difference was observed between high oleic and normal oleic in peanuts for seedling traits except germination percentage. In addition, marker-assisted backcrossing-introgression lines exhibited higher yield and resistance to foliar fungal diseases, i.e., late leaf spot and rust., Competing Interests: All authors have contributed equally to this research work and have no competing interests exist.

Chickpea is one of the most important food legumes that holds the key to meet rising global food and nutritional demand. In order to deploy molecular breeding approaches in crop improvement programs, user friendly and cost effective marker resources remain prerequisite. The advent of next generation sequencing (NGS) technology has resulted in the generation of several thousands of markers as part of several large scale genome sequencing and re-sequencing initiatives. Very recently, PCR based Insertion-deletions (InDels) are becoming a popular gel based genotyping solution because of their co-dominant, inexpensive, and highly polymorphic nature. With an objective to expand marker resources for genomics assisted breeding (GAB) in chickpea, whole genome re-sequencing data generated on five parental lines of one interspecific (ICC 4958 × PI 489777) and two intra-specific (ICC 283 × ICC 8261 and ICC 4958 × ICC 1882) mapping populations, were used for identification of InDels. A total of 231,658 InDels were identified using Dindel software with default parameters. Further, a total of 8,307 InDels with ≥20 bp size were selected for development of gel based markers, of which primers could be designed for 7,523 (90.56%) markers. On average, markers appeared at a frequency of 1,038 InDels/LG with a maximum number of markers on CaLG04 (1,952 InDels) and minimum on CaLG08 (360 InDels). In order to validate these InDels, a total of 423 primer pairs were randomly selected and tested on the selected parental lines. A high amplification rate of 80% was observed ranging from 46.06 to 58.01% polymorphism rate across parents on 3% agarose gel. This study clearly reflects the usefulness of available sequence data for the development of genome-wide InDels in chickpea that can further contribute and accelerate a wide range of genetic and molecular breeding activities in chickpea., Competing Interests: The authors have declared that no competing interests exist.

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

Common bean (Phaseolus vulgaris L.) is one of the most important grain legume crops in the world. The beans grown in north-western Himalayas possess huge diversity for seed color, shape and size but are mostly susceptible to Anthracnose disease caused by seed born fungus Colletotrichum lindemuthianum. Dozens of QTLs/genes have been already identified for this disease in common bean world-wide. However, this is the first report of gene/QTL discovery for Anthracnose using bean germplasm from north-western Himalayas of state Jammu & Kashmir, India. A core set of 96 bean lines comprising 54 indigenous local landraces from 11 hot-spots and 42 exotic lines from 10 different countries were phenotyped at two locations (SKUAST-Jammu and Bhaderwah, Jammu) for Anthracnose resistance. The core set was also genotyped with genome-wide (91) random and trait linked SSR markers. The study of marker-trait associations (MTAs) led to the identification of 10 QTLs/genes for Anthracnose resistance. Among the 10 QTLs/genes identified, two MTAs are stable (BM45 & BM211), two MTAs (PVctt1 & BM211) are major explaining more than 20% phenotypic variation for Anthracnose and one MTA (BM211) is both stable and major. Six (06) genomic regions are reported for the first time, while as four (04) genomic regions validated the already known QTL/gene regions/clusters for Anthracnose. The major, stable and validated markers reported during the present study associated with Anthracnose resistance will prove useful in common bean molecular breeding programs aimed at enhancing Anthracnose resistance of local bean landraces grown in north-western Himalayas of state Jammu and Kashmir.

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

Peanut is vulnerable to a range of diseases, such as Tomato spotted wilt virus (TSWV) and leaf spots which will cause significant yield loss. The most sustainable, economical and eco-friendly solution for managing peanut diseases is development of improved cultivars with high level of resistance. We developed a recombinant inbred line population from the cross between SunOleic 97R and NC94022, named as the S-population. An improved genetic linkage map was developed for the S-population with 248 marker loci and a marker density of 5.7 cM/loci. This genetic map was also compared with the physical map of diploid progenitors of tetraploid peanut, resulting in an overall co-linearity of about 60% with the average co-linearity of 68% for the A sub-genome and 47% for the B sub-genome. The analysis using the improved genetic map and multi-season (2010-2013) phenotypic data resulted in the identification of 48 quantitative trait loci (QTLs) with phenotypic variance explained (PVE) from 3.88 to 29.14%. Of the 48 QTLs, six QTLs were identified for resistance to TSWV, 22 QTLs for early leaf spot (ELS) and 20 QTLs for late leaf spot (LLS), which included four, six, and six major QTLs (PVE larger than 10%) for each disease, respectively. A total of six major genomic regions (MGR) were found to have QTLs controlling more than one disease resistance. The identified QTLs and resistance gene-rich MGRs will facilitate further discovery of resistance genes and development of molecular markers for these important diseases.

Rapid popularity and adaptation of next generation sequencing (NGS) approaches have generated huge volumes of data. High throughput platforms like Illumina HiSeq produce terabytes of raw data that requires quick processing. Quality control of the data is an important component prior to the downstream analyses. To address these issues, we have developed a quality control pipeline, NGS-QCbox that scales up to process hundreds or thousands of samples. Raspberry is an in-house tool, developed in C language utilizing HTSlib (v1.2.1) (http://htslib.org), for computing read/base level statistics. It can be used as stand-alone application and can process both compressed and uncompressed FASTQ format files. NGS-QCbox integrates Raspberry with other open-source tools for alignment (Bowtie2), SNP calling (SAMtools) and other utilities (bedtools) towards analyzing raw NGS data at higher efficiency and in high-throughput manner. The pipeline implements batch processing of jobs using Bpipe (https://github.com/ssadedin/bpipe) in parallel and internally, a fine grained task parallelization utilizing OpenMP. It reports read and base statistics along with genome coverage and variants in a user friendly format. The pipeline developed presents a simple menu driven interface and can be used in either quick or complete mode. In addition, the pipeline in quick mode outperforms in speed against other similar existing QC pipeline/tools. The NGS-QCbox pipeline, Raspberry tool and associated scripts are made available at the URL https://github.com/CEG-ICRISAT/NGS-QCbox and https://github.com/CEG-ICRISAT/Raspberry for rapid quality control analysis of large-scale next generation sequencing (Illumina) data.

Gene expression analysis using quantitative real-time PCR (qRT-PCR) is a very sensitive technique and its sensitivity depends on the stable performance of reference gene(s) used in the study. A number of housekeeping genes have been used in various expression studies in many crops however, their expression were found to be inconsistent under different stress conditions. As a result, species specific housekeeping genes have been recommended for different expression studies in several crop species. However, such specific housekeeping genes have not been reported in the case of pigeonpea (Cajanus cajan) despite the fact that genome sequence has become available for the crop. To identify the stable housekeeping genes in pigeonpea for expression analysis under drought stress conditions, the relative expression variations of 10 commonly used housekeeping genes (EF1α, UBQ10, GAPDH, 18SrRNA, 25SrRNA, TUB6, ACT1, IF4α, UBC and HSP90) were studied on root, stem and leaves tissues of Asha (ICPL 87119). Three statistical algorithms geNorm, NormFinder and BestKeeper were used to define the stability of candidate genes. geNorm analysis identified IF4α and TUB6 as the most stable housekeeping genes however, NormFinder analysis determined IF4α and HSP90 as the most stable housekeeping genes under drought stress conditions. Subsequently validation of the identified candidate genes was undertaken in qRT-PCR based gene expression analysis of uspA gene which plays an important role for drought stress conditions in pigeonpea. The relative quantification of the uspA gene varied according to the internal controls (stable and least stable genes), thus highlighting the importance of the choice of as well as validation of internal controls in such experiments. The identified stable and validated housekeeping genes will facilitate gene expression studies in pigeonpea especially under drought stress conditions.

Peanut, a high-oil crop with about 50% oil content, is either crushed for oil or used as edible products. Fatty acid composition determines the oil quality which has high relevance to consumer health, flavor, and shelf life of commercial products. In addition to the major fatty acids, oleic acid (C18:1) and linoleic acid (C18:2) accounting for about 80% of peanut oil, the six other fatty acids namely palmitic acid (C16:0), stearic acid (C18:0), arachidic acid (C20:0), gadoleic acid (C20:1), behenic acid (C22:0), and lignoceric acid (C24:0) are accounted for the rest 20%. To determine the genetic basis and to improve further understanding on effect of FAD2 genes on these fatty acids, two recombinant inbred line (RIL) populations namely S-population (high oleic line 'SunOleic 97R' × low oleic line 'NC94022') and T-population (normal oleic line 'Tifrunner' × low oleic line 'GT-C20') were developed. Genetic maps with 206 and 378 marker loci for the S- and the T-population, respectively were used for quantitative trait locus (QTL) analysis. As a result, a total of 164 main-effect (M-QTLs) and 27 epistatic (E-QTLs) QTLs associated with the minor fatty acids were identified with 0.16% to 40.56% phenotypic variation explained (PVE). Thirty four major QTLs (>10% of PVE) mapped on five linkage groups and 28 clusters containing more than three QTLs were also identified. These results suggest that the major QTLs with large additive effects would play an important role in controlling composition of these minor fatty acids in addition to the oleic and linoleic acids in peanut oil. The interrelationship among these fatty acids should be considered while breeding for improved peanut genotypes with good oil quality and desired fatty acid composition.

Peanut is an important and nutritious agricultural commodity and a livelihood of many small-holder farmers in the semi-arid tropics (SAT) of world which are facing serious production threats. Integration of genomics tools with on-going genetic improvement approaches is expected to facilitate accelerated development of improved cultivars. Therefore, high-resolution genotyping and multiple season phenotyping data for 50 important agronomic, disease and quality traits were generated on the 'reference set' of peanut. This study reports comprehensive analyses of allelic diversity, population structure, linkage disequilibrium (LD) decay and marker-trait association (MTA) in peanut. Distinctness of all the genotypes can be established by using either an unique allele detected by a single SSR or a combination of unique alleles by two or more than two SSR markers. As expected, DArT features (2.0 alleles/locus, 0.125 PIC) showed lower allele frequency and polymorphic information content (PIC) than SSRs (22.21 alleles /locus, 0.715 PIC). Both marker types clearly differentiated the genotypes of diploids from tetraploids. Multi-allelic SSRs identified three sub-groups (K = 3) while the LD simulation trend line based on squared-allele frequency correlations (r2) predicted LD decay of 15-20 cM in peanut genome. Detailed analysis identified a total of 524 highly significant MTAs (p value > 2.1 × 10-6) with wide phenotypic variance (PV) range (5.81-90.09%) for 36 traits. These MTAs after validation may be deployed in improving biotic resistance, oil/ seed/ nutritional quality, drought tolerance related traits, and yield/ yield components.

To estimate genetic diversity within and between 10 interfertile Cicer species (94 genotypes) from the primary, secondary and tertiary gene pool, we analysed 5,257 DArT markers and 651 KASPar SNP markers. Based on successful allele calling in the tertiary gene pool, 2,763 DArT and 624 SNP markers that are polymorphic between genotypes from the gene pools were analyzed further. STRUCTURE analyses were consistent with 3 cultivated populations, representing kabuli, desi and pea-shaped seed types, with substantial admixture among these groups, while two wild populations were observed using DArT markers. AMOVA was used to partition variance among hierarchical sets of landraces and wild species at both the geographical and species level, with 61% of the variation found between species, and 39% within species. Molecular variance among the wild species was high (39%) compared to the variation present in cultivated material (10%). Observed heterozygosity was higher in wild species than the cultivated species for each linkage group. Our results support the Fertile Crescent both as the center of domestication and diversification of chickpea. The collection used in the present study covers all the three regions of historical chickpea cultivation, with the highest diversity in the Fertile Crescent region. Shared alleles between different gene pools suggest the possibility of gene flow among these species or incomplete lineage sorting and could indicate complicated patterns of divergence and fusion of wild chickpea taxa in the past.

Open source single nucleotide polymorphism (SNP) discovery pipelines for next generation sequencing data commonly requires working knowledge of command line interface, massive computational resources and expertise which is a daunting task for biologists. Further, the SNP information generated may not be readily used for downstream processes such as genotyping. Hence, a comprehensive pipeline has been developed by integrating several open source next generation sequencing (NGS) tools along with a graphical user interface called Integrated SNP Mining and Utilization (ISMU) for SNP discovery and their utilization by developing genotyping assays. The pipeline features functionalities such as pre-processing of raw data, integration of open source alignment tools (Bowtie2, BWA, Maq, NovoAlign and SOAP2), SNP prediction (SAMtools/SOAPsnp/CNS2snp and CbCC) methods and interfaces for developing genotyping assays. The pipeline outputs a list of high quality SNPs between all pairwise combinations of genotypes analyzed, in addition to the reference genome/sequence. Visualization tools (Tablet and Flapjack) integrated into the pipeline enable inspection of the alignment and errors, if any. The pipeline also provides a confidence score or polymorphism information content value with flanking sequences for identified SNPs in standard format required for developing marker genotyping (KASP and Golden Gate) assays. The pipeline enables users to process a range of NGS datasets such as whole genome re-sequencing, restriction site associated DNA sequencing and transcriptome sequencing data at a fast speed. The pipeline is very useful for plant genetics and breeding community with no computational expertise in order to discover SNPs and utilize in genomics, genetics and breeding studies. The pipeline has been parallelized to process huge datasets of next generation sequencing. It has been developed in Java language and is available at http://hpc.icrisat.cgiar.org/ISMU as a standalone free software.

To understand the genetic basis of tolerance to drought and heat stresses in chickpea, a comprehensive association mapping approach has been undertaken. Phenotypic data were generated on the reference set (300 accessions, including 211 mini-core collection accessions) for drought tolerance related root traits, heat tolerance, yield and yield component traits from 1-7 seasons and 1-3 locations in India (Patancheru, Kanpur, Bangalore) and three locations in Africa (Nairobi, Egerton in Kenya and Debre Zeit in Ethiopia). Diversity Array Technology (DArT) markers equally distributed across chickpea genome were used to determine population structure and three sub-populations were identified using admixture model in STRUCTURE. The pairwise linkage disequilibrium (LD) estimated using the squared-allele frequency correlations (r2; when r2<0.20) was found to decay rapidly with the genetic distance of 5 cM. For establishing marker-trait associations (MTAs), both genome-wide and candidate gene-sequencing based association mapping approaches were conducted using 1,872 markers (1,072 DArTs, 651 single nucleotide polymorphisms [SNPs], 113 gene-based SNPs and 36 simple sequence repeats [SSRs]) and phenotyping data mentioned above employing mixed linear model (MLM) analysis with optimum compression with P3D method and kinship matrix. As a result, 312 significant MTAs were identified and a maximum number of MTAs (70) was identified for 100-seed weight. A total of 18 SNPs from 5 genes (ERECTA, 11 SNPs; ASR, 4 SNPs; DREB, 1 SNP; CAP2 promoter, 1 SNP and AMDH, 1SNP) were significantly associated with different traits. This study provides significant MTAs for drought and heat tolerance in chickpea that can be used, after validation, in molecular breeding for developing superior varieties with enhanced drought and heat tolerance.

Understanding genetic structure of Cajanus spp. is essential for achieving genetic improvement by quantitative trait loci (QTL) mapping or association studies and use of selected markers through genomic assisted breeding and genomic selection. After developing a comprehensive set of 1,616 single nucleotide polymorphism (SNPs) and their conversion into cost effective KASPar assays for pigeonpea (Cajanus cajan), we studied levels of genetic variability both within and between diverse set of Cajanus lines including 56 breeding lines, 21 landraces and 107 accessions from 18 wild species. These results revealed a high frequency of polymorphic SNPs and relatively high level of cross-species transferability. Indeed, 75.8% of successful SNP assays revealed polymorphism, and more than 95% of these assays could be successfully transferred to related wild species. To show regional patterns of variation, we used STRUCTURE and Analysis of Molecular Variance (AMOVA) to partition variance among hierarchical sets of landraces and wild species at either the continental scale or within India. STRUCTURE separated most of the domesticated germplasm from wild ecotypes, and separates Australian and Asian wild species as has been found previously. Among Indian regions and states within regions, we found 36% of the variation between regions, and 64% within landraces or wilds within states. The highest level of polymorphism in wild relatives and landraces was found in Madhya Pradesh and Andhra Pradesh provinces of India representing the centre of origin and domestication of pigeonpea respectively.

A comprehensive transcriptome assembly of chickpea has been developed using 134.95 million Illumina single-end reads, 7.12 million single-end FLX/454 reads and 139,214 Sanger expressed sequence tags (ESTs) from >17 genotypes. This hybrid transcriptome assembly, referred to as Cicer arietinumTranscriptome Assembly version 2 (CaTA v2, available at http://data.comparative-legumes.org/transcriptomes/cicar/lista&#95;cicar-201201), comprising 46,369 transcript assembly contigs (TACs) has an N50 length of 1,726 bp and a maximum contig size of 15,644 bp. Putative functions were determined for 32,869 (70.8%) of the TACs and gene ontology assignments were determined for 21,471 (46.3%). The new transcriptome assembly was compared with the previously available chickpea transcriptome assemblies as well as to the chickpea genome. Comparative analysis of CaTA v2 against transcriptomes of three legumes - Medicago, soybean and common bean, resulted in 27,771 TACs common to all three legumes indicating strong conservation of genes across legumes. CaTA v2 was also used for identification of simple sequence repeats (SSRs) and intron spanning regions (ISRs) for developing molecular markers. ISRs were identified by aligning TACs to the Medicago genome, and their putative mapping positions at chromosomal level were identified using transcript map of chickpea. Primer pairs were designed for 4,990 ISRs, each representing a single contig for which predicted positions are inferred and distributed across eight linkage groups. A subset of randomly selected ISRs representing all eight chickpea linkage groups were validated on five chickpea genotypes and showed 20% polymorphism with average polymorphic information content (PIC) of 0.27. In summary, the hybrid transcriptome assembly developed and novel markers identified can be used for a variety of applications such as gene discovery, marker-trait association, diversity analysis etc., to advance genetics research and breeding applications in chickpea and other related legumes.

Small heat shock protein 17.8 (HSP17.8) is produced abundantly in plant cells under heat and other stress conditions and may play an important role in plant tolerance to stress environments. However, HSP17.8 may be differentially expressed in different accessions of a crop species exposed to identical stress conditions. The ability of different genotypes to adapt to various stress conditions resides in their genetic diversity. Allelic variations are the most common forms of genetic variation in natural populations. In this study, single nucleotide polymorphisms (SNPs) of the HSP17.8 gene were investigated across 210 barley accessions collected from 30 countries using EcoTILLING technology. Eleven SNPs including 10 from the coding region of HSP17.8 were detected, which form nine distinguishable haplotypes in the barley collection. Among the 10 SNPs in the coding region, six are missense mutations and four are synonymous nucleotide changes. Five of the six missense changes are predicted to be deleterious to HSP17.8 function. The accessions from Middle East Asia showed the higher nucleotide diversity of HSP17.8 than those from other regions and wild barley (H. spontaneum) accessions exhibited greater diversity than the cultivated barley (H. vulgare) accessions. Four SNPs in HSP17.8 were found associated with at least one of the agronomic traits evaluated except for spike length, namely number of grains per spike, thousand kernel weight, plant height, flag leaf area and leaf color. The association between SNP and these agronomic traits may provide new insight for study of the gene's potential contribution to drought tolerance of barley.

Only a few genetic maps based on recombinant inbred line (RIL) and backcross (BC) populations have been developed for tetraploid groundnut. The marker density, however, is not very satisfactory especially in the context of large genome size (2800 Mb/1C) and 20 linkage groups (LGs). Therefore, using marker segregation data for 10 RILs and one BC population from the international groundnut community, with the help of common markers across different populations, a reference consensus genetic map has been developed. This map is comprised of 897 marker loci including 895 simple sequence repeat (SSR) and 2 cleaved amplified polymorphic sequence (CAPS) loci distributed on 20 LGs (a01-a10 and b01-b10) spanning a map distance of 3, 863.6 cM with an average map density of 4.4 cM. The highest numbers of markers (70) were integrated on a01 and the least number of markers (21) on b09. The marker density, however, was lowest (6.4 cM) on a08 and highest (2.5 cM) on a01. The reference consensus map has been divided into 20 cM long 203 BINs. These BINs carry 1 (a10&#95;02, a10&#95;08 and a10&#95;09) to 20 (a10&#95;04) loci with an average of 4 marker loci per BIN. Although the polymorphism information content (PIC) value was available for 526 markers in 190 BINs, 36 and 111 BINs have at least one marker with >0.70 and >0.50 PIC values, respectively. This information will be useful for selecting highly informative and uniformly distributed markers for developing new genetic maps, background selection and diversity analysis. Most importantly, this reference consensus map will serve as a reliable reference for aligning new genetic and physical maps, performing QTL analysis in a multi-populations design, evaluating the genetic background effect on QTL expression, and serving other genetic and molecular breeding activities in groundnut.

Pigeonpea (Cajanus cajan) is an annual or short-lived perennial food legume of acute regional importance, providing significant protein to the human diet in less developed regions of Asia and Africa. Due to its narrow genetic base, pigeonpea improvement is increasingly reliant on introgression of valuable traits from wild forms, a practice that would benefit from knowledge of its domestication history and relationships to wild species. Here we use 752 single nucleotide polymorphisms (SNPs) derived from 670 low copy orthologous genes to clarify the evolutionary history of pigeonpea (79 accessions) and its wild relatives (31 accessions). We identified three well-supported lineages that are geographically clustered and congruent with previous nuclear and plastid sequence-based phylogenies. Among all species analyzed Cajanus cajanifolius is the most probable progenitor of cultivated pigeonpea. Multiple lines of evidence suggest recent gene flow between cultivated and non-cultivated forms, as well as historical gene flow between diverged but sympatric species. Evidence supports that primary domestication occurred in India, with a second and more recent nested population bottleneck focused in tropical regions that is the likely consequence of pigeonpea breeding. We find abundant allelic variation and genetic diversity among the wild relatives, with the exception of wild species from Australia for which we report a third bottleneck unrelated to domestication within India. Domesticated C. cajan possess 75% less allelic diversity than the progenitor clade of wild Indian species, indicating a severe "domestication bottleneck" during pigeonpea domestication.

Chickpea (Cicer arietinum L.) is the third most important cool season food legume, cultivated in arid and semi-arid regions of the world. The goal of this study was to develop novel molecular markers such as microsatellite or simple sequence repeat (SSR) markers from bacterial artificial chromosome (BAC)-end sequences (BESs) and diversity arrays technology (DArT) markers, and to construct a high-density genetic map based on recombinant inbred line (RIL) population ICC 4958 (C. arietinum)×PI 489777 (C. reticulatum). A BAC-library comprising 55,680 clones was constructed and 46,270 BESs were generated. Mining of these BESs provided 6,845 SSRs, and primer pairs were designed for 1,344 SSRs. In parallel, DArT arrays with ca. 15,000 clones were developed, and 5,397 clones were found polymorphic among 94 genotypes tested. Screening of newly developed BES-SSR markers and DArT arrays on the parental genotypes of the RIL mapping population showed polymorphism with 253 BES-SSR markers and 675 DArT markers. Segregation data obtained for these polymorphic markers and 494 markers data compiled from published reports or collaborators were used for constructing the genetic map. As a result, a comprehensive genetic map comprising 1,291 markers on eight linkage groups (LGs) spanning a total of 845.56 cM distance was developed (http://cmap.icrisat.ac.in/cmap/sm/cp/thudi/). The number of markers per linkage group ranged from 68 (LG 8) to 218 (LG 3) with an average inter-marker distance of 0.65 cM. While the developed resource of molecular markers will be useful for genetic diversity, genetic mapping and molecular breeding applications, the comprehensive genetic map with integrated BES-SSR markers will facilitate its anchoring to the physical map (under construction) to accelerate map-based cloning of genes in chickpea and comparative genome evolution studies in legumes.