Allozyme variation in cork oak (Quercus suber L.): the role of phylogeography and genetic introgression by other Mediterranean oak species and human activities

 Genetic variation in the cork oak (Quercus suber L.) was investigated using 11 loci from seven enzyme systems in 40 populations sampled over the entire distribution of this species in the western Mediterranean Basin. Mean heterozygosity values over the polymorphic loci (Ho=0.283), the percentage of polymorphic populations (M=0.76), and the total genetic diversity (Ht=0.31) from which 11% was accounted for among-population variation, are among the highest recorded in oak species. In contrast to previous results in holm oak (Q. ilex L.), another evergreen species in the same area, cork oak possessed a smaller allele pool and a lower average number of alleles per locus and per population (A=2.0). More particularly, very few low-frequency alleles were observed in cork oak except for eight populations in which allozyme polymorphism at locus Pgi 1, diagnostic between both species, indicates that these low-frequency alleles are introgressed from holm oak. On the basis of the genetic distance estimated from allozyme frequencies, 32 of the 40 cork oak populations studied were classified into two very distinct sets which also corresponded to distinct geographic areas. One set gathered together the 18 populations from the Iberian peninsula and two adjacent areas in France, i.e. the centre of origin of cork oak, according to paleobotanical data. This set was characterized by a larger allele pool, a higher within-population genetic diversity and a lower differentiation between populations than was observed in the other set, which comprised the populations from North Africa, Sicily, Sardinia, Corsica, continental Italy and the region of Provence (southeastern France). In these more southern and eastern disjunct areas, cork oak migration from Iberia may have occurred at different periods since the end of the Tertiary. The possible effect of human activity on cork oak genetic structure, i.e. the selection of good-quality cork, acorn over-use for animal food, and even human nutrition, is discussed. Received: 3 March 1998 / Accepted: 19 March 1998     

How New Guinea rainforest flower resources vary in time and space: Implications for nectarivorous birds

Abstract Variability in spatial and temporal patterning of flowering by populations of rainforest trees fed upon by honeyeaters and flower-visiting parrots was studied for 2 years in lowland tropical hill forest in Papua New Guinea. All 2200 trees in a 3 ha plot were tagged, identified, mapped and monitored monthly. Of 274 tree species present, 86 flowered during the course of the study; during any given month, approximately 20% of the species flowering that month were visited by nectarivorous birds. Results showed that overall flower resources (total number of species, and number of bird-pollinated species, individuals and flowers) fluctuated during the year, decreased during the dry season and increased during the wet season. In addition, there was a wide range of temporal variation within and among tree species in length and timing of flowering period, percentage of each conspecific population flowering from year to year, and degree of synchrony among flowering conspecifics. Spatial dispersion of tree populations also varied, from clumps to scattered single individuals. Resident bird species were correlated with synchronously flowering trees, whereas nomadic bird species were correlated with asynchronously flowering trees. Resident birds were also associated with smaller blooming displays per tree, whereas nomadic birds were associated with trees that bloomed massively. There was no correlation between avian nomadism and spatial dispersion of tree populations. Thus nomadic birds seem to range in search of rich but unpredictable resources; resident birds may rely more on predictable, but smaller resources.     

Frequent cytoplasmic exchanges between oak species that are not closely related: Quercus suber and Q. ilex in Morocco

Chloroplast (cp) and mitochondrial (mt) DNA variation were studied in 97 populations of cork oak (Quercus suber) in Morocco; in 31 of these populations, holm oak (Quercus ilex), a clearly distinct species, also occurred and was compared with Q. suber. Three cpDNA and one mtDNA primer pairs were used in the survey, each in combination with one restriction enzyme. Six haplotypes belonging to two very divergent lineages were detected; one lineage predominates in each species, and is probably ancestral, as inferred from comparisons with other oak species. In the mixed-species populations, cytoplasmic genomes were frequently shared across species, as indicated by an introgression ratio of 0.63. This index is a new measure of the propensity of species to share locally genetic markers, varying from zero (complete differentiation) to one (no differentiation). By contrast, more closely related deciduous oak species (Q. robur, Q. petraea and Q. pubescens) have introgression ratios varying from 0.82 to 0.97. The introgression events appear to have been more frequent in the direction Q. ilex (female) x Q. suber (male), a finding which seems attributable to the flowering phenology of these two species. This asymmetry may have favoured immigration of Q. suber beyond its main range, in regions already colonized by Q. ilex. There, rare hybridization and further introgression through long distance pollen flow have established populations that are morphologically indistinguishable from Q. suber but that have cytoplasmic genomes originating from the local Q. ilex populations.     

