Search

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.

Sequencing of SARS-CoV-2 in wastewater provides a key opportunity to monitor the prevalence of variants spatiotemporally, potentially facilitating their detection simultaneously with, or even prior to, observation through clinical testing. However, there are multiple sequencing methodologies available. This study aimed to evaluate the performance of alternative protocols for detecting SARS-CoV-2 variants. We tested the detection of two synthetic RNA SARS-CoV-2 genomes in a wide range of ratios and at two concentrations representative of those found in wastewater using whole-genome andSpike-gene-only protocols utilising Illumina and Oxford Nanopore platforms. We developed a Bayesian hierarchical model to determine the predicted frequencies of variants and the error surrounding our predictions. We found that most of the sequencing protocols detected polymorphic nucleotide frequencies at a level that would allow accurate determination of the variants present at higher concentrations. Most methodologies, including theSpike-only approach, could also predict variant frequencies with a degree of accuracy in low-concentration samples but, as expected, with higher error around the estimates. All methods were additionally confirmed to detect the same prevalent variants in a set of wastewater samples. Our results provide the first quantitative statistical comparison of a range of alternative methods that can be used successfully in the surveillance of SARS-CoV-2 variant frequencies from wastewater.ImpactGenetic sequencing of SARS-CoV-2 in wastewater provides an ideal system for monitoring variant frequencies in the general population. The advantages over clinical data are that it is more cost efficient and has the potential to identify new variants before clinical testing. However, to date, there has been no direct comparison to determine which sequencing methodologies perform best at identifying the presence and prevalence of variants. Our study compares seven sequencing methods to determine which performs best. We also develop a Bayesian statistical methodology to estimate the confidence around variant frequency estimates. Our results will help monitor SARS-CoV-2 variants in wastewater, and the methodology could be adapted for other disease monitoring, including future pandemics.

Patagonian toothfish (Dissostichus eleginoides) is an economically and ecologically important fish species in the family Nototheniidae, found at depths between 70 and 2,500 meters on the southern shelves and slopes around the sub-Antarctic islands of the Southern Ocean. Genomic sequence data for this species is limited. Here, we report a high-quality assembly and annotation of theD. eleginoidesgenome, 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 ofD. eleginoidesdoes 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. The high-quality genome assembly of the Patagonian toothfish will provide a valuable genetic resource for ecological and evolutionary studies on this and other closely related species.

Wastewater-based epidemiology (WBE) has been used extensively throughout the COVID-19 pandemic to detect and monitor the spread and prevalence of SARS-CoV-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 (Alpha, Beta and Delta), as well as 13 wastewater samples collected in London between the 15–18 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.

Using the self-fertilizing mangrove killifish, we characterized two mutants, shorttail (stl) and balltail (btl). These mutants showed abnormalities in the posterior notochord and muscle development. Taking advantage of a highly inbred isogenic strain of the species, we rapidly identified the mutated genes, noto and msgn1 in the stl and btl mutants, respectively, using a single lane of RNA sequencing without the need of a reference genome or genetic mapping techniques. Next, we confirmed a conserved morphant phenotype in medaka and demonstrate a crucial role of noto and msgn1 in cell sorting between the axial and paraxial part of the tail mesoderm. This novel system could substantially accelerate future small-scale forward-genetic screening and identification of mutations. Therefore, the mangrove killifish could be used as a complementary system alongside existing models for future molecular genetic studies., Summary: The mangrove killifish, a self-fertilizing vertebrate model, allows rapid identification of the mutated genes noto and msgn1 from mutant strains and their genetic roles in tail bud development.

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 hemolymph, 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 hemocytes, resulting in hemolymph 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., Ministerio de Economía y Competitividad, Grant/Award Number: BFU2017- 90114- P; Horizon 2020 Framework Programme, Grant/Award Number: 747789; Department for Environment, Food and Rural Affairs, UK Government, Grant/Award Number: FB002A and FC1214; Hezkuntza, Hizkuntza Politika Eta Kultura Saila, Eusko Jaurlaritza, Grant/Award Number: PRE_2016_2_0124; Agencia Estatal de Investigación (AEI); Fondo Europeo de Desarrollo Regional (FEDER); Ayuda Juan de la Cierva—Incorporación, Grant/Award Number: IJC2018- 036657- I

