Dynamic transcriptome landscape of maize embryo and endosperm development

The biosynthesis of storage reserves and auxin is coordinated by a hierarchical regulatory network in maize endosperm

Grain filling in maize (Zea mays) is intricately linked to cell development, involving the regulation of genes responsible for the biosynthesis of storage reserves (starch, proteins, and lipids) and phytohormones. However, the regulatory network coordinating these biological functions remains unclear. In this study, we identified 1744 high-confidence target genes co-regulated by the transcription factors (TFs) ZmNAC128 and ZmNAC130 (ZmNAC128/130) through chromatin immunoprecipitation sequencing coupled with RNA-seq analysis in the zmnac128/130 loss-of-function mutants. We further constructed a hierarchical regulatory network using DNA affinity purification sequencing analysis of downstream TFs regulated by ZmNAC128/130. In addition to target genes involved in the biosynthesis of starch and zeins, we discovered novel target genes of ZmNAC128/130 involved in the biosynthesis of lipids and indole-3-acetic acid (IAA). Consistently, the number of oil bodies, as well as the contents of triacylglycerol, and IAA were significantly reduced in zmnac128/130. The hierarchical regulatory network centered by ZmNAC128/130 revealed a significant overlap between the direct target genes of ZmNAC128/130 and their downstream TFs, particularly in regulating the biosynthesis of storage reserves and IAA. Our results indicated that the biosynthesis of storage reserves and IAA is coordinated by a multi-TFs hierarchical regulatory network in maize endosperm.

Toward unveiling transcriptome dynamics and regulatory modules at the maternal/filial interface of developing maize kernel

The basal region of maize (Zea mays) kernels, which includes the pedicel, placenta-chalazal, and basal endosperm transfer layers, serves as the maternal/filial interface for nutrient transfer from the mother plant to the developing seed. However, transcriptome dynamics of this maternal/filial interface remain largely unexplored. To address this gap, we conducted high-temporal-resolution RNA sequencing of the basal and upper kernel regions between 4 and 32 days after pollination and deeply analyzed transcriptome dynamics of the maternal/filial interface. Utilizing 790 specifically and highly expressed genes in the basal region, we performed the gene ontology (GO) term and weighted gene co-expression network analyses. In the early-stage basal region, we identified five MADS-box transcription factors (TFs) as hubs. Their homologs have been demonstrated as pivotal regulators at the maternal/filial interface of rice or Arabidopsis, suggesting their potential roles in maize kernel development. In the filling-stage basal region, numerous GO terms associated with transcriptional regulation and transporters are significantly enriched. Furthermore, we investigated the molecular function of three hub TFs. Through genome-wide DNA affinity purification sequencing combined with promoter transactivation assays, we suggested that these three TFs act as regulators of 10 basal-specific transporter genes involved in the transfer of sugars, amino acids, and ions. This study provides insights into transcriptomic dynamic and regulatory modules of the maternal/filial interface. In the future, genetic investigation of these hub regulators must advance our understanding of maternal/filial interface development and function.

Discovery and characterization of dietary antigens in oral tolerance

Food antigens elicit immune tolerance through the action of regulatory T cells (Tregs) in the intestine. Although antigens that trigger common food allergies are known, the epitopes that mediate tolerance to most foods have not been described. Here, we identified murine T cell receptors specific for maize, wheat, and soy, and used expression cloning to de-orphan their cognate epitopes. All of the epitopes derive from seed storage proteins that are resistant to degradation and abundant in the edible portion of the plant. Multiple unrelated T cell clones were specific for an epitope at the C-terminus of 19 kDa alpha-zein, a protein from maize kernel. An MHC tetramer loaded with this antigen revealed that zein-specific T cells are predominantly Tregs localized to the intestine. These cells, which develop concurrently with weaning, constitute up to 2% of the peripheral Treg pool. Bulk and single-cell RNA sequencing revealed that these cells express higher levels of immunosuppressive markers and chemokines compared to other Tregs. These data suggest that immune tolerance to plant-derived foods is focused on a specific class of antigens with common features, and they reveal the functional properties of naturally occurring food-specific Tregs.

Establishing an optimized ATAC-seq protocol for the maize

The advent of next-generation sequencing in crop improvement offers unprecedented insights into the chromatin landscape closely linked to gene activity governing key traits in plant development and adaptation. Particularly in maize, its dynamic chromatin structure is found to collaborate with massive transcriptional variations across tissues and developmental stages, implying intricate regulatory mechanisms, which highlights the importance of integrating chromatin information into breeding strategies for precise gene controls. The depiction of maize chromatin architecture using Assay for Transposase Accessible Chromatin with high-throughput sequencing (ATAC-seq) provides great opportunities to investigate cis-regulatory elements, which is crucial for crop improvement. In this context, we developed an easy-to-implement ATAC-seq protocol for maize with fewer nuclei and simple equipment. We demonstrate a streamlined ATAC-seq protocol with four key steps for maize in which nuclei purification can be achieved without cell sorting and using only a standard bench-top centrifuge. Our protocol, coupled with the bioinformatic analysis, including validation by read length periodicity, key metrics, and correlation with transcript abundance, provides a precise and efficient assessment of the maize chromatin landscape. Beyond its application to maize, our testing design holds the potential to be applied to other crops or other tissues, especially for those with limited size and amount, establishing a robust foundation for chromatin structure studies in diverse crop species.

Central Roles of ZmNAC128 and ZmNAC130 in Nutrient Uptake and Storage during Maize Grain Filling

Grain filling is critical for determining yield and quality, raising the question of whether central coordinators exist to facilitate the uptake and storage of various substances from maternal to filial tissues. The duplicate NAC transcription factors ZmNAC128 and ZmNAC130 could potentially serve as central coordinators. By analyzing differentially expressed genes from zmnac128 zmnac130 mutants across different genetic backgrounds and growing years, we identified 243 highly and differentially expressed genes (hdEGs) as the core target genes. These 243 hdEGs were associated with storage metabolism and transporters. ZmNAC128 and ZmNAC130 play vital roles in storage metabolism, and this study revealed two additional starch metabolism-related genes, sugary enhancer1 and hexokinase1, as their direct targets. A key finding of this study was the inclusion of 17 transporter genes within the 243 hdEGs, with significant alterations in the levels of more than 10 elements/substances in mutant kernels. Among them, six out of the nine upregulated transporter genes were linked to the transport of heavy metals and metalloids (HMMs), which was consistent with the enrichment of cadmium, lead, and arsenic observed in mutant kernels. Interestingly, the levels of Mg and Zn, minerals important to biofortification efforts, were reduced in mutant kernels. In addition to their direct involvement in sugar transport, ZmNAC128 and ZmNAC130 also activate the expression of the endosperm-preferential nitrogen and phosphate transporters ZmNPF1.1 and ZmPHO1;2. This coordinated regulation limits the intake of HMMs, enhances biofortification, and facilitates the uptake and storage of essential nutrients.

Comparative transcriptomic analysis provides insights into the genetic networks regulating oil differential production in oil crops

BACKGROUND

Plants differ more than threefold in seed oil contents (SOCs). Soybean (Glycine max), cotton (Gossypium hirsutum), rapeseed (Brassica napus), and sesame (Sesamum indicum) are four important oil crops with markedly different SOCs and fatty acid compositions.

RESULTS

Compared to grain crops like maize and rice, expanded acyl-lipid metabolism genes and relatively higher expression levels of genes involved in seed oil synthesis (SOS) in the oil crops contributed to the oil accumulation in seeds. Here, we conducted comparative transcriptomics on oil crops with two different SOC materials. In common, DIHYDROLIPOAMIDE DEHYDROGENASE, STEAROYL-ACYL CARRIER PROTEIN DESATURASE, PHOSPHOLIPID:DIACYLGLYCEROL ACYLTRANSFERASE, and oil-body protein genes were both differentially expressed between the high- and low-oil materials of each crop. By comparing functional components of SOS networks, we found that the strong correlations between genes in "glycolysis/gluconeogenesis" and "fatty acid synthesis" were conserved in both grain and oil crops, with PYRUVATE KINASE being the common factor affecting starch and lipid accumulation. Network alignment also found a conserved clique among oil crops affecting seed oil accumulation, which has been validated in Arabidopsis. Differently, secondary and protein metabolism affected oil synthesis to different degrees in different crops, and high SOC was due to less competition of the same precursors. The comparison of Arabidopsis mutants and wild type showed that CINNAMYL ALCOHOL DEHYDROGENASE 9, the conserved regulator we identified, was a factor resulting in different relative contents of lignins to oil in seeds. The interconnection of lipids and proteins was common but in different ways among crops, which partly led to differential oil production.

