# Lacanobia oleracea nucleopolyhedrovirus (LaolNPV): A new European species of alphabaculovirus with a narrow host range

Oihane Simón, Martin A. Erlandson, Marie Frayssinet, Trevor Williams, David A. Theilmann, Anne Nathalie Volkoff, Primitivo Caballero

<figure><img src="image_1.png" /></figure>

## To cite this version:

Oihane Simón, Martin A. Erlandson, Marie Frayssinet, Trevor Williams, David A. Theilmann, et al. Lacanobia oleracea nucleopolyhedrovirus (LaolNPV): A new European species of alphabaculovirus with a narrow host range. PLoS ONE, 2017, 12 (4), 23 p. 10.1371/journal.pone.0176171. [http://10.1371/journal.pone.0176171](http://10.1371/journal.pone.0176171). [hal-01594567](https://hal.science/hal-01594567)

HAL Id: hal-01594567

https://hal.science/hal-01594567

Submitted on 26 Sep 2017

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
RESEARCH ARTICLE

*Lacanobia oleracea nucleopolyhedrovirus*
(LaoINPV): A new European species of
alphabaculovirus with a narrow host range

Oihane Simón¹*, Martin A. Erlandson², Marie Frayssinet³, Trevor Williams⁴, David A. Theilmann⁵, Anne-Nathalie Volkoff³, Primitivo Caballero¹,⁶

1 Bioinsecticidas Microbianos, Instituto de Agrobiotecnología, CSIC-UPNA, Gobierno de Navarra, Mutilva Baja, Navarra, Spain, 2 Sakatoon Research and Development Centre, Agriculture and Agri-Food Canada, Saskatoon, Saskatchewan, Canada, 3 UMR 1333 INRA—Université Montpellier 2, Diversité Génomique et Interactions Microorganismes-Hotes (DGIMI), Montpellier, France, 4 Instituto de Ecología AC, Xalapa, Veracruz, Mexico, 5 Summerland Research and Development Centre, Agriculture and Agri-Food Canada, Summerland, British Columbia, Canada, 6 Dpto. Producción Agraria, Universidad Pública de Navarra, Campus Arrosadía s/n, Pamplona, Spain

* oihane.simon@unavarra.es

<figure><img src="image_3.png" /></figure>

Abstract

During an insect sampling program in alfalfa crops near Montpellier, France in 2011, *Lacanobia oleracea* larvae were collected that died due to nucleopolyhedrovirus infection (LaoINPV). This virus was subjected to molecular and biological characterization. The virus was a multiple nucleocapsid NPV that showed similar restriction profiles to *Mamestra configurata* NPV-A (MacoNPV-A) but with significant differences. Polypeptide analysis demonstrated similar proteins in occlusion bodies and occlusion derived virions, to those observed in NPVs from *Mamestra* spp. Terminal sequencing revealed that the genome organization shared similarity with that of MacoNPV-A. The most homologous virus was MacoNPV-A 90/2 isolate (95.63% identity and 96.47% similarity), followed by MacoNPV-A 90/4 strain (95.37% and 96.26%), MacoNPV-B (89.21% and 93.53%) and *M. brassicae* MNPV (89.42% and 93.74%). Phylogenetic analysis performed with *lef-8*, *lef-9*, *polh* and a concatenated set of genes showed that LaolNPV and the *Mamestra* spp. NPVs clustered together with HaMNPV, but with a closer genetic distance to MacoNPV-A strains. The Kimura 2-parameter (K-2-P) distances of the complete genes were greater than 0.05 between LaolNPV and the MbMNPV/MacoNPV-B/HaMNPV complex, which indicates that LaolNPV is a distinct species. K-2-P distances were in the range 0.015–0.050 for comparisons of LaolNPV with MacoNPV-A strains, such that additional biological characteristics should be evaluated to determine species status. While MacoNPV-A was pathogenic to seven lepidopteran species tested, LaolNPV was only pathogenic to *Chrysodeixis chalcites*. Given these findings, *Lacanobia oleracea nucleopolyhedrovirus* should be considered as a new species in the *Alphabaculovirus* genus.

<figure><img src="image_4.png" /></figure>

OPEN ACCESS

**Citation:** Simón O, Erlandson MA, Frayssinet M, Williams T, Theilmann DA, Volkoff A-N, et al. (2017) *Lacanobia oleracea nucleopolyhedrovirus* (LaoINPV): A new European species of alphabaculovirus with a narrow host range. PLoS ONE 12(4): e0176171. https://doi.org/10.1371/journal.pone.0176171

**Editor:** Yi Li, Wuhan Bioengineering Institute, CHINA

**Received:** March 7, 2017

**Accepted:** April 6, 2017

**Published:** April 20, 2017

**Copyright:** © 2017 Simón et al. This is an open access article distributed under the terms of the [Creative Commons Attribution License](https://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

**Data Availability Statement:** All relevant data are within the paper.

**Funding:** This study received financial support from the National Institute of Agricultural Research from the Spanish Ministry project RTA2013-00114-C02-02 (Dr. Oihane Simón) and from the Department of Industry of the Government of Navarre project IIQ1406-RI1 (Dr. Primitivo Caballero). This study was supported in part by funding from Agriculture and Agri-Food Canada
Research Project J000197 (Dr. Martin A. Erlandson). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

**Competing interests:** The authors have declared that no competing interests exist.

## Introduction

Baculoviruses have been isolated from more than 700 insect species [1]. The genus *Alphabaculovirus* is the largest genus of baculoviruses, comprising over 90% of the presently known baculoviruses [2–4]. As well, considerable diversity has been observed among different geographical isolates of a single virus species and within baculovirus isolates [5–7], as exemplified by the isolation of multiple genotypes present within a single isolate [6, 8–10].

The alphabaculovirus populations that infect the bertha armyworm, *Mamestra configurata* (Lepidoptera: Noctuidae), represented a good example of this diversity. Populations of this pest are cyclic, with major regional outbreaks occurring every 6–8 years and lasting up to 3 years in western Canada [11, 12]. Epizootics of nucleopolyhedrovirus (NPV) are often associated with collapses of *M. configurata* larval populations and these viruses are major mortality factors that may dampen *M. configurata* outbreak cycles. A number of NPVs have been isolated from *M. configurata* larval populations [13–16]. Notably, two distinct species of alphabaculovirus that are closely related but distinguished by restriction endonuclease profiles, gene content and biological activity have been identified in *M. configurata*, namely *Mamestra configurata nucleopolyhedrovirus A* (MacoNPV-A) and *Mamestra configurata nucleopolyhedrovirus B* (MacoNPV-B) [15]. In addition, distinct geographical strains of MacoNPV-A, including 90/2 and 90/4, have been isolated from *M. configurata* populations. MacoNPV-A is much more prevalent than MacoNPV-B in western Canada, which might confer the higher infectivity due to host-pathogen evolution [17]. Significantly, MacoNPV-B may be a variant of the European *Mamestra brassicae* NPV (MbMNPV) [3, 4, 6, 11] and appears to have a wider host range than MacoNPV-A [6, 11, 14, 15]. Because of the differences in host range these two viruses could evolve divergently while infecting the same host, resulting in divergence in genomic and phenotypic characteristics [15].

According to the current definition by the International Committee on Taxonomy of Viruses (ICTV), a virus species is a polythetic class of viruses that constitute a replicating lineage and occupy a particular ecological niche [18]. A polythetic class is one whose members have several properties in common, although they do not necessarily all share a single common defining property. In the case of baculoviruses, Jehle et al. [3, 4] proposed the use of genome sequence-based phylogenies in addition to morphological and biological characteristics for classification.

Specifically, they suggested the use of a two-parameter model (K-2-P) proposed by Kimura et al. [19], to estimate the phylogenetic distance between two viruses. Using this method, the viruses in question should be considered as the same species when the K-2-P estimated distance between single genes or concatenated *polh*, *lef-8* and *lef-9* nucleotide sequences is less than 0.015. However, when this distance is larger than 0.05 the baculovirus strains should be considered as constituting different species. For viruses with K-2-P distances between 0.015 and 0.05, complementary information would be necessary for species demarcation. In this sense, morphological, pathological and ecological characteristics such as host species or geographical origin should be taken into account when classifying these viruses [4, 20].

During a sampling program in alfalfa crops near Montpellier in 2011 several lepidopteran larvae, presumably *Lacanobia oleracea* (synonymous of *Mamestra oleracea*), were collected with the typical signs of NPV infection. The present study aimed to characterize these isolates using molecular and phylogenetic tools in combination with host range studies to establish a key aspect of the ecology of the virus. Therefore, the discovery of NPV *L. oleracea* infected larvae in France offered an opportunity to apply the established baculovirus species demarcation criteria to determine whether these isolates represent a novel baculovirus species.
## Materials and methods

### Virus isolates, insect source and rearing

In April 2011, lepidopteran larvae were collected from alfalfa crops and the weeds present (mainly grasses and wild brassicas) near Montpellier, France. Two sites were examined for the presence of Lepidoptera; one located beside the village of Lattes (43°33'13.9"N, 3°56'10.0"E) and the other beside the village of Marsillargues (43°39'15.0"N 4°11'25.9"E) in southern France. These fields had been used for many years for the production of pasture and alfalfa. No specific permissions were required for access to the land and the field studies did not involve endangered or protected species.

The noctuids *Agrochola lychnidis*, *Lacanobia oleracea* and *Xylena exsoleta* were identified as the most abundant species present, whereas *Aporophyla australis* occurred in low numbers. In this region *L. oleracea* is common, especially during warm summers (M. Frayssenet, pers. obs.), although the larvae and adults of this species can be confused with *A. lychnidis*.

A total of 26 larvae were collected and individually reared on semisynthetic diet at 23-25°C. During laboratory rearing seven larvae showed the typical signs of lethal polyhedrosis disease. It was difficult to identify the host species; we believe these insects were *L. oleracea*, but we cannot exclude the possibility that one or more of them was *A. lychnidis*. Consequently, we will describe the novel isolate as Lacanobia oleracea NPV (LaolNPV) in the following text.

These isolates were compared with previously characterized viruses from our virus collections, namely Chrysodeixis chalcites nucleopolyhedrovirus (ChchSNPV [21]), Helicoverpa armigera nucleopololyhedrovirus (HearNPV [22]), Spodoptera littoralis nucleopolyhedrovirus (SpliNPV [23]), Spodoptera exigua multiple nucleopolyhedrovirus (SeMNPV [24]), Mamestra brassicae multiple nucleopolyhedrovirus (MbMNPV from the Mamestrin bioinsecticide [25]), Mamestra configurata NPV A (MacoNPV-A [14]) and Mamestra configurata NPV B (MacoNPV-B [15]). Once the phylogenetic relationships had been determined, the novel virus isolate was also compared at the biological level with closely related viruses.

The noctuid larvae used in the present study were *Mamestra brassicae*, *Chrysodeixis chalcites*, *Spodoptera littoralis*, *S. frugiperda*, *S. exigua* and *H. armigera*. These were obtained from laboratory colonies reared at the Universidad Pública de Navarra (UPNA) at 25±1°C, 70±5% relative humidity and 16:8 h day:night photoperiod on a semi-synthetic diet [26]. Bioassays involving *Trichoplusia ni* and *Mamestra configurata*, were performed in the Saskatoon Research and Development Centre, Agriculture and Agri-Food Canada at 21±1°C, 60±5% relative humidity and 16:8 h day:night photoperiod using a semi-synthetic diet [27].

### Molecular characterization

**OB purification, DNA extraction and restriction endonuclease analysis.** Occlusion bodies (OBs) of the different isolates were extracted from dead larvae by homogenizing the cadavers in 500 µl of 0.1% (wt/vol) sodium dodecyl sulfate (SDS), filtered through muslin and centrifuged at 2,500 x g for 5 min. Pellets were resuspended twice in 500 µl of 0.1% SDS and centrifuged for 5 min at 2,500 x g. The resulting pellets were washed twice in distilled water, resuspended in ~200 µl distilled water and OBs were stored at 4°C until required.

For DNA extraction, virions were released from OBs by dissolving the polyhedrin matrix by mixing 100 µl of purified OB suspensions, comprising ~10⁸ OBs, with 100 µl of 0.5 M Na₂CO₃ and 50 µl of 10% SDS in a final volume of 500 µl followed by incubation at 60°C for 10 min. Undissolved OBs and other debris were pelleted at 6,000 x g for 5 min. The virion-containing supernatant was transferred to sterile 1.5 ml vials and incubated at 50°C with 25 µl proteinase K (20 mg/ml) for 1 h. Viral DNA was extracted twice with phenol (pH 7.8):chloroform
(1:1). DNA was precipitated by addition of 10% (v/v) 3 M sodium acetate (pH 5.2) and 2.5 volumes of ice-cold 96% ethanol at 12,000 x g for 10 min. The DNA pellet was washed with 70% cold ethanol and centrifuged for 5 min. The ethanol was discarded and DNA was dried at room temperature for 5 min. Finally, DNA pellets were resuspended in 50–100 µl of 0.1X TE buffer (10 mM Tris, 1 mM EDTA), and kept at 4°C until use. For comparison, DNA was also extracted from OBs of ChchNPV, HearNPV, SpliNPV, SeMNPV, MbMNPV, MacoNPV-A and MacoNPV-B.

For restriction endonuclease (REN) digestion, a sample of viral DNA (1–2 µg) was incubated with one of the following enzymes BamHI, BglII, EcoRI or PstI (Takara) (10 U) at 37°C for 4–12 h. Each reaction was stopped by the addition of 4 µl of 6x loading buffer (0.25% w/v bromophenol blue and 40% w/v sucrose). Fragments were separated by electrophoresis using horizontal 1% agarose gels in TAE buffer (0.004 M Tris-acetate, 0.001 M EDTA, pH 8.0) at 16 V for 14 h. DNA fragments were stained with ethidium bromide and visualized on a UV transilluminator (GeneSnap, Syngene). DNA fragment sizes were estimated by comparison to a standard molecular weight marker (Tandem ladder, Lonza, Rockland, USA).

**Nucleocapsid packaging.** To determine whether the virus under study was a single or multiple nucleocapsid NPV, occlusion derived virions (ODVs) were harvested by treating 5x10⁸ OBs with an equal volume of 0.1 M Na₂CO₃, for 30 min at 28°C. Undissolved OBs and other debris were pelleted and discarded by low-speed centrifugation (2,500 x g, 2 min). The resulting suspensions containing the ODVs were placed on the top of continuous 30–60% (w/w) sucrose gradients and centrifuged at 30,000 x g for 1 h at 4°C in a Beckman Ti28 rotor. The banding pattern was visually inspected, photographed and compared with other viruses.

**Polypeptide analysis.** Structural polypeptides from the purified OBs and occlusion derived virions (ODVs) were analyzed on 11% SDS-polyacrylamide slab gels (SDS-PAGE). ODVs were harvested from OBs by mixing 10 µl of purified OBs at 10¹⁰ OBs/ml with an equal volume of 0.1 M Na₂CO₃ and incubating at 28°C for 30 min. Undissolved OBs and other debris were pelleted at 6,000 x g for 5 min. A 10 µl volume of purified OBs at a concentration of 10¹⁰ OBs/ml and 10 µl of harvested ODVs were solubilized with an equal volume of 2x Laemmli sample buffer (65.8 mM Tris-HCl, pH 6.8, 2.1% SDS, 26.3% (w/v) glycerol, 0.01% bromophenol blue, BioRad) by heating to 100°C for 5 min prior to electrophoresis. Electrophoresis was performed at 50 mA during 2 h. Finally, the gels were stained in Comassie Brilliant Blue R solution (0.1% Comassie R Brilliant Blue, 10% v/v acetic acid and 50% v/v ethanol) for 30 min and destained with a bleaching solution (9.45% v/v ethanol and 6.75% v/v acetic acid). The polypeptide profiles were visually inspected and photographed.

**DNA cloning, sequence analysis and phylogenetic analysis.** In order to obtain genomic sequence information four genomic fragments were terminally sequenced. For this, two genomic libraries were constructed in pUC19 (New England Biolabs) using the EcoRI or PstI digested DNAs. Briefly, 2 µg of viral DNA extracted from purified OBs were digested with EcoRI or PstI overnight at 37°C and then heated to 65°C during 15 min to inactivate enzymes. A 5 µg sample of pUC19 vector was digested with EcoRI or PstI in the same conditions. The digested vector (20 µl) was dephosphorylated during 2h at 37°C using an alkaline phosphatase (Roche Life Science, Basel, Switzerland). After dephosphorylation, the DNA (100 µl) was purified in low melting agarose gel using commercial kit for gel extraction (PCR clean-up, gel extraction, Macherey-Nagel, Düren, Germany) following manufacture's instructions. A 50 ng sample of the purified vector was ligated with 5 µl of the digested viral DNA (100 ng/µl of total fragmented DNA) using T4 DNA ligase (New England Biolabs) at 16°C overnight. After ligation, DH5α cells were transformed with the recombinant plasmids and plated on LB agar containing 100 µg/ml ampicillin, 1 µM IPTG and 80 µg/ml X-Gal. A total of 50 white colonies for
each ligation were amplified in LB broth containing 100 µg/ml ampicillin. Plasmid DNAs were purified by alkaline lysis and screened for the presence of EcoRI or PstI inserts by digestion with the respective enzymes followed by electrophoresis in 1% agarose gel. Inserts were authenticated by comparing their migration in agarose gels with the fragments of the viral DNA generated by the digestion with the same enzymes.

Two EcoRI and two PstI fragments (EcoRI-4Kb, EcoRI-6kb, PstI-5kb and PstI-6kb) were selected for terminal sequencing. Terminal nucleotide sequences were determined by Sanger Sequencing method performed by Stab Vida Company (Caparica, Portugal), employing standard M13 forward and M13 reverse primers. Sequence information was analyzed for the presence of open reading frames (ORFs) and for domain prediction using Clone Manager 9.0 (Scientific and Educational Software Server). Homology searches were performed both at the nucleotide and deduced amino acid levels, for all putative ORFs. DNA and protein comparisons with entries in the updated GenBank/EMBL, SWISS-PROT and PIR databases were performed using BLATn, BLASTp and FASTA programs [28, 29].

For phylogenetic analyses, the DNA sequences within the coding regions of three highly conserved genes, the late expression factor 8 (lef-8), late expression factor 9 (lef-9) and polyhedrin (polh), were used [4, 30, 31]. As recommended by Jehle et al. [4], for viruses with Kimura two-parameter (K-2-P) distances of more than 0.015 in the marker genes, the complete sequences of these three marker genes were determined and used in phylogenetic analyses. For this, primers were designed in the polh, lef-8 and lef-9 genes based on the genomic sequence of the most homologous NPV. The primer pairs used were MacoA-polh-F (ATGTATACCCGTTATAGTTA)-MacoA-polh-R (TTAGTAAGCCGGTCCGTTGTA), MacoA-lef8-F (ATGACGGACGTGATTGACGA)-MacoA-lef8-R (TCATCGAACCACTGTGTTGTG), and MacoA-lef9-F (ATGACCTTTAGCGGTCATTC)-MacoA-lef9-R (CTAGTCCAAAAACATGTCGA). The resulting fragments were cloned into pGEM-T Easy vector (Promega) following manufacturer's instructions and transformed as previously described. Two clones for each gene were selected and nucleotide sequences were determined by Sanger Sequencing method performed by Stab Vida Company (Caparica, Portugal), employing standard SP6 and T7 reverse primers.

The K-2-P distances were calculated for single and concatenated lef-8, lef-9 and polh genes. Concatenated sequences from the same virus were treated as a single sequence. Multiple sequence alignments were performed using MEGA6.06 software [32], and the K-2-P nucleotide substitution model was used for the analysis. Maximum parsimony (MP) phylogenetic trees (1,000 bootstrap replicates) were inferred from the nucleotide sequence alignments using MEGA 6.06. Introduced gaps were treated as missing data.

Genome sequences used in the comparative analysis were obtained from GenBank (accession number included as well as the publication): Autographa californica (Ac) MNPV (NC_001623/L22858 [33]), Agrotis ipsilon (Agip) NPV (NC_011345 [34]), Agrotis segetum (Agse) NPV (NC_007921 [35]), Bombyx mori (Bm) NPV (NC_001962 [36]), Chrysodeixis chalcites (Chch) NPV (NC_007151 [37]), Helicoverpa armigera multiple (HearM) NPV (NC_011615 [38]), Helicoverpa armigera single (HearS) NPV (AF303045 [39]), Helicoverpa zea (Hz) NPV (AF334030 [40]), Lymantria dispar (Ld) MNPV (AF081810 [41]), Mamestra brassicae (Mb) MNPV (NC_023681/JQ798165 [42]), Mamestra configurata A 90/2 (MacoA 90/2) NPV (U59461/AF467808 [14]), Mamestra configurata A 90/4 (MacoA 90/4) NPV (AF539999 [16]), Mamestra configurata B (MacoB) NPV (AY126275 [15]), Rachiplusia ou (Ro) MNPV (AY145471 [43]), Spodoptera exigua (Se) MNPV (AF169823 [44]), Spodoptera frugiperda (Sf) MNPV (HM595733 [45]), Spodoptera litura (Splt) NPV (NC_003102/AF325155 [46]) and Trichoplusia ni (Tn) NPV (NC_007383/DQ017380 [47]).
## Biological characterization

**Host range determination.** The host range of LaolNPV was determined by oral inoculation bioassays with the following lepidopteran species: *M. brassicae*, *M. configurata*, *T. ni*, *C. chalcites*, *S. littoralis*, *S. frugiperda*, *S. exigua* and *H. armigera*. The results of the oral inoculation studies were compared with those of MbMNPV, MacoNPV-A and MacoNPV-B. For this, second-instar larvae from the laboratory colonies were starved for 8–12 h at 26°C and then allowed to drink from an aqueous suspension containing 10% sucrose (w/v) and Fluorella blue and OBs at concentrations of 10⁵ and 10⁷ OBs/ml, following the droplet feeding method [48]. These concentrations were selected to permit determination of the susceptible host species based on preliminary tests (data not shown) and previous studies using MacoNPV [13, 15, 16] or MbMNPV [49, 50]. Larvae that ingested the suspension within 10 min, which were identified by the blue color of their intestine, were transferred to individual wells of a 25-well tissue culture plate with a cube (1 cm³) of semisynthetic diet [26]. Host range tests were performed using 25 larvae per virus concentration and 25 larvae treated with aqueous suspension containing sucrose and Fluorella blue without OBs as controls. Each assay comprised three replicates. Larvae were reared individually at 25°C and mortality was recorded every 24 h until the insects had either died or pupated. Dead larvae were examined microscopically to determine the presence of OBs. When present, OBs were purified from each group of dead larvae and viral DNA was extracted as previously described and subjected to REN analysis to determine the identity of the virus. The percentage of larval mortality was calculated for each virus and concentration and subjected to univariate analysis of variance (ANOVA) in SPSS v23 (IBM SPSS Statistics, Softtonic). The significance of treatments was determined for comparisons among the estimated means by Tukey test (p<0.05).

## Results

### Molecular characterization

**REN profiles indicated the proximity of LaolNPV to MacoNPV-A.** The BglII, EcoRI and PstI profiles of the isolates from each of the seven insects showed very similar patterns with the different enzymes tested. No submolar bands were observed in the different REN profiles, indicating that each isolate comprised a single or majority genotype (Fig 1A). Therefore, the first isolate was selected as the prototype and its DNA and OBs were used for subsequent studies.

When comparing the REN profiles with those of ChchNPV, HearNPV, SeMNPV or SpliNPV, marked differences were observed, suggesting that the novel virus was genetically distinct from the other virus species tested (Fig 1B). Additionally comparison with the REN profiles generated *in silico* with AcMNPV, LdMNPV or TnSNPV also indicated that LaolNPV was genetically distinct from these viruses (data not shown).

The REN profiles were then compared with those of MbMNPV, MacoNPV-A and MacoNPV-B. As observed in the BglII profiles (Fig 1C) the LaolNPV isolate from France might be more similar to MacoNPV-A than to MbMNPV or MacoNPV-B, although clear differences among these viruses were observed in the EcoRI and PstI profiles, suggesting that LaolNPV might represent a novel virus species.

**LaoINPV is a multiple nucleocapsid NPV.** The banding pattern observed following ODV density gradient centrifugation revealed that the ODVs of LaolNPV contained multiple nucleocapsids, as indicated by the multiple bands visible in the sucrose gradient (Fig 2). Based on visual inspection, the banding patterns of MbMNPV, MacoNPV-A and MacoNPV-B were quite similar, and comprised between 1 and 9 nucleocapsids/ODV in MbMNPV and
<figure><img src="image_5.png" /><figcaption>**Fig 1.** Restriction endonuclease analysis of LaoINPV DNA. Restriction endonuclease profiles following digestion of A) genomic DNAs of the seven isolates of LaoINPV collected in alfalfa crops nearby Montpellier with BgIII, EcoRI and PstI enzymes, B) genomic DNAs of Chrysodeixis chalcites NPV (Chch), Helicoverpa armigera SNPV (Hear), Spodoptera</figcaption></figure>
littoralis (Spli), Spodopera exigua MNPV (Se) and the virus under study (LaoINPV) with BamHI, BglII, EcoRI and PstI, and C) genomic DNAs of Mamestra brassicae MNPV (Mb), Mamestra configurata NPV-A (MacoA), Mamestra configurata NPV-B (MacoB) and LaoINPV with BglII, EcoRI and PstI.

