Genome-wide mapping of DNase I hypersensitive sites revealed differential chromatin accessibility and regulatory DNA elements under drought stress in rice cultivars.

Drought stress (DS) is one of the major constraints limiting yield in crop plants including rice. Gene regulation under DS is largely governed by accessibility of the transcription factors (TFs) to their cognate cis-regulatory elements (CREs). In this study, we used DNase I hypersensitive assays followed by sequencing to identify the accessible chromatin regions under DS in a drought-sensitive (IR64) and a drought-tolerant (N22) rice cultivar. Our results indicated that DNase I hypersensitive sites (DHSs) were highly enriched at transcription start sites (TSSs) and numerous DHSs were detected in the promoter regions. DHSs were concurrent with epigenetic marks and the genes harboring DHSs in their TSS and promoter regions were highly expressed. In addition, DS induced changes in DHSs (∆DHSs) in TSS and promoter regions were positively correlated with upregulation of several genes involved in drought/abiotic stress response, those encoding TFs and located within drought-associated quantitative trait loci, much preferentially in the drought-tolerant cultivar. The CREs representing the binding sites of TFs involved in DS response were detected within the ∆DHSs, suggesting differential accessibility of TFs to their cognate sites under DS in different rice cultivars, which may be further deployed for enhancing drought tolerance in rice.

Draft genome sequence of Indian mulberry (Morus indica) provides a resource for functional and translational genomics.

Mulberry is an important crop plant for the sericulture industry. Here, we report high-quality genome sequence of a cultivated Indian mulberry (Morus indica cv K2) obtained by combining data from four different technologies, including Illumina, single-molecule real-time sequencing, chromosome conformation capture and optical mapping, with a gene completeness of 96.5%. Based on the genome sequence, we identified 49.2% of repetitive DNA and 27,435 high-confidence protein-coding genes with >90% of them supported by transcript evidence. A comparative analysis with other plant genomes identified 4.8% of species-specific genes in the M. indica genome. Transcriptome profiling revealed tissue-specific and differential expression across multiple accessions of ~4.7% and 2-5% of protein-coding genes, respectively, implicated in diverse biological processes. Whole genome resequencing of 21 accessions/species revealed ~2.5 million single nucleotide polymorphisms and ~ 0.2 million insertions/deletions. These data and results provide a comprehensive resource to accelerate the genomics research in mulberry for its improvement.

Genome-wide analysis suggests the potential role of lncRNAs during seed development and seed size/weight determination in chickpea.

MAIN CONCLUSION: The integrated transcriptome data analyses suggested the plausible roles of lncRNAs during seed development in chickpea. The candidate lncRNAs associated with QTLs and those involved in miRNA-mediated seed size/weight determination in chickpea have been identified. Long non-coding RNAs (lncRNAs) are important regulators of various biological processes. Here, we identified lncRNAs at seven successive stages of seed development in small-seeded and large-seeded chickpea cultivars. In total, 4751 lncRNAs implicated in diverse biological processes were identified. Most of lncRNAs were conserved between the two cultivars, whereas only a few of them were conserved in other plants, suggesting their species-specificity. A large number of lncRNAs differentially expressed between the two chickpea cultivars associated with seed development-related processes were identified. The lncRNAs acting as precursors of miRNAs and those mimicking target protein-coding genes of miRNAs involved in seed size/weight determination, including HAIKU1, BIG SEEDS1, and SHB1, were also revealed. Further, lncRNAs located within seed size/weight associated quantitative trait loci were also detected. Overall, we present a comprehensive resource and identified candidate lncRNAs that may play important roles during seed development and seed size/weight determination in chickpea.

Genome-wide discovery of DNA polymorphisms via resequencing of chickpea cultivars with contrasting response to drought stress.

Drought stress limits plant growth, resulting in a significant yield loss in chickpea. The diversification in genome sequence and selective sweep of allele(s) in different genotypes of a crop plant may play an important role in the determination of agronomic traits, including drought stress response. We investigated, via whole genome resequencing, the DNA polymorphisms between two sets of chickpea genotypes with contrasting drought stress responses (3 drought-sensitive vs. 6 drought-tolerant). In total, 36,406 single nucleotide polymorphisms (SNPs) and 3407 insertions or deletions (InDels) differentiating drought-sensitive and drought-tolerant chickpea genotypes were identified. Interestingly, most (91%) of these DNA polymorphisms were located in chromosomes 1 and 4. The genes harboring DNA polymorphisms in their promoter and/or coding regions and exhibiting differential expression under control and/or drought stress conditions between/within the drought-sensitive and tolerant genotypes were found implicated in the stress response. Furthermore, we identified DNA polymorphisms within the cis-regulatory motifs in the promoter region of abiotic stress-related and QTL-associated genes, which might contribute to the differential expression of the candidate drought-responsive genes. In addition, we revealed the effect of nonsynonymous SNPs on mutational sensitivity and stability of the encoded proteins. Taken together, we identified DNA polymorphisms having relevance in drought stress response and revealed candidate genes to engineer drought tolerance in chickpea.

