Integrated multi-omics analysis reveals drought stress response mechanism in chickpea (Cicer arietinum L.)

Drought is one of the major constraints limiting chickpea productivity. To unravel complex mechanisms regulating drought response in chickpea, we generated transcriptomics, proteomics, and metabolomics datasets from root tissues of four contrasting drought-responsive chickpea genotypes: ICC 4958, JG 11, and JG 11+ (drought-tolerant), and ICC 1882 (drought-sensitive) under control and drought stress conditions. Integration of transcriptomics and proteomics data identified enriched hub proteins encoding isoflavone 4'-O-methyltransferase, UDP-d-glucose/UDP-d-galactose 4-epimerase, and delta-1-pyrroline-5-carboxylate synthetase. These proteins highlighted the involvement of pathways such as antibiotic biosynthesis, galactose metabolism, and isoflavonoid biosynthesis in activating drought stress response mechanisms. Subsequently, the integration of metabolomics data identified six metabolites (fructose, galactose, glucose, myoinositol, galactinol, and raffinose) that showed a significant correlation with galactose metabolism. Integration of root-omics data also revealed some key candidate genes underlying the drought-responsive "QTL-hotspot" region. These results provided key insights into complex molecular mechanisms underlying drought stress response in chickpea.

Genetic structure and ecological niche space of lentil's closest wild relative, Lens orientalis (Boiss.) Schmalh

Crops arose from wild ancestors and to understand their domestication it is essential to compare the cultivated species with their crop wild relatives. These represent an important source of further crop improvement, in particular in relation to climate change. Although there are about 58,000 Lens accessions held in genebanks, only 1% are wild. We examined the geographic distribution and genetic diversity of the lentil's immediate progenitor L. orientalis. We used Genotyping by Sequencing (GBS) to identify and characterize differentiation among accessions held at germplasm collections. We then determined whether genetically distinct clusters of accessions had been collected from climatically distinct locations. Of the 195 genotyped accessions, 124 were genuine L. orientalis with four identified genetic groups. Although an environmental distance matrix was significantly correlated with geographic distance in a Mantel test, the four identified genetic clusters were not found to occupy significantly different environmental space. Maxent modelling gave a distinct predicted distribution pattern centred in the Fertile Crescent, with intermediate probabilities of occurrence in parts of Turkey, Greece, Cyprus, Morocco, and the south of the Iberian Peninsula with NW Africa. Future projections did not show any dramatic alterations in the distribution according to the climate change scenarios tested. We have found considerable diversity in L. orientalis, some of which track climatic variability. The results of the study showed the genetic diversity of wild lentil and indicate the importance of ongoing collections and in situ conservation for our future capacity to harness the genetic variation of the lentil progenitor.

A genomic variation map provides insights into peanut diversity in China and associations with 28 agronomic traits

Peanut (Arachis hypogaea L.) is an important allotetraploid oil and food legume crop. China is one of the world's largest peanut producers and consumers. However, genomic variations underlying the migration and divergence of peanuts in China remain unclear. Here we reported a genome-wide variation map based on the resequencing of 390 peanut accessions, suggesting that peanuts might have been introduced into southern and northern China separately, forming two cultivation centers. Selective sweep analysis highlights asymmetric selection between the two subgenomes during peanut improvement. A classical pedigree from South China offers a context for the examination of the impact of artificial selection on peanut genome. Genome-wide association studies identified 22,309 significant associations with 28 agronomic traits, including candidate genes for plant architecture and oil biosynthesis. Our findings shed light on peanut migration and diversity in China and provide valuable genomic resources for peanut improvement.

Whole genome resequencing and phenotyping of MAGIC population for high resolution mapping of drought tolerance in chickpea

Terminal drought is one of the major constraints to crop production in chickpea (Cicer arietinum L.). In order to map drought tolerance related traits at high resolution, we sequenced multi-parent advanced generation intercross (MAGIC) population using whole genome resequencing approach and phenotyped it under drought stress environments for two consecutive years (2013-14 and 2014-15). A total of 52.02 billion clean reads containing 4.67 TB clean data were generated on the 1136 MAGIC lines and eight parental lines. Alignment of clean data on to the reference genome enabled identification of a total, 932,172 of SNPs, 35,973 insertions, and 35,726 deletions among the parental lines. A high-density genetic map was constructed using 57,180 SNPs spanning a map distance of 1606.69 cM. Using compressed mixed linear model, genome-wide association study (GWAS) enabled us to identify 737 markers significantly associated with days to 50% flowering, days to maturity, plant height, 100 seed weight, biomass, and harvest index. In addition to the GWAS approach, an identity-by-descent (IBD)-based mixed model approach was used to map quantitative trait loci (QTLs). The IBD-based mixed model approach detected major QTLs that were comparable to those from the GWAS analysis as well as some exclusive QTLs with smaller effects. The candidate genes like FRIGIDA and CaTIFY4b can be used for enhancing drought tolerance in chickpea. The genomic resources, genetic map, marker-trait associations, and QTLs identified in the study are valuable resources for the chickpea community for developing climate resilient chickpeas.

Celebrating Professor Rajeev K. Varshney's transformative research odyssey from genomics to the field on his induction as Fellow of the Royal Society

Professor Rajeev K. Varshney's transformative impact on crop genomics, genetics, and agriculture is the result of his passion, dedication, and unyielding commitment to harnessing the potential of genomics to address the most pressing challenges faced by the global agricultural community. Starting from a small town in India and reaching the global stage, Professor Varshney's academic and professional trajectory has inspired many scientists active in research today. His ground-breaking work, especially his effort to list orphan tropical crops to genomic resource-rich entities, has been transformative. Beyond his scientific achievements, Professor Varshney is recognized by his colleagues as an exemplary mentor, fostering the growth of future researchers, building institutional capacity, and strengthening scientific capability. His focus on translational genomics and strengthening seed system in developing countries for the improvement of agriculture has made a tangible impact on farmers' lives. His skills have been best utilized in roles at leading research centres where he has applied his expertise to deliver a new vision for crop improvement. These efforts have now been recognized by the Royal Society with the award of the Fellowship (FRS). As we mark this significant milestone in his career, we not only celebrate Professor Varshney's accomplishments but also his wider contributions that continue to transform the agricultural landscape.

A Single-Nucleus Resolution Atlas of Transcriptome and Chromatin Accessibility for Peanut (Arachis Hypogaea L.) Leaves

The peanut is an important worldwide cash-crop for edible oil and protein. However, the kinetic mechanisms that determine gene expression and chromatin accessibility during leaf development in peanut represented allotetraploid leguminous crops are poorly understood at single-cell resolution. Here, a single-nucleus atlas of peanut leaves is developed by simultaneously profiling the transcriptome and chromatin accessibility in the same individual-cell using fluorescence-activated sorted single-nuclei. In total, 5930 cells with 50 890 expressed genes are classified into 18 cell-clusters, and 5315 chromatin fragments are enriched with 26 083 target genes in the chromatin accessible landscape. The developmental trajectory analysis reveals the involvement of the ethylene-AP2 module in leaf cell differentiation, and cell-cycle analysis demonstrated that genome replication featured in distinct cell-types with circadian rhythms transcription factors (TFs). Furthermore, dual-omics illustrates that the fatty acid pathway modulates epidermal-guard cells differentiation and providescritical TFs interaction networks for understanding mesophyll development, and the cytokinin module (LHY/LOG) that regulates vascular growth. Additionally, an AT-hook protein AhAHL11 is identified that promotes leaf area expansion by modulating the auxin content increase. In summary, the simultaneous profiling of transcription and chromatin accessibility landscapes using snRNA/ATAC-seq provides novel biological insights into the dynamic processes of peanut leaf cell development at the cellular level.

Whole-genome sequencing based discovery of candidate genes and diagnostic markers for seed weight in groundnut

