Rana grylio virus 43R encodes an envelope protein involved in virus entry.

Rana grylio virus (RGV), a member of genus Ranavirus in the family Iridoviridae, is a viral pathogen infecting aquatic animal. RGV 43R has homologues only in Ranavirus and contains a transmembrane (TM) domain, but its role in RGV infection is unknown. In this study, 43R was determined to be associated with virion membrane. The transcripts encoding 43R and the protein itself appeared late in RGV-infected EPC cells and its expression was blocked by viral DNA replication inhibitor, indicating that 43R is a late expressed protein. Subcellular localization showed that 43R-EGFP fusion protein distributed in cytoplasm of EPC cells and that TM domain is essential for its distribution in cytoplasm. 43R-EGFP fusion protein colocalized with viral factories in RGV-infected cells. A recombinant RGV deleting 43R (Δ43R-RGV) was constructed by homologous recombination to investigate its role in virus infection. Compared with wild type RGV, the ability of Δ43R-RGV to induce the cytopathic effect and its virus titers were significantly reduced. Furthermore, it is revealed that 43R deletion significantly inhibited viral entry but did not influence viral DNA replication by measuring and comparing the DNA levels of RGV and Δ43R-RGV in the infected cells at the early stage of infection. RGV neutralization with anti-43R serum reduced the virus titer. Therefore, these data showed that RGV 43R is a late gene that encodes an envelope protein involved in RGV entry.

Discovery and characterization of conserved and novel microRNAs from blunt snout bream (Megalobrama amblycephala) by deep sequencing.

MicroRNAs (miRNAs) are short, single stranded RNA molecules with approximately 22 nts in length, which regulate the stability and translation of messenger RNAs in several organisms. To increase the repertoire of miRNAs characterized in M. amblycephala, we used the deep sequencing technology to sequence a small RNA library using pooled RNA sample isolated from the 4 different tissues of M. amblycephala. A total of 309 conserved miRNAs that originated from 131 miRNA families were detected. 15 novel candidates miRNA were identified. Randomly selected 6 miRNAs were analyzed by stem-loop qRT-PCR and differential expression patterns were observed in 6 different tissues of M. amblycephala. Furthermore, the potential targets were predicted. GO analysis showed that most of the targets were involved in a broad range of physiological functions including fish growth, development, metabolism, stress responses and so on. Overall, our results significantly increased the number of novel miRNAs in M. amblycephala, which should be useful for further investigation into the role of miRNAs in regulating diverse biological processes.

Protective Immunity Induced by DNA Vaccination against Ranavirus Infection in Chinese Giant Salamander Andrias davidianus.

Andrias davidianus ranavirus (ADRV) is an emerging viral pathogen that causes severe systemic hemorrhagic disease in Chinese giant salamanders. There is an urgent need for developing an effective vaccine against this fatal disease. In this study, DNA vaccines containing the ADRV 2L gene (pcDNA-2L) and the 58L gene (pcDNA-58L) were respectively constructed, and their immune protective effects were evaluated in Chinese giant salamanders. In vitro and in vivo expression of the vaccine plasmids were confirmed in transfected cells and muscle tissues of vaccinated Chinese giant salamanders by using immunoblot analysis or RT-PCR. Following ADRV challenge, the Chinese giant salamanders vaccinated with pcDNA-2L showed a relative percent survival (RPS) of 66.7%, which was significant higher than that in Chinese giant salamanders immunized with pcDNA-58L (RPS of 3.3%). Moreover, the specific antibody against ADRV was detected in Chinese giant salamanders vaccinated with pcDNA-2L at 14 and 21 days post-vaccination by indirect enzyme-linked immunosorbent assay (ELISA). Transcriptional analysis revealed that the expression levels of immune-related genes including type I interferon (IFN), myxovirus resistance (Mx), major histocompatibility complex class IA (MHCIA), and immunoglobulin M (IgM) were strongly up-regulated after vaccination with pcDNA-2L. Furthermore, vaccination with pcDNA-2L significantly suppressed the virus replication, which was seen by a low viral load in the spleen of Chinese giant salamander survivals after ADRV challenge. These results indicated that pcDNA-2L could induce a significant innate immune response and an adaptive immune response involving both humoral and cell-mediated immunity that conferred effective protection against ADRV infection, and might be a potential vaccine candidate for controlling ADRV disease in Chinese giant salamanders.

Transcriptome analysis of the Chinese giant salamander (Andrias davidianus) using RNA-sequencing.

