Application of a new real-time polymerase chain reaction assay for surveillance studies of lymphocystis disease virus in farmed gilthead seabream

Immune response of DNA vaccinated-gilthead seabream (Sparus aurata) against LCDV-Sa infection: relevance of the inflammatory process

Lymphocystis disease is one of the main viral pathologies affecting cultured gilthead seabream (Sparus aurata) in the Mediterranean region. Recently, we have developed a DNA vaccine based on the major capsid protein (MCP) of the Lymphocystis disease virus 3 (LCDV-Sa). The immune response triggered by either LCDV-Sa infection or vaccination have been previously studied and seem to be highly related to the modulation of the inflammatory and the IFN response. However, a comprehensive evaluation of immune-related gene expression in vaccinated fish after viral infection to identify immunogenes involved in vaccine-induced protection have not been carried out to date. The present study aimed to fulfill this objective by analyzing samples of head-kidney, spleen, intestine, and caudal fin from fish using an OpenArray® platform containing targets related to the immune response of gilthead seabream. The results obtained showed an increase of deregulated genes in the hematopoietic organs between vaccinated and non-vaccinated fish. However, in the intestine and fin, the results showed the opposite trend. The global effect of fish vaccination was a significant decrease (p<0.05) of viral replication in groups of fish previously vaccinated, and the expression of the following immune genes related to viral recognition (tlr9), humoral and cellular response (rag1 and cd48), inflammation (csf1r, elam, il1β, and il6), antiviral response (isg15, mx1, mx2, mx3), cell-mediated cytotoxicity (nccrp1), and apoptosis (prf1). The exclusive modulation of the immune response provoked by the vaccination seems to control the progression of the infection in the experimentally challenged gilthead seabream.

Advances in Viral Aquatic Animal Disease Knowledge: The Molecular Methods' Contribution

Aquaculture is the fastest-growing food-producing sector, with a global production of 122.6 million tonnes in 2020. Nonetheless, aquatic animal production can be hampered by the occurrence of viral diseases. Furthermore, intensive farming conditions and an increasing number of reared fish species have boosted the number of aquatic animals' pathogens that researchers have to deal with, requiring the quick development of new detection and study methods for novel unknown pathogens. In this respect, the molecular tools have significantly contributed to investigating thoroughly the structural constituents of fish viruses and providing efficient detection methods. For instance, next-generation sequencing has been crucial in reassignment to the correct taxonomic family, the sturgeon nucleo-cytoplasmic large DNA viruses, a group of viruses historically known, but mistakenly considered as iridoviruses. Further methods such as in situ hybridisation allowed objectifying the role played by the pathogen in the determinism of disease, as the cyprinid herpesvirus 2, ostreid herpesvirus 1 and betanodaviruses. Often, a combination of molecular techniques is crucial to understanding the viral role, especially when the virus is detected in a new aquatic animal species. With this paper, the authors would critically revise the scientific literature, dealing with the molecular techniques employed hitherto to study the most relevant finfish and shellfish viral pathogens.

New Insights into Lymphocystis Disease Virus Genome Diversity

Lymphocystis disease viruses (LCDVs) are viruses that infect bony fish which has been found in different locations across the globe. Four virus species have been classified by the International Committee on Taxonomy of Viruses (ICTV), despite remarkable discrepancies in genome size. Whole genome sequencing and phylogenetic analysis of LCDVs from wild fish from the North Sea and partial sequences from gilthead sea bream of an aquafarm located in the Aegean Sea in Turkey confirm that the LCDV1 genome at 100 kb is approximately half the size of the genomes of LCDV2-4. Since the fish species, of which LCDV1 was isolated, differ taxonomically at the order level, co-speciation can be excluded as the driver of the adaptation of the genome of this nucleocytoplasmic large DNA virus, but may represent an adaptation to the lifestyle of this demersal fish in the northeast Atlantic.

FELASA-AALAS Recommendations for Monitoring and Reporting of Laboratory Fish Diseases and Health Status, with an Emphasis on Zebrafish (Danio Rerio)

