Systemic immune responses do not affect significant immune responses in the skin.

Fish skin plays an important role in defending against pathogens in water, primarily through the secretion of skin mucus containing various immune-related factors. Local immune responses in the skin activate systemic immune responses by inflammatory cytokines. However, it remains unclear whether immune responses in the skin occur after systemic immune responses caused by pathogen invasion into the fish body. This study aimed to clarify the relationship between systemic immune responses and skin responses after intraperitoneal injection of formalin-killed cells (FKC) of Vibrio anguillarum. Although systemic inflammatory responses were observed in the spleen after injection, expression changes in the skin did not show significant differences. In contrast, expression of hemoglobin subunit genes significantly increased in the skin after FKC injection, suggesting that erythrocytes infiltrate extravascularly.

Investigating the impact of chlorine dioxide in shrimp-rearing water on stomach microbiome, gill transcriptome, and infection-related mortality in shrimp.

AIMS: This study aimed to assess the effects of chlorine dioxide (ClO2) in water on whiteleg shrimp Penaeus vannamei, evaluating its impact on stomach microbiota, gill transcriptome, and pathogens. METHODS AND RESULTS: ClO2 was added to the aquarium tanks containing the shrimp. The application of ClO2 to rearing water was lethal to shrimp at concentrations above 1.2 ppm. On the other hand, most of the shrimp survived at 1.0 ppm of ClO2. Microbiome analysis showed that ClO2 administration at 1.0 ppm significantly reduced the α-diversity of bacterial community composition in the shrimp stomach, and this condition persisted for at least 7 days. Transcriptome analysis of shrimp gill revealed that ClO2 treatment caused massive change of the gene expression profile including stress response genes. However, after 7 days of the treatment, the gene expression profile was similar to that of shrimp in the untreated control group, suggesting a recovery to the normal state. This 1.0-ppm ClO2 significantly reduced shrimp mortality in artificial challenges with an AHPND-causing Vibrio parahaemolyticus and white spot syndrome virus, which were added to rearing water. CONCLUSIONS: The use of ClO2 at appropriate concentrations effectively eliminates a significant portion of the bacteria in the shrimp stomach and pathogens in the water. The results of this study provide fundamental knowledge on the disinfection of pathogens in water using ClO2 and the creation of semi germ-free shrimp, which has significantly decreased microbiome in the stomach.

De Novo Assembly and Annotation of the Siganus fuscescens (Houttuyn, 1782) Genome: Marking a Pioneering Advance for the Siganidae Family.

This study presents the first draft genome of Siganus fuscescens, and thereby establishes the first whole-genome sequence for a species in the Siganidae family. Leveraging both long and short read sequencing technologies, i.e., Oxford Nanopore and Illumina sequencing, we successfully assembled a mitogenome spanning 16.494 Kb and a first haploid genome encompassing 498 Mb. The assembled genome accounted for a 99.6% of the estimated genome size and was organized into 164 contigs with an N50 of 7.2 Mb. This genome assembly showed a GC content of 42.9% and a high Benchmarking Universal Single-Copy Orthologue (BUSCO) completeness score of 99.5% using actinopterygii_odb10 lineage, thereby meeting stringent quality standards. In addition to its structural aspects, our study also examined the functional genomics of this species, including the intricate capacity to biosynthesize long-chain polyunsaturated fatty acids (LC-PUFAs) and secrete venom. Notably, our analyses revealed various repeats elements, which collectively constituted 17.43% of the genome. Moreover, annotation of 28,351 genes uncovered both shared genetic signatures and those that are unique to S. fuscescens. Our assembled genome also displayed a moderate prevalence of gene duplication compared to other fish species, which suggests that this species has a distinctive evolutionary trajectory and potentially unique functional constraints. Taken altogether, this genomic resource establishes a robust foundation for future research on the biology, evolution, and the aquaculture potential of S. fuscescens.

Molecular characterization of a short-chained pentraxin gene from kuruma shrimp Marsupenaeus japonicus hemocytes.

