Using giant scarlet runner bean embryos to uncover regulatory networks controlling suspensor gene activity

BSA-Seq and Transcriptomic Analysis Provide Candidate Genes Associated with Inflorescence Architecture and Kernel Orientation by Phytohormone Homeostasis in Maize

The developmental plasticity of the maize inflorescence depends on meristems, which directly affect reproductive potential and yield. However, the molecular roles of upper floral meristem (UFM) and lower floral meristem (LFM) in inflorescence and kernel development have not been fully elucidated. In this study, we characterized the reversed kernel1 (rk1) novel mutant, which contains kernels with giant embryos but shows normal vegetative growth like the wild type (WT). Total RNA was extracted from the inflorescence at three stages for transcriptomic analysis. A total of 250.16-Gb clean reads were generated, and 26,248 unigenes were assembled and annotated. Gene ontology analyses of differentially expressed genes (DEGs) detected in the sexual organ formation stage revealed that cell differentiation, organ development, phytohormonal responses and carbohydrate metabolism were enriched. The DEGs associated with the regulation of phytohormone levels and signaling were mainly expressed, including auxin (IAA), jasmonic acid (JA), gibberellins (GA), and abscisic acid (ABA). The transcriptome, hormone evaluation and immunohistochemistry observation revealed that phytohormone homeostasis were affected in rk1. BSA-Seq and transcriptomic analysis also provide candidate genes to regulate UFM and LFM development. These results provide novel insights for understanding the regulatory mechanism of UFM and LFM development in maize and other plants.

LCM and RNA-seq analyses revealed roles of cell cycle and translational regulation and homoeolog expression bias in cotton fiber cell initiation

Background

Cotton fibers provide a powerful model for studying cell differentiation and elongation. Each cotton fiber is a singular and elongated cell derived from epidermal-layer cells of a cotton seed. Efforts to understand this dramatic developmental shift have been impeded by the difficulty of separation between fiber and epidermal cells.

Results

Here we employed laser-capture microdissection (LCM) to separate these cell types. RNA-seq analysis revealed transitional differences between fiber and epidermal-layer cells at 0 or 2 days post anthesis. Specifically, down-regulation of putative cell cycle genes was coupled with upregulation of ribosome biosynthesis and translation-related genes, which may suggest their respective roles in fiber cell initiation. Indeed, the amount of fibers in cultured ovules was increased by cell cycle progression inhibitor, Roscovitine, and decreased by ribosome biosynthesis inhibitor, Rbin-1. Moreover, subfunctionalization of homoeologs was pervasive in fiber and epidermal cells, with expression bias towards 10% more D than A homoeologs of cell cycle related genes and 40–50% more D than A homoeologs of ribosomal protein subunit genes. Key cell cycle regulators were predicted to be epialleles in allotetraploid cotton. MYB-transcription factor genes displayed expression divergence between fibers and ovules. Notably, many phytohormone-related genes were upregulated in ovules and down-regulated in fibers, suggesting spatial-temporal effects on fiber cell development.

Conclusions

Fiber cell initiation is accompanied by cell cycle arrest coupled with active ribosome biosynthesis, spatial-temporal regulation of phytohormones and MYB transcription factors, and homoeolog expression bias of cell cycle and ribosome biosynthesis genes. These valuable genomic resources and molecular insights will help develop breeding and biotechnological tools to improve cotton fiber production.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07579-1.

Comparative analysis of embryo proper and suspensor transcriptomes in plant embryos with different morphologies

An important question is what genes govern the differentiation of plant embryos into suspensor and embryo proper regions following fertilization and division of the zygote. We compared embryo proper and suspensor transcriptomes of four plants that vary in embryo morphology within the suspensor region. We determined that genes encoding enzymes in several metabolic pathways leading to the formation of hormones, such as gibberellic acid, and other metabolites are up-regulated in giant scarlet runner bean and common bean suspensors. Genes involved in transport and Golgi body organization are up-regulated within the suspensors of these plants as well, strengthening the view that giant specialized suspensors serve as a hormone factory and a conduit for transferring substances to the developing embryo proper. By contrast, genes controlling transcriptional regulation, development, and cell division are up-regulated primarily within the embryo proper. Transcriptomes from less specialized soybean and Arabidopsis suspensors demonstrated that fewer genes encoding metabolic enzymes and hormones are up-regulated. Genes active in the embryo proper, however, are functionally similar to those active in scarlet runner bean and common bean embryo proper regions. We uncovered a set of suspensor- and embryo proper-specific transcription factors (TFs) that are shared by all embryos irrespective of morphology, suggesting that they are involved in early differentiation processes common to all plants. Chromatin immunoprecipitation sequencing (ChIP-Seq) experiments with scarlet runner bean and soybean WOX9, an up-regulated suspensor TF, gained entry into a regulatory network important for suspensor development irrespective of morphology.

