Effective population size as a driver for divergence of an antimicrobial peptide (Hymenoptaecin) in two common European bumblebee species

Social insects are the target of numerous pathogens. This is because the high density of closely‐related individuals frequently interacting with each other enhances the transmission and establishment of pathogens. This high selective pressure results in the rapid evolution of immune genes, which might be counteracted by a reduced effective population size (N_e) lowering the effectiveness of selection. We tested the effect of N_e on the evolutionary rate of an important immune gene for the antimicrobial peptide Hymenoptaecin in two common central European bumblebee species: Bombus terrestris and Bombus lapidarius. Both species are similar in their biology and are expected to be under similar selective pressures because pathogen prevalence does not differ between species. However, previous studies indicated a higher N_e in B. terrestris compared to B. lapidarius. We found high intraspecific variability in the coding sequence but low variability for silent polymorphisms in B. lapidarius. Estimates of long‐ and short‐term N_e were three‐ to four‐fold higher N_e in B. terrestris, although the species did not differ in census population sizes. The difference in N_e might result in less efficient selection and suboptimal adaptation of immune genes (e.g. hymenoptaecin) in B. lapidarius, and thus this species might become less resistant and more tolerant, turning into a superspreader of diseases.     

Bombus terrestris in a mass‐flowering pollinator‐dependent crop: A mutualistic relationship?

Bumblebees (Bombus spp.) rely on an abundant and diverse selection of floral resources to meet their nutritional requirements. In farmed landscapes, mass‐flowering crops can provide an important forage resource for bumblebees, with increased visitation from bumblebees into mass‐flowering crops having an additional benefit to growers who require pollination services. This study explores the mutualistic relationship between Bombus terrestris L. (buff‐tailed bumblebee), a common species in European farmland, and the mass‐flowering crop courgette (Cucurbita pepo L.) to see how effective B. terrestris is at pollinating courgette and in return how courgette may affect B. terrestris colony dynamics. By combining empirical data on nectar and pollen availability with model simulations using the novel bumblebee model Bumble‐BEEHAVE, we were able to quantify and simulate for the first time, the importance of courgette as a mass‐flowering forage resource for bumblebees. Courgette provides vast quantities of nectar to ensure a high visitation rate, which combined with abundant pollen grains, enables B. terrestris to have a high pollination potential. While B. terrestris showed a strong fidelity to courgette flowers for nectar, courgette pollen was not found in any pollen loads from returning foragers. Nonetheless, model simulations showed that early season courgette (nectar) increased the number of hibernating queens, colonies, and adult workers in the modeled landscapes. Synthesis and applications. Courgette has the potential to improve bumblebee population dynamics; however, the lack of evidence of the bees collecting courgette pollen in this study suggests that bees can only benefit from this transient nectar source if alternative floral resources, particularly pollen, are also available to fulfill bees’ nutritional requirements in space and time. Therefore, providing additional forage resources could simultaneously improve pollination services and bumblebee populations.     

Genetic differentiation of continental and island populations of Bombus terrestris (Hymenoptera: Apidae) in Europe

Ten microsatellite loci and a partial sequence of the COII mitochondrial gene were used to investigate genetic differentiation in B. terrestris, a bumble bee of interest for its high-value crop pollination. The analysis included eight populations from the European continent, five from Mediterranean islands (six subspecies altogether) and one from Tenerife (initially described as a colour form of B. terrestris but recently considered as a separate species, B. canariensis). Eight of the 10 microsatellite loci displayed high levels of polymorphism in most populations. In B. terrestris populations, the total number of alleles detected per polymorphic locus ranged from 3 to 16, with observed allelic diversity from 3.8 ± 0.5 to 6.5 ± 1.4 and average calculated heterozygosities from 0.41 ± 0.09 to 0.65 ± 0.07. B. canariensis showed a significantly lower average calculated heterozygosity (0.12 ± 0.08) and observed allelic diversity (1.5 ± 0.04) as compared to both continental and island populations of B. terrestris. No significant differentiation was found among populations of B. terrestris from the European continent. In contrast, island populations were all significantly and most of them strongly differentiated from continental populations. B. terrestris mitochondrial DNA is characterized by a low nucleotide diversity: 0.18%± 0.07%, 0.20%± 0.04% and 0.27%± 0.04% for the continental populations, the island populations and all populations together, respectively. The only haplotype found in the Tenerife population differs by a single nucleotide substitution from the most common continental haplotype of B. terrestris. This situation, identical to that of Tyrrhenian islands populations and quite different from that of B. lucorum (15 substitutions between terrestris and lucorum mtDNA) casts doubts on the species status of B. canariensis. The large genetic distance between the Tenerife and B. terrestris populations estimated from microsatellite data result, most probably, from a severe bottleneck in the Canary island population. Microsatellite and mitochondrial DNA data call for the protection of the island populations of B. terrestris against importation of bumble bees of foreign origin which are used as crop pollinators.     

