Speciation-by-depth on coral reefs: Sympatric divergence with gene flow or cryptic transient isolation?

The distributions of many sister species in the sea overlap geographically but are partitioned along depth gradients. The genetic changes leading to depth segregation may evolve in geographic isolation as a prerequisite to coexistence or may emerge during primary divergence leading to new species. These alternatives can now be distinguished via the power endowed by the thousands of scorable loci provided by second-generation sequence data. Here, we revisit the case of two depth-segregated, genetically isolated ecotypes of the nominal Caribbean candelabrum coral Eunicea flexuosa. Previous analyses based on a handful of markers could not distinguish between models of genetic exchange after a period of isolation (consistent with secondary contact) and divergence with gene flow (consistent with primary divergence). Analyses of the history of isolation, genetic exchange and population size based on 15,640 new SNP markers derived from RNAseq data best support models where divergence began 800K BP and include epochs of divergence with gene flow, but with an intermediate period of transient isolation. Results also supported the previous conclusion that recent exchange between the ecotypes occurs asymmetrically from the Shallow lineage to the Deep. Parallel analyses of data from two other corals with depth-segregated populations (Agaricia fragilis and Pocillopora damicornis) suggest divergence leading to depth-segregated populations may begin with a period of symmetric exchange, but that an epoch of population isolation precedes more complete isolation marked by asymmetric introgression. Thus, while divergence-with-gene flow may account for much of the differentiation that separates closely related, depth-segregated species, it remains to be seen whether any critical steps in the speciation process only occur when populations are isolated.

Palaeoclimate ocean conditions shaped the evolution of corals and their skeletons through deep time.

Identifying how past environmental conditions shaped the evolution of corals and their skeletal traits provides a framework for predicting their persistence and that of their non-calcifying relatives under impending global warming and ocean acidification. Here we show that ocean geochemistry, particularly aragonite-calcite seas, drives patterns of morphological evolution in anthozoans (corals, sea anemones) by examining skeletal traits in the context of a robust, time-calibrated phylogeny. The lability of skeletal composition among octocorals suggests a greater ability to adapt to changes in ocean chemistry compared with the homogeneity of the aragonitic skeleton of scleractinian corals. Pulses of diversification in anthozoans follow mass extinctions and reef crises, with sea anemones and proteinaceous corals filling empty niches as tropical reef builders went extinct. Changing environmental conditions will likely diminish aragonitic reef-building scleractinians, but the evolutionary history of the Anthozoa suggests other groups will persist and diversify in their wake.

Host-targeted RAD-Seq reveals genetic changes in the coral Oculina patagonica associated with range expansion along the Spanish Mediterranean coast.

Many organisms are expanding their ranges in response to changing environmental conditions. Understanding the patterns of genetic diversity and adaptation along an expansion front is crucial to assessing a species' long-term success. While next-generation sequencing techniques can reveal these changes in fine detail, ascribing them to a particular species can be difficult for organisms that live in close association with symbionts. Using a novel modified restriction site-associated DNA sequencing (RAD-Seq) protocol to target coral DNA, we collected 595 coral-specific single nucleotide polymorphisms from 189 colonies of the invasive coral Oculina patagonica from the Spanish Mediterranean coast, including established core populations and two expansion fronts. Surprisingly, populations from the recent northern expansion are genetically distinct from the westward expansion and core populations and also harbour greater genetic diversity. We found that temperature may have driven adaptation along the northern expansion, as genome scans for selection found three candidate loci associated with temperature in the north but none in the west. We found no genomic signature of selection associated with artificial substrate, which has been proposed for explaining the rapid spread of O. patagonica. This suggests that this coral is simply an opportunistic colonizer of free space made available by coastal habitat modifications. Our results suggest that unique genetic variation, possibly due to limited dispersal across the Ibiza Channel, an influx of individuals from different depths and/or adaptation to cooler temperatures along the northern expansion front may have facilitated the northward range expansion of O. patagonica in the western Mediterranean.

