Applications for next-generation sequencing in fish ecotoxicogenomics

Conjoint transcriptomics and metabolomics analyses provide insights into the toxicity of acetamiprid to Lethenteron reissneri larvae

The neonicotinoid pesticide acetamiprid has been widely used in agricultural pest control and was frequently detected in the water environment. There have been some studies of the toxic effects of acetamiprid on fish, but studies on aquatic lower vertebrates are still very limited. As a primitive jawless vertebrate, Lethenteron reissneri has a special position in evolution and is now listed as a national second level protected animal in China. The present study aimed to investigate the toxic effect of acetamiprid on the liver of L. reissneri larvae. A conjoint analysis of the transcriptomics and metabolomics was performed to determine the responses of L. reissneri larvae liver to acetamiprid at different concentrations (L for low concentration 25 mg/L and H for high concentration 100 mg/L). Even low concentrations of acetamiprid can cause significant liver damage to L. reissneri larvae in a short period. In omics analyses, 2141 differentially expressed genes (DEGs) and 183 differentially abundant metabolites (DAMs) were identified in the H/Control group, and 229 DEGs and 144 DAMs were identified in the L/C group. Correlation analyses revealed acetamiprid affected the metabolic pathways of L. reissneri larvae liver such as the glycerophospholipid metabolism and arachidonic acid metabolism. This study not only enriches the basis for understanding the toxic effect of acetamiprid exposure to L. reissneri larvae liver and provides more information on the breeding and conservation of L. reissneri, but also further causes attention on toxicity risk from acetamiprid to aquatic lower vertebrate species.

miR-124 mediates the expression of ccBax to regulate Cyprinid herpesvirus 2 (CyHV-2)-induced apoptosis and viral replication

Cyprinid herpesvirus 2 (CyHV-2), the etiological agent of herpesvirus haematopoietic necrosis (HVHN) in carp and goldfish, has caused significant economic losses in the aquaculture industry. During viral infection, the host initiates a series of active or passive defences to regulate the process of virus infection. Apoptosis is a key component of active cellular defence, and members of the Bcl-2 family have been shown to play a critical role in the apoptotic process. However, the mechanism of action of the Bcl-2 family in inducing apoptosis during CyHV-2 infection remains unclear. In this study, we revealed the molecular mechanism of miRNA-mediated silver crucian carp BAX (ccBax) in CyHV-2-induced apoptosis for the first time and demonstrated that the overexpression of miR-124 suppressed ccBax expression and significantly down-regulated apoptosis in caudal fin cells of Carassius auratus gibelio (GiCF), while miR-124 inhibitors were the opposite. These studies indicated that miR-124 inhibits CyHV-2-induced apoptosis by reducing the expression of ccBax. Furthermore, the fact that transfection of miR-124 mimics promoted CyHV-2 replication, whereas miR-124 inhibitors inhibited CyHV-2 replication, indicated that miR-124 inhibited CyHV-2-induced apoptosis and contributed to viral replication. All these results suggested that miR-124 suppresses virus-induced apoptosis and promotes viral replication by targeting and regulating ccBax expression.

Sex steroid profile during oocyte development and maturation in the intertidal worm Marphysa madrasi (Polychaeta: Eunicidae) from the east coast of India

Marphysa madrasi is a commercially valuable maturation diet in crustacean aquaculture. This study presents the first detailed investigation of oogenesis in the intertidal polychaete worm M. madrasi and reports the steroid profile during oocyte growth and development. Oogenesis is extraovarian type I, originating from coelomic epithelial cells, with four stages of development - primary growth, early vitellogenic, late vitellogenic, and maturation. The primary growth phase contains oogonial cells and previtellogenic oocyte clusters in the early, mid, and late stages of development form a dispersed ovary attached to blood vessels. The late previtellogenic oocytes detach from the ovary at the onset of vitellogenesis. The detached oocytes complete vitellogenesis and final maturation in the coelomic fluid as solitary free-floating cells without any connection with follicle cells. The worms display asynchronous reproduction with a heterogeneous population of developing oocytes. Steroid extracts from the polychaete homogenates in different stages of oogenesis were identified by HPLC and confirmed by LC-MS/MS. In M. madrasi, two vertebrate-type steroids, pregnenolone (P5) and 17α-hydroxyprogesterone (17-OHP) were detected and quantified. The P5 levels were low in immature worms but increased significantly by ∼ 8.3-fold in the previtellogenic stage and peaked during oocyte maturation. 17-OHP levels were low in immature worms but gradually increase as the oogenesis progress to the primary growth and early vitellogenic phase, with a significant increase (p < 0.001) during the late vitellogenic phase. Although an increase in the concentration of P5 and 17-OHP during vitellogenesis and maturation of oocytes points to a possible role in reproduction, the absence of other vertebrate-type steroids in the investigated polychaete signifies a plausible uptake of P5 and 17-OHP from the environment.

10 Years of Toxicogenomics section in Frontiers in Genetics: Past discoveries and Future Perspectives

