Search

Summary Single‐cell RNA‐seq (scRNA‐seq) has been highlighted as a powerful tool for the description of human cell transcriptome, but the technology has not been broadly applied in plant cells. Herein, we describe the successful development of a robust protoplast cell isolation system in the peanut leaf. A total of 6,815 single cells were divided into eight cell clusters based on reported marker genes by applying scRNA‐seq. Further, a pseudo‐time analysis was used to describe the developmental trajectory and interaction network of transcription factors (TFs) of distinct cell types during leaf growth. The trajectory enabled re‐investigation of the primordium‐driven development processes of the mesophyll and epidermis. These results suggest that palisade cells likely differentiate into spongy cells, while the epidermal cells originated earlier than the primordium. Subsequently, the developed method integrated multiple technologies to efficiently validate the scRNA‐seq result in a homogenous cell population. The expression levels of several TFs were strongly correlated with epidermal ontogeny in accordance with obtained scRNA‐seq values. Additionally, peanut AHL23 (AT‐HOOK MOTIF NUCLEAR LOCALIZED PROTEIN 23), which is localized in nucleus, promoted leaf growth when ectopically expressed in Arabidopsis by modulating the phytohormone pathway. Together, our study displays that application of scRNA‐seq can provide new hypotheses regarding cell differentiation in the leaf blade of Arachis hypogaea. We believe that this approach will enable significant advances in the functional study of leaf blade cells in the allotetraploid peanut and other plant species.

Celery (Apium graveolens L. 2n = 2x = 22), a member of the Apiaceae family, is among the most important and globally grown vegetables. Here, we report a high-quality genome sequence assembly, anchored to 11 chromosomes, with total length of 3.33 Gb and N50 scaffold length of 289.78 Mb. Most (92.91%) of the genome is composed of repetitive sequences, with 62.12% of 31 326 annotated genes confined to the terminal 20% of chromosomes. Simultaneous bursts of shared long-terminal repeats (LTRs) in different Apiaceae plants suggest inter-specific exchanges. Two ancestral polyploidizations were inferred, one shared by Apiales taxa and the other confined to Apiaceae. We reconstructed 8 Apiales proto-chromosomes, inferring their evolutionary trajectories from the eudicot common ancestor to extant plants. Transcriptome sequencing in three tissues (roots, leaves and petioles), and varieties with different-coloured petioles, revealed 4 and 2 key genes in pathways regulating anthocyanin and coumarin biosynthesis, respectively. A remarkable paucity of NBS disease-resistant genes in celery (62) and other Apiales was explained by extensive loss and limited production of these genes during the last ~10 million years, raising questions about their biotic defence mechanisms and motivating research into effects of chemicals, for example coumarins, that give off distinctive odours. Celery genome sequencing and annotation facilitates further research into important gene functions and breeding, and comparative genomic analyses in Apiales.

Haplotype-based breeding, a recent promising breeding approach to develop tailor-made crop varieties, deals with identification of superior haplotypes and their deployment in breeding programmes. In this context, whole genome re-sequencing data of 292 genotypes from pigeonpea reference set were mined to identify the superior haplotypes for 10 droughtresponsive candidate genes. A total of 83, 132 and 60 haplotypes were identified in breeding lines, landraces and wild species, respectively. Candidate gene-based association analysis of these 10 genes on a subset of 137 accessions of the pigeonpea reference set revealed 23 strong marker-trait associations (MTAs) in five genes influencing seven drought-responsive component traits. Haplo-pheno analysis for the strongly associated genes resulted in the identification of most promising haplotypes for three genes regulating five component drought traits. The haplotype C. cajan_23080-H2 for plant weight (PW), fresh weight (FW) and turgid weight (TW), the haplotype C. cajan_30211-H6 for PW, FW, TW and dry weight (DW), the haplotype C. cajan_26230-H11 for FW and DW and the haplotype C. cajan_26230-H5 for relative water content (RWC) were identified as superior haplotypes under drought stress condition. Furthermore, 17 accessions containing superior haplotypes for three drought-responsive genes were identified. The identified superior haplotypes and the accessions carrying these superior haplotypes will be very useful for deploying haplotype-based breeding to develop nextgeneration tailor-made better drought-responsive pigeonpea cultivars.