The exchange of fish for research may expose an aquatic laboratory to pathogen contamination as incoming fish can introduce bacteria, fungi, parasites, and viruses capable of affecting both experimental results and fish and personnel health and welfare. To develop risk mitigation strategies, FELASA and AALAS established a joint working group to recommend good practices for health monitoring of laboratory fish. The recommendations address all fish species used for research, with a particular focus on zebrafish (Danio rerio). First, the background of the working group and key definitions are provided. Next, fish diseases of high impact are described. Third, recommendations are made for health monitoring of laboratory fishes. The recommendations emphasize the importance of daily observation of the fish and strategies to determine fish colony health status. Finally, report templates are proposed for historical screening data and aquatic facility description to facilitate biohazard risk assessment when exchanging fish.

Initial Evidence That Gilthead Seabream (Sparus aurata L.) Is a Host for Lymphocystis Disease Virus Genotype I

Simple Summary

Nodular lesions were observed on the skin and fins of 95% of one and a half million juvenile gilthead seabreams cultured in Egypt, shortly after importation from Europe. We undertook a study to describe the clinical disease course, identify the causative agent, and investigate the origin of the causative agent. Preliminary diagnosis based on gross lesions and postmortem examination suggested lymphocystis disease caused by lymphocystis disease virus, Lymphocystivirus, Iridoviridae. Histopathological and ultrastructural pictures were typical of lymphocystis disease virus infections. Polymerase Chain Reaction followed by sequencing and phylogenetic analysis of the major capsid protein gene demonstrated the presence of lymphocystis disease virus genotype I, originally associated with lymphocystis disease in Northern European countries, with 99.7% and 100% nucleotide and deduced amino acid identity values, respectively. Lymphocystis disease virus genotype I has never been reported in this species or in the region. Regardless of whether it has maintained a previously undetected state of endemicity in Egypt or was introduced through importation or contamination of ship ballast water, the findings of this study add to existing knowledge about the lymphocystis disease’s ecology, and lymphocystis disease virus genotypes and their host range.

Abstract

Marine and brackish water aquacultures are rapidly expanding in the Mediterranean basin. In this context, Egypt recently received a shipment of a 1.5 million juvenile gilthead seabream (Sparus aurata L.) from European Mediterranean facility. Within a few weeks of their arrival, 95% of the imported fish developed nodules on their skin and fins that lasted for several months. This study was undertaken to describe the clinical disease course, to identify the causative agent, and to investigate its origin. Preliminary diagnosis based on gross lesions and postmortem examination suggested lymphocystis disease (LCD), caused by the lymphocystis disease virus (LCDV; genus Lymphocystivirus, family Iridoviridae). Histopathological and ultrastructural features were typical of LCDV infections. PCR followed by sequencing and phylogenetic analysis of a 306-bp fragment of the major capsid protein (MCP) gene demonstrated the presence of LCDV genotype I, originally associated with LCD in Northern European countries, with 99.7% and 100% nucleotide and deduced amino acid identity values, respectively. LCDV genotype I has neither been reported in this species nor in the region. Regardless of the source of infection, findings of this study add to existing knowledge about the ecology of LCDV genotype I and its host range.

Evaluation of Gilthead Seabream (Sparus aurata) Immune Response after LCDV-Sa DNA Vaccination

Simple Summary

Lymphocystis disease is the main viral pathology in gilthead seabream aquaculture. Currently, there are no treatments or vaccines to control this disease, thus our main goal was to construct a DNA vaccine that can be used in the future to stop the spread of this pathology in sea farms. The vaccine consisted of a plasmid DNA that contains a known viral gene. Once it was established that the vaccine drives the expression of the antigenic viral protein in fish, vaccination experiments were conducted to determine if the vaccinated fish become protected against the viral infection. In addition, the immune response triggered by the vaccine was also evaluated in order to understand the mechanisms underlying such protection. The obtained results showed that in vaccinated fish an activation of several genes relating to both the inflammatory process and the mucosal immunity were produced, as well as specific anti-viral antibodies. Although limited, our results deserve further investigation to assess the efficacy of the vaccine in bigger fish populations and to confirm the mode of action of the vaccine.

