First Report of a Ranavirus Associated with Morbidity and Mortality in Farmed Chinese Giant Salamanders (Andrias davidianus)

Nitrogen gas-bubble disease in two giant salamanders

Gas-bubble disease (GBD)-a non-infectious disease in aquatic organisms caused by supersaturated levels of total dissolved gases (oxygen and nitrogen) in water-is well known in various species, including fish and amphibians, but has not previously been reported in giant salamanders. In the present study, macroscopic and histopathological examinations of 2 mature Andrias spp. (kept with 293 fish in an aquarium) were performed to characterize GBD pathology. Bubbles developed on the body surfaces of the salamanders and fish, with erythema specifically noted in the salamanders. Within 3 d of the bubbles appearing, both salamanders and more than 270 fish had died. On Days 1 and 2, dissolved oxygen levels were 75.5 and 86.9%, respectively, while dissolved nitrogen gas levels were 90.6 and 103.1%, respectively. The 2 salamanders exhibited identical lesions characterized by erythema, congestion, and numerous bubbles in the major veins of the body cavity. Histopathologically, congestion and gas embolism-like dilatations were observed in the small vessels and capillaries. These lesions were found in the parenchymal and gastrointestinal organs, skin, eyeballs, and surrounding stromal tissue. Based on these findings and that GBD occurs at dissolved nitrogen gas and oxygen levels above 120 and 200%, respectively, the salamanders were diagnosed with nitrogen GBD. The exact etiology of this disease remains unconfirmed but likely involves circulatory system dysfunction within the aquarium environment, highlighting the importance of routine inspections and maintenance of equipment.

Chinese giant salamander Bcl-w: An inhibitory role in iridovirus-induced mitochondrial apoptosis and virus replication

B-cell lymphoma-2 (BCL-2) superfamily molecules play crucial roles in mitochondrial apoptosis induced by Chinese giant salamander iridovirus (GSIV). As an anti-apoptotic molecule in the BCL-2 family, the molecular mechanism of Bcl-w during GSIV infection remains unknown. In this study, we characterized for the first time an amphibian Bcl-w from Chinese giant salamander Andrias davidianus (AdBcl-w), and its function and regulatory mechanism during GSIV infection were investigated. AdBcl-w possesses the conserved structural features of Bcl-w and shares 35-54% sequence identities with other Bcl-w. mRNA expression of AdBcl-w was most abundant in liver and muscle. The AdBcl-w mRNA expression was regulated during GSIV infection. Western blotting assays revealed that the level of Bcl-w protein was downregulated markedly as the infection progresses. Confocal microscopy showed that overexpressed AdBcl-w was translocated to the mitochondria after infection with GSIV. Flow cytometry analysis demonstrated that compared with control, the apoptotic progress in cells transfected with AdBcl-w was reduced while that in cells transfected with AdBcl-w siRNA was enhanced. The number of virus major capsid protein gene copies was lower and protein synthesis was reduced in AdBcl-w overexpressing cells. In addition, AdBcl-w could bind directly to the pro-apoptotic molecule AdBak, while this interaction was weakened with GSIV infection. Moreover, p53 level was reduced and the mRNA expression levels of crucial regulatory molecules in the p53 pathway were regulated in AdBcl-w overexpressing cells during GSIV infection. These results suggested that AdBcl-w inhibit GSIV replication by regulating the virus induced mitochondrial apoptosis.

Fusing a TurboID tag with the Andrias davidianus ranavirus 2L reduced virus adsorption efficiency

Andrias davidianus ranavirus (ADRV) is a member of the genus ranavirus (family Iridoviridae). ADRV 2L is an envelope protein that could be essential in viral infection. In the present study, the function of ADRV 2L was investigated by fusion with the biotin ligase TurboID tag. A recombinant ADRV with a V5-TurboID tag fused in the N-terminal of 2L (ADRVT-2L) and a recombinant ADRV expressing V5-TurboID (ADRVT) were constructed, respectively. Infection of the recombinant viruses and wild-type ADRV (ADRVWT) in the Chinese giant salamander thymus cell line (GSTC) showed that ADRVT-2L had reduced cytopathic effect and lower virus titers than the other two viruses, indicating the fusion of a big tag affected ADRV infection. Analysis of the temporal expression profile showed that the expression of V5-TurboID-2L was delayed than wild-type 2L. However, electron microscopy found that the virion morphogenesis was not affected in ADRVT-2L-infected cells. Furthermore, the virus binding assay revealed that the adsorption efficiency of ADRVT-2L was considerably decreased compared to the other two viruses. Therefore, these data showed that linking the TurboID tag to ADRV 2L affected virus adsorption to the cell membrane, which suggested an important role of 2L in virus entry into cells.

Disease state associated with chronic toe lesions in hellbenders may alter anti-chytrid skin defenses

Hellbenders (Cryptobranchus alleganiensis) are large, aquatic salamanders from the eastern United States. Both subspecies, eastern and Ozark hellbenders, have experienced declines resulting in federal listing of Ozark hellbenders. The globally distributed chytrid fungus, Batrachochytrium dendrobatidis (Bd) has been detected in both subspecies, and Batrachochytrium salamandrivorans (Bsal) poses a new threat if introduced into North America. Ozark hellbenders also suffer a high prevalence of toe lesions of unknown etiology, with changes in host immunocompetence hypothesized to contribute. Antimicrobial peptides (AMPs) secreted from dermal granular glands may play a role in hellbender health. We collected skin secretions from free-ranging hellbenders and enriched them for small cationic peptides used for growth inhibition assays against Bd and Bsal. Generalized linear mixed models revealed the presence of active toe lesions as the strongest and only significant predictor of decreased Bd inhibition by skin peptides. We also found skin secretions were more inhibitory of Bsal than Bd. MALDI-TOF mass spectrometry revealed candidate peptides responsible for anti-chytrid activity. Results support the hypothesis that hellbender skin secretions are important for innate immunity against chytrid pathogens, and decreased production or release of skin peptides may be linked to other sub-lethal effects of disease associated with toe lesions.

The first identification of three AdIRAK2 genes from an evolutionarily important amphibian Andrias davidianus and their involvement in NF-κB activation and inflammatory responses

Interleukin-1 receptor associated kinases (IRAK) is the most important downstream kinases of TLRs/IL-1R signaling pathway for signal transduction and activation of inflammatory response against pathogen infections. However, the molecular identification and function characterization of IRAK2 homologs in lower vertebrate remains obscure. In this study, three IRAK2 genes (AdIRAK2a, AdIRAKb and AdIRAK2c) and their respective transcripts were identified from the Chinese giant salamander Andrias davidianus. This is the first evidence that three different IRAK2 genes exist in an ancient amphibian species, which has never been reported in other vertebrates. The complete open reading frames (ORFs) of AdIRAK2a, AdIRAK2b and AdIRAK2c were 2112 bp, 1917 bp and 816 bp encoding deduced proteins of 703 amino acids (aa), 628 aa and 271 aa, respectively. All three AdIRAK2 proteins contained two predicted conserved functional domains, including a death domain (DD) and a serine/threonine protein kinases domain (KD). Phylogenetic analysis showed that the three AdIRAK2s clustered together with other known IRAK2 of vertebrates. The three AdIRAK2s were ubiquitously expressed in all tested tissues with a similar tissues distribution pattern. After challenge of Aeromonas hydrophila (A. hydrophila), Staphylococcus aureus (S.aureus), giant salamander iridovirus (GSIV, belonging to the genus Ranavirus in the family Iridoviridae) and polyinosinic:polycytidylic acid (poly(I:C)), the expression levels of all AdIRAK2s in blood were significantly altered, however, they exhibited distinct response patterns. Moreover, the results of over-expression and RNAi of AdIRAK2s implied the involvement of AdIRAK2s in triggering NF-κB-mediated signaling pathways and inflammatory responses. This study might provide a better understanding of the presence and immune regulation function of IRAK2 in amphibians and even in vertebrates.

Is Xenopus laevis introduction linked with Ranavirus incursion, persistence and spread in Chile?