Seed weight in groundnut (Arachis hypogaea L.) has direct impact on yield as well as market price because of preference for bold seeds by consumers and industry, thereby making seed-size improvement as one of the most important objectives of groundnut breeding programs globally. Marker-based early generation selection can accelerate the process of breeding for developing large-seeded varieties. In this context, we deployed the quantitative trait locus-sequencing (QTL-seq) approach on a biparental mapping population (Chico × ICGV 02251) to identify candidate genes and develop markers for seed weight in groundnut. A total of 289.4-389.4 million reads sequencing data were generated from three libraries (ICGV 02251 and two extreme bulks) achieving 93.9-95.1% genome coverage and 8.34-9.29× average read depth. The analysis of sequencing data using QTL-seq pipeline identified five genomic regions (three on chromosome B06 and one each on chromosomes B08 and B09) for seed weight. Detailed analysis of above associated genomic regions detected 182 single-nucleotide polymorphisms (SNPs) in genic and intergenic regions, and 11 of these SNPs were nonsynonymous in the genomic regions of 10 candidate genes including Ulp proteases and BIG SEED locus genes. Kompetitive allele specific polymerase chain reaction (KASP) markers for 14 SNPs were developed, and four of these markers (snpAH0031, snpAH0033, snpAH0037, and snpAH0038) were successfully validated for deployment in breeding for large-seeded groundnut varieties.

Delving into the lifestyle of Sundarban Wetland resident, biofilm producing, halotolerant Salinicoccus roseus: a comparative genomics-based intervention

BACKGROUND

Microbial community played an essential role in ecosystem processes, be it mangrove wetland or other intertidal ecologies. Several enzymatic activities like hydrolases are effective ecological indicators of soil microbial function. So far, little is known on halophilic bacterial contribution and function on a genomic viewpoint of Indian Sundarban Wetland. Considering the above mentioned issues, the aims of this study was to understand the life style, metabolic functionalities and genomic features of the isolated bacterium, Salinicoccus roseus strain RF1H. A comparative genome-based study of S. roseus has not been reported yet. Henceforth, we have considered the inclusion of the intra-species genome comparison of S. roseus to gain insight into the high degree of variation in the genome of strain RF1H among others.

RESULTS

Salinicoccus roseus strain RF1H is a pink-red pigmented, Gram-positive and non-motile cocci. The bacterium exhibited high salt tolerance (up to 15% NaCl), antibiotic resistance, biofilm formation and secretion of extracellular hydrolytic enzymes. The circular genome was approximately 2.62978 Mb in size, encoding 574 predicted genes with GC content 49.5%. Presence of genomic elements (prophages, transposable elements, CRISPR-Cas system) represented bacterial virulence and multidrug-resistance. Furthermore, genes associated with salt tolerance, temperature adaptation and DNA repair system were distributed in 17 genomic islands. Genes related to hydrocarbon degradation manifested metabolic capability of the bacterium for potential biotechnological applications. A comparative pangenome analysis revealed two-component response regulator, modified C4-dicarboxylate transport system and osmotic stress regulated ATP-binding proteins. Presence of genes encoding arginine decarboxylase (ADC) enzyme being involved in biofilm formation was reported from the genome. In silico study revealed the protein is thermostable and made up with ~ 415 amino acids, and hydrophilic in nature. Three motifs appeared to be evolutionary conserved in all Salinicoccus sequences.

CONCLUSION

The first report of whole genome analysis of Salinicoccus roseus strain RF1H provided information of metabolic functionalities, biofilm formation, resistance mechanism and adaptation strategies to thrive in climate-change induced vulnerable spot like Sundarban. Comparative genome analysis highlighted the unique genome content that contributed the strain's adaptability. The biomolecules produced during metabolism are important sources of compounds with potential beneficial applications in pharmaceuticals.

Insight into the genome of an arsenic loving and plant growth-promoting strain of Micrococcus luteus isolated from arsenic contaminated groundwater

Contamination of arsenic in drinking water and foods is a threat for human beings. To achieve the goal for the reduction of arsenic availability, besides conventional technologies, arsenic bioremediation by using some potent bacteria is one of the hot topics for researchers. In this context, bacterium, AKS4c was isolated from arsenic contaminated water of Purbasthali, West Bengal, India, and through draft genome sequence; it was identified as a strain of Micrococcus luteus that comprised of 2.4 Mb genome with 73.1% GC content and 2256 protein coding genes. As the accessory genome, about 22 genomic islands (GIs) associated with many metal-resistant genes were identified. This strain was capable to tolerate more than 46,800 mg/L arsenate and 390 mg/L arsenite salts as well as found to be tolerable to multi-metals such as Fe, Pb, Mo, Mn, and Zn up to a certain limit of concentrations. Strain AKS4c was able to oxidize arsenite to less toxic arsenate, and its arsenic adsorption property was qualitatively confirmed through X-ray fluorescence (XRF) and Fourier transform infrared spectroscopy (FTIR) analysis. Quantitative estimation of plant growth-promoting attributes like Indole acetic acid (IAA), Gibberellic acid (GA), and proline production and enhancement of rice seedling growth in laboratory condition leads to its future applicability in arsenic bioremediation as a plant growth-promoting rhizobacteria (PGPR).

Haemoglobin diagnostic cut-offs for anaemia in Indian women of reproductive age

BACKGROUND

The persistent high prevalence of anaemia among Indian women of reproductive age (WRA) despite aggressive long-term iron supplementation could be related to over-diagnosis from an inappropriately high haemoglobin (Hb) diagnostic cut-off. To develop an appropriate cut-off for Indian WRA, we hypothesized that during iron-folic acid (IFA) supplementation to a mixed (anaemic/non-anaemic) WRA population, the positive slope of the Hb-plasma ferritin (PF) response in anaemic women would inflect into a plateau (zero-response) as a non-anaemic status is reached. The 2.5th percentile of the Hb distribution at this inflection point will be the diagnostic Hb cut-off for iron-responsive anaemia.

METHOD

A hierarchical mixed effects model, with a polynomial mean and variance model to account for intraclass correlation due to repeated measures, was used to estimate the response curve of Hb to PF, or body iron stores, in anaemic and non-anaemic WRA (without inflammation), who were receiving a 90-day IFA supplementation.

RESULTS

The Hb response curve at low PF values showed a steep increase, which inflected into a plateau at a PF of 10.1 µg/L and attained a steady state at a PF of 20.6 µg/L. The Hb distribution at the inflection was a normal probability distribution, with a mean of 12.3 g/dL. The 2.5th percentile value of this distribution, or the putative diagnostic Hb cut-off for anaemia, was 10.8 g/dL (~11 g/dL).

CONCLUSION

The derived Hb cut-off is lower than the current adult values of 12 g/dL and could partly explain the persistently high prevalence of anaemia.

Improved pearl millet genomes representing the global heterotic pool offer a framework for molecular breeding applications

High-quality reference genome assemblies, representative of global heterotic patterns, offer an ideal platform to accurately characterize and utilize genetic variation in the primary gene pool of hybrid crops. Here we report three platinum grade de-novo, near gap-free, chromosome-level reference genome assemblies from the active breeding germplasm in pearl millet with a high degree of contiguity, completeness, and accuracy. An improved Tift genome (Tift23D2B1-P1-P5) assembly has a contig N50 ~ 7,000-fold (126 Mb) compared to the previous version and better alignment in centromeric regions. Comparative genome analyses of these three lines clearly demonstrate a high level of collinearity and multiple structural variations, including inversions greater than 1 Mb. Differential genes in improved Tift genome are enriched for serine O-acetyltransferase and glycerol-3-phosphate metabolic process which play an important role in improving the nutritional quality of seed protein and disease resistance in plants, respectively. Multiple marker-trait associations are identified for a range of agronomic traits, including grain yield through genome-wide association study. Improved genome assemblies and marker resources developed in this study provide a comprehensive framework/platform for future applications such as marker-assisted selection of mono/oligogenic traits as well as whole-genome prediction and haplotype-based breeding of complex traits.