CONCLUSIONS

This study goes beyond the observations made in studies of individual species to provide new insights into which genes and networks may be fundamental to seed oil accumulation from a multispecies perspective.

Sucrose-associated SnRK1a1-mediated phosphorylation of Opaque2 modulates endosperm filling in maize

During maize endosperm filling, sucrose not only serves as a source of carbon skeletons for storage-reserve synthesis but also acts as a stimulus to promote this process. However, the molecular mechanisms underlying sucrose and endosperm filling are poorly understood. In this study, we found that sucrose promotes the expression of endosperm-filling hub gene Opaque2 (O2), coordinating with storage-reserve accumulation. We showed that the protein kinase SnRK1a1 can attenuate O2-mediated transactivation, but sucrose can release this suppression. Biochemical assays revealed that SnRK1a1 phosphorylates O2 at serine 41 (S41), negatively affecting its protein stability and transactivation ability. We observed that mutation of SnRK1a1 results in larger seeds with increased kernel weight and storage reserves, while overexpression of SnRK1a1 causes the opposite effect. Overexpression of the native O2 (O2-OE), phospho-dead (O2-SA), and phospho-mimetic (O2-SD) variants all increased 100-kernel weight. Although O2-SA seeds exhibit smaller kernel size, they have higher accumulation of starch and proteins, resulting in larger vitreous endosperm and increased test weight. O2-SD seeds display larger kernel size but unchanged levels of storage reserves and test weight. O2-OE seeds show elevated kernel dimensions and nutrient storage, like a mixture of O2-SA and O2-SD seeds. Collectively, our study discovers a novel regulatory mechanism of maize endosperm filling. Identification of S41 as a SnRK1-mediated phosphorylation site in O2 offers a potential engineering target for enhancing storage-reserve accumulation and yield in maize.

Two ABCI family transporters, OsABCI15 and OsABCI16, are involved in grain-filling in rice

Seed development is critical for plant reproduction and crop yield, with panicle seed-setting rate, grain-filling, and grain weight being key seed characteristics for yield improvement. However, few genes are known to regulate grain filling. Here, we identify two adenosine triphosphate (ATP)-binding cassette (ABC)I-type transporter genes, OsABCI15 and OsABCI16, involved in rice grain-filling. Both genes are highly expressed in developing seeds, and their proteins are localized to the plasma membrane and cytosol. Interestingly, knockout of OsABCI15 and OsABCI16 results in a significant reduction in seed-setting rate, caused predominantly by the severe empty pericarp phenotype, which differs from the previously reported low seed-setting phenotype resulting from failed pollination. Further analysis indicates that OsABCI15 and OsABCI16 participate in ion homeostasis and likely export ions between filial tissues and maternal tissues during grain filling. Importantly, overexpression of OsABCI15 and OsABCI16 enhances the seed-setting rate and grain yield in transgenic plants and decreases ion accumulation in brown rice. Moreover, the OsABCI15/16 orthologues in maize exhibit a similar role in kernel development, as demonstrated by their disruption in transgenic maize. Therefore, our findings reveal the important roles of two ABC transporters in cereal grain filling, highlighting their value in crop yield improvement.

ZmELP1, an Elongator complex subunit, is required for the maintenance of histone acetylation and RNA Pol II phosphorylation in maize kernels

The Elongator complex was originally identified as an interactor of hyperphosphorylated RNA polymerase II (RNAPII) in yeast and has histone acetyltransferase (HAT) activity. However, the genome-wide regulatory roles of Elongator on transcriptional elongation and histone acetylation remain unclear. We characterized a maize miniature seed mutant, mn7 and map-based cloning revealed that Mn7 encodes one of the subunits of the Elongator complex, ZmELP1. ZmELP1 deficiency causes marked reductions in the kernel size and weight. Molecular analyses showed that ZmELP1 interacts with ZmELP3, which is required for H3K14 acetylation (H3K14ac), and Elongator complex subunits interact with RNA polymerase II (RNAPII) C-terminal domain (CTD). Genome-wide analyses indicated that loss of ZmELP1 leads to a significant decrease in the deposition of H3K14ac and the CTD of phosphorylated RNAPII on Ser2 (Ser2P). These chromatin changes positively correlate with global transcriptomic changes. ZmELP1 mutation alters the expression of genes involved in transcriptional regulation and kernel development. We also showed that the decrease of Ser2P depends on the deposition of Elongator complex-mediated H3K14ac. Taken together, our results reveal an important role of ZmELP1 in the H3K14ac-dependent transcriptional elongation, which is critical for kernel development.

Identification and characterization of ACR gene family in maize for salt stress tolerance

Background

Members of the ACR gene family are commonly involved in various physiological processes, including amino acid metabolism and stress responses. In recent decades, significant progress has been made in the study of ACR genes in plants. However, little is known about their characteristics and function in maize.

Methods

In this study, ACR genes were identified from the maize genome, and their molecular characteristics, gene structure, gene evolution, gene collinearity analysis, cis-acting elements were analyzed. qRT-PCR technology was used to verify the expression patterns of the ZmACR gene family in different tissues under salt stress. In addition, Ectopic expression technique of ZmACR5 in Arabidopsis thaliana was utilized to identify its role in response to salt stress.

Results

A total of 28 ZmACR genes were identified, and their molecular characteristics were extensively described. Two gene pairs arising from segmented replication events were detected in maize, and 18 collinear gene pairs were detected between maize and 3 other species. Through phylogenetic analysis, three subgroups were revealed, demonstrating distinct divergence between monocotyledonous and dicotyledonous plants. Analysis of ZmACR cis-acting elements revealed the optional involvement of ZmACR genes in light response, hormone response and stress resistance. Expression analysis of 8 ZmACR genes under salt treatment clearly revealed their role in the response to salt stress. Ectopic overexpression of ZmACR5 in Arabidopsis notably reduced salt tolerance compared to that of the wild type under salt treatment, suggesting that ZmACR5 has a negative role in the response to salt stress.

Conclusion

Taken together, these findings confirmed the involvement of ZmACR genes in regulating salt stress and contributed significantly to our understanding of the molecular function of ACR genes in maize, facilitating further research in this field.

Combining association with linkage mapping to dissect the phenolamides metabolism of the maize kernel

Phenolamides are important secondary metabolites in plant species. They play important roles in plant defense responses against pathogens and insect herbivores, protection against UV irradiation and floral induction and development. However, the accumulation and variation in phenolamides content in diverse maize lines and the genes responsible for their biosynthesis remain largely unknown. Here, we combined genetic mapping, protein regulatory network and bioinformatics analysis to further enhance the understanding of maize phenolamides biosynthesis. Sixteen phenolamides were identified in multiple populations, and they were all significantly correlated with one or several of 19 phenotypic traits. By linkage mapping, 58, 58, 39 and 67 QTLs, with an average of 3.9, 3.6, 3.6 and 4.2 QTLs for each trait were mapped in BBE1, BBE2, ZYE1 and ZYE2, explaining 9.47%, 10.78%, 9.51% and 11.40% phenotypic variation for each QTL on average, respectively. By GWAS, 39 and 36 significant loci were detected in two different environments, 3.3 and 2.8 loci for each trait, explaining 10.00% and 9.97% phenotypic variation for each locus on average, respectively. Totally, 58 unique candidate genes were identified, 31% of them encoding enzymes involved in amine and derivative metabolic processes. Gene Ontology term analysis of the 358 protein-protein interrelated genes revealed significant enrichment in terms relating to cellular nitrogen metabolism, amine metabolism. GRMZM2G066142, GRMZM2G066049, GRMZM2G165390 and GRMZM2G159587 were further validated involvement in phenolamides biosynthesis. Our results provide insights into the genetic basis of phenolamides biosynthesis in maize kernels, understanding phenolamides biosynthesis and its nutritional content and ability to withstand biotic and abiotic stress.

Dynamic transcriptome landscape of maize pericarp development