Ranaviruses have been associated with amphibian, fish and reptile mortality events worldwide and with amphibian population declines in parts of Europe. Xenopus laevis is a widespread invasive amphibian species in Chile. Recently, Frog virus 3 (FV3), the type species of the Ranavirus genus, was detected in two wild populations of this frog near Santiago in Chile, however, the extent of ranavirus infection in this country remains unknown. To obtain more information about the origin of ranavirus in Chile, its distribution, species affected, and the role of invasive amphibians and freshwater fish in the epidemiology of ranavirus, a surveillance study comprising wild and farmed amphibians and wild fish over a large latitudinal gradient (2,500 km) was carried out in 2015-2017. In total, 1,752 amphibians and 496 fish were tested using a ranavirus-specific qPCR assay, and positive samples were analyzed for virus characterization through whole genome sequencing of viral DNA obtained from infected tissue. Ranavirus was detected at low viral loads in nine of 1,011 X. laevis from four populations in central Chile. No other amphibian or fish species tested were positive for ranavirus, suggesting ranavirus is not threatening native Chilean species yet. Phylogenetic analysis of partial ranavirus sequences showed 100% similarity with FV3. Our results show a restricted range of ranavirus infection in central Chile, coinciding with X. laevis presence, and suggest that FV3 may have entered the country through infected X. laevis, which appears to act as a competent reservoir host, and may contribute to the spread the virus locally as it invades new areas, and globally through the pet trade.

DNAJA4 Promotes the Replication of the Chinese Giant Salamander Iridovirus

The DNAJ family, a class of chaperone proteins involved in protein folding, assembly, and transport, plays an essential role in viral infections. However, the role of DNAJA4 (DnaJ Heat Shock Protein Family (Hsp40) Member A4) in the ranavirus infection has not been reported. This study demonstrates the function of the epithelial papilloma of carp (EPC) DNAJA4 in Chinese giant salamander (Andrias davidianus) iridovirus (CGSIV) replication. DNAJA4 consists of 1479 base pairs and encodes a 492 amino acid polypeptide. Sequence analysis has shown that EPC DNAJA4 contains a conserved J domain and shares 84% homology with Danio rerio DNAJA4 and 68% homology with Homo sapiens DNAJA4. EPC DNAJA4 was localized in the cytoplasm, and its expression was significantly upregulated after CGSIV infection. Overexpression of EPC DNAJA4 promotes CGSIV replication and CGSIV DNA replication. siRNA knockdown of DNAJA4 expression attenuates CGSIV replication and viral DNA replication. Overexpression and interference experiments have proved that EPC DNAJA4 is a pro-viral factor. Co-IP, GST-pulldown, and immunofluorescence confirmed the interaction between EPC DNAJA4 and CGSIV proliferating cell nuclear antigen (PCNA). Our results demonstrate for the first time that EPC DNAJA4 is involved in viral infection by promoting viral DNA replication and interacting with proteins associated with viral replication.

Glycosylphosphatidylinositol Mannosyltransferase Ⅰ Protects Chinese Giant Salamander, Andrias davidianus, against Iridovirus

Glycosylphosphatidylinositol mannosyltransferase I (GPI-MT-I) is an essential glycosyltransferase of glycosylphosphatidylinositol-anchor proteins (GPI-APs) that transfers the first of the four mannoses in GPI-AP precursors, which have multiple functions, including immune response and signal transduction. In this study, the GPI-MT-I gene that regulates GPI-AP biosynthesis in Andrias davidianus (AdGPI-MT-I) was characterized for the first time. The open reading frame (ORF) of AdGPI-MT-I is 1293 bp and encodes a protein of 430 amino acids that contains a conserved PMT2 superfamily domain. AdGPI-MT-I mRNA was widely expressed in the tissues of the Chinese giant salamander. The mRNA expression level of AdGPI-MT-I in the spleen, kidney, and muscle cell line (GSM cells) was significantly upregulated post Chinese giant salamander iridovirus (GSIV) infection. The mRNA expression of the virus major capsid protein (MCP) in AdGPI-MT-I-overexpressed cells was significantly reduced. Moreover, a lower level of virus MCP synthesis and gene copying in AdGPI-MT-I-overexpressed cells was confirmed by western blot and ddPCR. These results collectively suggest that GSIV replication in GSM cells was significantly reduced by the overexpression of the AdGPI-MT-I protein, which may contribute to a better understanding of the antiviral mechanism against iridovirus infection.

Comparative genomic analysis reveals new evidence of genus boundary for family Iridoviridae and explores qualified hallmark genes

Members of the family Iridoviridae (iridovirids) are globally distributed and trigger adverse economic and ecological impacts on aquaculture and wildlife. Iridovirids taxonomy has previously been studied based on a limited number of genomes, but this is not suitable for the current and future virological studies as more iridovirids are emerging. In our study, 57 representative iridovirids genomes were selected from a total of 179 whole genomes available on NCBI. Then 18 core genes were screened out for members of the family Iridoviridae. Average amino acid sequence identity (AAI) analysis indicated that a cut-off value of 70% is more suitable for the current iridovirids genome database than ICTV-defined 50% threshold to better clarify viral genus boundaries. In addition, more subgroups were divided at genus level with the AAI threshold of 70%. This observation was further confirmed by genomic synteny analysis, codon usage preference analysis, genome GC content and length analysis, and phylogenic analysis. According to the pairwise comparison analysis of core genes, 9 hallmark genes were screened out to conduct preliminary identification and investigation at the genus level of iridovirids in a more convenient and economical manner.

Prevalence of Ranavirus Infection in Three Anuran Species across South Korea

To cope with amphibian die-offs caused by ranavirus, it is important to know the underlying ranavirus prevalence in a region. We studied the ranavirus prevalence in tadpoles of two native and one introduced anuran species inhabiting agricultural and surrounding areas at 49 locations across eight provinces of South Korea by applying qPCR. The local ranavirus prevalence and the individual infection rates at infected locations were 32.6% and 16.1%, respectively, for Dryophytes japonicus (Japanese tree frog); 25.6% and 26.1% for Pelophylax nigromaculatus (Black-spotted pond frog); and 30.5% and 50.0% for Lithobates catesbeianus (American bullfrog). The individual infection rate of L. catesbeianus was significantly greater than that of D. japonicus. The individual infection rate of P. nigromaculatus was related to the site-specific precipitation and air temperature. The individual infection rate gradually increased from Gosner development stage 39, and intermittent infection was confirmed in the early and middle developmental stages. Our results show that ranavirus is widespread among wild amphibians living in agricultural areas of South Korea, and mass die-offs by ranavirus could occur at any time.

Bcl-xL Reduces Chinese Giant Salamander Iridovirus-Induced Mitochondrial Apoptosis by Interacting with Bak and Inhibiting the p53 Pathway

Chinese giant salamander iridovirus (GSIV) infection could lead to mitochondrial apoptosis in this animal, a process that involves B-cell lymphoma-2 (BCL-2) superfamily molecules. The mRNA expression level of Bcl-xL, a crucial antiapoptotic molecule in the BCL-2 family, was reduced in early infection and increased in late infection. However, the molecular mechanism remains unknown. In this study, the function and regulatory mechanisms of Chinese giant salamander (Andrias davidianus) Bcl-xL (AdBcl-xL) during GSIV infection were investigated. Western blotting assays revealed that the level of Bcl-xL protein was downregulated markedly as the infection progressed. Plasmids expressing AdBcl-xL or AdBcl-xL short interfering RNAs were separately constructed and transfected into Chinese giant salamander muscle cells. Confocal microscopy showed that overexpressed AdBcl-xL was translocated to the mitochondria after infection with GSIV. Additionally, flow cytometry analysis demonstrated that apoptotic progress was reduced in both AdBcl-xL-overexpressing cells compared with those in the control, while apoptotic progress was enhanced in cells silenced for AdBcl-xL. A lower number of copies of virus major capsid protein genes and a reduced protein synthesis were confirmed in AdBcl-xL-overexpressing cells. Moreover, AdBcl-xL could bind directly to the proapoptotic molecule AdBak with or without GSIV infection. In addition, the p53 level was inhibited and the mRNA expression levels of crucial regulatory molecules in the p53 pathway were regulated in AdBcl-xL-overexpressing cells during GSIV infection. These results suggest that AdBcl-xL plays negative roles in GSIV-induced mitochondrial apoptosis and virus replication by binding to AdBak and inhibiting p53 activation.