Switch between Morphospecies of Pocillopora Corals.

Pocillopora corals are the main reef builders in the eastern tropical Pacific. The validity of Pocillopora morphospecies remains under debate because of disagreements between morphological and genetic data. To evaluate the temporal stability of morphospecies in situ, we monitored the shapes of individual colonies in three communities in the southern Gulf of California for 44 months. Twenty-three percent of tagged colonies of Pocillopora damicornis changed to Pocillopora inflata morphology during this time. This switch in identity coincided with a shift to a higher frequency of storms and lower water turbidity (i.e., lower chlorophyll a levels). Seven months after the switch, P. inflata colonies were recovering their original P. damicornis morphology. All colonies of both morphospecies shared a common mitochondrial identity, but most P. damicornis colonies undergoing change were at a site with low-flow conditions. This is the first in situ study to document switching between described morphospecies, and it elucidates the influence of temporal shifts in environmental conditions on morphologically plastic responses.

The invasive coral Oculina patagonica has not been recently introduced to the Mediterranean from the western Atlantic.

BACKGROUND: Effective policies, management, and scientific research programs depend on the correct identification of invasive species as being either native or introduced. However, many species continue to be misidentified. Oculina patagonica, first recorded in the Mediterranean Sea in 1966, is believed to have been introduced in anthropogenic times and expanding in a west to east direction. However, its present identification and status as a recently introduced species remain to be explored. In this study, we used multi-locus genetic data to test whether O. patagonica in the Mediterranean has been recently introduced from the western North Atlantic. RESULTS: We found no genetic or historical demographic evidence to support a recent introduction of O. patagonica from the western North Atlantic or an expansion across the Mediterranean. Instead, Mediterranean and Atlantic populations are genetically distinct and appear to have begun diverging about 5 Mya. We also found evidence of a fossil record of Oculina spp. existing in the eastern North Atlantic millions of years before the present. CONCLUSIONS: Our results suggest that Mediterranean populations of O. patagonica have long been isolated from the western Atlantic, either in undetectable numbers or overlooked and undersampled sites and habitats, and have only recently been expanding to invasive levels as a result of environmental changes. Accurate identification of species' invasive statuses will enable more effective research programs aimed at better understanding the mechanisms promoting the invasive nature of species, which can then lead to the implementation of efficient management plans.

Discordance between morphological and molecular species boundaries among Caribbean species of the reef sponge Callyspongia.

Sponges are among the most species-rich and ecologically important taxa on coral reefs, yet documenting their diversity is difficult due to the simplicity and plasticity of their morphological characters. Genetic attempts to identify species are hampered by the slow rate of mitochondrial sequence evolution characteristic of sponges and some other basal metazoans. Here we determine species boundaries of the Caribbean coral reef sponge genus Callyspongia using a multilocus, model-based approach. Based on sequence data from one mitochondrial (COI), one ribosomal (28S), and two single-copy nuclear protein-coding genes, we found evolutionarily distinct lineages were not concordant with current species designations in Callyspongia. While C. fallax,C. tenerrima, and C. plicifera were reciprocally monophyletic, four taxa with different morphologies (C. armigera,C. longissima,C. eschrichtii, and C. vaginalis) formed a monophyletic group and genetic distances among these taxa overlapped distances within them. A model-based method of species delimitation supported collapsing these four into a single evolutionary lineage. Variation in spicule size among these four taxa was partitioned geographically, not by current species designations, indicating that in Callyspongia, these key taxonomic characters are poor indicators of genetic differentiation. Taken together, our results suggest a complex relationship between morphology and species boundaries in sponges.

Morphological variation and different branch modularity across contrasting flow conditions in dominant Pocillopora reef-building corals.

