Comparative Transcriptomics to Identify RNA Writers and Erasers in Microalgae.

Epitranscriptomics is considered as a new regulatory step in eukaryotes for developmental processes and stress responses. The aim of this study was, for the first time, to identify RNA methyltransferase (writers) and demethylase (erasers) in four investigated species, i.e., the dinoflagellates Alexandrium tamutum and Amphidinium carterae, the diatom Cylindrotheca closterium, and the green alga Tetraselmis suecica. As query sequences for the enzymatic classes of interest, we selected those ones that were previously detected in marine plants, evaluating their expression upon nutrient starvation stress exposure. The hypothesis was that upon stress exposure, the activation/deactivation of specific writers and erasers may occur. In microalgae, we found almost all plant writers and erasers (ALKBH9B, ALKBH10B, MTB, and FIP37), except for three writers (MTA, VIRILIZER, and HAKAI). A sequence similarity search by scanning the corresponding genomes confirmed their presence. Thus, we concluded that the three writer sequences were lacking from the studied transcriptomes probably because they were not expressed in those experimental conditions, rather than a real lack of these genes from their genomes. This study showed that some of them were expressed only in specific culturing conditions. We also investigated their expression in other culturing conditions (i.e., nitrogen depletion, phosphate depletion, and Zinc addition at two different concentrations) in A. carterae, giving new insights into their possible roles in regulating gene expression upon stress.

De novo transcriptome assembly of a lipoxygenase knock-down strain in the diatom Pseudo-nitzschia arenysensis.

Diatoms are microalgae that live in marine and freshwater environments and are responsible for about 20% of the world's carbon fixation. Population dynamics of these cells is finely regulated by intricate signal transduction systems, in which oxylipins are thought to play a relevant role. These are oxygenated fatty acids whose biosynthesis is initiated by a lipoxygenase enzyme (LOX) and are widely distributed in all phyla, including diatoms. Here, we present a de novo transcriptome obtained from the RNA-seq performed in the diatom species Pseudo-nitzschia arenysensis, using both a wild-type and a LOX-silenced strain, which will represent a reliable reference for comparative analyses within the Pseudo-nitzschia genus and at a broader taxonomic scale. Moreover, the RNA-seq data can be interrogated to go deeper into the oxylipins metabolic pathways.

Seagrass genomes reveal ancient polyploidy and adaptations to the marine environment.

We present chromosome-level genome assemblies from representative species of three independently evolved seagrass lineages: Posidonia oceanica, Cymodocea nodosa, Thalassia testudinum and Zostera marina. We also include a draft genome of Potamogeton acutifolius, belonging to a freshwater sister lineage to Zosteraceae. All seagrass species share an ancient whole-genome triplication, while additional whole-genome duplications were uncovered for C. nodosa, Z. marina and P. acutifolius. Comparative analysis of selected gene families suggests that the transition from submerged-freshwater to submerged-marine environments mainly involved fine-tuning of multiple processes (such as osmoregulation, salinity, light capture, carbon acquisition and temperature) that all had to happen in parallel, probably explaining why adaptation to a marine lifestyle has been exceedingly rare. Major gene losses related to stomata, volatiles, defence and lignification are probably a consequence of the return to the sea rather than the cause of it. These new genomes will accelerate functional studies and solutions, as continuing losses of the 'savannahs of the sea' are of major concern in times of climate change and loss of biodiversity.

Gene expression during the formation of resting spores induced by nitrogen starvation in the marine diatom Chaetoceros socialis.

BACKGROUND: Dormancy is widespread in both multicellular and unicellular organisms. Among diatoms, unicellular microalgae at the base of all aquatic food webs, several species produce dormant cells (spores or resting cells) that can withstand long periods of adverse environmental conditions. RESULTS: We present the first gene expression study during the process of spore formation induced by nitrogen depletion in the marine planktonic diatom Chaetoceros socialis. In this condition, genes related to photosynthesis and nitrate assimilation, including high-affinity nitrate transporters (NTRs), were downregulated. While the former result is a common reaction among diatoms under nitrogen stress, the latter seems to be exclusive of the spore-former C. socialis. The upregulation of catabolic pathways, such as tricarboxylic acid cycle, glyoxylate cycle and fatty acid beta-oxidation, suggests that this diatom could use lipids as a source of energy during the process of spore formation. Furthermore, the upregulation of a lipoxygenase and several aldehyde dehydrogenases (ALDHs) advocates the presence of oxylipin-mediated signaling, while the upregulation of genes involved in dormancy-related pathways conserved in other organisms (e.g. serine/threonine-protein kinases TOR and its inhibitor GATOR) provides interesting avenues for future explorations. CONCLUSIONS: Our results demonstrate that the transition from an active growth phase to a resting one is characterized by marked metabolic changes and provides evidence for the presence of signaling pathways related to intercellular communication.

