Comprehensive Stress-Based De Novo Transcriptome Assembly and Annotation of Guar (Cyamopsis tetragonoloba (L.) Taub.): An Important Industrial and Forage Crop

Development of de-novo transcriptome assembly and SSRs in allohexaploid Brassica with functional annotations and identification of heat-shock proteins for thermotolerance

Crop Brassicas contain monogenomic and digenomic species, with no evidence of a trigenomic Brassica in nature. Through somatic fusion (Sinapis alba + B. juncea), a novel allohexaploid trigenomic Brassica (H1 = AABBSS; 2n = 60) was produced and used for transcriptome analysis to uncover genes for thermotolerance, annotations, and microsatellite markers for future molecular breeding. Illumina Novaseq 6000 generated a total of 76,055,546 paired-end raw reads, which were used for de-novo assembly, resulting in the development of 486,066 transcripts. A total of 133,167 coding sequences (CDSs) were predicted from transcripts with a mean length of 507.12 bp and 46.15% GC content. The BLASTX search of CDSs against public protein databases showed a maximum of 126,131 (94.72%) and a minimum of 29,810 (22.39%) positive hits. Furthermore, 953,773 gene ontology (GO) terms were found in 77,613 (58.28%) CDSs, which were divided into biological processes (49.06%), cellular components (31.67%), and molecular functions (19.27%). CDSs were assigned to 144 pathways by a pathway study using the KEGG database and 1,551 pathways by a similar analysis using the Reactome database. Further investigation led to the discovery of genes encoding over 2,000 heat shock proteins (HSPs). The discovery of a large number of HSPs in allohexaploid Brassica validated our earlier findings for heat tolerance at seed maturity. A total of 15,736 SSRs have been found in 13,595 CDSs, with an average of one SSR per 4.29 kb length and an SSR frequency of 11.82%. The first transcriptome assembly of a meiotically stable allohexaploid Brassica has been given in this article, along with functional annotations and the presence of SSRs, which could aid future genetic and genomic studies.

Transcriptome Analysis Provides Insights into Potentilla bifurca Adaptation to High Altitude

Potentilla bifurca is widely distributed in Eurasia, including the Tibetan Plateau. It is a valuable medicinal plant in the Tibetan traditional medicine system, especially for the treatment of diabetes. This study investigated the functional gene profile of Potentilla bifurca at different altitudes by RNA-sequencing technology, including de novo assembly of 222,619 unigenes from 405 million clean reads, 57.64% of which were annotated in Nr, GO, KEGG, Pfam, and Swiss-Prot databases. The most significantly differentially expressed top 50 genes in the high-altitude samples were derived from plants that responded to abiotic stress, such as peroxidase, superoxide dismutase protein, and the ubiquitin-conjugating enzyme. Pathway analysis revealed that a large number of DEGs encode key enzymes involved in secondary metabolites, including phenylpropane and flavonoids. In addition, a total of 298 potential genomic SSRs were identified in this study, which provides information on the development of functional molecular markers for genetic diversity assessment. In conclusion, this study provides the first comprehensive assessment of the Potentilla bifurca transcriptome. This provides new insights into coping mechanisms for non-model organisms surviving in harsh environments at high altitudes, as well as molecular evidence for the selection of superior medicinal plants.

De novo Transcriptome Analysis of Drought-Adapted Cluster Bean (Cultivar RGC-1025) Reveals the Wax Regulatory Genes Involved in Drought Resistance

