Using NGS Technology and Association Mapping to Identify Candidate Genes Associated with Fusarium Stalk Rot Resistance

Stalk rot caused by Fusarium fungi is one of the most widespread and devastating diseases of maize, and the introduction of resistant genotypes is one of the most effective strategies for controlling the disease. Breeding genotypes with genetically determined resistance will also allow less use of crop protection products. The aim of the research was to identify molecular markers and associated candidate genes determining maize plant resistance to Fusarium stalk rot. The plant material for this study consisted of 122 maize hybrids. The experiment was conducted in two localities: Smolice and Kobierzyce. The Fusarium stalk rot values ranged from 1.65% (for genotype G01.10) to 31.18% (for genotype G03.07) in Kobierzyce and from 0.00% (for 58 genotypes) to 6.36% (G05.03) in Smolice. The analyzed genotypes were simultaneously subjected to next-generation sequencing using the Illumina platform. Illumina sequencing identified 60,436 SilicoDArT markers and 32,178 SNP markers (92,614 in total). For association mapping, 32,900 markers (26,234 SilicoDArT and 6666 SNP) meeting the criteria (MAF > 0.25 and the number of missing observations <10%) were used. The results of the observation of the degree of infection and sequencing were used for association mapping, which ultimately resulted in the selection of ten molecular markers important at both places. Among the identified markers, two SNP markers that are located inside candidate genes play an important role. Marker 4772836 is located inside the serine/threonine-protein kinase bsk3 gene, while marker 4765764 is located inside the histidine kinase 1 gene. Both genes can be associated with plant resistance to Fusarium stalk rot, and these genes can also be used in breeding programs to select resistant varieties.


Introduction
In the last decade, the increase in temperature, expansion of areas, intensification of cultivation, introduction of agrotechnical simplifications, and the appearance of new pest species have resulted in a significant increase in the threat to the amount and quality of maize yields [1].Crops of this cereal are frequently infected by pathogenic fungi, resulting in serious yield losses [2].Currently, fungi of the genus Fusarium can be classified as the most dangerous pathogens causing serious losses in maize plantations.The fungus Fusarium graminearum causes stalk rot in maize [3,4].This pathogen not only causes yield losses, but also poses a huge threat related to the presence of mycotoxins [5][6][7].Fungi of the genus Fusarium colonize soil and crop residues of maize and cereals, and are also present on the surface of grains.They reach plants directly from resting places and through wind, water, soil fragments (dust), insects, and animals.Maize is exposed to Fusarium fungi from the moment the grain is sown [8].As a result of infection, maize plants often lodge, leading to a decline in planting density.
Maize plant lodging caused by Fusarium fungus infection is often observed in fields where seed is sown untreated with a fungicide.This causes large losses in the quantity and quality of the crop.For some genotypes, these losses are even 100 percent.The development of the pathogen can also be favored by weather factors, such as a cool and rainy spring, heavy rainfall, or sowing too deeply in heavy soil [9].Infection can be prevented; several studies have shown that soil fertilization has a positive effect on plant tolerance to pathogenic fungi, particularly the use of fertilizers containing zinc or potassium phosphite [10][11][12].Stalk rot usually occurs when plant tissues are damaged by Ostrinia nubilalis.The feeding of this pest is largely responsible for the development of the disease, especially since O. nubilalis is also a vector of this fungus and, migrating between plants, transfers fungal spores from diseased plants to healthy ones [13].There are usually two generations of O. nubilalis larvae per year: the first generation attacks plants in the middle or end of the vegetative phase, and the second generation attacks plants in the reproductive phase (from early milk stage to maturity) [14].Strong decomposition of stem tissues, including nodes, leads to whole plants falling over, making harvesting difficult and sometimes impossible.
Fighting Fusarium diseases that appear at the end of the growing season is difficult, especially when no treatments are applied beforehand to directly or indirectly limit infection [15].Currently, the most effective method of combating this pathogen is the cultivation of varieties with genetically determined resistance to F. graminearum.A dynamically changing climate and an increase in demand for maize, driven by a natural increase in population, are catalysts for research into genes and genomic regions important to agronomy [16,17].Scientists have long used next-generation sequencing (NGS) and association mapping methods to identify markers associated with resistance genes.Genome-wide association studies (GWAS) are a useful tool for identifying candidate genes, especially when combined with QTL mapping to validate loci for quantitative traits.Zila et al. [18,19] performed GWAS tests in maize to detect SNPs associated with increased resistance to Fusarium.They identified 10 SNP markers significantly associated with resistance to this pathogen.Many previous studies have been conducted with GWAS and SNPs on Fusarium stalk rot.Liu et at.[20], based on 165 maize inbred lines and 4666 SNP markers, performed genetic diversity analysis, GWAS, linkage analysis, and phenotypic analysis through the inoculation of stalk rot pathogen on maize plants in the field.For identification and validation of genomic regions associated with Fusarium stalk rot resistance, Rashid et al. [21] used GWAS for 342 maize lines and 290,626 SNPs.The panel was screened for Fusarium stalk rot in three environments using standard artificial inoculation methodology.In GWAS conducted by Stagnati et al. [22] on 230 inbred lines, some of the SNPs identified for seedling rot co-localized within QTL intervals previously identified for Fusarium seed rot, Fusarium ear rot, and fumonisin accumulation in maize.The introduction of molecular analyses made it possible to develop a methodology for genetic marker-assisted selection (MAS).The scope of application of this methodology depends on the progress of knowledge of the genome of the species [23,24].The completion of the sequencing of the human genome has opened up a wide range of possibilities for understanding the genomic sequences of crop plants [25][26][27][28].For several years, maize breeding around the world has been supported by useful molecular markers.Many authors state in their publications that breeding supported by molecular marker accelerates yield increases not only in the United States, but also in other countries, offering great potential to increase the productivity and value of maize germplasm [29,30].
Therefore, the aim of this study was to identify new SNP and SilicoDArT markers linked to candidate genes for maize stalk rot resistance.

