Investigating biodiversity patterns across scales in freshwater mussels

Hypotheses relating genomic diversity to community attributes such as abundance and species diversity attract attention from a wide and varied audience because their applications are twofold. First, testing such hypotheses can further a theoretical-and hopefully generalizable-understanding of the forces that assemble communities and create observed patterns of biodiversity. Second, relationships that hold true could ease the burden of data collection for conservation or other urgent applications; for example, a strong correlation between species diversity and genetic diversity could make it possible to use one as a proxy for the other, and focus limited resources on measuring the easier of the two without sacrificing information gained. In a From the Cover article in this issue of Molecular Ecology, Bucholz et al. (2023) explore the relationships between within-species genomic diversity, community relative abundance and community species richness, testing three types of ecological hypotheses in the freshwater mussel communities of the southeastern United States. They find positive relationships between mussel density and species richness, and between genomic diversity within a species and density of that species, but no robust support for the expectation of correlated genomic and species diversity. Their analyses highlight the among-species variability in relationships among these different levels of organization and also the complex ways in which interactions with the broader ecosystem (i.e. unionid mussels require fish hosts for maturation) affect these quantitative relationships, nonetheless pushing forward into the important frontier of community-wide genomic assessment for theoretical and conservation applications.

The Complete Genome Sequences of 2 Species of Geukensia (Brachidontinae, Mytilidae)

We present the complete genome sequences of Geukensia demissa and Geukensia granosissima. Illumina sequencing was performed on genetic material from museum specimens. The reads were assembled using a de novo method followed by a finishing step. The raw and assembled data are available via Genbank.

Evidence of introgression, ecological divergence and adaptation in Asterias sea stars

Hybrid zones are important windows into the evolutionary dynamics of populations, revealing how processes like introgression and adaptation structure population genomic variation. Importantly, they are useful for understanding speciation and how species respond to their environments. Here, we investigate two closely related sea star species, Asterias rubens and A. forbesi, distributed along rocky European and North American coastlines of the North Atlantic, and use genome-wide molecular markers to infer the distribution of genomic variation within and between species in this group. Using genomic data and environmental niche modelling, we document hybridization occurring between northern New England and the southern Canadian Maritimes. We investigate the factors that maintain this hybrid zone, as well as the environmental variables that putatively drive selection within and between species. We find that the two species differ in their environmental niche breadth; Asterias forbesi displays a relatively narrow environmental niche while conversely, A. rubens has a wider niche breadth. Species distribution models accurately predict hybrids to occur within environmental niche overlap, thereby suggesting environmental selection plays an important role in the maintenance of the hybrid zone. Our results imply that the distribution of genomic variation in North Atlantic sea stars is influenced by the environment, which will be crucial to consider as the climate changes.

Twenty-year record of white pox disease in the Florida Keys: importance of environmental risk factors as drivers of coral health

Declining coral populations worldwide place a special premium on identifying risks and drivers that precipitate these declines. Understanding the relationship between disease outbreaks and their drivers can help to anticipate when the risk of a disease pandemic is high. Populations of the iconic branching Caribbean elkhorn coral Acropora palmata have collapsed in recent decades, in part due to white pox disease (WPX). To assess the role that biotic and abiotic factors play in modulating coral disease, we present a predictive model for WPX in A. palmata using 20 yr of disease surveys from the Florida Keys plus environmental information collected simultaneously in situ and via satellite. We found that colony size was the most influential predictor for WPX occurrence, with larger colonies being at higher risk. Water quality parameters of dissolved oxygen saturation, total organic carbon, dissolved inorganic nitrogen, and salinity were implicated in WPX likelihood. Both low and high wind speeds were identified as important environmental drivers of WPX. While high temperature has been identified as an important cause of coral mortality in both bleaching and disease scenarios, our model indicates that the relative influence of HotSpot (positive summertime temperature anomaly) was low and actually inversely related to WPX risk. The predictive model developed here can contribute to enabling targeted strategic management actions and disease surveillance, enabling managers to treat the disease or mitigate disease drivers, thereby suppressing the disease and supporting the persistence of corals in an era of myriad threats.

A Decade of Death and Other Dynamics: Deepening Perspectives on the Diversity and Distribution of Sea Stars and Wasting

AbstractMass mortality events provide valuable insight into biological extremes and also ecological interactions more generally. The sea star wasting epidemic that began in 2013 catalyzed study of the microbiome, genetics, population dynamics, and community ecology of several high-profile species inhabiting the northeastern Pacific but exposed a dearth of information on the diversity, distributions, and impacts of sea star wasting for many lesser-known sea stars and a need for integration across scales. Here, we combine datasets from single-site to coast-wide studies, across time lines from weeks to decades, for 65 species. We evaluated the impacts of abiotic characteristics hypothetically associated with sea star wasting (sea surface temperature, pelagic primary productivity, upwelling wind forcing, wave exposure, freshwater runoff) and species characteristics (depth distribution, developmental mode, diet, habitat, reproductive period). We find that the 2010s sea star wasting outbreak clearly affected a little over a dozen species, primarily intertidal and shallow subtidal taxa, causing instantaneous wasting prevalence rates of 5%-80%. Despite the collapse of some populations within weeks, environmental and species variation protracted the outbreak, which lasted 2-3 years from onset until declining to chronic background rates of ∼2% sea star wasting prevalence. Recruitment began immediately in many species, and in general, sea star assemblages trended toward recovery; however, recovery was heterogeneous, and a marine heatwave in 2019 raised concerns of a second decline. The abiotic stressors most associated with the 2010s sea star wasting outbreak were elevated sea surface temperature and low wave exposure, as well as freshwater discharge in the north. However, detailed data speaking directly to the biological, ecological, and environmental cause(s) and consequences of the sea star wasting outbreak remain limited in scope, unavoidably retrospective, and perhaps always indeterminate. Redressing this shortfall for the future will require a broad spectrum of monitoring studies not less than the taxonomically broad cross-scale framework we have modeled in this synthesis.

Latitudinal variation and plasticity in response to temperature in Geukensia demissa