The Frontiers Media family has over 200 journals, which are each headed by usually one Field Chief Editor, and several thousand specialty sections, which are each headed by one or more Specialty Chief Editors. The year 2021 was the 10th anniversary of the founding of the Frontiers in Genetics journal and the Frontiers in Toxicogenomics specialty section of this journal. In 2021, we also announce one of the newest of the Frontiers journals-Frontiers in Toxicology which is part of the Frontiers Media family of journals but independent of Frontiers in Genetics. Dr. Ruden is the founding, and currently sole, Specialty Chief Editor of Frontiers in Toxicogenomics and one of 9 Specialty Chief Editors of Frontiers in Toxicology. As of 2021, Frontiers in Toxicogenomics has published over 138 articles and has over 370 Editors including 90 Associate Editors and 280 Review Editors. The Frontiers in Genetics impact factor was initially approximately 2.5 when it was first listed in PubMed in 2015 and has risen steadily to its current value of 4.8, which is typical for the majority of the over 200 Frontiers journals that have established impact factors. In this overview of the first decade of Frontiers in Toxicogenomics, we discuss the top 5 articles with the highest Scopus citations, which were all written in the first few years of the journal. The article with the highest number of citations, with 353 Scopus over 600 Google Scholar citations, and the highest average number of citations (67) that steadily increased from 10 citations in 2013 to 119 citations in 2021, was written in 2012 by Dr. Ruden's laboratory and titled, "Using Drosophila melanogaster as a model for genotoxic chemical mutational studies with a new program, SnpSift." The five most influential authors who published in the journal in the past 10 years based on Scopus citations of a particular paper are Dr. Ruden's laboratory, with 353 Scopus citations for the SnpSift paper mentioned above; Drs. Brock Christensen and Carmen J. Marsit, with 86 Scopus citations for their review, "Epigenomics in environmental health"; Dr. Michael Aschner and colleagues, with 61 Scopus citations for their paper "Genetic factors and manganese-induced neurotoxicity"; and Dr. Sandra C. dos Santos and colleagues, with 59 Scopus citations for their paper, "Yeast toxicogenomics: genome-wide responses to chemical stresses with impact in environmental health, pharmacology, and biotechnology." While the top 5 papers were published in the early years of the journal, we will also discuss a more recent article published in 2018 on a comparison of RNA-seq and microarray methods by Dr. Michael Liguori's laboratory, "Comparison of RNA-Seq and Microarray Gene Expression Platforms for the Toxicogenomic Evaluation of Liver From Short-Term Rat Toxicity Studies," that far exceeds the number of downloads and views of all the other articles published in the first 10 years of the journal and will likely be a top cited paper in the second decade highlights of this journal. Finally, we discuss where the Frontiers in Toxicogenomics specialty journal and the Frontiers in Toxicology journal will go to advance the field of toxicogenomics, and more generally, toxicology, in the future.

Gene expression changes in Epinephelus marginatus (Teleostei, Serranidae) liver reveals candidate molecular biomarker of iron ore contamination

Wastes from iron ore mining activities are potentially damaging to adjacent aquatic ecosystems. We aimed to determine biomarkers of environmental exposure to this xenobiotic in the dusky grouper Epinephelus marginatus by differential gene expression analysis. For this, fish were exposed to iron ore (15.2 mg/L) and gene expression in liver was assessed by RNA-Seq and compared to the control group. A total of 124 differentially expressed genes were identified, from which 52 were upregulated and 72 were downregulated in response to iron ore. From these, ferritin (medium subunit), cytochrome b reductase and epoxide hydrolase genes were selected for validation by RT-qPCR that confirmed the upregulation of epoxide hydrolase in fish exposed to iron ore.

Transcriptomic and Histological Analysis of the Greentail Prawn (Metapenaeus bennettae) Following Light Crude Oil Exposure

Oil spills pose a significant threat to marine biodiversity. Crude oil can partition into sediments where it may be persistent, placing benthic species such as decapods at particular risk of exposure. Transcriptomic and histological tools are often used to investigate the effects of hydrocarbon exposure on marine organisms following oil spill events, allowing for the identification of metabolic pathways impacted by oil exposure. However, there is limited information available for decapod crustaceans, many of which carry significant economic value. In the present study, we assess the sublethal impacts of crude oil exposure in the commercially important Australian greentail prawn (Metapenaeus bennettae) using transcriptomic and histological analyses. Prawns exposed to light, unweathered crude oil "spiked" sediments for 90 h were transferred to clean sediments for a further 72 h to assess recovery. Chemical analyses indicated that polycyclic aromatic hydrocarbons increased by approximately 65% and 91% in prawn muscle following 24 and 90 h of exposure, respectively, and significantly decreased during 24- and 72-h recovery periods. Transcriptomic responses followed an exposure and recovery pattern with innate immunity and nutrient metabolism transcripts significantly lowered in abundance after 24 h of exposure and were higher in abundance after 72 h of recovery. In addition, transcription/translation, cellular responses, and DNA repair pathways were significantly impacted after 24 h of exposure and recovered after 72 h of recovery. However, histological alterations such as tubule atrophy indicated an increase in severity after 24 and 72 h of recovery. The present study provides new insights into the sublethal impacts of crude oil exposure in greentail prawns and identifies molecular pathways altered by exposure. We expect these findings to inform future management associated with oil extraction activity and spills. Environ Toxicol Chem 2022;41:2162-2180. © 2022 John Wiley & Sons Ltd. This article has been contributed to by U.S. Government employees and their work is in the public domain in the USA.

Comparative evaluation of fish larval preservation methods on microbiome profiles to aid in metagenomics research

Applications of microbiome research through metagenomics promise to generate microbiome manipulation strategies for improved larval survival in aquaculture. However, existing lacunae on the effects of sample preservation methods in metagenome profiles hinder the successful application of this technique. In this context, four preservation methods were scrutinized to identify reliable methods for fish larval microbiome research. The results showed that a total of ten metagenomics metrics, including DNA yield, taxonomic and functional microbiome profiles, and diversity measures, were significantly (P < 0.05) influenced by the preservation method. Activity ranking based on the performance and reproducibility showed that three methods, namely immediate direct freezing, room temperature preservation in absolute ethanol, and preservation at - 20 °C in lysis, storage, and transportation buffer, could be recommended for larval microbiome research. Furthermore, as there was an apparent deviation of the microbiome profiles of ethanol preserved samples at room temperature, the other methods are preferred. Detailed analysis showed that this deviation was due to the bias towards Vibrionales and Rhodobacterales. The microbial taxa responsible for the dissimilarity across different methods were identified. Altogether, the paper sheds light on the preservation protocols of fish larval microbiome research for the first time. The results can help in cross-comparison of future and past larval microbiome studies. Furthermore, this is the first report on the activity ranking of preservation methods based on metagenomics metrics. Apart from methodological perspectives, the paper provides for the first time certain insights into larval microbial profiles of Rachycentron canadum, a potential marine aquaculture species. KEY POINTS: • First report on effects of preservation methods on fish larval microbiome profiles. • First report on activity ranking of preservation methods based on metagenomics metrics. • Storage methods influenced DNA yield, taxonomic and functional microbiome profiles.

MicroRNA-124 Promotes Singapore Grouper Iridovirus Replication and Negatively Regulates Innate Immune Response

