Population genetics of Scirtothrips perseae: tracing the origin of a recently introduced exotic pest of Californian avocado orchards, using mitochondrial and microsatellite DNA markers

The Californian avocado industry has recently been impacted by the establishment of three exotic arthropod pests that are native to Mexico and Central America. Establishment of all three pests is thought to have resulted from illegal movement of host plants (avocado) into California. To identify likely sources and routes of entry of such pest invasions, we examined the population genetics of the most recent invader, avocado thrips [Scirtothrips perseae Nakahara (Thysanoptera: Thripidae)], using mitochondrial DNA (mtDNA) and microsatellite markers. The mtDNA sequences revealed three geographically distinct and divergent lineages, of which the mtDNA haplotypes of Californian individuals were most closely related to populations in the center of the pest's native range. Analysis of allele frequencies at four microsatellite loci added resolution, indicating Coatepec Harinas, Mexico, as the most likely source of the Californian population. Statistically, we did not detect any bottleneck in population size associated with the invasion of California. However, estimates of the effective population size of the invading population suggest that a severe bottleneck occurred, indicating that the quantity of host‐plant material entering California was small. Our findings implicate Coatepec Harinas, a large avocado germplasm and breeding center, as the most likely source of the introduced Californian population of S. perseae and as a likely source of previous and future avocado pest introductions. Efforts to identify natural enemies of S. perseae for biological control should focus on Coatepec Harinas and immediate surrounds. Moreover, identification of the source of invasive pests enables the establishment and enforcement of plant quarantine and free‐trade protocols.     

Intron size polymorphism of the Adh1 gene parallels the worldwide colonization history of the Mediterranean fruit fly, Ceratitis capitata

The exon-primed intron-crossing (EPIC) PCR technique was used to analyse the size variation at the first intron of the Ceratitis capitata Adh ₁ gene. A total of 27 samples from 16 natural populations was analysed from five geographical regions in the species range: Africa, Mediterranean Basin, Latin America, Hawaii and Australia. The Adh ₁ first intron varies extensively in length with at least 18 size variants ranging from 1400 bp to 3450 bp. These variants can be grouped into four distinct size categories: short, medium, long and very long. The majority of these variants are present only in the African populations. Only a subset of the ancestral variants appear to have succeeded in migrating from Africa during the medfly colonization process. The medfly population structure inferred from the intron size polymorphism is congruent with that observed from the analysis of allozyme variation. The geographical dispersal of the medfly from its source area is associated with a gradual and great reduction in intron variability which parallels the trend of decreasing variability evaluated at 26 biochemical loci. The intron phylogenetic tree is in agreement with allozyme data in portraying the dynamic population history of the medfly. Stochastic evolutionary forces such as drift, bottleneck effects and migration seem to have played the major roles in the dispersion pattern of Adh ₁ intron variation during the colonization of the medfly.     

Mitochondrial DNA variation in natural populations of the mushroom Agaricus bisporus

We investigated the patterns of mitochondrial DNA variation in the global population of the commercial mushroom Agaricus bisporus . Through the analysis of RFLP's among 441 isolates from nine countries in North America and Eurasia, we found a total of 140 mtDNA haplotypes. Based on population genetic analysis, there are four genetically distinct natural populations in this species, found in coastal California, desert California, France and Alberta (Canada). While 134 of the 140 mtDNA haplotypes were unique to single geographical regions, two mtDNA haplotypes, mt001 and mt002, were found in almost every population surveyed. These two mtDNA haplotypes also predominate among cultivars used throughout the world for at least the last two decades. These two mtDNA haplotypes are more similar to the cosmopolitan groups of mtDNA haplotypes than to the indigeneous clusters of mtDNA haplotypes from the two Californian regions.     

Female frequencies and fitness components between sex phenotypes among gynodioecious populations of the colonizing species Trifolium hirtum All. in California

