Population genetics of eight species of Trapezia (Brachyura: Xanthidae), symbionts of corals

Starch-gel electrophoresis was used to study gene-enzyme variation in thirteen populations of eight species of the genus Trapezia from Hawaii, Panamá, and Enewetak Atoll (Marshall Islands). Between 20 and 30 (mean = 27.8) gene-enzyme systems were resolved in each population, with 20 systems in common among all populations. The distribution of electrophenotypes was in agreement with Hardy-Weinberg-Castle expectations, except for T. digitalis, which consistently showed heterozygote deficiencies. Diagnostic loci among color forms support the hypothesis that color forms are distinct species. Low values of genetic distance among species suggest a recent radiation, perhaps during the Pleistocene. Genetic distance between the Hawaiian and Panamanian populations of T. ferruginea did not significantly differ from zero, indicating that the Eastern Pacific population of T. ferruginea has recently immigrated from the central Pacific, and/or that there is gene flow between the two areas. There were diagnostic loci between T. ferruginea and T. formosa from Enewetak and the populations of these species from Hawaii (T. ferruginea only) and Panamá (both species). Therefore, these geographic populations may represent separate species. The level and pattern of genetic variability in Trapezia spp. are in agreement with those observed in most other organisms.     

Lack of genetically-subdivided population structure in Bullia digitalis, a southern African marine gastropod with lecithotrophic development

Larval shell morphology in fossil and present-day gastropods is often used to infer modes of larval development and levels of dispersal. Dispersal ability influences not only genetic population structure, but also is thought to influence a species' geographical range and evolutionary duration. We tested these predictions in Bullia digitalis, a sandy-beach whelk, by examining genetic variability at 33 protein-coding loci in nine samples (N=739) taken in 1984 to 1985 at localities extending over about three-quarters of the geographical range of this species in southern Africa. Females of this species deposit eggs into benthic or brooded capsules in which larvae develop through the trochophore and veliger stages to emerge as crawling juveniles. Scanning electron micrographs confirmed a protoconch morphology typical for gastropods with lecithotrophic larval development. Contrary to expectations, subpopulations of B.␣digitalis had high levels of variability (H=0.102) and lacked a genetically-fragmented structure (=0.013). The lack of a genetically-subdivided population structure would not have been correctly inferred, if this species were known only from well-preserved fossil shells. Indirect estimates of migration between populations based on and the island model of migration, which assumes drift-mutation equilibrium, ranged between 19 and 23 individuals per generation. Either an undescribed mechanism of dispersal facilitates gene flow between populations, or the geographical range of this species has recently expanded to produce the appearance of high levels of gene flow. Gene-frequency distributions showed that relative to four other species of Bullia, populations of B. digitalis were in mutation-drift disequilibrium, with a significant excess of low-frequency alleles that is consistent with a recent rapid expansion from a small population. Also contrary to expectations, this species has a large geographical range (2 400 km) and an apparently long evolutionary history extending 5 to 20 million years, as estimated from an allozyme phylogeny with four other species of Bullia. Received: 15 January 1997 / Accepted: 28 January 1997     

Gpi genotypic effect on quantitative traits in the northern bay scallop,Argopecten irradians irradians

The relationships between glucose-6-phosphate isomerase (Gpi) genotype and quantitative variation related to reproduction and growth were explored over the lifespan of a single cohort of northern bay scallops, Argopecten irradians irradians (Lamarck), from September 1986 to April 1988 in the Niantic River estuary, Connecticut, USA. Analyses revealed that Gpi genotype explained an increasingly significant proportion of variation in scallop size (up to 15%) as the cohort aged. The genotype-specific effects were consistent across sampling dates and among measured traits. There was no evidence for heterosis with respect to size at this locus; however, rare genotypes contributed substantially to the relationship and showed some tendency to fall on the extremes of the phenotypic distribution. The strength and consistency of the genotypic effect on scallop size suggest that genetic variation for Gpi, or some locus in linkage disequilibrium with Gpi, may translate into biochemical and/or physiological variation and affect fitness in this species.     

Foster,but not genetic,father plumage coloration has a temperature-dependent effect on offspring quality