Elevational partitioning in species distribution,abundance and body size of Australian alpine grasshoppers (Kosciuscola)

Changes in the physical environment with elevation can influence species distributions and their morphological traits. In mountainous regions, steep temperature gradients can result in patterns of ecological partitioning among species that potentially increases their vulnerability to climate change. We collected data on species distributions, relative abundance and body size for three grasshopper species of the genus Kosciuscola (K. usitatus, K. tristis and K. cognatus) at three locations within the mountainous Kosciuszko National Park in Australia (Thredbo, Guthega and Jagungal). All three species showed differences in their distributions according to elevation, with K. usitatus ranging from 1400 to 2000 m asl, K. tristis from 1600 to 2000 m asl and K. cognatus from 1550 to 1900 m asl. Decreasing relative abundance with increasing elevation was found for K. usitatus, but the opposite pattern was found for K. tristis. The relative abundance of K. cognatus did not change with elevation but was negatively correlated with foliage cover. Body size decreased with elevation in both male and female K. usitatus, which was similarly observed in female K. tristis and male K. cognatus. Our results demonstrate spatial partitioning of species distributions and clines in body size in relation to elevational gradients. Species‐specific sensitivities to climatic gradients may help to predict the persistence of this grasshopper assemblage under climate change.     

Fine‐scale spatial genetic structure of common and declining bumble bees across an agricultural landscape

Land‐use changes have threatened populations of many insect pollinators, including bumble bees. Patterns of dispersal and gene flow are key determinants of species' ability to respond to land‐use change, but have been little investigated at a fine scale (<10 km) in bumble bees. Using microsatellite markers, we determined the fine‐scale spatial genetic structure of populations of four common Bombus species (B. terrestris, B. lapidarius, B. pascuorum and B. hortorum) and one declining species (B. ruderatus) in an agricultural landscape in Southern England, UK. The study landscape contained sown flower patches representing agri‐environment options for pollinators. We found that, as expected, the B. ruderatus population was characterized by relatively low heterozygosity, number of alleles and colony density. Across all species, inbreeding was absent or present but weak (F_IS = 0.01–0.02). Using queen genotypes reconstructed from worker sibships and colony locations estimated from the positions of workers within these sibships, we found that significant isolation by distance was absent in B. lapidarius, B. hortorum and B. ruderatus. In B. terrestris and B. pascuorum, it was present but weak; for example, in these two species, expected relatedness of queens founding colonies 1 m apart was 0.02. These results show that bumble bee populations exhibit low levels of spatial genetic structure at fine spatial scales, most likely because of ongoing gene flow via widespread queen dispersal. In addition, the results demonstrate the potential for agri‐environment scheme conservation measures to facilitate fine‐scale gene flow by creating a more even distribution of suitable habitats across landscapes.     

Neuronal Plasticity in the Mushroom‐Body Calyx of Bumble Bee Workers During Early Adult Development

Division of labor among workers is a key feature of social insects and frequently characterized by an age‐related transition between tasks, which is accompanied by considerable structural changes in higher brain centers. Bumble bees (Bombus terrestris), in contrast, exhibit a size‐related rather than an age‐related task allocation, and thus workers may already start foraging at two days of age. We ask how this early behavioral maturation and distinct size variation are represented at the neuronal level and focused our analysis on the mushroom bodies (MBs), brain centers associated with sensory integration, learning and memory. To test for structural neuronal changes related to age, light exposure, and body size, whole‐mount brains of age‐marked workers were dissected for synapsin immunolabeling. MB calyx volumes, densities, and absolute numbers of olfactory and visual projection neuron (PN) boutons were determined by confocal laser scanning microscopy and three‐dimensional image analyses. Dark‐reared bumble bee workers showed an early age‐related volume increase in olfactory and visual calyx subcompartments together with a decrease in PN‐bouton density during the first three days of adult life. A 12:12  h light‐dark cycle did not affect structural organization of the MB calyces compared to dark‐reared individuals. MB calyx volumes and bouton numbers positively correlated with body size, whereas bouton density was lower in larger workers. We conclude that, in comparison to the closely related honey bees, neuronal maturation in bumble bees is completed at a much earlier stage, suggesting a strong correlation between neuronal maturation time and lifestyle in both species.     

