Next Article in Journal
Analysis of DNA Methylation Differences during the JIII Formation of Bursaphelenchus xylophilus
Next Article in Special Issue
Overexpression of KvCHX Enhances Salt Tolerance in Arabidopsis thaliana Seedlings
Previous Article in Journal
Integrative Analysis of the Role of TP53 in Human Pan-Cancer
Previous Article in Special Issue
Exploring Plant Meiosis: Insights from the Kinetochore Perspective
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Identification of Multiple Genetic Loci Related to Low-Temperature Tolerance during Germination in Maize (Zea maize L.) through a Genome-Wide Association Study

1
Maize Research Institute of Heilongjiang Academy of Agricultural Sciences, Harbin 150086, China
2
Key Laboratory of Biology and Genetics Improvement of Maize in Northern Northeast Region, Ministry of Agriculture and Rural Affairs, Harbin 150086, China
3
Key Laboratory of Germplasm Resources Creation and Utilization of Maize, Harbin 150086, China
4
Keshan Branch of Heilongjiang Academy of Agricultural Sciences, Qiqihaer 161000, China
*
Authors to whom correspondence should be addressed.
These authors contributed equally to this work.
Curr. Issues Mol. Biol. 2023, 45(12), 9634-9655; https://doi.org/10.3390/cimb45120602
Submission received: 25 October 2023 / Revised: 13 November 2023 / Accepted: 21 November 2023 / Published: 29 November 2023
(This article belongs to the Special Issue Molecular Breeding and Genetics Research in Plants)

Abstract

:
Low-temperature stress during the germination stage is an important abiotic stress that affects the growth and development of northern spring maize and seriously restricts maize yield and quality. Although some quantitative trait locis (QTLs) related to low-temperature tolerance in maize have been detected, only a few can be commonly detected, and the QTL intervals are large, indicating that low-temperature tolerance is a complex trait that requires more in-depth research. In this study, 296 excellent inbred lines from domestic and foreign origins (America and Europe) were used as the study materials, and a low-coverage resequencing method was employed for genome sequencing. Five phenotypic traits related to low-temperature tolerance were used to assess the genetic diversity of maize through a genome-wide association study (GWAS). A total of 14 SNPs significantly associated with low-temperature tolerance were detected (−log10(P) > 4), and an SNP consistently linked to low-temperature tolerance in the field and indoors during germination was utilized as a marker. This SNP, 14,070, was located on chromosome 5 at position 2,205,723, which explained 4.84–9.68% of the phenotypic variation. The aim of this study was to enrich the genetic theory of low-temperature tolerance in maize and provide support for the innovation of low-temperature tolerance resources and the breeding of new varieties.

1. Introduction

The northern spring corn area is an important corn production area and commercial grain base in China, located at the northern end of China’s golden corn belt. However, due to the special geographical location and environmental conditions, low temperatures in spring are an important source of non-biotic stress that affects the seedling quality in this area, which seriously restricts the yield and quality of the corn produced. The low-temperature tolerance of maize belongs to a quantitative trait controlled by multiple genes. In recent years, with the development of molecular biology, scholars have carried out quantitative trait locis (QTL) analyses on its low-temperature tolerance, locating maize’s low-temperature tolerance on chromosomes 1–10. One QTL was located in an interval on chromosome 6, which was associated with three low-temperature tolerance traits and could explain 18.1–32.8% of the phenotypic variation [1]; twenty-six QTLs, associated with seed vigor, were detected under low temperatures during maize’s germination stage on chromosomes 2, 3, 5, and 9, alongside five meta-QTLs [2]. In the 176 IBM Syn10 doubled-haploid population from the B73 × Mo17 cross, there were thirteen QTLs associated with a low-temperature germination ability, three B73 upregulated genes, and five Mo17 upregulated genes found by combining the RNA-Seq technology and QTL analysis [3]. A recombinant inbred line population (IBM Syn4 RIL) from a B73 and Mo17 cross was used to identify QTLs and investigate the genetic architecture under low-temperature conditions at a young seedling stage: two QTLs (bin 1.02 and bin 5.05) with a high additive impact were detected, which were associated with cold tolerance [4]. A total of 406 recombinant inbred lines from a multi-parent, advanced-generation, intercross population were used and, as a result, many cold tolerance-related traits were recorded: the 858 SNPs were found that were significantly associated with all traits, which indicated that most QTLs are related to chlorophyll and Fv/Fm; the authors also located most of the QTLs in specific regions, particularly bin 10.04 [5]. An F2 population was constructed from the cross of IB030 and Mo17 to map QTLs associated with cold tolerance via QTL-seq and transcriptomic integrative analyses, and two positively regulated genes (ZmbZIP113 and ZmTSAH1) that control the low-temperature germination ability were identified [6]. Scholars performed QTL mapping on an IBM (intermated B73 × Mo17) Syn10 doubled-haploid (DH) population, and twenty-eight QTLs that significantly correlated with low-temperature germination were detected, and these QTLs explained 5.4–13.34% of the phenotypic variation. In addition, six QTL clusters were produced by fourteen overlapping QTLs on every chromosome, except for chromosomes 8 and 10 [7]. The identification of molecular marker loci associated with QTLs or genes can contribute to the study of the cold-tolerance mechanism of maize and could be further used for breeding cold-tolerant inbred lines or hybrids. The QTLs controlling low-temperature tolerance during the germination stage are distributed on chromosomes 1–10, and there are few QTLs that have been consistently identified using different methods and materials, with large intervals. At the same time, there is more than one main QTL interval, so it is necessary to continue to mine consistent, main QTLs and identify candidate genes.
In recent years, the construction of reference genomes, such as B73, MO17, W22, PH207, and CML247, has enabled the widespread application of high-throughput single-nucleotide sequence markers, greatly improving the accuracy and depth of maize’s whole genome sequencing and marker development. Relying on the progress of whole-genome sequencing technology and the development of whole-genome association analysis models and methods, and due to the higher level of genetic diversity in the mapping populations, GWAS has been used to analyze the variations of maize seedling and germination traits under low-temperature conditions. This is because GWAS offers increased mapping resolution and accuracy. A total of 338 cross experiments showed that some QTLs for four seedling cold-tolerance traits were detected using GWAS; thirty-two significant loci and thirty-six candidate genes related to stress tolerance were identified, suggesting that heterosis may be related to maize’s cold tolerance [8]. To identify and analyze cold-tolerance traits in 306 dent inbred lines and 206 European flint inbred lines from temperate regions, indirect cold-tolerance traits such as days from sowing to germination, relative chlorophyll content, and quantum yield of photosystem II were studied. Using the GWAS technology, 49,585 SNPs were used for genotyping, and associations between SNPs and cold-tolerance genes were located in both types. A total of 275 significant associated markers were found, and some candidate genes were consistent with current studies and previous reports [9]. A GWAS of 125 maize inbred lines was studied using 10 low-temperature tolerance traits during the seedling stage and the germination stage; finally, 43 SNPs were identified as being associated with low-temperature tolerance [10]. A study conducted a GWAS on 375 inbred lines grown outdoors and in an artificial climate chamber and identified 19 markers associated with low-temperature tolerance. These markers explained 5.7% to 52.5% of the phenotypic variation in the chlorophyll fluorescence parameters during the seedling stage. The candidate genes identified near the markers were related to ethylene signaling, brassinosteroid, and lignin synthesis [11]. A study employed two cold-tolerant inbred lines, 220 and P9-10, and two susceptible lines, Y1518 and PH4CV, to generate three F2:3 populations to detect QTLs associated with the low-temperature germination ability of seeds. Forty-three QTLs were detected, explaining 0.62% to 39.44% of the phenotypic variation. Among them, 17 QTLs explained more than 10% of the phenotypic variation, with 16 inheriting the favorable alleles from the tolerant lines. After constructing a linkage map, three meta-QTLs were identified, including at least two initial QTLs from different populations. mQTL1-1 includes seven initial QTLs for germination and seedling traits, with three explaining more than 30% of the phenotypic variation [12]. GWAS was used to conduct a germination test on 282 inbred lines and 17 loci associated with cold tolerance were identified [13]. GWAS and QTL mapping were performed on two populations; a total of four associated SNPs and twelve QTLs related to cold tolerance were identified, and the results showed that the Zm00001d002729 gene was a potential factor, with its overexpression being able to improve the cold tolerance of crops [14]. Using GWAS, a total of 30 SNPs were identified that were related to low-temperature tolerance during seed germination, and fourteen candidate genes were found to be directly related to these SNPs; in a further study of the linkage between these candidate genes and low-temperature tolerance, ten differentially expressed genes were identified via RNA-seq analysis [15]. Fifteen significant SNPs related to seed germination were identified via GWAS under cold stress in 300 inbred lines; among them, three genomic loci were repeatedly associated with multiple traits. In further candidate gene association analysis, Zm00001d010458, Zm00001d050021, Zm00001d010454, and Zm00001d010459 were identified as cold-tolerance germination-related candidate genes [16]. A total of 187 significant SNPs were identified via GWAS in 836 maize inbred lines, and there were 159 QTLs for emergence and traits related to early growth [17]. Many of the QTL and GWAS analyses have been widely used to express large variations in cold tolerance of maize, and these cited results open up new possibilities for improving cold tolerance and understanding the molecular and genetic mechanism of cold tolerance in maize. In addition, QTL mapping and GWAS can be applied as resources for conducting marker-assisted selection of cold-tolerant varieties, and we can use genomic selection technology to predict cold-tolerant varieties in large maize populations [18].
In this study, a population of 296 excellent inbred lines of maize from China and abroad was used as the study material, and their genotypes were analyzed via genome resequencing. The germination stage was then subjected to low-temperature tolerance identification in the field and laboratory, and indicators such as germination rate and germination index were detected. The TASSEL 5.0 method was used for GWAS to identify associated SNPs, aiming to provide theoretical support and material resources for the gene mining and breeding of low-temperature tolerance in maize.

