Transcriptomic profiling identifies differentially expressed genes associated with programmed cell death of nucellar cells in Ginkgo biloba L

Arabidopsis cell suspension culture and RNA sequencing reveal regulatory networks underlying plant-programmed cell death

Programmed cell death (PCD) facilitates selective, genetically controlled elimination of redundant, damaged, or infected cells. In plants, PCD is often an essential component of normal development and can mediate responses to abiotic and biotic stress stimuli. However, studying the transcriptional regulation of PCD is hindered by difficulties in sampling small groups of dying cells that are often buried within the bulk of living plant tissue. We addressed this challenge by using RNA sequencing and Arabidopsis thaliana suspension cells, a model system that allows precise monitoring of PCD rates. The use of three PCD-inducing treatments (salicylic acid, heat, and critical dilution), in combination with three cell death modulators (3-methyladenine, lanthanum chloride, and conditioned medium), enabled isolation of candidate core- and stimuli-specific PCD genes, inference of underlying regulatory networks and identification of putative transcriptional regulators of PCD in plants. This analysis underscored a disturbance of the cell cycle and mitochondrial retrograde signaling, and repression of pro-survival stress responses, as key elements of the PCD-associated transcriptional signature. Further, phenotyping of Arabidopsis T-DNA insertion mutants in selected candidate genes validated the potential of generated resources to identify novel genes involved in plant PCD pathways and/or stress tolerance.

A transcriptomic evaluation of the mechanism of programmed cell death of the replaceable bud in Chinese chestnut

Previous studies suggest that the senescence and death of the replaceable bud of the Chinese chestnut cultivar (cv.) "Tima Zhenzhu" involves programmed cell death (PCD). However, the molecular network regulating replaceable bud PCD is poorly characterized. Here, we performed transcriptomic profiling on the chestnut cv. "Tima Zhenzhu" replaceable bud before (S20), during (S25), and after (S30) PCD to unravel the molecular mechanism underlying the PCD process. A total of 5,779, 9,867, and 2,674 differentially expressed genes (DEGs) were discovered upon comparison of S20 vs S25, S20 vs S30, and S25 vs S30, respectively. Approximately 6,137 DEGs common to at least two comparisons were selected for gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses to interrogate the main corresponding biological functions and pathways. GO analysis showed that these common DEGs could be divided into three functional categories, including 15 cellular components, 14 molecular functions, and 19 biological processes. KEGG analysis found that "plant hormone signal transduction" included 93 DEGs. Overall, 441 DEGs were identified as related to the process of PCD. Most of these were found to be genes associated with ethylene signaling, as well as the initiation and execution of various PCD processes.

Characterization of the Calmodulin/Calmodulin-like Protein (CAM/CML) Family in Ginkgo biloba, and the Influence of an Ectopically Expressed GbCML Gene (Gb_30819) on Seedling and Fruit Development of Transgenic Arabidopsis

Calmodulins (CAMs) and calmodulin-like proteins (CMLs) can participate in the regulation of various physiological processes via sensing and decoding Ca²⁺ signals. To reveal the characteristics of the CAM/CML family in Ginkgo biloba, a comprehensive analysis was performed at the genome-wide level. A total of 26 CAMs/CMLs, consisting of 5 GbCAMs and 21 GbCMLs, was identified on 11 out of 12 chromosomes in G. biloba. They displayed a certain degree of multiplicity in their sequences, albeit with conserved EF hands. Collinearity analysis suggested that tandem rather than segmental or whole-genome duplications were likely to play roles in the evolution of the Ginkgo CAM/CML family. Furthermore, GbCAMs/GbCMLs were grouped into higher, lower, and moderate expression in magnitude. The cis-acting regulatory elements involved in phytohormone-responsiveness within GbCAM/GbCML promotors may explain their varied expression profiles. The ectopic expression of a GbCML gene (Gb_30819) in transgenic Arabidopsis led to phenotypes with significantly shortened root length and seedling height, and decreased yields of both pods and seeds. Moreover, an electrophoresis mobility shift assay demonstrated the Ca²⁺-binding activity of Gb_30819 in vitro. Altogether, these results contribute to insights into the characteristics of the evolution and expression of GbCAMs/GbCMLs, as well as evidence for Ca²⁺-CAM/CML pathways functioning within the ancient gymnosperm G. biloba.