Genome resequencing reveals DNA polymorphisms associated with seed size/weight determination in chickpea.

Crop productivity in legumes is determined by number and size/weight of seeds. To understand the genetic basis of seed size/weight in chickpea, we performed genome resequencing of 13 small- and 5 large-seeded genotypes using Illumina platform. Single nucleotide polymorphisms (SNPs) and insertions/deletions (InDels) differentiating small- and large-seeded genotypes were identified. A total of 17,902 SNPs and 2594 InDels located in promoter and/or coding regions that may contribute to seed size/weight were detected. Of these, 266 SNPs showed significant association with seed size/weight trait. Twenty-three genes including those involved in cell growth/division, encoding transcription factors and located within QTLs associated with seed size/weight harbored SNPs within transcription factor binding motif(s) and/or coding region. The non-synonymous SNPs were found to affect the mutational sensitivity and stability of the encoded proteins. Overall, we provided a high-quality SNP map for large-scale genotyping applications and identified candidate genes that determine seed size/weight in chickpea.

Discovery of DNA polymorphisms via whole genome resequencing and their functional relevance in salinity stress response in chickpea.

Salinity stress is one of the major constraints for plant growth and yield. The salinity stress response of different genotypes of crop plants may largely be governed by DNA polymorphisms. To determine the molecular genetic factors involved in salinity stress tolerance in chickpea, we performed a whole genome resequencing data analysis of three each of salinity-sensitive and salinity-tolerant genotypes. A total of 6173 single nucleotide polymorphisms and 920 insertions and deletions differentiating the chickpea genotypes with contrasting salinity stress responses were identified. Gene ontology analysis revealed the enrichment of functional terms related to stress response and development among the genes harboring DNA polymorphisms in their promoter and/or coding regions. DNA polymorphisms located within the cis-regulatory motifs of the quantitative trait loci (QTL)-associated and abiotic stress related genes were identified, which may influence salinity stress response via modulating binding affinity of the transcription factors. Several genes including QTL-associated and abiotic stress response related genes harboring DNA polymorphisms exhibited differential expression in response to salinity stress especially at the reproductive stage of development in the salinity-tolerant genotype. Furthermore, effects of non-synonymous DNA polymorphisms on mutational sensitivity and structural integrity of the encoded proteins by the candidate QTL-associated and abiotic stress response related genes were revealed. The results suggest that DNA polymorphisms may determine salinity stress response via influencing differential gene expression in genotype and/or stage-dependent manner. Altogether, we provide a high-quality set of DNA polymorphisms and candidate genes that may govern salinity stress tolerance in chickpea.

Draft genome and transcriptome analyses of halophyte rice Oryza coarctata provide resources for salinity and submergence stress response factors.

Oryza coarctata is a wild relative of rice that has adapted to diverse ecological environments, including high salinity and submergence. Thus, it can provide an important resource for discovering candidate genes/factors involved in tolerance to these stresses. Here, we report a draft genome assembly of 573 Mb comprised of 8877 scaffolds with N50 length of 205 kb. We predicted a total of 50,562 protein-coding genes, of which a significant fraction was found to be involved in secondary metabolite biosynthesis and hormone signal transduction pathways. Several salinity and submergence stress-responsive protein-coding and long noncoding RNAs involved in diverse biological processes were identified using RNA-sequencing data. Based on small RNA sequencing, we identified 168 unique miRNAs and 3219 target transcripts (coding and noncoding) involved in several biological processes, including abiotic stress responses. Further, whole genome bisulphite sequencing data analysis revealed at least 19%-48% methylcytosines in different sequence contexts and the influence of methylation status on gene expression. The genome assembly along with other datasets have been made publicly available at http://ccbb.jnu.ac.in/ory-coar. Altogether, we provide a comprehensive genomic resource for understanding the regulation of salinity and submergence stress responses and identification of candidate genes/factors involved for functional genomics studies.

