# Origin and invasion of the emerging infectious pathogen Sphaerothecum destruens

Salma Sana, Emilie A Hardouin, Rodolphe E Gozlan, Didem Ercan, Ali Serhan Tarkan, Tiantian Zhang, Demetra Andreou

► To cite this version:

Salma Sana, Emilie A Hardouin, Rodolphe E Gozlan, Didem Ercan, Ali Serhan Tarkan, et al.. Origin and invasion of the emerging infectious pathogen Sphaerothecum destruens. Emerging microbes & infections, 2017, 6, pp.e76. 10.1038/emi.2017.64. hal-01589799

HAL Id: hal-01589799

https://hal.sorbonne-universite.fr/hal-01589799

Submitted on 19 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_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
ORIGINAL ARTICLE

# Origin and invasion of the emerging infectious pathogen *Sphaerothecum destruens*

Salma Sana<fnref n="1" />, Emilie A Hardouin<fnref n="1" />, Rodolphe E Gozlan², Didem Ercan³, Ali Serhan Tarkan³, Tiantian Zhang<fnref n="1" /> and Demetra Andreou<fnref n="1" />

Non-native species are often linked to the introduction of novel pathogens with detrimental effects on native biodiversity. Since *Sphaerothecum destruens* was first discovered as a fish pathogen in the United Kingdom, it has been identified as a potential threat to European fish biodiversity. Despite this parasite's emergence and associated disease risk, there is still a poor understanding of its origin in Europe. Here, we provide the first evidence to support the hypothesis that *S. destruens* was accidentally introduced to Europe from China along with its reservoir host *Pseudorasbora parva* via the aquaculture trade. This is the first study to confirm the presence of *S. destruens* in China, and it has expanded the confirmed range of *S. destruens* to additional locations in Europe. The demographic analysis of *S. destruens* and its host *P. parva* in their native and invasive range further supported the close association of both species. This research has direct significance and management implications for *S. destruens* in Europe as a non-native parasite.

*Emerging Microbes & Infections* (2017) **6**, e76; doi:10.1038/emi.2017.64; published online 23 August 2017

**Keywords:** aquaculture; biological invasion; fungal pathogens; invasive

## INTRODUCTION

Fungal parasites have emerged as serious threats to biodiversity in the past 20 years, resulting in the worldwide decline of several taxa.<fnref n="1" /> Facilitated by the international trade in live animals, diseases such as chytridiomycosis and white-nose syndrome have resulted in high levels of mortalities in amphibian and North American bat populations,²⁻⁵ respectively. Emergent fungi are frequently generalists, and they can survive in a wide range of temperatures. They are often transmitted both directly and indirectly, enabling their transport and introduction to native communities across the globe.<fnref n="1" />,⁶,⁷ Wild population declines are difficult to detect, especially for aquatic species.⁸

Long-term monitoring of wild fish populations in combination with histological checks has provided the first empirical evidence that links the fungus-like parasite *Sphaerothecum destruens* to declines in endemic fish species in Turkey.⁹ This parasite was first identified in Europe in 2005 after it caused high fish mortalities in semi-natural and lab experiments.¹⁰ Its reservoir host, the invasive fish *Pseudorasbora parva*, is the suspected source of the parasite in Europe. *P. parva* was accidentally introduced to Europe from China and has invaded 32 countries in <40 years.¹¹ As it was first reported in Europe, *S. destruens* has been highlighted as a true generalist and an important pathogen of fishes, with over 14 species as its hosts (including important aquaculture species, such as salmon, carp and sea bass).⁹,¹²,¹³ Currently, there are two known geographically distinct strains of this parasite—one in North America and one in Eurasia.¹⁴

Disease outbreaks by *S. destruens* have been reported in both wild and cultured salmonid fishes on the west coast of the USA.¹⁵,¹⁶ It has caused high mortality in cultured Chinook salmon (*Oncorhynchus tshawytscha*)¹⁵,¹⁷ and in subadult Atlantic salmon (*Salmo salar*)¹⁶. Further work has confirmed the susceptibility of other salmonids (Coho salmon (*O. kisutch*), rainbow trout (*O. mykiss*) and brown trout (*S. trutta*)).¹⁷ The European strain was first detected in the cyprinid fishes sunbleak (*Leucaspius delineatus*) and fathead minnow (*Pimephales promelas*) with high mortalities.¹⁰ Experimental infections with the European strain have confirmed that it can cause the same pathology and morbidity in *S. salar* as the North American strain.¹³ Furthermore, the European strain can cause mortality in more cyprinid fishes, such as carp (*Cyprinus carpio*), roach (*Rutilus rutilus*), bream (*Abramis brama*), rudd (*Squalius fellowesii*)¹² and topmouth gudgeon (*Pseudorasbora parva*),¹⁸ as well as *Oxynoemacheilus* sp. (Family Nemacheilidae) and *Lepomis gibbosus* (Family Centrachidae).⁹ Despite the parasite's emergence and associated disease risk, there is still a poor understanding of its origin in Europe.

*S. destruens* is an obligate intracellular parasite belonging to the Rhinosporideaceae family and the Mesomycetozoea class.¹⁹ Within its fish host, *S. destruens* spores can infect multiple organs (liver, kidney, gonads, intestine, and gills⁹). The spores divide asexually within the cells, eventually leading to cell death, releasing the spores, and these spores then infect new cells. Spores are released in urine and seminal fluids and can cause infection to native fish through both direct transmission (contact with spores) and indirect transmission. Indirect

