Biodiversity Data Journal :
Research Article
|
Corresponding author: Michael Skvarla (mxs1578@psu.edu)
Academic editor: Therese Catanach
Received: 04 Dec 2019 | Accepted: 17 Jan 2020 | Published: 27 Jan 2020
This is an open access article distributed under the terms of the CC0 Public Domain Dedication.
Citation:
Skvarla M, Kramer M, Owen CL, Miller GL (2020) Reexamination of Rhopalosiphum (Hemiptera: Aphididae) using linear discriminant analysis to determine the validity of synonymized species, with some new synonymies and distribution data. Biodiversity Data Journal 8: e49102. https://doi.org/10.3897/BDJ.8.e49102
|
Although 17 species of Rhopalosiphum (Hemiptera: Aphididae) are currently recognized, 85 taxonomic names have been proposed historically. Some species are morphologically similar, especially alate individuals and most synonymies were proposed in catalogues without evidence. This has led to both confusion and difficulty in making accurate species-level identifications. In an attempt to address these issues, we developed a new approach to resolve synonymies based on linear discriminant analysis (LDA) and suggest that this approach may be useful for other taxonomic groups to reassess previously proposed synonymies. We compared 34 valid and synonymized species using 49 measurements and 20 ratios from 1,030 individual aphids. LDA was repeatedly applied to subsets of the data after removing clearly separated groups found in a previous iteration. We found our characters and technique worked well to distinguish among apterae. However, it separated well only those alatae with some distinctive traits, while those apterate which were morphologically similar were not well separated using LDA. Based on our morphological investigation, we transfer R. arundinariae (Tissot, 1933) to Melanaphis supported by details of the wing veination and other morphological traits and propose Melanaphis takahashii Skvarla and Miller as a replacement name for M. arundinariae (Takahashi, 1937); we also synonymize R. momo (Shinji, 1922) with R. nymphaeae (Linnaeus, 1761). Our analyses confirmed many of the proposed synonymies, which will help to stabilize the nomenclature and species concepts within Rhopalosiphum.
Aphidoidea, Aphididae, taxonomy, agriculture, species delimitation, linear discriminant analysis, Melanaphis arundinariae, Rhopalosiphum nymphaeae
The difficulty of aphid taxonomy and identification has been recognized as far back as the 18th Century by Carl Linnaeus (
Challenges in aphid taxonomy have resulted in different tools and methods to discern species. For example,
Both statistical and non-statistical morphological tools have been developed for use with slide-mounted specimens. For example, online interactive keys have been developed to distinguish among large numbers of aphid species, but they lack statistical grounds for identification (e.g.,
In this paper, we expand upon previous examples using discriminant analysis in aphid taxonomy to test whether multivariate statistics support currently recognized Rhopalosiphum species and their synonymies. Specifically, we use LDA to statistically compare 34 valid and synonymized Rhopalosiphum species to test the validity of historic synonomies. Valid species and synonymies were tested using iterative LDA analyses in which we removed the most distinct species clusters and reanalyzed the remaining species. We used LDA analyses with only valid species and applied the resulting discriminant functions to synonymized species to determine whether the synonymies are statistically correct using species specific LDA functions. The methods and analyses presented here are unique and repeatable with respect to previous aphid studies. No previous studies have tested taxonomic hypotheses in this manner. We use open source software and describe the analyses with sufficient detail to make them repeatable, which has been lacking from the literature. Lastly, our analyses are statistically robust with respect to LDA model assumptions as we thoroughly describe character transformations and missing data and their relationship to model assumptions.
Our example genus of aphids in need of comprehensive review and revision is Rhopalosiphum Koch, 1857 (Aphididae: Aphidinae: Aphidini: Rhopalosiphina) (
Rhopalosiphum species, associated original descriptions, type depositories, and number of specimens examined. Bold names are currently recognized as valid species with synonymized species listed immediately after. Citations marked with an asterisk (*) are not included in the literature cited as the authors were unable to locate them.
Species | Original description | Primary type depository | Secondary type depositories | Number of apterae measured | Number of alatae measured | Notes |
R. albigubernum |
|
- | - | Unable to locate types. | ||
R. arundinariae |
|
USNM | NHMUK | 3 | 5 | Moved to Melanaphis. |
R. cerasifoliae |
|
NHMUK | USNM | 0 | 10 | |
= Aphis furcata |
|
CNC | CNC | 0 | 6 | Lectotype and paralectotypes on one slide. |
= Aphis tahasa |
|
USNM | NHMUK | 27 | 10 | |
R. chusqueae |
|
CZULE | 2 | 0 | ||
R. dryopterae |
|
- | - | Unable to locate types. | ||
R. enigmae |
|
INHS | NHMUK, USNM | 120 | 32 | |
= R. laconae |
|
USNM | NCSU | 54 | 5 | Unable to locate holotype, apparently lost or never deposited. Neotype designated from NCSU paratypes. |
R. esculentum |
|
- | - | Unable to locate types; syn. possible de Aphis craccivora (Remaudiere & Remaudiere 1997) | ||
R. maidis |
|
USNM | 93 | 51 | ||
= Aphis adusta |
|
- | - | Unable to locate types. | ||
= Aphis africana |
|
NHMUK | NHMUK | 11 | 3 | |
= Aphis cookii |
|
EMEC | 0 | 13 | ||
= Stenaphis monticellii |
|
- | - | Unable to locate types. | ||
= Aphis obnoxia | Mordvilko (1916)* | - | - | Unable to locate types. | ||
= Schizaphis setariae | Rusanova (1962)* | - | - | Unable to locate types. | ||
= Aphis vulpiae |
|
- | - | Unable to locate types. | ||
= R. zeae | Rusanova (1960)* | - | - | Unable to locate types. | ||
R. musae |
|
NHMUK | NHMUK | 27 | 7 | Lecto- and paralectotypes. |
= R. scirpifolii |
|
USNM | NHMUK, USNM | 6 | 15 | |
R. nigrum |
|
CNC | CNC, NHMUK, USNM | 10 | 8 | |
R. nymphaeae |
|
58 | 14 | |||
= Aphis aquatica |
|
OSUC (?) | - | - | Unable to locate types. Type depository not given in original description. Jackson worked at The Ohio State and likely deposited his types there, however a search for the specimens did not locate them. | |
= Aphis butomi |
|
- | - | Unable to locate types. | ||
= Aphis infuscata |
|
- | - | Unable to locate types. | ||
= R. najadum |
|
- | - | Unable to locate types. | ||
= R. momo |
|
- | - | Unable to locate types, apparently lost. | ||
= Aphis prunaria |
|
NHMUK | 0 | 2 | ||
= Aphis prunorum |
|
Kiev Entomological Station (defunct) | - | - | Unable to locate types, apparently lost. | |
= Hyadaphis sparganii |
|
NHMUK | NHMUK | 6 | 0 | |
= R. yoksumi |
|
Entomology laboratory, Calcutta University | - | - | Types not examined. | |
R. oxyacanthae |
|
NHMUK | 0 | 9 | ||
= Aphis crataegella |
|
NHMUK | - | - | ||
= Aphis edentula |
|
NHMUK | NHMUK | - | - | Holotypes and 2 paratypes on one slide, holotype indicated. Ovipara. |
= Aphis fitchii |
|
USNM | NHMUK | 0 | 3 | |
= Aphis insertum |
|
USNM | 0 | 54 | ||
= Aphis macatata |
|
NHMUK | 0 | 1 | ||
= R. viridis |
|
CNC | CNC, NHMUK | 0 | 14 | |
= R. mali bivincta |
|
USNM | - | - | Unable to locate types, apparently lost. See also |
|
= R. mali fulviventris |
|
USNM | - | - | ||
= R. mali immaculata |
|
USNM | - | - | Unable to locate types, apparently lost. See also |
|
= R. mali nigricollis |
|
USNM | - | - | ||
= R. mali nigriventris |
|
USNM | - | - | ||
= R. mali obsoleta |
|
USNM | 0 | 1 | ||
= R. mali pallidicornis |
|
USNM | 0 | 1 | ||
= R. mali tergata |
|
USNM | 0 | 1 | ||
= R. mali thoracica |
|
USNM | 0 | 1 | ||
= R. mali triseriata |
|
USNM | 0 | 1 | ||
R. padi |
|
NHMUK | 49 | 44 | Holotype not examined. | |
= Siphocoryne acericola |
|
- | - | Unable to locate types. | ||
= Siphonaphis padi americana |
|
- | - | Unable to locate types. | ||
= Aphis avenaestivae |
|
- | - | Unable to locate types. | ||
= Siphocoryne donarium |
|
- | - | Unable to locate types. | ||
= Siphocoryne fraxinicola |
|
- | - | Unable to locate types. | ||
= Aphis holci |
|
- | - | Unable to locate types. | ||
= Aphis prunifoliae |
|
USNM | 24 | 24 | ||
= Aphis pseudoavenae |
|
CNC | 2 | 11 | ||
= Aphis tritici | Lawson (1866)* | - | - | Unable to locate types. | ||
= Aphis uwamizuskurae |
|
- | - | Unable to locate types. | ||
R. padiformis |
|
CNC | CNC, NHMUK, USNM | 4 | 0 | |
R. parvae |
|
USNM | 3 | 0 | ||
R. rufiabdominale |
|
USNM | USNM | 24 | 28 | |
= Cerosipha californica |
|
EMEC | US Bureau of Entomology and Plant Quarantine (defunct) | 0 | 1 | Unable to locate paratypes, apparently lost. |
= R. fucanoi | Moritsu (1947)* | Entomology Laboratory, Kyusyu Imperial University | - | - | Types not examined. | |
= R. gnaphalii |
|
USNM | 4 | 1 | ||
= Anuraphis mume |
|
- | - | Unable to locate types. | ||
= Yamataphis oryzae |
|
- | - | Unable to locate types. | ||
= Yamataphis papaveri | Taihoku agricultural experiment station (defunct) | - | - | Unable to locate types. | ||
= Pseudocerosipha pruni | Shinji (1932)* | - | - | Unable to locate types. | ||
= Aresha setigera |
|
- | - | Unable to locate types. | ||
= Aresha shelkovnikovi |
|
- | - | Unable to locate types. | ||
= Siphocoryne splendens |
|
NHMUK | 40 | 3 | ||
= R. subterraneum |
|
USNM | 17 | 19 | ||
R. rufulum |
|
CNC | CNC, NHMUK | 10 | 10 | |
R. sanguinarium |
|
NHMUK (cotype) | 0 | 1 | ||
sp. nov. "ex. Arisaema" | undescribed | USNM | 34 | 1 | ||
R. sp. T | undescribed, proposed by |
Crop and Food Research, Lincoln, New Zealand | - | - | ||
R. near insertum | undescribed, proposed by |
Crop and Food Research, Lincoln, New Zealand | - | - | ||
R. sp. x | undescribed, proposed by |
- | - |
Rhopalosiphum hosts and distributions. Host author names have been omitted for space. Unless otherwise noted, information is compiled from the original description (see Table
Taxon | Primary host | Secondary host | Distribution | Notes | Additional references |
R. albigubernum | Citrus | Baidicheng, Fengjie, Sichuan, China | Known only from the type series, which consists of two alate vivipara. | ||
R. arundinariae | Arundinaria tecta | Gainesville, FL | Known only from the type series, which were found as dense colonies on ventral side of younger A. tecta leaves that contained aptareae and alatae on 16 April 1930. | ||
R. cerasifoliae | Prunus pennsylvanica, P. virginiana | Cyperaceae, including Schoenoplectus, Scirpus, and Eleocharis. Also Juncus (Juncaceae) | North America, wherever host plants occur | Can persist throughout summer on primary host. | |
R. chusqueae | Chusquea tomentosa | Costa Rica | Lives close to the nodes and well protected by leaves of C. tomentosa, so not easily detectable. Alatae and lifecycle unknown | ||
R. dryopterae | Dryopteris filix-mas | Kyrgyzstan | |||
R. enigmae | Typha, especially. T. latifolia; also recorded from Sparganium | North America, wherever host plants occur, though apparently most common in the East (A. Jensen, pers. comm.). "R. laconae" is morphology restricted to southeastern coastal plain. | Autoecious on cattails, no primary host known. | ||
R. maidis | Prunus, including P. cornuta (Pakistan), P. mume and P. persica (Korea), and P. sargentii (Japan). | Zea, Sorghum, Hordeum,other Poaceae; occasionally Cyperaceae and Typhaceae. | Cosmopolitan, but cannot survive outdoors in regions with severe winter climates | Most lineages are autoecious anholocyclic, with males and ovipara occurring only rarely. Heteroecious, androcyclic populations are known from Asia, where the species originated. In the Northern Hemisphere, lineages with different chromosome counts are associated with different hosts: 2n = 10 colonize Hordeum and Echinochloa crus-galli, while 2n = 8 colonize Zea and Sorghum; lineages with 2n = 9 , 2n = 11, and heterozygous 2n = 8 are also known. In contrast, 2n = 8 and 2n = 9 lineages in Australia do not exhibit host preference. Morphological and early genetic investigations of North American lineages reported incomplete morphological separation and lack of genetic differences between the various parthenogenetic karyotype lineages of R. maidis, though future investigations may raise one or more lineages to species from within the complex. |
|
R. musae | Prunus beseyi, P. subcordata, P. fasciculata | Scirpus; also Musa, Ensete, Sterlitzia | WA to CO, MD (native); Europe, Middle East, Africa, and Australia (adventive) | Adventive populations outside of the native range are presumed to be anholocyclic. Specimen vouchers reportedly sent to the USNM are not present in the collection and were presumably never sent. |
|
R. nigrum | Crataegus | Avena sativa, Zizania aquatica, Alisma | ON, MB, AK; unconfirmed reports from OR, UT |
|
|
R. nymphaeae | Prunus | Aquatic and semi-aquatic plants, including Alisma, Juncus, Nelumbo, Nuphar, Sagittaria, Sparganium and Typha. Occasionally other hosts, including Canna, Glyceria, Lactuca, Triticum, and Tulipa | Cosmopolitan | ||
R. oxyacanthae | Various Rosaceae including Malus, Pyrus, Cotoneaster, Crataegus, Sorbus, and Prunus | Various Poaceae, including Agropyron, Agrostis, Festuca, Poa; occasionally cyperaceae and Juncaceae. | North America (probable native range), Europe, North Africa, Japan | The name R. insertum Walker was used widely for the species in North America, but |
|
R. padi | Primarily Prunus virginiana (North America) and P. padus (Europe), occasionally other Prunus. One record on Syringa vulgaris. | Primarily Poaceae, occasionally Asteraceae, Brassicaceae, Cyperaceae, Iridaceae, Juncaceae, Lilaceae, Typhaceae and other taxa. | Cosmopolitan | Can persist throughout summer on primary host. | |
R. padiformis | Poa, Triticum | BC, MT | Primary host and associated morphs not known; alate males have been obtained in culture. | ||
R. parvae | Carex (US), Scirpus lacustris (Italy) | Native to North America (IL); Italy (adventive) | Primary host not known. | ||
R. rufiabdominale | Primarily Prunus; also recorded from Malus, Chaenomeles, Pyrus, Rhodotypos, and Sorbus. | Underground parts of Poaceae (including cereals) and Cyperaceae. Can infest some dicots (e.g. Asteraceae, Solanaceae), especially in greenhouse and hydroponic situations. | Native to East Asia. Currently pan-tropical/subtropical and restricted to greenhouses in colder climates. | Heteroecious holocyclic in East Asia and Italy; autoecious anholocyclic in majority of introduced range. |
|
R. rufulum | Crataegus | Primarily Acorus, also recorded from Typha | North America and Europe, wherever hosts are found. | ||
R. sanguinarium | Crataegus mexicana | Unknown, but reared on various Poaceae in lab conditions | Mexico | ||
sp. nov. "ex. Arisaema" | Unknown | Arisaema | Maryland | ||
R. sp. T | Unknown | Cereals | New Zealand | ||
R. near insertum | Unknown | Cereals | New Zealand; Victoria, Australia | ||
R. sp. x | Unknown | Zea mays | Victoria, Australia |
The taxonomic history of Rhopalosiphum is complicated. Many taxa with slightly-swollen siphunculi were historically included within Rhopalosiphum, but are now placed in a different tribe, Macrosiphini (e.g. Hyadaphis Börner, Lipaphis Mordvilko, Rhopalomyzus Mordvilko), or other aphidine genera (e.g. Melanaphis van der Goot, Schizaphis Börner) (
Molecular studies (e.g.
Rhopalosiphum species concepts are further complicated by the presence of cryptic species (
Finally, while 17 species are currently recognized, at least 85 specific names (which includes misspellings of valid names) have been applied to Rhopalosiphum (
Herein, we confirm many historic synonymies using LDA, propose two new synonymies, and report new distribution data for R. rufulum Richards, 1960, based on material examined for the analyses.
Most specimens examined are housed in the National Museum of Natural History Aphidomorpha Collection (USNM), which is currently housed at the USDA-ARS Beltsville Agricultural Research Center in Beltsville, Maryland, U.S.A. Additionally, types and other material were borrowed from the Essig Museum of Entomology (EMEC), Berkeley, California, U.S.A.; the Florida State Collection of Arthropods (FSCA), Gainesville, Florida, U.S.A.; the Illinois Natural History Survey Insect Collection (INHS), Champaign, Illinois, U.S.A.; the North Carolina State University Insect Museum (NCSU), Raleigh, North Carolina, U.S.A.; The Ohio State University Triplehorn Insect Collection (OSUC), Columbus, Ohio, U.S.A; the personal collection of Andrew Jensen (AJ), Lakeview, OR, U.S.A.; The Canadian National Collection of Insects, Arachinds, and Nematodes (CNC), Ottawa, Ontario, Canada; the Natural History Museum (NHMUK), London, United Kingdom; and the Aphidological Collection of the University of León (CZULE). For a list of types, species, and number of specimens examined, see Table
Slides were labelled with individual, sequential numbers (MS 0001–0980) and specimens assigned a number that was appended to the slide number (e.g. MS 0001-1 for the first specimen on the first slide). Specimens were examined using a Zeiss Axio Imager M1 stereomicroscope; micrographs were taken and measurements of various morphological features of adult female apterous and alate specimens made using AxioVision 4.9.1 software (Carl Zeiss AG, Oberkochen, Germany) (Fig.
Measurement number | Measurement description or ratio |
1 | Antenna segment (AS) 1, length |
2 | Antenna segment (AS) 2, length |
3 | Antenna segment (AS) 3, length |
4 | Antenna segment (AS) 4, length |
5 | Antenna segment (AS) 5, length |
6 | Antenna segment (AS) 6 base, length |
7 | Antenna segment (AS) 6, process terminalis (pt), length |
8 | Head width across eyes |
9 | Ultimate rostral segment (RIV+V), length |
10 | Ultimate rostral segment (RIV+V), width |
11 | Hind femur, length |
12 | Hind tibia, length |
13 | Hind distitarsus, length |
14 | Siphunculus, length |
15 | Cauda, length |
16 | Abdominal segment 8 submedian seta, length |
17 | Body length (BL) |
18 | Siphunculus, width |
19 | RIV+V length: RIV+V width |
20 | RIV+V length: hind distitarsus length |
21 | Pt length: AS 6 base length |
22 | Siphunculus length: siphunculus width |
23 | Siphunculus length: cauda length |
24 | Siphunculus length: AS 3 length |
25 | RIV+V length: cauda length |
26 | BL: head width |
27 | AS 3 and 4 fused (0 = no, 1 = yes) |
28–33 | Angle, A1–A6 (see Fig. |
34–36 | Area, S1–S3 (see Fig. |
37–57 | Wing length, L1–L21 (see Fig. |
58 | L1:L2 (37:38) |
59 | L1:L3 (37:39) |
60 | L1:L4 (37:40) |
61 | L1:L5 (37:41) |
62 | L1:L6 (37:42) |
63 | L1:L7 (37:43) |
64 | L2:L3 (38:39) |
65 | L2:L4 (38:40) |
66 | L2:L6 (38:41) |
67 | L2:L7 (38:42) |
68 | L3:L4 (39:40) |
69 | L6:L7 (42:43) |
Morphological terms and structures were adapted from
We used linear discriminant analysis (LDA) in R (
The analysis proceeded in several steps, which were similar for apterae and alatae, the first of which was data cleaning, described below. LDA is one of the most commonly employed discriminant function analyses, used both to identify useful characters distinguishing specimens of valid species and groups of species, and to form discriminant functions that could then be applied to specimens of uncertain taxonomic status to determine if they cluster with valid species or cluster separately, which would suggest a new species. This LDA approach was applied in a systematic iterative manner; the initial linear discriminant functions separated the most distinctive species, leaving a large amorphous cluster of the other known species. In the next step, the distinctive valid species were removed in that iteration, and the method applied to the remaining valid species. This could be followed by another iteration, until all valid species were separated to the extent possible using available morphological characters.
Prior to the analyses, the dataset was checked for incorrectly entered data. This was performed for both the valid and synonymized by constructing boxplots for each trait, broken down by species (e.g. Fig.
We standardized the size measurements of each specimen, based on a ‘size’ variable (but also kept ‘size’ as a variable in the LDA, explained below). There are several reasons for this. One is that some species are more variable in size than others. If measures are not size adjusted, the LDA does not work as well; in fact, we implemented the size adjustments because size affected the usefulness of many of the measures, especially for taxa with a lot of size variability. Second, if one doesn't adjust for size, many of the measures are correlated through size, so less useful. Third, we wanted measures which made sense morphometrically, and these are often 'relative' measures, (e.g. antenna are relatively long for species A (in relation to its size) versus species B). If not adjusted for size, we would have to use ratios for more variables. Standardizing on size was very helpful for many traits and might be useful when developing criteria to separate species in other groups. Given that it is desirable to adjust for size, how does one best estimate 'size'? There is not a one-size-fits-all solution for this and, while we found one that seems to work well for these taxa, it may be improved upon or altered for other taxa. This variable was constructed by combining the body length, head width, and femur length using principal components (PC); head width and femur length were chosen as, among all measured characters, these were most highly correlated with body length. This is a dimension reduction technique, the idea being to create a single variable (the first PC) that best captures the variation in these three correlated measures. In cases where one or two of these measures were missing, the derived principal component measure (henceforth labelled ‘relative size’) for the specimen was imputed using linear regression based on the rest of the data set. Another method employed for developing a size measure is the use the geometric mean of the characters. We calculated the geometric means for data where we had all three characters and found the correlation between the geometric mean and first PC to be 0.9963; for this data set the two methods would yield essentially identical results.
Our relative size measure was retained as a character in the LDA analysis. It was also used to adjust all other size measures using linear regression, i.e. adjusted measures were residuals from regressing each of the size measures (dependent variable) on the relative size (independent variable). This resulted in smaller and larger individuals of the same species having similar adjusted measurements. Non-size measures (such as wing angles or ratios) were not adjusted using relative size, instead they were transformed by taking logs; this created measures that were closer to being normally distributed. All measures were then individually standardized (using all samples) to mean zero, standard deviation one (this helps one to interpret the coefficients of the linear discriminant analysis results). Missing values were then imputed by randomly sampling from the corresponding adjusted measure of other individuals of that same species. In the unusual case where this could not be done (e.g. all specimens were incomplete for this trait), the trait value was set to zero (which is the overall mean for each measure), removing its influence on the calculated discriminant functions. If all specimens of a species were naturally lacking a trait, a new character column was created for that trait, with either 0 (not missing) or 1 (missing). The end result was that a naturally missing trait could be used as a character when creating the linear discriminant function and that incomplete or aberrant individuals were not dropped from the analysis because of missing data.
A first LDA was performed on the cleaned dataset for only the valid species (this included the three variables used to create relative size, as well as the adjusted size measure) using 49 measured values and 20 ratios calculated from those values. The first three latent axes were sufficient to explain 80% or more of the variability. We looked at which variables loaded most heavily on each axis and compared that to the boxplots created at an earlier step as a check that the methodology was working as expected. The linear discriminant functions derived from the valid species were then applied to all specimens (valid and synonymized) so the specimens could be mapped to the two dimensional space created from pairs of the latent axes. Species were suitably coded so they could be distinguished on the plots. Since there was only one individual of R. sanguinarium Baker, 1934 for both apterae and alatae, it was analyzed with the synonymized, rather than the valid, species.
We then identified clusters of specimens comprising a valid species. Sometimes a valid species was well separated from other valid species and sometimes not. For those well-separated valid species, we looked to see if any of the synonymized species occupied the same space and manually outlined the group. When this happened, we concluded that the synonymized species and valid species were likely the same and we removed them from further analyses. If a synonymized species formed a distinct cluster, we concluded that it was not synonymized with any of the valid (or other synonymized) species and considered it a separate species and also removed it from further analyses.
With the remaining species, we repeated the previously described procedure, producing another LDA using only the valid species that did not separate well in the previous analysis, and using the same independent variables. With fewer species and the same number of independent variables, the software sometimes had difficulty due to insufficient degrees of freedom. When this occurred, we reduced the number of independent variables by removing those that were highly correlated with others in a stepwise manner until there were no independent variables that were highly correlated with other independent variables, and then produced an LDA with the reduced set of independent variables. For alatae, this methodology needed to be repeated a third time to finish separating all the known species.
Rhopalosiphum dryopterae Kan, 1986.
Rhopalosiphum esculentum Raychaudhuri and Roychoudhuri, 1978.
Rhopalosiphum momo Shinji, 1922.
(40) Rhopalosiphum momo SHINJI, n. sp. / pp. 791
Characteristics: Body green or pale. Antennae somewhat longer than body, III shorter than IV and V taken together, with 16–18 subcircular sensorial, IV and V subequal in length, the former with one sensorium in the middle and the latter with a subapical one; flagellum of VI about 3 times as long as the base. Antennae as a whole infuscated and each segment has a few hairs. Cornicles with a basal 2/3 part green or pale, the remaining 1/3 part swollen and infuscated.
Host plant: Prumus [sic] persica L.
Date of collection: 2 June 1920
Locality: Miyakonojo, Nakajo (Nagano Pref.)
Rhopalosiphum arundinariae (Tissot, 1933).
We were unable to locate apterae of R. cerasifoliae (Fitch, 1855) to include in the analyses. However, we included specimens of the junior synonym R. tahasa Hottes, 1950 in order to understand what might happen when synonymized species do not have a presumed conspecific valid species for comparison.
In total, 1,030 Rhopalosiphum specimens (625 apterae, 405 alatae), representing 34 valid and synonymized species, were measured (Table
Wing angle measurements for species of Melanaphis, Rhopalosiphum, and Schizaphis. Ranges are followed parenthetically by the average for each measurement. Generic measurements, which are listed in bold for easy comparisons, were calculated by adding together species in each genus. Rhopalosiphum arundinariae
Genus species | Original description | Specimens measured | A1 | A2 | A3 | A4 | A5 | A6 |
Melanaphis |
|
62 | 25.2–42.6 (30.6) | 92.9–118.7 (106.9) | 44.9–60.6 (52.6) | 37.9–53.1 (43.5) | 27–48.1 (36.8) | 36–53.1 (42.7) |
M. bambusae |
|
9 | 27–34.7 (31.6) | 92.9–104 (101) | 55.6–60.6 (57.9) | 41.1–53.1 (45.4) | 34.1–47.1 (41.3) | 41–53.1 (45.7) |
M. donacis |
|
4 | 27–30.8 (29.6) | 108.7–113.9 (111.1) | 44.9–49 (46.9) | 42.6–45 (43.4) | 27–30.2 (28.7) | 36–44.7 (41.5) |
M. japonica |
|
2 | 32–32 (32) | 111.9–111.9 (111.9) | 54.6–56.7 (55.7) | 42.6–45.2 (43.9) | 30.1–31.7 (30.9) | 40.6–43.3 (42) |
M. pyraria |
|
10 | 27.8–35.1 (32.1) | 106.9–118.7 (112.6) | 46–53.8 (49.7) | 42–47 (44.3) | 27.9–38.3 (33.4) | 38.9–48.6 (44.9) |
M. sacchari |
|
36 | 25.2–42.6 (29.7) | 94.8–112.1 (105.7) | 45.6–56.3 (52.4) | 37.9–48.2 (42.8) | 29.9–48.1 (37.9) | 37.4–49.7 (41.6) |
M. sorini |
|
1 | 29.5–29.5 (29.5) | 113.7–113.7 (113.7) | 48.2–48.2 (48.4) | 42.8–42.8 (42.8) | 33.4–33.4 (33.4) | 42.2–42.2 (42.2) |
R. arundinariae |
|
5 | 29.5–34 (31.5) | 102.2–114.1 (108.5) | 45.3–55.7 (48.8) | 43.8–45.6 (44.5) | 35.1–37.5 (36.5) | 44.7–49.1 (47.1) |
Rhopalosiphum |
|
342 | 26.6–44.6 (34.6) | 96.7–123.5 (113.1) | 28.8–44.6 (36.1) | 38.7–55.1 (47.1) | 15.2–47.4 (25.9) | 37.3–62.8 (50.3) |
R. cerasifoliae |
|
23 | 29.7–40.2 (33.6) | 106.3–122 (114.2) | 33.4–41.2 (37.8) | 41.8–48.8 (45.3) | 17.7–42.1 (26.7) | 43.7–52.6 (47.8) |
R. enigmae |
|
38 | 26.6–38.3 (32.3) | 110.3–123.5 (115.8) | 33–43 (37.6) | 39.9–51.8 (45.3) | 19.9–47.4 (30.4) | 37.3–53.7 (48.6) |
R. maidis |
|
53 | 28.1–44.6 (34.1) | 101.7–123.2 (112.7) | 28.9–40.5 (35.2) | 41.9–54.3 (46.3) | 15.6–36.9 (24.9) | 40.8–60.1 (49.7) |
R. musae |
|
15 | 30.6–39.2 (34.5) | 111.2–119.6 (116) | 31.1–36.3 (34.1) | 42.6–50.8 (46) | 15.2–25.3 (19.9) | 43.4–53.4 (49.2) |
R. niger |
|
8 | 29.2–35.2 (31.7) | 112.8–122.1 (116.4) | 36.6–42 (38.5) | 44.5–49.2 (46.9) | 21.5–32.6 (26.8) | 48.4–55.5 (51.2) |
R. nymphaeae |
|
14 | 29.6–40.8 (34.8) | 105.3–118 (112.4) | 35.3–44.4 (40.3) | 44.3–52 (47.6) | 19.6–46.7 (25.7) | 46.1–58 (49.1) |
R. oxyacanthae |
|
73 | 30.7–42.6 (35.9) | 105.5–120.9 (112.6) | 32.6–44.6 (36) | 42.4–51.5 (47.5) | 17–41.3 (25.6) | 42.3–58.6 (50.6) |
R. padi |
|
72 | 27–43.3 (35.7) | 96.7–123.1 (110.9) | 28.8–41.1 (35.1) | 38.7–52.6 (48) | 17.5–35.7 (23.6) | 42.9–62.8 (51.7) |
R. rufiabdominale |
|
43 | 28.5–39.7 (34.4) | 105.5–120.3 (113.8) | 29.6–40.1 (35) | 43.7–55.1 (48.9) | 18.6–36.1 (27.6) | 45.6–58.5 (51.6) |
R. sanguinarium |
|
1 | 40.9–40.9 (40.9) | 101.3–101.3 (101.3) | 39.8–39.8 (39.8) | 44.6–44.6 (44.6) | 29.6–29.6 (29.6) | 47.8–47.8 (47.8) |
Schizaphis |
|
99 | 26.3–41.5 (31.7) | 100–123.4 (112.4) | 29.5–44.9 (38.8) | 38.1–52.2 (44.2) | 18.5–53.8 (37.9) | 33.6–58 (46.3) |
S. caricis |
|
2 | 28.6–33.8 (31.2) | 111.1–117.9 (114.5) | 37–40.2 (38.6) | 46.3–46.3 (46.3) | 29.4–29.4 (29.4) | 51.4–51.4 (51.4) |
S. graminum |
|
41 | 28.4–35.5 (31.1) | 103.7–118.5 (111.9) | 36.3–44.9 (40.6) | 38.1–44.6 (41.4) | 32.7–53.8 (39) | 40.8–51.1 (44) |
S. minuta |
|
1 | 31.5–31.5 (31.5) | 115–115 (115) | 41.7–41.7 (41.7) | 42.3–42.3 (42.3) | 38.9–38.9 (28.9) | 45.8–45.8 (45.8) |
S. muhlenbergiae |
|
2 | 27.3–27.3 (27.3) | 109.6–109.6 (109.6) | 35.7–43.4 (39.6) | 38.25–40.9 (39.6) | 41.1–44.6 (42.9) | 43.7–49.3 (46.5) |
S. nigra |
|
44 | 26.3–41.5 (31.7) | 100–123.4 (113.5) | 29.5–42.1 (37.5) | 43–52.2 (46.5) | 18.5–46.8 (37.5) | 33.6–58 (47.9) |
S. palustris |
|
3 | 30–38.5 (34.1) | 105.3–112.1 (107.7) | 33.7–39.5 (37.1) | 46.5–49.2 (48.3) | 27.8–43.8 (35) | 46.8–53 (49.5) |
S. rotundiventris |
|
6 | 30.6–40.4 (35) | 107–115.8 (110.8) | 33.6–37.7 (36.3) | 44.5–50.6 (47.2) | 31.9–38.9 (35.5) | 46.2–53.6 (48.7) |
The type material of R. momo is apparently lost, so, left with only the description, we concur with
1) Rhopalosiphum nymphaeae is the only species in the genus with siphunculi that are green basally and dark and expanded apically, as was described in R. momo.
2) Seven species of Rhopalosiphum utilize Prunus as a primary host (Table
a) The body of R. momo is described as “green or pale”. Rhopalosiphum padi and R. rufiabdominale have a distinctive red patch between the siphunculi and R. oxyacanthae have dark green stripes that would probably have been noted in the description if present
b) The antennae of R. momo are “whol[ly] infuscated”. Antenna segment III of R. maidis is pale, rather than concolorous with the other dark segments, and can be pale or dark in R. padi.
3) One potential complication is that
The morphological characters suggested by
Pursuant to Article 57 of the International Code of Zoological Nomenclature (
While sorting undetermined Rhopalosiphum specimens in the USNM collection, three slides of specimens collected from Acorus in North America were discovered. The specimens match the measurements and brief descriptions of European R. rufulum apterous viviperae (
CANADA: 3 female apterae, locality unknown (label states “at HO-17214”), ex Acorus calamus, 29-IX-1952, J. Adams, USNM; UNITED STATES: Massachusetts: 2 female apterae, Hampshire Co., Amherst, ex. Acorus, 31-V-1964, M. Smith, USNM; New York: 3 nymphs, Suffolk Co., Greenport, ex. Acorus calamus, 25-VI-1963, R. Latham, USNM.
In the first LDA of apterae, R. maidis and its synonom R. africana (Theobald, 1914) formed a distinct cluster in the plots of LD1 x LD2 and LD1 x LD3 (Fig.
In the second linear discriminant analysis of apterae, R. padi and its synonyms R. prunifoliae (Fitch, 1855) and R. pseudoavenae (Patch, 1917) formed distinct clusters in both the LD1 x LD2 and LD1 x LD3 plots (Fig.
In the first linear discriminant analysis of alatae, R. maidis formed a cluster with its synonyms R. cookii (Essig, 1911) and R. africana and R. nymphaeae clustered with its synonym R. prunaria (Walker, 1848) in the plots of LD1 x LD2 and LD1 x LD3 (Fig.
In the second linear discriminant analysis of alatae, R. rufulum (along with one of the outlier R. africana not removed from the analyses) formed a distinct cluster in the LD1 x LD2 plot (Fig.
None of the remaining species formed distinct clusters in the third linear discriminant analysis (Fig.
The synonymization of R. momo with R. nymphaeae and movement of M. arundinariae (Tissot, 1933) to Melanaphis brings the total number of described Rhopalosiphum species to 17, two of which (R. dryopterae and R. esculentum) are questionably assigned to the genus, pending examination of the type material, with an additional four undescribed species known.
The discovery of R. rufulum on Acorus in North America and the presence of an undescribed Rhopalosiphum species on Arisaema in eastern North America echoes the sentiment of
Determining species boundaries and synonymies using molecular techniques has become the standard by which most modern taxonomy and systematics are measured. Indeed, with the continually lower costs and greater ease with which highly degraded DNA can be extracted and sequenced from historic museum specimens (
In general, the linear discriminant analyses were successful when applied to apterae. The analyses successfully grouped all synonymized species with their associated valid species over two iterations and we concluded that most of the synonymizations we tested are sound. This demonstrates that linear discriminant analyses can be used to test synonymizations when DNA is unavailable and provides a new method to examine and use historic, slide-mounted specimens. Additionally, R. sp. nov. "ex. Arisaema" formed a distinct group when it was included as a valid species. This supports its status as a valid, but undescribed, species and demonstrates that linear discriminant analysis can be used to distinguish potentially undescribed species from valid described species based on morphological similarities.
However, a few notable problems exist. First, R. parvae and R. rufulum consistently clustered together, but never as a distinct group away from other species. This may indicate that they are synonymous, although several factors indicate they are distinct species: they feed on different secondary hosts (Carex and Acorus, respectively), have non-overlapping ranges (Illinois versus New England and adjacent Canada), and have several morphological differences that separate them (Table
Morphological differences between R. parvae and R. rufulum. Measurement and ratio ranges are followed parenthetically by the mean and number of specimens measured.
Species | ant-6 | caud-l | bl | rost-l:rost-w | ant-pt:ant-6 | bl:head | bl:ant-pt |
R. parvae | 61.7–66.8(64.7, 3) | 91.8–103.1(98.4, 3) | 1410.2–1437.0(1427.9, 3) | 2.13–2.16(2.15, 2) | 3.37–3.77(3.56, 3) | 3.48–3.84(3.60, 3) | 5.80–6.90(6.22, 3) |
R. rufulum | 77.1–91.0(83.6, 9) | 94.9–131.6(119.7, 9) | 1829.0–2186.0(1969.3, 10) | 1.37–2.06(1.78, 9) | 2.42–3.45(3.14, 9) | 4.15–5.04(4.61, 10) | 6.65–9.71(7.50, 9) |
Unfortunately, we were unable to include verified apterous R. cerasifoliae in the analyses due to lack of specimens. We included specimens of R. tahasa, which is synonymized with R. cerasifoliae, in the analyses to see if they would still form a distinct cluster without valid R. cerasifoliae to act as a guide in the discriminant function. Instead, the R. tahasa specimens clustered well with R. padi, which it is not synonymized with, rather than forming a distinct cluster. This association has not previously been suggested in the literature, but without the inclusion of R. cerasifoliae, it is impossible to determine if R. tahasa should instead be synonymized with R. padi. This issue suggests that it is extremely important to include all valid species when creating the discriminant functions, so that synonymized species can be properly plotted. Additionally, assuming that R. tahasa is synonymous with R. cerasifoliae, it suggests that valid species analyzed with the synonymized species (either because they are incorrectly synonymized with a different species or because, as in this case, a synonymized species is included without specimens of its associated valid species) may not form distinct clusters.
The analyses of alatae were less decisive. Intuitively, species with the most distinctive apterae – R. maidis, R. nymphaeae, R. enigmae, and R. rufiabdominale – also had the most distinctive alatae and formed distinct clusters with their associated synonymized species in the first two LDA. However, the remaining species failed to form distinct clusters, even after a third LDA. Additionally, a number of synonymized species – R. furcata, R. fitchii, R. insertum, R. mactata, R. mali, and R. viridis – were only represented by alate specimens and did not cluster with the species with which they are synonymized. Based on our analyses, we cannot confirm that these synonymies are correct, although we also do not propose any be raised as valid species pending additional investigation.
Aphidologists have found that Rhopalosiphum alatae captured without host plant data (e.g., in suction traps) are difficult to identify to species due to similar, conserved morphologies. While published keys to alatae are available (i.e.
The loadings (coefficients) for each character on the discriminant functions could conceivably be used to create taxonomic keys. However, we do not suggest following that path for this group for a number of reasons: measures we used were adjusted by size, based on our data, so this same adjustment would need to be made for all measures (except wing venation angles); variables were then standardized to mean zero, variance one, so each variable would need to be standardized in the same way we did (again, the standardization we used was based on our data); the functions are built largely on continuous characters, rather than naturally dichotomous ones and all characters contribute to each discriminant function. However, since the variables were standardized, their importance (i.e. weighting) to the discriminant function can be evaluated, based on the absolute value of the coefficient. One could determine a cut-off for which variables were important. We looked at this for all the discriminant functions and found that, for any reasonable cut-off, some discriminant functions were largely determined by just a few characters (which would be useful for creating a key), but many were largely determined by at least 12 (which would be less useful). Since placement in the plane for any specimen is determined by two discriminant functions, many with large contributions from many characters, creating a usable key from these results would not be an easy exercise. One possible statistical approach to creating keys from continuous characters that should be explored is the use of regression trees.
Two species of Rhopalosiphum were moved or synonymized, bringing the total number of species in the genus to 17. While LDA has previously been used to distinguish between a limited number of species or ecotypes, the use of it to confirm previously proposed synonymies using historic slide-mounted specimens that lack DNA is novel and yielded promising results. In particular, the analyses confirmed most synonymizations when apterae were analyzed. However, while the most distinct alate Rhopalosiphum and associated synonymies were recovered in the LDA, many species and associated synonymies were not recovered as distinct. The failure of the analyses with some of the alatae using phenotypic traits mirrors problems previously documented for this morphologically similar group.
We thank the Essig Museum of Entomology, Florida State Collection of Arthropods, Illinois Natural History Survey Insect Collection, North Carolina State University Insect Museum, Canadian National collection of Insects, Arachnids, and Nematodes, Natural History Museum, Aphidological Collection of the University of León, and Andrew Jensen for kindly loaning specimens; Lu Musetti for searching the OSCU for the type specimens of Aphis aquatica. Mention of trade names or commercial products in this publication is solely for the purpose of providing specific information and does not imply recommendation or endorsement by the USDA; USDA is an equal opportunity provider and employer.
Length measurements (given in μm) and wing angles for alatae of valid Rhopalosiphum species.
Length measurements (given in μm) and wing angles for apterae of valid Rhopalosiphum species.
Length measurements (given in μm) and wing angles for alatae of synonymized Rhopalosiphum species.
Length measurements (given in μm) and wing angles for apterae of synonymized Rhopalosiphum species.