Search

Contaminants of emerging concern have been increasingly associated with the modulation of the epigenome, leading to potentially inherited and persistent impacts on apical endpoints. Here, we address the performance of the OECD Test No. 236 FET (fish embryo acute toxicity) in the identification of chemicals able to modulate the epigenome. Using zebrafish (Danio rerio) embryos, acute and chronic exposures were performed with the pharmaceutical, simvastatin (SIM), a widely prescribed hypocholesterolemic drug reported to induce inter and transgenerational effects. In the present study, the epigenetic effects of environmentally relevant concentrations of SIM (from 8 ng/L to 2000 ng/L) were addressed following (1) an acute embryo assay based on OECD Test No. 236 FET, (2) a chronic partial life-cycle exposure using adult zebrafish (90 days), and (3) F1 embryos obtained from parental exposed animals. Simvastatin induced significant effects in gene expression of key epigenetic biomarkers (DNA methylation and histone acetylation/deacetylation) in the gonads of exposed adult zebrafish and in 80 hpf zebrafish embryos (acute and chronic parental intergenerational exposure), albeit with distinct effect profiles between biological samples. In the chronic exposure, SIM impacted particularly DNA methyltransferase genes in males and female gonads, whereas in F1 embryos SIM affected mostly genes associated with histone acetylation/deacetylation. In the embryo acute direct exposure, SIM modulated the expression of both genes involved in DNA methylation and histone deacetylase. These findings further support the use of epigenetic biomarkers in zebrafish embryos in a high throughput approach to identify and prioritize epigenome-modulating chemicals.

The Mammalia skin exhibits a rich spectrum of evolutionary adaptations. The pilosebaceous unit, composed of the hair shaft, follicle, and the sebaceous gland, is the most striking synapomorphy. The evolutionary diversification of mammals across different ecological niches was paralleled by the appearance of an ample variety of skin modifications. Pangolins, order Pholidota, exhibit keratin-derived scales, one of the most iconic skin appendages. This formidable armor is intended to serve as a deterrent against predators. Surprisingly, while pangolins have hair on their abdomens, the occurrence of sebaceous and sweat glands is contentious. Here, we explore various molecular modules of skin physiology in four pangolin genomes, including that of sebum production. We show that genes driving wax monoester formation -Awat1/2, show patterns of inactivation in the stem pangolin branch; while the triacylglycerol synthesis geneDgat2l6seems independently eroded in the African and Asian clades. In contrast,Elovl3implicated in the formation of specific neutral lipids required for skin barrier function, is intact and expressed in the skin. An extended comparative analysis shows that genes involved in skin pathogen defense and structural integrity of keratinocyte layers also show inactivating mutations: associated with both ancestral and independent pseudogenization events. Finally, we deduce that the suggested absence of sweat glands is not paralleled by the inactivation of the ATP-binding cassette transporterAbcc11, as previously described in Cetacea. Our findings reveal the sophisticated, convergent and complex history of gene retention and loss as key mechanisms in the evolution of the highly modified mammalian skin phenotypes.

Vestigial organs are historical echoes of past phenotypes. Determining whether a specific organ constitutes a functional or vestigial structure can be a challenging task, given that distinct levels of atrophy may arise between and within lineages. The mammalian pineal gland, an endocrine organ involved in melatonin biorhythmicity, represents a classic example, often yielding contradicting anatomical observations. In Xenarthra (sloths, anteaters and armadillos), a peculiar mammalian order, the presence of a distinct pineal organ was clearly observed in some species (i.e. Linnaeus’s two-toed sloth) but undetected in other closely related species (i.e. brown-throated sloth). In the nine-banded armadillo, contradicting evidence supports either functional or vestigial scenarios. Thus, to untangle the physiological status of the pineal gland in Xenarthra, we used a genomic approach to investigate the evolution of the gene hub responsible for melatonin synthesis and signaling. We show that both synthesis and signaling compartments are eroded and were lost independently. Additionally, by expanding our analysis to 157 mammal genomes we offer a comprehensive view showing that species with very distinctive habitats and lifestyles have convergently evolved a similar phenotype: Cetacea, Pholidota, Dermoptera, Sirenia and Xenarthra. Our findings suggest that the recurrent inactivation of melatonin genes correlates with pineal atrophy, and endorse the use of genomic analyses to ascertain the physiological status of suspected vestigial structures.

