Influence of offshore oil and gas structures on seascape ecological connectivity.

Offshore platforms, subsea pipelines, wells and related fixed structures supporting the oil and gas (O&G) industry are prevalent in oceans across the globe, with many approaching the end of their operational life and requiring decommissioning. Although structures can possess high ecological diversity and productivity, information on how they interact with broader ecological processes remains unclear. Here, we review the current state of knowledge on the role of O&G infrastructure in maintaining, altering or enhancing ecological connectivity with natural marine habitats. There is a paucity of studies on the subject with only 33 papers specifically targeting connectivity and O&G structures, although other studies provide important related information. Evidence for O&G structures facilitating vertical and horizontal seascape connectivity exists for larvae and mobile adult invertebrates, fish and megafauna; including threatened and commercially important species. The degree to which these structures represent a beneficial or detrimental net impact remains unclear, is complex and ultimately needs more research to determine the extent to which natural connectivity networks are conserved, enhanced or disrupted. We discuss the potential impacts of different decommissioning approaches on seascape connectivity and identify, through expert elicitation, critical knowledge gaps that, if addressed, may further inform decision making for the life cycle of O&G infrastructure, with relevance for other industries (e.g. renewables). The most highly ranked critical knowledge gap was a need to understand how O&G structures modify and influence the movement patterns of mobile species and dispersal stages of sessile marine species. Understanding how different decommissioning options affect species survival and movement was also highly ranked, as was understanding the extent to which O&G structures contribute to extending species distributions by providing rest stops, foraging habitat, and stepping stones. These questions could be addressed with further dedicated studies of animal movement in relation to structures using telemetry, molecular techniques and movement models. Our review and these priority questions provide a roadmap for advancing research needed to support evidence-based decision making for decommissioning O&G infrastructure.

Marine environmental DNA biomonitoring reveals seasonal patterns in biodiversity and identifies ecosystem responses to anomalous climatic events.

Marine ecosystems are changing rapidly as the oceans warm and become more acidic. The physical factors and the changes to ocean chemistry that they drive can all be measured with great precision. Changes in the biological composition of communities in different ocean regions are far more challenging to measure because most biological monitoring methods focus on a limited taxonomic or size range. Environmental DNA (eDNA) analysis has the potential to solve this problem in biological oceanography, as it is capable of identifying a huge phylogenetic range of organisms to species level. Here we develop and apply a novel multi-gene molecular toolkit to eDNA isolated from bulk plankton samples collected over a five-year period from a single site. This temporal scale and level of detail is unprecedented in eDNA studies. We identified consistent seasonal assemblages of zooplankton species, which demonstrates the ability of our toolkit to audit community composition. We were also able to detect clear departures from the regular seasonal patterns that occurred during an extreme marine heatwave. The integration of eDNA analyses with existing biotic and abiotic surveys delivers a powerful new long-term approach to monitoring the health of our world's oceans in the context of a rapidly changing climate.

Setting priorities for conservation at the interface between ocean circulation, connectivity, and population dynamics.

Population persistence in the marine environment is driven by patterns of ocean circulation, larval dispersal, ecological interactions, and demographic rates. For habitat-forming organisms in particular, understanding the relationship between larval connectivity and meta-population dynamics aids in planning for marine spatial management. Here, we estimate networks of connectivity between fringing coral reefs in the northwest shelf of Australia by combining a particle tracking model based on shelf circulation with models of subpopulation dynamics of individual reefs. Coral cover data were used as a proxy for overall habitat quality, which can change as a result of natural processes, human-driven impacts, and management initiatives. We obtain three major results of conservation significance. First, the dynamics of the ecological network result from the interplay between network connectivity and ecological processes on individual reefs. The maximum coral cover a zone can sustain imposes a significant nonlinearity on the role an individual reef plays within the dynamics of the network, and thus on the impact of conservation interventions on specific reefs. Second, the role of an individual reef within these network dynamics changes considerably depending on the overall state of the system: a reef's role in sustaining the system's state can be different from the same reef's role in helping the system recover following major disturbance. Third, patterns of network connectivity change significantly as a function of yearly shelf circulation trends, and nonlinearity in network dynamics make mean connectivity a poor representation of yearly variations. From a management perspective, the priority list of reefs that are targets for management interventions depends crucially on what type of stressors (system-wide vs. localized) need addressing. This choice also depends not only on the ultimate purpose of management, but also on future oceanographic, climate change, and development scenarios that will determine the network connectivity and habitat quality.