Pentraxins (PTXs) are a family of pattern recognition proteins (PRPs) that play a role in pathogen recognition during infection via pathogen-associated molecular patterns (PAMPs). Here, we characterized a short-chained pentraxin isolated from kuruma shrimp (Marsupenaeus japonicus) hemocytes (MjPTX). MjPTX contains the pentraxin signature HxCxS/TWxS (where x can be any amino acid), although the second conserved residue of this signature differed slightly (L instead of C). In the phylogenetic analysis, MjPTX clustered closely with predicted sequences from crustaceans (shrimp, lobster, and crayfish) displaying high sequence identities exceeding 52.67 %. In contrast, MjPTX showed minimal sequence identity when compared to functionally similar proteins in other animals, with sequence identities ranging from 20.42 % (mouse) to 28.14 % (horseshoe crab). MjPTX mRNA transcript levels increased significantly after artificial infection with Vibrio parahaemolyticus (48 h), White Spot Syndrome Virus (72 h) and Yellow Head Virus (24 and 48 h). Assays done in vitro revealed that recombinant MjPTX (rMjPTX) has an ability to agglutinate Gram-negative and Gram-positive bacteria and to bind microbial polysaccharides and bacterial suspensions in the presence of Ca2+. Taken together, our results suggest that MjPTX functions as a classical pattern recognition protein in the presence of calcium ions, that is capable of binding to specific moieties present on the surface of microorganisms and facilitating their clearance.

Metagenome-assembled genomes of three Hepatoplasmataceae provide insights into isopod-mollicute symbiosis.

The digestive organs of terrestrial isopods harbour bacteria of the recently proposed mollicute family Hepatoplasmataceae. The only complete genome available so far for Hepatoplasmataceae is that of 'Candidatus Hepatoplasma crinochetorum'. The scarcity of genome sequences has hampered our understanding of the symbiotic relationship between isopods and mollicutes. Here, we present four complete metagenome-assembled genomes (MAGs) of uncultured Hepatoplasmataceae members identified from shotgun sequencing data of isopods. We propose genomospecies names for three MAGs that show substantial sequence divergence from any previously known Hepatoplamsataceae members: 'Candidatus Tyloplasma litorale' identified from the semiterrestrial isopod Tylos granuliferus, 'Candidatus Hepatoplasma vulgare' identified from the common pill bug Armadillidium vulgare, and 'Candidatus Hepatoplasma scabrum' identified from the common rough woodlouse Porcellio scaber. Phylogenomic analysis of 155 mollicutes confirmed that Hepatoplasmataceae is a sister clade of Metamycoplasmataceae in the order Mycoplasmoidales. The 16S ribosomal RNA gene sequences and phylogenomic analysis showed that 'Candidatus Tyloplasma litorale' and other semiterrestrial isopod-associated mollicutes represent the placeholder genus 'g_Bg2' in the r214 release of the Genome Taxonomy Database, warranting their assignment to a novel genus. Our analysis also revealed that Hepatoplasmataceae lack major metabolic pathways but has a likely intact type IIA CRISPR-Cas9 machinery. Although the localization of the Hepatoplasmatacae members have not been verified microscopically in this study, these genomic characteristics are compatible with the idea that these mollicutes have an ectosymbiotic lifestyle with high nutritional dependence on their host, as has been demonstrated for other members of the family. We could not find evidence that Hepatoplasmataceae encode polysaccharide-degrading enzymes that aid host digestion. If they are to provide nutritional benefits, it may be through extra-copy nucleases, peptidases, and a patatin-like lipase. Exploration of potential host-symbiont interaction-associated genes revealed large, repetitive open reading frames harbouring beta-sandwich domains, possibly involved with host cell adhesion. Overall, genomic analyses suggest that isopod-mollicute symbiosis is not characterized by carbohydrate degradation, and we speculate on their potential role as defensive symbionts through spatial competition with pathogens to prevent infection.

Integrase-associated niche differentiation of endogenous large DNA viruses in crustaceans.

IMPORTANCE: Crustacean genomes harbor sequences originating from a family of large DNA viruses called nimaviruses, but it is unclear why they are present. We show that endogenous nimaviruses selectively insert into repetitive sequences within the host genome, and this insertion specificity was correlated with different types of integrases, which are DNA recombination enzymes encoded by the nimaviruses themselves. This suggests that endogenous nimaviruses have colonized various genomic niches through the acquisition of integrases with different insertion specificities. Our results point to a novel survival strategy of endogenous large DNA viruses colonizing the host genomes. These findings may clarify the evolution and spread of nimaviruses in crustaceans and lead to measures to control and prevent the spread of pathogenic nimaviruses in aquaculture settings.