Summary Male sterility has been recently discovered in Californian populations of rose clover (Trifolium hirtum). This study describes the frequency distribution of male sterility in Turkish and Californian populations, and compares fitness components of hermaphrodites and females. As male-steriles were found in Turkey, it is likely that they were introduced to California during the 1940's with the original material derived from Turkey. The spread of male-sterile genotypes in California has given rise to an asymmetrical frequency distribution of male sterility with positive skewness. The frequency of females has not exceeded fifty percent and it appears to be temporally stable in most of the Californian populations. The hypothesis that female frequencies and fitness differences between phenotypes are correlated was tested by comparing sex phenotypes in seven populations with contrasting levels of male sterility. The analysis of those populations showed no evidence for such a correlation as no significant differences were found between sex phenotypes for fecundity and seed germination. The hypothesis that females are maintained due to fitness differences in the progeny of hermaphrodites and females was experimentally tested in the population with maximum frequency of male-steriles. The results showed no significant differences in the demographic performance of the progenies of hermaphrodites and females. The present results are discussed in terms of the possible mechanism of maintenance of gynodioecy in rose clover.     

Population-specific demography and invasion potential in medfly

Biological invasions are constantly gaining recognition as a significant component of global change. The Mediterranean fruit fly (medfly) constitutes an ideal model species for the study of biological invasions due to its (1) almost cosmopolitan geographic distribution, (2) huge economic importance, and (3) well-documented invasion history. Under a common garden experimental set up, we tested the hypothesis that medfly populations obtained from six global regions [Africa (Kenya), Pacific (Hawaii), Central America (Guatemala), South America (Brazil), Extra-Mediterranean (Portugal), and Mediterranean (Greece)] have diverged in important immature life-history traits such as preadult survival and developmental times. We also tested the hypothesis that medfly populations from the above regions exhibit different population growth rates. For this purpose, data on the life history of immatures were combined with adult survival and reproduction data derived from an earlier study in order to calculate population parameters for the above six populations. Our results clearly show that medfly populations worldwide exhibit significant differences in preadult survival, developmental rates of immatures and important population parameters such as the intrinsic rate of increase. Therefore, geographically isolated medfly populations may share different invasion potential, since population growth rates could influence basic population processes that operate mostly during the last two stages of an invasion event, such as establishment and spread. Our findings provide valuable information for designing population suppression measures and managing invasiveness of medfly populations worldwide.     

Inferring Geographic Isolation of Wolverines in California Using Historical DNA

ABSTRACT Delineating a species' geographic range using the spatial distribution of museum specimens or even contemporary detection-non-detection data can be difficult. This is particularly true at the periphery of a species range where species' distributions are often disjunct. Wolverines (Gulo gulo) are wide-ranging mammals with discontinuous and potentially isolated populations at the periphery of their range. One potentially disjunct population occurred in the Sierra Nevada Mountains, California, USA, and appears to have been extirpated by the 1930s. Many early 20th century naturalists believed that this population was connected to other populations occurring in the Cascade Range of northern California, Oregon, and Washington, USA, but a recent analysis of historical records suggests that California wolverines were isolated from other populations in North America. We used DNA extracted from museum specimens to examine whether California wolverines were isolated. Both nuclear and mitochondrial DNA data indicate that California wolverines were genetically distinct from extant populations, suggesting long-term isolation. We identified 2 new control region (mitochondrial DNA) haplotypes located only within California. We used these data and referenced sequences from the Rocky Mountains, USA, to make inferences regarding potential wolverine translocations into California. In addition, we used these genetic data to make inferences about wolverine conservation throughout western North America.     

Population structuring of the ubiquitous stingless bee Tetragonisca angustula in southern Brazil as revealed by microsatellite and mitochondrial markers

Tetragonisca angustula is one of the most widespread stingless bees in the Neotropics. This species swarms frequently and is extremely successful in urban environments. In addition, it is one of the most popular stingless bee species for beekeeping in Latin America, so nest transportation and trading is common. Nest transportation can change the genetic structure of the host population, reducing inbreeding and increasing homogenization. Here, we evaluate the genetic structure of 17 geographic populations of T. angustula in southern Brazil to quantify the level of genetic differentiation between populations. Analyses were conducted on partially sequenced mitochondrial genes and 11 microsatellite loci of 1002 workers from 457 sites distributed on the mainland and on 3 islands. Our results show that T. angustula populations are highly differentiated as demonstrated by mitochondrial DNA (mtDNA) and microsatellite markers. Of 73 haplotypes, 67 were population‐specific. MtDNA diversity was low in 9 populations but microsatellite diversity was moderate to high in all populations. Microsatellite data suggest 10 genetic clusters and low level of gene flow throughout the studied area. However, physical barriers, such as rivers and mountain ranges, or the presence or absence of forest appear to be unrelated to population clusters. Factors such as low dispersal, different ecological conditions, and isolation by distance are most likely shaping the population structure of this species. Thus far, nest transportation has not influenced the general population structure in the studied area. However, due to the genetic structure we found, we recommend that nest transportation should only occur within and between populations that are genetically similar.     