Spatially-explicit model for assessing wild dog control strategies in Western Australia.

Large predators can significantly impact livestock industries. In Australia, wild dogs (Canis lupus familiaris, Canis lupus dingo, and hybrids) cause economic losses of more than AUD$40M annually. Landscape-scale exclusion fencing coupled with lethal techniques is a widely practiced control method. In Western Australia, the State Barrier Fence encompasses approximately 260,000km2 of predominantly agricultural land, but its effectiveness in preventing wild dogs from entering the agricultural region is difficult to evaluate. We conducted a management strategy evaluation (MSE) based on spatially-explicit population models to forecast the effects of upgrades to the Western Australian State Barrier Fence and several control scenarios varying in intensity and spatial extent on wild dog populations in southwest Western Australia. The model results indicate that populations of wild dogs on both sides of the State Barrier Fence are self-sustaining and current control practices are not sufficient to effectively reduce their abundance in the agricultural region. Only when a combination of control techniques is applied on a large scale, intensively and continuously are wild dog numbers effectively controlled. This study identifies the requirement for addressing extant populations of predators within fenced areas to meet the objective of preventing wild dog expansion. This objective is only achieved when control is applied to the whole area where wild dogs are currently present within the fence plus an additional buffer of ~20 km. Our modelling focused on the use of baiting, trapping and shooting; however, we acknowledge that additional tools may also be applied. Finally, we recommend that a cost-benefit analysis be performed to evaluate the economic viability of an integrated control strategy.

Loggerhead sea turtle embryos (Caretta caretta) regulate expression of stress response and developmental genes when exposed to a biologically realistic heat stress.

Oviparous reptile embryos are expected to breach their critical thermal maxima if temperatures reach those predicted under current climate change models due to the lack of the maternal buffering processes and parental care. Heat-shock proteins (HSPs) are integral in the molecular response to thermal stress, and their expression is heritable, but the roles of other candidate families such as the heat-shock factors (HSFs) have not been determined in reptiles. Here, we subject embryonic sea turtles (Caretta caretta) to a biologically realistic thermal stress and employ de novo transcriptomic profiling of brain tissue to investigate the underlying molecular response. From a reference transcriptome of 302 293 transcripts, 179 were identified as differentially expressed between treatments. As anticipated, genes enriched in the heat-shock treatment were primarily associated with the Hsp families, or were genes whose products play similar protein editing and chaperone functions (e.g. bag3, MYOC and serpinh1). Unexpectedly, genes encoding the HSFs were not significantly upregulated under thermal stress, indicating their presence in unstressed cells in an inactive state. Genes that were downregulated under thermal stress were less well functionally defined but were associated with stress response, development and cellular organization, suggesting that developmental processes may be compromised at realistically high temperatures. These results confirm that genes from the Hsp families play vital roles in the thermal tolerance of developing reptile embryos and, in addition with a number of other genes, should be targets for evaluating the capacity of oviparous reptiles to respond adaptively to the effects of climate change.

Seascape genomics reveals fine-scale patterns of dispersal for a reef fish along the ecologically divergent coast of Northwestern Australia.

Understanding the drivers of dispersal among populations is a central topic in marine ecology and fundamental for spatially explicit management of marine resources. The extensive coast of Northwestern Australia provides an emerging frontier for implementing new genomic tools to comparatively identify patterns of dispersal across diverse and extreme environmental conditions. Here, we focused on the stripey snapper (Lutjanus carponotatus), which is important to recreational, charter-based and customary fishers throughout the Indo-West Pacific. We collected 1,016 L. carponotatus samples at 51 locations in the coastal waters of Northwestern Australia ranging from the Northern Territory to Shark Bay and adopted a genotype-by-sequencing approach to test whether realized connectivity (via larval dispersal) was related to extreme gradients in coastal hydrodynamics. Hydrodynamic simulations using CONNIE and a more detailed treatment in the Kimberley Bioregion provided null models for comparison. Based on 4,402 polymorphic single nucleotide polymorphism loci shared across all individuals, we demonstrated significant genetic subdivision between the Shark Bay Bioregion in the south and all locations within the remaining, more northern bioregions. More importantly, we identified a zone of admixture spanning a distance of 180 km at the border of the Kimberley and Canning bioregions, including the Buccaneer Archipelago and adjacent waters, which collectively experiences the largest tropical tidal range and some of the fastest tidal currents in the world. Further testing of the generality of this admixture zone in other shallow water species across broader geographic ranges will be critical for our understanding of the population dynamics and genetic structure of marine taxa in our tropical oceans.

