Species discrimination in Sisyrinchium (Iridaceae): assessment of DNA barcodes in a taxonomically challenging genus

DNA barcoding aims to develop an efficient tool for species identification based on short and standardized DNA sequences. In this study, the DNA barcode paradigm was tested among the genera of the tribe Sisyrinchieae (Iridoideae). Sisyrinchium, with more than 77% of the species richness in the tribe, is a taxonomically complex genus. A total of 185 samples belonging to 98 species of Sisyrinchium, Olsynium, Orthrosanthus and Solenomelus were tested using matK, trnH‐psbA and internal transcribed spacer (ITS). Candidate DNA barcodes were analysed either as single markers or in combination. Detection of a barcoding gap, similarity‐based methods and tree‐based analyses were used to assess the discrimination efficiency of DNA barcodes. The levels of species identification obtained from plastid barcodes were low and ranged from 17.35% to 20.41% for matK and 5.11% to 7.14% for trnH‐psbA. The ITS provided better results with 30.61–38.78% of species identified. The analyses of the combined data sets did not result in a significant improvement in the discrimination rate. Among the tree‐based methods, the best taxonomic resolution was obtained with Bayesian inference, particularly when the three data sets were combined. The study illustrates the difficulties for DNA barcoding to identify species in evolutionary complex lineages. Plastid markers are not recommended for barcoding Sisyrinchium due to the low discrimination power observed. ITS gave better results and may be used as a starting point for species identification.     

Comparative multilocus phylogeography of two Palaearctic spruce bark beetles: influence of contrasting ecological strategies on genetic variation

While phylogeographic patterns of organisms are often interpreted through past environmental disturbances, mediated by climate changes, and geographic barriers, they may also be strongly influenced by species‐specific traits. To investigate the impact of such traits, we focused on two Eurasian spruce bark beetles that share a similar geographic distribution, but differ in their ecology and reproduction. Ips typographus is an aggressive tree‐killing species characterized by strong dispersal, whereas Dendroctonus micans is a discrete inbreeding species (sib mating is the rule), parasite of living trees and a poor disperser. We compared genetic variation between the two species over both beetles’ entire range in Eurasia with five independent gene fragments, to evaluate whether their intrinsic differences could have an influence over their phylogeographic patterns. We highlighted widely divergent patterns of genetic variation for the two species and argue that the difference is indeed largely compatible with their contrasting dispersal strategies and modes of reproduction. In addition, genetic structure in I. typographus divides European populations in a northern and a southern group, as was previously observed for its host plant, and suggests past allopatric divergence. A long divergence time was estimated between East Asian and other populations of both species, indicating their long‐standing presence in Eurasia, prior to the last glacial maximum. Finally, the strong population structure observed in D. micans for the mitochondrial locus provides insights into the recent colonization history of this species, from its native European range to regions where it was recently introduced.     

Effects of Competition and Temporal Variation on the Evolutionary Potential of Two Native Bunchgrass Species

The capacity of restored plant populations to adapt to new environmental challenges depends on within‐population genetic variation. We examined how much genetic and environmentally based variation for fitness‐associated traits exists within populations of two native grasses commonly used for restoration in California. We were also interested in understanding how phenotypic expression of genetic variation for these traits varies with growth environment. Thirty maternal families of Elymus glaucus (Blue wild rye) and Nassella pulchra (Purple needlegrass) were sampled from both coastal and interior populations and reciprocally transplanted into three replicated common gardens with and without interspecific competition at each site. Reproductive output of families differed both among years and with competition treatments. Phenotypic expression of genetic variation in culm production differed among populations and was very low when families were grown with interspecific competition. Without interspecific competition, the degree of genetic determination peaked in year two in both species (8.4 and 15.1% in E. glaucus and N. pulchra, respectively). Significant genetic differences in reproduction and phenotypic plasticity were found among N. pulchra subpopulations sampled less than 3 km apart, further highlighting the importance of thoroughly sampling available genetic variation in populations used for restoration. The variable and generally low expression of genetic variation indicates that rates of adaptation in restored populations of these native grasses may vary temporally and may be especially slow within competitive environments.     

Large‐scale adaptive divergence in Boechera fecunda,an endangered wild relative of Arabidopsis