<footnote n="1">¹Bournemouth University, Faculty of Science and Technology, Fern Barrow, Talbot Campus, Poole, Dorset BH12 5BB, UK; ²UMR BOREA IRD-MNHN-Université Pierre et Marie Curie, Muséum National d'Histoire Naturelle, 47 Rue Cuvier, Paris, Cedex 5 75231, France and ³Faculty of Fisheries, Muğla Sitki Koçman University, Kötekli, Muğla 48000, Turkey</footnote>

<footnote>Correspondence: S Sana</footnote>

<footnote>E-mail: ssana@bournemouth.ac.uk</footnote>

<footnote>Received 8 February 2017; revised 7 June 2017; accepted 23 June 2017</footnote>
transmission is facilitated by zoospores, which are produced once the spores are incubated in freshwater.¹⁹ Zoospores and spores have a wide temperature tolerance (5–30 °C), potentially facilitating the spread of this parasite to new environments.⁶ Within tissues, *S. destruens* can cause disseminated and granulomatous histopathology.²⁰

In this paper, we provide the first evidence that *S. destruens* has invaded Europe from China along with its reservoir host, the highly invasive small cyprinid species *P. parva*. This result was achieved through extensive sampling of *P. parva* across its native and invasive range, coupled with pathogen detection in these populations. The close links of the reservoir host with aquaculture and the high susceptibility of aquaculture and native species to *S. destruens* raise serious threats to European freshwater biodiversity with important policy implications.

## MATERIALS AND METHODS

### *S. destruens* detection

A total of 420 *P. parva* from 21 populations was sampled across *P. parva*'s native and invasive range (Table 1; Figure 1). *P. parva* kidney DNA was extracted to investigate the presence of *S. destruens* in the fish tissue as the parasite is found at a high prevalence in this organ.²¹ All *P. parva* samples were collected during 2009–2010 and

<table><caption>**Table 1** Sampled populations of *Pseudorasbora parva* and the distribution of *Sphaerothecum destruens* in *P. parva* using molecular detection across *P. parva*'s native and non-native range</caption><thead><tr><th rowspan="2">Population</th><th rowspan="2">Code</th><th colspan="2">Coordinates</th><th rowspan="2">Prevalence of<br/>S. destruens</th><th colspan="3">Successfully amplified<br/>markers for every<br/>positive sample</th></tr><tr><th>X</th><th>Y</th><th>18S<br/>rRNA</th><th>Cyt-b</th><th>ITS 1</th></tr></thead><tbody><tr><td>China 1</td><td>S1</td><td>115.56</td><td>37.55</td><td>10% (2/20)</td><td>✓</td><td>✓</td><td></td></tr><tr><td>China 2</td><td>S2</td><td>117.12</td><td>34.81</td><td>5% (1/20)</td><td>✓</td><td></td><td></td></tr><tr><td>China 3</td><td>S3</td><td>118.59</td><td>33.19</td><td>5% (1/20)</td><td>✓</td><td>✓</td><td>✓</td></tr><tr><td>China 7</td><td>S7</td><td>110.32</td><td>25.27</td><td>0% (0/20)</td><td></td><td></td><td></td></tr><tr><td>China 9</td><td>S9</td><td>113.11</td><td>29.15</td><td>5% (1/20)</td><td>✓</td><td></td><td></td></tr><tr><td>China 11</td><td>S11</td><td>110.99</td><td>34.62</td><td>5% (1/20)</td><td>✓</td><td>✓</td><td>✓</td></tr><tr><td>China 12</td><td>S12</td><td>117</td><td>38.7</td><td>10% (2/20)</td><td>✓</td><td>✓</td><td></td></tr><tr><td>China 13</td><td>S13</td><td>122.52</td><td>40.1</td><td>5% (1/20)</td><td>✓</td><td>✓</td><td>✓</td></tr><tr><td>China 14</td><td>S14</td><td>124.99</td><td>45.03</td><td>10% (2/20)</td><td>✓</td><td></td><td></td></tr><tr><td>China 16</td><td>S16</td><td>118.27</td><td>40.9</td><td>5% (1/20)</td><td>✓</td><td></td><td></td></tr><tr><td>Austria</td><td>A</td><td>14.72</td><td>48.19</td><td>0% (0/20)</td><td></td><td></td><td></td></tr><tr><td>Bulgaria</td><td>BG</td><td>43</td><td>26</td><td>0% (0/20)</td><td></td><td></td><td></td></tr><tr><td>France</td><td>F</td><td>-1.73</td><td>47.1</td><td>0% (0/20)</td><td></td><td></td><td></td></tr><tr><td>Iran</td><td>IR</td><td>54.78</td><td>37.05</td><td>0% (0/20)</td><td></td><td></td><td></td></tr><tr><td>Italy</td><td>IT</td><td>10</td><td>44</td><td>0% (0/20)</td><td></td><td></td><td></td></tr><tr><td>Japan</td><td>JP</td><td>139.43</td><td>35.67</td><td>0% (0/20)</td><td></td><td></td><td></td></tr><tr><td>Morocco</td><td>M</td><td>32.11</td><td>2.89</td><td>0% (0/20)</td><td></td><td></td><td></td></tr><tr><td>Spain</td><td>SE</td><td>0.86</td><td>40.7</td><td>5% (1/20)</td><td>✓</td><td></td><td></td></tr><tr><td>Turkey</td><td>T</td><td>30.04</td><td>40.91</td><td>0% (0/20)</td><td></td><td></td><td></td></tr><tr><td>United Kingdom</td><td>UK</td><td>1</td><td>51</td><td>5% (1/20)</td><td>✓</td><td>✓</td><td></td></tr><tr><td>Hungary</td><td>H</td><td>18</td><td>46</td><td>0% (0/20)</td><td></td><td></td><td></td></tr></tbody></table>

