Variability in genome size of Trigonella foenum-graecum, Trigonella corniculata and Trigonella caerulea as estimated by flow cytometry indicates complex evolutionary history of fenugreek.

BACKGROUND: The determination of genome size is a fundamental step which provides a basis to initiate studies aimed at deciphering the genetic similarity of a species and to carry out other genomics based investigations. Fenugreek (Trigonella spp.) is an important spice crop which has numerous health promoting phytochemicals. Many species within this genus are known for their various health benefits owing to the presence of a wide diversity of important phytochemicals like diosgenin, trigonelline, fenugreekine, galactomannan, 4-hydroxy isoleucine, etc. It is a multipurpose crop being cultivated for food, animal feed and industrial purposes. Despite its importance, research on the genomics aspect of fenugreek remains scant. In the absence of sufficient genomic information, crop improvement in fenugreek is severely lagging. METHODS AND RESULTS: Estimation of genome size of a species is the preliminary step for initiation of any genomic studies and therefore in the present study we have estimated the genome size for fenugreek. Here, we have determined the genome sizes of three different Trigonella spp. namely T. foenum-graecum, T. corniculata and T. caerulea through flow cytometry (FC). The 2 C DNA content values were found to be 6.05 pg (T. foenum-graecum), 1.83 pg (T. corniculata) and 1.96 pg (T. caerulea). The genome size of T. foenum-graecum is approximately three times the genome size of T. corniculata and T. caerulea. This variation in genome size of more than three-fold indicates the level of genetic divergence among the three species, though within the same genus. CONCLUSIONS: The differences observed in the genome sizes of the three species provide conclusive evidence of their genetic divergence. Additionally, the information about the genome size would provide an impetus to the structural and functional genomics-based research in this crop.

Novel insights into drought-induced regulation of ribosomal genes through DNA methylation in chickpea.

Modifications within the epigenome of an organism in response to external environmental conditions allow it to withstand the hostile stress factors. Drought in chickpea is a severely limiting abiotic stress factor which is known to cause huge yield loss. To analyse the methylome of chickpea in response to drought stress conditions and how it affects gene expression, we performed whole-genome bisulfite sequencing (WGBS) and RNA-seq of two chickpea genotypes which contrast for drought tolerance. It was observed that the mCHH was most variable under drought stress and the drought tolerant (DT) genotype exhibited substantial genome-wide hypomethylation as compared to the drought sensitive (DS) genotype. Specifically, there was substantial difference in gene expression and methylation for the ribosomal genes for the tolerant and sensitive genotypes. The differential expression of these genes was in complete agreement with earlier reported transcriptomes in chickpea. Many of these genes were hypomethylated (q < 0.01) and downregulated under drought stress (p < 0.01) in the sensitive genotype. The gene RPS6 (ribosomal protein small subunit) was found to be downregulated and hypomethylated in the drought sensitive genotype which could possibly lead to reduced ribosomal biosynthesis. This study provides novel insights into regulation of drought-responsive genes in chickpea.

DNA methylation: an emerging paradigm of gene regulation under drought stress in plants.

Drought is an enormous threat to global crop production. In order to ensure food security for the burgeoning population, we must develop drought tolerant crop varieties. This necessitates the identification of drought-responsive genes and understanding the mechanisms involved in their regulation. DNA methylation is a widely studied mechanism of epigenetic regulation of gene expression, which is known to play vital role in conferring tolerance to various biotic and abiotic stress factors. The recent advances in next-generation sequencing (NGS) technologies, has allowed unprecedented access to genome-wide methylation marks, with single base resolution. The most important roles of DNA methylation have been studied in terms of gene body methylation (gbM), which is associated with regulation of both transcript abundance and its stability. The availability of mutants for the various genes encoding enzymes involved in methylation of DNA has allowed ascertainment of the biological significance of methylation. Even though a vast number of reports have emerged in the recent past, where both genome-wide methylation landscape and locus specific changes in DNA methylation have been studied, a conclusive picture with regards to the biological role of DNA methylation is still lacking. Compounding this, is the lack of sufficient evidence supporting the heritability of these epigenetic changes. Amongst the various epigenetic variations, the DNA methylation changes are observed to be the most stable. This review describes the drought-induced changes in DNA methylation identified across different plant species. We also briefly describe the stress memory contributed by these changes. The identification of heritable, drought-induced methylation marks would broaden the scope of crop improvement in the future.

