RNA sequencing of laser-capture microdissected compartments of the maize kernel identifies regulatory modules associated with endosperm cell differentiation

Integrated analysis of transcriptomics and defense-related phytohormones to discover hub genes conferring maize Gibberella ear rot caused by Fusarium Graminearum

BACKGROUND

Gibberella ear rot (GER) is one of the most devastating diseases in maize growing areas, which directly reduces grain yield and quality. However, the underlying defense response of maize to pathogens infection is largely unknown.

RESULTS

To gain a comprehensive understanding of the defense response in GER resistance, two contrasting inbred lines 'Nov-82' and 'H10' were used to explore transcriptomic profiles and defense-related phytohormonal alterations during Fusarium graminearum infection. Transcriptomic analysis revealed 4,417 and 4,313 differentially expressed genes (DEGs) from the Nov-82 and H10, respectively, and 647 common DEGs between the two lines. More DEGs were obviously enriched in phenylpropanoid biosynthesis, secondary metabolites biosynthesis, metabolic process and defense-related pathways. In addition, the concentration of the defense-related phytohormones, jasmonates (JAs) and salicylates (SAs), was greatly induced after the pathogen infection. The level of JAs in H10 was more higher than in Nov-82, whereas an opposite pattern for the SA between the both lines. Integrated analysis of the DEGs and the phytohormones revealed five vital modules based on co-expression network analysis according to their correlation. A total of 12 hub genes encoding fatty acid desaturase, subtilisin-like protease, ethylene-responsive transcription factor, 1-aminocyclopropane-1-carboxylate oxidase, and sugar transport protein were captured from the key modules, indicating that these genes might play unique roles in response to pathogen infection, CONCLUSIONS: Overall, our results indicate that large number DEGs related to plant disease resistance and different alteration of defensive phytohormones were activated during F. graminearum infection, providing new insight into the defense response against pathogen invasion, in addition to the identified hub genes that can be further investigated for enhancing maize GER resistance.

Developing a simple and rapid method for cell-specific transcriptome analysis through laser microdissection: insights from citrus rind with broader implications

BACKGROUND

With the rapid development of single-cell sequencing technology, histological studies are no longer limited to conventional homogenized tissues. Laser microdissection enables the accurate isolation of specific tissues or cells, and when combined with next-generation sequencing, it can reveal important biological processes at the cellular level. However, traditional laser microdissection techniques have often been complicated and time-consuming, and the quality of the RNA extracted from the collected samples has been inconsistent, limiting follow-up studies. Therefore, an improved, simple, and efficient laser microdissection method is urgently needed.

RESULTS

We omitted the sample fixation and cryoprotectant addition steps. Instead, fresh samples were embedded in Optimal Cutting Temperature medium within 1.5 ml centrifuge tube caps, rapidly frozen with liquid nitrogen, and immediately subjected to cryosectioning. A series of section thicknesses of citrus rind were tested for RNA extraction, which showed that 18 μm thickness yielded the highest quality RNA. By shortening the dehydration time to one minute per ethanol gradient and omitting the tissue clearing step, the resulting efficient dehydration and preserved morphology ensured high-quality RNA extraction. We also propose a set of laser microdissection parameters by adjusting the laser power to optimal values, reducing the aperture size, and lowering the pulse frequency. Both the epidermal and subepidermal cells from the citrus rind were collected, and RNA extraction was completed within nine hours. Using this efficient method, the transcriptome sequencing of the isolated tissues generated high-quality data with average Q30 values and mapping rates exceeding 91%. Moreover, the transcriptome analysis revealed significant differences between the cell layers, further confirming the effectiveness of our isolation approach.

CONCLUSIONS

We developed a simple and rapid laser microdissection method and demonstrated its effectiveness through a study based on citrus rind, from which we generated high-quality transcriptomic data. This fast and efficient method of cell isolation, combined with transcriptome sequencing not only contributes to precise histological studies at the cellular level in citrus but also provides a promising approach for cell-specific transcriptome analysis in a broader range of other plant tissues.

CaAOS as a hub gene based on physiological and transcriptomic analyses of cold-resistant and cold-sensitive pepper cultivars

The yield and quality of pepper are considerably influenced by the cold conditions. Herein, we performed morphological, physiological and transcriptomic analyses by using two pepper seedlings, '2379' (cold-resistant) and '2380' (cold-sensitive). Briefly, 60 samples from each cultivar were analyzed at four distinct time points (0, 6, 24 and 48 h) at 5 °C in darkness. The physiological indices and activities of enzymes exhibited marked differences between the two cultivars. Transcriptomic analysis indicated that, compared to the control group, 11,415 DEGs were identified in '2379' and '2380' at 24 h. In the early stage, the number of DEGs in '2379' was 5.68 times higher than that in '2380', potentially explaining the observed differences in tolerance to colds. Processes such as protein targeting to membranes, jasmonic acid (JA)-mediated signalling, cold response and abscisic acid-activated signalling were involved. Subsequently, we identified a hub gene, CaAOS, that is involved in JA biosynthesis, positively influences cold tolerance and is a target of CaMYC2. Variations in the GC-motif of the CaAOS's promoter may influence the expression levels of CaAOS under cold treatment. The result of this study may lead to the development of more effective strategies for enhancing cold tolerance, potentially benefitting pepper breeding in cold regions.

Transcriptome and metabolome analyses reveal regulatory networks associated with nutrition synthesis in sorghum seeds

Cereal seeds are vital for food, feed, and agricultural sustainability because they store and provide essential nutrients to human and animal food and feed systems. Unraveling molecular processes in seed development is crucial for enhancing cereal grain yield and quality. We analyze spatiotemporal transcriptome and metabolome profiles during sorghum seed development in the inbred line 'BTx623'. Morphological and molecular analyses identify the key stages of seed maturation, specifying starch biosynthesis onset at 5 days post-anthesis (dpa) and protein at 10 dpa. Transcriptome profiling from 1 to 25 dpa reveal dynamic gene expression pathways, shifting from cellular growth and embryo development (1-5 dpa) to cell division, fatty acid biosynthesis (5-25 dpa), and seed storage compounds synthesis in the endosperm (5-25 dpa). Network analysis identifies 361 and 207 hub genes linked to starch and protein synthesis in the endosperm, respectively, which will help breeders enhance sorghum grain quality. The availability of this data in the sorghum reference genome line establishes a baseline for future studies as new pangenomes emerge, which will consider copy number and presence-absence variation in functional food traits.

