Functional genomics of Brassica napus: Progresses, challenges, and perspectives.

Brassica napus, commonly known as rapeseed or canola, is a major oil crop contributing over 13% to the stable supply of edible vegetable oil worldwide. Identification and understanding the gene functions in the B. napus genome is crucial for genomic breeding. A group of genes controlling agronomic traits have been successfully cloned through functional genomics studies in B. napus. In this review, we present an overview of the progress made in the functional genomics of B. napus, including the availability of germplasm resources, omics databases and cloned functional genes. Based on the current progress, we also highlight the main challenges and perspectives in this field. The advances in the functional genomics of B. napus contribute to a better understanding of the genetic basis underlying the complex agronomic traits in B. napus and will expedite the breeding of high quality, high resistance and high yield in B. napus varieties.

Lipidomic and metabolomic analyses reveal changes of lipids and metabolites profile in rapeseed during nitrogen deficiency.

Nitrogen is one of the most essential macronutrients for plant growth and its availability in soil is vital for agricultural sustainability and productivity. However, excessive nitrogen application could reduce the nitrogen use efficiency and produce environmental pollution. Here, we systematically determined the response in lipidome and metabolome in rapeseed during nitrogen starvation. Plant growth was severely retarded during nitrogen deficiency, while the levels of most amino acids was significantly decreased. The levels of monogalactosyl diacyglycerol (MGDG) in leaves and roots was significantly decreased, while the level of digalactosyl diacylglycerol (DGDG) was significantly decreased in roots, resulting in significant reduction of MGDG/DGDG ratio during nitrogen starvation. Meanwhile, the levels of sulfoquinovosyl diacylglycerol, phosphatidylglycerol and glucuronosyl diacylglycerol was reduced to varying extents. Moreover, the levels of metabolites in the tricarboxylic acid cycle, Calvin cycle, and energy metabolism was changed during nitrogen deficiency. These findings show that nitrogen deprivation alters the membrane lipid metabolism and carbon metabolism, and our study provides valuable information to further understand the response of rapeseed to nitrogen deficiency at metabolism level.

Characterization of novel loci controlling seed oil content in Brassica napus by marker metabolite-based multi-omics analysis.

BACKGROUND: Seed oil content is an important agronomic trait of Brassica napus (B. napus), and metabolites are considered as the bridge between genotype and phenotype for physical traits. RESULTS: Using a widely targeted metabolomics analysis in a natural population of 388 B. napus inbred lines, we quantify 2172 metabolites in mature seeds by liquid chromatography mass spectrometry, in which 131 marker metabolites are identified to be correlated with seed oil content. These metabolites are then selected for further metabolite genome-wide association study and metabolite transcriptome-wide association study. Combined with weighted correlation network analysis, we construct a triple relationship network, which includes 21,000 edges and 4384 nodes among metabolites, metabolite quantitative trait loci, genes, and co-expression modules. We validate the function of BnaA03.TT4, BnaC02.TT4, and BnaC05.UK, three candidate genes predicted by multi-omics analysis, which show significant impacts on seed oil content through regulating flavonoid metabolism in B. napus. CONCLUSIONS: This study demonstrates the advantage of utilizing marker metabolites integrated with multi-omics analysis to dissect the genetic basis of agronomic traits in crops.

Regulation of seed oil accumulation by lncRNAs in Brassica napus.

BACKGROUND: Studies have indicated that long non-coding RNAs (lncRNAs) play important regulatory roles in many biological processes. However, the regulation of seed oil biosynthesis by lncRNAs remains largely unknown. RESULTS: We comprehensively identified and characterized the lncRNAs from seeds in three developing stages in two accessions of Brassica napus (B. napus), ZS11 (high oil content) and WH5557 (low oil content). Finally, 8094 expressed lncRNAs were identified. LncRNAs MSTRG.22563 and MSTRG.86004 were predicted to be related to seed oil accumulation. Experimental results show that the seed oil content is decreased by 3.1-3.9% in MSTRG.22563 overexpression plants, while increased about 2% in MSTRG.86004, compared to WT. Further study showed that most genes related to lipid metabolism had much lower expression, and the content of some metabolites in the processes of respiration and TCA (tricarboxylic acid) cycle was reduced in MSTRG.22563 transgenic seeds. The expression of genes involved in fatty acid synthesis and seed embryonic development (e.g., LEC1) was increased, but genes related to TAG assembly was decreased in MSTRG.86004 transgenic seeds. CONCLUSION: Our results suggest that MSTRG.22563 might impact seed oil content by affecting the respiration and TCA cycle, while MSTRG.86004 plays a role in prolonging the seed developmental time to increase seed oil accumulation.

