Next Article in Journal
Genome-Wide Identification and Expression Analysis of bZIP Family Genes in Stevia rebaudiana
Next Article in Special Issue
Advancements in Research on Prevention and Control Strategies for Maize White Spot Disease
Previous Article in Journal
The Small-Molecule E26-Transformation-Specific Inhibitor TK216 Attenuates the Oncogenic Properties of Pediatric Leukemia
Previous Article in Special Issue
Identification of Candidate QTLs and Genes for Ear Diameter by Multi-Parent Population in Maize
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Genes
Volume 14
Issue 10
10.3390/genes14101917
genes-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Transcriptome Analysis and QTL Mapping Identify Candidate Genes and Regulatory Mechanisms Related to Low-Temperature Germination Ability in Maize

by
Lei Du
3,†,
Xin Peng
2,†,
Hao Zhang
2,
Wangsen Xin
2,
Kejun Ma
2,
Yongzhong Liu
3 and
Guangcan Hu
1,*
1
Institute of Upland Food Crops, YiChang Academy of Agricultural Science, Yichang 443011, China
2
College of Plant Science and Technology, Huazhong Agricultural University, Wuhan 430070, China
3
Hubei Hongshan Laboratory, Wuhan 430070, China
*
Author to whom correspondence should be addressed.
†
These authors contributed equally to this work.
Genes 2023, 14(10), 1917; https://doi.org/10.3390/genes14101917
Submission received: 9 September 2023 / Revised: 29 September 2023 / Accepted: 5 October 2023 / Published: 8 October 2023
(This article belongs to the Special Issue Maize Molecular Genetics and Functional Genomics)
Browse Figures
Review Reports Versions Notes

Abstract

:
Low-temperature germination ability (LTGA) is an important characteristic for spring sowing maize. However, few maize genes related to LTGA were confirmed, and the regulatory mechanism is less clear. Here, maize-inbred lines Ye478 and Q1 with different LTGA were used to perform transcriptome analysis at multiple low-temperature germination stages, and a co-expression network was constructed by weighted gene co-expression network analysis (WGCNA). Data analysis showed that 7964 up- and 5010 down-regulated differentially expressed genes (DEGs) of Ye478 were identified at low-temperature germination stages, while 6060 up- and 2653 down-regulated DEGs of Q1 were identified. Gene ontology (GO) enrichment analysis revealed that ribosome synthesis and hydrogen peroxide metabolism were enhanced and mRNA metabolism was weakened under low-temperature stress for Ye478, while hydrogen peroxide metabolism was enhanced and mRNA metabolism was weakened for Q1. DEGs pairwise comparisons between the two genotypes found that Ye478 performed more ribosome synthesis at low temperatures compared with Q1. WGCNA analysis based on 24 transcriptomes identified 16 co-expressed modules. Of these, the MEbrown module was highly correlated with Ye478 at low-temperature stages and catalase and superoxide dismutase activity, and the MEred, MEgreen, and MEblack modules were highly correlated with Ye478 across low-temperature stages, which revealed a significant association between LTGA and these modules. GO enrichment analysis showed the MEbrown and MEred modules mainly functioned in ribosome synthesis and cell cycle, respectively. In addition, we conducted quantitative trait loci (QTL) analysis based on a doubled haploid (DH) population constructed by Ye478 and Q1 and identified a major QTL explanting 20.6% of phenotype variance on chromosome 1. In this QTL interval, we found three, four, and three hub genes in the MEbrown, MEred, and MEgreen modules, of which two hub genes (Zm00001d031951, Zm00001d031953) related to glutathione metabolism and one hub gene (Zm00001d031617) related to oxidoreductase activity could be the candidate genes for LTGA. These biological functions and candidate genes will be helpful in understanding the regulatory mechanism of LTGA and the directional improvement of maize varieties for LTGA.

1. Introduction

Maize (Zea mays L.) originated in tropical and subtropical regions and is inherently sensitive to low-temperature stress [1]. For maize, 25 °C is the optimum temperature for seed germination, while temperatures below 15 °C may cause low-temperature stress [2]. Low-temperature stress usually leads to seed emergence delay, slow seedling growth, and ultimately grain yield decreasing in many maize-production regions [3]. In China, more than 30 million hectares of spring sowing maize may encounter low-temperature stress during seed germination [4]. Therefore, it is necessary to improve low-temperature germination ability (LTGA) for cultivated varieties using breeding and genomic approaches.
Low-temperature tolerance is a complex quantitative trait involving numerous genes and mechanisms. In maize, several genes have been confirmed to be associated with this trait, including ZmDREB1/CBF, ZmRR1, ZmMPK8, bZIP68, and ZmICE1 [5,6,7,8]. Of these, the ZmDREB1/CBF transcription factor regulates the expression of the DRE (dehydration-responsive element) gene and results in higher low-temperature tolerance in seedlings [5]. ZmRR1 encoding a mitogen-activated protein kinase (MAPK) positively regulates the expression of ZmDREB1 to enhance low-temperature tolerance at the seedling stage, and ZmMPK8 is the negative regulator of ZmRR1 [6]. The bZIP68 transcription factor could repress the expression of ZmDREB1 and lead to poor low-temperature tolerance in seedlings [7]. ZmICE1 transcription factor positively regulates low-temperature tolerance both in seedling and seed germination, and the mechanism is blocking mitochondrial reactive oxygen species (ROS) bursting or directly regulating the expression of ZmDREB1s [8]. In addition, several candidate genes associated with LTGA were identified, for example, ZmbZIP113 and ZmTSAH1 [9], Zm00001d039219, and Zm00001d034319 [10]. These genes and mechanisms will be an important basis for genetic studies and directional variety improvement for LTGA in maize.
Except for that, these are other well-known regulatory mechanisms for low-temperature tolerance in plants. Plasma membrane functions, including fluidity and signal transduction, are widely considered to be correlated with plant low-temperature tolerance [11]. Membrane fluidity is affected by desaturated fatty acid content (DFC), and numerous genes have been found to modulate low-temperature responses by regulating DFC content [12,13,14]. Meanwhile, many signal channels in the plasma membrane are essential to transduce low-temperature sensing, and several genes involved in these channels have been confirmed to be associated with low-temperature tolerance [15,16,17]. In addition, low-temperature signal transmission intracellularly usually involves a variety of signal molecules and signal pathways. Plant hormones are well-known signal molecules for the regulation of low-temperature tolerance; for instance, increasing ABA biosynthesis could improve low-temperature tolerance [18], reducing free IAA content through overexpression of OsGH3-2 could increase low-temperature tolerance [19], and cytokinin [20,21], ethylene [22,23], and other hormones have been found to be correlated with low-temperature tolerance. Meanwhile, intracellular signal transmission usually relies on the MAPK cascades, which contain three kinases and amplify the signal step by step via phosphorylation. To date, numerous genes involved in MAPK function have been found to be associated with low-temperature tolerance in Arabidopsis [24,25], maize [6], and rice [26,27]. Transcription factors (TF) receive upstream signals and promote the transcription of downstream genes. Several TFs are reported to be involved in low-temperature tolerance, for example, maize IF genes ZmDREB1/CBF and bZIP68 [5,7], apple IF gene NAC [28], and rice IF genes bZIP73 and OsWRKY45–1 [29,30]. Except for these regulation mechanisms, there are a majority of other molecular mechanisms for low-temperature stress. Reactive oxygen (ROS) usually accumulates when plants encounter stress, such as low temperatures, and excessive ROS could lead to deleterious effects on plant cells, even cell death [31]. Therefore, the ROS-scavenging system was also an important mechanism for low-temperature tolerance. Of these, ROS-scavenging enzymes, including superoxide dismutase (SOD) and catalase (CAT), have been confirmed to increase low-temperature tolerance in plants [32,33]. Some chemicals, such as glutathione and ascorbic acid, could reduce ROS content to increase low-temperature tolerance in plants [34,35]. In addition, maintaining normal physiological function in cells under stress is important to withstand stress. Cell cycle is essential for the growth and development of plants, and some studies have found some genes regulate low-temperature tolerance by maintaining cell division in stress [36,37]. Similarly, ribosomes are essential for protein synthesis in living cells, and ribosomal synthesis and processing under low temperatures have been found to affect low-temperature tolerance, for example, ribosomal subunit synthesis-related genes RPS5 and CRASS [38,39] and rRNA processing-related gene STCH4/REIL2 [40]. All of these genes and regulatory mechanisms will be helpful in understanding the molecular basis and exploring excellent alleles for low-temperature tolerance.
Weighted gene co-expression network analysis (WGCNA) is usually used to research the gene expression patterns of multiple samples, which divide genes with similar expression patterns into the same modules. In addition, WGCNA can also analyze the relationship between modules and sample phenotypes and identify key regulatory modules and genes [41]. This method has been widely used to identify abiotic and biotic stress-responding genes in plants [41,42,43].
LTGA is an important characteristic for the early sowing of spring maize and high-altitude maize. The two inbred lines Ye478 and Q1 with different LTGAs were important maize germplasm sources for breeding in China. In this study, the objectives were to investigate the molecular basis of the LTGA difference between the two genotypes, explore the underling mechanism for LTGA, and identify key gene modules and candidate genes correlated with LTGA. To achieve these objectives, we conducted the transcriptome analysis for the two genotypes under low-temperature stress during seed germination, identified DEGs between the two genotypes and between the germination stages, and analyzed the biological function of these DEGs based on gene ontology (GO) enrichment analysis. In addition, we constructed a gene co-expression network based on WGCNA, identified key gene modules correlated with the samples and specific enzymatic activity, detected hub genes, and analyzed their underling biological function. In the end, we conducted QTL analysis based on the DH population constructed by Ye478 and Q1 and investigated candidate hub genes in the major QTL interval. These results will be used to understand the molecular basis of LTGA and further provide useful genes for directional breeding improvement.

