Decoupled dimensions of leaf economic and anti-herbivore defense strategies in a tropical canopy tree community.

Trade-offs among plant functional traits indicate diversity in plant strategies of growth and survival. The leaf economics spectrum (LES) reflects a trade-off between short-term carbon gain and long-term leaf persistence. A related trade-off, between foliar growth and anti-herbivore defense, occurs among plants growing in contrasting resource regimes, but it is unclear whether this trade-off is maintained within plant communities, where resource gradients are minimized. The LES and the growth-defense trade-off involve related traits, but the extent to which these trade-off dimensions are correlated is poorly understood. We assessed the relationship between leaf economic and anti-herbivore defense traits among sunlit foliage of 345 canopy trees in 83 species on Barro Colorado Island, Panama. We quantified ten traits related to resource allocation and defense, and identified patterns of trait co-variation using multivariate ordination. We tested whether traits and ordination axes were correlated with patterns of phylogenetic relatedness, juvenile demographic trade-offs, or topo-edaphic variation. Two independent axes described ~ 60% of the variation among canopy trees. Axis 1 revealed a trade-off between leaf nutritional and structural investment, consistent with the LES. Physical defense traits were largely oriented along this axis. Axis 2 revealed a trade-off between investments in phenolic defenses versus other foliar defenses, which we term the leaf defense spectrum. Phylogenetic relationships and topo-edaphic variation largely did not explain trait co-variation. Our results suggest that some trade-offs among the growth and defense traits of outer-canopy trees may be captured by the LES, while others may occur along additional resource allocation dimensions.

Persistence of Neighborhood Demographic Influences over Long Phylogenetic Distances May Help Drive Post-Speciation Adaptation in Tropical Forests.

Studies of forest dynamics plots (FDPs) have revealed a variety of negative density-dependent (NDD) demographic interactions, especially among conspecific trees. These interactions can affect growth rate, recruitment and mortality, and they play a central role in the maintenance of species diversity in these complex ecosystems. Here we use an equal area annulus (EAA) point-pattern method to comprehensively analyze data from two tropical FDPs, Barro Colorado Island in Panama and Sinharaja in Sri Lanka. We show that these NDD interactions also influence the continued evolutionary diversification of even distantly related tree species in these FDPs. We examine the details of a wide range of these interactions between individual trees and the trees that surround them. All these interactions, and their cumulative effects, are strongest among conspecific focal and surrounding tree species in both FDPs. They diminish in magnitude with increasing phylogenetic distance between heterospecific focal and surrounding trees, but do not disappear or change the pattern of their dependence on size, density, frequency or physical distance even among the most distantly related trees. The phylogenetic persistence of all these effects provides evidence that interactions between tree species that share an ecosystem may continue to promote adaptive divergence even after the species' gene pools have become separated. Adaptive divergence among taxa would operate in stark contrast to an alternative possibility that has previously been suggested, that distantly related species with dispersal-limited distributions and confronted with unpredictable neighbors will tend to converge on common strategies of resource use. In addition, we have also uncovered a positive density-dependent effect: growth rates of large trees are boosted in the presence of a smaller basal area of surrounding trees. We also show that many of the NDD interactions switch sign rapidly as focal trees grow in size, and that their cumulative effect can strongly influence the distributions and species composition of the trees that surround the focal trees during the focal trees' lifetimes.

Patterns of Phylogenetic Diversity of Subtropical Rainforest of the Great Sandy Region, Australia Indicate Long Term Climatic Refugia.

Australia's Great Sandy Region is of international significance containing two World Heritage areas and patches of rainforest growing on white sand. Previous broad-scale analysis found the Great Sandy biogeographic subregion contained a significantly more phylogenetically even subset of species than expected by chance contrasting with rainforest on white sand in Peru. This study aimed to test the patterns of rainforest diversity and relatedness at a finer scale and to investigate why we may find different patterns of phylogenetic evenness compared with rainforests on white sands in other parts of the world. This study focussed on rainforest sites within the Great Sandy and surrounding areas in South East Queensland (SEQ), Australia. We undertook field collections, expanded our three-marker DNA barcode library of SEQ rainforest plants and updated the phylogeny to 95% of the SEQ rainforest flora. We sampled species composition of rainforest in fixed area plots from 100 sites. We calculated phylogenetic diversity (PD) measures as well as species richness (SR) for each rainforest community. These combined with site variables such as geology, were used to evaluate patterns and relatedness. We found that many rainforest communities in the Great Sandy area were significantly phylogenetically even at the individual site level consistent with a broader subregion analysis. Sites from adjacent areas were either not significant or were significantly phylogenetically clustered. Some results in the neighbouring areas were consistent with historic range expansions. In contrast with expectations, sites located on the oldest substrates had significantly lower phylogenetic diversity (PD). Fraser Island was once connected to mainland Australia, our results are consistent with a region geologically old enough to have continuously supported rainforest in refugia. The interface of tropical and temperate floras in part also explains the significant phylogenetic evenness and higher than expected phylogenetic diversity.

