Small kernel 1 encodes a pentatricopeptide repeat protein required for mitochondrial nad7 transcript editing and seed development in maize (Zea mays) and rice (Oryza sativa)

The maize mTERF18 regulates transcriptional termination of the mitochondrial nad6 gene and is essential for kernel development

Mitochondria are semi-autonomous organelle present in eukaryotic cells, containing their own genome and transcriptional machinery. However, their functions are intricately linked to proteins encoded by the nuclear genome. Mitochondrial transcription termination factors (mTERFs) are nucleic acid-binding proteins involved in RNA splicing and transcription termination within plant mitochondria and chloroplasts. Despite their recognized importance, the specific roles of mTERF proteins in maize remain largely unexplored. Here, we clone and functionally characterize maize mTERF18 gene. Our findings reveal that mTERF18 mutations lead to severely undifferentiated embryos, resulting in abortive phenotypes. Early kernel exhibits abnormal basal endosperm transfer layer and a significant reduction in both starch and protein accumulation in mterf18. We identify the mTERF18 gene through mapping-based cloning and validate this gene through allelic tests. mTERF18 is widely expressed across various maize tissues and encodes a highly conserved mitochondrial protein. Transcriptome data reveal that mTERF18 mutations disrupt transcriptional termination of the nad6 gene, leading to undetectable levels of Nad6 protein and reduced complex I assembly and activity. Furthermore, transmission electron microscopy observation of mterf18 endosperm uncover severe mitochondrial defects. Collectively, these findings highlight the critical role of mTERF18 in mitochondrial gene transcription termination and its pivotal impact on maize kernel development.

Review: Recent advances in unraveling the genetic architecture of kernel row number in maize

Kernel row number (KRN) is an important trait in maize that significantly impacts maize yield. The high heritability of KRN underscores its significance in maize breeding programs. In this review, we summarize recent advances in understanding the mechanisms underlying the formation, differentiation, and regulation of KRN in maize. Specifically, we have discussed gene mapping studies, functional validation of KRN-associated genes, and the application of gene editing techniques to KRN in maize. We summarized the various methods used to map and fine-map QTLs controlling KRN and provide an overview of the current status of cloned KRN-regulating genes. Despite the identification of many genes associated with KRN, the complexity of its regulation-arising from multiple loci and intricate gene interactions-remains a challenge. Balancing KRN with kernel number per row (KNR) and kernel weight is critical for optimizing yield while ensuring stability across different environments. Furthermore, we analyzed the influence of environmental factors on KRN, noting that despite its high heritability, environmental conditions can significantly affect this trait. Combining genotype-phenotype relationships with environmental data using big data and artificial intelligence could enhance maize breeding efficiency and accelerate genetic gains. This review emphasizes the importance of balancing traits, integrating environmental factors, and leveraging advanced technologies in maize breeding to achieve optimal yield and stress tolerance. Finally, we outlined future research perspectives aimed at developing high-yielding maize varieties through advances in KRN-related research.

Integrated analysis of BSA-seq and RNA-seq identified the candidate genes for seed weight in Brassica juncea

Introduction

Brassica juncea is a major oilseed crop of Brassica. The seed weight is one of yield components in oilseed Brassica crops. Research on the genetic mechanism of seed weight is not only directly related to the yield and economic value of Brassica juncea but also can provide a theory foundation for studying other Brassica crops.

Methods

To map the genes for seed weight, the parental and F2 extreme bulks derived were constructed from the cross between the heavy-seeded accession 7981 and the light-seeded one Sichuan yellow (SY) of B. juncea, and used in bulk segregant sequencing (BSA-seq). Meanwhile, RNA-sequencing (RNA-seq) was performed for both parents at six seed development stages.

Results

Our results showed that a total of thirty five SNPs were identified in thirty two genes located on chromosomes A02 and A10, while fifty eight InDels in fifty one genes located on A01, A03, A05, A07, A09, A10, B01, B02 and B04. The 7,679 differentially expressed genes were identified in developing seeds between the parents. Furthermore, integrated analysis of BSA-seq and RNA-seq data revealed a cluster of nine genes on chromosome A10 and one gene on chromosome A05 that are putative candidate genes controlling seed weight in B. juncea.

Discussion

This study provides a new reference for research on Brassica seed weight and lays a solid foundation for the examination of seed in other Brassica crops.

Genome-Wide Identification of Pentatricopeptide Repeat (PPR) Gene Family and Multi-Omics Analysis Provide New Insights Into the Albinism Mechanism of Kandelia obovata Propagule Leaves

Pentatricopeptide repeat (PPR) gene family constitutes one of the largest gene families in plants, which mainly participate in RNA editing and RNA splicing of organellar RNAs, thereby affecting the organellar development. Recently, some evidence elucidated the important roles of PPR proteins in the albino process of plant leaves. However, the functions of PPR genes in the woody mangrove species have not been investigated. In this study, using a typical true mangrove Kandelia obovata, we systematically identified 298 PPR genes and characterized their general features and physicochemical properties, including evolutionary relationships, the subcellular localization, PPR motif type, the number of introns and PPR motifs, and isoelectric point, and so forth. Furthermore, we combined genome-wide association studies (GWAS) and transcriptome analysis to identify the genetic architecture and potential PPR genes associated with propagule leaves colour variations of K. obovata. As a result, we prioritized 16 PPR genes related to the albino phenotype using different strategies, including differentially expressed genes analysis and genetic diversity analysis. Further analysis discovered two genes of interest, namely Maker00002998 (PLS-type) and Maker00003187 (P-type), which were differentially expressed genes and causal genes detected by GWAS analysis. Moreover, we successfully predicted downstream target chloroplast genes (rps14, rpoC1 and rpoC2) bound by Maker00002998 PPR proteins. The experimental verification of RNA editing sites of rps14, rpoC1, and rpoC2 in our previous study and the verification of interaction between Maker00002998 and rps14 transcript using in vitro RNA pull-down assays revealed that Maker00002998 PPR protein might be involved in the post-transcriptional process of chloroplast genes. Our result provides new insights into the roles of PPR genes in the albinism mechanism of K. obovata propagule leaves.

Seeing the unseen in characterizing RNA editome during rice endosperm development

Rice (Oryza sativa L.) endosperm is essential to provide nutrients for seed germination and determine grain yield. RNA editing, a post-transcriptional modification essential for plant development, unfortunately, is not fully characterized during rice endosperm development. Here, we perform systematic analyses to characterize RNA editome during rice endosperm development. We find that most editing sites are C-to-U CDS-recoding in mitochondria, leading to increased hydrophobic amino acids and changed structures of mitochondrial proteins. Comparative analysis of RNA editome reveals that CDS-recoding sites present higher editing frequencies with lower variabilities and their resultant recoded amino acids tend to exhibit stronger evolutionary conservation across many land plants. Furthermore, we classify mitochondrial genes into three groups, presenting distinct patterns in terms of CDS-recoding events. Besides, we conduct genome-wide screening to detect pentatricopeptide repeat (PPR) proteins and construct PPR-RNA binding profiles, yielding candidate PPR editing factors related to rice endosperm development. Taken together, our findings provide valuable insights for deciphering fundamental mechanisms of rice endosperm development underlying RNA editing machinery.

Assembly and comparative analysis of the first complete mitochondrial genome of Salix psammophila, a good windbreak and sand fixation shrub

Salix psammophila, commonly known as the sandlive willow, is a vital shrub species within the Salicaceae family, particularly significant for its ecological role in regions susceptible to desertification and sandy soils. In this study, we assembled the complete S. psammophila mitochondrial genome using Pacbio HiFi third-generation sequencing data. The genome was found to be a typical single circular structure, with a total length of 715,555 bp and a GC content of 44.89%. We annotated 33 unique protein-coding genes (PCGs), which included 24 core mitochondrial genes and 9 variable genes, as well as 18 tRNA genes (5 of which were multicopy genes) and 3 rRNA genes. Comparative analysis of the PCGs from the mitochondrial genomes of S. psammophila, Populus deltoides, Populus simonii, Salix wilsonii, and Salix suchowensis revealed that these genes are relatively conserved within the Salicaceae family, with variability primarily occurring in the ribosomal protein genes. The absence of the rps14, which encodes a ribosomal protein, may have played a role in the evolution of stress tolerance in Salicaceae plants. Additionally, we identified 232 SSRs, 19 tandem repeat sequences, and 236 dispersed repeat sequences in the S. psammophila mitochondrial genome, with palindromic and forward repeats being the most abundant. The longest palindromic repeat measured 260 bp, while the longest forward repeat was 86,068 bp. Furthermore, 324 potential RNA editing sites were discovered, all involving C-to-U edits, with the nad4 having the highest number of edits. These findings provide valuable insights into the phylogenetic and genetic research of Salicaceae plants.

