Genotyping of a microsatellite locus to differentiate clinical Ostreid herpesvirus 1 specimens.

Ostreid herpesvirus 1 (OsHV-1) is a DNA virus belonging to the Malacoherpesviridae family from the Herpesvirales order. OsHV-1 has been associated with mortality outbreaks in different bivalve species including the Pacific cupped oyster, Crassostrea gigas. Since 2008, massive mortality events have been reported among C. gigas in Europe in relation to the detection of a variant of OsHV-1, called μVar. Since 2009, this variant has been mainly detected in France. These results raise questions about the emergence and the virulence of this variant. The search for association between specific virus genetic markers and clinical symptoms is of great interest and the characterization of the genetic variability of OsHV-1 specimens is an area of growing interest. Determination of nucleotide sequences of PCR-amplified virus DNA fragments has already been used to characterize OsHV-1 specimens and virus variants have thus been described. However, the virus DNA sequencing approach is time-consuming in the high-scale format. Identification and genotyping of highly polymorphic microsatellite loci appear as a suitable approach. The main objective of the present study was the development of a genotyping method in order to characterise clinical OsHV-1 specimens by targeting a particular microsatellite locus located in the ORF4 area. Genotyping results were compared to sequences already available. An excellent correlation was found between the detected genotypes and the corresponding sequences showing that the genotyping approach allowed an accuraté discrimination between virus specimens.

Introduction

Ostreid herpesvirus 1 (OsHV-1) is a DNA virus belonging to the Malacoherpesviridae family from the Herpesvirales order [1]. The virus has been purified from naturally infected Crassostrea gigas larvae [2] and its genome entirely sequenced [3]. The viral genome is a large linear duplex DNA molecule of 207 kb (GenBank accession number AY509253) that encodes at least 124 genes [3].

OsHV-1 has been associated with mortality outbreaks in different bivalve species including the Pacific cupped oyster, C. gigas[4,5]. Since 2008, massive mortality outbreaks have been reported among C. gigas spat in several farming areas in Europe [6-11] in relation to the detection of a newly described OsHV-1 variant called μVar [7]. Although the reference type (a viral specimen collected in France in 1995 during a mortality event affecting C. gigas larvae, GenBank accession number AY509253) and the variant μVar were detected in association with mortality outbreaks in 2008 in France, virus detection since 2009 has mainly concerned the μVar variant [7,10,12]. These results raise questions about the emergence and the virulence of the μVar variant. In this context, tools are needed in order to better describe OsHV-1 diversity in relation to virulence and geographical distribution. In light of the genetic diversity of OsHV-1, the search for associations between specific virus genetic markers and clinical symptoms is of great interest.

Determination of nucleotide sequences of PCR-amplified virus DNA fragments is the most accurate method for virus genotyping [13]. The DNA sequencing approach has been used to characterise OsHV-1 specimens and virus variants were thus reported [7,10,12,14-19]. The μVar variant [7] showed several differences in two genome areas when compared with the reference type (GenBank accession n° AY509523) and all these differences need to be observed to define a viral specimen as the μVar variant.

Virus DNA sequencing is, however, time-consuming in the high-scale format. The identification and genotyping of highly polymorphic microsatellite areas from vertebrate herpesviruses appears as a suitable approach. Microsatellites have been reported from different herpesviruses including human cytomegalovirus and they have been used as molecular markers to define virus polymorphism [20-24].

Since the μVar variant demonstrated a deletion of 12 bp in a microsatellite locus located up-stream of the ORF4 [7], the main objective of the present study was the development of a genotyping method. This method was used to characterise 47 clinical OsHV-1 specimens by targeting this microsatellite locus. DNA sequences already available were used to compare results obtained with both techniques. Sequencing and genotyping appeared to be equally useful to differentiate clinical OsHV-1 specimens.

Materials and methods

Oyster samples and C2/C6 sequences

Forty-seven samples of the Pacific cupped oyster, C. gigas, were selected in the present study in order to analyze them by genotyping. These included animals collected from 1993 to 2010 and covered different stages of development (larvae, spat and adults) (Table 1). Most of the samples (45) were collected in France during mortality outbreaks recorded by the national network for mollusc disease monitoring (Repamo, Ifremer) and were stored frozen at −20 °C. Two samples were of different geographical origins (Japan and USA) (Table 1).

Table 1

Collection of DNA extracted from samples: geographical origins, virus DNA amounts and related GenBank accession numbers

Isolate code	Origin	Nature of sample	Viral DNA quantity (genome copies ng - 1 of total DNA)	GenBank accession number
				C2/C6
