Comparative transcriptomic analysis reveals conserved programmes underpinning organogenesis and reproduction in land plants

The RALF signaling pathway regulates cell wall integrity during pollen tube growth in maize

Autocrine signaling pathways regulated by RAPID ALKALINIZATION FACTORs (RALFs) control cell wall integrity during pollen tube germination and growth in Arabidopsis (Arabidopsis thaliana). To investigate the role of pollen-specific RALFs in another plant species, we combined gene expression data with phylogenetic and biochemical studies to identify candidate orthologs in maize (Zea mays). We show that Clade IB ZmRALF2/3 mutations, but not Clade III ZmRALF1/5 mutations, cause cell wall instability in the sub-apical region of the growing pollen tube. ZmRALF2/3 are mainly located in the cell wall and are partially able to complement the pollen germination defect of their Arabidopsis orthologs AtRALF4/19. Mutations in ZmRALF2/3 compromise pectin distribution patterns leading to altered cell wall organization and thickness culminating in pollen tube burst. Clade IB, but not Clade III ZmRALFs, strongly interact as ligands with the pollen-specific Catharanthus roseus RLK1-like (CrRLK1L) receptor kinases Z. mays FERONIA-like (ZmFERL) 4/7/9, LORELEI-like glycosylphosphatidylinositol-anchor (LLG) proteins Z. mays LLG 1 and 2 (ZmLLG1/2), and Z. mays pollen extension-like (PEX) cell wall proteins ZmPEX2/4. Notably, ZmFERL4 outcompetes ZmLLG2 and ZmPEX2 outcompetes ZmFERL4 for ZmRALF2 binding. Based on these data, we suggest that Clade IB RALFs act in a dual role as cell wall components and extracellular sensors to regulate cell wall integrity and thickness during pollen tube growth in maize and probably other plants.

Decoding the leaf apical meristem of Guarea glabra Vahl (Meliaceae): insight into the evolution of indeterminate pinnate leaves

In seed plants, growth of shoots and roots is indeterminate, while leaves are typically determinate organs that cease to grow after a certain developmental stage. This is due to the characteristics of the leaf meristem, where cell proliferation activity is retained only for a limited period. However, several plants exhibit indeterminacy in their leaves, exemplified by the pinnate compound leaves of Guarea and Chisocheton genera in the Meliaceae family. In these plants, the leaf meristem at the tip of the leaf retains meristematic activity and produces leaflets over years, resulting in a single leaf that resembles a twig. The molecular mechanism underlying the indeterminate leaf meristem of these plants has not been examined. In this research, we used Guarea glabra as a model to investigate the development of indeterminate pinnate leaves. Transcriptome analyses revealed that the gene expression profile in leaf apex tissue differed from that in the shoot apex. However, a class 1 KNOTTED-LIKE HOMEOBOX (KNOX1) gene which is lost in Brassicaceae was highly expressed in both tissues. We established an in situ hybridisation system for this species using Technovit 9100 to analyse the spatial expression patterns of genes. We revealed that the leaf meristematic region of G. glabra expresses KNOX1, LEAFY and ANGUSTIFORIA3 simultaneously, suggesting the involvement of these genes in the indeterminacy of the leaf meristem.

How to utilize comparative transcriptomics to dissect morphological diversity in plants

Comparative transcriptomics has emerged as a powerful approach that allows us to unravel the genetic basis of organ morphogenesis and its diversification processes during evolution. However, the application of comparative transcriptomics in studying plant morphological diversity addresses challenges such as identifying homologous gene pairs, selecting appropriate developmental stages for comparison, and extracting biologically meaningful networks. Methods such as phylostratigraphy, clustering, and gene co-expression networks are explored to identify functionally equivalent genes, align developmental stages, and uncover gene regulatory relationships. In the current review, we highlight the importance of these approaches in overcoming the complexity of plant genomes, the impact of heterochrony on stage alignment, and the integration of gene networks with additional data for a comprehensive understanding of morphological evolution.

Design principles for synthetic control systems to engineer plants

KEY MESSAGE

Synthetic control systems have led to significant advancement in the study and engineering of unicellular organisms, but it has been challenging to apply these tools to multicellular organisms like plants. The ability to predictably engineer plants will enable the development of novel traits capable of alleviating global problems, such as climate change and food insecurity. Engineering predictable multicellular phenotypes will require the development of synthetic control systems that can precisely regulate how the information encoded in genomes is translated into phenotypes. Many efficient control systems have been developed for unicellular organisms. However, it remains challenging to use such tools to study or engineer multicellular organisms. Plants are a good chassis within which to develop strategies to overcome these challenges, thanks to their capacity to withstand large-scale reprogramming without lethality. Additionally, engineered plants have great potential for solving major societal problems. Here we briefly review the progress of control system development in unicellular organisms, and how that information can be leveraged to characterize control systems in plants. Further, we discuss strategies for developing control systems designed to regulate the expression of transgenes or endogenous loci and generate dosage-dependent or discrete traits. Finally, we discuss the utility that mathematical models of biological processes have for control system deployment.

Topological data analysis reveals a core gene expression backbone that defines form and function across flowering plants