Restorer of fertility like 30, encoding a mitochondrion-localized pentatricopeptide repeat protein, regulates wood formation in poplar

Nuclear-mitochondrial communication is crucial for plant growth, particularly in the context of cytoplasmic male sterility (CMS) repair mechanisms linked to mitochondrial genome mutations. The restorer of fertility-like (RFL) genes, known for their role in CMS restoration, remain largely unexplored in plant development. In this study, we focused on the evolutionary relationship of RFL family genes in poplar specifically within the dioecious Salicaceae plants. PtoRFL30 was identified to be preferentially expressed in stem vasculature, suggesting a distinct correlation with vascular cambium development. Transgenic poplar plants overexpressing PtoRFL30 exhibited a profound inhibition of vascular cambial activity and xylem development. Conversely, RNA interference-mediated knockdown of PtoRFL30 led to increased wood formation. Importantly, we revealed that PtoRFL30 plays a crucial role in maintaining mitochondrial functional homeostasis. Treatment with mitochondrial activity inhibitors delayed wood development in PtoRFL30-RNAi transgenic plants. Further investigations unveiled significant variations in auxin accumulation levels within vascular tissues of PtoRFL30-transgenic plants. Wood development anomalies resulting from PtoRFL30 overexpression and knockdown were rectified by NAA and NPA treatments, respectively. Our findings underscore the essential role of the PtoRFL30-mediated mitochondrion-auxin signaling module in wood formation, shedding light on the intricate nucleus-organelle communication during secondary vascular development.

Rice floury endosperm26 encoding a mitochondrial single-stranded DNA-binding protein is essential for RNA-splicing of mitochondrial genes and endosperm development

Endosperm, the major storage organ in cereal grains, determines the grain yield and quality. Mitochondria provide the energy for dry matter accumulation, in the endosperm development. Although mitochondrial single-stranded DNA-binding proteins (mtSSBs) play a canonical role in the maintenance of single-stranded mitochondrial DNA, their molecular functions in RNA processing and endosperm development remain obscure. Here, we report a defective rice endosperm mutant, floury endosperm26 (flo26), which develops abnormal starch grains in the endosperm. Map-based cloning and complementation experiments showed that FLO26 allele encodes a mitochondrial single-stranded DNA-binding protein, named as mtSSB1.1. Loss of function of mtSSB1.1 affects the transcriptional level of many mitochondrially-encoded genes and RNA splicing of nad1, a core component of respiratory chain complex I in mitochondria. As a result, dysfunctional mature nad1 led to dramatically decreased complex I activity, thereby reducing ATP production. Our results reveal that mtSSB1.1 plays an important role in the maintenance of mitochondrial function and endosperm development by stabilizing the splicing of mitochondrial RNA in rice.

Combined BSA-Seq and RNA-Seq to Identify Potential Genes Regulating Fruit Size in Bottle Gourd (Lagenaria siceraria L.)

Fruit size is a crucial agronomic trait in bottle gourd, impacting both yield and utility. Despite its significance, the regulatory mechanism governing fruit size in bottle gourd remains largely unknown. In this study, we used bottle gourd (small-fruited H28 and large-fruited H17) parent plants to measure the width and length of fruits at various developmental stages, revealing a single 'S' growth curve for fruit expansion. Paraffin section observations indicated that both cell number and size significantly influence bottle gourd fruit size. Through bulked segregant analysis and combined genotype-phenotype analysis, the candidate interval regulating fruit size was pinpointed to 17,747,353 bp-18,185,825 bp on chromosome 9, encompassing 0.44 Mb and including 44 genes. Parental fruits in the rapid expansion stage were subjected to RNA-seq, highlighting that differentially expressed genes were mainly enriched in pathways related to cell wall biosynthesis, sugar metabolism, and hormone signaling. Transcriptome and resequencing analysis, combined with gene function annotation, identified six genes within the localized region as potential regulators of fruit size. This study not only maps the candidate interval of genes influencing fruit size in bottle gourd through forward genetics, but also offers new insights into the potential molecular mechanisms underlying this trait through transcriptome analysis.

PPR proteins in plants: roles, mechanisms, and prospects for rice research

Pentatricopeptide repeat (PPR) proteins constitute one of the largest protein families in land plants, with over 300 members in various species. Nearly all PPR proteins are nuclear-encoded and targeted to the chloroplast and mitochondria, modulating organellar gene expression by participating in RNA metabolism, including mRNA stability, RNA editing, RNA splicing, and translation initiation. Organelle RNA metabolism significantly influences chloroplast and mitochondria functions, impacting plant photosynthesis, respiration, and environmental responses. Over the past decades, PPR proteins have emerged as a research focus in molecular biology due to their diverse roles throughout plant life. This review summarizes recent progress in understanding the roles and molecular mechanisms of PPR proteins, emphasizing their functions in fertility, abiotic and biotic stress, grain quality, and chloroplast development in rice. Furthermore, we discuss prospects for PPR family research in rice, aiming to provide a theoretical foundation for future investigations and applications.

DEK219 and HSF17 Collaboratively Regulate the Kernel Length in Maize

The kernel length is a crucial determinant of maize (Zea mays L.) yield; however, only a limited number of genes regulating kernel length have been validated, thus leaving our understanding of the mechanisms governing kernel length incomplete. We previously identified a maize kernel mutant, defective kernel219 (dek219), which encodes the DICER-LIKE1 protein that is essential for miRNA biogenesis. The present study revealed that dek219 consistently exhibits a stable phenotype characterized by a reduced kernel length. Further analysis indicated that dek219 may reduce the kernel length by inhibiting the expression of genes involved in regulating kernel length. By employing miRNA-target gene prediction, expression analysis, and correlation analysis, we successfully identified nine transcription factors that potentially participate in the regulation of kernel length under the control of DEK219. Among them, the upregulation fold change of HEAT SHOCK TRANSCRIPTION FACTOR17 (HSF17) expression was the highest, and the difference was most significant. The results of transient expression analysis and electrophoretic mobility shift assay (EMSA) indicated that HSF17 can inhibit the expression of DEFECTIVE ENDOSPERM18 (DE18), a gene involved in regulating kernel length. Furthermore, the hsf17 mutant exhibited a significant increase in kernel length, suggesting that HSF17 functions as a negative regulator of kernel length. The results of this study provide crucial evidence for further elucidating the molecular regulatory mechanism underlying maize kernel length and also offer valuable genetic resources for breeding high-yielding maize varieties.

Multiple factors interact in editing of PPR-E+-targeted sites in maize mitochondria and plastids

