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

Histologically resolved multiomics enables precise molecular profiling of human intratumor heterogeneity

Both the composition of cell types and their spatial distribution in a tissue play a critical role in cellular function, organ development, and disease progression. For example, intratumor heterogeneity and the distribution of transcriptional and genetic events in single cells drive the genesis and development of cancer. However, it can be challenging to fully characterize the molecular profile of cells in a tissue with high spatial resolution because microscopy has limited ability to extract comprehensive genomic information, and the spatial resolution of genomic techniques tends to be limited by dissection. There is a growing need for tools that can be used to explore the relationship between histological features, gene expression patterns, and spatially correlated genomic alterations in healthy and diseased tissue samples. Here, we present a technique that combines label-free histology with spatially resolved multiomics in unfixed and unstained tissue sections. This approach leverages stimulated Raman scattering microscopy to provide chemical contrast that reveals histological tissue architecture, allowing for high-resolution in situ laser microdissection of regions of interests. These microtissue samples are then processed for DNA and RNA sequencing to identify unique genetic profiles that correspond to distinct anatomical regions. We demonstrate the capabilities of this technique by mapping gene expression and copy number alterations to histologically defined regions in human oral squamous cell carcinoma (OSCC). Our approach provides complementary insights in tumorigenesis and offers an integrative tool for macroscale cancer tissues with spatial multiomics assessments.

The pho1;2a'-m1.1 allele of Phosphate1 conditions misregulation of the phosphorus starvation response in maize (Zea mays ssp. mays L.)

Plant PHO1 proteins play a central role in the translocation and sensing of inorganic phosphate. The maize (Zea mays ssp. mays) genome encodes two co-orthologs of the Arabidopsis PHO1 gene, designated ZmPho1;2a and ZmPho1;2b. Here, we report the characterization of the transposon footprint allele Zmpho1;2a'-m1.1, which we refer to hereafter as pho1;2a. The pho1;2a allele is a stable derivative formed by excision of an Activator transposable element from the ZmPho1;2a gene. The pho1;2a allele contains an 8-bp insertion at the point of transposon excision that disrupts the reading frame and is predicted to generate a premature translational stop. We show that the pho1;2a allele is linked to a dosage-dependent reduction in Pho1;2a transcript accumulation and a mild reduction in seedling growth. Characterization of shoot and root transcriptomes under full nutrient, low nitrogen, low phosphorus, and combined low nitrogen and low phosphorus conditions identified 1100 differentially expressed genes between wild-type plants and plants carrying the pho1;2a mutation. Of these 1100 genes, 966 were upregulated in plants carrying pho1;2a, indicating the wild-type PHO1;2a to predominantly impact negative gene regulation. Gene set enrichment analysis of the pho1;2a-misregulated genes revealed associations with phytohormone signaling and the phosphate starvation response. In roots, differential expression was broadly consistent across all nutrient conditions. In leaves, differential expression was largely specific to low phosphorus and combined low nitrogen and low phosphorus conditions. Of 276 genes upregulated in the leaves of pho1;2a mutants in the low phosphorus condition, 153 were themselves induced in wild-type plants with respect to the full nutrient condition. Our observations suggest that Pho1;2a functions in the fine-tuning of the transcriptional response to phosphate starvation through maintenance and/or sensing of plant phosphate status.

Physiological and omics analysis of maize inbred lines during late grain development

BACKGROUND

There were significant differences in the change of moisture content and grain composition at the late stage of grain development among different maize varieties, but the regulation mechanism is not clear.

OBJECTIVE

To explore the key genes causing the variation in physiological traits of two typical maize inbred lines in late grain development.

METHODS

The grains at different development stages were selected as materials to determine the content of water, sucrose, starch and ABA. Transcriptomic and proteomic analysis of the materials were performed to screen relevant genes.

RESULTS

The grain dehydration rate and the content of sucrose, starch and ABA were showed significant differences between two varieties in the late stage of grain development. The enrichment analysis of common differentially expressed genes (proteins) showed that most of the genes (proteins) were enriched in the extracellular region. The downregulated genes were mainly concentrated in carbohydrate metabolism and lipid metabolism, while the upregulated genes were mainly in response to stress. Furthermore, this study also identified many key candidate genes (dehydrin genes, pathogenesis-related genes, sucrose synthase and secondary metabolites related genes) related to late grain development of maize.

CONCLUSIONS

The suggested genes related to late grain development of maize can be candidates for further functional study.

Advances and applications of single-cell omics technologies in plant research

Single-cell sequencing approaches reveal the intracellular dynamics of individual cells and answer biological questions with high-dimensional catalogs of millions of cells, including genomics, transcriptomics, chromatin accessibility, epigenomics, and proteomics data across species. These emerging yet thriving technologies have been fully embraced by the field of plant biology, with a constantly expanding portfolio of applications. Here, we introduce the current technical advances used for single-cell omics, especially single-cell genome and transcriptome sequencing. Firstly, we overview methods for protoplast and nucleus isolation and genome and transcriptome amplification. Subsequently, we use well-executed benchmarking studies to highlight advances made through the application of single-cell omics techniques. Looking forward, we offer a glimpse of additional hurdles and future opportunities that will introduce broad adoption of single-cell sequencing with revolutionary perspectives in plant biology.

Deciphering the Molecular Signatures Associated With Resistance to Botrytis cinerea in Strawberry Flower by Comparative and Dynamic Transcriptome Analysis

Gray mold caused by Botrytis cinerea, which is considered to be the second most destructive necrotrophic fungus, leads to major economic losses in strawberry (Fragaria × ananassa) production. B. cinerea preferentially infects strawberry flowers and fruits, leading to flower blight and fruit rot. Compared with those of the fruit, the mechanisms of flower defense against B. cinerea remain largely unexplored. Therefore, in this study, we aimed to unveil the resistance mechanisms of strawberry flower through dynamic and comparative transcriptome analysis with resistant and susceptible strawberry cultivars. Our experimental data suggest that resistance to B. cinerea in the strawberry flower is probably regulated at the transcriptome level during the early stages of infection and strawberry flower has highly complex and dynamic regulatory networks controlling a multi-layered defense response to B. cinerea. First of all, the higher expression of disease-resistance genes but lower expression of cell wall degrading enzymes and peroxidases leads to higher resistance to B. cinerea in the resistant cultivar. Interestingly, CPKs, RBOHDs, CNGCs, and CMLs comprised a calcium signaling pathway especially play a crucial role in enhancing resistance by increasing their expression. Besides, six types of phytohormones forming a complex regulatory network mediated flower resistance, especially JA and auxin. Finally, the genes involved in the phenylpropanoid and amino acids biosynthesis pathways were gene sets specially expressed or different expression genes, both of them contribute to the flower resistance to B. cinerea. These data provide the foundation for a better understanding of strawberry gray mold, along with detailed genetic information and resistant materials to enable genetic improvement of strawberry plant resistance to gray mold.

Molecular mechanisms of maize endosperm transfer cell development

Endosperm transfer cells function as the nutrient transporter, antimicrobic barrier, and signal mediator between filial and maternal tissues. Sugar supply of maternal tissues, sugar demand of filial tissues, and requirement for defence against pathogens are three elemental factors inducing differentiation of endosperm transfer cells. Epigenetic factors, especially MEG1, moderate the key genetic factor ZmMRP-1 to activate endosperm transfer cell-specific genes that control the flange wall ingrowth formation and defensin-like protein secretion in maize. Auxin and cytokinin are primary hormones involved in development of maize endosperm transfer cells. Crosstalk between glucose and hormone signaling regulates endosperm transfer cell development via modifying ZmMRP-1 expression. This review summarizes the current knowledge on maize endosperm transfer cell development, and discusses its potential molecular mechanisms. It is expected to strengthen the theoretical basis for structural and functional optimization of endosperm transfer cells, and yield improvement of kernels in maize.

Maize RNA 3'-terminal phosphate cyclase-like protein promotes 18S pre-rRNA cleavage and is important for kernel development

