Skip to main content

Identification of candidate tolerance genes to low-temperature during maize germination by GWAS and RNA-seq approaches

Abstract

Background

Maize (Zea mays L.) is one of the main agricultural crops with the largest yield and acreage in the world. However, maize germplasm is very sensitive to low temperatures, mainly during germination, and low temperatures significantly affect plant growth and crop yield. Therefore, the identification of genes capable of increasing tolerance to low temperature has become necessary.

Results

In this study, fourteen phenotypic traits related to seed germination were used to assess the genetic diversity of maize through genome-wide association study (GWAS). A total of 30 single-nucleotide polymorphisms (SNPs) linked to low-temperature tolerance were detected (−log10(P) > 4), fourteen candidate genes were found to be directly related to the SNPs, further additional 68 genes were identified when the screen was extended to include a linkage disequilibrium (LD) decay distance of r2 ≥ 0.2 from the SNPs. RNA-sequencing (RNA-seq) analysis was then used to confirm the linkage between the candidate gene and low-temperature tolerance. A total of ten differentially expressed genes (DEGs) (|log2 fold change (FC)| ≥ 0.585, P < 0.05) were found within the set distance of LD decay (r2 ≥ 0.2). Among these genes, the expression of six DEGs was verified using qRT-PCR. Zm00001d039219 and Zm00001d034319 were putatively involved in ‘mitogen activated protein kinase (MAPK) signal transduction’ and ‘fatty acid metabolic process’, respectively, based on Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses. Thus, these genes appeared to be related to low-temperature signal transduction and cell membrane fluidity.

Conclusion

Overall, by integrating the results of our GWAS and DEG analysis of low-temperature tolerance during germination in maize, we were able to identify a total of 30 SNPs and 82 related candidate genes, including 10 DEGs, two of which were involved in the response to tolerance to low temperature. Functional analysis will provide valuable information for understanding the genetic mechanism of low-temperature tolerance during germination in maize.

Background

Maize (Zea mays L.) originated from tropical and subtropical regions and has a relatively high temperature threshold for germination [1]. As a consequence, maize is inherently sensitive to low temperatures [2], particularly during germination, and thus is seldom cultivated at higher latitudes or in mountainous regions. When maize is cultivated in cold zones, the plants grow more slowly and have a relatively short growing season, which leads to lower seedling vitality and reduced yields [3, 4]. It is known that the minimum temperature for maize seed germination is approximately 10 °C. When the temperature drops to 6 to 8 °C, irreversible damage to cells and tissues occurs. Under these conditions, seeds normally will not germinate, and the growth of seedlings will stop [1]. Thus, 10 °C is normally chosen as the growth temperature for the identification and screening of proper maize germplasm [5]. In recent years, the occurrence of cold weather has become more uncertain, and low temperature has been occurring more frequently, especially during the germination period and early developmental stage of seedlings. Low-temperature stress not only decrease the both emergence rate of maize seeds and seedling vigor but also increases the chance of pathogenic infection by soil bacteria, which can severely reduce maize yields. Therefore, it is urgent to identify the gene(s) that provide low-temperature-tolerant germination in maize.

At present, significant progress has been made in the identification of genetic loci associated with low-temperature tolerance in maize. Their corresponding genes are normally identified by gene mapping and GWAS. Some indices, such as the emergence rate, seedling emergence index, seedling dry weight, relative average germination time, and percentage of relative viability, have been used for genetic loci identification in maize inbred lines and populations [6, 7]. For instance, a major QTL for low-temperature tolerance of photosynthesis was detected on chromosome 6 using the ETH-DL3 × ETH-DH7 population [8]. Two other QTLs on chromosomes 3.01 and 6.03 were reported, which affect leaf color at low temperatures, and a candidate gene, luteus11, was identified [9]. By the use of a panel of European flint maize inbred lines, 47 lines were found to harbor favorable alleles for six significant QTLs [10]. In another study, 12 QTLs controlling the low-temperature germination rate and primary radicle length were detected on chromosomes 4, 5, 6, 7 and 9 by the use of 243 lines of the intermated B73 × Mo17 (IBM) Syn4 recombinant inbred line (RIL) population [11]. Later, researchers identified 43 QTLs that explained 0.62% ~ 39.44% of the phenotypic variance for low-temperature seed germination in maize by the use of three connected F2:3 populations with inbred lines of two low-temperature-tolerant inbred lines, 220 and P9–10, and two susceptible lines, Y1518 and PH4CV [5]. Several SNPs associated with low-temperature-tolerance traits at the seedling and seed germination stages were detected using GWAS. For example, 43 SNPs associated with 10 low-temperature-tolerance traits in maize seedlings or seed germination were found, although no overlapping SNPs were identified at either developmental stage [12]. A total of 19 markers related to low-temperature tolerance were identified through GWAS of 375 inbred lines in the field and in growth chambers; the loci at these markers explained 5.7% ~ 52.5% of the phenotypic genetic variation at the seedling growth stage and chlorophyll fluorescence parameters [13]. A further GWAS revealed 18 candidate genes from 17 genetic loci associated with low-temperature-tolerant germination, where 10 candidate genes were supported by previous QTL studies [14].

Numerous low-temperature-tolerance genes have been identified in many crop species, and their underlying molecular mechanisms have been studied. In rice, COLD1 interacts with the G-protein alpha-subunit, which in turn promotes GTPase activity, resulting in the activation of Ca2+ channels for low-temperature-induced responses [15]. During the reproductive growth period of rice, the protein kinase gene CTB4a is involved in maintaining high pollen fertility under low-temperature conditions, which increases the seed setting rate and yield. In addition, CTB4a interacts with AtpB to regulate the ATP content under low temperatures to enhance low-temperature tolerance [16]. The MADS-box family transcription factor OsMADS57 interacts with OsTB1 to regulate low-temperature tolerance in rice and helps balance growth and defense responses, which are dependent on OsWRKY94, their common target gene [17]. In addition to transcriptional control, phosphorylation of the Basic Transcription Factor 3 Like (BTF3L) protein by the protein kinase OST1 promotes the interaction between BTF3L and CBFs, which increases the stability of CBF proteins and enhances plant tolerance to low temperatures [18]. Similarly, numerous studies investigating the mechanisms of low-temperature tolerance in maize have been reported, but they have had a limited impact on maize breeding. Several maize genes, including ZmCDPK1, ZmSEC14p and ZmMPK5, with roles in low-temperature tolerance have been identified [19,20,21,22]. Some genes that control kernel weight or kernel number at low temperatures have also been cloned, including genes that were shown to induce the expression of bZIP-type and ERF/AP2-type transcription factors [23, 24].

There have been several studies on the low-temperature tolerance of maize during germination, but detailed evaluation indices are lacking. Thus, we established a standard evaluation system in the present study for genetic mapping, with particular consideration for the actual sowing conditions that occur during early spring in Northeast China. Low-temperature evaluation methods for maize often include both field-based and indoor identification methods. Compared with the indoor methods, the field-based methods are more objective and can be used to evaluate the low-temperature tolerance of maize during each growth period and to analyze genotype by environment interactions. However, such field studies are easily affected by changes in climatic conditions at different locations and different years, and environmental conditions are difficult to control. There are many confounding factors, and repeatability is generally poor. Indoor evaluation methods have several advantages. They are not limited by seasons, they have fewer confounding factors, the environmental conditions are easily controlled, and different temperature gradients can be applied during the experiment to accurately assess the tolerance of different maize varieties to low temperature; however, one disadvantage to this method is that it is generally unable to assess genotypes. Because of environmental interaction effects, two methods of evaluation, both indoor and outdoor, were used, and an improved technical system was established for the present study. Many performance indicators have been applied to evaluate the low-temperature tolerance of maize germination, such as the germination potential, germination rate, germination index, radicle length and germ length [5, 14]. For evaluation of seed germination ability at low temperatures, the ratio (relative value) of various traits measured at low temperature and normal temperature is usually used as an index to eliminate the differences in the genetic background of different materials [12, 14].