Establishment failure in biological invasions: a case history of Littorina littorea in California, USA

Background

The early stages of biological invasions are rarely observed, but can provide significant insight into the invasion process as well as the influence vectors have on invasion success or failure.

Methodology/Principal Findings

We characterized three newly discovered populations of an introduced gastropod, Littorina littorea (Linné, 1758), in California, USA, comparing them to potential source populations in native Europe and the North American East Coast, where the snail is also introduced. Demographic surveys were used to assess spatial distribution and sizes of the snail in San Francisco and Anaheim Bays, California. Mitochondrial DNA was sequenced and compared among these nascent populations, and various populations from the North American East Coast and Europe, to characterize the California populations and ascertain their likely source. Demographic and genetic data were considered together to deduce likely vectors for the California populations. We found that the three large California L. littorea populations contained only adult snails and had unexpectedly high genetic diversity rather than showing an extreme bottleneck as typically expected in recent introductions. Haplotype diversity in Californian populations was significantly reduced compared to European populations, but not compared to East Coast populations. Genetic analyses clearly suggested the East Coast as the source region for the California introductions.

Conclusions and Significance

The California L. littorea populations were at an early, non-established phase of invasion with no evidence of recruitment. The live seafood trade is the most likely invasion vector for these populations, as it preferentially transports large numbers of adult L. littorea, matching the demographic structure of the introduced California L. littorea populations. Our results highlight continued operation of live seafood trade vectors and the influence of vectors on the demographic and genetic structure of the resulting populations, especially early stages of the invasion process.     

Invasion Genetics of New World Medflies: Testing Alternative Colonization Scenarios

The Mediterranean fruit fly (Ceratitis capitata) is an invasive agricultural pest with a wide host range and a nearly global distribution. Efforts to forgo the medfly's spread into the United States are dependent on an understanding of population dynamics in newly established populations elsewhere. To explore the potential influence of demographic and historical parameters in six medfly populations distributed from Mexico to Peru, we created population genetic null models using Monte Carlo simulations. Null expectations for genetic differentiation (F _ST) were compared with actual sequence variation from four highly polymorphic nuclear loci. Four colonization scenarios that were modeled led to unique genetic signatures that could be used to interpret empirical data. Unless current gene flow across Latin America was assumed to be very high, we could reject colonizations consisting of multiple introductions, each of low genetic diversity. Further, if simulated populations were small (N _e = 5 × 10² individuals per population), small invasions from a single source consistently produced F _ST values comparable to those currently observed in Latin America. In contrast, only large invasions from diverse sources were compatible with the observed data for large populations (N _e 5 × 10³). This study demonstrates that alternative population genetic hypotheses can be tested empirically even when departures from equilibrium are extreme, and that population genetic theory can be used to explore the processes that underlie biological invasions.     

Multilocus genealogies reveal multiple cryptic species and biogeographical complexity in the California turret spider Antrodiaetus riversi (Mygalomorphae, Antrodiaetidae)

Antrodiaetus riversi (Araneae, Antrodiaetidae) is a dispersal-limited, habitat specialized mygalomorph spider species endemic to mesic woodlands of northern and central California. This species occupies a disjunct distribution, with populations in the Sierra Nevada and Coast Ranges, separated by the inhospitable Central Valley. Previous studies of morphological and allozyme variation have suggested that these populations may constitute cryptic species. We investigated the phylogeography of A. riversi using both nuclear and mitochondrial DNA sequences, collected for a comprehensive population sample. These data reveal the presence of at least five species in the A. riversi complex - these species are deeply diverged, and genealogically exclusive in both nuclear and mitochondrial genomes. Each of these species is characterized by extreme population subdivision and deep phylogeographical structuring, consistent with minimal gene flow across the dissected Californian landscape. Three species are restricted to the Coast Ranges, one to high altitudes of the central Sierran Nevada, and one species is found in both ranges. These species have allopatric distributions, although species parapatry is hypothesized to occur in several areas. Species diversification appears to have pulsed in the Late Miocene/Early Pliocene, a timing consistent with biogeographical reconstructions for many Californian taxa, and a time of turbulent geological activity in the region.     

