The effect of genetic and environmental variation on metabolic gene expression

What is the relationship between genetic or environmental variation and the variation in messenger RNA (mRNA) expression? To address this, microarrays were used to examine the effect of genetic and environmental variation on cardiac mRNA expression for metabolic genes in three groups of Fundulus heteroclitus: (i) individuals sampled in the field (field), (ii) field individuals acclimated for 6 months to laboratory conditions (acclimated), or (iii) individuals bred for 10 successive generations in a laboratory environment (G10). The G10 individuals have significantly less genetic variation than individuals obtained in the field and had a significantly lower variation in mRNA expression across all genes in comparison to the other two groups (P = 0.001). When examining the gene specific variation, 22 genes had variation in expression that was significantly different among groups with lower variation in G10 individuals than in acclimated individuals. Additionally, there were fewer genes with significant differences in expression among G10 individuals vs. either acclimated or field individuals: 66 genes have statistically different levels of expression vs. 107 or 97 for acclimated or field groups. Based on the permutation of the data, these differences in the number of genes with significant differences among individuals within a group are unlikely to occur by chance (P < 0.01). Surprisingly, variation in mRNA expression in field individuals is lower than in acclimated individuals. Relative to the variation among individual within a group, few genes have significant differences in expression among groups (seven, 2.3%) and none of these are different between acclimated and field individuals. The results support the concept that genetic variation affects variation in mRNA expression and also suggests that temporal environmental variation associated with estuarine environments does not increase the variation among individuals or add to the differences among groups.     

Multivariate analyses of morphometrical features from Apatemon gracilis (Rudolphi, 1819) Szidat, 1928 and A. annuligerum (v. Nordmann, 1832) (Digenea: Strigeidae) metacercariae

Only two species of the strigeid subgenus Apatemon (Apatemon) sensu Sudarikov (1959), Apatemon gracilis (Rudolhi, 1819) and A. annuligerum (von Nordmann, 1832), are known to infect British fishes. A. annuligerum is poorly characterised, with few life-cycle data available, since only the adult and metacercarial stages have been described. Both of these life-stages demonstrate only minor morphological differences from those of A. gracilis. The multidimensional approach of principal components analysis (PCA) was employed to determine the intra- and inter-specific morphological variation that exists for metacercariae of these species. Analyses revealed that metacercarial morphometrics cannot confidently discriminate between Apatemon gracilis and A. annuligerum specimens. Indeed, levels of intra-specific variation present between A. gracilis metacercariae originating from different piscine hosts were as marked for certain populations as those between the two nominal species. Individual populations of A. gracilis metacercariae differed most in the relative dimensions of their internal organs, whereas variation between the two nominal species was primarily in length-related variables. A. gracilis specimens originating from Welsh and Scottish stoneloach represented two populations with the most comparable total body lengths, yet clustered separately. Conversely, A. gracilis metacercariae infecting rainbow trout and salmon parr differed markedly in size, but demonstrated minimal composite morphological variation. Much, if not all, of the variation described between populations and species of parasites can be attributed to host or to metacercarial location within a particular host. These findings and the lack of strong life-cycle data for A. annuligerum place some doubt on its validity as a species discrete from A. gracilis.     

MOLECULAR GENETIC EVIDENCE FOR PARALLEL LIFE-HISTORY EVOLUTION WITHIN A PACIFIC SALMON (SOCKEYE SALMON AND KOKANEE,ONCORHYNCHUS NERKA)

