Abundance of active and inactive microcystin genotypes in populations of the toxic cyanobacterium Planktothrix spp

To investigate the abundance of active and inactive microcystin genotypes in populations of the filamentous cyanobacterium Planktothrix spp., individual filaments were grown as clonal strains in the laboratory and analysed for microcystin synthetase (mcy) genes and microcystin. Twenty-three green-pigmented strains of P. agardhii originating mostly from shallow water bodies fell into two groups, those possessing mcyA and those lacking mcyA. In contrast, all of the 49 strains that were assigned to the red-pigmented P. rubescens contained mcyA. One strain of P. agardhii and eight strains of P. rubescens contained the total microcystin synthetase gene cluster but were found inactive in microcystin synthesis. To investigate the natural abundance of inactive mcy genotypes in P. rubescens individual filaments sampled from Lake Irrsee and Lake Mondsee (Austria) were analysed directly for the presence of mcyA and microcystin by matrix-assisted laser desorption/ionisation time-of-flight mass spectrometry. All filaments assigned to P. rubescens contained mcyA. The proportion of inactive microcystin genotypes in populations with a low (Irrsee) or high density (Mondsee) of P. rubescens was 5% and 21%, each. The results of this study demonstrate that P. rubescens typically contain mcy genes whereas P. agardhii have a patchy distribution of mcy genes. In both species microcystin producers co-occur with non-microcystin producers due to the absence/inactivation of mcy genes.     

Strategies for the co-existence of zooplankton with the toxic cyanobacterium Planktothrix rubescens in Lake Zurich