Viral infections seriously affect the health of organisms including humans. Now, more and more researchers believe that microRNAs (miRNAs), one of the members of the non-coding RNA family, play significant roles in cell biological function, disease occurrence, and immunotherapy. However, the roles of miRNAs in virus infection (entry and replication) and cellular immune response remain poorly understood, especially in low vertebrate fish. In this study, based on the established virus-cell infection model, Singapore grouper iridovirus (SGIV)-infected cells were used to explore the roles of miR-124 of Epinephelus coioides, an economically mariculture fish in southern China and Southeast Asia, in viral infection and host immune responses. The expression level of E. coioides miR-124 was significantly upregulated after SGIV infection; miR-124 cannot significantly affect the entry of SGIV, but the upregulated miR-124 could significantly promote the SGIV-induced cytopathic effects (CPEs), the viral titer, and the expressions of viral genes. The target genes of miR-124 were JNK3/p38α mitogen-activated protein kinase (MAPK). Overexpression of miR-124 could dramatically inhibit the activation of NF-κB/activating protein-1 (AP-1), the transcription of proinflammatory factors, caspase-9/3, and the cell apoptosis. And opposite results happen when the expression of miR-124 was inhibited. The results suggest that E. coioides miR-124 could promote viral replication and negatively regulate host immune response by targeting JNK3/p38α MAPK, which furthers our understanding of virus and host immune interactions.

Immunotoxicological effects of insecticides in exposed fishes

Biologically active compounds used in agriculture that develop near aquatic environments easily spill into rivers or lakes. As a result, insecticides, herbicides and fungicides are observed worldwide in aquatic environments and accumulated in aquatic organism. Many insecticides, including organochlorine and organophosphate, have long been banned long ago because of their high persistence and non-target toxicity. However, previous studies have shown that persistent pesticides remain in aquatic organisms. The immune system is the first defense mechanism against exposure to persistent organic pollutants or pesticides that have been released into the aquatic environment. Many insecticides have been reported to cause immunotoxicity, which is represented by alteration of phagocytic and lysozyme activity. Recent studies show that immunotoxicity by insecticides exerts a more complex mechanism in fish. Insecticides induce immunotoxic effects, such as the release of inflammatory cytokines from head kidney macrophages and inhibition of immune cell proliferation in fish, which can lead to death in severe cases. Even currently used pesticides, such as pyrethroid, with low bioaccumulation have been shown to induce immunotoxicological effects in fish when exposed continuously. Therefore, this review describes the types and bioaccumulation of insecticides that cause immunotoxicity and detailed immunotoxicological mechanisms in fish tissues.

Gill and liver transcriptomic responses of Achirus lineatus (Neopterygii: Achiridae) exposed to water-accommodated fraction (WAF) of light crude oil reveal an onset of hypoxia-like condition

Crude oil is one of the most widespread pollutants released into the marine environment, and native species have provided useful information about the effect of crude oil pollution in marine ecosystems. We consider that the lined sole Achirus lineatus can be a useful monitor of the effect of crude oil in the Gulf of Mexico (GoM) because this flounder species has a wide distribution along the GoM, and its response to oil components is relevant. The objective of this study was to compare the transcriptomic changes in liver and gill of adults lined sole fish (Achirus lineatus) exposed to a sublethal acute concentration of water-accommodated fraction (WAF) of light crude oil for 48 h. RNA-Seq was performed to assess the transcriptional changes in both organs. A total of 1073 differentially expressed genes (DEGs) were detected in gills; 662 (61.69%) were upregulated, and 411 (38.30%) were downregulated whereas in liver, 515 DEGs; 306 (59.42%) were upregulated, and 209 (40.58%) were downregulated. Xenobiotic metabolism and redox metabolism, along with DNA repair mechanisms, were activated. The induction of hypoxia-regulated genes and the generalized regulation of multiple signaling pathways support the hypothesis that WAF exposition causes a hypoxia-like condition.

Ciguatera Mini Review: 21st Century Environmental Challenges and the Interdisciplinary Research Efforts Rising to Meet Them

Globally, the livelihoods of over a billion people are affected by changes to marine ecosystems, both structurally and systematically. Resources and ecosystem services, provided by the marine environment, contribute nutrition, income, and health benefits for communities. One threat to these securities is ciguatera poisoning; worldwide, the most commonly reported non-bacterial seafood-related illness. Ciguatera is caused by the consumption of (primarily) finfish contaminated with ciguatoxins, potent neurotoxins produced by benthic single-cell microalgae. When consumed, ciguatoxins are biotransformed and can bioaccumulate throughout the food-web via complex pathways. Ciguatera-derived food insecurity is particularly extreme for small island-nations, where fear of intoxication can lead to fishing restrictions by region, species, or size. Exacerbating these complexities are anthropogenic or natural changes occurring in global marine habitats, e.g., climate change, greenhouse-gas induced physical oceanic changes, overfishing, invasive species, and even the international seafood trade. Here we provide an overview of the challenges and opportunities of the 21st century regarding the many facets of ciguatera, including the complex nature of this illness, the biological/environmental factors affecting the causative organisms, their toxins, vectors, detection methods, human-health oriented responses, and ultimately an outlook towards the future. Ciguatera research efforts face many social and environmental challenges this century. However, several future-oriented goals are within reach, including digital solutions for seafood supply chains, identifying novel compounds and methods with the potential for advanced diagnostics, treatments, and prediction capabilities. The advances described herein provide confidence that the tools are now available to answer many of the remaining questions surrounding ciguatera and therefore protection measures can become more accurate and routine.

Transcriptomic response patterns of hornyhead turbot (Pleuronichthys verticalis) dosed with polychlorinated biphenyls and polybrominated diphenyl ethers

To evaluate the impact of environmental contaminants on aquatic health, extensive surveys of fish populations have been conducted using bioaccumulation as an indicator of impairment. While these studies have reported mixtures of chemicals in fish tissues, the relationship between specific contaminants and observed adverse impacts remains poorly understood. The present study aimed to characterize the toxicological responses induced by persistent organic pollutants in wild-caught hornyhead turbot (P. verticalis). To do so, hornyhead turbot were interperitoneally injected with a single dose of PCB or PBDE congeners prepared using environmentally realistic mixture proportions. After 96-hour exposure, the livers were excised and analyzed using transcriptomic approaches and analytical chemistry. Concentrations of PCBs and PBDEs measured in the livers indicated clear differences across treatments, and congener profiles closely mirrored our expectations. Distinct gene profiles were characterized for PCB and PBDE exposed fish, with significant differences observed in the expression of genes associated with immune responses, endocrine-related functions, and lipid metabolism. Our findings highlight the key role that transcriptomics can play in monitoring programs to assess chemical-induced toxicity in heterogeneous group of fish (mixed gender and life stage) as is typically found during field surveys. Altogether, the present study provides further evidence of the potential of transcriptomic tools to improve aquatic health assessment and identify causative agents.

