Genetic and phenotypic variation in a colourful treefrog across five geographic barriers

Aim Understanding the patterns and processes underlying phenotype in a polytypic species provides key insights into microevolutionary mechanisms of diversification. The red‐eyed treefrog, Agalychnis callidryas, exhibits strong regional differentiation in colour pattern, corresponding to five admixed mitochondrial DNA clades. We evaluated spatial diversity patterns across multiple, putative barriers to examine the fine‐scale processes that mediate phenotypic divergence between some regions while maintaining homogeneity between others. Location We examined patterns of phenotypic diversification among 17 sites that span five putative biogeographic barriers in lower Central America (Costa Rica and Panama). Methods We tested the extent to which genetic distance (F_ST) derived from six multilocus nuclear genotypes covaried with measures of phenotypic distance (leg coloration) within and between biogeographic regions. We used linear regression analyses to determine the role of geographic and genetic factors in structuring spatial patterns of phenotypic diversity. Results The factors that best explained patterns of phenotypic diversity varied among biogeographic regions. We identified one geographic barrier that impeded gene exchange and resulted in concordant phenotypic divergence across the Continental Divide, isolating Caribbean and Pacific populations. Across Caribbean Costa Rican populations, one barrier structured phenotypic but not genetic diversity patterns, indicating a role for selection. In other regions, the putative barriers had no determining effect on either genetic or leg colour structure. Main conclusions The processes mediating the distribution and diversification of colour pattern in this polytypic, wide‐ranging treefrog varied among biogeographic regions. Spatially varying selection combined with the isolating effects of geographic factors probably resulted in the patchy distribution of colour diversity across Costa Rican and Panamanian populations.     

To hatch and hatch not: similar selective trade-offs but different responses to egg predators in two closely related, syntopic treefrogs

Risk-sensitive hatching is adaptive for species facing a trade-off between egg-stage and post-hatching risks, and environmental variation in one or both stages. Such plasticity has been found in amphibians, fishes, reptiles and spiders, with red-eyed treefrogs (Agalychnis callidryas) being the best-studied case. We assessed hatching plasticity and egg- and larval-stage risks in a closely related, syntopic species, the gliding leaf-frog (Agalychnis spurrelli). We found a lower hatching response to egg-eating snakes in A. spurrelli (9–28% of embryos escaped) than in A. callidryas (59–80% escaped). Levels of snake predation were similarly high for clutches of both species monitored at a pond in Costa Rica, and in fish predation experiments early-hatched A. spurrelli tadpoles were more vulnerable than later hatchlings, as has been shown for A. callidryas. A. spurrelli thus face a risk trade-off similar to A. callidryas, and likely would benefit from predator-induced hatching; their lower responsiveness to snakes appears nonadaptive. A. spurrelli embryos showed a stronger hatching response (57% hatched in 1 h) to submergence underwater than to snake attacks even though submergence is a less frequent risk. This suggests they have a greater capacity for early hatching than is expressed in the context of snake attacks, but have much lower sensitivity to snake cues than to flooding cues. Development in A. spurrelli is accelerated compared to syntopic A. callidryas, and spontaneous hatching is earlier and more synchronous. This is congruent with predictions based on selection by egg predators in the absence of a strong escape hatching response.     

Lateralized Turning Biases in Two Neotropical Tadpoles

Lateralized turning behavior in startle responses and upon descent after surfacing for a breath of air has been documented for tadpoles in several anuran species. A left‐handed preference is most common and was previously thought to be linked to the asymmetry in spiracle location. Here, we investigate the presence of behavioral asymmetries in tadpoles of Agalychnis callidryas and Leptodactylus melanonotus in explosive turns after air‐breathing. Data were also collected on lateralized posture in the orientation of tails of embryonic A. callidryas within the egg case, as well as the startle response of free‐swimming tadpoles exposed to mechanical stimulation. A left‐curled tail bias was found among several clutches from Costa Rica, but this was not the case among clutches examined from Panama. Free‐living tadpoles of L. melanonotus displayed a distinct right‐handed preference during explosive turns. While some tadpoles of A. callidryas were at stages too early to detect any explosive turns when breathing, later‐staged individuals did display a left‐handed bias in startle response to mechanical stimulation. Additionally, it appears that the expression of behavioral lateralization of larvae (but not embryos) may predict whether or not the adults exhibit similar lateralization. Findings herein provide insight into the ontogeny and evolutionary origins of lateralized behavior in anurans.     