The Chinese giant salamander (Andrias davidianus) is an economically important animal on academic value. However, the genomic information of this species has been less studied. In our study, the transcripts of A. davidianus were obtained by RNA-seq to conduct a transcriptomic analysis. In total 132,912 unigenes were generated with an average length of 690 bp and N50 of 1263 bp by de novo assembly using Trinity software. Using a sequence similarity search against the nine public databases (CDD, KOG, NR, NT, PFAM, Swiss-prot, TrEMBL, GO and KEGG databases), a total of 24,049, 18,406, 36,711, 15,858, 20,500, 27,515, 36,705, 28,879 and 10,958 unigenes were annotated in databases, respectively. Of these, 6323 unigenes were annotated in all database and 39,672 unigenes were annotated in at least one database. Blasted with KEGG pathway, 10,958 unigenes were annotated, and it was divided into 343 categories according to different pathways. In addition, we also identified 29,790 SSRs. This study provided a valuable resource for understanding transcriptomic information of A. davidianus and laid a foundation for further research on functional gene cloning, genomics, genetic diversity analysis and molecular marker exploitation in A. davidianus.

Characterization of conserved and novel miRNAs using deep sequencing and prediction of miRNA targets in Crucian carp (Carassius auratus).

MicroRNAs (miRNAs) are a class of endogenous small non-coding RNAs of -22 nucleotides that can base pair with their target mRNAs, which represses their translation or induces their degradation in various biological processes. However, little is known about identification of miRNAs and their target genes in C. auratus. In the present study, a small RNA library from pooled tissue of C. auratus was constructed and sequenced using the deep sequencing. A total of 320 conserved miRNAs (belonging to 105 families) as well as 11 potentially novel miRNAs were identified. Stem-loop qRT-PCR analysis confirmed that both conserved and novel miRNAs were expressed in C. auratus, and some of them were preferentially expressed in certain tissues. Subsequently, a total of 1668 potential target genes were predicted for these identified miRNAs and GO analysis showed that most of the targets were involved in lots of physiological actions. This study represents a first large-scale identification and characterization of C. auratus miRNAs and their potential target genes. Taken together, our results add new information to existing data on C. auratus miRNAs and should be useful for investigating the biological functions of miRNAs in fishes and other aquatic species.

Identification and characterization of the Chinese giant salamander (Andrias davidianus) miRNAs by deep sequencing and predication of their targets.

MicroRNAs (miRNAs) are non-coding small RNA of approximately 22 nt in length. They are considered as key regulators for gene expression at the post-transcriptional level and play an important role in the regulation of many fundamental biological processes. Andrias davidianus, as one of the largest amphibian on earth, may represent the transitional type of animal from aquatic to terrestrial life, but so far, no miRNA has been identified in the species. In this study, Illumina deep sequencing was used for high-throughput analysis of miRNAs in a pooled small RNA library isolated from eight tissues sample of A. davidianus including the spleen, liver, muscle, kidney, skin, testis, gut, and heart. In total, 6,213,146 reads were obtained, 5,962,175 of which were related to 143 unique miRNAs, including 140 conserved and three novel A. davidianus-specific miRNAs. Among them, seven conserved miRNAs and one novel miRNA were selected to validate their expression pattern by stem-loop qRT-PCR. Moreover, 4700 potential target genes were predicted for 143 A. davidianus miRNAs; GO enrichment analysis and KEGG pathway showed that most of the targets are involved in diverse biological processes including ubiquitin-mediated proteolysis, FOXO signaling pathway, tight junction, and spliceosome. This study provides the first large-scale identification and characterization of A. davidianus miRNAs, and predicted their potential target genes; it will lay a valuable foundation for future understanding the role of these miRNAs on regulating diverse biological processes.

Identification and characterization of known and novel microRNAs in three tissues of Chinese giant salamander base on deep sequencing approach.

MicroRNAs (miRNAs) are a class of non-coding small RNAs that play important roles in the regulation of gene expression. Although animal miRNAs have been extensively studied in model systems, less is known in other animal with limited genome sequence data, including Chinese giant salamander (Andrias davidianus). The identification of the full repertoire of miRNAs expressed in the liver, spleen and muscle of A. davidianus would significantly increase our understanding for physiological function of A. davidianus, in this ancient and endangered urodele amphibian. In this study, three independent small RNA libraries were constructed from the liver, spleen and muscle of A. davidianus. The libraries were subjected to high-throughput sequencing by using the Illumina deep sequencing. As a result, a total of 12,831,239, 13,592,195 and 9,887,531 raw reads representing 2,240,771, 1,363,266 and 1,964,252 clean reads per library were obtained separately. Through bioinformatics analysis, we identified total of 553 known miRNAs and 44 putative novel miRNAs in our small RNA dataset from liver, spleen and muscle tissues. Five known miRNAs (gga-miR-10a-5p, pma-miR-29d-5p, aca-miR-338-3p, hsa-miR-455-3p and ssa-miR-2184-5p_R-1) and three novel miRNAs (PC-5p-891_1763, PC-5p-32538_50 and PC-3p-33645_48) showed different expression in eight different tissues as revealed by stem-loop qPCR analysis. This study characterized the miRNA of A. davidianus for the first time, which provides an opportunity for further understanding of miRNA regulation function in A. davidianus ranked as living fossils.