Opportunities and Challenges in Advancing Plant Research with Single-cell Omics

Plants possess diverse cell types and intricate regulatory mechanisms to adapt to the ever-changing environment of nature. Various strategies have been employed to study cell types and their developmental progressions, including single-cell sequencing methods which provide high-dimensional catalogs to address biological concerns. In recent years, single-cell sequencing technologies in transcriptomics, epigenomics, proteomics, metabolomics, and spatial transcriptomics have been increasingly used in plant science to reveal intricate biological relationships at the single-cell level. However, the application of single-cell technologies to plants is more limited due to the challenges posed by cell structure. This review outlines the advancements in single-cell omics technologies, their implications in plant systems, future research applications, and the challenges of single-cell omics in plant systems.

Recent advances in exploring transcriptional regulatory landscape of crops

Crop breeding entails developing and selecting plant varieties with improved agronomic traits. Modern molecular techniques, such as genome editing, enable more efficient manipulation of plant phenotype by altering the expression of particular regulatory or functional genes. Hence, it is essential to thoroughly comprehend the transcriptional regulatory mechanisms that underpin these traits. In the multi-omics era, a large amount of omics data has been generated for diverse crop species, including genomics, epigenomics, transcriptomics, proteomics, and single-cell omics. The abundant data resources and the emergence of advanced computational tools offer unprecedented opportunities for obtaining a holistic view and profound understanding of the regulatory processes linked to desirable traits. This review focuses on integrated network approaches that utilize multi-omics data to investigate gene expression regulation. Various types of regulatory networks and their inference methods are discussed, focusing on recent advancements in crop plants. The integration of multi-omics data has been proven to be crucial for the construction of high-confidence regulatory networks. With the refinement of these methodologies, they will significantly enhance crop breeding efforts and contribute to global food security.

Insights into dynamic coenocytic endosperm development: Unraveling molecular, cellular, and growth complexity

The endosperm, a product of double fertilization, is one of the keys to the evolution and success of angiosperms in conquering the land. While there are differences in endosperm development among flowering plants, the most common form is coenocytic growth, where the endosperm initially undergoes nuclear division without cytokinesis and eventually becomes cellularized. This complex process requires interplay among networks of transcription factors such as MADS-box, auxin response factors (ARFs), and phytohormones. The role of cytoskeletal elements in shaping the coenocytic endosperm and influencing seed growth also becomes evident. This review offers a recent understanding of the molecular and cellular dynamics in coenocytic endosperm development and their contributions to the final seed size.

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.

A MYB-related transcription factor ZmMYBR29 is involved in grain filling

BACKGROUND

The endosperm serves as the primary source of nutrients for maize (Zea mays L.) kernel embryo development and germination. Positioned at the base of the endosperm, the transfer cells (TCs) of the basal endosperm transfer layer (BETL) generate cell wall ingrowths, which enhance the connectivity between the maternal plant and the developing kernels. These TCs play a crucial role in nutrient transport and defense against pathogens. The molecular mechanism underlying BETL development in maize remains unraveled.

RESULTS

This study demonstrated that the MYB-related transcription factor ZmMYBR29, exhibited specific expression in the basal cellularized endosperm, as evidenced by in situ hybridization analysis. Utilizing the CRISPR/Cas9 system, we successfully generated a loss-of-function homozygous zmmybr29 mutant, which presented with smaller kernel size. Observation of histological sections revealed abnormal development and disrupted morphology of the cell wall ingrowths in the BETL. The average grain filling rate decreased significantly by 26.7% in zmmybr29 mutant in comparison to the wild type, which impacted the dry matter accumulation within the kernels and ultimately led to a decrease in grain weight. Analysis of RNA-seq data revealed downregulated expression of genes associated with starch synthesis and carbohydrate metabolism in the mutant. Furthermore, transcriptomic profiling identified 23 genes that expressed specifically in BETL, and the majority of these genes exhibited altered expression patterns in zmmybr29 mutant.

CONCLUSIONS

In summary, ZmMYBR29 encodes a MYB-related transcription factor that is expressed specifically in BETL, resulting in the downregulation of genes associated with kernel development. Furthermore, ZmMYBR29 influences kernels weight by affecting the grain filling rate, providing a new perspective for the complementation of the molecular regulatory network in maize endosperm development.

Physiological and transcriptomic comparisons shed light on the cold stress response mechanisms of Dendrobium spp

BACKGROUND

Dendrobium spp. comprise a group of tropical orchids with ornamental and medicinal value. Dendrobium spp. are sensitive to low temperature, and the underlying cold response regulatory mechanisms in this group are unclear. To understand how these plants respond to cold stress, we compared the transcriptomic responses of the cold-tolerant cultivar 'Hongxing' (HX) and the cold-sensitive cultivar 'Sonia Hiasakul' (SH) to cold stress.

RESULTS