RNA cytidine-to-uridine editing is essential for plant organellar gene expression. Pentatricopeptide repeat (PPR)-E+ proteins have been proposed to bind to target sites and recruit the cytidine deaminase AtDYW2, facilitated by AtNUWA. Here we analyze the function of ZmNUWA, ZmDYW2A, and ZmDYW2B and their relationships with other editing factors in maize. The zmdyw2a and zmdyw2b single mutants are normal, but the zmdyw2a::zmdyw2b and zmnuwa mutants are severely arrested in seed development. ZmNUWA, ZmDYW2A, and ZmDYW2B are dual localized in mitochondria and plastids. Loss of ZmNUWA decreases the editing at 99 mitochondrial sites and 8 plastid sites. Surprisingly, loss of ZmDYW2A:ZmDYW2B affects almost the same set of sites targeted by PPR-E+ proteins. ZmNUWA interacts with ZmDYW2A and ZmDYW2B, suggesting that ZmNUWA recruits ZmDYW2A/2B in the editing of PPR-E+-targeted sites in maize. Further protein interaction analyses show that ZmNUWA and ZmDYW2A/2B interact with ZmMORF1, ZmMORF8, ZmMORF2, and ZmMORF9 and that ZmOZ1 interacts with ZmORRM1, ZmDYW2A, ZmDYW2B, ZmMORF8, and ZmMORF9. These results suggest that the maize mitochondrial PPR-E+ editosome contains PPR-E+, ZmDYW2A/2B, ZmNUWA, and ZmMORF1/8, whereas the plastid PPR-E+ editosome is composed of PPR-E+, ZmDYW2A/2B, ZmNUWA, ZmMORF2/8/9, ZmORRM1, and ZmOZ1.

Assembly and comparative analysis of the initial complete mitochondrial genome of Primulina hunanensis (Gesneriaceae): a cave-dwelling endangered plant

BACKGROUND

Primulina hunanensis, a troglobitic plant within the Primulina genus of Gesneriaceae family, exhibits robust resilience to arid conditions and holds great horticultural potential as an ornamental plant. The work of chloroplast genome (cpDNA) has been recently accomplished, however, the mitochondrial genome (mtDNA) that is crucial for plant evolution has not been reported.

RESULTS

In this study, we sequenced and assembled the P. hunanensis complete mtDNA, and elucidated its evolutionary and phylogenetic relationships. The assembled mtDNA spans 575,242 bp with 43.54% GC content, encompassing 60 genes, including 37 protein-coding genes (PCGs), 20 tRNA genes, and 3 rRNA genes. Notably, high number of repetitive sequences in the mtDNA and substantial sequence translocation from chloroplasts to mitochondria were observed. To determine the evolutionary and taxonomic positioning of P. hunanensis, a phylogenetic tree was constructed using mitochondrial PCGs from P. hunanensis and 32 other taxa. Furthermore, an exploration of PCGs relative synonymous codon usage, identification of RNA editing events, and an investigation of collinearity with closely related species were conducted.

CONCLUSIONS

This study reports the initial assembly and annotation of P. hunanensis mtDNA, contributing to the limited mtDNA repository for Gesneriaceae plants and advancing our understanding of their evolution for improved utilization and conservation.

Unraveling the genetics of heat tolerance in chickpea landraces (Cicer arietinum L.) using genome-wide association studies

Chickpea, being an important grain legume crop, is often confronted with the adverse effects of high temperatures at the reproductive stage of crop growth, drastically affecting yield and overall productivity. The current study deals with an extensive evaluation of chickpea genotypes, focusing on the traits associated with yield and their response to heat stress. Notably, we observed significant variations for these traits under both normal and high-temperature conditions, forming a robust basis for genetic research and breeding initiatives. Furthermore, the study revealed that yield-related traits exhibited high heritability, suggesting their potential suitability for marker-assisted selection. We carried out single-nucleotide polymorphism (SNP) genotyping using the genotyping-by-sequencing (GBS) method for a genome-wide association study (GWAS). Overall, 27 marker-trait associations (MTAs) linked to yield-related traits, among which we identified five common MTAs displaying pleiotropic effects after applying a stringent Bonferroni-corrected p-value threshold of <0.05 [-log10(p) > 4.95] using the BLINK (Bayesian-information and linkage-disequilibrium iteratively nested keyway) model. Through an in-depth in silico analysis of these markers against the CDC Frontier v1 reference genome, we discovered that the majority of the SNPs were located at or in proximity to gene-coding regions. We further explored candidate genes situated near these MTAs, shedding light on the molecular mechanisms governing heat stress tolerance and yield enhancement in chickpeas such as indole-3-acetic acid-amido synthetase GH3.1 with GH3 auxin-responsive promoter and pentatricopeptide repeat-containing protein, etc. The harvest index (HI) trait was associated with marker Ca3:37444451 encoding aspartic proteinase ortholog sequence of Oryza sativa subsp. japonica and Medicago truncatula, which is known for contributing to heat stress tolerance. These identified MTAs and associated candidate genes may serve as valuable assets for breeding programs dedicated to tailoring chickpea varieties resilient to heat stress and climate change.

PPR596 Is Required for nad2 Intron Splicing and Complex I Biogenesis in Arabidopsis

Mitochondria are essential organelles that generate energy via oxidative phosphorylation. Plant mitochondrial genome encodes some of the respiratory complex subunits, and these transcripts require accurate processing, including C-to-U RNA editing and intron splicing. Pentatricopeptide repeats (PPR) proteins are involved in various organellar RNA processing events. PPR596, a P-type PPR protein, was previously identified to function in the C-to-U editing of mitochondrial rps3 transcripts in Arabidopsis. Here, we demonstrate that PPR596 functions in the cis-splicing of nad2 intron 3 in mitochondria. Loss of the PPR596 function affects the editing at rps3eU1344SS, impairs nad2 intron 3 splicing and reduces the mitochondrial complex I's assembly and activity, while inducing alternative oxidase (AOX) gene expression. This defect in nad2 intron splicing provides a plausible explanation for the slow growth of the ppr595 mutants. Although a few P-type PPR proteins are involved in RNA C-to-U editing, our results suggest that the primary function of PPR596 is intron splicing.

DEFECTIVE KERNEL 56 functions in mitochondrial RNA editing and maize seed development

Proper seed development is essential for achieving grain production, successful seed germination, and seedling establishment in maize (Zea mays). In the past few decades, pentatricopeptide repeat (PPR) proteins have been proven to play an essential role in regulating the development of maize kernels through posttranscriptional RNA modification of mitochondrial genes. However, the underlying mechanisms remain largely unknown. Here, we characterized a mutant of DEFECTIVE KERNEL 56 (DEK56) with defective kernels that exhibited arrested development of both the embryo and endosperm. Accordingly, we isolated DEK56 through a map-based cloning strategy and found that it encoded an E subgroup PPR protein located in the mitochondria. Dysfunction of DEK56 resulted in altered cytidine (C)-to-uridine (U) editing efficiency at 48 editing sites across 21 mitochondrial transcripts. Notably, the editing efficiency of the maturase-related (matR)-1124 site was substantially reduced or abolished in the dek56 mutant. Furthermore, we found that the splicing efficiency of NADH dehydrogenase subunit 4 (nad4) Introns 1 and 3 was substantially reduced in dek56 kernels, which might be a consequence of the defective MatR function. Through a protein-protein interaction test, we hypothesized that DEK56 carries out its function by recruiting the PPR-DYW protein PPR motif, coiled-coil, and DYW domain-containing protein 1 (PCW1). This interaction is facilitated by Multiple Organellar RNA Editing Factors (ZmMORFs) and Glutamine-Rich Protein 23 (ZmGRP23). Based on these findings, we developed a working model of PPR-mediated mitochondrial processing that plays an essential role in the development of maize kernels. The present study will further broaden our understanding of PPR-mediated seed development and provide a theoretical basis for maize improvement.

Maize Dek407 Encodes the Nitrate Transporter 1.5 and Is Required for Kernel Development

The kernel serves as the storage organ and harvestable component of maize, and it plays a crucial role in determining crop yield and quality. Understanding the molecular and genetic mechanisms of kernel development is of considerable importance for maize production. In this study, we obtained a mutant, which we designated defective kernel 407 (dek407), through ethyl methanesulfonate mutagenesis. The dek407 mutant exhibited reduced kernel size and kernel weight, as well as delayed grain filling compared with those of the wild type. Positional cloning and an allelism test revealed that Dek407 encodes a nitrate transporter 1/peptide transporter family (NPF) protein and is the allele of miniature 2 (mn2) that was responsible for a poorly filled defective kernel phenotype. A transcriptome analysis of the developing kernels showed that the mutation of Dek407 altered the expression of phytohormone-related genes, especially those genes associated with indole-3-acetic acid synthesis and signaling. Phytohormone measurements and analysis indicated that the endogenous indole-3-acetic acid content was significantly reduced by 66% in the dek407 kernels, which may be the primary cause of the defective phenotype. We further demonstrated that natural variation in Dek407 is associated with kernel weight and kernel size. Therefore, Dek407 is a potential target gene for improvement of maize yield.