https://doi.org/10.1371/journal.pone.0176171.g001

MacoNPV-B compared to between 1 and 8 nucleocapsids/ODV in MacoNPV-A and LaoINPV.

**OB and ODV structural polypeptides in LaolNPV.** The SDS-PAGE profiles of OB and ODV polypeptides were generally similar between LaolNPV and the other three NPVs (Fig 3). There were, however, clear differences in the molecular weights of specific major proteins and the presence or absence of particular proteins. In the OB polypeptide profiles, the virus under study presented a band at 75 kDa that was absent in MbMNPV, MacoNPV-A and MacoNPV-B. Additionally, the intensity of the band at 20 kDa was much reduced in LaolNPV compared to that of the other three viruses. Similarly, in the ODV profiles, LaolNPV presented a clear band of ~30 kDa that was absent in the other three viruses, while the band below the 25 kDa marker showed a lower molecular weight and the band of ~10 kDa had a stronger intensity in LaolNPV than in the other three viruses.

**Terminal sequencing and phylogenetic analysis indicates that LaolNPV evolved from the same ancestor as MacoNPV-A.** The terminal sequence information obtained from the four cloned fragments showed homology with 18 ORFs. However, although they appeared distinct in agarose gels, the two PstI fragments that were terminal sequenced were in fact identical fragments (corresponding to the 5,575 bp [nt 6,124 to 11,699] PstI fragment in the MacoNPV-A 90/2 genome [14]), showing the same ORFs. Therefore, just a total of 14 different ORFs were identified as being homologues of known baculovirus ORFs and all were assigned on the basis of their similarity to previously identified ORFs in the MacoNPV-A 90/2 genome [14]. The direction of transcription, its position in the sequenced fragment relative to that in the MacoNPV-A 90/2 genome, as well as its identity to MacoNPV-A 90/2 homologs are summarized in Table 1. Within the 5' end of PstI-5kb (M13-Rv sequence) two ORFs were identified, namely *orf5* and *odpv-6e* that showed 93% and 95% nt sequence identity with the homologous ORFs of MacoNPV-A 90/2, while the 3' end (M13-Fw sequence) included *efp* and *orf10*, with 97 and 99% nt identity to the corresponding homologs in MacoNPV-A 90/2, respectively (Table 1). Similarly, within the PstI-6kb 5'-terminal sequence (M13-Fw) *orf5* and *odpv-6e* were