Plant ribosomes contain four specialized ribonucleic acids, the 5S, 5.8S, 18S, and 25S ribosomal RNAs (rRNAs). Maturation of the latter three rRNAs requires cooperative processing of a single transcript by several endonucleases and exonucleases at specific sites. In maize (Zea mays), the exact nucleases and components required for rRNA processing remain poorly understood. Here, we characterized a conserved RNA 3'-terminal phosphate cyclase (RCL)-like protein, RCL1, that functions in 18S rRNA maturation. RCL1 is highly expressed in the embryo and endosperm during early seed development. Loss of RCL1 function resulted in lethality due to aborted embryo cell differentiation. We also observed pleiotropic defects in the rcl1 endosperm, including abnormal basal transfer cell layer growth and aleurone cell identity, and reduced storage reserve accumulation. The rcl1 seeds had lower levels of mature 18S rRNA and the related precursors were altered in abundance compared with wild type. Analysis of transcript levels and protein accumulation in rcl1 revealed that the observed lower levels of zein and starch synthesis enzymes mainly resulted from effects at the transcriptional and translational levels, respectively. These results demonstrate that RCL1-mediated 18S pre-rRNA processing is essential for ribosome function and messenger RNA translation during maize seed development.

Exploring the Developmental Progression of Endosperm Cavity Formation in Maize Grain and the Underlying Molecular Basis Using X-Ray Tomography and Genome Wide Association Study

Endosperm cavity (EC) in maize grain reduces yield and causes grain breakage during mechanical harvesting, hence representing a major problem in the maize industry. Despite this, little is known regarding the biological processes governing EC formation. Here, we attempted to address this issue by (i) determining the spatial and temporal progression of EC in a non-invasive manner and (ii) identifying candidate genes that may be linked to the formation of EC by using a genome wide association study (GWAS). Visualization and measurement using X-ray micro-computed tomography established that EC first appeared at the central starch endosperm at about 12 days after pollination (DAP) and became enlarged thereafter. GWAS-based screening of a panel of 299 inbred lines with a wide range of EC size identified nine candidate genes that showed significant association with EC formation. Most of the candidate genes exhibited a decrease at 12 DAP, coinciding with the timing of EC appearance. Among them, ZmMrp11 was annotated as a member encoding a multidrug resistance-associated protein that has been shown in other studies to sequestrate toxic metabolites from the cytosol to the vacuole, thereby detoxifying the cellular environment. This, together with the reduced expression of ZmMrp11 in maize grains from 12 DAP, prompted us to propose that the low expression of ZmMrp11 may block cellular detoxification in the maize endosperm cells, leading to cell death and ultimately the formation of EC.

A transcriptional landscape underlying sugar import for grain set in maize

Developing seed depends on sugar supply for its growth and yield formation. Maize (Zea mays L.) produces the largest grains among cereals. However, there is a lack of holistic understanding of the transcriptional landscape of genes controlling sucrose transport to, and utilization within, maize grains. By performing in-depth data mining of spatio-temporal transcriptomes coupled with histological and heterologous functional analyses, we identified transporter genes specifically expressed in the maternal-filial interface, including (i) ZmSWEET11/13b in the placento-chalazal zone, where sucrose is exported into the apoplasmic space, and (ii) ZmSTP3, ZmSWEET3a/4c (monosaccharide transporters), ZmSUT1, and ZmSWEET11/13a (sucrose transporters) in the basal endosperm transfer cells for retrieval of apoplasmic sucrose or hexoses after hydrolysis by extracellular invertase. In the embryo and its surrounding regions, an embryo-localized ZmSUT4 and a cohort of ZmSWEETs were specifically expressed. Interestingly, drought repressed those ZmSWEETs likely exporting sucrose but enhanced the expression of most transporter genes for uptake of apoplasmic sugars. Importantly, this drought-induced fluctuation in gene expression was largely attenuated by an increased C supply via controlled pollination, indicating that the altered gene expression is conditioned by C availability. Based on the analyses above, we proposed a holistic model on the spatio-temporal expression of genes that likely govern sugar transport and utilization across maize maternal and endosperm and embryo tissues during the critical stage of grain set. Collectively, the findings represent an advancement towards a holistic understanding of the transcriptional landscape underlying post-phloem sugar transport in maize grain and indicate that the drought-induced changes in gene expression are attributable to low C status.

Microautophagy Mediates Vacuolar Delivery of Storage Proteins in Maize Aleurone Cells

The molecular machinery orchestrating microautophagy, whereby eukaryotic cells sequester autophagic cargo by direct invagination of the vacuolar/lysosomal membrane, is still largely unknown, especially in plants. Here, we demonstrate microautophagy of storage proteins in the maize aleurone cells of the endosperm and analyzed proteins with potential regulatory roles in this process. Within the cereal endosperm, starchy endosperm cells accumulate storage proteins (mostly prolamins) and starch whereas the peripheral aleurone cells store oils, storage proteins, and specialized metabolites. Although both cell types synthesize prolamins, they employ different pathways for their subcellular trafficking. Starchy endosperm cells accumulate prolamins in protein bodies within the endoplasmic reticulum (ER), whereas aleurone cells deliver prolamins to vacuoles via an autophagic mechanism, which we show is by direct association of ER prolamin bodies with the tonoplast followed by engulfment via microautophagy. To identify candidate proteins regulating this process, we performed RNA-seq transcriptomic comparisons of aleurone and starchy endosperm tissues during seed development and proteomic analysis on tonoplast-enriched fractions of aleurone cells. From these datasets, we identified 10 candidate proteins with potential roles in membrane modification and/or microautophagy, including phospholipase-Dα5 and a possible EUL-like lectin. We found that both proteins increased the frequency of tonoplast invaginations when overexpressed in Arabidopsis leaf protoplasts and are highly enriched at the tonoplast surface surrounding ER protein bodies in maize aleurone cells, thus supporting their potential connections to microautophagy. Collectively, this candidate list now provides useful tools to study microautophagy in plants.

Maternal-filial transfer structures in endosperm: A nexus of nutritional dynamics and seed development

Although the ultimate purpose of a seed is the successful establishment of the next generation, seed development involves more than embryo growth. In angiosperms, seed development requires the intimate coordination of three distinct entities - maternal tissue and two offspring, embryo and embryo-nourishing endosperm. Although seeds are cornerstones of many terrestrial ecosystems and human diets, we are only beginning to understand the interactions among seed tissues and the molecular processes and genes that determine them. Recent studies of gene expression and function in distantly related angiosperms, combined with over 100 years of embryological research, have repeatedly highlighted the endosperm associated with maternal-filial boundaries as a central point in seed developmental dynamics. In this review, we highlight evidence that links this zone with nutritional dynamics, developmental signaling, and imprinted gene expression. We suggest that the underappreciated diversity of this specialized endosperm across angiosperms deserves further study from developmental, molecular, and genetic perspectives.

Gene regulation network analyses of pistil development in papaya

BACKGROUND

The pistil is an essential part of flowers that functions in the differentiation of the sexes and reproduction in plants. The stigma on the pistil can accept pollen to allow fertilization and seed development. Papaya (Carica papaya L.) is a dioecious plant, where female flowers exhibit normal pistil, while the male flowers exhibit aborted pistil at a late stage of pistil development.

RESULTS

The developmental stages of papaya pistil were analyzed after first dividing it into slices representing the primordium stage 1 (S1), the pre-meiotic stages S2, post-meiotic stage S3, and the mitotic stage S4. The SS scoring algorithm analysis of genes preferentially expressed at different stages revealed differentially expressed genes between male and female flowers. A transcription factor regulatory network for each stage based on the genes that are differentially expressed between male and female flowers was constructed. Some transcription factors related to pistil development were revealed based on the analysis of regulatory networks such as CpAGL11, CpHEC2, and CpSUPL. Based on the specific expression of genes, constructed a gene regulatory subnetwork with CpAGL11-CpSUPL-CpHEC2 functioning as the core. Analysis of the functionally enriched terms in this network reveals several differentially expressed genes related to auxin/ brassinosteroid signal transduction in the plant hormone signal transduction pathway. At the same time, significant differences in the expression of auxin and brassinosteroid synthesis-related genes between male and female flowers at different developmental stages were detected.

CONCLUSIONS

The pistil abortion of papaya might be caused by the lack of expression or decreased expression of some transcription factors and hormone-related genes, affecting hormone signal transduction or hormone biosynthesis. Analysis of aborted and normally developing pistil in papaya provided new insights into the molecular mechanism of pistil development and sex differentiation in dioecious papaya.