Bisulphite sequencing reveals dynamic DNA methylation under desiccation and salinity stresses in rice cultivars.

DNA methylation governs gene regulation in plants in response to environmental conditions. Here, we analyzed role of DNA methylation under desiccation and salinity stresses in three (IR64, stress-sensitive; Nagina 22, drought-tolerant and Pokkali, salinity-tolerant) rice cultivars via bisulphite sequencing. Methylation in CG context within gene body and methylation in CHH context in distal promoter regions were positively correlated with gene expression. Hypomethylation in Nagina 22 and hypermethylation in Pokkali in response to desiccation and salinity stresses, respectively, were correlated with higher expression of few abiotic stress response related genes. Most of the differentially methylated and differentially expressed genes (DMR-DEGs) were cultivar-specific, suggesting an important role of DNA methylation in abiotic stress responses in rice in cultivar-specific manner. DMR-DEGs harboring differentially methylated cytosines due to DNA polymorphisms between the sensitive and tolerant cultivars in their promoter regions and/or coding regions were identified, suggesting the role of epialleles in abiotic stress responses.

Genome-wide analysis of glutathione S-transferase gene family in chickpea suggests its role during seed development and abiotic stress.

Glutathione S-transferases (GSTs) are multifunctional proteins that help in oxidative stress metabolism and detoxification of xenobiotic compounds. Studies pertaining to GST gene family have been undertaken in various plant species, however no information is available with respect to GST genes in chickpea. In the current study, we identified a total of 51 GST encoding genes in chickpea (CaGST) genome. Phylogenetic analysis revealed that GST gene family can be divided into eleven distinct classes. Tau and phi were the major classes in chickpea and one third of the CaGST genes represented segmental duplication and purifying selection was common among these genes. Expression of many CaGST genes, in particular, members of tau class were found to be upregulated under abiotic stress conditions. In addition, CaGST genes displayed differential expression patterns across diverse organs/tissues, suggesting their roles in developmental processes. Many CaGST genes showed opposite expression pattern in small- and large-seeded chickpea cultivars during seed development. Higher expression of CaGST genes in small-seeded cultivar at maturation stages of seed development suggested their important role in seed development and seed size/weight determination in chickpea. Overall, these results provide a comprehensive information on GST gene family members in chickpea and is expected to provide a rational platform to explore versatile role of these genes in semi-arid legume crops.

PLncPRO for prediction of long non-coding RNAs (lncRNAs) in plants and its application for discovery of abiotic stress-responsive lncRNAs in rice and chickpea.

Long non-coding RNAs (lncRNAs) make up a significant portion of non-coding RNAs and are involved in a variety of biological processes. Accurate identification/annotation of lncRNAs is the primary step for gaining deeper insights into their functions. In this study, we report a novel tool, PLncPRO, for prediction of lncRNAs in plants using transcriptome data. PLncPRO is based on machine learning and uses random forest algorithm to classify coding and long non-coding transcripts. PLncPRO has better prediction accuracy as compared to other existing tools and is particularly well-suited for plants. We developed consensus models for dicots and monocots to facilitate prediction of lncRNAs in non-model/orphan plants. The performance of PLncPRO was quite better with vertebrate transcriptome data as well. Using PLncPRO, we discovered 3714 and 3457 high-confidence lncRNAs in rice and chickpea, respectively, under drought or salinity stress conditions. We investigated different characteristics and differential expression under drought/salinity stress conditions, and validated lncRNAs via RT-qPCR. Overall, we developed a new tool for the prediction of lncRNAs in plants and showed its utility via identification of lncRNAs in rice and chickpea.

Global transcriptome and coexpression network analyses reveal cultivar-specific molecular signatures associated with seed development and seed size/weight determination in chickpea.

Seed development is an intricate process regulated via a complex transcriptional regulatory network. To understand the molecular mechanisms governing seed development and seed size/weight in chickpea, we performed a comprehensive analysis of transcriptome dynamics during seed development in two cultivars with contrasting seed size/weight (small-seeded, Himchana 1 and large-seeded, JGK 3). Our analysis identified stage-specific expression for a significant proportion (>13%) of the genes in each cultivar. About one half of the total genes exhibited significant differential expression in JGK 3 as compared with Himchana 1. We found that different seed development stages can be delineated by modules of coexpressed genes. A comparative analysis revealed differential developmental stage specificity of some modules between the two cultivars. Furthermore, we constructed transcriptional regulatory networks and identified key components determining seed size/weight. The results suggested that extended period of cell division during embryogenesis and higher level of endoreduplication along with more accumulation of storage compounds during maturation determine large seed size/weight. Further, we identified quantitative trait loci-associated candidate genes harboring single nucleotide polymorphisms in the promoter sequences that differentiate small- and large-seeded chickpea cultivars. The results provide a valuable resource to dissect the role of candidate genes governing seed development and seed size/weight in chickpea.