Patterns of reproduction and habitat use in an assemblage of Neotropical hylid frogs

We censused pond-breeding hylid frogs in northeastern Costa Rica weekly for 15 months to deseribe patterns of reproduction, habitat use, and to establish baseline data on relative abundance for members of the assemblage. Reproduction in the Costa Rican assemblage was seasonal and occurred only during wet months. Some species called, but none reproduced, during the dry season. Three species (Agalychnis callidryas, Hyla ebraccata, and Scinax elaeochroa) accounted for more than 75% of the observations made during the study. The species overlapped broadly in time and space, but differed in substrate use and phenology. Two species of leaf-breeding frogs (A. callidryas and A. saltator) used perches that were significantly higher than those used by the other species. Some phenological differences were associated with different mating strategies. Explosive breeders (Scinax elaeochroa and Smilisca baudinii) were most common early in the wet seasons. Prolonged breeders (A. callidryas and H. ebraccata) were the most persistent members of this assemblage. Predation affects early and late life history stages of these hylids. Predation on arboreal egg masses by two snake species was observed. Ctenid spiders preyed on recently metamorphosed frogs and small adults. Our weekly samples were pooled into 21-day periods so that we could compare our results with those obtained for two communities of breeding anurans from South America. The patterns observed in the Costa Rican assemblage differed from those reported for South American pond-breeding frogs, but in all three assemblages reproduction was associated with wet periods.     

Melodic males and flashy females: Geographic variation in male and female reproductive behavior in red‐eyed treefrogs (Agalychnis callidryas)

Geographic variation in courtship behavior can affect reproductive success of divergent phenotypes via mate choice. Over time, this can lead to reproductive isolation and ultimately to speciation. The Neotropical red‐eyed treefrog (Agalychnis callidryas) exhibits high levels of phenotypic variation among populations in Costa Rica and Panama, including differences in color pattern, body size, and skin peptides. To test the extent of behavioral premating isolation among differentiated populations, we quantified male advertisement calls from six sites and female responses to male stimuli (acoustic and visual signals) from four sites. Our results show that both male advertisement calls and female behavior vary among populations: Discriminant function analyses can predict the population of origin for 99.3% ± 0.7 of males based on male call (dominant frequency and bandwidth) and 76.1% ± 6.6 of females based on female response behavior (frequency and duration of visual displays). Further, female mate choice trials (n = 69) showed that population divergence in male signals is coupled with female preference for local male stimuli. Combined, these results suggest that evolved differences among populations in male call properties and female response signals could have consequences for reproductive isolation. Finally, population variation in male and female behavior was not well explained by geographic or genetic distance, indicating a role for localized selection and/or drift. The interplay between male courtship and female responses may facilitate the evolution of local variants in courtship style, thus accelerating premating isolation via assortative mating.     

Sexual signal evolution and patterns of assortative mating across an intraspecific contact zone

Contact zones provide important insights into the evolutionary processes that underlie lineage divergence and speciation. Here, we use a contact zone to ascertain speciation potential in the red-eyed treefrog (Agalychnis callidryas), a brightly coloured and polymorphic frog that exhibits unusually high levels of intraspecific variation. Populations of A. callidryas differ in a number of traits, several of which are known sexual signals that mediate premating reproductive isolation in allopatric populations. Along the Caribbean coast of Costa Rica, a ~100 km contact zone, situated between two phenotypically and genetically divergent parent populations, contains multiple colour pattern phenotypes and late-generation hybrids. This contact zone provides the opportunity to examine processes that are important in the earliest stages of lineage divergence. We performed analyses of colour pattern variation in five contact zone sites and six parental sites and found complex, continuous colour variation along the contact zone. We found discordance between the geographic distribution of colour pattern and previously described genomic population structure. We then used a parental site and contact zone site to measure assortative mating and directional selection from naturally-occurring amplectant mating pairs. We found assortative mating in a parental population, but no assortative mating in the contact zone. Furthermore, we uncovered evidence of directional preference towards the adjacent parental phenotype in the contact zone population, but no directional preference in the parent population. Combined, these data provide insights into potential dynamics at the contact zone borders and indicate that incipient speciation between parent populations will be slowed.     

Distribution and diversity of begomoviruses in tomato and sweet pepper plants in Costa Rica

