Gene discovery and tissue-specific transcriptome analysis in chickpea with massively parallel pyrosequencing and web resource development

Drought priming modulates ABF, GRFs, related microRNAs and induce metabolic adjustment during heat stress in chickpea

Drought and high temperature stress may occur concomitantly or individually in succession causing cellular dysfunctions. Abscisic acid (ABA) is a key stress regulator, and its responsive genes are controlled by ABRE (Abscisic acid Responsive Element)-binding factors (ABFs)and G-Box Regulatory factors (GRFs). Here, we identify ABFs, GRFs and targeting miRNAs in desi and kabuli chickpea. To validate their role after drought priming and subsequent high temperature stress, two contrasting chickpea varieties (PBG1 and PBG5) were primed and exposed to 32 °C, 35 °C and 38 °C for 12, 6 and 2 h respectively and analyzed for Physio-biochemical, expression of ABFs, GRFs and MiRNAs, and GC-MS based metabolite analysis. To ascertain the ABF-GRF protein-protein interactions, docking studies were carried out between the ABF3 and GRF14. Genome-wide analysis identified total 9 & 11 ABFs, and 11 GRFsin desi and kabuli respectively. Their gene structure, and motif composition were conserved in all subfamilies and only 10 and 12 genes have undergone duplication in both desi and kabuli chickpea respectively. These genes were differentially expressed in-silico. MiR172 and miR396 were identified to target ABFs and GRFs respectively. Protein-protein interaction (ABF3 and GRF14) might be successful only when the ABF3 was phosphorylated. Drought priming downregulated miR172 and miR396 and eventually upregulated targeting ABFs, and GRFs. Metabolite profiling (GC-MS) revealed the accumulation of 87 metabolites in Primed (P) and Non-Primed (NP) Chickpea plants. Tolerant cultivar (PBG5) responded better in all respects however both severity of stress and exposure are important factors and can produce broadly similar cellular response.

Elucidating the Mesocarp Drupe Transcriptome of Açai (Euterpe oleracea Mart.): An Amazonian Tree Palm Producer of Bioactive Compounds

Euterpe oleracea palm, endemic to the Amazon region, is well known for açai, a fruit violet beverage with nutritional and medicinal properties. During E. oleracea fruit ripening, anthocyanin accumulation is not related to sugar production, contrarily to grape and blueberry. Ripened fruits have a high content of anthocyanins, isoprenoids, fibers, and proteins, and are poor in sugars. E. oleracea is proposed as a new genetic model for metabolism partitioning in the fruit. Approximately 255 million single-end-oriented reads were generated on an Ion Proton NGS platform combining fruit cDNA libraries at four ripening stages. The de novo transcriptome assembly was tested using six assemblers and 46 different combinations of parameters, a pre-processing and a post-processing step. The multiple k-mer approach with TransABySS as an assembler and Evidential Gene as a post-processer have shown the best results, with an N50 of 959 bp, a read coverage mean of 70x, a BUSCO complete sequence recovery of 36% and an RBMT of 61%. The fruit transcriptome dataset included 22,486 transcripts representing 18 Mbp, of which a proportion of 87% had significant homology with other plant sequences. Approximately 904 new EST-SSRs were described, and were common and transferable to Phoenix dactylifera and Elaeis guineensis, two other palm trees. The global GO classification of transcripts showed similar categories to that in P. dactylifera and E. guineensis fruit transcriptomes. For an accurate annotation and functional description of metabolism genes, a bioinformatic pipeline was developed to precisely identify orthologs, such as one-to-one orthologs between species, and to infer multigenic family evolution. The phylogenetic inference confirmed an occurrence of duplication events in the Arecaceae lineage and the presence of orphan genes in E. oleracea. Anthocyanin and tocopherol pathways were annotated entirely. Interestingly, the anthocyanin pathway showed a high number of paralogs, similar to in grape, whereas the tocopherol pathway exhibited a low and conserved gene number and the prediction of several splicing forms. The release of this exhaustively annotated molecular dataset of E. oleracea constitutes a valuable tool for further studies in metabolism partitioning and opens new great perspectives to study fruit physiology with açai as a model.

De Novo Transcriptome Profiling for the Generation and Validation of Microsatellite Markers, Transcription Factors, and Database Development for Andrographis paniculata

Andrographis paniculata belongs to the family Acanthaceae and is known for its medicinal properties owing to the presence of unique constituents belonging to the lactones, diterpenoids, diterpene glycosides, flavonoids, and flavonoid glycosides groups of chemicals. Andrographolide, a major therapeutic constituent of A. paniculata, is extracted primarily from the leaves of this plant and exhibits antimicrobial and anti-inflammatory activities. Using 454 GS-FLX pyrosequencing, we have generated a whole transcriptome profile of entire leaves of A. paniculata. A total of 22,402 high-quality transcripts were generated, with an average transcript length and N50 of 884 bp and 1007 bp, respectively. Functional annotation revealed that 19,264 (86%) of the total transcripts showed significant similarity with the NCBI-Nr database and were successfully annotated. Out of the 19,264 BLAST hits, 17,623 transcripts were assigned GO terms and distributed into three major functional categories: molecular function (44.62%), biological processes (29.19%), and cellular component (26.18%) based on BLAST2GO. Transcription factor analysis showed 6669 transcripts, belonging to 57 different transcription factor families. Fifteen TF genes that belong to the NAC, MYB, and bHLH TF categories were validated by RT PCR amplification. In silico analysis of gene families involved in the synthesis of biochemical compounds having medicinal values, such as cytochrome p450, protein kinases, heat shock proteins, and transporters, was completed and a total of 102 different transcripts encoding enzymes involved in the biosynthesis of terpenoids were predicted. Out of these, 33 transcripts belonged to terpenoid backbone biosynthesis. This study also identified 4254 EST-SSRs from 3661 transcripts, representing 16.34% of the total transcripts. Fifty-three novel EST-SSR markers generated from our EST dataset were used to assess the genetic diversity among eighteen A. paniculata accessions. The genetic diversity analysis revealed two distinct sub-clusters and all accessions based on the genetic similarity index were distinct from each other. A database based on EST transcripts, EST-SSR markers, and transcription factors has been developed using data generated from the present study combined with available transcriptomic resources from a public database using Meta transcriptome analysis to make genomic resources available in one place to the researchers working on this medicinal plant.

Underutilized legumes: nutrient status and advanced breeding approaches for qualitative and quantitative enhancement

Underutilized/orphan legumes provide food and nutritional security to resource-poor rural populations during periods of drought and extreme hunger, thus, saving millions of lives. The Leguminaceae, which is the third largest flowering plant family, has approximately 650 genera and 20,000 species and are distributed globally. There are various protein-rich accessible and edible legumes, such as soybean, cowpea, and others; nevertheless, their consumption rate is far higher than production, owing to ever-increasing demand. The growing global urge to switch from an animal-based protein diet to a vegetarian-based protein diet has also accelerated their demand. In this context, underutilized legumes offer significant potential for food security, nutritional requirements, and agricultural development. Many of the known legumes like Mucuna spp., Canavalia spp., Sesbania spp., Phaseolus spp., and others are reported to contain comparable amounts of protein, essential amino acids, polyunsaturated fatty acids (PUFAs), dietary fiber, essential minerals and vitamins along with other bioactive compounds. Keeping this in mind, the current review focuses on the potential of discovering underutilized legumes as a source of food, feed and pharmaceutically valuable chemicals, in order to provide baseline data for addressing malnutrition-related problems and sustaining pulse needs across the globe. There is a scarcity of information about underutilized legumes and is restricted to specific geographical zones with local or traditional significance. Around 700 genera and 20,000 species remain for domestication, improvement, and mainstreaming. Significant efforts in research, breeding, and development are required to transform existing local landraces of carefully selected, promising crops into types with broad adaptability and economic viability. Different breeding efforts and the use of biotechnological methods such as micro-propagation, molecular markers research and genetic transformation for the development of underutilized crops are offered to popularize lesser-known legume crops and help farmers diversify their agricultural systems and boost their profitability.