Molecular identification of cryptic bumblebee species from degraded samples using PCR–RFLP approach

The worldwide decline and local extinctions of bumblebees have raised a need for fast and accurate tools for species identification. Morphological characters are often not sufficient, and molecular methods have been increasingly used for reliable identification of bumblebee species. Molecular methods often require high‐quality DNA which makes them less suitable for analysis of low‐quality or older samples. We modified the PCR–RFLP protocol for an efficient and cost‐effective identification of four bumblebee species in the subgenus Bombus s. str. (B. lucorum, B. terrestris, B. magnus and B. cryptarum). We used a short partial mitochondrial COI fragment (446 bp) and three diagnostic restriction enzymes (Hinf I, Hinc II and Hae III) to identify species from degraded DNA material. This approach allowed us to efficiently determine the correct species from all degraded DNA samples, while only a subset of samples 64.6% (31 of 48) resulted in successful amplification of a longer COI fragment (1064 bp) using the previously described method. This protocol can be applied for conservation and management of bumblebees within this subgenus and is especially useful for fast species identification from degraded samples.     

Dietary opportunism,resource partitioning,and consumption of coffee berry borers by five species of migratory wood warblers (Parulidae) wintering in Jamaican shade coffee plantations

Diets reflect important ecological interactions, but are challenging to quantify for foliage‐gleaning birds. We used regurgitated stomach samples from five primarily insectivorous species of long‐distance migrant warblers (Parulidae) wintering in two moderate‐elevation shade coffee farms in Jamaica to assess both foraging opportunism and prey resource partitioning. Our results, based primarily on 6120 prey items in 80 stomach samples collected during a one‐week period in March 2000, confirm opportunism. The diets of all five warblers, including American Redstarts (Setophaga ruticilla), Black‐and‐White Warblers (Mniotilta varia), Black‐throated Blue Warblers (S. caerulescens), Northern Parulas (S. americana), and Prairie Warblers (S. discolor), overlapped strongly based on consumption of the same prey types, even many of the same prey species (4 of 10 interspecific overlaps >0.9, range = 0.74–0.97). Moreover, all five species fed on similarly small, often patchily distributed prey, including coffee berry borers (Hypothenemus hampei; Coleoptera, Curculionidae). Nonetheless, permutational multivariate analysis of variance also revealed that the diets of these species differed significantly, primarily with respect to prey mobility (winged vs. sessile); American Redstarts fed on the most mobile prey, and Northern Parulas on the least mobile prey and a relatively restricted set of prey taxa compared to the other four species of warblers. Overall, our results suggest both dietary opportunism consistent with a migratory life‐history, and interspecific resource partitioning consistent with differences in morphology and foraging behavior during a food‐limited season. Having provided evidence of the three necessary conditions, namely intraspecific competition, resource limitation, and interspecific overlap in resource use, the results of our study, in combination with those of other studies, also provide evidence of interspecific competition among wintering migrant insectivores. We thus argue that diffuse interspecific exploitative food competition may be more important than previously recognized.     

The effect of terrain and female density on survival of neonatal white‐tailed deer and mule deer fawns

Juvenile survival is a highly variable life‐history trait that is critical to population growth. Antipredator tactics, including an animal's use of its physical and social environment, are critical to juvenile survival. Here, we tested the hypothesis that habitat and social characteristics influence coyote (Canis latrans) predation on white‐tailed deer (Odocoileus virginianus) and mule deer (O. hemionus) fawns in similar ways during the neonatal period. This would contrast to winter when the habitat and social characteristics that provide the most safety for each species differ. We monitored seven cohorts of white‐tailed deer and mule deer fawns at a grassland study site in Alberta, Canada. We used logistic regression and a model selection procedure to determine how habitat characteristics, climatic conditions, and female density influenced fawn survival during the first 8 weeks of life. Fawn survival improved after springs with productive vegetation (high integrated Normalized Difference Vegetation Index values). Fawns that used steeper terrain were more likely to survive. Fawns of both species had improved survival in years with higher densities of mule deer females, but not with higher densities of white‐tailed deer females, as predicted if they benefit from protection by mule deer. Our results suggest that topographical variation is a critical resource for neonates of many ungulate species, even species like white‐tailed deer that use more gentle terrain when older. Further, our results raise the possibility that neonatal white‐tailed fawns may benefit from associating with mule deer females, which may contribute to the expansion of white‐tailed deer into areas occupied by mule deer.     