Transcriptome profiling reveals the expression and regulation of genes associated with Fusarium wilt resistance in chickpea (Cicer arietinum L.)

Fusarium wilt (FW) is one of the most significant biotic stresses limiting chickpea production worldwide. To dissect the molecular mechanism of FW resistance in chickpea, comparative transcriptome analyses of contrasting resistance sources of chickpea genotypes under control and Fusarium oxysporum f. sp. ciceris (Foc) inoculated conditions were performed. The high-throughput transcriptome sequencing generated about 1137 million sequencing reads from 24 samples representing two resistant genotypes, two susceptible genotypes, and two near-isogenic lines under control and stress conditions at two-time points (7th- and 12th-day post-inoculation). The analysis identified 5182 differentially expressed genes (DEGs) between different combinations of chickpea genotypes. Functional annotation of these genes indicated their involvement in various biological processes such as defense response, cell wall biogenesis, secondary metabolism, and disease resistance. A significant number (382) of transcription factor encoding genes exhibited differential expression patterns under stress. Further, a considerable number of the identified DEGs (287) co-localized with previously reported quantitative trait locus for FW resistance. Several resistance/susceptibility-related genes, such as SERINE/THREONINE PROTEIN KINASE, DIRIGENT, and MLO exhibiting contrasting expression patterns in resistant and susceptible genotypes upon Foc inoculation, were identified. The results presented in the study provide valuable insights into the transcriptional dynamics associated with FW stress response in chickpea and provide candidate genes for the development of disease-resistant chickpea cultivars.

Copper removal capability and genomic insight into the lifestyle of copper mine inhabiting Micrococcus yunnanensis GKSM13

Heavy metal pollution in mining areas is a serious environmental concern. The exploration of mine-inhabiting microbes, especially bacteria may use as an effective alternative for the remediation of mining hazards. A highly copper-tolerant strain GKSM13 was isolated from the soil of the Singhbhum copper mining area and characterized for significant copper (Cu) removal potential and tolerance to other heavy metals. The punctate, yellow-colored, coccoid strain GKSM13 was able to tolerate 500 mg L^-1 Cu²⁺. Whole-genome sequencing identified strain GKSM13 as Micrococcus yunnanensis, which has a 2.44 Mb genome with 2176 protein-coding genes. The presence of putative Cu homeostasis genes and other heavy metal transporters/response regulators or transcription factors may responsible for multi-metal resistance. The maximum Cu²⁺ removal of 89.2% was achieved at a pH of 7.5, a temperature of 35.5 °C, and an initial Cu²⁺ ion concentration of 31.5 mg L^-1. Alteration of the cell surface, deposition of Cu²⁺ in the bacterial cell, and the involvement of hydroxyl, carboxyl amide, and amine groups in Cu²⁺ removal were observed using microscopic and spectroscopic analysis. This study is the first to reveal a molecular-based approach for the multi-metal tolerance and copper homeostasis mechanism of M. yunnanensis GKSM13.

Developing future heat-resilient vegetable crops

Climate change seriously impacts global agriculture, with rising temperatures directly affecting the yield. Vegetables are an essential part of daily human consumption and thus have importance among all agricultural crops. The human population is increasing daily, so there is a need for alternative ways which can be helpful in maximizing the harvestable yield of vegetables. The increase in temperature directly affects the plants' biochemical and molecular processes; having a significant impact on quality and yield. Breeding for climate-resilient crops with good yields takes a long time and lots of breeding efforts. However, with the advent of new omics technologies, such as genomics, transcriptomics, proteomics, and metabolomics, the efficiency and efficacy of unearthing information on pathways associated with high-temperature stress resilience has improved in many of the vegetable crops. Besides omics, the use of genomics-assisted breeding and new breeding approaches such as gene editing and speed breeding allow creation of modern vegetable cultivars that are more resilient to high temperatures. Collectively, these approaches will shorten the time to create and release novel vegetable varieties to meet growing demands for productivity and quality. This review discusses the effects of heat stress on vegetables and highlights recent research with a focus on how omics and genome editing can produce temperature-resilient vegetables more efficiently and faster.

Comprehensive Transcriptome Profiling Uncovers Molecular Mechanisms and Potential Candidate Genes Associated with Heat Stress Response in Chickpea

Chickpea (Cicer arietinum L.) production is highly susceptible to heat stress (day/night temperatures above 32/20 °C). Identifying the molecular mechanisms and potential candidate genes underlying heat stress response is important for increasing chickpea productivity. Here, we used an RNA-seq approach to investigate the transcriptome dynamics of 48 samples which include the leaf and root tissues of six contrasting heat stress responsive chickpea genotypes at the vegetative and reproductive stages of plant development. A total of 14,544 unique, differentially expressed genes (DEGs) were identified across different combinations studied. These DEGs were mainly involved in metabolic processes, cell wall remodeling, calcium signaling, and photosynthesis. Pathway analysis revealed the enrichment of metabolic pathways, biosynthesis of secondary metabolites, and plant hormone signal transduction, under heat stress conditions. Furthermore, heat-responsive genes encoding bHLH, ERF, WRKY, and MYB transcription factors were differentially regulated in response to heat stress, and candidate genes underlying the quantitative trait loci (QTLs) for heat tolerance component traits, which showed differential gene expression across tolerant and sensitive genotypes, were identified. Our study provides an important resource for dissecting the role of candidate genes associated with heat stress response and also paves the way for developing climate-resilient chickpea varieties for the future.

Characterization of 'QTL-hotspot' introgression lines reveals physiological mechanisms and candidate genes associated with drought adaptation in chickpea

'QTL-hotspot' is a genomic region on linkage group 04 (CaLG04) in chickpea (Cicer arietinum) that harbours major-effect quantitative trait loci (QTLs) for multiple drought-adaptive traits, and it therefore represents a promising target for improving drought adaptation. To investigate the mechanisms underpinning the positive effects of 'QTL-hotspot' on seed yield under drought, we introgressed this region from the ICC 4958 genotype into five elite chickpea cultivars. The resulting introgression lines (ILs) and their parents were evaluated in multi-location field trials and semi-controlled conditions. The results showed that the 'QTL-hotspot' region improved seed yield under rainfed conditions by increasing seed weight, reducing the time to flowering, regulating traits related to canopy growth and early vigour, and enhancing transpiration efficiency. Whole-genome sequencing data analysis of the ILs and parents revealed four genes underlying the 'QTL-hotspot' region associated with drought adaptation. We validated diagnostic KASP markers closely linked to these genes using the ILs and their parents for future deployment in chickpea breeding programs. The CaTIFY4b-H2 haplotype of a potential candidate gene CaTIFY4b was identified as the superior haplotype for 100-seed weight. The candidate genes and superior haplotypes identified in this study have the potential to serve as direct targets for genetic manipulation and selection for chickpea improvement.

Genome-Wide Characterization of Ascorbate Peroxidase Gene Family in Peanut (Arachis hypogea L.) Revealed Their Crucial Role in Growth and Multiple Stress Tolerance

