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

The aleurone layer of cereal grains: Development, genetic regulation, and breeding applications

Cereal aleurone cells are differentiated from triploid endosperm cells and exhibit distinct cytological, physiological, and biochemical characteristics that distinguish them from the starchy endosperm cells of cereals. Aleurone cells maintain viability throughout seed development, whereas starchy endosperm cells undergo programmed cell death during maturation. Despite variations in aleurone-related traits among cereal species, the aleurone layer plays a crucial role in regulating many aspects of seed development, including the accumulation of storage reserves, the acquisition of dormancy, and germination. Given that many nutrients-such as lipids, dietary fibers, vitamins, and minerals like iron and zinc-are predominantly accumulated in the aleurone cells of cereal grains, this layer has attracted considerable attention aimed at improving the nutritional value of cereals. This review provides a comprehensive overview of the developmental, genetic, and molecular basis of aleurone cell differentiation and proliferation. It focuses on the improvement of aleurone-related traits informed by knowledge of the molecular networks governing aleurone development and presents a detailed discussion on the challenges and potential solutions associated with cereal improvement through the manipulation of aleurone-related traits.

Identification of candidate genes in sesquiterpenoid biosynthesis of Atractylodes lancea through combined metabolomic and transcriptomic analysis

Atractylodes lancea is a type of traditional Chinese medicine whose rhizome is rich in sesquiterpenoids. However, the mechanisms underlying the synthesis and metabolism of sesquiterpenoids in A. lancea remain poorly discovered. In the present research, the organs of A. lancea in different growth and development stages served as test materials. Initially, widely targeted metabolomic and transcriptomic analyses were integrated to elucidate sesquiterpenoid biosynthesis in A. lancea. Up to 177 differentially expressed sesquiterpene metabolites and 80 genes associated with the synthesis pathway of sesquiterpenes were identified. Of these, 13 terpene synthases (TPSs), including AlTPS17, AlTPS20, and AlTPS13, appeared to be directly involved in the synthesis of sesquiterpenes accumulated in the rhizomes. Both AlTPS13 and AlTPS20 used farnesyl diphosphate as a substrate to generate sesquiterpenoids, while AlTPS13 catalyzed the biosynthesis of γ-elemene. This study analyzed the quality-generation pattern of A. lancea under seasonal variations and thereby laid a strong foundation to further elucidate the mechanism of sesquiterpenoid biosynthesis.

Transcriptome analysis reveals potential medicinal ingredient synthesis in ornamental Dendrobium

Dendrobium is divided into ornamental and medicinal varieties due to ornamental and medicinal values. However, current research mainly focuses on medicinal Dendrobium, with less study on the medicinal value of ornamental Dendrobium. We analyzed the microstructures, active components of the stems from twelve ornamental Dendrobium, and explored the biosynthetic networks of these active components based on transcriptome sequencing. This study found the Dendrobium with the highest content of polysaccharide, alkaloid, and flavonoid was Dendrobium aphyllum (53.89 %), Dendrobium thyrsiflorum (2.11 %) and Dendrobium loddigesii (7.21 %). Further research revealed 9 DEGs associated with polysaccharide biosynthesis were highly expressed in D. aphyllum; 4 DEGs related to alkaloid biosynthesis were highly expressed in D. thyrsiflorum; 8 DEGs associated with flavonoid biosynthesis were highly expressed in D. loddigesii. This study revealed the potential medicinal value of ornamental Dendrobium and the synthetic mechanisms of its medicinal components, providing a foundation for the medical applications of ornamental Dendrobium.

The chromatin accessibility landscape during early maize seed development

Cis-regulatory elements (CREs) are enriched in accessible chromatin regions (ACRs) of eukaryotes. Despite extensive research on genome-wide ACRs in various plant tissues, the global impact of these changes on developmental processes in maize seeds remains poorly understood. In this study, we employed the assay for transposase-accessible chromatin sequencing (ATAC-seq) to reveal the chromatin accessibility profile throughout the genome during the early stages of maize seed development. We identified a total of 37 952 to 59 887 high-quality ACRs in maize seeds at 0 to 8 days after pollination (DAP). Furthermore, we examined the correlation between the identified ACRs and gene expression. We observed a positive correlation between the open degree of promoter-ACRs and the expression of most genes. Moreover, we identified binding footprints of numerous transcription factors (TFs) within chromatin accessibility regions and revealed key TF families involved in different stages. Through the footprints of accessible chromatin regions, we predicted transcription factor regulatory networks during early maize embryo development. Additionally, we discovered that DNA sequence diversity was notably reduced at ACRs, yet trait-associated SNPs were more likely to be located within ACRs. We edited the ACR containing the trait-associated SNP of NKD1. Both NKD1pro-1 and NKD1pro-2 showed phenotypes corresponding to the trait-associated SNP. Our results suggest that alterations in chromatin accessibility play a crucial role in maize seed development and highlight the potential contribution of open chromatin regions to advancements in maize breeding.