Breeding and Genomic Approaches towards Development of Fusarium Wilt Resistance in Chickpea

Chickpea is an important leguminous crop with potential to provide dietary proteins to both humans and animals. It also ameliorates soil nitrogen through biological nitrogen fixation. The crop is affected by an array of biotic and abiotic factors. Among different biotic stresses, a major fungal disease called Fusarium wilt, caused by Fusarium oxysporum f. sp. ciceris (FOC), is responsible for low productivity in chickpea. To date, eight pathogenic races of FOC (race 0, 1A, and 1B/C, 2-6) have been reported worldwide. The development of resistant cultivars using different conventional breeding methods is very time consuming and depends upon the environment. Modern technologies can improve conventional methods to solve these major constraints. Understanding the molecular response of chickpea to Fusarium wilt can help to provide effective management strategies. The identification of molecular markers closely linked to genes/QTLs has provided great potential for chickpea improvement programs. Moreover, omics approaches, including transcriptomics, metabolomics, and proteomics give scientists a vast viewpoint of functional genomics. In this review, we will discuss the integration of all available strategies and provide comprehensive knowledge about chickpea plant defense against Fusarium wilt.

Transcriptomic Analysis of Flowering Time Genes in Cultivated Chickpea and Wild Cicer

Chickpea (Cicer arietinum L.) is a major grain legume and a good source of plant-based protein. However, comprehensive knowledge of flowering time control in Cicer is lacking. In this study, we acquire high-throughput transcriptome sequencing data and analyze changes in gene expression during floral transition in the early flowering cultivar ICCV 96029, later flowering C. arietinum accessions, and two wild species, C. reticulatum and C. echinospermum. We identify Cicer orthologs of A. thaliana flowering time genes and analyze differential expression of 278 genes between four species/accessions, three tissue types, and two conditions. Our results show that the differences in gene expression between ICCV 96029 and other cultivated chickpea accessions are vernalization-dependent. In addition, we highlight the role of FTa3, an ortholog of FLOWERING LOCUS T in Arabidopsis, in the vernalization response of cultivated chickpea. A common set of differentially expressed genes was found for all comparisons between wild species and cultivars. The direction of expression change for different copies of the FT-INTERACTING PROTEIN 1 gene was variable in different comparisons, which suggests complex mechanisms of FT protein transport. Our study makes a contribution to the understanding of flowering time control in Cicer, and can provide genetic strategies to further improve this important agronomic trait.

Molecular Breeding and Drought Tolerance in Chickpea

Cicer arietinum L. is the third greatest widely planted imperative pulse crop worldwide, and it belongs to the Leguminosae family. Drought is the utmost common abiotic factor on plants, distressing their water status and limiting their growth and development. Chickpea genotypes have the natural ability to fight drought stress using certain strategies viz., escape, avoidance and tolerance. Assorted breeding methods, including hybridization, mutation, and marker-aided breeding, genome sequencing along with omics approaches, could be used to improve the chickpea germplasm lines(s) against drought stress. Root features, for instance depth and root biomass, have been recognized as the greatest beneficial morphological factors for managing terminal drought tolerance in the chickpea. Marker-aided selection, for example, is a genomics-assisted breeding (GAB) strategy that can considerably increase crop breeding accuracy and competence. These breeding technologies, notably marker-assisted breeding, omics, and plant physiology knowledge, underlined the importance of chickpea breeding and can be used in future crop improvement programmes to generate drought-tolerant cultivars(s).

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.

A full-length transcriptome and expression atlas of protein-coding genes and long non-coding RNAs is generated in chickpea. Components of transcriptional regulatory networks and candidate tissue-specific transcripts associated with quantitative trait loci are identified.

Herbal Arsenal against Skin Ailments: A Review Supported by In Silico Molecular Docking Studies

Maintaining healthy skin is important for a healthy body. At present, skin diseases are numerous, representing a major health problem affecting all ages from neonates to the elderly worldwide. Many people may develop diseases that affect the skin, including cancer, herpes, and cellulitis. Long-term conventional treatment creates complicated disorders in vital organs of the body. It also imposes socioeconomic burdens on patients. Natural treatment is cheap and claimed to be safe. The use of plants is as old as mankind. Many medicinal plants and their parts are frequently used to treat these diseases, and they are also suitable raw materials for the production of new synthetic agents. A review of some plant families, viz., Fabaceae, Asteraceae, Lamiaceae, etc., used in the treatment of skin diseases is provided with their most common compounds and in silico studies that summarize the recent data that have been collected in this area.

Aquaculture omics: An update on the current status of research and data analysis

Aquaculture is the fast-growing agricultural sector and has the ability to meet the growing demand for protein nutritional security for future population. In future aquaculture is going to be the major source of fish proteins as capture fisheries reached at its maximum. However, several challenges need to overcome such as lack of genetically improved strains/varieties, lack of species-specific feed/functional feed, round the year availability of quality fish seed, pollution of ecosystems and increased frequencies of disease occurrence etc. In recent years, the continuous development of high throughput sequencing technology has revolutionized the biological sciences and provided necessary tools. Application of 'omics' in aquaculture research have been successfully used to resolve several productive and reproductive issues and thus ensure its sustainability and profitability. To date, high quality draft genomes of over fifty fish species have been generated and successfully used to develop large number of single nucleotide polymorphism markers (SNPs), marker panels and other genomic resources etc in several aquaculture species. Similarly, transcriptome profiling and miRNAs analysis have been used in aquaculture research to identify key transcripts and expression analysis of candidate genes/miRNAs involved in reproduction, immunity, growth, development, stress toxicology and disease. Metagenome analysis emerged as a promising scientific tool to analyze the complex genomes contained within microbial communities. Metagenomics has been successfully used in the aquaculture sector to identify novel and potential pathogens, antibiotic resistance genes, microbial roles in microcosms, microbial communities forming biofloc, probiotics etc. In the current review, we discussed application of high-throughput technologies (NGS) in the aquaculture sector.

Chickpea Biofortification for Cytokinin Dehydrogenase via Genome Editing to Enhance Abiotic-Biotic Stress Tolerance and Food Security

Globally more than two billion people suffer from micronutrient malnutrition (also known as "hidden hunger"). Further, the pregnant women and children in developing nations are mainly affected by micronutrient deficiencies. One of the most important factors is food insecurity which can be mitigated by improving the nutritional values through biofortification using selective breeding and genetic enhancement techniques. Chickpea is the second most important legume with numerous economic and nutraceutical properties. Therefore, chickpea production needs to be increased from the current level. However, various kind of biotic and abiotic stresses hamper global chickpea production. The emerging popular targets for biofortification in agronomic crops include targeting cytokinin dehydrogenase (CKX). The CKXs play essential roles in both physiological and developmental processes and directly impact several agronomic parameters i.e., growth, development, and yield. Manipulation of CKX genes using genome editing tools in several crop plants reveal that CKXs are involved in regulation yield, shoot and root growth, and minerals nutrition. Therefore, CKXs have become popular targets for yield improvement, their overexpression and mutants can be directly correlated with the increased yield and tolerance to various stresses. Here, we provide detailed information on the different roles of CKX genes in chickpea. In the end, we discuss the utilization of genome editing tool clustered regularly interspaced short palindromic repeats/CRISPR associated protein 9 (CRISPR/Cas9) to engineer CKX genes that can facilitate trait improvement. Overall, recent advancements in CKX and their role in plant growth, stresses and nutrient accumulation are highlighted, which could be used for chickpea improvement.

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.

