Low mtDNA diversity among widespread Australian diamondback moth Plutella xylostella （L,） suggests isolation and a founder effect

Populations of Australian diamondback moth （DBM） Plutella xylostella （L.）, a serious pest of cruciferous crops, display extremely low levels of genetic differentiation across Australia and New Zealand sample locations, as determined previously using microsatellite markers. These data suggest high levels of contemporary gene flow that is consistent with Australian DBM being a vagile species. Here we examine Australian DBM samples for haplotype variation using the mitochondrial DNA sequences of a 257 bp fragment of the CO1 gene. We compare this variation to equivalent mtDNA sequence variation in samples from New Zealand, Kenya and Korea. Using 42 moths collected throughout Australia we show that Australian DBM have both low mtDNA haplotype and nucleotide diversities. The three Australian haplotypes detected are closely related and they cluster with the common haplotype group from Indonesia. In addition the Australian haplotype frequency distribution resembled more that from Indonesia than that from Kenya or Korea. These data are consistent with an original strong Australian/New Zealand founder effect, from a south-eastern Asian source, with subsequent continued isolation. In a single season, the frequency of PXMt01, the most common Australian haplotype, was estimated at 15 locations spread across southern Australia and New Zealand using a polymerase chain reaction BiPASA method. The PXMt01 haplotype frequency variation was heterogenous, suggesting a small degree of population isolation that was not detected using microsatellites. Differentiation was not a function of geographical distance. These data suggest transient and sporadic local colonisation events by small numbers of founding females.     

High microsatellite diversity and differential structuring among populations of the introduced common brushtail possum, Trichosurus vulpecula, in New Zealand

An understanding of genetic variation and structure of pest populations has the potential to improve the efficiency of measures to control them. Genetic analysis was undertaken at five microsatellite loci in four native Australian and 14 introduced New Zealand populations of the common brushtail possum Trichosurus vulpecula in order to document these parameters. Genetic variation in New Zealand populations, and phylogenetic relationships among Australian and New Zealand populations, were largely predicted by the recorded introduction history. Populations on the two main islands of New Zealand had only slightly lower genetic diversity than did Australian populations, except that allelic richness on the South Is. was significantly lower. Diversity was higher in North Is. than in South Is. populations (although not significantly so) and mainland New Zealand populations as a group were significantly more diverse than offshore islands that represented secondary population size bottlenecks. In phylogenetic analyses South Is. and offshore island populations grouped with Tasmania, while North Is. populations grouped either with mainland Australia or were intermediate between the two Australian sources. This scheme was supported by admixture coefficients showing that North and South Is./offshore island populations were largely mainland Australian and Tasmanian in origin, respectively. Population structure differed markedly between the North and South Islands: populations were typically more genetically differentiated on the former than the latter, which also showed significant isolation-by-distance. Substantial linkage disequilibrium in most sampled New Zealand but no Australian population between microsatellite loci Tv16 and Tv27 suggests they may be physically linked.     

Genetic differentiation among populations of the shortfinned eel Anguilla australis from East Australia and New Zealand

Variability at seven microsatellite loci was used to survey the genetic population structure of the shortfinned eel Anguilla australis . Samples were collected from six estuaries along the east coast of Australia and from three estuaries around New Zealand. Hierarchical analysis of molecular variance of the five loci with good fit to Hardy–Weinberg genotypic proportions detected highly significant differences among samples ( F _ST= 0·016, P < 0·001). The fixation index between countries ( F _CT= 0·012, P < 0·001) was more than double the index among samples within countries ( F _SC= 0·005, P < 0·05). An unweighted pair-group method with arithmetic mean (UPGMA) tree also supported the separation of Australian and New Zealand populations, as did assignment tests, which correctly assigned 80 and 84% of the individuals to Australia and New Zealand, respectively. Isolation-by-distance appeared among samples overall ( r = 0·807, P < 0·001), but not among samples within countries ( r = 0·027, P > 0·05 in Australia; r = 0·762, P > 0·05 in New Zealand). These findings indicate that populations of A. australis in East Australia and in New Zealand may be reproductively isolated from one another. Genetic differentiation among populations of A. australis was two- to 10-fold higher than that among populations of other temperate eels in the North Atlantic Ocean, suggesting that two group of A. australis may reflect sub-species. Anguilla australis in the two countries have different genetic structures and thus require separate management. Genetic isolation between Australian and New Zealand populations indicates that juveniles recruit independently into these two regions from geographically or temporally isolated spawning areas.     

