Linking Brassinosteroid and ABA Signaling in the Context of Stress Acclimation

Transcriptome Profiles Reveal Key Regulatory Networks during Single and Multifactorial Stresses Coupled with Melatonin Treatment in Pitaya (Selenicereus undatus L.)

In response to evolving climatic conditions, plants frequently confront multiple abiotic stresses, necessitating robust adaptive mechanisms. This study focuses on the responses of Selenicereus undatus L. to both individual stresses (cadmium; Cd, salt; S, and drought; D) and their combined applications, with an emphasis on evaluating the mitigating effects of (M) melatonin. Through transcriptome analysis, this study identifies significant gene expression changes and regulatory network activations. The results show that stress decreases pitaya growth rates by 30%, reduces stem and cladode development by 40%, and increases Cd uptake under single and combined stresses by 50% and 70%, respectively. Under stress conditions, enhanced activities of H2O2, POD, CAT, APX, and SOD and elevated proline content indicate strong antioxidant defenses. We identified 141 common DEGs related to stress tolerance, most of which were related to AtCBP, ALA, and CBP pathways. Interestingly, the production of genes related to signal transduction and hormones, including abscisic acid and auxin, was also significantly induced. Several calcium-dependent protein kinase genes were regulated during M and stress treatments. Functional enrichment analysis showed that most of the DEGs were enriched during metabolism, MAPK signaling, and photosynthesis. In addition, weighted gene co-expression network analysis (WGCNA) identified critical transcription factors (WRKYs, MYBs, bZIPs, bHLHs, and NACs) associated with antioxidant activities, particularly within the salmon module. This study provides morpho-physiological and transcriptome insights into pitaya's stress responses and suggests molecular breeding techniques with which to enhance plant resistance.

Transcriptomic analyses provide molecular insight into the cold stress response of cold-tolerant alfalfa

BACKGROUND

Daye No.3 is a novel cultivar of alfalfa (Medicago sativa L.) that is well suited for cultivation in high-altitude regions such as the Qinghai‒Tibet Plateau owing to its high yield and notable cold resistance. However, the limited availability of transcriptomic information has hindered our investigation into the potential mechanisms of cold tolerance in this cultivar. Consequently, we conducted de novo transcriptome assembly to overcome this limitation. Subsequently, we compared the patterns of gene expression in Daye No. 3 during cold acclimatization and exposure to cold stress at various time points.

RESULTS

A total of 15 alfalfa samples were included in the transcriptome assembly, resulting in 141.97 Gb of clean bases. A total of 441 DEGs were induced by cold acclimation, while 4525, 5016, and 8056 DEGs were identified at 12 h, 24 h, and 36 h after prolonged cold stress at 4 °C, respectively. The consistency between the RT‒qPCR and transcriptome data confirmed the accuracy and reliability of the transcriptomic data. KEGG enrichment analysis revealed that many genes related to photosynthesis were enriched under cold stress. STEM analysis demonstrated that genes involved in nitrogen metabolism and the TCA cycle were consistently upregulated under cold stress, while genes associated with photosynthesis, particularly antenna protein genes, were downregulated. PPI network analysis revealed that ubiquitination-related ribosomal proteins act as hub genes in response to cold stress. Additionally, the plant hormone signaling pathway was activated under cold stress, suggesting its vital role in the cold stress response of alfalfa.

CONCLUSIONS

Ubiquitination-related ribosomal proteins induced by cold acclimation play a crucial role in early cold signal transduction. As hub genes, these ubiquitination-related ribosomal proteins regulate a multitude of downstream genes in response to cold stress. The upregulation of genes related to nitrogen metabolism and the TCA cycle and the activation of the plant hormone signaling pathway contribute to the enhanced cold tolerance of alfalfa.

Chromatin modifications and memory in regulation of stress-related polyphenols: finding new ways to control flavonoid biosynthesis