The Pacific salmon Oncorhynchus nerka typically occurs as a sea-run form (sockeye salmon) or may reside permanently in lakes (kokanee) thoughout its native North Pacific. We tested whether such geographically extensive ecotypic variation resulted from parallel evolutionary divergence thoughout the North Pacific or whether the two forms are monophyletic groups by examining allelic variation between sockeye salmon and kokanee at two minisatellite DNA repeat loci and in mitochondrial DNA (mtDNA) Bgl II restriction sites. Our examination of over 750 fish from 24 populations, ranging from Kamchatka to the Columbia River, identified two major genetic groups of North Pacific O. nerka: a “northwestern” group consisting of fish from Kamchatka, western Alaska, and northwestern British Columbia, and a “southern” group consisting of sockeye salmon and kokanee populations from the Fraser and Columbia River systems. Maximum-likelihood analysis accompanied by bootstrapping provided strong support for these two genetic groups of O. nerka; the populations did not cluster by migratory form, but genetic affinities were organized more strongly by geographic proximity. The two major genetic groups resolved in our study probably stem from historical isolation and dispersal of O. nerka from two major Wisconsinan glacial refugia in the North Pacific. There were significant minisatellite DNA allele frequency differences between sockeye salmon and kokanee populations from different parts of the same watershed, between populations spawning in different tributaries of the same lake, and also between sympatric populations spawning in the same stream at the same time. MtDNA Bgl II restriction site variation was significant between sockeye salmon and kokanee spawning in different parts of the same major watershed but not between forms spawning in closer degrees of reproductive sympatry. Patterns of genetic affinity and allele sharing suggested that kokanee have arisen from sea-run sockeye salmon several times independently in the North Pacific. We conclude that sockeye salmon and kokanee are para- and polyphyletic, respectively, and that the present geographic distribution of the ecotypes results from parallel evolutionary origins of kokanee from sockeye (divergences between them) thoughout the North Pacific.     

MULTIVARIATE ADAPTATION BUT NO INCREASE IN COMPETITIVE ABILITY IN INVASIVE GERANIUM CAROLINIANUM L. (GERANIACEAE)

Adaptive evolution can affect the successful establishment of invasive species, but changes in selective pressures, loss of genetic variation in relevant traits, and/or altered trait correlations can make adaptation difficult to predict. We used a common‐garden experiment to assess trait correlations and patterns of adaptation in the invasive plant, Geranium carolinianum, sampled across 20 populations in its native (United States) and invasive (China) ranges. We used multivariate Q_ST – F_ST tests to determine if phenotypic differences between countries are attributable to adaptation. We also compared population‐level variation within each country to assess whether local adaptation resulted in similar multivariate phenotypes in the United States and China. Between countries, most phenotypic differences are indistinguishable from genetic drift, although we detected a signature of adaptation to the colder, drier winters in China. There was no evidence for increases in invasive traits in China. Within countries, strong multivariate adaptation appears to be driven by latitudinal climatic variation in the United States, but not in China. Additionally, adaptive trait combinations as well as their underlying correlations differ between the two countries, indicating that adaptation in invasive populations does not parallel patterns in native populations due to differences in selection pressures, genetic constraints, or both.     

Apochromic populations of Dactylorhiza incarnata s.l. (Orchidaceae): diversity and systematic significance as revealed by allozyme markers and morphology

Apochromic forms of the Eurasian Dactylorhiza incarnata s.l. were studied in northern Europe to reveal their genetic (allozyme) and morphological diversity and to assess their systematic significance. The study included eight localities with sympatric populations of plants with anthocyanin‐pigmented and apochromic flowers. Parallel samples of the two morphs were taken from each locality. Genetic variation was only found at the allozyme loci pgd, pgi and ugpp. Statistically significant differences in allele frequencies between the two colour morphs were found in two localities and demonstrate that the occurrence of apochromic individuals in D. incarnata s.l. is not always because of spontaneous mutation. At least in some localities the apochromic plants form distinct breeding groups (but local populations of different colour morphs may also be composed of several more or less distinct breeding groups). Based on molecular and morphometric data, it is proposed that the apochromic study populations from calcareous fens should be referred to D. incarnata var. ochroleuca, whereas the apochromic study populations from non‐calcareous fens are better treated as aberrant local populations of var. incarnata s.l. Possible evolutionary patterns and processes are discussed and guidelines for identification of var. ochroleuca from morphological features are given. © 2009 The Linnean Society of London, Botanical Journal of the Linnean Society, 2009, 159 , 396–407.     

Morphometric differentiation of the moor frog (Rana arvalis Nilss.) in Central Europe