Causal or spurious relationship? Climate and the distribution of Phelsuma geckos on Grand Comoro Island

The volcanic island of Grand Comoro, Malagasy biogeographic region, is inhabited by three species of Phelsuma day geckos; two island‐endemic taxa (Phelsuma comorensis and Phelsuma v‐nigra comoraegrandensis) and the introduced Phelsuma dubia. Phelsuma comorensis is restricted to elevations of greater than 150 m above sea level on the northern of the island's two volcanoes and is the only Phelsuma above 300 m. The other species are widespread at low elevations but also reach levels above 900 m at the southern volcano. To investigate these divergent distribution patterns, we used environmental niche models based on climate and habitat data and tested whether predicted climate change may influence species distributions. Analyses of niche overlap did not show significant differences between present‐day and predicted future potential distributions of any Phelsuma species studied, which could be seen as an indicator of resilience towards climate change. Climate models reflected the restricted distribution of P. comorensis with precipitation of the wettest month detected as most important variable, whereas habitat models predicted an island‐wide distribution. While climate appears to determine the distribution of P. comorensis, we propose isolation by migration barriers as an alternative and discuss the detection of causal versus spurious relationships in ecological niche models.     

Experimental viral spillover can harm Bombus terrestris workers under field conditions

1. Deformed wing virus (DWV), notorious for its virulence in the western honey bee (Apis mellifera) when vectored by the ectoparasitic mite Varroa destructor, is also widespread among wild bumble bee species, presumably through spillover from honey bees. Experimental studies on the virulence of DWV in Bombus spp. have provided equivocal results and have until now been confined to bumble bees under laboratory conditions.
2. Here, we inoculated commercially reared Bombus terrestris workers with DWV-A through feeding or injection and introduced them into experimental colonies placed in the field, thus exposing them to the environment and associated stressors. We monitored the survival of inoculated worker bumble bees and quantified their viral load at 10 days post inoculation.
3. Bombus terrestris workers injected with DWV-A supported high viral loads and exhibited significantly reduced median survival compared to controls. Bumble bees inoculated by feeding had low or zero detectable viral loads while their mortality did not differ from the control group.
4. Our results demonstrate that, although DWV-A is pathogenic for commercial B. terrestris, the risks for individual fitness from spillover of DWV-A during foraging on shared flowers appear limited.
5. The findings of this experiment also highlight the necessity to address the potential context-dependence of virulence when evaluating the impact of a pathogen in an alternative host.

     

American marten and fisher do not segregate in space and time during winter in a mixed‐forest system

Understanding the mechanisms of coexistence between ecologically similar species is an important issue in ecology. Carnivore coexistence may be facilitated by spatial segregation, temporal avoidance, and differential habitat selection. American martens Martes americana and fishers Pekania pennanti are medium‐sized mustelids that occur sympatrically across portions of North America, yet mechanisms of coexistence between the two species are not fully understood. We assessed spatial and temporal partitioning in martens and fishers in the Upper Peninsula of Michigan, USA, using camera trap data collected during winter 2013–2015. To investigate spatial segregation, we used a dynamic occupancy model to estimate species’ occupancy probabilities and probabilities of persistence and colonization as a function of covariates and yearly occupancy probability for the other species. Temporal segregation was assessed by estimating diel activity overlap between species. We found weak evidence of spatial or temporal niche partitioning of martens and fishers. There was high overlap in forest cover selection, and both marten and fisher occupancy were positively correlated with deciduous forests (excluding aspen [Populus tremuloides]). There was strong temporal overlap (; CI = 0.79–0.82) with both species exhibiting largely crepuscular activity patterns. Co‐occurrence of martens and fishers appears to be facilitated by mechanisms not investigated in this study, such as partitioning of snow features or diet. Our results add additional insights into resource partitioning of mesocarnivores, but further research is required to enhance our understanding of mechanisms that facilitate marten and fisher coexistence.     

Dietary overlap and seasonality in three species of mormoopid bats from a tropical dry forest