As a maternal tissue, the pericarp supports and protects for other components of seed, such as embryo and endosperm. Despite the importance of maize pericarp in seed, the genome-wide transcriptome pattern throughout maize pericarp development has not been well characterized. Here, we developed RNA-seq transcriptome atlas of B73 maize pericarp development based on 21 samples from 5 days before fertilization (DBP5) to 32 days after fertilization (DAP32). A total of 25 346 genes were detected in programming pericarp development, including 1887 transcription factors (TFs). Together with pericarp morphological changes, the global clustering of gene expression revealed four developmental stages: undeveloped, thickening, expansion and strengthening. Coexpression analysis provided further insights on key regulators in functional transition of four developmental stages. Combined with non-seed, embryo, endosperm, and nucellus transcriptome data, we identified 598 pericarp-specific genes, including 75 TFs, which could elucidate key mechanisms and regulatory networks of pericarp development. Cell wall related genes were identified that reflected their crucial role in the maize pericarp structure building. In addition, key maternal proteases or TFs related with programmed cell death (PCD) were proposed, suggesting PCD in the maize pericarp was mediated by vacuolar processing enzymes (VPE), and jasmonic acid (JA) and ethylene-related pathways. The dynamic transcriptome atlas provides a valuable resource for unraveling the genetic control of maize pericarp development.

Defective kernel 58 encodes an Rrp15p domain-containing protein essential to ribosome biogenesis and seed development in maize

Ribosome biogenesis is a highly dynamic and orchestrated process facilitated by hundreds of ribosomal biogenesis factors and small nucleolar RNAs. While many of the advances are derived from studies in yeast, ribosome biogenesis remains largely unknown in plants despite its importance to plant growth and development. Through characterizing the maize (Zea mays) defective kernel and embryo-lethal mutant dek58, we show that DEK58 encodes an Rrp15p domain-containing protein with 15.3% identity to yeast Rrp15. Over-expression of DEK58 rescues the mutant phenotype. DEK58 is localized in the nucleolus. Ribosome profiling and RNA gel blot analyses show that the absence of DEK58 reduces ribosome assembly and impedes pre-rRNA processing, accompanied by the accumulation of nearly all the pre-rRNA processing intermediates and the production of an aberrant processing product P-25S*. DEK58 interacts with ZmSSF1, a maize homolog of the yeast Ssf1 in the 60S processome. DEK58 and ZmSSF1 interact with ZmCK2α, a putative component of the yeast UTP-C complex involved in the small ribosomal subunit processome. These results demonstrate that DEK58 is essential to seed development in maize. It functions in the early stage of pre-rRNA processing in ribosome biogenesis, possibly through interacting with ZmSSF1 and ZmCK2α in maize.

Multi-omics-driven advances in the understanding of triacylglycerol biosynthesis in oil seeds

Vegetable oils are rich sources of polyunsaturated fatty acids and energy as well as valuable sources of human food, animal feed, and bioenergy. Triacylglycerols, which are comprised of three fatty acids attached to a glycerol backbone, are the main component of vegetable oils. Here, we review the development and application of multiple-level omics in major oilseeds and emphasize the progress in the analysis of the biological roles of key genes underlying seed oil content and quality in major oilseeds. Finally, we discuss future research directions in functional genomics research based on current omics and oil metabolic engineering strategies that aim to enhance seed oil content and quality, and specific fatty acids components according to either human health needs or industrial requirements.

Decoding the gene regulatory network of endosperm differentiation in maize

The persistent cereal endosperm constitutes the majority of the grain volume. Dissecting the gene regulatory network underlying cereal endosperm development will facilitate yield and quality improvement of cereal crops. Here, we use single-cell transcriptomics to analyze the developing maize (Zea mays) endosperm during cell differentiation. After obtaining transcriptomic data from 17,022 single cells, we identify 12 cell clusters corresponding to five endosperm cell types and revealing complex transcriptional heterogeneity. We delineate the temporal gene-expression pattern from 6 to 7 days after pollination. We profile the genomic DNA-binding sites of 161 transcription factors differentially expressed between cell clusters and constructed a gene regulatory network by combining the single-cell transcriptomic data with the direct DNA-binding profiles, identifying 181 regulons containing genes encoding transcription factors along with their high-confidence targets, Furthermore, we map the regulons to endosperm cell clusters, identify cell-cluster-specific essential regulators, and experimentally validated three predicted key regulators. This study provides a framework for understanding cereal endosperm development and function at single-cell resolution.

Dynamic transcriptome landscape of developing maize ear

Seed number and harvesting ability in maize (Zea mays L.) are primarily determined by the architecture of female inflorescence, namely the ear. Therefore, ear morphogenesis contributes to grain yield and as such is one of the key target traits during maize breeding. However, the molecular networks of this highly dynamic and complex grain-bearing inflorescence remain largely unclear. As a first step toward characterizing these networks, we performed a high-spatio-temporal-resolution investigation of transcriptomes using 130 ear samples collected from developing ears with length from 0.1 mm to 19.0 cm. Comparisons of these mRNA populations indicated that these spatio-temporal transcriptomes were clearly separated into four distinct stages stages I, II, III, and IV. A total of 23 793 genes including 1513 transcription factors (TFs) were identified in the investigated developing ears. During the stage I of ear morphogenesis, 425 genes were predicted to be involved in a co-expression network established by eight hub TFs. Moreover, 9714 ear-specific genes were identified in the seven kinds of meristems. Additionally, 527 genes including 59 TFs were identified as especially expressed in ear and displayed high temporal specificity. These results provide a high-resolution atlas of gene activity during ear development and help to unravel the regulatory modules associated with the differentiation of the ear in maize.

Spatial transcriptomics uncover sucrose post-phloem transport during maize kernel development

Maize kernels are complex biological systems composed of three genetic sources, namely maternal tissues, progeny embryos, and progeny endosperms. The lack of gene expression profiles with spatial information has limited the understanding of the specific functions of each cell population, and hindered the exploration of superior genes in kernels. In our study, we conduct microscopic sectioning and spatial transcriptomics analysis during the grain filling stage of maize kernels. This enables us to visualize the expression patterns of all genes through electronical RNA in situ hybridization, and identify 11 cell populations and 332 molecular marker genes. Furthermore, we systematically elucidate the spatial storage mechanisms of the three major substances in maize kernels: starch, protein, and oil. These findings provide valuable insights into the functional genes that control agronomic traits in maize kernels.

Endoreduplication in plant organogenesis: a means to boost fruit growth

Endoreduplication is the major source of somatic endopolyploidy in higher plants, and leads to variation in cell ploidy levels due to iterative rounds of DNA synthesis in the absence of mitosis. Despite its ubiquitous occurrence in many plant organs, tissues, and cells, the physiological meaning of endoreduplication is not fully understood, although several roles during plant development have been proposed, mostly related to cell growth, differentiation, and specialization via transcriptional and metabolic reprogramming. Here, we review recent advances in our knowledge of the molecular mechanisms and cellular characteristics of endoreduplicated cells, and provide an overview of the multi-scale effects of endoreduplication on supporting growth in plant development. In addition, the effects of endoreduplication in fruit development are discussed, since it is highly prominent during fruit organogenesis where it acts as a morphogenetic factor supporting rapid fruit growth, as illustrated by case of the model fleshy fruit, tomato (Solanum lycopersicum).

The transcription factors ZmNAC128 and ZmNAC130 coordinate with Opaque2 to promote endosperm filling in maize

Endosperm filling in maize (Zea mays), which involves nutrient uptake and biosynthesis of storage reserves, largely determines grain yield and quality. However, much remains unclear about the synchronization of these processes. Here, we comprehensively investigated the functions of duplicate NAM, ATAF1/2, and CUC2 (NAC)-type transcription factors, namely, ZmNAC128 and ZmNAC130, in endosperm filling. The gene-edited double mutant zmnac128 zmnac130 exhibits a poorly filled kernel phenotype such that the kernels have an inner cavity. RNA sequencing and protein abundance analysis revealed that the expression of many genes involved in the biosynthesis of zein and starch is reduced in the filling endosperm of zmnac128 zmnac130. Further, DNA affinity purification and sequencing combined with chromatin-immunoprecipitation quantitative PCR and promoter transactivation assays demonstrated that ZmNAC128 and ZmNAC130 are direct regulators of 3 (16-, 27-, and 50-kD) γ-zein genes and 6 important starch metabolism genes (Brittle2 [Bt2], pullulanase-type starch debranching enzyme [Zpu1], granule-bound starch synthase 1 [GBSS1], starch synthase 1 [SS1], starch synthase IIa [SSIIa], and sucrose synthase 1 [Sus1]). ZmNAC128 and ZmNAC130 recognize an additional cis-element in the Opaque2 (O2) promoter to regulate its expression. The triple mutant zmnac128 zmnac130 o2 exhibits extremely poor endosperm filling, which results in more than 70% of kernel weight loss. ZmNAC128 and ZmNAC130 regulate the expression of the transporter genes sugars that will eventually be exported transporter 4c (ZmSWEET4c), sucrose and glucose carrier 1 (ZmSUGCAR1), and yellow stripe-like2 (ZmYSL2) and in turn facilitate nutrient uptake, while O2 plays a supporting role. In conclusion, ZmNAC128 and ZmNAC130 cooperate with O2 to facilitate endosperm filling, which involves nutrient uptake in the basal endosperm transfer layer (BETL) and the synthesis of zeins and starch in the starchy endosperm (SE).

