Next Article in Journal
Role of Terpenophenolics in Modulating Inflammation and Apoptosis in Cardiovascular Diseases: A Review
Next Article in Special Issue
Genetic Analysis Based on CRISPR/Cas9 Technology in Plants
Previous Article in Journal
The Effect of Ketamine on Endoplasmic Reticulum Stress in Rats with Neuropathic Pain
Previous Article in Special Issue
Transcriptome and Metabolite Profiling of Tomato SGR-Knockout Null Lines Using the CRISPR/Cas9 System
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
IJMS
Volume 24
Issue 6
10.3390/ijms24065338
ijms-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Knockout Mutants of OsPUB7 Generated Using CRISPR/Cas9 Revealed Abiotic Stress Tolerance in Rice

by
Me-Sun Kim
1,†,
Seo-Rin Ko
1,†,
Yu Jin Jung
2,
Kwon-Kyoo Kang
2,
Yung-Jo Lee
3 and
Yong-Gu Cho
1,*
1
Department of Crop Science, College of Agriculture and Life & Environment Sciences, Chungbuk National University, Cheongju 28644, Republic of Korea
2
Division of Horticultural Biotechnology, Hankyong National University, Anseong 17579, Republic of Korea
3
Institute of Korean Prehistory, Cheongju 28763, Republic of Korea
*
Author to whom correspondence should be addressed.
†
These authors contributed equally to this work.
Int. J. Mol. Sci. 2023, 24(6), 5338; https://0-doi-org.brum.beds.ac.uk/10.3390/ijms24065338
Submission received: 10 January 2023 / Revised: 15 February 2023 / Accepted: 6 March 2023 / Published: 10 March 2023
(This article belongs to the Special Issue Genetic Analysis Based on CRISPR/Cas9 Technology)
Browse Figures
Review Reports Versions Notes

Abstract

:
Plants produce and accumulate stress-resistant substances when exposed to abiotic stress, which involves a protein conversion mechanism that breaks down stress-damaged proteins and supplies usable amino acids. Eukaryotic protein turnover is mostly driven by the ubiquitination pathway. Among the three enzymes required for protein degradation, E3 ubiquitin ligase plays a pivotal role in most cells, as it determines the specificity of ubiquitination and selects target proteins for degradation. In this study, to investigate the function of OsPUB7 (Plant U-box gene in Oryza sativa), we constructed a CRISPR/Cas9 vector, generated OsPUB7 gene-edited individuals, and evaluated resistance to abiotic stress using gene-edited lines. A stress-tolerant phenotype was observed as a result of drought and salinity stress treatment in the T2 OsPUB7 gene-edited null lines (PUB7-GE) lacking the T-DNA. In addition, although PUB7-GE did not show any significant change in mRNA expression analysis, it showed lower ion leakage and higher proline content than the wild type (WT). Protein–protein interaction analysis revealed that the expression of the genes (OsPUB23, OsPUB24, OsPUB66, and OsPUB67) known to be involved in stress increased in PUB7-GE and this, by forming a 1-node network with OsPUB66 and OsPUB7, acted as a negative regulator of drought and salinity stress. This result provides evidence that OsPUB7 will be a useful target for both breeding and future research on drought tolerance/abiotic stress in rice.

1. Introduction

Intracellular proteins play pivotal roles in signaling, transport, catalysis, membrane fusion, cellular protection, and regulation of biological processes [1]. However, proteins may be damaged in various ways; for example, their tertiary structure may be damaged, they may undergo non-specific aggregation, or there might be imbalance in their homeostasis due to exposure to various stresses [2]. The cyclical process of degrading damaged proteins to generate usable amino acids and resynthesizing them into necessary proteins is called protein turnover, and up to 80% of the proteins in living organisms are degraded by the ubiquitin-proteasome system [3]. The ubiquitin-proteasome system regulates the stability and activity of many proteins and affects various cellular processes in higher plants, including signal transduction, cell division, and responses to biotic and abiotic stresses [4]. It is also a specialized proteolytic system that plays an essential role in controlling protein degradation [5].
Ubiquitin is a small 8.6 kDa protein composed of 76 amino acids that has a C-terminal tail and seven lysine residues. In addition, human and yeast ubiquitins are highly conserved, with 96% sequence identity [6,7]. Ubiquitin plays the most central role in the ubiquitin-proteasome system; it binds to a target protein and acts as a labeling agent. Ubiquitin forms various types of poly-ubiquitin chains, which are detected and degraded by the proteasome [8,9]. The ubiquitin-proteasome system requires a cascade of three major enzymes: ubiquitin-activating enzyme (E1), ubiquitin-conjugating enzyme (E2), and ubiquitin-ligase (E3). The E1 enzyme activates ubiquitin via ATP-dependent reactions, which results in the formation of a high-energy thiol ester bond that connects the glycine at the C-terminal of ubiquitin and the cysteine in E1. The E2 enzyme receives the activated ubiquitin from E1, forms a thiol ester bond, and catalyzes the formation of an iso-peptide bond between ubiquitin and the target protein that is degraded by the E3 enzyme. Finally, the ubiquitinated protein is recognized and degraded by the 26S proteasome [10,11,12].
Compared to the 1 and 15 types of E1 and E2 enzymes known, respectively, hundreds of isoforms of the E3 enzyme have been identified, which detect target proteins for degradation and provide specificity to the ubiquitin-proteasome system [5,13]. E3 ligases are classified into three types: HECT (Homologous to E6-AP Carboxyl Terminus), RING finger (Really Interesting New Gene finger), and U-box, depending on the method via which ubiquitin is linked to the E2 binding domain or the presence of a zinc (Zn) residue [14]. Unlike the RING and U-box types, the HECT type directly attaches ubiquitin to the HECT domain and delivers it to the substrate proteins. The U-box and RING types allow ubiquitin to be directly transferred from E2 to the substrate protein and possess the same folded structure; however, the U-box type does not contain the zinc residue [15,16,17]. U-box proteins are widely present in eukaryotes, and 2, 64, and 77 U-box proteins have been identified in yeast, Arabidopsis, and rice, respectively [18,19,20]. The presence of more U-box proteins in highly advanced plants is possibly because the U-box E3 ligase plays an important role in plants; hence, it is subdivided into more types [21].
A recent study revealed that plant U-box (PUB) family members are involved in biotic and abiotic stress responses. In Arabidopsis, PUB12 and PUB13 interact with abscisic acid (ABA)-INTENSITIVE 1 (ABI1), a key PP2C protein involved in ABA signaling, to regulate the response to ABA [22]. The pub18 and pub19 double mutants showed resistance to drought stress by enhancing stomatal closure via interaction with Exo70B1, which is known to respond to ABA signaling when exposed to ABA [23]. CMPG1-V, a U-box type E3 ligase, contributes to defense responses against powdery mildew and broad-spectrum disease resistance due to an increase in the expression of salicylic acid-responsive genes and accumulation of H2O2 in plants [24]. The pub13 mutants have been shown to interact with LYK5 and induce plant immune responses via chitooctaose-induced rapid reactions such as reactive oxygen species production and mitogen-activated protein (MAP) kinase phosphorylation [25]. PUB genes in rice are classified according to the domain structure of the gene, and OsPUB genes with U-box/ARM domains are known to be tolerant to abiotic and biotic stress [26]. OsPUB15 transgenic rice caused severe growth delays and fatal phenotypes in seedling, and transcription levels of OsPUB15 knockout mutants were shown to increase with H2O2, salinity, and drought stress [27]. Knockout mutants of OsPUB16 generated by CRISPR/Cas9 showed improved water resistance by mediating ubiquitination and degradation of OsMADS23, a substrate of OSMOTIC STRESS/ABA-ACTIVATED PROTEIN KINASE 9 (SAPK9), and OsPUB16 over-expression lines showed an increased water deficiency susceptibility to stress [28].
The CRISPR/Cas9 technology is widely used in various organisms because of its ease of use and accuracy. Furthermore, CRISPR/Cas9-mediated target gene editing is used as an important tool for functional analysis of plant and crop genes [29,30,31,32].
In this study, we aimed to understand the physiological role of OsPUB7 with U-box/ARM domain in rice, such as its role in abiotic stress resistance (drought and salinity stress), by generating OsPUB7 knockout individuals using CRISPR/Cas9 technology. Our study will provide a useful genetic resource for breeding crops that are resistant to abiotic stress caused by climate change.

2. Results

2.1. Analysis of the Genetic Relationship of OsPUB7 among Plant Species

Phylogenetic analysis was performed to obtain information regarding the genetic relationship between OsPUB7 and genes encoding E3 ubiquitin ligases in other plant species. U-box type E3 ubiquitin ligases and other E3 ubiquitin ligases encoded by 26 genes obtained from the NCBI database were analyzed using the MEGA7.0 program (Supplementary Table S1). As shown in Figure 1, the E3 ubiquitin ligases were classified according to the conserved protein domains in plant species. Among the genes used in the analysis, one closely related to OsPUB7 was found to be well conserved in the rice genome. Additionally, PUB genes belonging to wheat (Triticum aestivum), barley (Hordeum vulgare), and millet (Setaria italica) were classified in the same group. These results suggested that rice OsPUB7, encoding a domain of the U-box type E3 ubiquitin ligase, exhibits considerably high homology in various plant species, irrespective of plant origin (Figure 1).