Limited tolerance by insects to high temperatures across tropical elevational gradients and the implications of global warming for extinction.

The critical thermal maximum (CTmax), the temperature at which motor control is lost in animals, has the potential to determine if species will tolerate global warming. For insects, tolerance to high temperatures decreases with latitude, suggesting that similar patterns may exist along elevational gradients as well. This study explored how CTmax varies among species and populations of a group of diverse tropical insect herbivores, the rolled-leaf beetles, across both broad and narrow elevational gradients. Data from 6,948 field observations and 8,700 museum specimens were used to map the elevational distributions of rolled-leaf beetles on two mountains in Costa Rica. CTmax was determined for 1,252 individual beetles representing all populations across the gradients. Initial morphological identifications suggested a total of 26 species with populations at different elevations displaying contrasting upper thermal limits. However, compared with morphological identifications, DNA barcodes (cytochrome oxidase I) revealed significant cryptic species diversity. DNA barcodes identified 42 species and haplotypes across 11 species complexes. These 42 species displayed much narrower elevational distributions and values of CTmax than the 26 morphologically defined species. In general, species found at middle elevations and on mountaintops are less tolerant to high temperatures than species restricted to lowland habitats. Species with broad elevational distributions display high CTmax throughout their ranges. We found no significant phylogenetic signal in CTmax, geography, or elevational range. The narrow variance in CTmax values for most rolled-leaf beetles, especially high-elevation species, suggests that the risk of extinction of insects may be substantial under some projected rates of global warming.

A new species of bromeliad-feeding Cephaloleia Chevrolat (Coleoptera, Chrysomelidae, Cassidinae) from Costa Rica: evidence from DNA barcodes, larval and adult morphology and insect diets.

The Neotropical genus Cephaloleia Chevrolat (Coleoptera: Chrysomelidae: Cassidinae) includes 214 species distributed from the south of Mexico to Argentina. Cephaloleia beetles feed mostly on plants from the order Zingiberales. The interactions between Cephaloleia beetles and their Zingiberales host plants is proposed as one of the oldest and most conservative associations. Here we describe a new species of Cephaloleia (Cephaloleiakuprewiczae sp. n.) that feeds on two species of bromeliads (Pitcairniaarcuata and Pitcairniabrittoniana, Bromeliaceae: Pitcairnioideae). Cephaloleiakuprewiczae was previously described as Cephaloleiahistrionica. This study includes evidence from DNA barcodes (COI), larval and adult morphology and insect diets that separates Cephaloleiakuprewiczae from Cephaloleiahistrionica as a new species.

A well-resolved phylogeny of the trees of Puerto Rico based on DNA barcode sequence data.

BACKGROUND: The use of phylogenetic information in community ecology and conservation has grown in recent years. Two key issues for community phylogenetics studies, however, are (i) low terminal phylogenetic resolution and (ii) arbitrarily defined species pools. METHODOLOGY/PRINCIPAL FINDINGS: We used three DNA barcodes (plastid DNA regions rbcL, matK, and trnH-psbA) to infer a phylogeny for 527 native and naturalized trees of Puerto Rico, representing the vast majority of the entire tree flora of the island (89%). We used a maximum likelihood (ML) approach with and without a constraint tree that enforced monophyly of recognized plant orders. Based on 50% consensus trees, the ML analyses improved phylogenetic resolution relative to a comparable phylogeny generated with Phylomatic (proportion of internal nodes resolved: constrained ML = 74%, unconstrained ML = 68%, Phylomatic = 52%). We quantified the phylogenetic composition of 15 protected forests in Puerto Rico using the constrained ML and Phylomatic phylogenies. We found some evidence that tree communities in areas of high water stress were relatively phylogenetically clustered. Reducing the scale at which the species pool was defined (from island to soil types) changed some of our results depending on which phylogeny (ML vs. Phylomatic) was used. Overall, the increased terminal resolution provided by the ML phylogeny revealed additional patterns that were not observed with a less-resolved phylogeny. CONCLUSIONS/SIGNIFICANCE: With the DNA barcode phylogeny presented here (based on an island-wide species pool), we show that a more fully resolved phylogeny increases power to detect nonrandom patterns of community composition in several Puerto Rican tree communities. Especially if combined with additional information on species functional traits and geographic distributions, this phylogeny will (i) facilitate stronger inferences about the role of historical processes in governing the assembly and composition of Puerto Rican forests, (ii) provide insight into Caribbean biogeography, and (iii) aid in incorporating evolutionary history into conservation planning.

