Aspergillus flavus pangenome (AflaPan) uncovers novel aflatoxin and secondary metabolite associated gene clusters.

BACKGROUND: Aspergillus flavus is an important agricultural and food safety threat due to its production of carcinogenic aflatoxins. It has high level of genetic diversity that is adapted to various environments. Recently, we reported two reference genomes of A. flavus isolates, AF13 (MAT1-2 and highly aflatoxigenic isolate) and NRRL3357 (MAT1-1 and moderate aflatoxin producer). Where, an insertion of 310 kb in AF13 included an aflatoxin producing gene bZIP transcription factor, named atfC. Observations of significant genomic variants between these isolates of contrasting phenotypes prompted an investigation into variation among other agricultural isolates of A. flavus with the goal of discovering novel genes potentially associated with aflatoxin production regulation. Present study was designed with three main objectives: (1) collection of large number of A. flavus isolates from diverse sources including maize plants and field soils; (2) whole genome sequencing of collected isolates and development of a pangenome; and (3) pangenome-wide association study (Pan-GWAS) to identify novel secondary metabolite cluster genes. RESULTS: Pangenome analysis of 346 A. flavus isolates identified a total of 17,855 unique orthologous gene clusters, with mere 41% (7,315) core genes and 59% (10,540) accessory genes indicating accumulation of high genomic diversity during domestication. 5,994 orthologous gene clusters in accessory genome not annotated in either the A. flavus AF13 or NRRL3357 reference genomes. Pan-genome wide association analysis of the genomic variations identified 391 significant associated pan-genes associated with aflatoxin production. Interestingly, most of the significantly associated pan-genes (94%; 369 associations) belonged to accessory genome indicating that genome expansion has resulted in the incorporation of new genes associated with aflatoxin and other secondary metabolites. CONCLUSION: In summary, this study provides complete pangenome framework for the species of Aspergillus flavus along with associated genes for pathogen survival and aflatoxin production. The large accessory genome indicated large genome diversity in the species A. flavus, however AflaPan is a closed pangenome represents optimum diversity of species A. flavus. Most importantly, the newly identified aflatoxin producing gene clusters will be a new source for seeking aflatoxin mitigation strategies and needs new attention in research.

Double-digest restriction-associated DNA sequencing-based genotyping and its applications in sesame germplasm management.

Sesame (Sesamum indicum L.) is an ancient oilseed crop belonging to the family Pedaliaceae and a globally cultivated crop for its use as oil and food. In this study, 2496 sesame accessions, being conserved at the National Genebank of ICAR-National Bureau of Plant Genetic Resources (NBPGR), were genotyped using genomics-assisted double-digest restriction-associated DNA sequencing (ddRAD-seq) approach. A total of 64,910 filtered single-nucleotide polymorphisms (SNPs) were utilized to assess the genome-scale diversity. Applications of this genome-scale information (reduced representation using restriction enzymes) are demonstrated through the development of a molecular core collection (CC) representing maximal SNP diversity. This information is also applied in developing a mid-density panel (MDP) comprising 2515 hyper-variable SNPs, representing almost equally the genic and non-genic regions. The sesame CC comprising 384 accessions, a representative set of accessions with maximal diversity, was identified using multiple criteria such as k-mer (subsequence of length "k" in a sequence read) diversity, observed heterozygosity, CoreHunter3, GenoCore, and genetic differentiation. The coreset constituted around 15% of the total accessions studied, and this small subset had captured >60% SNP diversity of the entire population. In the coreset, the admixture analysis shows reduced genetic complexity, increased nucleotide diversity (π), and is geographically distributed without any repetitiveness in the CC germplasm. Within the CC, India-originated accessions exhibit higher diversity (as expected based on the center of diversity concept), than those accessions that were procured from various other countries. The identified CC set and the MDP will be a valuable resource for genomics-assisted accelerated sesame improvement program.

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.

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.

Genetic mapping identified three hotspot genomic regions and candidate genes controlling heat tolerance-related traits in groundnut.