Ascorbate peroxidase (APX), an important antioxidant enzyme, plays a significant role in ROS scavenging by catalyzing the decrease of hydrogen peroxide under various environmental stresses. Nevertheless, information about the APX gene family and their evolutionary and functional attributes in peanut (Arachis hypogea L.) was not reported. Therefore, a comprehensive genome-wide study was performed to discover the APX genes in cultivated peanut genome. This study identified 166 AhAPX genes in the peanut genome, classified into 11 main groups. The gene duplication analysis showed that AhAPX genes had experienced segmental duplications and purifying selection pressure. Gene structure and motif investigation indicated that most of the AhAPX genes exhibited a comparatively well-preserved exon-intron pattern and motif configuration contained by the identical group. We discovered five phytohormones-, six abiotic stress-, and five growth and development-related cis-elements in the promoter regions of AhAPX. Fourteen putative ah-miRNAs from 12 families were identified, targeting 33 AhAPX genes. Furthermore, we identified 3,257 transcription factors from 38 families (including AP2, ARF, B3, bHLH, bZIP, ERF, MYB, NAC, WRKY, etc.) in 162 AhAPX genes. Gene ontology and KEGG enrichment analysis confirm the role of AhAPX genes in oxidoreductase activity, catalytic activity, cell junction, cellular response to stimulus and detoxification, biosynthesis of metabolites, and phenylpropanoid metabolism. Based on transcriptome datasets, some genes such as AhAPX4/7/17/77/82/86/130/133 and AhAPX160 showed significantly higher expression in diverse tissues/organs, i.e., flower, leaf, stem, roots, peg, testa, and cotyledon. Likewise, only a few genes, including AhAPX4/17/19/55/59/82/101/102/137 and AhAPX140, were significantly upregulated under abiotic (drought and cold), and phytohormones (ethylene, abscisic acid, paclobutrazol, brassinolide, and salicylic acid) treatments. qRT-PCR-based expression profiling presented the parallel expression trends as generated from transcriptome datasets. Our discoveries gave new visions into the evolution of APX genes and provided a base for further functional examinations of the AhAPX genes in peanut breeding programs.

Genome-wide association analysis to delineate high-quality SNPs for seed micronutrient density in chickpea (Cicer arietinum L.)

Chickpea is the most important nutrient-rich grain legume crop in the world. A diverse core set of 147 chickpea genotypes was genotyped with a Axiom(®)50K CicerSNP array and trait phenotyped in two different environments for four seed micronutrients (Zn, Cu, Fe and Mn). The trait data and high-throughput 50K SNP genotypic data were used for the genome-wide association study (GWAS). The study led to the discovery of genes/QTLs for seed Zn, Cu, Fe and Mn, concentrations in chickpea. The analysis of seed micronutrient data revealed significant differences for all four micronutrient concentrations (P ≤ 0.05). The mean concentrations of seed Zn, Cu, Fe and Mn pooled over the 2 years were 45.9 ppm, 63.8 ppm 146.1 ppm, and 27.0 ppm, respectively. The analysis of results led to the identification of 35 SNPs significantly associated with seed Zn, Cu, Fe and Mn concentrations. Among these 35 marker-trait associations (MTAs), 5 were stable (consistently identified in different environments), 6 were major (explaining more than 15% of the phenotypic variation for an individual trait) and 3 were both major and stable MTAs. A set of 6 MTAs, MTAs (3 for Mn, 2 for Fe, and 1 for Cu) reported by us during the present study have been also reported in the same/almost same genomic regions in earlier studies and therefore declared as validated MTAs. The stable, major and validated MTAs identified during the present study will prove useful in future chickpea molecular breeding programs aimed at enhancing the seed nutrient density of chickpea.

Genetic variation in CaTIFY4b contributes to drought adaptation in chickpea

Chickpea production is vulnerable to drought stress. Identifying the genetic components underlying drought adaptation is crucial for enhancing chickpea productivity. Here, we present the fine mapping and characterization of 'QTL-hotspot', a genomic region controlling chickpea growth with positive consequences on crop production under drought. We report that a non-synonymous substitution in the transcription factor CaTIFY4b regulates seed weight and organ size in chickpea. Ectopic expression of CaTIFY4b in Medicago truncatula enhances root growth under water deficit. Our results suggest that allelic variation in 'QTL-hotspot' improves pre-anthesis water use, transpiration efficiency, root architecture and canopy development, enabling high-yield performance under terminal drought conditions. Gene expression analysis indicated that CaTIFY4b may regulate organ size under water deficit by modulating the expression of GRF-INTERACTING FACTOR1 (GIF1), a transcriptional co-activator of Growth-Regulating Factors. Taken together, our study offers new insights into the role of CaTIFY4b and on diverse physiological and molecular mechanisms underpinning chickpea growth and production under specific drought scenarios.

Genome-wide association mapping of nutritional traits for designing superior chickpea varieties

Micronutrient malnutrition is a serious concern in many parts of the world; therefore, enhancing crop nutrient content is an important challenge. Chickpea (Cicer arietinum L.), a major food legume crop worldwide, is a vital source of protein and minerals in the vegetarian diet. This study evaluated a diverse set of 258 chickpea germplasm accessions for 12 key nutritional traits. A significant variation was observed for several nutritional traits, including crude protein (16.56-24.64/100 g), β-Carotene (0.003-0.104 mg/100 g), calcium (60.69-176.55 mg/100 g), and folate (0.413-6.537 mg/kg). These data, combined with the available whole-genome sequencing data for 318,644 SNPs, were used in genome-wide association studies comprising single-locus and multi-locus models. We also explored the effect of varying the minor allele frequency (MAF) levels and heterozygosity. We identified 62 significant marker-trait associations (MTAs) explaining up to 28.63% of the phenotypic variance (PV), of which nine were localized within genes regulating G protein-coupled receptor signaling pathway, proteasome assembly, intracellular signal transduction, and oxidation-reduction process, among others. The significant effect MTAs were located primarily on Ca1, Ca3, Ca4, and Ca6. Importantly, varying the level of heterozygosity was found to significantly affect the detection of associations contributing to traits of interest. We further identified seven promising accessions (ICC10399, ICC1392, ICC1710, ICC2263, ICC1431, ICC4182, and ICC16915) with superior agronomic performance and high nutritional content as potential donors for developing nutrient-rich, high-yielding chickpea varieties. Validation of the significant MTAs with higher PV could identify factors controlling the nutrient acquisition and facilitate the design of biofortified chickpeas for the future.

Development of High Yielding Fusarium Wilt Resistant Cultivar by Pyramiding of “Genes” Through Marker-Assisted Backcrossing in Chickpea (Cicer arietinum L.)

Pusa 391, a mega desi chickpea variety with medium maturity duration is extensively cultivated in the Central Zone of India. Of late, this variety has become susceptible to Fusarium wilt (FW), which has drastic impact on its yield. Presence of variability in the wilt causing pathogen, Fusarium oxysporum f.sp. ciceri (foc) across geographical locations necessitates the role of pyramiding for FW resistance for different races (foc 1,2,3,4 and 5). Subsequently, the introgression lines developed in Pusa 391 genetic background were subjected to foreground selection using three SSR markers (GA16, TA 27 and TA 96) while 48 SSR markers uniformly distributed on all chromosomes, were used for background selection to observe the recovery of recurrent parent genome (RPG). BC1F1 lines with 75–85% RPG recovery were used to generate BC2F1. The plants that showed more than 90% RPG recovery in BC2F1 were used for generating BC3F1. The plants that showed more than 96% RPG recovery were selected and selfed to generate BC3F3. Multi-location evaluation of advanced introgression lines (BC2F3) in six locations for grain yield (kg/ha), days to fifty percent flowering, days to maturity, 100 seed weight and disease incidence was done. In case of disease incidence, the genotype IL1 (BGM 20211) was highly resistant to FW in Junagarh, Indore, New Delhi, Badnapur and moderately resistant at Sehore and Nandyal. GGE biplot analysis revealed that IL1(BGM20211) was the most stable genotype at Junagadh, Sehore and Nandyal. GGE biplot analysis revealed that IL1(BGM 20211) and IL4(BGM 20212) were the top performers in yield and highly stable across six environments and were nominated for Advanced Varietal Trials (AVT) of AICRP (All India Coordinated Research Project on Chickpea) in 2018–19. BGM20211 and BGM 20212 recorded 29 and 28.5% average yield gain over the recurrent parent Pusa 391, in the AVT-1 and AVT-2 over five environments. Thus, BGM20211 was identified for release and notified as Pusa Manav/Pusa Chickpea 20211 for Madhya Pradesh, Gujarat and Maharashtra, Southern Rajasthan, Bundhelkhand region of Uttar Pradesh states by the Central Sub-Committees on Crop Standards, Notification and Release of Varieties of Agricultural Crops, Ministry of Agriculture and Farmers Welfare, Government of India, for commercial cultivation in India (Gazette notification number S.O.500 (E) dt. 29-1-2021).Such pyramided lines give resilience to multiple races of fusarium wilt with added yield advantage.