Long intergenic non-coding RNAs modulate proximal protein-coding gene expression and tolerance to Candidatus Liberibacter spp. in potatoes.

Long intergenic non-coding RNAs (lincRNAs) are emerging as regulators of protein-coding genes (PCGs) in many plant and animal developmental processes and stress responses. In this study, we characterize the genome-wide lincRNAs in potatoes responsive to a vascular bacterial disease presumably caused by Candidatus Liberibacter solanacearum (CLso). Approximately 4397 lincRNAs were detected in healthy and infected potato plants at various stages of zebra chip (ZC) disease progression. Of them, ~65% (2844) were novel lincRNAs, and less than 1% (9) were orthologs of Arabidopsis and rice based on reciprocal BLAST analysis, suggesting species-specific expansion. Among the proximal lincRNAs within 50 kbp from a PCG, ~49% were transcribed from the same strand, while ~39% and ~15% followed convergent (head-to-head) and divergent (tail-to-tail) orientations, respectively. Approximately 30% (1308) were differentially expressed following CLso infection, with substantial changes occurring 21 days after infection (DAI). Weighted Gene Co-expression Network Analysis (WGCNA) of lincRNAs and PCGs identified 46 highly correlated lincRNA-PCG pairs exhibiting co-up or co-downregulation. Furthermore, overexpression of selected lincRNAs in transgenic potato hairy roots resulted in perturbation of neighboring PCG expression and conferred tolerance to CLso infection. Our results provide novel insights into potato lincRNAs' identity, expression dynamics, and functional relevance to CLso infection.

DNA Methylation-Dependent Restriction of Tyrosine Hydroxylase Contributes to Pancreatic ß-Cell Heterogeneity.

The molecular and functional heterogeneity of pancreatic ß-cells is well recognized, but the underlying mechanisms remain unclear. Pancreatic islets harbor a subset of ß-cells that co-express tyrosine hydroxylase (TH), an enzyme involved in synthesis of catecholamines that repress insulin secretion. Restriction of the TH+ ß-cells within islets is essential for appropriate function in mice, such that a higher proportion of these cells corresponds to reduced insulin secretion. Here, we use these cells as a model to dissect the developmental control of ß-cell heterogeneity. We define the specific molecular and metabolic characteristics of TH+ ß-cells and show differences in their developmental restriction in mice and humans. We show that TH expression in ß-cells is restricted by DNA methylation during ß-cell differentiation. Ablation of de novo DNA methyltransferase Dnmt3a in the embryonic progenitors results in a dramatic increase in the proportion of TH+ ß-cells, whereas ß-cell-specific ablation of Dnmt3a does not. We demonstrate that maintenance of Th promoter methylation is essential for its continued restriction in postnatal ß-cells. Loss of Th promoter methylation in response to chronic overnutrition increases the number of TH+ ß-cells, corresponding to impaired ß-cell function. These results reveal a regulatory role of DNA methylation in determining ß-cell heterogeneity.

An integrated transcriptome mapping the regulatory network of coding and long non-coding RNAs provides a genomics resource in chickpea.

Large-scale transcriptome analysis can provide a systems-level understanding of biological processes. To accelerate functional genomic studies in chickpea, we perform a comprehensive transcriptome analysis to generate full-length transcriptome and expression atlas of protein-coding genes (PCGs) and long non-coding RNAs (lncRNAs) from 32 different tissues/organs via deep sequencing. The high-depth RNA-seq dataset reveal expression dynamics and tissue-specificity along with associated biological functions of PCGs and lncRNAs during development. The coexpression network analysis reveal modules associated with a particular tissue or a set of related tissues. The components of transcriptional regulatory networks (TRNs), including transcription factors, their cognate cis-regulatory motifs, and target PCGs/lncRNAs that determine developmental programs of different tissues/organs, are identified. Several candidate tissue-specific and abiotic stress-responsive transcripts associated with quantitative trait loci that determine important agronomic traits are also identified. These results provide an important resource to advance functional/translational genomic and genetic studies during chickpea development and environmental conditions.

