Expression Atlas of Selaginella moellendorffii Provides Insights into the Evolution of Vasculature, Secondary Metabolism, and Roots

Structural and biochemical characterization of SmoPG1, an exo-polygalacturonase from Selaginella moellendorffii

Polygalacturonases (PGs) can modulate chemistry and mechanical properties of the plant cell wall through the degradation of pectins, one of its major constituents. PGs are largely used in food, beverage, textile, and paper industries to increase processes' performances. To improve the use of PGs, knowledge of their biochemical, structural and functional features is of prime importance. Our study aims at characterizing SmoPG1, a polygalacturonase from Selaginella moellendorffii, that belongs to the lycophytes. Transcription data showed that SmoPG1 was mainly expressed in S. moellendorffii shoots while phylogenetic analyses suggested that SmoPG1 is an exo-PG, which was confirmed by the biochemical characterization following its expression in heterologous system. Indeed, LC-MS/MS oligoprofiling using various pectic substrates identified galacturonic acid (GalA) as the main hydrolysis product. We found that SmoPG1 was most active on polygalacturonic acid (PGA) at pH 5, and that its activity could be modulated by different cations (Ca2+, Cu2+, Fe2+, Mg2+, Mn2+, Na2+, Zn2+). In addition, SmoPG1 was inhibited by green tea catechins, including (-)-epigallocatechin-3-gallate (EGCG). Docking analyses and MD simulations showed in detail amino acids responsible for the SmoPG1-EGCG interaction. Considering its expression yield and activity, SmoPG1 appears as a prime candidate for the industrial production of GalA.

The Shoot Apical Meristem: An Evolutionary Molding of Higher Plants

The shoot apical meristem (SAM) gives rise to the aerial structure of plants by producing lateral organs and other meristems. The SAM is responsible for plant developmental patterns, thus determining plant morphology and, consequently, many agronomic traits such as the number and size of fruits and flowers and kernel yield. Our current understanding of SAM morphology and regulation is based on studies conducted mainly on some angiosperms, including economically important crops such as maize (Zea mays) and rice (Oryza sativa), and the model species Arabidopsis (Arabidopsis thaliana). However, studies in other plant species from the gymnosperms are scant, making difficult comparative analyses that help us understand SAM regulation in diverse plant species. This limitation prevents deciphering the mechanisms by which evolution gave rise to the multiple plant structures within the plant kingdom and determines the conserved mechanisms involved in SAM maintenance and operation. This review aims to integrate and analyze the current knowledge of SAM evolution by combining the morphological and molecular information recently reported from the plant kingdom.

Control of stem cell behavior by CLE-JINGASA signaling in the shoot apical meristem in Marchantia polymorpha

Land plants undergo indeterminate growth by the activity of meristems in both gametophyte (haploid) and sporophyte (diploid) generations. In the sporophyte of the flowering plant Arabidopsis thaliana, the apical meristems are located at the shoot and root tips in which a number of regulatory gene homologs are shared for their development, implying deep evolutionary origins. However, little is known about their functional conservation with gametophytic meristems in distantly related land plants such as bryophytes, even though genomic studies have revealed that the subfamily-level diversity of regulatory genes is mostly conserved throughout land plants. Here, we show that a NAM/ATAF/CUC (NAC) domain transcription factor, JINGASA (MpJIN), acts downstream of CLAVATA3 (CLV3)/ESR-related (CLE) peptide signaling and controls stem cell behavior in the gametophytic shoot apical meristem of the liverwort Marchantia polymorpha. In the meristem, strong MpJIN expression was associated with the periclinal cell division at the periphery of the stem cell zone (SCZ), whereas faint MpJIN expression was found at the center of the SCZ. Time course observation indicates that the MpJIN-negative cells are lost from the SCZ and respecified de novo at two separate positions during the dichotomous branching event. Consistently, the induction of MpJIN results in ectopic periclinal cell division in the SCZ and meristem termination. Based on the comparative expression data, we speculate that the function of JIN/FEZ subfamily genes was shared among the shoot apical meristems in the gametophyte and sporophyte generations in early land plants but was lost in certain lineages, including the flowering plant A. thaliana.

Genome and transcriptome of Selaginella kraussiana reveal evolution of root apical meristems in vascular plants