Inferring contemporary and historical genetic connectivity from juveniles.

Measuring population connectivity is a critical task in conservation biology. While genetic markers can provide reliable long-term historical estimates of population connectivity, scientists are still limited in their ability to determine contemporary patterns of gene flow, the most practical time frame for management. Here, we tackled this issue by developing a new approach that only requires juvenile sampling at a single time period. To demonstrate the usefulness of our method, we used the Speartooth shark (Glyphis glyphis), a critically endangered species of river shark found only in tropical northern Australia and southern Papua New Guinea. Contemporary adult and juvenile shark movements, estimated with the spatial distribution of kin pairs across and within three river systems, was contrasted with historical long-term connectivity patterns, estimated from mitogenomes and genome-wide SNP data. We found strong support for river fidelity in juveniles with the within-cohort relationship analysis. Male breeding movements were highlighted with the cross-cohort relationship analysis, and female reproductive philopatry to the river systems was revealed by the mitogenomic analysis. We show that accounting for juvenile river fidelity and female philopatry is important in population structure analysis and that targeted sampling in nurseries and juvenile aggregations should be included in the genomic toolbox of threatened species management.

Death by sex in an Australian icon: a continent-wide survey reveals extensive hybridization between dingoes and domestic dogs.

Hybridization between domesticated animals and their wild counterparts can disrupt adaptive gene combinations, reduce genetic diversity, extinguish wild populations and change ecosystem function. The dingo is a free-ranging dog that is an iconic apex predator and distributed throughout most of mainland Australia. Dingoes readily hybridize with domestic dogs, and in many Australian jurisdictions, distinct management strategies are dictated by hybrid status. Yet, the magnitude and spatial extent of domestic dog-dingo hybridization is poorly characterized. To address this, we performed a continent-wide analysis of hybridization throughout Australia based on 24 locus microsatellite DNA genotypes from 3637 free-ranging dogs. Although 46% of all free-ranging dogs were classified as pure dingoes, all regions exhibited some hybridization, and the magnitude varied substantially. The southeast of Australia was highly admixed, with 99% of animals being hybrids or feral domestic dogs, whereas only 13% of the animals from remote central Australia were hybrids. Almost all free-ranging dogs had some dingo ancestry, indicating that domestic dogs could have poor survivorship in nonurban Australian environments. Overall, wild pure dingoes remain the dominant predator over most of Australia, but the speed and extent to which hybridization has occurred in the approximately 220 years since the first introduction of domestic dogs indicate that the process may soon threaten the persistence of pure dingoes.

Y chromosome analysis of dingoes and southeast asian village dogs suggests a neolithic continental expansion from Southeast Asia followed by multiple Austronesian dispersals.

Dogs originated more than 14,000 BP, but the location(s) where they first arose is uncertain. The earliest archeological evidence of ancient dogs was discovered in Europe and the Middle East, some 5-7 millennia before that from Southeast Asia. However, mitochondrial DNA analyses suggest that most modern dogs derive from Southeast Asia, which has fueled the controversial hypothesis that dog domestication originated in this region despite the lack of supporting archeological evidence. We propose and investigate with Y chromosomes an alternative hypothesis for the proximate origins of dogs from Southeast Asia--a massive Neolithic expansion of dogs from this region that largely replaced more primitive dogs to the west and north. Previous attempts to test matrilineal findings with independent patrilineal markers have lacked the necessary genealogical resolution and mutation rate estimates. Here, we used Y chromosome genotypes, composed of 29 single-nucleotide polymorphism (SNPs) and 5 single tandem repeats (STRs), from 338 Australian dingoes, New Guinea singing dogs, and village dogs from Island Southeast Asia, along with modern European breed dogs, to estimate the evolutionary mutation rates of Y chromosome STRs based on calibration to the independently known age of the dingo population. Dingoes exhibited a unique haplogroup characterized by a single distinguishing SNP mutation and 14 STR haplotypes. The age of the European haplogroup was estimated to be only 1.7 times older than that of the dingo population, suggesting an origin during the Neolithic rather than the Paleolithic (as predicted by the Southeast Asian origins hypothesis). We hypothesize that isolation of Neolithic dogs from wolves in Southeast Asia was a key step accelerating their phenotypic transformation, enhancing their value in trade and as cargo, and enabling them to rapidly expand and replace more primitive dogs to the West. Our findings also suggest that dingoes could have arrived in Australia directly from Taiwan, independently of later dispersals of dogs through Thailand to Island Southeast Asia.

So long to genetic diversity, and thanks for all the fish.

The world faces a global fishing crisis. Wild marine fisheries comprise nearly 15% of all animal protein in the human diet, but, according to the U.N. Food and Agriculture Organization, nearly 60% of all commercially important marine fish stocks are overexploited, recovering, or depleted (FAO 2012; Fig. 1). Some authors have suggested that the large population sizes of harvested marine fish make even collapsed populations resistant to the loss of genetic variation by genetic drift (e.g. Beverton 1990). In contrast, others have argued that the loss of alleles because of overfishing may actually be more dramatic in large populations than in small ones (Ryman et al. 1995). In this issue, Pinsky & Palumbi (2014) report that overfished populations have approximately 2% lower heterozygosity and 12% lower allelic richness than populations that are not overfished. They also performed simulations which suggest that their estimates likely underestimate the actual loss of rare alleles by a factor of three or four. This important paper shows that the harvesting of marine fish can have genetic effects that threaten the long-term sustainability of this valuable resource.

Range-wide multilocus phylogeography of the red fox reveals ancient continental divergence, minimal genomic exchange and distinct demographic histories.

Widely distributed taxa provide an opportunity to compare biogeographic responses to climatic fluctuations on multiple continents and to investigate speciation. We conducted the most geographically and genomically comprehensive study to date of the red fox (Vulpes vulpes), the world's most widely distributed wild terrestrial carnivore. Analyses of 697 bp of mitochondrial sequence in ~1000 individuals suggested an ancient Middle Eastern origin for all extant red foxes and a 400 kya (SD = 139 kya) origin of the primary North American (Nearctic) clade. Demographic analyses indicated a major expansion in Eurasia during the last glaciation (~50 kya), coinciding with a previously described secondary transfer of a single matriline (Holarctic) to North America. In contrast, North American matrilines (including the transferred portion of Holarctic clade) exhibited no signatures of expansion until the end of the Pleistocene (~12 kya). Analyses of 11 autosomal loci from a subset of foxes supported the colonization time frame suggested by mtDNA (and the fossil record) but, in contrast, reflected no detectable secondary transfer, resulting in the most fundamental genomic division of red foxes at the Bering Strait. Endemic continental Y-chromosome clades further supported this pattern. Thus, intercontinental genomic exchange was overall very limited, consistent with long-term reproductive isolation since the initial colonization of North America. Based on continental divergence times in other carnivoran species pairs, our findings support a model of peripatric speciation and are consistent with the previous classification of the North American red fox as a distinct species, V. fulva.

Increased scalability and sequencing quality of an epigenetic age prediction assay.

Epigenetic ageing in a human context, has been used to better understand the relationship between age and factors such as lifestyle and genetics. In an ecological setting, it has been used to predict the age of individual animals for wildlife management. Despite the importance of epigenetic ageing in a range of research fields, the assays to measure epigenetic ageing are either expensive on a large scale or complex. In this study, we aimed to improve the efficiency and sequencing quality of an existing epigenetic ageing assay for the Australian Lungfish (Neoceratodus forsteri). We used an enzyme-based alternative to bisulfite conversion to reduce DNA fragmentation and evaluated its performance relative to bisulfite conversion. We found the sequencing quality to be 12% higher with the enzymatic alternative compared to bisulfite treatment (p-value < 0.01). This new enzymatic based approach, although currently double the cost of bisulfite treatment can increases the throughput and sequencing quality. We envisage this assay setup being adopted increasingly as the scope and scale of epigenetic ageing research continues to grow.

The genetic signature left by the range expansion of red foxes in Australia is detectable after more than 80 years of population stability.

Reconstructing biological invasions from historical sources can provide insights into how they occur but are difficult to do when invasions are poorly documented. Genetic signatures left by invaders can also offer insights into invasion routes, points of origin and general biology but often present conclusions that are contradictory to expectations. Here, we test the ability of continental-wide microsatellite genotype data from 29 loci and 3122 samples to reconstruct the well-documented invasion of red foxes Vulpes vulpes from the United Kingdom into Australia over 150 years ago, an invasion that has led to the extinction of many native species. Our analysis reveals several key signals of invasion evident in Australian foxes. They display lower levels of diversity than foxes sampled from the UK, exhibit clines in diversity from the point of introduction (south-east Australia) to the edge of their range, and show strong evidence of allele surfing in westerly and north-easterly directions. These characteristics are consistent with a single point of origin followed by rapid expansion in westerly and north-easterly directions as suggested by historical records. We also find little genetic structure in foxes across Australia with only the vast Nullarbor Plains and Great Victoria Desert region presenting a detectable barrier to their dispersal. As such, no mainland region within the current range of foxes can be considered genetically isolated and therefore appropriate for localised eradication efforts. Overall, our analyses demonstrate the ability of comprehensive population genetic studies to reconstruct invasion histories even after more than 80 years since colonisation was stabilised.

Calibrating epigenetic clocks with training data error.