Sequencing and de novo transcriptome assembly of the Chinese giant salamander (Andrias davidianus).

Next-generation technologies for determination of genomics and transcriptomics composition have a wide range of applications. Andrias davidianus, has become an endangered amphibian species of salamander endemic in China. However, there is a lack of the molecular information. In this study, we obtained the RNA-Seq data from a pool of A. davidianus tissue including spleen, liver, muscle, kidney, skin, testis, gut and heart using Illumina HiSeq 2500 platform. A total of 15,398,997,600 bp were obtained, corresponding to 102,659,984 raw reads. A total of 102,659,984 reads were filtered after removing low-quality reads and trimming the adapter sequences. The Trinity program was used to de novo assemble 132,912 unigenes with an average length of 690 bp and N50 of 1263 bp. Unigenes were annotated through number of databases. These transcriptomic data of A. davidianus should open the door to molecular evolution studies based on the entire transcriptome or targeted genes of interest to sequence. The raw data in this study can be available in NCBI SRA database with accession number of SRP099564.

Cloning, characterization and mRNA expression of interleukin-6 in blunt snout bream (Megalobrama amblycephala).

In the present study, the interleukin-6 gene (IL-6) cDNA in blunt snout bream (Megalobrama amblycephala) was identified and its expression profiles under ammonia stress and bacterial challenge were investigated. The IL-6 sequence consisted of 1045 bp, including a 696 bp ORF which translated into a 232 amino acid (AA) protein. The protein contained a putative signal peptide of 24 AA in length. IL-6 expression analysis showed that the it is differentially expressed in various tissues under normal conditions and the highest IL-6 level was observed in the intestine tissue, followed by the liver, and then in the gills. Under ammonia stress, the IL-6 mRNA level both in spleens and intestine increased significantly (P < 0.05), with the maximum levels attained at 6 h, 12 h (72, 10-fold, respectively). Thereafter, they all significantly decreased (P < 0.01) and returned to the basal value within 48 h. Whereas, in livers it slightly decreased at 3 h firstly (0.5-fold), and then significantly (P < 0.05) increased with the maximum level attained 12 h (3-fold). Further expression analysis showed that the mRNA level of IL-6 in spleens, intestine and livers of blunt snout bream all increased significantly (P < 0.05), with maximum values attained at 6 h, 3 h, 6 h (10, 6, 18-fold, respectively) after Aeromonas hydrophila (A. hydrophila) injection, and then decreased to the basal value within 24 h which suggested that IL-6 was involved in the immune response to A. hydrophila. The cloning and expression analysis of the IL-6 provide theoretical basis to further study the mechanism of anti-adverseness and expression characteristics under stress conditions in blunt snout bream.

Development and application of monoclonal antibodies for detection and analysis of aquareoviruses.

Monoclonal antibodies (mAbs) play an important role in detection of aquareoviruses. Three mAbs against grass carp reovirus (GCRV) were prepared. Isotyping revealed that all three mAbs were of subclass IgG2b. Western blot assay showed that all three mAbs reacted with GCRV 69 kDa protein (the putative VP5). In addition to the 69 kDa protein of GCRV, mAb 4B6 also recognize a 54 kDa protein. All three mAbs were used for detecting aquareovirus by Western blot assay and indirect immunofluorescence assay (IFA). All of them reacted with GCRV, and mAb 4A3 could also react with turbot Scophthalmus maximus reovirus (SMReV) and largemouth bass Microptererus salmonides reovirus (MsReV). Viral antigens were only observed in the cytoplasm of infected cells. Finally, syncytia formation was observed with light microscopy and fluorescence microscopy using fluorescein labelled 4A3 mAb at various times post-infection. Syncytia were observed at 36 hr post-infection (hpi) by light microscopy and at 12 hpi by fluorescence microscopy. The immunofluorescence based assay allowed earlier detection of virus than observation of virus-induced cytopathic effect (CPE) assay in inoculated cell cultures. The sensitivity and specificity of these mAbs may be useful for diagnosis and monitoring of aquareoviruses.