Begomoviruses (genus Begomovirus, family Geminiviridae) have emerged as important plant pathogens in tropical and subtropical regions worldwide. Although these viruses were reported during the 1970s in Costa Rica, they are still poorly known. Therefore, the objective of this study was to analyse the diversity and distribution of begomoviruses in commercial tomato and sweet pepper fields from different agricultural production systems of the major growing regions of Costa Rica. A total of 651 plants were randomly sampled from greenhouses and open field crops during 2011 and 2012 in three different geographical locations. The bipartite begomoviruses Tomato yellow mottle virus, Tomato leaf curl Sinaloa virus and Pepper golden mosaic virus, and the monopartite begomovirus Tomato yellow leaf curl virus were detected in the collected samples. The complete genome of isolates from each species was cloned and sequenced. The frequency of detection of these four begomoviruses in the analysed samples ranged from 0 to 9%, the presence, and the prevalent virus varied largely according to the geographical location, the host (tomato and pepper), and the production system (greenhouses or open fields). An association between geographical region and begomovirus species was observed suggesting that in Costa Rica the heterogeneity on climate, topography and agricultural system might influence the distribution of begomovirus species in the country. A broader survey needs to be conducted to confirm it, although these preliminary results may contribute to the management of begomoviruses in Costa Rica.     

Biogeographical analysis of two Polypodium species complexes (Polypodiaceae) in Mexico and Central America

We analyzed the geographical and elevational distributions of two Polypodium complexes from Mexico and Central America. Distribution data of nine species of the Polypodium colpodes complex and the Polypodium plesiosorum complex were obtained from almost 1500 herbarium specimens, field collections in Mexico and Costa Rica, and literature studies. The presence of each species was recorded for each Mesoamerican country, in 1° × 1° grid‐cells and biogeographical provinces. The rarity of species was also evaluated. Although the two complexes show extensive overlap, the P. colpodes complex is distributed mainly along the Pacific versant of Mexico and Central America, whereas the P. plesiosorum complex occurs mainly along the Atlantic versant. Those biogeographical provinces with maximum species diversity are Chiapas (seven species), Sierra Madre del Sur (six species), and the Trans‐Mexican Volcanic belt (six species). Grid‐cells with more species are located mainly in the mountains of central‐southern Mexico and northern Central America. Richness does not decrease or increase with latitude. Elevation distributions showed that most Polypodium species are concentrated in the montane interval and three species groups were recognized based on elevational preferences. Polypodium colpodes and P. plesiosorum are the most widely distributed species, whereas Polypodium castaneum and Polypodium flagellare are the only two species that possess the three attributes of rarity (narrow geographical distribution, high habitat specificity, and scarce local populations). Polypodium species of both complexes are present mainly in the montane regions of the study area and show some degree of geographical sympatry, especially in southern Mexico and northern Central America. This overlapping is explained by the elevation tolerance within montane systems and because most species inhabit three or more vegetation types. The distributional patterns of these complexes coincided with the three regional highlands of Mesoamerica, which are separated from each other by the Isthmus of Tehuantepec and by the lowlands of Nicaragua. © 2012 The Linnean Society of London, Biological Journal of the Linnean Society, 2012, ?? , ??–??.     

Interactions Between Competition and Predation Shape Early Growth and Survival of Two Neotropical Hylid Tadpoles

Experimental studies in temperate regions have revealed that competition and predation interact to shape aquatic communities. Predators typically reduce the effect of competition on growth and competitors provide alternative prey subjects, which may also alter predation. Here, we examine the independent and combined effects of competition and predation on the survival and growth of hatchling tadpoles of two widespread co‐occurring Neotropical hylid frogs (Agalychnis callidryas and Dendropsophus ebraccatus). Using 400 L mesocosms, we used a 2 × 3 factorial substitutive design, which crossed tadpole species composition with the presence or absence of a free‐roaming predator (Anax amazili dragonfly larva). Dragonflies were effective predators of both species, but had larger effects on A. callidryas survival. Both species had similar growth rates when alone, whereas A. callidryas grew 30 percent faster than D. ebraccatus when they co‐occurred, suggesting interspecific rather than intraspecific competition had relatively stronger effects on D. ebraccatus growth, while the opposite was true for A. callidryas. Predator presence dramatically reduced growth rates of both species and erased this asymmetry. Results suggest that the effects of predator induction (i.e., nonconsumptive effects) on growth were larger than both consumptive and competitive effects. Our study demonstrates that predators have strong effects on both survival and growth of prey, highlighting the potential importance of predators in shaping prey populations and tropical aquatic food web interactions. Abstract in Spanish is available at http://www.blackwell‐synergy.com/loi/btp .     