Transcriptional Profiling of Resistant and Susceptible Cultivars of Grapevine (Vitis L.) Reveals Hypersensitive Responses to Plasmopara viticola

Grapevine downy mildew is the most serious disease of grapevine cultivars that affects the rate of resistance/susceptibility to Plasmopara viticola. In this study, we used the susceptible cultivar "Zitian Seedless" and the resistant cultivar "Kober 5BB" as materials to determine the transcriptome differences and phenotypes of the leaves after inoculation with downy mildew. The differences in microstructures and molecular levels were compared and analyzed. Fluorescence staining and microscopic observations confirmed that hypersensitive cell death occurred around the stomata in "Kober 5BB" infected by downy mildew zoospores. Meanwhile, transcriptomic profiling indicated that there were 11,713 and 6,997 gene expression differences between the resistant and susceptible cultivars at 72 h after inoculation when compared to control (0 h), respectively. The differentially expressed genes of the two cultivars are significantly enriched in different pathways, including response to plant-pathogen interaction, mitogen-activated protein kinase (MAPK) signaling pathway, plant hormone signal transduction, phenylpropanoid, and flavonoid biosynthesis. Furthermore, the results of functional enrichment analysis showed that H₂O₂ metabolism, cell death, reactive oxygen response, and carbohydrate metabolism are also involved in the defense response of "Kober 5BB," wherein a total of 322 key genes have been identified. The protein interaction network showed that metacaspases (MCAs), vacuolar processing enzymes (VPEs), and Papain-like cysteine proteases (PLCPs) play an important role in the execution of hypersensitive responses (HR). In conclusion, we demonstrated that HR cell death is the key strategy in the process of grape defense against downy mildew, which may be mediated or activated by Caspase-like proteases.

Editing of the starch branching enzyme gene SBE2 generates high-amylose storage roots in cassava

The production of high-amylose cassava through CRISPR/Cas9-mediated mutagenesis of the starch branching enzyme gene SBE2 was firstly achieved. High-amylose cassava (Manihot esculenta Crantz) is desirable for starch industrial applications and production of healthier processed food for human consumption. In this study, we report the production of high-amylose cassava through CRISPR/Cas9-mediated mutagenesis of the starch branching enzyme 2 (SBE2). Mutations in two targeted exons of SBE2 were identified in all regenerated plants; these mutations, which included nucleotide insertions, and short or long deletions in the SBE2 gene, were classified into eight mutant lines. Three mutants, M6, M7 and M8, with long fragment deletions in the second exon of SBE2 showed no accumulation of SBE2 protein. After harvest from the field, significantly higher amylose (up to 56% in apparent amylose content) and resistant starch (up to 35%) was observed in these mutants compared with the wild type, leading to darker blue coloration of starch granules after quick iodine staining and altered starch viscosity with a higher pasting temperature and peak time. Further ¹H-NMR analysis revealed a significant reduction in the degree of starch branching, together with fewer short chains (degree of polymerization [DP] 15-25) and more long chains (DP>25 and especially DP>40) of amylopectin, which indicates that cassava SBE2 catalyzes short chain formation during amylopectin biosynthesis. Transition from A- to B-type crystallinity was also detected in the starches. Our study showed that CRISPR/Cas9-mediated mutagenesis of starch biosynthetic genes in cassava is an effective approach for generating novel varieties with valuable starch properties for food and industrial applications.

Nucleotide-binding leucine-rich repeat proteins: a missing link in controlling cell fate and plant adaptation to hostile environment?

Programmed cell death is a tightly regulated genetically controlled process that leads to cell suicide and eliminates cells that are either no longer needed or damaged/harmful. Nucleotide-binding leucine-rich repeat proteins have recently emerged as a novel class of Ca2+-permeable channels that operate in plant immune responses. This viewpoint argues that the unique structure of this channel, its permeability to other cations, and specificity of its operation make it an ideal candidate to mediate cell signaling and adaptive responses not only to pathogens but also to a broad range of abiotic stress factors.

Molecular mechanism of mulberry response to drought stress revealed by complementary transcriptomic and iTRAQ analyses

BACKGROUND

The use of mulberry leaves has long been limited to raising silkworms, but with the continuous improvement of mulberry (Morus alba) resource development and utilization, various mulberry leaf extension products have emerged. However, the fresh leaves of mulberry trees have a specific window of time for picking and are susceptible to adverse factors, such as drought stress. Therefore, exploring the molecular mechanism by which mulberry trees resist drought stress and clarifying the regulatory network of the mulberry drought response is the focus of the current work.

