Clone and Function Verification of the OPR gene in Brassica napus Related to Linoleic Acid Synthesis

Fatty acid composition and content affect rapeseed oil quality. Fatty acid synthesis-related genes in rapeseed have been studied globally by researchers. Nevertheless, rapeseed oil is mainly composed of seven different fatty acids (FA), and each fatty acid was regulated by different genes. Furthermore, different FA affect each other, which needs continuous and in-depth research to obtain more clear results in Brassica napus. In this paper, broad-scale miRNA expression profiles were constructed and 21 differentially expressed miRNAs were detected. GO enrichment analysis showed that most up-regulated proteins were involved in transcription factor activity and catalytic activity. KEGG pathway enrichment analysis indicated that 20 pathways involving 36 target genes were enriched, of which the bna00592 pathway may be involved in fatty acid metabolism. The results were verified using a quantitative real-time PCR (RT-qPCR) analysis, we found that the target gene of bna-miR156b > c > g was the OPR (12-oxo-phytodienoic acid reductase). Four copies of OPR gene were found, and the over-expression vectors (pCAMBIA1300-35 s-OPR and pCAMBIA1300-RNAi-OPR) were constructed to verify their functions. In T1 and T2 generation, the content of linoleic acid (LA) increased significantly in OE but deceased in OPRi. This is the first study to provide four copies of the OPR gene that regulates LA metabolism, can be used for the molecular mechanism of LA and optimizing fatty acid profiles in oilseed for breeding programs.


Background
Brassica campestris L. (rapeseed) is one of the most important oil crops in the world [1,2]. The quality of rapeseed oil mainly depends on its fatty acid composition, especially the proportion of three main unsaturated FA: oleic acid (OA, C 18: 1 ), LA (C 18: 2 ), and linolenic acid (C 18: 3 ) [3]. Many studies have shown that rapeseed oil with a high unsaturated fatty acid content will have better health effects and can prevent the occurrence of cardiovascular diseases [4]. LA can prevent or reduce the incidence of cardiovascular diseases and being used for the prevention and treatment of hypertension, hyperlipidemia, angina pectoris, coronary heart disease, atherosclerosis, and senile obesity [5,6]. LA could be hydrogenated to produceα-linolenic acid (ALA) andγlinolenic acid (GLA) [7]. miRNAs are short noncoding regulatory RNAs that regulate gene expression via post-transcriptional repression [8,9]. In plants, miRNAs are involved in various biological processes, including the regulation of plant development [10], architecture formation [11], photosynthesis [9], tolerance to biotic and abiotic stresses [12][13][14][15]. In recent years, several new miRNAs have been found in oilseed rape [16][17][18]. The process of seed development is a period that starts at embryo development and ends when dry seeds are mature; it is a key stage that affects the seed size, oil production, protein content, and antinutritional accumulation of rapeseed. Furthermore, gene expression changes during seed development [19][20][21].
Recently, miRNAs and their target transcripts involved in fatty acid and lipid metabolism have been studied in developing seed of Brassica napus, which is considered the third largest oil crop worldwide [22,23]. The generation of acetyl-CoA and carbon chain desaturase were also observed; nevertheless, the total number of known miRNAs and their functions in Brassica napus are still unknown [23].
The development of molecular biology technology has greatly promoted the research on LA breeding in rapeseed. (I) In terms of mining the function of genes, a study showed significant correlation between FAD3 and LA content [24], while another study revealed that the overexpression of SsDGAT1 significantly affected LA content (about 16%) [25]. McD6DES generates a double bond at the carboxyl end of LA, which reduces the LA content [26]; BJULFY increases LA content (approximately 5%) [27]. The single nucleotide mutation of FAD3 exon from G to A was screened from the rapeseed "Alboglabra", which treated with an Ethyl Methyl Sulfone (EMS) solution, and the mutant materials with high LA content and low alanine content (about 2.0%) were obtained [28]. (II) In terms of screening the genes, a genome wide association study (GWAS) discovered 53 and 24 SNP related to LA in 2013Cq and 2014Cq rapeseed materials, respectively [29]; when combined with RT-qPCR, 95 candidate genes were found to be highly related to LA and other fatty acid metabolism [30]. (III) In terms of locating the genes, FAD2 gene on chromosome A05 was found to be associated with LA content, and minor gene regulating LA content was found on chromosome A09 [31]; twenty QTLs related to LA were detected using two spring rapeseed varieties, "Polo" and "Topas", and subsequently distributed in seven linkage groups: A01, A02, A03, A05, C01, C03, and C09 [32]. In fact, many studies indicate that OA metabolism depend on FAD1 or FAD2 activity [33,34], and the have present the first evidence of Δ7 desaturation via the FADS1 gene product by Hui et al. [35] LA are used as precursors of ultra-long chain polyunsaturated FA, such as α-linolenic acid (ALA) and γ-linolenic acid (GLA), which are hydrogenated by desaturation enzymes [36]. Correlations between oleic versus linoleic and oleic versus linolenic were negative and highly significant, and correlations between linoleic versus linolenic were of lower magnitude by Kondra et al. [37]. According to Belo et al. [38], the impact of presenting an interesting study focusing on environmental effects on oil quality of oil crops, who showed that the differential effects of climate conditions had an effect on fatty acid composition. miRNA technology has been adopted to find several miRNAs in enzymes used for carbon chain desaturation, while being used in studies on fatty acid and lipid metabolism in rapeseed [22]. Employing this technology can be conducive to the breeding of rapeseed [23]. Therefore, we sequenced rapeseed material using miRNA technology. Then, we excavated four copies of OPR genes. The four copies of gene overexpression and RNAi vectors were transferred into A. thaliana. This is the first study to report that all copies of OPR gene could directly regulate the synthesis of LA. These findings will enhance the understanding of LA metabolism in rapeseed.