Genetic architecture in codling moth populations: comparison between microsatellite and insecticide resistance markers

The codling moth, Cydia pomonella, is renowned for developing resistance to insecticides and causing significant economic damage to pome fruits worldwide. In spite of its economic importance, little is known about the patterns of movement of this pest and the effects of insecticide treatment on the population genetic structure. Here, we investigated the genetic structure of the pest in 27 orchards from France, Italy, Armenia and Chile at seven microsatellite loci and two resistance markers [biochemical activity of cytochrome P450 oxidases and proportion of knockdown resistance (kdr) alleles in the sodium channel gene]. According to the microsatellite loci, we detected isolation by distance at the supranational scale but found no evidence of geographical structure among the 24 French orchards, which were mainly structured by the intensity of the insecticide treatments. Similarly, the highest levels of metabolic resistance associated with activity of the cytochrome P450 oxidases were detected in the most treated orchards. The kdr alleles were observed in southern France and Armenia where the pyrethroid insecticides were or have been intensively sprayed. The intensity of the insecticide treatments marginally affected the allelic richness in each orchard, but not the level of inbreeding. These results suggest important and high-distance gene flow among the codling moth populations, which were mainly structured according to the history of insecticide applications. Differences in mutation-migration-drift equilibrium among treated and untreated orchards also suggest that insecticide applications are the main force regulating the local dynamics of codling moth populations.     

Is there genetic structure in populations of Helicoverpa armigera from Australia?

Recent studies suggest that populations of the pest moth Helicoverpa armigera (Hübner) (Lepidoptera: Noctuidae) may be genetically differentiated over short distances and time periods within Queensland, Australia. To test for genetic structure in another region of Australia, we characterized population differentiation in Victorian samples of H. armigera using eight microsatellite loci. We found no evidence of genetic structure among samples from different locations or among samples collected at different times. Moreover, Victorian samples were not differentiated from other samples of H. armigera from Queensland and New Zealand. All samples showed substantial deviations from Hardy–Weinberg equilibrium, suggesting a high frequency of null alleles typically found in microsatellites of Lepidoptera. These results indicate that populations of H. armigera are not strongly structured among regions in south‐eastern Australia.     

ISSR-PCR: tool for discrimination and genetic structure analysis of Plutella xylostella populations native to different geographical areas

The diamondback moth (DBM), Plutella xylostella (L.) is considered as the most destructive pest of Brassicaceae crops world-wide. Its migratory capacities and development of insecticide resistance in many populations leads to more difficulties for population management. To control movement of populations and apparitions of resistance carried by resistant migrant individuals, populations must be identified using genetic markers. Here, seven different ISSR markers have been tested as a tool for population discrimination and genetic variations among 19 DBM populations from Canada, USA, Brazil, Martinique Island, France, Romania, Austria, Uzbekistan, Egypt, Benin, South Africa, Réunion Island, Hong Kong, Laos, Japan and four localities in Australia were assessed. Two classification methods were tested and compared: a common method of genetic distance analyses and a novel method based on an advanced statistical method of the Artificial Neural Networks' family, the Self-Organizing Map (SOM). The 188 loci selected revealed a very high variability between populations with a total polymorphism of 100% and a global coefficient of gene differentiation estimated by the Nei's index (Gst) of 0.238. Nevertheless, the largest part of variability was expressed among individuals within populations (AMOVA: 73.71% and mean polymorphism of 94% within populations). Genetic differentiation among the DBM populations did not reflect geographical distances between them. The two classification methods have given excellent results with less than 1.3% of misclassified individuals. The origin of the high genetic differentiation and efficiency of the two classification methods are discussed.     

Phylogeographic Structure in Penguin Ticks across an Ocean Basin Indicates Allopatric Divergence and Rare Trans-Oceanic Dispersal

