Transcriptome Analyses Provide Novel Insights into Heat Stress Responses in Chieh-Qua (Benincasa hispida Cogn. var. Chieh-Qua How)

Chromosome-level genome assembly and population genomics reveals crucial selection for subgynoecy development in chieh-qua

Chieh-qua is an important cucurbit crop and very popular in South China and Southeast Asia. Despite its significance, its genetic basis and domestication history are unclear. In this study, we have successfully generated a chromosome-level reference genome assembly for the chieh-qua 'A36' using a hybrid assembly strategy that combines PacBio long reads and Illumina short reads. The assembled genome of chieh-qua is approximately 953.3 Mb in size and is organized into 12 chromosomes, with contig N50 of 6.9 Mb and scaffold N50 of 68.2 Mb. Notably, the chieh-qua genome is comparable in size to the wax gourd genome. Through gene prediction analysis, we have identified a total of 24 593 protein-coding genes in the A36 genome. Additionally, approximately 56.6% (539.3 Mb) of the chieh-qua genome consists of repetitive sequences. Comparative genome analysis revealed that chieh-qua and wax gourd are closely related, indicating a close evolutionary relationship between the two species. Population genomic analysis, employing 129 chieh-qua accessions and 146 wax gourd accessions, demonstrated that chieh-qua exhibits greater genetic diversity compared to wax gourd. We also employed the GWAS method to identify related QTLs associated with subgynoecy, an interested and important trait in chieh-qua. The MYB59 (BhiCQ0880026447) exhibited relatively high expression levels in the shoot apex of four subgynoecious varieties compared with monoecious varieties. Overall, this research provides insights into the domestication history of chieh-qua and offers valuable genomic resources for further molecular research.

Unraveling the genetics of heat tolerance in chickpea landraces (Cicer arietinum L.) using genome-wide association studies

Chickpea, being an important grain legume crop, is often confronted with the adverse effects of high temperatures at the reproductive stage of crop growth, drastically affecting yield and overall productivity. The current study deals with an extensive evaluation of chickpea genotypes, focusing on the traits associated with yield and their response to heat stress. Notably, we observed significant variations for these traits under both normal and high-temperature conditions, forming a robust basis for genetic research and breeding initiatives. Furthermore, the study revealed that yield-related traits exhibited high heritability, suggesting their potential suitability for marker-assisted selection. We carried out single-nucleotide polymorphism (SNP) genotyping using the genotyping-by-sequencing (GBS) method for a genome-wide association study (GWAS). Overall, 27 marker-trait associations (MTAs) linked to yield-related traits, among which we identified five common MTAs displaying pleiotropic effects after applying a stringent Bonferroni-corrected p-value threshold of <0.05 [-log10(p) > 4.95] using the BLINK (Bayesian-information and linkage-disequilibrium iteratively nested keyway) model. Through an in-depth in silico analysis of these markers against the CDC Frontier v1 reference genome, we discovered that the majority of the SNPs were located at or in proximity to gene-coding regions. We further explored candidate genes situated near these MTAs, shedding light on the molecular mechanisms governing heat stress tolerance and yield enhancement in chickpeas such as indole-3-acetic acid-amido synthetase GH3.1 with GH3 auxin-responsive promoter and pentatricopeptide repeat-containing protein, etc. The harvest index (HI) trait was associated with marker Ca3:37444451 encoding aspartic proteinase ortholog sequence of Oryza sativa subsp. japonica and Medicago truncatula, which is known for contributing to heat stress tolerance. These identified MTAs and associated candidate genes may serve as valuable assets for breeding programs dedicated to tailoring chickpea varieties resilient to heat stress and climate change.

Genome-wide identification of bZIP gene family and expression analysis of BhbZIP58 under heat stress in wax gourd

BACKGROUND

The basic leucine zipper (bZIP) transcription factor family is one of the most abundant and evolutionarily conserved gene families in plants. It assumes crucial functions in the life cycle of plants, including pathogen defense, secondary metabolism, stress response, seed maturation, and flower development. Although the genome of wax gourd has been published, little is known about the functions, evolutionary background, and gene expression patterns of the bZIP gene family, which limits its utilization.

RESULTS