The influence of epigenetic factors on plant defense responses and the balance between growth and defense is becoming a central area in plant biology. It is believed that the biosynthesis of secondary metabolites can be regulated by epigenetic factors, but this is not associated with the formation of a "memory" to the previous biosynthetic status. This review shows that some epigenetic effects can result in epigenetic memory, which opens up new areas of research in secondary metabolites, in particular flavonoids. Plant-controlled chromatin modifications can lead to the generation of stress memory, a phenomenon through which information regarding past stress cues is retained, resulting in a modified response to recurring stress. How deeply are the mechanisms of chromatin modification and memory generation involved in the control of flavonoid biosynthesis? This article collects available information from the literature and interactome databases to address this issue. Visualization of the interaction of chromatin-modifying proteins with the flavonoid biosynthetic machinery is presented. Chromatin modifiers and "bookmarks" that may be involved in the regulation of flavonoid biosynthesis through memory have been identified. Through different mechanisms of chromatin modification, plants can harmonize flavonoid metabolism with: stress responses, developmental programs, light-dependent processes, flowering, and longevity programs. The available information points to the possibility of developing chromatin-modifying technologies to control flavonoid biosynthesis.

Genome-Wide Identification, Characterization, and Expression Analysis of the BES1 Family Genes under Abiotic Stresses in Phoebe bournei

The BRI1 EMS suppressor 1(BES1) transcription factor is a crucial regulator in the signaling pathway of Brassinosteroid (BR) and plays an important role in plant growth and response to abiotic stress. Although the identification and functional validation of BES1 genes have been extensively explored in various plant species, the understanding of their role in woody plants-particularly the endangered species Phoebe bournei (Hemsl.) Yang-remains limited. In this study, we identified nine members of the BES1 gene family in the genome of P. bournei; these nine members were unevenly distributed across four chromosomes. In our further evolutionary analysis of PbBES1, we discovered that PbBES1 can be divided into three subfamilies (Class I, Class II, and Class IV) based on the evolutionary tree constructed with Arabidopsis thaliana, Oryza sativa, and Solanum lycopersicum. Each subfamily contains 2-5 PbBES1 genes. There were nine pairs of homologous BES1 genes in the synteny analysis of PbBES1 and AtBES1. Three segmental replication events and one pair of tandem duplication events were present among the PbBES1 family members. Additionally, we conducted promoter cis-acting element analysis and discovered that PbBES1 contains binding sites for plant growth and development, cell cycle regulation, and response to abiotic stress. PbBES1.2 is highly expressed in root bark, stem bark, root xylem, and stem xylem. PbBES1.3 was expressed in five tissues. Moreover, we examined the expression profiles of five representative PbBES1 genes under heat and drought stress. These experiments preliminarily verified their responsiveness and functional roles in mediating responses to abiotic stress. This study provides important clues to elucidate the functional characteristics of the BES1 gene family, and at the same time provides new insights and valuable information for the regulation of resistance in P. bournei.

Whole-genome sequencing of tetraploid potato varieties reveals different strategies for drought tolerance

Climate changes leading to increasingly longer seasonal drought periods in large parts of the world increase the necessity for breeding drought-tolerant crops. Cultivated potato (Solanum tuberosum), the third most important vegetable crop worldwide, is regarded as drought-sensitive due to its shallow root architecture. Two German tetraploid potato cultivars differing in drought tolerance and their F1-progeny were evaluated under various drought scenarios. Bulked segregant analyses were combined with whole-genome sequencing (BSA-Seq) using contrasting bulks of drought-tolerant and drought-sensitive F1-clones. Applying QTLseqr, 15 QTLs comprising 588,983 single nucleotide polymorphisms (SNPs) in 2325 genes associated with drought stress tolerance were identified. SeqSNP analyses in an association panel of 34 mostly starch potato varieties using 1-8 SNPs for each of 188 selected genes narrowed the number of candidate genes down to 10. In addition, ent-kaurene synthase B was the only gene present under QTL 10. Eight of the identified genes (StABP1, StBRI1, StKS, StLEA, StPKSP1, StPKSP2, StYAB5, and StZOG1) address plant development, the other three genes (StFATA, StHGD and StSYP) contribute to plant protection under drought stress. Allelic variation in these genes might be explored in future breeding for drought-tolerant potato varieties.

The ORFans' tale: new insights in plant biology

Orphan genes (OGs) are protein-coding genes without a significant sequence similarity in closely related species. Despite their functional importance, very little is known about the underlying molecular mechanisms by which OGs participate in diverse biological processes. Here, we discuss the evolutionary mechanisms of OGs' emergence with relevance to species-specific adaptations. We also provide a mechanistic view of the involvement of OGs in multiple processes, including growth, development, reproduction, and carbon-metabolism-mediated immunity. We highlight the interconnection between OGs and the sucrose nonfermenting 1 (SNF1)-related protein kinases (SnRKs)-target of rapamycin (TOR) signaling axis for phytohormone signaling, nutrient metabolism, and stress responses. Finally, we propose a high-throughput pipeline for OGs' interspecies and intraspecies gene transfer through a transgenic approach for future biotechnological advances.

