Genome-wide analysis of C2H2 zinc finger family and their response to abiotic stresses in apple

C2H2-type zinc finger proteins are one of the most widely studied families in plants and play important roles in abiotic stress responses. In the present study, the physicochemical properties, chromosomal locations, evolutionary relationships, and gene structures of 54 C2H2 zinc finger protein (ZFP) family members were analyzed in apple. The MdC2H2-ZFP genes were phylogenetically clustered into seven subfamilies distributed in different densities on 16 chromosomes. The RNA-seq data from various tissues revealed that MdC2H2-ZFPs differentially expressed among root, stem, leaf, flower, and fruits. Quantitative analysis of its expression characteristics showed that the MdC2H2-ZFP genes were rapidly induced as exposure to abiotic stresses such as drought, salt and low temperature etc. Under drought stress, the expression of eight members was significantly up-regulated, and the highest was obtained from MdC2H2-17; as exposure to salt stress, nine MdC2H2-ZFPs was obviously up-regulated, with the highest expression of MdC2H2-13; and under low temperature stress, the expression of seven members was highly up-regulated, and MdC2H2-13 also demonstrated the highest expression which is same as the case under salt stress. Therefore, some members of MdC2H2-ZFP gene family considerably involve in the multiple abiotic stress responses, which may better understand the function of this family and facilitate the breeding of apple for stress tolerance.

Characterization and expression profiles of the ZmLBD gene family in Zea mays

BACKGROUND

The Lateral Organ Boundaries Domain (LBD) gene family is a family of plant-specific transcription factors (TFs) that are widely involved in processes such as lateral organ formation, stress response, and nutrient metabolism. However, the function of LBD genes in maize remains poorly understood.

METHODS AND RESULTS

In this study, a total of 49 ZmLBD genes were identified at the genome-wide level of maize, they were classified into nine branches based on phylogenetic relationships, and all of them were predicted to be nuclear localized. The 49 ZmLBD genes formed eight pairs of segmental duplicates, and members of the same branches' members had similar gene structure and conserved motif composition. The promoters of ZmLBD genes contain multiple types of cis-acting elements. In addition, by constructing the regulatory network of ZmLBD and other genes and miRNAs, 12 and 22 ZmLBDs were found to be involved in the gene regulatory network and miRNA regulatory network, respectively. The expression pattern analysis suggests that ZmLBD genes may be involved in different biological pathways, and drought stress induced the expressions of two inbred lines.

CONCLUSIONS

The findings enhance our comprehension of the potential roles of the ZmLBD gene family in maize growth and development, which is pivotal for genetic enhancement and breeding efforts pertaining to this significant crop.

Origin, evolution, and diversification of inositol 1,4,5-trisphosphate 3-kinases in plants and animals

BACKGROUND

In Eukaryotes, inositol polyphosphates (InsPs) represent a large family of secondary messengers and play crucial roes in various cellular processes. InsPs are synthesized through a series of pohophorylation reactions catalyzed by various InsP kinases in a sequential manner. Inositol 1,4,5-trisphosphate 3-kinase (IP3 3-kinase/IP3K), one member of InsP kinase, plays important regulation roles in InsPs metabolism by specifically phosphorylating inositol 1,4,5-trisphosphate (IP3) to inositol 1,3,4,5-tetrakisphosphate (IP4) in animal cells. IP3Ks were widespread in fungi, plants and animals. However, its evolutionary history and patterns have not been examined systematically.

RESULTS

A total of 104 and 31 IP3K orthologues were identified across 57 plant genomes and 13 animal genomes, respectively. Phylogenetic analyses indicate that IP3K originated in the common ancestor before the divergence of fungi, plants and animals. In most plants and animals, IP3K maintained low-copy numbers suggesting functional conservation during plant and animal evolution. In Brassicaceae and vertebrate, IP3K underwent one and two duplication events, respectively, resulting in multiple gene copies. Whole-genome duplication (WGD) was the main mechanism for IP3K duplications, and the IP3K duplicates have experienced functional divergence. Finally, a hypothetical evolutionary model for the IP3K proteins is proposed based on phylogenetic theory.

CONCLUSION

Our study reveals the evolutionary history of IP3K proteins and guides the future functions of animal, plant, and fungal IP3K proteins.

Genome-wide investigation of UDP-Glycosyltransferase family in Tartary buckwheat (Fagopyrum tataricum)

BACKGROUND

Tartary buckwheat (Fagopyrum tataricum) belongs to Polygonaceae family and has attracted increasing attention owing to its high nutritional value. UDP-glycosyltransferases (UGTs) glycosylate a variety of plant secondary metabolites to control many metabolic processes during plant growth and development. However, there have been no systematic reports of UGT superfamily in F. tataricum.

RESULTS

We identified 173 FtUGTs in F. tataricum based on their conserved UDPGT domain. Phylogenetic analysis of FtUGTs with 73 Arabidopsis UGTs clustered them into 21 families. FtUGTs from the same family usually had similar gene structure and motif compositions. Most of FtUGTs did not contain introns or had only one intron. Tandem repeats contributed more to FtUGTs amplification than segmental duplications. Expression analysis indicates that FtUGTs are widely expressed in various tissues and likely play important roles in plant growth and development. The gene expression analysis response to different abiotic stresses showed that some FtUGTs were involved in response to drought and cadmium stress. Our study provides useful information on the UGTs in F. tataricum, and will facilitate their further study to better understand their function.

CONCLUSIONS

Our results provide a theoretical basis for further exploration of the functional characteristics of FtUGTs and for understanding the growth, development, and metabolic model in F. tataricum.

Identification of WRKY gene family in Dioscorea opposita Thunb. reveals that DoWRKY71 enhanced the tolerance to cold and ABA stress

WRKY transcription factors constitute one of the largest plant-specific gene families, regulating various aspects of plant growth, development, physiological processes, and responses to abiotic stresses. This study aimed to comprehensively analyze the WRKY gene family of yam (Dioscorea opposita Thunb.), to understand their expression patterns during the growth and development process and their response to different treatments of yam and analyze the function of DoWRKY71 in detail. A total of 25 DoWRKY genes were identified from the transcriptome of yam, which were divided into six clades (I, IIa, IIc, IId, IIe, III) based on phylogenetic analysis. The analysis of conserved motifs revealed 10 motifs, varying in length from 16 to 50 amino acids. Based on real-time quantitative PCR (qRT-PCR) analysis, DoWRKY genes were expressed at different stages of growth and development and responded differentially to various abiotic stresses. The expression level of DoWRKY71 genes was up-regulated in the early stage and then down-regulated in tuber enlargement. This gene showed responsiveness to cold and abiotic stresses, such as abscisic acid (ABA) and methyl jasmonate (MeJA). Therefore, further study was conducted on this gene. Subcellular localization analysis revealed that the DoWRKY71 protein was localized in the nucleus. Moreover, the overexpression of DoWRKY71 enhanced the cold tolerance of transgenic tobacco and promoted ABA mediated stomatal closure. This study presents the first systematic analysis of the WRKY gene family in yam, offering new insights for studying WRKY transcription factors in yam. The functional study of DoWRKY71 lays theoretical foundation for further exploring the regulatory function of the DoWRKY71 gene in the growth and development related signaling pathway of yam.

Identification of the NOD-like receptor family of golden pompano and expression in response to bacterial and parasitic exposure reveal its key role in innate immunity

This study presents a genome-wide identification of NOD-like receptors (NLRs) in the golden pompano, key to its innate immunity. We identified 30 ToNLRs, analyzing their chromosomal positions, characteristics, evolutionary relationships, evidence of positive selection, and synteny with the yellowtail kingfish. Our findings categorize these NLRs into three main subgroups: NLRA, NLRC, and the distinct ToNLRX1. Post-exposure to Streptococcus agalactiae, most ToNLRs increased expression in the spleen, whereas NLRC3like13, NLRC3like16, and NLRC3like19 so in the kidneys. Upon Cryptocaryon irritans exposure, we categorized our groups based on the site of infection into the control group (BFS), the trophont-attached skin (TAS), and the nearby region skin (NRS). ToAPAF1 and ToNOD1 expressions rose in the NRS, in contrast to decreased expressions of ToNLRC5, ToNWD1 and ToCIITA. Other ToNLRs showed variable expressions in the TAS. Overall, this research lays the groundwork for further exploration of innate immunity in the golden pompano.

Genome-wide identification of the N6-methyladenosine regulatory genes reveals NtFIP37B increases drought resistance of tobacco (Nicotiana tabacum L.)

BACKGROUND

N6-methyladenosine (m6A) is one of the common internal RNA modifications found in eukaryotes. The m6A modification can regulate various biological processes in organisms through the modulation of alternative splicing, alternative polyadenylation, folding, translation, localization, transport, and decay of multiple types of RNA, without altering the nucleotide sequence. The three components involved in m6A modification, namely writer, eraser, and reader, mediate the abundance of RNA m6A modification through complex collaborative actions. Currently, research on m6A regulatory genes in plants is still in its infancy.

RESULTS