Pocillopora corals, the dominant reef-builders in the Eastern Tropical Pacific, exhibit a high level of phenotypic plasticity, making the interpretation of morphological variation and the identification of species challenging. To test the hypothesis that different coral morphospecies represent phenotypes that develop in different flow conditions, we compared branch characters in three Pocillopora morphospecies (P. damicornis, P. verrucosa, and P. meandrina) from two communities in the Gulf of California exposed to contrasting flow conditions. Morphological variation and branch modularity (i.e., the tendency of different sets of branch traits to vary in a coordinated way) were assessed in colonies classified as Pocillopora type 1 according to two mitochondrial regions. Our results can be summarized as follows. (1) Pocillopora type 1 morphospecies corresponded to a pattern of morphological variation in the Gulf of California. Overall, P. damicornis had the thinnest branches and its colonies the highest branch density, followed by P. verrucosa, and then by P. meandrina, which had the thickest branches and its colonies the lowest branch density. (2) The differentiation among morphospecies was promoted by different levels of modularity of traits. P. verrucosa had the highest coordination of traits, followed by P. damicornis, and P. meandrina. (3) The variation and modularity of branch traits were related to water flow condition. Morphology under the high-flow condition was more similar among morphospecies than under the low-flow condition and seemed to be related to mechanisms for coping with these conditions. Our results provide the first evidence that in scleractinian corals different levels of modularity can be promoted by different environmental conditions.

Historical and recent processes shaping the geographic range of a rocky intertidal gastropod: phylogeography, ecology, and habitat availability.

Factors shaping the geographic range of a species can be identified when phylogeographic patterns are combined with data on contemporary and historical geographic distribution, range-wide abundance, habitat/food availability, and through comparisons with codistributed taxa. Here, we evaluate range dynamism and phylogeography of the rocky intertidal gastropod Mexacanthina lugubris lugubris across its geographic range - the Pacific coast of the Baja peninsula and southern California. We sequenced mitochondrial DNA (CO1) from ten populations and compliment these data with museum records, habitat availability and range-wide field surveys of the distribution and abundance of M. l. lugubris and its primary prey (the barnacle Chthamalus fissus). The geographic range of M. l. lugubris can be characterized by three different events in its history: an old sundering in the mid-peninsular region of Baja (∼ 417,000 years ago) and more recent northern range expansion and southern range contraction. The mid-peninsular break is shared with many terrestrial and marine species, although M. l. lugubris represents the first mollusc to show it. This common break is often attributed to a hypothesized ancient seaway bisecting the peninsula, but for M. l. lugubris it may result from large habitat gaps in the southern clade. Northern clade populations, particularly near the historical northern limit (prior to the 1970s), have high local abundances and reside in a region with plentiful food and habitat - which makes its northern range conducive to expansion. The observed southern range contraction may result from the opposite scenario, with little food or habitat nearby. Our study highlights the importance of taking an integrative approach to understanding the processes that shape the geographic range of a species via combining range-wide phylogeography data with temporal geographic distributions and spatial patterns of habitat/food availability.

Strong natural selection on juveniles maintains a narrow adult hybrid zone in a broadcast spawner.

Natural selection can maintain and help form species across different habitats, even when dispersal is high. Selection against inferior migrants (immigrant inviability) acts when locally adapted populations suffer high mortality on dispersal to unsuitable habitats. Habitat-specific populations undergoing divergent selection via immigrant inviability should thus show (1) a change in the ratio of adapted to nonadapted individuals among age/size classes and (2) a cline (defined by the environmental gradient) as selection counterbalances migration. Here we examine the frequencies of two depth-segregated lineages in juveniles and adults of a Caribbean octocoral, Eunicea flexuosa. Distributions of the two lineages in both shallow and deep environments were more distinct when inferred from adults than juveniles. Despite broad larval dispersal, we also found an extremely narrow hybrid zone (<100 m), with coincident clines for molecular and morphological characters of the host coral and its algal symbiont. Effective dispersal estimates derived from the hybrid zone are remarkably small (<20 m) for a broadcast spawner. The large selection coefficient against mismatched genotypes derived from cohort data is consistent with that from field transplant experiments. Narrow hybrid zones and limited effective dispersal may be a common outcome of long periods of postsettlement, prereproductive selection across steep ecological gradients. Strong diversifying selection provides a mechanism to explain the prevalence of depth-segregated sibling species in the sea.