<figure><img src="image_6.png" /><figcaption>**Fig 2. ODV banding patterns.** ODV banding patterns of MbMNPV, MacoNPV-A, MacoNPV-B and LaoINPV after continuous sucrose gradient separation. Asterisks indicate the position of the different visible bands.</figcaption></figure>

https://doi.org/10.1371/journal.pone.0176171.g002
<figure><img src="image_7.png" /><figcaption>**Fig 3.** **OB and ODV structural polypeptides patterns.** OB and ODV structural polypeptides profiles of MbMNPV (Mb), MacoNPV-A (MacoA), MacoNPV-B (MacoB) and the virus under study (LaoINPV). OBs and harvested ODVs after solubilization were electrophoresed on an 11% SDS-PAGE gel. The molecular weights of the fragments of the Precision Plus Protein All Blue marker, from Bio-Rad, (M) are indicated above each band. The arrows indicated the differences observed in the presence and absence of proteins fragments and in differences in molecular weights.</figcaption></figure>

https://doi.org/10.1371/journal.pone.0176171.g003

identified, showing 93 and 96% identity to homologous ORFs in the MacoNPV-A 90/2 genome, respectively, while the 3'-terminal sequence (M13-Rv) included *efp* and *orf10* ORFs, showing 95 and 99% identity. EcoRI-4kb is the 3,899 bp fragment with homology to the MacoNPV-A 90/2 genome region located between nt 41,533 to 45,533. The M13-Fw sequence included *orf44*, *pkip-1* and *orf46* with 94, 95 and 97% of identity, respectively. The M13-Rv sequence included *pif2* and *pif1*, both with identities of 96%. Finally, the EcoRI-6kb is the 5,818 bp fragment homologous to nt 97,657 to 97,474 in the MacoNPV-A 90/2 genome. Within the M13-Fw sequence the *cg30* and *vp91* genes were located with identities of 95 and 98%, respectively. The M13-Rv sequence contained the *gp41*, *orf105* and *vlf-1* genes with identities of 98, 97 and 93%, to MacoNPV-A 90/2 homologues, respectively (Table 1).

The amino acid sequence identity and similarity of the 14 ORFs were also compared with those of other NPVs (Table 2). The most homologous NPV to LaoINPV was MacoNPV-A 90/2 with 95.93% identity and 96.64% similarity at the amino acid level, followed by MacoNPV-A 90/4 with 95.78% identity and 96.50% similarity. The identity and similarity values to MbMNPV (90.00 and 93.79%) and MacoNPV-B (89.78 and 93.64%), were lower, but were of a similar magnitude to the identity and similarity values (90.42 and 95.00%) for comparisons of MacoNPV-A 90/2 and MacoNPV-B or MbMNPV (data not shown). In contrast, the identity
<table><caption>**Table 1. ORFs identified in the terminal sequences of four restriction fragments.** The position and orientation of the 14 putative ORFs in the most homologous baculovirus MacoNPV-A 90/2 genome [14] along with the identity at nucleotide level and the percentage of gaps are also indicated.</caption><thead><tr><th rowspan="2">ORF num</th><th rowspan="2">Gene family</th><th rowspan="2">Size (nt)</th><th rowspan="2">Genomic fragment</th><th rowspan="2">Primer</th><th rowspan="2">Sequence size (nt)</th><th rowspan="2">Position in the fragment (nt)</th><th rowspan="2">Most homologous ORF</th><th rowspan="2">Position in MacoNPV-A 90/2</th><th colspan="2">MacoNPV-A 90/2</th><th colspan="2">MacoNPV-A 90/4</th></tr><tr><th>Identities</th><th>Gaps</th><th>Identities</th><th>Gaps</th></tr></thead><tbody><tr><td>1</td><td>orf5</td><td>176</td><td>PstI-5kb</td><td>M13-Rv</td><td>1,187</td><td>9&gt;185</td><td>Maco5</td><td>6,124&gt;6,312</td><td>175/189 (93%)</td><td>12/189 (6%)</td><td>175/189 (93%)</td><td>12/189 (6%)</td></tr><tr><td>2</td><td>odpv-6e</td><td>864</td><td>PstI-5kb</td><td>M13-Rv</td><td>1,187</td><td>266&gt;1,129</td><td>Maco6</td><td>6,400&gt;7,263</td><td>825/864 (95%)</td><td>0/864 (0%)</td><td>822/864 (95%)</td><td>3/864 (0%)</td></tr><tr><td>3</td><td>efp</td><td>845</td><td>PstI-5kb</td><td>M13-Fw</td><td>1,179</td><td>30&gt;874</td><td>Maco9</td><td>10,566&gt;11,410</td><td>819/845 (97%)</td><td>0/845 (0%)</td><td>821/845 (97%)</td><td>0/845 (0%)</td></tr><tr><td>4</td><td>orf10</td><td>189</td><td>PstI-5kb</td><td>M13-Fw</td><td>1,179</td><td>972&lt;1,160</td><td>Maco10</td><td>11,516&lt;11,704</td><td>187/189 (99%)</td><td>0/189 (0%)</td><td>187/189 (99%)</td><td>0/189 (0%)</td></tr><tr><td>1</td><td>orf5</td><td>177</td><td>PstI-6kb</td><td>M13-Fw</td><td>1,085</td><td>17&gt;193</td><td>Maco5</td><td>6,124&gt;6,312</td><td>175/189 (93%)</td><td>12/189 (6%)</td><td>175/189 (93%)</td><td>12/189 (6%)</td></tr><tr><td>2</td><td>odpv-6e</td><td>812</td><td>PstI-6kb</td><td>M13-Fw</td><td>1,085</td><td>274&gt;1085</td><td>Maco6</td><td>6,400&gt;7,211</td><td>777/812 (96%)</td><td>0/812 (0%)</td><td>774/812 (95%)</td><td>3/812 (0%)</td></tr><tr><td>3</td><td>efp</td><td>458</td><td>PstI-6kb</td><td>M13-Rv</td><td>1,051</td><td>1&gt;458</td><td>Maco9</td><td>10,953&gt;11,410</td><td>437/458 (95%)</td><td>0/458 (0%)</td><td>437/458 (95%)</td><td>0/458 (0%)</td></tr><tr><td>4</td><td>orf10</td><td>190</td><td>PstI-6kb</td><td>M13-Rv</td><td>1,051</td><td>556&lt;745</td><td>Maco10</td><td>11,516&lt;11,705</td><td>188/190 (99%)</td><td>0/190 (0%)</td><td>188/190 (99%)</td><td>0/190 (0%)</td></tr><tr><td>5</td><td>orf44</td><td>345</td><td>EcoRI-4kb</td><td>M13-Fw</td><td>1,177</td><td>21&gt;365</td><td>Maco44</td><td>41,562&gt;41,894</td><td>323/345 (94%)</td><td>12/345 (3%)</td><td>323/345 (94%)</td><td>12/345 (3%)</td></tr><tr><td>6</td><td>pkip-1</td><td>510</td><td>EcoRI-4kb</td><td>M13-Fw</td><td>1,177</td><td>388&gt;897</td><td>Maco45</td><td>41,917&gt;42,427</td><td>488/511 (95%)</td><td>1/511 (0%)</td><td>488/510 (96%)</td><td>0/510 (0%)</td></tr><tr><td>7</td><td>orf46</td><td>259</td><td>EcoRI-4kb</td><td>M13-Fw</td><td>1,177</td><td>919&lt;1,177</td><td>Maco46</td><td>42,449&lt;42,707</td><td>252/259 (97%)</td><td>0/259 (0%)</td><td>252/259 (97%)</td><td>0/259 (0%)</td></tr><tr><td>8</td><td>pif2</td><td>144</td><td>EcoRI-4kb</td><td>M13-Rv</td><td>1,011</td><td>3&gt;146</td><td>Maco48</td><td>44,540&gt;44,683</td><td>138/144 (96%)</td><td>0/144 (0%)</td><td>138/144 (96%)</td><td>0/144 (0%)</td></tr><tr><td>9</td><td>pif1</td><td>739</td><td>EcoRI-4kb</td><td>M13-Rv</td><td>1,011</td><td>161&gt;899</td><td>Maco49</td><td>44,698&gt;45,436</td><td>707/739 (96%)</td><td>0/739 (0%)</td><td>708/739 (96%)</td><td>0/739 (0%)</td></tr><tr><td>10</td><td>cg30</td><td>611</td><td>EcoRI-6kb</td><td>M13-Fw</td><td>1,206</td><td>19&gt;617</td><td>Maco100</td><td>91,689&gt;92,299</td><td>582/611 (95%)</td><td>12/611 (1%)</td><td>578/611 (95%)</td><td>12/611 (1%)</td></tr><tr><td>11</td><td>vp91</td><td>450</td><td>EcoRI-6kb</td><td>M13-Fw</td><td>1,206</td><td>679&lt;1,128</td><td>Maco101</td><td>92,662&lt;92,811</td><td>441/450 (98%)</td><td>0/450 (0%)</td><td>435/444 (98%)</td><td>0/444 (0%)</td></tr><tr><td>12</td><td>gp41</td><td>510</td><td>EcoRI-6kb</td><td>M13-Rv</td><td>1,203</td><td>2&gt;511</td><td>Maco104</td><td>96,354&gt;96,863</td><td>501/510 (98%)</td><td>0/510 (0%)</td><td>501/510 (98%)</td><td>0/510 (0%)</td></tr><tr><td>13</td><td>orf105</td><td>345</td><td>EcoRI-6kb</td><td>M13-Rv</td><td>1,203</td><td>508&gt;852</td><td>Maco105</td><td>96,860&gt;97,063</td><td>337/346 (97%)</td><td>1/346 (0%)</td><td>337/345 (98%)</td><td>0/345 (0%)</td></tr><tr><td>14</td><td>vlf-1</td><td>273</td><td>EcoRI-6kb</td><td>M13-Rv</td><td>1,203</td><td>854–1126</td><td>Maco106</td><td>97,207&gt;97,479</td><td>255/273 (93%)</td><td>2/273 (0%)</td><td>270/273 (99%)</td><td>0/273 (0%)</td></tr></tbody><tfoot><tr><td colspan="12">Mean</td><td>95.89%</td><td>0.88%</td><td>96.28%</td><td>0.88%</td></tr></tfoot></table>

