A whole-body transcriptome assembly of the annelid worm Hediste diversicolor.

The Annelida phylum is composed of a myriad of species exhibiting key phenotypic adaptations. They occupy key ecological niches in a variety of marine, freshwater and terrestrial ecosystems. Importantly, the increment of omic resources is rapidly modifying the taxonomic landscape and knowledge of species belonging to this phylum. Here, we comprehensively characterised and annotated a transcriptome of the common ragworm, Hediste diversicolor (OF Müller). This species belongs to the family Nereididae and inhabits estuarine and lagoon areas on the Atlantic coasts of Europe and North America. Ecologically, H. diversicolor plays an important role in benthic food webs. Given its commercial value, H. diversicolor is a promising candidate for aquaculture development and production in farming facilities, under a circular economy framework. We used Illumina next-generation sequencing technology, to produce a total of 105 million (M) paired-end (PE) raw reads and generate the first whole-body transcriptome assembly of H. diversicolor species. This high-quality transcriptome contains 69,335 transcripts with an N50 transcript length of 2313 bp and achieved a BUSCO gene completeness of 97.7% and 96% in Eukaryota and Metazoa lineage-specific profile libraries. Our findings offer a valuable resource for multiple biological applications using this species.

Unravelling the main immune repertoire of Paracentrotus lividus following Vibrio anguillarum bath challenge.

Paracentrotus lividus is the most abundant echinoid species in the North East Atlantic Ocean and Mediterranean Sea. Although there is abundant genomic information of the species, there is no deep characterisation of the genes involved in the immune response. Here, a reference transcriptome of male and female coelomocytes was produced. The generated P. lividus transcriptome assembly has 203,511 transcripts, N50 transcript length of 1079 bp, and more than 90% estimated gene completeness in Eukaryota and Metazoa BUSCO databases, respectively. Differential gene expression analyses showed 54 and 55 up-regulated genes in P. lividus female and male coelomocyte tissues, respectively. These results suggest a similar immune gene repertoire between sexes. To examine the immune response, P. lividus was challenged with Vibrio anguillarum, one of the candidate pathogens for bald disease. Immune parameters were evaluated at cell and humoral levels, as well as the expression analysis of immune related genes at an early response stage. No differences were found at cellular and humoral levels with the exception of the increase of nitric oxide in perivisceral fluid of challenged animals. At the gene expression level, a total of 2721 genes were upregulated in challenged animals, 13.6 times higher expression than control group. Our analysis revealed that four major KEGG pathways were enriched in challenged animals: Autophagy (KEGG:04140), Endocytosis (KEGG:04144), Phagosome (KEGG:04145) and Protein processing in endoplasmic reticulum (KEGG:04141). Several toll-like receptors (TLR), scavenger receptors cysteine-rich (SRCR) or nucleotide-binding oligomerisation domain like receptors (NLR) were identified as major family genes for pathogen recognition and immune defence. This study provides a valuable transcriptomic resource and unfolds the molecular basis of immune response to V. anguillarum exposure. Overall, our findings contribute to the conservation effort of the P. lividus populations, as well as its sustainable exploitation in an aquaculture context.

Dataset of the complete mitogenome of the deep-sea sailfin roughshark, Oxynotus paradoxus Frade, 1929.

Chondrichthyans comprise a diverse group of vertebrate species with extraordinary ecological relevance. Yet, multiple members of this evolutionary lineage are associated with significant extinction risk. The sailfin roughshark Oxynotus paradoxus is a deep-water benthic shark currently listed as vulnerable due to population declines in parts of its range. Here we provide the first complete mitochondrial genome of O. paradoxus, comprising also the first record for the genus and family Oxynotidae. These data can facilitate future monitoring of the genetic diversity in this and related species. Genomic DNA was extracted from O. paradoxus collected in the eastern North Atlantic off western Portugal (37.59°N, 9.51°W) and sent for Illumina Paired-End (2 × 150 bp) library construction and whole genome sequencing on a Novaseq6000 platform. Trimmomatic (version 0.38) was used to remove adapters and MitoZ (version 3.4) to assemble and annotate the mitogenome. This mitogenome with 17 100 bp has a total of 38 genes, 13 of which are protein-coding genes, 23 transfer RNA genes, and 2 ribosomal RNA genes. Eight transfer RNAs and 1 protein-coding gene (NADH dehydrogenase subunit 6, NAD6) are in the complementary strand. In the provided phylogenetic inference, with all available and verified Squalomorphii mitogenomes, the four orders are well separated, and as expected, O. paradoxus is placed in the Squaliformes order. This data reinforces the need for more genomic resources for the Oxynotidae family.

