Conservation and divergence of chemical defense system in the tunicate Oikopleura dioica revealed by genome wide response to two xenobiotics

Effects of petrogenic pollutants on North Atlantic and Arctic Calanus copepods: From molecular mechanisms to population impacts

Oil and gas industries in the Northern Atlantic Ocean have gradually moved closer to the Arctic areas, a process expected to be further facilitated by sea ice withdrawal caused by global warming. Copepods of the genus Calanus hold a key position in these cold-water food webs, providing an important energetic link between primary production and higher trophic levels. Due to their ecological importance, there is a concern about how accidental oil spills and produced water discharges may impact cold-water copepods. In this review, we summarize the current knowledge of the toxicity of petroleum on North Atlantic and Arctic Calanus copepods. We also review how recent development of high-quality transcriptomes from RNA-sequencing of copepods have identified genes regulating key biological processes, like molting, diapause and reproduction in Calanus copepods, to suggest linkages between exposure, molecular mechanisms and effects on higher levels of biological organization. We found that the available ecotoxicity threshold data for these copepods provide valuable information about their sensitivity to acute petrogenic exposures; however, there is still insufficient knowledge regarding underlying mechanisms of toxicity and the potential for long-term implications of relevance for copepod ecology and phenology. Copepod transcriptomics has expanded our understanding of how key biological processes are regulated in cold-water copepods. These advances can improve our understanding of how pollutants affect biological processes, and thus provide the basis for new knowledge frameworks spanning the effect continuum from molecular initiating events to adverse effects of regulatory relevance. Such efforts, guided by concepts such as adverse outcome pathways (AOPs), enable standardized and transparent characterization and evaluation of knowledge and identifies research gaps and priorities. This review suggests enhancing mechanistic understanding of exposure-effect relationships to better understand and link biomarker responses to adverse effects to improve risk assessments assessing ecological effects of pollutant mixtures, like crude oil, in Arctic areas.

Lufenuron treatment temporarily represses gene expression and affects the SUMO pathway in liver of Atlantic salmon

Lufenuron is a benzoylurea insecticide currently in use to combat sea lice infestation in salmon aquaculture in Chile. With pending approval in Norway, the aim of this work was to study the uptake and toxicity of lufenuron in liver tissue of Atlantic salmon. Juvenile salmon weighing 40 g were given a standard 7-day oral dose, and bioaccumulation and transcriptional responses in the liver were examined 1 day after the end-of-treatment (day 8) and after 1 week of elimination (day 14). Bioaccumulation levels of lufenuron were 29 ± 3 mg/kg at day 8 and 14 ± 1 mg/kg at day 14, indicating relatively rapid clearance. However, residues of lufenuron were still present in the liver after 513 days of depuration. The exposure gave a transient inhibition of transcription in the liver at day 8 (2437 significant DEGs, p-adj < .05), followed by a weaker compensatory response at day 14 (169 significant DEGs). Pathways associated with RNA metabolism such as the sumoylation pathway were most strongly affected at day 8, while the apelin pathway was most profoundly affected at day 14. In conclusion, this study shows that lufenuron easily bioaccumulates and that a standard 7-day oral dose induces a transient inhibition of transcription in liver of salmon.

Transcriptome responses in copepods Calanus finmarchicus, Calanus glacialis and Calanus hyperboreus exposed to phenanthrene and benzo[a]pyrene

