under drought stress

AP2 (APETALA2)/EREBP (ethylene-responsive element-binding protein), cytochrome c oxidase (COX) and nonspecific lipid transfer proteins (LTP) play important roles in the response to drought stress. This is the first study to identify the COX gene in Zea mays L. via genome-wide analysis. The qRT‒PCR results indicated that AP2/EREBP, COX and LTP were downregulated, with fold changes of 0.84, 0.53 and 0.31, respectively, after 12 h of drought stress. Genome-wide analysis identified 78 AP2/EREBP, 6 COX and 10 LTP genes in Z. mays L. Domain analysis confirmed the presence of the AP2 domain, Cyt_c_Oxidase_Vb domain and nsLTP1 in the AP2/EREBP, COX and LTP proteins, respectively. The AP2/EREBP protein family (AP2) includes five different domain types: the AP2/ERF domain, the EREBP-like factor (EREBP), the ethylene responsive factor (ERF), the dehydration responsive element binding protein (DREB) and the SHN SHINE. Synteny analysis of the AP2/EREBP, COX and LTP genes revealed collinearity orthologous relationships in O. sativa, H. vulgare and A. thaliana. AP2/EREBP genes were found on the 10 chromosomes of Z. mays L. COX genes were found on chromosomes 1, 3, 4, 5, 7 and 8. LTP genes were found on chromosomes 1, 3, 6, 8, 9 and 10. In the present study, the Ka/Ks ratios of the AP2/EREBP paralogous pairs indicated that the AP2/EREBP genes were influenced primarily by purifying selection, which indicated that the AP2/EREBP genes received strong environmental pressure during evolution. The Ka/Ks ratios of the COX-3/COX-4 paralogous pairs indicate that the COX-3/COX-4 genes were influenced primarily by Darwinian selection (driving change). For the LTP genes, the Ka/Ks ratios of the LTP-1/LTP-10, LTP-5/LTP-3 and LTP-4/LTP-8 paralogous pairs indicate that these genes were influenced primarily by purifying selection, while the Ka/Ks ratios of the LTP-2/LTP-6 paralogous pairs indicate that these genes were influenced primarily by Darwinian selection. The duplication time of the AP2/EREBP paralogous gene pairs in Z. mays L. ranged from approximately 9.364 to 100.935 Mya. The duplication time of the COX-3/COX-4 paralogous gene pair was approximately 5.217 Mya. The duplication time of the LTP paralogous gene pairs ranged from approximately 19.064 to 96.477 Mya. The major focus of research is to identify the genes that are responsible for drought stress tolerance to improve maize for drought stress tolerance. The results of the present study will improve the understanding of the functions of the AP2/EREBP, COX and LTP genes in response to drought stress.

Functional analysis of fasciclin-like arabinogalactan in carotenoid synthesis during tomato fruit ripening

Carotenoids are important pigmented nutrients synthesized by tomato fruits during ripening. To reveal the molecular mechanism underlying carotenoid synthesis during tomato fruit ripening, we analyzed carotenoid metabolites and transcriptomes in six development stages of tomato fruits. A total of thirty different carotenoids were detected and quantified in tomato fruits from 10 to 60 DPA. Based on differential gene expression profiles and WGCNA, we explored several genes that were highly significant and negatively correlated with lycopene, all of which encode fasciclin-like arabinogalactan proteins (FLAs). The FLAs are involved in plant signal transduction, however the functional role of these proteins has not been studied in tomato. Genome-wide analysis revealed that cultivated and wild tomato species contained 18 to 22 FLA family members, clustered into four groups, and mainly evolved by means of segmental duplication. The functional characterization of FLAs showed that silencing of SlFLA1, 5, and 13 were found to contribute to the early coloration of tomato fruits, and the expression of carotenoid synthesis-related genes was up-regulated in fruits that changed phenotypically, especially in SlFLA13-silenced plants. Furthermore, the content of multiple carotenoids (including (E/Z)-phytoene, lycopene, γ-carotene, and α-carotene) was significantly increased in SlFLA13-silenced fruits, suggesting that SlFLA13 has a potential inhibitory function in regulating carotenoid synthesis in tomato fruits. The results of the present study broaden the idea of analyzing the biological functions of tomato FLAs and preliminary evidence for the inhibitory role of SlFLA13 in carotenoid synthesis in fruit, providing the theoretical basis and a candidate for improving tomato fruit quality.

Genome-wide identification of the TIFY gene family in tobacco and expression analysis in response to Ralstonia solanacearum infection

The TIFY gene family plays an essential role in plant development and abiotic and biotic stress responses. In this study, genome-wide identification of TIFY members in tobacco and their expression pattern analysis in response to Ralstonia solanacearum infection were performed. A total of 33 TIFY genes were identified, including the TIFY, PPD, ZIM&ZML and JAZ subfamilies. Promoter analysis results indicated that a quantity of light-response, drought-response, SA-response and JA-response cis-elements exist in promoter regions. The TIFY gene family exhibited expansion and possessed gene redundancy resulting from tobacco ploidy change. In addition, most NtTIFYs equivalently expressed in roots, stems and leaves, while NtTIFY1, NtTIFY4, NtTIFY18 and NtTIFY30 preferentially expressed in roots. The JAZ III clade showed significant expression changes after inoculation with R. solanacearum, and the expression of NtTIFY7 in resistant varieties, compared with susceptible varieties, was more stably induced. Furthermore, NtTIFY7-silenced plants, compared with the control plants, were more susceptible to bacterial wilt. These results lay a foundation for exploring the evolutionary history of TIFY gene family and revealing gene function of NtTIFYs in tobacco bacterial wilt resistance.

Insights into sequence characteristics and evolutionary history of DGATs in arthropods

Triacylglycerol (TAG) is crucial in animal energy storage and membrane biogenesis. The conversion of diacylglycerol (DAG) to triacylglycerol (TAG) is catalyzed by diacylglycerol acyltransferase enzymes (DGATs), which are encoded by genes belonging to two distinct gene families. Although arthropods are known to possess DGATs activities and utilize the glycerol-3-phosphate pathway and MAG pathway for TAG biosynthesis, the sequence characterization and evolutionary history of DGATs in arthropods remains unclear. This study aimed to comparatively evaluate genomic analyses of DGATs in 13 arthropod species and 14 outgroup species. We found that arthropods lack SOAT2 genes within the DGAT1 family, while DGAT2, MOGAT3, AWAT1, and AWAT2 were absent from in DGAT2 family. Gene structure and phylogenetic analyses revealed that DGAT1 and DGAT2 genes come from different gene families. The expression patterns of these genes were further analyzed in crustaceans, demonstrating the importance of DGAT1 in TAG biosynthesis. Additionally, we identified the DGAT1 gene in Swimming crab (P. trituberculatus) undergoes a mutually exclusive alternative splicing event in the molt stages. Our newly determined DGAT inventory data provide a more complete scenario and insights into the evolutionary dynamics and functional diversification of DGATs in arthropods.

Genome-wide identification, structural characterization and gene expression analysis of the WRKY transcription factor family in pea (Pisum sativum L.)

BACKGROUND

The WRKY gene family is one of the largest families of transcription factors in higher plants, and WRKY transcription factors play important roles in plant growth and development as well as in response to abiotic stresses; however, the WRKY gene family in pea has not been systematically reported.

RESULTS

In this study, 89 pea WRKY genes were identified and named according to the random distribution of PsWRKY genes on seven chromosomes. The gene family was found to have nine pairs of tandem duplicates and 19 pairs of segment duplicates. Phylogenetic analyses of the PsWRKY and 60 Arabidopsis WRKY proteins were performed to determine their homology, and the PsWRKYs were classified into seven subfamilies. Analysis of the physicochemical properties, motif composition, and gene structure of pea WRKYs revealed significant differences in the physicochemical properties within the PsWRKY family; however, their gene structure and protein-conserved motifs were highly conserved among the subfamilies. To further investigate the evolutionary relationships of the PsWRKY family, we constructed comparative syntenic maps of pea with representative monocotyledonous and dicotyledonous plants and found that it was most recently homologous to the dicotyledonous WRKY gene families. Cis-acting element analysis of PsWRKY genes revealed that this gene family can respond to hormones, such as abscisic acid (ABA), indole-3-acetic acid (IAA), gibberellin (GA), methyl jasmonate (MeJA), and salicylic acid (SA). Further analysis of the expression of 14 PsWRKY genes from different subfamilies in different tissues and fruit developmental stages, as well as under five different hormone treatments, revealed differences in their expression patterns in the different tissues and fruit developmental stages, as well as under hormone treatments, suggesting that PsWRKY genes may have different physiological functions and respond to hormones.

CONCLUSIONS

In this study, we systematically identified WRKY genes in pea for the first time and further investigated their physicochemical properties, evolution, and expression patterns, providing a theoretical basis for future studies on the functional characterization of pea WRKY genes during plant growth and development.

Genome-wide identification of the Pyrus R2R3-MYB gene family and PhMYB62 regulation analysis in Pyrus hopeiensis flowers at low temperature

The R2R3-MYB gene family play an important role in plant growth, development and stress responses. In this study, a total of 122 PcoR2R3-MYB genes were identified and grouped into 26 clades in pear. And these PcoMYBs were unevenly distributed among 17 chromosomes. The sequence characteristics, conversed motifs, exon/intron structures, classification, duplication events and cis-acting elements were also investigated. The gene duplication events showed that segmental duplication may play key roles in expansion of the PcoMYB gene family. Pyrus hopeiensis, which is a valuable wild resource, has strong cold resistance. An integrative analyses of miRNA and mRNA showed that PhMYB62 was involved in regulating low-temperature stress in P. hopeiensis flower organs. Subcellular localization analysis showed that PhMYB62 protein was specifically localized to the nucleus. The result of DAP-seq showed that PhMYB62 responded to low-temperature stress in P. hopeiensis by regulating TFs, which were associated with plant stress resistance, and POD, GAUT12, AUX28 and CHS genes. Subsequently, yeast one-hybrid verified that PhMYB62 could bind and activate the promoter of POD gene. The current study would provide a comprehensive information for further functional research on the stress-responsive R2R3-MYB gene candidates in pear, and may help to identify the genes associated with cold resistance for the cultivation of cold-resistant pear varieties.