Secondary sexual characters have most likely evolved through sexual selection because such traits indicate the genetic or phenotypic quality of the bearer. While genetic variation in such fitness-related traits should be depleted by directional selection, there are many cases in which variation is higher than expected. One hypothesis explaining this variation is that different phenotypes within a population are adapted to different environmental conditions. We tested this hypothesis by comparing the offspring quality of male pied flycatchers, Ficedula hypoleuca, with different degrees of melanin-based dorsal plumage coloration under different environmental conditions. To create different environmental growth conditions and to be able to separate offspring genotype effects from paternal effects on offspring body mass, we used a partial cross-foster design where the brood size was reduced or enlarged by one chick. We also examined the interactive effects of temperature and male phenotype because previous correlative studies suggest such temperature-dependent effects in this species. We show that, while manipulated brood size did not interact with male phenotype to affect offspring quality, temperature during the nestling period influenced the offspring quality of dark and brown foster (but not genetic) fathers. When the temperature was relatively low during the nestling period, foster offspring of black males were lighter than those raised by brown males; the opposite was true if temperature was relatively high. These results add a new aspect to our understanding of how variation in the degree of melanin-based coloration is maintained in wild populations and how phenotypic variation may be maintained in general.     

Biochemical genetic variation and evidence of speciation in Chthamalus barnacles of the tropical eastern Pacific Ocean

Gene-enzyme variation was studied electrophoretically within and between barnacle populations of the genus Chthamalus from 8 intertidal stations from central California to the Pacific coast of Panamá. Horizontal starch-gel electrophoresis separated and resolved 18 enzymes from 512 individual barnacles. A maximum of 25 gene-enzyme systems was interpretable from the resulting zymograms. Electrophoretic phenotypes and patterns of phenotypic variation generally conformed with those observed in other organisms. The amount of genetic variation within barnacle populations varied; average heterozygosity, for example, ranged from a low of 3.4% in a Mexican population of the C. fissus group to a high of 10.4% in C. anisopoma from the Gulf of California. Observed and expected average heterozygosities agreed in all population samples, indicating that these species are outbreeding. In contrast to the prediction stemming from a hypothesis that trophic stability regulates the amount of genetic variation in marine species, average heterozygosity tends to be positively correlated with latitude. Data from these and other barnacle species may support a hybrid environmental heterogeneity-trophic diversity model recently proposed to explain genetic variation in decapod crustaceans. Juxtaposition of individuals from different localities revealed numerous genetic differences among populations of the C. fissus group. At least three partially sympatric sibling species are separated by genetic distances as large as those observed between the C. fissus group and the distinct species C. dalli and C. anisopoma. A cladistic analysis places C. anisopoma close to the Mexican and Panamanian sibling species, with C. fissus from San Diego and C. dalli successively farther removed.     

Biochemical genetic variation at a leucine aminopeptidase (LAP) locus in blue (Mytilus galloprovincialis) and greenshell (Perna canaliculus) mussel populations along a salinity gradient

In several marine bivalve species, biochemical genetic variation at a leucine aminopeptidase (LAP) locus is associated with environmental variability, primarily salinity fluctuation. Population genetic variation at a LAP locus was investigated here in two sympatric mussel species (Mytilus galloprovincialis and Perna canaliculus) from three locations along a salinity gradient in Wellington Harbour, New Zealand. The data for M. galloprovincialis and P. canaliculus do not support the hypothesis that the LAP polymorphism in either species is associated with salinity variation among adult mussels. Due to the absence of small mussels among the samples it is not possible to discount the hypothesis that selection acts primarily against juveniles, as it does for M. edulis in Long Island Sound, USA. Wellington Harbour populations of M. galloprovincialis exhibited large and often highly significant heterozygote excesses at the LAP locus, whereas populations of P. canaliculus from the same locations exhibited large and highly significant heterozygote deficiencies. The reason for this inter-specific difference in population structure is unknown. If it is the result of selection, this suggests that selection acts differentially upon the two species, because demographic attributes and reproductive biology are very similar in the two species. For both M. galloprovincialis and P. canaliculus, significant levels of population genetic heterogeneity were recorded among three locations separated by only 8 to 12 km. Neither species exhibited shell length-dependent genetic variation at the LAP locus, suggesting that for these two species the LAP polymorphism is not associated with variation in shell length. Received: 30 December 1996 / Accepted: 6 January 1997     

Lower genetic diversity in the limpet <Emphasis Type="Italic">Patella caerulea</Emphasis> on urban coastal structures compared to natural rocky habitats