High-throughput phenotyping-based quantitative trait loci mapping reveals the genetic architecture of the salt stress tolerance of Brassica napus.

Salt stress is a major limiting factor that severely affects the survival and growth of crops. It is important to understand the salt stress tolerance ability of Brassica napus and explore the underlying related genetic resources. We used a high-throughput phenotyping platform to quantify 2111 image-based traits (i-traits) of a natural population under three different salt stress conditions and an intervarietal substitution line (ISL) population under nine different stress conditions to monitor and evaluate the salt stress tolerance of B. napus over time. We finally identified 928 high-quality i-traits associated with the salt stress tolerance of B. napus. Moreover, we mapped the salt stress-related loci in the natural population via a genome-wide association study and performed a linkage analysis associated with the ISL population, respectively. These results revealed 234 candidate genes associated with salt stress response, and two novel candidate genes, BnCKX5 and BnERF3, were experimentally verified to regulate the salt stress tolerance of B. napus. This study demonstrates the feasibility of using high-throughput phenotyping-based quantitative trait loci mapping to accurately and comprehensively quantify i-traits associated with B. napus. The mapped loci could be used for genomics-assisted breeding to genetically improve the salt stress tolerance of B. napus.

Comprehensive transcriptional variability analysis reveals gene networks regulating seed oil content of Brassica napus.

BACKGROUND: Regulation of gene expression plays an essential role in controlling the phenotypes of plants. Brassica napus (B. napus) is an important source for the vegetable oil in the world, and the seed oil content is an important trait of B. napus. RESULTS: We perform a comprehensive analysis of the transcriptional variability in the seeds of B. napus at two developmental stages, 20 and 40 days after flowering (DAF). We detect 53,759 and 53,550 independent expression quantitative trait loci (eQTLs) for 79,605 and 76,713 expressed genes at 20 and 40 DAF, respectively. Among them, the local eQTLs are mapped to the adjacent genes more frequently. The adjacent gene pairs are regulated by local eQTLs with the same open chromatin state and show a stronger mode of expression piggybacking. Inter-subgenomic analysis indicates that there is a feedback regulation for the homoeologous gene pairs to maintain partial expression dosage. We also identify 141 eQTL hotspots and find that hotspot87-88 co-localizes with a QTL for the seed oil content. To further resolve the regulatory network of this eQTL hotspot, we construct the XGBoost model using 856 RNA-seq datasets and the Basenji model using 59 ATAC-seq datasets. Using these two models, we predict the mechanisms affecting the seed oil content regulated by hotspot87-88 and experimentally validate that the transcription factors, NAC13 and SCL31, positively regulate the seed oil content. CONCLUSIONS: We comprehensively characterize the gene regulatory features in the seeds of B. napus and reveal the gene networks regulating the seed oil content of B. napus.

A NAC-EXPANSIN module enhances maize kernel size by controlling nucellus elimination.

Maize early endosperm development is initiated in coordination with elimination of maternal nucellar tissues. However, the underlying mechanisms are largely unknown. Here, we characterize a major quantitative trait locus for maize kernel size and weight that encodes an EXPANSIN gene, ZmEXPB15. The encoded ß-expansin protein is expressed specifically in nucellus, and positively controls kernel size and weight by promoting nucellus elimination. We further show that two nucellus-enriched transcription factors (TFs), ZmNAC11 and ZmNAC29, activate ZmEXPB15 expression. Accordingly, these two TFs also promote kernel size and weight through nucellus elimination regulation, and genetic analyses support their interaction with ZmEXPB15. Importantly, hybrids derived from a ZmEXPB15 overexpression line have increased kernel weight, demonstrates its potential value in breeding. Together, we reveal a pathway modulating the cellular processes of maternal nucellus elimination and early endosperm development, and an approach to improve kernel weight.

Combining Physio-Biochemical Characterization and Transcriptome Analysis Reveal the Responses to Varying Degrees of Drought Stress in Brassica napus L.