Maize Endosperm Development: Tissues, Cells, Molecular Regulation and Grain Quality Improvement

Maize endosperm plays important roles in human diet, animal feed and industrial applications. Knowing the mechanisms that regulate maize endosperm development could facilitate the improvement of grain quality. This review provides a detailed account of maize endosperm development at the cellular and histological levels. It features the stages of early development as well as developmental patterns of the various individual tissues and cell types. It then covers molecular genetics, gene expression networks, and current understanding of key regulators as they affect the development of each tissue. The article then briefly considers key changes that have occurred in endosperm development during maize domestication. Finally, it considers prospects for how knowledge of the regulation of endosperm development could be utilized to enhance maize grain quality to improve agronomic performance, nutrition and economic value.

Integration of Transcriptome and Metabolome Provides New Insights to Flavonoids Biosynthesis in Dendrobium huoshanense

Dendrobium huoshanense is both a traditional herbal medicine and a plant of high ornamental and medicinal value. We used transcriptomics and metabolomics to investigate the effects of growth year on the secondary metabolites of D. huoshanense stems obtained from four different years of cultivation. In this study, a total of 428 differentially accumulated metabolites (DAMs) and 1802 differentially expressed genes (DEGs) were identified. The KEGG enrichment analysis of DEGs and DAMs revealed significant differences in "Flavonoid biosynthesis", "Phenylpropanoid biosynthesis" and "Flavone and flavonol biosynthesis". We summarize the biosynthesis pathway of flavonoids in D. huoshanense, providing new insights into the biosynthesis and regulation mechanisms of flavonoids in D. huoshanense. Additionally, we identified two candidate genes, FLS (LOC110107557) and F3'H (LOC110095936), which are highly involved in flavonoid biosynthesis pathway, by WGCNA analysis. The aim of this study is to investigate the effects of growth year on secondarily metabolites in the plant and provide a theoretical basis for determining a reasonable harvesting period for D. huoshanense.

ZmCTLP1 is required for the maintenance of lipid homeostasis and the basal endosperm transfer layer in maize kernels

Maize kernel weight is influenced by the unloading of nutrients from the maternal placenta and their passage through the transfer tissue of the basal endosperm transfer layer (BETL) and the basal intermediate zone (BIZ) to the upper part of the endosperm. Here, we show that Small kernel 10 (Smk10) encodes a choline transporter-like protein 1 (ZmCTLP1) that facilitates choline uptake and is located in the trans-Golgi network (TGN). Its loss of function results in reduced choline content, leading to smaller kernels with a lower starch content. Mutation of ZmCTLP1 disrupts membrane lipid homeostasis and the normal development of wall in-growths. Expression levels of Mn1 and ZmSWEET4c, two kernel filling-related genes, are downregulated in the smk10, which is likely to be one of the major causes of incompletely differentiated transfer cells. Mutation of ZmCTLP1 also reduces the number of plasmodesmata (PD) in transfer cells, indicating that the smk10 mutant is impaired in PD formation. Intriguingly, we also observed premature cell death in the BETL and BIZ of the smk10 mutant. Together, our results suggest that ZmCTLP1-mediated choline transport affects kernel development, highlighting its important role in lipid homeostasis, wall in-growth formation and PD development in transfer cells.

ZmGRAS11, transactivated by Opaque2, positively regulates kernel size in maize

Although the genetic basis for endosperm development in maize (Zea mays) has been well studied, the mechanism for coordinating grain filling with increasing kernel size remains elusive. Here, we report that increased kernel size was selected during modern breeding and identify a novel DELLA-like transcriptional regulator, ZmGRAS11, which positively regulates kernel size and kernel weight in maize. We find that Opaque2, a core transcription factor for zein protein and starch accumulation, transactivates the expression of ZmGRAS11. Our data suggest that the Opaque2-ZmGRAS11 module mediates synergistic endosperm enlargement with grain filling.

A systems genetics approach reveals PbrNSC as a regulator of lignin and cellulose biosynthesis in stone cells of pear fruit

BACKGROUND

Stone cells in fruits of pear (Pyrus pyrifolia) negatively influence fruit quality because their lignified cell walls impart a coarse and granular texture to the fruit flesh.

RESULTS

We generate RNA-seq data from the developing fruits of 206 pear cultivars with a wide range of stone cell contents and use a systems genetics approach to integrate co-expression networks and expression quantitative trait loci (eQTLs) to characterize the regulatory mechanisms controlling lignocellulose formation in the stone cells of pear fruits. Our data with a total of 35,897 expressed genes and 974,404 SNPs support the identification of seven stone cell formation modules and the detection of 139,515 eQTLs for 3229 genes in these modules. Focusing on regulatory factors and using a co-expression network comprising 39 structural genes, we identify PbrNSC as a candidate regulator of stone cell formation. We then verify the function of PbrNSC in regulating lignocellulose formation using both pear fruit and Arabidopsis plants and further show that PbrNSC can transcriptionally activate multiple target genes involved in secondary cell wall formation.

CONCLUSIONS

This study generates a large resource for studying stone cell formation and provides insights into gene regulatory networks controlling the formation of stone cell and lignocellulose.

Comparative transcriptomics and network analysis define gene coexpression modules that control maize aleurone development and auxin signaling

The naked endosperm1 (nkd1), naked endosperm2 (nkd2), and thick aleurone1 (thk1) genes are important regulators of maize (Zea mays L.) endosperm development. Double mutants of nkd1 and nkd2 (nkd1,2) show multiple aleurone (AL) cell layers with disrupted AL cell differentiation, whereas mutants of thk1 cause multiple cell layers of fully differentiated AL cells. Here, we conducted a comparative analysis of nkd1,2 and thk1 mutant endosperm transcriptomes to study how these factors regulate gene networks to control AL layer specification and cell differentiation. Weighted gene coexpression network analysis was incorporated with published laser capture microdissected transcriptome datasets to identify a coexpression module associated with AL development. In this module, both Nkd1,2+ and Thk1+ appear to regulate cell cycle and division, whereas Nkd1,2+, but not Thk1+, regulate auxin signaling. Further investigation of nkd1,2 differentially expressed genes combined with published putative targets of auxin response factors (ARFs) identified 61 AL-preferential genes that may be directly activated by NKD-modulated ARFs. All 61 genes were upregulated in nkd1,2 mutant and the enriched Gene Ontology terms suggested that they are associated with hormone crosstalk, lipid metabolism, and developmental growth. Expression of a transgenic DR5-red fluorescent protein auxin reporter was significantly higher in nkd1,2 mutant endosperm than in wild type, supporting the prediction that Nkd1,2+ negatively regulate auxin signaling in developing AL. Overall, these results suggest that Nkd1,2+ and Thk1+ may normally restrict AL development to a single cell layer by limiting cell division, and that Nkd1,2+ restrict auxin signaling in the AL to maintain normal cell patterning and differentiation processes.

Metabolome and transcriptome profiling of Theobroma cacao provides insights into the molecular basis of pod color variation

The Theobroma cacao presents a wide diversity in pod color among different cultivars. Although flavonoid biosynthesis has been studied in many plants, molecular mechanisms governing the diversity of coloration in cacao pods are largely unknown. The flavonoid metabolite profiles and flavonoid biosynthetic gene expression in the pod exocarps of light green pod 'TAS 410' (GW), green pod 'TAS 166' (GF), and mauve pod 'TAS 168' (PF) were determined. Changes in flavonoid metabolites, particularly the anthocyanins (cyanidin 3-O-galactoside, cyanidin 3-O-glucoside, and cyanidin O-syringic acid) were significantly up-accumulated in the mauve phenotype (PF) compared to the light green or green phenotypes, endowing the pod color change from light green or green to mauve. Consistently, the PF phenotype showed different expression patterns of flavonoid biosynthetic structural genes in comparison with GW/GF phenotypes. The expression level of LAR and ANR in GW/GF was significantly higher than PF, while the expression level of UFGT in GW/GF was lower than PF. These genes likely generated more anthocyanins in the exocarps samples of PF than that of GW/GF. Simultaneously, colorless flavan-3-ols (catechin, epicatechin and proanthocyanidin) content in the exocarp samples of PF was lower than GW/GF. Additionally, MYB (gene18079) and bHLH (gene5045 and gene21575) may participate in the regulation of the pod color. This study sheds light on the molecular basis of cacao pod color variation, which will contribute to breeding cacao varieties with enhanced flavonoid profiles for nutritional applications.