Autophagy during maize endosperm development dampens oxidative stress and promotes mitochondrial clearance

The selective turnover of macromolecules by autophagy provides a critical homeostatic mechanism for recycling cellular constituents and for removing superfluous and damaged organelles, membranes, and proteins. To better understand how autophagy impacts seed maturation and nutrient storage, we studied maize (Zea mays) endosperm in its early and middle developmental stages via an integrated multiomic approach using mutants impacting the core macroautophagy factor AUTOPHAGY (ATG)-12 required for autophagosome assembly. Surprisingly, the mutant endosperm in these developmental windows accumulated normal amounts of starch and Zein storage proteins. However, the tissue acquired a substantially altered metabolome, especially for compounds related to oxidative stress and sulfur metabolism, including increases in cystine, dehydroascorbate, cys-glutathione disulfide, glucarate, and galactarate, and decreases in peroxide and the antioxidant glutathione. While changes in the associated transcriptome were mild, the proteome was strongly altered in the atg12 endosperm, especially for increased levels of mitochondrial proteins without a concomitant increase in mRNA abundances. Although fewer mitochondria were seen cytologically, a heightened number appeared dysfunctional based on the accumulation of dilated cristae, consistent with attenuated mitophagy. Collectively, our results confirm that macroautophagy plays a minor role in the accumulation of starch and storage proteins during maize endosperm development but likely helps protect against oxidative stress and clears unneeded/dysfunctional mitochondria during tissue maturation.

Multiomics comparative analysis of the maize large grain mutant tc19 identified pathways related to kernel development

BACKGROUND

The mechanism of grain development in elite maize breeding lines has not been fully elucidated. Grain length, grain width and grain weight are key components of maize grain yield. Previously, using the Chinese elite maize breeding line Chang7-2 and its large grain mutant tc19, we characterized the grain size developmental difference between Chang7-2 and tc19 and performed transcriptomic analysis.

RESULTS

In this paper, using Chang7-2 and tc19, we performed comparative transcriptomic, proteomic and metabolomic analyses at different grain development stages. Through proteomics analyses, we found 2884, 505 and 126 differentially expressed proteins (DEPs) at 14, 21 and 28 days after pollination, respectively. Through metabolomics analysis, we identified 51, 32 and 36 differentially accumulated metabolites (DAMs) at 14, 21 and 28 days after pollination, respectively. Through multiomics comparative analysis, we showed that the phenylpropanoid pathways are influenced at transcriptomic, proteomic and metabolomic levels in all the three grain developmental stages.

CONCLUSION

We identified several genes in phenylpropanoid biosynthesis, which may be related to the large grain phenotype of tc19. In summary, our results provided new insights into maize grain development.

Large-scale long terminal repeat insertions produced a significant set of novel transcripts in cotton

Genomic analysis has revealed that the 1,637-Mb Gossypium arboreum genome contains approximately 81% transposable elements (TEs), while only 57% of the 735-Mb G. raimondii genome is occupied by TEs. In this study, we investigated whether there were unknown transcripts associated with TE or TE fragments and, if so, how these new transcripts were evolved and regulated. As sequence depths increased from 4 to 100 G, a total of 10,284 novel intergenic transcripts (intergenic genes) were discovered. On average, approximately 84% of these intergenic transcripts possibly overlapped with the long terminal repeat (LTR) insertions in the otherwise untranscribed intergenic regions and were expressed at relatively low levels. Most of these intergenic transcripts possessed no transcription activation markers, while the majority of the regular genic genes possessed at least one such marker. Genes without transcription activation markers formed their+1 and -1 nucleosomes more closely (only (117±1.4)bp apart), while twice as big spaces (approximately (403.5±46.0) bp apart) were detected for genes with the activation markers. The analysis of 183 previously assembled genomes across three different kingdoms demonstrated systematically that intergenic transcript numbers in a given genome correlated positively with its LTR content. Evolutionary analysis revealed that genic genes originated during one of the whole-genome duplication events around 137.7 million years ago (MYA) for all eudicot genomes or 13.7 MYA for the Gossypium family, respectively, while the intergenic transcripts evolved around 1.6 MYA, resultant of the last LTR insertion. The characterization of these low-transcribed intergenic transcripts can facilitate our understanding of the potential biological roles played by LTRs during speciation and diversifications.

Soybean gene co-expression network analysis identifies two co-regulated gene modules associated with nodule formation and development

Gene co-expression network analysis is an efficient systems biology approach for the discovery of novel gene functions and trait-associated gene modules. To identify clusters of functionally related genes involved in soybean nodule formation and development, we performed a weighted gene co-expression network analysis. Two nodule-specific modules (NSM-1 and NSM-2, containing 304 and 203 genes, respectively) were identified. The NSM-1 gene promoters were significantly enriched in cis-binding elements for ERF, MYB, and C2H2-type zinc transcription factors, whereas NSM-2 gene promoters were enriched in cis-binding elements for TCP, bZIP, and bHLH transcription factors, suggesting a role of these regulatory factors in the transcriptional activation of nodule co-expressed genes. The co-expressed gene modules included genes with potential novel roles in nodulation, including those involved in xylem development, transmembrane transport, the ethylene signalling pathway, cytoskeleton organization, cytokinesis and regulation of the cell cycle, regulation of meristem initiation and growth, transcriptional regulation, DNA methylation, and histone modifications. Functional analysis of two co-expressed genes using TILLING mutants provided novel insight into the involvement of unsaturated fatty acid biosynthesis and folate metabolism in nodule formation and development. The identified gene co-expression modules provide valuable resources for further functional genomics studies to dissect the genetic basis of nodule formation and development in soybean.

A dolabralexin-deficient mutant provides insight into specialized diterpenoid metabolism in maize

Two major groups of specialized metabolites in maize (Zea mays), termed kauralexins and dolabralexins, serve as known or predicted diterpenoid defenses against pathogens, herbivores, and other environmental stressors. To consider the physiological roles of the recently discovered dolabralexin pathway, we examined dolabralexin structural diversity, tissue-specificity, and stress-elicited production in a defined biosynthetic pathway mutant. Metabolomics analyses support a larger number of dolabralexin pathway products than previously known. We identified dolabradienol as a previously undetected pathway metabolite and characterized its enzymatic production. Transcript and metabolite profiling showed that dolabralexin biosynthesis and accumulation predominantly occur in primary roots and show quantitative variation across genetically diverse inbred lines. Generation and analysis of CRISPR-Cas9-derived loss-of-function Kaurene Synthase-Like 4 (Zmksl4) mutants demonstrated dolabralexin production deficiency, thus supporting ZmKSL4 as the diterpene synthase responsible for the conversion of geranylgeranyl pyrophosphate precursors into dolabradiene and downstream pathway products. Zmksl4 mutants further display altered root-to-shoot ratios and root architecture in response to water deficit. Collectively, these results demonstrate dolabralexin biosynthesis via ZmKSL4 as a committed pathway node biochemically separating kauralexin and dolabralexin metabolism, and suggest an interactive role of maize dolabralexins in plant vigor during abiotic stress.

Revisiting maize Brittle endosperm-2 reveals new insights in BETL development and starchy endosperm filling