Contemporary patterns in a historical context: phylogeographic history of the pipevine swallowtail,Battus philenor (Papilionidae)

Abstract We examined mitochondrial DNA (mtDNA) variation in pipevine swallowtail butterflies ( Battus philenor ) from throughout its extant range to provide a historical, phylogeographical context for ecological studies of the disjunct population in California. We evaluate current hypotheses regarding host plant use, behavior, and mimetic relationships of B. philenor populations and generate alternative hypotheses. Compared to populations throughout the rest of the species' range, California populations are ecologically distinct in that they lack mimics, lay significantly larger clutches of eggs, and exclusively use a unique, endemic larval host plant. Analysis of molecular variance, tests of population differentiation, and nested clade analysis of mtDNA variation indicate that, despite low levels of population genetic structure across the species' range, there is evidence of recent range expansion from presumed Pleistocene refuge(s) in southeastern North America. Colonization of California appears to have been a recent event. This phylogeographic investigation also suggests that the evolution of life-history adaptations to a novel larval host has occurred rapidly in California and the lack of mimics in California may be attributable to the recency of colonization.     

Geographic patterns of genome admixture in Latin American Mestizos

The large and diverse population of Latin America is potentially a powerful resource for elucidating the genetic basis of complex traits through admixture mapping. However, no genome-wide characterization of admixture across Latin America has yet been attempted. Here, we report an analysis of admixture in thirteen Mestizo populations (i.e. in regions of mainly European and Native settlement) from seven countries in Latin America based on data for 678 autosomal and 29 X-chromosome microsatellites. We found extensive variation in Native American and European ancestry (and generally low levels of African ancestry) among populations and individuals, and evidence that admixture across Latin America has often involved predominantly European men and both Native and African women. An admixture analysis allowing for Native American population subdivision revealed a differentiation of the Native American ancestry amongst Mestizos. This observation is consistent with the genetic structure of pre-Columbian populations and with admixture having involved Natives from the area where the Mestizo examined are located. Our findings agree with available information on the demographic history of Latin America and have a number of implications for the design of association studies in population from the region.     

Population Genetics of Ceratitis capitata in South Africa: Implications for Dispersal and Pest Management

The invasive Mediterranean fruit fly (medfly), Ceratitis capitata, is one of the major agricultural and economical pests globally. Understanding invasion risk and mitigation of medfly in agricultural landscapes requires knowledge of its population structure and dispersal patterns. Here, estimates of dispersal ability are provided in medfly from South Africa at three spatial scales using molecular approaches. Individuals were genotyped at 11 polymorphic microsatellite loci and a subset of individuals were also sequenced for the mitochondrial cytochrome oxidase subunit I gene. Our results show that South African medfly populations are generally characterized by high levels of genetic diversity and limited population differentiation at all spatial scales. This suggests high levels of gene flow among sampling locations. However, natural dispersal in C. capitata has been shown to rarely exceed 10 km. Therefore, documented levels of high gene flow in the present study, even between distant populations (>1600 km), are likely the result of human-mediated dispersal or at least some form of long-distance jump dispersal. These findings may have broad applicability to other global fruit production areas and have significant implications for ongoing pest management practices, such as the sterile insect technique.     

Admixture studies in Latin America: from the 20th to the 21st century