2.2. Production of OsPUB7-Edited Lines and Analysis of Mutation Types

2.2.1. CRISPR/Cas9-Mediated Editing of OsPUB7

To construct gene-edited rice using CRISPR/Cas9 technology, sgRNA was created based on the second exon of OsPUB7 (Supplementary Figure S3). Among them, sgRNAs with a GC content of 50%, out of frame score of 84.7, and no mismatches in the genome were selected, which can reduce the off-target probability. The selected sgRNA was located around 113 bp from the 5′-end of the second exon and was analyzed as a target in this region (Figure 2A,B). A Ti-plasmid vector for plant gene editing was constructed using Cas9 under the control of the 35S promoter, sgRNA under the control of the OsU3 promoter, and bar gene as the selection marker. The complete vector, named CRISPR/Cas9-OsPUB7 (Figure 2C), was inoculated into the callus obtained from the seeds of the wild type (WT) plant (Dongjin) using an Agrobacterium-mediated method, and finally 31 transformants were obtained (Figure 3A). The regenerated plants were transplanted into pots and grown. To confirm the introduction of the T-DNA, polymerase chain reaction (PCR) analysis was performed using the genomic DNA of the T0 generation plants (Figure 3B). The results confirmed that all 31 regenerated plants obtained after stable introduction of the CRISPR/Cas9-OsPUB7 vector and transplantation in soil were transgenic.

2.2.2. Analysis of Mutation Types Using Next-Generation Sequencing (NGS)

Genetic mutation types of 31 regenerated individuals were analyzed using NGS technology. Among the 31 transgenic individuals, 15 individuals harbored mutations in OsPUB7, and their gene editing efficiency was 48%. The genotypes of the 15 gene-edited individuals were homozygous in 4 individuals, bi-allelic in 6 individuals, and heterozygous in 5 individuals (Supplementary Table S2). In homozygous individuals, 1 base was deleted or inserted, and in bi-allelic individuals, 1 base was inserted or 2, 22, or 25 deletions were found. In heterozygous individuals, 1 or 2 bases were inserted (Figure 4; Table 1). Mutations that occurred during double-strand breaking repair in plants appeared mainly as deletions and 1bp insertions. These mutations are stably inherited in offspring in both T2 and T3 generations according to the classical Mendelian model. When both copies of a target gene are mutated, as in homozygous or bi-allelic editing plants, the genotype is stable and resistant to further editing by CRISPR/Cas9. Bi-allelic gene targeting in seed-propagated plants has the advantage of saving time in correcting genetic alterations [33,34]. T1 seeds were harvested by cultivating gene-edited individuals in which the gene mutation type was confirmed. In the T1 generation, null individuals lacking the T-DNA were selected using bar screening. First, null individuals in which bar was completely removed were selected after treatment with 40 ppm basta solution using the bar strip test and bar PCR (Figure 5). Next, null individuals with a segregation ratio of 3:1 according to the law of segregation (using the chi-square (χ2) test) and single copy introduction were confirmed (Table 2).

2.3. Selection of Tolerant Gene-Edited Lines and Examination of Physiological Characteristics According to Abiotic Stress Treatment

2.3.1. Screening of Drought and Salinity Stress

The T2 OsPUB7 gene-edited null lines (PUB7-GE), K343-8-1 and K343-17-8, along with control varieties Dongjin (WT), Sangnambat (T), and Gaya (S), were subjected to drought stress treatment (Figure 6A). As a result of drought stress treatment, all Dongjin, Sangnambat, and Gaya plants stopped growing completely and dried with a damage level of 4 on the 5th day after the treatment. Most plants were permanently dried and did not recover, even after watering. K343-7-5 plants had a similar phenotype to that before stress treatment until the 3rd day of drought stress treatment, and the damage level was 3 to 4 on the 5th day after drought stress treatment. The K343-8-1 and K343-17-8 plants stopped growing and dried; almost all plants showed a damage level of 4 on the 5th day after the drought treatment. However, most of them showed normal growth in the recovery stage and recovered to damage level 2.
Salinity stress treatment was performed on T2 PUB7-GE K343-7-5, K343-8-1, and K343-17-8, as well as on Dongjin (WT), Cheongho (T), and Nampyeong (S) (Figure 6B). As a result of salinity stress treatment, leaf curling and growth retardation were observed in Dongjin, Cheongho, and Nampyeong from the third day after treatment. On the 5th day after treatment, all plants with a damage level of 4 or 3 stopped growing completely and dried. In the recovery stage, the control varieties showed a damage level of 5, and all plants died without recovering. However, in K343-7-5, K343-8-1, and K343-17-8, damage level 2 was observed from the 3rd day after salinity stress treatment, although the damage was lower than that observed in Dongjin and Cheongho. Contrary to that observed in the control varieties in the recovery stage, K343-8-1 and K343-17-8 showed almost normal growth, with a damage level of 2; K343-7-5 recovered to damage level 3. Overall, the OsPUB7 knockout mutants exhibited an increase in the abiotic stress tolerance of rice.

2.3.2. Determination of Physiological Characteristics and Analysis of mRNA Expression

To confirm the degree of tolerance to stress treatment in terms of physiological traits, T2 PUB7-GE were subjected to ion leakage and proline content analyses. The degree of ion leakage before and after drought stress treatment was analyzed in three replicates using 100 mg leaves. Results showed that ion leakage with cell damage had the following order: Gaya (S, 62.37%), Dongjin (WT, 44.19%), Sangnambat (T, 42.19%), K343-8-1 (24.03%), and K343-17-8 (12.86%) on the 5th day after drought stress treatment. In addition, T2 PUB7-GE maintained a low level of ion leakage during recovery, with a tendency similar to that observed before recovery (Figure 7A). To confirm the degree of tolerance to drought stress treatment, proline content analysis was performed on the 3rd day, when the phenotype started to appear (Figure 7B). The results showed that the proline content of WT was 1.001 mg/g, whereas that of K343-7-5 was 1.108 mg/g, that of K343-8-1 was 2.160 mg/g, and that of K343-17-8 was 2.011 mg/g, which were approximately 10–110% more than that of the WT.
The ion leakage analysis of T2 PUB7-GE and that of the salinity stress control varieties, Cheongho (T) and Nampyeong (S), showed that ion leakage followed the order of Dongjin (WT, 64.03%), Nampyeong (S, 55.51%), Cheongho (T, 33.33%), K343-8-1 (24.33%), K343-7-5 (23.19%), and K343-17-8 (11.71%) on the 5th day after salinity stress treatment (Figure 7C). To observe the degree of tolerance to salinity stress treatment, proline content was analyzed on the 3rd day, when the stress phenotype started to appear (Figure 7D). The proline content of WT was 0.789 mg/g, whereas that of K343-7-5 was 1.093 mg/g, that of K343-8-1 was 0.963 mg/g, and that of K343-17-8 was 0.994 mg/g, which were 22–38% more than that of Dongjin. Overall, considering the results of ion leakage and proline content analyses, we concluded that T2 PUB7-GE showed tolerance to drought and salinity stress.
To survey the expression level of OsPUB7 during abiotic stress treatment, qRT-PCR analysis was performed using the leaves 0, 3, 6, 12, 24, and 36 h after drought and salinity stress treatment (Figure 8A). At 3 h after drought stress treatment, T2 PUB7-GE appeared to show the highest gene expression, although finally, the expression did not differ significantly between WT and PUB7-GE or between individuals. This suggested that although OsPUB7 harbored a mutation introduced using the CRISPR/Cas9 technology, the transcriptional regulatory region of the gene functioned normally. Transcription of OsPUB7 occurs normally, although translation is believed to be restricted due to the creation of an early stop codon (due to a codon frame shift caused by an InDel). Similar to that observed for drought stress treatment, RNA expression did not differ significantly between WT and PUB7-GE, nor was any specific trend in expression change observed after salinity stress treatment (Figure 8B).

2.4. Analysis of Protein Structure in PUB7-GE

Before confirming the codon frame shift in T2 PUB7-GE, the genotypes of the lines showing tolerance to abiotic stress were analyzed. Results confirmed that all bi-allelic genotypes in the T0 generation were fixed as homozygous genotypes in the T2 generation, and the mutation types were identified as 2-bp deletion and 1-bp insertion, respectively (Figure 9).
Analysis of the amino acid sequence of three T2 PUB7-GE, which showed tolerance to drought and salinity stress, revealed that an early stop codon was generated at 159 bp (Figure 10B,D) and 162 bp (Figure 10C) from the 5′ end of the second exon, due to a codon frame shift. We assumed that the function of OsPUB7 was lost because translation was halted by the generation of the early stop codon. The protein structure of T2 PUB7-GE with amino acid mutations was analyzed, and it was observed that, compared to the WT, 308 amino acids (Figure 11B,D) and 307 amino acids (Figure 11C) of OsPUB7 were not translated. Hence, we believe that OsPUB7-edited null lines lost the E3 ubiquitin ligase activity, as normal protein synthesis did not occur (Figure 11).