The association of ticks (Acarina) and seabirds provides an intriguing system for assessing the influence of long-distance dispersal on the evolution of parasitic species. Recent research has focused on host-parasite evolutionary relationships and dispersal capacity of ticks parasitising flighted seabirds. Evolutionary research on the ticks of non-flighted seabirds is, in contrast, scarce. We conducted the first phylogeographic investigation of a hard tick species (Ixodes eudyptidis) that parasitises the Little Blue Penguin (Eudyptula minor). Using one nuclear (28S) and two mitochondrial (COI and 16S) markers, we assessed genetic diversity among several populations in Australia and a single population on the South Island of New Zealand. Our results reveal two deeply divergent lineages, possibly representing different species: one comprising all New Zealand samples and some from Australia, and the other representing all other samples from Australian sites. No significant population differentiation was observed among any Australian sites from within each major clade, even those separated by hundreds of kilometres of coastline. In contrast, the New Zealand population was significantly different to all samples from Australia. Our phylogenetic results suggest that the New Zealand and Australian populations are effectively isolated from each other; although rare long-distance dispersal events must occur, these are insufficient to maintain trans-Tasman gene flow. Despite the evidence for limited dispersal of penguin ticks between Australia and New Zealand, we found no evidence to suggest that ticks are unable to disperse shorter distances at sea with their hosts, with no pattern of population differentiation found among Australian sites. Our results suggest that terrestrial seabird parasites may be quite capable of short-distance movements, but only sporadic longer-distance (trans-oceanic) dispersal.     

Geographical variation,population structure and gene flow between populations of Chrysophtharta agricola (Coleoptera: Chrysomelidae), a pest of Australian eucalypt plantations

Chrysophtharta agricola (Chapuis) is a pest of commercial eucalypt plantations in Tasmania and Victoria. Vagility of pest populations may result in difficulty predicting temporal and spatial pest outbreaks, and influence genetic resistance to chemical control. Gene flow in this pest species was estimated to assess predictability of attack, the potential efficacy of natural enemies, and the likelihood of resistance build-up. Ten geographic populations of C. agricola (six from Tasmania, one from the Australian Capital Territory, one from New South Wales and two from Victoria) were examined for genetic variation and gene flow using cellulose acetate allozyme electrophoresis. Six enzyme systems (PGI, PGD, PGM, IDH, HEX and MPI) were consistently polymorphic and scorable and were used to quantify estimated gene flow between populations. FST values and analysis of molecular variance indicated that gene flow was restricted between populations. Chrysophtharta agricola exhibited high levels of heterozygosity, probably because of high allelic diversity, and because all loci examined were polymorphic. The southern-most population was the most genetically different to other Tasmanian populations, and may also have been the most recently colonized. Limited gene flow implies that outbreaks of C. agricola should be spatially predictable and populations susceptible to control by natural enemies. Our results also imply that genetic resistance to chemical control may occur under frequent application of insecticide. However, testing population movement between plantations and native forest also needs to be conducted to assess gene flow between forest types.     

Patterns of mitochondrial haplotype diversity in the invasive pest Epiphyas postvittana (Lepidoptera: Tortricidae)

The light brown apple moth, Epiphyas postvittana (Walker) (Lepidoptera: Tortricidae), is a horticultural pest of Australia and New Zealand that has more recently invaded Hawaii, Europe, and California. A 2,216-bp region of the mitochondrial genome containing the cytochrome oxidase I and II genes was sequenced from 752 individuals. Haplotype network analyses revealed a major split between a predominantly Western Australian clade and all other samples, suggestive of either a deep genetic divergence or a cryptic species. Nucleotide and haplotype diversity were highest in the country of origin, Australia, and in New Zealand populations, with evidence of haplotype sharing between New Zealand and Tasmania. Nucleotide and haplotype diversity were higher in California than within the British Isles or Hawaii. From the total of 96 haplotypes, seven were found in California, of which four were private. Within California, there have been at least two introductions; based on genetic diversity we were unable to assign a likely source for a single moth found and eradicated in Los Angeles in 2007; however, our data suggest it is unlikely that Hawaii and the British Isles are sources of the major E. postvittana population found throughout the rest of the state since 2006.     

Identifying eradication units in an invasive mammalian pest species