Multiple enveloped viruses with rod-shaped nucleocapsids have been described, infecting the epithelial cell nuclei within the hepatopancreas tubules of crustaceans. These bacilliform viruses share the ultrastructural characteristics of nudiviruses, a specific clade of viruses infecting arthropods. Using histology, electron microscopy and high throughput sequencing, we characterise two further bacilliform viruses from aquatic hosts, the brown shrimp (Crangon crangon) and the European shore crab (Carcinus maenas). We assembled the full double stranded, circular DNA genome sequences of these viruses (~113 and 132 kbp, respectively). Comparative genomics and phylogenetic analyses confirm that both belong within the family Nudiviridae but in separate clades representing nudiviruses found in freshwater and marine environments. We show that the three thymidine kinase (tk) genes present in all sequenced nudivirus genomes, thus far, were absent in the Crangon crangon nudivirus, suggesting there are twenty-eight core genes shared by all nudiviruses. Furthermore, the phylogenetic data no longer support the subdivision of the family Nudiviridae into four genera (Alphanudivirus to Deltanudivirus), as recently adopted by the International Committee on Taxonomy of Viruses (ICTV), but rather shows two main branches of the family that are further subdivided. Our data support a recent proposal to create two subfamilies within the family Nudiviridae, each subdivided into several genera.

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.

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 naive 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.

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.

Recent genomic and metagenomic studies have led to a dramatic expansion of the known diversity of nucleocytoplasmic large DNA viruses (NCLDVs) of eukaryotes, which include giant viruses of protists and important pathogens of vertebrates, such as poxviruses. However, the characterization of viruses from nonmodel hosts still lags behind. We sequenced the complete genomes of three viruses infecting crustaceans, the Caribbean spiny lobster, demon shrimp, and European shore crab. These viruses have the smallest genomes among the known NCLDVs, with losses of many core genes, some of which are shared with iridoviruses. The deterioration of the transcription apparatus is compatible with microscopic and ultrastructural observations indicating that these viruses replicate in the nucleus of infected cells rather than in the cytoplasm. Phylogenomic analysis indicates that these viruses are sufficiently distinct from all other NCLDVs to justify the creation of a separate family, for which we propose the name “Mininucleoviridae” (i.e., small viruses reproducing in the cell nucleus)., Panulirus argus virus 1 (PaV1) is the only known virus infecting the Caribbean spiny lobster (Panulirus argus) from the Caribbean Sea. Recently, related viruses, Dikerogammarus haemobaphes virus 1 (DhV1) and Carcinus maenas virus 1 (CmV1), have been detected in the demon shrimp (Dikerogammarus haemobaphes) and the European shore crab (Carcinus maenas), respectively, from sites in the United Kingdom. The virion morphology of these crustacean viruses is similar to that of iridoviruses. However, unlike iridoviruses and other nucleocytoplasmic large DNA viruses (NCLDVs), these viruses complete their morphogenesis in the host cell nucleus rather than in the cytoplasm. To date, these crustacean viruses have remained unclassified due to a lack of genomic data. Using an Illumina MiSeq sequencer, we sequenced the complete genomes of PaV1, CmV1, and DhV1. Comparative genome analysis shows that these crustacean virus genomes encode the 10 hallmark proteins previously described for the NCLDVs of eukaryotes, strongly suggesting that they are members of this group. With a size range of 70 to 74 kb, these are the smallest NCLDV genomes identified to date. Extensive gene loss, divergence of gene sequences, and the accumulation of low-complexity sequences reflect the extreme degradation of the genomes of these “minimal” NCLDVs rather than any direct relationship with the NCLDV ancestor. Phylogenomic analysis supports the classification of these crustacean viruses as a distinct family, “Mininucleoviridae,” within the pitho-irido-Marseille branch of the NCLDVs.

Survival of marine fishes that are exposed to elevated near-future CO2 levels is threatened by their altered responses to sensory cues. Here we demonstrate a physiological and molecular mechanism in the olfactory system that helps to explain altered behaviour under elevated CO2. We combine electrophysiology measurements and transcriptomics with behavioural experiments to investigate how elevated CO2 affects the olfactory system of European sea bass (Dicentrarchus labrax). When exposed to elevated CO2 (approximately 1,000 µatm), fish must be up to 42% closer to an odour source for detection, compared with current CO2 levels (around 400 µatm), decreasing their chances of detecting food or predators. Compromised olfaction correlated with the suppression of the transcription of genes involved in synaptic strength, cell excitability and wiring of the olfactory system in response to sustained exposure to elevated CO2 levels. Our findings complement the previously proposed impairment of γ-aminobutyric acid receptors, and indicate that both the olfactory system and central brain function are compromised by elevated CO2 levels. Marine fishes exposed to elevated CO2 levels can have altered responses to sensory cues. Research now reveals a physiological and molecular mechanism in the olfactory system that helps to explain this altered behaviour under elevated CO2.

