Construction and Functional Evaluation of CRISPR/Cas9 Multiple Knockout Vectors of the FAD2 Gene Family

: Soybean oil is a traditional edible oil. Increasing the oleic acid content is an important direction of soybean breeding. The soybean FAD2 family consists of seven genes that regulate how oleic acid is converted into linoleic acid. Five genes of the soybean FAD2 gene family, GmFAD2-1B , GmFAD2-1A , GmFAD2-2B , GmFAD2-2C , and GmFAD2-2D , were taken as target genes in this study. Firstly, a multivalent CRISPR/Cas9 gene-editing vector was constructed to regulate FAD2 gene expression. Multiple knockout vectors were inserted into the soybean varieties JN38, T6098, and T7010 using Agrobacterium -mediated soybean cotyledon transformation. The functional analysis, agronomic character analysis, and comparison of the mutant lines of the offspring of different genotypes indicated that the JN38 mutant was signiﬁcantly taller in terms of plant height than the receptor JN38. The fatty acid content of the three groups showed the same trend. The fatty acid contents of mutant plants were higher than those of recipient plants, and the linoleic acid contents of mutant plants were lower than those of recipient plants. The best-performing among the three groups was the JN38 mutant, whose oleic acid content increased from 18.58% to 54.07% and whose linoleic acid content decreased from 57.79% to 26.17%. In conclusion, the knockout expression of multiple FAD2 genes increased the soybean oleic acid content and decreased the linoleic acid content in different receptors.


Introduction
Soybean (Glycine max (Linn.)Merr.) originated from China and is one of the most important legume and oil crops in the world [1].The index evaluation of soybean oil is mainly related to the content and composition ratios of five fatty acids: oleic acid (C18:1), linoleic acid (C18:2), linolenic acid (C18:3), stearic acid (C18:0), and palmitic acid (C16:0) [2].Linoleic acid accounts for the highest proportion of the total fatty acids, followed by oleic acid [3].Unsaturated fatty acids can not only prevent the incidence of cardiovascular diseases and reduce the human blood cholesterol content [4,5], but their proportion also affects the quality of edible oil and the nutritional value [6].
Linoleic acid is usually partly hydrogenated during manufacturing to improve its oxidative stability and shelf life [7].However, this process can easily produce trans-fatty acids, which pose a great danger to people's health [8].Oleic acid, a monounsaturated fatty acid containing only one unsaturated double bond, is more stable and has greater antioxidant capacity.The oleic acid/linoleic acid (O/L) ratio has a great influence on soybean food production, industrial storage, etc.Therefore, increasing the relative content of oleic acid and the ratio of oleic acid/linoleic acid (O/L) has become an important goal in soybean breeding [9].The Fatty Acid Desaturase 2 (FAD2) gene family plays a key role in regulating the conversion of oleic acid to linoleic acid [10,11].When plants synthesize linoleic acid and linolenic acid, the endoplasmic reticulum w-6 fatty acid desaturase (FAD2 enzyme) activity influences the synthesis of unsaturated fatty acids.In the metabolic process of fatty acids, FAD2 plays a significant role [12].Studies have shown that the FAD2 gene family can control oil storage [13], improve the oleic acid content, and reduce the linoleic acid content.For example, monounsaturated fatty acid (MUFA) levels were increased by 80% in the seeds of plants that had all three pairs of FAD2 homologues knocked out of hexaploid Camelina sativa.Additionally, when two of the converters of the CsFAD2 homologous gene were knocked out, the seed monounsaturated fatty acid (MUFA) content increased by 60% [14].Zhang et al. introduced the antisense FAD2-1 gene into soybean cotyledon nodes, and the results showed that the oleic acid content of the transsense FAD2-1 gene seed increased to 51.71% [15].Two genes in the soybean FAD2 gene family, GmFAD2-1B and GmFAD2-2C, were constructed with a single gene knockout vector and a double gene knockout vector for genetic transformation.The results showed that the regulation effect of multi-gene editing on the formation of linoleic acid was stronger than that of a single gene-editing body [16].
In recent years, CRISPR/Cas9 technology has been widely used in crop improvement.Bao et al. used CRISPR/Cas9 to mediate the mutation of four different combinations of GmSPL9 genes to change the plant structure of soybean [17].Li et al. used CRISPR/Cas9 to mutate the soybean seed storage protein gene, which provided germplasm resources for developing high-protein varieties [18].CRISPR/Cas9 gene-editing technology has been used to construct knock-out vectors of different genes of the FAD2 gene family in soybean for functional verification [16,19,20], but the effect of editing five genes of the soybean FAD2 family with knock-out vectors on oleic acid accumulation has not been reported.
Five genes with high homology were screened from the soybean FAD2 gene family.A multi-gene knockout expression vector was constructed using CRISPR/Cas9 gene editing technology, and five target genes of the FAD2 family were knocked out.The experimental receptor materials were screened from the natural population materials constructed in the previous period.The oleic acid contents of JN38, Z6098, and Z7010 were observed to be in the low range of 17-19% after many years of field planting in multiple sites.The higher linoleic acid contents ranged from 56% to 58%.During the growth cycle, all plants showed high adaptability, early maturation, and good lodging resistance.Therefore, we chose these three materials for use in the experiment.Three different soybean receptors, JN38, T6098, and T7010, were subsequently transformed by Agrobacterium-mediated genetic transformation.Mutants were screened, and the quality and agronomic traits of the transformed offspring of different genotypes were analyzed and compared to further elucidate the function of FAD2 gene-family genes.The results of this study show that when the CRISPR/Cas9 technique was used to introduce knockout vectors into the receptors JN38, T6098, and T7010, this clearly improved the grain qualities of the T3 gene-edited plants.In terms of agronomic traits, the plant height of the JN38 transgenic polyvalent vector was significantly higher than that of the receptor lines.The results of this experiment provide a germplasm basis for the molecular breeding of high-oleic-acid soybean.