Insights into the gene and protein structures of the CaSWEET family members in chickpea (Cicer arietinum), and their gene expression patterns in different organs under various stress and abscisic acid treatments

'Sugars Will Eventually be Exported Transporters' (SWEETs) are a group of sugar transporters that play crucial roles in various biological processes, particularly plant stress responses. However, no information is available yet for the CaSWEET family in chickpea. Here, we identified all putative CaSWEET members in chickpea, and obtained their major characteristics, including physicochemical patterns, chromosomal distribution, subcellular localization, gene organization, conserved motifs and three-dimensional protein structures. Subsequently, we explored available transcriptome data to compare spatiotemporal transcript abundance of CaSWEET genes in various major organs. Finally, we studied the changes in their transcript levels in leaves and/or roots following dehydration and exogenous abscisic acid treatments using RT-qPCR to obtain valuable information underlying their potential roles in chickpea responses to water-stress conditions. Our results provide the first insights into the characteristics of the CaSWEET family members and a foundation for further functional characterizations of selected candidate genes for genetic engineering of chickpea.

Transcriptomic Analysis of Early Flowering Signals in ‘Royal’ Flax

Canada is one of the world’s leading producers and exporters of flax seed, with most production occurring in the Prairie Provinces. However, reduced season length and risk of frost restricts production in the northern grain belt of the Canadian Prairies. To expand the growing region of flax and increase production in Canada, flax breeders need to develop earlier-flowering varieties capable of avoiding the risk of abiotic stress. A thorough understanding of flowering control of flax is essential for the efficient breeding of such lines. We identified 722 putative flax flowering genes that span all major flowering-time pathways. Frequently, we found multiple flax homologues for a single Arabidopsis flowering gene. We used RNA sequencing to quantify the expression of genes in the shoot apical meristem (SAM) at 10, 15, 19, and 29 days after planting (dap) using the ‘Royal’ cultivar. We observed the expression of 80% of putative flax flowering genes and the differential expression of only 30%; these included homologues of major flowering regulators, such as SOC1, FUL, and AP1. We also found enrichment of differentially expressed genes (DEGs) in transcription factor (TF) families involved in flowering. Finally, we identified the candidates’ novel flowering genes amongst the uncharacterized flax genes. Our transcriptomic dataset provides a useful resource for investigating the regulatory control of the transition to flowering in flax and for the breeding of northern-adapted varieties.

Transcriptome Analysis Reveals Key Pathways and Candidate Genes Controlling Seed Development and Size in Ricebean (Vigna umbellata)

Ricebean (Vigna umbellata) is a lesser known pulse with well-recognized potential. Recently, it has emerged as a legume with endowed nutritional potential because of high concentration of quality protein and other vital nutrients in its seeds. However, the genes and pathways involved in regulating seed development and size are not understood in this crop. In our study, we analyzed the transcriptome of two genotypes with contrasting grain size (IC426787: large seeded and IC552985: small seeded) at two different time points, namely, 5 and 10 days post-anthesis (DPA). The bold seeded genotype across the time points (B5_B10) revealed 6,928 differentially expressed genes (DEGs), whereas the small seeded genotype across the time point (S5_S10) contributed to 14,544 DEGs. We have also identified several candidate genes for seed development–related traits like seed size and 100-seed weight. On the basis of similarity search and domain analysis, some candidate genes (PHO1, cytokinin dehydrogenase, A-type cytokinin, and ARR response negative regulator) related to 100-seed weight and seed size showed downregulation in the small seeded genotype. The MapMan and KEGG analysis confirmed that auxin and cytokinin pathways varied in both the contrasting genotypes and can therefore be the regulators of the seed size and other seed development–related traits in ricebeans. A total of 51 genes encoding SCF^TIR1/AFB, Aux/IAA, ARFs, E3 ubiquitin transferase enzyme, and 26S proteasome showing distinct expression dynamics in bold and small genotypes were also identified. We have also validated randomly selected SSR markers in eight accessions of the Vigna species (V. umbellata: 6; Vigna radiata: 1; and Vigna mungo: 1). Cross-species transferability pattern of ricebean–derived SSR markers was higher in V. radiata (73.08%) than V. mungo (50%). To the best of our knowledge, this is the first transcriptomic study conducted in this crop to understand the molecular basis of any trait. It would provide us a comprehensive understanding of the complex transcriptome dynamics during the seed development and gene regulatory mechanism of the seed size determination in ricebeans.

Stage specific comparative transcriptomic analysis to reveal gene networks regulating iron and zinc content in pearl millet [Pennisetum glaucum (L.) R. Br.]

Pearl millet is an important staple food crop of poor people and excels all other cereals due to its unique features of resilience to adverse climatic conditions. It is rich in micronutrients like iron and zinc and amenable for focused breeding for these micronutrients along with high yield. Hence, this is a key to alleviate malnutrition and ensure nutritional security. This study was conducted to identify and validate candidate genes governing grain iron and zinc content enabling the desired modifications in the genotypes. Transcriptome sequencing using ION S5 Next Generation Sequencer generated 43.5 million sequence reads resulting in 83,721 transcripts with N₅₀ of 597 bp and 84.35% of transcripts matched with the pearl millet genome assembly. The genotypes having high iron and zinc showed differential gene expression during different stages. Of which, 155 were up-regulated and 251 were down-regulated while during flowering stage and milking stage 349 and 378 transcripts were differentially expressed, respectively. Gene annotation and GO term showed the presence of transcripts involved in metabolic activities associated with uptake and transport of iron and zinc. Information generated will help in gaining insights into iron and zinc metabolism and develop genotypes with high yield, grain iron and zinc content.

Identification of Pueraria spp. through DNA barcoding and comparative transcriptomics

BACKGROUND

Kudzu is a term used generically to describe members of the genus Pueraria. Kudzu roots have been used for centuries in traditional Chinese medicine in view of their high levels of beneficial isoflavones including the unique 8-C-glycoside of daidzein, puerarin. In the US, kudzu is seen as a noxious weed causing ecological and economic damage. However, not all kudzu species make puerarin or are equally invasive. Kudzu remains difficult to identify due to its diverse morphology and inconsistent nomenclature.

RESULTS

We have generated sequences for the internal transcribed spacer 2 (ITS2) and maturase K (matK) regions of Pueraria montana lobata, P. montana montana, and P. phaseoloides, and identified two accessions previously used for differential analysis of puerarin biosynthesis as P. lobata and P. phaseoloides. Additionally, we have generated root transcriptomes for the puerarin-producing P. m. lobata and the non-puerarin producing P. phaseoloides. Within the transcriptomes, microsatellites were identified to aid in species identification as well as population diversity.

CONCLUSIONS

The barcode sequences generated will aid in fast and efficient identification of the three kudzu species. Additionally, the microsatellites identified from the transcriptomes will aid in genetic analysis. The root transcriptomes also provide a molecular toolkit for comparative gene expression analysis towards elucidation of the biosynthesis of kudzu phytochemicals.

Genomics Associated Interventions for Heat Stress Tolerance in Cool Season Adapted Grain Legumes

Cool season grain legumes occupy an important place among the agricultural crops and essentially provide multiple benefits including food supply, nutrition security, soil fertility improvement and revenue for farmers all over the world. However, owing to climate change, the average temperature is steadily rising, which negatively affects crop performance and limits their yield. Terminal heat stress that mainly occurred during grain development phases severely harms grain quality and weight in legumes adapted to the cool season, such as lentils, faba beans, chickpeas, field peas, etc. Although, traditional breeding approaches with advanced screening procedures have been employed to identify heat tolerant legume cultivars. Unfortunately, traditional breeding pipelines alone are no longer enough to meet global demands. Genomics-assisted interventions including new-generation sequencing technologies and genotyping platforms have facilitated the development of high-resolution molecular maps, QTL/gene discovery and marker-assisted introgression, thereby improving the efficiency in legumes breeding to develop stress-resilient varieties. Based on the current scenario, we attempted to review the intervention of genomics to decipher different components of tolerance to heat stress and future possibilities of using newly developed genomics-based interventions in cool season adapted grain legumes.

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.