A total of 61 bZIP genes (BhbZIPs) were identified from wax gourd (Benincasa hispida) genome and divided into 12 subgroups. Whole-genome duplication (WGD) and dispersed duplication (DSD) were the main driving forces of bZIP gene family expansion in wax gourd, and this family may have undergone intense purifying selection pressure during the evolutionary process. We selected BhbZIP58, only one in the member of subgroup B, to study its expression patterns under different stresses, including heat, salt, drought, cold stress, and ABA treatment. Surprisingly, BhbZIP58 had a dramatic response under heat stress. BhbZIP58 showed the highest expression level in the root compared with leaves, stem, stamen, pistil, and ovary. In addition, BhbZIP58 protein was located in the nucleus and had transcriptional activation activity. Overexpression of BhbZIP58 in Arabidopsis enhanced their heat tolerance.

CONCLUSIONS

In this study, bZIP gene family is systematically bioinformatically in wax gourd for the first time. Particularly, BhbZIP58 may have an important role in heat stress. It will facilitate further research on the bZIP gene family regarding their evolutionary history and biological functions.

Developing future heat-resilient vegetable crops

Climate change seriously impacts global agriculture, with rising temperatures directly affecting the yield. Vegetables are an essential part of daily human consumption and thus have importance among all agricultural crops. The human population is increasing daily, so there is a need for alternative ways which can be helpful in maximizing the harvestable yield of vegetables. The increase in temperature directly affects the plants' biochemical and molecular processes; having a significant impact on quality and yield. Breeding for climate-resilient crops with good yields takes a long time and lots of breeding efforts. However, with the advent of new omics technologies, such as genomics, transcriptomics, proteomics, and metabolomics, the efficiency and efficacy of unearthing information on pathways associated with high-temperature stress resilience has improved in many of the vegetable crops. Besides omics, the use of genomics-assisted breeding and new breeding approaches such as gene editing and speed breeding allow creation of modern vegetable cultivars that are more resilient to high temperatures. Collectively, these approaches will shorten the time to create and release novel vegetable varieties to meet growing demands for productivity and quality. This review discusses the effects of heat stress on vegetables and highlights recent research with a focus on how omics and genome editing can produce temperature-resilient vegetables more efficiently and faster.

Genetic Dissection of Heat Stress Tolerance in Faba Bean (Vicia faba L.) Using GWAS

Heat waves are expected to become more frequent and intense, which will impact faba bean cultivation globally. Conventional breeding methods are effective but take considerable time to achieve breeding goals, and, therefore, the identification of molecular markers associated with key genes controlling heat tolerance can facilitate and accelerate efficient variety development. We phenotyped 134 accessions in six open field experiments during summer seasons at Terbol, Lebanon, at Hudeiba, Sudan, and at Central Ferry, WA, USA from 2015 to 2018. These accessions were genotyped using genotyping by sequencing (GBS), and 10,794 high quality single nucleotide polymorphisms (SNPs) were discovered. These accessions were clustered in one diverse large group, although several discrete groups may exist surrounding it. Fifteen lines belonging to different botanical groups were identified as tolerant to heat. SNPs associated with heat tolerance using single-trait (ST) and multi-trait (MT) genome-wide association studies (GWASs) showed 9 and 11 significant associations, respectively. Through the annotation of the discovered significant SNPs, we found that SNPs from transcription factor helix-loop-helix bHLH143-like S-adenosylmethionine carrier, putative pentatricopeptide repeat-containing protein At5g08310, protein NLP8-like, and photosystem II reaction center PSB28 proteins are associated with heat tolerance.

Comparative transcriptome analysis of heat stress responses of Clematis lanuginosa and Clematis crassifolia

BACKGROUND

Clematis species are attractive ornamental plants with a variety of flower colors and patterns. Heat stress is one of the main factors restricting the growth, development, and ornamental value of Clematis. Clematis lanuginosa and Clematis crassifolia are large-flowered and evergreen Clematis species, respectively, that show different tolerance to heat stress. We compared and analyzed the transcriptome of C. lanuginose and C. crassifolia under heat stress to determine the regulatory mechanism(s) of resistance.

RESULTS

A total of 1720 and 6178 differentially expressed genes were identified from C. lanuginose and C. crassifolia, respectively. The photosynthesis and oxidation-reduction processes of C. crassifolia were more sensitive than C. lanuginose under heat stress. Glycine/serine/threonine metabolism, glyoxylic metabolism, and thiamine metabolism were important pathways in response to heat stress in C. lanuginose, and flavonoid biosynthesis, phenylalanine metabolism, and arginine/proline metabolism were the key pathways in C. crassifolia. Six sHSPs (c176964_g1, c200771_g1, c204924_g1, c199407_g2, c201522_g2, c192936_g1), POD1 (c200317_g1), POD3 (c210145_g2), DREB2 (c182557_g1), and HSFA2 (c206233_g2) may be key genes in the response to heat stress in C. lanuginose and C. crassifolia.