Genome-Wide Association Study of Growth and Body-Shape-Related Traits in Large Yellow Croaker (Larimichthys crocea) Using ddRAD Sequencing

Large yellow croaker (Larimichthys crocea) is an economically important marine fish species of China. Due to overfishing and marine pollution, the wild stocks of this croaker have collapsed in the past decades. Meanwhile, the cultured croaker is facing the difficulties of reduced genetic diversity and low growth rate. To explore the molecular markers related to the growth traits of croaker and providing the related SNPs for the marker-assisted selection, we used double-digest restriction-site associated DNA (ddRAD) sequencing to dissect the genetic bases of growth traits in a cultured population and identify the SNPs that associated with important growth traits by GWAS. A total of 220 individuals were genotyped by ddRAD sequencing. After quality control, 27,227 SNPs were identified in 220 samples and used for GWAS analysis. We identified 13 genome-wide significant associated SNPs of growth traits on 8 chromosomes, and the beta P of these SNPs ranged from 0.01 to 0.86. Through the definition of candidate regions and gene annotation, candidate genes related to growth were identified, including important regulators such as fgf18, fgf1, nr3c1, cyp8b1, fabp2, cyp2r1, ppara, and ccm2l. We also identified SNPs and candidate genes that significantly associated with body shape, including bmp7, col1a1, col11a2, and col18a1, which are also economically important traits for large yellow croaker aquaculture. The results provided insights into the genetic basis of growth and body shape in large yellow croaker population and would provide reliable genetic markers for molecular marker-assisted selection in the future. Meanwhile, the result established a basis for our subsequent fine mapping and related gene study.

Transcriptional Basis of Copper-Induced Olfactory Impairment in the Sea Lamprey, a Primitive Invasive Fish

Olfaction mediates behaviors necessary for survival and reproduction in fishes. Anthropogenic inputs of contaminants into aquatic environments, specifically copper, are known to disrupt a broad range of olfactory-mediated behaviors and can cause long-lasting damage even at low concentrations that have profound impacts on the biology of aquatic organisms. The sea lamprey (Petromyzon marinus) is a primitive fish species invasive to the North American Great Lakes that relies on olfaction to navigate during natal homing and in mate choice during reproduction. To investigate effects of copper on sea lamprey olfaction and the potential for maintenance of olfactory function during copper exposure, we exposed juvenile sea lamprey to environmentally ecologically relevant copper concentrations (0, 5, 10 and 30 µg/L) for 24 hr and characterized gene transcription response in olfactory tissue (i.e., peripheral olfactory organ and olfactory bulb) and forebrain using whole transcriptome sequencing. Copper exposure induced a pattern of positive dose-dependent transcriptional response. Expression changes primarily reflected up-regulation of genes involved in apoptosis and wound healing. Unlike higher vertebrates, genes specifically related to the olfactory senses of the sea lamprey, e.g., olfactory receptors, exhibited little transcriptional response to copper exposure, suggesting the mechanism of copper-induced olfactory impairment is through necrosis of the olfactory bulb and not copper-selective inhibition of olfactory receptors. Fully two-thirds of the differentially expressed genes at higher doses of copper have no known function and thus represent important candidates for further study of the responses to copper-induced olfactory injury. Our results shed light on the evolution of vertebrate olfactory repair mechanisms and have important implications for the conservation and management of both invasive and native populations of lamprey.

Non-coding RNAs Function as Immune Regulators in Teleost Fish

Non-coding RNAs (ncRNAs) are functional RNA molecules that are transcribed from DNA but not translated into proteins. ncRNAs function as key regulators of gene expression and chromatin modification. Recently, the functional role of ncRNAs in teleost fish has been extensively studied. Teleost fish are a highly diverse group among the vertebrate lineage. Fish are also important in terms of aquatic ecosystem, food production and human life, being the source of animal proteins worldwide and models of biomedical research. However, teleost fish always suffer from the invasion of infectious pathogens including viruses and bacteria, which has resulted in a tremendous economic loss to the fishing industry worldwide. Emerging evidence suggests that ncRNAs, especially miRNAs and lncRNAs, may serve as important regulators in cytokine and chemokine signaling, antigen presentation, complement and coagulation cascades, and T cell response in teleost fish. In this review, we summarize current knowledge and understanding of the roles of both miRNAs and lncRNAs in immune regulation in teleost fish. Molecular mechanism insights into the function of ncRNAs in fish immune response may contribute to the development of potential biomarkers and therapeutic targets for the prevention and treatment of fish diseases.

Monitoring sublethal changes in fish physiology following exposure to a light, unweathered crude oil

Biomarkers are frequently used to determine the exposure of fish to petroleum hydrocarbons following an oil spill. These biomarkers must be chosen carefully if they are to be used to determine sublethal toxic impacts as well as oil exposure. Many commonly used biomarkers relate to the metabolism of high molecular weight, typically pyrogenic, polycyclic aromatic hydrocarbons (PAHs), which are not abundant in unweathered crude oil. The goal of this study was to compare the efficacy of different biomarkers, including histological examination and transcriptomic profiling, in showing exposure to oil and the potential for sublethal toxic impacts. To achieve these goals, subadults/adults of the spotted dragonet (Repomucenus calcaratus) were exposed to a representative light, unweathered Australian oil for 96 h, so that the physiological changes that occur with exposure could be documented. Fish were then transferred to clean sediment for 90 h to quantify recovery. Biomarker changes, including PAH metabolites, 7-ethoxyresorufin O-deethylase (EROD), and histopathology, are presented in this work. In addition, a de novo transcriptome for the spotted dragonet was assembled, and differential transcript abundance was determined for the gill and liver of petroleum-exposed fish relative to a control. Increased levels of some biliary phenanthrene metabolites were seen throughout the exposure period. EROD levels showed modest, but not significant, increases. Transcriptomic differences were noted in the abundances of transcripts with a role in inflammation, primary metabolism and cardiac function. The patterns of transcript abundance in the gill and the liver changed in a manner that reflected exposure and recovery. The histology showed elevated prevalence of lesions, most notably vacuolization in liver and heart tissue, multi-organ necrosis, and lamellar epithelial lifting and telangiectasia in the gill. These findings suggest that short-term exposures to low molecular weight PAHs could elicit changes in the health of fish that are well predicted by the transcriptome. Furthermore, when light oil is released into the environment, exposure and subsequent risk would be better estimated using phenanthrene metabolite levels rather than EROD. This study also adds to the weight of evidence that exposure to low molecular weight PAHs may cause cardiac problems in fish. Further study is needed to determine the impact of these changes on reproductive capacity, long-term survival, and other population specific parameters.