Patterns of infection, origins, and transmission of ranaviruses among the ectothermic vertebrates of Asia

Ranaviral infections, a malady of ectothermic vertebrates, are becoming frequent, severe, and widespread, causing mortality among both wild and cultured species, raising odds of species extinctions and economic losses. This increase in infection is possibly due to the broad host range of ranaviruses and the transmission of these pathogens through regional and international trade in Asia, where outbreaks have been increasingly reported over the past decade. Here, we focus attention on the origins, means of transmission, and patterns of spread of this infection within the region. Infections have been recorded in both cultured and wild populations in at least nine countries/administrative regions, together with mass die-offs in some regions. Despite the imminent seriousness of the disease in Asia, surveillance efforts are still incipient. Some of the viral strains within Asia may transmit across host-taxon barriers, posing a significant risk to native species. Factors such as rising temperatures due to global climate change seem to exacerbate ranaviral activity, as most known outbreaks have been recorded during summer; however, data are still inadequate to verify this pattern for Asia. Import risk analysis, using protocols such as Pandora+, pre-border pathogen screening, and effective biosecurity measures, can be used to mitigate introduction of ranaviruses to uninfected areas and curb transmission within Asia. Comprehensive surveillance using molecular diagnostic tools for ranavirus species and variants will help in understanding the prevalence and disease burden in the region. This is an important step toward conserving native biodiversity and safeguarding the aquaculture industry.

The Immune System and the Antiviral Responses in Chinese Giant Salamander, Andrias davidianus

The Chinese giant salamander, belonging to an ancient amphibian lineage, is the largest amphibian existing in the world, and is also an important animal for artificial cultivation in China. However, some aspects of the innate and adaptive immune system of the Chinese giant salamander are still unknown. The Chinese giant salamander iridovirus (GSIV), a member of the Ranavirus genus (family Iridoviridae), is a prominent pathogen causing high mortality and severe economic losses in Chinese giant salamander aquaculture. As a serious threat to amphibians worldwide, the etiology of ranaviruses has been mainly studied in model organisms, such as the Ambystoma tigrinum and Xenopus. Nevertheless, the immunity to ranavirus in Chinese giant salamander is distinct from other amphibians and less known. We review the unique immune system and antiviral responses of the Chinese giant salamander, in order to establish effective management of virus disease in Chinese giant salamander artificial cultivation.

Functional differences in the products of two TRAF3 genes in antiviral responses in the Chinese giant salamander, Andrias davidianus

Tumour necrosis factor receptor associated factor 3 (TRAF3) is a crucial transducing protein for linking upstream receptor signals and downstream antiviral signalling pathways. Previous studies mostly clarified the functions of TRAF3 in mammals, birds and fish, but little is known about the characterization and function of TRAF3 in amphibians. In this study, the molecular and functional identification of two TRAF3 genes, AdTRAF3A and AdTRAF3B, were investigated in the Chinese giant salamander Andrias davidianus. The complete open reading frames (ORFs) of AdTRAF3A and AdTRAF3B were 1698 bp and 1743 bp in length, encoding 565 and 580 amino acids, respectively. Both AdTRAF3A and AdTRAF3B deduced proteins contained a RING finger, two TRAF-type zinc fingers, a coiled-coil and a MATH domain. Phylogenetic analysis showed that the AdTRAF3 protein clustered together with other known TRAF3 proteins. Gene expression analysis showed that AdTRAF3s were broadly distributed in all examined tissues with similar distribution patterns. AdTRAF3s in the blood or spleen positively responded to Giant salamander iridovirus (GSIV) and poly (I:C) induction but exhibited distinct response patterns. Silencing AdTRAF3A/B remarkably suppressed the expression of IFN signalling pathway-related genes when leukocytes were treated with DNA virus and the viral RNA analogue. Moreover, overexpression of AdTRAF3A may induce the activation of the IFN-β promoter, and the zinc finger, coiled coil and MATH domains of AdTRAF3A were essential for IFN-β promoter activation. However, the overexpression of AdTRAF3B significantly suppressed IFN-β promoter activity, and its inhibitory effect was enhanced when the RING finger or MATH domain was deleted. Furthermore, AdTRAF3A rather than AdTRAF3B significantly induced NF-κB activation, implying that AdTRAF3A may function as an enhancer in both the IFN and NF-κB signalling pathways. Taken together, our results suggest that the two TRAF3 genes play different crucial regulatory roles in innate antiviral immunity in Chinese giant salamanders.

GEOGRAPHIC AND INDIVIDUAL DETERMINANTS OF IMPORTANT AMPHIBIAN PATHOGENS IN HELLBENDERS (CRYPTOBRANCHUS ALLEGANIENSIS) IN TENNESSEE AND ARKANSAS, USA

Wildlife diseases are a major threat for species conservation and there is a growing need to implement disease surveillance programs to protect species of concern. Globally, amphibian populations have suffered considerable losses from disease, particularly from chytrid fungi (Batrachochytrium dendrobatidis [Bd] and Batrachochytrium salamandrivorans [Bsal]) and ranavirus. Hellbenders (Cryptobranchus alleganiensis) are large riverine salamanders historically found throughout several watersheds of the eastern and midwestern US. Populations of both subspecies (Ozark hellbender, Cryptobranchus alleganiensis bishopi; eastern hellbender, Cryptobranchus alleganiensis alleganiensis) have experienced precipitous declines over at least the past five decades, and emerging pathogens are hypothesized to play a role. We surveyed Ozark hellbender populations in Arkansas (AR) and eastern hellbender populations in Middle Tennessee (MTN) and East Tennessee (ETN) for both chytrid fungi and ranavirus from swabs and tail tissue, respectively, from 2011 to 2017. Overall, we detected Bd on hellbenders from nine out of 15 rivers, with total prevalence of 26.7% (54/ 202) that varied regionally (AR: 33%, 28/86; MTN: 11%, 4/36; ETN: 28%, 22/80). Ranavirus prevalence (9.0%, 18/200) was comparatively lower than Bd, with less regional variation in prevalence (AR: 6%, 5/ 85; MTN: 11%, 4/36; ETN: 10%, 8/79). We did not detect Bsal in any hellbender populations. We detected a significant negative correlation between body condition score and probability of ranavirus infection (β=-0.13, SE=0.06, 95% confidence interval: -0.24, -0.02). Evaluation of infection load of positive individuals revealed different trends than prevalence alone for both ranavirus and Bd, with MTN having a significantly greater average ranaviral load than both other regions. We documented a variety of lesions that likely have multiple etiologies on hellbenders located within all geographic regions. Our data represent a multiyear pathogen dataset across several regions of C. alleganiensis, and we emphasize the need for continued pathogen surveillance.

The Insights of Genomic Synteny and Codon Usage Preference on Genera Demarcation of Iridoviridae Family

The members of the family Iridoviridae are large, double-stranded DNA viruses that infect various hosts, including both vertebrates and invertebrates. Although great progress has been made in genomic and phylogenetic analyses, the adequacy of the existing criteria for classification within the Iridoviridae family remains unknown. In this study, we redetermined 23 Iridoviridae core genes by re-annotation, core-pan analysis and local BLASTN search. The phylogenetic tree based on the 23 re-annotated core genes (Maximum Likelihood, ML-Tree) and amino acid sequences (composition vector, CV-Tree) were found to be consistent with previous reports. Furthermore, the information provided by synteny analysis and codon usage preference (relative synonymous codon usage, correspondence analysis, ENC-plot and Neutrality plot) also supports the phylogenetic relationship. Collectively, our results will be conducive to understanding the genera demarcation within the Iridoviridae family based on genomic synteny and component (codon usage preference) and contribute to the existing taxonomy methods for the Iridoviridae family.

Bid is involved in apoptosis induced by Chinese giant salamander iridovirus and contributes to the viral replication in an amphibian cell line

