Dissecting the Genetic Basis of Resistance to Soybean Cyst Nematode Combining Linkage and Association Mapping

TWAS facilitates gene-scale trait genetic dissection through gene expression, structural variations, and alternative splicing in soybean

A genome-wide association study (GWAS) identifies trait-associated loci, but identifying the causal genes can be a bottleneck, due in part to slow decay of linkage disequilibrium (LD). A transcriptome-wide association study (TWAS) addresses this issue by identifying gene expression-phenotype associations or integrating gene expression quantitative trait loci with GWAS results. Here, we used self-pollinated soybean (Glycine max [L.] Merr.) as a model to evaluate the application of TWAS to the genetic dissection of traits in plant species with slow LD decay. We generated RNA sequencing data for a soybean diversity panel and identified the genetic expression regulation of 29 286 soybean genes. Different TWAS solutions were less affected by LD and were robust to the source of expression, identifing known genes related to traits from different tissues and developmental stages. The novel pod-color gene L2 was identified via TWAS and functionally validated by genome editing. By introducing a new exon proportion feature, we significantly improved the detection of expression variations that resulted from structural variations and alternative splicing. As a result, the genes identified through our TWAS approach exhibited a diverse range of causal variations, including SNPs, insertions or deletions, gene fusion, copy number variations, and alternative splicing. Using this approach, we identified genes associated with flowering time, including both previously known genes and novel genes that had not previously been linked to this trait, providing insights complementary to those from GWAS. In summary, this study supports the application of TWAS for candidate gene identification in species with low rates of LD decay.

Application of SVR-Mediated GWAS for Identification of Durable Genetic Regions Associated with Soybean Seed Quality Traits

Soybean (Glycine max L.) is an important food-grade strategic crop worldwide because of its high seed protein and oil contents. Due to the negative correlation between seed protein and oil percentage, there is a dire need to detect reliable quantitative trait loci (QTL) underlying these traits in order to be used in marker-assisted selection (MAS) programs. Genome-wide association study (GWAS) is one of the most common genetic approaches that is regularly used for detecting QTL associated with quantitative traits. However, the current approaches are mainly focused on estimating the main effects of QTL, and, therefore, a substantial statistical improvement in GWAS is required to detect associated QTL considering their interactions with other QTL as well. This study aimed to compare the support vector regression (SVR) algorithm as a common machine learning method to fixed and random model circulating probability unification (FarmCPU), a common conventional GWAS method in detecting relevant QTL associated with soybean seed quality traits such as protein, oil, and 100-seed weight using 227 soybean genotypes. The results showed a significant negative correlation between soybean seed protein and oil concentrations, with heritability values of 0.69 and 0.67, respectively. In addition, SVR-mediated GWAS was able to identify more relevant QTL underlying the target traits than the FarmCPU method. Our findings demonstrate the potential use of machine learning algorithms in GWAS to detect durable QTL associated with soybean seed quality traits suitable for genomic-based breeding approaches. This study provides new insights into improving the accuracy and efficiency of GWAS and highlights the significance of using advanced computational methods in crop breeding research.

Large-scale data mining pipeline for identifying novel soybean genes involved in resistance against the soybean cyst nematode

The soybean cyst nematode (SCN) [Heterodera glycines Ichinohe] is a devastating pathogen of soybean [Glycine max (L.) Merr.] that is rapidly becoming a global economic issue. Two loci conferring SCN resistance have been identified in soybean, Rhg1 and Rhg4; however, they offer declining protection. Therefore, it is imperative that we identify additional mechanisms for SCN resistance. In this paper, we develop a bioinformatics pipeline to identify protein-protein interactions related to SCN resistance by data mining massive-scale datasets. The pipeline combines two leading sequence-based protein-protein interaction predictors, the Protein-protein Interaction Prediction Engine (PIPE), PIPE4, and Scoring PRotein INTeractions (SPRINT) to predict high-confidence interactomes. First, we predicted the top soy interacting protein partners of the Rhg1 and Rhg4 proteins. Both PIPE4 and SPRINT overlap in their predictions with 58 soybean interacting partners, 19 of which had GO terms related to defense. Beginning with the top predicted interactors of Rhg1 and Rhg4, we implement a "guilt by association" in silico proteome-wide approach to identify novel soybean genes that may be involved in SCN resistance. This pipeline identified 1,082 candidate genes whose local interactomes overlap significantly with the Rhg1 and Rhg4 interactomes. Using GO enrichment tools, we highlighted many important genes including five genes with GO terms related to response to the nematode (GO:0009624), namely, Glyma.18G029000, Glyma.11G228300, Glyma.08G120500, Glyma.17G152300, and Glyma.08G265700. This study is the first of its kind to predict interacting partners of known resistance proteins Rhg1 and Rhg4, forming an analysis pipeline that enables researchers to focus their search on high-confidence targets to identify novel SCN resistance genes in soybean.