RESULTS

In this study, natural and drought-treated mulberry grafted seedlings were used for transcriptomic and proteomic analyses (CK vs. DS9), aiming to clarify the molecular mechanism of the mulberry drought stress response. Through transcriptome and proteome sequencing, we identified 9889 DEGs and 1893 DEPs enriched in stress-responsive GO functional categories, such as signal transducer activity, antioxidant activity, and transcription regulator activity. KEGG enrichment analysis showed that a large number of codifferentially expressed genes were enriched in flavonoid biosynthesis pathways, hormone signalling pathways, lignin metabolism and other pathways. Through subsequent cooperation analysis, we identified 818 codifferentially expressed genes in the CK vs. DS9 comparison group, including peroxidase (POD), superoxide dismutase (SOD), aldehyde dehydrogenase (ALDHs), glutathione s-transferase (GST) and other genes closely related to the stress response. In addition, we determined that the mulberry gene MaWRKYIII8 (XP_010104968.1) underwent drought- and abscisic acid (ABA)-induced expression, indicating that it may play an important role in the mulberry response to drought stress.

CONCLUSIONS

Our research shows that mulberry can activate proline and ABA biosynthesis pathways and produce a large amount of proline and ABA, which improves the drought resistance of mulberry. MaWRKYIII8 was up-regulated and induced by drought and exogenous ABA, indicating that MaWRKYIII8 may be involved in the mulberry response to drought stress. These studies will help us to analyse the molecular mechanism underlying mulberry drought tolerance and provide important gene information and a theoretical basis for improving mulberry drought tolerance through molecular breeding in the future.

Structure, synthesis, biosynthesis, and activity of the characteristic compounds from Ginkgo biloba L

Covering: 1928-2021Ginkgo biloba L. is one of the most distinctive plants to have emerged on earth and has no close living relatives. Owing to its phylogenetic divergence from other plants, G. biloba contains many compounds with unique structures that have served to broaden the chemical diversity of herbal medicine. Examples of such compounds include terpene trilactones (ginkgolides), acylated flavonol glycosides (ginkgoghrelins), biflavones (ginkgetin), ginkgotides and ginkgolic acids. The extract of G. biloba leaf is used to prevent and/or treat cardiovascular diseases, while many ginkgo-derived compounds are currently at various stages of preclinical and clinical trials worldwide. The global annual sales of G. biloba products are estimated to total US$10 billion. However, the content and purity of the active compounds isolated by traditional methods are usually low and subject to varying environmental factors, making it difficult to meet the huge demand of the international market. This highlights the need to develop new strategies for the preparation of these characteristic compounds from G. biloba. In this review, we provide a detailed description of the structures and bioactivities of these compounds and summarize the recent research on the development of strategies for the synthesis, biosynthesis, and biotechnological production of the characteristic terpenoids, flavonoids, and alkylphenols/alkylphenolic acids of G. biloba. Our aim is to provide an important point of reference for all scientists who research ginkgo-related compounds for medicinal or other purposes.

Evolution and patterning of the ovule in seed plants

The ovule and its developmental successor, the seed, together represent a highly characteristic feature of seed plants that has strongly enhanced the reproductive and dispersal potential of this diverse group of taxa. Ovules encompass multiple tissues that perform various roles within a highly constrained space, requiring a complex cascade of genes that generate localized cell proliferation and programmed cell death during different developmental stages. Many heritable morphological differences among lineages reflect relative displacement of these tissues, but others, such as the second (outer) integuments of angiosperms and Gnetales, represent novel and apparently profound and independent innovations. Recent studies, mostly on model taxa, have considerably enhanced our understanding of gene expression in the ovule. However, understanding its evolutionary history requires a comparative and phylogenetic approach that is problematic when comparing extant angiosperms not only with phylogenetically distant extant gymnosperms but also with taxa known only from fossils. This paper reviews ovule characters across a phylogenetically broad range of seed plants in a dynamic developmental context. It discusses both well-established and recent theories of ovule and seed evolution and highlights potential gaps in comparative data that will usefully enhance our understanding of evolutionary transitions and developmental mechanisms.

Insights into Plant Programmed Cell Death Induced by Heavy Metals-Discovering a Terra Incognita