2. Materials and Methods

2.1. Plant Materials

We selected 296 representative inbred lines of maize (Zea maize L.) from both domestic and international sources, including 232 domestic lines, 36 US lines, and 28 European lines (see Appendix A). The seeds were provided by the maize research institute of Heilongjiang Academy of Agricultural Sciences, and the seed germination rate was above 95%.

2.2. Identification of Low-Temperature Tolerance during Germination in the Field

The experiment was conducted in the period of 2017–2019 at the experimental field of Heilongjiang Academy of Agricultural Sciences. The soil in the field was calcic soil, which is neutral, flat, and uniform. The experiment was conducted in two stages. In the first stage, seeds were sown as soon as the soil temperature at 5–10 cm depth reached and remained above 5 °C, while in the second stage, seeds were sown when the soil temperature at 5–10 cm depth remained stable at or above 10 °C. After sowing, timely irrigation was carried out. A randomized block design with two rows, each 5 m in length, with 20 cm between plants and 65 cm between rows, was employed with single-seed sowing, and three replicates were used. Daily records of soil temperature, maximum and minimum temperatures in the field, and daily average temperature were noted during the experiment. Natural low-temperature treatment was applied to the seeded plots, and the number of seedlings that germinated was recorded accurately every day. After the cessation of seedling germination, the field seedling germination rate was calculated, and the relative seedling germination rate and relative seedling germination index were determined as follows:
germination rate (%) = (number of germination seeds/total number of seeds) × 100
relative germination rate (%) = (germination rate of early sowing treatment/germination rate of appropriate sowing treatment) × 100
germination index = ∑Gt/Dt (Gt represents the number of germination seeds at time t, and Dt represents the corresponding days)
relative germination index (%) = (germination index of early sowing treatment/germination index of appropriate sowing treatment) × 100

2.3. Identification of Low-Temperature Tolerance during Germination in the Laboratory

Fifty plump seeds of each inbred line were selected, surface-sterilized with 0.5% sodium hypochlorite solution for 5 min, and then rinsed three times with sterile water. The sterilized seeds were transferred onto a culture dish lined with filter paper and covered with 3 cm thick vermiculite that was kept moist; these seeds were then allowed to germinate in a low-temperature incubator under dark conditions. Two low-temperature treatment stages were set up; the first included germination at 5 °C for 7 days, followed by 15 °C for 7 days and then 25 °C for another 7 days, whereas the control was germinated at 25 °C for 21 days. The germination of seedlings with germ breaking through the vermiculite was defined as germination, and the number of emerged seedlings was recorded daily. The experiment was carried out in three replicates. The germination rate was calculated, and the relative germination rate and relative germination index were determined to be 2.2.

2.4. Phenotypic Analysis

Data organization and analysis were performed using Microsoft Office Excel 2016 and R version 3.6.2 [19]. Basic statistical quantities were calculated using Microsoft Office Excel. ANOVA was performed using the aov function in the R language with a random blocking model [20]. Correlation analysis was conducted using the cor function in the R language.

2.5. Analysis of Genotype

2.5.1. Analysis of SNPs

The modified CTAB method [21] was used to extract the genomic DNA of 296 maize inbred lines, DNA quality was detected using a NanoDrop 2000 spectrophotometer (Thermo Fisher Scientific Inc., Kanagawa, Japan) and 0.80% agarose gel electrophoresis, and qualified DNA samples were used for SNP typing. In this study, the genotype of the association analysis population was analyzed by means of genotyping by sequencing (GBS). A combination of MseI, NlaIII, and EcoRI endonucleases was used to cleave the genomic DNA of the maize inbred lines, ligate the linker, construct the library, and sequence. After obtaining the original resequencing results, mutation detection was carried out using GATK (Genome Analysis Toolkit). Clean reads were compared to the reference sequence RefGen_v4 B73 using Bowtie2, and the resulting sam file was labeled, sorted, and removed through Picard to obtain a bam file for GATK. The indel around the bam file was re-aligned using GATK, and then SNP/INDEL analysis was performed using GATK’s HaplotyeCaller command. After merging all obtained vcf files, the SNP genotype data of the 296 maize inbred lines were finally obtained.

2.5.2. Analysis of Population Structure

Population structure analysis was performed using the LEA software package v3.3.2 in R [22]. First, the TASSEL 5.0 software [23] was used to remove SNP markers with rare allele frequencies (minor allele frequencies, MAFs) of less than 0.05, and the remaining SNP markers were exported in the ped format. The ped files were converted to geno- and lfmm-formatted genotyping data using the ped2geno and ped2lfmm functions of the LEA software package v3.3.2. Then, the snmf function of the software package was used to calculate the population structure. The number of subpopulations was set from 1 to 10, and each subpopulation was repeated 10 times. The cross-entropy criterion for subpopulation allocation was calculated using the cross-validation method built into the snmf function, and the appropriate number of subpopulations was selected based on this criterion. The Q-matrix was determined based on the maximum genetic similarity of each inbred line.

2.5.3. Analysis of LD

In the analysis of linkage disequilibrium (LD), the TASSEL 5.0 software [23] was first used to remove SNP markers with a minor allele frequency (MAF) of less than 0.05, and the markers were divided into 10 categories according to the 10 chromosomes and arranged based on their physical position from smallest to largest using B73 RefGen_v4 as a reference. Then, a sliding window approach was used to calculate the LD between these SNP markers, with each window consisting of 100 SNPs and sliding by 1 SNP at a time. The LD between markers was measured using r2 [24]. After obtaining the LD between pairs of SNP markers, an LD decay plot was generated as a function of the physical distance between the markers.

2.6. Analysis of Genome-Wide Association

Genome-wide association analysis was mainly performed using TASSEL 5.0 [23]. Based on the analysis of 296 maize inbred lines, high-quality SNP markers were selected for subsequent analysis by removing SNP markers with minor allele frequencies of less than 0.05 using the TASSEL software [23]. The kinship matrix was calculated using the TASSEL software to estimate the relatedness among the 296 maize inbred lines. The first 10 principal components were calculated using the TASSEL software’s PCA function as population structure parameters. The low-temperature tolerance indices of the 296 domestic and foreign elite maize inbred lines were used, together with SNP genotypes, population structures, and relatedness, to perform genome-wide association analysis using a mixed linear model in the TASSEL software. False positives resulting from multiple comparisons in the genome-wide association analysis results were controlled using the Benjamini and Hochberg method for controlling the false discovery rate, and the false discovery rate was set to 0.10 [25].

3. Results

3.1. Phenotypic Analysis of Low-Temperature Tolerance during Germination in the Field

A variance analysis was performed based on the relative germination index and field relative germination rate of 296 maize inbred lines (Table 1). The results show that there were highly significant differences in genotype, environment, and the interaction between genotype and environment for the relative germination index, with all results reaching a significance level of 0.001. For the field germination rate, there were also highly significant differences in genotype and environment, and both reached a significance level of 0.001. Overall, the relative germination index and field germination rate indicate significant differences in low-temperature tolerance among different maize inbred lines.
The phenotypic analysis results of relative germination rate and relative germination index of the 296 inbred lines under natural field conditions are shown in Table 2. The minimum, maximum, and average values of the field average germination rate for the inbred lines are 28.00%, 100.00%, and 66.62%, respectively. The number of inbred lines that fall on the right side of the mean is higher than those on the left side of the mean. In comparison, the relative seedling germination index in 2017 was similar to that in 2018, which was 86.80% and 83.15%, respectively. Overall, the distribution of the field-averaged relative seedling germination rate, the 2017 relative seedling germination index, and the 2018 relative seedling germination index vary widely, and the low-temperature tolerance variation in the inbred lines is relatively high, with a generally normal distribution.
Under suitable sowing and early sowing conditions, the average seedling germination rate of the 296 inbred lines was 86.87% and 62.48%, respectively. Low-temperature stress significantly reduced the germination rate of each inbred line. There were significant differences among the 296 inbred lines in their relative germination index, which could reduce genotypic differences among the inbred lines and better reflect the differences in their cold tolerance.