Gene co-expression network analysis for the selection of candidate early warning indicators of heat and nutrient stress in Posidonia oceanica.

The continuous worldwide seagrasses decline calls for immediate actions in order to preserve this precious marine ecosystem. The main stressors that have been linked with decline in seagrasses are 1) the increasing ocean temperature due to climate change and 2) the continuous inputs of nutrients (eutrophication) associated with coastal human activities. To avoid the loss of seagrass populations, an "early warning" system is needed. We used Weighed Gene Co-expression Network Analysis (WGCNA), a systems biology approach, to identify potential candidate genes that can provide an early warning signal of stress in the Mediterranean iconic seagrass Posidonia oceanica, anticipating plant mortality. Plants were collected from both eutrophic (EU) and oligotrophic (OL) environments and were exposed to thermal and nutrient stress in a dedicated mesocosm. By correlating the whole-genome gene expression after 2-weeks exposure with the shoot survival percentage after 5-weeks exposure to stressors, we were able to identify several transcripts that indicated an early activation of several biological processes (BP) including: protein metabolic process, RNA metabolic process, organonitrogen compound biosynthetic process, catabolic process and response to stimulus, which were shared among OL and EU plants and among leaf and shoot apical meristem (SAM), in response to excessive heat and nutrients. Our results suggest a more dynamic and specific response of the SAM compared to the leaf, especially the SAM from plants coming from a stressful environment appeared more dynamic than the SAM from a pristine environment. A vast list of potential molecular markers is also provided that can be used as targets to assess field samples.

Local environment modulates whole-transcriptome expression in the seagrass Posidonia oceanica under warming and nutrients excess.

The intensification of anomalous events of seawater warming and the co-occurrence with local anthropogenic stressors are threatening coastal marine habitats, including seagrasses, which form extensive underwater meadows. Eutrophication highly affects coastal environments, potentially summing up to the widespread effects of global climate changes. In the present study, we investigated for the first time in seagrasses, the transcriptional response of different plant organs (i.e., leaf and shoot apical meristem, SAM) of the Mediterranean seagrass Posidonia oceanica growing in environments with a different history of nutrient enrichment. To this end, a mesocosm experiment exposing plants to single (nutrient enrichment or temperature increase) and multiple stressors (nutrient enrichment plus temperature increase), was performed. Results revealed a differential transcriptome regulation of plants under single and multiple stressors, showing an organ-specific sensitivity depending on plants' origin. While leaf tissues were more responsive to nutrient stress, SAM revealed a higher sensitivity to temperature treatments, especially in plants already impacted in their native environment. The exposure to stress conditions induced the modulation of different biological processes. Plants living in an oligotrophic environment were more responsive to nutrients compared to plants from a eutrophic environment. Evidences that epigenetic mechanisms were involved in the regulation of transcriptional reprogramming were also observed in both plants' organs. These results represent a further step in the comprehension of seagrass response to abiotic stressors pointing out the importance of local pressures in a global warming scenario.

Silencing of a Pseudo-nitzschia arenysensis lipoxygenase transcript leads to reduced oxylipin production and impaired growth.

Because of their importance as chemical mediators, the presence of a rich and varied family of lipoxygenase (LOX) products, collectively named oxylipins, has been investigated thoroughly in diatoms, and the involvement of these products in important processes such as bloom regulation has been postulated. Nevertheless, little information is available on the enzymes and pathways operating in these protists. Exploiting transcriptome data, we identified and characterized a LOX gene, PaLOX, in Pseudo-nitzschia arenysensis, a marine diatom known to produce different species of oxylipins by stereo- and regio-selective oxidation of eicosapentaenoic acid (EPA) at C12 and C15. PaLOX RNA interference correlated with a decrease of the lipid-peroxidizing activity and oxylipin synthesis, as well as with a reduction of growth of P. arenysensis. In addition, sequence analysis and structure models of the C-terminal part of the predicted protein closely fitted with the data for established LOXs from other organisms. The presence in the genome of a single LOX gene, whose downregulation impairs both 12- and 15-oxylipins synthesis, together with the in silico 3D protein modelling suggest that PaLOX encodes for a 12/15S-LOX with a dual specificity, and provides additional support to the correlation between cell growth and oxylipin biosynthesis in diatoms.

