Annapurna Chitikineni

Draft genome of the peanut A-genome progenitor (Arachis duranensis) provides insights into geocarpy, oil biosynthesis, and allergens

Significance We present a draft genome of the peanut A-genome progenitor, Arachis duranensis, providing details on total genes present in the genome. Genome analysis suggests that the peanut lineage was affected by at least three polyploidizations since the origin of eudicots. Resequencing of synthetic Arachis tetraploids reveals extensive gene conversion since their formation by human hands. The A. duranensis genome provides a major source of candidate genes for fructification, oil biosynthesis, and allergens, expanding knowledge of understudied areas of plant biology and human health impacts of plants. This study also provides millions of structural variations that can be used as genetic markers for the development of improved peanut varieties through genomics-assisted breeding. Peanut or groundnut (Arachis hypogaea L.), a legume of South American origin, has high seed oil content (45–56%) and is a staple crop in semiarid tropical and subtropical regions, partially because of drought tolerance conferred by its geocarpic reproductive strategy. We present a draft genome of the peanut A-genome progenitor, Arachis duranensis, and 50,324 protein-coding gene models. Patterns of gene duplication suggest the peanut lineage has been affected by at least three polyploidizations since the origin of eudicots. Resequencing of synthetic Arachis tetraploids reveals extensive gene conversion in only three seed-to-seed generations since their formation by human hands, indicating that this process begins virtually immediately following polyploid formation. Expansion of some specific gene families suggests roles in the unusual subterranean fructification of Arachis. For example, the S1Fa-like transcription factor family has 126 Arachis members, in contrast to no more than five members in other examined plant species, and is more highly expressed in roots and etiolated seedlings than green leaves. The A. duranensis genome provides a major source of candidate genes for fructification, oil biosynthesis, and allergens, expanding knowledge of understudied areas of plant biology and human health impacts of plants, informing peanut genetic improvement and aiding deeper sequencing of Arachis diversity.

Draft genome sequence of adzuki bean, Vigna angularis

Adzuki bean (Vigna angularis var. angularis) is a dietary legume crop in East Asia. The presumed progenitor (Vigna angularis var. nipponensis) is widely found in East Asia, suggesting speciation and domestication in these temperate climate regions. Here, we report a draft genome sequence of adzuki bean. The genome assembly covers 75% of the estimated genome and was mapped to 11 pseudo-chromosomes. Gene prediction revealed 26,857 high confidence protein-coding genes evidenced by RNAseq of different tissues. Comparative gene expression analysis with V. radiata showed that the tissue specificity of orthologous genes was highly conserved. Additional re-sequencing of wild adzuki bean, V. angularis var. nipponensis, and V. nepalensis, was performed to analyze the variations between cultivated and wild adzuki bean. The determined divergence time of adzuki bean and the wild species predated archaeology-based domestication time. The present genome assembly will accelerate the genomics-assisted breeding of adzuki bean.

QTL-seq for rapid identification of candidate genes for 100-seed weight and root/total plant dry weight ratio under rainfed conditions in chickpea

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.

Translational Genomics in Agriculture: Some Examples in Grain Legumes

Recent advances in genomics and associated disciplines like bioinformatics have made it possible to develop genomic resources, such as large-scale sequence data for any crop species. While these datasets have been proven very useful for the understanding of genome architecture and dynamics as well as facilitating the discovery of genes, an obligation for, and challenge to the scientific community is to translate genome information to develop products, i.e. superior lines for trait(s) of interest. We call this approach, “translational genomics in agriculture” (TGA). TGA is currently in practice for cereal crops, such as maize (Zea mays) and rice (Oryza sativa), mainly in developed countries and by the private sector; progress has been slow for legume crops. Grown globally on 62.8 million ha with a production of 53.2 million tons and a value of nearly 24.2 billion dollars, the majority of these legumes have low crop productivity (<1 ton/ hectare) and are in the developing countries of sub Saharan Africa, Asia and South America. Interestingly, the last five years have seen enormous progress in genomics for these legume crops. Therefore, it is time to implement TGA in legume crops in order to enhance crop productivity and to ensure food security in developing countries. Prospects, as well as some success stories of TGA, in addition to advances in genomics, trait mapping and gene expression analysis are discussed for five leading legume crops, chickpea (Cicer arietinum), common bean (Phaseolus vulgaris), groundnut (Arachis hypogaea), pigeonpea (Cajanus cajan) and soybean (Glycine max). Some efforts have also been outlined to initiate/ accelerate TGA in three additional legume crops namely faba bean (Vicia faba), lentil (Lens culinaris) and pea (Pisum sativum).

