Bowman-Birk Inhibitor Mutants of Soybean Generated by CRISPR-Cas9 Reveal Drastic Reductions in Trypsin and Chymotrypsin Inhibitor Activities

1. Introduction

Soybean Glycine max (L.) Merr) is grown worldwide as a rotation partner for grain crops, particularly maize, due to an ability to fix atmospheric nitrogen, high seed yield potential, relatively few major biotic pathogens and relatively low input costs [1,2]. Soybean seeds contain far more protein than other legumes (on average 40%) and soybean meal is often the first choice for swine and poultry feed mixtures due to a balanced (albeit suboptimal) amino acid profile. Soybean meal, produced after extraction of seed oil, has a relatively low cost and is of high quality compared to other competing protein meal sources [3].

Despite the high quality of soybean protein, raw soybeans and soybean meal cannot be directly included into animal feed mixtures due to anti-nutritional compounds that interfere with digestion, cause pancreatic inflammation, and effectively reduce animal weight gain. Soybean seeds contain two major classes of proteinaceous anti-nutritional factors [4], which were named after the scientists that discovered them: 1) Kunitz-trypsin inhibitors [5] (KTi) inactivate the digestive enzyme trypsin; and 2) Bowman-Birk inhibitors [6] (BBi), which inactivate both trypsin and chymotrypsin animal digestive enzymes [7]. In soybean, three KTi isoforms (KTi1, KTi2 and KTi3) have been reported [8,9,10]. There are at least 11 different BBi encoding genes present in the soybean genome, however expression analysis of RNAseq data reveals only three of them are expressed at appreciable levels in developing seeds (Figure S1).

Experiments have conclusively shown that ingestion of high levels of active protease inhibitors are responsible for reduced animal weight gain, and can even result in hypoglycemia, pancreatic hypertrophy, or liver damage [11]. To avoid these issues and permit inclusion of soybean meal at high levels in animal feed, KTi and BBi anti-nutritional proteins in typical soybean seeds must be inactivated by heat-treatment, which has been shown to reduce, but not eliminate ANF activity [7]. Heat inactivation allows soybean meal to be included at moderately high levels in chicken and swine feed mixtures. While this step meets the needs of the industry, such treatment is energy intensive, adds expense, and requires a balance between sufficiently reducing ANF activity and negatively impacting the quality of seed proteins [3,11].

One success story in genetically reducing soybean’s trypsin inhibitor activity was the identification and analysis of a soybean accession (PI 157740) which features marked reduced (∼ 40%) in trypsin inhibitor activity. This genetic source has been examined through both genetic analysis and animal feeding trials. Raw extruded protein meal lacking KTi3 protein was found to result in superior animal weight gain compared with raw soybeans [12], but not a degree that could eliminate the need for heat inactivation. In our previous research, a soybean line (PI 68679) with reductions in trypsin inhibition due to a mutation affecting the KTi1 gene was identified [13]. Using hybridization, the KTi-1 mutation (PI 68679) was combined with the previously identified KTi-3 mutation (PI 542044), to yield soybean lines that lack both KTi1 and KTi3 protein. These soybean lines have lower trypsin inhibitor activity than the currently available KTi-3 mutant [13]. Unfortunately, all lines which have reduced KTi levels unexpectedly had unacceptable levels of trypsin inhibitor activity, due to overaccumulation of BBi proteins in seeds of all mutant lines (KTi1, KTi3, and KTi1/KTi3).

To develop soybean lines which require extremely low or no heat-inactivation processing, we hypothesize that both KTi and BBi gene expression in seeds must be reduced or ideally eliminated. Naturally occurring BBi mutants have been identified in perennial soybean collections [14]. However, to our knowledge no naturally occurring G. max mutant or germplasm entry has been reported with substantive reduced BBi accumulation or activity. However, one report indicated that transgenic expression of a modified protein can lower BBi activity and resulted in viable plants which were able to produce seed [15].