CONCLUSIONS

We compared important metabolic pathways and differentially expressed genes in response to heat stress between C. lanuginose and C. crassifolia. The results increase our understanding of the response mechanism and candidate genes of Clematis under heat stress. These data may contribute to the development of new Clematis varieties with greater heat tolerance.

Unfolding molecular switches in plant heat stress resistance: A comprehensive review

Plant heat stress response is a multi-factorial trait that is precisely regulated by the complex web of transcription factors from various families that modulate heat stress responsive gene expression. Global warming due to climate change affects plant growth and development throughout its life cycle. Adds to this, the frequent occurrence of heat waves is drastically reducing the global crop yield. Molecular plant scientists can help crop breeders by providing genetic markers associated with stress resistance. Plant heat stress response (HSR), however, is a multi-factorial trait and using a single stress resistance trait might not be ideal to develop thermotolerant crops. Transcription factors participate in regulation of plant biological processes and environmental stress responses. Recent studies have revealed that plant HSR is precisely regulated by the complex web of transcription factors from various families. These transcription factors enhance plant heat stress tolerance by regulating the expression level of several stress-responsive genes independently or in cross talk with different other transcription factors. This review explores how signaling pathways triggered by heat stress are regulated by multiple transcription factor families. To our knowledge, we for the first time analyze the role of major transcription factor families in plant HSR along with their regulatory mechanisms. In the end, we will also discuss the potential of emerging technologies to improve thermotolerance in plants.

Expression Profiling of Heat Shock Protein Genes as Putative Early Heat-Responsive Members in Lettuce

High temperatures due to global warming can cause harmful effects on the productivity of lettuce, a cool-season crop. To identify lettuce heat shock protein (HSP) genes that could be involved in early responses to heat stress in plants, we compared RNA transcriptomes between lettuce plants with and without heat treatment of 37 °C for 1 h. Using transcriptome sequencing analyses, a total of 7986 differentially expressed genes (DEGs) were identified including the top five, LsHSP70A, LsHSP70B, LsHSP17.3A, LsHSP17.9A and LsHSP17.9B, which were the most highly differentially expressed genes. In order to investigate the temporal expression patterns of 24 lettuce HSP genes with a fold-change greater than 100 under heat stress, the expression levels of the genes were measured by qRT-PCR at 0, 1, 4, 8, 14, and 24 h time points after heat treatment. The 24 LsHSP genes were classified into three groups based on the phylogenetic analysis and/or major domains available in each protein, and we provided a potential link between the phylogenetic relationships and expression patterns of the LsHSP genes. Our results showed putative early heat-responsive lettuce HSP genes that could be possible candidates as breeding guides for the development of heat-tolerant lettuce cultivars.

Differential expression gene/protein contribute to heat stress-responsive in Tetraena propinqua in Saudi Arabia

Within their natural habitat, plants are subjected to abiotic stresses that include heat stress. In the current study, the effect of 4 h, 24 h, and 48 h of heat stress on Tetraena propinqua ssp. migahidii seedling's protein profile and proteomic analyses were investigated. Total soluble protein SDS-PAGE (Sodium dodecyl sulfate-polyacrylamide gel electrophoresis) profile showed 18-protein bands, the newly synthesized protein band (with molecular weights 86.5, 30.2 and 31.4 KD) at 24 h of heat stress and 48 of normal conditions. Proteomic analysis showed that 81 and 930 targets are involved in gene and protein expression respectively. At 4 h, 57 genes and 110 proteins in C4 reached 56 genes and 173 proteins in T4. At 24 h, 63 genes and 180 proteins in C24 decreased to 54 genes and 151 protein in T24. After 48 h, 56 genes and 136 proteins in C48 increased to 64 genes and 180 proteins in T48. The genes and proteins involved in transcription, translation, photosynthesis, transport, and other unknown metabolic processes, were differentially expressed under treatments of heat stress. These findings provide insights into the molecular mechanisms related to heat stress, in addition to its influence on the physiological traits of T. propinqua seedlings. Heat stress-mediated differential regulation genes indicate a role in the development and stress response of T. propinqua. The candidate dual-specificity genes and proteins identified in this study paves way for more molecular analysis of up-and-down-regulation.