Human-made structures are increasingly found in marine coastal habitats. The aim of the present study was to explore whether urban coastal structures can affect the genetic variation of hard-bottom species. We conducted a population genetic analysis on the limpet Patella caerulea sampled in both natural and artificial habitats along the Adriatic coast. Five microsatellite loci were used to test for differences in genetic diversity and structure among samples. Three microsatellite loci showed strong Hardy–Weinberg disequilibrium likely linked with the presence of null alleles. Genetic diversity was significantly higher in natural habitat than in artificial habitat. A weak but significant differentiation over all limpet samples was observed, but not related to the type of habitat. While the exact causes of the differences in genetic diversity deserve further investigation, these results clearly point that the expansion of urban structures can lead to genetic diversity loss at regional scales.     

Wide-spread genetic variability and the paradox of effective population size in the gag,Mycteroperca microlepis,along the West Florida Shelf

Wide-ranging marine species are often described as having a low effective population size (N _e) to census size (N) ratio. This genetic phenomenon is typically attributed to large variation among individuals in reproductive success because of the high mortality rates and unpredictable environments associated with larval dispersal. In this study, we examined patterns of genetic variation in gag (Mycteroperca microlepis) on the West Florida Shelf across year classes of post-settlement juveniles and spawning adults. With no significant genetic differentiation among year classes despite varying recruitment dynamics, little evidence for chaotic genetic patchiness, and no truncation of adult genetic diversity in subsequent juvenile cohorts, there was little support for large variation among individual in reproductive success contributing to a low N _e/N ratio. In fact, the consistent lack of significant differences in annual recruitment classes indicated that reproductive success among individuals was resistant to skewing. Among the various evolutionary forces that may be affecting N _e, changes to demography due to fishing pressure are posited as a likely mechanism affecting current levels of genetic variation.     

Species boundaries in the starfish genus Linckia

 The genetic basis for species boundaries in the starfish genus Linckia was examined using variation observed in 613 base pairs (bp) of sequence from the cytochrome oxidase I gene of mtDNA and 16 allozyme loci. Five groups within Linckia were clearly genetically differentiated; L. columbiae, L. bouvieri, two clades within L. guildingi, and one clade with two sub-clades consisting of both L. laevigata and L. multifora. Genetic divergence among these groups is consistent with interspecific variation. The two clades within L. guildingi suggest the presence of a cryptic, partially sympatric, species. Genetic divergence between these two clades implies that they have been genetically distinct from each other for at least one million years. In contrast, genetic data suggest that L. laevigata and L. multifora are a single species, despite the fact that live individuals can be distinguished by their colour and colour pattern, number of madreporites and ratio of arm length to breadth. There are probably three closely related groups within the L. laevigata/L. multifora clade –L. multifora, and two groups in L. laevigata defined by biogeographic province. It is difficult to determine boundaries for these three entities, since genetic and morphological differences are complicated by phenotypic differences arising from both environmental variation and population genetic structure. The difficulties encountered in defining species boundaries in Linckia, particularly with respect to variation arising from the overlap of Indian and Pacific biogeographic provinces, may be a general issue for many marine organisms from this region. Received: 24 May 1999 / Accepted: 6 October 1999     

Length-weight relationships and coefficient of condition of Euphausia superba and Thysanoessa macrura (Crustacea: Euphausiacea) in southwest Indian Ocean during summer

Samples of antarctic euphausiids, Euphausia superba and Thysanoessa macrura were obtained during the MD 25 FIBEX expedition of the R. V. Marion Dufresne, in February 1981, to the southwest Indian Ocean. Individual and mean coefficients of condition (K _m) were calculated from variations in weight as a function of total length and of carapace length. In E. superba, no significant differences were found as a function of sex or development stage in either total length and weight or carapace length and weight. In T. macrura, a significant allometric difference emerged in the wet weight of juveniles, males and females as a function of total length, and in wet weight as a function of carapace length between adults and subadults. In both E. superba and T. macrura, wet weight as a function of total length differed significantly among stations, and also between species, with a greater weight increase for T. macrura in summer. The coefficients of condition calculated for all wet weight-total length relationships showed that the morphologically different Group II males were heavier than Group I males and mature females. K _m in T. macrura was higher for females than for juveniles or males, indicating a greater weight gain by the females. These differences probably reflect real differences in physiology and may affect the distribution of these two species.     

Interacting Effects of Phenotypic Plasticity and Evolution on Population Persistence in a Changing Climate