Arctic and sub-arctic pelagic organisms can be exposed to effluents and spills from offshore petroleum-related activities and thus it is important to understand how they respond to crude oil related contaminants such as polycyclic aromatic hydrocarbons (PAHs). The copepod species Calanus finmarchicus, Calanus glacialis and Calanus hyperboreus represent key links in the arctic marine food web. We performed a transcriptome analysis of the three species exposed to phenanthrene (Phe) and benzo[a]pyrene (BaP) representing low and high molecular weight PAHs, respectively. Differential expression of several genes involved in many cellular pathways was observed after 72 h exposure to Phe (0.1 μM) and BaP (0.1 μM). In C. finmarchicus and C. glacialis, the exposure resulted in up-regulation of genes encoding enzymes in xenobiotic biotransformation, particularly the phase II cytosolic sulfonation system that include 3'-phosphoadenosine 5'-phosphosulfate synthase (PAPSS) and sulfotransferases (SULTs). The sulfonation pathway genes were more strongly induced by BaP than Phe in C. finmarchicus and C. glacialis but were not affected in C. hyperboreus. However, a larger number of genes and pathways were modulated in C. hyperboreus by the PAHs including genes encoding xenobiotic biotransformation and lipid metabolism enzymes, suggesting stronger responses in this species. The results suggest that the cytosolic sulfonation is a major phase II conjugation pathway for PAHs in C. finmarchicus and C. glacialis. Some of the biotransformation systems affected are known to be involved in metabolism of endogenous compounds such as ecdysteroids, which may suggest potential interference with physiological and developmental processes of the copepod species.

The chemical defensome of five model teleost fish

How an organism copes with chemicals is largely determined by the genes and proteins that collectively function to defend against, detoxify and eliminate chemical stressors. This integrative network includes receptors and transcription factors, biotransformation enzymes, transporters, antioxidants, and metal- and heat-responsive genes, and is collectively known as the chemical defensome. Teleost fish is the largest group of vertebrate species and can provide valuable insights into the evolution and functional diversity of defensome genes. We have previously shown that the xenosensing pregnane x receptor (pxr, nr1i2) is lost in many teleost species, including Atlantic cod (Gadus morhua) and three-spined stickleback (Gasterosteus aculeatus), but it is not known if compensatory mechanisms or signaling pathways have evolved in its absence. In this study, we compared the genes comprising the chemical defensome of five fish species that span the teleosteii evolutionary branch often used as model species in toxicological studies and environmental monitoring programs: zebrafish (Danio rerio), medaka (Oryzias latipes), Atlantic killifish (Fundulus heteroclitus), Atlantic cod, and three-spined stickleback. Genome mining revealed evolved differences in the number and composition of defensome genes that can have implication for how these species sense and respond to environmental pollutants, but we did not observe any candidates of compensatory mechanisms or pathways in cod and stickleback in the absence of pxr. The results indicate that knowledge regarding the diversity and function of the defensome will be important for toxicological testing and risk assessment studies.

Cancer-related proteins in serum are altered in workers occupationally exposed to polycyclic aromatic hydrocarbons: a cross-sectional study

Exposure to some polycyclic aromatic hydrocarbons (PAH) increases the risk of cancer and is common particularly for workers in occupations such as chimney sweeping. In exposed workers, screening of early cancer-related markers provides important information to identify individuals at risk. Here, we aimed to elucidate the associations between PAH exposure and serum levels of cancer-related proteins in 118 chimney sweeps and 126 occupationally unexposed controls, all non-smoking males from Sweden. Monoydroxylated metabolites of pyrene, phenanthrene, benzo[a]pyrene and benzo[a]anthracene were measured in urine using liquid chromatography coupled to tandem mass spectrometry and 90 cancer-related proteins were measured in serum using a proximity extension assay. Linear regression analysis adjusted for age and body mass index, and false discovery rate (FDR) identified 17 serum proteins that were differentially expressed (16 upregulated and 1 downregulated) in chimney sweeps compared with controls (FDR < 0.05). Concentrations of the peptidase kallikrein 13 (KLK13) showed significant positive associations with urinary concentrations of the PAH metabolites 3-hydroxybenzo[a]pyrene (3-OH-BaP) [B, 95% confidence interval (CI): 0.042, 0.008–0.076] and 3-hydroxybenzo[a]anthracene (3-OH-BaA) (B, 95% CI: 0.068, 0.002–0.134). Moreover, dose–response relationships were observed between KLK13 and 3-OH-BaP (trend test P = 0.027) and 3-OH-BaA (P = 0.035). Pathway and gene ontology analyses showed that cell movement, cell adhesion and cell migration were the predominant molecular functions associated with the top differentially expressed proteins. In conclusion, we found a number of putative cancer-related proteins differentially expressed in workers exposed to PAH. This warrants effective measure to reduce PAH exposure among workers as well as further investigation to confirm these findings.

