Unraveling the habitat preferences of two closely related bumble bee species in Eastern Europe

Co‐occurrence of closely related species is often explained through resource partitioning, where key morphological or life‐history traits evolve under strong divergent selection. In bumble bees (genus Bombus), differences in tongue lengths, nest sites, and several life‐history traits are the principal factors in resource partitioning. However, the buff‐tailed and white‐tailed bumble bee (Bombus terrestris and B. lucorum respectively) are very similar in morphology and life history, but their ranges nevertheless partly overlap, raising the question how they are ecologically divergent. What little is known about the environmental factors determining their distributions stems from studies in Central and Western Europe, but even less information is available about their distributions in Eastern Europe, where different subspecies occur. Here, we aimed to disentangle the broad habitat requirements and associated distributions of these species in Romania and Bulgaria. First, we genetically identified sampled individuals from many sites across the study area. We then not only computed species distributions based on presence‐only data, but also expanded on these models using relative abundance data. We found that B. terrestris is a more generalist species than previously thought, but that B. lucorum is restricted to forested areas with colder and wetter climates, which in our study area are primarily found at higher elevations. Both vegetation parameters such as annual mean Leaf Area Index and canopy height, as well as climatic conditions, were important in explaining their distributions. Although our models based on presence‐only data suggest a large overlap in their respective distributions, results on their relative abundance suggest that the two species replace one another across an environmental gradient correlated to elevation. The inclusion of abundance enhances our understanding of the distribution of these species, supporting the emerging recognition of the importance of abundance data in species distribution modeling.     

Targeted agri‐environment schemes significantly improve the population size of common farmland bumblebee species

Changes in agricultural practice across Europe and North America have been associated with range contractions and local extinction of bumblebees (Bombus spp.). A number of agri‐environment schemes have been implemented to halt and reverse these declines, predominantly revolving around the provision of additional forage plants. Although it has been demonstrated that these schemes can attract substantial numbers of foraging bumblebees, it remains unclear to what extent they actually increase bumblebee populations. We used standardized transect walks and molecular techniques to compare the size of bumblebee populations between Higher Level Stewardship (HLS) farms implementing pollinator‐friendly schemes and Entry Level Stewardship (ELS) control farms. Bumblebee abundance on the transect walks was significantly higher on HLS farms than ELS farms. Molecular analysis suggested maximum foraging ranges of 566 m for Bombus hortorum, 714 m for B. lapidarius, 363 m for B. pascuorum and 799 m for B. terrestris. Substantial differences in maximum foraging range were found within bumblebee species between farm types. Accounting for foraging range differences, B. hortorum (47 vs 13 nests/km²) and B. lapidarius (45 vs 22 nests/km²) were found to nest at significantly greater densities on HLS farms than ELS farms. There were no significant differences between farm type for B. terrestris (88 vs 38 nests/km²) and B. pascuorum (32 vs 39 nests/km²). Across all bumblebee species, HLS management had a significantly positive effect on bumblebee nest density. These results show that targeted agri‐environment schemes that increase the availability of suitable forage can significantly increase the size of wild bumblebee populations.     

Volatility in the effective size of a freshwater gastropod population

Despite the utility of gastropod models for the study of evolutionary processes of great generality and importance, their effective population size has rarely been estimated in the field. Here, we report allele frequency variance at three allozyme‐encoding loci monitored over 7 years in a population of the invasive freshwater pulmonate snail Physa acuta (Draparnaud 1805), estimating effective population size with both single‐sample and two‐sample approaches. Estimated N_e declined from effectively infinite in 2009 to approximately 40–50 in 2012 and then rose back to infinity in 2015, corresponding to a striking fluctuation in the apparent census size of the population. Such volatility in N_e may reflect cryptic population subdivision.     

Food quality affects the expression of antimicrobial peptide genes upon simulated parasite attack in the larvae of greater wax moth

