Next Article in Journal
Identification of Myocardial Insulin Resistance by Using Liver Tests: A Simple Approach for Clinical Practice
Next Article in Special Issue
ATP-Binding Cassette G Transporters and Their Multiple Roles Especially for Male Fertility in Arabidopsis, Rice and Maize
Previous Article in Journal
Lidocaine Inhibited Tendon Cell Proliferation and Extracellular Matrix Production by Down Regulation of Cyclin A, CDK2, Type I and Type III Collagen Expression
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
IJMS
Volume 23
Issue 15
10.3390/ijms23158786
ijms-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

Combined QTL Mapping across Multiple Environments and Co-Expression Network Analysis Identified Key Genes for Embryogenic Callus Induction from Immature Maize Embryos

by
Yun Long
1,2,
Tianhu Liang
1,
Langlang Ma
1,
Peng Liu
1,
Yun Yang
3,
Xiaoling Zhang
1,
Chaoying Zou
1,
Minyan Zhang
1,
Fei Ge
1,
Guangsheng Yuan
1,
Thomas Lübberstedt
4,
Guangtang Pan
1 and
Yaou Shen
1,*
1
State Key Laboratory of Crop Gene Exploration and Utilization in Southwest China, Maize Research Institute, Sichuan Agricultural University, Chengdu 611130, China
2
Key Laboratory of Southwest China Wildlife Resources Conservation (Ministry of Education), College of Life Science, China West Normal University, Nanchong 637002, China
3
Nanchong Academy of Agricultural Sciences, Nanchong 637000, China
4
Department of Agronomy, Iowa State University, Ames, IA 50010, USA
*
Author to whom correspondence should be addressed.
Int. J. Mol. Sci. 2022, 23(15), 8786; https://doi.org/10.3390/ijms23158786
Submission received: 5 July 2022 / Revised: 5 August 2022 / Accepted: 5 August 2022 / Published: 7 August 2022
(This article belongs to the Special Issue Advances and Innovations in Maize Genetics and Breeding Research)
Browse Figures
Versions Notes

Abstract

:
The ability of immature embryos to induce embryogenic callus (EC) is crucial for genetic transformation in maize, which is highly genotype-dependent. To dissect the genetic basis of maize EC induction, we conducted QTL mapping for four EC induction-related traits, the rate of embryogenic callus induction (REC), rate of shoot formation (RSF), length of shoot (LS), and diameter of callus (DC) under three environments by using an IBM Syn10 DH population derived from a cross of B73 and Mo17. These EC induction traits showed high broad-sense heritability (>80%), and significantly negative correlations were observed between REC and each of the other traits across multiple environments. A total of 41 QTLs for EC induction were identified, among which 13, 12, 10, and 6 QTLs were responsible for DC, RSF, LS, and REC, respectively. Among them, three major QTLs accounted for >10% of the phenotypic variation, including qLS1-1 (11.54%), qLS1-3 (10.68%), and qREC4-1 (11.45%). Based on the expression data of the 215 candidate genes located in these QTL intervals, we performed a weighted gene co-expression network analysis (WGCNA). A combined use of KEGG pathway enrichment and eigengene-based connectivity (KME) values identified the EC induction-associated module and four hub genes (Zm00001d028477, Zm00001d047896, Zm00001d034388, and Zm00001d022542). Gene-based association analyses validated that the variations in Zm00001d028477 and Zm00001d034388, which were involved in tryptophan biosynthesis and metabolism, respectively, significantly affected EC induction ability among different inbred lines. Our study brings novel insights into the genetic and molecular mechanisms of EC induction and helps to promote marker-assisted selection of high-REC varieties in maize.

1. Introduction

Maize (Zea mays L.) is an important food and forage crop worldwide. Breeding superior varieties is the guarantee of high and stable yield in maize. Traditional genetic modification methods have a long cycle and are limited by natural variations [1], whereas transgenic breeding is able to break through these restrictions. In maize genetic transformation, the induction efficiency of embryogenic callus (EC) directly affects the transformation ratio of immature embryos [2]. Previous studies reported that EC induction from maize immature embryos was influenced by multiple factors, including genotypes, hormones, and compositions of induction media [3,4,5]. Notably, the capacity of EC formation showed serious genotype dependence when explants were incubated under the same condition [6]. Most elite lines in maize have low ratios of embryogenic callus induction (REC), which inhibits the development of maize transgenic breeding [7]. Identification of causal genes that control maize EC induction will contribute to the improvement of genetic transformation efficiency and transgenic breeding in maize.
Genetic architecture for EC induction and regeneration abilities have been studied in several crops. For instance, QTLs controlling callus induction and another culture ability have been detected in wheat, while QTLs for callus induction rate and plant regeneration ability have been identified in rice [8,9,10]. In maize, a high-REC maize inbred line, 18–599 (R), and a low-REC line, R15, were used to construct an F2 population [11]. Using this population, five QTLs controlling callus induction efficiency were mapped on chromosomes 1, 3, 7, and 8, respectively, with the explained phenotypic variation ranging from 5.25% to 23.4% [11]. Eight QTLs responsible for tissue culture capacities were determined by using 239 recombinant inbred lines (RILs) derived from the maize inbred lines Huangzao4 (with a lower REC) and Mo17 (with a higher REC) [12]. These QTLs were located on chromosomes 2, 3, 5, 6, 8, and 9, which explained 4.78–14.02% of phenotypic variation [12]. A B73 near-isogenic line with a high REC was used for fine mapping of the QTLs responsible for EC response, and finally, a QTL was located in a 3053 kb region in bin 3.06 of chromosome 3 [13]. Based on the physical positions, these reported EC induction-related QTLs were located in different genomic regions of maize. Moreover, owing to the large intervals of these identified QTLs, no causal genes were identified based on the above QTLs. Conventionally, to further identify the functional genes involved in EC induction, QTL fine mapping should be conducted based on enlarged segregation populations and enhanced marker resolutions [14].
In addition to QTLs, several genes were reported to influence EC formation in maize, including LEAFY COTYLEDON2 (LEC2), BABY BOOM (BBM), and WUSCHEL2 (WUS2). The expression of LEC2 was significantly upregulated during somatic embryogenesis in Arabidopsis explants induced in vitro, and a close link was observed between auxin and the LEC2 activity [15]. The overexpression of maize BBM and WUS2 genes produced higher transformation frequencies in many previously reported non-transformable maize inbred lines [16,17]. In a genome-wide association study of maize callus regeneration, the homologous gene of WUS2, WOX2, was confirmed as the candidate gene regulating the capacity of EC regeneration [18]. Our previous studies identified ZmMYB138 and ZmSAUR15, which modulated EC formation by participating in the hormone signal transduction pathway [2,19]. ZmSAUR15 encodes a small auxin-upregulated RNA, which negatively mediates the capability of maize EC formation by affecting IAA biosynthesis and cell division in immature embryo-derived callus [2]. ZmMYB138 promotes the formation of EC through GA signal transduction in maize immature embryos [19].
Recently, a weighted gene co-expression network analysis (WGCNA) was developed to identify the hub genes in gene co-expression modules according to the gene expression patterns derived from transcriptome data [20]. WGCNA has been proven to facilitate the excavation of causal genes in larger QTL intervals of target traits, as it circumvents the requirement of further QTL fine mapping [14]. In this study, an IBM Syn10 DH population [21], which was constructed from the cross of maize inbred lines B73 and Mo17 and has a high recombination ratio, was used to identify the QTLs controlling EC induction from immature embryos of maize. All the gene models located in these QTL intervals were used to perform a WGCNA based on the RNA-seq data of four maize lines during the EC induction process. The EC induction-associated modules and hub genes were identified according to the KEGG pathway enrichment and eigengene-based connectivity (KME) values. Finally, gene-based association studies were performed to detect the intragenic variations influencing EC induction. Our study provided insights into understanding genetic and molecular mechanisms underlying EC formation in maize and developing molecular markers for improving the REC of immature maize embryos.

2. Results

2.1. Phenotypic Performances of EC Induction Traits under Three Environments