Evolutionary and functional analysis of two-component system in chickpea reveals CaRR13, a TypeB RR, as positive regulator of symbiosis

The critical role of cytokinin in early nodulation in legumes is well known. In our study, exogenous cytokinin application to roots of the important crop legume, chickpea (Cicer arietinum L.), led to the formation of pseudo-nodules even in the absence of rhizobia. Hence, a genome-wide analysis of the cytokinin signalling, two-component system (TCS) genes, was conducted in chickpea, Medicago and Cajanus cajan. The integrated phylogenetic, evolutionary and expression analysis of the TCS genes was carried out, which revealed that histidine kinases (HKs) were highly conserved, whereas there was diversification leading to neofunctionalization at the level of response regulators (RRs) especially the TypeB RRs. Further, the functional role of the CaHKs in nodulation was established by complementation of the sln1Δ mutant of yeast and cre1 mutants of (Medicago) which led to restoration of the nodule-deficient phenotype. Additionally, the highest expressing TypeB RR of chickpea, CaRR13, was functionally characterized. Its localization in the nucleus and its Y1H assay-based interaction with the promoter of the early nodulation gene CaNSP2 indicated its role as a transcription factor regulating early nodulation. Overexpression, RNAi lines and complementation of cre1 mutants with CaRR13 revealed its critical involvement as an important signalling molecule regulating early events of nodule organogenesis in chickpea.

Genome-wide identification and functional prediction of salt- stress related long non-coding RNAs (lncRNAs) in chickpea (Cicer arietinum L.)

LncRNAs (long noncoding RNAs) are 200 bp length crucial RNA molecules, lacking coding potential and having important roles in regulating gene expression, particularly in response to abiotic stresses. In this study, we identified salt stress-induced lncRNAs in chickpea roots and predicted their intricate regulatory roles. A total of 3452 novel lncRNAs were identified to be distributed across all 08 chickpea chromosomes. On comparing salt-tolerant (ICCV 10, JG 11) and salt-sensitive cultivars (DCP 92-3, Pusa 256), 4446 differentially expressed lncRNAs were detected under various salt  treatments. We predicted 3373 lncRNAs to be regulating their target genes in cis regulating manner and 80 unique lncRNAs were observed as interacting with 136 different miRNAs, as eTMs (endogenous target mimic) targets of miRNAs and implicated them in the regulatory network of salt stress response. Functional analysis of these lncRNA revealed their association in targeting salt stress response-related genes like potassium transporter, transporter family genes, serine/threonine-protein kinase, aquaporins like TIP1-2, PIP2-5 and transcription factors like, AP2, NAC, bZIP, ERF, MYB and WRKY. Furthermore, about 614 lncRNA-SSRs (simple sequence repeats) were identified as a new generation of molecular markers with higher efficiency and specificity in chickpea. Overall, these findings will pave the understanding of comprehensive functional role of potential lncRNAs, which can help in providing insight into the molecular mechanism of salt tolerance in chickpea.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s12298-021-01093-0.

Fab Advances in Fabaceae for Abiotic Stress Resilience: From 'Omics' to Artificial Intelligence

Legumes are a better source of proteins and are richer in diverse micronutrients over the nutritional profile of widely consumed cereals. However, when exposed to a diverse range of abiotic stresses, their overall productivity and quality are hugely impacted. Our limited understanding of genetic determinants and novel variants associated with the abiotic stress response in food legume crops restricts its amelioration. Therefore, it is imperative to understand different molecular approaches in food legume crops that can be utilized in crop improvement programs to minimize the economic loss. 'Omics'-based molecular breeding provides better opportunities over conventional breeding for diversifying the natural germplasm together with improving yield and quality parameters. Due to molecular advancements, the technique is now equipped with novel 'omics' approaches such as ionomics, epigenomics, fluxomics, RNomics, glycomics, glycoproteomics, phosphoproteomics, lipidomics, regulomics, and secretomics. Pan-omics-which utilizes the molecular bases of the stress response to identify genes (genomics), mRNAs (transcriptomics), proteins (proteomics), and biomolecules (metabolomics) associated with stress regulation-has been widely used for abiotic stress amelioration in food legume crops. Integration of pan-omics with novel omics approaches will fast-track legume breeding programs. Moreover, artificial intelligence (AI)-based algorithms can be utilized for simulating crop yield under changing environments, which can help in predicting the genetic gain beforehand. Application of machine learning (ML) in quantitative trait loci (QTL) mining will further help in determining the genetic determinants of abiotic stress tolerance in pulses.

Systematic Identification and Validation of Housekeeping and Tissue-Specific Genes in Allotetraploid Chenopodium quinoa

Quinoa is a gluten-free food crop that contains all the essential amino acids and vitamins. The selection of proper housekeeping and tissue-specific genes is the crucial prerequisite for gene expression analysis using the common approach, real-time quantitative PCR (RT-qPCR). In this study, we identified 40 novel candidate housekeeping genes by the minimum transcript per million (TPM), coefficient of variation (CV) and maximum fold change (MFC) methods and 19 candidate tissue-specific genes by the co-expression network method based on an RNA-seq dataset that included 53 stem, leaf, flower and seed samples, as well as additional shoot and root samples under different stresses. The expression stability of 12 housekeeping and tissue-specific genes, as well as that of another two traditionally used housekeeping genes, was further evaluated using qPCR and ranked using NormFinder, BestKeeper and the comparative delta-Ct method. The results demonstrated that MIF, RGGA, VATE and UBA2B were ranked as the top four most stable candidate housekeeping genes. qPCR analysis also revealed three leaf-specific genes and five root-specific genes, but no stem-specific gene was identified. Gene Ontology (GO) enrichment analysis identified that housekeeping genes were mainly enriched in the small molecule metabolic process, organonitrogen compound metabolic process, NAD binding and ligase activity. In addition, tissue-specific genes are closely associated with the major functions of a specific tissue. Specifically, GO terms “photosynthesis” and “thylakoid” were most significantly overrepresented in candidate leaf-specific genes. The novel housekeeping and tissue-specific genes in our study will enable better normalization and quantification of transcript levels in quinoa.

The BBX gene family in Moso bamboo (Phyllostachys edulis): identification, characterization and expression profiles

BACKGROUND

The BBX (B-box) family are zinc finger protein (ZFP) transcription factors that play an essential role in plant growth, development and response to abiotic stresses. Although BBX genes have been characterized in many model organisms, genome-wide identification of the BBX family genes have not yet been reported in Moso bamboo (Phyllostachys edulis), and the biological functions of this family remain unknown.

RESULT

In the present study, we identified 27 BBX genes in the genome of Moso bamboo, and analysis of their conserved motifs and multiple sequence alignments revealed that they all shared highly similar structures. Additionally, phylogenetic and homology analyses indicated that PeBBX genes were divided into three clusters, with whole-genome duplication (WGD) events having facilitated the expansion of this gene family. Light-responsive and stress-related cis-elements were identified by analyzing cis-elements in the promoters of all PeBBX genes. Short time-series expression miner (STEM) analysis revealed that the PeBBX genes had spatiotemporal-specific expression patterns and were likely involved in the growth and development of bamboo shoots. We further explored the downstream target genes of PeBBXs, and GO/KEGG enrichment analysis predicted multiple functions of BBX target genes, including those encoding enzymes involved in plant photosynthesis, pyruvate metabolism and glycolysis/gluconeogenesis.

CONCLUSIONS

