Apricot aphid, Myzus mumecola (Matsumura), a new and important pest of apricot in Hungary

Apricot aphid, Myzus mumecola (Matsumura) (Homoptera: Aphididae), was recorded as a new pest of apricot in Hungary in the spring of 2020. Identification was based on morphological and genetical (mitochondrial COI region) characteristics. M. mumecola most likely arrived in Hungary in the last 2–3 years and has quickly become a widespread pest causing significant damage to young apricot trees. Colony development, damage and differences in susceptibility between cultivars are described. The presence of Plum pox virus in M. mumecola samples was detected, and all isolates belonged to the PPV-D subgroup. Illustrations of the most important diagnostic characters of M. mumecola are provided.


Introduction
Apricot (Prunus armeniaca L.) is an economically important crop in Europe, especially in the Mediterranean Region (FAOSTAT 2020). Although the presence of several aphid species was recorded in apricot orchards in Europe (Avinent et al. 1993;Mosco et al. 1997;Holman 2009Jeremovic et al. 2016, their most significant economic impact was not their direct damage caused on the yield, but the transmission of Plum pox potyvirus (PPV) (Avinent et al. 1993;Labonne et al. 1994;García et al. 2014). The most common aphid species landing on the foliage of apricot trees in Spain were Aphis gossypii Glover and A. spiraecola (Avinent et al. 1993). On the other hand, from the West Palearctic Region, the Hyalopterus pruni species complex, Myzus persicae Sulzer and Rhopalosiphum nymphaeae L. have been most commonly reported to form large colonies on apricot trees (Blackman and Eastop, 1984;Mosco et al. 1997;Lozier et al. 2008;Rakauskas et al. 2013). In Italy, Panini et al. (2017) recorded the appearance and rapid spread of an aphid species on apricot, which was identified as Myzus mumecola (Hemiptera: Aphididae), a new aphid species for Europe. Myzus mumecola originates from Eastern Asia where its main host plant is Prunus mume Siebold (Blackman and Eastop, 1984). The species has been recorded in Japan on P. ansu Maxim, P. armeniaca and P. mume (Takahashi, 1965;Miyazaki, 1971), in Taiwan (Mondor et al. 2007;Lai 2020) on Prunus sp., in China on P. armeniaca (Zhang et al. 1985), in East Russia on Prunus sp. (Pashchenko 1988), in India on Prunus cornuta Steudel (Basu and Raychaudhuri, 1976) and in the Himalayas on P. armeniaca and P. mume (Chakrabarti et al. 1970;Chakrabarti and Sarkar 2001). Based on the laboratory tests of Kimura et al. (2016), M. mumecola has the ability to transmit Plum pox virus (PPV) on P. mume, albeit it does not belong to the abundant aphid species in the P. mume orchards in Japan.

Morphological identification
Aphid samples were collected from apricot trees (P. armeniaca) of different ages at six different locations in Hungary in April and May 2020: Győr (lat. 47.699717 N;long. 17.748146 E) and Győrszentiván (lat. 47.676328 N;long. 17.662642 E) in Western Transdanubia, Balatonalmádi (lat. 47.024892 N;long. 18.008070 E) in Central Transdanubia, Budapest Budafok (lat. 47.422926 N;long. 19.032137 E) and Pomáz (lat. 47.652329 N;long. 19.012668 E) in Central Hungary, and Gönc (lat. 48.468253 N;long. 21.282107 E) in Northern Hungary ( Supplementary Fig. 1). The studied home garden trees and orchards received no pesticides. The aphid individuals were collected in 1.5-ml Eppendorf tubes with fine painting brush and were preserved in 70 V/V% ethanol for morphological and in 98 V/V % ethanol for molecular identification. From the collected aphids, 12-48 individuals per site were selected for identification. The aphid specimens were mounted after a 15-min soaking in hot (~ 95 °C) 20% KOH solution. The body content was dissolved by distilled water, and the clean cuticle of the aphids was fixed on microscope slides in glycerine jelly. The morphological identification of the aphid samples was performed by stereo-(Leica MZ6) and light (Euromex iScope 1153-PLi) microscopy, using the keys of Basu and Raychaudhuri (1976) and Blackman and Eastop (1984).

Molecular identification
The morphological identification was supported by the sequence analysis of the DNA mitochondrial cytochrome c oxidase subunit 1 gene (COI fragment) barcode (Hebert et al. 2003;Panini et al. 2017). Six aphid individuals (one specimen per site) were sequenced from the samples collected for morphological identification. Total genomic DNA was extracted from a single aphid with DNeasy blood and tissue kit (QIAGEN, Hilden, Germany) according to the manufacturer's instruction. Amplification of the COI barcode was performed in 20 μl reaction volume containing 20-80 ng DNA, 5 × Phire Reaction Buffer, 0.4 μl of Phire Hot Start II DNA Polymerase (Thermo Scientific, Hungary), 0.2 mM dNTP mix, 3% DMSO, 2.5 μmol of each 5′ and 3′ end primers (LCO 1490: 5′ GGT CAA CAA ATC ATA AAG ATA TTG G 3′ and HCO 2198: 5′ TAA ACT TCA GGG TGA CCA AAA AAT CA 3′ (Folmer et al. 1994)) and sterile distilled water. PCR was carried out in a Swift MaxPro thermocycler (ESCO Healthcare, Singapore). The cycling parameters were as follows: initial denaturation at 98 °C for 30 s, then 30 cycles of denaturation at 98 °C for 5 s, annealing at 49 °C for 5 s, extension at 72 °C for 15 s and a final extension for 1 min at 72 °C. Amplification was verified on a 1% (w/v) ethidium bromide-stained agarose gel in 1 × TBE buffer. Fragment size was estimated by comparison with the 1 kb DNA ladder (Fermentas, Waltham, MA, USA). The amplified fragments were purified using ExoSAP-IT Express purification kit (Thermo Scientific, Hungary) for direct sequencing. Sequencing was performed in an automated sequencer ABI PRISM 3100 Genetic Analyser (Applied Biosystems, Foster City, CA, USA). For each fragment, the nucleotide sequences were determined in both directions.
Forward and reverse sequences were edited and assembled; alignment and neighbour-joining analysis were conducted in MEGAX (Kumar et al. 2018). DNA sequences were verified using the BLASTN algorithm at NCBI.

Life cycle, damage and susceptibility
The development of the M. mumecola colonies was followed on individual trees in home gardens to assess the presence of different stages and the damage caused on the foliage of apricot trees. Susceptibility of different apricot cultivars to M. mumecola was assessed in five apricot orchards in Hungary in July 2020 ( Supplementary Fig. 1 cultivars 'Goldrich', 'Orange Red', 'Budapest', 'Ceglédi bíbor', 'Ceglédi kedves', 'Korai zamatos' and 'Pannónia'). All orchards were planted on cherry plum (Prunus cerasifera Ehrhart) rootstock. Cultivars within each orchard received the same horticultural practices and pest management (Supplementary Table 1). Number of shoots and number of M. mumecola-infested shoots were counted for each sampled tree. The extent of the damage was expressed as percentage of infested shoots on randomly selected trees (n ≥ 9 trees per cultivar) (Supplementary Table 1). Due to the non-normal distribution of the data and lack of homogeneity of variance, for intergroup comparisons the Brunner-Munzel test (two cultivars) or the adjusted rank Welch test followed by Brunner-Munzel test with Bonferroni adjustment (more than two cultivars) were used. As cultivars were grown in nonrandomised plots (spatial pseudoreplication), all statistical comparisons were tested at the significance level of 0.0001 and 0.05. Statistical analyses were performed with the software package ROPstat (Vargha et al. 2015).

Detection of PPV viral RNA
For detection of the presence of PPV in the aphids collected from apricot trees at the six aphid collection sites described above ( Supplementary Fig. 1), six aphid individuals (one specimen per sampling site) were tested. After releasing the viral RNA template from a single aphid with boiling technique (Kim et al., 2016), a two-step RT-PCR was carried out for the detection of PPV with a PPV-specific primer pair (PP3: 5′ TTA TCT CCA GGA (AG)TTG GAG C 3′ and PCI: 5′ TTG AGT CAA ATG G(AG)ACA GTT GG 3′ (Glasa et al. 2002)) amplifying a 836-bp-long fragment corresponding to the 5′-terminal region of the P3 gene (580 nt), the complete 6K1 (156 nt) and a 100 nt portion from the 3′-terminal part of the CI gene [3′P3-6K1-5′CI]. Amplification was performed with Phire Hot Start (Thermo Scientific, Hungary) polymerase as described above. The cycling parameters were as follows: initial denaturation at 98 °C for 30 s, then 30 cycles of denaturation at 98 °C for 5 s, annealing at 55 °C for 5 s, extension at 72 °C for 15 s and a final extension for 1 min at 72 °C. Amplification was verified on a 1% (w/v) ethidium bromide-stained agarose gel in 1 × TBE buffer. Fragment size was estimated by comparison with the 1 kb DNA ladder (Fermentas, Waltham, MA, USA).
For PPV subgroup typing, aliquots of the PCI/PP3 PCR products were subjected to CAPS analysis with restriction endonucleases DdeI, EcoRI and EcoRV (Thermo Scientific, Hungary), as suggested by the manufacturer. Based on the resulting digestion pattern, the three most common strains in Hungary (PPV-D, PPV-M, PPV-Rec) can be differentiated (Glasa et al. 2002;Ádám, et al. 2015).

Morphological identification
We confirmed the presence of M. mumecola in Hungary by morphological characteristics. Myzus mumecola appears to be widespread in Hungary, as we could detect it in all of the examined regions of Hungary: Győr (25 May 2020; number of identified specimens: 48) and Győrszentiván (23 April 2020; 24 specimens) in Western Transdanubia, Balatonalmádi (16 May 2020; 24 specimens) in Central Transdanubia, Budafok (2 May 2020; 12 specimens) and Pomáz (25 May 2020; 12 specimens) in Central Hungary, and Gönc (19 May 2020; 12 specimens) in Northern Hungary. The slide-mounted specimens were deposited in the Hungarian Natural History Museum, Budapest. These are the first records of M. mumecola from Central Europe. We propose the common name 'apricot aphid' for M. mumecola, referring to the host plant, apricot.

Morphological description of Myzus mumecola
Apterous (wingless) M. mumecola individuals in life are medium-sized (2-2.4 mm) (Panini et al. 2017), pale green in colour with pale cauda and siphunculi. The cauda is short and triangular, and the siphunculi are cylindrical, about 7 times longer than their width (Fig. 1a) (Panini et al. 2017). On the mounted specimens, key morphological characters of M. mumecola are clearly discernible: M. mumecola differs from Phorodon humuli Schrank in the lack of the finger-like processes on the antennal tubercles (Fig. 1b), while from Myzus persicae Sulzer in the divergent antennal tubercles (Fig. 1b), longer hairs of the antennae (Fig. 1c), and the narrowing shape of the siphunculi (Fig. 1d) (Basu and Raychaudhuri 1976;Blackman and Eastop 1984). Other aphid species recorded from apricot differ significantly from M. mumecola in size or colour (Blackman and Eastop 1984).

Molecular identification
Identification of M. mumecola was also confirmed by DNA sequence analysis. The COI barcode of six M. mumecola individuals (one specimen per site) was sequenced from the samples collected at six locations in Hungary (Győrszentiván, Győr, Balatonalmádi, Budapest Budafok, Pomáz, Gönc): the sequenced region covered 709 nucleotides. Forward and reverse sequences were aligned, and a consensus sequence was generated for each. All sequences are deposited in the NCBI GenBank under the accession number MT635054-MT635059. For sequence comparison, further M. mumecola (AB738876 and AB738877) and M. persicae (JF883917, HQ971262, GU668755, KR032469, KR039431, KY323048, KP759548, MN320357) COI sequences have been retrieved from NCBI GenBank. For NJ analysis, the sequences were trimmed; hence, finally 644 bases were used in the alignment and cluster analysis. The sequences both for M. mumecola and for M. persicae were uniform within the species, while comparing M. mumecola and M. persicae sequences 42 SNPs were found differentiating the two species. The neighbour-joining analysis (Fig. 2) revealed two independent clades with strong support (100%), one harbouring the M. mumecola individuals and the other the M. persicae individuals.

Life cycle, damage and susceptibility
Fundatrices were observed on the leaves as early as in the initial stages of leaf development, indicating that M. mumecola probably overwinters as eggs on apricot trees. The following 2-3 apterous generations were found to feed on the backside of the leaves (Fig. 3a) and caused significant malformation on young leaves: developing leaves were distorted, severely curled downwards (perpendicularly to the midrib) forming pseudo-galls as reported by Panini et al. (2017) (Fig. 3b). The apterae of the later spring generations were actively moving on the trees, colonising the tips of the growing shoots in the surrounding of the initial colonies. The growth of the infested shoots was affected, and the shoots became slightly distorted (Fig. 3c). In some cases, we observed delayed flowering (in mid-April) and fructification of the infested shoots. In some cases, aphid individuals were found feeding on the growing fruits. The infested leaves partially discoloured and fell prematurely within a couple of weeks. New shoot growth appeared at the base of the dead shoot tips in early June. The first alatae (winged) forms appeared in mid-May in the pseudo-galls (Fig. 3c), indicating that M. mumecola follows a heterocyclic life cycle in Hungary, although its summer host is still unknown (Panini et al. 2017). Alatae feeding on the winter host (apricot) were found by mid-June.

Detection of PPV viral RNA
The RT-PCR carried out with the PP3/PCI PPV-specific primer pair amplified fragments of the expected size (around Fig. 3 Damage of M. mumecola on apricot (a M. mumecola colony on the backside of a leaf; b malformed leaves on a growing shoot; c shoots disrupted in development with alateae forms) Table 1 Proportion of shoots (mean ± SD) infested by M. mumecola on different cultivars in five apricot orchards in Hungary Different capital and lowercase letters indicate significant differences among means on p < 0.0001 and p < 0.05 levels, respectively. Description of the orchards is given in Supplementary Table 1 Site 1

Discussion
In the past few years, Myzus mumecola became an important pest of apricot in Italy (Panini et al. 2017) and, as we found, in Hungary as well. Its rapid expansion across the apricotgrowing regions of Europe is highly probable. Since M. mumecola were not recorded in an extensive faunal survey of the apricot orchards in West Serbia between 2009 and 2011 (Jevremović et al. 2016), the new species had presumably been established in Central Europe in the past few years. As Northern Hungary is one of the northernmost apricot-growing regions of Europe (Bassi 1999;FAOSTAT 2020), it is less likely that continental climate can limit the distribution of the new species (Kaneko 1993). Aphids had been generally uncommon pest of apricot in Central Europe, and therefore insecticide sprays were not required for aphid control (Pénzes et al. 2003). As M. mumecola is now a significant pest of apricot in Central Europe, deeper understanding of its life cycle (including the range of summer hosts and the vector role in the virus and phytoplasma dissemination) and development of new effective control strategies are required. The morphological characteristics of M. mumecola specimens examined by us matched exactly with the description of Panini et al (2017) and the morphological keys of Takahashi (1965), Basu and Raychaudhuri (1976) and Blackman and Eastop (1984) (Fig. 1). The identity of the species was also confirmed by the COI barcode sequence analysis (Fig. 2), which fully matched with those of the Italian and Japanese samples of M. mumecola (AB738876 and AB738877) proving the reliability of the morphological descriptions in the distinction of the species from other aphid pests of apricot.
The 13 apricot cultivars evaluated varied considerably in their susceptibility to M. mumecola. While most of the cultivars showed low or moderate aphid incidence, 'Goldrich' and 'Spring Blush' appeared to be particularly susceptible ( Table 1). The abundance of aphids on fruit trees is primarily limited by the intensity of shoot growing, quality and quantity of plant soap, and the presence of natural enemies (Dixon, 1973). The susceptibility of different plant cultivars depends on several physiological and genetical factors (Dogimont et al. 2010;Pierson et al. 2011). Some cultivars compensate for the loss caused by aphids (Pierson et al. 2011) or restrain the feeding of aphids in a physical way (Giordanengo et al. 2010), while in other cases the expression of different resistance genes leads to some level of tolerance (Dogimont elt al. 2010). In the case of peach (Prunus persica L.), Pascal et al. (2002) found a single dominant resistance gene in the red leaf rootstock cultivar 'Rubira', which is responsible for induced resistance against M. persicae . Aphid resistance genes of apricot have not been identified yet.
Our findings, obtained with CAPS analysis of the PP3/ PCI amplified PCR product which is used for subgroup typing of PPV isolates (Glasa et al. 2002;Ádám et al. 2015), are in accordance with previously published results where the most prevalent strain of PPV in apricot in Hungary was PPV-D (Sihelská et al. 2017). Although M. mumecola has lower transmission rate (12%) than M. persicae (24.4%) (Kimura et al. 2016), it has the potential to be an important vector of PPV because of its high abundance and high rate of aphids carrying PPV.
Funding This research was supported by the National Research Development and Innovation Office of Hungary (K 129311) and by the Ministry for Innovation and Technology within the framework of the Thematic Excellence Programme 2020-Institutional Excellence Subprogram (TKP2020-IKA-12) for research on plant breeding and plant protection. Open access funding provided by Szent István University.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creat iveco mmons .org/licen ses/by/4.0/.