Differential regional response of the bushy‐tailed woodrat (Neotoma cinerea) to late Quaternary climate change

Aim To reconstruct the regional biogeographical history of the bushy‐tailed woodrat, Neotoma cinerea (Rodentia: Cricetidae), across its distribution using multiple sources of information, including genetic data, ecological niche models and the palaeorecord. Location Western North America. Methods We analysed complete cytochrome b gene (1143 bp) sequences from 182 specimens of N. cinerea using Bayesian and coalescent methods to infer phylogenetic relationships, time of major divergences, and recent demographic trends. For comparison, we developed clade‐specific ecological niche models for groups of interest and analysed spatial trends of N. cinerea in the palaeorecord as well as temporal frequency trends across strata of individual palaeomiddens. Results We found two largely allopatric clades within N. cinerea, with several regionally distinct subclades showing contrasting recent population dynamics. Niche models showed consistent habitat at the Last Glacial Maximum (LGM) and modern times in the Rocky Mountains and northern United States, while the Great Basin may have been markedly less suitable at the LGM than today. The palaeorecord showed great spatio‐temporal variability in the presence of N. cinerea, but documents broad‐scale patterns of occupancy and regional population trends. Main conclusions The Quaternary dynamics and evolutionary history of N. cinerea appear to have been shaped by both vicariant events associated with geographical barriers and the availability of suitable habitat through time. Divergence of the two major clades dates to the Pliocene–Pleistocene transition, with clades separated by the Green and Colorado rivers and northern Rocky Mountain Pleistocene glaciations. We found largely concordant genetic, niche model and palaeorecord patterns suggesting long‐term population stability in the Rocky Mountains, while extant clades in the Great Basin and far north appear to have expanded or re‐expanded into these regions relatively recently. Furthermore, disjunct haplotype distributions, regional demographic history, and historical distribution of suitable habitat suggest that the Great Basin has been a particularly dynamic region.     

Biogeography at the limits of life: Do extremophilic microbial communities show biogeographical regionalization?

Aim

Biogeographical regions are the fundamental geographical units for grouping Earth's biodiversity. Biogeographical regionalization has been demonstrated for many higher taxa, such as terrestrial plants and vertebrates, but not in microbial communities. Therefore, we sought to test empirically whether microbial communities, or taxa, show patterns consistent with biogeographical regionalization.

Location

Within halite (NaCl) crystals from coastal solar salterns of western Europe, the Mediterranean and east Africa.

Time period

Modern (2006–2013).

Major taxa studied

Archaea.

Methods

Using high‐throughput Illumina amplicon sequencing, we generated the most high‐resolution characterization of halite‐associated archaeal communities to date, using samples from 17 locations. We grouped communities into biogeographical clusters based on community turnover to test whether these communities show biogeographical regionalization. To examine whether individual taxa, rather than communities, show biogeographical patterns, we also tested whether the relative abundance of individual genera may be indicative of a community's biogeographical origins using machine learning methods, specifically random forest classification.

Results

We found that the rate of community turnover was greatest over subregional spatial scales (< 500 km), whereas at regional spatial scales the turnover was independent of geographical distance. Biogeographical clusters of communities were either not statistically robust or lacked spatial coherence, inconsistent with biogeographical regionalization. However, we identified several archaeal genera that were good indicators of biogeographical origin, providing classification error rates of < 10%.

Main conclusions

Overall, our results provide little support for the concept of biogeographical regions in these extremophilic microbial communities, despite the fact that some taxa do show biogeographical patterns. We suggest that variable dispersal ability among the halite‐associated Archaea may disrupt biogeographical patterns at the community level, preventing the formation of biogeographical regions. This means that the processes that lead to the formation of biogeographical regions operate differentially on individual microbial taxa rather than on entire communities.     

The biogeography of lower Mesoamerican freshwater fishes