2.5. Analysis of Protein–Protein Interactions

We aimed to identify the interaction between genes related to abiotic stress tolerance [35] and OsPUB7. Comparison of gene expression levels using T2 PUB7-GE mRNA derived after stress treatment revealed that the mRNA levels of the four genes, OsPUB23, OsPUB24, OsPUB66, and OsPUB67, increased (Figure 12). String DB-based protein interaction analysis between OsPUB7 and OsPUB families revealed that a 1-node network was formed with OsPUB66 (Figure 13). This indicated the functional association between OsPUB genes; it can be inferred that OsPUB7 acts as a negative regulator in WT and that OsPUB23, OsPUB24, OsPUB66, and OsPUB67 act to show tolerance to drought and salinity stress after knocking out OsPUB7.

2.6. Analysis of RNA-Sequencing Data of PUB7-GE

To identify the genes related to the abiotic stress tolerance pathway in PUB7-GE, we performed RNA-sequencing analysis in PUB7-GE K343-8-1 and compared it to that of the WT plants (Supplementary Table S3). We estimated the expression levels of transcripts using fragments per kilobase of exon per million mapped fragment (FPKM) values [36] (Supplementary Table S4). A gene was considered expressed in a sample if its FPKM was greater than 1. In total, 44,759 genes were expressed in at least one of four samples analyzed. Among the genes whose expression was confirmed, 27,831 genes were constitutively expressed in all samples, with a coefficient of variation <10%. Four samples were classified as harboring highly expressed genes under drought and salinity stress conditions. Differentially expressed gene (DEG) analysis revealed that the OsPUB7 knockout line contained 15,433 up-regulated genes and 17,162 down-regulated genes under drought stress treatment compared to that in the WT. Furthermore, K343-8-1 contained 15,654 up-regulated genes and 17,304 down-regulated genes under salinity stress treatment (Figure 14A). All of the differentially expressed genes were categorized into three main functional subgroups (Supplementary Figure S2), including Molecular Function (4789), Biological Process (4404), and Cellular Components (4401). For the Cellular Components category, the nucleus (GO:0005634) and integral component of membrane (GO:0016021) were the most highly represented categories. For the Molecular Function category, genes associated with metal ion binding (GO:0046872) and ATP binding (GO:0005524) have the highest fractions. Under the Biological Process category, proteins related to regulation of transcription (GO:0006355) and defense response (GO:0006952) were most frequent. After analyzing the down-regulated genes to determine the expression pattern of interacting genes in OsPUB7 and OsPUB gene families (Supplementary Figure S1), 5 out of 62 OsPUB genes were found to be down-regulated after abiotic stress treatment (Figure 14B). In addition, 6 of the 63 OsPUB genes were found to be up-regulated in OsPUB7-edited plants (Figure 14C). The up-regulation of OsPUB genes affected by the OsPUB7 knockout under drought and salinity stress treatment conditions suggests that OsPUB genes are closely related to abiotic stress responses. Based on these criteria, we finally identified six genes as the target candidates for the OsPUB7-related abiotic stress tolerance pathway.
To validate the results of sequencing analysis, qRT-PCR was performed using six OsPUB genes, which are abiotic stress-responsive genes selected from the in silico analysis of DEGs (Figure 15). Interestingly, RNA-sequencing analysis confirmed that, with the exception of OsPUB12, the remaining five genes showed significantly higher expression levels. The results of qRT-PCR of all six genes showed an up-regulated expression pattern in the OsPUB7-edited line.

2.7. Investigation of Agronomic Traits

To observe the growth of T3 PUB7-GE in a rice field, major agricultural traits such as plant height, culm length, panicle length, and number of tillers were investigated (Table 3, Figure 16). The plant height and number of tillers of K343-7-5 did not differ to those of WT, although culm length was approximately 4.8 cm shorter and panicle length was approximately 5 cm longer. K343-8-1 did not differ from the WT in terms of the four agricultural traits. The culm length, panicle length, and number of tillers of K343-17-8 did not differ from those of the WT, although plant height was approximately 4.8 cm longer than the WT. Despite some differences between OsPUB7-edited lines, overall, the agricultural traits were similar to those of WT, suggesting that gene editing using the CRISPR/Cas9 technology did not significantly affect genes other than OsPUB7.

3. Discussion

Target gene editing using the CRISPR/Cas9 technology has been applied to many crops for studying gene function and for improving agricultural traits such as resistance to various stress factors. CRISPR/Cas9-mediated gene editing is simple and easy to use [37], and consumers’ concerns regarding genetically modified organisms can also be avoided by removing foreign gene regions (T-DNA) introduced via selfing [38]. Hence, in this study, we have studied the physiological role of OsPUB7 in response to drought and salinity stress using CRISPR/Cas9-based genome-edited rice plants (PUB7-GE).
Drought and salinity stress treatment of these null lines revealed that in T2 PUB7-GE (K343-7-5 and K343-8-1), K343-17-8 showed tolerance to drought and salinity stress. In addition, analysis of the degree of abiotic stress and physiological tolerance performed by measuring ion leakage and proline content revealed that the tolerance pattern was the same as the observed phenotype. We observed that OsPUB7 acted as a negative regulator of drought and salinity stress and that the OsPUB7 knockout enhanced abiotic stress tolerance in rice. In the CRISPR/Cas9 system, translation of the target gene is halted when a premature termination codon (PTC) exists in the coding region of the gene to be corrected [39]. Sequencing analysis results showed that the mutation types in PUB7-GE used for drought and salinity stress treatment included a 2-bp deletion (K343-2-5 and K343-17-8) and 1-bp insertion (K343-8-1), respectively. Owing to a frame shift caused by an InDel in the coding region of OsPUB7, a PTC was formed at the beginning of the second exon sequence, as a result of which the PUB7-GE lines, K343-7-5, K343-18-8 and K343-8-1 were able to synthesize only 151 and 152 amino acids, respectively. PCT can occur during translation in the cytoplasm due to a codon frame shift caused by InDel in a gene. During this, when an mRNA binds to the ribosome and the stop codon is recognized, the ribosome degrades the mRNA via non-sense-mediated mRNA decay (NMD) to eliminate the potential risk associated with the accumulation of truncated proteins [40,41]. In the case of OsPUB7 gene editing lines, the frame shift-induced PCT is present at the beginning of the nucleotide sequence of the second exon of OsPUB7, which not only prevents subsequent amino acid synthesis, but also prevents production of the truncated protein. Protein–protein interaction analysis was performed with OsPUB23, OsPUB24, OsPUB66, OsPUB67, and OsPUB7, which are known to be involved in drought tolerance based on the results of transcriptome analysis under drought and salinity stress conditions. We confirmed that the expression of four genes increased in T2 PUB7-GE compared to that in the WT. The functions of Group II members in the PUB family of Arabidopsis have been extensively investigated in plant abiotic stress treatment. AtPUB22/AtPUB23 are negative regulators that mediate drought response in an ABA-independent pathway [42]. It has been reported that AtPUB30 negatively regulates salinity tolerance by promoting BRI1 KINASE INHIBITOR 1 (BKI1) degradation [43]. It can be assumed that the OsPUB7 gene in WT can act as a negative regulator of these four genes, and knockout of OsPUB7 affects the expression of several OsPUB genes, enhancing the drought and salinity tolerance of rice. To identify genes related to the abiotic stress tolerance pathway in the OsPUB7 knockout line, we performed RNA-sequencing analysis in PUB7-GE K343-8-1 and compared the results with those obtained using the WT plants. The up-regulated genes were analyzed to determine the expression pattern of interacting genes between the OsPUB7 and the OsPUB gene family. Results revealed that 6 out of 63 OsPUB genes were up-regulated after abiotic stress treatment. To validate the results of sequencing analysis, qRT-PCR was performed using six OsPUB genes, which are abiotic stress-responsive genes selected from the in silico analysis of DEGs. In qRT-PCR, five out of the six genes that showed up-regulated expression pattern in the OsPUB-edited line were found to be significantly highly expressed. This study provided clues regarding the interaction and linking role between OsPUB7 and a small number of OsPUB genes (OsPUB66, OsPUB23, OsPUB24, OsPUB67) under drought and salinity stress treatment. Among the up-regulated genes, two OsARM genes (OsPUB22 and OsPUB24) were differentially expressed in drought and salinity stress treatments. The ARM repeat (Armadillo repeat) gene has an ARM repeat domain consisting of one short α-helix and two long α-helices [20,44]. ARM genes have diverse functions such as protein degradation and signal transduction, nuclear transport, and cell adhesion [45]. In addition, it was confirmed that among up-regulated genes, genes encoding two protein kinases (OsPUB66 and OsPUB67) were over-represented. Most of these kinases are receptors such as Pelle kinase, and these receptors are the largest gene families in Arabidopsis and rice and respond to various abiotic and biotic stresses [46]. In rice, the receptor-like cytoplasmic kinase GROWTH UNDER DROUGHT KINASE (GUDK) has been shown to improve resistance to drought stress through activation of transcription factor APETALA2/ETHYLENE RESPONSE FACTOR OsAP37 [47]. In Arabidopsis, it was confirmed that water use efficiency is improved through the overexpression of Leucine-Rich Repeats Receptor-Like Kinase (LRR-RLK) genes [48]. Regulation of ubiquitination levels and stability in the changing expression of OsPUB7, and the signaling actions of sub-materials such as abscisic acid (ABA) and jasmonic acid (JA), should be explored in future work.
In conclusion, the editing of OsPUB7 using CRISPR/Cas9 technology resulted in the identification of several molecular mechanisms related to stress response that involve OsPUB gene-related pathways. These genes will be useful targets for both breeding and future research on drought tolerance/abiotic stress in rice. In addition, OsPUB7 can be expected to interact with the genes related to both biotic and abiotic stress, which may result in the degradation of damaged proteins to provide free amino acids for the synthesis of other proteins.