Cluster bean (Cyamopsis tetragonoloba L.) is one of the multipurpose underexplored crops grown as green vegetable and for gum production in dryland areas. Cluster bean is known as relatively tolerant to drought and salinity stress. To elucidate the molecular mechanisms involved in the drought tolerance of cluster bean cultivar RGC-1025, RNA sequencing (RNA-seq) of the drought-stressed and control samples was performed. De novo assembly of the reads resulted in 66,838 transcripts involving 203 pathways. Among these transcripts, differentially expressed gene (DEG) analysis resulted in some of the drought-responsive genes expressing alpha dioxygenase 2, low temperature-induced 65 kDa protein (LDI65), putative vacuolar amino acid transporter, and late embryogenesis abundant protein (LEA 3). The analysis also reported drought-responsive transcription factors (TFs), such as NAC, WRKY, GRAS, and MYB families. The relative expression of genes by qRT-PCR revealed consistency with the DEG analysis. Key genes involved in the wax biosynthesis pathway were mapped using the DEG data analysis. These results were positively correlated with epicuticular wax content and the wax depositions on the leaf surfaces, as evidenced by scanning electron microscope (SEM) image analysis. Further, these findings support the fact that enhanced wax deposits on the leaf surface had played a crucial role in combating the drought stress in cluster beans under drought stress conditions. In addition, this study provided a set of unknown genes and TFs that could be a source of engineering tolerance against drought stress in cluster beans.

Transcriptome Analysis and Identification of a Female-Specific SSR Marker in Pistacia chinensis Based on Illumina Paired-End RNA Sequencing

Pistacia chinensis Bunge (P. chinensis), a dioecious plant species, has been widely found in China. The female P. chinensis plants are more important than male plants in agricultural production, as their seeds can serve as an ideal feedstock for biodiesel. However, the sex of P. chinensis plants is hard to distinguish during the seedling stage due to the scarcity of available transcriptomic and genomic information. In this work, Illumina paired-end RNA sequencing assay was conducted to unravel the transcriptomic profiles of female and male P. chinensis flower buds. In total, 50,925,088 and 51,470,578 clean reads were obtained from the female and male cDNA libraries, respectively. After quality checks and de novo assembly, a total of 83,370 unigenes with a mean length of 1.3 kb were screened. Overall, 64,539 unigenes (77.48%) could be matched in at least one of the NR, NT, Swiss-Prot, COG, KEGG, and GO databases, 71 of which were putatively related to the floral development of P. chinensis. Additionally, 21,662 simple sequence repeat (SSR) motifs were identified in 17,028 unigenes of P. chinensis, and the mononucleotide motif was the most dominant type of repeats (52.59%) in P. chinensis, followed by dinucleotide (22.29%), trinucleotide (20.15%). The most abundant repeats were AG/CT (13.97%), followed by AAC/GTT (6.75%) and AT/TA (6.10%). Based on these SSR, 983 EST-SSR primers were designed, 151 of which were randomly chosen for validation. Of these validated EST-SSR markers, 25 SSR markers were found to be polymorphic between male and female plants. One SSR marker, namelyPCSSR55, displayed excellent specificity in female plants, which could clearly distinguish between male and female P. chinensis. Altogether, our findings not only reveal that the EST-SSR marker is extremely effective in distinguishing between male and female P. chinensis but also provide a solid framework for sex determination of plant seedlings.

Reference nodule transcriptomes for Melilotus officinalis and Medicago sativa cv. Algonquin

Host/symbiont compatibility is a hallmark of the symbiotic nitrogen-fixing interaction between rhizobia and legumes, mediated in part by plant-produced nodule-specific cysteine-rich (NCR) peptides and the bacterial BacA membrane protein that can act as a NCR peptide transporter. In addition, the genetic and metabolic properties supporting symbiotic nitrogen fixation often differ between compatible partners, including those sharing a common partner, highlighting the need for multiple study systems. Here, we report high-quality nodule transcriptome assemblies for Medicago sativa cv. Algonquin and Melilotus officinalis, two legumes able to form compatible symbioses with Sinorhizobium meliloti. The compressed M. sativa and M. officinalis assemblies consisted of 79,978 and 64,593 contigs, respectively, of which 33,341 and 28,278 were assigned putative annotations, respectively. As expected, the two transcriptomes showed broad similarity at a global level. We were particularly interested in the NCR peptide profiles of these plants, as these peptides drive bacterial differentiation during the symbiosis. A total of 412 and 308 NCR peptides were predicted from the M. sativa and M. officinalis transcriptomes, respectively, with approximately 9% of the transcriptome of both species consisting of NCR transcripts. Notably, transcripts encoding highly cationic NCR peptides (isoelectric point > 9.5), which are known to have antimicrobial properties, were ∼2-fold more abundant in M. sativa than in M. officinalis, and ∼27-fold more abundant when considering only NCR peptides in the six-cysteine class. We hypothesize that the difference in abundance of highly cationic NCR peptides explains our previous observation that some rhizobial bacA alleles which can support symbiosis with M. officinalis are unable to support symbiosis with M. sativa.