Summary Genotyping‐by‐sequencing has enabled approaches for genomic selection to improve yield, stress resistance and nutritional value. More and more resource studies are emerging providing 1000 and more genotypes and millions of SNPs for one species covering a hitherto inaccessible intraspecific genetic variation. The larger the databases are growing, the better statistical approaches for genomic selection will be available. However, there are clear limitations on the statistical but also on the biological part. Intraspecific genetic variation is able to explain a high proportion of the phenotypes, but a large part of phenotypic plasticity also stems from environmentally driven transcriptional, post‐transcriptional, translational, post‐translational, epigenetic and metabolic regulation. Moreover, regulation of the same gene can have different phenotypic outputs in different environments. Consequently, to explain and understand environment‐dependent phenotypic plasticity based on the available genotype variation we have to integrate the analysis of further molecular levels reflecting the complete information flow from the gene to metabolism to phenotype. Interestingly, metabolomics platforms are already more cost‐effective than NGS platforms and are decisive for the prediction of nutritional value or stress resistance. Here, we propose three fundamental pillars for future breeding strategies in the framework of Green Systems Biology: (i) combining genome selection with environment‐dependent PANOMICS analysis and deep learning to improve prediction accuracy for marker‐dependent trait performance; (ii) PANOMICS resolution at subtissue, cellular and subcellular level provides information about fundamental functions of selected markers; (iii) combining PANOMICS with genome editing and speed breeding tools to accelerate and enhance large‐scale functional validation of trait‐specific precision breeding.

Summary Spatio‐temporal and developmental stage‐specific transcriptome analysis plays a crucial role in systems biology‐based improvement of any species. In this context, we report here the Arachis hypogaea gene expression atlas (AhGEA) for the world's widest cultivated subsp. fastigiata based on RNA‐seq data using 20 diverse tissues across five key developmental stages. Approximately 480 million paired‐end filtered reads were generated followed by identification of 81 901 transcripts from an early‐maturing, high‐yielding, drought‐tolerant groundnut variety, ICGV 91114. Further, 57 344 genome‐wide transcripts were identified with ≥1 FPKM across different tissues and stages. Our in‐depth analysis of the global transcriptome sheds light into complex regulatory networks namely gravitropism and photomorphogenesis, seed development, allergens and oil biosynthesis in groundnut. Importantly, interesting insights into molecular basis of seed development and nodulation have immense potential for translational genomics research. We have also identified a set of stable expressing transcripts across the selected tissues, which could be utilized as internal controls in groundnut functional genomics studies. The AhGEA revealed potential transcripts associated with allergens, which upon appropriate validation could be deployed in the coming years to develop consumer‐friendly groundnut varieties. Taken together, the AhGEA touches upon various important and key features of cultivated groundnut and provides a reference for further functional, comparative and translational genomics research for various economically important traits.

Summary Pigeon pea (Cajanus cajan) is an important orphan crop mainly grown by smallholder farmers in India and Africa. Here, we present the first pigeon pea pangenome based on 89 accessions mainly from India and the Philippines, showing that there is significant genetic diversity in Philippine individuals that is not present in Indian individuals. Annotation of variable genes suggests that they are associated with self‐fertilization and response to disease. We identified 225 SNPs associated with nine agronomically important traits over three locations and two different time points, with SNPs associated with genes for transcription factors and kinases. These results will lead the way to an improved pigeon pea breeding programme.

Summary Hybrids are extensively used in agriculture to deliver an increase in yield, yet the molecular basis of heterosis is not well understood. Global DNA methylation analysis, transcriptome analysis and small RNA profiling were aimed to understand the epigenetic effect of the changes in gene expression level in the two hybrids and their parental lines. Increased DNA methylation was observed in both the hybrids as compared to their parents. This increased DNA methylation in hybrids showed that majority of the 24‐nt siRNA clusters had higher expression in hybrids than the parents. Transcriptome analysis revealed that various phytohormones (auxin and salicylic acid) responsive hybrid‐MPV DEGs were significantly altered in both the hybrids in comparison to MPV. DEGs associated with plant immunity and growth were overexpressed whereas DEGs associated with basal defence level were repressed. This antagonistic patterns of gene expression might contribute to the greater growth of the hybrids. It was also noticed that some common as well as unique changes in the regulatory pathways were associated with heterotic growth in both the hybrids. Approximately 70% and 67% of down‐regulated hybrid‐MPV DEGs were found to be differentially methylated in ICPH 2671 and ICPH 2740 hybrid, respectively. This reflected the association of epigenetic regulation in altered gene expressions. Our findings also revealed that miRNAs might play important roles in hybrid vigour in both the hybrids by regulating their target genes, especially in controlling plant growth and development, defence and stress response pathways. The above finding provides an insight into the molecular mechanism of pigeonpea heterosis.

We had the fortune of starting our scientific/research careers in the Molecular Biology and Crop Biotechnology Laboratory of Professor P.K. Gupta at Ch. Charan Singh University, Meerut, UP, India. Here, we describe the most important scientific contributions of our beloved mentor in the area of cytotaxonomy, cytogenetics, mutation breeding, quantitative genetics, molecular biology, crop biotechnology and plant genomics, on his 85th birthday. Important contributions made in the development and use of genomics resources including the development and use of different kinds of molecular markers, genetic and physical mapping, quantitative trait locus (QTL) interval mapping, genome-wide association mapping and molecular breeding including marker-assisted selection have been briefly summarized. Efforts have been also made to give readers a glimpse of important contributions in the study of cytology/apomixis of grasses, cytotaxonomic studies in asteraceae/fabaceae, nuclear/repetitive DNA content in Lolium, interspecific/intergeneric relationships involving the genus Hordeum and re-examining taxonomy of the tribe Triticeae.