Cryptic diversity hides host and habitat specialization in a gorgonian-algal symbiosis.

Shallow water anthozoans, the major builders of modern coral reefs, enhance their metabolic and calcification rates with algal symbionts. Controversy exists over whether these anthozoan-algae associations are flexible over the lifetimes of individual hosts, promoting acclimative plasticity, or are closely linked, such that hosts and symbionts co-evolve across generations. Given the diversity of algal symbionts and the morphological plasticity of many host species, cryptic variation within either partner could potentially confound studies of anthozoan-algal associations. Here, we used ribosomal, organelle and nuclear sequences, along with microsatellite variation, to study the relationship between lineages of a common Caribbean gorgonian and its algal symbionts. The gorgonian Eunicea flexuosa is a broadcast spawner, composed of two recently diverged, genetically distinct lineages largely segregated by depth. We sampled colonies of the two lineages across depth gradients at three Caribbean locations. We find that each host lineage is associated with a unique Symbiodinium B1/184 phylotype. This relationship between host and symbiont is maintained when host colonies are reciprocally transplanted, although cases of within phylotype switching were also observed. Even when the phylotypes of both partners are present at intermediate depths, the specificity between host and symbiont lineages remained absolute. Unrecognized cryptic diversity may mask host-symbiont specificity and change the inference of evolutionary processes in mutualistic associations. Symbiotic specificity thus likely contributes to the ecological divergence of the two partners, generating species diversity within coral reefs.

Unrecognized coral species diversity masks differences in functional ecology.

Porites corals are foundation species on Pacific reefs but a confused taxonomy hinders understanding of their ecosystem function and responses to climate change. Here, we show that what has been considered a single species in the eastern tropical Pacific, Porites lobata, includes a morphologically similar yet ecologically distinct species, Porites evermanni. While P. lobata reproduces mainly sexually, P. evermanni dominates in areas where triggerfish prey on bioeroding mussels living within the coral skeleton, thereby generating asexual coral fragments. These fragments proliferate in marginal habitat not colonized by P. lobata. The two Porites species also show a differential bleaching response despite hosting the same dominant symbiont subclade. Thus, hidden diversity within these reef-builders has until now obscured differences in trophic interactions, reproductive dynamics and bleaching susceptibility, indicative of differential responses when confronted with future climate change.

Evaluating summary statistics used to test for incomplete lineage sorting: mito-nuclear discordance in the reef sponge Callyspongia vaginalis.

Conflicting patterns of population differentiation between the mitochondrial and nuclear genomes (mito-nuclear discordance) have become increasingly evident as multilocus data sets have become easier to generate. Incomplete lineage sorting (ILS) of nucDNA is often implicated as the cause of such discordance, stemming from the large effective population size of nucDNA relative to mtDNA. However, selection, sex-biased dispersal and historical demography can also lead to mito-nuclear discordance. Here, we compare patterns of genetic diversity and subdivision for six nuclear protein-coding gene regions to those for mtDNA in a common Caribbean coral reef sponge, Callyspongia vaginalis, along the Florida reef tract. We also evaluated a suite of summary statistics to determine which are effective metrics for comparing empirical and simulated data when testing drivers of mito-nuclear discordance in a statistical framework. While earlier work revealed three divergent and geographically subdivided mtDNACOI haplotypes separated by 2.4% sequence divergence, nuclear alleles were admixed with respect to mitochondrial clade and geography. Bayesian analysis showed that substitution rates for the nuclear loci were up to 7 times faster than for mitochondrial COI. Coalescent simulations and neutrality tests suggested that mito-nuclear discordance in C. vaginalis is not the result of ILS in the nucDNA or selection on the mtDNA but is more likely caused by changes in population size. Sperm-mediated gene flow may also influence patterns of population subdivision in the nucDNA.