https://doi.org/10.1371/journal.pone.0176171.t001

and similarity values for comparisons of the MacoNPV-A 90/2 and MacoNPV-A 90/4 isolates, using the same sequences, increased to 99.53 and 99.74%, respectively (data not shown). Similar results were obtained at the nucleotide level: the nt sequence identity between LaolNPV and the MacoNPV-A isolates was intermediate (95.89%) between that found between the same species (isolates of MacoNPV-A; 99.68%), and different species (MacoNPV-A 90/2 and MacoNPV-B or MbMNPV; 90.42 and 88.95%, respectively).

To further determine the relationship between LaolNPV and the other NPVs and to establish whether this is a novel virus or a variant of MacoNPV-A, a phylogenetic study was performed using a total of 19 NPV sequences and the complete nucleotide sequences of single or concatenated *lef-8*, *lef-9* and *polh* using maximum parsimony (MP) and MEGA6.06 distance
<table><caption>Table 2. (Continued)</caption><thead><tr><th rowspan="3">Protein family</th><th rowspan="3">Size (aa)</th><th colspan="10">Protein position (aa)</th></tr><tr><th colspan="2">Identity (%)</th><th colspan="2">Positive (%)</th><th colspan="5">Gaps (%)</th></tr><tr><th>MacoNPV-A (90/2)</th><th>MacoNPV-A (90/4)</th><th>MacoNPV-B</th><th>MbNPV</th><th>Heart/NPV</th><th>AgsNPV</th><th>STNPV</th><th>ChchNPV</th><th>ToSNPV</th><th>LdNPV</th><th>AdINPV</th></tr></thead><tbody><tr><td rowspan="2">VP91</td><td>194/202 (96%)</td><td>195/202 (96%)</td><td>171/206 (83%)</td><td>173/206 (83%)</td><td>171/206 (83%)</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td>4/202 (96%)</td><td>4/202 (96%)</td><td>17/206 (8%)</td><td>17/206 (8%)</td><td>17/206 (8%)</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td rowspan="2">GP41</td><td>661 to 809 aa</td><td>662 to 810 aa</td><td>664 to 812 aa</td><td>664 to 812 aa</td><td>664 to 812 aa</td><td>630 to 778 aa</td><td>672 to 813 aa</td><td>645 to 769 aa</td><td>675 to 819 aa</td><td>702 to 849 aa</td><td>428 to 571 aa</td></tr><tr><td>146/149 (98%)</td><td>146/149 (98%)</td><td>147/149 (99%)</td><td>147/149 (99%)</td><td>146/149 (98%)</td><td>94/150 (63%)</td><td>92/143 (64%)</td><td>75/145 (52%)</td><td>78/1457 (53%)</td><td>72/149 (48%)</td><td>72/145 (50%)</td></tr><tr><td rowspan="2">ORF105</td><td>148/149 (99%)</td><td>148/149 (99%)</td><td>149/149 (100%)</td><td>149/149 (100%)</td><td>148/149 (98%)</td><td>116/150 (77%)</td><td>121/143 (84%)</td><td>166/145 (73%)</td><td>104/147 (70%)</td><td>99/149 (66%)</td><td>106/145 (73%)</td></tr><tr><td>0/149 (0%)</td><td>0/149 (0%)</td><td>0/149 (0%)</td><td>0/149 (0%)</td><td>0/149 (0%)</td><td>3/150 (2%)</td><td>1/143 (0%)</td><td>3/145 (2%)</td><td>6/147 (4%)</td><td>4/149 (2%)</td><td>4/145 (2%)</td></tr><tr><td rowspan="2">VLF-1</td><td>165 to 333 aa</td><td>165 to 333 aa</td><td>165 to 333 aa</td><td>165 to 333 aa</td><td>165 to 333 aa</td><td>174 to 332 aa</td><td>171 to 329 aa</td><td>140 to 297 aa</td><td>164 to 321 aa</td><td>158 to 318 aa</td><td>232 to 387 aa</td></tr><tr><td>169/169 (100%)</td><td>169/169 (100%)</td><td>167/169 (99%)</td><td>167/169 (9%)</td><td>167/169 (9%)</td><td>130/159 (62%)</td><td>130/159 (62%)</td><td>97/161 (60%)</td><td>98/161 (61%)</td><td>93/165 (56%)</td><td>92/160 (58%)</td></tr><tr><td rowspan="2">Mean Identity</td><td>169/169 (100%)</td><td>169/169 (100%)</td><td>167/169 (98%)</td><td>167/169 (98%)</td><td>167/169 (98%)</td><td>147/159 (92%)</td><td>147/159 (92%)</td><td>120/161 (74%)</td><td>122/161 (75%)</td><td>118/165 (71%)</td><td>118/160 (73%)</td></tr><tr><td>0/169 (0%)</td><td>0/169 (0%)</td><td>0/169 (0%)</td><td>0/169 (0%)</td><td>0/159 (0%)</td><td>0/159 (0%)</td><td>0/159 (0%)</td><td>5/161 (3%)</td><td>5/161 (3%)</td><td>9/165 (5%)</td><td>7/160 (4%)</td></tr><tr><td rowspan="2">Mean Similarity</td><td>1 to 65 aa</td><td>1 to 114 aa</td><td>40 to 151 aa</td><td>40 to 151 aa</td><td>1 to 124 aa</td><td>1 to 108 aa</td><td>1 to 116 aa</td><td>1 to 119 aa</td><td>9 to 102 aa</td><td>1 to 92 aa</td><td>39/104 (38%)</td></tr><tr><td>64/65 (98%)</td><td>112/114 (98%)</td><td>105/114 (92%)</td><td>105/114 (92%)</td><td>63/124 (51%)</td><td>61/109 (56%)</td><td>51/118 (43%)</td><td>50/121 (41%)</td><td>49/104 (47%)</td><td>38/104 (38%)</td><td>57/104 (54%)</td></tr><tr><td rowspan="2">Mean Identity</td><td>64/65 (98%)</td><td>112/114 (98%)</td><td>110/114 (96%)</td><td>110/114 (96%)</td><td>79/124 (63%)</td><td>72/109 (66%)</td><td>79/118 (66%)</td><td>72/121 (61%)</td><td>68/104 (65%)</td><td>57/104 (54%)</td><td>57/104 (54%)</td></tr><tr><td>0/65 (0%)</td><td>0/114 (0%)</td><td>2/114 (1%)</td><td>2/114 (1%)</td><td>2/114 (1%)</td><td>11/124 (8%)</td><td>3/109 (2%)</td><td>11/118 (9%)</td><td>14/121 (11%)</td><td>11/104 (10%)</td><td>12/104 (11%)</td></tr><tr><td rowspan="2">Mean Similarity</td><td>1 to 91</td><td>1 to 91</td><td>1 to 91 aa</td><td>1 to 91 aa</td><td>1 to 91 aa</td><td>1 to 91 aa</td><td>1 to 91 aa</td><td>2 to 90 aa</td><td>2 to 90 aa</td><td>21 to 100 aa</td><td>5 to 88 aa</td></tr><tr><td>89/91 (98%)</td><td>89/91 (98%)</td><td>87/91 (96%)</td><td>87/91 (96%)</td><td>87/91 (96%)</td><td>85/94 (90%)</td><td>85/91 (93%)</td><td>58/89 (65%)</td><td>57/89 (64%)</td><td>52/80 (65%)</td><td>55/84 (65%)</td><td>68/84 (80%)</td></tr><tr><td rowspan="2">Mean Identity</td><td>89/91 (97%)</td><td>89/91 (97%)</td><td>87/91 (95%)</td><td>87/91 (95%)</td><td>87/91 (95%)</td><td>86/94 (91%)</td><td>87/91 (95%)</td><td>71/189 (79%)</td><td>72/89 (80%)</td><td>68/80 (85%)</td><td>68/84 (80%)</td><td>68/84 (80%)</td></tr><tr><td>0/91 (0%)</td><td>0/91 (0%)</td><td>0/91 (0%)</td><td>0/91 (0%)</td><td>0/91 (0%)</td><td>0/94 (0%)</td><td>0/91 (0%)</td><td>1/89 (1%)</td><td>1/80 (1%)</td><td>49.00%</td><td>50.00%</td><td>65.62%</td></tr><tr><td rowspan="2">Mean Similarity</td><td>95.93%</td><td>95.78%</td><td>89.78%</td><td>90.00%</td><td>90.07%</td><td>61.42%</td><td>61.50%</td><td>51.00%</td><td>52.12%</td><td>49.00%</td><td>50.00%</td><td>65.62%</td></tr><tr><td>96.64%</td><td>96.50%</td><td>93.64%</td><td>93.79%</td><td>93.64%</td><td>73.25%</td><td>75.10%</td><td>68.78%</td><td>69.88%</td><td>65.11%</td><td>65.62%</td><td></td></tr></tbody></table>

https://doi.org/10.1371/journal.pone.0176171.t002
methods [4]. The K-2-P model was used to calculate the genetic distances between the different viruses.

The phylogenetic trees using complete genes sequences showed that in the four cases *Mamestra* NPVs formed a separate branch together with HearMNPV and LaolNPV, although with different bootstrap values between branches depending on the gene (Fig 4). MbMNPV shares close sequence identity with MacoNPV-B and HearMNPV, whereas MacoNPV-A strains formed a separate branch. In all cases, except the *lef-8* gene, LaolNPV clustered on the same branch as MacoNPV-A isolates, clearly separated from the MbMNPV/MacoNPV-B/HaMNPV complex. MacoNPV-A 90/2 and MacoNPV-A 90/4 invariably clustered together, whereas LaolNPV, although clearly proximal, was separated on another branch.