were fixed in 92% ethanol upon collection. The DNA was extracted using Qiagen DNeasy Blood and Tissue kit following the manufacturer's protocol. *S. destruens* 18S rRNA was amplified using the method described in Mendonca and Arkush,²² which is currently the most sensitive DNA-based detection method for this parasite. All samples identified as positive for *S. destruens* DNA using 18S rRNA were also sequenced for the ITS 1 and mitochondrial DNA Cyt-b regions (total length of 700 bp, spanning the Cyt-b gene (629 bp), the Cyt-b-COI intergenic region (60 bp) and the COI gene (11 bp)) in order to investigate the putative origin of *S. destruens*. Disease prevalence was calculated as: (number of *S. destruens* positive fish/ total number of fish tested) × 100.

The nested PCR for the ITS 1 region (623 bp) used the primers Sdes2F (5'-CTT CGG ATT GGC CCT GTA C-3')¹⁴ and NC 2 (5'-TTA GTT TCT TTT CCT CCG CT-3')²³ in the first step-PCR, and Sdes2F and SD-ITS R1 (5'- CGA TGC ACG AGC CAA GAG-3')⁹ in the second step-PCR. The PCR conditions were 1× Promega (Madison, WI, USA) Flexi buffer, 1.5 mM MgCl₂, 0.2 mM dNTPs, 0.3 µM primers and 0.5 U Taq polymerase (Promega). PCR cycling conditions included an initial denaturation at 95 °C for 3 min, followed by 35 cycles of 30 s at 95 °C, 45 s at X °C and 90 s at 62 °C, and final extension at 62 °C for 7 min (X= annealing temperatures were 60 °C and 59 °C for the first and the second PCR, respectively). The amplified product was gel extracted before sequencing. Cyt-b (700 bp) sequences were obtained for *S. destruens* using nested PCR. The primers used were Nt-Cytb-F1 (5-ATG AGT TTA TGG GGA GCG) coupled with Nt-Cytb-R1 (5-GCT CCA GCC AAC ACA GGT AAG GAT AAT AAC) in the first step-PCR. The second step-PCR used the primer Nt-Cytb-F2 (5-GGA GGG TTT AGT GTG GAT AAT GC) coupled with Nt-Cytb-R1 (5-TCA TCG TCA AAT CCA ACT CAC C). The PCR conditions were 1× Promega Flexi buffer, 1.5 mM MgCl₂, 0.2 mM dNTPs, 0.3 µM forward and reverse primer and 0.5 U Taq polymerase (Promega). Cycling conditions included an initial denaturation at 95 °C for 2 min, followed by 35 cycles with 95 °C for 40 s, X °C for 40 s, and 72 °C for 60 s, and a final extension step at 72 °C for 5 min (X= annealing temperatures were 56 °C and 58 °C for the first and the second PCR step, respectively).