Red sea bream iridovirus infection downregulates inflammation-related genes in the spleen of rock bream (Oplegnathus fasciatus).

This study investigated the kinetics of red sea bream iridovirus and host gene expression during infection in rock bream (Oplegnathus fasciatus), a species highly sensitive to the virus. After intraperitoneal injection of the viral solution at 104 TCID50/fish, the viral genome copy number in the spleen was 104.7 ± 0.2 and 105.9 ± 0.4 copies/µg DNA at 3 and 5 days post-injection (dpi), respectively. Using transcriptomic analyses via MiSeq, viral gene transcripts were detected at 3 and 5 dpi. Six genes including RING-finger domain-containing protein and laminin-type epidermal growth factor-like domain genes were significantly expressed at 5 dpi. Further, 334 host genes were differentially expressed compared with those before infection. Genes were clustered into four groups based on their expression profiles. Interferon-stimulated genes were more prevalent in groups showing upregulation at 5 dpi and 3 and 5 dpi. In contrast, the group showing downregulation at 3 dpi included inflammation-related genes, such as granzyme and eosinophil peroxidase genes. Downregulation of certain inflammation-related genes may contribute to the susceptibility of this fish to the virus.

A novel liver-specific immunoglobulin heavy chain-like gene in a cartilaginous fish.

We identified a novel immunoglobulin (Ig) heavy chain-like gene (tsIgH) expressed in the liver of the banded houndshark Triakis scyllium by preliminary transcriptomic analysis. The tsIgH gene showed less than 30% of amino acid identities to Ig genes of the shark. The gene encodes one variable domain (VH) and three conserved domains (CH1-CH3) with a predicted signal peptide. Interestingly, this protein has only one cysteine residue in a linker region between VH and CH1 other than those required for the formation of the immunoglobulin domain. Genome sequencing revealed that each of the domains was encoded by a corresponding single exon, and the exon-intron structures of the homologues are conserved in the other cartilaginous fishes. By RT-qPCR analysis, the transcript of the tsIgH gene was observed only in the liver, while that of the IgM was mainly detected in the epigonal organ, liver, and spleen. The novel Ig-heavy chain-like gene in cartilaginous fish may provide new clues to the evolution of immunoglobulin genes.

Genomic characterization and identification of virulence-related genes in Vibrio nigripulchritudo isolated from white leg shrimp Penaeus vannamei.

Vibrio nigripulchritudo causes vibriosis in penaeid shrimps. Here, we used Illumina and Nanopore sequencing technologies to sequence the genomes of three of its strains (TUMSAT-V. nig1, TUMSAT-V. nig2, and TUMSAT-V. nig3) to explore opportunities for disease management. Putative virulence factors and mobile genetic elements were detected while evaluating the phylogenetic relationship of each isolated strain. The genomes consisted of two circular chromosomes (I and II) plus one or two plasmids. The size of chromosome I ranged from 4.02 to 4.07 Mb with an average GC content of 46%, while the number of predicted CDSs ranged from 3563 to 3644. The size of chromosome II ranged from 2.16 to 2.18 Mb, with an average GC content of 45.5%, and the number of predicted CDSs ranged from 1970 to 1987. Numerous virulence genes were identified related to adherence, antiphagocytosis, chemotaxis, motility, iron uptake, quorum sensing, secretion systems, and toxins in all three genomes. Higher numbers of prophages and genomic islands found in TUMSAT-V. nig1 suggest that the strain has experienced numerous horizontal gene transfer events. The presence of antimicrobial resistance genes suggests that the strains have multidrug resistance. Comparative genomic analysis showed that all three strains belonged to the same clade.

Genome sequence of Chionoecetes opilio bacilliform virus, a nimavirus infecting the snow crab Chionoecetes opilio.

