Phylogeographic analysis of invasive Asian earthworms (Amynthas) in the northeast United States

Phylogeographic studies are useful in reconstructing the history of species invasions, and in some instances can elucidate cryptic diversity of invading taxa. This can help in predicting or managing the spread of invasive species. Among terrestrial invasive species in North America, earthworms can have profound ecological effects. We are familiar with the centuries‐old invasions of European earthworms (Lumbricidae) and their impacts on nutrient cycling in soils. More recent invasions by Asian earthworms of the family Megascolecidae are less fully understood. We used data for two mitochondrial gene fragments, cytochrome oxidase I (COI) and 16S rRNA, to examine the relationships among populations of Asian earthworms in the megascolecid genus Amynthas in the northeast United States. Recent reports have indicated that one species in particular, Amynthas agrestis, is having detrimental effects in mixed forest ecosystems, and we were interested in understanding the invasion history for this species. We were surprised to discover three divergent mitochondrial lineages of Amynthas occurring sympatrically in upstate New York. Given the gap between intra‐ and inter‐lineage sequence divergences, we propose that these three lineages represent cryptic species of Amynthas, one of which is A. agrestis. For each of the three lineages of Amynthas, we observed shared haplotypes across broad geographic distances. This may reflect common origins for populations in each lineage, either by direct routes from native ranges or through post‐introduction spread by natural dispersal or human‐mediated transport within North America. Management efforts focused on horticultural imports from Asia, commercial nurseries within the USA, and on prohibition of bait disposal may help to reduce the further invasion success of Amynthas.     

New insights into the mitochondrial phylogeny of the whitefly Bemisia tabaci (Hemiptera: Aleyrodidae) in the Mediterranean Basin

The whitefly Bemisia tabaci is vector of plant infecting viruses and it is considered as one of the most important agricultural pests around the Mediterranean basin. At present, five biotypes of B. tabaci have been reported in the Mediterranean Basin: B, Q, S, T and M. To establish the phylogeographic relationship of these Mediterranean biotypes with others, 54 samples collected in Europe and Africa were analysed by sequencing the mitochondrial cytochrome oxidase I gene (mtCOI). The phylogeny showed that Spanish samples corresponding to the biotype S were related to the haplotype Uganda 2 of the African clade, associated with recent epidemic upsurges of cassava mosaic virus (CMD) in that country. This phylogeographic relationship gave support to a distinct subgroup revealed within the African clade. Bemisia tabaci collected from Euphorbia plants in Italy (biotype T) formed one of the three distinct subgroups existing within the Southeast/Far East Asian clade, while samples from Turkey (biotype M) clustered together with reference mitochondrial sequences from whiteflies from Pakistan and Thailand. Recent reports indicate that Bemisia populations corresponding to the biotypes S and T are distributed in areas larger than those initially delimited. Other results indicated that samples collected in Sudan grouped within the Mediterranean–North Africa clade together with reference sequences of the biotype Q corresponding to insects collected in Spain and Morocco. Mitochondrial haplotypes of B. tabaci samples collected on sweet potato in Ghana clustered with reference sequences of samples from Cameroon corresponding to one of the five Sub-Saharan subgroups already described in the African clade. These data extends the phylogenetic information of the B. tabaci species complex and present new questions to be investigated.     

Genetic implications of phylogeographical patterns in the conservation of the boreal wetland butterfly Colias palaeno (Pieridae)

The boreo‐montane wetland butterfly species Colias palaeno has a European distribution from the Alps to northern Fennoscandia. Within its European range, the species’ populations have shrunk dramatically in recent historical times. Therefore, detailed baseline knowledge of the genetic makeup of the species is pivotal in planning potential conservation strategies. We collected 523 individuals from 21 populations across the entire European range and analyzed nuclear (20 allozyme loci) and mitochondrial (600 bp of the cytochrome c oxidase subunit I gene) genetic markers. The markers revealed contrasting levels of genetic diversity and divergence: higher in allozymes and lower in mitochondrial sequences. Five main groups were identified by allozymes: Alps, two Czech groups, Baltic countries, Fennoscandia, and Poland. The haplotype mitochondrial network indicates a recent range expansion. The most parsimonious interpretation for our results is the existence of a continuous Würm glacial distribution in Central Europe, with secondary disjunction during the Last Glacial Maximum into a south‐western and a north‐eastern fragment and subsequent moderate differentiation. Both groups present signs of postglacial intermixing in the Czech Republic. However, even a complete extinction in this region would not considerably affect the species’ genetic basis, as long as the source populations in the Alps and in northern Europe, comprising the most relevant evolutionary units for conservation, are surviving.     

