Reconstructing the Evolutionary History of Pinna nobilis: New Genetic Signals from the Past of a Species on the Brink of Extinction

Pinna nobilis, commonly known as the noble pen shell, is a marine bivalve endemic to the Mediterranean Sea. Unfortunately, due to a multifactorial disease that began affecting its populations in 2016, the species is currently facing the threat of extinction. To gain insights into the evolutionary history of P. nobilis before the mass mortality event (MME), and to obtain a comprehensive understanding of how evolutionary processes led to the adaptation of the species into the Mediterranean Sea, phylogenetic and phylogeographic analyses were carried out. The dataset analysed includes 469 sequences of COI gene fragment both from GenBank and the present study (100). The analysis performed evidenced that P. nobilis diverged about 2.5 mya, after the entrance of its ancestor into the Mediterranean Sea following the Zanclean flood (5.33 mya). Moreover, our results suggest that the starting point of colonisation was the central part of the western Mediterranean basin, with the eastern basin being populated subsequently. From a conservational viewpoint, these results provide important hints for present and future restocking plans, helping to reconstruct the pre-existing genetic variability in sites where the species became extinct.

Plankton and marine aggregates as transmission vectors for V. aestuarianus 02/041 infecting the pacific oyster Crassostrea gigas

Vibrio aestuarianus is a bacterium related to mass mortality outbreaks of the Pacific oyster, Crassostrea gigas in Europe. In this study, the role of different planktonic substrates (phytoplankton cells, marine aggregates and chitin fragments) in mediating V. aestuarianus 02/041 infection of oysters was evaluated by controlled infection experiments. It was shown that phytoplankton cells and, to a greater extent, marine aggregates, significantly promote V. aestuarianus 02/041 intake by C. gigas maintained under stressful conditions in the laboratory. Such intake is associated with higher concentration of the pathogen in the bivalve hemolymph and compromised health status of infected oysters. In contrast, chitin particles do not play a significant role as transmission vector for V. aestuarianus 02/041 infecting its bivalve host. Interestingly, incorporation into marine aggregates foster extracellular proteases (ECPs) activity and a higher expression of bacterial virulence genes, that are potentially involved in bivalve infection. Results from this study contribute to elucidate transmission patterns of V. aestuarianus 02/041 to C. gigas that may be useful for the development of efficient measures to prevent and control oyster disease outbreaks.

Taxonomic Distribution and Molecular Evolution of Mytilectins

R-type lectins are a widespread group of sugar-binding proteins found in nearly all domains of life, characterized by the presence of a carbohydrate-binding domain that adopts a β-trefoil fold. Mytilectins represent a recently described subgroup of β-trefoil lectins, which have been functionally characterized in a few mussel species (Mollusca, Bivalvia) and display attractive properties, which may fuel the development of artificial lectins with different biotechnological applications. The detection of different paralogous genes in mussels, together with the description of orthologous sequences in brachiopods, supports the formal description of mytilectins as a gene family. However, to date, an investigation of the taxonomic distribution of these lectins and their molecular diversification and evolution was still lacking. Here, we provide a comprehensive overview of the evolutionary history of mytilectins, revealing an ancient monophyletic evolutionary origin and a very broad but highly discontinuous taxonomic distribution, ranging from heteroscleromorphan sponges to ophiuroid and crinoid echinoderms. Moreover, the overwhelming majority of mytilectins display a chimera-like architecture, which combines the β-trefoil carbohydrate recognition domain with a C-terminal pore-forming domain, suggesting that the simpler structure of most functionally characterized mytilectins derives from a secondary domain loss.

Viral Diversity in Benthic Abyssal Ecosystems: Ecological and Methodological Considerations

Viruses are the most abundant 'biological entities' in the world's oceans. However, technical and methodological constraints limit our understanding of their diversity, particularly in benthic abyssal ecosystems (>4000 m depth). To verify advantages and limitations of analyzing virome DNA subjected either to random amplification or unamplified, we applied shotgun sequencing-by-synthesis to two sample pairs obtained from benthic abyssal sites located in the North-eastern Atlantic Ocean at ca. 4700 m depth. One amplified DNA sample was also subjected to single-molecule long-read sequencing for comparative purposes. Overall, we identified 24,828 viral Operational Taxonomic Units (vOTUs), belonging to 22 viral families. Viral reads were more abundant in the amplified DNA samples (38.5-49.9%) compared to the unamplified ones (4.4-5.8%), with the latter showing a greater viral diversity and 11-16% of dsDNA viruses almost undetectable in the amplified samples. From a procedural point of view, the viromes obtained by direct sequencing (without amplification step) provided a broader overview of both ss and dsDNA viral diversity. Nevertheless, our results suggest that the contextual use of random amplification of the same sample and long-read technology can improve the assessment of viral assemblages by reducing off-target reads.

Gene expression profiling in white blood cells reveals new insights into the molecular mechanisms of thalidomide in children with inflammatory bowel disease

Thalidomide has emerged as an effective immunomodulator in the treatment of pediatric patients with inflammatory bowel disease (IBD) refractory to standard therapies. Cereblon (CRBN), a component of E3 protein ligase complex that mediates ubiquitination and proteasomal degradation of target proteins, has been identified as the primary target of thalidomide. CRBN plays a crucial role in thalidomide teratogenicity, however it is unclear whether it is also involved in the therapeutic effects in IBD patients. This study aimed at identifying the molecular mechanisms underpinning thalidomide action in pediatric IBD. In this study, ten IBD pediatric patients responsive to thalidomide were prospectively enrolled. RNA-sequencing (RNA-seq) analysis and functional enrichment analysis were carried out on peripheral blood mononuclear cells (PBMC) obtained before and after twelve weeks of treatment with thalidomide. RNA-seq analysis revealed 378 differentially expressed genes before and after treatment with thalidomide. The most deregulated pathways were cytosolic calcium ion concentration, cAMP-mediated signaling, eicosanoid signaling and inhibition of matrix metalloproteinases. Neuronal signaling mechanisms such as CREB signaling in neurons and axonal guidance signaling also emerged. Connectivity Map analysis revealed that thalidomide gene expression changes were similar to those exposed to MLN4924, an inhibitor of NEDD8 activating enzyme, suggesting that thalidomide exerts its immunomodulatory effects by acting on the ubiquitin-proteasome pathway. In vitro experiments on cell lines confirmed the effect of thalidomide on candidate altered pathways observed in patients. These results represent a unique resource for enhanced understanding of thalidomide mechanism in pediatric patients with IBD, providing novel potential targets associated with drug response.

Unveiling the Impact of Gene Presence/Absence Variation in Driving Inter-Individual Sequence Diversity within the CRP-I Gene Family in Mytilus spp

Mussels (Mytilus spp.) tolerate infections much better than other species living in the same marine coastal environment thanks to a highly efficient innate immune system, which exploits a remarkable diversification of effector molecules involved in mucosal and humoral responses. Among these, antimicrobial peptides (AMPs) are subjected to massive gene presence/absence variation (PAV), endowing each individual with a potentially unique repertoire of defense molecules. The unavailability of a chromosome-scale assembly has so far prevented a comprehensive evaluation of the genomic arrangement of AMP-encoding loci, preventing an accurate ascertainment of the orthology/paralogy relationships among sequence variants. Here, we characterized the CRP-I gene cluster in the blue mussel Mytilus edulis, which includes about 50 paralogous genes and pseudogenes, mostly packed in a small genomic region within chromosome 5. We further reported the occurrence of widespread PAV within this family in the Mytilus species complex and provided evidence that CRP-I peptides likely adopt a knottin fold. We functionally characterized the synthetic peptide sCRP-I H1, assessing the presence of biological activities consistent with other knottins, revealing that mussel CRP-I peptides are unlikely to act as antimicrobial agents or protease inhibitors, even though they may be used as defense molecules against infections from eukaryotic parasites.

Expansion and neofunctionalization of actinoporin-like genes in Mediterranean mussel (Mytilus galloprovincialis)