Evolution and expression analysis of carotenoid cleavage oxygenase gene family in Chinese mitten crab Eriocheir sinensis

Carotenoid cleavage oxygenase (CCO) plays a pivotal role in various biological activities, including antioxidant and immune functions in animals. This paper investigates the evolution and expression of CCO genes based on three chordates and 27 arthropods. Aquatic animals exhibit a higher abundance of CCO genes. Despite this, research on CCO in crustaceans has been notably limited, with a complete absence of any previous studies on the CCO genes for the Chinese mitten crab (Eriocheir sinensis). In this study, six CCO genes were identified in the E. sinensis genome database. Results reveal that the evolution of the CCO gene family in Crustacea is primarily characterized by purifying selection, with a preference for employing similar codons. EsCCO1 and EsCCO3 were mainly expressed in the epidermal layer, and EsCCO4 was mainly expressed in the hindgut. Meanwhile, EsCCO5 and EsCCO6 were mainly expressed in the hepatopancreas and endometrium. A notable detail that different EsCCO genes demonstrate distinct expression patterns within various tissues of E. sinensis. The findings of this study offer fundamental insights that could serve as a basis for further exploration into the functions and regulatory mechanisms of CCO genes in crustacean species.

Genome-wide identification and functional analysis of the cellulose synthase-like gene superfamily in common oat (Avena sativa L.)

Hemicelluloses constitute approximately one-third of the plant cell wall and can be used as a dietary fiber and food additive, and as raw materials for biofuels. Although genes involved in hemicelluloses synthesis have been investigated in some model plants, no comprehensive analysis has been conducted in common oat at present. In this study, we identified and systematically analyzed the cellulose synthase-like gene (Csl) family members in common oat and investigated them using various bioinformatics tools. The results showed that there are 76 members of the oat Csl gene family distributed on 17 chromosomes, and phylogenetic analysis indicated that the 76 Csl genes belong to the CslA, CslC, CslD, CslE, CslF, CslH, and CslJ subfamilies. A total of 14 classes of cis-acting elements were identified in the promoter regions, including hormone response, light response, cell development, and defense stress elements. The collinearity analysis identified 28 pairs of segmentally duplicated genes, most of which were found on chromosomes 2D and 6A. Expression pattern analysis showed that oat Csl genes display strong tissue-specific expression; of the 76 Csl genes, 33 were significantly up-regulated in stems and 30 were up-regulated in immature seeds. The expression of most members of the AsCsl gene family is repressed by abiotic stress, while the expression of some members is up-regulated by light. Immunoelectron microscopy shows that the product of AsCsl61, a member of CslF subfamily, mediates (1,3; 1,4)-β-D-glucan synthesis in transgenic Arabidopsis. These findings provide a fundamental understanding of the structural, functional, and evolutionary features of the oat Csl genes and may contribute to our general understanding of hemicellulose biosynthesis. Moreover, this information will be helpful in designing experiments for genetic manipulation of mixed-linkage glucan (MLG) synthesis with the goal of quality improvement in oat.

Genome-wide identification and expression characterization of the GH3 gene family of tea plant (Camellia sinensis)

To comprehensively understand the characteristics of the GH3 gene family in tea plants (Camellia sinensis), we identified 17 CsGH3 genes and analyzed their physicochemical properties, phylogenetic relationships, gene structures, promoters, and expression patterns in different tissues. The study showed that the 17 CsGH3 genes are distributed on 9 chromosomes, and based on evolutionary analysis, the CsGH3 members were divided into three subgroups. Gene duplication analysis revealed that segmental duplications have a significant impact on the amplification of CsGH3 genes. In addition, we identified and classified cis-elements in the CsGH3 gene promoters and detected elements related to plant hormone responses and non-biotic stress responses. Through expression pattern analysis, we observed tissue-specific expression of CsGH3.3 and CsGH3.10 in flower buds and roots. Moreover, based on predictive analysis of upstream regulatory transcription factors of CsGH3, we identified the potential transcriptional regulatory role of gibberellin response factor CsDELLA in CsGH3.14 and CsGH3.15. In this study, we found that CsGH3 genes are involved in a wide range of activities, such as growth and development, stress response, and transcription. This is the first report on CsGH3 genes and their potential roles in tea plants. In conclusion, these results provide a theoretical basis for elucidating the role of GH3 genes in the development of perennial woody plants and offer new insights into the synergistic effects of multiple hormones on plant growth and development in tea plants.

Genome-wide identification of the GRF family in sweet orange (Citrus sinensis) and functional analysis of the CsGRF04 in response to multiple abiotic stresses

BACKGROUND

Citrus is one of the most valuable fruits worldwide and an economic pillar industry in southern China. Nevertheless, it frequently suffers from undesirable environmental stresses during the growth cycle, which severely restricts the growth, development and yield of citrus. In plants, the growth-regulating factor (GRF) family of transcription factors (TF) is extensively distributed and plays an vital part in plant growth and development, hormone response, as well as stress adaptation. However, the systematic identification and functional analysis of GRF TFs in citrus have not been reported.

RESULTS

Here, a genome-wide identification of GRF TFs was performed in Citrus sinensis, 9 members of CsGRFs were systematically identified and discovered to be scattered throughout 5 chromosomes. Subsequently, physical and chemical properties, phylogenetic relationships, structural characteristics, gene duplication events, collinearity and cis-elements of promoter were elaborately analyzed. In particular, the expression patterns of the CsGRF genes in response to multiple phytohormone and abiotic stress treatments were investigated. Predicated on this result, CsGRF04, which exhibited the most differential expression pattern under multiple phytohormone and abiotic stress treatments was screened out. Virus-induced gene silencing (VIGS) technology was utilized to obtain gene silenced plants for CsGRF04 successfully. After the three stress treatments of high salinity, low temperature and drought, the CsGRF04-VIGS lines showed significantly reduced resistance to high salinity and low temperature stresses, but extremely increased resistance to drought stress.

CONCLUSIONS

Taken together, our findings systematically analyzed the genomic characterization of GRF family in Citrus sinensis, and excavated a CsGRF04 with potential functions under multiple abiotic stresses. Our study lay a foundation for further study on the function of CsGRFs in abiotic stress and hormone signaling response.

Genome-wide identification and analysis of the GGCT gene family in wheat

BACKGROUND

γ-glutamylcyclotransferase (GGCT), an enzyme to maintain glutathione homeostasis, plays a vital role in the response to plant growth and development as well as the adaptation to various stresses. Although the GGCT gene family analysis has been conducted in Arabidopsis and rice, the family genes have not yet been well identified and analyzed at the genome-wide level in wheat (Triticum aestivum L.).

RESULTS

In the present study, 20 TaGGCT genes were identified in the wheat genome and widely distributed on chromosomes 2A, 2B, 2D, 3A, 4A, 5A, 5B, 5D, 6A, 6B, 6D, 7A, 7B, and 7D. Phylogenetic and structural analyses showed that these TaGGCT genes could be classified into three subfamilies: ChaC, GGGACT, and GGCT-PS. They exhibited similar motif compositions and distribution patterns in the same subgroup. Gene duplication analysis suggested that the expansion of TaGGCT family genes was facilitated by segmental duplications and tandem repeats in the wheat evolutionary events. Identification of diverse cis-acting response elements in TaGGCT promoters indicated their potential fundamental roles in response to plant development and abiotic stresses. The analysis of transcriptome data combined with RT-qPCR results revealed that the TaGGCTs genes exhibited ubiquitous expression across plant organs, with highly expressed in roots, stems, and developing grains. Most TaGGCT genes were up-regulated after 6 h under 20% PEG6000 and ABA treatments. Association analysis revealed that two haplotypes of TaGGCT20 gene displayed significantly different Thousand-kernel weight (TKW), Kernel length (KL), and Kernel width (KW) in wheat. The geographical and annual distribution of the two haplotypes of TaGGCT20 gene further revealed that the frequency of the favorable haplotype TaGGCT20-Hap-I was positively selected in the historical breeding process of wheat.

CONCLUSION

This study investigated the genome-wide identification, structure, evolution, and expression analysis of TaGGCT genes in wheat. The motifs of TaGGCTs were highly conserved throughout the evolutionary history of wheat. Most TaGGCT genes were highly expressed in roots, stems, and developing grains, and involved in the response to drought stresses. Two haplotypes were developed in the TaGGCT20 gene, where TaGGCT20-Hap-I, as a favorable haplotype, was significantly associated with higher TKW, KL, and KW in wheat, suggesting that the haplotype is used as a function marker for the selection in grain yield in wheat breeding.

Insights into ZmWAKL in maize kernel development: genome-wide investigation and GA-mediated transcription

BACKGROUND

The functional roles of the Wall Associated Kinase (WAK) and Wall Associated Kinase Like (WAKL) families in cellular expansion and developmental processes have been well-established. However, the molecular regulation of these kinases in maize development is limited due to the absence of comprehensive genome-wide studies.

RESULTS

Through an in-depth analysis, we identified 58 maize WAKL genes, and classified them into three distinct phylogenetic clusters. Moreover, structural prediction analysis showed functional conservation among WAKLs across maize. Promoter analysis uncovered the existence of cis-acting elements associated with the transcriptional regulation of ZmWAKL genes by Gibberellic acid (GA). To further elucidate the role of WAKL genes in maize kernels, we focused on three highly expressed genes, viz ZmWAKL38, ZmWAKL42 and ZmWAKL52. Co-expression analyses revealed that their expression patterns exhibited a remarkable correlation with GA-responsive transcription factors (TF) TF5, TF6, and TF8, which displayed preferential expression in kernels. RT-qPCR analysis validated the upregulation of ZmWAKL38, ZmWAKL42, ZmWAKL52, TF5, TF6, and TF8 following GA treatment. Additionally, ZmWAKL52 showed significant increase of transcription in the present of TF8, with ZmWAKL52 localizing in both the plasma membrane and cell wall. TF5 positively regulated ZmWAKL38, while TF6 positively regulated ZmWAKL42.