Sugar maple (Acer saccharum Marsh.) establishment in boreal forest: results of a transplantation experiment

Aim To evaluate whether seedlings of sugar maple (Acer saccharum Marsh.) can establish beyond the species northern range limit in adjacent boreal forest. Location The hardwood–boreal forest transition zone on the north‐east shore of Lake Superior, Ontario, Canada. Methods Seed fall of sugar maple was monitored for 5 years in a stand of this species at its northern range limit, and seed from this stand was transplanted to five micro‐habitat types in an adjacent boreal forest. The establishment and survival of sugar maple seedlings there, and in the seed‐source stand, was monitored for the following 7–11 years. Soil‐surface light levels were measured in both forest types. Results Most seed fell in the final year of monitoring, when c. 250 seeds m^?2 were recorded. First‐year seedling establishment rates in the maple stand, deriving from this mast seed year, was approximately double that deriving from seed transplanted to the boreal forest sites; this is tentatively attributed to seed predator satiation in the maple stand. However, subsequent seedling survivorship in the boreal forest was greater than that in the maple stand, resulting in comparable seedling densities by the end of 6 years. This difference is tentatively attributed to better illumination in the boreal forest sites, and canopy‐opening disturbances appear to be especially facilitative of seedling survival. Main conclusions There is no fundamental impediment to sugar maple seedlings establishing in boreal forest communities if climate warming occurs and seed is available. If management intervention is needed to accelerate seed availability in a rapidly warming boreal forest, then diffuse seed application to disturbed boreal forest sites during mast years of local boreal tree species is recommended as the most effective way of avoiding seed predation and increasing seedling survival.     

Population genetic structure of the newly invasive Q biotype of Bemisia tabaci in Taiwan

The sweetpotato whitefly, Bemisia tabaci (Gennadius) (Hemiptera: Aleyrodidae), is among the top 100 invasive pests in the world, and it causes serious agricultural damage and economic losses in many countries. More than 24 biotypes of the sweetpotato whitefly have been detected worldwide, of which the Q biotype has recently been reported to be a new invasive pest spreading throughout the world via trade in poinsettias, Euphorbia pulcherrima Willd. ex Klotzsch (Euphorbiaceae). In 2006, the Q biotype was first recorded in Taiwan in greenhouses, but not in the field, suggesting that the invasion of this biotype might be at an early stage in that country. The mitochondrial cytochrome oxidase I (COI) gene and 12 microsatellite loci were used to investigate the genetic structure of multiple B. tabaci Q biotype populations. The presence of only a few COI haplotypes and a low number of nucleotide differences suggest high genetic similarity among these populations. Microsatellite analyses also revealed low genetic differentiation and frequent gene flow among greenhouses. The molecular evidence supports the occurrence of a recent genetic bottleneck in the B. tabaci Q biotype. Bayesian cluster analyses indicated that at least two invasion events have occurred in Taiwan. Phylogenetic analyses of microsatellites support Q biotype migration among greenhouses, which was likely facilitated by the frequent movement of poinsettias between greenhouses. Future management strategies should focus on developing plantlet trade regulations to avoid further anthropogenic dispersal of the B. tabaci Q biotype among greenhouses in Taiwan.     

Spatial expansion and population structure of the neotropical malaria vector, Anopheles darlingi (Diptera: Culicidae)