Programmed cell death (PCD) is a process that plays a fundamental role in plant development and responses to biotic and abiotic stresses. Knowledge of plant PCD mechanisms is still very scarce and is incomparable to the large number of studies on PCD mechanisms in animals. Quick and accurate assays, e.g., the TUNEL assay, comet assay, and analysis of caspase-like enzyme activity, enable the differentiation of PCD from necrosis. Two main types of plant PCD, developmental (dPCD) regulated by internal factors, and environmental (ePCD) induced by external stimuli, are distinguished based on the differences in the expression of the conserved PCD-inducing genes. Abiotic stress factors, including heavy metals, induce necrosis or ePCD. Heavy metals induce PCD by triggering oxidative stress via reactive oxygen species (ROS) overproduction. ROS that are mainly produced by mitochondria modulate phytotoxicity mechanisms induced by heavy metals. Complex crosstalk between ROS, hormones (ethylene), nitric oxide (NO), and calcium ions evokes PCD, with proteases with caspase-like activity executing PCD in plant cells exposed to heavy metals. This pathway leads to very similar cytological hallmarks of heavy metal induced PCD to PCD induced by other abiotic factors. The forms, hallmarks, mechanisms, and genetic regulation of plant ePCD induced by abiotic stress are reviewed here in detail, with an emphasis on plant cell culture as a suitable model for PCD studies. The similarities and differences between plant and animal PCD are also discussed.

The regulatory framework of developmentally programmed cell death in floral organs: A review

Developmentally programmed cell death (dPCD) is a tightly controlled biological process. In recent years, vital roles of dPCD on regulating floral organ growth and development have been reported. It is well known that flower is an essential organ for reproduction and a turning point of plants' life cycle. Hence, uncovering the complex molecular networks which regulates dPCD processes in floral organs is utmost important. So far, our understanding of dPCD on floral organ growth and development is just starting. Herein, we summarize the important factors that involved in the tapetal degeneration, pollen tube rupture, receptive synergid cell death, nucellar degradation, and antipodal cell degradation. Meanwhile, the known factors that involved in transmitting tract formation and self-incompatibility-induced PCD were also introduced. Furthermore, the genes that associated with anther dehiscence and petal senescence and abscission were reviewed as well. The functions of various types of factors involved in floral dPCD processes are highlighted principally. The regulatory panorama described here can provide us some insights about flower-specific dPCD process.

Genome-Wide Identification and Coexpression Network Analysis of DNA Methylation Pathway Genes and Their Differentiated Functions in Ginkgo biloba L

DNA methylation plays a vital role in diverse biological processes. DNA methyltransferases (DNMTs) genes and RNA-directed DNA methylation (RdDM)-related genes are key genes responsible for establishing and maintaining genome DNA methylation in plants. In the present study, we systematically identified nine GbDNMTs in Ginkgo biloba, including the three common families of GbMET1a/1b, GbCMT2, and GbDRMa/b/2a/2b/2c, and a fourth family—GbDNMT3—which is absent in most angiosperms. We also identified twenty RdDM-related genes, including four GbDCLs, six GbAGOs, and ten GbRDRs. Expression analysis of the genes showed the different patterns of individual genes, and 15 of 29 genes displayed expression change under five types of abiotic stress. Gene coexpression analysis and weighted gene co-expression network analysis (WGCNA) using 126 public transcriptomic datasets revealed that these genes were clustered into two groups. In group I, genes covered members from all six families which were preferentially expressed in the ovulate strobile and fruit. A gene ontology (GO) enrichment analysis of WGCNA modules indicated that group I genes were most correlated with the biological process of cell proliferation. Group II only consisted of RdDM-related genes, including GbDRMs, GbAGOs, and GbRDRs, but no GbDCLs, and these genes were specifically expressed in the cambium, suggesting that they may function in a dicer-like (DCL)-independent RdDM pathway in specific tissues. The gene module related to group II was most enriched in signal transduction, cell communication, and the response to the stimulus. These results demonstrate that gene family members could be conserved or diverged across species, and multi-member families in the same pathway may cluster into different modules to function differentially. The study provides insight into the DNA methylation genes and their possible functions in G. biloba, laying a foundation for the further study of DNA methylation in gymnosperms.

Programmed cell death (PCD) control in plants: New insights from the Arabidopsis thaliana deathosome

