DUF538 protein superfamily is predicted to be chlorophyll hydrolyzing enzymes in plants

Domain of unknown function (DUF) proteins in plants: function and perspective

Domains of unknown function (DUFs), which are deposited in the protein family database (Pfam), are protein domains with conserved amino acid sequences and uncharacterized functions. Proteins with the same DUF were classified as DUF families. Although DUF families are generally not essential for the survival of plants, they play roles in plant development and adaptation. Characterizing the functions of DUFs is important for deciphering biological puzzles. DUFs were generally studied through forward and reverse genetics. Some novelty approaches, especially the determination of crystal structures and interaction partners of the DUFs, should attract more attention. This review described the identification of DUF genes by genome-wide and transcriptome-wide analyses, summarized the function of DUF-containing proteins, and addressed the prospects for future studies in DUFs in plants.

Genome-wide association mapping combined with gene-based haplotype analysis identify a novel gene for shoot length in rice (Oryza sativa L.)

KEY MESSAGE

Genome-wide association mapping revealed a novel QTL for shoot length across multiple environments. Its causal gene, LOC_Os01g68500, was identified firstly through gene-based haplotype analysis, gene expression and knockout transgenic verification. Strong seedling vigor is an important breeding target for rice varieties used in direct seeding. Shoot length (SL) is one of the important traits associated with seedling vigor characterized by rapid growth of seedling, which enhance seedling emergence. Therefore, mining genes for SL and conducting molecular breeding help to develop varieties for direct seeding. However, few QTLs for SL have been fine mapped or cloned so far. In this study, a genome-wide association study of SL was performed in a diverse rice collection consisting of 391 accessions in two years, using phenotypes generated by different cultivation methods according to the production practice, and a total of twenty-four QTLs for SL were identified. Among them, the novel QTL qSL-1f on chromosome 1 could be stably detected across all three cultivation methods in the whole population and indica subpopulation. Through gene-based haplotype analysis of the annotated genes within the putative region of qSL-1f, and validated by gene expression and knockout transgenic experiments, LOC_Os01g68500 (i.e., Os01g0913100 in RAP-DB) was identified as the causal gene for SL, which has a single-base variation (C-to-A transversion) in its CDS region, resulting in the significant difference in SL of rice. LOC_Os01g68500 encodes a DUF538 (Domain of unknown function) containing protein, and the function of DUF538 protein gene on rice seedling growth is firstly reported in this study. These results provide a new clue for exploring the molecular mechanism regulating SL, and promising gene source for the molecular breeding in rice.

Antagonistic Regulation of ABA Responses by Duplicated Tandemly Repeated DUF538 Protein Genes in Arabidopsis

The plant hormone ABA (abscisic acid) regulates plant responses to abiotic stresses by regulating the expression of ABA response genes. However, the functions of a large portion of ABA response genes have remained unclear. We report in this study the identification of ASDs (ABA-inducible signal peptide-containing DUF538 proteins), a subgroup of DUF538 proteins with a signal peptide, as the regulators of plant responses to ABA in Arabidopsis. ASDs are encoded by four closely related DUF538 genes, with ASD1/ASD2 and ASD3/ASD4 being two pairs of duplicated tandemly repeated genes. The quantitative RT-PCR (qRT-PCR) results showed that the expression levels of ASDs increased significantly in response to ABA as well as NaCl and mannitol treatments, with the exception that the expression level of ASD2 remained largely unchanged in response to NaCl treatment. The results of Arabidopsis protoplast transient transfection assays showed that ASDs were localized on the plasma membrane and in the cytosol and nucleus. When recruited to the promoter of the reporter gene via a fused GD domain, ASDs were able to slightly repress the expression of the co-transfected reporter gene. Seed germination and cotyledon greening assays showed that ABA sensitivity was increased in the transgenic plants that were over-expressing ASD1 or ASD3 but decreased in the transgenic plants that were over-expressing ASD2 or ASD4. On the other hand, ABA sensitivity was increased in the CRISPR/Cas9 gene-edited asd2 single mutants but decreased in the asd3 single mutants. A transcriptome analysis showed that differentially expressed genes in the 35S:ASD2 transgenic plant seedlings were enriched in several different processes, including in plant growth and development, the secondary metabolism, and plant hormone signaling. In summary, our results show that ASDs are ABA response genes and that ASDs are involved in the regulation of plant responses to ABA in Arabidopsis; however, ASD1/ASD3 and ASD2/ASD4 have opposite functions.