Genetic variability of fragmented stands of pedunculate oak (Quercus robur) in Finland

The genetic structure of 33 natural Quercus robur stands in Finland was studied using 13 allozyme loci to analyze the effects of fragmentation in a wind-pollinated tree species. The present fragmented and discontinuous distribution of oak is a result of both short-term human impact and long-term climatic and geological change, including post-glacial land uplift. In accordance with general expectations, genetic diversity in small populations was lower than that in large populations, and differentiation among small populations was higher than that among large populations. Heterozygote deficiency was more pronounced in large populations, which is proposed to be a Wahlund effect created by either spatial sub-structuring or the existence of synchronized flowering lineages. Also genetic differentiation was higher and diversity lower in Finland than the estimates reported for Central Europe. There were differences in the genetic structure on sites of different geological age. We suggest that on most geologically old sites drift has a prominent effect whereas on younger sites also founder effects may be important.     

Seeing the trees for the leaves – oaks as mosaics for a host-specific moth

From the perspective of a specialist herbivore, how homogenous are individual tree crowns as patches of habitat? We partitioned variation in physical and chemical host leaf traits and in the abundance and performance of a specialist oak leaf miner, Tischeria ekebladella, into variation at different hierarchical levels. For the phenolic contents of the leaves, we examined variation among oak stands, among trees within stands and among branches within trees. For leaf size and water content, we assessed variation among trees within a single stand, among shoots within trees, and among leaves within shoots. For moth abundance and performance, we examined variation across all levels: among oak stands, among trees within stands, among branches within trees, among shoots within branches and among leaves and insect individuals. For measures of phenolic contents, we found little variation among stands but substantial variation among individual trees. Yet, a tree particularly rich in a given compound was often comparatively poor in another. At a finer spatial scale, the phenolic composition of individual parts of a single tree was quite consistent, whereas leaf weight and water content varied widely within individual tree crowns. Moth abundances varied more among shoots within branches than at any other spatial level, whereas moth survival showed equal levels of variation within individual shoots as among separate oak stands. Likewise, for four other measures of larval performance (assessed at the level of trees and lower), we found more variation within than between trees. In conclusion, the large variation observed in the performance of a specialist moth and in the physical traits of the leaves among different parts of single tree crowns refutes the image of an oak tree as an ‘island’ of internally homogeneous quality. Hence, we may expect little evolutionary adaptation of T. ekebladella at the scale of individual trees. The moths may instead evolve to behaviourally select their resource at a very fine scale. Large variation within trees also calls for extensive replication within trees in ecological sampling designs and/or the sampling of maximally similar leaves.     

Assessing the resilience of Mediterranean holm oaks to disturbances using selective thinning

Climate change will increase the frequency and the intensity of droughts in the Mediterranean region, likely reducing growth and increasing mortality of holm oaks (Quercus ilex), one of the most abundant species of Mediterranean forests. In water-limited systems such as those of the Mediterranean, carbon allocation patterns strongly favour belowground accumulation, especially in large subterranean structures called lignotubers. The resilience of these forests depends largely on the replenishment rate of these carbon reserves after disturbances. An experimental thinning, with two intensities (removal of 40% and 80% of basal area), was performed in 1992 in a holm oak forest at the Prades Experimental Complex of Catchments (NE Spain). In 2002, a second thinning was carried out in subplots within the former experimental 0.5 ha plots. Samples from the lignotubers of holm oak trees were analyzed for starch, and both mobile and immobile chemical components, in order to assess the resilience of holm oaks to repeated disturbances. Our results show that after 10 years, starch stocks in the lignotubers have only recovered to half their former values. Removing 40% of the basal area instead of 80% is suggested to be the better managing option for this kind of forests.     

Genetic differences within natural and planted stands of <Emphasis Type="Italic">Quercus petraea</Emphasis>

