The unveiling of long non-coding RNAs as important gene regulators in many biological contexts has increased the demand for efficient and robust computational methods to identify novel long non-coding RNAs from transcripts assembled with high throughput RNA-seq data. Several classes of sequence-based features have been proposed to distinguish between coding and non-coding transcripts. Among them, open reading frame, conservation scores, nucleotide arrangements, and RNA secondary structure have been used with success in literature to recognize intergenic long non-coding RNAs, a particular subclass of non-coding RNAs.

Neurotrophins (NTF) are a family of secreted nerve growth factors with affinity for tyrosine kinase (Ntrk) and p75 receptors. To fully understand the variety of developmental roles played by NTFs, it is critical to know when and where genes encoding individual ligands and receptors are transcribed. Identification of ntf and ntrk transcripts in zebrafish development remains to be fully characterized for further uncovering the potential function(s) of the NTF signal transduction pathway. Here, we conducted a systematic analysis of the expression profiles of four ntf and five ntrk genes during zebrafish development using whole-mount in situ hybridization. Our study unveils new expression domains in the developing embryo, confirms those previously known, and shows that ntf and ntrk genes have different degrees of cell- and tissue-type specificity. The unique and overlapping expression patterns here depicted indicate the coordination of the redundant and divergent functions of NTFs and represent valuable tools for deciphering the molecular pathways involved in the specification and function of embryonic cell types.

Alternative splicing increases neuronal transcriptomic complexity throughout animal phylogeny. To delve into the mechanisms controlling the assembly and evolution of this regulatory layer, we characterized the neuronal microexon program in Drosophila and compared it with that of mammals. In nonvertebrate bilaterians, this splicing program is restricted to neurons by the posttranscriptional processing of the enhancer of microexons (eMIC) domain in Srrm234. In Drosophila, this processing is dependent on regulation by Elav/Fne. eMIC deficiency or misexpression leads to widespread neurological alterations largely emerging from impaired neuronal activity, as revealed by a combination of neuronal imaging experiments and cell type-specific rescues. These defects are associated with the genome-wide skipping of short neural exons, which are strongly enriched in ion channels. We found no overlap of eMIC-regulated exons between flies and mice, illustrating how ancient posttranscriptional programs can evolve independently in different phyla to affect distinct cellular modules while maintaining cell-type specificity.

The anti-proliferative effects of diatoms, described for the first time in copepods, have also been demonstrated in benthic invertebrates such as polychaetes, sea urchins and tunicates. In these organisms PUAs (polyunsaturated aldehydes) induce the disruption of gametogenesis, gamete functionality, fertilization, embryonic mitosis, and larval fitness and competence. These inhibitory effects are due to the PUAs, produced by diatoms in response to physical damage as occurs during copepod grazing. The cell targets of these compounds remain largely unknown. Here we identify some of the genes targeted by the diatom PUA 2-trans-4-trans-decadienal (DD) using the tunicate Ciona intestinalis. The tools, techniques and genomic resources available for Ciona, as well as the suitability of Ciona embryos for medium-to high-throughput strategies, are key to their employment as model organisms in different fields, including the investigation of toxic agents that could interfere with developmental processes. We demonstrate that DD can induce developmental aberrations in Ciona larvae in a dose-dependent manner. Moreover, through a preliminary analysis, DD is shown to affect the expression level of genes involved in stress response and developmental processes.

Epithelial-mesenchymal interactions are crucial for the development of numerous animal structures. Thus, unraveling how molecular tools are recruited in different lineages to control interplays between these tissues is key to understanding morphogenetic evolution. Here, we study Esrp genes, which regulate extensive splicing programs and are essential for mammalian organogenesis. We find that Esrp homologs have been independently recruited for the development of multiple structures across deuterostomes. Although Esrp is involved in a wide variety of ontogenetic processes, our results suggest ancient roles in non-neural ectoderm and regulating specific mesenchymal-to-epithelial transitions in deuterostome ancestors. However, consistent with the extensive rewiring of Esrp-dependent splicing programs between phyla, most developmental defects observed in vertebrate mutants are related to other types of morphogenetic processes. This is likely connected to the origin of an event in Fgfr, which was recruited as an Esrp target in stem chordates and subsequently co-opted into the development of many novel traits in vertebrates.

