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.

Molecular mechanisms of embryonic tail development in the self-fertilizing mangrove killifish Kryptolebias marmoratus.

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.

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

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.

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.

Molecular Characterization of an Endozoicomonas-Like Organism Causing Infection in the King Scallop (Pecten maximus L.).

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.

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.

'Candidatus Aquirickettsiella gammari' (Gammaproteobacteria: Legionellales: Coxiellaceae): A bacterial pathogen of the freshwater crustacean Gammarus fossarum (Malacostraca: Amphipoda).

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.

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.

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.

Genome Characterization and Phylogenetic Analysis of Scale Drop Disease Virus Isolated from Asian Seabass (Lates calcarifer).

Scale drop disease virus (SDDV), a double-stranded DNA virus in the family Iridoviridae, has been reported widely in southeast Asian countries as a causative agent of scale drop syndrome (SDS) in Asian seabass. SDS has resulted in high mortality and significant economic losses to the aquaculture industry. This study demonstrated the use of metagenomic methods to investigate bacterial and viral communities present in infected fish tissues and recover a complete genome of the causative agent named SDDV TH7_2019. Characterization of the TH7_2019 genome revealed a genome size of 131 kb with 134 putative ORFs encoding viral proteins potentially associated with host apoptosis manipulation. A comparative genome analysis showed a high degree of amino acid identity across SDDV strains, with variations in number of repeat sequences and mutations within core genes. Phylogenetic analyses indicate a close relationship among SDDV genomes. This research enhances our understanding of the genetic diversity and evolutionary relationship of SDDV, contributing valuable insights for further development of effective control strategies of SDDV.

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.

Evolutionary dynamics of the successful expansion of pandemic Vibrio parahaemolyticus ST3 in Latin America.

The underlying evolutionary mechanisms driving global expansions of pathogen strains are poorly understood. Vibrio parahaemolyticus is one of only two marine pathogens where variants have emerged in distinct climates globally. The success of a Vibrio parahaemolyticus clone (VpST3) in Latin America- the first spread identified outside its endemic region of tropical Asia- provided an invaluable opportunity to investigate mechanisms of VpST3 expansion into a distinct marine climate. A global collection of VpST3 isolates and novel Latin American isolates were used for evolutionary population genomics, pangenome analysis and combined with oceanic climate data. We found a VpST3 population (LatAm-VpST3) introduced in Latin America well before the emergence of this clone in India, previously considered the onset of the VpST3 epidemic. LatAm-VpST3 underwent successful adaptation to local conditions over its evolutionary divergence from Asian VpST3 isolates, to become dominant in Latin America. Selection signatures were found in genes providing resilience to the distinct marine climate. Core genome mutations and accessory gene presences that promoted survival over long dispersals or increased environmental fitness were associated with environmental conditions. These results provide novel insights into the global expansion of this successful V. parahaemolyticus clone into regions with different climate scenarios.

Realising a global One Health disease surveillance approach: insights from wastewater and beyond.

One Health is a recognition of the shared environment inhabited by humans, animals and plants, and the impact of their interactions on the health of all organisms. The COVID-19 pandemic highlighted the need for a framework of pathogen surveillance in a tractable One Health paradigm to allow timely detection and response to threats to human and animal health. We present case studies centered around the recent global approach to tackle antimicrobial resistance and the current interest in wastewater testing, with the concept of "one sample many analyses" to be further explored as the most appropriate means of initiating this endeavor.

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.

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.

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.

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.

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.