Something Fishy about Siamese Fighting Fish (Betta splendens) Sex: Polygenic Sex Determination or a Newly Emerged Sex-Determining Region?

Fishes provide a unique and intriguing model system for studying the genomic origin and evolutionary mechanisms underlying sex determination and high sex-chromosome turnover. In this study, the mode of sex determination was investigated in Siamese fighting fish, a species of commercial importance. Genome-wide SNP analyses were performed on 75 individuals (40 males and 35 females) across commercial populations to determine candidate sex-specific/sex-linked loci. In total, 73 male-specific loci were identified and mapped to a 5.6 kb region on chromosome 9, suggesting a putative male-determining region (pMDR) containing localized dmrt1 and znrf3 functional sex developmental genes. Repeat annotations of the pMDR revealed an abundance of transposable elements, particularly Ty3/Gypsy and novel repeats. Remarkably, two out of the 73 male-specific loci were located on chromosomes 7 and 19, implying the existence of polygenic sex determination. Besides male-specific loci, five female-specific loci on chromosome 9 were also observed in certain populations, indicating the possibility of a female-determining region and the polygenic nature of sex determination. An alternative explanation is that male-specific loci derived from other chromosomes or female-specific loci in Siamese fighting fish recently emerged as new sex-determining loci during domestication and repeated hybridization.

Comparative transcriptomic analysis and antioxidant defense mechanisms in clusterbean (Cyamopsis tetragonoloba (L.) Taub.) genotypes with contrasting drought tolerance

To understand drought tolerance mechanism(s) in clusterbean (Cyamopsis tetragonoloba), we conducted physiological, biochemical, and de novo comparative transcriptome analysis of drought-tolerant (RGC-1002) and drought-sensitive (RGC-1066) genotypes subjected to 30 days of drought stress. Relative water content (RWC) was maintained in tolerant genotype but was reduced in sensitive genotype. Leaf pigment concentrations were higher in tolerant genotype. Net photosynthesis was significantly decreased in sensitive genotype but insignificant reduction was found in tolerant genotype. Enzymatic antioxidant (GR, APX, DHAR) activities were enhanced in tolerant genotype, while there were insignificant changes in these enzymes in sensitive genotype. The ratios of antioxidant molecules (ASC/DHA and GSH/GSSG) were higher in tolerant genotype as compared to sensitive genotype. In sensitive genotype, 6625 differentially expressed genes (DEGs) were upregulated and 5365 genes were downregulated. In tolerant genotype, 5206 genes were upregulated and 2793 genes were downregulated. In tolerant genotype, transketolase family protein, phosphoenolpyruvate carboxylase 3, temperature-induced lipocalin, and cytochrome oxidase were highly upregulated. Moreover, according to Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis, the drought tolerance may be attributed to upregulated starch and sucrose metabolism-related genes in tolerant genotype. Finally, quantitative real-time PCR confirmed the reproducibility of the RNA-seq data.

Progress of Genomics-Driven Approaches for Sustaining Underutilized Legume Crops in the Post-Genomic Era