Nimaviridae (class Naldaviricetes) is a family of double-stranded DNA viruses infecting crustaceans, with the only officially recognized representative being white spot syndrome virus (WSSV). Chionoecetes opilio bacilliform virus (CoBV) was isolated as the causative agent of milky hemolymph disease in the snow crab Chionoecetes opilio, an economically important crustacean in the northwestern Pacific. Here, we present the complete genome sequence of CoBV and show that it is unambiguously a nimavirus. The CoBV genome is a 240-kb circular DNA molecule with 40% GC content that encodes 105 proteins, including 76 WSSV orthologs. Phylogenetic analysis based on eight naldaviral core genes established that CoBV is a member of the family Nimaviridae. The availability of the CoBV genome sequence provides a deeper understanding of CoBV pathogenicity and nimavirus evolution.

Identification of a rabbit Ig light chain recombinant protein bound to serum immunoglobulins from different marine fish species.

The structures of fish serum immunoglobulin differ among different fish species. In this study, we accidently isolated a rabbit immunoglobulin (Ig) light chain bound to serum immunoglobulin from different marine fish species using phage display. Fish Ig was separated using a protein A column. The phage library was generated from variable regions of rabbit spleen B cells immunized with bluefin tuna Thunnus orientalis Ig. Fish Ig-specific phages were enriched using two rounds of bio-panning with yellowtail Seriola quinqueradiata serum Ig, followed by two rounds of bio-panning with red seabream Pagrus major serum Ig. The enriched phages demonstrated an increase in binding specificity to the tuna, yellowtail, and red seabream Igs compared to the phages listed in the unpanned library. A recombinant protein of a single clonal phage, which encodes the rabbit Ig light chain, was produced, and the binding specificities to fish Igs were analyzed using enzyme-linked immunosorbent assay (ELISA) and western blotting. The recombinant protein exhibited binding properties to fish Igs in the ELISA. However, the recombinant protein that bound to serum protein(s), but not IgM, was detected via western blotting. The recombinant protein may provide a novel information on the common structural feature in the fish immunoglobulins.

Whole-genome analysis of red sea bream iridovirus spread in 2021 in Japan provided epidemiological and viral traits insight.

Red sea bream iridovirus (RSIV) is the pathogen that causes red sea bream iridoviral disease. It causes a huge loss to the Japanese aquaculture industry. In 2021, outbreaks of red sea bream iridovirus occurred in South Japan. This study analysed nine whole-genome sequences of RSIV isolated in Oita and Ehime Prefectures in 2021 using a short-read next-generation sequencer. Nine isolates had highly uniform sequences, and there was no variant depending on locations or host species. Phylogenetic analyses with other reported megalocytivirus isolates showed that RSIV isolated in 2021 was genetically different from RSIV previously isolated in Oita and Ehime Prefectures in 2017-2019. These results suggest that RSIV isolated in Oita and Ehime Prefectures in 2021 might spread from a common ancestor different from the recent one. Additionally, it was found that RSIV isolated in 2021 had sequence mutations on protein-coding sequences that may be involved in viral pathogenicity and infectivity.

Comparative genome analyses of five Vibrio penaeicida strains provide insights into their virulence-related factors.

Vibrio penaeicida (family Vibrionaceae) is an important bacterial pathogen that affects Japanese shrimp aquaculture. Only two whole-genome sequences of V. penaeicida are publicly available, which has hampered our understanding of the pathogenesis of shrimp vibriosis caused by this bacterium. To gain insight into the genetic features, evolution and pathogenicity of V. penaeicida, we sequenced five V. penaeicida strains (IFO 15640T, IFO 15641, IFO 15642, TUMSAT-OK1 and TUMSAT-OK2) and performed comparative genomic analyses. Virulence factors and mobile genetic elements were detected. Furthermore, average nucleotide identities (ANIs), clusters of orthologous groups and phylogenetic relationships were evaluated. The V. penaeicida genome consists of two circular chromosomes. Chromosome I sizes ranged from 4.1 to 4.3 Mb, the GC content ranged from 43.9 to 44.1 %, and the number of predicted protein-coding sequences (CDSs) ranged from 3620 to 3782. Chromosome II sizes ranged from 2.2 to 2.4 Mb, the GC content ranged from 43.5 to 43.8 %, and the number of predicted CDSs ranged from 1992 to 2273. All strains except IFO 15641 harboured one plasmid, having sizes that ranged from 150 to 285 kb. All five genomes had typical virulence factors, including adherence, anti-phagocytosis, flagella-related proteins and toxins (repeats-in-toxin and thermolabile haemolysin). The genomes also contained factors responsible for iron uptake and the type II, IV and VI secretion systems. The genome of strain TUMSAT-OK2 tended to encode more prophage regions than the other strains, whereas the genome of strain IFO 15640T had the highest number of regions encoding genomic islands. For comparative genome analysis, we used V. penaeicida (strain CAIM 285T) as a reference strain. ANIs between strain CAIM 285T and the five V. penaeicida strains were >95 %, which indicated that these strains belong to the same species. Orthology cluster analysis showed that strains TUMSAT-OK1 and TUMSAT-OK2 had the greatest number of shared gene clusters, followed by strains CAIM 285T and IFO 15640T. These strains were also the most closely related to each other in a phylogenetic analysis. This study presents the first comparative genome analysis of V. penaeicida and these results will be useful for understanding the pathogenesis of this bacterium.