Unravelling the due importance of pseudogenes and their resurrection in plants.

The complexities of a genome are underpinned to the vast expanses of the intergenic region, which constitutes ∼97-98% of the genome. This region is essentially composed of what is colloquially referred to as the "junk DNA" and is composed of various elements like transposons, repeats, pseudogenes, etc. The latter have long been considered as dead elements merely contributing to transcriptional noise in the genome. Many studies now describe the previously unknown regulatory functions of these genes. Recent advances in the Next-generation sequencing (NGS) technologies have allowed unprecedented access to these regions. With the availability of whole genome sequences of more than 788 different plant species in past 20 years, genome annotation has become feasible like never before. Different bioinformatic pipelines are available for the identification of pseudogenes. However, still little is known about their biological functions. The functional validation of these genes remains challenging and research in this area is still in infancy, particularly in plants. CRISPR/Cas-based genome editing could provide solutions to understand the biological roles of these genes by allowing creation of precise edits within these genes. The possibility of pseudogene reactivation or resurrection as has been demonstrated in a few studies might open new avenues of genetic manipulation to yield a desirable phenotype. This review aims at comprehensively summarizing the progress made with regards to the identification of pseudogenes and understanding their biological functions in plants.

Genetic mapping of quantitative trait loci associated with drought tolerance in chickpea (Cicer arietinum L.).

Elucidation of the genetic basis of drought tolerance is vital for genomics-assisted breeding of drought tolerant crop varieties. Here, we used genotyping-by-sequencing (GBS) to identify single nucleotide polymorphisms (SNPs) in recombinant inbred lines (RILs) derived from a cross between a drought tolerant chickpea variety, Pusa 362 and a drought sensitive variety, SBD 377. The GBS identified a total of 35,502 SNPs and subsequent filtering of these resulted in 3237 high-quality SNPs included in the eight linkage groups. Fifty-one percent of these SNPs were located in the genic regions distributed throughout the genome. The high density linkage map has total map length of 1069 cm with an average marker interval of 0.33 cm. The linkage map was used to identify 9 robust and consistent QTLs for four drought related traits viz. membrane stability index, relative water content, seed weight and yield under drought, with percent variance explained within the range of 6.29%-90.68% and LOD scores of 2.64 to 6.38, which were located on five of the eight linkage groups. A genomic region on LG 7 harbors quantitative trait loci (QTLs) explaining > 90% phenotypic variance for membrane stability index, and > 10% PVE for yield. This study also provides the first report of major QTLs for physiological traits such as membrane stability index and relative water content for drought stress in chickpea. A total of 369 putative candidate genes were identified in the 6.6 Mb genomic region spanning these QTLs. In-silico expression profiling based on the available transcriptome data revealed that 326 of these genes were differentially expressed under drought stress. KEGG analysis resulted in reduction of candidate genes from 369 to 99, revealing enrichment in various signaling pathways. Haplotype analysis confirmed 5 QTLs among the initially identified 9 QTLs. Two QTLs, qRWC1.1 and qYLD7.1, were chosen based on high SNP density. Candidate gene-based analysis revealed distinct haplotypes in qYLD7.1 associated with significant phenotypic differences, potentially linked to pathways for secondary metabolite biosynthesis. These identified candidate genes bolster defenses through flavonoids and phenylalanine-derived compounds, aiding UV protection, pathogen resistance, and plant structure.The study provides novel genomic regions and candidate genes which can be utilized in genomics-assisted breeding of superior drought tolerant chickpea cultivars.

Genome wide identification and characterization of the amino acid transporter (AAT) genes regulating seed protein content in chickpea (Cicer arietinum L.).

