Adaptive divergence despite strong genetic drift: genomic analysis of the evolutionary mechanisms causing genetic differentiation in the island fox (Urocyon littoralis)

The evolutionary mechanisms generating the tremendous biodiversity of islands have long fascinated evolutionary biologists. Genetic drift and divergent selection are predicted to be strong on islands and both could drive population divergence and speciation. Alternatively, strong genetic drift may preclude adaptation. We conducted a genomic analysis to test the roles of genetic drift and divergent selection in causing genetic differentiation among populations of the island fox (Urocyon littoralis). This species consists of six subspecies, each of which occupies a different California Channel Island. Analysis of 5293 SNP loci generated using Restriction‐site Associated DNA (RAD) sequencing found support for genetic drift as the dominant evolutionary mechanism driving population divergence among island fox populations. In particular, populations had exceptionally low genetic variation, small N_e (range = 2.1–89.7; median = 19.4), and significant genetic signatures of bottlenecks. Moreover, islands with the lowest genetic variation (and, by inference, the strongest historical genetic drift) were most genetically differentiated from mainland grey foxes, and vice versa, indicating genetic drift drives genome‐wide divergence. Nonetheless, outlier tests identified 3.6–6.6% of loci as high F_ST outliers, suggesting that despite strong genetic drift, divergent selection contributes to population divergence. Patterns of similarity among populations based on high F_ST outliers mirrored patterns based on morphology, providing additional evidence that outliers reflect adaptive divergence. Extremely low genetic variation and small N_e in some island fox populations, particularly on San Nicolas Island, suggest that they may be vulnerable to fixation of deleterious alleles, decreased fitness and reduced adaptive potential.     

Differential introgression from a sister species explains high FST outlier loci within a mussel species

Scanning genomes for loci with high levels of population differentiation has become a standard of population genetics. F_ST outlier loci are most often interpreted as signatures of local selection, but outliers might arise for many other reasons too often left unexplored. Here, we tried to identify further the history and genetic basis underlying strong differentiation at F_ST outlier loci in a marine mussel. A genome scan of genetic differentiation has been conducted between Atlantic and Mediterranean populations of Mytilus galloprovincialis. The differentiation was low overall (F_ST = 0.03), but seven loci (2%) were strong F_ST outliers. We then analysed DNA sequence polymorphism at two outlier loci. The genetic structure proved to be the consequence of differential introgression of alleles from the sister‐hybridizing species Mytilus edulis. Surprisingly, the Mediterranean population was the most introgressed at these two loci, although the contact zone between the two species is nowadays localized along the Atlantic coasts of France and the British Isles. A historical contact between M. edulis and Mediterranean M. galloprovincialis should have happened during glacial periods. It proved difficult to disentangle two hypotheses: (i) introgression was adaptive, implying edulis alleles have been favoured in Mediterranean populations, or (ii) the genetic architecture of the barrier to edulis gene flow is different between the two M. galloprovincialis backgrounds. Five of the seven outliers between M. galloprovincialis populations were also outliers between M. edulis and Atlantic M. galloprovincialis, which would support the latter hypothesis. Differential introgression across semi‐permeable barriers to gene flow is a neglected scenario to interpret outlying loci that may prove more widespread than anticipated.     

Phylogeography of the California sheephead,Semicossyphus pulcher: the role of deep reefs as stepping stones and pathways to antitropicality

In the past decade, the study of dispersal of marine organisms has shifted from focusing predominantly on the larval stage to a recent interest in adult movement. Antitropical distributions provide a unique system to assess vagility and dispersal. In this study, we have focused on an antitropical wrasse genus, Semicossyphus, which includes the California sheephead, S. pulcher, and Darwin's sheephead, S. darwini. Using a phylogenetic approach based on mitochondrial and nuclear markers, and a population genetic approach based on mitochondrial control region sequences and 10 microsatellite loci, we compared the phylogenetic relationships of these two species, as well as the population genetic characteristics within S. pulcher. While S. pulcher and S. darwini are found in the temperate eastern Pacific regions of the northern and southern hemispheres, respectively, their genetic divergence was very small (estimated to have occurred between 200 and 600 kya). Within S. pulcher, genetic structuring was generally weak, especially along mainland California, but showed weak differentiation between Sea of Cortez and California, and between mainland California and Channel Islands. We highlight the congruence of weak genetic differentiation both within and between species and discuss possible causes for maintenance of high gene flow. In particular, we argue that deep and cooler water refugia are used as stepping stones to connect distant populations, resulting in low levels of genetic differentiation.     