3.2. Phenotypic Analysis of Low-Temperature Tolerance during Germination in the Laboratory

The results of the variance analysis of indoor relative germination rate are shown in Table 3. The differences between genotypes, environments, and the interaction between genotype and environment were highly significant, reaching a significance level of 0.001. The differences between genotypes and blocks were also highly significant, reaching a significance level of 0.001. Overall, the indoor germination rate indicates significant differences in cold tolerance among different inbred lines.
The phenotypic analysis results of the indoor relative germination rate of the inbred lines are shown in Table 4. From the table, it can be seen that the average relative germination rates in 2018 and 2019 were 79.51% and 84.60%, respectively. Overall, the distribution range of the indoor relative germination rates in 2018 and 2019 was relatively large, indicating a high variation in cold tolerance among different inbred lines.

3.3. Correlation Analysis of Low-Temperature Tolerance during Germination between Field and Indoor

The correlation analysis showed that the indoor relative germination rate in 2018 was significantly correlated with the indoor relative germination rate in 2019 and the relative germination rate in the field, with correlation coefficients of 0.67 and 0.18, respectively, and both reached a significant level of 0.001 (Table 5). The indoor relative germination rate in 2019 was significantly correlated with the relative germination rate in the field, with a correlation coefficient of 0.20, which reached a significant level of 0.001. The relative germination rate in the field was significantly correlated with the field relative germination indices in 2017 and 2018, with correlation coefficients of 0.50 and 0.49, respectively, and both reached a significant level of 0.001. Among the significantly correlated indicators, the correlation coefficient between the indoor relative germination rate in 2018 and that in 2019 was the highest, reaching 0.67, while the correlation coefficient between the indoor relative germination rate in 2018 and the relative germination rate in the field was the lowest, at 0.18.

3.4. Analysis of Genotype

3.4.1. Analysis of SNPs

A total of 24,042 high-quality SNP markers were identified across the entire maize genome (Table 6). The identified SNPs were distributed relatively evenly across the ten chromosomes, with the highest number found on chromosome 1 (3687 SNPs), followed by chromosome 2 (3217 SNPs). The lowest number of SNPs was found on chromosome 10 (2976 SNPs). Of the 24,042 SNPs identified, 98% had a minor allele frequency greater than 0.05, and 36% had a minor allele frequency greater than 0.1. Additionally, SNPs with different minor allele frequencies were distributed relatively evenly across the ten chromosomes.

3.4.2. Analysis of Population Structure

We used the snmf function with a cross-validation technique to calculate and select the appropriate number of subpopulations based on the standard, and we classified the individuals based on their maximum genetic similarity (Q-matrix). When the number of subpopulations is set from 1 to 10, the cross-entropy criterion for assigning subpopulations gradually decreases, but no obvious turning point is observed (Figure 1). When the number of subpopulations is varied, the clustering of the inbred lines is clearly distinguished (Figure 2). When the number of subpopulations is set to five, the 296 excellent inbred lines are divided into five subpopulations. Subpopulations A, B, C, D, and E include 21, 22, 178, 10, and 65 inbred lines, respectively.

3.4.3. Analysis of LD

Using 23,497 SNP markers with an MAF greater than 0.05 for LD analysis, the LD r2 between these SNP markers basically decreases with an increase in genetic distance between the markers, and all values are distributed between 0.000 and 1.000 (Table 7, Figure 3). The average value of LD between the SNP markers on chromosome 9 is the highest, at 0.122, while the average value on chromosome 7 is the lowest, at 0.052 (Table 7).

3.5. Analysis of Genome-Wide Association

3.5.1. Genome-Wide Association Analysis of Low-Temperature Tolerance in the Field

Using the MLM model in the TASSEL software at a significance threshold of p < 1 × 10−4, seven SNP markers associated with low-temperature tolerance were detected based on the average relative germination rate and the relative germination indices in the field in 2017 and 2018. These markers are located on chromosomes 5, 6, 7, and 10, and their phenotypic contributions range from 5.03% to 9.68% (Table 8, Figure 4). Among them, five significantly associated SNP markers were identified based on the relative germination index in 2018, including marker.17002, marker.17003, marker.17009, and marker.17105 located on chromosome 6 and marker.19874 located on chromosome 7, which explained 6.55%, 6.55%, 5.86%, 8.46%, and 7.30% of the phenotypic variation, respectively. No associated SNP loci were identified for the year 2017, which might be due to the lack of effective low-temperature stress between early sowing and sowing at the optimum time in that year.

3.5.2. Genome-Wide Association Analysis of Low-Temperature Tolerance in the Laboratory

Using the relative germination rate in the laboratory during the germination stage as the indicator, a total of 14 SNP loci associated with low-temperature tolerance during germination were detected; these SNP markers are located on chromosomes 1, 3, 4, 5, and 10, explaining 4.84% to 9.68% of the phenotypic variance. Among them, six significantly associated SNP markers were identified using the relative germination rate in 2018 (Table 9, Figure 5), and eight significantly associated SNP markers were identified using the relative germination rate in 2019 (Table 9, Figure 6).

3.5.3. Consistency Analysis of SNP Markers Associated with Low-Temperature Tolerance

Eight significantly associated SNP markers were identified using the indoor relative germination rate in 2019, mainly distributed on chromosomes 1, 3, 4, and 10. When using the indoor relative germination rate in 2018 and the field relative germination index in 2018, five significantly associated SNP markers were identified for each. No significantly associated SNP markers were identified based on the germination index in 2017. Overall, significantly associated SNP markers were distributed on chromosomes 1, 3, 4, 5, 6, 7, and 10, with most markers on chromosome 1 (up to nine), and no significantly associated SNP markers were found on chromosomes 2, 8, and 9. The −Lg(p) values of significantly associated SNP markers ranged from 4.00 to 6.86, with an average of 4.65. The phenotypic variation explained by a single SNP marker ranged from 4.84% to 9.68%, with an average of 6.13%.
Significantly associated SNP markers also showed clustered distribution. Using the indoor relative germination rates in 2018 and 2019, four significantly associated SNP markers (marker.1723, marker.1724, marker.1726, and marker.1729) were identified in the interval of 31,809,859–31,954,983 on chromosome 1, with an average distance of 36.28 Kb between markers. Using the indoor relative germination rates in 2018 and 2019, a significantly associated SNP marker (marker.8339) was identified on chromosome 3 at position 6,292,001, explaining up to 5.61% of the phenotypic variation. Using the indoor relative germination rate in 2018 and the field relative germination rate in 2019, a significantly associated SNP marker (marker.8340) was identified on chromosome 3 at position 6,292,053, explaining up to 6.87% of the phenotypic variation. Using the indoor relative germination rate in 2019 and the field relative germination rate in 2019, a significantly associated SNP marker (marker.14070) was identified on chromosome 5 at position 2,205,723, explaining up to 9.68% of the phenotypic variation.

4. Discussion