StBIN2 Positively Regulates Potato Formation through Hormone and Sugar Signaling

Potato is an important food crop worldwide. Brassinosteroids (BRs) are widely involved in plant growth and development, and BIN2 (brassinosteroid insensitive 2) is the negative regulator of their signal transduction. However, the function of BIN2 in the formation of potato tubers remains unclear. In this study, transgenic methods were used to regulate the expression level of StBIN2 in plants, and tuber related phenotypes were analyzed. The overexpression of StBIN2 significantly increased the number of potatoes formed per plant and the weight of potatoes in transgenic plants. In order to further explore the effect of StBIN2 on the formation of potato tubers, this study analyzed BRs, ABA hormone signal transduction, sucrose starch synthase activity, the expression levels of related genes, and interacting proteins. The results show that the overexpression of StBIN2 enhanced the downstream transmission of ABA signals. At the same time, the enzyme activity of the sugar transporter and the expression of synthetic genes were increased in potato plants overexpressing StBIN2, which also demonstrated the upregulation of sucrose and the expression of the starch synthesis gene. Apparently, StBIN2 affected the conversion and utilization of key substances such as glucose, sucrose, and starch in the process of potato formation so as to provide a material basis and energy preparation for forming potatoes. In addition, StBIN2 also promoted the expression of the tuber formation factors StSP6A and StS6K. Altogether, this investigation enriches the study on the mechanism through which StBIN2 regulates potato tuber formation and provides a theoretical basis for achieving a high and stable yield of potato.

A novel TF molecular switch-mechanism found in two contrasting ecotypes of a psammophyte, Agriophyllum squarrosum, in regulating transcriptional drought memory

BACKGROUND

Prior drought stress may change plants response patterns and subsequently increase their tolerance to the same condition, which can be referred to as "drought memory" and proved essential for plants well-being. However, the mechanism of transcriptional drought memory in psammophytes remains unclear. Agriophyllum squarrosum, a pioneer species on mobile dunes, is widely spread in Northern China's vast desert areas with outstanding ability of water use efficiency. Here we conducted dehydration-rehydration treatment on A. squarrosum semi-arid land ecotype AEX and arid land ecotype WW to dissect the drought memory mechanism of A. squarrosum, and to determine the discrepancy in drought memory of two contrasting ecotypes that had long adapted to water heterogeneity.

RESULT

Physiological traits monitoring unveiled the stronger ability and longer duration in drought memory of WW than that of AEX. A total of 1,642 and 1,339 drought memory genes (DMGs) were identified in ecotype AEX and WW, respectively. Furthermore, shared DMGs among A. squarrosum and the previously studied species depicted that drought memory commonalities in higher plants embraced pathways like primary and secondary metabolisms; while drought memory characteristics in A. squarrosum were mainly related to response to heat, high light intensity, hydrogen peroxide, and dehydration, which might be due to local adaptation to desert circumstances. Heat shock proteins (HSPs) occupied the center of the protein-protein interaction (PPI) network in drought memory transcription factors (TF), thus playing a key regulatory role in A. squarrosum drought memory. Co-expression analysis of drought memory TFs and DMGs uncovered a novel regulating module, whereby pairs of TFs might function as molecular switches in regulating DMG transforming between high and low expression levels, thus promoting drought memory reset.

CONCLUSION

Based on the co-expression analysis, protein-protein interaction prediction, and drought memory metabolic network construction, a novel regulatory module of transcriptional drought memory in A. squarrosum was hypothesized here, whereby recurrent drought signal is activated by primary TF switches, then amplified by secondary amplifiers, and thus regulates downstream complicated metabolic networks. The present research provided valuable molecular resources on plants' stress-resistance basis and shed light on drought memory in A. squarrosum.

Phytohormones regulate the non-redundant response of ω-3 fatty acid desaturases to low temperatures in Chorispora bungeana