With the above established evaluation system, a maize population panel of 222 diverse inbred lines was used for the analysis of traits correlated with low-temperature tolerance via GWAS. Candidate genes were predicted based on RNA-seq data to achieve the following objectives: (1) identify potential SNPs responsible for low-temperature tolerance during seed germination, (2) identify maize inbred lines with extremely low-temperature tolerance, and (3) predict and identify the involved candidate genes for future studies and agricultural applications.

Results

Low-temperature germination ability of the maize lines

The germination rate (Fig. 1a) and seedling performance (Fig. 1b, c, d, e) of the 222 maize lines were evaluated. The following 14 traits were measured: the relative germination rate (RGR), relative germ length (RGL), relative radicle length (RRL), relative radicle surface area (RRSA), relative radicle volume (RRV), relative germination index (RGI), relative vitality index (RVI), relative simple vitality index (RSVI), XiangYang relative germination rate (XYRGR), XiangYang relative germ length (XYRGL), XiangYang relative simple vitality index (XYRSVI), KeShan relative germination rate (KSRGR), KeShan relative germ length (KSRGL), and KeShan relative simple vitality index (KSRSVI). Descriptive statistics of the relative values of the germination traits under low-temperature conditions and normal conditions, including the numbers of observations (n), means, medians, standard deviations (SDs), range, kurtosis and skewness, were calculated. For most of the traits, considerable phenotypic variation was detected among the lines, with medians that ranged from 0.089 (for the RVI) to 0.788 (for the RGR). The RGR varied from 0.027 to 1.000, with a mean of 0.737, and the RVI varied from 0.001 to 0.485, with a mean of 0.101. In particular, the RGL, RRL, RRSA, and RRV of several inbred lines were larger than 1.000, suggesting that the germ and radicle growth of these inbred lines was stimulated by low-temperature treatment. In general, the segregation of most traits in the panel fit a normal distribution except for that of the RGR, RRSA, RRV, RGI, and RVI (Table 1). ANOVA revealed highly significant differences (P < 0.001) among the genotypes for all the relative traits. The broad-sense heritability (H2) of the 14 relative traits ranged from 87.24% (for the XYRGL) to 99.21% (for the RVI) in the association panel (Additional file 1: Table S1).

Fig. 1
figure1

Seed germination of the four maize inbred lines under low-temperature conditions. a Germination rate comparison among four maize inbred lines in the low-temperature treatment (10 °C) and the control treatment (25 °C); b Comparison of the germ length and radicle length of maize inbred line Mo17; c inbred line Ji 853; d inbred line K10; and e inbred line Shen5003

Table 1 Descriptive statistics for germination traits under low-temperature and normal conditions

In total, 14 traits, which were monomodally distributed, were measured (Fig. 2). Correlation analysis was carried out on the relative values of each index under low-temperature stress. The traits during the germination period were closely related to the regulation of maize growth (P < 0.05), except for the RGR and RGL. For the indoor traits, the RSVI and RGL (0.83), RVI and RGI (0.81), and RRSA and RRL (0.78) showed strong correlations (P < 0.001). There were significant positive correlations between the other traits, indicating that the eight traits selected could be used as important indicators for low-temperature tolerance in maize germination. Except for the XYRGL and KSRGL, most traits in the field reached a significant level of correlation (P < 0.05). More significant positive correlations were found between the other traits in the field, such as the XYRGR and XYRSVI (0.94) and the KSRGR and KSRSVI (0.89) (P < 0.001). However, correlations were weak between most traits under indoor and field conditions, such as the RSVI and XYRGL (0.11) (P < 0.05) (Fig. 2). Pairwise correlation coefficients among the eight indoor traits were significant and positive for the low-temperature conditions (r = 0.264–0.948), and pairwise correlation coefficients among the six field traits were also significant and positive for the low-temperature conditions (r = 0.310–0.970). Under normal conditions, most indoor traits and field traits showed a weak correlation, which may be due to the different effects of these environments (Additional file 2: Table S2).

Fig. 2
figure2

Distributions of and correlations between 14 relative phenotypic traits. The frequency distribution of each trait is shown on a central diagonal in the form of a histogram. Scatter plots of between every pair of traits are shown in the areas below the diagonal, and numerical correlation coefficients between every pair of traits are shown in the areas above in the diagonal. *, ** and *** indicate significance at P < 0.05, P < 0.01 and P < 0.001, respectively

Population structure and LD decay

The population structure of the 222 inbred lines was calculated using STRUCTURE version 2.3 [25]. ΔK attains the peak value when the K = 4, it is indicated that the 222 inbred lines could be divided into the four subpopulations (Additional file 3: Fig. S1). A total of 40,697 high-quality SNPs were used in LD analysis for all the chromosomes. LD varied along each chromosome. When r2 = 0.1, the mean LD decay was 395–1125 kb across all chromosomes. The distribution of r2 values along with the physical distance for each chromosome is shown in Additional file 4: Table S3 and Additional file 5: Fig. S2.

Associated SNPs for GWAS

GWAS was conducted on the 14 relative traits (RGR, RGL, RRL, RRSA, RRV, RGI, RVI, RSVI, XYRGR, XYRGL, XYRSVI, KSRGR, KSRGL, KSRSVI) for the 222 maize inbred lines. The GWAS of the SNP markers and traits was performed using the mixed-linear model (MLM) combined with population structure and kinship using TASSEL 5.0 software. A total of 30 SNPs showed a highly significant association with 13 traits except the KSRGL, which may be due to the different effects of these environments (the P values ranged from 1.70E-07 to 9.95E-05; Fig. 3).

Fig. 3
figure3

Manhattan plots of GWAS results showing the significant SNPs associated with 14 relative phenotypic traits. The traits included the following: the relative germination rate (RGR), relative germ length (RGL), relative radicle length (RRL), relative radicle surface area (RRSA), relative radicle volume (RRV), relative germination index (RGI), relative vitality index (RVI), relative simple vitality index (RSVI), XiangYang relative germination rate (XYRGR), XiangYang relative germ length (XYRGL), XiangYang relative simple vitality index (XYRSVI), KeShan relative germination rate (KSRGR), KeShan relative germ length (KSRGL), KeShan relative simple vitality index (KSRSVI)

Among the 30 associated SNPs, three SNPs (SYN30173, PZE-105034900, PZE-109053558) were associated with two traits (RGL and RSVI, RRV and RVI, XYRGR and XYRSVI). These SNPs were distributed on all 10 maize chromosomes, with the highest number on chromosome 2, which contained five SNPs. For each trait, such as the RGR, RRL, XYRGR, XYRGL, XYRSVI, KSRGR and KSRSVI, one SNP was found. Two SNPs were detected for the RGL, and the most significant association was for SYN30173, with a P value of 3.75E-05. Ten SNPs were associated with the RRV and are located on chromosomes 1 (PZE-101192647), 2 (SYN21841), 4 (PZE-104024779), 5 (PZE-105034900), and 7 (PZE-107036384). Three SNPs were located on chromosome 6 (PZE-106003222, PZE-106129965 and PZE-106130106), and two SNPs (PZE-108064544 and SYN26538) were found on chromosome 8. Four SNPs were found for the RRSA (SYN29778, PZE-108068725, PZE-103072583, PZE-110029252), RGI (PZE-101120376, PZE-102099570, PZE-106097864, SYN5516) and RVI (SYN25470, PZE-102100684, PZE-105034900, PZE-110057591). Two SNPs were associated with the RSVI—one on chromosome 5 (SYN30173) and the other on chromosome 7 (SYN4961). Thirteen relative traits contained 30 SNPs explained 12.0–28.6% of the total phenotypic variation (Table 2).