Construction and Validation of CRISPR/Cas9 Multiplex Knockout Expression Vector
For comparison, the constructed vectors were sequenced.The results show that the target site sequences in the CRISPR/Cas9 multiple knockout vectors were 100% consistent with the target site sequences in the reference sequence, indicating that the multiple knockout vectors were successfully constructed (Table 1).The findings in respect of the vector structure (Figure 1) and sequencing comparison are displayed in Figure 2 below.

Gene
Target Sequence + PAM SG Sequence

GmFAD2-2C
TGGCAGGAGGAACAGCAGTT CGG GGGCAGGAGGAACAGCAGT PEER REVIEW 3 of 14 Table 1.Target design information, the yellow part is the PAM structure.

Gene Target Sequence + PAM SG Sequence
TGGCAGGAGGAACAGCAGTTCGG GGGCAGGAGGAACAGCAGT

PEER REVIEW 3 of 14
Table 1.Target design information, the yellow part is the PAM structure.

Gene Target Sequence + PAM SG Sequence
TGGCAGGAGGAACAGCAGTTCGG GGGCAGGAGGAACAGCAGT

Detection of Positive Soybean
The T0 generation produced transgenic offspring of soybean.The findings of several of the assays are displayed in Figure 4 below.

Detection of T1 Generation Mutant Plants of CRISPR/Cas9 Multiple Knockout Vector
The seeds harvested from the T0 generation were bred and planted in an artificial chamber.Seeds with full and well-developed kernels were sown in soil and cultivated in an artificial culture room (25 °C, 16 h light/8 h dark).Seven plants were randomly selected from each of the three different recipient plants from generation T1, and the target segments of a total of 21 positive soybean plants were sequenced randomly.The editing sites of the CRISPR/Cas 9 system were detected and compared with the sequences of the unmodified plants.

Detection of Positive Soybean
The T0 generation produced transgenic offspring of soybean.The findings of several of the assays are displayed in Figure 4 below.

Detection of Positive Soybean
The T0 generation produced transgenic offspring of soybean.The findings of severa of the assays are displayed in Figure 4 below.