Fast-forward breeding for a food-secure world

Crop production systems need to expand their outputs sustainably to feed a burgeoning human population. Advances in genome sequencing technologies combined with efficient trait mapping procedures accelerate the availability of beneficial alleles for breeding and research. Enhanced interoperability between different omics and phenotyping platforms, leveraged by evolving machine learning tools, will help provide mechanistic explanations for complex plant traits. Targeted and rapid assembly of beneficial alleles using optimized breeding strategies and precise genome editing techniques could deliver ideal crops for the future. Realizing desired productivity gains in the field is imperative for securing an adequate future food supply for 10 billion people.

Novel Genes and Genetic Loci Associated With Root Morphological Traits, Phosphorus-Acquisition Efficiency and Phosphorus-Use Efficiency in Chickpea

Chickpea-the second most important grain legume worldwide-is cultivated mainly on marginal soils. Phosphorus (P) deficiency often restricts chickpea yields. Understanding the genetics of traits encoding P-acquisition efficiency and P-use efficiency will help develop strategies to reduce P-fertilizer application. A genome-wide association mapping approach was used to determine loci and genes associated with root architecture, root traits associated with P-acquisition efficiency and P-use efficiency, and any associated proxy traits. Using three statistical models-a generalized linear model (GLM), a mixed linear model (MLM), and a fixed and random model circulating probability unification (FarmCPU) -10, 51, and 40 marker-trait associations (MTAs), respectively were identified. A single nucleotide polymorphism (SNP) locus (Ca1_12310101) on Ca1 associated with three traits, i.e., physiological P-use efficiency, shoot dry weight, and shoot P content was identified. Genes related to shoot P concentration (NAD kinase 2, dynamin-related protein 1C), physiological P-use efficiency (fasciclin-like arabinogalactan protein), specific root length (4-coumarate-CoA ligase 1) and manganese concentration in mature leaves (ABC1 family protein) were identified. The MTAs and novel genes identified in this study can be used to improve P-use efficiency in chickpea.

Comparative Transcriptome Analysis Identified Candidate Genes for Late Leaf Spot Resistance and Cause of Defoliation in Groundnut

Late leaf spot (LLS) caused by fungus Nothopassalora personata in groundnut is responsible for up to 50% yield loss. To dissect the complex nature of LLS resistance, comparative transcriptome analysis was performed using resistant (GPBD 4), susceptible (TAG 24) and a resistant introgression line (ICGV 13208) and identified a total of 12,164 and 9954 DEGs (differentially expressed genes) respectively in A- and B-subgenomes of tetraploid groundnut. There were 135 and 136 unique pathways triggered in A- and B-subgenomes, respectively, upon N. personata infection. Highly upregulated putative disease resistance genes, an RPP-13 like (Aradu.P20JR) and a NBS-LRR (Aradu.Z87JB) were identified on chromosome A02 and A03, respectively, for LLS resistance. Mildew resistance Locus (MLOs)-like proteins, heavy metal transport proteins, and ubiquitin protein ligase showed trend of upregulation in susceptible genotypes, while tetratricopeptide repeats (TPR), pentatricopeptide repeat (PPR), chitinases, glutathione S-transferases, purple acid phosphatases showed upregulation in resistant genotypes. However, the highly expressed ethylene responsive factor (ERF) and ethylene responsive nuclear protein (ERF2), and early responsive dehydration gene (ERD) might be related to the possible causes of defoliation in susceptible genotypes. The identified disease resistance genes can be deployed in genomics-assisted breeding for development of LLS resistant cultivars to reduce the yield loss in groundnut.

Introgression of "QTL-hotspot" region enhances drought tolerance and grain yield in three elite chickpea cultivars

With an aim of enhancing drought tolerance using a marker-assisted backcrossing (MABC) approach, we introgressed the "QTL-hotspot" region from ICC 4958 accession that harbors quantitative trait loci (QTLs) for several drought-tolerance related traits into three elite Indian chickpea (Cicer arietinum L.) cultivars: Pusa 372, Pusa 362, and DCP 92-3. Of eight simple sequence repeat (SSR) markers in the QTL-hotspot region, two to three polymorphic markers were used for foreground selection with respective cross-combinations. A total of 47, 53, and 46 SSRs were used for background selection in case of introgression lines (ILs) developed in genetic backgrounds of Pusa 372, Pusa 362, and DCP 92-3, respectively. In total, 61 ILs (20 BC₃ F₃ in Pusa 372; 20 BC₂ F₃ in Pusa 362, and 21 BC₃ F₃ in DCP 92-3), with >90% recurrent parent genome recovery were developed. Six improved lines in different genetic backgrounds (e.g. BGM 10216 in Pusa 372; BG 3097 and BG 4005 in Pusa 362; IPC(L4-14), IPC(L4-16), and IPC(L19-1) in DCP 92-3) showed better performance than their respective recurrent parents. BGM 10216, with 16% yield gain over Pusa 372, has been released as Pusa Chickpea 10216 by the Central Sub-Committees on Crop Standards, Notification and Release of Varieties of Agricultural Crops, Ministry of Agriculture and Farmers Welfare, Government of India, for commercial cultivation in India. In summary, this study reports introgression of the QTL-hotspot for enhancing yield under rainfed conditions, development of several introgression lines, and release of Pusa Chickpea 10216 developed through molecular breeding in India.

Natural polymorphisms in a pair of NSP2 homoeologs can cause loss of nodulation in peanut

Microbial symbiosis in legumes is achieved through nitrogen-fixing root nodules, and these are important for sustainable agriculture. The molecular mechanisms underlying development of root nodules in polyploid legume crops are largely understudied. Through map-based cloning and QTL-seq approaches, we identified a pair of homoeologous GRAS transcription factor genes, Nodulation Signaling Pathway 2 (AhNSP2-B07 or Nb) and AhNSP2-A08 (Na), controlling nodulation in cultivated peanut (Arachis hypogaea L.), an allotetraploid legume crop, which exhibited non-Mendelian and Mendelian inheritance, respectively. The segregation of nodulation in the progeny of Nananbnb genotypes followed a 3:1 Mendelian ratio, in contrast to the 5:3~1:1 non-Mendelian ratio for nanaNbnb genotypes. Additionally, a much higher frequency of the nb allele (13%) than the na allele (4%) exists in the peanut germplasm collection, suggesting that Nb is less essential than Na in nodule organogenesis. Our findings reveal the genetic basis of naturally occurred non-nodulating peanut plants, which can be potentially used for nitrogen fixation improvement in peanut. Furthermore, the results have implications for and provide insights into the evolution of homoeologous genes in allopolyploid species.

Genomic resources in plant breeding for sustainable agriculture