Genome-Wide Analysis and Expression Profiling of DUF668 Genes in Glycine max under Salt Stress

The DUF668 gene performs a critical role in mitigating the impact of abiotic stress factors. In this study, we identified 30 DUF668 genes in a soybean genome, distributed across fifteen chromosomes. The phylogenetic analysis classified the DUF668 genes into three groups (group I, group II, and group III). Interestingly, gene structure analysis illustrated that several GmDUF668 genes were without introns. Furthermore, the subcellular localization results suggested that GmDUF668 proteins were present in the nucleus, mitochondria, cytoplasm, and plasma membrane. GmDUF668 promoters were analyzed in silico to gain insight into the presence of regulatory sequences for TFs binding. The expression profiling illustrated that GmDUF668 genes showed expression in leaves, roots, nodules, and flowers. To investigate their response to salt stress, we utilized the RNA sequencing data of GmDUF668 genes. The results unveiled that GmDUF668-8, GmDUF668-20, and GmDUF668-30 genes were upregulated against salt stress treatment. We further validated these findings using qRT-PCR analysis. These findings provide a scientific basis to explore the functions of GmDUF668 genes against different stress conditions.

PeCLH2 Gene Positively Regulate Salt Tolerance in Transgenic Populus alba × Populus glandulosa

Salt is an important environmental stress factor, which seriously affects the growth, development and distribution of plants. Chlorophyllase plays an important role in stress response. Nevertheless, little is known about the physiological and molecular mechanism of chlorophyll (Chlase, CLH) genes in plants. We cloned PeCLH2 from Populus euphratica and found that PeCLH2 was differentially expressed in different tissues, especially in the leaves of P. euphratica. To further study the role of PeCLH2 in salt tolerance, PeCLH2 overexpression and RNA interference transgenic lines were established in Populus alba × Populus glandulosa and used for salt stress treatment and physiologic indexes studies. Overexpressing lines significantly improved tolerance to salt treatment and reduced reactive oxygen species production. RNA interference lines showed the opposite. Transcriptome analysis was performed on leaves of control and transgenic lines under normal growth conditions and salt stress to predict genes regulated during salt stress. This provides a basis for elucidating the molecular regulation mechanism of PeCLH2 in response to salt stress and improving the tolerance of poplar under salt stress.

Unraveling the Diverse Roles of Neglected Genes Containing Domains of Unknown Function (DUFs): Progress and Perspective

Domain of unknown function (DUF) is a general term for many uncharacterized domains with two distinct features: relatively conservative amino acid sequence and unknown function of the domain. In the Pfam 35.0 database, 4795 (24%) gene families belong to the DUF type, yet, their functions remain to be explored. This review summarizes the characteristics of the DUF protein families and their functions in regulating plant growth and development, generating responses to biotic and abiotic stress, and other regulatory roles in plant life. Though very limited information is available about these proteins yet, by taking advantage of emerging omics and bioinformatic tools, functional studies of DUF proteins could be utilized in future molecular studies.

A Multi-Year, Multi-Cultivar Approach to Differential Expression Analysis of High- and Low-Protein Soybean (Glycine max)