CONCLUSIONS

Collectively, these findings provide novel insights into the characterization and regulatory mechanisms of specific ZmWAKL genes involved in maize kernel development, offering prospects for their utilization in maize breeding programs.

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.

Genome-wide analysis of PIN genes in cultivated peanuts (Arachis hypogaea L.): identification, subcellular localization, evolution, and expression patterns

BACKGROUND

Auxin is an important hormone in plants and the PIN-FORMED (PIN) genes are essential to auxin distribution in growth and developmental processes of plants. Peanut is an influential cash crop, but research into PIN genes in peanuts remains limited.

RESULTS

In this study, 16 PIN genes were identified in the genome of cultivated peanut, resolving into four subfamilies. All PIN genes were predicted to be located in the plasma membrane and a subcellular location experiment confirmed this prediction for eight of them. The gene structure, cis-elements in the promoter, and evolutionary relationships were elucidated, facilitating our understanding of peanut PINs and their evolution. In addition, the expression patterns of these PINs in various tissues were analyzed according to a previously published transcriptome dataset and qRT-PCR, which gave us a clear understanding of the temporal and spatial expression of PIN genes in different growth stages and different tissues. The expression trend of homologous genes was similar. AhPIN2A and AhPIN2B exhibited predominant expression in roots. AhPIN1A-1 and AhPIN1B-1 displayed significant upregulation following peg penetration, suggesting a potential close association with peanut pod development. Furthermore, we presented the gene network and gene ontology enrichment of these PINs. Notably, AhABCB19 exhibited a co-expression relationship with AhPIN1A and AhPIN1B-1, with all three genes displaying higher expression levels in peanut pegs and pods. These findings reinforce their potential role in peanut pod development.

CONCLUSIONS

This study details a comprehensive analysis of PIN genes in cultivated peanuts and lays the foundation for subsequent studies of peanut gene function and phenotype.

Function identification and characterization of Oryza sativa ZRT and IRT-like proteins computationally for nutrition and biofortification in rice

Zinc plays a very critical role and function in all organisms. Its deficiency can cause a serious issue. In Oryza sativa, the ZRT/IRT transporter-like proteins play a role in the zinc metal uptake and transport. Few OsZIPs genes have been validated and characterized for their biological functions and most of OsZIPs are not well physiologically, biochemically and phenotypically characterized. In the current study, they analyzed for their function through subcellular localization, phylogenetic analysis, homology modeling, expression analysis, protein-protein interaction (PPI) network prediction, and prediction of their binding sites. Hierarchical clustering of OsZIP genes based on different anatomical parts and developmental stages also orthologs prediction was identified. The presence of SNPs, SSRs, ESTs, FSTs, MPSS, and SAGE tags were analyzed for useful development of markers. SNPs were identified in all OsZIPs genes and each gene was further classified based on their number and position in the 3'UTR and 5'UTR regions of the gene-specific sequences. Binding clusters and their location on the protein sequences were predicted. We found Changing in residues number and position which were due to partial overlapping and sequence alignment, but they share the same mechanism of binding and transporting Zinc. A wide range of CRISPR Cas9 gRNAs was designed based on single nucleotide polymorphism (SNP) for each OsZIP transporter gene for well-function identification and characterization with genome-wide association studies. Hence this study would provide useful information, understanding, and predicting molecular insights for the future studies that will help for improvement of nutritional quality of rice varieties.Communicated by Ramaswamy H. Sarma.

Genome-wide identification and analysis of the WRKY gene family and low-temperature stress response in Prunus sibirica

BACKGROUND

WRKY transcription factors are a prominent gene family in plants, playing a crucial role in various biological processes including development, metabolism, defense, differentiation, and stress response. Although the WRKY gene family has been extensively studied and analysed in numerous plant species, research on Prunus sibirica's WRKY genes (PsWRKY) remains lacking.

RESULTS

This study analysed the basic physicochemical properties, phylogeny, gene structure, cis-acting elements, and Gene ontology (GO) annotation of PsWRKY gene family members using bioinformatics methods based on the whole-genome data of P. sibirica. In total, 55 WRKYs were identified in P. sibirica and were heterogeneously distributed on eight chromosomes. Based on the phylogenetic analysis, these WRKYs were classified into three major groups: Group I, Group II (II-a, II-b, II-c, II-d, II-e), and Group III. Members of different subfamilies have different cis-acting elements, conserved motifs, and intron-exon structures, indicating functional heterogeneity of the WRKY family. Prediction of subcellular localisation indicated that PsWRKYs were mainly located in the nucleus. Twenty pairs of duplicated genes were identified, and segmental duplication events may play an important role in PsWRKY gene family expansion. Analysis of the Ka/Ks ratio showed that the PsWRKY family's homologous genes were primarily purified by selection. Additionally, GO annotation analysis showed that the WRKY gene family was mainly involved in responses to stimuli, immune system processes, and reproductive processes. Furthermore, quantitative real-time PCR (qRT-PCR) analysis showed that 23 PsWRKYs were highly expressed in one or more tissues (pistils and roots) and PsWRKYs showed specific expression patterns under different low-temperature stress conditions.

CONCLUSIONS

Our results provide a scientific basis for the further exploration and functional validation of WRKYs in P. sibirica.

Genome-wide identification and expression analysis of the polygalacturonase gene family in sweetpotato

BACKGROUND

Polygalacturonase (PG), a crucial enzyme involved in pectin degradation, is associated with various plants' developmental and physiological processes such as seed germination, fruit ripening, fruit softening and plant organ abscission. However, the members of PG gene family in sweetpotato (Ipomoea batatas) have not been extensively identified.

RESULTS

In this study, there were 103 PG genes identified in sweetpotato genome, which were phylogenetically clustered into divergent six clades. The gene structure characteristics of each clade were basically conserved. Subsequently, we renamed these PGs according to their locations of the chromosomes. The investigation of collinearity between the PGs in sweetpotato and other four species, contained Arabidopsis thaliana, Solanum lycopersicum, Malus domestica and Ziziphus jujuba, revealed important clues about the potential evolution of the PG family in sweetpotato. Gene duplication analysis showed that IbPGs with collinearity relationships were all derived from segmental duplications, and these genes were under purifying selection. In addition, each promoter region of IbPG proteins contained cis-acting elements related to plant growth and development processes, environmental stress responses and hormone responses. Furthermore, the 103 IbPGs were differentially expressed in various tissues (leaf, stem, proximal end, distal end, root body, root stalk, initiative storage root and fibrous root) and under different abiotic stresses (salt, drought, cold, SA, MeJa and ABA treatment). IbPG038 and IbPG039 were down-regulated with salt, SA and MeJa treatment. According to the further investigation, we found that IbPG006, IbPG034 and IbPG099 had different patterns under the drought and salt stress in fibrous root of sweetpotato, which provided insights into functional differences among these genes.

CONCLUSION

A total of 103 IbPGs were identified and classified into six clades from sweetpotato genome. The results of RNA-Seq and qRT-PCR suggested that IbPG006, IbPG034 and IbPG099 might play a significant role in tissue specificity as well as drought and salt stress responses, which showed valuable information for further functional characterization and application of the IbPGs.

Global dissection of R2R3-MYB in Pogostemon cablin uncovers a species-specific R2R3-MYB clade

MYB family is one of the largest transcription factor families in plants and plays a crucial role in regulating plant biochemical and physiological processes. However, R2R3-MYBs in patchouli have not been systematically investigated. Here, based on the gene annotation of patchouli genome sequence, 484 R2R3-MYB transcripts were detected. Further in-depth analysis of the gene structure and expression of R2R3-MYBs supported the tetraploid hybrid origin of patchouli. When combined with R2R3-MYBs from Arabidopsis, a phylogenetic tree of patchouli R2R3-MYBs was constructed and divided into 31 clades. Interestingly, a patchouli-specific R2R3-MYB clade was found and confirmed by homologous from other Lamiaceae species. The syntenic analysis demonstrated that tandem duplication contributed to its evolution. This study systematically analysed the R2R3-MYB family in patchouli, providing information on its gene characterization, functional prediction, and species evolution.

Comprehensive analysis of OpHD-ZIP transcription factors related to the regulation of camptothecin biosynthesis in Ophiorrhiza pumila

Ophiorrhiza pumila, as a folk herb belonging to the Rubiaceae family, has become a potential source of camptothecin (CPT), which is a monoterpenoid indole alkaloid with good antitumor property. However, the camptothecin content in this herb is low, and is far from meeting the increasing clinical demand. Understanding the transcriptional regulation of camptothecin biosynthesis provides an effective strategy for improvement of camptothecin yield. Previous studies have demonstrated several transcription factors that are related to camptothecin biosynthesis, while the functions of HD-ZIP members in O. pumila have not been investigated yet. In this study, 32 OpHD-ZIP transcription factor members were genome-wide identified. Phylogenetic tree showed that these OpHD-ZIP proteins are divided into four subfamilies. Based on the transcriptome data, nine OpHD-ZIP genes were shown to be predominantly expressed in O. pumila roots, which were in line with the camptothecin biosynthetic genes. Co-expression analysis showed that OpHD-ZIP7 and OpHD-ZIP20 were potentially related to the modulation of camptothecin biosynthesis. Dual-luciferase reporter assays (Dual-LUC) showed that both OpHD-ZIP7 and OpHD-ZIP20 could activate the expression of camptothecin biosynthetic genes OpIO and OpTDC. In conclusion, this study offered the promising data for exploring the roles of OpHD-ZIP transcription factors in regulating camptothecin biosynthesis.

Genome-wide identification and expression profile analysis of metal tolerance protein gene family in Eucalyptus grandis under metal stresses