Transcriptomic, lipid, and histological profiles suggest changes in health in fish from a pesticide hot spot

Barramundi (Lates calcarifer) were collected at the beginning (1st sampling) and end (2nd sampling) of the wet season from Sandy Creek, an agriculturally impacted catchment in the Mackay Whitsundays region of the Great Barrier Reef catchment area, and from Repulse Creek, located approximately 100 km north in Conway National Park, to assess the impacts of pesticide exposure. Gill and liver histology, lipid class composition in muscle, and the hepatic transcriptome were examined. The first sample of Repulse Creek fish showed little tissue damage and low transcript levels of xenobiotic metabolism enzymes. Sandy Creek fish showed altered transcriptomic patterns, including those that regulate lipid metabolism, xenobiotic metabolism, and immune response; gross histological alterations including lipidosis; and differences in some lipid classes. The second sampling of Repulse Creek fish showed similar alterations in hepatic transcriptome and tissue structure as fish from Sandy Creek. These changes may indicate a decrease in health of pesticide exposed fish.

Application of molecular tools to elucidate the microbiota of seafood

The aim of this review is to present the methodologies currently applied to identify microbiota and pathogens transmitted to humans through seafood consumption, focusing on molecular techniques and pointing out their importance, advantages, disadvantages and applicability. Knowledge of available techniques allows researchers to identify which technique best fits their expectations. With such discernment, it will be possible to infer which disadvantages will be present and, therefore, not interfering with the final result. Two methodologies can be employed for this purpose, dependent and independent cultures. However, the dependent culture has certain limitations that can be solved through the independent cultivation techniques, such as PCR, PFGE and NGS, especially through the sequencing of the 16S rRNA region, providing a complete view of microbial diversity. These have revolutionized microbiological knowledge, mainly because they allow for the identification of uncultivable micro-organisms, which represent a substantial portion of total micro-organisms, making it possible to elucidate not yet described taxa which may display pathogenic potential, besides quantifying microbial communities, microbiota genetics, translated proteins and produced metabolites. In addition, transcriptomic and metabolomic techniques also allow for the evaluation of possible impacts that microbial communities may create in their environment, as well as the determination of potential pathogenicity to humans.

The genome of the marine medaka Oryzias melastigma

Marine medaka (Oryzias melastigma) is considered to be a useful fish model for marine and estuarine ecotoxicology studies and has good potential for field-based population genomics because of its geographical distribution in Asian estuarine and coastal areas. In this study, we present the first whole-genome draft of O. melastigma. The genome assembly consists of 8,602 scaffolds (N50 = 23.737 Mb) and a total genome length of 779.4 Mb. A total of 23,528 genes were predicted, and 12,670 gene families shared with three teleost species (Japanese medaka, mangrove killifish and zebrafish) were identified. Genome analyses revealed that the O. melastigma genome is highly heterozygous and contains a large number of repeat sequences. This assembly represents a useful genomic resource for fish scientists.

Transcriptome Analysis Based on RNA-Seq in Understanding Pathogenic Mechanisms of Diseases and the Immune System of Fish: A Comprehensive Review

In recent years, with the advent of next-generation sequencing along with the development of various bioinformatics tools, RNA sequencing (RNA-Seq)-based transcriptome analysis has become much more affordable in the field of biological research. This technique has even opened up avenues to explore the transcriptome of non-model organisms for which a reference genome is not available. This has made fish health researchers march towards this technology to understand pathogenic processes and immune reactions in fish during the event of infection. Recent studies using this technology have altered and updated the previous understanding of many diseases in fish. RNA-Seq has been employed in the understanding of fish pathogens like bacteria, virus, parasites, and oomycetes. Also, it has been helpful in unraveling the immune mechanisms in fish. Additionally, RNA-Seq technology has made its way for future works, such as genetic linkage mapping, quantitative trait analysis, disease-resistant strain or broodstock selection, and the development of effective vaccines and therapies. Until now, there are no reviews that comprehensively summarize the studies which made use of RNA-Seq to explore the mechanisms of infection of pathogens and the defense strategies of fish hosts. This review aims to summarize the contemporary understanding and findings with regard to infectious pathogens and the immune system of fish that have been achieved through RNA-Seq technology.

Marine copepod cytochrome P450 genes and their applications for molecular ecotoxicological studies in response to oil pollution

Recently, accidental spills of heavy oil have caused adverse effects in marine organisms. Oil pollution can induce damages on development and reproduction, linking with detrimental effects on diverse molecular levels of genes and proteins in plankton and fish. However, most information was mainly focused on marine vertebrates and consequently, limited information was available in marine invertebrates. Furthermore, there is still a lack of knowledge bridging in vivo endpoints with the functional regulation of cytochrome P450 (CYP) genes in response to oil spill pollution in marine invertebrates. In this paper, adverse effects of oil spill pollution in marine invertebrates are summarized with the importance of CYP genes as a potential biomarker, applying for environmental monitoring to detect oil spill using marine copepods.

Hepatic transcriptomic profiles from barramundi, Lates calcarifer, as a means of assessing organism health and identifying stressors in rivers in northern Queensland

Resource managers need to differentiate between sites with and without contaminants and those where contaminants cause impacts. Potentially, transcriptomes could be used to evaluate sites where contaminant-induced effects may occur, to identify causative stressors of effects and potential adverse outcomes. To test this hypothesis, the hepatic transcriptomes in Barramundi, a perciforme teleost fish, (Lates calcarifer) from two reference sites, two agriculturally impacted sites sampled during the dry season, and an impacted site sampled during the wet season were compared. The hepatic transcriptome was profiled using RNA-Seq. Multivariate analysis showed that transcriptomes were clustered based on site and by inference water quality, but not sampling time. The largest differences in transcriptomic profile were between reference sites and a site sampled during high run-off, showing that impacted sites can be identified via RNA-Seq. Transcripts with altered abundance were linked to xenobiotic metabolism, peroxisome proliferation and stress responses, indicating putative stressors with the potential for adverse outcomes in barramundi.

Characterisation of the transcriptome of male and female wild-type guppy brains with RNA-Seq and consequences of exposure to the pharmaceutical pollutant, 17α-ethinyl estradiol

Waterways are increasingly being contaminated by chemical compounds that can disrupt the endocrinology of organisms. One such compound is 17α-ethinyl estradiol (EE2), a synthetic estrogen used in the contraceptive pill. Despite considerable research interest in the effects of EE2 on reproduction and gene expression, surprisingly, only a few studies have capitalised on technologies, such as next-generation sequencing (NGS), to uncover the molecular pathways related to EE2 exposure. Accordingly, using high-throughput sequencing technologies, the aim of our study was to explore the effects of EE2 on brain transcriptome in wild-type male and female guppy (Poecilia reticulata). We conducted two sets of experiments, where fish were exposed to EE2 (measured concentrations: 8ng/L and 38ng/L) in a flow-through system for 21days. The effects on the brain transcriptome on both males and females were assessed using Illumina sequencing (MiSeq and HiSeq) platform followed by bioinformatics analysis (edgeR, DESeq2). Here, we report that exposure to EE2 caused both up- and downregulation of specific transcript abundances, and affected transcript abundance in a sex-specific manner. Specifically, we found 773 transcripts, of which 60 were male-specific, 61 female-specific and 285 treatment-specific. EE2 affected expression of 165 transcripts in males, with 88 downregulated and 77 upregulated, while in females, 120 transcripts were affected with 62 downregulated and 58 upregulated. Finally, RT-qPCR validation demonstrated that expression of transcripts related to transposable elements, neuroserpin and heat shock protein were significantly affected by EE2-exposure. Our study is the first to report brain transcriptome libraries for guppies exposed to EE2. Not only does our study provide a valuable resource, it offers insights into the mechanisms underlying the feminizing effects on the brains of organisms exposed to environmentally realistic concentrations of EE2.

Natural environmental impacts on teleost immune function

The environment in which teleosts exist can experience considerable change. Short-term changes can occur in relation to tidal movements or adverse weather events. Long-term changes can be caused by anthropogenic impacts such as climate change, which can result in changes to temperature, acidity, salinity and oxygen capacity of aquatic environments. These changes can have important impacts on the physiology of an animal, including its immune system. This can have consequences on the well-being of the animal and its ability to protect against pathogens. This review will look at recent investigations of these types of environmental change on the immune response in teleosts.

Omics of the marine medaka (Oryzias melastigma) and its relevance to marine environmental research

In recent years, the marine medaka (Oryzias melastigma), also known as the Indian medaka or brackish medaka, has been recognized as a model fish species for ecotoxicology and environmental research in the Asian region. O. melastigma has several promising features for research, which include a short generation period (3-4 months), daily spawning, small size (3-4 cm), transparent embryos, sexual dimorphism, and ease of mass culture in the laboratory. There have been extensive transcriptome and genome studies on the marine medaka in the past decade. Such omics data can be useful in understanding the signal transduction pathways of small teleosts in response to environmental stressors. An omics-integrated approach in the study of the marine medaka is important for strengthening its role as a small fish model for marine environmental studies. In this review, we present current omics information about the marine medaka and discuss its potential applications in the study of various molecular pathways that can be targets of marine environmental stressors, such as chemical pollutants. We believe that this review will encourage the use of this small fish as a model species in marine environmental research.

Defining the role of omics in assessing ecosystem health: Perspectives from the Canadian environmental monitoring program

Scientific reviews and studies continue to describe omics technologies as the next generation of tools for environmental monitoring, while cautioning that there are limitations and obstacles to overcome. However, omics has not yet transitioned into national environmental monitoring programs designed to assess ecosystem health. Using the example of the Canadian Environmental Effects Monitoring (EEM) program, the authors describe the steps that would be required for omics technologies to be included in such an established program. These steps include baseline collection of omics endpoints across different species and sites to generate a range of what is biologically normal within a particular ecosystem. Natural individual variability in the omes is not adequately characterized and is often not measured in the field, but is a key component to an environmental monitoring program, to determine the critical effect size or action threshold for management. Omics endpoints must develop a level of standardization, consistency, and rigor that will allow interpretation of the relevance of changes across broader scales. To date, population-level consequences of routinely measured endpoints such as reduced gonad size or intersex in fish is not entirely clear, and the significance of genome-wide molecular, proteome, or metabolic changes on organism or population health is further removed from the levels of ecological change traditionally managed. The present review is not intended to dismiss the idea that omics will play a future role in large-scale environmental monitoring studies, but rather outlines the necessary actions for its inclusion in regulatory monitoring programs focused on assessing ecosystem health.

Draft genome sequence of the silver pomfret fish, Pampus argenteus

Silver pomfret, Pampus argenteus, is a fish species from coastal waters. Despite its high commercial value, this edible fish has not been sequenced. Hence, its genetic and genomic studies have been limited. We report the first draft genome sequence of the silver pomfret obtained using a Next Generation Sequencing (NGS) technology. We assembled 38.7 Gb of nucleotides into scaffolds of 350 Mb with N50 of about 1.5 kb, using high quality paired end reads. These scaffolds represent 63.7% of the estimated silver pomfret genome length. The newly sequenced and assembled genome has 11.06% repetitive DNA regions, and this percentage is comparable to that of the tilapia genome. The genome analysis predicted 16 322 genes. About 91% of these genes showed homology with known proteins. Many gene clusters were annotated to protein and fatty-acid metabolism pathways that may be important in the context of the meat texture and immune system developmental processes. The reference genome can pave the way for the identification of many other genomic features that could improve breeding and population-management strategies, and it can also help characterize the genetic diversity of P. argenteus.

Comprehensive identification and profiling of host miRNAs in response to Singapore grouper iridovirus (SGIV) infection in grouper (Epinephelus coioides)