Rerouting the starch biosynthesis pathway in maize can generate specialty types, like sweet corn and waxy corn, with a drastically increasing global demand. Hence, a fine-tuning of starch metabolism is relevant to create diverse maize cultivars for end-use applications. Here, we characterized a new maize brittle endosperm mutant, referred to as bt1774, which exhibited decreased starch content but a dramatic increase of soluble sugars at maturity. Both endosperm and embryo development was impaired in bt1774 relative to the wild-type (WT), with a prominently arrested basal endosperm transfer layer (BETL). Map-based cloning revealed that BRITTLE ENDOSPERM2 (Bt2), which encodes a small subunit of ADP-glucose pyrophosphorylase (AGPase), is the causal gene for bt1774. A MuA2 element was found to be inserted into intron 2 of Bt2, leading to a severe decrease of its expression, in bt1774. This is in line with the irregular and loosely packed starch granules in the mutant. Transcriptome of endosperm at grain filling stage identified 1, 013 differentially expressed genes in bt1774, which were notably enriched in the BETL compartment, including ZmMRP1, Miniature1, MEG1, and BETLs. Gene expression of the canonical starch biosynthesis pathway was marginally disturbed in Bt1774. Combined with the residual 60% of starch in this nearly null mutant of Bt2, this data strongly suggests that an AGPase-independent pathway compensates for starch synthesis in the endosperm. Consistent with the BETL defects, zein accumulation was impaired in bt1774. Co-expression network analysis revealed that Bt2 probably has a role in intracellular signal transduction, besides starch synthesis. Altogether, we propose that Bt2 is likely involved in carbohydrate flux and balance, thus regulating both the BETL development and the starchy endosperm filling.

Genetic dissection of embryo size and weight related traits for enhancement of kernel oil in maize

Embryo is a key determinant of kernel-oil in maize. Higher calorific value of maize kernel is attributed to increment in kernel-oil and it stores in specialised structure called embryo. Understanding the genetic behaviour of embryo size and weight related-traits is inevitable task for genetic improvement of kernel-oil. Here, the six-basic generations (P₁, P₂, F₁, F₂, BC₁P₁ and BC₁P₂) of three crosses (CRPBIO-962 × EC932601, CRPBIO-973 × CRPBIO-966 and CRPBIO-966 × CRPBIO-979) between contrasting embryo-sized maize inbreds were field evaluated at three locations to decipher the genetics of twenty embryo, kernel and embryo-to-kernel related-traits through generation-mean-analysis (GMA). Combined ANOVA revealed the significance of all the traits among generations; however, location and generation × location were found to be non-significant (P > 0.05) for most of the traits. Significance (P < 0.05) of scaling and joint-scaling tests revealed the presence of non-allelic interactions. Elucidation of six-parameters disclosed the predominance of dominance main-effect (h) and dominance × dominance interaction-effect (l) for most of traits. The signs of (h) and (l) indicated the prevalence of duplicate-epistasis type across crosses and locations. Thus, the population improvement approaches along with heterosis breeding method could be effective for improvement of these traits. Quantitative inheritance pattern was observed for all the traits with high broad-sense heritability and better-stability across locations. The study also predicted one to three major-gene blocks/QTLs for embryo-traits and up to 11 major-gene blocks/QTLs for embryo-to-kernel traits. These findings could provide deep insights to strategize extensive breeding methods to improve embryo traits for enhancing kernel-oil in sustainable manner.

Spatiotemporal transcriptomics and metabolic profiling provide insights into gene regulatory networks during lentil seed development

Lentil (Lens culinaris Medik.) is a nutritious legume with seeds rich in protein, minerals and an array of diverse specialized metabolites. The formation of a seed requires regulation and tight coordination of developmental programs to form the embryo, endosperm and seed coat compartments, which determines the structure and composition of mature seed and thus its end-use quality. Understanding the molecular and cellular events and metabolic processes of seed development is essential for improving lentil yield and seed nutritional value. However, such information remains largely unknown, especially at the seed compartment level. In this study, we generated high-resolution spatiotemporal gene expression profiles in lentil embryo, seed coat and whole seeds from fertilization through maturation. Apart from anatomic differences between the embryo and seed coat, comparative transcriptomics and weighted gene co-expression network analysis revealed embryo- and seed coat-specific genes and gene modules predominant in specific tissues and stages, which highlights distinct genetic programming. Furthermore, we investigated the dynamic profiles of flavonoid, isoflavone, phytic acid and saponin in seed compartments across seed development. Coupled with transcriptome data, we identified sets of candidate genes involved in the biosynthesis of these metabolites. The global view of the transcriptional and metabolic changes of lentil seed tissues throughout development provides a valuable resource for dissecting the genetic control of secondary metabolism and development of molecular tools for improving seed nutritional quality.

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.

Genome-wide association analysis of kernel nutritional quality in two natural maize populations

As one of the three staple crops, nutritional traits in maize are important for human and animal nutrition. Grain quality-related traits are closely related to grain commercial value. Understanding the genetic basis of quality-related traits in maize would be helpful for breeding high-quality maize varieties. In this study, two association panels (AM122 and AM180) were subjected to genome-wide association analysis of grain quality-related traits, including protein content, oil content, starch content, and fiber content. In total, 98 SNPs (P < 1 × 10-4) were identified to be significantly associated with these four grain quality-related traits. By integrating two sets of public transcriptome data, 31 genes located in 200 kb regions flanking the associated SNP showed high expression during kernel development and were differentially expressed in two maize inbred lines, KA225 and KB035, with significantly different quality. These genes might regulate maize grain quality by participating in plant hormone processes, autophagy processes, and others. All these results could provide important reference information for breeding high‑quality maize varieties.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-023-01360-w.

Maize tubulin folding cofactor B is required for cell division and cell growth through modulating microtubule homeostasis

Tubulin folding cofactors (TFCs) are required for tubulin folding, α/β tubulin heterodimer formation, and microtubule (MT) dynamics in yeast and mammals. However, the functions of their plant counterparts remain to be characterized. We identified a natural maize crumpled kernel mutant, crk2, which exhibits reductions in endosperm cell number and size, as well as embryo/seedling lethality. Map-based cloning and functional complementation confirmed that ZmTFCB is causal for the mutation. ZmTFCB is targeted mainly to the cytosol. It facilitates α-tubulin folding and heterodimer formation through sequential interactions with the cytosolic chaperonin-containing TCP-1 ε subunit ZmCCT5 and ZmTFCE, thus affecting the organization of both the spindle and phragmoplast MT array and the cortical MT polymerization and array formation, which consequently mediated cell division and cell growth. We detected a physical association between ZmTFCB and the maize MT plus-end binding protein END-BINDING1 (ZmEB1), indicating that ZmTFCB1 may modulate MT dynamics by sequestering ZmEB1. Our data demonstrate that ZmTFCB is required for cell division and cell growth through modulating MT homeostasis, an evolutionarily conserved machinery with some species-specific divergence.

High temporal-resolution transcriptome landscapes of maize embryo sac and ovule during early seed development

Here we provided a high temporal-resolution transcriptome atlas of maize embryo sac and ovule to reveal the gene activity dynamic during early seed development. The early maize (Zea mays) seed development is initiated from double fertilization in the embryo sac and needs to undergo a highly dynamic and complex development process to form the differentiated embryo and endosperm. Despite the importance of maize seed for food, feed, and biofuel, many regulators responsible for controlling its early development are not known yet. Here, we reported a high temporal-resolution transcriptome atlas of embryo sac and ovule based on 44 time point samples collected within the first four days of seed development. A total of 25,187 genes including 1598 transcription factors (TFs) involved in early seed development were detected. Global comparisons of the expressions of these genes revealed five distinct development stages of early seed, which are mainly related to double fertilization, asymmetric cell division of the zygote, as well as coenocyte formation, cellularization and differentiation in endosperm. We identified 3327 seed-specific genes, which more than one thousand seed-specific genes with main expressions during early seed development were newly identified here, including 859 and 186 genes predominantly expressed in the embryo sac and ovule, respectively. Combined with the published transcriptome data of seed, we uncovered the dominant auxin biosynthesis, transport and signaling related genes at different development stages and subregions of seed. These results are helpful for understanding the genetic control of early seed development.

Transcriptomic profiling of wheat stem during meiosis in response to freezing stress