Two novel CYP3A isoforms in marine mussel Mytilus coruscus: Identification and response to cadmium and benzo[a]pyrene

CYP3A enzymes play a crucial role in metabolic clearance of a variety of xenobiotics. However, their genetic information and function remain unclear in molluscs. In the present study, two novel CYP3A genes i.e. McCYP3A-1 and McCYP3A-2 were identified and characterized from the thick shell mussel Mytilus coruscus, and their tissue distribution as well as the response to cadmium (Cd) and benzo[a]pyrene (B[α]P) exposure were addressed using real time quantitative RT-PCR (qRT-PCR) and erythromycin N-demethylase (ERND) assay. McCYP3A-1 and McCYP3A-2 possess typically domains of CYP family such as helix-C, helix-I, helix-K, PERF and the heme binding domain as well as the characteristic domains of CYP3s including six SRS motifs. McCYP3A-1 and McCYP3A-2 transcripts were constitutively expressed in all examined tissues with high expression level in digestive glands, hepatopancreas and gonads. Upon B[α]P exposure, McCYP3A-1 and McCYP3A-2 mRNA expression in digestive glands showed a pattern of up-regulation followed by down-regulation, while under Cd exposure, showed a time-dependent induction profile. In addition, ERND activity, generally used as an indicator of CYP3, increased in a time-dependent manner after exposure to Cd and B[α]P. These results collectively indicated that McCYP3A-1 and McCYP3A-2 are CYP3A family member and may play a potential role in metabolic clearance of xenobiotics. Meanwhile, the current results may provide some baseline data to support McCYP3A-1 and McCYP3A-2 as candidate biomarkers for monitoring of PAHs and heavy metal pollution.

The UDP-Glycosyltransferase (UGT) Superfamily: New Members, New Functions, and Novel Paradigms

UDP-glycosyltransferases (UGTs) catalyze the covalent addition of sugars to a broad range of lipophilic molecules. This biotransformation plays a critical role in elimination of a broad range of exogenous chemicals and by-products of endogenous metabolism, and also controls the levels and distribution of many endogenous signaling molecules. In mammals, the superfamily comprises four families: UGT1, UGT2, UGT3, and UGT8. UGT1 and UGT2 enzymes have important roles in pharmacology and toxicology including contributing to interindividual differences in drug disposition as well as to cancer risk. These UGTs are highly expressed in organs of detoxification (e.g., liver, kidney, intestine) and can be induced by pathways that sense demand for detoxification and for modulation of endobiotic signaling molecules. The functions of the UGT3 and UGT8 family enzymes have only been characterized relatively recently; these enzymes show different UDP-sugar preferences to that of UGT1 and UGT2 enzymes, and to date, their contributions to drug metabolism appear to be relatively minor. This review summarizes and provides critical analysis of the current state of research into all four families of UGT enzymes. Key areas discussed include the roles of UGTs in drug metabolism, cancer risk, and regulation of signaling, as well as the transcriptional and posttranscriptional control of UGT expression and function. The latter part of this review provides an in-depth analysis of the known and predicted functions of UGT3 and UGT8 enzymes, focused on their likely roles in modulation of levels of endogenous signaling pathways.

Oikopleura dioica: An Emergent Chordate Model to Study the Impact of Gene Loss on the Evolution of the Mechanisms of Development

The urochordate Oikopleura dioica is emerging as a nonclassical animal model in the field of evolutionary developmental biology (a.k.a. evo-devo) especially attractive for investigating the impact of gene loss on the evolution of mechanisms of development. This is because this organism fulfills the requirements of an animal model (i.e., has a simple and accessible morphology, a short generation time and life span, and affordable culture in the laboratory and amenable experimental manipulation), but also because O. dioica occupies a key phylogenetic position to understand the diversification and origin of our own phylum, the chordates. During its evolution, O. dioica genome has suffered a drastic process of compaction, becoming the smallest known chordate genome, a process that has been accompanied by exacerbating amount of gene losses. Interestingly, however, despite the extensive gene losses, including entire regulatory pathways essential for the embryonic development of other chordates, O. dioica retains the typical chordate body plan. This unexpected situation led to the formulation of the so-called inverse paradox of evo-devo, that is, when a genetic diversity is able to maintain a phenotypic unity. This chapter reviews the biological features of O. dioica as a model animal, along with the current data on the evolution of its genes and genome. We pay special attention to the numerous examples of gene losses that have taken place during the evolution of this unique animal model, which is helping us to understand to which the limits of evo-devo can be pushed off.