Tropical forest fragmentation limits pollination of a keystone understory herb.

Loss of native vegetation cover is thought to be a major driver of declines in pollination success worldwide. However, it is not well known whether reducing the fragmentation of remaining vegetation can ameliorate these negative effects. We tested the independent effects of composition vs. configuration on the reproductive success of a keystone tropical forest herb (Heliconia tortuosa). To do this we designed a large-scale mensurative experiment that independently varied connected forest-patch size (configuration) and surrounding amount of forest (composition). In each patch, we tested whether pollen tubes, fruit, and seed set were associated with these landscape variables. We also captured hummingbirds as an indication of pollinator availability in a subset of patches according to the same design. We found evidence for an effect of configuration on seed set of H. tortuosa, but not on other aspects of plant reproduction; proportion of seeds produced increased 40% across the gradient in patch size we observed (0.64 to > 1300 ha), independent of the amount of forest in the surrounding landscape at both local and landscape scales. We also found that the availability of pollinators was dependent upon forest configuration; hummingbird capture rates increased three and one-half times across the patch size gradient, independent of forest amount. Finally, pollinator availability was strongly positively correlated with seed set. We hypothesize that the effects of configuration on plant fitness that we observed are due to reduced pollen quality resulting from altered hummingbird availability and/or movement behavior. Our results suggest that prioritizing larger patches of tropical forest may be particularly important for conservation of this species.

Wood nitrogen concentrations in tropical trees: phylogenetic patterns and ecological correlates.

In tropical and temperate trees, wood chemical traits are hypothesized to covary with species' life-history strategy along a 'wood economics spectrum' (WES), but evidence supporting these expected patterns remains scarce. Due to its role in nutrient storage, we hypothesize that wood nitrogen (N) concentration will covary along the WES, being higher in slow-growing species with high wood density (WD), and lower in fast-growing species with low WD. In order to test this hypothesis we quantified wood N concentrations in 59 Panamanian hardwood species, and used this dataset to examine ecological correlates and phylogenetic patterns of wood N. Wood N varied > 14-fold among species between 0.04 and 0.59%; closely related species were more similar in wood N than expected by chance. Wood N was positively correlated with WD, and negatively correlated with log-transformed relative growth rates, although these relationships were relatively weak. We found evidence for co-evolution between wood N and both WD and log-transformed mortality rates. Our study provides evidence that wood N covaries with tree life-history parameters, and that these patterns consistently co-evolve in tropical hardwoods. These results provide some support for the hypothesized WES, and suggest that wood is an increasingly important N pool through tropical forest succession.

Low plant density enhances gene dispersal in the Amazonian understory herb Heliconia acuminata.

In theory, conservation genetics predicts that forest fragmentation will reduce gene dispersal, but in practice, genetic and ecological processes are also dependent on other population characteristics. We used Bayesian genetic analyses to characterize parentage and propagule dispersal in Heliconia acuminata L. C. Richard (Heliconiaceae), a common Amazonian understory plant that is pollinated and dispersed by birds. We studied these processes in two continuous forest sites and three 1-ha fragments in Brazil's Biological Dynamics of Forest Fragments Project. These sites showed variation in the density of H. acuminata. Ten microsatellite markers were used to genotype flowering adults and seedling recruits and to quantify realized pollen and seed dispersal distances, immigration of propagules from outside populations, and reproductive dominance among parents. We tested whether gene dispersal is more dependent on fragmentation or density of reproductive plants. Low plant densities were associated with elevated immigration rates and greater propagule dispersal distances. Reproductive dominance among inside-plot parents was higher for low-density than for high-density populations. Elevated local flower and fruit availability is probably leading to spatially more proximal bird foraging and propagule dispersal in areas with high density of reproductive plants. Nevertheless, genetic diversity, inbreeding coefficients and fine-scale spatial genetic structure were similar across populations, despite differences in gene dispersal. This result may indicate that the opposing processes of longer dispersal events in low-density populations vs. higher diversity of contributing parents in high-density populations balance the resulting genetic outcomes and prevent genetic erosion in small populations and fragments.