In recent years, with the rapid development of sequencing technology and statistical algorithms, GWAS has become one of the most effective methods for identifying genetic variants associated with important agronomic traits in crops [26,27,28,29]. Compared to the traditional linkage analysis, GWAS can use natural populations as materials directly; it can detect more QTLs than traditional QTL mapping by using biparental populations because it uses a larger number of molecular markers and datasets from hundreds of maize inbred lines, which have a rich allelic diversity [30,31]. Moreover, GWAS can analyze multiple phenotypic traits in multiple environments and across multiple time points at the same time. Its high-throughput sequencing and high precision have greatly improved the efficiency of crop breeding [32,33]. Currently, GWAS has greatly advanced genetic research on maize functional genomics [34]; many agronomic traits such as flowering time, leaf angle, leaf size, and disease resistance have been identified in maize. For example, using 368 maize inbred lines and approximately 1 million SNPs, a GWAS analysis successfully detected 74 loci associated with seed oil content and fatty acid composition in maize [35,36]. The US-NAM population was used to detect maize flowering time variants, and a total of 90 flowering time regions were identified in the whole genome via GWAS; among them, one third of regions were associated with the environmental sensitivity of maize flowering time [37]. In another study, 513 inbred lines were used to identify 678 SNPs associated with 17 agronomic traits via GWAS, such as plant height, seed morphology, and flowering time; the results found that 54.3% of these SNPs were associated with at least two or more agronomic traits [38]. A total of 217 inbred lines were genotyped using the GBS technology, and 39 SNPs were identified to be significantly associated with fumonisin resistance in maize kernels based on GWAS analysis [39]. A panel of 143 elite lines were genotyped by using the MaizeSNP50 chip, combined with GWAS and transcriptome analysis; the results showed that 15 common quantitative trait nucleotides were associated with maize white spot, and SYN10137-PZA00131.14 was identified as a key genetic region for improving resistance to MWS; in this region, three candidate genes were identified [40].
Maize can grow in cool-temperate climates but is often exposed to cold temperatures in spring, which can affect seedling growth. Currently, although studies have shown that the growth and development of maize plants are closely related to low temperatures, the genetics of low-temperature tolerance in maize is not well understood. For example, low-temperature stress can increase the expression of related genes, resulting in the accumulation of folate in maize plants [41]. Cold stress can result in a series of physiological responses, such as the expression of osmotic stress-related genes, accumulation of ROS, activities of antioxidant enzymes, and levels of plant hormones and MDA production [42,43,44,45]; thus, plants need to stabilize cell membranes and biologically active proteins in order to survive under low-temperature conditions. However, low-temperature tolerance in maize is a complex trait because the identification and evaluation of low-temperature tolerance traits are complex and have not been standardized. Classic quantitative genetics studies have shown that low-temperature tolerance is controlled by multiple genes and is easily affected by environmental conditions. Quantitative genetic analyses of cold tolerance have shown that genotype, additive effects, growth stage, heterosis, and reciprocal and environmental factors are all involved in the expression of cold tolerance in maize [46]. Six maize lines were used to evaluate the expression of CAT, APX, SOD, and other genes; the results showed that there was heterosis for germination under cold stress, and non-additive genes were more important [47]. The studies cited above show that the genetics mechanisms of low-temperature tolerance in maize are very complicated.
Maize has rich genetic variability, a fast LD decay rate, and abundant information on SNP loci, so maize is an ideal model crop for GWAS analysis [48]. In this study, we used GWAS to identify the genetic loci associated with five traits related to low-temperature tolerance during germination. We identified 30 markers significantly associated with low-temperature tolerance, which were located on chromosomes 1, 2, 3, 4, 5, 6, 7, and 10. Two markers (marker.17002 and marker.17009) significantly associated with the relative germination index in the field in 2018 were located in bin 2.05. This interval has been mapped to several traits under different temperature conditions, including SPAD values [49], antioxidant activity under cold treatment, chlorophyll b, chlorophyll a + b, and Fv/Fm under different temperatures and sowing times [50]. Marker.7569, significantly associated with the relative germination index in the field in 2018, was located in bin 2.06 and was involved in the photosynthetic traits of the third leaf under 15 °C conditions, including CO2 assimilation rate and ΦPSII [1]. Marker.19874, located in bin 2.08, was associated with hundred-grain weight under 14 °C/10 °C (day/night) conditions [51]. Some of these associated markers are consistent with previous studies on low-temperature tolerance, although the traits they are associated with may differ, which may be due to pleiotropy.
Eight significant SNPs related to relative germination rate were detected using the indoor relative germination rate in 2019, and five significant SNPs related to low-temperature tolerance were detected using both the indoor relative germination rate in 2018 and the field relative germination index in 2018. Overall, these significant SNPs related to low-temperature tolerance were distributed on chromosomes 1, 3, 4, 5, 6, 7, and 10, with most SNPs distributed on chromosome 1 (nine SNPs). Previous studies have also shown that SNPs that are associated with seedling-related traits in maize under cold stress are concentrated on chromosomes 1, 2, 3, 5, 6, 8, and 10 [13,52,53]. The above research results further indicate that the cold tolerance of maize is a polygenic quantitative trait controlled by multiple genes. Using polygenic aggregation or multiple molecular markers for the genetic improvement of cold tolerance in maize is an effective strategy. No significant SNPs related to relative germination index were detected in 2017 using GWAS analysis, indicating that low-temperature tolerance in maize may be easily affected by environmental conditions, particularly climate conditions during the growing season. Some scholars have reported that the gene expression related to cold tolerance was affected by the environment of maize [54]. Under controlled conditions, the highest number of significant SNPs related to relative germination rate was detected using the indoor relative germination rate in 2019 (nine SNPs), further demonstrating that controlling low-temperature environmental conditions is important for identifying variations in low-temperature tolerance among different maize inbred lines. Therefore, future studies on identifying genes for low-temperature tolerance and improving maize varieties should focus more on phenotypic evaluations under artificial controlled low-temperature conditions.
Using the indoor relative germination rates in 2018 and 2019, four SNP markers significantly associated with low-temperature tolerance were identified in the range of 31,809,859–31,954,983 on chromosome 1, with an average distance of 36.28 Kb between markers. On chromosome 3, two SNP markers significantly associated with low-temperature tolerance were identified, namely marker.8339 and marker.8340. These markers can be directly applied or developed into easily detectable molecular markers for use in the marker-assisted selection for low-temperature tolerance in maize. In addition, some markers were found to be significantly associated with multiple traits, such as marker.1723 and marker.1724, which are only 43 bp apart and were associated with both indoor germination rates in 2018 and 2019. Marker.8339 and marker.8340 were associated with both indoor relative germination rates in 2018 and 2019, and marker.14070 was associated with both the field average relative germination rate and the 2019 indoor relative germination rate. The reason for the same markers being associated with multiple traits may be due to the strong correlation between traits, and it also indicates that the traits identified in this study are all effective for low temperature-tolerance identification. Another reason may be due to the pleiotropy of genes, where genes not only directly control the expression of a trait through the action of enzymes, but also affect many other traits through the modification of a particular trait. This requires further investigation into the function of candidate genes to avoid any negative effects of gene expression on maize breeding.

5. Conclusions

In the present study, GWAS was performed with 296 maize inbred lines, and a total of 14 SNPs significantly associated with low-temperature tolerance were detected. The SNP consistently linked to low-temperature tolerance in the field and indoors during germination was marker.14070, located on chromosome 5 at position 2,205,723, which explained 4.84–9.68% of the phenotypic variation.

Author Contributions

Conceptualization, T.Y., J.Z. and J.C.; data curation, J.Z., X.L. and Q.C.; formal analysis, T.Y., J.Z., J.C., S.L. (Shujun Li) and C.H.; funding acquisition, J.Z. and J.C.; investigation, T.Y., J.Z., C.H. and X.L.; methodology, J.Z., S.L. (Shujun Li), Y.L. and S.L. (Sinan Li); project administration, J.C.; resources, S.L. (Shujun Li) and Y.L.; software, T.Y., X.L. and C.H.; supervision, J.Z.; visualization, S.L. (Sinan Li) and Q.C.; writing—original draft preparation, T.Y.; writing—review and editing, S.L. (Sinan Li), J.C. and X.M. All authors have read and agreed to the published version of the manuscript.

Funding