Transcriptome profiling suggests roles of innate immunity and digestion metabolism in purplish Washington clam

The purplish Washington clam (Saxidomus purpuratus) in the family Veneridae is distributed widely along the intertidal zones of northeast Asia and is increasingly being utilized as a commercially important food resource. Bivalves maintain homeostasis by regulating their food intake and digestion, innate immunity, and biotransformation in a mollusk-specific organ, the digestive gland. To understand digestive gland-specific pathways, we generated a high-quality de novo assembly of the digestive gland transcriptome of this clam using the Illumina Miseq platform. A total of 9.9 million raw reads were obtained and assembled using the Oases assembly platform, resulting in 27,358 contigs with an N50 of 433 bp. Functional gene annotations were performed using Gene Ontology, Eukaryotic Orthologous Groups, and Kyoto Encyclopedia of Genes and Genomes pathway analyses. In the transcriptome, many crucial genes involved in innate immunity and digestion metabolism were detected. A number of enzymes associated with drug metabolism were annotated, as much as that identified from the whole transcriptome of the Pacific oyster Crassostrea gigas. We provide valuable sequence information of S. purpuratus to predict functional understandings of the bivalve-specific digestive gland. This resource will be valuable for researchers comparing gene compositions and their expression levels in the digestive glands of bivalves.

Response to metals treatment of Fra1, a member of the AP-1 transcription factor family, in P. lividus sea urchin embryos

Lithium (Li), Nickel (Ni), and Zinc (Zn) are metals normally present in the seawater, although they can have adverse effects on the marine ecosystem at high concentrations by interfering with many biological processes. These metals are toxic for sea urchin embryos, affecting their morphology and developmental pathways. In particular, they perturb differently the correct organization of the embryonic axes (animal-vegetal, dorso-ventral): Li is a vegetalizing agent and Ni disrupts the dorso-ventral axis, while Zn has an animalizing effect. To deeply address the response of Paracentrotus lividus embryos to these metals, we studied the expression profiling of Pl-Fra transcription factor (TF), relating it to Pl-jun, a potential partner for AP-1 complex formation, and to Pl-MT, known to be an AP-1 target and to have a protective role against heavy metals. The AP-1 TFs are found throughout the animal kingdom and are involved in many cellular events, i.e. cell proliferation and differentiation, immune and stress responses, cancer growth. Here, we isolated the complete Pl-Fra cDNA and showed that Pl-Fra transcript, already present in the unfertilized eggs, was newly synthesized from the blastula stage, while its spatial distribution was mainly observed in skeletogenic cells, similarly to Pl-jun. Interestingly, Pl-Fra expression was induced by the different metals and the induction kinetics revealed its persistent expression during treatments. Moreover, its temporal and spatial behavior in response to the three metals was comparable to that of Pl-jun and Pl-MT. The understanding of AP-1 functions in invertebrates may provide new knowledge about the mechanisms of response to metal injuries, as well as it might lead to acknowledge the TFs as new type of biomarkers for the evaluation of hazards in polluted environment.

Increased fitness of a key appendicularian zooplankton species under warmer, acidified seawater conditions

Ocean warming and acidification (OA) may alter the fitness of species in marine pelagic ecosystems through community effects or direct physiological impacts. We used the zooplanktonic appendicularian, Oikopleura dioica, to assess temperature and pH effects at mesocosm and microcosm scales. In mesocosms, both OA and warming positively impacted O. dioica abundance over successive generations. In microcosms, the positive impact of OA, was observed to result from increased fecundity. In contrast, increased pH, observed for example during phytoplankton blooms, reduced fecundity. Oocyte fertility and juvenile development were equivalent under all pH conditions, indicating that the positive effect of lower pH on O. dioica abundance was principally due to increased egg number. This effect was influenced by food quantity and quality, supporting possible improved digestion and assimilation at lowered pH. Higher temperature resulted in more rapid growth, faster maturation and earlier reproduction. Thus, increased temperature and reduced pH had significant positive impacts on O. dioica fitness through increased fecundity and shortened generation time, suggesting that predicted future ocean conditions may favour this zooplankton species.