An Ancestral Major Histocompatibility Complex Organization in Cartilaginous Fish: Reconstructing MHC Origin and Evolution.

Cartilaginous fish (sharks, rays, and chimeras) comprise the oldest living jawed vertebrates with a mammalian-like adaptive immune system based on immunoglobulins (Ig), T-cell receptors (TCRs), and the major histocompatibility complex (MHC). Here, we show that the cartilaginous fish "adaptive MHC" is highly regimented and compact, containing (i) a classical MHC class Ia (MHC-Ia) region containing antigen processing (antigen peptide transporters and immunoproteasome) and presenting (MHC-Ia) genes, (ii) an MHC class II (MHC-II) region (with alpha and beta genes) with linkage to beta-2-microglobulin (ß2m) and bromodomain-containing 2, (iii) nonclassical MHC class Ib (MHC-Ib) regions with 450 million-year-old lineages, and (iv) a complement C4 associated with the MHC-Ia region. No MHC-Ib genes were found outside of the elasmobranch MHC. Our data suggest that both MHC-I and MHC-II genes arose after the second round of whole-genome duplication (2R) on a human chromosome (huchr) 6 precursor. Further analysis of MHC paralogous regions across early branching taxa from all jawed vertebrate lineages revealed that Ig/TCR genes likely arose on a precursor of the huchr9/12/14 MHC paralog. The ß2m gene is linked to the Ig/TCR genes in some vertebrates suggesting that it was present at 1R, perhaps as the donor of C1 domain to the primordial MHC gene. In sum, extant cartilaginous fish exhibit a conserved and prototypical MHC genomic organization with features found in various vertebrates, reflecting the ancestral arrangement for the jawed vertebrates.

Mitochondrial replication's role in vertebrate mtDNA strand asymmetry.

Mitogenomes are defined as compact and structurally stable over aeons. This perception results from a vertebrate-centric vision, where few types of mtDNA rearrangements are described. Here, we bring a new light to the involvement of mitochondrial replication in the strand asymmetry of the vertebrate mtDNA. Using several species of deep-sea hatchetfish (Sternoptychidae) displaying distinct mtDNA structural arrangements, we unravel the inversion of the coding direction of protein-coding genes (PCGs). This unexpected change is coupled with a strand asymmetry nucleotide composition reversal and is shown to be directly related to the strand location of the Control Region (CR). An analysis of the fourfold redundant sites of the PCGs (greater than 6000 vertebrates), revealed the rarity of this phenomenon, found in nine fish species (five deep-sea hatchetfish). Curiously, in Antarctic notothenioid fishes (Trematominae), where a single PCG inversion (the only other record in fish) is coupled with the inversion of the CR, the standard asymmetry is disrupted for the remaining PCGs but not yet reversed, suggesting a transitory state. Our results hint that a relaxation of the classic vertebrate mitochondrial structural stasis promotes disruption of the natural balance of asymmetry of the mtDNA. These findings support the long-lasting hypothesis that replication is the main molecular mechanism promoting the strand-specific compositional bias of this unique and indispensable molecule.

Extensive gene loss parallels kidney aglomerulism in Syngnathidae.

The eccentric seahorses, seadragons, pipehorses and pipefishes (Syngnathidae) have an aglomerular kidney1. Here, we show that nephron genes2 conserved in Bilateria are secondarily eroded/deleted in Syngnathidae genomes. A transcriptome enrichment analysis suggests the predominance of excretion processes in the Syngnathidae kidney. In a lineage where crypsis and idleness are tightly associated, we propose that aglomerulism evolved as an energy-saving strategy.

Examination of gammarid transcriptomes reveals a widespread occurrence of key metabolic genes from epibiont bdelloid rotifers in freshwater species.

Previous data revealed the unexpected presence of genes encoding for long-chain polyunsaturated fatty acid (LC-PUFA) biosynthetic enzymes in transcriptomes from freshwater gammarids but not in marine species, even though closely related species were compared. This study aimed to clarify the origin and occurrence of selected LC-PUFA biosynthesis gene markers across all published gammarid transcriptomes. Through systematic searches, we confirmed the widespread occurrence of sequences from seven elongases and desaturases involved in LC-PUFA biosynthesis, in transcriptomes from freshwater gammarids but not marine species, and clarified that such occurrence is independent from the gammarid species and geographical origin. The phylogenetic analysis established that the retrieved elongase and desaturase sequences were closely related to bdelloid rotifers, confirming that multiple transcriptomes from freshwater gammarids contain contaminating rotifers' genetic material. Using the Adineta steineri genome, we investigated the genomic location and exon-intron organization of the elongase and desaturase genes, establishing they are all genome-anchored and, importantly, identifying instances of horizontal gene transfer. Finally, we provide compelling evidence demonstrating Bdelloidea desaturases and elongases enable these organisms to perform all the reactions for de novo biosynthesis of PUFA and, from them, LC-PUFA, an advantageous trait when considering the low abundance of these essential nutrients in freshwater environments.