Aim This paper examines the importance of regional processes in determining the patterns of distribution and diversity of lower Mesoamerican freshwater fishes. Location We focused our analyses on the lower Mesoamerican region, which we define to include all the rivers of Panama and Costa Rica. The geographic boundaries are the Colombian Choco to the south and Lake Nicaragua to the north. Methods We described the biogeographical provinces of lower Mesoamerica (LMA) using presence/absence data of primary and secondary LMA freshwater fishes. We conducted subsequent analyses at the spatial resolution of the biogeographical provinces and described patterns of community composition, species richness, endemism, range size, and the permeability of dispersal barriers between biogeographical provinces. Results This study represents the first attempt since that of W. A. Bussing in 1976 to investigate the biogeographical regions of Mesoamerica, and our analyses demonstrate increased regional complexity in biodiversity patterns relative to previous studies. Changes in community composition across LMA clearly highlight the importance of both extrinsic geological processes and intrinsic biological differences among freshwater fish species in shaping the dispersal and diversification histories of the LMA freshwater fish fauna. The influence of biology and geology is also exemplified by patterns of endemism and turnover between biogeographical provinces, which suggests that the relative importance of regional speciation and dispersal varies spatially across the LMA landscape. Finally, it would seem to follow that secondary freshwater fishes will have larger range sizes than primary fishes as a result of the increased salinity tolerance posited for the former group, and thus the increased probability of dispersal along coastlines. We did not, however, find a significant difference between the average range size of primary and secondary freshwater fishes, indicating that the putative differences in physiological tolerance to seawater between the two groups are not reflected in their distribution patterns at the scale of LMA. The geometric distribution of range size of LMA freshwater fishes suggests that dispersal of both primary and secondary freshwater fishes along coastlines must be infrequent. Main conclusion The observation that regional processes exerted a strong influence on the assembly and maintenance of LMA freshwater fish communities has important consequences for both theory and conservation. We suggest that large‐scale biogeographical analyses are required to illuminate the backdrop upon which local interactions play themselves out, supporting a top‐down approach to the study of biological diversity. Our results also identify areas of high conservation priority, providing a baseline for informing conservation strategies for freshwater fishes in LMA. We conclude by calling for conservation planning and action that acknowledges the importance that regional processes play in determining patterns of organismal diversity, and that incorporates these processes in strategies to conserve remnant biological diversity.     

Biogeographically dynamic genetic structure bridging two continents in the monotypic Central American rodent Ototylomys phyllotis

Central America is an ideal region in which to study patterns of historical divergence and population genetic differentiation, because of its extraordinarily dynamic biogeographical, tectonic, and climatic history. The rodent Ototylomys phyllotis is the only extant species of the genus Ototylomys and is distributed within this region from the Isthmus of Tehuantepec, Mexico, to central Costa Rica, offering an excellent opportunity to study spatial and temporal patterns of population structure of the species and to explain the ecological and evolutionary processes responsible for those patterns. We estimated the genetic diversity and structure within and between populations of O. phyllotis, times of divergence, and migration patterns using mitochondrial DNA and a comprehensive combination of phylogenetic and phylogeographical computational analyses. Our results support monophyly of the genus Ototylomys. We identified three major phylogeographical lineages within O. phyllotis that are linked to its diversification and coincide with the main geological features that shaped Middle America. The origin of the genus was before 3.35 Mya, prior to the Great American Biotic Interchange (GABI), and its initial occurrence was near the centre of its current distribution (Honduras/El Salvador), from which it later spread (3.20–2.84 Mya) following a series of GABI pulses. The species showed an initial northward dispersal to the Chiapas and Guatemala highlands (2.27 Mya) followed by diversification. A later dispersal (1.82 Mya) occurred toward both the south (Nicaragua, Costa Rica) and the north (Belize). The Yucatan peninsula was colonized (0.8 Mya) by individuals from Belize. Extremely high radiation and range expansion occurred throughout the entire range, the highest of which was in the Yucatan peninsula (0.125 Mya). © 2012 The Linnean Society of London, Biological Journal of the Linnean Society, 2012, 107 , 593–610.     

Discordant patterns of evolutionary differentiation in two Neotropical treefrogs