Chemometric results showed that the physiological response of SH in the later stages of cold stress is similar to that of HX throughout the cold treatment. Orthogonal partial least squares discriminant analysis (OPLS-DA) revealed that soluble protein content and peroxidase activity are key physiological parameters for assessing the cold tolerance of these two Dendrobium spp. cultivars. Additionally, weighted gene co-expression network analysis (WGCNA) results showed that many cold response genes and metabolic pathways significantly associated with the physiological indices were enriched in the 12 detected modules. The Kyoto Encyclopedia of Genes and Genomes (KEGG) and gene ontology (GO) enrichment analyses of the 105 hub genes showed that Dendrobium spp. adapt to cold stress by regulating signal transduction, phytohormones, transcription factors, protein translation and modification, functional proteins, biosynthesis and metabolism, cell structure, light, and the circadian clock. Hub genes of the cold stress response network included the remorin gene pp34, the abscisic acid signaling pathway-related genes PROTEIN PHOSPATASE 2 C (PP2C), SNF1-RELATED PROTEIN KINASE 2 (SnRK2), ABRE-BINDING FACTOR 1 (ABF1) and SKI-INTERACTING PROTEIN 17 (SKIP17), the Ca2+ signaling-related GTP diphosphokinase gene CRSH1, the carbohydrate-related gene STARCH SYNTHASE 2 (SS2), the cell wall biosynthesis gene CINNAMYL ALCOHOL DEHYDROGENASE (CAD7), and the endocytosis-related gene VACUOLAR PROTEIN SORTING-ASSOCIATED PROTEIN 52 A (VPS52A).

CONCLUSIONS

The cold-responsive genes and metabolic pathways of Dendrobium spp. revealed in this study provide important insight to enable the genetic enhancement of cold tolerance in Dendrobium spp., and to facilitate cold tolerance breeding in related plants.

Reducing Acrylamide Formation Potential by Targeting Free Asparagine Accumulation in Seeds

Acrylamide is a probable carcinogen in humans and is formed when reducing sugars react with free asparagine (Asn) during thermal processing of food. Although breeding for low reducing sugars worked well in potatoes, it is less successful in cereals. However, reducing free Asn in cereals has great potential for reducing acrylamide formation, despite the role that Asn plays in nitrogen transport and amino acid biosynthesis. In this perspective, we summarize the efforts aimed at reducing free Asn in cereal grains and discuss the potentials and challenges associated with targeting this essential amino acid, especially in a seed-specific manner.

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.

The conservation of allelic DNA methylation and its relationship with imprinting in maize

Genomic imprinting refers to allele-specific expression of genes depending on parental origin, and it is regulated by epigenetic modifications. Intraspecific allelic variation for imprinting has been detected; however, the intraspecific genome-wide allelic epigenetic variation in maize and its correlation with imprinting variants remain unclear. Here, three reciprocal hybrids were generated by crossing Zea mays inbred lines CAU5, B73, and Mo17 in order to examine the intraspecific conservation of the imprinted genes in the kernel. The majority of imprinted genes exhibited intraspecific conservation, and these genes also exhibited interspecific conservation (rice, sorghum, and Arabidopsis) and were enriched in some specific pathways. By comparing intraspecific allelic DNA methylation in the endosperm, we found that nearly 15% of DNA methylation existed as allelic variants. The intraspecific whole-genome correlation between DNA methylation and imprinted genes indicated that DNA methylation variants play an important role in imprinting variants. Disruption of two conserved imprinted genes using CRISPR/Cas9 editing resulted in a smaller kernel phenotype. Our results shed light on the intraspecific correlation of DNA methylation variants and variation for imprinting in maize, and show that imprinted genes play an important role in kernel development.

United by Conflict: Convergent Signatures of Parental Conflict in Angiosperms and Placental Mammals

Endosperm in angiosperms and placenta in eutherians are convergent innovations for efficient embryonic nutrient transfer. Despite advantages, this reproductive strategy incurs metabolic costs that maternal parents disproportionately shoulder, leading to potential inter-parental conflict over optimal offspring investment. Genomic imprinting-parent-of-origin-biased gene expression-is fundamental for endosperm and placenta development and has convergently evolved in angiosperms and mammals, in part, to resolve parental conflict. Here, we review the mechanisms of genomic imprinting in these taxa. Despite differences in the timing and spatial extent of imprinting, these taxa exhibit remarkable convergence in the molecular machinery and genes governing imprinting. We then assess the role of parental conflict in shaping evolution within angiosperms and eutherians using four criteria: (1) Do differences in the extent of sibling relatedness cause differences in the inferred strength of parental conflict? (2) Do reciprocal crosses between taxa with different inferred histories of parental conflict exhibit parent-of-origin growth effects? (3) Are these parent-of-origin growth effects caused by dosage-sensitive mechanisms and do these loci exhibit signals of positive selection? (4) Can normal development be restored by genomic perturbations that restore stoichiometric balance in the endosperm/placenta? Although we find evidence for all criteria in angiosperms and eutherians, suggesting that parental conflict may help shape their evolution, many questions remain. Additionally, myriad differences between the two taxa suggest that their respective biologies may shape how/when/where/to what extent parental conflict manifests. Lastly, we discuss outstanding questions, highlighting the power of comparative work in quantifying the role of parental conflict in evolution.

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.

NAKED ENDOSPERM1, NAKED ENDOSPERM2, and OPAQUE2 interact to regulate gene networks in maize endosperm development

NAKED ENDOSPERM1 (NKD1), NKD2, and OPAQUE2 (O2) are transcription factors important for cell patterning and nutrient storage in maize (Zea mays) endosperm. To study the complex regulatory interrelationships among these 3 factors in coregulating gene networks, we developed a set of nkd1, nkd2, and o2 homozygous lines, including all combinations of mutant and wild-type genes. Among the 8 genotypes tested, we observed diverse phenotypes and gene interactions affecting cell patterning, starch content, and storage proteins. From ∼8 to ∼16 d after pollination, maize endosperm undergoes a transition from cellular development to nutrient accumulation for grain filling. Gene network analysis showed that NKD1, NKD2, and O2 dynamically regulate a hierarchical gene network during this period, directing cellular development early and then transitioning to constrain cellular development while promoting the biosynthesis and storage of starch, proteins, and lipids. Genetic interactions regulating this network are also dynamic. The assay for transposase-accessible chromatin using sequencing (ATAC-seq) showed that O2 influences the global regulatory landscape, decreasing NKD1 and NKD2 target site accessibility, while NKD1 and NKD2 increase O2 target site accessibility. In summary, interactions of NKD1, NKD2, and O2 dynamically affect the hierarchical gene network and regulatory landscape during the transition from cellular development to grain filling in maize endosperm.