Morphometric variation of 11 characters was studied in 14 samples of Rana arvalis from allopatric population groups in Central Europe, representing two subspecies formerly recognised: R. a. arvalis and R. a. wolterstorffi. All samples from Poland (nine) were collected from the area, which is believed to be populated by the nominal form. In addition, the isolated population of R. arvalis from the Eastern Carpathians of Romania was classified to that form. All the Hungarian samples fall into the range of R. a. wolterstorffi. No significant differences in body size between sexes and among the populations from the three geographic regions (Poland, Hungary and Romanian Eastern Carpathians) were found. Multivariate analysis of the differences in body shape among groups showed that the populations from Hungary differed from the nominal form mainly in the relative hindlimb length. However, substantial overlap in overall shape differences resulted in a relatively low percentage of correct classifications to the respective geographic groups in the Discriminant Analysis. Large intrapopulation variation in the values of the two ratios (SVL/TL and TL/IMT) formerly used for differentiate the two European subspecies of R. arvalis makes them unreliable characters for distinguishing the two forms. The pattern of morphometric variation does not match the overall genetic divergence of the R. arvalis populations, which suggests that the body shape differences in this species result from the phenotypic plasticity correlated with local climatic factors.     

Little evidence for genetic variation associated with susceptibility to sea star wasting syndrome in the keystone species Pisaster ochraceus

The keystone species Pisaster ochraceus suffered mass mortalities along the northeast Pacific Ocean from Sea Star Wasting Syndrome (SSWS) outbreaks in 2013–2016. SSWS causation remains of debate, leading to concerns as to whether outbreaks will continue to impact this species. Considering the apparent link between ocean temperature and SSWS, the future of this species and intertidal communities remains uncertain. Surveys of co-occurring apparently normal and wasting P. ochraceus along the central Oregon coast in 2016 allowed us to address whether variation in disease status showed genetic variation that may be associated with differences in susceptibility to SSWS. We performed restriction site-associated DNA sequencing (2bRAD-seq) to genotype ~72,000 single nucleotide polymorphism (SNP) loci across apparently normal and wasting sea stars. Locus-specific analyses of differentiation (F_ST) between disease-status groups revealed no signal of genetic differences separating the two groups. Using a multivariate approach, we observed weak separation between the groups, but identified 18 SNP loci showing highest discriminatory power between the groups and scanned the genome annotation for linked genes. A total of 34 protein-coding genes were found to be located within 15 kb (measured by linkage disequilibrium decay) of at least one of the 18 SNPs, and 30 of these genes had homologies to annotated protein databases. Our results suggest that the likelihood of developing SSWS symptoms does not have a strong genetic basis. The few genomic regions highlighted had only modest levels of differentiation, but the genes associated with these regions may form the basis for functional studies aiming to understand disease progression.     

Variation in genetic structure among populations of the caddisfly Helicopsyche borealis from three streams in northern California, U.S.A.