Legume crops, belonging to the Fabaceae family, are of immense importance for sustaining global food security. Many legumes are profitable crops for smallholder farmers due to their unique ability to fix atmospheric nitrogen and their intrinsic ability to thrive on marginal land with minimum inputs and low cultivation costs. Recent progress in genomics shows promise for future genetic gains in major grain legumes. Still it remains limited in minor legumes/underutilized legumes, including adzuki bean, cluster bean, horse gram, lathyrus, red clover, urd bean, and winged bean. In the last decade, unprecedented progress in completing genome assemblies of various legume crops and resequencing efforts of large germplasm collections has helped to identify the underlying gene(s) for various traits of breeding importance for enhancing genetic gain and contributing to developing climate-resilient cultivars. This review discusses the progress of genomic resource development, including genome-wide molecular markers, key breakthroughs in genome sequencing, genetic linkage maps, and trait mapping for facilitating yield improvement in underutilized legumes. We focus on 1) the progress in genomic-assisted breeding, 2) the role of whole-genome resequencing, pangenomes for underpinning the novel genomic variants underlying trait gene(s), 3) how adaptive traits of wild underutilized legumes could be harnessed to develop climate-resilient cultivars, 4) the progress and status of functional genomics resources, deciphering the underlying trait candidate genes with putative function in underutilized legumes 5) and prospects of novel breeding technologies, such as speed breeding, genomic selection, and genome editing. We conclude the review by discussing the scope for genomic resources developed in underutilized legumes to enhance their production and play a critical role in achieving the "zero hunger" sustainable development goal by 2030 set by the United Nations.

Deciphering Molecular Mechanisms Involved in Salinity Tolerance in Guar (Cyamopsis tetragonoloba (L.) Taub.) Using Transcriptome Analyses

Guar is a commercially important legume crop known for guar gum. Guar is tolerant to various abiotic stresses, but the mechanisms involved in its salinity tolerance are not well established. This study aimed to understand molecular mechanisms of salinity tolerance in guar. RNA sequencing (RNA-Seq) was employed to study the leaf and root transcriptomes of salt-tolerant (Matador) and salt-sensitive (PI 340261) guar genotypes under control and salinity. Our analyses identified a total of 296,114 unigenes assembled from 527 million clean reads. Transcriptome analysis revealed that the gene expression differences were more pronounced between salinity treatments than between genotypes. Differentially expressed genes associated with stress-signaling pathways, transporters, chromatin remodeling, microRNA biogenesis, and translational machinery play critical roles in guar salinity tolerance. Genes associated with several transporter families that were differentially expressed during salinity included ABC, MFS, GPH, and P-ATPase. Furthermore, genes encoding transcription factors/regulators belonging to several families, including SNF2, C₂H₂, bHLH, C3H, and MYB were differentially expressed in response to salinity. This study revealed the importance of various biological pathways during salinity stress and identified several candidate genes that may be used to develop salt-tolerant guar genotypes that might be suitable for cultivation in marginal soils with moderate to high salinity or using degraded water.

Identification of Key Metabolic Pathways and Biomarkers Underlying Flowering Time of Guar (Cyamopsis tetragonoloba (L.) Taub.) via Integrated Transcriptome-Metabolome Analysis

Guar (Cyamopsis tetragonoloba (L.) Taub.) is an annual legume crop native to India and Pakistan. Seeds of the plant serve as a source of galactomannan polysaccharide (guar gum) used in the food industry as a stabilizer (E412) and as a gelling agent in oil and gas fracturing fluids. There were several attempts to introduce this crop to countries of more northern latitudes. However, guar is a plant of a short photoperiod, therefore, its introduction, for example, to Russia is complicated by a long day length during the growing season. Breeding of new guar varieties insensitive to photoperiod slowed down due to the lack of information on functional molecular markers, which, in turn, requires information on guar genome. Modern breeding strategies, e.g., genomic predictions, benefit from integration of multi-omics approaches such as transcriptome, proteome and metabolome assays. Here we present an attempt to use transcriptome-metabolome integration to understand the genetic determination of flowering time variation among guar plants that differ in their photoperiod sensitivity. This study was performed on nine early- and six delayed-flowering guar varieties with the goal to find a connection between 63 metabolites and 1,067 differentially expressed transcripts using Shiny GAM approach. For the key biomarker of flowering in guar myo-inositol we also evaluated the KEGG biochemical pathway maps available for Arabidopsis thaliana. We found that the phosphatidylinositol signaling pathway is initiated in guar plants that are ready for flowering through the activation of the phospholipase C (PLC) gene, resulting in an exponential increase in the amount of myo-inositol in its free form observed on GC-MS chromatograms. The signaling pathway is performed by suppression of myo-inositol phosphate kinases (phosphorylation) and alternative overexpression of phosphatases (dephosphorylation). Our study suggests that metabolome and transcriptome information taken together, provide valuable information about biomarkers that can be used as a tool for marker-assisted breeding, metabolomics and functional genomics of this important legume crop.

