De novo assembly and annotation of the Patagonian toothfish (Dissostichus eleginoides) genome

BACKGROUND

Patagonian toothfish (Dissostichus eleginoides) is an economically and ecologically important fish species in the family Nototheniidae. Juveniles occupy progressively deeper waters as they mature and grow, and adults have been caught as deep as 2500 m, living on or in just above the southern shelves and slopes around the sub-Antarctic islands of the Southern Ocean. As apex predators, they are a key part of the food web, feeding on a variety of prey, including krill, squid, and other fish. Despite its importance, genomic sequence data, which could be used for more accurate dating of the divergence between Patagonian and Antarctic toothfish, or establish whether it shares adaptations to temperature with fish living in more polar or equatorial climes, has so far been limited.

RESULTS

A high-quality D. eleginoides genome was generated using a combination of Illumina, PacBio and Omni-C sequencing technologies. To aid the genome annotation, the transcriptome derived from a variety of toothfish tissues was also generated using both short and long read sequencing methods. The final genome assembly was 797.8 Mb with a N50 scaffold length of 3.5 Mb. Approximately 31.7% of the genome consisted of repetitive elements. A total of 35,543 putative protein-coding regions were identified, of which 50% have been functionally annotated. Transcriptomics analysis showed that approximately 64% of the predicted genes (22,617 genes) were found to be expressed in the tissues sampled. Comparative genomics analysis revealed that the anti-freeze glycoprotein (AFGP) locus of D. eleginoides does not contain any AFGP proteins compared to the same locus in the Antarctic toothfish (Dissostichus mawsoni). This is in agreement with previously published results looking at hybridization signals and confirms that Patagonian toothfish do not possess AFGP coding sequences in their genome.

CONCLUSIONS

We have assembled and annotated the Patagonian toothfish genome, which will provide a valuable genetic resource for ecological and evolutionary studies on this and other closely related species.

Transcriptomic Responses to Koi Herpesvirus in Isolated Blood Leukocytes from Infected Common Carp

Koi herpesvirus (KHV, CyHV-3) causes severe economic losses in carp farms. Its eradication is challenging due to the establishment of latency in blood leukocytes and other tissues. To understand the molecular mechanisms leading to KHV infection in leukocytes, common carp were bath-exposed to KHV at 17 °C. After confirming the presence of viral transcripts in blood leukocytes at ten days post infection, RNA-Seq was performed on peripheral blood leukocytes on the Illumina NovaSeq. KHV infection triggered a robust immune response mediated by pattern recognition receptors, mainly toll-like receptors (tlr2, tlr5, tlr7, and tlr13), urokinase plasminogen activator surface receptor-like, galectin proteins, and lipid mediators such as leukotriene B4 receptor 1. Enriched pathways showed increased mitochondria oxidative phosphorylation and the activation of signalling pathways such as mitogen-activated protein kinases (MAPKs) and vascular endothelial growth factor (VEGF). KHV-infected leukocytes showed low production of reactive oxygen species (ROS) and glutathione metabolism, high iron export and phagocytosis activity, and low autophagy. Macrophage polarization was deduced from the up-regulation of genes such as arginase non-hepatic 1-like, macrophage mannose receptor-1, crem, il-10, and il-13 receptors, while markers for cytotoxic T cells were observed to be down-regulated. Further work is required to characterise these leukocyte subsets and the molecular events leading to KHV latency in blood leukocytes.

Utility of wastewater genomic surveillance compared to clinical surveillance to track the spread of the SARS-CoV-2 Omicron variant across England

The world has moved into a new stage of managing the SARS-CoV-2 pandemic with minimal restrictions and reduced testing in the population, leading to reduced genomic surveillance of virus variants in individuals. Wastewater-based epidemiology (WBE) can provide an alternative means of tracking virus variants in the population but decision-makers require confidence that it can be applied to a national scale and is comparable to individual testing data. We analysed 19,911 samples from 524 wastewater sites across England at least twice a week between November 2021 and February 2022, capturing sewage from >70% of the English population. We used amplicon-based sequencing and the phylogeny based de-mixing tool Freyja to estimate SARS-CoV-2 variant frequencies and compared these to the variant dynamics observed in individual testing data from clinical and community settings. We show that wastewater data can reconstruct the spread of the Omicron variant across England since November 2021 in close detail and aligns closely with epidemiological estimates from individual testing data. We also show the temporal and spatial spread of Omicron within London. Our wastewater data further reliably track the transition between Omicron subvariants BA1 and BA2 in February 2022 at regional and national levels. Our demonstration that WBE can track the fast-paced dynamics of SARS-CoV-2 variant frequencies at a national scale and closely match individual testing data in time shows that WBE can reliably fill the monitoring gap left by reduced individual testing in a more affordable way.

Comparison of metagenomic and targeted methods for sequencing human pathogenic viruses from wastewater

Wastewater-based epidemiology is a powerful tool for monitoring the emergence and spread of viral pathogens at the population scale. Typical polymerase chain reaction (PCR)-based methods of quantitative and genomic monitoring of viruses in wastewater provide high sensitivity and specificity. However, these methods are limited to the surveillance of target viruses in a single assay and require prior knowledge of the target genome(s). Metagenomic sequencing methods may represent a target-agnostic approach to viral wastewater monitoring, allowing for the detection of a broad range of target viruses, including potentially novel and emerging pathogens. In this study, targeted and untargeted metagenomic sequencing methods were compared with tiled-PCR sequencing for the detection and genotyping of viral pathogens in wastewater samples. Deep shotgun metagenomic sequencing was unable to generate sufficient genome coverage of human pathogenic viruses for robust genomic epidemiology, with samples dominated by bacteria. Hybrid-capture enrichment of shotgun libraries for respiratory viruses led to significant increases in genome coverage for a range of targets. Tiled-PCR sequencing led to further improvements in genome coverage compared to hybrid capture for severe acute respiratory syndrome coronavirus 2, enterovirus D68, norovirus GII, and human adenovirus F41 in wastewater samples. In conclusion, untargeted shotgun sequencing was unsuitable for genomic monitoring of the low virus concentrations in wastewater samples analyzed in this study. Hybrid-capture enrichment represented a viable method for simultaneous genomic epidemiology of a range of viral pathogens, while tiled-PCR sequencing provided the optimal genome coverage for individual viruses with the minimum sequencing depth. IMPORTANCE Most public health initiatives that monitor viruses in wastewater have utilized quantitative polymerase chain reaction (PCR) and whole genome PCR sequencing, mirroring techniques used for viral epidemiology in individuals. These techniques require prior knowledge of the target viral genome and are limited to monitoring individual or small groups of viruses. Metagenomic sequencing may offer an alternative strategy for monitoring a broad spectrum of viruses in wastewater, including novel and emerging pathogens. In this study, while amplicon sequencing gave high viral genome coverage, untargeted shotgun sequencing of total nucleic acid samples was unable to detect human pathogenic viruses with enough sensitivity for use in genomic epidemiology. Enrichment of shotgun libraries for respiratory viruses using hybrid-capture technology provided genotypic information on a range of viruses simultaneously, indicating strong potential for wastewater surveillance. This type of targeted metagenomics could be used for monitoring diverse targets, such as pathogens or antimicrobial resistance genes, in environmental samples.

Optimised protocol for monitoring SARS-CoV-2 in wastewater using reverse complement PCR-based whole-genome sequencing