To evaluate the phenotypic performances of EC induction capability among the IBM Syn10 DH population, we investigated the REC, DC, RSF, and LS for this population and its parents, B73 and Mo17, across three environments. The EC had a loose and fragile structure and formed globular nodules on the callus surface (Figure 1a,c), whereas non-embryogenic callus (NEC) presents a soft, transparent, and waterlogged structure (Figure 1b,d). After 30 d of incubation, the immature embryos of the parent B73 only formed NEC, with the REC = 0; approximately 5.79% (average across the three environments) of the Mo17 embryos formed EC (REC = 5.79%) (Table 1). However, the mean values of the rate of shoot formation (RSF), length of shoot (LS), and diameter of callus (DC) across the three environments were significantly greater in B73 than those in Mo17.
Among the IBM Syn10 DH population, the REC ranged from 0 to 52.17%, 0 to 58.61%, and 0 to 56.48% under the three environments of Xishuangbanna (XSBN), Ya’an (YA), and Chongzhou (CZ), respectively, with the coefficients of variation (CVs) being 2.00, 2.90, and 2.67. However, most of these recombinant lines cannot form EC (with the REC = 0). The skewness and kurtosis of REC were > 1 under each environment (Table 1), suggesting that REC was an atypical quantitative trait (Figure 2). Interestingly, five recombinant lines displayed an REC > 10% (56.48%, 42.88%, 44.91%, 58.61%, and 12.04%) (Table 1 and Figure 1), exceeding their two parents, B73 and Mo17, in REC. It indicated that higher REC-associated alleles from B73 and Mo17 have been integrated into these recombinant lines. Meanwhile, the RSF, LS, and DC were continuously distributed among this population, with the respective CVs ranging from 0.35 to 0.37, 0.42 to 0.48, and 0.16 to 0.23 across different environments (Table 1). The absolute skewness and kurtosis values of these three traits were all smaller than 1 except for the RSF in the CZ environment, verifying that the RSF, LS, and DC generally conformed to normal distributions and were typical quantitative traits (Table 1 and Figure 2).
Significantly negative correlations were observed between REC and RSF, REC and LS, and REC and DC in the three environments except for REC and RSF in CZ (Table 2). The correlation coefficients (CCs) of the three pairs of traits ranged from −0.105 to −0.342, −0.134 to −0.263, and −0.177 to −0.318 across the three environments, respectively. Furthermore, higher phenotypic correlations between DC and RSF, DC and LS, and LS and RSF were found across different environments, and the CCs ranged from 0.482 to 0.571, 0.612 to 0.715, and 0.692 to 0.737, respectively (p < 0.01). The broad-sense heritability (H2) of the four EC induction-related traits varied between 83.17 and 88.27% (Table 3), indicating that genetic factors played dominant roles in the formation of EC. Analysis of variance (ANOVA) revealed that the variances from genotype (G), environment (E), and genotype × environmental interactions (G × E) were all significant for the four traits (p < 0.01, Table 3).

2.2. QTL Responsible for EC Induction

To detect the QTLs responsible for maize EC induction, we conducted linkage mapping by using the phenotype values collected from each environment and best linear unbiased prediction (BLUP) values. Finally, a total of 41 QTLs were identified for the four EC induction-related traits under the three environments and by using the BLUP values. These QTLs were distributed on 10 chromosomes, and the phenotypic variation explained (PVE) by a single QTL ranged from 4.01 (qLS10-1 using BLUP) to 11.54% (qLS1-1 in XSBN) (Table 4 and Figure 3). Three QTLs had the PVE > 10% and were thus considered as the major QTLs for EC induction ability. Furthermore, seven QTLs were individually simultaneously identified across different environments or across environment(s) and BLUP values, including three for DC, one for RSF, and three for LS. Among the 41 QTLs, 18 and 23 exhibited positively and negatively additive effects, respectively. The detailed information for these QTLs is presented below.

2.2.1. DC

Thirteen QTLs controlling DC were mapped to chromosomes 1, 3, 4, 7, 9, and 10, respectively. These individually explained 4.61–9.98% of the phenotypic variation and accounted for 28.15%, 30.32%, 22.32%, and 28.19% of the total DC variation under XSBN, YA, CZ, and using the BLUP values, respectively (Table 4). Among them, 11 QTLs presented positive additive effects, indicating that alleles from B73 contributed to the increase in DC (Table 4). Notably, qDC1-2 and qDC9-3 were separately co-detected across CZ and the BLUP model, whereas qDC7-1 was repeatedly identified under YA and by using the BLUP values (Table 4 and Figure 3).

2.2.2. RSF

A total of 12 QTLs located on chromosomes 1, 2, 3, 4, 6, 8, 9, and 10 were identified for RSF (Table 4 and Figure 3). The PVE of a single QTL ranged from 5.33 to 9.97%, and the logarithm of odds (LOD) scores of these QTLs varied between 2.68 and 5.71 (Table 4). The total phenotypic variation explained by these RSF QTLs was 19.50%, 29.48%, 29.36%, and 14.03%, in XSBN, YA, CZ, and using BLUP, respectively (Table 4). One QTL (qRSF9-1) was co-identified across XSBN and BLUP, whereas the others were environment-specific QTLs (Table 4). In addition, eight QTLs on chromosomes 1, 2, 4, 8, and 10 showed negative effects, implying that alleles from Mo17 contributed to the increased RSF (Table 4).

2.2.3. LS

Ten QTLs affecting LS were detected on chromosomes 1, 2, 3, 9, and 10 across multiple environments and using BLUP. These QTLs explained 4.01–11.54% of the phenotypic variation, among which two major QTLs (qLS1-1 and qLS1-3) had a PVE > 10%. Moreover, three QTLs were separately repeatedly detected under multiple environments or across environment(s) and BLUP (Table 4 and Figure 3). Two common QTLs (qLS1-1 and qLS3-2) were individually repeatedly detected in XSBN and using BLUP, which explained 74–11.54% and 7.34–9.20% of phenotypic variation across different environments, respectively (Table 4). The other common QTL (qLS10-1) was simultaneously identified across XSBN, YA, and BLUP, with the PVE ranging from 4.01 to 6.98%. Moreover, nine QTLs showed negative additive effects, suggesting that the increased LS was mainly contributed by the alleles from Mo17 (Table 4).

2.2.4. REC

Across the three environments and BLUP values, six QTLs for REC were mapped to chromosomes 1, 4, 5, 6, and 7, contributing 4.01–11.45% to the phenotypic variation (Table 4 and Figure 3). No QTL was repeatedly detected across multiple environments or across environment(s) and BLUP, suggesting that these identified QTLs were environment-sensitive QTLs (Table 4 and Figure 3). Notably, a major QTL (qREC4-1) was detected in XSBN, which explained 11.45% of the phenotypic variation (Table 4). In addition, among the six QTLs controlling REC, four QTLs on chromosomes 5, 6, and 7 displayed negative addition effects, indicating that the increased REC were mainly contributed by the alleles from Mo17 (Table 4).

2.3. QTL Clusters for Embryogenic Callus Induction-Related Traits

By comparing the physical positions among different QTLs identified in the present study, 12 were classified into four QTL clusters, which were located on chromosomes 1, 9 and 10, respectively (Table 5). Cluster A contained the QTLs controlling DC, LS, and RSF. Clusters B and D involved the QTLs for RSF and LS, whereas Cluster C included RSF and DC QTLs. Notably, each of the clusters contained the QTLs responsible for RSF, probably due to the high phenotypic correlation between RSF and each of DC and LS (Table 2 and Table 5). These universally observed clusters among DC, LS, and RSF reflected the pleiotropic effects of these QTLs on DC, LS, and RSF during the process of EC induction in immature embryos. However, the QTLs affecting REC were not included in any of clusters, which was consistent with the observation of the negative correlation between REC and each of DC, LS, and RSF. Collectively, these results suggested that the casual genes responsible for REC were independent from those for DC, LS, and RSF (Table 2 and Table 5).

2.4. Validation of QTL Intervals

Reliability of one of the detected QTLs (qREC6-1) was validated by evaluating co-segregation between genetic variations in the QTL interval and REC phenotype using six high-REC lines (Mo17 and five other lines with REC >10%) and six low-REC lines (REC approximately 0). A total of 10 fragments covering 1035–3300 nucleotides in qREC6-1 were amplified by PCR among these 12 lines. In fragment 1 (F1), only one SNP (G/A) was found; all six high-REC lines contained the G-allele, and the six low-REC ones contained the A-allele (Figure S1). In F2, three SNPs (G/A, T/C, T/C) were identified, and they formed two haplotypes (GTT and ACC). The six high-REC lines had the haplotype GTT, whereas the six low-REC lines had the haplotype ACC (Figure S2). Moreover, we identified four SNPs in F3, one SNP and one InDel in F4, four SNPs in F6, and three SNPs in F7. These 12 lines displayed co-segregations between the REC phenotype and genotype of each fragment (Figures S3, S4, S6, and S7). For each of F5, F8, F9, and F10, the high-REC lines possessed the same allele/haplotype, whereas the low-REC lines had the other allele/haplotype, except for the line IBM182 (a high-REC line), which contained the same alleles as the six low-REC lines (Figures S5 and S8–S10). Combined, these results verified the high credibility of the QTLs detected in the present study. Moreover, these SNPs/InDels can be used to develop the molecular markers for improvement of REC in immature maize embryos.

2.5. Candidate Genes and Co-Expression Networks