The K-2-P values between LaolNPV and the MbMNPV/MacoNPV-B/HearMNPV complex for *lef-9*, *polh* and concatenated sequences were greater than 0.05 in all cases (underlined in Table 3 and Table 4). Similarly, the K-2-P values for MacoNPV-A and MacoNPV-B, considered as distinct species, also exceeded 0.05 (dashed line underlined in Table 3 and Table 4). In contrast, K-2-P distances between the two viruses previously considered to be the same species, MbMNPV and MacoNPV-B, were consistent below 0.015 (shown in bold in Table 3 and Table 4). In contrast, when LaolNPV was compared with MacoNPV-A strains the K-2-P distances were intermediate, between 0.05 and 0.015, being of 0.049, 0.028, 0.026 and 0.039 for *lef-8*, *lef-9*, *polh* and the concatenated sequences, respectively (in bold and underlined in Table 3 and Table 4). In this case other characteristics, such as the host range need to be evaluated to establish species status.

## Biological characterization

**LaoINPV has a narrow host range compared to MacoNPV-A.** The host range of LaolNPV was compared with that of MacoNPV-A, MacoNPV-B and MbMNPV at high (1x10⁷ OBs/ml) (Fig 5A) and moderate (1x10⁵ OBs/ml) (Fig 5B) concentrations in eight lepidopteran species; *M. brassicae*, *M. configurata*, *T. ni*, *C. chalcites*, *S. littoralis*, *S. exigua*, *S. frugiperda* and *H. armigera*. In all cases dead larvae died due to the inoculated virus as showed in the REN profiles obtained with the OBs from cadavers (data not shown). The REN profile results indicate that cross-contamination did not occur during host range bioassays and OB treatments did not activate sublethal infections by homologous NPVs.

Mortality differed significantly among the insect species at high (F<sub>7,159</sub> = 114.026; p<0.0001) and moderate (F<sub>7,159</sub> = 74.955; p<0.0001) OB concentrations. Mortality also varied significantly according to identity of the virus at high (F<sub>3,159</sub> = 393.261; p<0.0001) and moderate (F<sub>3,159</sub> = 78.882; p<0.0001) OB concentrations.

MacoNPV-A OBs produced variable mortalities that differed significantly across the different host species (F<sub>7,32</sub> = 60.253; p<0.0001). The lower concentration of MacoNPV-A OB produced mortalities between 0%, in the least susceptible species such as *S. littoralis* or *H. armigera*, to 12-72% mortality in the other species. MacoNPV-A was significantly more pathogenic (in terms of the mortality caused by a particular concentration of OBs) to *M. brassicae* (72% mortality) than to its homologous species, *M. configurata* (12% mortality) (Tukey test p<0.05). At the higher OB concentration MacoNPV-A also caused mortality in *S. littoralis* (17%) and *H. armigera* (14%), and mortality differed significantly among the different species examined (F<sub>7,32</sub> = 42.279; p<0.0001). Similarly, at the highest concentration MacoNPV-A was more pathogenic to *M. brassicae* (90%) than to *M. configurata* (59%) (Tukey test p<0.05), followed by *S. exigua* (73%) *C. chalcites* (72%), and *T. ni* (58%) (Tukey test p>0.05) and was least pathogenic to *S. frugiperda* (44%).

MacoNPV-B also showed significant variation in pathogenicity to the different species (F<sub>7,32</sub> = 55.436; p<0.0001). MacoNPV-B was more pathogenic to *M. brassicae* (85%), *T. ni* (55%).
<figure><img src="image_8.png" /><figcaption>**Fig 4. Phylogenetic analysis.** Molecular phylogenetic analysis by maximum likelihood method obtained with the complete gene sequences of single or concatenated *lef-8*, *lef-9*, and *polh* genes. The phylogeny was inferred using the maximum likelihood method based on the Kimura two-parameter model. The tree with the highest log likelihood value (-41803.7062) is shown. The percentage of trees in which the associated taxa clustered together is shown next to the branches. Initial tree(s) for the heuristic search were obtained automatically by applying Neighbor-Join and BioNJ algorithms to a matrix of pairwise distances estimated using the Maximum Composite Likelihood (MCL) approach, and then selecting the topology with the superior log likelihood value. The tree is drawn to scale, with branch lengths measured in the number of substitutions per site (shown as values next to the branches). The analysis involved 19 nucleotide sequences. All positions containing gaps and missing data were eliminated. There were a total of 4,702 positions in the final dataset. Analyses were conducted in MEGA6.</figcaption></figure>

https://doi.org/10.1371/journal.pone.0176171.g004
<table><caption>**Table 3. K-2-P values.** Pairwise distances of the nucleotide sequences of *lef*-9 and *lef*- complete genes. The distances of the virus under study (LaoINPV) to MbMNPV/MacoNPV-B/ HearMNPV complex are underlined, the distances to MacoNPV-A strains are in given in bold and underlined. The K-2-P values for MacoNPV-A and MacoNPV-B, two different species, are underlined with a dashed line, while those values between MbMNPV and MacoNPV-B, two viruses of the same viral species, are in bold.</caption><thead><tr><th></th><th>1</th><th>2</th><th>3</th><th>4</th><th>5</th><th>6</th><th>7</th><th>8</th><th>9</th><th>10</th><th>11</th><th>12</th><th>13</th><th>14</th><th>15</th><th>16</th><th>17</th><th>18</th><th>19</th></tr></thead><tbody><tr><td><u>lef-9</u></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>1. AcMNPV</td><td>-</td><td>0.503</td><td>0.492</td><td>0.044</td><td>0.461</td><td>0.457</td><td>0.487</td><td>0.483</td><td>0.492</td><td>0.457</td><td>0.468</td><td>0.468</td><td>0.455</td><td>0.470</td><td>0.005</td><td>0.468</td><td>0.488</td><td>0.464</td><td>0.428</td></tr><tr><td>2. AgipNPV</td><td>0.455</td><td>-</td><td>0.190</td><td>0.500</td><td>0.369</td><td>0.282</td><td>0.446</td><td>0.445</td><td>0.353</td><td>0.280</td><td>0.254</td><td>0.254</td><td>0.283</td><td>0.262</td><td>0.503</td><td>0.235</td><td>0.255</td><td>0.213</td><td>0.379</td></tr><tr><td>3. AgseNPV</td><td>0.457</td><td>0.193</td><td>-</td><td>0.482</td><td>0.369</td><td>0.278</td><td>0.468</td><td>0.469</td><td>0.382</td><td>0.276</td><td>0.277</td><td>0.276</td><td>0.275</td><td>0.285</td><td>0.485</td><td>0.220</td><td>0.244</td><td>0.218</td><td>0.388</td></tr><tr><td>4. BmMNPV</td><td>0.029</td><td>0.466</td><td>0.457</td><td>-</td><td>0.450</td><td>0.448</td><td>0.461</td><td>0.458</td><td>0.478</td><td>0.450</td><td>0.458</td><td>0.458</td><td>0.450</td><td>0.456</td><td>0.045</td><td>0.468</td><td>0.486</td><td>0.458</td><td>0.430</td></tr><tr><td>5. ChchSNPV</td><td>0.471</td><td>0.415</td><td>0.415</td><td>0.462</td><td>-</td><td>0.363</td><td>0.403</td><td>0.406</td><td>0.482</td><td>0.361</td><td>0.355</td><td>0.356</td><td>0.362</td><td>0.357</td><td>0.456</td><td>0.362</td><td>0.367</td><td>0.369</td><td>0.110</td></tr><tr><td>6. HearMNPV</td><td>0.458</td><td>0.315</td><td>0.318</td><td>0.456</td><td>0.395</td><td>-</td><td>0.422</td><td>0.426</td><td>0.421</td><td>0.011</td><td>0.092</td><td>0.092</td><td>0.008</td><td>0.095</td><td>0.545</td><td>0.309</td><td>0.289</td><td>0.288</td><td>0.375</td></tr><tr><td>7. HearNPV</td><td>0.453</td><td>0.387</td><td>0.395</td><td>0.456</td><td>0.403</td><td>0.401</td><td>-</td><td>0.005</td><td>0.456</td><td>0.416</td><td>0.411</td><td>0.411</td><td>0.420</td><td>0.413</td><td>0.484</td><td>0.438</td><td>0.422</td><td>0.444</td><td>0.406</td></tr><tr><td>8. HzSNPV</td><td>0.448</td><td>0.386</td><td>0.393</td><td>0.452</td><td>0.398</td><td>0.397</td><td>0.010</td><td>-</td><td>0.456</td><td>0.420</td><td>0.415</td><td>0.415</td><td>0.424</td><td>0.417</td><td>0.480</td><td>0.437</td><td>0.425</td><td>0.443</td><td>0.408</td></tr><tr><td>9. LdMNPV</td><td>0.512</td><td>0.329</td><td>0.374</td><td>0.525</td><td>0.484</td><td>0.477</td><td>0.474</td><td>0.468</td><td>-</td><td>0.419</td><td>0.413</td><td>0.412</td><td>0.420</td><td>0.412</td><td>0.487</td><td>0.400</td><td>0.447</td><td>0.392</td><td>0.477</td></tr><tr><td>10. MbMNPV</td><td>0.458</td><td>0.317</td><td>0.317</td><td>0.455</td><td>0.398</td><td>0.006</td><td>0.404</td><td>0.400</td><td>0.478</td><td>-</td><td>0.091</td><td>0.091</td><td><strong>0.005</strong></td><td>0.093</td><td>0.454</td><td>0.304</td><td>0.286</td><td>0.284</td><td>0.371</td></tr><tr><td>11. MacoNPV-A 90/2</td><td>0.460</td><td>0.302</td><td>0.309</td><td>0.462</td><td>0.396</td><td>0.090</td><td>0.406</td><td>0.403</td><td>0.460</td><td>0.089</td><td>-</td><td>0.003</td><td>0.089</td><td><strong>0.028</strong></td><td>0.465</td><td>0.290</td><td>0.295</td><td>0.279</td><td>0.353</td></tr><tr><td>12. MacoNPV A 90/4</td><td>0.458</td><td>0.302</td><td>0.309</td><td>0.460</td><td>0.396</td><td>0.090</td><td>0.407</td><td>0.403</td><td>0.457</td><td>0.090</td><td>0.003</td><td>-</td><td>0.089</td><td><strong>0.029</strong></td><td>0.465</td><td>0.293</td><td>0.297</td><td>0.280</td><td>0.353</td></tr><tr><td>13. MacoNPV-B</td><td>0.460</td><td>0.313</td><td>0.317</td><td>0.458</td><td>0.398</td><td>0.005</td><td>0.400</td><td>0.398</td><td>0.477</td><td>0.051</td><td><strong>0.049</strong></td><td>0.091</td><td>0.091</td><td>-</td><td>0.092</td><td>0.453</td><td>0.306</td><td>0.285</td><td>0.287</td><td>0.370</td></tr><tr><td>14. LaoINPV</td><td>0.464</td><td>0.311</td><td>0.308</td><td>0.468</td><td>0.395</td><td>0.055</td><td>0.392</td><td>0.388</td><td>0.457</td><td>0.051</td><td><strong>0.049</strong></td><td>0.091</td><td>0.091</td><td>-</td><td>0.467</td><td>0.292</td><td>0.295</td><td>0.295</td><td>0.279</td><td>0.279</td></tr><tr><td>15. RoMNPV</td><td>0.040</td><td>0.479</td><td>0.467</td><td>0.061</td><td>0.471</td><td>0.462</td><td>0.463</td><td>0.457</td><td>0.524</td><td>0.463</td><td>0.469</td><td>0.468</td><td>0.466</td><td>0.472</td><td>-</td><td>0.468</td><td>0.490</td><td>0.464</td><td>0.427</td><td>0.427</td></tr><tr><td>16. SeMNPV</td><td>0.454</td><td>0.252</td><td>0.264</td><td>0.465</td><td>0.424</td><td>0.346</td><td>0.404</td><td>0.400</td><td>0.407</td><td>0.348</td><td>0.328</td><td>0.328</td><td>0.345</td><td>0.336</td><td>0.454</td><td>-</td><td>0.193</td><td>0.124</td><td>0.366</td><td>0.366</td></tr><tr><td>17. SfMNPV</td><td>0.462</td><td>0.312</td><td>0.308</td><td>0.464</td><td>0.402</td><td>0.357</td><td>0.407</td><td>0.409</td><td>0.471</td><td>0.359</td><td>0.353</td><td>0.353</td><td>0.355</td><td>0.353</td><td>0.458</td><td>0.234</td><td>-</td><td>0.194</td><td>0.372</td><td>0.372</td></tr><tr><td>18. SplitNPV</td><td>0.462</td><td>0.237</td><td>0.272</td><td>0.479</td><td>0.429</td><td>0.349</td><td>0.400</td><td>0.396</td><td>0.389</td><td>0.354</td><td>0.353</td><td>0.333</td><td>0.351</td><td>0.340</td><td>0.470</td><td>0.132</td><td>0.211</td><td>-</td><td>0.368</td><td>0.368</td></tr><tr><td>19. TnSNPV</td><td>0.471</td><td>0.406</td><td>0.396</td><td>0.463</td><td>0.127</td><td>0.391</td><td>0.402</td><td>0.400</td><td>0.476</td><td>0.395</td><td>0.391</td><td>0.393</td><td>0.393</td><td>0.397</td><td>0.465</td><td>0.404</td><td>0.386</td><td>0.413</td><td>-</td><td>-</td></tr></tbody></table>