Bid is a pro-apoptotic BH3-only member of the Bcl-2 superfamily that functions to link the extrinsic apoptotic pathway and the mitochondrial amplification loop of the intrinsic pathway. In this study, the expression and functions of Chinese giant salamander (Andrias davidianus) Bid (AdBid) were investigated. The AdBid cDNA sequence contains an open reading frame (ORF) of 576 nucleotides, encoding a putative protein of 191 aa. AdBid possesses the conserved BH3 interacting domain and shared 34-52% sequence identities with other amphibian Bid. mRNA expression of AdBid was most abundant in muscle. The expression level of AdBid in Chinese giant salamander muscle, kidney and spleen significantly increased after Chinese giant salamander iridovirus (GSIV) infection. Additionally, a plasmid expressing AdBid was constructed and transfected into the Chinese giant salamander muscle cell line (GSM cells). The morphology and cytopathic effect (CPE) and apoptotic process in AdBid over-expressed GSM cells was significantly enhanced during GSIV infection compared with that in control cells. Moreover, a higher level of the virus major capsid protein (MCP) gene copies and protein synthesis was confirmed in the AdBid over-expressed cells. These results indicated that AdBid played a positive role in GSIV induced apoptosis and the viral replication. This study may contribute to the better understanding on the infection mechanism of iridovirus-induced apoptosis.

Range-wide decline of Chinese giant salamanders Andrias spp. from suitable habitat

Over recent decades, Chinese giant salamanders Andrias spp. have declined dramatically across much of their range. Overexploitation and habitat degradation have been widely cited as the cause of these declines. To investigate the relative contribution of each of these factors in driving the declines, we carried out standardized ecological and questionnaire surveys at 98 sites across the range of giant salamanders in China. We did not find any statistically significant differences between water parameters (temperature, dissolved oxygen, ammonia, nitrite, nitrate, salinity, alkalinity, hardness and flow rate) recorded at sites where giant salamanders were detected by survey teams and/or had been recently seen by local respondents, and sites where they were not detected and/or from which they had recently been extirpated. Additionally, we found direct and indirect evidence that the extraction of giant salamanders from the wild is ongoing, including within protected areas. Our results support the hypothesis that the decline of giant salamanders across China has been primarily driven by overexploitation. Data on water parameters may be informative for the establishment of conservation breeding programmes, an initiative recommended for the conservation of these species.

Differential liver histopathological responses to amphibian chytrid infection

Amphibians have been facing a pandemic caused by the deadly fungus Batrachochytrium dendrobatidis (Bd). Although studies have elucidated cutaneous and homeostatic disturbances, it is still unknown if the hepatic function can be affected or if hepatic effects differ among host species. Thus, we evaluated the effects of an experimental Bd infection on the liver (histopathology and the hepatosomatic index) of 2 anuran species (Xenopus laevis and Physalaemus albonotatus) with different susceptibilities to Bd infection and compared them to uninfected controls. Bd infection increased the melanomacrophage cell area and induced leukocyte infiltration in both species. The effects were more pronounced in the sensitive species, P. albonotatus, which showed severe reduction in glycogen stores and liver atrophy, due to energetic imbalance. Hepatocytes of P. albonotatus also showed ballooning degeneration (vacuolization), which could lead to cell death and liver failure. Our results provide evidence that although the sensitive species showed more severe effects, the tolerant species also had hepatic responses to the infection. These findings indicate that hepatic function can play an important role in detoxification and in immune responses to chytridiomycosis, and that it may be used as a new biomarker of health status in chytrid infections.

Selected Emerging Infectious Diseases of Amphibians

This article updates the understanding of two extirpation-driving infectious diseases, Batrachochytrium dendrobatidis and Batrachochytrium salamandrivorans, and Ranavirus. Experimental studies and dynamic, multifactorial population modeling have outlined the epidemiology and future population impacts of B dendrobatidis, B salamandrivorans, and Ranavirus. New genomic findings on divergent fungal and viral pathogens can help optimize control and disease management strategies. Although there have been major advances in knowledge of amphibian pathogens, controlled studies are needed to guide population recovery to elucidate and evaluate transmission routes for several pathogens, examine environmental control, and validate new diagnostic tools to confirm the presence of disease.

FTR67, a member of the fish-specific finTRIM family, triggers IFN pathway and against spring viremia of carp virus

Tripartite motif (TRIM) proteins have attracted particular research interest because of their multiple functions in the antiviral innate immune response. TRIM proteins perform different functions during virus infection, some play a role in inhibiting while others play a role in promoting. In this study, we described a species-specific TRIM gene named ftr67. Analysis of tissue distribution showed that ftr67 was mainly expressed in the gill and liver in five examined tissues of zebrafish. The phylogenic analysis showed that ftr67 was closest to the grass carp TRIM67. Overexpression of ftr67 resulted in a significantly decreased SVCV entry and impaired SVCV replication in FHM cells. Furthermore, overexpression of ftr67 could significantly induce the upregulation of molecular sensor RIG-I, IRF3/7, IFN and ISGs. In addition, RING domain of ftr67 was a required part essential for the antiviral effect. In summary, our results demonstrated that the important role of ftr67 in regulating SVCV infection, which offers a potential target for development of anti-SVCV therapies.

Zebrafish NIK Mediates IFN Induction by Regulating Activation of IRF3 and NF-κB

Type I IFN mediates the innate immune system to provide defense against viral infections. NF-κB-inducing kinase (NIK) potentiates the basal activation of endogenous STING, which facilitates the recruitment of TBK1 with the ectopically expressed IRF3 to induce IFN production. Moreover, NIK phosphorylates IKKα and confers its ability to phosphorylate p100 (also known as NF-κB2) in mammals. Our study demonstrated that NIK plays a critical role in IFN production in teleost fish. It was found that NIK interacts with IKKα in the cytoplasm and that IKKα phosphorylates the NIK at the residue Thr432, which is different from the mammals. Overexpression of NIK caused the activation of IRF3 and NF-κB, which in turn led to the production of IFN and IFN-stimulated genes (ISGs). Furthermore, the ectopic expression of NIK was observed to be associated with a reduced replication of the fish virus, whereas silencing of endogenous NIK had an opposite effect in vitro. Furthermore, NIK knockdown significantly reduced the expression of IFN and key ISGs in zebrafish larvae after spring viremia of carp virus infection. Additionally, the replication of spring viremia of carp virus was enhanced in NIK knockdown zebrafish larvae, leading to a lower survival rate. In summary, our findings revealed a previously undescribed function of NIK in activating IFN and ISGs as a host antiviral response. These findings may facilitate the establishment of antiviral therapy to combat fish viruses.

Evaluation of Severity and Factors Contributing to Foot Lesions in Endangered Ozark Hellbenders, Cryptobranchus alleganiensis bishopi

Arkansas populations of Ozark Hellbenders, Cryptobranchus alleganiensis bishopi have declined precipitously over the past few decades and are now limited to a single river. Biologists have also observed an increase of distal limb lesions with unidentified etiology and unknown role in morbidity and mortality of the species in this location. We documented lesions and collected associated individual size class data and pathogen samples in Ozark Hellbenders of Arkansas (n = 73) from 2011 to 2014 with the following two objectives: (1) document spatiotemporal patterns and severity of lesions present in this last remaining Arkansas Ozark Hellbender population, and (2) determine if host factors and infection status are associated with lesion severity. A scoring system was created from 0 to 7 based on lesion observations. Linear mixed model regressions followed by AICc model evaluation were used to determine associations among infection status for amphibian pathogens Batrachochytrium dendrobatidis (Bd) and Ranavirus as well as individual biometrics on lesion score. We discovered 93.2% of Hellbenders had lesions characterized by digit swelling that often progressed toward toe-tip ulceration. In severe cases we observed digital necrosis progressing to digit loss. Any recaptured individuals had the same or worse lesion score from previous captures. The top predictive model for lesion severity included individual mass and Bd infection status with a significant, positive association of Bd with increased lesion severity (β = 0.87 ± 0.39 S.E., C.I.: 0.11, 1.63). Our findings highlight a widespread and progressive disease that is an important factor to consider for the future of Ozark Hellbenders. This syndrome is presumptively multifactorial, and future studies will benefit from investigating several factors of host, infectious agents, and environment and their roles in disease manifestation for the purpose of developing effective, multi-faceted conservation strategies. A summary of potential etiologies and mechanisms is provided that may explain observed lesion distribution and that will be applicable to future disease and epidemiological investigations.