In conclusion, we analyzed the PeBBX genes at multiple different levels, which will contribute to further studies of the BBX family and provide valuable information for the functional validation of this family.

Full-Length Transcriptome Analysis of Four Different Tissues of Cephalotaxus oliveri

Cephalotaxus oliveri is a tertiary relict conifer endemic to China, regarded as a national second-level protected plant in China. This species has experienced severe changes in temperature and precipitation in the past millions of years, adapting well to harsh environments. In view of global climate change and its endangered conditions, it is crucial to study how it responds to changes in temperature and precipitation for its conservation work. In this study, single-molecule real-time (SMRT) sequencing and Illumina RNA sequencing were combined to generate the complete transcriptome of C. oliveri. Using the RNA-seq data to correct the SMRT sequencing data, the four tissues obtained 63,831 (root), 58,108 (stem), 33,013 (leaf) and 62,436 (male cone) full-length unigenes, with a N50 length of 2523, 3480, 3181, and 3267 bp, respectively. Additionally, 35,887, 11,306, 36,422, and 25,439 SSRs were detected for the male cone, leaf, root, and stem, respectively. The number of long non-coding RNAs predicted from the root was the largest (11,113), and the other tissues were 3408 (stem), 3193 (leaf), and 3107 (male cone), respectively. Functional annotation and enrichment analysis of tissue-specific expressed genes revealed the special roles in response to environmental stress and adaptability in the different four tissues. We also characterized the gene families and pathways related to abiotic factors. This work provides a comprehensive transcriptome resource for C. oliveri, and this resource will facilitate further studies on the functional genomics and adaptive evolution of C. oliveri.

Screening of mungbean for drought tolerance and transcriptome profiling between drought-tolerant and susceptible genotype in response to drought stress

Mungbean, is a widely cultivated pulse crop in India, experiences severe drought stress during the cultivation period. The mechanism of drought tolerance in mungbean is not well understood. In this presents study we screened 7 widely cultivated mungbean genotypes towards their drought sensitivity at seedling stage and transcriptome sequencing of drought-tolerant and susceptible genotype to understand the drought tolerance mechanism. Our physiological data such as increase in root length, shoot length, fresh weight, dry weight, relative water content (RWC), proline content, MDA content and molecular data in terms of quantitative expression of drought stress responsive genes under 3-d drought stress in mungbean suggests that, K851 seems to be most drought tolerant and PDM-139 as drought susceptible genotype. The transcriptomic study between K-851 and PDM-139 revealed 22,882 differentially expressed genes (DEGs) which were identified under drought stress, and they were mainly mapped to phytohormone signal transduction, carbon metabolism and flavonoid biosynthesis. Out of these, 10,235 genes were up-regulated and 12,647 genes were down-regulated. Furthermore, we found that, the DEGs related to, phytohormone signal transduction, carbon metabolism and flavonoid biosynthesis and they were more induced in K-851. Our data suggested that, the drought tolerant genotype K-851, scavenges the damage of drought stress by producing more amount of osmolytes, ROS scavengers and sugar biosynthesis.

Root-specific expression of chickpea cytokinin oxidase/dehydrogenase 6 leads to enhanced root growth, drought tolerance and yield without compromising nodulation

Cytokinin group of phytohormones regulate root elongation and branching during post-embryonic development. Cytokinin-degrading enzymes cytokinin oxidases/dehydrogenases (CKXs) have been deployed to investigate biological activities of cytokinin and to engineer root growth. We expressed chickpea cytokinin oxidase 6 (CaCKX6) under the control of a chickpea root-specific promoter of CaWRKY31 in Arabidopsis thaliana and chickpea having determinate and indeterminate growth patterns, respectively, to study the effect of cytokinin depletion on root growth and drought tolerance. Root-specific expression of CaCKX6 led to a significant increase in lateral root number and root biomass in Arabidopsis and chickpea without any penalty to vegetative and reproductive growth of shoot. Transgenic chickpea lines showed increased CKX activity in root. Soil-grown advanced chickpea transgenic lines exhibited higher root-to-shoot biomass ratio and enhanced long-term drought tolerance. These chickpea lines were not compromised in root nodulation and nitrogen fixation. The seed yield in some lines was up to 25% higher with no penalty in protein content. Transgenic chickpea seeds possessed higher levels of zinc, iron, potassium and copper. Our results demonstrated the potential of cytokinin level manipulation in increasing lateral root number and root biomass for agronomic trait improvement in an edible legume crop with indeterminate growth habit.

Identification of Alternative Mitochondrial Electron Transport Pathway Components in Chickpea Indicates a Differential Response to Salinity Stress between Cultivars

All plants contain an alternative electron transport pathway (AP) in their mitochondria, consisting of the alternative oxidase (AOX) and type 2 NAD(P)H dehydrogenase (ND) families, that are thought to play a role in controlling oxidative stress responses at the cellular level. These alternative electron transport components have been extensively studied in plants like Arabidopsis and stress inducible isoforms identified, but we know very little about them in the important crop plant chickpea. Here we identify AP components in chickpea (Cicer arietinum) and explore their response to stress at the transcript level. Based on sequence similarity with the functionally characterized proteins of Arabidopsis thaliana, five putative internal (matrix)-facing NAD(P)H dehydrogenases (CaNDA1-4 and CaNDC1) and four putative external (inter-membrane space)-facing NAD(P)H dehydrogenases (CaNDB1-4) were identified in chickpea. The corresponding activities were demonstrated for the first time in purified mitochondria of chickpea leaves and roots. Oxidation of matrix NADH generated from malate or glycine in the presence of the Complex I inhibitor rotenone was high compared to other plant species, as was oxidation of exogenous NAD(P)H. In leaf mitochondria, external NADH oxidation was stimulated by exogenous calcium and external NADPH oxidation was essentially calcium dependent. However, in roots these activities were low and largely calcium independent. A salinity experiment with six chickpea cultivars was used to identify salt-responsive alternative oxidase and NAD(P)H dehydrogenase gene transcripts in leaves from a three-point time series. An analysis of the Na:K ratio and Na content separated these cultivars into high and low Na accumulators. In the high Na accumulators, there was a significant up-regulation of CaAOX1, CaNDB2, CaNDB4, CaNDA3 and CaNDC1 in leaf tissue under long term stress, suggesting the formation of a stress-modified form of the mitochondrial electron transport chain (mETC) in leaves of these cultivars. In particular, stress-induced expression of the CaNDB2 gene showed a striking positive correlation with that of CaAOX1 across all genotypes and time points. The coordinated salinity-induced up-regulation of CaAOX1 and CaNDB2 suggests that the mitochondrial alternative pathway of respiration is an important facet of the stress response in chickpea, in high Na accumulators in particular, despite high capacities for both of these activities in leaf mitochondria of non-stressed chickpeas.

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.

Legume genomics and transcriptomics: From classic breeding to modern technologies

Legumes are essential and play a significant role in maintaining food standards and augmenting physiochemical soil properties through the biological nitrogen fixation process. Biotic and abiotic factors are the main factors limiting legume production. Classical breeding methodologies have been explored extensively about the problem of truncated yield in legumes but have not succeeded at the desired rate. Conventional breeding improved legume genotypes but with more resources and time. Recently, the invention of next-generation sequencing (NGS) and high-throughput methods for genotyping have opened new avenues for research and developments in legume studies. During the last decade, genome sequencing for many legume crops documented. Sequencing and re-sequencing of important legume species have made structural variation and functional genomics conceivable. NGS and other molecular techniques such as the development of markers; genotyping; high density genetic linkage maps; quantitative trait loci (QTLs) identification, expressed sequence tags (ESTs), single nucleotide polymorphisms (SNPs); and transcription factors incorporated into existing breeding technologies have made possible the accurate and accelerated delivery of information for researchers. The application of genome sequencing, RNA sequencing (transcriptome sequencing), and DNA sequencing (re-sequencing) provide considerable insights for legume development and improvement programs.

