Systematic characterization of Brassica napus UBC13 genes involved in DNA-damage response and K63-linked polyubiquitination
BMC Plant Biology volume 23, Article number: 24 (2023)
Ubc13 is the only known ubiquitin conjugating enzyme (Ubc/E2) dedicated to promoting Lys (K)63-linked polyubiquitination, and this process requires a Ubc/E2 variant (UEV). Unlike conventional K48-linked polyubiquitination that targets proteins for degradation, K63-linked polyubiquitination, which is involved in several cellular processes, does not target proteins for degradation but alter their activities.
In this study we report the identification and functional characterization of 12 Brassica napus UBC13 genes. All the cloned UBC13 gene products were able to physically interact with AtUev1D, an Arabidopsis UEV, to form stable complexes that are capable of catalyzing K63-linked polyubiquitination in vitro. Furthermore, BnUBC13 genes functionally complemented the yeast ubc13 null mutant defects in spontaneous mutagenesis and DNA-damage responses, suggesting that BnUBC13s can replace yeast UBC13 in mediating K63-linked polyubiquitination and error-free DNA-damage tolerance.
Collectively, this study provides convincing data to support notions that B. napus Ubc13s promote K63-linked polyubiquitination and are probably required for abiotic stress response. Since plant Ubc13-UEV are also implicated in other developmental and stress responses, this systematic study sets a milestone in exploring roles of K63-linked polyubiquitination in this agriculturally important crop.
Ubiquitin (Ub) is an abundant and one of most highly conserved proteins in eukaryotes, from unicellular yeast to human. Ub is able to conjugate with other proteins by forming an isopeptide bond using its C-terminal glycine (G76) residue and eventually produces a monoubiquitinated substrate . Ub has two key features. One is that its G76 can form an isopeptide bond with a substrate lysine (K) residue. Another is that seven lysine residues in Ub, namely K6, K11, K27, K29, K33, K48 and K63 can be potentially used to form distinct types of poly-Ub chains; at least five of which have been observed in vitro or in vivo .
Ubiquitination, a process of Ub attachment to targeted proteins, regulates diverse cellular activities such as proteasomal and lysosomal degradation, subcellular localization , DNA damage response [4, 5], ribosomal biogenesis , cell cycle progression , apoptosis , mitochondrial inheritance  and transcriptional regulation . Ub conjugated to the target protein can alter the protein stability, localization or activity . The ubiquitination requires three basic enzymatic activities to work in concert to transfer Ub to client substrates and form poly-Ub chains. Firstly, high energy thioester bond is formed with the C terminus of Ub by a Ub activating enzyme (Uba or E1). The Ub molecule is then transferred to the active-site Cys of a Ub-conjugating enzyme (Ubc or E2). The subsequent transfer of Ub to the ε-amino group of a Lys side chain in the substrate is catalyzed by the Ub-charged E2 after binding to a Ub ligase (E3). It is well established that distinct structural and functional information is conveyed by poly-Ub chains bearing different linkages, in which K48-linked Ub chains mainly target protein for degradation by the 26 S proteosome. In contrast, K63-linked chains confer non-proteolytic functions involved in well-characterized pathways including DNA damage response (DDR) [12, 13], NF-κB activation [14, 15], mitochondrial inheritance , plasma membrane protein endocytosis , ribosome function  and cell-cycle checkpoints , in which target proteins conjugated with K63-linked poly-Ub chains affect protein-protein interaction and/or protein localization . Thus far, Ubc13 is the only known E2 specialized in catalyzing K63-linked polyubiquitination, and this activity absolutely requires a Ubc/E2 variant (UEV) as a cofactor [12, 19]. It is well established that K48- and K63-linked poly-Ub chains represent two distinct conformations, in which the 26 S proteosome can only recognize the K48-linked zig-zag conformation, whereas proteins containing specialized tandem Ub-binding domains can recognize K63-linked stretched Ub chains [1, 20].
The model plant Arabidopsis thaliana contains two highly conserved UBC13 genes capable of functionally complementing the yeast ubc13 null mutant in terms of spontaneous mutagenesis and sensitivity to DNA-damaging agents . Both AtUbc13s can interact with yeast and human UEVs  as well as Arabidopsis UEVs . Subsequent studies reveal that AtUBC13s function in apical dominance , DDR , iron metabolism , auxin signaling , low temperature stress response  and plant immunity [26, 27]. In addition, systematic studies in Arabidopsis revealed additional target proteins modified by K63-linked polyubiquitination [28, 29], indicating that plant Ubc13s play multiple roles in plant development and stress responses.
Arabidopsis belongs to the Brassicaceae family, which includes large number of economic crops like rapeseed (Brassica napus) and its value-added breeding derivative canola . Given the involvement of Ubc13 and K63-linked polyubiquitination in agriculturally related traits, we wish to investigate B. napus UBC13 genes. Here we report molecular cloning and functional characterization of UBC13 genes from B. napus, following strategies as outlined in Supplementary Fig. S1. All analyzed BnUbc13s can interact with AtUev1D and catalyze K63-linked poly-Ub chain assembly. Since UBC13 genes have been shown to play critical roles in plant development and responses to both abiotic and biotic stresses including agriculturally important features, this study sets a corner stone in improving the oilseed crop through genetic manipulations.
Identification and bioinformatics analyses of Brassica napus UBC13 genes
To identify B. napus UBC13 genes, the Arabidopsis Ubc13 protein sequence was used to Blast the B. napus protein database (https://plants.ensembl.org/Brassica_napus/Tools/Blast). Eight predicted highly similar proteins (E-values at or below 10− 56), namely BnaA06g11360D (GenBank accession OP380669), BnaC08g38130D (OP380670), BnaAnng13030D (OP380671), BnaC08g17090D (OP380672), BnaA07g34450D (OP380673), BnaA08g23450D (OP380674), BnaC05g12900D (OP380675) and BnaC06g39290D (OP380676), were found and named as BnUbc13A-H, respectively. In addition, four incomplete proteins in the database, BnaC02g25260D (OP380677), BnaA07g38410D (OP380679), BnaC06g20310D (OP380678) and BnaA02g19070D (OP380680), have about 99% sequence identity with AtUbc13. After cDNA amplification and sequencing, these four encoded proteins were also found to be highly homologous to AtUbc13 and subsequently named BnUbc13I-L. Features of the 12 identified BnUBC13 genes and gene products, including number of exons, gene location, protein domain and length, molecular weight (MW), isoelectric point (pI) and grand average of hydropathicity index (GRAVY), are summarized in Supplementary Table S1. All 12 BnUbc13 proteins contain 153 amino acids and their alignment with AtUbc13s is shown in Fig. 1A. AtUbc13A and AtUbc13B differ by only two conserved amino acids . Within BnUbc13s, BnUbc13A, B, I, J and L are identical to AtUbc13B; BnUbc13C and D are identical to each other but differ from AtUbc13B by one amino acid; BnUbc13E and K are identical to each other but differ from AtUbc13B by two amino acids; BnUbc13G, F and H differs from AtUbc13B by 1, 2 and 3 amino acids, respectively. These amino acid substitutions do not affect known functional residues, including M66 (blue asterisk) involved in physical interaction with E3s , E57, F59 and R72  required for the interaction with UEV, and the active site C89 to conjugate Ub . All BnUbc13s appear to be closer to AtUbc13B in sequence than to AtUbc13A, as the two amino acid variations in AtUbc13A were not found in any BnUbc13s (Fig. 1A).
The phylogenetic analysis of BnUbc13s to AtUbc13s was performed. As shown in Fig. 1B, all 12 BnUbc13s can be grouped into two clades separately derived from AtUbc13A and AtUbc13B, and the bootstraps in the clade of AtUbc13B are higher than that of AtUbc13A. To further assess the conservation of BnUBC13 genes in Brassicaceae, we performed synteny analysis among four related species: A. thaliana, B. napus, B. rapa and B. oleracea. As anticipated, the B. napus genome was split into two subgenomes A and C, corresponding to B. rapa and B. oleracea genomes, respectively (Supplementary Fig. S2).
Due to amino acid sequence redundancy, eight BnUBC13 genes (A to H) representing all BnUbc13 variants were used for the subsequent studies.
Functional complementation of yeast ubc13 null mutants by BnUBC13s
Budding yeast UBC13 functions in the error-free DNA-damage tolerance (DDT) pathway [33,34,35]. To ask whether BnUBC13s have same functions as yeast UBC13, pGBT-BnUBC13 plasmids were transformed into the yeast ubc13∆ mutant and cell survival in the presence of various DNA-damaging agents were examined by a serial dilution assay. As previously reported , the yeast ubc13 mutant displayed an increased sensitivity to a variety of DNA-damaging agents including methyl methanesulfonate (MMS), 4-nitroquinoline 1-oxide (4NQO) and ultraviolet (UV) irradiation, while expression of BnUBC13s could functionally complement the yeast ubc13 DDT defect (Fig. 2A). Similarly, in a gradient plate assay, deletion of UBC13 resulted in an increased sensitivity to MMS, while pGBT-BnUBC13s, but not the pGBT9 empty vector, were able to restore cellular resistance to MMS (Fig. 2B).
Yeast UBC13 functions in an error-free DDT pathway in parallel to the translesion DNA synthesis (TLS) pathway in response to lesions that block replication . While mutations in each pathway cause moderate sensitivity to DNA damaging agents, double mutations defective in both pathways result in synergistic effects . Figure 3 shows that in the presence of 0.001% MMS, neither ubc13∆ nor rev3∆ single mutants displayed increased sensitivity; however, the ubc13∆ rev3∆ double mutant did not grow at all in the gradient plat assay. Under the same experimental conditions, expression of any of the eight GAL4BD-BnUBC13 genes, but not GAL4BD alone, was able to rescue the double mutant to the level indistinguishable from the wild-type or single mutants (Fig. 3). The above observations collectively demonstrate that B. napus UBC13 genes can functionally complement the yeast ubc13∆ mutant from killing by DNA-damaging agents.
BnUBC13 genes protect ubc13 cells from spontaneous mutagenesis
Yeast UBC13 is a member of the error-free DDT pathway and plays an important role in protecting yeast cells from spontaneous mutagenesis [33, 34]. Therefore, a spontaneous mutagenesis assay was performed to determine whether BnUBC13s could functionally complement the error-free DDT defect in the yeast. It was apparent from Fig. 4 that inactivation of UBC13 in wild-type yeast cells caused nearly 27-fold increase in the spontaneous mutation rate. The large increase in spontaneous mutagenesis supports a notion that UBC13 plays a vital role in maintaining host genome stability. In contrast, when ubc13Δs cells were transformed with BnUBC13s, the spontaneous mutation rate dropped to near wild-type levels. These results indicate that BnUbc13s can replace Ubc13 in yeast cells to limit endogenous DNA-damage stress.
Physical interactions of BnUbc13s with AtUev1D
One condition for the yeast Ubc13 function is that it has to form a stable complex with Mms2, the only known UEV in yeast cells . The fact that BnUbc13s can replace yeast Ubc13 in spontaneous and DNA-damage induced functions indicates that it must interact with Mms2 in yeast cells. To ask whether BnUbc13s are also able to interact with plant UEVs, a yeast two hybrid (Y2H) assay  was employed to analyze the protein-protein interaction between BnUbc13s and AtUev1D. All eight Gal4BD-BnUbc13s gave positive results with Gal4AD-AtUev1D under high stringency (SD-His + 3-AT and SD-Ade) conditions in comparison to negative controls including Gal4BD-BnUbc13s with Gal4AD and Gal4BD with Gal4AD-AtUev1D (Fig. 5). All the above interactions are deemed robust and strong, and no difference in the interaction strength among BnUbc13s was observed. Indeed, these interactions seem to be specific between BnUbc13s and AtUev1D as neither of the proteins alone was able to activate reporter genes in the Y2H assay. Therefore, results from the Y2H assay indicate that all BnUbc13s are able to physically interact with AtUev1D.
The physical interaction between BnUbc13 and AtUev1D was further confirmed independently in vitro by a GST-affinity pulldown assay. In this experiment, bacterial cells were transformed with both His6-tagged BnUbc13s and GST-tagged AtUev1D, and the produced proteins in bacterial cells were co-purified by adding to a column containing glutathione beads. After incubation, washing and elution, all eight His6-BnUbc13s were found to be co-eluted with GST-AtUev1D, but not with GST (Fig. 6). Hence, all BnUbc13s can form stable heterodimers with AtUev1D in vitro.
Dual rescue of yeast ubc13 mms2 by BnUBC13s and AtUEV1D
To assess in vivo complex formation and functions between BnUbc13s and AtUev1D, the ubc13 mms2 double mutant was created and co-transformed with pGBT-BnUBC13s and pGAD-AtUEV1D, or with their respective empty vectors. When the double mutant cells were transformed with only pGBT-BnUBC13s or pGAD-AtUEV1D, the transformed cells did not display enhanced resistance to MMS (Fig. 7), implying that both Ubc13 and UEV are required for the DDT function. In contrast, double mutant cells carrying both pGBT-BnUBC13 and pGAD-AtUEV1D plasmids displayed MMS resistance comparable to the wild-type level (Fig. 7). Since BnUBC13s and AtUEV1D can jointly complement ubc13Δ and mms2Δ defects in yeast, one can envisage that BnUbc13s must be able to bind AtUev1D in yeast cells to form a functional E2 complex and promote K63-linked polyubiquitination on PCNA , which is a highly conserved process within eukaryotes including plants .
BnUbc13 mediates K63-linked polyubiquitination with AtUev1 in vitro
To directly examine whether BnUbc13s can promote K63-linked poly-Ub chain assembly, we performed an in vitro Ub conjugation assay with selected BnUbc13s. To date, Ubc13 is the only known E2 capable of catalyzing K63-linked polyubiquitination; however, Ubc13 alone is insufficient and requires a UEV as a cofactor . As shown in Fig. 8, in our reconstituted Ub conjugation assay, neither BnUbc13s alone (Fig. 8A, lanes 2, 6 and 10; Fig. 8B, lanes 3 and 7) nor AtUev1D alone (Fig. 8A, lanes 3, 7 and 11; and Fig. 8B, lanes 3 and 7) was able to catalyze the poly-Ub chain formation. When both BnUbc13s and AtUev1D are present (Fig. 8A, lanes 1, 5 and 9; Fig. 8B, lanes 1 and 5), free poly-Ub chains are formed. These poly-Ub chains are deemed to be K63-linked, since when Ub was replaced by Ub-K63R, the poly-Ub chain formation was completely abolished (Fig. 8A, lanes 4, 8 and 12; Fig. 8B, lanes 4 and 8). These results demonstrate that BnUbc13s and AtUev1D can jointly form K63-linked poly-Ub chains in vitro.
In this study, we identified 12 UBC13 genes from the B. napus genome. Most non-plant eukaryotic organisms contain only one UBC13 gene . Characterized diploid plant genomes contain either one  or two [21, 41, 42] UBC13 genes. Hence, it is surprising that the amphidiploid B. napus contains 12 UBC13 genes. Although their DNA sequences vary, the encoded BnUbc13s are highly conserved with no more than 3 amino acid differences. Furthermore, five BnUbc13s are identical in sequence with AtUbc13B, while the two amino acid variations found in AtUbc13A were not shared by any BnUbc13s, supporting a previous speculation through phylogenetic analysis that AtUbc13A was derived from gene duplication after Arabidopsis speciation . Similarly, the remaining seven BnUBC13 genes encoding BnUbc13 variants were likely evolved within the Brassica genus. It remains interesting what evolutionary pressure caused sequential UBC13 gene duplications in B. rapa and B. oleracea, and then in B. napus.
In this study, we cloned BnUBC13 genes encoding all BnUbc13 variants and characterized their functions in budding yeast. The only known function of yeast UBC13 to date is its involvement in error-free DDT, which is achieved by promoting K63-linked polyubiquitination of PCNA at the K164 residue  that promotes template switch to bypass replication-blocking lesions . This study demonstrated that BnUBC13s can functionally complement the yeast ubc13 null defect in error-free DDT. Firstly, heterologous expression of BnUBC13s rescued ubc13∆ and ubc13∆ rev3∆ cells from killing by DNA-damaging agents, indicating that they function in DDR. Secondly, BnUBC13s could also limit the increased spontaneous mutagenesis in ubc13∆ cells, suggesting that they protect cells from genome instability under non-stressed conditions. This study also demonstrated physical interactions between BnUbc13s and AtUev1 by Y2H and GST pulldown assays, as well as by a dual rescuing experiment in which BnUbc13s must interact with AtUev1 in yeast cells to restore the DDT function, hence confirming that they form stable heterodimers in vitro and in vivo. Finally, selected BnUbc13 proteins were used in an in vitro Ub conjugation assay, which demonstrated that only in the presence of both BnUbc13 and AtUev1 could free poly-Ub chains be formed, and that these chains were K63 linked. In addition, the dual rescue experiment of yeast ubc13∆ mms2∆ mutant cells by BnUBC13s and AtUEV1D also implies that the BnUbc13-AtUev1D heterodimer must promote PCNA polyubiquitination to confer error-free DDT functions. Although this study did not use a UEV from B. napus, the facts that, like Ubc13s, eukaryotic UEVs are also highly conserved and that BnUbc13s could physically interact with both AtUev1 and yeast Mms2 to promote K63-linked polyubiquitination support a notion that B. napus Ubc13s and UEVs also form stable heterodimers to promote K63-linked polyubiquitination and are required for plant DDR.
All assays employed in this study did not reveal qualitative or quantitative difference among all BnUbc13s, which is not surprising, since they are highly conserved in sequence with very few amino acid variations and that all amino acid residues with functional implications are identical. Indeed, plant UBC13 appears to be a housekeeping gene, particularly when the plant genome contains only one UBC13 gene . However, the similarity in biochemical activity does not necessarily mean that all BnUBC13 genes function equally and redundantly, as one cannot rule out possibilities that these genes contain different promoter sequences and display different spatial and temporal expression patterns, particularly in response to environmental stresses.
Plant Ubc13-Uev1 complexes are not only required for DDR, but also for other agriculturally related processes. It appears that in multicellular organisms, UEVs serve as regulatory subunits to alter substrate specificity. In human cells, two UEVs, Uev1A and Mms2 contain different N-terminal extension and function in NF-kB activation and DDR, respectively . Plants also contain two classes of Uev1s with different C-terminal extension and cellular activities . We are particularly interested in roles of plant Ubc13 and its mediated K63-linked polyubiquitination in auxin signaling [23, 25, 45], nutrient metabolism  and innate immunity against pathogen infections [26, 27, 46], as alteration of these processes may improve quality and yield of the oilseed crop.
All methods were performed in accordance with the relevant guidelines and regulations.
Plant materials and yeast cell culture
Brassica napus L. cv Westar was developed at the Saskatoon Research and Development Centre, Agriculture and Agri-Food Canada (AAFC), and registered by AAFC in Canada in 1982 under the registration # 2238. The doubled haploid seed of Westar used in the study was produced at Nutrien Ag Solutions Inc. (Saskatoon, SK, Canada), and provided to AAFC for research uses, including this study. The seeds were planted in Sunshine #3 soil (SunGro Horticulture, Vancouver, BC) in 4-inch pots, and placed in a growth chamber at 22/16°C (day/night) with a 16 h photoperiod until sample collection.
Saccharomyces cerevisiae haploid strains used in this study are listed in Supplementary Table S2. Either rich YPD or a synthetic dextrose (SD) medium (0.67% Bacto-yeast nitrogen base without amino acids, 2% glucose) were used to grow yeast cells at 30 oC with necessary nutrients as recommended . When making solid plates, 2% agar was added to either YPD or SD medium prior to autoclaving. Yeast transformation was carried out by dimethyl sulfoxide (DMSO)-enhanced method as described .
The genome structure and location of BnUBC13 genes were retrieved from Ensemble Plants (https://plants.ensembl.org/index.html), illustrating the annotated protein-coding regions. The protein domain was found in the NCBI Conserved Domain Search (https://www.ncbi.nlm.nih.gov/Structure/cdd/wrpsb.cgi), while the subcellular location was predicted by UniProt (https://www.uniprot.org/). The PI and MW were calculated in Expasy (https://web.expasy.org/compute_pi/) by entering the sequences of each gene. The SMS Protein GRAVY software (https://www.bioinformatics.org/sms2/protein_gravy.html) was used to calculate the GRAVY value from the entered FASTA sequences.
A phylogenetic tree based on B. napus and Arabidopsis Ubc13 family amino acid sequences was constructed by using a maximum likelihood (ML) estimation method of a JTT matrix-based model of MEGA7.0.26. The bootstrap consensus tree inferred from 1,000 replicates. Values < 50% are omitted (parameters: Gaps/Missing data treatment -Partial deletion; Site coverage cutoff − 95%). The phylogenetic tree was exported as Newick format and opened in the Interactive Tree Of Life (iTOL), a web-based tool for phylogenetic tree manipulation and annotation .
To identify homologs of UBC13 family genes and the conservation of these homologous genes in Brassicaceae, the synteny analysis was performed among four species: A. thaliana, B. napus, B. rapa and B. oleracea. The genome assemblies for these four species were obtained from EnsemblPlants (https://plants.ensembl.org/info/data/ftp, version 55). The synteny analysis was performed by MCScanX  with default parameters from top-five BlastP hits. Two Arabidopsis UBC13 genes (AT1G16890 and AT1G78870) were used to search for gene pairs in three Brassica species. When performing the synteny analysis among Brassica species, the genome of B. napus was split into two subgenomes (A and C), with gene pairs detected between the A subgenome and B. rapa, and between the C subgenome and B. oleracea, correspondingly. Circos plot for synteny analysis was further performed through http://circos.ca/.
Cloning B. napus UBC13 cDNAs and plasmid construction
Total RNA was extracted from B. napus seedlings with a TRIzol reagent and used for reverse transcriptase (RT)-PCR with the SuperScript RT-PCR III system (Invitrogen) according to manufacturer’s instructions. Full-length BnUBC13 coding sequences were cloned from the above resulting cDNA library. All primers used for the PCR are shown in Table S3 and include a SalI restriction site in forward primers and a BamH1 restriction site in reverse primers. The PCR products were cleaved by SalI and BamHI, and subsequently cloned into Y2H vectors pGBT9E and pGAD424E, which were derived from pGBT9 (Gal4BD) and pGAD424 (Gal4AD) , respectively, with a frameshift at the multiple cloning site. The identity of the cloned inserts was confirmed by DNA sequencing.
Yeast survival assays
A yeast strain HK578-10D and its isogenic mutant derivatives (Table S2) were transformed or co-transformed with pGBT9 and/or pGAD424 based plasmids as indicated, and at least five colonies were selected after 3-day incubation on the selective plate (SD-Leu plates for the single transformation and SD-Leu-Trp for double transformations) and streaked onto the same selective plates. The cell survival assays were as previously described . Briefly, for a gradient plate assay, transformed yeast cells were used to inoculate liquid SD minimal medium. Equal number of cells were taken after overnight incubation and imprinted onto YPD alone or YPD gradient plates containing 0.025% MMS gradient in the bottom layer. 0.1 mL of overnight culture was mixed with 0.4 mL of sterile distilled water and 0.5 mL of molten agar on a sterile microscopic slide and the mixture was printed on the plates across the gradient with another microscopic slide. All the plates were incubated at 30 oC for 3 days before photography.
For a serial dilution assay, overnight cultures were adjusted with equal cell density, made tenfold serial dilutions with sterile distilled water, and 5 µL of each sample was spotted to the plates containing desired concentrations of MMS or 4NQO. For the UV irradiation, samples were spotted to the YPD plates and exposed to predetermined UV irradiation doses in a 254 nM UV crosslinker (UVP). After the liquid was absorbed, the plates were incubated at 30 oC for 3 days before photography.
Spontaneous mutagenesis assay
A yeast strain DBY747 and its isogenic ubc13∆ derivative WXY849 (Table S2) were used in the spontaneous mutagenesis assay. These strains bear a trp1-289 amber mutation that can be reverted to Trp+ by different mutation events . WXY849 was transformed with pGAD-BnUbc13s and the transformants were selected on SD-Leu plates. The overnight yeast cell culture was used to inoculate 5 mL of YPD liquid medium with a final concentration of 20 cells/mL, which was incubated at 30 oC for 3 days. Cells were collected by centrifugation at 4,000 rpm, resuspended in sterile distilled water and plated on SD-Leu to count for the total colony forming units and on SD-Leu-Trp to count for Trp+ revertant. Spontaneous mutation rates were calculated as described .
Yeast two-hybrid analysis
The Y2H strain PJ69-4A (Table S2) was co-transformed with Gal4BD and Gal4AD constructs. SD-Leu-Trp plates were used to select the co-transformed colonies. At least five independent colonies from each co-transformation were printed onto SD-Leu-Trp plates and SD-Leu-Trp-His selective plates with different concentrations of 1,2,4 amino triazole (3 A-T) to test activation of the PGAL1-HIS3 reporter gene and onto Leu-Trp-Ade plates to test activation of the PGAL2-ADE2 reporter gene. The above plates were incubated at 30 oC for the indicated time before photography.
Recombinant protein production and extraction
The BnUBC13 open-reading frames (ORFs) were isolated from pGBT-BnUbc13s and cloned into pET30a (Novogene) as a His6 fusion. The resulting pET-BnUbc13s were co-transformed with pGEX-AtUev1D  into Epicurian coli BL21-CodonPLus (DE3)-RIL strain (Thermo Sci. 2,287,225). The transformed cells were cultured overnight in LB + Amp + Kan, diluted 1:50 into fresh culture and incubated until OD600 nm of 0.6–0.8. The His6-BnUbc13s and GST-AtUev1D fusion proteins were induced by adding isopropyl-b-D-thiogalactopyranoside (IPTG) to the final concentration of 0.2 mM and the incubation continued for 6 h. The cells were harvested by 8,000 rpm centrifugation in a Beckman Coulter Avanti JA17 rotor for 1 h at 4 °C, resuspended in phosphate-buffered saline (PBS, 140 mM NaCl, 2.7 mM KCl, 10 mM Na2HPO4, 1.8 mM KH2PO4, pH 7.3), and passed through Constant Systems one shot cell disrupter at 25 PSI. The resulting crude extract was centrifuged at 17,000 rpm in the same JA17 rotor for 30 min at 4 oC, and the soluble fraction was used for the GST pulldown assay.
GST pulldown assay
The GST-pulldown assay was carried out using Glutathione Sepharose 4B Microspin™ GST purification columns as previously described . Purified GST-AtUev1D and His6-BnUbc13s were added to individual columns with glutathione beads and incubated for 1 h at 4 oC. Finally, wash buffer (1x PBS buffer with 350 mM NaCl) was used to wash the beads by running through the column. The eluted samples were subjected to 15% SDS-PAGE and visualized by Coomassie blue staining.
Ub conjugation reaction
E. coli BL21-RIL cells separately transformed with pET-BnUbc13s and pGEX-AtUev1D were incubated and induced by IPTG as described above, and the cell extracts were used for protein purification with Bio-Rad poly-prep chromatography columns (731–1550) containing Ni-NTA and glutathione beads, respectively. After washing beads with respective lysis buffers, affinity-purified recombinant proteins were eluted from the columns by using 20 mL of His6 elution buffer for His6-BnUbc13s and GST elution buffer for GST-AtUev1D. The ubiquitination assay kit containing Ub thioester/conjugation initiation reagents was purchased from Abcam (ab139467). A 20-mL reaction mixture contained E1, Ub, MgATP, reaction buffer from the kit, plus His6-BnUbc13 and GST-AtUev1D prepared from this study. The Ub-K63R protein was purchased from Abcam (UM-K63R). Conjugation reactions were performed at 37 °C for 4 h according to manufacturer’s instructions, followed by running of 15% SDS-PAGE and western blotting using polyclonal goat anti-Ub antibodies (Bio-Rad).
Availability of data and materials
The BnUBC13 sequences are deposited to the GenBank with accession numbers OP380669-OP380680. The datasets acquired and/or analyzed during the current study are available from the corresponding author on reasonable request.
Komander D, Rape M. The ubiquitin code. Annu Rev Biochem. 2012;81:203–29.
Volk S, Wang M, Pickart CM. Chemical and genetic strategies for manipulating polyubiquitin chain structure. Methods Enzymol. 2005;399:3–20.
Kerscher O, Felberbaum R, Hochstrasser M. Modification of proteins by ubiquitin and ubiquitin-like proteins. Annu Rev Cell Dev Biol. 2006;22:159–80.
Jentsch S, McGrath JP, Varshavsky A. The yeast DNA repair gene RAD6 encodes a ubiquitin-conjugating enzyme. Nature. 1987;329(6135):131–4.
Pastushok L, Xiao W. DNA postreplication repair modulated by ubiquitination and sumoylation. Adv Protein Chem. 2004;69:279–306.
Finley D, Bartel B, Varshavsky A. The tails of ubiquitin precursors are ribosomal proteins whose fusion to ubiquitin facilitates ribosome biogenesis. Nature. 1989;338(6214):394–401.
Harper JW. A phosphorylation-driven ubiquitination switch for cell-cycle control. Trends Cell Biol. 2002;12(3):104–7.
Zhang HG, Wang J, Yang X, Hsu HC, Mountz JD. Regulation of apoptosis proteins in cancer cells by ubiquitin. Oncogene. 2004;23(11):2009–15.
Fisk HA, Yaffe MP. A role for ubiquitination in mitochondrial inheritance in Saccharomyces cerevisiae. J Cell Biol. 1999;145(6):1199–208.
Kao KC, Yang YL, Boscolo R, Sabatti C, Roychowdhury V, Liao JC. Transcriptome-based determination of multiple transcription regulator activities in Escherichia coli by using network component analysis. Proc Natl Acad Sci USA. 2004;101(2):641–6.
Dorval V, Fraser PE. SUMO on the road to neurodegeneration. Biochim Biophys Acta. 2007;1773(6):694–706.
Hofmann RM, Pickart CM. Noncanonical MMS2-encoded ubiquitin-conjugating enzyme functions in assembly of novel polyubiquitin chains for DNA repair. Cell. 1999;96(5):645–53.
Hoege C, Pfander B, Moldovan GL, Pyrowolakis G, Jentsch S. RAD6 -dependent DNA repair is linked to modification of PCNA by ubiquitin and SUMO. Nature. 2002;419(6903):135–41.
Deng L, Wang C, Spencer E, Yang L, Braun A, You J, Slaughter C, Pickart C, Chen ZJ. Activation of the IkB kinase complex by TRAF6 requires a dimeric ubiquitin-conjugating enzyme complex and a unique polyubiquitin chain. Cell. 2000;103(2):351–61.
Zhou H, Wertz I, O’Rourke K, Ultsch M, Seshagiri S, Eby M, Xiao W, Dixit VM. Bcl10 activates the NF-kB pathway through ubiquitination of NEMO. Nature. 2004;427(6970):167–71.
Galan JM, Haguenauer-Tsapis R. Ubiquitin lys63 is involved in ubiquitination of a yeast plasma membrane protein. EMBO J. 1997;16(19):5847–54.
Spence J, Gali RR, Dittmar G, Sherman F, Karin M, Finley D. Cell cycle-regulated modification of the ribosome by a variant multiubiquitin chain. Cell. 2000;102(1):67–76.
Bothos J, Summers MK, Venere M, Scolnick DM, Halazonetis TD. The Chfr mitotic checkpoint protein functions with Ubc13-Mms2 to form Lys63-linked polyubiquitin chains. Oncogene. 2003;22(46):7101–7.
McKenna S, Spyracopoulos L, Moraes T, Pastushok L, Ptak C, Xiao W, Ellison MJ. Noncovalent interaction between ubiquitin and the human DNA repair protein Mms2 is required for Ubc13-mediated polyubiquitination. J Biol Chem. 2001;276(43):40120–6.
Sims JJ, Cohen RE. Linkage-specific avidity defines the lysine 63-linked polyubiquitin-binding preference of Rap80. Mol Cell. 2009;33(6):775–83.
Wen R, Newton L, Li G, Wang H, Xiao W. Arabidopsis thaliana UBC13: implication of error-free DNA damage tolerance and Lys63-linked polyubiquitylation in plants. Plant Mol Biol. 2006;61(1–2):241–53.
Wen R, Torres-Acosta JA, Pastushok L, Lai X, Pelzer L, Wang H, Xiao W. Arabidopsis UEV1D promotes lysine-63-linked polyubiquitination and is involved in DNA damage response. Plant Cell. 2008;20(1):213–27.
Yin XJ, Volk S, Ljung K, Mehlmer N, Dolezal K, Ditengou F, Hanano S, Davis SJ, Schmelzer E, Sandberg G, et al. Ubiquitin lysine 63 chain forming ligases regulate apical dominance in Arabidopsis. Plant Cell. 2007;19(6):1898–911.
Li W, Schmidt W. A lysine-63-linked ubiquitin chain-forming conjugase, UBC13, promotes the developmental responses to iron deficiency in Arabidopsis roots. Plant J. 2010;62(2):330–43.
Wen R, Wang S, Xiang D, Venglat P, Shi X, Zang Y, Datla R, Xiao W, Wang H. UBC13, an E2 enzyme for Lys63-linked ubiquitination, functions in root development by affecting auxin signaling and Aux/IAA protein stability. Plant J. 2014;80(3):424–36.
Wang L, Wen R, Wang J, Xiang D, Wang Q, Zang Y, Wang Z, Huang S, Li X, Datla R, et al. Arabidopsis UBC13 differentially regulates two programmed cell death pathways in responses to pathogen and low-temperature stress. New Phytol. 2019;221(2):919–34.
Yao D, Arguez MA, He P, Bent AF, Song J. Coordinated regulation of plant immunity by poly(ADP-ribosyl)ation and K63-linked ubiquitination. Mol Plant. 2021;14:208–103.
Yang K, Xiao W. Functions and mechanisms of the Ubc13-UEV complex and lysine 63-linked polyubiquitination in plants. J Exp Bot. 2022;73(16):5372–87.
Romero-Barrios N, Monachello D, Dolde U, Wong A, San Clemente H, Cayrel A, Johnson A, Lurin C, Vert G. Advanced cataloging of lysine-63 polyubiquitin networks by genomic, interactome, and sensor-based proteomic analyses. Plant Cell. 2020;32(1):123–38.
Franzke A, German D, Al-Shehbaz IA, Mummenhoff K. Arabidopsis family ties: molecular phylogeny and age estimates in Brassicaceae. Taxon. 2009;58(2):425–37.
Wooff J, Pastushok L, Hanna M, Fu Y, Xiao W. The TRAF6 RING finger domain mediates physical interaction with Ubc13. FEBS Lett. 2004;566(1–3):229–33.
Pastushok L, Moraes TF, Ellison MJ, Xiao W. A single Mms2 “key” residue insertion into a Ubc13 pocket determines the interface specificity of a human Lys63 ubiquitin conjugation complex. J Biol Chem. 2005;280(18):17891–900.
Broomfield S, Chow BL, Xiao W. MMS2, encoding a ubiquitin-conjugating-enzyme-like protein, is a member of the yeast error-free postreplication repair pathway. Proc Natl Acad Sci USA. 1998;95(10):5678–83.
Brusky J, Zhu Y, Xiao W. UBC13, a DNA-damage-inducible gene, is a member of the error-free postreplication repair pathway in Saccharomyces cerevisiae. Curr Genet. 2000;37(3):168–74.
Fan L, Bi T, Wang L, Xiao W. DNA-damage tolerance through PCNA ubiquitination and sumoylation. Biochem J. 2020;477(14):2655–77.
Xiao W, Chow BL, Fontanie T, Ma L, Bacchetti S, Hryciw T, Broomfield S. Genetic interactions between error-prone and error-free postreplication repair pathways in Saccharomyces cerevisiae. Mutat Res. 1999;435(1):1–11.
Fields S, Song O. A novel genetic system to detect protein-protein interactions. Nature. 1989;340(6230):245–6.
Yang K, Xiao W. Ubc13-UEV mediated K63-linked polyubiquitination in plants. J Exp Bot. 2022. In press.
Bai Z, Li Z, Xiao W. Drosophila bendless catalyzes K63-linked polyubiquitination and is involved in the response to DNA damage. Mutat Res. 2018;808:39–47.
Zang Y, Wang Q, Xue C, Li M, Wen R, Xiao W. Rice UBC13, a candidate housekeeping gene, is required for K63-linked polyubiquitination and tolerance to DNA damage. Rice. 2012;5(1):24.
Guo H, Wang L, Hu R, He Y, Xiao W. Molecular cloning and functional characterization of Physcomitrella patens UBC13-UEV1 genes required for Lys63-linked polyubiquitination. Plant Sci. 2020;297:110518.
Guo H, Wen R, Liu Z, Datla R, Xiao W. Molecular cloning and functional characterization of two Brachypodium distachyon UBC13 genes whose products promote K63-linked polyubiquitination. Front Plant Sci. 2016;6:1222.
Andersen PL, Zhou H, Pastushok L, Moraes T, McKenna S, Ziola B, Ellison MJ, Dixit VM, Xiao W. Distinct regulation of Ubc13 functions by the two ubiquitin-conjugating enzyme variants Mms2 and Uev1A. J Cell Biol. 2005;170(5):745–55.
Wang Q, Liu M, Zang Y, Xiao W. The C-terminal extension of Arabidopsis Uev1A/B with putative prenylation site plays critical roles in protein interaction, subcellular distribution and membrane association. Plant Sci. 2020;291:110324.
Leitner J, Petrášek J, Tomanov K, Retzer K, Pařezová M, Korbei B, Bachmair A, Zažímalová E, Luschnig C. Lysine63-linked ubiquitylation of PIN2 auxin carrier protein governs hormonally controlled adaptation of Arabidopsis root growth. Proc Natl Acad Sci USA. 2012;109(21):8322–7.
Mural RV, Liu Y, Rosebrock TR, Brady JJ, Hamera S, Connor RA, Martin GB, Zeng L. The tomato Fni3 lysine-63-specific ubiquitin-conjugating enzyme and suv ubiquitin E2 variant positively regulate plant immunity. Plant Cell. 2013;25(9):3615–31.
Xiao W. Yeast protocols, Second Edition. Vol. 313. Totowa: Humana Press; 2006.
Hill J, Donald KA, Griffiths DE. DMSO-enhanced whole cell yeast transformation. Nucleic Acids Res. 1991;19(20):5791.
Letunic I, Bork P. Interactive tree of life (iTOL) v5: an online tool for phylogenetic tree display and annotation. Nucleic Acids Res. 2021;49(W1):W293–6.
Wang Y, Tang H, Debarry JD, Tan X, Li J, Wang X, Lee TH, Jin H, Marler B, Guo H, et al. MCScanX: a toolkit for detection and evolutionary analysis of gene synteny and collinearity. Nucleic Acids Res. 2012;40(7):e49.
Bartel PL, Fields S. Analyzing protein-protein interactions using two-hybrid system. Methods Enzymol. 1995;254:241–63.
Xu X, Lambrecht AD, Xiao W. Yeast survival and growth assays. In: Xiao W, editor. Yeast Protocols. Vol. 1163. 3rd ed. New York: Humana Press; 2014. p. 183–91.
Xiao W, Samson L. In vivo evidence for endogenous DNA alkylation damage as a source of spontaneous mutation in eukaryotic cells. Proc Natl Acad Sci USA. 1993;90(6):2117–21.
Williamson MS, Game JC, Fogel S. Meiotic gene conversion mutants in Saccharomyces cerevisiae. I. Isolation and characterization of pms1-1 and pms1-2. Genetics. 1985;110(4):609–46.
The authors wish to thank the Xiao Laboratory members for technical assistance and helpful discussion.
This work was supported by the Saskatchewan Agriculture Development and Saskatchewan Canola Development Commission Fund No. 20180243 to WX.
Ethics approval and consent to participate
All methods were performed in accordance with the relevant guidelines and regulations.
Consent for publication
The authors declare that they have no known competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
The characteristics of identified Brassica napus UBC13 genes and gene products. Table S2. Saccharomyces cerevisiae strains. Table S3. Primers used to amplify BnUBC13 genes. Fig. S1. A flowchart outlining experimental designs of this study. Fig. S2. Synteny analysis of UBC13 family genes among Arabidopsis thaliana, Brassica napus, B. rapa and B. oleracea. Identified gene pairs were displayed by circos plots (http://circos.ca/). All chromosomes from four species are indicated by blocks scaled by chromosome length with chromosome names. The UBC13 family genes and chromosomes from different species are indicated by different colors: red, A. thaliana; yellow, B. napus; purple, B. rapa and green, B. oleracea. Synteny gene pairs are linked by lines: pink, genes with a close relationship to AtUBC13A; brown, genes with a close relationship to AtUBC13B. For gene ID from B. rapa, the species name “_BraROA” is omitted to reduce the length for plotting. Fig. S3. Original images of (A) Fig. 6A and (B) Fig. 6B. Fig. S4. Original images of (A) Fig. 8A and (B) Fig. 8B.
About this article
Cite this article
Kumasaruge, I., Wen, R., Wang, L. et al. Systematic characterization of Brassica napus UBC13 genes involved in DNA-damage response and K63-linked polyubiquitination. BMC Plant Biol 23, 24 (2023). https://doi.org/10.1186/s12870-023-04035-y