PHB3 Is Required for the Assembly and Activity of Mitochondrial ATP Synthase in Arabidopsis

Mitochondrial ATP synthase is a multiprotein complex, which consists of a matrix-localized F1 domain (F1-ATPase) and an inner membrane-embedded Fo domain (Fo-ATPase). The assembly process of mitochondrial ATP synthase is complex and requires the function of many assembly factors. Although extensive studies on mitochondrial ATP synthase assembly have been conducted on yeast, much less study has been performed on plants. Here, we revealed the function of Arabidopsis prohibitin 3 (PHB3) in mitochondrial ATP synthase assembly by characterizing the phb3 mutant. The blue native PAGE (BN-PAGE) and in-gel activity staining assays showed that the activities of ATP synthase and F1-ATPase were significantly decreased in the phb3 mutant. The absence of PHB3 resulted in the accumulation of the Fo-ATPase and F1-ATPase intermediates, whereas the abundance of the Fo-ATPase subunit a was decreased in the ATP synthase monomer. Furthermore, we showed that PHB3 could interact with the F1-ATPase subunits β and δ in the yeast two-hybrid system (Y2H) and luciferase complementation imaging (LCI) assay and with Fo-ATPase subunit c in the LCI assay. These results indicate that PHB3 acts as an assembly factor required for the assembly and activity of mitochondrial ATP synthase.

Genetic Localization and Homologous Genes Mining for Barley Grain Size

Grain size is an important agronomic trait determining barley yield and quality. An increasing number of QTLs (quantitative trait loci) for grain size have been reported due to the improvement in genome sequencing and mapping. Elucidating the molecular mechanisms underpinning barley grain size is vital for producing elite cultivars and accelerating breeding processes. In this review, we summarize the achievements in the molecular mapping of barley grain size over the past two decades, highlighting the results of QTL linkage analysis and genome-wide association studies. We discuss the QTL hotspots and predict candidate genes in detail. Moreover, reported homologs that determine the seed size clustered into several signaling pathways in model plants are also listed, providing the theoretical basis for mining genetic resources and regulatory networks of barley grain size.

Rice FLOURY ENDOSPERM22, encoding a pentatricopeptide repeat protein, is involved in both mitochondrial RNA splicing and editing and is crucial for endosperm development

Most of the reported P-type pentatricopeptide repeat (PPR) proteins play roles in organelle RNA stabilization and splicing. However, P-type PPRs involved in both RNA splicing and editing have rarely been reported, and their underlying mechanism remains largely unknown. Here, we report a rice floury endosperm22 (flo22) mutant with delayed amyloplast development in endosperm cells. Map-based cloning and complementation tests demonstrated that FLO22 encodes a mitochondrion-localized P-type PPR protein. Mutation of FLO22 resulting in defective trans-splicing of mitochondrial nad1 intron 1 and perhaps causing instability of mature transcripts affected assembly and activity of complex Ⅰ, and mitochondrial morphology and function. RNA-seq analysis showed that expression levels of many genes involved in starch and sucrose metabolism were significantly down-regulated in the flo22 mutant compared with the wild type, whereas genes related to oxidative phosphorylation and the tricarboxylic acid cycle were significantly up-regulated. In addition to involvement in splicing as a P-type PPR protein, we found that FLO22 interacted with DYW3, a DYW-type PPR protein, and they may function synergistically in mitochondrial RNA editing. The present work indicated that FLO22 plays an important role in endosperm development and plant growth by participating in nad1 maturation and multi-site editing of mitochondrial messager RNA.

Characterization and Transcriptome Analysis of Maize Small-Kernel Mutant smk7a in Different Development Stages

The kernel serves as a storage organ for various nutrients and determines the yield and quality of maize. Understanding the mechanisms regulating kernel development is important for maize production. In this study, a small-kernel mutant smk7a of maize was characterized. Cytological observation suggested that the development of the endosperm and embryo was arrested in smk7a in the early development stage. Biochemical tests revealed that the starch, zein protein, and indole-3-acetic acid (IAA) contents were significantly lower in smk7a compared with wild-type (WT). Consistent with the defective development phenotype, transcriptome analysis of the kernels 12 and 20 days after pollination (DAP) revealed that the starch, zein, and auxin biosynthesis-related genes were dramatically downregulated in smk7a. Genetic mapping indicated that the mutant was controlled by a recessive gene located on chromosome 2. Our results suggest that disrupted nutrition accumulation and auxin synthesis cause the defective endosperm and embryo development of smk7a.

A Systemic Investigation of Genetic Architecture and Gene Resources Controlling Kernel Size-Related Traits in Maize

Grain yield is the most critical and complex quantitative trait in maize. Kernel length (KL), kernel width (KW), kernel thickness (KT) and hundred-kernel weight (HKW) associated with kernel size are essential components of yield-related traits in maize. With the extensive use of quantitative trait locus (QTL) mapping and genome-wide association study (GWAS) analyses, thousands of QTLs and quantitative trait nucleotides (QTNs) have been discovered for controlling these traits. However, only some of them have been cloned and successfully utilized in breeding programs. In this study, we exhaustively collected reported genes, QTLs and QTNs associated with the four traits, performed cluster identification of QTLs and QTNs, then combined QTL and QTN clusters to detect consensus hotspot regions. In total, 31 hotspots were identified for kernel size-related traits. Their candidate genes were predicted to be related to well-known pathways regulating the kernel developmental process. The identified hotspots can be further explored for fine mapping and candidate gene validation. Finally, we provided a strategy for high yield and quality maize. This study will not only facilitate causal genes cloning, but also guide the breeding practice for maize.

Maize PPR-E proteins mediate RNA C-to-U editing in mitochondria by recruiting the trans deaminase PCW1

RNA C-to-U editing in organelles is essential for plant growth and development; however, the underlying mechanism is not fully understood. Here, we report that pentatricopeptide repeat (PPR)-E subclass proteins carry out RNA C-to-U editing by recruiting the trans deaminase PPR motifs, coiled-coil, and DYW domain-containing protein 1 (PCW1) in maize (Zea mays) mitochondria. Loss-of-function of bZIP and coiled-coil domain-containing PPR 1 (bCCP1) or PCW1 arrests seed development in maize. bCCP1 encodes a bZIP and coiled-coil domain-containing PPR protein, and PCW1 encodes an atypical PPR-DYW protein. bCCP1 is required for editing at 66 sites in mitochondria and PCW1 is required for editing at 102 sites, including the 66 sites that require bCCP1. The PCW1-mediated editing sites are exclusively associated with PPR-E proteins. bCCP1 interacts with PCW1 and the PPR-E protein Empty pericarp7 (EMP7). Two multiple organellar RNA editing factor (MORF) proteins, ZmMORF1 and ZmMORF8, interact with PCW1, EMP7, and bCCP1. ZmMORF8 enhanced the EMP7-PCW1 interaction in a yeast three-hybrid assay. C-to-U editing at the ccmFN-1553 site in maize required EMP7, bCCP1, and PCW1. These results suggest that PPR-E proteins function in RNA editing by recruiting the trans deaminase PCW1 and bCCP1, and MORF1/8 assist this recruitment through protein-protein interactions.

Ensemble classification based feature selection: a case of identification on plant pentatricopeptide repeat proteins

In order to identify plant pentatricopeptide repeat (PPR) proteins, a framework of variable selection has been proposed. In fact, it is an effective feature selection strategy that focuses on the performance of classification. Random forest has been used as the classifier with certain variables automatically selected for discrimination between PPR functional and non-functional proteins. However, it is found that samples regarded as PPR functional proteins are wrongly classified in a high rate. In this paper, we plan to improve the framework in order to achieve better classification results. Modifications are made on the framework for better identifying PPR functional proteins. Instead of random forest, a hybrid ensemble classifier is built with its base classifiers derived from six different classification methods. Besides, an incremental strategy and a clustering by search in descending order are alternatively used for feature selection, which can effectively select the most representative variables for identification on PPR proteins. In addition, it can be found that different base classifiers alternately play an important role in the ensemble classifier with feature dimension increasing. The experimental results demonstrate the effectiveness of our improvements.

