|Genetic evidence links invasive monk parakeet populations in the United States to the international pet trade.|
|Jump to Full Text|
|PMID: 18652686 Owner: NLM Status: MEDLINE|
|BACKGROUND: Severe ecological and economic impacts caused by some invasive species make it imperative to understand the attributes that permit them to spread. A notorious crop pest across its native range in South America, the monk parakeet (Myiopsitta monachus) has become established on four other continents, including growing populations in the United States. As a critical first step to studying mechanisms of invasion success in this species, here we elucidated the geographical and taxonomic history of the North American invasions of the monk parakeet. Specifically, we conducted a genetic assessment of current monk parakeet taxonomy based on mitochondrial DNA control region sequences from 73 museum specimens. These data supported comparative analyses of mtDNA lineage diversity in the native and naturalized ranges of the monk parakeet and allowed for identification of putative source populations.
RESULTS: There was no molecular character support for the M. m. calita, M. m. cotorra, and M. m. monachus subspecies, while the Bolivian M. m. luchsi was monophyletic and diagnosably distinct. Three haplotypes sampled in the native range were detected within invasive populations in Florida, Connecticut, New Jersey and Rhode Island, the two most common of which were unique to M. m. monachus samples from eastern Argentina and bordering areas in Brazil and Uruguay.
CONCLUSION: The lack of discrete morphological character differences in tandem with the results presented here suggest that M. m. calita, M. m. cotorra and M. m. monachus are in need of formal taxonomic revision. The genetic distinctiveness of M. m. luchsi is consistent with previous recommendations of allospecies status for this taxon. The geographic origins of haplotypes sampled in the four U.S. populations are concordant with trapping records from the mid-20th century and suggest that propagule pressure exerted by the international pet bird trade contributed to the establishment of invasive populations in the United States.
|Michael A Russello; Michael L Avery; Timothy F Wright|
Related Documents :
|24203186 - Genetic characterization and relationships of populus alba, p. tremula, and p. x canesc...
24976076 - The evolutionary history of darwin's finches: speciation, gene flow, and introgression ...
23379646 - No evidence for nuclear introgression despite complete mtdna replacement in the carpath...
24085436 - Short communication microsatellite markers in paulownia kawakamii (scrophulariaceae) an...
24446526 - Association of il1b -31c/t and il1ra variable number of an 86-bp tandem repeat with den...
19290876 - Mutation scanning using high-resolution melting.
|Type: Journal Article; Research Support, Non-U.S. Gov't Date: 2008-07-24|
|Title: BMC evolutionary biology Volume: 8 ISSN: 1471-2148 ISO Abbreviation: BMC Evol. Biol. Publication Date: 2008|
|Created Date: 2008-08-18 Completed Date: 2008-08-25 Revised Date: 2013-06-05|
Medline Journal Info:
|Nlm Unique ID: 100966975 Medline TA: BMC Evol Biol Country: England|
|Languages: eng Pagination: 217 Citation Subset: IM|
|Department of Biology and Physical Geography, University of British Columbia Okanagan, Kelowna, British Columbia V1V 1V7, Canada. firstname.lastname@example.org|
|APA/MLA Format Download EndNote Download BibTex|
Animals, Wild / genetics*
DNA, Mitochondrial / genetics*
Parakeets / genetics*
Sequence Analysis, DNA
Journal ID (nlm-ta): BMC Evol Biol
Publisher: BioMed Central
Copyright © 2008 Russello et al; licensee BioMed Central Ltd.
open-access: This is an Open Access article distributed under the terms of the Creative Commons Attribution License (), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
Received Day: 18 Month: 1 Year: 2008
Accepted Day: 24 Month: 7 Year: 2008
collection publication date: Year: 2008
Electronic publication date: Day: 24 Month: 7 Year: 2008
Volume: 8First Page: 217 Last Page: 217
Publisher Id: 1471-2148-8-217
PubMed Id: 18652686
|Genetic evidence links invasive monk parakeet populations in the United States to the international pet trade|
|Michael A Russello12||Email: email@example.com|
|Michael L Avery3||Email: Michael.L.Avery@aphis.usda.gov|
|Timothy F Wright4||Email: firstname.lastname@example.org|
1Department of Biology and Physical Geography, University of British Columbia Okanagan, Kelowna, British Columbia V1V 1V7, Canada
2Centre for Species at Risk and Habitat Studies, University of British Columbia Okanagan, Kelowna, British Columbia V1V 1V7, Canada
3USDA Wildlife Services, National Wildlife Research Center, Gainesville, Florida 32641, USA
4Department of Biology, New Mexico State University, Las Cruces, NM 88003, USA
The introduction of exotic species into native ecosystems has modified habitats, reduced species diversity and adversely altered ecosystem functioning across the globe . In the United States only habitat degradation poses a higher threat to endangered taxa . In other regions around the world, however, as many as 80% of endangered species are threatened due to pressures from non-native species . From an economic perspective, the environmental damage caused by the approximately 50,000 alien-invasive species in the United States, coupled with the costs of controlling these species, exceeds $120 billion per year . The severe ecological and economic impacts of invasive species render it imperative to understand the attributes that permit them to establish and spread within their expanded ranges.
The monk parakeet (Myiopsitta monachus) is one of the most successful parrot invaders . It has been historically regarded as an agricultural pest in its native range in South America, as noted by Charles Darwin during his voyage on the H.M.S. Beagle:
"A small green parrot (Conurus murinus; early synonym of M. monachus), with a grey breast, appears to prefer the tall trees on the islands to any other situation for its building-place. A number of nests are placed so close together as to form one great mass of sticks. These parrots always live in flocks, and commit great ravages on the corn-fields. I was told that near Colonia 2500 were killed in the course of one year." [, Chapter VII, p. 101]
Over the past century, the widespread introduction of Eucalyptus has facilitated the expansion of M. monachus populations in its native range [6-8]. In Argentina, this rapid increase in native population sizes has been implicated in the loss of 2–15% of sunflower and corn yields, crop damages estimated by some sources to be as high as US$1 billion per year [9,10]. In addition to rapid population growth within their endemic range, monk parakeets have become broadly established on four other continents, presumably due to their widespread presence in the international pet bird trade. Despite the rapid spread of monk parakeets around the globe and their potential as an agricultural pest, little is known about the geographical history of the invasions. Such information may provide important insights into the mechanisms of invasion success and potential for future range expansions.
The monk parakeet is distributed in its native range across the lowlands of South America, east of the Andes from Bolivia to Patagonia [, Figure 1]. Four subspecies are currently recognized based on geographical variation in wing length, bill size, body mass and plumage coloration . The nominate, M. m. monachus, is the largest of the four subspecies and is found in extreme southeast Brazil (Rio Grande do Sul), Uruguay and northeastern Argentina (provinces of Entre Rios, Santa Fe, Córdoba, south to northern Rio Negro). M. m. calita is distributed in western Argentina from Salta province south to Rio Negro and described as having bluer wings and a darker gray head. M. m. cotorra is distributed in southeast Bolivia (department of Tarija), Paraguay, southern Brazil (Mato Grosso do Sul), south to northern Argentina (provinces of Formosa and Chaco). M. m. cotorra has been reported as brighter green on the upper parts and less yellowish on the abdomen than M. m. calita , yet their general lack of distinctiveness in these characters and their similarity in size has brought their status as separate taxonomic entities into question . Lastly, M. m. luchsi is geographically and altitudinally isolated from the other subspecies, restricted to the arid, intermontane valleys of the east Andes in Bolivia, from southern Cochabamba to northern Chuquisaca . In addition, M. m. luchsi exhibits distinctive plumage coloration, reported as generally brighter than the other subspecies, with a bright yellow lower breast, paler underwings, a dark area at the base of the upper mandible, and a breast entirely pale grey without the barred effect observed in the other three subspecies . In contrast to the colonial, tree-nesting behavior of all other monk parakeets, M. m. luchsi build single-chambered nests on cliffs. These behavioral and morphological differences led del Hoyo  to elevate this group to allospecies status (Myiopsitta luchsi), a designation that is not widely recognized.
In addition to the South American populations, naturalized breeding populations of M. monachus have been established in such disparate regions as the United Kingdom, Puerto Rico, Kenya, Japan, Spain, Italy, Belgium, Czech Republic, and throughout the United States including growing populations in Florida, Texas and Connecticut [4,15-18]. The origin of initial invaders in the U.S. has been traced back to purposeful and accidental releases of individuals from the pet trade for which approximately 64,225 monk parakeets were imported between 1968–1972 alone . In general, U.S. naturalized populations are a collection of disjunct colonies, most common in southern and coastal regions, with an estimated 6,000 to 200,000 individuals in residence nationally . Once feared as a potentially devastating crop pest, M. monachus is still generally considered a moderate threat as populations continue to grow exponentially [15,19]. A less publicized impact of the monk parakeet invasion has been their preference for power structures as nesting substrates. In 2001, an estimated 1,027 power outages in south Florida were attributed to monk parakeet activities at an approximate cost of $585,000 . Moreover, the cost of nest removal alone in south Florida was estimated at $1.3 to $4.7 million over the past five years (2003–2007) . In addition to the financial impacts to energy providers and the communities they serve, monk parakeets may have ecological effects within the local ecosystems. Overall, these real and potential impacts have resulted in the presence of statewide controls or bans in over 15 states .
The objectives of this study were to identify the taxonomic and geographic source(s) of the invasive populations along the eastern seaboard of the United States. As the accuracy of monk parakeet taxonomy has been questioned, we initially conducted a genetic assessment of the biological validity of the four currently recognized subspecies of M. monachus by way of historical DNA analysis from museum specimens collected throughout the native range. This broad sampling of mitochondrial DNA lineage diversity in the endemic range provides a reference database by which to infer the origin and extent of mtDNA haplotype diversity in invasive populations in Florida, Connecticut, New Jersey and Rhode Island. Our results suggest that these invasive populations are derived from a localized area in eastern Argentina and bordering areas in Brazil and Uruguay within the described range of M. m. monachus, the most commonly exported subspecies for the international pet trade.
Toepad tissue was obtained from 73 museum specimens of Myiopsitta monachus representing all four subspecies (M. m. calita, n = 9; M. m. cotorra, n = 16; M. m. luchsi, n = 14; M. m. monachus, n = 38) courtesy of the American Museum of Natural History (AMNH). Blood samples from four individuals of M. m. monachus collected in Entre Rios, Argentina were also included in the sampling. In addition, feather or tissue samples were obtained from 64 individuals from four localities across the naturalized range of M. monachus in the eastern United States. Specifically, muscle tissue was obtained from individuals culled in Miami, Florida (n = 43) as part of a management program by the local electric utility company. Feathers were obtained from colonies sampled from different trees and areas in Bridgeport, Connecticut (n = 9) and Edgewater, New Jersey (n = 11). In addition, we sampled a single museum specimen collected in Kent County, Rhode Island (AMNH832643). Table 1 includes detailed collection information for all samples while Figure 1 plots the individual sampling localities within the native range of M. monachus.
DNA was extracted from blood, tissue, and feather samples using the DNeasy Tissue kit and manufacturer protocols (Qiagen, Inc.). Museum specimens were handled in a dedicated ancient DNA facility using a modified Qiagen DNeasy Tissue kit protocol . Other necessary precautions were taken to prevent and detect contamination by contemporary specimens, including use of extraction and PCR negative controls, PCR amplification of short, overlapping fragments (see below), and confirmation of all unique haplotype sequences by way of cloning .
A 558 basepair segment of the mitochondrial DNA (mtDNA) control region (CR) was amplified as a single fragment using external primers LGlu and CR522Rb  for the DNA extractions from blood and feather samples or, in the case of the DNA extractions from museum specimens, as a set of four overlapping fragments not exceeding 180 basepairs in length each [Lglu/MyiopCR1B (TGCCAATGGTTGCCCTAATAA); MyiopCR2A (GACATTGCATGCTCGTCCTA)/MyiopCR2B (TGGAATTGGAGAGGAGTGTTTT); MyiopCR3A (AGCAACTAAACCGAATGATCC)/MyiopCR3B (TGGGCCTGAAGCTAGTAACG); MyiopCR4A (CCACTCACGAGAAACCATCA)/CR522Rb]. All PCR reactions were carried out on an MJ Research DNA Engine thermal cycler in 25 μl reactions containing: ~20–50 ng of DNA, 10 mM Tris-HCl (pH 8.3), 50 mM KCl, 1.5 mM MgCl2, 200 μM dNTPs, 0.5 μM of each primer and 0.5 U of AmpliTaq Gold DNA polymerase (Applied Biosystems). Cycling conditions for all primer pairs consisted of 95°C for 10 minutes, 35 cycles of 95°C for 30 seconds, 50°C for 30 seconds, 72°C for 30 seconds, and a final extension of 72°C for 7 minutes. Double-stranded PCR products were sequenced using Big Dye 3.1 terminators on an ABI 3730 DNA sequencer (Applied Biosystems).
Previous work revealed duplication and concerted evolution of the control region in Amazona and Pionus parrots . Subsequent surveys of mtDNA gene order across the entire order of parrots via PCR across selected gene junctions has revealed that this duplication is absent in many parrot species, including M. monachus (Schirtzinger E, Gonzalez L, Eberhard JR, Graves G, Wright TF, unpublished data). Furthermore, long-range PCR followed by sequencing of the entire mtDNA genome of M. monachus has shown that it conforms to a typical avian gene order with a single control region (Schirtzinger E, Eberhard JR, Wright TF, unpublished data).
Haplotypic (h)  and nucleotide (π)  diversity estimates were calculated based on mtDNA CR sequences as executed in ARLEQUIN . Pairwise genetic distances were calculated in PAUP*4.0b10  assuming the HKY+G model of nucleotide substitution as selected according to the Akaike information criterion as implemented in Modeltest . Levels of genetic divergence between samples were calculated with the fixation index (PhiST)  as executed in ARLEQUIN . Because the HKY model is not implemented in ARLEQUIN the more inclusive Tamura-Nei (TrN)  model with the same parameters for ti/tv rate and α was used. Significance of PhiST for all possible pairwise population comparisons was assessed using 2,000 permutations. Tests for significant geographic structure among subspecies sampled across the native range were conducted using analysis of molecular variance (AMOVA) . MtDNA CR sequence alignments for all four subspecies were further employed to identify diagnostic nucleotide sites by means of population aggregation analysis . The presence of characters fixed within and differing among populations was used as evidence to diagnose distinct units.
Sequences were unambiguously aligned in Clustal X  employing default settings for gap opening and extension costs. Genealogical relationships among all sampled haplotypes throughout the native range were reconstructed as a haplotype network using the statistical parsimony method of Templeton et al.  as implemented in TCS, version 1.06 . Gaps were treated as a 5th character state. Networks are especially appropriate for inferring intraspecific gene genealogies because of the potential for extant ancestral nodes and multifurcating relationships .
A Bayesian haplotype tree was reconstructed using MrBayes 3.1  assuming the HKY+G model of nucleotide substitution as selected by Modeltest  as described above. The orange-chinned parakeet (Brotogeris jugularis) was used as an outgroup to root the tree, as previous phylogenetic studies have revealed species from this genus to be sister to M. monachus [37,38]. The Bayesian phylogenetic analysis ran four simultaneous chains for 2.0 × 106 total generations, each using a random tree as a starting point, the default heating scheme, and saving a tree every 100 generations for a total 20,000 trees. The first 2,000 trees were discarded as burn-in samples and the remaining 18,000 trees were used to construct a majority-rule consensus tree and derive posterior probability values. Violation of a criterion of monophyly was used to indicate incorrect taxonomic assignment.
A total of 17 mtDNA CR haplotypes were recovered among the 77 individuals sampled from across the native range of the four described subspecies of M. monachus (GenBank Accession No. EU545521-EU545537). The number of haplotypes identified ranged from four (M. m. calita) to eight (M. m. monachus), with levels of haplotypic and nucleotide diversity relatively consistent across the subspecies (Table 2). Of the 17 detected haplotypes, three were shared among a combination of M. m. calita, M. m. cotorra, and M. m. monachus. One shared haplotype was widely distributed, sampled in individuals from all three of these subspecies in disparate localities ranging from northern (Tucumán province) and central (Entre Rios and Mendoza provinces) Argentina, to Concepción, Paraguay and Mato Grosso, Brazil. All five haplotypes recovered for M. m. luchsi in Bolivia were unique to that subspecies. Overall, sequence divergence among M. monachus haplotypes recovered from the four subspecies ranged from 0.20% to 1.66% (luchsi01/calita02) based on HKY+G distances.
Genetic variation across the samples was highly structured with significant levels of genetic variation distributed among, rather than within, the four M. monachus subspecies (p < 0.0001; Table 3a). When the Bolivian M. m. luchsi was removed from the AMOVA, the results were reversed, with the vast majority of variation distributed within (96.41%) rather than among (3.59) subspecies (Table 3b). A similar pattern was revealed by the fixation indices, with all pairwise comparisons involving M. m. luchsi highly significant (Table 3c). None of the pairwise comparisons of M. m. calita, M. m. cotorra and M. m. monachus approached significance. Likewise, M. m. luchsi was diagnosably distinct from each of the other three subspecies, with the number of diagnostic characters detected ranging from three (M. m. cotorra, M. m. monachus) to five (M. m. calita) across the 558 basepairs of the mtDNA CR (Table 3c).
A single haplotype network was reconstructed within which all haplotypes had a 95% probability of being parsimoniously connected (Figure 2). Overall, the network was characterized by reticulation and little structure to the recovered relationships (Figure 2). The only distinct clustering was of the four M. m. luchsi haplotypes sampled in Bolivia, which were three to four steps different than the nearest M. m. monachus or M. m. cotorra haplotypes (Figure 2). The remaining haplotypes constituted a mixed assemblage, exhibiting neither geographic structure nor clustering patterns consistent with currently described subspecies boundaries. Results of a Bayesian phylogenetic analysis mirrored those of the haplotype network, reconstructing a well-supported (posterior probability = 90; Figure 3), monophyletic M. m. luchsi with the remaining three subspecies forming a paraphyletic assemblage.
Three haplotypes were recovered from the 64 individuals sampled in populations in the eastern United States, all of which were identical to haplotypes detected in the native range of M. monachus. The most common haplotype, detected in the naturalized range at a frequency of 0.63, was one previously found unique to M. m. monachus (monachus01; Figure 3). Initially sampled in Entre Rios, Argentina and Rio Grande do Sul, Brazil, this haplotype was fixed in the Bridgeport, CT (n = 9) and Kent County, RI (n = 1) samplings and was likewise detected at high frequencies in the Miami, FL (0.56) and Edgewater, NJ (0.55) populations. A second high frequency haplotype unique to the M. m. monachus subspecies (monachus02; Figure 3) was found in the Miami, FL (0.40) and Edgewater, NJ (0.45) populations. In the native range, the monachus02 haplotype was recovered over a wide geographic area, found in Rio Grande do Sul, Brazil, Soriano, Uruguay, and throughout sampling localities in northern and central Argentina (Figure 3). Lastly, a haplotype shared by M. m. monachus, M. m. calita and M. m. cotorra in the native range (shared01; Figure 3) was also found at very low frequency in the Miami, FL population (0.03).
In this study we employed extensive geographic sampling and historical DNA analysis to describe patterns of genetic variation among subspecies of the invasive monk parakeet in its native range, assess current taxonomic designations, and infer the source(s) of introduced populations in the United States.
The biological relevance of subspecies has been widely debated since the 1950s . Ernst Mayr, who wrote the most influential book on speciation analysis , also grappled with the concept of subspecies. Although in early writings he clearly assigned evolutionary status to subspecies , later work directly acknowledged the subjectivity associated with this level of taxonomic classification, explicitly stating that subspecies are not units of evolution . Continentally distributed avian subspecies are a prime example, with a recent survey finding that 97% lack the population genetic structure indicative of historically independent units .
We share the view of many that an accurate taxonomy should reflect evolutionary history. The phylogenetic species concept (PSC) offers such an approach, directly linking patterns of evolution with species status . Under the PSC sensu Cracraft , a species is the smallest diagnosable cluster of individual organisms within which there is a parental pattern of ancestry and descent. Applied to the monk parakeets, three (M. m. calita, M. m. cotorra, M. m. monachus) of four M. monachus subspecies lacked diagnostic character support and violated a criterion of monophyly. The absence of genetic distinctiveness of these three taxa reflects the uncertainty surrounding the differences in size and plumage characteristics upon which they were initially described [11,13]. The lack of discrete morphological character differences in tandem with the results presented here suggest that M. m. calita, M. m. cotorra and M. m. monachus are in need of formal taxonomic revision.
In contrast to the uncertainty associated with the descriptions of the other three subspecies, the controversy surrounding M. m. luchsi has been related to its relative distinctiveness, and proposed elevation to allospecies status . Restricted to the intermontane valleys in Bolivia, M. m. luchsi is morphologically distinct from the other subspecies , including M. m. cotorra, despite the fact that their known ranges come within 175 km of each other . Moreover, M. m. luchsi is altitudinally distributed between 1300–3000 m, in sharp contrast to other monk parakeet taxa, which are routinely found below 1000 m. Another unique characteristic involves the cliff-nesting behavior of M. m. luchsi, which contrasts with the colonial, tree-nesting exhibited across the remainder of the range of M. monachus. This assorted evidence has been used to elevate M. luchsi to allospecies status, forming a superspecies with the remaining taxa of M. monachus . Although this taxonomic revision is not generally recognized, the results of the current study further highlight the uniqueness of this taxon. In addition to displaying between three and five diagnostic molecular characters relative to M. m. calita, M. m. cotorra, and M. m. monachus (Table 2), the Bolivian luchsi formed a well-supported, monophyletic group based on the mtDNA control region sequence data (Figure 3). Collectively, the morphological, behavioral and genetic data support M. luchsi as a distinct, phylogenetic species [32,43] and suggest that a formal taxonomic revision is in order.
Over the past 35 years, monk parakeets have been recorded on U.S. Christmas Bird Counts in 14 states: Connecticut, Delaware, Florida, Georgia, Illinois, Massachusetts, Nebraska, New Jersey, New York, Ohio, Oregon, Pennsylvania, Texas, and Washington D.C./Virginia . Other states where monk parakeet nesting has been observed include Alabama , California , Louisiana , North Carolina , South Carolina , and Rhode Island . The United States Fish and Wildlife Service conducted an eradication campaign from 1970 to 1975 that effectively eliminated populations in California and reduced the naturalized range of monk parakeets in the U.S. to seven localities in five states . Since 1975, M. monachus populations in the U.S. grew exponentially and spread throughout the country to its present distribution . Currently, two of the largest naturalized populations of M. monachus reside in Florida and southern Connecticut. Both appear to be expanding in size and geographic distribution. The Florida population, in particular, continues to increase at an exponential rate, with a recent study by Pruett-Jones et al.  estimating a statewide population size of 18,025 to 32,044.
Despite multiple introductions and the widespread distribution of M. monachus in the U.S., we detected a low level of haplotype diversity across four sampling localities in Connecticut, Florida, New Jersey and Rhode Island. Only three different haplotypes were recovered, the most common of which (monachus01; 0.63) was found in all four localities and was identical to a haplotype sampled from M. m. monachus in a localized area in eastern Argentina in Entre Rios to Rio Grande do Sul, Brazil on the Uruguayan border. The other high-frequency haplotype in the naturalized range (monachus02; 0.34) was also specific to M. m. monachus and was likewise sampled in Rio Grande do Sul, Brazil, Soriano, Uruguay and a number of localities in central and northern Argentina. These results are consistent with preliminary morphometric analyses (M. Avery, unpublished) as well as trapping records that indicate that the vast majority of birds captured for the pet trade were M. m. monachus exported from eastern Argentina and Uruguay . The concordance between the trapping records and our genetic results support the idea that the invasion of monk parakeets has been facilitated, at least initially, by their widespread presence in the international pet bird trade. As a likely source of large and repeated release events, the international trade in parrots may have historically exerted significant propagule pressure, generally a key determinant of invasion success in birds and other taxa [52,53]. Nuclear data and additional geographical sampling may provide important sources of historical information for further testing this hypothesis.
The Wild Bird Conservation Act of 1992 prohibits the importation of monk parakeets into the United States, reducing the chances of future introductions of wild-caught individuals . However, monk parakeets, known as Quaker parakeets in the pet trade, remain one of the most popular cage birds and are widely bred and sold by aviculturists in the U.S. At the state level, local laws vary, with some states banning the possession of monk parakeets while others placing no restrictions on them. Nevertheless, this domestic trade in monk parakeets remains the most likely source of introductions into states not currently reporting self-sustaining breeding populations .
By and large, it is likely that non-native populations of monk parakeets will continue to grow. One reason is that the popularity of monk parakeets has in recent years extended to introduced populations. Efforts to remove birds and nests from electric utility structures in Connecticut, Illinois, Florida, Washington and New Jersey have often met with substantial resistance by a vocal subset of the local communities. In addition, current control strategies have not effectively prevented the establishment and continued growth of naturalized populations in the U.S. and Western Europe, most notably in Spain . Unlike other psittacines, monk parakeets are not constrained by the availability of nesting cavities. Rather, they construct nests of sticks and branches and they tend to select man-made structures as nesting substrates . Furthermore, by exploiting feeding opportunities provided by humans, monk parakeets persist in even cold temperate winters . Population growth and expansion seems assured, as impractically large management efforts would be needed to reverse the trend . Consequently, broader understanding of the mechanisms of monk parakeet invasion success and local adaptation constitute important areas for future basic and applied research.
MAR designed the study, collected samples (museum collections, Connecticut and Rhode Island invasive populations), carried out the molecular studies, performed data analyses, and drafted the manuscript. MLA facilitated sample collection in the U.S. populations, aided in interpretation of results, and helped draft the manuscript. TFW participated in the design of the study, aided in interpretation of results, and helped draft the manuscript. All authors read and approved the final manuscript.
This work was funded by an American Philosophical Society grant to MR. We thank the American Museum of Natural History for granting access to their collections and P. Sweet for his assistance sampling the specimens. J. Eberhard provided field-collected samples from Argentina and offered helpful comments on the manuscript. J. Wright, D. Hoffmeier, and J. Lindsay with Florida Power and Light Company provided samples from south Florida. The work further benefited from discussions with G. Amato and A. Caccone. C. Hyseni, S. Glaberman and E. Moseman aided in data and sample collection. An internal grant awarded by UBC Okanagan (MR) assisted in covering publication costs.
|Pimentel D,Zuniga R,Morrison D. Update on the environmental and economic costs associated with alien-invasive species in the United StatesEcological Economics 2005;52:273–288.|
|Wilcove DS,Rothstein D,Dubow J,Phillips A,Losos E. Quantifying threats to imperiled species in the United StatesBioscience 1998;48:607–615. [doi: 10.2307/1313420]|
|Armstrong S. Rare Plants Protect Capes Water-SuppliesNew Scientist 1995;145:8–8. [pmid: 11659935]|
|Lever C. Naturalized birds of the world (3rd). 19873rd. London , Longman Scientific & Technical;|
|Darwin C. Journal of researches into the natural history and geology of the countries visited during the voyage of H.M.S. Beagle. 1839London , Henry Colburn;|
|Gibson E. Further ornithological notes from the neighbourhood of Cape San Antonio, Buenos Ayres. Part II.Ibis 1919;11:495–537.|
|Forshaw JM. Parrots of the world. 2006New Jersey , Princeton University Press;|
|Navarro JL,Martella MB,Bucher EH. Breeding-Season and Productivity of Monk Parakeets in Cordoba, ArgentinaWilson Bulletin 1992;104:413–424.|
|Davis LR. Johnson WV, Marsh REThe monk parakeet: a potential threat to agriculture: University of California, Davis. 1974:253–256.|
|Iriarte JA,Lobos GA,Jaksic FM. Invasive vertebrate species in Chile and their control and monitoring by governmental agenciesRevista Chilena De Historia Natural 2005;78:143–154.|
|Juniper T,Parr M. Parrots: a guide to parrots of the world. 1998New Haven , Yale University Press;|
|Peters J. Checklist of birds of the world 1937;IIICambridge , Harvard University Press;|
|Spreyer MF,Bucher EH. Poole A, Gill FMonk Parakeet (Myiopsitta monachus)The Birds of North America, No 322. 1998Philadelphia , The Birds of North America, Inc.;|
|del Hoyo J,Elliott A,Sargatal J. Handbook of the birds of the world, Volume 4. Sandgrouse to Cuckoos. 1997Cambridge , Birdlife International; [pmid: 9021763]|
|Pruett-Jones S,Tarvin KA. Baker RO, Crabb ACMonk Parakeets in the United States: population growth and regional patterns of distribution: University of California, Davis. 1998:55–58.|
|Butler CJ. Feral parrots in the continental United States and United Kingdom: Past, present, and futureJournal of Avian Medicine and Surgery 2005;19:142–149. [doi: 10.1647/183]|
|Munoz AR,Real R. Assessing the potential range expansion of the exotic monk parakeet in SpainDiversity and Distributions 2006;12:656–665. [doi: 10.1111/j.1472-4642.2006.00272.x]|
|Buhrman-Deever SC,Rappaport AR,Bradbury JW. Geographic variation in contact calls of feral North American populations of the Monk ParakeetCondor 2007;109:389–398. [doi: 10.1650/0010-5422(2007)109[389:GVICCO]2.0.CO;2]|
|Van Bael S,Pruett-Jones S. Exponential population growth of Monk Parakeets in the United StatesWilson Bulletin 1996;108:584–588.|
|Avery ML,Greiner EC,Lindsay JR,Newman JR,Pruett-Jones S. Timm RM, Schmidt RHMonk parakeet management at electric utility facilities in south Florida: University of California, Davis. 2002:140–145.|
|Avery ML,Yoder CA,Tilman EA. Diazacon inhibits reproduction in invasive monk parakeet populationsJournal of Wildlife Management In press|
|Stafford T. Pest risk assessment for the monk parakeet in Oregon. 2003 Oregon Department of Agriculture;|
|Russello MA,Brazaitis P,Gratten J,Watkins-Colwell GJ,Caccone A. Molecular assessment of the genetic integrity, distinctiveness and phylogeographic context of the Saltwater crocodile (Crocodylus porosus) on PalauConservation Genetics 2007;8:777–787. [doi: 10.1007/s10592-006-9225-7]|
|Wandeler P,Hoeck PEA,Keller LF. Back to the future: museum specimens in population geneticsTrends in Ecology & Evolution 2007;22:634–642. [pmid: 17988758] [doi: 10.1016/j.tree.2007.08.017]|
|Eberhard JR,Wright TF,Bermingham E. Duplication and concerted evolution of the mitochondrial control region in the parrot genus AmazonaMolecular Biology and Evolution 2001;18:1330–1342. [pmid: 11420371]|
|Nei M. Molecular Evolutionary Genetics. 1987New York , Columbia University Press;|
|Schneider S,Roessli D,Excoffier L. Arlequin v2.000: A software for population genetics data analysis. 2000Switerland , Genetics and Biometry Laboratory, University of Geneva;|
|Swofford DL. PAUP*: Phylogenetic Analysis Using Parsimony, v.4.0b10 . 2002Illinois Natural History Survey, Champaign ; [pmid: 12504223]|
|Posada D,Crandall KA. MODELTEST: testing the model of DNA substitutionBioinformatics 1998;14:817–818. [pmid: 9918953] [doi: 10.1093/bioinformatics/14.9.817]|
|Excoffier L,Smouse PE,Quattro JM. Analysis of Molecular Variance Inferred from Metric Distances among DNA Haplotypes - Application to Human Mitochondrial-DNA Restriction DataGenetics 1992;131:479–491. [pmid: 1644282]|
|Tamura K,Nei M. Estimation of the Number of Nucleotide Substitutions in the Control Region of Mitochondrial-DNA in Humans and ChimpanzeesMolecular Biology and Evolution 1993;10:512–526. [pmid: 8336541]|
|Davis JI,Nixon KC. Populations, Genetic-Variation, and the Delimitation of Phylogenetic SpeciesSystematic Biology 1992;41:421–435. [doi: 10.2307/2992584]|
|Thompson JD,Gibson TJ,Plewniak F,Jeanmougin F,Higgins DG. The CLUSTAL_X windows interface: flexible strategies for multiple sequence alignment aided by quality analysis toolsNucleic Acids Research 1997;25:4876–4882. [pmid: 9396791] [doi: 10.1093/nar/25.24.4876]|
|Templeton AR,Crandall KA,Sing CF. A Cladistic-Analysis of Phenotypic Associations with Haplotypes Inferred from Restriction Endonuclease Mapping and DNA-Sequence Data .3. Cladogram EstimationGenetics 1992;132:619–633. [pmid: 1385266]|
|Clement M,Posada D,Crandall KA. TCS: a computer program to estimate gene genealogiesMolecular Ecology 2000;9:1657–1659. [pmid: 11050560] [doi: 10.1046/j.1365-294x.2000.01020.x]|
|Ronquist F,Huelsenbeck JP. MrBayes 3: Bayesian phylogenetic inference under mixed modelsBioinformatics 2003;19:1572–1574. [pmid: 12912839] [doi: 10.1093/bioinformatics/btg180]|
|de Kloet RS,de Kloet SR. The evolution of the spindlin gene in birds: Sequence analysis of an intron of the spindlin W and Z gene reveals four major divisions of the PsittaciformesMolecular Phylogenetics and Evolution 2005;36:706–721. [pmid: 16099384] [doi: 10.1016/j.ympev.2005.03.013]|
|Tavares ES,Baker AJ,Pereira SL,Miyaki CY. Phylogenetic relationships and historical biogeography of Neotropical parrots (Psittaciformes : Psittacidae : Arini) inferred from mitochondrial and nuclear DNA sequencesSystematic Biology 2006;55:454–470. [pmid: 16861209] [doi: 10.1080/10635150600697390]|
|Wilson EO,Brown WL. The subspecies concept and its taxonomic applicationSystematic Zoology 1953;2:97–111. [doi: 10.2307/2411818]|
|Mayr E. Systematics and the Origin of Species. 1942New York , Columbia University Press;|
|Mayr E. Animal Species and Evolution. 1963Cambridge , Harvard University Press;|
|Zink RM. The role of subspecies in obscuring avian biological diversity and misleading conservation policyProceedings of the Royal Society of London Series B-Biological Sciences 2004;271:561–564. [doi: 10.1098/rspb.2003.2617]|
|Cracraft J. Species concepts and speciation analysisCurrent Ornithology 1983;1:159–187.|
|Lanning DV. Distribution and Nest Sites of the Monk Parakeet in BoliviaWilson Bulletin 1991;103:366–372.|
|Summerour B. Monk parakeets (Myiopsitta monachus) nesting in HuntsvilleAlabama Birdlife 1990;37:9–10.|
|Hardy JW. Feral Exotic Birds in Southern-CaliforniaWilson Bulletin 1973;85:506–512.|
|Simpson MB,Ruiz RC. Monk-Parakeets Breeding in Buncombe-County, North-CarolinaWilson Bulletin 1974;86:171–172.|
|McNair DB,Post W. Supplement to status and distribution of South Carolina birdsCharleston Museum Ornithological Contribution No 8. 1993|
|Enser RW. The atlas of breeding birds in Rhode Island (1982-1987). 1992 Rhode Island Department of Environmental Management;|
|Neidermeyer WJ,Hickey JJ. The monk parakeet in the United States, 1970-75American Birds 1977;31:273–278.|
|Pruett-Jones S,Newman JR,Newman CM,Lindsay JR. Population growth of monk parakeets in FloridaFlorida Field Naturalist 2005;33:1–14.|
|Blackburn TM,Duncan RP. Determinants of establishment success in introduced birdsNature 2001;414:195–197. [pmid: 11700555] [doi: 10.1038/35102557]|
|Lockwood JL,Cassey P,Blackburn T. The role of propagule pressure in explaining species invasionsTrends in Ecology & Evolution 2005;20:223–228. [pmid: 16701373] [doi: 10.1016/j.tree.2005.02.004]|
|Wild Bird Conservation Act. 1992U.S.A. ;|
|Avery ML,Lindsay JR,Newman JR,Pruett-Jones S,Tillman EA. Feare CJ, Cowan DPReducing monk parakeet impacts to electric utility facilities in south FloridaAdvances in Vertebrate Pest Management 2006;4Fürth , Filander Verlag; :125–136.|
|South JM,Pruett-Jones S. Patterns of flock size, diet, and vigilance of naturalized Monk Parakeets in Hyde Park, ChicagoCondor 2000;102:848–854. [doi: 10.1650/0010-5422(2000)102[0848:POFSDA]2.0.CO;2]|
|Pruett-Jones S,Newman JR,Newman CM,Avery ML,Lindsay JR. Population viability analysis of monk parakeets in the United States and examination of alternative management strategiesHuman-Wildlife Conflicts 2007;1:35–44.|
[Figure ID: F1]
Distribution of Myiopsitta monachus across its native range in South America. Alternative shading denotes the individual ranges of the four subspecies [redrawn from ] including M. m. monachus (light gray), M. m. calita (black), M. m. cotorra (dark gray), and M. m. luchsi (striped). Localities of specimens sampled for this study are indicated by dots, with associated abbreviations following Table 1.
[Figure ID: F2]
Network showing genealogical relationships among Myiopsitta monachus haplotypes sampled in the native range. Haplotypes are connected with a 95% confidence limit. The size of each oval is proportional to the frequency of the haplotype in the analysis. White dots represent mutational steps separating the observed haplotypes. Different shades represent the proportion of individuals of each subspecies exhibiting that particular haplotype (colors as in Figure 1).
[Figure ID: F3]
Bayesian haplotype tree depicting relationships among sampled Myiopsitta monachus haplotypes relative to their geographic and taxonomic distributions. The names of each haplotype are as in Figure 2. Bayesian posterior probabilities (> 50%) are indicated above the branches. Each column in the associated table is a locality sorted by country with abbreviations following Table 1. Each row is a haplotype according to its placement in the tree on the left; the number of individuals at that sampling locality exhibiting that particular haplotype is indicated in each cell. Shading represents the subspecies designation for the distribution of haplotypes according to Figure 1. Bolded italicized numbers indicate the distribution of individuals collected in the naturalized range in the United States. Total number of sampled individuals exhibiting each haplotype (N) is denoted in the last column. For illustration purposes, accurate branch lengths leading to the outgroup are not shown (indicated by dashed line).
Sampling of Myiopsitta monachus subspecies in native and naturalized ranges
|M. m. calita||3||Argentina||Santiago del Estero||SE||140653, 474808–474809|
|M. m. cotorra||5||Brazil||Mato Grosso||MG||127356–127359, 474802|
|11||Paraguay||Concepción||PRG||149404, 320771–320777, 748687–748688, 811356|
|M. m. luchsi||12||Bolivia||Chuquisaca||CH||139094–139096,139098–139106|
|M. m. monachus||2||Argentina||Santiago del Estero||SE||140649, 140651|
|14*||Argentina||Entre Rios||ER||779017–779019, 779025, 779037, 779059–779061, 779065, 779083|
|13||Argentina||Corrientes||CR||793580–793581, 793586–793587, 793589, 793597–793598, 793603, 793605, 793616, 793631, 793633, 793640|
|6||Brazil||Rio Grande do Sul||RGS||321247–321249, 321560–321562|
|11||United States||New Jersey||NJ||n/a|
|1||United States||Rhode Island||RI||832643|
*Sampling includes four field-collected samples.
† Accession numbers of specimens sampled in the collections at the American Museum of Natural History or, in the case of Florida individuals, at the USDA National Wildlife Research Center.
†† Individuals sampled in the naturalized populations are of unknown taxonomic affinity.
Genetic variation within Myiopsitta monachus subspecies
|Haplotypes†||Diversity, h||Diversity, π|
|M. m. calita||9||4||0.58||0.0031|
|M. m. cotorra||16||5||0.73||0.0020|
|M. m. monachus||38||8||0.77||0.0028|
|M. m. luchsi||14||5||0.66||0.0015|
†Results based on 558 base pairs of the mtDNA control region. All haplotypes recovered for each subspecies are considered. Three haplotypes were shared among M. m. calita, M. m. cotorra, and M. m. monachus. All haplotypes sampled in U.S.A. were also recovered in native range (see text).
‡Values in parentheses are the standard errors for h and π.
Genetic divergence among Myiopsitta monachus subspecies
|a. Analysis of molecular variance including all subspecies|
|Subspecies||Source of||d.f.||% of||P-value|
|M. m. calita||Among||3||61.23||<0.0001|
|M. m. cotorra||Within||74||38.77|
|M. m. monachus||Total||77|
|M. m. luchsi|
|b. Analysis of molecular variance excluding M. m. luchsi|
|Subspecies||Source of||d.f.||% of||P-value|
|M. m. calita||Among||2||3.59||0.1369|
|M. m. cotorra||Within||61||96.41|
|M. m. monachus||Total||63|
|c. Diagnostic characters and fixation indices*|
|Subspecies||M. m. calita||M. m. cotorra||M. m. monachus||M. m. luchsi|
|M. m. calita||-||-0.0370||0.0500||0.8062**|
|M. m. cotorra||0||-||0.0422||0.8266**|
|M. m. monachus||0||0||-||0.7757**|
|M. m. luchsi||5||3||3||-|
‡ Among populations, within populations or total.
* Number of diagnostic characters (below diagonal) and PhiST (above diagonal) based on mtDNA control region sequence data.
** Indicates statistical significance (p < 0.001).
Previous Document: Evaluation of the bacterial diversity in the feces of cattle using 16S rDNA bacterial tag-encoded FL...
Next Document: Deciphering downstream gene targets of PI3K/mTOR/p70S6K pathway in breast cancer.