Pore-forming toxins are an important component of the venom of many animals. Actinoporins are potent cytolysins that were first detected in the venom of sea anemones; however, they are occasionally found in animals other than cnidarians and are expanded in a few predatory gastropods. Here, we report the presence of 27 unique actinoporin-like genes with monophyletic origin in Mytilus galloprovincialis, which we have termed mytiporins. These mytiporins exhibited a remarkable level of molecular diversity and gene presence–absence variation, which warranted further studies aimed at elucidating their functional role. We structurally and functionally characterized mytiporin-1 and found significant differences from the archetypal actinoporin fragaceatoxin C. Mytiporin-1 showed weaker permeabilization activity, no specificity towards sphingomyelin, and weak activity in model lipid systems with negatively charged lipids. In contrast to fragaceatoxin C, which forms octameric pores, functional mytiporin-1 pores on negatively charged lipid membranes were hexameric. Similar hexameric pores were observed for coluporin-26 from Cumia reticulata and a conoporin from Conus andremenezi. This indicates that also other molluscan actinoporin-like proteins differ from fragaceatoxin C. Although the functional role of mytiporins in the context of molluscan physiology remains to be elucidated, the lineage-specific gene family expansion event that characterizes mytiporins indicates that strong selective forces acted on their molecular diversification. Given the tissue distribution of mytiporins, this process may have broadened the taxonomic breadth of their biological targets, which would have important implications for digestive processes or mucosal immunity.

Amplicon sequence variant-based meiofaunal community composition revealed by DADA2 tool is compatible with species composition

The present study is aimed at implementing the morphological identification-free amplicon sequence variant (ASV) concept for describing meiofaunal species composition, while strongly indicating reasonable compatibility with the underlying species. A primer pair was constructed and demonstrated to PCR amplify a 470-490 bp 18S barcode from a variety of meiofaunal taxa, high throughput sequenced using the Illumina 300 × 2 bps platform. Sixteen 18S multi-species HTS assemblies were created from meiofaunal samples and merged to one assembly of ~2,150,000 reads. Five quality scores (q = 35, 30, 25, 20, 15) were implemented to filter five 18S barcode assemblies, which served as inputs for the DADA2 software, ending with five reference ASV libraries. Each of these libraries was clustered, applying 3% dissimilarity threshold, revealed an average number of 1.38 ± 0.078 ASVs / cluster. Hence, demonstrating high level of ASV uniqueness. The libraries which were based on q ≤ 25 reached a near-asymptote number of ASVs which together with the low average number of ASVs / cluster, strongly indicated fair representation of the actual number of the underlying species. Hence, the q = 25 library was selected to be used as metabarcoding reference library. It contained 461 ASVs and 342-3% clusters with average number of 1.34 ± 1.036 ASV / cluster and their BLASTN annotation elucidated a variety of expected meiofaunal taxa. The sixteen assemblies of sample-specific paired reads were mapped to this reference library and sample ASV profiles, namely the list of ASVs and their proportional copy numbers were created and clustered.

The role of outdoor and indoor air quality in the spread of SARS-CoV-2: Overview and recommendations by the research group on COVID-19 and particulate matter (RESCOP commission)

There are important questions surrounding the potential contribution of outdoor and indoor air quality in the transmission of SARS-CoV-2 and perpetuation of COVID-19 epidemic waves. Environmental health may be a critical component of COVID-19 prevention. The public health community and health agencies should consider the evolving evidence in their recommendations and statements, and work to issue occupational guidelines. Evidence coming from the current epidemiological and experimental research is expected to add knowledge about virus diffusion, COVID-19 severity in most polluted areas, inter-personal distance requirements and need for wearing face masks in indoor or outdoor environments. The COVID-19 pandemic has highlighted the need for maintaining particulate matter concentrations at low levels for multiple health-related reasons, which may also include the spread of SARS-CoV-2. Indoor environments represent even a more crucial challenge to cope with, as it is easier for the SARS-COV2 to spread, remain vital and infect other subjects in closed spaces in the presence of already infected asymptomatic or mildly symptomatic people. The potential merits of preventive measures, such as CO₂ monitoring associated with natural or controlled mechanical ventilation and air purification, for schools, indoor public places (restaurants, offices, hotels, museums, theatres/cinemas etc.) and transportations need to be carefully considered. Hospital settings and nursing/retirement homes as well as emergency rooms, infectious diseases divisions and ambulances represent higher risk indoor environments and may require additional monitoring and specific decontamination strategies based on mechanical ventilation or air purification.

Anisaxins, helical antimicrobial peptides from marine parasites, kill resistant bacteria by lipid extraction and membrane disruption

An infecting and propagating parasite relies on its innate defense system to evade the host's immune response and to survive challenges from commensal bacteria. More so for the nematode Anisakis, a marine parasite that during its life cycle encounters both vertebrate and invertebrate hosts and their highly diverse microbiotas. Although much is still unknown about how the nematode mitigates the effects of these microbiota, its antimicrobial peptides likely play an important role in its survival. We identified anisaxins, the first cecropin-like helical antimicrobial peptides originating from a marine parasite, by mining available genomic and transcriptomic data for Anisakis spp. These peptides are potent bactericidal agents in vitro, selectively active against Gram-negative bacteria, including multi-drug resistant strains, at sub-micromolar concentrations. Their interaction with bacterial membranes was confirmed by solid state NMR (ssNMR) and is highly dependent on the peptide concentration as well as peptide to lipid ratio, as evidenced by molecular dynamics (MD) simulations. MD results indicated that an initial step in the membranolytic mode of action involves membrane bulging and lipid extraction; a novel mechanism which may underline the peptides' potency. Subsequent steps include membrane permeabilization leading to leakage of molecules and eventually cell death, but without visible macroscopic damage, as shown by atomic force microscopy and flow cytometry. This membranolytic antibacterial activity does not translate to cytotoxicity towards human peripheral blood mononuclear cells (HPBMCs), which was minimal at well above bactericidal concentrations, making anisaxins promising candidates for further drug development. STATEMENT OF SIGNIFICANCE: Witnessing the rapid spread of antibiotic resistance resulting in millions of infected and dozens of thousands dying worldwide every year, we identified anisaxins, antimicrobial peptides (AMPs) from marine parasites, Anisakis spp., with potent bactericidal activity and selectivity towards multi-drug resistant Gram-negative bacteria. Anisaxins are membrane-active peptides, whose activity, very sensitive to local peptide concentrations, involves membrane bulging and lipid extraction, leading to membrane permeabilization and bacterial cell death. At the same time, their toxicity towards host cells is negligible, which is often not the case for membrane-active AMPs, therefore making them suitable drug candidates. Membrane bulging and lipid extraction are novel concepts that broaden our understanding of peptide interactions with bacterial functional structures, essential for future design of such biomaterials.

Comparative analysis of novel and common reference genes in adult tissues of the mussel Mytilus galloprovincialis

Background

Real-time quantitative PCR is a widely used method for gene expression analyses in various organisms. Its accuracy mainly relies on the correct selection of reference genes. Any experimental plan involving real-time PCR needs to evaluate the characteristics of the samples to be examined and the relative stability of reference genes. Most studies in mollusks rely on reference genes commonly used in vertebrates.

Results

In this study, we focused on the transcriptome of the bivalve mollusk Mytilus galloprovincialis in physiological state to identify suitable reference genes in several adult tissues. Candidate genes with highly stable expression across 51 RNA-seq datasets from multiple tissues were selected through genome-wide bioinformatics analysis. This approach led to the identification of three genes (Rpl14, Rpl32 and Rpl34), whose suitability was evaluated together with 7 other reference genes commonly reported in literature (Act, Cyp-A, Ef1α, Gapdh, 18S, 28S and Rps4). The stability analyses performed with geNorm, NormFinder and Bestkeeper identified specific either single or pairs of genes suitable as references for gene expression analyses in specific tissues and revealed the Act/Cyp-A pair as the most appropriate to analyze gene expression across different tissues.

Conclusion

Mytilus galloprovincialis is a model system increasingly used in ecotoxicology and molecular studies. Our transcriptome-wide approach represents the first comprehensive investigation aimed at the identification of suitable reference genes for expression studies in this species.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08553-1.

Spectroscopic investigation of faeces with surface-enhanced Raman scattering: a case study with coeliac patients on gluten-free diet

Surface-enhanced Raman scattering (SERS) spectra of faecal samples can be obtained by adding AuNP to their methanol extracts according to the reported protocol, and display bands that are due to bilirubin-like species but also to xanthine and hypoxanthine, two metabolic products secreted by gut bacteria. A total of 27 faecal samples from three different groups, i.e. coeliac patients (n = 9), coeliac patients on gluten-free diet (n = 10) and a control group (n = 8), were characterized with both SERS spectroscopy and 16S rRNA sequencing analysis. Significant differences are present between SERS spectra of coeliac patients and those on gluten-free diet, with a marked increase in the relative intensity of both xanthine and hypoxanthine for the latter. Interestingly, these differences do not correlate with bacterial composition as derived from 16S rRNA sequencing.

Recent advances in bivalve-microbiota interactions for disease prevention in aquaculture