Abstract

Lymphocystis disease is the main viral pathology reported in gilthead seabream. Its etiological agent is Lymphocystis disease virus 3 (LCDV-Sa), genus Lymphocystivirus, family Iridoviridae. There are no effective treatments or vaccines for LCDV control, thus the main aim of this study was to develop a DNA vaccine, and to evaluate both the protection conferred against LCDV-Sa infection and the immune response in vaccinated fish. The vaccine was constructed by cloning the mcp gene (ORF LCDVSa062R) into pcDNA3.1/NT-GFP-TOPO. Two independent vaccination trials were conducted. In the first one, 5–7 g fish were intramuscularly injected with the vaccine (pcDNA-MCP) or the empty-plasmid, and the distribution and expression of the vaccine was investigated. Furthermore, vaccinated fish were challenged with LCDV-Sa in order to access the protective capacity of the vaccine. In the second trial, 70–100 g fish were vaccinated as specified, and the immune response was evaluated analyzing the expression of 23 immune-related genes and the production of specific antibodies. The results showed that the vaccine triggers an immune response characterized by the overexpression of genes relating to the inflammatory process, but not the innate antiviral immunity relating to type I IFN (interferon), and also induces the production of specific neutralizing antibodies, which could explain the protection against LCDV-Sa in vaccinated fish.

Nano-scale applications in aquaculture: Opportunities for improved production and disease control

Aquaculture is the fastest growing food-production sector and is vital to food security, habitat restoration and endangered species conservation. One of the continued challenges to the industry is our ability to manage aquatic disease agents that can rapidly decimate operations and are a constant threat to sustainability. Such threats also evolve as microbes acquire resistance and/or new pathogens emerge. The advent of nanotechnology has transformed our approach to fisheries disease management with advances in water disinfection, food conversion, fish health and management systems. In this review, several nano-enabled technology successes will be discussed as they relate to the challenges associated with disease management in the aquaculture sector, with a particular focus on fishes. Future perspectives on how nanotechnology can offer functional approaches for improving disinfection and innovating at the practical space of early warning systems will be discussed. Finally, the importance of "safety by design" approaches to the development of novel commercial nano-enabled products will be emphasized.

Microalgal extracts induce larval programming and modify growth and the immune response to bioactive treatments and LCDV in Senegalese sole post-larvae

Immunostimulants are key molecules in aquaculture since they heighten defensive responses and protection against pathogens. The present study investigated the treatment of Senegalese sole larvae with a whole-cell crude extract of the microalgae Nannochloropsis gaditana (Nanno) and programming of growth and the immune system. Larvae at hatch were treated with the Nanno extracts for 2 h and thereafter were cultivated for 32 days post-hatch (dph) in parallel with an untreated control group (CN). Dry weight and length at 21 days post-hatch (dph) were higher in post-larvae of the Nanno than CN group. These differences in weight were later confirmed at 32 dph. To evaluate changes in the immune response associated with Nanno-programming treatments, the Nanno and CN post-larvae were supplied with two bioactive compounds yeast β-glucan (Y) and a microalga extract from the diatom Phaeodactylum tricornutum (MAe). The bioactive treatments were administrated to the treatment groups through the live prey (artemia metanauplii, 200 artemia mL^-1) enriched for 30 min with MAe or Y (at 2 mg mL^-1 SW) or untreated prey in the case of the negative control (SW). The effect of the treatments was assessed by monitoring gene expression, enzyme activity and mortality over 48 h. The post-larvae sole supplied with the bioactive compounds Y and MAe had increased mortality at 48 h compared to the SW group. Moreover, mortality was higher in Nanno-programmed than CN post-larvae. Lysozyme and total anti-protease enzymatic activities at 6 and 24 h after the start of the trial were significantly higher in the Nanno and MAe supplied post-larvae compared to their corresponding control (CN and SW, respectively). Immune gene transcripts revealed that il1b, cxc10 and mx mRNAs were significantly different between Nanno and CN post-larvae at 6 and 24 h. Moreover, the expression of il1b, tnfa, cxc10, irf3, irf7 and mx was modified by bioactive treatments but with temporal differences. At 48 h after bioactive treatments, Y and SW post-larvae were challenged with the lymphocystis disease virus (LCDV). No difference existed in viral copy number between programming or bioactive treatment groups at 3, 6 and 24 h after LCDV challenge although the total number of copies reduced with time. Gene expression profiles in the LCDV-challenged group indicated that post-larvae triggered a wide defensive response compared to SWC 24 h after challenge, which was modulated by programming and bioactive compound treatments. Cluster analysis of expressed genes separated the SW and Y groups indicating long-lasting effects of yeast β-glucan treatment in larvae. A noteworthy interaction between Nanno-programming and Y-treatment on the regulation of antiviral genes was observed. Overall, the data demonstrate the capacity of microalgal crude extracts to modify sole larval plasticity with long-term effects on larval growth and the immune responses.

