Cryptic diversity and spatial genetic variation in the coral Acropora tenuis and its endosymbionts across the Great Barrier Reef.

Genomic studies are uncovering extensive cryptic diversity within reef-building corals, suggesting that evolutionarily and ecologically relevant diversity is highly underestimated in the very organisms that structure coral reefs. Furthermore, endosymbiotic algae within coral host species can confer adaptive responses to environmental stress and may represent additional axes of coral genetic variation that are not constrained by taxonomic divergence of the cnidarian host. Here, we examine genetic variation in a common and widespread, reef-building coral, Acropora tenuis, and its associated endosymbiotic algae along the entire expanse of the Great Barrier Reef (GBR). We use SNPs derived from genome-wide sequencing to characterize the cnidarian coral host and organelles from zooxanthellate endosymbionts (genus Cladocopium). We discover three distinct and sympatric genetic clusters of coral hosts, whose distributions appear associated with latitude and inshore-offshore reef position. Demographic modelling suggests that the divergence history of the three distinct host taxa ranges from 0.5 to 1.5 million years ago, preceding the GBR's formation, and has been characterized by low-to-moderate ongoing inter-taxon gene flow, consistent with occasional hybridization and introgression typifying coral evolution. Despite this differentiation in the cnidarian host, A. tenuis taxa share a common symbiont pool, dominated by the genus Cladocopium (Clade C). Cladocopium plastid diversity is not strongly associated with host identity but varies with reef location relative to shore: inshore colonies contain lower symbiont diversity on average but have greater differences between colonies as compared with symbiont communities from offshore colonies. Spatial genetic patterns of symbiont communities could reflect local selective pressures maintaining coral holobiont differentiation across an inshore-offshore environmental gradient. The strong influence of environment (but not host identity) on symbiont community composition supports the notion that symbiont community composition responds to habitat and may assist in the adaptation of corals to future environmental change.

Decadal erosion of coral assemblages by multiple disturbances in the Palm Islands, central Great Barrier Reef.

Increases in the frequency of perturbations that drive coral community structure, such as severe thermal anomalies and high intensity storms, highlight the need to understand how coral communities recover following multiple disturbances. We describe the dynamics of cover and assemblage composition of corals on exposed inshore reefs in the Palm Islands, central Great Barrier Reef, over 19 years encapsulating major disturbance events such as the severe bleaching event in 1998 and Cyclone Yasi in 2011, along with other minor storm and heat stress events. Over this time, 47.8% of hard coral cover was lost, with a concomitant shift in coral assemblage composition due to taxon-specific rates of mortality during the disturbances, and asymmetric recovery in the aftermath thereof. High recruitment rates of some broadcast-spawning corals, particularly corymbose Acropora spp., even in the absence of adult colonies, indicate that a strong external larval supply replenished the stocks. Conversely, the time required for recovery of slow-growing coral morphologies and life histories was longer than the recurrence times of major disturbances. With interludes between bleaching and cyclones predicted to decrease, the probability of another severe disturbance event before coral cover and assemblage composition approximates historical levels suggests that reefs will continue to erode.

Congruent patterns of connectivity can inform management for broadcast spawning corals on the Great Barrier Reef.

Connectivity underpins the persistence and recovery of marine ecosystems. The Great Barrier Reef (GBR) is the world's largest coral reef ecosystem and managed by an extensive network of no-take zones; however, information about connectivity was not available to optimize the network's configuration. We use multivariate analyses, Bayesian clustering algorithms and assignment tests of the largest population genetic data set for any organism on the GBR to date (Acropora tenuis, >2500 colonies; >50 reefs, genotyped for ten microsatellite loci) to demonstrate highly congruent patterns of connectivity between this common broadcast spawning reef-building coral and its congener Acropora millepora (~950 colonies; 20 reefs, genotyped for 12 microsatellite loci). For both species, there is a genetic divide at around 19°S latitude, most probably reflecting allopatric differentiation during the Pleistocene. GBR reefs north of 19°S are essentially panmictic whereas southern reefs are genetically distinct with higher levels of genetic diversity and population structure, most notably genetic subdivision between inshore and offshore reefs south of 19°S. These broadly congruent patterns of higher genetic diversities found on southern GBR reefs most likely represent the accumulation of alleles via the southward flowing East Australia Current. In addition, signatures of genetic admixture between the Coral Sea and outer-shelf reefs in the northern, central and southern GBR provide evidence of recent gene flow. Our connectivity results are consistent with predictions from recently published larval dispersal models for broadcast spawning corals on the GBR, thereby providing robust connectivity information about the dominant reef-building genus Acropora for coral reef managers.