Metal tolerance proteins (MTPs) as Me²⁺/H⁺(K⁺) antiporters participate in the transport of divalent cations, leading to heavy metal stress resistance and mineral utilization in plants. In the present study, to obtain better knowledge of the biological functions of the MTPs family, 20 potential EgMTPs genes were identified in Eucalyptus grandis and classified into seven groups belonging to three cation diffusion facilitator groups (Mn-CDFs, Zn/Fe-CDFs, and Zn-CDFs) and seven groups. EgMTP-encoded amino acids ranged from 315 to 884, and most of them contained 4-6 recognized transmembrane domains and were clearly prognosticated to localize into the cell vacuole. Almost all EgMTP genes experienced gene duplication events, in which some might be uniformly distributed in the genome. The numbers of cation efflux and the zinc transporter dimerization domain were highest in EgMTP proteins. The promoter regions of EgMTP genes have different cis-regulatory elements, indicating that the transcription rate of EgMTP genes can be a controlled response to different stimuli in multiple pathways. Our findings provide accurate perception on the role of the predicted miRNAs and the presence of SSR marker in the Eucalyptus genome and clarify their functions in metal tolerance regulation and marker-assisted selection, respectively. Gene expression profiling based on previous RNA-seq data indicates a probable function for EgMTP genes during development and responses to biotic stress. Additionally, the upregulation of EgMTP6, EgMTP5, and EgMTP11.1 to excess Cd²⁺ and Cu²⁺ exposure might be responsible for metal translocation from roots to leaves.

Identification and transcriptome data analysis of ARF family genes in five Orchidaceae species

The Orchidaceae is a large family of perennial herbs especially noted for the exceptional diversity of specialized flowers. Elucidating the genetic regulation of flowering and seed development of orchids is an important research goal with potential utility in orchid breeding programs. Auxin Response Factor (ARF) genes encode auxin-responsive transcription factors, which are involved in the regulation of diverse morphogenetic processes, including flowering and seed development. However, limited information on the ARF gene family in the Orchidaceae is available. In this study, 112 ARF genes were identified in the genomes of 5 orchid species (Apostasia shenzhenica, Dendrobium catenatum, Phalaenopsis aphrodite, Phalaenopsis equestris and Vanilla planifolia,). These genes were grouped into 7 subfamilies based on their phylogenetic relationships. Compared with the ARF family in model plants, such as Arabidopsis thaliana and Oryza sativa, one group of ARF genes involved in pollen wall synthesis has been lost during evolution of the Orchidaceae. This loss corresponds with absence of the exine in the pollinia. Through mining of the published genomic and transcriptomic data for the 5 orchid species: the ARF genes of subfamily 4 may play an important role in flower formation and plant growth, whereas those of subfamily 3 are potentially involved in pollen wall development. the study results provide novel insights into the genetic regulation of unique morphogenetic phenomena of orchids, which lay a foundation for further analysis of the regulatory mechanisms and functions of sexual reproduction-related genes in orchids.

Genome-wide identification and molecular characterization of CRK gene family in cucumber (Cucumis sativus L.) under cold stress and sclerotium rolfsii infection

BACKGROUND

The plant cysteine-rich receptor-like kinases (CRKs) are a large family having multiple roles, including defense responses under both biotic and abiotic stress. However, the CRK family in cucumbers (Cucumis sativus L.) has been explored to a limited extent. In this study, a genome-wide characterization of the CRK family has been performed to investigate the structural and functional attributes of the cucumber CRKs under cold and fungal pathogen stress.

RESULTS

A total of 15 C. sativus CRKs (CsCRKs) have been characterized in the cucumber genome. Chromosome mapping of the CsCRKs revealed that 15 genes are distributed in cucumber chromosomes. Additionally, the gene duplication analysis of the CsCRKs yielded information on their divergence and expansion in cucumbers. Phylogenetic analysis divided the CsCRKs into two clades along with other plant CRKs. Functional predictions of the CsCRKs suggested their role in signaling and defense response in cucumbers. The expression analysis of the CsCRKs by using transcriptome data and via qRT-PCR indicated their involvement in both biotic and abiotic stress responses. Under the cucumber neck rot pathogen, Sclerotium rolfsii infection, multiple CsCRKs exhibited induced expressions at early, late, and both stages. Finally, the protein interaction network prediction results identified some key possible interacting partners of the CsCRKs in regulating cucumber physiological processes.

CONCLUSIONS

The results of this study identified and characterized the CRK gene family in cucumbers. Functional predictions and validation via expression analysis confirmed the involvement of the CsCRKs in cucumber defense response, especially against S. rolfsii. Moreover, current findings provide better insights into the cucumber CRKs and their involvement in defense responses.

Genome-wide identification of the fibrillin gene family in chickpea (Cicer arietinum L.) and its response to drought stress

Fibrillin family members play multiple roles in growth, development, and protection against abiotic stress. In this study, we identified 12 potential CaFBNs that are ranging from 25 kDa-42.92 kDa and are mostly basic. These proteins were hydrophilic in nature and generally resided in the chloroplast. The CaFBN genes were located on different chromosomes like 1, 4, 5, and 7. All FBNs shared conserved motifs and possessed a higher number of stress-responsive elements. For evolutionary analysis, a phylogenetic tree of CaFBNs with other plants' FBNs was constructed and clustered into 11 FBN subgroups. For expression analysis, 21 day old chickpea seedling was exposed to dehydration stress by withholding water. We also performed various physiological and biochemical analyses to check that plant changes at the physiological and cellular levels while undergoing stress conditions. The transcript expression of CaFBNs was higher in aerial parts, especially in stems and leaves. Dehydration-specific transcriptome and qPCR analysis showed that FBN-1, FBN-2, and FBN-6 were highly expressed. In addition, our study provides a comprehensive overview of the FBN protein family and their importance during the dehydration stress condition in Cicer arietinum.

Genome-wide identification and comparative expression profiling of the WRKY transcription factor family in two Citrus species with different Candidatus Liberibacter asiaticus susceptibility

BACKGROUND

Salicylic Acid (SA) is a pivotal phytohormone in plant innate immunity enhancement of triggered by various pathogens, such as Candidatus Liberibacter asiaticus (CLas), the causal agent of Huanglongbing (HLB). WRKY is a plant specific transcription factor (TF) family, which plays crucial roles in plant response to biotic stresses. So far, the evolutionary history, functions, and expression patterns under SA treatment and CLas infection of WRKY family are poorly understood in Citrus, despite the release of the genome of several Citrus species. A comprehensive genomic and expressional analysis is worth to conduct for this family.

RESULTS

Here, a genome-wide identification of WRKY TFs was performed in two Citrus species: Citrus sinensis (HLB-sensitive) and Poncirus trifoliata (HLB-tolerant). In total, 52 CsWRKYs and 51 PtrWRKYs were identified, whose physical and chemical properties, chromosome locations, phylogenetic relationships and structural characteristics were comparatively analyzed. Especially, expression patterns of these WRKY genes before and after SA treatment and CLas infection were compared. Based on this result, seven pairs of orthologous WRKY genes showing opposite expression patterns in two Citrus species were screened out. Moreover, two pairs of orthologous WRKY genes with significant differences in the number or type of stress-responsive cis-elements in the promoter regions were discovered. Subcellular localization and transcriptional activation activity assays revealed that these two pairs of orthologous genes are classic WRKY TFs localize in the nucleus and could function as transcriptional activators.

CONCLUSION

In this study, we systematically analyzed the genomic characterization of WRKY family in two Citrus species, together with the analyses of expression patterns under SA signaling and CLas infection. Our study laid a foundation for further study on the function of WRKY TFs in HLB response and SA signaling of Citrus.

Identification of HubHLH family and key role of HubHLH159 in betalain biosynthesis by activating the transcription of HuADH1, HuCYP76AD1-1, and HuDODA1 in pitaya

Basic helix-loop-helix (bHLH) proteins are dimeric transcription factors (TFs) involved in various plant physiological and biological processes. Despite this, little is known about the molecular properties and roles of bHLH TFs in pitaya betalain biosynthesis. Here we report the identification of 165 HubHLH genes in H. undantus genome, their chromosomal distribution, physiochemical characteristics, conserved motifs, gene structure, phylogeny and synteny of HubHLH genes. Based on phylogenetic relationship analysis, the 165 HubHLHs were divided into 26 subfamilies and unequally distributed on the 11 chromosomes of pitaya. Based on the pitaya transcriptome data, a candidate gene HubHLH159 was obtained, and the real-time quantitative PCR analysis confirmed that HubHLH159 showed a high expression level in 'Guanhuahong' pitaya (red-pulp) at mature stage, indicating its role in betalain biosynthesis. HubHLH159 is a Group II protein and contains a bHLH domain. It is a nuclear protein with transcriptional activation activity. Dual luciferase reporter assays and virus-induced gene silencing (VIGS) experiments showed that HubHLH159 promotes betalain biosynthesis by activating the expression of HuADH1, HuCYP76AD1-1, and HuDODA1. The results of the present study lay a new theoretical reference for the regulation of pitaya betalain biosynthesis and also provides as essential basis for the future analysis of the functions of HubHLH gene family.

Genome-wide identification glutathione-S-transferase gene superfamily in Daphnia pulex and its transcriptional response to nanoplastics

Glutathione S-transferases (GSTs) are key multifunctional phase II detoxification enzymes involved in the regulation of growth, development, and stress responses. However, the knowledge of GSTs in the model invertebrate organism Daphnia pulex at the genomic level remains limited. In the present study, 35 GST genes were identified in D. pulex (Dp-GST), belonging to eight subfamilies, with the sigma, mu, and delta/epsilon subfamilies constituting approximately 29 %, 20 %, and 20 % of the GST superfamily, respectively. Chromosome tandem duplication of genes within the same subfamily was observed, which may be the main force driving GST expansion in D. pulex. DpGST genes showed different expression patterns in response to nanoplastic exposure for 96 h and 21 days. Some homologous GST genes in D. pulex showed similar expression patterns in response to nanoplastic exposure, likely owing to their unique motifs. For example, motif 9 is found in all delta/epsilon GST genes, whereas motifs 1, 2, 3, 5, and 7 are highly conserved in sigma GST genes. The characterization of D. pulex GSTs extending the knowledge of GST-mediated environmental contaminants, especially nanoplastics.