Genome and transcriptome assemblies of the kuruma shrimp, Marsupenaeus japonicus.

The kuruma shrimp Marsupenaeus japonicus (order Decapoda, family Penaeidae) is an economically important crustacean that occurs in shallow, warm seas across the Indo-Pacific. Here, using a combination of Illumina and Oxford Nanopore Technologies platforms, we produced a draft genome assembly of M. japonicus (1.70 Gbp; 18,210 scaffolds; scaffold N50 = 234.9 kbp; 34.38% GC, 93.4% BUSCO completeness) and a complete mitochondrial genome sequence (15,969 bp). As with other penaeid shrimp genomes, the M. japonicus genome is extremely rich in simple repeats, which occupies 27.4% of the assembly. A total of 26,381 protein-coding gene models (94.7% BUSCO completeness) were predicted, of which 18,005 genes (68.2%) were assigned functional description by at least one method. We also produced an Illumina-based transcriptome shotgun assembly (40,991 entries; 93.0% BUSCO completeness) and a PacBio Iso-Seq transcriptome assembly (25,415 entries; 67.5% BUSCO completeness). We envision that the M. japonicus genome and transcriptome assemblies will serve as useful resources for the basic research, fisheries management, and breeding programs of M. japonicus.

Genome Sequence of Lymphocystis Disease Virus 2 LCDV-JP_Oita_2018, Isolated from a Diseased Japanese Flounder (Paralichthys olivaceus) in Japan.

Here, we present the genome sequence of lymphocystis disease virus 2 LCDV-JP_Oita_2018 (genus Lymphocystivirus, family Iridoviridae), which was isolated from a diseased Japanese flounder (Paralichthys olivaceus) in Japan. The LCDV-JP_Oita_2018 genome was assembled into a circular contig of 186,627 bp, with 140 predicted protein-coding genes and a GC content of 27%.

Development of single nucleotide polymorphism (SNP) application for detection and genotyping of RSIV-type megalocytiviruses.

Red sea bream iridovirus (RSIV) belonging to the genus Megalocytivirus of the family Iridoviridae is the cause of serious mass mortality of cultured marine fishes. RSIV-type megalocytiviruses show extremely high nucleotide sequence identities. Thus, epidemiological studies on this virus are limited. This study developed two primer sets amplifying the regions possessing single nucleotide polymorphism (SNP) to determine the relationships and divergence of RSIV-type megalocytiviruses isolated from cultured marine fishes in Japan. The two regions were designed according to the genome sequences of the representative RSIV genotype II of megalocytivirus members in GenBank. The SNP 1 and 2 regions have sequences homologous to hypothetical protein ORF 24 and ORF 31, respectively, of RSIV (accession no. AP017456.1). By sequencing the regions, 53 polymorphic sites were identified. The phylogenetic analysis of 25 RSIV-type megalocytivirus isolates, classified into RSIV cluster, was clustered into eight haplotypes (seven haplotypes from Oita, two haplotypes from Ehime, and one haplotype shared between Oita and Ehime). These findings suggested that SNP in the RSIV genome is a powerful application for the detection and identification of RSIV-type megalocytiviruses.

Genome Sequence of Vibrio nigripulchritudo Strain TUMSAT-TG-2018, Isolated from Diseased Pacific White Shrimp, Litopenaeus vannamei.

The Gram-negative bacterium Vibrio nigripulchritudo is an important shrimp pathogen. Here, we present the genome sequence of Vibrio nigripulchritudo TUMSAT-TG-2018, which was isolated from a diseased Pacific white shrimp (Litopenaeus vannamei). The assembly totaled 6.8 Mbp, consisting of two chromosomes and four plasmids.