Since they emerged approximately 125 million years ago, flowering plants have evolved to dominate the terrestrial landscape and survive in the most inhospitable environments on earth. At their core, these adaptations have been shaped by changes in numerous, interconnected pathways and genes that collectively give rise to emergent biological phenomena. Linking gene expression to morphological outcomes remains a grand challenge in biology, and new approaches are needed to begin to address this gap. Here, we implemented topological data analysis (TDA) to summarize the high dimensionality and noisiness of gene expression data using lens functions that delineate plant tissue and stress responses. Using this framework, we created a topological representation of the shape of gene expression across plant evolution, development, and environment for the phylogenetically diverse flowering plants. The TDA-based Mapper graphs form a well-defined gradient of tissues from leaves to seeds, or from healthy to stressed samples, depending on the lens function. This suggests that there are distinct and conserved expression patterns across angiosperms that delineate different tissue types or responses to biotic and abiotic stresses. Genes that correlate with the tissue lens function are enriched in central processes such as photosynthetic, growth and development, housekeeping, or stress responses. Together, our results highlight the power of TDA for analyzing complex biological data and reveal a core expression backbone that defines plant form and function.

A bHLH heterodimer regulates germ cell differentiation in land plant gametophytes

Land plants are a monophyletic group of photosynthetic eukaryotes that diverged from streptophyte algae about 470 million years ago. During both the alternating haploid and diploid stages of the life cycle, land plants form multicellular bodies.1,2,3,4 The haploid multicellular body (gametophyte) produces progenitor cells that give rise to gametes and the reproductive organs.5,6,7,8 In the liverwort Marchantia polymorpha, differentiation of the initial cells of gamete-producing organs (gametangia) from the gametophyte is regulated by MpBONOBO (MpBNB), a member of the basic helix-loop-helix (bHLH) transcription factor subfamily VIIIa. In Arabidopsis thaliana, specification of generative cells in developing male gametophytes (pollen) requires redundant action of BNB1 and BNB2.9 Subfamily XI bHLHs, such as LOTUS JAPONICUS ROOTHAIRLESS LIKE1 (LRL1)/DEFECTIVE REGION OF POLLEN1 (DROP1) and LRL2/DROP2 in A. thaliana and the single LRL/DROP protein MpLRL in M. polymorpha, are the evolutionarily conserved regulators of rooting system development.10 Although the role of LRL1/DROP1 and LRL2/DROP2 in gametogenesis remains unclear, their loss leads to the formation of abnormal pollen devoid of sperm cells.11 Here, we show that BNBs and LRL/DROPs co-localize to gametophytic cell nuclei and form heterodimers. LRL1/DROP1 and LRL2/DROP2 act redundantly to regulate BNB expression for generative cell specification in A. thaliana after asymmetric division of the haploid microspore. MpLRL is required for differentiation of MpBNB-expressing gametangium initial cells in M. polymorpha gametophytes. Our findings suggest that broadly expressed LRL/DROP stabilizes BNB expression, leading to the formation of an evolutionarily conserved bHLH heterodimer, which regulates germ cell differentiation in the haploid gametophyte of land plants.

Conserved and non-conserved RNA-target modules in plants: lessons for a better understanding of Marchantia development

A wide variety of functional regulatory non-coding RNAs (ncRNAs) have been identified as essential regulators of plant growth and development. Depending on their category, ncRNAs are not only involved in modulating target gene expression at the transcriptional and post-transcriptional levels but also are involved in processes like RNA splicing and RNA-directed DNA methylation. To fulfill their molecular roles properly, ncRNAs must be precisely processed by multiprotein complexes. In the case of small RNAs, DICER-LIKE (DCL) proteins play critical roles in the production of mature molecules. Land plant genomes contain at least four distinct classes of DCL family proteins (DCL1-DCL4), of which DCL1, DCL3 and DCL4 are also present in the genomes of bryophytes, indicating the early divergence of these genes. The liverwort Marchantia polymorpha has become an attractive model species for investigating the evolutionary history of regulatory ncRNAs and proteins that are responsible for ncRNA biogenesis. Recent studies on Marchantia have started to uncover the similarities and differences in ncRNA production and function between the basal lineage of bryophytes and other land plants. In this review, we summarize findings on the essential role of regulatory ncRNAs in Marchantia development. We provide a comprehensive overview of conserved ncRNA-target modules among M. polymorpha, the moss Physcomitrium patens and the dicot Arabidopsis thaliana, as well as Marchantia-specific modules. Based on functional studies and data from the literature, we propose new connections between regulatory pathways involved in Marchantia's vegetative and reproductive development and emphasize the need for further functional studies to understand the molecular mechanisms that control ncRNA-directed developmental processes.

oggmap: a Python package to extract gene ages per orthogroup and link them with single-cell RNA data

SUMMARY

For model species, single-cell RNA-based cell atlases are available. A good cell atlas includes all major stages in a species' ontogeny, and soon, they will be standard even for nonmodel species. Here, we propose a Python package called oggmap, which allows for the easy extraction of an orthomap (gene ages per orthogroup) for any given query species from OrthoFinder and other gene family data resources, like homologous groups from eggNOG or PLAZA. oggmap provides extracted gene ages for more than thousand eukaryotic species which can be further used to calculate gene age-weighted expression data from scRNA sequencing objects using the Python Scanpy toolkit. Not limited to one transcriptome evolutionary index, oggmap can visualize the individual gene category (e.g. age class, nucleotide diversity bin) and their corresponding expression profiles to investigate scRNA-based cell type assignments in an evolutionary context.