Many biological species are threatened with extinction because of a number of factors such as climate change and habitat loss, and their preservation depends on an accurate understanding of the extent of their genetic variability within and among populations. In this study, we assessed the genetic divergence of five quantitative traits in 10 populations of an endangered cruciferous species, Boechera fecunda, found in only several populations in each of two geographic regions (WEST and EAST) in southwestern Montana. We analyzed variation in quantitative traits, neutral molecular markers, and environmental factors and provided evidence that despite the restricted geographical distribution of this species, it exhibits a high level of genetic variation and regional adaptation. Conservation efforts therefore should be directed to the preservation of populations in each of these two regions without attempting transplantation between regions. Heritabilities and genetic coefficients of variation estimated from nested ANOVAs were generally high for leaf and rosette traits, although lower (and not significantly different from 0) for water‐use efficiency. Measures of quantitative genetic differentiation, Q_ST, were calculated for each trait from each pair of populations. For three of the five traits, these values were significantly higher between regions compared with those within regions (after adjustment for neutral genetic variation, F_ST). This suggested that natural selection has played an important role in producing regional divergence in this species. Our analysis also revealed that the B. fecunda populations appear to be locally adapted due, at least in part, to differences in environmental conditions in the EAST and WEST regions.     

Diploid hybrid origin of Hippophaë gyantsensis (Elaeagnaceae) in the western Qinghai–Tibet Plateau

Homoploid hybrid speciation, the origin of a hybrid species without change in chromosome number, is currently considered to be a rare form of speciation. In the present study, we examined the phylogenetic origin of Hippophaë gyantsensis, a diploid species occurring in the western Qinghai–Tibet Plateau. Some of its morphological and molecular traits suggest a close relationship to H. rhamnoides ssp. yunnanensis while others indicate H. neurocarpa. We conducted phylogenetic analyses of sequence data of two maternally inherited chloroplast (cp) DNA fragments and the bi‐parentally inherited nuclear ribosomal internal transcribed spacer (ITS) from 17 populations of H. gyantsensis, 15 populations of H. rhamnoides ssp. yunnanensis and 27 populations of H. neurocarpa across their distributional ranges, and modelled the niche differentiation of the three taxa. Multiple lines of evidence suggested that H. gyantsensis is a morphologically stable, genetically independent and ecologically distinct species. The inconsistent phylogenetic placements of the H. gyantsensis clade that comprised the dominant cpDNA haplotypes and ITS ribotypes suggested a probable diploid hybrid origin from multiple crosses between H. rhamnoides ssp. yunnanensis and H. neurocarpa. This tentative hypothesis is more parsimonious than alternative explanations according to the data available, although more evidence based on further testing is needed.     

Host‐associated genomic differentiation in congeneric butterflies: now you see it,now you do not

Ecotypic variation among populations may become associated with widespread genomic differentiation, but theory predicts that this should happen only under particular conditions of gene flow, selection and population size. In closely related species, we might expect the strength of host‐associated genomic differentiation (HAD) to be correlated with the degree of phenotypic differentiation in host‐adaptive traits. Using microsatellite and Amplified Fragment Length Polymorphism (AFLP) markers, and controlling for isolation by distance between populations, we sought HAD in two congeneric species of butterflies with different degrees of host plant specialization. Prior work on Euphydryas editha had shown strong interpopulation differentiation in host‐adapted traits, resulting in incipient reproductive isolation among host‐associated ecotypes. We show here that Euphydryas aurinia had much weaker host‐associated phenotypic differentiation. Contrary to our expectations, we detected HAD in Euphydryas aurinia, but not in E. editha. Even within an E. aurinia population that fed on both hosts, we found weak but significant sympatric HAD that persisted in samples taken 9 years apart. The finding of significantly stronger HAD in the system with less phenotypic differentiation may seem paradoxical. Our findings can be explained by multiple factors, ranging from differences in dispersal or effective population size, to spatial variation in genomic or phenotypic traits and to structure induced by past histories of host‐adapted populations. Other infrequently measured factors, such as differences in recombination rates, may also play a role. Our result adds to recent work as a further caution against assumptions of simple relationships between genomic and adaptive phenotypic differentiation.     

The comparative reproductive biology of a tetraploid species,Hedychium villosum,and its diploid progenitor H. tenuiflorum (Zingiberaceae)