Evolutionary Adaptation of the Thyroid Hormone Signaling Toolkit in Chordates.

The specification of the endostyle in non-vertebrate chordates and of the thyroid gland in vertebrates are fundamental steps in the evolution of the thyroid hormone (TH) signaling to coordinate development and body physiology in response to a range of environmental signals. The physiology and biology of TH signaling in vertebrates have been studied in the past, but a complete understanding of such a complex system is still lacking. Non-model species from non-vertebrate chordates may greatly improve our understanding of the evolution of this complex endocrine pathway. Adaptation of already existing proteins in order to perform new roles is a common feature observed during the course of evolution. Through sequence similarity approaches, we investigated the presence of bona fide thyroid peroxidase (TPO), iodothyronine deiodinase (DIO), and thyroid hormone receptors (THRs) in non-vertebrate and vertebrate chordates. Additionally, we determined both the conservation and divergence degrees of functional domains at the protein level. This study supports the hypothesis that non-vertebrate chordates have a functional thyroid hormone signaling system and provides additional information about its possible evolutionary adaptation.

The Ascidia Ciona robusta Provides Novel Insights on the Evolution of the AP-1 Transcriptional Complex.

The Activator Protein-1 transcription factor family (AP-1) transcriptional complex is historically defined as an early response group of transcription factors formed by dimeric complexes of the Jun, Fos, Atf, and Maf bZIP proteins that control cell proliferation and differentiation by regulating gene expression. It has been greatly investigated in many model organisms across metazoan evolution. Nevertheless, its complexity and variability of action made its multiple functions difficult to be defined. Here, we place the foundations for understanding the complexity of AP-1 transcriptional members in tunicates. We investigated the gene members of this family in the ascidian Ciona robusta and identified single copies of Jun, Fos, Atf3, Atf2/7, and Maf bZIP-related factors that could have a role in the formation of the AP-1 complex. We highlight that mesenchyme is a common cellular population where all these factors are expressed during embryonic development, and that, moreover, Fos shows a wider pattern of expression including also notochord and neural cells. By ectopic expression in transgenic embryos of Jun and Fos genes alone or in combination, we investigated the phenotypic alterations induced by these factors and highlighted a degree of functional conservation of the AP-1 complex between Ciona and vertebrates. The lack of gene redundancy and the first pieces of evidence of conserved functions in the control of cell movements and structural organization exerted by these factors open the way for using Ciona as a helpful model system to uncover the multiple potentialities of this highly complex family of bZIP transcription factors.

Insights into phenolic compounds from microalgae: structural variety and complex beneficial activities from health to nutraceutics.

Phenolic compounds (PCs) are a family of secondary metabolites with recognized biological activities making them attractive for the biomedical "red" biotechnology. The development of the eco-sustainable production of natural bioactive metabolites requires using easy cultivable organisms, such as microalgae, which represents one of the most promising sources for biotechnological applications. Microalgae are photosynthetic organisms inhabiting aquatic systems, displaying high levels of biological and functional diversities, and are well-known producers of fatty acids and carotenoids. They are also rich in other families of bioactive molecules e.g. phenolic compounds. Microalgal PCs however are less investigated than other molecular components. This study aims to provide a state-of-art picture of the actual knowledge on microalgal phenolic compounds, reviewing information on the PC content variety and chemodiversity in microalgae, their environmental modulation, and we aim to report discuss data on PC biosynthetic pathways. We report the challenges of promoting microalgae as a relevant source of natural PCs, further enhancing the interests of microalgal "biofactories" for biotechnological applications (i.e. nutraceutical, pharmacological, or cosmeceutical products).

First De Novo Transcriptome of the Copepod Rhincalanus gigas from Antarctic Waters.