In bivalves, no clear-cut functional role of microbiota has yet been identified, although many publications suggest that they could be involved in nutrition or immunity of their host. In the context of climate change, integrative approaches at the crossroads of disciplines have been developed to explore the environment-host-pathogen-microbiota system. Here, we attempt to synthesize work on (1) the current methodologies to analyse bivalve microbiota, (2) the comparison of microbiota between species, between host compartments and their surrounding habitat, (3) how the bivalve microbiota are governed by environmental factors and host genetics and (4) how host-associated microorganisms act as a buffer against pathogens and/or promote recovery, and could thereby play a role in the prevention of disease or mortalities.

Environmental DNA reveals diversity and abundance of Alternaria species in neighbouring heterogeneous landscapes in Worcester, UK

Alternaria is a pathogenic and allergenic fungus affecting 400 plant species and 334 million people globally. This study aimed at assessing the diversity of Alternaria species in airborne samples collected from closely located (7 km apart) and heterogeneous sites (rural, urban and unmanaged grassland) in Worcester and Lakeside, the UK. A secondary objective was to examine how the ITS1 subregion varies from ITS2 in Alternaria species diversity and composition. Airborne spores were collected using Burkard 7-day and multi-vial Cyclone samplers for the period 5 July 2016-9 October 2019. Air samples from the Cyclone were amplified using the ITS1and ITS2 subregions and sequenced using Illumina MiSeq platform whereas those from the Burkard sampler were identified and quantified using optical microscopy. Optical microscopy and eDNA revealed a high abundance of Alternaria in the rural, urban and unmanaged sites. ITS1 and ITS2 detected five and seven different Alternaria species at the three sampling sites, respectively. A. dactylidicola, A. metachromatica and A. infectoria were the most abundant. The rural, urban and unmanaged grassland sites had similar diversity (PERMANOVA) of the species due to similarity in land use and proximity of the sites. Overall, the study showed that heterogeneous and neighbouring sites with similar land uses can have similar Alternaria species. It also demonstrated that an eDNA approach can complement the classical optical microscopy method in providing more precise information on fungal species diversity in an environment for targeted management. Similar studies can be replicated for other allergenic and pathogenic fungi.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s10453-022-09760-9.

The mitochondrial protein Opa1 promotes adipocyte browning that is dependent on urea cycle metabolites

White to brown/beige adipocytes conversion is a possible therapeutic strategy to tackle the current obesity epidemics. While mitochondria are key for energy dissipation in brown fat, it is unknown if they can drive adipocyte browning. Here, we show that the mitochondrial cristae biogenesis protein optic atrophy 1 (Opa1) facilitates cell-autonomous adipocyte browning. In two cohorts of patients with obesity, including weight discordant monozygotic twin pairs, adipose tissue OPA1 levels are reduced. In the mouse, Opa1 overexpression favours white adipose tissue expandability as well as browning, ultimately improving glucose tolerance and insulin sensitivity. Transcriptomics and metabolomics analyses identify the Jumanji family chromatin remodelling protein Kdm3a and urea cycle metabolites, including fumarate, as effectors of Opa1-dependent browning. Mechanistically, the higher cyclic adenosine monophosphate (cAMP) levels in Opa1 pre-adipocytes activate cAMP-responsive element binding protein (CREB), which transcribes urea cycle enzymes. Flux analyses in pre-adipocytes indicate that Opa1-dependent fumarate accumulation depends on the urea cycle. Conversely, adipocyte-specific Opa1 deletion curtails urea cycle and beige differentiation of pre-adipocytes, and is rescued by fumarate supplementation. Thus, the urea cycle links the mitochondrial dynamics protein Opa1 to white adipocyte browning.

High-Throughput Sequencing of Environmental DNA as a Tool for Monitoring Eukaryotic Communities and Potential Pathogens in a Coastal Upwelling Ecosystem

The marine environment includes diverse microeukaryotic organisms that play important functional roles in the ecosystem. With molecular approaches, eukaryotic taxonomy has been improved, complementing classical analysis. In this study, DNA metabarcoding was performed to describe putative pathogenic eukaryotic microorganisms in sediment and marine water fractions collected in Galicia (NW Spain) from 2016 to 2018. The composition of eukaryotic communities was distinct between sediment and water fractions. Protists were the most diverse group, with the clade TSAR (Stramenopiles, Alveolata, Rhizaria, and Telonemida) as the primary representative organisms in the environment. Harmful algae and invasive species were frequently detected. Potential pathogens, invasive pathogenic organisms as well as the causative agents of harmful phytoplanktonic blooms were identified in this marine ecosystem. Most of the identified pathogens have a crucial impact on the aquacultural sector or affect to relevant species in the marine ecosystem, such as diatoms. Moreover, pathogens with medical and veterinary importance worldwide were also found, as well as pathogens that affect diatoms. The evaluation of the health of a marine ecosystem that directly affects the aquacultural sector with a zoonotic concern was performed with the metabarcoding assay.

Improved high throughput protocol for targeting eukaryotic symbionts in metazoan and eDNA samples

Eukaryote symbionts of animals are major drivers of ecosystems not only because of their diversity and host interactions from variable pathogenicity but also through different key roles such as commensalism and to different types of interdependence. However, molecular investigations of metazoan eukaryomes require minimising coamplification of homologous host genes. In this study we (1) identified a previously published “antimetazoan” reverse primer to theoretically enable amplification of a wider range of microeukaryotic symbionts, including more evolutionarily divergent sequence types, (2) evaluated in silico several antimetazoan primer combinations, and (3) optimised the application of the best performing primer pair for high throughput sequencing (HTS) by comparing one‐step and two‐step PCR amplification approaches, testing different annealing temperatures and evaluating the taxonomic profiles produced by HTS and data analysis. The primer combination 574*F – UNonMet_DB tested in silico showed the largest diversity of nonmetazoan sequence types in the SILVA database and was also the shortest available primer combination for broadly‐targeting antimetazoan amplification across the 18S rRNA gene V4 region. We demonstrate that the one‐step PCR approach used for library preparation produces significantly lower proportions of metazoan reads, and a more comprehensive coverage of host‐associated microeukaryote reads than the two‐step approach. Using higher PCR annealing temperatures further increased the proportion of nonmetazoan reads in all sample types tested. The resulting V4 region amplicons were taxonomically informative even when only the forward read is analysed. This region also revealed a diversity of known and putatively parasitic lineages and a wider diversity of host‐associated eukaryotes.

The proline-rich myticalins from Mytilus galloprovincialis display a membrane-permeabilizing antimicrobial mode of action

Bivalve mollusks are continuously exposed to potentially pathogenic microorganisms living in the marine environment. Not surprisingly, these filter-feeders developed a robust innate immunity to protect themselves, which includes a broad panel of antimicrobial peptides. Among these, myticalins represent a recently discovered family of linear cationic peptides expressed in the gills of Mytilus galloprovincialis. Even though myticalins and insect and mammalian proline-rich antimicrobial peptides (PrAMPs) share a similar amino acid composition, we here show that none of the tested mussel peptides use a non-lytic mode of action relying on the bacterial transporter SbmA. On the other hand, all the tested myticalins perturbed and permeabilized the membranes of E. coli BW25113, as shown by flow-cytometry and atomic force microscopy. Circular dichroism spectra revealed that most myticalins did not adopt recognizable secondary structures in the presence of amphipathic environments, such as biological membranes. To explore possible uses of myticalins for biotech, we assessed their biocompatibility with a human cell line. Non-negligible cytotoxic effects displayed by myticalins indicate that their optimization would be required before their further use as lead compounds in the development of new antibiotics.

Suitability of a dual COI marker for marine zooplankton DNA metabarcoding

As DNA metabarcoding has become an emerging tool for surveying biodiversity, including its application in legally binding assessments, reliable and efficient barcodes are requested, especially for the highly diverse group of zooplankton. This study focuses on comparing the efficiency of two mitochondrial COI barcodes based on the internal primers mlCOIintF and mlCOIintR utilizing mesozooplankton samples collected in a Mediterranean lagoon. Our results indicate that after a slight adjustment, the mlCOIintR primer performs in combination with jdgLCO1490 (herein) very comparably to the much more widely used primer system mlCOIintF/jgHCO2198+dgHCO2198, in terms of level of taxonomic resolution, species detection and their relative abundance in terms of numbers of reads. As for some groups, like Ctenophora, this barcode is not suitable; a combination of them may be the best option to rely on the Folmer region in its entirety without the risk of losing information for a limited primer match.

Evolutionary History of DNA Methylation Related Genes in Bivalvia: New Insights From Mytilus galloprovincialis

