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.

Fine mapping and gene cloning in the post-NGS era: advances and prospects.

Improvement in traits of agronomic importance is the top breeding priority of crop improvement programs. Majority of these agronomic traits show complex quantitative inheritance. Identification of quantitative trait loci (QTLs) followed by fine mapping QTLs and cloning of candidate genes/QTLs is central to trait analysis. Advances in genomic technologies revolutionized our understanding of genetics of complex traits, and genomic regions associated with traits were employed in marker-assisted breeding or cloning of QTLs/genes. Next-generation sequencing (NGS) technologies have enabled genome-wide methodologies for the development of ultra-high-density genetic linkage maps in different crops, thus allowing placement of candidate loci within few kbs in genomes. In this review, we compare the marker systems used for fine mapping and QTL cloning in the pre- and post-NGS era. We then discuss how different NGS platforms in combination with advanced experimental designs have improved trait analysis and fine mapping. We opine that efficient genotyping/sequencing assays may circumvent the need for cumbersome procedures that were earlier used for fine mapping. A deeper understanding of the trait architectures of agricultural significance will be crucial to accelerate crop improvement.

Plant vigour QTLs co-map with an earlier reported QTL hotspot for drought tolerance while water saving QTLs map in other regions of the chickpea genome.

BACKGROUND: Terminal drought stress leads to substantial annual yield losses in chickpea (Cicer arietinum L.). Adaptation to water limitation is a matter of matching water supply to water demand by the crop. Therefore, harnessing the genetics of traits contributing to plant water use, i.e. transpiration rate and canopy development dynamics, is important to design crop ideotypes suited to a varying range of water limited environments. With an aim of identifying genomic regions for plant vigour (growth and canopy size) and canopy conductance traits, 232 recombinant inbred lines derived from a cross between ICC 4958 and ICC 1882, were phenotyped at vegetative stage under well-watered conditions using a high throughput phenotyping platform (LeasyScan). RESULTS: Twenty one major quantitative trait loci (M-QTLs) were identified for plant vigour and canopy conductance traits using an ultra-high density bin map. Plant vigour traits had 13 M-QTLs on CaLG04, with favourable alleles from high vigour parent ICC 4958. Most of them co-mapped with a previously fine mapped major drought tolerance "QTL-hotspot" region on CaLG04. One M-QTL was found for canopy conductance on CaLG03 with the ultra-high density bin map. Comparative analysis of the QTLs found across different density genetic maps revealed that QTL size reduced considerably and % of phenotypic variation increased as marker density increased. CONCLUSION: Earlier reported drought tolerance hotspot is a vigour locus. The fact that canopy conductance traits, i.e. the other important determinant of plant water use, mapped on CaLG03 provides an opportunity to manipulate these loci to tailor recombinants having low/high transpiration rate and plant vigour, fitted to specific drought stress scenarios in chickpea.

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

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

Genotyping-by-sequencing based intra-specific genetic map refines a ''QTL-hotspot" region for drought tolerance in chickpea.

To enhance the marker density in the "QTL-hotspot" region, harboring several QTLs for drought tolerance-related traits identified on linkage group 04 (CaLG04) in chickpea recombinant inbred line (RIL) mapping population ICC 4958 × ICC 1882, a genotyping-by-sequencing approach was adopted. In total, 6.24 Gb data from ICC 4958, 5.65 Gb data from ICC 1882 and 59.03 Gb data from RILs were generated, which identified 828 novel single-nucleotide polymorphisms (SNPs) for genetic mapping. Together with these new markers, a high-density intra-specific genetic map was developed that comprised 1,007 marker loci spanning a distance of 727.29 cM. QTL analysis using the extended genetic map along with precise phenotyping data for 20 traits collected over one to seven seasons identified 49 SNP markers in the "QTL-hotspot" region. These efforts have refined the "QTL-hotspot" region to 14 cM. In total, 164 main-effect QTLs including 24 novel QTLs were identified. In addition, 49 SNPs integrated in the "QTL-hotspot" region were converted into cleaved amplified polymorphic sequence (CAPS) and derived CAPS (dCAPS) markers which can be used in marker-assisted breeding.

Integrated physical, genetic and genome map of chickpea (Cicer arietinum L.).

