Spatial autocorrelation in two Iris pumila populations estimated on morphological data from natural clones and their samples grown in two different habitats

Morphological data from two Iris pumila populations (measured on native clones, on their replants into the same habitat, and on their transplants into alternative habitat) were combined with native clones spatial position and spatial autocorrelations (SA) were calculated. Naturally growing I. pumila clones revealed significant SA that were positive on small distances and negative on medium ones in both open Hillock and shaded Woodland populations. No significant SA were detected when calculated with original clone positions, but with morphometric data from replants into the experimental plot in the same habitat. Some significant SA were, however, detected when morphometric data from transplants to alternative habitat were used. Detected SA on I. pumila clones were primarily a consequence of spatial structuring of environmental factors but also, in a lesser degree, a result of genetic spatial arrangements (most probably due to patterns of gene flow). The text was submitted by the authors in English.     

Ontogeny of flower parts on naturally growing Iris pumila clones: implications for population differentiation and phenotypic plasticity studies

Previous studies revealed significant phenotypic plasticity, genetic variability and population differentiation of flower morphometric traits on dwarf bearded iris Iris pumila. Also, study of I. pumila flowering phenology revealed significant impact of habitat type as well as population differentiation for flowering time. Since the flowering time can influence other flower traits, we performed this analysis of flower morphometric traits in three time points during the flower bud ontogenic development in two habitat types (open vs. shaded). Analysis revealed that for most of the traits greater trait values were recorded for open habitat but only on latter time points. For most of the analyzed traits direction of differences in bud stage was the opposite to the direction of differences in mature flower stage detected in previous studies. However, length of the stem, a trait that showed the greatest variability between habitats and populations and therefore greatest genetic differentiation and phenotypic plasticity, was significantly greater in the samples from the late flowering shaded habitat in all time samples, indicating that in case of this trait different mechanisms were involved. Those findings have implications for design of the future studies on I. pumila.     

Variation in survival of Iris pumila L. juvenile plants grown in two nutrient regimes

I analyzed variability of early life survival on 881 Iris pumila juvenile plants in their prereproductive period. Seedlings were produced by applying half-sib mating scheme on clones originating from two Deliblato Sand populations that were accommodated to common garden conditions. Nineteen clones served as pollen donors (sires), and the other 90 clones served as pollen recipients (dames). Seedlings were grown in two nutrient levels (full strength and 1/10 of full-strength Hoagland solution). While I failed to detect significant mean nutrient level effect as well as significant between population differentiation for prereproductive survival, a statistically significant additive genetic variability for this ultimate prereproductive measure of fitness was detected.     

Phenotypic plasticity and between population differentiation in <Emphasis Type="Italic">Iris pumila</Emphasis> transplants between native open and anthropogenic shade habitats

Response to environmental heterogeneity caused by human activity was analyzed on Iris pumila reciprocal transplants between native steppe and anthropogenic (planted pine forest) habitats that were monitored during several growing seasons in a protected area of Deliblato Sand. Morphometric traits exhibited significant plastic responses to the environmental variability between native and anthropogenic habitats that differed in light quantity and quality, as well as in some other ecological indices. Significant differentiation between populations occupying those habitats was also detected. Plastic responses and population differences were substantially related to the variation in general size and had the same direction, plastic responses being larger in magnitude. Estimates of reproductive and vegetative performance of reciprocal transplants detected home site advantage in the native open but not in the secondary shade habitat created under anthropogenic influence.     

Flowering phenology in natural populations of Iris pumila

In the present study, I examined flowering phenology of dwarf bearded ins, Iris pumila , using naturally occurring clones (three sites, four microsites) and clones transplanted between two sites representing two habitats Naturally occurring clones in contrasting habitats and microsites differed significantly in phenology, with exposition of the site or microsite being the most important Genotypes from more exposed habitats flowered earlier in both habitats and these differences were statistically significant Patterns of between-habitat and between-population differences were stable over years even though years differed markedly in the flowering onset Within-population genetic variability for flowering phenology was also found to be significant Although clones with dark-colored flowers tended to flower earlier in all habitats and microsites I failed to detect statistically significant differences in flowering time among color morphs Flowering phenology in I pumila is highly susceptible to environmental variability, but this factor has not prevented population divergence in flowering time Between-habitat differences in flowering time turned out to be a result of both phenotypic plasticity of individuals and genetic differentiation of populations     

Dung beetle movements at two spatial scales