DNA methylation is an essential epigenetic mechanism influencing gene expression in all organisms. In metazoans, the pattern of DNA methylation changes during embryogenesis and adult life. Consequently, differentiated cells develop a stable and unique DNA methylation pattern that finely regulates mRNA transcription during development and determines tissue-specific gene expression. Currently, DNA methylation remains poorly investigated in mollusks and completely unexplored in Mytilus galloprovincialis. To shed light on this process in this ecologically and economically important bivalve, we screened its genome, detecting sequences homologous to DNA methyltransferases (DNMTs), methyl-CpG-binding domain (MBD) proteins and Ten-eleven translocation methylcytosine dioxygenase (TET) previously described in other organisms. We characterized the gene architecture and protein domains of the mussel sequences and studied their phylogenetic relationships with the ortholog sequences from other bivalve species. We then comparatively investigated their expression levels across different adult tissues in mussel and other bivalves, using previously published transcriptome datasets. This study provides the first insights on DNA methylation regulators in M. galloprovincialis, which may provide fundamental information to better understand the complex role played by this mechanism in regulating genome activity in bivalves.

Molecular detection of SARS-CoV-2 from indoor air samples in environmental monitoring needs adequate temporal coverage and infectivity assessment

The relevance of airborne exposure to SARS-CoV-2 in indoor environments is a matter of research and debate, with special importance for healthcare low-risk settings. Experimental approaches to the bioaerosol sampling are neither standardized nor optimized yet, leading in some cases to limited representativity of the temporal and spatial variability of viral presence in aerosols. Airborne viral viability moreover needs to be assessed. A study has been conducted collecting five 24-h PM10 samples in a COVID-19 geriatric ward in late June 2020, and detecting E and RdRp genes by RT-qPCR with a Ct between 36 and 39. The viral RNA detection at Ct = 36 was related to the maximal numerosity of infected patients hosted in the ward. Lacking a direct infectivity assessment for the collected samples an experimental model has been defined, by seeding twelve nasopharyngeal swab extracts from COVID-19 positive patients on Vero E6 cells; only the four extracts with a viral load above E+10 viral copies (approximately Ct<24) have been able to establish a persistent infection in vitro. Therefore, the cytopathic effect, a key feature of residual infectivity, could be considered unlikely for the environmental PM₁₀ samples showing amplification of viral RNA at Ct = 36 or higher. A standardization of airborne SARS-CoV-2 long-term monitoring and of environmental infectivity assessment is urgently needed.

Innate immune-gene expression during experimental amyloodiniosis in European seabass (Dicentrarchus labrax)

The ectoparasite protozoan Amyloodinium ocellatum (AO) is the causative agent of amyloodiniosis in European seabass (ESB, Dicentrarchus labrax). There is a lack of information about basic molecular immune response mechanisms of ESB during AO infestation. Therefore, to compare gene expression between experimental AO-infested ESB tissues and uninfested ESB tissues (gills and head kidney) RNA-seq was adopted. The RNA-seq revealed multiple differentially expressed genes (DEG), namely 679 upregulated genes and 360 downregulated genes in the gills, and 206 upregulated genes and 170 downregulated genes in head kidney. In gills, genes related to the immune system (perforin, CC1) and protein binding were upregulated. Several genes involved in IFN related pathways were upregulated in the head kidney. Subsequently, to validate the DEG from amyloodiniosis, 26 ESB (mean weight 14 g) per tank in triplicate were bath challenged for 2 h with AO (3.5 × 10⁶/tank; 70 dinospores/mL) under controlled conditions (26-28 °C and 34‰ salinity). As a control group (non-infested), 26 ESB per tank in triplicate were also used. Changes in the expression of innate immune genes in gills and head kidney at 2, 3, 5, 7 and 23 dpi were analysed using real-time PCR. The results indicated that the expression of cytokines (CC1, IL-8) and antimicrobial peptide (Hep) were strongly stimulated and reached a peak at 5 dpi in the early infestation stage, followed by a gradual reduction in the recovery stage (23 dpi). Noticeably, the immunoglobulin (IgM) expression was higher at 23 dpi compared to 7 dpi. Furthermore, in-situ hybridization showed positive signals of CC1 mRNA in AO infested gills compared to the control group. Altogether, chemokines were involved in the immune process under AO infestation and this evidence allows a better understanding of the immune response in European seabass during amyloodiniosis.

Cold Adaptation in Antarctic Notothenioids: Comparative Transcriptomics Reveals Novel Insights in the Peculiar Role of Gills and Highlights Signatures of Cobalamin Deficiency

Far from being devoid of life, Antarctic waters are home to Cryonotothenioidea, which represent one of the fascinating cases of evolutionary adaptation to extreme environmental conditions in vertebrates. Thanks to a series of unique morphological and physiological peculiarities, which include the paradigmatic case of loss of hemoglobin in the family Channichthyidae, these fish survive and thrive at sub-zero temperatures. While some of the distinctive features of such adaptations have been known for decades, our knowledge of their genetic and molecular bases is still limited. We generated a reference de novo assembly of the icefish Chionodraco hamatus transcriptome and used this resource for a large-scale comparative analysis among five red-blooded Cryonotothenioidea, the sub-Antarctic notothenioid Eleginops maclovinus and seven temperate teleost species. Our investigation targeted the gills, a tissue of primary importance for gaseous exchange, osmoregulation, ammonia excretion, and its role in fish immunity. One hundred and twenty genes were identified as significantly up-regulated in Antarctic species and surprisingly shared by red- and white-blooded notothenioids, unveiling several previously unreported molecular players that might have contributed to the evolutionary success of Cryonotothenioidea in Antarctica. In particular, we detected cobalamin deficiency signatures and discussed the possible biological implications of this condition concerning hematological alterations and the heavy parasitic loads typically observed in all Cryonotothenioidea.

Gene expression response of the alga Fucus virsoides (Fucales, Ochrophyta) to glyphosate solution exposure

Fucus virsoides is an ecologically important canopy-forming brown algae endemic to the Adriatic Sea. Once widespread in marine coastal areas, this species underwent a rapid population decline and is now confined to small residual areas. Although the reasons behind this progressive disappearance are still a matter of debate, F. virsoides may suffer, like other macroalgae, from the potential toxic effects of glyphosate-based herbicides. Here, through a transcriptomic approach, we investigate the molecular basis of the high susceptibility of this species to glyphosate solution, previously observed at the morphological and eco-physiological levels. By simulating runoff event in a factorial experiment, we exposed F. virsoides to glyphosate (Roundup® 2.0), either alone or in association with nutrient enrichment, highlighting significant alterations of gene expression profiles that were already visible after three days of exposure. In particular, glyphosate exposure determined the near-complete expression shutdown of several genes involved in photosynthesis, protein synthesis and stress response molecular pathways. Curiously, these detrimental effects were partially mitigated by nutrient supplementation, which may explain the survival of relict population in confined areas with high nutrient inputs.

PLANiTS: a curated sequence reference dataset for plant ITS DNA metabarcoding

DNA metabarcoding combines DNA barcoding with high-throughput sequencing to identify different taxa within environmental communities. The ITS has already been proposed and widely used as universal barcode marker for plants, but a comprehensive, updated and accurate reference dataset of plant ITS sequences has not been available so far. Here, we constructed reference datasets of Viridiplantae ITS1, ITS2 and entire ITS sequences including both Chlorophyta and Streptophyta. The sequences were retrieved from NCBI, and the ITS region was extracted. The sequences underwent identity check to remove misidentified records and were clustered at 99% identity to reduce redundancy and computational effort. For this step, we developed a script called 'better clustering for QIIME' (bc4q) to ensure that the representative sequences are chosen according to the composition of the cluster at a different taxonomic level. The three datasets obtained with the bc4q script are PLANiTS1 (100 224 sequences), PLANiTS2 (96 771 sequences) and PLANiTS (97 550 sequences), and all are pre-formatted for QIIME, being this the most used bioinformatic pipeline for metabarcoding analysis. Being curated and updated reference databases, PLANiTS1, PLANiTS2 and PLANiTS are proposed as a reliable, pivotal first step for a general standardization of plant DNA metabarcoding studies. The bc4q script is presented as a new tool useful in each research dealing with sequences clustering. Database URL: https://github.com/apallavicini/bc4q; https://github.com/apallavicini/PLANiTS.

Massive gene presence-absence variation shapes an open pan-genome in the Mediterranean mussel

BACKGROUND

The Mediterranean mussel Mytilus galloprovincialis is an ecologically and economically relevant edible marine bivalve, highly invasive and resilient to biotic and abiotic stressors causing recurrent massive mortalities in other bivalves. Although these traits have been recently linked with the maintenance of a high genetic variation within natural populations, the factors underlying the evolutionary success of this species remain unclear.