Moreover, RNA-Seq helps to characterize genes, including differentially expressed genes, and can be applied for functional genomics studies, especially when there is limited information available for the studied genomes. Genome-based crop development studies and the availability of genomics data as well as decision-making gears look be specific for breeding programs. This review mainly presents an overview of the path from classical breeding to new emerging genomics tools, which will trigger and accelerate genomics-assisted breeding for recognition of novel genes for yield and quality characters for sustainable legume crop production.

Transcript profiling of chickpea pod wall revealed the expression of floral homeotic gene AGAMOUS-like X2 (CaAGLX2)

The differentially expressed genes in the chickpea pod wall have been identified for the first time using a forward suppression subtractive hybridization (SSH) library. In all, 226 clones of SSH library were sequenced and analyzed. A total of 179 high-quality expressed sequence tags (ESTs) were generated and based on the CAP3 assembly of these ESTs, 126 genes (97 singletons and 29 contigs) were computationally annotated. The mapping of 88.26% ESTs by gene ontology (GO) annotation distributed them into 751 GO terms of three categories, cellular location, molecular function, and biological process. The KEGG pathway analysis revealed 45 ESTs are involved in 49 different biological pathways. Also, 67 ESTs encodes four different classes of enzymes such as oxidoreductases (29), transferase (20), hydrolases (16) and isomerase (2). Six genes were selected and subjected to qPCR analysis, of these, two genes (FHG Floral homeotic AGAMOUS-like isoform X2, MADS1 MADS-box transcription factor) showed significant up-regulation in the pod wall compared to leaves. Surprisingly, one of the MADS1 box gene, FHG (CaAGLX2), responsible for flower development expressed in the pod wall. Therefore, understanding its specific role in the pod wall could be interesting. Thus, the transcript dynamics of the chickpea pod wall revealed differentially expressed genes in the pod wall, which may be participating in the metabolic build-up of both pod wall and seeds.

Transcriptome sequencing, de novo assembly, characterisation of wild accession of blackgram (Vigna mungo var. silvestris) as a rich resource for development of molecular markers and validation of SNPs by high resolution melting (HRM) analysis

BACKGROUND

Blackgram [Vigna mungo (L.) Hepper], is an important legume crop of Asia with limited genomic resources. We report a comprehensive set of genic simple sequence repeat (SSR) and single nucleotide polymorphism (SNPs) markers using Illumina MiSeq sequencing of transcriptome and its application in genetic variation analysis and mapping.

RESULTS

Transcriptome sequencing of immature seeds of wild blackgram, V. mungo var. silvestris by Illumina MiSeq technology generated 1.9 × 10⁷ reads, which were assembled into 40,178 transcripts (TCS) with an average length of 446 bp covering 2.97 GB of the genome. A total of 38,753 CDS (Coding sequences) were predicted from 40,178 TCS and 28,984 CDS were annotated through BLASTX and mapped to GO and KEGG database resulting in 140 unique pathways. The tri-nucleotides were most abundant (39.9%) followed by di-nucleotide (30.2%). About 60.3 and 37.6% of SSR motifs were present in the coding sequences (CDS) and untranslated regions (UTRs) respectively. Among SNPs, the most abundant substitution type were transitions (Ts) (61%) followed by transversions (Tv) type (39%), with a Ts/Tv ratio of 1.58. A total of 2306 DEGs were identified by RNA Seq between wild and cultivar and validation was done by quantitative reverse transcription polymerase chain reaction. In this study, we genotyped SNPs with a validation rate of 78.87% by High Resolution Melting (HRM) Assay.

CONCLUSION

In the present study, 1621genic-SSR and 1844 SNP markers were developed from immature seed transcriptome sequence of blackgram and 31 genic-SSR markers were used to study genetic variations among different blackgram accessions. Above developed markers contribute towards enriching available genomic resources for blackgram and aid in breeding programmes.

An ethnobotanical study of medicinal plants used to treat skin diseases in northern Pakistan

Background

Skin diseases are a major health concern especially in association with human immune deficiency syndrome and acquired an immune deficiency. The aim of this study was to document the ethnomedicinal information of plants used to treat skin diseases in Northern Pakistan. This is the first quantitative ethnobotanical study of therapeutic herbs utilized by the indigenous people of Northern Pakistan for skin diseases.

Methods

Interviews were taken to obtain information from 180 participants. Quantitative methods including fidelity level (FL), Frequency of citation (FC), Use-value (UV), Jaccard indices (JI), Family importance value (FIV), Relative frequency of citation (RFC) and Chi-square test were applied. Medicinal plants uses are also compared with 50 national and international publications.

Results

In this study, we recorded 106 plant species belonged to 56 floral families for treatment of skin ailments. The dominant life form reported was herb while the preferred method of utilization was powder, along with leaf as the most used plant part. RFC ranges from 0.07 to 0.25% whereas the highest FIV was recorded for family Pteridaceae. FL values range from 36.8 to 100%. The study reported 88% of new plant reports for the treatment of skin diseases.

Conclusion

The present study revealed the importance of several plants used to treat skin diseases by the local communities of Northern Pakistan. The available literature supported the evidence of plant dermatological properties. Plants having high UV and RFC can be considered for further scientific analysis. There is dire need to create awareness among local, government and scientific communities for the preservation of medicinal species and ethnomedicinal knowledge in Northern Pakistan.

Electronic supplementary material

The online version of this article (10.1186/s12906-019-2605-6) contains supplementary material, which is available to authorized users.

The Transcriptome Landscape of Walnut Interspecies Hybrid (Juglans hindsii × Juglans regia) and Regulation of Cambial Activity in Relation to Grafting

Walnuts (Juglans, Juglandaceae) are known throughout the world as economically important trees that provide fat, protein, vitamins, and minerals as a food source, and produce high-quality timber. We have amended the purpose section to say "However," the omics resources are limited, which hampered the elucidation of molecular mechanisms resulting in their economically important traits (such as yield, fertility alternation, oil synthesis, and wood formation). To enrich the omics database of walnut, there is great need for analyses of its genomic and transcriptomic characteristics. In this study, we reported for the first time of the transcriptome landscape of six important organs or tissues in walnut interspecies hybrid using next-generation sequencing technology. Over 338 million clean reads were obtained. This yielded 74,072 unigenes with an average length of 782.71 bp. To develop an understanding of gene functions and regulatory pathways, 66,355 of the unigenes were identified as homologs of annotated genes and classified into three general categories with 61 functional subcategories. 2,288 out of 2,549 unmapped unigenes had at least one BLAST hit against the public databases. A total of 1,237 transcription factor-encoding genes (TFs) and 2,297 tissue-specific unigenes were identified. Interestingly, in the new shoot between an adult seedling and a grafted tree, the expression of 9,494 unigenes were significantly different, among which 4,388 were up-regulated and 5,106 were down-regulated. Of these, 195, 177, 232, 75, 114, and 68 unigenes were related to transcription factors, cell wall, defense response, transport, plant hormone biosynthesis, and other cambial activity-related functions, respectively. The obtained sequences and putative functional data constitute a resource for future functional analyses in walnut and other woody plants. These findings will be useful in further studies addressing the molecular mechanisms underlying grafting-related cambial activity.

Genomic structure and expression of alternative oxidase genes in legumes