As global temperatures warm, species must adapt to a changing climate or transition to a different location suitable for their survival. Understanding the extent to which species are able to do so, particularly keystone species, is imperative to ensuring the survival of key ecosystems. The ribbed mussel Geukensia demissa is an integral part of salt marshes along the Atlantic coast of North America. Spatial patterns of genomic and phenotypic divergence have been previously documented, although their link with coastal environmental variation is unknown. Here, we study how populations of G. demissa in the northern (Massachusetts) and southern (Georgia) portions of the species range respond to changes in temperature. We combine assays of variation in oxygen consumption and RNA transcriptomic data with genomic divergence analyses to identify how separate populations of G. demissa may vary in distinct thermal environments. Our results show differences in constitutive oxygen consumption between mussels from Georgia and Massachusetts, as well as shared and disparate patterns of gene expression across temperature profiles. We also find that metabolic genes seem to be a strong component of divergence between these two populations. Our analysis highlights the importance of studying integrative patterns of genomic and phenotypic variation in species that are key for particular ecosystems, and how they might respond to further changes in climate.

Minor Genetic Consequences of a Major Mass Mortality: Short-Term Effects in Pisaster ochraceus

AbstractMass mortality events are increasing globally in frequency and magnitude, largely as a result of human-induced change. The effects of these mass mortality events, in both the long and short term, are of imminent concern because of their ecosystem impacts. Genomic data can be used to reveal some of the population-level changes associated with mass mortality events. Here, we use reduced-representation sequencing to identify potential short-term genetic impacts of a mass mortality event associated with a sea star wasting outbreak. We tested for changes in the population for genetic differentiation, diversity, and effective population size between pre-sea star wasting and post-sea star wasting populations of Pisaster ochraceus-a species that suffered high sea star wasting-associated mortality (75%-100% at 80% of sites). We detected no significant population-based genetic differentiation over the spatial scale sampled; however, the post-sea star wasting population tended toward more differentiation across sites than the pre-sea star wasting population. Genetic estimates of effective population size did not detectably change, consistent with theoretical expectations; however, rare alleles were lost. While we were unable to detect significant population-based genetic differentiation or changes in effective population size over this short time period, the genetic burden of this mass mortality event may be borne by future generations, unless widespread recruitment mitigates the population decline. Prior results from P. ochraceus indicated that natural selection played a role in altering allele frequencies following this mass mortality event. In addition to the role of selection found in a previous study on the genomic impacts of sea star wasting on P. ochraceus, our current study highlights the potential role the stochastic loss of many individuals plays in altering how genetic variation is structured across the landscape. Future genetic monitoring is needed to determine long-term genetic impacts in this long-lived species. Given the increased frequency of mass mortality events, it is important to implement demographic and genetic monitoring strategies that capture baselines and background dynamics to better contextualize species' responses to large perturbations.

Is It in the Stars? Exploring the Relationships between Species' Traits and Sea Star Wasting Disease

AbstractAn explanation for variation in impacts of sea star wasting disease across asteroid species remains elusive. Although various traits have been suggested to play a potential role in sea star wasting susceptibility, currently we lack a thorough comparison that explores how life-history and natural history traits shape responses to mass mortality across diverse asteroid taxa. To explore how asteroid traits may relate to sea star wasting, using available data and recognizing the potential for biological correlations to be driven by phylogeny, we generated a supertree, tested traits for phylogenetic association, and evaluated associations between traits and sea star wasting impact. Our analyses show no evidence for a phylogenetic association with sea star wasting impact, but there does appear to be phylogenetic association for a subset of asteroid life-history traits, including diet, substrate, and reproductive season. We found no relationship between sea star wasting and developmental mode, diet, pelagic larval duration, or substrate but did find a relationship with minimum depth, reproductive season, and rugosity (or surface complexity). Species with the greatest sea star wasting impacts tend to have shallower minimum depth distributions, they tend to have their median reproductive period 1.5 months earlier, and they tend to have higher rugosities relative to species less affected by sea star wasting. Fully understanding sea star wasting remains challenging, in part because dramatic gaps still exist in our understanding of the basic biology and phylogeny of asteroids. Future studies would benefit from a more robust phylogenetic understanding of sea stars, as well as leveraging intra- and interspecific comparative transcriptomics and genomics to elucidate the molecular pathways responding to sea star wasting.

A Review of Asteroid Biology in the Context of Sea Star Wasting: Possible Causes and Consequences

AbstractSea star wasting-marked in a variety of sea star species as varying degrees of skin lesions followed by disintegration-recently caused one of the largest marine die-offs ever recorded on the west coast of North America, killing billions of sea stars. Despite the important ramifications this mortality had for coastal benthic ecosystems, such as increased abundance of prey, little is known about the causes of the disease or the mechanisms of its progression. Although there have been studies indicating a range of causal mechanisms, including viruses and environmental effects, the broad spatial and depth range of affected populations leaves many questions remaining about either infectious or non-infectious mechanisms. Wasting appears to start with degradation of mutable connective tissue in the body wall, leading to disintegration of the epidermis. Here, we briefly review basic sea star biology in the context of sea star wasting and present our current knowledge and hypotheses related to the symptoms, the microbiome, the viruses, and the associated environmental stressors. We also highlight throughout the article knowledge gaps and the data needed to better understand sea star wasting mechanistically, its causes, and potential management.

Gene family amplification facilitates adaptation in freshwater unionid bivalve Megalonaias nervosa

Freshwater unionid bivalves currently face severe anthropogenic challenges. Over 70% of species in the United States are threatened, endangered or extinct due to pollution, damming of waterways and overfishing. These species are notable for their unusual life history strategy, parasite-host co-evolution and biparental mitochondrial inheritance. Among this clade, the washboard mussel Megalonaias nervosa is one species that remains prevalent across the Southeastern United States, with robust population sizes. We have created a reference genome for M. nervosa to determine how genome content has evolved in the face of these widespread environmental challenges. We observe dynamic changes in genome content, with a burst of recent transposable element proliferation causing a 382 Mb expansion in genome content. Birth-death models suggest rapid expansions among gene families, with a mutation rate of 1.16 × 10^-8 duplications per gene per generation. Cytochrome P450 gene families have experienced exceptional recent amplification beyond expectations based on genome-wide birth-death processes. These genes are associated with increased rates of amino acid changes, a signature of selection driving evolution of detox genes. Fitting evolutionary models of adaptation from standing genetic variation, we can compare adaptive potential across species and mutation types. The large population size in M. nervosa suggests a 4.7-fold advantage in the ability to adapt from standing genetic variation compared with a low diversity endemic E. hopetonensis. Estimates suggest that gene family evolution may offer an exceptional substrate of genetic variation in M. nervosa, with P_sgv  = 0.185 compared with P_sgv  = 0.067 for single nucleotide changes. Hence, we suggest that gene family evolution is a source of 'hopeful monsters' within the genome that may facilitate adaptation when selective pressures shift. These results suggest that gene family expansion is a key driver of adaptive evolution in this key species of freshwater Unionidae that is currently facing widespread environmental challenges. This work has clear implications for conservation genomics on freshwater bivalves as well as evolutionary theory. This genome represents a first step to facilitate reverse ecological genomics in Unionidae and identify the genetic underpinnings of phenotypic diversity.