Dynamic network biomarker analysis discovers IbNAC083 in the initiation and regulation of sweet potato root tuberization

The initiation and development of storage roots (SRs) are intricately regulated by a transcriptional regulatory network. One key challenge is to accurately pinpoint the tipping point during the transition from pre-swelling to SRs and to identify the core regulators governing such a critical transition. To solve this problem, we performed a dynamic network biomarker (DNB) analysis of transcriptomic dynamics during root development in Ipomoea batatas (sweet potato). First, our analysis identified stage-specific expression patterns for a significant proportion (>9%) of the sweet potato genes and unraveled the chronology of events that happen at the early and later stages of root development. Then, the results showed that different root developmental stages can be depicted by co-expressed modules of sweet potato genes. Moreover, we identified the key components and transcriptional regulatory network that determine root development. Furthermore, through DNB analysis an early stage, with a root diameter of 3.5 mm, was identified as the critical period of SR swelling initiation, which is consistent with morphological and metabolic changes. In particular, we identified a NAM/ATAF/CUC (NAC) domain transcription factor, IbNAC083, as a core regulator of this initiation in the DNB-associated network. Further analyses and experiments showed that IbNAC083, along with its associated differentially expressed genes, induced dysfunction of metabolism processes, including the biosynthesis of lignin, flavonol and starch, thus leading to the transition to swelling roots.

Genome-wide identification, phylogeny and expression analysis of AP2/ERF transcription factors family in sweet potato

Background

In recent years, much attention has been given to AP2/ERF transcription factors because they play indispensable roles in many biological processes, such as plant development and biotic and abiotic stress responses. Although AP2/ERFs have been thoroughly characterised in many plant species, the knowledge about this family in the sweet potato, which is a vital edible and medicinal crop, is still limited. In this study, a comprehensive genome-wide investigation was conducted to characterise the AP2/ERF gene family in the sweet potato.

Results

Here, 198 IbAP2/ERF transcription factors were obtained. Phylogenetic analysis classified the members of the IbAP2/ERF family into three groups, namely, ERF (172 members), AP2 (21 members) and RAV (5 members), which was consistent with the analysis of gene structure and conserved protein domains. The evolutionary characteristics of these IbAP2/ERF genes were systematically investigated by analysing chromosome location, conserved protein motifs and gene duplication events, indicating that the expansion of the IbAP2/ERF gene family may have been caused by tandem duplication. Furthermore, the analysis of cis-acting elements in IbAP2/ERF gene promoters implied that these genes may play crucial roles in plant growth, development and stress responses. Additionally, the available RNA-seq data and quantitative real-time PCR (qRT-PCR) were used to investigate the expression patterns of IbAP2/ERF genes during sweet potato root development as well as under multiple forms of abiotic stress, and we identified several developmental stage-specific and stress-responsive IbAP2/ERF genes. Furthermore, g59127 was differentially expressed under various stress conditions and was identified as a nuclear protein, which was in line with predicted subcellular localization results.

Conclusions

This study originally revealed the characteristics of the IbAP2/ERF superfamily and provides valuable resources for further evolutionary and functional investigations of IbAP2/ERF genes in the sweet potato.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-08043-w.

Rejuvenation remodels transcriptional network to improve rhizogenesis in mature Juglans tree

Adventitious rooting of walnut species (Juglans L.) is known to be rather difficult, especially for mature trees. The adventitious root formation (ARF) capacities of mature trees can be significantly improved by rejuvenation. However, the underlying gene regulatory networks (GRNs) of rejuvenation remain largely unknown. To characterize such regulatory networks, we carried out the transcriptomic study using RNA samples of the cambia and peripheral tissues on the bottom of rejuvenated and mature walnut (Juglans hindsii × J. regia) cuttings during the ARF. The RNA sequencing data suggested that zeatin biosynthesis, energy metabolism and substance metabolism were activated by rejuvenation, whereas photosynthesis, fatty acid biosynthesis and the synthesis pathways for secondary metabolites were inhibited. The inter- and intra-module GRNs were constructed using differentially expressed genes. We identified 35 hub genes involved in five modules associated with ARF. Among these hub genes, particularly, beta-glucosidase-like (BGLs) family members involved in auxin metabolism were overexpressed at the early stage of the ARF. Furthermore, BGL12 from the cuttings of Juglans was overexpressed in Populus alba × P. glandulosa. Accelerated ARF and increased number of ARs were observed in the transgenic poplars. These results provide a high-resolution atlas of gene activity during ARF and help to uncover the regulatory modules associated with the ARF promoted by rejuvenation.

Weighted gene co-expression network analysis unveils gene networks regulating folate biosynthesis in maize endosperm

Folates are essential elements for human growth and development, and their deficiency can lead to serious disorders. Waxy maize is a rich source of folates; however, the regulatory mechanism underlying folate biosynthesis in the endosperm remains unclear. Here, we examined changes in the folate content of maize endosperm collected at 15, 18, 21, 24, and 27 days after pollination (DAP) using liquid chromatograph-mass spectrometry and identified genes related to folate biosynthesis using transcriptome sequencing data. The results showed that 5-methyl-tetrahydrofolate and 5,10-methylene tetrahydrofolate were the main storage forms of folates in the endosperm, and their contents were relatively high at 21-24 days. We also identified 569, 3183, 4365, and 5513 differentially expressed genes (DEGs) in different days around milk stage. Functional annotation revealed 518 transcription factors (TFs) belonging to 33 families exhibiting specific expression in at least one sampling time. The key hub genes involved in folate biosynthesis were identified by weighted gene co-expression network analysis. In total, 24,976 genes were used to construct a co-expression network with 29 co-expression modules, among which the brown and purple modules were highly related to folate biosynthesis. Further, 187 transcription factors in the brown and purple modules were considered potential transcription factors related to endosperm folate biosynthesis. These results may improve the understanding of the molecular mechanism underlying folate biosynthesis in waxy maize and lead to the development of nutritionally fortified varieties.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s13205-021-02974-7.

High-resolution spatiotemporal transcriptome analyses during cellularization of rice endosperm unveil the earliest gene regulation critical for aleurone and starchy endosperm cell fate specification

The major tissues of the cereal endosperm are the starchy endosperm (SE) in the inner and the aleurone layer (AL) at the outer periphery. The fates of the cells that comprise these tissues are determined according to positional information; however, our understanding of the underlying molecular mechanisms remains limited. Here, we conducted a high-resolution spatiotemporal analysis of the rice endosperm transcriptome during early cellularization. In rice, endosperm cellularization proceeds in a concentric pattern from a primary alveolus cell layer, such that developmental progression can be defined by the number of cell layers. Using laser-capture microdissection to obtain precise tissue sections, transcriptomic changes were followed through five histologically defined stages of cellularization from the syncytial to 3-cell layer (3 L) stage. In addition, transcriptomes were compared between the inner and the outermost peripheral cell layers. Large differences in the transcriptomes between stages and between the inner and the peripheral cells were found. SE attributes were expressed at the alveolus-cell-layer stage but were preferentially activated in the inner cell layers that resulted from periclinal division of the alveolus cell layer. Similarly, AL attributes started to be expressed only after the 2 L stage and were localized to the outermost peripheral cell layer. These results indicate that the first periclinal division of the alveolus cell layer is asymmetric at the transcriptome level, and that the cell-fate-specifying positional cues and their perception system are already operating before the first periclinal division. Several genes related to epidermal identity (i.e., type IV homeodomain-leucine zipper genes and wax biosynthetic genes) were also found to be expressed at the syncytial stage, but their expression was localized to the outermost peripheral cell layer from the 2 L stage onward. We believe that our findings significantly enhance our knowledge of the mechanisms underlying cell fate specification in rice endosperm.

Maize unstable factor for orange1 is essential for endosperm development and carbohydrate accumulation