To construct the co-expression networks involved in EC induction, we carried out a WGCNA for candidate genes responsible for EC induction capability. A total of 215 gene models were located in the confidential intervals of detected QTLs, based on the B73 v4 genome. Among them, 74, 58, 87, and 14 genes were responsible for DC, RSF, LS, and REC, respectively (Table S1). Our previous study performed transcriptome sequencing for four maize lines (CN9802, ZM28, JS0251, and YA3237) during the EC induction process. According to the transcriptome data, 122 of these 215 genes were differentially expressed at 5 d, 10 d, and/or 15 d after EC induction culture in at least one line, with the |log2(fold change)| > 1, and Q < 0.05. We then conducted a WGCNA using the expression data of the 122 genes. Finally, these 122 genes were grouped into six co-expression modules, namely green (14 genes), turquoise (32 genes), brown (24 genes), yellow (15 genes), blue (26 genes), and red (11 genes) (Table S2). Remarkably, KEGG enrichment analysis showed that the genes in the blue module were enriched in the phenylalanine, tyrosine and tryptophan biosynthesis pathway (Zm00001d028477 and Zm00001d047896), tryptophan metabolism pathway (Zm00001d034388), and plant hormone signal transduction pathway (Zm00001d022542) (Figure 4a,b). Since plant callus induction was extensively reported to correlate with tryptophan biosynthesis and metabolism, as well as hormone signal transduction [22,23,24], the blue module was thus considered as the EC formation-associated module in the present study. The above four genes were accordingly taken as the hub genes affecting EC induction, which showed high KME values (0.87–0.92) in the blue module (Figure 4c and Table S2). Generally, the four hub genes displayed increasing expression at each EC induction stage when compared with 0 h, verifying their response to EC induction (Figure 4d).

2.6. Hub Gene-Based Association Mapping

To further investigate whether these variations within the hub genes affected the REC among different maize lines, we performed hub gene-based association mapping by combining the phenotypes and genotypes of 75 lines with distinct genetic backgrounds. The maize lines consisted of the Tropical, Stiff Stalk (SS), and non-Stiff Stalk (NSS) germplasms previously genotyped using a maize 56K SNP Array (Zhang et al., 2016). Cluster analysis based on the genotypes showed that these 75 lines were classified into six different subgroups with 1–19 lines in each subgroup (Figure S11). The DC, RSF, LS, and REC of the lines were collected under three environments in our previous study [2]. The genotypes of these lines were obtained using PCR amplification and SNP/InDel calling for the four genes in the present study. After removing variants with missing rate > 0.2 or heterozygous proportion > 0.3, we identified 42 variants (42 SNPs) in Zm00001d028477, 10 variants (10 SNPs) in Zm00001d047896, 50 variants (48 SNPs and 2 InDels) in Zm00001d034388, and 44 variants (41 SNPs and 3 InDels) in Zm00001d022542 (Table S3). Association mapping revealed that six SNPs were significantly (p < 0.05) associated with the REC, DC, RSF, and LS in Zm00001d028477, and 20 SNPs were significantly (p < 0.05) associated with the four traits in Zm00001d034388 (Table S4). Meanwhile, two separate SNPs, one in Zm00001d022542 and the other in Zm00001d047896, were significantly (p < 0.05) associated with LS (Table S4). Among the six significant SNPs in Zm00001d028477, three SNPs (LS: S1_36255280 and S1_36256420; RSF: S1_36256406) were intron-splicing or missense variants, whereas the other three SNPs (REC: S1_36247233 and S1_36249669; DC: S1_36247066) were intron or synonymous variants (Table S4). For Zm00001d034388, nine SNPs (REC: S1_291663139; DC: S1_291664657; LS: S1_291666764 and S1_291663835; RSF: S1_291663009, S1_291665448, S1_291665505, S1_291664082, S1_291664657, and S1_291665851) were missense, 3′-UTR, or intron-splicing variants, whereas the other eleven SNPs were synonymous or intron variants (Table S4). Two SNPs located in Zm00001d022542 (LS: S7_179752120) and Zm00001d047896 (LS: S9_144619659) were also intron variants (Table S4).
It is widely accepted that the missense and intron-splicing variants and those located in the UTR and promoter regions tend to generate alterations in amino acid or mRNA abundance [25,26]. We thus divided these 75 lines into different haplotypes by combining the types of significant variants (called effective variants in the present study) in each gene for a single associated trait. For Zm00001d034388, the effective REC-associated SNP classified these lines into two groups, and the group with the C-allele showed a significantly (p < 0.05) lower REC (5.50%) than that (22.60%) of the group with the T-allele (Figure 5a,b). Meanwhile, the effective DC-associated SNP divided these lines into two groups, and the group with the A-allele displayed a significantly (p < 0.05) higher DC (0.35 cm) relative to that (0.21 cm) of the group with the G-allele (Figure 5a,c). Furthermore, the six effective RSF-associated SNPs categorized these lines into three major haplotypes (Hap1: CGAACA; Hap2: CGGGGC; Hap3: TGGGGC), among which Hap2 had a significantly (p < 0.05) higher RSF (87.72%) than that (74.90%) in Hap3 (Figure 5a,d). According to the two effective LS-associated variants, these lines were divided into two major haplotypes (Hap1: CT; Hap2: TT), and Hap1 displayed a significantly (p < 0.001) larger LS value (1.04 cm) than that (0.59 cm) of Hap2 (Figure 5a,e). As for Zm00001d028477, the two effective LS-associated SNPs formed two major haplotypes (Hap1: AG; Hap2: AT), and the LS (0.73 cm) of Hap1 lines was significantly (p < 0.05) higher than that (0.47 cm) of Hap2 lines (Figure 5f,g). The effective RSF-associated SNP divided these 75 lines into two genotypes (A and G), and the RSF (51.23%) of the lines containing the A-allele was significantly (p < 0.01) lower than that (85.99%) of those containing the G-allele (Figure 5f,h). Collectively, these findings illustrated that the variants in Zm00001d028477 and Zm00001d034388 influenced the formation of EC among different maize lines.

3. Discussion

3.1. Use of IBM Syn10 DH Population for Mapping Embryogenic Callus Induction-Related Traits

Abundant phenotypic variations and high heritability estimates are both important to perform QTL mapping of complex traits [27]. In the present study, the parent B73 cannot form EC, whereas the REC of parent Mo17 ranged from 3.39–7.05% across three environments (Table 1). Notably, the REC of the IBM Syn10 DH population showed large variation, with the phenotypic value ranging from 0 to 52.17%, 0 to 58.61%, and 0 to 56.48% in XSBN, YA, and CZ, respectively (Table 1 and Figure 2). The broad-sense heritability estimates were larger than 80% for the investigated four traits, which contributed to the identification of the environmentally stable QTLs. A total of forty-one EC induction-related QTLs were detected, among which seven common QTLs were individually simultaneously detected across multiple environments or across the environment(s) and BLUP values, including three for DC (qDC1-2, qDC7-1, and qDC9-3), one for RSF (qRSF9-1), and three for LS (qLS1-1, qLS3-2, and qLS10-1). Moreover, high recombination rates and high-density markers in segregation populations facilitate fine mapping of the genetic loci controlling target traits [27]. The IBM Syn 10 DH population was obtained via six additional generations of open pollination on the IBM Syn4 RIL population and genotyped using Next-Generation Sequencing, resulting in a higher ratio of genetic recombination and a smaller genetic distance between two adjacent markers [18]. In the present study, the majority of the QTLs were mapped within < 0.2 Mb physical intervals, such as all the REC QTLs (qREC1-1, qREC4-1, qREC5-1, qREC6-1, qREC7-1, and qREC7-2). Compared to the larger intervals reported in previous studies on QTL mapping of EC induction [9,12,13], the smaller QTL intervals detected in our study facilitate further excavation of causal genes responsible for maize EC induction.

3.2. Maize EC Induction Is Probably Controlled by a Few Major Genes

In contrast to LS, RSF, and DC, REC did not conform to a normal distribution among the IBM Syn 10 DH population (Figure 2), which was consistent with our observations in an association panel in our previous study [2]. It is widely accepted that quantitative traits show normal distributions only under the polygenic hypothesis [27]. However, the phenotype values will display a skewness among populations when the number of QTLs is small and there are a few QTLs with large genetic effects [28,29]. In the present study, only six QTLs were detected across different environments and using BLUP, including one major QTL with a PVE > 10% and five minor QTLs (Table 4 and Figure 3). Presumably, EC induction in immature maize embryos is controlled by a few major genes, which is helpful in the exploration of causal genes responsible for EC induction.

3.3. QTLs for Embryogenic Callus Induction-Related Traits