Soybean (Glycine max (L.) Merr.) is among the most valuable crops based on its nutritious seed protein and oil. Protein quality, evaluated as the ratio of glycinin (11S) to β-conglycinin (7S), can play a role in food and feed quality. To help uncover the underlying differences between high and low protein soybean varieties, we performed differential expression analysis on high and low total protein soybean varieties and high and low 11S soybean varieties grown in four locations across Eastern and Western Canada over three years (2018-2020). Simultaneously, ten individual differential expression datasets for high vs. low total protein soybeans and ten individual differential expression datasets for high vs. low 11S soybeans were assessed, for a total of 20 datasets. The top 15 most upregulated and the 15 most downregulated genes were extracted from each differential expression dataset and cross-examination was conducted to create shortlists of the most consistently differentially expressed genes. Shortlisted genes were assessed for gene ontology to gain a global appreciation of the commonly differentially expressed genes. Genes with roles in the lipid metabolic pathway and carbohydrate metabolic pathway were differentially expressed in high total protein and high 11S soybeans in comparison to their low total protein and low 11S counterparts. Expression differences were consistent between East and West locations with the exception of one, Glyma.03G054100. These data are important for uncovering the genes and biological pathways responsible for the difference in seed protein between high and low total protein or 11S cultivars.

Genome-Wide Identification and Expression Profiling of DUF221 Gene Family Provides New Insights Into Abiotic Stress Responses in Potato

The domain of the unknown function 221 proteins regulate several processes in plants, including development, growth, hormone transduction mechanism, and abiotic stress response. Therefore, a comprehensive analysis of the potato genome was conducted to identify the deafness-dystonia peptide (DDP) proteins' role in potatoes. In the present study, we performed a genome-wide analysis of the potato domain of the unknown function 221 (DUF221) genes, including phylogenetic inferences, chromosomal locations, gene duplications, gene structures, and expression analysis. In our results, we identified 10 DDP genes in the potato genome. The phylogenetic analysis results indicated that StDDPs genes were distributed in all four clades, and clade IV was the largest clade. The gene duplication under selection pressure analysis indicated various positive and purifying selections in StDDP genes. The putative stu-miRNAs from different families targeting StDDPs were also predicted in the present study. Promoter regions of StDDP genes contain different cis-acting components involved in multiple stress responses, such as phytohormones and abiotic stress-responsive factors. The analysis of the tissue-specific expression profiling indicated the StDDPs gene expression in stem, root, and leaf tissues. We subsequently observed that StDDP4, StDDP5, and StDDP8 showed higher expressions in roots, stems, and leaves. StDDP5 exhibited high expression against heat stress response, and StDDP7 showed high transcript abundance against salt stress in potatoes. Under abscisic acid (ABA) and indole acetic acid (IAA) treatments, seven StDDP genes' expressions indicated that ABA and IAA performed important roles in immunity response. The expression profiling and real-time qPCR of stems, roots, and leaves revealed StDDPs' significant role in growth and development. These expression results of DDPs are primary functional analysis and present basic information for other economically important crops.

Overexpression of cassava RSZ21b enhances drought tolerance in Arabidopsis

Drought is one of the major environmental constraints affecting crop productivity. Plants have to adjust their developmental and physiological processes to cope with drought. We previously identified 18 cassava serine/arginine-rich (SR) proteins that had a pivotal role in alternative splicing in response to environmental stress. However, functional characterization of SR proteins is rarely explored. Here, we characterized the RSZ subfamily gene MeRSZ21b in cassava. The RSZ21b belongs to the RSZ subfamily, which was widely distributed in major crops and was highly conserved. Quantitative RT-PCR assay showed that the expression of MeRSZ21b was significantly induced by drought. Moreover, overexpression of MeRSZ21b in Arabidopsis was hypersensitive to abscisic acid (ABA) in the phases of seed germination and post-germination seedling growth. Meantime, MeRSZ21b overexpression lines were resistant to sorbitol treatment, and quickly closed the stomata when compared with Col-0 under drought condition. Importantly, overexpression of MeRSZ21b resulted in improved drought tolerance through modulating ABA-dependent signaling. Therefore, our findings refine our knowledge of the SR protein-coding genes and provide novel insights for enhancing plant resistance to environmental stress.

