Next Article in Journal
Parameter Optimization of Finger Clip Plate Garlic Seed-Metering Device
Previous Article in Journal
In Vitro Regeneration of Chrysanthemum from Ovaries and Ovules Treated with Thermal and Chemical Stimuli: Morphogenic and Cytogenetic Effects
 
 
Correction published on 31 May 2024, see Agriculture 2024, 14(6), 880.
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Identification of a Novel KTi-1 Allele Associated with Reduced Trypsin Inhibitor Activity in Soybean Accessions

1
Department of Applied Plant Science, Chonnam National University, Gwangju 61186, Republic of Korea
2
BK21 Interdisciplinary Program in IT-Bio Convergence System, Chonnam National University, Gwangju 61186, Republic of Korea
3
Department of Applied Biosciences, Kyungpook National University, Daegu 41566, Republic of Korea
4
Department of Crop Science, College of Agriculture and Life Sciences, Chungnam National University, Daejeon 34134, Republic of Korea
5
National Agrobiodiversity Center, National Institute of Agricultural Sciences, RDA, Jeonju 54874, Republic of Korea
*
Authors to whom correspondence should be addressed.
These authors contributed equally to this work.
Agriculture 2023, 13(11), 2070; https://doi.org/10.3390/agriculture13112070
Submission received: 4 October 2023 / Revised: 26 October 2023 / Accepted: 27 October 2023 / Published: 28 October 2023 / Corrected: 31 May 2024
(This article belongs to the Section Crop Genetics, Genomics and Breeding)

Abstract

:
Trypsin inhibitors (TIs) in soybean seeds reduce the availability of processed soybean foods and animal feed. This study aimed to evaluate the trypsin inhibitor activity (TIA) in 999 Korean soybean accessions and conduct molecular characterization of soybean accessions with low TIA. TIA was evaluated using colorimetric analysis through a substrate–enzyme reaction. The average TIA of the 999 soybean accessions was 90.31%, ranging from 43.70% to 99.51%. Kunitz trypsin inhibitor 3 (KTi-3) gene target sequencing analysis was performed on seven soybean germplasm accessions (IT105782, IT170889, IT273590, IT274513, IT274515, IT276197, and IT022891) showing less than 60% TIA. Four soybean accessions (IT274513, IT274515, IT276197, and IT022891) exhibited the same mutations (a G-to-T transversion and an AG deletion) in the KTi-3 gene (Glyma.08g341500) in PI542044. In addition, whole-genome re-sequencing was performed on three soybean accessions with no mutations in the KTi-3 gene. Compared with the reference soybean genome, an in-frame insertion and five missense mutations were identified in the coding sequencing of the KTi-1 gene (Glyma.01g09500) in IT105782. A RT-qPCR analysis showed that the mRNA expression level of KTi-1 was reduced by approximately 17% in IT105782 during seed development. In this study, we identified a previously unreported mutation in the KTi-1 gene and developed a KASP marker using this new KTi-1 variant.

1. Introduction

Soybean seeds (Glycine max (L.) Merr.) are valuable sources of nutrition for both humans and animals. Because of the high demand for protein meals and vegetable oils, global soybean production is expected to increase by 6% year-on-year, reaching 361 million tonnes. Brazil (139 million tonnes) and the United States (120.8 million tonnes) were the world’s largest soybean producers in 2021–2022 [1]. Soybean seeds comprise approximately 40% protein, 30% carbohydrates, and 20% fat, with the remaining 10% containing various functional substances [2]. In addition, soybeans contain non-nutritive elements such as isoflavones, which reduce blood cholesterol levels and increase bone strength [3]. Because of their health benefits, soybeans are widely consumed in various processed forms, such as soybean paste, soy sauce, tofu, and soy milk, in Korea, China, and Japan.
However, soybeans also contain anti-nutritional elements such as lipoxygenase, tannins, trypsin inhibitor (TI), and phytates [4,5]. Especially, TI interferes with the amino acid hydrolysis reactions of trypsin and chymotrypsin, strongly inhibiting protein digestion [6]. Therefore, the presence of TI in animal feed causes digestive disturbances, pancreatic enlargement, and liver damage in monogastric animals [7,8,9]. TI makes up about 6% of the total soybean protein and can be divided into two major types: the Kunitz trypsin inhibitor (KTi) and the Bowman–Brik inhibitor (BBi) [10]. In the context of trypsin inhibition, KTi plays a more substantial role compared to BBi [11,12]. When evaluating trypsin inhibition on an equal-mass basis, KTi exerts a more pronounced influence than BBi [13]. This suggests that KTi exhibits specificity in TI. KTi is present in higher concentrations than BBi, and KTi content varies widely among different soybean genotypes, ranging from 4.28 to 15.0 mg per gram of soy flour [14,15]. Soybean KTi is a protein with a 21 kDa molecular weight that exists as a monomer, consisting of 181 amino acids linked by two disulfide cysteine residues, while BBi is a low-molecular-weight protein (7 kDa) consisting of 71 amino acids with seven disulfide cysteine residues [16,17]. The disulfide bonds in TIs increase their stability at high temperatures and decrease the availability of soybean processed foods [18]. To decrease the activity of TIs, high-temperature treatment at 100 °C for 40 min is sufficient, but it causes nutrient denaturation and requires significant effort [18].
Three genes (KTi-1, KTi-2, and KTi-3) were cloned from developing soybean seeds [19,20]. The KTi-1 gene (Glyma.01g095000) is located on chromosome 1, while the KTi-3 gene (Glyma.08g341500) is located on chromosome 8. The KTi-2 gene, on the other hand, is not included in the most recent soybean genome annotation [21]. The KTi-1 and KTi-3 genes are composed of a single CDS domain and exhibit about 80% sequence homology [20].
The mRNA expression of KTi is mainly observed in the endosperm and cotyledon of seedlings, with rare occurrences in other regions [20]. In soybean seeds, the expression level of KTi protein is first detected at a seed length of 3–4 mm and reaches its peak in 11–12 mm seeds, gradually decreasing thereafter [22]. The KTi-3 gene accounts for more than 80% of total KTi expression [23].
Several soybean accessions with low KTi activity have been found in soybean germplasm. According to Orf and Hymowitz, KTi activity is reduced by 40% to 50% in PI157440 and PI196168 in the USDA homogeneous germplasm collection [24]. These two accessions contain the same mutations (one G-to-T transversion and two nucleotide sequence deletions) in the KTi-3 gene, resulting in a truncated protein [25]. In addition, Gillman et al. (2015) identified PI68679 with reduced KTi content among 520 soybean accessions in the USDA soybean germplasm collection. PI68679 has a single-base deletion in the KTi-1 gene (GGG356 → GG358), resulting in a frameshift mutation [21]. The KTi-1 null allele shows a reduction in trypsin inhibitor activity (TIA) similar to that of the KTi-3 null allele, and the combination of the KTi-1 and KTi-3 null alleles results in a greater reduction in TIA [21]. Recently, a low KTi genotype, Punjab1, was identified in Indian soybean germplasm [14]. However, the specific genetic variation contributing to the low KTi phenotype in Punjab1 has not yet been determined.
Several molecular markers have been developed to assist in the selection of null-KTi soybean lines without the need for biochemical analysis [26,27]. Since Kim et al. [28] identified three simple sequence repeat (SSR) markers (Satt228, Satt409, and Satt429) closely linked to the KTi-3 null allele within a range of 0–10 cM, these SSR markers have been broadly used to introgress the KTi-3 null allele in soybean breeding programs [26,27]. Recently, the single-nucleotide polymorphism (SNP) loci of the KTi-3 null allele have been employed to develop functional markers based on high-throughput and automated SNP genotyping platforms such as the Simple Probe assay and the Kompetitive allele-specific polymerase chain reaction (KASP) genotyping assay [21,29].
In general, sodium dodecyl sulfate-polyacrylamide gel electrophoresis (SDS-PAGE) has been utilized to identify KTi in soybeans [30]. On the other hand, the Lab-on-a-Chip, ELISA, two-dimensional liquid chromatography, and high-performance liquid chromatography (HPLC) methods have been developed to quantify KTi [21,31,32,33,34]. However, these methods need specialized equipment. Currently, a colorimetric enzyme assay based on the substrate–enzyme reaction mode is widely used for large-scale analysis [25,35]. DL-BAPA (Na-Benzoyl-D, L-arginine 4-nitroanilide hydrochloride) is used as a synthetic substrate, and the bovine trypsin enzyme cleaves it to release a yellow-colored p-Nitroaniline, which shows strong UV/V absorption at 416 nm [35]. A high amount of TI in the sample extract reduces TI, leading to a decrease in p-Nitroaniline, while a low amount of TI in the sample extract leads to an increase in p-Nitroaniline [35]. Substrate–enzyme reaction methods are classified into the E-last test and the S-last test. The E-last test involves reacting the inhibitor and the substrate first and then adding the enzyme to confirm the colorimetric reaction. In contrast, the S-last test involves reacting the inhibitor and the enzyme first and then adding the substrate [36]. The stability of the experimental results differs depending on the order in which the substrate and the enzyme bind to the TI. In a recent study, it was confirmed that the E-last test shows consistent results between pH 3.0 and pH 9.0, while the S-last test shows higher-than-average measured values above pH 7.0 [37].
In this study, we adapted Liu’s method to analyze 999 soybean accessions rapidly at a low cost, with the goal of screening for low TIA and identifying novel KTi-null mutations for the development of molecular markers [37]. We successfully screened low TIA, identified the KTi-1 mutation, and developed a KASP marker. The newly developed screening marker is efficient in selecting for low TIA traits, which will aid in the breeding of KTi-null soybeans in the future.