Long Non-Coding RNAs: Discoveries, Mechanisms, and Research Strategies in Seeds

Seeds provide nutrients for the embryo and allow for dormancy in stressed environments to better adapt the plant to its environment. In addition, seeds are an essential source of food for human survival and are the basis for the formation of food production and quality. Therefore, the research on the genetic mechanism of seed development and germination will provide a theoretical basis and technical support for the improvement of crop yield and quality. Recent studies have shown that long non-coding RNAs (lncRNAs) occupy a pivotal position in seed development and germination. In this review, we describe the key processes in seed biology and examine discoveries and insights made in seed lncRNA, with emphasis on lncRNAs that regulate seed biology through multiple mechanisms. Given that thousands of lncRNAs are present in the seed transcriptome, characterization has lagged far behind identification. We provide an overview of research strategies and approaches including some exciting new techniques that may uncover the function of lncRNAs in seed. Finally, we discuss the challenges facing the field and the opening questions. All in all, we hope to provide a clear perspective on discoveries of seed lncRNA by linking discoveries, mechanisms, and technologies.

Auxin and light-mediated regulation of growth, morphogenesis, and alkaloid biosynthesis in Crinum x powellii 'Album' callus

Crinum x powellii 'Album' belongs to the Amaryllidaceae medicinal plant family that produces a range of structurally diverse alkaloids with potential therapeutic properties. The optimal conditions for in vitro tissue growth, morphogenesis, and alkaloid biosynthesis remain unclear. Auxin and light play critical roles in regulating plant growth, development, and alkaloid biosynthesis in several Amaryllidaceae plants. Here, we have succeeded in showing, for the first time, that the combination of auxin and light significantly influence C. x powellii "Album" in vitro tissue growth, survival, and morphogenesis compared to individual treatments. Furthermore, this combination also upregulates the expression of alkaloid biosynthetic genes and led to an increase in the content of certain alkaloids, suggesting a positive impact on the defense and therapeutic potential of the calli. Our findings provide insights into the regulation of genes involved in alkaloid biosynthesis in C. x powellii "Album" callus and underline the potential of auxin and light as tools for enhancing their production in plants. This study provides a foundation for further exploration of C. x powellii "Album" calli as a sustainable source of bioactive alkaloids for pharmaceutical and agricultural applications. Furthermore, this study paves the way to the discovery of the biosynthetic pathway of specialized metabolites from C. x powellii "Album", such as cherylline and lycorine.

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.

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).

Defective kernel 66 encodes a GTPase essential for kernel development in maize

The mitochondrion is a semi-autonomous organelle that provides energy for cell activities through oxidative phosphorylation. In this study, we identified a defective kernel 66 (dek66)-mutant maize with defective kernels. We characterized a candidate gene, DEK66, encoding a ribosomal assembly factor located in mitochondria and possessing GTPase activity (which belongs to the ribosome biogenesis GTPase A family). In the dek66 mutant, impairment of mitochondrial structure and function led to the accumulation of reactive oxygen species and promoted programmed cell death in endosperm cells. Furthermore, the transcript levels of most of the key genes associated with nutrient storage, mitochondrial respiratory chain complex, and mitochondrial ribosomes in the dek66 mutant were significantly altered. Collectively, the results suggest that DEK66 is essential for the development of maize kernels by affecting mitochondrial function. This study provides a reference for understanding the impact of a mitochondrial ribosomal assembly factor in maize kernel development.

Global Transcriptome and Co-Expression Network Analyses Revealed Hub Genes Controlling Seed Size/Weight and/or Oil Content in Peanut

Cultivated peanut (Arachis hypogaea L.) is an important economic and oilseed crop worldwide, providing high-quality edible oil and high protein content. Seed size/weight and oil content are two important determinants of yield and quality in peanut breeding. To identify key regulators controlling these two traits, two peanut cultivars with contrasting phenotypes were compared to each other, one having a larger seed size and higher oil content (Zhonghua16, ZH16 for short), while the second cultivar had smaller-sized seeds and lower oil content (Zhonghua6, ZH6). Whole transcriptome analyses were performed on these two cultivars at four stages of seed development. The results showed that ~40% of the expressed genes were stage-specific in each cultivar during seed development, especially at the early stage of development. In addition, we identified a total of 5356 differentially expressed genes (DEGs) between ZH16 and ZH6 across four development stages. Weighted gene co-expression network analysis (WGCNA) based on DEGs revealed multiple hub genes with potential roles in seed size/weight and/or oil content. These hub genes were mainly involved in transcription factors (TFs), phytohormones, the ubiquitin-proteasome pathway, and fatty acid synthesis. Overall, the candidate genes and co-expression networks detected in this study could be a valuable resource for genetic breeding to improve seed yield and quality traits in peanut.

Laser capture microdissection-based spatiotemporal transcriptomes uncover regulatory networks during seed abortion in seedless Ponkan (Citrus reticulata)

Seed abortion is an important process in the formation of seedless characteristics in citrus fruits. However, the molecular regulatory mechanism underlying citrus seed abortion is poorly understood. Laser capture microdissection-based RNA-seq combined with Pacbio-seq was used to profile seed development in the Ponkan cultivars 'Huagan No. 4' (seedless Ponkan) (Citrus reticulata) and 'E'gan No. 1' (seeded Ponkan) (C. reticulata) in two types of seed tissue across three developmental stages. Through comparative transcriptome and dynamic phytohormone analyses, plant hormone signal, cell division and nutrient metabolism-related processes were revealed to play critical roles in the seed abortion of 'Huagan No. 4'. Moreover, several genes may play indispensable roles in seed abortion of 'Huagan No. 4', such as CrWRKY74, CrWRKY48 and CrMYB3R4. Overexpression of CrWRKY74 in Arabidopsis resulted in severe seed abortion. By analyzing the downstream regulatory network, we further determined that CrWRKY74 participated in seed abortion regulation by inducing abnormal programmed cell death. Of particular importance is that a preliminary model was proposed to depict the regulatory networks underlying seed abortion in citrus. The results of this study provide novel insights into the molecular mechanism across citrus seed development, and reveal the master role of CrWRKY74 in seed abortion of 'Huagan No. 4'.