Tyrosinases, widely distributed among animals, plants and fungi, are involved in the biosynthesis of melanin, a pigment that has been exploited, in the course of evolution, to serve different functions. We conducted a deep evolutionary analysis of tyrosinase family amongst metazoa, thanks to the availability of new sequenced genomes, assessing that tyrosinases (tyr) represent a distinctive feature of all the organisms included in our study and, interestingly, they show an independent expansion in most of the analyzed phyla. Tyrosinase-related proteins (tyrp), which derive from tyr but show distinct key residues in the catalytic domain, constitute an invention of chordate lineage. In addition we here reported a detailed study of the expression territories of the ascidian Ciona intestinalis tyr and tyrps. Furthermore, we put efforts in the identification of the regulatory sequences responsible for their expression in pigment cell lineage. Collectively, the results reported here enlarge our knowledge about the tyrosinase gene family as valuable resource for understanding the genetic components involved in pigment cells evolution and development.

Modern high-throughput technologies enable the processing of a large number of samples simultaneously, while also providing rapid and accurate procedures. In recent years, automated liquid handling workstations have emerged as an established technology for reproducible sample preparation. They offer flexibility, making them suitable for an expanding range of applications. Commonly, such approaches are well-developed for experimental procedures primarily designed for cell-line processing and xenobiotics testing. Conversely, little attention is focused on the application of automated liquid handlers in the analysis of whole organisms, which often involves time-consuming laboratory procedures.

Although splicing occurs largely co-transcriptionally, the order by which introns are removed does not necessarily follow the order in which they are transcribed. Whereas several genomic features are known to influence whether or not an intron is spliced before its downstream neighbor, multiple questions related to adjacent introns' splicing order (AISO) remain unanswered. Here, we present Insplico, the first standalone software for quantifying AISO that works with both short and long read sequencing technologies. We first demonstrate its applicability and effectiveness using simulated reads and by recapitulating previously reported AISO patterns, which unveiled overlooked biases associated with long read sequencing. We next show that AISO around individual exons is remarkably constant across cell and tissue types and even upon major spliceosomal disruption, and it is evolutionarily conserved between human and mouse brains. We also establish a set of universal features associated with AISO patterns across various animal and plant species. Finally, we used Insplico to investigate AISO in the context of tissue-specific exons, particularly focusing on SRRM4-dependent microexons. We found that the majority of such microexons have non-canonical AISO, in which the downstream intron is spliced first, and we suggest two potential modes of SRRM4 regulation of microexons related to their AISO and various splicing-related features. Insplico is available on gitlab.com/aghr/insplico.

Long non-coding RNAs (lncRNAs) represent a novel class of non-coding RNAs having a crucial role in many biological processes. The identification of long non-coding homologs among different species is essential to investigate such roles in model organisms as homologous genes tend to retain similar molecular and biological functions. Alignment-based metrics are able to effectively capture the conservation of transcribed coding sequences and then the homology of protein coding genes. However, unlike protein coding genes the poor sequence conservation of long non-coding genes makes the identification of their homologs a challenging task.

The ascidian neural plate consists of a defined number of identifiable cells organized in a grid of rows and columns, representing a useful model to investigate the molecular mechanisms controlling neural patterning in chordates. Distinct anterior brain lineages are specified via unique combinatorial inputs of signalling pathways with Nodal and Delta-Notch signals patterning along the medial-lateral axis and FGF/MEK/ERK signals patterning along the anterior-posterior axis of the neural plate. The Ciona Gsx gene is specifically expressed in the a9.33 cells in the row III/column 2 position of anterior brain lineages, characterised by a combinatorial input of Nodal-OFF, Notch-ON and FGF-ON. Here, we identify the minimal cis-regulatory element (CRE) of 376 bp, which can recapitulate the early activation of Gsx. We show that this minimal CRE responds in the same way as the endogenous Gsx gene to manipulation of FGF- and Notch-signalling pathways and to overexpression of Snail, a mediator of Nodal signals, and Six3/6, which is required to demarcate the anterior boundary of Gsx expression at the late neurula stage. We reveal that sequences proximal to the transcription start site include a temporal regulatory element required for the precise transcriptional onset of gene expression. We conclude that sufficient spatial and temporal information for Gsx expression is integrated in 376 bp of non-coding cis-regulatory sequences.