To explore the contributions of ω-3 fatty acid desaturases (FADs) to cold stress response in a special cryophyte, Chorispora bungeana, two plastidial ω-3 desaturase genes (CbFAD7, CbFAD8) were cloned and verified in an Arabidopsis fad7fad8 mutant, before being compared with the microsomal ω-3 desaturase gene (CbFAD3). Though these genes were expressed in all tested tissues of C. bungeana, CbFAD7 and CbFAD8 have the highest expression in leaves, while CbFAD3 was mostly expressed in suspension-cultured cells. Low temperatures resulted in significant increases in trienoic fatty acids (TAs), corresponding to the cooperation of CbFAD3 and CbFAD8 in cultured cells, and the coordination of CbFAD7 and CbFAD8 in leaves. Furthermore, the cold induction of CbFAD8 in the two systems were increased with decreasing temperature and independently contributed to TAs accumulation at subfreezing temperature. A series of experiments revealed that jasmonie acid and brassinosteroids participated in the cold-responsive expression of ω-3 CbFAD genes in both C. bungeana cells and leaves, while the phytohormone regulation in leaves was complex with the participation of abscisic acid and gibberellin. These results point to the hormone-regulated non-redundant contributions of ω-3 CbFADs to maintain appropriate level of TAs under low temperatures, which help C. bungeana survive in cold environments.

Advances in Receptor-like Protein Kinases in Balancing Plant Growth and Stress Responses

Accompanying the process of growth and development, plants are exposed to ever-changing environments, which consequently trigger abiotic or biotic stress responses. The large protein family known as receptor-like protein kinases (RLKs) is involved in the regulation of plant growth and development, as well as in the response to various stresses. Understanding the biological function and molecular mechanism of RLKs is helpful for crop breeding. Research on the role and mechanism of RLKs has recently received considerable attention regarding the balance between plant growth and environmental adaptability. In this paper, we systematically review the classification of RLKs, the regulatory roles of RLKs in plant development (meristem activity, leaf morphology and reproduction) and in stress responses (disease resistance and environmental adaptation). This review focuses on recent findings revealing that RLKs simultaneously regulate plant growth and stress adaptation, which may pave the way for the better understanding of their function in crop improvement. Although the exact crosstalk between growth constraint and plant adaptation remains elusive, a profound study on the adaptive mechanisms for decoupling the developmental processes would be a promising direction for the future research.

Genome-wide identification, characterization and gene expression of BES1 transcription factor family in grapevine (Vitis vinifera L.)

BES1, as the most important transcription factor responsible for brassinolide (BR) signaling, has been confirmed to play a significant role in regulating plant growth and the improvement of stress resistance. The transcriptional regulatory mechanism of BES1 has been well elucidated in several plants, such as Arabidopsis thaliana (A. thaliana), Triticum aestivum L. (T. aestivum), and Oryza sativa L. (O. sativa). Nevertheless, the genome-wide analysis of the BES1 family in Vitis vinifera L. (V. vinifera). has not been comprehensively carried out. Thus, we have conducted a detailed analysis and identification of the BES1 transcription factors family in V. vinifera; a total of eight VvBES1 genes was predicted, and the phylogenetic relationships, gene structures, and Cis-acting element in their promoters were also analyzed. BES1 genes have been divided into three groups (I, II and III) based on phylogenetic relationship analysis, and most of VvBES1 genes were in group III. Also, we found that VvBES1 genes was located at seven of the total nineteen chromosomes, whereas VvBES1-2 (Vitvi04g01234) and VvBES1-5 (Vitvi18g00924) had a collinearity relationship, and their three copies are well preserved. In addition, the intron-exon model of VvBES1 genes were mostly conserved, and there existed several Cis-acting elements related to stress resistance responsive and phytohormones responsive in BES1s genes promoter. Moreover, the BES1 expressions were different in different V. vinifera organs, and BES1 expressions were different in different V. vinifera varieties under saline-alkali stress and heat stress, the expression of VvBES1 also changed with the prolongation of saline-alkali stress treatment time. The above findings could not only lay a primary foundation for the further validation of VvBES1 function, but could also provide a reference for molecular breeding in V. vinifera.

Physiological and Transcriptomic Analysis provide Molecular Insight into 24-Epibrassinolide mediated Cr(VI)-Toxicity Tolerance in Pepper Plants