Deciphering the genetic basis of resistance to soybean cyst nematode combining IBD and association mapping

IBD analysis clarified the dynamics of chromosomal recombination during the ZP pedigree breeding process and identified ten genomic regions resistant to SCN race3 combining association mapping. Soybean cyst nematode (SCN, Heterodera glycines Ichinohe) is one of the most devastating pathogens for soybean production worldwide. The cultivar Zhongpin03-5373 (ZP), derived from SCN-resistant progenitor parents, Peking, PI 437654 and Huipizhi Heidou, is an elite line with high resistance to SCN race3. In the current study, a pedigree variation map was generated for ZP and its ten progenitors using 3,025,264 high-quality SNPs identified from an average of 16.2 × re-sequencing for each genome. Through identity by decent (IBD) tracking, we showed the dynamic change of genome and detected important IBD fragments, which revealed the comprehensively artificial selection of important traits during ZP breeding process. A total of 2,353 IBD fragments related to SCN resistance including SCN-resistant genes rhg1, rhg4 and NSF_RAN07 were identified based on the resistant-related genetic paths. Moreover, 23 genomic regions underlying resistance to SCN race3 were identified by genome-wide association study (GWAS) in 481 re-sequenced cultivated soybeans. Ten common loci were found by both IBD tracking and GWAS analysis. Haplotype analysis of 16 potential candidate genes suggested a causative SNP (C/T, - 1065) located in the promoter of Glyma.08G096500 and encoding a predicted TIFY5b-related protein on chr8 was highly correlated with SCN race3 resistance. Our results more thoroughly elucidated the dynamics of genomic fragments during ZP pedigree breeding and the genetic basis of SCN resistance, which will provide useful information for gene cloning and the development of resistant soybean cultivars using a marker-assisted selection approach.

Breeding for disease resistance in soybean: a global perspective

KEY MESSAGE

This review provides a comprehensive atlas of QTLs, genes, and alleles conferring resistance to 28 important diseases in all major soybean production regions in the world. Breeding disease-resistant soybean [Glycine max (L.) Merr.] varieties is a common goal for soybean breeding programs to ensure the sustainability and growth of soybean production worldwide. However, due to global climate change, soybean breeders are facing strong challenges to defeat diseases. Marker-assisted selection and genomic selection have been demonstrated to be successful methods in quickly integrating vertical resistance or horizontal resistance into improved soybean varieties, where vertical resistance refers to R genes and major effect QTLs, and horizontal resistance is a combination of major and minor effect genes or QTLs. This review summarized more than 800 resistant loci/alleles and their tightly linked markers for 28 soybean diseases worldwide, caused by nematodes, oomycetes, fungi, bacteria, and viruses. The major breakthroughs in the discovery of disease resistance gene atlas of soybean were also emphasized which include: (1) identification and characterization of vertical resistance genes reside rhg1 and Rhg4 for soybean cyst nematode, and exploration of the underlying regulation mechanisms through copy number variation and (2) map-based cloning and characterization of Rps11 conferring resistance to 80% isolates of Phytophthora sojae across the USA. In this review, we also highlight the validated QTLs in overlapping genomic regions from at least two studies and applied a consistent naming nomenclature for these QTLs. Our review provides a comprehensive summary of important resistant genes/QTLs and can be used as a toolbox for soybean improvement. Finally, the summarized genetic knowledge sheds light on future directions of accelerated soybean breeding and translational genomics studies.

The GmSNAP11 Contributes to Resistance to Soybean Cyst Nematode Race 4 in Glycine max