Comparison of microsatellite distribution in the genomes of Pteropus vampyrus and Miniopterus natalensis (Chiroptera)

BACKGROUND

Microsatellites are a ubiquitous occurrence in prokaryotic and eukaryotic genomes. Microsatellites have become one of the most popular classes of genetic markers due to their high reproducibility, multi-allelic nature, co-dominant mode of inheritance, abundance and wide genome coverage. We characterised microsatellites in the genomes and genes of two bat species, Pteropus vampyrus and Miniopterus natalensis. This characterisation was used for gene ontology analysis and the Kyoto Encyclopedia of Genes and Genomes pathway enrichment of coding sequences (CDS).

RESULTS

Compared to M. natalensis, the genome size of P. vampyrus is larger and contains more microsatellites, but the total diversity of both species is similar. Mononucleotide and dinucleotide repeats were the most diverse in the genome of the two species. In each bat species, the microsatellite bias was obvious. The microsatellites with the largest number of repeat motifs in P. vampyrus from mononucleotide to hexanucleotide were (A)_n, (AC)_n, (CAA)_n, (AAAC)_n, (AACAA)_n and (AAACAA)_n, with frequencies of 97.94%, 58.75%, 30.53%, 22.82%, 54.68% and 22.87%, respectively, while in M. natalensis were (A)_n, (AC)_n, (TAT)_n, (TTTA)_n, (AACAA)_n and (GAGAGG)_n, with of 92.00%, 34.08%, 40.36%, 21.83%, 25.42% and 12.79%, respectively. In both species, the diversity of microsatellites was highest in intergenic regions, followed by intronic, untranslated and exonic regions and lowest in coding regions. Location analysis indicated that microsatellites were mainly concentrated at both ends of the genes. Microsatellites in the CDS are thus subject to higher selective pressure. In the GO analysis, two unique GO terms were found only in P. vampyrus and M. natalensis, respectively. In KEGG enriched pathway, the biosynthesis of other secondary metabolites and metabolism of other amino acids in metabolism pathways were present only in M. natalensis. The combined biological process, cellular components and molecular function ontology are reflected in the GO analysis and six functional enrichments in KEGG annotation, suggesting advantageous mutations during species evolution.

CONCLUSIONS

Our study gives a comparative characterization of the genomes of microsatellites composition in the two bat species. And also allow further study on the effect of microsatellites on gene function as well as provide an insight into the molecular basis for species adaptation to new and changing environments.

Genome-wide identification and characterization of polycomb repressive complex 2 core components in upland cotton (Gossypium hirsutum L.)

BACKGROUND

The evolutionarily conserved Polycomb Repressive Complex 2 (PRC2) plays a vital role in epigenetic gene repression by depositing tri-methylation on lysine residue K27 of histone H3 (H3K27me3) at the target loci, thus participating in diverse biological processes. However, few reports about PRC2 are available in plant species with large and complicated genomes, like cotton.

RESULTS

Here, we performed a genome-wide identification and comprehensive analysis of cotton PRC2 core components, especially in upland cotton (Gossypium hirsutum). Firstly, a total of 8 and 16 PRC2 core components were identified in diploid and tetraploid cotton species, respectively. These components were classified into four groups, E(z), Su(z)12, ESC and p55, and the members in the same group displayed good collinearity, similar gene structure and domain organization. Next, we cloned G. hirsutum PRC2 (GhPRC2) core components, and found that most of GhPRC2 proteins were localized in the nucleus, and interacted with each other to form multi-subunit complexes. Moreover, we analyzed the expression profile of GhPRC2 genes. The transcriptome data and quantitative real-time PCR (qRT-PCR) assays indicated that GhPRC2 genes were ubiquitously but differentially expressed in various tissues, with high expression levels in reproductive organs like petals, stamens and pistils. And the expressions of several GhPRC2 genes, especially E(z) group genes, were responsive to various abiotic and biotic stresses, including drought, salinity, extreme temperature, and Verticillium dahliae (Vd) infection.

CONCLUSION

We identified PRC2 core components in upland cotton, and systematically investigated their classifications, phylogenetic and synteny relationships, gene structures, domain organizations, subcellular localizations, protein interactions, tissue-specific and stresses-responsive expression patterns. Our results will provide insights into the evolution and composition of cotton PRC2, and lay the foundation for further investigation of their biological functions and regulatory mechanisms.

Zinc supply influenced the growth, yield, zinc content, and expression of ZIP family transporters in sorghum

Zinc deficiency altered shoot and root growth, plant biomass, yield, and ZIP family transporter gene expression in sorghum. Zinc (Zn) deficiency affects several crop plants' growth and yield, including sorghum. We have evaluated the sorghum under various concentrations of Zn supply for phenotypic changes, Zn content, and expression of Zn-regulated, iron-regulated transporter-like proteins (ZIP) family genes. Zn deficiency reduced the shoot and root growth, plant biomass, and yield by > 50%. The length and number of lateral roots were increased by more than 50% under deficient Zn compared to sufficient Zn. Ten SbZIP family transporter genes showed dynamic expression in shoot and root tissues of sorghum under deficient and sufficient Zn. SbZIP2, 5, 6, 7, and 8 were expressed in all tissues under deficient and sufficient Zn. SbZIP2, 4, 5, 6, 7, 8, and 10 were highly induced in shoot tissues by deficient Zn. The expression level of SbZIP6, 7, 8, and 9 was higher in root tissues under deficient Zn. The phylogenetic analysis revealed that most SbZIP family proteins are closely associated with the ZmZIP family of maize. The functional residues His177 and Gly182 are fully conserved in all SbZIP family transporters, as revealed by homology modeling and multiple sequence alignment. This study may provide a foundation for improving the Zn-use efficiency of sorghum.

The LEA gene family in tomato and its wild relatives: genome-wide identification, structural characterization, expression profiling, and role of SlLEA6 in drought stress

BACKGROUND

Late embryogenesis abundant (LEA) proteins are widely distributed in higher plants and play crucial roles in regulating plant growth and development processes and resisting abiotic stress. Cultivated tomato (Solanum lycopersicum) is an important vegetable crop worldwide; however, its growth, development, yield, and quality are currently severely constrained by abiotic stressors. In contrast, wild tomato species are more tolerant to abiotic stress and can grow normally in extreme environments. The main objective of this study was to identify, characterize, and perform gene expression analysis of LEA protein families from cultivated and wild tomato species to mine candidate genes and determine their potential role in abiotic stress tolerance in tomatoes.

RESULTS

Total 60, 69, 65, and 60 LEA genes were identified in S. lycopersicum, Solanum pimpinellifolium, Solanum pennellii, and Solanum lycopersicoides, respectively. Characterization results showed that these genes could be divided into eight clusters, with the LEA_2 cluster having the most members. Most LEA genes had few introns and were non-randomly distributed on chromosomes; the promoter regions contained numerous cis-acting regulatory elements related to abiotic stress tolerance and phytohormone responses. Evolutionary analysis showed that LEA genes were highly conserved and that the segmental duplication event played an important role in evolution of the LEA gene family. Transcription and expression pattern analyses revealed different regulatory patterns of LEA genes between cultivated and wild tomato species under normal conditions. Certain S. lycopersicum LEA (SlLEA) genes showed similar expression patterns and played specific roles under different abiotic stress and phytohormone treatments. Gene ontology and protein interaction analyses showed that most LEA genes acted in response to abiotic stimuli and water deficit. Five SlLEA proteins were found to interact with 11 S. lycopersicum WRKY proteins involved in development or resistance to stress. Virus-induced gene silencing of SlLEA6 affected the antioxidant and reactive oxygen species defense systems, increased the degree of cellular damage, and reduced drought resistance in S. lycopersicum.

CONCLUSION

These findings provide comprehensive information on LEA proteins in cultivated and wild tomato species and their possible functions under different abiotic and phytohormone stresses. The study systematically broadens our current understanding of LEA proteins and candidate genes and provides a theoretical basis for future functional studies aimed at improving stress resistance in tomato.

Genome-wide identification and gene expression analysis of the 14-3-3 gene family in potato (Solanum tuberosum L.)

BACKGROUND

14-3-3 proteins are essential in regulating various biological processes and abiotic stress responses in plants. Although 14-3-3 proteins have been studied in model plants such as Arabidopsis thaliana and Oryza sativa, there is a lack of research on the 14-3-3 gene family in potatoes (Solanum tuberosum L.).

RESULTS

A total of 18 14-3-3 genes encoding proteins containing a typical conserved PF00244 domain were identified by genome-wide analysis in potatoes. The St14-3-3 gene family members were unevenly distributed across the chromosomes, and gene structure analysis showed that gene length and intron number varied greatly among the members. Phylogenetic analysis of 14-3-3 proteins in potatoes and other plant species showed that they could be divided into two distinct groups (ε and non-ε). Members in the ε group tended to have similar exon-intron structures and conserved motif patterns. Promoter sequence analysis showed that the St14-3-3 gene promoters contained multiple hormone-, stress-, and light-responsive cis-regulatory elements. Synteny analysis suggested that segmental duplication events contributed to the expansion of the St14-3-3 gene family in potatoes. The observed syntenic relationships between some 14-3-3 genes from potato, Arabidopsis, and tomato suggest that they evolved from a common ancestor. RNA-seq data showed that St14-3-3 genes were expressed in all tissues of potatoes but that their expression patterns were different. qRT-PCR assays revealed that the expression levels of nearly all tested St14-3-3 genes were affected by drought, salt, and low-temperature stresses and that different St14-3-3 genes had different responses to these stresses.

CONCLUSIONS