Chlorophyllase, a Common Plant Hydrolase Enzyme with a Long History, Is Still a Puzzle

Chlorophyllase (Chlase, CLH) is one of the earliest discovered enzymes present in plants and green algae. It was long considered to be the first enzyme involved in chlorophyll (Chl) degradation, while strong evidence showed that it is not involved in Chl breakdown during leaf senescence. On the other hand, it is possible that CLH is involved in Chl breakdown during fruit ripening. Recently, it was discovered that Arabidopsis CLH1 is located in developing chloroplasts but not in mature chloroplasts, and it plays a role in protecting young leaves from long-term photodamage by catalysing Chl turnover in the photosystem II (PSII) repair cycle. However, there remain other important questions related to CLH. In this article, we briefly reviewed the research progress on CLH and listed the main unanswered questions related to CLH for further study.

A pair of DUF538 domain-containing proteins modulates plant growth and trichome development through the transcriptional regulation of GLABRA1 in Arabidopsis thaliana

SMALLER TRICHOMES WITH VARIABLE BRANCHES (SVB) is an emerging plant growth regulator in trichome development, endoplasmic reticulum stress response, and phosphoinositide signaling, and belongs to the land plant-specific DUF538 domain-containing protein family. Despite its multifaceted roles, the functions of this protein family are poorly understood in plant growth and development. Here, we report that SVB-like (SVBL), the closest homolog of SVB, modulates plant growth and trichome development with SVB in Arabidopsis thaliana. Although none of the single mutants showed an obvious growth defect, the double mutants of svb svbl exhibited dwarfed plant growth. In trichome development, the defects in svb mutant were greatly enhanced by the additional mutation in SVBL, despite the single knockout of SVBL showing the mild defects. The double mutation reduced the transcript level of one of the central hub genes for trichome development, GLABRA1 (GL1), which in turn affects the other downstream genes, GLABRA2 (GL2), TRANSPARENT TESTA GLABRA2 (TTG2), TRIPTYCHON (TRY), CAPRICE (CPC), and ENHANCER OF TRY AND CPC1 (ETC1). In situ translational reporter assays showed that SVB and SVBL share highly similar localization patterns both at tissue and subcellular levels. The present study suggests that SVB and SVBL play a pivotal role in plant growth and trichome development by affecting a specific subset of known trichome developmental regulators, highlighting the importance of the DUF538 protein family in higher plants.

Genome-wide identification of the DUF668 gene family in cotton and expression profiling analysis of GhDUF668 in Gossypium hirsutum under adverse stress

Background

Domain of unknown function 668 (DUF668) may play a crucial role in the plant growth and developmental response to adverse stress. However, our knowledge of the function of the DUF668 gene family is limited.

Results

Our study was conducted based on the DUF668 gene family identified from cotton genome sequencing. Phylogenetic analysis showed that the DUF668 family genes can be classified into four subgroups in cotton. We identified 32 DUF668 genes, which are distributed on 17 chromosomes and most of them located in the nucleus of Gossypium hirsutum. Gene structure and motif analyses revealed that the members of the DUF668 gene family can be clustered in G. hirsutum into two broad groups, which are relatively evolutionarily conserved. Transcriptome data analysis showed that the GhDUF668 genes are differentially expressed in different tissues under various stresses (cold, heat, drought, salt, and Verticillium dahliae), and expression is generally increased in roots and stems. Promoter and expression analyses indicated that Gh_DUF668–05, Gh_DUF668–08, Gh_DUF668–11, Gh_DUF668–23 and Gh_DUF668–28 in G. hirsutum might have evolved resistance to adverse stress. Additionally, qRT-PCR revealed that these 5 genes in four cotton lines, KK1543 (drought resistant), Xinluzao 26 (drought sensitive), Zhongzhimian 2 (disease resistant) and Simian 3 (susceptible), under drought and Verticillium wilt stress were all significantly induced. Roots had the highest expression of these 5 genes before and after the treatment. Among them, the expression levels of Gh_DUF668–08 and Gh_DUF668–23 increased sharply at 6 h and reached a maximum at 12 h under biotic and abiotic stress, which showed that they might be involved in the process of adverse stress resistance in cotton.