Diversity among peripheral populations: genetic and evolutionary differentiation of Salamandra atra at the southern edge of the Alps

Separate populations at the edge of a species range are receiving great attention and have been shown to be often different from populations in the core area. However, it has rarely been tested whether neighboring peripheral populations are genetically and evolutionarily similar to each other, as expected for their geographical proximity and similar ecological conditions, or differ due to historical contingency. We investigated isolation and differentiation, within‐population genetic diversity and evolutionary relationships among multiple peripheral populations of a cold‐adapted terrestrial salamander, Salamandra atra, at the southern edge of the species core range. We carried out population genetic, phylogeographic, and phylogenetic analyses on various molecular markers (10 autosomal microsatellite loci, three mitochondrial loci with total length >2,100 bp, two protein‐coding nuclear genes) sampled from more than 100 individuals from 13 sites along the southern Prealps. We found at least seven isolated peripheral populations, all highly differentiated from the remaining populations and differentiated from each other at various levels. The within‐population genetic diversity was variable in the peripheral populations, but consistently lower than in the remaining populations. All peripheral populations along the southern Prealps belong to an ancient lineage that is also found in the Dinarides but did not contribute to the postglacial recolonization of the inner and northern Alps. All fully melanistic populations from the Orobian mountains to the southern Dinarides represent a single clade, to the exclusion of the two yellow‐patched populations inhabiting the Pasubio massif and the Sette Comuni plateau, which are distinguished as S. atra pasubiensis and S. atra aurorae, respectively. In conclusion, multiple populations of S. atra at the southern edge of the species core area have different levels of differentiation, different amount of within‐population genetic diversity, and different evolutionary origin. Therefore, they should be regarded as complementary conservation targets to preserve the overall genetic and evolutionary diversity of the species.     

Genetic rescue of an inbred captive population of the critically endangered Puerto Rican crested toad (<Emphasis Type="Italic">Peltophryne lemur</Emphasis>) by mixing lineages

The Puerto Rican crested toad (Peltophryne lemur) is currently composed of a single wild population on the south coast of Puerto Rico and two captive populations founded by animals from the northern and southern coasts. The main factors contributing to its decline are habitat loss, inundation of breeding ponds during storms, and impacts of invasive species. Recovery efforts have been extensive, involving captive breeding and reintroductions, habitat restoration, construction of breeding ponds, and public education. To guide future conservation efforts, genetic variation and differentiation were assessed for the two captive colonies and the remaining wild population using the mitochondrial control region and six novel microsatellite loci. Only two moderately divergent mitochondrial haplotypes were found, with one fixed in each of the southern and northern lineages. Moderate genetic variation exists for microsatellite loci in all three groups. The captive southern population has not diverged substantially from the wild population at microsatellite loci (F _ST = 0.03), whereas there is little allelic overlap between the northern and southern lineages at five of six loci (F _ST > 0.3). Despite this differentiation, they are no more divergent than many populations of other amphibian species. As the northern breeding colony may not remain viable due to its small size and inbred nature, it is recommended that a third breeding colony be established in which northern and southern individuals are combined. This will preserve any northern adaptive traits that may exist, and provide animals for release in the event that the pure northern lineage becomes extirpated.     

Genetic and morphological variation in the venus clam Cyclina?sinensis along the coast of China