Fish reovirus GCReV-109 VP33 protein elicits protective immunity in rare minnows.

Grass carp reovirus strain 109 (GCReV-109) was previously isolated from a grass carp (Ctenopharyngodon idellus) with hemorrhagic disease, and its complete genome has been sequenced. However, the infectivity of GCReV-109 has not been studied, and the viral protein VP33, encoded on genome segment S11, had no detectable sequence homology to other known reovirus proteins. In this study, we characterized GCReV-109 infections in vivo and in vitro, as well as the VP33 protein. Infectivity analysis showed that GCReV-109 caused severe hemorrhagic disease and 100% mortality at dilutions up to 10(-4) in rare minnows (Gobiocypris rarus) by 8 days postinfection, but no visible cytopathic effect was observed in GCReV-109-infected subcultured grass carp muscle (GCM) cells. To confirm that GCReV-109 could be propagated in GCM cells, three virus genome segments were detected by RT-PCR, and large numbers of virus particles were observed by transmission electron microscopy in samples from the infected GCM cells. The suspension of GCReV-109-infected GCM cells was pathogenic to rare minnows. VP33 protein was expressed and purified for generation of an anti-VP33 antiserum. In western blot analysis of purified GCReV-109 particles, the antiserum specifically recognized a protein band (approximately 33 kDa). This revealed that VP33 is a major structural protein of GCReV-109 that might have immunogenic properties. The protective efficacy of the anti-VP33 antiserum against GCReV-109 infection was tested. The death of infected fish was delayed and the mortality fell to 10% when fish were treated with the anti-VP33 antiserum, suggesting that it might be useful for the prevention and control of fish reoviral disease.

Genome analysis and gene nblA identification of Microcystis aeruginosa myovirus (MaMV-DC) reveal the evidence for horizontal gene transfer events between cyanomyovirus and host.

The genome sequence, genetic characterization and nblA gene function of Microcystis aeruginosa myovirus isolated from Lake Dianchi in China (MaMV-DC) have been analysed. The genome DNA is 169 223 bp long, with 170 predicted protein-coding genes (001L­170L) and a tRNA gene. About one-sixth of these genes have homologues in the host cyanobacteria M. aeruginosa. The genome carries a gene homologous to host nblA, which encodes a protein involved in the degradation of cyanobacterial phycobilisome. Its expression during MaMV-DC infection was confirmed by reverse transcriptase PCR and Western blot detection and abundant expression was companied by the significant decline of phycocyanin content and massive release of progeny MaMV-DC. In addition, expressing MaMV-DC nblA reduced the phycocyanin peak and the phycocyanin to chlorophyll ratio in model cyanobacteria. These results confirm that horizontal gene transfer events have occurred between cyanobacterial host and cyanomyovirus and suggest that MaMV-DC carrying host-derived genes (such as 005L, that codes for NblA) is responsible for more efficient expression of cyanophage genes and release of progeny cyanophage. This study provides novel insight into the horizontal gene transfer in cyanophage and the interactions between cyanophage and their host.

Whole-Genome Analysis of a Novel Fish Reovirus (MsReV) Discloses Aquareovirus Genomic Structure Relationship with Host in Saline Environments.

Aquareoviruses are serious pathogens of aquatic animals. Here, genome characterization and functional gene analysis of a novel aquareovirus, largemouth bass Micropterus salmoides reovirus (MsReV), was described. It comprises 11 dsRNA segments (S1-S11) covering 24,024 bp, and encodes 12 putative proteins including the inclusion forming-related protein NS87 and the fusion-associated small transmembrane (FAST) protein NS22. The function of NS22 was confirmed by expression in fish cells. Subsequently, MsReV was compared with two representative aquareoviruses, saltwater fish turbot Scophthalmus maximus reovirus (SMReV) and freshwater fish grass carp reovirus strain 109 (GCReV-109). MsReV NS87 and NS22 genes have the same structure and function with those of SMReV, whereas GCReV-109 is either missing the coiled-coil region in NS79 or the gene-encoding NS22. Significant similarities are also revealed among equivalent genome segments between MsReV and SMReV, but a difference is found between MsReV and GCReV-109. Furthermore, phylogenetic analysis showed that 13 aquareoviruses could be divided into freshwater and saline environments subgroups, and MsReV was closely related to SMReV in saline environments. Consequently, these viruses from hosts in saline environments have more genomic structural similarities than the viruses from hosts in freshwater. This is the first study of the relationships between aquareovirus genomic structure and their host environments.