Predator‐prey interactions are an important evolutionary force affecting the immunity of the prey. Parasitoids and mites pierce the cuticle of their prey, which respond by activating their immune system against predatory attacks. Immunity is a costly function for the organism, as it often competes with other life‐history traits for limited nutrients. We tested whether the expression of antimicrobial peptides (AMP) of the larvae of the greater wax moth Galleria mellonella (L.) (Lepidoptera: Pyralidae) changes as a consequence of insertion of a nylon monofilament, which acts like a synthetic parasite. The treatment was done for larvae grown on a high‐quality vs. a low‐quality diet. The expression of Gloverin and 6‐tox were upregulated in response to the insertion of the nylon monofilament. The expression of 6‐tox, Cecropin‐D, and Gallerimycin were significantly higher in the ‘low‐quality diet’ group than in the ‘high‐quality diet’ group. As food quality seems to affect AMP gene expression in G. mellonella larvae, it should always be controlled for in studies on bacterial and fungal infections in G. mellonella.     

Relationships between intra‐specific variation in seed size and recruitment in four species in two contrasting habitats

Large seeds contain more stored resources, and seedlings germinating from large seeds generally cope better with environmental stresses such as shading, competition and thick litter layers, than seedlings germinating from small seeds. A pattern with small‐seeded species being associated with open habitats and large‐seeded species being associated with closed (shaded) habitats has been suggested and supported by comparative studies. However, few studies have assessed the intra‐specific relationship between seed size and recruitment, comparing plant communities differing in canopy cover. Here, seeds from four plant species commonly occurring in ecotones between open and closed habitats (Convallaria majalis, Frangula alnus, Prunus padus and Prunus spinosa) were weighed and sown individually (3200 seeds per species) in open and closed‐canopy sites, and seedling emergence and survival recorded over 3 years. Our results show a generally positive, albeit weak, relationship between seed size and recruitment. In only one of the species, C. majalis, was there an association between closed canopy habitat and a positive seed size effect on recruitment. We conclude that there is a weak selection gradient favouring larger seeds, but that this selection gradient is not clearly related to habitat.     

Extremely low effective population sizes, genetic structuring and reduced genetic diversity in a threatened bumblebee species, Bombus sylvarum (Hymenoptera: Apidae)

Habitat fragmentation may severely affect survival of social insect populations as the number of nests per population, not the number of individuals, represents population size, hence they may be particularly prone to loss of genetic diversity. Erosion of genetic diversity may be particularly significant among social Hymenoptera such as bumblebees (Bombus spp.), as this group may be susceptible to diploid male production, a suggested direct cost of inbreeding. Here, for the first time, we assess genetic diversity and population structuring of a threatened bumblebee species (Bombus sylvarum) which exists in highly fragmented habitat (rather than oceanic) islands. Effective population sizes, estimated from identified sisterhoods, were very low (range 21-72) suggesting that isolated populations will be vulnerable to loss of genetic variation through drift. Evidence of significant genetic structuring between populations (theta = 0.084) was found, but evidence of a bottleneck was detected in only one population. Comparison across highly fragmented UK populations and a continental population (where this species is more widespread) revealed significant differences in allelic richness attributable to a high degree of genetic diversity in the continental population. While not directly related to population size, this is perhaps explained by the high degree of isolation between UK populations relative to continental populations. We suggest that populations now existing on isolated habitat islands were probably linked by stepping-stone populations prior to recent habitat loss.     

The genetic effective and adult census size of an Australian population of tiger prawns (Penaeus esculentus)

This study compares estimates of the census size of the spawning population with genetic estimates of effective current and long-term population size for an abundant and commercially important marine invertebrate, the brown tiger prawn (Penaeus esculentus). Our aim was to focus on the relationship between genetic effective and census size that may provide a source of information for viability analyses of naturally occurring populations. Samples were taken in 2001, 2002 and 2003 from a population on the east coast of Australia and temporal allelic variation was measured at eight polymorphic microsatellite loci. Moments-based and maximum-likelihood estimates of current genetic effective population size ranged from 797 to 1304. The mean long-term genetic effective population size was 9968. Although small for a large population, the effective population size estimates were above the threshold where genetic diversity is lost at neutral alleles through drift or inbreeding. Simulation studies correctly predicted that under these experimental conditions the genetic estimates would have non-infinite upper confidence limits and revealed they might be overestimates of the true size. We also show that estimates of mortality and variance in family size may be derived from data on average fecundity, current genetic effective and census spawning population size, assuming effective population size is equivalent to the number of breeders. This work confirms that it is feasible to obtain accurate estimates of current genetic effective population size for abundant Type III species using existing genetic marker technology.     

Pollen limitation and Allee effect related to population size and sex ratio in the endangered Ottelia acuminata (Hydrocharitaceae): implications for conservation and reintroduction