Detection of T1 Generation Mutant Plants of CRISPR/Cas9 Multiple Knockout Vector
The seeds harvested from the T0 generation were bred and planted in an artificia chamber.Seeds with full and well-developed kernels were sown in soil and cultivated i an artificial culture room (25 °C, 16 h light/8 h dark).Seven plants were randomly selecte from each of the three different recipient plants from generation T1, and the target seg ments of a total of 21 positive soybean plants were sequenced randomly.The editing site of the CRISPR/Cas 9 system were detected and compared with the sequences of the un modified plants.

Detection of T1 Generation Mutant Plants of CRISPR/Cas9 Multiple Knockout Vector
The seeds harvested from the T0 generation were bred and planted in an artificial chamber.Seeds with full and well-developed kernels were sown in soil and cultivated in an artificial culture room (25 • C, 16 h light/8 h dark).Seven plants were randomly selected from each of the three different recipient plants from generation T1, and the target segments of a total of 21 positive soybean plants were sequenced randomly.The editing sites of the CRISPR/Cas 9 system were detected and compared with the sequences of the unmodified plants.
A total of 7 mutant plants were obtained, and 5 FAD2 gene mutations were found in the offspring plants of transgenic soybean, with an average mutation rate of 52.4%.Seven mutants were edited at all four targets and at different sites (Table 2).Among them, there were three mutants of the JN38 variety (Figure 5).The main types of editing were basic insertions, deletions, and substitutions.The results show that the target gene was successfully knocked out.
Table 2. Detection of target mutation in CRISPR/Cas9 polyvalent vector-positive soybean, the yellow part is the PAM structure.

Plant Genotype
Target Spot Target Sequence Mutation Type

PCR Detection of Mutant Plants
DNA was extracted from the T0 (Figure 6) and T1 (Figure 7) generation plants.The Bar and Cas9 genes in the DNA of the modified plants were detected using PCR, and specific bands in accordance with the length of the target segments were visible (marker bar band: 552 bp, marker Cas9 band: 663 bp).Neither blank nor negative water controls showed bands.This indicates that the editing vector was successfully transferred into positive soybean plants.

PCR Detection of Mutant Plants
DNA was extracted from the T0 (Figure 6) and T1 (Figure 7) generation plants.The Bar and Cas9 genes in the DNA of the modified plants were detected using PCR, and specific bands in accordance with the length of the target segments were visible (marker bar band: 552 bp, marker Cas9 band: 663 bp).Neither blank nor negative water controls showed bands.This indicates that the editing vector was successfully transferred into positive soybean plants.

PCR Detection of Mutant Plants
DNA was extracted from the T0 (Figure 6) and T1 (Figure 7) generation plants.The Bar and Cas9 genes in the DNA of the modified plants were detected using PCR, and specific bands in accordance with the length of the target segments were visible (marker bar band: 552 bp, marker Cas9 band: 663 bp).Neither blank nor negative water controls showed bands.This indicates that the editing vector was successfully transferred into positive soybean plants.

Southern Blot Analysis of T2 Transgenic Plants
In order to determine whether the foreign gene was integrated into the genome of T2 generation plants, we used Bar as a probe to perform Southern blot detection on some of the T2 generation edited plants.The results show that the untransformed negative plants had no hybridization signal, while the edited plants showed that the exogenous gene was integrated into the genome in the form of a single copy, and the hybridization signals between individuals were slightly different (Figure 8).

Southern Blot Analysis of T2 Transgenic Plants
In order to determine whether the foreign gene was integrated into the genome of T2 generation plants, we used Bar as a probe to perform Southern blot detection on some of the T2 generation edited plants.The results show that the untransformed negative plants had no hybridization signal, while the edited plants showed that the exogenous gene was integrated into the genome in the form of a single copy, and the hybridization signals between individuals were slightly different (Figure 8).