Low temperature injury in spring has seriously destabilized the production and grain quality of common wheat. However, the molecular mechanisms underlying spring frost tolerance remain elusive. In this study, we investigated the response of a frost-tolerant wheat variety Zhongmai8444 to freezing stress at the meiotic stage. Transcriptome profiles over a time course were subsequently generated by high-throughput sequencing. Our results revealed that the prolonged freezing temperature led to the significant reductions in plant height and seed setting rate. Cell wall thickening in the vascular tissue was also observed in the stems. RNA-seq analyses demonstrated the identification of 1010 up-regulated and 230 down-regulated genes shared by all time points of freezing treatment. Enrichment analysis revealed that gene activity related to hormone signal transduction and cell wall biosynthesis was significantly modulated under freezing. In addition, among the identified differentially expressed genes, 111 transcription factors belonging to multiple gene families exhibited dynamic expression pattern. This study provided valuable gene resources beneficial for the breeding of wheat varieties with improved spring frost tolerance.

Regulation of seed storage protein synthesis in monocot and dicot plants: A comparative review

Seeds are a major source of nutrients for humans and animal livestock worldwide. With improved living standards, high nutritional quality has become one of the main targets for breeding. Storage protein content in seeds, which is highly variable depending on plant species, serves as a pivotal criterion of seed nutritional quality. In the last few decades, our understanding of the molecular genetics and regulatory mechanisms of storage protein synthesis has greatly advanced. Here, we systematically and comprehensively summarize breakthroughs on the conservation and divergence of storage protein synthesis in dicot and monocot plants. With regard to storage protein accumulation, we discuss evolutionary origins, developmental processes, characteristics of main storage protein fractions, regulatory networks, and genetic modifications. In addition, we discuss potential breeding strategies to improve storage protein accumulation and provide perspectives on some key unanswered problems that need to be addressed.

Nitrogen-dependent binding of the transcription factor PBF1 contributes to the balance of protein and carbohydrate storage in maize endosperm

Fluctuations in nitrogen (N) availability influence protein and starch levels in maize (Zea mays) seeds, yet the underlying mechanism is not well understood. Here, we report that N limitation impacted the expression of many key genes in N and carbon (C) metabolism in the developing endosperm of maize. Notably, the promoter regions of those genes were enriched for P-box sequences, the binding motif of the transcription factor prolamin-box binding factor 1 (PBF1). Loss of PBF1 altered accumulation of starch and proteins in endosperm. Under different N conditions, PBF1 protein levels remained stable but PBF1 bound different sets of target genes, especially genes related to the biosynthesis and accumulation of N and C storage products. Upon N-starvation, the absence of PBF1 from the promoters of some zein genes coincided with their reduced expression, suggesting that PBF1 promotes zein accumulation in the endosperm. In addition, PBF1 repressed the expression of sugary1 (Su1) and starch branching enzyme 2b (Sbe2b) under normal N supply, suggesting that, under N-deficiency, PBF1 redirects the flow of C skeletons for zein toward the formation of C compounds. Overall, our study demonstrates that PBF1 modulates C and N metabolism during endosperm development in an N-dependent manner.

Conservation Study of Imprinted Genes in Maize Triparental Heterozygotic Kernels

Genomic imprinting is a classic epigenetic phenomenon related to the uniparental expression of genes. Imprinting variability exists in seeds and can contribute to observed parent-of-origin effects on seed development. Here, we conducted allelic expression of the embryo and endosperm from four crosses at 11 days after pollination (DAP). First, the F1 progeny of B73(♀) × Mo17(♂) and the inducer line CAU5 were used as parents to obtain reciprocal crosses of BM-C/C-BM. Additionally, the F1 progeny of Mo17(♀) × B73(♂) and CAU5 were used as parents to obtain reciprocal crosses of MB-C/C-MB. In total, 192 and 181 imprinted genes were identified in the BM-C/C-BM and MB-C/C-MB crosses, respectively. Then, by comparing the allelic expression of these imprinted genes in the reciprocal crosses of B73 and CAU5 (BC/CB), fifty-one Mo17-added non-conserved genes were identified as exhibiting imprinting variability. Fifty-one B73-added non-conserved genes were also identified by comparing the allelic expression of imprinted genes identified in BM-C/C-BM, MB-C/C-MB and MC/CM crosses. Specific Gene Ontology (GO) terms were not enriched in B73-added/Mo17-added non-conserved genes. Interestingly, the imprinting status of these genes was less conserved across other species. The cis-element distribution, tissue expression and subcellular location were similar between the B73-added/Mo17-added conserved and B73-added/Mo17-added non-conserved imprinted genes. Finally, genotypic and phenotypic analysis of one non-conserved gene showed that the mutation and overexpression of this gene may affect embryo and kernel size, which indicates that these non-conserved genes may also play an important role in kernel development. The findings of this study will be helpful for elucidating the imprinting mechanism of genes involved in maize kernel development.

THP9 enhances seed protein content and nitrogen-use efficiency in maize

Teosinte, the wild ancestor of maize (Zea mays subsp. mays), has three times the seed protein content of most modern inbreds and hybrids, but the mechanisms that are responsible for this trait are unknown^1,2. Here we use trio binning to create a contiguous haplotype DNA sequence of a teosinte (Zea mays subsp. parviglumis) and, through map-based cloning, identify a major high-protein quantitative trait locus, TEOSINTE HIGH PROTEIN 9 (THP9), on chromosome 9. THP9 encodes an asparagine synthetase 4 enzyme that is highly expressed in teosinte, but not in the B73 inbred, in which a deletion in the tenth intron of THP9-B73 causes incorrect splicing of THP9-B73 transcripts. Transgenic expression of THP9-teosinte in B73 significantly increased the seed protein content. Introgression of THP9-teosinte into modern maize inbreds and hybrids greatly enhanced the accumulation of free amino acids, especially asparagine, throughout the plant, and increased seed protein content without affecting yield. THP9-teosinte seems to increase nitrogen-use efficiency, which is important for promoting a high yield under low-nitrogen conditions.

Maize sterility gene DRP1 encodes a desiccation-related protein that is critical for Ubisch bodies and pollen exine development

Desiccation tolerance is a remarkable feature of pollen, seeds, and resurrection-type plants. Exposure to desiccation stress can cause sporophytic defects, resulting in male sterility. Here, we report the novel maize sterility gene DRP1 (Desiccation-Related Protein 1), which was identified by bulked-segregant analysis sequencing and encodes a desiccation-related protein. Loss of function of DRP1 results in abnormal Ubisch bodies, defective tectum of the pollen exine, and complete male sterility. Our results suggest that DRP1 may facilitate anther dehydration to maintain appropriate water status. DRP1 is a secretory protein that is specifically expressed in the tapetum and microspore from the tetrad to the uninucleate microspore stage. Differentially expressed genes in drp1 are enriched in Gene Ontology terms for pollen exine formation, polysaccharide catabolic process, extracellular region, and response to heat. In addition, DRP1 is a target of selection that appears to have played an important role in the spread of maize from tropical/subtropical to temperate regions. Taken together, our results suggest that DRP1 encodes a desiccation-related protein whose loss of function causes male sterility. Our findings provide a potential genetic resource that may be used to design crops for heterosis utilization.

The sugar transporter ZmSUGCAR1 of the nitrate transporter 1/peptide transporter family is critical for maize grain filling

Maternal-to-filial nutrition transfer is central to grain development and yield. nitrate transporter 1/peptide transporter (NRT1-PTR)-type transporters typically transport nitrate, peptides, and ions. Here, we report the identification of a maize (Zea mays) NRT1-PTR-type transporter that transports sucrose and glucose. The activity of this sugar transporter, named Sucrose and Glucose Carrier 1 (SUGCAR1), was systematically verified by tracer-labeled sugar uptake and serial electrophysiological studies including two-electrode voltage-clamp, non-invasive microelectrode ion flux estimation assays in Xenopus laevis oocytes and patch clamping in HEK293T cells. ZmSUGCAR1 is specifically expressed in the basal endosperm transfer layer and loss-of-function mutation of ZmSUGCAR1 caused significantly decreased sucrose and glucose contents and subsequent shrinkage of maize kernels. Notably, the ZmSUGCAR1 orthologs SbSUGCAR1 (from Sorghum bicolor) and TaSUGCAR1 (from Triticum aestivum) displayed similar sugar transport activities in oocytes, supporting the functional conservation of SUGCAR1 in closely related cereal species. Thus, the discovery of ZmSUGCAR1 uncovers a type of sugar transporter essential for grain development and opens potential avenues for genetic improvement of seed-filling and yield in maize and other grain crops.

Paternal imprinting of dosage-effect defective1 contributes to seed weight xenia in maize