Long prereproductive selection and divergence by depth in a Caribbean candelabrum coral.

Long-lived corals, the foundation of modern reefs, often follow ecological gradients, so that populations or sister species segregate by habitat. Adaptive divergence maintains sympatric congeners after secondary contact or may even generate species by natural selection in the face of gene flow. Such ecological divergence, initially between alternative phenotypes within populations, may be aided by immigrant inviability, especially when a long period separates larval dispersal and the onset of reproduction, during which selection can sort lineages to match different habitats. Here, we evaluate the strength of one ecological factor (depth) to isolate populations by comparing the genes and morphologies of pairs of depth-segregated populations of the candelabrum coral Eunicea flexuosa across the Caribbean. Eunicea is endemic to the Caribbean and all sister species co-occur. Eunicea flexuosa is widespread both geographically and across reef habitats. Our genetic analysis revealed two depth-segregated lineages. Field survivorship data, combined with estimates of selection coefficients based on transplant experiments, suggest that selection is strong enough to segregate these two lineages. Genetic exchange between the Shallow and Deep lineages occurred either immediately after divergence or the two have diverged with gene flow. Migration occurs asymmetrically from the Shallow to Deep lineage. Limited recruitment to reproductive age, even under weak annual selection advantage, is sufficient to generate habitat segregation because of the cumulative prolonged prereproductive selection. Ecological factors associated with depth can act as filters generating strong barriers to gene flow, altering morphologies, and contributing to the potential for speciation in the sea.

No gene flow across the Eastern Pacific Barrier in the reef-building coral Porites lobata.

The expanse of deep water between the central Pacific islands and the continental shelf of the Eastern Tropical Pacific is regarded as the world's most potent marine biogeographic barrier. During recurrent climatic fluctuations (ENSO, El Niño Southern Oscillation), however, changes in water temperature and the speed and direction of currents become favourable for trans-oceanic dispersal of larvae from central Pacific to marginal eastern Pacific reefs. Here, we investigate the population connectivity of the reef-building coral Porites lobata across the Eastern Pacific Barrier (EPB). Patterns of recent gene flow in samples (n = 1173) from the central Pacific and the Eastern Tropical Pacific (ETP) were analysed with 12 microsatellite loci. Results indicated that P. lobata from the ETP are strongly isolated from those in the central Pacific and Hawaii (F(ct) ' = 0.509; P < 0.001). However, samples from Clipperton Atoll, an oceanic island on the eastern side of the EPB, grouped with the central Pacific. Within the central Pacific, Hawaiian populations were strongly isolated from three co-occurring clusters found throughout the remainder of the central Pacific. No further substructure was evident in the ETP. Changes in oceanographic conditions during ENSO over the past several thousand years thus appear insufficient to support larval deliveries from the central Pacific to the ETP or strong postsettlement selection acts on ETP settlers from the central Pacific. Recovery of P. lobata populations in the frequently disturbed ETP thus must depend on local larval sources.

The Tegula tango: a coevolutionary dance of interacting, positively selected sperm and egg proteins.

Reproductive proteins commonly show signs of rapid divergence driven by positive selection. The mechanisms driving these changes have remained ambiguous in part because interacting male and female proteins have rarely been examined. We isolate an egg protein the vitelline envelope receptor for lysin (VERL) from Tegula, a genus of free-spawning marine snails. Like VERL from abalone, Tegula VERL is a major component of the VE surrounding the egg, includes a conserved zona pellucida (ZP) domain at its C-terminus, and possesses a unique, negatively charged domain of about 150 amino acids implicated in interactions with the positively charged lysin. Unlike for abalone VERL, where this unique VERL domain occurs in a tandem array of 22 repeats, Tegula VERL has just one such domain. Interspecific comparisons show that both lysin and the VERL domain diverge via positive selection, whereas the ZP domain evolves neutrally. Rates of nonsynonymous substitution are correlated between lysin and the VERL domain, consistent with sexual antagonism, although lineage-specific effects, perhaps owing to different ecologies, may alter the relative evolutionary rates of sperm- and egg-borne proteins.