Development of Embryo Suspensors for Five Genera of Crassulaceae with Special Emphasis on Plasmodesmata Distribution and Ultrastructure

The suspensor in the majority of angiosperms is an evolutionally conserved embryonic structure functioning as a conduit that connects ovule tissues with the embryo proper for nutrients and growth factors flux. This is the first study serving the purpose of investigating the correlation between suspensor types and plasmodesmata (PD), by the ultrastructure of this organ in respect of its full development. The special attention is paid to PD in representatives of Crassulaceae genera: Sedum, Aeonium, Monanthes, Aichryson and Echeveria. The contribution of the suspensor in transporting nutrients to the embryo was confirmed by the basal cell structure of the suspensor which produced, on the micropylar side of all genera investigated, a branched haustorium protruding into the surrounding ovular tissue and with wall ingrowths typically associated with cell transfer. The cytoplasm of the basal cell was rich in endoplasmic reticulum, mitochondria, dictyosomes, specialized plastids, microtubules, microbodies and lipid droplets. The basal cell sustained a symplasmic connection with endosperm and neighboring suspensor cells. Our results indicated the dependence of PD ultrastructure on the type of suspensor development: (i) simple PD are assigned to an uniseriate filamentous suspensor and (ii) PD with an electron-dense material are formed in a multiseriate suspensor. The occurrence of only one or both types of PD seems to be specific for the species but not for the genus. Indeed, in the two tested species of Sedum (with the distinct uniseriate/multiseriate suspensors), a diversity in the structure of PD depends on the developmental pattern of the suspensor. In all other genera (with the multiseriate type of development of the suspensor), the one type of electron-dense PD was observed.

Using Giant Scarlet Runner Bean (Phaseolus coccineus) Embryos to Dissect the Early Events in Plant Embryogenesis

The giant embryo of the scarlet runner bean (Phaseolus coccineus) has been used historically to investigate the molecular and developmental processes that control the early events of plant embryo development. In more recent years, our laboratory has been using scarlet runner bean embryos to uncover the genes and regulatory events that control embryo proper and suspensor region differentiation shortly after fertilization. In this chapter we describe methods that we have developed to isolate scarlet runner bean embryos at the globular stage of development, and capture embryo proper and suspensor regions by either hand dissection or laser capture microdissection (LCM) for use in downstream genomic analysis. These methods are also applicable for use in investigating the early events of common bean (Phaseolus vulgaris) embryo development, a close relative of scarlet runner bean, which also has a giant embryo in addition to a sequenced genome.

Transcriptome-derived investigation of biosynthesis of quinolizidine alkaloids in narrow-leafed lupin (Lupinus angustifolius L.) highlights candidate genes linked to iucundus locus

Unravelling the biosynthetic pathway of quinolizidine alkaloids (QAs), regarded as antinutritional compounds of narrow-leafed lupin (NLL) seeds, is fundamental to best exploit NLL as food or feed. We investigated 12 candidate genes connected to QA biosynthesis, selecting them by transcriptomic and genomic approaches, from the landscape of genes differentially expressed in leaves of the high- and low-alkaloid NLL accessions. Linkage analysis enabled the assessment of the location of the candidate genes in relation to iucundus, a major locus of unknown identity, that confers reduced QA content in seeds. The key finding was the identification of APETALA2/ethylene response transcription factor, RAP2-7, cosegregating with the iucundus locus and located within a region with highly significant QTLs that affect QA composition. We additionally identified a 4-hydroxy-tetrahydrodipicolinate synthase (DHDPS) gene involved in L-lysine biosynthesis as being closely linked to iucundus. The distributed location of other remaining candidates (including previously known QA genes) across different linkage groups, also indirectly supports the transcription factor as a possible regulator of lupin alkaloid biosynthesis. Our findings provide crucial insight into QA biosynthesis in NLL. Additionally, we evaluated and selected appropriate reference genes for qRT-PCRs to analyse the expression levels of QA genes in NLL.

Endoreplication and its consequences in the suspensor of Pisum sativum

KEY MESSAGE