Monitoring the spread of viral pathogens in the population during epidemics is crucial for mounting an effective public health response. Understanding the viral lineages that constitute the infections in a population can uncover the origins and transmission patterns of outbreaks and detect the emergence of novel variants that may impact the course of an epidemic. Population-level surveillance of viruses through genomic sequencing of wastewater captures unbiased lineage data, including cryptic asymptomatic and undiagnosed infections, and has been shown to detect infection outbreaks and novel variant emergence before detection in clinical samples. Here, we present an optimised protocol for quantification and sequencing of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) in influent wastewater, used for high-throughput genomic surveillance in England during the COVID-19 pandemic. This protocol utilises reverse compliment PCR for library preparation, enabling tiled amplification across the whole viral genome and sequencing adapter addition in a single step to enhance efficiency. Sequencing of synthetic SARS-CoV-2 RNA provided evidence validating the efficacy of this protocol, while data from high-throughput sequencing of wastewater samples demonstrated the sensitivity of this method. We also provided guidance on the quality control steps required during library preparation and data analysis. Overall, this represents an effective method for high-throughput sequencing of SARS-CoV-2 in wastewater which can be applied to other viruses and pathogens of humans and animals.

An integrated eco-evolutionary framework to predict population-level responses of climate-sensitive pathogens

It is critical to gain insight into how climate change impacts evolutionary responses within climate-sensitive pathogen populations, such as increased resilience, opportunistic responses and the emergence of dominant variants from highly variable genomic backgrounds and subsequent global dispersal. This review proposes a framework to support such analysis, by combining genomic evolutionary analysis with climate time-series data in a novel spatiotemporal dataframe for use within machine learning applications, to understand past and future evolutionary pathogen responses to climate change. Recommendations are presented to increase the feasibility of interdisciplinary applications, including the importance of robust spatiotemporal metadata accompanying genome submission to databases. Such workflows will inform accessible public health tools and early-warning systems, to aid decision-making and mitigate future human health threats.

Evaluation of variant calling algorithms for wastewater-based epidemiology using mixed populations of SARS-CoV-2 variants in synthetic and wastewater samples

Wastewater-based epidemiology has been used extensively throughout the COVID-19 (coronavirus disease 19) pandemic to detect and monitor the spread and prevalence of SARS-CoV-2 (severe acute respiratory syndrome coronavirus 2) and its variants. It has proven an excellent, complementary tool to clinical sequencing, supporting the insights gained and helping to make informed public-health decisions. Consequently, many groups globally have developed bioinformatics pipelines to analyse sequencing data from wastewater. Accurate calling of mutations is critical in this process and in the assignment of circulating variants; yet, to date, the performance of variant-calling algorithms in wastewater samples has not been investigated. To address this, we compared the performance of six variant callers (VarScan, iVar, GATK, FreeBayes, LoFreq and BCFtools), used widely in bioinformatics pipelines, on 19 synthetic samples with known ratios of three different SARS-CoV-2 variants of concern (VOCs) (Alpha, Beta and Delta), as well as 13 wastewater samples collected in London between the 15th and 18th December 2021. We used the fundamental parameters of recall (sensitivity) and precision (specificity) to confirm the presence of mutational profiles defining specific variants across the six variant callers. Our results show that BCFtools, FreeBayes and VarScan found the expected variants with higher precision and recall than GATK or iVar, although the latter identified more expected defining mutations than other callers. LoFreq gave the least reliable results due to the high number of false-positive mutations detected, resulting in lower precision. Similar results were obtained for both the synthetic and wastewater samples.

Resistance to white spot syndrome virus in the European shore crab is associated with suppressed virion trafficking and heightened immune responses

Introduction

All decapod crustaceans are considered potentially susceptible to White Spot Syndrome Virus (WSSV) infection, but the degree of White Spot Disease (WSD) susceptibility varies widely between species. The European shore crab Carcinus maenas can be infected with the virus for long periods of time without signs of disease. Given the high mortality rate of susceptible species, the differential susceptibility of these resistant hosts offers an opportunity to investigate mechanisms of disease resistance.

Methods

Here, the temporal transcriptional responses (mRNA and miRNA) of C. maenas following WSSV injection were analysed and compared to a previously published dataset for the highly WSSV susceptible Penaeus vannamei to identify key genes, processes and pathways contributing to increased WSD resistance.

Results

We show that, in contrast to P. vannamei, the transcriptional response during the first 2 days following WSSV injection in C. maenas is limited. During the later time points (7 days onwards), two groups of crabs were identified, a recalcitrant group where no replication of the virus occurred, and a group where significant viral replication occurred, with the transcriptional profiles of the latter group resembling those of WSSV-susceptible species. We identify key differences in the molecular responses of these groups to WSSV injection.

Discussion

We propose that increased WSD resistance in C. maenas may result from impaired WSSV endocytosis due to the inhibition of internal vesicle budding by dynamin-1, and a delay in movement to the nucleus caused by the downregulation of cytoskeletal transcripts required for WSSV cytoskeleton docking, during early stages of the infection. This response allows resistant hosts greater time to fine-tune immune responses associated with miRNA expression, apoptosis and the melanisation cascade to defend against, and clear, invading WSSV. These findings suggest that the initial stages of infection are key to resistance to WSSV in the crab and highlight possible pathways that could be targeted in farmed crustacean to enhance resistance to WSD.

Txikispora philomaios n. sp., n. g., a Micro-Eukaryotic Pathogen of Amphipods, Reveals Parasitism and Hidden Diversity in Class Filasterea

This study provides a morphological, ultrastructural, and phylogenetic characterization of a novel micro-eukaryotic parasite (2.3-2.6 µm) infecting amphipod genera Echinogammarus and Orchestia. Longitudinal studies across two years revealed that infection prevalence peaked in late April and May, reaching 64% in Echinogammarus sp. and 15% in Orchestia sp., but was seldom detected during the rest of the year. The parasite infected predominantly haemolymph, connective tissue, tegument, and gonad, although hepatopancreas and nervous tissue were affected in heavier infections, eliciting melanization and granuloma formation. Cell division occurred inside walled parasitic cysts, often within host haemocytes, resulting in haemolymph congestion. Small subunit (18S) rRNA gene phylogenies including related environmental sequences placed the novel parasite as a highly divergent lineage within Class Filasterea, which together with Choanoflagellatea represent the closest protistan relatives of Metazoa. We describe the new parasite as Txikispora philomaios n. sp. n. g., the first confirmed parasitic filasterean lineage, which otherwise comprises four free-living flagellates and a rarely observed endosymbiont of snails. Lineage-specific PCR probing of other hosts and surrounding environments only detected T. philomaios in the platyhelminth Procerodes sp. We expand the known diversity of Filasterea by targeted searches of metagenomic datasets, resulting in 13 previously unknown lineages from environmental samples.

Social partner cooperativeness influences brain oxytocin transcription in Trinidadian guppies (Poecilia reticulata)