Quantitative proteomic analysis to capture the role of heat-accumulated proteins in moss plant acquired thermotolerance

At dawn of a scorching summer day, land plants must anticipate upcoming extreme midday temperatures by timely establishing molecular defences that can keep heat-labile membranes and proteins functional. A gradual morning pre-exposure to increasing sub-damaging temperatures induces heat-shock proteins (HSPs) that are central to the onset of plant acquired thermotolerance (AT). To gain knowledge on the mechanisms of AT in the model land plant Physcomitrium patens, we used label-free LC-MS/MS proteomics to quantify the accumulated and depleted proteins before and following a mild heat-priming treatment. High protein crowding is thought to promote protein aggregation, whereas molecular chaperones prevent and actively revert aggregation. Yet, we found that heat priming (HP) did not accumulate HSP chaperones in chloroplasts, although protein crowding was six times higher than in the cytosol. In contrast, several HSP20s strongly accumulated in the cytosol, yet contributing merely 4% of the net mass increase of heat-accumulated proteins. This is in poor concordance with their presumed role at preventing the aggregation of heat-labile proteins. The data suggests that under mild HP unlikely to affect protein stability. Accumulating HSP20s leading to AT, regulate the activity of rare and specific signalling proteins, thereby preventing cell death under noxious heat stress.

Comparative transcriptome analysis uncovers different heat stress responses in heat-resistant and heat-sensitive jujube cultivars

Jujube (Ziziphus jujuba Mill.) is an economically and agriculturally significant fruit crop and is widely cultivated throughout the world. Heat stress has recently become a primary abiotic stressor limiting the productivity and growth of jujube, as well as other crops. There are few studies, however, that have performed transcriptome profiling of jujube when it is exposed to heat stress. In this study, we observed the physiochemical changes and analyzed gene expression profiles in resistant jujube cultivar ‘HR’ and sensitive cultivar ‘HS’ subjected to heat stress for 0, 1, 3, and 5d. Twenty-four cDNA libraries from ‘HR’ and ‘HS’ leaves were built with a transcriptome assay. A total of 6887 and 5077 differentially expressed genes were identified in ‘HR’ and ‘HS’ after 1d, 3d, and 5d of heat stress compared with the control treatment, GO and KEGG enrichment analysis revealed that some of the genes were highly enriched in oxidation-reduction process, response to stress, response to water deprivation, response to heat, carbon metabolism, protein processing in endoplasmic reticulum, and plant hormone signal transduction and may play vital roles in the heat stress response in jujube plants. Differentially expressed genes were identified in the two cultivars, including heat shock proteins, transcriptional factors, and ubiquitin-protein ligase genes. And the expression pattern of nine genes was also validated by qRT-PCR. These results will provide useful information for elucidating the molecular mechanism underlying heat stress in different jujube cultivars.

Genome-wide analysis of H3K4me3 and H3K27me3 modifications due to Lr28 for leaf rust resistance in bread wheat (Triticum aestivum)

Present study revealed a complex relationship among histone H3 methylation (examined using H3K4/K27me3 marks), cytosine DNA methylation and differential gene expression during Lr28 mediated leaf rust resistance in wheat. During the present study, genome-wide histone modifications were examined in a pair of near isogenic lines (NILs) (with and without Lr28 in the background of cv. HD2329). The two histone marks used included H3K4me3 (an activation mark) and H3K27me3 (a repression mark). The results were compared with levels of expression (using RNA-seq) and DNA methylation (MeDIP) data obtained using the same pair of NILs. Some of the salient features of the present study include the following: (i) large scale differential binding sites (DBS) were available for only H3K4me3 in the susceptible cultivar, but for both H3K4me3 and H3K27me3 in its resistant NIL; (ii) DBSs for H3K27me3 mark were more abundant (> 80%) in intergenic regions, whereas DBSs for H3K4me3 were distributed in all genomic regions including exons, introns, intergenic, TTS (transcription termination sites) and promoters; (iii) fourteen (14) genes associated with DBSs showed co-localization for both the marks; (iv) only a small fraction (7% for H3K4me3 and 12% for H3K27me3) of genes associated with DBSs matched with the levels of gene expression inferred from RNA-seq data; (v) validation studies using qRT-PCR were conducted on 26 selected representative genes; results for only 11 genes could be validated. The proteins encoded by important genes involved in promoting infection included domains generally carried by R gene proteins such as Mlo like protein, protein kinases and purple acid phosphatase. Similarly, proteins encoded by genes involved in resistance included those carrying domains for lectin kinase, R gene, aspartyl protease, etc. Overall, the results suggest a very complex network of downstream genes that are expressed during compatible and incompatible interactions; some of the genes identified during the present study may be used in future validation studies involving RNAi/overexpression approaches.