Maize (Zea mays L.) Ufo1-1 is a spontaneous dominant mutation of the unstable factor for orange1 (ufo1). We recently cloned ufo1, which is a Poaceae-specific gene highly expressed during seed development in maize. Here, we have characterized Ufo1-1 and a loss-of-function Ds insertion allele (ufo1-Dsg) to decipher the role of ufo1 in maize. We found that both ufo1 mutant alleles impact sugars and hormones, and have defects in the basal endosperm transfer layer (BETL) and adjacent cell types. The Ufo1-1 BETL had reduced cell elongation and cell wall ingrowth, resulting in cuboidal shaped transfer cells. In contrast, the ufo1-Dsg BETL cells showed a reduced overall size with abnormal wall ingrowth. Expression analysis identified the impact of ufo1 on several genes essential for BETL development. The overexpression of Ufo1-1 in various tissues leads to ectopic phenotypes, including abnormal cell organization and stomata subsidiary cell defects. Interestingly, pericarp and leaf transcriptomes also showed that as compared with wild type, Ufo1-1 had ectopic expression of endosperm development-specific genes. This study shows that Ufo1-1 impacts the expression patterns of a wide range of genes involved in various developmental processes.

Localised expression of OsIAA29 suggests a key role for auxin in regulating development of the dorsal aleurone of early rice grains

Non-canonical AUX/IAA protein, OsIAA29, and ZmMPR-1 homologues, OsMRPLs, are part of an auxin-related signalling cascade operating in the dorsal aleurone during early rice grain development. Endosperm of rice and other cereals accumulates high concentrations of the predominant in planta auxin, indole-3-acetic acid (IAA) during early grain development. However, IAA signalling and function during endosperm development are poorly understood. Here, we report that OsYUC12 (an auxin biosynthesis gene) and OsIAA29 (encoding a non-canonical AUX/IAA) are both expressed exclusively in grains, reaching a maximum 5-6 days after pollination. OsYUC12 expression is localised in the aleurone, sub-aleurone and embryo, whereas OsIAA29 expression is restricted to a narrow strip in the dorsal aleurone, directly under the vascular bundle. Although rice has been reported to lack endosperm transfer cells (ETCs), this region of the aleurone is enriched with sugar transporters and is likely to play a key role in apoplastic nutrient transfer, analogous to ETCs in other cereals. OsIAA29 has orthologues only in grass species; expression of which is also specific to early grain development. OsYUC12 and OsIAA29 are temporally co-expressed with two genes (AL1 and OsPR602) previously linked to the development of dorsal aleurone or ETCs. Also up-regulated at the same time is a cluster of MYB-related genes (designated OsMRPLs) homologous to ZmMRP-1, which regulates maize ETC development. Wheat homologues of ZmMRP-1 are similarly expressed in ETCs. Although previous work has suggested that other cereals do not have orthologues of ZmMRP-1, our work suggests OsIAA29 and OsMRPLs and their homologues in other grasses are part of an auxin-regulated, conserved signalling network involved in the differentiation of cells with ETC-like function in developing cereal grains.

Comparative transcriptome analysis reveals regulatory network and regulators associated with proanthocyanidin accumulation in persimmon

BACKGROUND

Proanthocyanidins (PAs) are important plant secondary metabolites that confer flavor, nutritional value, and resistance to pathogens. Persimmon is one of the PA richest crops. Mature fruits can be inedible because of the astringency caused by high PA levels and need to go through a de-astringency treatment before consumption. The molecular basis for PA accumulation is poorly known, particularly transcriptional regulators. We characterised three genotypes ('Luotiantianshi' (LT), 'Mopanshi' (MP), and 'Youhou' (YH)) with different PA accumulation patterns using an approach that combined PacBio full-length sequencing and Illumina-based RNA sequencing to build high-quality full-length transcriptomes. Additionally, we analysed transcriptome dynamics of the three genotypes (LT, MP, and YH) at four key fruit developmental stages.

RESULTS

A total of 96,463 transcripts were obtained. We identified 80,075 protein-coding sequences (CDSs), 71,137 simple sequence repeats (SSRs), and 27,845 long noncoding RNAs (lncRNAs). Pearson correlation coefficient (PCC), principal component analysis (PCA), and differentially expressed transcripts (DETs) analyses indicated that the four different developmental stages within a genotype exhibited similar transcriptome activities. A total of 2,164 transcripts specific to each fruit developmental stage were detected. The transcripts specific to early stages were attributed to phenylpropanoid and flavonoid biosynthesis. Co-expression network analyses revealed MEbrown and MEblue modules were strongly associated to PA accumulation. From these two modules, 20 hub TFs are potential regulators for PA accumulation. Among them, Cluster_78388 (SBP protein), Cluster_63454 (bZIP protein), and Cluster_66595 (MYB protein) appear to involve in the PA biosynthesis in Chinese genotypes.

CONCLUSIONS

This is the first high-quality reference transcriptome for commercial persimmon. Our work provides insights into the molecular pathways underlying PA accumulation and enhances our global understanding of transcriptome dynamics throughout fruit development.

Transcriptome analysis of differentially expressed genes in chrysanthemum MET1 RNA interference lines

DNA methylation is the most important epigenetic modification involved in many essential biological processes. MET1 is one of DNA methyltransferases that affect the level of methylation in the entire genome. To explore the effect of MET1 gene silencing on gene expression profile of Chrysanthemum × morifolium 'Zijingling'. The stem section and leaves at the young stage were taken for transcriptome sequencing. MET1-_RNAi leaves had 8 differentially expressed genes while 156 differentially expressed genes were observed in MET1-_RNAi stem compared with control leaves and stem. These genes encode many key proteins in plant biological processes, such as transcription factors, signal transduction mechanisms, secondary metabolite synthesis, transport and catabolism and interaction. In general, 34.58% of the differentially expressed genes in leaves and stems were affected by the reduction of the MET1 gene. The differentially expressed genes in stem and leaves of transgenic plants went through significant changes. We found adequate amount of candidate genes associated with flowering, however, the number of genes with significant differences between transgenic and control lines was not too high. Several flowering related genes were screened out for gene expression verification and all of them were obseved as consistent with transcriptome data. These candidate genes may play important role in flowering variation of chrysanthemum. This study reveals the mechanism of CmMET1 interference on the growth and development of chrysanthemum at the transcriptional level, which provides the basis for further research on the epigenetic regulation mechanism in flower induction and development.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s12298-021-01022-1.

High-spatiotemporal-resolution transcriptomes provide insights into fruit development and ripening in Citrus sinensis

Citrus fruit has a unique structure with soft leathery peel and pulp containing vascular bundles and several segments with many juice sacs. The function and morphology of each fruit tissue are different. Therefore, analysis at the organ-wide or mixed-tissue level inevitably obscures many tissue-specific phenomena. High-throughput RNA sequencing was used to profile Citrus sinensis fruit development based on four fruit tissue types and six development stages from young fruits to ripe fruits. Using a coexpression network analysis, modules of coexpressed genes and hub genes of tissue-specific networks were identified. Of particular, importance is the discovery of the regulatory network of phytohormones during citrus fruit development and ripening. A model was proposed to illustrate how ABF2 mediates the ABA signalling involved in sucrose transport, chlorophyll degradation, auxin homoeostasis, carotenoid and ABA biosynthesis, and cell wall metabolism during citrus fruit development. Moreover, we depicted the detailed spatiotemporal expression patterns of the genes involved in sucrose and citric acid metabolism in citrus fruit and identified several key genes that may play crucial roles in sucrose and citric acid accumulation in the juice sac, such as SWEET15 and CsPH8. The high spatial and temporal resolution of our data provides important insights into the molecular networks underlying citrus fruit development and ripening.

Transcriptomic and Weighted Gene Co-expression Correlation Network Analysis Reveal Resveratrol Biosynthesis Mechanisms Caused by Bud Sport in Grape Berry

Resveratrol is a natural polyphenol compound produced in response to biotic and abiotic stresses in grape berries. However, changes in resveratrol caused by bud sport in grapes are scarcely reported. In this study, trans-resveratrol and cis-resveratrol were identified and quantified in the grape berries of 'Summer Black' and its bud sport 'Nantaihutezao' from the veraison to ripening stages using ultra performance liquid chromatography-high resolution tandem mass spectrometry (UPLC-HRMS). We found that bud sport accumulates the trans-resveratrol earlier and increases the contents of cis-resveratrol in the earlier stages but decreases its contents in the later stages. Simultaneously, we used RNA-Seq to identify 51 transcripts involved in the stilbene pathways. In particular, we further identified 124 and 19 transcripts that negatively correlated with the contents of trans-resveratrol and cis-resveratrol, respectively, and four transcripts encoding F3'5'H that positively correlated with the contents of trans-resveratrol by weighted gene co-expression network analysis (WGCNA). These transcripts may play important roles in relation to the synergistic regulation of metabolisms of resveratrol. The results of this study can provide a theoretical basis for the genetic improvement of grapes.