Eradication of invasive species is an important component for species conservation and ecosystem restoration. Success of eradications is dependent on knowledge of population connectivity in order to determine reinvasion pathways, and hence populations requiring simultaneous eradication (eradication units). The common brushtail possum (Trichosurus vulpecula) was introduced into New Zealand from Australia and Tasmania, and now occupies a wide range of habitats across the majority of New Zealand. Possums are one of the most destructive invaders within New Zealand, with extensive control operations occurring throughout the country. Understanding the population connectivity of possums on mainland New Zealand (North and South Islands) will enhance the success of planned eradications. We examined the genetic population structure of invasive possums to identify gene flow, thus reinvasion pathways, between seven populations around Dunedin and on the Otago Peninsula where eradication of possums is occurring. Genetic variation at 12 microsatellite loci was comparable between all sampled populations and exhibited a significant isolation by distance pattern. Bayesian clustering methods supported the existence of two population clusters, indicating the presence of a reinvasion pathway onto the Otago Peninsula from urban areas at the Southern end of the Peninsula. To avoid recolonisation, all possums on the Otago Peninsula should be eradicated simultaneously, with the implementation and ongoing maintenance and monitoring of an urban buffer zone. We recommend pre-eradication genetic analyses be adopted by all pest managers to define appropriate eradication units, thereby maximising eradication success and avoiding costly failures.     

Microsatellites isolated from diamondback moth,Plutella xylostella (L.), for studies of dispersal in Australian populations

Diamondback moth, Plutella xylostella (L.), is a worldwide agricultural pest that has developed resistance to many insecticides used for its control. Population structure and gene flow are yet to be determined for P. xylostella in Australia, but are important factors for the design of effective control strategies. We have isolated six polymorphic microsatellite markers: three from a partial genomic library, two from an Expressed Sequence Tagged library and one from an aminopeptidase intron of P. xylostella. These microsatellites will be used to determine population structure and gene flow in Australian populations of P. xylostella to improve insecticide resistance management.     

Independent histories of human Y chromosomes from Melanesia and Australia

To investigate the origins and relationships of Australian and Melanesian populations, 611 males from 18 populations from Australia, Melanesia, and eastern/southeastern Asia were typed for eight single-nucleotide polymorphism (SNP) loci and seven short tandem-repeat loci on the Y chromosome. A unique haplotype, DYS390.1del/RPS4Y711T, was found at a frequency of 53%-69% in Australian populations, whereas the major haplotypes found in Melanesian populations (M4G/M5T/M9G and DYS390.3del/RPS4Y711T) are absent from the Australian populations. The Y-chromosome data thus indicate independent histories for Australians and Melanesians, a finding that is in agreement with evidence from mtDNA but that contradicts some analyses of autosomal loci, which show a close relationship between Australian and Melanesian (specifically, highland Papua New Guinean) populations. Since the Australian and New Guinean landmasses were connected when first colonized by humans > or =50,000 years ago but separated some 8,000 years ago, a possible way to reconcile all the genetic data is to infer that the Y-chromosome and mtDNA results reflect the past 8,000 years of independent history for Australia and New Guinea, whereas the autosomal loci reflect the long preceding period of common origin and shared history. Two Y-chromosome haplotypes (M119C/M9G and M122C/M9G) that originated in eastern/southeastern Asia are present in coastal and island Melanesia but are rare or absent in both Australia and highland Papua New Guinea. This distribution, along with demographic analyses indicating that population expansions for both haplotypes began approximately 4,000-6,000 years ago, suggests that these haplotypes were brought to Melanesia by the Austronesian expansion. Most of the populations in this study were previously typed for mtDNA SNPs; population differentiation is greater for the Y chromosome than for mtDNA and is significantly correlated with geographic distance, a finding in agreement with results of similar analyses of European populations.     

The effects on genetic variability following a recent colonization event: the Australian sheep blowfly, Lucilia cuprina arrives in New Zealand

The Australian sheep blowfly, Lucilia cuprina, was first identified in New Zealand in 1988 and is now found to have spread throughout many sheep-farming regions. L. cuprina is estimated to have been present in New Zealand < 20 years, while in Australia L. cuprina has been estimated present > 100 years. The aim of this study was to determine the genetic effects of colonization of L. cuprina and to compare populations of L. cuprina from these two countries in terms of genetic variability and differentiation. Allozyme electrophoresis was used which revealed variability at six loci. 1680 blowflies were examined from 56 sites throughout L. cuprina's range in both countries. Genetic variability at each locus in terms of allele composition was found to be high and genetic differentiation varied considerably in New Zealand in comparison to Australia. Temporal sampling in New Zealand suggests seasonal fluctuations of population size in the recently colonized region of the South Island.     