Small populations may suffer more severe pollen limitation and result in Allee effects. Sex ratio may also affect pollination and reproduction success in dioecious species, which is always overlooked when performing conservation and reintroduction tasks. In this study, we investigated whether and how population size and sex ratio affected pollen limitation and reproduction in the endangered Ottelia acuminata, a dioecious submerged species. We established experimental plots with increasing population size and male sex ratio. We observed insect visitation, estimated pollen limitation by hand‐pollinations and counted fruit set and seed production per fruit. Fruit set and seed production decreased significantly in small populations due to pollinator scarcity and thus suffered more severe pollen limitation. Although frequently visited, female‐biased larger populations also suffered severe pollen limitation due to few effective visits and insufficient pollen availability. Rising male ratio enhanced pollination service and hence reproduction. Unexpectedly, pollinator preferences did not cause reduced reproduction in male‐biased populations because of high pollen availability. However, reproductive outputs showed more variability in severe male‐biased populations. Our results revealed two component Allee effects in fruit set and seed production, mediated by pollen limitation in O. acuminata. Moreover, reproduction decreased significantly in larger female‐biased populations, increasing the risk of an Allee effect.     

The population genomic basis of geographic differentiation in North American common ragweed (Ambrosia artemisiifolia L.)

Common ragweed (Ambrosia artemisiifolia L.) is an invasive, wind‐pollinated plant nearly ubiquitous in disturbed sites in its eastern North American native range and present across growing portions of Europe, Africa, Asia, and Australia. Phenotypic divergence between European and native‐range populations has been described as rapid evolution. However, a recent study demonstrated major human‐mediated shifts in ragweed genetic structure before introduction to Europe and suggested that native‐range genetic structure and local adaptation might fully explain accelerated growth and other invasive characteristics of introduced populations. Genomic differentiation that potentially influenced this structure has not yet been investigated, and it remains unclear whether substantial admixture during historical disturbance of the native range contributed to the development of invasiveness in introduced European ragweed populations. To investigate fine‐scale population genetic structure across the species' native range, we characterized diallelic SNP loci via a reduced‐representation genotyping‐by‐sequencing (GBS) approach. We corroborate phylogeographic domains previously discovered using traditional sequencing methods, while demonstrating increased power to resolve weak genetic structure in this highly admixed plant species. By identifying exome polymorphisms underlying genetic differentiation, we suggest that geographic differentiation of this important invasive species has occurred more often within pathways that regulate growth and response to defense and stress, which may be associated with survival in North America's diverse climatic regions.     

Contemporary effective population and metapopulation size (Ne and meta‐Ne): comparison among three salmonids inhabiting a fragmented system and differing in gene flow and its asymmetries

We estimated local and metapopulation effective sizes ( and meta‐) for three coexisting salmonid species (Salmo salar, Salvelinus fontinalis, Salvelinus alpinus) inhabiting a freshwater system comprising seven interconnected lakes. First, we hypothesized that might be inversely related to within‐species population divergence as reported in an earlier study (i.e., F_ST: S. salar> S. fontinalis> S. alpinus). Using the approximate Bayesian computation method implemented in ONeSAMP, we found significant differences in () between species, consistent with a hierarchy of adult population sizes (). Using another method based on a measure of linkage disequilibrium (LDNE: ), we found more finite values for S. salar than for the other two salmonids, in line with the results above that indicate that S. salar exhibits the lowest among the three species. Considering subpopulations as open to migration (i.e., removing putative immigrants) led to only marginal and non‐significant changes in , suggesting that migration may be at equilibrium between genetically similar sources. Second, we hypothesized that meta‐ might be significantly smaller than the sum of local s (null model) if gene flow is asymmetric, varies among subpopulations, and is driven by common landscape features such as waterfalls. One ‘bottom‐up’ or numerical approach that explicitly incorporates variable and asymmetric migration rates showed this very pattern, while a number of analytical models provided meta‐ estimates that were not significantly different from the null model or from each other. Our study of three species inhabiting a shared environment highlights the importance and utility of differentiating species‐specific and landscape effects, not only on dispersal but also in the demography of wild populations as assessed through local s and meta‐s and their relevance in ecology, evolution and conservation.     