Transcriptome analysis reveals the mechanism associated with dynamic changes in fatty acid and phytosterol content in foxtail millet (Setaria italica) during seed development

Foxtail millet (Setaria italica) is an excellent source of beneficial natural fatty acids and phytosterols. However, the mechanisms underlying the dynamic changes of fatty acids and phytosterols during seed development are unknown. In this study, a comprehensive dynamic change analysis of the bioactive compounds during seed development was conducted in two cultivars with different crude fat content (high-fat, JG 35 [5.40%]; and low-fat, JG 39 [2.90%]). GC-FID/MS analysis showed that the proportion of unsaturated fatty acids (UFAs) were higher than the saturated fatty acids (SFAs). UFA content first increased, then decreased during seed development, while SFA content showed the opposite trend. Oil contents continuously increased with seed development, especially at the S2 stage. Phytosterol contents initially increased, then decreased with seed development. Transcriptome analysis revealed that 152 genes were associated with fatty acid metabolism and phytosterol biosynthesis, of which 46 and 62 were related to UFA and phytosterol biosynthesis, respectively. Furthermore, the key genes involved in fatty acid synthesis (ACCase and FATA/B), triacylglycerol biosynthesis (LACS, GPAT, and DGAT), and phytosterols synthesis (CAS1, STM1, EGR6, and DWF1) were overexpressed. This led to maximum UFA, oil, and phytosterol accumulation in JG 35 at the S2 stage. This study reveals the mechanism behind the dynamic changes of fatty acid and phytosterol contents in foxtail millet during seed development.

Whitefly Network Analysis Reveals Gene Modules Involved in Host Plant Selection, Development and Evolution

Whiteflies are Hemipterans that typically feed on the undersides of plant leaves. They cause severe damage by direct feeding as well as transmitting plant viruses to a wide range of plants. However, it remains largely unknown which genes play a key role in development and host selection. In this study, weighted gene co-expression network analysis was applied to construct gene co-expression networks in whitefly. Nineteen gene co-expression modules were detected from 15560 expressed genes of whitefly. Combined with the transcriptome data of salivary glands and midgut, we identified three gene co-expression modules related to host plant selection. These three modules contain genes related to host-plant recognition, such as detoxification genes, chemosensory genes and some salivary gland-associated genes. Results of Gene Ontology and Kyoto Encyclopedia of Genes and Genomes analyses elucidated the following pathways involved in these modules: lysosome, metabolic and detoxification pathways. The modules related to the development contain two co-expression modules; moreover, the genes were annotated to the development of chitin-based cuticle. This analysis provides a basis for future functional analysis of genes involved in host-plant recognition.

Miniature Seed6, encoding an endoplasmic reticulum signal peptidase, is critical in seed development

Endoplasmic reticulum (ER) type I signal peptidases (ER SPases I) are vital proteases that cleave signal peptides from secreted proteins. However, the specific function of ER SPase I in plants has not been genetically characterized, and the substrate is largely unknown. Here, we report the identification of a maize (Zea mays) miniature seed6 (mn6) mutant. The loss-of-function mn6 mutant exhibited severely reduced endosperm size. Map-based cloning and molecular characterization indicated that Mn6 is an S26-family ER SPase I, with Gly102 (box E) in Mn6 critical for protein function during processing. Mass spectrometric and immunoprecipitation analyses revealed that Mn6 is predominantly involved in processing carbohydrate synthesis-related proteins, including the cell wall invertase miniature seed1 (Mn1), which is specifically expressed in the basal endosperm transfer layer. RNA and protein expression levels of Mn1 were both significantly downregulated in the mn6 mutant. Due to the significant reduction in cell wall invertase activity in the transfer cell layer, mutation of Mn6 caused dramatic defects in endosperm development. These results suggest that proper maturation of Mn1 by Mn6 may be a crucial step for proper seed filling and maize development.

Maize endosperm development

Recent breakthroughs in transcriptome analysis and gene characterization have provided valuable resources and information about the maize endosperm developmental program. The high temporal-resolution transcriptome analysis has yielded unprecedented access to information about the genetic control of seed development. Detailed spatial transcriptome analysis using laser-capture microdissection has revealed the expression patterns of specific populations of genes in the four major endosperm compartments: the basal endosperm transfer layer (BETL), aleurone layer (AL), starchy endosperm (SE), and embryo-surrounding region (ESR). Although the overall picture of the transcriptional regulatory network of endosperm development remains fragmentary, there have been some exciting advances, such as the identification of OPAQUE11 (O11) as a central hub of the maize endosperm regulatory network connecting endosperm development, nutrient metabolism, and stress responses, and the discovery that the endosperm adjacent to scutellum (EAS) serves as a dynamic interface for endosperm-embryo crosstalk. In addition, several genes that function in BETL development, AL differentiation, and the endosperm cell cycle have been identified, such as ZmSWEET4c, Thk1, and Dek15, respectively. Here, we focus on current advances in understanding the molecular factors involved in BETL, AL, SE, ESR, and EAS development, including the specific transcriptional regulatory networks that function in each compartment during endosperm development.

Conserved co-functional network between maize and Arabidopsis aid in the identification of seed defective genes in maize

BACKGROUND

The biological pathways related to Arabidopsis seed development have been well studied and functional genes involved in it have been discovered. However, functional studies about maize seed development were more limited compared to Arabidopsis.

OBJECTIVE

Therefore, transferring knowledge from Arabidopsis into maize would facilitate functional studies about maize seed development.

METHOD

In this study, public transcriptome data of the two species related to seed development were obtained. Co-expression network in each species was compared by integrating orthology information.

RESULTS

This conserved co-functional network contained 4510 maize and 4808 Arabidopsis genes, respectively. Most of these genes were expressed in throughout embryo, early or later endosperm/seed. These conserved co-functional genes were significantly enriched for members of PPR protein family, which was consistent with that PPR proteins play an important role in maize seed development. Spatial-temporally co-functional genes were discovered in the seed coat and embryo. Furthermore, 66 well-studied genes involved in Arabidopsis seed development were co-functional with 319 maize genes and one maize gene (GRMZM2G036050) was further confirmed using an EMS-induced seed defective mutant by bulked segregating RNA sequencing (BSR) analysis.

CONCLUSIONS

Altogether, these results showed the potential of this approach to support functional studies in maize seed development by transferring knowledge from Arabidopsis.

Transcriptome analysis reveals significant difference in gene expression and pathways between two peanut cultivars under Al stress

Aluminum (Al) toxicity is an important factor in limiting peanut growth on acidic soil. The molecular mechanisms underlying peanut responses to Al stress are largely unknown. In this study, we performed transcriptome analysis of the root tips (0-1 cm) of peanut cultivar ZH2 (Al-sensitive) and 99-1507 (Al-tolerant) respectively. Root tips of peanuts that treated with 100 μM Al for 8 h and 24 h were analyzed by RNA-Seq, and a total of 8,587 differentially expressed genes (DEGs) were identified. GO and KEGG pathway analysis excavated a group of important Al-responsive genes related to organic acid transport, metal cation transport, transcription regulation and programmed cell death (PCD). These homologs were promising targets to modulate Al tolerance in peanuts. It was found that the rapid transcriptomic response to Al stress in 99-1507 helped to activate effective Al tolerance mechanisms. Protein and protein interaction analysis indicated that MAPK signal transduction played important roles in the early response to Al stress in peanuts. Moreover, weighted correlation network analysis (WGCNA) identified a predicted EIL (EIN3-like) gene with greatly increased expression as an Al-associated gene, and revealed a link between ethylene signaling transduction and Al resistance related genes in peanut, which suggested the enhanced signal transduction mediated the rapid transcriptomic responses. Our results revealed key pathways and genes associated with Al stress, and improved the understanding of Al response in peanut.