This work was funded by the Innovation Project of Heilongjiang Academy of Agricultural Sciences (CX23JQ04, CX23ZD05); the Natural Science Foundation of Heilongjiang Province (LH2022C096); the National Key Research and Development Program of China (2021YFD1201001, 2021YFD1201001-2); and the Key Research and Development Program of Heilongjiang Province (JD22A010).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data used in this study are available within the text.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Table A1. The population distribution characteristics when the number of subgroups is 5.
Table A1. The population distribution characteristics when the number of subgroups is 5.
NO.NameAncestry MatrixSubgroup
ABCDE
1LX460.8600.0500.09A
2UH3020.92000.080A
3LX560.570.14000.29A
4LX570.590.14000.27A
5LX650.520.020.370.050.04A
6LX760.390.090.3200.2A
7LX980.90.1000A
8LX990.670.33000A
938P05f0.7700.0700.16A
1038P05m10000A
11LX11310000A
12YA1M0.900.040.030.03A
13LX1470.8500.060.030.06A
14LX1480.8800.050.030.04A
15LX1510.740.09000.17A
16LX1560.670.33000A
17LX1620.3700.370.160.1A
18LX1710.530.030.060.20.18A
19UH3030.9100.050.040A
20LX1940.450.150.2700.13A
21688F0.940.06000A
22LX9400.480.230.010.28B
23LX1150.010.560.0200.41B
24xy335f00.970.0300B
25LX1320.010.440.1500.4B
26LX1370.040.470.1100.38B
27LY88M00.780.0600.16B
28LX16400.650.0700.28B
29LX16600.68000.32B
30LX16700.680.040.020.26B
31LX16800.540.0100.45B
32LX16900.610.0500.34B
33LX1720.030.55000.42B
34LX18000.880.10.020B
35LX18600.890.090.020B
36LX20500.510.2600.23B
37LX2060.040.480.240.030.21B
38LY99M01000B
39MEIXI0.030.430.170.040.33B
40101M0.010.53000.46B
41420F01000B
42738M01000B
43820F0.050.550.330.070B
44LX10.050.040.490.090.33C
45LX20.030.050.620.160.14C
46LX30.010.060.640.20.09C
47LX400.130.530.170.17C
48LX50.020.020.830.130C
49LX600.020.920.040.02C
50LX80.020.020.590.040.33C
51LX900100C
52LX100.030.030.760.050.13C
53LuYuan920.020.030.720.080.15C
54LX120.2500.450.220.08C
55LX130.020.050.530.050.35C
56LX140.020.050.590.040.3C
57LX150.030.090.690.050.14C
58LX170.020.030.690.040.22C
59LX180.10.030.440.30.13C
60H26100.020.740.050.19C
61LX190.050.050.610.050.24C
62LX2000.010.610.050.33C
63LX210.190.080.370.070.29C
64Q3190.040.060.710.050.14C
65LX230.020.050.730.180.02C
66LX250.10.080.60.070.15C
67LX9500.050.740.210C
68LX27000.8600.14C
69LX290.040.040.660.070.19C
70LX3000.080.560.290.07C
71LX32000.880.120C
72LX3300.070.570.310.05C
73LX3400.020.750.160.07C
74LX350.10.050.450.260.14C
75LX360.270.010.40.210.11C
76LX370.040.060.460.280.16C
77LX38000.9700.03C
78LX39000.820.180C
79LX400.0300.7100.26C
80LX41000.540.150.31C
81LX420.010.050.470.030.44C
82LX43000.870.130C
83LX450.020.030.780.080.09C
84B31700.010.840.080.07C
85LX480.010.020.770.110.09C
86B1440.020.040.470.040.43C
87Si287000.840.160C
887-004000.6600.34C
89SD1900.010.050.780.080.08C
90LX490.020.030.730.010.21C
91LX50000.8900.11C
92LX520.0500.690.050.21C
93LX53000.90.050.05C
94LX550.010.050.620.160.16C
95LX590.060.020.720.050.15C
96LX600.020.040.620.080.24C
97LX610.030.030.730.120.09C
98LX62000.760.240C
99LX6300.080.790.10.03C
100LX6400.010.750.20.04C
101LX660.2300.420.130.22C
102LX670.050.070.530.080.27C
103LX680.020.040.330.290.32C
104LX690.020.010.550.290.13C
105LX710.0200.870.040.07C
106LX7400.210.5800.21C
107LX750.040.180.410.050.32C
108LX820.030.010.480.220.26C
109LX83000.740.220.04C
110LX84000.80.190.01C
111LX850.030.060.490.020.4C
112LX86000.930.070C
113LX87000.950.050C
114LX900.020.040.470.060.41C
115LX92000.740.260C
116LX93000.660.30.04C
117LX100000.530.290.18C
118LX10200.030.860.090.02C
119LM33M0.030.050.610.070.24C
120LX1050.0500.880.050.02C
121LX1060.010.10.810.050.03C
122LX1070.0600.890.050C
123LX1120.070.140.470.090.23C
124LX1140.060.120.40.020.4C
125Huangzao4000.860.140C
126Longxi530.040.040.410.130.38C
1277060.080.050.710.110.05C
128LX1170.010.010.960.020C
129LX1180.030.10.720.060.09C
130LX11900.030.780.140.05C
131K1000.090.470.350.09C
132Chang30.020.030.60.130.22C
133Zhong7490-920.010.040.760.180.01C
134785990.0300.910.060C
135478000.8900.11C
136He34400.010.750.180.06C
137Longkang1100100C
138MO170.040.170.560.030.2C
1394F100100C
140330000.90.10C
1415003000.980.020C
142B730.030.010.9600C
143LX1200.070.030.660.190.05C
144LX1210.020.020.870.070.02C
145LX12200.040.690.270C
146LX1230.0300.830.110.03C
147LX124000.990.010C
14806S021000.9700.03C
14906S03200100C
15006S0340.130.010.780.050.03C
15106S052000.530.010.46C
15206S06000.010.8400.15C
15306S068000.9200.08C
15406S075000.770.140.09C
15506S09500.040.710.10.15C
156R1170.0100.9100.08C
157FLAF0.130.010.480.220.16C
158YAM0.070.080.660.090.1C
159LX12500.050.810.140C
160Jidan27♂0.010.010.490.040.45C
161LX126000.590.010.4C
162LX12700.030.750.040.18C
163LX12800.020.830.130.02C
164LX129000.980.020C
165LX1300.0100.670.010.31C
166LX1310.030.130.610.080.15C
167Zheng580.0100.8700.12C
168Chang7-200.010.80.140.05C
169HY6M0.1800.490.190.14C
170HY6F0.040.030.770.060.1C
171LX1340.2600.640.080.02C
172LX140000.8700.13C
173LX1440.030.090.510.060.31C
174Co117-20.10.040.530.190.14C
175Co2200.090.020.470.30.12C
176Co2280.150.030.40.30.12C
177Co2660.090.040.440.310.12C
178Co2740.070.090.550.150.14C
179Co2850.070.020.510.080.32C
180Co3580.130.020.480.210.16C
181Co3720.070.020.490.270.15C
182Co3730.10.060.470.280.09C
183Co3800.080.030.460.190.24C
184LY88F000.80.10.1C
185LX149000.940.060C
186LX150000.810.080.11C
187LX1520.0200.9300.05C
188LX1530.060.020.710.090.12C
189LX1540.0100.870.120C
190LX15500.020.830.150C
191H127RE0.040.040.640.120.16C
192LX1630.050.070.8200.06C
193DY100.010.880.110C
194DY100.010.030.5300.43C
195DY210.070.120.440.030.34C
196DY530.0100.860.130C
197DY710.020.020.80.160C
198DY9700.010.6300.36C
199DY990.030.030.50.220.22C
200Dy13-170.140.030.460.20.17C
201LX1760.030.020.9500C
202LX17700.060.90.040C
203LX1780.030.060.570.060.28C
204LX18100.070.590.060.28C
205LX1840.010.030.670.010.28C
206LX1850.040.220.480.040.22C
207LX1880.030.260.370.020.32C
208LX1900.090.170.40.020.32C
209LX19300.020.950.030C
210LX19500.010.910.080C
211LX1960.070.010.770.030.12C
212LX19700.010.750.010.23C
213LX19800.010.770.010.21C
214LX1990.0200.7600.22C
215LX2010.040.020.490.020.43C
216LX2020.0500.440.080.43C
217LX2090.130.060.670.080.06C
218LX2100.130.130.540.060.14C
219698F0.040.040.680.010.23C
220820M0.130.140.370.050.31C
221YA2M0.270.090.470.060.11C
222LX2600.030.170.690.11D
223UH00400010D
224FLAM00.040.240.580.14D
225YAF0000.970.03D
226Co3740.210.030.20.450.11D
227YA1F00.0200.940.04D
228LX1590.36000.640D
229LX1600.47000.530D
230LX1610.210.10.060.520.11D
231UH00900.0100.930.06D
232Ji871000.390.10.51E
233LX110.040.030.420.080.43E
234LX1600.270.250.020.46E
235LX22000.120.030.85E
236LX2400.22000.78E
237LX310.030.030.410.040.49E
238LX510.0100.360.050.58E
239LX54000.2100.79E
240LX720.340.030.0900.54E
241LX730.050.260.250.070.37E
242LX770.160.110.260.030.44E
243LX780.140.180.0600.62E
244LX790.150.190.0400.62E
245LX800.060.110.210.010.61E
246LX880.190.110.2800.42E
247LX890.030.020.340.020.59E
248LX91000.090.110.8E
249LX1010.060.250.270.040.38E
250LX104000.130.010.86E
251LX10800.1000.9E
252LX1090.020.170.240.020.55E
253LX1100.130.010.030.040.79E
254LX1110.020.050.330.050.55E
255LX11600.020.0100.97E
256Lv2800001E
257LX1240.380.010.2100.4E
258xy335m00.090.0200.89E
259LX13300.410.160.010.42E
260LX13500.330.090.010.57E
261LX13600.430.120.010.44E
262LX13800.340.230.020.41E
263LX1390.060.230.230.030.45E
264LX1410.020.080.370.060.47E
265LX14200001E
266LX14300001E
267LX1450000.020.98E
268LX1460.010.020.140.010.82E
269Co3710.070.090.280.170.39E
270LX1650.30.040.30.010.35E
271LX1700.020.41000.57E
272DY700001E
273DY290.01000.010.98E
274DY360.050.030.140.020.76E
275DY4900.120.200.68E
276LX17300.070.090.020.82E
277LX17400001E
278LX175000.0400.96E
279LX17900.090.340.050.52E
280LX1820.420.0900.020.47E
281LX1830.060.10.240.20.4E
282LX1870.010.030.360.070.53E
283LX1910.040.230.270.040.42E
284LX1920.180.260.240.010.31E
285LX200000.450.040.51E
286LX20300001E
287LX20400001E
288LX20700001E
289LX20800.070.180.020.73E
290252M0.010.060.210.010.71E
291335MG00.04000.96E
292688M00.04000.96E
293XZD27600.1000.9E
294XZD17000.08000.92E
295XZD1710.020.130.310.020.52E
296YA2F00.33000.67E