Competing hypotheses explaining species’ use of resources have been advanced. Resource limitations in habitat and/or food are factors that affect assemblages of species. These limitations could drive the evolution of morphological and/or behavioural specialization, permitting the coexistence of closely related species through resource partitioning and niche differentiation. Alternatively, when resources are unlimited, fluctuations in resources availability will cause concomitant shifts in resource use regardless of species identity. Here, we used next‐generation sequencing to test these hypotheses and characterize the diversity, overlap and seasonal variation in the diet of three species of insectivorous bats of the genus Pteronotus. We identified 465 prey (MOTUs) in the guano of 192 individuals. Lepidoptera and Diptera represented the most consumed insect orders. Diet of bats exhibited a moderate level of overlap, with the highest value between Pteronotus parnellii and Pteronotus personatus in the wet season. We found higher dietary overlap between species during the same seasons than within any single species across seasons. This suggests that diets of the three species are driven more by prey availability than by any particular predator‐specific characteristic. P. davyi and P. personatus increased their dietary breadth during the dry season, whereas P. parnellii diet was broader and had the highest effective number of prey species in all seasons. This supports the existence of dietary flexibility in generalist bats and dietary niche overlapping among groups of closely related species in highly seasonal ecosystems. Moreover, the abundance and availability of insect prey may drive the diet of insectivores.     

Alternate life history phases of a common seaweed have distinct microbial surface communities

Macroalgal life histories are complex, often involving the alternation of distinct free‐living life history phases that differ in morphology, longevity and ploidy. The surfaces of marine macroalgae support diverse microbial biofilms, yet the degree of microbial variation between alternate phases is unknown. We quantified bacterial (16S rRNA gene) and microeukaryote (18S rRNA gene) communities on the surface of the common intertidal seaweed, Mastocarpus spp., which alternates between gametophyte (foliose, haploid) and sporophyte (encrusting, diploid) life history phases. A large portion (97%) of bacterial taxa on the surface Mastocarpus was also present in samples from the environment, indicating that macroalgal surface communities are largely assembled from the surrounding seawater. Still, changes in the relative abundance of bacterial taxa result in significantly different communities on alternate Mastocarpus life history phases, rocky substrate and seawater at all intertidal elevations. For microeukaryote assemblages, only high intertidal samples had significant differences between life history phases although sporophytes were not different from the rocky substrate at this elevation; gametophytes and sporophytes did not differ in microeukaryote communities in the mid and low zones. By sequencing three host genes, we identified three cryptic species of Mastocarpus in our data set, which co‐occur in the mid‐to‐low intertidal zone. In these samples, M. alaskensis sporophytes harboured distinct bacterial communities compared to M. agardhii and M. intermedius sporophytes, which were not distinguishable. Conversely, microeukaryote communities did not differ among species.     

Realized trophic niche driven by apparent competition: an example with marsupials

According to apparent competition theory, the co‐occurrence of two species that share the same predators appears to affect each other's population growth and abundance. However, due to habitat loss and over‐hunting, top predators are being made rare worldwide. Considering that apparent competitors share similar resources, we would expect the absence of top predators to reflect in changes on prey realized trophic niches. To test our hypothesis, we developed a model to predict the abundance ratio of apparent competitor species based on changes in their realized trophic niches. We tested our model against field data on the Neotropical marsupials Didelphis aurita and Metachirus nudicaudatus. Our results revealed that D. aurita and M. nudicaudatus are two species under apparent competition and their realized trophic niche and diet overlap change according to the presence of top predators. The model was able to predict the actual relative abundances of D. aurita and M. nudicaudatus in the three empirical studies analyzed. Our study presents quantitative support to the apparent competition theory; however, the model's applications to other groups still need to be verified. Additionally, our study shows that the lack of top predators has consequences on the realized trophic niche of their prey, and therefore, we reinforce that conservation plans need to focus on the effects of top predator loss on ecosystems.     

Population genomics reveals a candidate gene involved in bumble bee pigmentation

Variation in bumble bee color patterns is well‐documented within and between species. Identifying the genetic mechanisms underlying such variation may be useful in revealing evolutionary forces shaping rapid phenotypic diversification. The widespread North American species Bombus bifarius exhibits regional variation in abdominal color forms, ranging from red‐banded to black‐banded phenotypes and including geographically and phenotypically intermediate forms. Identifying genomic regions linked to this variation has been complicated by strong, near species level, genome‐wide differentiation between red‐ and black‐banded forms. Here, we instead focus on the closely related black‐banded and intermediate forms that both belong to the subspecies B. bifarius nearcticus. We analyze an RNA sequencing (RNAseq) data set and identify a cluster of single nucleotide polymorphisms (SNPs) within one gene, Xanthine dehydrogenase/oxidase‐like, that exhibit highly unusual differentiation compared to the rest of the sequenced genome. Homologs of this gene contribute to pigmentation in other insects, and results thus represent a strong candidate for investigating the genetic basis of pigment variation in B. bifarius and other bumble bee mimicry complexes.     

