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

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

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

Transcriptome and Metabolome Analysis of Color Changes during Fruit Development of Pepper (Capsicum baccatum)

Fruit color is one of the most critical characteristics of pepper. In this study, pepper (Capsicum baccatum L.) fruits with four trans-coloring periods were used as experimental materials to explore the color conversion mechanism of pepper fruit. By transcriptome and metabolome analysis, we identified a total of 307 flavonoid metabolites, 68 carotenoid metabolites, 29 DEGs associated with flavonoid biosynthesis, and 30 DEGs related to carotenoid biosynthesis. Through WGCNA (weighted gene co-expression network analysis) analysis, positively correlated modules with flavonoids and carotenoids were identified, and hub genes associated with flavonoid and carotenoid synthesis and transport were anticipated. We identified Pinobanksin, Naringenin Chalcone, and Naringenin as key metabolites in the flavonoid biosynthetic pathway catalyzed by the key genes chalcone synthase (CHS CQW23_29123, CQW23_29380, CQW23_12748), cinnamic acid 4-hydroxylase (C4H CQW23_16085, CQW23_16084), cytochrome P450 (CYP450 CQW23_19845, CQW23_24900). In addition, phytoene synthase (PSY CQW23_09483), phytoene dehydrogenase (PDS CQW23_11317), zeta-carotene desaturase (ZDS CQW23_19986), lycopene beta cyclase (LYC CQW23_09027), zeaxanthin epoxidase (ZEP CQW23_05387), 9-cis-epoxycarotenoid dioxygenase (NCED CQW23_17736), capsanthin/capsorubin synthase (CCS CQW23_30321) are key genes in the carotenoid biosynthetic pathway, catalyzing the synthesis of key metabolites such as Phytoene, Lycopene, β-carotene and ε-carotene. We also found that transcription factor families such as p450 and NBARC could play important roles in the biosynthesis of flavonoids and carotenoids in pepper fruits. These results provide new insights into the interaction mechanisms of genes and metabolites involved in the biosynthesis of flavonoids and carotenoids in pepper fruit leading to color changes in pepper fruit.

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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.