Soybean cyst nematode (SCN) has devastating effects on soybean production, making it crucial to identify genes conferring SCN resistance. Here we employed next-generation sequencing-based bulked segregant analysis (BSA) to discover genomic regions, candidate genes, and diagnostic markers for resistance to SCN race 4 (SCN4) in soybean. Phenotypic analysis revealed highly significant differences among the reactions of 145 recombinant inbred lines (RILs) to SCN4. In combination with euclidean distance (ED) and Δsingle-nucleotide polymorphism (SNP)-index analyses, we identified a genomic region on Gm11 (designated as rhg1-paralog) associated with SCN4 resistance. Overexpression and RNA interference analyzes of the two candidate genes identified in this region (GmPLAC8 and GmSNAP11) revealed that only GmSNAP11 significantly contributes to SCN4 resistance. We developed a diagnostic marker for GmSNAP11. Using this marker, together with previously developed markers for SCN-resistant loci, rhg1 and Rhg4, we evaluated the relationship between genotypes and SCN4 resistance in 145 RILs and 30 soybean accessions. The results showed that all the SCN4-resistant lines harbored all the three loci, however, some lines harboring the three loci were still susceptible to SCN4. This suggests that these three loci are necessary for the resistance to SCN4, but they alone cannot confer full resistance. The GmSNAP11 and the diagnostic markers developed could be used in genomic-assisted breeding to develop soybean varieties with increased resistance to SCN4.

Machine-Learning-Based Genome-Wide Association Studies for Uncovering QTL Underlying Soybean Yield and Its Components

A genome-wide association study (GWAS) is currently one of the most recommended approaches for discovering marker-trait associations (MTAs) for complex traits in plant species. Insufficient statistical power is a limiting factor, especially in narrow genetic basis species, that conventional GWAS methods are suffering from. Using sophisticated mathematical methods such as machine learning (ML) algorithms may address this issue and advance the implication of this valuable genetic method in applied plant-breeding programs. In this study, we evaluated the potential use of two ML algorithms, support-vector machine (SVR) and random forest (RF), in a GWAS and compared them with two conventional methods of mixed linear models (MLM) and fixed and random model circulating probability unification (FarmCPU), for identifying MTAs for soybean-yield components. In this study, important soybean-yield component traits, including the number of reproductive nodes (RNP), non-reproductive nodes (NRNP), total nodes (NP), and total pods (PP) per plant along with yield and maturity, were assessed using a panel of 227 soybean genotypes evaluated at two locations over two years (four environments). Using the SVR-mediated GWAS method, we were able to discover MTAs colocalized with previously reported quantitative trait loci (QTL) with potential causal effects on the target traits, supported by the functional annotation of candidate gene analyses. This study demonstrated the potential benefit of using sophisticated mathematical approaches, such as SVR, in a GWAS to complement conventional GWAS methods for identifying MTAs that can improve the efficiency of genomic-based soybean-breeding programs.

Identification of Candidate Genes Controlling Soybean Cyst Nematode Resistance in "Handou 10" Based on Genome and Transcriptome Analyzes

Soybean cyst nematode (SCN; Heterodera glycines Ichinohe) is a highly destructive pathogen for soybean production worldwide. The use of resistant varieties is the most effective way of preventing yield loss. Handou 10 is a commercial soybean variety with desirable agronomic traits and SCN resistance, however genes underlying the SCN resistance in the variety are unknown. An F_2:8 recombinant inbred line (RIL) population derived from a cross between Zheng 9525 (susceptible) and Handou 10 was developed and its resistance to SCN HG type 2.5.7 (race 1) and 1.2.5.7 (race 2) was identified. We identified seven quantitative trait loci (QTLs) with additive effects. Among these, three QTLs on Chromosomes 7, 8, and 18 were resistant to both races. These QTLs could explain 1.91-7.73% of the phenotypic variation of SCN's female index. The QTLs on chromosomes 8 and 18 have already been reported and were most likely overlapped with rhg1 and Rhg4 loci, respectively. However, the QTL on chromosome 7 was novel. Candidate genes for the three QTLs were predicted through genes functional analysis and transcriptome analysis of infected roots of Handou 10 vs. Zheng 9525. Transcriptome analysis performed also indicated that the plant-pathogen interaction played an important role in the SCN resistance for Handou 10. The information will facilitate SCN-resistant gene cloning, and the novel resistant gene will be a source for improving soybeans' resistance to SCN.

SoySNP618K array: A high-resolution single nucleotide polymorphism platform as a valuable genomic resource for soybean genetics and breeding