Comparative transcriptomic analysis of transcription factors and hormones during flower bud differentiation in 'Red Globe' grape under red‒blue light

Grape is a globally significant fruit-bearing crop, and the grape flower bud differentiation essential to fruit production is closely related to light quality. To investigate the regulatory mechanism of grape flower bud differentiation under red‒blue light, the transcriptome and hormone content were determined at four stages of flower bud differentiation. The levels of indole-3-acetic acid (IAA) and abscisic acid (ABA) in grape flower buds at all stages of differentiation under red‒blue light were higher than those in the control. However, the levels of cytokinins (CKs) and gibberellic acid (giberellins, GAs) fluctuated continuously over the course of flower bud differentiation. Moreover, many differentially expressed genes were involved in auxin, CK, GA, and the ABA signal transduction pathways. There were significant differences in the AUX/IAA, SAUR, A-RR, and ABF gene expression levels between the red‒blue light treatment and the control buds, especially in regard to the ABF genes, the expression levels of which were completely different between the two groups. The expression of GBF4 and AI5L2 in the control was always low, while the expression under red‒blue light increased. AI5L7 and AI5L5 expression levels showed an upwards trend in the control plant buds and gradually decreased in red‒blue light treatment plant buds. Through weighted gene coexpression network analysis, we determined that the transcription factors WRK48 (WRKY family), EF110 (ERF family), ABR1, CAMTA3 (CAMTA family), and HSFA3 (HSF family) may be involved in the regulation of the GBF4 gene. This study lays a foundation for further analysis of grape flower bud differentiation regulation under red‒blue light.

Insights into the Functional Components in Wheat Grain: Spatial Pattern, Underlying Mechanism and Cultivation Regulation

Wheat is a staple crop; its production must achieve both high yield and good quality due to worldwide demands for food security and better quality of life. It has been found that the grain qualities vary greatly within the different layers of wheat kernels. In this paper, the spatial distributions of protein and its components, starch, dietary fiber, and microelements are summarized in detail. The underlying mechanisms regarding the formation of protein and starch, as well as spatial distribution, are discussed from the views of substrate supply and the protein and starch synthesis capacity. The regulating effects of cultivation practices on gradients in composition are identified. Finally, breakthrough solutions for exploring the underlying mechanisms of the spatial gradients of functional components are presented. This paper will provide research perspectives for producing wheat that is both high in yield and of good quality.

The translational landscape of bread wheat during grain development

The dynamics of gene expression in crop grains has typically been investigated at the transcriptional level. However, this approach neglects translational regulation, a widespread mechanism that rapidly modulates gene expression to increase the plasticity of organisms. Here, we performed ribosome profiling and polysome profiling to obtain a comprehensive translatome data set of developing bread wheat (Triticum aestivum) grains. We further investigated the genome-wide translational dynamics during grain development, revealing that the translation of many functional genes is modulated in a stage-specific manner. The unbalanced translation between subgenomes is pervasive, which increases the expression flexibility of allohexaploid wheat. In addition, we uncovered widespread previously unannotated translation events, including upstream open reading frames (uORFs), downstream open reading frames (dORFs), and open reading frames (ORFs) in long noncoding RNAs, and characterized the temporal expression dynamics of small ORFs. We demonstrated that uORFs act as cis-regulatory elements that can repress or even enhance the translation of mRNAs. Gene translation may be combinatorially modulated by uORFs, dORFs, and microRNAs. In summary, our study presents a translatomic resource that provides a comprehensive and detailed overview of the translational regulation in developing bread wheat grains. This resource will facilitate future crop improvements for optimal yield and quality.

ZmDRR206 Regulates Nutrient Accumulation in Endosperm through Its Role in Cell Wall Biogenesis during Maize Kernel Development

Dirigent proteins (DIRs) contribute to plant fitness by dynamically reorganizing the cell wall and/or by generating defense compounds during plant growth, development, and interactions with environmental stresses. ZmDRR206 is a maize DIR, it plays a role in maintaining cell wall integrity during seedling growth and defense response in maize, but its role in regulating maize kernel development is unclear. Association analysis of candidate genes indicated that the natural variations of ZmDRR206 were significantly associated with maize hundred-kernel weight (HKW). ZmDRR206 plays a dominant role in storage nutrient accumulation in endosperm during maize kernel development, ZmDRR206 overexpression resulted in small and shrunken maize kernel with significantly reduced starch content and significantly decreased HKW. Cytological characterization of the developing maize kernels revealed that ZmDRR206 overexpression induced dysfunctional basal endosperm transfer layer (BETL) cells, which were shorter with less wall ingrowth, and defense response was constitutively activated in developing maize kernel at 15 and 18 DAP by ZmDRR206 overexpression. The BETL-development-related genes and auxin signal-related genes were down-regulated, while cell wall biogenesis-related genes were up-regulated in developing BETL of the ZmDRR206-overexpressing kernel. Moreover, the developing ZmDRR206-overexpressing kernel had significantly reduced contents of the cell wall components such as cellulose and acid soluble lignin. These results suggest that ZmDRR206 may play a regulatory role in coordinating cell development, storage nutrient metabolism, and stress responses during maize kernel development through its role in cell wall biogenesis and defense response, and provides new insights into understanding the mechanisms of kernel development in maize.

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.

Dissection of Developmental Programs and Regulatory Modules Directing Endosperm Transfer Cell and Aleurone Identity in the Syncytial Endosperm of Barley