2. Materials and Methods

2.1. Plant Materials

Two maize elite inbred lines were used to conduct transcriptome analysis in this study. The inbred line Ye478 was tolerant to low-temperature stress during seed germination, which has been widely applied over the past two decades. The inbred line Q1 was sensitive to low temperatures during seed germination and was bred from the cross of the commercial hybrid variety S7913 and the tropical landrace Tuxpeno. In addition, a DH population consisting of 123 DH lines constructed by Ye478 and Q1 as described previously [44] was used to perform QTL analysis for low-temperature germination rate.

2.2. Phenotypic Evaluation of Low-Temperature Tolerance during Germination

To ensure the quality and quantity of the seeds, the two inbred lines and the DH population were propagated in the spring of 2019 in Wuhan city. For the two inbred lines, the germination rate of the two inbred lines was evaluated with three replicates at low (12 °C) and optimum (25 °C) temperatures in the growth chamber, respectively. In brief, the seeds of each genotype were soaked in a 0.5% (v/v) hydrogen peroxide solution for 20 h at 25 °C. The seedling substrate was prepared using sterilized vermiculite with a 20% (w/w) distilled water content and put in a transparent germination box (19 × 13 × 12 cm). In each germination box, the thirty soaked seeds were uniformly placed on top of 2.0 cm of the seeding substrate and then covered with 1.0 cm. These prepared germination boxes were, respectively, located in the low and optimum temperature growth chambers, and 3000 lx light intensity was maintained through the whole germination stage. The germination rate for each germination box was investigated after emergence and repeatedly surveyed until the germination rate did not change. Seedlings are considered to emerge when they exceed 0.5 cm of the seeding bed. The germination rate for each germination box was the ratio of the seedling numbers to the seed numbers. For the DH population, the germination rates were evaluated with three replicates only at low temperatures in the growth chamber, and the germination rate was investigated 15 days after sowing. The phenotypic data of the two inbred lines and the DH population were organized in Excel 2010 and plotted and depicted using R software (version 4.3.1).

2.3. Measurement of Enzyme Activity

We measured the catalase (CAT) and superoxide dismutase (SOD) activity for the Ye478 and Q1 at four germination stages under low-temperature stress. The first sampling was taken after the seed soaking described in Section 2.2. The second, third, and fourth samplings were taken at the first, third, and fifth days after sowing, when the two genotypes were germinating at the low temperature conditions described in Section 2.2. At each sampling, the embryos of the seeds were dissected on ice and immediately frozen in a −80 °C refrigerator until enzyme activity was measured. Each sample consisted of two biological replicates, and each replicate contained ten embryos. The activities of CAT and SOD were measured with a commercial Boxbio kit (Beijing Boxbio Science Technology, Beijing, China).

2.4. RNA-Sequencing and Data Analysis

To explore the molecular mechanism of low-temperature tolerance during germination, the transcriptome analysis for Ye478 and Q1 was conducted at four stages under low-temperature stress. In brief, the seeds of the two genotypes were soaked in 0.5% (v/v) hydrogen peroxide solution for 20 h at 25 °C, and then the first sampling was conducted as the control, namely the 0d stage. The other soaked seeds of the two genotypes were, respectively, sown in transparent germination boxes described in Section 2.2, and then the germination boxes were put in a growth chamber at 12 °C and 3000 lx. Subsequently, the second, third, and fourth samplings were taken at the first, third, and fifth days after sowing and named as 1d, 3d, and 5d stages, respectively. At each sampling, the embryos of the seeds were dissected on ice and immediately frozen in a −80 °C refrigerator until RNA isolation. Each sample consisted of three biological replicates, and each replicate contained ten embryos.
Finally, a total of 24 samples from the two genotypes at four germination stages with three replicates were used for RNA sequencing. RNA was extracted using TRIzol reagent (Invitrogen, Life Technologies) following the standard protocol. The quality and quantity of the RNA were checked for integrity on an Agilent Technologies 2100 Bioanalyzer (Agilent Technologies, Santa Clara, CA, USA) according to their RNA Integrity Number (RIN) value. RNA-sequencing library construction was performed following the TruSeq RNA Sample Prep v2 protocol. The libraries were sequenced on an Illumina HiSeq 4000 platform with 150 bp paired-end reads at the Huazhong Agricultural University (Wuhan, China). The data analysis was conducted according to standard procedures; briefly, the raw data from the high-throughput sequencing were first filtered using cutadapt (v1.13) [45], and then clean reads were aligned to the maize reference genome (B73 AGPv4) using the HISAT2 (v2.1.0) tool [46]. Transcripts per kilobase of exon model per million mapped reads (TPM) values were calculated to represent the gene expression level. Differential expression analysis was performed between pairwise samples using the DESeq (2022) R package. The raw p-values were adjusted for multiple testing using the BH method [47], and genes were declared differentially expressed between two samples if the adjusted p-values were less than 0.05 and the fold-change was greater than 2. Finally, gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis was performed using the packages ‘clusterProfiler’ [48] in R software [49], and the GO term was declared significantly enriched if the q-value were less than 0.05.

2.5. WGCNA Co-Expression Network Construction Module Identification