4. Materials and Methods

4.1. Plant Materials

In this study, Dongjinbyeo, a rice variety that is not tolerant to drought and salinity stress, was used as plant material for the generation of transgenic rice using CRISPR/Cas9-OsPUB7 vector. Seedlings were transplanted in a glass greenhouse and experimental farm at 30 × 15 cm spacing, with one seedling per hill arranged in an incomplete block design. The fertilizer N-P2O5-K2O were applied at a rate of 90-45-47 kg/ha. Cultivation management was performed following the rice cultivation standards adapted to the experimental area of Chungbuk National University.

4.2. Analysis of Phylogenetic Tree

To infer the relationship between PUB7 in rice and Arabidopsis and other crops, we calculated evolutionary divergence estimates and constructed a phylogenetic tree. The OsPUB7 gene was queried in RAP-DB(https://rapdb.dna.affrc.go.jp/index.html accessed on 3 August 2020), and the amino acid sequences of other plant gene homologues related to OsPUB7 were collected from NCBI (https://blast.ncbi.nlm.nih.gov accessed on 3 August 2020) using the OsPUB7 protein sequence. Corresponding genomic sequences and coding region sequences were blast searched in Gramene (http://www.gramene.org accessed on 3 August 2020), Rice Genome Annotation Project (rice.plantbiolog.msu.edu accessed on 3 August 2020), and NCBI (http://0-www-ncbi-nlm-nih-gov.brum.beds.ac.uk/gene accessed on 3 August 2020). The evolutionary history was inferred using the neighbor-joining method [49], and the phylogenetic tree was constructed and viewed with MEGA 7 software [50].

4.3. Selection of Target Sequences and Vector Construction

Target sites and sgRNAs for the 1000 bp exon region of OsPUB7 (Os04g0348400) adjacent to a protospacer-adjacent motif (PAM) were designed using the CRISPR RGEN tool (http://www.rgenome.net/ accessed on 2 September 2020) developed by the Hanyang University ([51]. According to the principles for designing target sequences in the CRISPR/Cas9 system, 20 base pairs excluding the PAM sequence of the selected sgRNA were designed as candidate target sequences as follows: target-up, 5′-ggcaGTGTGCACCACTGTGAGATCATGG-3′; target-down, 5′-aaacCCATGATCTCACAGTGGTGCACAC-3′. A vector for rice transformation was constructed to introduce the created OsPUB7 RGEN into the WT line (Oryza sativa L. cv. Dongjin-byeo). The pPZP-3′PINII-Bar vector was used as the parent vector for plant transformation. Cas9, under the CaMV 35S promoter, and sgRNA, under the OsU3 promoter, were ligated to this vector. Bar, under the control of the CaMV 35S promoter, was used as the selection marker; it was first linked to p35S at the HindIII/SacI site of the pPZP-3′PINII-Bar vector, which was then digested with BamHI/SacI and ligated with Cas9. Finally, AarI and XhoI sites were added to connect the sgRNA to the OsU3 promoter (Figure 2). The vector constructed was named CRISPR/Cas9-OsPUB7. The binary constructs were then introduced into Agrobacterium tumefaciens strain EHA105 via electroporation [52].

4.4. Agrobacterium-Mediated Transformation and Analysis of Mutation Types Using NGS

Agrobacterium-mediated transformation of embryogenic callus was performed by Lee et al. [52]. After 4 weeks of rooting, the regenerated rice plants were transferred to pots in a greenhouse maintained at 30 °C during the day and at 23 °C at night. Genomic DNA was extracted from the leaves of T0 individuals confirmed to harbor the T-DNA for gene mutation type analysis and mutant selection. Deep-sequencing was performed to amplify the genomic region containing the CRISPR/Cas9 target sites using specific primers adjacent to the designed target site [53]. PCR amplicons were sequenced by forming paired-end reads using MiniSeq (Illumina, San Diego, CA, USA). The NGS data obtained were analyzed using Cas-Analyzer (https://www.rgenome.net/cas-analyzer accessed on 5 April 2021) [54]. Reads that occurred only once were excluded to remove errors associated with amplification and sequencing. Insertion and deletion mutations were considered as mutations induced by Cas9. To select null lines without T-DNA, 40 ppm of basta solution was used to treat individuals that generated more than six leaves after the T1 generation. The basta solution was embedded in a gauze and placed on a part of the seedling leaves grown to the six-leaf stage, which were observed after 7 days to distinguish between resistance and sensitivity. These plants were also used to reconfirm whether they were resistant or sensitive using the bar strip test. T-DNA-free plants were selected on the basis of the death of the leaves after basta treatment and the absence of bands in the bar strip experiments. Finally, PCR analysis was performed on null lines clearly lacking the T-DNA using a bar primer, and individuals with no band in 1% agarose gel electrophoresis were selected as the final null T1 lines.

4.5. Screening for Drought and Salt Stress Resistance

For screening drought stress tolerance, 4-week-old T2 PUB7-GE lines, along with drought-tolerant Sangnambat and susceptible Gaya, were grown in the greenhouse and transplanted to soil with 5% moisture content. Then the screening was performed without watering until the appearance of the drought phenotype. The soil moisture content was measured using a soil moisture meter (ProCheck, Armidale, NSW). The screening was stopped when the susceptible Gaya line showed clear signs of damage, following which it was irrigated to recover for 2 weeks.
For the screening of salinity stress tolerance of the GE null lines, 4-week-old plants of salt-tolerant Cheongho and susceptible Nampyeong were treated with 250 mM sodium chloride (NaCl). NaCl solution was replaced once every two days to prevent the reduction of salt concentration due to absorption and evaporation during treatment. The screening was stopped when the susceptible Nampyeong showed clear signs of damage. The plant was then irrigated with clean water and allowed to recover for 2 weeks. Leaf rolling and leaf drying were evaluated [55] when the susceptible variety showed symptoms under stress. The seedlings were evaluated using the scales described in IRRI [56]: 1 = normal growth and no leaf symptoms, 2 = nearly normal growth although tips of few leaves are discolored and rolled, 3 = growth is severely retarded with most leaves rolled and only a few elongating, 4 = complete cessation of growth along with drying of some leaves and death of some plants, and 5 = almost all plants died (Figure 17).

4.6. Measurement of the Physiological Parameters of Transgenic Rice

The proline content in leaves was estimated according to the method used by Abdula et al. [57], with some modifications. Briefly, 0.1 g of rice leaves was ground with 1 mL MCW buffer (MeOH:chloroform:water = 12:5:1), and the homogenate was centrifuged at 4 °C and 13,000 rpm for 10 min. The extract was mixed with 200 μL of 2% ninhydrin reagent and 200 μL of acetic acid. The mixture was mixed using a vortex mixer and boiled at 100 °C for 45 min. The boiled sample was then frozen on ice, combined with 1000 μL toluene, and then left to stand for 5–10 min. Absorbance of the reddish pink upper phase was recorded at 520 nm against a toluene blank. L-proline was used as the standard solution.
Ion leakage was measured before and after autoclaving using 0.1 g leaves (1 cm in diameter). The leaf discs were placed in closed tubes containing 10 mL distilled water and incubated at 32 °C for 2 h. Subsequently, the initial electrical conductivity of the solution (EC1) was determined using a multi-range EC meter (HANNA Instruments, Woonsocket, RI, USA). The samples were autoclaved to release all electrolytes and cooled down to 25 °C, and their final electrical conductivity (EC2) was measured. The electrolyte leakage (EL) was calculated as EL = (EC1/EC2) × 100 (%).

4.7. mRNA Expression Analysis

Leaf samples were collected at 0, 3, 6, 12, and 24 h after treatment and at 72 h after recovery. Total RNA was extracted from leaf tissues using the RNeasy plant mini kit (QIAGEN, Redwood City, USA) according to the manufacturer’s instructions. The relative purity and concentration of RNA were estimated using NanoDrop One (Thermo Fisher Scientific, Wilmington, USA) and stored in a −80 °C freezer. The first-strand cDNAs were synthesized using Oligo (dT)20 primer and ReverTra AceTM qPCR RT master mix (TOYOBO, Osaka, Japan). Actin primers were used as the internal control to normalize the results of the real-time qRT-PCR. All samples were analyzed in triplicate to increase the accuracy of the experiment.

4.8. Analysis of OsPUB7-GE

Codon changes caused by frame shift mutations were confirmed from the nucleotide and amino acid sequence analysis of PUB7-GE-harboring InDels. For the amino acid sequence, the translator tool (https://web.expasy.org/translate/ accessed on 5 December 2021) provided by Expasy was used. Based on the amino acid sequence of PUB7-GE, the protein structure was confirmed using SBI’s modeling program (https://swissmodel.expasy.org/interactive accessed on 5 December 2021).

4.9. Protein–Protein Interaction Network Analysis

Protein–protein interactions were predicted using the STRING database (https://string-db.org accessed on 5 August 2021) to investigate interaction between the U-box type E3 ubiquitin ligase genes and OsPUB7 in rice. The OsPUB7 amino acid sequence obtained from RAP-DB (https://rapdb.dna.affrc.go.jp accessed on 5 August 2021) was used in the amino acid sequence column, and orthologs of Oryza sativa were selected as references. After completing the BLAST step, the network was constructed using the highest score gene.

4.10. RNA-Sequencing Analysis

Leaves of PUB7-GE and WT lines were collected after drought stress treatment. RNA libraries were sequenced using the Illumina High-Seq 2500 platform, provided by a commercial service provider (Theragen Bio, Seongnam, Korea). The low-quality bases (Q < 15) were trimmed at both ends of the sequence using a customized program, and the adapter was trimmed using Cutadapt. The high-quality reads were subsequently aligned to the IRGSP-1.0 reference genome sequence using Bowtie and TopHat. The expression levels of each transcript were expressed as FPKM values calculated based on the number of mapped reads. All DEGs were determined using Cufflinks program v.2.0.1 (http://cufflinks.cbcb.umd.edu accessed on 1 January 2022) [58]. All DEGs were identified using DEGseq and categorized according to the Gene Ontology framework using DAVID bioinformatics resources (ver. 6.8) (https://david,ncifcrf.gov/tools.jsp accessed on 1 January 2022).

4.11. Statistical Analysis

Data were analyzed by one-way analysis of variance (ANOVA) using Statistical Analysis System (SAS version 9.4). Values are mean ± SE (n = 3) and statistical significance was set to p < 0.05 according to Duncan’s multiple range test.

Supplementary Materials

The following supporting information can be downloaded at: https://0-www-mdpi-com.brum.beds.ac.uk/article/10.3390/ijms24065338/s1.

Author Contributions

Experimental analysis and writing of the article, M.-S.K. and S.-R.K.; review and editing, K.-K.K., Y.J.J. and Y.-J.L.; supervision, idea development, and editing, Y.-G.C. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by Chungbuk National University Korea National University Development Project (2022), Republic of Korea.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

This research was supported by Chungbuk National University Korea National University Development Project (2022), Republic of Korea.

Conflicts of Interest

The authors have no conflict of interest relevant to this study to disclose.

References

  1. Amm, I.; Sommer, T.; Wolf, D.H. Protein quality control and elimination of protein waste: The role of the ubiqui-tin-proteasome system. Biochim. Biophys. Acta Mol. Cell Res. 2014, 1843, 182–196. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Doyle, S.M.; Genest, O.; Wickner, S. Protein rescue from aggregates by powerful molecular chaperone ma-chines. Nat. Rev. Mol. Cell Biol. 2013, 14, 617–629. [Google Scholar] [CrossRef] [PubMed]
  3. Nelson, C.; Millar, A. Protein turnover in plant biology. Nat. Plants 2015, 1, 15017. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Santner, A.; Estelle, M. The ubiquitin-proteasome system regulates plant hormone signaling. Plant J. 2010, 61, 1029–1040. [Google Scholar] [CrossRef] [PubMed]
  5. Ciechanover, A. The ubiquitin-proteasome proteolytic pathway. Cell 1994, 79, 13–21. [Google Scholar] [CrossRef] [PubMed]
  6. Peng, J.; Schwartz, D.; Elias, J.E.; Thoreen, C.C.; Cheng, D.; Marsischky, G.; Roelofs, J.; Finley, D.; Gygi, S.P. A proteomics approach to understanding protein ubiquitination. Nat. Biotechnol. 2003, 21, 921–926. [Google Scholar] [CrossRef]
  7. Sadowski, M.; Suryadinata, R.; Tan, A.R.; Roesley, S.N.A.; Sarcevic, B. Protein monoubiquitination and polyubiquitination generate structural diversity to control distinct biological processes. IUBMB Life 2012, 64, 136–142. [Google Scholar] [CrossRef]
  8. Kirisako, T.; Kamei, K.; Murata, S.; Kato, M.; Fukumoto, H.; Kanie, M.; Sano, S.; Tokunaga, F.; Tanaka, K.; Iwai, K. A ubiquitin ligase complex assembles linear polyubiquitin chains. EMBO J. 2016, 25, 4877–4887. [Google Scholar] [CrossRef]
  9. Berndsen, C.; Wolberger, C. New insights into ubiquitin E3 ligase mechanism. Nat. Struct. Mol. Biol. 2014, 21, 301–307. [Google Scholar] [CrossRef]
  10. Ciechanover, A. The ubiquitin–proteasome pathway: On protein death and cell life. EMBO J. 1998, 17, 7151–7160. [Google Scholar] [CrossRef] [Green Version]
  11. Pickart, C.M.; Eddins, M.J. Ubiquitin: Structures, functions, mechanisms. Biochim. Biophys. Acta-Mol. Cell Res. 2004, 1695, 55–72. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Lecker, S.H.; Goldberg, A.L.; Mitch, W.E. Protein degradation by the ubiquitin-proteasome pathway in normal and disease states. J. Am. Soc. Nephrol. 2006, 177, 1807–1819. [Google Scholar] [CrossRef] [PubMed]
  13. Bae, H.; Kim, W.T. Classification and interaction modes of 40 rice E2 ubiquitin-conjugating enzymes with 17 rice ARM-U-box E3 ubiquitin ligases. Biochem. Biophys. Res. Commun. 2014, 444, 575–580. [Google Scholar] [CrossRef] [PubMed]
  14. Patterson, C. A new gun in town: The U box is a ubiquitin ligase domain. Sci. STKE 2002, 2002, PE4. [Google Scholar] [CrossRef] [PubMed]
  15. Aravind, L.; Koonin, E.V. The U box is a modified RING finger: A common domain in ubiquitination. Curr. Biol. 2000, 10, R132–R134. [Google Scholar] [CrossRef] [Green Version]
  16. Smalle, J.; Vierstra, R.D. The ubiquitin 26S proteasome proteolytic pathway. Annu. Rev. Plant Biol. 2004, 55, 555–559. [Google Scholar] [CrossRef] [PubMed]
  17. Stone, S.L.; Callis, J. Ubiquitin ligases mediate growth and development by promoting protein death. Curr. Opin. Plant Biol. 2007, 10, 624–632. [Google Scholar] [CrossRef]
  18. Koegl, M.; Hoppe, T.; Schlenker, S.; Ulrich, H.D.; Mayer, T.U.; Jentsch, S. A novel ubiquitination factor, E4, is involved in multiubiquitin chain assembly. Cell 1999, 96, 635–644. [Google Scholar] [CrossRef] [Green Version]
  19. Hatakeyama, S.; Yada, M.; Matsumoto, M.; Ishida, N.; Nakayama, K.I. U Box Proteins as a New Family of Ubiquitin-Protein Ligases. J. Biol. Chem. 2001, 276, 33111–33120. [Google Scholar] [CrossRef] [Green Version]
  20. Mudgil, Y.; Shiu, S.H.; Stone, S.L.; Salt, J.N.; Goring, D.R. A Large Complement of the Predicted Arabidopsis ARM Repeat Proteins Are Members of the U-Box E3 Ubiquitin Ligase Family. Plant Physiol. 2004, 134, 59–66. [Google Scholar] [CrossRef] [Green Version]
  21. Yee, D.; Goring, D.R. The diversity of plant U-box E3 ubiquitin ligases: From upstream activators to down-stream target substrates. J. Exp. Bot. 2009, 60, 1109–1121. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Kong, L.Y.; Cheng, J.K.; Zhu, Y.J.; Ding, Y.L.; Meng, J.J.; Chen, Z.Z.; Xie, Q.; Guo, Y.; Li, J.; Yang, S.; et al. Degradation of the ABA co-receptor ABI1 by PUB12/13 U-box E3ligases. Nat. Commun. 2015, 6, 8630–8642. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Seo, D.H.; Ahn, M.Y.; Park, K.Y.; Kim, E.Y.; Kim, W.T. TheN-terminal UND motif of the Arabidopsis U-Box E3 Ligase PUB18 Is critical forthe negative regulation of ABA-mediated stomatal movement anddetermines its ubiquitination specificity for Exocyst Subunit Exo70B1. Plant Cell 2016, 28, 2952–2973. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Zhu, Y.F.; Li, Y.B.; Fei, F.; Wang, Z.K.; Wang, W.; Cao, A.Z.; Liu, Y.; Han, S.; Xing, L.P.; Wang, H.Y.; et al. E3 ubiquitin ligase gene CMPG1-V from Haynaldia villosa L. contributes to powdery mildew resistance in common wheat (Triticum aestivum L.). Plant J. 2015, 84, 154–168. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Liao, D.; Cao, Y.; Sun, X.; Espinoza, C.; Nguyen, C.T.; Liang, Y.; Stacey, G. Arabidopsis E3 ubiquitin ligase PLANT U-BOX13 (PUB13) regulates chitin receptor Lysin motif receptor kinase5 (LYK5) protein abundance. New Phytol. 2017, 214, 1646–1656. [Google Scholar] [CrossRef] [Green Version]
  26. Zeng, L.R.; Park, C.H.; Venu, R.C.; Gough, J.; Wang, G.L. Classification, expression pattern, and E3 ligase activity assay of rice U-box-containing proteins. Mol. Plant 2008, 1, 800–815. [Google Scholar] [CrossRef] [Green Version]
  27. Park, J.J.; Yi, J.; Yoon, J.; Cho, L.H.; Ping, J.; Jeong, H.J.; Cho, S.K.; Kim, W.T.; An, G. OsPUB15, an E3 ubiquitin ligase, functions to reduce cellular oxidative stress during seedling establishment. Plant J. Cell Mol. Biol. 2011, 65, 194–205. [Google Scholar] [CrossRef]
  28. Lv, Q.; Li, X.; Jin, X.; Sun, Y.; Wu, Y.; Wang, W.; Huang, J. Rice OsPUB16 modulates the ‘SAPK9-OsMADS23-OsAOC’ pathway to reduce plant water-deficit tolerance by repressing ABA and JA biosynthesis. PLoS Genet. 2022, 18, e1010520. [Google Scholar] [CrossRef]
  29. Cong, J.; Ran, F.A.; Cox, D.; Lin, S.; Barretto, R.; Habib, N.; Hsu, P.D.; Wu, X.; Jiang, W.; Marraffini, L.A.; et al. Multiplex Genome Engineering Using CRISPR/Cas Systems. Physiol. Behav. 2013, 176, 139–148. [Google Scholar] [CrossRef] [Green Version]
  30. Doudna, J.A.; Charpentier, E. The new frontier of genome engineering with CRISPR-Cas9. Science 2014, 346, 1258096. [Google Scholar] [CrossRef]
  31. Jinek, M.; Jiang, F.; Taylor, D.W.; Sternberg, S.H.; Kaya, E.; Ma, E.; Anders, C.; Hauer, M.; Zhou, K.; Lin, S.; et al. Structures of Cas9 endonucleases reveal RNA-mediated conformational activation. Science 2014, 343, 1247997. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  32. Sedeek, K.E.M.; Mahas, A.; Mahfouz, M. Plant genome engineering for targeted improvement of crop traits. Front. Plant Sci. 2019, 10, 1–16. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Feng, Z.; Mao, Y.; Xu, N.; Zhang, B.; Wei, P.; Yang, D.L.; Wang, Z.; Zhang, Z.; Zheng, R.; Yang, L.; et al. Multigeneration analysis reveals the inheritance, specificity, and patterns of CRISPR/Cas-induced gene modifications in Arabidopsis. Proc. Natl. Acad. Sci. USA 2014, 111, 4632–4637. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Endo, M.; Mikami, M.; Toki, S. Biallelic Gene Targeting in Rice. Plant Physiol. 2016, 170, 667–677. [Google Scholar] [CrossRef] [Green Version]
  35. Kim, M.S.; Kang, K.K.; Cho, Y.G. Molecular and Functional Analysis of U-box E3 Ubiquitin Ligase Gene Family in Rice (Oryza sativa). Int. J. Mol. Sci. 2021, 22, 12088. [Google Scholar] [CrossRef]
  36. Trapnell, C.; Williams, B.A.; Pertea, G.; Mortazavi, A.; Kwan, G.; van Baren, M.J.; Salzberg, S.L.; Barbara, J.W.; Lior, P. Transcript as-sembly and quantification by RNA-Seq reveals unannotated transcripts and isoform switching during cell differentiation. Nat. Biotechnol. 2010, 28, 511–515. [Google Scholar] [CrossRef] [Green Version]
  37. Adli, M. The CRISPR tool kit for genome editing and beyond. Nat. Commun. 2018, 9, 1911. [Google Scholar] [CrossRef] [Green Version]
  38. Waltz, E. Gene-edited CRISPR mushroom escapes US regulation. Nature 2016, 532, 293. [Google Scholar] [CrossRef] [Green Version]
  39. Kapahnke, M.; Banning, A.; Tikkanen, R. Random splicing of several exons caused by a single base change in the target exon of CRISPR/cas9 mediated gene knockout. Cells 2016, 5, 45. [Google Scholar] [CrossRef] [Green Version]
  40. Maquat, L.E. Nonsense-mediated mRNA decay: Splicing, translation and mRNP dynamics. Nat. Rev. Mol. Cell Biol. 2004, 5, 89–99. [Google Scholar] [CrossRef]
  41. Peccarelli, M.; Kebaara, B.W. Regulation of natural mRNAs by the nonsense-mediated mRNA decay pathway. Eukaryot. Cell 2014, 13, 1126–1135. [Google Scholar] [CrossRef] [PubMed]
  42. Seo, D.H.; Ryu, M.Y.; Jammes, F.; Hwang, J.H.; Turek, M.; Kang, B.G.; Kwak, J.M.; Kim, W.T. Roles of four Arabidopsis U-box E3 ubiquitin ligases in negative regulation of abscisic acid-mediated drought stress responses. Plant Physiol. 2012, 160, 556–568. [Google Scholar] [CrossRef] [Green Version]
  43. Zhang, M.; Zhao, J.; Li, L.; Gao, Y.; Zhao, L.; Patil, S.B.; Fang, J.; Zhang, W.; Yang, Y.; Li, M.; et al. The Arabidopsis U-box E3 ubiquitin ligase PUB30 negatively regulates salt tolerance by facilitating BRI1 kinase inhibitor 1 (BKI1) degradation. Plant Cell Environ. 2017, 40, 2831–2843. [Google Scholar] [CrossRef] [PubMed]
  44. Sharma, M.; Pandey, G.K. Expansion and function of repeat domain proteins during stress and development in plants. Front. Plant Sci. 2016, 6, 1218. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Sharma, M.; Singh, A.; Shankar, A.; Pandey, A.; Baranwal, V.; Kapoor, S.; Tyagi, A.K.; Pandey, G.K. Comprehensive expression analysis of rice Armadillo gene family during abiotic stress and development. DNA Res. Int. J. Rapid Publ. Rep. Genes Genomes 2014, 21, 267–283. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  46. Lehti-Shiu, M.; Zou, C.; Hanada, K.; Shiu, S.H. Evolutionary history and stress regulation of plant receptor-like kinase/pelle genes. Plant Physiol. 2009, 150, 12–26. [Google Scholar] [CrossRef] [Green Version]
  47. Ramegowda, V.; Krishnan, A.; Pereira, A. Rice GROWTH UNDER DROUGHT KINASE is required for drought tolerance and grain yield under normal and drought stress conditions. Plant Physiol. 2014, 166, 1634–1645. [Google Scholar] [CrossRef] [Green Version]
  48. Xing, H.T.; Guo, P.; Xia, X.L.; Yin, W.L. PdERECTA, a leucine-rich repeat receptor-like kinase of poplar, confers enhanced water use efficiency in Arabidopsis. Planta 2011, 234, 229–241. [Google Scholar] [CrossRef]
  49. Saitou, N.; Nei, M. The neighbor-joining method: A new method for reconstructing phylogenetic trees. Mol. Biol. Evol. 1987, 4, 406–425. [Google Scholar]
  50. Kumar, S.; Stecher, G.; Tamura, K. MEGA7: Molecular Evolutionary Genetics Analysis version 7.0 for bigger datasets. Mol. Biol. Evol. 2016, 33, 1870–1874. [Google Scholar] [CrossRef] [Green Version]
  51. Park, J.; Bae, S.; Kim, J. Cas-Designer: A web-based tool for choice of CRISPR-Cas9 target sites. Bioinformatics 2015, 31, 4014–4016. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  52. Lee, H.J.; Abdula, S.E.; Jee, M.G.; Jang, D.W.; Cho, Y.G. High-efciency and Rapid Agrobacterium-mediated genetic trans-formation method using germinating rice seeds. J. Plant Biotechnol. 2011, 38, 251–257. [Google Scholar] [CrossRef]
  53. Kim, H.; Choi, J.; Won, K. A stable DNA-free screening system for CRISPR/RNPs-mediated gene editing in hot and sweet cultivars of Capsicum annuum. BMC Plant Biol. 2020, 20, 449. [Google Scholar] [CrossRef] [PubMed]
  54. Park, J.; Lim, K.; Kim, J.; Bae, S. Cas-analyzer: An online tool for assessing genome editing results using NGS data. Bioinformatics 2017, 33, 286–288. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  55. International Rice Research Institute (IRRI). Standard Evaluation System for Rice; IRRI: Los Banos, Philippines, 2002. [Google Scholar]
  56. IRRI. Standard Evaluation System for Rice 4th Edition; IRRI, INGER Genetic Resources Center: Manila, Philippines, 1996. [Google Scholar]
  57. Abdula, S.E.; Lee, H.; Ryu, H.; Kang, K.; Nou, I.S.; Sorrells, M.E.; Cho, Y. Overexpression of BrCIPK1 Gene En-hances Abiotic Stress Tolerance by Increasing Proline Biosynthesis in Rice. Plant Mol. Biol. Report. 2015, 34, 501–511. [Google Scholar] [CrossRef]
  58. Trapnell, C.; Roberts, A.; Goff, L.; Pertea, G.; Kim, D.; Kelley, D.R.; Pimentel, H.; Salzberg, S.L.; Rinn, J.L.; Pachter, L. Differential gene and transcript expression analysis of RNA-seq experiments with TopHat and Cufflinks. Nat. Protoc. 2012, 7, 562. [Google Scholar] [CrossRef] [Green Version]
Ijms 24 05338 g001 550
Figure 1. Phylogenetic analysis of OsPUB7 in various plant species. The phylogenetic tree was constructed using the neighbor-joining method of MEGAX. The numbers represent bootstrap values from 1000 replicates.
Figure 1. Phylogenetic analysis of OsPUB7 in various plant species. The phylogenetic tree was constructed using the neighbor-joining method of MEGAX. The numbers represent bootstrap values from 1000 replicates.
Ijms 24 05338 g001
Ijms 24 05338 g002 550
Figure 2. Designing of sgRNA based on CRISPR RGEN Tool (http://www.rgenome.net/ accessed on 2 September 2020). (A) Selection of the target region of OsPUB7 using the CRISPR RGEN tool. (B) The position of the sgRNA target sites in the OsPUB7 sequence. (C) Construction of the Ti-plasmid vector of the sgRNA region for CRISPR/Cas9-mediated mutagenesis of OsPUB7 in rice. LB, left border; RB, right border; Nos, nopaline synthase terminator; Bar, basta resistance gene; sgRNA, single guide RNA; OsU3-P, Oryza sativa U3 promoter driven sgRNA cassette; Cas9, human codon-optimized Cas9 expressing cassette.
Figure 2. Designing of sgRNA based on CRISPR RGEN Tool (http://www.rgenome.net/ accessed on 2 September 2020). (A) Selection of the target region of OsPUB7 using the CRISPR RGEN tool. (B) The position of the sgRNA target sites in the OsPUB7 sequence. (C) Construction of the Ti-plasmid vector of the sgRNA region for CRISPR/Cas9-mediated mutagenesis of OsPUB7 in rice. LB, left border; RB, right border; Nos, nopaline synthase terminator; Bar, basta resistance gene; sgRNA, single guide RNA; OsU3-P, Oryza sativa U3 promoter driven sgRNA cassette; Cas9, human codon-optimized Cas9 expressing cassette.
Ijms 24 05338 g002
Ijms 24 05338 g003 550
Figure 3. Development of transgenic rice plants with CRISPR/Cas9-OsPUB7 vector using Agrobacterium-mediated transformation. (A) Dongjin (WT) transformed with pPZP-3′PinII-Bar vector containing the CRISPR/Cas9:sgRNA gene: (a) seeds plated on 2N6 media; (b) callus formation; (c) Agrobacterium co-culture in 2N6 AS media; (d,e) multi-shoot differentiation; (f) regenerated plants in rooting medium; (g,h) acclimation in tissue culture room and planted in soil. (B) Detection of sgRNA/Cas9-mediated DNA modifications using PCR analysis of bar and nos terminator regions. L, DNA ladder; N, negative control (sterile water); P, positive control (plasmid vector).
Figure 3. Development of transgenic rice plants with CRISPR/Cas9-OsPUB7 vector using Agrobacterium-mediated transformation. (A) Dongjin (WT) transformed with pPZP-3′PinII-Bar vector containing the CRISPR/Cas9:sgRNA gene: (a) seeds plated on 2N6 media; (b) callus formation; (c) Agrobacterium co-culture in 2N6 AS media; (d,e) multi-shoot differentiation; (f) regenerated plants in rooting medium; (g,h) acclimation in tissue culture room and planted in soil. (B) Detection of sgRNA/Cas9-mediated DNA modifications using PCR analysis of bar and nos terminator regions. L, DNA ladder; N, negative control (sterile water); P, positive control (plasmid vector).
Ijms 24 05338 g003
Ijms 24 05338 g004 550
Figure 4. Identification of mutation genotypes in OsPUB7 generated using CRISPR/Cas9 technology using NGS analysis in Dongjin (WT) and T0 plants. Green indicates the target sequence, and blue indicates the PAM region. The InDel sequence of the target locus is shown in red. Deletion in the analysis region is indicated by “-”, and insertion is indicated by an insertion of nucleotide such as “T” or “G”. The minus (−) and plus (+) signs indicate the number of nucleotides deleted and inserted at the target sites, respectively. WT, wild type sequence with no mutations.
Figure 4. Identification of mutation genotypes in OsPUB7 generated using CRISPR/Cas9 technology using NGS analysis in Dongjin (WT) and T0 plants. Green indicates the target sequence, and blue indicates the PAM region. The InDel sequence of the target locus is shown in red. Deletion in the analysis region is indicated by “-”, and insertion is indicated by an insertion of nucleotide such as “T” or “G”. The minus (−) and plus (+) signs indicate the number of nucleotides deleted and inserted at the target sites, respectively. WT, wild type sequence with no mutations.
Ijms 24 05338 g004
Ijms 24 05338 g005 550
Figure 5. Selection of T-DNA-free null lines. Bar screening after glufosinate/basta treatment on leaf tip, and bar strip analysis in T1 generation. (A) Screening of OsPUB7 transgenic T0 mutants after treatment with 40 ppm basta. (B) Detection of the selectable marker gene using PAT/bar test strip and bar PCR: (a) selection of OsPUB7 null lines after PAT/bar strip test; (b) selection of OsPUB7 transgenic T1 plants after bar PCR. The minus (−) and plus (+) signs indicate resistance to bar.
Figure 5. Selection of T-DNA-free null lines. Bar screening after glufosinate/basta treatment on leaf tip, and bar strip analysis in T1 generation. (A) Screening of OsPUB7 transgenic T0 mutants after treatment with 40 ppm basta. (B) Detection of the selectable marker gene using PAT/bar test strip and bar PCR: (a) selection of OsPUB7 null lines after PAT/bar strip test; (b) selection of OsPUB7 transgenic T1 plants after bar PCR. The minus (−) and plus (+) signs indicate resistance to bar.
Ijms 24 05338 g005
Ijms 24 05338 g006 550
Figure 6. Phenotype of responses induced by drought and salinity stress treatment in T2 PUB7-GE. (A) Phenotype of PUB7-GE in the presence of 5% soil moisture content. WT, Dongjin; T, Sangnambat; S, Gaya. (B) Phenotype of PUB7-GE in the presence of 250 mM NaCl. WT, Dongjin; T, Cheongho; S, Nampyeong.
Figure 6. Phenotype of responses induced by drought and salinity stress treatment in T2 PUB7-GE. (A) Phenotype of PUB7-GE in the presence of 5% soil moisture content. WT, Dongjin; T, Sangnambat; S, Gaya. (B) Phenotype of PUB7-GE in the presence of 250 mM NaCl. WT, Dongjin; T, Cheongho; S, Nampyeong.
Ijms 24 05338 g006
Ijms 24 05338 g007 550
Figure 7. Determination of physiological traits in T2 PUB7-GE during drought and salinity stress treatment. (A) Ion leakage levels (%) were measured after 5 days of drought stress treatment and after 3 weeks of recovery in WT and PUB7-GE. (B) Proline contents (mg/g FW) of WT and PUB7-GE after drought treatment (C) Ion leakage levels (%) were measured after 5 days of salinity stress treatment and after 2 weeks of recovery in WT and PUB7-GE. (D) Proline contents (mg/g) of WT and PUB7-GE after salinity treatment. Three replicates were used, and data are shown with SE.
Figure 7. Determination of physiological traits in T2 PUB7-GE during drought and salinity stress treatment. (A) Ion leakage levels (%) were measured after 5 days of drought stress treatment and after 3 weeks of recovery in WT and PUB7-GE. (B) Proline contents (mg/g FW) of WT and PUB7-GE after drought treatment (C) Ion leakage levels (%) were measured after 5 days of salinity stress treatment and after 2 weeks of recovery in WT and PUB7-GE. (D) Proline contents (mg/g) of WT and PUB7-GE after salinity treatment. Three replicates were used, and data are shown with SE.
Ijms 24 05338 g007
Ijms 24 05338 g008 550
Figure 8. RNA expression patterns of OsPUB7 in response to drought and salinity stress treatment analyzed using qRT-PCR. Relative gene expression level was normalized to that of the rice actin gene, ACT1. Values are presented as mean ± standard error. Three independent biological replicates were used. Bars represent standard deviations (SD) of three technical replicates. Significant differences according to t-test are indicated using asterisks (*, p < 0.05; **, p < 0.01).
Figure 8. RNA expression patterns of OsPUB7 in response to drought and salinity stress treatment analyzed using qRT-PCR. Relative gene expression level was normalized to that of the rice actin gene, ACT1. Values are presented as mean ± standard error. Three independent biological replicates were used. Bars represent standard deviations (SD) of three technical replicates. Significant differences according to t-test are indicated using asterisks (*, p < 0.05; **, p < 0.01).
Ijms 24 05338 g008
Ijms 24 05338 g009 550
Figure 9. Identification of genotypes using NGS analysis in T2 PUB7-GE. The minus (−) and plus (+) signs indicate the number of nucleotides deleted and inserted at the target sites, respectively. WT, wild type sequence with no mutations.
Figure 9. Identification of genotypes using NGS analysis in T2 PUB7-GE. The minus (−) and plus (+) signs indicate the number of nucleotides deleted and inserted at the target sites, respectively. WT, wild type sequence with no mutations.
Ijms 24 05338 g009
Ijms 24 05338 g010 550
Figure 10. Amino acid sequence based on nucleotide sequence of PUB7-GE. Green indicates the start codon, red indicates the stop codon, and gray indicates the InDel mutation type. (A) OsPUB7 wild type. (B) PUB7-GE K343-7-5. (C) PUB7-GE K343-8-1. (D) PUB7-GE K343-17-8.
Figure 10. Amino acid sequence based on nucleotide sequence of PUB7-GE. Green indicates the start codon, red indicates the stop codon, and gray indicates the InDel mutation type. (A) OsPUB7 wild type. (B) PUB7-GE K343-7-5. (C) PUB7-GE K343-8-1. (D) PUB7-GE K343-17-8.
Ijms 24 05338 g010
Ijms 24 05338 g011 550
Figure 11. 3D protein structure based on amino acid sequences of the gene-edited PUB7-GE. (A) OsPUB7 wild type. (B) PUB7-GE K343-7-5. (C) PUB7-GE K343-8-1. (D) PUB7-GE K343-17-8.
Figure 11. 3D protein structure based on amino acid sequences of the gene-edited PUB7-GE. (A) OsPUB7 wild type. (B) PUB7-GE K343-7-5. (C) PUB7-GE K343-8-1. (D) PUB7-GE K343-17-8.
Ijms 24 05338 g011
Ijms 24 05338 g012 550
Figure 12. Expression analysis of OsPUB23, OsPUB24, OsPUB66, and OsPUB67 after drought and salinity stress treatment in T2 PUB7-GE. (A) Confirmation of the mRNA expression of stress-related genes in PUB7-GE 12 h after stress treatment using agarose gel electrophoresis. (B) Band intensity was quantified using ImageJ.
Figure 12. Expression analysis of OsPUB23, OsPUB24, OsPUB66, and OsPUB67 after drought and salinity stress treatment in T2 PUB7-GE. (A) Confirmation of the mRNA expression of stress-related genes in PUB7-GE 12 h after stress treatment using agarose gel electrophoresis. (B) Band intensity was quantified using ImageJ.
Ijms 24 05338 g012
Ijms 24 05338 g013 550
Figure 13. Result of protein–protein interaction analysis obtained using String DB analysis.
Figure 13. Result of protein–protein interaction analysis obtained using String DB analysis.
Ijms 24 05338 g013
Ijms 24 05338 g014 550
Figure 14. Comparative analysis of differentially expressed genes (DEGs) after drought and salinity stress treatment. (A) Number of down-regulated and up-regulated DEGs in samples in response to stress treatment. (B) Venn diagram of down-regulated genes in PUB7-GE (K343-8-1) vs. WT (Dongjin). (C) Venn diagram of up-regulated genes in PUB7-GE (K343-8-1) vs. WT (Dongjin).
Figure 14. Comparative analysis of differentially expressed genes (DEGs) after drought and salinity stress treatment. (A) Number of down-regulated and up-regulated DEGs in samples in response to stress treatment. (B) Venn diagram of down-regulated genes in PUB7-GE (K343-8-1) vs. WT (Dongjin). (C) Venn diagram of up-regulated genes in PUB7-GE (K343-8-1) vs. WT (Dongjin).
Ijms 24 05338 g014
Ijms 24 05338 g015 550
Figure 15. Relative gene expression determined using qRT-PCR. (A) Gene expression of six candidate differentially expressed genes (DEGs) in WT and PUB7-GE (K343-8-1) after drought stress treatment. (B) Expression of six candidate DEGs in WT and PUB7-GE (K343-8-1) after salinity stress treatment. Error bars indicate SEs (n = 3, three biological replicates). Values for WT and PUB7-GE differed significantly. *, p < 0.05; ***, p < 0.001.
Figure 15. Relative gene expression determined using qRT-PCR. (A) Gene expression of six candidate differentially expressed genes (DEGs) in WT and PUB7-GE (K343-8-1) after drought stress treatment. (B) Expression of six candidate DEGs in WT and PUB7-GE (K343-8-1) after salinity stress treatment. Error bars indicate SEs (n = 3, three biological replicates). Values for WT and PUB7-GE differed significantly. *, p < 0.05; ***, p < 0.001.
Ijms 24 05338 g015
Ijms 24 05338 g016 550
Figure 16. Phenotypes of selected PUB7-GE in the T3 generation.
Figure 16. Phenotypes of selected PUB7-GE in the T3 generation.
Ijms 24 05338 g016
Ijms 24 05338 g017 550
Figure 17. Description of scoring conditions for plant phenotypes after abiotic stress treatments. Description of the symptoms by each stress level.
Figure 17. Description of scoring conditions for plant phenotypes after abiotic stress treatments. Description of the symptoms by each stress level.
Ijms 24 05338 g017
Table
Table 1. Analysis of OsPUB7 transformation efficiency and mutation type ratio at the target site in T0 mutant rice plants.
Table 1. Analysis of OsPUB7 transformation efficiency and mutation type ratio at the target site in T0 mutant rice plants.
Target GeneNo. of Plants
Examined		No. of Plants with
Mutations		Mutation Rate
(%)		Putative Homozygous
Mutations		Putative Bi-Allelic
Mutations		Putative Heterozygous
Mutations
No. of Plants%No. of Plants%No. of Plants%
OsPUB7311548426.7640.0533.3
Table
Table 2. Chi-square (χ2) analysis of the genotypes with bar screening from OsPUB7 T1 plants.
Table 2. Chi-square (χ2) analysis of the genotypes with bar screening from OsPUB7 T1 plants.
GeneLineGenotypeNo. of Resistant PlantsNo. of Susceptible PlantsNo. of All Edited Plantsχ2Significance
(α)
OsPUB7K343-7Bi-allelic2675432111.4490.003
K343-8Bi-allelic130421720.0310.07
K343-17Bi-allelic2316395.3420.001
Table
Table 3. Agronomic traits of PUB7-GE in the T3 generation.
Table 3. Agronomic traits of PUB7-GE in the T3 generation.
GeneNo. of LinePlant Height
(cm)		Culm Length
(cm)		Panicle Length
(cm)		No. of Tiller
DongjinWT113.3 ± 1.394.5 ± 1.616.6 ± 0.912.4 ± 1.7
OsPUB7-GEK343-7-5 ns115.1 ± 1.7 *89.7 ± 0.6 *21.6 ± 0.6 *9.7 ± 0.5 ns
K343-8-1 ns111.9 ± 2.3 ns92.6 ± 2.7 ns19.0 ± 1.7 ns15.7 ± 1.0 ns
K343-17-8 *118.5 ± 1.6 ns94.3 ± 1.2 ns18.5 ± 0.9 ns15.0 ± 0.8 ns
*, p < 0.05; ns, not significant.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Kim, M.-S.; Ko, S.-R.; Jung, Y.J.; Kang, K.-K.; Lee, Y.-J.; Cho, Y.-G. Knockout Mutants of OsPUB7 Generated Using CRISPR/Cas9 Revealed Abiotic Stress Tolerance in Rice. Int. J. Mol. Sci. 2023, 24, 5338. https://0-doi-org.brum.beds.ac.uk/10.3390/ijms24065338

AMA Style

Kim M-S, Ko S-R, Jung YJ, Kang K-K, Lee Y-J, Cho Y-G. Knockout Mutants of OsPUB7 Generated Using CRISPR/Cas9 Revealed Abiotic Stress Tolerance in Rice. International Journal of Molecular Sciences. 2023; 24(6):5338. https://0-doi-org.brum.beds.ac.uk/10.3390/ijms24065338

Chicago/Turabian Style

Kim, Me-Sun, Seo-Rin Ko, Yu Jin Jung, Kwon-Kyoo Kang, Yung-Jo Lee, and Yong-Gu Cho. 2023. "Knockout Mutants of OsPUB7 Generated Using CRISPR/Cas9 Revealed Abiotic Stress Tolerance in Rice" International Journal of Molecular Sciences 24, no. 6: 5338. https://0-doi-org.brum.beds.ac.uk/10.3390/ijms24065338

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop