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.

A multi-tissue de novo transcriptome assembly and relative gene expression of the vulnerable freshwater salmonid Thymallus ligericus.

Freshwater ecosystems are among the most endangered ecosystems worldwide. While numerous taxa are on the verge of extinction as a result of global changes and direct or indirect anthropogenic activity, genomic and transcriptomic resources represent a key tool for comprehending species' adaptability and serve as the foundation for conservation initiatives. The Loire grayling, Thymallus ligericus, is a freshwater European salmonid endemic to the upper Loire River basin. The species is comprised of fragmented populations that are dispersed over a small area and it has been identified as a vulnerable species. Here, we provide a multi-tissue de novo transcriptome assembly of T. ligericus. The completeness and integrity of the transcriptome were assessed before and after redundancy removal with lineage-specific libraries from Eukaryota, Metazoa, Vertebrata, and Actinopterygii. Relative gene expression was assessed for each of the analyzed tissues, using the de novo assembled transcriptome and a genome-based analysis using the available T. thymallus genome as a reference. The final assembly, with a contig N50 of 1221 and Benchmarking Universal Single-Copy Orthologs (BUSCO) scores above 94%, is made accessible along with structural and functional annotations and relative gene expression of the five tissues (NCBI SRA and FigShare databases). This is the first transcriptomic resource for this species, which provides a foundation for future research on this and other salmonid species that are increasingly exposed to environmental stressors.

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.

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.

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.

A Highly Complex, MHC-Linked, 350 Million-Year-Old Shark Nonclassical Class I Lineage.

Cartilaginous fish, or Chondrichthyes, are the oldest extant vertebrates to possess the MHC and the Ig superfamily-based Ag receptors, the defining genes of the gnathostome adaptive immune system. In this work, we have identified a novel MHC lineage, UEA, a complex multigene nonclassical class I family found in sharks (division Selachii) but not detected in chimaeras (subclass Holocephali) or rays (division Batoidea). This new lineage is distantly related to the previously reported nonclassical class I lineage UCA, which appears to be present only in dogfish sharks (order Squaliformes). UEA lacks conservation of the nine invariant residues in the peptide (ligand)-binding regions (PBR) that bind to the N and C termini of bound peptide in most vertebrate classical class I proteins, which are replaced by relatively hydrophobic residues compared with the classical UAA. In fact, UEA and UCA proteins have the most hydrophobic-predicted PBR of all identified chondrichthyan class I molecules. UEA genes detected in the whale shark and bamboo shark genome projects are MHC linked. Consistent with UEA comprising a very large gene family, we detected weak expression in different tissues of the nurse shark via Northern blotting and RNA sequencing. UEA genes fall into three sublineages with unique characteristics in the PBR. UEA shares structural and genetic features with certain nonclassical class I genes in other vertebrates, such as the highly complex XNC nonclassical class I genes in Xenopus, and we anticipate that each shark gene, or at least each sublineage, will have a unique function, perhaps in bacterial defense.

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.

PseudoChecker: an integrated online platform for gene inactivation inference.

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.

Convergent Cortistatin losses parallel modifications in circadian rhythmicity and energy homeostasis in Cetacea and other mammalian lineages.

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.

A network-based approach to identify protein kinases critical for regulating srebf1 in lipid deposition causing obesity.

Obesity is a rapidly growing health pandemic, underlying a wide variety of disease conditions leading to increases in global mortality. It is known that the phosphorylation of various proteins regulates sterol regulatory element-binding transcription factors 1 (srebf1), a key lipogenic transcription factor, to cause the development of obesity. To detect the key protein kinases for regulating srebf1 in lipid deposition, we established the srebf1 knockout model in zebrafish (KO, srebf1-/-) by CRISPR/Cas9. The KO zebrafish exhibited a significant reduction of total free fatty acid content (fell 60.5%) and lipid deposition decrease compared with wild-type (WT) zebrafish. Meanwhile, srebf1 deletion in zebrafish eliminated lipid deposition induced by high-fat diet feeding. Compared with WT zebrafish, a total of 697 differentially expressed proteins and 316 differentially expressed phosphoproteins with 439 sites were identified in KO by differential proteomic and phosphoproteomic analyses. A significant number of proteins identified were involved in lipid and glucose metabolism. Moreover, some protein kinases critical for regulating srebf1 in lipid deposition, including Cdk2, Pkc, Prkceb, mTORC1, Mapk12, and Wnk1, were determined by network analyses. An in vitro study was performed to verify the network analysis results. Our findings provide potential targets (kinases) for human obesity treatments.

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.

Complete Inactivation of Sebum-Producing Genes Parallels the Loss of Sebaceous Glands in Cetacea.

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.

Cartilaginous fishes offer unique insights into the evolution of the nuclear receptor gene repertoire in gnathostomes.

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 evolution of S100A7: an unusual gene expansion in Myotis bats.

BACKGROUND: The S100A7 gene, also called psoriasin, was first described as an upregulated protein in psoriatic skin. For the past years, the importance of this protein as a key effector of innate immunity has been clearly established, not only due to its importance protecting against bacteria skin insult in humans, but also because of its important role in amplifying inflammatory processes. Given the importance of S100A7 in host defense, S100A7 genes have been mostly studied in humans. Here we provide a detailed analysis of the evolution of the gene family encoding for the S100A7 protein in mammals. RESULTS: Examination of several mammalian genomes revealed an unexpected variation in the copy number of S100A7. Among the most representative mammalian groups, we report that multiple events of duplication, gene loss and high mutation rates are shaping the evolution of this gene family. An unexpected result comes from Myotis species (order Chiroptera), where we found an outstanding S100A7 gene radiation, resulting in more than 10 copies in M. lucifugus and 5 copies in M. brandtii. These findings suggest a unique adaptive road in these species and are suggestive of special role of this protein in their immune system. CONCLUSIONS: We found different evolutionary histories among different mammalian groups. Overall, our results suggest that this gene family is evolving under the birth-and-death model of evolution. To our knowledge, this work represents the first detailed analysis of phylogenetic relationships of S100A7 within mammals and therefore will pave the way to further clarify their unique function in the immune system.