The ever-increasing industrial activities over the decades have generated high toxic metals such as chromium (Cr) that hampers plant growth and development. To counter Cr-toxicity, plants have evolved complex defensive systems including hormonal crosstalk with various signaling pathways. 24-epibrassinolide (24-EBR) lowers oxidative stress and alleviates Cr(VI)-toxicity in plants. In this study, the concealed BR-mediated influences on Cr(VI)-stress tolerance were explored by transcriptome analysis in the Capsicum annuum. Results revealed a linkage between plant development under Cr(VI)-stress and the mitigating effect of 24-epibrassinolide and brassinazole. Growth inhibition, chlorophyll degradation, and a significant rise of malondialdehyde (MDA) were observed after 40 mg/L Cr(VI) treatment in Brz supplemented seedlings, whereas 24-EBR supplemented seedlings exhibited commendatory effect. Comparative transcriptome analysis showed that the expression levels of 6687 genes changed (3846 up-regulated and 2841 downregulated) under Cr(VI)-stress with Brz supplementation. Whereas the expression levels of only 1872 genes changed under Cr(VI)-stress with 24-EBR supplementation (1223 up-regulated and 649 downregulated). The functional categories of the differentially expressed genes (DEGs) by gene ontology (GO) revealed that drug transport, defense responses, and drug catabolic process were the considerable enrichments between 24-EBR and Brz supplemented seedlings under Cr(VI)-stress. Furthermore, auxin signaling, glutathione metabolism, ABC transporters, MAPK pathway, and 36 heavy metal-related genes were significantly differentially expressed components between Cr(VI)-stress, 24-EBR, and Brz supplemented seedlings. Overall, our data demonstrate that employing 24-EBR can commendably act as a growth stimulant in plants subjected to Cr(VI)-stress by modulating the physiological and defense regulatory system.

Brassinosteroid Accelerates Wound Healing of Potato Tubers by Activation of Reactive Oxygen Metabolism and Phenylpropanoid Metabolism

Wound healing could effectively reduce the decay rate of potato tubers after harvest, but it took a long time to form typical and complete healing structures. Brassinosteroid (BR), as a sterol hormone, is important for enhancing plant resistance to abiotic and biotic stresses. However, it has not been reported that if BR affects wound healing of potato tubers. In the present study, we observed that BR played a positive role in the accumulation of lignin and suberin polyphenolic (SPP) at the wounds, and effectively reduced the weight loss and disease index of potato tubers (cv. Atlantic) during healing. At the end of healing, the weight loss and disease index of BR group was 30.8% and 23.1% lower than the control, respectively. Furthermore, BR activated the expression of StPAL, St4CL, StCAD genes and related enzyme activities in phenylpropanoid metabolism, and promoted the synthesis of lignin precursors and phenolic acids at the wound site, mainly by inducing the synthesis of caffeic acid, sinapic acid and cinnamyl alcohol. Meanwhile, the expression of StNOX was induced and the production of O2- and H2O2 was promoted, which mediated oxidative crosslinking of above phenolic acids and lignin precursors to form SPP and lignin. In addition, the expression level of StPOD was partially increased. In contrast, the inhibitor brassinazole inhibited phenylpropanoid metabolism and reactive oxygen metabolism, and demonstrated the function of BR hormone in healing in reverse. Taken together, the activation of reactive oxygen metabolism and phenylpropanoid metabolism by BR could accelerate the wound healing of potato tubers.

Genome-wide survey of the dehydrin genes in bread wheat (Triticum aestivum L.) and its relatives: identification, evolution and expression profiling under various abiotic stresses

BACKGROUND

Bread wheat (Triticum aestivum) is an important staple cereal grain worldwide. The ever-increasing environmental stress makes it very important to mine stress-resistant genes for wheat breeding programs. Therefore, dehydrin (DHN) genes can be considered primary candidates for such programs, since they respond to multiple stressors.

RESULTS

In this study, we performed a genome-wide analysis of the DHN gene family in the genomes of wheat and its three relatives. We found 55 DHN genes in T. aestivum, 31 in T. dicoccoides, 15 in T. urartu, and 16 in Aegilops tauschii. The phylogenetic, synteny, and sequence analyses showed we can divide the DHN genes into five groups. Genes in the same group shared similar conserved motifs and potential function. The tandem TaDHN genes responded strongly to drought, cold, and high salinity stresses, while the non-tandem genes respond poorly to all stress conditions. According to the interaction network analysis, the cooperation of multiple DHN proteins was vital for plants in combating abiotic stress.