In this study, we identified 52 candidate m6A regulatory genes in common tobacco (Nicotiana tabacum L.). Gene structure, conserved domains, and motif analysis showed structural and functional diversity among different subgroups of tobacco m6A regulatory genes. The amplification of m6A regulatory genes were mainly driven by polyploidization and dispersed duplication, and duplicated genes evolved through purified selection. Based on the potential regulatory network and expression pattern analysis of m6A regulatory genes, a significant number of m6A regulatory genes might play important roles in growth, development, and stress response processes. Furthermore, we have confirmed the critical role of NtFIP37B, an m6A writer gene in tobacco, in enhancing drought resistance.

CONCLUSIONS

This study provides useful information for better understanding the evolution of m6A regulatory genes and the role of m6A modification in tobacco stress response, and lays the foundation for further elucidating the function of m6A regulatory genes in tobacco.

Regulation of blue infertile flower pigmentation by WD40 transcription factor HmWDR68 in Hydrangea macrophylla 'forever summer'

BACKGROUND

WD40 transcription factors are crucial in plant growth and developmental, significantly impacting plant growth regulation. This study investigates the WD40 transcription factor HmWDR68's role in developing the distinctive blue infertile flower colors in Hydrangea macrophylla 'Forever Summer'.

METHODS AND RESULTS

The HmWDR68 gene was isolated by PCR, revealing an open reading frame of 1026 base pairs, which encodes 341 amino acids. Characterized by four WD40 motifs, HmWDR68 is a member of the WD40 family. Phylogenetic analysis indicates that HmWDR68 shares high homology with PsWD40 in Camellia sinensis and CsWD40 in Paeonia suffruticosa, both of which are integral in anthocyanin synthesis regulation. Quantitative real-time PCR (qRT-PCR) analysis demonstrated that HmWDR68 expression in the blue infertile flowers of 'Forever Summer' hydrangea was significantly higher compared to other tissues and organs. Additionally, in various hydrangea varieties with differently colored infertile flowers, HmWDR68 expression was markedly elevated in comparison to other hydrangea varieties, correlating with the development of blue infertile flowers. Pearson correlation analysis revealed a significant association between HmWDR68 expression and the concentration of delphinidin 3-O-glucoside, as well as key genes involved in anthocyanin biosynthesis (HmF3H, HmC3'5'H, HmDFR, and HmANS) in the blue infertile flowers of 'Forever Summer' hydrangea (P < 0.01).

CONCLUSION

These findings suggest HmWDR68 may specifically regulate blue infertile flower formation in hydrangea by enhancing delphinidin-3-O-glucoside synthesis, modulating expression of HmF3H, HmC3'5'H, HmDFR and HmANS. This study provides insights into HmWDR68's role in hydrangea's blue flowers development, offering a foundation for further research in this field.

Genome-wide identification, stress- and hormone-responsive expression characteristics, and regulatory pattern analysis of Scutellaria baicalensis SbSPLs

SQUAMOSA PROMOTER BINDING PROTEIN-LIKEs (SPLs) encode plant-specific transcription factors that regulate plant growth and development, stress response, and metabolite accumulation. However, there is limited information on Scutellaria baicalensis SPLs. In this study, 14 SbSPLs were identified and divided into 8 groups based on phylogenetic relationships. SbSPLs in the same group had similar structures. Abscisic acid-responsive (ABRE) and MYB binding site (MBS) cis-acting elements were found in the promoters of 8 and 6 SbSPLs. Segmental duplications and transposable duplications were the main causes of SbSPL expansion. Expression analysis based on transcriptional profiling showed that SbSPL1, SbSPL10, and SbSPL13 were highly expressed in roots, stems, and flowers, respectively. Expression analysis based on quantitative real-time polymerase chain reaction (RT‒qPCR) showed that most SbSPLs responded to low temperature, drought, abscisic acid (ABA) and salicylic acid (SA), among which the expression levels of SbSPL7/9/10/12 were significantly upregulated in response to abiotic stress. These results indicate that SbSPLs are involved in the growth, development and stress response of S. baicalensis. In addition, 8 Sba-miR156/157 s were identified, and SbSPL1-5 was a potential target of Sba-miR156/157 s. The results of target gene prediction and coexpression analysis together indicated that SbSPLs may be involved in the regulation of L-phenylalanine (L-Phe), lignin and jasmonic acid (JA) biosynthesis. In summary, the identification and characterization of the SbSPL gene family lays the foundation for functional research and provides a reference for improved breeding of S. baicalensis stress resistance and quality traits.

Comprehensive bioinformation analysis of homeodomain-leucine zipper gene family and expression pattern of HD-Zip I under abiotic stress in Salix suchowensis

BACKGROUND

Homeodomain-leucine zipper (HD-Zip) transcription factors are plant-specific and play important roles in plant defense against environmental stresses. Identification and functional studies have been carried out in model plants such as rice, Arabidopsis thaliana, and poplar, but comprehensive analysis on the HD-Zip family of Salix suchowensis have not been reported.

RESULTS

A total of 55 HD-Zip genes were identified in the willow genome, unevenly distributed on 18 chromosomes except for chromosome 19. And segmental duplication events containing SsHD-Zip were detected on all chromosomes except chromosomes 13 and 19. The SsHD-Zip were classified into 4 subfamilies subfamilies (I-IV) according to the evolutionary analysis, and members of each subfamily shared similar domain structure and gene structure. The combination of GO annotation and promoter analysis showed that SsHD-Zip genes responded to multiple abiotic stresses. Furthermore, the results of qPCR analysis showed that the SsHD-Zip I gene exhibited different degrees of expression under salt stress, PEG treatment and heat treatment. Moreover, there was a synergistic effect between SsHD-Zip I genes under stress conditions based on coregulatory networks analysis.

CONCLUSIONS

In this study, HD-Zip transcription factors were systematically identified and analyzed at the whole genome level. These results preliminarily clarified the structural characteristics and related functions of willow HD-Zip family members, and it was found that SsHox34, SsHox36 and SsHox51 genes were significantly involved in the response to various stresses. Together, these findings laid the foundation for further research on the resistance functions of willow HD-Zip genes.

Genome-wide identification and analysis of WD40 proteins reveal that NtTTG1 enhances drought tolerance in tobacco (Nicotiana tabacum)

BACKGROUND

WD40 proteins, which are highly prevalent in eukaryotes, play important roles in plant development and stress responses. However, systematic identification and exploration of WD40 proteins in tobacco have not yet been conducted.

RESULTS

In this study, a total of 399 WD40 regulatory genes were identified in common tobacco (Nicotiana tabacum). Gene structure and motif analysis revealed structural and functional diversity among different clades of tobacco WD40 regulatory genes. The expansion of tobacco WD40 regulatory genes was mainly driven by segmental duplication and purifying selection. A potential regulatory network of NtWD40s suggested that NtWD40s might be regulated by miRNAs and transcription factors in various biological processes. Expression pattern analysis via transcriptome analysis and qRT-PCR revealed that many NtWD40s exhibited tissue-specific expression patterns and might be involved in various biotic and abiotic stresses. Furthermore, we have validated the critical role of NtTTG1, which was located in the nuclei of trichome cells, in enhancing the drought tolerance of tobacco plants.

CONCLUSIONS

Our study provides comprehensive information to better understand the evolution of WD40 regulatory genes and their roles in different stress responses in tobacco.

Clinicopathological Significance and Expression Pattern of Bcl2 in Breast Cancer: A Comprehensive in silico and in vitro Study

B-cell lymphoma/leukemia gene-2 (Bcl-2) is the primary proto-oncogene that has been shown to work by preventing apoptosis/programmed cell death. Bcl-2 combines a variety of cell-generated signals associated to the survival and death of cells. In glioma, lung, and breast cancer, Bcl-2 over-expression has been linked to an increase in invasion and migration. Many treatment regimens that target Bcl2 have been established and approved, and thus increasing the survival rates of the patients. The primary goal of this research was to recognize new therapeutic compounds that target Bcl2 and assess Bcl2 expression pattern in BC patients. We used various bioinformatic tools as well as several in vitro assays to look out the expression and inhibition of Bcl2 in BC. Our study depicted that Bcl2 had a strong connection with tumour stroma, notably with suppressor cells originating from myeloid tissues. Moreover, in vitro and in silico research identified Paclitaxel as a promising natural substance that targets Bcl2. Overall, this work shows that Bcl2 overexpression accelerates the development of BC, and that targeting Bcl2 in combination with other drugs will dramatically improve BC patient's response to treatment and prevent the emergence of drug resistance.

Immunolocalization patterns of histone-deacetylases in salivary glands of mice during postnatal development

OBJECTIVE

Histone-deacetylases (HDACs) are epigenetic modulators involved in the control of gene expression. No data are available on the expression or subcellular localization of HDACs in salivary glands. The present study aims to examine the subcellular distribution of HDACs in salivary glands during postnatal development.

DESIGN

The major salivary glands of C57/BL6 mice were separately removed at 10, 25, 30,60 and 90 days after birth. Hematoxylin-eosin (H&E) and immunohistochemical staining were performed for HDACs. Gene Expression of HDACs in C57BL/6. NOD-Aec1Aec2 mice salivary glands during the development of Sjögren's syndrome-like illness were also analyzed by using the gene expression datasets (GSE 15640).

RESULTS

