Next generation conservation genetics at sea: detecting and conserving adaptive potential


There is a long tradition in conservation genetics of using neutral genetic markers to assess demographic patterns and processes relevant to conservation. For example, the identification of the units of conservation requires information on population genetic structure and patterns of connectivity. These analyses assume differentiation by neutral processes. Neutral markers also provide information on effective population size and historical population dynamics, and reflect inbreeding (and can therefore be correlated with fitness; Reed & Frankham 2003). At the same time, a primary objective (and almost a mantra) of conservation genetics has been the conservation of adaptive potential. This is necessarily associated with phenotype and local adaptation, each of which will be determined by functional genetic loci, and gene-environment interactions. The objective is to conserve the ability of wildlife species to respond to changing environments, essential in the context of anthropogenic global climate change. However, this has been difficult to assess and support. Although there is weak correlation between quantitative trait diversity (QST, reflecting phenotypic diversity) and Wright’s inbreeding coefficient measure of population structure (FST, based on neutral markers; McKay & Latta 2002), a clear relationship has remained elusive. Variation at adaptive markers could be expected to show a more direct relationship, but the relationship between phenotype and genotype is poorly understood. A meta-analyses of plant studies showed a stronger relationship between adaptive diversity and ecological and community structure than for neutral markers in plant species using indirect inference for adaptation (e.g. ecotypes and cultivars; Whitlock 2014), but to better understand the broader relationship we can now use genomic data.
In the marine environment there are few obvious barriers to gene flow, and in fact marine species of management concern (such as pelagic fish species) often show little evidence of structure at neutral genetic markers. This was the case for several deep sea fish species studied in our lab, where 15-20 polymorphic microsatellite DNA markers showed no evidence of structure across a broad geographic range (e.g. White et al. 2009). However, for each of these species, genome sampling (by restriction-associated DNA analysis: RADseq), showed evidence of weak structure at neutral loci, consistent with isolation by distance in most cases. However, RADseq also permits an analysis of structure at putative functional loci, detected as outliers from predicted divergence metrics under neutral theory. At those loci, a different pattern of structure emerged, highlighting especially stand out differentiation for populations south of the sub-polar front (in warmer waters, influenced by a different current system; e.g. see Figure 1; Goncalves da Silva et al. 2019). Some similar patterns associated instead with ecotype were found for marine mammal species (e.g. Moura et al. 2015). The differential pattern of structure suggests local adaptation, but the more specific mechanisms associated with specific loci or genetic pathways could not easily be determined from the RADseq data (where 1000’s of single nucleotide polymorphic (SNP) sites are revealed, but millions may be required).

Aims & novelty:
To test hypotheses about the more specific mechanisms that generate distinct patterns of diversity at functional loci across environmental gradients using high resolution data and working with reference genomes. The key objectives will be to identify loci or gene systems that reflect diversity requiring conservation management that are not detectable using conventional methods based on neutral loci. At the same time, these data will enhance our understanding of the process of adaptive evolution, providing novel inference about the evolution of diversity in the marine ecosystem.

Click on an image to expand

Image Captions

Redesigned from figure in Goncalves da Silva et al. 2019


Extensive existing sample sets across relevant environmental boundaries will be investigated using high-resolution re-sequencing methodologies on an Illumina platform (available in Durham), in comparison with new datasets as required and appropriate. Analyses will be based on comparison against available annotated reference sequences using GWAS and similar bioinformatic methods. Neutral patterns of population structure and diversity will be compared with signals of local adaptation and admixture using dense SNP databases and coalescent, site frequency spectra, ordination, redundancy analyses, and other methods that consider genetic diversity in the context of ecological and environmental variables.

Project Timeline

Year 1

The student will research the background literature, be trained in relevant methodologies, manage collections for study, extract DNA as needed and begin DNA sequencing work.

Year 2

Focus on next generation sequencing and bioinformatics, together with database searches and initial analytical work.

Year 3

Further bioinformatics, statistical analyses, data interpretation and thesis writing.

Year 3.5

Complete thesis writing and work on generating published papers.

& Skills

In addition to training in Durham as part of the Biosciences program, together with elective training courses through the Durham University, the student will undertake a course on next generation sequencing bioinformatics through the Edinburgh Genomics facility.

References & further reading

Goncalves da Silva, A., Barendse, B., Kijas, J., England, P.R., Hoelzel, A.R. (2019) Genomic data suggest environmental drivers of fish population structure in the deep sea; a case study for the orange roughy (Hoplostethus atlanticus) J. Applied Ecol. In press.

McKay, J.K., Latta, R.G. (2002) Adaptive population divergence: markers, QTL and traits. TREE 17, 285-291.

Moura, A.E., Kenny, J.G., Chaudhuri, R., Hughes, M.A., Welch, A., Reisinger, R.R., de Bruyn, P.J.N., Dahlheim, M.E., Hall, N.,

Hoelzel, A.R. 2014. Population genomics of the killer whale indicates ecotype evolution in sympatry involving both selection and drift. Mol. Ecol. 23, 5179-5192

Reed, D.H., Frankham, R. (2003) Correlation between fitness and genetic diversity. Cons. Biol. 17, 230-237.

White, T.A., Stefanni, S., Stamford, J. & Hoelzel, A.R. 2009. Ocean basin panmixia in a long-lived, deep-sea fish with well defined habitat dependence and relatively low fecundity. Mol. Ecol. 18: 2563-2573.

Whitlock, R. (2014) Relationships between adaptive and neutral genetic diversity and ecological structure and functioning: a meta-analysis. J. Ecol. 102, 857-872.

Further Information

Supervisor contact details: Prof. Rus Hoelzel, Department of Biosciences, Durham University (; Dr. Per Berggren, Natural and Environmental Sciences, Newcastle University (

Apply Now