RESULTS

Here, after the assembly of a 1.28-Gb reference genome and the resequencing of 14 individuals from two independent populations, we reveal a complex pan-genomic architecture in M. galloprovincialis, with a core set of 45,000 genes plus a strikingly high number of dispensable genes (20,000) subject to presence-absence variation, which may be entirely missing in several individuals. We show that dispensable genes are associated with hemizygous genomic regions affected by structural variants, which overall account for nearly 580 Mb of DNA sequence not included in the reference genome assembly. As such, this is the first study to report the widespread occurrence of gene presence-absence variation at a whole-genome scale in the animal kingdom.

CONCLUSIONS

Dispensable genes usually belong to young and recently expanded gene families enriched in survival functions, which might be the key to explain the resilience and invasiveness of this species. This unique pan-genome architecture is characterized by dispensable genes in accessory genomic regions that exceed by orders of magnitude those observed in other metazoans, including humans, and closely mirror the open pan-genomes found in prokaryotes and in a few non-metazoan eukaryotes.

Insight into the microbial communities associated with first larval stages of Mytilus galloprovincialis: Possible interference by estrogenic compounds

The microbiota, the host-associated community of microbes, play important roles in health status and whole body homeostasis of all organisms, including marine species. In bivalves, the microbiota composition has been mainly investigated in adults, whereas little information is available during development. In this work, the microbiota composition of the first larval stages of Mytilus galloprovincialis was evaluated by 16S rRNA gene-based profiling, at 24 and 48 hours post fertilization in comparison with those of eggs and sperm. The main genera detected in both larvae (Vibrio, Pseudoalteromonas, Psychrobium, Colwellia) derived from eggs. However, a clear shift in microbiota was observed in developing larvae compared to eggs, both in terms of core microbiome and relative abundance of different genera. The results provide a first insight into the composition of the microbial communities associated with gametes and early larvae of mussels. Moreover, the impact on larval microbiome of estrogenic chemicals that potentially affect Mytilus early development, 17βestradiol-E₂, Bisphenol A-BPA and Bisphenol F-BPF (10 μg/L), was investigated. Exposure to estrogenic chemicals leads to changes in abundance of different genera, with distinct and common effects depending on the compound and larval stage. Both potential pathogens (Vibrio, Arcobacter, Tenacibaculum) and genera involved in xenobiotic biotransformation (Oleispira, Shewanella) were affected. The effects of estrogenic compounds on larval microbiome were not related to their developmental effects: however, the results address the importance of evaluating the impact of emerging contaminants on the microbiota of marine invertebrates, including larval stages, that are most sensitive to environmental perturbations.

Environmental DNA assessment of airborne plant and fungal seasonal diversity

Environmental DNA (eDNA) metabarcoding and metagenomics analyses can improve taxonomic resolution in biodiversity studies. Only recently, these techniques have been applied in aerobiology, to target bacteria, fungi and plants in airborne samples. Here, we present a nine-month aerobiological study applying eDNA metabarcoding in which we analyzed simultaneously airborne diversity and variation of fungi and plants across five locations in North and Central Italy. We correlated species composition with the ecological characteristics of the sites and the seasons. The most abundant taxa among all sites and seasons were the fungal genera Cladosporium, Alternaria, and Epicoccum and the plant genera Brassica, Corylus, Cupressus and Linum, the latter being much more variable among sites. PERMANOVA and indicator species analyses showed that the plant diversity from air samples is significantly correlated with seasons, while that of fungi varied according to the interaction between seasons and sites. The results consolidate the performance of a new eDNA metabarcoding pipeline for the simultaneous amplification and analysis of airborne plant and fungal particles. They also highlight the promising complementarity of this approach with more traditional biomonitoring frameworks and routine reports of air quality provided by environmental agencies.

Identification of an IgD/IgT chimera in the European sea bass (Dicentrarchus labrax L.)

Three classes of immunoglobulins have been identified in Teleosts: IgM, IgT/Z and IgD. They are fundamental for fish immune responses and, therefore, their functional activities are heavily investigated. In this paper, we describe the identification of a new IgD/IgT chimera in sea bass (Dicentrarchus labrax) from a gills transcriptome. This transcript joined the first six constant domains of the IgD chain with the two terminal constant domains of IgT, generating a long in-frame coding sequence with a junction between the canonical δ6 exon splicing donor site and the τ3 exon splicing acceptor site. Studies performed on genomic DNA confirmed the presence of the sequence and identifies and intronic region of 656 bp within this joining region. The basal expression of the IgD/IgT chimera was investigated both in silico and in vivo: high level of expression was found in gills, gut and head kidney. Moreover, IgD/IgT transcripts were up-regulated after in vitro stimulation of sea bass HK leukocytes with LPS. The IgD/IgT chimera was found also in two congener species, Morone saxatilis and Morone chrysops. It is not possible to have a precise idea on the evolutionary scenario that lead to the appearance of this sequence due to the lack of genomic information, but we could speculate that an ancestral duplication of the entire IgH locus was followed by the chimerization of Cδ/Cτ in one of the two loci. Finally, the IgD/IgT high basal expression in tissues and organs fundamental for sea bass immune response and its modulation after LPS stimulation provide a very preliminary indication that this unusual Ig variant could have a functional activity.

SARS-Cov-2RNA found on particulate matter of Bergamo in Northern Italy: First evidence

BACKGROUND

The burden of COVID-19 was extremely severe in Northern Italy, an area characterized by high concentrations of particulate matter (PM), which is known to negatively affect human health. Consistently with evidence already available for other viruses, we initially hypothesized the possibility of SARS-CoV-2 presence on PM, and we performed a first experiment specifically aimed at confirming or excluding this research hyphotesys.

METHODS

We have collected 34 PM10 samples in Bergamo area (the epicenter of the Italian COVID-19 epidemic) by using two air samplers over a continuous 3-weeks period. Filters were properly stored and underwent RNA extraction and amplification according to WHO protocols in two parallel blind analyses performed by two different authorized laboratories. Up to three highly specific molecular marker genes (E, N, and RdRP) were used to test the presence of SARS-CoV-2 RNA on particulate matter.

RESULTS

The first test showed positive results for gene E in 15 out of 16 samples, simultaneously displaying positivity also for RdRP gene in 4 samples. The second blind test got 5 additional positive results for at least one of the three marker genes. Overall, we tested 34 RNA extractions for the E, N and RdRP genes, reporting 20 positive results for at least one of the three marker genes, with positivity separately confirmed for all the three markers. Control tests to exclude false positivities were successfully accomplished.

CONCLUSION

This is the first evidence that SARS-CoV-2 RNA can be present on PM, thus suggesting a possible use as indicator of epidemic recurrence.

SARS-Cov-2RNA found on particulate matter of Bergamo in Northern Italy: First evidence

BACKGROUND

The burden of COVID-19 was extremely severe in Northern Italy, an area characterized by high concentrations of particulate matter (PM), which is known to negatively affect human health. Consistently with evidence already available for other viruses, we initially hypothesized the possibility of SARS-CoV-2 presence on PM, and we performed a first experiment specifically aimed at confirming or excluding this research hyphotesys.

METHODS

We have collected 34 PM10 samples in Bergamo area (the epicenter of the Italian COVID-19 epidemic) by using two air samplers over a continuous 3-weeks period. Filters were properly stored and underwent RNA extraction and amplification according to WHO protocols in two parallel blind analyses performed by two different authorized laboratories. Up to three highly specific molecular marker genes (E, N, and RdRP) were used to test the presence of SARS-CoV-2 RNA on particulate matter.

RESULTS

The first test showed positive results for gene E in 15 out of 16 samples, simultaneously displaying positivity also for RdRP gene in 4 samples. The second blind test got 5 additional positive results for at least one of the three marker genes. Overall, we tested 34 RNA extractions for the E, N and RdRP genes, reporting 20 positive results for at least one of the three marker genes, with positivity separately confirmed for all the three markers. Control tests to exclude false positivities were successfully accomplished.

CONCLUSION

This is the first evidence that SARS-CoV-2 RNA can be present on PM, thus suggesting a possible use as indicator of epidemic recurrence.

DNA metabarcoding and morphological analysis - Assessment of zooplankton biodiversity in transitional waters