Six geographic samples of Cyclina sinensis were collected from the coast along China and analyzed to reveal morphological and genetic variation by using nine allozyme loci and 11 morphological variables. The discriminant function analysis (DFA) of morphology suggested a clear separation between the southern and northern populations. Polymorphism was detected at nine loci across all six populations. The mean allele number ranged from 2.44 to 2.78, and the mean observed heterozygosity ranged from 0.218 to 0.296. High level of genetic differentiation was found between three northern populations and three southern populations. The marked genetic differentiation can be explained by the upwelling of the Zhejiang province and the freshwater outflow of Yangtze River. The results obtained in this study indicate that the northern and southern populations of C. sinensis should be treated as separate units for conservation management.     

Testing hypotheses driving genetic structure in the cooperatively breeding Brown‐headed Nuthatch Sitta pusilla

Habitat fragmentation has often been implicated in the decline of many species. For habitat specialists and/or sedentary species, loss of habitat can result in population isolation and lead to negative genetic effects. However, factors other than fragmentation can often be important and also need to be considered when assessing the genetic structure of a species. We genotyped individuals from 13 populations of the cooperatively breeding Brown‐headed Nuthatch Sitta pusilla in Florida to test three alternative hypotheses regarding the effects that habitat fragmentation might have on genetic structure. A map of potential habitat developed from recent satellite imagery suggested that Brown‐headed Nuthatch populations in southern Florida occupied smaller and more isolated habitat patches (i.e. were more fragmented) than populations in northern Florida. We also genotyped individuals from a small, isolated Brown‐headed Nuthatch population on Grand Bahama Island. We found that populations associated with more fragmented habitat in southern Florida had lower allelic richness than populations in northern Florida (P = 0.02), although there were no differences in heterozygosity. Although pairwise estimates of F_ST were low overall, values among southern populations were generally higher than northern populations. Population assignment tests identified K = 3 clusters corresponding to a northern cluster, a southern cluster and a unique population in southeast Florida; using sampling localities as prior information revealed K = 7 clusters, with greater structure only among southern Florida populations. The Bahamas population showed moderate to high differentiation compared with Florida populations. Overall, our results suggest that fragmentation could affect gene flow in Brown‐headed Nuthatch populations and is likely to become more pronounced over time.     

Genomic signatures of host‐associated divergence and adaptation in a coral‐eating snail,Coralliophila violacea (Kiener, 1836)

The fluid nature of the ocean, combined with planktonic dispersal of marine larvae, lowers physical barriers to gene flow. However, divergence can still occur despite gene flow if strong selection acts on populations occupying different ecological niches. Here, we examined the population genomics of an ectoparasitic snail, Coralliophila violacea (Kiener 1836), that specializes on Porites corals in the Indo‐Pacific. Previous genetic analyses revealed two sympatric lineages associated with different coral hosts. In this study, we examined the mechanisms promoting and maintaining the snails’ adaptation to their coral hosts. Genome‐wide single nucleotide polymorphism (SNP) data from type II restriction site‐associated DNA (2b‐RAD) sequencing revealed two differentiated clusters of C. violacea that were largely concordant with coral host, consistent with previous genetic results. However, the presence of some admixed genotypes indicates gene flow from one lineage to the other. Combined, these results suggest that differentiation between host‐associated lineages of C. violacea is occurring in the face of ongoing gene flow, requiring strong selection. Indeed, 2.7% of all SNP loci were outlier loci (73/2,718), indicative of divergence with gene flow, driven by adaptation of each C. violacea lineage to their specific coral hosts.     

Different genetic structures revealed resident populations of a specialist parasitoid wasp in contrast to its migratory host

Genetic comparisons of parasitoids and their hosts are expected to reflect ecological and evolutionary processes that influence the interactions between species. The parasitoid wasp, Cotesia vestalis, and its host diamondback moth (DBM), Plutella xylostella, provide opportunities to test whether the specialist natural enemy migrates seasonally with its host or occurs as resident population. We genotyped 17 microsatellite loci and two mitochondrial genes for 158 female adults of C. vestalis collected from 12 geographical populations, as well as nine microsatellite loci for 127 DBM larvae from six separate sites. The samplings covered both the likely source (southern) and immigrant (northern) areas of DBM from China. Populations of C. vestalis fell into three groups, pointing to isolation in northwestern and southwestern China and strong genetic differentiation of these populations from others in central and eastern China. In contrast, DBM showed much weaker genetic differentiation and high rates of gene flow. TESS analysis identified the immigrant populations of DBM as showing admixture in northern China. Genetic disconnect between C. vestalis and its host suggests that the parasitoid did not migrate yearly with its host but likely consisted of resident populations in places where its host could not survive in winter.     