Extended phenotypes on coral reefs: cryptic phenotypes modulate coral-vermetid interactions

Phenotypic variation can lead to variation in the strength and outcome of species interactions. Variation in phenotypic traits can arise due to plastic responses to environmental stimuli, underlying genetic variation, or both, and may reflect differences in the focal organism or aspects of the extended phenotype (e.g., associated microbes). We used a reciprocal transplant experiment of Porites corals to evaluate the role of plasticity vs. heritable diversity on phenotypic traits and performance of corals that varied in their prior exposure to vermetid gastropods, an organism known to reduce coral growth and survival. We measured a suite of phenotypic traits associated with coral performance, many of which showed a plastic response to vermetid exposure. Vermetids decreased calcification of corals, increased microbial diversity, and shifted microbial composition. Most traits also showed a signature of previous exposure environment that persisted even when exposure was reversed: i.e., under the same conditions, corals naïve to vermetids had slower calcification rates, thicker tissues, higher Symbiodiniaceae densities, and different microbiomes than corals previously exposed to vermetids. We suggest the phenotypic differences are heritable, as reefs with and without vermetids were comprised of two different mitotypes, that revealed high, consistent genetic variation. Vermetids were only found on the fast-growing coral mitotype that was characterized by thin tissue, and that likely had a history of disturbance. As extended phenotypes can have community impacts, we suggest vermetid, in addition to microbes, are part of the extended community phenotype of these corals. Coral genotypes can establish different reef trajectories, with thin-tissue types more prone to disturbance and subsequent colonization by other species, like vermetids, which can further facilitate the degradation of coral reefs. The effects of the extended phenotype of species likely influence heterogeneity across landscapes as well as temporal differences in community composition.

Evidence That Microorganisms at the Animal-Water Interface Drive Sea Star Wasting Disease

Sea star wasting (SSW) disease describes a condition affecting asteroids that resulted in significant Northeastern Pacific population decline following a mass mortality event in 2013. The etiology of SSW is unresolved. We hypothesized that SSW is a sequela of microbial organic matter remineralization near respiratory surfaces, one consequence of which may be limited O₂ availability at the animal-water interface. Microbial assemblages inhabiting tissues and at the asteroid-water interface bore signatures of copiotroph proliferation before SSW onset, followed by the appearance of putatively facultative and strictly anaerobic taxa at the time of lesion genesis and as animals died. SSW lesions were induced in Pisaster ochraceus by enrichment with a variety of organic matter (OM) sources. These results together illustrate that depleted O₂ conditions at the animal-water interface may be established by heterotrophic microbial activity in response to organic matter loading. SSW was also induced by modestly (∼39%) depleted O₂ conditions in aquaria, suggesting that small perturbations in dissolved O₂ may exacerbate the condition. SSW susceptibility between species was significantly and positively correlated with surface rugosity, a key determinant of diffusive boundary layer thickness. Tissues of SSW-affected individuals collected in 2013–2014 bore δ¹⁵N signatures reflecting anaerobic processes, which suggests that this phenomenon may have affected asteroids during mass mortality at the time. The impacts of enhanced microbial activity and subsequent O₂ diffusion limitation may be more pronounced under higher temperatures due to lower O₂ solubility, in more rugose asteroid species due to restricted hydrodynamic flow, and in larger specimens due to their lower surface area to volume ratios which affects diffusive respiratory potential.

An initial comparative genomic autopsy of wasting disease in sea stars

Beginning in 2013, sea stars throughout the Eastern North Pacific were decimated by wasting disease, also known as "asteroid idiopathic wasting syndrome" (AIWS) due to its elusive aetiology. The geographic extent and taxonomic scale of AIWS meant events leading up to the outbreak were heterogeneous, multifaceted, and oftentimes unobserved; progression from morbidity to death was rapid, leaving few tell-tale symptoms. Here, we take a forensic genomic approach to discover candidate genes that may help explain sea star wasting syndrome. We report the first genome and annotation for Pisaster ochraceus, along with differential gene expression (DGE) analyses in four size classes, three tissue types, and in symptomatic and asymptomatic individuals. We integrate nucleotide polymorphisms associated with survivors of the wasting disease outbreak, DGE associated with temperature treatments in P. ochraceus, and DGE associated with wasting in another asteroid Pycnopodia helianthoides. In P. ochraceus, we found DGE across all tissues, among size classes, and between asymptomatic and symptomatic individuals; the strongest wasting-associated DGE signal was in pyloric caecum. We also found previously identified outlier loci co-occur with differentially expressed genes. In cross-species comparisons of symptomatic and asymptomatic individuals, consistent responses distinguish genes associated with invertebrate innate immunity and chemical defence, consistent with context-dependent stress responses, defensive apoptosis, and tissue degradation. Our analyses thus highlight genomic constituents that may link suspected environmental drivers (elevated temperature) with intrinsic differences among individuals (age/size, alleles associated with susceptibility) that elicit organismal responses (e.g., coelomocyte proliferation) and manifest as sea star wasting mass mortality.

Maintenance of a Genetic Cline in the Barnacle Balanus glandula

The barnacle Balanus glandula is a broadly distributed species in the temperate northeastern Pacific that is notable for a robust genetic cline between about 36° and 40° N latitude. Prior work established the evolutionary origins of this pattern and proposed that it is maintained by environmental selection. In recent years, "climate velocity" studies in marine habitats have shown dramatic distributional shifts for many species as they track their preferred temperature range in a warming ocean. We re-sampled B. glandula across its entire geographic range to determine whether there has been any shift in this genetic distribution, a development signaling that temperature or other climate factors are maintaining this genetic cline. Additionally, we asked whether the spatially distributed mitochondrial lineages also vary in reproductive output with latitude, using location as a proxy for temperature and other coastal environmental factors. Here we show that although the distribution of the genetic cline has not appreciably changed, there is a notable association of decreased reproductive output at lower latitudes of the distribution in the "northern" lineage of B. glandula.

RNA expression and disease tolerance are associated with a "keystone mutation" in the ochre sea star Pisaster ochraceus