The introduction of new genome editing tools such as ZFNs, TALENs and, more recently, the CRISPR/Cas9 system, has greatly expanded the ability to knock-out genes in different animal models, including zebrafish. However, time and costs required for the screening of a huge number of animals, aimed to identify first founder fishes (F0), and then carriers (F1) are still a bottleneck. Currently, high-resolution melting (HRM) analysis is the most efficient technology for large-scale InDels detection, but the very expensive equipment demanded for its application may represent a limitation for research laboratories. Here, we propose a rapid and cheap method for high-throughput genotyping that displays efficiency rate similar to the HRM. In fact, using a common ViiA™7 real-time PCR system and optimizing the parameters of the melting analysis, we demonstrated that it is possible to discriminate between the mutant and the wild type melting curves. Due to its simplicity, rapidity and cheapness, our method can be used as a preliminary one-step approach for massive screening, in order to restrict the scope at a limited number of embryos and to focus merely on them for the next sequencing step, necessary for the exact sequence identification of the induced mutation. Moreover, thanks to its versatility, this simple approach can be readily adapted to the detection of any kind of genome editing approach directed to genes or regulatory regions and can be applied to many other animal models.

RNA splicing is widely dysregulated in cancer, frequently due to altered expression or activity of splicing factors (SFs). Microexons are extremely small exons (3-27 nucleotides long) that are highly evolutionarily conserved and play critical roles in promoting neuronal differentiation and development. Inclusion of microexons in mRNA transcripts is mediated by the SF Serine/Arginine Repetitive Matrix 4 (SRRM4), whose expression is largely restricted to neural tissues. However, microexons have been largely overlooked in prior analyses of splicing in cancer, as their small size necessitates specialized computational approaches for their detection. Here, we demonstrate that despite having low expression in normal nonneural tissues, SRRM4 is further silenced in tumors, resulting in the suppression of normal microexon inclusion. Remarkably, SRRM4 is the most consistently silenced SF across all tumor types analyzed, implying a general advantage of microexon down-regulation in cancer independent of its tissue of origin. We show that this silencing is favorable for tumor growth, as decreased SRRM4 expression in tumors is correlated with an increase in mitotic gene expression, and up-regulation of SRRM4 in cancer cell lines dose-dependently inhibits proliferation in vitro and in a mouse xenograft model. Further, this proliferation inhibition is accompanied by induction of neural-like expression and splicing patterns in cancer cells, suggesting that SRRM4 expression shifts the cell state away from proliferation and toward differentiation. We therefore conclude that SRRM4 acts as a proliferation brake, and tumors gain a selective advantage by cutting off this brake.

Plastic pollution is recognized as a global environmental threat and concern is increasing regarding the potential interactions of the smallest fragments, nanoplastics (1 µm), with either physical and chemical entities encountered in the natural environment, including toxic pollutants. The smallest size of nanoplastics (<100nm) rebounds to their safety associated with remarkable biological, chemical and physical reactivity that allow them to interact with cellular machinery by crossing biological barriers and causing damage to living beings. Recent findings on nanoplastic occurrence in marine coastal waters, including the Mediterranean Sea, leave open the question on their ability to act as a vector of other contaminants of emerging concerns (CECs) concomitantly released by wastewater treatment plants and reaching marine coastal waters. Here, we assess for the first time the role of non-functionalized polystyrene nanoparticles (PS NPs, 20 nm) as a proxy for nanoplastics (1 and 10 µg/mL) alone and in combination with bisphenol A (BPA) (4.5 and 10 µM) on Ciona robusta embryos (22 h post fertilization, hpf) by looking at embryotoxicity through phenotypic alterations. We confirmed the ability of BPA to impact ascidian C. robusta embryo development, by affecting sensory organs pigmentation, either alone and in combination with PS NPs. Our findings suggest that no interactions are taking place between PS NPs and BPA in filtered sea water (FSW) probably due to the high ionic strength of seawater able to trigger the sorption surface properties of PS NPs. Further studies are needed to elucidate such peculiarities and define the risk posed by combined exposure to BPA and PS NPs in marine coastal waters.