https://doi.org/10.1371/journal.pone.0176171.t003
<table><caption>**Table 4. K-2-P values.** Pairwise distances of the nucleotide sequences of polh complete gene and concatenated polh/lef-8/lef-9 genes. The distances of the virus under study (LaoINPV) to MbMNPV/MacoNPV-B/HeaMNPV complex are underlined, the distances to MacoNPV-A strains are in given in bold and underlined. The K-2-P values for MacoNPV-A and MacoNPV-B, two different species, are underlined with a dashed line, while those values between MbMNPV and MacoNPV-B, two viruses of the same viral species, are in bold.</caption><thead><tr><th>polh<br>polh/lef-8/lef-9</th><th>1</th><th>2</th><th>3</th><th>4</th><th>5</th><th>6</th><th>7</th><th>8</th><th>9</th><th>10</th><th>11</th><th>12</th><th>13</th><th>14</th><th>15</th><th>16</th><th>17</th><th>18</th><th>19</th></tr></thead><tbody><tr><td>1. AcMNPV</td><td>-</td><td>0.243</td><td>0.266</td><td>0.285</td><td>0.211</td><td>0.218</td><td>0.331</td><td>0.337</td><td>0.300</td><td>0.218</td><td>0.199</td><td>0.199</td><td>0.199</td><td>0.220</td><td>0.198</td><td>0.230</td><td>0.286</td><td>0.236</td><td>0.259</td><td>0.186</td></tr><tr><td>2. AgipNPV</td><td>0.432</td><td>-</td><td>0.163</td><td>0.281</td><td>0.238</td><td>0.196</td><td>0.325</td><td>0.318</td><td>0.308</td><td>0.230</td><td>0.189</td><td>0.189</td><td>0.189</td><td>0.203</td><td>0.183</td><td>0.294</td><td>0.183</td><td>0.180</td><td>0.169</td><td>0.246</td></tr><tr><td>3. AgseNPV</td><td>0.434</td><td>0.188</td><td>-</td><td>0.303</td><td>0.248</td><td>0.227</td><td>0.325</td><td>0.320</td><td>0.317</td><td>0.237</td><td>0.212</td><td>0.212</td><td>0.212</td><td>0.233</td><td>0.209</td><td>0.322</td><td>0.204</td><td>0.209</td><td>0.209</td><td>0.269</td></tr><tr><td>4. BmMNPV</td><td>0.064</td><td>0.444</td><td>0.440</td><td>-</td><td>0.304</td><td>0.293</td><td>0.358</td><td>0.353</td><td>0.315</td><td>0.298</td><td>0.302</td><td>0.302</td><td>0.294</td><td>0.292</td><td>0.080</td><td>0.307</td><td>0.304</td><td>0.298</td><td>0.303</td><td>0.303</td></tr><tr><td>5. ChchSNPV</td><td>0.419</td><td>0.374</td><td>0.374</td><td>0.430</td><td>-</td><td>0.225</td><td>0.280</td><td>0.281</td><td>0.319</td><td>0.227</td><td>0.231</td><td>0.231</td><td>0.227</td><td>0.230</td><td>0.298</td><td>0.263</td><td>0.213</td><td>0.235</td><td>0.113</td><td>0.113</td></tr><tr><td>6. HearMNPV</td><td>0.415</td><td>0.284</td><td>0.287</td><td>0.424</td><td>0.356</td><td>-</td><td>0.301</td><td>0.295</td><td>0.300</td><td>0.011</td><td>0.080</td><td>0.080</td><td>0.008</td><td>0.072</td><td>0.304</td><td>0.208</td><td>0.200</td><td>0.219</td><td>0.208</td><td>0.208</td></tr><tr><td>7. HearNPV</td><td>0.441</td><td>0.396</td><td>0.401</td><td>0.440</td><td>0.384</td><td>0.390</td><td>-</td><td>0.012</td><td>0.393</td><td>0.296</td><td>0.315</td><td>0.315</td><td>0.296</td><td>0.317</td><td>0.354</td><td>0.290</td><td>0.308</td><td>0.333</td><td>0.291</td><td>0.291</td></tr><tr><td>8. HzSNPV</td><td>0.437</td><td>0.394</td><td>0.400</td><td>0.435</td><td>0.382</td><td>0.388</td><td>0.009</td><td>-</td><td>0.390</td><td>0.289</td><td>0.309</td><td>0.309</td><td>0.290</td><td>0.311</td><td>0.347</td><td>0.299</td><td>0.306</td><td>0.340</td><td>0.291</td><td>0.291</td></tr><tr><td>9. LdMNPV</td><td>0.470</td><td>0.334</td><td>0.365</td><td>0.474</td><td>0.458</td><td>0.432</td><td>0.455</td><td>0.452</td><td>-</td><td>0.308</td><td>0.276</td><td>0.276</td><td>0.304</td><td>0.278</td><td>0.331</td><td>0.334</td><td>0.285</td><td>0.280</td><td>0.300</td><td>0.300</td></tr><tr><td>10. MbMNPV</td><td>0.415</td><td>0.286</td><td>0.288</td><td>0.426</td><td>0.358</td><td>0.008</td><td>0.389</td><td>0.387</td><td>0.433</td><td>-</td><td>0.089</td><td>0.089</td><td><strong>0.003</strong></td><td>0.079</td><td>0.304</td><td>0.207</td><td>0.205</td><td>0.227</td><td>0.207</td><td>0.207</td></tr><tr><td>11. MacoNPV-A 90/2</td><td>0.415</td><td>0.269</td><td>0.282</td><td>0.433</td><td>0.357</td><td>0.088</td><td>0.394</td><td>0.392</td><td>0.414</td><td>0.089</td><td>-</td><td>0.000</td><td>0.089</td><td><strong>0.026</strong></td><td>0.302</td><td>0.207</td><td>0.205</td><td>0.227</td><td>0.207</td><td>0.207</td></tr><tr><td>12. MacoNPV A 90/4</td><td>0.414</td><td>0.269</td><td>0.281</td><td>0.432</td><td>0.357</td><td>0.088</td><td>0.395</td><td>0.393</td><td>0.412</td><td>0.089</td><td>0.002</td><td>-</td><td>0.089</td><td><strong>0.026</strong></td><td>0.302</td><td>0.207</td><td>0.187</td><td>0.204</td><td>0.218</td><td>0.218</td></tr><tr><td>13. MacoNPV-B</td><td>0.416</td><td>0.285</td><td>0.287</td><td>0.426</td><td>0.358</td><td>0.006</td><td>0.389</td><td>0.387</td><td>0.433</td><td><strong>0.004</strong></td><td>0.089</td><td>0.089</td><td>-</td><td>0.080</td><td>0.300</td><td>0.208</td><td>0.206</td><td>0.227</td><td>0.210</td><td>0.210</td></tr><tr><td>14. LaoINPV</td><td>0.417</td><td>0.274</td><td>0.284</td><td>0.433</td><td>0.356</td><td>0.069</td><td>0.387</td><td>0.385</td><td>0.412</td><td>0.067</td><td><strong>0.039</strong></td><td><strong>0.039</strong></td><td>0.068</td><td>-</td><td>0.298</td><td>0.220</td><td>0.187</td><td>0.217</td><td>0.210</td><td>0.210</td></tr><tr><td>15. RoMNPV</td><td>0.055</td><td>0.453</td><td>0.449</td><td>0.058</td><td>0.436</td><td>0.432</td><td>0.448</td><td>0.442</td><td>0.480</td><td>0.432</td><td>0.438</td><td>0.437</td><td>0.432</td><td>0.439</td><td>-</td><td>0.333</td><td>0.300</td><td>0.305</td><td>0.294</td><td>0.294</td></tr><tr><td>16. SeMNPV</td><td>0.428</td><td>0.235</td><td>0.239</td><td>0.438</td><td>0.380</td><td>0.310</td><td>0.397</td><td>0.395</td><td>0.393</td><td>0.310</td><td>0.294</td><td>0.295</td><td>0.309</td><td>0.302</td><td>0.436</td><td>-</td><td>0.156</td><td>0.095</td><td>0.272</td><td>0.272</td></tr><tr><td>17. SfMNPV</td><td>0.430</td><td>0.272</td><td>0.272</td><td>0.443</td><td>0.359</td><td>0.308</td><td>0.395</td><td>0.397</td><td>0.431</td><td>0.309</td><td>0.304</td><td>0.305</td><td>0.307</td><td>0.305</td><td>0.441</td><td>0.208</td><td>-</td><td>0.152</td><td>0.232</td><td>0.232</td></tr><tr><td>18. SplINPV</td><td>0.428</td><td>0.219</td><td>0.244</td><td>0.442</td><td>0.377</td><td>0.308</td><td>0.403</td><td>0.401</td><td>0.370</td><td>0.311</td><td>0.295</td><td>0.295</td><td>0.310</td><td>0.301</td><td>0.301</td><td>0.440</td><td>0.123</td><td>0.195</td><td>-</td><td>0.251</td></tr><tr><td>19. TnSNPV</td><td>0.405</td><td>0.373</td><td>0.372</td><td>0.425</td><td>0.120</td><td>0.354</td><td>0.386</td><td>0.386</td><td>0.448</td><td>0.355</td><td>0.350</td><td>0.351</td><td>0.354</td><td>0.354</td><td>0.424</td><td>0.372</td><td>0.356</td><td>0.374</td><td>-</td><td>-</td></tr></tbody></table>

https://doi.org/10.1371/journal.pone.0176171.t004
A)

<figure><img src="image_9.png" /></figure>

B)

<figure><img src="image_10.png" /><figcaption>**Fig 5.** Host range of LaoINPV. Percentage of mortality produced by MacoNPV-A, MacoNPV-B, MbMNPV and LaoINPV OBs inoculated at concentrations of A) 1x10⁵ OBs/ml (≈33 OBs/larva) and B) 1x10⁷ OBs/ml (≈3,300 OBs/larva) against *M. brassicae*,</figcaption></figure>
*M. configurata*, *T. ni*, *C. chalcites*, *S. littoralis*, *S. exigua*, *S. frugiperda* and *H. armigera* larvae. Vertical bars indicate standard errors. Values above bars indicate means. Values followed by different letters are significantly different (Tukey test p<0.05).

https://doi.org/10.1371/journal.pone.0176171.g005

(67%), *C. chalcites* (69%) and *S. exigua* (38%) than to the homologous species *M. configurata* (11%) (Tukey test p<0.05). Similarly, at the moderate OB concentration MacoNPV-B did not cause lethal disease in *S. littoralis* or *H. armigera*, whereas at the higher OB concentration, MacoNPV-B caused mortality in both *S. littoralis* (19%) and *H. armigera* (23%). MacoNPV-B was markedly more pathogenic to heterologous hosts *M. brassicae* (100%), *T. ni* (88%), *C. chalcites* (88%) and *S. exigua* (92%), than to its homologous host *M. configurata* (62%).