Comparative studies of codistributed taxa test the degree to which historical processes have shaped contemporary population structure. Discordant patterns of lineage divergence among taxa indicate that species differ in their response to common historical processes. The complex geologic landscape of the Isthmus of Central America provides an ideal setting to test the effects of vicariance and other biogeographic factors on population history. We compared divergence patterns between two codistributed Neotropical frogs ( Dendropsophus ebraccatus and Agalychnis callidryas ) that exhibit colour pattern polymorphisms among populations, and found significant differences between them in phenotypic and genetic divergence among populations. Colour pattern in D. ebraccatus did not vary with genetic or geographic distance, while colour pattern co-varied with patterns of gene flow in A. callidryas . In addition, we detected significant species differences in the phylogenetic history of populations, gene flow among them, and the extent to which historical diversification and recent gene flow have been restricted by five biogeographic barriers in Costa Rica and Panama. We inferred that alternate microevolutionary processes explain the unique patterns of diversification in each taxon. Our study underscores how differences in selective regimes and species-typical ecological and life-history traits maintain spatial patterns of diversification.     

Biogeography of the livebearing fish Poecilia gillii in Costa Rica: are phylogeographical breaks congruent with fish community boundaries?

One of the original goals of phylogeography was to use genetic data to identify historical events that might contribute to breaks among communities. In this study, we examine the phylogeography of a common livebearing fish ( Poecilia gillii ) from Costa Rica. Our goal was to determine if phylogeographical breaks in this species were congruent with previously defined boundaries among four fish community provinces. We hypothesized that if abiotic factors influence both community boundaries and genetic structuring in P. gillii then we might find four clades within our focal species that were geographically separated along community boundary lines. Similarly, we expected to find most of the genetic variation in P. gillii partitioned among these four geographical regions. We generated DNA sequence data (mitochondrial cytochrome b and nuclear S7 small ribosomal subunit) for 260 individuals from 42 populations distributed across Costa Rica. We analysed these data using phylogenetic (parsimony and likelihood) and coalescent approaches to estimate phylogenetic relationships among haplotypes, patterns of gene flow and effective population size . Contrary to our expectations, we did not find four monophyletic groups that mapped cleanly to our geographical community provinces. However, one of our clades was restricted to a single province, suggesting that common earth history events could be responsible for both genetic structuring in P. gillii and fish community composition in this area. However, our results show a complex pattern of gene flow throughout other regions in Costa Rica where genetic structuring is not predicted by community province boundaries.     

Intraspecific application of the mid-domain effect model: spatial and temporal nest distributions of green turtles, Chelonia mydas, at Tortuguero, Costa Rica

The mid‐domain effect (MDE) model was developed to evaluate patterns of species richness. We applied the MDE model to intraspecific distribution patterns – the spatial and temporal nest distributions of green turtles, Chelonia mydas, at Tortuguero, Costa Rica, from 1972 to 2000. Spatial and temporal distributions of green turtle nests at Tortuguero did not exhibit significant annual variation over this time period. The spatial and temporal distribution of nests largely conformed to the predictions of the MDE model, although the spatial model has a better fit. Environmental factors that may cause deviations from the MDE model are discussed. The model also indirectly provided a first estimate of the mean spatial nesting range of individual green turtles at Tortuguero: 10.1 km (SD 8.7 km). The MDE model provides insight into intraspecific as well as interspecific distribution patterns.     

Evolution of colour patterns in East African cichlid fish

African cichlid fishes have undergone outbursts of explosive speciation in several lakes, accompanied by rapid radiations in coloration and ecology. Little is known about the evolutionary forces that triggered these events but a hypothesis, published by Wallace Dominey in 1984, has figured prominently. It states that the evolution of colour patterns is driven by sexual selection and that these colour patterns are important in interspecific mate choice, a combination which holds the potential for rapid speciation. Here we present phylogenetic analyses that describe major events in colour evolution and test predictions yielded by Dominey's hypothesis. We assembled information on stripe patterns and the presence or absence of nuptial coloration from more than 700 cichlid species representing more than 90 taxa for which molecular phylogenetic hypotheses were available. We show that sexual selection is most likely the selection force that made male nuptial coloration arise and evolve quickly. In contrast, stripe patterns, though phylogenetically not conserved either, are constrained ecologically. The evolution of vertical bar patterns is associated with structurally complex habitats, such as rocky substrates or vegetation. The evolution of a horizontal stripe is associated with a piscivorous feeding mode. Horizontal stripes are also associated with shoaling behaviour. Strength of sexual selection, measured in terms of the mating system (weak in monogamous, strong in promiscuous species), has no detectable effects on stripe pattern evolution. In promiscuous species the frequency of difference between sister species in nuptial hue is higher than in pair bonding and harem forming species, but the frequency of difference in stripe pattern is lower. We argue that differences between the two components of coloration in their exposure to natural selection explain their very different evolutionary behaviour. Finally, we suggest that habitat-mediated selection upon chromomotor flexibility, a special form of phenotypic plasticity found in the river-dwelling outgroups of the lake-dwelling cichlids, explains the rapid and recurrent ecology-associated radiation of stripe patterns in lake environments, a new hypothesis that yields experimentally testable predictions.     