Sequences were cleaned and aligned using Clustal W in BioEdit Version 5.0.9.²⁴ A phylogenetic network was generated for the 18S rRNA marker using DnaSP version 5.10,²⁵ and Network and Network Publisher²⁶ (Available at: http://www.fluxus-engineering.com) using the sequences generated in the present study and all the published sequences available for this marker (FN996945.1, AY267344.1, AY267345.1, and AY267346.1¹³,¹⁹). A phylogenetic tree was drawn for the ITS 1 marker using MrBayes²⁷ with the sequences generated in the present study and the sequences available in GenBank (FJ440707.1, FJ440708.1, FJ440709.1, FJ440702.1, FJ440703.1, FJ440704.1 and KT361608.1⁹,¹⁴). The best model fitting our data was Hasegawa–Kishino–Yano (HKY), and it was investigated with jModel test v2.1. 4.²⁸ The number of haplotypes and haplotype diversity were calculated using DnaSP.²⁵ All the sequences obtained through this work have been deposited in GenBank (MF138119, MF062546-MF062560 and MF101749-MF101755).

### *Pseudorasbora parva* Cyt-b sequencing

PCR amplification of the Cyt-b gene was carried out on *P. parva* individuals that had tested positive for *S. destruens* using the primers L15267 (5-AAT GAC TTG AAG AAC CAC CGT-3') and H15891Ph (5-GTT TGA TCC CGT TTC GTG TA-3')²⁹ with an amplicon size of ~600 bp. The reaction conditions included a 50-µL reaction volume with 100 ng of template DNA, 1 × Promega Flexi buffer, 2 mM MgCl₂,
<figure><img src="image_2.png" /><figcaption>**Figure 1** *Pseudorasbora parva* sites screened in this study for *Sphaerothecum destruens* presence in the host's native range (China) and invasive range (Europe, North Africa and Eurasia). The red sites were positive for *S. destruens* and the black sites were negative using the 18S rRNA marker. Europe- Austria, A; Bulgaria, BG; France, F; Hungary; H, Italy, IT; Spain, SP; Turkey, T; United Kingdom, UK. North Africa- M, Morocco; Eurasia- I, Iran; Asia- China (S1, S2, S3, S7, S9, S11, S12, S13, S14 and S16); Japan, JP.</figcaption></figure>

0.2 mM dNTPs, 0.3 µM forward and reverse primer and 0.5 U Taq polymerase. The cycling conditions included an initial denaturation at 95 °C for 15 min followed by 35 cycles of 30 s at 95 °C, 90 s at 60 °C and 60 s at 72 °C, with a final elongation step at 72 °C for 15 min.

The Cyt-b sequences for *P. parva* that had tested positive for *S. destruens* were aligned with all available *P. parva* Cyt-b sequences across populations in its native range and invasive range³⁰ (JF489575–JF489887, and KR074432–KR074994). These populations included the populations tested for *S. destruens* in this study. The sequences were aligned by Clustal W in BioEdit.²⁴ Haplotype diversity was calculated in DnaSP version 5.10.²⁵ A phylogenetic tree was constructed to identify the *P. parva* haplotypes associated with the presence of *S. destruens*. The phylogenetic analysis was performed using MrBayes,²⁷ and posterior probabilities were obtained after 2 500 000 generations with a burn-in of 25%. The tree was calculated using a Hasegawa–Kishino–Yano model with a Gamma distribution (HKY +G) model³¹ determined with jModel test v 2.1.4.²⁸ The Cyt-b gene sequences from *Ictiobus bubalus* (JF799443.1), *Hypentelium nigricans* (JF799441.1) and *Danio rerio* (JN234356.1) were used as outgroups.

## Population demographic analysis

The association between *S. destruens* and *P. parva* in China was further investigated by inferring their demographic history. A strong association between the two species would result in them sharing a similar demographic history. This analysis was achieved by determining the mismatch distribution using the 18S rRNA gene for *S. destruens* and the Cyt-b gene for *P. parva*. This analysis plots the distribution of nucleotide differences between each pair of sequences and compares it to the expected values for a model of population expansion. A unimodal distribution is indicative of a population expansion in the recent past whereas a bimodal/multimodal distribution indicates that a population is at demographic equilibrium.³² Demographic changes were analyzed by calculating Harpending's raggedness index (*Hri*), which quantifies the smoothness of the observed mismatch distribution³³ and the sum of squared deviations between the observed and expected mismatch for the nucleotide differences³⁴ in Arlequin version 3.5.³⁵

Three statistical methods were used to test for the population expansion of *P. parva*: Fu's *F*<sub>s</sub> test,³⁶ Tajima's *D* test³⁷ and Ramos-

Onsins & Rozas' *R*<sub>2</sub> test.³⁸ The *R*<sub>2</sub> test was performed by DnaSP,²⁵ and Fu's *F*<sub>s</sub> and Tajima's *D* tests were carried out in Arlequin version 3.5.³⁵ Only Tajima's *D* test and the *R*<sub>2</sub> test were performed for *S. destruens* as Fu's *F*<sub>s</sub> test is not suitable for small sample sizes.³⁸ For Tajima's *D* and *F*<sub>s</sub>, *P*-values were calculated based on a coalescent simulation algorithm, and for the *R*<sub>2</sub> test, the *P*-values were based on parametric boot strapping with coalescence simulations.

## RESULTS

### Prevalence

Of the 10 Chinese populations tested in this study, 9 were found to be positive for the presence of *S. destruens* (Table 1 and Figure 1). The prevalence of this parasite in the Chinese populations ranged from 0 to 10%. The overall prevalence of *S. destruens* across all Chinese populations was 6% (12/200). *S. destruens* was also found in two European populations: Spain and the United Kingdom (Figure 1 and Table 1), with a prevalence of 5% (1/20) in both populations. Overall, the prevalence in Europe was 1.4% (2/140). *S. destruens* was not detected in Morocco, Iran or Japan.

### Genetic diversity of the parasite *S. destruens*

A sequence of 397 bp of 18S rRNA of *S. destruens* was obtained for 14 *P. parva* individuals from 11 different populations (Table 1), and it was aligned with the sequences published in GenBank (FN996945.1AY267344.1, AY267345.1 and AY267346.1¹³,¹⁹). The haplotype diversity of *S. destruens* across all populations (native and invasive) was 0.22 with four identified haplotypes. Only three individuals, one from the USA (host: *S. salar*) and two from China (host: *P. parva*), were found to display different haplotypes (Figure 2).

We managed to sequence the ITS 1 region from only one of the 11 *S. destruens*-positive populations (Chinese population S3). However, when combined with the published ITS 1 sequences, the overall haplotypic diversity was high (Hd = 0.97, Figure 3). Three clades were identified in our study. Individuals originating from the UK and China clustered together, the Turkish samples grouped by themselves, and the North American samples constituted the third clade (Figure 3). The Turkish samples were more closely related to the UK and China (S3) samples than the USA strains (Figure 3), indicating that the two European strains are closely related to the Chinese strain. In addition,
<figure><img src="image_3.png" /><figcaption>**Figure 2** Minimum spanning network based on the 18S rRNA segment (397 bp) of *Sphaerothecum destruens*. The size of the different circles represents the frequency of each respective haplotype. The numbers on the branches indicate the number of mutations between the nodes. Black circles indicate branch splits. The color code indicates *S. destruens* individuals from different localities; Red: China (*n*=11), Yellow: Spain (*n*=1), Light blue: Turkey (*n*=1), Brown: UK (*n*=2) and Green: USA (*n*=3). Haplotypes H_1, H_2, H_3 and H_4 represent *S. destruens* 18S rRNA haplotypes. The Cyt-b *P. parva* haplotype for individuals positive for *S. destruens* (with *S. destruens* haplotypes H_3 and H_4) was Haplotype 6.</figcaption></figure>

<figure><img src="image_4.png" /><figcaption>**Figure 3** Unrooted phylogenetic tree resulting from Bayesian inference method in MrBayes²⁷ based on Hasegawa-Kishino-Yano model³¹ analysis of the ribosomal ITS 1 sequences of *Sphaerothecum destruens*. Isolate origins and GenBank accession numbers are RA 1-3 (FJ440707.1), RA 3-1 (FJ440708.1), RA 3-2 (FJ440709.1), RA 3-3 (FJ440710.1), RA 4-1 (FJ440702.1), RA 4-3 (FJ440703.1), RA 4-4 (FJ440704.1), RA-Turkey (KT361608.1) and RA-China (MF138119).</figcaption></figure>

a 700-bp fragment of *S. destruens* Cyt-b was successfully amplified in six *P. parva* populations (China (S1, S3, S11, S12 and S13) and the UK). No Cyt-b DNA sequence variation was found between the Chinese and the UK *S. destruens* samples (Supplementary Figure S1), which supported the IT S1 results. In combination, the close clustering between *S. destruens* isolated from Europe and China supports the hypothesis that the parasite was introduced to Europe from China.

## Genetic diversity of the host *P. parva*

A total of 91 haplotypes from 949 *P. parva* individuals were identified in the dataset of Cyt-b sequences of *P. parva* populations (a total of 62 *P. parva* populations) across its invasive and native range. The two main haplogroups identified in China were Haplogroups A and B

(Figure 4). Six Cyt-b haplotypes in *P. parva* were found to be associated with the presence of *S. destruens* (*Hap_1*, *Hap_4*, *Hap_6*, *Hap_7*, *Hap_12* and *Hap_55*; Figure 4). The highest number of *P. parva* individuals (*n* = 7) that were positive for *S. destruens* had the Cyt-b haplotype *Hap_6*. The remaining haplotypes each had one *P. parva* individual positive for *S. destruens*.

## Demographic analysis of *P. parva* and *S. destruens*

The demographic analysis of the host, *P. parva*, and the parasite, *S. destruens*, suggested a potential recent population expansion in both species, supporting the close relationship between the two species. The sum of squared differences and Harpending's raggedness index (*Hri*) were not significant for the pathogen and its host, indicating that the data are a
<figure><img src="image_5.png" /><figcaption>**Figure 4** Molecular phylogenetic analysis of Cyt-b haplotypes of *Pseudorasbora parva* populations across its invasive and native range. The tree was inferred from the Bayesian inference method based on the Hasegawa-Kishino-Yano model³¹ with Gamma distribution in MrBayes.²⁷ The colored circles indicate the countries that each haplotype has been found in and the colored stars indicate *S. destruens* positive haplotypes in that country.</figcaption></figure>

relatively good fit with population expansion (Figure 5). Both species also
had significant negative values for the Tajima's *D* test, further supporting
population expansion. *P. parva's* population expansion was further

supported by the R<sub>2</sub> test. The R<sub>2</sub> test, however, was not significant for
the *S. destruens* population, which was in contrast to the negative value of
Tajima's *D* and the unimodal mismatch distribution.
<figure><img src="image_6.png" /><figcaption>**Figure 5** Frequency distributions of the number of pairwise nucleotide differences (mismatch) between (A) 18S rRNA haplotypes of *Sphaerothecum destruens* populations in China; (B) cytochrome b haplotypes of *Pseudorasbora parva* populations in China; (C) *P. parva* haplogroup A, and (D) *P. parva* haplogroup B. The solid line is the theoretical distribution under the hypothesis of population expansion. Sum of squared differences, Harpending's raggedness index (*Hri*), Fu's *F*<sub>s</sub>, Tajima's *D* and Ramos-Onsins and Rozas' *R*<sub>2</sub> statistics are listed next to each dataset. The *P*-values for each statistical test can be found in parentheses. Significance was set at a *P*-value of 0.05, except for *F*<sub>s</sub>, for which it was set at 0.01.</figcaption></figure>

The observed mismatch distribution for *P. parva* when all the populations are considered (Figure 5B) is bimodal. The bimodality of the mismatch distribution for *P. parva* could be due to the presence of different haplogroups.<sup>39</sup> In China, there are two main established haplogroups, A and B (Figure 4), despite the statistical test supporting that the population as a whole has undergone a recent population expansion. When the data were split by haplogroups, the mismatch distribution was unimodal (Figures 5C and 5D).

DISCUSSION

Our results support that *S. destruens* was introduced to Europe from China via its reservoir host *P. parva*, designating *S. destruens* as a non-native parasite to Europe. As a non-native parasite to Europe, *S. destruens* will need to be risk assessed<sup>40</sup> with respective fish movement limitations on fishes that test positive for the parasite. This is especially important in regions where *S. destruens* has been associated with population declines of native freshwater fishes.<sup>9,10</sup> The healthy reservoir of *S. destruens*, *P. parva*, has been accidentally
introduced to Europe through the aquaculture trade of Asian carp (Hypophthalmichthys molitri, Ctenopharyngodonidella).¹¹ It is often challenging to determine whether a parasite was introduced, as parasitological surveys are neither regular nor exhaustive.⁴¹ The disseminated form of disease caused by *S. destruens* is difficult to detect through traditional histo-parasitological surveys, which focus on gross pathological changes in the organs. Thus, we used molecular tools to test the hypothesis that *S. destruens* was introduced to Europe using an extensive survey of *P. parva* across its native and invasive range.

No geographical isolation was identified between the Chinese and the European *S. destruens* populations using the phylogenetically informative marker ITS 1. The results indicate that there are two geographical clades for *S. destruens*, one encompassing the samples from North America and one including the samples from China, the UK and Turkey. However, within the European clade, the UK and China sequences were more closely related compared with the sequence from Turkey. This similarity can be explained by the invasion history of the parasite's reservoir host, *P. parva*. The UK populations of *P. parva* genetically group with *P. parva* populations found north of the River Yangtze (Haplogroup A; Figure 4). However, the *P. parva* found in Turkey is genetically more similar to populations found to the south of the River Yangtze⁴² (Haplogroup B; Figure 4), which could explain why the Turkish *S. destruens* isolate (Figure 3) grouped on its own within the overall European and China clade. Thus, the genetic diversity of the host species, *P. parva*, is likely to reflect the observed diversity of *S. destruens*, further supporting the hypothesis that the *S. destruens* populations found in Europe have been introduced from China.

Using two nuclear and one mitochondrial marker, we demonstrated that *S. destruens* is widely distributed in China (present in 90% of the sampled Chinese *P. parva* populations). A higher presence of the parasite in its native range would be expected because the process of species introduction can lead to the loss of associated parasites through stochastic effects.⁴³,⁴⁴ The wide distribution of the *P. parva* populations that are positive for *S. destruens* across China suggests that the two species could share a long co-evolutionary history. Similarly, the emergence and association of fungal parasites with reservoir hosts over a long evolutionary history has been recently demonstrated for the chytrid fungus *Batrachochytrium salamandrivorans*.⁴⁵ *B. salamandrivorans* has likely originated and coexisted with its reservoir hosts for millions of years in Asia before being introduced across the world with the trade of its reservoir hosts.⁴⁵ Similar to *B. salamandrivorans*, *S. destruens* appears to have been introduced to Europe via the accidental introduction of its reservoir host.

The demographic analysis of *P. parva* and *S. destruens* partially suggests that both species have undergone a recent population expansion. The partial congruence between the demographic history of the two species is surprising, especially in light of the true generalist nature of *S. destruens* and its ability to use a number of different hosts. Recent work, however, has indicated that, following the establishment of a generalist parasite in a community, its population dynamics are driven via intra-host transmission rather than by inter-host transmission.⁴⁶,⁴⁷ This could explain the observed similarity in the population demographic history of the two species.

This study represents the first screening of the native and invasive *P. parva* populations for the presence of *S. destruens*. It is important to note that the prevalence values recorded in positive populations are very likely to be underestimates of the true prevalence of this parasite, as only the kidney was sampled. *S. destruens* infects multiple organs and does so unequally.²¹ This lack of infection localization makes this

organism harder to detect. Populations that have been detected as negative for *S. destruens* in this study must be treated with caution as it cannot be excluded that the parasite might still be present in these and other populations in the country. For example, the Turkish population sampled in this study was found to be negative for *S. destruens*, but the parasite has been detected in another *P. parva* population in Turkey.⁹ The screening recommendations for *S. destruens* include sampling multiple populations, and where possible, a minimum sample size of 30 fish should be tested and multiple tissues (kidney, liver, testis/ova) should be tested to increase the probability of detection.⁴⁰

Currently, there are 14 known species that are susceptible to *S. destruens*, including valuable aquaculture species (salmon, carp and sea bass)⁹,¹²,¹³ and endemic fishes to Europe that are of high conservation value (*L. delineatus*, *Oxynoemacheilus* sp., *Petroleuciscus smyrnaeus*, *S. fellowesii* and *Dicentrarchus labrax*).⁹,²¹ The close association of the reservoir host, *P. parva*, with aquaculture facilities (due to its accidental introduction along with carp from China)¹¹ and the ability of the parasite to establish in local freshwater communities within a year of its introduction⁴⁶ increases the risk of disease to native fishes. In the last 20 years, aquaculture production has increased exponentially to support economic growth, with its expansion being highly reliant on non-native species.⁴⁸ The introduction of non-native species can be detrimental both to ecosystem services and naive communities.⁴⁹ The potential threats associated with aquaculture production and the resultant fish movements highlight the importance of risk assessments to identify emergent parasites. Horizon scanning for potential emergent diseases will be critical in informing strict biosecurity controls in order to prevent disease introduction.

## ACKNOWLEDGEMENTS

We would like to thank all the colleagues who have helped us collect samples across the native and invasive populations.

1. Fisher MC, Henk DA, Briggs CJ *et al*. Emerging fungal threats to animal, plant and ecosystem health. *Nature* 2012; **484**: 186–194.

2. Frick WF, Pollock JF, Hicks AC *et al*. An emerging disease causes regional population collapse of a common North American bat species. *Science* 2010; **329**: 679–682.

3. Gilbert M, Bickford D, Clark L *et al*. Amphibian pathogens in Southeast Asian frog trade. *EcoHealth* 2012; **9**: 386–398.

4. Kolby JE, Ramirez SD, Berger L *et al*. Terrestrial dispersal and potential environmental transmission of the amphibian chytrid fungus (*Batrachochytrium dendrobatidis*). *PLoS ONE* 2015; **10**: e0125386.

5. Schloegel LM, Toledo LF *et al*. Novel, panzootic and hybrid genotypes of amphibian chytridiomycosis associated with the bullfrog trade. *Mol Ecol* 2012; **21**: 5162–5177.

6. Andreou D, Gozlan R, Paley R. Temperature influence on production and longevity of *Sphaerothecum destruens* zoospores. *J Parasitol* 2009; **95**: 1539–1541.

7. Mitchell KM, Churcher TS, Garner TW *et al*. Persistence of the emerging pathogen *Batrachochytrium dendrobatidis* outside the amphibian host greatly increases the probability of host extinction. *Proc R Soc Lond B Biol Sci* 2008; **275**: 329–334.

8. Gozlan R. Pathogens: Monitoring fungal infections in fish. *Nature* 2012; **485**: 446–446.

9. Ercan D, Andreou D, Sana S *et al*. Evidence of threat to European economy and biodiversity following the introduction of an alien pathogen on the fungal-animal boundary. *Emerg Microbes Infect* 2015; **4**: e52.

10. Gozlan RE, St-Hilaire S, Feist SW *et al*. Biodiversity: disease threat to European fish. *Nature* 2005; **435**: 1046–1046.

11. Gozlan R, Andreou D, Asaeda T *et al*. Pan-continental invasion of *Pseudorasbora parva*: towards a better understanding of freshwater fish invasions. *Fish Fish* 2010; **11**: 315–340.

12. Andreou D, Arkush KD, Guegan JF *et al*. Introduced pathogens and native freshwater biodiversity: a case study of *Sphaerothecum destruens*. *PLoS ONE* 2012; **7**: e36998.

13. Paley RK, Andreou D, Bateman K *et al*. Isolation and culture of *Sphaerothecum destruens* from Sunbleak (*Leucaspis delineatus*) in the UK and pathogenicity experiments in Atlantic salmon (*Salmo salar*). *Parasitology* 2012; **139**: 904.

14. Gozlan RE, Whipps CM, Andreou D *et al*. Identification of a rosette-like agent as *Sphaerothecum destruens*, a multi-host fish pathogen. *Int J Parasitol* 2009; **39**: 1055–1058.
15 Harrell LW, Elston RA, Scott TM *et al*. A significant new systemic disease of net-pen reared chinook salmon (*Oncorhynchus tshawytscha*) brood stock. *Aquaculture* 1986; **55**: 249-262.

16 Hedrick R, Friedman C, Modin J. Systemic infection in Atlantic salmon Salmo salar with a Dermocystidium-like species. *Dis Aquat Organ* 1989; **7**: 171-177.

17 Arkush K, Frasca Jr S, Hedrick R. Pathology associated with the rosette agent, a systemic protist infecting salmonid fishes. *J Aquat Anim Health* 1998; **10**: 1-11.

18 Spikmans F, van Tongeren T, van Alen TA *et al*. High prevalence of the parasite *Sphaerothecum destruens* in the invasive topmouth gudgeon *Pseudorasbora parva* in the Netherlands, a potential threat to native freshwater fish. *Aquat Invasions* 2013; **8**: 355-360.

19 Arkush KD, Mendoza L, Adkison MA *et al*. Observations on the life stages of *Sphaerothecum destruens* ng, n. sp., a mesomycetozoean fish pathogen formally referred to as the rosette agent. *J Eukaryot Microbiol* 2003; **50**: 430-438.

20 Andreou D, Gozlan RE, Stone D *et al*. *Sphaerothecum destruens* pathology in cyprinids. *Dis Aquat Organ* 2011; **95**: 145-152.

21 Andreou D, Hussey M, Griffiths SW *et al*. Influence of host reproductive state on *Sphaerothecum destruens* prevalence and infection level. *Parasitology* 2011; **138**: 26-34.

22 Mendonca HL, Arkush KD. Development of PCR-based methods for detection of *Sphaerothecum destruens* in fish tissues. *Dis Aquat Organ* 2004; **61**: 187-197.

23 Gasser RB, Chilton NB, Hoste H *et al*. Rapid sequencing of rDNA from single worms and eggs of parasitic helminths. *Nucleic Acids Res* 1993; **21**: 2525.

24 Hall TA. BioEdit: a user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT. *Nucleic Acids Symp Ser* 1999; **41**: 95-98.

25 Librado P, Rozas J. DnaSP v5: a software for comprehensive analysis of DNA polymorphism data. *Bioinformatics* 2009; **25**: 1451-1452.

26 Bandelt H-J, Forster P, Röhl A. Median-joining networks for inferring intraspecific phylogenies. *Mol Biol Evol* 1999; **16**: 37-48.

27 Ronquist F, Teslenko M, van der Mark P *et al*. MrBayes 3.2: efficient Bayesian phylogenetic inference and model choice across a large model space. *Syst Biol* 2012; **61**: 539-542.

28 Darriba D, Taboada GL, Doallo R *et al*. *jModelTest 2*: more models, new heuristics and parallel computing. *Nat Methods* 2012; **9**: 772-772.

29 Briolay J, Galtier N, Brito RM *et al*. Molecular phylogeny of cyprinidae inferred fromcytochrome bDNA Sequences. *Mol Phylogenet Evol* 1998; **9**: 100-108.

30 Simon A, Britton R, Gozlan R *et al*. Invasive cyprinid fish in Europe originate from the single introduction of an admixed source population followed by a complex pattern of spread. *PLoS ONE* 2011; **6**: e18560.

31 Hasegawa M, Kishino H, Yano T-A. Dating of the human-ape splitting by a molecular clock of mitochondrial DNA. *J Mol Evol* 1985; **22**: 160-174.

32 Rogers AR, Harpending H. Population growth makes waves in the distribution of pairwise genetic differences. *Mol Biol Evol* 1992; **9**: 552-569.

33 Harpending HC, Sherry ST, Rogers AR *et al*. The genetic structure of ancient human populations. *Curr Anthropol* 1993; **34**: 483-496.

34 Schneider S, Excoffier L. Estimation of past demographic parameters from the distribution of pairwise differences when the mutation rates vary among sites: application to human mitochondrial DNA. *Genetics* 1999; **152**: 1079-1089.

35 Excoffier L, Lischer HE. Arlequin suite ver 3.5: a new series of programs to perform population genetics analyses under Linux and Windows. *Mol Ecol Resour* 2010; **10**: 564-567.

36 Fu Y-X. Statistical tests of neutrality of mutations against population growth, hitchhiking and background selection. *Genetics* 1997; **147**: 915-925.

37 Tajima F. Statistical method for testing the neutral mutation hypothesis by DNA polymorphism. *Genetics* 1989; **123**: 585-595.

38 Ramos-Onsins SE, Rozas J. Statistical properties of new neutrality tests against population growth. *Mol Biol Evol* 2002; **19**: 2092-2100.

39 Maltagliati F, Di Giuseppe G, Barbieri M *et al*. Phylogeography and genetic structure of the edible sea urchin *Paracentrotus lividus* (Echinodermata: Echinoida) inferred from the mitochondrial cytochrome *b* gene. *Biol J Linn Soc* 2010; **100**: 910-923.

40 Andreou D, Gozlan RE. Associated disease risk from the introduced generalist pathogen *Sphaerothecum destruens*: management and policy implications. *Parasitology* 2016; **143**: 1204-1210.

41 Lymbery AJ, Morine M, Kanani HG *et al*. Co-invaders: the effects of alien parasites on native hosts. *Int J Parasitol Parasites* 2014; **3**: 171-177.

42 Simon A. Genetic factors affecting establishment during invasions: the introduction of the topmouth gudgeon (*Pseudorasbora parva*) and the rainbow trout (*Oncorhynchus mykiss*) in Europe. University of Hull 2012.

43 Sheath DJ, Williams CF, Reading AJ *et al*. Parasites of non-native freshwater fishes introduced into England and Wales suggest enemy release and parasite acquisition. *Biol Invasions* 2015; **17**: 2235-2246.

44 Torchin ME, Lafferty KD, Dobson AP *et al*. Introduced species and their missing parasites. *Nature* 2003; **421**: 628-630.

45 Martel A, Blooi M, Adriaensen C *et al*. Recent introduction of a chytrid fungus endangers Western Palearctic salamanders. *Science* 2014; **346**: 630-631.

46 Al-Shorbaji F, Roche B, Gozlan R *et al*. The consequences of reservoir host eradication on disease epidemiology in animal communities. *Emerg Microbes Infect* 2016; **5**: e46.

47 Fenton A, Streicker DG, Petchey OL *et al*. Are all Hosts created equal? Partitioning host species contributions to parasite persistence in multihost communities. *Am Nat* 2015; **186**: 610-622.

48 Peeler EJ, Oidtmann BC, Middlyng PJ *et al*. Non-native aquatic animals introductions have driven disease emergence in Europe. *Biol Invasions* 2011; **13**: 1291-1303.

49 Pelicice FM, Vitule JRS, Lima Junior DP *et al*. A serious new threat to Brazilian freshwater ecosystems: the naturalization of nonnative fish by decree. *Conserv Lett* 2014; **7**: 55-60.

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

This work is licensed under a Creative Commons Attribution 4.0 International License. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in the credit line; if the material is not included under the Creative Commons license, users will need to obtain permission from the license holder to reproduce the material. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/

© The Author(s) 2017

Supplementary Information for this article can be found on the *Emerging Microbes & Infections* website (http://www.nature.com/emi)