The evolution of roots allowed vascular plants to adapt to land environments. Fossil evidence indicates that roots evolved independently in euphyllophytes (ferns and seed plants) and lycophytes, the two lineages of extant vascular plants. Based on a high-quality genome assembly, mRNA sequencing (mRNA-seq) data, and single-cell RNA-seq data for the lycophyte Selaginella kraussiana, we show that the two root origin events in lycophytes and euphyllophytes adopted partially similar molecular modules in the regulation of root apical meristem (RAM) development. In S. kraussiana, the RAM initiates from the rhizophore primordium guided by auxin and duplicates itself by dichotomous branching. The auxin signaling pathway directly upregulates euAINTEGUMENTAb (SkeuANTb), and then SkeuANTb directly promotes the expression of SkeuANTa and the WUSCHEL-RELATED HOMEOBOX13b (SkWOX13b) for RAM maintenance, partially similar to the molecular pathway involving the euANT-branch PLETHORA (AtPLT) genes and AtWOX5 in root initiation in the seed plant Arabidopsis thaliana. Other molecular modules, e.g., SHORT-ROOT and SCARECROW, also have partially similar expression patterns in the RAMs of S. kraussiana and A. thaliana. Overall, our study not only provides genome and transcriptome tools of S. kraussiana but also indicates the employment of some common molecular modules in RAMs during root origins in lycophytes and euphyllophytes.

Diverse branching forms regulated by a core auxin transport mechanism in plants

Diverse branching forms have evolved multiple times across the tree of life to facilitate resource acquisition and exchange with the environment. In the vascular plant group, the ancestral pattern of branching involves dichotomy of a parent shoot apex to form two new daughter apices. The molecular basis of axillary branching in Arabidopsis is well understood, but few regulators of dichotomous branching are known. Through analyses of dichotomous branching in the lycophyte, Selaginella kraussiana, we identify PIN-mediated auxin transport as an ancestral branch regulator of vascular plants. We show that short-range auxin transport out of the apices promotes dichotomy and that branch dominance is globally coordinated by long-range auxin transport. Uniquely in Selaginella, angle meristems initiate at each dichotomy, and these can develop into rhizophores or branching angle shoots. We show that long-range auxin transport and a transitory drop in PIN expression are involved in angle shoot development. We conclude that PIN-mediated auxin transport is an ancestral mechanism for vascular plant branching that was independently recruited into Selaginella angle shoot development and seed plant axillary branching during evolution.

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.

Cellular and gene expression patterns associated with root bifurcation in Selaginella

The roots of lycophytes branch through dichotomy or bifurcation, during which the root apex splits into two daughter roots. This is morphologically distinct from lateral root (LR) branching in the extant euphyllophytes, with LRs developing along the root axis at different distances from the apex. Although the process of root bifurcation is poorly understood, such knowledge can be important, because it may represent an evolutionarily ancient strategy that roots recruited to form new stem cells or meristems. In this study, we examined root bifurcation in the lycophyte Selaginella moellendorffii. We characterized an in vitro developmental time frame based on repetitive apex bifurcations, allowing us to sample different stages of dichotomous root branching and analyze the root meristem and root branching in S. moellendorffii at the microscopic and transcriptomic level. Our results showed that, in contrast to previous assumptions, initial cells (ICs) in the root meristem are mostly not tetrahedral but rather show an irregular shape. Tracking down the early stages of root branching argues for the occurrence of a symmetric division of the single IC, resulting in two apical stem cells that initiate root meristem bifurcation. Moreover, we generated a S. moellendorffii root branching transcriptome that resulted in the delineation of a subset of core meristem genes. The occurrence of multiple putative orthologs of meristem genes in this dataset suggests the presence of conserved pathways in the control of meristem and root stem cell establishment or maintenance.

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.

Transcriptional analysis of Ceratopteris richardii young sporophyte reveals conservation of stem cell factors in the root apical meristem

Gene expression in roots has been assessed in different plant species in studies ranging from complete organs to specific cell layers, and more recently at the single cell level. While certain genes or functional categories are expressed in the root of all or most plant species, lineage-specific genes have also been discovered. An increasing amount of transcriptomic data is available for angiosperms, while a limited amount of data is available for ferns, and few studies have focused on fern roots. Here, we present a de novo transcriptome assembly from three different parts of the Ceratopteris richardii young sporophyte. Differential gene expression analysis of the root tip transcriptional program showed an enrichment of functional categories related to histogenesis and cell division, indicating an active apical meristem. Analysis of a diverse set of orthologous genes revealed conserved expression in the root meristem, suggesting a preserved role for different developmental roles in this tissue, including stem cell maintenance. The reconstruction of evolutionary trajectories for ground tissue specification genes suggests a high degree of conservation in vascular plants, but not for genes involved in root cap development, showing that certain genes are absent in Ceratopteris or have intricate evolutionary paths difficult to track. Overall, our results suggest different processes of conservation and divergence of genes involved in root development.

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.