Extensive population structuring is known to occur in Anopheles darlingi , the primary malaria vector of the Neotropics. We analysed the phylogeographic structure of the species using the mitochondrial cytochrome oxidase I marker. Diversity is divided into six main population groups in South America: Colombia, central Amazonia, southern Brazil, south-eastern Brazil, and two groups in north-east Brazil. The ancestral distribution of the taxon is hypothesized to be central Amazonia, and there is evidence of expansion from this region during the late Pleistocene. The expansion was not a homogeneous front, however, with at least four subgroups being formed due to geographic barriers. As the species spread, populations became isolated from each other by the Amazon River and the coastal mountain ranges of south-eastern Brazil and the Andes. Analyses incorporating distances around these barriers suggest that the entire South American range of An. darlingi is at mutation–dispersal–drift equilibrium. Because the species is distributed throughout such a broad area, the limited dispersal across some landscape types promotes differentiation between otherwise proximate populations. Moreover, samples from the An. darlingi holotype location in Rio de Janeiro State are substantially derived from all other populations, implying that there may be additional genetic differences of epidemiological relevance. The results obtained contribute to our understanding of gene flow in this species and allow the formulation of human mosquito health protocols in light of the potential population differences in vector capacity or tolerance to control strategies.  © 2009 The Linnean Society of London, Biological Journal of the Linnean Society , 2009, 97 , 854–866.     

Genetic diversity of Bemisia tabaci (Gennadius) (Hemiptera: Aleyrodidae) populations and presence of the B biotype and a non-B biotype that can induce silverleaf symptoms in squash, in Uganda

The extent of genetic variability and host‐plant distribution of Bemisia tabaci (Gennadius) genotypes colonising cultivated and uncultivated plant species occurring adjacent to cassava fields in selected cassava‐producing areas of Uganda in 2003/04 were investigated using the mitochondrial cytochrome oxidase I (mtCOI) gene as the molecular marker. Eight genotype clusters, Ug1–Ug8, which are supported by high bootstrap values (≥80), at 3–18% nt divergence, were revealed among the collective Ugandan B. tabaci populations. Ug1 and Ug2 (both cassava‐associated) and Ug8 (sweetpotato‐associated) have been reported previously in Uganda. Ug3 was genetically dissimilar to B. tabaci described elsewhere and colonised a single species, Ocimum gratissimum. Ug4–Ug7 formed four closely related subclusters (93–97% nt identity) and diverged by 15–18% from Ug1, Ug2, Ug3 and Ug8, respectively. Ug4 had as its closest relatives (at 97–99% nt identity) the Ivory Coast okra biotype, whereas genotypes Ug5 and Ug6 had as their closest relatives (at 95–99% and 99% nt identity, respectively) the Mediterranean–North Africa–Middle East (MED‐NAFR‐ME) biotypes, which also include the well‐studied B and Q biotypes. Ug7 was closely related (at 98–99% nt identity) to biotype Ms from the Reunion Island in the Indian Ocean. Ug4 colonised Cucurbita pepo, Cucurbita sativus, Leonotis nepetifolia and Pavonia urens, while Ug7 colonised Commelina benghalensis, Gossypium hirsutum and Phaseolus vulgaris. Ug6, the B‐biotype‐like genotype colonised Abelmoschus esculentus and C. benghalensis only. None of Ug4–Ug7 genotypes was found associated with, or colonising, cassava or sweetpotato plants. In addition to colonising sweetpotato, the Ug8 genotypes colonised Lycopersicon esculentum and L. nepetifolia. Ug6 and Ug7, both members of the B biotype/B‐like cluster, induced silverleaf symptoms on Cucurbita sp. The discovery of five previously identified B. tabaci genotype clusters, Ug3–Ug7, in Uganda, among which are some of the world's most economically important biotypes, namely B and Q, is particularly significant in the spread of geminiviruses with devastating effects to crop production in Africa.     

Genetic structure of Japanese populations of an ambrosia beetle,Xylosandrus germanus (Curculionidae: Scolytinae)