Groundnut productivity and quality have been impeded by rising temperatures in semi-arid environments. Hence, understanding the effects and molecular mechanisms of heat stress tolerance will aid in tackling yield losses. In this context, a recombinant inbred line (RIL) population was developed and phenotyped for eight seasons at three locations for agronomic, phenological, and physiological traits under heat stress. A genetic map was constructed using genotyping-by-sequencing with 478 single-nucleotide polymorphism (SNP) loci spanning a map distance of 1,961.39 cM. Quantitative trait locus (QTL) analysis using phenotypic and genotypic data identified 45 major main-effect QTLs for 21 traits. Intriguingly, three QTL clusters (Cluster-1-Ah03, Cluster-2-Ah12, and Cluster-3-Ah20) harbor more than half of the major QTLs (30/45, 66.6%) for various heat tolerant traits, explaining 10.4%-38.6%, 10.6%-44.6%, and 10.1%-49.5% of phenotypic variance, respectively. Furthermore, important candidate genes encoding DHHC-type zinc finger family protein (arahy.J0Y6Y5), peptide transporter 1 (arahy.8ZMT0C), pentatricopeptide repeat-containing protein (arahy.4A4JE9), Ulp1 protease family (arahy.X568GS), Kelch repeat F-box protein (arahy.I7X4PC), FRIGIDA-like protein (arahy.0C3V8Z), and post-illumination chlorophyll fluorescence increase (arahy.92ZGJC) were the underlying three QTL clusters. The putative functions of these genes suggested their involvement in seed development, regulating plant architecture, yield, genesis and growth of plants, flowering time regulation, and photosynthesis. Our results could provide a platform for further fine mapping, gene discovery, and developing markers for genomics-assisted breeding to develop heat-tolerant groundnut varieties.

A pilot-scale comparison between single and double-digest RAD markers generated using GBS strategy in sesame (Sesamum indicum L.).

To reduce the genome sequence representation, restriction site-associated DNA sequencing (RAD-seq) protocols is being widely used either with single-digest or double-digest methods. In this study, we genotyped the sesame population (48 sample size) in a pilot scale to compare single and double-digest RAD-seq (sd and ddRAD-seq) methods. We analysed the resulting short-read data generated from both protocols and assessed their performance impacting the downstream analysis using various parameters. The distinct k-mer count and gene presence absence variation (PAV) showed a significant difference between the sesame samples studied. Additionally, the variant calling from both datasets (sdRAD-seq and ddRAD-seq) exhibits a significant difference between them. The combined variants from both datasets helped in identifying the most diverse samples and possible sub-groups in the sesame population. The most diverse samples identified from each analysis (k-mer, gene PAV, SNP count, Heterozygosity, NJ and PCA) can possibly be representative samples holding major diversity of the small sesame population used in this study. The best possible strategies with suggested inputs for modifications to utilize the RAD-seq strategy efficiently on a large dataset containing thousands of samples to be subjected to molecular analysis like diversity, population structure and core development studies were discussed.

Genome-wide association studies reveal novel loci for resistance to groundnut rosette disease in the African core groundnut collection.

KEY MESSAGE: We identified markers associated with GRD resistance after screening an Africa-wide core collection across three seasons in Uganda Groundnut is cultivated in several African countries where it is a major source of food, feed and income. One of the major constraints to groundnut production in Africa is groundnut rosette disease (GRD), which is caused by a complex of three agents: groundnut rosette assistor luteovirus, groundnut rosette umbravirus and its satellite RNA. Despite several years of breeding for GRD resistance, the genetics of the disease is not fully understood. The objective of the current study was to use the African core collection to establish the level of genetic variation in their response to GRD, and to map genomic regions responsible for the observed resistance. The African groundnut core genotypes were screened across two GRD hotspot locations in Uganda (Nakabango and Serere) for 3 seasons. The Area Under Disease Progress Curve combined with 7523 high quality SNPs were analyzed to establish marker-trait associations (MTAs). Genome-Wide Association Studies based on Enriched Compressed Mixed Linear Model detected 32 MTAs at Nakabango: 21 on chromosome A04, 10 on B04 and 1 on B08. Two of the significant markers were localised on the exons of a putative TIR-NBS-LRR disease resistance gene on chromosome A04. Our results suggest the likely involvement of major genes in the resistance to GRD but will need to be further validated with more comprehensive phenotypic and genotypic datasets. The markers identified in the current study will be developed into routine assays and validated for future genomics-assisted selection for GRD resistance in groundnut.

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.

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.