Rice sheath blight disease, caused by the basidiomycetous necrotroph Rhizoctonia solani, became one of the major threats to the rice cultivation worldwide, especially after the adoption of high-yielding varieties. The pathogen is challenging to manage because of its extensively broad host range and high genetic variability and also due to the inability to find any satisfactory level of natural resistance from the available rice germplasm. It is high time to find remedies to combat the pathogen for reducing rice yield losses and subsequently to minimize the threat to global food security. The development of genetic resistance is one of the alternative means to avoid the use of hazardous chemical fungicides. This review mainly focuses on the effort of better understanding the host–pathogen relationship, finding the gene loci/markers imparting resistance response and modifying the host genome through transgenic development. The latest development and trend in the R. solani–rice pathosystem research with gap analysis are provided.

Summary Multiparental genetic mapping populations such as nested‐association mapping (NAM) have great potential for investigating quantitative traits and associated genomic regions leading to rapid discovery of candidate genes and markers. To demonstrate the utility and power of this approach, two NAM populations, NAM_Tifrunner and NAM_Florida‐07, were used for dissecting genetic control of 100‐pod weight (PW) and 100‐seed weight (SW) in peanut. Two high‐density SNP‐based genetic maps were constructed with 3341 loci and 2668 loci for NAM_Tifrunner and NAM_Florida‐07, respectively. The quantitative trait locus (QTL) analysis identified 12 and 8 major effect QTLs for PW and SW, respectively, in NAM_Tifrunner, and 13 and 11 major effect QTLs for PW and SW, respectively, in NAM_Florida‐07. Most of the QTLs associated with PW and SW were mapped on the chromosomes A05, A06, B05 and B06. A genomewide association study (GWAS) analysis identified 19 and 28 highly significant SNP–trait associations (STAs) in NAM_Tifrunner and 11 and 17 STAs in NAM_Florida‐07 for PW and SW, respectively. These significant STAs were co‐localized, suggesting that PW and SW are co‐regulated by several candidate genes identified on chromosomes A05, A06, B05, and B06. This study demonstrates the utility of NAM population for genetic dissection of complex traits and performing high‐resolution trait mapping in peanut.

Summary The transcriptome connects genome to the gene function and ultimate phenome in biology. So far, transcriptomic approach was not used in peanut for performing trait mapping in bi‐parental populations. In this research, we sequenced the whole transcriptome in immature seeds in a peanut recombinant inbred line (RIL) population and explored thoroughly the landscape of transcriptomic variations and its genetic basis. The comprehensive analysis identified total 49 691 genes in RIL population, of which 92 genes followed a paramutation‐like expression pattern. Expression quantitative trait locus (eQTL) analysis identified 1207 local eQTLs and 15 837 distant eQTLs contributing to the whole‐genome transcriptomic variation in peanut. There were 94 eQTL hot spot regions detected across the genome with the dominance of distant eQTL. By integrating transcriptomic profile and annotation analyses, we unveiled a putative candidate gene and developed a linked marker InDel02 underlying a major QTL responsible for purple testa colour in peanut. Our result provided a first understanding of genetic basis of whole‐genome transcriptomic variation in peanut and illustrates the potential of the transcriptome‐aid approach in dissecting important traits in non‐model plants.

Summary Silique number is the most important component of yield in rapeseed (Brassica napus L.). To dissect the mechanism underlying the natural variation of silique number in rapeseed germplasm, a series of studies were performed. A panel of 331 core lines was employed to genome‐wide association study (GWAS), and 27 loci (including 20 novel loci) were identified. The silique number difference between the more‐ and fewer‐silique lines can be attributed to the accumulative differences in flower number and silique setting rate. Each of them accounted for 75.2% and 24.8%, respectively. The silique number was highly associated with the total photosynthesis and biomass. Microscopic analysis showed that the difference between extremely more‐ and fewer‐silique lines normally occurred at the amount of flower bud but not morphology. Transcriptome analysis of shoot apical meristem (SAM) suggested that most of enriched groups were associated with the auxin biosynthesis/metabolism, vegetative growth and nutrition/energy accumulation. By integrating GWAS and RNA‐seq results, six promising candidate genes were identified, and some of them were related to biomass accumulation. In conclusion, the natural variation of silique number is largely affected by the biomass and nutrition accumulation, which essentially reflects the positive regulatory relationship between the source and sink. Our study provides a comprehensive and systematic explanation for natural variation of silique number in rapeseed, which provides a foundation for its improvement.