A population genetic assessment of coral recovery on highly disturbed reefs of the Keppel Island archipelago in the southern Great Barrier Reef.

Coral reefs surrounding the islands lying close to the coast are unique to the Great Barrier Reef (GBR) in that they are frequently exposed to disturbance events including floods caused by cyclonic rainfall, strong winds and occasional periods of prolonged above-average temperatures during summer. In one such group of islands in the southern GBR, the Keppel Island archipelago, climate-driven disturbances frequently result in major coral mortality. Whilst these island reefs have clearly survived such dramatic disturbances in the past, the consequences of extreme mortality events may include the loss of genetic diversity, and hence adaptive potential, and a reduction in fitness due to inbreeding, especially if new recruitment from external sources is limited. Here we examined the level of isolation of the Keppel Island group as well as patterns of gene flow within the Keppel Islands using 10 microsatellite markers in nine populations of the coral, Acropora millepora. Bayesian cluster analysis and assignment tests indicated gene flow is restricted, but not absent, between the outer and inner Keppel Island groups, and that extensive gene flow exists within each of these island groups. Comparison of the Keppel Island data with results from a previous GBR-wide study that included a single Keppel Island population, confirmed that A. millepora in the Keppel Islands is genetically distinct from populations elsewhere on the GBR, with exception of the nearby inshore High Peak Reef just north of the Keppel Islands. We compared patterns of genetic diversity in the Keppel Island populations with those from other GBR populations and found them to be slightly, but significantly lower, consistent with the archipelago being geographically isolated, but there was no evidence for recent bottlenecks or deviation from mutation-drift equilibrium. A high incidence of private alleles in the Keppel Islands, particularly in the outer islands, supports their relative isolation and contributes to the conservation value of the archipelago. The lack of evidence for genetic erosion, in combination with our observation that the North Keppel Island population samples collected in 2002 and 2008, respectively, exhibited a pairwise genetic distance of zero, supports previous published work indicating that, following bleaching, Acropora corals in the Keppel Islands predominantly recover from regrowth of small amounts of remaining live tissue in apparently dead coral colonies. This is likely supplemented by recruitment of larvae from genetically similar, less disturbed populations at nearby reefs, particularly following extreme flood events.

The importance of coral larval recruitment for the recovery of reefs impacted by cyclone Yasi in the central Great Barrier Reef.

Cyclone Yasi, one of the most severe tropical storms on record, crossed the central Great Barrier Reef (GBR) in February 2011, bringing wind speeds of up to 285 km hr⁻¹ and wave heights of at least 10 m, and causing massive destruction to exposed reefs in the Palm Island Group. Following the cyclone, mean (± S.E.) hard coral cover ranged from just 2.1 (0.2) % to 5.3 (0.4) % on exposed reefs and no reproductively mature colonies of any species of Acropora remained. Although no fragments of Acropora were found at impacted exposed sites following the cyclone, small juvenile colonies of Acropora (<10 cm diameter) were present, suggesting that their small size and compact morphologies enabled them to survive the cyclone. By contrast, sheltered reefs appeared to be unaffected by the cyclone. Mean (± S.E.) hard coral cover ranged from 18.2 (2.4) % to 30.0 (1.0) % and a large proportion of colonies of Acropora were reproductively mature. Macroalgae accounted for 8 to 16% of benthic cover at exposed sites impacted by cyclone Yasi but were absent at sheltered sites. Mean (± S.E.) recruitment of acroporids to settlement tiles declined from 25.3 (4.8) recruits tile⁻¹ in the pre-cyclone spawning event (2010) to 15.4 (2.2) recruits tile⁻¹ in the first post-cyclone spawning event (2011). Yet, post-cyclone recruitment did not differ between exposed (15.2±2.1 S.E.) and sheltered sites (15.6±2.2 S.E.), despite the loss of reproductive colonies at the exposed sites, indicating larval input from external sources. Spatial variation in impacts, the survival of small colonies, and larval replenishment to impacted reefs suggest that populations of Acropora have the potential to recover from this severe disturbance, provided that the Palm Islands are not impacted by acute disturbances or suffer additional chronic stressors in the near future.