AVAILABILITY AND IMPLEMENTATION

oggmap source code is available at https://github.com/kullrich/oggmap, documentation is available at https://oggmap.readthedocs.io/en/latest/. oggmap can be installed via PyPi or directly used via a docker container.

Comparative phylotranscriptomics reveals ancestral and derived root nodule symbiosis programmes

Symbiotic interactions such as the nitrogen-fixing root nodule symbiosis (RNS) have structured ecosystems during the evolution of life. Here we aimed at reconstructing ancestral and intermediate steps that shaped RNS observed in extant flowering plants. We compared the symbiotic transcriptomic responses of nine host plants, including the mimosoid legume Mimosa pudica for which we assembled a chromosome-level genome. We reconstructed the ancestral RNS transcriptome composed of most known symbiotic genes together with hundreds of novel candidates. Cross-referencing with transcriptomic data in response to experimentally evolved bacterial strains with gradual symbiotic proficiencies, we found the response to bacterial signals, nodule infection, nodule organogenesis and nitrogen fixation to be ancestral. By contrast, the release of symbiosomes was associated with recently evolved genes encoding small proteins in each lineage. We demonstrate that the symbiotic response was mostly in place in the most recent common ancestor of the RNS-forming species more than 90 million years ago.

Cell-type-specific alternative splicing in the Arabidopsis germline

During sexual reproduction in flowering plants, the two haploid sperm cells (SCs) embedded within the cytoplasm of a growing pollen tube are carried to the embryo sac for double fertilization. Pollen development in flowering plants is a dynamic process that encompasses changes at transcriptome and epigenome levels. While the transcriptome of pollen and SCs in Arabidopsis (Arabidopsis thaliana) is well documented, previous analyses have mostly been based on gene-level expression. In-depth transcriptome analysis, particularly the extent of alternative splicing (AS) at the resolution of SC and vegetative nucleus (VN), is still lacking. Therefore, we performed RNA-seq analysis to generate a spliceome map of Arabidopsis SCs and VN isolated from mature pollen grains. Based on our de novo transcriptome assembly, we identified 58,039 transcripts, including 9,681 novel transcripts, of which 2,091 were expressed in SCs and 3,600 in VN. Four hundred and sixty-eight genes were regulated both at gene and splicing levels, with many having functions in mRNA splicing, chromatin modification, and protein localization. Moreover, a comparison with egg cell RNA-seq data uncovered sex-specific regulation of transcription and splicing factors. Our study provides insights into a gamete-specific AS landscape at unprecedented resolution.

Thermosensitivity of pollen: a molecular perspective

A current trend in climate comprises adverse weather anomalies with more frequent and intense temperature events. Heatwaves are a serious threat to global food security because of the susceptibility of crop plants to high temperatures. Among various developmental stages of plants, even a slight rise in temperature during reproductive development proves detrimental, thus making sexual reproduction heat vulnerable. In this context, male gametophyte or pollen development stages are the most sensitive ones. High-temperature exposure induces pollen abortion, reducing pollen viability and germination rate with a concomitant effect on seed yield. This review summarizes the ultrastructural, morphological, biochemical, and molecular changes underpinning high temperature-induced aberrations in male gametophytes. Specifically, we highlight the temperature sensing cascade operating in pollen, involving reactive oxygen species (ROS), heat shock factors (HSFs), a hormones and transcriptional regulatory network. We also emphasize integrating various omics approaches to decipher the molecular events triggered by heat stress in pollen. The knowledge of genes, proteins, and metabolites conferring thermotolerance in reproductive tissues can be utilized to breed/engineer thermotolerant crops to ensure food security.

50 years of sperm cell isolations: from structural to omic studies

The fusion of male and female gametes is a fundamental process in the perpetuation and diversification of species. During the last 50 years, significant efforts have been made to isolate and characterize sperm cells from flowering plants, and to identify how these cells interact with female gametes to achieve double fertilization. The first techniques and analytical approaches not only provided structural and biochemical characterizations of plant sperm cells but also paved the way for in vitro fertilization studies. Further technological advances then led to unique insights into sperm biology at transcriptomic, proteomic and epigenetic level. Starting with a historical overview of sperm cell isolation techniques, we provide examples of how these contributed to create our current knowledge of sperm cell biology, and point out remaining challenges.

RUVBL proteins are involved in plant gametophyte development

The proper development of male and female gametophytes is critical for successful sexual reproduction and requires a carefully regulated series of events orchestrated by a suite of various proteins. RUVBL1 and RUVBL2, plant orthologues of human Pontin and Reptin, respectively, belong to the evolutionarily highly conserved AAA⁺ family linked to a wide range of cellular processes. Previously, we found that RUVBL1 and RUVBL2A mutations are homozygous lethal in Arabidopsis. Here, we report that RUVBL1 and RUVBL2A play roles in reproductive development. We show that mutant plants produce embryo sacs with an abnormal structure or with various numbers of nuclei. Although pollen grains of heterozygous mutant plants exhibit reduced viability and reduced pollen tube growth in vitro, some of the ruvbl pollen tubes are capable of targeting ovules in vivo. Similarly, some ruvbl ovules retain the ability to attract wild-type pollen tubes but fail to develop further. The activity of the RUVBL1 and RUVBL2A promoters was observed in the embryo sac, pollen grains, and tapetum cells and, for RUVBL2A, also in developing ovules. In summary, we show that the RUVBL proteins are essential for the proper development of both male and particularly female gametophytes in Arabidopsis.