Few genetic differences between Victorian and Western Australian blue penguins,Eudyptula minor

Abstract

Blue penguins, Eudyptula minor, breeding on Penguin Island, Western Australia are considerably larger than other blue penguins in Australia. If genetic isolation is the cause, it may have implications for the conservation status of some blue penguin populations. We compared the sequences of two mitochondrial gene regions (cytochrome‐b and the control region) from Western Australian blue penguins with other populations of blue penguins from Australia and New Zealand. We found few differences between sequences from Western Australia, Phillip Island, Victoria and Otago, New Zealand, although all three differed considerably from other New Zealand blue penguins. Sequences for the control region from the Western Australian blue penguins and 30 more birds breeding at various Australasian sites provided further support for two major clades within Eudyptula; an Australian clade (including Otago) and a New Zealand clade.     

Reconstructing recent divergence: evaluating nonequilibrium population structure in New Zealand chinook salmon

Newly established or perturbed populations are often the focus of conservation concerns but they pose special challenges for population genetics because drift?migration equilibrium is unlikely. To advance our understanding of the evolution of such populations, we investigated structure and gene flow among populations of chinook salmon that formed via natural straying following introduction to New Zealand in the early 1900s. We examined 11 microsatellite loci from samples collected in several sites and years to address two questions: (i) what population differentiation has arisen in the ≈ 30 generations since salmon were introduced to New Zealand, relative to temporal variation within populations; and (ii) what are the approximate effective population sizes and amounts of gene flow in these populations? These questions are routinely addressed in studies of indigenous populations, but less often in the case of new populations and rarely with consideration of equilibrium assumptions. We show that despite the recent introduction, continued gene flow and high temporal variability among samples, detectable population structure has arisen among the New Zealand populations, consistent with their colonization pattern and isolation by geographical distance. Furthermore, we use simple individual‐based simulations and estimates of effective population sizes to estimate the effective gene flow among drainages under likely nonequilibrium conditions. Similar methodology may be broadly applicable to other studies of population structure and phenotypic evolution under similar nonequilibrium, high gene flow conditions.     

Phylogeographic analysis of the brown alga Cutleria multifida (Tilopteridales,Phaeophyceae) suggests a complicated introduction history

In depth genetic comparisons of populations of Cutleria multifida (Tilopteridales, Phaeophyceae) collected from Europe, the northwestern Pacific Ocean, Australia and New Zealand using the DNA sequences of four gene regions (the mitochondrial cox2 and cox3 genes, the intergeneric spacer region adjacent to cox3, and the open reading frame) suggested that the northwestern European and Japanese populations were considerably greater in terms of their genetic divergence than Mediterranean, Australian or New Zealand populations. The haplotypes of the populations in northwestern European (distribution range including the type locality, seven haplotypes) and Japanese populations (seven haplotypes) were unique except for one shared haplotype. There were weak but positive correlations between the geographical distance and the genetic divergence among northwestern European and Japanese populations. Moreover, both female and male gametophytes occurred in eight of the nine Japanese localities, suggesting Japanese populations showed normal sexual heteromorphic life history of the species. In light of these results, it appears that Japanese populations were native to the area despite earlier hypothesis. In contrast, Australian and New Zealand populations were composed of only one haplotype that is very close to those found in northwestern Europe and Japan, suggesting a recent introduction history from Europe (or from northeastern Asia via Europe) by ship transport to Australia and New Zealand. The Mediterranean populations included two haplotypes identical to those found in northwestern Europe and Japan, and it is suggestive of transoceanic introductions of some populations between Mediterranean and Japanese coasts.     

Population genetic structure of Cydia pomonella: a review and case study comparing spatiotemporal variation

Analysis of population genetic structure is a key aspect to understand insect pest population dynamics in agricultural scenarios. Here the role of geography, hosts and time on the population genetic structure of codling moth Cydia pomonella (Linnaeus) (Lep., Tortricidae) populations is described. Temporal variation was examined in two French orchards among each of three adult flights during two successive years. Analyses were conducted using two insecticide resistance markers (variation at the sodium channel gene and enzymatic activity of cytochrome P450 oxidases) and three microsatellite loci. Levels of genetic variation among temporal populations were not significant based on variation in the sodium channel gene and microsatellite loci. However, P450 oxidase activity differed significantly during both flights and years, decreasing during the three flights of the first year and increasing during the second. These results suggest that phytosanitary measures are among the factors shaping the genetic structure of C. pomonella populations over temporal and geographical scales. We discuss the relative importance of natural and passive dispersal related to anthropogenic activities affecting C. pomonella population genetics and highlight population genetic research needs in order to design more efficient pest management practices.     