1993/002	France	Larvae	5.15E + 04	JN80065
1994/005	France	< 1 year	6.54E + 03	JN80068
1994/006	France	< 1 year	4.04E + 04	JN80069
1994/011	France	< 1 year	6.60E + 04	JN80070
1994/012	France	Larvae	3.79E + 05	JN80071
1995/020	France	Larvae	1.35E + 05	JN80072
2005/001	France	< 1 year	7.05E + 04	JN80081
2005/005	France	< 1 year	1.23E + 05	JN80082
2005/012	France	< 1 year	9.96E + 04	JN80084
2006/005	France	Larvae	3.21E + 03	JN80087
2006/009	France	< 1 year	2.21E + 04	JN80088
2006/018	France	< 1 year	1.31E + 04	JN80090
2007/004	France	Larvae	1.03E + 06	JN80091
2007/026	France	< 1 year	1.41E + 04	JN80094
2007/028	France	< 1 year	1.00E + 04	JN80095
2007/029	France	< 1 year	1.04E + 05	JN80096
2007/030	France	< 1 year	1.21E + 04	JN80097
2007/034	France	1-2 years	1.18E + 04	JN80098
2007/035	France	< 1 year	8.93E + 03	JN80099
2008/017	France	< 1 year	4.30E + 02	JN800100
2008/019	France	1-2 years	8.43E + 06	JN800101
2008/021	France	Larvae	2.45E + 06	JN800102
2008/023	France	Adult	2.65E + 05	JN800103
2008/020	France	Larvae	1.66E + 03	-
2008/025	France	Adult	4.33E + 05	JN800104
2008/045	France	1-2 years	8.60E + 02	JN800107
2008/050	France	< 1 year	4.68E + 02	JN800108
2008/055	France	< 1 year	2.56E + 05	JN800109
2008/059	France	< 1 year	2.60E + 05	JN800110
2008/073	France	1-2 years	4.55E + 04	JN800111
2008/079	France	< 1 year	1.45E + 04	JN800112
2008/083	France	< 1 year	2.74E + 04	JN800113
2008/092	France	< 1 year	2.24E + 03	JN800114
2009/002	France	< 1 year	2.67E + 03	JN800115
2009/021	France	< 1 year	7.17E + 03	JN800116
2009/022	France	< 1 year	1.04E + 04	JN800117
2009/027	France	< 1 year	3.31E + 03	JN800118
2009/035	France	< 1 year	2.19E + 03	JN800119
2010/002	France	< 1 year	4.44E + 04	JN800120
2010/008	France	< 1 year	2.57E + 04	JN800121
2010/012	France	< 1 year	1.58E + 05	JN800122
2010/013	France	< 1 year	1.61E + 05	JN800123
2010/021	France	< 1 year	8.38E + 03	JN800124
2010/023	France	< 1 year	6.42E + 05	JN800125
2010/026	France	< 1 year	5.98E + 05	JN800126
2007/07-CB2	USA	-	5.02E + 04	JN800128
2010/158-144	Japan	< 1 year	3.00E + 02	JN800133

Total nucleic acids were previously extracted from oyster samples using the QIAamp DNA Mini Kit (Qiagen, Courtaboeuf, France) [25] and the quantity of viral DNA was estimated by real time PCR using the primer pair C9/C10 [26] for the purpose of a previous study [10]. Sequences of C2/C6 [27] PCR products from the sample set were previously defined in the laboratory [10] and deposited in GenBank (Table 1).

Genotyping of a microsatellite locus

A microsatellite locus (called H10) was identified up-stream of ORF4 based on an analysis of the whole genome of OsHV-1 (GenBank accession n° AY509523) using the MsatFinder algorithm [28] and selected for genotyping. The number of repeat units of this trinucleotide microsatellite was 8 in the reference sequence (GenBank accession n° AY509523). This microsatellite (H10) was located in inverted repeated regions (4440–4463 positions and 178547–178547 positions). A primer pair (H10F/H10R) was designed with primer3 [29] (H10F: gtgatggctttggtcaaggt and H10R: ggcgcgatttgtcagtttag). The expected size of the PCR product was 151 bp.