Rice Chalky Grain 5 regulates natural variation for grain quality under heat stress

Heat stress occurring during rice (Oryza sativa) grain development reduces grain quality, which often manifests as increased grain chalkiness. Although the impact of heat stress on grain yield is well-studied, the genetic basis of rice grain quality under heat stress is less explored as quantifying grain quality is less tractable than grain yield. To address this, we used an image-based colorimetric assay (Red, R; and Green, G) for genome-wide association analysis to identify genetic loci underlying the phenotypic variation in rice grains exposed to heat stress. We found the R to G pixel ratio (RG) derived from mature grain images to be effective in distinguishing chalky grains from translucent grains derived from control (28/24°C) and heat stressed (36/32°C) plants. Our analysis yielded a novel gene, rice Chalky Grain 5 (OsCG5) that regulates natural variation for grain chalkiness under heat stress. OsCG5 encodes a grain-specific, expressed protein of unknown function. Accessions with lower transcript abundance of OsCG5 exhibit higher chalkiness, which correlates with higher RG values under stress. These findings are supported by increased chalkiness of OsCG5 knock-out (KO) mutants relative to wildtype (WT) under heat stress. Grains from plants overexpressing OsCG5 are less chalky than KOs but comparable to WT under heat stress. Compared to WT and OE, KO mutants exhibit greater heat sensitivity for grain size and weight relative to controls. Collectively, these results show that the natural variation at OsCG5 may contribute towards rice grain quality under heat stress.

Genome-wide analysis suggests the potential role of lncRNAs during seed development and seed size/weight determination in chickpea

The integrated transcriptome data analyses suggested the plausible roles of lncRNAs during seed development in chickpea. The candidate lncRNAs associated with QTLs and those involved in miRNA-mediated seed size/weight determination in chickpea have been identified. Long non-coding RNAs (lncRNAs) are important regulators of various biological processes. Here, we identified lncRNAs at seven successive stages of seed development in small-seeded and large-seeded chickpea cultivars. In total, 4751 lncRNAs implicated in diverse biological processes were identified. Most of lncRNAs were conserved between the two cultivars, whereas only a few of them were conserved in other plants, suggesting their species-specificity. A large number of lncRNAs differentially expressed between the two chickpea cultivars associated with seed development-related processes were identified. The lncRNAs acting as precursors of miRNAs and those mimicking target protein-coding genes of miRNAs involved in seed size/weight determination, including HAIKU1, BIG SEEDS1, and SHB1, were also revealed. Further, lncRNAs located within seed size/weight associated quantitative trait loci were also detected. Overall, we present a comprehensive resource and identified candidate lncRNAs that may play important roles during seed development and seed size/weight determination in chickpea.

Identification of major genomic regions for soybean seed weight by genome-wide association study

The hundred-seed weight (HSW) is an important yield component and one of the principal breeding traits in soybean. More than 250 quantitative trait loci (QTL) for soybean HSW have been identified. However, most of them have a large genomic region or are environmentally sensitive, which provide limited information for improving the phenotype in marker-assisted selection (MAS) and identifying the candidate genes. Here, we utilized 281 soybean accessions with 58,112 single nucleotide polymorphisms (SNPs) to dissect the genetic basis of HSW in across years in the northern Shaanxi province of China through one single-locus (SL) and three multi-locus (ML) genome-wide association study (GWAS) models. As a result, one hundred and fifty-four SNPs were detected to be significantly associated with HSW in at least one environment via SL-GWAS model, and 27 of these 154 SNPs were detected in all (three) environments and located within 7 linkage disequilibrium (LD) block regions with the distance of each block ranged from 40 to 610 Kb. A total of 15 quantitative trait nucleotides (QTNs) were identified by three ML-GWAS models. Combined with the results of different GWAS models, the 7 LD block regions associated with HSW detected by SL-GWAS model could be verified directly or indirectly by the results of ML-GWAS models. Eleven candidate genes underlying the stable loci that may regulate seed weight in soybean were predicted. The significantly associated SNPs and the stable loci as well as predicted candidate genes may be of great importance for marker-assisted breeding, polymerization breeding, and gene discovery for HSW in soybean.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-022-01310-y.

Cereal Endosperms: Development and Storage Product Accumulation

The persistent triploid endosperms of cereal crops are the most important source of human food and animal feed. The development of cereal endosperms progresses through coenocytic nuclear division, cellularization, aleurone and starchy endosperm differentiation, and storage product accumulation. In the past few decades, the cell biological processes involved in endosperm formation in most cereals have been described. Molecular genetic studies performed in recent years led to the identification of the genes underlying endosperm differentiation, regulatory network governing storage product accumulation, and epigenetic mechanism underlying imprinted gene expression. In this article, we outline recent progress in this area and propose hypothetical models to illustrate machineries that control aleurone and starchy endosperm differentiation, sugar loading, and storage product accumulations. A future challenge in this area is to decipher the molecular mechanisms underlying coenocytic nuclear division, endosperm cellularization, and programmed cell death.

EMP80 mediates the C-to-U editing of nad7 and atp4 and interacts with ZmDYW2 in maize mitochondria

RNA C-to-U editing is important to the expression and function of organellar genes in plants. Although several families of proteins have been identified to participate in this process, the underlying mechanism is not fully understood. Here we report the function of EMP80 in the C-to-U editing at the nad7-769 and atp4-118 sites, and the potential recruitment of ZmDYW2 as a trans deaminase in maize (Zea mays) mitochondria. Loss of EMP80 function arrests embryogenesis and endosperm development in maize. EMP80 is a PPR-E+ protein localised to mitochondria. An absence of EMP80 abolishes the C-to-U RNA editing at nad7-769 and atp4-118 sites, resulting in a cysteine-to-arginine (Cys→Arg) change in Nad7 and Atp4 in the emp80 mutant. The amino acid change consequently reduces the assembly of complexes I and V, leading to an accumulation of the F₁ subcomplex of complex V. EMP80 was found to interact with atypical DYW-type PPR protein ZmDYW2, which interacts with ZmNUWA. Co-expression of ZmNUWA enhances the interaction between EMP80 and ZmDYW2, suggesting that EMP80 potentially recruits ZmDYW2 as a trans deaminase through protein-protein interaction, and ZmNUWA may function as an enhancer of this interaction.

Genome-wide association studies dissect the genetic architecture of seed shape and size in common bean

Seed weight and size are important yield components. Thus, selecting for large seeds has been a key objective in crop domestication and breeding. In common bean, seed shape is also important since it influences industrial processing and plays a vital role in determining the choices of consumers and farmers. In this study, we performed genome-wide association studies on a core collection of common bean accessions to dissect the genetic architecture and identify genomic regions associated with seed morphological traits related to weight, size, and shape. Phenotypic data were collected by high-throughput image-based approaches, and utilized to test associations with 10,362 single-nucleotide polymorphism markers using multilocus mixed models. We searched within genome-associated regions for candidate genes putatively involved in seed phenotypic variation. The collection exhibited high variability for the entire set of seed traits, and the Andean gene pool was found to produce larger, heavier seeds than the Mesoamerican gene pool. Strong pairwise correlations were verified for most seed traits. Genome-wide association studies identified marker-trait associations accounting for a considerable amount of phenotypic variation in length, width, projected area, perimeter, and circularity in 4 distinct genomic regions. Promising candidate genes were identified, e.g. those encoding an AT-hook motif nuclear-localized protein 8, type 2C protein phosphatases, and a protein Mei2-like 4 isoform, known to be associated with seed size and weight regulation. Moreover, the genes that were pinpointed are also good candidates for functional analysis to validate their influence on seed shape and size in common bean and other related crops.

DEK48 Is Required for RNA Editing at Multiple Mitochondrial Sites and Seed Development in Maize