Historically, xenia effects were hypothesized to be unique genetic contributions of pollen to seed phenotype, but most examples represent standard complementation of Mendelian traits. We identified the imprinted dosage-effect defective1 (ded1) locus in maize (Zea mays) as a paternal regulator of seed size and development. Hypomorphic alleles show a 5–10% seed weight reduction when ded1 is transmitted through the male, while homozygous mutants are defective with a 70–90% seed weight reduction. Ded1 encodes an R2R3-MYB transcription factor expressed specifically during early endosperm development with paternal allele bias. DED1 directly activates early endosperm genes and endosperm adjacent to scutellum cell layer genes, while directly repressing late grain-fill genes. These results demonstrate xenia as originally defined: Imprinting of Ded1 causes the paternal allele to set the pace of endosperm development thereby influencing grain set and size.

Xenia effects describe the genetic contribution of pollen to seed phenotypes. Here the authors show that paternal imprinting of Ded1 contributes to the xenia effect in maize by setting the pace of endosperm development.

Transcriptome and metabolome changes in Chinese cedar during cold acclimation reveal the roles of flavonoids in needle discoloration and cold resistance

Cryptomeria fortunei growth and development are usually affected by low temperatures. Despite the evergreen nature of this species, most needles turn yellowish-brown in cold winters. The underlying discoloration mechanisms that cause this phenomenon in response to cold acclimation remain poorly understood. Here, we measured the pigment content and ultrastructure of normal wild-type (Wt) and evergreen mutant (GM) C. fortunei needles and performed integrated transcriptomic and metabolomic analyses to explore potential discoloration mechanisms. The results showed that the needle chlorophyll content of these two genotypes decreased in winter. Wt needles showed greater decrease in the chlorophyll content and local destruction of chloroplast ultrastructure and contained larger amounts of flavonoids than GM needles, as shown by metabolomics analysis. We subsequently identified key differentially expressed genes in the flavonoid biosynthesis pathway and observed significantly upregulated flavonol synthase expression in Wt needles compared with GM needles that significantly increased the anthoxanthin (flavones and flavonols) content, which is likely a key factor underlying the difference in needle color between these two genotypes. Therefore, flavonoid metabolism may play important roles in the cold resistance and needle discoloration of C. fortunei, and our results provide an excellent foundation for the molecular mechanism of C. fortunei in response to cold stress.

easyMF: A Web Platform for Matrix Factorization-Based Gene Discovery from Large-scale Transcriptome Data