Comprehensive analysis of 1R- and 2R-MYBs reveals novel genic and protein features, complex organisation, selective expansion and insights into evolutionary tendencies

Myeloblastosis (MYB) family, the largest plant transcription factor family, has been subcategorised based on the number and type of repeats in the MYB domain. In spite of several reports, evolution of MYB genes and repeats remains enigmatic. Brassicaceae members are endowed with complex genomes, including dysploidy because of its unique history with multiple rounds of polyploidisation, genomic fractionations and rearrangements. The present study is an attempt to gain insights into the complexities of MYB family diversity, understand impacts of genome evolution on gene families and develop an evolutionary framework to understand the origin of various subcategories of MYB gene family. We identified and analysed 1129 MYBs that included 1R-, 2R-, 3R- and atypical-MYBs across sixteen species representing protists, fungi, animals and plants and exclude MYB identified from Brassicaceae except Arabidopsis thaliana; in addition, a total of 1137 2R-MYB genes from six Brassicaceae species were also analysed. Comparative analysis revealed predominance of 1R-MYBs in protists, fungi, animals and lower plants. Phylogenetic reconstruction and analysis of selection pressure suggested ancestral nature of R1-type repeat containing 1R-MYBs that might have undergone intragenic duplication to form multi-repeat MYBs. Distinct differences in gene structure between 1R-MYB and 2R-MYBs were observed regarding intron number, the ratio of gene length to coding DNA sequence (CDS) length and the length of exons encoding the MYB domain. Conserved as well as novel and lineage-specific intron phases were identified. Analyses of physicochemical properties revealed drastic differences indicating functional diversification in MYBs. Phylogenetic reconstruction of 1R- and 2R-MYB genes revealed a shared structure-function relationship in clades which was supported when transcriptome data was analysed in silico. Comparative genomics to study distribution pattern and mapping of 2R-MYBs revealed congruency and greater degree of synteny and collinearity among closely related species. Micro-synteny analysis of genomic segments revealed high conservation of genes that are immediately flanking the surrounding tandemly organised 2R-MYBs along with instances of local duplication, reorganisations and genome fractionation. In summary, polyploidy, dysploidy, reshuffling and genome fractionation were found to cause loss or gain of 2R-MYB genes. The findings need to be supported with functional validation to understand gene structure-function relationship along the evolutionary lineage and adaptive strategies based on comparative functional genomics in plants.

Identification of genes differentially expressed between prostrate shoots and erect shoots in the lycophyte Selaginella nipponica using an RNA-seq approach

Lycophytes are the earliest vascular plants and Selaginella is the most studied genus among them. Prostrate shoots are produced during early growth and erect shoots emerge later in S. nipponica, thus providing an opportunity for exploring the evolution of the mechanism underlying the transition between growth phases. Six libraries were sequenced for the prostrate and the erect shoots, and a total of 206 768 genes were identified. Some genes were differentially expressed in prostate and erect shoot, with relatively high expression in the prostate shoots being related to hormone responses and defence reactions, while higher expression in the erect shoots was related to spore formation and shoot development. Some SPL genes possessed a miR156 binding site and were highly expressed in the erect shoots, while AP2-like genes were more highly expressed in the prostrate shoots but simultaneously lacked any miR172 binding site. MiR156 was detected at a higher concentration in the prostrate shoots. Thus, the mechanism for the vegetative to reproductive transition of sporophytes probably originated in the common ancestor of vascular plants and must have experienced stepwise development during evolution.

Oxylipin biosynthesis in spikemoss Selaginella moellendorffii: Identification of allene oxide synthase (CYP74L2) and hydroperoxide lyase (CYP74L1)