An overdominant mutation in an intron of the elongation factor 1-α (EF1A) gene in the sea star Pisaster ochraceus has shown itself to mediate tolerance to "sea star wasting disease", a pandemic that has significantly reduced sea star populations on the Pacific coast of North America. Here we use RNA sequencing of healthy individuals to identify differences in constitutive expression of gene regions that may help explain this tolerance phenotype. Our results show that individuals carrying this mutation have lower expression at a large contingent of gene regions. Individuals without this mutation also appear to have a greater cellular response to temperature stress, which has been implicated in the outbreak of sea star wasting disease. Given the ecological significance of P. ochraceus, these results may be useful in predicting the evolutionary and demographic future for Pacific intertidal communities.

Parallel Patterns of Host-Specific Morphology and Genetic Admixture in Sister Lineages of a Commensal Barnacle

Symbiotic relationships are often species specific, allowing symbionts to adapt to their host environments. Host generalists, on the other hand, have to cope with diverse environments. One coping strategy is phenotypic plasticity, defined by the presence of host-specific phenotypes in the absence of genetic differentiation. Recent work indicates that such host-specific phenotypic plasticity is present in the West Pacific lineage of the commensal barnacle Chelonibia testudinaria (Linnaeus, 1758). We investigated genetic and morphological host-specific structure in the genetically distinct Atlantic sister lineage of C. testudinaria. We collected adult C. testudinaria from loggerhead sea turtles, horseshoe crabs, and blue crabs along the eastern U.S. coast between Delaware and Florida and in the Gulf of Mexico off Mississippi. We find that shell morphology, especially shell thickness, is host specific and comparable in similar host species between the Atlantic and West Pacific lineages. We did not detect significant genetic differentiation related to host species when analyzing data from 11 nuclear microsatellite loci and mitochondrial sequence data, which is comparable to findings for the Pacific lineage. The most parsimonious explanation for these parallel patterns between distinct lineages of C. testudinaria is that C. testudinaria maintained phenotypic plasticity since the lineages diverged 4-5 mya.

Mitochondrial diversity in Gonionemus (Trachylina:Hydrozoa) and its implications for understanding the origins of clinging jellyfish in the Northwest Atlantic Ocean

Determining whether a population is introduced or native to a region can be challenging due to inadequate taxonomy, the presence of cryptic lineages, and poor historical documentation. For taxa with resting stages that bloom episodically, determining origin can be especially challenging as an environmentally-triggered abrupt appearance of the taxa may be confused with an anthropogenic introduction. Here, we assess diversity in mitochondrial cytochrome oxidase I sequences obtained from multiple Atlantic and Pacific locations, and discuss the implications of our findings for understanding the origin of clinging jellyfish Gonionemus in the Northwest Atlantic. Clinging jellyfish are known for clinging to seagrasses and seaweeds, and have complex life cycles that include resting stages. They are especially notorious as some, although not all, populations are associated with severe sting reactions. The worldwide distribution of Gonionemus has been aptly called a “zoogeographic puzzle” and our results refine rather than resolve the puzzle. We find a relatively deep divergence that may indicate cryptic speciation between Gonionemus from the Northeast Pacific and Northwest Pacific/Northwest Atlantic. Within the Northwest Pacific/Northwest Atlantic clade, we find haplotypes unique to each region. We also find one haplotype that is shared between highly toxic Vladivostok-area populations and some Northwest Atlantic populations. Our results are consistent with multiple scenarios that involve both native and anthropogenic processes. We evaluate each scenario and discuss critical directions for future research, including improving the resolution of population genetic structure, identifying possible lineage admixture, and better characterizing and quantifying the toxicity phenotype.

Large-scale gene flow in the barnacle Jehlius cirratus and contrasts with other broadly-distributed taxa along the Chilean coast

We evaluate the population genetic structure of the intertidal barnacle Jehlius cirratus across a broad portion of its geographic distribution using data from the mitochondrial cytochrome oxidase I (COI) gene region. Despite sampling diversity from over 3,000 km of the linear range of this species, there is only slight regional structure indicated, with overall Φ CT of 0.036 (p < 0.001) yet no support for isolation by distance. While these results suggest greater structure than previous studies of J. cirratus had indicated, the pattern of diversity is still far more subtle than in other similarly-distributed species with similar larval and life history traits. We compare these data and results with recent findings in four other intertidal species that have planktotrophic larvae. There are no clear patterns among these taxa that can be associated with intertidal depth or other known life history traits.

The oceanic concordance of phylogeography and biogeography: a case study in Notochthamalus

Dispersal and adaptation are the two primary mechanisms that set the range distributions for a population or species. As such, understanding how these mechanisms interact in marine organisms in particular - with capacity for long-range dispersal and a poor understanding of what selective environments species are responding to - can provide useful insights for the exploration of biogeographic patterns. Previously, the barnacle Notochthamalus scabrosus has revealed two evolutionarily distinct lineages with a joint distribution that suggests an association with one of the two major biogeographic boundaries (~30°S) along the coast of Chile. However, spatial and genomic sampling of this system has been limited until now. We hypothesized that given the strong oceanographic and environmental shifts associated with the other major biogeographic boundary (~42°S) for Chilean coastal invertebrates, the southern mitochondrial lineage would dominate or go to fixation in locations further to the south. We also evaluated nuclear polymorphism data from 130 single nucleotide polymorphisms to evaluate the concordance of the signal from the nuclear genome with that of the mitochondrial sample. Through the application of standard population genetic approaches along with a Lagrangian ocean connectivity model, we describe the codistribution of these lineages through a simultaneous evaluation of coastal lineage frequencies, an approximation of larval behavior, and current-driven dispersal. Our results show that this pattern could not persist without the two lineages having distinct environmental optima. We suggest that a more thorough integration of larval dynamics, explicit dispersal models, and near-shore environmental analysis can explain much of the coastal biogeography of Chile.

Microsatellite loci discovery from next-generation sequencing data and loci characterization in the epizoic barnacle Chelonibia testudinaria (Linnaeus, 1758)

Microsatellite markers remain an important tool for ecological and evolutionary research, but are unavailable for many non-model organisms. One such organism with rare ecological and evolutionary features is the epizoic barnacle Chelonibia testudinaria (Linnaeus, 1758). Chelonibia testudinaria appears to be a host generalist, and has an unusual sexual system, androdioecy. Genetic studies on host specificity and mating behavior are impeded by the lack of fine-scale, highly variable markers, such as microsatellite markers. In the present study, we discovered thousands of new microsatellite loci from next-generation sequencing data, and characterized 12 loci thoroughly. We conclude that 11 of these loci will be useful markers in future ecological and evolutionary studies on C. testudinaria.

What doesn't kill them makes them stronger: an association between elongation factor 1-α overdominance in the sea star Pisaster ochraceus and "sea star wasting disease"