De novo transcriptome sequencing of triton shell Charonia lampas sauliae: Identification of genes related to neurotoxins and discovery of genetic markers

Charonia lampas sauliae (triton snails, triton shells or tritons; Mollusca, Caenogastropoda, Littorinimorpha, Ranellidae) is a marine species with a wide distribution. In Korea, this species is listed as vulnerable and is regionally protected as an endangered species. Here, we report the first comprehensive transcriptome dataset of C. lampas sauliae obtained using the Illumina HiSeq 2500 platform. In total, 97.68% of raw read sequences were processed as clean reads. Of the 577,478 contigs obtained, 146,026 sequences were predicted to contain coding regions. About 89.34% of all annotated unigene sequences showed homologous matches to protein sequences in PANM DB (Protostome database). Further, about one-third of the unigene sequences were annotated using the UniGene, Swiss-Prot, Clusters of Orthologous Groups (COG) and Gene Ontology (GO) databases. In total, 190 enzymes were predicted under key metabolic pathways under stood through Kyoto Encyclopedia of Genes and Genomes (KEGG) database annotation. Repetitive elements such as long terminal repeats (LTRs), short interspersed nuclear elements (SINEs), long interspersed nuclear elements (LINEs), and DNA elements were enriched in the unigene sequences. Among the identified transcripts were the channel proteins, some of which were blocked by tetrodotoxin, which is thought to be synthesized by symbiotic bacteria inhabiting the shells. In addition, conotoxin superfamily peptides, such as B-conotoxin, conotoxin superfamily T and alpha-conotoxin, were identified, which may have relevance to biomedical and evolutionary research. A transcriptome-wide search for polymorphic loci identified 21,568 simple sequence repeats (SSRs) in the unigene sequences. Most SSRs were dinucleotides, among which AC/GT was the dominant SSR type. The molecular and genetic resources revealed in this study could be utilized for investigations on the fitness of the species in the marine environment and sustainability in a changing habitat.

Transcriptome annotation in the cloud: complexity, best practices, and cost

BACKGROUND

The NIH Science and Technology Research Infrastructure for Discovery, Experimentation, and Sustainability (STRIDES) initiative provides NIH-funded researchers cost-effective access to commercial cloud providers, such as Amazon Web Services (AWS) and Google Cloud Platform (GCP). These cloud providers represent an alternative for the execution of large computational biology experiments like transcriptome annotation, which is a complex analytical process that requires the interrogation of multiple biological databases with several advanced computational tools. The core components of annotation pipelines published since 2012 are BLAST sequence alignments using annotated databases of both nucleotide or protein sequences almost exclusively with networked on-premises compute systems.

FINDINGS

We compare multiple BLAST sequence alignments using AWS and GCP. We prepared several Jupyter Notebooks with all the code required to submit computing jobs to the batch system on each cloud provider. We consider the consequence of the number of query transcripts in input files and the effect on cost and processing time. We tested compute instances with 16, 32, and 64 vCPUs on each cloud provider. Four classes of timing results were collected: the total run time, the time for transferring the BLAST databases to the instance local solid-state disk drive, the time to execute the CWL script, and the time for the creation, set-up, and release of an instance. This study aims to establish an estimate of the cost and compute time needed for the execution of multiple BLAST runs in a cloud environment.

CONCLUSIONS