Brassica napus L. has become one of the most important oil-bearing crops, and drought stress severely influences its yield and quality. By combining physio-biochemical characterization and transcriptome analysis, we studied the response of B. napus plants to different degrees of drought stress. Some physio-biochemical traits, such as fresh weight (FW), dry weight (DW), abscisic acid (ABA) content, net photosynthetic rate (Pn), stomatal conductance (gs), and transpiration rate (Tr), were measured, and the total content of the epidermal wax/cutin, as well as their compositions, was determined. The results suggest that both stomatal transpiration and cuticular transpiration are affected when B. napus plants are subjected to varying degrees of drought stress. A total of 795 up-regulated genes and 1050 down-regulated genes were identified under severe drought stress by transcriptome analysis. Gene ontology (GO) enrichment analysis of differentially expressed genes (DEGs) revealed that the up-regulated genes were mainly enriched in the stress response processes, such as response to water deprivation and abscisic acid, while the down-regulated genes were mainly enriched in the chloroplast-related parts affecting photosynthesis. Moreover, overexpression of BnaA01.CIPK6, an up-regulated DEG, was found to confer drought tolerance in B. napus. Our study lays a foundation for a better understanding of the molecular mechanisms underlying drought tolerance in B. napus.

Genome-Wide Association Studies of Salt-Alkali Tolerance at Seedling and Mature Stages in Brassica napus.

Most plants are sensitive to salt-alkali stress, and the degree of tolerance to salt-alkali stress varies from different species and varieties. In order to explore the salt-alkali stress adaptability of Brassica napus, we collected the phenotypic data of 505 B. napus accessions at seedling and mature stages under control, low and high salt-alkali soil stress conditions in Inner Mongolia of China. Six resistant and 5 sensitive materials, respectively, have been identified both in Inner Mongolia and Xinjiang Uygur Autonomous Region of China. Genome-wide association studies (GWAS) for 15 absolute values and 10 tolerance coefficients (TCs) of growth and agronomic traits were applied to investigate the genetic basis of salt-alkali tolerance of B. napus. We finally mapped 9 significant QTLs related to salt-alkali stress response and predicted 20 candidate genes related to salt-alkali stress tolerance. Some important candidate genes, including BnABA4, BnBBX14, BnVTI12, BnPYL8, and BnCRR1, were identified by combining sequence variation annotation and expression differences. The identified valuable loci and germplasms could be useful for breeding salt-alkali-tolerant B.napus varieties. This study laid a foundation for understanding molecular mechanism of salt-alkali stress adaptation and provides rich genetic resources for the large-scale production of B. napus on salt-alkali land in the future.

Multi-omics analysis dissects the genetic architecture of seed coat content in Brassica napus.

BACKGROUND: Brassica napus is an important vegetable oil source worldwide. Seed coat content is a complex quantitative trait that negatively correlates with the seed oil content in B. napus. RESULTS: Here we provide insights into the genetic basis of natural variation of seed coat content by transcriptome-wide association studies (TWAS) and genome-wide association studies (GWAS) using 382 B. napus accessions. By population transcriptomic analysis, we identify more than 700 genes and four gene modules that are significantly associated with seed coat content. We also characterize three reliable quantitative trait loci (QTLs) controlling seed coat content by GWAS. Combining TWAS and correlation networks of seed coat content-related gene modules, we find that BnaC07.CCR-LIKE (CCRL) and BnaTT8s play key roles in the determination of the trait by modulating lignin biosynthesis. By expression GWAS analysis, we identify a regulatory hotspot on chromosome A09, which is involved in controlling seed coat content through BnaC07.CCRL and BnaTT8s. We then predict the downstream genes regulated by BnaTT8s using multi-omics datasets. We further experimentally validate that BnaCCRL and BnaTT8 positively regulate seed coat content and lignin content. BnaCCRL represents a novel identified gene involved in seed coat development. Furthermore, we also predict the key genes regulating carbon allocation between phenylpropane compounds and oil during seed development in B. napus. CONCLUSIONS: This study helps us to better understand the complex machinery of seed coat development and provides a genetic resource for genetic improvement of seed coat content in B. napus breeding.

Comparative transcriptome analysis of genes and metabolic pathways involved in sporulation in Ganoderma lingzhi.