Table 2 SNPs for low-temperature germination-related traits detected in 222 maize inbred lines

Identification of candidate genes in maize

A total of 14 candidate genes were found to be directly associated with SNPs using the B73 RefGen_v4 Maize Gene Database (http://www.maizegdb.org/) (Table 3). Two of these genes (Zm00001d018864 and Zm00001d013415) which are directly associated with SNPs, were correlated with the RSVI, while Zm00001d013415 was also correlated with the RGL. Three candidate genes (Zm00001d041438, Zm00001d010670 andZm00001d002676) were associated with the RRSA, and three genes (Zm00001d049442, Zm00001d010497 and Zm00001d013794) were associated with the RRV. Two genes (Zm00001d029193 and Zm00001d013794) were associated with the RVI, while Zm00001d013794 was correlated with both the RRV and RVI. For the RGR (Zm00001d002243), RGI (Zm00001d023538), XYRGL (Zm00001d025380), KSRGR (Zm00001d034313), and KSRSVI (Zm00001d008224), only one candidate gene was associated.

Table 3 Functions of candidate genes that are directly associated with germination traits

The mean LD decay of the 10 chromosomes across the 222 maize lines was 60 kb, and the LD decay of chromosomes 1 to 10 was 55, 60, 110, 100, 60, 45, 80, 70, 60, and 110 kb (r2 = 0.2), respectively. A 55 kb window was selected to fall within the estimated window of LD decay on chromosome 1; 60 kb, chromosome 2; 110 kb, chromosome 3; and so on. Furthermore, an additional 68 candidate genes were identified within the LD decay regions of 26 associated SNPs (Additional file 6: Table S4).

Differential gene expression and determination of candidate genes

To help identify candidate genes for the identified SNPs, the low-temperature-resistant maize line (Zao 8–3, referred to as 55) and the low-temperature-sensitive line (Ji 853, referred to as 102) were selected from among the 222 inbred lines, and RNA-seq analysis was performed to evaluate genome-wide gene expression levels. A total of 4982 DEGs were identified as downregulated in the CT_102vsCT_55 comparison group, 5550 DEGs were identified at upregulated in the CT_102vsCT_55 comparison group, 5477 DEGs were identified as downregulated in the LT_102vsLT_55 comparison group, and 5661 DEGs were identified as upregulated in the LT_102vsLT_55 comparison group. The results showed that 1022 DEGs were downregulated only in the LT_102vsLT_55 comparison group and that 985 DEGs were upregulated only in the LT_102vsLT_55 comparison group. Four DEGs (|log2FC| ≥ 0.585, P < 0.05) were downregulated in the LT_102vsLT_55 comparison group but upregulated in the CT_102vsCT_55 comparison group. Three DEGs (|log2FC | ≥ 0.585, P < 0.05) were upregulated in the LT_102vsLT_55 comparison groups but were downregulated in the CT_102vsCT_55 comparison group (Fig. 4).

Fig. 4
figure4

Venn diagram of DEGs distributed in maize inbred lines 102 and 55. The treatments included the normal control treatment (CT) and the low-temperature treatment (LT). The low-temperature-tolerant line (Zao 8–3) and low-temperature-sensitive line (Ji 853) are labeled “55” and “102”, respectively. The five treatment-line comparisons groups included the following: low-temperature treatment line 102 vs low-temperature treatment 55, upregulated expression (LT_102vsLT_55.up); low-temperature treatment line 102 vs low-temperature treatment 55, downregulated expression (LT_102vsLT_55.down); normal control treatment line 102 vs normal control treatment 55, upregulated expression (CT_102vsCT_55.up); normal control treatment line 102 vs normal control treatment 55, downregulated expression (CT_102vsCT_55.down); and candidate genes

A total of 82 candidate genes were associated with SNPs, and ten important DEGs were identified, including six candidate genes (Zm00001d034319, Zm00001d025379, Zm00001d021653, Zm00001d039219, Zm00001d002241, and Zm00001d002677) that were downregulated only in the LT_102vsLT_55 comparison group and four candidate genes (Zm00001d002676, Zm00001d038373, Zm00001d029193, and Zm00001d012512) that were upregulated only in the LT_102vsLT_55 comparison group. The putative identities of the proteins encoded by these genes are listed in Table 4. Zm00001d034319 encodes an inositol phosphoryl ceramide-B C-26 hydroxylase, Zm00001d025379 encodes a photoperiod responsive protein, Zm00001d021653 encodes glucose-6-phosphate/phosphate translocator 2, Zm00001d039219 encodes a pleckstrin homology (PH) domain-containing protein, Zm00001d002241 encodes an AT hook motif family protein, Zm00001d002677 encodes coiled-coil domain-containing protein 124, Zm00001d002676 encodes a pseudouridine synthase family protein, Zm00001d038373 encodes a fasciclin-like arabinogalactan family protein (SOS5), and Zm00001d029193 encodes a monocopper oxidase.

Table 4 Candidate DEGs that are inconsistent between the CT_102vsCT_55 and LT_102vsLT_55 comparison groups

Functional prediction of candidate genes

A total of 26 GO terms were found to be associated with the eight candidate genes, and three KEGG pathways were found to be associated with one candidate gene (Tables 5, 6; Additional file 7: Table S5). These GO terms are associated with multiple functions. The first type of function is related to biological processes, including carbohydrate transmembrane transport (Zm00001d021653), fatty acid metabolic processes (Zm00001d034319), protein ubiquitination (Zm00001d025379), fucose metabolic process, protein glycosylation, and cell adhesion (Zm00001d038373). The second type is related to molecular functions, which involve copper ion binding and oxidoreductase activity (Zm00001d029193), fatty acid alpha-hydroxylase activity (Zm00001d034319), ubiquitin-protein transferase activity (Zm00001d025379), transferase activity, transferring glycosyl groups (Zm00001d038373), and nucleoside diphosphate kinase activity (Zm00001d039219). The third type of function is related to cellular components, which involve the endoplasmic reticulum membrane (Zm00001d034319), chloroplasts (Zm00001d039219), plant-type cell wall, plasmodesmata and mitochondria (Zm00001d029193), the endoplasmic reticulum (Zm00001d034319), the Golgi apparatus (Zm00001d038373), the cytoplasm (Zm00001d025379 and Zm00001d002677), the nucleus (Zm00001d039219 and Zm00001d002676), and integral components of membrane (Zm00001d038373 and Zm00001d029193) (Table 5; Fig. 5a). Two KEGG pathways were identified. The first type of function is related to MAPK signal transduction, and the second type of function is related to nucleotide metabolism, including both purine metabolism and pyrimidine metabolism (Table 6; Fig. 5b).

Table 5 GO analysis of candidate genes related to the RNA-seq results
Table 6 Distribution of genes and pathways related to the RNA-seq results
Fig. 5
figure5

GO annotations and KEGG pathways of DEGs. a Bubble chart of GO classifications of the DEGs. The green circle represents the GO annotations common to the Zm00001d039219 and Zm00001d002676 genes. The purple circle represents the GO annotations common to the Zm00001d029193 and Zm00001d038373 genes. The blue square represents the GO annotations common to the Zm00001d034319 gene. A total of 8 DEGs were annotated; b Histogram of the KEGG pathways of the annotated DEGs

Validation of qRT-PCR for differentially expressed genes

Six candidate genes that showed significant differential expression levels according to the RNA-seq were chosen for further qRT-PCR analysis: Zm00001d039219, Zm00001d029193, Zm00001d002676, Zm00001d021653, Zm00001d034319 and Zm00001d025379. A comparison of the qRT-PCR results and RNA-seq data showed consistent expression trends for all six genes. Four genes (Zm00001d039219, Zm00001d021653, Zm00001d034319, and Zm00001d025379) showed significantly higher expression levels in 55 than in 102 under 2 h and 4 h of treatment with low temperatures, and two genes (Zm00001d029193 and Zm00001d002676) showed significantly lower expression levels in 55 than in 102 (Fig. 6).

Fig. 6
figure6

qRT-PCR validation of the GWAS and RNA-seq results. Expression of six candidate genes in the low-temperature-resistant line Zao 8–3 (referred to as 55) and the low-temperature-sensitive line Ji 853 (referred to as 102). Expression analysis was conducted on embryos that were collected at 2 h and 4 h under low (10 °C) and normal (25 °C) temperature conditions, respectively. * and ** indicate significance at P < 0.05 and P < 0.01, respectively

Two candidate genes are of particular interest: Zm00001d039219 and Zm00001d034319. They showed significantly upregulated expression in 55 and showed significant downregulated expression in 102 under all conditions (both low temperatures for 2 h and 4 h) and are involved in functions related to low-temperature resistance, such as MAPK signaling and fatty acid metabolic processes.

Discussion

Relationship between phenotypic characteristics and low-temperature tolerance of maize

The entire growth process of maize is affected by low-temperature stress beginning at the seed germination stage. The most important indicator of low-temperature tolerance during germination is root emergence [1], which is a critical factor for plant development and yield. Low temperature stress decreases root activity, shortens the root length, and results in fewer lateral roots in affected plants [26]. In accordance with previous studies, root and shoot growth were evaluated in controlled growth chambers under a range of temperature regimes [12]. Low-temperature stress also affects maize seedlings after they germinate. It was found that the relative water content, leaf area and leaf dry weight, plant height, root length, stem length and dry weight, and whole-plant fresh weight can be affected [12, 13, 27]. In our study, low-temperature conditions were defined as 10 °C and, for normal conditions, 25 °C. The germination traits assessed were radicle length, radicle surface area, and radicle volume. We observed a large phenotypic variation in radicle length among the 222 maize inbred lines under low-temperature conditions and a strong correlation between radicle length and germination rate. However, their genetic loci were different, indicating that different mechanisms are involved. Our study focused on germination, which was defined as radicle emergence from seeds, under low-temperature conditions. The 30 associated SNPs and two candidate genes identified in this study provide valuable resources for future studies to improve the understanding of the genetic underpinnings of low-temperature tolerance in maize and to improve maize varieties through breeding.

Integrating GWAS data and RNA-seq data for candidate gene prediction

GWAS has been widely used to identify potential candidate genes for important abiotic stress traits in maize [28], but there are still the problems such as the false positives and so on. RNA-seq has become the preferred technique for detecting genome-wide gene expression patterns [29]. However, it is difficult to identify potential key candidate genes as large amounts of DEG are usually obtained through RNA-seq. Recently, the method of integrating GWAS and RNA-seq has been widely adopted to predict candidate genes. For example, a method combining GWAS results with linked DEGs and co-expression network analysis was used to identify seven candidate genes that respond to drought stress in maize [28]. We identified 10 candidate genes related to seed germination at low temperatures through this method. Among them, the two candidate genes, Zm00001d039219 and Zm00001d034319 could provide valuable information for understanding the genetic mechanism of low-temperature tolerance during maize germination.

Consistent SNPs in previous reports

To evaluate the reliability of SNPs detected in this work, we compared the 30 SNPs identified in the present study to those identified in several related publications. Three SNPs overlapped with the physical positions of published QTLs (Additional file 8: Table S6). QTL-8 for ФPSII at the seedling stage [30] was consistent with PUT-163a-149,007,696-748. Two SNPs (PZE-102099570 and PZE-102100684) on chromosome 2 were located in QTL regions associated with straw dry weight [31]; leaf greenness (SPAD) and trapping efficiency of PSII (F′v/F′m)at 15 °C [8]. In addition, a candidate gene (Zm00001d010671) that was strongly correlated with a nearby SNP (PZE-108068725) was supported by a previously identified candidate gene [5]. Thus, our analysis successfully revealed the SNPs known to be associated with low-temperature tolerance, indicating that the identified SNPs from the present study are highly reliable for use in gene cloning and maize breeding.

MAPK signaling pathways in response to low-temperature stress

MAPKs are serine–threonine kinases that mediate intracellular signaling and play vital roles in regulating plant growth, development, and stress responses. At present, many protein kinase genes, including those encoding MAPKs, have been shown to mediate the transduction of abiotic stress response signals [32]. In plants, the accumulation of permeants and antioxidants can be induced by low temperature, drought and salt stress, which is mediated by MAPK pathways in yeast and animals [33]. These MAPK pathways are activated by different stimuli via receptors such as protein tyrosine kinases, G-protein-coupled receptors and two-component histidine kinases. Arabidopsis has approximately 60 MAPKKKs, 10 MAPKKs, and 20 MAPKs. These kinases can be activated by low temperature and other abiotic stresses and are thought to be important components of abiotic stress signaling [34]. In Medicago sativa, low-temperature treatment results in the activation of a MAPK within ten min [35]. Similarly, in tobacco cells, a MAPK and another protein kinase were shown to be activated by osmotic stress in response to Ca2+ or in an ABA-independent manner [36]. The maize gene ZmMAPK5 showed increased expression in response to specific low-temperature treatments [22]. To date, different mechanisms underlying the low-temperature response have been proposed, and relevant coordinated regulatory networks have been analyzed in rice and Arabidopsis. In particular, the signal transduction pathways between MPK activation and ICE1 stability at low temperatures have been elucidated. This marks an important breakthrough in the field of plant regulation in response to low temperature [37,38,39] and reveals the important role of MAPK cascade signals. A candidate gene (Zm00001d039219) reported in this study is related to the MAPK signaling pathway, which may be related to low-temperature tolerance in maize.

Functional analysis of candidate genes

Two candidate genes (Zm00001d039219 and Zm00001d034319) have putative functions related to low-temperature resistance, such as MAPK signaling [22] and fatty acid hydroxylase activity [40], in maize and other species. Zm00001d039219 encodes a pleckstrin homology (PH) domain-containing protein, which is homologous to Arabidopsis AT4G23895 and rice LOC_Os05g51710.1. Pleckstrin homology (PH) domains are typically involved in targeting proteins to the appropriate cellular location and in protein-protein interactions. Despite minimal sequence conservation, they share a common electrostatically polarized fold. Some (< 10%) PH domains bind phosphoinositide phosphates (PIPs) with high specificity and affinity. They are found in a wide range of cellular signaling proteins, including serine/threonine kinases, adaptors, cytoskeletal-associated molecules, lipid-associated enzymes, tyrosine kinases, regulators of G-proteins, and endocytotic GTPases [41]. The putative protein encoded by Zm00001d034319 is an inositol phosphorylceramide-B C-26 hydroxylase that belongs to the fatty acid hydroxylase superfamily and shares sequence homology with the fatty acid hydroxylase in Arabidopsis (AT2G34770) and in rice (LOC_Os03g56820.1). The fatty acid hydroxylase superfamily includes both fatty acid and carotene hydroxylases and sterol desaturases. Beta-carotene hydroxylase hydroxylates beta-carotene in zeaxanthin synthesis and may be involved in other pathways. Other family members include C-5 sterol desaturases and C-4 sterol methyl oxidases. The family members containing two copies of the HXHH motif are involved in cholesterol biosynthesis and biosynthesis of plant cuticular wax; these members are typically integral membrane proteins [41]. Maize cell membrane fluidity significantly decreases under low-temperature stress, and the normal physiological function of membrane-bound proteins is lost. Alteration of the composition of membrane lipid fatty acids through genetic manipulation was shown to improve the low-temperature tolerance of plants [40].

In conclusion, two candidate genes (Zm00001d039219 and Zm00001d034319) were found to have significantly different gene expression levels under low-temperature treatment in resistant and sensitive maize lines. Homologs of these genes in Arabidopsis and rice have functions related to low-temperature resistance to stress, so these genes are attractive candidate genes for involvement in low-temperature tolerance in maize.

Conclusion

In the present study, we contributed to the understanding of the genetic control of low-temperature tolerance in maize at the germination stage by the use of a panel of 222 maize inbred lines. Through GWAS, RNA-seq analysis, and validation via qRT-PCR, the two candidate genes (Zm00001d039219 and Zm00001d034319) were found to have significantly different gene expression levels under low-temperature treatment in resistant and sensitive maize lines.

Methods

Plant materials

A panel consisting of 222 maize inbred lines was used for the association study (This seeds were originally acquired from researcher Xinhai Li and associate researcher Jianfeng Weng, Institute of Crop Science, Chinese Academy of Agricultural Sciences) [25, 42, 43]. These inbred lines were shown to have a wide range of variation in yield components and biotic stress tolerance. They are generally grown in the Yellow and Huai River valley regions, which are in northeastern and southwestern China [42]. Seeds were produced in the spring of 2016 by manual self-pollination of each plant in at Harbin city, Heilongjiang Province, under temperate climatic conditions. Plants were well managed and were kept free of any disease, insect or weed issues during the whole growing season. Harvested seeds were fully dried and then stored at 4 °C. The 222 maize inbred lines used for the GWAS had germination rates greater than 90% at 25 °C in previous studies in our laboratory (Additional file 9: Table S7).

Two maize lines that were resistant (Zao 8–3, referred to as 55) and sensitive (Ji 853, referred to as 102) to low temperature were selected from among the 222 inbred lines, and RNA-seq was performed to evaluate whole-genome gene expression levels.

Seed germination and field experiments

The surface of the seeds was disinfected with 1% sodium hypochlorite (NaOCl) for 5 min, and the seeds were rinsed three times with sterile distilled water for germination experiments. After 6 h of imbibition at normal temperature (25 °C), in accordance with the International Seed Testing Association (https://www.seedtest.org/en/home.html, ISTA) protocol, germination experiments were performed on germination paper in a dark chamber at 10 °C for low-temperature conditions (treatment) and at 25 °C for normal conditions (control). The germination paper was wetted, and 50 sterilized seeds were sown on one piece, which was then covered by another piece of wet paper towel [5]. Seed germination was defined as an observed radicle emergence of 0.5 cm. After a treatment incubation at 10 °C for 31 days and recovery at 15 °C for 7 days or a control incubation at 25 °C for 6 days, radicle traits were measured by an Epson Perfection V800 scanner, and Regent WinRHIZO (Canada) software was used for the data analysis. For each condition (treatment and control), three independent experiments were conducted per line.

The experimental field sites were located in XiangYang, Harbin and KeShan of Qiqihar in 2018. The seeds of the low-temperature treatment group were sown when the ground temperature was consistently approximately 6 °C. The normal sowing time was used as a control. The experiment involved a random block design with three replications, and the seed germination number and seedling length were determined.

Collection of phenotypic data

The germination rate (GR) was expressed as the percentage of germinating plants out of the total number of seeds used. The germ traits and radicle traits were measured using an Epson Perfection V800 scanner. The germ trait was the germ length (GL), and the radicle traits included the radicle length (RL), radicle surface area (RAS) and radicle volume (RV). All traits were measured as the mean of 10 seedlings. The germination index of plants was calculated from 19 to 31 days at 3-day intervals at 10 °C and was calculated daily, from 3 to 6 days at 25 °C. The calculated traits were included the germination index (GI), vitality index (VI) and simple vitality index (SVI), as follows:

$$ \mathrm{GI}=\sum \frac{\mathrm{Gt}}{\mathrm{Dt}} $$

where Gt is the number of germinating plants on a given day (Dt, days after sowing). VI = TL × GI, where TL is the total length of the seedling (including both the germ length and radicle length) on last day and GI is the germination index. SVI = GL × GR, where GR is the germination rate and GL is the germ length.

The field experiment statistics involving the germination rate were calculated on the last day of the XiangYang (XYGR) and KeShan (KSGR) tests, the germ length for XiangYang (XYGL) and KeShan (KSGL) was measured by a ruler, and a simple vitality index was calculated for both locations (XYRSVI and KSRSVI). Each experimental replicate consisted of 50 seedlings, and the mean of 10 seedlings was used for the XYGL and KSGL.

The relative performance of the 14 traits (RGR, RGL, RRL, RRSA, RRV, RGI, RVI, RSVI, XYRGR, XYRGL, XYRSVI, KSRGR, KSRGL, and KSRSVI) was calculated simply as the ratios of the mean values of measurements (n = 3) taken under low-temperature stress conditions and control conditions, and these ratios were used as indicators for low-temperature tolerance.

IBM SPSS Statistics version 20.0 software (https://www.ibm.com/support/pages/node/230551) was used to analyze the phenotypic data, including the range, mean, median, standard deviation, kurtosis and skewness of each relative trait. Phenotypic correlations were analyzed using IBM SPSS Statistics version 20.0 software and with the ‘PerformanceAnalytics’ package in R software. ANOVA for the 14 relative traits was performed for each association panel using the ‘lme4’ function in the R package ‘lme4’. The broad-sense heritability (H2) at each location was estimated using the following formula: H2 = δ2g/(δ2g + δ2/r). The Hb2 in cross-locations was estimated using the following formula: \( {H}_b^2={\updelta}^2\mathrm{g}/\left({\updelta}^2\mathrm{g}+{\updelta}^2\mathrm{g}\mathrm{e}/\mathrm{n}+{\updelta}^2/\mathrm{n}\mathrm{r}\right) \). δ2g, δ2ge, and δ2 are estimates of the genetic, G × E and error variances, and n and r are the number of environments and replications per environment, respectively [44].

Genotypic data and GWAS

Genotyping was carried out on the association panel using an Illumina Maize SNP50 BeadChip, which revealed 56,110 SNPs in the population and was filtered such that SNPs with a missing percentage > 20%, SNPs with a minor allele frequency (MAF) < 0.05, and SNPs with a heterozygosity > 20% were removed [25, 42, 43]. In total, 40,697 SNPs were used for the association analysis, with a MAF of > 0.05 in the population. Using the STRUCTURE 2.3 software, 7742 distributed SNP datasets were assessed for structural parameters [45]. ∆K was calculated using StructureHarvester [46]. The kinship information for 222 inbred lines was estimated using the software TASSEL 5.0.

The GWAS was performed in accordance with the MLM in TASSEL 5.0 [47], and the following 14 traits were used for association analysis: the RGR, RGL, RRL, RRSA, RRV, RGI, RVI, RSVI, XYRGR, XYRGL, XYRSVI, KSRGR, KSRGL and KSRSVI. GEC software (http://grass.cgs.hku.hk/gec/estimateB.php?function=Bonferroni) was used to calculate the effective number of markers (Ne) and to calculate the recommended threshold (0.05/Ne) as the basis for whether the 14 trait values were significantly correlated with a given SNP. Because a Bonferroni correction (0.05 / 23,398 = 2.140e-6) was too conservative (there were very few SNPs significantly associated with the 14 traits), a less stringent threshold of -log10(P) > 4 was used to detect significant association signals [25, 43, 48]. Manhattan plots were subsequently generated by the CMplot package in R software.

The linkage disequilibrium measurement parameter r2 was used to estimate the LD between all SNPs with less than 25% missing data for each chromosome via the software PopLDdecay 3.30 (https://github.com/BGI-shenzhen/PopLDdecay). All significantly associated SNPs on the same chromosome whose physical distance was less than the LD decay distance were defined as one site, and the range of each LD decay upstream and downstream of the SNP of the -log10(P) peak for each site was used to mine candidate genes. The B73 RefGen_V4 gene model from the maizeGDB website (http://www.maizegdb.org/) was used to map the loci and to retrieve genetic information.

RNA-seq and analysis of DEGs

The low-temperature-resistant line Zao 8–3 (referred to as 55) and a low-temperature-sensitive line Ji 853 (referred to as 102) were used in this study. Before imbibition, the seeds were sterilized and rinsed as described above. The seeds imbibed at normal temperature (25 °C) for 6 h and then were sampled at 2 h after imbibition on the germination paper under the ‘normal’ condition (25 °C) for CK_102 and CK_55 and under low-temperature conditions (10 °C) for LT_102 and LT_55. Embryos were isolated from the seeds of each sample, frozen in liquid nitrogen, and stored at − 80 °C for RNA extraction. There were three biological replicates for each treatment. The maize RNA samples were sent to LC-BIO (Hangzhou, China) for library construction and sequencing. Briefly, for each total RNA sample, the genomic DNA was removed by DNase I treatment, after which the mRNA was enriched using oligo (dT) magnetic beads and subsequently fragmented. Random hexamer primers were used for the synthesis of double-strand cDNA, which was further modified by end repair and 3′-polyadenylation. The fragments were amplified for sequencing via an Illumina NovaSeq 6000. The maize genome database ZmB73RefGenv4 (ftp://ftp.ncbi.nlm.nih.gov/genomes/all/GCF/000/005/005/GCF_000005005.2_B73_RefGen_v4) was used as the reference genome. The read count for each gene was obtained from the mapping results and normalized to fragments per kilobase of transcript per million mapped reads (FPKM). The default threshold for a significant difference in gene expression was set to |log2FC| ≥ 0.585 and p < 0.05. The raw sequence data are available in the NCBI GEO database under accession number GSE146666 (GSM4403922- GSM4403933).

Identification and annotation of candidate genes

GO enrichment analysis was used to identify all the GO terms of the genes whose expression was significantly enriched compared to the background from among the lists of DEGs and to filter the DEGs corresponding to specific biological functions. All DEGs were categorized and grouped based on their GO terms using the publicly available database at http://www.geneontology.org/. The gene numbers were calculated for every term. A hypergeometric test and GO Term Finder (http://www.yeastgenome.org/help/analyze/go-term-finder) were used to identify significantly enriched GO terms. The three categories, ‘biological process’, ‘cellular component’, and ‘molecular function’, were overrepresented and were further filtered using the Bonferroni multitest adjustment method and Fisher’s exact test [49, 50]. To further understand the biological functions of the DEGs, the significantly enriched metabolic or signal transduction pathways were identified using KEGG pathway enrichment analysis [51] (http://www.genome.ad.jp/kegg/kegg2.html). Genes were annotated using the maize GDB (http://www.maizegdb.org) and NCBI databases (http://www.ncbi.nlm.nih.gov/).

qRT-PCR-based validation of candidate genes

The total RNA samples used for transcriptome sequencing of the low-temperature-resistant line Zao 8–3 (referred to as 55) and the low-temperature-sensitive line Ji 853 (referred to as 102) were also subjected to qRT-PCR. For each sample, 500 ng of RNA was used for reverse transcription with a TransScript® One-Step gDNA Removal Kit in conjunction with cDNA Synthesis SuperMix (TransGen Biotech, Beijing, China). qRT-PCR was carried out using an Analytik Jena Real Time PCR Detection System, with the maize actin gene serving as a control. Gene-specific primers were designed using Primer 5.0 software (Additional file 10: Table S8). The 20 μl reaction mixtures consisted of 2 μl of diluted cDNA, 0.4 μl of reverse and forward primers, 7.2 μl of ddH2O and 10 μl of 2× TransStart® Tip Green qPCR SuperMix (TransGen Biotech, Beijing, China). The gene expression level was calculated using the relative 2^-∆∆CT analytical method [52]. For each sample, three biological replicates were included, with each biological replicate consisting of three technical repeats. For gene expression, the mean value of three replicates was used.

Availability of data and materials

The raw RNA-seq data are available in the NCBI GEO database under accession number GSE146666 (https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE146666).

Abbreviations

GWAS:

Genome-wide association study

SNP(s):

Single nucleotide polymorphism(s)

RNA-seq:

RNA-sequencing

LD:

Linkage disequilibrium

DEGs:

Differentially expressed genes

FC:

Fold change

GO:

Gene ontology

KEGG:

Kyoto encyclopedia of genes and genomes

qRT-PCR:

Quantitative real-time PCR

RGR:

Relative germination rate

RGL:

Relative germ length

RRL:

Relative radicle length

RRSA:

Relative radicle surface area

RRV:

Relative radicle volume

RGI:

Relative germination index

RVI:

Relative vitality index

RSVI:

Relative simple vitality index

XYRGR:

XiangYang relative germination rate

XYRGL:

XiangYang relative germ length

XYRSVI:

XiangYang relative simple vitality index

KSRGR:

KeShan relative germination rate

KSRGL:

KeShan relative germ length

KSRSVI:

KeShan relative simple vitality index

MLM:

Mixed-linear model

FPKM:

Fragments per kilobase of transcript per million mapped reads

SD:

Standard deviation

MAPK:

Mitogen-activated protein kinase

PH:

Pleckstrin homology

PIPs:

Phosphoinositide phosphates

PM:

Plasma membrane

References

  1. 1.

    Greaves JA. Improving suboptimal temperature tolerance in maize the search for variation. J Exp Bot. 1996;47(3):307–23.

    CAS  Google Scholar 

  2. 2.

    Rodríguez VM, Romay MC, Ordás A, Revilla P. Evaluation of European maize (Zea mays L.) germplasm under cold conditions. Genet Resour Crop Ev. 2010;57(3):329–35.

    Google Scholar 

  3. 3.

    Verheul MJ, Picatto C, Stamp P. Growth and development of maize (Zea mays L.) seedlings under chilling conditions in the field. Eur J Agron. 1996;5:31–43.

    Google Scholar 

  4. 4.

    Hope HJ, Maamari R, Séguin S, Hamilton RI, Dwyer LM, White RP. Low temperature emergence potential of short season corn hybrids grown under controlled environment and plot conditions. Can J Plant Sci. 1992;72(1):83–91.

    Google Scholar 

  5. 5.

    Li XH, Wang GH, Fu JJ, Li L, Jia GY, Ren LS, Lubberstedt T, Wang GY, Wang JH, Gu RL. 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–75.

    PubMed  PubMed Central  Google Scholar 

  6. 6.

    Mock JJ, Mcneill MJ. Cold tolerance of maize inbred lines adapted to various latitudes in North America. Crop Sci. 1979;19(2):239–42.

    Google Scholar 

  7. 7.

    Crosbie TM, Mock JJ, Smith OS. Comparison of gains predicted by several selection methods for cold tolerance traits of two maize populations. Crop Sci. 1980;20(5):649–55.

    Google Scholar 

  8. 8.

    Fracheboud Y, Jompuk C, Ribaut JM, Stamp P, Leipner J. Genetic analysis of cold-tolerance of photosynthesis in maize. Plant Mol Biol. 2004;56(2):241–53.

    CAS  PubMed  Google Scholar 

  9. 9.

    Rodríguez VM, Butrón A, Malvar RA, Ordas A, Revilla P. Quantitative trait loci for cold tolerance in the maize IBM population. Int J Plant Sci. 2008;169(4):551–6.

    Google Scholar 

  10. 10.

    Revilla P, Rodríguez VM, Ordás A, Rincent R, Charcosset A, Giauffret C, Melchinger AE, Schön C-C, Bauer E, Altmann T, et al. Association mapping for cold tolerance in two large maize inbred panels. BMC Plant Biol. 2016;16(1):127–36.

    PubMed  PubMed Central  Google Scholar 

  11. 11.

    Hu SD, Lübberstedt T, Zhao GW, Lee M. QTL mapping of low-temperature germination ability in the maize IBM Syn4 RIL population. PLoS One. 2016;11(3):e0152795.

    PubMed  PubMed Central  Google Scholar 

  12. 12.

    Huang J, Zhang JH, Li WZ, Hu W, Duan LC, Feng Y, Qiu FZ, Yue B. Genome-wide association analysis of ten cold tolerance indices at the germination and seeding stages in maize. J Integr Plant Biol. 2013;55(8):735–44.

    CAS  PubMed  Google Scholar 

  13. 13.

    Strigens A, Freitag NM, Gilbert X, Grieder C, Riedelsheimer C, Schrag TA, Messmer R, Melchinger AE. Association mapping for cold tolerance in elite flint and dent maize inbred lines evaluated in growth chamber and field experiments. Plant Cell Environ. 2013;36(10):1871–87.

    CAS  PubMed  Google Scholar 

  14. 14.

    Hu GH, Li Z, Lu YC, Li CX, Gong SC, Yan SQ, Li GL, Wang MQ, Ren HL, Guan HT, et al. Genome-wide association study identified multiple genetic loci on chilling resistance during germination in maize. Sci Rep. 2017;7:10840.

    PubMed  PubMed Central  Google Scholar 

  15. 15.

    Ma Y, Dai XY, Xu YY, Luo W, Zheng XM, Zeng DL, Pan YJ, Lin XL, Liu HH, Zhang DJ, et al. COLD1 confers chilling tolerance in rice. Cell. 2015;160(6):1209–21.

    CAS  PubMed  Google Scholar 

  16. 16.

    Zhang ZY, Li JJ, Pan YH, Li JL, Zhou L, Shi HL, Zeng YW, Guo HF, Yang SM, Zheng WW, et al. Natural variation in CTB4a enhances rice adaptation to cold habitats. Nat Commun. 2017;8:14788.

    CAS  PubMed  PubMed Central  Google Scholar 

  17. 17.

    Chen LP, Zhao Y, Xu SJ, Zhang ZY, Xu YY, Zhang JY, Chong K. OsMADS57 together with OsTB1 coordinates transcription of its target OsWRKY94 and D14 to switch its organogenesis to defense for cold adaptation in rice. New Phytol. 2018;218(1):219–31.

    CAS  PubMed  PubMed Central  Google Scholar 

  18. 18.

    Ding YL, Jia YX, Shi YT, Zhang XY, Song CP, Gong ZZ, Yang SH. OST1-mediated BTF3L phosphorylation positively regulates CBFs during plant cold responses. EMBO J. 2018;37(8):e98228.

    PubMed  PubMed Central  Google Scholar 

  19. 19.

    Berberich T, Kusano T. Cycloheximide induces a subset of low temperature-inducible in maize. Mol Gen Genet. 1997;254(3):275–83.

    CAS  PubMed  Google Scholar 

  20. 20.

    Kong XP, Lv W, Jiang SS, Zhang D, Cai GH, Pan JW, Li DQ. Genome-wide identification and expression analysis of calcium-dependent protein kinase in maize. BMC Genomics. 2013;14(1):433–47.

    CAS  PubMed  PubMed Central  Google Scholar 

  21. 21.

    Wang XY, Shan XH, Xue CM, Wu Y, Su SZ, Li SP, Liu HK, Jiang Y, Zhang YF, Yuan YP. Isolation and functional characterization of a cold responsive phosphatidylinositol transfer-associated protein, ZmSEC14p, from maize (Zea may L.). Plant Cell Rep. 2016;35(8):1671–86.

    CAS  PubMed  Google Scholar 

  22. 22.

    Berberich T, Sano H, Kusano T. Involvement of a MAP kinase, ZmMPK5, in senescence and recovery from low-temperature stress in maize. Mol Gen Genet. 1999;262(3):534–42.

    CAS  PubMed  Google Scholar 

  23. 23.

    Kusano T, Berberich T, Harada M, Suzuki N. A maize DNA-binding factor with a bZIP motif is induced by low temperature. Mol Gen Genet. 1995;248(5):507–17.

    CAS  PubMed  Google Scholar 

  24. 24.

    Qin F, Sakuma Y, Li J, Liu Q, Li YQ, Shinozaki K, Yamagushi-Shinozaki K. Cloning and functional analysis of a novel DREB1/CBF transcription factor involved in cold-responsive gene expression in Zea mays L. Plant Cell Physiol. 2004;45(8):1042–52.

    CAS  PubMed  Google Scholar 

  25. 25.

    Weng JF, Xie CX, Hao ZF, Wang JJ, Liu CL, Li MS, Zhang DG, Bai L, Zhang SH, Li XH. Genome-wide association study identifies candidate genes that affect plant height in Chinese elite maize (Zea mays L.) inbred lines. PLoS One. 2011;6(12):e29229.

    CAS  PubMed  PubMed Central  Google Scholar 

  26. 26.

    Nagel KA, Kastenholz B, Jahnke S, Dusschoten D, Aach T, Mühlich M, Truhn D, Scharr H, Terjung S, Walter A, et al. Temperature responses of roots: impact on growth, root system architecture and implications for phenotyping. Funct Plant Biol. 2009;36(11):947–59.

    CAS  Google Scholar 

  27. 27.

    Farooq M, Aziz T, Wahid A, Lee DJ, Siddique KH. Chilling tolerance in maize: agronomic and physiological approaches. Crop Pasture Sci. 2009;60(6):501–16.

    Google Scholar 

  28. 28.

    Guo J, Li CH, Zhang XQ, Li YX, Zhang DF, Shi YS, Song YC, Li Y, Yang DG, Wang TY. Transcriptome and GWAS analyses reveal candidate gene for seminal root length of maize seedlings under drought stress. Plant Sci. 2020;292:110380.

    CAS  PubMed  Google Scholar 

  29. 29.

    Tai HH, Lu X, Opitz N, Marcon C, Paschold A, Lithio A, Nettleton D, Hochholdinger F. Transcriptomic and anatomical complexity of primary, seminal, and crown roots highlight root type-specific functional diversity in maize (Zea mays L.). J Exp Bot. 2016;67(4):1123–35.

    CAS  PubMed  Google Scholar 

  30. 30.

    Rodríguez VM, Butrón A, Rady MOA, Soengas P, Revilla P. Identification of quantitative trait loci involved in the response to cold stress in maize (Zea mays L.). Mol Breed. 2014;33(2):363–71.

    Google Scholar 

  31. 31.

    Leipner J, Jompuk C, Camp KH, Stamp P, Fracheboud Y. QTL studies reveal little relevance of chilling-related seedling traits for yield in maize. Theor Appl Genet. 2008;116(4):555–62.

    PubMed  Google Scholar 

  32. 32.

    Ichimura K, Mizoguchi T, Yoshida R, Yuasa T, Shinozaki K. Various abiotic stresses rapidly activate Arabidopsis MAP kinases ATMPK4 and ATMPK6. Plant J. 2000;24(5):655–65.

    CAS  PubMed  Google Scholar 

  33. 33.

    Hasegawa PM, Bressan RA, Zhu JK, Bohnert HJ. Plant cellular and molecular responses to high salinity. Annu Rev Plant Physiol Plant Mol Biol. 2000;51:463–99.

    CAS  PubMed  Google Scholar 

  34. 34.

    Chinnusamy V, Schumaker K, Zhu JK. Molecular genetic perspectives on cross-talk and specificity in abiotic stress signalling in plants. J Exp Bot. 2004;55(395):225–36.

    CAS  PubMed  Google Scholar 

  35. 35.

    Jonak C, Kiegerl S, Ligterink W, Barker PJ, Huskisson NS, Hirt H. Stress signaling in plants: a mitogen-activated protein kinase pathway is activated by cold and drought. Proc Natl Acad Sci U S A. 1996;93(20):11274–9.

    CAS  PubMed  PubMed Central  Google Scholar 

  36. 36.

    Hoyos ME, Zhang SQ. Calcium-independent activation of salicylic acid-induced protein kinase and a 40-kilodalton protein kinase by hyperosmotic stress. Plant Physiol. 2000;122(4):1355–63.

    CAS  PubMed  PubMed Central  Google Scholar 

  37. 37.

    Li H, Ding YL, Shi YT, Zhang XY, Zhang SQ, Gong ZZ, Yang SH. MPK3-and MPK6-mediated ICE1 phosphorylation negatively regulates ICE1 stability and freezing tolerance in Arabidopsis. Dev Cell. 2017;43(5):630–42.

    CAS  PubMed  Google Scholar 

  38. 38.

    Zhao CZ, Wang PC, Si T, Hsu C-C, Wang L, Zayed O, Yu ZP, Zhu YF, Dong J, Tao WA, et al. MAP kinase cascades regulate the cold response by modulating ICE1 protein stability. Dev Cell. 2017;43(5):618–29.

    CAS  PubMed  PubMed Central  Google Scholar 

  39. 39.

    Zhang ZY, Li JH, Li F, Liu HH, Yang WS, Chong K, Xu YY. OsMAPK3 phosphorylates OsbHLH002/OsICE1 and inhibits its ubiquitination to activate OsTPP1 and enhances rice chilling tolerance. Dev Cell. 2017;43(6):731–43.

    CAS  PubMed  Google Scholar 

  40. 40.

    Palta JP, Whitaker BD, Weiss LS. Plasma membrane lipids associated with genetic variability in freezing tolerance and cold acclimation of solanum species. Plant Physiol. 1993;103(3):793–803.

    CAS  PubMed  PubMed Central  Google Scholar 

  41. 41.

    Marchler-Bauer A, Bo Y, Han L, He J, Lanczycki CJ, Lu S, Chitsaz F, Derbyshire MK, Geer RC, Gonzales NR, et al. CDD/SPARCLE: functional classification of proteins via subfamily domain architectures. Nucleic Acids Res. 2017;45(D1):D200–3.

    CAS  PubMed  Google Scholar 

  42. 42.

    Liu CL, Hao ZF, Zhang DG, Xie CX, Li MS, Zhang XC, Yong HJ, Zhang SH, Weng JF, Li XH. Genetic properties of 240 maize inbred lines and identity-by-descent segments revealed by high-density SNP markers. Mol Breed. 2015;35(7):146–57.

    CAS  Google Scholar 

  43. 43.

    Zhang CS, Zhou ZQ, Yong HJ, Zhang XC, Hao ZF, Zhang FJ, Li MS, Zhang DG, Li XH, Wang ZH, et al. Analysis of the genetic architecture of maize ear and grain morphological traits by combined linkage and association mapping. Theor Appl Genet. 2017;130(5):1011–29.

    CAS  PubMed  Google Scholar 

  44. 44.

    Knapp SJ, Stroup WW, Ross WM. Exact confidence intervals for heritability on a progeny mean basis. Crop Sci. 1985;25(1):192–4.

    Google Scholar 

  45. 45.

    Evanno G, Regnaut S, Goudet J. Detecting the number of clusters of individuals using the software STRUCTURE: a simulation study. Mol Ecol. 2005;14(8):2611–20.

    CAS  PubMed  PubMed Central  Google Scholar 

  46. 46.

    Earl DA. STRUCTURE HARVESTER: a website and program for visualizing STRUCTURE output and implementing the Evanno method. Conserv Genet Resour. 2012;4(2):359–61.

    Google Scholar 

  47. 47.

    Bradbury PJ, Zhang ZW, Kroon DE, Casstevens TM, Ramdoss Y, Buckler ES. TASSEL: software for association mapping of complex traits in diverse samples. Bioinformatics. 2007;23(19):2633–5.

    CAS  PubMed  Google Scholar 

  48. 48.

    Jia T, Wang L, Li J, Ma J, Cao Y, Lübberstedt T, Li HY. Integrating a genome-wide association study with transcriptomic analysis to detect genes controlling grain drying rate in maize (Zea may, L.). Theor Appl Genet. 2020;133(2):623–34.

    CAS  PubMed  Google Scholar 

  49. 49.

    Du Z, Zhou X, Ling Y, Zhang ZH, Su Z. agriGO: a GO analysis toolkit for the agricultural community. Nucleic Acids Res. 2010;38:W64–70.

    CAS  PubMed  PubMed Central  Google Scholar 

  50. 50.

    Han YC, Lv P, Hou SL, Li SY, Ji GS, Ma X, Du RH, Liu GQ. Combining next generation sequencing with bulked segregant analysis to fine map a stem moisture locus in sorghum (Sorghum bicolor L. Moench). PLoS One. 2015;10(5):e0127065.

    PubMed  PubMed Central  Google Scholar 

  51. 51.

    Kanehisa M, Araki M, Goto S, Hattori M, Hirakawa M, Itoh M, Katayama T, Kawashima S, Okuda S, Tokimatsu T, et al. KEGG for linking genomes to life and the environment. Nucleic Acids Res. 2008;36:D480–4.

    CAS  PubMed  Google Scholar 

  52. 52.

    Kenneth J, Schmittgen TD. Analysis of relative gene expression data using realtime quantitative PCR and the 2−ΔΔCT method. Methods. 2001;25:402–8.

    Google Scholar 

Download references

Acknowledgments

We would like to thank Prof. Xinhai Li and Jianfeng Weng for providing the plant materials.

Funding

This study was supported by the University Nursing Program for Young Scholars with Creative Talents in Heilongjiang Province (UNPYSCT-2017011). The funding agency had no role in the design of the study and collection, analysis and interpretation of data or in writing the manuscript.

Author information

Affiliations

Authors

Contributions

HZ and JZ performed the experiments and wrote the paper. QX, DW, HD, JH, XY, ZW, LZ, and LD performed the experiments. ZW and YZ designed the experiments and revised the paper. All authors in the “Contributions” section have read and approved the manuscript.

Corresponding authors

Correspondence to Zhenhua Wang or Yu Zhou.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary information

Additional file 1 Table S1.

ANOVA results and heritability of 14 traits of the lines composing the association panel.

Additional file 2 Table S2.

Correlation coefficients (r) between 14 phenotypic traits of 222 maize inbred lines under normal and low-temperature conditions.

Additional file 3 Figure S1.

Estimated population structure of 222 maize inbred lines.

Additional file 4 Table S3.

Mean of LD decay distance of the 10 chromosomes for r2 values which are equal to greater than 0.1 and 0.2.

Additional file 5 Figure S2.

Illustration of whole-genome LD in a panel of 222 maize lines.

Additional file 6 Table S4.

Functions of an additional 68 candidate genes identified within a LD decay of 26 associated SNPs.

Additional file 7 Table S5.

Twenty-six GO terms associated with eight candidate genes.

Additional file 8 Table S6.

SNPs or candidate genes that overlap with published QTLs.

Additional file 9 Table S7.

Germination rates of 222 maize inbred lines at 25 °C.

Additional file 10 Table S8.

Primer sequences used for quantitative real-time PCR.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Zhang, H., Zhang, J., Xu, Q. et al. Identification of candidate tolerance genes to low-temperature during maize germination by GWAS and RNA-seq approaches. BMC Plant Biol 20, 333 (2020). https://doi.org/10.1186/s12870-020-02543-9

Download citation

Keywords

  • Maize
  • Low-temperature
  • Germination
  • Genome-wide association study
  • Candidate genes