Invasive and non-native species can pose risks to vulnerable ecosystems by co-introducing bacterial pathogens. Alternatively, co-introduced bacterial pathogens may regulate invasive population size and invasive traits. We describe a novel candidate genus and species of bacteria (‘Candidatus Aquirickettsiella gammari’) found to infect Gammarus fossarum, from its native range in Poland. The bacterium develops intracellularly within the haemocytes and cells of the musculature, hepatopancreas, connective tissues, nervous system and gonad of the host. The developmental cycle of ‘Candidatus Aquirickettsiella gammari’ includes an elementary body (496.73 nm ± 37.56 nm in length, and 176.89 nm ± 36.29 nm in width), an elliptical, condensed spherical stage (737.61 nm ± 44.51 nm in length and 300.07 nm ± 44.02 nm in width), a divisional stage, and a spherical initial body (1397.59 nm ± 21.26 nm in diameter). We provide a partial genome for ‘Candidatus Aquirickettsiella gammari’, which clades phylogenetically alongside environmental 16S rRNA sequences from aquatic habitats, and bacterial symbionts from aquatic isopods (Asellus aquaticus), grouping separately from the Rickettsiella, a genus that includes bacterial pathogens of terrestrial insects and isopods. Increased understanding of the diversity of symbionts carried by G. fossarum identifies those that might regulate host population size, or those that could pose a risk to native species in the invasive range. Identification of ‘Candidatus Aquirickettsiella gammari’ and its potential for adaptation as a biological control agent is explored.

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.

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.

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.

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.