PCR was performed in 20 μL reaction volumes in a thermocycler (Applied Biosystems, Villebon-sur-Yvette, France). Each reaction contained 11.76 μL dH2O, 4.0 μL buffer 5X (Promega, Charbonnières-les-Bains, France), 0.8 μL MgCl2 (25 mM), 2 μL dNTP (2 mM), 0.18 μL of each primer (20 μM), and 0.08 μL GoTaq®DNA polymerase (5U/μL, Promega). PCR cycling conditions were 96 °C for 5 min; then 30 cycles of 95 °C for 30 s, 57 °C for 30 s and 72 °C for 30 s; and finally a step of 72 °C for 2 min. PCR products were verified by electrophoresis (1% agarose gel).

PCR products were mixed with formamide and GeneScan 500-ROX size standard (Applied Biosystems) respectively according to the manufacturer’s recommendations (1.5 μL PCR products, 0.25 μL Rox size standard and 13.25 μL formamide). After 5 min denaturation followed by rapid cooling, PCR products were detected using an ABI 3130xl Genetic Analyzer (Applied Biosystems), and the fragment length was estimated through the GeneMapper 3.7 software.

Phylogenetic analysis

Phylogenetic analysis was performed on C2/C6 [27] sequences [10] using 3 computational approaches. Information concerning genotyping results (length of the fragment) was included in specimen codes. For the first approach, a phylogenetic tree was created from sequence alignments using the Neighbor-Joining (NJ) method [30]. The significance of the branching orders was assessed by bootstrap resampling of 1000 replicates. The second approach was based on phylogeny inference according to the Maximum Likelihood method based on the Tamura-Nei model [31]. Bootstrap data sets (1000 replicates) were generated. The Maximum Parsimony method was also used as the third approach. All approaches were implemented using the MEGA5 program [32].

Results

Polymorphism of C2/C6 PCR products

Comparing sequences of C2/C6 PCR products (Table 1) demonstrated a high polymorphism with 82 positions of a 482 nucleotide sequence (17%) showing a substitution/deletion/insertion defining 9 virus specimen types from the analysed samples (data not shown). H10F/H10R sequences revealed only a variability in the number of repeat units at the targeted short tandem repeat (H10) defining 7 virus specimen types from the analysed samples (Figure 1). An additional type corresponded to acute viral necrosis virus (AVNV) infecting cultured scallops, Chlamys farreri, in China [33] (Figure 1). The minimum and maximum numbers of repeat units of the trinucleotide motif were 3 (AVNV) and 13 (2006/18/France), respectively (Figure 1).

Figure 1

H10F/H10R sequence alignments between virus specimens. Partial C2-C6 (H10F/H10R) sequence alignments between virus specimens demonstrating variability at the microsatellite locus (H10) and mutation points. Locations of H10F/H10R primers are identified as surligned. OsHV-1 reference type, the variant μVar and AVNV sequences are highlited in grey. Stars represent identity at a particular nucleotide position.

Microsatellite genotyping

The 47 samples were genotyped for the H10 microsatellite locus which is located up-stream of the ORF4. In this region, a 13 bp deletion is one of the characteristics of the μVar variant [7] and was hence chosen to achieve its interest as a diagnostic tool. Protocol optimisation focussed on the concentration of the labelled primers and the DNA in the PCR mix, but also the annealing temperature, time of elongation, and finally the dilution factor of the PCR products in the formamide before fragment length analysis in the Genetic Analyzer.

Fragments were successfully amplified for the 47 samples (Figure 2). Seven different genotypes were detected corresponding to different fragment lengths estimated through the GeneMapper 3.7 software (135, 152, 155, 159, 161, 165 and 168 bp; Figure 2 and Table 2). Two genotypes (135 and 152) were more frequent with 20 and 16 of the samples, respectively (Table 2). Five more genotypes were detected in 11 samples (Table 2). When comparing the detected genotypes and the corresponding H10F/H10R sequences (Table 2), a good correlation was reported (Figure 3) showing that the genotypes identified by genotyping of the microsatellite reflected the sequences and allowed a clear discrimination between them. Moreover, sequencing showed that specimens presenting a fragment length estimated at 135 bp and 152 bp corresponded to the reference type and the μVar variant, respectively. Although all the samples collected in France in 1993, 1994 and 1995 demonstrated a fragment length estimated at 152 bp, French samples collected in 2009 and 2010 showed a single pattern at 135 bp. (Table 2). Different genotypes were detected each year for samples collected in France from 2005 to 2008 (Table 2).

Figure 2

Electrophoregrams from virus specimens. The two traces represent separate reactions on 2 DNA extracts (2008–017 and 2008–020). The boxed numbers under the peaks of the traces are the fragment sizes in base pairs assigned by comparison with the standard curve generated with the internal size standard.

Table 2