Over the past decades, Ganoderma lingzhi spores have received considerable attention as a great potential pharmaceutical resource. However, the genetic regulation of sporulation is not well understood. In this study, a comparative transcriptome analysis of the low-sporing HZ203 and high-sporing YW-1 was performed to characterize the mechanism underlying sporulation. A total of 917 differentially expressed genes were identified in HZ203 and 1,450 differentially expressed genes in YW-1. Differentially expressed genes involved in sporulation were identified, which included HOP1, Mek1, MSH4, MSH5, and Spo5 in meiosis. Positive regulatory pathways of sporulation were proposed as 2 transcriptional factors had high connectivity with MSH4 and Spo5. Furthermore, we found that the pathways associated with energy production were enriched in the high-sporing genotype, such as the glyoxylate and dicarboxylate metabolism, starch and sucrose metabolism. Finally, we performed a weighted gene coexpression network analysis and found that the hub genes of the module which exhibit strong positive relationship with the high-sporing phase purportedly participate in signal transduction, carbohydrate transport and metabolism. The dissection of differentially expressed genes during sporulation extends our knowledge about the genetic and molecular networks mediating spore morphogenesis and sheds light on the importance of energy source during sporulation.

Genome- and transcriptome-wide association studies reveal the genetic basis and the breeding history of seed glucosinolate content in Brassica napus.

A high content of seed glucosinolates and their degradation products imposes anti-nutritional effects on livestock; therefore, persistent efforts are made to reduce the seed GSL content to increase the commercial value of rapeseed meal. Here, we dissected the genetic structure of SGC by genome-wide association studies (GWAS) combined with transcriptome-wide association studies (TWAS). Fifteen reliable quantitative trait loci (QTLs) were identified to be associated with the reduced SGC in modern B. napus cultivars by GWAS. Analysis of the selection strength and haplotypes at these QTLs revealed that low SGC was predominantly generated by the co-selection of qGSL.A02.2, qGSL.C02.1, qGSL.A09.2, and qGSL.C09.1. Integration of the results from TWAS, comprehensive bioinformatics, and POCKET algorithm analyses indicated that BnaC02.GTR2 (BnaC02g42260D) is a candidate gene underlying qGSL.C02.1. Using CRISPR/Cas9-derived Bna.gtr2s knockout mutants, we experimentally verified that both BnaC02.GTR2 and its three paralogs positively regulate seed GSL accumulation but negatively regulated vegetative tissue GSL contents. In addition, we observed smaller seeds with higher seed oil content in these Bna.gtr2 mutants. Furthermore, both RNA-seq and correlation analyses suggested that Bna.GTR2s might play a comprehensive role in seed development, such as amino acid accumulation, GSL synthesis, sugar assimilation, and oil accumulation. This study unravels the breeding selection history of low-SGC improvement and provides new insights into the molecular function of Bna.GTR2s in both seed GSL accumulation and seed development in B. napus.

Transcriptome Analysis Reveals Genes of Flooding-Tolerant and Flooding-Sensitive Rapeseeds Differentially Respond to Flooding at the Germination Stage.

Flooding results in significant crop yield losses due to exposure of plants to hypoxic stress. Various studies have reported the effect of flooding stress at seedling establishment or later stages. However, the molecular mechanism prevailing at the germination stage under flooding stress remains enigmatic. The present study highlights the comparative transcriptome analysis in two rapeseed lines, i.e., flooding-tolerant (Santana) and -sensitive (23651) lines under control and 6-h flooding treatments at the germination stage. A total of 1840 up-regulated and 1301 down-regulated genes were shared by both lines in response to flooding. There were 4410 differentially expressed genes (DEGs) with increased expression and 4271 DEGs with reduced expression shared in both control and flooding conditions. Gene ontology (GO) enrichment analysis revealed that "transcription regulation", "structural constituent of cell wall", "reactive oxygen species metabolic", "peroxidase", oxidoreductase", and "antioxidant activity" were the common processes in rapeseed flooding response. In addition, the processes such as "hormone-mediated signaling pathway", "response to organic substance response", "motor activity", and "microtubule-based process" are likely to confer rapeseed flooding resistance. Mclust analysis clustered DEGs into nine modules; genes in each module shared similar expression patterns and many of these genes overlapped with the top 20 DEGs in some groups. This work provides a comprehensive insight into gene responses and the regulatory network in rapeseed flooding stress and provides guidelines for probing the underlying molecular mechanisms in flooding resistance.

Fine mapping of a kernel length-related gene with potential value for maize breeding.