Five sessile oak [Quercus petraea (Matt.) Liebl.] stands from the Czech Republic were studied to learn about the impact of different types of forest management regimes on the genetic differences among tree populations and on population structures. One population had not been markedly affected by human activity, two populations represented unplanted stands that were extensively managed for a long period of time using the coppice system, and two populations were planted stands. Approximately 100 trees from each stand were mapped and subsequently genotyped using 10 nuclear microsatellite loci. We determined the spatial genetic structure of each population and the genetic differentiation among the populations. We found that: (i) the populations were genetically differentiated, but the differences between the unplanted and planted stands were not markedly significant; (ii) the genetic differentiation among the populations depended on the geographical distribution of the populations; (iii) within unplanted stands, a strong spatial genetic structure was seen; and (iv) within planted stands, no spatial genetic structure was observed. Our findings implies that the analysis of spatial genetic structure of the sessile oak forest stand can help reveal and determine its origin.     

Reproductive ecology of a cleistogamous annual,Impatiens noli-tangere L., occurring under different environmental conditions

Demographic and reproductive schedules were compared among five populations of a cleistogamous annual,Impatients noli-tangere L., occurring in habitats with contrasting moisture and/or light environments. In all populations, flowering extended for 2–3 months during which light environment and mortality changed. Seasonal patterns of growth and mortality were significantly different among the five populations studied. The beginning and duration of flowering and the ratio of chasmogamous flowers to cleistogamous ones were also significantly different among populations. An experiment was conducted under different light conditions (open and closed) to separate the genetic and environmental components of the variation in reproductive traits observed among populations ofI. noli-tangere. Transplanted plants showed significant among-population variation in flowering time, as is observed in natural populations, suggesting genetic differentiation among populations of this species. On the other hand, the ratio of chasmogamous flowers did not differ among plants transplanted from three populations. Based on these results, the authors suggest that facultative cleistogamy is a conditional strategy under seasonally changing environments.     

Inter-populational variation,but no-annual variation within populations,in terms of reproductive size and genetic structure in a monocarpic perennial herb,Cardiocrinum cordatum var. glehnii

Cardiocrinum cordatum var. glehnii (Liliaceae) is a monocarpic perennial herb living in temperate broad-leaved deciduous forests. In the present study we examined the sizes (basal diameter of the stem and flower numbers) of flowering individuals and genetic diversity using microsatellite loci of 23 populations of C. cordatum var. glehnii in Hokkaido, Japan, over 2 years (2009 and 2010). As a result, we found both the basal stem diameter and the number of flowers varied widely among the populations. However, although the sizes of flowering individuals differed among the populations, these were very stable in each population and in each year. In addition, for genetic diversity, the same trends (i.e. wide variation among the populations but non-annual variation) were detected.     

The Balearic Islands: a reservoir of cpDNA genetic variation for evergreen oaks

Aim To analyse the role of the Balearic Islands as a refuge area for evergreen Quercus (cork oak: Quercus suber L., holm oak: Q. ilex L., kermes oak: Q. coccifera L.), by using molecular, historical and palaeobotanical data. Location The Western Mediterranean Basin (Balearic Islands, eastern Iberia, Provence, Sardinia, Corsica, Sicily, Malta, Italy, Northern Africa). Methods We sampled 108 populations and used the PCR‐RFLP technique with five universal cpDNA primers to define haplotypes in the sampled populations. Diversity, differentiation parameters and spatial analysis of the populations, using a spatial version of amova , were linked to the geological history of the Western Mediterranean Basin in order to explain the present spatial pattern of the evergreen Quercus populations in the Balearics. Results Evergreen Quercus cpDNA shows a complex structure, with remnants of ancient diversity in the Balearics. Balearic populations of holm oak are related to Iberian populations, while for cork and kermes oaks, we found both Tyrrhenian and Iberian haplotypes. Main conclusions The complex spatial patterns of cpDNA in Balearic evergreen Quercus appears explicable in terms of a combination of physical (vicariance and long distance dispersal) and biological (introgressive hybridization) factors. The Balearics constitute a glacial refuge area and a reservoir of genetic variation with traces of ancient diversity from Messinian–Pliocene stages.     

Efficient mitigation of founder effects during the establishment of a leading-edge oak population

