1. Introduction
All living organisms are exposed to unfavorable stimuli for shorter or longer periods of time. Unlike animals, plants are stationary life forms that cannot escape stressful factors by moving away. Nonetheless, plants successfully survive because they have developed numerous defense mechanisms and tolerance to environmental stresses [
1]. Any harmful biotic or abiotic factor has the potential to induce oxidative stress in plant cells [
2]. Conditions for oxidative stress occur when the cell enters a state of imbalance between the production of reactive oxygen species (ROS) and its ability to detoxify these reactive intermediates. Highly reactive molecules can cause serious damage to biomolecules such as proteins, lipids, and DNA, leading to a disturbance in cellular functions and even cell death [
3]. Plants respond to oxidative stress through a complex antioxidant defense system that includes both enzymatic and non-enzymatic mechanisms. Superoxide dismutase, catalase, and ascorbate peroxidase are enzymatic antioxidants responsible for transforming ROS into less harmful molecules [
4]. Non-enzymatic antioxidants are different compounds, such as ascorbate, glutathione, and tocopherol, that have the ability to neutralize ROS [
3]. At the foundation of plant defense against oxidative stress lie intricate signaling networks. These networks are pivotal in regulating the production of defensive secondary metabolites and coordinating alterations in gene expression. Understanding these mechanisms is important for developing strategies to enhance plant resistance to stress. New insights can enhance crop yields and overcome agriculture issues arising from climate change.
Along with the well-known components of the antioxidant system, members of a large family of naturally disordered proteins (IDPs—intrinsically disordered proteins) are also essential for plants to react effectively to stress. IDPs are highly represented in the proteome of eukaryotes. Plant IDPs play a key role in plant biology, acting as a connection between different intracellular regulatory signals and signals from the external environment. Research has shown that these proteins play a significant role in abiotic stress, transcriptional regulation, light perception, and plant development [
5]. Late embryogenesis abundant (LEA) proteins, gibberellic acid-insensitive repressor of ga1-3 scarecrow (GRAS) proteins, and cryptochromes (CRYs) are well-known protein families in plants that link intracellular regulatory processes with external environmental signals [
5]. Among these families, LEA proteins are the most prominent and are typically induced by various environmental signals such as cold, frost, heat, drought, salinity, and increased evaporation [
6]. These proteins serve diverse functions. For instance, LEA proteins can interact with partially denatured polypeptides, preventing protein aggregation in the cytoplasm during dehydration. Additionally, they can serve as both chaperones and transcription factors [
7]. As transcription factors, LEA proteins connect external signals to specific gene expression patterns. One IDP, specifically the DSS1 (deletion of split hand/split foot 1) protein, could be one of the many factors involved in the plant stress response, operating within the proteasome machinery alongside various other factors. It was shown that the absence of DSS1 proteins in
Arabidopsis thaliana mutant lines exposed to hydrogen peroxide (H
2O
2) leads to the accumulation of oxidized proteins and disruption of protein homeostasis [
8].
DSS1 is a small, acidic, and highly conserved eukaryotic protein. It belongs to a large family of naturally disordered proteins lacking three-dimensional structure. DSS1 proteins also lack well-defined secondary elements. Even when DSS1 is bound to other molecules, large regions of the protein remain disordered, and only a limited secondary helical structure is observed at the C-terminus of the protein [
9,
10,
11]. Despite the absence of a 3D structure, IDPs participate in cell functioning and help maintain the stability of biological complexes, primarily as regulatory components [
12]. DSS1 is a multifunctional protein that participates in various biological processes through its interaction with various components of numerous protein complexes. When this flexible protein binds to partner proteins, it acquires final 3D conformation. The molecular weights of DSS1 proteins range from 7–9 kDa, with lengths of polypeptide chains ranging from 70 to 90 amino acids in different organisms. Multiple alignments of amino acid sequences of DSS1 orthologues exhibit high levels of homology between fungi, plants, animals, and humans [
11]. The fully conserved region of 15 amino acids presents a functionally significant domain for interaction with other proteins [
13]. The secondary structure of DSS1 contains a transient α-helix at the end of the C-terminus and two conserved acidic regions [
13,
14].
One of the first well-defined roles of DSS1 is its involvement in the regulation of homologous recombination (HR), the most precise mechanism of DNA double-strand break repair [
15]. In most eukaryotes, the core of HR machinery is comprised of recombinase RAD51, its mediator BRCA2, and regulator DSS1. The occurrence of DNA double-strand breaks induces the recruitment of a number of proteins that process the breaks, generating long 3′ DNA tails that are then coated with RPA homofilament. Recombinase RAD51 uses this single-stranded DNA to search for and invade homologous DNA. However, displacement of RPA by RAD51 requires the BRCA2 protein, which is conserved in higher eukaryotes and also found in fungi from the Basidiomycota division, such as
Ustilago maydis [
16]. The interplay of DSS1 and BRCA2 is essential for the properly regulated HR activity of the BRCA2 machinery, contributing to the maintenance of genome integrity [
13,
14,
15,
16,
17,
18,
19,
20,
21].
The suppressor of exocyst mutation 1, Sem1 (or Dss1), has been hypothesized to take part in two structurally related complexes with distinct functions in mRNA processing in yeast. One of these is the ternary complex Thp3-Csn12-Sem1, which regulates transcription and facilitates pre-mRNA splicing [
14,
22,
23,
24]. In another larger complex in which Sem1/Dss1 contributes to the assembly of the nuclear pore (NCP)—transcription-export complex-2 (TREX-2), it plays a role in mRNA export.
As the smallest regulatory and structural subunit of the 26S proteasome, besides ubiquitin-binding activity, Dss1 also participates in the correct proteasome assembly. Dss1 contributes to the stability of the proteasome by effectively recruiting its subunits Rpn3 and Rpn7, which normally have a low affinity for each other [
23,
25,
26]. Although complete assembly of the proteasome is possible in the absence of Dss1, the proteolytic activity of this proteasome is incomplete. During proteolysis, Dss1 is necessary for the fine modulation of Rpn7 in the process of efficient ATP-dependent substrate unwinding [
25].
Our studies have revealed that the
A. thaliana dss1(V) mutant has increased sensitivity to oxidative stress when compared to
dss1(I) or wild-type (WT) plants. The lack of DSS1(V) leads to an excessive accumulation of oxidized proteins, likely caused by the impaired functioning of the 26S proteasome. Conversely, the absence of DSS1(I) in
dss1(I) mutant plants caused only moderate accumulation of damaged proteins under H
2O
2-induced oxidative stress [
8].
In this paper, we further explored the potential differences in the functions of two highly homologous DSS1 proteins using knockdown dss1 mutant plants. Transcriptome analysis revealed that the absence of either DSS1 gene caused alterations in the expression of different genes under stress conditions. In terms of transcriptional regulation, major differences between the DSS1(I) and DSS1(V) genes were not detected. However, it is important to note that the lack of significant differences does not rule out subtle variations in their regulatory functions or interactions with other cellular components. Our results imply that DSS1s act as negative regulators of the genes involved in the synthesis of secondary metabolites that play a vital role in plant defense against stress. Based on the results of the functional complementation test, we showed that the DSS1(I) gene may make a more substantial contribution to the HR process. Further research on the interactions of individual plant DSS1 isoforms with other proteins involved in DNA repair or studies on their interaction with DNA regulatory sequences of genes that showed altered expression in transcriptome analysis could reveal their specific contributions to cellular processes.
3. Discussion
Research on
AtDSS1 genes points to the possibility of functional divergence between two highly homologous
DSS1 genes as regards their differential contributions to plant phenotypes and developmental dynamics [
8,
28]. Also, we demonstrated that
dss1(I) and
dss1(V) mutant plants have different sensitivities to oxidative stress compared to WT plants and to each other. This can be attributed to the different roles of DSS1(I) and DSS1(V) in maintaining protein homeostasis [
8]. While DSS1 proteins exhibit high similarity, the Q44L substitution was hypothesized to induce a discernible variance between the two protein isoforms, which could also influence the observed phenotypic differences. These DSS1 proteins could have different functions in the 26S proteasome, including the elimination of damaged proteins. The absence of DSS1(V) protein leads to a significant accumulation of damaged proteins due to the inefficient removal of oxidized proteins, while the absence of DSS1(I) protein moderately influences oxidized protein accumulation [
8]. Current research based on the functional complementation test and comparison of mutant and WT transcriptomes further confirms the differences between the two
A. thaliana DSS1 genes.
Using a semiquantitative approach, the functional complementation test demonstrated that the DSS1(I) protein of
A. thaliana can restore the WT phenotype of the
U. maydis Δdss1 mutant. Overexpression of
A. thaliana DSS1(I) in the
Δdss1 mutant results in significant resistance to sublethal doses of genotoxic agents. The survival rate of transformants
Δdss1 + AtDSS1(I) is 1000 times higher than that of
Δdss1, but lower than the rate achieved in transformants
Δdss1 + UmDss1. DSS1 and BRCA2 are highly conserved proteins that interact via helix domains (HD). Literature data indicate a high degree of similarity in the HD domain, as 38% similarity is noted between
U. maydis (Brh2) and human BRCA2, and 44% similarity between
A. thaliana and human BRCA2 [
15]. Although UmDss1 and AtDSS1 differ in their amino acid sequences, the regions responsible for interactions with partner proteins are highly conserved [
15]. According to the presented data, we hypothesize that besides the differences in the amino acid sequences of the protein, DSS1(I) successfully interacts with the Brh2 complex of
U. maydis due to evolutionary conserved domains that interact with each other. This is supported by the observation that in
A. thaliana, DSS1(I) can bind both AtBRCA2(IV) and AtBRCA2(V) isoforms [
29]. This result indicates that two variants of the
A. thaliana DSS1 protein have different roles in recombinant DNA repair lesions. However, the absence of functional complementation of the
U. maydis Δdss1 by
A. thaliana DSS1(V) can be attributed to the distinct structural organization of these two
A. thaliana DSS1 proteins, which are linked to their primary amino acid sequences. Genetic redundancy in plants is the result of polyploidy, and this phenomenon allows for the accumulation of mutations, leading to the appearance of new alleles and the expansion of the gene family [
30]. The evolutionary model of neofunctionalization suggests that due to gene duplication, the “ancestral” gene retains its functions while its paralog acquires new ones [
31]. Given the assumption that the
AtDSS1(V) gene is the result of gene duplication, as we discussed in our previous publication [
8], and considering new results from the functional complementation tests and transcriptome analysis of the mutant lines, it is plausible to suggest that changes in protein structure during evolution contribute to the appearance of functional differences among
DSS1 genes.
To explore the potential differences in functions and uncover novel roles of the DSS1(I) and DSS1(V) genes in plant defense against oxidative stress, we analyzed the transcriptomes of dss1 mutant seedlings and compared them not only to each other but also to the WT. Significant differences in gene expression between the CRISPR/Cas9 dss1 mutants and WT seedlings under oxidative stress were observed. Genes showing more significant expression changes in the mutated lines compared to WT plants were subsequently classified into metabolic pathways using KEGG bioinformatics analysis.
Dss1 mutant transcriptome analysis revealed that the biosynthesis of a broad range of protective compounds with antioxidant properties could be interrupted by their decreased gene expression due to the absence of DSS1 protein. However, there is still a high similarity between the transcriptomes of
dss1(I) and
dss1(V) plants regarding these metabolic pathways. Reduced expression of genes involved in the biosynthesis of secondary metabolites such as flavonoids, polyphenolic compounds, and phenylpropanoids was noted in both mutants; thus, it can be assumed that DSS1 proteins have an important role in the transcriptional regulation of the genes involved in their biosynthesis. Generally, these secondary metabolites help plants remove cytotoxic products and ROS molecules [
32]. Although numerous transcriptional regulators of secondary metabolite biosynthesis have been identified and characterized, a lack of information hinders our comprehensive understanding of metabolic pathway regulation in response to environmental signals [
33]. If DSS1 proteins are involved in transcriptional complexes, they could potentially play a role in regulating the synthesis of secondary metabolites through transcriptional regulation [
14,
23]. However, during transcriptome analysis of the mutants, it was found that some genes involved in the synthesis of phenylpropanoids have elevated expression in
dss1(I) mutants compared to the
dss1(V) line. This points to an unambiguous function of the DSS1(V) protein in response to oxidative stress. DSS1(V) could be a modulator of the protein network involved in initiating the expression of enzymes involved in the biosynthesis of phenylpropanoids. We assume that the absence of DSS1 proteins could make the transcriptional machinery less stable. Reducing the expression of genes that encode enzymes involved in the biosynthesis of protective compounds could render mutated plants more sensitive to oxidative stress. Additionally, we confirmed that
CYP75B and
LDOX are notably downregulated only in
dss1(I) mutants. CYP75B contributes to the synthesis of plant antioxidants and thus participates in stress tolerance mechanisms, helping plants adapt to adverse conditions [
34]. LDOX encodes a 2-oxoglutarate-dependent dioxygenase called anthocyanidin synthase, which converts leucoanthocyanidins to anthocyanidins, the precursors of anthocyanins. This enzyme is essential to the production of anthocyanins, which are also responsible for plant protection against UV radiation and cold stress [
35,
36]. In this paper, we also demonstrated that
dss1(V) mutant plants have lower levels of AOCS6,
PR1, and
BGLU28 transcripts when exposed to oxidative stress. Both
BGLU28 and
PR1 genes are involved in the plant defense response against biotic and abiotic stresses. BGLU28 is a beta-glucosidase enzyme that hydrolyzes glucosinolates, a class of sulfur-containing compounds that have anti-microbial and anti-herbivore activities [
37]. PR1 is a pathogenesis-related protein that is induced by salicylic acid and confers resistance to fungal and bacterial pathogens [
38], while AOCS6 is a MAPK substrate and a participant in ethylene biosynthesis [
39]. DSS1(V) could negatively regulate various genes involved in plant protection against stress. This explains our previous results indicating that the
dss1(V) mutant is more sensitive to stress than the
dss1(I). A significant difference in expression was observed in the expression profile of the
PRX70 gene among the genes involved in the process of biogenesis of secondary metabolites. It is considerably downregulated in
dss1(V) mutants. According to different studies, this gene is involved in the removal of H
2O
2, the biosynthesis and degradation of lignin, the response to environmental stresses, pathogen attack, and oxidative stress [
40]. However, it remains unproven whether the regulation of PRX70 enzyme transcription could be directly influenced by DSS1(V).
Transcriptome analysis revealed that the disruption of
DSS1(V) in
dss1(V) mutants leads to decreased expression of genes involved in the biosynthesis of plant hormones when compared to WT plants. Plant hormones also play a crucial role in plant defense against oxidative stress. Disruption of biosynthetic pathways due to DSS1 absence could cause mutant sensitivity to oxidative stress. In silico analysis of the promoter region of the
DSS1(V) gene revealed the presence of
cis-activating regulatory elements involved in the response to phytohormones such as jasmonates and salicylic acid [
28]. In the complex and dynamic processes of transcriptional regulation in plant stress response, DSS1(V) and phytohormones could regulate each other’s expression through positive feedback mechanisms.
In comparison to WT seedlings,
dss1(I) and
dss1(V) mutants showed reduced expression of the genes involved in the MAPK signaling pathway under oxidative stress. MAPK has different roles in cell signal transmission, especially in the transduction of extracellular signals such as biotic and abiotic stresses, as well as in the transduction of stimuli in processes during development, differentiation, proliferation, and cell death [
41]. Participants in the MAPK signaling pathway with decreased levels of transcripts also have a significant role in defense against pathogens and ROS molecules, and it is not surprising that
dss1(I) and
dss1(V) mutants have increased sensitivity to oxidative stress [
42,
43]. The downregulation of 17 genes within the MAPK cascade pathway in plants lacking the DSS1 gene may certainly disrupt plant immunity and impact their stress response. However, the specific level and mechanism by which DSS1 influences the transcriptional regulation of these genes remain unclear.
Reduced expression of certain participants in glutathione metabolism was also observed in both
dss1 mutant lines. Glutathione is a tripeptide with many vital functions, one of which is the detoxification of toxic metabolites produced under oxidative stress. In plants, the metabolism of GSH involves enzymes for the biosynthesis and degradation of GSH, as well as enzymes that utilize GSH as a substrate [
44]. According to transcriptome analysis, we hypothesize that the lack of DSS1 protein also inhibits the expression of the genes involved in glutathione metabolism. It is conceivable that this inhibitory effect also contributes to the mutated
dss1 line’s inability to overcome the oxidative imbalance. On the whole, perturbations in MAPK signaling and glutathione metabolism induced by the absence of
DSS1 genes could have broader impacts on plant stress responses, cellular health, and overall adaptation to environmental fluctuations.
In contrast, both
dss1 mutant lines show increased expression of many genes involved in HR. DSS1 is one of the components of the BRCA2 complex and is a part of the HR machinery [
22,
29]. In this paper, several HR genes that exhibited the highest expression were considered for validation. Interestingly, the HR genes,
GIS5,
RECG,
RPA70,
RECQ4A, and
EME1B, were upregulated in both mutants as compared to WT plants, while
BRCA1,
CLSY3, and
RECQ4B showed increased expression in the transcriptome and real-time analysis, but only in the
dss1(I) mutant. The GIS5 protein has a catalytic role within DNA polymerase δ, which participates in DNA replication and repair [
45]. BRCA1, RECG, RPA70, and EME1B are involved in DNA repair [
46,
47,
48,
49], CLSY3 has ATP-dependent chromatin remodeling activity [
50], while RECQ4A and RECQ4B have ATP-dependent DNA unwinding activity [
51] (
Table S2). DSS1(I) may represent a specific negative regulator of the transcription of these genes, considering that the absence of DSS1(I) leads to elevated levels of these genes. Although it has not been experimentally demonstrated so far, DSS1(I) may be an unusual member of the complexes along with BRCA1, CLSY3, or RECQ4B. Existing literature provides evidence that proteins belonging to the same complex could be co-regulated at the transcriptional level [
52,
53]. We hypothesize that the expression of DSS1 and other protein partners involved in homologous recombination may likewise be co-regulated. Consequently, if one protein is absent—as in the case of the
dss1 mutants, where DSS1 is absent—the cell may compensate for the deficiency through negative feedback. The compensatory mechanism may result in the up-regulation of other co-regulated protein partners within the complex, potentially explaining the elevated expression of the analyzed homologous recombination genes.
Dss1(I) mutants showed significantly elevated expression of genes that encode ABC transporters. ABC transporters are the largest family of proteins in all living organisms [
45]. Compared to animals, the plant genome contains twice the number of ABC proteins. These proteins play a role in DNA repair, RNA translocation, and the active transport of numerous compounds through different types of membranes [
54]. What has been established to date is that DSS1 has a potential role in proteasome regulation and RNA export. Under conditions of oxidative stress,
dss1(I) mutant plants experience disruptions in the maintenance of cellular homeostasis. This could trigger enhanced expression of several ABC transporters to accelerate RNA transport intracellularly and eliminate toxic molecules extracellularly [
54]. However, more research is needed to fully understand the molecular mechanisms and biological functions of this interplay between DSS1 proteins and ABC transporters.
The characterization of
dss1(I) or
dss1(V) mutant transcriptomes reveals a new potential role for DSS1 in either activating or repressing numerous genes, leading to differential gene expression in mutant lines compared to the WT. Much the same, many other intrinsically disordered proteins, DSS1 proteins have ubiquitous and diverse functions in transcriptional regulation and cell signaling [
55,
56]. These results and our previous findings suggest that DSS1(I) and DSS1(V) proteins play a role in various cellular processes. DSS1(I) has a role in HR, as demonstrated by the ability of AtDSS(I) to reconstitute the wild-type DNA repair phenotype of the
U. maydis Δdss1 mutant. In the absence of DSS1(I), plants showed slower developmental dynamics and the occurrence of abortive seeds. DSS1(V) provides a greater contribution to plant defense against oxidative stress and maintenance of protein homeostasis, which is evidenced by the increased sensitivity of
dss1(V) mutant plants to oxidative stress as well as the increased accumulation of oxidized proteins compared to WT plants, attributed to impaired function of the 26S proteasome [
8,
28]. These findings underscore the significance of understanding the molecular mechanisms governing plant responses to oxidative stress. Understanding the role of the DSS1 protein in plant defense mechanisms could be vital for crops, especially given the prevalent factors of climate change and oxidative stress.
4. Materials and Methods
4.1. Plant Material, Cultivation, and Treatment
WT plants (Columbia-0),
dss1(I), and
dss1(V) mutants of
Arabidopsis thaliana were grown in sterile Murashige and Skoog (MS) agar medium containing 1% sucrose [
57]. The homozygous mutant lines with large indels (
dss1(I)del25 and
dss1(V)ins18) were obtained using CRISPR/Cas9 technology. The CRISPR/Cas9 gene editing technique, Arabidopsis transformation, and mutant screening have been described previously [
8]. Fourteen-day-old Arabidopsis seedlings were grown for three days on agar plates containing 10 mM hydrogen peroxide. On the third day of treatment, the seedlings were harvested, immediately placed, and stored at −80 °C.
4.2. Ustilago Maydis Strains, Methods, and Stress Treatment
Manipulations with
Ustilago maydis, culture methods, transformations, gene transfer, treatments, and survival after DNA damage were carried out as described [
27,
58].
U. maydis’ UCM350 and Δ
dss1 strains were used for the functional complementation test. UCM350 is the nominal WT strain (
nar1-6 pan1-1 a1b1); the Δ
dss1 strain was generated by disruption of the
UmDss1 gene in UCM350 [
15,
27].
The self-replicating pCM955 vector was used to express plant AtDSS1 genes in the U. maydis Δdss1 mutant. pCM955 contains a promoter for glyceraldehyde 3-phosphate dehydrogenase (gap) and a positive selectable marker conferring resistance to hygromycin (the HPH cassette). Transformation of the U. maydis Δdss1 mutant was also carried out with the pCM1019 vector, which contains the HPH marker and the complementary gene for UmDss1, driven by the gap promoter.
The degree of resistance of U. maydis to genotoxic agents was determined semi-quantitatively. Cell counting from the overnight culture was done under a microscope (BA80, Motic, Xiamen, China) using a hemocytometer. U. maydis cells were grown in a-liquid rich YEPS medium (1% yeast extract, 2% sucrose, and 2% peptone) at 30 °C with agitation at 200 rpm. Spot assays were conducted by making five serial 10-fold dilutions from the initial cell suspension of 2 × 107/mL. Then, 10 μL portions of each dilution were placed in order from left to right onto solid YEPS medium. Plates were incubated for three days at 30 °C for colonies to develop. To test the resistance of cells to genotoxic agents, the medium contained either 0.015% MMS or 0.0075% DEB, or an open plate was irradiated with a 254 nm germicidal UV lamp (OSRAM, Germicidal lamp, Puritec 30 W, G13 G30T8/OF, Moscow, Russia) at a dose rate of 100 J/m2.
4.3. The Functional Complementation Test of U. maydis Δdss1 Mutant
For the functional complementation test, we conducted cloning and subcloning series to analyze the genetic complementation of the
U. maydis Δ
dss1 mutant sensitive to UV stress with the AtDSS1(I) and AtDSS1(V) proteins of
Arabidopsis thaliana. The open reading frames of the plant
AtDSS1(I) and
AtDSS1(V) genes were amplified from the total cDNA of
Arabidopsis thaliana using gene-specific primers through PCR and then subcloned into two separate pCRII plasmids (Invitrogen™, Thermo Fisher Scientific™, Waltham, MA, USA) (
Table S1). cDNA was obtained from total RNA extracted from wild-type seedlings of
A. thaliana grown on plates. The
AtDSS1(I) cDNA sequence was amplified using the UmNdeI_F/3UTR DI_R primers, and the
AtDSS1(V) cDNA sequence was amplified using the UmNdeI_F/3UTR DV_R primers. The recognition site for the restriction enzyme NdeI was engineered upstream of the translation start site of the
AtDSS1 gene using the primer UmNdeI_F. The PCR product was ligated into the pCRII vector using the TA Cloning™ Kit, Dual Promoter Kit (Invitrogen™, Thermo Fisher Scientific™, Waltham, MA, USA), following the manufacturer’s protocol. Electrocompetent cells MegaX DH10B T1R, transformed by the ligation mixture pCRII-
AtDSS1(I) and pCRII-
AtDSS1(V), were selected in a medium with 50 μg/mL kanamycin.
pCM1027, pCRII-AtDSS1(I), and pCRII-AtDSS1(V) vectors were digested with NdeI and HindIII (Thermo Scientific Scientific™, Waltham, MA, USA). The digestion fragments were isolated by gel purification using the PureLink™ Quick Gel Extraction Kit (Invitrogen™, Thermo Fisher Scientific™, Waltham, MA, USA) according to the manufacturer’s instructions. Fragments with AtDSS1(I) or AtDSS1(V) were inserted into the linearized pCM1027 vector downstream of the gap promoter for expression in U. maydis. The transformation of the DH10B E. coli strain was performed with ligation products, and the transformants were selected in a medium containing 100 μg/mL ampicillin.
pCM955, pCM1027-AtDSS1(I), and pCM1027-AtDSS1(V) vectors were digested with Acc65I and BamHI (Thermo Scientific™, Waltham, MA, USA). The fragments gap + DSS1(I) and gap + DSS1(V) were purified and ligated into the linearized expression pCM955 shuttle vector. The success of every cloning was validated by sequencing using universal M13 primers as well as verifying the accuracy of the gene sequences (Macrogen Europe BV, Amsterdam, The Netherlands).
4.4. RNA Extraction and cDNA Synthesis
For the analysis of the AtDSS1 transcripts in U. maydis strains, UCM350 and Δdss1 strains were grown overnight, and cells were washed twice. The cells were collected, and the total RNA of each sample was extracted using the GeneJET RNA Purification Kit (Thermo Scientific Scientific™, Waltham, MA, USA) following the manufacturer’s protocol with some modifications. To increase the amount of extracted RNA, the collected cells were resuspended in yeast lysis buffer (1 M sorbitol, 0.1 M EDTA, pH 7.4) containing 20 mM dithiothreitol (DTT) and lysis enzyme from Trichoderma harzianum (Sigma, St. Louis, MO, USA). Suspensions were incubated for 30 min at 30 °C with agitation. Further steps were performed according to the manufacturer’s instructions. RNA purity was checked using the BioSpec-nano spectrophotometer (Shimadzu, Kyoto, Japan). Total RNA was treated with the DNA-free™ DNA Removal Kit (Invitrogen™, Thermo Fisher Scientific™, Waltham, MA, USA) to remove DNA contamination before cDNA synthesis. For cDNA synthesis, reverse transcription was performed using the RevertAid First Strand cDNA Synthesis Kit (Thermo Fisher Scientific, Waltham, MA, USA).
The examination of the transcriptome in this study involved the use of previously described 14-day-old
dss1 mutants and WT seedlings of
A. thaliana, from which total RNA was extracted and cDNA synthesized by reverse transcription as described previously [
8].
4.5. PCR and Gene Expression
The primer pairs used in this study are listed in
Table S1. They were designed using Primer3web (accessed in August and September 2023). The primer pairs RECG and BRCA1 utilized herein were sourced from previous research [
47,
59]. Each 25 μL PCR reaction contained 12.5-μL of the DreamTaq Green PCR Master Mix (2×) (Thermo Fisher Scientific™, Waltham, MA, USA), 10 ng of the template, and 200nM of the primer set (
Table S1). PCR reactions were incubated in the Tprofessional Thermocycler (Biometra, Jena, Germany) under the following cycling program: 5 min at 95 °C, 35 cycles at 95 °C for 30 s, at 62 °C for 30 s, at 72 °C for 30 s, and for 5 min at 72 °C. The PCR products were resolved on 2% agarose gels, stained with ethidium bromide, and visualized under UV light.
For gene expression analysis, before the SYBR Green assay, total cDNAs were diluted 1:16 with nuclease-free water. Reactions were performed in 25 µL containing 300 nM of each primer (
Table S1) and Power SYBR™ Green PCR Master Mix (Applied Biosystems™, Thermo Fisher Scientific™, Waltham, MA, USA). Real-time PCR was conducted on the Magnetic Induction Cycler (Mic qPCR, Bio Molecular Systems, Queensland, Australia) under the following cycles: 10 min at 95° C, and 40 cycles at 95° C for 15 s, and 60°C for 1 min. Each PCR reaction was performed in triplicate, and no template controls were included. Amplification of PCR products was detected in real-time, and
Arabidopsis ACT2 (AT3G18780) was used as an internal normalization control. The results were analyzed with MicPCR software v2.12.7 (Bio Molecular Systems, Brisbane, Queensland, Australia) and presented as 2
−ΔCt.
4.6. RNA Library Construction and Sequencing
To carry out RNA sequencing of three different biological samples (WT,
dss1(I), and
dss1(V) mutants under oxidative stress), 400 ng of the isolated RNA per sample were sent to Novogene Bioinformatics Technology Co., Ltd. (Cambridge, UK) at room temperature stabilized by the DNA/RNA shield buffer (Zymo Research, Irvine, CA, USA). RNA material from each sample was collected from three biologically independent experiments [
8].
Total RNA was used as input material for the RNA sample preparations. RNA purity, concentration, and quality were assessed using the RNA Nano 6000 Assay Kit of the Bioanalyzer 2100 system (Agilent Technologies, Palo Alto, CA, USA). After poly-A-based mRNA enrichment and cDNA library preparation, 150-bp pair-end RNA sequencing was conducted on the Illumina NovaSeq 6000 platform (Illumina, Inc., San Diego, CA, USA).
4.7. Bioinformatic Analysis
Basic data analysis with a well-annotated reference genome was carried out by Novogene (Cambridge, UK). The initial step involved processing the raw reads in FASTQ format, followed by the calculation of Q20, Q30, and the GC content for the clean reads. The raw reads were subjected to filtering, excluding any reads with adapter contamination, reads containing more than 10% of uncertain nucleotides, and low-quality reads (with a Qscore < 5 for more than 50% of the bases) using the FASTP tool (
https://github.com/OpenGene/fastp). All the downstream analyses were based on clean, high-quality data.
Mapping paired-end clean reads to the Arabidopsis thaliana reference genome (ensemblplants_arabidopsis_thaliana_tair10_gca_000001735_1) was accomplished by HISAT2 (v2.0.5). featureCounts (v1.5.0-p3) was used to count the number of reads mapped to each gene. The FPKM (Fragments Per Kilobase of transcript per Million mapped reads) of each gene was calculated based on the length of the gene and the read count mapped to this gene.
Before analyzing differential gene expression, the read counts for each sequenced library were adjusted using the edgeR program package by applying a scaling normalization factor. Differential expression analysis of two groups was performed using the edgeR (without biological replicates) R package (v3.22.5). The p-values were adjusted using the Benjamini-Hochberg method to control the false discovery rate. Genes with a corrected p-value (padj) < 0.05 and an absolute fold change of 2 (|log2(Fold Change)| ≥ 1) identified by edgeR were classified as differentially expressed.
Enrichment analysis was done using the clusterProfiler (v3.8.1) R package. R package (with correction of gene length bias and threshold of corrected
p-value < 0.05) was considered to test the statistical enrichment of differentially expressed genes in GO terms (
http://www.geneontology.org/; accessed on December 2022) and KEGG pathways (
http://www.genome.jp/kegg/; accessed on December 2022). The data discussed in this publication have been deposited in the NCBI Gene Expression Omnibus [
60] and can be accessed through GEO Series accession number GSE264581 (
https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE264581).