RT-qPCR Assay of Gene-Edited Progeny
The β-lactin gene is an internal reference gene.The relative expression analysis of the gene shows that in the seeds of T2 generation genome-edited vector transgenic lines, the relative expression level of five target genes was significantly lower than that of the control line.In addition, the relative expression level of five target genes in the leaves of T2 generation genome-edited vector transgenic lines was also significantly lower than that of the control lines (Figure 9).This indicates that the relative expression levels of the five target genes were affected in both leaves and seeds.

RT-qPCR Assay of Gene-Edited Progeny
The β-lactin gene is an internal reference gene.The relative expression analysis of the gene shows that in the seeds of T2 generation genome-edited vector transgenic lines, the relative expression level of five target genes was significantly lower than that of the control line.In addition, the relative expression level of five target genes in the leaves of T2 generation genome-edited vector transgenic lines was also significantly lower than that of the control lines (Figure 9).This indicates that the relative expression levels of the five target genes were affected in both leaves and seeds.
gene shows that in the seeds of T2 generation genome-edited vector transgenic lines, the relative expression level of five target genes was significantly lower than that of the control line.In addition, the relative expression level of five target genes in the leaves of T2 generation genome-edited vector transgenic lines was also significantly lower than that of the control lines (Figure 9).This indicates that the relative expression levels of the five target genes were affected in both leaves and seeds.

Determination of Fatty Acid Content of Different Editing Vectors in T2 Plants
All fatty acid components of T2 positive soybean were changed compared with their respective receptors, with a significant increase in oleic acid content and a significant decrease in linoleic acid content (Figure 10).In contrast, palmitic acid content decreased as oleic acid levels increased, followed by an increase in stearic acid content (Table 3).The effect of different receptors on oleic acid content was different.Plants with simultaneous knockout of two or three FAD2 target genes of the same receptor were less able to regulate oleic acid accumulation than plants with knockout of five target genes.The high oleic acid material mutant with the best fatty acid trait screened from the T2 transgenic line was JN38-FAD2-1A-1B-2B-2C-2D-mutant-1.Its oleic acid content in soybean seeds was 54.07%, while its linoleic acid content was reduced to 26.17% and its linolenic acid content was reduced to 7.67%.The findings suggest that the combined action of five target genes in the oleic acid metabolic pathway prevented the metabolism of oleic acid to linoleic acid and linolenic acid, which led to an increase in oleic acid content and a decrease in linoleic acid and linolenic acid contents.

R PEER REVIEW
8 of 14 levels in leaf; (B) the relative expression levels in seed.** indicate significant differences in the studied traits of the gene expression at p ≤ 0.01.

Determination of Fatty Acid Content of Different Editing Vectors in T2 Plants
All fatty acid components of T2 positive soybean were changed compared with their respective receptors, with a significant increase in oleic acid content and a significant decrease in linoleic acid content (Figure 10).In contrast, palmitic acid content decreased as oleic acid levels increased, followed by an increase in stearic acid content (Table 3).The effect of different receptors on oleic acid content was different.Plants with simultaneous knockout of two or three FAD2 target genes of the same receptor were less able to regulate oleic acid accumulation than plants with knockout of five target genes.The high oleic acid material mutant with the best fatty acid trait screened from the T2 transgenic line was JN38-FAD2-1A-1B-2B-2C-2D-mutant-1.Its oleic acid content in soybean seeds was 54.07%, while its linoleic acid content was reduced to 26.17% and its linolenic acid content was reduced to 7.67%.The findings suggest that the combined action of five target genes in the oleic acid metabolic pathway prevented the metabolism of oleic acid to linoleic acid and linolenic acid, which led to an increase in oleic acid content and a decrease in linoleic acid and linolenic acid contents.

Comparison of Agronomic Traits
Under field management, the agronomic traits of T2 generation positive mutant plants and control JN38, T6098, and T7010 materials were compared.Some of the mutant plants with JN38 as the recipient differed significantly from the recipient plants in plant height and no other traits (Table 4).Note: Means followed by the same letter within a column are not significantly different according to the least significant difference (LSD) test at p ≤ 0.05.