Climate change during the last 40 years has had a serious impact on agriculture and threatens global food and nutritional security. From over half a million plant species, cereals and legumes are the most important for food and nutritional security. Although systematic plant breeding has a relatively short history, conventional breeding coupled with advances in technology and crop management strategies has increased crop yields by 56 % globally between 1965-85, referred to as the Green Revolution. Nevertheless, increased demand for food, feed, fiber, and fuel necessitates the need to break existing yield barriers in many crop plants. In the first decade of the 21st century we witnessed rapid discovery, transformative technological development and declining costs of genomics technologies. In the second decade, the field turned towards making sense of the vast amount of genomic information and subsequently moved towards accurately predicting gene-to-phenotype associations and tailoring plants for climate resilience and global food security. In this review we focus on genomic resources, genome and germplasm sequencing, sequencing-based trait mapping, and genomics-assisted breeding approaches aimed at developing biotic stress resistant, abiotic stress tolerant and high nutrition varieties in six major cereals (rice, maize, wheat, barley, sorghum and pearl millet), and six major legumes (soybean, groundnut, cowpea, common bean, chickpea and pigeonpea). We further provide a perspective and way forward to use genomic breeding approaches including marker-assisted selection, marker-assisted backcrossing, haplotype based breeding and genomic prediction approaches coupled with machine learning and artificial intelligence, to speed breeding approaches. The overall goal is to accelerate genetic gains and deliver climate resilient and high nutrition crop varieties for sustainable agriculture.

Allelic Diversity, Structural Analysis, and Genome-Wide Association Study (GWAS) for Yield and Related Traits Using Unexplored Common Bean (Phaseolus vulgaris L.) Germplasm From Western Himalayas

The north-western Indian Himalayas possesses vast diversity in common bean germplasm due to several years of natural adaptation and farmer's selection. Systematic efforts have been made for the first time for the characterization and use of this huge diversity for the identification of genes/quantitative trait loci (QTLs) for yield and yield-contributing traits in common bean in India. A core set of 96 diverse common bean genotypes was characterized using 91 genome-wide genomic and genic simple sequence repeat (SSR) markers. The study of genetic diversity led to the identification of 691 alleles ranging from 2 to 21 with an average of 7.59 alleles/locus. The gene diversity (expected heterozygosity, He) varied from 0.31 to 0.93 with an average of 0.73. As expected, the genic SSR markers detected less allelic diversity than the random genomic SSR markers. The traditional clustering and Bayesian clustering (structural analysis) analyses led to a clear cut separation of a core set of 96 genotypes into two distinct groups based on their gene pools (Mesoamerican and Andean genotypes). Genome-wide association mapping for pods/plant, seeds/pod, seed weight, and yield/plant led to the identification of 39 significant marker-trait associations (MTAs) including 15 major, 15 stable, and 13 both major and stable MTAs. Out of 39 MTAs detected, 29 were new MTAs reported for the first time, whereas the remaining 10 MTAs were already identified in earlier studies and therefore declared as validation of earlier results. A set of seven markers was such, which were found to be associated with multiple (two to four) different traits. The important MTAs will be used for common bean molecular breeding programs worldwide for enhancing common bean yield.

MutMap Approach Enables Rapid Identification of Candidate Genes and Development of Markers Associated With Early Flowering and Enhanced Seed Size in Chickpea (Cicer arietinum L.)

Globally terminal drought is one of the major constraints to chickpea (Cicer arietinum L.) production. Early flowering genotypes escape terminal drought, and the increase in seed size compensates for yield losses arising from terminal drought. A MutMap population for early flowering and large seed size was developed by crossing the mutant line ICC4958-M3-2828 with wild-type ICC 4958. Based on the phenotyping of MutMap population, extreme bulks for days to flowering and 100-seed weight were sequenced using Hi-Seq2500 at 10X coverage. On aligning 47.41 million filtered reads to the CDC Frontier reference genome, 31.41 million reads were mapped and 332,395 single nucleotide polymorphisms (SNPs) were called. A reference genome assembly for ICC 4958 was developed replacing these SNPs in particular positions of the CDC Frontier genome. SNPs specific for each mutant bulk ranged from 3,993 to 5,771. We report a single unique genomic region on Ca6 (between 9.76 and 12.96 Mb) harboring 31, 22, 17, and 32 SNPs with a peak of SNP index = 1 for low bulk for flowering time, high bulk for flowering time, high bulk for 100-seed weight, and low bulk for 100-seed weight, respectively. Among these, 22 SNPs are present in 20 candidate genes and had a moderate allelic impact on the genes. Two markers, Ca6EF10509893 for early flowering and Ca6HSDW10099486 for 100-seed weight, were developed and validated using the candidate SNPs. Thus, the associated genes, candidate SNPs, and markers developed in this study are useful for breeding chickpea varieties that mitigate yield losses under drought stress.

Major QTLs and Potential Candidate Genes for Heat Stress Tolerance Identified in Chickpea (Cicer arietinum L.)

In the context of climate change, heat stress during the reproductive stages of chickpea (Cicer arietinum L.) leads to significant yield losses. In order to identify the genomic regions responsible for heat stress tolerance, a recombinant inbred line population derived from DCP 92-3 (heat sensitive) and ICCV 92944 (heat tolerant) was genotyped using the genotyping-by-sequencing approach and evaluated for two consecutive years (2017 and 2018) under normal and late sown or heat stress environments. A high-density genetic map comprising 788 single-nucleotide polymorphism markers spanning 1,125 cM was constructed. Using composite interval mapping, a total of 77 QTLs (37 major and 40 minor) were identified for 12 of 13 traits. A genomic region on CaLG07 harbors quantitative trait loci (QTLs) explaining >30% phenotypic variation for days to pod initiation, 100 seed weight, and for nitrogen balance index explaining >10% PVE. In addition, we also reported for the first time major QTLs for proxy traits (physiological traits such as chlorophyll content, nitrogen balance index, normalized difference vegetative index, and cell membrane stability). Furthermore, 32 candidate genes in the QTL regions that encode the heat shock protein genes, heat shock transcription factors, are involved in flowering time regulation as well as pollen-specific genes. The major QTLs reported in this study, after validation, may be useful in molecular breeding for developing heat-tolerant superior lines or varieties.

Superior haplotypes for haplotype-based breeding for drought tolerance in pigeonpea (Cajanus cajan L.)

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 drought-responsive 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 next-generation tailor-made better drought-responsive pigeonpea cultivars.

Development of a dense genetic map and QTL analysis for pod borer Helicoverpa armigera (Hübner) resistance component traits in chickpea (Cicer arietinum L.)

Genetic enhancement for resistance against the pod borer, Helicoverpa armigera is crucial for enhancing production and productivity of chickpea. Here we provide some novel insights into the genetic architecture of natural variation in H. armigera resistance in chickpea, an important legume, which plays a major role in food and nutritional security. An interspecific recombinant inbred line (RIL) population developed from a cross between H. armigera susceptible accession ICC 4958 (Cicer arietinum) and resistant accession PI 489777 (Cicer reticulatum) was evaluated for H. armigera resistance component traits using detached leaf assay and under field conditions. A high-throughput AxiomCicerSNP array was utilized to construct a dense linkage map comprising of 3,873 loci and spanning a distance of 949.27 cM. Comprehensive analyses of extensive genotyping and phenotyping data identified nine main-effect QTLs and 955 epistatic QTLs explaining up to 42.49% and 38.05% phenotypic variance, respectively, for H. armigera resistance component traits. The main-effect QTLs identified in this RIL population were linked with previously described genes, known to modulate resistance against lepidopteran insects in crop plants. One QTL cluster harbouring main-effect QTLs for three H. armigera resistance component traits and explaining up to 42.49% of the phenotypic variance, was identified on CaLG03. This genomic region, after validation, may be useful to improve H. armigera resistance component traits in elite chickpea cultivars.

Transcriptional responses of toxigenic and atoxigenic isolates of Aspergillus flavus to oxidative stress in aflatoxin-conducive and non-conducive media