For non-kin cooperation to be maintained, individuals need to respond adaptively to the cooperative behaviour of their social partners. Currently, however, little is known about the biological responses of individuals to experiencing cooperation. Here, we quantify the neuroregulatory response of Trinidadian guppies (Poecilia reticulata) experiencing cooperation or defection by examining the transcriptional response of the oxytocin gene (oxt; also known as isotocin), which has been implicated in cooperative decision-making. We exposed wild-caught females to social environments where partners either cooperated or defected during predator inspection, or to a control (non-predator inspection) context, and quantified the relative transcription of the oxt gene. We tested an experimental group, originating from a site where individuals are under high predation threat and have previous experience of large aquatic predators (HP), and a control group, where individuals are under low predation threat and naïve to large aquatic predators (LP). LP, but not HP, fish showed different behavioural responses to the behaviour of their social environment, cooperating with cooperative partners and defecting when paired with defecting ones. In HP, but not LP, fish brain mid-section oxt relative transcription varied depending on social partner behaviour. HP fish experiencing cooperation during predator inspection had lower oxt transcription than those experiencing defection. This effect was not present in the control population or in the control context, where the behaviour of social partners did not affect oxt transcription. Our findings provide insight into the neuromodulation underpinning behavioural responses to social experiences, and ultimately to the proximate mechanisms underlying social decision-making.

Global mRNA and miRNA Analysis Reveal Key Processes in the Initial Response to Infection with WSSV in the Pacific Whiteleg Shrimp

White Spot Disease (WSD) presents a major barrier to penaeid shrimp production. Mechanisms underlying White Spot Syndrome Virus (WSSV) susceptibility in penaeids are poorly understood due to limited information related to early infection. We investigated mRNA and miRNA transcription in Penaeus vannamei over 36 h following infection. Over this time course, 6192 transcripts and 27 miRNAs were differentially expressed—with limited differential expression from 3–12 h post injection (hpi) and a more significant transcriptional response associated with the onset of disease symptoms (24 hpi). During early infection, regulated processes included cytoskeletal remodelling and alterations in phagocytic activity that may assist WSSV entry and translocation, novel miRNA-induced metabolic shifts, and the downregulation of ATP-dependent proton transporter subunits that may impair cellular recycling. During later infection, uncoupling of the electron transport chain may drive cellular dysfunction and lead to high mortalities in infected penaeids. We propose that post-transcriptional silencing of the immune priming gene Dscam (downregulated following infections) by a novel shrimp miRNA (Pva-pmiR-78; upregulated) as a potential mechanism preventing future recognition of WSSV that may be suppressed in surviving shrimp. Our findings improve our understanding of WSD pathogenesis in P. vannamei and provide potential avenues for future development of prophylactics and treatments.

Evidence of Transcriptional Shutoff by Pathogenic Viral Haemorrhagic Septicaemia Virus in Rainbow Trout

The basis of pathogenicity of viral haemorrhagic septicaemia virus (VHSV) was analysed in the transcriptome of a rainbow trout cell line inoculated with pathogenic and non-pathogenic VHSV isolates. Although both VHSV isolates showed similar viral replication patterns, the number of differentially expressed genes was 42-fold higher in cells inoculated with the non-pathogenic VHSV at 3 h post inoculation (hpi). Infection with the non-pathogenic isolate resulted in Gene Ontologies (GO) enrichment of terms such as immune response, cytokine-mediated signalling pathway, regulation of translational initiation, unfolded protein binding, and protein folding, and induced an over-representation of the p53, PPAR, and TGF-β signalling pathways. Inoculation with the pathogenic isolate resulted in the GO enrichment of terms related to lipid metabolism and the salmonella infection KEGG pathway involved in the rearrangement of the cytoskeleton. Antiviral response was evident at 12hpi in cells infected with the pathogenic isolate. Overall, the data showed a delay in the response of genes involved in immune responses and viral sensing in cells inoculated with the pathogenic isolate and suggest transcriptional shutoff and immune avoidance as a critical mechanism of pathogenicity in VHSV. These pathways offer opportunities to further understand and manage VHSV pathogenicity in rainbow trout.

A Novel RNA Virus, Macrobrachium rosenbergii Golda Virus (MrGV), Linked to Mass Mortalities of the Larval Giant Freshwater Prawn in Bangladesh

Mass mortalities of the larval stage of the giant freshwater prawn, Macrobrachium rosenbergii, have been occurring in Bangladesh since 2011. Mortalities can reach 100% and have resulted in an 80% decline in the number of hatcheries actively producing M. rosenbergii. To investigate a causative agent for the mortalities, a disease challenge was carried out using infected material from a hatchery experiencing mortalities. Moribund larvae from the challenge were prepared for metatranscriptomic sequencing. De novo virus assembly revealed a 29 kb single‑stranded positive-sense RNA virus with similarities in key protein motif sequences to yellow head virus (YHV), an RNA virus that causes mass mortalities in marine shrimp aquaculture, and other viruses in the Nidovirales order. Primers were designed against the novel virus and used to screen cDNA from larvae sampled from hatcheries in the South of Bangladesh from two consecutive years. Larvae from all hatcheries screened from both years were positive by PCR for the novel virus, including larvae from a hatchery that at the point of sampling appeared healthy, but later experienced mortalities. These screens suggest that the virus is widespread in M. rosenbergii hatchery culture in southern Bangladesh, and that early detection of the virus can be achieved by PCR. The hypothesised protein motifs of Macrobrachium rosenbergii golda virus (MrGV) suggest that it is likely to be a new species within the Nidovirales order. Biosecurity measures should be taken in order to mitigate global spread through the movement of post-larvae within and between countries, which has previously been linked to other virus outbreaks in crustacean aquaculture.

Whole Genome Sequencing of Hepatitis A Virus Using a PCR-Free Single-Molecule Nanopore Sequencing Approach

Hepatitis A virus (HAV) is one of the most common causes of acute viral hepatitis in humans. Although HAV has a relatively small genome, there are several factors limiting whole genome sequencing such as PCR amplification artefacts and ambiguities in de novo assembly. The recently developed Oxford Nanopore technologies (ONT) allows single-molecule sequencing of long-size fragments of DNA or RNA using PCR-free strategies. We have sequenced the whole genome of HAV using a PCR-free approach by direct reverse-transcribed sequencing. We were able to sequence HAV cDNA and obtain reads over 7 kilobases in length containing almost the whole genome of the virus. The comparison of these raw long nanopore reads with the HAV reference wild type revealed a nucleotide sequence identity between 81.1 and 96.6%. By de novo assembly of all HAV reads we obtained a consensus sequence of 7362 bases, with a nucleotide sequence identity of 99.0% with the genome of the HAV strain pHM175/18f. When the assembly was performed using as reference the HAV strain pHM175/18f a consensus with a sequence similarity of 99.8 % was obtained. We have also used an ONT amplicon-based assay to sequence two fragments of the VP3 and VP1 regions which showed a sequence similarity of 100% with matching regions of the consensus sequence obtained using the direct cDNA sequencing approach. This study showed the applicability of ONT sequencing technologies to obtain the whole genome of HAV by direct cDNA nanopore sequencing, highlighting the utility of this PCR-free approach for HAV characterization and potentially other viruses of the Picornaviridae family.

How do abiotic environmental conditions influence shrimp susceptibility to disease? A critical analysis focussed on White Spot Disease