Nonclassical P450s of the CYP74 family catalyse the secondary conversions of fatty acid hydroperoxides to bioactive oxylipins in plants. The model organism, spikemoss Selaginella moellendorffii Hieron, possesses at least ten CYP74 genes of novel J, K, L, and M subfamilies. The cloning of three CYP74L genes and catalytic properties of recombinant proteins are described in the present work. The CYP74L1 possessed mainly hydroperoxide lyase (HPL) activity towards the 13(S)-hydroperoxide of α-linolenic acids (13-HPOT) and nearly equal HPL and allene oxide synthase (AOS) activities towards the 13(S)-hydroperoxide of linoleic acids (13-HPOD). The 9-hydroperoxides were poor substrates for CYP74L1 and led to the production of mainly the α-ketols (AOS products) and minorities of HPL and epoxyalcohol synthase (EAS) products. The CYP74L2 possessed the AOS activity towards all tested hydroperoxides. CYP74L3 possessed low HPL/EAS activity. Besides, the aerial parts of S. moellendorffii plants possessed complex oxylipins patterns including divinyl ethers, epoxyalcohols, and 12-oxo-phytodienoic acid. Characterization of the CYP74L enzymes and oxylipin pattern updates the knowledge on the complex oxylipin biosynthetic machinery in the surviving oldest taxa of vascular plants.

Green Leaf Volatile-Burst in Selaginella moellendorffii

Green leaf volatiles (GLVs) consist of six-carbon volatile aldehydes, alcohols, and their esters. They are formed from polyunsaturated fatty acids and are involved in the defense of plants against herbivores and pathogens. GLVs generally have low concentrations in intact healthy plant tissues, but the biosynthetic pathway to form GLVs is quickly activated by mechanical damage to tissues, an event called the GLV-burst. Most seed plants have the ability to implement GLV-burst; however, this potential in non-seed plants has not been extensively researched. In this study, we examined the GLV-burst capacity of monilophytes, lycophytes, and bryophytes, and confirmed that monilophytes and lycophytes showed substantial GLV-burst ability, while bryophytes did not, with a few exceptions. When the genome sequence of a model lycophyte, Selaginella moellendorffii was reviewed, 10 genes were found that showed high similarity with the non-canonical cytochrome P450 enzymes, CYP74s, specialized in oxylipin formation. Recombinant proteins expressed with Escherichia coli showed that one of them had the ability to encode allene oxide synthase, and another encoded hydroperoxide lyase (HPL), preferring linolenic acid 13-hydroperoxide, and it was inferred that this gene was responsible for GLV-burst in S. moellendorffii. Based on the phylogenetic tree constructed with CYP74s of non-seed and seed plants, we hypothesized that HPL was acquired independently in the lycophyte and seed plants through diversification of CYP74 genes.

LSTrAP-Kingdom: an automated pipeline to generate annotated gene expression atlases for kingdoms of life

MOTIVATION

There are now more than two million RNA sequencing experiments for plants, animals, bacteria and fungi publicly available, allowing us to study gene expression within and across species and kingdoms. However, the tools allowing the download, quality control and annotation of this data for more than one species at a time are currently missing.

RESULTS

To remedy this, we present the Large-Scale Transcriptomic Analysis Pipeline in Kingdom of Life (LSTrAP-Kingdom) pipeline, which we used to process 134 521 RNA-seq samples, achieving ∼12 000 processed samples per day. Our pipeline generated quality-controlled, annotated gene expression matrices that rival the manually curated gene expression data in identifying functionally related genes.

AVAILABILITY AND IMPLEMENTATION

LSTrAP-Kingdom is available from: https://github.com/wirriamm/plants-pipeline and is fully implemented in Python and Bash.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Early "Rootprints" of Plant Terrestrialization: Selaginella Root Development Sheds Light on Root Evolution in Vascular Plants

Roots provide multiple key functions for plants, including anchorage and capturing of water and nutrients. Evolutionarily, roots represent a crucial innovation that enabled plants to migrate from aquatic to terrestrial environment and to grow in height. Based on fossil evidence, roots evolved at least twice independently, once in the lycophyte clade and once in the euphyllophyte (ferns and seed plants) clade. In lycophytes, roots originated in a stepwise manner. Despite their pivotal position in root evolution, it remains unclear how root development is controlled in lycophytes. Getting more insight into lycophyte root development might shed light on how genetic players controlling the root meristem and root developmental processes have evolved. Unfortunately, genetic studies in lycophytes are lagging behind, lacking advanced biotechnological tools, partially caused by the limited economic value of this clade. The technology of RNA sequencing (RNA-seq) at least enabled transcriptome studies, which could enhance the understanding or discovery of genes involved in the root development of this sister group of euphyllophytes. Here, we provide an overview of the current knowledge on root evolution followed by a survey of root developmental events and how these are genetically and hormonally controlled, starting from insights obtained in the model seed plant Arabidopsis and where possible making a comparison with lycophyte root development. Second, we suggest possible key genetic regulators in root development of lycophytes mainly based on their expression profiles in Selaginella moellendorffii and phylogenetics. Finally, we point out challenges and possible future directions for research on root evolution.

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