Overview of sRNA sequencing results
We have already upload data and ensure the deposited data is made public (

Differentially expressed miRNAs identification and their function analysis
The clustering analysis method was used to investigate the similarity between samples by calculating the differential miRNA distance between LOAR (A) and HOAR (B) (Fig. 1B). As shown, three samples of high or low OA materials were found in a cluster, suggesting that these miRNAs might have similar biological functions. After applying a P-value < 0.05 and an absolute value of log2 (treatment/control: LOAR/HOAR) greater than 1.5 to identify differentially expressed miRNAs, 21 differentially expressed miRNAs (8.39% of the total) were detected ( Table 2, Supplemental  Table 2); their frequencies were calculated using TPM (Table 3). Among them, nine genes (42.86%) were upregulated and 12 (57.14%) were down-regulated in HOAR materials (Fig. 1C). Briefly, 21 miRNAs (including .1_17408*, NC_027769.1_17408, NC_027768.1_15573, NC_027761.1_6665, NC_027760.1_5272*, NC_027760.1_5272, NC_027757.1_266) and 358 target genes were obtained. Target finder software was used to analyze the target genes of differential miRNAs, refer to Fahlgren and Carrington (2010) [62]. In a word, 21 miRNAs were found to target 358 genes. Subsequently, according to p-value < 0.05, 133 putative target genes were obtained by screening significantly enriched GO and KEGG pathways from 358 target genes.
The target genes were then subjected to GO functional and KEGG Pathway analyses. In many cases, multiple terms were assigned to the same miRNA. Thus, 133 putative target genes were associated with 21 differentially expressed miRNAs and distributed into the following subcategories: 62 "Biological process", 32 "Cellular component", and 39 "Molecular function" (Fig. 1D). Under "Biological process", most of the target genes were related to "transcription", and "regulation of transcription". Within the "Cellular component" category, "nucleus" and "cytosol" were observed as much as "cytoplasm". Among genes in the "Molecular function" category, most potential functions were related to "transcription factor activity", "DNA binding", and "Catalytic activity". The distribution of target genes indicated that rapeseed underwent active metabolization.

Expression pattern of bna-miR156b > c > g gene was detected by RT-qPCR
To confirm the results of the miRNA sequence analysis, 21 annotated differentially expressed miRNAs were compared to the Brassica napus genome using BLAST [39] (Fig. 2B).
Most expression trends of the RT-qPCR analysis results agreed with the miRNA sequencing data (NC_027757.1_266, NC_027760.1_5272, NC_027760.1_5272*, NC_027761.1_6665, NC_027768.1_15573, NC_027769.1_17408, NC_027769.1_17408*, NC_027774.1_24533, NC_027774.1_24533*, NW_013650328.1_26640*, bna-miR162a, bna-miR167a > b, bna-miR169m, bna-miR172d, bna-miR396a, bna-miR824.), In addition, among the 13 miRNAs with significant difference, 9 had which may be the novel miRNAs related to FA. Moreover, the expressions of fatty acid metabolism related to differential miRNAs, such as bna-miR396a, bna-miR156b > c > g, and their target genes, were studied in different developmental stages (Fig. 2C). The bna-miR396a has opposite expression pattern with its target gene, at first, the bna-miR396 had up-regulated expression, until the bud stage reached the peak, and the expression decreased with the growth stage; bna-miR156b > c > g had an opposite expression pattern with its target gene, it had down-regulated expression in the whole growth stages, and the expression decreased with the growth stages. Differentially expressed miRNAs and their target genes were related to fatty acid metabolism in bna-mi156b > c > g at different developmental stages (Fig. 1C). In contrast to the expression pattern of bna-miR156b > c > g and its target gene, the expression level was down-regulated throughout the whole growth stage of the plant, and the expression level gradually decreased with the growth process of rapeseed.

Cloning of OPR genes in rapeseed and bioinformatic analysis
Target gene: bna-miR156b > c > g was cloned by miRNA sequencing, and four copies were detected: GSBRNA2T00012422001, GSBRNA2T00135385001, GSBRNA2T00082938001, and GSBR-NA2T00094910001, which were named OPR1, OPR2, OPR3, and OPR4, respectively. Both OPR1 and OPR3 were 1119 bp, OPR2 and OPR4 were 1125 bp and 1122 bp, respectively (Fig. 3A).  (Fig. 3B), homology was more than 99% with the published sequences (Fig. 3C). OPR1, OPR3, OPR4 had no base difference with the published sequence, there were 10 base differences between OPR2 and published sequence and the homology was 99.11%. Preliminary identification of OPR2 and OPR3 were located in A genome and OPR4 and OPR1 in C genome was conducted.
The number of four copies of OPR gene amino acids ranged from 372 to 374 bp, with a molecular weight of about 41 ku; The encoded amino acids were acidic (< 7), unstable (< 40), exhibited a fat coefficient of about 75, and belonged to fat-soluble proteins. Predicted subcellular localization of proteins were encoded by different copies of OPR genes and we found that these four proteins were located in the cytoplasm. We found that the four copies were all extracellular proteins without a transmembrane structure. Predicted protein secondary structures were summarized in Table 4.
The predicted tertiary structure model of the protein shows that the tertiary structures of OPR1, OPR2, OPR3, and OPR4 were all adapted to the 12-O-plant dienoate reductase model (integrated with the crystal structure of At1g76680 protein in A. thaliana), but their  conformations were slightly different (Fig. 3D). It has been reported that the protein structure A. thaliana of At1g76680 is similar to that of yeast ScOYE1 [40].
The recombinant vector was transformed into A. thaliana and the obtained A. thaliana was detected. The T 1 transgenic A. thaliana seeds were screened using hygromycin (Fig. 5A) and the results of hygromycin primer identification (Fig. 5B) showed that each copy of A. thaliana OPR had been successfully transformed, and the target plasmid T-DNA had been inserted into the genome of A. thaliana.

Analysis of fatty acid composition and fatty acid content
The transformation methods in A. thaliana, reference Clough and Bent (1998) [41]. fatty acid composition was detected by gas chromatography [42]. we obtained fatty acid composition results of the T 1 and T 2 generations. Fatty acid composition of A. thaliana T-DNA insertion lines in Table 5, and transformation materials fatty acid composition in the contrast T 1 /T 2 . The contents of OA and stearic acid in OPR1i were significantly increased and the LA content decreased significantly; OPR1-OE will lead to a significant increase in palmitic acid and LA content; OPR2i significantly increased stearic acid content and decreased LA content; and OPR2-OE increased LA content. The RT-qPCR results of 21 differential miRNAs. The results of low oleic acid rapeseed materials were used as control. (* and ** signify the difference level of P < 0.05 and P < 0.01, respectively); (C) The expressions of miRNAs (miR396a, miR1566b > c > g) and their target genes. The expression of miRNA and its target gene were used as controls The LA content in OPR3i decreased significantly; OPR3-OE significantly increased LA content; OPR4i significantly decreased LA content; and OPR4-OE increased LA content significantly.
Each copy of OPR-OE increased LA content, with an average increase of 12.56% in T 1 generation and 7.185% in T 2 generation. Subsequently, LA content in OPRi gene was significantly decreased, with an average decrease of 5.98% in T 1 generation and 0.86% in T 2 generation.
As shown in Table 5, oleic, linolenic, arachidonic, and erucic acids with the same variation trend as that of the fatty acid composition were selected for variance analysis. The results (Table 6) showed that the linolenic acid content in OPR1i significantly increased, while OPR2i significantly increased the linolenic acid content. Both OPR3-OE and OPR4-OE affected the content of arachidonic acid, which decreased significantly. In addition, OPR4i had no significant effect on the arachidonic acid content. We correlation analysis of FA in transgenic Arabidopsis seeds, according to the Table 7 shows that: In T 1 , the four OPR copies showed significant negative correlation between the OA, LA and linolenic acid content, OA mass fraction is higher, LA and linolenic acid mass fraction, the lower the relative; There was a positive correlation between LA and linolenic acid. Overall variation trend of LA and linolenic acid was basically the same. The mass fraction of linolenic acid was positively correlated with that of arachidonic acid. The stearic acid massed fraction and LA mass fraction had the same trend, and a larger trend with that of OA and linolenic acid. The mass fraction of palmitic acid was positively correlated with stearic acid, negatively correlated with OA (except OPR1), and positively correlated with LA and linolenic acid. The variation trend of FA in T 1 and T 2 generations is almost the same, there are also differences, and here was a significant negative correlation between LA and linolenic acid.

miRNA expression and enrichment analysis
In recent years, the number of known miRNAs has increased continually, such as A. thaliana [43] and Oryza sativa [44,45]. Brassica napus L. has a relatively high genome size and complexity [39], and the number and function of miRNAs in Brassica napus L. have not been adequately studied; which suggests that many miRNAs have not yet been discovered, especially in seeds.
In this study, miRNA libraries were constructed from self-pollinated seeds that were collected 20-35 d after pollination of rapeseed with high OA content in the near-isogenic lines. Clean reads exhibited 87.45% (A) and 88.26% (B) homology with the referenced genome sequences ( Table 1). The small RNAs (24 nt) were most abundant in all the samples (Fig. 1A). The clustering analysis results showed that three samples of high or low OA contents were found in a cluster, revealing that these miRNAs might have similar biological functions (Fig. 1B). A total of 21 differentially expressed miRNAs were detected ( Table 2, Supplemental Table 2), including 9 (42.86%) up-regulated and 12 (57.14%) down-regulated genes (Fig. 1C) using GO and KEGG pathways. In our study, bna-miR156b > c > g may be involved in fatty acid metabolism. Genes related to Cd stress have previously been discovered using GO and KEGG pathways, such as BNPCS1, BNGSTU12, and BNGSTU5 [46]. Jian et al. (2018) found a total of 13 differentially expressed miR-NAs were confirmed by RT-qPCR, and a hypothetical model of cadmium response mechanism in B.napus was proposed on this basis [47].

Multi copy phenomenon of OPR genes
We found that OPR gene has four copies in rapeseed, with OPR1 and OPR4 located in the C genome and OPR2 and OPR3 located in the A genome. The OPR gene has been identified in several species and there are often multiple copies. Three OPRs were found in A. thaliana and tomato [48], 6 in peas [49], and 8 in corn [50]. Meanwhile, rice comprised 13 OPRs [51] and wheat had 48 OPRs [52]. Multi-copy genes are ubiquitous in plants and play an important role in maintaining plant genetic stability; However, they have hindered molecular breeding research. The loss of a few copies of gene function often does not cause phenotypic changes, and the probability of simultaneous mutation of multiple copy sites is too low to create a gene family or multiple copies of genes change simultaneously [53]. Rapeseed is an allotetraploid with many multi-copy genes. Conventional molecular breeding research methods are difficult to obtain phenotypic multi-copy gene mutants. Handa (2003) found that the main DNA sequence of the protein coding region was highly conserved between rapeseed and A. thaliana [54]. Transformation with A. thaliana as a receptor is helpful to study the function of multi-copy genes in rapeseed. In this study, we found four copies of OPR gene in rapeseed.
Based on these multi-copy genes, we transformed four copies of OPR genes in A. thaliana.

Regulation function of OPR genes to fatty acid
In this study, we found that OPR genes may affect the metabolism of LA and each copy was transferred separately into A. thaliana. The LA content of OPR-OE transgenic plants was significantly increased (T 1 12.56%, T 2 7.185%), while the LA content of OPRi transgenic plants was significantly decreased (T 1 5.98%, T 2 0.86%). These results have rarely been reported before. However, it has been described that the ClOPR genes, particularly ClOPR2 and ClOPR4, significantly upregulated by exogenous jasmonic acid, salicylic acid, and ethylene treatments in watermelon [55]. Virus-induced gene silencing (VIGS) analysis suggested that knockdown of GhOPR9 could increase the susceptibility of cotton to V. dahliae infection [56]. OPR gene was cloned from Oryza sativa; the overexpression of OPR genes was found to enhance the stress resistance of tobacco to heavy metal Cd 2+ [57]. Previous studies concluded that OPR genes were widely involved in abiotic stress processes [58]. In addition, OPR genes involved in fatty acid β oxidation, cilinolenic acid reduction, and the octadecanoic acid metabolic pathway [59]. However, there were few reports on OPR genes regulating OA, LA, or saturated fatty acid synthesis. In our study, we found that OPR genes directly affect the synthesis of LA and indirectly affect the content of other FA (Fig. 6), which is consistent with the theoretical pathway wherein OPR genes regulate jasmonic acid synthesis using alpha-linolenic acid (18: 3) [59,60]. According to the Table 7, The four OPR copies showed significant negative correlation between the OA, LA and linolenic acid content, which was consistent with Zhao J (2008) [61] and Yan (2012) [62]; The stearic acid massed fraction and LA mass fraction had the same trend, and a larger trend with that of OA and linolenic acid. This is consistent with findings of Shang et al. [63]. Linolenic acid content is not only affected by OA and LA, but also affected by another independent gene system, which was consistent with Kondra (1975) [64]. The variation trend of FA in T 1 and T 2 generations is almost the same, there are also differences, and here was a significant negative correlation between LA and linolenic acid. In the meantime, there are many researchers believed that FA content of rapeseed was controlled by maternal genotype, and there were interaction, additive and nonadditive gene effects between genotype and environment [65][66][67].
A. thaliana and rapeseed are both cruciferous plants and current studies have shown that their gene functions are basically the same, it is of practical significance to study in model plants [54]. We found the optimizing quality FA in A. thaliana can be changed by regulating OPR genes. OPR may be involved in regulating LA synthesis and improving fatty acid composition in rapeseed. This is the first study which discovered that the OPR gene can regulate LA metabolism. Therefore, this study is a good reference for studies researching the molecular mechanism of LA synthesis and molecular breeding in rapeseed.
In this study, 20 pathways were enriched using the KEGG pathway through high-throughput sequencing, of which 15 may be involved in the regulation of fatty acid metabolism. The reliability of the results was verified by performing the RT-qPCR analysis, which provided a basis for subsequent functional verification. We excavated a target gene OPR, bna-miR156b > c > g, from rapeseed that may be related to fatty acid synthesis and identified the function of OPR genes through transformation of A. thaliana. The LA content of OPR-OE transgenic plants significantly increased (T 1 12.56%, T 2 7.185%), and the LA content of OPRi transgenic plants significantly decreased  Table 5 Analysis of fatty acid composition of T1 and T2 seeds Note:*. P < 0.05; **. 0.01 < P < 0.05. The same as (T 1 5.98%, T 2 0.86%). In addition, by performing a bioinformatics analysis, we found four copies of the OPR gene in the cytoplasm that were located on chromosomes A and C. In this study, by detecting the fatty acid content of different generations of transgenic A. thaliana. OA and LA, linolenic acid content showed significant relationship, LA can be affected by FAD-related gene regulation and the environment. We found the four copies of OPR gene that can directly affect LA content and indirectly affect other high quality FA were discovered for the first time. These results can be used in breeding programs aimed at optimizing fatty acid profiles in rapeseed [68].  Table 7 FA correlation analysis of T1 and T2 seeds Note: Lower triangle is r,upper triangle is p;"**" highly significant correlation(P < 0.01),"*" significant correlation(P < 0.05)

Plant materials and growth conditions
The near-isogenic rapeseed lines with high (81.4%) and low (56.2%) OA contents were provided by Rapeseed Molecular Breeding of Hunan Agricultural University, high OA rapeseed materials became known as HOAR and low OA rapeseed materials became known as LOAR, the strain was originally cultivated by National Oil Improvement Center of Hunan Agricultural University. The materials exhibited stable traits and were planted in the  [69]. The seeds 20 d after pollination were quickly frozen in liquid nitrogen and stored at -80 °C for subsequent studies [70]. The sample treatment method adopted in this study is the same with our previous study [71].
The seeds of wild-type (WT) A. thaliana (ecotype: Columbia) were bought from Think Gene Biological Technology Co., LTD, Shanghai, China. Plants were grown under greenhouse conditions: 24 °C with a photoperiod of 18.5 h-light/5.5 h-dark, with a light intensity of 6500 lx.

sRNA library construction and high throughput sequencing
Total RNA was extracted from frozen seeds using the Trizol reagent (Sigma Aldrich, St. Louis, MO, USA), according to the manufacturer's instructions. The quality and quantity of the purified RNA were assessed using an Agilent 2100 Bioanalyzer (Agilent Technologies, Palo Alto, Santa Clara, USA) and RNA 6000 nanokit (Agilent Technologies, Palo Alto, Santa Clara, USA). sRNAs with lengths of 18-30 nt were separated and purified using 15% denaturing polyacrylamide gel electrophoresis. Consequently, sRNAs fractions were ligated to the 5′ lectrophoresters using T4 RNA ligase (Epicentre, America). The adapter-ligated fragments were then reverse transcribed and amplified by performing PCR with a pair of adapter complementary primers. These PCR products were purified and sequenced using IlluminaHiseq XTEN (Illumina, USA). Construction of the sRNA libraries and deep sequencing were carried out by Oebiotech Genomics (Shanghai, China).

Bioinformatic analysis of the sRNAs sequencing data
Clean reads were generated after eliminating the lowquality reads, poly As, reads smaller than 18 nt, and gener adaptor contaminants and subsequently inserting nulls. Using bowtie software [72], the clean reads were aligned against the NCBI Gen Bank [73], Brassica napus oilseed genome [39], and Rfam databases (version 10.0). Reads annotated into the noncoding RNA categories, including rRNA, tRNA, snRNA, and snoRNA were filtered. The remaining sRNA sequences were aligned against the mi RBase21.0 [74]. The nearly matched sequences (less than two mismatches) were considered to be the known miRNAs. Target miRAN will computationally predict small RNA binding sites on target transcripts from a sequence database, this is done by aligning the input small RNA sequence against all transcripts, followed by site scoring using a position-weighted scoring matrix [75]. First, used 'targetfinder' software, set threshold -C < = 4, predicting the miRNA target genes, we got 64 target genes with miRNA156; and then, to GO and KEGG analysis (P-value ≤ 0.05) of 64 target genes, we found there are six of them related to the metabolism of FA. The unannotated sRNAs were further analyzed to predict the novel miRNAs using the Mirdeep 2 software. The secondary structures of premiRNAs were also predicted using the RNAfold software.
To reveal the continuous changes in expression of miRNAs during the biosynthesis process, the variation in expression was analyzed in immature seed libraries of HOAR. The frequency of miRNAs was normalized as transcripts per million (TPM) for further analysis. In addition, miRNAs were assessed using the negative binomial distribution test, with P-value ≤ 0.05 and absolute value of log2 (treatment/control: LOAR/HOAR) > 1.5 being considered as differentially expressed. Moreover, the similarity between samples was investigated by the clustering method. The Blast 2 GO software with default parameters was applied to determine the functional annotation and categorization of the target genes [76]. The KEGG (Kyoto Encyclopedia of Genes and Genomes) and InterPro databases were also searched with an evalue of 1e −10 [77][78][79].

Real time quantitative polymerase chain reaction (RT-qPCR) validation of differentially expressed miRNAs
To validate the expression of differentially expressed miRNAs, all annotated miRNAs were selected for RT-qPCR validation using the poly (T) adaptor RT-qPCR method [80]. The RT-qPCR amplifications were performed according to an established procedure [81]. 5S rRNA was used as an internal reference gene of miRNA and two primers were used: F: 5'-CTC GGC AAC GGA TAT CTC G-3' and R: 5'-CTA ATG GCT TGG GGCG-3' . The internal reference gene of the miRNA target gene was UBC9 and two primers were used: F: 5'-TCC ATC CGA CAG CCC TTA CTCT-3' and R: 5'-ACA CTT TGG TCC TAA AAG CCACC-3' . All reactions were performed in triplicate. The RT-qPCR was conducted on the Step One Plus RT-qPCR System (ABI, America). The relative expression was calculated using the 2 − ΔΔCT method [82]. Statistical analyses consisted of analysis of variance and Fisher'spost-hoc tests. A p-value < 0.05 was considered statistically significant. Details of the RT-qPCR primers were provided in Table 2.

Cloning of OPR genes in rapeseed and bioinformatic analysis
Four copies of OPR genes were found using miRNA sequencing homologous cloning. We used material HOAR (high OA rapeseed materials materials, B) to clone, and the total RNA was extracted from the rapeseed leaf samples using the CTAB method adopted by Niu et al. (2018) and Meng et al. (2006) [83,84], with minor modifications. cDNA was produced by implementing reverse transcription, which was followed by PCR amplification (Table 8). According to each copy sequence of OPR, the physical and chemical properties of corresponding proteins were analyzed by using online websites such as the ExPASY-Protparam Tool and modifying the methods described by Li et al. (2019) and Sun et al. (2016) [85,86].

Vector construction and transformation of A. thaliana
Specific RNAi primer (Table 9) cDNA was used as a template. The base before the 5'-end enzyme digestion site was used for the construction of recombinant plasmids (the same below), according to the method adopted by Li et al. (2021) with minor modifications [87]. The cloned OPR sequence was digested with BamH I and Sal I ( Table 8) and recombinant with pCAMBIA1300-35 s vector (Kangyan Corporation, China) to construct pCAM-BIA1300-35 s-OPR vector.
Similarly, the pCAMBIA1300-RNAi-OPR vector (Kangyan Corporation, China) was reconstructed using a double digestion process in vitro (Table 9) [6,29]. pCAM-BIA1300-35 s-OPR and pCAMBIA1300-RNAi-OPR vectors were transformed into A. thaliana because of inflorescence infection to explore the function of each copy by using the inflorescence infection method of Clough and Bent (1998) [41], with minor modifications.

Fatty acid detection and statistical methods
Using Agilent Technologies 7890B gas chromatograph, MSD detector, 5977A chromatographic workstation, and DB-23 capillary chromatographic column, sample size Table 8 Primers for OPR genes amplification Note: the underline indicates the restriction site. the same as Table 9 Target ID Primer name Primer sequence (