To date, most research on maize EC has focused on REC, callus weight, callus color, callus texture, callus proliferation ability, callus browning tendency, and callus regeneration ability [13,30,31,32]. However, adventitious shoots growing on immature embryos are considered to negatively influence the formation of EC in embryos [33]. Consistently, LS and RSF were significantly negatively associated with the REC across different environments in the present study (Table 2). Meanwhile, DC reflects the growth rate of EC [34]. Therefore, RSF, LS, and DC were taken as the novel traits for detecting genetic architecture of maize EC formation in the present study.
The results of QTL mapping are susceptible to environmental influences [35], and QTL loci repeatedly identified in multiple environments are considered more convincing [36]. In this study, seven common QTLs (DC: three; RSF: one; LS: three) were separately simultaneously detected across multiple environments or across the environment(s) and BLUP (Table 4 and Figure 3). The QTLs controlling REC were separately mapped on chromosomes 1, 4, 5, 6, and 7, whereas none was repeatedly detected in multiple environments in this study (Table 4 and Figure 3). This implied a high environmental effect on REC in maize. A previous study reported two QTLs associated with the formation of embryo-like structures, which were located on the chromosomes 2 and 8, respectively [37]. The regions on chromosomes 1, 2, and 3 harbored the QTL controlling somatic embryogenesis in maize [38]. Five QTLs on chromosomes 1, 3, 7, and 8 were found to affect efficiency callus induction [11], whereas one QTL on chromosome 5 conferred EC initiation in maize [12]. Furthermore, a 23.9 Mb region on chromosome 3 was found to regulate embryogenic and regenerable tissue culture response in maize [13]. We compared the physical positions between the REC QTL identified in the present study and those reported previously; however, no overlap was found. One possible reason is that the mapping populations are different between our study and the previous studies, resulting in distinct variation effects on REC.

3.4. Candidate Genes Involved in Embryogenic Callus Induction

Up to now, only a few genes have been proven to affect EC induction in maize, such as ZmSAUR15, WUS2, BBM, LEC, SERK, and ZmMYB138 [15,16,17,19,39]. However, none of the above causal genes were located in the QTLs identified in the present study. The possible explanation was that no EC induction-associated variants existed within the genes among this IBM Syn10 DH population. In the present study, a combination of QTL mapping and WGCNA suggested that Zm00001d028477, Zm00001d047896, Zm00001d034388, and Zm00001d022542 were the hub genes in the EC induction-associated module. Among them, Zm00001d028477 and Zm00001d047896 encode chorismite synthases and participate in tryptophan biosynthesis. Zm00001d034388 was annotated as aldehyde oxidase 4, which is involved in tryptophan metabolism. Tryptophan acts as the precursor in auxin biosynthesis of plants [40], and the biosynthesis and metabolism of tryptophan exert impacts on auxin signal transduction [41]. A number of studies demonstrated the significance of auxin signal transduction on EC induction [42,43,44]. Our previous study revealed that the knockdown and overexpression of a key gene (ZmSAUR15) involved in auxin signal transduction improved and inhibited the induction of maize EC [2], respectively. Zm00001d022542 encodes a TGA transcription factor, a key factor in the pathway of salicylic signal transduction. Salicylic acid was extensively reported to affect callus culture and regeneration in different plant species [45,46,47]. The four hub genes showed increasing expression during the process of EC induction in the four lines, with the exception of Zm00001d047896 in the line CN9802. Gene-based associations revealed that significant variations within Zm00001d028477 and Zm00001d034388 affected EC induction-related traits among different inbred lines. These significant loci can be used to develop functional markers for improving the REC of immature maize embryos.

4. Materials and Methods

4.1. Plant Materials and Field Trials

The IBM Syn10 DH population used for QTL mapping was derived from a cross between the lines B73 and Mo17, as described in our previous study [21]. In brief, this population was constructed through six additional generations of intercrossing on the IBM Syn4 population followed by haploid induction and chromosome doubling. Therefore, the IBM Syn10 DH population has a high genetic recombination frequency [22]. A total of 210 doubled haploid (DH) lines and their parents (B73 and Mo17) were investigated for the phenotypes of REC, RSF, LS, and DC. These lines were grown under three environments, namely XSBN (Yunnan province; 22.00° N,100.79° E; elevation, 553 m) in 2015, YA (Sichuan province; 29.59° N, 102.57° E; elevation, 516 m) in 2015, and CZ (Sichuan province, 30.30° N, 103.07° E; elevation, 560 m) in 2016. XSBN has a typical monsoon climate with a mean annual temperature of 21.3 °C in 2015, whereas YA and CZ have a monsoon-influenced humid subtropical climate with mean annual temperatures of 17.3 °C (YA, 2015) and 16.8 °C (CZ, 2016). A completely randomized block design was carried out with three replicates for each environment. Each line was grown in three rows with 14 plants per row. The row length was 3 m, and the space between two adjacent rows was 0.75 m. In each location, this population was managed using a standard cultivation practice for corn.

4.2. Immature Embryo Culture and Callus Induction

Three immature ears were harvested from the middle row per line at 12 d after self-pollination for each replicate. For each line, 108 immature embryos with 1.2–1.5 mm in length were collected from the middle sections of the three immature ears and evenly incubated among three Petri dishes with a modified N6 inducting medium [2]. These embryos were aseptically cultured at 28 °C in the dark for 30 days, with the scutella facing upward. A total of nine ears from three field repetitions were accordingly classified into three replicates in tissue culture.

4.3. Phenotype Investigation

The RSF, which was calculated as (number of embryos forming shoots/total number of inoculated embryos) × 100%, was investigated after 7 d of incubation. The LS (mm) and DC (mm) were measured after 10 d of incubation. The REC, which was equal to (number of embryos forming embryogenic calli/total number of inoculated embryos) × 100%, was scored after 30 d of incubation. The criteria for identification of EC were in accordance with the previous description [48].

4.4. Phenotypic Data Analysis

The average value of each trait among three replicates was taken as the phenotype value for each line. The SPSS software (version 26.0, IBM Corp., Armonk, NY, USA) was utilized to analyze the phenotypic data, including the range, mean value, CV, skewness, kurtosis, and correlation of each trait across three environments. Variance analysis incorporating genotype, environment, and genotype–environment interaction was conducted with the general linear model (GLM) in the SPSS 26.0 software. The broad-sense heritability (H2) of each trait was calculated following the previous description [49]. To reduce environmental effects on genotype, the BLUP values for the four traits were calculated using the lme4 package in the R software. The model is as follows: Phenotype~ (1|genotype) + (1|repeat % in % environment) + (1|genotypeand environment).

4.5. QTL Analysis

A bin map with high-density markers was constructed for the IBM Syn10 DH population in our previous study [50]. This bin map was approximately 11,198.5 cM in length, which contained 6618 bin markers with the average genetic distance between bin markers being 1.7 cM. The QTL IciMapping software version 3.0 [51] was applied to detect QTLs under three environments and using the BLUP values, based on inclusive composite interval mapping (ICIM). The testing window size of scan configuration was set to 10 cM, with a 2 cM walk speed. The significance threshold of LOD for QTL detection was set to 2.5, as described in the previous study [52]. Moreover, the QTLs explaining phenotypic variation > 10% were considered as major QTLs regulating EC induction. For each trait, QTLs with overlapped confidence intervals were assumed to be the same QTL [51]. The rule of QTL naming was as follows: q + trait + serial number of chromosome-serial number of the identified QTL. For example, in “qREC6-1”, “q” stands for QTL, “REC” is the abbreviation of induction rate of embryogenic callus, “6” represents chromosome 6, and “1” means the first QTL identified on chromosome 6.

4.6. DNA Extraction and Variation Validation

