Next Article in Journal
Genome-Wide Identification and Expression Analysis of Heat Shock Protein 20 (HSP20) Gene Family in Response to High-Temperature Stress in Chickpeas (Cicer arietinum L.)
Next Article in Special Issue
Uncovering Fusarium Species Associated with Fusarium Wilt in Chickpeas (Cicer arietinum L.) and the Identification of Significant Marker–Trait Associations for Resistance in the International Center for Agricultural Research in the Dry Areas’ Chickpea Collection Using SSR Markers
Previous Article in Journal
Effect of γ-Irradiation on the Growth and Yield Response of Three Varieties of Pea (Pisum spp.)
Previous Article in Special Issue
Exploring the Efficacy of Four Essential Oils as Potential Insecticides against Thrips flavus
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Identification of QTLs and Candidate Genes for Red Crown Rot Resistance in Two Recombinant Inbred Line Populations of Soybean [Glycine max (L.) Merr.]

1
Key Laboratory of Biology and Genetics Improvement of Soybean, Ministry of Agriculture, Zhongshan Biological Breeding Laboratory (ZSBBL), National Innovation Platform for Soybean Breeding and Industry-Education Integration, State Key Laboratory of Crop Genetics & Germplasm Enhancement and Utilization, College of Agriculture, Nanjing Agricultural University, Nanjing 210095, China
2
Council for Scientific and Industrial Research-Crops Research Institute (CSIR-CRI), Fumesua, Kumasi P.O. Box 3785, Ghana
3
College of Plant Protection, Nanjing Agricultural University, Nanjing 210095, China
*
Authors to whom correspondence should be addressed.
Agronomy 2024, 14(8), 1693; https://doi.org/10.3390/agronomy14081693
Submission received: 27 June 2024 / Revised: 28 July 2024 / Accepted: 30 July 2024 / Published: 1 August 2024
(This article belongs to the Special Issue Recent Advances in Legume Crop Protection)

Abstract

:
With the rapid emergence and distribution of red crown rot (RCR) across countries, durable sources of resistance against Calonectria ilicicola in soybean [Glycine max (L.) Merrill] is required to control the disease. We employed two RIL populations for the experiment. We identified 15 and 14 QTLs associated with RCR resistance in ZM6 and MN populations, respectively, totaling 29 QTLs. Six and eight QTLs had phenotypic variation above 10% in ZM6 and MN populations, respectively. We identified six (6) “QTL hotspots” for resistance to RCR from the ZM6 and MN RIL populations on chromosomes 1, 7, 10, 11, 13, and 18. Gene annotations, gene ontology enhancement, and RNA sequencing assessment detected 23 genes located within six “QTL Hotspots” as potential candidate genes that could govern RCR resistance in soybeans. Our data will generally assist breeders in rapidly and effectively incorporating RCR resistance into high-yielding accession through marker-assisted selection.

1. Introduction