CMTV-like ranavirus infection associated with high mortality in captive catfish-like loach, Triplophysa siluorides, in China

Ranaviruses are important emerging pathogens of ectothermic vertebrates that threaten aquaculture and wildlife worldwide. A mortality event occurred in a cultured population of catfish-like loach (Triplophysa siluorides) in Sichuan Province, China. Gross clinical signs of the affected fish included skin lesions and haemorrhagic ulcers, which are often associated with ranaviruses. Inoculation of liver, kidney and spleen tissue homogenates in epithelioma papulosum cyprini (EPC) cells at 25°C resulted in cytopathic effect within 24 hr. Transmission electron microscopy of infected EPC cells revealed hexagonal viral arrays in the cytoplasm and icosahedral geometry of the virions. Following exposure of T. siluroides to the isolated virus, similar clinical signs were observed and the fish experienced 40% and 90% mortality after 21 days at 10^3.58 and 10^7.8 TCID₅₀ /0.1 ml doses, respectively, providing evidence the isolated virus was the main causative agent of the mortality event. Diagnostic PCR of the major capsid protein gene of ranavirus showed that all samples of diseased fish and isolated virus were positive. Phylogenetic analysis revealed that the isolated virus, designated as FYLl40220, was associated with the Common Midwife Toad Virus (CMTV)-like ranavirus clade. To our knowledge, this case represents the first report of CMTV-associated mortality in a fish species. Collectively, these results suggest that the host range of CMTV-like ranaviruses is greater than previously thought, and this clade of ranaviruses could have significant economic and biodiversity impacts.

An examination of the Iridovirus core genes for reconstructing Ranavirus phylogenies

Ranaviruses are globally emerging infections of poikilothermic vertebrates and belong to the viral family Iridoviridae. The six species of ranaviruses are responsible for unknown numbers of infections and disease and mortality events around the world in amphibians, fish, and reptiles. Genomic investigations have shown that there are 24 core genes shared by all iridoviruses. In this study, we examine the utility of each of these genes in reconstructing phylogenetic relationships across six species of Ranavirus. We also performed dot-plot analysis for the 17 isolates in the study. For large-scale differentiation, using the major capsid protein gene creates a tree similar to the whole genome tree. Other comparable genes include open reading frame (ORF) 19R (a serine–theonine protein kinase) and ORF 88R (Erv I/Alr Family protein). The poorest candidate for phylogenetic reconstruction, due to high homology, was ORF 1R (a putative replication factor and (or) DNA binding-packing protein). There are a plethora of genes that may be useful to examine phylogenies at smaller scales (e.g., to examine local adaptation); however, they do not necessarily belong to the set of highly conserved core genes.

When Micro Drives the Macro: A Fresh Look at Disease and its Massive Contributions in the Hindu Kush-Himalaya

The outbreaks of emerging and reemerging diseases have a high impact on the human and animal health because they are the underlying causes of disability, death, and long-term illness. For many regions those details are not, or just poorly known. Here we present on the morbidity and mortality in faunal diversities including domestic and wild species caused by various viral, bacterial, parasitic, and fungal diseases prevalent in Nepal and relevant for the wider Hindu Kush Himalaya. In addition, we provide details how antibiotic resistivity, vectors, and zoonosis have resulted on a landscape-scale in the huge public and veterinary health problem has been dealt with in the context of Nepal and the wider region.

Chinese Giant Salamander (Andrias davidianus) Iridovirus Infection Leads to Apoptotic Cell Death through Mitochondrial Damage, Caspases Activation, and Expression of Apoptotic-Related Genes

Chinese giant salamander iridovirus (GSIV) is the causative pathogen of Chinese giant salamander (Andrias davidianus) iridovirosis, leading to severe infectious disease and huge economic losses. However, the infection mechanism by GSIV is far from clear. In this study, a Chinese giant salamander muscle (GSM) cell line is used to investigate the mechanism of cell death during GSIV infection. Microscopy observation and DNA ladder analysis revealed that DNA fragmentation happens during GSIV infection. Flow cytometry analysis showed that apoptotic cells in GSIV-infected cells were significantly higher than that in control cells. Caspase 8, 9, and 3 were activated in GSIV-infected cells compared with the uninfected cells. Consistently, mitochondria membrane potential (MMP) was significantly reduced, and cytochrome c was released into cytosol during GSIV infection. p53 expression increased at an early stage of GSIV infection and then slightly decreased late in infection. Furthermore, mRNA expression levels of pro-apoptotic genes participating in the extrinsic and intrinsic pathway were significantly up-regulated during GSIV infection, while those of anti-apoptotic genes were restrained in early infection and then rose in late infection. These results collectively indicate that GSIV induces GSM apoptotic cell death involving mitochondrial damage, caspases activation, p53 expression, and pro-apoptotic molecules up-regulation.

A transmembrane domain of Andrias davidianus ranavirus 13R is crucial for co-localization to endoplasmic reticulum and viromatrix

13R, a core gene of Andrias davidianus ranavirus (ADRV), encoded a protein containing a transmembrane domain (TMD) and a restriction endonuclease-like domain. However, the characterization and function of 13R and the protein it encodes remain unclear. In this study, Chinese giant salamander thymus cell (GSTC) was used to investigate the function of 13R. The results showed that the 13R transcripts were detected first at 8 h post-infection (hpi) by RT-PCR and the protein was detected first at 24 hpi by western blot, but the transcription was inhibited by cycloheximide and cytosine arabinofuranoside, indicating that 13R is a viral late gene. Subcellular localization showed that the 13R was co-localized with endoplasmic reticulum (ER) in the cytoplasm, while 13R deleting TMD (13RΔTM) was distributed in cytoplasm and nucleus. During ADRV infection, 13R was observed first in the cytoplasm and nucleus, and later aggregated into the viromatrix, whereas 13RΔTM remain dispersed in cytoplasm and nucleus. Western blot analysis suggested that 13R was a viral non-structural protein and its overexpression did not affect the viral titer in GSTC. All these indicated that the TMD of 13R is crucial for the co-localization into the ER and the viromatrix.

Rapid detection of giant salamander iridovirus by cross-priming amplification

Giant salamander iridovirus (GSIV) belongs to the epizootic genus Ranavirus, and is the cause of epidemic diseases associated with high mortality and great losses to artificial breeding and farming. Here, we established a simple, accurate, and reliable cross-priming amplification (CPA) method to detect GSIV. The CPA assay targets the major caspid protein gene of the GSIV genome to design crossing primer pairs, and the reaction conditions were optimized, including optimal concentrations of the primers, betaine, dNTPs, Mg²⁺, and Bst DNA polymerase, and reaction conditions. The sensitivity was shown to be 10 times higher than that of conventional polymerase chain reaction (PCR), and the specificity was 100%. The results were identified on nucleic acid strips within 3-5 min. Application of the CPA and PCR to 54 samples of giant salamander showed a positive rate of 72.22% and 74.07%, respectively, demonstrating high coincidence (94.44%, kappa = 8.7, P < 0.0001). The sensitivity of the CPA assay was 97.50% and the specificity was 92.86%. Thus, the CPA assay is as effective as conventional PCR, but with added practical advantages of simplicity and an almost instrument-free platform, which will be useful for both laboratories and giant salamander farms.

Low intraspecific variation of Frog virus 3 with evidence for novel FV3-like isolates in central and northwestern Canada