Physical map of chickpea was developed for the reference chickpea genotype (ICC 4958) using bacterial artificial chromosome (BAC) libraries targeting 71,094 clones (~12× coverage). High information content fingerprinting (HICF) of these clones gave high-quality fingerprinting data for 67,483 clones, and 1,174 contigs comprising 46,112 clones and 3,256 singletons were defined. In brief, 574 Mb genome size was assembled in 1,174 contigs with an average of 0.49 Mb per contig and 3,256 singletons represent 407 Mb genome. The physical map was linked with two genetic maps with the help of 245 BAC-end sequence (BES)-derived simple sequence repeat (SSR) markers. This allowed locating some of the BACs in the vicinity of some important quantitative trait loci (QTLs) for drought tolerance and reistance to Fusarium wilt and Ascochyta blight. In addition, fingerprinted contig (FPC) assembly was also integrated with the draft genome sequence of chickpea. As a result, ~965 BACs including 163 minimum tilling path (MTP) clones could be mapped on eight pseudo-molecules of chickpea forming 491 hypothetical contigs representing 54,013,992 bp (~54 Mb) of the draft genome. Comprehensive analysis of markers in abiotic and biotic stress tolerance QTL regions led to identification of 654, 306 and 23 genes in drought tolerance "QTL-hotspot" region, Ascochyta blight resistance QTL region and Fusarium wilt resistance QTL region, respectively. Integrated physical, genetic and genome map should provide a foundation for cloning and isolation of QTLs/genes for molecular dissection of traits as well as markers for molecular breeding for chickpea improvement.

Genetic dissection of drought tolerance in chickpea (Cicer arietinum L.).