Corrigendum: Uncovering a 500 million year old history and evidence of pseudogenization for TLR15.

[This corrects the article DOI: 10.3389/fimmu.2022.1020601.].

Decay of Skin-Specific Gene Modules in Pangolins.

The mammalian 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 gene Dgat2l6 seems independently eroded in the African and Asian clades. In contrast, Elovl3 implicated in the formation of specific neutral lipids required for skin barrier function is intact and expressed in the pangolin 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 transporter Abcc11, as previously described in Cetacea. Our findings reveal the sophisticated and complex history of gene retention and loss as key mechanisms in the evolution of the highly modified mammalian skin phenotypes.

The chromosome-level genome and key genes associated with mud-dwelling behavior and adaptations of hypoxia and noxious environments in loach (Misgurnus anguillicaudatus).

BACKGROUND: The loach (Misgurnus anguillicaudatus), the most widely distributed species of the family Cobitidae, displays a mud-dwelling behavior and intestinal air-breathing, inhabiting the muddy bottom of extensive freshwater habitats. However, lack of high-quality reference genome seriously limits the interpretation of the genetic basis of specialized adaptations of the loach to the adverse environments including but not limited to the extreme water temperature, hypoxic and noxious mud environment. RESULTS: This study generated a 1.10-Gb high-quality, chromosome-anchored genome assembly, with a contig N50 of 3.83 Mb. Multiple comparative genomic analyses found that proto-oncogene c-Fos (fos), a regulator of bone development, is positively selected in loach. Knockout of fos (ID: Mis0086400.1) led to severe osteopetrosis and movement difficulties, combined with the comparison results of bone mineral density, supporting the hypothesis that fos is associated with loach mud-dwelling behavior. Based on genomic and transcriptomic analysis, we identified two key elements involved in the intestinal air-breathing of loach: a novel gene (ID: mis0158000.1) and heat shock protein beta-1 (hspb1). The flavin-containing monooxygenase 5 (fmo5) genes, central to xenobiotic metabolism, undergone expansion in loach and were identified as differentially expressed genes in a drug stress trial. A fmo5-/- (ID: Mis0185930.1) loach displayed liver and intestine injury, indicating the importance of this gene to the adaptation of the loach to the noxious mud. CONCLUSIONS: Our work provides valuable insights into the genetic basis of biological adaptation to adverse environments.

The gill transcriptome of threatened European freshwater mussels.

Genomic tools applied to non-model organisms are critical to design successful conservation strategies of particularly threatened groups. Freshwater mussels of the Unionida order are among the most vulnerable taxa and yet almost no genetic resources are available. Here, we present the gill transcriptomes of five European freshwater mussels with high conservation concern: Margaritifera margaritifera, Unio crassus, Unio pictorum, Unio mancus and Unio delphinus. The final assemblies, with N50 values ranging from 1069-1895 bp and total BUSCO scores above 90% (Eukaryote and Metazoan databases), were structurally and functionally annotated, and made available. The transcriptomes here produced represent a valuable resource for future studies on these species' biology and ultimately guide their conservation.

The repertoire of the elongation of very long-chain fatty acids (Elovl) protein family is conserved in tambaqui (Colossoma macropomum): Gene expression profiles offer insights into the sexual differentiation process.

Elongation of very long-chain fatty acids (Elovl) proteins are critical players in the regulation of the length of a fatty acid. At present, eight members of the Elovl family (Elovl1-8), displaying a characteristic fatty acid substrate specificity, have been identified in vertebrates, including teleost fish. In general, Elovl1, Elovl3, Elovl6 and Elovl7 exhibit a substrate preference for saturated and monounsaturated fatty acids, while Elovl2, Elovl4, Elovl5 and Elovl8 use polyunsaturated fatty acids (PUFA) as substrates. PUFA elongases have received considerable attention in aquatic animals due to their involvement in the conversion of C18 PUFAs to long-chain polyunsaturated fatty acids (LC-PUFA). Here, we identified the full repertoire of elovl genes in the tambaqui Colossoma macropomum genome. A detailed phylogenetic and synteny analysis suggests a conservation of these genes among teleosts. Furthermore, based on RNAseq gene expression data, we discovered a gender bias expression of elovl genes during sex differentiation of tambaqui, toward future males. Our findings suggest a role of Elovl enzymes and fatty acid metabolism in tambaqui sexual differentiation.