The appearance of plant organs mediated the explosive radiation of land plants, which shaped the biosphere and allowed the establishment of terrestrial animal life. The evolution of organs and immobile gametes required the coordinated acquisition of novel gene functions, the co-option of existing genes and the development of novel regulatory programmes. However, no large-scale analyses of genomic and transcriptomic data have been performed for land plants. To remedy this, we generated gene expression atlases for various organs and gametes of ten plant species comprising bryophytes, vascular plants, gymnosperms and flowering plants. A comparative analysis of the atlases identified hundreds of organ- and gamete-specific orthogroups and revealed that most of the specific transcriptomes are significantly conserved. Interestingly, our results suggest that co-option of existing genes is the main mechanism for evolving new organs. In contrast to female gametes, male gametes showed a high number and conservation of specific genes, which indicates that male reproduction is highly specialized. The expression atlas capturing pollen development revealed numerous transcription factors and kinases essential for pollen biogenesis and function.

Tissue and regional expression patterns of dicistronic tRNA-mRNA transcripts in grapevine (Vitis vinifera) and their evolutionary co-appearance with vasculature in land plants

Transfer RNAs (tRNA) are crucial adaptor molecules between messenger RNA (mRNA) and amino acids. Recent evidence in plants suggests that dicistronic tRNA-like structures also act as mobile signals for mRNA transcripts to move between distant tissues. Co-transcription is not a common feature in the plant nuclear genome and, in the few cases where polycistronic transcripts have been found, they include non-coding RNA species, such as small nucleolar RNAs and microRNAs. It is not known, however, the extent to which dicistronic transcripts of tRNA and mRNAs are expressed in field-grown plants, or the factors contributing to their expression. We analysed tRNA-mRNA dicistronic transcripts in the major horticultural crop grapevine (Vitis vinifera) using a novel pipeline developed to identify dicistronic transcripts from high-throughput RNA-sequencing data. We identified dicistronic tRNA-mRNA in leaf and berry samples from 22 commercial vineyards. Of the 124 tRNA genes that were expressed in both tissues, 18 tRNA were expressed forming part of 19 dicistronic tRNA-mRNAs. The presence and abundance of dicistronic molecules was tissue and geographic sub-region specific. In leaves, the expression patterns of dicistronic tRNA-mRNAs significantly correlated with tRNA expression, suggesting that their transcriptional regulation might be linked. We also found evidence of syntenic genomic arrangements of tRNAs and protein-coding genes between grapevine and Arabidopsis thaliana, and widespread prevalence of dicistronic tRNA-mRNA transcripts among vascular land plants but no evidence of these transcripts in non-vascular lineages. This suggests that the appearance of plant vasculature and tRNA-mRNA occurred concurrently during the evolution of land plants.

Using Gene Expression to Study Specialized Metabolism-A Practical Guide

Plants produce a vast array of chemical compounds that we use as medicines and flavors, but these compounds' biosynthetic pathways are still poorly understood. This paucity precludes us from modifying, improving, and mass-producing these specialized metabolites in suitable bioreactors. Many of the specialized metabolites are expressed in a narrow range of organs, tissues, and cell types, suggesting a tight regulation of the responsible biosynthetic pathways. Fortunately, with unprecedented ease of generating gene expression data and with >200,000 publicly available RNA sequencing samples, we are now able to study the expression of genes from hundreds of plant species. This review demonstrates how gene expression can elucidate the biosynthetic pathways by mining organ-specific genes, gene expression clusters, and applying various types of co-expression analyses. To empower biologists to perform these analyses, we showcase these analyses using recently published, user-friendly tools. Finally, we analyze the performance of co-expression networks and show that they are a valuable addition to elucidating multiple the biosynthetic pathways of specialized metabolism.

Genome-wide identification and characterization of COMT gene family during the development of blueberry fruit

BACKGROUND

Caffeic acid O-methyltransferases (COMTs) play an important role in the diversification of natural products, especially in the phenylalanine metabolic pathway of plant. The content of COMT genes in blueberry and relationship between their expression patterns and the lignin content during fruit development have not clearly investigated by now.