Apolygus lucorum genome provides insights into omnivorousness and mesophyll feeding

Apolygus lucorum (Miridae) is an omnivorous pest that occurs worldwide and is notorious for the serious damage it causes to various crops and substantial economic losses. Although some studies have examined the biological characteristics of the mirid bug, no reference genome is available in Miridae, limiting in‐depth studies of this pest. Here, we present a chromosome‐scale reference genome of A. lucorum, the first sequenced Miridae species. The assembled genome size was 1.02 Gb with a contig N50 of 785 kb. With Hi‐C scaffolding, 1,016 Mb contig sequences were clustered, ordered and assembled into 17 large scaffolds with scaffold N50 length 68 Mb, each corresponding to a natural chromosome. Numerous transposable elements occur in this genome and contribute to the large genome size. Expansions of genes associated with omnivorousness and mesophyll feeding such as those related to digestion, chemosensory perception, and detoxification were observed in A. lucorum, suggesting that gene expansion contributed to its strong environmental adaptability and severe harm to crops. We clarified that a salivary enzyme polygalacturonase is unique in mirid bugs and has significantly expanded in A. lucorum, which may contribute to leaf damage from this pest. The reference genome of A. lucorum not only facilitates biological studies of Hemiptera as well as an understanding of the damage mechanism of mesophyll feeding, but also provides a basis on which to develop efficient control technologies for mirid bugs.     

Cryptic species diversity in a widespread bumble bee complex revealed using mitochondrial DNA RFLPs

Cryptic species diversity is thought to be common within the class Insecta, posing problems for basic ecological and population genetic studies and conservation management. Within the temperate bumble bee (Bombus spp.) fauna, members of the subgenus Bombus sensu stricto are amongst the most abundant and widespread. However, their species diversity is controversial due to the extreme difficulty or inability morphologically to identify the majority of individuals to species. Our character-based phylogenetic analyses of partial CO1 (700 bp) from 39 individuals spread across their sympatric European ranges provided unequivocal support for five taxa (3–22 diagnostic DNA base pair sites per species). Inclusion of 20 Irish specimens to the dataset revealed ≥2.3% sequence divergence between taxa and ≤1.3% within taxa. We developed a PCR-RFLP based method for unequivocally distinguishing amongst the four cryptic European taxa of this subgenus, B. cryptarum, B. lucorum, B. magnus and B. terrestris, and used it to analyse 391 females of the former three species collected across Ireland, all of which could be unambiguously assigned to species. Bombus lucorum was the most widely distributed and abundant of the cryptarum–lucorum–magnus species complex, comprising 56% of individuals, though it was significantly less abundant at higher altitudes (>200 m) whilst B. cryptarum was relatively more abundant at higher altitudes. Bombus magnus was rarely encountered at urban sites. Both B. lucorum and B. terrestris are nowadays reared commercially for pollination and transported globally. Our RFLP approach to identify native fauna can underpin ecological studies of these important cryptic species as well as the impact of commercial bumble bees on them.     

Remarkable life history polymorphism may be evolving under divergent selection in the silverleaf sunflower

Substantial intraspecific variation in life history is rare and potentially a signal of incipient ecological speciation, if variation is driven by geographically heterogenous natural selection. We present the first report of extensive life history polymorphism in Helianthus argophyllus, the silverleaf sunflower, and examine evidence for its evolution by divergent selection. In 18 populations sampled from across the species range and grown in a common garden, most quantitative traits covaried such that individuals could be assigned to two distinct life history syndromes: tall and late flowering with small initial flowerheads, or short and early flowering with larger initial flowerheads. Helianthus argophyllus exhibits regional genetic structure, but this population structure does not closely correspond with patterns of phenotypic variation. The early‐flowering syndrome is primarily observed in populations from coastal barrier islands, while populations from the nearby mainland coast, although geographically and genetically close, are primarily late flowering. Additionally, several traits are more differentiated among regions than expected based on neutral genetic divergence (Q_ST > F_ST), including the first principal component score corresponding with life history syndrome. This discordance between patterns of phenotypic and genetic variation suggests that divergent selection is driving genetic differences in life history across the species range. If so, the silverleaf sunflower may be in early stages of ecological speciation.