Endosperm development in barley starts with the formation of a multinucleate syncytium, followed by cellularization in the ventral part of the syncytium generating endosperm transfer cells (ETCs) as first differentiating subdomain, whereas aleurone (AL) cells will originate from the periphery of the enclosing syncytium. Positional signaling in the syncytial stage determines cell identity in the cereal endosperm. Here, we performed a morphological analysis and employed laser capture microdissection (LCM)-based RNA-seq of the ETC region and the peripheral syncytium at the onset of cellularization to dissect developmental and regulatory programs directing cell specification in the early endosperm. Transcriptome data revealed domain-specific characteristics and identified two-component signaling (TCS) and hormone activities (auxin, ABA, ethylene) with associated transcription factors (TFs) as the main regulatory links for ETC specification. On the contrary, differential hormone signaling (canonical auxin, gibberellins, cytokinin) and interacting TFs control the duration of the syncytial phase and timing of cellularization of AL initials. Domain-specific expression of candidate genes was validated by in situ hybridization and putative protein-protein interactions were confirmed by split-YFP assays. This is the first transcriptome analysis dissecting syncytial subdomains of cereal seeds and provides an essential framework for initial endosperm differentiation in barley, which is likely also valuable for comparative studies with other cereal crops.

Multi-omic characterization of the maize GPI synthesis mutant gwt1 with defects in kernel development

BACKGROUND

Glycosylphosphatidylinositol (GPI) and GPI-anchored proteins (GAPs) are important for cell wall formation and reproductive development in Arabidopsis. However, monocot counterparts that function in kernel endosperm development have yet to be discovered. Here, we performed a multi-omic analysis to explore the function of GPI related genes on kernel development in maize.

RESULTS

In maize, 48 counterparts of human GPI synthesis and lipid remodeling genes were identified, in which null mutation of the glucosaminyl-phosphatidylinositol O-acyltransferase1 gene, ZmGWT1, caused a kernel mutant (named gwt1) with defects in the basal endosperm transport layer (BETL). We performed plasma membrane (PM) proteomics to characterize the potential GAPs involved in kernel development. In total, 4,981 proteins were successfully identified in 10-DAP gwt1 kernels of mutant and wild-type (WT), including 1,638 membrane-anchored proteins with different posttranslational modifications. Forty-seven of the 256 predicted GAPs were differentially accumulated between gwt1 and WT. Two predicted BETL-specific GAPs (Zm00001d018837 and Zm00001d049834), which kept similar abundance at general proteome but with significantly decreased abundance at membrane proteome in gwt1 were highlighted.

CONCLUSIONS

Our results show the importance of GPI and GAPs for endosperm development and provide candidate genes for further investigation of the regulatory network in which ZmGWT1 participates.

A translatome-transcriptome multi-omics gene regulatory network reveals the complicated functional landscape of maize

BACKGROUND

Maize (Zea mays L.) is one of the most important crops worldwide. Although sophisticated maize gene regulatory networks (GRNs) have been constructed for functional genomics and phenotypic dissection, a multi-omics GRN connecting the translatome and transcriptome is lacking, hampering our understanding and exploration of the maize regulatome.

RESULTS

We collect spatio-temporal translatome and transcriptome data and systematically explore the landscape of gene transcription and translation across 33 tissues or developmental stages of maize. Using this comprehensive transcriptome and translatome atlas, we construct a multi-omics GRN integrating mRNAs and translated mRNAs, demonstrating that translatome-related GRNs outperform GRNs solely using transcriptomic data and inter-omics GRNs outperform intra-omics GRNs in most cases. With the aid of the multi-omics GRN, we reconcile some known regulatory networks. We identify a novel transcription factor, ZmGRF6, which is associated with growth. Furthermore, we characterize a function related to drought response for the classic transcription factor ZmMYB31.

CONCLUSIONS

Our findings provide insights into spatio-temporal changes across maize development at both the transcriptome and translatome levels. Multi-omics GRNs represent a useful resource for dissection of the regulatory mechanisms underlying phenotypic variation.

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.

Comparative analysis of the chrysanthemum transcriptome with DNA methylation inhibitors treatment and silencing MET1 lines

BACKGROUND

As one of the ten most famous flowers in China, the chrysanthemum has rich germplasm with a variety of flowering induction pathways, most of which are photoperiod-induced. After treatment with DNA methylation inhibitors, it was found that DNA methylation plays an important role in flowering regulation, but the mechanism of action remains unclear. Therefore, in this study, curcumin, 5-azaC, their mixed treatment, and MET1_-RNAi lines were used for transcriptome sequencing to find out how different treatments affected gene expression in chrysanthemums at different stages of flowering.

RESULTS

Genomic DNA methylation levels were measured using HPLC technology. The methylation level of the whole genome in the vegetative growth stage was higher than that in the flowering stage. The methylation level of DNA in the vegetative growth stage was the lowest in the curcumin and mixed treatment, and the methylation level of DNA in the transgenic line, mixed treatment, and curcumin treatment was the lowest in the flowering stage. The flowering rate of mixed treatment and curcumin treatment was the lowest. Analysis of differentially expressed genes in transcriptomes showed that 5-azaC treatment had the most differentially expressed genes, followed by curcumin and transgenic lines, and mixed treatment had the fewest. In addition, 5-azaC treatment resulted in the differential expression of multiple DNA methylation transferases, which led to the differential expression of many genes. Analysis of differentially expressed genes in different treatments revealed that different treatments had gene specificity. However, the down-regulated GO pathway in all 4 treatments was involved in the negative regulation of the reproductive process, and post-embryonic development, and regulation of flower development. Several genes associated with DNA methylation and flowering regulation showed differential expression in response to various treatments.

CONCLUSIONS

Both DNA methylase reagent treatment and targeted silencing of the MET1 gene can cause differential expression of the genes. The operation of the exogenous application is simple, but the affected genes are exceedingly diverse and untargeted. Therefore, it is possible to construct populations with DNA methylation phenotypic diversity and to screen genes for DNA methylation regulation.

A genome-wide association study identifies a transporter for zinc uploading to maize kernels