The evolutionary advantages of polyploidy may result from a number of changes in floral traits and breeding system, which may enable polyploids to exploit new habitats and become widespread. In this study, we comparatively investigated the floral biology of the tetraploid species Hedychium villosum and its diploid progenitor H. tenuiflorum, to assess reproductive divergence between the two species. The results showed that flowers of the tetraploid species last longer and produce more nectar than did diploid species. The flowering times of the two species did not overlap at all. Observations of floral visitors in natural populations demonstrated that butterflies and hawkmoths were effective pollinators of both species, but there was a significant difference in butterfly and hawkmoth assemblages between the two species. The hand‐pollination experiments and pollen tube growth experiments suggested that diploid H. tenuiflorum was self‐incompatible, while tetraploid H. villosum was completely self‐compatible. H. villosum has a much wider distribution range and occupies more diverse habitats than H. tenuiflorum. Polyploidisation may enable tetraploid H. villosum to exploit new habitats previously unavailable to diploid H. tenuiflorum.     

Variation in floral morphology and plant reproductive success in four Ipomoea species (Convolvulaceae) with contrasting breeding systems

• This study tested the hypothesis that self‐compatibility would be associated with floral traits that facilitate autonomous self‐pollination to ensure reproduction under low pollinator visitation. In a comparison of two pairs of Ipomoea species with contrasting breeding systems, we predicted that self‐compatible (SC) species would have smaller, less variable flowers, reduced herkogamy, lower pollinator visitation and higher reproductive success than their self‐incompatible (SI) congeners.
• We studied sympatric species pairs, I. hederacea (SC)– I. mitchellae (SI) and I. purpurea (SC)–I. indica (SI), in Mexico, over two years. We quantified variation in floral traits and nectar production, documented pollinator visitation, and determined natural fruit and seed set. Hand‐pollination and bagging experiments were conducted to determine potential for autonomous self‐pollination and apomixis.
• Self‐compatible Ipomoea species had smaller flowers and lower nectar production than SI species; however, floral variation and integration did not vary according to breeding system. Bees were primary pollinators of all species, but visitation rates were seven times lower in SC than SI species. SC species had a high capacity for autonomous self‐pollination due to reduced herkogamy at the highest anther levels. Self‐compatible species had two to six times higher fruit set than SI species.
• Results generally support the hypothesis that self‐compatibility and autonomous self‐pollination ensure reproduction under low pollinator visitation. However, high variation in morphological traits of SC Ipomoea species suggests they maintain variation through outcrossing. Furthermore, reduced herkogamy was associated with high potential for autonomous self‐pollination, providing a reproductive advantage that possibly underlies transitions to self‐compatibility in Ipomoea.

     

Microgeographic morphological variation across larval wood frog populations associated with environment despite gene flow

Gene flow has historically been thought to constrain local adaptation; yet, recent research suggests that populations can diverge despite exchanging genes. Here I use a common garden experiment to assess the combined effects of gene flow and natural selection on morphological variation of 16 wood frog (Rana sylvatica) populations, a species known to experience divergent selection pressures in open‐ and closed‐canopy ponds across relatively small geographic scales. Wood frog tadpoles from different ponds showed significant morphological variation associated with canopy type with a trade‐off between tail length and body depth consistent with previous research. In contrast, neutral genetic differentiation of nine microsatellite loci as measured by Jost's D was not associated with canopy type, indicating no pattern of isolation by environment. Genetic structure analyses indicated some substructure across the 16 ponds (K = 4); however, three out of four assigned clusters included both open‐ and closed‐canopy ponds. Together, these results suggest that morphological divergence among these wood frog populations is occurring despite gene flow and that selection within these environments is strong. Furthermore, morphological variation among ponds differed across two sampling periods during larval development, demonstrating the importance of evaluating phenotypic divergence over multiple time periods and at a time relevant to the processes being studied.     

Geographic variation of reproductive traits and competition for pollinators in a bird‐pollinated plant