Weighted gene co-expression network analysis was performed using the WGCNA package (Version 4.0.2) in R software [50]. This analysis is based on the TPM value of all 24 samples from the two genotypes during germination and filtering out the DEGs with a low expression (TPM < 1 in all samples). The basic process of WGCNA analysis was as follows: First, the soft thresholding power was determined using the R function pickSoftThreshold when the fitted curve first approached 0.9. Second, gene modules were detected using a hierarchical clustering tree according to the calculated soft threshold. Third, correlation analysis between modules and samples, between modules and activities of CAT and SOD, was calculated to explore specific modules. Finally, hub genes in target modules were extracted according to the criteria of |KME| > 0.9 (module eigengene-based connectivity) and |GS| > 0.2 (gene significance). Function annotations and GO annotations of the target hub genes were obtained from the NCBI website (https://www.ncbi.nlm.nih.gov/, accessed on 15 August 2023) and the Gramene website (https://maize-pangenome.gramene.org/, accessed on 15 August 2023).

2.6. QTL Analysis and Investigate Effects of the Target Genes

The DH population was genotyped using genotyping by sequencing (GBS) technology, and an ultra-high-density linkage map of the DH population was developed as previously described [44]. In brief, the ultra-high-density linkage maps contained 64,553 high-quality SNP and INDEL markers and were integrated into 1101 bins based on the linkage and crossover among these markers. QTL analysis was performed using the composite interval mapping method (CIM) in Windows QTL Cartographer v2.5 [51]. The average germination rate at 15 days after sowing under low temperatures was used for QTL analysis, and the LOD threshold was set to 3.0. The confidence interval (CI) of QTLs was defined as the 2-LOD interval flanking the QTL peak. Additionally, the effects of the target genes were investigated based on flanking markers of the gene locus, and a single-tailed student t-test was performed to check the difference. Significant marker-trait associations were declared at p < 0.05.

3. Results

3.1. Germination Rate Analysis of the Two Genotypes

Germination rates of Q1 and Ye478 at low and optimum temperatures were used to study low-temperature tolerance during germination. At optimum temperature conditions, the germination rate of both the Ye478 and Q1 could reach 100%, but the Ye478 was earlier to reach (Figure 1A,C). In the low-temperature condition, the germination rate of Ye478 reached 100% on the fifth day after sowing, while Q1 reached only 60% by the 14th day after sowing and no longer changed (Figure 1B,D). These results indicated that the Ye478 could be more tolerant to low-temperature stress during germination compared with the Q1.

3.2. RNA-Sequencing Analysis for the Two Genotypes during Seed Germination

To explore the molecular basis of LTGA, transcriptome analysis was conducted for Ye478 and Q1 at multiple low-temperature germination stages. In total, 24 libraries from the two genotypes at four germination stages with three biological replicates were constructed and analyzed. Based on RNA-sequencing analysis, 686.7 million clean paired-end (PE) reads were obtained, and the average number of PE reads per library was 28.6 million (Table S1). Genes with normalized reads lower than 0.5 TPM were removed from this analysis. Finally, 21632, 21699, 22173, and 23351 transcripts were found to be expressed at the 0d, 1d, 3d, and 5d stages of Ye478, respectively. In a similar vein, 21776, 21935, 22568, and 23506 transcripts were also identified in the samples from the respective stages of Q1. Approximately 35.4% of expressed genes were in the 0.5–5 TPM range, and 56.5% of expressed genes were in the range 5–100 TPM (Figure 2A).
Principal component analysis for the 24 samples was performed based on expression quantity, and the results are depicted in Figure 2B. The results showed that the biological replicates of each sample could be clearly clustered together, suggesting the reliability of RNA-sequencing data. Significant differences between Ye478 and Q1 were observed at each of the 0d, 1d, 3d, and 5d stages, and the variation tendency of the two genotypes was consistent through the germination stages, suggesting that the overall transcriptome profiling is determined by both genotype and germination stage.

3.3. Identification of DEGs at the Low-Temperature Germination Stage for Ye478

To identify genes with altered expression levels under low-temperature stress during germination for Ye478, we conducted pairwise comparisons between the control stage (0d) and the three low-temperature stress stages (1d, 3d, and 5d), respectively. In total, 6466 DEGs (4068 up- and 2398 down-regulated), 10519 DEGs (6227 up- and 4292 down-regulated), and 12527 DEGs (7532 up- and 4995 down-regulated) were identified at the 1d, 3d, and 5d stages compared with the 0d stage, respectively (Figure 3A,B). In order to further investigate the molecular mechanism of Ye748 responding to low-temperature stress, we combined these DEGs into one gene set and then conducted GO enrichment analysis. A total of 7964 up-regulated genes in this gene set were enriched in 98 biological process (BP) terms, 28 molecular function (MF) terms, and 50 biological process (CC) terms (q-value < 0.05), and the top 30 terms are shown in Figure 4A. These top terms are mainly involved in ribosome synthesis (“ribosome”, “structural constituent of ribosome”, “cytosolic ribosome”, “ribosomal subunit”, “large ribosomal subunit”, “cytosolic large ribosomal subunit”, “cytosolic small ribosomal subunit”, “small ribosomal subunit”), and hydrogen peroxide metabolism and function (“hydrolase activity”, “response to oxidative stress”, “hydrogen peroxide metabolic process”, “hydrogen peroxide catabolic process”, “reactive oxygen species metabolic process”, “hydrolase activity”, “acting on peroxide as acceptor”, “peroxidase activity”). In addition, 5010 down-regulated DEGs in this gene set were enriched in 39 BP terms, 14 MF terms, and 20 CC terms, most of which were related to the mRNA metabolism (Figure 4B). These results implied that ribosome synthesis and hydrogen peroxide metabolism were enhanced and mRNA metabolism was weakened at the low-temperature germination stage for Ye478.

3.4. Identification of DEGs at the Low-Temperature Germination Stage for the Q1

Similarly, we studied the DEGs between the 0d stage and the three low-temperature stress stages for Q1 and identified 2061 DEGs (1594 up- and 467 down-regulated), 5775 DEGs (3866 up- and 1909 down-regulated), and 6673 DEGs (4487 up- and 2186 down-regulated) at the 1d, 3d, and 5d stages compared with the 0d stage, respectively (Figure 3A,C). These DEGs were also combined into one gene set, which contained 6060 up-regulated genes and 2653 down-regulated genes. GO enrichment analysis showed that the up-regulated genes were enriched in 150 BP terms, 34 MF terms, and 28 CC terms, and the top 30 terms are shown in Figure 4C. Based on these results, the most significant terms involved in hydrogen peroxide metabolism and response to abiotic stimulus (“response to oxidative stress”, “hydrogen peroxide metabolic process”, “hydrogen peroxide catabolic process”, “reactive oxygen species metabolic process”, “acting on peroxide as acceptor”, “peroxidase activity”, “antioxidant activity”, “monooxygenase activity”, “response to oxygen-containing compounds”, “response to temperature stimulus”, “response to temperature stimulus”). Additionally, GO enrichment analysis was also performed for the down-regulated DEGs and significantly enriched three BP terms (mRNA metabolic process, mRNA processing, regulation of cell morphogenesis), three MF terms (cysteine dioxygenase activity, helicase activity, chromatin binding), and one CC term (plasma membrane protein complex) (Figure 4D). These results implied that hydrogen peroxide metabolism and response to abiotic stimulus were enhanced, and mRNA metabolism was weakened at the low-temperature germination stage for Q1.

3.5. Identification of DEGs between the Two Genotypes

Pairwise comparisons were performed to identify the DEGs at the same low-temperature germination stage between the two genotypes, and the numbers of DEGs are shown in Figure 3A. At the 0d stage, 4338 and 5340 genes were found to be significantly up- or down-regulated, respectively, in samples of the Ye478 compared with Q1. Of these, the up-regulated DEGs were significantly enriched in 22 BP terms and two MF terms (Figure 5A). In these BP terms, “response to heat”, “response to temperature stimulus”, “response to abiotic stimulus” and “cellular response to heat” were directly associated with temperature response, and “response to hydrogen peroxide”, “response to oxygen-containing compounds”, “response to abscisic acid”, and “response to alcohol” could be associated with tolerance. These results suggested differences in biological function related to temperature responses between the two genotypes. In addition, GO enrichment analysis was also conducted for down-regulated DEGs at the T0 stage. Based on the results, 111 BP terms, 49 MF terms, and 14 CC terms were significantly enriched, and the top 30 terms are shown in Figure 5B. These top terms are mainly involved in MF in terms of hydrolase and oxidoreductase activity and BP in terms of amino acid and glycometabolism, which could be involved in stress response.
Similarly, at the 1d stage, 3168 and 3976 genes were significantly up- or down-regulated in Ye478, respectively. At the 3d stage, the up- or down-regulated genes in Ye478 were 2727 and 3934, respectively. At the 5d stage, the up- or down-regulated genes in Ye478 were 2643 and 2396, respectively (Figure 3A). We combined DEGs from these three stages into one gene set, which included 4655 up-regulated DEGs and 5722 down-regulated DEGs. This gene set was considered to be the DEGs between the two genotypes under low-temperature stress during germination, and GO enrichment analysis was carried out. Based on the results, the up-regulated DEGs were significantly enriched in ten BP terms and three CC terms (Figure 5C). Most of these BP terms were associated with ribosome biogenesis and processing, and the three CC terms were “nucleolus”, “preribosome” and “nuclear lumen”, which suggested that the Ye478 has more ribosome synthesis under low-temperature stress during germination. For those down-regulated DEGs, a total of 106 BP terms, 58 MF terms, and nine CC terms were significantly enriched, and the top 30 terms are depicted in Figure 5D. Based on the results, most of the top BP and MF terms were related to response to oxidative stress and secondary metabolic processes, for example “response to oxidative stress”, “hydrogen peroxide metabolic process”, “reactive oxygen species metabolic process”, “oxidoreductase activity”, “peroxidase activity”, “secondary metabolic process”, “lipid biosynthetic process” and “response to external stimulus”. These results indicated that the Ye478 conducted more ribosome functions under low-temperature stress, while the Q1 conducted more stress responses.

3.6. WGCNA Analysis

In order to study the specific gene modules and genes that are highly associated with LTGA, we performed a WGCNA using the expression values (TPM) based on all 24 samples. After filtering out the DEGs with low expression (TPM < 1 in all samples), a total of 19753 genes were retained for the WGCNA analysis. The soft thresholding powers with β = 10 (R2 = 0.91) were selected for further cluster dendrogram analysis (Figure S1). Finally, these genes were divided into 16 modules (labeled with different colors), as shown in Figure 6A. Subsequently, the correlation analysis of these modules with samples was conducted (Figure 6B). Based on the results, we found that the MEbrown and MEred modules were highly positively correlated with the Ye478 at the low-temperature germination stage. And the MEgreen module was highly positively correlated with the Q1 across low-temperature stress, while the MEblack module was highly negatively correlated with the Q1. Similarly, some modules were significantly correlated with specific samples, for example, between the MEpink module with the sample Q0d (r = 0.78), between the MEyellow module with the sample Q0d (r = 0.79), between the MEtan module with the sample Y0d (r = −0.71), between the MEturquoise module with the sample Y0d (r = 0.95), between the MEsalmon module with the sample Q1d (r = 0.72), between the MEpurple module with the sample Q3d (r = 0.80), and between the MEblue module with the sample Q5d (r = 0.77). In addition, we also analyzed the relationship between these modules and oxidoreductase activity and found that the MEbrown module was significantly correlated with CAT and SOD activity (Figure 6B). All of these specific modules could be associated with LTGA.
Notably, the MEbrown module was highly correlated with the Ye478 at the low-temperature germination stage and was significantly correlated with CAT and SOD activity, which could be the key gene module for LTGA. GO enrichment analysis showed that a total of 2398 genes in the MEbrown module were significantly (q < 0.05) enriched in 26 BP terms, 11 MF terms, and six CC terms (Table 1). Of them, the BP terms were mainly involved in the regulation of cell cycle and phosphorylation, and the MF terms were mainly associated with trans-membrane transporters and cyclin-dependent protein kinases. KEGG enrichment analysis found that nine pathways were significantly enriched in the MEbrown module and mainly involved in the biosynthesis of plant secondary metabolites (Table S3). In addition, the MEred module was highly related to the Ye478 across the three low-temperature germination stages, which could be important to LGTA. GO and KEGG enrichment analyses found that 1090 genes were significantly enriched in 61 GO terms (33 BP terms, eight MF terms, and 20 CC terms) and two pathways, and almost all of them were related to ribosome synthesis (Table 1 and Table S3). Meanwhile, the MEgreen and MEblack modules had a continually negative and positive correlation with Ye478 at the four germination stages, which indicates that the genes in these two modules were continually down- and up-regulated in Ye478. GO and KEGG enrichment analyses found that the genes in the MEgreen module are mainly enriched in GO terms of ions trans-membrane transport and pathways of oxidative phosphorylation, carbon metabolism, and amino acid metabolism. For the MEblack module, we did not identify significantly enriched GO term and KEGG pathway.
Meanwhile, we analyzed hub genes in the MEbrown, MEred, MEgreen, and MEblack modules based on the principle that |KME| > 0.9 and |GS| > 0.3 and identified 414, 188, 239, and 137 hub genes, respectively. These four gene modules and hub genes could be essential to low-temperature tolerance during seed germination.

3.7. QTL Analysis and Investigate Effects of the Target Genes

QTL analysis was performed for the germination rate on the 15th day after sowing at low temperatures using the DH population. We identified a major QTL with LOD scores of 7.2, which explanted 20.6% of phenotypic variation and spanned a physical region from 191.3 to 209.5 Mb on chromosome 1 (Table 2). And this QTL could increase germination rate under low-temperature stress when the allele is derived from the Ye478 (Table 2). In this major QTL interval, we, respectively, found three, four, three, and zero hub genes from MEbrown, MEred, MEgreen, and MEred modules (Table 3), of which six genes (Zm00001d031951, Zm00001d031953, Zm00001d031651, Zm00001d031667, Zm00001d031617, and Zm00001d031992) have significant expression differences between the Ye478 and Q1 (Table S2). In addition, correlation analysis between flanking marker and low-temperature germination rate found significant differences at eight gene locus (Zm00001d031951, Zm00001d031953, Zm00001d031560, Zm00001d031640, Zm00001d031651, Zm00001d031667, Zm00001d031617, and Zm00001d031992), especially at the Zm00001d031951 and Zm00001d031953 locus (Figure 7). Of them, Zm00001d031951 and Zm00001d031953 are involved in the glutathione metabolic process; Zm00001d031651 is involved in the mitochondrial intermembrane space; Zm00001d031617 is involved in oxidoreductase activity; and Zm00001d031667 and Zm00001d031992 do not have function annotations (Table 3). All of these hub genes could be the candidate genes for LTGA in this study.

4. Discussion

Low-temperature germination tolerance is necessary for spring sowing maize, especially in high-altitude and high-latitude regions [52]. However, conventional genetic breeding of maize mainly focuses on increasing yield potential and usually results in the loss of specific alleles. Understanding molecular mechanisms and discovering excellent alleles will be helpful for the directional breeding improvement of cultivated varieties. In this study, Q1 is an important inbred line for hybrid corn production in Southwest and Northwest China, but it has poor germination ability at low-temperature stress, while another elite inbred line, Ye478 retains the characteristic of low-temperature germination tolerance (Figure 1). Therefore, exploring low-temperature tolerance alleles in Ye478 and transferring these alleles into Q1 is an effective strategy to create low-temperature tolerance germplasm resources and directionally improve the LTGA of maize varieties.
Low-temperature tolerance is a complex biological process that usually involves sophisticated regulatory mechanisms. Of these, ribosome synthesis and processing have been found to correlate with low-temperature tolerance in many plants, including Arabidopsis [53,54], maize [38,39,40], and rice [55,56]. In our results, at the low-temperature germination stage, a large number of ribosome synthesis and processing biological processes were found only in Ye478 but not in Q1 (Figure 4A). More importantly, based on the WGCNA and correlation analysis, we also found that one of the most significant modules, Mered, mainly functioned ribosome synthesis and processing (Table 1). These results suggested the importance of normal ribosome synthesis for adapting to low-temperature stress. In addition, numerous genes involved in the cell cycle were considered to be closely associated with low-temperature responses in many plants, such as maize [37] and rice [36,57]. In this study, we identified a significant module, MEbrown, using WGCNA analysis, which was highly correlated with Ye478 at the low-temperature germination stage and CAT and SOD activity. And GO analysis for the MEbrown module found most biological processes related to cell cycle (Table 1), which indicated that the normal cell cycle function of Ye478 could be essential to withstand low-temperature stress, and this function could be influenced by CAT and SOD activity. Except for these potential molecular mechanisms, ROS-scavenging enzyme systems, including SOD and CAT, were also important mechanisms for low-temperature tolerance [34]. Here, we identified numerous biological processes related to ROS metabolism in both Ye478 and Q1 when they were exposed to low temperature during germination, which suggests ROS could be a response molecule to low-temperature stress.
Exploiting low-temperature tolerance genes and understanding the molecular basis of LTGA will be helpful for directional breeding improvement. In maize, many genes have been demonstrated to confer low-temperature tolerance. For example, two IF genes (ZmDREB1/CB, bZIP68) and two MAPK genes (ZmRR1, ZmMPK8) regulate low-temperature tolerance only at the seedling stage, and another IF gene ZmICE1 regulate low-temperature tolerance both in the seedling and germination stages [5,6,7,8]. However, in DEGs between Ye478 and Q1, we did not identify significant enriched GO terms related to IF and MAPK. In contrast, we found ribosome synthesis, ROS-scavenging enzyme systems, and the cell cycle could be the main mechanisms for low-temperature tolerance in this study. Meanwhile, WGCNA and correlation analysis also found many hub genes related to ribosome synthesis and cell cycle (Table 1), and these hub genes have a higher expression level in the Ye478 at the low-temperature germination stage, which indicated the importance of these hub genes to withstand low-temperature stress. Additionally, we identified a major QTL associated with LTGA that spanned a physical region from 191.3 to 209.5 Mb on chromosome 1 based on a DH population. In this major QTL interval, we found three and four hub genes from MEbrown and MEred modules, respectively (Table 3). Of these, Zm00001d031951 and Zm00001d031953 are involved in the glutathione metabolic process, and glutathione is a well-known chemical for scavenging ROS [34,35], which further indicates the importance of these genes to low-temperature stress. The Zm00001d031617 is involved in oxidoreductase activity, which could be related to the ROS-scavenging system, which was an important mechanism for low-temperature tolerance [32,33]. Therefore, these genes could be the key candidate genes to withstand low temperatures in this study.

5. Conclusions

LTGA is essential to the earlier sowing of spring maize and high-altitude maize. To investigate the molecular mechanism and explore candidate genes for this trait, we performed transcriptome analysis for inbred lines Ye478 and Q1 at low-temperature germination stages and constructed a co-expression network with WGCNA. The results showed that 7964 up- and 5010 down-regulated DEGs of Ye478 were identified after low-temperature treatment, and 6060 up- and 2653 down-regulated DEGs of Q1 were identified. GO enrichment analysis found ribosome synthesis and hydrogen peroxide metabolism of Ye478 were enhanced under low-temperature stress, while hydrogen peroxide metabolism of Q1 was enhanced. DEGs pairwise comparisons between the two genotypes found Ye478 performed more ribosome functions at low-temperature stress compared with Q1. WGCNA analysis MEbrown and MEred modules could be target modules related to LTGA, and they are mainly involved in ribosome synthesis and cell cycle, respectively. In QTL analysis based on the DH population constructed by Ye478 and Q1, we identified a major QTL explaining 20.7% of the phenotype variance. In the QTL interval, ten hub genes in the MEbrown, MEred, MEgreen, and MEblack modules were significantly associated with low-temperature germination rate, of which two hub genes related to glutathione metabolism and one hub gene related to oxidoreductase activity could be the candidate genes for LTGA. These biological functions and key genes will be helpful in understanding the regulatory mechanism and directional improvement variety for LTGA.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/genes14101917/s1, Table S1: Base data for RNA sequencing; Table S2: Expression level of candidate genes in two genotypes at low-temperature germination stages; Table S3: The significantly enriched KEGG pathways of the target modules.

Author Contributions

Validation, X.P.; Formal analysis, X.P. and L.D.; Investigation, L.D., H.Z. and W.X.; Resources, K.M.; Writing—original draft, L.D. and X.P.; Supervision, Y.L. and G.H. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by the Major Project of Hubei Hongshan Laboratory (2022hszd003).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. He, F.; Shen, H.; Lin, C.; Fu, H.; Sheteiwy, M.S.; Guan, Y.; Huang, Y.; Hu, J. Transcriptome analysis of chilling-imbibed embryo revealed membrane recovery related genes in maize. Front. Plant Sci. 2016, 7, 1978. [Google Scholar] [CrossRef] [PubMed]
  2. Kołodziejczyk, I.; Kazmierczak, A.; Posmyk, M.M. Melatonin application modifies antioxidant defense and induces endoreplication in maize seeds exposed to chilling stress. Int. J. Mol. Sci. 2021, 22, 8628. [Google Scholar] [CrossRef] [PubMed]
  3. Huang, J.; Zhang, J.; Li, W.; Hu, W.; Duan, L.; Feng, Y.; Qiu, F.; 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, 735–744. [Google Scholar] [CrossRef]
  4. China Statistics Press. China Statistical Yearbook; China Statistics Press: Beijing, China, 2022. [Google Scholar]
  5. Qin, F.; Sakuma, Y.; Li, J.; Liu, Q.; Li, Y.Q.; Shinozaki, K.; Yamaguchi-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, 1042–1052. [Google Scholar] [CrossRef] [PubMed]
  6. Zeng, R.; Li, Z.; Shi, Y.; Fu, D.; Yin, P.; Cheng, J.K.; Jiang, C.F.; Yang, S.H. Natural variation in a type-A response regulator confers maize chilling tolerance. Nat. Commun. 2021, 12, 4713. [Google Scholar] [CrossRef] [PubMed]
  7. Li, Z.Y.; Fu, D.Y.; Wang, X.; Zeng, R.; Zhang, X.; Tian, J.G.; Zhang, S.S.; Yang, X.H.; Tian, F.; Lai, J.S.; et al. The transcription factor bZIP68 negatively regulates cold tolerance in maize. Plant Cell 2022, 34, 2833–2851. [Google Scholar] [CrossRef] [PubMed]
  8. Jiang, H.; Shi, Y.; Liu, J.; Li, Z.; Fu, D.Y.; Wu, S.F.; Li, M.Z.; Yang, Z.J.; Shi, Y.L.; Lai, J.S.; et al. Natural polymorphism of ZmICE1 contributes to amino acid metabolism that impacts cold tolerance in maize. Nat. Plants 2022, 8, 1176–1190. [Google Scholar] [CrossRef]
  9. 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, 24, 2629. [Google Scholar] [CrossRef] [PubMed]
  10. He, R.Y.; Zheng, J.J.; Chen, Y.; Pan, Z.Y.; Yang, T.; Zhou, Y.; Li, X.F.; Li, Y.Z.; Cheng, M.J.; Li, Y.; 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]
  11. Takahashi, D.; Li, B.; Nakayama, T.; Kawamura, Y.; Uemura, M. Plant plasma membrane proteomics for improving cold tolerance. Front. Plant Sci. 2013, 4, 90. [Google Scholar] [CrossRef]
  12. Martiniere, A.; Shvedunova, M.; Thomson, A.J.; Evans, N.H.; Penfield, S.; Runions, J.; McWatters, H.G. Homeostasis of plasma membrane viscosity in fluctuating temperatures. New Phytol. 2011, 192, 328–337. [Google Scholar] [CrossRef] [PubMed]
  13. Shi, J.; Cao, Y.; Fan, X.R.; Li, M.; Wang, Y.F.; Ming, F. A rice microsomal delta-12 fatty acid desaturase can enhance resistance to cold stress in yeast and Oryza sativa. Mol. Breed. 2012, 29, 743–757. [Google Scholar] [CrossRef]
  14. Wang, X.; Yu, C.; Liu, Y.; Yang, L.; Li, Y.; Yao, W.; Cai, Y.; Yan, X.; Li, S.; Cai, Y.; et al. GmFAD3A, A ω-3 fatty acid desaturase gene, enhances cold tolerance and seed germination rate under low temperature in rice. Int. J. Mol. Sci. 2019, 20, 3796. [Google Scholar] [CrossRef] [PubMed]
  15. Ma, Y.; Dai, X.; Xu, Y.; Luo, W.; Zheng, X.; Zeng, D.; Pan, Y.; Lin, X.; Liu, H.; Zhang, D.; et al. COLD1 confers chilling tolerance in Rice. Cell 2015, 160, 1209–1221. [Google Scholar] [CrossRef] [PubMed]
  16. Mori, K.; Renhu, N.; Naito, M.; Nakamura, A.; Shiba, H.; Yamamoto, T.; Suzaki, T.; Iida, H.; Miura, K. Ca2+-permeable mechanosensitive channels MCA1 and MCA2 mediate cold-induced cytosolic Ca2+ increase and cold tolerance in Arabidopsis. Sci. Rep. 2018, 8, 550. [Google Scholar] [CrossRef]
  17. Liu, Q.; Ding, Y.; Shi, Y.; Ma, L.; Wang, Y.; Song, C.; Wilkins, K.A.; Davies, J.M.; Knight, H.; Knight, M.R.; et al. The calcium transporter ANNEXIN1 mediates cold-induced calcium signaling and freezing tolerance in plants. EMBO J. 2021, 40, e104559. [Google Scholar] [CrossRef]
  18. Mega, R.; Meguro-Maoka, A.; Endo, A.; Shimosaka, E.; Murayama, S.; Nambara, E.; Seo, M.; Kanno, Y.; Abrams, S.R.; Sato, Y. Sustained low abscisic acid levels increase seedling vigor under cold stress in rice (Oryza sativa L.). Sci. Rep. 2015, 5, 13819. [Google Scholar] [CrossRef]
  19. Du, H.; Wu, N.; Fu, J.; Wang, S.; Li, X.; Xiao, J.; Xiong, L. A GH3 family member, OsGH3-2, modulates auxin and abscisic acid levels and differentially affects drought and cold tolerance in rice. J. Exp. Bot. 2012, 63, 6467–6480. [Google Scholar] [CrossRef]
  20. Zhu, J.; Zhang, K.X.; Wang, W.S.; Gong, W.; Liu, W.C.; Chen, H.G.; Xu, H.H.; Lu, Y.T. Low temperature inhibits root growth by reducing auxin accumulation via ARR1/12. Plant Cell Physiol. 2015, 56, 727–736. [Google Scholar] [CrossRef]
  21. Zwack, P.J.; Compton, M.A.; Adams, C.I.; Rashotte, A.M. Cytokinin response factor 4 (CRF4) is induced by cold and involved in freezing tolerance. Plant Cell Rep. 2016, 35, 573–584. [Google Scholar] [CrossRef] [PubMed]
  22. Bolt, S.; Zuther, E.; Zintl, S.; Hincha, D.K.; Schmulling, T. ERF105 is a transcription factor gene of Arabidopsis thaliana required for freezing tolerance and cold acclimation. Plant Cell Environ. 2017, 40, 108–120. [Google Scholar] [CrossRef] [PubMed]
  23. Catala, R.; Salinas, J. The Arabidopsis ethylene overproducer mutant eto1-3 displays enhanced freezing tolerance. Plant Signal. Behav. 2015, 10, e989768. [Google Scholar] [CrossRef] [PubMed]
  24. Li, H.; Ding, Y.; Shi, Y.; Zhang, X.; Zhang, S.; Gong, Z.; Yang, S. MPK3- and MPK6-Mediated ICE1 Phosphorylation Negatively Regulates ICE1 Stability and Freezing Tolerance in Arabidopsis. Dev. Cell 2017, 43, 630–642. [Google Scholar] [CrossRef] [PubMed]
  25. Zhao, C.; Wang, P.; Si, T.; Hsu, C.C.; Wang, L.; Zayed, O.; Yu, Z.; Zhu, Y.; Dong, J.; Tao, W.A.; et al. MAP Kinase Cascades Regulate the Cold Response by Modulating ICE1 Protein Stability. Dev. Cell 2017, 43, 618–629. [Google Scholar] [PubMed]
  26. Zhang, Z.; Li, J.H.; Li, F.; Liu, H.H.; Yang, W.S.; Chong, K.; Cu, Y.Y. OsMAPK3 phosphorylates OsbHLH002/OsICE1 and inhibits its ubiquitination to activate OsTPP1 and enhances rice chilling tolerance. Dev. Cell 2017, 43, 731–743. [Google Scholar] [CrossRef] [PubMed]
  27. Jia, M.R.; Meng, X.B.; Song, X.G.; Zhang, D.H.; Kou, L.Q.; Zhang, J.H.; Jing, Y.H.; Liu, G.F.; Liu, H.H.; Huang, X.H.; et al. Chilling-induced phosphorylation of IPA1 by OsSAPK6 activates chilling tolerance responses in rice. Cell Discov. 2022, 8, 71. [Google Scholar] [CrossRef] [PubMed]
  28. An, J.P.; Li, R.; Qu, F.J.; You, C.X.; Wang, X.F.; Hao, Y.J. An apple NAC transcription factor negatively regulates cold tolerance via CBF-dependent pathway. J. Plant Physiol. 2018, 221, 74–80. [Google Scholar] [CrossRef]
  29. Liu, C.T.; Schläppi, M.R.; Mao, B.G.; Wang, W.; Wang, A.J.; Chu, C.C. The bZIP73 transcription factor controls rice cold tolerance at the reproductive stage. Plant Biotechnol. J. 2019, 17, 1834–1849. [Google Scholar] [CrossRef] [PubMed]
  30. Zhang, H.J.; Wu, T.; Li, Z.; Huang, K.; Kim, N.E.; Ma, Z.M.; Kwon, S.W.; Jiang, W.Z.; Du, X.L. OsGATA16, a GATA Transcription Factor, Confers Cold Tolerance by Repressing OsWRKY45–1 at the Seedling Stage in Rice. Rice 2021, 14, 42. [Google Scholar] [CrossRef] [PubMed]
  31. You, J.; Chan, Z. ROS regulation during abiotic stress responses in crop plants. Front. Plant Sci. 2015, 6, 1092. [Google Scholar] [CrossRef] [PubMed]
  32. Sato, Y.; Masuta, Y.; Saito, K.; Ozawa, K. Enhanced chilling tolerance at the booting stage in rice by transgenic overexpression of the ascorbate peroxidase gene, OsAPXa. Plant Cell Rep. 2011, 30, 399–406. [Google Scholar] [PubMed]
  33. Xu, J.; Duan, X.G.; Yang, J.; Beeching, J.R.; Zhang, P. Coupled expression of Cu/Zn-superoxide dismutase and catalase in cassava improves tolerance against cold and drought stresses. Plant Signal. Behav. 2013, 8, e24525. [Google Scholar] [CrossRef] [PubMed]
  34. Kocsy, C.; Galiba, C.; Brunold, C. Role of glutathione in adaptation and signalling during chilling and cold acclimation in plants. Physiol. Plant. 2008, 113, 158–164. [Google Scholar] [CrossRef] [PubMed]
  35. Zhang, Y.; Ming, R.H.; Khan, M.; Wang, Y.; Dahro, B.; Xiao, W.; Li, C.L.; Liu, J.H. ERF9 of Poncirus trifoliata (L.) Raf. undergoes feedback regulation by ethylene and modulates cold tolerance via regulating a glutathione S-transferase U17 gene. Plant Biotechnol. J. 2022, 20, 183–200. [Google Scholar] [CrossRef]
  36. Ma, Q.B.; Dai, X.Y.; Xu, Y.Y.; Guo, J.; Liu, Y.J.; Chen, N.; Xiao, J.; Zhang, D.J.; Xu, Z.H.; Zhang, X.S.; et al. Enhanced tolerance to chilling stress in OsMYB3R-2 transgenic rice is mediated by alteration in cell cycle and ectopic expression of stress genes. Plant Physiol. 2009, 150, 244–256. [Google Scholar] [CrossRef] [PubMed]
  37. Rymen, B.; Fiorani, F.; Kartal, F.; Vandepoele, K.; Inzé, D.; Beemster, G.T. Cold nights impair leaf growth and cell cycle progression in maize through transcriptional changes of cell cycle genes. Plant Physiol. 2007, 143, 1429–1438. [Google Scholar] [CrossRef]
  38. Zhang, J.X.; Yuan, H.; Yang, Y.; Fish, T.; Lyi, S.M.; Thannhauser, T.M.; Zhang, L.; Li, L. Plastid ribosomal protein S5 is involved in photosynthesis, plant development, and cold stress tolerance in Arabidopsis. J. Exp. Bot. 2016, 67, 2731–2744. [Google Scholar] [CrossRef] [PubMed]
  39. Pulido, P.; Zagari, N.; Manavski, N.; Gawronski, P.; Matthes, A.; Scharff, L.B.; Meurer, J.; Leister, D. CHLOROPLAST RIBOSOME ASSOCIATED Supports Translation under Stress and Interacts with the Ribosomal 30S Subunit. Plant Physiol. 2018, 177, 1539–1554. [Google Scholar] [CrossRef] [PubMed]
  40. Yu, H.S.; Kong, X.F.; Huang, H.; Wu, W.W.; Park, J.; Yun, D.J.; Lee, B.H.; Shi, H.Z.; Zhu, J.K. STCH4/REIL2 Confers Cold Stress Tolerance in Arabidopsis by Promoting rRNA Processing and CBF Protein Translation. Cell Rep. 2020, 30, 229–242.e5. [Google Scholar] [CrossRef] [PubMed]
  41. Ma, L.; Zhang, M.; Chen, J.; Qing, C.; He, S.; Zou, C.; Yuan, G.; Yang, C.; Peng, H.; Pan, G. GWAS and WGCNA uncover hub genes controlling salt tolerance in maize (Zea mays L.) seedlings. Theor. Appl. Genet. 2021, 134, 3305–3318. [Google Scholar] [CrossRef] [PubMed]
  42. Zhu, M.; Xie, H.; Wei, X.; Dossa, K.; Yu, Y.; Hui, S.; Tang, G.; Zeng, X.; Yu, Y.; Hu, P.; et al. WGCNA analysis of salt-responsive core transcriptome identifies novel hub genes in rice. Genes 2019, 10, 719. [Google Scholar] [CrossRef] [PubMed]
  43. Yu, T.; Zhang, J.; Cao, J.; Ma, X.; Li, W.; Yang, G. Hub gene mining and co-expression network construction of low-temperature response in maize of seedling by WGCNA. Genes 2023, 14, 1598. [Google Scholar] [CrossRef] [PubMed]
  44. Du, L.; Yu, F.; Zhang, H.; Wang, B.; Ma, K.; Yu, C.; Xin, W.; Huang, X.; Liu, Y.; Liu, K. Genetic mapping of quantitative trait loci and a major locus for resistance to grey leaf spot in maize. Theor. Appl. Genet. 2020, 133, 2521–2533. [Google Scholar] [CrossRef] [PubMed]
  45. Martin, M. Cutadapt removes adapter sequences from high-throughput sequencing reads. EMBnet. J. 2011, 17, 10–12. [Google Scholar] [CrossRef]
  46. Kim, D.; Langmead, B.; Salzberg, S.L. HISAT: A fast spliced aligner with low memory requirements. Nat. Methods 2015, 12, 357–360. [Google Scholar] [CrossRef]
  47. Love, M.I.; Huber, W.; Anders, S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014, 15, 550. [Google Scholar] [CrossRef] [PubMed]
  48. Yu, G.C.; Wang, L.G.; Han, Y.Y.; He, Q.Y. clusterProfiler: An R package for comparing biological themes among gene clusters. OMICS J. Integr. Biol. 2012, 16, 284–287. [Google Scholar] [CrossRef]
  49. R Core Team. R: A Language and Environment for Statistical Computing; R Foundation for Statistical Vienna, Austria. 2014. Available online: http://www.R-project.org/ (accessed on 16 June 2023).
  50. Langfelder, P.; Horvath, S. WGCNA: An R package for weighted correlation network analysis. BMC Bioinform. 2008, 9, 559. [Google Scholar] [CrossRef]
  51. Wang, S.; Basten, C.; Zeng, Z. Windows QTL Cartographer 2.5; Department of Statistics, North Carolina State University: Raleigh, NC, USA, 2012. [Google Scholar]
  52. Liu, Z.J.; Hubbard, K.G.; Lin, X.M.; Yang, X.G. Negative effects of climate warming on maize yield are reversed by the changing of sowing date and cultivar selection in Northeast China. Glob. Chang. Biol. 2013, 19, 3481–3492. [Google Scholar] [CrossRef]
  53. Martinez-Seidel, F.; Suwanchaikasem, P.; Gentry-Torfer, D.; Rajarathinam, R.; Ebert, A.; Erban, A.; Firmino, A.A.P.; Nie, S.; Leeming, M.G.; Williamson, N.A.; et al. Remodelled ribosomes synthesise a specific proteome in proliferating plant tissue during cold. bioRxiv 2022. [Google Scholar] [CrossRef]
  54. Martinez-Seidel, F.; Beine-Golovchuk, O.; Hsieh, Y.-C.; Eshraky, K.E.; Gorka, M.; Cheong, B.-E.; Jimenez-Posada, E.V.; Walther, D.; Skirycz, A.; Roessner, U.; et al. Spatially enriched paralog rearrangements argue functionally diverse ribosomes arise during cold acclimation in Arabidopsis. Int. J. Mol. Sci. 2021, 22, 6160. [Google Scholar] [CrossRef]
  55. Wang, Z.; Sun, J.; Zu, X.F.; Gong, J.; Deng, H.J.; Hang, R.; Zhang, X.F.; Liu, C.Y.; Deng, X.; Luo, L.L.; et al. Pseudouridylation of chloroplast ribosomal RNA contributes to low temperature acclimation in rice. New Phytol. 2022, 236, 1708–1720. [Google Scholar] [CrossRef] [PubMed]
  56. Wang, Y.L.; Wang, C.M.; Zheng, M.; Lyu, J.; Xu, Y.; Li, X.H.; Niu, M.; Long, W.H.; Wang, D.; Wang, H.Y.; et al. WHITE PANICLE1, a Val-tRNA synthetase regulating chloroplast ribosome biogenesis in rice, is essential for early chloroplast development. Plant Physiol. 2016, 170, 2110–2123. [Google Scholar] [CrossRef] [PubMed]
  57. Xu, P.P.; Cai, W.M. RAN1 is involved in plant cold resistance and development in rice (Oryza sativa). J. Exp. Bot. 2014, 65, 3277–3287. [Google Scholar] [CrossRef] [PubMed]