Aflatoxin production by isolates of Aspergillus flavus varies, ranging from highly toxigenic to completely atoxigenic. Several mechanisms have been identified which regulate aflatoxin production including medium carbon source and oxidative stress. In recent studies, aflatoxin production has been implicated in partially ameliorating oxidative stress in A. flavus. To better understand the role of aflatoxin production in oxidative stress responses, a selection of toxigenic and atoxigenic isolates of A. flavus with moderate to high oxidative stress tolerance were exposed to increasing concentrations of H2O2 in both aflatoxin-conducive and non-conducive media. Mycelial mats were collected for global transcriptome sequencing followed by differential expression, functional prediction, and weighted co-expression analyses. Oxidative stress and medium carbon source had a significant effect on the expression of several secondary metabolite gene clusters including those for aflatoxin, aflatrem, aflavarin, cyclopiazonic acid, and kojic acid. Atoxigenic biological control isolates showed less differential expression under stress than other atoxigenic isolates suggesting expression profiles may be useful in screening. Increasing stress also resulted in regulation of SakA/Hog1 and MpkA MAP kinase signalling pathways pointing to their potential roles in regulating oxidative stress responses. Their expression was also influenced by medium carbon source. These results suggest that aflatoxin production along with that of other mycotoxins may occur as part of a concerted coping mechanism for oxidative stress and its effects in the environment. This mechanism is also regulated by availability of simple sugars and glycolytic compounds for their biosynthesis.

Transcriptional responses of toxigenic and atoxigenic isolates of Aspergillus flavus to oxidative stress in aflatoxin-conducive and non-conducive media

Aflatoxin production by isolates of Aspergillus flavus varies, ranging from highly toxigenic to completely atoxigenic. Several mechanisms have been identified which regulate aflatoxin production including medium carbon source and oxidative stress. In recent studies, aflatoxin production has been implicated in partially ameliorating oxidative stress in A. flavus. To better understand the role of aflatoxin production in oxidative stress responses, a selection of toxigenic and atoxigenic isolates of A. flavus with moderate to high oxidative stress tolerance were exposed to increasing concentrations of H2O2 in both aflatoxin-conducive and non-conducive media. Mycelial mats were collected for global transcriptome sequencing followed by differential expression, functional prediction, and weighted co-expression analyses. Oxidative stress and medium carbon source had a significant effect on the expression of several secondary metabolite gene clusters including those for aflatoxin, aflatrem, aflavarin, cyclopiazonic acid, and kojic acid. Atoxigenic biological control isolates showed less differential expression under stress than other atoxigenic isolates suggesting expression profiles may be useful in screening. Increasing stress also resulted in regulation of SakA/Hog1 and MpkA MAP kinase signalling pathways pointing to their potential roles in regulating oxidative stress responses. Their expression was also influenced by medium carbon source. These results suggest that aflatoxin production along with that of other mycotoxins may occur as part of a concerted coping mechanism for oxidative stress and its effects in the environment. This mechanism is also regulated by availability of simple sugars and glycolytic compounds for their biosynthesis.

Arachis hypogaea gene expression atlas for fastigiata subspecies of cultivated groundnut to accelerate functional and translational genomics applications

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.

Genome-Wide SNP Discovery and Mapping QTLs for Seed Iron and Zinc Concentrations in Chickpea (Cicer arietinum L.)

Biofortification through plant breeding is a cost-effective and sustainable approach towards addressing micronutrient malnutrition prevailing across the globe. Screening cultivars for micronutrient content and identification of quantitative trait loci (QTLs)/genes and markers help in the development of biofortified varieties in chickpea (Cicer arietinum L.). With the aim of identifying the genomic regions controlling seed Fe and Zn concentrations, the F2:3 population derived from a cross between MNK-1 and Annigeri 1 was genotyped using genotyping by sequencing approach and evaluated for Fe and Zn concentration. An intraspecific genetic linkage map comprising 839 single nucleotide polymorphisms (SNPs) spanning a total distance of 1,088.04 cM with an average marker density of 1.30 cM was constructed. By integrating the linkage map data with the phenotypic data of the F2:3 population, a total of 11 QTLs were detected for seed Fe concentration on CaLG03, CaLG04, and CaLG05, with phenotypic variation explained ranging from 7.2% (CaqFe3.4) to 13.4% (CaqFe4.2). For seed Zn concentration, eight QTLs were identified on CaLG04, CaLG05, and CaLG08. The QTLs individually explained phenotypic variations ranging between 5.7% (CaqZn8.1) and 13.7% (CaqZn4.3). Three QTLs for seed Fe and Zn concentrations (CaqFe4.4, CaqFe4.5, and CaqZn4.1) were colocated in the "QTL-hotspot" region on CaLG04 that harbors several drought tolerance-related QTLs. We identified genes in the QTL regions that encode iron-sulfur metabolism and zinc-dependent alcohol dehydrogenase activity on CaLG03, iron ion binding oxidoreductase on CaLG04, and zinc-induced facilitator-like protein and ZIP zinc/iron transport family protein on CaLG05. These genomic regions and the associated markers can be used in marker-assisted selection to increase seed Fe and Zn concentrations in agronomically superior chickpea varieties.

Unraveling candidate genomic regions responsible for delayed leaf senescence in rice

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 BC₂F₄ 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).

Molecular and Physiological Alterations in Chickpea under Elevated CO2 Concentrations

The present study reports profiling of the elevated carbon dioxide (CO2) concentration responsive global transcriptome in chickpea, along with a combinatorial approach for exploring interlinks between physiological and transcriptional changes, important for the climate change scenario. Various physiological parameters were recorded in two chickpea cultivars (JG 11 and KAK 2) grown in open top chambers under ambient [380 parts per million (ppm)] and two stressed/elevated CO2 concentrations (550 and 700 ppm), at different stages of plant growth. The elevated CO2 concentrations altered shoot and root length, nodulation (number of nodules), total chlorophyll content and nitrogen balance index, significantly. RNA-Seq from 12 tissues representing vegetative and reproductive growth stages of both cultivars under ambient and elevated CO2 concentrations identified 18,644 differentially expressed genes including 9,687 transcription factors (TF). The differential regulations in genes, gene networks and quantitative real-time polymerase chain reaction (qRT-PCR) -derived expression dynamics of stress-responsive TFs were observed in both cultivars studied. A total of 138 pathways, mainly involved in sugar/starch metabolism, chlorophyll and secondary metabolites biosynthesis, deciphered the crosstalk operating behind the responses of chickpea to elevated CO2 concentration.

Genome-wide analysis of epigenetic and transcriptional changes associated with heterosis in pigeonpea

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.

Multiomics approach unravels fertility transition in a pigeonpea line for a two-line hybrid system

Pigeonpea [Cajanus cajan (L.) Millsp.] is a pulse crop cultivated in the semi-arid regions of Asia and Africa. It is a rich source of protein and capable of alleviating malnutrition, improving soil health and the livelihoods of small-holder farmers. Hybrid breeding has provided remarkable improvements for pigeonpea productivity, but owing to a tedious and costly seed production system, an alternative two-line hybrid technology is being explored. In this regard, an environment-sensitive male sterile line has been characterized as a thermosensitive male sterile line in pigeonpea precisely responding to day temperature. The male sterile and fertile anthers from five developmental stages were studied by integrating transcriptomics, proteomics and metabolomics supported by precise phenotyping and scanning electron microscopic study. Spatio-temporal analysis of anther transcriptome and proteome revealed 17 repressed DEGs/DEPs in sterile anthers that play a critical role in normal cell wall morphogenesis and tapetal cell development. The male fertility to sterility transition was mainly due to a perturbation in auxin homeostasis, leading to impaired cell wall modification and sugar transport. Limited nutrient utilization thus leads to microspore starvation in response to moderately elevated day temperature which could be restored with auxin-treatment in the male sterile line. Our findings outline a molecular mechanism that underpins fertility transition responses thereby providing a process-oriented two-line hybrid breeding framework for pigeonpea.

