Sean Hoban

Computer simulations: tools for population and evolutionary genetics

Computer simulations are excellent tools for understanding the evolutionary and genetic consequences of complex processes whose interactions cannot be analytically predicted. Simulations have traditionally been used in population genetics by a fairly small community with programming expertise, but the recent availability of dozens of sophisticated, customizable software packages for simulation now makes simulation an accessible option for researchers in many fields. The in silico genetic data produced by simulations, along with greater availability of population-genomics data, are transforming genetic epidemiology, anthropology, evolutionary and population genetics and conservation. In this Review of the state-of-the-art of simulation software, we identify applications of simulations, evaluate simulator capabilities, provide a guide for their use and summarize future directions.

Breaking RAD: an evaluation of the utility of restriction site-associated DNA sequencing for genome scans of adaptation

Understanding how and why populations evolve is of fundamental importance to molecular ecology. Restriction site‐associated DNA sequencing (RADseq), a popular reduced representation method, has ushered in a new era of genome‐scale research for assessing population structure, hybridization, demographic history, phylogeography and migration. RADseq has also been widely used to conduct genome scans to detect loci involved in adaptive divergence among natural populations. Here, we examine the capacity of those RADseq‐based genome scan studies to detect loci involved in local adaptation. To understand what proportion of the genome is missed by RADseq studies, we developed a simple model using different numbers of RAD‐tags, genome sizes and extents of linkage disequilibrium (length of haplotype blocks). Under the best‐case modelling scenario, we found that RADseq using six‐ or eight‐base pair cutting restriction enzymes would fail to sample many regions of the genome, especially for species with short linkage disequilibrium. We then surveyed recent studies that have used RADseq for genome scans and found that the median density of markers across these studies was 4.08 RAD‐tag markers per megabase (one marker per 245 kb). The length of linkage disequilibrium for many species is one to three orders of magnitude less than density of the typical recent RADseq study. Thus, we conclude that genome scans based on RADseq data alone, while useful for studies of neutral genetic variation and genetic population structure, will likely miss many loci under selection in studies of local adaptation.

Bringing genetic diversity to the forefront of conservation policy and management

In this essay we explore questions on how to increase the visibility and utility of genetic information for biodiversity managers and policy makers. This is discussed in the light of Aichi CBD Target 13, which for the first time impels signatories to minimise genetic erosion and safeguard genetic diversity. Drawing on qualitative results from a questionnaire sent to European conservation professionals by the ConGRESS Framework 7 Support Action (www.congressgenetics.eu), we summarise our preliminary findings on the attitudes and experiences of European conservation professionals in using genetics. We then discuss the implications of these findings for academics involved in conservation genetics and suggest that a much closer partnership between academic conservation geneticists and conservation practitioners is necessary if the full potential of genetic tools in conservation is to be realised.

The number of markers and samples needed for detecting bottlenecks under realistic scenarios, with and without recovery: a simulation-based study

Detecting bottlenecks is a common task in molecular ecology. While several bottleneck detection methods exist, evaluations of their power have focused only on severe bottlenecks (e.g. to Ne ~10). As a component of a recent review, Peery et al. ( ) analysed the power of two approaches, the M‐ratio and heterozygote excess tests, to detect moderate bottlenecks (e.g. to Ne ~100), which is realistic for many conservation situations. In this Comment, we address three important points relevant to but not considered in Peery et al. Under moderate bottleneck scenarios, we test the (i) relative advantage of sampling more markers vs. more individuals, (ii) potential power to detect the bottleneck when utilizing dozens of microsatellites (a realistic possibility for contemporary studies) and (iii) reduction in power when postbottleneck recovery has occurred. For the realistic situations examined, we show that (i) doubling the number of loci shows equal or better power than tripling the number of individuals, (ii) increasing the number of markers (up to 100) results in continued additive gains in power, and (iii) recovery after a moderate amount of time or gradual change in size reduces power, by up to one‐half. Our results provide a practical supplement to Peery et al. and encourage the continued use of bottleneck detection methods in the genomic age, but also emphasize that the power under different sampling schemes should be estimated, using simulation modelling, as a routine component of molecular ecology studies.

Responsible RAD: Striving for best practices in population genomic studies of adaptation

Two recent articles were written in response to our paper “Breaking RAD: An evaluation of the utility of restriction site associated DNA sequencing scans of adaptation.” While we agree with some of the comments made by the authors of these two response papers, we still believe caution should be employed in RADseq studies that aim to detect loci that contribute to adaptation. In this rebuttal, we evaluate the key points made in these papers, attempt to identify a middle ground and make suggestions for responsibly conducting future studies to understand the genomewide mechanisms of adaptation.

High variance in reproductive success generates a false signature of a genetic bottleneck in populations of constant size: a simulation study

BackgroundDemographic bottlenecks can severely reduce the genetic variation of a population or a species. Establishing whether low genetic variation is caused by a bottleneck or a constantly low effective number of individuals is important to understand a species’ ecology and evolution, and it has implications for conservation management. Recent studies have evaluated the power of several statistical methods developed to identify bottlenecks. However, the false positive rate, i.e. the rate with which a bottleneck signal is misidentified in demographically stable populations, has received little attention. We analyse this type of error (type I) in forward computer simulations of stable populations having greater than Poisson variance in reproductive success (i.e., variance in family sizes). The assumption of Poisson variance underlies bottleneck tests, yet it is commonly violated in species with high fecundity.ResultsWith large variance in reproductive success (Vk ≥ 40, corresponding to a ratio between effective and census size smaller than 0.1), tests based on allele frequencies, allelic sizes, and DNA sequence polymorphisms (heterozygosity excess, M-ratio, and Tajima’s D test) tend to show erroneous signals of a bottleneck. Similarly, strong evidence of population decline is erroneously detected when ancestral and current population sizes are estimated with the model based method MSVAR.ConclusionsOur results suggest caution when interpreting the results of bottleneck tests in species showing high variance in reproductive success. Particularly in species with high fecundity, computer simulations are recommended to confirm the occurrence of a population bottleneck.