Isolation and characterization of 11 polymorphic microsatellite markers for the marine gastropod Concholepas concholepas (Brugière, 1789)

Because of its long‐lived planktonic stage, the marine gastropod Concholepas concholepas is expected to exchange larvae over large distances. However, discrepancies between expected and realized dispersal have been documented in marine invertebrates. To investigate relationships between potential and effective (i.e. gene flow) dispersal, we developed 11 microsatellite markers and investigate their usefulness by analysing two populations distant by c. 4000 km. The 11 loci were found to be highly polymorphic in both populations, with 12–51 alleles according to the locus. This polymorphism is strong enough to allow fine‐scale population analyses including larval studies and paternity analyses.     

SNP panel development for genetic management of wild and domesticated white bass (Morone chrysops)

White bass (Morone chrysops), striped bass and their interspecific hybrid are important game fishes, whereas the hybrid striped bass is an important aquaculture species in the US. Numerous state, federal and private hatcheries, therefore, rear these species for stocking purposes as well as for food fish. Although striped bass populations (both wild and domesticated) have been extensively evaluated, relatively little effort has been directed toward the study and improvement of white bass. In this study, we developed SNP resources to examine the genetic relationships among a long‐term domesticated white bass line and five potential founder stocks for selective breeding collected from drainages in Arkansas, Texas and Alabama. Using genotyping‐by‐sequencing, we generated 13 872 genome‐wide SNP loci across the six populations. Stringent filtering of SNP‐calling parameters identified 426 informative SNP loci. Population genetic and structure analyses using these loci revealed only moderate genetic differentiation between populations (global F_st = 0.083) and indicated two major genetic clusters. A final 57‐SNP assay was successfully designed and validated using the MassARRAY system. The developed SNP panel assigned 96 additional genotyped individuals to their population of origin with 100% accuracy. The SNP resources developed in this study should facilitate ongoing efforts in selective breeding and conservation of white bass.     

High degree of genetic differentiation in marine three‐spined sticklebacks (Gasterosteus aculeatus)

Populations of widespread marine organisms are typically characterized by a low degree of genetic differentiation in neutral genetic markers, but much less is known about differentiation in genes whose functional roles are associated with specific selection regimes. To uncover possible adaptive population divergence and heterogeneous genomic differentiation in marine three‐spined sticklebacks (Gasterosteus aculeatus), we used a candidate gene‐based genome‐scan approach to analyse variability in 138 microsatellite loci located within/close to (<6 kb) functionally important genes in samples collected from ten geographic locations. The degree of genetic differentiation in markers classified as neutral or under balancing selection—as determined with several outlier detection methods—was low (F_ST = 0.033 or 0.011, respectively), whereas average F_ST for directionally selected markers was significantly higher (F_ST = 0.097). Clustering analyses provided support for genomic and geographic heterogeneity in selection: six genetic clusters were identified based on allele frequency differences in the directionally selected loci, whereas four were identified with the neutral loci. Allelic variation in several loci exhibited significant associations with environmental variables, supporting the conjecture that temperature and salinity, but not optic conditions, are important drivers of adaptive divergence among populations. In general, these results suggest that in spite of the high degree of physical connectivity and gene flow as inferred from neutral marker genes, marine stickleback populations are strongly genetically structured in loci associated with functionally relevant genes.     

Geographic structure in the Southern Ocean circumpolar brittle star Ophionotus victoriae (Ophiuridae) revealed from mtDNA and single‐nucleotide polymorphism data