Geographic variation in the reproductive traits of animal‐pollinated plants can be shaped by spatially variable selection imposed by differences in the local pollination environment. We investigated this process in Babiana ringens (Iridaceae), an enigmatic species from the Western Cape region of South Africa. B. ringens has evolved a specialized perch facilitating cross‐pollination by sunbirds and displays striking geographic variation in perch size and floral traits. Here, we investigate whether this variation can be explained by geographic differences in the pollinator communities. We measured floral and inflorescence traits, and abiotic variables (N, P, C, and rainfall) and made observations of sunbirds in populations spanning the range of B. ringens. In each population, we recorded sunbird species identity and measured visitation rates, interfloral pollen transfer, and whether the seed set of flowers was pollen limited. To evaluate whether competition from co‐occurring sunbird‐pollinated species might reduce visitation, we quantified nectar rewards in B. ringens and of other co‐flowering bird‐pollinated species in local communities in which populations occurred. Variation in abiotic variables was not associated with geographical variation of traits in B. ringens. Malachite sunbirds were the dominant visitor (97% of visits) and populations with larger‐sized traits exhibited higher visitation rates, more between‐flower pollen transfer and set more seed. No sunbirds were observed in four populations, all with smaller‐sized traits. Sunbird visitation to B. ringens was not associated with local sunbird activity in communities, but sunbird visitation was negatively associated with the amount of B. ringens sugar relative to the availability of alternative nectar sources. Our study provides evidence that B. ringens populations with larger floral traits are visited more frequently by sunbirds, and we propose that visitation rates to B. ringens may be influenced, in part, by competition with other sunbird‐pollinated species.     

Right here,right now: Populations of Actinotus helianthi differ in their early performance traits and interactions

Populations across the geographical distribution of a species are shaped by different local environments to produce distinctive patterns of variation in plant traits. Among‐population variation is, therefore, important for understanding potential shifts in distributions under changing environments, but is often not included in studies. In particular, critical data on the suitability of local environments for plant traits expressed at different life stages are lacking. To address this we performed two experiments to disentangle the influence of the local environment on multiple plant traits for populations of Actinotus helianthi from across its latitudinal range. A common environment experiment was used to compare early plant traits of germination, early seedling growth and survival for 17 populations of A. helianthi. To examine how biotic interactions vary across populations, we evaluated whether plant traits, including height and number of pseudanthia, influence visitor diversity and abundance, and if insect visitor abundance or diversity was associated with seed set success. We found that populations varied in germination success between 0.2 ± 0.1% and 64.2 ± 2.3%. Seedling growth and early survival varied among populations by as much as a factor of two and 44 respectively. We recorded variation in plant traits across hierarchical spatial scales from the maternal plant to biogeographical regions. The abundance and diversity of insect visitors also varied among populations and seed set was found to be site specific. There was a trend for populations with taller plants and larger floral display sizes to be more frequently visited by pollinators. We also identified a positive linear relationship between the number of visits by flies and seed set success. These results suggest that the local environment has a strong role in directly and indirectly influencing variation in plant traits within populations of A. helianthi, and potentially other perennial species.     

Maintaining their genetic distance: Little evidence for introgression between widely hybridizing species of Geum with contrasting mating systems

Within the plant kingdom, many genera contain sister lineages with contrasting outcrossing and inbreeding mating systems that are known to hybridize. The evolutionary fate of these sister lineages is likely to be influenced by the extent to which they exchange genes. We measured gene flow between outcrossing Geum rivale and selfing Geum urbanum, sister species that hybridize in contemporary populations. We generated and used a draft genome of G. urbanum to develop dd‐RAD data scorable in both species. Coalescent analysis of RAD data from allopatric populations indicated that the species diverged 2–3 Mya, and that historical gene flow between them was extremely low (1 migrant every 25 generations). Comparison of genetic divergence between species in sympatry and allopatry, together with an analysis of allele frequencies in potential parental and hybrid populations, provided no evidence of contemporary introgression in sympatric populations. Cluster‐ and species‐specific marker analyses revealed that, apart from four early‐generation hybrids, individuals in sympatric populations fell into two genetically distinct groups that corresponded exactly to their morphological species classification with maximum individual admixture estimates of only 1–3%. However, we did observe joint segregation of four putatively introgressed SNPs across two scaffolds in the G. urbanum population that was associated with significant morphological variation, interpreted as tentative evidence for rare, recent interspecific gene flow. Overall, our results indicate that despite the presence of hybrids in contemporary populations, genetic exchange between G. rivale and G. urbanum has been extremely limited throughout their evolutionary history.     

Expression of anatomical leaf traits in homoploid hybrids between deciduous and evergreen species of Vaccinium