Amino acid transporters (AATs), besides, being a crucial component for nutrient partitioning system are also vital for growth and development of the plants and stress resilience. In order to understand the role of AAT genes in seed quality proteins, a comprehensive analysis of AAT gene family was carried out in chickpea leading to identification of 109 AAT genes, representing 10 subfamilies with random distribution across the chickpea genome. Several important stress responsive cis-regulatory elements like Myb, ABRE, ERE were detected in the promoter region of these CaAAT genes. Most of the genes belonging to the same sub-families shared the intron-exon distribution pattern owing to their conserved nature. Random distribution of these CaAAT genes was observed on plasma membrane, vacuolar membrane, Endoplasmic reticulum and Golgi membranes, which may be associated to distinct biochemical pathways. In total 92 out 109 CaAAT genes arise as result of duplication, among which segmental duplication was more prominent over tandem duplication. As expected, the phylogenetic tree was divided into 2 major clades, and further sub-divided into different sub-families. Among the 109 CaAAT genes, 25 were found to be interacting with 25 miRNAs, many miRNAs like miR156, miR159 and miR164 were interacting only with single AAT genes. Tissues specific expression pattern of many CaAAT genes was observed like CaAAP7 and CaAVT18 in nodules, CaAAP17, CaAVT5 and CaCAT9 in vegetative tissues while CaCAT10 and CaAAP23 in seed related tissues as per the expression analysis. Mature seed transcriptome data revealed that genotypes having high protein content (ICC 8397, ICC 13461) showed low CaAATs expression as compared to the genotypes having low protein content (FG 212, BG 3054). Amino acid profiling of these genotypes revealed a significant difference in amount of essential and non-essential amino acids, probably due to differential expression of CaAATs. Thus, the present study provides insights into the biological role of AAT genes in chickpea, which will facilitate their functional characterization and role in various developmental stages, stress responses and involvement in nutritional quality enhancement.

Small cardamom genome: development and utilization of microsatellite markers from a draft genome sequence of Elettaria cardamomum Maton.

Small cardamom (Elettaria cardamomum Maton), the queen of spices, is the third most expensive spice in the world after saffron and vanilla, valued highly for its aroma and taste. This perennial herbaceous plant is a native of coastal parts of Southern India and displays a significant amount of morphological diversity. Its genetic potential has not been exploited due to lack of genomic resources limiting our understanding of the genome and important metabolic pathways which give it the economic advantage in the spice industry. Here, we report upon the de novo assembled, draft whole genome sequence of cardamom variety, Njallani Green Gold. We used a hybrid assembly strategy using the reads from the Oxford Nanopore, Illumina and 10x Genomics GemCode sequencing chemistries. The assembled genome length was 1.06 Gb (gigabases) which is close to the estimated genome size of cardamom. More than 75% of the genome was captured in 8000 scaffolds with a N50 of 0.15 Mb. The genome appears to have a high repeat content and 68055 gene models were predicted. The genome is close to Musa species and displays an expansion and contraction in different gene families. The draft assembly was used for in silico mining of simple sequence repeats (SSRs). A total of 2,50,571 SSRs were identified of which 2,18,270 were perfect SSRs and 32,301 were compound SSRs. Among the perfect SSRs, trinucleotides were most abundant (1,25,329) and hexanucleotide repeats appear least (2,380). From the 2,50,571 SSRs mined, 2,27,808 primer pairs were designed based on flanking sequence information. Wet lab validation was performed for 246 SSR loci and based on their amplification profiles, 60 SSR markers were used for diversity analysis of a set of 60 diverse cardamom accessions. The average number of alleles detected per locus were 14.57 with a minimum of 4 and maximum of 30 alleles. Population structure analysis revealed the presence of high degree of admixtures which could primarily be due to cross-pollination prevalent in this species. The SSR markers identified would help in the development of gene or trait-linked markers which can be subsequently used for marker-assisted breeding for crop improvement in cardamom. The information on utilization of the SSR loci for generation of markers has been developed into a public database, 'cardamomSSRdb' that is freely available for use by the cardamom community.

The SPL transcription factor genes are potential targets for epigenetic regulation in response to drought stress in chickpea (C. arietinum L.).

BACKGROUND: Crop improvement for tolerance to various biotic and abiotic stress factors necessitates understanding the key gene regulatory mechanisms. One such mechanism of gene regulation involves changes in cytosine methylation at the gene body and flanking regulatory sequences. The present study was undertaken to identify genes which might be potential targets of drought-induced DNA methylation in chickpea. METHODS AND RESULTS: Two chickpea genotypes, which contrast for drought tolerance, were subjected to drought stress conditions and their differential response was studied by analysing different morpho-physiological traits. Utilizing the in-house, high throughput sequencing data, the SQUAMOSA promoter-binding (SBP) protein-like (SPL) transcription factor genes were identified to be differentially methylated and expressed amongst the two genotypes, in response to drought stress. The methylation status of one of these genes was examined and validated through bisulfite PCR (BS-PCR). The identified genes could be possible homologs to known epialleles and can therefore serve as potential epialleles which can be utilized for crop improvement in chickpea. CONCLUSION: The SPL TF genes are potential targets of epigenetic regulation in response to drought stress in chickpea. Since these are TFs, they might play important roles in controlling the expression of other genes, thus contributing to differential drought response of the two genotypes.