Heritability Estimates and Genetic Correlation for Growth Traits and LCDV Susceptibility in Gilthead Sea Bream (Sparus aurata)

The lymphocystis disease (LCD) is a viral infection with a high economic impact in gilthead sea bream aquaculture. In this study, genetic estimates associated with lymphocystis disease virus (LCDV) disease susceptibility and growth were determined in sea bream juveniles. Two fish batches (named batch 1 and batch 2) were built from mass spawning and reared under industrial conditions until disease outbreak. At the moment of the sampling (n = 500 specimens for each batch), all animals had the typical LCDV lesions in the skin. For phenotyping, animals were weighted and photographed for image analysis (surface covered and lesion intensity). LCDV DNA copies were quantified in the liver by qPCR. Batch 1 had a higher surface covered and lesion intensity than batch 2, and the body caudal region was the lowest affected region in both batches. The average LCDV DNA copies in liver were higher in the batch 1 than batch 2, and they were positively correlated with severity index (SI) categories (r2 = 0.90–0.94). The total number of families evaluated were 150 and 128 for batch 1 and batch 2, respectively, with a high bias in offspring contribution by family and broodstock. Heritabilities for weight and length were 0.18 and 0.14 in batch 1 and 0.06 and 0.05 in batch 2, respectively. Heritability for the number of viral DNA copies was low (<0.08) in both batches. Heritabilities for SI in binary scale were 0.32/0.33 and 0.21/0.24 (underlying liability/Bayesian approach) for batch 1 and batch 2, respectively. Genetic correlations were very high and positive when growth traits (weight and length) or disease traits (LCDV DNA copies and SI) were compared. In contrast, the genetic correlations between growth and disease traits were moderate–high and positive in the batch 1 but negative in batch 2. These results indicate the genetic selection for LCDV susceptibility and growth is feasible in sea bream juveniles, although estimates are highly dependent on the age. The information provided is relevant to designing selective breeding programs in sea bream.

Immune gene expression in gilthead seabream (Sparus aurata) after Lymphocystis disease virus (LCDV-Sa) challenge resulting in asymptomatic infection

AIMS

To determine the immune gene expression response of gilthead seabream (Sparus aurata) that is experimentally infected with the lymphocystivirus LCDV-Sa.

METHODS AND RESULTS

Viral DNA and transcripts were detected by qPCR in all samples from fish injected with LCDV-Sa, demonstrating that the virus establish a systemic and asymptomatic infection. The expression of 23 immune-related genes was also analysed by RT-qPCR in the head kidney (HK) and intestine at several times post-infection (dpi). In HK, the expression of five type I interferon (IFN)-related genes (ifn, irf3, mx2, mx3 and isg15), il10 and ck10 was upregulated at 1-3 dpi, while genes related to the inflammation process (tnfα, il1ß, il6, casp1) were not differentially expressed or even downregulated. The expression profile in the intestine was different regarding type I INF-related genes. An upregulated c3 and ighm expression was observed in both HK and intestine at 3-8 dpi. Finally, the transcription of nccrp1 and mhcIIα was induced in HK, whereas tcrβ expression was downregulated in both organs.