Sea snakes rarely venture far from home.

The extent to which populations are connected by dispersal influences all aspects of their biology and informs the spatial scale of optimal conservation strategies. Obtaining direct estimates of dispersal is challenging, particularly in marine systems, with studies typically relying on indirect approaches to evaluate connectivity. To overcome this challenge, we combine information from an eight-year mark-recapture study with high-resolution genetic data to demonstrate extremely low dispersal and restricted gene flow at small spatial scales for a large, potentially mobile marine vertebrate, the turtleheaded sea snake (Emydocephalus annulatus). Our mark-recapture study indicated that adjacent bays in New Caledonia (<1.15 km apart) contain virtually separate sea snake populations. Sea snakes could easily swim between bays but rarely do so. Of 817 recaptures of marked snakes, only two snakes had moved between bays. We genotyped 136 snakes for 11 polymorphic microsatellite loci and found statistically significant genetic divergence between the two bays (F(ST)= 0.008, P < 0.01). Bayesian clustering analyses detected low mixed ancestry within bays and genetic relatedness coefficients were higher, on average, within than between bays. Our results indicate that turtleheaded sea snakes rarely venture far from home, which has strong implications for their ecology, evolution, and conservation.

Evaluating fossil calibrations for dating phylogenies in light of rates of molecular evolution: a comparison of three approaches.

Evolutionary and biogeographic studies increasingly rely on calibrated molecular clocks to date key events. Although there has been significant recent progress in development of the techniques used for molecular dating, many issues remain. In particular, controversies abound over the appropriate use and placement of fossils for calibrating molecular clocks. Several methods have been proposed for evaluating candidate fossils; however, few studies have compared the results obtained by different approaches. Moreover, no previous study has incorporated the effects of nucleotide saturation from different data types in the evaluation of candidate fossils. In order to address these issues, we compared three approaches for evaluating fossil calibrations: the single-fossil cross-validation method of Near, Meylan, and Shaffer (2005. Assessing concordance of fossil calibration points in molecular clock studies: an example using turtles. Am. Nat. 165:137-146), the empirical fossil coverage method of Marshall (2008. A simple method for bracketing absolute divergence times on molecular phylogenies using multiple fossil calibration points. Am. Nat. 171:726-742), and the Bayesian multicalibration method of Sanders and Lee (2007. Evaluating molecular clock calibrations using Bayesian analyses with soft and hard bounds. Biol. Lett. 3:275-279) and explicitly incorporate the effects of data type (nuclear vs. mitochondrial DNA) for identifying the most reliable or congruent fossil calibrations. We used advanced (Caenophidian) snakes as a case study; however, our results are applicable to any taxonomic group with multiple candidate fossils, provided appropriate taxon sampling and sufficient molecular sequence data are available. We found that data type strongly influenced which fossil calibrations were identified as outliers, regardless of which method was used. Despite the use of complex partitioned models of sequence evolution and multiple calibrations throughout the tree, saturation severely compressed basal branch lengths obtained from mitochondrial DNA compared with nuclear DNA. The effects of mitochondrial saturation were not ameliorated by analyzing a combined nuclear and mitochondrial data set. Although removing the third codon positions from the mitochondrial coding regions did not ameliorate saturation effects in the single-fossil cross-validations, it did in the Bayesian multicalibration analyses. Saturation significantly influenced the fossils that were selected as most reliable for all three methods evaluated. Our findings highlight the need to critically evaluate the fossils selected by data with different rates of nucleotide substitution and how data with different evolutionary rates affect the results of each method for evaluating fossils. Our empirical evaluation demonstrates that the advantages of using multiple independent fossil calibrations significantly outweigh any disadvantages.