DNA replication and continuous process of transcription during ongoing amitotic division accelerate the development of four-celled pea suspensor containing nuclei which create transient gradient of polyploidy necessary for correct embryo development. A suspensor, the link between embryo proper and surrounding tissues, differs significantly in size, morphology, and degree of polyploidy among the species. The suspensor of Pisum sativum consists of four polynuclear cells (two hemispherical and two elongated) formed in two layers. Their nuclei undergo endoreplication reaching, respectively, up to 256C and 128-256C DNA levels in its hemispherical and elongated parts. Our study shows that endoreplication first appears in the spherical part of the suspensor, and, subsequently, in the elongated one. At the next stages of suspensor development, the increase in DNA content takes place also in a similar order. Thus, despite simple construction of the suspensor, its development, supported by endoreplication, creates a certain gradient of polyploidy, which occurs in more extensive suspensors. Moreover, the rapid development of suspensor is supported both by the initiation of DNA replication prior to the completion of amitotic division of its polyploidal nuclei and by a continuous process of transcription, which is silenced by chromatin condensation throughout mitosis. Furthermore, the increase in DNA content correlates with the greater amount of transcripts; however, the multiplication of DNA copies does not entail an increase (but fluctuation) in the mean transcriptional activity of a particular nucleus during the next stages of suspensor development.

A protocol combining multiphoton microscopy and propidium iodide for deep 3D root meristem imaging in rice: application for the screening and identification of tissue-specific enhancer trap lines

BACKGROUND

The clear visualization of 3D organization at the cellular level in plant tissues is needed to fully understand plant development processes. Imaging tools allow the visualization of the main fluorophores and in vivo growth monitoring. Confocal microscopy coupled with the use of propidium iodide (PI) counter-staining is one of the most popular tools used to characterize the structure of root meristems in A. thaliana. However, such an approach is relatively ineffective in species with more complex and thicker root systems.

RESULTS

We adapted a PI counter-staining protocol to visualize the internal 3D architecture of rice root meristems using multiphoton microscopy. This protocol is simple and compatible with the main fluorophores (CFP, GFP and mCherry). The efficiency and applicability of this protocol were demonstrated by screening a population of 57 enhancer trap lines. We successfully characterized GFP expression in all of the lines and identified 5 lines with tissue-specific expression.

CONCLUSIONS

All of these resources are now available for the rice community and represent critical tools for future studies of root development.

A shared cis-regulatory module activates transcription in the suspensor of plant embryos

The mechanisms controlling the transcription of gene sets in specific regions of a plant embryo shortly after fertilization remain unknown. Previously, we showed that G564 mRNA, encoding a protein of unknown function, accumulates to high levels in the giant suspensor of both Scarlet Runner Bean (SRB) and Common Bean embryos, and a cis-regulatory module containing three unique DNA sequences, designated as the 10-bp, Region 2, and Fifth motifs, is required for G564 suspensor-specific transcription [Henry KF, et al. (2015) Plant Mol Biol 88:207-217; Kawashima T, et al. (2009) Proc Natl Acad Sci USA 106:3627-3632]. We tested the hypothesis that these motifs are also required for transcription of the SRB GA 20-oxidase gene, which encodes a gibberellic acid hormone biosynthesis enzyme and is coexpressed with G564 at a high level in giant bean suspensors. We used deletion and gain-of-function experiments in transgenic tobacco embryos to show that two GA 20-oxidase DNA regions are required for suspensor-specific transcription, one in the 5' UTR (+119 to +205) and another in the 5' upstream region (-341 to -316). Mutagenesis of sequences in these two regions determined that the cis-regulatory motifs required for G564 suspensor transcription are also required for GA 20-oxidase transcription within the suspensor, although the motif arrangement differs. Our results demonstrate the flexibility of motif positioning within a cis-regulatory module that activates gene transcription within giant bean suspensors and suggest that G564 and GA 20-oxidase comprise part of a suspensor gene regulatory network.

BLISS: Shining a light on lignification in plants

Lignin is a polyphenolic polymer of the plant cell wall formed by the oxidative polymerization of 3 main monomers called monolignols that give rise to the lignin H-, G- and S-units. Together with cellulose and hemicelluloses, lignin is a major component of plant biomass that is widely exploited by humans in numerous industrial processes. Despite recent advances in our understanding of monolignol biosynthesis, our current understanding of the spatio-temporal regulation of their transport and polymerization is more limited. In a recent publication, we have reported the development of an original Bioorthogonal Labeling Imaging Sequential Strategy (BLISS) that allows us to visualize the simultaneous incorporation dynamics of H and G monolignol reporters into lignifying cell walls of the flax stem. ¹¹ Here, we extend the application of this strategy to other plant organs such as roots and rapidly discuss some of the contributions and perspectives of this new technique for improving our understanding of the lignification process in plants.

Gene expression and metabolite profiling of gibberellin biosynthesis during induction of somatic embryogenesis in Medicago truncatula Gaertn