In recent years, a massive mortality event has killed millions of sea stars, of many different species, along the Pacific coast of North America. This disease event, known as 'sea star wasting disease' (SSWD), is linked to viral infection. In one affected sea star (Pisaster ochraceus), previous work had identified that the elongation factor 1-α locus (EF1A) harbored an intronic insertion allele that is lethal when homozygous yet appears to be maintained at moderate frequency in populations through increased fitness for heterozygotes. The environmental conditions supporting this increased fitness are unknown, but overdominance is often associated with disease. Here, we evaluate populations of P. ochraceus to identify the relationship between SSWD and EF1A genotype. Our data suggest that there may be significantly decreased occurrence of SSWD in individuals that are heterozygous at this locus. These results suggest further studies are warranted to understand the functional relationship between diversity at EF1A and survival in P. ochraceus.

Spatial Homogeneity of Bacterial Communities Associated with the Surface Mucus Layer of the Reef-Building Coral Acropora palmata

Coral surface mucus layer (SML) microbiota are critical components of the coral holobiont and play important roles in nutrient cycling and defense against pathogens. We sequenced 16S rRNA amplicons to examine the structure of the SML microbiome within and between colonies of the threatened Caribbean reef-building coral Acropora palmata in the Florida Keys. Samples were taken from three spatially distinct colony regions—uppermost (high irradiance), underside (low irradiance), and the colony base—representing microhabitats that vary in irradiance and water flow. Phylogenetic diversity (PD) values of coral SML bacteria communities were greater than surrounding seawater and lower than adjacent sediment. Bacterial diversity and community composition was consistent among the three microhabitats. Cyanobacteria, Bacteroidetes, Alphaproteobacteria, and Proteobacteria, respectively were the most abundant phyla represented in the samples. This is the first time spatial variability of the surface mucus layer of A. palmata has been studied. Homogeneity in the microbiome of A. palmata contrasts with SML heterogeneity found in other Caribbean corals. These findings suggest that, during non-stressful conditions, host regulation of SML microbiota may override diverse physiochemical influences induced by the topographical complexity of A. palmata. Documenting the spatial distribution of SML microbes is essential to understanding the functional roles these microorganisms play in coral health and adaptability to environmental perturbations.

Phylogeographic structure and northward range expansion in the barnacle Chthamalus fragilis

The barnacle Chthamalus fragilis is found along the US Atlantic seaboard historically from the Chesapeake Bay southward, and in the Gulf of Mexico. It appeared in New England circa 1900 coincident with warming temperatures, and is now a conspicuous member of rocky intertidal communities extending through the northern shore of Cape Cod, Massachusetts. The origin of northern C. fragilis is debated. It may have spread to New England from the northern end of its historic range through larval transport by ocean currents, possibly mediated by the construction of piers, marinas, and other anthropogenic structures that provided new hard substrate habitat. Alternatively, it may have been introduced by fouling on ships originating farther south in its historic distribution. Here we examine mitochondrial cytochrome c oxidase I sequence diversity and the distribution of mitochondrial haplotypes of C. fragilis from 11 localities ranging from Cape Cod, to Tampa Bay, Florida. We found significant genetic structure between northern and southern populations. Phylogenetic analysis revealed three well-supported reciprocally monophyletic haplogroups, including one haplogroup that is restricted to New England and Virginia populations. While the distances between clades do not suggest cryptic speciation, selection and dispersal barriers may be driving the observed structure. Our data are consistent with an expansion of C. fragilis from the northern end of its mid-19th century range into Massachusetts.

Mitochondrial cytochrome b sequence data are not an improvement for species identification in scleractinian corals

There are well-known difficulties in using the cytochrome oxidase I (COI) mitochondrial gene region for population genetics and DNA barcoding in corals. A recent study of species divergence in the endemic Caribbean genus Agaricia reinforced such knowledge. However, the growing availability of whole mitochondrial genomes may help indicate more promising gene regions for species delineation. I assembled the whole mitochondrial genome for Agaricia fragilis from Illumina single-end 250 bp reads and compared this sequence to that of the congener A. humilis. Although these data suggest that the cytochrome b (CYB) gene region is more promising, comparison of available CYB sequence data from scleractinian and other reef-building corals indicates that multilocus approaches are still probably necessary for phylogenetic and population genetic analysis of recently-diverged coral taxa.

Circulation constrains the evolution of larval development modes and life histories in the coastal ocean

The evolutionary pressures that drive long larval planktonic durations in some coastal marine organisms, while allowing direct development in others, have been vigorously debated. We introduce into the argument the asymmetric dispersal of larvae by coastal currents and find that the strength of the currents helps determine which dispersal strategies are evolutionarily stable. In a spatially and temporally uniform coastal ocean of finite extent, direct development is always evolutionarily stable. For passively drifting larvae, long planktonic durations are stable when the ratio of mean to fluctuating currents is small and the rate at which larvae increase in size in the plankton is greater than the mortality rate (both in units of per time). However, larval behavior that reduces downstream larval dispersal for a given time in plankton will be selected for, consistent with widespread observations of behaviors that reduce dispersal of marine larvae. Larvae with long planktonic durations are shown to be favored not for the additional dispersal they allow, but for the additional fecundity that larval feeding in the plankton enables. We analyzed the spatial distribution of larval life histories in a large database of coastal marine benthic invertebrates and documented a link between ocean circulation and the frequency of planktotrophy in the coastal ocean. The spatial variation in the frequency of species with planktotrophic larvae is largely consistent with our theory; increases in mean currents lead to a decrease in the fraction of species with planktotrophic larvae over a broad range of temperatures.

Selection and demographic history shape the molecular evolution of the gamete compatibility protein bindin in Pisaster sea stars

Reproductive compatibility proteins have been shown to evolve rapidly under positive selection leading to reproductive isolation, despite the potential homogenizing effects of gene flow. This process has been implicated in both primary divergence among conspecific populations and reinforcement during secondary contact; however, these two selective regimes can be difficult to discriminate from each other. Here, we describe the gene that encodes the gamete compatibility protein bindin for three sea star species in the genus Pisaster. First, we compare the full-length bindin-coding sequence among all three species and analyze the evolutionary relationships between the repetitive domains of the variable second bindin exon. The comparison suggests that concerted evolution of repetitive domains has an effect on bindin divergence among species and bindin variation within species. Second, we characterize population variation in the second bindin exon of two species: We show that positive selection acts on bindin variation in Pisaster ochraceus but not in Pisaster brevispinus, which is consistent with higher polyspermy risk in P. ochraceus. Third, we show that there is no significant genetic differentiation among populations and no apparent effect of sympatry with congeners that would suggest selection based on reinforcement. Fourth, we combine bindin and cytochrome c oxidase 1 data in isolation-with-migration models to estimate gene flow parameter values and explore the historical demographic context of our positive selection results. Our findings suggest that positive selection on bindin divergence among P. ochraceus alleles can be accounted for in part by relatively recent northward population expansions that may be coupled with the potential homogenizing effects of concerted evolution.