Metallothioneins of the urochordate Oikopleura dioica have Cys-rich tandem repeats, large size and cadmium-binding preference

Oikopleura dioica has the longest metallothionein described so far, made of repeats generated by a modular and step-wise evolution.

The choreography of the chemical defensome response to insecticide stress: insights into the Anopheles stephensi transcriptome using RNA-Seq

Animals respond to chemical stress with an array of gene families and pathways termed “chemical defensome”. In arthropods, despite many defensome genes have been detected, how their activation is arranged during toxic exposure remains poorly understood. Here, we sequenced the transcriptome of Anopheles stephensi larvae exposed for six, 24 and 48 hours to the LD₅₀ dose of the insecticide permethrin to monitor transcriptional changes of defensome genes across time. A total of 177 genes involved in insecticide defense were differentially expressed (DE) in at least one time-point, including genes encoding for Phase 0, I, II, III and antioxidant enzymes and for Heat Shock and Cuticular Proteins. Three major patterns emerged throughout time. First, most of DE genes were down-regulated at all time-points, suggesting a reallocation of energetic resources during insecticide stress. Second, single genes and clusters of genes turn off and on from six to 48 hours of treatment, showing a modulated response across time. Third, the number of up-regulated genes peaked at six hours and then decreased during exposure. Our results give a first picture of how defensome gene families respond against toxicants and provide a valuable resource for understanding how defensome genes work together during insecticide stress.

Rapid Evolutionary Rates and Unique Genomic Signatures Discovered in the First Reference Genome for the Southern Ocean Salp, Salpa thompsoni (Urochordata, Thaliacea)

A preliminary genome sequence has been assembled for the Southern Ocean salp, Salpa thompsoni (Urochordata, Thaliacea). Despite the ecological importance of this species in Antarctic pelagic food webs and its potential role as an indicator of changing Southern Ocean ecosystems in response to climate change, no genomic resources are available for S. thompsoni or any closely related urochordate species. Using a multiple-platform, multiple-individual approach, we have produced a 318,767,936-bp genome sequence, covering >50% of the estimated 602 Mb (±173 Mb) genome size for S. thompsoni. Using a nonredundant set of predicted proteins, >50% (16,823) of sequences showed significant homology to known proteins and ∼38% (12,151) of the total protein predictions were associated with Gene Ontology functional information. We have generated 109,958 SNP variant and 9,782 indel predictions for this species, serving as a resource for future phylogenomic and population genetic studies. Comparing the salp genome to available assemblies for four other urochordates, Botryllus schlosseri, Ciona intestinalis, Ciona savignyi and Oikopleura dioica, we found that S. thompsoni shares the previously estimated rapid rates of evolution for these species. High mutation rates are thus independent of genome size, suggesting that rates of evolution >1.5 times that observed for vertebrates are a broad taxonomic characteristic of urochordates. Tests for positive selection implemented in PAML revealed a small number of genes with sites undergoing rapid evolution, including genes involved in ribosome biogenesis and metabolic and immune process that may be reflective of both adaptation to polar, planktonic environments as well as the complex life history of the salps. Finally, we performed an initial survey of small RNAs, revealing the presence of known, conserved miRNAs, as well as novel miRNA genes; unique piRNAs; and mature miRNA signatures for varying developmental stages. Collectively, these resources provide a genomic foundation supporting S. thompsoni as a model species for further examination of the exceptional rates and patterns of genomic evolution shown by urochordates. Additionally, genomic data will allow for the development of molecular indicators of key life history events and processes and afford new understandings and predictions of impacts of climate change on this key species of Antarctic pelagic ecosystems.

Evolution by gene loss