Discussion
Endoplasmic reticulum omega-6 fatty acid desaturases (FAD2 enzymes) are a family of enzymes that play an important role in regulating the oleic acid content in plants [21,22].The FAD2 enzyme enhances the content of linoleic acid by inserting a double bond at the delta-12 (omega-6) position of oleic acid [23].The FAD2 gene is the coding gene that regulates the FAD2 enzyme [24], which has been studied in a variety of plants with high oil content.For example, the down-regulated expression of two homologous FAD2 genes in peanut, based on the CRISPR/Cas9 method, resulted in increased oleic acid levels in stably transformed seeds, compared to the wild type [20].Soybean seeds with SACPD-C and FAD2-1A mutant alleles showed stearic acid contents from 14% to 21%, while the oleic acid content was increased [25].CRISPR/Cas9 was used to knock out the GmFATB1 gene, or one or both of them, to produce the mutant.The individual mutants FATB1A and FATB1B reduced the palmitic acid and stearic acid contents in leaves by 11% to 21%.In the double mutant FATB1A:1B, weak acid and stearic acid levels were reduced by 42% and 35%, respectively [26].
In previous studies, it was found that the bivalent and polyvalent CRISPR/Cas9 knockout of the same family gene often had a higher regulatory effect on oleic acid than the single gene knockout [16,19].The soybean FAD2 gene family exists in two forms: Fad2-1 and Fad2-2 [27].Two FAD2-1 members, GmFAD2-1B and GmFAD2-1A, were selected in this experiment.Three FAD2-2 members, GmFAD2-2B, GmFAD2-2C, and GmFAD2-2D, were studied.CRISPR/Cas9 gene-editing technology was used to construct a multi-gene knockout expression vector that knocked out five target genes of the FAD2 family.The fatty acid content of the mutant progeny showed that the oleic acid content of the mutant plants with different receptors was higher than that of the receptor plants, and that the linoleic acid content was lower than that of the receptor plants.The polyvalent knockout vector JN38 increased the oleic acid content from 18.97% to 51.14%, and the linoleic acid content was reduced from 57.79% to 43.36%.
The fatty acid content of soybean Is a quantitative character, which is influenced by many factors such as the environment and heredity.Editing different FAD2 genes also increases the oleic acid content to varying degrees.In a study of the soybean FAD2 gene family, the oleic acid improvement degree was higher than that of a double knockout vector test with knockout soybeans GmFAD2-1B and GmFAD2-2C [3], and the transformed offspring of the FAD2-1 gene were silenced using RNAi interference [28].We believe that the FAD-2 family genes may interact with each other, and that regulation of the fatty acid ratio is more obvious when the expression of the FAD-2 gene is regulated.The new allele FAD2-1A was induced by Jo et al. in soybean mutant PE529; they showed that oleic acid content in soybean oil can be increased to 80%.The results were significantly higher than those of this experiment, which may be due to the influence of different soybean lines on genetic factors and different experimental environments [29].Due to the different levels of expression, and the inhibition of the FAD2 gene family in different plants [30], the degree of increase in the oleic acid content varies.However, the trend of increasing the oleic acid content and decreasing the linoleic acid content is consistent.Different crops can also express genes in the FAD2 gene family.By reducing the expression, the oleic acid content is increased and the linoleic acid content is decreased.The results of this study are consistent with those of previous studies, and they provide theoretical support for the FAD2 gene family regulating the oleic acid synthesis pathway and new soybean varieties with high oleic acid contents.
The FAD2 gene family has a wide range of functions.In addition to the obvious regulation of oleic acid function, some studies have shown that FAD2 genes may also affect plant growth and development.The soybean FAD2 gene is expressed differently in different plant tissues, which improves its ability to resist cold and salt stress by regulating the contents of diene fatty acids [20].It has also been observed that combinations of polyvalent mutants do not seem to affect other agronomic characteristics, such as plant morphology and other quality traits [16,25].In this study, various agronomic traits of plants of different genotypes were investigated and analyzed during the whole growth period, and it was found that the plant height of the JN38-mutant was significantly higher than that of receptor plants, which might be due to differences in varieties.