Antarctic waters are the largest almost untapped diversified resource of our planet. Molecular resources for Antarctic organisms are very limited and mostly represented by sequences used for species genotyping. In this study, we present the first transcriptome for the copepod Rhincalanus gigas, one of the predominant zooplankton species of Antarctic waters. This transcriptome represents also the first molecular resource for an eucalanoid copepod. The transcriptome is of high quality and completeness. The presence of three predicted genes encoding antifreeze proteins and gene duplication within the glutathione metabolism pathway are suggested as possible adaptations to cope with this harsh environment. The R. gigas transcriptome represents a powerful new resource for investigating the molecular basis associated with polar biological processes and ecology.

m6A RNA Methylation in Marine Plants: First Insights and Relevance for Biological Rhythms.

Circadian regulations are essential for enabling organisms to synchronize physiology with environmental light-dark cycles. Post-transcriptional RNA modifications still represent an understudied level of gene expression regulation in plants, although they could play crucial roles in environmental adaptation. N6-methyl-adenosine (m6A) is the most prevalent mRNA modification, established by "writer" and "eraser" proteins. It influences the clockwork in several taxa, but only few studies have been conducted in plants and none in marine plants. Here, we provided a first inventory of m6A-related genes in seagrasses and investigated daily changes in the global RNA methylation and transcript levels of writers and erasers in Cymodocea nodosa and Zostera marina. Both species showed methylation peaks during the dark period under the same photoperiod, despite exhibiting asynchronous changes in the m6A profile and related gene expression during a 24-h cycle. At contrasting latitudes, Z. marina populations displayed overlapping daily patterns of the m6A level and related gene expression. The observed rhythms are characteristic for each species and similar in populations of the same species with different photoperiods, suggesting the existence of an endogenous circadian control. Globally, our results indicate that m6A RNA methylation could widely contribute to circadian regulation in seagrasses, potentially affecting the photo-biological behaviour of these plants.

Bioinformatics Resources for Plant Abiotic Stress Responses: State of the Art and Opportunities in the Fast Evolving -Omics Era.

Abiotic stresses are among the principal limiting factors for productivity in agriculture. In the current era of continuous climate changes, the understanding of the molecular aspects involved in abiotic stress response in plants is a priority. The rise of -omics approaches provides key strategies to promote effective research in the field, facilitating the investigations from reference models to an increasing number of species, tolerant and sensitive genotypes. Integrated multilevel approaches, based on molecular investigations at genomics, transcriptomics, proteomics and metabolomics levels, are now feasible, expanding the opportunities to clarify key molecular aspects involved in responses to abiotic stresses. To this aim, bioinformatics has become fundamental for data production, mining and integration, and necessary for extracting valuable information and for comparative efforts, paving the way to the modeling of the involved processes. We provide here an overview of bioinformatics resources for research on plant abiotic stresses, describing collections from -omics efforts in the field, ranging from raw data to complete databases or platforms, highlighting opportunities and still open challenges in abiotic stress research based on -omics technologies.

De novo transcriptome of the diatom Cylindrotheca closterium identifies genes involved in the metabolism of anti-inflammatory compounds.

Diatoms are the most diverse and abundant group of phytoplankton species and represent a huge reservoir of marine natural products with possible application for human health. Several diatoms are known to have anticancer, anti-inflammatory, antioxidant and anti-microbial properties, but the compounds responsible of these activities are often still unknown. The diatom Cylindrotheca closterium showed anti-inflammatory properties inhibiting TNFα release in human monocytic leukemia cells. In this study, we present the full transcriptome of C. closterium, and used an -omic approach to identify transcripts coding enzymes that can be involved in the synthesis/degradation of anti-inflammatory compounds. This approach allowed to identify phosphatidylinositol-3-phosphatase, phosphatidylinositol 3-kinase catalytic subunit type 3, phosphatidylinositol N-acetylglucosaminyltransferase subunit A, monogalactosyldiacylglycerol synthase and violaxanthin de-epoxidase, which are known to be involved in anti-inflammatory compound metabolism. When C. closterium was cultured in silica-starvation conditions, selected as stress condition to potentially trigger the synthesis of bioactive metabolites, anti-inflammatory activity was lost and expression levels of the analyzed transcripts were reduced. These data suggested that the control culturing condition was the most active. This study used for the first time a transcriptomic-guided approach to identify enzymes involved in anti-inflammatory compound metabolism, directing future discoveries of marine natural products in microalgae.