Establishment of three cell lines from Chinese giant salamander and their sensitivities to the wild-type and recombinant ranavirus.

Known as lethal pathogens, Ranaviruses have been identified in diseased fish, amphibians (including Chinese giant salamander Andrias davidianus, the world's largest amphibian) and reptiles, causing organ necrosis and systemic hemorrhage. Here, three Chinese giant salamander cell lines, thymus cell line (GSTC), spleen cell line (GSSC) and kidney cell line (GSKC) were initially established. Their sensitivities to ranaviruses, wild-type Andrias davidianus ranavirus (ADRV) and recombinant Rana grylio virus carrying EGFP gene (rRGV-EGFP) were tested. Temporal transcription pattern of ranavirus major capsid protein (MCP), fluorescence and electron microscopy observations showed that both the wild-type and recombinant ranavirus could replicate in the cell lines.

Unraveling the genome structure of cyanobacterial podovirus A-4L with long direct terminal repeats.

The freshwater cyanobacterial virus (cyanophage) A-4L, a podovirus, can infect the model cyanobacterium Anabaena sp. strain PCC 7120 resulting in a high burst size and forming concentric plaques on its lawns. The complete genome sequence of A-4L was determined by the combination of high-throughput sequencing, terminal transferase-mediated polymerase chain reaction and restriction mapping. It contains 41,750 bp with 810 bp direct terminal repeats and 38 potential open reading frames. As compared with other cyanobacterial podoviruses in diverse ecosystems, the A-4L has the longest terminal repeat and shares similar genome organizations with freshwater members. Furthermore, phylogenetic analysis based on concatenated sequences of eight core proteins indicated that freshwater cyanobacterial podoviruses were clustered together and distinct from marine counterparts, suggesting a clear divergence in the cyanobacterial podovirus lineage between freshwater and marine ecosystems. Our findings uncover the unique genome structure of A-4L which contains long direct terminal repeats, and create the first model system to address knowledge gaps in understanding cyanobacterial virus-host interactions at the molecular level.

Rana grylio virus (RGV) envelope protein 2L: subcellular localization and essential roles in virus infectivity revealed by conditional lethal mutant.

Rana grylio virus (RGV) is a pathogenic iridovirus that has resulted in high mortality in cultured frog. Here, an envelope protein gene, 2L, was identified from RGV and its possible role in virus infection was investigated. Database searches found that RGV 2L had homologues in all sequenced iridoviruses and is a core gene of iridoviruses. Western blotting detection of purified RGV virions confirmed that 2L protein was associated with virion membrane. Fluorescence localization revealed that 2L protein co-localized with viral factories in RGV infected cells. In co-transfected cells, 2L protein co-localized with two other viral envelope proteins, 22R and 53R. However, 2L protein did not co-localize with the major capsid protein of RGV in co-transfected cells. Meanwhile, fluorescence observation showed that 2L protein co-localized with endoplasmic reticulum, but did not co-localize with mitochondria and Golgi apparatus. Moreover, a conditional lethal mutant virus containing the lac repressor/operator system was constructed to investigate the role of RGV 2L in virus infection. The ability to form plaques and the virus titres were strongly reduced when expression of 2L was repressed. Therefore, the current data showed that 2L protein is essential for virus infection. Our study is the first report, to our knowledge, of co-localization between envelope proteins in iridovirus and provides new insights into the understanding of envelope proteins in iridovirus.

A conditional lethal mutation in Rana grylio virus ORF 53R resulted in a marked reduction in virion formation.

Rana grylio virus (RGV) is a pathogenic iridovirus that has resulted in high mortality in cultured frog. Here, a recombinant RGV (i53R-RGV-lacIO) containing the inducible lac repressor/operator system was constructed. i53R-RGV-lacIO was a conditional lethal mutant in which the expression of envelope protein 53R was regulated by IPTG. i53R-RGV-lacIO shared characteristics similar to RGV in the presence of IPTG. However, the expression level of 53R, the ability of plaques formation, and the virus titers were strongly reduced in the absence of IPTG. Electron microscopy showed that the number of progeny virus produced by i53R-RGV-lacIO was remarkably reduced without IPTG. Furthermore, over-expression of 53R in vitro could increase titers of i53R-RGV-lacIO in the absence of IPTG. Therefore, the current data suggested that the lac repressor/operator system could regulate gene expression in the recombinant iridovirus. Our study was thought to be the first report of the system in aquatic virus.