Mitochondria isolated from chickpea (Cicer arietinum) possess substantial alternative oxidase (AOX) activity, even in non-stressed plants, and one or two AOX protein bands were detected immunologically, depending on the organ. Four different AOX isoforms were identified in the chickpea genome: CaAOX1 and CaAOX2A, B and D. CaAOX2A was the most highly expressed form and was strongly expressed in photosynthetic tissues, whereas CaAOX2D was found in all organs examined. These results are very similar to those of previous studies with soybean and siratro. Searches of available databases showed that this pattern of AOX genes and their expression was common to at least 16 different legume species. The evolution of the legume AOX gene family is discussed, as is the in vivo impact of an inherently high AOX capacity in legumes on growth and responses to environmental stresses.

Analysis of genes encoding seed storage proteins (SSPs) in chickpea (Cicer arietinum L.) reveals co-expressing transcription factors and a seed-specific promoter

Improvement of the quality and quantity of chickpea seed protein can be greatly facilitated by an understanding of the genic organization and the genetic architecture of the genes encoding seed storage proteins (SSPs). The aim of this study was to provide a comprehensive analysis of the chickpea SSP genes, putative co-expressing transcription factors (TFs), and to identify a seed-specific SSP gene promoter. A genome-wide identification of SSP genes in chickpea led to the identification of 21 non-redundant SSP encoding genes located on 6 chromosomes. Phylogenetic analysis grouped SSP genes into 3 subgroups where members within the same clade demonstrated similar motif composition and intron-exon organization. Tandem duplications were identified to be the major contributors to the expansion of the SSP gene family in chickpea. Co-expression analysis revealed 14 TFs having expression profiles similar to the SSP genes that included members of important TF families that are known to regulate seed development. Expression analysis of SSP genes and TFs revealed significantly higher expression in late stages of seed development as well as in high seed protein content (HPC) genotypes. In silico analysis of the promoter regions of the SSP encoding genes revealed several seed-specific cis-regulatory elements such as RY repeats, ACGT motifs, CAANTG, and GCN4. A candidate promoter was analyzed for seed specificity by generating stable transgenics in Arabidopsis. Overall, this study provides a useful resource to explore the regulatory networks involved in SSP synthesis and/or accumulation for utilization in developing nutritionally improved chickpea genotypes.

Identification, Structural Characterization and Gene Expression Analysis of Members of the Nuclear Factor-Y Family in Chickpea (Cicer arietinum L.) under Dehydration and Abscisic Acid Treatments

In plants, the Nuclear Factor-Y (NF-Y) transcription factors (TFs), which include three distinct types of NF-YA, NF-YB, and NF-YC TFs, have been identified to play key roles in the regulation of various plant growth and developmental processes under both normal and environmental stress conditions. In this work, a total of 40 CaNF-Y-encoding genes, including eight CaNF-YAs, 21 CaNF-YBs, and 11 CaNF-YCs, were identified in chickpea, and their major gene and protein characteristics were subsequently obtained using various web-based tools. Of our interest, a phylogenetically-based analysis predicted 18 CaNF-Ys (eight CaNF-YAs, seven CaNF-YBs, and three CaNF-YCs) that potentially play roles in chickpea responses to dehydration according to their close relationship with the well-characterized GmNF-Ys in soybean. These results were in good agreement with the enrichment of drought-responsive cis-regulatory motifs and expression patterns obtained from in silico analyses using publically available transcriptome data. Most of the phylogenetically predicted drought-responsive CaNF-Y genes (15 of 18) were quantitatively validated to significantly respond to dehydration treatment in leaves and/or roots, further supporting the results of in silico analyses. Among these CaNF-Y genes, the transcript levels of CaNF-YA01 and CaNF-YC10 were the most highly accumulated in leaves (by approximately eight-fold) and roots (by approximately 18-fold), respectively, by dehydration. Furthermore, 12 of the 18 CaNF-Y genes were found to be responsive to the most well-known stress hormone, namely abscisic acid (ABA), in leaves and/or roots, suggesting that these genes may act in chickpea response to dehydration in ABA-dependent manner. Taken together, our study has provided a comprehensive and fundamental information for further functional analyses of selected CaNF-Y candidate genes, ultimately leading to the improvement of chickpea growth under water-limited conditions.

Genome-wide transcriptome profiling of transgenic hop (Humulus lupulus L.) constitutively overexpressing HlWRKY1 and HlWDR1 transcription factors

Background

The hop plant (Humulus lupulus L.) is a valuable source of several secondary metabolites, such as flavonoids, bitter acids, and essential oils. These compounds are widely implicated in the beer brewing industry and are having potential biomedical applications. Several independent breeding programs around the world have been initiated to develop new cultivars with enriched lupulin and secondary metabolite contents but met with limited success due to several constraints. In the present work, a pioneering attempt has been made to overexpress master regulator binary transcription factor complex formed by HlWRKY1 and HlWDR1 using a plant expression vector to enhance the level of prenylflavonoid and bitter acid content in the hop. Subsequently, we performed transcriptional profiling using high-throughput RNA-Seq technology in leaves of resultant transformants and wild-type hop to gain in-depth information about the genome-wide functional changes induced by HlWRKY1 and HlWDR1 overexpression.

Results

The transgenic WW-lines exhibited an elevated expression of structural and regulatory genes involved in prenylflavonoid and bitter acid biosynthesis pathways. In addition, the comparative transcriptome analysis revealed a total of 522 transcripts involved in 30 pathways, including lipids and amino acids biosynthesis, primary carbon metabolism, phytohormone signaling and stress responses were differentially expressed in WW-transformants. It was apparent from the whole transcriptome sequencing that modulation of primary carbon metabolism and other pathways by HlWRKY1 and HlWDR1 overexpression resulted in enhanced substrate flux towards secondary metabolites pathway. The detailed analyses suggested that none of the pathways or genes, which have a detrimental effect on physiology, growth and development processes, were induced on a genome-wide scale in WW-transgenic lines.

Conclusions

Taken together, our results suggest that HlWRKY1 and HlWDR1 simultaneous overexpression positively regulates the prenylflavonoid and bitter acid biosynthesis pathways in the hop and thus these transgenes are presented as prospective candidates for achieving enhanced secondary metabolite content in the hop.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5125-8) contains supplementary material, which is available to authorized users.

Identification of genes associated with stress tolerance in moth bean [Vigna aconitifolia (Jacq.) Marechal], a stress hardy crop

Moth bean is the most drought and heat tolerant cultigens among Asian Vigna. We performed comparative transcriptome analysis of moth bean cultivar "Marumoth" under control and stress condition. De novo transcriptome assembly was carried out by using Velvet followed by Oases softwares. Differential expression analyses, SSR identification and validation and mapping of pathways and transcription factors were conducted. A total of 179,979 and 201,888 reads were generated on Roche 454 platform and 48,617,205 and 45,449,053 reads were generated on ABI Solid platform for the control and stressed samples. Combined assembly from Roche and ABI Solid platforms generated 16,090 and 15,096 transcripts for control and stressed samples. We found 1287 SSRs and 5606 transcripts involved in 179 pathways. The 55 transcription factor families represented 19.42% of total mothbean transcripts. In expression profiling, ten transcripts were found to be up-regulated and 41 down-regulated while 490 showed no major change under moisture stress condition. Stress inducible genes like Catalase, Cyt P450 monooxygenase, heat shock proteins (HSP 90 and HSP 70), oxidoreductase, protein kinases, dehydration responsive protein (DRP), universal stress protein and ferridoxin NADH oxidoreductase genes were up-regulated in stressed sample. Genes which might be involved in moisture stress tolerance in moth bean were identified and these might be useful for stress tolerance breeding in moth bean and other related crops.

Current advances in chickpea genomics: applications and future perspectives