Innovations in genomics have enabled the development of low-cost, high-resolution, single nucleotide polymorphism (SNP) genotyping arrays that accelerate breeding progress and support basic research in crop science. Here, we developed and validated the SoySNP618K array (618,888 SNPs) for the important crop soybean. The SNPs were selected from whole-genome resequencing data containing 2,214 diverse soybean accessions; 29.34% of the SNPs mapped to genic regions representing 86.85% of the 56,044 annotated high-confidence genes. Identity-by-state analyses of 318 soybeans revealed 17 redundant accessions, highlighting the potential of the SoySNP618K array in supporting gene bank management. The patterns of population stratification and genomic regions enriched through domestication were highly consistent with previous findings based on resequencing data, suggesting that the ascertainment bias in the SoySNP618K array was largely compensated for. Genome-wide association mapping in combination with reported quantitative trait loci enabled fine-mapping of genes known to influence flowering time, E2 and GmPRR3b, and of a new candidate gene, GmVIP5. Moreover, genomic prediction of flowering and maturity time in 502 recombinant inbred lines was highly accurate (>0.65). Thus, the SoySNP618K array is a valuable genomic tool that can be used to address many questions in applied breeding, germplasm management, and basic crop research.

Genome-Wide Association Studies of Soybean Yield-Related Hyperspectral Reflectance Bands Using Machine Learning-Mediated Data Integration Methods

In conjunction with big data analysis methods, plant omics technologies have provided scientists with cost-effective and promising tools for discovering genetic architectures of complex agronomic traits using large breeding populations. In recent years, there has been significant progress in plant phenomics and genomics approaches for generating reliable large datasets. However, selecting an appropriate data integration and analysis method to improve the efficiency of phenome-phenome and phenome-genome association studies is still a bottleneck. This study proposes a hyperspectral wide association study (HypWAS) approach as a phenome-phenome association analysis through a hierarchical data integration strategy to estimate the prediction power of hyperspectral reflectance bands in predicting soybean seed yield. Using HypWAS, five important hyperspectral reflectance bands in visible, red-edge, and near-infrared regions were identified significantly associated with seed yield. The phenome-genome association analysis of each tested hyperspectral reflectance band was performed using two conventional genome-wide association studies (GWAS) methods and a machine learning mediated GWAS based on the support vector regression (SVR) method. Using SVR-mediated GWAS, more relevant QTL with the physiological background of the tested hyperspectral reflectance bands were detected, supported by the functional annotation of candidate gene analyses. The results of this study have indicated the advantages of using hierarchical data integration strategy and advanced mathematical methods coupled with phenome-phenome and phenome-genome association analyses for a better understanding of the biology and genetic backgrounds of hyperspectral reflectance bands affecting soybean yield formation. The identified yield-related hyperspectral reflectance bands using HypWAS can be used as indirect selection criteria for selecting superior genotypes with improved yield genetic gains in large breeding populations.

Linkage and association study discovered loci and candidate genes for glycinin and β-conglycinin in soybean (Glycine max L. Merr.)

KEY MESSAGE

Linkage mapping and GWAS identified 67 QTLs related to soybean glycinin, β-conglycinin and relevant traits. Polymorphisms of the candidate gene Gy1 promoter were associated with the glycinin content in soybean. The major components of storage proteins in soybean seeds are glycinin and β-conglycinin, which play important roles in determining protein nutrition and soy food processing properties. Increasing the protein content while improving the ratio of glycinin to β-conglycinin is substantially important for soybean protein improvement. To investigate the genetic mechanism of storage protein subunits, 184 recombinant inbred lines (RILs) derived from a cross of Kefeng No. 1 and Nannong 1138-2 and 211 diverse soybean cultivars were used to detect loci related to glycinin (11S), β-conglycinin (7S), the sum of glycinin and β-conglycinin (SGC), and the ratio of glycinin to β-conglycinin (RGC). Sixty-seven QTLs and 11 hot genomic regions were identified as affecting the four traits. One genetic region (q10-1) on chromosome 10 was associated with multiple traits by both linkage and association analysis. Eight genes in 11 hot genomic regions might be related to soybean protein subunit. The candidate gene analysis showed that polymorphisms in Gy1 promoters were significantly correlated with the 11S content. The QTLs and candidate genes identified in the present study allow for further understanding the genetic basis of 11S and 7S regulation and provide useful information for marker-assisted selection (MAS) in soybean quality improvement.