Draft Genome Sequences of Vibrio atypicus Strains DSM 25292T and TUMSAT1.

Vibrio atypicus is a Gram-negative bacterium associated with the digestive tract of penaeid shrimp. Here, we present the draft genome sequences of two V. atypicus strains, DSM 25292T and TUMSAT1. Both assemblies were approximately 4.8 Mbp long, with GC contents of around 44%.

Anti-PirA-like toxin immunoglobulin (IgY) in feeds passively immunizes shrimp against acute hepatopancreatic necrosis disease.

Acute hepatopancreatic necrosis disease (AHPND), caused by a toxin-producing Vibrio parahaemolyticus strain, has become a serious threat to shrimp aquaculture. The need to regulate antibiotic use prompted the development of alternative ways to treat infections in aquaculture including the use of chicken egg yolk immunoglobulin (IgY) for passive immunization. This study evaluated the protective effect of IgY against AHPND infection in Litopenaeus vannamei (Boone). IgY was isolated from eggs laid by hens immunized with recombinant PirA-like (rPirA) and PirB-like (rPirB) toxins. Whole-egg powders having IgY specific to rPirA (anti-PirA-IgY) and rPirB (anti-PirB-IgY) and IgY from non-immunized hen (control-IgY) were mixed with basal diets at 20% concentrations and used to prefeed shrimp 3 days before the bacterial challenge test. Survival rates of the challenged shrimp fed the anti-PirA-IgY, anti-PirB-IgY and control-IgY diets were 86%, 14% and 0%, respectively. Only the feed containing anti-PirA-IgY protected shrimp against AHPND. Increasing the concentration of rPirA antigen to immunize hens and lowering the amount of egg powder in feeds to 10% consistently showed higher survival rates in shrimp fed with anti-PirA-IgY (87%) compared with the control (12%). These results confirm that addition of anti-PirA-IgY in feeds could be an effective prophylactic method against AHPND infection in shrimp.

Crustacean Genome Exploration Reveals the Evolutionary Origin of White Spot Syndrome Virus.

White spot syndrome virus (WSSV) is a crustacean-infecting, double-stranded DNA virus and is the most serious viral pathogen in the global shrimp industry. WSSV is the sole recognized member of the family Nimaviridae, and the lack of genomic data on other nimaviruses has obscured the evolutionary history of WSSV. Here, we investigated the evolutionary history of WSSV by characterizing WSSV relatives hidden in host genomic data. We surveyed 14 host crustacean genomes and identified five novel nimaviral genomes. Comparative genomic analysis of Nimaviridae identified 28 "core genes" that are ubiquitously conserved in Nimaviridae; unexpected conservation of 13 uncharacterized proteins highlighted yet-unknown essential functions underlying the nimavirus replication cycle. The ancestral Nimaviridae gene set contained five baculoviral per os infectivity factor homologs and a sulfhydryl oxidase homolog, suggesting a shared phylogenetic origin of Nimaviridae and insect-associated double-stranded DNA viruses. Moreover, we show that novel gene acquisition and subsequent amplification reinforced the unique accessory gene repertoire of WSSV. Expansion of unique envelope protein and nonstructural virulence-associated genes may have been the key genomic event that made WSSV such a deadly pathogen.IMPORTANCE WSSV is the deadliest viral pathogen threatening global shrimp aquaculture. The evolutionary history of WSSV has remained a mystery, because few WSSV relatives, or nimaviruses, had been reported. Our aim was to trace the history of WSSV using the genomes of novel nimaviruses hidden in host genome data. We demonstrate that WSSV emerged from a diverse family of crustacean-infecting large DNA viruses. By comparing the genomes of WSSV and its relatives, we show that WSSV possesses an expanded set of unique host-virus interaction-related genes. This extensive gene gain may have been the key genomic event that made WSSV such a deadly pathogen. Moreover, conservation of insect-infecting virus protein homologs suggests a common phylogenetic origin of crustacean-infecting Nimaviridae and other insect-infecting DNA viruses. Our work redefines the previously poorly characterized crustacean virus family and reveals the ancient genomic events that preordained the emergence of a devastating shrimp pathogen.