Genealogical approaches to the temporal origins of the Central American Gap: speciation and divergence in pacific Chthamalus (Sessilia: Chthamalidae)

A large section of the tropical Eastern Pacific coastline is nearly devoid of reef or consolidated habitat, and is known as the Central American Gap as it is associated with a biogeographic transition in fish and invertebrate species. We analyze phylogeographic data for intertidal barnacles (Chthamalus) to identify relevant temporal patterns that describe the origins of this biogeographic transition (the Mexican-Panamic Transition Zone). These contrasts of populations on either side of the transition zone include two pairs of closely related species (C. panamensis and C. hedgecocki; C. southwardorum and a Southern form of C. southwardorum), as well as gene flow data within one species (C. panamensis) that currently is found on both sides of the boundary between provinces. Using sequence data from a prior phylogenetic study, we used traditional (net nucleotide divergence) measures as well as coalescent analyses that incorporate the isolation-migration model to identify the likely time of separation between Northern and Southern taxa in two species pairs. A total of 67 individuals were sequenced at two mitochondrial (cytochrome c oxidase I, 16S) and one nuclear (elongation factor 1-alpha) gene regions. Our analyses indicate that the regional isolation of these intertidal barnacles occurred approximately 315-400kya, with subsequent expansion of C. panamensis from the Southern region into the North much more recently. There are insufficient survey data to conclusively document the absence of species from this group within the Central American Gap region near the Gulf of Tehuantepec. However, appropriate habitat is quite sparse in this region and other environmental factors, including upwelling and water temperature, are likely to be associated with isolation of many species in the Mexican and Panamic provinces sensu stricto. Some taxa may maintain gene flow across this region, but very few genetic studies have been completed on such taxa. Until further work is done, distinguishing between prior hypotheses of a faunal gap, or a faunal transition zone, is somewhat speculative. Additional taxonomic revision will be necessary in Chthamalus but is beyond the scope of this paper.

Duplication and population dynamics shape historic patterns of selection and genetic variation at the major histocompatibility complex in rodents

Genetic variation at the major histocompatibility complex (MHC) is vitally important for wildlife populations to respond to pathogen threats. As natural populations can fluctuate greatly in size, a key issue concerns how population cycles and bottlenecks that could reduce genetic diversity will influence MHC genes. Using 454 sequencing, we characterized genetic diversity at the DRB Class II locus in montane voles (Microtus montanus), a North American rodent that regularly undergoes high-amplitude fluctuations in population size. We tested for evidence of historic balancing selection, recombination, and gene duplication to identify mechanisms maintaining allelic diversity. Counter to our expectations, we found strong evidence of purifying selection acting on the DRB locus in montane voles. We speculate that the interplay between population fluctuations and gene duplication might be responsible for the weak evidence of historic balancing selection and strong evidence of purifying selection detected. To further explore this idea, we conducted a phylogenetically controlled comparative analysis across 16 rodent species with varying demographic histories and MHC duplication events (based on the maximum number of alleles detected per individual). On the basis of phylogenetic generalized linear model-averaging, we found evidence that the estimated number of duplicated loci was positively related to allelic diversity and, surprisingly, to the strength of purifying selection at the DRB locus. Our analyses also revealed that species that had undergone population bottlenecks had lower allelic richness than stable species. This study highlights the need to consider demographic history and genetic structure alongside patterns of natural selection to understand resulting patterns of genetic variation at the MHC.

Genetic diversity of fluorescent proteins in Caribbean agariciid corals

The fluorescent protein (FP) gene family is a highly diverse group of proteins whose expression govern color diversity in corals. Here, we examine the genetic diversity of FPs and the extent to which it can be used to assess phylogenetic relationships within the coral genus Agaricia. Tissue samples were collected throughout the Florida Keys from a wide range of phenotypes within the genus Agaricia (A. agaricites [n = 7], A. fragilis [n = 13], and A. lamarcki [n = 2]), as well as the confamilial species Helioseris cucullata (n = 3). Primers were developed from published cDNA sequences to amplify a region of coding and noncoding sequences of FPs. Cloning reactions were performed to capture the multiple copies of FPs and allele diversity. In the resulting 116 cloned sequences, we identified a 179-bp coding region for phylogenetic analysis. Three distinct clades were found in all 3 species of Agaricia, potentially representing 3 copies of the FP gene. Of the 3 gene copies, 2 contain distinct subclades that display reciprocal monophyly between A. agaricites and A. fragilis, whereas A. lamarcki is polyphyletic. Further resolution of the species phylogeny is necessary to fully understand how genetic diversity within this gene family is distributed among taxa and habitats.

Approximate Bayesian estimation of extinction rate in the Finnish Daphnia magna metapopulation

Approximate Bayesian computation (ABC) is useful for parameterizing complex models in population genetics. In this study, ABC was applied to simultaneously estimate parameter values for a model of metapopulation coalescence and test two alternatives to a strict metapopulation model in the well-studied network of Daphnia magna populations in Finland. The models shared four free parameters: the subpopulation genetic diversity (θS), the rate of gene flow among patches (4Nm), the founding population size (N0) and the metapopulation extinction rate (e) but differed in the distribution of extinction rates across habitat patches in the system. The three models had either a constant extinction rate in all populations (strict metapopulation), one population that was protected from local extinction (i.e. a persistent source), or habitat-specific extinction rates drawn from a distribution with specified mean and variance. Our model selection analysis favoured the model including a persistent source population over the two alternative models. Of the closest 750,000 data sets in Euclidean space, 78% were simulated under the persistent source model (estimated posterior probability = 0.769). This fraction increased to more than 85% when only the closest 150,000 data sets were considered (estimated posterior probability = 0.774). Approximate Bayesian computation was then used to estimate parameter values that might produce the observed set of summary statistics. Our analysis provided posterior distributions for e that included the point estimate obtained from previous data from the Finnish D. magna metapopulation. Our results support the use of ABC and population genetic data for testing the strict metapopulation model and parameterizing complex models of demography.

Extinction hazards in experimental Daphnia magna populations: effects of genotype diversity and environmental variation