Programmed cell death (PCD) is a genetically controlled process that leads to cell suicide in both eukaryotic and prokaryotic organisms. In plants PCD occurs during development, defence response and when exposed to adverse conditions. PCD acts controlling the number of cells by eliminating damaged, old, or unnecessary cells to maintain cellular homeostasis. Unlike in animals, the knowledge about PCD in plants is limited. The molecular network that controls plant PCD is poorly understood. Here we present a review of the current mechanisms involved with the genetic control of PCD in plants. We also present an updated version of the AtLSD1 deathosome, which was previously proposed as a network controlling HR-mediated cell death in Arabidopsis thaliana. Finally, we discuss the unclear points and open questions related to the AtLSD1 deathosome.

Full-length sequencing of ginkgo transcriptomes for an in-depth understanding of flavonoid and terpenoid trilactone biosynthesis

Ginkgo biloba L. is regarded as the most ancient living tree, and its kernel has been used as a traditional Chinese medicine for more than 2,000 years. The leaf extracts of this tree have been among the bestselling herbal remedies in Western countries since the last century. To understand the biosynthesis of the pharmacologically active ingredients in G. biloba, flavonoids and terpenoid trilactones (TTLs), we sequenced the transcriptomes of G. biloba leaves, kernels and testae with Iso-Seq and RNA-Seq technologies and obtained 152,524 clean consensus reads. When these reads were used to improve the annotation of the G. biloba genome, 4,856 novel genes, 25,583 new isoforms of previously annotated genes and 4,363 lncRNAs were discovered. Gene ontology and Kyoto Encyclopedia of Genes and Genomes analyses indicated that genes involved in growth, regulation and response to stress were more likely to be regulated by alternative splicing (AS) or alternative polyadenylation (APA), which represent the two most important posttranscriptional regulation mechanisms. It was found that some of the characterized genes involved in the biosynthesis of flavonoids and TTLs were also possibly regulated by AS and APA. Using phylogenetic and gene expression pattern analyses, some candidate genes for the biosynthesis of flavonoids and TTLs were screened. After qRT-PCR validation, the final candidate genes for flavonoid biosynthesis included three UDP-glycosyltransferases and one MYB transcription factor, while the candidate genes for TTL biosynthesis included two cytochrome P450 and one WRKY transcription factor. Our study suggested that Iso-Seq may play an important role in improving genome annotation, elucidating AS and APA mechanisms and discovering candidate genes involved in the biosynthesis of some secondary metabolites.

Transcriptomic analysis of early fruit development in Chinese white pear (Pyrus bretschneideri Rehd.) and functional identification of PbCCR1 in lignin biosynthesis

BACKGROUND

The content of stone cells and lignin is one of the key factors affecting the quality of pear fruit. In a previous study, we determined the developmental regularity of stone cells and lignin in 'Dangshan Su' pear fruit 15-145 days after pollination (DAP). However, the development of fruit stone cells and lignin before 15 DAP has not been heavily researched.

RESULTS

In this study, we found that primordial stone cells began to appear at 7 DAP and that the fruit had formed a large number of stone cells at 15 DAP. Subsequently, transcriptome sequencing was performed on fruits at 0, 7, and 15 DAP and identified 3834 (0 vs. 7 DAP), 4049 (7 vs. 15 DAP) and 5763 (0 vs. 15 DAP) DEGs. During the 7-15 DAP period, a large number of key enzyme genes essential for lignin biosynthesis are gradually up-regulated, and their expression pattern is consistent with the accumulation of lignin in this period. Further analysis found that the biosynthesis of S-type lignin in 'Dangshan Su' pear does not depend on the catalytic activity of PbSAD but is primarily generated by the catalytic activity of caffeoyl-CoA through CCoAOMT, CCR, F5H, and CAD. We cloned PbCCR1, 2 and analysed their functions in Chinese white pear lignin biosynthesis. PbCCR1 and 2 have a degree of functional redundancy; both demonstrate the ability to participate in lignin biosynthesis. However, PbCCR1 may be the major gene for lignin biosynthesis, while PbCCR2 has little effect on lignin biosynthesis.

CONCLUSIONS

Our results revealed that 'Dangshan Su' pear began to form a large number of stone cells and produce lignin after 7 DAP and mainly accumulated materials from 0 to 7 DAP. PbCCR1 is mainly involved in the biosynthesis of lignin in 'Dangshan Su' pear and plays a positive role in lignin biosynthesis.