Although the mortality across species produced by MbMNPV was quite similar to those produced by MacoNPV-A and MacoNPV-B at high and moderate OB concentrations, MbMNPV showed a greater host range as it was pathogenic at the moderate OB concentration to species such as *H. armigera*. In general, MbMNPV was more pathogenic to the homologous host *M. brassicae* at high (100%) and moderate (89%) concentrations than to the heterologous hosts (Tukey, p<0.05). At the moderate OB concentration MbMNPV produced 100% mortality in *M. brassicae* and mortalities lower than 100% in all heterologous species, although these differences were not significant in the case of *C. chalcites* or *S. exigua* (Tukey, p>0.05).

LaolNPV showed a markedly different host range compared with MacoNPV-A. LaolNPV was not infective to *M. brassicae* or *M. configurata* even at the high OB concentration. LaolNPV was only pathogenic to *C. chalcites* at moderate and high concentrations, which resulted in 19 and 50% mortality, respectively. The host range of LaolNPV was completely different, therefore, to that of MacoNPV-A, indicating that these viruses occupy distinct ecological niches. Regrettably, we were unable to determine the pathogenicity of LaolNPV OBs in *L. oleracea* or *A. australis* as laboratory or field populations have not been available for testing.

## Discussion

The present study aimed to characterize a novel NPV from diseased larvae of *L. oleracea* collected in 2011 from a field of alfalfa near Montpellier, France, although because of difficulties in differentiation of noctuid species in the larval stage we could not be certain of the identity of each and every diseased larva collected.

Restriction endonuclease analyses suggested that the baculovirus under study, LaolNPV, and MacoNPV-A diverged to a greater degree than MbMNPV and MacoNPV-B, that were previously considered to be the same viral species isolated from different hosts [4, 15]. The differences in REN profiles were of a similar magnitude to those found between two different viruses, namely MacoNPV-A and MacoNPV-B. Sucrose density gradients demonstrated that, like MbMNPV, MacoNPV-A and MacoNPV-B [16], LaolNPV is a multiple nucleocapsid virus with a banding pattern quite similar to that of MacoNPV-A, suggesting a similar distribution of multiple nucleocapsids within the ODVs. Additionally, OB and ODV structural polypeptide profiles were also similar to MbMNPV, MacoNPV-A and MacoNPV-B, but with clear differences in the presence and molecular weights of several proteins. These results were consistent with the hypothesis that LaolNPV was closely related to, but distinct from the other NPVs that infect *Mamestra* spp.

Terminal sequencing information revealed that the most homologous NPV to the virus under study was MacoNPV-A 90/2 [14] closely followed by MacoNPV-A 90/4 [16]. The percentage of amino acid sequence identity between LaolNPV and MacoNPV-A was intermediate to that found between strains of MacoNPV-A or MbMNPV/MacoNPV-B, viruses of the same viral species, or between viruses of different species, such as MacoNPV-A and MbMNPV/MacoNPV-B. When comparing with MbMNPV, the identity clearly decreased to a values
corresponding to different species. Consequently, it appears that LaolNPV is more similar to the New World MacoNPV-A than to MbMNPV, a species from the Old World (Europe).

Phylogenetic analysis demonstrated that LaolNPV formed a separate branch within a MacoNPV-A clade, with high bootstrap support, as a distinct lineage from MbMNPV, MacoNPV-B and HearMNPV. Additionally, the baculovirus under study differed with respect to the MbMNPV/MacoNPV-B/HaMNPV group with a K-2-P distance greater than 0.05 in all cases, supporting the concept that LaolNPV is phylogenetically different from MbMNPV, HaMNPV or MacoNPV-B. In contrast, the MbMNPV/MacoNPV-B/HaMNPV group had K-2-P distances among the different viruses of less than 0.015.

In many cases the same virus isolated from different hosts has been assigned different names. This was the case with *Anagrapha falcifera* NPV (AnfaNPV) and *Rachioplusia ou* MNPV (RoMNPV) [43, 51] and HearSNPV and HzSNPV [40]. Alternatively, different viruses isolated from the same species have received the same name [52]. Although both MacoNPV-A and MacoNPV-B were isolated from the same host and their genomes are closely related, MacoNPV-B represents a separately evolving virus from MacoNPV-A, as the K-2-P distances estimated in the present study were consistently higher than 0.05 (0.088), and exceed the genetic distance between RoMNPV and BmNPV (0.055), which are also thought to be different species [4, 43]. MacoNPV-A is much more prevalent and infective to *M. configurata* populations than MacoNPV-B, however MacoNPV-B has higher pathogenicity to a wider range of host species. It is likely that there are alternative hosts for MacoNPV-B, so that infections in *M. configurata* are only commonly observed during outbreak years of high population density. It is intriguing that these two viruses could evolve divergently while infecting the same host, *M. configurata* [4, 14, 15]. This type of situation highlights the present drawbacks in baculovirus nomenclature as virus names are based on the host from which the virus was first isolated. However, it is also clear that host range can prove to be an informative aspect among the ecological characteristics that contribute to the baculovirus species definition criteria.

Baculovirus isolates from separate populations of a particular host species, or from closely related species, often differ in pathogenicity and the origin of a particular NPV is not a useful predictor of its virulence in a specific host, although local isolates tend to be more effective in controlling local insect populations than geographically distant isolates [6, 13, 21, 53]. Surprisingly, although the baculovirus under study was phylogenetically close to MacoNPV-A, LaolNPV showed particular host range differences and was not pathogenic to *M. brassicae* or *M. configurata*, even at moderately high OB concentrations. In contrast, MacoNPV-A produced mortality across a variety of species tested. Therefore the effective host range of MacoNPV-A appears to be notably broader than that of LaolNPV. This suggests that these viruses occupy different ecological niches, and therefore, in combination with the phylogenetic analyses, are likely to represent different species. Surprisingly, the divergence in pathogenicity between viruses within related insect populations was also observed with the MacoNPV viruses. MacoNPV-B was more pathogenic than MacoNPV-A for several species. Additionally both MacoNPV-A and -B viruses were more pathogenic to *M. brassicae* than to the homologous host, *M. configurata*. This surprised us, as NPVs tend to be more pathogenic to homologous hosts rather than heterologous species [21, 54, 55]. However, in the present study we observed that MbMNPV was more pathogenic to *M. configurata* than MacoNPV-A or -B, which was also reported by Erlandson [6]. MbMNPV and MacoNPV-B represent variants of the same viral species, with distributions that include southern Canada and Europe [15]. Both viruses infect *M. configurata*, and as MbMNPV is common in Europe, it is possible that MacoNPV-B-like viruses may be found worldwide.

We conclude that *Lacanobia oleracea* *nucleopolyhedrovirus* (LaolNPV) should be considered as a new species of multiple nucleocapsid NPV in the *Alphabaculovirus* genus on the basis of REN profiles, gene sequence features and host range properties, in addition to the original
insect host species and geographical point of isolation. Future studies should examine intraspecific variation in LaolNPV and establish the dose or concentration mortality response of this virus in its homologous host and sympatric noctuid species, such as *Agrochola lychnidis*, that are present in the alfalfa agroecosystem of southern Europe. Such information could prove valuable in evaluating the virus' potential as a biological insecticide for control of *L. oleracea* and related noctuid pests in this region.

## Acknowledgments

We thank Noelia Gorría and Itxaso Ibáñez for insect rearing and Leire Bardaji and Eduardo Aguirre for technical assistance. We thank Doug Baldwin for assistance with bioassays in *M. configurata* and *T. ni*. TW received logistical support from Gabriel Mercado (INECOL). This study received financial support from the National Institute of Agricultural Research from the Spanish Ministry project RTA2013-00114-C02-02 (Dr. Oihane Simón) and from the Department of Industry of the Government of Navarre project IIQ1406-RI1 (Dr. Primitivo Cabalero). This study was supported in part by funding from Agriculture and Agri-Food Canada Research Project J000197 (Dr. Martin A. Erlandson).

## Author Contributions

**Conceptualization:** OS MAE PC NV.

**Data curation:** MF OS.

**Formal analysis:** OS ME TW.

**Funding acquisition:** OS PC MAE DAT.

**Investigation:** OS MF MAE DAT.

**Methodology:** OS PC MAE TW.

**Project administration:** OS PC MAE TW NV.

**Resources:** OS PC MAE NV.

**Supervision:** PC MAE DAT TW NV.

**Validation:** OS MAE.

**Visualization:** OS.

**Writing – original draft:** OS MAE MF.

**Writing – review & editing:** PC MAE DAT TW NV.

## References

1. Herniou EA, Olszewski JA, Cory JS, O'Reilly DR. The genome sequence and evolution of baculoviruses. Annu Rev Entomol. 2003; 48: 211–234. https://doi.org/10.1146/annurev.ento.48.091801.112756 PMID: 12414741

2. Herniou EA, Arif BM, Becnel JJ, Blissard GW, Bonning B, Harrison R, et al. Virus taxonomy: classification and nomenclature of viruses: Ninth Report of the International Committee on Taxonomy of Viruses. Amsterdam: Elsevier; 2012. pp. 167–174.

3. Jehle JA, Blissard GW, Bonning BC, Cory JS, Herniou EA, Rohrmann GF, et al. On the classification and nomenclature of baculoviruses: a proposal for revision. Arch Virol. 2006; 151: 257–266.

4. Jehle JA, Lange M, Wang H, Hu Z, Wang Y, Hauschild R. Molecular identification and phylogenetic analysis of baculoviruses from Lepidoptera. Virology. 2006; 346: 180–193. https://doi.org/10.1016/j.virol.2005.10.032 PMID: 16313938
5. Craveiro SR, Melo FL, Ribeiro ZMA, Ribeiro BM,áo SN, et al. *Pseudoplusia includes* single nucleopo-
lyhedrovirus: genetic diversity, phylogeny and hypervariability of the *pif-2* gene. J Invertebr Pathol.
2013; 114: 258–267. https://doi.org/10.1016/j.jip.2013.08.005 PMID: 24012501

6. Erlandson MA. Genetic variation in field populations of baculoviruses: mechanisms for generating variation and its potential role in baculovirus epizootiology. Virol Sin. 2009; 24: 458-469.

7. Harrison RL, Keena MA, Rowley DL. Classification, genetic variation and pathogenicity of *Lymantria dispar* nucleopolyhedrovirus isolates from Asia, Europe and North America. J Invertebr Pathol. 2014; 116: 27–35. https://doi.org/10.1016/j.jip.2013.12.005 PMID: 24370838

8. Cory JS, Green BM, Paul RK, Hunter-Fujita F. Genotypic and phenotypic diversity of a baculovirus pop-
ulation within an individual insect host. J Invertebr Pathol. 2005; 89: 101–111. https://doi.org/10.1016/j.
jip.2005.03.008 PMID: 15876438

9. Redman EM, Wilson K, Grzywacz D, Cory JS. High levels of genetic diversity in *Spodoptera exempta* NPV from Tanzania. J Invertebr Pathol. 2010; 105: 190–193. https://doi.org/10.1016/j.jip.2010.06.008 PMID: 20600096

10. Redman EM, Wilson K, Cory JS. Trade-offs and mixed infections in an obligate-killing insect pathogen. J Anim Ecol. 2016; 85: 1200–12009. https://doi.org/10.1111/1365-2656.12547 PMID: 27155461

11. Erlandson MA. *Mamestra configurata* Walker, bertha armyworm (Lepidoptera: Noctuidae). In: Mason PG, Gillespie DR editors. Biological Control Programmes in Canada 2001–2012. Wallingford: CABI; 2013. pp. 228–237.

12. Mason PG, Arthur AP, Olfert OO, Erlandson MA. The bertha armyworm (*Mamestra configurata*) (Lepidoptera: Noctuidae) in western Canada. Can Entomol. 1998; 130: 321–336.