In flowering plants, C-to-U RNA editing can be critical to normal functions of mitochondrion-encoded proteins. Mitochondrial C-to-U RNA editing is facilitated by many factors from diverse protein families, of which the pentatricopeptide repeat (PPR) proteins play an important role. Owing to their large number and frequent embryo lethality in mutants, functions of many PPRs remain unknown. In this study, we characterized a mitochondrion-localized DYW-type PPR protein, DEK48, functioning in the C-to-U RNA editing at multiple mitochondrial transcripts in maize. Null mutation of Dek48 severely arrests embryo and endosperm development, causing a defective kernel (dek) phenotype, named dek48. DEK48 loss of function abolishes the C-to-U editing at nad3-185, -215, and nad4-376, -977 sites and decreases the editing at 11 other sites, resulting in the alteration of the corresponding amino acids. Consequently, the absence of editing caused reduced assembly and activity of complex I in dek48. Interestingly, we identified a point mutation in dek48-3 causing a deletion of the Tryptophan (W) residue in the DYW motif that abolishes the editing function. In sum, this study reveals the function of DEK48 in the C-to-U editing in mitochondrial transcripts and seed development in maize, and it demonstrates a critical role of the W residue in the DYW triplet motif of DEK48 for the C-to-U editing function in vivo.

Maize PPR278 Functions in Mitochondrial RNA Splicing and Editing

Pentatricopeptide repeat (PPR) proteins are a large protein family in higher plants and play important roles during seed development. Most reported PPR proteins function in mitochondria. However, some PPR proteins localize to more than one organelle; functional characterization of these proteins remains limited in maize (Zea mays L.). Here, we cloned and analyzed the function of a P-subfamily PPR protein, PPR278. Loss-function of PPR278 led to a lower germination rate and other defects at the seedling stage, as well as smaller kernels compared to the wild type. PPR278 was expressed in all investigated tissues. Furthermore, we determined that PPR278 is involved in the splicing of two mitochondrial transcripts (nad2 intron 4 and nad5 introns 1 and 4), as well as RNA editing of C-to-U sites in 10 mitochondrial transcripts. PPR278 localized to the nucleus, implying that it may function as a transcriptional regulator during seed development. Our data indicate that PPR278 is involved in maize seed development via intron splicing and RNA editing in mitochondria and has potential regulatory roles in the nucleus.

Dek504 Encodes a Mitochondrion-Targeted E+-Type Pentatricopeptide Repeat Protein Essential for RNA Editing and Seed Development in Maize

In flowering plants, RNA editing is a post-transcriptional process that selectively deaminates cytidines (C) to uridines (U) in organellar transcripts. Pentatricopeptide repeat (PPR) proteins have been identified as site-specific recognition factors for RNA editing. Here, we report the map-based cloning and molecular characterization of the defective kernel mutant dek504 in maize. Loss of Dek504 function leads to delayed embryogenesis and endosperm development, which produce small and collapsed kernels. Dek504 encodes an E+-type PPR protein targeted to the mitochondria, which is required for RNA editing of mitochondrial NADH dehydrogenase 3 at the nad3-317 and nad3-44 sites. Biochemical analysis of mitochondrial protein complexes revealed a significant reduction in the mitochondrial NADH dehydrogenase complex I activity, indicating that the alteration of the amino acid sequence at nad3-44 and nad3-317 through RNA editing is essential for NAD3 function. Moreover, the amino acids are highly conserved in monocots and eudicots, whereas the events of C-to-U editing are not conserved in flowering plants. Thus, our results indicate that Dek504 is essential for RNA editing of nad3, which is critical for NAD3 function, mitochondrial complex I stability, and seed development in maize.

Genetic Architecture of Grain Yield-Related Traits in Sorghum and Maize

Grain size, grain number per panicle, and grain weight are crucial determinants of yield-related traits in cereals. Understanding the genetic basis of grain yield-related traits has been the main research object and nodal in crop science. Sorghum and maize, as very close C4 crops with high photosynthetic rates, stress tolerance and large biomass characteristics, are extensively used to produce food, feed, and biofuels worldwide. In this review, we comprehensively summarize a large number of quantitative trait loci (QTLs) associated with grain yield in sorghum and maize. We placed great emphasis on discussing 22 fine-mapped QTLs and 30 functionally characterized genes, which greatly hinders our deep understanding at the molecular mechanism level. This review provides a general overview of the comprehensive findings on grain yield QTLs and discusses the emerging trend in molecular marker-assisted breeding with these QTLs.

Identification and Allele Combination Analysis of Rice Grain Shape-Related Genes by Genome-Wide Association Study

Grain shape is an important agronomic character of rice, which affects the appearance, processing, and the edible quality. Screening and identifying more new genes associated with grain shape is beneficial to further understanding the genetic basis of grain shape and provides more gene resources for genetic breeding. This study has a natural population containing 623 indica rice cultivars. Genome-wide association studies/GWAS of several traits related to grain shape (grain length/GL, grain width/GW, grain length to width ratio/GLWR, grain circumferences/GC, and grain size/grain area/GS) were conducted by combining phenotypic data from four environments and the second-generation resequencing data, which have identified 39 important Quantitative trait locus/QTLs. We analyzed the 39 QTLs using three methods: gene-based association analysis, haplotype analysis, and functional annotation and identified three cloned genes (GS3, GW5, OsDER1) and seven new candidate genes in the candidate interval. At the same time, to effectively utilize the genes in the grain shape-related gene bank, we have also analyzed the allelic combinations of the three cloned genes. Finally, the extreme allele combination corresponding to each trait was found through statistical analysis. This study’s novel candidate genes and allele combinations will provide a valuable reference for future breeding work.

A Structural Perspective on the RNA Editing of Plant Respiratory Complexes

The last steps of respiration, a core energy-harvesting process, are carried out by a chain of multi-subunit complexes in the inner mitochondrial membrane. Several essential subunits of the respiratory complexes are RNA-edited in plants, frequently leading to changes in the encoded amino acids. While the impact of RNA editing is clear at the sequence and phenotypic levels, the underlying biochemical explanations for these effects have remained obscure. Here, we used the structures of plant respiratory complex I, complex III₂ and complex IV to analyze the impact of the amino acid changes of RNA editing in terms of their location and biochemical features. Through specific examples, we demonstrate how the structural information can explain the phenotypes of RNA-editing mutants. This work shows how the structural perspective can bridge the gap between sequence and phenotype and provides a framework for the continued analysis of RNA-editing mutants in plant mitochondria and, by extension, in chloroplasts.

Research Progress of PPR Proteins in RNA Editing, Stress Response, Plant Growth and Development

RNA editing is a posttranscriptional phenomenon that includes gene processing and modification at specific nucleotide sites. RNA editing mainly occurs in the genomes of mitochondria and chloroplasts in higher plants. In recent years, pentatricopeptide repeat (PPR) proteins, which may act as trans-acting factors of RNA editing have been identified, and the study of PPR proteins has become a research focus in molecular biology. The molecular functions of these proteins and their physiological roles throughout plant growth and development are widely studied. In this minireview, we summarize the current knowledge of the PPR family, hoping to provide some theoretical reference for future research and applications.

Functions of PPR Proteins in Plant Growth and Development

Pentatricopeptide repeat (PPR) proteins form a large protein family in land plants, with hundreds of different members in angiosperms. In the last decade, a number of studies have shown that PPR proteins are sequence-specific RNA-binding proteins involved in multiple aspects of plant organellar RNA processing, and perform numerous functions in plants throughout their life cycle. Recently, computational and structural studies have provided new insights into the working mechanisms of PPR proteins in RNA recognition and cytidine deamination. In this review, we summarized the research progress on the functions of PPR proteins in plant growth and development, with a particular focus on their effects on cytoplasmic male sterility, stress responses, and seed development. We also documented the molecular mechanisms of PPR proteins in mediating RNA processing in plant mitochondria and chloroplasts.

GmPGL2, Encoding a Pentatricopeptide Repeat Protein, Is Essential for Chloroplast RNA Editing and Biogenesis in Soybean