Materials
The JN38, Z6098, and Z7010 soybean lines selected in this study belong to soybean varieties with excellent agronomic traits.As the receptor material of this study, they were provided by the Biotechnology Center of Jilin Agricultural University.

Construction of the CRISPR/Cas9 Vector
The high-resolution gRNA sequence was selected from the CRISPR-P website at OOV in September 2000.The target sequence mainly included 20 bases and a PAM structure, making Cas9 more convenient to bind to gRNA.The off-target effects of the two target sequences were screened using the BLAST program on the NCBI website, and the results show that there was no matching off-target sequence.U6 was designed as a promoter, and the CRISPR/Cas vector construction kit provided by Jiangsu Weimi Biotechnology Company was used to construct the vector.
Amounts of 2 µL CRISPR/CAS9 vector, 2 µL oligo dimer, 1 µL enzyme mixture, and ddH 2 O were mixed on ice to a total volume of 10 µL, then mixed and reacted at room temperature (20 • C) for 1 h and transformed into Escherichia coli.DNA was extracted from the recombinant plasmid for PCR verification (Table S1) and Sanger sequencing, and then compared using DNAMAN software (V6.0) to verify whether the gRNA was successfully connected to the vector.

Genetic Transformation
Using the cotyledon segments of soybean JN38, Z6098, and Z7010 as receptor materials, the editing vector was transferred to the receptor through the Agrobacterium tumefaciensmediated method.The tissue of the soybean cotyledon segments was cut 2-3 times and placed in a nutrient solution containing Agrobacterium tumefaciens for germination, screening, differentiation, rooting, and other cultures.It was also inoculated on a medium rich in other plant hormones to promote its growth.The genomic DNA of T2 soybean was amplified by PCR using primers Cas9-1F/Cas9-1R and Bar-1F/Bar-1R (Table S1).Complete PCR amplification was undertaken to a final reaction volume of 20 µL PCR conditions: 94 • C, 5 min; 30 cycles, 94 • C for 30 s, 60 • C for 30 s, 72 • C for 10 min.

Bar-Positive Plant Detection Test Strip
Bar-positive plant detection test strips were obtained from Wuhan Shangcheng Biotechnology Co., Ltd.(Wuhan, China); a total of eight T0 generation positive soybean plants were discovered using Bar test strips to test T0 generation plants.When a control-line appears, this is a negative control, indicating that the transformation was not successful.When both the control-line and T-line appear, this indicates that the sample has been successfully transformed into soybean.