White Spot Syndrome Virus (WSSV) causes White Spot Disease (WSD) and is historically the most devastating disease in the shrimp industry. Global losses from this disease have previously exceeded $3 bn annually, having a major impact on a global industry worth US$19 bn per annum. Shrimp are cultured predominantly in enclosed ponds that are subject to considerable fluctuations in abiotic conditions and WSD outbreaks are increasingly linked to periods of extreme weather, which may cause major fluctuations in pond culture conditions. Combined with the intensity of production in these systems, the resulting suboptimal physicochemical conditions have a major bearing on the susceptibility of shrimp to infection and disease. Current knowledge indicates that pond temperature and salinity are major factors determining outbreak severity. WSSV appears to be most virulent in water temperatures between 25 and 28 °C and salinities far removed from the isoosmotic point of shrimp. Elevated temperatures (>30 °C) may protect against WSD, depending on the stage of infection, however the mechanisms mediating this effect have not been well established. Other factors relating to water quality that may play key roles in determining outbreak severity include dissolved oxygen concentration, nitrogenous compound concentration, partial pressure of carbon dioxide and pH, but data on their impacts on WSSV susceptibility in cultured shrimps is scarce. This illustrates a major research gap in our understanding of the influence of environmental conditions on disease. For example, it is not clear whether temperature manipulations can be used effectively to prevent or mitigate WSD in cultured shrimp. Therefore, developing our understanding of the impact of environmental conditions on shrimp susceptibility to WSSV may provide insight for WSD mitigation when, even after decades of research, there is no effective practical prophylaxis or treatment.

The Segment Matters: Probable Reassortment of Tilapia Lake Virus (TiLV) Complicates Phylogenetic Analysis and Inference of Geographical Origin of New Isolate from Bangladesh

Tilapia lake virus (TiLV), a negative sense RNA virus with a 10 segment genome, is an emerging threat to tilapia aquaculture worldwide, with outbreaks causing over 90% mortality reported on several continents since 2014. Following a severe tilapia mortality event in July 2017, we confirmed the presence of TiLV in Bangladesh and obtained the near-complete genome of this isolate, BD-2017. Phylogenetic analysis of the concatenated 10 segment coding regions placed BD-2017 in a clade with the two isolates from Thailand, separate from the Israeli and South American isolates. However, phylogenetic analysis of individual segments gave conflicting results, sometimes clustering BD-2017 with one of the Israeli isolates, and splitting pairs of isolates from the same region. By comparing patterns of topological difference among segments of quartets of isolates, we showed that TiLV likely has a history of reassortment. Segments 5 and 6, in particular, appear to have undergone a relatively recent reassortment event involving Ecuador isolate EC-2012 and Israel isolate Til-4-2011. The phylogeny of TiLV isolates therefore depends on the segment sequenced. Our findings illustrate the need to exercise caution when using phylogenetic analysis to infer geographic origin and track the movement of TiLV, and we recommend using whole genomes wherever possible.

Computational methods for 16S metabarcoding studies using Nanopore sequencing data

Assessment of bacterial diversity through sequencing of 16S ribosomal RNA (16S rRNA) genes has been an approach widely used in environmental microbiology, particularly since the advent of high-throughput sequencing technologies. An additional innovation introduced by these technologies was the need of developing new strategies to manage and investigate the massive amount of sequencing data generated. This situation stimulated the rapid expansion of the field of bioinformatics with the release of new tools to be applied to the downstream analysis and interpretation of sequencing data mainly generated using Illumina technology. In recent years, a third generation of sequencing technologies has been developed and have been applied in parallel and complementarily to the former sequencing strategies. In particular, Oxford Nanopore Technologies (ONT) introduced nanopore sequencing which has become very popular among molecular ecologists. Nanopore technology offers a low price, portability and fast sequencing throughput. This powerful technology has been recently tested for 16S rRNA analyses showing promising results. However, compared with previous technologies, there is a scarcity of bioinformatic tools and protocols designed specifically for the analysis of Nanopore 16S sequences. Due its notable characteristics, researchers have recently started performing assessments regarding the suitability MinION on 16S rRNA sequencing studies, and have obtained remarkable results. Here we present a review of the state-of-the-art of MinION technology applied to microbiome studies, the current possible application and main challenges for its use on 16S rRNA metabarcoding.

Heart Regeneration in the Mexican Cavefish

Although Astyanax mexicanus surface fish regenerate their hearts after injury, their Pachón cave-dwelling counterparts cannot and, instead, form a permanent fibrotic scar, similar to the human heart. Myocardial proliferation peaks at similar levels in both surface fish and Pachón 1 week after injury. However, in Pachón, this peak coincides with a strong scarring and immune response, and ultimately, cavefish cardiomyocytes fail to replace the scar. We identified lrrc10 to be upregulated in surface fish compared with Pachón after injury. Similar to cavefish, knockout of lrrc10 in zebrafish impairs heart regeneration without affecting wound cardiomyocyte proliferation. Furthermore, using quantitative trait locus (QTL) analysis, we have linked the degree of heart regeneration to three loci in the genome, identifying candidate genes fundamental to the difference between scarring and regeneration. Our study provides evidence that successful heart regeneration entails a delicate interplay between cardiomyocyte proliferation and scarring.

•

Astyanax mexicanus surface fish regenerate their hearts, but Pachón cavefish cannot

•

Successful regeneration is a delicate interplay between proliferation and scarring

•

lrrc10 is required for heart regeneration

•

Using QTL analysis, we have identified three genomic loci linked to heart regeneration

De novo transcriptome assembly of the Qatari pearl oyster Pinctada imbricata radiata

The pearl oyster Pinctada imbricata radiata is an iconic species in Qatar, representing an integral part of the nation's cultural heritage and one of the main economic foundations upon which the nation developed. During the early part of the 20th century, nearly half the Qatar population was involved in the pearl oyster industry. However, the fishery has undergone steady decline since the 1930s, and the species is now under threat due to multiple confounding pressures. This manuscript presents the first de novo transcriptome of the Qatari pearl oyster assembled into 30,739 non-redundant coding sequences and with a BUSCO completeness score of 98.4%. Analysis of the transcriptome reveals the close evolutionary distance to the conspecific animal Pinctada imbricata fucata but also highlights differences in immune genes and the presence of distinctive transposon families, suggesting recent adaptive divergence. This data is made available for all to utilise in future studies on the species.

Health assessment of the cleaner fish ballan wrasse Labrus bergylta from the British south-west coast

Wild-caught ballan wrasse Labrus bergylta are translocated en masse from the British south-west coast to Scotland for use as cleaner fish to tackle Atlantic salmon Salmo salar sea lice infestations; however, very little is known about the background health status of this species. This is the first health assessment of wild ballan wrasse from the British south-west. Wild-caught ballan wrasse (n = 75) from coastal populations off Dorset and Cornwall were subjected to a full health screen for viral, bacterial and parasitic infections and associated pathology. A range of metazoan and protozoan parasites were observed in histological sections, including copepods (sea lice Caligus centrodonti), nematodes, cestodes, digenean metacercariae, Cryptocaryon-like ciliates and an intestinal coccidian (Eimeria sp.) observed in 26.6% of the samples. The mycoplasma Acholeplasma laidlawii was associated with cytopathic effect in cell culture inoculated with tissue homogenates. The opportunistic pathogen Photobacterium damselae damselae was isolated from a single fish with a systemic infection. The isolate was confirmed to possess the virulence factors hlyAch and plpV, previously associated with cell toxicity and pathogenicity to fish. There are no immediate concerns for the continued mass translation of ballan wrasse, however careful monitoring of the population is recommended.