Marine systems have traditionally been thought of as “open” with few barriers to gene flow. In particular, many marine organisms in the Southern Ocean purportedly possess circumpolar distributions that have rarely been well verified. Here, we use the highly abundant and endemic Southern Ocean brittle star Ophionotus victoriae to examine genetic structure and determine whether barriers to gene flow have existed around the Antarctic continent. Ophionotus victoriae possesses feeding planktotrophic larvae with presumed high dispersal capability, but a previous study revealed genetic structure along the Antarctic Peninsula. To test the extent of genetic differentiation within O. victoriae, we sampled from the Ross Sea through the eastern Weddell Sea. Whereas two mitochondrial DNA markers (16S rDNA and COI) were employed to allow comparison to earlier work, a 2b‐RAD single‐nucleotide polymorphism (SNP) approach allowed sampling of loci across the genome. Mitochondrial data from 414 individuals suggested three major lineages, but 2b‐RAD data generated 1,999 biallelic loci that identified four geographically distinct groups from 89 samples. Given the greater resolution by SNP data, O. victoriae can be divided into geographically distinct populations likely representing multiple species. Specific historical scenarios that explain current population structure were examined with approximate Bayesian computation (ABC) analyses. Although the Bransfield Strait region shows high diversity possibly due to mixing, our results suggest that within the recent past, dispersal processes due to strong currents such as the Antarctic Circumpolar Current have not overcome genetic subdivision presumably due to historical isolation, questioning the idea of large open circumpolar populations in the Southern Ocean.     

Outlier SNP markers reveal fine‐scale genetic structuring across European hake populations (Merluccius merluccius)

Shallow population structure is generally reported for most marine fish and explained as a consequence of high dispersal, connectivity and large population size. Targeted gene analyses and more recently genome‐wide studies have challenged such view, suggesting that adaptive divergence might occur even when neutral markers provide genetic homogeneity across populations. Here, 381 SNPs located in transcribed regions were used to assess large‐ and fine‐scale population structure in the European hake (Merluccius merluccius), a widely distributed demersal species of high priority for the European fishery. Analysis of 850 individuals from 19 locations across the entire distribution range showed evidence for several outlier loci, with significantly higher resolving power. While 299 putatively neutral SNPs confirmed the genetic break between basins (F_CT = 0.016) and weak differentiation within basins, outlier loci revealed a dramatic divergence between Atlantic and Mediterranean populations (F_CT range 0.275–0.705) and fine‐scale significant population structure. Outlier loci separated North Sea and Northern Portugal populations from all other Atlantic samples and revealed a strong differentiation among Western, Central and Eastern Mediterranean geographical samples. Significant correlation of allele frequencies at outlier loci with seawater surface temperature and salinity supported the hypothesis that populations might be adapted to local conditions. Such evidence highlights the importance of integrating information from neutral and adaptive evolutionary patterns towards a better assessment of genetic diversity. Accordingly, the generated outlier SNP data could be used for tackling illegal practices in hake fishing and commercialization as well as to develop explicit spatial models for defining management units and stock boundaries.     

Drifting fronds and drifting alleles: range dynamics,local dispersal and habitat isolation shape the population structure of the estuarine seaweed Fucus ceranoides

Aim The seaweed Fucus ceranoides is restricted to spatially discrete estuarine habitats and lacks planktonic dispersal phases; it is therefore expected to exhibit strong population differentiation. Its cold‐temperate affinities and mtDNA variation imply that the northern part of the species’ range, where F. ceranoides is now ubiquitous, was recently colonized after the onset of the last deglaciation, potentially resulting in areas with greater genetic homogeneity. Here we examine the population structure of F. ceranoides to test these predictions, emphasizing the contrasting genetic signatures of limited dispersal in refugial versus recently colonized regions. Location North‐eastern Atlantic estuaries from Portugal to Norway. Methods A total of 504 individuals from 21 estuarine sites spanning the entire range of F. ceranoides were sampled and genotyped for nine microsatellite loci. Population structure was inferred from several genotypic and allele‐frequency analyses. Geographical patterns of genetic diversity were used to reconstruct the historical biogeography of the species. Results Genetic diversity and regional population differentiation showed a consistent decline with increasing latitude. Southernmost populations harboured most of the endemic variation, whereas the northern populations (> 55° N) were almost fixed for the same alleles across loci. In southern and central regions of its distribution, F. ceranoides showed striking population subdivision, with many of the sampled estuaries corresponding to coherent genetic units that were easily discriminated from one another with standard clustering methods. Main conclusions The geographical pattern of genetic diversity supports the long‐term refugial status of Iberia and a post‐glacial range expansion of F. ceranoides into previously glaciated latitudes. Despite the species’ capacity to colonize newly available habitats, the genetic structure of F. ceranoides outside the recently (re)colonized range reveals that gene flow between populations is extremely low. This study provides a remarkable example of how infrequent and spatially limited dispersal can have contrasting effects at the scales of meta‐population (connectivity) versus range dynamics (habitat tracking), and of how dispersal restrictions can result in either genetic divergence or homogeneity depending on the maturity and demographic conditions of the populations.     