Insights into the genome structure of four acetogenic bacteria with specific reference to the Wood-Ljungdahl pathway.

Acetogenic bacteria are obligate anaerobes with the ability of converting carbon dioxide and other one-carbon substrates into acetate through the Wood-Ljungdahl (WL) pathway. These substrates are becoming increasingly important feedstock in industrial microbiology. The main potential industrial application of acetogenic bacteria is the production of metabolites that constitute renewable energy sources (biofuel); such bacteria are of particular interest for this purpose thanks to their low energy requirements for large-scale cultivation. Here, we report new genome sequences for four species, three of them are reported for the first time, namely Acetobacterium paludosum DSM 8237, Acetobacterium tundrae DSM 917, Acetobacterium bakii DSM 8239, and Alkalibaculum bacchi DSM 221123. We performed a comparative genomic analysis focused on the WL pathway's genes and their encoded proteins, using Acetobacterium woodii as a reference genome. The Average Nucleotide Identity (ANI) values ranged from 70% to 95% over an alignment length of 5.4-6.5 Mbp. The core genome consisted of 363 genes, whereas the number of unique genes in a single genome ranged from 486 in A. tundrae to 2360 in A.bacchi. No significant rearrangements were detected in the gene order for the Wood-Ljungdahl pathway however, two species showed variations in genes involved in formate metabolism: A. paludosum harbor two copies of fhs1, and A. bakii a truncated fdhF1. The analysis of protein networks highlighted the expansion of protein orthologues in A. woodii compared to A. bacchi, whereas protein networks involved in the WL pathway were more conserved. This study has increased our understanding on the evolution of the WL pathway in acetogenic bacteria.

Bioinformatics for Marine Products: An Overview of Resources, Bottlenecks, and Perspectives.

The sea represents a major source of biodiversity. It exhibits many different ecosystems in a huge variety of environmental conditions where marine organisms have evolved with extensive diversification of structures and functions, making the marine environment a treasure trove of molecules with potential for biotechnological applications and innovation in many different areas. Rapid progress of the omics sciences has revealed novel opportunities to advance the knowledge of biological systems, paving the way for an unprecedented revolution in the field and expanding marine research from model organisms to an increasing number of marine species. Multi-level approaches based on molecular investigations at genomic, metagenomic, transcriptomic, metatranscriptomic, proteomic, and metabolomic levels are essential to discover marine resources and further explore key molecular processes involved in their production and action. As a consequence, omics approaches, accompanied by the associated bioinformatic resources and computational tools for molecular analyses and modeling, are boosting the rapid advancement of biotechnologies. In this review, we provide an overview of the most relevant bioinformatic resources and major approaches, highlighting perspectives and bottlenecks for an appropriate exploitation of these opportunities for biotechnology applications from marine resources.

Functional conserved non-coding elements among tunicates and chordates.

Non-coding regions with dozens to several hundred base pairs of extreme conservation have been found in all metazoan genomes. The distribution of these conserved non-coding elements (CNE) within and across genomes has suggested that many of them may have roles as transcriptional regulatory elements. A combination of bioinformatics and experimental approaches can be used to identify CNEs with regulatory activity in phylogenetically distant species. Nevertheless, the high divergent rate of genomic sequences of several organisms, such as tunicates, complicates the characterization of these conserved elements and very few examples really may prove their functional activity. We used a comparative approach to facilitate the identification of CNEs among distantly related or highly divergent species and experimentally demonstrated the functional significance of these novel CNEs. We first experimentally tested, in C. robusta and D. rerio transgenic embryos, the regulatory activity of conserved elements associated to genes involved in developmental control among different chordates (Homo sapiens and Danio rerio for vertebrates, Ciona robusta and Ciona savignyi for tunicates and Branchiostoma floridae for cephalochordates). Once demonstrated the cross-species functional conservation of these CNEs, the same gene loci were used as references to locate homologous regions and possible CNEs in available tunicate genomes. Comparison of tunicate-specific and chordate-specific CNEs revealed absence of conservation of the regulatory elements in spite of conservation of regulatory patterns, likely due to evolutionary specification of the respective developmental networks. This result highlights the importance of an integrative in-silico/in-vivo approach to CNEs investigation, encompassing both bioinformatics, essential for putative CNEs identification, and laboratory experiments, pivotal for the understanding of CNEs functionality.