In the mice salivary gland, HDACs were found to have different localization patterns at various stages of development (10, 25, 30, 60, and 90 days). Apart from HDAC6, ductal cells of salivary glands were the primary sites for HDAC localization. HDAC2, 8, 5, 10 and 11 were expressed at high levels in the salivary gland after birth while HDAC6 showed no expression during postnatal development. This suggests that these HDAC subtypes may have different roles in salivary gland function. In the context of Sjögren's syndrome-like illness, HDAC 2, 8 and 10 showed low expression while HDAC1, 6,5,3 and 11 had relatively high expression in the salivary gland.

CONCLUSIONS

This study has provided an important reference for understanding the spatiotemporal-specific expression of HDACs in the salivary gland. These results offer new clues for the experimenters and hold promise for developing innovative therapeutic strategies for salivary gland-related diseases.

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.

GRAS gene family in rye (Secale cereale L.): genome-wide identification, phylogeny, evolutionary expansion and expression analyses

BACKGROUND

The GRAS transcription factor family plays a crucial role in various biological processes in different plants, such as tissue development, fruit maturation, and environmental stress. However, the GRAS family in rye has not been systematically analyzed yet.

RESULTS

In this study, 67 GRAS genes in S. cereale were identified and named based on the chromosomal location. The gene structures, conserved motifs, cis-acting elements, gene replications, and expression patterns were further analyzed. These 67 ScGRAS members are divided into 13 subfamilies. All members include the LHR I, VHIID, LHR II, PFYRE, and SAW domains, and some nonpolar hydrophobic amino acid residues may undergo cross-substitution in the VHIID region. Interested, tandem duplications may have a more important contribution, which distinguishes them from other monocotyledonous plants. To further investigate the evolutionary relationship of the GRAS family, we constructed six comparative genomic maps of homologous genes between rye and different representative monocotyledonous and dicotyledonous plants. The response characteristics of 19 ScGRAS members from different subfamilies to different tissues, grains at filling stages, and different abiotic stresses of rye were systematically analyzed. Paclobutrazol, a triazole-based plant growth regulator, controls plant tissue and grain development by inhibiting gibberellic acid (GA) biosynthesis through the regulation of DELLA proteins. Exogenous spraying of paclobutrazol significantly reduced the plant height but was beneficial for increasing the weight of 1000 grains of rye. Treatment with paclobutrazol, significantly reduced gibberellin levels in grain in the filling period, caused significant alteration in the expression of the DELLA subfamily gene members. Furthermore, our findings with respect to genes, ScGRAS46 and ScGRAS60, suggest that these two family members could be further used for functional characterization studies in basic research and in breeding programmes for crop improvement.

CONCLUSIONS

We identified 67 ScGRAS genes in rye and further analysed the evolution and expression patterns of the encoded proteins. This study will be helpful for further analysing the functional characteristics of ScGRAS genes.

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.

Genome-wide identification and expression analysis of the glycosyl hydrolase family 1 genes in Medicago sativa revealed their potential roles in response to multiple abiotic stresses

Glycoside hydrolase family 1 (GH1) β-glucosidases (BGLUs), are encoded by a large number of genes, which participate in the development and stress response of plants, particularly under biotic and abiotic stresses through the activation of phytohormones. However, there are few studies systematically analyzing stress or hormone-responsive BGLU genes in alfalfa. In this study, a total of 179 BGLU genes of the glycoside hydrolase family 1 were identified in the genome of alfalfa, and then were classified into five distinct clusters. Sequence alignments revealed several conserved and unique motifs among these MsBGLU proteins. Many cis-acting elements related to abiotic stresses and phytohormones were identified in the promoter of some MsBGLUs. Moreover, RNA-seq and RT-qPCR analyses showed that these MsBGLU genes exhibited distinct expression patterns in response to different abiotic stress and hormonal treatments. In summary, this study suggests that MsBGLU genes play crucial roles in response to various abiotic stresses and hormonal responses, and provides candidate genes for stress tolerance breeding in alfalfa.

Genome-wide identification of the Pinus tabuliformis CONSTANS-like gene family and their potential roles in reproductive cone development

The CONSTANS-like (COL) genes, as a core transcription factor in the photoperiod regulation pathway, play a key role in plant reproduction development. However, their molecular characterization has rarely been studied in Pinus tabuliformis. Here, 10 PtCOL genes were identified in the P. tabuliformis genome and multiple sequence alignments have indicated that the PtCOL proteins contained highly conserved B-BOX1 and CCT domains. Sequence similarity analysis showed that PtCOL1 and PtCOL3 had the higher similarity with Norway spruce COLs (PaCOL2 and PaCOL1) and Arabidopsis COLs (AtCOL3, 4 and 5), respectively. Phylogeny and gene structure analyses revealed that PtCOLs were divided into three subgroups, each with identical or similar distributions of exons, introns, and motifs. Moreover, 10 PtCOLs were distributed on 6 chromosomes and PtCOL9 has syntenic gene pairs in both Ginkgo biloba and Sequoiadendron giganteum. Interestingly, in transcriptome profiles, most PtCOLs exhibited a diurnal oscillation pattern under both long (LD) and short (SD) day conditions. Additionally, PtCOLs were highly expressed in needles and female cones, and showed different spatial expression patterns. Among the ten PtCOLs, PtCOL1/3 heterologous overexpression Arabidopsis displayed a delayed-flowering phenotype under SD, indicating that they are likely to play a crucial role in the reproductive development. Additionally, PtCOL1 and PtCOL3 were not only capable of interacting with each other, but they were each capable of interacting with themselves. Furthermore, PtCOL1 and PtCOL3 were also involved in the MADS-box protein-protein interaction (PPI) network in P. tabuliformis cone development. Direct interactions of PtDAL11 with PtCOL1/3 impeded PtCOL1/3 translocation into the nucleus. In summary, this study provided comprehensive understanding for the functions of the PtCOL gene family and revealed their biological roles in the photoperiod-dependent P. tabuliformis cone development.

Genome-wide characterization of the xyloglucan endotransglucosylase/hydrolase family genes and their response to plant hormone in sugar beet

Xyloglucan endotransglucosylase/hydrolases (XTHs) play a crucial role in plant growth and development. However, their functional response to phytohormone in sugar beet still remains obscure. In this study, we identified 30 putative BvXTH genes in the sugar beet genome. Phylogenetic and evolutionary relationship analysis revealed that they were clustered into three groups and have gone through eight tandem duplication events under purifying selection. Gene structure and motif composition analysis demonstrated that they were highly conserved and all contained one conserved glycoside hydrolase family 16 domain (Glyco_hydro_16) and one xyloglucan endotransglycosylase C-terminus (XET_C) domain. Transcriptional expression analysis exhibited that all BvXTHs were ubiquitously expressed in leaves, root hairs and tuberous roots, and most of them were up-regulated by brassinolide (BR), jasmonic acid (JA), abscisic acid (ABA) and gibberellic acid (GA3). Further mutant complementary experiment demonstrated that expression of BvXTH17 rescued the retarded growth phenotype of xth22, an Arabidopsis knock out mutant of AtXTH22. The findings in our work provide fundamental information on the structure and evolutionary relationship of the XTH family genes in sugar beet, and reveal the potential function of BvXTH17 in plant growth and hormone response.

The potential roles of acid invertase family in Dendrobium huoshanense: Identification, evolution, and expression analyses under abiotic stress

Dendrobium huoshanense, a traditional Chinese medicine prized for its horticultural and medicinal properties, thrives in an unfavorable climate and is exposed to several adverse environmental conditions. Acid invertase (AINV), a widely distributed enzyme that has been demonstrated to play a significant role in response to environmental stresses. However, the identification of the AINV gene family in D. huoshanense, the collinearity between relative species, and the expression pattern under external stress have yet to be resolved. We systematically retrieved the D. huoshanense genome and screened out four DhAINV genes, which were further classified into two subfamilies by the phylogenetic analysis. The evolutionary history of AINV genes in D. huoshanense was uncovered by comparative genomics investigations. The subcellular localization predicted that the DhVINV genes may be located in the vacuole, while the DhCWINV genes may be located in the cell wall. The exon/intron structures and conserved motifs of DhAINV genes were found to be highly conserved in two subclades. The conserved amino acids and catalytic motifs in DhAINV proteins were determined to be critical to their function. Notably, the cis-acting elements in all DhAINV genes were mainly relevant to abiotic stresses and light response. In addition, the expression profile coupled with qRT-PCR revealed the typical expression patterns of DhAINV in response to diverse abiotic stresses. Our findings could be beneficial to the characterization and further investigation of AINV functions in Dendrobium plants.

Genome- and transcriptome-wide identification of trehalose-6-phosphate phosphatases (TPP) gene family and their expression patterns under abiotic stress and exogenous trehalose in soybean

BACKGROUND

Trehalose-6-phosphate phosphatase (TPP) is an essential enzyme catalyzing trehalose synthesis, an important regulatory factor for plant development and stress response in higher plants. However, the TPP gene family in soybean has not been reported.

RESULTS