Chloroplast biogenesis and development are highly complex processes requiring interactions between plastids and nuclear genomic products. Pentatricopeptide repeat (PPR) proteins play an essential role in the development of chloroplasts; however, it remains unclear how RNA editing factors influence soybean development. In this study, a Glycine max pale green leaf 2 mutant (Gmpgl2) was identified with decreased chlorophyll contents. Genetic mapping revealed that a single-nucleotide deletion at position 1949 bp in the Glyma.05g132700 gene in the Gmpgl2 mutant, resulting in a truncated GmPGL2 protein. The nuclear-encoded GmPGL2 is a PLS-type PPR protein that localizes to the chloroplasts. The C-to-U editing efficiencies of rps16, rps18, ndhB, ndhD, ndhE, and ndhF were reduced in the Gmpgl2 mutant. RNA electrophoresis mobility shift assay (REMSA) analysis further revealed that GmPGL2 binds to the immediate upstream sequences at RNA editing sites of rps16 and ndhB in vitro, respectively. In addition, GmPGL2 was found to interact with GmMORF8, GmMORF9, and GmORRM6. These results suggest that GmPGL2 participates in C-to-U RNA editing via the formation of a complex RNA editosome in soybean chloroplasts.

PPR-DYW Protein EMP17 Is Required for Mitochondrial RNA Editing, Complex III Biogenesis, and Seed Development in Maize

The conversion of cytidines to uridines (C-to-U) at specific sites in mitochondrial and plastid transcripts is a post-transcriptional processing event that is important to the expression of organellar genes. Pentatricopeptide repeat (PPR) proteins are involved in this process. In this study, we report the function of a previously uncharacterized PPR-DYW protein, Empty pericarp17 (EMP17), in the C-to-U editing and kernel development in maize. EMP17 is targeted to mitochondria. The loss-function of EMP17 arrests maize kernel development, abolishes the editing at ccmF _C -799 and nad2-677 sites, and reduces the editing at ccmF _C -906 and -966 sites. The absence of editing causes amino acid residue changes in CcmF_C-267 (Ser to Pro) and Nad2-226 (Phe to Ser), respectively. As CcmF_C functions in cytochrome c (Cytc) maturation, the amount of Cytc and Cytc ₁ protein is drastically reduced in emp17, suggesting that the CcmF_C-267 (Ser to Pro) change impairs the CcmF_C function. As a result, the assembly of complex III is strikingly decreased in emp17. In contrast, the assembly of complex I appears less affected, suggesting that the Nad2-226 (Phe to Ser) change may have less impact on Nad2 function. Together, these results indicate that EMP17 is required for the C-to-U editing at several sites in mitochondrial transcripts, complex III biogenesis, and seed development in maize.

Deepred-Mt: Deep representation learning for predicting C-to-U RNA editing in plant mitochondria

In land plant mitochondria, C-to-U RNA editing converts cytidines into uridines at highly specific RNA positions called editing sites. This editing step is essential for the correct functioning of mitochondrial proteins. When using sequence homology information, edited positions can be computationally predicted with high precision. However, predictions based on the sequence contexts of such edited positions often result in lower precision, which is limiting further advances on novel genetic engineering techniques for RNA regulation. Here, a deep convolutional neural network called Deepred-Mt is proposed. It predicts C-to-U editing events based on the 40 nucleotides flanking a given cytidine. Unlike existing methods, Deepred-Mt was optimized by using editing extent information, novel strategies of data augmentation, and a large-scale training dataset, constructed with deep RNA sequencing data of 21 plant mitochondrial genomes. In comparison to predictive methods based on sequence homology, Deepred-Mt attains significantly better predictive performance, in terms of average precision as well as F1 score. In addition, our approach is able to recognize well-known sequence motifs linked to RNA editing, and shows that the local RNA structure surrounding editing sites may be a relevant factor regulating their editing. These results demonstrate that Deepred-Mt is an effective tool for predicting C-to-U RNA editing in plant mitochondria. Source code, datasets, and detailed use cases are freely available at https://github.com/aedera/deepredmt.

Construction of a high-density genetic map and QTL analysis for yield, yield components and agronomic traits in chickpea (Cicer arietinum L.)

Unravelling the genetic architecture underlying yield components and agronomic traits is important for enhancing crop productivity. Here, a recombinant inbred line (RIL) population, developed from ICC 4958 and DCP 92–3 cross, was used for constructing linkage map and QTL mapping analysis. The RIL population was genotyped using a high-throughput Axiom^®CicerSNP array, which enabled the development of a high-density genetic map consisting of 3,818 SNP markers and spanning a distance of 1064.14 cM. Analysis of phenotyping data for yield, yield components and agronomic traits measured across three years together with genetic mapping data led to the identification of 10 major-effect QTLs and six minor-effect QTLs explaining up to 59.70% phenotypic variance. The major-effect QTLs identified for 100-seed weight, and plant height possessed key genes, such as C3HC4 RING finger protein, pentatricopeptide repeat (PPR) protein, sugar transporter, leucine zipper protein and NADH dehydrogenase, amongst others. The gene ontology studies highlighted the role of these genes in regulating seed weight and plant height in crop plants. The identified genomic regions for yield, yield components, and agronomic traits, and the closely linked markers will help advance genetics research and breeding programs in chickpea.

A Major and Stable Quantitative Trait Locus qSS2 for Seed Size and Shape Traits in a Soybean RIL Population

Seed size and shape traits are important determinants of seed yield and appearance quality in soybean [Glycine max (L.) Merr.]. Understanding the genetic architecture of these traits is important to enable their genetic improvement through efficient and targeted selection in soybean breeding, and for the identification of underlying causal genes. To map seed size and shape traits in soybean, a recombinant inbred line (RIL) population developed from K099 (small seed size) × Fendou 16 (large seed size), was phenotyped in three growing seasons. A genetic map of the RIL population was developed using 1,485 genotyping by random amplicon sequencing-direct (GRAS-Di) and 177 SSR markers. Quantitative trait locus (QTL) mapping was conducted by inclusive composite interval mapping. As a result, 53 significant QTLs for seed size traits and 27 significant QTLs for seed shape traits were identified. Six of these QTLs (qSW8.1, qSW16.1, qSLW2.1, qSLT2.1, qSWT1.2, and qSWT4.3) were identified with LOD scores of 3.80-14.0 and R ² of 2.36%-39.49% in at least two growing seasons. Among the above significant QTLs, 24 QTLs were grouped into 11 QTL clusters, such as, three major QTLs (qSL2.3, qSLW2.1, and qSLT2.1) were clustered into a major QTL on Chr.02, named as qSS2. The effect of qSS2 was validated in a pair of near isogenic lines, and its candidate genes (Glyma.02G269400, Glyma.02G272100, Glyma.02G274900, Glyma.02G277200, and Glyma.02G277600) were mined. The results of this study will assist in the breeding programs aiming at improvement of seed size and shape traits in soybean.

Rice FLOURY SHRUNKEN ENDOSPERM 5 Encodes a Putative Plant Organelle RNA Recognition Protein that Is Required for cis-Splicing of Mitochondrial nad4 Intron 1

BACKGROUND

The sequences of several important mitochondrion-encoded genes involved in respiration in higher plants are interrupted by introns. Many nuclear-encoded factors are involved in splicing these introns, but the mechanisms underlying this splicing remain unknown.

RESULTS

We isolated and characterized a rice mutant named floury shrunken endosperm 5 (fse5). In addition to having floury shrunken endosperm, the fse5 seeds either failed to germinate or produced seedlings which grew slowly and died ultimately. Fse5 encodes a putative plant organelle RNA recognition (PORR) protein targeted to mitochondria. Mutation of Fse5 hindered the splicing of the first intron of nad4, which encodes an essential subunit of mitochondrial NADH dehydrogenase complex I. The assembly and NADH dehydrogenase activity of complex I were subsequently disrupted by this mutation, and the structure of the mitochondria was abnormal in the fse5 mutant. The FSE5 protein was shown to interact with mitochondrial intron splicing factor 68 (MISF68), which is also a splicing factor for nad4 intron 1 identified previously via yeast two-hybrid (Y2H) assays.

CONCLUSION