KEY MESSAGE: A key candidate gene for maize kernel length was fine mapped to an interval of 942 kb; the locus significantly increases kernel length (KL) and hundred-kernel weight (HKW). Kernel size is a major determinant of yield in cereals. Kernel length, one of the determining factors of kernel size, is a target trait for both domestication and artificial breeding. However, there are few reports of fine mapping and quantitative trait loci (QTLs)/cloned genes for kernel length in maize. In this project, a novel major QTL, named qKL9, controlling maize kernel length was identified. We verified the authenticity and stability of qKL9 via BC2F2 and BC3F1 populations, respectively, and ultimately mapped qKL9 to an ~ 942-kb genomic interval by testing the progenies of recombination events derived from BC3F2 and BC4F2 populations in multiple environments. Additionally, one new line (McqKL9-A) containing the ~ 942-kb segment was screened from the BC4F2 population. Combining transcriptome analysis between McqKL9-A and Mc at 6, 9 and 14 days after pollination and candidate regional association mapping, Zm00001d046723 was preliminarily identified as the key candidate gene for qKL9. Importantly, the replacement in the Mc line of the Mc's alleles by the V671's alleles in the qKL9 region improved the performances of single-cross hybrids obtained with elite lines, illustrating the potential value of this QTL for the genetic improvement in maize kernel-related traits. These findings facilitate molecular breeding for kernel size and cloning of the gene underlying qKL9, shedding light on the genetic basis of kernel size in maize.

Genome-Wide Association Studies of Salt Tolerance at Seed Germination and Seedling Stages in Brassica napus.

Most crops are sensitive to salt stress, but their degree of susceptibility varies among species and cultivars. In order to understand the salt stress adaptability of Brassica napus to salt stress, we collected the phenotypic data of 505 B. napus accessions at the germination stage under 150 or 215 mM sodium chloride (NaCl) and at the seedling stage under 215 mM NaCl. Genome-wide association studies (GWAS) of 16 salt tolerance coefficients (STCs) were applied to investigate the genetic basis of salt stress tolerance of B. napus. In this study, we mapped 31 salts stress-related QTLs and identified 177 and 228 candidate genes related to salt stress tolerance were detected at germination and seedling stages, respectively. Overexpression of two candidate genes, BnCKX5 and BnERF3 overexpression, were found to increase the sensitivity to salt and mannitol stresses at the germination stage. This study demonstrated that it is a feasible method to dissect the genetic basis of salt stress tolerance at germination and seedling stages in B. napus by GWAS, which provides valuable loci for improving the salt stress tolerance of B. napus. Moreover, these candidate genes are rich genetic resources for the following exploration of molecular mechanisms in adaptation to salt stress in B. napus.

ZmSRL5 is involved in drought tolerance by maintaining cuticular wax structure in maize.

Cuticular wax is a natural barrier on terrestrial plant organs, which protects plants from damages caused by a variety of stresses. Here, we report the identification and functional characterization of a cuticular-wax-related gene, Zea mays L. SEMI-ROLLED LEAF 5 (ZmSRL5). The loss-of-function mutant srl5, which was created by a 3,745 bp insertion in the first intron that led to the premature transcript, exhibited abnormal wax crystal morphology and distribution, which, in turn, caused the pleiotropic phenotypes including increased chlorophyll leaching and water loss rate, decreased leaf temperature, sensitivity to drought, as well as semi-rolled mature leaves. However, total wax amounts showed no significant difference between wild type and semi-rolled leaf5 (srl5) mutant. The phenotype of srl5 was confirmed through the generation of two allelic mutants using CRISPR/Cas9. ZmSRL5 encodes a CASPARIAN-STRIP-MEMBRANE-DOMAIN-LIKE (CASPL) protein located in plasma membrane, and highly expressed in developing leaves. Further analysis showed that the expressions of the most wax related genes were not affected or slightly altered in srl5. This study, thus, primarily uncovers that ZmSRL5 is required for the structure formation of the cuticular wax and could increase the drought tolerance by maintaining the proper cuticular wax structure in maize.

A Comprehensive Transcriptomics Analysis Reveals Long Non-Coding RNA to be Involved in the Key Metabolic Pathway in Response to Waterlogging Stress in Maize.