Single Nucleotide Polymorphism Genotyping for Breeding and Genetics Applications in Chickpea and Pigeonpea using the BeadXpress Platform

Single nucleotide polymorphisms (SNPs) are ideal molecular markers due to their higher abundance. Although several types of genotyping platforms for assaying large number of SNPs are available, in cases such as marker‐assisted selection, where few markers are required for genotyping a set of potential lines, high‐throughput SNP genotyping platforms (e.g., iScan or Infinium) may not be cost effective. In this scenario, GoldenGate assays based on VeraCode technology using Illumina BeadXpress seems to be the most cost‐effective platform. The objective of this study was to develop cost‐effective SNP genotyping platforms in chickpea (Cicer arietinum L.) and pigeonpea (Cajanus cajan L.). Two sets of SNPs, one each for chickpea (96 SNPs) and pigeonpea (48 SNPs), were developed and tested by genotyping 288 diverse genotypes from respective reference sets. The SNPs selected for the oligo pool assays had high transferability to crop wild relative species. The mean polymorphism information content value of assayed SNP markers was 0.31 and 0.32 in chickpea and pigeonpea, respectively. No unique pattern was observed in the chickpea reference set whereas two major groups were observed in the case of the pigeonpea reference set. The Illumina BeadXpress platform assays developed for chickpea and pigeonpea are highly informative and cost effective for undertaking genetic studies in these legume species.

Development and Evaluation of a High Density Genotyping 'Axiom_Arachis' Array with 58 K SNPs for Accelerating Genetics and Breeding in Groundnut.

Single nucleotide polymorphisms (SNPs) are the most abundant DNA sequence variation in the genomes which can be used to associate genotypic variation to the phenotype. Therefore, availability of a high-density SNP array with uniform genome coverage can advance genetic studies and breeding applications. Here we report the development of a high-density SNP array ‘Axiom_Arachis’ with 58 K SNPs and its utility in groundnut genetic diversity study. In this context, from a total of 163,782 SNPs derived from DNA resequencing and RNA-sequencing of 41 groundnut accessions and wild diploid ancestors, a total of 58,233 unique and informative SNPs were selected for developing the array. In addition to cultivated groundnuts (Arachis hypogaea), fair representation was kept for other diploids (A. duranensis, A. stenosperma, A. cardenasii, A. magna and A. batizocoi). Genotyping of the groundnut ‘Reference Set’ containing 300 genotypes identified 44,424 polymorphic SNPs and genetic diversity analysis provided in-depth insights into the genetic architecture of this material. The availability of the high-density SNP array ‘Axiom_Arachis’ with 58 K SNPs will accelerate the process of high resolution trait genetics and molecular breeding in cultivated groundnut.

A chromosomal genomics approach to assess and validate the desi and kabuli draft chickpea genome assemblies

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.

Next‐generation sequencing for identification of candidate genes for Fusarium wilt and sterility mosaic disease in pigeonpea (Cajanus cajan)

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.

QTL-seq approach identified genomic regions and diagnostic markers for rust and late leaf spot resistance in groundnut (Arachis hypogaea L.)

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.

Pearl millet genome sequence provides a resource to improve agronomic traits in arid environments

Pearl millet [Cenchrus americanus (L.) Morrone] is a staple food for more than 90 million farmers in arid and semi-arid regions of sub-Saharan Africa, India and South Asia. We report the ∼1.79 Gb draft whole genome sequence of reference genotype Tift 23D2B1-P1-P5, which contains an estimated 38,579 genes. We highlight the substantial enrichment for wax biosynthesis genes, which may contribute to heat and drought tolerance in this crop. We resequenced and analyzed 994 pearl millet lines, enabling insights into population structure, genetic diversity and domestication. We use these resequencing data to establish marker trait associations for genomic selection, to define heterotic pools, and to predict hybrid performance. We believe that these resources should empower researchers and breeders to improve this important staple crop.