CONCLUSIONS

LCDV-Sa seems to trigger an immune response in gilthead seabream characterized by a partial activation of type I IFN system and a lack of systemic inflammatory response which may be related to viral persistence.

SIGNIFICANCE AND IMPACT OF THE STUDY

The immune response observed in gilthead seabream infected by LCDV-Sa could be implicated in the establishment of an asymptomatic persistent infection.

Feed and immersion challenges with lymphocystis disease virus (LCDV) reveals specific mechanisms for horizontal transmission and immune response in senegalese sole post-larvae

The horizontal transmission of lymphocystis disease virus (LCDV) through contaminated water and feed (using artemia as vehicle) and the associated immune gene expression profiles in Senegalese sole post-larvae were investigated. All specimens analyzed were positive for LCDV DNA detection at 1-day post-challenge (1 dpc) with the highest viral levels in specimens infected through the immersion route. However, the percentage of LCDV-positive animals and number of viral DNA copies dropped progressively at 2 and 7 dpc. The histological analysis identified structural changes in the skin, muscle and gills of sole post-larvae LCDV-challenged by immersion. In situ hybridization confirmed a wide distribution of LCDV in the skin, gut, surrounding vessels in trunk muscle and head kidney in the immersion route, while the signals were restricted to the liver and lamina propria in the feeding treatment. Expression analysis using a set of 22 genes related to innate immune defense system demonstrated clear differences in the time-course response to LCDV as function of the infection route. Most antiviral defense genes, the proinflammatory cytokines, the complement c3, g-type lysozyme and T-cell markers cd4 and cd8a were rapidly induced in the feeding-infected post-larvae, and they were remained activated at 2 dpc. In contrast, in the immersion-infected post-larvae the induction of most defensive genes was delayed, with a low intensity at 2 dpc. All these data demonstrate that LCDV can horizontally infect Senegalese sole post-larvae through the water or feed although with different patterns of histopathological disorders, virus distribution and route-specific expression profiles.

Artemia spp., a Susceptible Host and Vector for Lymphocystis Disease Virus

Different developmental stages of Artemia spp. (metanauplii, juveniles and adults) were bath-challenged with two isolates of the Lymphocystis disease virus (LCDV), namely, LCDV SA25 (belonging to the species Lymphocystis disease virus 3) and ATCC VR-342 (an unclassified member of the genus Lymphocystivirus). Viral quantification and gene expression were analyzed by qPCR at different times post-inoculation (pi). In addition, infectious titres were determined at 8 dpi by integrated cell culture (ICC)-RT-PCR, an assay that detects viral mRNA in inoculated cell cultures. In LCDV-challenged Artemia, the viral load increased by 2–3 orders of magnitude (depending on developmental stage and viral isolate) during the first 8–12 dpi, with viral titres up to 2.3 × 10² Most Probable Number of Infectious Units (MPNIU)/mg. Viral transcripts were detected in the infected Artemia, relative expression values showed a similar temporal evolution in the different experimental groups. Moreover, gilthead seabream (Sparus aurata) fingerlings were challenged by feeding on LCDV-infected metanauplii. Although no Lymphocystis symptoms were observed in the fish, the number of viral DNA copies was significantly higher at the end of the experimental trial and major capsid protein (mcp) gene expression was consistently detected. The results obtained support that LCDV infects Artemia spp., establishing an asymptomatic productive infection at least under the experimental conditions tested, and that the infected metanauplii are a vector for LCDV transmission to gilthead seabream.

Development and validation of a reverse transcription quantitative polymerase chain reaction for tilapia lake virus detection in clinical samples and experimentally challenged fish