TaDA1, a conserved negative regulator of kernel size,has an additive effect with TaGW2 in common wheat (Triticum aestivum L.)

Kernel size is an important trait determining cereal yields. In this study, we cloned and characterized TaDA1, a conserved negative regulator of kernel size in wheat (Triticum aestivum). The overexpression of TaDA1 decreased the size and weight of wheat kernels, while its down‐regulation using RNA interference (RNAi) had the opposite effect. Three TaDA1‐A haplotypes were identified in Chinese wheat core collections, and a haplotype association analysis showed that TaDA1‐A‐HapI was significantly correlated with the production of larger kernels and higher kernel weights in modern Chinese cultivars. The haplotype effect resulted from a difference in TaDA1‐A expression levels between genotypes, with TaDA1‐A‐HapI resulting in lower TaDA1‐A expression levels. This favourable haplotype was found having been positively selected during wheat breeding over the last century. Furthermore, we demonstrated that TaDA1‐A physically interacts with TaGW2‐B. The additive effects of TaDA1‐A and TaGW2‐B on kernel weight were confirmed not only by the phenotypic enhancement arising from the simultaneous down‐regulation of TaDA1 and TaGW2 expression, but also by the combinational haplotype effects estimated from multi‐environment field data from 348 wheat cultivars. A comparative proteome analysis of developing transgenic and wild‐type grains indicated that TaDA1 and TaGW2 are involved in partially overlapping but relatively independent protein regulatory networks. Thus, we have identified an important gene controlling kernel size in wheat and determined its interaction with other genes regulating kernel weight, which could have beneficial applications in wheat breeding.     

Brucin,an antibacterial peptide derived from fruit protein of Fructus Bruceae,Brucea javanica (L.) Merr

A novel antibacterial peptide specific to Streptococcus pyogenes was produced from dried fruit protein of Brucea javanica (L.) Merr. A mixture of active peptides from the fruit protein was produced in vitro by pepsin hydrolysis. The hydrolysate was purified by reverse‐phase HPLC, and antimicrobial peptides active against Gram‐negative and Gram‐positive bacteria were analysed using SDS‐PAGE and nanoLC‐MS/MS. Here, four possible peptides were obtained and chemically synthesized for comparative study of the growth inhibition of Strep. pyogenes. One chemically synthesized peptide with a molecular mass of 1168·31 Da, His‐Thr‐Leu‐Cys‐Met‐Asp‐Gly‐Gly‐Ala‐Thr‐Tyr, showed the most potent antibacterial activity against Strep. pyogenes. This 11‐amino acid peptide was named Brucin. Its bacterial inhibitory activity was 16‐fold and 12·5‐fold higher than penicillin G and chloramphenicol, respectively, with a MIC value of 20 μmol l^?1. The results suggest that Brucin, a potent antibiotic peptide, may be developed as an alternative drug for the treatment of the disease caused by Strep. pyogenes.

Significance and Impact of the Study

An antibacterial peptide, named Brucin with specificity for Streptococcus pyogenes, was produced in vitro from dried fruit protein of Brucea javanica (L.) Merr. by pepsin‐catalysed hydrolysis. Its inhibitory activity towards the Gram‐positive bacteria was higher than penicillin G and chloramphenicol. The result suggested that Brucin may be applied for the treatment of the disease caused by Strep. pyogenes^*.     

RAD sequencing reveals genomewide divergence between independent invasions of the European green crab (Carcinus maenas) in the Northwest Atlantic

Genomic studies of invasive species can reveal both invasive pathways and functional differences underpinning patterns of colonization success. The European green crab (Carcinus maenas) was initially introduced to eastern North America nearly 200 years ago where it expanded northwards to eastern Nova Scotia. A subsequent invasion to Nova Scotia from a northern European source allowed further range expansion, providing a unique opportunity to study the invasion genomics of a species with multiple invasions. Here, we use restriction‐site‐associated DNA sequencing‐derived SNPs to explore fine‐scale genomewide differentiation between these two invasions. We identified 9137 loci from green crab sampled from 11 locations along eastern North America and compared spatial variation to mitochondrial COI sequence variation used previously to characterize these invasions. Overall spatial divergence among invasions was high (pairwise F_ST ~0.001 to 0.15) and spread across many loci, with a mean F_ST ~0.052 and 52% of loci examined characterized by F_ST values >0.05. The majority of the most divergent loci (i.e., outliers, ~1.2%) displayed latitudinal clines in allele frequency highlighting extensive genomic divergence among the invasions. Discriminant analysis of principal components (both neutral and outlier loci) clearly resolved the two invasions spatially and was highly correlated with mitochondrial divergence. Our results reveal extensive cryptic intraspecific genomic diversity associated with differing patterns of colonization success and demonstrates clear utility for genomic approaches to delineating the distribution and colonization success of aquatic invasive species.     