CONCLUSIONS

Conserved, duplicated DHN genes may be important for wheat being adaptable to a different stress conditions, thus contributing to its worldwide distribution as a staple food. This study not only highlights the role of DHN genes help the Triticeae species against abiotic stresses, but also provides vital information for the future functional studies in these crops.

Control of ABA Signaling and Crosstalk with Other Hormones by the Selective Degradation of Pathway Components

A rapid and appropriate genetic and metabolic acclimation, which is crucial for plants' survival in a changing environment, is maintained due to the coordinated action of plant hormones and cellular degradation mechanisms influencing proteostasis. The plant hormone abscisic acid (ABA) rapidly accumulates in plants in response to environmental stress and plays a pivotal role in the reaction to various stimuli. Increasing evidence demonstrates a significant role of autophagy in controlling ABA signaling. This field has been extensively investigated and new discoveries are constantly being provided. We present updated information on the components of the ABA signaling pathway, particularly on transcription factors modified by different E3 ligases. Then, we focus on the role of selective autophagy in ABA pathway control and review novel evidence on the involvement of autophagy in different parts of the ABA signaling pathway that are important for crosstalk with other hormones, particularly cytokinins and brassinosteroids.

Integration of mRNA and miRNA analysis reveals the molecular mechanism of potato (Solanum tuberosum L.) response to alkali stress

The continuing increase in the global saline-alkali land area has made saline-alkali stress the principal abiotic stress limiting plant growth. Potato is the most important non-grain crop, and its production is also severely limited by saline-alkali stress. However, few studies have addressed the mechanism of saline-alkali tolerance of potato with a focus on its response to neutral salt NaCl stress, or its response to alkali stress. Recently, miRNA-mRNA analyses have helped advance our understanding of how plants respond to stress. Here, we have characterized the morphological, physiological, and transcriptome changes of tissue culture seedlings of potato variety "Qingshu No. 9" treated with NaHCO₃ (for 0, 2, 6, and 24 h). We found that the leaves of tissue culture seedlings wilted and withered under alkali stress, and the contents of ABA, BRs, trehalose, and lignin in roots increased significantly. The contents of GAs decreased significantly. Subsequently, miRNA-seq analysis results identified 168 differentially expressed miRNAs (DEMIs) under alkali stress, including 21 exist miRNAs and 37 known miRNAs from 47 families and 110 novel miRNAs. The mRNA-seq results identified 5731 differentially expressed mRNAs (DEMs) under alkali stress. By miRNA-mRNA integrated analysis, were obtained 33 miRNA-target gene pairs composed of 20 DEMIs and 33 DEMs. Next, we identified the "phenylpropanoid biosynthesis", "plant hormone signal transduction", and "starch and sucrose metabolism" pathways as necessary for potato to respond to alkali stress. miR4243-x and novel-m064-5p were involved in the response of potato to alkali stress by their negative regulatory effects on shikimate O-hydroxycinnamoyltransferase (HCT) and sucrose-phosphate synthase (SPS) genes, respectively. The expression results of miRNA and mRNA were verified by quantitative real-time PCR (qRT-PCR). Our results clarify the mechanism of potato response to alkali stress at the miRNA level, providing new insights into the molecular mechanisms of potato's response to alkali stress. We report many candidate miRNAs and mRNAs for molecular-assisted screening and salt-alkali resistance breeding.

Basic Protein Modules Combining Abscisic Acid and Light Signaling in Arabidopsis

It is generally accepted that plants use the complex signaling system regulated by light and abscisic acid (ABA) signaling components to optimize growth and development in different situations. The role of ABA-light interactions is evident in the coupling of stress defense reactions with seed germination and root development, maintaining of stem cell identity and stem cell specification, stem elongation and leaf development, flowering and fruit formation, senescence, and shade avoidance. All these processes are regulated jointly by the ABA-light signaling system. Although a lot of work has been devoted to ABA-light signal interactions, there is still no systematic description of central signaling components and protein modules, which jointly regulate plant development. New data have emerged to promote understanding of how ABA and light signals are integrated at the molecular level, representing an extensively growing area of research. This work is intended to fill existing gaps by using literature data combined with bioinformatics analysis.