The Zn content in cereal seeds is an important trait for crop production as well as for human health. However, little is known about how Zn is loaded to plant seeds. Here, through a genome-wide association study (GWAS), we identify the Zn-NA (nicotianamine) transporter gene ZmYSL2 that is responsible for loading Zn to maize kernels. High promoter sequence variation in ZmYSL2 most likely drives the natural variation in Zn concentrations in maize kernels. ZmYSL2 is specifically localized on the plasma membrane facing the maternal tissue of the basal endosperm transfer cell layer (BETL) and functions in loading Zn-NA into the BETL. Overexpression of ZmYSL2 increases the Zn concentration in the kernels by 31.6%, which achieves the goal of Zn biofortification of maize. These findings resolve the mystery underlying the loading of Zn into plant seeds, providing an efficient strategy for breeding or engineering maize varieties with enriched Zn nutrition.

Coexpression Network Construction and Visualization from Transcriptomes Underlying ER Stress Responses

Dynamic gene expression changes are primary cellular reactions in response to most stresses and developmental cues in all organisms, including plants. With the ever-decreasing cost and increasing access, high-throughput transcriptome analyses have become a significant research tool to understand a wide spectrum of complex gene regulatory mechanisms. However, it is still challenging to understand the complete picture of gene responses because of the interactive and dynamic nature of gene expression in biological networks. Coexpression network analyses followed by network mapping are being increasingly applied to overcome this challenge. In this chapter, we will introduce detailed instructions for performing a weighted coexpression network analysis (WGCNA) and network visualization using a transcriptome dataset obtained during recovery from endoplasmic reticulum (ER) stress in Arabidopsis thaliana. The streamlined workflow described here allows biologists to identify and visualize coexpression interactions among genes, accessing a comprehensive landscape of dynamic gene expression changes for further downstream analyses using their datasets.

Protoplast Isolation for Plant Single-Cell RNA-seq

The growth and development of plants depends on diversified gene expression in different cell types. Compared to traditional bulk RNA sequencing, droplet-based single-cell RNA sequencing (scRNA-seq) allows for transcriptome profiling of individual cells within heterogeneous tissues. scRNA-seq provides a high-resolution atlas of cellular characterization and vastly improves our understandings of the interactions between individual cells and the microenvironment. However, the difficulty in protoplast isolation has limited the application of single-cell sequencing technology in plant research. Here we describe a high-efficiency protoplast isolation protocol for scRNA-seq.

Coexpression network and trans-activation analyses of maize reproductive phasiRNA loci

The anther-enriched phased, small interfering RNAs (phasiRNAs) play vital roles in sustaining male fertility in grass species. Their long non-coding precursors are synthesized by RNA polymerase II and are likely regulated by transcription factors (TFs). A few putative transcriptional regulators of the 21- or 24-nucleotide phasiRNA loci (referred to as 21- or 24-PHAS loci) have been identified in maize (Zea mays), but whether any of the individual TFs or TF combinations suffice to activate any PHAS locus is unclear. Here, we identified the temporal gene coexpression networks (modules) associated with maize anther development, including two modules highly enriched for the 21- or 24-PHAS loci. Comparisons of these coexpression modules and gene sets dysregulated in several reported male sterile TF mutants provided insights into TF timing with regard to phasiRNA biogenesis, including antagonistic roles for OUTER CELL LAYER4 and MALE STERILE23. Trans-activation assays in maize protoplasts of individual TFs using bulk-protoplast RNA-sequencing showed that two of the TFs coexpressed with 21-PHAS loci could activate several 21-nucleotide phasiRNA pathway genes but not transcription of 21-PHAS loci. Screens for combinatorial activities of these TFs and, separately, the recently reported putative transcriptional regulators of 24-PHAS loci using single-cell (protoplast) RNA-sequencing, did not detect reproducible activation of either 21-PHAS or 24-PHAS loci. Collectively, our results suggest that the endogenous transcriptional machineries and/or chromatin states in the anthers are necessary to activate reproductive PHAS loci.

Comparative Transcriptome and Co-Expression Network Analyses Reveal the Molecular Mechanism of Calcium-Deficiency-Triggered Tipburn in Chinese Cabbage (Brassica rapa L. ssp. Pekinensis)

Chinese cabbage tipburn is characterized by the formation of necrotic lesions on the margin of leaves, including on the insides of the leafy head. This physiological disorder is associated with a localized calcium deficiency during leaf development. However, little information is available regarding the molecular mechanisms governing Ca-deficiency-triggered tipburn. This study comprehensively analysed the transcriptomic comparison between control and calcium treatments (CK and 0 mM Ca) in Chinese cabbage to determine its molecular mechanism in tipburn. Our analysis identified that the most enriched gene ontology (GO) categories are photosynthesis, thylakoid and cofactor binding. Moreover, the KEGG pathway was most enriched in photosynthesis, carbon metabolism and carbon fixation. We also analyzed the co-expression network by functional categories and identified ten critical hub differentially expressed genes (DEGs) in each gene regulatory network (GRN). These DEGs might involve abiotic stresses, developmental processes, cell wall metabolism, calcium distribution, transcription factors, plant hormone biosynthesis and signal transduction pathways. Under calcium deficiency, CNX1, calmodulin-binding proteins and CMLs family proteins were downregulated compared to CK. In addition, plant hormones such as GA, JA, BR, Auxin and ABA biosynthesis pathways genes were downregulated under calcium treatment. Likewise, HATs, ARLs and TCP transcription factors were reported as inactive under calcium deficiency, and potentially involved in the developmental process. This work explores the specific DEGs' significantly different expression levels in 0 mM Ca and the control involved in plant hormones, cell wall developments, a light response such as chlorophylls and photosynthesis, transport metabolism and defence mechanism and redox. Our results provide critical evidence of the potential roles of the calcium signal transduction pathway and candidate genes governing Ca-deficiency-triggered tipburn in Chinese cabbage.

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.

Developmental evidence for parental conflict in driving Mimulus species barriers