We investigated the anatomical expression of leaf traits in hybrids between evergreen Vaccinium vitis‐idaea and deciduous V. myrtillus. We compared parents from four populations with their respective F1 hybrids and tested whether (i) transgression can be the source of novel anatomical traits in hybrids; (ii) expression of transgressive traits is more probable for traits with similar values in parents and intermediate for more distinct values, as predicted by theory; and (iii) independent origin of hybrids leads to identical trait expression profiles among populations. We found that anatomical leaf traits can be divided into four categories based on their similarity to parents: intermediate, parental‐like, transgressive and non‐significant. Contrary to the common view, parental‐like trait values were equally important in shaping the hybrid profile, as were intermediate traits. Transgression was revealed in 17/144 cases and concerned mainly cell and tissue sizes. As predicted by theory, we observed transgressive segregation more often when there was little phenotypic divergence, but intermediate values when parental traits were differentiated. It is likely that cell and tissue sizes are phylogenetically more conserved due to stabilising selection, whereas traits such as leaf thickness and volume fraction of the intercellular spaces, showing a consistent intermediate pattern across populations, are more susceptible to directional selection. Hybrid populations showed little similarity in expression profile, with only three traits identically expressed across all populations. Thus local adaptation of parental species and specific genetic background may be of importance.     

Variability and adaptability of Miscanthus species evaluated for energy crop domestication

A growing body of evidence indicates that second‐generation energy crops can play an important role in the development of renewable energy and the mitigation of climate change. However, dedicated energy crops have yet to be domesticated in order to fully realize their productive potential under unfavorable soil and climatic conditions. To explore the possibility of domesticating Miscanthus crops in northern China where marginal and degraded land is abundant, we conducted common garden experiments at multiple locations to evaluate variation and adaptation of three Miscanthus species that are likely to serve as the wild progenitors of the energy crops. A total of 93 populations of Miscanthus sinensis, Miscanthus sacchariflorus, and Miscanthus lutarioriparius were collected across their natural distributional ranges in China and grown in three locations that represent temperate grassland with cold winter, the semiarid Loess Plateau, and relatively warm and wet central China. Evaluated with growth traits such as plant height, tiller number, tiller diameter, and flowering time, the Miscanthus species showed high levels of genetic variation within and between species. There were significant site × population interactions for almost all traits of M. sacchariflorus and M. sinensis, but not M. lutarioriparius. The northern populations of M. sacchariflorus had the highest establishment rates at the most northern site owing to their strong cold tolerance. An endemic species in central China, M. lutarioriparius, produced not only the highest biomass of the three species but also higher biomass at the Loess Plateau than the southern site near its native habitats. These results demonstrated that the wild species harbored a high level of genetic variation underlying traits important for crop establishment and production at sites that are colder and drier than their native habitats. The natural variation and adaptive plasticity found in the Miscanthus species indicated that they could provide valuable resources for the development of second‐generation energy crops.     

Hybridization and introgression among three Aconitum (Ranunculaceae) species of different ploidy levels in the Tatra Mountains (Western Carpathians)

Hybridization among species of A conitum effects their morphology and ecology. In this study the hybridization between the diploid 2n(2x) = 16 ( A . lasiocarpum and A . variegatum) and tetraploid 2n(4x) = 32 ( A . firmum) species was documented in the Tatra Mountains (Western Carpathians) in a small, local population in a semi‐natural site. The hybrid species were: homoploid A . × pawlowskii ( A . lasiocarpum × A . variegatum), and triploid A . × berdaui ( A . firmum × A . variegatum, 2n(3x) = 24). Chloroplast DNA (cpDNA) alleles formed two distinct haplotypes, one typical for the tetraploid and another for diploid lines, shared between the tetraploid, triploid and diploid groups, indicating introgressive hybridization. The presumed gene flow was from the tetraploid to diploid species via the triploid bridge. The only two specimens of A . × pawlowskii that harbored tetraploid ( A . firmum) type cpDNA possessed bracteoles of A . firmum‐type. The remaining introgressed (cpDNA and Inter Simple Sequence Repeats (ISSR)) specimens ( A . variegatum) were morphologically pure, implying cryptic introgression. ISSR loci shared between the tetraploid A . firmum and diploid A . variegatum support the hypothesis of an ancient allopolyploid origin of A . firmum and the diploid species of A . variegatum‐type as one of its parent.     

Cytogeography and chromosomal variation of the endemic East Asian herb Lycoris radiata

