RNA polymerase common subunit ZmRPABC5b is transcriptionally activated by Opaque2 and essential for endosperm development in maize.

Maize (Zea mays) kernel size is an important factor determining grain yield; although numerous genes regulate kernel development, the roles of RNA polymerases in this process are largely unclear. Here, we characterized the defective kernel 701 (dek701) mutant that displays delayed endosperm development but normal vegetative growth and flowering transition, compared to its wild type. We cloned Dek701, which encoded ZmRPABC5b, a common subunit to RNA polymerases I, II and III. Loss-of-function mutation of Dek701 impaired the function of all three RNA polymerases and altered the transcription of genes related to RNA biosynthesis, phytohormone response and starch accumulation. Consistent with this observation, loss-of-function mutation of Dek701 affected cell proliferation and phytohormone homeostasis in maize endosperm. Dek701 was transcriptionally regulated in the endosperm by the transcription factor Opaque2 through binding to the GCN4 motif within the Dek701 promoter, which was subjected to strong artificial selection during maize domestication. Further investigation revealed that DEK701 interacts with the other common RNA polymerase subunit ZmRPABC2. The results of this study provide substantial insight into the Opaque2-ZmRPABC5b transcriptional regulatory network as a central hub for regulating endosperm development in maize.

Tetratricopeptide-containing SMALL KERNEL 11 is essential for the assembly of cytochrome c oxidase in maize mitochondria.

Assembly of the functional complexes of the mitochondrial respiratory chain requires sophisticated and efficient regulatory mechanisms. In plants, the subunit composition and assembly factors involved in the biogenesis of cytochrome c oxidase (complex IV) are substantially less defined than in mammals and yeast. In this study, we cloned maize (Zea mays) Small kernel 11 (Smk11) via map-based cloning. Smk11 encodes a mitochondria-localized tetratricopeptide repeat protein. Disruption of Smk11 severely affected the assembly and activity of mitochondrial complex IV, leading to delayed plant growth and seed development. Protein interactions studies revealed that SMK11 might interact with four putative complex IV assembly factors, Inner membrane peptidase 1A (ZmIMP1A), MYB domain protein 3R3 (ZmMYB3R-3), cytochrome c oxidase 23 (ZmCOX23), and mitochondrial ferredoxin 1 (ZmMFDX1), among which ZmMFDX1 might interact with subunits ZmCOX6a and ZmCOX-X1; ZmMYB3R-3 might also interact with ZmCOX6a. The mutation of SMK11 perturbed the normal assembly of these subunits, leading to the inactivation of complex IV. The results of this study revealed that SMK11 serves as an accessory assembly factor required for the normal assembly of subunits into complex IV, which will accelerate the elucidation of the assembly of complex IV in plant mitochondria.

cis-Regulatory variation affecting gene expression contributes to the improvement of maize kernel size.

cis-Regulatory variations contribute to trait evolution and adaptation during crop domestication and improvement. As the most important harvested organ in maize (Zea mays L.), kernel size has undergone intensive selection for size. However, the associations between maize kernel size and cis-regulatory variations remain unclear. We chose two independent association populations to dissect the genetic architecture of maize kernel size together with transcriptomic and genotypic data. The resulting phenotypes reflected a strong influence of population structure on kernel size. Compared with genome-wide association studies (GWASs), which accounted for population structure and relatedness, GWAS based on a naïve or simple linear model revealed additional associated single-nucleotide polymorphisms significantly involved in the conserved pathways controlling seed size in plants. Regulation analyses through expression quantitative trait locus mapping revealed that cis-regulatory variations likely control kernel size by fine-tuning the expression of proximal genes, among which ZmKL1 (GRMZM2G098305) was transgenically validated. We also proved that the pyramiding of the favorable cis-regulatory variations has contributed to the improvement of maize kernel size. Collectively, our results demonstrate that cis-regulatory variations, together with their regulatory genes, provide excellent targets for future maize improvement.

Combined population transcriptomic and genomic analysis reveals cis-regulatory differentiation of non-coding RNAs in maize.