Summary Cultivated peanut (Arachis hypogaea L.) is an important grain legume providing high‐quality cooking oil, rich proteins and other nutrients. Shelling percentage (SP) is the 2nd most important agronomic trait after pod yield and this trait significantly affects the economic value of peanut in the market. Deployment of diagnostic markers through genomics‐assisted breeding (GAB) can accelerate the process of developing improved varieties with enhanced SP. In this context, we deployed the QTL‐seq approach to identify genomic regions and candidate genes controlling SP in a recombinant inbred line population (Yuanza 9102 × Xuzhou 68‐4). Four libraries (two parents and two extreme bulks) were constructed and sequenced, generating 456.89–790.32 million reads and achieving 91.85%–93.18% genome coverage and 14.04–21.37 mean read depth. Comprehensive analysis of two sets of data (Yuanza 9102/two bulks and Xuzhou 68‐4/two bulks) using the QTL‐seq pipeline resulted in discovery of two overlapped genomic regions (2.75 Mb on A09 and 1.1 Mb on B02). Nine candidate genes affected by 10 SNPs with non‐synonymous effects or in UTRs were identified in these regions for SP. Cost‐effective KASP (Kompetitive Allele‐Specific PCR) markers were developed for one SNP from A09 and three SNPs from B02 chromosome. Genotyping of the mapping population with these newly developed KASP markers confirmed the major control and stable expressions of these genomic regions across five environments. The identified candidate genomic regions and genes for SP further provide opportunity for gene cloning and deployment of diagnostic markers in molecular breeding for achieving high SP in improved varieties.

Summary Agricultural production faces a Herculean challenge to feed the increasing global population. Food production systems need to deliver more with finite land and water resources while exerting the least negative influence on the ecosystem. The unpredictability of climate change and consequent changes in pests/pathogens dynamics aggravate the enormity of the challenge. Crop improvement has made significant contributions towards food security, and breeding climate‐smart cultivars are considered the most sustainable way to accelerate food production. However, a fundamental change is needed in the conventional breeding framework in order to respond adequately to the growing food demands. Progress in genomics has provided new concepts and tools that hold promise to make plant breeding procedures more precise and efficient. For instance, reference genome assemblies in combination with germplasm sequencing delineate breeding targets that could contribute to securing future food supply. In this review, we highlight key breakthroughs in plant genome sequencing and explain how the presence of these genome resources in combination with gene editing techniques has revolutionized the procedures of trait discovery and manipulation. Adoption of new approaches such as speed breeding, genomic selection and haplotype‐based breeding could overcome several limitations of conventional breeding. We advocate that strengthening varietal release and seed distribution systems will play a more determining role in delivering genetic gains at farmer’s field. A holistic approach outlined here would be crucial to deliver steady stream of climate‐smart crop cultivars for sustainable agriculture.

Summary Jatropha curcas (physic nut), a non‐edible oilseed crop, represents one of the most promising alternative energy sources due to its high seed oil content, rapid growth and adaptability to various environments. We report ~339 Mbp draft whole genome sequence of J. curcas var. Chai Nat using both the PacBio and Illumina sequencing platforms. We identified and categorized differentially expressed genes related to biosynthesis of lipid and toxic compound among four stages of seed development. Triacylglycerol (TAG), the major component of seed storage oil, is mainly synthesized by phospholipid:diacylglycerol acyltransferase in Jatropha, and continuous high expression of homologs of oleosin over seed development contributes to accumulation of high level of oil in kernels by preventing the breakdown of TAG. A physical cluster of genes for diterpenoid biosynthetic enzymes, including casbene synthases highly responsible for a toxic compound, phorbol ester, in seed cake, was syntenically highly conserved between Jatropha and castor bean. Transcriptomic analysis of female and male flowers revealed the up‐regulation of a dozen family of TFs in female flower. Additionally, we constructed a robust species tree enabling estimation of divergence times among nine Jatropha species and five commercial crops in Malpighiales order. Our results will help researchers and breeders increase energy efficiency of this important oil seed crop by improving yield and oil content, and eliminating toxic compound in seed cake for animal feed.

Summary The cultivated [Glycine max (L) Merr.] and wild [Glycine soja Siebold & Zucc.] soybean species comprise wide variation in seed composition traits. Compared to wild soybean, cultivated soybean contains low protein, high oil, and high sucrose. In this study, an interspecific population was derived from a cross between G. max (Williams 82) and G. soja (PI 483460B). This recombinant inbred line (RIL) population of 188 lines was sequenced at 0.3× depth. Based on 91 342 single nucleotide polymorphisms (SNPs), recombination events in RILs were defined, and a high‐resolution bin map was developed (4070 bins). In addition to bin mapping, quantitative trait loci (QTL) analysis for protein, oil, and sucrose was performed using 3343 polymorphic SNPs (3K‐SNP), derived from Illumina Infinium BeadChip sequencing platform. The QTL regions from both platforms were compared, and a significant concordance was observed between bin and 3K‐SNP markers. Importantly, the bin map derived from next‐generation sequencing technology enhanced mapping resolution (from 1325 to 50 Kb). A total of five, nine, and four QTLs were identified for protein, oil, and sucrose content, respectively, and some of the QTLs coincided with soybean domestication‐related genomic loci. The major QTL for protein and oil were mapped on Chr. 20 (qPro_20) and suggested negative correlation between oil and protein. In terms of sucrose content, a novel and major QTL were identified on Chr. 8 (qSuc_08) and harbours putative genes involved in sugar transport. In addition, genome‐wide association using 91 342 SNPs confirmed the genomic loci derived from QTL mapping. A QTL‐based haplotype using whole‐genome resequencing of 106 diverse soybean lines identified unique allelic variation in wild soybean that could be utilized to widen the genetic base in cultivated soybean.