The skin immune response of rainbow trout, Oncorhynchus mykiss (Walbaum), associated with puffy skin disease (PSD)

Puffy skin disease (PSD) is an emerging skin condition which affects rainbow trout, Oncorhynchus mykiss (Walbaum). The transmission pattern of PSD suggests an infectious aetiology, however, the actual causative infectious agent(s) remain(s) unknown. In the present study, the rainbow trout epidermal immune response to PSD was characterised. Skin samples from infected fish were analysed and classified as mild, moderate or severe PSD by gross pathology and histological assessment. The level of expression of 26 immune-associated genes including cytokines, immunoglobulins and cell markers were examined by TaqMan qPCR assays. A significant up-regulation of the gene expression of C3, lysozyme, IL-1β and T-bet and down-regulation of TGFβ and TLR3 was observed in PSD fish compared to control fish. MHCI gene expression was up-regulated only in severe PSD lesions. Histological examinations of the epidermis showed a significant increase in the number of eosinophil cells and dendritic melanocytes in PSD fish. In severe lesions, mild diffuse lymphocyte infiltration was observed. IgT and CD8 positive cells were detected locally in the skin of PSD fish by in situ hybridisation (ISH), however, the gene expression of those genes was not different from control fish. Total IgM in serum of diseased animals was not different from control fish, measured by a sandwich ELISA, nor was significant up regulation of IgM gene expression in PSD lesions observed. Taken together, these results show activation of the complement pathway, up-regulation of a Th17 type response and eosinophilia during PSD. This is typical of a response to extracellular pathogens (i.e. bacteria and parasites) and allergens, commonly associated with acute dermatitis.

Membrane Trafficking Modulation during Entamoeba Encystation

Entamoeba histolytica is an intestinal parasite that infects 50-100 million people and causes up to 55,000 deaths annually. The transmissive form of E. histolytica is the cyst, with a single infected individual passing up to 45 million cysts per day, making cyst production an attractive target for infection control. Lectins and chitin are secreted to form the cyst wall, although little is known about the underlying membrane trafficking processes supporting encystation. As E. histolytica does not readily form cysts in vitro, we assessed membrane trafficking gene expression during encystation in the closely related model Entamoeba invadens. Genes involved in secretion are up-regulated during cyst formation, as are some trans-Golgi network-to-endosome trafficking genes. Furthermore, endocytic and general trafficking genes are up-regulated in the mature cyst, potentially preserved as mRNA in preparation for excystation. Two divergent dynamin-related proteins found in Entamoeba are predominantly expressed during cyst formation. Phylogenetic analyses indicate that they are paralogous to, but quite distinct from, classical dynamins found in human, suggesting that they may be potential drug targets to block encystation. The membrane-trafficking machinery is clearly regulated during encystation, providing an additional facet to understanding this crucial parasitic process.

Bioavailability and Kidney Responses to Diclofenac in the Fathead Minnow (Pimephales promelas)

Diclofenac is one of the most widely prescribed nonsteroidal anti-inflammatory drugs worldwide. It is frequently detected in surface waters; however, whether this pharmaceutical poses a risk to aquatic organisms is debated. Here we quantified the uptake of diclofenac by the fathead minnow (Pimephales promelas) following aqueous exposure (0.2-25.0 μg L^-1) for 21 days, and evaluated the tissue and biomolecular responses in the kidney. Diclofenac accumulated in a concentration- and time-dependent manner in the plasma of exposed fish. The highest plasma concentration observed (for fish exposed to 25 μg L^-1 diclofenac) was within the therapeutic range for humans. There was a strong positive correlation between exposure concentration and the number of developing nephrons observed in the posterior kidney. Diclofenac was not found to modulate the expression of genes in the kidney associated with its primary mode of action in mammals (prostaglandin-endoperoxide synthases) but modulated genes associated with kidney repair and regeneration. There were no significant adverse effects following 21 days exposure to concentrations typical of surface waters. The combination of diclofenac's uptake potential, effects on kidney nephrons and relatively small safety margin for some surface waters may warrant a longer term chronic health effects analysis for diclofenac in fish.

Puffy Skin Disease Is an Emerging Transmissible Condition in Rainbow Trout Oncorhynchus mykiss Walbaum

The transmission of puffy skin disease (PSD) to rainbow trout Oncorhynchus mykiss Walbaum was tested in the laboratory by conducting co-habitation challenges with puffy skin (PS)-affected fish (Trojans) collected from the field. Two separate challenges were conducted using Trojans sourced from two different sites and diploid (first trial) or triploid (second trial) naïve fish. PSD-specific clinical signs were observed in both groups of naïve fish, with 66% of the fish sampled during the challenges showing signs of varying severity. The first clinical features of PSD were presented as white oval skin patches on one or both flanks 15–21 days post-challenge (dpc). The extent of the lesions ranged from 10 to 90% of the body surface, depending on the severity of the lesion. Both the severity and number of affected fish increased during the challenge. Macroscopically, oedema of the skin and multifocal petechial haemorrhaging were observed towards the end of the trials. Abnormal fish behaviour consisting of “flashing” and excessive mucous production was noted from 15 dpc onwards. Fish with severe PSD lesions also displayed inappetence and associated emaciation. Rodlet cells were observed in 41% of the fresh skin scrapes analysed from the second trial. Histologically epidermal oedema was observed in 31% of the naive fish showing gross pathology, with additional 12% displaying epidermal hyperplasia, mostly observed at the end of the challenge. Other concomitant features of the PSD lesions in challenged fish were epithelial erosion and sloughing, and occasionally mild or focal inflammation. No consistent pathology of internal organs was observed. The parasites Ichthyophthirius multifiliis and Ichthyobodo necator were observed in skin samples of a proportion of naïve challenged fish and in Trojans but not in control fish. The presence of these and other known fish pathogens in the skin of PSD-fish was confirmed by high-throughput sequencing analysis. In summary, we have demonstrated that PSD is a transmissible condition. However, even though a number of known fish pathogens were identified in the skin tissues of PSD-fish, the actual causative infectious agent(s) remain(s) unknown.

In vivo virulence of viral haemorrhagic septicaemia virus (VHSV) in rainbow trout Oncorhynchus mykiss correlates inversely with in vitro Mx gene expression

The in vitro replication of viral haemorrhagic septicaemia virus (VHSV) isolates from each VHSV genotype and the associated cellular host Mx gene expression were analysed. All the isolates were able to infect RTG-2 cells and induce increased Mx gene expression (generic assay detecting isoforms 1 and 3 [Mx1/3]). A trout pathogenic, genotype Ia isolate (J167), showing high replication in RTG-2 cells (by infective titre and N gene expression) induced lower Mx1/3 gene expression than observed in VHSV isolates known to be non-pathogenic to rainbow trout: 96-43/8, 96-43/10 (Ib); 1p49, 1p53 (II); and MI03 (IVb). Paired co-inoculation assays were analysed using equal number of plaque forming units per ml (PFU) of J167 (Ia genotype) with other less pathogenic VHSV genotypes. In these co-inoculations, the Mx1/3 gene expression was significantly lower than for the non-pathogenic isolate alone. Of the three rainbow trout Mx isoforms, J167 did not induce Mx1 up-regulation in RTG-2 or RTgill-W1 cells. Co-inoculating isolates resulted in greater inhibition of Mx in both rainbow trout cell lines studied. Up-regulation of sea bream Mx in SAF-1 cells induced by 96-43/8 was also lower in co-inoculation assays with J167. The RTG-P1 cell line, expressing luciferase under the control of the interferon-induced Mx rainbow trout gene promoter, showed low luciferase activity when inoculated with pathogenic strains: J167, DK-5131 (Ic), NO-A-163/68 (Id), TR-206239-1, TR-22207111 (Ie), 99-292 (IVa), and CA-NB00-01 (IVc). Co-inoculation assays showed a J167-dose dependent inhibition of the luciferase activity. The data suggest that virulent VHSV isolates may interfere in the interferon pathways, potentially determining higher pathogenicity.