We examined the genetic structures of 13 Japanese populations of an ambrosia beetle, Xylosandrus germanus (Curculionidae: Scolytinae), to understand the effects of geographical barriers on the colonization dynamics of this species. The genetic structure was studied using portions of the mitochondrial cytochrome oxidase I (COI) gene. A phylogenetic analysis revealed three distinct lineages (clades A, B and C) within X. germanus. Clade A contained 21 haplotypes from all 13 populations; whereas clade B contained eight haplotypes from Hokkaido (Sapporo and Furano), Iwate and Nagano populations; and clade C contained only a single a haplotype from the Hokkaido (Furano) population. In the analysis of molecular variance (amova ), the greatest amount of genetic variation was detected between populations in Hokkaido and those in Honshu and other southern islands. Between these two groups of populations, all the values of the coefficient of gene differentiation were significantly larger than zero, except for the Hokkaido (Sapporo) versus Nagano comparison. Our results confirm that for X. germanus, gene flow has been interrupted between Hokkaido and Honshu since the last glacial maximum.     

Aboveground impacts of a belowground invader: how invasive earthworms alter aboveground arthropod communities in a northern North American forest

Declining arthropod communities have recently gained a lot of attention, with climate and land-use change among the most frequently discussed drivers. Here, we focus on a seemingly underrepresented driver of arthropod community decline: biological invasions. For approximately 12 000 years, earthworms have been absent from wide parts of northern North America, but they have been re-introduced with dramatic consequences. Most studies investigating earthworm-invasion impacts focus on the belowground world, resulting in limited knowledge on aboveground-community changes. We present observational data on earthworm, plant and aboveground arthropod communities in 60 plots, distributed across areas with increasing invasion status (low, medium and high) in a Canadian forest. We analysed how earthworm-invasion status and biomass impact aboveground arthropod community abundance, biomass and species richness, and how earthworm impacts cascade across trophic levels. We sampled approximately 13 000 arthropods, dominated by Hemiptera, Diptera, Araneae, Thysanoptera and Hymenoptera. Total arthropod abundance, biomass and species richness declined significantly from areas of low to those with high invasion status, with reductions of 61, 27 and 18%, respectively. Structural equation models suggest that earthworms directly and indirectly impact arthropods across trophic levels. We show that earthworm invasion can alter aboveground multi-trophic arthropod communities and suggest that belowground invasions might be underappreciated drivers of aboveground arthropod decline.     

Genetic diversity of genus Vespa including an invaded species of V. velutina (Hymenoptera: Vespidae) in Korea inferred from DNA barcoding data

With about 5000 known species, the Vespidae is a large family belongs to order Hymenoptera. The genus Vespa with 22 species is one of the four genera of the subfamily Vespinae. In Korea, 10 species and subspecies are recognized. Because of their social behavior, their treat to human health and their impact in apiculture, the reliable and sometimes automated identification of these insects to species level are important. To test the efficacy of DNA barcoding method for identification of species of the genus Vespa in Korea, 30 samples of eight Korean species of genus Vespa were collected and mitochondrial DNAs of 658 bp fragment cytochrome oxidase subunit 1 (CO1) region were sequenced. A Bayesian Inference based on COI gene of the Korean Vespa species was constructed. The phylogenetic tree shoed that identification of all specimens is possible based on COI gene and we found strong relation between the sequences of the collected species from different localities in South Korea which clustered together with 100% support with sequences of the same species in GenBank. The results demonstrate that DNA barcoding is a useful technique for rapid and accurate species recognition in Korean Vespa species. The DNA barcode part of COI for V. binghami is provided for the first time that can help for identification of this species through DNA barcoding. Also, the genetic diversity among Korean Vespa velutina was zero suggests that the invasion might have occurred in a single event with small number of founders.     

Contrasting patterns of population structure in the montane caddisflies Hydropsyche tenuis and Drusus discolor in the Central European highlands