CRISPR/Cas9 gene editing [16,17], is now a routine practice for inducing targeted loss of function mutations in soybean [17,18,19]. This method could be leveraged to induce mutations of soybean seed BBi genes and decrease or eliminate their accumulation in seeds. In this study, we report on the creation of an allelic series of loss-of-function mutations affecting the major BBi seed expressed genes in soybean, and mutation of two genes was shown to lead to substantial reductions in both trypsin and chymotrypsin inhibitor activities due to loss of function mutations affecting BBi genes.

2. Results

2.1. Generation of BBi Mutants by CRISPR/Cas9 Gene Editing

At least 11 BBi encoding genes exist in the soybean genome, however transcriptome analysis implicates only one major and two minor seed-specific BBi genes in developing soybean seeds (Supplementary Figure S1) [20]. Because BBi genes share extensive sequence homology we designed target sites that enabled simultaneous knockout of all three seed-expressed genes, and two CRISPR/Cas9 expression cassettes were designed, which differed in use of either single or double guideRNA (Supplementary Figure S2). Utilizing these two CRISPR/Cas9 expression cassettes, we generated several independent transgenic events by Agrobacterium-mediated soybean transformation (Figure 1A). Six transgenic events were grown in the greenhouse and leaf tissue from these plants were utilized for the isolation of genomic DNA. Sanger sequencing of PCR products generated by utilizing primers that are specific for the most abundantly expressed BBi gene (Glyma.16g208900) identified several mutations in this gene. Figure 1A shows a schematic summary of mutations identified in these six independent events leading to deletions, insertions, and frameshift mutations on amino acid sequence. (Figure 1B, Figure S3). Similarly, we also sequenced PCR products generated by utilizing primers that are specific for the second most abundantly expressed BBi gene (Glyma.14G117700). This analysis revealed several mutations that introduced deletions, insertions, and frameshift mutations on amino acid sequence (Figure 2A and 2B; Figure S3)

2.2. SDS-PAGE Analysis of Seed Proteins of CRISPR/Cas9 Transgenic Plants

First, we wanted to compare the total seed protein profile between the wildtype and CRISPR/Cas9 transgenic plants. For this purpose, we first isolated total seed proteins from the wild-type and the six independent transgenic events and resolved them by SDS-PAGE (Figure 3A). A visual examination of the total seed protein profile between wild type and transgenic events did not reveal any drastic alterations between them (Figure 3A). All soybean seeds examined accumulated the two most abundant seed storage proteins (7S β-conglycinin and 11S glycinin). However, the accumulation of the β-subunit of β-conglycinin, a 52 kDa glycoprotein, in all transgenic events were noticeably lower when compared to the wild type (Figure 3A). To substantiate these changes, we performed immunoblot analysis using antibodies raised against the 7S β-conglycinin and 11S glycinin antibodies. The 7S β-conglycinin antibody reacted against all three subunits. The accumulation pattern of the α’ and α subunits of β-conglycinin was similar between the wild type and the six CRISPR/Cas9 transgenic events (Figure 3B). In contrast, the accumulation of the β-subunit of β-conglycinin was much lower in all transgenic events when compared to the wild type, and for one transgenic event (Figure 3B, lane 7) only an extremely faint band was detected. We also performed immunoblot analysis to quantify any changes in the 11S glycinin accumulation between wild type and CRISPR/Cas9 transgenic events. The glycinin antibody recognized three glycinin species (A3; A1a, A1b, A2, and A4; A5). In contrast to the β-subunit of β-conglycinin, the accumulation of 11S glycinin showed little variation in their abundance between the wild type and the CRISPR/Cas9 transgenic events (Figure 3C).

2.3. Immunoblot Analysis Reveal the Absence of BBi in the Seeds of Transgenic BBi Mutants

To examine the accumulation of BBi in seeds of the CRISPR/Cas9 transgenic events we performed western blot analysis utilizing BBi-specific peptide-antibodies [21]. A protein fraction enriched for protease inhibitors obtained by calcium fractionation method [22] was utilized for the immunoblot analysis and transferred to a nitrocellulose membrane. The nitrocellulose membrane when challenged with soybean BBi-specific antibodies failed to detect the accumulation of this protein in all transgenic events. In contrast, the BBi antibody reacted against a 11-12 kDa protein in the wild-type plant (Figure 4B). We also examined the accumulation of KTi, which also contributes significantly to the overall protease inhibitor of soybean seed, by immunoblot analysis, via antibodies raised against soybean KTi [9]. Immunoblot analysis clearly demonstrated normal accumulation of KTi proteins in all the plants examined in our study (Figure 4C).