The genomic DNA was extracted from the young leaves of B73, Mo17, 5 high-REC IBM Syn10 DH lines (IBM003, IBM097, IBM182, IBM270, and IBM304), and 5 low-REC IBM Syn10 DH lines (IBM009, IBM062, IBM144, IBM234, and IBM327) using the CTAB method [53]. In the qREC6-1 interval, 10 fragments with 1035–3300 nt in length were randomly selected and amplified using PCR. The amplified fragments were sequenced with Sanger sequencing. The sequencing results were aligned and analyzed with the SnapGene software (version 2.3.2, Insightful Science, USA, https://www.snapgene.com/).

4.7. WGCNA

Expression data of the gene models located in the detected QTL intervals were obtained from our RNA-Seq data of four maize lines (CN9802, ZM28, JS0251, and YA3237) during EC induction (Accession number: CRA007381, https://ngdc.cncb.ac.cn/gsa, accessed on 1 July 2022). Briefly, transcriptomes of the embryos were analyzed at 0 d, 5 d, 10 d, and 15 d after callus induction culture in the four maize lines. The clean reads were obtained by filtering the raw sequencing reads using the fastp software (version 0.23.1, Shifu Chen, Shenzhen, China) [54]. Gene expression values were calculated based on the clean reads mapped to the B73 RefGen_V4 genome and were subsequently normalized to transcripts per kilobase million (TPM) using an in-house script. Among the candidate genes obtained by QTL mapping, the differentially expressed genes (DEGs) were identified according to the absolute value of log2(fold change) between each induction stage (5, 10, and 15 d) and 0 d, and the threshold of DEGs was set to 1. The TPM of the DEGs were then subjected to a WGCNA using the WGCNA package in R software [20]. To detect the co-expression modules, WGCNA program parameters were set as follows: expression of variance data > 0; soft threshold = 11; max block size = 122; deep split = 2; min module size = 10; merge cut height = 0.2. In each of the co-expression modules, the Kyoto encyclopedia of genes and genomes (KEGG) pathway was analyzed on the GENE DENOVO platform (https://www.omicshare.com/tools/Home/Soft/pathwaygseasenior, accessed on 8 may 2021, GENE DENOVO Biotechnology Co., Ltd., Guangzhou, China). The co-expression network of the hub genes was drawn using the Cytoscape software (verision 3.9.1, Paul Shannon, Seattle, USA) [55].

4.8. Association Analysis

For each hub gene, the sequences of the gene body and its 2000 bp upstream were amplified by PCR in 80 randomly selected lines from our previously reported association panel [2]. PCR primers were designed in Primer 3 (version 0.4.0, Triinu Koressaar, Tartu, Estonia) based on the gene sequences in B73 RefGen_V4. Sequence alignment between the PCR-amplified sequences and the B73 genome sequence was conducted using the DNAMAN software (version 5.2.2, Lynnon Corp., Quebec, Canada) to detect sequence polymorphisms including single nucleotide polymorphisms (SNPs) and insertion/deletions (InDels). A GLM + PCA model in TASSEL software (version 4.0, Peter J. Bradbury, New York, USA) was used to detect associations between the phenotype values and polymorphism loci, with the significance threshold set as p = 0.05 [2]. Linkage disequilibrium (LD) decay between the markers was evaluated using HaploView software (Jeffrey Barrett, Cambridge, MA, USA). Haplotype identification was performed using the significant markers located in each hub gene.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/ijms23158786/s1.

Author Contributions

Conceptualization, Y.S.; Data curation, Y.L., L.M. and P.L.; Formal analysis, T.L. (Thomas Lübberstedt); Funding acquisition, Y.S., L.M. and F.G.; Investigation, Y.L., X.Z., C.Z., M.Z., F.G. and G.Y.; Methodology, Y.S.; Project administration, Y.S.; Resources, T.L. (Tianhu Liang); Software, Y.L., L.M. and P.L.; Supervision, G.P.; Validation, Y.L. and Y.Y.; Visualization, Y.L.; Writing—original draft, Y.L.; Writing—review and editing, Y.S. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Foundation of China (32072073, 31871637, and 32001500) and Sichuan Science and Technology Program (2021JDTD0004 and 2021YJ0476).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Transcriptome data used in this study are available in the Genome Sequence Archive (GSA) in National Genomics Data Center (NGDC) database with the accession number CRA007381.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Ahmar, S.; Gill, R.A.; Jung, K.-H.; Faheem, A.; Qasim, M.U.; Mubeen, M.; Zhou, W. Conventional and molecular techniques from simple breeding to speed breeding in crop plants: Recent advances and future outlook. Int. J. Mol. Sci. 2020, 21, 2590. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Wang, Y.; He, S.; Long, Y.; Zhang, X.; Zhang, X.; Hu, H.; Li, Z.; Hou, F.; Ge, F.; Gao, S. Genetic variations in ZmSAUR15 contribute to the formation of immature embryo-derived embryonic calluses in maize. Plant J. 2022, 109, 980–991. [Google Scholar] [CrossRef]
  3. Satish, L.; Rathinapriya, P.; Ceasar, S.A.; Rency, A.S.; Pandian, S.; Rameshkumar, R.; Subramanian, A.; Ramesh, M. Effects of cefotaxime, amino acids and carbon source on somatic embryogenesis and plant regeneration in four Indian genotypes of foxtail millet (Setaria italica L.). In Vitr. Cell. Dev. Biol. Plant 2015, 52, 140–153. [Google Scholar] [CrossRef]
  4. Saha, S.; Islam, Z.; Islam, S.; Hassan, M.F.; Hossain, M.S.; Islam, S.M.S. Enhancement of Somatic Embryogenesis by Mature and Immature Seeds in Wheat (Triticum aestivum L.). J. Biol. Life Sci. 2017, 8, 20. [Google Scholar] [CrossRef] [Green Version]
  5. Ana, L.M.L.; Ivone, B.d.O.E.; Carlos, A.S.; Claudete, A.M.; Maria, d.F.P.S.M. Somatic embryogenesis and plant regeneration in popcorn (Zea mays L.) inbred lines. Afr. J. Biotechnol. 2017, 16, 1738–1742. [Google Scholar] [CrossRef] [Green Version]
  6. Abe, T. Genotypic variability for callus formation and plant regeneration in rice. Theor. Appl. Genet. 1986, 72, 3–10. [Google Scholar] [CrossRef]
  7. Ge, F.; Qu, J.; Liu, P.; Pan, L.; Zou, C.; Yuan, G.; Yang, C.; Pan, G.; Huang, J.; Ma, L. Genome assembly of the maize inbred line A188 provides a new reference genome for functional genomics. Crop J. 2022, 10, 47–55. [Google Scholar] [CrossRef]
  8. Jia, H.; Yi, D.; Yu, J.; Xue, S.; Ma, Z. Mapping QTLs for tissue culture response of mature wheat embryos. Mol. Cells 2007, 23, 323–330. [Google Scholar]
  9. Zhao, L.; Manli, Y.U.; Guo, H.; Jiayu, G.U.; Xie, Y.; Zhao, S.; Xiong, H.; Liu, L. Preliminary Gene Mapping for Anther Culture Ability Based on a New Wheat Germplasm with High Anther Culture Ability. J. Nucl. Agric. Sci. 2018, 32, 2305–2312. [Google Scholar]
  10. Kwon, Y.S.; Kim, K.M.; Eun, M.Y.; Sohn, J.K. Quantitative trait loci mapping associated with plant regeneration ability from seed derived calli in rice (Oryza sativa L.). Mol. Cells 2001, 11, 64–67. [Google Scholar]
  11. Pan, G.; Zhang, Z.; Wei, X.; Song, Y.; Zhao, M.; Xia, Y.; Rong, T. QTL Analysis of Maize (Zea may L.) Embryo Culturing Capacity. Acta Agron. Sin. 2006, 32, 7–13. [Google Scholar]
  12. Zhang, H.; Liu, Y.; Guo, X.; Zhang, F.; Li, J.; Tan, Z.; Chen, G.; Sun, D. QTL mapping for callus induction and plant regeneration in maize immature embryos. Acta Agron. Sin. 2006, 32, 385–389. [Google Scholar]
  13. Salvo, S.; Cook, J.; Carlson, A.R.; Hirsch, C.N.; Kaeppler, S.M.; Kaeppler, H.F. Genetic Fine-Mapping of a Quantitative Trait Locus (QTL) Associated with Embryogenic Tissue Culture Response and Plant Regeneration Ability in Maize (Zea mays L.). Plant Genome 2018, 11, 170111. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Ma, L.; An, R.; Jiang, L.; Zhang, C.; Li, Z.; Zou, C.; Yang, C.; Pan, G.; Lübberstedt, T.; Shen, Y. Effects of ZmHIPP on lead tolerance in maize seedlings: Novel ideas for soil bioremediation. J. Hazard. Mater. 2022, 430, 128457. [Google Scholar] [CrossRef]
  15. Ledwoń, A.; Gaj, M.D. LEAFY COTYLEDON2 gene expression and auxin treatment in relation to embryogenic capacity of Arabidopsis somatic cells. Plant Cell Rep. 2009, 28, 1677. [Google Scholar] [CrossRef] [PubMed]
  16. Lowe, K.; Wu, E.; Wang, N.; Hoerster, G.; Hastings, C.; Cho, M.-J.; Scelonge, C.; Lenderts, B.; Chamberlin, M.; Cushatt, J.; et al. Morphogenic Regulators Baby boom and Wuschel Improve Monocot Transformation. Plant Cell 2016, 28, 1998–2015. [Google Scholar] [CrossRef] [Green Version]
  17. Zhang, T.-Q.; Lian, H.; Zhou, C.-M.; Xu, L.; Jiao, Y.; Wang, J.-W. A Two-Step Model for de Novo Activation of WUSCHEL during Plant Shoot Regeneration. Plant Cell 2017, 29, 1073–1087. [Google Scholar] [CrossRef] [Green Version]
  18. Ma, L.; Liu, M.; Yan, Y.; Qing, C.; Zhang, X.; Zhang, Y.; Long, Y.; Wang, L.; Pan, L.; Zou, C. Genetic dissection of maize embryonic callus regenerative capacity using multi-locus genome-wide association studies. Front. Plant Sci. 2018, 9, 561. [Google Scholar] [CrossRef] [Green Version]
  19. Ge, F.; Luo, X.; Huang, X.; Zhang, Y.; He, X.; Liu, M.; Lin, H.; Peng, H.; Li, L.; Zhang, Z. Genome-wide analysis of transcription factors involved in maize embryonic callus formation. Physiol. Plant. 2016, 158, 452–462. [Google Scholar] [CrossRef]
  20. Langfelder, P.; Horvath, S. WGCNA: An R package for weighted correlation network analysis. BMC Bioinform. 2008, 9, 1–13. [Google Scholar] [CrossRef] [Green Version]
  21. Jansen, C.; Zhang, Y.; Liu, H.; Gonzalez-Portilla, P.J.; Lauter, N.; Kumar, B.; Trucillo-Silva, I.; Martin, J.P.S.; Lee, M.; Simcox, K.; et al. Genetic and agronomic assessment of cob traits in corn under low and normal nitrogen management conditions. Theor. Appl. Genet. 2015, 128, 1231–1242. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Morino, K.; Matsuda, F.; Miyazawa, H.; Sukegawa, A.; Miyagawa, H.; Wakasa, K. Metabolic profiling of tryptophan-overproducing rice calli that express a feedback-insensitive α subunit of anthranilate synthase. Plant Cell Physiol. 2005, 46, 514–521. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Ge, F.; Hu, H.; Huang, X.; Zhang, Y.; Wang, Y.; Li, Z.; Zou, C.; Peng, H.; Li, L.; Gao, S. Metabolomic and proteomic analysis of maize embryonic callus induced from immature embryo. Sci. Rep. 2017, 7, 1–16. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Du, X.; Fang, T.; Liu, Y.; Huang, L.; Zang, M.; Wang, G.; Liu, Y.; Fu, J. Transcriptome profiling predicts new genes to promote maize callus formation and transformation. Front. Plant Sci. 2019, 10, 1633. [Google Scholar] [CrossRef]
  25. Roos, D.; de Boer, M. Mutations in cis that affect mRNA synthesis, processing and translation. Biochim. Biophys. Acta (BBA) Mol. Basis Dis. 2021, 1867, 166166. [Google Scholar] [CrossRef]
  26. Haddad, L.A. Protein-Coding Genes. In Human Genome Structure, Function and Clinical Considerations; Haddad, L.A., Ed.; Springer International Publishing: Cham, Switzerland, 2021; pp. 93–138. [Google Scholar]
  27. Li, H.; Zhang, L.; Wang, J. Analysis and answers to frequently asked questions in quantitative trait locus mapping. Acta Agron. Sin. 2010, 36, 918–931. [Google Scholar] [CrossRef] [Green Version]
  28. Wang, J.-K. Inclusive composite interval mapping of quantitative trait genes. Acta Agron. Sin. 2009, 35, 239–245. [Google Scholar] [CrossRef] [Green Version]
  29. Zhang, L.; Li, H.; Li, Z.; Wang, J. Interactions Between Markers Can Be Caused by the Dominance Effect of Quantitative Trait Loci. Genetics 2008, 180, 1177–1190. [Google Scholar] [CrossRef] [Green Version]
  30. Krakowsky, M.D.; Lee, M.; Garay, L.; Woodman-Clikeman, W.; Long, M.J.; Sharopova, N.; Frame, B.; Wang, K. Quantitative trait loci for callus initiation and totipotency in maize (Zea mays L.). Theor. Appl. Genet. 2006, 113, 821–830. [Google Scholar] [CrossRef] [Green Version]
  31. Taguchi-Shiobara, F.; Yamamoto, T.; Yano, M.; Oka, S. Mapping QTLs that control the performance of rice tissue culture and evaluation of derived near-isogenic lines. Theor. Appl. Genet. 2006, 112, 968–976. [Google Scholar] [CrossRef]
  32. Li, S.; Yan, S.; Wang, A.h.; Zou, G.; Huang, X.; Han, B.; Qian, Q.; Tao, Y. Identification of QTLs associated with tissue culture response through sequencing-based genotyping of RILs derived from 93-11 × Nipponbare in rice (Oryza sativa). Plant Cell Rep. 2013, 32, 103–116. [Google Scholar] [CrossRef] [PubMed]
  33. Fehér, A. Callus, Dedifferentiation, Totipotency, Somatic Embryogenesis: What These Terms Mean in the Era of Molecular Plant Biology? Front. Plant Sci. 2019, 10, 536. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Etedali, F.; Khossroshahli, M.; Motallebi-Azar, A.; Valizadeh, M.; Kazemiani, S. Gene action and heterosis for callus induction and growth from mature embryo culture of rapeseed (Brassica napus L.). Int. Res. J. Agric. Sci. Soil Sci. 2011, 1, 402–407. [Google Scholar]
  35. Chang, L.; He, K.; Cui, T.; Xue, J.; Liu, J. QTL mapping and QTL × environment interaction analysis of kernel ratio in maize (Zea mays). J. Agric. Biotechnol. 2017, 25, 517–525. [Google Scholar]
  36. van Eeuwijk, F.A.; Bink, M.C.A.M.; Chenu, K.; Chapman, S.C. Detection and use of QTL for complex traits in multiple environments. Curr. Opin. Plant Biol. 2010, 13, 193–205. [Google Scholar] [CrossRef] [PubMed]
  37. Wan, Y.; Rocheford, T.R.; Widholm, J.M. RFLP analysis to identify putative chromosomal regions involved in the anther culture response and callus formation of maize. Theor. Appl. Genet. 1992, 85, 360–365. [Google Scholar] [CrossRef] [PubMed]
  38. Armstrong, C.L.; Romero-Severson, J.; Hodges, T.K. Improved tissue culture response of an elite maize inbred through backcross breeding, and identification of chromosomal regions important for regeneration by RFLP analysis. Theor. Appl. Genet. 1992, 84, 755–762. [Google Scholar] [CrossRef]
  39. Zhang, S.; Liu, X.; Lin, Y.; Xie, G.; Fu, F.; Liu, H.; Wang, J.; Gao, S.; Lan, H.; Rong, T. Characterization of a ZmSERK gene and its relationship to somatic embryogenesis in a maize culture. Plant Cell Tissue Organ Cult. 2011, 105, 29–37. [Google Scholar] [CrossRef]
  40. Glawischnig, E.; Tomas, A.; Eisenreich, W.; Spiteller, P.; Bacher, A.; Gierl, A. Auxin biosynthesis in maize kernels. Plant Physiol. 2000, 123, 1109–1120. [Google Scholar] [CrossRef] [Green Version]
  41. Guilfoyle, T.; Hagen, G. Auxin signal transduction. Essays Biochem. 2015, 58, 1–12. [Google Scholar] [CrossRef]
  42. Su, Y.H.; Zhao, X.Y.; Liu, Y.B.; Zhang, C.L.; O’Neill, S.D.; Zhang, X.S. Auxin-induced WUS expression is essential for embryonic stem cell renewal during somatic embryogenesis in Arabidopsis. Plant J. 2009, 59, 448–460. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. Nic-Can, G.I.; Loyola-Vargas, V.M. The Role of the Auxins During Somatic Embryogenesis. In Somatic Embryogenesis: Fundamental Aspects and Applications; Loyola-Vargas, V.M., Ochoa-Alejo, N., Eds.; Springer International Publishing: Cham, Switzerland, 2016; pp. 171–182. [Google Scholar]
  44. Wójcik, A.M.; Wójcikowska, B.; Gaj, M.D. Current Perspectives on the Auxin-Mediated Genetic Network that Controls the Induction of Somatic Embryogenesis in Plants. Int. J. Mol. Sci. 2020, 21, 1333. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Srivastava, S.; Dwivedi, U.N. Plant regeneration from callus of Cuscuta reflexa–an angiospermic parasite and modulation of catalase and peroxidase activity by salicylic acid and naphthalene acetic acid. Plant Physiol. Biochem. 2001, 39, 529–538. [Google Scholar] [CrossRef]
  46. Hao, L.; Zhou, L.; Xu, X.; Cao, J.; Xi, T. The role of salicylic acid and carrot embryogenic callus extracts in somatic embryogenesis of naked oat (Avena nuda). Plant Cell Tissue Organ Cult. 2006, 85, 109–113. [Google Scholar] [CrossRef]
  47. Gadzovska, S.; Maury, S.; Delaunay, A.; Spasenoski, M.; Hagège, D.; Courtois, D.; Joseph, C. The influence of salicylic acid elicitation of shoots, callus, and cell suspension cultures on production of naphtodianthrones and phenylpropanoids in Hypericum perforatum L. Plant Cell Tissue Organ Cult. 2013, 113, 25–39. [Google Scholar] [CrossRef]
  48. Armstrong, C.; Green, C. Establishment and maintenance of friable, embryogenic maize callus and the involvement of L-proline. Planta 1985, 164, 207–214. [Google Scholar] [CrossRef] [PubMed]
  49. Knapp, S.; Stroup, W.; Ross, W. Exact confidence intervals for heritability on a progeny mean basis 1. Crop Sci. 1985, 25, 192–194. [Google Scholar] [CrossRef]
  50. Liu, H.; Niu, Y.; Gonzalez-Portilla, P.J.; Zhou, H.; Wang, L.; Zuo, T.; Qin, C.; Tai, S.; Jansen, C.; Shen, Y.; et al. An ultra-high-density map as a community resource for discerning the genetic basis of quantitative traits in maize. BMC Genom. 2015, 16, 1078. [Google Scholar] [CrossRef] [Green Version]
  51. Meng, L.; Li, H.; Zhang, L.; Wang, J. QTL IciMapping: Integrated software for genetic linkage map construction and quantitative trait locus mapping in biparental populations. Crop J. 2015, 3, 269–283. [Google Scholar] [CrossRef] [Green Version]
  52. Sun, Z.; Li, H.; Zhang, L.; Wang, J. Properties of the test statistic under null hypothesis and the calculation of LOD threshold in quantitative trait loci (QTL) mapping. Acta Agron. Sin. 2013, 39, 1–11. [Google Scholar] [CrossRef]
  53. Saghai-Maroof, M.A.; Soliman, K.M.; Jorgensen, R.A.; Allard, R. Ribosomal DNA spacer-length polymorphisms in barley: Mendelian inheritance, chromosomal location, and population dynamics. Proc. Natl. Acad. Sci. USA 1984, 81, 8014–8018. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  54. Chen, S.; Zhou, Y.; Chen, Y.; Gu, J. fastp: An ultra-fast all-in-one FASTQ preprocessor. Bioinformatics 2018, 34, i884–i890. [Google Scholar] [CrossRef] [PubMed]
  55. Shannon, P.; Markiel, A.; Ozier, O.; Baliga, N.S.; Wang, J.T.; Ramage, D.; Amin, N.; Schwikowski, B.; Ideker, T. Cytoscape: A software environment for integrated models of biomolecular interaction networks. Genom. Res. 2003, 13, 2498–2504. [Google Scholar] [CrossRef] [PubMed]
Ijms 23 08786 g001 550
Figure 1. The calli induced from different doubled haploid lines in the IBM Syn10 DH population. (a,c) show the embryogenic callus (EC). (b,d) exhibit the non-embryogenic callus (NEC). Bar scale is 1 cm in (a,b) and 3 mm in (c,d).
Figure 1. The calli induced from different doubled haploid lines in the IBM Syn10 DH population. (a,c) show the embryogenic callus (EC). (b,d) exhibit the non-embryogenic callus (NEC). Bar scale is 1 cm in (a,b) and 3 mm in (c,d).
Ijms 23 08786 g001
Ijms 23 08786 g002 550
Figure 2. Phenotypic frequency distribution of the IBM Syn10 DH population across three environments. RSF, LS, DC, and REC represent the rate of shoot formation, length of shoot, diameter of callus, and rate of embryogenic callus induction. XSBN, YA, and CZ denote the environments Xishuangbanna, Ya’an, and Chongzhou, respectively. Mo17 and B73 show the two parents of the population.
Figure 2. Phenotypic frequency distribution of the IBM Syn10 DH population across three environments. RSF, LS, DC, and REC represent the rate of shoot formation, length of shoot, diameter of callus, and rate of embryogenic callus induction. XSBN, YA, and CZ denote the environments Xishuangbanna, Ya’an, and Chongzhou, respectively. Mo17 and B73 show the two parents of the population.
Ijms 23 08786 g002
Ijms 23 08786 g003 550
Figure 3. Chromosomal distributions of identified QTLs under three environments in this study. The red, blue, green, and purple represent Xishuangbanna (XSBN), Ya’an (YA), Chongzhou (CZ), and the best linear unbiased prediction (BLUP), respectively. Squares, circles, triangles, and pentagrams represent diameter of callus (DC), rate of shoot formation (RSF), length of shoot (LS), and rate of embryogenic callus induction (REC), respectively.
Figure 3. Chromosomal distributions of identified QTLs under three environments in this study. The red, blue, green, and purple represent Xishuangbanna (XSBN), Ya’an (YA), Chongzhou (CZ), and the best linear unbiased prediction (BLUP), respectively. Squares, circles, triangles, and pentagrams represent diameter of callus (DC), rate of shoot formation (RSF), length of shoot (LS), and rate of embryogenic callus induction (REC), respectively.
Ijms 23 08786 g003
Ijms 23 08786 g004 550
Figure 4. Gene co-expression network and hub genes responsible for maize embryogenic callus induction. (a) Eigengene expression pattern of the embryogenic callus induction-associated module (blue). Upper panel shows the gene expression heat map; lower panel shows the eigengene expression histogram. (b) Top 20 enriched pathways in blue module. (c) Co-expression network of hub genes. The four hub genes are shown in pink and orange. (d) The expression patterns of the four hub genes during embryogenic callus induction. CN9802, ZM28, JS0251, and YA3237 represent four different maize inbred lines. The significance of difference between each pair of these lines was estimated by using a t-test (2-tailed), with n = 3 (three biological replicates). *, **, and *** represent significant difference at p = 0.05, p = 0.01, and p = 0.001 levels, respectively, at least between one pair of the four lines, for each stage. ns = not significant. The bars show the standard deviation (SD).
Figure 4. Gene co-expression network and hub genes responsible for maize embryogenic callus induction. (a) Eigengene expression pattern of the embryogenic callus induction-associated module (blue). Upper panel shows the gene expression heat map; lower panel shows the eigengene expression histogram. (b) Top 20 enriched pathways in blue module. (c) Co-expression network of hub genes. The four hub genes are shown in pink and orange. (d) The expression patterns of the four hub genes during embryogenic callus induction. CN9802, ZM28, JS0251, and YA3237 represent four different maize inbred lines. The significance of difference between each pair of these lines was estimated by using a t-test (2-tailed), with n = 3 (three biological replicates). *, **, and *** represent significant difference at p = 0.05, p = 0.01, and p = 0.001 levels, respectively, at least between one pair of the four lines, for each stage. ns = not significant. The bars show the standard deviation (SD).
Ijms 23 08786 g004
Ijms 23 08786 g005 550
Figure 5. Hub gene-based association mapping. (a,f) Manhattan plots (top) and LD heat maps (bottom) of Zm00001d034388 and Zm00001d028477. The significance threshold is shown as a red dashed line (p < 0.05). The gene structure is shown in the middle. The exons are represented by filled red boxes. (b–e) Differences in (b) REC (rate of embryogenic callus), (c) DC (diameter of callus), (d) RSF (rate of shoot formation), and (e) LS (length of shoot) between alleles or haplotypes of Zm00001d034388. (g,h) Differences in (g) LS and (h) RSF between alleles or haplotypes of Zm00001d028477. Statistical significance is determined by two-tailed t-test: *, **, and *** denote significance at p = 0.05, p = 0.01, and p = 0.001 levels, respectively.
Figure 5. Hub gene-based association mapping. (a,f) Manhattan plots (top) and LD heat maps (bottom) of Zm00001d034388 and Zm00001d028477. The significance threshold is shown as a red dashed line (p < 0.05). The gene structure is shown in the middle. The exons are represented by filled red boxes. (b–e) Differences in (b) REC (rate of embryogenic callus), (c) DC (diameter of callus), (d) RSF (rate of shoot formation), and (e) LS (length of shoot) between alleles or haplotypes of Zm00001d034388. (g,h) Differences in (g) LS and (h) RSF between alleles or haplotypes of Zm00001d028477. Statistical significance is determined by two-tailed t-test: *, **, and *** denote significance at p = 0.05, p = 0.01, and p = 0.001 levels, respectively.
Ijms 23 08786 g005
Table
Table 1. Phenotypic performances of four traits in the IBM Syn10 DH population and its parents.
Table 1. Phenotypic performances of four traits in the IBM Syn10 DH population and its parents.
Trait aEnvironment bParents (n = 6, t-Test)IBM Syn10 DH Population
Mo17 ± SDB73 ± SDRangeAverageSDCV%SkewnessKurtosis
RSF (%)XSBN81.34 ± 4.0095.77 ± 0.64 **0.00–99.0767.6424.990.37−0.78−0.42
YA81.82 ± 6.2496.67 ± 1.96 **0.00–100.0070.8425.530.36−0.86−0.25
CZ79.06 ± 2.3594.44 ± 5.01 **0.00–100.0070.2724.870.35−1.020.26
LS (mm)XSBN4.0 ± 0.010.3 ± 0.6 **0.0–15.05.52.650.480.700.20
YA4.3 ± 0.311.0 ± 1.0 **0.0–12.06.02.540.420.26−0.32
CZ5.0 ± 0.010.2 ± 1.8 **0.0–11.75.32.290.430.380.13
DC (mm)XSBN2.7 ± 0.33.8 ± 0.3 **1.5–4.32.90.660.230.34−0.85
YA3.0 ± 0.03.8 ± 0.3 **2.0–5.03.00.560.190.33−0.50
CZ3.3 ± 0.33.8 ± 0.32.0–4.02.80.460.160.150.39
REC (%)XSBN6.94 ± 1.390.00 ± 0.00 **0.00–52.174.398.792.003.0810.54
YA7.05 ± 1.910.00 ± 0.00 **0.00–58.612.747.952.904.4822.56
CZ3.39 ± 1.230.00 ± 0.00 **0.00–56.482.276.032.674.8432.70
SD, standard deviation, CV, coefficient of variation. ** The difference is significant at the p = 0.01 level between B73 and Mo17. a RSF, LS, DC, and REC represent the rate of shoot formation, length of shoot, diameter of callus, and rate of embryogenic callus induction. b XSBN, YA, and CZ represent the three environments, namely Xishuangbanna, Ya’an, and Chongzhou, respectively.
Table
Table 2. Phenotypic correlation between the traits across the three environments.
Table 2. Phenotypic correlation between the traits across the three environments.
EnvironmentTraitREC
(n = 210)		SC
(n = 210)		LS
(n = 210)
XSBNRSF−0.310 **0.531 **0.692 **
LS−0.263 **0.725 **
DC−0.304 **
YARSF−0.342 **0.482 **0.693 **
LS−0.243 **0.612 **
DC−0.318 **
CZRSF−0.1050.571 **0.737 **
LS−0.134 *0.639 **
DC−0.177 *
** and * represent significant correlation at the 0.01 and 0.05 levels by a t-test (2-tailed).
Table
Table 3. Analysis of variance (ANOVA) and broad-sense heritability (H2) for four traits under three environments.
Table 3. Analysis of variance (ANOVA) and broad-sense heritability (H2) for four traits under three environments.
TraitSource of VariationdfMean SquareSignificanceH2 (%)
RSFGenotype(G)2092092.792<0.01 **84.02
Environment(E)21136.991<0.01 **
G × E4181147.495<0.01 **
Error1260140.281
LSGenotype(G)20922.363<0.01 **85.64
Environment(E)240.006<0.01 **
G × E41810.892<0.01 **
Error12601.067
DCGenotype(G)2091.076<0.01 **83.17
Environment(E)23.584<0.01 **
G × E4180.633<0.01 **
Error12600.061
RECGenotype(G)209194.825<0.01 **88.27
Environment(E)2369.202<0.01 **
G × E41873.082<0.01 **
Error126013.689
RSF, LS, DC, and REC represent the rate of shoot formation, length of shoot, diameter of callus, and rate of embryogenic callus induction. H2 represents the broad-sense heritability. ** denotes significance at p = 0.01 level by a t-test (2-tailed).
Table
Table 4. QTLs controlling embryogenic callus induction-related traits across three environments.
Table 4. QTLs controlling embryogenic callus induction-related traits across three environments.
TraitNameEnv.Chr.Genetic Position (cM)Physical Position (Mb)LOD aPVE bADD c
DCqDC1-1CZ1125.84191.425–191.6003.027.10−0.1174
qDC1-2CZ1217.21286.200–286.3002.615.730.1089
qDC1-2BLUP1217.49286.400–286.5503.554.640.0398
qDC3-1YA39.852.700–2.8003.196.59−0.1322
qDC4-1YA4101.97170.300–170.6753.407.130.1407
qDC4-2YA4116.7180.300–180.4004.769.98−0.1705
qDC4-3XSBN4138.78210.875–211.0002.855.300.1734
qDC7-1YA771.29132.275–132.5253.216.620.1437
qDC7-1BLUP772.01132.800–132.9004.846.400.0501
qDC7-2BLUP7140.95174.175–174.3003.424.61−0.0391
qDC9-1BLUP986.05133.775–133.9005.297.130.0504
qDC9-2XSBN999.62141.700–141.9504.147.890.1858
qDC9-3CZ9131.38152.200–152.3003.929.490.1512
qDC9-3BLUP9132.23152.600–152.7004.065.400.047
qDC10-1XSBN1027.827.400–7.7754.007.790.1936
qDC10-2XSBN1041.1367.250–68.9003.817.17−0.184
RSFqRSF1-1CZ1127.16192.525–193.6504.297.37−6.7656
qRSF1-2YA1211.45280.975–281.1003.467.59−6.9403
qRSF1-3XSBN1230.48290.700–290.8003.216.49−6.7815
qRSF2-1YA275.7932.650–33.1752.685.33−5.8722
qRSF2-2YA2190.64232.500–232.6004.759.95−8.0007
qRSF3-1CZ381.31155.675–155.6755.719.978.2628
qRSF4-1CZ45.092.500–2.6003.776.36−6.4367
qRSF6-1CZ617.1313.600–13.7253.315.656.0506
qRSF8-1XSBN80.060.100–0.3502.825.556.1132
qRSF8-2YA8115.14169.300–169.4003.306.61−6.8164
qRSF9-1BLUP986.7133.900–134.0004.496.892.655
qRSF9-1XSBN987.33134.000–134.1003.737.467.3017
qRSF10-1BLUP1095.16145.200–145.3004.307.13−2.8252
LSqLS1-1BLUP158.2623.675–23.9256.839.740.3231
qLS1-1XSBN160.7724.975–25.2257.7011.540.946
qLS1-2XSBN171.2735.375–35.6253.154.44−0.5845
qLS1-3CZ1127.16192.525–193.6505.5910.68−0.8065
qLS1-4XSBN1231.2290.800–290.9003.515.03−0.6832
qLS2-1BLUP2115.59185.350–185.5503.544.900.2473
qLS3-1YA317.744.100–4.2002.986.290.6625
qLS3-2BLUP352.4114.300–14.4006.339.20−0.8419
qLS3-2XSBN352.6814.650–14.8505.147.34−0.2805
qLS3-3BLUP3161.37216.325–216.5503.154.420.2312
qLS9-1CZ977.71117.825–118.2003.276.05−0.6074
qLS10-1BLUP1087.58143.675–143.8002.904.01−0.2087
qLS10-1XSBN1091.16144.600–144.7003.965.90−0.7004
qLS10-1YA1096.74145.400–145.5003.246.98−0.741
RECqREC1-1BLUP189.6166.800–66.9003.836.251.0291
qREC4-1XSBN414.294.400–4.5003.1811.452.3818
qREC5-1YA558.421.200–21.3003.137.73−2.2782
qREC6-1CZ6132.21164.300–164.4002.556.00−1.5787
qREC7-1BLUP723.375.600–5.7002.514.01−0.8092
qREC7-2BLUP7120.07167.200–167.3003.966.36−0.9828
a LOD represents the logarithm of odds score. b PVE denotes the phenotypic variance explained by each QTL. c ADD, additive effect values of the QTLs. Positive and negative values indicate that the increasing effect of phenotypic value is derived from the alleles of B73 and Mo17, respectively. Env. and Chr. represent environment and chromosome, respectively. XSBN, YA, and CZ denote the environments, Xishuangbanna, Ya’an, and Chongzhou, respectively. BLUP, the best linear unbiased prediction.
Table
Table 5. QTL clusters for the embryogenic callus induction traits in the IBM Syn10 DH population across three environments.
Table 5. QTL clusters for the embryogenic callus induction traits in the IBM Syn10 DH population across three environments.
QTL Cluster NumberChromosomeTraits aQTL NamesPosition (Mb)Positive AllelesRange of Explained Phenotypic Variation (%)
a1DC + LS + RSFqDC1-1; qRSF1-1; qLS1-3191.425–193.650Mo17 + Mo17 + Mo177.10–10.68
b1RSF + LSqRSF1-3; qLS1-4290.700–290.900Mo17 + Mo175.03–6.49
c9DC + RSFqDC9-1; qRSF9-1133.775–134.100B73 + B73 + B736.89–7.46
d10RSF + LSqRSF10-1; qLS10-1143.675–145.500Mo17 + Mo17 + Mo17 + Mo174.01–7.13
a RSF, LS, and DC represent the rate of shoot formation, length of shoot, and diameter of callus, respectively.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Long, Y.; Liang, T.; Ma, L.; Liu, P.; Yang, Y.; Zhang, X.; Zou, C.; Zhang, M.; Ge, F.; Yuan, G.; et al. Combined QTL Mapping across Multiple Environments and Co-Expression Network Analysis Identified Key Genes for Embryogenic Callus Induction from Immature Maize Embryos. Int. J. Mol. Sci. 2022, 23, 8786. https://doi.org/10.3390/ijms23158786

AMA Style

Long Y, Liang T, Ma L, Liu P, Yang Y, Zhang X, Zou C, Zhang M, Ge F, Yuan G, et al. Combined QTL Mapping across Multiple Environments and Co-Expression Network Analysis Identified Key Genes for Embryogenic Callus Induction from Immature Maize Embryos. International Journal of Molecular Sciences. 2022; 23(15):8786. https://doi.org/10.3390/ijms23158786

Chicago/Turabian Style

Long, Yun, Tianhu Liang, Langlang Ma, Peng Liu, Yun Yang, Xiaoling Zhang, Chaoying Zou, Minyan Zhang, Fei Ge, Guangsheng Yuan, and et al. 2022. "Combined QTL Mapping across Multiple Environments and Co-Expression Network Analysis Identified Key Genes for Embryogenic Callus Induction from Immature Maize Embryos" International Journal of Molecular Sciences 23, no. 15: 8786. https://doi.org/10.3390/ijms23158786

APA Style

Long, Y., Liang, T., Ma, L., Liu, P., Yang, Y., Zhang, X., Zou, C., Zhang, M., Ge, F., Yuan, G., Lübberstedt, T., Pan, G., & Shen, Y. (2022). Combined QTL Mapping across Multiple Environments and Co-Expression Network Analysis Identified Key Genes for Embryogenic Callus Induction from Immature Maize Embryos. International Journal of Molecular Sciences, 23(15), 8786. https://doi.org/10.3390/ijms23158786

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