No evidence for adaptation to local rhizobial mutualists in the legume Medicago lupulina

Local adaptation is a common but not ubiquitous feature of species interactions, and understanding the circumstances under which it evolves illuminates the factors that influence adaptive population divergence. Antagonistic species interactions dominate the local adaptation literature relative to mutualistic ones, preventing an overall assessment of adaptation within interspecific interactions. Here, we tested whether the legume Medicago lupulina is adapted to the locally abundant species of mutualistic nitrogen‐fixing rhizobial bacteria that vary in frequency across its eastern North American range. We reciprocally inoculated northern and southern M. lupulina genotypes with the northern (Ensifer medicae) or southern bacterium (E. meliloti) in a greenhouse experiment. Despite producing different numbers of root nodules (the structures in which the plants house the bacteria), neither northern nor southern plants produced more seeds, flowered earlier, or were more likely to flower when inoculated with their local rhizobia. We then used a pre‐existing dataset to perform a genome scan for loci that showed elevated differentiation between field‐collected plants that hosted different bacteria. None of the loci we identified belonged to the well‐characterized suite of legume–rhizobia symbiosis genes, suggesting that the rhizobia do not drive genetic divergence between M. lupulina populations. Our results demonstrate that symbiont local adaptation has not evolved in this mutualism despite large‐scale geographic variation in the identity of the interacting species.     

Genome‐wide analyses suggest parallel selection for universal traits may eclipse local environmental selection in a highly mobile carnivore

Ecological and environmental heterogeneity can produce genetic differentiation in highly mobile species. Accordingly, local adaptation may be expected across comparatively short distances in the presence of marked environmental gradients. Within the European continent, wolves (Canis lupus) exhibit distinct north–south population differentiation. We investigated more than 67‐K single nucleotide polymorphism (SNP) loci for signatures of local adaptation in 59 unrelated wolves from four previously identified population clusters (northcentral Europe n = 32, Carpathian Mountains n = 7, Dinaric‐Balkan n = 9, Ukrainian Steppe n = 11). Our analyses combined identification of outlier loci with findings from genome‐wide association study of individual genomic profiles and 12 environmental variables. We identified 353 candidate SNP loci. We examined the SNP position and neighboring megabase (1 Mb, one million bases) regions in the dog (C. lupus familiaris) genome for genes potentially under selection, including homologue genes in other vertebrates. These regions included functional genes for, for example, temperature regulation that may indicate local adaptation and genes controlling for functions universally important for wolves, including olfaction, hearing, vision, and cognitive functions. We also observed strong outliers not associated with any of the investigated variables, which could suggest selective pressures associated with other unmeasured environmental variables and/or demographic factors. These patterns are further supported by the examination of spatial distributions of the SNPs associated with universally important traits, which typically show marked differences in allele frequencies among population clusters. Accordingly, parallel selection for features important to all wolves may eclipse local environmental selection and implies long‐term separation among population clusters.     

Latitudinal divergence in a widespread amphibian: Contrasting patterns of neutral and adaptive genomic variation