A mitochondrial genome assembly of the opal chimaera, Chimaera opalescens Luchetti, Iglésias et Sellos 2011, using PacBio HiFi long reads.

Chondrichthyans (sharks, rays and chimeras) are a fascinating and highly vulnerable group of early branching gnathostomes. However, they remain comparatively poorly sampled from the point of view of molecular resources, with deep water taxa being particularly data deficient. The development of long-read sequencing technologies enables the analysis of phylogenetic relationships through a precise and reliable assembly of complete mtDNA genomes. The sequencing and characterization of the complete mitogenome of the opal chimera Chimera opalescens Luchetti, Iglésias et Sellos 2011, using the long-read technique PacBio HiFi is presented. The entire mitogenome was 23,411 bp long and shows the same overall content, i.e. 13 protein-coding genes, 22 transfer RNA and 2 ribosomal RNA genes, as all other examined Chondrichthyan mitogenomes. Phylogenetic reconstructions using all available Chondrichthyan mitogenomes, including 11 Holocephali (chimeras and ratfishes), places C. opalescens within the Chimaeridae family. Furthermore, the results reinforce previous findings, showing the genus Chimera as paraphyletic and thus highlighting the need to expand molecular approaches in this group of cartilaginous fishes.

A zebrafish pparγ gene deletion reveals a protein kinase network associated with defective lipid metabolism.

Peroxisome proliferator-activated receptor γ (Pparγ) is a master regulator of adipogenesis. Chronic pathologies such as obesity, cardiovascular diseases, and diabetes involve the dysfunction of this transcription factor. Here, we generated a zebrafish mutant in pparγ (KO) with CRISPR/Cas9 technology and revealed its regulatory network. We uncovered the hepatic phenotypes of these male and female KO, and then the male wild-type zebrafish (WT) and KO were fed with a high-fat (HF) or standard diet (SD). We next conducted an integrated analyze of the proteomics and phosphoproteomics profiles. Compared with WT, the KO showed remarkable hyalinization and congestion lesions in the liver of males. Strikingly, pparγ deletion protected against the influence of high-fat diet feeding on lipid deposition in zebrafish. Some protein kinases critical for lipid metabolism, including serine/threonine-protein kinase TOR (mTOR), ribosomal protein S6 kinase (Rps6kb1b), and mitogen-activated protein kinase 14A (Mapk14a), were identified to be highly phosphorylated in KO based on differential proteome and phosphoproteome analysis. Our study supplies a pparγ deletion animal model and provides a comprehensive description of pparγ-induced expression level alterations of proteins and their phosphorylation, which are vital to understand the defective lipid metabolism risks posed to human health.

Neuroendocrine pathways at risk? Simvastatin induces inter and transgenerational disruption in the keystone amphipod Gammarus locusta.

The primary focus of environmental toxicological studies is to address the direct effects of chemicals on exposed organisms (parental generation - F0), mostly overlooking effects on subsequent non-exposed generations (F1 and F2 - intergenerational and F3 transgenerational, respectively). Here, we addressed the effects of simvastatin (SIM), one of the most widely prescribed human pharmaceuticals for the primary treatment of hypercholesterolemia, using the keystone crustacean Gammarus locusta. We demonstrate that SIM, at environmentally relevant concentrations, has significant inter and transgenerational (F1 and F3) effects in key signaling pathways involved in crustaceans' neuroendocrine regulation (Ecdysteroids, Catecholamines, NO/cGMP/PKG, GABAergic and Cholinergic signaling pathways), concomitantly with changes in apical endpoints, such as depressed reproduction and growth. These findings are an essential step to improve hazard and risk assessment of biological active compounds, such as SIM, and highlight the importance of studying the transgenerational effects of environmental chemicals in animals' neuroendocrine regulation.

A genome assembly of the Atlantic chub mackerel (Scomber colias): a valuable teleost fishing resource.

The Atlantic chub mackerel, Scomber colias (Gmelin, 1789), is a medium-sized 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 Mbp, distributed into 2,028 scaffolds and 2,093 contigs with an N50 length of 4.19 and 3.34 Mbp, 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 genome assembly ranks fourth in high-quality genome assemblies within the order Scombriformes and first in the genus Scomber.