Conclusion

The significant changes in GhDUF668 expression in the roots after adverse stress indicate that GhDUF668 is likely to increase plant resistance to stress. This study provides an important theoretical basis for further research on the function of the DUF668 gene family and the molecular mechanism of adverse stress resistance in cotton.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07716-w.

Identification of cassava alternative splicing-related genes and functional characterization of MeSCL30 involvement in drought stress

Alternative splicing (AS) is an important post-transcriptional regulation strategy that can increase the proteome diversity and regulate mRNA level in eukaryote. Multi-exon genes can be alternative spliced to generate two or more transcripts, thereby increasing the adaptation to the external stress conditions in planta. However, AS-related proteins were less explored in cassava which is an important staple crop in the tropical area. A total of 365 genes encoding AS-related proteins were identified and renamed in the cassava genome, and the transcriptional and splicing changes of 15 randomly selected genes were systematically investigated in the tissues under diverse abiotic stress conditions. 13 out of 15 genes undergo AS in the tissues and under diverse environmental stress condition. Importantly, the greatest changes of splicing patterns were found in the leaf or in response to temperature stress, indicating that AS-related proteins had their tissue-specific regulation patterns and might be participated in the plant adaptation to temperature stress. We then found that overexpression of MeSCL30 in Arabidopsis enhanced the tolerance to drought stress through maintaining reactive oxygen species (ROS) homeostasis and increasing the expression of drought-responsive genes. Therefore, these findings refined the AS-related protein-coding genes and provided novel insights for manipulation of AS-related genes in order to enhance the resistance to environmental stress in plant.

The DUF221 domain-containing (DDP) genes identification and expression analysis in tomato under abiotic and phytohormone stress

The domain of unknown function (DUF221 domain-containing) proteins regulates various aspects of plant growth, development, responses to abiotic stresses, and hormone transduction pathways. To understand the role of DDP proteins in tomato, a comprehensive genome-wide analysis was performed in the tomato genome. A total of 12 DDP genes were identified and distributed in 8 chromosomes in the tomato genome. Phylogenetically all SlDDPs were clustered into four clades, subsequently supported by their gene structure and conserved motifs distribution. The SlDDPs contained various cis-acting elements involved in plant responses to abiotic and various phytohormones stresses. The tissue-specific expression profile analysis revealed the constitutive expression of SlDDPs in roots, leaves, and developmental phases of fruit. It was found that SlDDP1, SlDDP3, SlDDP4, SlDDP9, SlDDP10, and SlDDP12 exhibited high expression levels in fruits at different development stages. Of these genes, SlDDP12 contained ethylene (ERE) responsive elements in their promoter regions, suggesting its role in ethylene-dependent fruit ripening. It was found that a single SlDDP induced by two or more abiotic and phytohormone stresses. These include, SlDDP1, SlDDP2, SlDDP3, SlDDP4, SlDDP7, SlDDP8, and SlDDP10 was induced under salt, drought, ABA, and IAA stresses. Moreover, tomato SlDDPs were targeted by multiple miRNA gene families as well. In conclusion, this study predicted that the putative DDP genes might help improve abiotic and phytohormone tolerance in plants, particularly tomato, rice, and other economically important crop plant species.

Cell Wall Proteome of Wheat Grain Endosperm and Outer Layers at Two Key Stages of Early Development