Genetic diversity and population structure of pigeonpea (Cajanus cajan [L.] Millspaugh) landraces grown in Benin revealed by Genotyping-By-Sequencing.

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

Root exudation of contrasting drought-stressed pearl millet genotypes conveys varying biological nitrification inhibition (BNI) activity.

Roots secrete a vast array of low molecular weight compounds into the soil broadly referred to as root exudates. It is a key mechanism by which plants and soil microbes interact in the rhizosphere. The effect of drought stress on the exudation process and composition is rarely studied, especially in cereal crops. This study focuses on comparative metabolic profiling of the exudates from sensitive and tolerant genotypes of pearl millet after a period of drought stress. We employed a combined platform of gas and liquid chromatography coupled to mass spectrometry to cover both primary and secondary metabolites. The results obtained demonstrate that both genotype and drought stress have a significant impact on the concentration and composition of root exudates. The complexity and function of these differential root exudates are discussed. To reveal the potential effect of root exudates on the soil microbial community after a period of drought stress, we also tested for biological nitrification inhibition (BNI) activity. The analysis revealed a genotype-dependent enhancement of BNI activity after a defined period of drought stress. In parallel, we observed a genotype-specific relation of elongated root growth and root exudation under drought stress. These data suggest that the drought stress-dependent change in root exudation can manipulate the microbial soil communities to adapt and survive under harsh conditions. Supplementary Information: The online version contains supplementary material available at 10.1007/s00374-021-01578-w.

A chickpea genetic variation map based on the sequencing of 3,366 genomes.

Zero hunger and good health could be realized by 2030 through effective conservation, characterization and utilization of germplasm resources1. So far, few chickpea (Cicer arietinum) germplasm accessions have been characterized at the genome sequence level2. Here we present a detailed map of variation in 3,171 cultivated and 195 wild accessions to provide publicly available resources for chickpea genomics research and breeding. We constructed a chickpea pan-genome to describe genomic diversity across cultivated chickpea and its wild progenitor accessions. A divergence tree using genes present in around 80% of individuals in one species allowed us to estimate the divergence of Cicer over the last 21 million years. Our analysis found chromosomal segments and genes that show signatures of selection during domestication, migration and improvement. The chromosomal locations of deleterious mutations responsible for limited genetic diversity and decreased fitness were identified in elite germplasm. We identified superior haplotypes for improvement-related traits in landraces that can be introgressed into elite breeding lines through haplotype-based breeding, and found targets for purging deleterious alleles through genomics-assisted breeding and/or gene editing. Finally, we propose three crop breeding strategies based on genomic prediction to enhance crop productivity for 16 traits while avoiding the erosion of genetic diversity through optimal contribution selection (OCS)-based pre-breeding. The predicted performance for 100-seed weight, an important yield-related trait, increased by up to 23% and 12% with OCS- and haplotype-based genomic approaches, respectively.

Transcriptome analysis reveals key genes associated with root-lesion nematode Pratylenchus thornei resistance in chickpea.