Molecular Mechanisms of White Spot Syndrome Virus Infection and Perspectives on Treatments

Since its emergence in the 1990s, White Spot Disease (WSD) has had major economic and societal impact in the crustacean aquaculture sector. Over the years shrimp farming alone has experienced billion dollar losses through WSD. The disease is caused by the White Spot Syndrome Virus (WSSV), a large dsDNA virus and the only member of the Nimaviridae family. Susceptibility to WSSV in a wide range of crustacean hosts makes it a major risk factor in the translocation of live animals and in commodity products. Currently there are no effective treatments for this disease. Understanding the molecular basis of disease processes has contributed significantly to the treatment of many human and animal pathogens, and with a similar aim considerable efforts have been directed towards understanding host-pathogen molecular interactions for WSD. Work on the molecular mechanisms of pathogenesis in aquatic crustaceans has been restricted by a lack of sequenced and annotated genomes for host species. Nevertheless, some of the key host-pathogen interactions have been established: between viral envelope proteins and host cell receptors at initiation of infection, involvement of various immune system pathways in response to WSSV, and the roles of various host and virus miRNAs in mitigation or progression of disease. Despite these advances, many fundamental knowledge gaps remain; for example, the roles of the majority of WSSV proteins are still unknown. In this review we assess current knowledge of how WSSV infects and replicates in its host, and critique strategies for WSD treatment.

De novo assembly of the Carcinus maenas transcriptome and characterization of innate immune system pathways

Background

The European shore crab, Carcinus maenas, is used widely in biomonitoring, ecotoxicology and for studies into host-pathogen interactions. It is also an important invasive species in numerous global locations. However, the genomic resources for this organism are still sparse, limiting research progress in these fields. To address this resource shortfall we produced a C. maenas transcriptome, enabled by the progress in next-generation sequencing technologies, and applied this to assemble information on the innate immune system in this species.

Results

We isolated and pooled RNA for twelve different tissues and organs from C. maenas individuals and sequenced the RNA using next generation sequencing on an Illumina HiSeq 2500 platform. After de novo assembly a transcriptome was generated encompassing 212,427 transcripts (153,699 loci). The transcripts were filtered, annotated and characterised using a variety of tools (including BLAST, MEGAN and RSEM) and databases (including NCBI, Gene Ontology and KEGG). There were differential patterns of expression for between 1,223 and 2,741 transcripts across tissues and organs with over-represented Gene Ontology terms relating to their specific function. Based on sequence homology to immune system components in other organisms, we show both the presence of transcripts for a series of known pathogen recognition receptors and response proteins that form part of the innate immune system, and transcripts representing the RNAi, Toll-like receptor signalling, IMD and JAK/STAT pathways.

Conclusions

We have produced an assembled transcriptome for C. maenas that provides a significant molecular resource for wide ranging studies in this species. Analysis of the transcriptome has revealed the presence of a series of known targets and functional pathways that form part of their innate immune system and illustrate tissue specific differences in their expression patterns.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1667-1) contains supplementary material, which is available to authorized users.

Bmp suppression in mangrove killifish embryos causes a split in the body axis

Bone morphogenetic proteins (Bmp) are major players in the formation of the vertebrate body plan due to their crucial role in patterning of the dorsal-ventral (DV) axis. Despite the highly conserved nature of Bmp signalling in vertebrates, the consequences of changing this pathway can be species-specific. Here, we report that Bmp plays an important role in epiboly, yolk syncytial layer (YSL) movements, and anterior-posterior (AP) axis formation in embryos of the self-fertilizing mangrove killifish, Kryptolebias marmoratus. Stage and dose specific exposures of embryos to the Bmp inhibitor dorsomorphin (DM) produced three distinctive morphologies, with the most extreme condition creating the splitbody phenotype, characterised by an extremely short AP axis where the neural tube, somites, and notochord were bilaterally split. In addition, parts of caudal neural tissues were separated from the main body and formed cell islands in the posterior region of the embryo. This splitbody phenotype, which has not been reported in other animals, shows that modification of Bmp may lead to significantly different consequences during development in other vertebrate species.

Molecular mechanisms of toxicity of silver nanoparticles in zebrafish embryos

Silver nanoparticles cause toxicity in exposed organisms and are an environmental health concern. The mechanisms of silver nanoparticle toxicity, however, remain unclear. We examined the effects of exposure to silver in nano-, bulk-, and ionic forms on zebrafish embryos (Danio rerio) using a Next Generation Sequencing approach in an Illumina platform (High-Throughput SuperSAGE). Significant alterations in gene expression were found for all treatments and many of the gene pathways affected, most notably those associated with oxidative phosphorylation and protein synthesis, overlapped strongly between the three treatments indicating similar mechanisms of toxicity for the three forms of silver studied. Changes in oxidative phosphorylation indicated a down-regulation of this pathway at 24 h of exposure, but with a recovery at 48 h. This finding was consistent with a dose-dependent decrease in oxygen consumption at 24 h, but not at 48 h, following exposure to silver ions. Overall, our data provide support for the hypothesis that the toxicity caused by silver nanoparticles is principally associated with bioavailable silver ions in exposed zebrafish embryos. These findings are important in the evaluation of the risk that silver particles may pose to exposed vertebrate organisms.

Bioavailability of nanoscale metal oxides TiO(2), CeO(2), and ZnO to fish

Nanoparticles (NPs) are reported to be a potential environmental health hazard. For organisms living in the aquatic environment, there is uncertainty on exposure because of a lack of understanding and data regarding the fate, behavior, and bioavailability of the nanomaterials in the water column. This paper reports on a series of integrative biological and physicochemical studies on the uptake of unmodified commercial nanoscale metal oxides, zinc oxide (ZnO), cerium dioxide (CeO(2)), and titanium dioxide (TiO(2)), from the water and diet to determine their potential ecotoxicological impacts on fish as a function of concentration. Particle characterizations were performed and tissue concentrations were measured by a wide range of analytical methods. Definitive uptake from the water column and localization of TiO(2) NPs in gills was demonstrated for the first time by use of coherent anti-Stokes Raman scattering (CARS) microscopy. Significant uptake of nanomaterials was found only for cerium in the liver of zebrafish exposed via the water and ionic titanium in the gut of trout exposed via the diet. For the aqueous exposures undertaken, formation of large NP aggregates (up to 3 mum) occurred and it is likely that this resulted in limited bioavailability of the unmodified metal oxide NPs in fish.

Identifying health impacts of exposure to copper using transcriptomics and metabolomics in a fish model