Spatial structuring of the population genetics of a European subterranean termite species

In population genetics studies, detecting and quantifying the distribution of genetic variation can help elucidate ecological and evolutionary processes. In social insects, the distribution of population‐level genetic variability is generally linked to colony‐level genetic structure. It is thus especially crucial to conduct complementary analyses on such organisms to examine how spatial and social constraints interact to shape patterns of intraspecific diversity. In this study, we sequenced the mitochondrial COII gene for 52 colonies of the subterranean termite Reticulitermes grassei (Isoptera: Rhinotermitidae), sampled from a population in southwestern France. Three haplotypes were detected, one of which was found exclusively in the southern part of the study area (near the Pyrenees). After genotyping 6 microsatellite loci for 512 individual termites, we detected a significant degree of isolation by distance among individuals over the entire range; however, the cline of genetic differentiation was not continuous, suggesting the existence of differentiated populations. A spatial principal component analysis based on allele frequency data revealed significant spatial autocorrelation among genotypes: the northern and southern groups were strongly differentiated. This finding was corroborated by clustering analyses; depending on the randomized data set, two or three clusters, exhibiting significant degrees of differentiation, were identified. An examination of colony breeding systems showed that colonies containing related neotenic reproductives were prevalent, suggesting that inbreeding may contribute to the high level of homozygosity observed and thus enhance genetic contrasts among colonies. We discuss the effect of evolutionary and environmental factors as well as reproductive and dispersal modes on population genetic structure.     

Post‐fragmentation population structure in a cooperative breeding Afrotropical cloud forest bird: emergence of a source‐sink population network

The impact of demographic parameters on the genetic population structure and viability of organisms is a long‐standing issue in the study of fragmented populations. Demographic and genetic tools are now readily available to estimate census and effective population sizes and migration and gene flow rates with increasing precision. Here we analysed the demography and genetic population structure over a recent 15‐year time span in five remnant populations of Cabanis's greenbul (Phyllastrephus cabanisi), a cooperative breeding bird in a severely fragmented cloud forest habitat. Contrary to our expectation, genetic admixture and effective population sizes slightly increased, rather than decreased between our two sampling periods. In spite of small effective population sizes in tiny forest remnants, none of the populations showed evidence of a recent population bottleneck. Approximate Bayesian modelling, however, suggested that differentiation of the populations coincided at least partially with an episode of habitat fragmentation. The ratio of meta‐N_e to meta‐N_c was relatively low for birds, which is expected for cooperative breeding species, while N_e/N_c ratios strongly varied among local populations. While the overall trend of increasing population sizes and genetic admixture may suggest that Cabanis's greenbuls increasingly cope with fragmentation, the time period over which these trends were documented is rather short relative to the average longevity of tropical species. Furthermore, the critically low N_c in the small forest remnants keep the species prone to demographic and environmental stochasticity, and it remains open if, and to what extent, its cooperative breeding behaviour helps to buffer such effects.     

Successional loss of two key food tree species best explains decline in group size of Panamanian howler monkeys (Alouatta palliata)

Negative impacts of discrete, short‐term disturbances to wildlife populations are well‐documented. The consequences of more gradual environmental change are less apparent and harder to study because they play out over longer periods and are often indirect in their action. Yet, they can drive the decline of wildlife populations even in seemingly pristine and currently well‐protected habitats. One such environmental change is a successional shift in a community's species composition as it regenerates from disturbance caused by past human land use. Early and middle successional tree species often provide key foods to folivores and frugivores, but the abundance of these resources drops as the forest matures, with adverse repercussions for these consumers. Our 44‐year record (1974–2018) of howler monkey (Alouatta palliata) group sizes and demographic composition from Barro Colorado Island, Panama, a protected reserve, documents an example of this phenomenon. After 70 years of relative stability, the mean size of howler monkey groups exhibited a marked decline, beginning in 2003. This downward trajectory in group size has continued through the most recent census in 2018. The composition of howler groups also changed significantly during the study period, with the patterns of decline differing among age/sex classes. There is no evidence that these changes were caused by increased rates of emigration, group fission, predation, parasitism, or disease. Rather, they are best explained by an island‐wide, succession‐driven decline in the densities of two species of free‐standing fig trees, Ficus yoponensis and F. insipida, which together were providing ~36% of BCI howlers’ annual diet. Abstract in Spanish is available with online material.     