The roots of diversity: below ground species richness and rooting distributions in a tropical forest revealed by DNA barcodes and inverse modeling.

BACKGROUND: Plants interact with each other, nutrients, and microbial communities in soils through extensive root networks. Understanding these below ground interactions has been difficult in natural systems, particularly those with high plant species diversity where morphological identification of fine roots is difficult. We combine DNA-based root identification with a DNA barcode database and above ground stem locations in a floristically diverse lowland tropical wet forest on Barro Colorado Island, Panama, where all trees and lianas >1 cm diameter have been mapped to investigate richness patterns below ground and model rooting distributions. METHODOLOGY/PRINCIPAL FINDINGS: DNA barcode loci, particularly the cpDNA locus trnH-psba, can be used to identify fine and small coarse roots to species. We recovered 33 species of roots from 117 fragments sequenced from 12 soil cores. Despite limited sampling, we recovered a high proportion of the known species in the focal hectare, representing approximately 14% of the measured woody plant richness. This high value is emphasized by the fact that we would need to sample on average 13 m(2) at the seedling layer and 45 m(2) for woody plants >1 cm diameter to obtain the same number of species above ground. Results from inverse models parameterized with the locations and sizes of adults and the species identifications of roots and sampling locations indicates a high potential for distal underground interactions among plants. CONCLUSIONS: DNA barcoding techniques coupled with modeling approaches should be broadly applicable to studying root distributions in any mapped vegetation plot. We discuss the implications of our results and outline how second-generation sequencing technology and environmental sampling can be combined to increase our understanding of how root distributions influence the potential for plant interactions in natural ecosystems.

Plant-pollinator interactions and floral convergence in two species of Heliconia from the Caribbean Islands.

Variation in interspecific interactions across geographic space is a potential driver of diversification and local adaptation. This study quantitatively examined variation in floral phenotypes and pollinator service of Heliconia bihai and H. caribaea across three Antillean islands. The prediction was that floral characters would correspond to the major pollinators of these species on each island. Analysis of floral phenotypes revealed convergence among species and populations of Heliconia from the Greater Antilles. All populations of H. caribaea were similar, characterized by long nectar chambers and short corolla tubes. In contrast, H. bihai populations were strongly divergent: on Dominica, H. bihai had flowers with short nectar chambers and long corollas, whereas on Hispaniola, H. bihai flowers resembled those of H. caribaea with longer nectar chambers and shorter corolla tubes. Morphological variation in floral traits corresponded with geographic differences or similarities in the major pollinators on each island. The Hispaniolan mango, Anthracothorax dominicus, is the principal pollinator of both H. bihai and H. caribaea on Hispaniola; thus, the similarity of floral phenotypes between Heliconia species suggests parallel selective regimes imposed by the principal pollinator. Likewise, divergence between H. bihai populations from Dominica and Hispaniola corresponded with differences in the pollinators visiting this species on the two islands. The study highlights the putative importance of pollinator-mediated selection as driving floral convergence and the evolution of locally-adapted plant variants across a geographic mosaic of pollinator species.

What makes a leaf tough? Patterns of correlated evolution between leaf toughness traits and demographic rates among 197 shade-tolerant woody species in a neotropical forest.

Slow-growing juveniles of shade-tolerant plant species are predicted to have tough leaves because of the high cost of leaf replacement in shade relative to potential carbon gain. We assessed the degree of correlated evolution among eight traits associated with leaf toughness and the relationships of those traits with the growth and mortality rates of 197 tree and shrub species from the understory of the 50-ha forest dynamics plot on Barro Colorado Island, Panama. Path analysis with phylogenetically independent contrasts revealed that leaves attained material toughness (resistance to fracture per unit fracture area) through increases in tissue density, percent cellulose per unit dry mass, and vein fracture toughness. Lamina density and cellulose content evolved independently and thus represent different paths to material toughness. Structural toughness (resistance to fracture per unit fracture length) depended on material toughness and lamina thickness. Mortality rates of individuals 1-10 cm in stem diameter were negatively correlated with material toughness and lamina density but were independent of structural toughness and cell wall fiber contents. Leaf toughness traits were uncorrelated with relative growth rates. Results imply that material toughness enhances resistance to natural enemies, which increases survival and offsets the biomass allocation cost of producing tough leaves in the shaded understory.