Genome-wide comprehensive characterization and expression analysis of TLP gene family revealed its responses to hormonal and abiotic stresses in watermelon (Citrullus lanatus).

Thaumatin-like protein (TLP) is the well-known sweetest protein which plays a crucial role in diverse developmental processes and different stress conditions in plants, fungi and animals. The TLP gene family is extensively studied in different plant species including crop plants. Watermelon (Citrullus lanatus) is an important cucurbit crop cultivated worldwide; however, the comprehensive information about the TLP gene family is not available in watermelon. In the present study, we identified the 29 TLP genes as gene family members in watermelon using various computational methods to understand its role in different developmental processes and stress conditions. ClaTLP gene family members were not uniformly distributed on 22 chromosomes. Phylogenetic analysis revealed that the ClaTLP gene family members were grouped into 10 sub-groups. Further, gene duplication analysis showed thirteen gene duplication events which included one tandem and twelve segmental duplications. Amino acid sequence alignment has shown that ClaTLP proteins shared 16 conserved cysteine residues in their THN domain. Furthermore, cis-acting regulatory elements analysis also displayed that ClaTLP gene family members contain diverse phytohormone, various defense, and stress-responsive elements in their promoter region. The expression profile of the ClaTLP gene family revealed the differential expression of gene family members in different tissues and abiotic stresses conditions. Moreover, the expression profile of ClaTLP genes was further validated by semi-quantitative reverse transcriptase PCR. Taken together, these results indicate that ClaTLP genes might play an important role in developmental processes and diverse stress conditions. Therefore, the outcome of this study brings forth the valuable information for further interpret the precise role of ClaTLP gene family members in watermelon.

Genome-wide identification, in silico characterization and expression analysis of the RNA helicase gene family in chickpea (C. arietinum L.).

The RNA helicases are an important class of enzymes which are known to influence almost every aspect of RNA metabolism. The majority of RNA helicases belong to the SF2 (superfamily 2) superfamily, members of which are further categorized into three separate subfamilies i.e., the DEAD, DEAH and DExD/H-box subfamilies. In chickpea, these RNA helicases have not been characterized until now. A genome-wide analysis across the chickpea genome led to the identification of a total of 150 RNA helicase genes which included 50 DEAD, 33 DEAH and 67 DExD/H-box genes. These were distributed across all the eight chromosomes, with highest number on chromosome 4 (26) and least on chromosome 8 (8). Gene duplication analysis resulted in identification of 15 paralogous gene pairs with Ka/Ks values < 1, indicating towards the genes being under purifying selection during the course of evolution. The promoter regions of the RNA helicase genes were enriched in cis-acting elements like the light and ABA-responsive elements. The drought responsiveness of the genes was analysed by studying the expression profiles of few of these genes, in two different genotypes, the cultivated variety ICC 8261 (kabuli, C. arietinum) and the wild accession ILWC 292 (C. reticulatum), through qRT-PCR. These genotypes were selected based on their drought responsiveness in a field experiment, where it was observed that the percentage (%) reduction in relative water content (RWC) and membrane stability index (MSI) for the drought stressed plants after withholding water for 24 days, over the control or well-watered plants, was least for both the genotypes. The genes CaDEAD50 and CaDExD/H66 were identified as drought-responsive RNA helicase genes in chickpea. The protein encoded by the CaDExD/H66 gene shares a high degree of homology with one of the CLSY (CLASSY) proteins of A. thaliana. We hypothesize that this gene could possibly be involved in regulation of DNA methylation levels in chickpea by regulating siRNA production, in conjunction with other proteins like the Argonaute, RNA dependent RNA polymerases and Dicer-like proteins.

Genetically modified foods: bibliometric analysis on consumer perception and preference.