Abstract: Climate change affects individual organisms by altering development, physiology, behavior, and fitness, and populations by altering genetic and phenotypic composition, vital rates, and dynamics. We sought to clarify how selection, phenotypic plasticity, and demography are linked in the context of climate change. On the basis of theory and results of recent empirical studies of plants and animals, we believe the ecological and evolutionary issues relevant to population persistence as climate changes are the rate, type, magnitude, and spatial pattern of climate‐induced abiotic and biotic change; generation time and life history of the organism; extent and type of phenotypic plasticity; amount and distribution of adaptive genetic variation across space and time; dispersal potential; and size and connectivity of subpopulations. An understanding of limits to plasticity and evolutionary potential across traits, populations, and species and feedbacks between adaptive and demographic responses is lacking. Integrated knowledge of coupled ecological and evolutionary mechanisms will increase understanding of the resilience and probabilities of persistence of populations and species.     

High levels of gene flow in the surf bivalve Donax deltoides (Bivalvia: Donacidae) on the east coast of Australia

We surveyed patterns of allelic variation within twelve samples of the pipi Donax deltoides Lamarck from beaches separated by up to 1200 km but connected to varying degrees by the East Australian Current. We used these data to test the prediction that the irregular patterns of water movement would cause genetic differentiation in pipis, so that there would be more genetic variation within and among the more southern regions than the northern regions. We found that six loci were at least moderately variable within all samples, and there were no clear geographic patterns in allelic frequencies. In general, genotype frequencies within samples were consistent with predictions for an outcrossed, sexually reproducing species, and we detected no evidence of population subdivision. Within samples, with the exception of the peptidase loci, single-locus genotype frequencies were in close agreement with expectations for Hardy–Weinberg equilibrium. We observed no significant linkage disequilibrium for any pairwise comparison of loci in any sample. Our hierarchical analysis of genetic variation revealed little variation among all samples (F st = 0.009). Loci showed consistently low levels of subdivision (F st from 0.003 to 0.018). We found almost no variation among the four geographic regions sampled (F rt = 0.001). All variation was therefore attributable to variation among samples within regions (F sr = 0.010). These data imply that larvae are moving between regions and that levels of present or recent gene flow are high, and support the conclusions of other studies which have inferred widespread gene flow for animals dispersing via planktonic, outcrossed larvae in parts of this region. This implies that the East Australian Current is sufficient to produce strong larval connections despite its intermittent nature. If existing levels of population subdivision reflect current levels of gene flow, then these data imply that D. deltoides represents a single fishery on the east coast of Australia. Received: 16 September 1996 / Accepted: 25 September 1996     

Invasion patterns inferred from cytochrome oxidase I sequences in three bryozoans, Bugula neritina, Watersipora subtorquata, and Watersipora arcuata

Nucleotide variation in cytochrome c oxidase subunit I (COI) was used to examine population structure in three invasive bryozoans: Bugula neritina (Linnaeus, 1758), Watersipora subtorquata (d’Orbigny, 1852), and W. arcuata (Banta, 1969). These species are found on ship hulls and have a short (≤2 days) larval phase. Samples were collected from 1998–2001 at multiple sites in Australia, and in Hong Kong, New Zealand, Hawaii, California, Curaçao, and England. B. neritina is known to include three cryptic species, including species Type S (Davidson and Haygood in Biol Bull 196:273–280, 1999) which occurs on the east and west coasts of the USA. One haplotype recorded previously in the USA, S1, was found to be widespread, occurring throughout Australia and in Hong Kong, Curaçao, Hawaii, and England. W. subtorquata, a Caribbean–Atlantic species which has invaded southern Australia, New Zealand, and California, had low nucleotide diversity in these areas (π=0.0016±0.0014), consisting of three haplotypes connected by one or two nucleotide mutations. W. arcuata, an Eastern-Pacific native, had comparatively high diversity (π=0.0221±0.0115) in introduced populations from Australia and Hawaii. In each species, identical haplotypes were identified on separate coastlines providing evidence of widespread, rather than genetically independent, introductions. The major contrast in nucleotide diversity suggests that different propagule-source models explain introductions.     

Phylogeography of the Indo-Pacific parrotfish Scarus psittacus: isolation generates distinctive peripheral populations in two oceans