ZmSKS13, a cupredoxin domain-containing protein, is required for maize kernel development via modulation of redox homeostasis

The SKU5 similar (SKS) genes encode a family of multi-copper-oxidase-like proteins with cupredoxin domains similar to those in laccase and ascorbate oxidase. Although SKS proteins are known to function in root growth and cotyledon vascular patterning in Arabidopsis, their role in plant reproductive processes is poorly understood. Here, we identified a seed mutant of maize (Zea mays), generated by ethyl methane sulfonate (EMS) mutagenesis, that we designated defective kernel-zk1 (dek-zk1). The mutant produced small, shriveled kernels with an aberrant basal endosperm transfer layer (BETL) and placento-chalazal (PC) layer and irregular starch granules. Map-based cloning revealed that Dek-zk1 encodes an SKU5 similar 13 (GenBank: ONM36900.1), so it was named ZmSKS13. ZmSKS13 comprises a paralogous pair with Zm00001d012524, but the transcript abundance of ZmSKS13 in developing kernels is 15 times higher than that of Zm00001d012524, resulting in dek-zk1 mutation conveying a distinct kernel phenotype. ZmSKS13 loss of function led to overaccumulation of reactive oxygen species (ROS) and severe DNA damage in the nucellus and BETL and PC layer cells, and exogenous antioxidants significantly alleviated the defects of the mutant kernels. Our results thus demonstrate that ZmSKS13 is a novel regulator that plays a crucial role in kernel development in maize through the modulation of ROS homeostasis.

Identification of Key Tissue-Specific, Biological Processes by Integrating Enhancer Information in Maize Gene Regulatory Networks

Enhancers are key players in the spatio-temporal coordination of gene expression during numerous crucial processes, including tissue differentiation across development. Characterizing the transcription factors (TFs) and genes they connect, and the molecular functions underpinned is important to better characterize developmental processes. In plants, the recent molecular characterization of enhancers revealed their capacity to activate the expression of several target genes. Nevertheless, identifying these target genes at a genome-wide level is challenging, particularly for large-genome species, where enhancers and target genes can be hundreds of kilobases away. Therefore, the contribution of enhancers to plant regulatory networks remains poorly understood. Here, we investigate the enhancer-driven regulatory network of two maize tissues at different stages: leaves at seedling stage (V2-IST) and husks (bracts) at flowering. Using systems biology, we integrate genomic, epigenomic, and transcriptomic data to model the regulatory relationships between TFs and their potential target genes, and identify regulatory modules specific to husk and V2-IST. We show that leaves at the V2-IST stage are characterized by the response to hormones and macromolecules biogenesis and assembly, which are regulated by the BBR/BPC and AP2/ERF TF families, respectively. In contrast, husks are characterized by cell wall modification and response to abiotic stresses, which are, respectively, orchestrated by the C2C2/DOF and AP2/EREB families. Analysis of the corresponding enhancer sequences reveals that two different transposable element families (TIR transposon Mutator and MITE Pif/Harbinger) have shaped part of the regulatory network in each tissue, and that MITEs have provided potential new TF binding sites involved in husk tissue-specificity.

Maize kernel development

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

Cereal asparagine synthetase genes

Asparagine synthetase catalyses the transfer of an amino group from glutamine to aspartate to form glutamate and asparagine. The accumulation of free (nonprotein) asparagine in crops has implications for food safety because free asparagine is the precursor for acrylamide, a carcinogenic contaminant that forms during high-temperature cooking and processing. Here we review publicly available genome data for asparagine synthetase genes from species of the Pooideae subfamily, including bread wheat and related wheat species (Triticum and Aegilops spp.), barley (Hordeum vulgare) and rye (Secale cereale) of the Triticeae tribe. Also from the Pooideae subfamily: brachypodium (Brachypodium dIstachyon) of the Brachypodiae tribe. More diverse species are also included, comprising sorghum (Sorghum bicolor) and maize (Zea mays) of the Panicoideae subfamily and rice (Oryza sativa) of the Ehrhartoideae subfamily. The asparagine synthetase gene families of the Triticeae species each comprise five genes per genome, with the genes assigned to four groups: 1, 2, 3 (subdivided into 3.1 and 3.2) and 4. Each species has a single gene per genome in each group, except that some bread wheat varieties (genomes AABBDD) and emmer wheat (Triticum dicoccoides; genomes AABB) lack a group 2 gene in the B genome. This raises questions about the ancestry of cultivated pasta wheat and the B genome donor of bread wheat, suggesting that the hybridisation event that gave rise to hexaploid bread wheat occurred more than once. In phylogenetic analyses, genes from the other species cluster with the Triticeae genes, but brachypodium, sorghum and maize lack a group 2 gene, while rice has only two genes, one group 3 and one group 4. This means that TaASN2, the most highly expressed asparagine synthetase gene in wheat grain, has no equivalent in maize, rice, sorghum or brachypodium. An evolutionary pathway is proposed in which a series of gene duplications gave rise to the five genes found in modern Triticeae species.

Phenotypic Characteristics and Transcriptome of Cucumber Male Flower Development Under Heat Stress

Cucumber (Cucumis sativus L.) is an important vegetable crop, which is thermophilic not heat resistant. High-temperature stress always results in sterility at reproductive stage. In the present study, we evaluate the male flower developmental changes under normal (CK) and heat stress (HS) condition. After HS, the activities of peroxidase (POD) and superoxide dismutase (SOD) and the contents of malondialdehyde (MDA) were increased. In addition, the pollen fertility was significantly decreased; and abnormal tapetum and microspore were observed by paraffin section. Transcriptome analysis results presented that total of 5828 differentially expressed genes (DEGs) were identified after HS. Among these DEGs, 20 DEGs were found at four stages, including DNA binding transcription factor, glycosyltransferase, and wound-responsive family protein. The gene ontology term of carbohydrate metabolic process was significantly enriched in all anther stages, and many saccharides and starch synthase-related genes, such as invertase, sucrose synthase, and starch branching enzyme, were significantly different expressed in HS compared with CK. Furthermore, co-expression network analysis showed a module (midnightblue) strongly consistent with HS, and two hub genes (CsaV3_6G004180 and CsaV3_5G034860) were found with a high degree of connectivity to other genes. Our results provide comprehensive understandings on male flower development in cucumber under HS.

Single-Cell Genomics and Epigenomics: Technologies and Applications in Plants

The development of genomics and epigenomics has allowed rapid advances in our understanding of plant biology. However, conventional bulk analysis dilutes cell-specific information by providing only average information, thereby limiting the resolution of genomic and functional genomic studies. Recent advances in single-cell sequencing technology concerning genomics and epigenomics open new avenues to dissect cell heterogeneity in multiple biological processes. Recent applications of these approaches to plants have provided exciting insights into diverse biological questions. We highlight the methodologies underlying the current techniques of single-cell genomics and epigenomics before covering their recent applications, potential significance, and future perspectives in plant biology.

The thick aleurone1 Gene Encodes a NOT1 Subunit of the CCR4-NOT Complex and Regulates Cell Patterning in Endosperm

Maize (Zea mays) thick aleurone1 (thk1-R) mutants form multiple aleurone layers in the endosperm and have arrested embryogenesis. Prior studies suggest that thk1 functions downstream of defective kernel1 (dek1) in a regulatory pathway that controls aleurone cell fate and other endosperm traits. The original thk1-R mutant contained an ∼2-Mb multigene deletion, which precluded identification of the causal gene. Here, ethyl methanesulfonate mutagenesis produced additional alleles, and RNA sequencing from developing endosperm was used to identify a candidate gene based on differential expression compared with the wild-type progenitor. Gene editing confirmed the gene identity by producing mutant alleles that failed to complement existing thk1 mutants and that produced multiple-aleurone homozygous phenotypes. Thk1 encodes a homolog of NEGATIVE ON TATA-LESS1, a protein that acts as a scaffold for the CARBON CATABOLITE REPRESSION4-NEGATIVE ON TATA-LESS complex. This complex is highly conserved and essential in all eukaryotes for regulating a wide array of gene expression and cellular activities. Maize also harbors a duplicate locus, thick aleurone-like1, which likely accounts for the ability of thk1 mutants to form viable cells. Transcriptomic analysis indicated that THK1 regulates activities involving cell division, signaling, differentiation, and metabolism. Identification of thk1 provides an important new component of the DEK1 regulatory system that patterns cell fate in endosperm.