Convergent inactivation of the skin-specific C-C motif chemokine ligand 27 in mammalian evolution.

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, and start and stop codon mutations independently occurred in Cetacea, Pholidota, Sirenia, Chiroptera, and Rodentia, totalizing 18 species. The diverse habitat settings and lifestyles 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.

Tributyltin Affects Retinoid X Receptor-Mediated Lipid Metabolism in the Marine Rotifer Brachionus koreanus.

To examine how tributyltin (TBT), a model obesogen, affects the lipid metabolism in the marine rotifer Brachionus koreanus, we carried out life-cycle studies and determined the in vitro and in silico interactions of retinoid X receptor (RXR) with TBT, the transcriptional levels of RXR and lipid metabolic genes, and the fatty acid content. The lethal concentration 10% (LC10) was determined to be 5.12 µg/L TBT, and negative effects on ecologically relevant end points (e.g., decreased lifespan and fecundity) were detected at 5 µg/L TBT. On the basis of these findings, subsequent experiments were conducted below 1 µg/L TBT, which did not show any negative effects on ecologically relevant end points in B. koreanus. Nile red staining analysis showed that after exposure to 1 µg/L TBT, B. koreanus stored neutral lipids and had significantly increased transcriptional levels of RXR and lipid metabolism-related genes compared to the control. However, the content of total fatty acids did not significantly change at any exposure level. In the single fatty acids profile, a significant increase in saturated fatty acids (SFAs) 14:0 and 20:0 was observed, but the contents of omega-3 and omega-6 fatty acids were significantly decreased. Also, a transactivation assay of TBT with RXR showed that TBT is an agonist of Bk-RXR with a similar fold-induction to the positive control. Taken together, these results demonstrate that TBT-modulated RXR signaling leads to increase in transcriptional levels of lipid metabolism-related genes and the synthesis of SFAs but decreases the content of polyunsaturated fatty acids (PUFAs). Our findings support a wider taxonomic scope of lipid perturbation due to xenobiotic exposure that occurs via NRs in aquatic animals.

Linking chemical exposure to lipid homeostasis: A municipal waste water treatment plant influent is obesogenic for zebrafish larvae.

Obesity, a risk factor for the development of type-2 diabetes, hypertension, cardiovascular disease, hepatic steatosis and some cancers, has been ranked in the top 10 health risk in the world by the World Health Organization. Despite the growing body of literature evidencing an association between the obesity epidemic and specific chemical exposure across a wide range of animal taxa, very few studies assessed the effects of chemical mixtures and environmental samples on lipid homeostasis. Additionally, the mode of action of several chemicals reported to alter lipid homeostasis is still poorly understood. Aiming to fill some of these gaps, we combined an in vivo assay with the model species zebrafish (Danio rerio) to screen lipid accumulation and evaluate expression changes of key genes involved in lipid homeostasis, alongside with an in vitro transactivation assay using human and zebrafish nuclear receptors, retinoid X receptor α and peroxisome proliferator-activated receptor γ. Zebrafish larvae were exposed from 4 th day post-fertilization until the end of the experiment (day 18), to six different treatments: experimental control, solvent control, tributyltin at 100 ng/L Sn and 200 ng/L Sn (positive control), and wastewater treatment plant influent at 1.25% and 2.5%. Exposure to tributyltin and to 2.5% influent led to a significant accumulation of lipids, with white adipose tissue deposits concentrating in the perivisceral area. The highest in vitro tested influent concentration (10%) was able to significantly transactivate the human heterodimer PPARγ/RXRα, thus suggesting the presence in the influent of HsPPARγ/RXRα agonists. Our results demonstrate, for the first time, the ability of complex environmental samples from a municipal waste water treatment plant influent to induce lipid accumulation in zebrafish larvae.

Evolutionary Plasticity in Detoxification Gene Modules: The Preservation and Loss of the Pregnane X Receptor in Chondrichthyes Lineages.

To appraise how evolutionary processes, such as gene duplication and loss, influence an organism's xenobiotic sensitivity is a critical question in toxicology. Of particular importance are gene families involved in the mediation of detoxification responses, such as members of the nuclear receptor subfamily 1 group I (NR1I), the pregnane X receptor (PXR), and the constitutive androstane receptor (CAR). While documented in multiple vertebrate genomes, PXR and CAR display an intriguing gene distribution. PXR is absent in birds and reptiles, while CAR shows a tetrapod-specific occurrence. More elusive is the presence of PXR and CAR gene orthologs in early branching and ecologically-important Chondrichthyes (chimaeras, sharks and rays). Therefore, we investigated various genome projects and use them to provide the first identification and functional characterization of a Chondrichthyan PXR from the chimaera elephant shark (Callorhinchus milii, Holocephali). Additionally, we substantiate the targeted PXR gene loss in Elasmobranchii (sharks and rays). Compared to other vertebrate groups, the chimaera PXR ortholog displays a diverse expression pattern (skin and gills) and a unique activation profile by classical xenobiotic ligands. Our findings provide insights into the molecular landscape of detoxification mechanisms and suggest lineage-specific adaptations in response to xenobiotics in gnathostome evolution.

Retention of fatty acyl desaturase 1 (fads1) in Elopomorpha and Cyclostomata provides novel insights into the evolution of long-chain polyunsaturated fatty acid biosynthesis in vertebrates.

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.

Cetacea are natural knockouts for IL20.

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.