KEY MESSAGE: Long intergenic non-coding RNA (lincRNA), cis-acting expression quantitative trait locus (cis-eQTL), maize, regulatory evolution. The law of genetic variation during domestication explains the evolutionary mechanism and provides a theoretical basis for improving existing varieties of maize. Previous studies focused on exploiting regulatory variations controlling the expression of protein-coding genes rather than of non-protein-coding genes. Here, we examined the genetic and evolutionary features of long non-coding RNAs from intergenic regions (long intergenic non-coding RNAs, lincRNAs) using population-scale transcriptome data and identified 1168 lincRNAs with cis-acting expression quantitative trait loci (cis-eQTLs). We found that lincRNAs are more likely to be regulated by cis-eQTLs, which exert stronger effects than the protein-coding genes. During maize domestication and improvement, upregulated alleles of lincRNAs, which originated from both standing variation and new mutation, accumulate more frequently and show larger effect sizes than the coding genes. A stronger signature of genetic differentiation was observed in their regulatory regions compared to those of randomly sampled lincRNAs. In addition, we found that cis-regulatory differentiation of lincRNAs is related to the sequence conservation of lincRNA transcripts. Non-conserved lincRNAs more tend to gain upregulated alleles and show a stronger relationship with selected traits than conserved lincRNAs between maize and its wild relatives. Our findings in maize improve the understanding of cis-regulatory variation in lincRNA genes during domestication and improvement and provide an effective approach for prioritizing candidates for further investigation.

Characterization and fine mapping of a maize lesion mimic mutant (Les8) with enhanced resistance to Curvularia leaf spot and southern leaf blight.

KEY MESSAGE: A novel light-dependent dominant lesion mimic mutant with enhanced multiple disease resistance was physiologically, biochemically, and genetically characterized; the causal gene was fine mapped to a 909 kb interval containing 38 genes. Identification of genes that confer multiple disease resistance (MDR) is crucial for the improvement of maize disease resistance. However, very limited genes are identified as MDR genes in maize. In this study, we characterized a dominant disease lesion mimics 8 (Les8) mutant that had chlorotic lesions on the leaves and showed enhanced resistance to both curvularia leaf spot and southern leaf blight. Major agronomic traits were not obviously altered, while decreased chlorophyll content was observed in the mutant, and the genetic effect of the Les8 mutation was stable in different genetic backgrounds. By BSR-seq analysis and map-based cloning, the LES8 gene was mapped into a 909 kb region containing 38 candidate genes on chromosome 9 wherein no lesion mimic or disease-resistance genes were previously reported. Using transcriptomics analysis, we found that genes involved in defense responses and secondary metabolite biosynthesis were enriched in the significantly up-regulated genes, while genes involved in photosynthesis and carbohydrate-related pathways were enriched in the significantly down-regulated genes in Les8. In addition, there was an overaccumulation of jasmonic acid and lignin but not salicylic acid in Les8. Taken together, this study revealed candidate genes and potential mechanism underlying Les8-conferred MDR in maize.

An ARF24-ZmArf2 module influences kernel size in different maize haplotypes.

Members of the ADP-ribosylation factor family, which are GTP-binding proteins, are involved in metabolite transport, cell division, and expansion. Although there has been a significant amount of research on small GTP-binding proteins, their roles and functions in regulating maize kernel size remain elusive. Here, we identified ZmArf2 as a maize ADP-ribosylation factor-like family member that is highly conserved during evolution. Maize zmarf2 mutants showed a characteristic smaller kernel size. Conversely, ZmArf2 overexpression increased maize kernel size. Furthermore, heterologous expression of ZmArf2 dramatically elevated Arabidopsis and yeast growth by promoting cell division. Using expression quantitative trait loci (eQTL) analysis, we determined that ZmArf2 expression levels in various lines were mainly associated with variation at the gene locus. The promoters of ZmArf2 genes could be divided into two types, pS and pL, that were significantly associated with both ZmArf2 expression levels and kernel size. In yeast-one-hybrid screening, maize Auxin Response Factor 24 (ARF24) is directly bound to the ZmArf2 promoter region and negatively regulated ZmArf2 expression. Notably, the pS and pL promoter types each contained an ARF24 binding element: an auxin response element (AuxRE) in pS and an auxin response region (AuxRR) in pL, respectively. ARF24 binding affinity to AuxRR was much higher compared with AuxRE. Overall, our results establish that the small G-protein ZmArf2 positively regulates maize kernel size and reveals the mechanism of its expression regulation.

Underlying mechanism of accelerated cell death and multiple disease resistance in a maize lethal leaf spot 1 allele.