Frog virus 3 (FV3) and FV3-like ranaviruses can infect a variety of cold-blooded aquatic species and present a primary threat to amphibians across the globe. Previous studies of FV3-like viruses have largely investigated higher-level phylogenetic distinctions of these pathogens via portions of the conserved major capsid protein (MCP), and the putative virulence gene vIF-2α. Few studies, however, have investigated the spatial distribution of FV3 variants at the population level3-data that can be used to further understand the spatial epidemiology of this disease. In this study, we sequenced the MCP and vIF-2α of 127 FV3-positive amphibians sampled from Canadian water bodies in Ontario, northeastern Alberta, and southern Northwest Territories to explore whether intraspecific genetic variation exists within FV3. There was a lack of variation at the 2 markers across these regions, suggesting that there is a lack of FV3 sequence diversity in Canada, which may hint at a single source of infection that has spread. However, an undocumented variant termed Wood Buffalo ranavirus (WBRV) was detected in samples from 3 sites in Alberta and Northwest Territories that clustered within the FV3-like lineage with 99.3% sequence homology for MCP. For vIF-2α, all sequences were the expected truncated variant except for 6 samples in Ontario. These latter sequences were suggestive of recombination with common midwife toad virus (CMTV). The lack of variation suggests that higher-resolution genome analyses will be required to further explore the spatial spread and intraspecific variation of the disease.

Common midwife toad ranaviruses replicate first in the oral cavity of smooth newts (Lissotriton vulgaris) and show distinct strain-associated pathogenicity

Ranavirus is the second most common infectious cause of amphibian mortality. These viruses affect caudates, an order in which information regarding Ranavirus pathogenesis is scarce. In the Netherlands, two strains (CMTV-NL I and III) were suspected to possess distinct pathogenicity based on field data. To investigate susceptibility and disease progression in urodeles and determine differences in pathogenicity between strains, 45 adult smooth newts (Lissotriton vulgaris) were challenged via bath exposure with these ranaviruses and their detection in organs and feces followed over time by PCR, immunohistochemistry and in situ hybridization. Ranavirus was first detected at 3 days post infection (p.i.) in the oral cavity and upper respiratory mucosa. At 6 days p.i, virus was found in connective tissues and vasculature of the gastrointestinal tract. Finally, from 9 days p.i onwards there was widespread Ranavirus disease in various organs including skin, kidneys and gonads. Higher pathogenicity of the CMTV-NL I strain was confirmed by higher correlation coefficient of experimental group and mortality of challenged animals. Ranavirus-exposed smooth newts shed virus in feces intermittently and infection was seen in the absence of lesions or clinical signs, indicating that this species can harbor subclinical infections and potentially serve as disease reservoirs.

Pathogen Risk Analysis for Wild Amphibian Populations Following the First Report of a Ranavirus Outbreak in Farmed American Bullfrogs (Lithobates catesbeianus) from Northern Mexico

Ranaviruses are the second deadliest pathogens for amphibian populations throughout the world. Despite their wide distribution in America, these viruses have never been reported in Mexico, the country with the fifth highest amphibian diversity in the world. This paper is the first to address an outbreak of ranavirus in captive American bullfrogs (Lithobates catesbeianus) from Sinaloa, Mexico. The farm experienced high mortality in an undetermined number of juveniles and sub-adult bullfrogs. Affected animals displayed clinical signs and gross lesions such as lethargy, edema, skin ulcers, and hemorrhages consistent with ranavirus infection. The main microscopic lesions included mild renal tubular necrosis and moderate congestion in several organs. Immunohistochemical analyses revealed scant infected hepatocytes and renal tubular epithelial cells. Phylogenetic analysis of five partial ranavirus genes showed that the causative agent clustered within the Frog virus 3 clade. Risk assessment with the Pandora⁺ protocol demonstrated a high risk for the pathogen to affect amphibians from neighboring regions (overall Pandora risk score: 0.619). Given the risk of American bullfrogs escaping and spreading the disease to wild amphibians, efforts should focus on implementing effective containment strategies and surveillance programs for ranavirus at facilities undertaking intensive farming of amphibians.

Fish-specific finTRIM FTR36 triggers IFN pathway and mediates inhibition of viral replication

The tripartite motif (TRIM) family involves many cellular processes, including fundamental functions in antiviral immunity. Antiviral activities of TRIMs are reported in a variety of patterns, and one of the most significant channels is related to the activation of the type-I interferon (IFN) pathway. In this study, we described a fintrim (ftr) gene named ftr36, which is mainly expressed in the gills, skin, and intestines. This study shows that ftr36 encodes a protein affording a potent antiviral effect. In vitro, overexpression of FTR36 mediated an upregulated pattern of recognition receptor retinoic acid-inducible gene I (RIG-I), interferon regulatory factor 3/7(IRF3/7), IFN, and IFN-stimulated genes (ISGs) expression. Thereby, FTR36 expression could afford host defense against the spring viremia of carp virus (SVCV) and the giant salamander iridovirus (GSIV). With the deletion of the RING domain or B30.2 domain separately, the antiviral ability of FTR36 was abolished partially and almost lost its ability to activate the IFN-pathway. These findings indicate that both RING and B30.2 domains are indispensable for the antiviral activity of FTR36. Altogether, this study described a finTRIM FTR36, which can activate IFN-pathways and stimulate ISGs to provide host defense against viral infections.

Rag1 and rag2 gene expressions identify lymphopoietic tissues in juvenile and adult Chinese giant salamander (Andrias davidianus)

Rag1 and rag2 are two closely linked recombination activating genes required for V(D)J recombination of antigen receptors in immature lymphocytes, whose expression can serve as marker to identify the lymphopoietic tissues. To study the development of lymphopoietic tissues in Chinese giant salamander (Andrias davidianus), the Chinese giant salamander rag1 and rag2 coding sequences were cloned and determined. High transcript levels of rag1 and rag2 were co-detected in the thymus before 14 months of age, whereas levels were lower in spleen, liver and kidney at all stage of development. The spatial expression patterns of rag1 and rag2 were studied in combination with igY and tcrβ gene expression using in situ hybridization. Significant transcript signals for rag1, rag2, tcrβ and igY were detected not only in the thymus and spleen but also the liver and kidney of juvenile and adult Chinese giant salamanders, which suggests that cells of lymphocyte lineage are present in multiple tissues of the Chinese giant salamander. This implies that lymphopoiesis may take place in these tissues. The tissue morphology of thymus suggested that the branched thymic primordium developed into mature organ with the development of thymocyte from juvenile to adult. These results not only confirm that as expected the thymus and spleen are primordial lymphopoietic tissues but also suggest that the liver and kidney provide site of lymphocyte differentiation in Chinese giant salamander.

Partial validation of a TaqMan real-time quantitative PCR for the detection of ranaviruses

Ranaviruses are globally emerging pathogens negatively impacting wild and cultured fish, amphibians, and reptiles. Although conventional and diagnostic real-time PCR (qPCR) assays have been developed to detect ranaviruses, these assays often have not been tested against the known diversity of ranaviruses. Here we report the development and partial validation of a TaqMan real-time qPCR assay. The primers and TaqMan probe targeted a conserved region of the major capsid protein (MCP) gene. A series of experiments using a 10-fold dilution series of Frog virus 3 (FV3) MCP plasmid DNA revealed linearity over a range of 7 orders of magnitude (107-101), a mean correlation coefficient (R2) of >0.99, and a mean efficiency of 96%. The coefficient of variation of intra- and inter-assay variability ranged from <0.1-3.5% and from 1.1-2.3%, respectively. The analytical sensitivity was determined to be 10 plasmid copies of FV3 DNA. The qPCR assay detected a panel of 33 different ranaviral isolates originating from fish, amphibian, and reptile hosts from all continents excluding Africa and Antarctica, thereby representing the global diversity of ranaviruses. The assay did not amplify highly divergent ranaviruses, members of other iridovirus genera, or members of the alloherpesvirus genus Cyprinivirus. DNA from fish tissue homogenates previously determined to be positive or negative for the ranavirus Epizootic hematopoietic necrosis virus by virus isolation demonstrated a diagnostic sensitivity of 95% and a diagnostic specificity of 100%. The reported qPCR assay provides an improved expedient diagnostic tool and can be used to elucidate important aspects of ranaviral pathogenesis and epidemiology in clinically and sublinically affected fish, amphibians, and reptiles.

Ranavirus genotypes in the Netherlands and their potential association with virulence in water frogs (Pelophylax spp.)