Summary Whole‐genome resequencing (WGRS) of mapping populations has facilitated development of high‐density genetic maps essential for fine mapping and candidate gene discovery for traits of interest in crop species. Leaf spots, including early leaf spot (ELS) and late leaf spot (LLS), and Tomato spotted wilt virus (TSWV) are devastating diseases in peanut causing significant yield loss. We generated WGRS data on a recombinant inbred line population, developed a SNP‐based high‐density genetic map, and conducted fine mapping, candidate gene discovery and marker validation for ELS, LLS and TSWV. The first sequence‐based high‐density map was constructed with 8869 SNPs assigned to 20 linkage groups, representing 20 chromosomes, for the ‘T’ population (Tifrunner × GT‐C20) with a map length of 3120 cM and an average distance of 1.45 cM. The quantitative trait locus (QTL) analysis using high‐density genetic map and multiple season phenotyping data identified 35 main‐effect QTLs with phenotypic variation explained (PVE) from 6.32% to 47.63%. Among major‐effect QTLs mapped, there were two QTLs for ELS on B05 with 47.42% PVE and B03 with 47.38% PVE, two QTLs for LLS on A05 with 47.63% and B03 with 34.03% PVE and one QTL for TSWV on B09 with 40.71% PVE. The epistasis and environment interaction analyses identified significant environmental effects on these traits. The identified QTL regions had disease resistance genes including R‐genes and transcription factors. KASP markers were developed for major QTLs and validated in the population and are ready for further deployment in genomics‐assisted breeding in peanut.

Identification of candidate genomic regions associated with target traits using conventional mapping methods is challenging and time‐consuming. In recent years, a number of single nucleotide polymorphism (SNP)‐based mapping approaches have been developed and used for identification of candidate/putative genomic regions. However, in the majority of these studies, insertion–deletion (Indel) were largely ignored. For efficient use of Indels in mapping target traits, we propose Indel‐seq approach, which is a combination of whole‐genome resequencing (WGRS) and bulked segregant analysis (BSA) and relies on the Indel frequencies in extreme bulks. Deployment of Indel‐seq approach for identification of candidate genomic regions associated with fusarium wilt (FW) and sterility mosaic disease (SMD) resistance in pigeonpea has identified 16 Indels affecting 26 putative candidate genes. Of these 26 affected putative candidate genes, 24 genes showed effect in the upstream/downstream of the genic region and two genes showed effect in the genes. Validation of these 16 candidate Indels in other FW‐ and SMD‐resistant and FW‐ and SMD‐susceptible genotypes revealed a significant association of five Indels (three for FW and two for SMD resistance). Comparative analysis of Indel‐seq with other genetic mapping approaches highlighted the importance of the approach in identification of significant genomic regions associated with target traits. Therefore, the Indel‐seq approach can be used for quick and precise identification of candidate genomic regions for any target traits in any crop species.

Summary Rust and late leaf spot (LLS) are the two major foliar fungal diseases in groundnut, and their co‐occurrence leads to significant yield loss in addition to the deterioration of fodder quality. To identify candidate genomic regions controlling resistance to rust and LLS, whole‐genome resequencing (WGRS)‐based approach referred as ‘QTL‐seq’ was deployed. A total of 231.67 Gb raw and 192.10 Gb of clean sequence data were generated through WGRS of resistant parent and the resistant and susceptible bulks for rust and LLS. Sequence analysis of bulks for rust and LLS with reference‐guided resistant parent assembly identified 3136 single‐nucleotide polymorphisms (SNPs) for rust and 66 SNPs for LLS with the read depth of ≥7 in the identified genomic region on pseudomolecule A03. Detailed analysis identified 30 nonsynonymous SNPs affecting 25 candidate genes for rust resistance, while 14 intronic and three synonymous SNPs affecting nine candidate genes for LLS resistance. Subsequently, allele‐specific diagnostic markers were identified for three SNPs for rust resistance and one SNP for LLS resistance. Genotyping of one RIL population (TAG 24 × GPBD 4) with these four diagnostic markers revealed higher phenotypic variation for these two diseases. These results suggest usefulness of QTL‐seq approach in precise and rapid identification of candidate genomic regions and development of diagnostic markers for breeding applications.