The present study is a review of admixture studies in Latin America, an interesting subject because of the unique history of the area, in which populations from 3 different origins had contact and intercrossed. The most often used methods of analysis of admixture in Latin America and some problems related to them, such as the determination of the parental populations and selection of genetic markers, are briefly reviewed. Several sources of data for admixture studies (surnames, quantitative traits, proteins, and molecular information) are summarized. The results obtained using protein systems and blood groups, the most often used markers in Latin America, are considered. They are classified according to their application in 3 groups of populations: urban centers, native Americans, and African-descended subjects. The data show that almost every population is dihybrid or trihybrid, and when African influence is not detected, it is probably due more to the method than to an absence of that contribution. A special section is dedicated to the direction of gene flow, and results about directional mating based on mtDNA, Y-chromosome, and nuclear DNA or proteins are also given. From these studies it is possible to conclude that Amerindian admixture came mainly from female lineages, but it is difficult to establish what happened with the African contribution. A last subject considered is the relation between interethnic crosses and diseases; it is easy to analyze that relation when the pathological condition is related to a unique allele, but when complex diseases are considered, the results are not as clear because of the influence of nongenetic factors. Finally, the perspectives for admixture studies in the 21st century are considered, and some attempts to predict their future in Latin America are made.     

Sunshine versus gold: The effect of population age on genetic structure of an invasive mosquito

The genetic diversity and structure of invasive species are affected by the time since invasion, but it is not well understood how. We compare likely the oldest populations of Aedes aegypti in continental North America with some of the newest to illuminate the range of genetic diversity and structure that can be found within the invasive range of this important disease vector. Aedes aegypti populations in Florida have probably persisted since the 1600‐1700s, while populations in southern California derive from new invasions that occurred in the last 10 years. For this comparison, we genotyped 1,193 individuals from 28 sites at 12 highly variable microsatellites and a subset of these individuals at 23,961 single nucleotide polymorphisms (SNPs). This is the largest sample analyzed for genetic structure for either region, and it doubles the number of southern California populations previously analyzed. As predicted, the older populations (Florida) showed fewer indicators of recent founder effect and bottlenecks; in particular, these populations have dramatically higher genetic diversity and lower genetic structure. Geographic distance and driving distance were not good predictors of genetic distance in either region, especially southern California. Additionally, southern California had higher levels of genetic differentiation than any comparably sized documented region throughout the worldwide distribution of the species. Although population age and demographic history are likely driving these differences, differences in climate and transportation practices could also play a role.     

Population genetics of black abalone, Haliotis cracherodii, along the central California coast

Over the past three decades, the black abalone, Haliotis cracherodii, has experienced precipitous declines in abundance over portions of its range in southern and central California. The potential for recovery of these populations is dependent in part on dispersal processes; that is, can distant populations serve as sources of recruits to locales that no longer harbor H. cracherodii? Here we use population genetic analysis to assess levels of population subdivision and infer recruitment processes. Epipodial tissue samples were obtained from over 400 black abalone from seven geographic sites between Santa Cruz and Santa Barbara counties in central California. Allelic frequencies were determined for each population at three polymorphic enzyme-encoding loci (GPI, AAT-1 and PGM). Significant allelic frequency differentiation among sites was observed at all three loci. Genetic distance was found to be independent of geographic distance over the approximately 300-km sampling range. In addition, a limited number of DNA sequences (total N=51) were obtained for the mitochondrial cytochrome oxidase subunit I gene (COI) from five of the populations. Since the same common COI haplotype dominated each population, this analysis had little statistical power and failed to detect population structure. The observed level of population differentiation at allozyme loci was three-fold higher than that observed in California red abalone, H. rufescens. The species differ in their breeding period and it is suggested that the relatively short, summer breeding season of black abalone limits dispersal because larvae experience reduced variance in oceanographic conditions relative to red abalone that spawn year-round. Based on these results, rates of recolonization and recovery of locally depressed or extirpated black abalone populations are likely to be slow despite harvest restrictions.     

Patterns of introduction and adaptation during the invasion of Aegilops triuncialis (Poaceae) into Californian serpentine soils

Multiple introductions can play a prominent role in explaining the success of biological invasions. One often cited mechanism is that multiple introductions of invasive species prevent genetic bottlenecks by parallel introductions of several distinct genotypes that, in turn, provide heritable variation necessary for local adaptation. Here, we show that the invasion of Aegilops triuncialis into California, USA, involved multiple introductions that may have facilitated invasion into serpentine habitats. Using microsatellite markers, we compared the polymorphism and genetic structure of populations of Ae. triuncialis invading serpentine soils in California to that of accessions from its native range. In a glasshouse study, we also compared phenotypic variation in phenological and fitness traits between invasive and native populations grown on loam soil and under serpentine edaphic conditions. Molecular analysis of invasive populations revealed that Californian populations cluster into three independent introductions (i.e. invasive lineages). Our glasshouse common garden experiment found that all Californian populations exhibited higher fitness under serpentine conditions. However, the three invasive lineages appear to represent independent pathways of adaptation to serpentine soil. Our results suggest that the rapid invasion of serpentine habitats in California may have been facilitated by the existence of colonizing Eurasian genotypes pre‐adapted to serpentine soils.     