13. Erlandson MA. Biological and biochemical comparison of *Mamestra configurata* and *Mamestra brassicae* nuclear polyhedrosis virus isolates pathogenic for the bertha armyworm, *Mamestra configuarata* (Lepidoptera: Noctuidae). J Invertebr Pathol. 1990; 56: 47–56.

14. Li Q, Donly C, Li L, Willis LG, Theilmann DA, Erlandson MA. Sequence organization of the *Mamestra configurata* nucleopolyhedrovirus genome. Virology. 2002; 294: 106–121. https://doi.org/10.1006/viro. 2001.1313 PMID: 11886270

15. Li L, Donly C, Li Q, Willis LG, Keddie BA, Erlandson MA, et al. Identification and genomic analysis of a second species of nucleopolyhedrovirus isolated from *Mamestra configurata*. Virology. 2002; 297: 226–244. PMID: 12083822

16. Li L, Li Q, Willis LG, Erlandson M, Theilmann DA, Donly C. Complete comparative genomic analysis of two field isolates of *Mamestra configurata* nucleopolyhedrovirus-A. J Gen Virol. 2005; 86: 91–105. https://doi.org/10.1099/vir.0.80488-0 PMID: 15604435

17. Herniou EA, Olszewski JA, O'Reilly DR, Cory JS. Ancient coevolution of baculoviruses and their insect hosts. J Virol. 2004; 78: 3244–3251. https://doi.org/10.1128/JVI.78.7.3244-3251.2004 PMID: 15016845

18. King AMQ, Adams MJ, Carstens EB, Lefkowitz EJ. Virus taxonomy: ninth report of the International Committee on Taxonomy of Viruses. Amsterdam: Elsevier; 2012.

19. Kimura M, A simple method for estimating evolutionary rates of base substitutions through comparative studies of nucleotide sequences. J Mol Evol. 1980; 16: 111–120. PMID: 7463489

20. Cory JS, Myers JH. The ecology and evolution of insect baculoviruses. Annu Rev Ecol Evol Syst. 2003; 34: 239–272.

21. Bernal A, Carnero A, Hernández-Suárez E, Williams T. Caballero P, Simón O. A native variant of *Chrysodeixis chalcites* nucleopolyhedrovirus: the basis for a promising bioinsecticide for control of *C. chalcites* in Canary Islands' banana crops. Biol. Control 2013; 67: 101–110.

22. Arrizubeta M, Williams T, Caballero P, Simón O. Selection of a nucleopolyhedrovirus isolate from *Helicoverpa armigera* as the basis for a biological insecticide. Pest Manag Sci. 2013; 70: 967–976. https://doi.org/10.1002/ps.3637 PMID: 23983128

23. Kikhno I, Gutierrez S, Croizier L, Croizier G, López-Ferber M. Characterization of *pif*, a gene required for the *per* os infectivity of *Spodoptera littoralis* nucleopolyhedrovirus. J. Gen Virol. 2002; 83: 3013–3022. https://doi.org/10.1099/0022-1317-83-12-3013 PMID: 12466478

24. Muñoz D, Murillo R, Krell PJ, Vlak JM. Caballero P. Four genotypic variants of a *Spodoptera exigua* nucleopolyhedrovirus (Se-SP2) are distinguishable by a hypervariable genomic region. Virus Res. 1999; 59: 61–74. PMID: 10854166

25. Murillo R, Muñoz D, Lipa JJ, Caballero P. Biochemical characterization of three nucleopolyhedrovirus isolates of *Spodoptera exigua* and *Mamestra brassicae*. J Appl Ent. 2001; 125: 267–270.

26. Greene GL, Leppa NC, Dickerson WA. Velvetbean caterpillar: a rearing procedure and artificial medium. J. Econ. Entomol. 1976; 69, 487–488.
27. Bucher DE, Bracken GK. The bertha armyworm, *Mamestra configurata* (Lepidoptera: Noctuidae). Artificial diet and rearing technique. Can Entomol. 1976; 108: 1327-1338.

28. Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ. Basic local alignment search tool. J Mol Biol. 1990; 215: 403-410. https://doi.org/10.1016/S0022-2836(05)80360-2 PMID: 2231712

29. Pearson WR. Rapid and sensitive sequence comparison with FASTP and FASTA. Methods Enzymol. 1990; 183: 63-98. PMID: 2156132

30. Eberle K, Sayed S, Rezapanah M, Shojai-Estabragh S, Jehle J. Diversity and evolution of the *Cydia pomonella* granulovirus. J Gen Virol. 2009; 90: 662-671. https://doi.org/10.1099/vir.0.006999-0 PMID: 19218212

31. Lange M, Wang H, Zhihong H, Jehle J. Towards a molecular identification and classification system of lepidopteran-specific baculoviruses. Virology. 2004; 325: 36-47. https://doi.org/10.1016/j.virol.2004.04.023 PMID: 15231384

32. Tamura K, Stecher G, Peterson D, Filipski A, Kumar S. MEGA6: Molecular Evolutionary Genetics Analysis version 6.0. Mol Biol Evol. 2013; 30: 2725-2729. https://doi.org/10.1093/molbev/mst197 PMID: 24132122

33. Ayres MD, Howard SC, Kuzio J, López-Ferber M, Possee RD. The complete DNA sequence of *Autographa californica* nuclear polyhedrosis virus. Virology. 1994; 202: 586-605. https://doi.org/10.1006/viro.1994.1380 PMID: 8030224

34. Harrison RL. Genomic sequence analysis of the Illinois strain of the *Agrotis ipsilon* multiple nucleopolyhedrovirus. Virus Genes. 2009; 38: 155-170. https://doi.org/10.1007/s11262-008-0297-y PMID: 19015970

35. Jakubowska AK, Peters SA, Ziemnicka J, Vlak JM, van Oers MM. Genome sequence of an enhancing gene-rich nucleopolyhedrovirus (NPV) from *Agrotis segetum*: collinearity with *Spodoptera exigua* multiple NPV. J Gen Virol. 2006; 87: 537-551. https://doi.org/10.1099/vir.0.81461-0 PMID: 16476975

36. Gomi S, Majima K, Maeda S. Sequence analysis of the genome of *Bombyx mori* nucleopolyhedrovirus. J Gen Virol. 1999; 80: 1323-1337. https://doi.org/10.1099/0022-1317-80-5-1323 PMID: 10355780

37. van Oers MM, Abma-Henkens MH, Herniou EA, de Groot JC, Peters S, Vlak JM. Genome sequence of *Chrysodeixis chalcites* nucleopolyhedrovirus, a baculovirus with two DNA photolyase genes. J Gen Virol. 2005; 86: 2069-2080. https://doi.org/10.1099/vir.0.80964-0 PMID: 15958686

38. Tang P, Zhang H, Li Y, Han B, Wang G, Qin Q, et al. Genomic sequencing and analyses of HearMNPV-a new multinucleocapsid nucleopolyhedrovirus isolated from *Helicoverpa armigera*. Virol J. 2012; 9: 168. https://doi.org/10.1186/1743-422X-9-168 PMID: 22913743

39. Zhang CX, Ma XC, Guo ZJ. Comparison of complete genome sequence between C1 and G4 isolates of the *Helicoverpa armigera* single nucleocapsid nucleopolyhedrovirus. Virology. 2005; 333: 190-199. https://doi.org/10.1016/j.virol.2004.12.028 PMID: 15708604

40. Chen X, Zhang WJ, Wong J, Chun G, Lu A, McCutchen BF, et al. Comparative analysis of the complete genome sequences of *Helicoverpa zea* and *Helicoverpa armigera* single-nucleocapsid nucleopolyhedroviruses. J Gen Virol. 2002; 83: 673-684. https://doi.org/10.1099/0022-1317-83-3-673 PMID: 11842262

41. Kuzio JM, Pearson MN, Harwood SH, Funk CJ, Evans JT, Slavicek JM, et al. Sequence analysis of the genome of a baculovirus pathogenic for *Lymantria dispar*. Virology. 1999; 253: 17-34. https://doi.org/10.1006/viro.1998.9469 PMID: 9887315

42. Choi JB, Heo WI, Shin TY, Bae SM, Kim WJ, Kim JI, et al. Complete genomic sequences and comparative analysis of *Mamestra brassicae* nucleopolyhedrovirus isolated in Korea. Virus Genes. 2013; 47: 133-151. https://doi.org/10.1007/s11262-013-0922-2 PMID: 23712441

43. Harrison RL, Bonning BC. Comparative analysis of the genomes of *Rachiplusia ou* and *Autographa californica* multiple nucleopolyhedroviruses. J Gen Virol. 2003; 84: 1827-1842. https://doi.org/10.1099/vir.0.19146-0 PMID: 12810877

44. Ijkel WF, van Estrien EA, Heldens JG, Broer JG, Zuidema D, Goldbach RW, et al. Sequence and organization of the *Spodoptera exigua* multicapsid nucelopolyhedrovirus genome. J Gen Virol. 1999; 80: 3289-3304. https://doi.org/10.1099/0022-1317-80-12-3289 PMID: 10567663

45. Simón O, Palma L, Beperet I, Muñoz D, López-Ferber M, Caballero P, et al. Sequence comparison between three geographically distinct *Spodoptera frugiperda* multiple nucleopolyhedrovirus isolates: Detecting positively selected genes. J Invertebr Pathol. 2011; 107: 33-42. https://doi.org/10.1016/j.jip.2011.01.002 PMID: 21238456

46. Pang Y, Yu J, Wang L, Hu X, Bao W, Li G, et al. Sequence analysis of the *Spodoptera litura* multicapsid nucelopolyhedrovirus genome. Virology. 2001; 287: 391404.
47. Willis LG, Seipp R, Stewart TM, Erlandson MA, Theilmann DA. Sequence analysis of the complete genome of *Trichoplusia ni* single nucleopolyhedrovirus and the identification of a baculoviral photolyase gene. Virology. 2005; 338: 209–226. https://doi.org/10.1016/j.virol.2005.04.041 PMID: 15951000

48. Hughes PR, Wood HA. A synchronous peroral technique for the bioassay of insect viruses. J Invertebr Pathol. 1981; 37: 154–159.

49. Mukawa S, Goto C. *In vivo* characterization of a group II nucleopolyhedrovirus isolated from *Mamestra brassicae* (Lepidoptera: Noctuidae) in Japan. J Virol. 2006; 87: 1491–1500.

50. Rovesti L, Crook NE, Winstabley D. Biological and biochemical relationships between the nucleopolyhedroviruses of *Mamestra brassicae* and *Heliothis armigera*. J Invertebr Pathol. 2000; 75: 2–8. https://doi.org/10.1006/jipa.1999.4903 PMID: 10631051

51. Harrison RL, Bonning BC. The nucleopolyhedroviruses of *Rachoplusia ou* and *Anagrapha falcifera* are isolates of the same virus. J Gen Virol. 1999; 80: 2793–2798. https://doi.org/10.1099/0022-1317-80-10-2793 PMID: 10573177

52. Kislev N, Edelman M. DNA restriction pattern differences from geographic isolates of *Spodoptera littora-lis* nuclear polyhedrosis virus. Virology. 1982; 119: 219–222. PMID: 18635146

53. Barrera G, Simón O, Villamizar L, Williams T, Caballero P. *Spodoptera frugiperda* multiple nucleopolyhedrovirus as a potential biological insecticide: genetic and phenotypic comparison of field isolates from Colombia. Biol Control. 2011; 58: 113–120.

54. Shapiro M. Enhancement in activity of homologous and heterologous baculoviruses infectious to beet armyworm (Lepidoptera: Noctuidae) by optical brightener. J Econ Entomol. 2000; 93: 572–576. PMID: 10902301

55. Shapiro M, Vaughn JL. Enhancement in activity of homologous and heterologous baculoviruses infectious to cotton bollworm (Lepidoptera: Noctuidae) by optical brightener. J Econ Entomol. 1995; 88: 265–269.