Genomic data detect corresponding signatures of population size change on an ecological time scale in two salamander species

Understanding the demography of species over recent history (e.g. <100 years) is critical in studies of ecology and evolution, but records of population history are rarely available. Surveying genetic variation is a potential alternative to census‐based estimates of population size, and can yield insight into the demography of a population. However, to assess the performance of genetic methods, it is important to compare their estimates of population history to known demography. Here, we leveraged the exceptional resources from a wetland with 37 years of amphibian mark–recapture data to study the utility of genetically based demographic inference on salamander species with documented population declines (Ambystoma talpoideum) and expansions (A. opacum), patterns that have been shown to be correlated with changes in wetland hydroperiod. We generated ddRAD data from two temporally sampled populations of A. opacum (1993, 2013) and A. talpoideum (1984, 2011) and used coalescent‐based demographic inference to compare alternate evolutionary models. For both species, demographic model inference supported population size changes that corroborated mark–recapture data. Parameter estimation in A. talpoideum was robust to our variations in analytical approach, while estimates for A. opacum were highly inconsistent, tempering our confidence in detecting a demographic trend in this species. Overall, our robust results in A. talpoideum suggest that genome‐based demographic inference has utility on an ecological scale, but researchers should also be cognizant that these methods may not work in all systems and evolutionary scenarios. Demographic inference may be an important tool for population monitoring and conservation management planning.     

The solution structure of horseshoe crab antimicrobial peptide tachystatin B with an inhibitory cystine-knot motif.

Tachystatin B is an antimicrobial and a chitin-binding peptide isolated from the Japanese horseshoe crab (Tachypleus tridentatus) consisting of two isopeptides called tachystatin B1 and B2. We have determined their solution structures using NMR experiments and distance geometry calculations. The 20 best converged structures of tachystatin B1 and B2 exhibited root mean square deviations of 0.46 and 0.49 A, respectively, for the backbone atoms in Cys(4)-Arg(40). Both structures have identical conformations, and they contain a short antiparallel beta-sheet with an inhibitory cystine-knot (ICK) motif that is distributed widely in the antagonists for voltage-gated ion channels, although tachystatin B does not have neurotoxic activity. The structural homology search provided several peptides with structures similar to that of tachystatin B. However, most of them have the advanced functions such as insecticidal activity, suggesting that tachystatin B may be a kind of ancestor of antimicrobial peptide in the molecular evolutionary history. Tachystatin B also displays a significant structural similarity to tachystatin A, which is member of the tachystatin family. The structural comparison of both tachystatins indicated that Tyr(14) and Arg(17) in the long loop between the first and second strands might be the essential residues for binding to chitin.     

Habitat patch size alters the importance of dispersal for species diversity in an experimental freshwater community

Increased dispersal of individuals among discrete habitat patches should increase the average number of species present in each local habitat patch. However, experimental studies have found variable effects of dispersal on local species richness. Priority effects, predators, and habitat heterogeneity have been proposed as mechanisms that limit the effect of dispersal on species richness. However, the size of a habitat patch could affect how dispersal regulates the number of species able to persist. We investigated whether habitat size interacted with dispersal rate to affect the number of species present in local habitats. We hypothesized that increased dispersal rates would positively affect local species richness more in small habitats than in large habitats, because rare species would be protected from demographic extinction. To test the interaction between dispersal rate and habitat size, we factorially manipulated the size of experimental ponds and dispersal rates, using a model community of freshwater zooplankton. We found that high‐dispersal rates enhanced local species richness in small experimental ponds, but had no effect in large experimental ponds. Our results suggest that there is a trade‐off between patch connectivity (a mediator of dispersal rates) and patch size, providing context for understanding the variability observed in dispersal effects among natural communities, as well as for developing conservation and management plans in an increasingly fragmented world.