Complete genome sequence of sixteen plant growth promoting Streptomyces strains

The genome sequences of 16 Streptomyces strains, showing potential for plant growth-promotion (PGP) activities in rice, sorghum, chickpea and pigeonpea, isolated from herbal vermicompost, have been decoded. The genome assemblies of the 16 Streptomyces strains ranged from 6.8 Mb to 8.31 Mb, with a GC content of 72 to 73%. The extent of sequence similarity (in terms of shared ortholog) in 16 Streptomyces strains showed 70 to 85% common genes to the closest publicly available Streptomyces genomes. It was possible to identify ~1,850 molecular functions across these 16 strains, of which close to 50% were conserved across the genomes of Streptomyces strains, whereas, ~10% were strain specific and the rest were present in various combinations. Genome assemblies of the 16 Streptomyces strains have also provided genes involved in key pathways related to PGP and biocontrol traits such as siderophores, auxin, hydrocyanic acid, chitinase and cellulase. Further, the genome assemblies provided better understanding of genetic similarity among target strains and with the publically available Streptomyces strains.

The genome sequence of segmental allotetraploid peanut Arachis hypogaea

Like many other crops, the cultivated peanut (Arachis hypogaea L.) is of hybrid origin and has a polyploid genome that contains essentially complete sets of chromosomes from two ancestral species. Here we report the genome sequence of peanut and show that after its polyploid origin, the genome has evolved through mobile-element activity, deletions and by the flow of genetic information between corresponding ancestral chromosomes (that is, homeologous recombination). Uniformity of patterns of homeologous recombination at the ends of chromosomes favors a single origin for cultivated peanut and its wild counterpart A. monticola. However, through much of the genome, homeologous recombination has created diversity. Using new polyploid hybrids made from the ancestral species, we show how this can generate phenotypic changes such as spontaneous changes in the color of the flowers. We suggest that diversity generated by these genetic mechanisms helped to favor the domestication of the polyploid A. hypogaea over other diploid Arachis species cultivated by humans.

Integrated transcriptome, small RNA and degradome sequencing approaches provide insights into Ascochyta blight resistance in chickpea

Ascochyta blight (AB) is one of the major biotic stresses known to limit the chickpea production worldwide. To dissect the complex mechanisms of AB resistance in chickpea, three approaches, namely, transcriptome, small RNA and degradome sequencing were used. The transcriptome sequencing of 20 samples including two resistant genotypes, two susceptible genotypes and one introgression line under control and stress conditions at two time points (3rd and 7th day post inoculation) identified a total of 6767 differentially expressed genes (DEGs). These DEGs were mainly related to pathogenesis-related proteins, disease resistance genes like NBS-LRR, cell wall biosynthesis and various secondary metabolite synthesis genes. The small RNA sequencing of the samples resulted in the identification of 651 miRNAs which included 478 known and 173 novel miRNAs. A total of 297 miRNAs were differentially expressed between different genotypes, conditions and time points. Using degradome sequencing and in silico approaches, 2131 targets were predicted for 629 miRNAs. The combined analysis of both small RNA and transcriptome datasets identified 12 miRNA-mRNA interaction pairs that exhibited contrasting expression in resistant and susceptible genotypes and also, a subset of genes that might be post-transcriptionally silenced during AB infection. The comprehensive integrated analysis in the study provides better insights into the transcriptome dynamics and regulatory network components associated with AB stress in chickpea and, also offers candidate genes for chickpea improvement.

Genotyping-by-sequencing based genetic mapping reveals large number of epistatic interactions for stem rot resistance in groundnut

Genetic mapping identified large number of epistatic interactions indicating the complex genetic architecture for stem rot disease resistance. Groundnut (Arachis hypogaea) is an important global crop commodity and serves as a major source of cooking oil, diverse confectionery preparations and livestock feed. Stem rot disease caused by Sclerotium rolfsii is the most devastating disease of groundnut and can cause up to 100% yield loss. Genomic-assisted breeding (GAB) has potential for accelerated development of stem rot resistance varieties in short period with more precision. In this context, linkage analysis and quantitative trait locus (QTL) mapping for resistance to stem rot disease was performed in a bi-parental recombinant inbred line population developed from TG37A (susceptible) × NRCG-CS85 (resistant) comprising of 270 individuals. Genotyping-by-sequencing approach was deployed to generate single nucleotide polymorphism (SNP) genotyping data leading to development of a genetic map with 585 SNP loci spanning map distance of 2430 cM. QTL analysis using multi-season phenotyping and genotyping data could not detect any major main-effect QTL but identified 44 major epistatic QTLs with phenotypic variation explained ranging from 14.32 to 67.95%. Large number interactions indicate the complexity of genetic architecture of resistance to stem rot disease. A QTL of physical map length 5.2 Mb identified on B04 comprising 170 different genes especially leucine reach repeats, zinc finger motifs and ethyleneresponsive factors, etc., was identified. The identified genomic regions and candidate genes will further validate and facilitate marker development to deploy GAB for developing stem rot disease resistance groundnut varieties.

Genome sequence of Jatropha curcas L., a non-edible biodiesel plant, provides a resource to improve seed-related traits

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.

Development and Application of High-Density Axiom Cajanus SNP Array with 56K SNPs to Understand the Genome Architecture of Released Cultivars and Founder Genotypes

As one of the major outputs of next-generation sequencing (NGS), a large number of genome-wide single-nucleotide polymorphisms (SNPs) have been developed in pigeonpea [ (L.) Huth.]. However, SNPs require a genotyping platform or assay to be used in different evolutionary studies or in crop improvement programs. Therefore, we developed an Axiom SNP array with 56K SNPs uniformly distributed across the genome and assessed its utility in a genetic diversity study. From the whole-genome resequencing (WGRS) data on 104 pigeonpea lines, ∼2 million sequence variations (SNPs and insertion-deletions [InDels]) were identified, from which a subset of 56,512 unique and informative sequence variations were selected to develop the array. The Axiom SNP array developed was used for genotyping 103 pigeonpea lines encompassing 63 cultivars released between 1960 and 2014 and 40 breeding, germplasm, and founder lines. Genotyping data thus generated on 103 pigeonpea lines provided 51,201 polymorphic SNPs and InDels. Genetic diversity analysis provided in-depth insights into the genetic architecture and trends in temporal diversity in pigeonpea cultivars. Therefore, the continuous use of the high-density Axiom SNP array developed will accelerate high-resolution trait mapping, marker-assisted breeding, and genomic selection efforts in pigeonpea.

The RNA-Seq-based high resolution gene expression atlas of chickpea (Cicer arietinum L.) reveals dynamic spatio-temporal changes associated with growth and development

Chickpea is one of the world's largest cultivated food legumes and is an excellent source of high-quality protein to the human diet. Plant growth and development are controlled by programmed expression of a suite of genes at the given time, stage, and tissue. Understanding how the underlying genome sequence translates into specific plant phenotypes at key developmental stages, information on gene expression patterns is crucial. Here, we present a comprehensive Cicer arietinum Gene Expression Atlas (CaGEA) across different plant developmental stages and organs covering the entire life cycle of chickpea. One of the widely used drought tolerant cultivars, ICC 4958 has been used to generate RNA-Seq data from 27 samples at 5 major developmental stages of the plant. A total of 816 million raw reads were generated and of these, 794 million filtered reads after quality control (QC) were subjected to downstream analysis. A total of 15,947 unique number of differentially expressed genes across different pairwise tissue combinations were identified. Significant differences in gene expression patterns contributing in the process of flowering, nodulation, and seed and root development were inferred in this study. Furthermore, differentially expressed candidate genes from "QTL-hotspot" region associated with drought stress response in chickpea were validated.