Multiplexed transcriptomic analyzes of the plant embryonic hourglass

Zoologists have adduced morphological convergence among embryonic stages of closely related taxa, which has been called the phylotypic stage of embryogenesis. Transcriptomic analyzes reveal an hourglass pattern of gene expression during plant and animal embryogenesis, characterized by the accumulation of evolutionarily older and conserved transcripts during mid-embryogenesis, whereas younger less-conserved transcripts predominate at earlier and later embryonic stages. In contrast, comparisons of embryonic gene expression among different animal phyla describe an inverse hourglass pattern, where expression is correlated during early and late stages but not during mid-embryo development. Here, multiplexed spatial-transcriptomic analyzes is used to investigate embryogenesis and homology in maize, which has grass-specific morphology. A set of shared, co-expressed genes is identified during initiation of maize embryonic organs, replete for ancient/conserved genes manifesting an hourglass pattern during mid-embryogenesis. Transcriptomic comparisons of maize and Arabidopsis embryogenesis with that of the moss Physcomitrium patens identify an inverse hourglass pattern across plant phyla, as in animals. The data suggest that the phylotypic stages in plants and animals are characterized by expression of ancient and conserved genes during histogenesis, organization of embryonic axes, and initial morphogenesis. We propose a mechanism for gene evolution during the innovation of morphological novelty.

Dissecting the cellular architecture and genetic circuitry of the soybean seed

Seeds are complex structures composed of three regions, embryo, endosperm, and seed coat, with each further divided into subregions that consist of tissues, cell layers, and cell types. Although the seed is well characterized anatomically, much less is known about the genetic circuitry that dictates its spatial complexity. To address this issue, we profiled mRNAs from anatomically distinct seed subregions at several developmental stages. Analyses of these profiles showed that all subregions express similar diverse gene numbers and that the small gene numbers expressed subregion specifically provide information about the biological processes that occur in these seed compartments. In parallel, we profiled RNAs in individual nuclei and identified nuclei clusters representing distinct cell identities. Integrating single-nucleus RNA and subregion mRNA transcriptomes allowed most cell identities to be assigned to specific subregions and cell types and/or cell states. The number of cell identities exceeds the number of anatomically distinguishable cell types, emphasizing the spatial complexity of seeds. We defined gene coexpression networks that underlie distinct biological processes during seed development. We showed that network distribution among subregions and cell identities is highly variable. Some networks operate in single subregions and/or cell identities, and many coexpression networks operate in multiple subregions and/or cell identities. We also showed that single subregions and cell identities possess several networks. Together, our studies provide unique insights into the biological processes and genetic circuitry that underlie the spatial landscape of the seed.

Maize unstable factor for orange1 encodes a nuclear protein that affects redox accumulation during kernel development

The basal endosperm transfer layer (BETL) of the maize (Zea mays L.) kernel is composed of transfer cells for nutrient transport to nourish the developing kernel. To understand the spatiotemporal processes required for BETL development, we characterized 2 unstable factor for orange1 (Zmufo1) mutant alleles. The BETL defects in these mutants were associated with high levels of reactive oxygen species, oxidative DNA damage, and cell death. Interestingly, antioxidant supplementation in in vitro cultured kernels alleviated the cellular defects in mutants. Transcriptome analysis of the loss-of-function Zmufo1 allele showed differential expression of tricarboxylic acid cycle, redox homeostasis, and BETL-related genes. The basal endosperms of the mutant alleles had high levels of acetyl-CoA and elevated histone acetyltransferase activity. The BETL cell nuclei showed reduced electron-dense regions, indicating sparse heterochromatin distribution in the mutants compared with wild-type. Zmufo1 overexpression further reduced histone methylation marks in the enhancer and gene body regions of the pericarp color1 (Zmp1) reporter gene. Zmufo1 encodes an intrinsically disordered nuclear protein with very low sequence similarity to known proteins. Yeast two-hybrid and luciferase complementation assays established that ZmUFO1 interacts with proteins that play a role in chromatin remodeling, nuclear transport, and transcriptional regulation. This study establishes the critical function of Zmufo1 during basal endosperm development in maize kernels.