Numerous plant species are shifting their range polewards in response to ongoing climate change. Range shifts typically involve the repeated establishment and growth of leading-edge populations well ahead of the main species range. How these populations recover from founder events and associated diversity loss remains poorly understood. To help fill this gap, we exhaustively investigated a newly established population of holm oak (Quercus ilex) growing more than 30 km ahead of the nearest larger stands. Pedigree reconstructions showed that plants belong to two non-overlapping generations and that the whole population originates from only two founder trees. The four first-generation trees that have reached maturity showed disparate mating patterns despite being full-sibs. Long-distance pollen immigration was notable despite the strong isolation of the stand: 6 per cent gene flow events in acorns collected on the trees (n = 255), and as much as 27 per cent among their established offspring (n = 33). Our results show that isolated leading-edge populations of wind-pollinated forest trees can rapidly restore their genetic diversity through the interacting effects of efficient long-distance pollen flow and purging of inbred individuals during recruitment. They imply that range expansions of these species are primarily constrained by initial propagule arrival rather than by subsequent gene flow.     

Growth and yield of young Quercus ilex coppice stands in the Tafferte forest (Morocco)

The coppice method is the silvicultural system commonly applied to Quercus ilex (holm oak) forests for the production of firewood and charcoal. Although the holm oak is the species that covers the largest area in Morocco (1400 000 hectares), little is known about its growth and yield. The objective of this study is to quantify the stand growth and yield of young coppice holm oak stands in the Tafferete forest in the eastern part of the Middle Atlas (Morocco).Equations relating dominant height to age, and basal area, volume, and number of trees to dominant height were developed and used to derive young holm oak coppice yield estimates.The growth and yield information developed in this study may have a greater utility in defining cultural measures that should be applied in the early stage of holm-oak coppice stands.     

Natural hybridisation between kermes (Quercus coccifera L.) and holm oaks (Q. ilex L.) revealed by microsatellite markers

Hybridisation between species of the genus Quercus is a common phenomenon as a result of weak reproductive isolation mechanisms between phylogenetically close species that frequently co-occur in mixed stands. In this study, we use microsatellite markers to analyse introgression between kermes ( Quercus coccifera L.) and holm ( Q. ilex L.) oak, two closely related taxa that frequently dominate the landscape in extensive areas in the Mediterranean region. All tested microsatellites amplified and were polymorphic in both kermes and holm oaks. Bayesian admixture analyses showed a good correspondence between each species and one of the two inferred genetic clusters. Five sampled individuals were a priori tentatively identified as hybrids on the basis of intermediate morphological characteristics, and it was confirmed that they also presented mixed genotypes. However, we also detected different levels of genetic introgression among morphologically pure individuals, suggesting that successful backcrossing and/or reduced phenotypic expression of genetic variance in certain individuals may have resulted in strong convergence towards a single species phenotype.     

Substantial genetic diversity in cultivated Moroccan olive despite a single major cultivar: a paradoxical situation evidenced by the use of SSR loci

To assess the genetic diversity in Moroccan cultivated olive, Olea europaea L. subsp. europaea, we performed molecular analysis of olive trees sampled in four geographic zones representing all areas of traditional olive culture. The analysis of 215 trees using 15 simple sequence repeat (SSR) loci revealed 105 alleles distributed among 60 SSR profiles. The analysis of chloroplast deoxyribonucleic acid polymorphism for these 60 olive genotypes allowed to identify four chlorotypes: 42 CE1, one CE2, nine COM1 and eight CCK. Among the 60 SSR profiles, 52 corresponded to cultivated olive trees for which neither denomination nor characterisation is available. These local olive genotypes displayed a spatial genetic structuring over the four Moroccan geographic zones (northwest, north centre, Atlas and southwest), as pairwise Fst values ranged from 0.0394 to 0.1383 and varied according to geographic distance. As single alleles detected in local olive were also observed in Moroccan oleaster populations, results suggest that plant material was mainly selected from indigenous populations. The assumption that Picholine marocaine cultivar is a multi-clonal cultivar was not supported by our data because we found a single genotype for 112 olive trees representing 31 to 93% of the olives sampled locally in the 14 different areas. Picholine marocaine and the few other named cultivars do not seem to belong to the same gene pools as the unnamed genotypes cultivated only locally. The situation is paradoxical: a substantial genetic diversity in Moroccan olive germplasm, probably resulting from much local domestication, but a single cultivar is predominant.     

Studying genetics of adaptive variation in model organisms: flowering time variation in Arabidopsis lyrata

Arabidopsis thaliana has emerged as a model organism for plant developmental genetics, but it is also now being widely used for population genetic studies. Outcrossing relatives of A. thaliana are likely to provide suitable additional or alternative species for studies of evolutionary and population genetics. We have examined patterns of adaptive flowering time variation in the outcrossing, perennial A. lyrata. In addition, we examine the distribution of variation at marker genes in populations form North America and Europe. The probability of flowering in this species differs between southern and northern populations. Northern populations are much less likely to flower in short than in long days. A significant daylength by region interaction shows that the northern and southern populations respond differently to the daylength. The timing of flowering also differs between populations, and is made shorter by long days, and in some populations, by vernalization. North American and European populations show consistent genetic differentiation over microsatellite and isozyme loci and alcohol dehydrogenase sequences. Thus, the patterns of variation are quite different from those in A. thaliana, where flowering time differences show little relationship to latitude of origin and the genealogical trees of accessions vary depending on the genomic region studied. The genetic architecture of adaptation can be compared in these species with different life histories.     

Development and feeding efficiency of Malacosoma neustrium larvae reared with Quercus spp. leaves

Observations on larval development of Malacosoma neustrium were conducted both in a cork oak stand and in the laboratory by using leaves of different host trees (cork oak, holm oak and downy oak) as food source. Instars were determined using head capsule and frass measurements. In the field the larvae progressed up to the fifth instar before pupating, and the increase in head capsule width followed Dyars Rule with a rate of increase (R.I.) value of 1.74. The same number of instars was determined for the larvae reared with cork oak (R.I. = 1.73) and holm oak (R.I. = 1.70) leaves. The caterpillars reared with downy oak foliage completed larval development in five, six and seven instars and the R.I. values obtained were 1.60, 1.52 and 1.44 respectively. A lower mortality was recorded for the larvae reared on holm oak. Growth and feeding indices were determined for the larvae from the third up to the last instar. The highest leaf consumption was detected for the fifth instar larvae reared on holm oak. For the caterpillars which completed five instars before pupating, the relative consumption rate (RCR) decreased from the third up to the fifth instar: from 4.8 to 1.7 (cork oak), from 7.4 to 3.3 (holm oak) and from 14.3 to 2.1 (downy oak). The relative growth rate (RGR) was highest during the fourth stadium (0.24, 0.27 and 0.33 for larvae reared with cork oak, holm oak and downy oak leaves respectively) and decreased in the fifth instar (0.09, 0.14 and 0.14 for larvae reared with cork oak, holm oak and downy oak leaves respectively), probably because of greater expense of energy due to the approach of maturity. Feeding and growth indices could be useful to define a defoliation prediction model.     

Estimation of Mediterranean forest transpiration and photosynthesis through the use of an ecosystem simulation model driven by remotely sensed data

Aim This paper investigates the use of an ecosystem simulation model, FOREST‐BGC, to estimate the main ecophysiological processes (transpiration and photosynthesis) of Mediterranean coastal forest areas using remotely sensed data. Location Model testing was carried out at two protected forest sites in central Italy, one of which was covered by Turkey oak (Circeo National Park) and the other by holm‐oak (Castelporziano Estate). Methods At both sites, transpiration and photosynthesis measurements were collected in the field during the growing seasons over a four‐year period (1999 and 2001 for the Turkey oak; 1997, 1999 and 2000 for the holm‐oak). Calibration of the model was obtained through combining information derived from ground measurements and remotely sensed data. In particular, remote sensing estimates of the Leaf Area Index derived from 1 × 1‐km NOAA AVHRR Normalized Difference Vegetation Index data were used to improve the adaptation of the model to local forest conditions. Results The results indicated different strategies regarding water use efficiency, ‘water spending’ for Turkey oak and ‘water saving’ for holm‐oak. The water use efficiency for the holm‐oak was consistently higher than that for the Turkey oak and the relationship between VPD and WUE for the holm‐oak showed a higher coefficient of determination (R² = 0.9238). Comparisons made between the field measurements of transpiration and photosynthesis and the model estimates showed that the integration procedure used for the deciduous oak forest was effective, but that there is a need for further studies regarding the sclerophyllous evergreen forest. In particular, for Turkey oak the simulations of transpiration yielded very good results, with errors lower than 0.3 mm H₂O/day, while the simulation accuracy for photosynthesis was lower. In the case of holm‐oak, transpiration was markedly overestimated for all days considered, while the simulations of photosynthesis were very accurate. Main conclusions Overall, the approach offers interesting operational possibilities for the monitoring of Mediterranean forest ecosystems, particularly in view of the availability of new satellite sensors with a higher spatial and temporal resolution, which have been launched in recent years.     

Compensatory growth in shoot populations of young white pine trees

Summary White pine (Pinus strobus L.) trees have shoot populations composed of subpopulations of terminal and lateral shoots. I tested whether the subpopulations would show compensatory (increased) growth when separated from each other. Ten-year-old white pine (Pinus strobus L.) trees growing under an oak (Quercus) overstory were untreated or treated in winter by removing either all terminal, or all lateral buds (10 trees per treatment). Growth was compared between control and treated shoot subpopulations. In the 1st year, shoot-length frequency distributions were similar between control and treated subpopulations. There was significant compensatory shoot elongation (mean of 1.5 cm per shoot) in both treated subpopulations. In the 2nd year each subpopulation produced both terminal and lateral shoots. Shoot-length frequency distributions were similar, but shifted toward longer shoots in treated populations. Shoot number, mean length and total shoot length were greater in treated populations. The increased growth in treated subpopulations was due both to differences in parent shoot length and to compensatory shoot production and elongation.     

Escape from floral herbivory by early flowering in Arabidopsis halleri subsp. gemmifera

Natural selection on flowering phenology has been studied primarily in terms of plant–pollinator interactions and effects of abiotic conditions. Little is known, however, about geographic variation in other biotic factors such as herbivores and its consequence for differential selection on flowering phenology among populations. Here, we examine selection by floral herbivores on the flowering phenology of Arabidopsis halleri subsp. gemmifera using two adjacent populations with contrasting herbivory regimes. Intensive floral herbivory by the leaf beetle Phaedon brassicae occurs in one population, while the beetle is absent in another population. We tested the hypothesis that the two populations experience differential selection on flowering time that is attributable to the presence or absence of floral herbivory. A two-year field study showed that early flowering was favoured in the population under intensive floral herbivory, whereas selection for early flowering was not found in one year in the population where floral herbivory was absent. Selection for early flowering disappeared when the abundance of floral herbivores was artificially decreased in a field experiment. Thus, the heterogeneous distribution of P. brassicae was a major agent for differential selection on flowering time. However, flowering time did not differ between the two populations when plants were grown in the laboratory. The lack of genetic differentiation in flowering time may be explained by ongoing gene flow or recent invasion of P. brassicae. This study illustrates that the role of floral herbivory in shaping geographic variation in selection on flowering phenology may be more important than previously thought.     

Environmentally-based maternal effects: a hidden force in insect population dynamics?

Summary The nutritional environment of the parental generation of the polyphagous gypsy moth, Lymantria dispar, can significantly influence the growth and reproductive potential of the next generation through environmentally-based maternal effects. In the first experiment, members of the parental generation were reared on red oak trees (Quercus rubra L.) with known defoliation and phenolic levels. Diet in the offspring generation was homogeneous (synthetic diet). With genetic effects accounted for 1) offspring attained greater pupal weights when their mothers fed on trees with higher leaf damage levels, 2) daughters had a shorter prefeeding stage, a trait associated with dispersal tendency, when their mothers experienced higher condensed tannin levels, and 3) sons had lower pupal weights when their mothers experienced greater condensed tannin levels. In the second experiment, members of the parental generation were reared on either red or black oak (Q. velutina) trees. Offspring of each mother were divided among four diets: red oak, chestnut oak (Q. prinus L.), a standard synthetic diet, and a low-protein synthetic diet. The parental host species accounted for 24% of the variation in daughters' development time; offspring diet accounted for 52%. When mothers were reared on black oak rather than red oak, their offspring developed significantly faster when the F1 diet was chestnut oak. Environmentally-based maternal effects can significantly influence the expression of offspring dispersal potential, growth rate, and offspring fecundity. These traits contribute to natality and survival in natural populations and, hence, to population growth potential. Theoretical models of insect population dynamics demonstrate that the presence of a time delay in a density dependent response can induce destabilization. Maternal effects act on a time delay and may participate in the destabilization characteristic of outbreak species.