Guidelines for planning genomic assessment and monitoring of locally adaptive variation to inform species conservation

Identifying and monitoring locally adaptive genetic variation can have direct utility for conserving species at risk, especially when management may include actions such as translocations for restoration, genetic rescue, or assisted gene flow. However, genomic studies of local adaptation require careful planning to be successful, and in some cases may not be a worthwhile use of resources. Here, we offer an adaptive management framework to help conservation biologists and managers decide when genomics is likely to be effective in detecting local adaptation, and how to plan assessment and monitoring of adaptive variation to address conservation objectives. Studies of adaptive variation using genomic tools will inform conservation actions in many cases, including applications such as assisted gene flow and identifying conservation units. In others, assessing genetic diversity, inbreeding, and demographics using selectively neutral genetic markers may be most useful. And in some cases, local adaptation may be assessed more efficiently using alternative approaches such as common garden experiments. Here, we identify key considerations of genomics studies of locally adaptive variation, provide a road map for successful collaborations with genomics experts including key issues for study design and data analysis, and offer guidelines for interpreting and using results from genomic assessments to inform monitoring programs and conservation actions.

Quantitative evaluation of hybridization and the impact on biodiversity conservation

Abstract Anthropogenic hybridization is an increasing conservation threat worldwide. In South Africa, recent hybridization is threatening numerous ungulate taxa. For example, the genetic integrity of the near‐threatened bontebok (Damaliscus pygargus pygargus) is threatened by hybridization with the more common blesbok (D. p. phillipsi). Identifying nonadmixed parental and admixed individuals is challenging based on the morphological traits alone; however, molecular analyses may allow for accurate detection. Once hybrids are identified, population simulation software may assist in determining the optimal conservation management strategy, although quantitative evaluation of hybrid management is rarely performed. In this study, our objectives were to describe species‐wide and localized rates of hybridization in nearly 3,000 individuals based on 12 microsatellite loci, quantify the accuracy of hybrid assignment software (STRUCTURE and NEWHYBRIDS), and determine an optimal threshold of bontebok ancestry for management purposes. According to multiple methods, we identified 2,051 bontebok, 657 hybrids, and 29 blesbok. More than two‐thirds of locations contained at least some hybrid individuals, with populations varying in the degree of introgression. HYBRIDLAB was used to simulate four generations of coexistence between bontebok and blesbok, and to optimize a threshold of ancestry, where most hybrids will be detected and removed, and the fewest nonadmixed bontebok individuals misclassified as hybrids. Overall, a threshold Q‐value (admixture coefficient) of 0.90 would remove 94% of hybrid animals, while a threshold of 0.95 would remove 98% of hybrid animals but also 8% of nonadmixed bontebok. To this end, a threshold of 0.90 was identified as optimal and has since been implemented in formal policy by a provincial nature conservation agency. Due to widespread hybridization, effective conservation plans should be established and enforced to conserve native populations that are genetically unique.

Chloroplast haplotypes suggest preglacial differentiation and separate postglacial migration paths for the threatened North American forest tree Juglans cinerea L.

Postglacial migration paths for most of the tree species of Eastern North America remain unknown. The presence of no-analogue forest communities prior to the last glacial advance suggests that individual trees species in Eastern North America may respond differently as climate changes and human impacts increase. In this study, we examined chloroplast haplotypes from natural populations of Juglans cinerea L., a North American forest tree, to infer postglacial migration patterns. Sequences from eight different regions of the chloroplast genome in 197 trees distributed across the range revealed 10 haplotypes. A minimum spanning network, phylogenetic analysis and haplotype distributions revealed that the three most common haplotypes were geographically disjunct and not closely related. Haplotype 6 (73 trees) occurred only in western populations, haplotype 10 (83 trees) occurred only in eastern populations and haplotype 7 (21 trees) occurred only at the southern edge of the native range. The southernmost population contained the most haplotype diversity but included no eastern haplotypes. Haplotype phylogeography suggested geographical differentiation prior to the last glacial advance in eastern populations and separate postglacial migration paths for eastern and western populations. As migration of J. cinerea to Atlantic Canada from southern refugia does not appear possible given known seed dispersal mechanisms, the possibility of northern refugia or dispersal by extinct megafauna merits serious consideration. Differences among species in preglacial history, ecological niche preferences and seed dispersal mechanisms suggest that response to long-term climate change and acute human disturbance may be highly species specific.

Challenges in global biodiversity conservation and solutions that cross sociology, politics, economics and ecology.

The study and practice of conservation biology is inherently interdisciplinary, addresses short and long time-scales and occurs within complex human–natural interfaces. Zoos and aquaria, in partnership with researchers, other non-government organizations, government, industry and educators, are combining knowledge of species and ecosystems with economics, psychology and law to create solutions for conserving biodiversity. From 22 to 25 May, the Conservation Forum of the European Association of Zoos and Aquaria was a venue for discussing conservation research, education and interventions, from the scale of villages to global policy.