A comprehensive atlas of endogenous peptides in maize

In this study, we present a comprehensive peptidomic atlas of 13 maize tissues, covering both vegetative and reproductive phases. Using a three-frame translation of canonical coding sequences, we identified 6100 nonredundant endogenous peptides, significantly expanding the known plant peptide repertoire. By integrating peptidomic coexpression profiles with previously reported proteomic profiles, we found that the peptide abundance did not consistently correlate with the abundance of their source proteins, suggesting the presence of complex regulatory mechanisms. This integrated peptidomic and proteomic map can serve as a valuable resource for exploring the functional roles of endogenous peptides in maize development and facilitates the investigation of the functional relationship among genes, peptides, and proteins across various biological contexts.

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.

Dissecting the temporal genetic networks programming soybean embryo development from embryonic morphogenesis to post-germination

KEY MESSAGE

Desiccation-stage transcription factors perform similar functions, with early ones focused on desiccation tolerance and later ones on development. Gene networks governing late embryo development diverge between soybean and Arabidopsis. To understand gene activities programming seed embryo development, we profiled the soybean embryo transcriptome across embryonic morphogenesis through post-germination. Transcriptomic landscapes across embryo development feature highly prevalent transcripts, categorized into early and late groups, with shared and distinct functions. During the mid-storage reserve accumulation stage, the upregulated genes are enriched with regulatory tasks at both the transcriptional and chromatin levels, including DNA methylation and chromatin remodeling. The epigenetic-related functions also dominate in the upregulated genes during germination, involving core histone variants and histone chaperones. Gene network analysis reveals both stage-specific modules and modules active across multiple stages. The desiccation-associated gene module integrates diverse transcription factors (TFs) that are sequentially active during different desiccation stages, transitioning from abiotic stress functions early on to developmental functions later. Two TFs, active during the early and mid-desiccation stages were functionally assessed in Arabidopsis overexpression lines to uncover their potential roles in desiccation processes. Interestingly, nearly half of the Arabidopsis orthologs of soybean TFs active in the desiccation-associated module are inactive during Arabidopsis desiccation. Our results reveal that chromatin and transcriptional regulation coordinate during mid-storage reserve accumulation, while distinct epigenetic mechanisms drive germination. Additionally, gene modules either perform stage-specific functions or are required across multiple stages, and gene networks during late embryogenesis diverge between soybean and Arabidopsis. Our studies provide new information on the biological processes and gene networks underlying development from embryonic morphogenesis to post-germination.

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.

Transcriptomic and coexpression network analyses revealed the regulatory mechanism of Cydia pomonella infestation on the synthesis of phytohormones in walnut husks

The codling moth (Cydia pomonella) has a major effect on the quality and yield of walnut fruit. Plant defences respond to insect infestation by activating hormonal signalling and the flavonoid biosynthetic pathway. However, little is known about the role of walnut husk hormones and flavonoid biosynthesis in response to C. pomonella infestation. The phytohormone content assay revealed that the contents of salicylic acid (SA), abscisic acid (ABA), jasmonic acid (JA), jasmonic acid-isoleucine conjugate (JA-ILE), jasmonic acid-valine (JA-Val) and methyl jasmonate (MeJA) increased after feeding at different time points (0, 12, 24, 36, 48, and 72 h) of walnut husk. RNA-seq analysis of walnut husks following C. pomonella feeding revealed a temporal pattern in differentially expressed genes (DEGs), with the number increasing from 3,988 at 12 h to 5,929 at 72 h postfeeding compared with the control at 0 h postfeeding. Walnut husks exhibited significant upregulation of genes involved in various defence pathways, including flavonoid biosynthesis (PAL, CYP73A, 4CL, CHS, CHI, F3H, ANS, and LAR), SA (PAL), ABA (ZEP and ABA2), and JA (AOS, AOC, OPR, JAZ, and MYC2) pathways. Three gene coexpression networks that had a significant positive association with these hormonal changes were constructed based on the basis of weighted gene coexpression network analysis (WGCNA). We identified several hub transcription factors, including the turquoise module (AIL6, MYB4, PRE6, WRKY71, WRKY31, ERF003, and WRKY75), the green module (bHLH79, PCL1, APRR5, ABI5, and ILR3), and the magenta module (ERF27, bHLH35, bHLH18, TIFY5A, WRKY31, and MYB44). Taken together, these findings provide useful genetic resources for exploring the defence response mediated by phytohormones in walnut husks.

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.

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.

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.

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.

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.