The development of Next-Generation Sequencing (NGS) technologies has revolutionized multiple fields of Biology. The ability to sequence DNA and RNA in an automated, parallel and low-cost approach, was key to boost the applicability and acquisition efficiency of omics resources. In this context, the availability of these tools has become indispensable towards our understanding of biodiversity, deepening our knowledge on the distinct complexity levels, from cells to ecosystems. In this study we comprehensively characterised and annotated a whole-body transcriptome of Hediste diversicolor. This annelid worm species belongs to the family Nereididae and inhabits estuarine and lagoon areas on the Atlantic coasts of Europe and North America. Ecologically, this species plays an important role in benthic food webs. The ability of rework sediments through bioturbation activity makes this species essential to the estuarine mechanisms. Here, we used Illumina next-generation sequencing technology, to sequence a total 105 million (M) paired-end (PE) raw reads and produce the first multi-tissue transcriptome assembly of an adult of H. diversicolor. This transcriptome contains 69,335 transcripts with a N50 transcript length of 2,313 bp and a BUSCO gene completeness of 97.7% and 96% (S: 88.2%; D: 7.8%) in Eukaryota and Metazoa lineage-specific profile libraries. Our findings offer a valuable resource for multiple biological studies using this species., This research was partially supported by the ERA-NET BlueBio COFUND Project SIDESTREAM [Grant ID 68], co-funded through national funds provided by FCT [BLUEBIO/0005/2019]; Agencia Estatal de Investigación [PCI2020-111960], NRC – Norwegian Research Council [#311701]. A.G.S. was funded by the Portuguese Foundation for Science and Technology (FCT) under the Grant [SFRH/BD/137935/2018]. Additional strategic funding was provided by FCT [UIDP/04423/2020].

The rapid expansion of high-quality genome assemblies, exemplified by ongoing initiatives such as the Genome-10K and i5k, demands novel automated methods to approach comparative genomics. Of these, the study of inactivating mutations in the coding region of genes, or pseudogenization, as a source of evolutionary novelty is mostly overlooked. Thus, to address such evolutionary/genomic events, a systematic, accurate and computationally automated approach is required. Here, we present PseudoChecker, the first integrated online platform for gene inactivation inference. Unlike the few existing methods, our comparative genomics-based approach displays full automation, a built-in graphical user interface and a novel index, PseudoIndex, for an empirical evaluation of the gene coding status. As a multi-platform online service, PseudoChecker simplifies access and usability, allowing a fast identification of disruptive mutations. An analysis of 30 genes previously reported to be eroded in mammals, and 30 viable genes from the same lineages, demonstrated that PseudoChecker was able to correctly infer 97% of loss events and 95% of functional genes, confirming its reliability. PseudoChecker is freely available, without login required, at http://pseudochecker.ciimar.up.pt.

Three peroxisome proliferator-activated receptor paralogues (PPARα, -β and -γ) are currently recognized in vertebrate genomes. PPARγ is known to modulate nutrition, adipogenesis and immunity in vertebrates. Natural ligands of PPARγ have been proposed; however, the receptor also binds synthetic ligands such as endocrine disruptors. Two paralogues of PPARα and PPARβ have been documented in teleost species, a consequence of the 3R WGD. Recently, two PPARγ paralogue genes were also identified in Astyanax mexicanus. We aimed to determine whether the presence of two PPARγ paralogues is prevalent in other teleost genomes, through genomic and phylogenetic analysis. Our results showed that besides Characiformes, two PPARγ paralogous genes were also identified in other teleost taxa, coinciding with the teleost-specific, whole-genome duplication and with the retention of both genes prior to the separation of the Clupeocephala. To functionally characterize these genes, we used the European sardine (Sardina pilchardus) as a model. PPARγA and PPARγB display a different tissue distribution, despite the similarity of their functional profiles: they are unresponsive to tested fatty acids and other human PPARγ ligands yet yield a transcriptional response in the presence of tributyltin (TBT). This observation puts forward the relevance of comparative analysis to decipher alternative binding architectures and broadens the disruptive potential of man-made chemicals for aquatic species.

The Atlantic chub mackerel, Scomber colias Gmelin, 1789, is a medium-size pelagic fish with substantial importance in the fisheries of the Atlantic Ocean and the Mediterranean Sea. Over the past decade, this species has gained special relevance being one of the main targets of pelagic fisheries in the NE Atlantic. Here, we sequenced and annotated the first high-quality draft genome assembly of S. colias, produced with Pacbio HiFi long reads and Illumina Paired-End short reads. The estimated genome size is 814 Mb distributed into 2,028 scaffolds and 2,093 contigs with an N50 length of 4,19 and 3,34 Mb, respectively. We annotated 27,675 protein-coding genes and the BUSCO analyses indicated high completeness, with 97.3 % of the single-copy orthologs in the Actinopterygii library profile. The present genome assembly represents a valuable resource to address the biology and management of this relevant fishery. Finally, this is the fourth high-quality genome assembly within the Order Scombriformes and the first in the genus Scomber.

18 Pág. Departamento de Reproducción Animal.(INIA), This research was funded by COMPETE 2020 ; Portugal 2020 ; the European Union through the ERDF, grant number 031342 ; FCT through national funds ( PTDC/CTA-AMB/31342/2017 ); and grant RTI2018-093548-BI00 from the Spanish Ministry of Science and Innovation . R.V. is funded by the FCT under the grant SFRH/BD/144786/2019 . The study was also supported by trainers and veterinarians of Oceanogràfic Valencia, especially Carlos Barros-García, for the obtaining and transfer of dolphin ejaculates. S.R.D. is supported by a “Doctorados Industriales 2018” fellowship of Comunidad de Madrid ( IND2018/BIO-9610 ).

Environmental chemicals have been reported to greatly disturb the endocrine and metabolic systems of multiple animal species. A recent example involves the exploitation of the nuclear receptor (NR) heterodimeric pair composed by PPAR/RXR (peroxisome proliferator-activated receptor/retinoid X receptor), which shows lipid perturbation in mammalian species. While gene orthologues of both of these receptors have been described outside vertebrates, no functional characterization of PPAR has been carried in protostome lineages. We provide the first functional analysis of PPAR in Patella sp. (Mollusca), using model obesogens such as tributyltin (TBT), triphenyltin (TPT), and proposed natural ligands (fatty acid molecules). To gain further insights, we used site-directed mutagenesis to PPAR and replaced the tyrosine 277 by a cysteine (the human homologous amino acid and TBT anchor residue) and an alanine. Additionally, we explored the alterations in the fatty acid profiles after an exposure to the model obesogen TBT, in vivo. Our results show that TBT and TPT behave as an antagonist of Patella sp. PPAR/RXR and that the tyrosine 277 is important, but not essential in the response to TBT. Overall, these results suggest a relation between the response of the mollusc PPAR-RXR to TBT and the lipid profile alterations reported at environmentally relevant concentrations. Our findings highlight the importance of comparative analysis between protostome and deuterostome lineages to decipher the differential impact of environmental chemicals.

Globally persistent man-made chemicals display ever-growing ecosystemic consequences, a hallmark of the Anthropocene epoch. In this context, the assessment of how lineage-specific gene repertoires influence organism sensitivity toward endocrine disruptors is a central question in toxicology. A striking example highlights the role of a group of compounds known as obesogens. In mammals, most examples involve the modulation of the nuclear receptor peroxisome proliferator-activated receptor γ (PPARγ). To address the structural and biological determinants of PPARγ exploitation by a model obesogen, tributyltin (TBT), in chordates, we employed comparative genomics, transactivation and ligand binding assays, homology modeling, and site-directed-mutagenesis. We show that the emergence of multiple PPARs (α, β and γ) in vertebrate ancestry coincides with the acquisition of TBT agonist affinity, as can be deduced from the conserved transactivation and binding affinity of the chondrichthyan and mammalian PPARγ. The amphioxus single-copy PPAR is irresponsive to TBT; as well as the investigated teleosts, this is a probable consequence of a specific mutational remodeling of the ligand binding pocket. Our findings endorse the modulatory ability of man-made chemicals and suggest an evolutionarily diverse setting, with impacts for environmental risk assessment.

The hypothesis that exposure to certain environmental chemicals during early life stages may disrupt reproduction across multiple non-exposed generations has significant implications for understanding disease etiology and adverse outcomes. We demonstrate here reproductive multi and transgenerational effects, at environmentally relevant levels, of one of the most prescribed human pharmaceuticals, simvastatin, in a keystone species, the amphipod Gammarus locusta. The transgenerational findings has major implications for hazard and risk assessment of pharmaceuticals and other contaminants of emerging concern given that transgenerational effects of environmental chemicals are not addressed in current hazard and risk assessment schemes. Considering that the mevalonate synthesis, one of the key metabolic pathways targeted by simvastatin, is highly conserved among metazoans, these results may also shed light on the potential transgenerational effects of simvastatin on other animals, including humans. This research was funded by COMPETE 2020, Portugal 2020, the European Union through the ERDF and the Portuguese Foundation for Science and Technology – FCT (Transobesogen project – Trans-phyletic obesogenic responses: from epigenetic modules to transgenerational environmental impacts, reference: PTDC/CTA-AMB/31544/2017 – NORTE-01-0145-FEDER-031544). This research was also supported by the National Funds through FCT under the projects (UIDB/04423/2020; UIDP/04423/2020), by the Spanish Agencia Estatal de Investigación (CTM2017-84763-C3-2-R), and by the Galician Council of Culture, Education and Universities (ED431C2017/36), cofounded by ERDF. A PhD grant awarded to Susana Barros (PD/BD/143090/2018) was funded by the FCT SI

The ancestors of Cetacea underwent profound morpho-physiological alterations. By displaying an exclusive aquatic existence, cetaceans evolved unique patterns of locomotor activity, vigilant behaviour, thermoregulation and circadian rhythmicity. Deciphering the molecular landscape governing many of these adaptations is key to understand the evolution of phenotypes. Here, we investigate Cortistatin (CORT), a neuropeptide displaying an important role in mammalian biorhythm regulation. This neuropeptide is a known neuroendocrine factor, stimulating slow-wave sleep, but also involved in the regulation of energy metabolism and hypomotility inducement. We assessed the functional status of CORT in 359 mammalian genomes (25 orders), including 30 species of Cetacea. Our findings indicate that cetaceans and other mammals with atypical biorhythms, thermal constraints and/or energy metabolism, have accumulated deleterious mutations in CORT. In light of the pleiotropic action of this neuropeptide, we suggest that this inactivation contributed to a plethora of phenotypic adjustments to accommodate adaptive solutions to specific ecological niches.

Biological control of key processes, such as development and reproduction, is largely ascribed to a superfamily of ligand-dependent and independent transcription factors named Nuclear Receptors (NRs). Given their ability to accommodate ligands, NRs are prime targets of man-made compounds that mimic or antagonise the action of endogenous ligands. Accordingly, NRs occupy a prominent role in OECD and EPA guidelines for testing and assessment of Endocrine disrupting chemicals (EDCs). Although NR assays are already a key instrument in the OECD Conceptual Framework for Testing and Assessment of EDCs, the focus is mostly on vertebrate NRs. Here, we address the chief knowledge gaps in the field. More specifically, we (1) verify the growing availability of genomes/transcriptome projects, (2) highlight gaps in the identification and characterization of metazoan NR and in the establishment of (3) life cycle and (4) toxicity testing protocols. An overall bias towards vertebrates and selected invertebrate groups, notably Arthropoda, Annelida and Mollusca, was observed. Hence, if we aim to improve risk assessment of EDCs and emerging pollutants at an ecosystems scale, and understand their mode of action (MOA), we must establish a framework to include a broad phylogenetic sampling of Metazoans.

Signalling molecules and their cognate receptors are central components of the Metazoa endocrine system. Defining their presence or absence in extant animal lineages is critical to accurately devise evolutionary patterns, physiological shifts and the impact of endocrine disrupting chemicals. Here, we address the evolution of retinoic acid (RA) signalling in the Priapulida worm, Priapulus caudatus Lamarck, 1816, an Ecdysozoa. RA signalling has been shown to be central to chordate endocrine homeostasis, participating in multiple developmental and physiological processes. Priapulids, with their slow rate of molecular evolution and phylogenetic position, represent a key taxon to investigate the early phases of Ecdysozoa evolution. By exploring a draft genome assembly, we show, by means of phylogenetics and functional assays, that an orthologue of the nuclear receptor retinoic acid receptor (RAR) subfamily, a central mediator of RA signalling, is present in Ecdysozoa, contrary to previous perception. We further demonstrate that the Priapulida RAR displays low-affinity for retinoids (similar to annelids), and is not responsive to common endocrine disruptors acting via RAR. Our findings provide a timeline for RA signalling evolution in the Bilateria and give support to the hypothesis that the increase in RA affinity towards RAR is a late acquisition in the evolution of the Metazoa.

The expression pattern of genes coding for enzymes of the retinoic acid (RA) synthetic and degradation pathways was characterized in adult female zebrafish Danio rerio. Females were conditioned until maturation and post-spawn expression dynamics were determined. A striking upregulation of cyp26b1, but not cyp26a1, was observed following egg deposition, decreasing to initial levels during recovery. A similar, yet lower, fluctuation was observed for aldh1a2 and rdh10a, the enzymes participating in the two-step RA biosynthesis cascade. The present work highlights the dynamics of the adult D. rerio oogenesis and uncovers novel, yet elusive, metabolic contributors. Possible compartmentalized roles for the different gene paralogue isoforms are discussed.

Recent reports raise the concern that exposure to several environmental chemicals may induce persistent changes that go beyond the exposed organisms, being transferred to subsequent generations even in the absence of the original chemical insult. These changes in subsequent non-exposed generations have been related to epigenetic changes. Although highly relevant for hazard and risk assessment, biomarkers of epigenetic modifications that can be associated with adversity, are still not integrated into hazard assessment frameworks. Here, in order to validate new biomarkers of epigenetic modifications in a popular animal model, zebrafish embryos were exposed to different concentrations of Bisphenol A (0.01, 0.1, 1 and 10 mg/L) and Valproic Acid (0.8, 4, 20 and 100 mg/L), two chemicals reported to alter the modulation of the epigenome. Morphological abnormalities and epigenetic changes were assessed at 80 hours-post fertilization, including DNA global methylation and gene expression of both DNA and histone epigenetic modifications. Gene expression changes were detected at concentrations below those inducing morphological abnormalities. These results further support the importance of combining epigenetic biomarkers with apical endpoints to improve guidelines for chemical testing and hazard assessment, and favour the integration of new biomarkers of epigenetic modifications into the standardized OECD test guideline 236 with zebrafish embryos.

In the present work we provide the first isolation and functional characterization of a RXR orthologue in an annelid species, the Platynereis dumerilii. Using an in vitro luciferase reporter assay we evaluated the annelid receptor ability to respond to ligand 9-cis-retinoic acid, TBT and TPT. Our results show that the annelid RXR responds to 9-cis-retinoic acid and to the organotins by activating reporter gene transcription. The findings suggest a conserved mode of action of the receptor in response to a common signaling molecule and modulation by organotin compounds between vertebrates and Lophotrochozoans. © 2017 Elsevier This work was supported by the Funda??o para a Ci?ncia e a Tecnologia (FCT) projects (PTDC/MAR/105199/2008 and EXPL/MAR/EST/1540/2012) and by Norte 2020 and FEDER (Coral - Sustainable Ocean Exploitation - Norte-01-0145-FEDER-000036). A Postdoctoral grant awarded to Raquel Ruivo (SFRH/BPD/72519/2010), a Research Technician grant awarded to In?s P?scoa (EXPL/MAR/EST/1540/2012/BTI/2013/015), and a PhD grant awarded to Ana Andr? (SFRH/BD/81243/2011) and to Ana Capit?o (SFRH/BD/90664/2012).

The appearance of mammalian-specific skin features was a key evolutionary event contributing for the elaboration of physiological processes such as thermoregulation, adequate hydration, locomotion and inflammation. Skin inflammatory and autoimmune processes engage a population of skin-infiltrating T cells expressing a specific C-C chemokine receptor (CCR10), which interacts with an epidermal CC chemokine, the skin-specific C-C motif chemokine ligand 27 (CCL27). CCL27 is selectively produced in the skin by keratinocytes, particularly upon inflammation, mediating the adhesion and homing of skin-infiltrating T cells. Here, we examined the evolution and coding condition of Ccl27 in 112 placental mammalian species. Our findings reveal that a number of open reading frame inactivation events such as insertions, deletions, start and stop codon mutations, independently occurred in Cetacea, Pholidota, Sirenia, Chiroptera, and Rodentia, totalizing 18 species. The diverse habitat settings and life-styles of Ccl27-eroded lineages probably implied distinct evolutionary triggers rendering this gene unessential. For example, in Cetacea the rapid renewal of skin layers minimizes the need for an elaborate inflammatory mechanism, mirrored by the absence of epidermal scabs. Our findings suggest that the convergent and independent loss of Ccl27 in mammalian evolution concurred with unique adaptive roads for skin physiology.

Publisher's version (útgefin grein), Genomes are dynamic biological units, with processes of gene duplication and loss triggering evolutionary novelty. The mammalian skin provides a remarkable case study on the occurrence of adaptive morphological innovations. Skin sebaceous glands (SGs), for instance, emerged in the ancestor of mammals serving pivotal roles, such as lubrication, waterproofing, immunity, and thermoregulation, through the secretion of sebum, a complex mixture of various neutral lipids such as triacylglycerol, free fatty acids, wax esters, cholesterol, and squalene. Remarkably, SGs are absent in a few mammalian lineages, including the iconic Cetacea. We investigated the evolution of the key molecular components responsible for skin sebum production: Dgat2l6, Awat1, Awat2, Elovl3, Mogat3, and Fabp9. We show that all analyzed genes have been rendered nonfunctional in Cetacea species (toothed and baleen whales). Transcriptomic analysis, including a novel skin transcriptome from blue whale, supports gene inactivation. The conserved mutational pattern found in most analyzed genes, indicates that pseudogenization events took place prior to the diversification of modern Cetacea lineages. Genome and skin transcriptome analysis of the common hippopotamus highlighted the convergent loss of a subset of sebum-producing genes, notably Awat1 and Mogat3. Partial loss profiles were also detected in non-Cetacea aquatic mammals, such as the Florida manatee, and in terrestrial mammals displaying specialized skin phenotypes such as the African elephant, white rhinoceros and pig. Our findings reveal a unique landscape of “gene vestiges” in the Cetacea sebum-producing compartment, with limited gene loss observed in other mammalian lineages: suggestive of specific adaptations or specializations of skin lipids., This work was supported by Project No. 031342 cofinanced by COMPETE 2020, Portugal 2020 and the European Union through the ERDF, and by Fundac¸a~o para a Cie^ncia e a Tecnologia through national funds. R.R.F. thanks the Danish National Research Foundation for its support of the Center for Macroecology, Evolution, and Climate (grant DNRF96). We acknowledge the various Cetacea genome consortiums for genome sequencing and assemblies. We also thank Gısli Vikingsson at the Marine and Freshwater Research Institute in Iceland for lending us the Larsen gun and to North Sailing whale watching for the use of their zodiac.

The use of transcriptomics data brings new insights and works as a powerful tool to explore the molecular mode of action (MoA) of transgenerational inheritance effects of contaminants of emerging concern. Therefore, in this dataset, we present the transcriptomic data of the transgenerational effects of environmentally relevant simvastatin levels, one of the most prescribed human pharmaceuticals, in the keystone amphipod species Gammarus locusta. In summary, G. locusta juveniles were maintained under simvastatin exposure up to adulthood (exposed group - F0E) and the offspring of F0E were transferred to control water for the three subsequent generations (transgenerational group - F1T, F2T and F3T). To gain insights into the biological functions and canonical pathways transgenerationally disrupted by simvastatin, a G. locusta de novo transcriptome assembly was produced and the transcriptomic profiles of three individual G. locusta females, per group, over the four generations (F0 to F3) - solvent control groups (F0.C, F1.C, F2.C and F3.C), F0 320 ng/L simvastatin exposed group (F0.320E) and F1 to F3 320 transgenerational group (F1.320T; F2.320T and F3.320T) - were analyzed. Briefly, Illumina HiSeq™ 2500 platform was used to perform RNA sequencing, and due to the unavailability of G. locusta genome, the RNA-seq datasets were assembled de novo using Trinity and annotated with Trinotate software. After assembly and post-processing steps, 106093 transcripts with N50 of 2371 bp and mean sequence length of 1343.98 bp was produced. BUSCO analyses showed a transcriptome with gene completeness of 97.5 % Arthropoda library profile. The Bowtie2, RSEM and edgeR tools were used for the differential gene expression (DEGs) analyses that allowed the identification of a high quantity of genes differentially expressed in all generations. Finally, to identify the main metabolic pathways affected by the transgenerational effects of SIM across all generations, the DGEs genes were blasted onto KEGG pathways database using the KAAS webserver. The data furnished in this article allows a better molecular understanding of the transgenerational effects produced by simvastatin in the keystone amphipod G. locusta and has major implications for hazard and risk assessment of pharmaceuticals and other emerging contaminants. This article is related to the research article entitled “Transgenerational inheritance of chemical-induced signature: a case study with simvastatin This article was developed under the Transobesogen project - Trans-phyletic obesogenic responses: from epigenetic modules to transgenerational environmental impacts (reference PTDC/CTA-AMB/31544/2017 - NORTE-01-0145-FEDER-031544), cofunded by Portugal 2020, the European Union through the ERDF and by the Portuguese Foundation for Science and Technology – FCT. This article was also supported by FCT through national funds (UIDB/04423/2020; UIDP/04423/2020), by the Spanish Agencia Estatal de Investigación (CTM2017-84763-C3-2-R) and by the Galician Council of Culture, Education and Universities (ED431C2017/36), cofounded by ERDF. A PhD grant awarded to Susana Barros acknowledges the doctoral grant attributed by FCT with reference PD/BD/143090/2018 SI

The origin and diversification of the metazoan endocrine systems represents a fundamental research issue in biology. Nuclear receptors are critical components of these systems. A particular group named VDR/PXR/CAR (NR1I/J) is central in the mediation of detoxification responses. While orthologues have been thoroughly characterized in vertebrates, a sparse representation is currently available for invertebrates. Here, we provide the first isolation and characterization of a lophotrochozoan protostome VDR/PXR/CAR nuclear receptor (NR1J), in the estuarine bivalve the peppery furrow shell (Scrobicularia plana). Using a reporter gene assay, we evaluated the xenobiotic receptor plasticity comparing the human PXR with the S. plana NR1Jβ. Our results show that the molluscan receptor responds to a natural toxin (okadaic acid) in a similar fashion to that reported for other invertebrates. In contrast, the pesticide esfenvalerate displayed a unique response, since it down regulated transactivation at higher concentrations, while for triclosan no response was observed. Additionally, we uncovered lineage specific gene duplications and gene loss in the gene group encoding NRs in protostomes with likely impacts on the complexity of detoxification mechanisms across different phyla. Our findings pave the way for the development of multi-specific sensor tools to screen xenobiotic compounds acting via the NR1I/J group.

Nuclear receptors (NRs) are key transcription factors that originated in the common ancestor of metazoans. The vast majority of NRs are triggered by binding to either endogenous (e.g. retinoic acid) or exogenous (e.g. xenobiotics) ligands, and their evolution and expansion is tightly linked to the function of endocrine systems. Importantly, they represent classic targets of physiological exploitation by endocrine disrupting chemicals. The NR gene repertoire in different lineages has been shaped by gene loss, duplication and mutation, denoting a dynamic evolutionary route. As the earliest diverging class of gnathostomes (jawed vertebrates), cartilaginous fishes offer an exceptional opportunity to address the early diversification of NR gene families and the evolution of the endocrine system in jawed vertebrates. Here we provide an exhaustive analysis into the NR gene composition in five elasmobranch (sharks and rays) and two holocephalan (chimaeras) species. For this purpose, we generated also a low coverage draft genome assembly of the chimaera small-eyed rabbitfish, Hydrolagus affinis. We show that cartilaginous fish retain an archetypal NR gene repertoire, similar to that of mammals and coincident with the two rounds of whole genome duplication that occurred in the gnathostome ancestor. Furthermore, novel gene members of the non-canonical NR0B receptors were found in the genomes of this lineage. Our findings provide an essential view into the early diversification of NRs in gnathostomes, paving the way for functional studies.

The wide ecological relevance of lipid homeostasis modulators in the environment has been increasingly acknowledged. Tributyltin (TBT), for instance, was shown to cause lipid modulation, not only in mammals, but also in fish, molluscs, arthropods and rotifers. In vertebrates, TBT is known to interact with a nuclear receptor heterodimer module, formed by the retinoid X receptor (RXR) and the peroxisome proliferator-activated receptor (PPAR). These modulate the expression of genes involved in lipid homeostasis. In the present work, we isolated for the first time the complete coding region of the Echinodermata (Paracentrotus lividus) gene orthologues of PPAR and RXR and evaluated the ability of a model lipid homeostasis modulator, TBT, to interfere with the lipid metabolism in this species. Our results demonstrate that TBT alters the gonadal fatty acid composition and gene expression patterns: yielding sex-specific responses in fatty acid levels, including the decrease of eicosapentaenoic acid (C20:5 n-3, EPA) in males, and increase of arachidonic acid (20:4n-6, ARA) in females, and upregulation of long-chain acyl-CoA synthetase (acsl), ppar and rxr. Furthermore, an in vitro test using COS-1 cells as host and chimeric receptors with the ligand binding domain (LBD) of P. lividus PPAR and RXR shows that organotins (TBT and TPT (Triphenyltin)) suppressed activity of the heterodimer PPAR/RXR in a concentration-dependent manner. Together, these results suggest that TBT acts as a lipid homeostasis modulator at environmentally relevant concentrations in Echinodermata and highlight a possible conserved mode of action via the PPAR/RXR heterodimer.

Clupeiformes, such as sardines and herrings, represent an important share of worldwide fisheries. Among those, the European sardine (Sardina pilchardus, Walbaum 1792) exhibits significant commercial relevance. While the last decade showed a steady and sharp decline in capture levels, recent advances in culture husbandry represent promising research avenues. Yet, the complete absence of genomic resources from sardine imposes a severe bottleneck to understand its physiological and ecological requirements. We generated 69 Gbp of paired-end reads using Illumina HiSeq X Ten and assembled a draft genome assembly with an N50 scaffold length of 25579 bp and BUSCO completeness of 82.1% (Actinopterygii). The estimated size of the genome ranges between 655 and 850 Mb. Additionally, we generated a relatively high-level liver transcriptome. To deliver a proof of principle of the value of this dataset, we established the presence and function of enzymes (elovl2, elovl5 and fads2) that have pivotal roles in the biosynthesis of long chain polyunsaturated fatty acids, essential nutrients particularly abundant in oily fish such as sardines. Our study provides the first omics dataset from a valuable economic marine teleost species, the European sardine, an essential resource for their effective conservation, management and sustainable exploitation.

Retinoid acid receptor (RAR)-dependent signalling pathways are essential for the regulation and maintenance of essential biological functions and are recognized targets of disruptive anthropogenic compounds. Recent studies put forward the inability of mollusc RARs to bind and respond to the canonical vertebrate ligand, retinoic acid: a feature that seems to have been lost during evolution. Yet, these studies were carried out in a limited number of molluscs. Therefore, using an in vitro transactivation assay, the present work aimed to characterize phylogenetically relevant mollusc RARs, as monomers or as functional units with RXR, not only in the presence of vertebrate bone fine ligands but also known endocrine disruptors, described to modulate retinoid-dependent pathways. In general, none of the tested mollusc RARs were able to activate reporter gene transcription when exposed to retinoic acid isomers, suggesting that the ability to respond to retinoic acid was lost across molluscs. Similarly, the analysed mollusc RAR were unresponsive towards organochloride pesticides. In contrast, transcriptional repressions were observed with the RAR/RXR unit upon exposure to retinoids or RXR-specific ligands. Loss-of-function and gain-of-function mutations further corroborate the obtained results and suggest that the repressive behaviour, observed with mollusc and human RAR/RXR heterodimers, is possibly mediated by ligand biding to RXR.

To compare genelociconsidering a phylogenetic framework is a promising approach to uncover the genetic basis of human diseases. Imbalance of dopaminergic systems is suspected to underlie some emerging neurological disorders. The physiological functions of dopamine are transduced via G-protein-coupled receptors, including DRD5which displays a relatively higher affinity towards dopamine. Importantly, DRD5knockout mice are hypertense, a condition emerging from an increase in sympathetic tone. We investigated the evolution of DRD5, a high affinity receptor for dopamine, in mammals. Surprisingly, among 124 investigated mammalian genomes, we found that Cetacea lineages (Mysticeti and Odontoceti) have independently lost this gene, as well as the burrowingChrysochloris asiatica (Cape golden mole).We suggest that DRD5inactivation parallels hypoxia-induced adaptations, such as peripheral vasoconstriction required for deep-diving in Cetacea, in accordance with the convergent evolution of vasoconstrictor genes in hypoxia-exposed animals. Our findings indicate that Cetacea are natural knockouts for DRD5and might offer valuable insights into the mechanisms of some forms of vasoconstriction responses and hypertension in humans.

The Cetacea infraorder comprises a very unique group within the mammalian lineage. While sharing common ancestors with terrestrial mammals, their exclusive dependence on aquatic environments makes them attractive models to explore the landscape of molecular shifts in radical habitat transitions. Among their diverse anatomical and physiological solutions, we find detectable genetic remodeling of the immune system. In agreement, here we show that the gene sequence of interleukin-20 (IL20) displays unambiguous signs of inactivation with several disruptive mutations, including stop codons, insertions, and a conserved trans-species mutation abolishing a canonical splice site, in nine analyzed cetacean genomes. Considering the suggested role of IL20 in skin immunity processes, including inflammation, epithelization, and remodeling, we propose that gene inactivation follows specific adaptations of cetacean skin to the aquatic environment, in frame with the less-is-more hypothesis.

Fatty acids (FAs) constitute a considerable fraction of all lipid molecules with a fundamental role in numerous physiological processes. In animals, the majority of complex lipid molecules are derived from the transformation of FAs through several biochemical pathways. Yet, for FAs to enroll in these pathways they require an activation step. FA activation is catalyzed by the rate limiting action of Acyl-CoA synthases. Several Acyl-CoA enzyme families have been previously described and classified according to the chain length of FAs they process. Here, we address the evolutionary history of the ACSBG gene family which activates, FAs with >16 carbons. Currently, two different ACSBG gene families, ACSBG1 and ACSBG2, are recognized in vertebrates. We provide evidence that a wider and unequal ACSBG gene repertoire is present in vertebrate lineages. We identify a novel ACSBG-like gene lineage which occurs specifically in amphibians, ray finned fishes, coelacanths and cartilaginous fishes named ACSBG3. Also, we show that the ACSBG2 gene lineage duplicated in the Theria ancestor. Our findings, thus offer a far richer understanding on FA activation in vertebrates and provide key insights into the relevance of comparative and functional analysis to perceive physiological differences, namely those related with lipid metabolic pathways. © 2018 Elsevier This work was supported by Project INNOVMAR—Innovation and Sustainability in the Management and Exploitation of Marine Resources (reference NORTE-01-0145-FEDER-000035 , within Research Line NOVELMAR/INSEAFOOD/ECOSERVICES), supported by North Portugal Regional Operational Programme ( NORTE 2020 ), under the PORTUGAL 2020 Partnership Agreement, through the European Regional Development Fund (ERDF). We acknowledge Fundação para a Ciência e a Tecnologia for the support to EF (SFRH/BD/79305/2011).

Nucleobase-Ascorbate Transporter (NAT) family includes ascorbic acid, nucleobases and uric acid transporters, with a broad evolutionary distribution. In vertebrates, four members have been previously recognized, the ascorbate transporters Slc23a1 and Slc3a2, the nucleobase transporter Slc23a4 and an orphan transporter SLC23A3. Here we identify a fifth member of the vertebrateslc23complement (slc23a5), expressed in neopterygians (gars and teleosts) and amphibians, and clarify the evolutionary relationships between the novel gene and knownslc23genes. Further comparative analysis puts forward uric acid as the preferred substrate for Slc23a5. Gene expression quantification suggests kidney and testis as major expression sites inXenopus tropicalis(western clawed frog) andDanio rerio(zebrafish). Additional expression in brain was detected inD. rerio, while in the NeoteleosteiOryzias latipes(medaka)slc23a5expression is restricted to brain. The biological relevance of the retention of an extra transporter in fish and amphibians is examined: with respect to the (1) antioxidant role of uric and ascorbic acid in seminal fluid and brain, (2) the ability to endogenously synthesize ascorbic acid and (3) the morphological adaptations of the male urogenital system.

[Background] Provision of long-chain polyunsaturated fatty acids (LC-PUFA) in vertebrates occurs through the diet or via endogenous production from C18 precursors through consecutive elongations and desaturations. It has been postulated that the abundance of LC-PUFA in the marine environment has remarkably modulated the gene complement and function of Fads in marine teleosts. In vertebrates two fatty acyl desaturases, namely Fads1 and Fads2, encode ∆5 and ∆6 desaturases, respectively. To fully clarify the evolutionary history of LC-PUFA biosynthesis in vertebrates, we investigated the gene repertoire and function of Fads from species placed at key evolutionary nodes., [Results] We demonstrate that functional Fads1Δ5 and Fads2∆6 arose from a tandem gene duplication in the ancestor of vertebrates, since they are present in the Arctic lamprey. Additionally, we show that a similar condition was retained in ray-finned fish such as the Senegal bichir and spotted gar, with the identification of fads1 genes in these lineages. Functional characterisation of the isolated desaturases reveals the first case of a Fads1 enzyme with ∆5 desaturase activity in the Teleostei lineage, the Elopomorpha. In contrast, in Osteoglossomorpha genomes, while no fads1 was identified, two separate fads2 duplicates with ∆6 and ∆5 desaturase activities respectively were uncovered., [Conclusions] We conclude that, while the essential genetic components involved LC-PUFA biosynthesis evolved in the vertebrate ancestor, the full completion of the LC-PUFA biosynthesis pathway arose uniquely in gnathostomes., We acknowledge Fundação para a Ciência e a Tecnologia for the support to M.L.-M (SFRH/BD/84238/2012). LFCC research is supported by Norte2020 and FEDER (Coral - Sustainable Ocean Exploitation - Norte-01-0145-FEDER-000036) a project supported by the North Portugal Regional Operational Program (NORTE 2020), under the PORTUGAL 2020 Partnership Agreement, through the European Regional Development Fund (ERDF). The funding bodies played no role in the design of the study and collection, analysis, interpretation of data, or in writing the manuscript.

Recent genomic and biochemical evidence in invertebrate species pushes back the origin of the retinoid metabolic and signaling modules to the last common ancestor of all bilaterians. However, the evolution of retinoid pathways are far from fully understood. In the majority of non-chordate invertebrate lineages, the ongoing functional characterization of retinoid-related genes (metabolism and signaling pathways), as well as the characterization of the endogenous retinoid content (precursors and active retinoids), is still incomplete. Despite limited, the available data supports the presence of biologically active retinoid pathways in invertebrates. Yet, the mechanisms controlling the spatial and temporal distribution of retinoids as well as their physiological significance share similarities and differences with vertebrates. For instance, retinol storage in the form of retinyl esters, a key feature for the maintenance of retinoid homeostatic balance in vertebrates, was only recently demonstrated in some mollusk species, suggesting that such ability is older than previously anticipated. In contrast, the enzymatic repertoire involved in this process is probably unlike that of vertebrates. The suggested ancestry of active retinoid pathways implies that many more metazoan species might be potential targets for endocrine disrupting chemicals. Here, we review the current knowledge about the occurrence and functionality of retinoid metabolic and signaling pathways in invertebrate lineages, paying special attention to the evolutionary origin of retinoid storage mechanisms. Additionally, we summarize existing information on the endocrine disruption of invertebrate retinoid modules by environmental chemicals. Research priorities in the field are highlighted.

Gene duplication and loss are powerful drivers of evolutionary change. The role of loss in phenotypic diversification is notably illustrated by the variable enzymatic repertoire involved in vertebrate protein digestion. Among these we find the pepsin family of aspartic proteinases, including chymosin (Cmy). Previous studies demonstrated that Cmy, a neo-natal digestive pepsin, is inactivated in some primates, including humans. This pseudogenization event was hypothesized to result from the acquisition of maternal immune immunoglobulin G (IgG) transfer. By investigating 94 mammalian subgenomes we reveal an unprecedented level of Cmy erosion in placental mammals, with numerous independent events of gene loss taking place in Primates, Dermoptera, Rodentia, Cetacea and Perissodactyla. Our findings strongly suggest that the recurrent inactivation of Cmy correlates with the evolution of the passive transfer of IgG and uncovers a noteworthy case of evolutionary cross-talk between the digestive and the immune system, modulated by gene loss. © 2017 Elsevier We acknowledge Fundação para a Ciência e a Tecnologia for the support to M.L.-M (SFRH/BD/84238/2012), to E.F. (SFRH/BD/79305/2011) and to R.R. (SFRH/BPD/72519/2010). LFCC research is supported by the MarInfo – Integrated Platform for Marine Data Acquisition and Analysis (reference NORTE-01-0145-FEDER-000031), a project supported by the North Portugal Regional Operational Program (NORTE 2020), under the PORTUGAL 2020 Partnership Agreement, through the European Regional Development Fund (ERDF) and by the European Regional Development Fund (ERDF) through the COMPETE - Operational Competitiveness Programme and POPH - Operational Human Potential Programme and national funds through FCT – Foundation for Science and Technology (UID/Multi/04423/2013). We would like to thank Liliana Gonçalves for the embryo illustration in the graphical abstract. Appendix A

Nuclear receptors (NRs) regulate numerous aspects of the endocrine system. They mediate endogenous and exogenous cues, ensuring a homeostatic control of development and metabolism. Gene duplication, loss and mutation have shaped the repertoire and function of NRs in metazoans. Here, we examine the evolution of a pivotal orchestrator of cholesterol metabolism in vertebrates, the liver X receptors (LXRs). Previous studies suggested that LXRα and LXRβ genes emerged in themammalian ancestor. However, we show through genomean alysis and functional assay that bona fide LXRα and LXRβ orthologues are present in reptiles, coelacanth and chondrichthyans but not in cyclostomes. These findings show that LXR duplicated before gnathostome radiation, followed by asymmetric paralogue loss in some lineages. We suggest that a tighter control of cholesterol levels in vertebrates was achieved through the exploitation of a wider range of oxysterols, an ability contingent on ligand-binding pocket remodeling. © 2016 The Author. We thank Neelakanteswar Aluru and Pedro Morais for the samples of little skate and amphioxus, respectively. This work was supported by Norte2020 and FEDER (Coral-Sustainable Ocean Exploitation-Norte-01-0145-FEDER- 000036), EXPL/MAR-EST/1540/2012. SFRH/BD/79305/2011 to E.F., SFRH/BD/84238/2012 to M.L.-M. and SFRH/BPD/72519/2010 to R.R. L.F.C.C. and M.M.S. are supported by the European Regional Development Fund (ERDF) through the COMPETE-Operational Competitiveness Programme and POPH-Operational Human Potential Programme and national funds through FCT-Foundation for Science and Technology (UID/Multi/04423/2013). This work was supported by the Branchiostoma lanceolatum genome consortium that provided access to the Branchiostoma lanceolatum genome sequence. We thank three anonymous reviewers for their comments and suggestions.