Metabolome and transcriptome analyses reveal chlorophyll and anthocyanin metabolism pathway associated with cucumber fruit skin color

BACKGROUND

Fruit skin color play important role in commercial value of cucumber, which is mainly determined by the content and composition of chlorophyll and anthocyanins. Therefore, understanding the related genes and metabolomics involved in composition of fruit skin color is essential for cucumber quality and commodity value.

RESULTS

The results showed that chlorophyll a, chlorophyll b and carotenoid content in fruit skin were higher in Lv (dark green skin) than Bai (light green skin) on fruit skin. Cytological observation showed more chloroplast existed in fruit skin cells of Lv. A total of 162 significantly different metabolites were found between the fruit skin of the two genotypes by metabolome analysis, including 40 flavones, 9 flavanones, 8 flavonols, 6 anthocyanins, and other compounds. Crucial anthocyanins and flavonols for fruit skin color, were detected significantly decreased in fruit skin of Bai compared with Lv. By RNA-seq assay, 4516 differentially expressed genes (DEGs) were identified between two cultivars. Further analyses suggested that low expression level of chlorophyll biosynthetic genes, such as chlM, por and NOL caused less chlorophylls or chloroplast in fruit skin of Bai. Meanwhile, a predicted regulatory network of anthocyanin biosynthesis was established to illustrate involving many DEGs, especially 4CL, CHS and UFGT.

CONCLUSIONS

This study uncovered significant differences between two cucumber genotypes with different fruit color using metabolome and RNA-seq analysis. We lay a foundation to understand molecular regulation mechanism on formation of cucumber skin color, by exploring valuable genes, which is helpful for cucumber breeding and improvement on fruit skin color.

Transcriptome profiling and gene expression analyses of eggplant (Solanum melongena L.) under heat stress

Global warming induces heat stress in eggplant, seriously affecting its quality and yield. The response to heat stress is a complex regulatory process; however, the exact mechanism in eggplant is unknown. We analyzed the transcriptome of eggplant under different high-temperature treatments using RNA-Seq technology. Three libraries treated at high temperatures were generated and sequenced. There were 40,733,667, 40,833,852, and 40,301,285 clean reads with 83.98%, 79.69%, and 84.42% of sequences mapped to the eggplant reference genome in groups exposed to 28°C (CK), 38°C (T38), and 43°C (T43), respectively. There were 3,067 and 1,456 DEGs in T38 vs CK and T43 vs CK groups, respectively. In these two DEG groups, 315 and 342 genes were up- and down-regulated, respectively, in common. Differential expression patterns of DEGs in antioxidant enzyme systems, detoxication, phytohormones, and transcription factors under heat stress were investigated. We screened heat stress-related genes for further validation by qRT-PCR. Regulation mechanisms may differ under different temperature treatments, in which heat shock proteins and heat stress transcription factors play vital roles. These results provide insight into the molecular mechanisms of the heat stress response in eggplant and may be useful in crop breeding.

Understanding the heat resistance of cucumber through leaf transcriptomics

Heat stress is a major environmental factor limiting plant productivity and quality in agriculture. Cucumber, one of the most important vegetables among cucurbitaceae, prefers to grow in a warm environment. Until now the molecular knowledge of heat stress in cucumber remained unclear. In this study, we performed transcriptome analysis using two diverse genetic cucumber cultivars, L-9 and A-16 grown under normal and heat stress. L-9 displayed heat-tolerance phenotype with higher superoxide dismutase enzyme (SOD) enzyme activity and lower malondialdehyde (MDA) content than A-16 under heat stress. RNA-sequencing revealed that a total of 963 and 2778 genes are differentially expressed between L-9 and A-16 under normal and heat stress respectively. In addition, we found that differentially expressed genes (DEGs) associated with plant hormones signally pathway, transcription factors, and secondary metabolites showed significantly change in expression level after heat stress, which were confirmed by quantitative real-time PCR assay. Our results not only explored several crucial genes involved in cucumber heat resistance, but also provide a new insight into studying heat stress.