Summary Terminal drought is a major constraint to chickpea productivity. Two component traits responsible for reduction in yield under drought stress include reduction in seeds size and root length/root density. QTL‐seq approach, therefore, was used to identify candidate genomic regions for 100‐seed weight (100SDW) and total dry root weight to total plant dry weight ratio (RTR) under rainfed conditions. Genomewide SNP profiling of extreme phenotypic bulks from the ICC 4958 × ICC 1882 population identified two significant genomic regions, one on CaLG01 (1.08 Mb) and another on CaLG04 (2.7 Mb) linkage groups for 100SDW. Similarly, one significant genomic region on CaLG04 (1.10 Mb) was identified for RTR. Comprehensive analysis revealed four and five putative candidate genes associated with 100SDW and RTR, respectively. Subsequently, two genes (Ca_04364 and Ca_04607) for 100SDW and one gene (Ca_04586) for RTR were validated using CAPS/dCAPS markers. Identified candidate genomic regions and genes may be useful for molecular breeding for chickpea improvement.

Summary Bacterial wilt, caused by Ralstonia solanacearum, is a devastating disease affecting over 350 plant species. A few peanut cultivars were found to possess stable and durable bacterial wilt resistance (BWR). Genomics‐assisted breeding can accelerate the process of developing resistant cultivars by using diagnostic markers. Here, we deployed sequencing‐based trait mapping approach, QTL‐seq, to discover genomic regions, candidate genes and diagnostic markers for BWR in a recombination inbred line population (195 progenies) of peanut. The QTL‐seq analysis identified one candidate genomic region on chromosome B02 significantly associated with BWR. Mapping of newly developed single nucleotide polymorphism (SNP) markers narrowed down the region to 2.07 Mb and confirmed its major effects and stable expressions across three environments. This candidate genomic region had 49 nonsynonymous SNPs affecting 19 putative candidate genes including seven putative resistance genes (R‐genes). Two diagnostic markers were successfully validated in diverse breeding lines and cultivars and could be deployed in genomics‐assisted breeding of varieties with enhanced BWR.

Summary Whole genome re‐sequencing (WGRS) was conducted on a panel of 20 Cajanus spp. accessions (crossing parentals of recombinant inbred lines, introgression lines, multiparent advanced generation intercross and nested association mapping population) comprising of two wild species and 18 cultivated species accessions. A total of 791.77 million paired‐end reads were generated with an effective mapping depth of ~12X per accession. Analysis of WGRS data provided 5 465 676 genome‐wide variations including 4 686 422 SNPs and 779 254 InDels across the accessions. Large structural variations in the form of copy number variations (2598) and presence and absence variations (970) were also identified. Additionally, 2 630 904 accession‐specific variations comprising of 2 278 571 SNPs (86.6%), 166 243 deletions (6.3%) and 186 090 insertions (7.1%) were also reported. Identified polymorphic sites in this study provide the first‐generation HapMap in Cajanus spp. which will be useful in mapping the genomic regions responsible for important traits.

Summary APETALA2/ethylene response factor (AP2/ERF) and heat‐shock protein 90 (HSP90) are two significant classes of transcription factor and molecular chaperone proteins which are known to be implicated under abiotic and biotic stresses. Comprehensive survey identified a total of 147 AP2/ERF genes in chickpea, 176 in pigeonpea, 131 in Medicago, 179 in common bean and 140 in Lotus, whereas the number of HSP90 genes ranged from 5 to 7 in five legumes. Sequence alignment and phylogenetic analyses distinguished AP2, ERF, DREB, RAV and soloist proteins, while HSP90 proteins segregated on the basis of their cellular localization. Deeper insights into the gene structure allowed ERF proteins to be classified into AP2s based on DNA‐binding domains, intron arrangements and phylogenetic grouping. RNA‐seq and quantitative real‐time PCR (qRT‐PCR) analyses in heat‐stressed chickpea as well as Fusarium wilt (FW)‐ and sterility mosaic disease (SMD)‐stressed pigeonpea provided insights into the modus operandi of AP2/ERF and HSP90 genes. This study identified potential candidate genes in response to heat stress in chickpea while for FW and SMD stresses in pigeonpea. For instance, two DREB genes (Ca_02170 and Ca_16631) and three HSP90 genes (Ca_23016, Ca_09743 and Ca_25602) in chickpea can be targeted as potential candidate genes. Similarly, in pigeonpea, a HSP90 gene, C.cajan_27949, was highly responsive to SMD in the resistant genotype ICPL 20096, can be recommended for further functional validation. Also, two DREB genes, C.cajan_41905 and C.cajan_41951, were identified as leads for further investigation in response to FW stress in pigeonpea.