Extinction is ubiquitous in natural systems and the ultimate fate of all biological populations. However, the factors that contribute to population extinction are still poorly understood, particularly genetic diversity and composition. A laboratory experiment was conducted to examine the influences of environmental variation and genotype diversity on persistence in experimental Daphnia magna populations. Populations were initiated in two blocks with one, two, three, or six randomly selected and equally represented genotypes, fed and checked for extinction daily, and censused twice weekly over a period of 170 days. Our results show no evidence for an effect of the number of genotypes in a population on extinction hazard. Environmental variation had a strong effect on hazards in both experimental blocks, but the direction of the effect differed between blocks. In the first block, variable environments hastened extinction, while in the second block, hazards were reduced under variable food input. This occurred despite greater fluctuations in population size in variable environments in the second block of our experiment. Our results conflict with previous studies, where environmental variation consistently increased extinction risk. They are also at odds with previous studies in other systems that documented significant effects of genetic diversity on population persistence. We speculate that the lack of sexual reproduction, or the phenotypic similarity among our experimental lines, might underlie the lack of a significant effect of genotype diversity in our study.

Genetic estimates of population age in the water flea, Daphnia magna

Genetic datasets can be used to date evolutionary events, even on recent time scales if sufficient data are available. We used statistics calculated from multilocus microsatellite datasets to estimate population ages in data generated through coalescent simulations and in samples from populations of known age in a metapopulation of Daphnia magna in Finland. Our simulation results show that age estimates improve with additional loci and define a time frame over which these statistics are most useful. On the most recent time scales, assumptions regarding the model of mutation (infinite sites vs. stepwise mutation) have little influence on estimated ages. In older populations, size homoplasy among microsatellite alleles results in a downwards bias for estimates based on the infinite sites model (ISM). In the Finnish D. magna metapopulation, our genetically derived estimated ages were biased upwards. Potential sources of this bias include the underlying model of mutation, gene flow, founder size, and the possibility of persistent source populations in the system. Our simulated data show that genetic age estimation is possible, even for very young populations, but our empirical data highlight the importance of factors such as migration when these statistics are applied in natural populations.

Consequences of a poecilogonous life history for genetic structure in coastal populations of the polychaete Streblospio benedicti

In many species, alternative developmental pathways lead to the production of two distinct phenotypes, promoting the evolution of morphological novelty and diversification. Offspring type in marine invertebrates influences transport time by ocean currents, which dictate dispersal potential and gene flow, and thus has sweeping evolutionary effects on the potential for local adaptation and on rates of speciation, extinction and molecular evolution. Here, we use the polychaete Streblospio benedicti to investigate the effects of dimorphic offspring type on gene flow and genetic structure in coastal populations. We use 84 single nucleotide polymorphism (SNP) markers for this species to assay populations on the East and West Coasts of the United States. Using these markers, we found that in their native East Coast distribution, populations of S. benedicti have high-population genetic structure, but this structure is associated primarily with geographic separation rather than developmental differences. Interestingly, very little genetic differentiation is recovered between individuals of different development types when they occur in the same or nearby populations, further supporting that this is a true case of poecilogony. In addition, we were able to demonstrate that the recently introduced (~100 ya) West Coast populations probably originated from a lecithotrophic population near Delaware.

Mitochondrial lineages in Notochthamalus scabrosus as indicators of coastal recruitment and interactions

A significant genetic cline has previously been identified along the Chilean coast in the barnacle Notochthamalus scabrosus. Samples from the previous study, spanning 800 km, were not able to show whether the southern lineage ultimately goes to fixation at higher latitudes. In addition to expanding the geographic sampling of this species, locations that were sampled approximately four to five generations ago were resampled for this study, enabling a temporal comparison of the location and strength of the observed cline. Here, we show that although the cline persists as expected, the tremendous changes in observed lineage frequencies near the clinal boundary are indicative of source-sink dynamics that may be associated with a codistributed biogeographic transition zone. We also find that the southern lineage does not increase in frequency in more southern populations, suggesting that this lineage is maintained through a combination of density-dependent interactions and a coastal fitness gradient.

Examining an outlier: molecular diversity in the cirripedia

Despite the typical assumption in studies of mitochondrial diversity that such data are useful for approximating population size and demography, studies of sequence diversity in mitochondrial DNA across the Metazoa have shown a surprising excess of rare alleles, a pattern associated either with strong selection or population growth. Previous work has shown that this bias toward an excess of rare alleles is typical across the Crustacea, and in particular, in the Cirripedia (barnacles). Here, we directly evaluate sequence data from studies of barnacle populations to ensure that inclusion of cryptic species is not the cause of this pattern. The results shown here reinforce previous studies that suggest caution in interpreting such patterns of allele frequencies, as they are likely to be influenced both by demographic changes and selection.

Transcriptome analysis and SNP development can resolve population differentiation of Streblospio benedicti, a developmentally dimorphic marine annelid

Next-generation sequencing technology is now frequently being used to develop genomic tools for non-model organisms, which are generally important for advancing studies of evolutionary ecology. One such species, the marine annelid Streblospio benedicti, is an ideal system to study the evolutionary consequences of larval life history mode because the species displays a rare offspring dimorphism termed poecilogony, where females can produce either many small offspring or a few large ones. To further develop S. benedicti as a model system for studies of life history evolution, we apply 454 sequencing to characterize the transcriptome for embryos, larvae, and juveniles of this species, for which no genomic resources are currently available. Here we performed a de novo alignment of 336,715 reads generated by a quarter GS-FLX (Roche 454) run, which produced 7,222 contigs. We developed a novel approach for evaluating the site frequency spectrum across the transcriptome to identify potential signatures of selection. We also developed 84 novel single nucleotide polymorphism (SNP) markers for this species that are used to distinguish coastal populations of S. benedicti. We validated the SNPs by genotyping individuals of different developmental modes using the BeadXPress Golden Gate assay (Illumina). This allowed us to evaluate markers that may be associated with life-history mode.

Population genetics of a trochid gastropod broadens picture of Caribbean Sea connectivity

Background

Regional genetic connectivity models are critical for successful conservation and management of marine species. Even though rocky shore invertebrates have been used as model systems to understand genetic structure in some marine environments, our understanding of connectivity in Caribbean communities is based overwhelmingly on studies of tropical fishes and corals. In this study, we investigate population connectivity and diversity of Cittarium pica, an abundant rocky shore trochid gastropod that is commercially harvested across its natural range, from the Bahamas to Venezuela.

Methodology/Principal Findings