Reconciling molecules and morphology: molecular systematics and biogeography of Neotropical blennies (Acanthemblemaria).

Neotropical reef fish communities are species-poor compared to those of the Indo-West Pacific. An exception to that pattern is the blenny clade Chaenopsidae, one of only three rocky and coral reef fish families largely endemic to the Neotropics. Within the chaenopsids, the genus Acanthemblemaria is the most species-rich and is characterized by elaborate spinous processes on the skull. Here we construct a species tree using five nuclear markers and compare the results to those from Bayesian and parsimony phylogenetic analyses of 60 morphological characters. The sequence-based species tree conflicted with the morphological phylogenies for Acanthemblemaria, primarily due to the convergence of a suite of characters describing the distribution of spines on the head. However, we were able to resolve some of these conflicts by performing phylogenetic analyses on suites of characters not associated with head spines. By using the species tree as a guide, we used a quantitative method to identify suites of correlated morphological characters that, together, produce the distinctive skull phenotypes found in these fishes. A time calibrated phylogeny with nearly complete taxon sampling provided divergence time estimates that recovered a mid-Miocene origin for the genus, with a temporally and geographically complex pattern of speciation both before and after the closure of the Isthmus of Panama. Some sister taxa are broadly sympatric, but many occur in allopatry. The ability to infer the geography of speciation in Acanthemblemaria is complicated by extinctions, incomplete knowledge of their present geographic ranges and by wide-spread taxa that likely represent cryptic species complexes.

A new lysozyme from the eastern oyster, Crassostrea virginica, and a possible evolutionary pathway for i-type lysozymes in bivalves from host defense to digestion.

BACKGROUND: Lysozymes are enzymes that lyse bacterial cell walls, an activity widely used for host defense but also modified in some instances for digestion. The biochemical and evolutionary changes between these different functional forms has been well-studied in the c-type lysozymes of vertebrates, but less so in the i-type lysozymes prevalent in most invertebrate animals. Some bivalve molluscs possess both defensive and digestive lysozymes. RESULTS: We report a third lysozyme from the oyster Crassostrea virginica, cv-lysozyme 3. The chemical properties of cv-lysozyme 3 (including molecular weight, isoelectric point, basic amino acid residue number, and predicted protease cutting sites) suggest it represents a transitional form between lysozymes used for digestion and immunity. The cv-lysozyme 3 protein inhibited the growth of bacteria (consistent with a defensive function), but semi-quantitative RT-PCR suggested the gene was expressed mainly in digestive glands. Purified cv-lysozyme 3 expressed maximum muramidase activity within a range of pH (7.0 and 8.0) and ionic strength (I = 0.005-0.01) unfavorable for either cv-lysozyme 1 or cv-lysozyme 2 activities. The topology of a phylogenetic analysis of cv-lysozyme 3 cDNA (full length 663 bp, encoding an open reading frame of 187 amino acids) is also consistent with a transitional condition, as cv-lysozyme 3 falls at the base of a monophyletic clade of bivalve lysozymes identified from digestive glands. Rates of nonsynonymous substitution are significantly high at the base of this clade, consistent with an episode of positive selection associated with the functional transition from defense to digestion. CONCLUSION: The pattern of molecular evolution accompanying the shift from defensive to digestive function in the i-type lysozymes of bivalves parallels those seen for c-type lysozymes in mammals and suggests that the lysozyme paralogs that enhance the range of physiological conditions for lysozyme activity may provide stepping stones between defensive and digestive forms.

Ecological partitioning among parapatric cryptic species.