Chickpea genomics promises to illuminate our understanding of genome organization, structural variations, evolutionary and domestication-related insights and fundamental biology of legume crops. Unprecedented advancements of next generation sequencing (NGS) technologies have enabled in decoding of multiple chickpea genome sequences and generating huge genomic resources in chickpea both at functional and structural level. This review is aimed to update the current progress of chickpea genomics ranging from high density linkage map development, genome-wide association studies (GWAS), functional genomics resources for various traits, emerging role of abiotic stress responsive coding and non-coding RNAs after the completion of draft chickpea genome sequences. Additionally, the current efforts of whole genome re-sequencing (WGRS) approach of global chickpea germplasm to capture the global genetic diversity existing in the historically released varieties across the world and increasing the resolution of the previously identified candidate gene(s) of breeding importance have been discussed. Thus, the outcomes of these genomics resources will assist in genomics-assisted selection and facilitate breeding of climate-resilient chickpea cultivars for sustainable agriculture.

Genome-wide bisulphite-sequencing reveals organ-specific methylation patterns in chickpea

DNA methylation is widely known to regulate gene expression in eukaryotes. Here, we unraveled DNA methylation patterns in cultivated chickpea to understand the regulation of gene expression in different organs. We analyzed the methylation pattern in leaf tissue of wild chickpea too, and compared it with cultivated chickpea. Our analysis indicated abundant CG methylation within gene-body and CHH methylation in intergenic regions of the chickpea genome in all the organs examined. Analysis of differentially methylated regions (DMRs) demonstrated a higher number of CG context DMRs in wild chickpea and CHH context DMRs in cultivated chickpea. We observed increased preponderance of hypermethylated DMRs in the promoter regions and hypomethylated DMRs in the genic regions in cultivated chickpea. Genomic location and context of the DMRs correlated well with expression of proximal genes. Our results put forth a positive correlation of promoter hypermethylation with increased transcript abundance via identification of DMR-associated genes involved in flower development in cultivated chickpea. The atypical correlation observed between promoter hypermethylation and increased transcript abundance might be dependent on 24-nt small RNAs and transcription factors binding to the promoter region. This study provides novel insights into DNA methylation patterns in chickpea and their role in regulation of gene expression.

Comprehensive metabolic and transcriptomic profiling of various tissues provide insights for saponin biosynthesis in the medicinally important Asparagus racemosus

Asparagus racemosus (Shatavari), belongs to the family Asparagaceae and is known as a “curer of hundred diseases” since ancient time. This plant has been exploited as a food supplement to enhance immune system and regarded as a highly valued medicinal plant in Ayurvedic medicine system for the treatment of various ailments such as gastric ulcers, dyspepsia, cardiovascular diseases, neurodegenerative diseases, cancer, as a galactogogue and against several other diseases. In depth metabolic fingerprinting of various parts of the plant led to the identification of 13 monoterpenoids exclusively present in roots. LC-MS profiling led to the identification of a significant number of steroidal saponins (33). However, we have also identified 16 triterpene saponins for the first time in A. racemosus. In order to understand the molecular basis of biosynthesis of major components, transcriptome sequencing from three different tissues (root, leaf and fruit) was carried out. Functional annotation of A. racemosus transcriptome resulted in the identification of 153 transcripts involved in steroidal saponin biosynthesis, 45 transcripts in triterpene saponin biosynthesis, 44 transcripts in monoterpenoid biosynthesis and 79 transcripts in flavonoid biosynthesis. These findings will pave the way for better understanding of the molecular basis of steroidal saponin, triterpene saponin, monoterpenoids and flavonoid biosynthesis in A. racemosus.

The aquatic animals' transcriptome resource for comparative functional analysis

Background

Aquatic animals have great economic and ecological importance. Among them, non-model organisms have been studied regarding eco-toxicity, stress biology, and environmental adaptation. Due to recent advances in next-generation sequencing techniques, large amounts of RNA-seq data for aquatic animals are publicly available. However, currently there is no comprehensive resource exist for the analysis, unification, and integration of these datasets. This study utilizes computational approaches to build a new resource of transcriptomic maps for aquatic animals. This aquatic animal transcriptome map database dbATM provides de novo assembly of transcriptome, gene annotation and comparative analysis of more than twenty aquatic organisms without draft genome.

Results

To improve the assembly quality, three computational tools (Trinity, Oases and SOAPdenovo-Trans) were employed to enhance individual transcriptome assembly, and CAP3 and CD-HIT-EST software were then used to merge these three assembled transcriptomes. In addition, functional annotation analysis provides valuable clues to gene characteristics, including full-length transcript coding regions, conserved domains, gene ontology and KEGG pathways. Furthermore, all aquatic animal genes are essential for comparative genomics tasks such as constructing homologous gene groups and blast databases and phylogenetic analysis.

Conclusion

In conclusion, we establish a resource for non model organism aquatic animals, which is great economic and ecological importance and provide transcriptomic information including functional annotation and comparative transcriptome analysis. The database is now publically accessible through the URL http://dbATM.mbc.nctu.edu.tw/.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4463-x) contains supplementary material, which is available to authorized users.

Genome-wide analysis identifies chickpea (Cicer arietinum) heat stress transcription factors (Hsfs) responsive to heat stress at the pod development stage

The heat stress transcription factors (Hsfs) play a prominent role in thermotolerance and eliciting the heat stress response in plants. Identification and expression analysis of Hsfs gene family members in chickpea would provide valuable information on heat stress responsive Hsfs. A genome-wide analysis of Hsfs gene family resulted in the identification of 22 Hsf genes in chickpea in both desi and kabuli genome. Phylogenetic analysis distinctly separated 12 A, 9 B, and 1 C class Hsfs, respectively. An analysis of cis-regulatory elements in the upstream region of the genes identified many stress responsive elements such as heat stress elements (HSE), abscisic acid responsive element (ABRE) etc. In silico expression analysis showed nine and three Hsfs were also expressed in drought and salinity stresses, respectively. Q-PCR expression analysis of Hsfs under heat stress at pod development and at 15 days old seedling stage showed that CarHsfA2, A6, and B2 were significantly upregulated in both the stages of crop growth and other four Hsfs (CarHsfA2, A6a, A6c, B2a) showed early transcriptional upregulation for heat stress at seedling stage of chickpea. These subclasses of Hsfs identified in this study can be further evaluated as candidate genes in the characterization of heat stress response in chickpea.

Transcriptome profiling of genes involving in carotenoid biosynthesis and accumulation between leaf and root of carrot (Daucus carota L.)

Carrot provides abundant carotenoid for human diet and is one of the most widely cultivated root vegetables in the world. However, the absence of the tissue-specific transcriptome of carrots hampers the investigation of the association of secondary metabolic mechanism with the different tissue types. In this study, we obtained 46,119,008/48,414,508 raw reads and 45,394,846/47,887,648 clean reads from the carrot leaf and root, respectively. Moreover, α- and β-carotene were found to accumulate in both tissues. Then, using Trinity assembly into contigs and mapped back to contigs, these reads were assembled to 56,267 and 62,427 leaf and root unigenes, respectively, after Ns removal and paired-end extension. In addition, a total of 18,354 DEGs were found between the carrot leaf and root unigenes, and 99 of these DEGs were found to be involved in carotenoid biosynthesis as revealed by integrated function annotation. In the carotenoid pathway DEGs, DcPSY1, DcZ-ISO, DcCISO2, DcLBCY, DcLECY, DcZEP1, DcZEP2, DcVDE1, DcVDE2, DcNSY1, DcNSY2, DcA8H-CYP707A1.2, DcAAO3a, DcCCD4, and DcMAX1 were expressed dramatically in the carrot leaf compared with in the root. This result was consistent with the results from the quantitative real-time PCR analysis of DEG expression profiles. Moreover, 67 more carotenoid biosynthesis-related genes were found in this transcriptome database. Most of these DEGs were up-regulated in the carrot leaf compared with those in the root. The expression of DEGs may be related to the higher carotenoid pathway flux in the carrot leaf than in the root. These results will help to further understand the carotenoid biosynthesis in carrot.