Zooplankton biodiversity assessment is a crucial element in monitoring marine ecosystem processes and community responses to environmental alterations. In order to evaluate the suitability of metabarcoding for zooplankton biodiversity assessment and biomonitoring as a fast and more cost-effective method, seasonal zooplankton sampling was carried out in the Venice Lagoon and the nearby coastal area (Northern Adriatic Sea). The molecular analysis showed higher taxa richness compared to the classical morphological method (224 vs. 88 taxa), discriminating better the meroplanktonic component, morphologically identified only up to order level. Both methods revealed a similar spatio-temporal distribution pattern and the sequence abundances and individual counts were significantly correlated for various taxonomic groups. These results indicate that DNA metabarcoding is an efficient tool for biodiversity assessments in ecosystems with high spatial and temporal variability, where high sampling effort is required as well as fast alert systems for non-native species (NIS).

SARS-Cov-2RNA found on particulate matter of Bergamo in Northern Italy: First evidence

BACKGROUND

The burden of COVID-19 was extremely severe in Northern Italy, an area characterized by high concentrations of particulate matter (PM), which is known to negatively affect human health. Consistently with evidence already available for other viruses, we initially hypothesized the possibility of SARS-CoV-2 presence on PM, and we performed a first experiment specifically aimed at confirming or excluding this research hyphotesys.

METHODS

We have collected 34 PM10 samples in Bergamo area (the epicenter of the Italian COVID-19 epidemic) by using two air samplers over a continuous 3-weeks period. Filters were properly stored and underwent RNA extraction and amplification according to WHO protocols in two parallel blind analyses performed by two different authorized laboratories. Up to three highly specific molecular marker genes (E, N, and RdRP) were used to test the presence of SARS-CoV-2 RNA on particulate matter.

RESULTS

The first test showed positive results for gene E in 15 out of 16 samples, simultaneously displaying positivity also for RdRP gene in 4 samples. The second blind test got 5 additional positive results for at least one of the three marker genes. Overall, we tested 34 RNA extractions for the E, N and RdRP genes, reporting 20 positive results for at least one of the three marker genes, with positivity separately confirmed for all the three markers. Control tests to exclude false positivities were successfully accomplished.

CONCLUSION

This is the first evidence that SARS-CoV-2 RNA can be present on PM, thus suggesting a possible use as indicator of epidemic recurrence.

Vibrio parahaemolyticus-specific Halobacteriovorax From Seawater of a Mussel Harvesting Area in the Adriatic Sea: Abundance, Diversity, Efficiency and Relationship With the Prey Natural Level

This research aimed to study the abundance and molecular diversity of Vibrio parahaemolyticus-specific Halobacteriovorax strains isolated from seawater of the Adriatic Sea and the relationship between predator and prey abundances. Moreover, predator efficiency of the Halobacteriovorax isolates toward V. parahaemolyticus and Vibrio cholerae non-O1/O139 strains was tested. V. parahaemolyticus NCTC 10885 was used as primary host for the isolation of Halobacteriovorax from seawater by the plaque assay. Molecular identification was performed by PCR detection of a fragment of the 16S rRNA gene of the Halobacteriovoraceae family members. Moreover, 700 bp PCR products were sequenced and compared between them and to clones described for other sampling sites. Vibrio counts were performed on TCBS agar from 100 ml of filtered water samples and presumptive colonies were confirmed by standard methods. Predatory efficiency of Halobacteriovorax isolates was tested by monitoring abilities of 3-day enrichments to form clear lytic halos on a lawn of Vibrio preys, by the plaque assay. Out of 12 seawater samples monthly collected from June 2017 to May 2018, 10 were positive for V. parahaemolyticus specific Halobacteriovorax with counts ranging from 4 to 1.4 × 10³ PFU per 7.5 ml. No significant relationship was found between Halobacteriovorax and Vibrio abundances. The 16SrRNA sequences of our Halobacteriovorax strains, one for each positive sample, were divided into three lineages. Within the lineages, some sequences had 100% similarity. Sequence similarity between lineages was always <94.5% suggesting that they may therefore well belong to three different species. All Halobacteriovorax isolates had the ability to prey all tested Vibrio strains. Additional research is necessary to assess whether stable strains of Halobacteriovorax are present in the Adriatic Sea and to understand the mechanisms by which Halobacteriovorax may modulate the abundance of V. parahaemolyticus and other vibrios in a complex marine ecosystem.

Impact of nanoplastics on hemolymph immune parameters and microbiota composition in Mytilus galloprovincialis

Ocean contamination by micro- and nanoplastics represents a potential threat to marine biota, from bacterial communities to higher organisms. In this work, the effect of in vivo exposure of Mytilus galloprovincialis to amino modified nanopolystyrene (PS-NH₂) (10 μg/L, 96 h) on hemolymph immune parameters and microbiota composition were investigated. Nanoplastics significantly affected immune parameters (decreased phagocytosis, increased ROS and lysozyme activity, inhibition of NO production). These changes were associated with a shift in hemolymph microbiota composition, with increase in some genera (Arcobacter-like, Psychrobium, Vibrio), and decreases in others (Shewanella, Mycoplasma). The results indicate that exposure to nanoplastics can impact on the microbiome of marine bivalves, and suggest that downregulation of immune defences induced by PS-NH₂ may favour potentially pathogenic bacteria. These data underline how exposure to nanoplastics may represent a potential threat to the complex interplay between innate immunity and host microbiota, thus affecting the homeostatic processes involved in maintenance of organism health.

Searching for SARS-COV-2 on Particulate Matter: A Possible Early Indicator of COVID-19 Epidemic Recurrence

A number of nations were forced to declare a total shutdown due to COVID-19 infection, as extreme measure to cope with dramatic impact of the pandemic, with remarkable consequences both in terms of negative health outcomes and economic loses. However, in many countries a "Phase-2" is approaching and many activities will re-open soon, although with some differences depending on the severity of the outbreak experienced and SARS-COV-2 estimated diffusion in the general population. At the present, possible relapses of the epidemic cannot be excluded until effective vaccines or immunoprophylaxis with human recombinant antibodies will be properly set up and commercialized. COVD-19-related quarantines have triggered serious social challenges, so that decision makers are concerned about the risk of wasting all the sacrifices imposed to the people in these months of quarantine. The availability of possible early predictive indicators of future epidemic relapses would be very useful for public health purposes, and could potentially prevent the suspension of entire national economic systems. On 16 March, a Position Paper launched by the Italian Society of Environmental Medicine (SIMA) hypothesized for the first time a possible link between the dramatic impact of COVID-19 outbreak in Northern Italy and the high concentrations of particulate matter (PM10 and PM2.5) that characterize this area, along with its well-known specific climatic conditions. Thereafter, a survey carried out in the U.S. by the Harvard School of Public Health suggested a strong association between increases in particulate matter concentration and mortality rates due to COVID-19. The presence of SARS-COV-2 RNA on the particulate matter of Bergamo, which is not far from Milan and represents the epicenter of the Italian epidemic, seems to confirm (at least in case of atmospheric stability and high PM concentrations, as it usually occurs in Northern Italy) that the virus can create clusters with the particles and be carried and detected on PM10. Although no assumptions can be made concerning the link between this first experimental finding and COVID-19 outbreak progression or severity, the presence of SARS-COV-2 RNA on PM10 of outdoor air samples in any city of the world could represent a potential early indicator of COVID-19 diffusion. Searching for the viral genome on particulate matter could therefore be explored among the possible strategies for adopting all the necessary preventive measures before future epidemics start.

Identification and functional characterization of the astacidin family of proline-rich host defence peptides (PcAst) from the red swamp crayfish (Procambarus clarkii, Girard 1852)

This study reports the identification of four novel proline-rich antimicrobial peptides (PR-AMP) from the transcriptome of the red swamp crayfish Procambarus clarkii. The newly identified putative peptides (PcAst-1b, -1c, -2 and -3), which are related with the previously identified hemocyte-specific PR-AMP astacidin-1, are encoded by the multi-genic astacidin gene family. The screening of available and proprietary transcriptomes allowed to define the taxonomical range of distribution of this gene family to Astacoidea and Parastacoidea. The antimicrobial properties of three synthetic PcAst peptides (PcAst-1a, -1b/c and -2), were characterized against reference bacteria or multidrug resistant clinical isolates, and their cytotoxicity was evaluated towards human transformed cell lines. The antimicrobial activity ranged from potent and broad-spectrum, in low-salt medium, to poor, whereas it was generally low in full nutrient broth. No significant toxic effects were observed on cultured human cells. RNA-seq data from 12 different tissues indicated a strong specificity for haemocytes under naïve physiological condition, with moderate expression (5-fold lower) in gills. Quantitative real time PCR revealed a rapid (within 2 h) and significant up-regulation of PcAst-1a (Astacidin 1) and PcAst-2 expression in response to LPS injection. Due to the variation in antimicrobial potency and inducibility, the roles of the other astacidins (PcAst-1b, -1c and -3) need to be further investigated to determine their significance to the immune responses of the red swamp crayfish.