In this study, we present the bibliometric trends emerging from research outputs on consumer perception and preference for genetically modified (GM) foods and policy prescriptions for enabling the consumption using VOSviewer visualization software. Consumers' positive response is largely influenced by the decision of the governments to ban or approve the GM crops cultivation. Similarly, the public support increases when the potential benefits of the technology are well articulated, consumption increases with a price discount, people's trust on the government and belief in science increases with a positive influence by the media. Europe and the USA are the first region and country, respectively, in terms of the number of active institutions per research output, per-capita GDP publication and citations. We suggest research-, agri-food industries-, and society-oriented policies to be implemented by the stakeholders to ensure the safety of GM foods, encourage consumer-based studies, and increase public awareness toward these food products.

Genome-wide identification and expression analysis of the GRAS gene family in response to drought stress in chickpea (Cicer arietinum L.).

The GRAS (gibberellic acid insensitive, repressor of GAI and scarecrow) transcription factors (TFs) regulate diverse biological processes involved in plant growth and development. These TFs are also known to regulate gene expression in response to various abiotic stress factors like cold, drought, etc. In chickpea one of the most devastating abiotic stress factors is terminal drought. The GRAS TF family has not been characterized in chickpea (Cicer arietinum L.) until now. In this study, we report 46 GRAS TF genes (CaGRAS genes) in the chickpea genome. The CaGRAS proteins were categorized into nine subfamilies based on their phylogenetic relationship with known GRAS members of Arabidopsis and soybean. The PAT subfamily was the largest consisting of ten CaGRAS members whereas the LAS subfamily was the smallest with only one member. Gene duplication analysis revealed that segmental duplication was the primary reason for the expansion of this gene family within the chickpea genome. The gene expression levels of CaGRAS genes were analysed using two different chickpea varieties contrasting for drought tolerance trait, i.e., ICC 4958 (drought tolerant) and ICC 1882 (drought sensitive). On exposure to drought stress, the two chickpea genotypes, exhibited differential drought response, which was quantified and estimated in terms of differences in leaf relative water content (RWC). The well-watered or control plants of the drought tolerant variety were able to maintain a higher leaf RWC by the end of the drought stress period, whereas the control plants of the drought sensitive variety continued to show a decline in leaf RWC. The two genotypes also differed in their root morphologies, under well-watered and drought stress conditions. The gene expression analysis revealed a potential role of PAT, SCR, SCL3 and SHR GRAS members in the regulation of differential response to drought, in the root tissues, for both the genotypes. CaGRAS 12 (SCR) was identified as a drought-responsive GRAS TF gene, which could serve as a potential candidate gene for utilization in developing chickpea varieties with improved drought tolerance. This study demonstrates the drought-responsive expression of CaGRAS genes in chickpea and also describes the morpho-physiological response of chickpea plants to drought stress conditions. SUPPLEMENTARY INFORMATION: The online version contains supplementary material available at 10.1007/s13205-021-03104-z.

Assessment of genetic diversity in Indian rice germplasm (Oryza sativa L.): use of random versus trait-linked microsatellite markers.

Assessment of genetic diversity in a crop germplasm is a vital part of plant breeding. DNA markers such as microsatellite or simple sequence repeat markers have been widely used to estimate the genetic diversity in rice. The present study was carried out to decipher the pattern of genetic diversity in terms of both phenotypic and genotypic variability, and to assess the efficiency of random vis-á-vis QTL linked/gene based simple sequence repeat markers in diversity estimation. A set of 88 rice accessions that included landraces, farmer's varieties and popular Basmati lines were evaluated for agronomic traits and molecular diversity. The random set of SSR markers included 50 diversity panel markers developed under IRRI's Generation Challenge Programme (GCP) and the trait-linked/gene based markers comprised of 50 SSR markers reportedly linked to yield and related components. For agronomic traits, significant variability was observed, ranging between the maximum for grains/panicle and the minimum for panicle length. The molecular diversity based grouping indicated that varieties from a common centre were genetically similar, with few exceptions. The trait-linked markers gave an average genetic dissimilarity of 0.45 as against that of 0.37 by random markers, along with an average polymorphic information constant value of 0.48 and 0.41 respectively. The correlation between the kinship matrix generated by trait-linked markers and the phenotype based distance matrix (0.29) was higher than that of random markers (0.19). This establishes the robustness of trait-linked markers over random markers in estimating genetic diversity of rice germplasm.