Genetic monandry in 6 viviparous species of true sea snakes.

Using a suite of 10 highly variable microsatellite loci, we conducted genetic paternity analyses for 76 embryos in the broods of 12 pregnant females representing 6 viviparous species of true sea snakes (Hydrophis clade) in the family Elapidae. To our surprise, we uncovered no evidence for multiple paternity within any of the clutches despite the fact that the genetic markers showed high intraspecific heterozygosities and as many as 20 conspecific alleles per locus. This outcome stands in sharp contrast to the rather high (but also variable) frequency of multiple paternity previously reported in many other reptilian species. However, because our study appears to be the first assessment of genetic parentage for any sea snake species (and indeed for any member of the elapid clade), we can only speculate as to whether this apparent monandry by females is a broader phylogenetic characteristic of elapid snakes or whether it might relate somehow to the sea snakes' peculiar lifestyle that uniquely combines viviparity with a marine existence.

Phylogeography of the Mekong mud snake (Enhydris subtaeniata): the biogeographic importance of dynamic river drainages and fluctuating sea levels for semiaquatic taxa in Indochina.

During the Cenozoic, Southeast Asia was profoundly affected by plate tectonic events, dynamic river systems, fluctuating sea levels, shifting coastlines, and climatic variation, which have influenced the ecological and evolutionary trajectories of the Southeast Asian flora and fauna. We examined the role of these paleogeographic factors on shaping phylogeographic patterns focusing on a species of semiaquatic snake, Enhydris subtaeniata (Serpentes: Homalopsidae) using sequence data from three mitochondrial fragments (cytochrome b, ND4, and ATPase-2785 bp). We sampled E. subtaeniata from seven locations in three river drainage basins that encompassed most of this species' range. Genetic diversities were typically low within locations but high across locations. Moreover, each location had a unique suite of haplotypes not shared among locations, and pairwise φ(ST) values (0.713-0.998) were highly significant between all location pairs. Relationships among phylogroups were well resolved and analysis of molecular variance (AMOVA) revealed strong geographical partitioning of genetic variance among the three river drainage basins surveyed. The genetic differences observed among the populations of E. subtaeniata were likely shaped by the Quaternary landscapes of Indochina and the Sunda Shelf. Historically, the middle and lower Mekong consisted of strongly dissected river valleys separated by low mountain ranges and much of the Sunda Shelf consisted of lowland river valleys that served to connect faunas associated with major regional rivers. It is thus likely that the contemporary genetic patterns observed among populations of E. subtaeniata are the result of their histories in a complex terrain that created abundant opportunities for genetic isolation and divergence yet also provided lowland connections across now drowned river valleys.

Contamination level of mercury in red meat products from cetaceans available from South Korea markets.

Levels of total mercury (T-Hg) were surveyed in red meat (n=73) and liver (n=3) from toothed whales, dolphins and porpoises (odontocetes) sold for human consumption in the coastal cities of South Korea. High concentrations of T-Hg were found in the liver products of finless porpoises (18.7 and 156 microg/wet g) and common dolphins (13.2 microg/wet g). The T-Hg concentrations in red meat products were highest in the false killer whale (9.66+/-12.3 microg/wet g, n=9), bottlenose dolphin (10.6+/-12.6 microg/wet g, n=3) and killer whale (13.3 microg/wet g, n=1), and lowest in Cuvier's beaked whale and the harbour porpoise (0.4-0.5 microg/wet g). Thus, most of the products that originated from odontocetes exceeded the safety limit of 0.5 microg/wet g for T-Hg set by the South Korean health authorities for the fishery industry. Pregnant women and other vulnerable sectors of the population living in South Korea should therefore limit their consumption of odontocete products.