References

  1. Fracheboud, Y.; Ribaut, J.M.; Vaigas, M.; Messmer, R.; Stamp, P. Identification of quantitative trait loci for cold tolerance of photosynthesis in maize. J. Exp. Bot. 2002, 376, 1967–1977. [Google Scholar] [CrossRef]
  2. Shi, Y.; Li, G.H.; Tian, Z.Q.; Wang, Z.Y.; Wang, X.B.; Zhu, Y.G.; Chen, Y.H.; Guo, S.L.; Qi, J.S.; Zhang, X.; et al. Genetic dissection of seed vigour traits in maize (Zea mays L.) under low-temperature conditions. J. Genet. 2016, 95, 1017–1022. [Google Scholar] [CrossRef] [PubMed]
  3. Han, Q.; Zhu, Q.; Shen, Y.; Lee, M.; Lübberstedt, T.; Zhao, G. QTL Mapping Low-Temperature Germination Ability in the Maize IBM Syn10 DH Population. Plants 2022, 14, 214. [Google Scholar] [CrossRef] [PubMed]
  4. Goering, R.; Larsen, S.; Tan, J.; Whelan, J.; Makarevitch, I. QTL mapping of seedling tolerance to exposure to low temperature in the maize IBM RIL population. PLoS ONE 2021, 16, e0254437. [Google Scholar] [CrossRef]
  5. Yi, Q.; Malvar, R.A.; Álvarez-Iglesias, L.; Ordás, B.; Revilla, P. Dissecting the genetics of cold tolerance in a multiparental maize population. Theor. Appl. Genet. 2020, 133, 503–516. [Google Scholar] [CrossRef] [PubMed]
  6. He, R.Y.; Zheng, J.J.; Chen, Y.; Pan, Z.Y.; Yang, T.; Zhou, Y.; Li, X.F.; Nan, X.; Li, Y.Z.; Cheng, M.J.; et al. QTL-seq and transcriptomic integrative analyses reveal two positively regulated genes that control the low-temperature germination ability of MTP-maize introgression lines. Theor. Appl. Genet. 2023, 136, 116. [Google Scholar] [CrossRef] [PubMed]
  7. Xu, Q.; Wang, X.; Wang, Y.; Zhang, H.; Zhang, H.; Di, H.; Zhang, L.; Dong, L.; Zeng, X.; Liu, X.; et al. Combined QTL mapping and RNA-Seq pro-filing reveal candidate genes related to low-temperature tolerance in maize. Mol. Breed. 2022, 42, 33. [Google Scholar] [CrossRef] [PubMed]
  8. Yan, J.; Wu, Y.; Li, W.; Qin, X.; Yue, B. Genetic mapping with testcrossing associations and F2:3 populations reveals the importance of heterosis in chilling tolerance at maize seedling stage. Sci. Rep. 2017, 7, 427–428. [Google Scholar] [CrossRef]
  9. Revilla, P.; Rodríguez, V.M.; Ordás, A.; Rincent, R.; Charcosset, A.; Giauffret, C. Association mapping for cold tolerance in two large maize inbred panels. BMC Plant Biol. 2016, 16, 127. [Google Scholar] [CrossRef]
  10. Huang, J.; Zhang, J.; Li, W.; Hu, W.; Duan, L.; Feng, Y.; Qiu, F.Z.; Yue, B. Genome-wide association analysis of ten chilling tolerance indices at the germination and seedling stages in maize. J. Integr. Plant Biol. 2013, 55, 10. [Google Scholar] [CrossRef]
  11. Strigens, A.; Grieder, C.; Haussmann, B.I.G.; Melchinger, A.E. Genetic variation among inbred lines and testcrosses of maize for early growth parameters and their relationship to final dry matter yield. Crop Sci. 2012, 52, 1084. [Google Scholar] [CrossRef]
  12. Li, X.H.; Wang, G.H.; Fu, J.J.; Li, L.; Jia, G.Y.; Ren, L.S.; Lubberstedt, T.; Wang, G.Y.; Wang, J.H.; Gu, R.L. QTL mapping in three connected populations reveals a set of consensus genomic regions for low temperature germination ability in Zea mays L. Front. Plant Sci. 2018, 9, 65–72. [Google Scholar] [CrossRef]
  13. Hu, G.H.; Li, Z.; Lu, Y.C.; Li, C.X.; Gong, S.C.; Yan, S.Q.; Li, G.L.; Wang, M.Q.; Ren, H.L.; Guan, H.T.; et al. Genome-wide association study identified multiple genetic loci on chilling resistance during germination in maize. Sci. Rep. 2017, 7, 10840. [Google Scholar] [CrossRef] [PubMed]
  14. Jin, Y.; Zhang, Z.; Xi, Y.; Yang, Z.; Xiao, Z.; Guan, S.; Qu, J.; Wang, P.; Zhao, R. Identification and Functional Verification of Cold Tolerance Genes in Spring Maize Seedlings Based on a Genome-Wide Association Study and Quantitative Trait Locus Mapping. Front. Plant Sci. 2021, 12, 776972. [Google Scholar] [CrossRef]
  15. Zhang, H.; Zhang, J.; Xu, Q.; Wang, D.; Di, H.; Huang, J.; Yang, X.; Wang, Z.; Zhang, L.; Dong, L.; et al. Identification of candidate tolerance genes to low-temperature during maize germination by GWAS and RNA-seqapproaches. BMC Plant Biol. 2020, 20, 333. [Google Scholar] [CrossRef] [PubMed]
  16. Zhang, Y.; Liu, P.; Wang, C.; Zhang, N.; Zhu, Y.; Zou, C.; Yuan, G.; Yang, C.; Gao, S.; Pan, G.; et al. Genome-wide association study uncovers new genetic loci and candidate genes underlying seed chilling-germination in maize. PeerJ 2021, 9, e11707. [Google Scholar] [CrossRef]
  17. Yi, Q.; Álvarez-Iglesias, L.; Malvar, R.A.; Romay, M.C.; Revilla, P.A. worldwide maize panel revealed new genetic variation for cold tolerance. Theor. Appl. Genet. 2021, 4, 1083–1094. [Google Scholar] [CrossRef]
  18. Zhou, X.; Muhammad, I.; Lan, H.; Xia, C. Recent Advances in the Analysis of Cold Tolerance in Maize. Front. Plant Sci. 2022, 13, 866034. [Google Scholar] [CrossRef]
  19. R Core Team. R: A language and environment for statistical computing. Computing 2011, 1, 12–21. [Google Scholar]
  20. Kabacoff, R.I. R in Action: Data Analysis and Graphics with R; Simon and Schuster: New York, NY, USA, 2015. [Google Scholar]
  21. Murray, M.; Thompson, W.F. Rapid isolation of high molecular weight plant DNA. Nucleic Acids Res. 1980, 8, 4321–4326. [Google Scholar] [CrossRef]
  22. Frichot, E.F.O. Lea: An R package for landscape and ecological association studies. Methods Ecol. Evol. 2015, 6, 925–929. [Google Scholar] [CrossRef]
  23. Bradbury, P.J.; Zhang, Z.; Kroon, D.E.; Casstevens, T.M.; Buckler, E. TASSEL: Software for association mapping of complex traits in diverse samples. Bioinformatics 2007, 23, 2633–2635. [Google Scholar] [CrossRef]
  24. Gaut, B.S.; Long, A.D. The lowdown on linkage disequilibrium. Plant Cell 2003, 15, 1502–1506. [Google Scholar] [CrossRef] [PubMed]
  25. Benjamini, Y.; Hochberg, Y. Controlling the false discovery rate: A practical and powerful approach to multiple testing. J. R. Stat. Soc. 1995, 57, 289–300. [Google Scholar] [CrossRef]
  26. Wang, Q.; Tang, J.; Han, B.; Huang, X. Advances in genome-wide association studies of complex traits in rice. Theor. Appl. Genet. 2020, 3, 1415–1425. [Google Scholar] [CrossRef]
  27. Alseekh, S.; Kostova, D.; Bulut, M.; Fernie, A.R. Genome-wide association studies: Assessing trait characteristics in model and crop plants. Cell Mol. Life Sci. 2021, 8, 5743–5754. [Google Scholar] [CrossRef] [PubMed]
  28. Saini, D.K.; Chopra, Y.; Singh, J.; Sandhu, K.S.; Kumar, A.; Bazzer, S.; Srivastava, P. Comprehensive evaluation of mapping complex traits in wheat using genome-wide association studies. Mol. Breed. 2021, 42, 1. [Google Scholar] [CrossRef] [PubMed]
  29. Xiao, Y.; Liu, H.; Wu, L.; Warburton, M.; Yan, J. Genome-wide Association Studies in Maize: Praise and Stargaze. Mol. Plant 2017, 10, 359–374. [Google Scholar] [CrossRef] [PubMed]
  30. Weng, J.F.; Xie, C.X.; Hao, Z.F.; Wang, J.J.; Liu, C.L.; Li, M.S.; Zhang, D.G.; Bai, L.; Zhang, S.H.; Li, X.H. Genome-wide association study identifies candidate genes that affect plant height in Chinese elite maize (Zea mays L.) inbred lines. PLoS ONE 2011, 6, e29229. [Google Scholar] [CrossRef]
  31. Huang, X.H.; Zhao, Y.; Wei, X.H.; Li, C.Y.; Wang, A.H.; Zhao, Q.; Li, W.J.; Guo, Y.L.; Deng, L.W.; Zhu, C.R.; et al. Genome-wide association study of flowering time and grain yield traits in a worldwide collection of rice germplasm. Nat. Genet. 2012, 44, 32–39. [Google Scholar] [CrossRef]
  32. Xiao, Q.; Bai, X.; Zhang, C.; He, Y. Advanced high-throughput plant phenotyping techniques for genome-wide association studies: A review. J. Adv. Res. 2021, 35, 215–230. [Google Scholar] [CrossRef]
  33. Gupta, P.K.; Kulwal, P.L.; Jaiswal, V. Association mapping in plants in the post-GWAS genomics era. Adv. Genet. 2019, 104, 75–154. [Google Scholar]
  34. Wang, B.; Lin, Z.; Li, X.; Zhao, Y.; Zhao, B.; Wu, G.; Ma, X.; Wang, H.; Xie, Y.; Li, Q.; et al. Genome-wide selection and genetic improvement during modern maize breeding. Nat. Genet. 2020, 52, 565–571. [Google Scholar] [CrossRef]
  35. Li, H.; Peng, Z.; Yang, X.; Wang, W.; Fu, J.; Wang, J.; Han, Y.; Chai, Y.; Guo, T.; Yang, N.; et al. Genome-wide association study dissects the genetic architecture of oil biosynthesis in maize kernels. Nat. Genet. 2013, 45, 43. [Google Scholar] [CrossRef]
  36. Jia, G.; Huang, X.; Zhi, H.; Zhao, Y.; Zhao, Q.; Li, W.J.; Chai, Y.; Yang, L.F.; Liu, K.Y.; Lu, H.Y.; et al. A haplotype map of genomic variations and genome-wide association studies of agronomic traits in foxtail millet (Setaria italica). Nat. Genet. 2013, 45, 957–961. [Google Scholar] [CrossRef]
  37. Li, Y.X.; Li, C.; Bradbury, P.J.; Liu, X.; Lu, F.; Romay, C.M.; Glaubitz, J.C.; Wu, X.; Peng, B.; Shi, Y.; et al. Identification of genetic variants associated with maize flowering time using an extremely large multi-genetic background population. Plant J. 2016, 86, 391–402. [Google Scholar] [CrossRef]
  38. Yang, N.; Lu, Y.; Yang, X.; Huang, J.; Zhou, Y.; Ali, F.; Wen, W.; Liu, J.; Li, J.; Yan, J. Genome wide association studies using a new nonparametric model reveal the genetic architecture of 17 agronomic traits in an enlarged maize association panel. PLoS Genet. 2014, 10, e1004573. [Google Scholar] [CrossRef] [PubMed]
  39. Samayoa, L.F.; Cao, A.; Santiago, R.; Malvar, R.A.; Butrón, A. Genome-wide association analysis for fumonisin content in maize kernels. BMC Plant Biol. 2019, 19, 166. [Google Scholar] [CrossRef] [PubMed]
  40. Wang, D.; He, Y.; Nie, L.; Guo, S.; Tu, L.; Guo, X.; Wang, A.; Liu, P.; Zhu, Y.; Wu, X.; et al. Integrated IBD Analysis, GWAS Analysis and Transcriptome Analysis to Identify the Candidate Genes for White Spot Disease in Maize. Int. J. Mol. Sci. 2023, 24, 10005. [Google Scholar] [CrossRef]
  41. Xiang, N.; Hu, J.G.; Wen, T.X.; Brennan, M.A.; Brennan, C.S.; Guo, X.B. Effects of temperature stress on the accumulation of ascorbic acid and folates in sweet corn (Zea mays L.) seedlings. J. Sci. Food Agric. 2020, 100, 1694–1701. [Google Scholar] [CrossRef] [PubMed]
  42. Ritonga, F.N.; Chen, S. Physiological and Molecular Mechanism Involved in Cold Stress Tolerance in Plants. Plants 2020, 9, 560. [Google Scholar] [CrossRef] [PubMed]
  43. Farhangi-Abriz, S.; Torabian, S. Antioxidant enzyme and osmotic adjustment changes in bean seedlings as affected by biochar under salt stress. Ecotoxicol. Environ. Saf. 2017, 137, 64–70. [Google Scholar] [CrossRef] [PubMed]
  44. Nadarajah, K.K. ROS Homeostasis in Abiotic Stress Tolerance in Plants. Int. J. Mol. Sci. 2020, 21, 5208. [Google Scholar] [CrossRef] [PubMed]
  45. Qin, P.; Zhang, G.; Hu, B.; Wu, J.; Chen, W.; Ren, Z.; Liu, Y.; Xie, J.; Yuan, H.; Tu, B.; et al. Leaf-derived ABA regulates rice seed development via a transporter-mediated and temperature-sensitive mechanism. Sci. Adv. 2021, 7, eabc8873. [Google Scholar] [CrossRef] [PubMed]
  46. Ma, Y. emical Responses’ and QTL/Genes’ Identification. Plants 2022, 11, 2082. [Google Scholar] [CrossRef] [PubMed]
  47. Costa Silva Neta, I.; Vilela de Resende Von Pinho, É.; de Abreu, V.M.; Rezende Vilela, D.; Santos, M.C.; Oliveira Dos Santos, H.; Diniz Cabral Ferreira, R.A.; Garcia Von Pinho, R.; Coelho de Castro Vasconcellos, R. Gene expression and genetic control to cold tolerance during maize seed germination. BMC Plant Biol. 2020, 20, 188. [Google Scholar] [CrossRef]
  48. Liu, M.; Tan, X.; Yang, Y.; Liu, P.; Zhang, X.; Zhang, Y.; Wang, L.; Hu, Y.; Ma, L.; Li, Z.; et al. Analysis of the genetic architecture of maize kernel size traits by combined linkage and association mapping. Plant Biotechnol. J. 2020, 18, 207–221. [Google Scholar] [CrossRef]
  49. Peraza, O.G.; Leipner, J.; Reimer, R.; Nguyen, H.T.; Fracheboud, Y. Temperature at night affects the genetic control of acclimation to cold in maize seedlings. Maydica 2012, 56, 367–377. [Google Scholar]
  50. Jompuk, C. Mapping of quantitative trait loci associated with chilling tolerance in maize (Zea mays L.) seedlings grown under field conditions. J. Exp. Bot. 2005, 56, 1153–1163. [Google Scholar] [CrossRef]
  51. Allam, M.; Revilla, P.; Djemel, A.; Tracy, W.F.; Ordás, B. Identification of QTLs involved in cold tolerance in sweetfield corn. Euphytica 2016, 208, 353–365. [Google Scholar] [CrossRef]
  52. Zhou, Y.; Lu, Q.; Ma, J.; Wang, D.; Li, X.; Di, H.; Zhang, L.; Hu, X.; Dong, L.; Liu, X.; et al. Using a high density bin map to analyze quantitative trait locis of germination ability of maize at low temperatures. Front. Plant Sci. 2022, 22, 978941. [Google Scholar] [CrossRef] [PubMed]
  53. He, R.Y.; Yang, T.; Zheng, J.J.; Pan, Z.Y.; Chen, Y.; Zhou, Y.; Li, X.F.; Li, Y.Z.; Iqbal, M.Z.; Yang, C.Y.; et al. QTL Mapping and a Transcriptome Integrative Analysis Uncover the Candidate Genes That Control the Cold Tolerance of Maize Introgression Lines at the Seedling Stage. Int. J. Mol. Sci. 2023, 30, 2629. [Google Scholar] [CrossRef] [PubMed]
  54. Presterl, T.; Ouzunova, M.; Schmidt, W.; Möller, E.M.; Röber, F.K.; Knaak, C.; Ernst, K.; Westhoff, P.; Geiger, H.H. Quantitative trait loci for early plant vigour of maize grown in chilly environments. Theor. Appl. Genet. 2007, 114, 1059–1070. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Variation in the clustering standard cross-entropy as the number of subgroups K increases.