Geographic range differences among species may result from differences in their physiological tolerances. In the intertidal zone, marine and terrestrial environments intersect to create a unique habitat, across which physiological tolerance strongly influences range. Traits to cope with environmental extremes are particularly important here because many species live near their physiological limits and environmental gradients can be steep. The snail Melampus bidentatus occurs in coastal salt marshes in the western Atlantic and the Gulf of Mexico. We used sequence data from one mitochondrial (COI) and two nuclear markers (histone H3 and a mitochondrial carrier protein, MCP) to identify three cryptic species within this broad-ranging nominal species, two of which have partially overlapping geographic ranges. High genetic diversity, low population structure, and high levels of migration within these two overlapping species suggest that historical range limitations do not entirely explain their different ranges. To identify microhabitat differences between these two species, we modelled their distributions using data from both marine and terrestrial environments. Although temperature was the largest factor setting range limits, other environmental components explained features of the ranges that temperature alone could not. In particular, the interaction of precipitation and salinity likely sets physiological limits that lead to range differences between these two cryptic species. This suggests that the response to climatic change in these snails will be mediated by changes to multiple environmental factors, and not just to temperature alone.

Nuclear and mitochondrial sequence data reveal and conceal different demographic histories and population genetic processes in Caribbean reef fishes.

Mitochondrial and nuclear sequence data should recover historical demographic events at different temporal scales due to differences in their effective population sizes and substitution rates. This expectation was tested for two closely related coral reef fish, the tube blennies Acanthemblemaria aspera and A. spinosa. These two have similar life histories and dispersal potentials, and co-occur throughout the Caribbean. Sequence data for one mitochondrial and two nuclear markers were collected for 168 individuals across the species' Caribbean ranges. Although both species shared a similar pattern of genetic subdivision, A. spinosa had 20-25 times greater nucleotide sequence divergence among populations than A. aspera at all three markers. Substitution rates estimated using a relaxed clock approach revealed that mitochondrial COI is evolving at 11.2% pairwise sequence divergence per million years. This rapid mitochondrial rate had obscured the signal of old population expansions for both species, which were only recovered using the more slowly evolving nuclear markers. However, the rapid COI rate allowed the recovery of a recent expansion in A. aspera corresponding to a period of increased habitat availability. Only by combining both nuclear and mitochondrial data were we able to recover the complex demographic history of these fish.

High amino acid diversity and positive selection at a putative coral immunity gene (tachylectin-2).

BACKGROUND: Genes involved in immune functions, including pathogen recognition and the activation of innate defense pathways, are among the most genetically variable known, and the proteins that they encode are often characterized by high rates of amino acid substitutions, a hallmark of positive selection. The high levels of variation characteristic of immunity genes make them useful tools for conservation genetics. To date, highly variable immunity genes have yet to be found in corals, keystone organisms of the world's most diverse marine ecosystem, the coral reef. Here, we examine variation in and selection on a putative innate immunity gene from Oculina, a coral genus previously used as a model for studies of coral disease and bleaching. RESULTS: In a survey of 244 Oculina alleles, we find high nonsynonymous variation and a signature of positive selection, consistent with a putative role in immunity. Using computational protein structure prediction, we generate a structural model of the Oculina protein that closely matches the known structure of tachylectin-2 from the Japanese horseshoe crab (Tachypleus tridentatus), a protein with demonstrated function in microbial recognition and agglutination. We also demonstrate that at least three other genera of anthozoan cnidarians (Acropora, Montastrea and Nematostella) possess proteins structurally similar to tachylectin-2. CONCLUSIONS: Taken together, the evidence of high amino acid diversity, positive selection and structural correspondence to the horseshoe crab tachylectin-2 suggests that this protein is 1) part of Oculina's innate immunity repertoire, and 2) evolving adaptively, possibly under selective pressure from coral-associated microorganisms. Tachylectin-2 may serve as a candidate locus to screen coral populations for their capacity to respond adaptively to future environmental change.