With the development of high-throughput experimental technologies, large-scale RNA sequencing (RNA-Seq) data have been and continue to be produced, but have led to challenges in extracting relevant biological knowledge hidden in the produced high-dimensional gene expression matrices. Here, we develop easyMF ( https://github.com/cma2015/easyMF ), a web platform that can facilitate functional gene discovery from large-scale transcriptome data using matrix factorization (MF) algorithms. Compared with existing MF-based software packages, easyMF exhibits several promising features, such as greater functionality, flexibility and ease of use. The easyMF platform is equipped using the Big-Data-supported Galaxy system with user-friendly graphic user interfaces, allowing users with little programming experience to streamline transcriptome analysis from raw reads to gene expression, carry out multiple-scenario MF analysis, and perform multiple-way MF-based gene discovery. easyMF is also powered with the advanced packing technology to enhance ease of use under different operating systems and computational environments. We illustrated the application of easyMF for seed gene discovery from temporal, spatial, and integrated RNA-Seq datasets of maize (Zea mays L.), resulting in the identification of 3,167 seed stage-specific, 1,849 seed compartment-specific, and 774 seed-specific genes, respectively. The present results also indicated that easyMF can prioritize seed-related genes with superior prediction performance over the state-of-art network-based gene prioritization system MaizeNet. As a modular, containerized and open-source platform, easyMF can be further customized to satisfy users' specific demands of functional gene discovery and deployed as a web service for broad applications.

Transcriptomes across fertilization and seed development in the water lily Nymphaea thermarum (Nymphaeales): evidence for epigenetic patterning during reproduction

The first record of gene expression during seed development within the Nymphaeales provides evidence for a variety of biological processes, including dynamic epigenetic patterning during sexual reproduction in the water lily Nymphaea thermarum. Studies of gene expression during seed development have been performed for a growing collection of species from a phylogenetically broad sampling of flowering plants (angiosperms). However, angiosperm lineages whose origins predate the divergence of monocots and eudicots have been largely overlooked. In order to provide a new resource for understanding the early evolution of seed development in flowering plants, we sequenced transcriptomes of whole ovules and seeds from three key stages of reproductive development in the waterlily Nymphaea thermarum, an experimentally tractable member of the Nymphaeales. We first explore patterns of gene expression, beginning with mature ovules and continuing through fertilization into early- and mid-stages of seed development. We find patterns of gene expression that corroborate histological/morphological observations of seed development in this species, such as expression of genes involved in starch synthesis and transcription factors that have been associated with embryo and endosperm development in other species. We also find evidence for processes that were previously not known to be occurring during seed development in this species, such as epigenetic modification. We then examine the expression of genes associated with patterning DNA and histone methylation-processes that are essential for seed development in distantly related and structurally diverse monocots and eudicots. Around 89% of transcripts putatively homologous to DNA and histone methylation modifiers are expressed during seed development in N. thermarum, including homologs of genes known to pattern imprinting-related epigenetic modifications. Our results suggest that dynamic epigenetic patterning is a deeply conserved aspect of angiosperm seed development.

Dynamic transcriptome analysis suggests the key genes regulating seed development and filling in Tartary buckwheat (Fagopyrum tataricum Garetn.)

Tartary buckwheat is highly attractive for the richness of nutrients and quality, yet post-embryonic seed abortion greatly halts the yield. Seed development is crucial for determining grain yield, whereas the molecular basis and regulatory network of Tartary buckwheat seed development and filling is not well understood at present. Here, we assessed the transcriptional dynamics of filling stage Tartary buckwheat seeds at three developmental stages by RNA sequencing. Among the 4249 differentially expressed genes (DEGs), genes related to seed development were identified. Specifically, 88 phytohormone biosynthesis signaling genes, 309 TFs, and 16 expansin genes participating in cell enlargement, 37 structural genes involved in starch biosynthesis represented significant variation and were candidate key seed development genes. Cis-element enrichment analysis indicated that the promoters of differentially expressed expansin genes and starch biosynthesis genes are rich of hormone-responsive (ABA-, AUX-, ET-, and JA-), and seed growth-related (MYB, MYC and WRKY) binding sites. The expansin DEGs showed strong correlations with DEGs in phytohormone pathways and transcription factors (TFs). In total, phytohormone ABA, AUX, ET, BR and CTK, and related TFs could substantially regulate seed development in Tartary buckwheat through targeting downstream expansin genes and structural starch biosynthetic genes. This transcriptome data could provide a theoretical basis for improving yield of Tartary buckwheat.

Epigenetic Mutation in a Tubulin-Folding Cofactor B (ZmTFCB) Gene Arrests Kernel Development in Maize

Epialleles, the heritable epigenetic variants that are not caused by changes in DNA sequences, can broaden genetic and phenotypic diversity and benefit to crop breeding, but very few epialleles related to agricultural traits have been identified in maize. Here, we cloned a small kernel mutant, smk-wl10, from maize, which encoded a tubulin-folding cofactor B (ZmTFCB) protein. Expression of the ZmTFCB gene decreased in the smk-wl10 mutant, which arrested embryo, endosperm and basal endosperm transfer layer developments. Overexpression of ZmTFCB could complement the defective phenotype of smk-wl10. No nucleotide sequence variation in ZmTFCB could be found between smk-wl10 and wild type (WT). Instead, we detected hypermethylation of nucleotide CHG (where H is A, C or T nucleotide) sequence contexts and increased level of histone H3K9me2 methylation in the upstream sequence of ZmTFCB in smk-wl10 compared with WT, which might respond to the attenuating transcription of ZmTFCB. In addition, yeast two-hybrid and bimolecular fluorescence complementation assays identified a strong interaction between ZmTFCB and its homolog ZmTFCE. Thus, our work identifies a novel epiallele of the maize ZmTFCB gene, which might represent a common phenomenon in the epigenetic regulation of important traits such as kernel development in maize.

Identification of two new QTLs of maize (Zea mays L.) underlying kernel row number using the HNAU-NAM1 population

BACKGROUND

Maize kernel row number (KRN) is one of the most important yield traits and has changed greatly during maize domestication and selection. Elucidating the genetic basis of KRN will be helpful to improve grain yield in maize.

RESULTS

Here, we measured KRN in four environments using a nested association mapping (NAM) population named HNAU-NAM1 with 1,617 recombinant inbred lines (RILs) that were derived from 12 maize inbred lines with a common parent, GEMS41. Then, five consensus quantitative trait loci (QTLs) distributing on four chromosomes were identified in at least three environments along with the best linear unbiased prediction (BLUP) values by the joint linkage mapping (JLM) method. These QTLs were further validated by the separate linkage mapping (SLM) and genome-wide association study (GWAS) methods. Three KRN genes cloned through the QTL assay were found in three of the five consensus QTLs, including qKRN1.1, qKRN2.1 and qKRN4.1. Two new QTLs of KRN, qKRN4.2 and qKRN9.1, were also identified. On the basis of public RNA-seq and genome annotation data, five genes highly expressed in ear tissue were considered candidate genes contributing to KRN.

CONCLUSIONS

This study carried out a comprehensive analysis of the genetic architecture of KRN by using a new NAM population under multiple environments. The present results provide solid information for understanding the genetic components underlying KRN and candidate genes in qKRN4.2 and qKRN9.1. Single-nucleotide polymorphisms (SNPs) closely linked to qKRN4.2 and qKRN9.1 could be used to improve inbred yield during molecular breeding in maize.

Distinct Hormone Signalling-Modulation Activities Characterize Two Maize Endosperm-Specific Type-A Response Regulators

ZmTCRR1 and 2 are type-A response regulators expressed in the maize endosperm transfer cells (TC). While type-B response regulators transcriptionally control canonical type-A response regulators, as part of the cytokinin signal transduction mechanism, the ZmTCRRs are regulated by ZmMRP1, a master regulator of TC identity. In addition, the corresponding proteins are not detected in the TC, accumulating in the inner endosperm cells instead. These features suggest these molecules are not involved in classical, cell-autonomous, cytokinin signalling pathways. Using transgenic Arabidopsis plants ectopically expressing these genes, we have shown that ZmTCRR1 and 2 can modulate auxin and cytokinin signalling, respectively. In Arabidopsis, the ectopic expression of ZmTCRR2 blocked, almost completely, cytokinin perception. Given the conservation of these signalling pathways at the molecular level, our results suggest that the ZmTCRRs modulate cytokinin and auxin perception in the inner endosperm cells.

Genome-Wide Characterization and Phylogenetic Analysis of GSK Genes in Maize and Elucidation of Their General Role in Interaction with BZR1

Glycogen synthase kinase-3 (GSK-3) is a nonreceptor serine/threonine protein kinase that is involved in diverse processes, including cell development, photomorphogenesis, biotic and abiotic stress responses, and hormone signaling. In contrast with the deeply researched GSK family in Arabidopsis and rice, maize GSKs’ common bioinformatic features and protein functions are poorly understood. In this study, we identified 11 GSK genes in the maize (Zea mays L.) genome via homologous alignment, which we named Zeama;GSKs (ZmGSKs). The results of ZmGSK protein sequences, conserved motifs, and gene structures showed high similarities with each other. The phylogenetic analyses showed that a total of 11 genes from maize were divided into four clades. Furthermore, semi-quantitative RT-PCR analysis of the GSKs genes showed that ZmGSK1, ZmGSK2, ZmGSK4, ZmGSK5, ZmGSK8, ZmGSK9, ZmGSK10, and ZmGSK11 were expressed in all tissues; ZmGSK3, ZmGSK6, and ZmGSK7 were expressed in a specific organization. In addition, GSK expression profiles under hormone treatments demonstrated that the ZmGSK genes were induced under BR conditions, except for ZmGSK2 and ZmGSK5. ZmGSK genes were regulated under ABA conditions, except for ZmGSK1 and ZmGSK8. Finally, using the yeast two-hybrid and BiFC assay, we determined that clads II (ZmGSK1, ZmGSK4, ZmGSK7, ZmGSK8, and ZmGSK11) could interact with ZmBZR1. The results suggest that clade II of ZmGSKs is important for BR signaling and that ZmGSK1 may play a dominant role in BR signaling as the counterpart to BIN2. This study provides a foundation for the further study of GSK3 functions and could be helpful in devising strategies for improving maize.

Meta-QTL analysis explores the key genes, especially hormone related genes, involved in the regulation of grain water content and grain dehydration rate in maize

BACKGROUND

Low grain water content (GWC) at harvest of maize (Zea mays L.) is essential for mechanical harvesting, transportation and storage. Grain drying rate (GDR) is a key determinant of GWC. Many quantitative trait locus (QTLs) related to GDR and GWC have been reported, however, the confidence interval (CI) of these QTLs are too large and few QTLs has been fine-mapped or even been cloned. Meta-QTL (MQTL) analysis is an effective method to integrate QTLs information in independent populations, which helps to understand the genetic structure of quantitative traits.

RESULTS

In this study, MQTL analysis was performed using 282 QTLs from 25 experiments related GDR and GWC. Totally, 11 and 34 MQTLs were found to be associated with GDR and GWC, respectively. The average CI of GDR and GWC MQTLs was 24.44 and 22.13 cM which reduced the 57 and 65% compared to the average QTL interval for initial GDR and GWC QTL, respectively. Finally, 1494 and 5011 candidate genes related to GDR and GWC were identified in MQTL intervals, respectively. Among these genes, there are 48 genes related to hormone metabolism.

CONCLUSIONS

Our studies combined traditional QTL analyses, genome-wide association study and RNA-seq to analysis major locus for regulating GWC in maize.

Exploring key developmental phases and phase-specific genes across the entirety of anther development in maize

Anther development from stamen primordium to pollen dispersal is complex and essential to sexual reproduction. How this highly dynamic and complex developmental process is controlled genetically is not well understood, especially for genes involved in specific key developmental phases. Here we generated RNA sequencing libraries spanning 10 key stages across the entirety of anther development in maize (Zea mays). Global transcriptome analyses revealed distinct phases of cell division and expansion, meiosis, pollen maturation, and mature pollen, for which we detected 50, 245, 42, and 414 phase-specific marker genes, respectively. Phase-specific transcription factor genes were significantly enriched in the phase of meiosis. The phase-specific expression of these marker genes was highly conserved among the maize lines Chang7-2 and W23, indicating they might have important roles in anther development. We explored a desiccation-related protein gene, ZmDRP1, which was exclusively expressed in the tapetum from the tetrad to the uninucleate microspore stage, by generating knockout mutants. Notably, mutants in ZmDRP1 were completely male-sterile, with abnormal Ubisch bodies and defective pollen exine. Our work provides a glimpse into the gene expression dynamics and a valuable resource for exploring the roles of key phase-specific genes that regulate anther development.

Starch metabolism in potato <i>Solanum tuberosum</i> L

Starch is a major storage carbohydrate in plants. It is an important source of calories in the human and animal diet. Also, it is widely used in various industries. Native starch consists of water-insoluble semicrystalline granules formed by natural glucose polymers amylose and amylopectin. The physicochemical properties of starch are determined by the amylose:amylopectin ratio in the granule and degrees of their polymerization and phosphorylation. Potato Solanum tuberosum L. is one of the main starch-producing crops. Growing industrial needs necessitate the breeding of plant varieties with increased starch content and specified starch properties. This task demands detailed information on starch metabolism in the producing plant. It is a complex process, requiring the orchestrated work of many enzymes, transporter and targeting proteins, transcription factors, and other regulators. Two types of starch are recognized with regard to their biological functions. Transitory starch is synthesized in chloroplasts of photosynthetic organs and degraded in the absence of light, providing carbohydrates for cell needs. Storage starch is synthesized and stored in amyloplasts of storage organs: grains and tubers. The main enzymatic reactions of starch biosynthesis and degradation, as well as carbohydrate transport and metabolism, are well known in the case of transitory starch of the model plant Arabidopsis thaliana. Less is known about features of starch metabolism in storage organs, in particular, potato tubers. Several issues remain obscure: the roles of enzyme isoforms and different regulatory factors in tissues at various plant developmental stages and under different environmental conditions; alternative enzymatic processes; targeting and transport proteins. In this review, the key enzymatic reactions of plant carbohydrate metabolism, transitory and storage starch biosynthesis, and starch degradation are discussed, and features specific for potato are outlined. Attention is also paid to the known regulatory factors affecting starch metabolism