A characteristic feature of peanut is the subterranean fructification, geocarpy, in which the gynophore ('peg'), a specialized organ that transitions from upward growth habit to downward outgrowth upon fertilization, drives the developing pod into the soil for subsequent development underground. As a step towards understanding this phenomenon, we explore the developmental dynamics of the peanut pod transcriptome at 11 successive stages. We identified 110 217 transcripts across developmental stages and quantified their abundance along a pod developmental gradient in pod wall. We found that the majority of transcripts were differentially expressed along the developmental gradient as well as identified temporal programs of gene expression, including hundreds of transcription factors. Thought to be an adaptation to particularly harsh subterranean environments, both up- and down-regulated gene sets in pod wall were enriched for response to a broad array of stimuli, like gravity, light and subterranean environmental factors. We also identified hundreds of transcripts associated with gravitropism and photomorphogenesis, which may be involved in the geocarpy. Collectively, this study forms a transcriptional baseline for geocarpy in peanut as well as provides a considerable body of evidence that transcriptional regulation in peanut aerial and subterranean fruits is complex.

Summary To map resistance genes for Fusarium wilt (FW) and sterility mosaic disease (SMD) in pigeonpea, sequencing‐based bulked segregant analysis (Seq‐BSA) was used. Resistant (R) and susceptible (S) bulks from the extreme recombinant inbred lines of ICPL 20096 × ICPL 332 were sequenced. Subsequently, SNP index was calculated between R‐ and S‐bulks with the help of draft genome sequence and reference‐guided assembly of ICPL 20096 (resistant parent). Seq‐BSA has provided seven candidate SNPs for FW and SMD resistance in pigeonpea. In parallel, four additional genotypes were re‐sequenced and their combined analysis with R‐ and S‐bulks has provided a total of 8362 nonsynonymous (ns) SNPs. Of 8362 nsSNPs, 60 were found within the 2‐Mb flanking regions of seven candidate SNPs identified through Seq‐BSA. Haplotype analysis narrowed down to eight nsSNPs in seven genes. These eight nsSNPs were further validated by re‐sequencing 11 genotypes that are resistant and susceptible to FW and SMD. This analysis revealed association of four candidate nsSNPs in four genes with FW resistance and four candidate nsSNPs in three genes with SMD resistance. Further, In silico protein analysis and expression profiling identified two most promising candidate genes namely C.cajan_01839 for SMD resistance and C.cajan_03203 for FW resistance. Identified candidate genomic regions/SNPs will be useful for genomics‐assisted breeding in pigeonpea.

Calcium is a universal signal in the regulation of wide aspects in biology, but few are known about the function of calcium in the control of early embryo development. Ca(2+) deficiency in soil induces early embryo abortion in peanut, producing empty pods, which is a general problem; however, the underlying mechanism remains unclear. In this study, embryo abortion was characterized to be caused by apoptosis marked with cell wall degradation. Using a method of SSH cDNA libraries associated with library lift (SSHaLL), 62 differentially expressed genes were isolated from young peanut embryos. These genes were classified to be stress responses, catabolic process, carbohydrate and lipid metabolism, embryo morphogenesis, regulation, etc. The cell retardation with cell wall degradation was caused by up-regulated cell wall hydrolases and down-regulated cellular synthases genes. HsfA4a, which was characterized to be important to embryo development, was significantly down-regulated under Ca(2+) -deficient conditions from 15 days after pegging (DAP) to 30 DAP. Two AhCYP707A4 genes, encoding abscisic acid (ABA) 8'-hydroxylases, key enzymes for ABA catabolism, were up-regulated by 21-fold under Ca(2+) -deficient conditions upstream of HsfA4a, reducing the ABA level in early embryos. Over-expression of AhCYP707A4 in Nicotiana benthamiana showed a phenotype of low ABA content with high numbers of aborted embryos, small pods and less seeds, which confirms that AhCYP707A4 is a key player in regulation of Ca(2+) deficiency-induced embryo abortion via ABA-mediated apoptosis. The results elucidated the mechanism of low Ca(2+) -induced embryo abortion and described the method for other fields of study.

Summary To accelerate genomics research and molecular breeding applications in chickpea, a high‐throughput SNP genotyping platform ‘Axiom® CicerSNP Array’ has been designed, developed and validated. Screening of whole‐genome resequencing data from 429 chickpea lines identified 4.9 million SNPs, from which a subset of 70 463 high‐quality nonredundant SNPs was selected using different stringent filter criteria. This was further narrowed down to 61 174 SNPs based on p‐convert score ≥0.3, of which 50 590 SNPs could be tiled on array. Among these tiled SNPs, a total of 11 245 SNPs (22.23%) were from the coding regions of 3673 different genes. The developed Axiom® CicerSNP Array was used for genotyping two recombinant inbred line populations, namely ICCRIL03 (ICC 4958 × ICC 1882) and ICCRIL04 (ICC 283 × ICC 8261). Genotyping data reflected high success and polymorphic rate, with 15 140 (29.93%; ICCRIL03) and 20 018 (39.57%; ICCRIL04) polymorphic SNPs. High‐density genetic maps comprising 13 679 SNPs spanning 1033.67 cM and 7769 SNPs spanning 1076.35 cM were developed for ICCRIL03 and ICCRIL04 populations, respectively. QTL analysis using multilocation, multiseason phenotyping data on these RILs identified 70 (ICCRIL03) and 120 (ICCRIL04) main‐effect QTLs on genetic map. Higher precision and potential of this array is expected to advance chickpea genetics and breeding applications.