Uncovering a 500 million year old history and evidence of pseudogenization for TLR15.

Introduction: Toll like receptors (TLRs) are at the front line of pathogen recognition and host immune response. Many TLR genes have been described to date with some being found across metazoans while others are restricted to specific lineages. A cryptic member of the TLR gene family, TLR15, has a unique phylogenetic distribution. Initially described in extant species of birds and reptiles, an ortholog has been reported for cartilaginous fish. Methods: Here, we significantly expanded the evolutionary analysis of TLR15 gene evolution, taking advantage of large genomic and transcriptomic resources available from different lineages of vertebrates. Additionally, we objectively search for TLR15 in lobe-finned and ray-finned fish, as well as in cartilaginous fish and jawless vertebrates. Results and discussion: We confirm the presence of TLR15 in early branching jawed vertebrates - the cartilaginous fish, as well as in basal Sarcopterygii - in lungfish. However, within cartilaginous fish, the gene is present in Holocephalans (all three families) but not in Elasmobranchs (its sister-lineage). Holocephalans have long TLR15 protein sequences that disrupt the typical TLR structure, and some species display a pseudogene sequence due to the presence of frameshift mutations and early stop codons. Additionally, TLR15 has low expression levels in holocephalans when compared with other TLR genes. In turn, lungfish also have long TLR15 protein sequences but the protein structure is not compromised. Finally, TLR15 presents several sites under negative selection. Overall, these results suggest that TLR15 is an ancient TLR gene and is experiencing ongoing pseudogenization in early-branching vertebrates.

Convergent Loss of the Necroptosis Pathway in Disparate Mammalian Lineages Shapes Viruses Countermeasures.

Programmed cell death is a vital process in the life cycle of organisms. Necroptosis, an evolutionary form of programmed necrosis, contributes to the innate immune response by killing pathogen-infected cells. This virus-host interaction pathway is organized around two components: the receptor-interacting protein kinase 3 (RIPK3), which recruits and phosphorylates the mixed lineage kinase-like protein (MLKL), inducing cellular plasma membrane rupture and cell death. Critically, the presence of necroptotic inhibitors in viral genomes validates necroptosis as an important host defense mechanism. Here, we show, counterintuitively, that in different mammalian lineages, central components of necroptosis, such as RIPK3 and MLKL, are deleted or display inactivating mutations. Frameshifts or premature stop codons are observed in all the studied species of cetaceans and leporids. In carnivores' genomes, the MLKL gene is deleted, while in a small number of species from afrotheria and rodentia premature stop codons are observed in RIPK3 and/or MLKL. Interestingly, we also found a strong correlation between the disruption of necroptosis in leporids and cetaceans and the absence of the N-terminal domain of E3-like homologs (responsible for necroptosis inhibition) in their naturally infecting poxviruses. Overall, our study provides the first comprehensive picture of the molecular evolution of necroptosis in mammals. The loss of necroptosis multiple times during mammalian evolution highlights the importance of gene/pathway loss for species adaptation and suggests that necroptosis is not required for normal mammalian development. Moreover, this study highlights a co-evolutionary relationship between poxviruses and their hosts, emphasizing the role of host adaptation in shaping virus evolution.

Complete mitochondrial genome of the ragworm annelid Hediste diversicolor (of Müller, 1776) (Annelida: Nereididae).

Marine annelids are a globally distributed and species-rich group, performing important ecological roles in macrobenthic communities. Yet, the availability of molecular resources to study these organisms is scarcer, comparatively with other phyla. Here, we present the first complete mitogenome of the Atlantic ragworm Hediste diversicolor (OF Muller, 1776). The mitogenome (15,904 bp long) contains 13 protein-coding genes, 22 transfer RNA, and two ribosomal RNA genes, all encoded in the same strand. Gene arrangement and composition are identical to those observed in two available congeneric species, Hediste diadroma and Hediste japonica. The phylogenetic analysis using both maximum-likelihood and Bayesian inference methods reveal a well-supported monophyly of genus Hediste and the already reported paraphyletic relationships within the subfamilies Nereidinae and Gymnonereidinae. Our results highlight the relevance of increasing the molecular sampling within this diverse group of marine fauna.

Functional or Vestigial? The Genomics of the Pineal Gland in Xenarthra.

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 probably lost independently among Xenarthra orders. 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 endorses the use of genomic analyses to ascertain the physiological status of suspected vestigial structures.