Multilevel comparative bioinformatics to investigate evolutionary relationships and specificities in gene annotations: an example for tomato and grapevine.

BACKGROUND: "Omics" approaches may provide useful information for a deeper understanding of speciation events, diversification and function innovation. This can be achieved by investigating the molecular similarities at sequence level between species, allowing the definition of ortholog and paralog genes. However, the spreading of sequenced genome, often endowed with still preliminary annotations, requires suitable bioinformatics to be appropriately exploited in this framework. RESULTS: We presented here a multilevel comparative approach to investigate on genome evolutionary relationships and peculiarities of two fleshy fruit species of relevant agronomic interest, Solanum lycopersicum (tomato) and Vitis vinifera (grapevine). We defined 17,823 orthology relationships between tomato and grapevine reference gene annotations. The resulting orthologs are associated with the detected paralogs in each species, permitting the definition of gene networks, useful to investigate the different relationships. The reconciliation of the compared collections in terms of an updating of the functional descriptions was also exploited. All the results were made accessible in ComParaLogs, a dedicated bioinformatics platform available at http://biosrv.cab.unina.it/comparalogs/gene/search . CONCLUSIONS: The aim of the work was to suggest a reliable approach to detect all similarities of gene loci between two species based on the integration of results from different levels of information, such as the gene, the transcript and the protein sequences, overcoming possible limits due to exclusive protein versus protein comparisons. This to define reliable ortholog and paralog genes, as well as species specific gene loci in the two species, overcoming limits due to the possible draft nature of preliminary gene annotations. Moreover, reconciled functional descriptions, as well as common or peculiar enzymatic classes and protein domains from tomato and grapevine, together with the definition of species-specific gene sets after the pairwise comparisons, contributed a comprehensive set of information useful to comparatively exploit the two species gene annotations and investigate on differences between species with climacteric and non-climacteric fruits. In addition, the definition of networks of ortholog genes and of associated paralogs, and the organization of web-based interfaces for the exploration of the results, defined a friendly computational bench-work in support of comparative analyses between two species.

Bioinformatics Resources for Plant Genomics: Opportunities and Bottlenecks in The -omics Era.

The sudden exponential increase of biological data concerning genome structure and functionalities, also fostered by the advent of Next Generation Sequencing (NGS) technologies, while expanding the opportunity to highlight still uncovered molecular aspects, challenges bioinformatics in several repects. Data management, processing, updating, dissemination and integration are the major areas of concern. The rapid increase in various omics technologies causes two major issues, which may even appear contrasting: the dissemination of poorly curated datasets, still in the form of raw collections or preliminary draft results, and the fast updating of information that, as a consequence, affects the establishment of stable reliable resources. These issues are mainly caused by the lower rate of bioinformatics in extracting added value information from the large amount of data, when compared to the faster technologies involved in data production. This review describes main bioinformatics resources for plants genomics to underline the heterogeneity of the available collections, coherent with the multifaceted complexity of plant sciences. It aims to provide an in-depth report highlighting bottlenecks that may significantly affect a fluent progress in the field and attempts to suggest possible solutions to the various issues.

Transcriptologs: A Transcriptome-Based Approach to Predict Orthology Relationships.

The detection of orthologs is a key approach in genomics, useful to understand gene evolution and phylogenetic relationships and essential for gene function prediction. However, a reliable annotation of the encoded protein regions is still a limiting aspect in genomics, mainly due to the lack of confirmatory experimental evidence at proteome level. Nevertheless, the current ortholog collections are generally based on protein sequence comparisons, in addition to the availability of large transcriptome sequence collections. We developed Transcriptologs, a method for the prediction of orthologs based on similarities of translated fragments from messenger RNAs of 2 species. We implemented a procedure to extend BLAST-based alignments and to define orthologs based on the Bidirectional Best Hit approach. Results from a test case on Arabidopsis thaliana and Sorghum bicolor transcript collections revealed in some cases outperformance of Transcriptologs in comparison with a classical protein-based analysis in terms of alignment quality, revealing similarities otherwise not detectable.