DNA methylation reprogramming during seed development and its functional relevance in seed size/weight determination in chickpea.

Seed development is orchestrated via complex gene regulatory networks and pathways. Epigenetic factors may also govern seed development and seed size/weight. Here, we analyzed DNA methylation in a large-seeded chickpea cultivar (JGK 3) during seed development stages. Progressive gain of CHH context DNA methylation in transposable elements (TEs) and higher frequency of small RNAs in hypermethylated TEs during seed development suggested a role of the RNA-dependent DNA methylation pathway. Frequency of intragenic TEs was higher in CHH context differentially methylated region (DMR) associated differentially expressed genes (DEGs). CG context hyper/hypomethylation within the gene body was observed for most of DMR-associated DEGs in JGK 3 as compared to small-seeded chickpea cultivar (Himchana 1). We identified candidate genes involved in seed size/weight determination exhibiting CG context hypermethylation within the gene body and higher expression in JGK 3. This study provides insights into the role of DNA methylation in seed development and seed size/weight determination in chickpea.

Updates on Genomic Resources in Chickpea for Crop Improvement.

In recent years, rapid advancement has been done in generation of genomic resources for the important legume crop chickpea. Here, we provide an update on important advancements made on availability of genomic resources for this crop. The availability of reference genome and transcriptome sequences, and resequencing of several accessions have enabled the discovery of gene space and molecular markers in chickpea. These resources have helped in elucidating evolutionary relationship and identification of quantitative trait loci for important agronomic traits. Gene expression in different tissues/organs during development and under abiotic/biotic stresses has been interrogated. In addition, single-base resolution DNA methylation patterns in different organs have been analyzed to understand gene regulation. Overall, we provide a consolidated overview of available genomic resources of chickpea that may help in fulfilling the promises for improvement of this important crop.

Genome-wide discovery of DNA polymorphisms among chickpea cultivars with contrasting seed size/weight and their functional relevance.

Seed size/weight is a major agronomic trait which determine crop productivity in legumes. To understand the genetic basis of seed size determination, we sought to identify DNA polymorphisms between two small (Himchana 1 and Pusa 362) and two large-seeded (JGK 3 and PG 0515) chickpea cultivars via whole genome resequencing. We identified a total of 75535 single nucleotide polymorphisms (SNPs), 6486 insertions and deletions (InDels), 1938 multi-nucleotide polymorphisms (MNPs) and 5025 complex variants between the two small and two large-seeded chickpea cultivars. Our analysis revealed 814, 244 and 72 seed-specific genes harboring DNA polymorphisms in promoter or non-synonymous and large-effect DNA polymorphisms, respectively. Gene ontology analysis revealed enrichment of cell growth and division related terms in these genes. Among them, at least 22 genes associated with quantitative trait loci, and those involved in cell growth and division and encoding transcription factors harbored promoter and/or large-effect/non-synonymous DNA polymorphisms. These also showed higher expression at late-embryogenesis and/or mid-maturation stages of seed development in the large-seeded cultivar, suggesting their role in seed size/weight determination in chickpea. Altogether, this study provided a valuable resource for large-scale genotyping applications and a few putative candidate genes that might play crucial role in governing seed size/weight in chickpea.

Genome-wide identification and co-expression network analysis provide insights into the roles of auxin response factor gene family in chickpea.

Auxin response factors (ARFs) are the transcription factors that regulate auxin responses in various aspects of plant growth and development. Although genome-wide analysis of ARF gene family has been done in some species, no information is available regarding ARF genes in chickpea. In this study, we identified 28 ARF genes (CaARF) in the chickpea genome. Phylogenetic analysis revealed that CaARFs can be divided into four different groups. Duplication analysis revealed that 50% of CaARF genes arose from duplication events. We analyzed expression pattern of CaARFs in various developmental stages. CaARF16.3, CaARF17.1 and CaARF17.2 showed highest expression at initial stages of flower bud development, while CaARF6.2 had higher expression at later stages of flower development. Further, CaARF4.2, CaARF9.2, CaARF16.2 and CaARF7.1 exhibited differential expression under different abiotic stress conditions, suggesting their role in abiotic stress responses. Co-expression network analysis among CaARF, CaIAA and CaGH3 genes enabled us to recognize components involved in the regulatory network associated with CaARFs. Further, we identified microRNAs that target CaARFs and TAS3 locus that trigger production of trans-acting siRNAs targeting CaARFs. The analyses presented here provide comprehensive information on ARF family members and will help in elucidating their exact function in chickpea.