Since the cyanobacterium Planktothrix rubescens, which dominates the phytoplankton community in Lake Zurich, is generally considered toxic to zooplankton, we addressed the question whether co-occurring zooplankton species have developed adaptive responses. Artificially shortened filaments (<30 m in length) of P.rubescens significantly reduced survival of Thamnocephalus platyurus (Crustacea, Branchiopoda, Anostraca) naturally occurring in temporary ponds. In contrast to Thamnocephalus, the survival of co-existing zooplankton was unaffected (Eudiaptomus gracilis (or enhanced (Daphnia hyalina and Cyclops abyssorum). High sensitivity to the microcystins of Planktothrix was coupled to strict food avoidance in Eudiaptomus, but not in Thamnocephalus. Daphnia and Cyclops exhibited higher physiological resistance to cyanobacterial toxins, and ingested Planktothrix. For the lake zooplankton species, the feeding rates on high-quality algae were not significantly reduced in the presence of Planktothrix. In order to separate the effects of mechanical interference (filament length) versus toxins, clearance rates on Planktothrix filaments were compared to clearance rates on filaments subjected to toxin extraction. The results show that microcystins are important feeding deterrents against grazing by Daphnia since feeding rates on Planktothrix increased significantly after an aqueous-methanolic extraction of the major part of microcystins. On the other hand, copepods persisted in food avoidance, but exhibited high clearance rates on Planktothrix after a more lipophilic extraction was applied. Both microcystins and a lipophilic, unidentified toxin may contribute to the avoidance behaviour of copepods. For both Daphnia and copepods, the grazing resistance of Planktothrix is mediated by chemical defences rather than by the large size and the rigidity of the filaments.      

Diversity of microcystin genotypes among populations of the filamentous cyanobacteria Planktothrix rubescens and Planktothrix agardhii

Microcystins (MCs) are toxic heptapeptides that are produced by filamentous cyanobacteria Planktothrix rubescens and Planktothrix agardhii via nonribosomal peptide synthesis. MCs share a common structure cyclo (-D-Alanine(1)-L-X(2)- D-erythro-beta-iso-aspartic acid(3)-L-Z(4)-Adda(5)-D-Glutamate(6)- N-methyl-dehydroalanine(7)) where X(2) and Z(2) are variable L-amino acids in positions 2, 4 of the molecule. Part of the mcyB gene (1,451 bp) that is involved in the activation of the X(2) amino acid during MC synthesis was sequenced in 49 strains containing different proportions of arginine, homotyrosine, and leucine in position 2 of the MC molecule. Twenty-five genotypes were found that consisted of eight genotype groups (A-H, comprising 2-11 strains) and 17 unique genotypes. P. rubescens and P. agardhii partly consisted of the same mcyB genotypes. The occurrence of numerous putative recombination events that affected all of the genotypes can explain the conflict between taxonomy and mcyB genotype distribution. Genotypes B (homotyrosine and leucine in X(2)) and C (arginine in X(2)) showed higher nonsynonymous/synonymous (d(N)/d(S)) substitution ratios implying a relaxation of selective constraints. In contrast, other genotypes (arginine, leucine, homotyrosine) showed lowest d(N)/d(S) ratios implying purifying selection. Restriction fragment length polymorphism (RFLP) revealed the unambiguous identification of mcyB genotypes, which are indicative of variable X(2) amino acids in eight populations of P. rubescens in the Alps (Austria, Germany, and Switzerland). The populations were found to differ significantly in the proportion of specific genotypes and the number of genotypes that occurred over several years. It is concluded that spatial isolation might favour the genetic divergence of microcystin synthesis in Planktothrix spp.     

Modular structure of genes encoding multifunctional peptide synthetases required for non-ribosomal peptide synthesis

Abstract Peptide synthetases are large multienzyme complexes that catalyze the non-ribosomal synthesis of a structurally diverse family of bioactive peptides. They possess a multidomain structure and employ the thiotemplate mechanism to activate, modify and link together by amide or ester bonds the constituent amino acids of the peptide product. The domains, which represent the functional building units of peptide synthetases, appear to act as independent enzymes whose specific linkage order forms the protein-template that defines the sequence of the incorporated amino acids. Two types of domains have been characterized in peptide synthetases of bacterial and fungal origin: type I comprises about 600 amino acids and contains at least two modules involved in substrate recognition, adenylation and thioester formation, whereas type II domains carry in addition an insertion of about 430 amino acids that may function as a N-methyltransferase module. The role of other genes associated with bacterial opérons encoding peptide synthetases is also discussed.     

Structural and expressional divergence of genes encoding O-methyltransferase in wheat

Enzymatic methylation, which is catalyzed by the large number of O-methyltransferases (OMTs), is one of the important reactions in the flow of primary and (or) secondary metabolism. In a previous study, the gene TaOMT1 was induced by Hessian fly infestation of a wheat-rye translocation line. In this study we considered other wheat OMT genes - TaOMT3, TaOMT4, and TaOMT5 - using a bioinformatics approach and examined the TaOMT genes for their genomic organization, tissue-specific expression, responses to abiotic stresses and hormones, and cis-elements. There appeared to be a homoeologous relationship between TaOMT4 (6DS) and TaOMT5 (6BS), whereas TaOMT1 and TaOMT3 were placed on chromosome arms 7BL and 5DL, respectively. Differences in the tissue-specific, constitutive, and stress-inducible expression patterns among the TaOMT genes were found in both healthy and stressed plants. A number of cis-elements, which are potentially correlated with the responses of the TaOMT genes, were detected in the analysis of the TaOMT promoter sequences. In addition, evolutionary perspectives of the TaOMT genes are discussed. The nucleotide sequences have been deposited in the GenBank database under accession Nos. AAP23942 (TaCOMT1), EF423610 (TaOMT5), EF423611 (TaOMT4), EF423612 (TaOMT3), EU831287 (5' upstream of TaOMT1), EU831288 (5' upstream of TaOMT3), EU831289 (5' upstream of TaOMT4), and EU831290 (5' upstream of TaOMT5).     

Spatial synchrony of wader populations in inland lakes of the Iberian Peninsula

Spatial synchronization refers to similarity in temporal variations between spatially separated populations. Three mechanisms have been associated with the spatial synchrony of populations: Moran effect, dispersal and trophic interactions. In this study, we explored the degree of spatial synchrony of three wader species populations (Pied Avocet, Black-winged Stilt and Kentish Plover) using monthly estimates of their abundance in inland lakes of the Iberian Peninsula. The effect of several types of wetland variables (structural, hydroperiod and landscape) on spatial synchronization was explored. Groups of lakes with significant synchronization were identified for all three species. The lakes with wastewater input presented longer hydroperiods than those that did not receive these effluents, and this factor was positively related to the spatial synchrony of the Pied Avocet and Kentish Plover populations. The distance between lakes (used as an indicator of the dispersal effect on synchronization) was significant only in Pied Avocet. No structural or landscape variables were related to spatial synchronization in any species. It was impossible to identify any variable related to the spatial synchronization of Black-winged Stilt abundance as a possible result of the high ecological plasticity of this species. Our data provides the first evidence for mechanisms that act on the spatial synchronizing of wader populations in temporary continental lakes in central Spain, and show that the hydroperiod of lakes acts as an important factor in the spatial synchronization of aquatic species and that its effect is mediated by the reception of urban wastewater.     

Genetic divergence among geographical populations of the migratory locust in China

The random amplified polymorphic DNA (RAPD) technique was used to examine genetic divergence and interrelations of 11 geographical populations of the migratory locust in China, and the role of spatial separation in the population differentiations. AMOVA analysis of genetic variations in all the populations indicated greater within- (79.55%) than among-population variability (20.45%), and that there were significant differentiations among the populations; 11 populations were divided into four regional groups, with significantly greater variability within (82.99%) than among the groups (17.01%), and there existed apparent regional differentiations. Paired comparisons showed significantly greater variability within-than between-groups, indicating significant differentiations between populations of different regional groups. Of all the pairwise comparisons, Hainan and Tibetan groups displayed the greatest differentiation, with the difference between the two groups being seven folds of that between populations within the groups; the least differentiations were exhibited between the groups of Hainan, Xinjiang, and Inner Mongolia, with the differences between groups being only half of the differences between populations within the groups. Mantel tests of the genetic and spatial distances showed that the two matrices were significantly correlated (p<0.01), indicating that the geographical isolation played an important role in the differentiations of the geographical populations of the migratory locusts. Cluster analysis divided all populations into four major groups: Xinjiang and Inner Mongolia group, the Great Plains of North China (the Yellow River and Huai River Plains) group, Hainan group, and Tibet group. Principal component analysis (PCA) supported the division of populations based on the cluster analysis. However, analysis of individuals clustered the locusts into five populations: Xinjiang and Inner Mongolia, Hami in Xinjiang, the Great Plains of North China, Hainan, and Tibet. The locust populations in eastern China displayed apparently continous and gradient variations; as such authors consider that there were no necessity and valid reasons for further division of subspecies. The subspecific status for the main geographical populations of the migratory locusts in China was discussed.     

Molecular characterization of the genes encoding acetohydroxy acid synthase in the cyanobacterium Spirulina platensis.

The enzyme acetohydroxy acid synthase (AHS), which catalyses the first common step in the biosynthesis of isoleucine, leucine and valine, has been demonstrated to be present in Spirulina platensis in two isoenzymic forms. The complete nucleotide sequences of the genes ilvX and ilvW encoding these two enzymes have been determined. Sequence analysis revealed the presence of two open reading frames, of 1836 and 1737 nucleotides for ilvX and ilvW, respectively. The predicted amino acid sequences of the two isoenzymes, compared with the Synechococcus PCC 7942 AHS enzyme and the large subunits of the Escherichia coli AHSI, II, III isoenzymes, revealed a notable degree of similarity. A small subunit has not been identified for either of the S. platensis AHS isoenzymes. Analysis by Northern blot hybridization demonstrated that the ilvX and ilvW genes are transcribed to give mRNA species of approximately 2.15 kb and 1.95 kb, respectively.     

Amino acid availability determines the ratio of microcystin variants in the cyanobacterium Planktothrix agardhii

Cyanobacteria are capable of producing multiple microcystin variants simultaneously. The mechanisms that determine the composition of microcystin variants in cyanobacteria are still debated. [Asp³]microcystin-RR contains arginine at the position where the more toxic [Asp³]microcystin-LR incorporates leucine. We cultured the filamentous cyanobacterium Planktothrix agardhii strain 126/3 with and without external addition of leucine and arginine. Addition of leucine to the growth medium resulted in a strong increase of the [Asp³]microcystin LR/RR ratio, while addition of arginine resulted in a decrease. This demonstrates that amino acid availability plays a role in the synthesis of different microcystin variants. Environmental changes affecting cell metabolism may cause differences in the intracellular availability of leucine and arginine, which can thus affect the production of microcystin variants. Because leucine contains one nitrogen atom while arginine contains four nitrogen atoms, we hypothesized that low nitrogen availability might shift the amino acid composition in favor of leucine, which might explain seasonal increases in the [Asp³]microcystin LR/RR ratio in natural populations. However, when a continuous culture of P. agardhii was shifted from nitrogen-saturated to a nitrogen-limited mineral medium, leucine and arginine concentrations decreased, but the leucine/arginine ratio did not change. Accordingly, while the total microcystin concentration of the cells decreased, we did not observe changes in the [Asp³]microcystin LR/RR ratio in response to nitrogen limitation.     

Localization of polyketide synthase encoding genes to the toxic dinoflagellate Karenia brevis

Karenia brevis is a toxic marine dinoflagellate endemic to the Gulf of Mexico. Blooms of this harmful alga cause fish kills, marine mammal mortalities and neurotoxic shellfish poisonings. These harmful effects are attributed to a suite of polyketide secondary metabolites known as the brevetoxins. The carbon framework of all polyketides is assembled by a polyketide synthase (PKS). Previously, PKS encoding genes were amplified from K. brevis culture and their similarity to a PKS gene from the closely related protist, Cryptosporidium parvum, suggested that these genes originate from the dinoflagellate. However, K. brevis has not been grown axenically. The associated bacteria might be the source of the toxins or the PKS genes. Herein we report the localization of PKS encoding genes by a combination of flow cytometry/PCR and fluorescence in situ hybridization (FISH). Two genes localized exclusively to K. brevis cells while a third localized to both K. brevis and associated bacteria. While these genes have not yet been linked to toxin production, the work describes the first definitive evidence of resident PKS genes in any dinoflagellate.     

Spatial and temporal divergence of expression in duplicated barley germin-like protein-encoding genes

Subfunctionalization is the process by which a pair of duplicated genes, or paralogs, experiences a reduction of individual expression patterns or function while still reproducing the complete expression pattern and function of the ancestral gene. Two germin-like protein (GLP)-encoding genes, GerB and GerF, are paralogs that belong to a small gene family in barley (Hordeum vulgare). Both genes share high nucleotide sequence similarity in coding and noncoding regions and encode identical apoplastic proteins. The use of RNA gel blots, coupled with single-stranded conformation polymorphism (SSCP) analysis of RT-PCR products, elucidated the developmental and tissue-specific expression patterns of each gene. Individual expression patterns provided evidence of both overlapping redundancy and early subfunctionalization. GerB is predominantly expressed in developing shoots, while GerF is predominantly expressed in seedling roots, developing spikes, and pericarp/testa. GerF promoter deletion studies located a region (-356/-97) responsible for high promoter activity and showed the ability of GerB and GerF upstream regions to drive gfp expression in coleoptiles, epicarps, and lemma/palea of developing spikes. The observed expression patterns are consistent with proposed roles in plant development and defense mechanisms for this gene family. These roles may explain why redundancy has been selectively maintained in this duplicate gene pair.     

Spatial and temporal variability in filament length of a toxic cyanobacterium (Anabaena affinis)

1. This study examines the distribution of Anabaena affinis filament lengths under natural conditions as a function of depth and season, and in the laboratory as a function of growth phase. Because Anabaena affinis is only toxic when consumed, both its filament length and position in the water column are important determinants of its potential impact on zooplankton populations. 2. Star Lake (Norwich, Vermont, U.S.A.), a natural, eutrophic pond, remained thermally stratified throughout the Anabaena bloom. Filament number and length differed significantly with both sampling date and water depth. Most filaments occurred at 0.5 m, particularly at the height of the bloom. Throughout the entire water column average filament length decreased from approximately 0.53 mm in May to 0.14 mm in July. The shortest filaments occurred at the 2.5 m depth. Filament length distributions (combined for all depths) for 29 May, 12 June and 3 July, corresponding to the beginning, middle and end of the bloom, respectively, differed significantly among the three dates. These patterns most likely reflect variable growth conditions, both during the season and in the water column. 3. In the laboratory, Anabaena filament length was affected by medium composition and growth phase. Filaments were significantly longer when grown in MBL than in ASM medium. Also, the average length of Anabaena filaments grown in MBL changed significantly as cultures aged; by day 13 filament length (2.01 ± 0.38 mm, mean ± SD) was twice that on day 0 (0.97 ± 0.71 mm). As cell concentration continued to increase, mean filament length gradually decreased.     

Characterization of the locus of genes encoding enzymes producing heptadepsipeptide micropeptin in the unicellular cyanobacterium Microcystis

The gene cluster involved in producing the cyclic heptadepsipeptide micropeptin was cloned from the genome of the unicellular cyanobacterium Microcystis aeruginosa K-139. Sequencing revealed four genes encoding non-ribosomal peptide synthetases (NRPSs) that are highly similar to the gene cluster involved in cyanopeptolins biosynthesis. According to predictions based on the non-ribosomal consensus code, the order of the mcnABCE NPRS modules was well consistent with that of the biosynthetic assembly of cyclic peptides. The biochemical analysis of a McnB(K-139) adenylation domain and the knock-out of mcnC in a micropeptin-producing strain, M. viridis S-70, revealed that the mcn gene clusters were responsible for the production of heptadepsipeptide micropeptins. A detailed comparison of nucleotide sequences also showed that the regions between the mcnC and mcnE genes of M. aeruginosa K-139 retained short stretches of DNA homologous to halogenase genes involved in the synthesis of halogenated cyclic peptides of the cyanopeptolin class including anabaenopeptilides. This suggests that the mcn clusters of M. aeruginosa K-139 have lost the halogenase genes during evolution. Finally, a comparative bioinformatics analysis of the congenial gene cluster for depsipetide biosynthesis suggested the diversification and propagation of the NRPS genes in cyanobacteria.     

Multiplicity and regulation of genes encoding peptide transporters in Saccharomyces cerevisiae

The model eukaryote Saccharomyces cerevisiae has two distinct peptide transport mechanisms, one for di-/tripeptides (the PTR system) and another for tetra-/pentapeptides (the OPT system). The PTR system consists of three genes, PTR1, PTR2 and PTR3. The transporter (Ptr2p), encoded by the gene PTR2, is a 12 transmembrane domain (TMD) integral membrane protein that translocates di-/tripeptides. Homologues to Ptr2p have been identified in virtually all organisms examined to date and comprise the PTR family of transport proteins. In S. cerevisiae, the expression of PTR2 is highly regulated at the cellular level by complex interactions of many genes, including PTR1, PTR3, CUP9 and SSY1. Oligopeptides, consisting of four to five amino acids, are transported by the 12-14 TMD integral membrane protein Opt1p. Unlike Ptr2p, distribution of this protein appears limited to fungi and plants, and there appears to be three paralogues in S. cerevisiae. This transporter has an affinity for enkephalin, an endogenous mammalian pentapeptide, as well as for glutathione. Although it is known that OPT1 is normally expressed only during sporulation, to date little is known about the genes and proteins involved in the regulation of OPT1 expression.     

Multiplicity and regulation of genes encoding peptide transporters in Saccharomyces cerevisiae

The model eukaryote Saccharomyces cerevisiae has two distinct peptide transport mechanisms, one for di-/tripeptides (the PTR system) and another for tetra-/pentapeptides (the OPT system). The PTR system consists of three genes, PTR1, PTR2 and PTR3. The transporter (Ptr2p), encoded by the gene PTR2, is a 12 transmembrane domain (TMD) integral membrane protein that translocates di-/tripeptides. Homologues to Ptr2p have been identified in virtually all organisms examined to date and comprise the PTR family of transport proteins. In S. cerevisiae, the expression of PTR2 is highly regulated at the cellular level by complex interactions of many genes, including PTR1, PTR3, CUP9 and SSY1. Oligopeptides, consisting of four to five amino acids, are transported by the 12 - 14 TMD integral membrane protein Opt1p. Unlike Ptr2p, distribution of this protein appears limited to fungi and plants, and there appears to be three paralogues in S. cerevisiae. This transporter has an affinity for enkephalin, an endogenous mammalian pentapeptide, as well as for glutathione. Although it is known that OPT1 is normally expressed only during sporulation, to date little is known about the genes and proteins involved in the regulation of OPT1 expression.     

Planktothrix populations in subalpine lakes: selection for strains with strong gas vesicles as a function of lake depth,morphometry and circulation

1. The genus Planktothrix (Cyanobacteria) usually produces concentrated populations of filaments in the summer metalimnion of thermally stratifying lakes. This has been associated with the action of gas vesicles, cellular structures providing positive buoyancy. At the end of the summer, filaments are carried by convective mixing deeper into the water column where some gas vesicles collapse as a result of high hydrostatic pressure. They then lose their buoyancy, sink and are lost from the euphotic zone. 2. The resistance of gas vesicles to hydrostatic pressures is critical for the survival of Planktothrix in deep lakes. However, comparative observations on populations from lakes of a range of depths and hydrodynamic regimes are still needed to examine the relationships between the adaptive trait (i.e. the ‘critical’ pressure at which each gas vesicle collapses) with the environmental factor (i.e. the maximum hydrostatic pressure). 3. To explore the adaptation of Planktothrix populations to the depth of winter circulation in different systems, we collected 276 strains of P. cf. rubescens from eight lakes (z_max = 24–410 m) in Northern Italy during summer 2009 and we analysed the multicopy gene gvpC coding for a protein that crucially influences the critical pressure. 4. The strains analysed clustered into two main groups having gas vesicles with a mean critical pressure of 1.1 and 0.9 MPa, respectively. The proportion of the stronger strains was generally positively related to lake depth, although the overall pattern was complicated by individual lake morphology and hydrology. The relative frequency of stronger filaments was (i) greatest in deep basins with concave slopes and (ii) least in one deep, but permanently stratified lake. 5. The simultaneous presence of ‘weaker’ and ‘stronger’ filaments could allow for a rapid adaptive response to changes in hydrostatic pressures, related to changes in the amplitude of vertical circulation characterising deep lakes.     

Spatial and temporal patterns of light attenuation among lakes of the Mackenzie Delta

SUMMARY 1. The seasonal dynamics of light attenuation, and the relative roles of total suspended solids (TSS), dissolved organic carbon (DOC) and chlorophyll as light attenuators among two sets of lakes in the Mackenzie Delta, were assessed during the open‐water periods of 1998 and 1999. 2. The first set consisted of 40 spatially discrete lakes where the frequency of flooding with river water was controlled by sill height (‘sill‐set lakes’). The second set consisted of a chain of six lakes connected to a main river channel (frequently flooded, all with same frequency), but where riverine influence was controlled by the distance from the channel connection point (‘chain‐set lakes’). 3. As the flooding frequency of lakes decreased (sill‐set), and as the distance from the channel connection point increased (chain‐set), lake water became increasingly transparent and the stability (decreasing temporal variability) of underwater light increased. 4. The effect of flooding on transparency was greater in years with a high minimum summer water level. However, the effect of river flooding on lake water transparency was damped more by an increase in the frequency and duration of flooding than by an increase in distance from the channel connection point. 5. The index of scattering was linearly related to TSS over the common range of concentrations in both sets of lakes. The specific attenuation coefficient for TSS (and scattering) increased substantially from the most turbid to the most transparent waters. 6. During the summer, DOC provided an approximate index of water colour in the sill‐set lakes but not in the chain‐set lakes, where the gradient of DOC ran counter to the gradient of water colour. The specific attenuation coefficient for water colour was roughly constant among both sets of lakes. 7. Calculations of partial attenuation show that, during the spring flood peak, TSS is the dominant attenuator among most lakes, other than those with high sills or positioned far from channel connection points. During the lengthy summer period of open water, however, water colour appeared to be the most important light attenuator among almost all of the lakes in the central delta, with chlorophyll a of only minor importance. 8. Lakes of the Mackenzie Delta may be quite sensitive to changes in climate and ultraviolet‐b (UV‐b) radiation in the circumpolar arctic because of the role of DOC as an attenuator of photosynthetically active radiation and UV‐b irradiance and as an energy source for microbial foodwebs in this system.     

Local selection modifies phenotypic divergence among Rana temporaria populations in the presence of gene flow

In ectotherms, variation in life history traits among populations is common and suggests local adaptation. However, geographic variation itself is not a proof for local adaptation, as genetic drift and gene flow may also shape patterns of quantitative variation. We studied local and regional variation in means and phenotypic plasticity of larval life history traits in the common frog Rana temporaria using six populations from central Sweden, breeding in either open‐canopy or partially closed‐canopy ponds. To separate local adaptation from genetic drift, we compared differentiation in quantitative genetic traits (Q_ST) obtained from a common garden experiment with differentiation in presumably neutral microsatellite markers (F_ST). We found that R. temporaria populations differ in means and plasticities of life history traits in different temperatures at local, and in F_ST at regional scale. Comparisons of differentiation in quantitative traits and in molecular markers suggested that natural selection was responsible for the divergence in growth and development rates as well as in temperature‐induced plasticity, indicating local adaptation. However, at low temperature, the role of genetic drift could not be separated from selection. Phenotypes were correlated with forest canopy closure, but not with geographical or genetic distance. These results indicate that local adaptation can evolve in the presence of ongoing gene flow among the populations, and that natural selection is strong in this system.