Multiple disease resistance (MDR) in maize has attracted increasing attention. However, the interplay between cell death and metabolite changes and their contributions to MDR remains elusive in maize. In this study, we identified a mutant named as lesion mimic 30 (les30) that showed 'suicidal' lesion formation in the absence of disease and had enhanced resistance to the fungal pathogen Curvularia lunata. Using map-based cloning, we identified the causal gene encoding pheophorbide a oxidase (PAO), which is known to be involved in chlorophyll degradation and MDR, and is encoded by LETHAL LEAF SPOT1 (LLS1). LLS1 was found to be induced by both biotic and abiotic stresses. Transcriptomics analysis showed that genes involved in defense responses and secondary metabolite biosynthesis were mildly activated in leaves of the les30 mutant without lesions, whilst they were strongly activated in leaves with lesions. In addition, in les30 leaves with lesions, there was overaccumulation of defense-associated phytohormones including jasmonic acid and salicylic acid, and of phytoalexins including phenylpropanoids, lignin, and flavonoids, suggesting that their biosynthesis was activated in a lesion-dependent manner. Taken together, our study implies the existence of an interactive amplification loop of interrupted chlorophyll degradation, cell death, expression of defense-related genes, and metabolite changes that results in suicidal lesion formation and MDR, and this has the potential to be exploited by genetic manipulation to improve maize disease resistance.

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.

Natural variation of CT2 affects the embryo/kernel weight ratio in maize.

Embryo size is a critical trait determining not only grain yield but also the nutrition of the maize kernel. Up to the present, only a few genes have been characterized affecting the maize embryo/kernel ratio. Here, we identify 63 genes significantly associated with maize embryo/kernel weight ratio using a genome-wide association study (GWAS). The peak GWAS signal shows that the natural variation in Zea mays COMPACT PLANT2 (CT2), encoding the heterotrimeric G protein α subunit, is significantly associated with the Embryo/Kernel Weight Ratio (EKWR). Further analyses show that a missense mutation of CT2 increases its enzyme activity and associates with EKWR. The function of CT2 on affecting embryo/kernel weight ratio is further validated by the characterization of two ct2 mutants, for which EKWR is significantly decreased. Subsequently, the key downstream genes of CT2 are identified by combining the differential expression analysis (DEG) of the ct2 mutant and quantitative trait transcript analysis in the GWAS population. In addition, the allele frequency spectrum shows that CT2 was under selective pressure during maize domestication. This study provides important genetic insights into the natural variation of maize embryo/kernel weight ratio, which could be applied in future maize breeding programs to improve grain yield and nutritional content.

PNGSeqR: An R Package for Rapid Candidate Gene Selection through Pooled Next-Generation Sequencing.

Although bulked segregant analysis (BSA) has been used extensively in genetic mapping, user-friendly tools which can integrate current algorithms for researchers with no background in bioinformatics are scarce. To address this issue, we developed an R package, PNGSeqR, which takes single-nucleotide polymorphism (SNP) markers from next-generation sequencing (NGS) data in variant call format (VCF) as the input file, provides four BSA algorithms to indicate the magnitude of genome-wide signals, and rapidly defines the candidate region through the permutation test and fractile quantile. Users can choose the analysis methods according to their data and experimental design. In addition, it also supports differential expression gene analysis (DEG) and gene ontology analysis (GO) to prioritize the target gene. Once the analysis is completed, the plots can conveniently be exported.

Global profiling of N6 -methyladenosine methylation in maize callus induction.

Callus induction is a dedifferentiation process that accompanies a cell fate transition, and epigenetic regulation plays a crucial role in the process. N6 -methyladenosine (m6A) methylation is an important mechanism in post-transcriptional epigenetic regulation and functions in cell reprogramming. However, the function of m6A methylation during callus induction is still unknown. Here, we performed transcriptome-wide m6A-seq on immature maize embryos after culturing for 2, 4, or 8 days with or without the auxin analogue 2,4-D. A total of 26,794 unique m6A peaks were detected from 17,456 maize genes; and 2,338 specific, 2,4-D-induced m6A peaks (D-specific m6A) were detected only in embryos cultured with 2,4-D. Furthermore, a positive correlation between m6A methylation and mRNA abundance was discovered in the genes with D-specific m6A deposition, especially at the beginning of callus induction. Key genes involved in callus induction, i.e. BABY BOOM and LBD transcription factors, underwent m6A methylation, increasing their transcript levels, thus improving callus induction. These results revealed the importance of m6A methylation during the early stage of callus induction and provided new insights into the molecular mechanism of callus induction at an epitranscriptomic level.