To understand the dynamics of spatially structured populations, we need to know the level of movements at different spatial scales. This paper reports on Aphodius dung beetle movements at two scales: movements between dung pats within pastures, and movements between pastures. First, I test an assumption common to many recent models of spatially structured populations – that the probability of an individual moving between habitat patches decreases exponentially with distance. For dung beetles, I find sufficient evidence to reject this assumption. The distribution of dispersal distances was clearly leptokurtic, with more individuals moving short and long distances than expected on the basis of an exponential function. In contrast, the data were well described by a power function. I conclude that dung beetle movements include an element of non-randomness not captured by the simplistic exponential model. The power function offers a promising alternative, but the actual mechanisms behind the pattern need to be clarified. Second, I compare several species of Aphodius to each other. Although these species occur in the same network of habitat patches, their movement patterns are different enough to result in a mixture of different spatial population structures. Movements between pastures were more frequent the larger the species, the more specific its occurrence in relation to pat age, and the more specialized it is on cow dung and open pasture habitats. Within pastures, all species form "patchy" populations, with much movement among individual pats.     

Genet age in marginal populations of two clonal Carex species in the Siberian Arctic

During a Swedish-Russian expedition to northern Siberia 1994, we sampled two marginal populations of two Carex species at two high arctic sites ( C, stans Drej. on Faddeyevsky Island and C. ensifolia V. Krecz ssp. arctisibirica Jurtz. at north-eastern Taymyr Peninsula), both north of previously documented localities in that areas for the two species. These populations were composed of a few distinct patches of ramet colonies, some of them shaped like fairy rings with dead centres. We measured the size of all colonies and collected samples for detailed morphometric analyses of rhizome growth. By using RAPD (random amplified polymorphic DNA) analysis we established that the largest colony at each site consisted of a single genet, based on 41 polymorphic bands amplified with three primers. Pooled samples from each of two additional colonies of C. stans on Faddeyevsky Island were analysed and showed that clones of the same species at the same site were relatively dissimilar (Dice's similarity index 0.26 0,43), We then assumed that each ramet colony represented a single genet. Based on the morphometric data, we developed a deterministic growth model that simulates the clonal growth of these species and enabled estimates of the time since establishment of the genets. The estimated age of the five C. Stans clones varied from 17 to 154 yr and the age of the two C. ensifolia ssp, arctisibirica clones was well over 3000 yr.     

Clonal and fine-scale genetic structure in populations of a restricted Korean endemic, Hosta jonesii (Liliaceae) and the implications for conservation

BACKGROUND AND AIMS: In plant populations the magnitude of spatial genetic structure of apparent individuals (including clonal ramets) can be different from that of sexual individuals (genets). Thus, distinguishing the effects of clonal versus sexual individuals in population genetic analyses could provide important insights for evolutionary biology and conservation. To investigate the effects of clonal spread on the fine-scale spatial genetic structure within plant populations, Hosta jonesii (Liliaceae), an endemic species to Korea, was chosen as a study species. METHODS: Using allozymes as genetic markers, spatial autocorrelation analysis of ramets and of genets was conducted to quantify the spatial scale of clonal spread and genotype distribution in two populations of H. jonesii. KEY RESULTS: Join-count statistics revealed that most clones are significantly aggregated at < 3-m interplant distance. Spatial autocorrelation analysis of all individuals resulted in significantly higher Moran's I values at 0-3-m interplant distance than analyses of population samples in which clones were excluded. However, significant fine-scale genetic structure was still observed when clones were excluded. CONCLUSIONS: These results suggest that clones enhance the magnitude of spatial autocorrelation due to localized clonal spread. The significant fine-scale genetic structure detected in samples excluding clones is consistent with the biological and ecological traits exhibited by H. jonesii including bee pollination and limited seed dispersal. For conservation purposes, genetic diversity would be maximized in local populations of H. jonesii by collecting or preserving individuals that are spaced at least 5 m apart.     

Ontogeny of flower parts on naturally growing <Emphasis Type="Italic">Iris pumila</Emphasis> clones: Implications for population differentiation and phenotypic plasticity studies

Previous studies revealed significant phenotypic plasticity, genetic variability and population differentiation of flower morphometric traits on dwarf bearded iris Iris pumila. Also, study of I. pumila flowering phenology revealed significant impact of habitat type as well as population differentiation for flowering time. Since the flowering time can influence other flower traits, we performed this analysis of flower morphometric traits in three time points during the flower bud ontogenic development in two habitat types (open vs. shaded). Analysis revealed that for most of the traits greater trait values were recorded for open habitat but only on latter time points. For most of the analyzed traits direction of differences in bud stage was the opposite to the direction of differences in mature flower stage detected in previous studies. However, length of the stem, a trait that showed the greatest variability between habitats and populations and therefore greatest genetic differentiation and phenotypic plasticity, was significantly greater in the samples from the late flowering shaded habitat in all time samples, indicating that in case of this trait different mechanisms were involved. Those findings have implications for design of the future studies on I. pumila.     

Further Studies on Differentiation of Photosynthetic Properties in Sun and Shade Ecotypes of Solidago virgaurea

Measurements of the fraction of the incident light absorbed by diverse Solidago leaves revealed that differences in light harvesting capacity cannot explain the differences in efficiency of utilization of weak light in photosynthesis that have previously been shown to exist between sun and shade ecotypes when these have been grown in strong light and between identical clones of shade ecotypes when grown at different light intensities. Photosynthesis measurements at low and normal oxygen concentrations, provided no evidence that a different degree of inhibition of photo-synthetic CO₂ uptake by atmospheric oxygen is responsible for the observed differences in photosynthetic efficiency, at low or high light intensities. These results support the conclusion that the markedly less efficient use of weak light by shaded habitat clones grown in strong as compared with weak light is caused primarily by damage to the photosystems, or to a site close to them. Measurements of Emerson enhancement and of light-induced absorbance changes provide some evidence that photoreaction II is more affected than I. Enzyme extracts prepared from clones native to an exposed habitat were found to contain considerably higher activities of carboxydismutase (ribulosc-l,5-diphos-phate carboxylase) than from clones native to a shaded habitat when the plants were previously grown at a moderately high light intensity. Exposed habitat clones apparently have a genetically determined, higher capacity to produce the carboxyla-tion enzyme than shaded habitat clones. The high degree of correlation found when the light-saturated rate of CO₂ uptake in vivo of a number of individual Solidago leaves is plotted against the carboxydismutase activities found in the extracts of these same leaves suggests that low carboxydismutase activity is one of the intrinsic properties responsible for the low capacity for light-saturated photosynthesis of clones from shaded habitats. It is concluded from this and other investigations that differentiation between plants from habitats with contrasting light intensities, whether unrelated species or ecotypos of the same species, probably involves the capacity of several component steps of the photosynthetic process.     

Assessment of variation in bacterial composition among microhabitats in a mangrove environment using DGGE fingerprints and barcoded pyrosequencing

Here, we use DGGE fingerprinting and barcoded pyrosequencing data, at six cut-off levels (85-100%), of all bacteria, Alphaproteobacteria and Betaproteobacteria to assess composition in the rhizosphere of nursery plants and nursery-raised transplants, native plants and bulk sediment in a mangrove habitat. When comparing compositional data based on DGGE fingerprinting and barcoded pyrosequencing at different cut-off levels, all revealed highly significant differences in composition among microhabitats. Procrustes superimposition revealed that ordination results using cut-off levels from 85-100% and DGGE fingerprint data were highly congruent with the standard 97% cut-off level. The various approaches revealed a primary gradient in composition from nursery to mangrove samples. The affinity between the nursery and transplants was greatest when using Betaproteobacteria followed by Alphaproteobacteria data. There was a distinct secondary gradient in composition from transplants to bulk sediment with native plants intermediate, which was most prevalent using all bacteria at intermediate cut-off levels (92-97%). Our results show that PCR-DGGE provides a robust and cost effective exploratory approach and is effective in distinguishing among a priori defined groups.     

Local population success in heterogeneous habitats: reciprocal transplant experiments completed on a desert spider

Reciprocal transplant experiments were completed to test for selection against the mixing of behavioural phenotypes in a desert spider. Most Agelenopsis aperta populations experience low prey abundances and competition for web‐sites that provide shelter from thermal extremes. These conditions favour aggressiveness towards both prey and conspecifics (an ‘arid‐land behavioural phenotype’). The spider also occupies narrow stretches of riparian habitat bordering spring‐fed streams and rivers. Here it is released from competition for prey and foraging sites, but is subject to predation by birds. A less aggressive/more fearful behaviour is selected for in these riparian habitats (a ‘riparian behavioural phenotype’). Previous work with this spider indicates that there is genetic differentiation between arid‐land and riparian populations. However, the degree to which genetic differentiation is achieved may be limited by gene flow. Reciprocal sets of enclosures were established in: (1) a dry evergreen woodland site (arid‐land phenotype) and (2) a neighbouring riparian site (riparian phenotype) in south‐eastern Arizona. Equal numbers of field collected, early instar A. aperta were introduced into native and transplant enclosures in each habitat. After 6 months of site‐imposed selection, survivorship was determined and growth estimates and behavioural trials completed on spiders remaining in the different enclosures. The same behavioural test was subsequently applied to lab‐reared offspring of the spiders surviving the respective selection regimes. Riparian transplants showed both poor survival and retarded growth in the dry woodland habitat when compared with both arid‐land and riparian natives. Arid‐land transplants that survived, however, grew equally well in riparian habitat as did dry woodland and riparian natives. Behavioural assays conducted on test subjects after selection and on their offspring reared in a controlled laboratory environment indicate that phenotypes that were inappropriate to the respective habitats were selected against in the transplant experiments. The frequency distribution of transplant spider behaviour on a continuum from fearful to aggressive was intermediate between that exhibited by respective native riparian and dry woodland spiders. It is concluded that while arid‐land and riparian behavioural ecotypes do exist, directional gene flow of arid‐land phenotypes into riparian habitat limits population subdivision.     

Incorporation of DPP6a and DPP6K variants in ternary Kv4 channel complex reconstitutes properties of A-type K current in rat cerebellar granule cells

Dipeptidyl peptidase-like protein 6 (DPP6) proteins co-assemble with Kv4 channel α-subunits and Kv channel-interacting proteins (KChIPs) to form channel protein complexes underlying neuronal somatodendritic A-type potassium current (I(SA)). DPP6 proteins are expressed as N-terminal variants (DPP6a, DPP6K, DPP6S, DPP6L) that result from alternative mRNA initiation and exhibit overlapping expression patterns. Here, we study the role DPP6 variants play in shaping the functional properties of I(SA) found in cerebellar granule (CG) cells using quantitative RT-PCR and voltage-clamp recordings of whole-cell currents from reconstituted channel complexes and native I(SA) channels. Differential expression of DPP6 variants was detected in rat CG cells, with DPP6K (41 ± 3%)>DPP6a (33 ± 3%)>DPP6S (18 ± 2%)>DPP6L (8 ± 3%). To better understand how DPP6 variants shape native neuronal I(SA), we focused on studying interactions between the two dominant variants, DPP6K and DPP6a. Although previous studies did not identify unique functional effects of DPP6K, we find that the unique N-terminus of DPP6K modulates the effects of KChIP proteins, slowing recovery and producing a negative shift in the steady-state inactivation curve. By contrast, DPP6a uses its distinct N-terminus to directly confer rapid N-type inactivation independently of KChIP3a. When DPP6a and DPP6K are co-expressed in ratios similar to those found in CG cells, their distinct effects compete in modulating channel function. The more rapid inactivation from DPP6a dominates during strong depolarization; however, DPP6K produces a negative shift in the steady-state inactivation curve and introduces a slow phase of recovery from inactivation. A direct comparison to the native CG cell I(SA) shows that these mixed effects are present in the native channels. Our results support the hypothesis that the precise expression and co-assembly of different auxiliary subunit variants are important factors in shaping the I(SA) functional properties in specific neuronal populations.     

不同年龄巴山木竹种群克隆结构

很多竹类植物是典型的克隆植物,也是大熊猫的食物。研究典型竹子种群克隆结构的形成和发展对竹林的生产和抚育具有理论和实践意义,可为预测该竹林群落的演替趋势和大熊猫保护提供科学依据。利用SSR标记研究不同年龄A(7龄)、B(30龄)和C(60龄)巴山木竹种群的克隆结构和多样性,探讨小尺度范围内不同年龄巴山木竹种群的克隆结构及斑块的建立和发展。8对SSR引物共扩增出了118个位点,3个种群样地的256个样本共检测到了49个克隆(基因型),A、B和C种群分别检测出31、10个和8个克隆。随着种群年龄的增长,巴山木竹克隆面积增加,克隆数量减少;A和B样地各克隆分布格局为团块状,而C样地克隆既有团块状又有离散状。这一结果显示出在幼苗定居的初期,基株可能以短距离的克隆延伸为主从而呈现出团块状;而随着年龄的增长,克隆面积不断扩大,当复轴混生型的巴山木竹克隆受到强大的压迫时,基株可能会进行较多的单轴和长距离克隆延伸,呈现出离散状。Mantel检测和空间自相关分析都支持3个样地在小尺度范围内存在明显的克隆空间遗传结构。3个样地在10 m等级下显著的正相关空间遗传结构距离为3.1、28、48 m,X-轴截距为9.051、30.698和50.536,空间自相关系数的范围分别为0.1—0.167、0.008—0.703和0.006—0.735。由此可推断,随着年龄的增长,巴山木竹克隆斑块的规模在不断地扩大,同一克隆的分株数量增加,在均匀取样情况下,正相关空间遗传结构距离范围内取到具有相同基因型的可能性越大。A、B和C 3样地的基因型比率(G/N)为1、0.14和0.055,Simpson多样性指数(D)分别为1、0.876和0.744。这说明巴山木竹幼苗期基因型比例远远高于成年的竹林,随着年龄的增长巴山木竹克隆多样性虽有所降低,但由于有性繁殖的作用仍然保持了较高的多样性。聚类和主坐标分析均表明总体上各样地的克隆被聚为一类,但不同样地少数克隆的基因型有重叠和聚集,可推断出不同巴山木竹种群之间可能存在着基因流动和近似的克隆起源。     

Genetic variation and population differentiation in the lichen-forming ascomycete Xanthoria parietina on the island Storfosna, central Norway

Genetic diversity and fine-scale population structure in the lichen-forming ascomycete Xanthoria parietina was investigated using sequence variation in part of the intergenic spacer (IGS) and the complete internal transcribed spacer (ITS) regions of the nuclear ribosomal DNA. Sampling included 213 and 225 individuals, respectively, from seven populations in two different habitats, bark and rock, on the island Storfosna off the central west coast of Norway. Both markers revealed significant variation and a total of 10 IGS and 16 ITS haplotypes were found. There were no signs of significant positive spatial autocorrelation at any spatial size class down to 10% of transect length, nor did we find significant deviations from neutrality or signs of historical population expansion. Analysis of molecular variance (amova) indicated that most of the genetic variance observed was within populations, but when populations were grouped according to habitat, more than a quarter of the variance was explained among groups. Pairwise comparisons of populations (F(ST), exact tests of population differentiation) revealed significant differentiation between populations in different habitats (on bark or rock), but not between populations in the same habitat. Haplotype networks show that internal and presumably old haplotypes are shared between habitats, whereas terminal haplotypes tend to be unique to a habitat, mostly bark. We interpret the observed pattern to mean that there is no evidence of restricted gene flow between populations in the same habitat at the present spatial scale (interpopulation distances one or a few kilometres). On the other hand, differentiation between habitats is considerable, which we attribute to restricted gene flow between habitats (habitat isolation). Evidence suggests that the observed differentiation did not evolve locally. Estimates of divergence time between populations in the respective habitats indicate that an ancestral population started to diverge at least 34,000 years ago but probably much further back in time.     

Genetic structure of a native cyprinid in a reservoir‐altered stream network

1. Reservoirs modify riverine ecosystems worldwide, and often with deleterious impacts on native biota. The immediate effects of reservoirs on native fish species below dams and in impounded reaches have received considerable attention, but it is unclear how reservoirs may affect fish species at larger spatial and temporal scales. Documented declines of stream fish populations in direct tributaries of reservoirs suggest reservoir pools may reduce gene flow among historically connected populations. 2. Because of increased predator densities in reservoirs and the extent of habitat alteration in impounded reaches, I predicted reservoir habitats would reduce gene flow among small‐bodied fish populations separated by reservoir habitat. I used microsatellite markers to assess the spatial genetic structure of populations of the red shiner (Cyprinella lutrensis), in a reservoir‐fragmented stream network (Lake Texoma, U.S.A.). I also tested the prediction that populations in two direct tributaries that have experienced population declines would have low genetic diversity. Individuals were collected from six sites upstream of the reservoir, three sites in the reservoir and three sites in direct tributaries of the reservoir during 2008 and 2009. 3. Results indicate that most populations were isolated by distance with little divergence among populations. In one direct tributary population, however, there was substantial genetic divergence, and genetic diversity was significantly lower than in other populations. Gene flow also seemed to be lower in reservoir habitats than in intact stream habitats, suggesting reservoir habitats may be reducing gene flow among the reservoir‐separated populations. These results indicate that reservoirs may reduce gene flow among reservoir‐fragmented stream fish populations, altering the evolutionary trajectories of fragmented populations.     

Limits to local adaptation in six populations of the annual plant Diodia teres

* Local adaptation is common, but tests for adaptive differentiation frequently compare populations from strongly divergent environments, making it unlikely that any influence of stochastic processes such as drift or mutation on local adaptation will be detected. Here, the hypothesis that local adaptation is more likely to develop when the native environments of populations are more distinct than when they are similar was tested. * A reciprocal transplant experiment including two populations from each of three habitats was conducted to determine the pattern of local adaptation. In addition to testing for local adaptation at the population level, the hypothesis was tested that local adaptation is more common between populations from different habitats than between populations from the same habitat. * Local adaptation was not common, but more evidence was found of local adaptation between populations from different habitats than between populations from the same habitat. Two instances of foreign genotype fitness advantage confirm that stochastic processes such as drift can limit local adaptation. * These results are consistent with the hypothesis that stochastic processes can inhibit local adaptation but are more likely to be overwhelmed by natural selection when populations occur in divergent environments.     

Genetic and morphological divergence among Gravel Bank Grasshoppers, <Emphasis Type="BoldItalic">Chorthippus pullus</Emphasis> (Acrididae), from contrasting environments

Gravel Bank Grasshopper (Chorthippus pullus) populations inhabit two contrasting environments, pebbly gravel banks with scarce vegetation cover in mountainous areas along the Alps and lowland grasslands dominated by Common Heather (Calluna vulgaris). Heath populations of C. pullus have been rediscovered only recently, and show a distribution scattered across Central Europe. The wings are reduced in this species; thus, it has low potential for long-distance dispersal. We used sequence data on a newly developed non-coding nuclear marker from three gravel-bank and four heath populations to test whether grasshoppers from the two environments represent distinct lineages. Gravel-bank populations were studied in southern Germany (Bavaria), heath populations in eastern Germany (Brandenburg and Saxony) and Ukraine. We compared those genetic data with an analysis of variation in a suite of morphometric traits. Finally, we combined genetic and morphometric data to reconstruct a plausible scenario for the ecological shift observed in C. pullus. Our newly developed marker did not sort populations from contrasting environments in two monophyletic lineages. Nevertheless, we found a general lack of gene flow between the gravel-bank and heath populations. There was pronounced variation among populations in morphometric traits. That variation was partially partitioned by habitat type, and populations from the same habitat tended to be more similar than those from different habitats. Our data suggest that heath populations originated through northward expansion from multiple southern European refugia, and that the gravel-bank populations represent one of these sources. Patterns of genetic and morphometric divergence suggest that gravel-bank and heath populations may be in the process of incipient speciation.     

Variation in early developmental stages in two populations of an intertidal crab, Neohelice (Chasmagnathus) granulata

Duration of embryonic development, egg size, larval size at hatching, and starvation tolerance of the first zoeal stage were studied in an intertidal crab from the southwestern Atlantic, Neohelice (formerly Chasmagnathus) granulata. These reproductive traits were quantified comparing (a) two populations living in ecologically contrasting coastal habitats in Argentina, a brackish lagoon, Mar Chiquita, MC vs. an open marine habitat near San Antonio, Patagonia, SA, (b) beginning vs. end of the reproductive season, and (c) two temperatures during egg development (18 vs. 27°C). Eggs in an early stage of embryonic development were in both populations larger at the beginning than at the end of the season, and were consistently larger in the SA population. These size differences persisted through larval hatching, independent of the temperature during embryogenesis. At 18°C, eggs produced at the beginning of the season developed in both populations more rapidly than those from the end of the reproductive season, while the opposite trend was observed at 27°C. The stage duration of the zoea I was in both populations shorter at the beginning as compared to the end of the season. The nutritional flexibility of the zoea I stage was compared using as indices the point-of-reserve-saturation (PRS₅₀) and the point-of-no-return (PNR₅₀). The PRS₅₀ was consistently lower in larvae from SA than in those from MC. In the MC population, this index was lower at the beginning than at the end of the season, while no significant seasonal difference was observed in larvae from SA. The PNR₅₀ varied between temperatures of embryonic development and populations, showing also significant interactions between all three factors. The PRS₅₀ was on average lower, and the PNR₅₀ was higher, than values previously reported for N. granulata, suggesting a stronger nutritional flexibility in the larvae used in the present study. Our results indicate significant intraspecific variability among separate populations, seasonal variation, and carry-over effects of environmental conditions prevailing during the embryonic phase, all of which may affect the performance of the larval phase.