Transcriptional profile of red seabream iridovirus in a fish model as revealed by viral DNA microarrays

Partial validation of a TaqMan quantitative polymerase chain reaction for the detection of the three genotypes of Infectious spleen and kidney necrosis virus

Megalocytiviruses (MCVs) are double-stranded DNA viruses known to infect important freshwater and marine fish species in the aquaculture, food, and ornamental fish industries worldwide. Infectious spleen and kidney necrosis virus (ISKNV) is the type species within the genus Megalocytivirus that causes red seabream iridoviral disease (RSIVD) which is a reportable disease to the World Animal Health Organization (WOAH). To better control the transboundary spread of this virus and support WOAH reporting requirements, we developed and partially validated a TaqMan real-time qPCR assay (ISKNV104R) to detect all three genotypes of ISKNV, including the two genotypes that cause RSIVD. Parameters averaged across 48 experiments used a 10-fold dilution series of linearized plasmid DNA (107-101 copies), carrying a fragment of the three-spot gourami iridovirus (TSGIV) hypothetical protein revealed that the assay was linear over 7 orders of magnitude (107-101), a mean efficiency of 99.97 ± 2.92%, a mean correlation coefficient of 1.000 ± 0.001, and a limit of detection (analytical sensitivity) of ≤10 copies of TSGIV DNA. The diagnostic sensitivity and specificity for the ISKNV104R qPCR assay was evaluated and compared to other published assays using a panel of 397 samples from 21 source populations with different prevalence of ISKNV infection (0-100%). The diagnostic sensitivity and specificity for the ISKNV104R qPCR assay was 91.99% (87.28-95.6; 95% CI) and 89.8% (83.53-94.84). The latent class analysis showed that the ISKNV104R qPCR assay had similar diagnostic sensitivities and specificities with overlapping confidence limits compared to a second TaqMan qPCR assay and a SYBR green assay. This newly developed TaqMan assay represents a partially validated qPCR assay for the detection of the three genotypes of the species ISKNV. The ISKNV104R qPCR assay once fully validated, will serve as an improved diagnostic tool that can be used for ISKNV surveillance efforts and diagnosis in subclinical fish to prevent further spread of MCVs throughout the aquaculture and ornamental fish industries.

Analysis of the transcriptomic profiles of Mandarin fish (Siniperca chuatsi) infected with red sea bream iridovirus (RSIV)