In summary, genome-wide identification, evolutionary, and expression analyses of the 14-3-3 gene family in potato were conducted. These results provide important information for further studies on the function and regulation of St14-3-3 gene family members in potatoes.

Beta-galactosidase gene family genome-wide identification and expression analysis of members related to fruit softening in melon (Cucumis melo L.)

BACKGROUND

Texture quality is impotent for melon (Cucumis melo L.) fruit. β-galactosidase (β-Gal, EC 3.2.1.23) is an important cell wall glycosyl hydrolase involved in fruit softening, However, the β-Gal gene (BGALs) family hasn't been identified genome-wide in melon. Thus, it's necessary to conduct an in-depth bioinformatic analysis on melon BGALs family and to seek out the key members who participated in melon fruit softening.

RESULTS

A total of 21 BGALs members designated as CmBGAL1-CmBGAL21 were identified genome-wide in melon, clustered into A-G seven clades. Among them, three duplications CmBGAL1:CmBGAL3, CmBGAL19:CmBGAL21, and CmBGAL20:CmBGAL21 happened. For conserved domains, besides the Glyco_hydro_35 domain (PF01301), all the members also contained the GHD domain (PF17834) except for CmBGAL12, and the Gal_Lectin (PF02140) domain existed in most CmBGALs at the C-termini. Motifs, protein secondary and tertiary structure analysis showed that the CmBGAL12 is a unique member. Moreover, protein-protein association network analysis showed that the CmBGAL12 is the only node protein. Furthermore, spatiotemporal expression pattern analysis by quantitative real-time PCR (qRT-PCR) suggested that most of CmBGALs expressed in tissues with vigorous cell wall remodeling/disassembly. In addition, cis-acting regulatory elements analysis in promoters inferred that CmBGALs might participate in diverse responsiveness to phytohormone, biotic and abiotic signaling.

CONCLUSIONS

A novel clade of CmBGAL members (Clade F) related to melon fruit softening was discovered, since their expression showed a specific surge in the mature fruit of 'HPM' with mealy texture (softening sharply), but not in 'HDB' with crisp texture (softening bluntly). The homologous CmBGAL7-11 in Clade F exhibited identical spatiotemporal expression patterns may multiple genes leading to melon fruit softening.

Genome-wide identification and expression profiling of WUSCHEL-related homeobox (WOX) genes confer their roles in somatic embryogenesis, growth and abiotic stresses in banana

Plant-specific WUSCHEL-related homeobox (WOX) transcription factors are known to be involved in plant developmental processes, especially in embryogenesis. In this study, a total of thirteen WOX members were identified in the banana (Musa acuminata) genome (MaWOX) and characterized for in-silico analysis. Phylogenetic analysis revealed that these genes were divided into three clades (ancient, intermediate and modern) which reflected the evolutionary history of WOX families. Furthermore, modern clade members have shown higher variations in gene structural features and carried unique conserved motifs (motif 3 and motif 4) when compared to the members of other clades. The differential expression of all 13 MaWOX was observed in early (embryogenic cell suspension (ECS), multiplying ECS, germinating embryos, young leaflet and node of germinated plantlets) and late (unripe fruit peel and pulp, ripe fruit peel and pulp) developmental stages of banana cultivar Grand Naine. The maximum expression of MaWOX6 (18 fold) and MaWOX13 (120 fold) was found during somatic embryogenesis and in unripe fruit pulp, respectively. Moreover, numerous cis-elements responsive to drought, cold, ethylene, methyl jasmonate (MeJA), abscisic acid (ABA) and gibberellic acid (GA) were observed in all MaWOX promoter regions. The subsequent expression analysis under various abiotic stresses (cold, drought and salt) revealed maximum expression of the MaWOX3 (830 fold), MaWOX8a (30 fold) and MaWOX11b (105 fold) in salt stress. It gives evidence about their possible role in salt stress tolerance in banana. Hence, the present study provides precise information on the MaWOX gene family and their expression in various tissues and stressful environmental conditions that may help to develop climate-resilient banana plants.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-022-03387-w.

Identification and characteristics of SnRK genes and cold stress-induced expression profiles in Liriodendron chinense

Background

The sucrose non-fermenting 1 (SNF1)-related protein kinases (SnRKs) play a vivid role in regulating plant metabolism and stress response, providing a pathway for regulation between metabolism and stress signals. Conducting identification and stress response studies on SnRKs in plants contributes to the development of strategies for tree species that are more tolerant to stress conditions.

Results

In the present study, a total of 30 LcSnRKs were identified in Liriodendron chinense (L. chinense) genome, which was distributed across 15 chromosomes and 4 scaffolds. It could be divided into three subfamilies: SnRK1, SnRK2, and SnRK3 based on phylogenetic analysis and domain types. The LcSnRK of the three subfamilies shared the same Ser/Thr kinase structure in gene structure and motif composition, while the functional domains, except for the kinase domain, showed significant differences. A total of 13 collinear gene pairs were detected in L. chinense and Arabidopsis thaliana (A. thaliana), and 18 pairs were detected in L. chinense and rice, suggesting that the LcSnRK family genes may be evolutionarily more closely related to rice. Cis-regulation element analysis showed that LcSnRKs were LTR and TC-rich, which could respond to different environmental stresses. Furthermore, the expression patterns of LcSnRKs are different at different times under low-temperature stress. LcSnRK1s expression tended to be down-regulated under low-temperature stress. The expression of LcSnRK2s tended to be up-regulated under low-temperature stress. The expression trend of LcSnRK3s under low-temperature stress was mainly up-or down-regulated.

Conclusion

The results of this study will provide valuable information for the functional identification of the LcSnRK gene in the future.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08902-0.

Genome-Wide Identification and Expression Profile Analysis of the SnRK2 Gene Family in Nicotiana tabacum

SnRK2 protein kinase family plays an important role in plant response to abiotic stress and has been identified in various plants. This study aimed to identify SnRK2 genes in tobacco and systematically analyze their expression under abscisic acid treatment and abiotic stress. We identified 22 NtSnRK2 members, which were divided into three groups and located on 13 chromosomes, mainly at both ends of the chromosomes; additionally, 11 duplicated NtSnRK2 gene pairs were observed. Phylogenetic analysis showed that these SnRK2 members were divided into three groups in tobacco. The motifs of NtSnRK2 proteins in the same group were highly similar. Subcellular localization indicated that NtSnRK2s in Group3 were present in the nucleus, cytomembrane, and cytoplasm. Gene expression pattern analysis revealed that NtSnRK2 genes played a role in the responses to several abiotic stresses (salt, drought, and low-temperature stress), indicating that they are widely involved in the adaptation of tobacco to adverse environmental conditions.

Unraveling the genetics of polyamine metabolism in barley for senescence-related crop improvement

We explored the polyamine (PA) metabolic pathway genes in barley (Hv) to understand plant development and stress adaptation in Gramineae crops with emphasis on leaf senescence. Bioinformatics and functional genomics tools were utilized for genome-wide identification, comprehensive gene features, evolution, development and stress effects on the expression of the polyamine metabolic pathway gene families (PMGs). Three S-adenosylmethionine decarboxylases (HvSAMDCs), two ornithine decarboxylases (HvODCs), one arginine decarboxylase (HvADC), one spermidine synthase (HvSPDS), two spermine synthases (HvSPMSs), five copper amine oxidases (HvCuAOs) and seven polyamine oxidases (HvPAOs) members of PMGs were identified and characterized in barley. All the HvPMG genes were found to be distributed on all chromosomes of barley. The phylogenetic and comparative assessment revealed that PA metabolic pathway is highly conserved in plants and the prediction of nine H. vulgare miRNAs (hvu-miR) target sites, 18 protein-protein interactions and 961 putative CREs in the promoter region were discerned. Gene expression of HvSAMDC3, HvCuAO7, HvPAO4 and HvSPMS1 was apparent at every developmental stage. SPDS/SPMS gene family was found to be the most responsive to induced leaf senescence. This study provides a reference for the functional investigation of the molecular mechanism(s) that regulate polyamine metabolism in plants as a tool for future breeding decision management systems.

Genome-wide analysis of the DREB family genes and functional identification of the involvement of BrDREB2B in abiotic stress in wucai (Brassica campestris L.)

Dehydration responsive element binding protein (DREB) is a significant transcription factor class known to be implicated in abiotic stresses. In this study, we systematically conducted a genome-wide identification and expression analysis of the DREB gene family, including gene structures, evolutionary relationships, chromosome distribution, conserved domains, and expression patterns. A total of 65 DREB family gene members were identified in Chinese cabbage (Brassica rapa L.) and were classified into five subgroups based on phylogenetic analysis. Through analysis of the conserved domains of BrDREB family genes, only one exon existed in the gene structure. Through the analysis of cis-acting elements, these genes were mainly involved in hormone regulation and adversity stress. In order to identify the function of BrDREB2B, overexpressed transgenic Arabidopsis was constructed. After different stress treatments, the germination rate, root growth, survival rate, and various plant physiological indicators were measured. The results showed that transgenic Arabidopsis thaliana plants overexpressing BrDREB2B exhibited enhanced tolerance to salt, heat and drought stresses. Taken together, our results are the first to report the BrDREB2B gene response to drought and heat stresses in Chinese cabbage and provide a basis for further studies to determine the function of BrDREBs in response to abiotic stresses.

Different responses of banana classical AGP genes and cell wall AGP components to low-temperature between chilling sensitive and tolerant cultivars