Results obtained for each analysed specimen in terms of genotyping and DNA sequencing were grouped per year of collection

Specimen code	Genotyping	Number of repeats (CTA)	Detection frequency (%)	DNA sequencing
1993/002	152,19	8	34	151
1994/005	152,19	8	34	151
1994/006	152,18	8	34	151
1994/011	152,18	8	34	151
1994/012	152,19	8	34	151
1995/020	152,16	8	34	151
2005/001	152,23	8	34	151
2005/005	161,71	11	4.3	160
2005/012	152,27	8	34	151
2006/005	152,27	8	34	151
2006/009	164,88	12	2	163
2006/018	167,91	13	2	166
2007/004	152,35	8	34	151
2007/026	158,69	10	12.7	157
2007/028	161,82	11	4.3	160
2007/029	158,67	10	12.7	157
2007/030	152,18	8	34	151
2007/034	152,32	8	34	151
2007/035	152,32	8	34	151
2008/017	135,28	4	42.5	136
2008/019	158,56	10	12.7	157
2008/020	151,9	8	34	151
2008/021	158,72	10	12.7	157
2008/023	158,68	10	12.7	157
2008/025	158,65	10	12.7	157
2008/045	135,25	4	42.5	136
2008/050	135,18	4	42.5	136
2008/055	135,3	4	42.5	136
2008/059	135,17	4	42.5	136
2008/083	135,18	4	42.5	136
2008/092	135,21	4	42.5	136
2008/073	152,3	8	34	150
2008/079	152,32	8	34	151
2009/002	135,22	4	42.5	136
2009/021	135,14	4	42.5	136
2009/022	135,28	4	42.5	136
2009/027	135,29	4	42.5	136
2009/035	135,29	4	42.5	136
2010/002	135,29	4	42.5	136
2010/008	135,23	4	42.5	136
2010/012	135,24	4	42.5	136
2010/013	135,18	4	42.5	136
2010/021	135,31	4	42.5	136
2010/023	135,17	4	42.5	136
2010/026	135,24	4	42.5	136
2007/07-CB2	155,47	9	2	154
2010/158-144	135,17	4	42.5	136

Although several genotypes were reported the same year in 2005, 2006, 2007 and 2008, a single one was noticed the other years.

Figure 3

Correlation between genotype corresponding sequence lengths. Correlation between the genotypes detected by genotyping and the corresponding sequence lengths of the fragment obtained by sequencing for virus specimens.

The phylogenetic trees built from the C2/C6 amplicon sequences (ORF4 and its related up stream area) using 3 different approaches allowed identification of 2 major groups from the 47 analysed virus specimens (Figure 4). A first group contained French specimens collected from 1993 to 2008 including the reference type (OsHV-1, GenBank n° accession AY509253) and the sample collected in the USA (California) in 2007. Although a main genotype (152 bp) was represented in this group, several other genotypes were also observed (155, 159, 162, 165 and 168 bp) (Figure 4). The second large group comprised French specimens collected from 2008 to 2010. It included the sequence of the μVar variant deposited in GenBank (accession n° HQ842610) and the sequence of the sample from Japan (Figure 4). All samples grouped with the μVar variant demonstrated a similar genotype at 135 bp (Figure 4). Two additional groups were defined. Each of them included a single member: a virus specimen collected in France in 1993 on the one hand and AVNV on the other hand.

Figure 4

Phylogenetic tree representing the relationships of virus specimens. Phylogenetic tree representing the relationships of 47 virus specimens (fragment lengths obtained by genotyping were included in specimen codes) and 3 reference sequences (OsHV-1, the variant μVar and AVNV) based on a fragment of the ORF4 and its up-stream zone (460 nts). Fragment lengths were included in specimen codes. The analysis involved 50 nucleotide sequences. Evolutionary analysis was conducted in MEGA5. The tree was generated by the Maximum Likelihood method.

Discussion

This study reports for the first time the use of a microsatellite locus (H10) present in the OsHV-1 genome to analyze the virus diversity using 47 OsHV-1 specimens. Microsatellites are short tandem repeats that occur in eukaryote, prokaryote, and also some virus genomes. They are highly DNA mutable sequences and represent hot spots of length mutation. Replication slippage errors are considered as the main cause of insertions and deletions at microsatellite loci. Microsatellites have thus been extensively used as molecular markers in numerous genetic diversity and genome mapping studies. Short microsatellite polymorphisms have already been used to describe genetic polymorphism of different vertebrate herpesviruses [20-24]. As an example, Deback et al. [24] used the microsatellite technology to determine genetic relationship between HSV-1 strains and showed that each patient was characterized by its own HSV-1 microsatellite haplotype.