Copper (Cu) is a micronutrient essential for the biochemical functioning of numerous processes in vertebrates but is also often present in the aquatic environment at concentrations able to cause adverse health effects in aquatic organisms. This study investigated the signaling pathways mediating the effects of exposure to Cu using a toxicogenomic approach in a fish model, the stickleback ( Gasterosteus aculeatus ). Freshwater-acclimated male fish were exposed via the water to Cu, including at environmentally relevant concentrations (3.2-128 microg of Cu/L for 4 days), and the biological responses explored through analyses of the hepatic transcriptome and metabolome and phenotypic end points, including assessment of DNA damage in blood cells. The Cu exposures resulted in DNA strand breaks in blood cells at all exposure concentrations and alterations in hepatic gene expression and metabolite concentrations in a concentration-dependent manner (from 10 microg of Cu/L). Genes associated with the cholesterol biosynthesis pathway were significantly over-represented and consistently down-regulated (at 128 microg of Cu/L), similar to that occurring in a mouse model for Wilson's disease. Additionally, inductions in metallothionein and catalase were also observed. The concentrations of NAD(+) and lactate increased significantly with the Cu exposure, consistent with a shift toward anaerobic metabolism, and these aligned closely with changes observed in gene expression. The pathways of Cu toxicity identified in our study support the conserved mechanisms of Cu toxicity from lower vertebrates to mammals, provide novel insights into the deleterious effects of Cu in fish, and further demonstrate the utility of fish as environmental sentinels for chemical impacts on both environmental and human health.

The kisspeptin/gonadotropin-releasing hormone pathway and molecular signaling of puberty in fish

The mechanisms underlying the initiation of puberty in fish are poorly understood, and whether the Kiss1 receptor (Kiss1r; previously designated G protein-coupled receptor 54; GPR54) and its ligands, kisspeptins, play a significant role, as has been established in mammals, is not yet known. We determined (via real-time PCR) temporal patterns of expression in the brain of kiss1r, gnrh2, and gnrh3 and a suite of related genes in the hypothalamo-pituitary-gonadal (HPG) axis and analyzed them against the timing of gonadal germ cell development in male and female fathead minnow (Pimephales promelas). Full- or partial-length cDNAs for kiss1r (736 bp), gnrh2 (698 bp), and gnrh3 (804 bp) cloned from fathead minnow were found to be expressed only in the brain, testis, and ovary of adult fish. Localization of kiss1r, gnrh2, and gnrh3 within the brain provided evidence for their physiological roles and a likely hypophysiotropic role for GnRH3 in this species (which, like other cyprinids, does not appear to express gnrh1). In both sexes, kiss1r expression in the brain increased at the onset of puberty and reached maximal expression in males when spermatagonia type B appeared in the testis and in females when cortical alveolus-stage oocytes first appeared in the ovary, the timings of which differed for the two sexes. However, kiss1r expression was considerably lower during more advanced stages of spermatogenesis and oogenesis. The expression of kiss1r closely aligned with that of the gnrh genes (gnrh3 in particular), suggesting the Kiss1r/kisspeptin system in fish has a similar role in puberty to that occurring in mammals, and this hypothesis was supported by the induction of gnrh3 (2.25-fold) and kiss1r (1.5-fold) in early-mid pubertal fish injected with mammalian kisspeptin-10 (2 nmol/g wet weight). An intriguing finding, and contrasting that in mammals, was an elevated expression of esr1, ar, and cyp19a2 (genes involved in sex steroid signaling) in the brain at the onset of puberty, and in females slightly in advance of the elevation in the expression of kiss1r.

Development and validation of a direct homologous quantitative sandwich ELISA for fathead minnow (Pimephales promelas) vitellogenin

Vitellogenin (Vtg) is an established and sensitive endpoint for analysis of exposure to (anti-)oestrogens and their mimics in fish [Sumpter, J.P., 1995. Feminized responses in fish to environmental estrogens. Toxicol. Lett. 82, 737-742; Arukwe, A., Goksøyr, A., 2003. Eggshell and egg yolk proteins in fish: hepatic proteins for the next generation: oogenetic, population, and evolutionary implications of endocrine disruption. Comp. Hepatol. 2, 4. ]. In some instances, links have been drawn between high level induction of Vtg and adverse health effects in fish [Herman, R.L., Kincaide, H.L., 1988. Pathological effects of orally administered estradiol to rainbow trout. Aquaculture 72, 165-172; Schwaiger, J., Spieser, O.H., Bauer, C., Ferling, H., Mallow, U., Kalbfus, W., Negele, R.D., 2000. Chronic toxicity of nonylphenol and ethinyloestraiol: haematological and histopathological effects in juvenile common carp (Cyprinus carpio). Aquat. Toxicol. 51, 69-78]. The widespread use of Vtg as a biomarker has led to the development of a variety of assays to quantitatively measure Vtg concentrations in tissue samples from fish, and hence a need for a standardization of the performance criteria and validation of such assays [Goksøyr, A., Eidem, J.K., Kristiansen, S.I., Nilsen, B.M., 2003. On the need for a standardized set-up for validation studies of fish vitellogenin assays as an endpoint in endocrine disruptor testing and screening-a proposal. ]. One of the most popular test fish species for assessing chemical effects is the fathead minnow (Pimephales promelas), which is now used widely for studies into endocrine disruption [Panter, G.H., Hutchinson, T.H., Lange, R., Lye, C.M., Sumpter, J.P., Zerulla, M., Tyler, C.R., 2002. Utility of a juvenile fathead minnow screening assay for detecting (anti)estrogenic substances. Environ. Toxicol. Chem. 21, 319-326; Hutchinson, T.H., Yokota, H., Hagino, S., Ozato, K., 2003. Development of fish tests for endocrine disruptors. Pure Appl. Chem. 75, 2343-2353]. This paper describes the development and validation of a new, homologous enzyme-linked immunosorbent assay (ELISA) for quantification of Vtg in this fish species.

COMPRENDO: Focus and approach

Tens of thousands of man-made chemicals are in regular use and discharged into the environment. Many of them are known to interfere with the hormonal systems in humans and wildlife. Given the complexity of endocrine systems, there are many ways in which endocrine-disrupting chemicals (EDCs) can affect the body's signaling system, and this makes unraveling the mechanisms of action of these chemicals difficult. A major concern is that some of these EDCs appear to be biologically active at extremely low concentrations. There is growing evidence to indicate that the guiding principle of traditional toxicology that "the dose makes the poison" may not always be the case because some EDCs do not induce the classical dose-response relationships. The European Union project COMPRENDO (Comparative Research on Endocrine Disrupters--Phylogenetic Approach and Common Principles focussing on Androgenic/Antiandrogenic Compounds) therefore aims to develop an understanding of potential health problems posed by androgenic and antiandrogenic compounds (AACs) to wildlife and humans by focusing on the commonalities and differences in responses to AACs across the animal kingdom (from invertebrates to vertebrates) .

Predicted exposures to steroid estrogens in U.K. rivers correlate with widespread sexual disruption in wild fish populations