Figure 1. Variation in the clustering standard cross-entropy as the number of subgroups K increases.
Cimb 45 00602 g001
Figure 2. The population structure of subgroups from K = 4 to K = 6. Note: (a) K = 4; (b) K = 5; and (c) K = 6.
Figure 2. The population structure of subgroups from K = 4 to K = 6. Note: (a) K = 4; (b) K = 5; and (c) K = 6.
Cimb 45 00602 g002
Figure 3. The attenuation in maize chromosome LD with an increase in SNP distance.
Figure 3. The attenuation in maize chromosome LD with an increase in SNP distance.
Cimb 45 00602 g003
Figure 4. Correlation analysis based on the relative germination index in 2018.
Figure 4. Correlation analysis based on the relative germination index in 2018.
Cimb 45 00602 g004
Figure 5. Correlation analysis based on the indoor relative germination index in 2018.
Figure 5. Correlation analysis based on the indoor relative germination index in 2018.
Cimb 45 00602 g005
Figure 6. Correlation analysis based on the indoor relative germination index in 2019.
Figure 6. Correlation analysis based on the indoor relative germination index in 2019.
Cimb 45 00602 g006
Table 1. The ANOVA analysis of traits related to low-temperature tolerance in the field from germination to the seedling stage.
Table 1. The ANOVA analysis of traits related to low-temperature tolerance in the field from germination to the seedling stage.
TraitsSource of VariationDegrees of FreedomSum of SquaresMean SquaresF-Valuep-Value
error56572,604129--
Relative germination indexGenotype292191,85965787.98<2 × 10−10 ***
Environment142434243568.15<2 × 10−10 ***
Block1440.520.47
Genotype × environment281195,40869593.12<2 × 10−10 ***
Error57442877--
Relative germination rateGenotype293113,8563891.010.47
Environment112,85812,85833.33<2 × 10−8 ***
Error285109,946386--
Note: ‘***’ indicates significance at the 0.001 level.
Table 2. The low-temperature tolerance phenotype statistics from germination to the seedling stage.
Table 2. The low-temperature tolerance phenotype statistics from germination to the seedling stage.
TraitsNumberMin.Max.MeanMedianSDKurtosisSkewness
Relative germination rate29328.00100.0066.6267.0014.10−0.39−0.25
Relative germination index in 201729120.69100.0086.8092.6816.652.78−1.66
Relative germination index in 20182844.08100.0083.1589.7020.020.85−1.17
Table 3. The ANOVA analysis of traits related to low-temperature tolerance in the laboratory during the germination stage.
Table 3. The ANOVA analysis of traits related to low-temperature tolerance in the laboratory during the germination stage.
TraitsSource of VariationDegrees of FreedomSum of SquaresMean SquaresF-Valuep-Value
Relative germination rateGenotype287478,817166812.98<2 × 10−10 ***
Environment170,89370,893551.68<2 × 10−10 ***
Block1220.0120.91
Genotype × environment28098,0283502.72<2 × 10−10 ***
Error56572,604129--
Note: ‘***’ indicates significance at the 0.001 level.
Table 4. The low-temperature tolerance phenotype statistics during the germination stage.
Table 4. The low-temperature tolerance phenotype statistics during the germination stage.
TraitsNumberMinMaxMeanMedianSDKurtosisSkewness
Relative germination index in 20182810.00100.0079.5176.0025.250.09−0.10
Relative germination index in 20192880.00100.0084.6092.0020.094.68−2.12
Table 5. The correlation analysis of low-temperature tolerance during the germination stage.
Table 5. The correlation analysis of low-temperature tolerance during the germination stage.
TraitsCorrelation Coefficient
Indoor Relative Germination Rate (2018)Indoor Relative Germination Rate (2019)Field Relative Germination RateField Relative Germination Index (2017)Field Relative Germination Index (2018)
Indoor relative germination rate (2018)10.67 ***0.18 ***0.070.02
Indoor relative germination rate (2019)0.67 ***10.20 ***0.020.11
Field-relative germination rate0.18 ***0.20 ***10.50 ***0.49 ***
Field-relative germination index (2017)0.070.020.50 ***1−0.06
Field-relative germination index (2018)0.020.110.49 ***−0.061
Note: ‘***’ indicates significance at the 0.001 level.
Table 6. The allele frequency characteristics of SNP markers.
Table 6. The allele frequency characteristics of SNP markers.
Chr.SNPMinor Allele Frequency (MAF)
>0.05Percent (%)>0.1Percent (%)>0.2Percent (%)
136873599981383381153
23217315998117436893
3289628339893032913
4288228029794733953
5238423349883835944
6189818689871538905
7157715329758037574
8184817939763534623
92294226499854371115
10135913139752939363
Total24,04223,497988585368403
Table 7. LD and LD attenuation of maize chromosomes.
Table 7. LD and LD attenuation of maize chromosomes.
Chr.LD Decay
(r2 < 0.2)
LD Decay
(r2 < 0.1)
MinMaxMeanMedianSDKurtosisSkewness
11204100.0001.0000.0620.0080.15213.8713.605
232010000.0001.0000.0880.0090.2016.1982.69
31306200.0001.0000.0670.0070.16910.4113.264
41005000.0001.0000.0610.0080.15213.9993.648
52709900.0001.0000.0820.0080.1877.8872.903
6853500.0001.0000.0640.0080.1611.8333.435
7904400.0001.0000.0520.0080.14117.754.069
82208200.0001.0000.0680.0080.16911.3353.378
934019900.0001.0000.1220.0110.242.9032.069
10703400.0001.0000.070.010.16511.3263.351
Note: chr.1 to chr.10 represent chromosome 1 to chromosome 10 in maize, respectively.
Table 8. SNPs of maize with significant correlation with low-temperature tolerance.
Table 8. SNPs of maize with significant correlation with low-temperature tolerance.
TraitsMarkerChr.Physical
Position
Lg (p)Contribution
(%)
Field-relative germination ratemarker.1407052,205,7236.869.68
Relative germination index (2018)marker.17002664,236,7754.566.55
Relative germination index (2018)marker.17003664,236,7814.566.55
Relative germination index (2018)marker.17009664,298,5664.155.86
Relative germination index (2018)marker.17105672,142,7515.698.46
Relative germination index (2018)marker.198747180,326,3885.017.30
Field-relative germination ratemarker.5361063,529,7694.065.03
Note: The physical position of SNP markers was determined in reference to B73 RefGen-v4.
Table 9. SNPs of maize with significant correlations with low-temperature tolerance.
Table 9. SNPs of maize with significant correlations with low-temperature tolerance.
TraitsMarkerChr.Physical
Position
Lg (p)Contribution
(%)
Relative germination rate (2019)marker.1723131,809,859 5.01 6.54
Relative germination rate (2018)marker.1724131,809,902 4.07 5.46
Relative germination rate (2019)marker.1726131,897,277 4.18 5.31
Relative germination rate (2019)marker.1729131,954,983 4.95 6.41
Relative germination rate (2018)marker.833936,292,001 4.16 5.61
Relative germination rate (2019)marker.833936,292,001 5.42 7.12
Relative germination rate (2018)marker.834036,292,053 5.01 6.87
Relative germination rate (2019)marker.834036,292,053 4.87 6.34
Relative germination rate (2019)marker.128164140,575,088 4.545.87
Relative germination rate (2018)marker.1407052,205,723 5.176.26
Relative germination rate (2018)marker.1901022,696,941 4.13 5.56
Relative germination rate (2019)marker.7531090,874,322 4.75 6.17
Relative germination rate (2019)marker.1407052,205,7234.004.84
Relative germination rate (2018)marker.84310100,622,715 4.526.10
Note: The physical position of SNP markers was determined in reference to B73 RefGen-v4.
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