microRNAs (miRNAs) are an evolutionarily conserved class of non-coding RNA molecules that participate in various biological processes. Employment of high-throughput screening strategies greatly prompts the investigation and profiling of miRNAs in diverse species. In recent years, grouper (Epinephelus spp.) aquaculture was severely affected by iridoviral diseases. However, knowledge regarding the host immune responses to viral infection, especially the miRNA-mediated immune regulatory roles, is rather limited. In this study, by employing Solexa deep sequencing approach, we identified 116 grouper miRNAs from grouper spleen-derived cells (GS). As expected, these miRNAs shared high sequence similarity with miRNAs identified in zebrafish (Danio rerio), pufferfish (Fugu rubripes), and other higher vertebrates. In the process of Singapore grouper iridovirus (SGIV) infection, 45 and 43 miRNAs with altered expression (>1.5-fold) were identified by miRNA microarray assays in grouper spleen tissues and GS cells, respectively. Furthermore, target prediction revealed 189 putative targets of these grouper miRNAs.

Omics for aquatic ecotoxicology: control of extraneous variability to enhance the analysis of environmental effects

There are multiple sources of biological and technical variation in a typical ecotoxicology study that may not be revealed by traditional endpoints but that become apparent in an omics dataset. As researchers increasingly apply omics technologies to environmental studies, it will be necessary to understand and control the main source(s) of variability to facilitate meaningful interpretation of such data. For instance, can variability in omics studies be addressed by changing the approach to study design and data analysis? Are there statistical methods that can be employed to correctly interpret omics data and make use of unattributed, inherent variability? The present study presents a review of experimental design and statistical considerations applicable to the use of omics methods in systems toxicology studies. In addition to highlighting potential sources that contribute to experimental variability, this review suggests strategies with which to reduce and/or control such variability so as to improve reliability, reproducibility, and ultimately the application of omics data for systems toxicology.

454 pyrosequencing-based analysis of gene expression profiles in the amphipod Melita plumulosa: transcriptome assembly and toxicant induced changes

Next generation sequencing using Roche's 454 pyrosequencing platform can be used to generate genomic information for non-model organisms, although there are bioinformatic challenges associated with these studies. These challenges are compounded by a lack of a standardized protocol to either assemble data or to evaluate the quality of a de novo transcriptome. This study presents an assembly of the control and toxicant responsive transcriptome of Melita plumulosa, an Australian amphipod commonly used in ecotoxicological studies. RNA was harvested from control amphipods, juvenile amphipods, and from amphipods exposed to either metal or diesel contaminated sediments. This RNA was used as the basis for a 454 based transcriptome sequencing effort. Sequencing generated 1.3 million reads from control, juvenile, metal-exposed and diesel-exposed amphipods. Different read filtering and assembly protocols were evaluated to generate an assembly that (i) had an optimal number of contigs; (ii) had long contigs; (iii) contained a suitable representation of conserved genes; and (iv) had long ortholog alignment lengths relative to the length of each contig. A final assembly, generated using fixed-length trimming based on the sequence quality scores, followed by assembly using the MIRA algorithm, produced the best results. The 26,625 contigs generated via this approach were annotated using Blast2GO, and the differential expression between treatments and control was determined by mapping with BWA followed by DESeq. Although the mapping generated low coverage, many differentially expressed contigs, including some with known developmental or toxicological function, were identified. This study demonstrated that 454 pyrosequencing is an effective means of generating reference transcriptome information for organisms, such as the amphipod M. plumulosa, that have no genomic information available in databases or in closely related sequenced species. It also demonstrated how optimization of read filtering protocols and assembly approaches changes the utility of results obtained from next generation sequencing studies, and establishes criteria to determine the quality of a de novo assembly in species lacking a reference genome. This new transcriptomic knowledge provides the genomic foundation for the creation of microarray and qPCR assays, serving as a reference transcriptome in future RNAseq studies, and allowing both the biology and ecotoxicology of this organism to be better understood. This approach will allow genomics-based methodology to be applied to a wider range of environmentally relevant species.

Assessing mechanisms of toxicant response in the amphipod Melita plumulosa through transcriptomic profiling

This study describes the function of transcripts with altered abundance in the epibenthic amphipod, Melita plumulosa, following whole-sediment exposure to a series of common environmental contaminants. M. plumulosa were exposed for 48 h to sediments spiked and equilibrated with the following contaminants at concentrations predicted to cause sublethal effects to reproduction: porewater ammonia 30 mg L(-1); bifenthrin at 100 μg kg(-1); fipronil at 50 μg kg(-1); 0.6% diesel; 0.3% crude oil; 250 mg Cu kg(-1); 400 mg Ni kg(-1); and 400 mg Zn kg(-1). RNA was extracted and hybridized against a custom Agilent microarray developed for this species. Although the microarray represented a partial transcriptome and not all features on the array could be annotated, unique transcriptomic profiles were generated for each of the contaminant exposures. Hierarchical clustering grouped the expression profiles together by contaminant class, with copper and zinc, the petroleum products and nickel, and the pesticides each forming a distinct cluster. Many of the transcriptional changes observed were consistent with patterns previously described in other crustaceans. The changes in the transcriptome demonstrated that contaminant exposure caused changes in digestive function, growth and moulting, and the cytoskeleton following metal exposure, whereas exposure to petroleum products caused changes in carbohydrate metabolism, xenobiotic metabolism and hormone cycling. Functional analysis of these gene expression profiles can provide a better understanding of modes of toxic action and permits the prediction of mixture effects within contaminated ecosystems.

Assessing gene network stability and individual variability in the fathead minnow (Pimephales promelas) transcriptome

Transcriptomics is increasingly used to assess biological responses to environmental stimuli and stressors such as aquatic pollutants. However, fundamental studies characterizing individual variability in mRNA levels are lacking, which currently limits the use of transcriptomics in environmental monitoring assessments. To address individual variability in transcript abundance, we performed a meta-analysis on 231 microarrays that were conducted in the fathead minnow (FHM), a widely used toxicological model. The mean variability for gene probes was ranked from most to least variable based upon the coefficient of variation. Transcripts that were the most variable in individual tissues included NADH dehydrogenase flavoprotein 1, GTPase IMAP family member 7-like and v-set domain-containing T-cell activation inhibitor 1-like while genes encoding ribosomal proteins (rpl24 and rpl36), basic transcription factor 3, and nascent polypeptide-associated complex alpha subunit were the least variable in individuals across a range of microarray experiments. Gene networks that showed high variability (based upon the variation in expression of individual members within the network) included cell proliferation, metabolism (steroid, lipids, and glucose), cell adhesion, vascularization, and regeneration while those that showed low variability (more stability) included mRNA and rRNA processing, regulation of translational fidelity, RNA splicing, and ribosome biogenesis. Real-time PCR was conducted on a subset of genes for comparison of variability collected from the microarrays. There was a significant positive relationship between the two methods when measuring individual variability, suggesting that variability detected in microarray data can be used to guide decisions on sample sizes for measuring transcripts in real-time PCR experiments. A power analysis revealed that measuring estrogen receptor ba (esrba) requires fewer biological replicates than that of estrogen receptor bb (esrbb) in the gonad and samples sizes required to detect a 50% change for reproductive-related transcripts is between 12 and 20. Characterizing individual variability at the molecular level will prove necessary as efforts are made toward integrating molecular tools into environmental risk assessments.