Field Experiment
The experiment was established in plots of 10 m 2 , in three repetitions, in a randomized complete block design, at two localities (Smolice (51

Weather Conditions
In 2022, the average rainfall in Smolice was 34.45 mm, which was 13.82 mm lower than the multi-year average.The highest rainfall was recorded in July (55 mm) and the lowest in March (15 mm).The average air temperature this year in Smolice was 9.54 • C, higher than the long-term average temperature by 0.8 • C. The warmest month in 2004 was August (20 • C), while the lowest temperature was recorded in December (1.1 • C).In 2022, the amount of rainfall and temperature were unfavorable during the initial development of maize plants.Despite the early sowing date, maize remained at the 2-3 leaf phase for a long time, during which purple discoloration was visible due to difficulties in taking up phosphorus from the soil.May had an abundance of rainfall, which had a positive effect on maize development.In Kobierzyce, the average rainfall in 2022 was 51.52 mm, 3.22 mm higher than the average multi-year rainfall.The highest rainfall was recorded in May, with 95.8 mm, and the lowest in March, with 21.6 mm.The average air temperature in Kobierzyce was 11.46 • C, which was higher than the multi-year average air temperature by 2.58 • C. The warmest month in 2022 in Kobierzyce, as in Smolice, was August (22 • C), and the coldest was December (1.5 • C).

Determining the Degree of Fusarium Infection
Infection of maize plants by fungi Fusarium verticillioides was observed in natural conditions, in field experiments that were conducted in two localities in three repetitions.The plants were not artificially inoculated.In the stage of full maize maturity (BBCH 86), shortly before harvest, the number of plants with Fusarium wilt symptoms was determined, distinguishing the base of the stalk rot.The number of plants with symptoms of stalk rot was determined by lightly hitting the base of the stem with a foot, which caused fractures of the infected plants [32].The degree of plant infection with stalk rot caused by fungi of the genus Fusarium spp. was determined during maize harvest based on visual inspection of the longitudinal section of the three lower internodes of the plant.The observation was made on ten plants.A 9-point scale was used [33], adopting the following parameters for individual degrees (even numbers mean indirect infection): 1-no visible symptoms (healthy plant), 3-one or two nodes turn brown and decay, 5-one or two internodes brownish, decaying, 7-more than two internodes brown and rotting, but pith tissue still exists around the vascular bundles, 9-all tissues decomposed, only conductive bundles visible.Next, for each repetition, based on observations of ten plants, the average percentage of stalk rot infected with Fusarium species was calculated.

DNA Isolation
DNA isolation from 122 maize hybrids was carried out using a ready-made set of reagents from Symbios.The isolated DNA samples were subjected to next-generation sequencing.The concentration and purity of the isolated DNA samples were determined using a DS-11 spectrophotometer from DeNovix.The isolated DNA matrix was adjusted to a uniform concentration of 100 ng µL −1 by dilution with deionized distilled water.

Genotyping
DArTseq technology, based on next-generation sequencing, was used for genotyping.Isolated DNA from the 122 maize hybrids tested at 50 µL at a concentration of 100 ng was transferred on two 96-well Eppendorf plates for analyses to identify SilicoDArT and SNP polymorphisms.The analyses were performed at Diversity Arrays Technology, University of Canberra, Australia.The methods used are described in detail on the Diversity Arrays Technology website: (https://www.diversityarrays.com/technology-and-resources/dartseq/, accessed on 15 January 2024).

Association Mapping Using GWAS Analysis
Based on the results obtained from genotyping and phenotyping, association mapping was performed using GWAS analysis.Associations were determined for 122 analyzed maize hybrids.Genotypic data were obtained from DArTseq analysis, while phenotypic data were determined from observations of the percentage of stem infection and lodging of maize.Based on GWAS analysis, SilicoDArT and SNP markers with the highest level of significance, i.e., those that were most strongly associated with plant resistance to Fusarium stalk rot, were selected for further research.For the association analysis, only SilicoDArT and SNP sequences meeting the following criteria were selected: one SilicoDArT and/or SNP within a given sequence (69 nt), minor allele frequency (MAF) > 0.25, and the missing observation fractions <10%.

Statistical Analysis and Association Mapping
Data were analyzed using analysis of variance in a model with fixed effects of location and random effects of genotypes and genotype × location interaction.Association mapping, based on SilicoDArT and SNP data and average trait values, was conducted separately for data from the two locations using the method based on the mixed linear model with a population structure estimated by eigenanalysis and modelled by random effects [34,35].All analyses and visualizations of the results were performed in GenStat 23.1 [36], using procedure QSASSOCIATION.QSASSOCIATION performs a mixed model marker-trait association analysis (also known as linkage disequilibrium mapping) with data from a single-environment trial.To avoid false positives in association mapping studies, some form of control is necessary for the genetic relatedness.The model used was specified by the RELATIONSHIPMODEL=eigenanalysis option, which infers the underlying genetic substructure in the population by retaining the most significant principal components from the molecular marker matrix [37]; the scores of the significant axes are used as covariables in the mixed model, which effectively is an approximation to the structuring of the genetic variance covariance matrix by a coefficient of coancestry matrix (kinship matrix).Significance of association between Fusarium stalk rot and SilicoDArT and SNP markers was assessed using p values corrected for multiple testing using the Benjamini-Hochberg method.

Physical Mapping
SilicoDArT and SNP marker sequences selected based on GWAS analysis were subjected to basic local alignment search tool (BLAST) analysis, which involves searching databases to find sequences with high homology to the selected silicoDArT and SNP marker sequences.Publicly available web browsers were used for this purpose, specifically CEPH Genotype database http://www.cephb.fr/en/cephdb/(accessed on 15 January 2024), NCBI Map Viewer http://www.ncbi.nlm.nih.gov/projects/mapview/(accessed on 15 January 2024), UCSC Genome Browser http://genome.ucsc.edu/(accessed on 15 January 2024), and Ensembl Map View https://www.ensembl.org/index.html(accessed on 15 January 2024).The programs used were applied to indicate the chromosomal locations of the searched sequences, similar to the sequences analyzed, and to determine their physical location.The sequences of all genes located in the designated area on the chromosome were further analyzed.

Functional Analysis of Gene Sequences
Functional analysis was carried out using Blast2GO 6.0 software.The sequences of all genes located in the chromosome regions determined by BLAST analysis were analyzed.The goal was to obtain information on the biological function of gene sequences located in the designated chromosome region.

Design of Primers for Identified SilicoDArT and SNP Polymorphisms Related to Yield and Its Characteristics
The Primer 3 program was used to design the primers.The designed primers were used for a conventional PCR reaction; therefore, the PCR product size range from 50 to 500 bp was selected in the program.The melting temperature of the primers was selected as follows: min.57.0 • C, opt.60.0 • C, max.63.0 • C, max.Tm difference 3. The other settings were not changed.

Symptoms of Fusarium Stalk Rot
Typical symptoms of Fusarium stalk rot were observed on individual plants: rotting at roots, crown, and lower internodes; wilting of the plants, taking on a grayish-green hue, then turning tan; and light pink to tan discoloration, indicating splitting of the inner stalk (Figure 1).These symptoms were observed to varying degrees depending on the hybrid.Significantly higher Fusarium stalk rot infestation was observed in Kobierzyce than in Smolice.

Phenotyping
Detailed information on phenotypes such as mean, skewness, kurtosis, standard error, and coefficient of variation observed in each locality, as well as broad-sense heritability, is presented in Table 1.Significantly higher Fusarium stalk rot infestation was observed in Kobierzyce than in Smolice (Table 1).The large number of completely resistant

Phenotyping
Detailed information on phenotypes such as mean, skewness, kurtosis, standard error, and coefficient of variation observed in each locality, as well as broad-sense heritability, is presented in Table 1.Significantly higher Fusarium stalk rot infestation was observed in Kobierzyce than in Smolice (Table 1).The large number of completely resistant plants (Fusarium stalk rot values equal to zero) in the experiment conducted in Smolice resulted in higher skewness, kurtosis, and coefficients of variation than in the experiment conducted in Kobierzyce.Broad-sense heritability for Fusarium stalk rot estimated from environments was 66.1 (Table 1).The observed Fusarium stalk rot was characterized by a normal distribution about both localities, which made it possible to carry out association mapping.Analysis of variance indicated that the main effects of genotype and location, as well as genotype × location interaction were significant for Fusarium stalk rot (Table 2).The differences between the locations were large (11.39% in Kobierzyce and 1.20% in Smolice), which was the result of large differences in weather conditions between the two locations (Figure 2).The Fusarium stalk rot values ranged from 1.65% (for genotype G01.10) to 31.18% (for genotype G03.07) in Kobierzyce and from 0.00% (for 58 genotypes) to 6.36% (G05.03) in Smolice.On average for both localities, the Fusarium stalk rot values ranged from 0.82% (for genotype G01.10) to 17.08% (for genotype G03.07) (Table S1).The variability of the observed values (expressed in standard deviations) also differed between locations, from 0% to 16.65% in Kobierzyce and from 0% to 9.45% in Smolice (Table S1).

Genotype
Despite very large differences between Fusarium stalk rot values, their correlation was observed between Kobierzyce and Smolice (r = 0.36, p < 0.001).

Genotyping
Illumina sequencing identified 60,436 SilicoDArT markers and 32,178 SNP markers (92,614 in total).For association mapping, 32,900 markers (26,234 SilicoDArT and 6666 SNP) meeting the criteria (MAF > 0.25 and the number of missing observations <10%) were used.Based on the identified SilicoDArT and SNP molecular markers, a dendrogram of genetic similarity between the 122 hybrids analyzed was prepared.The dendrogram distinguished three groups of similarity.Group I contained 53 hybrids, Group II contained 50 hybrids, and Group III contained 19 hybrids (Figure 3).
SNP) meeting the criteria (MAF > 0.25 and the number of missing observations <10%) w used.Based on the identified SilicoDArT and SNP molecular markers, a dendrogram genetic similarity between the 122 hybrids analyzed was prepared.The dendrogram d tinguished three groups of similarity.Group I contained 53 hybrids, Group II contain 50 hybrids, and Group III contained 19 hybrids (Figure 3).The association mapping identified 11,282 markers (8405 SilicoDArT and 2877 SN in Smolice (Table 3, Figure 4) and 9439 markers (7148 SilicoDArT and 2291 SNP) Kobierzyce (Table 3, Figure 5) that were significant at a level of 0.05.The percentage variation explained by each marker ranged from 2.4% to 34.4% in Smolice and from 2 to 27.0% in Kobierzyce (Table 3).A total of 8702 markers were statistically significant s ultaneously in both localities.Of these, 91 were highly statistically significant (LOD >  Of the 8702 physically mapped statistically significant markers in both localities (Smolice and Kobierzyce), 10 were significantly associated with plant resistance to Fusarium stalk rot at both locations (Table 5).These were the markers with the highest statistical significance in both localities (LOD > 9.5 in Kobierzyce, LOD > 7.5 in Smolice) and the highest coefficients of determination (R 2 •100 > 27.4% in Kobierzyce, R 2 •100 > 22.1% in Smolice).Using the BLAST database, the location of the selected markers was determined and the associated candidate genes were provided.The next step was to design primers that would be used to identify the 10 selected markers.The sequences of the primers are shown in Table S2.The primers were designed to identify the analyzed markers using the conventional PCR method.Among the identified markers, two SNP markers located inside candidate genes play an important role.Marker 4772836 is located inside the serine/threonine-protein kinase bsk3 gene, while marker 4765764 is located inside the histidine kinase 1 gene.Both genes may be associated with plant resistance to Fusarium stalk rot.

Discussion
Maize is currently a leading crop in the global grain market, behind only rice and wheat [38].Its largest producers are the United States and China, with 2018 production of 392 million tons and 257 million tons, respectively [39].In terms of yield growth rate in production, maize dominates among cereal plants, with yield potential reaching 12-15 tons of dry matter per 1 ha.The high fertility of maize is the result of many favorable physiological and morphological traits of this species [40].However, it should be remembered that yield levels are affected by weather conditions during a given growing season and local environmental conditions [41].Maize has successfully acclimatized to Polish soil and climatic conditions [42,43].The reason for this phenomenon is the introduction of linear hybrids into maize cultivation, providing access to varieties with an appropriate growing season.Maize, as a crop, is of great importance in the world, both in terms of utility and economy.This is due to the possibility of using virtually all of the plant's biomass as feed, food, or industrial raw material [39].
Plants, including maize, are constantly exposed to many stress factors, both abiotic and biotic, that can affect their growth and development [41].Biotic stress factors include viruses, bacteria, fungi, and insect pests.Often the consequence of stress factors is the occurrence of pathogenic diseases, which can interfere with normal plant growth, development, and function, and cause significant yield reduction [44].Maize diseases occurring worldwide include common smut of maize, maize small leaf spot, maize sheath blight, common rust of maize, maize Helminthosporium, and maize stalk rot [45].
Fusarium stalk rot of maize is caused by fungi of the genus Fusarium.Damage caused by European maize stalk rot promotes the development of the pathogenic fungus [46].Characteristic traits of the disease are wilting and drying of leaf blades, weakening of the stem, and often lodging [47].In addition, Fusarium fungi increase the content of mycotoxins, which in turn disqualify grain from further processing [48][49][50].
In our own research, the field experiment made it possible to perform and analyze measurements of the degree of Fusarium infection of maize stalks.The observed trait was characterized by a normal distribution at both localities, which made it possible to perform association mapping.Analysis of variance indicated that the main effects of genotype and location, as well as genotype and location interaction, were significant for Fusarium stalk rot.
The most important mycotoxins include aflatoxins, fumonisins, ochratoxin A, deoxynivalenol, zearalenone, and ergot alkaloids, which are mainly produced by the genera Aspergillus, Penicillium, Fusarium, and Claviceps [6].Selection of resistant genotypes based on their phenotype can be unreliable and is time-and labor-consuming.Phenotyping maize for Fusarium resistance requires field trials conducted in several environments due to the high interaction of genotype and environment [51].In addition, proper management of inoculation timing is required due to the diversity of genotypes in terms of flowering [52].Currently, the most effective and fastest methods to support selection are transcriptomic and genomic association studies (GWAS).These are useful tools for identifying candidate genes, especially when combined with QTL mapping to map and validate loci for quantitative traits [53].Combining these methods has allowed us to overcome their individual limitations [54,55].
Our study identified 60,436 SilicoDArT markers and 32,178 SNP markers (92,614 in total).For association mapping, 32,900 markers (26,234 silicoDArT and 6666 SNP) were selected, meeting assumptions about the distribution of values (MAF > 0.25) and the number of missing observations (<10%).The differentiation of hybrids at the genetic level makes the association mapping results obtained even more reliable and will be solid selection tools in maize breeding.The association mapping identified 11,282 markers (8405 silicoDArT and 2877 SNP) in the Smolice experiment and 9439 markers (7148 silicoDArT and 2291 SNP) in the Kobierzyce experiment.Of the 8702 markers associated with physical mapping, 10 were selected that were statistically significant in both localities (Smolice and Kobierzyce).The identification of ten markers coupled to specific genes with specific functions provides information about the genetic basis, as well as pointing to the molecular mechanisms controlling resistance to Fusarium stalk rot.The use of maize hybrids as an association population makes it possible to create favorable conditions for increasing maize resistance to F. verticillioides through genomics-assisted breeding [56,57].
In recent years, many authors have attempted to identify molecular markers associated with functionally important traits in maize.Using next-generation sequencing, Sobiech et al. [58] identified markers associated with Fusarium resistance in maize plants.Bocianowski et al. [59] used NGS technology and association mapping to identify markers associated with the heterosis effect in maize.A similar study was conducted by Tomkowiak et al. [60], who identified six SNP markers (1818, 14506, 2317, 3233, 11657, and 12812) located within genes, on chromosomes 8, 9, 7, 3, 5, and 1, respectively, associated with maize yield.Other authors [61] identified four genes (the isoform × 2 gene of sucrose synthase 4, the isoform × 1 gene of phosphoinositide phosphatase, the putative isoform × 1 gene of the SET domain-containing protein family, and the grx_c8-glutaredoxin subgroup iii gene) which significantly regulate the level of vigor and germination of maize seeds.
Stalk rot is a destructive disease in maize and occurs worldwide.The genetic mechanism of maize stalk rot is complex.It is caused by one or more fungal pathogen(s).The causative agents, including F. verticillioides, Fusarium moniliforme, and Fusarium proliferatum, are soil-borne fungus, meaning that it is difficult to control stalk rot.Thus, it is necessary to understand the genetic basis of stalk rot and to breed resistant varieties.The stalk rot resistance in maize is controlled by polygene and the additive genetic effects among the genes [62,63].Until now, a number of quantitative trait loci (QTLs) associated with maize stalk rot have been identified [64][65][66].In the study presented here, of the ten markers selected, two of them are located inside genes.SNP marker 4772836, located on chromosome 2, is located inside the serine/threonine-protein kinase bsk3 gene, while SNP marker 4765764, located on chromosome 5, is inside the histidine kinase 1 gene.Histidine kinases are typically trans-membrane proteins of the transferase class that play a role in signal transduction across the cellular membrane.The vast majority of histidine kinases (HKs) are homodimers that exhibit autokinase, phosphotransfer, and phosphatase activity [67].Unlike other classes of protein kinases, HKs are usually parts of a two-component signal transduction mechanism, in which the HK transfers a phosphate group from ATP to a histidine residue inside the kinase, then to an aspartate residue in the receiver domain of the response regulator protein (or sometimes on the kinase itself).Histidine kinases have been implicated in cytokinin and ethylene signaling, in regulation responses to environmental stress [68].It can be speculated that they are also involved in the immune response associated with Fusarium infection of plants.
The second gene analyzed, a probable serine/threonine kinase, acts as a positive regulator of brassinosteroid (BR) signaling.Brassinosteroids, another class of growth-promoting regulator, regulate many aspects of plant growth and development, including cell expansion and elongation (in combination with auxin) and vascular differentiation [69][70][71][72][73].The necessity of brassinosteroids in plant growth and development has been demonstrated by the identification of many BR-deficient dwarf mutants and subsequent BRs treating experiments [74].Brassinosteroid-signaling kinases (BSKs) are critical in BRs signal transduction.BSK3 is the only BSK member involved in BR-mediated plant root growth [75].
To confirm the involvement of the two candidate genes in the immune response to infection of maize plants by Fusarium fungi, their expression will be analyzed.If the involvement of these genes in the immune reaction is confirmed, SNP markers associated with them can be used to select resistant varieties.

Conclusions
In maize cultivation, it is crucial to reduce the severity of diseases that cause a drastic decrease in yields.Excessive use of chemical protection negatively affects the natural environment, so the best method of crop protection is to use varieties with genetically determined resistance in breeding.Understanding the genetic basis and molecular mechanisms controlling Fusarium resistance is a key requirement for developing maize varieties with enhanced resistance.Because resistance to Fusarium fungi is a quantitative trait based on a distributed architecture of multiple genes with small effectors, the best approach for future molecular breeding is MAS.Using next-generation sequencing (NGS), two significant markers were identified in candidate genes associated with Fusarium resistance in maize plants.SNP marker 4772836, located on chromosome 2, is inside the serine/threonineprotein kinase bsk3 gene, while SNP marker 4765764, located on chromosome 5, is inside the histidine kinase 1 gene.If the involvement of these genes in the immune response is confirmed, SNP markers associated with them can be used to select resistant varieties.

Figure 2 .
Figure 2. Density charts showing the distribution of the Fusarium stalk rot.

Figure 2 .
Figure 2. Density charts showing the distribution of the Fusarium stalk rot.

Figure 3 .
Figure 3. Dendrogram showing the genetic similarity between the 122 analyzed hybrids construc based on the observation of 32,900 markers (26,234 SilicoDArT and 6666 SNP).

Figure 3 .
Figure 3. Dendrogram showing the genetic similarity between the 122 analyzed hybrids constructed based on the observation of 32,900 markers (26,234 SilicoDArT and 6666 SNP).

Table 1 .
Statistical summary of phenotypes of 122 maize hybrids for Fusarium stalk rot.

Table 2 .
Results of two-way analysis of variance for Fusarium stalk rot.

Table 3 .
SilicoDArT and SNP molecular markers significantly associated with maize plant resistance to Fusarium stalk rot in Smolice and Kobierzyce (significant associations selected at p < 0.001 with correction for Benjamini-Hochberg multiple testing).

Table 5 .
Characteristics and location of markers significantly related to plant resistance to Fusarium stalk rot.