Phylogenetic, phylogeographic, population genetic and coalescence analyses were combined to examine the recent evolutionary history of the widespread Indo-Pacific parrotfish, Scarus psittacus, over a geographic range spanning three marine biogeographic realms. We sequenced 164 individuals from 12 locations spanning 17,000 km, from 55oE to 143oW, using 322 base pairs of mitochondrial control region (D-loop). S. psittacus displayed high haplotype (h = 0.83–0.98), but low nucleotide (<1%) diversity. Most (>83%) genetic variation was within populations. AMOVA revealed significant partitioning and identified five geographic groups. These included one central population and four populations peripheral to the centre. The central population occupied reefs from Western Australia to Tahiti and represented the central Indo-Pacific biogeographic realm. Cocos Keeling was distinct from central and western Indo-Pacific biogeographic realms occupying a position intermediate to these. Peripheral populations (Hawaii, Marquesas) represented the eastern Indo-Pacific biogeographic realm, while Seychelles represented the western Indo-Pacific biogeographic realm. All but the central population expanded (<163 kya). Whilst all populations experienced major sea level and SST changes associated with Pleistocene glaciation cycles, the genetic structure of the central population was relatively homogenous unlike the remaining genetically distinctive populations.     

Subpopulations of the Hawaiian striped mullet Mugil cephalus: analysis of variations of nuclear eye-lens protein electropherograms and nuclear eye-lens weights

The method of electrophoretic separation of nuclear eye-lens proteins was applied to a study of the population structure of the striped mullet Mugil cephalus L. in Hawaii. Electropherograms produced 5 pattern types, with most variations due to non-developmental genetic polymorphism. Some variations, however, were associated with size, and an ontogenetic basis was considered. The frequencies of the pattern types were significantly (p<0.005) independent among the three geographical subpopulations examined: Kaneohe Bay, Oahu; Diamond Head, Oahu; and Kawaihae Bay, Hawaii. Growth rates of mullet eye-lens nuclei demonstrated significant (p<0.05) differences between sexes and localities, except for Kaneohe Bay males versus Diamond Head males, which were nearly significant (p 0.075). These variations in eye-lens nuclei growth rates are probably regulated by genetic factors. Nearly complete genetic isolation is evident for the three subpopulations of mullet in the Hawaiian Islands. The occurrence of genetic interchange in apparently less frequent between the Oahu and Hawaii populations than between the two Oahu populations, in accordance with their respective geographical distances.Oceanic Institute Contribution No. 72.     

Genetic consequences of long larval life in the starfish Linckia laevigata (Echinodermata: Asteroidea) on the Great Barrier Reef

Gene flow between populations of the asteroid Linckia laevigata (Linnaeus) was investigated by examining over 1000 individuals collected from ten reefs throughout the Great Barrier Reef (GBR), Australia, for genetic variation at seven polymorphic enzyme loci. Despite geographic separations in excess of 1000 km, Nei's unbiased genetic distance (0 to 0.003) and standardised genetic variation between populations (F _ST) values (mean 0.0011) were small and not significant. Genetic homogeneity among L. laevigata populations is consistent with the long-distance dispersal capability of its 28 d planktonic larval phase, and is greater than that observed for other asteroid species, including another high-dispersal species, Acanthaster planci, which has a 14 d larval phase. Variation within populations was also higher than previously recorded for asteroids (mean heterozygosity=0.384; number of alleles per locus ranged from 5.1 to 6.0 in each population). Among asteroids, dispersal ability is positively correlated with gene flow and levels of variation, and negatively correlated with levels of differentiation.     

Genetic evidence for the cryptic species pair, <Emphasis Type="Italic">Lottia digitalis </Emphasis>and <Emphasis Type="Italic">Lottia austrodigitalis</Emphasis> and microhabitat partitioning in sympatry

It has been proposed that the common West Coast limpet, Lottia digitalis, is actually the northern counterpart of a cryptic species duo including, Lottia austrodigitalis. Allele frequency differences between southern and northern populations at two polymorphic enzyme loci provided the basis for this claim. Due to lack of further evidence, L. austrodigitalis is still largely unrecognized in the literature. Seven additional enzyme loci were examined from populations in proposed zones of allopatry and sympatry to determine the existence of L. austrodigitalis as a sibling species to L. digitalis. Significant allele frequency differences were found at five enzyme loci between populations in Laguna Beach, southern California, and Bodega Bay, northern California; strongly supporting the existence of separate species. Both species exhibit two microhabitat morphotypes, a gooseneck barnacle morph in the mid-intertidal zone and a rock morph in the high-intertidal zone. In sympatry, L. austrodigitalis was more abundant higher in the intertidal on rocks, whereas L. digitalis was more abundant lower in the intertidal on barnacles. This finding supports earlier claims of microhabitat partitioning in this sibling species pair. In addition to this finding, the transition zone between the species was found to have shifted substantially northward in only two decades, from Monterey Peninsula, CA to near Pigeon Point, CA, where L. digitalis previously dominated.     