A comprehensive analysis of the TPP gene family identified 18 GmTPPs classified into eight groups based on the phylogenetic relationships and the conservation of protein in six monocot and eudicot plants. The closely linked subfamilies had similar motifs and intron/exon numbers. Segmental duplication was the main driving force of soybean GmTPPs expansion. In addition, analysis of the cis-regulatory elements and promoter regions of GmTPPs revealed that GmTPPs regulated the response to several abiotic stresses. Moreover, RNA-seq and qRT-PCR analysis of the tissue-specific GmTPPs under different abiotic stresses revealed that most GmTPPs were associated with response to different stresses, including cold, drought, saline-alkali, and exogenous trehalose. Notably, exogenous trehalose treatment up-regulated the expression of most TPP genes under saline-alkali conditions while increasing the carbohydrate and trehalose levels and reducing reactive oxygen species (ROS) accumulation in soybean sprouts, especially in the saline-alkali tolerant genotype. Furthermore, the interaction network and miRNA target prediction revealed that GmTPPs interacted with abiotic stress response-related transcription factors.

CONCLUSIONS

The findings in this study lay a foundation for further functional studies on TPP-based breeding to improve soybean development and stress tolerance.

Transcriptome analysis and expression profiles of odorant binding proteins and chemosensory proteins in Orius sauteri

The flower bug Orius sauteri (Heteroptera: Anthocoridae), is a polyphagous predator and a natural enemy widely used in biological pest control to micro-pests including aphids, spider mites, thrips and so on. In the present study, the transcriptome analysis of adult heads in O. sauteri were performed and identified a total of 38 chemosensory genes including 24 odorant binding proteins (OBPs) and 14 chemosensory proteins (CSPs). Subsequently, we conducted quantitative real-time PCR to detect the tissue expression level of 18 OBPs and 8 CSPs. The results showed that almost all OsauOBPs and OsauCSPs have a high expression level in the adult heads of both sexes. In addition, 5 OsauOBPs (OBP1, OBP2, OBP3, OBP4 and OBP14) have a significantly higher expressed in male heads than female, indicating that these chemosensory proteins might be involved in the male-specific behaviors such as pheromone reception and mate-seeking. This study will provide helpful reference for subsequent understanding of chemoreception mechanism in O. sauteri.

Genome-wide identification of GH9 gene family and the assessment of its role during fruit abscission zone formation in Vaccinium ashei

KEY MESSAGE

The genomic location and stage-specific expression pattern of GH9 genes reveal their critical roles during fruit abscission zone formation in Vaccinium ashei. Glycosyl hydrolase family 9 (GH9) cellulases play a crucial role in both cellulose synthesis and hydrolysis during plant growth and development. Despite this importance, there is currently no study on the involvement of GH9-encoding genes, specifically VaGH9s, in abscission zone formation of rabbiteye blueberries (Vaccinium ashei). In this study, we identified a total of 61 VaGH9s in the genome, which can be classified into 3 subclasses based on conserved motifs and domains, gene structures, and phylogenetic analyses. Our synteny analysis revealed that VaGH9s are more closely related to the GH9s of Populus L. than to those of Arabidopsis, Vitis vinifera, and Citrus sinensis. In silico structural analysis predicted that most of VaGH9s are hydrophilic, and localized in cell membrane and/or cell wall, and the variable sets of cis-acting regulatory elements and functional diversity with four categories of stress response, hormone regulation, growth and development, and transcription factor-related elements are present in the promoter sequence of VaGH9s genes. Transcriptomic analysis showed that there were 22 differentially expressed VaGH9s in fruit abscission zone tissue at the veraison stage, and the expression of VaGH9B2 and VaGH9C10 was continuously increased during fruit maturation, which were in parallel with the increasing levels of cellulase activity and oxidative stress indicators, suggesting that they are involved in the separation stage of fruit abscission in Vaccinium ashei. Our work identified 22 VaGH9s potentially involved in different stages of fruit abscission and would aid further investigation into the molecular regulation of abscission in rabbiteye blueberries fruit.

Lack of the CCT domain changes the ability of mango MiCOL14A to resist salt and drought stress in Arabidopsis

CONSTANS (CO) is the key gene in the photoperiodic pathway that regulates flowering in plants. In this paper, a CONSTANS-like 14A (COL14A) gene was obtained from mango, and its expression patterns and functions were characterized. Sequence analysis shows that MiCOL14A-JH has an additional A base, which leads to code shifting in subsequent coding boxes and loss of the CCT domain. The MiCOL14A-JH and MiCOL14A-GQ genes both belonged to group Ⅲ of the CO/COL gene family. Analysis of tissue expression patterns showed that MiCOL14A was expressed in all tissues, with the highest expression in the leaves of seedling, followed by lower expression levels in the flowers and stems of adult leaves. However, there was no significant difference between different mango varieties. At different development stages of flowering, the expression level of MiCOL14A-GQ was the highest in the leaves before floral induction period, and the lowest at flowering stage, while the highest expression level of MiCOL14A-JH appeared in the leaves at flowering stage. The transgenic functional analysis showed that both MiCOL14A-GQ and MiCOL14A-JH induced delayed flowering of transgenic Arabidopsis. In addition, MiCOL14A-JH enhanced the resistance of transgenic Arabidopsis to drought stress, while MiCOL14A-GQ increased the sensitivity of transgenic Arabidopsis to salt stress. Further proteinprotein interaction analysis showed that MiCOL14A-JH directly interacted with MYB30-INTERACTING E3 LIGASE 1 (MiMIEL1), CBL-interacting protein kinase 9 (MiCIPK9) and zinc-finger protein 4 (MiZFP4), but MiCOL14A-GQ could not interact with these three stress-related proteins. Together, our results demonstrated that MiCOL14A-JH and MiCOL14A-GQ not only regulate flowering but also play a role in the abiotic stress response in mango, and the lack of the CCT domain affects the proteinprotein interaction, thus affecting the gene response to stress. The insertion of an A base can provide a possible detection site for mango resistance breeding.

Genome-wide analysis of the Pleurotus eryngii laccase gene (PeLac) family and functional identification of PeLac5

The laccase gene family encodes multiple isozymes that are crucial for the degradation of substrates and the regulation of developmental processes in fungi. Pleurotus eryngii is an important edible and medicinal fungus belonging to the Basidiomycota phylum and can grow on a variety of natural substrates. In the present study, genome-wide profiling of P. eryngii identified 10 genes encoding its laccase isoenzymes. Conservative sequence analysis demonstrated that all PeLacs possess classical laccase structural domains. Phylogenetic analysis yielded four major subgroups, the members of which are similar with respect to conserved gene organization, protein domain architecture, and consensus motifs. The 10 PeLacs formed three groups together with 12 PoLacs in Pleurotus ostreatus, indicating that they share a high level of evolutionary homology. Cis-responsive element analysis implied that PeLacs genes play a role in growth and development and lignocellulose degradation. Targeted overexpression of PeLac5 reduced the time to primordia formation and their development to fruiting bodies. Gene expression patterns in the presence of different lignocellulosic substrates indicate that three PeLacs genes (2, 4, and 9) are key to lignocellulose degradation. This work presents the first inventory of laccase genes in P. eryngii and preliminarily explores their functions, which may help to uncover the manner by which these proteins utilize substrates.

Genome-wide analysis of the heat shock transcription factor family reveals saline-alkali stress responses in Xanthoceras sorbifolium

The heat shock transcription factor (HSF) family is involved in regulating growth, development, and abiotic stress. The characteristics and biological functions of HSF family member in X. sorbifolium, an important oil and ornamental plant, have never been reported. In this study, 21 XsHSF genes were identified from the genome of X. sorbifolium and named XsHSF1-XsHSF21 based on their chromosomal positions. Those genes were divided into three groups, A, B, and C, containing 12, one, and eight genes, respectively. Among them, 20 XsHSF genes are located on 11 chromosomes. Protein structure analysis suggested that XsHSF proteins were conserved, displaying typical DNA binding domains (DBD) and oligomerization domains (OD). Moreover, HSF proteins within the same group contain specific motifs, such as motif 5 in the HSFC group. All XsHSF genes have one intron in the CDS region, except XsHSF1 which has two introns. Promoter analysis revealed that in addition to defense and stress responsiveness elements, some promoters also contained a MYB binding site and elements involved in multiple hormones responsiveness and anaerobic induction. Duplication analysis revealed that XsHSF1 and XsHSF4 genes were segmentally duplicated while XsHSF2, XsHSF9, and XsHSF13 genes might have arisen from transposition. Expression pattern analysis of leaves and roots following salt-alkali treatment using qRT-PCR indicated that five XsHSF genes were upregulated and one XsHSF gene was downregulated in leaves upon NaCl treatment suggesting these genes may play important roles in salt response. Additionally, the expression levels of most XsHSFs were decreased in leaves and roots following alkali-induced stress, indicating that those XsHSFs may function as negative regulators in alkali tolerance. MicroRNA target site prediction indicated that 16 of the XsHSF genes may be regulated by multiple microRNAs, for example XsHSF2 might be regulated by miR156, miR394, miR395, miR408, miR7129, and miR854. And miR164 may effect the mRNA levels of XsHSF3 and XsHSF17, XsHSF9 gene may be regulated by miR172. The expression trends of miR172 and miR164 in leaves and roots on salt treatments were opposite to the expression trend of XsHSF9 and XsHSF3 genes, respectively. Promoter analysis showed that XsHSFs might be involved in light and hormone responses, plant development, as well as abiotic stress responses. Our results thus provide an overview of the HSF family in X. sorbifolium and lay a foundation for future functional studies to reveal its roles in saline-alkali response.