Genome-wide association study and genomic selection for tolerance of soybean biomass to soybean cyst nematode infestation

Soybean cyst nematode (SCN), Heterodera glycines Ichinohe, is one of the most devastating pathogens affecting soybean production in the U.S. and worldwide. The use of SCN-resistant soybean cultivars is one of the most affordable strategies to cope with SCN infestation. Because of the limited sources of SCN resistance and changes in SCN virulence phenotypes, host resistance in current cultivars has increasingly been overcome by the pathogen. Host tolerance has been recognized as an additional tool to manage the SCN. The objectives of this study were to conduct a genome-wide association study (GWAS), to identify single nucleotide polymorphism (SNP) markers, and to perform a genomic selection (GS) study for SCN tolerance in soybean based on reduction in biomass. A total of 234 soybean genotypes (lines) were evaluated for their tolerance to SCN in greenhouse using four replicates. The tolerance index (TI = 100 × Biomass of a line in SCN infested / Biomass of the line without SCN) was used as phenotypic data of SCN tolerance. GWAS was conducted using a total of 3,782 high quality SNPs. GS was performed based upon the whole set of SNPs and the GWAS-derived SNPs, respectively. Results showed that (1) a large variation in soybean TI to SCN infection among the soybean genotypes was identified; (2) a total of 35, 21, and 6 SNPs were found to be associated with SCN tolerance using the models SMR, GLM (PCA), and MLM (PCA+K) with 6 SNPs overlapping between models; (3) GS accuracy was SNP set-, model-, and training population size-dependent; and (4) genes around Glyma.06G134900, Glyma.15G097500.1, Glyma.15G100900.3, Glyma.15G105400, Glyma.15G107200, and Glyma.19G121200.1 (Table 4). Glyma.06G134900, Glyma.15G097500.1, Glyma.15G100900.3, Glyma.15G105400, and Glyma.19G121200.1 are best candidates. To the best of our knowledge, this is the first report highlighting SNP markers associated with tolerance index based on biomass reduction under SCN infestation in soybean. This research opens a new approach to use SCN tolerance in soybean breeding and the SNP markers will provide a tool for breeders to select for SCN tolerance.

Molecular Basis of Soybean Resistance to Soybean Aphids and Soybean Cyst Nematodes

Soybean aphid (SBA; Aphis glycines Matsumura) and soybean cyst nematode (SCN; Heterodera glycines Ichninohe) are major pests of the soybean (Glycine max [L.] Merr.). Substantial progress has been made in identifying the genetic basis of limiting these pests in both model and non-model plant systems. Classical linkage mapping and genome-wide association studies (GWAS) have identified major and minor quantitative trait loci (QTLs) in soybean. Studies on interactions of SBA and SCN effectors with host proteins have identified molecular cues in various signaling pathways, including those involved in plant disease resistance and phytohormone regulations. In this paper, we review the molecular basis of soybean resistance to SBA and SCN, and we provide a synthesis of recent studies of soybean QTLs/genes that could mitigate the effects of virulent SBA and SCN populations. We also review relevant studies of aphid–nematode interactions, particularly in the soybean–SBA–SCN system.

Genome-Wide Association Analysis Pinpoints Additional Major Genomic Regions Conferring Resistance to Soybean Cyst Nematode (Heterodera glycines Ichinohe)