The ruddy duck Oxyura jamaicensis in Europe: natural colonization or human introduction?

Native to North America, ruddy ducks Oxyura jamaicensis now occur in 21 countries in the western Palaearctic (including Iceland) and their expanding population threatens the native white-headed duck, Oxyura leucocephala, through hybridization and possibly competition for food and nest sites. We used mitochondrial DNA sequences and nuclear microsatellites to test whether the European ruddy duck population is descended solely from the captive population in the UK, which traces to seven individuals imported from the USA in 1948, or, alternatively, has been augmented by natural dispersal of birds from North America. Limited genetic diversity in the European population is consistent with a founder population as small as seven birds. In addition, shifts in allele frequencies at several loci, presumably due to genetic drift in the founding population, result in significant differentiation between the European and North American populations. Despite the recent separation of these populations, almost all individuals could be unambiguously assigned based on their composite genotypes, to one of two distinct populations, one comprising all of the European ruddy ducks we sampled (including those from Iceland and captive birds in the UK) and the other comprising all North American samples. Our results confirm that the European ruddy duck population is likely to derive solely from the captive population in the UK and we find no evidence of recent arrivals from North America or of admixture between ruddy ducks from Europe and North America.     

Why do Californian striders fly?

Discerning the adaptive significance of migratory strategies poses significant challenges, not the least of which is measuring migratory capability in natural populations. We take advantage of a visible migratory dimorphism to study variation in migratory capability in the stream-dwelling water strider, Aquarius remigis. Theory predicts loss of migratory capability in this species because streams have been viewed as stable and persistent habitats. As expected, A. remigis lack wings throughout most of North America. However, Californian populations are noted for unexpectedly high frequencies of winged, migratory morphs. To deduce the adaptive significance of this anomalous regional variation, we compare proportion winged among 37 Californian populations. We discover a strong, positive correlation with altitude, but no correlations with latitude, rainfall or stream size. A common garden experiment reveals that both proportion winged and its reaction norm to temperature differ genetically among populations, and a half-sibling experiment demonstrates that wing morph has high heritability, moderate genetic correlations across environments and a significant genotype by environment interaction. These results support the hypothesis that proportion winged and its reaction norm to temperature have diverged genetically in California. We conclude that high migratory capability is an evolutionary adaptation to the unusual harshness and instability of Californian stream habitats, and particularly to the high elevational gradients and extreme seasonal variation characteristic of montane streams.     

Patterns of Admixture and Population Structure in Native Populations of Northwest North America

The initial contact of European populations with indigenous populations of the Americas produced diverse admixture processes across North, Central, and South America. Recent studies have examined the genetic structure of indigenous populations of Latin America and the Caribbean and their admixed descendants, reporting on the genomic impact of the history of admixture with colonizing populations of European and African ancestry. However, relatively little genomic research has been conducted on admixture in indigenous North American populations. In this study, we analyze genomic data at 475,109 single-nucleotide polymorphisms sampled in indigenous peoples of the Pacific Northwest in British Columbia and Southeast Alaska, populations with a well-documented history of contact with European and Asian traders, fishermen, and contract laborers. We find that the indigenous populations of the Pacific Northwest have higher gene diversity than Latin American indigenous populations. Among the Pacific Northwest populations, interior groups provide more evidence for East Asian admixture, whereas coastal groups have higher levels of European admixture. In contrast with many Latin American indigenous populations, the variance of admixture is high in each of the Pacific Northwest indigenous populations, as expected for recent and ongoing admixture processes. The results reveal some similarities but notable differences between admixture patterns in the Pacific Northwest and those in Latin America, contributing to a more detailed understanding of the genomic consequences of European colonization events throughout the Americas.