Information on the spatial distribution of cytotypes and karyotype variation in plants is critical for studies of the origin and evolution of polyploid complexes. Here, the spatial distribution of cytological races and intraspecific variation in the karyotype of Lycoris radiata, an endemic species to East Asia, is investigated. Conventional karyotype analysis methods were used to determine ploidy level and karyotypical characteristics in 2,420 individuals from 114 populations of L. radiata nearly covering the whole distribution areas in China. Of 114 populations studied, 52 (45.61%), 58 (50.88%), and 4 (3.51%) are diploid, triploid, and mixoploid populations, respectively, with 1,224, 1,195, and 1 individuals being diploid, triploid, and tetraploid, respectively. The triploid possesses a much wider distribution range than the diploid, with the former almost occupying the entire range of this complex species in East Asia and the latter distributing in the middle and east regions of China. Triploids tend to occur at high altitudes, and the relationship between the ploidy and altitude is significantly positive but low (r² = 0.103, p < 0.01). About 98.6% of examined bulbs have a common karyotype consisting of 22 or 33 acrocentric (A) chromosomes. Some aberrant chromosomes which should be generated from A‐type chromosome have been found including metacentrics (m), small metacentrics (m′), and B‐type chromosome. The results can provide a fundamental cytogeographic data for further studies on the evolutionary origins and adaptive divergences of polyploids, especially the triploid, within L. radiata using molecular and/or ecological methods in the future.     

Genomic variation in cline shape across a hybrid zone

Hybrid zones are unique biological interfaces that reveal both population level and species level evolutionary processes. A genome‐scale approach to assess gene flow across hybrid zones is vital, and now possible. In Mexican towhees (genus Pipilo), several morphological hybrid gradients exist. We completed a genome survey across one such gradient (9 populations, 140 birds) using mitochondrial DNA, 28 isozyme, and 377 AFLP markers. To assess variation in introgression among loci, cline parameters (i.e., width, center) for the 61 clinally varying loci were estimated and compiled into genomic distributions for tests against three empirical models spanning the range of observed cline shape. No single model accounts for observed variation in cline shape among loci. Numerous backcross individuals near the gradient center confirm a hybrid origin for these populations, contrary to a previous hypothesis based on social mimicry and character displacement. In addition, the observed variation does not bin into well‐defined categories of locus types (e.g., neutral vs. highly selected). Our multi‐locus analysis reveals cross‐genomic variation in selective constraints on gene flow and locus‐specific flexibility in the permeability of the interspecies membrane.     

Subspecies delineation amid phenotypic,geographic and genetic discordance in a songbird

Understanding the processes that drive divergence within and among species is a long‐standing goal in evolutionary biology. Traditional approaches to assessing differentiation rely on phenotypes to identify intra‐ and interspecific variation, but many species express subtle morphological gradients in which boundaries among forms are unclear. This intraspecific variation may be driven by differential adaptation to local conditions and may thereby reflect the evolutionary potential within a species. Here, we combine genetic and morphological data to evaluate intraspecific variation within the Nelson's (Ammodramus nelsoni) and salt marsh (Ammodramus caudacutus) sparrow complex, a group with populations that span considerable geographic distributions and a habitat gradient. We evaluated genetic structure among and within five putative subspecies of A. nelsoni and A. caudacutus using a reduced‐representation sequencing approach to generate a panel of 1929 SNPs among 69 individuals. Although we detected morphological differences among some groups, individuals sorted along a continuous phenotypic gradient. In contrast, the genetic data identified three distinct clusters corresponding to populations that inhabit coastal salt marsh, interior freshwater marsh and coastal brackish–water marsh habitats. These patterns support the current species‐level recognition but do not match the subspecies‐level taxonomy within each species—a finding which may have important conservation implications. We identified loci exhibiting patterns of elevated divergence among and within these species, indicating a role for local selective pressures in driving patterns of differentiation across the complex. We conclude that this evidence for adaptive variation among subspecies warrants the consideration of evolutionary potential and genetic novelty when identifying conservation units for this group.     

Intraspecific trait variation and trade‐offs within and across populations of a toxic dinoflagellate

Intraspecific trait diversity can promote the success of a species, as complementarity of functional traits within populations may enhance its competitive success and facilitates resilience to changing environmental conditions. Here, we experimentally determined the variation and relationships between traits in 15 strains of the toxic dinoflagellate Alexandrium ostenfeldii derived from two populations. Measured traits included growth rate, cell size, elemental composition, nitrogen uptake kinetics, toxin production and allelochemical potency. Our results demonstrate substantial variation in all analysed traits both within and across populations, particularly in nitrogen affinity, which was even comparable to interspecific variation across phytoplankton species. We found distinct trade‐offs between maximum nitrogen uptake rate and affinity, and between defensive and competitive traits. Furthermore, we identified differences in trait variation between the genetically similar populations. The observed high trait variation may facilitate development and resilience of harmful algal blooms under dynamic environmental conditions.