RESULTS

Ninety-two VcCOMTs were identified in Vaccinium corymbosum. According to phylogenetic analyses, the 92 VcCOMTs were divided into 2 groups. The gene structure and conserved motifs within groups were similar which supported the reliability of the phylogenetic structure groupings. Dispersed duplication (DSD) and whole-genome duplication (WGD) were determined to be the major forces in VcCOMTs evolution. The results showed that the results of qRT-PCR and lignin content for 22 VcCOMTs, VcCOMT40 and VcCOMT92 were related to lignin content at different stages of fruit development of blueberry.

CONCLUSION

We identified COMT gene family in blueberry, and performed comparative analyses of the phylogenetic relationships in the 15 species of land plant, and gene duplication patterns of COMT genes in 5 of the 15 species. We found 2 VcCOMTs were highly expressed and their relative contents were similar to the variation trend of lignin content during the development of blueberry fruit. These results provide a clue for further study on the roles of VcCOMTs in the development of blueberry fruit and could promisingly be foundations for breeding blueberry clutivals with higher fruit firmness and longer shelf life.

Long-distance stress and developmental signals associated with abscisic acid signaling in environmental responses

Flowering plants consist of highly differentiated organs, including roots, leaves, shoots and flowers, which have specific roles: root system for water and nutrient uptake, leaves for photosynthesis and gas exchange and reproductive organs for seed production. The communication between organs through the vascular system, by which water, nutrient and signaling molecules are transported, is essential for coordinated growth and development of the whole plant, particularly under adverse conditions. Here, we highlight recent progress in understanding how signaling pathways of plant hormones are associated with long-distance stress and developmental signals, with particular focus on environmental stress responses. In addition to the root-to-shoot peptide signal that induces abscisic acid accumulation in leaves under drought stress conditions, we summarize the diverse stress-responsive peptide signals reported to date to play a role in environmental responses.

Plant TDP1 (Tyrosyl-DNA Phosphodiesterase 1): A Phylogenetic Perspective and Gene Expression Data Mining

The TDP1 (tyrosyl-DNA phosphodiesterase 1) enzyme removes the non-specific covalent intermediates between topoisomerase I and DNA, thus playing a crucial role in preventing DNA damage. While mammals possess only one TDP1 gene, in plants two genes (TDP1α and TDP1β) are present constituting a small gene subfamily. These display a different domain structure and appear to perform non-overlapping functions in the maintenance of genome integrity. Namely, the HIRAN domain identified in TDP1β is involved in the interaction with DNA during the recognition of stalled replication forks. The availability of transcriptomic databases in a growing variety of experimental systems provides new opportunities to fill the current gaps of knowledge concerning the evolutionary origin and the specialized roles of TDP1 genes in plants. Whereas a phylogenetic approach has been used to track the evolution of plant TDP1 protein, transcriptomic data from a selection of representative lycophyte, eudicots, and monocots have been implemented to explore the transcriptomic dynamics in different tissues and a variety of biotic and abiotic stress conditions. While the phylogenetic analysis indicates that TDP1α is of non-plant origin and TDP1β is plant-specific originating in ancient vascular plants, the gene expression data mining comparative analysis pinpoints for tissue- and stress-specific responses.

Development and Cell Cycle Activity of the Root Apical Meristem in the Fern Ceratopteris richardii

Ferns are a representative clade in plant evolution although underestimated in the genomic era. Ceratopteris richardii is an emergent model for developmental processes in ferns, yet a complete scheme of the different growth stages is necessary. Here, we present a developmental analysis, at the tissue and cellular levels, of the first shoot-borne root of Ceratopteris. We followed early stages and emergence of the root meristem in sporelings. While assessing root growth, the first shoot-borne root ceases its elongation between the emergence of the fifth and sixth roots, suggesting Ceratopteris roots follow a determinate developmental program. We report cell division frequencies in the stem cell niche after detecting labeled nuclei in the root apical cell (RAC) and derivatives after 8 h of exposure. These results demonstrate the RAC has a continuous mitotic activity during root development. Detection of cell cycle activity in the RAC at early times suggests this cell acts as a non-quiescent organizing center. Overall, our results provide a framework to study root function and development in ferns and to better understand the evolutionary history of this organ.

Fungi.guru: Comparative genomic and transcriptomic resource for the fungi kingdom