Animal age data are valuable for management of wildlife populations. Yet, for most species, there is no practical method for determining the age of unknown individuals. However, epigenetic clocks, a molecular-based method, are capable of age prediction by sampling specific tissue types and measuring DNA methylation levels at specific loci. Developing an epigenetic clock requires a large number of samples from animals of known ages. For most species, there are no individuals whose exact ages are known, making epigenetic clock calibration inaccurate or impossible. For many epigenetic clocks, calibration samples with inaccurate age estimates introduce a degree of error to epigenetic clock calibration. In this study, we investigated how much error in the training data set of an epigenetic clock can be tolerated before it resulted in an unacceptable increase in error for age prediction. Using four publicly available data sets, we artificially increased the training data age error by iterations of 1% and then tested the model against an independent set of known ages. A small effect size increase (Cohen's d >0.2) was detected when the error in age was higher than 22%. The effect size increased linearly with age error. This threshold was independent of sample size. Downstream applications for age data may have a more important role in deciding how much error can be tolerated for age prediction. If highly precise age estimates are required, then it may be futile to embark on the development of an epigenetic clock when there is no accurately aged calibration population to work with. However, for other problems, such as determining the relative age order of pairs of individuals, a lower-quality calibration data set may be adequate.

Fish species lifespan prediction from promoter cytosine-phosphate-guanine density.

Lifespan is a key attribute of a species' life cycle and varies extensively among major lineages of animals. In fish, lifespan varies by several orders of magnitude, with reported values ranging from less than 1 year to approximately 400 years. Lifespan information is particularly useful for species management, as it can be used to estimate invasion potential, extinction risk and sustainable harvest rates. Despite its utility, lifespan is unknown for most fish species. This is due to the difficulties associated with accurately identifying the oldest individual(s) of a given species, and/or deriving lifespan estimates that are representative for an entire species. Recently it has been shown that CpG density in gene promoter regions can be used to predict lifespan in mammals and other vertebrates, with variable accuracy across taxa. To improve accuracy of lifespan prediction in a non-mammalian vertebrate group, here we develop a fish-specific genomic lifespan predictor. Our new model includes more than eight times the number of fish species included in the previous vertebrate model (n = 442) and uses fish-specific gene promoters as reference sequences. The model predicts fish lifespan from genomic CpG density alone (measured as CpG observed/expected ratio), explaining 64% of the variance between known and predicted lifespans. The predictions are highly robust to variation in genome quality and are applicable to all classes of fish; a taxonomically diverse and speciose group. The results demonstrate the value of promoter CpG density as a universal predictor of fish lifespan that can applied where empirical data are unavailable, or impracticable to obtain.

Ocean connectivity and habitat characteristics predict population genetic structure of seagrass in an extreme tropical setting.

Understanding patterns of gene flow and processes driving genetic differentiation is important for a broad range of conservation practices. In marine organisms, genetic differentiation among populations is influenced by a range of spatial, oceanographic, and environmental factors that are attributed to the seascape. The relative influences of these factors may vary in different locations and can be measured using seascape genetic approaches. Here, we applied a seascape genetic approach to populations of the seagrass, Thalassia hemprichii, at a fine spatial scale (~80 km) in the Kimberley coast, western Australia, a complex seascape with strong, multidirectional currents greatly influenced by extreme tidal ranges (up to 11 m, the world's largest tropical tides). We incorporated genetic data from a panel of 16 microsatellite markers, overwater distance, oceanographic data derived from predicted passive dispersal on a 2 km-resolution hydrodynamic model, and habitat characteristics from each meadow sampled. We detected significant spatial genetic structure and asymmetric gene flow, in which meadows 12-14 km apart were less connected than ones 30-50 km apart. This pattern was explained by oceanographic connectivity and differences in habitat characteristics, suggesting a combined scenario of dispersal limitation and facilitation by ocean current with local adaptation. Our findings add to the growing evidence for the key role of seascape attributes in driving spatial patterns of gene flow. Despite the potential for long-distance dispersal, there was significant genetic structuring over small spatial scales implicating dispersal and recruitment bottlenecks and highlighting the importance of implementing local-scale conservation and management measures.

Understanding age-specific dispersal in fishes through hydrodynamic modelling, genetic simulations and microsatellite DNA analysis.

Many marine species have vastly different capacities for dispersal during larval, juvenile and adult life stages, and this has the potential to complicate the identification of population boundaries and the implementation of effective management strategies such as marine protected areas. Genetic studies of population structure and dispersal rarely disentangle these differences and usually provide only lifetime-averaged information that can be considered by managers. We address this limitation by combining age-specific autocorrelation analysis of microsatellite genotypes, hydrodynamic modelling and genetic simulations to reveal changes in the extent of dispersal during the lifetime of a marine fish. We focus on an exploited coral reef species, Lethrinus nebulosus, which has a circum-tropical distribution and is a key component of a multispecies fishery in northwestern Australia. Conventional population genetic analyses revealed extensive gene flow in this species over vast distances (up to 1,500 km). Yet, when realistic adult dispersal behaviours were modelled, they could not account for these observations, implying adult dispersal does not dominate gene flow. Instead, hydrodynamic modelling showed that larval L. nebulosus are likely to be transported hundreds of kilometres, easily accounting for the observed gene flow. Despite the vast scale of larval transport, juvenile L. nebulosus exhibited fine-scale genetic autocorrelation, which declined with age. This implies both larval cohesion and extremely limited juvenile dispersal prior to maturity. The multidisciplinary approach adopted in this study provides a uniquely comprehensive insight into spatial processes in this marine fish.

Comparison of materials for rapid passive collection of environmental DNA.

Passive collection is an emerging sampling method for environmental DNA (eDNA) in aquatic systems. Passive eDNA collection is inexpensive and efficient, and requires minimal equipment, making it suited to high-density sampling and remote deployment. Here, we compare the effectiveness of nine membrane materials for passively collecting fish eDNA from a 3-million-litre marine mesocosm. We submerged materials (cellulose, cellulose with 1% and 3% chitosan, cellulose overlayed with electrospun nanofibres and 1% chitosan, cotton fibres, hemp fibres, and sponge with either zeolite or active carbon) for intervals between 5 and 1080 min. We show that for most materials, with as little as 5 min of submersion, mitochondrial fish eDNA measured with qPCR, and fish species richness measured with metabarcoding, was comparable to that collected by conventional filtering. Furthermore, PCR template DNA concentrations and species richness were generally not improved significantly by longer submersion. Species richness detected for all materials ranged between 11 and 37 species, with a median of 27, which was comparable to the range for filtered eDNA (19-32). Using scanning electron microscopy, we visualized biological matter adhering to the surface of materials, rather than entrapped, with images also revealing a diversity in size and structure of putative eDNA particles. eDNA can be collected rapidly from seawater with a passive approach and using a variety of materials. This will suit cost- and time-sensitive biological surveys, and where access to equipment is limited.

Age prediction of green turtles with an epigenetic clock.

Age is a fundamental life history attribute that is used to understand the dynamics of wild animal populations. Unfortunately, most animals do not have a practical or nonlethal method to determine age. This makes it difficult for wildlife managers to carry out population assessments, particularly for elusive and long-lived fauna such as marine turtles. In this study, we present an epigenetic clock that predicts the age of marine turtles from skin biopsies. The model was developed and validated using DNA from known-age green turtles (Chelonia mydas) from two captive populations, and mark-recapture wild turtles with known time intervals between captures. Our method, based on DNA methylation levels at 18 CpG sites, was highly accurate with a median absolute error of 2.1 years (4.7% of maximum age in data set). This is the first epigenetic clock developed for a reptile and illustrates their broad applicability across a broad variety of vertebrate species. It has the potential to transform marine turtle management through a nonlethal and inexpensive method to provide key life history information.