Tilapia lake virus (TiLV) is an emerging pathogen associated with high mortalities of wild and farm-raised tilapia in different countries. In this study, a SYBR green-based reverse transcription quantitative polymerase chain reaction (RT-qPCR) assay targeting segment three of the virus was developed to detect and quantify TiLV in clinical samples and experimentally challenged fish. All 30 field samples with clinical signs and history consistent with TiLV infection were positive for TiLV as detected by the developed RT-qPCR method. The RT-qPCR technique provided 100 and 10,000 times more sensitive for virus detection than those offered by the RT-PCR and virus isolation in cell culture methods, respectively. The detection limit of the RT-qPCR method was as low as two viral copies/μl. Moreover, the RT-qPCR technique could be applied for TiLV detection in various fish tissues including gills, liver, brain, heart, anterior kidney and spleen. Significantly, this study delivered an accurate and reliable method for rapid detection of TiLV viruses that facilitates active surveillance programme and disease containment.

Gene expression profiles associated with lymphocystis disease virus (LCDV) in experimentally infected Senegalese sole (Solea senegalensis)

In the present study, the pathogenesis of lymphocystis disease virus (LCDV) and the immune gene expression patterns associated with this viral infection were determined in the flatfish Senegalese sole. The results indicate that LCDV spreads rapidly from the peritoneal cavity through the bloodstream to reach target organs such as kidney, gut, liver, and skin/fin. The viral load was highest in kidney and reduced progressively thorough the experiment in spite of the viral major capsid protein gene was transcribed. The LCDV injection activated a similar set of differentially expressed transcripts in kidney and intestine although with some differences in the intensity and time-course response. This set included antiviral-related transcripts (including the mx and interferon-related factors irf1, irf2, irf3, irf7, irf8, irf9, irf10), cytokines (il1b, il6, il8, il12 and tnfa) and their receptors (il1r, il8r, il10r, il15ra, il17r), chemokines (CXC-type, CC-type and IL-8), prostaglandins (cox-2), g-type lysozymes, hepcidin, complement fractions (c2, c4-1 and c4-2) and the antigen differentiation factors cd4, cd8a, and cd8b. The expression profile observed indicated that the host triggered a systemic defensive response including inflammation able to cope with the viral challenge.

Target organs for lymphocystis disease virus replication in gilthead seabream (Sparus aurata)

The lymphocystis disease (LCD), the main viral pathology described in cultured gilthead seabream (Sparus aurata), is a self-limiting condition characterized by the appearance of hypertrophied fibroblasts (named lymphocysts) in the connective tissue of fish, primarily in the skin and fins. The causative agent of the disease is the Lymphocystis disease virus (LCDV), a member of the Iridoviridae family. In the present study, LCDV genome and transcripts were detected by real-time PCR in caudal fin, as well as in several internal organs, such as intestine, liver, spleen, kidney and brain, from asymptomatic, diseased and recovered gilthead seabream juveniles. These results indicate that the LCDV has a broad range tissue tropism, and can establish a systemic infection, even in subclinically infected fish. As showed by in situ hybridization, the permissive cells for LCDV infection seem to be fibroblasts, hepatocytes and cells of the mononuclear phagocyte system. Histopathological alterations associated with LCD were observed in all the organs analysed, including necrotic changes in liver and kidney, inflammatory response in the intestine submucosa or brain haemorrhage, although lymphocysts were only detected in the dermis of the caudal fin. Nevertheless, these histological changes were reverted in recovered animals.

Rapid and Sensitive Detection of Lymphocystis Disease Virus Genotype VII by Loop-Mediated Isothermal Amplification