A set of 2486 single nucleotide polymorphisms (SNPs) were compiled in chickpea using four approaches, namely (i) Solexa/Illumina sequencing (1409), (ii) amplicon sequencing of tentative orthologous genes (TOGs) (604), (iii) mining of expressed sequence tags (ESTs) (286) and (iv) sequencing of candidate genes (187). Conversion of these SNPs to the cost-effective and flexible throughput Competitive Allele Specific PCR (KASPar) assays generated successful assays for 2005 SNPs. These marker assays have been designated as Chickpea KASPar Assay Markers (CKAMs). Screening of 70 genotypes including 58 diverse chickpea accessions and 12 BC(3) F(2) lines showed 1341 CKAMs as being polymorphic. Genetic analysis of these data clustered chickpea accessions based on geographical origin. Genotyping data generated for 671 CKAMs on the reference mapping population (Cicer arietinum ICC 4958 × Cicer reticulatum PI 489777) were compiled with 317 unpublished TOG-SNPs and 396 published markers for developing the genetic map. As a result, a second-generation genetic map comprising 1328 marker loci including novel 625 CKAMs, 314 TOG-SNPs and 389 published marker loci with an average inter-marker distance of 0.59 cM was constructed. Detailed analyses of 1064 mapped loci of this second-generation chickpea genetic map showed a higher degree of synteny with genome of Medicago truncatula, followed by Glycine max, Lotus japonicus and least with Vigna unguiculata. Development of these cost-effective CKAMs for SNP genotyping will be useful not only for genetics research and breeding applications in chickpea, but also for utilizing genome information from other sequenced or model legumes.

With the expansion of next-generation sequencing technology and advanced bioinformatics, there has been a rapid growth of genome sequencing projects. However, while this technology enables the rapid and cost-effective assembly of draft genomes, the quality of these assemblies usually falls short of gold standard genome assemblies produced using the more traditional BAC by BAC and Sanger sequencing approaches. Assembly validation is often performed by the physical anchoring of genetically mapped markers, but this is prone to errors and the resolution is usually low, especially towards centromeric regions where recombination is limited. New approaches are required to validate reference genome assemblies. The ability to isolate individual chromosomes combined with next-generation sequencing permits the validation of genome assemblies at the chromosome level. We demonstrate this approach by the assessment of the recently published chickpea kabuli and desi genomes. While previous genetic analysis suggests that these genomes should be very similar, a comparison of their chromosome sizes and published assemblies highlights significant differences. Our chromosomal genomics analysis highlights short defined regions that appear to have been misassembled in the kabuli genome and identifies large-scale misassembly in the draft desi genome. The integration of chromosomal genomics tools within genome sequencing projects has the potential to significantly improve the construction and validation of genome assemblies. The approach could be applied both for new genome assemblies as well as published assemblies, and complements currently applied genome assembly strategies.

Chickpea (Cicer arietinum L.) is an important legume crop in the semi-arid regions of Asia and Africa. Gains in crop productivity have been low however, particularly because of biotic and abiotic stresses. To help enhance crop productivity using molecular breeding techniques, next generation sequencing technologies such as Roche/454 and Illumina/Solexa were used to determine the sequence of most gene transcripts and to identify drought-responsive genes and gene-based molecular markers. A total of 103 215 tentative unique sequences (TUSs) have been produced from 435 018 Roche/454 reads and 21 491 Sanger expressed sequence tags (ESTs). Putative functions were determined for 49 437 (47.8%) of the TUSs, and gene ontology assignments were determined for 20 634 (41.7%) of the TUSs. Comparison of the chickpea TUSs with the Medicago truncatula genome assembly (Mt 3.5.1 build) resulted in 42 141 aligned TUSs with putative gene structures (including 39 281 predicted intron/splice junctions). Alignment of ∼37 million Illumina/Solexa tags generated from drought-challenged root tissues of two chickpea genotypes against the TUSs identified 44 639 differentially expressed TUSs. The TUSs were also used to identify a diverse set of markers, including 728 simple sequence repeats (SSRs), 495 single nucleotide polymorphisms (SNPs), 387 conserved orthologous sequence (COS) markers, and 2088 intron-spanning region (ISR) markers. This resource will be useful for basic and applied research for genome analysis and crop improvement in chickpea.