I Choose You: Selecting Accurate Reference Genes for qPCR Expression Analysis in Reproductive Tissues in Arabidopsis thaliana

Quantitative real-time polymerase chain reaction (qPCR) is a widely used method to analyse the gene expression pattern in the reproductive tissues along with detecting gene levels in mutant backgrounds. This technique requires stable reference genes to normalise the expression level of target genes. Nonetheless, a considerable number of publications continue to present qPCR results normalised to a single reference gene and, to our knowledge, no comparative evaluation of multiple reference genes has been carried out in specific reproductive tissues of Arabidopsis thaliana. Herein, we assessed the expression stability levels of ten candidate reference genes (UBC9, ACT7, GAPC-2, RCE1, PP2AA3, TUA2, SAC52, YLS8, SAMDC and HIS3.3) in two conditional sets: one across flower development and the other using inflorescences from different genotypes. The stability analysis was performed using the RefFinder tool, which combines four statistical algorithms (geNorm, NormFinder, BestKeeper and the comparative ΔCt method). Our results showed that RCE1, SAC52 and TUA2 had the most stable expression in different flower developmental stages while YLS8, HIS3.3 and ACT7 were the top-ranking reference genes for normalisation in mutant studies. Furthermore, we validated our results by analysing the expression pattern of genes involved in reproduction and examining the expression of these genes in published mutant backgrounds. Overall, we provided a pool of appropriate reference genes for expression studies in reproductive tissues of A. thaliana, which will facilitate further gene expression studies in this context. More importantly, we presented a framework that will promote a consistent and accurate analysis of gene expression in any scientific field. Simultaneously, we highlighted the relevance of clearly defining and describing the experimental conditions associated with qPCR to improve scientific reproducibility.

Cross-stress gene expression atlas of Marchantia polymorpha reveals the hierarchy and regulatory principles of abiotic stress responses