The cell wall is an important compartment in grain cells that fulfills both structural and functional roles. It has a dynamic structure that is constantly modified during development and in response to biotic and abiotic stresses. Non-structural cell wall proteins (CWPs) are key players in the remodeling of the cell wall during events that punctuate the plant life. Here, a subcellular and quantitative proteomic approach was carried out to identify CWPs possibly involved in changes in cell wall metabolism at two key stages of wheat grain development: the end of the cellularization step and the beginning of storage accumulation. Endosperm and outer layers of wheat grain were analyzed separately as they have different origins (maternal and seed) and functions in grains. Altogether, 734 proteins with predicted signal peptides were identified (CWPs). Functional annotation of CWPs pointed out a large number of proteins potentially involved in cell wall polysaccharide remodeling. In the grain outer layers, numerous proteins involved in cutin formation or lignin polymerization were found, while an unexpected abundance of proteins annotated as plant invertase/pectin methyl esterase inhibitors were identified in the endosperm. In addition, numerous CWPs were accumulating in the endosperm at the grain filling stage, thus revealing strong metabolic activities in the cell wall during endosperm cell differentiation, while protein accumulation was more intense at the earlier stage of development in outer layers. Altogether, our work gives important information on cell wall metabolism during early grain development in both parts of the grain, namely the endosperm and outer layers. The wheat cell wall proteome is the largest cell wall proteome of a monocot species found so far.

Characterization and Functional Divergence of a Novel DUF668 Gene Family in Rice Based on Comprehensive Expression Patterns

The domain of unknown function (DUF) superfamily encodes proteins of unknown functions in plants. Among them, DUF668 family members in plants possess a 29 amino-acid conserved domain, and this family has not been described previously. Here, we report this plant-specific novel DUF668 gene family containing 12 OsDUF668 genes in rice (Oryza sativa) and 91 DUF668s for the other seven plant species. In our study, DUF668 genes were present in both dicot and monocot plants, indicating that DUF668 is a conserved gene family that originated by predating the dicot–monocot divergence. Based on the gene structure and motif composition, the DUF668 family consists of two distinct clades, I and II in the phylogenetic tree. Remarkably, OsDUF668 genes clustered on the chromosomes merely show close phylogenetic relationships, suggesting that gene duplications or collinearity seldom happened. Cis-elements prediction display that over 80% of DUF668s contain phytohormone and light responsiveness factors. Further comprehensive experimental analyses of the OsDUF668 family are implemented in 22 different tissues, five hormone treatments, seven environmental factor stresses, and two pathogen-defense related stresses. The OsDUF668 genes express ubiquitously in analyzed rice tissues, and seven genes show tissue-specific high expression profiles. All OsDUF668s respond to drought, and some of Avr9/Cf-9 rapidly elicited genes resist to salt, wound, and rice blast with rapidly altered expression patterns. These findings imply that OsDUF668 is essential for drought-enduring and plant defense. Together, our results bring the important role of the DUF668 gene family in rice development and fitness to the fore.

Maize glossy6 is involved in cuticular wax deposition and drought tolerance

Cuticular waxes, long-chain hydrocarbon compounds, form the outermost layer of plant surfaces in most terrestrial plants. The presence of cuticular waxes protects plants from water loss and other environmental stresses. Cloning and characterization of genes involved in the regulation, biosynthesis, and extracellular transport of cuticular waxes onto the surface of epidermal cells have revealed the molecular basis of cuticular wax accumulation. However, intracellular trafficking of synthesized waxes to the plasma membrane for cellular secretion is poorly understood. Here, we characterized a maize glossy (gl6) mutant that exhibited decreased epicuticular wax load, increased cuticle permeability, and reduced seedling drought tolerance relative to wild-type. We combined an RNA-sequencing-based mapping approach (BSR-Seq) and chromosome walking to identify the gl6 candidate gene, which was confirmed via the analysis of multiple independent mutant alleles. The gl6 gene represents a novel maize glossy gene containing a conserved, but uncharacterized, DUF538 domain. This study suggests that the GL6 protein may be involved in the intracellular trafficking of cuticular waxes, opening the door to elucidating the poorly understood process by which cuticular wax is transported from its site of biosynthesis to the plasma membrane.