Ranaviruses are pathogenic viruses for poikilothermic vertebrates worldwide. The identification of a common midwife toad virus (CMTV) associated with massive die-offs in water frogs (Pelophylax spp.) in the Netherlands has increased awareness for emerging viruses in amphibians in the country. Complete genome sequencing of 13 ranavirus isolates collected from ten different sites in the period 2011–2016 revealed three CMTV groups present in distinct geographical areas in the Netherlands. Phylogenetic analysis showed that emerging viruses from the northern part of the Netherlands belonged to CMTV-NL group I. Group II and III viruses were derived from the animals located in the center-east and south of the country, and shared a more recent common ancestor to CMTV-amphibian associated ranaviruses reported in China, Italy, Denmark, and Switzerland. Field monitoring revealed differences in water frog host abundance at sites where distinct ranavirus groups occur; with ranavirus-associated deaths, host counts decreasing progressively, and few juveniles found in the north where CMTV-NL group I occurs but not in the south with CMTV-NL group III. Investigation of tandem repeats of coding genes gave no conclusive information about phylo-geographical clustering, while genetic analysis of the genomes revealed truncations in 17 genes across CMTV-NL groups II and III compared to group I. Further studies are needed to elucidate the contribution of these genes as well as environmental variables to explain the observed differences in host abundance.

Divergent transcriptomic responses underlying the ranaviruses-amphibian interaction processes on interspecies infection of Chinese giant salamander

Background

Ranaviruses (family Iridoviridae, nucleocytoplasmic large DNA viruses) have been reported as promiscuous pathogens of cold-blooded vertebrates. Rana grylio virus (RGV, a ranavirus), from diseased frog Rana grylio with a genome of 105.79 kb and Andrias davidianus ranavirus (ADRV), from diseased Chinese giant salamander (CGS) with a genome of 106.73 kb, contains 99% homologous genes.

Results

To uncover the differences in virus replication and host responses under interspecies infection, we analyzed transcriptomes of CGS challenged with RGV and ADRV in different time points (1d, 7d) for the first time. A total of 128,533 unigenes were obtained from 820,858,128 clean reads. Transcriptome analysis revealed stronger gene expression of RGV than ADRV at 1 d post infection (dpi), which was supported by infection in vitro. RGV replicated faster and had higher titers than ADRV in cultured CGS cell line. RT-qPCR revealed the RGV genes including the immediate early gene (RGV-89R) had higher expression level than that of ADRV at 1 dpi. It further verified the acute infection of RGV in interspecies infection.

The number of differentially expressed genes and enriched pathways from RGV were lower than that from ADRV, which reflected the variant host responses at transcriptional level. No obvious changes of key components in pathway “Antigen processing and presentation” were detected for RGV at 1 dpi. Contrarily, ADRV infection down-regulated the expression levels of MHC I and CD8. The divergent host immune responses revealed the differences between interspecies and natural infection, which may resulted in different fates of the two viruses. Altogether, these results revealed the differences in transcriptome responses among ranavirus interspecies infection of amphibian and new insights in DNA virus-host interactions in interspecies infection.

Conclusion

The DNA virus (RGV) not only expressed self-genes and replicated quickly after entry into host under interspecies infection, but also avoided the over-activation of host responses. The strategy could gain time for the survival of interspecies pathogen, and may provide opportunity for its adaptive evolution and interspecies transmission.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4596-y) contains supplementary material, which is available to authorized users.

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.

Pathogen richness and abundance predict patterns of adaptive major histocompatibility complex variation in insular amphibians

The identification of the factors responsible for genetic variation and differentiation at adaptive loci can provide important insights into the evolutionary process and is crucial for the effective management of threatened species. We studied the impact of environmental viral richness and abundance on functional diversity and differentiation of the MHC class Ia locus in populations of the black-spotted pond frog (Pelophylax nigromaculatus), an IUCN-listed species, on 24 land-bridge islands of the Zhoushan Archipelago and three nearby mainland sites. We found a high proportion of private MHC alleles in mainland and insular populations, corresponding to 32 distinct functional supertypes, and strong positive selection on MHC antigen-binding sites in all populations. Viral pathogen diversity and abundance were reduced at island sites relative to the mainland, and islands housed distinctive viral communities. Standardized MHC diversity at island sites exceeded that found at neutral microsatellites, and the representation of key functional supertypes was positively correlated with the abundance of specific viruses in the environment (Frog virus 3 and Ambystoma tigrinum virus). These results indicate that pathogen-driven diversifying selection can play an important role in maintaining functionally important MHC variation following island isolation, highlighting the importance of considering functionally important genetic variation and host-pathogen associations in conservation planning and management.

Ranavirus phylogenomics: Signatures of recombination and inversions among bullfrog ranaculture isolates

Ranaviruses are emerging pathogens of fish, amphibians, and reptiles that threaten aquatic animal industries and wildlife worldwide. Our objective was to genetically characterize ranaviruses isolated during separate bullfrog Lithobates catesbeianus die-offs that occurred eight years apart on the same North American farm. The earlier outbreak was due to a highly pathogenic strain of common midwife toad virus (CMTV) previously known only from Europe and China. The later outbreak was due to a chimeric ranavirus that displayed a novel genome arrangement and a DNA backbone typical for Frog virus 3 (FV3) strains except for interspersed fragments acquired through recombination with the CMTV isolated earlier. Both bullfrog ranaviruses are more pathogenic than wild-type FV3 suggesting recombination may have resulted in the increased pathogenicity observed in the ranavirus isolated in the later outbreak. Our study underscores the role international trade in farmed bullfrogs may have played in the global dissemination of highly pathogenic ranaviruses.

Vaccination with Recombinant Baculovirus Expressing Ranavirus Major Capsid Protein Induces Protective Immunity in Chinese Giant Salamander, Andrias davidianus

The Chinese giant salamander iridovirus (CGSIV), belonging to the genus Ranavirus in the family Iridoviridae, is the causative agent of an emerging infectious disease causing high mortality of more than 90% and economic losses in Chinese giant salamanders in China. In this study, a recombinant baculovirus-based vaccine expressing the CGSIV major capsid protein (MCP) was developed and its protective immunity in Chinese giant salamanders was evaluated. The recombinant Autographacalifornica nucleopolyhedrosis virus (AcNPV), expressing CGSIV MCP, designated as AcNPV-MCP, was generated with the highest titers of 1 × 10⁸ plaque forming units/mL (PFU/mL) and confirmed by Western blot and indirect immunofluorescence (IIF) assays. Western blot analysis revealed that the expressed MCP reacted with mouse anti-MCP monoclonal antibodies at the band of about 53 kDa. The results of IIF indicated that the MCP was expressed in the infected Spodoptera frugiperda 9 (Sf9) cells with the recombinant baculovirus, and the Chinese giant salamander muscle cells also transduced with the AcNPV-MCP. Immunization with the recombinant baculovirus of AcNPV-MCP elicited robust specific humoral immune responses detected by ELISA and neutralization assays and potent cellular immune responses in Chinese giant salamanders. Importantly, the effective immunization conferred highly protective immunity for Chinese giant salamanders against CGSIV challenge and produced a relative percent of survival rate of 84%. Thus, the recombinant baculovirus expressing CGSIV MCP can induce significant immune responses involving both humoral and cell-mediated immunity in Chinese giant salamanders and might represent a potential baculovirus based vaccine candidate for Chinese giant salamanders against CGSIV.

Mycolactone-producing Mycobacterium marinum infection in captive Hong Kong warty newts and pathological evidence of impaired host immune function

A mass mortality event of captive Hong Kong warty newts Paramesotriton hongkongensis with non-granulomatous necrotic lesions occurred in Taipei Zoo, Taiwan, in 2014. Clinically, the sick newts were lethargic and often covered with water mold Saprolegnia sp. on the skin of the body trunk or extremities. Predominant pathological findings were multifocal non-granulomatous necrotic lesions in the liver, spleen, and kidneys and severe skin infection with Saprolegnia sp., with deep invasion and involvement of underlying muscles. The possibility of ranavirus infection was ruled out by negative PCR results. Unexpectedly, abundant intralesional acid-fast positive bacilli were found in the necrotic lesions of the liver, spleen, and kidney in all 14 sick newts. PCR targeting the hsp65, ITS region, and partial 16S rRNA genes was performed, and the sequence identity from amplified amplicons of hsp65 and partial 16S rRNA genes was 100% identical to that of the corresponding gene fragment of Mycobacterium marinum. Further molecular investigations demonstrated that the current M. marinum was a mycolactone-producing mycobacterium with the presence of esxA/esxB genes. Mycolactone is a plasmid-encoded, immunosuppressive, and cytotoxic toxin. The possible immunosuppression phenomenon characterized by systemic non-granulomatous necrotic lesions caused by M. marinum and the unusual deep invasive infection caused by water mold might be associated with the immunosuppressive effect of mycolactone. Therefore, it should be noted that non-granulomatous necrotic lesions in amphibians can be caused not only by ranavirus infection but also by mycobacteriosis.