KEY MESSAGE: Analysis of phenotypic data for 20 drought tolerance traits in 1-7 seasons at 1-5 locations together with genetic mapping data for two mapping populations provided 9 QTL clusters of which one present on CaLG04 has a high potential to enhance drought tolerance in chickpea improvement. Chickpea (Cicer arietinum L.) is the second most important grain legume cultivated by resource poor farmers in the arid and semi-arid regions of the world. Drought is one of the major constraints leading up to 50% production losses in chickpea. In order to dissect the complex nature of drought tolerance and to use genomics tools for enhancing yield of chickpea under drought conditions, two mapping populations-ICCRIL03 (ICC 4958 × ICC 1882) and ICCRIL04 (ICC 283 × ICC 8261) segregating for drought tolerance-related root traits were phenotyped for a total of 20 drought component traits in 1-7 seasons at 1-5 locations in India. Individual genetic maps comprising 241 loci and 168 loci for ICCRIL03 and ICCRIL04, respectively, and a consensus genetic map comprising 352 loci were constructed ( http://cmap.icrisat.ac.in/cmap/sm/cp/varshney/). Analysis of extensive genotypic and precise phenotypic data revealed 45 robust main-effect QTLs (M-QTLs) explaining up to 58.20% phenotypic variation and 973 epistatic QTLs (E-QTLs) explaining up to 92.19% phenotypic variation for several target traits. Nine QTL clusters containing QTLs for several drought tolerance traits have been identified that can be targeted for molecular breeding. Among these clusters, one cluster harboring 48% robust M-QTLs for 12 traits and explaining about 58.20% phenotypic variation present on CaLG04 has been referred as "QTL-hotspot". This genomic region contains seven SSR markers (ICCM0249, NCPGR127, TAA170, NCPGR21, TR11, GA24 and STMS11). Introgression of this region into elite cultivars is expected to enhance drought tolerance in chickpea.

A conserved SNP variation in the pre-miR396c flanking region in Oryza sativa indica landraces correlates with mature miRNA abundance.

Plant precursor miRNAs (pre-miRNA) have conserved evolutionary footprints that correlate with mode of miRNA biogenesis. In plants, base to loop and loop to base modes of biogenesis have been reported. Conserved structural element(s) in pre-miRNA play a major role in turn over and abundance of mature miRNA. Pre-miR396c sequences and secondary structural characteristics across Oryza species are presented. Based on secondary structure, twelve Oryza pre-miR396c sequences are divided into three groups, with the precursor from halophytic Oryza coarctata forming a distinct group. The miRNA-miRNA* duplex region is completely conserved across eleven Oryza species as are other structural elements in the pre-miRNA, suggestive of an evolutionarily conserved base-to-loop mode of miRNA biogenesis. SNPs within O. coarctata mature miR396c sequence and miRNA* region have the potential to alter target specificity and association with the RNA-induced silencing complex. A conserved SNP variation, rs10234287911 (G/A), identified in O. sativa pre-miR396c sequences alters base pairing above the miRNA-miRNA* duplex. The more stable structure conferred by the 'A10234287911' allele may promote better processing vis-à-vis the structure conferred by 'G10234287911' allele. We also examine pri- and pre-miR396c expression in cultivated rice under heat and salinity and their correlation with miR396c expression.

Morpho-physiological, biochemical and molecular characterization of coastal rice landraces to identify novel genetic sources of salinity tolerance.

Soil salinity is a leading cause for yield losses in rice, affecting nearly 6% of global rice cultivable area. India is host to a rich diversity of coastal rice landraces that are naturally tolerant to salinity and an untapped source to identify novel determinants of salinity tolerance. In the present study, we have assessed the relative salinity tolerance of 43 previously genotyped rice landraces at seedling stage, using thirteen morpho-physiological and biochemical parameters using a hydroponics system. Among 43 rice varieties, 25 were tolerant, 15 were moderately tolerant, 1 was moderately susceptible and 2 sensitive checks were found to be highly susceptible based on standard salinity scoring methods. In addition to previously known saline tolerant genotypes (Pokkali, FL478 and Nona Bokra), the present study has novel genotypes such as Katrangi, Orkyma, Aduisen 1, Orumundakan 1, Hoogla, and Talmugur 2 as potential sources of salinity tolerance through measurement of morpho-physiological and biochemical parameters including Na+, K+ estimations and Na+/K+ ratios. Further, Pallipuram Pokkali may be an important source of the tissue tolerance trait under salinity. Four marker trait associations (RM455-root Na+; RM161-shoot and root Na+/K+ ratios; RM237-salinity tolerance index) accounted for phenotypic variations in the range of 20.97-39.82%. A significant increase in root endodermal and exodermal suberization was observed in selected rice landraces under salinity. For the first time, variation in the number of suberized sclerenchymatous layers as well as passage cells is reported, in addition to expression level changes in suberin biosynthetic genes (CYP86A2, CYP81B1, CYP86A8 and PERL).

Creation of novel alleles of fragrance gene OsBADH2 in rice through CRISPR/Cas9 mediated gene editing.

Fragrance in rice grains is a key quality trait determining its acceptability and marketability. Intensive research on rice aroma identified mutations in betaine aldehyde dehydrogenase (OsBADH2) leading to production of aroma in rice. Gene editing technologies like CRISPR/Cas9 system has opened new avenues for accelerated improvement of rice grain quality through targeted mutagenesis. In this study, we have employed CRISPR/Cas9 tool to create novel alleles of OsBADH2 leading to introduction of aroma into an elite non-aromatic rice variety ASD16. PCR analysis of putative transformants using primers targeting the flanking regions of sgRNA in the 7th exon of OsBADH2 identified 37.5% potential multi-allelic mutations in T0 generation. Sensory evaluation test in the leaves of T0 lines identified thirteen lines belonging to five independent events producing aroma. Sequence analysis of these aromatic T0 lines identified 22 different types of mutations located within -17 bp to +15bp of sgRNA region. The -1/-2 bp deletion in the line # 8-19 and -8/-5 bp deletion in the line # 2-16 produced strong aroma and the phenotype was stably inherited in the T1 generation. Comparative volatile profiling detected novel aromatic compounds viz., pyrrolidine, pyridine, pyrazine, pyradazine and pyrozole in the grains of T1 progenies of line # 8-19. This study has demonstrated the use of CRISPR/Cas9 in creating novel alleles of OsBADH2 to introduce aroma into any non-aromatic rice varieties.

Analysis of genetic diversity and population structure using SSR markers and validation of a Cleavage Amplified Polymorphic Sequences (CAPS) marker involving the sodium transporter OsHKT1;5 in saline tolerant rice (Oryza sativa L.) landraces.

Naturally evolved saline tolerant rice landraces found along the coastline of India are a valuable genomic resource to explore the complex, polygenic nature of salinity tolerance. In the present study, a set of 28 genome wide SSR markers, 11 salt responsive genic SSR markers and 8 Saltol QTL linked SSR markers were used to estimate genetic relatedness and population structure within a collection of 47 rice landraces (including a tolerant and 2 sensitive checks) originating from geographically divergent coastal regions of India. All three marker types identified substantial genetic variation among the landraces, as evident from their higher PIC values (0.53 for genomic SSRs, 0.43 for Genic SSRs and 0.59 for Saltol SSRs). The markers RM431, RM484 (Genomic SSRs), OsCAX (D), OsCAX (T) (Genic SSRs) and RM562 (Saltol SSR) were identified as good candidates to be used in breeding programs for improving salinity tolerance in rice. STRUCTURE analysis divided the landraces into five distinct populations, with classification correlating with their geographical locations. Principal coordinate and hierarchical cluster analyses (UPGMA and neighbor joining) are in close agreement with STRUCTURE results. AMOVA analysis indicated a higher magnitude of genetic differentiation within individuals of groups (58%), than among groups (42%). We also report the development and validation of a new Cleavage Amplified Polymorphic Sequence (CAPS) marker (OsHKT1;5V395) that targets a codon in the sodium transporter gene OsHKT1;5 (Saltol/SKC1 locus) that is associated with sodium transport rates in the above rice landraces. The CAPS marker was found to be present in all landraces except in IR29, Kamini, Gheus, Matla 1 and Matla 2. Significant molecular genetic diversity established among the analyzed salt tolerant rice landraces will aid in future association mapping; the CAPS marker, OsHKT1;5V395 can be used to map rice landraces for the presence of the SNP (Single Nucleotide Polymorphism) associated with increased sodium transport rates and concomitant salinity tolerance in rice.

CRISPR for Crop Improvement: An Update Review.

The availability of genome sequences for several crops and advances in genome editing approaches has opened up possibilities to breed for almost any given desirable trait. Advancements in genome editing technologies such as zinc finger nucleases (ZFNs), transcription activator-like effector nucleases (TALENs) has made it possible for molecular biologists to more precisely target any gene of interest. However, these methodologies are expensive and time-consuming as they involve complicated steps that require protein engineering. Unlike first-generation genome editing tools, CRISPR/Cas9 genome editing involves simple designing and cloning methods, with the same Cas9 being potentially available for use with different guide RNAs targeting multiple sites in the genome. After proof-of-concept demonstrations in crop plants involving the primary CRISPR-Cas9 module, several modified Cas9 cassettes have been utilized in crop plants for improving target specificity and reducing off-target cleavage (e.g., Nmcas9, Sacas9, and Stcas9). Further, the availability of Cas9 enzymes from additional bacterial species has made available options to enhance specificity and efficiency of gene editing methodologies. This review summarizes the options available to plant biotechnologists to bring about crop improvement using CRISPR/Cas9 based genome editing tools and also presents studies where CRISPR/Cas9 has been used for enhancing biotic and abiotic stress tolerance. Application of these techniques will result in the development of non-genetically modified (Non-GMO) crops with the desired trait that can contribute to increased yield potential under biotic and abiotic stress conditions.

Prioritization of candidate genes in "QTL-hotspot" region for drought tolerance in chickpea (Cicer arietinum L.).

A combination of two approaches, namely QTL analysis and gene enrichment analysis were used to identify candidate genes in the "QTL-hotspot" region for drought tolerance present on the Ca4 pseudomolecule in chickpea. In the first approach, a high-density bin map was developed using 53,223 single nucleotide polymorphisms (SNPs) identified in the recombinant inbred line (RIL) population of ICC 4958 (drought tolerant) and ICC 1882 (drought sensitive) cross. QTL analysis using recombination bins as markers along with the phenotyping data for 17 drought tolerance related traits obtained over 1-5 seasons and 1-5 locations split the "QTL-hotspot" region into two subregions namely "QTL-hotspot_a" (15 genes) and "QTL-hotspot_b" (11 genes). In the second approach, gene enrichment analysis using significant marker trait associations based on SNPs from the Ca4 pseudomolecule with the above mentioned phenotyping data, and the candidate genes from the refined "QTL-hotspot" region showed enrichment for 23 genes. Twelve genes were found common in both approaches. Functional validation using quantitative real-time PCR (qRT-PCR) indicated four promising candidate genes having functional implications on the effect of "QTL-hotspot" for drought tolerance in chickpea.