Fse5 which encodes a PORR domain-containing protein, is essential for the splicing of nad4 intron 1, and loss of Fse5 function affects seed development and seedling growth.

Nuclear-Encoded Maturase Protein 3 Is Required for the Splicing of Various Group II Introns in Mitochondria during Maize (Zea mays L.) Seed Development

Splicing of plant organellar group II introns from precursor-RNA transcripts requires the assistance of nuclear-encoded splicing factors. Maturase (nMAT) is one such factor, as its three homologs (nMAT1, 2 and 4) have been identified as being required for the splicing of various mitochondrial introns in Arabidopsis. However, the function of nMAT in maize (Zea mays L.) is unknown. In this study, we identified a seed development mutant, empty pericarp 2441 (emp2441) from maize, which showed severely arrested embryogenesis and endosperm development. Positional cloning and transgenic complementation assays revealed that Emp2441 encodes a maturase-related protein, ZmnMAT3. ZmnMAT3 is highly expressed during seed development and its protein locates to the mitochondria. The loss of function of ZmnMAT3 resulted in the reduced splicing efficiency of various mitochondrial group II introns, particularly of the trans-splicing of nad1 introns 1, 3 and 4, which consequently abolished the transcript of nad1 and severely impaired the assembly and activity of mitochondrial complex I. Moreover, the Zmnmat3 mutant showed defective mitochondrial structure and exhibited expression and activity of alternative oxidases. These results indicate that ZmnMAT3 is essential for mitochondrial complex I assembly during kernel development in maize.

Maize kernel development

Maize (Zea mays) is a leading cereal crop in the world. The maize kernel is the storage organ and the harvest portion of this crop and is closely related to its yield and quality. The development of maize kernel is initiated by the double fertilization event, leading to the formation of a diploid embryo and a triploid endosperm. The embryo and endosperm are then undergone independent developmental programs, resulting in a mature maize kernel which is comprised of a persistent endosperm, a large embryo, and a maternal pericarp. Due to the well-characterized morphogenesis and powerful genetics, maize kernel has long been an excellent model for the study of cereal kernel development. In recent years, with the release of the maize reference genome and the development of new genomic technologies, there has been an explosive expansion of new knowledge for maize kernel development. In this review, we overviewed recent progress in the study of maize kernel development, with an emphasis on genetic mapping of kernel traits, transcriptome analysis during kernel development, functional gene cloning of kernel mutants, and genetic engineering of kernel traits.

Association Mapping and Transcriptome Analysis Reveal the Genetic Architecture of Maize Kernel Size

Kernel length, kernel width, and kernel thickness are important traits affecting grain yield and product quality. Here, the genetic architecture of the three kernel size traits was dissected in an association panel of 309 maize inbred lines using four statistical methods. Forty-two significant single nucleotide polymorphisms (SNPs; p < 1.72E-05) and 70 genes for the three traits were identified under five environments. One and eight SNPs were co-detected in two environments and by at least two methods, respectively, and they explained 5.87-9.59% of the phenotypic variation. Comparing the transcriptomes of two inbred lines with contrasting seed size, three and eight genes identified in the association panel showed significantly differential expression between the two genotypes at 15 and 39 days after pollination, respectively. Ten and 17 genes identified by a genome-wide association study were significantly differentially expressed between the two development stages in the two genotypes. Combining environment-/method-stable SNPs and differential expression analysis, ribosomal protein L7, jasmonate-regulated gene 21, serine/threonine-protein kinase RUNKEL, AP2-EREBP-transcription factor 16, and Zm00001d035222 (cell wall protein IFF6-like) were important candidate genes for maize kernel size and development.

Maize Defective Kernel605 Encodes a Canonical DYW-Type PPR Protein that Edits a Conserved Site of nad1 and Is Essential for Seed Nutritional Quality

Pentatricopeptide repeat (PPR) proteins involved in mitochondrial RNA cytidine (C)-to-uridine (U) editing mostly result in stagnant embryo and endosperm development upon loss of function. However, less is known about PPRs that are involved in farinaceous endosperm formation and maize quality. Here, we cloned a maize DYW-type PPR Defective Kernel605 (Dek605). Mutation of Dek605 delayed seed and seedling development. Mitochondrial transcript analysis of dek605 revealed that loss of DEK605 impaired C-to-U editing at the nad1-608 site and fails to alter Ser203 to Phe203 in NAD1 (dehydrogenase complex I), disrupting complex I assembly and reducing NADH dehydrogenase activity. Meanwhile, complexes III and IV in the cytochrome pathway, as well as AOX2 in the alternative respiratory pathway, are dramatically increased. Interestingly, the dek605 mutation resulted in opaque endosperm and increased levels of the free amino acids alanine, aspartic acid and phenylalanine. The down- and upregulated genes mainly involved in stress response-related and seed dormancy-related pathways, respectively, were observed after transcriptome analysis of dek605 at 12 d after pollination. Collectively, these results indicate that Dek605 specifically affects the single nad1-608 site and is required for normal seed development and resulted in nutritional quality relevant amino acid accumulation.

The novel E-subgroup pentatricopeptide repeat protein DEK55 is responsible for RNA editing at multiple sites and for the splicing of nad1 and nad4 in maize

BACKGROUND

Pentatricopeptide repeat (PPR) proteins compose a large protein family whose members are involved in both RNA processing in organelles and plant growth. Previous reports have shown that E-subgroup PPR proteins are involved in RNA editing. However, the additional functions and roles of the E-subgroup PPR proteins are unknown.

RESULTS

In this study, we developed and identified a new maize kernel mutant with arrested embryo and endosperm development, i.e., defective kernel (dek) 55 (dek55). Genetic and molecular evidence suggested that the defective kernels resulted from a mononucleotide alteration (C to T) at + 449 bp within the open reading frame (ORF) of Zm00001d014471 (hereafter referred to as DEK55). DEK55 encodes an E-subgroup PPR protein within the mitochondria. Molecular analyses showed that the editing percentage of 24 RNA editing sites decreased and that of seven RNA editing sites increased in dek55 kernels, the sites of which were distributed across 14 mitochondrial gene transcripts. Moreover, the splicing efficiency of nad1 introns 1 and 4 and nad4 intron 1 significantly decreased in dek55 compared with the wild type (WT). These results indicate that DEK55 plays a crucial role in RNA editing at multiple sites as well as in the splicing of nad1 and nad4 introns. Mutation in the DEK55 gene led to the dysfunction of mitochondrial complex I. Moreover, yeast two-hybrid assays showed that DEK55 interacts with two multiple organellar RNA-editing factors (MORFs), i.e., ZmMORF1 (Zm00001d049043) and ZmMORF8 (Zm00001d048291).

CONCLUSIONS

Our results demonstrated that a mutation in the DEK55 gene affects the mitochondrial function essential for maize kernel development. Our results also provide novel insight into the molecular functions of E-subgroup PPR proteins involved in plant organellar RNA processing.

Maize pentatricopeptide repeat protein DEK53 is required for mitochondrial RNA editing at multiple sites and seed development

Pentatricopeptide repeat (PPR) proteins were identified as site-specific recognition factors for RNA editing in plant mitochondria and plastids. In this study, we characterized maize (Zea mays) kernel mutant defective kernel 53 (dek53), which has an embryo lethal and collapsed endosperm phenotype. Dek53 encodes an E-subgroup PPR protein, which possesses a short PLS repeat region of only seven repeats. Subcellular localization analysis indicated that DEK53 is localized in the mitochondrion. Strand- and transcript-specific RNA-seq analysis showed that the dek53 mutation affected C-to-U RNA editing at more than 60 mitochondrial C targets. Biochemical analysis of mitochondrial protein complexes revealed a significant reduction in the assembly of mitochondrial complex III in dek53. Transmission electron microscopic examination showed severe morphological defects of mitochondria in dek53 endosperm cells. In addition, yeast two-hybrid and luciferase complementation imaging assays indicated that DEK53 can interact with the mitochondrion-targeted non-PPR RNA editing factor ZmMORF1, suggesting that DEK53 might be a functional component of the organellar RNA editosome.