Abstract
Previous studies of coral viruses have employed either microscopy or metagenomics, but few have attempted to comprehensively link the presence of a virus-like particle (VLP) to a genomic sequence. We conducted transmission electron microscopy imaging and virome analysis in tandem to characterize the most conspicuous viral types found within the dominant Pacific reef-building coral genus Acropora. Collections for this study inadvertently captured what we interpret as a natural outbreak of viral infection driven by aerial exposure of the reef flat coincident with heavy rainfall and concomitant mass bleaching. All experimental corals in this study had high titers of viral particles. Three of the dominant VLPs identified were observed in all tissue layers and budding out from the epidermis, including viruses that were ∼70, ∼120, and ∼150 nm in diameter; these VLPs all contained electron dense cores. These morphological traits are reminiscent of retroviruses, herpesviruses, and nucleocytoplasmic large DNA viruses (NCLDVs), respectively. Some 300–500 nm megavirus-like VLPs also were observed within and associated with dinoflagellate algal endosymbiont (Symbiodinium) cells. Abundant sequence similarities to a gammaretrovirus, herpesviruses, and members of the NCLDVs, based on a virome generated from five Acropora aspera colonies, corroborated these morphology-based identifications. Additionally sequence similarities to two diagnostic genes, a MutS and (based on re-annotation of sequences from another study) a DNA polymerase B gene, most closely resembled Pyramimonas orientalis virus, demonstrating the association of a cosmopolitan megavirus with Symbiodinium. We also identified several other virus-like particles in host tissues, along with sequences phylogenetically similar to circoviruses, phages, and filamentous viruses. This study suggests that viral outbreaks may be a common but previously undocumented component of natural bleaching events, particularly following repeated episodes of multiple environmental stressors.
Introduction
Viruses (phages and eukaryotic viruses) are abundant and diverse residents of stony coral colonies (reviewed in ). These viruses likely play multiple, parasitic and commensal roles in the health of coral reefs (e.g., ; ; ; ; ). Research interrogating the impact of viruses on colony fitness and survival under different environmental contexts is of particular importance, given anthropogenic climate forcing and other impacts (). For example, abiotic conditions that stress coral colonies, such as elevated seawater temperatures or UV exposure, may trigger viral infections that contribute to coral bleaching and disease (; ; ). Identifying potential mechanisms of coral reef decline is increasingly important given accelerations in this process during recent decades (e.g., ; ), and the current global mass bleaching event1.
Although the field of coral virology remains in its infancy, several groups have applied microscopy or genomics to examine the diversity and roles of viruses in coral holobionts. Microscopy studies have presented evidence that virus-like particles (VLPs) are present in all tissue layers of apparently healthy and diseased corals: the gastrodermis, mesoglea, and epidermis, as well as in the coral surface microlayer (CSM; e.g., ; ; ; ; ). The physical structure of VLPs also has been examined within cultured Symbiodinium (; ; ). Some of these observed VLPs likely represent particles produced during the lytic replication phase of previously latent or endogenous infections of the coral animal, its dinoflagellate algae, or its microbiota (; ). , for example, were able to distinguish 17 sub-groups of VLPs associated with the CSM of four species of Australian corals based on morphological similarities. The role of each of these groups of viruses is uncertain, especially considering that the density of some VLPs within the CSM is relatively low. However, in other cases, transmission electron microscopy has revealed structures within corals that are highly indicative of massive viral infection (e.g., crystalline arrays, viral factories; ).
Yet standalone transmission electron microscopy (TEM) images can pose interpretive challenges. A set of TEM images may contain VLPs that present only some of the diagnostic morphological characteristics of a viral group, or characteristics that appear representative of many described viral groups (e.g., Figure 3 in ). Further, since viral Families can encompass a range of capsid sizes and shapes and may overlap in these characteristics, microscopy-based studies may not fully resolve a group of VLPs. VLPs reminiscent of a large group of phylogenetically related viruses, the nucleocytoplasmic large DNA viruses (NCLDVs), exemplify this issue (e.g., ). The NCLDV group includes the Phycodnaviridae, Iridoviridae, Poxviridae, Mimiviridae, and Ascoviridae, as well as the recently described giant viruses, marseillevirus and lausannevirus (; ; Yutin and Koonin, 2012). VLPs that are within the cytoplasm, larger than 120 nm, and generally icosahedral in shape are often interpreted as NCLDV-like but several exceptions to this rule remain, such as the poxviruses, and pandoravirus, which are NCLDVs that exhibit very different morphological characteristics.
Genomic and proteomic-based studies have identified patterns in the diversity and abundance of genomic sequences similar to described viruses within healthy and diseased tropical corals (; ; ; ; ; ) and cultured Symbiodinium (; ; ; ; ; ; ,), as well as in cold water corals (; ). For example, a strong correlation between specific viral markers in bleached, diseased, and healthy Orbicella corals was used to establish a role for small circular ssDNA viruses (SCSDVs) in white plague disease (). Although this work was somewhat substantiated by a TEM study on another coral white disease (), viral metagenome studies often contain sequence similarities to many viral groups, most of which have not been corroborated by microscopy-based studies. For example, numerous studies have found sequences similar to mimiviruses and baculoviruses, and yet no TEM study has confirmed these annotations (; ; ). This is likely due to several issues. A single sequence read or contig may have significant similarity to multiple viral groups because many related viruses share some gene homology. Alternatively, reads or contigs may have only a few sequence similarities with relatively high associated e-values. It can also be difficult to rule out contamination as a source of error for sequence similarities represented by few reads within a metagenome, or for reads similar to cosmopolitan or host-associated viral remnants (e.g., retrotransposons, retroelements). Thus, ambiguity remains when metagenomics is the sole approach applied to characterize the viral consortia associated with corals.
A comparative analysis of several metagenome studies recently addressed some of these challenges and cataloged a cosmopolitan set of viruses in corals and their symbionts (). This meta-analysis showed that, based on presence in >90% of 35 surveyed metagenomes, coral holobionts contain signatures of nine major Families in the dsDNA Group I viral lineages. These families include all of the major Caudovirales (Sipho-, Podo-, and Myoviridae) and many eukaryotic viruses. Within the eukaryotic viruses, the Herpesvirales Order, as well as five Families (Phycodnaviridae, Iridoviridae, Poxviridae, Mimiviridae, and Ascoviridae) within the NCLDVs are members of this 90% carriage cosmopolitan virome. Genomic signatures from some ssDNA (e.g., Circoviridae) and ssRNA (e.g., Retroviridae and Caulimoviridae) viral lineages are also well represented, but fall below a 90% threshold, likely due to biological and technical differences among studies (for discussion, see ).
Regardless of the approach applied, a body of evidence indicates that herpes-like viruses and one or more NCLDVs associate with coral holobionts. Yet neither microscopy nor metagenomics alone has fully resolved the identity of either viral group within corals. Thus, this study sought finer taxonomic resolution for one or more groups by characterizing the viruses associated with fragments of dominant reef-building Pacific acroporid corals using morphological and sequence-based approaches in tandem. The aims of this work were to: (1) identify and compare specific VLPs in coral tissues; and (2) improve delineation of the core coral virome through the use of visual descriptors of viral taxonomy in conjunction with metagenomics analysis.
Materials and Methods
Overview of Environmental Setting, Experimental Setup, and Design
The experiments reported here were conducted using acroporid colonies (N = 5 for A. aspera; 4 for A. millepora) collected from the tidal reef flat off of Heron Island, Queensland, Australia (23°26′39.63′′S, 151°54′46.70′′E, Figure 1) in March of 2011. To evaluate the ambient environmental conditions at the study site, we used the Integrated Marine Observing System2 run by the Australian Institute of Marine Science to characterize the following parameters: rain intensity and timing, air temperature, water temperature (on the reef flat and at ∼8 m depth on the reef slope), and water height variations. Prior to and during our study, this reef flat experienced a period of low tides that caused repeated aerial exposure of reef flat colonies (e.g., Figure 1 photograph) and increased residence time of water on the reef flat. Heavy rainfall was coincident with some low tides and aerial exposure of colonies. The temperature range experienced on the reef flat (measured at 1.1 m depth) was neither unique for the season nor as extreme as previous months, however (Figure 2, Environmental Setting in the Supplementary Material).
FIGURE 1
FIGURE 2
Corals were collected for the A. aspera experiment on March 18, 2011 (Figure 2). The experiment was initiated on March 19, 2011 and ran for 6 days. Corals were collected for the A. millepora experiment on March 20, 2011. The A. millepora experiment was initiated on March 22, 2011 and ran for ∼4.75 days. We hypothesized that aerial exposure and rainfall stressors prior to collection might have triggered bleaching in our acroporid colonies, and primed them for viral production. Therefore, we placed our coral colonies in flow-through (∼3 L min–1) seawater tanks and characterized their health states for at least 24 h prior to initiating the experimental periods. All colonies collected for use in this study exhibited normal pigmentation and were “apparently healthy” from the time of collection through the start of the experimental treatments (Supplementary Figure S1A). On March 20, 2011 and through the remainder of the experimental period, we observed a coral bleaching event on the reef (Figure 1 photograph). Figure 2 integrates details of the environmental setting of this study with observations of in situ bleaching on the reef flat and our experimental design.
Once acclimated to flow-through tanks, coral branches were subjected to a variety of experimental injection treatments (e.g., viral inoculate, heat-killed viral inoculate, not injected) and some A. aspera were additionally exposed to a thermal stress treatment (∼2°C above ambient) and/or virus-free seawater inoculation; these treatments are described in the Supplementary Material. Following this, every 24 h, all coral branches were photographed and monitored for signs of stress and/or disease. The visual appearance of each coral branch, as well as the presence/absence of mucus, bleaching, and lesions were recorded. On March 25, 2011, 40 A. aspera specimens were sacrificed for TEM and five control saline-injected A. aspera samples were frozen at –80°C for virome analysis. Different fragments were used for microscopic and genomic analysis, but these fragments were from the same parent colonies and experienced identical treatments. On March 27, 2011, 12 A. millepora specimens were sacrificed for TEM.
Transmission Electron Microscopy
Approximately 1 to 5 mm3 per coral specimen was immersed in a TEM fixative (2% EM-grade glutaraldehyde in virus-free 3x PBS) and stored in the dark at 4°C until processing. Samples were processed following the methods of . In brief, decalcified coral samples (20% ETDA) were washed in 0.2 M cacodylate buffer (pH 7.0) and post-fixed with 1% osmium tetroxide in 0.1 M cacodylate buffer (pH 7.0). The samples were subsequently washed in distilled water, dehydrated through immersion in a series of ethanol and propylene oxide baths, and embedded in resin. Ultra-thin sections were generated using a Leica ultracut UC6 microtome and diamond knife. Sections were stained on 200 μm copper grids for 3 minutes with 5% aqueous uranyl acetate followed by a 1-min stain with lead citrate. Replicated tissue sections were collected from each resin-blocked sample until all coral tissue regions and layers could be assessed. Multiple sections from each sample were visualized; each sample was viewed for equal time (approximately 1 h) on a JEM-1010 Transmission Electron Microscope at the University of Queensland Centre for Microscopy and Microanalysis.
Virome Generation and Analysis
A single virome was constructed from viral DNA from five control saline-injected A. aspera specimens to corroborate TEM data. All fragments appeared healthy (i.e., did not exhibit evidence of paling or tissue sloughing) when frozen at the end of the experiment (Supplementary Figure S1B). With some minor modifications, we used our standard protocol for isolating viral particles from coral tissues (). Briefly, each coral specimen was defrosted and tissue removed with an airbrush containing virus-free 1x PBS (pH 7.4). This combined slurry was centrifuged at 3220 × g for 20 min at room temperature. Samples were decanted and the supernatant passaged through a 0.8 μm nucleopore impact filter from Whatman. Filtrate was placed on a CsCl density gradient containing four densities (1.7, 1.5, 1.35, 1.2 g ml–1) and spun at 22,000 rpm for 2 h at 4°C. Abundant viral particles were identified from the 1.2 g ml–1 CsCl density layer () and recovered. This viral isolate was then filtered through a 0.45 μm Sterivex. Particle DNA was extracted with a formamide procedure () and then amplified using the phi29 polymerase multiple displacement method (; Methodological Considerations provided in Supplementary Materials). The Genomiphi kit V2 (GE Healthcare Life Sciences) was used according the manufacturer’s recommended procedure.
Approximately 16.5 ng of this combined A. aspera control saline-injected viral DNA was prepared for sequencing using the Nextera XT kit. Illumina MiSeq 150 bp paired-end sequencing generated reads of approximately 300 bp. Sequences were quality filtered (phred = 30) and trimmed. High quality paired-end reads were merged using PEAR (Zhang et al., 2014). Merged and singleton sequences were then combined into a single file for further analysis. Sequences were screened for host and human contamination by using BLASTn (e-value ≤ 10–20) against the Acropora digitifera and human genomes, respectively. Sequences were further filtered with BLASTn (e-value ≤ 10–20) to the entire Refseq NCBI database to remove any reads that annotated as potential cellular organisms (i.e., sequence contaminants). Reads with BLASTn annotations to viruses were then assembled with Velvet using a kmer size of 71 (Zerbino and Birney, 2008). A tBLASTx analysis of the contigs was then performed against the NCBI RefSeq viral database (e-value ≤ 10–7). Viral taxonomy was assigned using NCBI’s taxonomy tree and in-house python scripts. Similarities were then parsed at the viral Family level. Reads were archived at the European Nucleotide Archive (ENA; Accession #PRJEB12107).
Results
At the end of the experimental period, some coral fragments (particularly those exposed to heat treatments) exhibited visible signs of stress including paling, bleaching, and tissue sloughing. These signs generally coincided with the onset of an in situ bleaching event on Heron Island. No potential effects from the control saline injection were evident in any A. aspera samples during the experimental period based on daily photographs of each coral fragment and visual inspection of fragments at regular intervals. Regardless of our experimental challenges, all A. aspera and A. millepora fragments (including non-injected controls) showed microscopic evidence consistent with a massive viral infection. Given this, we chose to generate a single virome from the control saline-injected A. aspera samples (N = 5 fragments) and we here interpret all TEM data jointly (independent of challenge). Based on these congruent morphological and genomic data, we show that this outbreak consisted of four major viruses: an atypical herpes-like virus, a retrovirus similar to gamma-retroviruses, and two NCLDVs: one 150–180 nm VLP most similar to phycodnaviruses and associated with the host coral and another ∼300–500 nm NCLDV in the candidate Family Megaviridae and associated with resident Symbiodinium.
Virome Analysis Reveals Dominance by Diverse Eukaryotic Viruses and Phages
Using viral particle purification and the Illumina MiSeq platform, a 5,069,340 sequence virome library was generated that had a mean sequence length of 301 bp, of which 829,330 (16%) passed quality control and pre-screening for similarities to non-viral targets. Velvet assembly resulted in a total of 70,807 contigs, with a mean length of 1,160 bp and a maximum of 3,070 bp. Of these, 18,432 contigs (26%) were highly similar (i.e., e-values ≤ 10–7) to a completed viral genome. These contigs contained ∼42,000 viral gene annotations, of which 2,587 were unique. Of these, 841 unique contigs fell within viral Families.
To determine the types of viruses present in Acropora aspera, unique contig similarities were binned hierarchically based on annotations into 19 viral Families (Figure 3). Overall, these contigs were similar to phages and three groups of eukaryotic viruses (dsDNA, ssDNA, and RNA genomes). A majority of the annotations fell within the Group I classification of dsDNA eukaryotic viruses (Figure 3, red bars), retroviruses (light blue bars) and phages (green and dark blue bars). The five dominant families were: Siphoviridae, Myoviridae, Retroviridae, Herpesviridae, and Phycodnaviridae, in decreasing order of relative abundance. Family-level analyses based on the (1) top five similarities to a given contig, (2) top read hits, and (3) top five similarities to a given read produced results concordant with the top contig results described here (Supplementary Tables S1A,B). Three ssDNA viral Families, the Circoviridae, Inoviridae, and Microviridae, produced many more similarities to reads than to contigs (see Methodological Considerations in Supplementary Material, Supplementary Table S1A).
FIGURE 3
Herpesvirus-Like Viral Particles are Abundant in Acropora aspera and Acropora millepora Corals
Based on TEM evidence, the most commonly identified VLP in both coral species was composed of an enveloped and circular capsid ranging from 120 to 150 nm in diameter and containing an electron dense core (Figure 4), a morphology highly reminiscent of herpesviruses. These herpes-like VLPs were present individually and as large clusters (Figures 4A,B,F) that ranged from ∼5 to 40 VLPs within host coral epidermal and gastrodermal cells. Acropora aspera contained most of the large inclusions of this VLP. In some instances, clusters of these VLPs were found in cellular vacuoles (Figures 4F–H) similar to those commonly found in herpes infections (). In vacuoles, herpes-like VLPs often were associated with other VLPs (Figures 4G,H) and constituted the dominant structure within host cells (Figure 4F). Herpes-like VLPs were not observed in host cell nuclei.
FIGURE 4
Results from the A. aspera virome were consistent with the observation of herpes-like VLPs in TEM images. The third most abundant unique eukaryotic virus similarities (72 unique contig and 816 read annotations) were to the Herpesvirales Order (Herpesviridae and Alloherpesviridae; Figure 3). Similarities to several important functional genes were characterized, including: a uracil DNA glycosylase and a DNA polymerase (Table 1). Further, when considering all contig annotations (not just unique ones), a total of 15,083 contigs were similar to a single betaherpesvirinae genome, Human herpesvirus 6A. Phylogenetic analysis of a DNA polymerase-like contig generated in this study indicates that it originates from an undescribed virus within the Herpesviridae that is most similar to mammalian gammaherpesviruses (methods and detailed results provided in the Supplementary Material, Supplementary Figure S2).
Table 1
| Viral family annotation | ENA Accession No. | Contig coverage | nt Size | E-Value | aa % Identity | Viral genome similarity | Gene annotation |
|---|---|---|---|---|---|---|---|
| Herpesviridae | LT009377 | 6.43 | 511 | 2e – 44 | 52% | Elephant endotheliotropic herpesvirus 6 | Uracil DNA glycosylase |
| Herpesviridae | LT009378 | 1.18 | 497 | 2e – 17 | 34% | Gorilla lymphocryptovirus 2 | DNA polymerase |
| NCLDV/Mega | LT009379 | 2.5 | 635 | 2e – 11 | 28% | Pyramimonas orientalis virus | MutS protein |
| NCLDV/Mega | LT009380 | 1.34 | 321 | 4e – 53 | 64% | Lausannevirus | Eukaryotic peptide chain release factor subunit 1 |
| NCLDV/Mega | LT009381 | 2 | 907 | 2e – 08 | 30% | Acanthamoeba polyphaga mimivirus | Putative serine/threonine-protein kinase/receptor |
| NCLDV/Mega | LT009382 | 5.1 | 896 | 1e – 34 | 36% | Megavirus chilensis | UDP-N-acetylglucosamine 2-epimerase |
| NCLDV/Mega | LT009383 | 1.94 | 736 | 2e – 10 | 49% | Cafeteria roenbergensis virus BV-PW1 | Superfamily II helicase/eIF-4AIII |
| NCLDV/Mega | LT009384 | 2.45 | 751 | 2e – 22 | 45% | Cafeteria roenbergensis virus BV-PW1 | Putative photolyase |
| NCLDV/Mega | LT009385 | 1.58 | 259 | 7e – 21 | 67% | Cafeteria roenbergensis virus BV-PW1 | Putative DnaK/Hsp70 |
| NCLDV/Phyco | LT009386 | 2.37 | 1209 | 3e – 13 | 31% | Bathycoccus sp. RCC1105 virus | UDP-N-acetylglucosamine O-acyltransferase |
| NCLDV/Phyco | LT009387 | 5.75 | 720 | 7e – 09 | 47% | Paramecium bursaria Chlorella virus | ATP-dependent helicase |
| NCLDV/Pox | LT009388 | 1.76 | 880 | 3e – 09 | 31% | Pigeonpox virus | Hypothetical protein fep_013 |
| Retroviridae | LT009389 | 1.1 | 666 | 3e – 07 | 26% | Walleye epidermal hyperplasia virus 1 | RT_ZFREV_like reverse transcriptases |
| Retroviridae | LT009390 | 2.38 | 474 | 6e – 19 | 40% | Walleye epidermal hyperplasia virus 2 | Polymerase protein |
| Retroviridae | LT009391 | 4.7 | 856 | 5e – 15 | 32% | Reticuloendotheliosis virus | gag protein |
| Retroviridae | LT009392 | 3.2 | 474 | 6e – 13 | 46% | Reticuloendotheliosis virus | Envelope glycoprotein |
| Retroviridae | LT009393 | 4.86 | 461 | 2e – 14 | 48% | Reticuloendotheliosis virus | Protease |
| Retroviridae | LT009394 | 2.18 | 466 | 2e – 11 | 28% | Simian foamy virus | Integrase core domain |
Examples of gene annotations from a single control Acropora aspera virome (generated from five coral fragments), including those from Herpesviridae, the Nucleocytoplasmic Large DNA Viruses (NCLDVs), and Retroviridae.
Annotations are based on tBLASTx to the NCBI viral database (e-value threshold of <10–7). Contig sequences can be obtained through the European Nucleotide Archive (ENA).
NCLDV-Like Viruses in Host Tissues
The next most common VLP type observed in both A. aspera and A. millepora epidermal cells fell within the NCLDVs (Figure 5). In TEM images where this type was observed, 10 to 59 VLPs were typically present. This VLP morphology was icosahedral, electron-dense and enveloped. However, these VLPs were consistently larger in diameter (150–180 nm) than the 120–150 nm herpes-like virus. Very dark membranes, a more angular appearance, and a wider space between the electron-dense core and the envelope membrane (Figures 5B,I) also distinguished these VLPs, relative to the herpes-like viruses.
FIGURE 5
Within the A. aspera virome, the second most abundant group of eukaryotic virus annotations was to the NCLDVs (79 unique contigs and 733 read annotations). These unique similarities fell among the Families that make up the candidate NCLDV cluster: Phycodnaviridae (n = 41 best contigs), Mimiviridae (n = 17), Poxviridae (n = 13), Iridoviridae (n = 5), Marseillevirus (n = 2), and Ascoviridae (n = 1). Importantly, these annotations contain multiple phylogenetically informative protein encoding genes including: an ATP-dependent helicase, a UDP-N-acetylglucosamine O-acyltransferase, a peptide release factor, a DnaK/Hsp70, and a putative photolyase (Table 1). When analyzed based on all contig annotations (not just unique ones) the most abundant similarities to this classification were to four genomes within the candidate Family Megaviridae and one Phycodnaviridae genome: Acanthamoeba polyphaga mimivirus (n = 50), Megavirus chilensis (n = 28), Cafeteria roenbergensis virus (n = 22), Phaeocystis globosa virus (n = 14), and Paramecium bursaria Chlorella virus (n = 50), respectively.
Symbiodinium-Associated Megaviruses
In addition to the ∼150 nm NCLDV described above, another NCLDV-like icosahedral VLP was observed within and adjacent to probable symbiosomes (e.g., Figure 6C). These Symbiodinium-associated VLPs had mean diameters from 300 to 500 nm but varied in capsid morphology with some being rounded (Figures 6A,B) and others polyhedral (Figure 6C). VLPs in this size range have been referred to as megaviruses, and several of the NCLDV viral Families (e.g., Mimiviridae, Marseillevirus) fall within this size category (; ). Given their proximity to or location within Symbiodinium cells, these putative viruses are a potentially distinct nucleocytoplasmic large DNA viral type perhaps within Phycodnaviridae or candidate Family Megaviridae. For example, although numerous gene annotations were found across these Families, one phylogenetically informative sequence, MutS, indicated that this putative virus is most similar to Pyramimonas orientalis virus (Table 1), a genome recently reclassified as a megavirus (; ). Similarly, phylogenetic analysis of this contig indicated that it originated from an undescribed virus within the Megaviridae (for methods and detailed results, see Supplementary Material; Supplementary Figure S3).
FIGURE 6
Morphological and Metagenomic Evidence of Additional Acroporid Virus Diversity
Additional VLP morphologies were identified in all corals, but appeared to be present at lower prevalence than the herpes-like and NCLDV-like VLP types described above. These other putative viral types included small (∼17–40 nm) star/pentagonal shaped viruses similar to astro-, circo-, parvo-, or nano-like virus particles (Figures 7A,B), amorphous or egg-shaped retrovirus-like VLPs (∼70 nm in diameter; Figures 7C,D), filamentous VLPs within large viral factories or inclusion bodies (Figures 7E,F), and phage-like VLPs within bacterial aggregates enclosed in peri-algal spaces of the host gastroderm (Figure 7G).
FIGURE 7
In the A. aspera virome, phage sequences were dominated by the Sipho- and Myoviridae (Figure 3) with the top two most abundant contig annotations being similar to a Streptococcus myophage (n = 616 total) and a Streptococcus siphophage (n = 298). These large numbers of sequences are notable because the next most abundant group of phage annotations was to a Pelagibacter phage, which only had 28 total similarities across all contigs. Sequences with similarity to the Retroviridae were highly abundant in the A. aspera virome; Circoviridae were also notable (<5%; Figure 3).
Discussion
Using a modification and combination of established methods for coral viral metagenomics and TEM, we found highly congruent molecular (1000s of high quality sequences and contigs) and physical (1000s of VLPs) evidence to support the hypothesis that although a diverse group of viruses are present in acroporid corals, four main or core eukaryotic viruses dominate these species: two NCLDVs including a new megavirus, a somewhat atypical herpes-like virus, and a gamma-retrovirus. It should be noted that genomic support for these four core viral groups are based on metagenomic analyses of control saline-injected fragments of A. aspera that showed a high abundance of diverse VLPs in the TEM images. Although sequencing of additional samples may have produced some novel evidence of viral diversity in terms of similarities to diagnostic viral genes or the construction of longer contigs, it is unlikely that it would have profoundly influenced the core viral groups identified from these two acroporid species. We infer this based on the fact that TEM images generated from all experimental corals were examined, and all major VLP morphologies from the total dataset were linked to predominant viral sequence similarities in the metagenome. Thus, to a certain extent, joint application of microscopy and metagenomics (even to a single treatment) improved our confidence in the identification of the dominant or core viral types within an environmental consortium.
Multiple Nucleocytoplasmic Large DNA Viruses Infect Coral Holobionts
Nucleocytoplasmic large DNA virus sequences have been found in every coral virus metagenomic study undertaken thus far () and often are the most abundant similarities found in corals and in cultures of their dinoflagellate algal symbionts (; ; ). Two distinct variants reminiscent of NCLDVs were commonly identified in this study. One variant was a ∼150 nm NCLDV located within host tissues (Figure 5); the other variant was significantly larger at ∼300 nm and located near or within Symbiodinium (Figure 6). Congruently, NCLDV contig similarities were highly abundant in the virome we generated. Work on NCLDVs is rapidly advancing, and new viruses are now found routinely (e.g., pithovirus, ; faustovirus, ). NCLDVs span a large range of particle sizes (∼140–1,100 nm diameter capsids) and genome lengths (∼100–2.5 Mb; ; ). Thus ascribing gene sequences within the systematic framework of this group of viruses remains problematic due to the vagaries of viral morphologies and genome sequences. However, the microscopic and Family-level genomic data generated here (Figure 3), in combination with a review of the literature and an in-depth phylogenetic analysis of sequences from this virome and other previously annotated viromes (e.g., ), provide us with sufficient data to clearly delineate the taxonomic identity of the ∼150 nm versus the ∼300 nm coral-associated NCLDVs.
The identity of the larger (∼300 nm) NCLDV characterized here is almost certainly a new relative of the megaviruses. The most common annotations in the dataset were to members of this candidate Family including: Acanthamoeba polyphaga mimivirus, Megavirus chilensis, Cafeteria roenbergensis virus, and Phaeocystis globosa virus strain 16T virus genomes. Importantly, one of the genes identified in our virome was a MutS homolog (Table 1); this gene has become diagnostic for megaviruses and has led to a reorganization of this clade of the NCLDVs (; ). In our unrooted tree, this A. aspera MutS contig is placed as sister taxa to a clade containing two megaviruses and a phycodnavirus with strong bootstrap support (Supplementary Figure S3). Further, re-annotation of our previously identified viral mRNA sequences (GenBank # JX026066.1; ) from both a coral and a Symbiodinium culture showed that another diagnostic gene for this candidate Family, DNA Polymerase B, also annotates as Pyramimonas orientalis virus (BLASTx 66% identity, e-value ≤ 8e–58). Given that gene similarities to these viruses have now been found in corals from the Atlantic and Pacific, as well as in Symbiodinium cultures, they likely represent a cosmopolitan Symbiodinium-infecting megavirus.
Given the physical size of these VLPs imaged in and adjacent to Symbiodinium (Figure 6) and the identified MutS gene similarity in our A. aspera virome, we hypothesize that coral holobionts harbor megaviruses most likely related to the Pyramimonas orientalis virus () or perhaps Cafeteria roenbergensis virus (). However, megaviruses vary in their physical structure both in size and shape (e.g., the presence/absence of projections; ). Interestingly, Cafeteria roenbergensis virus, like many mimiviruses, contains a unique star-like structure at one end of the virion (), whereas Pyramimonas orientalis virus does not (). As mentioned above, after re-analysis of viral mRNA sequences from our previous work on NCLDVs in corals and Symbiodinium, we found that another diagnostic gene for this candidate Family (), DNA polymerase B, also annotates to Pyramimonas orientalis virus. Therefore, given that: (1) our NCLDV lacks a star-shaped structure and is genetically similar to the Pyramimonas orientalis clade of megaviruses based on phylogenetic reconstruction of a diagnostic gene (MutS), and (2) similarities to this viral clade have previously been recovered from both corals and Symbiodinium cultures (), the most parsimonious interpretation of these physical and genetic data is that the Symbiodinium-associated megavirus in this study is a cosmopolitan relative of the Pyramimonas orientalis megaviruses.
In contrast, the identity of the 150 nm NCLDV particles described here is less straightforward. Importantly, these putative viruses are not reminiscent of the originally described mimiviruses (; C. Desnues, personal communication), but more similar to phycodnaviruses and iridoviruses (), as well as megaviruses () in that they lack the characteristic capsid hair-like projections of mimiviruses (Figure 5). Many phylogenetically relevant genes within the Phycodnaviridae were annotated in this work (Table 1). For example, 50 unique contig similarities were found to a single member, Paramecium bursaria Chlorella virus NY2A, of this phycovirus Family. Although phycodnaviruses are not generally thought to associate with multicellular eukaryotes, recent evidence has shown that they can infect non-algal hosts, including humans (Yolken et al., 2014). Thus, we suggest that these Phycodnaviridae similarities represent a phycovirus relative that infects the coral host. An alternative, but unlikely explanation, is that these 150 nm viruses do infect Symbiodinium as they are similar to those in but were only visualized moving through the host en route to the free-living environment.
Atypical Herpes-Like Viruses in Corals
A major interest and perplexing aspect of this work, was the high prevalence of the ∼120 nm VLPs that were morphologically similar to and biologically different from herpesviruses. In this study, VLPs of this average size were visually indistinguishable from many herpesviruses in terms of their capsid size and envelope. However, described herpesviruses predominantly replicate in the nuclei of cells (), whereas herpes-like particles in this study were never identified in the nuclei of coral cells. Thus, we hypothesize that the ∼120 nm VLPs described here are not true herpesviruses, but perhaps something distinct.
Further, although we (; ; ) and others (; ) have found metagenomic evidence that herpes-like annotations can dominate dsDNA viral taxa in corals, these similarities align with only a few coding regions of herpesvirus genomes, have low percent identity to known herpesviruses, and rarely span large portions of the herpesvirus genomes (see ). In the unrooted DNA polymerase gene tree, the A. aspera contig is placed in a clade containing all primate taxa in the phylogeny except the squirrel monkey, Saimiri sciureus (Supplementary Figure S2). The exact position of the A. aspera DNA polymerase contig within this clade could not be resolved, but the sequences generated in this study appear to be relatively distinct from previously sequenced herpesvirus DNA polymerase genes. Thus we hypothesize that these viruses are ‘atypical’ in that they are highly morphologically reminiscent of herpesviruses, but only marginally similar to herpesviruses at the genomic and cell cycle levels. What these novel viruses truly are and how they affect their host corals remains an intriguing question. Future investigations should aim to evaluate the genomes of these viruses, perhaps by using deeper and longer sequencing approaches (e.g., PacBio), cell culture-based work, or size selection-based flow sorting ().
Multiple Lines of Evidence for Retroviruses in Acroporid Corals
Retrovirus-like sequence similarities have previously been characterized from stony corals. For example, they comprised 6.8 and 10%, respectively, of sequence similarities obtained from control and heat-stressed viromes generated from the Caribbean coral, Montastraea cavernosa (). The number of retrovirus-like sequence similarities recovered in this study from Acropora aspera (12.7%) is thus comparable to previous findings. Although VLPs physically similar to retroviruses were observed in the Acropora aspera samples used for virome generation, they were encountered relatively infrequently (Figures 7C,D). These particles were ∼70 nm in diameter with a somewhat amorphous or egg-shaped morphology reminiscent of amphotropic retroviruses (; Figure 1 on page 477 of ; ). Yet, the most abundant set of unique eukaryotic virus similarities were to the Retroviridae (n = 112 best contig annotations, Figure 3). These annotations included all the structural genes specifically important to this group including gag, pro, pol, and env, as well as non-structural genes, such as an integrase (Table 1). A majority (73%) of these annotations were from gammaretroviruses most similar to the Reticuloendotheliosis virus and Porcine type-C oncovirus genomes. The next most abundant (12.6%) of these annotations were to a spumavirus, Simian foamy virus. Several of these annotations were to a special group of reverse transcriptases in these viruses that contain RT_ZFREV_like domains (Table 1) that are only found in true retroviruses and not retro-elements, confirming their viral origin ().
With regards to RNA viral diversity, it should be noted that we only enriched for DNA viruses; any RNA virus that does not have an intermediate DNA stage would have been missed. If the small (∼17–40 nm) and abundant filamentous viruses cataloged in a few of the TEMs (Figures 7E,F) were RNA-based as hypothesized previously (; ; ; ), then the virome data would tell us little about their genomic identity. In a previous effort to enrich for RNA genomes (), we identified five sequence similarities to Heterocapsa circularisquama RNA virus (HcRNAV), a +ssRNA virus that infects free-living dinoflagellates. However, sequence similarities to HcRNAV were not observed in this study and, overall, there remains little molecular data to suggest that RNA viruses similar to previously described groups are a major component of the acroporid coral virome. Additional work should be conducted in this capacity.
Viral Outbreak in Corals Associated with a Reef Flat Bleaching Event
This study characterizes the viral consortia from acroporid colonies that: (1) experienced aerial exposure and hyposmotic stress in situ on the reef flat, (2) were collected and acclimated to flow-through experimental tanks, and then (3) experienced experimental injection and heat (A. aspera only) treatments. Microscopic and genomic data were then interpreted from these samples. No experimental corals were bleached at the time of collection or at the start of the experimental treatments, yet the collected coral colonies did experience abiotic stressors on the reef flat prior to collection for this experiment (Figure 2, Supplementary Material). These stressors could have primed or triggered viral lytic cycles in the experimental corals that were only evident at the end of the experimental period (in TEM and genomic investigations). VLP morphologies representative of all predominant groups identified from the metagenome were observed in all experimental treatments (including non-injected controls), and no potential effects from injection treatments were visually evident in any coral fragments. Thus, we infer that the documented viral outbreak most parsimoniously reflects a common response among all experimental corals to aerial exposure and hyposaline (rainfall) stressors on the reef prior to collection for this study, and is not likely driven by our experimental treatments.
We interpret this event as an outbreak for several reasons. First, we identified a diversity of viruses from a relatively small number of samples and sequencing effort, yet observed VLP abundances that were two to three times higher than those previously reported from Acropora muricata at Heron Island and Lizard Island, Great Barrier Reef (i.e., 2 to 20 VLPs per cell or membrane-bound vacuole for all VLP size ranges; ). Further, similar viruses were detected from all experimental fragments of two different acroporid species (based on TEM). Relative to experimental corals, conspecific acroporids that remained on the reef flat during the study period experienced additional episodes of aerial exposure and coincident intense rainfall (on March 21st–22nd), which likely triggered the observed massive bleaching (e.g., ). Thus, we hypothesize that our experimental corals and many Heron Island reef flat acroporids had high viral loads simultaneously.
Conclusion
The combined application of physical and genomic-based methods in this study provided some significant benefits over using either approach in isolation. Since the TEMs and virome generated in this experiment contained evidence of numerous and diverse putative viruses, we suspect that environmental conditions (low tide-driven aerial exposure coincident with hyposaline conditions due to heavy rainfall) on the reef flat led to mass bleaching on the Heron Island reef flat in March of 2011, which was associated with a viral outbreak. This suggests that stressful environmental conditions can rapidly trigger the onset of viral infection by multiple etiological agents (e.g., atypical herpes-like virus, NCLDV-like viruses including megaviruses) concurrently, and highlights our uncertainty regarding the disease signs exhibited in coral viral infections. Future studies should explore whether increased viral loads are ubiquitous in bleached corals regardless of the stressor triggering bleaching.
Statements
Author contributions
AC, AT, and RV designed and implemented the study, and collected data in the field. TA performed the microscopy, and AC, TA and RV interpreted the microscopy results. CB created the viral metagenome; SR and RV analyzed it; and AC, SR, and RV interpreted the metagenome results. AC, TA, RS, and RV wrote the manuscript. All authors edited the manuscript.
Funding
. This research was supported by a National Science Foundation grant (OCE-0960937) to RT.
Acknowledgments
We thank W. Leggat for comprehensive logistical support and M. van Oppen for providing us with SYBR Gold. We are grateful to members of the D. Miller lab for their intellectual input and for assistance with field logistics, S. Dove for confirming our coral colony identifications, and D. Kline for his assistance with A. aspera aquaria setup. We thank the Heron Island Research Station Staff bsts for their logistical support of our experiments. We also thank three anonymous reviewers for their insightful suggestions on an earlier version of this manuscript.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Supplementary material
The Supplementary Material for this article can be found online at: http://journal.frontiersin.org/article/10.3389/fmicb.2016.00127
Footnotes
1.^“NOAA Declares Third Ever Global Coral Bleaching Event.” National Oceanic and Atmospheric Administration. NOAA, October 8, 2015. Web. January 12, 2016.
References
1
BakerA. C.CunningR. (2015). “Coral “bleaching” as a generalized stress response to environmental disturbance,” in Diseases of Coral, edsWoodleyC. M.DownsC. A.BrucknerA. W.PorterJ. W.GallowayS. B. (Hoboken, NJ:John Wiley & Sons, Inc.), 396–409.
2
BettarelY.BouvierT.NguyenH. K.ThuP. T. (2015). The versatile nature of coral-associated viruses.Environ. Microbiol.173433–3439. 10.1111/1462-2920.12579
3
BettarelY.ThuyN. T.HuyT. Q.HoangP. K.BouvierT. (2013). Observation of virus-like particles in thin sections of the bleaching scleractinian coral Acropora cytherea.J. Mar. Biol. Assoc.93909–912. 10.1017/S0025315411002062
4
ClaverieJ. M.GrzelaR.LartigueA.BernadacA.NitscheS.VaceletJ.et al (2009). Mimivirus and Mimiviridae: giant viruses with an increasing number of potential hosts, including corals and sponges.J. Invertr. Pathol.101172–180. 10.1016/j.jip.2009.03.011
5
ColsonP.De LamballerieX.YutinN.AsgariS.BigotY.BideshiD. K.et al (2013). “Megavirales,” a proposed new order for eukaryotic nucleocytoplasmic large DNA viruses.Arch. Virol.1582517–2521. 10.1007/s00705-013-1768-6
6
ColsonP.GimenezG.BoyerM.FournousG.RaoultD. (2011). The giant Cafeteria roenbergensis virus that infects a widespread marine phagocytic protist is a new member of the fourth domain of life.PLoS ONE6:e18935. 10.1371/journal.pone.0018935
7
CorreaA. M. S.WelshR. M.ThurberR. L. V. (2013). Unique nucleocytoplasmic dsDNA and +ssRNA viruses are associated with the dinoflagellate endosymbionts of corals.ISME J.713–27. 10.1038/ismej.2012.75
8
DavyJ. E.PattenN. L. (2007). Morphological diversity of virus-like particles within the surface microlayer of scleractinian corals.Aquat. Microb. Ecol.4737–44. 10.3354/ame047037
9
DeanF. B.NelsonJ. R.GieslerT. L.LaskenR. S. (2001). Rapid amplification of plasmid and phage DNA using phi29 DNA polymerase and multiply-primed rolling circle amplification.Genome Res.111095–1099. 10.1101/gr.180501
10
De’AthG.FabriciusK. E.SweatmanH.PuotinenM. (2012). The 27-year decline of coral cover on the great barrier reef and its causes.Proc. Natl. Acad. Sci. U.S.A.10917995–17999. 10.1073/pnas.1208909109
11
FischerM. G.AllenM. J.WilsonW. H.SuttleC. A. (2010). Giant virus with a remarkable complement of genes infects marine zooplankton.Proc. Natl. Acad. Sci. U.S.A.10719508–19513. 10.1073/pnas.1007615107
12
GardnerT.CoteI.GillJ.GrantA.WatkinsonA. (2003). Long-term region-wide declines in Caribbean corals.Science301958–960. 10.1126/science.1086050
13
IyerL. A.BalajiS.KooninE. V.AravindL. (2006). Evolutionary genomics of nucleo-cytoplasmic large DNA viruses.Virus Res.117156–184. 10.1016/j.virusres.2006.01.009
14
KingA. M. Q.AdamsM. J.CarstensE. B.LefkowitzE. J. (2012). Virus Taxonomy: Classification and Nomenclature of Viruses: Ninth Report of the International Committee on Taxonomy of Viruses.San Diego, CA:Elsevier Academic Press.
15
La ScolaB.AudicS.RobertC.JungangL.de LamballerieX.DrancourtM.et al (2003). A giant virus in amoebae.Science299:2033. 10.1126/science.1081867
16
LawrenceS. A.WilkinsonS. P.DavyJ. E.ArlidgeW. N. S.WilliamsG. J.WilsonW. H.et al (2015). Influence of local environmental variables on the viral consortia associated with the coral Montipora capitata from Kaneohe Bay, Hawaii, USA.Aquat. Microb. Ecol.74251–262. 10.3354/ame01743
17
LawrenceS. A.WilsonW. H.DavyJ. E.DavyS. K. (2014). Latent virus-like infections are present in a diverse range of Symbiodinium spp. (Dinophyta).J. Phycol.50977–997. 10.1111/jpy.12242
18
LegendreM.BartoliJ.ShmakovaL.JeudyS.LabadieK.AdraitA.et al (2014). Thirty-thousand-year-old distant relative of giant icosahedral DNA viruses with a pandoravirus morphology.Proc. Natl. Acad. Sci. U.S.A.1114274–4279. 10.1073/pnas.1320670111
19
LerusteA.BouvierT.BettarelY. (2012). Enumerating viruses in coral mucus.Appl. Environ. Microbiol.786377–6379. 10.1128/AEM.01141-12
20
Le TissierM. D. A. (1991). The nature of the skeleton and skeletogenic tissues in the Cnidaria.Hydrobiologia21397–402. 10.1007/BF00026492
21
LittmanR. A.WillisB. L.BourneD. G. (2011). Metagenomic analysis of the coral holobiont during a natural bleaching event on the Great Barrier Reef.Environ. Microbiol. Rep.3651–660. 10.1111/j.1758-2229.2010.00234.x
22
LohrJ.MunnC. B.WilsonW. H. (2007). Characterization of a latent virus-like infection of symbiotic zooxanthellae.Appl. Environ. Microbiol.732976–2981. 10.1128/AEM.02449-06
23
MaierC.De KluijverA.AgisM.BrussaardC. P. D.Van DuylF. C.WeinbauerM. G. (2011). Dynamics of nutrients, total organic carbon, prokaryotes and viruses in onboard incubations of cold-water corals.Biogeosciences82609–2620. 10.5194/bg-8-2609-2011
24
MarhaverK. L.EdwardsR. A.RohwerF. (2008). Viral communities associated with healthy and bleaching corals.Environ. Microbiol.102277–2286. 10.1111/j.1462-2920.2008.01652.x
25
MartinezM. J.SwanB. K.WilsonW. H. (2014). Marine viruses, a genetic reservoir revealed by targeted viromics.ISME J.81079–1088. 10.1038/ismej.2013.214
26
MonierA.LarsenJ. B.SandaaR. A.BratbakG.ClaverieJ. M.OgataH. (2008). Marine Mimivirus relatives are probably large algal viruses.Virol. J.5:12. 10.1186/1743-422X-5-12
27
MoniruzzamanM.LecleirG. R.BrownC. M.GoblerC. J.BidleK. D.WilsonW. H.et al (2014). Genome of brown tide virus (AaV), the little giant of the Megaviridae, elucidates NCLDV genome expansion and host-virus coevolution.Virology46660–70. 10.1016/j.virol.2014.06.031
28
Nguyen-KimH.BouvierT.BouvierC.HaiD. N.LamN. N.Rochelle-NewallE.et al (2014). High occurrence of viruses in the mucus layer of scleractinian corals.Environ. Microbiol. Rep.6675–682. 10.1111/1758-2229.12185
29
NobleR. T.FuhrmanJ. A. (1998). Use of SYBR Green I for rapid epifluorescence counts of marine viruses and bacteria.Aquat. Microb. Ecol.14113–118. 10.3354/ame014113
30
PattenN. L.HarrisonP. L.MitchellJ. G. (2008). Prevalence of virus-like particles within a staghorn scleractinian coral (Acropora muricata) from the Great Barrier Reef.Coral Reefs27569–580. 10.1007/s00338-008-0356-9
31
PollockF. J.Wood-CharlsonE. M.Van OppenM. J. H.BourneD. G.WillisB. L.WeynbergK. D. (2014). Abundance and morphology of virus-like particles associated with the coral Acropora hyacinthus differ between healthy and white syndrome-infected states.Mar. Ecol. Prog. Ser.51039–43. 10.3354/meps10927
32
RetenoD. G.BenamarS.KhalilJ. B.AndreaniJ.ArmstrongN.KloseT.et al (2015). Faustovirus, an Asfarvirus-related new lineage of giant viruses infecting amoebae.J. Virol.896585–6594. 10.1128/JVI.00115-15
33
RosarioK.SchenckR. O.HarbeitnerR. C.LawlerS. N.BreitbartM. (2015). Novel circular single-stranded DNA viruses identified in marine invertebrates reveal high sequence diversity and consistent predicted intrinsic disorder patterns within putative structural proteins.Front. Microbiol.6:696. 10.3389/fmicb.2015.00696
34
RosenbergE.Zilber-RosenbergI. (2014). “Microbiotas are part of holobiont fitness,” in The Hologenome Concept: Human, Animal, and Plant Microbiota, edsRosenbergE.Zilber-RosenbergI. (New York, NY:Springer International Publishing), 55–80.
35
SandaaR. A.HeldalM.CastbergT.ThyrhaugR.BratbakG. (2001). Isolation and characterization of two viruses with large genome size infecting Chrysochromulina ericina (Prymnesiophyceae) and Pyramimonas orientalis (Prasinophyceae).Virology290272–280. 10.1006/viro.2001.1161
36
SchlabergR.ChoeD. J.BrownK. R.ThakerH. M.SinghI. R. (2009). XMRV is present in malignant prostatic epithelium and is associated with prostate cancer, especially high-grade tumors.Proc. Natl. Acad. Sci. U.S.A.10616351–16356. 10.1073/pnas.0906922106
37
SchmidM.SpeisederT.DobnerT.GonzalezR. A. (2014). DNA virus replication compartments.J. Virol.881404–1420. 10.1128/JVI.02046-13
38
SharmaV.ColsonP.GiorgiR.PontarottiP.RaoultD. (2014). DNA-dependent RNA polymerase detects hidden giant viruses in published databanks.Genome Biol. Evol.61603–1610. 10.1093/gbe/evu128
39
ShenC. H.SteinerL. A. (2004). Genome structure and thymic expression of an endogenous retrovirus in zebrafish.J. Virol.78899–911. 10.1128/JVI.78.2.899-911.2004
40
SofferN.BrandtM. E.CorreaA. M. S.SmithT. B.Vega ThurberR. L. (2014a). Potential role of viruses in white plague coral disease.ISME J.8271–283. 10.1038/ismej.2013.137
41
SofferN.ZaneveldJ.Vega ThurberR. L. (2014b). Phage-bacteria network analysis and its implication for the understanding of coral disease.Environ. Microbiol.171203–1218. 10.1111/1462-2920.12553
42
TakemuraM.YokoboriS.OgataH. (2015). Evolution of eukaryotic DNA polymerases via interaction between cells and large DNA viruses.J. Mol. Evol.8124–33. 10.1007/s00239-015-9690-z
43
ThurberR. L. V.HaynesM.BreitbartM.WegleyL.RohwerF. (2009). Laboratory procedures to generate viral metagenomes.Nat. Protoc.4470–483. 10.1038/nprot.2009.10
44
van OppenM. J. H.LeongJ. A.GatesR. D. (2009). Coral-virus interactions: a double-edged sword?Symbiosis471–8. 10.1007/BF03179964
45
Vega ThurberR. L.BarottK. L.HallD.LiuH.Rodriguez-MuellerB.DesnuesC.et al (2008). Metagenomic analysis indicates that stressors induce production of herpes-like viruses in the coral Porites compressa.Proc. Natl. Acad. Sci. U.S.A.10518413–18418. 10.1073/pnas.0808985105
46
Vega ThurberR. L.CorreaA. M. S. (2011). Viruses of reef-building scleractinian corals.J. Exp. Mar. Biol. Ecol.408102–113. 10.1016/j.jembe.2011.07.030
47
WegleyL.EdwardsR.BeltranR.-B.HongL.RohwerF. (2007). Metagenomic analysis of the microbial community associated with the coral Porites astreoides.Environ. Microbiol.92707–2719. 10.1111/j.1462-2920.2007.01383.x
48
WestonA. J.DunlapW. C.ShickJ. M.KlueterA.IglicK.VukelicA.et al (2012). A profile of an endosymbiont-enriched fraction of the coral Stylophora pistillata reveals proteins relevant to microbe-host interactions.Mol. Cell. Proteom.11:M111.015487. 10.1074/mcp.M111.015487
49
WeynbergK. D.VoolstraC. R.NeaveM. J.BuergerP.Van OppenM. J. H. (2015). From cholera to corals: viruses as drivers of virulence in a major coral bacterial pathogen.Sci. Rep.5:17889. 10.1038/srep17889
50
WeynbergK. D.Wood-CharlsonE. M.SuttleC. A.Van OppenM. J. H. (2014). Generating viral metagenomes from the coral holobiont.Front. Microbiol.5:206. 10.3389/fmicb.2014.00206
51
WilsonW. H. (2011). “Coral viruses,” in Studies in Viral Ecology, Volume Two: Animal Host Systems, ed. HurstC. J. (Hoboken, NJ:John Wiley & Sons, Inc.), 141–149.
52
WilsonW. H.DaleA. L.DavyJ. E.DavyS. K. (2005). An enemy within? Observations of virus-like particles in reef corals.Coral Reefs24145–148. 10.1007/s00338-004-0448-0
53
WilsonW. H.FrancisI.RyanK.DavyS. K. (2001). Temperature induction of viruses in symbiotic dinoflagellates.Aquat. Microb. Ecol.2599–102. 10.3354/ame025099
54
WilsonW. H.GilgI. C.DuarteA.OgataH. (2014). Development of DNA mismatch repair gene, MutS, as a diagnostic marker for detection and phylogenetic analysis of algal Megaviruses.Virology466123–128. 10.1016/j.virol.2014.07.001
55
Wood-CharlsonE. M.WeynbergK. D.SuttleC. A.RouxS.Van OppenM. J. H. (2015). Metagenomic characterization of viral communities in corals: mining biological signal from methodological noise.Environ. Microbiol.173440–3449. 10.1111/1462-2920.12803
56
YamadaT. (2011). Giant viruses in the environment: their origins and evolution.Curr. Opin. Virol.158–62. 10.1016/j.coviro.2011.05.008
57
YolkenR. H.Jones-BrandoL.DuniganD. D.KannanG.DickersonF.SeveranceE.et al (2014). Chlorovirus ATCV-1 is part of the human oropharyngeal virome and is associated with changes in cognitive functions in humans and mice.Proc. Natl. Acad. Sci. U.S.A.11116106–16111. 10.1073/pnas.1418895111
58
YutinN.KooninE. V. (2012). Hidden evolutionary complexity of Nucleo-Cytoplasmic large DNA viruses of eukaryotes.Virol. J.91–18. 10.1186/1743-422X-9-161
59
ZerbinoD. R.BirneyE. (2008). Velvet: algorithms for de novo short read assembly using de Bruijn graphs.Genome Res.18821–829. 10.1101/gr.074492.107
60
ZhangJ.KobertK.FlouriT.StamatakisA. (2014). PEAR: a fast and accurate Illumina paired-end reAd mergeR.Bioinformatics30614–620. 10.1093/bioinformatics/btt593
Summary
Keywords
virome, tropical coral reef, virus-like particle (VLP), herpesvirus, megavirus, nucleocytoplasmic large DNA virus (NCLDV)
Citation
Correa AMS, Ainsworth TD, Rosales SM, Thurber AR, Butler CR and Vega Thurber RL (2016) Viral Outbreak in Corals Associated with an In Situ Bleaching Event: Atypical Herpes-Like Viruses and a New Megavirus Infecting Symbiodinium. Front. Microbiol. 7:127. doi: 10.3389/fmicb.2016.00127
Received
11 September 2015
Accepted
25 January 2016
Published
22 February 2016
Volume
7 - 2016
Edited by
Ian Hewson, Cornell University, USA
Reviewed by
Fabiano Thompson, Federal University of Rio de Janeiro, Brazil; Stéphan Jacquet, Institut National de la Recherche Agronomique, France; Mya Breitbart, University of South Florida, USA
Updates
Copyright
© 2016 Correa, Ainsworth, Rosales, Thurber, Butler and Vega Thurber.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Rebecca L. Vega Thurber, athurber@coas.oregonstate.edu
This article was submitted to Aquatic Microbiology, a section of the journal Frontiers in Microbiology
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.