The fungi kingdom is composed of eukaryotic heterotrophs, which are responsible for balancing the ecosystem and play a major role as decomposers. They also produce a vast diversity of secondary metabolites, which have antibiotic or pharmacological properties. However, our lack of knowledge of gene function in fungi precludes us from tailoring them to our needs and tapping into their metabolic diversity. To help remedy this, we gathered genomic and gene expression data of 19 most widely-researched fungi to build an online tool, fungi.guru, which contains tools for cross-species identification of conserved pathways, functional gene modules, and gene families. We exemplify how our tool can elucidate the molecular function, biological process and cellular component of genes involved in various biological processes, by identifying a secondary metabolite pathway producing gliotoxin in Aspergillus fumigatus, the catabolic pathway of cellulose in Coprinopsis cinerea and the conserved DNA replication pathway in Fusarium graminearum and Pyricularia oryzae. The tool is available at www.fungi.guru.

What can lycophytes teach us about plant evolution and development? Modern perspectives on an ancient lineage

All Evo-Devo studies rely on representative sampling across the tree of interest to elucidate evolutionary trajectories through time. In land plants, genetic resources are well established in model species representing lineages including bryophytes (mosses, liverworts, and hornworts), monilophytes (ferns and allies), and seed plants (gymnosperms and flowering plants), but few resources are available for lycophytes (club mosses, spike mosses, and quillworts). Living lycophytes are a sister group to the euphyllophytes (the fern and seed plant clade), and have retained several ancestral morphological traits despite divergence from a common ancestor of vascular plants around 420 million years ago. This sister relationship offers a unique opportunity to study the conservation of traits such as sporophyte branching, vasculature, and indeterminacy, as well as the convergent evolution of traits such as leaves and roots which have evolved independently in each vascular plant lineage. To elucidate the evolution of vascular development and leaf formation, molecular studies using RNA Seq, quantitative reverse transcription polymerase chain reaction, in situ hybridisation and phylogenetics have revealed the diversification and expression patterns of KNOX, ARP, HD-ZIP, KANADI, and WOX gene families in lycophytes. However, the molecular basis of further trait evolution is not known. Here we describe morphological traits of living lycophytes and their extinct relatives, consider the molecular underpinnings of trait evolution and discuss future research required in lycophytes to understand the key evolutionary innovations enabling the growth and development of all vascular plants.

Characterization of In Vivo Function(s) of Members of the Plant Mitochondrial Carrier Family

Although structurally related, mitochondrial carrier family (MCF) proteins catalyze the specific transport of a range of diverse substrates including nucleotides, amino acids, dicarboxylates, tricarboxylates, cofactors, vitamins, phosphate and H⁺. Despite their name, they do not, however, always localize to the mitochondria, with plasma membrane, peroxisomal, chloroplast and thylakoid and endoplasmic reticulum localizations also being reported. The existence of plastid-specific MCF proteins is suggestive that the evolution of these proteins occurred after the separation of the green lineage. That said, plant-specific MCF proteins are not all plastid-localized, with members also situated at the endoplasmic reticulum and plasma membrane. While by no means yet comprehensive, the in vivo function of a wide range of these transporters is carried out here, and we discuss the employment of genetic variants of the MCF as a means to provide insight into their in vivo function complementary to that obtained from studies following their reconstitution into liposomes.

LSTrAP-Cloud: A User-Friendly Cloud Computing Pipeline to Infer Coexpression Networks

As genomes become more and more available, gene function prediction presents itself as one of the major hurdles in our quest to extract meaningful information on the biological processes genes participate in. In order to facilitate gene function prediction, we show how our user-friendly pipeline, the Large-Scale Transcriptomic Analysis Pipeline in Cloud (LSTrAP-Cloud), can be useful in helping biologists make a shortlist of genes involved in a biological process that they might be interested in, by using a single gene of interest as bait. The LSTrAP-Cloud is based on Google Colaboratory, and provides user-friendly tools that process quality-control RNA sequencing data streamed from the European Nucleotide Archive. The LSTRAP-Cloud outputs a gene coexpression network that can be used to identify functionally related genes for any organism with a sequenced genome and publicly available RNA sequencing data. Here, we used the biosynthesis pathway of Nicotiana tabacum as a case study to demonstrate how enzymes, transporters, and transcription factors involved in the synthesis, transport, and regulation of nicotine can be identified using our pipeline.