Stochastic effects from demographic processes and selection are expected to shape the distribution of genetic variation in spatially heterogeneous environments. As the amount of genetic variation is central for long‐term persistence of populations, understanding how these processes affect variation over large‐scale geographical gradients is pivotal. We investigated the distribution of neutral and putatively adaptive genetic variation, and reconstructed demographic history in the moor frog (Rana arvalis) using 136 individuals from 15 populations along a 1,700‐km latitudinal gradient from northern Germany to northern Sweden. Using double digest restriction‐site associated DNA sequencing we obtained 27,590 single nucleotide polymorphisms (SNPs), and identified differentiation outliers and SNPs associated with growing season length. The populations grouped into a southern and a northern cluster, representing two phylogeographical lineages from different post‐glacial colonization routes. Hybrid index estimation and demographic model selection showed strong support for a southern and northern lineage and evidence of gene flow between regions located on each side of a contact zone. However, patterns of past gene flow over the contact zone differed between neutral and putatively adaptive SNPs. While neutral nucleotide diversity was higher along the southern than the northern part of the gradient, nucleotide diversity in differentiation outliers showed the opposite pattern, suggesting differences in the relative strength of selection and drift along the gradient. Variation associated with growing season length decreased with latitude along the southern part of the gradient, but not along the northern part where variation was lower, suggesting stronger climate‐mediated selection in the north. Outlier SNPs included loci involved in immunity and developmental processes.     

PATTERNS OF REPRODUCTION,GENETIC DIVERSITY,AND GENETIC DIFFERENTIATION IN CALIFORNIA POPULATIONS OF THE GENICULATE CORALLINE ALGA LITHOTHRIX ASPERGILLUM (RHODOPHYTA)1

The reproductive composition and genetic diversity of populations of the red seaweed Lithothrix aspergillum Gray (O. Corallinales) were studied at three southern California sites (Shaw's Cove and Treasure Island, Laguna Beach; Indian Rock, Santa Catalina Island) and at a fourth site (Bodega Bay) located in northern California. Sexually reproducing populations were confined to southern California. Diploid individuals were numerically dominant over haploid (gametophytic) individuals at all sites. Intertidal and subtidal subpopulations from Shaw's Cove differed in their reproductive profiles. Most intertidal specimens found on emersed surfaces were densely branched, turf-forming, and bore tetrasporangial (68.6%), carposporangial (11.4%), or spermatangial (5.7%) conceptacles, reflecting a sexual life history; none produced asexual bispores. In contrast, 74.3% of the larger, loosely branched subtidal specimens bore bisporangial conceptacles indicative of asexual reproduction. Nearly 70% of the Indian Rock thalli showed no evidence of conceptacle formation. Only asexual, diploid bispore-producing thalli were obtained from the Bodega Bay site. Genetic diversity (mean number of alleles per locus, percent of polymorphic loci, and average expected heterozygosity) of diploid L. aspergillum populations varied with life-history characteristics and geographic location. A total of 30 alleles was inferred from zymograms of 16 loci examined by starch-gel electrophoresis; of these loci, 11 were polymorphic. The genetic diversity of sexual, diploid populations of L. aspergillum (alleles per locus [A/L] = 1.4-1.5; percent polymorphic loci [%P] = 37.5-50.0) was relatively high compared with other red seaweeds. Lowest diversity (A/L = 1.0; %P = 0.0) occurred in the exclusively asexual Bodega Bay population which consisted of genetic clones. All sexual L. aspergillum populations deviated significantly from Hardy-Wein-berg expectations due to lower than expected heterozygosity. Genetic differentiation (Wright's F_statistic [F_ST]; Nei's Genetic Distance [D]) among sexually reproducing southern California populations was low (F_ST= 0.030) on a local scale (ca. 5 km), suggesting high levels of gene flow, but high genetic differention (F_ST= 0.390 and 0.406) occurred among southern California populations separated by ca. 70 km. Very high genetic differentiation (F_ST= 0.583–0.683) was obtained between northern and southern California populations separated by 700–760 km. Our genetic and reproductive data suggest that the L. aspergillum population from Bodega Bay is sustained by perennation, vegetative propagation, or asexual reproduction by bispores and may represent an isolated remnant or a population established by a founder event.