We tested for genetic structure using DNA sequence variation at the mitochondrial COI and 16S loci, AMOVA and distance-based methods. We found substantial differentiation among Caribbean sites. Yet, genetic differentiation was associated only with larger geographic scales within the Caribbean, and the pattern of differentiation only partially matched previous assessments of Caribbean connectivity, including those based on larval dispersal from hydrodynamic models. For instance, the Bahamas, considered an independent region by previous hydrodynamic studies, showed strong association with Eastern Caribbean sites in our study. Further, Bonaire (located in the east and close to the meridional division of the Caribbean basin) seems to be isolated from other Eastern sites.

Conclusions/Significance

The significant genetic structure and observed in C. pica has some commonalities in pattern with more commonly sampled taxa, but presents features, such as the differentiation of Bonaire, that appear unique. Further, the level of differentiation, together with regional patterns of diversity, has important implications for the application of conservation and management strategies in this commercially harvested species.

A "shallow phylogeny" of shallow barnacles (chthamalus)

BACKGROUND

We present a multi-locus phylogenetic analysis of the shallow water (high intertidal) barnacle genus Chthamalus, focusing on member species in the western hemisphere. Understanding the phylogeny of this group improves interpretation of classical ecological work on competition, distributional changes associated with climate change, and the morphological evolution of complex cirripede phenotypes.

METHODOLOGY AND FINDINGS

We use traditional and Bayesian phylogenetic and 'deep coalescent' approaches to identify a phylogeny that supports the monophyly of the mostly American 'fissus group' of Chthamalus, but that also supports a need for taxonomic revision of Chthamalus and Microeuraphia. Two deep phylogeographic breaks were also found within the range of two tropical American taxa (C. angustitergum and C. southwardorum) as well.

CONCLUSIONS

Our data, which include two novel gene regions for phylogenetic analysis of cirripedes, suggest that much more evaluation of the morphological evolutionary history and taxonomy of Chthamalid barnacles is necessary. These data and associated analyses also indicate that the radiation of species in the late Pliocene and Pleistocene was very rapid, and may provide new insights toward speciation via transient allopatry or ecological barriers.

Evolutionary dynamics of transferrin in Notropis

This study reports on the development of comparative data for the transferrin (TF) gene in cyprinid fishes, focusing on the genus Notropis. While previous studies have suggested varied patterns of adaptation influencing the diversity at this gene locus both within and among species, sequence data for (TF) in Notropis exhibit limited evidence for selection. However, there are significant effects of sequence variation associated with Notropis lutipinnis, suggesting some form of diversifying selection acting among populations of this species. Overall, the gene performs well as a second locus for phylogenetic and biogeographic inference and may help improve description of the pattern and the process of diversification in Notropis.

Overdominant maintenance of diversity in the sea star Pisaster ochraceus

When individuals have higher evolutionary fitness because of being heterozygous at a given gene region, it is known as overdominance. Although overdominant selection could represent an important mechanism for maintaining genetic variation within species, the prevalence of this mode of selection appears to be relatively low. Identification of cases of true single-locus heterozygote advantage are thus useful reference points in our overall understanding of how various forms of balancing selection influence and maintain genetic variation in natural populations. Here we report the apparent long-term maintenance of diversity via overdominant selection with homozygous lethality at an elongation factor locus in the sea star Pisaster ochraceus. Observing this pattern in a gene with such major effects on protein assembly indicates that overdominant selection could be a more prevalent factor in maintaining allelic diversity in the wild than previously recognized.

The Apicomplexan whole-genome phylogeny: an analysis of incongruence among gene trees

The protistan phylum Apicomplexa contains many important pathogens and is the subject of intense genome sequencing efforts. Based upon the genome sequences from seven apicomplexan species and a ciliate outgroup, we identified 268 single-copy genes suitable for phylogenetic inference. Both concatenation and consensus approaches inferred the same species tree topology. This topology is consistent with most prior conceptions of apicomplexan evolution based upon ultrastructural and developmental characters, that is, the piroplasm genera Theileria and Babesia form the sister group to the Plasmodium species, the coccidian genera Eimeria and Toxoplasma are monophyletic and are the sister group to the Plasmodium species and piroplasm genera, and Cryptosporidium forms the sister group to the above mentioned with the ciliate Tetrahymena as the outgroup. The level of incongruence among gene trees appears to be high at first glance; only 19% of the genes support the species tree, and a total of 48 different gene-tree topologies are observed. Detailed investigations suggest that the low signal-to-noise ratio in many genes may be the main source of incongruence. The probability of being consistent with the species tree increases as a function of the minimum bootstrap support observed at tree nodes for a given gene tree. Moreover, gene sequences that generate high bootstrap support are robust to the changes in alignment parameters or phylogenetic method used. However, caution should be taken in that some genes can infer a "wrong" tree with strong support because of paralogy, model violations, or other causes. The importance of examining multiple, unlinked genes that possess a strong phylogenetic signal cannot be overstated.

Population divergence in plant species reflects latitudinal biodiversity gradients

The trend for increasing biodiversity from the poles to the tropics is one of the best-known patterns in nature. This latitudinal biodiversity gradient has primarily been documented so far with extant species as the measure of biodiversity. Here, we evaluate the global pattern in biodiversity across latitudes based on the magnitude of genetic population divergence within plant species, using a robust spatial design to compare published allozyme datasets. Like the pattern of plant species richness across latitudes, we expected the divergence among populations of current plant species would have a similar pattern and direction. We found that lower latitudinal populations showed greater genetic differentiation within species after controlling for geographical distance. Our analyses are consistent with previous population-level studies in animals, suggesting a high possibility of tropical peaks in speciation rates associated with observed levels of species richness.

Drift by drift: effective population size is limited by advection

BACKGROUND

Genetic estimates of effective population size often generate surprising results, including dramatically low ratios of effective population size to census size. This is particularly true for many marine species, and this effect has been associated with hypotheses of "sweepstakes" reproduction and selective hitchhiking.

RESULTS

Here we show that in advective environments such as oceans and rivers, the mean asymmetric transport of passively dispersed reproductive propagules will act to limit the effective population size in species with a drifting developmental stage. As advection increases, effective population size becomes decoupled from census size as the persistence of novel genetic lineages is restricted to those that arise in a small upstream portion of the species domain.

CONCLUSION

This result leads to predictions about the maintenance of diversity in advective systems, and complements the "sweepstakes" hypothesis and other hypotheses proposed to explain cases of low allelic diversity in species with high fecundity. We describe the spatial extent of the species domain in which novel allelic diversity will be retained, thus determining how large an appropriately placed marine reserve must be to allow the persistence of endemic allelic diversity.