Soybean [Glycine max (L.) Merr.] is a significant legume crop globally because of its high nutritional values [1]. In China, soybeans are essential legumes for vegetable oil and animal feed are concerned, ranking fourth in its production globally in 2022 (FAOSTAT, https://www.fao.org/faostat, accessed on 31 December 2023). Economic analysis of soybean production shows a high return contributing to the economic and social impact, especially on smallholder farmers [2,3]. Hence, some researchers suggest the need to establish interventions to protect China’s economic, environmental, and sustainable practices for soybean production [4].
Unfortunately, soybean yield parameters are usually confronted with a number of stresses ranging from biotic (pests, diseases, and weeds) to abiotic (drought, salinity, and heavy metals), regulating the quantity and quality of soybean seeds and their nutritional components [5,6]. One of the key constraints impeding soybean production arises from fungal pathogens, which cause 26 different diseases [7]. Particularly in China, Calonectria ilicicola Boedijn and Reitsma are among the economically important pathogens as they infect several plants such as soybean, blueberry, stout camphor tree, groundnuts, alfalfa, sassafras, and ginger [8,9,10]. Now, soybean red crown rot (RCR) disease caused by C. ilicicola is a threat to major soybean regions, especially in southern Asia. Soybean RCR negatively impacts the quantity of seed in a pod, its seed weight, and quality [11,12,13], thus limiting yield by 25–30% [14,15,16], and it is projected to reach 50–100% disease incidence on susceptible cultivars [11,14,17]. Morphologically, strains of C. ilicicola are characterized by anamorphic structures (macroconidia, vesicles, conidiophores, and conidia) and teleomorphic structures (perithecia, asci, and ascospores) [18,19]. C. ilicicola has also been identified via internal transcribed spacer (ITS) regions, calmodulin, histone3, and translation elongation factor 1-α [19,20]. Recently, novel Calonectria species and new hosts have been recorded in Southern China [21]. The fungus C. ilicicola infects the roots of soybeans, resulting in root rot advancing to pencil-like roots, damping off of younger seedlings, and untimely defoliation [22,23,24]. Currently, neither an effective fungicide nor a resistant cultivar controls RCR; however, agronomic practices are adopted to manage the disease [15,25,26]. The generation of RCR-resistant varieties continues as the foundation for disease management.
Research on soybean diseases is of growing interest owing to its consequences. Hence, several researchers have investigated the inherent features of disease traits via quantitative trait loci (QTL) through linkage analysis and quantitative trait nucleotides (QTNs) via genome-wide association studies (GWAS) [27,28]. The disease traits in soybeans are associated with the growing season, isolates and their diversity, and the soybean cultivar [7,29]. Thus, disease traits are complex and are controlled by several genes [30,31]. Therefore, it is crucial to identify key genes that could provide resistance to soybean crops. Molecular breeding provides an inclusive comprehension of the genetic basis underlying disease traits. Previous researchers have identified QTLs/QTNs utilizing either GWAS or linkage mapping analysis to unravel the genetic basis of soybean disease traits [31,32]. The success of QTL/QTN identification largely depends on the mapping strategies used. However, few soybean germplasm accessions show RCR resistance [25]. Hence, the application of GWAS becomes comparatively difficult as rare variants could be omitted due to the minor allele frequency (MAF) cut-off, leading to skewed phenotypic distribution toward the susceptibility genotypes [33]. On the contrary, biparental population RCR-resistance sources could be studied by utilizing the genetic background of the resistant parent. Linkage mapping analysis has proven fruitful in detecting markers linked with various soybean disease traits. For instance, linkage mapping has successfully been applied in soybean diseases such as rust [34,35], bacterial leaf pustule [32,36], mold white [37,38], Phytophthora root rot [39,40,41], and frog eye leafspot [42,43], resulting in detecting QTLs. Conversely, bi-parental mapping populations have yet to be employed in the search for QTL related to RCR resistance in soybeans. Therefore, the current investigation used a high-density linkage map of two RIL populations, namely ZM6 and MN, and evaluated two different times in an environment to identify QTLs and mine for potential candidate genes for RCR resistance in soybeans. We sought to identify QTLs and potential genes governing regulating RCR resistance in soybeans. This study deepens our understanding of the genetics of RCR resistance in soybeans, thus contributing to the identification, development, and application of markers related to RCR resistance for marker-assisted selection toward developing resistant cultivars.

2. Materials and Methods

2.1. Plant Material

After preliminary RCR resistance evaluation of parents of six recombinant inbred line (RIL) populations (LM6, ZM6, M6T, M6T, MN, ZN), the parents of ZM6 and MN RIL populations showed significant resistance differences, so these two populations were selected for mapping RCR resistance QTL. Both RIL populations were developed through the single-seed descent (SSD) method by advancing to the seventh generation. The ZM6, consisting of 122 lines, was developed by crossing a disease-sensitive variety, Zhengyang, and a disease-tolerant cultivar, Meng8206 to C. ilicicola. Also, MN, consisting of 98 lines, was developed using the soybean cvs. M8180 and NN1138-2, which are resistant and susceptible to C. ilicicola strain Y62, respectively. All the RIL population accessions were sourced from the Soybean Germplasm Gene Bank, located at the National Centre for Soybean Improvement (Ministry of Agriculture), Nanjing Agricultural University, Nanjing, China.

2.2. Pathogen Culture and Inoculation

The C. ilicicola strain Y62 was provided by the College of Plant Protection, Nanjing Agricultural University, Nanjing, China. The fungi were maintained on V8 media plates (90 mm) at 26 °C for short-term use by subculturing. Pathogen inoculation was carried out following the protocol of Jiang et al. (2020) [25].

2.3. Growth Conditions, Experimental Conditions, and Design

The planting media consisted of a mixture of vermiculite and nutritive soil in equal volumes (1:1, v/v). The media was autoclaved. After cooling at ambient temperature for two days, the media was combined with an inoculum–soil mixture to achieve a concentration of 2% (w/v) and then filled into the plastic pot as the growing media as described by Jiang et al. (2020) [25]. The seeds for planting consist of the two RIL populations, along with their parents. Each pot was planted with ten seeds, and the surface was then covered with a two-millimeter layer of the growing medium. Pots were positioned within a container, and water was directed to the pots through their drainage openings to make sure they were completely saturated in a greenhouse maintained at a temperature of 26 °C and a relative humidity of 50%. Water was supplied to the container to maintain soil moisture throughout the assay. To mitigate against potential location-related impacts, the plastic pots in the greenhouse were changed every two days. The study used a completely randomized design (CRD) with three replications.

2.4. Phenotypic Data

Phenotypic data were taken on the emergence rate (ER), survival rate (SR), and disease severity (DS). The data on the traits were taken twice. The first data set was taken between February and April 2023, termed as the first (1st) data set. The second time data set was taken between June and August 2023 as the second (2nd) data set.

2.4.1. Assessment of the Rate of Emergence and Survival for Resistance to Calonectria ilicicola

The soybean accessions were scored for emergence rate on the 5th day after planting (DAP). The emergence rate was expressed as the total number of seeds that emerged out of the total number of seeds planted expressed in percentage. The survival rate was taken on the 12th DAP. Survival rate was calculated as the total number of plants alive out of the total number of plants that emerged expressed in percentage.

2.4.2. RIL Populations Lines Response to Calonectria ilicicola Infection

The seedlings were evaluated for C. ilicicola resistance on the 14 DAP according to Jiang et al. (2020) [25]. The seedlings were gently removed from the pot, and the roots were washed to ensure that the roots were free from soil for visual examination. The resistance level of each RIL was determined using DS in addition to their SR and ER traits. A threshold of SR > 0.90 and ER > 0.85 was used as a standard to identify accessions with resistance. The SR helped to determine the mortality rate of the seedling, which was incorporated into the disease resistance rating scale. The DS employs a numerical scale ranging from 0 to 5, as indicated in Table 1. These were employed to categorize RILs according to their responses to C. ilicicola infection. The average of the three repeats of each RIL was further used to classify them using y ≤ 1.5, 1.5 ≥ y ≤ 2.5, 2.5 ≥ y ≤ 3.5, 3.5 ≥ y ≤ 4.5, and y > 4.5 as highly resistant, resistant, moderately susceptible, susceptible, and highly susceptible to C. ilicicola, respectively.

2.4.3. Analysis of Phenotypic Data

The GenStat software, version 12 (VSN International Ltd., Hemel Hempstead, UK), was utilized to generate descriptive statistics for the ER, SR, and DS traits in both RIL populations as well as their parents. These statistics included the mean, range, degree of variation (CV, %), skewness, and kurtosis. A Pearson correlational study was conducted to examine the relationship between ER, SR, and DS. The results were shown using the Corrplot package in R (version 0.84, Auckland, New Zealand). The level of significance was set at p < 0.05 [44].

2.5. SNP Genotyping and Linkage Map Construction for ZM6 and MN RIL Populations

We extracted the genomic DNA using young and fresh leaves of the entire populations of ZM6 and MN RILs in addition to the parents, following the standard established by [45]. The DNA library construction, sequencing, SNP acquisition, and integration of Bin/SLAF (specific-locus amplified fragment) markers were performed for the ZM6 and MN populations. For the ZM6 population, bin markers were used for genotyping [46,47]. Also, the MN population was genotyped using SLAF markers [46,47]. The high genetic maps of the ZM6 and MN populations comprised 2601 bin markers and 2062 SLAF markers, respectively. The combined length of the ZM6 and MN maps was 2630.22 cM and 2054.5 cM, having a mean distance of 1.02 cM and 1.36 cM between neighboring markers, respectively (Table S1). The mean length of each linkage group was 131.51 cM for the ZM6 linkage map and 102.73 cM for the MN linkage map. The average number of markers per chromosome was 130.05 and 103.1 from the ZM6 and MN populations, respectively (Table S1).

2.6. QTL Mapping for Emergence Rate, Survival Rate, and Disease Severity

The WinQTLCart (version 2.5, Raleigh, NC, USA) application was used to locate QTLs via composite interval mapping (CIM) [48]. The window size, working speed, and control were configured to be 10 cM, 1 cM, and less than 5 cM, respectively [49]. The LOD thresholds for identifying significant QTLs were established at 2.5, using 1000 permutations for each test and a type I error rate of 5% [50]. The precise physical positions of individual QTLs on all chromosomes were determined using the MapChart (version 2.1, Wageningen, Netherlands) tool [51]. Names were assigned to the QTLs by adhering to the accepted nomenclature system [52]. For instance, for the QTL denoted as qER-7-1zm6, the letter q stands for QTL preceded every QTL name, which was followed by a two-letter descriptor of the ER trait, 7 refers to the specific chromosome on which the QTL was identified, 1 displays the number of each QTL trait on the chromosome, and zm6 represents the ZM6-RIL population where the QTL existed.

2.7. Mining of Candidate Genes for Major QTLs

We obtained the complete genomic data from the main “QTL hotspots” by utilizing the internet-based resource libraries from Phytozome (http://phytozome.jgi.doe.gov (accessed on 13 January 2024)) and SoyBase (http://www.soybase.org (accessed on 13 January 2024)), which also assisted in candidate gene identification along with the data obtained from previously published publications. We used the Gene Ontology (GO) online platform to analyze the identified potential genes. GO enrichment analysis was performed for all the QTL hotspots using agriGo (version 2.0, Beijing, China) [53].
The available RNA sequencing dataset by Kobayashi et al. (2022) [54] and our unpublished data were applied to investigate the transcription of identified candidate genes in different soybean tissues and phases of development. A heat map was generated using TBtools version 1.6 to visually represent the changes in gene expression trends for the anticipated candidate genes [55].

3. Results

3.1. Variability Characteristics of the Mapping Traits for the RIL Populations

The use of parents with varied phenotypes is necessary for the creation of mapping populations, and this is crucial for maximizing the effectiveness of QTL discovery [56]. For this investigation, two mapping RIL populations were utilized: ZM6 RIL, which originated from the cross between Zhengyang (♀) and Meng8206 (♂), and MN, which originated from the cross between M8180 (♀) and NN1138-2 (♂) (Table 2). The ZM6 and MN RIL populations display the mean, range (minimum and maximum value), standard deviation, skewness, and kurtosis of the ER, SR, and DS for the two RIL populations (ZM6 and MN) and their parents (Table 2). These measurements were examined for two conservative times (first and second). The parental genotypes exhibited relatively large variations for ER, SR, and DS (Table 2). The RILs mostly showed a transgressive segregation. The skewness and kurtosis were mostly less than 1, which is typical for quantitative attributes, implying that the populations are suitable for conducting QTL mapping (Table 2). Furthermore, the measured variables (ER, SR, and DS) exhibited a distinct range of distribution across both RIL populations (Figure 1). Based on the three traits screened, there were significant positive correlation coefficients between SR and ER in both mapping populations, while significant negative correlation coefficients existed between DS against ER and SR (Figure 2).

3.2. QTL Detected in the ZM6 Population for Red Crown Rot Resistance in Soybean

The composite interval mapping (CIM) strategy was employed using the 2601 bin markers obtained among the ZM6 RIL population and was polymorphic between the parents. A total of 15 QTL for the three traits in ZM6 RIL population were located on eight different chromosomes (Chr06, Chr07, Chr08, Chr10, Chr11, Chr13, Chr17, and Chr18), of which 5, 6, and 4 were for ER, SR, and DS, respectively, and each QTL had phenotypic variance explained (PVE) ranging between 8.40 and 14.87% (Table 3; Figure 3). The 5 QTL for ER consisted of qER-7-1zm6 (Chr07), qER-8-1zm6 (Chr08), qER-10-1zm6 and qER-10-2zm6 (Chr10), and (qER-11-1zm6) (Chr11) with their PVE ranging from 8.40 to 10.70%. With the exception of the 2 QTL on Chr10 (qER-10-1zm6 and qER-10-2zm6), which inherited its alleles from the male parent, Zhengyang, the remaining three QTLs (qER-7-1zm6, qER-8-1zm6, and qER-11-1zm6) had their alleles from the female parent, Meng8206 (Table 3).
Six QTLs were mapped for SR across five chromosomes: Chr06 (qSR-6-1zm6), Chr07 (qSR-7-1zm6, and qSR-7-2zm6), Chr10 (qSR-10-1zm6), Chr11 (qSR-11-1zm6), and Chr17 (qSR-17-1zm6) (Table 3). Out of these, qSR-7-2zm6 had the most prominent QTL with a LOD value of 5.19, accounting for 14.87% of the phenotypic variation, which was followed by qSR-7-1zm6 with a LOD score of 4.55 and a difference in phenotype of 12.71% (Table 3). Four QTLs influencing DS were detected across three chromosomes (Chr11, Chr13, and Chr18), accounting for an average of 10.04% of PVE (Table 3). The Chr011 and Chr13 harbored one QTL each (qDS-11-1zm6 and qDS-13-1zm6, respectively), while Chr18 had two QTL (qDS-12-1zm6 and qDS-12-2zm6). Figure 4 shows the location of QTLs on the linkage genetic map of the ZM6 RIL population.

3.3. QTL Detected in MN Population for Red Crown Rot Resistance in Soybean

The MN RIL population detected 14 QTLs across eight chromosomes (Chr01, Chr02, Chr04, Chr08, Chr10, Chr13, Chr15, and Chr17) mapped for ER, SR, and DS comprising 5, 4, and 5 QTLs, respectively, with their phenotypic variation ranging from 8.69% to 25.24% (Figure 5; Table 4). For the ER, the five QTLs, namely, qER-1-1mn, qER-10-1mn, qER-8-1mn, qER-8-2mn, and qER-15-1mn had their LOD values in the range of 2.63 to 7.61 (Table 4). The QTLs for the SR were located on Chr01 (qSR-1-1mn, and qSR-1-2mn), Chr02 (qSR-2-1mn), and Chr17 (qSR-17-1mn) and their PVE ranged from 8.68% to 10.7% and their alleles were inherited from NN1138-2 (♂) except qSR-2-1mn (Table 4). Five QTLs, namely, qDS-1-1mn, qDS-1-2mn, qDS-1-3mn, qDS-4-1mn, and qDS-13-1mn, were mapped for DS with PVE ranging from 10.51% to 12.28% (Table 4). Two QTLs (qER-1-1mn and qSR-1-2mn) were detected on Chr01 across the second environment, accounting for 33.93% of the phenotypic variance, suggesting it is relative stability (Table 4). These two QTLs colocalized within the genomic interval 37,147,226–38,075,381 bp. Figure 6 shows the location of QTLs on the genetic linkage map of the MN RIL population.

3.4. Colocalization of QTLs in the QTL Hotspot

The QTL hotspot refers to the concentrated area on a chromosome that hosts many QTLs related to distinct traits. The stability of QTL across mapping populations and environments is essential for their use in practical plant breeding [57]. Hence, our stable QTLs were detected either linked to at least two of the traits (ER, SR, or DS) across an environment (first and second time of screening) within a population or across a population (found in ZM6 and MN). The QTLs qER-1-1mn and qSR-1-2mn were detected on Chr01 for two conservative times during the screening; hence, the stable QTL was termed as Hotspot A (Table 5). The QTLs were stable and are colocalized within the genomic interval of 8,823,531–44,479,895 bp. The mapping results of the ZM6 population found two pairs of QTLs, namely, qSR-7-2zm6 colocalized within the genomic region of 14,134,797–15,903,280 bp across the two environments and renamed as QTL Hotspot B (Table 5). Also, qER-10-1zm6 and qSR-10-1zm6 colocalized within genomic interval 1,603,735–2,732,880 bp, termed as QTL Hotspot C. Hotspot D refers to QTL qER-11-1zm6 in addition to qSR-11-1zm6 colocalized within the genomic interval of 14,962,695–16,816,800. Interestingly, a stable QTL qDS-13-1zm6 and qDS-13-1mn colocalized within genomic interval 68,8713-5,375,644 bp was detected across the MN and ZM6 population and named QTL Hotspot E (Table 3, Table 4 and Table 5). The QTLs qSR-7-2zm6 and qDs-18-2zm6 were constantly found in ZM6 in the first and second screenings, indicating that they were stable; they were thereafter classified as QTL hotspots B and F, respectively (Table 5).

3.5. Candidate Gene Mining within Major “QTL Hotspots”

One of the goals of QTL mapping research is to identify possible candidate genes. This is made possible by the presence of complete genetic information and detailed gene descriptions. Therefore, model genes within the QTL hotspot and their corresponding gene annotations were obtained via SoyBase (https://www.soybase.org, retrieved on 26 December 2023) using soybean genome version 2. A total of 1279 gene models were found in the genomic region for each of the significant “QTL hotspots,” of which 529, 116, 118, 122, 274, and 120 genes were found within hotspots of A, B, C, D, E, and F, respectively (Table S3).
The gene ontology (GO) analysis of enrichment was conducted using the web-based toolbox agriGO (version 2.0, Beijing, China) [53]. It enabled the depiction of the main categorizations of biological process (BP), molecular function (MF), and cellular component (CC) (Figures S1–S3). We observed the highest number of genes in CC followed by BP, and MF was the least for all the QTL hotspots, as expected for QTL hotspot C, for which the BP was higher than the CC, and the MF recorded the lowest number of genes (Figures S1–S3). Analysis of the GO showed that the majority of the genes associated with the six “QTL hotspots” in CC engage in the cell, cell part, and organelle processes; the BP partake largely in processes involving cellular, metabolic, and single-organism processes; and the MF genes mainly engage in the activities of catalytic, binding, and transport (Figures S1–S3).
Using gene annotation descriptions, existing research findings, and GO enrichment examination, our study identified 32, 12, 15, 12, 25, and 15 potential genes from “QTL Hotspot A”, “QTL Hotspot B”, “QTL Hotspot C”, “QTL Hotspot D”, “QTL Hotspot E”, and “QTL Hotspot F”, totaling 111 predicted genes (Table S4). These genes play a role in controlling disease infection development, including developing resistance and innate immune response. They are involved in various processes such as signal transduction; biosynthetic processes involving salicylic acid, lignin, and flavonoid; phenylpropanoid metabolic processes; and cell death. To enhance the accuracy of the earlier-predicted list of candidate genes, we used an available online database by Kobayashi et al. (2022) [54] (Table S5) and that of our RNA-Seq data (unpublished, though is available upon request). According to the RNA-seq data, 23 genes from a total of 111 projected candidate genes exhibited a significantly greater level of transcription log2 (FC) of ≥1 in the root, stem, R2-1d, R2-2d, S-1d, and S-2d (Figure S4A–F; Table 6). These include five genes (Glyma.01G127100, Glyma.01G126600, Glyma.01G112300, Glyma.01G127200, and Glyma.01G127700); two genes (Glyma.07G133900, and Glyma.07G134100); five genes (Glyma.10G019900, Glyma.10G023300, Glyma.10G023400, Glyma.10G028200, and Glyma.10G025700); three genes (Glyma.11G181200, Glyma.11G192300, and Glyma.11G193600); three genes (Glyma.13G066100, Glyma.13G069200, and Glyma.13G076200); and five genes (Glyma.18G287000, Glyma.18G289100, Glyma.18G289600, Glyma.18G293200, and Glyma.18G293300) from QTL Hotspots A, B, C, D, E, and F, respectively (Figure S4; Table 6). Therefore, the 23 genes are potential candidates for regulating C. ilicicola resistance in soybeans. Nevertheless, additional functional validation is required to verify their specific functions in modulating resistance to soybean RCR disease.

4. Discussion

The objective of this study was to employ a high level of intraspecific mapping of ZM6 and MN populations, analyzed for two different periods, to find stable QTLs (QTL Hotspots) and to determine probable candidate genes associated with soybean’s resistance to C. ilicicola. We found six QTL hotspots and predicted 23 potential candidate genes (PCGs).

4.1. Examining the Resistance to RCR in Two RIL Populations

Due to its fast spreading and severity in major production regions, soybean pathologists and breeders have aimed to identify and incorporate RCR resistance into high-quality soybean varieties. Earlier works have developed screening protocols for evaluating soybean germplasm for resistant accession identification [25,58,59,60]. There are a few sources of incomplete resistance, primarily found in wild relatives or non-elite germplasm. This requires a process of crossing and careful selection, which is made more difficult by environmental factors that affect the pathogen’s lifestyle. The ANOVA analysis showed a substantial variation across all the RILs in the ZM6 and MN populations for all three traits, with the environment playing a substantial role (Table S2). The research found that the paternal accessions Zhengyang and NN1138-2 showed a significant resistance level, while Meng8206 and M8180 revealed some level of susceptibility. The disease severity index analysis revealed that the resistant accessions Zhengyang and NN1138-2 potentially harbor distinct resistance genes. The high correlation coefficients among the traits used for screening RCR resistance imply an independent and interdependent relationship among the traits, suggesting some common loci may be responsible for at least two of the traits studied (Figure 2).

4.2. QTL Detected for Red Crown Rot Resistance in ZM6 and MN Population

The availability of genetic markers is crucial to assist QTL identification. Furthermore, combining QTLs that resist C. ilicicola infection within a cultivar is an additional advantage of identifying QTLs connected with soybean RCR resistance. Although some progress has been made in establishing resistance to RCR, data regarding DNA markers associated with resistance to RCR are scarce. Linkage mapping has the potential to facilitate mutual verification and enhance the accuracy of produced results, promoting its application in several studies involving disease resistance gene identification in soybeans using RIL populations [34]. The high-density genetic maps for ZM6 were genotyped by RAD-seq technology [32], and the MN population was also genotyped by SLAF technology [61]. The markers on both linkage maps, ZM6 and MN, were incorporated into all the 20 linkage groups and covered a length of 2630.22 cM and 2054.50 cM, respectively. In ZM6, the average distance among adjacent markers was 1.01 cM, while in MN, it was 1.00 cM. By combining these high-density genetic maps, there is potential to accurately locate QTL that has significant and subtle effects on resistance to RCR in the mapping populations.
The capacity to detect QTL determines the extent of the genomic region, which is decided by the gene variation between parental individuals and is influenced by population structure and size [62,63]. The two RIL populations identified 15 and 14 QTLs associated with RCR resistance in ZM6 and MN populations, respectively, totaling 29 QTLs. For the ZM6, the most substantial effect was on chromosome 7, qSR-7-2zm6, with LOD 5.19, while for MN, the highest was on chromosome 1, qER-1-1mn with a LOD of 7.61 (Table 3 and Table 4). This confirms that the phenotypic data gathered in different settings exhibit statistical differences [64]. The identification of QTL for resistance to RCR provides evidence that quantitative traits are impacted by multiple QTLs, each with relatively small individual effects. Hence, we hypothesize that resistance to RCR is a complex genetic characteristic regulated by several resistance gene loci. For instance, in soybeans, several QTLs are identified for S. sclerotiorum and P. sojae resistance, confirming that their genetic architecture is influenced by several genes [7,65,66].
The QTLs qER-7-1zm6, qER-10-2zm6, qER-11-1zm6, qSR-7-1zm6, qSR-7-2zm6, and qSR-11-1zm6 in the ZM6 accounted for phenotypic variation above 10% and were on Chr07, Chr10, or Chr11. Therefore, the desired alleles of the genetic markers can significantly decrease the disease traits of RIL progenies and alter the response characteristics of C. ilicicola. Also, the MN population had eight QTLs with a PVE of more than 10%, of which qER-1-1mn had its PVE of 25.24% on Chr01. Thus, our findings demonstrated the dependability of QTL mapping. Furthermore, QTLs could be utilized as main targets for discovering candidate genes and implementing marker-assisted selection in future studies. The results prove that RIL has segregated homozygous alternative alleles and can detect QTL due to the presence of half alleles in the RIL [67].

4.3. Identification of QTL Hotspot for Resistance to RCR in Soybean

The durability of the QTL is essential for its application in a breeding strategy. This study additionally showed the localization of QTLs on chromosomes for many variables (ER, SR, and DS) associated with RCR resistance. For instance, in the two environments, we confirmed the presence of six (6) QTL hotspots for resistance to RCR from the ZM6 and MN RIL populations on chromosomes 1, 7, 10, 11, 13, and 18 (Table 5). Hotspot E consisted of qDS-13-1zm6 and qDS-13-1mn, which were obtained from both populations colocalized in the same physical interval. Our investigation found that the qSR-7-2zm6 and qDS-18-2zm6 genes in the ZM6 population were localized in the same area on chromosome 7 for two consecutive years (Table 5). These findings indicate that there may be a shared resistance mechanism to several RCRs in soybeans, and a specific gene on chromosome 7 could be responsible for providing resistance to RCR. We have a firm conviction that these QTLs have the potential to be significant QTLs for either gene cloning or marker-assisted selection.
Determining the precise gene(s) accountable for resistance loci is an essential stage in the ongoing endeavor to examine the molecular processes and biological basis of quantitative resistance to diseases. The gene ontology (GO) study revealed that the majority of genes in the six hotspots were mainly involved in cellular functions such as the cell, cell component, organelle, cellular process, metabolism process, catalytic activity, binding, and transporter activity. These substances are crucial for the growth and development of plants, and multiple studies have verified their role in safeguarding plants from diseases caused by pathogens. To illustrate, the cells, cell parts, and organelles of plants participate in the hypersensitive response in plants due to gene-for-gene resistance. For instance, the cells located around the area where the fungus enters experience a process of programmed cell death (PCD) [68,69], cellular processes, metabolic processes [70] and catalytic activity, binding, and transporter activity [71,72,73] and are involved in plant immunity in order to impede the progression of the disease.
At the physical position of the six “QTL Hotspots”, 1279 gene models were determined. Among these, 111 genes were identified as PCGs based on the study of GO enrichment, gene function, and existing literature. The actions of these 111 genes engage in activating the plant’s immune response during a pathogen attack. Also, some of the genes activate secondary metabolite processes (lignin biosynthetic, flavonoid biosynthetic, and phenylpropanoid metabolic) as well as signal transduction. Of the 111 predicted PCGs, 23 are noted to provide resistance against C. ilicicola invasion in soybeans based on their gene expression data obtained from RNA-seq. Among the 23 PCGs, seven genes, namely, Glyma.01G112300, Glyma.01G126600, Glyma.01G127100, Glyma.01G127200, Glyma.01G127700, Glyma.13G076200, and Glyma.18G287000, are reported to offer resistance to crops upon pathogen infection either by activating defense response, signal transduction, or leucine-rich repeat responses/processes. For instance, in soybean, a gene annotated as defense response governs tolerance to Phytophthora root rot disease [74], while leucine-rich repeat has the highest number of documented disease resistance genes [75]. Similarly, Glyma.07G133900 and Glyma.07G134100 encode multicopper oxidase. In wheat, ascorbate oxidase and skewed5-similar proteins resulting from multicopper oxidase partake in the plant’s immunity, and their silencing promotes the plant’s immunity to Verticillium wilt [76]. The soybean GmHSP40 results in cell death and enhances its immunity to mosaic virus, as its silencing results in making soybean more sensitive to the virus [77]. Similar roles are confirmed in pepper by overexpressing CaHSP70a [78]. Hence, Glyma.18G289100 and Glyma.18G289600 encoding the Hsp70 protein could possibly function to regulate RCR resistance.
We anticipated that the 23 predicted candidate genes may have a potential role in regulating RCR resistance in soybeans. However, additional testing and verification are necessary to demonstrate their specific involvement in soybean resistance to RCR, as well as their possible application in disease breeding programs. Also, the identified loci and QTL hotspot will aid in validating C. ilicicola resistance in soybeans through systematic breeding. Additionally, these genomic areas can serve as targets for enhancing our understanding of the RCR mechanism and improving soybean resistance to C. ilicicola. The identified harbor QTLs demonstrating promise will be employed for precise mapping and molecular cloning of crucial loci in the future. These regions can then be utilized to improve soybean plant immunity to C. ilicicola invasion.

5. Conclusions

The study examined the genetic architecture of soybean resistance to C. ilicicola invasion. A total of 29 QTLs were identified via linkage mapping using two different RIL populations. All the detected QTLs were grouped into six significant “QTL hotspots” and represent the major and consistent genomic areas that control the inheritance of soybean resistance to RCR. Among them, QTLs qER-1-1mn and qSR-7-2zm6 accounted for 25.24% and 14.87% of the overall phenotypic variation, respectively. A common QTL qDS-13-1zm6 and qDS-13-1mn were obtained from both populations colocalized in the same physical interval. We spotted 23 genes likely to regulate RCR resistance based on the six genomic areas known as “QTL hotspots”. The identified genes serve as a ground for future research in developing genetic resources for enhancing soybean resistance. Similarly, additional functional studies are required for the validation and cloning of functional genes of the proposed candidate genes to determine their exact functions in regulating RCR resistance. These findings establish a crucial basis for generating fungal-resistant QTL for RCR in soybeans through cloning.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/agronomy14081693/s1, Table S1: Distribution of bin and SLAF markers mapped on soybean chromosomes/linkage groups in ZM6 and MN Populations. Table S2: Combined analysis of variance (ANOVA) for RCR evaluation trait (ER, SR, and DS) in ZM6 & MN RIL Population across two time screening (1st and 2nd time screening). Table S3: Model genes within hotspot A, hotspot B, hotspot C, hotspot D and hotspot E regions in both RIL Populations (ZM6 & MN). Table S4: Predicted candidate genes within hotspot A, hotspot B, hotspot C, hotspot D and hotspot E regions in both RIL Populations (ZM6 & MN) based on known functional annotation. Table S5: DEGs following infection by C. ilicicola Kobayashi et al., 2022 [54]. Figure S1: AgriGO annotation Palatino Linotypefor the QTL Hotspot A and B. (A) QTL Hotspot A, (B) QTL Hotspot B. Figure S2: AgriGO annotation information for the QTL Hotspot C and D. (C) QTL Hotspot C, (D) QTL Hotspot D. Figure S3: AgriGO annotation information for the QTL Hotspot E and F. (E) QTL Hotspot E, and (F) QTL Hotspot F. Figure S4: Heat map displaying the expression patterns of 23 selected genes in various soybean tissues from six QTL hotspots.

Author Contributions

A.A.-B., T.Z. and J.F. conceived and conducted the experiments; A.A., J.L. N.G., and C.C. performed part of the data collection; C.Z. and S.J. performed part of the data analysis work. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the National Natural Science Foundation of China (Grant Nos. 31871646, 32171965), the Core Technology Development for Breeding Program of Jiangsu Province (JBGS-2021-014), the Program of Collaborative Innovation Center for Modern Crop Production co-sponsored by Province and Ministry (CIC-MCP).

Data Availability Statement

The original contributions presented in the study are included in the article/Supplementary Materials, further inquiries can be directed to the corresponding author.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Henchion, M.; Hayes, M.; Mullen, A.M.; Fenelon, M.; Tiwari, B. Future protein supply and demand: Strategies and factors influencing a sustainable equilibrium. Foods 2017, 6, 53. [Google Scholar] [CrossRef] [PubMed]
  2. Cattelan, A.J.; Dall’Agnol, A. The rapid soybean growth in Brazil. OCL 2018, 25, D102. [Google Scholar] [CrossRef]
  3. Ayalew, B.; Bekele, A.; Mazengia, Y. Analysis of cost and return of soybean production under small holder farmers in Pawe District, North Western Ethiopia. J. Nat. Sci. Res. 2018, 8, 28–34. [Google Scholar]
  4. Liao, Z.; Pei, S.; Bai, Z.; Lai, Z.; Wen, L.; Zhang, F.; Li, Z.; Fan, J. Economic Evaluation and Risk Premium Estimation of Rainfed Soybean under Various Planting Practices in a Semi-Humid Drought-Prone Region of Northwest China. Agronomy 2023, 13, 2840. [Google Scholar] [CrossRef]
  5. Sabagh, A.E.; Hossain, A.; Islam, M.S.; Barutçular, C.; Ratnasekera, D.; Kumar, N.; Meena, R.S.; Gharib, H.S.; Saneoka, H.; da Silva, J.A.T. Sustainable soybean production and abiotic stress management in saline environments: A critical review. Aust. J. Crop Sci. 2019, 13, 228–236. [Google Scholar] [CrossRef]
  6. Miransari, M. Environmental Stresses in Soybean Production: Soybean Production Volume 2; Academic Press: Cambridge, MA, USA, 2016; Volume 2. [Google Scholar]
  7. Lin, F.; Chhapekar, S.S.; Vieira, C.C.; Da Silva, M.P.; Rojas, A.; Lee, D.; Liu, N.; Pardo, E.M.; Lee, Y.-C.; Dong, Z. Breeding for disease resistance in soybean: A global perspective. Theor. Appl. Genet. 2022, 135, 3773–3872. [Google Scholar] [CrossRef]
  8. Gai, Y.; Deng, Q.; Chen, X.; Guan, M.; Xiao, X.; Xu, D.; Deng, M.; Pan, R. Phylogenetic diversity of Calonectria ilicicola causing Cylindrocladium black rot of peanut and red crown rot of soybean in southern China. J. Gen. Plant Pathol. 2017, 83, 273–282. [Google Scholar] [CrossRef]
  9. Fei, N.; Qi, Y.; Meng, T.; Fu, J.; Yan, X. First report of root rot caused by Calonectria ilicicola on blueberry in Yunnan Province, China. Plant Dis. 2018, 102, 1036. [Google Scholar] [CrossRef]
  10. Zhang, Q.; Zhou, D.; Jiang, W.; Zhu, H.; Deng, S.; Wei, L. First report of soft rot of ginger caused by Calonectria ilicicola in Guangxi Province, China. Plant Dis. 2020, 104, 993. [Google Scholar] [CrossRef]
  11. Roy, K.; McLean, K.; Lawrence, G.; Patel, M.; Moore, W. First report of red crown rot on soybeans in Mississippi. Plant Dis. 1989, 73, 273. [Google Scholar] [CrossRef]
  12. Ochi, S.; Mimuro, G.; Kishi, S. Effect of red crown rot of soybean on occurrence of wrinkled seeds. J. Gen. Plant Pathol. 2022, 88, 232–238. [Google Scholar] [CrossRef]
  13. Yamamoto, S.; Nomoto, S.; Hashimoto, N.; Maki, M.; Hongo, C.; Shiraiwa, T.; Homma, K. Monitoring spatial and time-series variations in red crown rot damage of soybean in farmer fields based on UAV remote sensing. Plant Prod. Sci. 2023, 26, 36–47. [Google Scholar] [CrossRef]
  14. Kleczewski, N.; Plewa, D.; Kangas, C.; Phillippi, E.; Kleczewski, V. First report of red crown rot of soybeans caused by Calonectria ilicicola (anamorph: Cylindrocladium parasiticum) in Illinois. Plant Dis. 2019, 103, 1777. [Google Scholar] [CrossRef]
  15. Akamatsu, H.; Fujii, N.; Saito, T.; Sayama, A.; Matsuda, H.; Kato, M.; Kowada, R.; Yasuta, Y.; Igarashi, Y.; Komori, H. Factors affecting red crown rot caused by Calonectria ilicicola in soybean cultivation. J. Gen. Plant Pathol. 2020, 86, 363–375. [Google Scholar] [CrossRef]
  16. Win, K.T.; Kobayashi, M.; Tanaka, F.; Takeuchi, K.; Oo, A.Z.; Jiang, C.-J. Identification of Pseudomonas strains for the biological control of soybean red crown root rot. Sci. Rep. 2022, 12, 14510. [Google Scholar] [CrossRef] [PubMed]
  17. Berner, D.; Berggren, G.; Pace, M.; White, E.; Gershey, J.; Freedman, J. Red crown rot: Now a major disease of soybeans. La Agric. 1986, 29, 4–5. [Google Scholar]
  18. Lombard, L.; Crous, P.W.; Wingfield, B.D.; Wingfield, M.J. Phylogeny and systematics of the genus Calonectria. Stud. Mycol. 2010, 66, 31–69. [Google Scholar] [CrossRef] [PubMed]
  19. Wang, H.; Wu, J.; Fang, L.; Xie, Y.; Wang, L. First Report of Calonectria ilicicola Causing Fruit Rot on Postharvest Prunus persica in Zhejiang Province, China. Plant Dis. 2023, 107, 3313. [Google Scholar] [CrossRef] [PubMed]
  20. Liu, H.; Shen, Y.; Chang, H.; Tseng, M.; Lin, Y. First report of soybean red crown rot caused by Calonectria ilicicola in Taiwan. Plant Dis. 2020, 104, 979. [Google Scholar] [CrossRef]
  21. Zhang, Y.; Chen, C.; Chen, C.; Chen, J.; Xiang, M.; Wanasinghe, D.N.; Hsiang, T.; Hyde, K.D.; Manawasinghe, I.S. Identification and characterization of Calonectria species associated with plant diseases in Southern China. J. Fungi 2022, 8, 719. [Google Scholar] [CrossRef]
  22. Tazawa, J. Occurrence and control of red crown rot of soybean. Plant Prot. 2013, 67, 46–49. [Google Scholar]
  23. Ochi, S. Studies on the epidemiology and control of red crown rot of soybean. J. Gen. Plant Pathol. 2017, 83, 427–428. [Google Scholar] [CrossRef]
  24. Guan, M.; Pan, R.; Gao, X.; Xu, D.; Deng, Q.; Deng, M. First report of red crown rot caused by Cylindrocladium parasiticum on soybean in Guangdong, Southern China. Plant Dis. 2010, 94, 485. [Google Scholar] [CrossRef] [PubMed]
  25. Jiang, C.-J.; Sugano, S.; Ochi, S.; Kaga, A.; Ishimoto, M. Evaluation of Glycine max and Glycine soja for resistance to Calonectria ilicicola. Agronomy 2020, 10, 887. [Google Scholar] [CrossRef]
  26. Shafique, H.A.; Sultana, V.; Ehteshamul-Haque, S.; Athar, M. Management of soil-borne diseases of organic vegetables. J. Plant Prot. Res. 2016, 56, 221–230. [Google Scholar] [CrossRef]
  27. Jiang, H.; Lv, S.; Zhou, C.; Qu, S.; Liu, F.; Sun, H.; Zhao, X.; Han, Y. Identification of QTL, QTL-by-environment interactions, and their candidate genes for resistance HG Type 0 and HG Type 1.2. 3.5. 7 in soybean using 3VmrMLM. Front. Plant Sci. 2023, 14, 1177345. [Google Scholar] [CrossRef]
  28. Kantartzi, S.K.; McAllister, K.; Lee, Y.-C. Targeting and mapping resistance to Cercospora sojina in two elite soybean (Glycine max L.) populations. J. Plant Breed. Crop Sci. 2023, 15, 118–128. [Google Scholar]
  29. Ferreira, E.G.C.; Marcelino-Guimarães, F.C. Mapping Major Disease Diseases Resistance Resistance Genes in Soybean by Genome-Wide Association Studies. Methods Mol. Biol. 2022, 2481, 313–340. [Google Scholar] [PubMed]
  30. Xiong, H.; Chen, Y.; Pan, Y.-B.; Wang, J.; Lu, W.; Shi, A. A genome-wide association study and genomic prediction for Phakopsora pachyrhizi resistance in soybean. Front. Plant Sci. 2023, 14, 1179357. [Google Scholar] [CrossRef]
  31. Sang, Y.; Liu, X.; Yuan, C.; Yao, T.; Li, Y.; Wang, D.; Zhao, H.; Wang, Y. Genome-wide association study on resistance of cultivated soybean to Fusarium oxysporum root rot in Northeast China. BMC Plant Biol. 2023, 23, 625. [Google Scholar] [CrossRef]
  32. Zhao, F.; Cheng, W.; Wang, Y.; Gao, X.; Huang, D.; Kong, J.; Antwi-Boasiako, A.; Zheng, L.; Yan, W.; Chang, F. Identification of novel genomic regions for bacterial leaf pustule (BLP) resistance in soybean (Glycine max L.) via integrating linkage mapping and association analysis. Int. J. Mol. Sci. 2022, 23, 2113. [Google Scholar] [CrossRef]
  33. Fournier, T.; Abou Saada, O.; Hou, J.; Peter, J.; Caudal, E.; Schacherer, J. Extensive impact of low-frequency variants on the phenotypic landscape at population-scale. eLife 2019, 8, e49258. [Google Scholar] [CrossRef] [PubMed]
  34. Chanchu, T.; Yimram, T.; Chankaew, S.; Kaga, A.; Somta, P. Mapping QTLs Controlling Soybean Rust Disease Resistance in Chiang Mai 5, an Induced Mutant Cultivar. Genes 2022, 14, 19. [Google Scholar] [CrossRef]
  35. Yamanaka, N.; Aoyagi, L.N.; Hossain, M.M.; Aoyagi, M.B.; Muraki, Y. Genetic Mapping of Seven Kinds of Locus for Resistance to Asian Soybean Rust. Plants 2023, 12, 2263. [Google Scholar] [CrossRef]
  36. Kim, K.H.; Park, J.-H.; Kim, M.Y.; Heu, S.; Lee, S.-H. Genetic mapping of novel symptom in response to soybean bacterial leaf pustule in PI 96188. J. Crop Sci. Biotechnol. 2011, 14, 119–123. [Google Scholar] [CrossRef]
  37. Li, D.; Sun, M.; Han, Y.; Teng, W.; Li, W. Identification of QTL underlying soluble pigment content in soybean stems related to resistance to soybean white mold (Sclerotinia sclerotiorum). Euphytica 2010, 172, 49–57. [Google Scholar] [CrossRef]
  38. Huynh, T.; Bastien, M.; Iquira, E.; Turcotte, P.; Belzile, F. Identification of QTLs associated with partial resistance to white mold in soybean using field-based inoculation. Crop Sci. 2010, 50, 969–979. [Google Scholar] [CrossRef]
  39. Han, Y.; Teng, W.; Yu, K.; Poysa, V.; Anderson, T.; Qiu, L.; Lightfoot, D.A.; Li, W. Mapping QTL tolerance to Phytophthora root rot in soybean using microsatellite and RAPD/SCAR derived markers. Euphytica 2008, 162, 231–239. [Google Scholar] [CrossRef]
  40. Zhang, Z.; Hao, J.; Yuan, J.; Song, Q.; Hyten, D.L.; Cregan, P.B.; Zhang, G.; Gu, C.; Li, M.; Wang, D. Phytophthora root rot resistance in soybean E00003. Crop Sci. 2014, 54, 492–499. [Google Scholar] [CrossRef]
  41. Niu, J.; Guo, N.; Sun, J.; Li, L.; Cao, Y.; Li, S.; Huang, J.; Zhao, J.; Zhao, T.; Xing, H. Fine mapping of a resistance gene RpsHN that controls Phytophthora sojae using recombinant inbred lines and secondary populations. Front. Plant Sci. 2017, 8, 538. [Google Scholar] [CrossRef]
  42. McAllister, K.R.; Lee, Y.-C.; Kantartzi, S.K. QTL mapping for resistance to Cercospora sojina in Essex Forrest soybean (Glycine max L.) lines. J. Plant Breed. Crop Sci. 2021, 13, 14–22. [Google Scholar]
  43. Pham, A.-T.; Harris, D.K.; Buck, J.; Hoskins, A.; Serrano, J.; Abdel-Haleem, H.; Cregan, P.; Song, Q.; Boerma, H.R.; Li, Z. Fine mapping and characterization of candidate genes that control resistance to Cercospora sojina K. Hara in two soybean germplasm accessions. PLoS ONE 2015, 10, e0126753. [Google Scholar] [CrossRef] [PubMed]
  44. Wei, T.; Simko, V. R package “Corrplot”: Visualization of a Correlation Matrix (Version 0.84). 2017. Available online: https://github.com/taiyun/corrplot (accessed on 7 January 2023).
  45. Zhang, X.; Hina, A.; Song, S.; Kong, J.; Bhat, J.A.; Zhao, T. Whole-genome mapping identified novel “QTL hotspots regions” for seed storability in soybean (Glycine max L.). BMC Genom. 2019, 20, 499. [Google Scholar] [CrossRef] [PubMed]
  46. Huang, X.; Feng, Q.; Qian, Q.; Zhao, Q.; Wang, L.; Wang, A.; Guan, J.; Fan, D.; Weng, Q.; Huang, T. High-throughput genotyping by whole-genome resequencing. Genome Res. 2009, 19, 1068–1076. [Google Scholar] [CrossRef] [PubMed]
  47. Cao, Y.; Li, S.; Chen, G.; Wang, Y.; Bhat, J.A.; Karikari, B.; Kong, J.; Gai, J.; Zhao, T. Deciphering the genetic architecture of plant height in soybean using two RIL populations sharing a common M8206 parent. Plants 2019, 8, 373. [Google Scholar] [CrossRef] [PubMed]
  48. Wang, S. Windows QTL Cartographer 2.5; Department of Statistics, North Carolina State University: Raleigh, NC, USA, 2007; Available online: https://brcwebportal.cos.ncsu.edu/qtlcart/WQTLCart.htm (accessed on 1 January 2019).
  49. Song, Q.; Marek, L.; Shoemaker, R.; Lark, K.; Concibido, V.; Delannay, X.; Specht, J.E.; Cregan, P. A new integrated genetic linkage map of the soybean. Theor. Appl. Genet. 2004, 109, 122–128. [Google Scholar] [CrossRef] [PubMed]
  50. Churchill, G.A.; Doerge, R.W. Empirical threshold values for quantitative trait mapping. Genetics 1994, 138, 963–971. [Google Scholar] [CrossRef] [PubMed]
  51. Voorrips, R. MapChart: Software for the graphical presentation of linkage maps and QTLs. J. Hered. 2002, 93, 77–78. [Google Scholar] [CrossRef] [PubMed]
  52. McCouch, S.R.; Chen, X.; Panaud, O.; Temnykh, S.; Xu, Y.; Cho, Y.G.; Huang, N.; Ishii, T.; Blair, M. Microsatellite marker development, mapping and applications in rice genetics and breeding. In Oryza: From Molecule to Plant; Sasaki, T., Moore, G., Eds.; Springer: Dordrecht, The Netherlands, 1997. [Google Scholar] [CrossRef]
  53. Tian, T.; Liu, Y.; Yan, H.; You, Q.; Yi, X.; Du, Z.; Xu, W.; Su, Z. agriGO v2. 0: A GO analysis toolkit for the agricultural community, 2017 update. Nucleic Acids Res. 2017, 45, W122–W129. [Google Scholar] [CrossRef]
  54. Kobayashi, M.; Win, K.T.; Jiang, C.-J. Soybean hypocotyls prevent Calonectria ilicicola invasion by multi-layered defenses. Front. Plant Sci. 2022, 12, 813578. [Google Scholar] [CrossRef]
  55. Chen, C.; Chen, H.; He, Y.; Xia, R. TBtools, a toolkit for biologists integrating various biological data handling tools with a user-friendly interface. BioRxiv 2018, 1, 289660. [Google Scholar] [CrossRef]
  56. Hung, H.Y.; Browne, C.; Guill, K.; Coles, N.; Eller, M.; Garcia, A.; Lepak, N.; Melia-Hancock, S.; Oropeza-Rosas, M.; Salvo, S.; et al. The relationship between parental genetic or phenotypic divergence and progeny variation in the maize nested association mapping population. Heredity 2012, 108, 490–499. [Google Scholar] [CrossRef] [PubMed]
  57. Gelli, M.; Konda, A.R.; Liu, K.; Zhang, C.; Clemente, T.E.; Holding, D.R.; Dweikat, I.M. Validation of QTL mapping and transcriptome profiling for identification of candidate genes associated with nitrogen stress tolerance in sorghum. BMC Plant Biol. 2017, 17, 123. [Google Scholar] [CrossRef] [PubMed]
  58. Win, K.T.; Jiang, C.-J. A fresh weight-based method for evaluating soybean resistance to red crown rot. Breed. Sci. 2021, 71, 384–389. [Google Scholar] [CrossRef] [PubMed]
  59. Nakajima, T.; Sakai, S.; Gomi, T.; Kikuchi, A. Development of methods for assessing resistance to black root rot caused by Calonectria crotalariae in soybean [Glycine max] and screening for resistant germplasm. Bull. Tohoku Natl. Agric. Exp. Stn. 1994, 88, 39–56. [Google Scholar]
  60. Antwi-Boasiako, A.; Jia, S.; Liu, J.; Guo, N.; Chen, C.; Karikari, B.; Feng, J.; Zhao, T. Identification and Genetic Dissection of Resistance to Red Crown Rot Disease in a Diverse Soybean Germplasm Population. Plants 2024, 13, 940. [Google Scholar] [CrossRef] [PubMed]
  61. Cao, Y.; Li, S.; Wang, Z.; Chang, F.; Kong, J.; Gai, J.; Zhao, T. Identification of Major Quantitative Trait Loci for Seed Oil Content in Soybeans by Combining Linkage and Genome-Wide Association Mapping. Front. Plant Sci. 2017, 8, 1222. [Google Scholar] [CrossRef] [PubMed]
  62. Li, Z.; Xu, Y. Bulk segregation analysis in the NGS era: A review of its teenage years. Plant J. 2022, 109, 1355–1374. [Google Scholar] [CrossRef] [PubMed]
  63. Mackay, I.; Caligari, P. Efficiencies of F2 and backcross generations for bulked segregant analysis using dominant markers. Crop Sci. 2000, 40, 626–630. [Google Scholar] [CrossRef]
  64. Jansen, R.; Van Ooijen, J.; Stam, P.; Lister, C.; Dean, C. Genotype-by-environment interaction in genetic mapping of multiple quantitative trait loci. Theor. Appl. Genet. 1995, 91, 33–37. [Google Scholar] [CrossRef]
  65. Han, Y.; Zhang, Y.; Wu, D.; Zhao, X.; Teng, W.; Li, D.; Li, W. Identification of novel quantitative trait loci associated with tolerance to Phytophthora root rot in the soybean cultivar Hefeng 25 using two recombinant inbred line populations. Can. J. Plant Sci. 2017, 97, 827–834. [Google Scholar] [CrossRef]
  66. Antwi-Boasiako, A.; Zheng, L.; Begum, N.; Amoah, S.; Zhao, T. Progress towards germplasm evaluation and genetic improvement for resistance to Sclerotinia white mold in soybean. Euphytica 2021, 217, 178. [Google Scholar] [CrossRef]
  67. Huang, X.; Han, B. Natural variations and genome-wide association studies in crop plants. Annu. Rev. Plant Biol. 2014, 65, 531–551. [Google Scholar] [CrossRef] [PubMed]
  68. Ding, L.-N.; Li, Y.-T.; Wu, Y.-Z.; Li, T.; Geng, R.; Cao, J.; Zhang, W.; Tan, X.-L. Plant disease resistance-related signaling pathways: Recent progress and future prospects. Int. J. Mol. Sci. 2022, 23, 16200. [Google Scholar] [CrossRef] [PubMed]
  69. Ngou, B.P.M.; Ding, P.; Jones, J.D. Thirty years of resistance: Zig-zag through the plant immune system. Plant Cell 2022, 34, 1447–1478. [Google Scholar] [CrossRef] [PubMed]
  70. Sertsuvalkul, N.; DeMell, A.; Dinesh-Kumar, S.P. The complex roles of autophagy in plant immunity. FEBS Lett. 2022, 596, 2163–2171. [Google Scholar] [CrossRef] [PubMed]
  71. Devanna, B.N.; Jaswal, R.; Singh, P.K.; Kapoor, R.; Jain, P.; Kumar, G.; Sharma, Y.; Samantaray, S.; Sharma, T.R. Role of transporters in plant disease resistance. Physiol. Plant. 2021, 171, 849–867. [Google Scholar] [CrossRef] [PubMed]
  72. Strobel, G.A. A mechanism of disease resistance in plants. Sci. Am. 1975, 232, 80–89. [Google Scholar] [CrossRef]
  73. Xi, Y.; Cesari, S.; Kroj, T. Insight into the structure and molecular mode of action of plant paired NLR immune receptors. Essays Biochem. 2022, 66, 513–526. [Google Scholar]
  74. Zhou, L.; Deng, S.; Xuan, H.; Fan, X.; Sun, R.; Zhao, J.; Wang, H.; Guo, N.; Xing, H. A novel TIR-NBS-LRR gene regulates immune response to Phytophthora root rot in soybean. Crop J. 2022, 10, 1644–1653. [Google Scholar] [CrossRef]
  75. Ameline-Torregrosa, C.; Wang, B.-B.; O’Bleness, M.S.; Deshpande, S.; Zhu, H.; Roe, B.; Young, N.D.; Cannon, S.B. Identification and characterization of nucleotide-binding site-leucine-rich repeat genes in the model plant Medicago truncatula. Plant Physiol. 2008, 146, 5–21. [Google Scholar] [CrossRef] [PubMed]
  76. Wang, F.; Lu, T.; Zhu, L.; Cao, A.; Xie, S.; Chen, X.; Shen, H.; Xie, Q.; Li, R.; Zhu, J. Multicopper oxidases GbAO and GbSKS are involved in the Verticillium dahliae resistance in Gossypium barbadense. J. Plant Physiol. 2023, 280, 153887. [Google Scholar] [CrossRef] [PubMed]
  77. Liu, J.Z.; Whitham, S.A. Overexpression of a soybean nuclear localized type–III DnaJ domain-containing HSP40 reveals its roles in cell death and disease resistance. Plant J. 2013, 74, 110–121. [Google Scholar] [CrossRef] [PubMed]
  78. Kim, N.H.; Hwang, B.K. Pepper heat shock protein 70a interacts with the type III effector AvrBsT and triggers plant cell death and immunity. Plant Physiol. 2015, 167, 307–322. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Frequency distribution of emergence rate, survival rate, and disease severity in ZM6 and MN RIL populations. The phenotype was averaged from the first and second time screenings in each RIL population. The trend shows the moving average. The arrows show the mean value of corresponding. The horizontal and vertical axis represent trait value and number of genotypes, respectively.
Figure 1. Frequency distribution of emergence rate, survival rate, and disease severity in ZM6 and MN RIL populations. The phenotype was averaged from the first and second time screenings in each RIL population. The trend shows the moving average. The arrows show the mean value of corresponding. The horizontal and vertical axis represent trait value and number of genotypes, respectively.
Agronomy 14 01693 g001
Figure 2. Pearson correlation matrix among the emergence rate (ER), survival rate (SR), and disease severity (DS) used to screen for red crown rot resistance in ZM6 and MN RIL populations (ZM6, upper part, and MN, lower part of the matrix). The phenotype was averaged from the first and second time screenings of the study in each RIL population. The asterisk shows significance at p < 0.05. The color scale of each correlation coefficient is shown in the color legend on the right-hand side of the correlation matrix.
Figure 2. Pearson correlation matrix among the emergence rate (ER), survival rate (SR), and disease severity (DS) used to screen for red crown rot resistance in ZM6 and MN RIL populations (ZM6, upper part, and MN, lower part of the matrix). The phenotype was averaged from the first and second time screenings of the study in each RIL population. The asterisk shows significance at p < 0.05. The color scale of each correlation coefficient is shown in the color legend on the right-hand side of the correlation matrix.
Agronomy 14 01693 g002
Figure 3. Genome-wide view of QTL detected in the ZM6 RIL population for soybean root crown rot disease via the CIM strategy. The phenotypic data were obtained from two time screenings (first time and second time). The disease resistance was screened by three traits: emergence rate (ER), survival rate (SR), and disease severity (DS). The y-axis of the top figure represents a logarithm of odd value with a threshold of 2.5 (sea blue line), while the x-axis represents the chromosome (Chr) number in the soybean genome. The bottom figure denotes the additivity of alleles and their origin from either of the parents: positive values represent allele from Meng8206 (♀), while negative values emanate from Zhengyang (♂). The color legend is shown on the right-hand side of the top figure.
Figure 3. Genome-wide view of QTL detected in the ZM6 RIL population for soybean root crown rot disease via the CIM strategy. The phenotypic data were obtained from two time screenings (first time and second time). The disease resistance was screened by three traits: emergence rate (ER), survival rate (SR), and disease severity (DS). The y-axis of the top figure represents a logarithm of odd value with a threshold of 2.5 (sea blue line), while the x-axis represents the chromosome (Chr) number in the soybean genome. The bottom figure denotes the additivity of alleles and their origin from either of the parents: positive values represent allele from Meng8206 (♀), while negative values emanate from Zhengyang (♂). The color legend is shown on the right-hand side of the top figure.
Agronomy 14 01693 g003
Figure 4. Location of QTLs on the genetic linkage map of ZM6 RIL population. Distances among markers are indicated using the physical location to the left of the populations and the QTL names located at their position on the right side. Chr-Chromosome. Colored bars represent different QTLs.
Figure 4. Location of QTLs on the genetic linkage map of ZM6 RIL population. Distances among markers are indicated using the physical location to the left of the populations and the QTL names located at their position on the right side. Chr-Chromosome. Colored bars represent different QTLs.
Agronomy 14 01693 g004
Figure 5. Genome-wide view of QTL detected in MN RIL population for root crown rot resistance via the CIM strategy. The phenotypic data were obtained from two time screenings (first time and second time). The disease resistance was screened by three traits: emergence rate (ER), survival rate (SR), and disease severity (DS). The y-axis of the top figure represents a logarithm of odd value with a threshold of 2.5 (sea blue line), while the x-axis represents chromosome (Chr) number in the soybean genome. The bottom figure denotes the additivity of alleles and their origin from either of the parents: positive values represent alleles from M8180 (♀), while negative values emanate from NN1138-2 (♂). The color legend is shown on the right-hand side of the top figure.
Figure 5. Genome-wide view of QTL detected in MN RIL population for root crown rot resistance via the CIM strategy. The phenotypic data were obtained from two time screenings (first time and second time). The disease resistance was screened by three traits: emergence rate (ER), survival rate (SR), and disease severity (DS). The y-axis of the top figure represents a logarithm of odd value with a threshold of 2.5 (sea blue line), while the x-axis represents chromosome (Chr) number in the soybean genome. The bottom figure denotes the additivity of alleles and their origin from either of the parents: positive values represent alleles from M8180 (♀), while negative values emanate from NN1138-2 (♂). The color legend is shown on the right-hand side of the top figure.
Agronomy 14 01693 g005
Figure 6. Location of the QTLs on the genetic linkage map of MN RIL populations. Distances among markers are indicated using the physical location to the left of the populations and the QTL names located at their position on the right side. Colored bars represent different QTLs.
Figure 6. Location of the QTLs on the genetic linkage map of MN RIL populations. Distances among markers are indicated using the physical location to the left of the populations and the QTL names located at their position on the right side. Colored bars represent different QTLs.
Agronomy 14 01693 g006
Table 1. Classification for the disease severity of red crown rot.
Table 1. Classification for the disease severity of red crown rot.
ScaleDamage DegreeResistance Degree
0No visible sign of necrotic lesions on the rootImmune
1Only the tap root exhibits small necrotic lesions without obvious changes in its form and shapeHighly Resistant
2Necrotic lesions spread to the crown and root system and seedling mortality less than 10%Resistant
3Roots show serious necrotic lesions with less than 50% loss by rot and seedling mortality of 10–20%Moderately Susceptible
4Roots show serious necrotic lesions with more than 50% root loss by rot and seedling mortality of 21–50%Susceptible
5Over 50% of root loss by rot with seedling mortality of more than 50%Highly Susceptible
Table 2. Response of two soybean RIL populations and their parents to RCR.
Table 2. Response of two soybean RIL populations and their parents to RCR.
Population/
Environment a
Trait bMean of Parent (%) cRIL Populations d
Mean ± SEMin.Max.KurtosisSkewness
ZM6/2ER73.33 ± 3.3383.33 ± 3.3383.11 ± 0.0063.33100.00−0.860.13
SR82.14 ± 3.57100 ± 0.0086.07 ± 0.0054.76100.00−0.09−0.63
DS3.33 ± 0.331.00 ± 0.002.83 ± 1.001.004.00−0.50−0.10
ZM6/1ER76.67 ± 3.3393.33 ± 3.3380.44 ± 1.0056.67100.00−0.30−0.32
SR82.74 ± 3.90100 ± 0.0081.79 ± 1.0039.68100.00−0.05−0.59
DS4.33 ± 0.332.00 ± 0.003.46 ± 1.001.004.33−0.50−0.10
MN/2ER100.00 ± 0.0080.00 ± 0.0082.50 ± 5.2035.67100.001.31−1.08
SR100.00 ± 0.0083.33 ± 6.9396.45 ± 0.0474.17100.003.86−1.93
DS1.67 ± 34.643.67 ± 5.752.46 ± 6.091.004.00−0.660.30
MN/1ER96.67 ± 5.9770.00 ± 0.0078.44 ± 3.9260.0090.003.67−1.36
SR100.00 ± 0.0076.67 ± 8.0087.86 ± 3.0052.86100.001.28−1.03
DS2.00 ± 0.003.6 ± 6.003.00 ± 6.001.334.00−0.53−0.52
a Environment: the first and second time screenings (1 and 2) in each RIL population. b ER-emergence rate, SR—survival rate, and DS—disease severity. c Zhengyang (♀)and Meng8206 (♂) for ZM6 RILs, while M8180 (♀)and NN1138-2 (♂) for MN RILs. d SE—standard error of mean, min—minimum, and max—maximum.
Table 3. QTLs detected for RCR disease resistance via ER, SR, and DS in the ZM6 RIL population via the CIM strategy.
Table 3. QTLs detected for RCR disease resistance via ER, SR, and DS in the ZM6 RIL population via the CIM strategy.
QTL aPos (cM) bLOD cAdd dPVE (%) eCI (cM) fPhysical Region (bp)Env g
qER-7-1zm664.513.825.8110.3761.4–69.715,375,767–17,747,9731
qER-8-1zm6150.813.225.048.40147.7–154.741,485,100–42,915,2552
qER-10-1zm616.613.25−5.248.5612.8–19.61,694,367–2,679,8371
qER-10-2zm6106.414.04−6.0010.70105.1–109.943,900,754–44,741,9602
qER-11-1zm673.813.867.2010.4071.2–75.514,962,695–15,949,2961
qSR-6-1zm69.712.975.327.786.7–15.61,813,130–3,196,5552
qSR-7-1zm648.814.557.0812.7146.2–50.19,304,376–10,428,5322
qSR-7-2zm656.015.197.7114.8753–6214,134,797–15,903,2802
59.213.216.008.6252.7–61.414,134,797–15,452,7981
qSR-10-1zm616.613.37−6.008.8810.3–19.61,603,735–2,732,8801
qSR-11-1zm676.414.197.6611.3475.5–81.815,676,274–16,816,8001
qSR-17-1zm646.413.32−5.938.7432.8–53.76,777,393–9,645,3252
qDS-11-1zm6139.113.81−0.269.55136–140.837,603,249–38,850,6962
qDS-13-1zm68.913.690.279.263.1–11.24,552,834–5,592,4481
qDS-18-1zm6120.113.63−0.289.57118–122.759,218,992–60,685,6751
qDS-18-2zm6128.515.04−0.3212.95124.8–129.561,300,197–62,014,7061
128.513.57−0.258.89124.2–129.560,909,812–62,014,7062
a QTLs detected in different environments at the same, adjacent, or overlapping marker intervals were considered the same QTL. b Pos—Position of QTL in centiMorgan. c LOD—Logarithm of odds. d Add—Addictive, indicating the origin of beneficial alleles; positive values represent allele from Meng8206 (♀), while negative values emanated from Zhengyang (♂)). e PEV—phenotypic variance explained (%) expressed by the QTL. f CI—Confidence interval. g Env—Environment, the phenotypic data obtained from two time screenings (first time and first time).
Table 4. QTLs detected for RCR disease resistance via ER, SR, and DS in the MN RIL population via the CIM strategy.
Table 4. QTLs detected for RCR disease resistance via ER, SR, and DS in the MN RIL population via the CIM strategy.
QTL aPos(cM) bLOD cAdd dPVE (%) eCI (cM) fPhysical Region (bp)Env g
qER-1-1mn27.117.61−7.2425.2425–298,823,531–22,021,3581
qER-10-1mn0.013.38−3.8110.720–82,275,280–4,174,7911
qER-8-1mn0.012.633.889.590–1014,650,727–11,805,2462
qER-8-2mn48.713.35−4.5112.9233.2–57.98,218,976–18,160,0782
qER-15-1mn16.113.444.3112.6613.5–22.35,519,255–113,832,932
qSR-2-1mn111.512.672.979.32105.5–125.67,361,306–15,293,2251
qSR-17-1mn34.012.91−3.049.8331.4–34.618,508,753–33,427,2031
qSR-1-1mn19.812.64−1.8310.714.2–23.33,617,559–8,823,7882
qSR-1-2mn29.012.73−1.688.6927.1–3521,174,218–44,479,8952
qDS-1-1mn37.912.740.2510.632.8–41.310,404,837–43,932,9071
qDS-1-2mn46.613.220.2711.744.8–52.638,197,263–49,324,4051
qDS-1-3mn56.713.67−0.2611.8451.4–58.241,313,930–50,832,2922
qDS-4-1mn154.413.28−0.2510.51152.4–16411,631,171–36,956,7692
qDS-13-1mn3.413.65−0.2712.280–7.8688,713–5,592,4482
a QTLs detected in different environments at the same, adjacent, or overlapping marker intervals were considered the same QTL. b Pos—Position of QTL in centiMorgan. c LOD—Logarithm of odds. d Add—Addictive, indicating the origin of beneficial alleles; positive values represent alleles from Meng8206 (♀), while negative values emanate from Zhengyang (♂). e phenotypic variance explained (%) expressed by the QTL. f CI—Confidence interval of the QTL. g Env—Environment.
Table 5. The six QTL hotspots identified in the ZM6 and MN RIL populations.
Table 5. The six QTL hotspots identified in the ZM6 and MN RIL populations.
QTL Hotspot NameQTL NameLOD aAdd bPVE (%) cPhysical Region (bp)
Hotspot AqER-1-1mn7.61−7.2425.248,823,531–44,479,895
qSR-1-2mn2.73−1.688.69
Hotspot BqSR-7-2zm65.197.7114.8714,134,797–15,903,280
3.206.008.62
Hotspot CqER-10-1zm63.25−5.248.561,603,735–2,732,880
qSR-10-1zm63.37−6.008.88
Hotspot DqER-11-1zm63.867.2010.4014,962,695–16,816,800
qSR-11-1zm64.197.6611.34
Hotspot EqDS-13-1zm63.81−0.269.55688,713–5,592,448
qDS-13-1mn3.65−0.2712.28
Hotspot FqDs-18-2zm65.04−0.3212.9560,909,812–62,014,706
3.57−0.258.89
a LOD—Logarithm of odds. b Add—Adductive indicating the origin of beneficial alleles; positive values represent alleles from Meng8206 (♀), while negative values emanate from Zhengyang (♂). c PVE—phenotypic variance explained (%) explained by the QTL.
Table 6. Predictive gene annotation information within the six QTL hotspots identified in the ZM6 and MN RIL populations.
Table 6. Predictive gene annotation information within the six QTL hotspots identified in the ZM6 and MN RIL populations.
QTL HotspotGene Mapped IDs aAnnotation Descriptions b
AGlyma.01G127100Disease resistance-responsive; Dirigent-like protein
Glyma.01G126600Disease resistance-responsive; Dirigent-like protein
Glyma.01G112300Signal transduction; Leucine Rich Repeat
Glyma.01G127200Disease resistance-responsive; Dirigent-like protein
Glyma.01G127700Signal transduction; Defense response
BGlyma.07G133900Lignin catabolic process; multicopper oxidase
Glyma.07G134100Lignin catabolic process; multicopper oxidase
CGlyma.10G019900Glutathione metabolic process; glutathione S-transferase
Glyma.10G023300Protein phosphorylation; serine/threonine protein kinase
Glyma.10G023400Protein phosphorylation; serine/threonine protein kinase
Glyma.10G028200Bifunctional inhibitor/lipid-transfer protein/seed storage 2S albumin superfamily protein
Glyma.10G025700F-box family protein; protein binding
DGlyma.11G181200F-box family protein; protein binding
Glyma.11G192300Oxidation-reduction process; pheophorbide an oxygenase
Glyma.11G193600Cellular glucan metabolic process; cell wall biogenesis
EGlyma.13G066100DNA repair; ATP binding
Glyma.13G069200Zinc finger (AN1-like) family protein
Glyma.13G076200Defense response; Signal transduction; Leucine Rich Repeat
FGlyma.18G287000Defense response; signal transduction
Glyma.18G289100ATP binding, hsp70 protein
Glyma.18G289600ATP binding, hsp70 protein
Glyma.18G293200Drug transmembrane transport; multidrug resistance protein
Glyma.18G293300Drug transmembrane transport; multidrug resistance protein
a model gene retrieved on reference genome V2 from SoyBase (https://www.soybase.org/ accessed on 4 April 2024). b Gene ontology of the retrieved genes from SoyBase (https://www.soybase.org/ accessed on 3 April 2024).
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

Antwi-Boasiako, A.; Zhang, C.; Almakas, A.; Liu, J.; Jia, S.; Guo, N.; Chen, C.; Zhao, T.; Feng, J. Identification of QTLs and Candidate Genes for Red Crown Rot Resistance in Two Recombinant Inbred Line Populations of Soybean [Glycine max (L.) Merr.]. Agronomy 2024, 14, 1693. https://doi.org/10.3390/agronomy14081693

AMA Style

Antwi-Boasiako A, Zhang C, Almakas A, Liu J, Jia S, Guo N, Chen C, Zhao T, Feng J. Identification of QTLs and Candidate Genes for Red Crown Rot Resistance in Two Recombinant Inbred Line Populations of Soybean [Glycine max (L.) Merr.]. Agronomy. 2024; 14(8):1693. https://doi.org/10.3390/agronomy14081693

Chicago/Turabian Style

Antwi-Boasiako, Augustine, Chunting Zhang, Aisha Almakas, Jiale Liu, Shihao Jia, Na Guo, Changjun Chen, Tuanjie Zhao, and Jianying Feng. 2024. "Identification of QTLs and Candidate Genes for Red Crown Rot Resistance in Two Recombinant Inbred Line Populations of Soybean [Glycine max (L.) Merr.]" Agronomy 14, no. 8: 1693. https://doi.org/10.3390/agronomy14081693

APA Style

Antwi-Boasiako, A., Zhang, C., Almakas, A., Liu, J., Jia, S., Guo, N., Chen, C., Zhao, T., & Feng, J. (2024). Identification of QTLs and Candidate Genes for Red Crown Rot Resistance in Two Recombinant Inbred Line Populations of Soybean [Glycine max (L.) Merr.]. Agronomy, 14(8), 1693. https://doi.org/10.3390/agronomy14081693

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