We demonstrate that public cloud providers are a practical alternative for the execution of advanced computational biology experiments at low cost. Using our cloud recipes, the BLAST alignments required to annotate a transcriptome with ∼500,000 transcripts can be processed in <2 hours with a compute cost of ∼$200-$250. In our opinion, for BLAST-based workflows, the choice of cloud platform is not dependent on the workflow but, rather, on the specific details and requirements of the cloud provider. These choices include the accessibility for institutional use, the technical knowledge required for effective use of the platform services, and the availability of open source frameworks such as APIs to deploy the workflow.

Identification of Differentially Expressed Drought-Responsive Genes in Guar [Cyamopsis tetragonoloba (L.) Taub]

Drought remains one of the most serious environmental stresses because of the continuous reduction in soil moisture, which requires the improvement of crops with features such as drought tolerance. Guar [Cyamopsis tetragonoloba (L.) Taub], a forage and industrial crop, is a nonthirsty plant. However, the information on the transcriptome changes that occur under drought stress in guar is very limited; therefore, a gene expression analysis is necessary in this context. Here, we studied the differentially expressed genes (DEGs) in response to drought stress and their metabolic pathways. RNA-Seq via an expectation-maximization algorithm was used to estimate gene abundance. Subsequently, an Empirical Analysis of Digital Gene Expression Data in the R Bioconductor package was used to identify DEGs. Blast2GO, InterProScan, and the Kyoto Encyclopedia of Genes and Genomes were used to explore functional annotation, protein analysis, enzymes, and metabolic pathways. Transcription factors were identified using the PlantTFDB database. Our study identified 499 upregulated and 191 downregulated genes in response to drought stress. Of those, 32 upregulated and six downregulated genes were deemed as novel genes exclusive to guar. An aggregate of 137 protein families, 306 domains, 12 repeats, and two sites were upregulated. The proton-dependent oligopeptide transporter family and transferase, aquaporin transporter, calcium/calmodulin-dependent/calcium-dependent protein kinase, aspartic peptidase A1 family, UDP-glucuronosyl/UDP-glucosyltransferase, and major intrinsic protein were the most upregulated protein families. The upregulated unigenes were associated with 88 enzymes and 77 KEGG pathways. Finally, the MYB-related, MYB, and ERF transcription factor families were upregulated. These data may be useful for understanding the plant molecular response to drought stress.

Identification of Heat-Responsive Genes in Guar [Cyamopsis tetragonoloba (L.) Taub]

The threat of heat stress on crop production increased dramatically due to global warming leading to the rise on the demand of heat-tolerant crops and understanding their tolerance. The leguminous forage crop Guar [Cyamopsis tetragonoloba (L.) Taub] is a high-temperature tolerant plant with numerous works on its tolerance at morph-physiological levels but lack on molecular thermotolerance level. In the current study, the differential gene expression and the underlying metabolic pathways induced by heat treatment were investigated. An RNA-Seq study on Guar leaves was carried out to estimate gene abundance and identify genes involved in heat tolerance to better understand the response mechanisms to heat stress. The results uncovered 1551 up- and 1466 downregulated genes, from which 200 and 72 genes with unknown function could be considered as new genes specific to guar. The upregulated unigenes were associated with 158 enzymes and 102 KEGG pathways. Blast2GO, InterProScan, and Kyoto Encyclopaedia of Genes and Genomes packages were utilized to search the functional annotation, protein analysis, enzymes, and metabolic pathways and revealed hormone signal transduction were enriched during heat stress tolerance. A total of 301 protein families, 551 domains, 15 repeats, and 3 sites were upregulated and matched to those unigenes. A batch of heat-regulated transcription factor transcripts were identified using the PlantTFDB database, which may play roles in heat response in Guar. Interestingly, several heat shock protein families were expressed in response to exposure to stressful conditions for instance small HSP20, heat shock transcription factor family, heat shock protein Hsp90 family, and heat shock protein 70 family. Our results revealed the expressional changes associated with heat tolerance and identified potential key genes in the regulation of this process. These results will provide a good start to dissect the molecular behaviour of plants induced by heat stress and could identify the key genes in stress response for marker-assisted selection in Guar and reveal their roles in stress adaptation in plants.