Characterization of the Cutaneous Bacterial Communities of Two Giant Salamander Subspecies

Pathogens currently threaten the existence of many amphibian species. In efforts to combat global declines, researchers have characterized the amphibian cutaneous microbiome as a resource for disease management. Characterization of microbial communities has become useful in studying the links between organismal health and the host microbiome. Hellbender salamanders (Cryptobranchus alleganiensis) provide an ideal system to explore the cutaneous microbiome as this species requires extensive conservation management across its range. In addition, the Ozark hellbender subspecies (Cryptobranchus alleganiensis bishopi) exhibits chronic wounds hypothesized to be caused by bacterial infections, whereas the eastern hellbender (Cryptobranchus alleganiensis alleganiensis) does not. We assessed the cutaneous bacterial microbiome of both subspecies at two locations in the state of Missouri, USA. Through 16S rRNA gene-based amplicon sequencing, we detected more than 1000 distinct operational taxonomic units (OTUs) in the cutaneous and environmental bacterial microbiome. Phylogenetic and abundance-based dissimilarity matrices identified differences in the bacterial communities between the two subspecies, but only the abundance-based dissimilarity matrix identified differences between wounds and healthy skin on Ozark hellbenders. The higher abundance of OTUs on Ozark wounds suggests that commensal bacteria present on the skin and environment may be opportunistically colonizing the wounds. This brief exploration of the hellbender cutaneous bacterial microbiome provides foundational support for future studies seeking to understand the hellbender cutaneous bacterial microbiome and the role of the bacterial microbiota on chronic wounds of Ozark hellbenders.

In vitro antiviral efficacy of moroxydine hydrochloride and ribavirin against grass carp reovirus and giant salamander iridovirus

Moroxydine hydrochloride (Mor) and ribavirin (Rib) have been reported to exhibit multi-antiviral activities against DNA and RNA viruses, but their antiviral activities and pharmacologies have seldom been studied in aquaculture. This paper has selected 3 aquatic viruses including a double-stranded RNA virus (grass carp reovirus, GCRV), a single-stranded RNA virus (spring viraemia of carp virus, SVCV) and a DNA virus (giant salamander iridovirus, GSIV) for antiviral testing. The results showed that Mor and Rib can effectively control the infection of GCRV and GSIV in respective host cells. Further study was undertaken to explore the antivirus efficiencies and pharmacological mechanisms of Mor and Rib on GCRV and GSIV in vitro. Briefly, compounds showed over 50% protective effects at 15.9 µg ml-1 except for the group of GSIV-infected epithelioma papulosum cyprinid (EPC) cells treated with Mor. Moreover, Mor and Rib blocked the virus-induced cytopathic effects and apoptosis in host cells to keep the normal cellular structure. The expression of VP1 (GCRV) and major capsid protein (MCP; GSIV) gene was also significantly inhibited in the virus-infected cells when treated with Mor and Rib. Cytotoxicity assay verified the 2 compounds had no toxic effects on grass carp ovary (GCO) cells and EPC cells at ≤96 µg ml-1. In conclusion, these results indicated that exposing GCRV-infected GCO cells and GSIV-infected EPC cells to Mor and Rib could elicit significant antiviral responses, and the 2 compounds have been shown to be promising agents for viral control in the aquaculture industry.

Reconstructing the emergence of a lethal infectious disease of wildlife supports a key role for spread through translocations by humans

There have been few reconstructions of wildlife disease emergences, despite their extensive impact on biodiversity and human health. This is in large part attributable to the lack of structured and robust spatio-temporal datasets. We overcame logistical problems of obtaining suitable information by using data from a citizen science project and formulating spatio-temporal models of the spread of a wildlife pathogen (genus Ranavirus, infecting amphibians). We evaluated three main hypotheses for the rapid increase in disease reports in the UK: that outbreaks were being reported more frequently, that climate change had altered the interaction between hosts and a previously widespread pathogen, and that disease was emerging due to spatial spread of a novel pathogen. Our analysis characterized localized spread from nearby ponds, consistent with amphibian dispersal, but also revealed a highly significant trend for elevated rates of additional outbreaks in localities with higher human population density—pointing to human activities in also spreading the virus. Phylogenetic analyses of pathogen genomes support the inference of at least two independent introductions into the UK. Together these results point strongly to humans repeatedly translocating ranaviruses into the UK from other countries and between UK ponds, and therefore suggest potential control measures.

Autophagy and apoptosis induced by Chinese giant salamander (Andrias davidianus) iridovirus (CGSIV)

The outbreak of Chinese Giant Salamander (Andrias davidianus, CGS) Iridovirus (CGSIV) caused massive death of CGSs. However, some CGSs with low level of CGSIV usually could survive. In our study, major capsid protein (MCP) DNA replicates of CGSIV in shedding skin were employed to assess the relative content of CGSIV in the living CGSs by qPCR. Furthermore, the examinations of autophagy and apoptosis in CGSs in vivo and in the primary renal cells in vitro were performed, respectively. The results showed that the relative contents of CGSIV in the shedding skin could reflect those in liver, spleen, and kidney of the CGSs. In these tissues of the CGSs with low-level replicates of CGSIV, there were not obviously macroscopic lesions. But the irregularly-shaped vesicles perhaps involving in autophagosome were observed by transmission electron microscopy (TEM). The LC3B protein displayed uneven distribution by Immunohistochemistry and the level mRNA of Atg5 was higher in these tissues than that in the tissues of healthy CGSs using qRT-PCR. Meanwhile, the apoptosis also appeared in these tissues by TUNEL staining and higher level mRNA of type I IFN were detected in these tissues using qRT-PCR. Further, both the expression level of LC3B II protein and Atg5 mRNA increased significantly at 2h after the virus infected the primary renal cells from the health CGSs in vitro. In addition, apoptosis and type I IFN mRNA began to increase significantly at 4h after the virus infected the renal cells. It was suggested that autophagy may be a pivotal role for survival of CGSIV in the CGSs during early infection and the rapid proliferation of CGSIV could be inhibited by innate immune response and apoptosis.

Reintroduction and Post-Release Survival of a Living Fossil: The Chinese Giant Salamander

Captive rearing and reintroduction / translocation are increasingly used as tools to supplement wild populations of threatened species. Reintroducing captive-reared Chinese giant salamanders may help to augment the declining wild populations and conserve this critically endangered amphibian. We released 31 captive-reared juvenile giant salamanders implanted with VHF radio transmitters at the Heihe River (n = 15) and the Donghe River (n = 16) in the Qinling Mountains of central China. Salamanders were monitored every day for survival from April 28^th 2013 to September 3^rd 2014. We attempted to recapture all living individuals by the end of the study, measured their body mass and total body length, and checked for abnormalities and presence of external parasites. Two salamanders at the Heihe River and 10 animals at the Donghe River survived through the project timeline. Nine salamanders were confirmed dead, while the status of the other 10 animals was undetermined. The annual survival rate of giant salamanders at the Donghe River (0.702) was 1.7-fold higher than that at the Heihe River (0.405). Survival increased as individuals were held longer following surgery, whereas body mass did not have a significant impact on survival rate. All salamanders recaptured from the Donghe River (n = 8) increased in mass (0.50 ± 0.13 kg) and length (5.5 ± 1.5 cm) after approximately 11 months in the wild, and they were only 7% lighter than wild animals of the same length (mean residual = -0.033 ± 0.025). Our results indicate that captive-reared Chinese giant salamanders can survive in the wild one year after release and adequate surgical recovery time is extremely important to post-release survival. Future projects may reintroduce older juveniles to achieve better survival and longer monitoring duration.