Steroidal estrogens, originating principally from human excretion, are likely to play a major role in causing widespread endocrine disruption in wild populations of the roach (Rutilus rutilus), a common cyprinid fish, in rivers contaminated by treated sewage effluents. Given the extent of this problem, risk assessment models are needed to predict the location and severity of endocrine disruption in river catchments and to identify areas where regulation of sewage discharges to remove these contaminants is necessary. In this study we attempted to correlate the extent of endocrine disruption in roach in British rivers, with their predicted exposure to steroid estrogens derived from the human population. The predictions of steroid estrogen exposure at each river site were determined by combining the modeled concentrations of the individual steroid estrogens [17beta-estradiol (E2), estrone (E1), and 17alpha-ethinylestradiol (EE2)] in each sewage effluent with their predicted dilution in the immediate receiving water. This model was applied to 45 sites on 39 rivers throughout the United Kingdom. Each site studied was then categorized as either high, medium, or low "risk" on the basis of the assumed additive potency of the three steroid estrogens calculated from data derived from published studies in various cyprinid fish species. We sampled 1,438 wild roach from the predicted high-, medium-, and low-risk river sites and examined them for evidence and severity of endocrine disruption. Both the incidence and the severity of intersex in wild roach were significantly correlated with the predicted concentrations of the natural estrogens (E1 and E2) and the synthetic contraceptive pill estrogen (EE2) present. Predicted steroid estrogen exposure was, however, less well correlated with the plasma vitellogenin concentration measured in the same fish. Moreover, we found no correlation between any of the end points measured in the roach and the proportion of industrial effluents entering the rivers we studied. Overall, our results provide further and substantive evidence to support the hypothesis that steroidal estrogens play a major role in causing intersex in wild freshwater fish in rivers in the United Kingdom and clearly show that the location and severity of these endocrine-disrupting effects can be predicted.

Endocrine (sexual) disruption is not a prominent feature in the pike (Esox lucius), a top predator, living in English waters

The high incidence of intersex roach (Rutilus rutilus) in some United Kingdom rivers that has been associated with exposure to sewage treatment works (STWs) effluent led us to hypothesize that top predator fish also may be affected by estrogenic chemicals, because they are likely to bioaccumulate lipophilic compounds through a predator-prey relationship. To investigate this possibility, pike (Esox lucius) were sampled both upstream and downstream of STWs and then examined for total estrogenic activity of their bile, as measured using a yeast-based estrogen assay to determine the degree of recent exposure of the pike to estrogens and vitellogenin induction, and for possible disruption of sexual development, as measured using histological analysis of the gonads. No evidence of severe disruption was found in the sampled fish, which came from 16 sampling sites that were representative of English rivers. However, 14% of pike were intersex, of which 15 of 16 showed patches of male germ cells among predominantly female gonadal tissue. The incidence of masculinization was independent of whether the pike had been sampled upstream or downstream of STWs. Although pike are gonochoristic, it is not known if this masculinization of presumptive female pike is normal or, instead, indicative of endocrine disruption. Vitellogenin concentrations were not elevated in male pike at sites either upstream or downstream of STWs. The results suggest that sexual disruption is not common in pike, a fish at the top of the food chain in the freshwaters of England.

Effects of atrazine on sex steroid dynamics, plasma vitellogenin concentration and gonad development in adult goldfish (Carassius auratus)

Sexually mature goldfish (Carassius auratus) of both sexes were exposed to two doses (100 and 1000 microg/l) of the widely used herbicide atrazine (2-chloro-4-ethylamino-6-isopropylamino-s-triazine) for a period of 21 days and effects on the concentrations of gonad and plasma sex steroids (testosterone (T), 17beta-estradiol (E2) and 11-ketotestosterone (11-KT)), plasma vitellogenin (VTG) and gonad histo-morphology assessed. Atrazine did not show any obvious estrogenic effect in males, as determined by a lack of vitellogenin induction. There were, however, effects of atrazine on plasma androgen concentrations (androgen dynamics) and tissue (plasma and gonad) estrogen concentrations in male goldfish; exposure to 1000 microg/l atrazine induced suppression in both testosterone and 11-ketotestosterone, and resulted in elevated 17beta-estradiol, after 21 day of exposure. Further, these suppressive effects on plasma androgens and the induction in estrogen were dose- and time-related. The highest atrazine exposure dose induced structural disruption in the testis and both 100 and 1000 microg/l induced elevated levels of atresia in ovaries.

Window of sensitivity for the estrogenic effects of ethinylestradiol in early life-stages of fathead minnow, Pimephales promelas

Sexual differentiation in fish occurs after hatching during early life-stages and is believed to be a time when the gonad has a heightened sensitivity to disruption by chemicals that mimic hormones. In this study fathead minnows (Pimephales promelas) were exposed to an environmentally relevant concentration of ethinylestradiol (EE2) for short intervals in fish early life-stages and vitellogenic and gonadal responses were measured at 30 and 100 dph (sexual maturity), respectively. All EE2 exposure regimes induced vitellogenin (VTG) synthesis and disruption in duct development (a feminization) in males, with a window of enhanced sensitivity between 10 and 15 dph (where 60% of the males had feminized ducts). There was an altered pattern in sex cell development in males (inhibition of spermatogenesis) in the solvent controls (ethanol 0.1 ml/l) and all EE2 treatments when compared with the dilution water controls. Furthermore, fewer spermatozoa were observed in the testis of males exposed to EE2 from 15 to 20 dph and fertilized eggs (<24 h post-fertilization)-20 dph, compared with both the solvent and dilution water controls. These data show that short exposures of embryos/very early life-stage fathead minnows to an environmentally relevant concentration of EE2 lead to alterations in gonadal development that potentially have reproductive consequences and thus population level effects.

Monoclonal antibody enzyme-linked immunosorbent assay to quantify vitellogenin for studies on environmental estrogens in the rainbow trout (Oncorhynchus mykiss)

Vitellogenin (VTG) induction has proved to be a valuable biomarker for assessing exposure to environmental estrogens in fish. The widespread use of VTG in this regard has lead to the need for standardized assays to quantify VTG, and monoclonal antibodies have the potential to help accomplish this. A VTG enzyme-linked immunosorbent assay (ELISA) was developed using a monoclonal antibody prepared against Atlantic salmon (Salmo salar) VTG (MAb BN-5) and its ability to quantify VTG in the rainbow trout (Oncorhynchus mykiss) compared with a rainbow trout vitellogenin (rt-VTG) ELISA that employed homologous polyclonal antibodies (PAb). In routine protocols, the working range of the homologous rt-PAb VTG ELISA was between 9 ng/ml and 70 ng/ml (80- 20% relative maximum binding [B/Bo]) with a 50% B/Bo of 25+/-0.9 ng/ml and inter- and intraassay variations at 50% B/Bo of 7% (n = 7) and 8% (n = 15), respectively. The working range of the MAb BN-5 VTG ELISA was between 60 ng/ml and 850 ng/ml (80-20% B/Bo) with a 50% B/Bo of 227+/-22 ng/ml and inter- and intraassay variations at 50% B/Bo of 5% (n = 10) and 9% (n = 12), respectively. In the routine protocols, detection limits for measurement of plasma VTG in rainbow trout (at 80% B/Bo; and given the requirement to dilute plasma to a minimum of 1:10 for the assays) were 90 ng/ml for the polyclonal rt-VTG assay and approximately 600 ng/ml in the monoclonal antibody assay. In juvenile female rainbow trout exposed to a series of doses of estradiol-17beta (E2) and 4-tert nonylphenol (4-NP), there were no differences in the vitellogenic responses measured in the PAb and MAb BN-5 VTG ELISAs. The monoclonal MAb BN-5 VTG ELISA is likely to be of considerable value for studies on environmental estrogens in juvenile female rainbow trout in standardized tests.