Waterlogging stress (WS) in a dynamic environment seriously limits plant growth, development, and yield. The regulatory mechanism underlying WS conditions at an early stage in maize seedlings is largely unknown. In the present study, the primary root tips of B73 seedlings were sampled before (0 h) and after (2 h, 4 h, 6 h, 8 h, 10 h, and 12 h) WS and then subjected to transcriptome sequencing, resulting in the identification of differentially expressed protein-coding genes (DEpcGs) and long non-coding RNAs (DElncRs) in response to WS. These DEpcGs were classified into nine clusters, which were significantly enriched in several metabolic pathways, such as glycolysis and methionine metabolism. Several transcription factor families, including AP2-EREBP, bZIP, NAC, bHLH, and MYB, were also significantly enriched. In total, 6099 lncRNAs were identified, of which 3190 were DElncRs. A co-expression analysis revealed lncRNAs to be involved in 11 transcription modules, 10 of which were significantly associated with WS. The DEpcGs in the four modules were enriched in the hypoxia response pathways, including phenylpropanoid biosynthesis, MAPK signaling, and carotenoid biosynthesis, in which 137 DElncRs were also co-expressed. Most of the co-expressed DElncRs were co-localized with previously identified quantitative trait loci associated with waterlogging tolerance. A quantitative reverse transcription-polymerase chain reaction analysis of DEpcG and DElncR expression among the 32 maize genotypes after 4 h of WS verified significant expression correlations between them as well as significant correlation with the phenotype of waterlogging tolerance. Moreover, the high proportion of hypoxia response elements in the promoter region increased the reliability of the DElncRs identified in this study. These results provide a comprehensive transcriptome in response to WS at an early stage of maize seedlings and expand our understanding of the regulatory network involved in hypoxia in plants.

ZmSMK9, a pentatricopeptide repeat protein, is involved in the cis-splicing of nad5, kernel development and plant architecture in maize.

Maize kernel size and weight are essential contributors to its yield. So the identification of the genes controlling kernel size and weight can give us a chance to gain the yield. Here, we identified a small kernel mutant, Zea mays small kernel 9 (Zmsmk9), in maize. Cytological observation showed that the development of the endosperm and embryo was delayed in Zmsmk9 mutants at the early stages, resulting in a small kernel phenotype. Interestingly, despite substantial variation in kernel size, the germination of Zmsmk9 seeds was comparable to that of WT, and could develop into normal plants with upright leaf architecture. We cloned Zmsmk9 via map-based cloning. ZmSMK9 encodes a P-type pentatricopeptide repeat protein that targets to mitochondria, and is involved in RNA splicing in mitochondrial NADH dehydrogenase5 (nad5) intron-1 and intron-4. Consistent with the delayed development phenotype, transcriptome analysis of 12-DAP endosperm showed that starch and zeins biosynthesis related genes were dramatically down regulated in Zmsmk9, while cell cycle and cell growth related genes were dramatically increased. As a result, ZmSMK9 is a novel gene required for the splicing of nad5 intron-1 and intron-4, kernel development, and plant architecture in maize.

A Mitochondrial Transcription Termination Factor, ZmSmk3, Is Required for nad1 Intron4 and nad4 Intron1 Splicing and Kernel Development in Maize.

The expression systems of the mitochondrial genes are derived from their bacterial ancestors, but have evolved many new features in their eukaryotic hosts. Mitochondrial RNA splicing is a complex process regulated by families of nucleus-encoded RNA-binding proteins, few of which have been characterized in maize (Zea mays L.). Here, we identified the Zea mays small kernel 3 (Zmsmk3) candidate gene, which encodes a mitochondrial transcription termination factor (mTERF) containing two mTERF motifs, which is conserved in monocotyledon; and the target introns were also quite conserved during evolution between monocotyledons and dicotyledons. The mutations of Zmsmk3 led to arrested embryo and endosperm development, resulting in small kernels. A transcriptome of 12 days after pollination endosperm analysis revealed that the starch biosynthetic pathway and the zein gene family were down-regulated in the Zmsmk3 mutant kernels. ZmSMK3 is localized in mitochondria. The reduced expression of ZmSmk3 in the mutant resulted in the splicing deficiency of mitochondrial nad4 intron1 and nad1 intron4, causing a reduction in complex I assembly and activity, impairing mitochondria structure and activating the alternative respiratory pathway. So, the results suggest that ZmSMK3 is required for the splicing of nad4 intron 1 and nad1 intron 4, complex I assembly and kernel development in maize.