Abiotic stresses negatively impact ecosystems and the yield of crops, and climate change will increase their frequency and intensity. Despite progress in understanding how plants respond to individual stresses, our knowledge of plant acclimatization to combined stresses typically occurring in nature is still lacking. Here, we used a plant with minimal regulatory network redundancy, Marchantia polymorpha, to study how seven abiotic stresses, alone and in 19 pairwise combinations, affect the phenotype, gene expression, and activity of cellular pathways. While the transcriptomic responses show a conserved differential gene expression between Arabidopsis and Marchantia, we also observe a strong functional and transcriptional divergence between the two species. The reconstructed high-confidence gene regulatory network demonstrates that the response to specific stresses dominates those of others by relying on a large ensemble of transcription factors. We also show that a regression model could accurately predict the gene expression under combined stresses, indicating that Marchantia performs arithmetic multiplication to respond to multiple stresses. Lastly, two online resources ( https://conekt.plant.tools and http://bar.utoronto.ca/efp_marchantia/cgi-bin/efpWeb.cgi ) are provided to facilitate the study of gene expression in Marchantia exposed to abiotic stresses.

Toward kingdom-wide analyses of gene expression

Gene expression data for Archaeplastida are accumulating exponentially, with more than 300 000 RNA-sequencing (RNA-seq) experiments available for hundreds of species. The gene expression data stem from thousands of experiments that capture gene expression in various organs, tissues, cell types, (a)biotic perturbations, and genotypes. Advances in software tools make it possible to process all these data in a matter of weeks on modern office computers, giving us the possibility to study gene expression in a kingdom-wide manner for the first time. We discuss how the expression data can be accessed and processed and outline analyses that take advantage of cross-species analyses, allowing us to generate powerful and robust hypotheses about gene function and evolution.

Tomato POLLEN DEFICIENT 2 encodes a G-type lectin receptor kinase required for viable pollen grain formation

Pollen development is a crucial biological process indispensable for seed set in flowering plants and for successful crop breeding. However, little is known about the molecular mechanisms regulating pollen development in crop species. This study reports a novel male-sterile tomato mutant, pollen deficient 2 (pod2), characterized by the production of non-viable pollen grains and resulting in the development of small parthenocarpic fruits. A combined strategy of mapping-by-sequencing and RNA interference-mediated gene silencing was used to prove that the pod2 phenotype is caused by the loss of Solanum lycopersicum G-type lectin receptor kinase II.9 (SlG-LecRK-II.9) activity. In situ hybridization of floral buds showed that POD2/SlG-LecRK-II.9 is specifically expressed in tapetal cells and microspores at the late tetrad stage. Accordingly, abnormalities in meiosis and tapetum programmed cell death in pod2 occurred during microsporogenesis, resulting in the formation of four dysfunctional microspores leading to an aberrant microgametogenesis process. RNA-seq analyses supported the existence of alterations at the final stage of microsporogenesis, since we found tomato deregulated genes whose counterparts in Arabidopsis are essential for the normal progression of male meiosis and cytokinesis. Collectively, our results revealed the essential role of POD2/SlG-LecRK-II.9 in regulating tomato pollen development.

How plants conquered land: evolution of terrestrial adaptation

The transition of plants from water to land is considered one of the most significant events in the evolution of life on Earth. The colonization of land by plants, accompanied by their morphological, physiological and developmental changes, resulted in plant biodiversity. Besides significantly influencing oxygen levels in the air and on land, plants manufacture organic matter from CO₂ and water with the help of sunlight, paving the way for the diversification of nonplant lineages ranging from microscopic organisms to animals. Land plants regulate the climate by adjusting total biomass and energy flow. At the genetic level, these innovations are achieved through the rearrangement of pre-existing genetic information. Advances in genome sequencing technology are revamping our understanding of plant evolution. This study highlights the morphological and genomic innovations that allow plants to integrate life on Earth.

RNA-seq analysis of synchronized developing pollen isolated from a single anther

Pollen development, from unicellular microspores to anthesis, is a complex process involving the coordinated specification, differentiation and functions of different cell types. Key to understanding this development is identifying the genes expressed at precise stages of development. However, transcriptomic studies on pollen prior to anthesis are complicated by the inaccessible nature of pollen developing in the anther and the resistant pollen wall. To assist with understanding gene expression during pollen development we have developed a protocol to perform RNA-Seq on pollen isolated from a single anther (SA RNA-Seq). The protocol involves removing pollen from a single anther for analysis and viewing the remaining pollen to determine the developmental stage. The isolated pollen is chemically lysed and mRNA isolated from the lysate using an oligo-dT column before library preparation. Here, we report on the development and testing of our method and the generation of a transcriptome for three stages of pollen development from Arabidopsis (Arabidopsis thaliana) and two stages from male kiwifruit (Actinidia chinensis). This protocol enables the transcriptome of precise developmental stages of pollen to be analyzed, and uses a small number of plants, potentially facilitating studies that require a range of treatments or the analysis of the first generation of transgenic plants.

Evolutionary analysis of the LORELEI gene family in plants reveals regulatory subfunctionalization

A signaling complex comprising members of the LORELEI (LRE)-LIKE GPI-anchored protein (LLG) and Catharanthus roseus RECEPTOR-LIKE KINASE 1-LIKE (CrRLK1L) families perceive RAPID ALKALINIZATION FACTOR (RALF) peptides and regulate growth, reproduction, immunity, and stress responses in Arabidopsis (Arabidopsis thaliana). Genes encoding these proteins are members of multigene families in most angiosperms and could generate thousands of signaling complex variants. However, the links between expansion of these gene families and the functional diversification of this critical signaling complex as well as the evolutionary factors underlying the maintenance of gene duplicates remain unknown. Here, we investigated LLG gene family evolution by sampling land plant genomes and explored the function and expression of angiosperm LLGs. We found that LLG diversity within major land plant lineages is primarily due to lineage-specific duplication events, and that these duplications occurred both early in the history of these lineages and more recently. Our complementation and expression analyses showed that expression divergence (i.e. regulatory subfunctionalization), rather than functional divergence, explains the retention of LLG paralogs. Interestingly, all but one monocot and all eudicot species examined had an LLG copy with preferential expression in male reproductive tissues, while the other duplicate copies showed highest levels of expression in female or vegetative tissues. The single LLG copy in Amborella trichopoda is expressed vastly higher in male compared to in female reproductive or vegetative tissues. We propose that expression divergence plays an important role in retention of LLG duplicates in angiosperms.

MarpolBase Expression: A Web-Based, Comprehensive Platform for Visualization and Analysis of Transcriptomes in the Liverwort Marchantia polymorpha

The liverwort Marchantia polymorpha is equipped with a wide range of molecular and genetic tools and resources that have led to its wide use to explore the evo-devo aspects of land plants. Although its diverse transcriptome data are rapidly accumulating, there is no extensive yet user-friendly tool to exploit such a compilation of data and to summarize results with the latest annotations. Here, we have developed a web-based suite of tools, MarpolBase Expression (MBEX, https://marchantia.info/mbex/), where users can visualize gene expression profiles, identify differentially expressed genes, perform co-expression and functional enrichment analyses and summarize their comprehensive output in various portable formats. Using oil body biogenesis as an example, we demonstrated that the results generated by MBEX were consistent with the published experimental evidence and also revealed a novel transcriptional network in this process. MBEX should facilitate the exploration and discovery of the genetic and functional networks behind various biological processes in M. polymorpha and promote our understanding of the evolution of land plants.

Regulatory dynamics of gene expression in the developing male gametophyte of Arabidopsis

Sexual reproduction in angiosperms requires the production and delivery of two male gametes by a three-celled haploid male gametophyte. This demands synchronized gene expression in a short developmental window to ensure double fertilization and seed set. While transcriptomic changes in developing pollen are known for Arabidopsis, no studies have integrated RNA and proteomic data in this model. Further, the role of alternative splicing has not been fully addressed, yet post-transcriptional and post-translational regulation may have a key role in gene expression dynamics during microgametogenesis. We have refined and substantially updated global transcriptomic and proteomic changes in developing pollen for two Arabidopsis accessions. Despite the superiority of RNA-seq over microarray-based platforms, we demonstrate high reproducibility and comparability. We identify thousands of long non-coding RNAs as potential regulators of pollen development, hundreds of changes in alternative splicing and provide insight into mRNA translation rate and storage in developing pollen. Our analysis delivers an integrated perspective of gene expression dynamics in developing Arabidopsis pollen and a foundation for studying the role of alternative splicing in this model.

A genome resource for Acacia, Australia's largest plant genus

Acacia (Leguminosae, Caesalpinioideae, mimosoid clade) is the largest and most widespread genus of plants in the Australian flora, occupying and dominating a diverse range of environments, with an equally diverse range of forms. For a genus of its size and importance, Acacia currently has surprisingly few genomic resources. Acacia pycnantha, the golden wattle, is a woody shrub or tree occurring in south-eastern Australia and is the country's floral emblem. To assemble a genome for A. pycnantha, we generated long-read sequences using Oxford Nanopore Technology, 10x Genomics Chromium linked reads, and short-read Illumina sequences, and produced an assembly spanning 814 Mb, with a scaffold N50 of 2.8 Mb, and 98.3% of complete Embryophyta BUSCOs. Genome annotation predicted 47,624 protein-coding genes, with 62.3% of the genome predicted to comprise transposable elements. Evolutionary analyses indicated a shared genome duplication event in the Caesalpinioideae, and conflict in the relationships between Cercis (subfamily Cercidoideae) and subfamilies Caesalpinioideae and Papilionoideae (pea-flowered legumes). Comparative genomics identified a suite of expanded and contracted gene families in A. pycnantha, and these were annotated with both GO terms and KEGG functional categories. One expanded gene family of particular interest is involved in flowering time and may be associated with the characteristic synchronous flowering of Acacia. This genome assembly and annotation will be a valuable resource for all studies involving Acacia, including the evolution, conservation, breeding, invasiveness, and physiology of the genus, and for comparative studies of legumes.

Diversification of heat shock transcription factors expanded thermal stress responses during early plant evolution

The copy numbers of many plant transcription factor (TF) genes substantially increased during terrestrialization. This allowed TFs to acquire new specificities and thus create gene regulatory networks (GRNs) with new biological functions to help plants adapt to terrestrial environments. Through characterizing heat shock factor (HSF) genes MpHSFA1 and MpHSFB1 in the liverwort Marchantia polymorpha, we explored how heat-responsive GRNs widened their functions in M. polymorpha and Arabidopsis thaliana. An interspecies comparison of heat-induced transcriptomes and the evolutionary rates of HSFs demonstrated the emergence and subsequent rapid evolution of HSFB prior to terrestrialization. Transcriptome and metabolome analyses of M. polymorpha HSF-null mutants revealed that MpHSFA1 controls canonical heat responses such as thermotolerance and metabolic changes. MpHSFB1 also plays essential roles in heat responses, as well as regulating developmental processes including meristem branching and antheridiophore formation. Analysis of cis-regulatory elements revealed development- and stress-related TFs that function directly or indirectly downstream of HSFB. Male gametophytes of M. polymorpha showed higher levels of thermotolerance than female gametophytes, which could be explained by different expression levels of MpHSFA1U and MpHSFA1V on sex chromosome. We propose that the diversification of HSFs is linked to the expansion of HS responses, which enabled coordinated multicellular reactions in land plants.

Feature importance network reveals novel functional relationships between biological features in Arabidopsis thaliana

Understanding how the different cellular components are working together to form a living cell requires multidisciplinary approaches combining molecular and computational biology. Machine learning shows great potential in life sciences, as it can find novel relationships between biological features. Here, we constructed a dataset of 11,801 gene features for 31,522 Arabidopsis thaliana genes and developed a machine learning workflow to identify linked features. The detected linked features are visualised as a Feature Important Network (FIN), which can be mined to reveal a variety of novel biological insights pertaining to gene function. We demonstrate how FIN can be used to generate novel insights into gene function. To make this network easily accessible to the scientific community, we present the FINder database, available at finder.plant.tools.

High-energy-level metabolism and transport occur at the transition from closed to open flowers

During the maturation phase of flower development, the onset of anthesis visibly marks the transition from buds to open flowers, during which petals stretch out, nectar secretion commences, and pollination occurs. Analysis of the metabolic changes occurring during this developmental transition has primarily focused on specific classes of metabolites, such as pigments and scent emission, and far less on the whole network of primary and secondary metabolites. To investigate the metabolic changes occurring at anthesis, we performed multi-platform metabolomics alongside RNA sequencing in individual florets harvested from the main inflorescence of Arabidopsis (Arabidopsis thaliana) ecotype Col-0. To trace metabolic fluxes at the level of the whole inflorescence and individual florets, we further integrated these studies with radiolabeled experiments. These extensive analyses revealed high-energy-level metabolism and transport of carbohydrates and amino acids, supporting intense metabolic rearrangements occurring at the time of this floral transition. These comprehensive data are discussed in the context of our current understanding of the metabolic shifts underlying flower opening. We envision that this analysis will facilitate the introgression of floral metabolic traits promoting pollination in crop species for which a comprehensive knowledge of flower metabolism is still limited.

An overview of bioinformatics, genomics, and transcriptomics resources for bryophytes

Bryophytes are useful models for the study of plant evolution, development, plant-fungal symbiosis, stress responses, and gametogenesis. Additionally, their dominant haploid gametophytic phase makes them great models for functional genomics research, allowing straightforward genome editing and gene knockout via CRISPR or homologous recombination. Until 2016, however, the only bryophyte genome sequence published was that of Physcomitrium patens. Throughout recent years, several other bryophyte genomes and transcriptome datasets became available, enabling better comparative genomics in evolutionary studies. The increase in the number of bryophyte genome and transcriptome resources available has yielded a plethora of annotations, databases, and bioinformatics tools to access the new data, which covers the large diversity of this clade and whose biology comprises features such as association with arbuscular mycorrhiza fungi, sex chromosomes, low gene redundancy, or loss of RNA editing genes for organellar transcripts. Here we provide a guide to resources available for bryophytes with regards to genome and transcriptome databases and bioinformatics tools.

Exploiting plant transcriptomic databases: Resources, tools, and approaches

There are now more than 300 000 RNA sequencing samples available, stemming from thousands of experiments capturing gene expression in organs, tissues, developmental stages, and experimental treatments for hundreds of plant species. The expression data have great value, as they can be re-analyzed by others to ask and answer questions that go beyond the aims of the study that generated the data. Because gene expression provides essential clues to where and when a gene is active, the data provide powerful tools for predicting gene function, and comparative analyses allow us to study plant evolution from a new perspective. This review describes how we can gain new knowledge from gene expression profiles, expression specificities, co-expression networks, differential gene expression, and experiment correlation. We also introduce and demonstrate databases that provide user-friendly access to these tools.

The adaptive evolution of Euryale ferox to the aquatic environment through paleo-hexaploidization

Occupation of living space is one of the main driving forces of adaptive evolution, especially for aquatic plants whose leaves float on the water surface and thus have limited living space. Euryale ferox, from the angiosperm basal family Nymphaeaceae, develops large, rapidly expanding leaves to compete for space on the water surface. Microscopic observation found that the cell proliferation of leaves is almost completed underwater, while the cell expansion occurs rapidly after they grow above water. To explore the mechanism underlying the specific development of leaves, we performed sequences assembly and analyzed the genome and transcriptome dynamics of E. ferox. Through reconstruction of the three sub-genomes generated from the paleo-hexaploidization event in E. ferox, we revealed that one sub-genome was phylogenetically closer to Victoria cruziana, which also exhibits gigantic floating leaves. Further analysis revealed that while all three sub-genomes promoted the evolution of the specific leaf development in E. ferox, the genes from the sub-genome closer to V. cruziana contributed more to this adaptive evolution. Moreover, we found that genes involved in cell proliferation and expansion, photosynthesis, and energy transportation were over-retained and showed strong expression association with the leaf development stages, such as the expression divergence of SWEET orthologs as energy uploaders and unloaders in the sink and source leaf organs of E. ferox. These findings provide novel insights into the genome evolution through polyploidization, as well as the adaptive evolution regarding the leaf development accomplished through biased gene retention and expression sub-functionalization of multi-copy genes in E. ferox.

A Plastid-Bound Ankyrin Repeat Protein Controls Gametophyte and Early Embryo Development in Arabidopsis thaliana

Proplastids are essential precursors for multi-fate plastid biogenesis, including chloroplast differentiation, a powerhouse for photosynthesis in plants. Arabidopsis ankyrin repeat protein (AKRP, AT5G66055) is a plastid-localized protein with a putative function in plastid differentiation and morphogenesis. Loss of function of akrp leads to embryo developmental arrest. Whether AKRP is critical pre-fertilization has remained unresolved. Here, using reverse genetics, we report a new allele, akrp-3, that exhibited a reduced frequency of mutant embryos (<13%) compared to previously reported alleles. akrp-3 affected both male and female gametophytes resulting in reduced viability, incompetence in pollen tube attraction, altered gametic cell fate, and embryo arrest that were depleted of chlorophyll. AKRP is widely expressed, and the AKRP-GFP fusion localized to plastids of both gametophytes, in isolated chloroplast and co-localized with a plastid marker in pollen and pollen tubes. Cell-type-specific complementation of akrp-3 hinted at the developmental timing at which AKRP might play an essential role. Our findings provide a plausible insight into the crucial role of AKRP in the differentiation of both gametophytes and coupling embryo development with chlorophyll synthesis.

The Physcomitrium patens egg cell expresses several distinct epigenetic components and utilizes homologues of BONOBO genes for cell specification

Although land plant germ cells have received much attention, knowledge about their specification is still limited. We thus identified transcripts enriched in egg cells of the bryophyte model species Physcomitrium patens, compared the results with angiosperm egg cells, and selected important candidate genes for functional analysis. We used laser-assisted microdissection to perform a cell-type-specific transcriptome analysis on egg cells for comparison with available expression profiles of vegetative tissues and male reproductive organs. We made reporter lines and knockout mutants of the two BONOBO (PbBNB) genes and studied their role in reproduction. We observed an overlap in gene activity between bryophyte and angiosperm egg cells, but also clear differences. Strikingly, several processes that are male-germline specific in Arabidopsis are active in the P. patens egg cell. Among those were the moss PbBNB genes, which control proliferation and identity of both female and male germlines. Pathways shared between male and female germlines were most likely present in the common ancestors of land plants, besides sex-specifying factors. A set of genes may also be involved in the switches between the diploid and haploid moss generations. Nonangiosperm gene networks also contribute to the specification of the P. patens egg cell.

Gynoecium and fruit development in Arabidopsis

Flowering plants produce flowers and one of the most complex floral structures is the pistil or the gynoecium. All the floral organs differentiate from the floral meristem. Various reviews exist on molecular mechanisms controlling reproductive development, but most focus on a short time window and there has been no recent review on the complete developmental time frame of gynoecium and fruit formation. Here, we highlight recent discoveries, including the players, interactions and mechanisms that govern gynoecium and fruit development in Arabidopsis. We also present the currently known gene regulatory networks from gynoecium initiation until fruit maturation.

Nanotechnology Strategies for Plant Genetic Engineering

Plant genetic engineering is essential for improving crop yield, quality, and resistance to abiotic/biotic stresses for sustainable agriculture. Agrobacterium-, biolistic bombardment-, electroporation-, and poly(ethylene glycol) (PEG)-mediated genetic-transformation systems are extensively used in plant genetic engineering. However, these systems have limitations, including species dependency, destruction of plant tissues, low transformation efficiency, and high cost. Recently, nanotechnology-based gene-delivery methods have been developed for plant genetic transformation. This nanostrategy shows excellent transformation efficiency, good biocompatibility, adequate protection of exogenous nucleic acids, and the potential for plant regeneration. However, the nanomaterial-mediated gene-delivery system in plants is still in its infancy, and there are many challenges for its broad applications. Herein, the conventional genetic transformation techniques used in plants are briefly discussed. After that, the progress in the development of nanomaterial-based gene-delivery systems is considered. CRISPR-Cas-mediated genome editing and its combined applications with plant nanotechnology are also discussed. The conceptual innovations, methods, and practical applications of nanomaterial-mediated genetic transformation summarized herein will be beneficial for promoting plant genetic engineering in modern agriculture.

Multi-Omics-Based Discovery of Plant Signaling Molecules

Plants produce numerous structurally and functionally diverse signaling metabolites, yet only relatively small fractions of which have been discovered. Multi-omics has greatly expedited the discovery as evidenced by increasing recent works reporting new plant signaling molecules and relevant functions via integrated multi-omics techniques. The effective application of multi-omics tools is the key to uncovering unknown plant signaling molecules. This review covers the features of multi-omics in the context of plant signaling metabolite discovery, highlighting how multi-omics addresses relevant aspects of the challenges as follows: (a) unknown functions of known metabolites; (b) unknown metabolites with known functions; (c) unknown metabolites and unknown functions. Based on the problem-oriented overview of the theoretical and application aspects of multi-omics, current limitations and future development of multi-omics in discovering plant signaling metabolites are also discussed.

Genomic Signatures of Sexual Selection on Pollen-Expressed Genes in Arabis alpina

Fertilization in angiosperms involves the germination of pollen on the stigma, followed by the extrusion of a pollen tube that elongates through the style and delivers two sperm cells to the embryo sac. Sexual selection could occur throughout this process when male gametophytes compete for fertilization. The strength of sexual selection during pollen competition should be affected by the number of genotypes deposited on the stigma. As increased self-fertilization reduces the number of mating partners, and the genetic diversity and heterozygosity of populations, it should thereby reduce the intensity of sexual selection during pollen competition. Despite the prevalence of mating system shifts, few studies have directly compared the molecular signatures of sexual selection during pollen competition in populations with different mating systems. Here we analyzed whole-genome sequences from natural populations of Arabis alpina, a species showing mating system variation across its distribution, to test whether shifts from cross- to self-fertilization result in molecular signatures consistent with sexual selection on genes involved in pollen competition. We found evidence for efficient purifying selection on genes expressed in vegetative pollen, and overall weaker selection on sperm-expressed genes. This pattern was robust when controlling for gene expression level and specificity. In agreement with the expectation that sexual selection intensifies under cross-fertilization, we found that the efficacy of purifying selection on male gametophyte-expressed genes was significantly stronger in genetically more diverse and outbred populations. Our results show that intra-sexual competition shapes the evolution of pollen-expressed genes, and that its strength fades with increasing self-fertilization rates.

The 3D architecture and molecular foundations of de novo centriole assembly via bicentrioles

Centrioles are structurally conserved organelles, composing both centrosomes and cilia. In animal cycling cells, centrioles often form through a highly characterized process termed canonical duplication. However, a large diversity of eukaryotes assemble centrioles de novo through uncharacterized pathways. This unexplored diversity is key to understanding centriole assembly mechanisms and how they evolved to assist specific cellular functions. Here, we show that, during spermatogenesis of the bryophyte Physcomitrium patens, centrioles are born as a co-axially oriented centriole pair united by a cartwheel. Interestingly, we observe that these centrioles are twisted in opposite orientations. Microtubules emanate from the bicentrioles, which localize to the spindle poles during cell division. After their separation, the two resulting sister centrioles mature asymmetrically, elongating specific microtubule triplets and a naked cartwheel. Subsequently, two motile cilia are assembled that appear to alternate between different motility patterns. We further show that centriolar components SAS6, Bld10, and POC1, which are conserved across eukaryotes, are expressed during spermatogenesis and required for this de novo biogenesis pathway. Our work supports a scenario where centriole biogenesis, while driven by conserved molecular modules, is more diverse than previously thought.