Gibberellins (GAs) are involved in the regulation of numerous developmental processes in plants including zygotic embryogenesis, but their biosynthesis and role during somatic embryogenesis (SE) is mostly unknown. In this study we show that during three week- long induction phase, when cells of leaf explants from non-embryogenic genotype (M9) and embryogenic variant (M9-10a) were forming the callus, all the bioactive gibberellins from non-13-hydroxylation (GA4, GA7) and 13-hydroxylation (GA1, GA5, GA3, GA6) pathways were present, but the contents of only a few of them differed between the tested lines. The GA53 and GA19 substrates synthesized by the 13-hydroxylation pathway accumulated specifically in the M9-10a line after the first week of induction; subsequently, among the bioactive gibberellins detected, only the content of GA3 increased and appeared to be connected with acquisition of embryogenic competence. We fully annotated 20 Medicago truncatula orthologous genes coding the enzymes which catalyze all the known reactions of gibberellin biosynthesis. Our results indicate that, within all the genes tested, expression of only three: MtCPS, MtGA3ox1 and MtGA3ox2, was specific to embryogenic explants and reflected the changes observed in GA53, GA19 and GA3 contents. Moreover, by analyzing expression of MtBBM, SE marker gene, we confirmed the inhibitory effect of manipulation in GAs metabolism, applying exogenous GA3, which not only impaired the production of somatic embryos, but also significantly decreased expression of this gene.

Gene expression and metabolite profiling of gibberellin biosynthesis during induction of somatic embryogenesis in Medicago truncatula Gaertn

Gibberellins (GAs) are involved in the regulation of numerous developmental processes in plants including zygotic embryogenesis, but their biosynthesis and role during somatic embryogenesis (SE) is mostly unknown. In this study we show that during three week- long induction phase, when cells of leaf explants from non-embryogenic genotype (M9) and embryogenic variant (M9-10a) were forming the callus, all the bioactive gibberellins from non-13-hydroxylation (GA4, GA7) and 13-hydroxylation (GA1, GA5, GA3, GA6) pathways were present, but the contents of only a few of them differed between the tested lines. The GA53 and GA19 substrates synthesized by the 13-hydroxylation pathway accumulated specifically in the M9-10a line after the first week of induction; subsequently, among the bioactive gibberellins detected, only the content of GA3 increased and appeared to be connected with acquisition of embryogenic competence. We fully annotated 20 Medicago truncatula orthologous genes coding the enzymes which catalyze all the known reactions of gibberellin biosynthesis. Our results indicate that, within all the genes tested, expression of only three: MtCPS, MtGA3ox1 and MtGA3ox2, was specific to embryogenic explants and reflected the changes observed in GA53, GA19 and GA3 contents. Moreover, by analyzing expression of MtBBM, SE marker gene, we confirmed the inhibitory effect of manipulation in GAs metabolism, applying exogenous GA3, which not only impaired the production of somatic embryos, but also significantly decreased expression of this gene.

A cis-regulatory module activating transcription in the suspensor contains five cis-regulatory elements

Little is known about the molecular mechanisms by which the embryo proper and suspensor of plant embryos activate specific gene sets shortly after fertilization. We analyzed the upstream region of the Scarlet Runner Bean (Phaseolus coccineus) G564 gene in order to understand how genes are activated specifically in the suspensor during early embryo development. Previously, we showed that a 54-bp fragment of the G564 upstream region is sufficient for suspensor transcription and contains at least three required cis-regulatory sequences, including the 10-bp motif (5'-GAAAAGCGAA-3'), the 10 bp-like motif (5'-GAAAAACGAA-3'), and Region 2 motif (partial sequence 5'-TTGGT-3'). Here, we use site-directed mutagenesis experiments in transgenic tobacco globular-stage embryos to identify two additional cis-regulatory elements within the 54-bp cis-regulatory module that are required for G564 suspensor transcription: the Fifth motif (5'-GAGTTA-3') and a third 10-bp-related sequence (5'-GAAAACCACA-3'). Further deletion of the 54-bp fragment revealed that a 47-bp fragment containing the five motifs (the 10-bp, 10-bp-like, 10-bp-related, Region 2 and Fifth motifs) is sufficient for suspensor transcription, and represents a cis-regulatory module. A consensus sequence for each type of motif was determined by comparing motif sequences shown to activate suspensor transcription. Phylogenetic analyses suggest that the regulation of G564 is evolutionarily conserved. A homologous cis-regulatory module was found upstream of the G564 ortholog in the Common Bean (Phaseolus vulgaris), indicating that the regulation of G564 is evolutionarily conserved in closely related bean species.