© 2019 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. This publication was made possible by the NPRP award [NPRP9-394-1-090 “The Pearl Oyster: from national icon to guardian of Qatar's marine environment”] from the Qatar National Research Fund (a member of Qatar Foundation). The findings herein reflect the work, and are solely the responsibility, of the authors. In addition, the authors would like to thank Prof. Eduarda Santos for advice regarding sequencing requirements and Mr. Mark Chatting and Mr. Reyniel Gasang for their support in the collection of the samples.

Viral diseases of crustaceans are increasingly recognised as challenges to shellfish farms and fisheries. Here we describe the first naturally-occurring virus reported in any clawed lobster species. Hypertrophied nuclei with emarginated chromatin, characteristic histopathological lesions of DNA virus infection, were observed within the hepatopancreatic epithelial cells of juvenile European lobsters (Homarus gammarus). Transmission electron microscopy revealed infection with a bacilliform virus containing a rod shaped nucleocapsid enveloped in an elliptical membrane. Assembly of PCR-free shotgun metagenomic sequencing produced a circular genome of 107,063 bp containing 97 open reading frames, the majority of which share sequence similarity with a virus infecting the black tiger shrimp: Penaeus monodon nudivirus (PmNV). Multiple phylogenetic analyses confirm the new virus to be a novel member of the Nudiviridae: Homarus gammarus nudivirus (HgNV). Evidence of occlusion body formation, characteristic of PmNV and its closest relatives, was not observed, questioning the horizontal transmission strategy of HgNV outside of the host. We discuss the potential impacts of HgNV on juvenile lobster growth and mortality and present HgNV-specific primers to serve as a diagnostic tool for monitoring the virus in wild and farmed lobster stocks.

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.

1AbstractMass 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 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.

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.

One of the fastest growing fisheries in the UK is the king scallop ( Pecten maximus L.), also currently rated as the second most valuable fishery. Mass mortality events in scallops have been reported worldwide, often with the causative agent(s) remaining uncharacterized. In May 2013 and 2014, two mass mortality events affecting king scallops were recorded in the Lyme Bay marine protected area (MPA) in Southwest England. Histopathological examination showed gill epithelial tissues infected with intracellular microcolonies (IMCs) of bacteria resembling Rickettsia -like organisms (RLOs), often with bacteria released in vascular spaces. Large colonies were associated with cellular and tissue disruption of the gills. Ultrastructural examination confirmed the intracellular location of these organisms in affected epithelial cells. The 16S rRNA gene sequences of the putative IMCs obtained from infected king scallop gill samples, collected from both mortality events, were identical and had a 99.4% identity to 16S rRNA gene sequences obtained from “ Candidatus Endonucleobacter bathymodioli” and 95% with Endozoicomonas species. In situ hybridization assays using 16S rRNA gene probes confirmed the presence of the sequenced IMC gene in the gill tissues. Additional DNA sequences of the bacterium were obtained using high-throughput (Illumina) sequencing, and bioinformatic analysis identified over 1,000 genes with high similarity to protein sequences from Endozoicomonas spp. (ranging from 77 to 87% identity). Specific PCR assays were developed and applied to screen for the presence of IMC 16S rRNA gene sequences in king scallop gill tissues collected at the Lyme Bay MPA during 2015 and 2016. There was 100% prevalence of the IMCs in these gill tissues, and the 16S rRNA gene sequences identified were identical to the sequence found during the previous mortality event. IMPORTANCE Molluscan mass mortalities associated with IMCs have been reported worldwide for many years; however, apart from histological and ultrastructural characterization, characterization of the etiological agents is limited. In the present work, we provide detailed molecular characterization of an Endozoicomonas -like organism (ELO) associated with an important commercial scallop species.

Vibrio parahaemolyticus is the leading cause of seafood-related infections with illnesses undergoing a geographic expansion. In this process of expansion, the most fundamental change has been the transition from infections caused by local strains to the surge of pandemic clonal types. Pandemic clone sequence type 3 (ST3) was the only example of transcontinental spreading until 2012, when ST36 was detected outside the region where it is endemic in the U.S. Pacific Northwest causing infections along the U.S. northeast coast and Spain. Here, we used genome-wide analyses to reconstruct the evolutionary history of the V. parahaemolyticus ST36 clone over the course of its geographic expansion during the previous 25 years. The origin of this lineage was estimated to be in ~1985. By 1995, a new variant emerged in the region and quickly replaced the old clone, which has not been detected since 2000. The new Pacific Northwest (PNW) lineage was responsible for the first cases associated with this clone outside the Pacific Northwest region. After several introductions into the northeast coast, the new PNW clone differentiated into a highly dynamic group that continues to cause illness on the northeast coast of the United States. Surprisingly, the strains detected in Europe in 2012 diverged from this ancestral group around 2000 and have conserved genetic features present only in the old PNW lineage. Recombination was identified as the major driver of diversification, with some preliminary observations suggesting a trend toward a more specialized lifestyle, which may represent a critical element in the expansion of epidemics under scenarios of coastal warming. IMPORTANCE Vibrio parahaemolyticus and Vibrio cholerae represent the only two instances of pandemic expansions of human pathogens originating in the marine environment. However, while the current pandemic of V. cholerae emerged more than 50 years ago, the global expansion of V. parahaemolyticus is a recent phenomenon. These modern expansions provide an exceptional opportunity to study the evolutionary process of these pathogens at first hand and gain an understanding of the mechanisms shaping the epidemic dynamics of these diseases, in particular, the emergence, dispersal, and successful introduction in new regions facilitating global spreading of infections. In this study, we used genomic analysis to examine the evolutionary divergence that has occurred over the course of the most recent transcontinental expansion of a pathogenic Vibrio , the spreading of the V. parahaemolyticus sequence type 36 clone from the region where it is endemic on the Pacific coast of North America to the east coast of the United States and finally to the west coast of Europe.

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

Background Schizophrenia (SZ) is a severe neuropsychiatric disorder characterized by episodic psychosis and altered cognitive function. About 20% of SZ patients are resistant to the commonly-prescribed antipsychotic medications used to treat disease symptoms. Clozapine is an atypical antipsychotic drug often prescribed to treatment-resistant SZ patients, although the functional pathways mediating its action are not well understood. Importantly, 40–70% of patients treated with clozapine show an inadequate response and there are sever side-effects. Understanding the molecular pathways involved in antipsychotic response will help in the development of new improved therapeutics that act on pathogenicity rather than just treating the acute manifestations of SZ. Methods We exposed wild type zebrafish to high and low doses of the antipsychotic drug clozapine over 72 hours. Relevant controls were included at all stages of the experiment. We recorded the animals and analysed changes in behavior associated with clozapine exposure. We then performed RNA-sequencing in the brain and identified changes in gene expression in response to clozapine exposure. Network analysis was performed to identify modules of co-expressed genes sensitive to clozapine exposure and pathway analysis was employed to identify overrepresentation of GO terms in each of the clozapine exposure-associated modules Results The fish showed a dose-sensitive shift in behaviour in response to clozapine. We identified five genes showing significant gene expression changes in response to clozapine, including ODC1, a gene previously implicated in SZ which shows changes in expression during brain development. Using network analysis, we identified fifteen modules of co-expressed genes which present a striking change in module connectivity in response to the drug. These modules show overrepresentation of key gene ontology terms, implicating processes such as calcium ion transport, cell-cell signalling, G-protein coupled receptor activity, ribosomal activity, and the regulation of transcription. Discussion Our study highlights the utility of zebrafish as a model for assessing the molecular consequences of antipsychotic medications. Our data show marked behavioural effects induced by clozapine and significant transcriptomic alterations in important functional pathways in the brain.