Soybean cyst nematode (Heterodera glycines Ichinohe) (SCN) is the most destructive pest affecting soybeans [Glycine max (L.) Merr.] in the U.S. To date, only two major SCN resistance alleles, rhg1 and Rhg4, identified in PI 88788 (rhg1) and Peking (rhg1/Rhg4), residing on chromosomes (Chr) 18 and 8, respectively, have been widely used to develop SCN resistant cultivars in the U.S. Thus, some SCN populations have evolved to overcome the PI 88788 and Peking derived resistance, making it a priority for breeders to identify new alleles and sources of SCN resistance. Toward that end, 461 soybean accessions from various origins were screened using a greenhouse SCN bioassay and genotyped with Illumina SoySNP50K iSelect BeadChips and three KASP SNP markers developed at the Rhg1 and Rhg4 loci to perform a genome-wide association study (GWAS) and a haplotype analysis at the Rhg1 and Rhg4 loci. In total, 35,820 SNPs were used for GWAS, which identified 12 SNPs at four genomic regions on Chrs 7, 8, 10, and 18 that were significantly associated with SCN resistance (P < 0.001). Of those, three SNPs were located at Rhg1 and Rhg4, and 24 predicted genes were found near the significant SNPs on Chrs 7 and 10. KASP SNP genotyping results of the 462 accessions at the Rhg1 and Rhg4 loci identified 30 that carried PI 88788-type resistance, 50 that carried Peking-type resistance, and 58 that carried neither the Peking-type nor the PI 88788-type resistance alleles, indicating they may possess novel SCN resistance alleles. By using two subsets of SNPs near the Rhg1 and Rhg4 loci obtained from SoySNP iSelect BeadChips, a haplotype analysis of 461 accessions grouped those 58 accessions differently from the accessions carrying Peking or PI 88788 derived resistance, thereby validating the genotyping results at Rhg1 and Rhg4. The significant SNPs, candidate genes, and newly characterized SCN resistant accessions will be beneficial for the development of DNA markers to be used for marker-assisted breeding and developing soybean cultivars carrying novel sources of SCN resistance.

Genome wide association study discovers genomic regions involved in resistance to soybean cyst nematode (Heterodera glycines) in common bean

Common bean (Phaseolus vulgaris L.) is an important high protein crop grown worldwide. North Dakota and Minnesota are the largest producers of common beans in the USA, but crop production is threatened by soybean cyst nematode (SCN; Heterodera glycines Ichinohe) because most current cultivars are susceptible. Greenhouse screening data using SCN HG type 0 from 317 plant introductions (PI's) from the USDA core collection was used to conduct a genome wide association study (GWAS). These lines were divided into two subpopulations based on principal component analysis (Middle American vs. Andean). Phenotypic results based on the female index showed that accessions could be classified as highly resistant (21% and 27%), moderately resistant (51% and 48%), moderately susceptible (27% and 22%) and highly susceptible (1% and 3%) for Middle American and Andean gene pools, respectively. Mixed models with two principal components (PCs) and kinship matrix for Middle American genotypes and Andean genotypes were used in the GWAS analysis using 3,985 and 4,811 single nucleotide polymorphic (SNP) markers, respectively which were evenly distributed across all 11 chromosomes. Significant peaks on Pv07, and Pv11 in Middle American and on Pv07, Pv08, Pv09 and Pv11 in Andean group were found to be associated with SCN resistance. Homologs of soybean rhg1, a locus which confers resistance to SCN in soybean, were identified on chromosomes Pv01 and Pv08 in the Middle American and Andean gene pools, respectively. These genomic regions may be the key to develop SCN-resistant common bean cultivars.

Genome-Wide Analysis of Tar Spot Complex Resistance in Maize Using Genotyping-by-Sequencing SNPs and Whole-Genome Prediction

Tar spot complex (TSC) is one of the most destructive foliar diseases of maize ( L.) in tropical and subtropical areas of Central and South America, causing significant grain yield losses when weather conditions are conducive. To dissect the genetic architecture of TSC resistance in maize, association mapping, in conjunction with linkage mapping, was conducted on an association-mapping panel and three biparental doubled-haploid (DH) populations using genotyping-by-sequencing (GBS) single-nucleotide polymorphisms (SNPs). Association mapping revealed four quantitative trait loci (QTL) on chromosome 2, 3, 7, and 8. All the QTL, except for the one on chromosome 3, were further validated by linkage mapping in different genetic backgrounds. Additional QTL were identified by linkage mapping alone. A major QTL located on bin 8.03 was consistently detected with the largest phenotypic explained variation: 13% in association-mapping analysis and 13.18 to 43.31% in linkage-mapping analysis. These results indicated that TSC resistance in maize was controlled by a major QTL located on bin 8.03 and several minor QTL with smaller effects on other chromosomes. Genomic prediction results showed moderate-to-high prediction accuracies in different populations using various training population sizes and marker densities. Prediction accuracy of TSC resistance was >0.50 when half of the population was included into the training set and 500 to 1,000 SNPs were used for prediction. Information obtained from this study can be used for developing functional molecular markers for marker-assisted selection (MAS) and for implementing genomic selection (GS) to improve TSC resistance in tropical maize.