1. Allozyme electrophoresis was used to describe the genetic structure of Helicopsyche borealis caddisflies collected from three sites in each of three streams in northern California, U.S.A.: the Rice Fork of the Eel River; Big Sulphur Creek; and Alameda Creek. Between 7 and 11km separated adjacent sites within these three streams. Helicopsyche borealis from three additional streams in eastern North America (Christiana Creek, Indiana; Byrd's Mill Creek, Oklahoma; Susquehanna River, Pennsylvania) were also analysed electrophoretically to address taxonomic questions that arose during the study. 2. Four genetically distinct groups of individuals were identified (i.e. Helicopsyche types A, B, C, and D). Lack of interbreeding between sympatric groups (as evidenced by fixed allelic differences) and large genetic differences (mean Nei's genetic distances = 0.396–0.693) indicate that these four groups of Helicopsyche were actually reproductively isolated species rather than genetic variants of a single species. 3. Occurrence of Helicopsyche type A at multiple sites permitted an analysis of spatial variation in genetic structure. Within a drainage basin, small differences in allele frequencies were observed among sites in the Rice Fork and Big Sulphur Creek, but not in Alameda Creek. Larger genetic differences were found among sites in separate drainage basins. Genetic distances (Nei's) between Helicopsyche type A from California and from eastern North America sites (mean = 0.236) were greater than interpopulation differences commonly observed for insects, which suggests that Helicopsyche type A from California may represent a different species than Helicopsyche type A from eastern North America. 4. Geographical and taxonomic differences observed in this study underscore the importance of understanding both population structure and genetic relationships among populations in the design and interpretation of stream faunal studies.     

Experimental studies of adaptation in Clarkia xantiana. I. Sources of trait variation across a subspecies border

Abstract Both genetic differentiation and phenotypic plasticity might be expected to affect the location of geographic range limits. Co‐gradient variation (CoGV), plasticity that is congruent with genetic differentiation, may enhance performance at range margins, whereas its opposite, counter‐gradient variation (CnGV) may hinder performance. Here we report findings of reciprocal transplant experiments intended to tease apart the roles of differentiation and plasticity in producing phenotypic variation across a geographic border between two plant subspecies. Clarkia xantiana ssp. xantiana and C. xantiana ssp. parviflora are California‐endemic annuals that replace each other along a west‐east gradient of declining precipitation. We analyzed variation in floral traits, phenological traits, and vegetative morphological and developmental traits by sowing seeds of 18 populations (six of ssp. xantiana and 12 of ssp. parviflora) at three sites (one in each subspecies' exclusive range and one in the subspecies' contact zone), in two growing seasons (an exceptionally wet El Niño winter and a much drier La Niña winter). Significant genetic differences between subspecies appeared in 11 of 12 traits, and differences were of the same sign as in nature. These findings are consistent with the hypothesis that selection is responsible for subspecies differences. Geographic variation within subspecies over part of the spatial gradient mirrored between‐subspecies differences present at a larger scale. All traits showed significant plasticity in response to spatial and temporal environmental variation. Plasticity patterns ranged from spatial and temporal CoGV (e.g., in node of first flower), to spatial CnGV (e.g., in flowering time), to patterns that were neither CoGV nor CnGV (the majority of traits). Instances of CoGV may reflect adaptive plasticity and may serve to increase performance under year‐to‐year environmental variation and at sites near the subspecies border. However, the presence of spatial CnGV in some critical traits suggests that subspecies ranges may also be constrained by patterns of plasticity.     

Further evidence for geographic differentiation in R. appendiculatus (Acari: Ixodidae) from Eastern and Southern provinces of Zambia

Studies in the biology, ecology and behaviour of R. appendiculatus in Zambia have shown considerable variation within and between populations often associated with their geographical origin. We studied variation in the mitochondrial COI (mtCOI) gene of adult R. appendiculatus ticks originating from the Eastern and Southern provinces of Zambia. Rhipicephalus appendiculatus ticks from the two provinces were placed into two groups on the mtCOI sequence data tree. One group comprised all haplotypes of specimens from the Eastern province plateau districts of Chipata and Petauke. The second group consisted of a single haplotype of specimens from the Southern province districts and Nyimba, an Eastern province district on the fringes of the valley. This variation provides additional evidence to the earlier observations in the 12S rDNA and ITS2 data for the geographic subdivision of R. appendiculatus from Southern province and Eastern province plateau. The geographic subdivision further corresponds with differences in body size and diapause between R. appendiculatus from these geographic areas. The possible implications of these findings on the epidemiology of East Coast fever (ECF) the disease for which R. appendiculatus is one of the vectors are discussed. Nucleotide sequence data reported in this paper are available in the GenBank^TM under the Accession Numbers DQ859259–DQ859266     

VARIATION IN GENETIC ARCHITECTURE OF CALLING SONG AMONG POPULATIONS OF ALLONEMOBIUS SOCIUS,A. FASCIATUS,AND A HYBRID POPULATION: DRIFT OR SELECTION?

Predictions using quantitative genetic models generally assume that the variance-covariance matrices remain constant over time. This assumption is based on the supposition that selection is generally weak and hence variation lost through selection can be replaced by new mutations. Whether this is generally true can only be ascertained from empirical studies. Ideally for such a study we should be able to make a prediction concerning the relative strength of selection versus genetic drift. If the latter force is prevalent then the variance-covariances matrices should be proportional to each other. Previous studies have indicated that females in the two sibling cricket species Allonemobius socius and A. fasciatus do not discriminate between males of the two species by their calling song. Therefore, differences between the calling song of the two males most likely result from drift rather than sexual selection. We test this hypothesis by comparing the genetic architecture of calling song of three populations of A. fasciatus with two populations of A. socius. We found no differences among populations within species, but significant differences in the G (genetic) and P (phenotypic) matrices between species, with the matrices being proportional as predicted under the hypothesis of genetic drift. Because of the proportional change in the (co)variances no differences between species are evident in the heritabilities or genetic correlations. Comparison of the two species with a hybrid population from a zone of overlap showed highly significant nonproportional variation in genetic architecture. This variation is consistent with a general mixture of two separate genomes or selection. Qualitative conclusions reached using the phenotypic matrices are the same as those reached using the genetic matrices supporting the hypothesis that the former may be used as surrogate measures of the latter.     

MOLECULAR AND MORPHOLOGICAL ANALYSIS OF ENTEROMORPHA INTESTINALIS AND E. COMPRESSA (CHLOROPHYTA) IN THE BRITISH ISLES

The very common green seaweeds Enteromorpha intestinalis (L.) Nees and E. compressa (L.) Nees are important fouling organisms and have commonly been used as indicators of eutrophication, but their taxonomic status is problematic. The genus presents extreme difficulties because there is wide intraspecific variation in morphology, but morphological differences between species are small and difficult to detect. In this study, molecular data were used in parallel with morphological characters to resolve the taxonomic problems. Phylogenetic analysis of sequences of the internal transcribed spacers ITS1 and ITS2 and the 5.8S gene distinguished two groups of samples, which were identified by morphological characters as E. compressa (branched) and E. intestinalis (normally unbranched). There was a low level of sequence divergence within each group of samples, but divergence between groups was as great as that between either of the two species and the outgroup E. prolifera. Clades representing E. compressa and E. intestinalis were also found in analyses of an independent molecular data set, chloroplast DNA restriction fragment length polymorphisms (RFLPs). Enteromorpha intestinalis and E. compressa represent two distinct, genetically divergent species. Reinterpretation of published studies shows that these species are reproductively isolated. However, E. compressa and E. intestinalis are sometimes very difficult to distinguish from each other and could be regarded as cryptic species. The presence or absence of branching was the most useful character distinguishing these two species, but there was an element of ambiguity because low salinity or salinity shock can induce branching in E. intestinalis. If environmental factors such as salinity are taken into account, branching can be used to identify the great majority of thalli correctly. This study therefore provides a basis for identifying the two most important marine fouling macroalgae and for their use in environmental monitoring and experimentation. Typification of these two Linnaean species showed that current usage of the names accords with the lectotype and protologue of both species. Samples that resembled E. usneoides did not form a clade in any of the trees, and constraining the data to support the monophyly of this group incurred a penalty. Enteromorpha usneoides appears to be an ecotype of E. compressa.     

Numerical morphological analysis of Lettuce cultivars and species (Lactuca sect.Lactuca,Asteraceae)

L. sativa L. cultivars were compared morphologically with accessions ofL. saligna L.,L. serriola L., andL. virosa L. in order to define the infra- and interspecific variation. Multivariate analysis (principal component analysis) distinguished four groups corresponding with the four species. Infraspecific variation (e.g., colour, leaf shape) withinL. sativa is the largest, because breeding has caused a large variation in vegetative features; infraspecific variation withinL. serriola is also large. Three plants of one accession were placed betweenL. sativa andL. serriola because they showed morphological characteristics of both species. They are still classified asL. sativa. The morphological differences betweenL. sativa andL. serriola are too large to consider these two as one species. The descent of the cultivar groups and the four species is described.     

Intraspecific and interspecific variation of female genitalia in two species of watersnake

The morphological differences in female genitalia within and between species are little studied and poorly understood, yet understanding patterns of variation in female genitalia can provide insights into mechanisms of genital evolution. The present study aimed to explore the patterns of intraspecific and interspecific variation in female genitalia in two sister taxa of watersnake (Nerodia sipedon and Nerodia fasciata) that have similar genital shape. We used a geometric morphometric (GM) approach to study variation in shape of the vagina between and within two sister species. We examined genital shape in female watersnakes ranging from small, sexually immature females to large reproductive females that had recently given birth. We found that shape variation of genitalia is strongly correlated with body size, where larger but not smaller females have a bifurcation in the vagina. However, we also found significant shape variation in the structure of the vagina between the two species, where N. fasciata has narrower genitalia with more prominent bifurcation, whereas N. sipedon has wider genitalia with less marked bifurcation. Using GM allowed us to detect significant differences in genital shape that were not apparent upon visual examination alone, suggesting that shape variation in female genitalia may be greater than previously assumed. Additional study of morphological differences in male reproductive organs for these species would help to determine whether there has been genital co‐evolution, and potentially mechanical reproductive isolation, in these two closely‐related and occasionally sympatric species. © 2013 The Linnean Society of London, Biological Journal of the Linnean Society, 2014, 111, 183–191.     

Morphological divergence of lake and stream Phoxinus of Northern Italy and the Danube basin based on geometric morphometric analysis

Minnows of the genus Phoxinus are promising candidates to investigate adaptive divergence, as they inhabit both still and running waters of a variety of altitudes and climatic zones in Europe. We used landmark‐based geometric morphometric methods to quantify the level of morphological variability in Phoxinus populations from streams and lakes of Northern Italy and the Danube basin. We analyzed body shape differences of populations in the dorsal, lateral, and ventral planes, using a large array of landmarks and semilandmarks. As the species identification of Phoxinus on morphological characters is ambiguous, we used two mitochondrial genes to determine the genetic background of the samples and to ensure we are comparing homogenous groups. We have found significant body shape differences between habitats: Minnow populations inhabiting streams had a deeper body and caudal peduncle and more laterally inserted pectoral fins than minnows inhabiting lakes. We have also found significant body shape differences between genetic groups: Italian minnows had deeper bodies, deeper and shorter caudal peduncles, and a shorter and wider gape than both groups from the Danube. Our results show that the morphology of Phoxinus is highly influenced by habitat and that body shape variation between habitats was within the same range as between genetic groups. These morphological differences are possibly linked to different modes of swimming and foraging in the respective habitats and are likely results of phenotypic plasticity. However, differences in shape and interlandmark distances between the groups suggest that some (though few) morphometric characters might be useful for separating Phoxinus species.     

Quantitative genetic analysis of cranial morphology in the cotton-top (Saguinus oedipus) and saddle-back (S. fuscicollis) tamarins

Patterns of morphological variation play an important role in evolutionary diversification and are critical to an informed interpretation of interspecific differences. When patterns of genetic variation have not diverged substantially, it is possible to reconstruct the differences in selection which gave rise to morphological differences among extant species. Morphological variation patterns are compared between two tamarin species, the cotton-top tamarin (Saguinus oedipus) and the saddle-back tamarin (S. fuscicollis illigeri). Genetic, phenotypic, and environmental variance/covariance and correlation matrices were obtained for a series of 39 cranial characters in each species (cotton-top tamarin, N = 328; saddle-back tamarin, N = 209) and for the species combined using crania from individuals of known genealogical relationship. After accounting for the effects of estimation error on measures of matrix similarity, patterns of phenotypic, genetic, and environmental variation and correlation were found to be very similar across species and among the types of variance within species. Taking the saddle-back tamarins as the standard, cotton-top tamarins have been selected for an enlarged anterior temporalis attachment area and increased facial prognathism. In primates, an enlarged anterior temporalis muscle is associated with incisive food preparation, especially at wide gape.     

Population dynamics and comparative demographics in sympatric populations of the round-leaved orchids Platanthera macrophylla and P. orbiculata

Orchids (Orchidaceae) are a family of flowering plants with a high proportion of threatened taxa making them an important focus of plant conservation. Orchid conservation efforts are most effective when informed by reliable demographic research. We utilized transition matrix models to examine the population dynamics and demography within sympatric populations of a species pair of terrestrial round-leaved orchids, Platanthera macrophylla and P. orbiculata. The models were parameterized from a large data set spanning 9 years from field observations of over 1,000 orchids. Life table response experiments (LTRE) were used to identify which life history transitions, and which vital rates within those transitions, most contributed to observed differences between the two species and most contributed to interannual variation within each species. Results from mean transition matrices projected finite rates of population growth that were not significantly different between the two species, with P. macrophylla near the replacement rate and P. orbiculata below it. LTRE revealed that the difference in population growth rates between the two species was mostly due to differences in fecundity (flowering adult to protocorm transition) driven by differences in fruit set and seed germination into protocorm, which were much greater for P. macrophylla. The two primary contributors to interannual variation in population growth rates for both orchids were adult survival and fruit set, respectively. These findings indicate that any environmental disturbances harming adult survival or fecundity will have a disproportionately negative effect on the orchid populations.     

Genetic variation within and among North Atlantic and Baltic populations of the benthic alga Phycodrys rubens (Rhodophyta)

Genetic variation was examined within and among North Atlantic, North Sea and Baltic populations of the benthic red alga Phycodrys rubens using allozymes and random amplified polymorphic DNA (RAPD) markers. On western and eastern North Atlantic coasts distinct allozyme types were found, with the exception of western Newfoundland where East and West Atlantic types co-occur. Along the European coasts, two genetic groups were distinguished by fixed allelic differences: an outer oceanic group and a North Sea/Baltic group. The two genetic types co-occur in the Skagerrak and Kattegat region. Reproductive isolation between the two types is suggested by the lack of hybrids in the overlap zones, and they may therefore represent sibling species. Unexpectedly, an analysis of RAPD variation was unable to recover the two cryptic species identified using allozymes. Within-population RAPD variation was similar to or greater than between-population variation. The lack of structure in the RAPD data cannot be attributed solely to technical artefacts of the method but appears to reflect real biological variability. Within-population genomic polymorphisms caused by frequent mutational events are discussed, as are high amounts of genetic drift and possible disruptive selection brought about by stressed habitats. Finally, Baltic and extra-Baltic salinity ecotypes are known to exist in P. rubens. However, no correlation between ecotypic variation and allozyme groups was detected.     

GENETIC, MORPHOLOGICAL, AND TOXICOLOGICAL VARIATION AMONG GLOBALLY DISTRIBUTED STRAINS OF NODULARIA (CYANOBACTERIA)

Morphological, toxicological, and genetic variation was examined among 19 strains of Nodularia. The strains examined could be morphologically discriminated into four groups corresponding to N. spumigena Mertens, N. sphaerocarpa Bornet et Flahault, and two strains that did not clearly correspond to currently accepted Nodularia species. Genetic variation was examined using nucleotide sequencing of the phycocyanin intergenic spacer region (cpcBA-IGS) and RAPD-PCR. The PCR-RFLP of the cpcBA-IGS differentiated four genotypes corresponding to the four morphological groups. However, nucleotide sequencing of 598 bp of the 690-bp fragment showed that one of the three strains corresponding to N. sphaerocarpa (PCC 7804) was genetically divergent from the other two, suggesting that it constitutes a distinct species. Nucleotide variation within the morphospecies groups was limited (<1%), and all 14 Australian strains of N. spumigena possessed identical cpcBA-IGS sequences. The RAPD-PCR differentiated the same groups as the cpcBA sequencing and discriminated each of the seven different Australian populations of N. spumigena. Strains from within a bloom appeared genetically identical; however, strains isolated from different blooms could be separated into either a western or a southeastern Australian cluster, with one strain from western Australia showing considerable genetic divergence. The pattern of variation suggests that individual blooms of N. spumigena are clonal but also that Australian N. spumigena populations are genetically distinct from each other. Examination of genetic distance within and between blooms and within and between morphological groups showed clear genetic dicontinuities that, in combination with the cpcBA-IGS data, suggest that Nodularia contains genetically distinct morphospecies rather than a continuous cline of genetic variation. Furthermore, these morphospecies are genetically variable, exhibiting hierarchical patterns of genetic variation on regional and global scales. Production of the hepatotoxin nodularin was not restricted to one genetic lineage but was distributed across three of the five genotypic groups. A strain of N. spumigena from a nontoxic Australian population was found to fall within the range of genetic variation for other toxic Australian strains and appears to be a unique nontoxic strain that might have arisen by loss of toxin production capacity.