The microsatellite selected in the present study (H10) is found in a noncoding region. This microsatellite was selected since numerous sequences are already available for this region demonstrating a high level of length polymorphism [7,10,12,19].

In the present study analysis of C2-C6 sequences [10] was first carried out to identify polymorphisms among the selected OsHV-1 specimens and to prepare genotyping. The ORF4 area with 82 substitutions/deletions/insertions appeared highly polymorphic presenting variability in the number of repeat units at the targeted short tandem repeat and a variety of point mutations defining 9 virus types. Sequence alignment allowed identification of polymorphisms among virus specimens interpreted as being the reference type (GenBank AY509253). Several French samples collected from 1993 to 2008 demonstrated 100% identity with the reference type and as such could be identified as OsHV-1 [3]. Other samples collected in France from 2003 to 2008 showed some differences in comparison with the reference type. Finally, a French virus specimen collected in 1993 presented high homologies with the variant OsHV-1 Var [14,16,34]. These results showed that different OsHV-1 variants are represented in the sample set selected for the present study. Acute Viral Necrosis Virus (AVNV) [33] was included in comparing C2-C6 PCR product sequences sinced its complete genome is available in GenBank and it presents the shorter sequence for the H10 microsatellite. The number of sequences from countries other than France used in this study was low. Complementary analysis of additional specimens is ongoing in the laboratory and detailed comparison of sequences would present further epidemiological information on OsHV-1.

Among the 47 samples analysed, 7 different genotypes were detected with 2 more frequent ones. They respectively included specimens interpreted as the reference type and the μVar variant. Five more genotypes were also detected. When comparing the genotypes detected and the corresponding C2/C6 sequences, a good correlation was reported showing that the detected genotypes reflected the sequences and allowed a clear discrimination between specimens. However, the number of virus specimen types (9) obtained by sequencing of C2-C6 PCR products remains higher than the number of genotypes defined by genotyping (7). Although analysis of variation in length through genotyping offers a first order of discrimination, sequencing of alleles and viral length variants adds a second level. Sequencing is a necessary step to obtain maximum resolution between viral specimens also revealing SNP.

The polymorphism reported for the selected microsatellite in the present study confirms the interest of such analysis to describe OsHV-1 genome diversity. Moreover, a multiplex genotyping based on analysis of several microsatellites needs to be developed for OsHV-1. The in silico analysis of the OsHV-1 genome using the MsatFinder algorithm demonstrated the presence of 12 short repeat sites including 4 mononucleotide units, 5 dinucleotide repeats, and 3 trinucleotide repeats (data not shown). Most of the identified microsatellites were localized in noncoding parts of the OsHV-1 genome, except for 3 of them located in ORF 66, 77 and 106, respectively. The number of repeat units of dinucleotide or trinucleotide microsatellites was 5 or 8. The longest mononucleotide sequence was 18 bases and 3 microsatellites were located in inverted repeated regions. As most of OsHV-1 repeats are found in noncoding areas, they can be considered as evolutionarily neutral or nearly so and therefore as suitable markers for epidemiology studies. Such a technique targeting several microsatellite loci may provide a rapid and accurate tool that can be used to compare OsHV-1 specimens and to study the epidemiology of viral infections. Finally, polymorphism of microsatellites may also be used to study viral strain virulence. Perdue et al. [35] reported that the increased virulence of a particular strain of the avian influenzae virus is related to the increase in the length of a microsatellite at the hemagglutinin cleavage site.

In conclusion, genotyping based on microsatellite loci appears as a powerful tool to study OsHV-1 polymorphism and can offer a first level of discrimination between specimens in order to select best candidates for complete genome sequencing. Futhermore comparative diversity studies between the host, Crassostrea gigas, and OsHV-1 can be easily performed using oyster mircosatellite markers [36] and to characterize coevolution in this recently introduced oyster species in Europe [37].

Competing interests

This work was supported by Ifremer (Institut Français pour l’Exploitation de la Mer). The authors declare that they have no competing interests.

Authors’ contributions

This study is the result of a collective work. TR and SL conceived this study, and participated in its design. TR, GT and SL carried out the design of the protocol, with the help of the other authors. GT, NF, PM and AS carried out sample analyses. TR and SL drafted the manuscript. VBS participated in the discussion and modification of the manuscript. All authors read, corrected, and approved the final manuscript.