2.4. Chymotrypsin and Trypsin Inhibitor Activities Are Drastically Reduced in CRISPR/Cas9 gene Edited Transgenic BBi Mutants

Since our immunoblot analysis clearly showed the absence of BBi protein accumulation in the CRISPR/Cas9 gene edited transgenic BBi mutants we also wanted to examine protease inhibitor activities in these plants. The chymotrypsin inhibitor activity (CIU) was measured using N-Succinyl-Ala-Ala-Pro-Phe p-nitroanilide as a substrate and is expressed as CIU /mg seed powder. The wild-type soybean seed contained 17.2 CIU while in transgenic seeds the activity ranged from 2.89-15.9 CIU. Out of the six transgenic events examined in this study, four had drastic reductions; we saw a 2.8-6.0-fold reduction in CIU compared to wild type (Table 1). Since soybean BBi proteins can inhibit both chymotrypsin and trypsin, we also measured the trypsin inhibitor activity (TIU) using the same seed extracts (Table 1), with Nα-Benzoyl-DL-arginine 4-nitroanilide hydrochloride as a substrate. The trypsin inhibitor activity in the wild-type soybean seeds averaged 57.25 TIU/mg seed powder, whereas we saw reduced activities in mutant lines ranging from a low of 14.46 (Event 8) to a high of 35.7 TIU (Event 5) for different transgenic BBi events. This corresponds to a 1.6-4.0-fold reduction in TIU activities. TIU activity in all examined transgenic BBi mutants were significantly lower than in wild type seeds. Moreover, in four of the transgenic events (Events 2, 4, 7 and 8; Table 1) we observed drastic reductions in the TIU activities ranging from 25.3-41.9% of the TIU activities found in wild type seeds. Interestingly, the reduction in TIU activity in these BBi mutants is more significant than we observed for a naturally occurring soybean frameshift KTi-3 mutant (Table 1). CIU and TIU activities in seeds of mutant and the control line were fully concordant when activities were expressed as TIU/mg seed powder or as TIU/mg seed protein (Table 1).

2.5. qRT-PCR Analysis Confirms the Absence of Expression of the Abundant BBi Genes in CRISPR/Cas9 Gene Edited Transgenic BBi Mutants