Identification of Dioscorea opposite Thunb. CDPK gene family reveals that DoCDPK20 is related to heat resistance

Temperature affects the growth and yield of yam (Dioscorea opposite Thunb.), and calcium-dependent protein kinases (CDPKs) play an important role in the plant stress response. However, there has been a lack of system analyses of yam's CDPK gene family. In this study, 29 CDPK transcriptome sequences with complete open reading frames (ORFs) were identified from yam RNA sequencing data. The sequences were classified into four groups (I-VI) using phylogenetic analysis. Two DoCDPK genes were randomly selected from each group and the gene patterns of yam leaves were determined using quantitative real-time PCR (qRT-PCR) under high and low temperature stress in order to show their unique functions in mediating specific responses. Among them, DoCDPK20 was significantly induced in high temperatures. The pPZP221-DoCDPK20 was transformed into tobacco leaves using an agrobacterium-mediated method. Under high temperature stress, DoCDPK20 overexpression reduced photosynthesis and improved heat tolerance in transgenic tobacco. Our research offers meaningful perspectives into CDPK genes and new avenues for the genetic engineering and molecular breeding of yam.

Characterization and in silico analysis of the domain unknown function DUF568-containing gene family in rice (Oryza sativa L.)

BACKGROUND

Domains of unknown function (DUF) proteins are a number of uncharacterized and highly conserved protein families in eukaryotes. In plants, some DUFs have been predicted to play important roles in development and response to abiotic stress. Among them, DUF568-containing protein family is plant-specific and has not been described previously. A basic analysis and expression profiling was performed, and the co-expression and interaction networks were constructed to explore the functions of DUF568 family in rice.

RESULTS

The phylogenetic tree showed that the 8, 9 and 11 DUF568 family members from rice, Arabidopsis and maize were divided into three groups. The evolutionary relationship between DUF568 members in rice and maize was close, while the genes in Arabidopsis were more distantly related. The cis-elements prediction showed that over 82% of the elements upstream of OsDUF568 genes were responsive to light and phytohormones. Gene expression profile prediction and RT-qPCR experiments revealed that OsDUF568 genes were highly expressed in leaves, stems and roots of rice seedling. The expression of some OsDUF568 genes varied in response to plant hormones (abscisic acid, 6-benzylaminopurine) and abiotic stress (drought and chilling). Further analysis of the co-expression and protein-protein interaction networks using gene ontology showed that OsDUF568 - related genes were enriched in cellular transports, metabolism and processes.

CONCLUSIONS

In summary, our findings suggest that the OsDUF568 family may be a vital gene family for the development of rice roots, leaves and stems. In addition, the OsDUF568 family may participate in abscisic acid and cytokinin signaling pathways, and may be related to abiotic stress resistance in these vegetative tissues of rice.

Genome-wide identification of the NAC gene family and its functional analysis in Liriodendron

As one of the largest plant specific transcription factor families, NAC family members play an important role in plant growth, development and stress resistance. To investigate the function of NAC transcription factors during abiotic stress, as well as during somatic embryogenesis, we identified and characterized the NAC gene family in Liriodendron chinense. We found that most LcNAC members contain more than three exons, with a relatively conserved gene and motif structure, especially at the N-terminus. Interspecies collinearity analysis revealed a closer relationship between the L. chinense NACs and the P. trichocarpa NACs. We analyzed the expression of LcNAC in different tissues and under three abiotic stresses. We found that 12 genes were highly expressed during the ES3 and ES4 stages of somatic embryos, suggesting that they are involved in the development of somatic embryos. 6 LcNAC genes are highly expressed in flower organs. The expression pattern analysis of LcNACs based on transcriptome data and RT-qPCR obtained from L. chinense leaves indicated differential expression responses to drought, cold, and heat stress. Genes in the NAM subfamily expressed differently during abiotic stress, and LcNAC6/18/41/65 might be the key genes in response to abiotic stress. LcNAC6/18/41/65 were cloned and transiently transformed into Liriodendron protoplasts, where LcNAC18/65 was localized in cytoplasm and nucleus, and LcNAC6/41 was localized only in nucleus. Overall, our findings suggest a role of the NAC gene family during environmental stresses in L. chinense. This research provides a basis for further study of NAC genes in Liriodendron chinense.

Discovery of unique CYP716C oxidase involved in pentacyclic triterpene biosynthesis from Camptotheca acuminata

Dozens of triterpenes have been isolated from Camptotheca acuminata, however, triterpene metabolism in this plant remains poorly understood. The common C28 carboxy located in the oleanane-type and ursane-type triterpenes indicates the existence of a functionally active triterpene, C28 oxidase, in this plant. Thorough mining and screening of the CYP716 genes were initiated using the multi-omics database for C. acuminata. Two CYP716A (CYP716A394 and CYP716A395) and three CYP716C (CYP716C80-CYP716C82) were identified based on conserved domain analyses and hierarchical cluster analyses. CYP716 microsomal proteins were prepared and their enzymatic activities were evaluated in vitro. The CYP716 classified into the CYP716C subfamily displays β-amyrin oxidation activity, and CYP716A displays α-amyrin and lupeol oxidation activity, based on gas chromatography-mass spectrometry analyses. The oxidation products were determined based on their mass and nuclear magnetic resonance spectrums. The optimum reaction conditions and kinetic parameters for CYP716C were determined, and functions were verified in Nicotiana benthaminana. Relative quantitative analyses revealed that these CYP716C genes were enriched in the leaves of C. acuminata plantlets after 60 d. These results indicate that CYP716C plays a dominant role in oleanane-type triterpene metabolism in the leaves of C. acuminata via a substrate-specific manner, and CYP716A is responsible for ursane- and lupane-type triterpene metabolism in fruit. This study provides valuable insights into the unique CYP716C-mediated oxidation step of pentacyclic triterpene biosynthesis in C. acuminata.

The TRAF gene family in turbot (Scophthalmus maximus): Identification, characterization, molecular evolution and expression patterns analysis

Tumor necrosis factor receptor-associated factor (TRAF) is an important structural protein, which can bind to TNF receptors and participate in the regulation of TNF signaling pathway. Nonetheless, few studies have been conducted to investigate the systematic identification of TRAF gene family in teleost and role in innate immunity of turbot (Scophthalmus maximus). In this study, eight TRAF genes, namely SmTRAF2aa, SmTRAF2ab, SmTRAF2b, SmTRAF3, SmTRAF4a, SmTRAF5, SmTRAF6 and SmTRAF7, were identified and annotated in turbot by using bioinformatics methods. Analysis of the phylogenetic, syntenic and molecular evolution demonstrated that all SmTRAF members were evolutionarily conserved in teleost. Domain analysis showed all SmTRAF proteins contained a typical conserved N-terminal RING finger domain. Most SmTRAF proteins contained a MATH domain at the C-terminal, while SmTRAF7 contains seven duplicate WD40 domains. In addition, quantitative real-time PCR was performed to detect the expression patterns of SmTRAFs in tissues from healthy and Vibrio anguillarum infected turbots. The results indicated SmTRAFs had diverse tissue expression patterns and the expression of TRAF gene changed significantly after V. anguillarum infection. This study provided a basis for understanding the roles of TRAFs in the innate immune response of turbot.

Genome-wide analysis of the MADS-box gene family in Lonicera japonica and a proposed floral organ identity model

BACKGROUND

Lonicera japonica Thunb. is widely used in traditional Chinese medicine. Medicinal L. japonica mainly consists of dried flower buds and partially opened flowers, thus flowers are an important quality indicator. MADS-box genes encode transcription factors that regulate flower development. However, little is known about these genes in L. japonica.

RESULTS

In this study, 48 MADS-box genes were identified in L. japonica, including 20 Type-I genes (8 Mα, 2 Mβ, and 10 Mγ) and 28 Type-II genes (26 MIKCc and 2 MIKC*). The Type-I and Type-II genes differed significantly in gene structure, conserved domains, protein structure, chromosomal distribution, phylogenesis, and expression pattern. Type-I genes had a simpler gene structure, lacked the K domain, had low protein structure conservation, were tandemly distributed on the chromosomes, had more frequent lineage-specific duplications, and were expressed at low levels. In contrast, Type-II genes had a more complex gene structure; contained conserved M, I, K, and C domains; had highly conserved protein structure; and were expressed at high levels throughout the flowering period. Eleven floral homeotic MADS-box genes that are orthologous to the proposed Arabidopsis ABCDE model of floral organ identity determination, were identified in L. japonica. By integrating expression pattern and protein interaction data for these genes, we developed a possible model for floral organ identity determination.

CONCLUSION

This study genome-widely identified and characterized the MADS-box gene family in L. japonica. Eleven floral homeotic MADS-box genes were identified and a possible model for floral organ identity determination was also developed. This study contributes to our understanding of the MADS-box gene family and its possible involvement in floral organ development in L. japonica.

Genome-wide identification and expression analyses of the pectate lyase (PL) gene family in Fragaria vesca

BACKGROUND