Experimental analysis of the influence of pest management practice on the efficacy of an endemic arthropod natural enemy complex of the diamondback moth

Maximizing the contribution of endemic natural enemies to integrated pest management (IPM), programs requires a detailed knowledge of their interactions with the target pest. This experimental field study evaluated the impact of the endemic natural enemy complex of Plutella xylostella (L.) (Lepidoptera: Yponomeutidae) on pest populations in commercial cabbage crops in southeastern Queensland, Australia. Management data were used to score pest management practices at experimental sites on independent Brassica farms practicing a range of pest management strategies, and mechanical methods of natural enemy exclusion were used to assess the impact of natural enemies on introduced cohorts of P. xylostella at each site. Natural enemy impact was greatest at sites adopting IPM and least at sites practicing conventional pest management strategies. At IPM sites, the contribution of natural enemies to P. xylostella mortality permitted the cultivation of marketable crops with no yield loss but with a substantial reduction in insecticide inputs. Three species of larval parasitoids (Diadegma semiclausum Hellén [Hymenoptera: Ichneumonidae], Apanteles ippeus Nixon [Hymenoptera: Braconidae], and Oomyzus sokolowskii Kurdjumov [Hymenoptera: Eulophidae]) and one species of pupal parasitoid Diadromus collaris Gravenhorst (Hymenoptera: Ichneumonidae) attacked immature P. xylostella. The most abundant groups of predatory arthropods caught in pitfall traps were Araneae (Lycosidae) > Coleoptera (Carabidae, Coccinelidae, Staphylinidae) > Neuroptera (Chrysopidae) > Formicidae, whereas on crop foliage Araneae (Clubionidae, Oxyopidae) > Coleoptera (Coccinelidae) > Neuroptera (Chrysopidae) were most common. The abundance and diversity of natural enemies was greatest at sites that adopted IPM, correlating greater P. xylostella mortality at these sites. The efficacy of the natural enemy complex to pest mortality under different pest management regimes and appropriate strategies to optimize this important natural resource are discussed.     

Spatial subdivision among assemblages of Spanish mackerel, Scomberomorus commerson (Pisces: Scombridae) across northern Australia: implications for fisheries management

Aim This study investigated the use of stable δ¹³C and δ¹⁸O isotopes in the sagittal otolith carbonate of narrow‐barred Spanish mackerel, Scomberomorus commerson, as indicators of population structure across Australia. Location Samples were collected from 25 locations extending from the lower west coast of Western Australia (30°), across northern Australian waters, and to the east coast of Australia (18°) covering a coastline length of approximately 9500 km, including samples from Indonesia. Methods The stable δ¹³C and δ¹⁸O isotopes in the sagittal otolith carbonate of S. commerson were analysed using standard mass spectrometric techniques. The isotope ratios across northern Australian subregions were subjected to an agglomerative hierarchical cluster analysis to define subregions. Isotope ratios within each of the subregions were compared to assess population structure across Australia. Results Cluster analysis separated samples into four subregions: central Western Australia, north Western Australia, northern Australia and the Gulf of Carpentaria and eastern Australia. Isotope signatures for fish from a number of sampling sites from across Australia and Indonesia were significantly different, indicating population separation. No significant differences were found in otolith isotope ratios between sampling times (no temporal variation). Main conclusions Significant differences in the isotopic signatures of S. commerson demonstrate that there is unlikely to be any substantial movement of fish among these spatially discrete adult assemblages. The lack of temporal variation among otolith isotope ratios indicates that S. commerson populations do not undergo longshore spatial shifts in distribution during their life history. The temporal persistence of spatially explicit stable isotopic signatures indicates that, at these spatial scales, the population units sampled comprise functionally distinct management units or separate ‘stocks’ for many of the purposes of fisheries management. The spatial subdivision evident among populations of S. commerson across northern and western Australia indicates that it may be advantageous to consider S. commerson population dynamics and fisheries management from a metapopulation perspective (at least at the regional level).