2. Materials and Methods

2.1. Plant Materials

In this study, a total of 999 soybean accessions were obtained from the Genebank of the National Agrobiodiversity Center at the Rural Development Administration (RDA) in the Republic of Korea. All accessions were grown in 2019 in the experimental field of the RDA Genebank in Jeonju, Republic of Korea (located at 35°49′53″ N, 127°03′50″ E). After sowing two plants per hole in a 90 cm × 15 cm spaced plot, one plant was removed, and three plants were harvested. We acquired 30 seeds for each accession from the Genebank, all of which were harvested in 2019, and used them directly for TIA assay. The controls consisted of two germplasms: Daepung 2 (IT269977) (KTi-1 and KTi-3) as the wild type and PI542044 with the null KTi-3 mutation. PI542044, a soybean accession with a null KTi-3 allele, was created through five backcrossing generations involving ‘Williams 82’ and PI157440 at both the Illinois Agricultural Experimentation Station and the USDA-ARS in the USA. PI542044 was obtained from the USDA Germplasm Resources Information Network (GRIN) (http://www.ars-grin.gov, accessed on 26 October 2023). On 21 June 2021, one wild-type soybean, null KTi-3 soybean and eight low-TIA soybean accessions were planted in Chonnam National University Field (35°10′32″ N, 126°54′25″ E). Two plants were placed in one hole at 20 cm intervals on 0.7 m × 1.5 m plots. These plants were used for DNA extraction, RNA extraction, and agronomic trait evaluation.

2.2. TIA Assay

This experiment modified Keshun Liu’s TIA to test 15 samples simultaneously in three replications, streamlining large-scale analysis with minimal resource utilization [37]. The changes included reducing the volume to 1 mL and performing the reactions in a dry bath at a constant temperature of 40 °C. In this experiment, 10 mg of seed powder, sieved using a 500 μm mesh sieve, was mixed with 1.5 mL of 10 mM NaOH buffer. The protein extract was obtained after centrifugation at 10,000 rpm for 30 min. In a 2 mL deep-well plate, a specified volume of 100 mM Tris-HCl buffer (pH 8.2, containing 20 mM CaCl2), protein extract, and 40 °C BAPA Working Solution were added and allowed to react for 10 min. Trypsin solution was added to react at 40 °C for 20 min, and the reaction was stopped with a 30% acetic acid solution.

2.3. Calculation of TIA

The TIA (%) was calculated by measuring at 410 nm after the completion of the reaction. The data were displayed as the average ± standard deviation (SD). TIA indicates the activity of TI as a percentage. To obtain TIA (%), absorbance values for four conditions (Reference, Reference Blank, Sample, and Sample Blank) were calculated and substituted into the formula below [38,39]. TIA was measured in triplicate, and absorbance was evaluated at 410 nm using an Epoch Microplate Spectrometer (BioTek instrument, Inc. Winooski, VT, USA).
T I A % = A R A R B A s A S B A R A R B × 100 %
AR, absorbance of Reference;
ARB, absorbance of Reference Blank;
AS, absorbance of Sample;
ASB, absorbance of Sample Blank.

2.4. DNA Extraction and Kompetitive Allele-Specific PCR (KASP) Assay

Fresh leaves were harvested from one wild-type soybean plant and eight low-TIA soybean plants in 2021. DNA was isolated using a DNeasy Plant Kit (Qiagen, Hilden, Germany), and its concentration and purity were analyzed using a Nanodrop ND 2000 spectrophotometer (Thermo Fisher Scientific, Waltham, MA, USA). KASP primers were designed using candidate gene sequences to distinguish the types of KTis in soybean mutants (Table 1). KASP assays were conducted in a final reaction volume of 10 μL, consisting of 5 μL of KASP master mix, 0.14 μL of primer mix, and 5 μL of 50 ng/μL DNA. The KASP PCR cycling involved an initial pre-reading cycle at 30 °C for 1 min, a holding stage at 94 °C for 15 min, 10 touchdown cycles at 94 °C for 20 s starting from 61 °C and reducing by 0.6 °C per cycle down to 55 °C for 1 min, 35 cycles at 94 °C for 20 s and 55 °C for 1 min, and a post-PCR reading stage at 30 °C for 1 min, with subsequent analysis of end-point fluorescence data using StepOne software (Applied Biosystems).

2.5. Whole-Genome Re-Sequencing

Re-sequencing analysis was performed on low-TIA soybeans for which no mutations were found by target gene sequencing. Using an Illumina DNA library preparation kit, we constructed three low-TIA soybean sequencing libraries and sequenced them on an Illumina NovaSeq 6000 platform. The trimmed reads were aligned to the soybean reference genome [40] using the BWA-MEM [41], resulting in the generation of a SAM format file for SNP analysis. After filtering, the SNP matrix was generated using the SAM tools (v.0.1.16) [42]. SNP annotation was performed based on the genetic location information of the soybean reference genome. Reads with over 90% identity were categorized as “homozygous SNPs”; those with similarity ranging from 40% to 60% were classified as “heterozygous SNPs”. In addition, “N” was utilized to classify mixed cases, while “n” was utilized to indicate cases where read mapping was omitted.

2.6. RT-qPCR Assay

Total RNA was isolated from two seed stages: mid maturation (200–300 mg, seed fresh weight) and late maturation (400–500 mg, seed fresh weight). Total RNA extraction was conducted as per the manufacturer’s guidelines using the RibospinTM Plant (GeneAll, Republic of Korea). RNA quality was determined by measuring the absorbance at 260 nm using a Nanodrop 2000 spectrophotometer (Thermo Scientific, Waltham, MA, USA), and the RNA concentration of the samples was adjusted to 50 ng/µL. Complementary DNA (cDNA) synthesis was carried out with the SuperScriptTM III First-strand Synthesis SuperMix kit (Invitrogen, Carlsbad, CA, USA).
Quantitative PCR was performed on the ABI StepOnePlus system (Applied Biosystems, Waltham, MA, USA) using iQTM SYBR Green Supermix (Bio-RAD). The qPCR reaction was carried out with 2 µL of cDNA (100 ng), 10 µL of cDNA (100 ng), 10 µL of 2X SYBR® Green Supermix, 3 µL of KTi-1 primer (Table 2; 300 nM), and 3 µL of D.W. in a total reaction volume of 20 µL. Quantitative PCR assays were performed under the following conditions: preheating for 1 min at 95 °C, denaturation for 3 min at 95 °C, and 40 cycles of incubation at 95 °C for 15 s and at 60 °C for 1 min for amplification. Following that, a melting curve analysis was carried out to validate both the absence of a product and the formation of dimers among the primers. The relative expression levels of target genes were calculated using the previously reported 2−ΔΔCT method [40], and cons7 (insulin-degrading enzyme, metalloprotease) was used as a housekeeping gene [43,44].

2.7. Evaluation of Agronomic Traits

We investigated the agronomic traits of three accessions that did not have the KTi-3 mutation. In 2021, agronomic trait evaluations were conducted for three accessions, following the criteria introduced by the RDA’s GenBank for soybean genetic resource characteristic evaluations. These evaluations aimed to assess various agronomic traits, such as plant height, number of branches, number of nodes, pod count, seed yield per pod, 100-seed weight, flowering date, and maturity date. Each evaluation was repeated five times to ensure accuracy.

2.8. Statistical Analysis

The statistical analyses were conducted using IBM SPSS Statistics 27 (IBM, Armonk, NY, USA) [45] Student’s t-tests were employed to compare the impact of genotypes on gene expression. The t-test results were indicated with an asterisk (*) for a significant difference with p < 0.05, with two asterisks (**) for a significant difference with p < 0.01, and with “ns” for non-significant differences. Agronomic trait results were analyzed using a one-way analysis of variance (ANOVA) and the Duncan test.

3. Results

3.1. Analysis of TIA in 999 Soybean Accessions

A substrate enzyme analysis was performed to measure the activity of TIs in 999 soybean accessions. The TIA values for these accessions ranged from 43.70% to 99.51%, with an average of 90.31% (Figure 1 and Table S1). The majority of soybean accessions (625) exhibited a high level of TIA (>90%), while 263 accessions showed a moderate level of TIA (80–90%). However, seven soybean accessions, namely, IT274515, IT273590, IT274513, IT170889, IT276197, IT022891, and IT105782, showed less than 70% TIA (Table 3). Among the 999 soybean germplasms, the lowest TIA was detected in IT274515. In addition, five accessions, excluding IT022891 and IT105782, exhibited lower TIA levels compared to PI542044 (null KTi-3).

3.2. Results of KASP Assay Screening for KTi-3 Allele

A comparison of the base sequences between Daepung 2 and PI542044 confirmed the presence of two mutations: the 436G → T missense mutation and the 440GAG → G frameshift mutation. Based on the null KTi-3 mutation found in PI542044, we developed a KASP marker targeting the 436G → T mutation in the KTI-3 gene. To distinguish KTi-3 mutations, we analyzed seven accessions, including the KTI-null PI542044 and the wild type Daepung 2. The analysis showed that IT274513, IT274515, IT276197, and IT02289 had the same mutation as PI542044 (Figure 2). Furthermore, Daepung 2, IT105782, IT170889, and IT273590 were identified as having the wild-type allele in the KTi-3 gene.

3.3. Whole-Genome Re-Sequencing

A resequencing analysis was performed on three soybean accessions (IT105782, IT170889, and IT273590) that showed no KTi-3 gene mutations in the KASP analysis. Genome sequencing of the three accessions generated over 587.86 million reads (88.77 Gb), with an average of 195.96 million reads (29.59 Gb) per accession (Table 4). After the cleaning step, 95.75% to 95.81% of the reads were mapped to the reference soybean genome. On average, each soybean accession covered 95.79% of the reference genome. The mean sequencing depth ranged from 26.32× in IT273590 to 27.38× in IT170889, with an average of 26.90×. A comparison of the three soybean accessions with the reference sequence revealed a total of 5,619,625 variants, including 4,649,462 SNPs and 970,163 indels. Specifically, IT105782 had 1,583,974 SNPs, 333,575 indels, and 1,917,549 total variations. Furthermore, IT170889 had 1,561,108 SNPs and 320,010 indels, while IT273590 had 1,504,380 SNPs and 316,578 indels.

3.4. Novel Allele in KTi-1 Gene

Based on the NGS data, no significant variations were found in the KTi genes of IT170889 and IT273590. However, we found that IT105782 had five base-substitution mutations and an in-frame insertion of 3 bp in the KTi-1 gene (Figure 3 and Figure S1). The five nucleotide substitutions resulted in missense mutations, specifically causing the alterations G → A at position 154, changing Gly to Ser; A → C at position 392, changing Lys to Thr; C → G at position 484, changing Gln to Glu; G → A at position 529, changing Asp to Asn; and C → G at position 580, changing Gln to Glu (Figure S2). These missense mutations were all non-conservative. Additionally, an in-frame insertion of CTA at position 461 resulted in the generation of tyrosine. Overall, the results of genome sequencing indicate that the KTi-1 mutations in IT105782 were different from those in PI68679 (GGG356 → GG358, relative to the start codon; Figure 3A).

3.5. Analysis of mRNA Expression Levels of KTi-1 Genes

To evaluate the role of the KTi-1 gene in soybean seeds, we analyzed its expression levels at two different stages of seed development and compared them to both IT105782 and the wild type genotype. The analysis revealed a significant decrease in the expression of the KTi-1 gene in IT105782 (a KTi-1 mutant) during both the mid- and late stages (Figure 4). In the mid-stage, IT105782 showed a significant decrease of 63% compared to Daepung 2. Additionally, a significant decrease of 17% was observed in the late stage compared to Daepung 2.

3.6. Detection of Mutations in KTi-1 Candidate Genes Using KASP Markers

The KASP marker was based on the KTi-1 variant candidate gene Glyma01g095000 and was designed using the missense mutation at location 392. The newly developed KASP marker effectively distinguished between the KTi-1 mutant and the non-variant alleles (Figure 5).

3.7. Evaluation of Agronomic Traits

The agricultural characteristics of four soybean accessions harvested in 2021, namely, IT105782, IT170889, IT273590, and Daepung 2, were evaluated. As shown in Figure 6 and Table 5, the three soybean accessions (IT105782, IT170889, and IT273590) exhibited higher plant height compared to Daepung 2, accompanied by more pods and branches. Daepung 2 showed a flowering date that was 10.6 days earlier than the average of the three accessions. IT105782 showed the quickest maturation, while IT273590 exhibited the slowest. The average 100-seed weight (16.9 g) of the three soybean accessions was half that of Daepung 2 (30.88 g), and the seed numbers per pod were lower than that of Daepung 2. In general, the agronomic traits of the three soybean accessions were deemed unsuitable for use as soybean cultivars.

4. Discussion

Soybean TIs interfere with protein digestion in animals and reduce food availability [46]. Although the TIA of soybeans can be decreased by various processing methods, such as heat treatment, starch treatment, electromagnetic wave treatment, and cryogenic treatment, these treatments are expensive and may induce damage to other proteins [47,48,49]. Therefore, the most effective way to inactivate this nutritional element in soybean meals is to develop KTi-free soybean cultivars. Currently, three soybean accessions (PI157440, PI196168, and PI68679) lacking the KTi have been identified among the USDA soybean germplasm collection and extensively utilized in soybean breeding programs to introduce the null allele of the TI [26,50,51]. In this study, we screened many genotypes in Korean soybean germplasm collections to identify new soybean genotypes with low KTi contents and characterize novel genetic variations in KTi genes.
In this study, TIA analysis was conducted on a large number of soybean accessions using a colorimetric enzyme assay that relies on enzyme–substrate reactions [37]. However, enzyme–substrate reactions are significantly influenced by temperature and incubation time. Lowering the temperature of the reagent or reducing the reaction time may lower the TIA content. In this experiment, to maintain a consistent temperature, the assay buffer was pre-stored in a dry bath at 40 °C to create the BAPA working solution. Additionally, the protein extract was refrigerated at 2 °C. As the temperature of the protein extract approached room temperature, lower values were observed in soybean varieties with low KTI content. However, there was no significant reduction in varieties with high KTI content. This suggests that the modified assay in this study effectively screens for KTi and allows for the accurate and rapid screening of low TIA levels.
Among the 999 soybean accessions, more than 88% of genotypes exhibited a high level of TIA, ranging from 90% to 100% (Figure 1). Additionally, we identified six soybean accessions with low TIA levels (below 60%) similar to the KTi-free genotypes of PI542044 (Table 2). Among these six genotypes, three accessions (IT274515, IT273590, and IT274513) originated in Korea (Table 2). Interestingly, two KTi-3 null lines, PI157740 and PI196168, were also introduced into the USA from South Korea [52]. To determine if the low TIA levels in these six accessions were associated with KTi-3 null lines, the coding regions of the KTi-3 gene (Glyma.08g341500) were compared with PI542044. A KASP analysis revealed that IT274513, IT274515, IT276197, and IT022891 carried the same mutation as PI542044. Since these three accessions were recently deposited at the RDA GenBank in 2013–2014, the low KTi trait in IT274513, IT274515, and IT276197 might have originated from breeding programs utilizing these KTi-3 null lines.
Re-sequencing analysis is the most effective method for identifying genetic variations in candidate genes. Jofuku and Goldberg reported a minimum of 10 KTi genes in the soybean genome. [20]. In this study, KTi genes were identified from the Williams82 soybean reference genome (Wm82.a2.v1.genome). Using re-sequencing analysis, we examined the genetic variations of the gene in three soybean accessions (IT105782, IT170889, and IT273590) containing no KTi-3 gene mutation. Remarkably, we found that IT105782 contained a three-base-insertion mutation (A461 → ACTA) and five missense mutations in the coding region of the KTi-1 gene (Glyma.01g09500) (Figure 3). The insertion of an additional three-base pairs led to the insertion of one more amino acid (Y154 → YY154), without resulting in a frameshift change. The previously reported KTi-1 null line, PI68679, had a single-base deletion (GGG → 358GG) at position 356 bp, resulting in a frame shift and a truncated protein [21]. Therefore, the mutation found in IT105782 differed from PI69679, suggesting that it is a new mutant allele of the KTi-1 gene.
To evaluate the impact of the novel KTi-1 gene mutation on mRNA levels, we conducted RT-qPCR analysis on mature seeds. The results showed a significant reduction (17–63%) in the expression level of the KTi-1 gene in IT105782 compared to the wild type Daepung 2 during the mid-maturation and late-maturation stages (Figure 4). This reduction suggests a direct correlation between the reduced KTi-1 mRNA level and the decline in TI protein production. Moreover, it is crucial to consider the suitability of these three soybean varieties for cultivation. An assessment of various agricultural characteristics, including average seed weight, seed count, and other quantitative aspects, indicated that these varieties may not be ideal choices for cultivation. In the development of genetic markers related to KTi-1, some efforts have faced challenges and been unsuccessful [21]. Despite these challenges, we successfully developed the KTi-1 KASP marker, demonstrating its effectiveness in distinguishing between the mutation type and the wild type.

5. Conclusions

In this study, we identified several soybean accessions with low TIA from a pool of 999 soybean germplasm accessions. Among them, IT105782 had an in-frame insertion (A461 → ACTA) mutation and five missense mutations in the Glyma01g095000 gene. This mutation is a novel KTi-1 gene mutation that has not been previously reported. RT-qPCR analysis revealed that the KTi-1 mRNA expression level in IT105782 significantly decreased by 17%-63% compared to Daepung 2 (an elite soybean variety) during the two seed stages. IT105782 showed early maturation characteristics similar to those of null-KTi soybean, although its other agronomic traits were not superior to those of Daepung 2 (an elite cultivar soybean). This study identified a new KTi-1 mutation and developed a new KTi-1 molecular marker. Breeder-friendly KASP markers for selecting KTi-null plants will soon be readily available in numerous breeding programs.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/2077-0472/13/11/2070/s1. Figure S1: Base sequence of IT105782 in KTi-1 Glyma01g095000 gene. Sequence variations are shown in square brackets, with the front-base sequence representing reference and the back-base sequence representing IT105782; Figure S2: Amino acid differences between IT105782 and reference in KTi-1 Glyma01g095000 gene; Table S1: Trypsin inhibitor activity (TIA) in 999 soybean accessions.

Author Contributions

Writing—original draft preparation, A.P. and S.-H.K.; methodology, B.-H.K., S.C., S.-Y.S. and W.-H.L.; resources, J.-D.L., S.L. and Y.-M.C.; funding acquisition, J.-D.L., S.L., Y.-M.C. and B.-K.H.; supervision, Y.-M.C. and B.-K.H. All authors have read and agreed to the published version of the manuscript.

Funding

This work was carried out with the support of the “Cooperative Research Program for Agriculture Science and Technology Development (Project No. PJ014168/RS-2019-RD008897)” of the Rural Development Administration, Republic of Korea.

Institutional Review Board Statement

Not applicable.

Data Availability Statement

The datasets generated during the current study are available from the corresponding author upon reasonable request.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Dahl, E.; Nelson, R. Oilseeds. Agric. Commod. 2021, 11, 53–61. [Google Scholar]
  2. Liu, K. Chemistry and nutritional value of soybean components. In Soybeans; Springer: Boston, MA, USA, 1997; pp. 25–113. [Google Scholar]
  3. O’Keefe, S.F.; Bianchi, L.; Sharman, J. Soybean nutrition. SM J. Nutr. Metab. 2015, 1, 1006. [Google Scholar]
  4. Adeyemo, S.; Onilude, A. Enzymatic reduction of anti-nutritional factors in fermenting soybeans by Lactobacillus plantarum isolates from fermenting cereals. Niger. Food J. 2013, 31, 84–90. [Google Scholar] [CrossRef]
  5. Hwang, C.; Lee, S.; Kang, J.; Kwon, M.; Kwon, H.; Chung, J.; Sung, N. Physicochemical characteristics and antioxidant activity of Kanjang made from soybean cultivars lacking lipoxygenase and kunitz trypsin inhibitor protein. J. Agric. Life Sci. 2012, 46, 111–125. [Google Scholar]
  6. Kunitz, M. Crystalline soybean trypsin inhibitor II. General properties. J. Gen. Physiol. 1947, 30, 291–310. [Google Scholar] [CrossRef]
  7. Liener, I.E.; Goodale, R.L.; Deshmukh, A.; Satterberg, T.L.; Ward, G.; DiPietro, C.M.; Bankey, P.E.; Borner, J.W. Effect of a trypsin inhibitor from soybeans (Bowman-Birk) on the secretory activity of the human pancreas. Gastroenterology 1988, 94, 419–427. [Google Scholar] [CrossRef] [PubMed]
  8. Palacios, M.; Easter, R.; Soltwedel, K.; Parsons, C.; Douglas, M.; Hymowitz, T.; Pettigrew, J. Effect of soybean variety and processing on growth performance of young chicks and pigs. J. Anim. Sci. 2004, 82, 1108–1114. [Google Scholar] [CrossRef]
  9. Smith, J.; Wilson, F.; Allen, P.; Berry, D. Hypertrophy and hyperplasia of the rat pancreas produced by short-term dietary administration of soya-derived protein and soybean trypsin inhibitor. J. Appl. Toxicol. 1989, 9, 175–179. [Google Scholar] [CrossRef]
  10. DiPietro, C.; Liener, I. Soybean protease inhibitors in foods. J. Food Sci. 1989, 54, 606–609. [Google Scholar] [CrossRef]
  11. DiPietro, C.M.; Liener, I.E. Heat inactivation of the Kunitz and Bowman-Birk soybean protease inhibitors. J. Agric. Food Chem. 1989, 37, 39–44. [Google Scholar] [CrossRef]
  12. Friedman, M.; Brandon, D.L.; Bates, A.H.; Hymowitz, T. Comparison of a commercial soybean cultivar and an isoline lacking the Kunitz trypsin inhibitor: Composition, nutritional value, and effects of heating. J. Agric. Food Chem. 1991, 39, 327–335. [Google Scholar] [CrossRef]
  13. Yang, Y.; Chang, S.K.C.; Zhang, Y. Determination of protease inhibitors, glycinin, and beta-conglycinin in soybeans and their relationships. J. Food Sci. 2022, 87, 1082–1095. [Google Scholar] [CrossRef] [PubMed]
  14. Kumar, V.; Rani, A.; Mittal, P.; Shuaib, M. Kunitz trypsin inhibitor in soybean: Contribution to total trypsin inhibitor activity as a function of genotype and fate during processing. J. Food Meas. Charact. 2019, 13, 1583–1590. [Google Scholar] [CrossRef]
  15. Pesic, M.B.; Vucelic-Radovic, B.V.; Barac, M.B.; Stanojevic, S.P.; Nedovic, V.A. Influence of different genotypes on trypsin inhibitor levels and activity in soybeans. Sensors 2007, 7, 67–74. [Google Scholar] [CrossRef]
  16. Sweet, R.; Wright, H.; Janin, J.; Chothia, C.; Blow, D.A. Crystal structure of the complex of porcine trypsin with soybean trypsin inhibitor (Kunitz) at 2.6 Å resolution. Biochemistry 1974, 13, 4212–4228. [Google Scholar] [CrossRef]
  17. Voss, R.H.; Ermler, U.; Essen, L.O.; Wenzl, G.; Kim, Y.M.; Flecker, P. Crystal Structure of the Bifunctional Soybean Bowman-Birk Inhibitor at 0.28-nm Resolution: Structural Peculiarities in a Folded Protein Conformation. Eur. J. Biochem. 1996, 242, 122–131. [Google Scholar] [CrossRef]
  18. Chen, Y.; Xu, Z.; Zhang, C.; Kong, X.; Hua, Y. Heat-induced inactivation mechanisms of Kunitz trypsin inhibitor and Bowman-Birk inhibitor in soymilk processing. Food Chem. 2014, 154, 108–116. [Google Scholar] [CrossRef] [PubMed]
  19. Jofuku, K.D.; Schipper, R.D.; Goldberg, R.B. A frameshift mutation prevents Kunitz trypsin inhibitor mRNA accumulation in soybean embryos. Plant Cell 1989, 1, 427–435. [Google Scholar]
  20. Jofuku, K.D.; Goldberg, R.B. Kunitz trypsin inhibitor genes are differentially expressed during the soybean life cycle and in transformed tobacco plants. Plant Cell 1989, 1, 1079–1093. [Google Scholar]
  21. Gillman, J.D.; Kim, W.-S.; Krishnan, H.B. Identification of a New Soybean Kunitz Trypsin Inhibitor Mutation and Its Effect on Bowman− Birk Protease Inhibitor Content in Soybean Seed. J. Agric. Food Chem. 2015, 63, 1352–1359. [Google Scholar] [CrossRef]
  22. McGrain, A.K.; Chen, J.C.; Wilson, K.A.; Tan-Wilson, A.L. Proteases catalysing processing and degradation of Kunitz soybean trypsin inhibitor during seed maturation. Phytochemistry 1992, 31, 421–426. [Google Scholar] [CrossRef]
  23. De Moraes, R.M.A.; Soares, T.C.B.; Colombo, L.R.; Salla, M.F.S.; de Almeida Barros, J.G.; Piovesan, N.D.; De Barros, E.G.; Moreira, M.A. Assisted selection by specific DNA markers for genetic elimination of the kunitz trypsin inhibitor and lectin in soybean seeds. Euphytica 2006, 149, 221–226. [Google Scholar] [CrossRef]
  24. Orf, J.; Hymowitz, T. Inheritance of the Absence of the Kunitz Trypsin Inhibitor in Seed Protein of Soybeans 1. Crop Sci. 1979, 19, 107–109. [Google Scholar] [CrossRef]
  25. Liu, K.; Markakis, P. An improved colorimetric method for determining antitryptic activity in soybean products. Cereal Chem. 1989, 66, 415–422. [Google Scholar]
  26. Maranna, S.; Verma, K.; Talukdar, A.; Lal, S.K.; Kumar, A.; Mukherjee, K. Introgression of null allele of Kunitz trypsin inhibitor through marker-assisted backcross breeding in soybean (Glycine max L. Merr.). BMC Genet. 2016, 17, 1–9. [Google Scholar] [CrossRef] [PubMed]
  27. Rani, A.; Kumar, V.; Shukla, S.; Jha, P.; Tayalkar, T.; Mittal, P. Changes in storage protein composition on genetic removal of Kunitz trypsin inhibitor maintain protein content in soybean (Glycine max). J. Agric. Food Res. 2020, 2, 100065. [Google Scholar] [CrossRef]
  28. Kim, M.S.; Park, M.J.; Jeong, W.H.; Nam, K.C.; Chung, J.I. SSR marker tightly linked to the Ti locus in soybean [Glycine max (L.) Merr.]. Euphytica 2006, 152, 361–366. [Google Scholar] [CrossRef]
  29. Rosso, M.L.; Shang, C.; Song, Q.; Escamilla, D.; Gillenwater, J.; Zhang, B. Development of Breeder-Friendly KASP Markers for Low Concentration of Kunitz Trypsin Inhibitor in Soybean Seeds. Int. J. Mol. Sci. 2021, 22, 2675. [Google Scholar] [CrossRef]
  30. Duranti, M.; Barbiroli, A.; Scarafoni, A.; Tedeschi, G.; Morazzoni, P. One-step purification of Kunitz soybean trypsin inhibitor. Protein Expr. Purif. 2003, 30, 167–170. [Google Scholar] [CrossRef]
  31. Brandon, D.L.; Bates, A.H.; Friedman, M. ELISA analysis of soybean trypsin inhibitors in processed foods. In Nutritional and Toxicological Consequences of Food Processing; Springer: Boston, MA, USA, 1991; pp. 321–337. [Google Scholar]
  32. Rosso, M.L.; Shang, C.; Correa, E.; Zhang, B. An Efficient HPLC Approach to Quantify Kunitz Trypsin Inhibitor in Soybean Seeds. Crop Sci. 2018, 58, 1616–1623. [Google Scholar] [CrossRef]
  33. Rouquié, D.; Capt, A.; Eby, W.H.; Sekar, V.; Hérouet-Guicheney, C. Investigation of endogenous soybean food allergens by using a 2-dimensional gel electrophoresis approach. Regul. Toxicol. Pharmacol. 2010, 58, S47–S53. [Google Scholar] [CrossRef]
  34. Zhou, T.; Han, S.; Li, Z.; He, P. Purification and quantification of Kunitz trypsin inhibitor in soybean using two-dimensional liquid chromatography. Food Anal. Methods 2017, 10, 3350–3360. [Google Scholar] [CrossRef]
  35. Call, L.; Reiter, E.V.; Wenger-Oehn, G.; Strnad, I.; Grausgruber, H.; Schoenlechner, R.; D’Amico, S. Development of an enzymatic assay for the quantitative determination of trypsin inhibitory activity in wheat. Food Chem. 2019, 299, 125038. [Google Scholar] [CrossRef] [PubMed]
  36. Kakade, M.; Rackis, J.; McGhee, J.; Puski, G. Determination of trypsin inhibitor activity of soy products: A collaborative analysis of an improved procedure. Cereal Chem. 1974, 51, 376–381. [Google Scholar]
  37. Liu, K. Soybean trypsin inhibitor assay: Further improvement of the standard method approved and reapproved by American oil Chemists’ Society and American Association of Cereal Chemists International. J. Am. Oil Chem. Soc. 2019, 96, 635–645. [Google Scholar] [CrossRef]
  38. Liu, K. Trypsin Inhibitor Assay: Expressing, Calculating, and Standardizing Inhibitor Activity in Absolute Amounts of Trypsin Inhibited or Trypsin Inhibitors. J. Am. Oil Chem. Soc. 2021, 98, 355–373. [Google Scholar] [CrossRef]
  39. Van Eys, J.; Offner, A.; Bach, A. Manual of quality analyses for soybean products in the feed industry. Am. Soybean Assoc. Bruss. Belg. 2004. [Google Scholar]
  40. Schmutz, J.; Cannon, S.; Schlueter, J.; Ma, J.; Mitros, T.; Nelson, W.; Hyten, D. Genome sequence of the palaeopolyploid soybean. Nature 2010, 463. [Google Scholar] [CrossRef]
  41. Li, H. Aligning sequence reads, clone sequences and assembly contigs with BWA-MEM. arXiv 2013, arXiv:1303.3997. [Google Scholar]
  42. Li, H.; Handsaker, B.; Wysoker, A.; Fennell, T.; Ruan, J.; Homer, N.; Marth, G.; Abecasis, G.; Durbin, R.; Subgroup, G.P.D.P. The sequence alignment/map format and SAMtools. Bioinformatics 2009, 25, 2078–2079. [Google Scholar] [CrossRef]
  43. Livak, K.J.; Schmittgen, T.D. Analysis of relative gene expression data using real-time quantitative PCR and the 2−ΔΔCT method. Methods 2001, 25, 402–408. [Google Scholar] [CrossRef] [PubMed]
  44. Libault, M.; Thibivilliers, S.; Bilgin, D.; Radwan, O.; Benitez, M.; Clough, S.; Stacey, G. Identification of four soybean reference genes for gene expression normalization. Plant Genome 2008, 1, 44–54. [Google Scholar] [CrossRef]
  45. Levesque, R. SPSS Programming and Data Management. A Guide for SPSS and SAS Users; SPSS Inc.: Chicago, IL, USA, 2007. [Google Scholar]
  46. Makkar, H.P.; Siddhuraju, P.; Becker, K. Trypsin Inhibitor. In Plant Secondary Metabolites; Springer: Boston, MA, USA, 2007; pp. 1–6. [Google Scholar]
  47. Baker, E.; Mustakas, G. Heat inactivation of trypsin inhibitor, lipoxygenase and urease in soybeans: Effect of acid and base additives. J. Am. Oil Chem. Soc. 1973, 50, 137–141. [Google Scholar] [CrossRef]
  48. Barać, M.; Stanojević, S. The effect of microwave roasting on soybean protein composition and components with trypsin inhibitor activity. Acta Aliment. 2005, 34, 23–31. [Google Scholar] [CrossRef]
  49. Vagadia, B.H.; Vanga, S.K.; Raghavan, V. Inactivation methods of soybean trypsin inhibitor–A review. Trends Food Sci. Technol. 2017, 64, 115–125. [Google Scholar] [CrossRef]
  50. Kumar, V.; Rani, A.; Rawal, R.; Mourya, V. Marker assisted accelerated introgression of null allele of kunitz trypsin inhibitor in soybean. Breed. Sci. 2015, 65, 447–452. [Google Scholar] [CrossRef]
  51. Schmidt, M.A.; Hymowitz, T.; Herman, E.M. Breeding and characterization of soybean Triple Null; a stack of recessive alleles of Kunitz Trypsin Inhibitor, Soybean Agglutinin, and P34 allergen nulls. Plant Breed. 2015, 134, 310–315. [Google Scholar] [CrossRef]
  52. Krishnan, H.B. Characterization of a soybean [Glycine max (L.) Merr.] mutant with reduced levels of Kunitz trypsin inhibitor. Plant Sci. 2001, 160, 979–986. [Google Scholar] [CrossRef]
Figure 1. Distribution of trypsin inhibitor activity (TIA) in 999 soybean germplasms.
Figure 1. Distribution of trypsin inhibitor activity (TIA) in 999 soybean germplasms.
Agriculture 13 02070 g001
Figure 2. KASP assay screening for KTi-3 using HEX and FAM fluorescent fluorophores. Shown in red are the soybean lines with KTi-3 mutation presented by HEX fluorescence, and in blue are the soybean lines with normal KTi-3 allele presented by FAM fluorescence.
Figure 2. KASP assay screening for KTi-3 using HEX and FAM fluorescent fluorophores. Shown in red are the soybean lines with KTi-3 mutation presented by HEX fluorescence, and in blue are the soybean lines with normal KTi-3 allele presented by FAM fluorescence.
Agriculture 13 02070 g002
Figure 3. (A) Overall visual description of the reference sequence, KTi-1 mutation, and the novel KTi-1 mutation. (B) Detailed visual description of the novel KTi-1 mutation.
Figure 3. (A) Overall visual description of the reference sequence, KTi-1 mutation, and the novel KTi-1 mutation. (B) Detailed visual description of the novel KTi-1 mutation.
Agriculture 13 02070 g003
Figure 4. KTi-1 mRNA expression levels of three soybean accessions with low Kunitz trypsin inhibitor (KTi) activity. **, p < 0.01.
Figure 4. KTi-1 mRNA expression levels of three soybean accessions with low Kunitz trypsin inhibitor (KTi) activity. **, p < 0.01.
Agriculture 13 02070 g004
Figure 5. KASP assay for KTi-1 using HEX and FAM fluorescent fluorophores. Shown in red are the soybean lines with KTi-1 mutation concentration accumulation, indicated by HEX fluorescence, and in blue are the differences from KTi-1 cultivars, indicated by FAM fluorescence.
Figure 5. KASP assay for KTi-1 using HEX and FAM fluorescent fluorophores. Shown in red are the soybean lines with KTi-1 mutation concentration accumulation, indicated by HEX fluorescence, and in blue are the differences from KTi-1 cultivars, indicated by FAM fluorescence.
Agriculture 13 02070 g005
Figure 6. Plants (A) and seeds (B) of three soybean accessions with low Kunitz trypsin inhibitor (KTi) activity.
Figure 6. Plants (A) and seeds (B) of three soybean accessions with low Kunitz trypsin inhibitor (KTi) activity.
Agriculture 13 02070 g006
Table 1. List of primers used in KASP.
Table 1. List of primers used in KASP.
MarkerGene ID
(Glyma.Wm82.gnm2)
Fluorescent PrimerSequence
KTi-1Glyma01g095000FAM_primerCTACTGTATCGCGTGCAGCAAGTT
HEX_primerTACTGTATCGCGTGCAGCAAGTG
Common reverse primerGAGAGAGAGGGTCTACAAGCTGTTA
KTi-3Glyma08g341500FAM_primerTGCAATGGATGGTTGGTTTAGACTTG
HEX_primerATGCAATGGATGGTTGGTTTAGACTTT
Common reverse primerCTGTGGACAGAACACAAGCTTATAGTTAT
Table 2. List of primer sequences used for target sequencing and RT-qPCR.
Table 2. List of primer sequences used for target sequencing and RT-qPCR.
GeneGene ID (Glyma.Wm82.gnm2)Forward Primer (5′ → 3′)Reverse Primer (5′ → 3′)Product SizePurpose
KTi-1Glyma01g09500TGGAGAGAGAGGGTCTACAAGCTTAGACAGCACTAAACGCCTGA196 bpRT-qPCR
cons7NoneATGAATGACGGTTCCCATGTAGGCATTAAGGCAGCTCACTCT114 bpRT-qPCR
Table 3. Trypsin inhibitor activity (TIA) percentage values in seven soybean accessions and reference soybeans (mean ± standard deviation, n = 3).
Table 3. Trypsin inhibitor activity (TIA) percentage values in seven soybean accessions and reference soybeans (mean ± standard deviation, n = 3).
Accession NameGenotype NameTIA (%)Origin
IT27451501450243.70 ± 0.21KOR
IT273590MNG-PSARI-1998-1153.10 ± 2.67KOR
IT27451301449954.50 ±4.38KOR
IT170889Hood 7562.90 ± 0.57USA
IT276197GSI 01409964.15 ± 5.71Unknown
IT022891PI 7486664.75 ± 1.34Unknown
IT105782Kongnamul Kong68.98 ± 1.89KOR
IT269977Daepung 2 (normal KTi)89.15 ± 0.22KOR
PI542044Kunitz (null KTi-3)66.60 ± 1.41USA
Table 4. Summary of genome re-sequencing and DNA sequence variants identified in three soybean accessions with low Kunitz trypsin inhibitor (KTi) activity.
Table 4. Summary of genome re-sequencing and DNA sequence variants identified in three soybean accessions with low Kunitz trypsin inhibitor (KTi) activity.
IT105782IT170889IT273590MeanTotal
Total reads196,617,960199,490,702191,750,912195,953,191587,859,574
Average sequence length (bp)151151151151151
Total size29,689,311,96030,123,096,00228,954,387,71229,588,931,89188,766,795,674
Total size (Mb)29,68930,12328,95429,588.788,766
Sequencing depth (×)26.9927.3826.3226.90-
Mapped reads186,915,887189,267,121182,758,171186,313,726558,941,179
Genome coverage (%)95.81%95.75%95.80%95.79%-
Number of SNPs1,583,9741,561,1081,504,3801,549,8214,649,462
Number of indels333,575320,010316,578323,388970,163
Total variation1,917,5491,881,1181,820,9581,873,2085,619,625
Table 5. Major agronomic traits of three soybean accessions with low Kunitz trypsin inhibitor (KTi) activity and the control cultivar.
Table 5. Major agronomic traits of three soybean accessions with low Kunitz trypsin inhibitor (KTi) activity and the control cultivar.
CultivarPlant Height (cm)Number of BranchesNumber of NodesNumber of PodsSeeds per Pod100-Seed Weight (g)Flowering DateMaturity Date
Daepung 250.8 ± 3.27 a†2.6 ± 0.55 a12.4 ± 2.07 a57.0 ± 10.44 a330.9 ± 1.60 d29 July27 October
IT10578270.0 ± 5.05 b11.0 ± 1.87 b15.4 ± 1.14 b120.4 ± 14.84 b212.7 ± 0.74 a2 August15 October
IT17088988.6 ± 1.34 c13.0 ± 2.35 c15.8 ± 0.84 b178.4 ± 12.78 c322.9 ± 0.70 c6 August 6 November
IT273590133.2 ± 9.83 d9.4 ± 2.70 b22.0 ± 1.22 c122.0 ± 38.94 b215.1 ± 0.61 b18 August 15 November
The different letters indicate a significant difference at p < 0.05 by ANOVA and Duncan’s test.
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

Park, A.; Kang, S.-H.; Kang, B.-H.; Chowdhury, S.; Shin, S.-Y.; Lee, W.-H.; Lee, J.-D.; Lee, S.; Choi, Y.-M.; Ha, B.-K. Identification of a Novel KTi-1 Allele Associated with Reduced Trypsin Inhibitor Activity in Soybean Accessions. Agriculture 2023, 13, 2070. https://doi.org/10.3390/agriculture13112070

AMA Style

Park A, Kang S-H, Kang B-H, Chowdhury S, Shin S-Y, Lee W-H, Lee J-D, Lee S, Choi Y-M, Ha B-K. Identification of a Novel KTi-1 Allele Associated with Reduced Trypsin Inhibitor Activity in Soybean Accessions. Agriculture. 2023; 13(11):2070. https://doi.org/10.3390/agriculture13112070

Chicago/Turabian Style

Park, Aron, Se-Hee Kang, Byeong-Hee Kang, Sreeparna Chowdhury, Seo-Young Shin, Won-Ho Lee, Jeong-Dong Lee, Sungwoo Lee, Yu-Mi Choi, and Bo-Keun Ha. 2023. "Identification of a Novel KTi-1 Allele Associated with Reduced Trypsin Inhibitor Activity in Soybean Accessions" Agriculture 13, no. 11: 2070. https://doi.org/10.3390/agriculture13112070

APA Style

Park, A., Kang, S. -H., Kang, B. -H., Chowdhury, S., Shin, S. -Y., Lee, W. -H., Lee, J. -D., Lee, S., Choi, Y. -M., & Ha, B. -K. (2023). Identification of a Novel KTi-1 Allele Associated with Reduced Trypsin Inhibitor Activity in Soybean Accessions. Agriculture, 13(11), 2070. https://doi.org/10.3390/agriculture13112070

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