Pectate lyase (PL, EC 4.2.2.2), as an endo-acting depolymerizing enzyme, cleaves α-1,4-glycosidic linkages in esterified pectin and involves a broad range of cell wall modifications. However, the knowledge concerning the genome-wide analysis of the PL gene family in Fragaria vesca has not been thoroughly elucidated.

RESULTS

In this study, sixteen PLs members in F. vesca were identified based on a genome-wide investigation. Substantial divergences existed among FvePLs in gene duplication, cis-acting elements, and tissue expression patterns. Four clusters were classified according to phylogenetic analysis. FvePL6, 8 and 13 in cluster II significantly contributed to the significant expansions during evolution by comparing orthologous PL genes from Malus domestica, Solanum lycopersicum, Arabidopsis thaliana, and Fragaria×ananassa. The cis-acting elements implicated in the abscisic acid signaling pathway were abundant in the regions of FvePLs promoters. The RNA-seq data and in situ hybridization revealed that FvePL1, 4, and 7 exhibited maximum expression in fruits at twenty days after pollination, whereas FvePL8 and FvePL13 were preferentially and prominently expressed in mature anthers and pollens. Additionally, the co-expression networks displayed that FvePLs had tight correlations with transcription factors and genes implicated in plant development, abiotic/biotic stresses, ions/Ca2+, and hormones, suggesting the potential roles of FvePLs during strawberry development. Besides, histological observations suggested that FvePL1, 4 and 7 enhanced cell division and expansion of the cortex, thus negatively influencing fruit firmness. Finally, FvePL1-RNAi reduced leaf size, altered petal architectures, disrupted normal pollen development, and rendered partial male sterility.

CONCLUSION

These results provide valuable information for characterizing the evolution, expansion, expression patterns and functional analysis, which help to understand the molecular mechanisms of the FvePLs in the development of strawberries.

Magnesium transporter CsMGT10 of tea plants plays a key role in chlorosis leaf vein greening

Magnesium (Mg2+), as the central atom of chlorophyll, is the most abundant divalent cation for plant growth and development in living cells. MRS2/MGT magnesium transporters play important roles in coping with magnesium stress, chloroplast development and photosynthesis. However, the molecular mechanism of MGT influencing tea plant leaf vein color remains unknown. Here, we demonstrate that CsMGT10 may be a potential transporter influencing leaf vein color. CsMGT10 belongs to Clade A member of MRS2/MGT family. CsMGT10 has the highest expression level in leaves of tea plants. And it is mainly expressed in aboveground parts, especially in vascular bundles. Moreover, CsMGT10 localizes to the chloroplast envelope of tea plants with a high affinity to Mg2+. And the GMN motif is required for its magnesium transport function. Ectopic expression of CsMGT10 in Arabidopsis leaf variegation mutant var5-1 can restore green color of chlorosis leaf veins, and the contents of chlorophyll and carotenoid change significantly, proving its essential role in leaf vein greening. Furthermore, the chlorophyll and carotenoid of tea leaves treated with CsMGT10 antisense oligonucleotides also decrease significantly. Our findings indicate that CsMGT10 mainly acts as Mg2+ transporter in chloroplast envelope of leaf veins, which may play a key role in leaf vein greening of tea plants.

Image-based quantification of histological features as a function of spatial location using the Tissue Positioning System

BACKGROUND

Tissues such as the liver lobule, kidney nephron, and intestinal gland exhibit intricate patterns of zonated gene expression corresponding to distinct cell types and functions. To quantitatively understand zonation, it is important to measure cellular or genetic features as a function of position along a zonal axis. While it is possible to manually count, characterize, and locate features in relation to the zonal axis, it is labor-intensive and difficult to do manually while maintaining precision and accuracy.

METHODS

We addressed this challenge by developing a deep-learning-based quantification method called the "Tissue Positioning System" (TPS), which can automatically analyze zonation in the liver lobule as a model system.

FINDINGS

By using algorithms that identified vessels, classified vessels, and segmented zones based on the relative position along the portal vein to central vein axis, TPS was able to spatially quantify gene expression in mice with zone specific reporters.

INTERPRETATION

TPS could discern expression differences between zonal reporter strains, ages, and disease states. TPS could also reveal the zonal distribution of cells previously thought to be positioned randomly. The design principles of TPS could be generalized to other tissues to explore the biology of zonation.

FUNDING