The root-lesion nematode, Pratylenchus thornei, is one of the major plant-parasitic nematode species causing significant yield losses in chickpea (Cicer arietinum). In order to identify the underlying mechanisms of resistance to P. thornei, the transcriptomes of control and inoculated roots of three chickpea genotypes viz. D05253 > F3TMWR2AB001 (resistant advanced breeding line), PBA HatTrick (moderately resistant cultivar), and Kyabra (susceptible cultivar) were studied at 20 and 50 days post inoculation using the RNA-seq approach. On analyzing the 633.3 million reads generated, 962 differentially expressed genes (DEGs) were identified. Comparative analysis revealed that the majority of DEGs upregulated in the resistant genotype were downregulated in the moderately resistant and susceptible genotypes. Transcription factor families WRKY and bZIP were uniquely expressed in the resistant genotype. The genes Cysteine-rich receptor-like protein kinase 10, Protein lifeguard-like, Protein detoxification, Bidirectional sugar transporter Sugars Will Eventually be Exported Transporters1 (SWEET1), and Subtilisin-like protease were found to play cross-functional roles in the resistant chickpea genotype against P. thornei. The identified candidate genes for resistance to P. thornei in chickpea can be explored further to develop markers and accelerate the introgression of P. thornei resistance into elite chickpea cultivars.

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.

Spatial distribution of proteins and metabolites in developing wheat grain and their differential regulatory response during the grain filling process.

Grain filling and grain development are essential biological processes in the plant's life cycle, eventually contributing to the final seed yield and quality in all cereal crops. Studies of how the different wheat (Triticum aestivum L.) grain components contribute to the overall development of the seed are very scarce. We performed a proteomics and metabolomics analysis in four different developing components of the wheat grain (seed coat, embryo, endosperm, and cavity fluid) to characterize molecular processes during early and late grain development. In-gel shotgun proteomics analysis at 12, 15, 20, and 26 days after anthesis (DAA) revealed 15 484 identified and quantified proteins, out of which 410 differentially expressed proteins were identified in the seed coat, 815 in the embryo, 372 in the endosperm, and 492 in the cavity fluid. The abundance of selected protein candidates revealed spatially and temporally resolved protein functions associated with development and grain filling. Multiple wheat protein isoforms involved in starch synthesis such as sucrose synthases, starch phosphorylase, granule-bound and soluble starch synthase, pyruvate phosphate dikinase, 14-3-3 proteins as well as sugar precursors undergo a major tissue-dependent change in abundance during wheat grain development suggesting an intimate interplay of starch biosynthesis control. Different isoforms of the protein disulfide isomerase family as well as glutamine levels, both involved in the glutenin macropolymer pattern, showed distinct spatial and temporal abundance, revealing their specific role as indicators of wheat gluten quality. Proteins binned into the functional category of cell growth/division and protein synthesis/degradation were more abundant in the early stages (12 and 15 DAA). At the metabolome level all tissues and especially the cavity fluid showed highly distinct metabolite profiles. The tissue-specific data are integrated with biochemical networks to generate a comprehensive map of molecular processes during grain filling and developmental processes.

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.

Global Transcriptome Profiling Identified Transcription Factors, Biological Process, and Associated Pathways for Pre-Harvest Aflatoxin Contamination in Groundnut.

Pre-harvest aflatoxin contamination (PAC) in groundnut is a serious quality concern globally, and drought stress before harvest further exacerbate its intensity, leading to the deterioration of produce quality. Understanding the host-pathogen interaction and identifying the candidate genes responsible for resistance to PAC will provide insights into the defense mechanism of the groundnut. In this context, about 971.63 million reads have been generated from 16 RNA samples under controlled and Aspergillus flavus infected conditions, from one susceptible and seven resistant genotypes. The RNA-seq analysis identified 45,336 genome-wide transcripts under control and infected conditions. This study identified 57 transcription factor (TF) families with major contributions from 6570 genes coding for bHLH (719), MYB-related (479), NAC (437), FAR1 family protein (320), and a few other families. In the host (groundnut), defense-related genes such as senescence-associated proteins, resveratrol synthase, seed linoleate, pathogenesis-related proteins, peroxidases, glutathione-S-transferases, chalcone synthase, ABA-responsive gene, and chitinases were found to be differentially expressed among resistant genotypes as compared to susceptible genotypes. This study also indicated the vital role of ABA-responsive ABR17, which co-regulates the genes of ABA responsive elements during drought stress, while providing resistance against A. flavus infection. It belongs to the PR-10 class and is also present in several plant-pathogen interactions.

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.