Seventeen classical MaAGPs and 9 MbAGPs were identified and analyzed. MaAGP1/2/6/9/16/17, the antigens of JIM13 and LM2 antibodies are likely to be involved in banana chilling tolerance. Classical arabinogalactan proteins (AGPs) belong to glycosylphosphatidylinositol-anchored proteins, which are proved to be involved in signaling and cell wall metabolism upon stresses. However, rare information is available on the roles of classical AGPs in low temperature (LT) tolerance. Cultivation of banana in tropical and subtropical region is seriously threatened by LT stress. In the present study, 17 classical MaAGPs and nine MbAGPs in banana A and B genome were identified and characterized, respectively. Great diversity was present among different classical MaAGP/MbAGP members while five members (AGP3/6/11/13/14) showed 100% identity between these two gene families. We further investigated different responses of classical AGPs to LT between a chilling sensitive (CS) and tolerant (CT) banana cultivars. In addition, different changes in the temporal and spatial distribution of cell wall AGP components under LTs between these two cultivars were compared using immunofluorescence labeling. Seven classical MbAGPs were upregulated by LT(s) in the CT cultivar. Classical MaAGP4/6 was induced by LT(s) in both cultivars while MaAGP1/2/9/16/17 only in the CT cultivar. Moreover, these genes showed significantly higher transcription abundance in the CT cultivar than the CS one under LT(s) except classical MaAGP4. Similar results were observed with the epitopes of JIM13 and LM2 antibodies. The antigens of these antibodies and classical MaAGP1/2/6/9/16/17 might be related to LT tolerance of banana. These results provide additional information about plant classical AGPs and their involvement in LT tolerance, as well as their potential as candidate genes to be targeted when breeding CT banana.

Genome-wide analysis and functional characterization of CHYR gene family associated with abiotic stress tolerance in bread wheat (Triticum aestivum L.)

BACKGROUND

CHY zinc-finger and RING finger (CHYR) proteins have been functionally characterized in plant growth, development and various stress responses. However, the genome-wide analysis was not performed in wheat.

RESULTS

In this study, a total of 18 TaCHYR genes were identified in wheat and classified into three groups. All TaCHYR genes contained CHY-zinc finger, C3H2C3-type RING finger and zinc ribbon domains, and group III members included 1-3 hemerythrin domains in the N-terminus regions. TaCHYR genes in each group shared similar conserved domains distribution. Chromosomal location, synteny and cis-elements analysis of TaCHYRs were also analyzed. Real-time PCR results indicated that most of selected 9 TaCHYR genes exhibited higher expression levels in leaves during wheat seedling stage. All these TaCHYR genes were up-regulated after PEG treatment, and these TaCHYRs exhibited differential expression patterns in response to salt, cold and heat stress in seedling leaves. The growth of yeast cells expressing TaCHYR2.1, TaCHYR9.2 and TaCHYR11.1 were inhibited under salt and dehydration stress. Moreover, gene ontology (GO) annotation, protein interaction and miRNA regulatory network of TaCHYR genes were analyzed.

CONCLUSIONS

These results increase our understanding of CHYR genes and provide robust candidate genes for further functional investigations aimed at crop improvement.

Identification and characterization of the BEL1-like genes reveal their potential roles in plant growth and abiotic stress response in tomato

BEL1-like (BELL) transcription factors, belonging to three-amino acid-loop-extension (TALE) superfamily, are ubiquitous in plants. BELLs regulate a wide range of plant biological processes, but the understanding of the BELL family in tomato (Solanum lycopersicum) remains fragmentary. In this study, a total of 14 members of the SlBELL family were identified in tomato. SlBELL proteins contained the conserved BELL and SKY domains that served as typical structures of the BELL family. Syntenic analysis indicated that the BELL orthologs between tomato and other dicots had close evolutionary relationships. Furthermore, the promoters of SlBELLs contained numerous cis-elements related to plant growth, development, and stress response. The SlBELL genes exhibited different tissue-specific expression profiles and responded to cold, heat, and drought stresses, implying their potential functions in regulating multiple aspects of plant growth, as well as in response to abiotic stresses. Through the interaction network prediction, we found that most SlBELL proteins displayed probable interactions with the KNOTTED1-like (KNOX) proteins, another kind of transcription factor in the TALE superfamily. These findings laid foundations for further dissection of the functions of SlBELL genes in tomato, as well as for exploration of the evolutionary relationships of BELL homologs among different plant species.

Genome-wide identification and expression profiling of thes MAPK, MAPKK, and MAPKKK gene families in Ophiocordyceps sinensis

Mitogen-activated protein kinase (MAPK) cascades have a universal cell signaling mechanism in eukaryotes. A typical MAPK signal transduction module comprises three kinds of sequentially phosphorylated protein kinases: MAPK, Mitogen-activated protein kinase kinase (MAPKK), and Mitogen-activated protein kinase kinase kinase (MAPKKK). However, little is known regarding the genes involved in MAPK cascades in Ophiocordyceps sinensis. Nine genes (three MAPK, three MAPKK, and three MAPKKK) were identified in this study. The MAPK, MAPKK, and MAPKKK genes were divided into three subfamilies, according to the phylogenetic analysis. TEY and TGY represented the activation domains of the MAPKs; the corresponding domains in MAPKKs were SDIWS and SDVWS, and those in the MAPKKs were GSVFYWMAPEV and GTPMYMSPEV. Transcription data analysis and quantitative real-time polymerase chain reaction showed that the MAPK cascade was related to the growth of the fruiting body. This is the first study to report a genome-wide identification of the MAPK, MAPKK, and MAPKK gene families in O. sinensis.

Genome-wide characterization on MT family and their expression in response to environmental cues in upland cotton (Gossypium hirsutum L.)

Metallothioneins (MTs) are believed as key metal chelators and scavengers of reactive oxygen species (ROS), which are involved in tolerance and de-toxicity to multiple environmental stresses in plants. The MT gene family was characterized from upland cotton (Gossypium hirsutum L.), compared with its putative genome donors G. arboretum and raimondii. Subsequently, gene functions were predicted by promoter analysis. Moreover, gene expressions subjecting to exogenous stimuli, as well as in terms of developments, were studied. The main findings were shown as follows: 1) 19 GhMTs were identified from G. hirsutum, and the family completely included all four sub-types, namely p1, p2, p3, and pec. Sub-type p2 GhMTs were most conservative in protein motif compositions, gene structures, phylogenic relationships, and group numbers, while p3 GhMTs demonstrated much more diversiform and distant genetic relationships. 2) The GhMT family experienced apparent gene expansion, and the members from the D sub-genome were subjected to stronger environmental selection. 3) GhMTs played differential and overlapped roles in response to environmental cues. 4) GhMT6, GhMT8, and GhMT14 were involved in both vegetative and reproductive developments. These findings must provide valuable insights into understanding the plant MT gene family and novel gene resources for cotton breeding for environmental stresses, phytoremediation, and beyond.

Genome-wide identification of PIP5K in wheat and its relationship with anther male sterility induced by high temperature

BACKGROUND

Phosphatidylinositol 4 phosphate 5-kinase (PIP5K) plays a key enzyme role in the inositol signal transduction system and has essential functions in plants in terms of growth, development, and stress responses. However, systematic studies on the wheat PIP5K gene family and its relation to male sterility have not been reported yet.

RESULTS

Sixty-four TaPIP5K genes were identified. The TaPIP5K genes contained similar gene structures and conserved motifs on the same branches of the evolutionary tree, and their cis-regulatory elements were related to MeJA-responsiveness. Furthermore, 49 pairs of collinearity genes were identified and mainly subjected to purification selection during evolution. Synteny analyses showed that some PIP5K genes in wheat and the other four species shared a relatively conserved evolutionary process. The expression levels of many conservative TaPIP5K genes in HT-ms anthers were significantly lower than that in Normal anthers. In addition, HT-ms anthers have no dehiscence, and levels of OPDA and JA-ILE are significantly lower at the trinucleus stage.

CONCLUSION

These results indicate that the PIP5K gene family may be associated with male sterility induced by HT, and the reduction of JA-ILE levels and the abnormal levels of these genes expression may be one reason for the HT-ms anthers having no dehiscence, ultimately leading to the abortion of the anthers.

Genome-wide identification of PME genes, evolution and expression analyses in soybean (Glycine max L.)

BACKGROUND

Pectin methylesterase (PME) is one of pectin-modifying enzyme that affects the pectin homeostasis in cell wall and regulates plant growth and diverse biological processes. The PME genes have been well explored and characterized in different plants. Nevertheless, systematic research on the soybean (Glycine max L.) PME genes remain lacking.

RESULTS

We identified 127 Glycine max PME genes (GmPME) from the soybean Wm82.a2.v1 genome, which unevenly distributed on 20 soybean chromosomes. Phylogenetic analysis classified the GmPME genes into four clades (Group I, Group II, Group III and Group IV). GmPME gene members in the same clades displayed similar gene structures and motif patterns. The gene family expansion analysis demonstrated that segmental duplication was the major driving force to acquire novel GmPME genes compared to the tandem duplication events. Further synteny and evolution analyses showed that the GmPME gene family experienced strong purifying selective pressures during evolution. The cis-element analyses together with the expression patterns of the GmPME genes in various tissues suggested that the GmPME genes broadly participate in distinct biological processes and regulate soybean developments. Importantly, based on the transcriptome data and quantitative RT-PCR validations, we examined the potential roles of the GmPME genes in regulating soybean flower bud development and seed germination.

CONCLUSION

In conclusion, we provided a comprehensive characterization of the PME genes in soybean, and our work laid a foundation for the functional study of GmPME genes in the future.

Comparative genomic analysis of expansin superfamily gene members in zucchini and cucumber and their expression profiles under different abiotic stresses

Zucchini and cucumber belong to the Cucurbitaceae family, a group of economical and nutritious food plants that is consumed worldwide. Expansin superfamily proteins are generally localized in the cell wall of plants and are known to possess an effect on cell wall modification by causing the expansion of this region. Although the whole genome sequences of cucumber and zucchini plants have been resolved, the determination and characterization of expansin superfamily members in these plants using whole genomic data have not been implemented yet. In the current study, a genome-wide analysis of zucchini (Cucurbita pepo) and cucumber (Cucumis sativus) genomes was performed to determine the expansin superfamily genes. In total, 49 and 41 expansin genes were identified in zucchini and cucumber genomes, respectively. All expansin superfamily members were subjected to further bioinformatics analysis including gene and protein structure, ontology of the proteins, phylogenetic relations and conserved motifs, orthologous relations with other plants, targeting miRNAs of those genes and in silico gene expression profiles. In addition, various abiotic stress responses of zucchini and cucumber expansin genes were examined to determine their roles in stress tolerance. CsEXPB-04 and CsEXPA-11 from cucumber and CpEXPA-20 and CpEXPLA-14 from zucchini can be candidate genes for abiotic stress response and tolerance in addition to their roles in the normal developmental processes, which are supported by the gene expression analysis. This work can provide new perspectives for the roles of expansin superfamily genes and offers comprehensive knowledge for future studies investigating the modes of action of expansin proteins.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s12298-021-01108-w.