CPRIT (RP190208, RP220614, RP230330) and NIH (P30CA142543, R01AA028791, R01CA251928, R01DK1253961, R01GM140012, 1R01GM141519, 1R01DE030656, 1U01CA249245). The Pollack Foundation, Simmons Comprehensive Cancer Center Cancer & Obesity Translational Pilot Award, and the Emerging Leader Award from the Mark Foundation For Cancer Research (#21-003-ELA).

Genome-wide identification, characterization and transcriptional profile of the SWEET gene family in Dendrobium officinale

BACKGROUND

Dendrobium officinale Kimura et Migo (D. officinale) is a well-known traditional Chinese medicine with high content polysaccharides in stems. The SWEET (Sugars Will Eventually be Exported Transporters) family is a novel class of sugar transporters mediating sugar translocation among adjacent cells of plants. The expression patterns of SWEETs and whether they are associated with stress response in D. officinale remains uncovered.

RESULTS

Here, 25 SWEET genes were screened out from D. officinale genome, most of which typically contained seven transmembrane domains (TMs) and harbored two conserved MtN3/saliva domains. Using multi-omics data and bioinformatic approaches, the evolutionary relationship, conserved motifs, chromosomal location, expression patterns, correlationship and interaction network were further analyzed. DoSWEETs were intensively located in nine chromosomes. Phylogenetic analysis revealed that DoSWEETs were divided into four clades, and conserved motif 3 specifically existed in DoSWEETs from clade II. Different tissue-specific expression patterns of DoSWEETs suggested the division of their roles in sugar transport. In particular, DoSWEET5b, 5c, and 7d displayed relatively high expression levels in stems. DoSWEET2b and 16 were significantly regulated under cold, drought, and MeJA treatment, which were further verified using RT-qPCR. Correlation analysis and interaction network prediction discovered the internal relationship of DoSWEET family.

CONCLUSIONS

Taken together, the identification and analysis of the 25 DoSWEETs in this study provide basic information for further functional verification in D. officinale.

Genome-wide identification and immune response analysis of mitogen-activated protein kinase cascades in tea geometrid, Ectropis grisescens Warren (Geometridae, Lepidoptera)

BACKGROUND

Tea geometrid Ectropis grisescens (Geometridae: Lepidoptera), is one of the most destructive defoliators in tea plantations in China. The MAPK cascade is known to be an evolutionarily conserved signaling module, acting as pivotal cores of host-pathogen interactions. Although the chromosome-level reference genome of E. grisescens was published, the whole MAPK cascade gene family has not been fully identified yet, especially the expression patterns of MAPK cascade gene family members upon an ecological biopesticide, Metarhizium anisopliae, remains to be understood.

RESULTS

In this study, we have identified 19 MAPK cascade gene family members in E. grisescens, including 5 MAPKs, 4 MAP2Ks, 8 MAP3Ks, and 2 MAP4Ks. The molecular evolution characteristics of the whole Eg-MAPK cascade gene family, including gene structures, protein structural organization, chromosomal localization, orthologs construction and gene duplication, were systematically investigated. Our results showed that the members of Eg-MAPK cascade gene family were unevenly distributed in 13 chromosomes, and the clustered members in each group shared similar structures of the genes and proteins. Gene expression data revealed that MAPK cascade genes were expressed in all four developmental stages of E. grisescens and were fairly and evenly distributed in four different larva tissues. Importantly, most of the MAPK cascade genes were induced or constitutively expressed upon M. anisopliae infection.

CONCLUSIONS

In summary, the present study was one of few studies on MAPK cascade gene in E. grisescens. The characterization and expression profiles of Eg-MAPK cascades genes might help develop new ecofriendly biological insecticides to protect tea trees.

Genome-wide identification and characterization of Glutathione S-Transferases (GSTs) and their expression profile under abiotic stresses in tobacco (Nicotiana tabacum L.)

BACKGROUND

Glutathione S-transferases (GSTs) are large and multifunctional proteases that play an important role in detoxification, protection against biotic and abiotic stresses, and secondary metabolite transportation which is essential for plant growth and development. However, there is limited research on the identification and function of NtGSTs.

RESULTS

This study uses K326 and other six tobacco varieties (Hongda, HG, GDH11, Va116, VG, and GDH88) as materials to conduct comprehensive genome-wide identification and functional characterization of the GST gene in tobacco. A total of 59 NtGSTs were identified and classified into seven subfamilies via the whole-genome sequence analysis, with the Tau type serving as the major subfamily. The NtGSTs in the same branch of the evolutionary tree had similar exon/intron structure and motif constitution. There were more than 42 collinear blocks between tobacco and pepper, tomato, and potato, indicating high homology conservation between them. Twelve segmental duplicated gene pairs and one tandem duplication may have had a substantial impact on the evolution and expansion of the tobacco GST gene family. The RT-qPCR results showed that the expression patterns of NtGSTs varied significantly among tissues, varieties, and multiple abiotic stresses, suggesting that NtGST genes may widely respond to various abiotic stresses and hormones in tobacco, including NtGSTF4, NtGSTL1, NtGSTZ1, and NtGSTU40.

CONCLUSIONS

This study provides a comprehensive analysis of the NtGST gene family, including structures and functions. Many NtGSTs play a critical regulatory role in tobacco growth and development, and responses to abiotic stresses. These findings offer novel and valuable insights for understanding the biological function of NtGSTs and the reference materials for cultivating highly resistant varieties and enhancing the yield and quality of crops.

Whole-genome identification and expression profiling of growth-regulating factor (GRF) and GRF-interacting factor (GIF) gene families in Panax ginseng

BACKGROUND

Panax ginseng is a perennial herb and one of the most widely used traditional medicines in China. During its long growth period, it is affected by various environmental factors. Past studies have shown that growth-regulating factors (GRFs) and GRF-interacting factors (GIFs) are involved in regulating plant growth and development, responding to environmental stress, and responding to the induction of exogenous hormones. However, GRF and GIF transcription factors in ginseng have not been reported.

RESULTS

In this study, 20 GRF gene members of ginseng were systematically identified and found to be distributed on 13 chromosomes. The ginseng GIF gene family has only ten members, which are distributed on ten chromosomes. Phylogenetic analysis divided these PgGRFs into six clades and PgGIFs into two clades. In total, 18 of the 20 PgGRFs and eight of the ten PgGIFs are segmental duplications. Most PgGRF and PgGIF gene promoters contain some hormone- and stress- related cis-regulatory elements. Based on the available public RNA-Seq data, the expression patterns of PgGRF and PgGIF genes were analysed from 14 different tissues. The responses of the PgGRF gene to different hormones (6-BA, ABA, GA3, IAA) and abiotic stresses (cold, heat, drought, and salt) were studied. The expression of the PgGRF gene was significantly upregulated under GA3 induction and three weeks of heat treatment. The expression level of the PgGIF gene changed only slightly after one week of heat treatment.

CONCLUSIONS

The results of this study may be helpful for further study of the function of PgGRF and PgGIF genes and lay a foundation for further study of their role in the growth and development of Panax ginseng.

Genome-wide identification and expression of monoacylglycerol lipase (MAGL) gene family in peanut (Arachis hypogaea L.) and functional analysis of AhMGATs in neutral lipid metabolism

Monoacylglycerol lipase (MAGL) involved in regulating plant growth and development and stress responses, hydrolyzes monoacylglycerol (MAG) into free fatty acid and glycerol, which is the last step of triacylglycerol (TAG) breakdown. Here, a genome-wide characterization of MAGL gene family from cultivated peanut (Arachis hypogaea L.) was performed. In total, 24 MAGL genes were identified and unevenly distributed on 14 chromosomes, encoding 229-414 amino acids with molecular weights ranging from 25.91 to 47.01 kDa. Spatiotemporal and stress-induced expression was analyzed by qRT-PCR. Multiple sequence alignment revealed that AhMAGL1a/b and AhMAGL3a/b were the only four bifunctional enzymes with conserved regions of hydrolase and acyltransferase, which could also be named as AhMGATs. GUS histochemical assay showed that AhMAGL1a and -1b were strongly expressed in all tissues of the plants; whereas both AhMAGL3a and -3b were weakly expressed in plants. Subcellular localization analysis indicated that AhMGATs were localized in the endoplasmic reticulum and/or Golgi complex. Seed-specific overexpression of AhMGATs in Arabidopsis decreased the oil content of the seeds and altered the fatty acid compositions, indicating that AhMGATs were involved in TAG breakdown but not TAG biosynthesis in plant seeds. This study lays the foundation for better understanding AhMAGL genes biological function in planta.

Selection and verification of reliable internal reference genes in stem development of herbaceous peony (Paeonia lactiflora Pall.)

Herbaceous peony (Paeonia lactiflora Pall.) has emerged in the cut flower market due to its beautiful appearance. The bending flower stems caused by a lack of mechanical strength is the main problem restricting the development of the cut P. lactiflora industry. So it is of great worth to reveal the basis of stem development changes in P. lactiflora to improve its cut flower quality. Quantitative research on gene expression characteristics can provide clues for understanding their biological functions, and the screening of relatively stable expression genes is a prerequisite for the quantitative study of gene expression characteristics. Thus, it is necessary to find appropriate genes during stem development so as to analyze the qRT‒PCR results. In this study, 10 genes were screened, and these genes expressed stably in stems of different stem strengths at three different developmental stages. Then, their expressions were evaluated by RefFinder, BestKeeper, NormFinder, and GeNorm programs. The results demonstrated that γ-tubulin (γ-TUB) was the most suitable gene, followed by α-tubulin (α-TUB) and β-D-glucosidase (β-GUS), whereas histone H3 (His) was the least suitable gene. Besides, the temporal and spatial expression characteristics of PlCOMT1, the key gene concerned with the synthesis of cell wall fillers in P. lactiflora, were also used to evaluate the suitability of genes. Consequently, γ-TUB and α-TUB are the two best combinations during stem development, and their combination can be used for the stem development of P. lactiflora. These findings will provide a reference for the selection of genes related to stem development and the study of molecular mechanisms related to stem development in P. lactiflora.

Supplementary Information

The online version contains supplementary material available at 10.1007/s12298-023-01325-5.

The heavy metal-associated isoprenylated plant protein (HIPP) gene family plays a crucial role in cadmium resistance and accumulation in the tea plant (Camellia sinensis L.)

Heavy metal-associated isoprenylated plant proteins (HIPPs) are only distributed in vascular plants, and are essential for the detoxification and vascular transport of heavy metals in plants. However, the HIPP gene family has not been thoroughly explored in the tea plant (Camellia sinensis). In this study, we systematically identified 56C. sinensis CsHIPP genes from five groups and characterized their phylogeny, structures, and the features of the encoded proteins. The expression patterns of CsHIPP genes in various tissues of C. sinensis were investigated based on a previous RNA-seq data analysis. The expression patterns of CsHIPP genes were explored in cadmium (Cd)-treated C. sinensis roots using our RNA-seq data. Three CsHIPP genes (CsHIPP22, CsHIPP24, and CsHIPP36) with high expression levels in Cd-treated C. sinensis roots were selected as candidate genes associated with Cd tolerance. Overexpression of CsHIPP22, CsHIPP24, and CsHIPP36 in a yeast mutant (ycf1) rescued Cd-sensitive ycf1 yeast and increased the yeast resistance to Cd stress, implying that these three CsHIPPs might be involved in Cd tolerance. These findings will enable the roles of HIPPs in Cd absorption and detoxification to be better understood as well as improving our understanding of the Cd-resistance and Cd-accumulation mechanisms in tea plant.

Dissecting the features of TGA gene family in Saccharum and the functions of ScTGA1 under biotic stresses

Sugarcane is an important sugar and energy crop and smut disease caused by Sporisorium scitamineum is a major fungal disease which can seriously reduce the yield and quality of sugarcane. In plants, TGACG motif binding (TGA) transcription factors are involved in the regulation of salicylic acid (SA) and methyl jasmonate (MeJA) signaling pathways, as well as in response to various biotic and abiotic stresses. However, no TGA-related transcription factor has been reported in Saccharum. In the present study, 44 SsTGA genes were identified from Saccharum spontaneum, and were assorted into three clades (I, II, III). Cis-regulatory elements (CREs) analysis revealed that SsTGA genes may be involved in hormone and stress response. RNA-seq data and RT-qPCR analysis indicated that SsTGAs were constitutively expressed in different tissues and induced by S. scitamineum stress. In addition, a ScTGA1 gene (GenBank accession number ON416997) was cloned from the sugarcane cultivar ROC22, which was homologous to SsTGA1e in S. spontaneum and encoded a nucleus protein. It was constitutively expressed in sugarcane tissues and up-regulated by SA, MeJA and S. scitamineum stresses. Furthermore, transient overexpression of ScTGA1 in Nicotiana benthamiana could enhance its resistance to the infection of Ralstonia solanacearum and Fusarium solani var. coeruleum, by regulating the expression of immune genes related to hypersensitive response (HR), ethylene (ET), SA and jasmonic acid (JA) pathways. This study should contribute to our understanding on the evolution and function of the SsTGA gene family in Saccharum, and provide a basis for the functional identification of ScTGA1 under biotic stresses.

Genome-wide analysis of UDP-glycosyltransferases family and identification of UGT genes involved in abiotic stress and flavonol biosynthesis in Nicotiana tabacum

BACKGROUND

Uridine disphosphate (UDP) glycosyltransferases (UGTs) act upon a huge variety of highly diverse and complex substrates, such as phytohormones and specialized metabolites, to regulate plant growth, development, disease resistance, and environmental interactions. However, a comprehensive investigation of UGT genes in tobacco has not been conducted.

RESULTS

In this study, we carried out a genome-wide analysis of family-1 UDP glycosyltransferases in Nicotiana tabacum. We predicted 276 NtUGT genes, which were classified into 18 major phylogenetic subgroups. The NtUGT genes were invariably distributed among all the 24 chromosomes with structural diversity in exon/intron structure, conserved motifs, and cis-acting elements of promoters. Three groups of proteins which involved in flavonoid biosynthesis, plant growth and development, transportation and modification were identified that interact with NtUGT proteins using the PPI analysis. Expression analysis of NtUGT genes in cold stress, drought stress and different flower color using both online RNA-Seq data and the realtime PCR analysis, suggested the distinct role of NtUGT genes in resistance of cold, drought and in flavonoid biosynthesis. The enzymatic activities of seven NtUGT proteins that potentially involved in flavonoid glycosylation were analyzed, and found that all seven exhibited activity on myricetin; six (NtUGT108, NtUGT123, NtUGT141, NtUGT155, NtUGT179, and NtUGT195) showed activity on cyanidin; and three (NtUGT108, NtUGT195, and NtUGT217) were active on the flavonol aglycones kaempferol and quercetin, which catalyzing the substrates (myricetin, cyanidin or flavonol) to form new products. We further investigated the enzymatic products and enzymatic properties of NtUGT108, NtUGT195, and NtUGT217, suggested their diverse enzymatic activity toward flavonol, and NtUGT217 showed the highest catalyzed efficient toward quercetin. Overexpression of NtUGT217 significantly increase the content levels of the quercetin-3-O-glucoside, quercetin-3-O-rutinoside and kaempferol-3-O-rutinoside in transgenic tobacco leaves.

CONCLUSION

We identified 276 UGT genes in Nicotiana tabacum. Our study uncovered valuable information about the phylogenetic structure, distribution, genomic characters, expression patterns and enzymatic activity of NtUGT genes in tobacco. We further identified three NtUGT genes involved in flavonoid biosynthesis, and overexpressed NtUGT217 to validate its function in catalyze quercetin. The results provide key candidate NtUGT genes for future breeding of cold and drought resistance and for potential metabolic engineering of flavonoid compounds.

Global expression and functional analysis of human piRNAs during HSV-1 infection

Piwi-interacting RNAs (piRNAs) are a class of non-coding RNAs that play a key role in spermatogenesis. However, little is known about their expression characterization and role in somatic cells infected with herpes simplex virus type 1 (HSV-1). In this study, we systematically investigated the cellular piRNA expression profiles of HSV-1-infected human lung fibroblasts. Compared with the control group, 69 differentially expressed piRNAs were identified in the infection group, among which 52 were up-regulated and 17 were down-regulated. The changes in the expression of 8 piRNAs were further verified by RT-qPCR with a similar expression trend. Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis showed that the target genes of piRNAs were mainly involved in antiviral immunity and various human disease-related signaling pathways. Furthermore, we tested the effects of four up-regulated piRNAs on viral replication by transfecting piRNA mimics. The results showed that the virus titers of the group transfected with piRNA-hsa-28,382 (alias piR-36,233) mimic decreased significantly, and that of the group transfected piRNA-hsa-28,190 (alias piR-36,041) mimic significantly increased. Overall, our results revealed the expression characteristics of piRNAs in HSV-1-infected cells. We also screened two piRNAs that potentially regulate HSV-1 replication. These results may promote a better understanding of the regulatory mechanism of pathophysiological changes induced by HSV-1 infection.

Spatiotemporal expression pattern of the chicken glucocorticoid receptor during early embryonic development

Glucocorticoids - commonly known as stress hormones - belong to the family of steroid hormones and regulate numerous life essential physiological processes. As lipophilic molecules, glucocorticoids are known to cross the placental barrier in mammals, which - applied for therapeutic reasons or arising from environmental influences - illustrates the role of prenatal stress during embryonic developmental processes. The hormones employ their functions by binding to the glucocorticoid receptor (GR) and thus are involved in regulating the transcription of thousands of genes. Therefore, the aim of this study was to investigate the spatiotemporal expression pattern of the GR during early embryonic vertebrate development, using the chicken embryo as a model organism. The results should contribute to enhance and expand the current understanding of glucocorticoid signaling. By performing in-situ hybridization on whole mount chicken embryos from stage HH10 to HH29 and analyzing vibratome sections of hybridized embryos, we described the spatiotemporal expression pattern of the GR during early embryogenesis. Moreover, we compared the expression pattern of the GR with other developmental markers such as Pax7, Desmin, MyoD and HNK-1 using double in-situ hybridization and immunohistochemistry. We were able to determine the first emergence of GR expression in stage HH13 of chicken development in the cranial area, especially in the muscle anlagen of the branchial arches and of non-somitic neck muscles. Furthermore, we monitored the extension of GR expression pattern throughout later stages and found transcripts of GR during somitogenesis, limb development, myogenesis, neurulation and neural differentiation and moreover during organogenesis of the gastrointestinal organs, the heart, the kidneys and the lungs. Toward later stages, GR expression transitioned from more distinct areas of expression to an increasingly ubiquitous expression pattern. Our results support the notion of an enormous relevance of glucocorticoid signaling during vertebrate embryonic development and contribute to a better understanding of the consequences of prenatal stress and the clinical administration of prenatal glucocorticoids.

Genome-wide analysis of the GLP gene family and overexpression of GLP1-5-1 to promote lignin accumulation during early somatic embryo development in Dimocarpus longan

Longan (Dimocarpus longan Lour.) is an economically important subtropical fruit tree. Its fruit quality and yield are affected by embryo development. As a plant seed germination marker gene, the germin-like protein (GLP) gene plays an important role in embryo development. However, the mechanism underlying the role of the GLP gene in somatic embryos is still unclear. Therefore, we conducted genome-wide identification of the longan GLP (DlGLP) gene and preliminarily verified the function of DlGLP1-5-1. Thirty-five genes were identified as longan GLP genes and divided into 8 subfamilies. Based on transcriptome data and qRT‒PCR results, DlGLP genes exhibited the highest expression levels in the root, and the expression of most DlGLPs was upregulated during the early somatic embryogenesis (SE) in longan and responded to high temperature stress and 2,4-D treatment; eight DlGLP genes were upregulated under MeJA treatment, and four of them were downregulated under ABA treatment. Subcellular localization showed that DlGLP5-8-2 and DlGLP1-5-1 were located in the cytoplasm and extracellular stroma/chloroplast, respectively. Overexpression of DIGLP1-5-1 in the globular embryos (GEs) of longan promoted the accumulation of lignin and decreased the H₂O₂ content by regulating the activities of ROS-related enzymes. The results provide a reference for the functional analysis of DlGLPs and related research on improving lignin accumulation in the agricultural industry through genetic engineering.

Genome-wide identification and comparative analyses of key genes involved in C4 photosynthesis in five main gramineous crops

Compared to C₃ species, C₄ plants showed higher photosynthetic capacity as well as water and nitrogen use efficiency due to the presence of the C₄ photosynthetic pathway. Previous studies have shown that all genes required for the C₄ photosynthetic pathway exist in the genomes of C₃ species and are expressed. In this study, the genes encoding six key C₄ photosynthetic pathway enzymes (β-CA, PEPC, ME, MDH, RbcS, and PPDK) in the genomes of five important gramineous crops (C₄: maize, foxtail millet, and sorghum; C₃: rice and wheat) were systematically identified and compared. Based on sequence characteristics and evolutionary relationships, their C₄ functional gene copies were distinguished from non-photosynthetic functional gene copies. Furthermore, multiple sequence alignment revealed important sites affecting the activities of PEPC and RbcS between the C₃ and C₄ species. Comparisons of expression characteristics confirmed that the expression patterns of non-photosynthetic gene copies were relatively conserved among species, while C₄ gene copies in C₄ species acquired new tissue expression patterns during evolution. Additionally, multiple sequence features that may affect C₄ gene expression and subcellular localization were found in the coding and promoter regions. Our work emphasized the diversity of the evolution of different genes in the C₄ photosynthetic pathway and confirmed that the specific high expression in the leaf and appropriate intracellular distribution were the keys to the evolution of C₄ photosynthesis. The results of this study will help determine the evolutionary mechanism of the C₄ photosynthetic pathway in Gramineae and provide references for the transformation of C₄ photosynthetic pathways in wheat, rice, and other major C₃ cereal crops.

Genome-wide characterization, phylogenetic and expression analysis of Histone gene family in cucumber (Cucumis sativus L.)

Histones are essential components of chromatin and play an important role in regulating gene transcription and participating in DNA replication. Here, we performed a comprehensive analysis of this gene family. In this study, we identified 37 CsHistones that were classified into five groups (H1, H2A, H2B, H3 and H4). The closely linked subfamilies exhibited more similarity in terms of motifs and intron/exon numbers. Segmental duplication (SD) is the main driving force of cucumber CsHistones expansion. Analysis of cis-regulatory elements in the promoter region of CsHistones showed that CsHistones can respond to a variety of stresses. RNA-Seq analysis indicated that the expression of most CsHistones was associated with different stresses, including downy mildew, powdery mildew, wilt, heat, cold, salt stress, and waterlogging. Expression analysis showed that several genes of H3 group were highly expressed in different reproductive organs. Notably, CsCENH3 (CsHistone30) has the characteristics of a variant histone, and we demonstrated that CsCENH3 was localized on the nucleus and its proteins were expressed in centromere region. These findings provide valuable information for the identification and potential functions of Histone genes and ideas for the cultivation of CENH3-mediated haploid induction lines in cucumber.