1. In this study, we compared mitochondrial sequence data (cytochrome oxidase I) to infer the population structure of the two montane caddisflies Hydropsyche tenuis and Drusus discolor. The two species are contrasting examples of montane aquatic insects with insular distributions: D. discolor is restricted to altitudes above 600 m, H. tenuis is limited to the same mountain ranges in Central Europe but inhabits lower altitudes. 2. In particular, we ask whether these two species with similar regional distributions show similar patterns of population structure and haplotype diversity, and whether any differences can be attributed to population history and/or autecology. 3. To determine the population structure of both species, we applied conventional population genetics analyses to mitochondrial sequence data. We collected and sampled 121 specimens of H. tenuis from 29 sites in 10 different regions of the Central European highlands and 138 individuals of D. discolor from 40 sites in 11 different regions. 4. Nine unique haplotypes were identified for H. tenuis and 34 for D. discolor. There were eight variable positions in H. tenuis and 41 in D. discolor. The maximum difference between haplotypes was 0.8% (4 bp) for H. tenuis and 4.2% (21 bp) for D. discolor. We observed haplotype overlap between geographic regions for both species. Analysis of molecular variance showed that two‐thirds of the total variance in H. tenuis was found among regions while in D. discolor, a larger portion of variance was found within regions and populations due to a higher number of haplotypes observed within regions. Mantel test showed a significant relationship between genetic and geographic distance in D. discolor, but no significant relationship in H. tenuis. 5. Our analyses show that, despite their very similar overall distribution pattern in Europe, the two species exhibit distinct population structures, which may reflect differences in phylogeographic history, dispersal capabilities, habitat specifity or within‐region geographic occurrence.     

From introduction to equilibrium: reconstructing the invasive pathways of the Argentine ant in a Mediterranean region

Determining the geographical range of invasive species is an important component of formulating effective management strategies. In the absence of detailed distributional data, species distribution models can provide estimates of an invasion range and increase our understanding of the ecological processes acting at various spatial scales. We used two complementary approaches to evaluate the influence of historical and environmental factors in shaping the distribution of the Argentine ant ( Linepithema humile ), a widespread, highly invasive species native to South America. Occurrence data were combined with environmental data at incremental spatial scales (extent and resolution) to predict the suitable range of the ant invasion using ecological niche models. In addition, we also used a spread model that simulated the jump dispersal of the species to identify the most plausible scenarios of arrival of L. humile in the NE Iberian Peninsula at local scales. Based on the results of both modelling practices, we suggest that L. humile might have reached its maximum geographic range at regional scales in the NE Iberian Peninsula. However, the species does not appear in equilibrium with the environment at small spatial scales, and further expansions are expected along coastal and inland localities of the Costa Brava. Long-distance jumps are ultimately responsible for the spread of the Argentine ant in the area. Overall, our study shows the utility of combining niche based models with spread models to understand the dynamics of species' invasions.     

Genetic diversity and population structure of the commercially harvested sea urchin Paracentrotus lividus (Echinodermata, Echinoidea)

The population structure of the edible Atlanto-Mediterranean sea urchin Paracentrotus lividus is described by analysing sequence variation in a fragment of the mitochondrial gene cytochrome c oxidase subunit I in 127 individuals from 12 localities across south-west Europe. The study revealed high levels of genetic diversity but low levels of genetic structure, suggesting a large degree of gene flow between populations and panmixis within each, the Mediterranean and Atlantic basins. However, we found significant genetic differentiation between the two basins probably due to restricted gene flow across the geographical boundary imposed by the area of the Strait of Gibraltar. Populations of P. lividus appeared to have experienced a recent demographic expansion in the late Pleistocene. We provide new evidence on the population structure of this commercial species, predicting a healthy stock of this sea urchin on the Mediterranean and Atlantic coasts.     

Spatial patterns of genetic structure among populations of a stone‐cased caddis (Trichoptera: Tasimiidae) in south‐east Queensland,Australia

1. Gene flow and dispersal among populations of a stone‐cased caddis (Tasimiidae: Tasimia palpata) were estimated indirectly using a 460 bp region of the cytochrome oxidase I gene of mitochondrial DNA. 2. There was no significant differentiation at the largest spatial scale (between catchments) and no correlation between genetic distance and geographic distance. These results are consistent with widespread adult dispersal. 3. Conversely, significant genetic differentiation was detected at the smallest spatial scale examined (among reaches within streams). This pattern was primarily because of significant F_ST values in a single stream (Bundaroo Creek). 4. Bundaroo Creek also had the lowest mean number of haplotypes per population (n = 7) suggesting that a limited number of females may be responsible for recruitment at these sites. Significant F_ST's at the reach scale may be a result of this ‘patchy’ recruitment. However, additional evidence regarding the long‐range dispersal ability and fecundity of T. palpata females is needed to test this hypothesis fully.     