We also examined the expression of the two most abundantly expressed BBi genes (Glyma.16G208900, Glyma.14G117700) along with few other BBi genes (Glyma.09G158500, Glyma.09G158600, Glyma.09G158700) in the CRISPR/Cas9 gene edited transgenic BBi mutants. Total RNA was isolated from mid-stage (35 days after flowering) developing T1 soybean seeds and used for qRT-PCR analysis. The gene specific primers used in our study are shown in Table 3. In the wild-type soybean seeds, we could detect the expression of most of these BBi genes. At this stage of seed development, we observed high level expression of Glyma.16G208900, Glyma.14G117700 and Glyma.09G158700 (Figure 5). In contrast, the expression of these abundantly expressed BBi genes was essentially not detected in two of the transgenic events (event # 2 and 4), while several fold lower expression was detected in one of the mutants (event # 8).

3. Discussion

In this study, we have generated several stable transgenic soybean lines that exhibit the lowest levels of trypsin and chymotrypsin inhibitor activities that have been reported in the literature. Several of our CRISPR-Cas9 BBi lines have much lower trypsin inhibitor activity (62% to 75% lower) when compared to control plants. The relative abundance of KTi proteins in soybean seeds is higher than BBi [23]. It has been reported that BBi proteins are more heat stable than KTi protein since the BBi contains multiple cysteine residues leading to conformational change of BBi during heating [24]. Since our CRISPR-Cas9 BBi lines contain undetectable levels of BBi proteins, as evidenced by immunoblot analysis, one would expect these newly developed BBi mutant lines may require only minimal heat treatment.

Until now, no soybean null mutants with reduced or ablated BBi proteins have been identified in G. max germplasm. However, BBi mutants have been identified in several perennial Glycine spp. [14]. Unfortunately, substantial hybridization barriers exist between perennial Glycine spp. and G. max which has hampered the development of G. max BBi mutants by traditional breeding. Attempts have been made to lower the protease inhibitor activity in soybean seeds through overexpression of a mutated BBi gene [25]. The transgenic seeds generated in that study revealed a 20-50% reduction in protease inhibitory activity. In contrast, we have developed several CRISPR-Cas9 mutated BBi lines that showed remarkable reduction (ranging from 63-84%) in chymotrypsin inhibitor activity when compared to the wild-type seeds. Even though our CRISPR-Cas9 constructs were targeted to knockout seed specific BBi genes, its pronounced effect on lowering the trypsin inhibitor is striking. This drastic reduction in trypsin inhibitor activity is not mediated by a reduction in the accumulation of KTi proteins in the CRISPR-Cas9 BBi lines. Our immunoblot analysis (Figure 4) clearly demonstrates that KTi proteins accumulate to the similar levels in BBi mutant lines as in wild type seeds. Thus, despite the presence of functional versions of all three KTi seed expressed genes, we saw a significant reduction in the trypsin inhibitor activity. This reduction is solely because of BBi mutations, and unexpectedly revealed that BBi contributes significantly to overall trypsin inhibitor activity in soybean seeds.

Our immunoblot analysis (Figure 4) and qRT-PCR analysis (Figure 5) confirmed that the gRNA we employed in our study to knockout the seed specific BBi genes was very effective. In addition to eliminating the expression of the two most abundant BBi genes (Glyma.16G208900 and Glyma.14G117700) other BBi genes (Glyma.09G158500, Glyma.09G158600, Glyma.09G158700) were also eliminated in some the CRISPR-Cas9 BBi induced mutant lines. This was an expected result as BBi genes share extensive sequence homology and constructs were designed to target multiple BBi genes. The full number of mutated genes is unknown and will likely require creation and analysis of segregating populations to discern the number of mutated genes and their relative contributions to trypsin and chymotrypsin inhibition activities. One intriguing prospect would be the combination of null mutations for both seed expressed KTi genes and our newly identified null mutations affecting all seed expressed BBi genes.

Our study also reveals proteome alteration due to knock out of BBi genes. One-dimensional gel analysis of total seed proteins from BBi mutant lines reveals lower accumulation of the β-subunit of β-conglycinin. The accumulation of this protein is influenced by nitrogen and sulfur status of the plant [26,27]. The accumulation of β subunit of β-conglycinin, a protein that is deficient in sulfur containing amino acids, is elevated when the plants are grown under limited sulfur supply [28,29,30]. In the current study, there is a drastic reduction in the accumulation of BBi in all CRISPR-Cas9 mutated BBi lines. BBi contains 14 cysteine residues and acts as a main sink for sulfur storage in the seed [31,32]. In the absence of BBi accumulation, cysteine that would be normally incorporated in BBi will be freely available. We speculate that the reduction in the accumulation of the β-subunit of β-conglycinin could be related to an increase in the availability of cysteine in the seeds, which would be incorporated into other, normally lower abundance seed expressed proteins. Our suggestion for this possibility is strengthened by previous reports that has shown an increased availability of sulfur amino acids in transgenic soybean plants reduces the accumulation of β subunit of β-conglycinin [33]. However, this hypothesis remains to be rigorously tested through greenhouse and/or field studies. Although we clearly demonstrate the creation of an allelic series of BBi null mutations through CRISPR Cas9 mutagenesis, we acknowledge we do not currently know what impact, if any, such mutations (or a combination of mutations) will have on field agronomic performance or seed yield. Similar work with KTi-3 mutant lines did not identify any deleterious impacts on whole plant agronomic traits or seed yield [34].

Our results open the possibility that breeding lines could be developed that require substantially less heat-treatment, which would save both money and energy which is used to heat treat soybeans. One study performed with KTi-3 and lectin free seed demonstrated that partial loss of trypsin/chymotrypsin inhibitors resulted in superior animal weight gain performance as compared to conventional soybean meal, however both were inferior to heat processed soybean meal [11]. Given our finding that our BBi mutants had reduced trypsin inhibitor activity equivalent to the KTi-3 mutation (Table 1), the BBi allelic series we have created may open the door to even greater reductions on trypsin and chymotrypsin inhibitor content in animal feed mixtures and for human food products, such as tofu and soymilk, while also reducing the energy required for heat processing.

4. Materials and Methods

4.1. Chemicals

Most chemicals and reagents used in this study were of analytical grade. Chemicals for electrophoresis including acrylamide, bis acrylamide, SDS, TEMED, ammonium persulfate were purchased from GE healthcare (Piscataway, NJ). Tris-HCL, trypsin, chymotrypsin, β-mercaptoethanol, Nα-Benzoyl-DL-arginine p-nitroanilide hydrochloride, and N-Succinyl-Ala-Ala-Pro-Phe p-nitroanilide were purchased from Sigma-Aldrich (St. Louis, MO).

4.2. Analysis of Expression of BBi Genes Using Public Data Repositories

Soybase.org is a repository for data from multiple experimental datasets for soybean. We identified protein coding genes in the ‘Williams 82’ genome soybean with the annotation PFAM PF00228 (Bowman-Birk serine protease inhibitor family) as targets. Previously collected data from the ‘RNA-Seq Atlas’ [20] was used to prioritize target BBi genes for CRISPR-Cas9 construct, gene sequencing and RT-PCR analysis. In developing seeds, two genes, Glyma.14g117700 and Glyma.16g208900 represent between 17.0-31.3% and 68.5-82.8% of the total detectable BBi gene expression from seeds samples between 21and 42 days after flowering.

4.3. Generation of CRISPR/CAS9 induced BBi Knockout Soybean Mutants

For CRISPR/Cas9-mediated mutagenesis of the BBi genes in soybean, two single gRNAs were designed to target the two major expressed BBi genes and three other BBi genes. gRNA1: (5’-GAACAACATGGTGGTGCTAAAGG-3’) was a common target in five BBi genes (Glyma.09G158500; Glyma.14G117700; Glyma.09G158600; Glyma.09G158700; Glyma.16G202800) and gRNA2 (5’-CACATGCAGAGATCACAGCATGG-3’, reverse complimentary of the sense strand sequence 5’- CCATGCTGTGATCTCTGCATGTG-3’) was a common target in only 3 BBi genes (Glyma.09G158600; Glyma.09G158700; Glyma.16G202800). Two plasmid vectors were constructed: the first plasmid contained only gRNA1 and the second plasmid contained both gRNA1 and gRNA2. For cloning of the BBi gRNAs in the T-DNA vector, a two-step Golden Gate Assembly was performed firstly in the module vector pMOD_B2103 and then subsequently in the destination T-DNA vector pTRANS230d as described previously [18,35]. Recombinant pMOD_B2103_BBi gRNAs plasmids harboring BBi gRNAs were verified by Sanger sequencing of the gRNA cassette with primers flanking the gRNAs: the forward primer TC320 (5’- CTAGAAGTAGTCAAGGCGGC-3’) and reverse primer TC089R (5’-GGAACCCTAATTCCCTTATCTGG-3’), respectively. The second Gloden Gate Assembly reaction was composed of the recombinant pMOD_B2103_BBi gRNAs, pMOD_A0521, pMOD_C2906, and a destination T-DNA vector pTRANS230d. Final CRISPR/Cas9_BBi gRNA constructs in the destination T-DNA vector were designated as pTRANS230/BBi_1gRNA and pTRANS230/BBi_2gRNA, respectively. The CRISPR constructs were mobilized into a disarmed Agrobacterium tumefaciens strain EHA105 by heat shock procedure.

4.4. Soybean Transformation

The soybean genotype ‘Maverick’ [36] was transformed using Agrobacterium-mediated transformation of half-seed explants. The transformation protocol was modified from Zeng et al., [37]. The process began by germinating sterilized seeds for 16 hours in Petri dishes on germination media. Then, half seed explants were made from the seeds and treated with an Agrobacterium culture [38]. After 5 days of co-cultivation, the explants were transferred to shoot induction medium. Subsequently, the explants underwent two rounds of subcultures before being moved into shoot elongation medium. Elongated shoots were then cut from the explants and used for rooting. After 2 weeks, the rooted plants were transferred into soil. To identify putative transgenic soybean plants, leaf paginating assay was performed with 100 mg L^− 1 glufosinate-ammonium solution on leaves at three stages. Herbicide resistant plants were cultivated in greenhouses under controlled conditions, with a 16-h day at 28°C and 8-h night at 24°C until T1 seeds were harvested.

4.5. Protein Extraction and 1-D Electrophoresis

Total soybean seed protein extraction and subsequent separation of proteins by 1-D SDS-PAGE were carried out as previously described [39]. Briefly, dry soybean seeds were ground to a fine power with a mortar and pestle and extracted with 1 mL of sodium dodecyl sulfate (SDS)-sample buffer (60 mM Tris-HCl, pH 6.8, 2% SDS, 10% glycerol and 5% 2-mercaptoethanol). Tubes were vortexed (Vortex Genie, Scientific Industries, San Diego, CA) for 10 min at room temperature, and then heated to 100° C for 5 min, followed by centrifugation (Eppendorf, Beckman Coulter, Inc., Indianapolis, IN) at 15800xg for 5 min. The supernatant is the total seed protein fraction. This fraction was resolved by 13.5% SDS-PAGE gels using a Hoefer SE 260 minigel apparatus (Amersham Biosciences, Piscataway, NJ). Electrophoretic separation of proteins was achieved with 20 mA per gel (constant current). SDS-PAGE separated proteins were visualized using Coomassie Blue R-250 solution (0.3 % Coomassie Brilliant Blue R-250, 45 % methanol, and 10 % glacial acetic acid).

4.6. Real-Time qRT-PCR Analysis of BBi Genes in Developing Soybean Seeds

Total RNA from 35-day after anthesis soybean seeds were extracted using TRIzol reagent (Invitrogen) according to the manufacturer’s instructions. Real-time qRT-PCR was conducted using a LightCycler 480 II instrument (Roche, Indianapolis, IN). Each qRT-PCR reactions were conducted with two replications in 20 µL reaction using the QuantiNova SYBR Green RT-PCR Kit (Qiagen, Hilden, Germany) following the manufacturer’s instructions. The cycling condition of absolute quantification analysis were 10 min at 50 °C then 2 min 95°C for cDNA synthesize followed by 40 cycles of 20 sec at 95 °C and 40 sec at 60 °C, with the data acquired at 60 °C in green channel. The reference gene (F-box protein, Glyma.12G051100) was amplified by RT-PCR from Maverick soybean seed total RNA with gene specific primer set (Supplementary Table S1). The amplified reference gene RT-PCR fragment was purified from agarose gel and standard was calculated based on purified fragment concentration and molecular weight. The reference gene was subjected to a 10-fold serial dilution (1.00E+1 to 1.00E-10) with DEPC. For absolute quantification analysis, five standard dilutions from 1.00E-5 to 1.00E-9 for standard curve were used for quantifying both target and reference gene expression in cDNA samples.

4.7. Screening of Putative Transgenic Plants and Identification of CRISPR-Cas9 induced Mutations

Twenty seeds from 10 individual T0 BBi events were planted in the Sears Plant greenhouse in 2-gallon pots containing Pro-Mix^® HP Mycorrhizae™ (Premier Tech, Quakertown, PA). The plants were fertilized with Osmocote Plus (Scotts) according to the manufacturer’s recommendations. Greenhouse settings were 16-h daylength and 30° C/18° C day/night temperatures. For every plant which germinated, fresh tissue was collected, and genomic DNA isolated from ~20-30 mg lyophilized leaf tissue using a Maxwell Promega robot with the Maxwell RSC PureFood GMO and Authentication Kit (Promega).

PCR primer pairs were designed to bear at least two gene-specific SNP differences as compared to homeologous genes. Primers used for PCR amplification and sequencing are detailed in supplementary Table S1. We blasted all primers against the unmasked Glycine max ‘Williams 82’ genomic sequence (E-value cut off = 10.0) to ensure primer specificity (www.soybase.org). PCR amplification was performed using Titanium taq according to manufacturer’s recommendations (Takara) in a PTC-200 thermocycler (MJ Research/Bio-Rad), using these conditions: 95°C for an initial 5 minute denaturation, followed by 40 cycles of 95°C for 30 seconds, followed by 60°C for 30 seconds, and an extension step at 72°C for 1 minute/kilobase of target sequence. PCR products were separated on a 1% agarose gel and compared to DNA ladder size standards to ensure appropriate size. PCR products were purified through use of a QIAquick PCR purification kit (Qiagen). After PCR purification, products were Sanger sequenced at the DNA Core Facility of the University of Missouri-Columbia.

PCR product sequencing traces were imported into Geneious Primer software, version 11.0.11+9 (Invitrogen), trimmed, aligned to gene sequences from the G275_W82.a2.v1 assembly (http://www.phytozome.net/soybean) using the Geneious Primer software, version 11.0.11+9 (Invitrogen). Disagreements were manually assessed in comparison to the G275_W82.a2.v1 assembly [40]. Putative indels and polymorphisms were verified by a minimum of two independent PCR reactions.

4.8. Immunoblot Analysis

Total seed proteins isolated from wild-type soybean cultivar ‘Maverick’ and T1:2 seeds of BBi mutated CRISPR/CAS9 lines were resolved on a 13.5% SDS-PAGE gels. Resolved proteins were transferred to nitrocellulose membranes and incubated with TBS (10 mM Tris-HCl, pH 7.5, 500 mM NaCl) containing 5% non-fat dry milk for 1 h at room temperature. Following this step, the nitrocellulose membranes were washed 3X with TBST (15 min each) and incubated over-night with 1:10,000 diluted antibodies raised either against purified soybean Kunitz trypsin or Bowman-Birk protease inhibitors. 7S β-conglycinin and 11S glycinin specific antibodies were used at 1:50,000 dilution. Specifically bound primary antibodies with detected using anti-rabbit IgG-horseradish peroxidase conjugates and a SuperSignal West Pico kit (Pierce, Rockford, IL, U.S.A.) as previously described [21].

4.9. Trypsin and Chymotrypsin Inhibitor Assay

Trypsin inhibitor activity in different BBi mutated CRISPR/CAS9 lines was measured following the procedure described by [41]. Briefly, to 20 mg of dry soybean seed powder was placed in a 2 mL Eppendorf centrifuge. To this 1 mL of 10 mM NaOH solution was added and vigorously agitated on a vortex mixer for 10 min at room temperature. The resulting slurry was clarified by centrifugation at 16,000xg for 10 min. Clear supernatant was utilized for measuring KTi activity. Trypsin inhibitor activity was expressed as trypsin units inhibited (TUI) per mg of sample, and results were expressed as the mean ± SD from three biological replicates.

Chymotrypsin inhibitor activity in soybean seed extracts were measured N-Succinyl-Ala-Ala-Pro-Phe p-nitroanilide (AAPF) as the substrate. Soybean seed extract, which resulted in 40 to 60% chymotrypsin inhibition range, was added to a 1.5 mL Eppendorf tube containing 900 µL of assay buffer (100 mM-Tris-HCL, pH 8.0), 8 µL of α-chymotrypsin (Sigma-Aldrich Company) dissolved in 1 mM HCl solution (0.1 mg/mL) and seed extract (8 to 10 µL) and incubated at 37 C for 5 min. Following this step, 80 µL of AAPF (1 mg/mL) was added and incubated for additional 10 min at 37 C. The assay was stopped by adding 500 µl of 30% acetic acid. Absorbance at 410 nm was measured, and chymotrypsin inhibitor units were calculated as the amount of inhibitor that reduced the absorbance at 410 nm by 0.1 optical density. The results are expressed as the mean ± SD from three biological replicates. Trypsin and Chymotrypsin inhibitor assay data was visualized and compared using JMP statistical software 16.0 (SAS Institute Inc., Cary, NC, USA). One-way ANOVA was performed and significant differences (α=0.05) between means were determined using the Tukey-Kramer HSD test.