Transfer cells: what regulates the development of their intricate wall labyrinths?

Transfer cells (TCs) support high nutrient rates into, or at symplasmic discontinuities within, the plant body. Their transport capacity is conferred by an amplified plasma membrane surface area, enriched in nutrient transporters, supported on an intricately invaginated wall labyrinth (WL). Thus, development of the WL is at the heart of TC function. Enquiry has shifted from describing WL architecture and formation to discovering mechanisms regulating WL assembly. Experimental systems used to examine these phenomena are critiqued. Considerable progress has been made in identifying master regulators that commit stem cells to a TC fate (e.g. the maize Myeloblastosis (MYB)-related R1-type transcription factor) and signals that induce differentiated cells to undergo trans-differentiation to a TC phenotype (e.g. sugar, auxin and ethylene). In addition, signals that provide positional information for assembly of the WL include apoplasmic hydrogen peroxide and cytosolic Ca²⁺ plumes. The former switches on, and specifies the intracellular site for WL construction, while the latter creates subdomains to direct assembly of WL invaginations. Less is known about macromolecule species and their spatial organization essential for WL assembly. Emerging evidence points to a dependency on methyl-esterified homogalacturonan accumulation, unique patterns of cellulose and callose deposition and spatial positioning of arabinogalactan proteins.

Genome Wide Analysis of the Transcriptional Profiles in Different Regions of the Developing Rice Grains

BACKGROUND

Rice is an important food source for humans worldwide. Because of its nutritional and agricultural significance, a number of studies addressed various aspects of rice grain development and grain filling. Nevertheless, the molecular processes underlying grain filling and development, and in particular the contributions of different grain tissues to these processes, are not understood.

MAIN TEXT

Using RNA-sequencing, we profiled gene expression activity in grain tissues comprised of cross cells (CC), the nucellar epidermis (NE), ovular vascular trace (OVT), endosperm (EN) and the aleurone layer (AL). These tissues were dissected using laser capture microdissection (LCM) at three distinct grain development stages. The mRNA expression datasets offer comprehensive and new insights into the gene expression patterns in different rice grain tissues and their contributions to grain development. Comparative analysis of the different tissues revealed their similar and/or unique functions, as well as the spatio-temporal regulation of common and tissue-specific genes. The expression patterns of genes encoding hormones and transporters indicate an important role of the OVT tissue in metabolite transport during grain development. Gene co-expression network prediction on OVT-specific genes identified several distinct and common development-specific transcription factors. Further analysis of enriched DNA sequence motifs proximal to OVT-specific genes revealed known and novel DNA sequence motifs relevant to rice grain development.

CONCLUSION

Together, the dataset of gene expression in rice grain tissues is a novel and useful resource for further work to dissect the molecular and metabolic processes during rice grain development.

Loss of Function of an RNA Polymerase III Subunit Leads to Impaired Maize Kernel Development

Kernel size is an important factor determining grain yield. Although a number of genes affecting kernel development in maize (Zea mays) have been identified by analyzing kernel mutants, most of the corresponding mutants cannot be used in maize breeding programs due to low germination or incomplete seed development. Here, we characterized small kernel7, a recessive small-kernel mutant with a mutation in the gene encoding the second-largest subunit of RNA polymerase III (RNAPΙΙΙ; NRPC2). A frame shift in ZmNRPC2 leads to a premature stop codon, resulting in significantly reduced levels of transfer RNAs and 5S ribosomal RNA, which are transcribed by RNAPΙΙΙ. Loss-of-function nrpc2 mutants created by CRISPR/CAS9 showed significantly reduced kernel size due to altered endosperm cell size and number. ZmNRPC2 affects RNAPIII activity and the expression of genes involved in cell proliferation and endoreduplication to control kernel development via physically interacting with RNAPIII subunits RPC53 and AC40, transcription factor class C1 and Floury3. Notably, unlike the semidominant negative mutant floury3, which has defects in starchy endosperm, small kernel7 only affects kernel size but not the composition of kernel storage proteins. Our findings provide novel insights into the molecular network underlying maize kernel size, which could facilitate the genetic improvement of maize in the future.

Spatiotemporal Resolved Leaf Angle Establishment Improves Rice Grain Yield via Controlling Population Density

Leaf angle is mainly determined by the lamina joint (LJ) and contributes to ideal crop architecture for high yield. Here, we dissected five successive stages with distinct cytological features of LJs spanning organogenesis to leaf angle formation and obtained the underlying stage-specific mRNAs and small RNAs, which well explained the cytological dynamics during LJ organogenesis and leaf angle plasticity. Combining the gene coexpression correlation with high-throughput promoter analysis, we identified a set of transcription factors (TFs) determining the stage- and/or cytological structure-specific profiles. The functional studies of these TFs demonstrated that cytological dynamics determined leaf angle and that the knockout rice of these TFs with erect leaves significantly enhanced yield by maintaining the proper tiller number under dense planting. This work revealed the high-resolution mechanisms of how the cytological dynamics of LJ determined leaf erectness and served as a valuable resource to remodel rice architecture for high yield by controlling population density.

Transcriptomic Profiling of Young Cotyledons Response to Chilling Stress in Two Contrasting Cotton (Gossypium hirsutum L.) Genotypes at the Seedling Stage

Young cotyledons of cotton seedlings are most susceptible to chilling stress. To gain insight into the potential mechanism of cold tolerance of young cotton cotyledons, we conducted physiological and comparative transcriptome analysis of two varieties with contrasting phenotypes. The evaluation of chilling injury of young cotyledons among 74 cotton varieties revealed that H559 was the most tolerant and YM21 was the most sensitive. The physiological analysis found that the ROS scavenging ability was lower, and cell membrane damage was more severe in the cotyledons of YM21 than that of H559 under chilling stress. RNA-seq analysis identified a total of 44,998 expressed genes and 19,982 differentially expressed genes (DEGs) in young cotyledons of the two varieties under chilling stress. Weighted gene coexpression network analysis (WGCNA) of all DEGs revealed four significant modules with close correlation with specific samples. The GO-term enrichment analysis found that lots of genes in H559-specific modules were involved in plant resistance to abiotic stress. Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis revealed that pathways such as plant hormone signal transduction, MAPK signaling, and plant–pathogen interaction were related to chilling stress response. A total of 574 transcription factors and 936 hub genes in these modules were identified. Twenty hub genes were selected for qRT-PCR verification, revealing the reliability and accuracy of transcriptome data. These findings will lay a foundation for future research on the molecular mechanism of cold tolerance in cotyledons of cotton.

The dynamic transcriptome of waxy maize (Zea mays L. sinensis Kulesh) during seed development

BACKGROUND

Waxy maize (Zea mays L. sinensis Kulesh) is a mutant of maize (Zea mays L.) with a mutation at Waxy1 (Wx1) gene locus. The seed of waxy maize has higher viscosity compared to regular maize. By now, we know little about the expression patterns of genes that involved in the seed development of waxy maize.

OBJECTIVE

By analyzing the transcriptome data during waxy maize seed development, we attempt to dig out the genes that may influence the seed development of waxy maize.

METHODS

The seeds of waxy maize inbred line SWL01 from six phases after pollination were used to do RNA-seq. Bioinformatics methods were used to analyze the expression patterns of the expressed genes, to identify the genes involved in waxy maize seed development.

RESULTS

A total of 24,546 genes including 1611 transcription factors (TFs) were detected during waxy maize seed development. Coexpression analysis of expressed genes revealed the dynamic processes of waxy maize seed development. Particularly, 2457 genes including 177 TFs were specially expressed in waxy maize seed, some of which mainly involved in the process of seed dormancy and maturation. In addition, 2681, 5686, 4491, 4386, 3669 and 4624 genes were identified to be differential expressed genes (DEGs) at six phases compared to regular maize B73, and 113 DEGs among them may be key genes that lead the difference of seed development between waxy and regular maizes in milk stage.

CONCLUSION

In summary, we elucidated the expression patterns of expressed genes during waxy maize seed development globally. A series of genes that associated with seed development were identified in our research, which may provide an important resource for functional study of waxy maize seed development to help molecular assisted breeding.