A single polyploidization event at the origin of the tetraploid genome of Coffea arabica is responsible for the extremely low genetic variation in wild and cultivated germplasm

The genome of the allotetraploid species Coffea arabica L. was sequenced to assemble independently the two component subgenomes (putatively deriving from C. canephora and C. eugenioides) and to perform a genome-wide analysis of the genetic diversity in cultivated coffee germplasm and in wild populations growing in the center of origin of the species. We assembled a total length of 1.536 Gbp, 444 Mb and 527 Mb of which were assigned to the canephora and eugenioides subgenomes, respectively, and predicted 46,562 gene models, 21,254 and 22,888 of which were assigned to the canephora and to the eugeniodes subgenome, respectively. Through a genome-wide SNP genotyping of 736 C. arabica accessions, we analyzed the genetic diversity in the species and its relationship with geographic distribution and historical records. We observed a weak population structure due to low-frequency derived alleles and highly negative values of Taijma’s D, suggesting a recent and severe bottleneck, most likely resulting from a single event of polyploidization, not only for the cultivated germplasm but also for the entire species. This conclusion is strongly supported by forward simulations of mutation accumulation. However, PCA revealed a cline of genetic diversity reflecting a west-to-east geographical distribution from the center of origin in East Africa to the Arabian Peninsula. The extremely low levels of variation observed in the species, as a consequence of the polyploidization event, make the exploitation of diversity within the species for breeding purposes less interesting than in most crop species and stress the need for introgression of new variability from the diploid progenitors.

First Report of Clinostomum complanatum (Trematoda: Digenea) in European Perch (Perca fluviatilis) from an Italian Subalpine Lake: A Risk for Public Health?

Clinostomum complanatum, a digenean trematode of the Clinostomidae family, is a fish-borne zoonotic parasite responsible for Halzoun syndrome in humans and is transmitted through the consumption of raw or undercooked freshwater fish. Of the total of 112 specimens of European perch (Perca fluviatilis) sampled from a subalpine lake (Lake Endine) in North Italy in 2019, 21 (18.75%) tested positive for encysted metacercariae in the fillet. This study reports the first isolation of C. complanatum in P. fluviatilis and highlights the possible zoonotic risk for consumers, since P. fluviatilis is a food fish used in the traditional local cuisine.

Differential Effects of Dietary Supplementation of Krill Meal, Soybean Meal, Butyrate, and Bactocell® on the Gene Expression of Atlantic Salmon Head Kidney

The head kidney is a key organ that plays a fundamental role in the regulation of the fish immune response and in the maintenance of endocrine homeostasis. Previous studies indicate that the supplementation of exogenous dietary components, such as krill meal (KM), soybean meal (SM), Bactocell® (BA), and butyrate (BU), can have a significant effect on the immune function of the head kidney. The aim of this study was to investigate the differential effect of these four dietary ingredients on the transcriptional profiles of the head kidney of the Atlantic salmon. This study revealed that just a small number of genes were responsive to the feeding regime after a long-term (12 weeks) treatment, and evidenced that the most significant alterations, both in terms of the number of affected genes and magnitude of changes in gene expression, were detectable in the BU- and KM-fed groups compared with controls, while the SM diet had a nearly negligible effect, and BA had no significant effects at all. Most of the differentially expressed genes were involved in the immune response and, in line with data previously obtained from pyloric caeca, major components of the complement system were significantly affected. These alterations were accompanied by an increase in the density of melanomacrophage centers in the KM- and SM-fed group and their reduction in the BU-fed group. While three types of dietary supplements (BU, KM, and SM) were able to produce a significant modulation of some molecular players of the immune system, the butyrate-rich diet was revealed as the one with the most relevant immune-stimulating properties in the head kidney. These preliminary results suggest that further investigations should be aimed towards the elucidation of the potential beneficial effects of butyrate and krill meal supplementation on farmed salmon health and growth performance.

Molecular Diversity of Mytilin-Like Defense Peptides in Mytilidae (Mollusca, Bivalvia)

The CS-αβ architecture is a structural scaffold shared by a high number of small, cationic, cysteine-rich defense peptides, found in nearly all the major branches of the tree of life. Although several CS-αβ peptides involved in innate immune response have been described so far in bivalve mollusks, a clear-cut definition of their molecular diversity is still lacking, leaving the evolutionary relationship among defensins, mytilins, myticins and other structurally similar antimicrobial peptides still unclear. In this study, we performed a comprehensive bioinformatic screening of the genomes and transcriptomes available for marine mussels (Mytilida), redefining the distribution of mytilin-like CS-αβ peptides, which in spite of limited primary sequence similarity maintain in all cases a well-conserved backbone, stabilized by four disulfide bonds. Variations in the size of the alpha-helix and the two antiparallel beta strand region, as well as the positioning of the cysteine residues involved in the formation of the C1-C5 disulfide bond might allow a certain degree of structural flexibility, whose functional implications remain to be investigated. The identification of mytilins in Trichomya and Perna spp. revealed that many additional CS-αβ AMPs remain to be formally described and functionally characterized in Mytilidae, and suggest that a more robust scheme should be used for the future classification of such peptides with respect with their evolutionary origin.

A GM1b/asialo-GM1 oligosaccharide-binding R-type lectin from purplish bifurcate mussels Mytilisepta virgata and its effect on MAP kinases

A 15‐kDa lectin, termed SeviL, was isolated from Mytilisepta virgata (purplish bifurcate mussel). SeviL forms a noncovalent dimer that binds strongly to ganglio‐series GM1b oligosaccharide (Neu5Acɑ2‐3Galβ1‐3GalNAcβ1‐4Galβ1‐4Glc) and its precursor, asialo‐GM1 (Galβ1‐3GalNAcβ1‐4Galβ1‐4Glc). SeviL also interacts weakly with the glycan moiety of SSEA‐4 hexaose (Neu5Acα2‐3Galβ1‐3GalNAcβ1‐3Galα1‐4Galβ1‐4Glc). A partial protein sequence of the lectin was determined by mass spectrometry, and the complete sequence was identified from transcriptomic analysis. SeviL, consisting of 129 amino acids, was classified as an R(icin B)‐type lectin, based on the presence of the QxW motif characteristic of this fold. SeviL mRNA is highly expressed in gills and, in particular, mantle rim tissues. Orthologue sequences were identified in other species of the family Mytilidae, including Mytilus galloprovincialis, from which lectin MytiLec‐1 was isolated and characterized in our previous studies. Thus, mytilid species contain lectins belonging to at least two distinct families (R‐type lectins and mytilectins) that have a common β‐trefoil fold structure but differing glycan‐binding specificities. SeviL displayed notable cytotoxic (apoptotic) effects against various cultured cell lines (human breast, ovarian, and colonic cancer; dog kidney) that possess asialo‐GM1 oligosaccharide at the cell surface. This cytotoxic effect was inhibited by the presence of anti‐asialo‐GM1 oligosaccharide antibodies. With HeLa ovarian cancer cells, SeviL showed dose‐ and time‐dependent activation of kinase MKK3/6, p38 MAPK, and caspase‐3/9. The transduction pathways activated by SeviL via the glycosphingolipid oligosaccharide were triggered apoptosis.

Database

Nucleotide sequence data have been deposited in the GenBank database under accession numbers MK434191, MK434192, MK434193, MK434194, MK434195, MK434196, MK434197, MK434198, MK434199, MK434200, and MK434201.

Dynamics of the Pacific oyster pathobiota during mortality episodes in Europe assessed by 16S rRNA gene profiling and a new target enrichment next-generation sequencing strategy

Infectious agents such as the bacteria Vibrio aestuarianus or Ostreid herpesvirus 1 have been repeatedly associated with dramatic disease outbreaks of Crassostrea gigas beds in Europe. Beside roles played by these pathogens, microbial infections in C. gigas may derive from the contribution of a larger number of microorganisms than previously thought, according to an emerging view supporting the polymicrobial nature of bivalve diseases. In this study, the microbial communities associated with a large number of C. gigas samples collected during recurrent mortality episodes at different European sites were investigated by real-time PCR and 16SrRNA gene-based microbial profiling. A new target enrichment next-generation sequencing protocol for selective capturing of 884 phylogenetic and virulence markers of the potential microbial pathogenic community in oyster tissue was developed allowing high taxonomic resolution analysis of the bivalve pathobiota. Comparative analysis of contrasting C. gigas samples conducted using these methods revealed that oyster experiencing mortality outbreaks displayed signs of microbiota disruption associated with the presence of previously undetected potential pathogenic microbial species mostly belonging to genus Vibrio and Arcobacter. The role of these species and their consortia should be targeted by future studies aiming to shed light on mechanisms underlying polymicrobial infections in C. gigas.