Genes 14 01917 g001
Figure 1. Germination rates of Ye478 and Q1 under optimum and low temperatures during germination. (A) germination rates under optimum temperature. (B) germination rates under low temperature. (C,D) the final germination performance under optimum and low temperatures, respectively.
Figure 1. Germination rates of Ye478 and Q1 under optimum and low temperatures during germination. (A) germination rates under optimum temperature. (B) germination rates under low temperature. (C,D) the final germination performance under optimum and low temperatures, respectively.
Genes 14 01917 g001
Genes 14 01917 g002
Figure 2. Global gene expression profiling of the two genotypes in all germination stages. (A) Numbers of detected transcripts in each sample. (B) Principal component analysis of the RNA-sequencing data. And Q0d, Q1d, Q3d, and Q5d represent germination stages of the Q1 at 0d, 1d, 3d, and 5d stages, respectively; Y0d, Y1d, Y3d, and Y5d represent germination stages of the Ye478 at 0d, 1d, 3d, and 5d stages, respectively.
Figure 2. Global gene expression profiling of the two genotypes in all germination stages. (A) Numbers of detected transcripts in each sample. (B) Principal component analysis of the RNA-sequencing data. And Q0d, Q1d, Q3d, and Q5d represent germination stages of the Q1 at 0d, 1d, 3d, and 5d stages, respectively; Y0d, Y1d, Y3d, and Y5d represent germination stages of the Ye478 at 0d, 1d, 3d, and 5d stages, respectively.
Genes 14 01917 g002
Genes 14 01917 g003
Figure 3. Numbers of DEGs of the two genotypes at each germination stage. (A) Numbers of DEGs between the two genotypes at the same germination stage (on the diagonal line), among different stages for Q1 (above the diagonal line), and among different stages for Ye478 (below the diagonal line), and the red and black numbers are the up- and down-regulated DEGs when the sample in the later stage is compared with the sample in the earlier stage. (B,C) the overlapping numbers of the up-and-down-regulated DEGs between Ye478 and Q1 at the same germination stages, respectively.
Figure 3. Numbers of DEGs of the two genotypes at each germination stage. (A) Numbers of DEGs between the two genotypes at the same germination stage (on the diagonal line), among different stages for Q1 (above the diagonal line), and among different stages for Ye478 (below the diagonal line), and the red and black numbers are the up- and down-regulated DEGs when the sample in the later stage is compared with the sample in the earlier stage. (B,C) the overlapping numbers of the up-and-down-regulated DEGs between Ye478 and Q1 at the same germination stages, respectively.
Genes 14 01917 g003
Genes 14 01917 g004
Figure 4. GO enrichment analysis of DEGs under low-temperature stress compared with the control stage. (A) The top 30 GO terms of up-regulated DEGs for Ye478. (B) The top 30 GO terms of down-regulated DEGs for Ye478. (C) The top 30 GO terms of up-regulated DEGs for Q1. (D) The significant enriched GO terms of down-regulated DEGs for Q1.
Figure 4. GO enrichment analysis of DEGs under low-temperature stress compared with the control stage. (A) The top 30 GO terms of up-regulated DEGs for Ye478. (B) The top 30 GO terms of down-regulated DEGs for Ye478. (C) The top 30 GO terms of up-regulated DEGs for Q1. (D) The significant enriched GO terms of down-regulated DEGs for Q1.
Genes 14 01917 g004
Genes 14 01917 g005
Figure 5. GO enrichment analysis of DEGs between the two genotypes at the same germination stage. (A) The significant enriched GO terms of up-regulated DEGs at T0 stage. (B) The top 30 GO terms of down-regulated DEGs at the T0 stage. (C) The significant enriched GO terms of up-regulated DEGs at the low-temperature stages. (D) The top 30 GO terms of down-regulated DEGs at low-temperature stages.
Figure 5. GO enrichment analysis of DEGs between the two genotypes at the same germination stage. (A) The significant enriched GO terms of up-regulated DEGs at T0 stage. (B) The top 30 GO terms of down-regulated DEGs at the T0 stage. (C) The significant enriched GO terms of up-regulated DEGs at the low-temperature stages. (D) The top 30 GO terms of down-regulated DEGs at low-temperature stages.
Genes 14 01917 g005
Genes 14 01917 g006
Figure 6. WGCNA of genes in Ye478 and Q1 at each germination stage. (A) Hierarchical dendrogram reveals co-expression modules identified by WGCNA. The major tree branches constitute 16 modules, labeled with different colors. (B) The correlation coefficient and correlation significance between each module with sample and oxidoreductase activity.
Figure 6. WGCNA of genes in Ye478 and Q1 at each germination stage. (A) Hierarchical dendrogram reveals co-expression modules identified by WGCNA. The major tree branches constitute 16 modules, labeled with different colors. (B) The correlation coefficient and correlation significance between each module with sample and oxidoreductase activity.
Genes 14 01917 g006
Genes 14 01917 g007
Figure 7. Low-temperature germination effects of candidate genes in the DH population.
Figure 7. Low-temperature germination effects of candidate genes in the DH population.
Genes 14 01917 g007
Table 1. Top GO terms of the target modules.
Table 1. Top GO terms of the target modules.
OntologyIDDescriptionq-Value
MEbrown:
BPGO:0043549regulation of kinase activity1.5 × 105
BPGO:0045859regulation of protein kinase activity1.5 × 105
BPGO:0000079regulation of cyclin-dependent protein serine/threonine kinase activity1.5 × 105
BPGO:1904029regulation of cyclin-dependent protein kinase activity1.5 × 105
BPGO:0001932regulation of protein phosphorylation4.1 × 105
BPGO:0042325regulation of phosphorylation4.1 × 105
BPGO:0071900regulation of protein serine/threonine kinase activity4.1 × 105
MFGO:0046983protein dimerization activity9.2 × 105
BPGO:0007049cell cycle1.8 × 104
BPGO:0051338regulation of transferase activity2.3 × 104
BPGO:0051726regulation of cell cycle2.3 × 104
BPGO:0040008regulation of growth1.0 × 103
BPGO:0000278mitotic cell cycle1.2 × 103
BPGO:0048638regulation of developmental growth1.3 × 103
BPGO:0044770cell cycle phase transition1.3 × 103
MEred:
CCGO:0005840ribosome2.8 × 10140
MFGO:0003735structural constituent of ribosome8.4 × 10137
MFGO:0005198structural molecule activity4.8 × 10124
CCGO:0022626cytosolic ribosome9.8 × 1089
CCGO:0044391ribosomal subunit2.5 × 1085
CCGO:1990904ribonucleoprotein complex3.2 × 1072
CCGO:0022625cytosolic large ribosomal subunit6.5 × 1056
CCGO:0015934large ribosomal subunit6.1 × 1054
CCGO:0022627cytosolic small ribosomal subunit1.3 × 1032
CCGO:0015935small ribosomal subunit7.0 × 1030
BPGO:0042254ribosome biogenesis1.5 × 1025
BPGO:0022613ribonucleoprotein complex biogenesis3.0 × 1023
BPGO:0042273ribosomal large subunit biogenesis1.5 × 1017
BPGO:0042255ribosome assembly7.8 × 1017
BPGO:0002181cytoplasmic translation6.7 × 1016
MEgreen:
CCGO:0016469proton-transporting two-sector ATPase complex8.6 × 108
CCGO:0033178proton-transporting two-sector ATPase complex7.4 × 107
CCGO:0033176proton-transporting V-type ATPase complex9.7 × 107
CCGO:0033180proton-transporting V-type ATPase, V1 domain1.3 × 105
CCGO:0098796membrane protein complex1.6 × 103
CCGO:1990204oxidoreductase complex4.4 × 102
MFGO:0044769ATPase activity, coupled to transmembrane movement of ions1.7 × 107
MFGO:0046961proton-transporting ATPase activity, rotational mechanism1.7 × 107
MFGO:0019829ATPase-coupled cation transmembrane transporter activity1.3 × 104
MFGO:0042625ATPase-coupled ion transmembrane transporter activity1.3 × 104
MFGO:0015078proton transmembrane transporter activity8.3 × 104
MFGO:0016627oxidoreductase activity3.9 × 102
Table 2. QTLs for germination rate at 15th days after sowing.
Table 2. QTLs for germination rate at 15th days after sowing.
Chr.Peak (cM)LODAddR2 (%)2 LOD (cM)Location (Mb)
1105.05.0−23.713.893.2–106.1118.5–174.8
1128.87.219.020.6117.8–139.7191.3–209.5
2114.23.49.48.9109.3–119.1223.6–229.7
Table 3. The hub genes of target modules in the major QTL interval.
Table 3. The hub genes of target modules in the major QTL interval.
GeneIDKMEGSNR AnnotationGO TermLocation (Chr. 1)
MEbrown:
Zm00001d0316550.900.41uncharacterized protein LOC100281072regulation of transcription197,545,631–197,549,052
Zm00001d0319510.980.63uncharacterized protein LOC100272476glutathione metabolic process207,697,908–207,854,723
Zm00001d0319530.920.69uncharacterized protein LOC100272476glutathione metabolic process207,932,227–207,935,640
MEred:
Zm00001d0314290.960.84pfkB-like carbohydrate kinase family proteinstructural constituent of ribosome189,856,343–189,862,099
Zm00001d0315600.950.8160S ribosomal protein L32structural constituent of ribosome194,266,307–194,270,465
Zm00001d0316400.900.66hypothetical protein SETIT_2G267400v2*197,108,702–197,111,394
Zm00001d0316510.920.88Mitochondrial import inner membrane translocase subunit Tim8mitochondrial intermembrane space197,393,450–197,395,554
MEgreen:
Zm00001d0316670.92−0.42hypothetical protein Zm00014a_031611*196,224,574–196,228,047
Zm00001d0316170.92−0.50Discolored-paralog2oxidoreductase activity198,119,074–198,130,785
Zm00001d0319920.93−0.49indeterminate domain p1*208,942,917–208,950,782
“*” represent that the genes have no GO annotation in the Gramene website.
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