Molecular systematics of the Chrysoperla carnea group (Neuroptera: Chrysopidae) in Europe

The green lacewing Chrysoperla carnea is a complex of cryptic species whose differentiation has been so far based upon morphology, ecophysiology, behaviour and preliminary mitochondrial DNA data using cytochrome oxidase subunit II (COII) and NADH dehydrogenase subunit II. In this work we extended the DNA data by screening nucleotide sequences of COII, cytochrome oxidase I, cytochrome b and the large ribosomal subunit of the mtDNA. These new data suggest that C. carnea s.s. is a well-supported, separate taxon, but that other taxa of the complex are not consistently differentiated by the current DNA data.     

Species status and population structure of the Australian Eucalyptus pest Paropsis atomaria Olivier (Coleoptera: Chrysomelidae)

1 Paropsis atomaria Olivier represents an emergent pest of Eucalyptus plantations in Queensland and New South Wales, Australia. Most prior studies on the biology and control of P. atomaria have centred on populations from Canberra in the Australian Capital Territory, but the biological relationship between beetles from Canberra and those from up to 1500 km further north are unknown. 2 DNA markers were used to determine whether P. atomaria from Canberra are the same biological species as those from Eucalyptus forestry plantations in northern New South Wales and Queensland, where the beetle has become an important pest. Using the mitochondrial gene, cytochrome c oxidase I (COI), individuals collected from across the distribution of P. atomaria were investigated for haplotype diversity and levels of mitochondrial divergence. 3 Within P. atomaria, genetic distance averaged 0.5% across 23 unique haplotypes for 93 individuals, with an average of 14% difference between P. atomaria and the outgroup species, Paropsis obsoleta. Significant genetic structure was observed relative to geographical distribution, but not with respect to host plant species of origin. Greatest divergence was between the southern‐most sample site (Canberra) and northern sites in New South Wales and Queensland, indicating reduced gene flow between these regions. 4 Individuals from across eastern Australia belong to the same genetic species with population substructuring evident. Consequently, there is no evidence to suggest cryptic species complexes exist within the currently defined taxon. Continued implementation of control strategies for P. atomaria across its distribution is appropriate.     

Crystal structure of the Mycobacterium tuberculosis dUTPase: insights into the catalytic mechanism

The structure of Mycobacterium tuberculosis dUTP nucleotidohydrolase (dUTPase) has been determined at 1.3 Angstrom resolution in complex with magnesium ion and the non-hydrolyzable substrate analog, alpha,beta-imido dUTP. dUTPase is an enzyme essential for depleting potentially toxic concentrations of dUTP in the cell. Given the importance of its biological role, it has been proposed that inhibiting M.tuberculosis dUTPase might be an effective means to treat tuberculosis infection in humans. The crystal structure presented here offers some insight into the potential for designing a specific inhibitor of the M.tuberculosis dUTPase enzyme. The structure also offers new insights into the mechanism of dUTP hydrolysis by providing an accurate representation of the enzyme-substrate complex in which both the metal ion and dUTP analog are included. The structure suggests that inclusion of a magnesium ion is important for stabilizing the position of the alpha-phosphorus for an in-line nucleophilic attack. In the absence of magnesium, the alpha-phosphate of dUTP can have either of the two positions which differ by 4.5 Angstrom. A transiently ordered C-terminal loop further assists catalysis by shielding the general base, Asp83, from solvent thus elevating its pK(a) so that it might in turn activate a tightly bound water molecule for nucleophilic attack. The metal ion coordinates alpha, beta, and gamma phosphate groups with tridentate geometry identical with that observed in the crystal structure of DNA polymerase beta complexed with magnesium and dNTP analog, revealing some common features in catalytic mechanism.