Genome-Wide Microsatellite Survey and Introgression of Null KTI Allele in Charcoal Rot Resistant Soybean
- 1. Department of Agricultural Botany, Post Graduate Institute, Dr. Panjabrao Deshmukh Krishi Vidyapeeth, India
- 2. Regional Research Centre, Amravati, Dr. Panjabrao Deshmukh Krishi Vidyapeeth, India
- 3. Gen Scan Inc Global Consultancy, Houston, Texas, USA and Namangan Engineering and Technology Institute, Uzbekistan
- 4. Voronezh State Agrarian University, Voronezh, Russia
Abstract
Soybean though being a vegetable crop have not been used for table purpose due to various quality constraints. Kunitz trypsin inhibitor is a major ANF in soybean which interfere protein digestion and absorption. The objective of the experiment is formulated to address these issues with soybean that are preventing high-quality soybean production is necessary. For parental polymorphism, 388 SSR markers across all 20 linkage groups were examined in four parental combinations comprising four soybean genotypes: AMS MB-5-18, AMS MB-5-19, NRC 101, and NRC 127. Seventy-eight molecular markers showed polymorphism with 37 of them being highly polymorphic. Allelic frequency of 3.84, 4.3, 4.7 and 3.95 per marker was observed respectively for parental combinations 1,2,3 and 4 PIC value ranging from 0.11 to 0.75, 0.15 to 0.75, 0.25 to 0.75 and 0.13 to 0.75 respectively. Using a marker specific to the Null KTI gene, foreground selection was applied to the F1 progenies acquired from the experiment. Future soybean breeding programs will benefit from the obtained polymorphic markers.
Keywords
• Soybean; KTI
• Antinutritional factor
• Charcoal rot
• M. phaseolina
• Polymorphic information content
CITATION
Jadhav PV, Shinde AN, Thakur PP, Shinde UD, Mahobia A, et al. (2023) Genome-Wide Microsatellite Survey and Introgression of Null KTI Allele in Charcoal Rot Resistant Soybean. J Hum Nutr Food Sci 11(4): 1172.
INTRODUCTION
With high protein content, soybean is also referred to as functional food and is currently being highlighted for its potential to end malnutrition. According to the USDA, the worldwide soy food business was worth $44.7 billion in 2021. The soy food business is primarily driven by consumers’ changing preferences for protein-rich diets over high-calorie diets. The soy food items including textured vegetable proteins, tofu, and soy milk, have been in high demand on the international market.
However, its popularity has been limited because of antinutritional components like Kunitz trypsin inhibitor (KTI), beany flavor, and the presence of allergens. The anti-nutritional substances in soybean seed, especially from the two families of Kunitz trypsin inhibitors (KTI) and Bowman-Birk inhibitors, prevent the activity of digesting proteases (BBI). High quantities of raw soybean meal in feed mixes can result in poor body weight and pancreatic anomalies by inactivating trypsin/chymotrypsin enzymes. Heat-treating soybean protein meal in order to disable inhibitors is common practice that can deplete some essential amino acids, is expensive and energy-intensive [1] and 80% of the overall trypsin inhibitor activity is accounted for by KTI [2]. Despite the fact that KTI is heat labile and heat treatment costs money and reduces protein solubility by about 20%. Similarly, beans must be boiled before being ground with wheat (1:9) to make chapatti flour with soy supplements [3].
Genotypes with a null allele of KTI have been identified in soybean germplasm. KTI is a monomeric protein containing 181 amino acid residues [4]. Ten independent differentially expressed KTI genes have been identified with KTI3 being seed specific [5]. Thirteen KTI3 protein isozymes have also been identified which are encoded by multiple KTI alleles [6]. A single dominant gene (Ti) controls the presence of KTI, whereas a recessive gene (titi) controls its absence which encodes for a truncated protein. One substitution and two deletions are present in recessive null kti allele which alters the translation process resulting in lower levels of KTI in seed embryos [7]. Genetic research indicates that numerous alleles, notably Ti a, Ti b, Ti c, and Ti d, are combined at a single locus controlling the KTI trait in soybean. These four soybean KTI electrophoretic forms are controlled by a variety of allelic series of co-dominance (Ti a, Ti b, Ti c and Ti d). A final version ‘ti’ that lacks Kunitz trypsin inhibitor activity is controlled by a recessive gene [8]. The gene has been located in linkage group (LG) A2 on the genetic linkage map of soybean [9]. Three SSR markers namely Satt409, Satt228 and Satt429 are closely linked with the ti locus at distances of 4.5, 0-3.7, and 5.1 cM, respectively [10].
Abiotic and biotic stresses pose significant challenges to soybean production. Charcoal rot, a broad-spectrum biotic stress caused by the soil-borne polyphagous fungus Macrophomina phaseolina that reduces seed quality and yield, is the second most important yield-reducing disease in soybean-growing countries worldwide after brown spot (Septoria glycines). This major disease is responsible for 80% of soybean yield loss. The fungus infects the roots and lower stems of soybeans, causing general root rot. Charcoal rot disease appearance looked like roots dipped in charcoal and hence its name. It infects wide range of hosts from sorghum, soybean to cucurbits and various weeds. This was first noticed in the United States in 1949.Every part of the sensitive soybean plant becomes infected as it travels from the roots up through the entire plant. The yield loss in extreme circumstances might be as high as 80%. It is drawing more attention from the breeders in present time because of changing global climatic conditions. The prevalence of recurrent droughts or drought-like circumstances, in particular, increases the susceptibility of soybean to this disease [11]. This issue must be addressed as soon as possible. Therefore, it is anticipated that the creation of KTI-free charcoal rot-resistant soybean cultivars, which are currently unavailable in India, will increase the use of soybean in food items [12].
Marker assisted selection has been widely used for the development of disease and insect resistance in a variety of crops [13]. Nutritional quality traits such as the absence of protease inhibitor or specific fatty acid profile) which cannot be scored visually at the field level, but can be scored effectively by the MAS (Marker assisted selection) technique over morphological traits. Marker-assisted backcross selection (MABS) can make it possible to successfully transfer the ‘ti’ allele responsible for NULL KTI character to superior soybean varieties from genotypes of soybean resistant to charcoal rot disease. In the current study, NRC101 and NRC 127, two superior soybean varieties with a null allele of KTI, were crossed with AMS-MB-5-18 and AMSMB-5-19 to introduce the ti allele via marker-assisted selection. SSR markers were used to determine the purity of the F1 hybrids and the parental polymorphism. To do this, a null allele-specific marker must be utilized in the foreground selection approach to choose the positive plants. The objective of the current investigation was to assess the level of polymorphic SSR loci between the donor and recurrent parents.
MATERIALS AND METHODS
Plant material and fungal isolate
For the investigation’s parental polymorphism survey, four genotypes were examined. The two soybean genotypes resistant to charcoal rot, AMS-MB-5-18 and AMS-MB-5-19, obtained from Dr. Panjabrao Deshmukh Krishi Vidyapeeth, Akola were used as recurrent parents. The other two genotypes, NRC-101 and NRC-127 obtained from Indian Institute for soybean research, Indore provided a null KTI allele and were employed as a donor. The experiment was done in the year 2021-22 and the breeding population was maintained at the field of the Biotechnology Center, Dr. Panjabrao Deshmukh Krishi Vidyapeeth, Akola (Latitude: 200 42’34.45”N, Longitude : 76 0 59’ 53.16” E ) and all molecular work was carried out in the Molecular Breeding laboratory of the Biotechnology Centre, Post Graduate Institute. The experiments were carried out in medium textured loamy soil with a pH of 6.5 to 7.5 and low to medium organic carbon content. For raising the crop in the experiment, the recommended package of practices was followed.
Fungal isolates of Macrophomina phaseolina were collected from the charcoal rot infected soybean plants as well as from soil found in the experimental fields of Regional Research Center, Amravati (Latitude: 20.9320° N, Longitude: 77.7523° E) which is known hot spot for charcoal rot disease). The samples were fetched to the laboratory of Department of Plant Pathology, Dr. PDKV, Akola. The isolates brought were cultured and purified and then re-infected to soybean plants to confirm its identity and virulence. Identified pure culture of M. phaseolina was maintained in controlled conditions and used for in-vitro screening to inoculate the target genotypes.
Performance of parental genotypes in response to Macrophomina phaseolina under in-vitro condition
The lab has already developed methods for screening genotypes against charcoal rot. Out of which pot and invivo showed similar responses. So, the experimental trial for estimating the performance of parental genotypes in response to M. phaseolina was set up using earthen pots. A set of earthen pots was raised to screen the genotypes through artificial inoculation. The experiment was conducted in controlled conditions at the Department of plant pathology, Dr PDKV, Akola. The detailed procedure followed is given below (Figure 1).
Figure 1: Procedure for screening of genotypes in controlled conditions
The scoring of the plants was done using 0-9 scale described by [14], where 0 and 9 indicated immune and highly susceptible, respectively (Table 1).
Confirmation of presence and absence of KTI peptide in selected genotypes
The parental seeds of soybean were biochemically tested through native PAGE for confirming absence of KTI polypeptides. For this purpose, finely ground seed flour (100 mg) was incubated in 1 ml Tris–HCl buffer (pH 8.0) for 30 min. in 700 C and then centrifuged. After this step, prepare 5 % stacking gel and pour and place comb immediately. After solidifying of gel, load the sample prepared having equal volumes of supernatant was added and 5× sample buffer and run at 35 mA for 2 h 30 min. Gels were stained overnight in 0.25 % aqueous solution of coomassie brilliant blue (R-250) in methanol, water and glacial acetic acid (45:45:10) followed by de-staining in methanol, water and glacial acetic acid (45:45:10) solution. Standard trypsin inhibitor protein (21.0 kDa) Thermo Fischer make was loaded in a separate lane for identification of KTI polypeptide.
DNA Extraction from leaf tissue
The chosen genotypes’ young leaves were harvested for the extraction of genomic DNA using the CTAB (Cetyltrimethyl ammonium bromide) technique [15]. The DNA extraction solution (2% CTAB, 100 mM Tris-HCl, 50 mM ethylene diamine tetra acetate (EDTA), and 100 mM NaCl) was applied to crush 15 to 20-day-old soybean leaves [16]. A nanodrop spectrophotometer (IMPLEN make Munchen, Germany) was used to measure optical density for the quantity and purity of the DNA at 260 and 280 nm and agarose gel electrophoresis was also used to assess the quality of DNA. The dilution of genomic DNA was prepared at 40ng/µL for further experiment.
Source of genomic microsatellite sequences
The primer sequence and SSR sites were retrieved from the soybean database [17]. A total of 388 SSR markers from 20 linkage groups (LG) were used, with an average density of one SSR marker per 5 cM, to assess the level of polymorphism between different parental combinations for background selection (Figure 2). Null allele specific markers were used for screening of positive F1 plants.
Figure 2: Histogram showing SSR markers throughout the linkage groups of soybean used under study
SSR Profiling of donor and recurrent parents
The polymerase chain reaction was performed using 325 primers in an Eppendorf thermal cycler. In a 10 μL reaction mixture with 2.0 μL DNA (20 ng/l), 1.0 μL 10X buffer, 1.1 μLMgCl2 , 0.1 μL dNTPs (10 mM), 0.4 μL each of forward and reverse SSR primers, 0.068 μL Taq DNA polymerase (5U/ μL), and 4.932 μL nuclease-free water was subjected to PCR amplification. The template DNAs were amplified in a thermo cycler (Applied Bioscience by Thermo Fischer, Waltham, Massachusetts, United states) using the PCR profile: Initial denaturation at 94°C for 2 minutes. Followed by 32 cycles of denaturation at 94°C for 1 min, primer annealing at 49 to 52 °C for 2 minutes, and primer elongation at 72 °C for 3 minutes each are then performed.
The amplicons’ final elongation was allowed to run its course at 72 °C for 10 min before being put on hold at 4 °C. After completion of amplification, PCR products were stored at –20°C and the amplified products were analyzed by electrophoresis using 8 % Polyacrylamide gel and visualized by silver stain. The silver nitrate staining solution was prepared at 0.18% concentration and the gel was exposed for 10 minutes with continuous shaking and then washed foe 30s in distilled water. Then PAGE gel was treated with a developer solution consisting of 3 % sodium hydroxide with 0.3 % formaldehyde until the bands appeared clearly. Gel documentation unit (Biorad make, Hercules, California, United States) available at Biotechnology Centre, Dr. PDKV., Akola was used for the visualization of resolved PCR products.
Analysis of SSR polymorphism and Detection of informative markers
To assess four genotypes, genotyping by SSR marker was scored based on the presence or absence of the bands was used. The data was then organized as discrete variables into a binary matrix, where score ‘1’ denoted the existence of the band and ‘0’ denoted its absence. This data matrix was subsequently subjected to analysis [18]. The polymorphism information content (PIC) for each chromosomal SSR was calculated using the formula to assess the markers’ in formativeness [19].
PIC=1-(Pi2)
Where Pi is the frequency of the ith plus allele among the group of 4 parental genotypes under investigation, and I is the total number of alleles discovered for the SSR marker. The SSR marker’s discriminating power is estimated by the PIC value [20]. The PIC values are divided into three classes: PIC> 0.7 = highly informative, then 0.25> PIC> 0.7= moderate, and PIC < 0.25 = low. Markers with PIC value above 0.5 are informative in genetic studies and are extremely useful in determining polymorphism rate of specific locus [21]. The evaluation of genetic markers is based on DNA bands of PCR amplification results.
The polymorphism percentage was also calculated between the parents with the formula illustrated below:
Polymorphismpercentage=
The polymorphism percentage of markers will be calculated
Polymorphismpercentageof markers=
The higher the polymorphism percentage the more informative and discriminative it will be.
RESULTS AND DISUSSION
Response of parental genotypes against Macrophomina phaseolina under in vitro The isolates of Macrophomina phaseolina were collected from the charcoal rot infected soybean plants as well as from soil found in the experimental fields of RRC, Amravati. The isolates were cultured, purified and re-infected the soybean plants to confirm its identity and virulence. Pure culture of M. phaseolina was maintained in controlled condition and used for inoculating the target genotypes.
A set of four genotypes was subjected to screening in pots under artificial inoculated conditions as proposed by [14]. The genotypes were selected based on the experiments conducted earlier. The experiment was conducted as per procedure in the net-house of Department of Plant Pathology, PGI, Dr. PDKV Akola (Table 2).
To permit disease development in the potted plants, limited supply of nutrients and water was done. Based on the disease score, 2 genotypes viz., AMS MB 5-19 and AMS MB 5-18 were identified as resistant ones. Similarly, one genotypes NRC 101 was found to be moderately resistant and NRC 127 was found highly susceptible (Figure 3).
Figure 3: Results obtained by screening of the soybean genotypes through screening in pots along with the graphical representation of the survival % of the plants after 45 days.
Therefore, these genotypes were used for developing the mapping population for studying the genetics of charcoal rot resistance.
Confirmation of presence and absence of KTI peptide in selected genotypes
The parental seeds were checked biochemically through NATIVE PAGE for confirming the absence of KTI polypeptides. The samples were loaded in NATIVE PAGE gel and it has been observed that the bands similar to KTI standard were present in genotypes AMS MB 5-18 and AMS MB 5-19 but NRC 101 and NRC 127 does not have. This indicates the absence of KTI allele in donor genotypes NRC 101 and NRC 127 (Figure 4).
Figure 4: Confirmation for presence of Null KTI peptide
Microsatellite marker-based profiling and polymorphism
78 of the 388 SSR markers that were tested were found to be polymorphic collectively for four parental combinations of selected genotypes under investigation among them 46 markers were highly polymorphic (Table 3,4). Three markers Satt144, Satt543, and Satt 301 present on the F, D2, and L linkage groups respectively, were found highly polymorphic in all the parental combinations with the highest PIC value of > 0.75. The representative gel picture showing parental polymorphism of SSR markers is depicted in Figure 5 and monomorphic markers in Figure 6.
Figure 5: Representative images of polymorphic marker for selected genotypes L- Ladder 100kb, P1 - AMS MB 5-18, P2 =AMS-MB-5-19, P3=NRC-101 and P4=NRC-127
Figure 6: Representative images of monomorphic markers for selected genotypes L- Ladder 100kb P1 - AMS MB 5-18, P2 =AMS-MB-5-19, P3=NRC-101 and P4=NRC-127
The level of polymorphism (PIC) is required to select markers that can differentiate between the lines/hybrid used. Markers that produce fewer alleles have a smaller ability to distinguish samples tested. Other investigations have likewise shown low polymorphism [22,23]. In order to produce a polymorphism, it is encouraged to use genetically diverse genotypes in hybridization programs.
The Polymorphic Information Content (pic) value of 31 polymorphic markers for parental combination 1 (AMS MB 5-18 X NRC 101) ranged from 0.11 to 0.75 with 100 polymorphic alleles and average polymorphic percentage of 63.7% among identified polymorphic markers for parental combination 1. The polymorphism percentage for different polymorphic markers ranged from 25 to 100 % and seven markers showed 100 % polymorphism with 0.75 PIC value for parental combination 1. Similarly, 48 markers for parental combination 2 (AMS MB 5-18 X NRC 127) were found polymorphic with PIC value ranging from 0.15 to 0.75. The polymorphism percentage of markers ranged from 20 to 100 % with 220 polymorphic alleles detected out of total 339. 57 polymorphic markers for parental combination 3 (AMS MB 5-19 X NRC 101) has been observed with PIC value ranging from 0.13 to 0.75 and 91 polymorphic alleles out of total 143 alleles with average polymorphism percentage of 63.6 %. In addition, 67 markers have been found polymorphic for parental combination 4 (AMS MB 5-19 X NRC 127) whose PIC value ranged from 0.25 to 0.75 and polymorphism percentage ranged from 25 to 100% with 226 polymorphic alleles (Table 5).
The graph showing the PIC value of all the parental combinations have been shown in Figure 7. The location of each SSR marker found polymorphic has been depicted in Figure 8. The SSR markers which are more informative with greater than 0.30 pic value will be selected for further screening in the breeding program.
Figure 7: PIC value of different polymorphic markers in parental combinations PC1, PC2, PC3 and PC4
Figure 8: Chromosome wise distribution of polymorphic SSR markers of soybean.
Among highly polymorphic SSR markers, three markers namely, Satt 144, Satt 543 and Satt_301 were found polymorphic in all the four parental combinations (Figure 9)
Figure 9: Histogram depicting number of highly polymorphic markers for different parental combinations with highest PIC value
with high PIC value of 0.75 and these markers are located on chromosome 13 (Satt 144) and 17 (Satt 543 and Satt_301).
The average allele frequency for parental combination 1 (AMS MB 5-18 X NRC 101) obtained was 3.84 per marker, for parental combination 2 (AMS MB 5-18 X NRC 127) allelic frequency was 4.3 per marker. Similarly, an allelic frequency of 3.95 per marker was observed for parental combination 3 (AMS MB 5-19 X NRC 101) and parental combination 4 (AMS MB 5-19 X NRC 127) have an allelic frequency of 4.7 per marker.
The number of alleles per marker was found to range from 2 to 5, only. It reflected the poor variability among the soybean genotypes used in the study. The kind of motif (di, tri, tetra, etc.) and the number of repeats of the motif i.e. n=5, 10, 15, etc., in the SSRs have a significant impact on the detection of polymorphism. In the present investigation, it was shown that tri-nucleotide repeat SSR motifs were more polymorphic than other SSR motifs. Ten of the 19 highly polymorphic markers had tri-nucleotide repeat motifs, and seven had di-nucleotide repeat patterns. The SSR motif (TA) 13-58, which has 13–58 repeats, seems to be more variable than other di-nucleotide repeats. In the case of tri-nucleotide repeats, the four cross combinations revealed increased variation in motifs with 5-32 repeats (Table 6, Figure 10).
Figure 10: Histogram depicting number of highly polymorphic markers with different number of motifs in different parental combinations
The markers which have been found to be having highest PIC value and polymorphism percentage in all the four parental combinations have trinucleotide motif (TAT) and (ATT) and only one marker have dinucleotide repeat of (AT) motif with varying number of repeats. In parental combination 1, 4 out of 7 markers have trinucleotide motif, similarly 7 out of 11 in case of parental combination 2, 4 out of 6 in PC3 and 7 out of 12 in PC4 have trinucleotide motif which detected more polymorphism among the four cross combinations. Similar results have already been reported in soybean [24].
The exploitation of informatics markers for hybrid confirmation
There are two possible approaches for adding extra SSR markers to the genome’s sparsely polymorphic regions. The first strategy would include choosing SSR markers from previously existing resources, such as a map released in 2008 by the USDA, and testing them for polymorphism. The second strategy would involve creating new SSRs using the clonal sequences of those locations. In this instance, clone identification and the development of new primers would be made easier by the physical map and the soybean genome’s whole sequence. There have already been reports of similar outcomes in soybean. Null KTI allele-specific marker was used for foreground selection of hybrids (Figure 9) and the selected progenies will be used for further generation advancement. In the future, this information will be helpful for selecting and/or creating SSR primers.
CONCLUSION
Soybean improvement basically revolves around the development of food grade soybean which involves reduction of potent compounds affecting food quality and yield increase by addressing biotic and abiotic stress. This experiment helps in inferring some set of markers polymorphic to specific four soybean genotypes which can be involved in the development of KTI free soybean with charcoal rot disease resistance. Precision molecular breeding in soybeans would benefit from the use of these techniques.
Significance statement
Soybean being a crop with high nutritional value and potential to eradicate malnutrition, improvement is of great significance. Agronomical traits for soybean improvement have been widely studied in soybean and now research is focused to improve its quality. This research experiment was based on addressing the charcoal rot disease caused by M. phaseolina which accounts for 80% yield loss and improving food quality by removing antinutritional factor (KTI) from soybean. The soybean genotypes selected for the research work included two charcoal rot resistant genotypes along with two genotypes with Null KTI allele. The purpose was to introgressed Null KTI allele in the background of charcoal rot resistant genotype. This research article represents first key findings of the experiment in which in vitro screening was done for detecting resistant genotypes and polymorphism detection of all the four genotypes with different parental combinations. The results showed 78 markers out of 388 to be highly polymorphic which will be used for screening in further generation advancement.
ACKNOLWLEDGMENT
The Department of Biotechnology, Government of India, provided financial support for this work as part of the project “Accelerated development of Kunitz trypsin inhibitor free soybean genotype with charcoal rot resistance through marker assisted selection” sponsored by the Accelerated Programme on Crop Improvement. Sanction No. BT/PR30841/AGIII/103/1112/2019
REFERENCES
17. Anonymous. 2023.
24. Kumar B, Talukdar A, Verma K. Diversity analysis in soybean genotypes using SSR markers.In: Book of Abstracts, National Seminar on “Contemporary approaches to crop improvement”, 2011, UAS, Bangalore. 126. 2011.