Comparative population genetic structure of marine gastropods (Littorina spp.) with and without pelagic larval dispersal

 Population genetic theory predicts that marine animal species with planktonic larvae will have less genetic structure than those with direct development. We compared the genetic structure of four species of littorinid snails – two with planktonic egg capsules that hatch as planktonic larvae and two with benthic egg masses that hatch as crawl-away juveniles. We used DNA sequencing and single stranded conformational polymorphism (SSCP) to assess sequence variation in a 480 bp fragment of the mitochondrial cytochrome b gene and then used an analysis of molecular variance (AMOVA) to estimate Φ_st for populations from the northeastern Pacific coast. One of the two direct-developing species, Littorina subrotundata, had a moderate amount of population structure (Φ_st=0.209) as expected but the other direct-developing species, L. sitkana, was nearly fixed for a single haplotype that made it impossible to precisely estimate Φ_st. One of the two planktonic-developing species, L. scutulata, did not show any significant population structure (Φ_st=0.004). In contrast to our expectations, the other planktonic-developing species, L. plena, showed some weak but statistically significant population structure (Φ_st=0.052). We discuss how differences in population genetic structure between species with the same type of development may reflect differences in their historical demography. Received: 22 December 1999 / Accepted: 24 July 2000     

Genetic differentiation between populations of Talitrus saltator and Talorchestia deshayesii (Crustacea: Amphipoda) from coastal areas of the north-western European continent

Electrophoretic studies of gene-enzyme variation in the littoral talitrids Talitrus saltator (Montagu) and Talorchestia deshayesii (Audouin) were undertaken to estimate the amount of divergence among geographically separated populations. Samples of both species were taken from sandy beaches over a transect of approximately 3 500 km along the coast of the European continent including Baltic, North Sea and Atlantic locations. A total of 22 T. saltator and 15 T. deshayesii populations were analysed for genetic variation at various enzyme loci. Both amphipods revealed relatively low levels of polymorphism and heterozygosity. Among the loci studied, phosphoglucose isomerase (Pgi) and phosphoglucomutase (Pgm) were highly polymorphic. Patterns of micro- and macrogeographic variation in terms of distributions of allele frequencies at these particular loci are compared. Interpopulation allozymic variation was shown to be lower in T. deshayesii than in T. saltator. As demonstrated by T. saltator populations sampled in coastal sites ranging from Denmark to western France, clinal variation in frequencies of two alleles became evident at the PGI locus, exhibiting a steady increase in the level of polymorphism with decreasing latitudes. It is argued that limited gene flow and, to some extent, random genetic drift may account for the gene pool structure of the talitrid species investigated.     

Long-term food modification by Acartia clausi: A preliminary view

After acclimation, the copepod Acartia clausi was allowed to graze for 5 days in a mixed suspension of two discrete size classes (species) of the chaining diatoms Thalassiosira spp. derived from continuous culture. Total particle numbers and particle size distributions of Thalassiosira spp. were stable throughout the 5 days, indicating that the effects of algal removal and modification due to grazing were balanced by algal growth. Grazer ingestion is the predominant process affecting all size classes of the smaller diatom population (T. nordenskioldii); however, both ingestion and chain modification are observed with the larger diatom (T. gravida). Although the greatest percentage removal of algal volume occurs in the largest algal size classes for each algal species, the greatest volume removal occurs at the modal peak (T. nordenskioldii) or just to the right of the modal peak (T. gravida). Flask-to-flask replicability of experiments was tested with the single-celled T. fluviatilis, and these tests were compared to the long-term experiments with T. nordenskioldii and T. gravida. Net particle removal occurs on both the large and small algal species in mixed suspension, not just on the larger-cell-side of the T. gravida distribution. Although 80% of the particles ingested are from the T. nordenskioldii distribution, 80% of the volume ingested is from the T. gravida distribution. The apparent multi-peak selection observed in our data is discussed in reference to two separate hypotheses and in light of other recent work pertaining to selective grazing by copepods.