Advances in the use of DNA barcodes to build a community phylogeny for tropical trees in a Puerto Rican forest dynamics plot.

BACKGROUND: Species number, functional traits, and phylogenetic history all contribute to characterizing the biological diversity in plant communities. The phylogenetic component of diversity has been particularly difficult to quantify in species-rich tropical tree assemblages. The compilation of previously published (and often incomplete) data on evolutionary relationships of species into a composite phylogeny of the taxa in a forest, through such programs as Phylomatic, has proven useful in building community phylogenies although often of limited resolution. Recently, DNA barcodes have been used to construct a robust community phylogeny for nearly 300 tree species in a forest dynamics plot in Panama using a supermatrix method. In that study sequence data from three barcode loci were used to generate a well-resolved species-level phylogeny. METHODOLOGY/PRINCIPAL FINDINGS: Here we expand upon this earlier investigation and present results on the use of a phylogenetic constraint tree to generate a community phylogeny for a diverse, tropical forest dynamics plot in Puerto Rico. This enhanced method of phylogenetic reconstruction insures the congruence of the barcode phylogeny with broadly accepted hypotheses on the phylogeny of flowering plants (i.e., APG III) regardless of the number and taxonomic breadth of the taxa sampled. We also compare maximum parsimony versus maximum likelihood estimates of community phylogenetic relationships as well as evaluate the effectiveness of one- versus two- versus three-gene barcodes in resolving community evolutionary history. CONCLUSIONS/SIGNIFICANCE: As first demonstrated in the Panamanian forest dynamics plot, the results for the Puerto Rican plot illustrate that highly resolved phylogenies derived from DNA barcode sequence data combined with a constraint tree based on APG III are particularly useful in comparative analysis of phylogenetic diversity and will enhance research on the interface between community ecology and evolution.

Plant DNA barcodes and a community phylogeny of a tropical forest dynamics plot in Panama.

The assembly of DNA barcode libraries is particularly relevant within species-rich natural communities for which accurate species identifications will enable detailed ecological forensic studies. In addition, well-resolved molecular phylogenies derived from these DNA barcode sequences have the potential to improve investigations of the mechanisms underlying community assembly and functional trait evolution. To date, no studies have effectively applied DNA barcodes sensu strictu in this manner. In this report, we demonstrate that a three-locus DNA barcode when applied to 296 species of woody trees, shrubs, and palms found within the 50-ha Forest Dynamics Plot on Barro Colorado Island (BCI), Panama, resulted in >98% correct identifications. These DNA barcode sequences are also used to reconstruct a robust community phylogeny employing a supermatrix method for 281 of the 296 plant species in the plot. The three-locus barcode data were sufficient to reliably reconstruct evolutionary relationships among the plant taxa in the plot that are congruent with the broadly accepted phylogeny of flowering plants (APG II). Earlier work on the phylogenetic structure of the BCI forest dynamics plot employing less resolved phylogenies reveals significant differences in evolutionary and ecological inferences compared with our data and suggests that unresolved community phylogenies may have increased type I and type II errors. These results illustrate how highly resolved phylogenies based on DNA barcode sequence data will enhance research focused on the interface between community ecology and evolution.

Heliconia acuminata reproductive success is independent of local floral density

Plantas reprodutivas em florestas tropicas são distribuidas em manchas, com algumas em grandes agregações coespecíficas e outras relativamente isoladas. A hipótese é que esta variação na densidade de flores em um local tem um grande efeito no sucesso reprodutivo de plantas, já que indivíduos em agregações maiores poderiam atrair mais polinizadores ou polinizadores de melhor qualidade. Esta hipótese foi testada na Amazônia central com populações da erva de sub-bosque Heliconia acuminata. Foram criadas parcelas em que a densidade de plantas reprodutivas simula a densidade natural, e medimos componentes de sucesso reprodutivo em plantas focais em cada parcela. Não houve diferença significativa entre nenhum dos tratamentos de densidade de flores em termos da taxa de produção de frutos, de produção total de sementes por planta, ou de produção de sementes por fruto. Visitas por polinizadores eram extremamente raras, e muitas plantas não receberam nenhuma visita. Isto poderia ser porque os beija-flores que polinizam Heliconia acuminata são muito ineficientes, mas isto parece não ser o caso. A densidade de flores pode simplesmente estar abaixo do limiar em que efeitos de densidade local são importantes, até mesmo em agregações de maior densidade. Limitação de nutrientes, aborto seletivo de frutos, e reprodução via função masculina e não feminina também poderiam ser responsáveis. Sugere-se que a ausência de efeitos de densidade local na reprodução de plantas pode ser um fenômeno geral nas florestas da Amazônia central, embora experimentos adicionais com outros sistemas são necessários para determinar se esta hipótese é valida.