Yu, T.; Zhang, J.; Cao, J.; Li, S.; Cai, Q.; Li, X.; Li, S.; Li, Y.; He, C.; Ma, X. Identification of Multiple Genetic Loci Related to Low-Temperature Tolerance during Germination in Maize (Zea maize L.) through a Genome-Wide Association Study. Curr. Issues Mol. Biol. 2023, 45, 9634-9655. https://doi.org/10.3390/cimb45120602

AMA Style

Yu T, Zhang J, Cao J, Li S, Cai Q, Li X, Li S, Li Y, He C, Ma X. Identification of Multiple Genetic Loci Related to Low-Temperature Tolerance during Germination in Maize (Zea maize L.) through a Genome-Wide Association Study. Current Issues in Molecular Biology. 2023; 45(12):9634-9655. https://doi.org/10.3390/cimb45120602

Chicago/Turabian Style

Yu, Tao, Jianguo Zhang, Jingsheng Cao, Shujun Li, Quan Cai, Xin Li, Sinan Li, Yunlong Li, Changan He, and Xuena Ma. 2023. "Identification of Multiple Genetic Loci Related to Low-Temperature Tolerance during Germination in Maize (Zea maize L.) through a Genome-Wide Association Study" Current Issues in Molecular Biology 45, no. 12: 9634-9655. https://doi.org/10.3390/cimb45120602

APA Style

Yu, T., Zhang, J., Cao, J., Li, S., Cai, Q., Li, X., Li, S., Li, Y., He, C., & Ma, X. (2023). Identification of Multiple Genetic Loci Related to Low-Temperature Tolerance during Germination in Maize (Zea maize L.) through a Genome-Wide Association Study. Current Issues in Molecular Biology, 45(12), 9634-9655. https://doi.org/10.3390/cimb45120602

Article Metrics

Back to TopTop