The endosperm, a tissue that nourishes the embryo in the seeds of flowering plants, is often disrupted in inviable hybrid seeds of closely related species. A key question is whether parental conflict is a major driver of this common form of reproductive isolation. Here, we performed reciprocal crosses between pairs of three monkeyflower species (Mimulus caespitosa, Mimulus tilingii, and Mimulus guttatus). The severity of hybrid seed inviability varies among these crosses, which we inferred to be due to species divergence in effective ploidy. By performing a time series experiment of seed development, we discovered parent-of-origin phenotypes that provide strong evidence for parental conflict in shaping endosperm evolution. We found that the chalazal haustorium, a tissue within the endosperm that is found at the maternal-filial boundary, shows pronounced differences between reciprocal hybrid seeds formed from Mimulus species that differ in effective ploidy. These parent-of-origin effects suggest that the chalazal haustorium might act as a mediator of parental conflict, potentially by controlling sucrose movement from the maternal parent into the endosperm. Our study suggests that parental conflict in the endosperm may function as a driver of speciation by targeting regions and developmental stages critical for resource allocation and thus proper seed development.

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.

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

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

Current studies of maize kernel development mostly focus on endosperm and embryo development, and little is known about the control of the nucellus tissue. Here, the authors report a NAC-EXPB15 pathway that regulate maize kernel size by modulating nucellus elimination and early endosperm development.

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.

Rearrangement with the nkd2 promoter contributed to allelic diversity of the r1 gene in maize (Zea mays)

The maize red1 (r1) locus regulates anthocyanin accumulation and is a classic model for allelic diversity; changes in regulatory regions are responsible for most of the variation in gene expression patterns. Here, an intrachromosomal rearrangement between the distal upstream region of r1 and the region of naked endosperm 2 (nkd2) upstream to the third exon generated a nkd2 null allele lacking the first three exons, and the R1-st (stippled) allele with a novel r1 5' promoter region homologous to 5' regions from nkd2-B73. R1-sc:124 (an R1-st derivative) shows increased and earlier expression than a standard R1-g allele, as well as ectopic expression in the starchy endosperm compartment. Laser capture microdissection and RNA sequencing indicated that ectopic R1-sc:124 expression impacted expression of genes associated with RNA modification. The expression of R1-sc:124 resembled nkd2-W22 expression, suggesting that nkd2 regulatory sequences may influence the expression of R1-sc:124. The r1-sc:m3 allele is derived from R1-sc:124 by an insertion of a Ds6 transposon in intron 4. This insertion blocks anthocyanin regulation by causing mis-splicing that eliminates exon 5 from the mRNA. This allele serves as an important launch site for Ac/Ds mutagenesis studies, and two Ds6 insertions believed to be associated with nkd2 mutant alleles were actually located in the r1 5' region. Among annotated genomes of teosinte and maize varieties, the nkd2 and r1 loci showed conserved overall gene structures, similar to the B73 reference genome, suggesting that the nkd2-r1 rearrangement may be a recent event.

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.

Regulators of early maize leaf development inferred from transcriptomes of laser capture microdissection (LCM)-isolated embryonic leaf cells

The superior photosynthetic efficiency of C₄ leaves over C₃ leaves is owing to their unique Kranz anatomy, in which the vein is surrounded by one layer of bundle sheath (BS) cells and one layer of mesophyll (M) cells. Kranz anatomy development starts from three contiguous ground meristem (GM) cells, but its regulators and underlying molecular mechanism are largely unknown. To identify the regulators, we obtained the transcriptomes of 11 maize embryonic leaf cell types from five stages of pre-Kranz cells starting from median GM cells and six stages of pre-M cells starting from undifferentiated cells. Principal component and clustering analyses of transcriptomic data revealed rapid pre-Kranz cell differentiation in the first two stages but slow differentiation in the last three stages, suggesting early Kranz cell fate determination. In contrast, pre-M cells exhibit a more prolonged transcriptional differentiation process. Differential gene expression and coexpression analyses identified gene coexpression modules, one of which included 3 auxin transporter and 18 transcription factor (TF) genes, including known regulators of Kranz anatomy and/or vascular development. In situ hybridization of 11 TF genes validated their expression in early Kranz development. We determined the binding motifs of 15 TFs, predicted TF target gene relationships among the 18 TF and 3 auxin transporter genes, and validated 67 predictions by electrophoresis mobility shift assay. From these data, we constructed a gene regulatory network for Kranz development. Our study sheds light on the regulation of early maize leaf development and provides candidate leaf development regulators for future study.

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.

Transcriptomic analysis of genes related to alkaloid biosynthesis and the regulation mechanism under precursor and methyl jasmonate treatment in Dendrobium officinale

Dendrobium officinale is both a traditional herbal medicine and a plant of high ornamental and medicinal value. Alkaloids, especially terpenoid indole alkaloids (TIAs), with pharmacological activities are present in the tissues of D. officinale. A number of genes involved in alkaloid biosynthetic pathways have been identified. However, the regulatory mechanisms underlying the precursor and methyl jasmonate (MeJA)-induced accumulation of alkaloids in D. officinale are poorly understood. In this study, we collected D. officinale protocorm-like bodies (PLBs) and treated them with TIA precursors (tryptophan and secologanin) and MeJA for 0 (T0), 4 (T4) and 24 h (T24); we also established control samples (C4 and C24). Then, we measured the total alkaloid content of the PLBs and performed transcriptome sequencing using the Illumina HiSeq 2,500 system. The total alkaloid content increased significantly after 4 h of treatment. Go and KEGG analysis suggested that genes from the TIA, isoquinoline alkaloid, tropane alkaloid and jasmonate (JA) biosynthetic pathways were significantly enriched. Weighted gene coexpression network analysis (WGCNA) uncovered brown module related to alkaloid content. Six and seven genes related to alkaloid and JA bisosynthetic pathways, respectively, might encode the key enzymes involved in alkaloid biosynthesis of D. officinale. Moreover, 13 transcription factors (TFs), which mostly belong to AP2/ERF, WRKY, and MYB gene families, were predicted to regulate alkaloid biosynthesis. Our data provide insight for studying the regulatory mechanism underlying TIA precursor and MeJA-induced accumulation of three types of alkaloids in D. officinale.