Heliconia acuminata reproductive success is independent of local floral density

Reproductive plants in tropical forests are patchily distributed, with some in large aggregations of reproductive consepecifics while others are relatively isolated. This variation in floral density is hypothesized to have a major effect on plant reproductive success, since individuals in higher density neighborhoods can attract more or higher quality pollinators. We experimentally tested this hypothesis with populations of the understory herb Heliconia acuminata in central Amazonia. We created replicated plots in which reproductive plant density spanned the range of naturally occurring floral neighborhood size, then measured three surrogates of plant fitness in focal plants in each array. There was no significant difference between any of the three floral neighborhood treatments in total seed production, fruit set, or the number of seeds produced per fruit. Pollinator visitation rates to plants in all treatments were extremely low, with many plants not visited at all during the observation period. This could be because H. acuminata's hummingbird pollinators are unable to find the widely scattered reproductive plants, however this hypothesis appears unlikely. Instead, natural flowering plant densities may simply be below the threshold value at which neighborhood effects become important, even in "high-density" aggregations. Nutrient limitation, selective fruit abortion, and reproduction via male rather than female function may also be playing a role. We argue the absence of neighborhood effects may be a general phenomenon in central Amazonian forests, though additional experiments with other plant-pollinator systems are needed to determine the extent to which this hypothesis is supported.

Plantas reprodutivas em florestas tropicas são distribuidas em manchas, com algumas em grandes agregações coespecíficas e outras relativamente isoladas. A hipótese é que esta variação na densidade de flores em um local tem um grande efeito no sucesso reprodutivo de plantas, já que indivíduos em agregações maiores poderiam atrair mais polinizadores ou polinizadores de melhor qualidade. Esta hipótese foi testada na Amazônia central com populações da erva de sub-bosque Heliconia acuminata. Foram criadas parcelas em que a densidade de plantas reprodutivas simula a densidade natural, e medimos componentes de sucesso reprodutivo em plantas focais em cada parcela. Não houve diferença significativa entre nenhum dos tratamentos de densidade de flores em termos da taxa de produção de frutos, de produção total de sementes por planta, ou de produção de sementes por fruto. Visitas por polinizadores eram extremamente raras, e muitas plantas não receberam nenhuma visita. Isto poderia ser porque os beija-flores que polinizam Heliconia acuminata são muito ineficientes, mas isto parece não ser o caso. A densidade de flores pode simplesmente estar abaixo do limiar em que efeitos de densidade local são importantes, até mesmo em agregações de maior densidade. Limitação de nutrientes, aborto seletivo de frutos, e reprodução via função masculina e não feminina também poderiam ser responsáveis. Sugere-se que a ausência de efeitos de densidade local na reprodução de plantas pode ser um fenômeno geral nas florestas da Amazônia central, embora experimentos adicionais com outros sistemas são necessários para determinar se esta hipótese é valida.

Adaptation in a plant-hummingbird association.

Sexual dimorphism in bill morphology and body size of the Caribbean purple-throated carib hummingbird is associated with a reversal in floral dimorphism of its Heliconia food plants. This hummingbird is the sole pollinator of H. caribaea and H. bihai, with flowers of the former corresponding to the short, straight bills of males, the larger sex, and flowers of the latter corresponding to the long, curved bills of females. On St. Lucia, H. bihai compensates for the rarity of H. caribaea by evolving a second color morph with flowers that match the bills of males, whereas on Dominica, H. caribaea evolves a second color morph with flowers that match the bills of females. The nectar rewards of all Heliconia morphs are consistent with each sex's choice of the morph that corresponds to its bill morphology and energy requirements, supporting the hypothesis that feeding preferences have driven their coadaptation.