Transcriptome-wide identification and characterization of the Procambarus clarkii microRNAs potentially related to immunity against Spiroplasma eriocheiris infection

We used the Illumina/Solexa deep sequencing technology to sequence two small RNA libraries prepared from hemocytes of Procambarus clarkii under normal and infection conditions. The high-throughput sequencing approach resulted in approximately 12,801,827 and 8,410,455 raw reads corresponding to 10,949,754 and 6,648,161 high-quality mappable reads for the normal and infected hemocyte samples, respectively. Bioinformatic analyses identified 195 unique miRNAs, including 30 that are conserved in crustaceans, 48 that are novel to crayfish but are present in other arthropods (PN-type), and 117 that are completely new (PC-type). Thirty-three miRNAs displayed significant differential expressions between the two hemocyte samples (p < 0.0001). Of these, 15 (45.5%) were significantly up-regulated and 18 (54.5%) were significantly down-regulated upon challenge with Spiroplasma eriocheiris. Integrating comparative genomic and bibliomic analysis, of the 33 significant miRNAs identified, 19 were conserved and immune-related in P. clarkii and Eriocheir sinensis infected with S. eriocheiris infection; 24 were conserved and immune-related in P. clarkii and Marsupenaeus japonicus immune response to S. eriocheiris or white spot syndrome virus (WSSV) infection. Function annotation of target genes revealed a broad range of biological processes and signal transduction pathways that regulated by crayfish miRNAs. Thereinto, pcl-miR-34, pcl-miR-7, PN-pcl-let-7, pcl-miR-1, and pcl-miR-2b are highly conserved in vertebrates and invertebrates and function in the similar pathways. To our knowledge, this is the first report of comprehensive identification of P. clarkii miRNAs and of expression analysis of P. clarkii miRNAs after exposure to S. eriocheiris in crayfish, and many miRNAs were differentially regulated under normal and infection conditions. Our results should help develop new control strategies for efficient immune protection against S. eriocheiris infections in crustaceans.

Induction of gene responses in St. Lawrence River northern pike (Esox lucius) environmentally exposed to perfluorinated compounds

Municipal waste water effluents (MWWEs) are important sources of chemical contamination for aquatic environments. This study investigated the presence and effects of perfluorinated compounds (PFCs) in environmentally exposed northern pike (Esox lucius) collected upstream and downstream a major municipal waste water treatment plant (WWTP) in the St. Lawrence River, Canada. Twelve PFCs, including the newly detected perfluoroethylcyclohexane sulfonate (PFECHS), were quantified in fish muscle, liver, and plasma. Additionally, the expression of eight genes and the activity of three biomarkers were analyzed in fish tissues at both sites. Mean ∑PFC concentration in fish plasma collected upstream the WWTP was 185ng/g w.w. compared to 545ng/g w.w. downstream the point of release. PFECHS was quantified for the first time in St. Lawrence River fish (mean plasma concentration in MWWE fish: 5.07±4.72ng/g w.w.). Results of transcriptomic responses were tissue-specific and indicated significant up-regulation for metallothionein (MT) in blood and MT, glutathion-S-transferase (GST), superoxide dismutase (SOD), and cytochromes P450 1A1 (CYP1A1) in gill tissue of fish collected in the MWWE suggesting greater stress responses for organisms at this location. Significant relationships were found between vitellogenin (Vtg) gene expression in liver, Vtg activity in plasma and perfluorotridecanoic acid (PFTrA), perfluorotetradecanoic acid (PFTeA), and perfluorodecane sulfonate (PFDS) plasma concentrations. The possible endocrine effects of these PFCs should be further investigated.

Genome wide analysis of Silurana (Xenopus) tropicalis development reveals dynamic expression using network enrichment analysis

Development involves precise timing of gene expression and coordinated pathways for organogenesis and morphogenesis. Functional and sub-network enrichment analysis provides an integrated approach for identifying networks underlying development. The objectives of this study were to characterize early gene regulatory networks over Silurana tropicalis development from NF stage 2 to 46 using a custom Agilent 4×44K microarray. There were >8000 unique gene probes that were differentially expressed between Nieuwkoop-Faber (NF) stage 2 and stage 16, and >2000 gene probes differentially expressed between NF 34 and 46. Gene ontology revealed that genes involved in nucleosome assembly, cell division, pattern specification, neurotransmission, and general metabolism were increasingly regulated throughout development, consistent with active development. Sub-network enrichment analysis revealed that processes such as membrane hyperpolarisation, retinoic acid, cholesterol, and dopamine metabolic gene networks were activated/inhibited over time. This study identifies RNA transcripts that are potentially maternally inherited in an anuran species, provides evidence that the expression of genes involved in retinoic acid receptor signaling may increase prior to those involved in thyroid receptor signaling, and characterizes novel gene expression networks preceding organogenesis which increases understanding of the spatiotemporal embryonic development in frogs.

Promises and problems for the new paradigm for risk assessment and an alternative approach involving predictive systems models

The need for cost-effective risk assessment of chemicals is leading to the development of a reductionist paradigm that tries to assess impacts on humans and ecosystems from molecular changes. However, the biggest challenge for this paradigm comes from the emergence of properties that arise out of the interactions of the parts that are not included and yet which are key for assessing likely impacts. Although identifying key events and adverse outcome pathways can shed light on the involvement of important metabolic processes in toxicity, this does not mean that particular molecular initiating events are likely to be robust or accurate predictors of impacts that matter. There are even greater challenges for the new paradigm applied to ecological systems than to human health because of the need to link across more levels of biological organization. The present study argues for a predictive systems approach that makes the linkages through systems models in a mechanistic way that allows for emergence and that also has the potential for reducing the costs and use of animals in ecological risk assessments.