Species distribution models predict the geographic expansion of an enzootic amphibian pathogen

Globally, numerous amphibian species have declined due to the introduction of the chytrid fungus Batrachochytrium dendrobatidis (Bd). However, the understanding of the spatiotemporal dynamics remains incomplete. Therefore, estimating the current geographic distribution of Bd is urgently needed, especially in countries like Costa Rica, where susceptible species are still recovering from Bd‐driven declines. We conducted model tuning and spatial analysis to compare the habitat suitability for epizootic and enzootic Bd in Costa Rica and to identify data‐deficient regions, opportunistic sampling, and Bd hotspots. Our dataset combined two methods of detection (histology and PCR methods) for a total of 451 Bd‐positive records from 34 localities. We found that the distribution of enzootic Bd in Costa Rica increased 60% since previous estimates in the early 2000s and extended to highlands and dry lowlands that were considered unsuitable for Bd. We also found that Bd is common across protected lands (80%) and within the herpetological provinces containing the highest amphibian richness and endemism in Costa Rica. Opportunistic sampling of Bd has focused on sites where epizootics occurred with the strongest intensity, leading to deficient or absent sampling across the Talamanca Range, the Nicoya Peninsula, and the northern lowlands. Our results showed that PCR increased the power of Bd detection in lowlands and favored the identification of Bd hotspots across the Caribbean side of Costa Rica. Our results add to the understanding of disease spread during enzootics and can be used to identify new hotspots for disease to mitigate future outbreaks of this pathogen. Abstract in Spanish is available with online material     

Patterns and processes in complex landscapes: testing alternative biogeographical hypotheses through integrated analysis of phylogeography and community ecology in Hawai'i

The Island of Hawai'i is a dynamic assemblage of five volcanoes with wet forest habitat currently existing in four distinct natural regions that vary in area, age and geographical isolation. In this complex landscape, alternative assumptions of the relative importance of specific habitat characteristics on evolutionary and ecological processes predict strikingly different general patterns of local diversity and regional similarity. In this study, we compare alternative a priori hypotheses against observed patterns within two distinct biological systems and scales: community composition of wet forest vascular plant species and mitochondrial and nuclear genes of Drosophila sproati, a wet‐forest‐restricted endemic. All observed patterns display strong and similar regional structuring, with the greatest local diversity found in Kohala and the windward side of Mauna Loa, the least in Ka'ū and Kona, and a distinctive pattern of regional similarity that probably reflects the historical development of this habitat on the island. These observations largely corroborate a biogeographical model that integrates multiple lines of evidence, including climatic reconstruction, over those relying on single measures, such as current habitat configuration or substrate age. This method of testing alternative hypotheses across biological systems and scales is an innovative approach for understanding complex landscapes and should prove valuable in diverse biogeographical systems.     

Continental and regional ranges of North American mammals: Rapoport's rule in real and null worlds

Aim To assess the relationship between species richness and distribution within regions arranged along a latitudinal gradient we use the North American mammalian fauna as a study case for testing theoretical models. Location North America. Methods We propose a conceptual framework based on a fully stochastic mid‐domain model to explore geographical patterns of range size and species richness that emerge when the size and position of species ranges along a one‐dimensional latitudinal gradient are randomly generated. We also analyse patterns for the mammal fauna of North America by comparing empirical results from a biogeographical data base with predictions based on randomization null models. Results We confirmed the validity of Rapoport's rule for the mammals of North America by documenting gradients in the size of the continental ranges of species. Additionally, we demonstrated gradients of mean regional range size that parallel those of continental range. Our data also demonstrated that mean range size, measured both as a continental or a regional variable, is significantly correlated with the geographical pattern in species richness. All these patterns deviated sharply from null models. Main conclusions Rapoport's statement of an areographic relationship between species distribution and richness is highly relevant in modern discussions about ecological patterns at the geographical scale.