Evaluation of Transformation Efficiency and Gene-Editing Efficiency
Specific primers were designed to include target sites.Genomic DNA was extracted from positive soybean for PCR amplification.The amplicons were sequenced by Changchun Kumei (Changchun, China).DNAMAN (v6.0) was used to compare and analyze the generated sequences.The base changes in the gRNA sequence indicated that the target gene was successfully edited.Statistical data in respect of soybean transformation and gene-editing efficiency were obtained.The amino acid sequence of transgenic positive soybean obtained via Sanger sequencing was used as a template.The SWISS-MODEL website (http://swissmodel.expasy.org/,accessed on 22 January 2021) was visited, and the amino acid sequence was submitted for template identification.

Southern Blot Analysis of the T2 Transgenic Plants
Genomic DNA was extracted from the leaves of T2 generation transgenic soybeans and digested using restriction endonucleases Hind III and BamH I. We used exogenous gene Cas9 as a probe, DIG-HIGH Prime DNA labeling and detection kit I (C58-11745832910, Roche, Shanghai, China) for detection, and conducted Southern blot analysis experiments.

Real-Time Quantitative PCR (RT-qPCR)
The total RNA of soybean seeds and leaves that were positive according to PCR testing of the T2 generation was extracted via RNAiso-Pluskits (TaKaRa, Beijing, China).A fluorescent quantitative reverse transcription kit (TakaraBio, RR036A, Beijing, China) was used to reverse-transcribe into cDNA, which was diluted 5-fold.The reaction procedure was 95 • C for 30 s, 95 • C for 10 s, 60 • C for 30 s, and 40 cycles.The melting point curve and data analysis were performed using the 2 −∆∆CT method [31].The soybean β-actin gene (GenBank accession number: NM_001252731.2) was selected as the internal reference gene and the online design primer website (Primer Quest Tool) (https://sg.idtdna.com/pages/tools (accessed on 10 September 2022)) (Table S1) was used.To assess the reproducibility and stability of our experiment, we selected 3 biological replicates.For the selection of RNA concentration, we utilized a concentration of 431.7 ng/µL.In order to eliminate the possibility of DNA contamination, we added DNase I (70 µL Reaction Buffer + 10 µL DNaseI) before using RNA for qPCR experiments.To evaluate RNA purity, we employed a Nano Drop 2000 microspectrophotometer (Thermo Fisher, Waltham, MA, USA) instrument for detection.Prime-ScriptTM RT Master Mix (Fermentas, Waltham, MA, USA) reagent was used for the amplification of qRT-PCR.The Mx3000P fluorescence quantitative PCR instrument (Agilent Technologies, Lexington, MA, USA) was used for analysis.The relative expression levels of the target gene in the leaves and the mature seed of the transformed plants were detected, and the relative expression levels were calculated using the formula 2 −∆∆CT .

Determination of Main Oil Content
The seeds of each plant from the T2 generation were collected and analyzed using gas chromatography (GC).The method of Yang [21] and Zhang [15] was slightly modified and tested.First, the soybean seeds were ground into a powder and dried overnight in a 90 • C oven.An amount of 2 mL sodium methoxide (0.5 M L −1 ) was added to each flask, and the sample was placed in a 30 • C water bath for 30 min.After transesterification with sodium methoxy, the lipids were extracted with n-hexane.The samples were then separated for 45 min and the fatty acid methyl ester was separated using HP-6890 gas chromatography with a hydrogen flame ionization detector and an HP INNOWAX capillary column (30 m; 0.25 mm id) with an N 2 carrier at 2 mL/min.The chromatographic conditions were as follows: initial temperature 100 • C; final temperature, 250 • C; injection temperature 250 • C; rates, 15 • C min −1 to 250 • C. Samples were maintained at 250 • C to compare the retention time of fatty acid methyl ester (FAME) with that of the standard (37-component FAME mix: Supelco 47885-U).Relative fatty acid compositions were then calculated as the percentage that each fatty acid represented of the total measured fatty acid.

Analysis of the Main Agronomic Traits
The standards for yield traits include the number of pods per plant, grain weight per plant, and 100 grain weight.The number of pods was determined by measuring the average of three plants per row.We determined the weight of grains by measuring the average of three plants, and we determined the weight of 100 grains by measuring the weight of a single plant of uniform size.We repeated the measurements (g) three times.Other agronomic traits evaluated included plant height, number of branches, and number of nodes.We measured the plant height (cm) of positive soybean plants in the laboratory.We calculated the number of branches and nodes on the main stem.These data were analyzed using IBM SPSS statistical software.

Statistical Analysis
All of the experiments in this study were replicated three times.For statistical analysis of experimental measurement data, SPSS 23.0 software was utilized, and unidirectional ANOVA was performed to confirm the variability of results between treatments.Nonsignificance (ns) was set at * p ≤ 0.05 and ** p ≤ 0.01.

Figure 4 .
Figure 4. Positive detection of some T0 generation soybean plants.WT: wild-type plant; 1-3: T generation transformed plants to be tested.

Table 1 .
Target design information, the yellow part is the PAM structure.

Table 3 .
Determination of fatty acid content of different editing vectors in T2 plants.
Note: Means followed by the same letter within a column are not different according to the least significant difference (LSD) test at p ≤ 0.05.

Table 4 .
Comparison of the agronomic characters of T2 generation plants.