Red sea bream iridovirus (RSIV) belongs to the family Iridoviridae, genus Megalocytivirus, which could widely infect marine fish, causing diseases and huge economic losses. Now it has been reported that RSIV was also detected in diseased mandarin fish. Transmission electron microscopy and immunohistochemistry showed that spleen was the main target organ in mandarin fish infected with RSIV. To investigate the immune response mechanism of mandarin fish to RSIV infection, transcriptomics of RSIV-infected mandarin fish was analyzed. A total of 53,040 unigenes were obtained, and there were 21,576 and 17,904 unigenes had significant hit the Nr and SwissProt databases, respectively. In RSIV-infected and non-infected spleen tissues, there were 309 differentially expressed genes (DEGs), including 100 up-regulated genes and 209 down-regulated genes. Gene Ontology database (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways analysis were performed to reveal the function information and give a better understanding of the signal transduction pathways of DEGs. Further analysis of the cytokine-cytokine receptor interactions pathway exhibited that the expression of cytokines was widely activated after viral infection. In addition, ten DEGs were randomly selected and verified by quantitative real-time PCR, which revealed a similar expression tendency as the high-throughput sequencing data. These findings present valuable information that will benefit for better understanding of RSIV infection in mandarin fish.

Conserved motifs in the invertebrate iridescent virus 6 (IIV6) genome regulate virus transcription

Invertebrate iridescent virus 6 (IIV6) is the type species of the Iridovirus genus in the Betairidovirinae subfamily of the Iridoviridae family. Transcription of the 215 predicted IIV6 genes is temporally regulated, dividing the genes into three kinetic classes: immediate-early (IE), delayed-early (DE), and late (L). So far, the transcriptional class has been determined for a selection of virion protein genes and only for three genes the potential promoter regions have been analyzed in detail. In this study, we investigated the transcriptional class of all IIV6 genes that had not been classified until now. RT-PCR analysis of total RNA isolated from virus-infected insect cells in the presence or absence of protein and DNA synthesis inhibitors, placed 113, 23 and 22 of the newly analyzed viral ORFs into the IE, DE and L gene classes, respectively. Afterwards, in silico analysis was performed to the upstream regions (200 bp) of all viral ORFs using the MEME Suite Software. The AA(A/T)(T/A)TG(A/G)A and (T/A/C)(T/G/C)T(T/A)ATGG motifs were identified in the upstream region of IE and DE genes, respectively. These motifs were validated by luciferase reporter assays as crucial sequences for promoter activity. For the L genes two conserved motifs were identified for all analyzed genes: (T/G)(C/T)(A/C)A(T/G/C)(T/C)T(T/C) and (C/G/T)(G/A/C)(T/A)(T/G) (G/T)(T/C). However, the presence of these two motifs did not influence promoter activity. Conversely, the presence of these two sequences upstream of the reporter decreased its expression. Single nucleotide mutations in the highly conserved nucleotides at the end of the second motif (TTGT) showed that this motif acted as a repressor sequence for late genes in the IIV6 genome. Next, upstream sequences of IIV6 L genes from which we removed this second motif in silico, were re-analyzed for the presence of potential conserved promoter sequences. Two additional motifs were identified in this way for L genes: (T/A)(A/T)(A/T/G)(A/T)(T/C)(A/G)(A/C)(A/C) and (C/G)(T/C)(T/A/C)C(A/T)(A/T)T(T/G) (T/G)(T/G/A). Independent mutations in either motif caused a severe decrease in luciferase expression. Information on temporal classes and upstream regulatory sequences will contribute to our understanding of the transcriptional mechanisms in IIV6.

Transcriptional programs of infectious spleen and kidney necrosis virus (ISKNV) in vitro and in vivo

Infectious spleen and kidney necrosis virus (ISKNV), causing serious infectious diseases to marine and freshwater fishes, is the type species of the genus Megalocytivirus, family Iridoviridae. In this study, the transcriptional programs of ISKNV in vitro (MFF-1 cells) and in vivo (spleens from mandarin fish) were investigated using real-time PCR. Transcription of all the putative open reading frames (ORFs) of ISKNV was verified. The temporal expression patterns of ISKNV ORFs in vitro and in vivo, including peak expression times (PETs) and relative maximal expression levels, were determined and compared. The K-means clustering with Spearman rank correlation was generated in heat maps constructed based on ISKNV ORF expression profiles in vivo and in vitro. The current study may provide a global picture of ISKNV infection at the transcription level and help better understand the molecular pathogenic mechanism of megalocytiviruses.

Transcriptomic analysis of the host response to an iridovirus infection in Chinese giant salamander, Andrias davidianus

The emergence of an infectious viral disease caused by the Chinese giant salamander iridovirus (GSIV) has led to substantial economic losses. However, no more molecular information is available for the understanding of the mechanisms associated with virus–host interaction. In this study, de novo sequencing was used to obtain abundant high-quality ESTs and investigate differentially-expressed genes in the spleen of Chinese giant salamanders that were either infected or mock infected with GSIV. Comparative expression analysis indicated that 293 genes were down-regulated and 220 genes were up-regulated. Further enrichment analysis showed that the most enriched pathway is “complement and coagulation cascades”, and significantly enriched diseases include “inherited thrombophilia”, “immune system diseases”, “primary immunodeficiency”, “complement regulatory protein defects”, and “disorders of nucleotide excision repair”. Additionally, 30 678 simple sequence repeats (SSRs) from all spleen samples, 26 355 single nucleotide polymorphisms (SNPs) from the spleens of uninfected animals and 36 070 SNPs from the spleens of infected animals were detected. The large amount of variation was specific for the Chinese giant salamanders that were infected with GSIV. The results reported herein provided significant and new EST information that could contribute greatly in investigations into the molecular functions of immune genes in the Chinese giant salamander.

Rana grylio virus as a vector for foreign gene expression in fish cells

In the present study, Rana grylio virus (RGV, an iridovirus) thymidine kinase (TK) gene and viral envelope protein 53R gene were chosen as targets for foreign gene insertion. ΔTK-RGV and Δ53R-RGV, two recombinant RGV, expressing enhanced green fluorescence protein (EGFP) were constructed and analyzed in Epithelioma papulosum cyprinid (EPC) cells. The EGFP gene which fused to the virus major capsid protein (MCP) promoter p50 was inserted into TK and 53R gene loci of RGV, respectively. Cells infected with these two recombinant viruses not only displayed plaques, but also emitted strong green fluorescence under fluorescence microscope, providing a simple method for selection and purification of recombinant viruses. ΔTK-RGV was purified by seven successive rounds of plaque isolation and could be stably propagated in EPC cells. All of the plaques produced by the purified recombinant virus emitted green fluorescence. However, Δ53R-RGV was hard to be purified even through twenty rounds of plaque isolation. The purified recombinant virus ΔTK-RGV was verified by PCR analysis and Western blotting. These results showed EGFP was expressed in ΔTK-RGV infected cells. Furthermore, one-step growth curves and electron microscopy revealed that infection with recombinant ΔTK-RGV and wild-type RGV are similar. Therefore, RGV was demonstrated could be as a viral vector for foreign gene expression in fish cells.

Microarray technology is an effective tool for identifying genes related to the aquacultural improvement of Japanese flounder, Paralichthys olivaceus

Molecular techniques are now essential for discovering new functional genes for the improvement of cultured marine organisms. Such techniques can be used to obtain expressed sequence tags (EST) and, through the use of microarrays, the entire transcriptome. For example, microarrays can be used to reveal biomarkers of health conditions, nutrient changes and immune response in fish and shellfish. EST-based microarray chips were constructed for cultured fish and shellfish species including Japanese flounder (Paralichthys olivaceus). Using the flounder microarray chip, the efficacy of two DNA vaccines derived from pathogenic viruses [hirame rhabdovirus (HRV) and viral hemorrhagic septicemia virus (VHSV)] was evaluated through gene expression profiles. The results suggest that both DNA vaccines were effective in protecting the flounder from HIRRV and VHSV. The flounder microarray was also used to compare gene expression patterns in fish that are susceptible and resistant to Edwardsiella tarda. At 3days post infection, the gene expression patterns between the two groups were dramatically changed. Thus, microarray analysis is a very powerful tool to understand gene expression profiles in fish.

Stress transcriptomics in fish: a role for genomic cortisol signaling

The physiological responses to stressors, including hormonal profiles and associated tissue responsiveness have been extensively studied in teleosts, but the molecular mechanisms associated with this adaptive response are not well understood. The advent of cDNA microarray technology has transformed the field of functional genomics by revealing global gene expression changes in response to stressor exposures even in non-mammalian vertebrates, including fish. A unifying response in studies related to stressor exposure is activation of the hypothalamus-pituitary-interrenal (HPI) axis in fish, leading to cortisol release into the circulation. Here we will discuss the implications of some of the gene expression changes observed in response to acute stress in fish, while highlighting a role for cortisol in this adaptive stress response. As liver is a key organ for metabolic adjustments to stressors and also is a major target for cortisol action, the genomic studies pertaining to stress and glucocorticoid regulation have focused mainly on this tissue. The studies have identified several genes that are altered transiently after an acute stressor exposure in fish. A number of these stress-responsive genes were also modulated by glucocorticoid receptor activation, suggesting that elevation in cortisol levels during stressor exposure may be a key signal for target tissue molecular programming, essential for stress adaptation. The identification of regulatory gene networks that are stress activated, and modulated by cortisol, both in hepatic and extra-hepatic tissues, including gonads, brain, immune cells and gills, will provide a mechanistic framework to characterize the multifaceted role of cortisol during stress adaptation.

Transcriptome analysis of Frog virus 3, the type species of the genus Ranavirus, family Iridoviridae

Frog virus 3 is the best characterized species within the genus Ranavirus, family Iridoviridae. FV3's large ( approximately 105 kbp) dsDNA genome encodes 98 putative open reading frames (ORFs) that are expressed in a coordinated fashion leading to the sequential appearance of immediate early (IE), delayed early (DE) and late (L) viral transcripts. As a step toward elucidating molecular events in FV3 replication, we sought to identify the temporal class of viral messages. To accomplish this objective an oligonucleotide microarray containing 70-mer probes corresponding to each of the 98 FV3 ORFs was designed and used to examine viral gene expression. Viral transcription was initially monitored during the course of a productive replication cycle at 2, 4 and 9 h after infection. To confirm results of the time course assay, viral gene expression was also monitored in the presence of cycloheximide (CHX), which limits expression to only IE genes, and following infection with a temperature-sensitive (ts) mutant which at non-permissive temperatures is defective in viral DNA synthesis and blocked in late gene expression. Subsequently, microarray analyses were validated by RT-PCR and qRT-PCR. Using these approaches we identified 33 IE genes, 22 DE genes and 36 L viral genes. The temporal class of the 7 remaining genes could not be determined. Comparison of protein function with temporal class indicated that, in general, genes encoding putative regulatory factors, or proteins that played a part in nucleic acid metabolism and immune evasion, were classified as IE and DE genes, whereas those involved in DNA packaging and virion assembly were considered L genes. Information on temporal class will provide the basis for determining whether members of the same temporal class contain common upstream regulatory regions and perhaps allow us to identify virion-associated and virus-induced proteins that control viral gene expression.

Phenotypic diversity of infectious red sea bream iridovirus isolates from cultured fish in Japan

Megalocytivirus is causing economically serious mass mortality by infecting fish in and around the Pacific region of Asia. The recent emergence of many new iridoviruses has drawn attention to the marked taxonomic variation within this virus family. Most studies of these viruses have not included extensive study of these emergent species. We explored the emergence of red sea bream iridovirus (RSIV) on a fish farm in Japan, and we specifically endeavored to quantify genetic and phenotypic differences between RSIV isolates using in vitro and in vivo methods. The three isolates had identical major capsid protein sequences, and they were closely related to Korean RSIV isolates. In vitro studies revealed that the isolates differed in replication rate, which was determined by real-time quantitative PCR of viral genomes in infected cells and cell culture supernatant, and in cell viability, estimated by the MTT (3-[4,5-dimethylthiazol-2-yl]-2,5-diphenyltetrazolium bromide) assay for infected cells. In vivo studies showed that the isolates exhibit different virulence characteristics: infected red sea bream showed either acute death or subacute death according to infection with different isolates. Significant differences were seen in the antigenicity of isolates by a formalin-inactivated vaccine test. These results revealed that variant characteristics exist in the same phylogenetic location in emergent iridoviruses. We suggest that this strain variation would expand the host range in iridoviral epidemics.

Whole-genome transcriptional profiles of a novel marine fish iridovirus, Singapore grouper iridovirus (SGIV) in virus-infected grouper spleen cell cultures and in orange-spotted grouper, Epinephulus coioides

A DNA microarray containing all Singapore grouper iridovirus (SGIV) open reading frames (ORFs) was constructed to map the viral gene transcriptional profiles in virus-infected grouper spleen (GS) cells and in spleen tissues of virus-infected grouper. The results showed that viral genes started to be transcribed as early as 1 h postinfection (p.i.), and followed by a rapid increasing gene expression along with virus infection in cell cultures. The three temporal kinetic classes (15 immediate-early, 89 early and 53 late transcripts) were classified during an in vitro infection by their dependence on de novo protein synthesis and viral DNA replication inhibitors. In SGIV-infected grouper, Epinephulus coioides, most of the viral genes were expressed between 1 and 4 d p.i., and the number and expression levels started to decrease after 5 d p.i. These data were confirmed by real-time RT-PCR. This study provides an experimental basis for investigation of virus-host interactions and the development of control strategies against SGIV infection.