The complex evolutionary history of sulfoxide synthase in ovothiol biosynthesis

Sulfoxide synthases are enzymes involved in the biosynthesis of small sulfur-containing natural products. Their enzymatic activity represents a unique sulfur transfer strategy in nature that is the insertion of a sulfur atom on the imidazole ring of histidine. To date, only two enzymes are known to carry out this function: the sulfoxide synthase EgtB, involved in the biosynthesis of ergothioneine in fungi and bacteria, and the 5-histidylcysteine sulfoxide synthase OvoA, involved in the biosynthesis of ovothiols, found in the eggs and biological fluids of marine invertebrates, some proteobacteria and protists. In particular, ovothiols, thanks to their unique redox properties, are probably the most intriguing marine sulfur-containing molecules. Although they have long been considered as cellular protective molecules, new evidence suggest that their biological activities and ecological role might be more complex than originally thought. Here, we investigate the evolutionary history of OvoA in Metazoa, reporting its monophyletic ancient origins, which could be traced back to the latest common ancestor of Choanozoa. Nevertheless, we show that OvoA is missing in several major extant taxa and we discuss this patchy distribution in the light of the massive genome reduction events documented in Metazoa. We also highlight two interesting cases of secondary acquisition through horizontal gene transfer, which occurred in hydrozoans and bdelloid rotifers. The evolutionary success of this metabolic pathway is probably ascribable to its role in the maintenance of cellular redox homeostasis, which enables organisms to survive in different environmental niches.

Identification, molecular characterization and functional analysis of interleukin (IL)-2 and IL-2like (IL-2L) cytokines in sea bass (Dicentrarchus labrax L.)

In mammals, interleukin (IL)-2, initially known as a T-cell grow factor, is an immunomodulatory cytokine involved in the proliferation of T cells upon antigen activation. In bony fish, some IL-2 orthologs have been identified, but, recently, an additional IL-2like (IL-2L) gene has been found. In this paper, we report the presence of these two divergent IL-2 isoforms in sea bass (Dicentrarchus labrax L.). Genomic analyses revealed that they originated from a gene duplication event, as happened in most percomorphs. These two IL-2 paralogs show differences in the amino acid sequence and in the exon 4 size, and these features could be an indication that they bind preferentially to different specific IL-2 receptors. Sea bass IL-2 paralogs are highly expressed in gut and spleen, which are tissues and organs involved in fish T cell immune functions, and the two cytokines could be up-regulated by both PHA stimulation and vaccination with a bacterial vaccine, with IL-2L being more inducible. To investigate the functional activities of sea bass IL-2 and IL-2L we produced the corresponding recombinant molecules in E. coli and used them to in vitro stimulate HK and spleen leukocytes. IL-2L is able to up-regulate the expression of markers related to different T cell subsets (Th1, Th2 and Th17) and to Treg cells in HK, whereas it has little effect in spleen. IL-2 is not active on these markers in HK, but shows an effect on Th1 markers in spleen. Finally, the stimulation with recombinant IL-2 and IL-2L is also able to induce in vitro proliferation of HK- and spleen-derived leukocytes. In conclusion, we have demonstrated that sea bass possess two IL-2 paralogs that likely have an important role in regulating T cell development in this species and that show distinct bioactivities.

De novo gonad transcriptome analysis of the common littoral shrimp Palaemon serratus: novel insights into sex-related genes

BACKGROUND

The common littoral shrimp Palaemon serratus is an economically important decapod resource in some European communities. Aquaculture practices prevent the genetic deterioration of wild stocks caused by overfishing and at the same time enhance the production. The biotechnological manipulation of sex-related genes has the proved potential to improve the aquaculture production but the scarcity of genomic data about P. serratus hinders these applications. RNA-Seq analysis has been performed on ovary and testis samples to generate a reference gonadal transcriptome. Differential expression analyses were conducted between three ovary and three testis samples sequenced by Illumina HiSeq 4000 PE100 to reveal sex-related genes with sex-biased or sex-specific expression patterns.

RESULTS

A total of 224.5 and 281.1 million paired-end reads were produced from ovary and testis samples, respectively. De novo assembly of ovary and testis trimmed reads yielded a transcriptome with 39,186 transcripts. The 29.57% of the transcriptome retrieved at least one annotation and 11,087 differentially expressed genes (DEGs) were detected between ovary and testis replicates. Six thousand two hundred seven genes were up-regulated in ovaries meanwhile 4880 genes were up-regulated in testes. Candidate genes to be involved in sexual development and gonadal development processes were retrieved from the transcriptome. These sex-related genes were discussed taking into account whether they were up-regulated in ovary, up-regulated in testis or not differentially expressed between gonads and in the framework of previous findings in other crustacean species.

CONCLUSIONS

This is the first transcriptome analysis of P. serratus gonads using RNA-Seq technology. Interesting findings about sex-related genes from an evolutionary perspective (such as Dmrt1) and for putative future aquaculture applications (Iag or vitellogenesis genes) are reported here. We provide a valuable dataset that will facilitate further research into the reproductive biology of this shrimp.

Extensive Tandem Duplication Events Drive the Expansion of the C1q-Domain-Containing Gene Family in Bivalves

C1q-domain-containing (C1qDC) proteins are rapidly emerging as key players in the innate immune response of bivalve mollusks. Growing experimental evidence suggests that these highly abundant secretory proteins are involved in the recognition of microbe-associated molecular patterns, serving as lectin-like molecules in the bivalve proto-complement system. While a large amount of functional data concerning the binding specificity of the globular head C1q domain and on the regulation of these molecules in response to infection are quickly accumulating, the genetic mechanisms that have led to the extraordinary lineage-specific expansion of the C1qDC gene family in bivalves are still largely unknown. The analysis of the chromosome-scale genome assembly of the Eastern oyster Crassostrea virginica revealed that the 476 oyster C1qDC genes, far from being uniformly distributed along the genome, are located in large clusters of tandemly duplicated paralogs, mostly found on chromosomes 7 and 8. Our observations point out that the evolutionary process behind the development of a large arsenal of C1qDC lectin-like molecules in marine bivalves is still ongoing and likely based on an unequal crossing over.

Expansion and loss events characterized the occurrence of MIF-like genes in bivalves

Macrophage migration inhibitory factor (MIF) dynamically connects innate and adaptive immune systems in vertebrate animals, allowing highly orchestrated systemic responses to various insults. The occurrence of MIF-like genes in non-vertebrate organisms suggests its origin from an ancestral metazoan gene, whose function is still a matter of debate. In the present work, by analyzing available genomic and transcriptomic data from bivalve mollusks, we identified 137 MIF-like sequences, which were classified into three types, based on phylogeny and conservation of key residues: MIF, D-DT, and the lineage-specific type MDL. Comparative genomics revealed syntenic conservation of homologous genes at the family level, the loss of D-DT in the Ostreidae family as well as the expansion of MIF-like genes in the Mytilidae family, possibly underpinning the neofunctionalization of duplicated gene copies. In M. galloprovincialis, MIF and one D-DT were mostly expressed in haemocytes and mantle rim of untreated animals, while D-DT paralogs often showed very limited expression, suggesting an accessory role or their persistence as relict genes.

Immune Tolerance in Mytilus galloprovincialis Hemocytes After Repeated Contact With Vibrio splendidus

Mediterranean mussels (Mytilus galloprovincialis) are sessile filter feeders that live in close contact with numerous marine microorganisms. As is the case in all invertebrates, mussels lack an adaptive immune system, but they respond to pathogens, injuries or environmental stress in a very efficient manner. However, it is not known if they are able to modify their immune response when they reencounter the same pathogen. In this work, we studied the transcriptomic response of mussel hemocytes before and after two consecutive sublethal challenges with Vibrio splendidus. The first exposure significantly regulated genes related to inflammation, migration and response to bacteria. However, after the second exposure, the differentially expressed genes were related to the control and inhibition of ROS production and the resolution of the inflammatory response. Our results also show that the second injection with V. splendidus led to changes at the transcriptional (control of the expression of pro-inflammatory transcripts), cellular (shift in the hemocyte population distribution), and functional levels (inhibition of ROS production). These results suggest that a modified immune response after the second challenge allowed the mussels to tolerate rather than fight the infection, which minimized tissue damage.