The recent increase in genomic data is revealing an unexpected perspective of gene loss as a pervasive source of genetic variation that can cause adaptive phenotypic diversity. This novel perspective of gene loss is raising new fundamental questions. How relevant has gene loss been in the divergence of phyla? How do genes change from being essential to dispensable and finally to being lost? Is gene loss mostly neutral, or can it be an effective way of adaptation? These questions are addressed, and insights are discussed from genomic studies of gene loss in populations and their relevance in evolutionary biology and biomedicine.

Identification of a CYP3A-like gene and CYPs mRNA expression modulation following exposure to benzo[a]pyrene in the bivalve mollusk Chlamys farreri

In this study, we isolated a CYP3A-like gene from ovary of the scallop (Chlamys farreri). High levels of CYP3A-like gene expression occur in the digestive gland and gonad, which suggested their role in the metabolism of steroids and xenobiotics. Scallops were exposed to a polycyclic aromatic hydrocarbons (PAHs), benzo[a]pyrene (B[a]P) for 10 days. The CYP4 and CYP3A-like gene can be up-regulated by B[a]P in a dose-dependent manner after 10 days exposure. But no induction of the CYP3A-like was observed in 10 μg/L B[a]P group. The CYP1A-like expression can only be induced by 0.025 μg/L B[a]P. 0.5 and 10 μg/L B[a]P caused significant DNA damage and 10 μg/L B[a]P can also lead to oxidative damage. These results demonstrate that the mollusk CYPs can be modulated by environmental pollutant, and the blocked induction of CYP3A-like and CYP1A-like expression probably results from the high genotoxicity and oxidative damage partly.

Expression pattern of entire cytochrome P450 genes and response of defensomes in the benzo[a]pyrene-exposed monogonont rotifer Brachionus koreanus

Cytochrome P450 (CYP) proteins are involved in the first line of detoxification mechanism against diverse polycyclic aromatic hydrocarbons (PAHs) including benzo[a]pyrene (B[a]P). In aquatic invertebrates, there is still a lack of knowledge on the CYP genes involved in the molecular response to B[a]P exposure due to limited gene information. In this study, we cloned the entire 25 CYP genes in the monogonont rotifer Brachionus koreanus with the aid of next generation sequencing (NGS) technologies and analyzed their transcript profiles with a real-time RT-PCR array to better understand B[a]P-triggered molecular response over different time courses. As a result, B[a]P exposure induced CYP2/3-involved detoxification mechanisms and defensome, including phase II detoxification and antioxidant systems with a modulation of the chaperone heat shock protein (hsp) expression but did not change expression of other CYP clans in B. koreanus . Therefore, we found that B[a]P induced a strong detoxification mechanism to overcome detrimental effects of B[a]P associated with B[a]P-induced growth retardation as a trade-off in fitness costs. Also, this approach revealed that the entire CYP profiling can be a way of providing a better understanding on the mode of action of B[a]P in B. koreanus with respect to molecular defense metabolism.

Expression of three novel cytochrome P450 (CYP) and antioxidative genes from the polychaete, Perinereis nuntia exposed to water accommodated fraction (WAF) of Iranian crude oil and benzo[a]pyrene

To report a novel CYP genes and to evaluate its potency as a biomarker for oil pollution, we cloned three CYP genes and measured their expression profiles under controlled lab conditions using real-time reverse transcription PCR (real-time RT-PCR) after exposure of the water accommodated fraction (WAF) of Iranian crude oil and benzo[α]pyrene (B[α]P) as a positive control. Of these, CYP432A1 (CYP3 clan) gene was significantly induced by B[α]P exposure, indicating that the CYP3 clan gene would play an important role in polycyclic aromatic hydrocarbon (PAH) metabolisms, particularly for B[α]P in this species. However, the Perinereis nuntia CYP431A1 mRNA, a CYP2 clan gene, was sensitively expressed to WAF exposure with other two CYP genes. As one of Phase II detoxification enzymes, the glutathione S-transferase (GST) genes also upregulated with other antioxidant genes (SOD and CAT), indicating that WAF-exposed P. nuntia was properly responding to this kind of chemical stress. Thus, three CYP genes from the polychaete, P. nuntia have a potential as a biomarker in monitoring of the marine sediment after an oil spill accident.