Comprehensive genome wide identification and expression analysis of MTP gene family in tomato (Solanum lycopersicum) under multiple heavy metal stress

Plant metal tolerance proteins (MTPs) play major roles in enhancing resistance to heavy metal tolerance and homeostasis. However, the role of MTPs genes in tomato, which is one of the most popular crops, is still largely limited. Hence, we investigated genome-wide study of tomato MTPs, including phylogenetic, duplication, gene structure, gene ontology and previous transcriptomic data analysis. Moreover, the MTPs expression behaviour under various heavy metals stress has rarely been investigated. In the current study, eleven MTP candidate genes were genome-wide identified and classified into three major groups; Mn-cation diffusion facilitators (CDFs), Fe/Zn-CDFs, and Zn-CDFs based on the phylogeny. Structural analysis of SlMTPs showed high gene similarity within the same group with cation_efflux or ZT_dimerdomains. Evolutionary analysis revealed that segmental duplication contributed to the expansion of the SlMTP family. Gene ontology further showed the vital roles of MTPs in metal-related processes. Tissue-specific expression profiling exhibited similar expression patterns in the same group, whereas gene expression varied among groups. The MTPs expression was evaluated after tomato treatments by five divalent heavy metals (Cd2+, Co2+, Mn2+, Zn2+, and Fe2+). SlMTP genes displayed differential responses in either plant leaves or roots under heavy metals treatments. Nine and ten SlMTPs responded to at least one metal ion treatment in leaves and roots, respectively. In addition SlMTP1, SlMTP3, SlMTP4, SlMTP8, SlMTP10 and SlMTP11 exhibited the highest expression responses in most of heavy metals treatments. Overall, our findings presented a standpoint on the evolution of MTPs and their evolution in tomato and paved the way for additional functional characterization under heavy metal toxicity.

Genome-wide identification and expression profiling of invertase gene family for abiotic stresses tolerance in Poncirus trifoliata

BACKGROUND

Sucrose (Suc) hydrolysis is directly associated with plants tolerance to multiple abiotic stresses. Invertase (INV) enzymes irreversibly catalyze Suc degradation to produce glucose (Glc) and fructose (Frc). However, genome-wide identification and function of individual members of the INV gene family in Poncirus trifoliata or its Citrus relatives in response to abiotic stresses are not fully understood.

RESULTS

In this report, fourteen non-redundant PtrINV family members were identified in P. trifoliata including seven alkaline/neutral INV genes (PtrA/NINV1-7), two vacuolar INV genes (PtrVINV1-2), and five cell wall INV isoforms (PtrCWINV1-5). A comprehensive analysis based on the biochemical characteristics, the chromosomal location, the exon-intron structures and the evolutionary relationships demonstrated the conservation and the divergence of PtrINVs. In addition, expression analysis of INV genes during several abiotic stresses in various tissues indicated the central role of A/NINV7 among INV family members in response to abiotic stresses. Furthermore, our data demonstrated that high accumulation of Suc, Glc, Frc and total sugar contents were directly correlated with the elevated activities of soluble INV enzymes in the cold-tolerant P. trifoliata, C. ichangensis and C. sinensis, demonstrating the potential role of soluble INV enzymes for the cold tolerance of Citrus.

CONCLUSIONS

This work offered a framework for understanding the physiological role of INV genes and laid a foundation for future functional studies of these genes in response to abiotic stresses.

Comparative genomic analysis of the tricarboxylic acid cycle members in four Solanaceae vegetable crops and expression pattern analysis in Solanum tuberosum

Background

The tricarboxylic acid (TCA) cycle is crucial for energy supply in animal, plant, and microbial cells. It is not only the main pathway of carbohydrate catabolism but also the final pathway of lipid and protein catabolism. Some TCA genes have been found to play important roles in the growth and development of tomato and potato, but no comprehensive study of TCA cycle genes in Solanaceae crops has been reported.

Results

In this study, we analyzed TCA cycle genes in four important Solanaceae vegetable crops (potato (Solanum tuberosum), tomato (Solanum lycopersicum), eggplant (Solanum melongena), and pepper (Capsicum annuum)) based on comparative genomics. The four Solanaceae crops had a total of 180 TCA cycle genes: 43 in potato, 44 in tomato, 40 in eggplant, and 53 in pepper. Phylogenetic analysis, collinearity analysis, and tissue expression patterns revealed the conservation of and differences in TCA cycle genes between the four Solanaceae crops and found that there were unique subgroup members in Solanaceae crops that were independent of Arabidopsis genes. The expression analysis of potato TCA cycle genes showed that (1) they were widely expressed in various tissues, and some transcripts like Soltu.DM.01G003320.1(SCoAL) and Soltu.DM.04G021520.1 (SDH) mainly accumulate in vegetative organs, and some transcripts such as Soltu.DM.12G005620.3 (SDH) and Soltu.DM.02G007400.4 (MDH) are preferentially expressed in reproductive organs; (2) several transcripts can be significantly induced by hormones, such as Soltu.DM.08G023870.2 (IDH) and Soltu.DM.06G029290.1 (SDH) under ABA treatment, and Soltu.DM.07G021850.2 (CSY) and Soltu.DM.09G026740.1 (MDH) under BAP treatment, and Soltu.DM.02G000940.1 (IDH) and Soltu.DM.01G031350.4 (MDH) under GA treatment; (3) Soltu.DM.11G024650.1 (SDH) can be upregulated by the three disease resistance inducers including Phytophthora infestans, acibenzolar-S-methyl (BTH), and DL-β-amino-n-butyric acid (BABA); and (4) the levels of Soltu.DM.01G045790.1 (MDH), Soltu.DM.01G028520.3 (CSY), and Soltu.DM.12G028700.1 (CSY) can be activated by both NaCl and mannitol. The subcellular localization results of three potato citrate synthases showed that Soltu.DM.01G028520.3 was localized in mitochondria, while Soltu.DM.12G028700.1 and Soltu.DM.07G021850.1 were localized in the cytoplasm.

Conclusions

This study provides a scientific foundation for the comprehensive understanding and functional studies of TCA cycle genes in Solanaceae crops and reveals their potential roles in potato growth, development, and stress response.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-08109-9.

Characterization of NAC family genes in Salvia miltiorrhiza and NAC2 potentially involved in the biosynthesis of tanshinones

The NAC (NAM, ATAF, and CUC) family members are specific transcription factors in plants. The large family is involved in many plant growth and developmental processes, as well as in abiotic/biotic stress responses. It has been well studied in the genomes of various plants, including Arabidopsis thaliana, tomato, and quinoa. However, identification and functional studies of NAC family members in medicinal Salvia miltiorrhiza are limited. Here, we systematically identified 84 NAC genes and named them according to their gene IDs in the recently sequenced genome. The phylogeny of NAC family protein sequences was analyzed using bioinformatics methods, which divided them into nine subfamilies. Then, their chromosomal locations, gene structures and conserved domains were analyzed comprehensively. To further investigate the regulatory functions of NACs in S. miltiorrhiza, we analyzed the response of 10 selected NAC genes to methyl jasmonate and used NAC2 for transgenic experiments. The overexpression of Sm-NAC2 decreased the tanshinone I and IIA contents by 56% and 62%, respectively. However, Sm-NAC2-RNAi promoted the accumulation of four tanshinones, tanshinone I, tanshinone IIA, cryptotanshinone, and dihydrotanshinone I, which increased 3.68-, 4.1-, 3.13- and 5.9- fold, respectively, compared with wild type. In the tanshinone biosynthetic pathways, the overexpression of Sm-NAC2 down-regulated CYP76AH1, and the silencing of Sm-NAC2 up-regulated the expression levels of HMGR1, DXS2, KSL2, and CYP76AH1. This study provides information on the evolution of Sm-NAC genes and their possible functions, and it lays a foundation for further research into the NAC family-associated regulation of tanshinone biosynthesis.

Comparative genomic analysis of MYB transcription factors for cuticular wax biosynthesis and drought stress tolerance in Helianthus annuus L

Sunflower is an important oil-seed crop in Pakistan, it is mainly cultivated in the spring season. It is severely affected by drought stress resulting in lower yield. Cuticular wax acts as the first defense line to protect plants from drought stress condition. It seals the aerial parts of plants and reduce the water loss from leaf surfaces. Various myeloblastosis (MYB) transcription factors (TFs) are involved in biosynthesis of epicuticular waxes under drought-stress. However, less information is available for MYB, TFs in drought stress and wax biosynthesis in sunflower. We used different computational tools to compare the Arabidopsis MYB, TFs involved in cuticular wax biosynthesis and drought stress tolerance with sunflower genome. We identified three putative MYB genes (MYB16, MYB94 and MYB96) in sunflower along with their seven homologs in Arabidopsis. Phylogenetic association of MYB TFs in Arabidopsis and sunflower indicated strong conservation of TFs in plant species. From gene structure analysis, it was observed that intron and exon organization was family-specific. MYB TFs were unevenly distributed on sunflower chromosomes. Evolutionary analysis indicated the segmental duplication of the MYB gene family in sunflower. Quantitative Real-Time PCR revealed the up-regulation of three MYB genes under drought stress. The gene expression of MYB16, MYB94 and MYB96 were found many folds higher in experimental plants than control. The present study provided the first insight into MYB TFs family's characterization in sunflower under drought stress conditions and wax biosynthesis TFs.