Du, L.; Peng, X.; Zhang, H.; Xin, W.; Ma, K.; Liu, Y.; Hu, G. Transcriptome Analysis and QTL Mapping Identify Candidate Genes and Regulatory Mechanisms Related to Low-Temperature Germination Ability in Maize. Genes 2023, 14, 1917. https://doi.org/10.3390/genes14101917

AMA Style

Du L, Peng X, Zhang H, Xin W, Ma K, Liu Y, Hu G. Transcriptome Analysis and QTL Mapping Identify Candidate Genes and Regulatory Mechanisms Related to Low-Temperature Germination Ability in Maize. Genes. 2023; 14(10):1917. https://doi.org/10.3390/genes14101917

Chicago/Turabian Style

Du, Lei, Xin Peng, Hao Zhang, Wangsen Xin, Kejun Ma, Yongzhong Liu, and Guangcan Hu. 2023. "Transcriptome Analysis and QTL Mapping Identify Candidate Genes and Regulatory Mechanisms Related to Low-Temperature Germination Ability in Maize" Genes 14, no. 10: 1917. https://doi.org/10.3390/genes14101917

APA Style

Du, L., Peng, X., Zhang, H., Xin, W., Ma, K., Liu, Y., & Hu, G. (2023). Transcriptome Analysis and QTL Mapping Identify Candidate Genes and Regulatory Mechanisms Related to Low-Temperature Germination Ability in Maize. Genes, 14(10), 1917. https://doi.org/10.3390/genes14101917

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop