The Difficulty of Predicting Evolutionary Change in Response to Novel Ecological Interactions: A Field Experiment with Anolis Lizards.

AbstractDetermining whether and how evolution is predictable is an important goal, particularly as anthropogenic disturbances lead to novel species interactions that could modify selective pressures. Here, we use a multigeneration field experiment with brown anole lizards (Anolis sagrei) to test hypotheses about the predictability of evolution. We manipulated the presence/absence of predators and competitors of A. sagrei across 16 islands in the Bahamas that had preexisting brown anole populations. Before the experiment and again after roughly five generations, we measured traits related to locomotor performance and habitat use by brown anoles and used double-digest restriction enzyme-associated DNA sequencing to estimate genome-wide changes in allele frequencies. Although previous work showed that predators and competitors had characteristic effects on brown anole behavior, diet, and population sizes, we found that evolutionary change at both phenotypic and genomic levels was difficult to forecast. Phenotypic changes were contingent on sex and habitat use, whereas genetic change was unpredictable and not measurably correlated with phenotypic changes, experimental treatments, or other environmental factors. Our work shows how differences in ecological context can alter evolutionary outcomes over short timescales and underscores the difficulty of forecasting evolutionary responses to multispecies interactions in natural conditions, even in a well-studied system with ample supporting ecological information.

Predator-induced collapse of niche structure and species coexistence.

Biological invasions are both a pressing environmental challenge and an opportunity to investigate fundamental ecological processes, such as the role of top predators in regulating biodiversity and food-web structure. In whole-ecosystem manipulations of small Caribbean islands on which brown anole lizards (Anolis sagrei) were the native top predator, we experimentally staged invasions by competitors (green anoles, Anolis smaragdinus) and/or new top predators (curly-tailed lizards, Leiocephalus carinatus). We show that curly-tailed lizards destabilized the coexistence of competing prey species, contrary to the classic idea of keystone predation. Fear-driven avoidance of predators collapsed the spatial and dietary niche structure that otherwise stabilized coexistence, which intensified interspecific competition within predator-free refuges and contributed to the extinction of green-anole populations on two islands. Moreover, whereas adding either green anoles or curly-tailed lizards lengthened food chains on the islands, adding both species reversed this effect-in part because the apex predators were trophic omnivores. Our results underscore the importance of top-down control in ecological communities, but show that its outcomes depend on prey behaviour, spatial structure, and omnivory. Diversity-enhancing effects of top predators cannot be assumed, and non-consumptive effects of predation risk may be a widespread constraint on species coexistence.

Phylogeography of the Neotropical epiphytic orchid, Brassavola nodosa: evidence for a secondary contact zone in northwestern Costa Rica.

Spatial patterns of genetic variation can reveal otherwise cryptic evolutionary and landscape processes. In northwestern Costa Rica, an approximately concordant genetic discontinuity occurs among populations of several plant species. We conducted phylogeographic analyses of an epiphytic orchid, Brassavola nodosa, to test for genetic discontinuity and to explore its underlying causes. We genotyped 18 populations with 19 nuclear loci and two non-coding chloroplast sequence regions. We estimated genetic diversity and structure, relative importance of pollen and seed dispersal, and divergence time to understand how genetic diversity was spatially partitioned. Nuclear genetic diversity was high with little differentiation among populations (GSTn = 0.065). In contrast, chloroplast haplotypes were highly structured (GSTc = 0.570) and reveal a discontinuity between northwestern and southeastern populations within Costa Rica. Haplotype differences suggest two formerly isolated lineages that diverged ~10,000-100,000 YBP. Haplotype mixing and greater genetic diversity occur in an intermediate transition zone. Patterns of nuclear and chloroplast data were consistent. Different levels of genetic differentiation for the two genomes reflect the relative effectiveness of biotic versus abiotic dispersers of pollen and seeds, respectively. Isolation of the two lineages likely resulted from the complex environmental and geophysical history of the region. Our results suggest a recent cryptic seed dispersal barrier and/or zone of secondary contact. We hypothesize that powerful northeasterly trade winds hinder movement of wind-borne seeds between the two regions, while the multi-directional dispersal of pollen by strong-flying sphinx moths resulted in lower differentiation of nuclear loci.

Molecular detection of invertebrate prey in vertebrate diets: trophic ecology of Caribbean island lizards.

Understanding community assembly and population dynamics frequently requires detailed knowledge of food web structure. For many consumers, obtaining precise information about diet composition has traditionally required sacrificing animals or other highly invasive procedures, generating tension between maintaining intact study populations and knowing what they eat. We developed 16S mitochondrial DNA sequencing methods to identify arthropods in the diets of generalist vertebrate predators without requiring a blocking primer. We demonstrate the utility of these methods for a common Caribbean lizard that has been intensively studied in the context of small island food webs: Anolis sagrei (a semi-arboreal 'trunk-ground' anole ecomorph). Novel PCR primers were identified in silico and tested in vitro. Illumina sequencing successfully characterized the arthropod component of 168 faecal DNA samples collected during three field trips spanning 12 months, revealing 217 molecular operational taxonomic units (mOTUs) from at least nine arthropod orders (including Araneae, Blattodea, Coleoptera, Hemiptera, Hymenoptera, Isoptera, Lepidoptera and Orthoptera). Three mOTUs (one beetle, one cockroach and one ant) were particularly frequent, occurring in ≥50% of samples, but the majority of mOTUs were infrequent (180, or 83%, occurred in ≤5% of samples). Species accumulation curves showed that dietary richness and composition were similar between size-dimorphic sexes; however, female lizards had greater per-sample dietary richness than males. Overall diet composition (but not richness) was significantly different across seasons, and we found more pronounced interindividual variation in December than in May. These methods will be generally useful in characterizing the diets of diverse insectivorous vertebrates.

Relative importance of pollen and seed dispersal across a Neotropical mountain landscape for an epiphytic orchid.

Populations of many species are isolated within narrow elevation bands of Neotropical mountain habitat, and how well dispersal maintains genetic connectivity is unknown. We asked whether genetic structure of an epiphytic orchid, Epidendrum firmum, corresponds to gaps between Costa Rican mountain ranges, and how these gaps influence pollen and seed flow. We predicted that significant genetic structure exists among mountain ranges due to different colonization histories and limited gene flow. Furthermore, we predicted that pollen movement contributes more to gene flow than seeds because seeds are released into strong winds perpendicular to the narrow northwest-southeast species distribution, while the likely pollinators are strong fliers. Individuals from 12 populations and three mountain ranges were genotyped with nuclear microsatellites (nDNA) and chloroplast sequences (cpDNA). Genetic diversity was high for both markers, while nDNA genetic structure was low (FSTn  = 0.020) and cpDNA structure was moderate (FSTc  = 0.443). Significant cpDNA barriers occurred within and among mountain ranges, but nDNA barriers were not significant after accounting for geographic distance. Consistent with these contrasting patterns of genetic structure, pollen contributes substantially more to gene flow among populations than seed (mp /ms  = 46). Pollinators mediated extensive gene flow, eroding nDNA colonization footprints, while seed flow was comparatively limited, possibly due to directional prevailing winds across linearly distributed populations. Dispersal traits alone may not accurately inform predictions about gene flow or genetic structure, supporting the need for research into the potentially crucial role of pollinators and landscape context in gene flow among isolated populations.

Highly diverse and spatially heterogeneous mycorrhizal symbiosis in a rare epiphyte is unrelated to broad biogeographic or environmental features.

Symbiotic interactions are common in nature. In dynamic or degraded environments, the ability to associate with multiple partners (i.e. broad specificity) may enable species to persist through fluctuations in the availability of any particular partner. Understanding how species interactions vary across landscapes is necessary to anticipate direct and indirect consequences of environmental degradation on species conservation. We asked whether mycorrhizal symbiosis by populations of a rare epiphytic orchid (Epidendrum firmum) is related to geographic or environmental heterogeneity. The latter would suggest that interactions are governed by environmental conditions rather than historic isolation of populations and/or mycorrhizal fungi. We used DNA-based methods to identify mycorrhizal fungi from eleven E. firmum populations in Costa Rica. We used molecular and phylogenetic analyses to compare associations. Epidendrum firmum exhibited broad specificity, associating with diverse mycorrhizal fungi, including six Tulasnellaceae molecular operational taxonomic units (MOTUs), five Sebacinales MOTUs and others. Notably, diverse mycorrhizal symbioses formed in disturbed pasture and roadside habitats. Mycorrhizal fungi exhibited significant similarity within populations (spatial and phylogenetic autocorrelation) and significant differences among populations (phylogenetic community dissimilarity). However, mycorrhizal symbioses were not significantly associated with biogeographic or environmental features. Such unexpected heterogeneity among populations may result from complex combinations of fine-scale environmental factors and macro-evolutionary patterns of change in mycorrhizal specificity. Thus, E. firmum exhibits broad specificity and the potential for opportunistic associations with diverse fungi. We suggest that these characteristics could confer symbiotic assurance when mycorrhizal fungi are stochastically available, which may be crucial in dynamic or disturbed habitats such as tropical forest canopies.

Genetic inference of epiphytic orchid colonization; it may only take one.

Colonization of vacant habitat is a fundamental ecological process that affects the ability of species to persist and undergo range modifications in continually shifting landscapes. Thus, understanding factors that affect and limit colonization has important ecological and conservation implications. Epiphytic orchids are increasingly threatened by various factors, including anthropogenic habitat disturbance. As cleared areas (e.g. pastures) are recolonized by suitable host trees, the establishment and genetic composition of epiphytic orchid populations are likely a function of their colonization patterns. We used genetic analyses to infer the prevailing colonization pattern of the epiphytic orchid, Brassavola nodosa. Samples from three populations (i.e. individuals within a tree) from each of five pastures in the dry forest of Costa Rica were genotyped with neutral nuclear and chloroplast markers. Spatial autocorrelation and hierarchical genetic structure analyses were used to assess the relatedness of individuals within populations, among populations within pastures and among populations in different pastures. The results showed significant relatedness within populations (mean r=0.166) and significant but lower relatedness among populations within a pasture (mean r=0.058). Our data suggest that colonization of available habitats is by few individuals with subsequent population expansion resulting from in situ reproduction, and that individuals within a tree are not a random sample of the regional seed pool. Furthermore, populations within a pasture were likely colonized by seeds produced by founders of a neighbouring population within that pasture. These results have important ramifications for understanding conservation measures needed for this species and other epiphytic orchids.

Microsatellite primers for the neotropical epiphyte Epidendrum firmum (Orchidaceae).

PREMISE OF THE STUDY: Ten microsatellite loci were isolated and characterized for the neotropical epiphytic orchid Epidendrum firmum to examine levels of genetic diversity and genetic structure at multiple spatial scales. METHODS AND RESULTS: We screened loci in 12-25 individuals from each of two populations in Costa Rica and identified 10 polymorphic loci. The number of alleles per locus ranged from one to 15 while observed heterozygosity for polymorphic loci ranged from 0.360 to 0.960. CONCLUSIONS: Primers for these informative genetic markers will be useful for quantifying genetic diversity, spatial genetic structure, and gene flow in E. firmum.