| Literature DB >> 31861567 |
Sa Yang1,2, Philippe Gayral3, Hongxia Zhao4, Yaojun Wu5, Xuejian Jiang5, Yanyan Wu1,2, Diane Bigot3, Xinling Wang1,2, Dahe Yang1,2, Elisabeth A Herniou3, Shuai Deng1,2, Fei Li1,2, Qingyun Diao1,2, Eric Darrouzet3, Chunsheng Hou1,2.
Abstract
Since the discovery that honey bee viruses play a role in colony decline, researchers have made major breakthroughs in understanding viral pathology and infection processes in honey bees. Work on virus transmission patterns and virus vectors, such as the mite Varroa destructor, has prompted intense efforts to manage honey bee health. However, little is known about the occurrence of honey bee viruses in bee predators, such as vespids. In this study, we characterized the occurrence of 11 honey bee viruses in five vespid species and one wasp from four provinces in China and two vespid species from four locations in France. The results showed that all the species from China carried certain honey bee viruses, notably Apis mellifera filamentous virus (AmFV), Deformed wing virus (DWV), and Israeli acute paralysis virus (IAPV); furthermore, in some vespid colonies, more than three different viruses were identified. In France, DWV was the most common virus; Sacbrood virus (SBV) and Black queen cell virus (BQCV) were observed in one and two samples, respectively. Phylogenetic analyses of IAPV and BQCV sequences indicated that most of the IAPV sequences belonged to a single group, while the BQCV sequences belonged to several groups. Additionally, our study is the first to detect Lake Sinai virus (LSV) in a hornet from China. Our findings can guide further research into the origin and transmission of honey bee viruses in Vespidae, a taxon of ecological, and potentially epidemiological, relevance.Entities:
Keywords: AmFV; DWV; IAPV; honeybee viruses; hornets
Mesh:
Year: 2019 PMID: 31861567 PMCID: PMC7019919 DOI: 10.3390/v12010006
Source DB: PubMed Journal: Viruses ISSN: 1999-4915 Impact factor: 5.048
Figure 1Sampling locations in (a) China and (b) France. Sampled species are indicated by names and photos.
Presence of major honey bee viruses in various vespid species.
| Country | Developmental Stage | Colony ID/Province (Species) | AmFV | IAPV | DWV-A | DWV-B | DWV-A/KV | CSBV | SBV | BQCV | LSV |
|---|---|---|---|---|---|---|---|---|---|---|---|
| China | Adult | B1/Beijing ( | nt | nt | + | nt | + | nt | nt | nt | nt |
| Pupa | + | + | + | nt | + | nt | nt | nt | nt | ||
| Adult | B2/Beijing ( | nt | nt | + | nt | nt | nt | nt | + | nt | |
| Pupa | + | + | + | nt | + | + | nt | + | + | ||
| Adult | B3/Beijing ( | nt | nt | nt | nt | nt | nt | nt | nt | nt | |
| Pupa | nt | nt | nt | nt | nt | nt | nt | nt | nt | ||
| Adult | B4/Beijing ( | + | nt | nt | nt | nt | nt | nt | nt | nt | |
| Pupa | + | nt | nt | nt | nt | nt | nt | + | nt | ||
| Adult | B5/Beijing ( | nt | nt | nt | nt | nt | nt | nt | nt | nt | |
| Pupa | + | + | nt | nt | nt | nt | nt | nt | nt | ||
| Adult | Q1/Qinghai | nt | nt | nt | nt | nt | nt | nt | nt | nt | |
| Adult | GD1/Guangdong ( | nt | nt | + | nt | nt | nt | nt | nt | nt | |
| Pupa | nt | + | + | nt | nt | nt | nt | nt | nt | ||
| Adult | GD2/Guangdong ( | nt | nt | nt | nt | nt | nt | nt | nt | nt | |
| Pupa | nt | + | nt | nt | nt | nt | nt | nt | nt | ||
| Adult | GD3/Guangdong ( | nt | nt | nt | nt | nt | nt | nt | nt | nt | |
| Pupa | nt | + | nt | nt | nt | nt | nt | nt | nt | ||
| Adult | GD4/Guangdong ( | nt | nt | nt | nt | nt | nt | nt | nt | nt | |
| Pupa | nt | + | + | + | + | nt | nt | + | nt | ||
| Adult | GD5/Guangdong ( | nt | + | + | nt | nt | nt | nt | nt | nt | |
| Pupa | nt | + | + | nt | nt | nt | nt | nt | nt | ||
| Adult | JX/Jiangxi | nt | nt | nt | nt | nt | nt | nt | nt | nt | |
| Pupa | nt | nt | nt | nt | nt | nt | nt | nt | nt | ||
| France | Adult | Fondettes ( | nt | nt | + | nt | nt | nt | nt | nt | nt |
| Adult | Saché ( | nt | nt | + | nt | nt | nt | nt | nt | nt | |
| Adult | Saché ( | nt | nt | + | nt | nt | nt | nt | nt | nt | |
| Adult | Lussault-sur-Loire ( | nt | nt | + | nt | nt | nt | nt | nt | nt | |
| Adult | Lussault-sur-Loire ( | nt | nt | nt | nt | nt | nt | + | nt | nt | |
| Adult | Fondettes ( | nt | nt | nt | nt | nt | nt | nt | + | nt | |
| Adult | Tours ( | nt | nt | nt | nt | nt | nt | nt | + | nt |
nt: not tested; only the positive samples are included in the table.
Co-occurrence of virus species in vespid colonies in China.
| Colony ID | Vespid Species | Developmental Stage | Viruses Present |
|---|---|---|---|
| B2 |
| Adult | DWV+BQCV |
| GD5 |
| Adult | IAPV+DWV |
| B4 |
| Pupa | AmFV+BQCV |
| B5 |
| Pupa | AmFV+IAPV |
| GD1 |
| Pupa | DWV+IAPV |
| B1 |
| Pupa | AmFV+IAPV+DWV |
| GD4 | Pupa | IAPV+DWV+BQCV | |
| B2 |
| Pupa | AmFV+IAPV+DWV+CSBV+BQCV+LSV |
Figure 2Maximum Likelihood phylogenetic trees of new AmFV, BQCV, IAPV, and SBV sequences found in China and France. (a) Apis mellifera filamentous virus phylogeny (524 sites, TN93+G model), (b) Black queen cell virus phylogeny (634 sites, GTR+I+G model), (c) Israeli acute bee paralysis virus (439 sites, TN93+I model) and (d) Sacbrood virus phylogeny (277 sites, GTR+G model). Names of countries and species for each known and new sequence were added in each leaf of the tree, joined with the corresponding GenBank accession number. The scale bars represent the substitution rate per site. The values correspond to the values of nodes supported by aLRT statistics.
Figure 3Concatenate Maximum Likelihood phylogenetic tree of Deformed wing virus sequences found in China and France. This phylogeny was built by concatenated four sequences’ alignment regions, due to multiple PCR detection primers used in protocols (Region 1: 2170–2438; region 2: 2843–3182; Region 3: 6111–6310; Region 4: 8377–8753) (1004 sites, GTR+G model). To tree topology of the concatenation was verified using each independent region in their respective phylogenetics tree (Supplementary Figure S2a–d). Country and species names for each known and new sequence were added in each leaf of the tree, joined with the corresponding GenBank accession number. The scale bars represent the substitution rate per site. The values correspond to the values of nodes supporting by aLRT statistics.