Lymphocystis disease virus (LCDV) infections have been described in gilthead seabream (Sparus aurata L.) and Senegalese sole (Solea senegalensis, Kaup), two of the most important marine fish species in the Mediterranean aquaculture. In this study, a rapid, specific, and sensitive detection method for LCDV genotype VII based on loop-mediated isothermal amplification (LAMP) was developed. The LAMP assay, performed using an apparatus with real-time amplification monitoring, was able to specifically detect LCDV genotype VII from clinically positive samples in less than 12 min. In addition, the assay allowed the detection of LCDV in all asymptomatic carrier fish analysed, identified by qPCR, showing an analytical sensitivity of ten copies of viral DNA per reaction. The LCDV LAMP assay has proven to be a promising diagnostic method that can be used easily in fish farms to detect the presence and spread of this iridovirus.

Evaluation of an integrated cell culture RT-PCR assay to detect and quantify infectious lymphocystis disease virus

The lymphocystis disease virus (LCDV), a member of the Iridoviridae family, infects a wide range of fish species including gilthead seabream (Sparus aurata L.), the most important species cultured in the Mediterranean. LCDV is difficult to propagate in cell culture and does not produce clear and consistent cytopathic effects (CPE), especially in samples collected from subclinically infected fish. An integrated cell culture reverse transcription-polymerase chain reaction (ICC-RT-PCR) assay, followed by dot-blot hybridization of the RT-PCR products, was developed to improve the detection of infectious LCDV. The sensitivity of the ICC-RT-PCR assay, which can be performed in 7 d, was at least 100-fold higher than viral diagnosis obtained by CPE development. The developed assay thus allows the determination of infectious titres in samples with low viral loads, including those from asymptomatic carrier fish, in which no CPE was recorded after a 14-d incubation period. The ICC-RT-PCR assay enables rapid, specific and sensitive detection and quantification of infectious LCDV, and may be a valuable tool in the study of aspects of LCDV infection including transmission or epizootiology.

Concurrence of Iridovirus, Polyomavirus, and a Unique Member of a New Group of Fish Papillomaviruses in Lymphocystis Disease-Affected Gilthead Sea Bream

Lymphocystis disease is a geographically widespread disease affecting more than 150 different species of marine and freshwater fish. The disease, provoked by the iridovirus lymphocystis disease virus (LCDV), is characterized by the appearance of papillomalike lesions on the skin of affected animals that usually self-resolve over time. Development of the disease is usually associated with several environmental factors and, more frequently, with stress conditions provoked by the intensive culture conditions present in fish farms. In gilthead sea bream (Sparus aurata), an economically important cultured fish species in the Mediterranean area, a distinct LCDV has been identified but not yet completely characterized. We have used direct sequencing of the virome of lymphocystis lesions from affected S. aurata fish to obtain the complete genome of a new LCDV-Sa species that is the largest vertebrate iridovirus sequenced to date. Importantly, this approach allowed us to assemble the full-length circular genome sequence of two previously unknown viruses belonging to the papillomaviruses and polyomaviruses, termed Sparus aurata papillomavirus 1 (SaPV1) and Sparus aurata polyomavirus 1 (SaPyV1), respectively. Epidemiological surveys showed that lymphocystis disease was frequently associated with the concurrent appearance of one or both of the new viruses. SaPV1 has unique characteristics, such as an intron within the L1 gene, and as the first member of the Papillomaviridae family described in fish, provides evidence for a more ancient origin of this family than previously thought.

IMPORTANCE

Lymphocystis disease affects marine and freshwater fish species worldwide. It is characterized by the appearance of papillomalike lesions on the skin that contain heavily enlarged cells (lymphocysts). The causative agent is the lymphocystis disease virus (LCDV), a large icosahedral virus of the family Iridoviridae In the Mediterranean area, the gilthead sea bream (Sparus aurata), an important farmed fish, is frequently affected. Using next-generation sequencing, we have identified within S. aurata lymphocystis lesions the concurrent presence of an additional LCDV species (LCDV-Sa) as well as two novel viruses. These are members of polyomavirus and papillomavirus families, and here we report them to be frequently associated with the presence of lymphocysts in affected fish. Because papillomaviruses have not been described in fish before, these findings support a more ancient origin of this virus family than previously thought and evolutionary implications are discussed.