Brendan N. Reid

Reliability of genetic bottleneck tests for detecting recent population declines

The identification of population bottlenecks is critical in conservation because populations that have experienced significant reductions in abundance are subject to a variety of genetic and demographic processes that can hasten extinction. Genetic bottleneck tests constitute an appealing and popular approach for determining if a population decline has occurred because they only require sampling at a single point in time, yet reflect demographic history over multiple generations. However, a review of the published literature indicates that, as typically applied, microsatellite‐based bottleneck tests often do not detect bottlenecks in vertebrate populations known to have experienced declines. This observation was supported by simulations that revealed that bottleneck tests can have limited statistical power to detect bottlenecks largely as a result of limited sample sizes typically used in published studies. Moreover, commonly assumed values for mutation model parameters do not appear to encompass variation in microsatellite evolution observed in vertebrates and, on average, the proportion of multi‐step mutations is underestimated by a factor of approximately two. As a result, bottleneck tests can have a higher probability of ‘detecting’ bottlenecks in stable populations than expected based on the nominal significance level. We provide recommendations that could add rigor to inferences drawn from future bottleneck tests and highlight new directions for the characterization of demographic history.

Comparing and combining distance-based and character-based approaches for barcoding turtles.

Molecular barcoding can serve as a powerful tool in wildlife forensics and may prove to be a vital aid in conserving organisms that are threatened by illegal wildlife trade, such as turtles (Order Testudines). We produced cytochrome oxidase subunit one (COI) sequences (650 bp) for 174 turtle species and combined these with publicly available sequences for 50 species to produce a data set representative of the breadth of the order. Variability within the barcode region was assessed, and the utility of both distance‐based and character‐based methods for species identification was evaluated. For species in which genetic material from more than one individual was available (n = 69), intraspecific divergences were 1.3% on average, although divergences greater than the customary 2% barcode threshold occurred within 15 species. High intraspecific divergences could indicate species with a high degree of internal genetic structure or possibly even cryptic species, although introgression is also probable in some of these taxa. Divergences between species of the same genus were 6.4% on average; however, 49 species were <2% divergent from congeners. Low levels of interspecific divergence could be caused by recent evolutionary radiations coupled with the low rates of mtDNA evolution previously observed in turtles. Complementing distance‐based barcoding with character‐based methods for identifying diagnostic sets of nucleotides provided better resolution in several cases where distance‐based methods failed to distinguish species. An online identification engine was created to provide character‐based identifications. This study constitutes the first comprehensive barcoding effort for this seriously threatened order.

DNA barcodes for globally threatened marine turtles: a registry approach to documenting biodiversity.

DNA barcoding is a global initiative that provides a standardized and efficient tool to catalogue and inventory biodiversity, with significant conservation applications. Despite progress across taxonomic realms, globally threatened marine turtles remain underrepresented in this effort. To obtain DNA barcodes of marine turtles, we sequenced a segment of the cytochrome c oxidase subunit I (COI) gene from all seven species in the Atlantic and Pacific Ocean basins (815 bp; n = 80). To further investigate intraspecific variation, we sequenced green turtles (Chelonia mydas) from nine additional Atlantic/Mediterranean nesting areas (n = 164) and from the Eastern Pacific (n = 5). We established character‐based DNA barcodes for each species using unique combinations of character states at 76 nucleotide positions. We found that no haplotypes were shared among species and the mean of interspecific variation ranged from 1.68% to 13.0%, and the mean of intraspecific variability was relatively low (0–0.90%). The Eastern Pacific green turtle sequence was identical to an Australian haplotype, suggesting that this marker is not appropriate for identifying these phenotypically distinguishable populations. Analysis of COI revealed a north–south gradient in green turtles of Western Atlantic/Mediterranean nesting areas, supporting a hypothesis of recent dispersal from near equatorial glacial refugia. DNA barcoding of marine turtles is a powerful tool for species identification and wildlife forensics, which also provides complementary data for conservation genetic research.

Finding the right coverage: the impact of coverage and sequence quality on single nucleotide polymorphism genotyping error rates

Restriction‐enzyme‐based sequencing methods enable the genotyping of thousands of single nucleotide polymorphism (SNP) loci in nonmodel organisms. However, in contrast to traditional genetic markers, genotyping error rates in SNPs derived from restriction‐enzyme‐based methods remain largely unknown. Here, we estimated genotyping error rates in SNPs genotyped with double digest RAD sequencing from Mendelian incompatibilities in known mother–offspring dyads of Hoffmans two‐toed sloth (Choloepus hoffmanni) across a range of coverage and sequence quality criteria, for both reference‐aligned and de novo‐assembled data sets. Genotyping error rates were more sensitive to coverage than sequence quality and low coverage yielded high error rates, particularly in de novo‐assembled data sets. For example, coverage ≥5 yielded median genotyping error rates of ≥0.03 and ≥0.11 in reference‐aligned and de novo‐assembled data sets, respectively. Genotyping error rates declined to ≤0.01 in reference‐aligned data sets with a coverage ≥30, but remained ≥0.04 in the de novo‐assembled data sets. We observed approximately 10‐ and 13‐fold declines in the number of loci sampled in the reference‐aligned and de novo‐assembled data sets when coverage was increased from ≥5 to ≥30 at quality score ≥30, respectively. Finally, we assessed the effects of genotyping coverage on a common population genetic application, parentage assignments, and showed that the proportion of incorrectly assigned maternities was relatively high at low coverage. Overall, our results suggest that the trade‐off between sample size and genotyping error rates be considered prior to building sequencing libraries, reporting genotyping error rates become standard practice, and that effects of genotyping errors on inference be evaluated in restriction‐enzyme‐based SNP studies.

Resolving the phylogenetic history of the short-necked turtles, genera Elseya and Myuchelys (Testudines: Chelidae) from Australia and New Guinea

Phylogenetic relationships and taxonomy of the short-necked turtles of the genera Elseya, Myuchelys, and Emydura in Australia and New Guinea have long been debated as a result of conflicting hypotheses supported by different data sets and phylogenetic analyses. To resolve this contentious issue, we analyzed sequences from two mitochondrial genes (cytochrome b and ND4) and one nuclear intron gene (R35) from all species of the genera Elseya, Myuchelys, Emydura, and their relatives. Phylogenetic analyses using three methods (maximum parsimony, maximum likelihood, and Bayesian inference) produce a single, well resolved, and strongly corroborated hypothesis, which provides support for the three genera, with the exception that the genus Myuchelys is paraphyletic - Myuchelys purvisi is the sister taxon to the remaining Elseya, Myuchelys and Emydura. A new genus is proposed for the species Myuchelys purvisi to address this paraphyletic relationship. Time-calibration analysis suggests that diversification of the group in Australia coincides with periods of aridification in the late Eocene and between the mid-Miocene and early Pliocene. Other speciation events occurred during the faunal exchange between Australia and the island of New Guinea during the late Miocene and early Pliocene. Lineages distributed in New Guinea are likely influenced by the complex geologic history of the island, and include cryptic species diversity.

Genetic effects of landscape, habitat preference and demography on three co-occurring turtle species

Expanding the scope of landscape genetics beyond the level of single species can help to reveal how species traits influence responses to environmental change. Multispecies studies are particularly valuable in highly threatened taxa, such as turtles, in which the impacts of anthropogenic change are strongly influenced by interspecific differences in life history strategies, habitat preferences and mobility. We sampled approximately 1500 individuals of three co‐occurring turtle species across a gradient of habitat change (including varying loss of wetlands and agricultural conversion of upland habitats) in the Midwestern USA. We used genetic clustering and multiple regression methods to identify associations between genetic structure and permanent landscape features, past landscape composition and landscape change in each species. Two aquatic generalists (the painted turtle, Chrysemys picta, and the snapping turtle Chelydra serpentina) both exhibited population genetic structure consistent with isolation by distance, modulated by aquatic landscape features. Genetic divergence for the more terrestrial Blandings turtle (Emydoidea blandingii), on the other hand, was not strongly associated with geographic distance or aquatic features, and Bayesian clustering analysis indicated that many Emydoidea populations were genetically isolated. Despite long generation times, all three species exhibited associations between genetic structure and postsettlement habitat change, indicating that long generation times may not be sufficient to delay genetic drift resulting from recent habitat fragmentation. The concordances in genetic structure observed between aquatic species, as well as isolation in the endangered, long‐lived Emydoidea, reinforce the need to consider both landscape composition and demographic factors in assessing differential responses to habitat change in co‐occurring species.

More precisely biased: increasing the number of markers is not a silver bullet in genetic bottleneck testing

In response to our review of the use of genetic bottleneck tests in the conservation literature (Peery et al. 2012, Molecular Ecology, 21, 3403–3418), Hoban et al. (2013, Molecular Ecology, in press) conducted population genetic simulations to show that the statistical power of genetic bottleneck tests can be increased substantially by sampling large numbers of microsatellite loci, as they suggest is now possible in the age of genomics. While we agree with Hoban and co‐workers in principle, sampling large numbers of microsatellite loci can dramatically increase the probability of committing type 1 errors (i.e. detecting a bottleneck in a stable population) when the mutation model is incorrectly assumed. Using conservative values for mutation model parameters can reduce the probability of committing type 1 errors, but doing so can result in significant losses in statistical power. Moreover, we believe that practical limitations associated with developing large numbers of high‐quality microsatellite loci continue to constrain sample sizes, a belief supported by a literature review of recent studies using next generation sequencing methods to develop microsatellite libraries. conclusion, we maintain that researchers employing genetic bottleneck tests should proceed with caution and carefully assess both statistical power and type 1 error rates associated with their study design.

Individual reproductive strategies shape the mating system of tree sloths

While polygyny is the dominant mating system in mammals, it is increasingly recognized that promiscuity occurs in most species. Using a long-term genetic and space-use data set, we documented the mating system for 2 sedentary and uniparous species of tree sloths, brown-throated three-toed (Bradypus variegatus) and Hoffmanns two-toed (Choloepus hoffmanni) sloths. We predicted that the life history of these species facilitates female strategies that promote mating with multiple males across breeding seasons, and shape central features of the mating system in tree sloths. We found that many female sloths mated with different males during our study: 70% of female B. variegatus and 50% of female C. hoffmanni switched mates among years at least once during our study. Our observations of individual movements suggested that females employed 2 strategies that appeared to influence mate switching across breeding seasons: 1) selecting a male from a pool of males in their activity center, and 2) mating with different males by shifting their home ranges during estrus. Collectively, our findings suggest that individual variation in female reproductive strategies contributes to shaping the mating systems for a sedentary mammal like sloths, and highlights the need for long-term studies to effectively capture the mating systems of mammals with slow life histories. Aunque la poligamia es el sistema de apareamiento dominante en los mamíferos, es cada vez más reconocido que las hembras en la mayoría de las especies de mamíferos se aparean con múltiples machos. Usando un set de datos genéticos y de uso de hábitat a largo plazo, se documentó el sistema de apareamiento de dos especies sedentarias y uníparas, el perezoso grisáceo (Bradypus variegatus) y el perico ligero (Choloepus hoffmanni). Se predijo que las características de estas especies, facilitan que las hembras presenten estrategias que promuevan el apareamiento con múltiples machos, determinando finalmente el sistema de apareamiento en estas especies. Se encontró que muchas hembras se aparearon con múltiples machos: 70% de las hembras de B. variegatus y el 50% de las de C. hoffmanni presentaron múltiples parejas por lo menos una vez durante esta investigación. Nuestras observaciones sobre movimientos individuales revelaron que las hembras cambian de pareja bajo diferentes escenarios: (1) apareamiento con machos diferentes en cada estro, de los que se superponen en su ámbito hogareño (2), apareamiento con machos diferentes al cambiar su ámbito hogareño durante el estro. Estas evidencias sugieren que las estrategias reproductivas individuales de las hembras ayudan a modelar el sistema reproductivo de especies sedentarias como los perezosos; finalmente, es importante destacar la necesidad de investigaciones a largo plazo para poder entender los sistemas de apareamiento de especies con lenta historia de vida.

Finding the right coverage : The impact of coverage and sequence quality on SNP genotyping error rates

Restriction‐enzyme‐based sequencing methods enable the genotyping of thousands of single nucleotide polymorphism (SNP) loci in nonmodel organisms. However, in contrast to traditional genetic markers, genotyping error rates in SNPs derived from restriction‐enzyme‐based methods remain largely unknown. Here, we estimated genotyping error rates in SNPs genotyped with double digest RAD sequencing from Mendelian incompatibilities in known mother–offspring dyads of Hoffmans two‐toed sloth (Choloepus hoffmanni) across a range of coverage and sequence quality criteria, for both reference‐aligned and de novo‐assembled data sets. Genotyping error rates were more sensitive to coverage than sequence quality and low coverage yielded high error rates, particularly in de novo‐assembled data sets. For example, coverage ≥5 yielded median genotyping error rates of ≥0.03 and ≥0.11 in reference‐aligned and de novo‐assembled data sets, respectively. Genotyping error rates declined to ≤0.01 in reference‐aligned data sets with a coverage ≥30, but remained ≥0.04 in the de novo‐assembled data sets. We observed approximately 10‐ and 13‐fold declines in the number of loci sampled in the reference‐aligned and de novo‐assembled data sets when coverage was increased from ≥5 to ≥30 at quality score ≥30, respectively. Finally, we assessed the effects of genotyping coverage on a common population genetic application, parentage assignments, and showed that the proportion of incorrectly assigned maternities was relatively high at low coverage. Overall, our results suggest that the trade‐off between sample size and genotyping error rates be considered prior to building sequencing libraries, reporting genotyping error rates become standard practice, and that effects of genotyping errors on inference be evaluated in restriction‐enzyme‐based SNP studies.

Finding the right coverage

Restriction‐enzyme‐based sequencing methods enable the genotyping of thousands of single nucleotide polymorphism (SNP) loci in nonmodel organisms. However, in contrast to traditional genetic markers, genotyping error rates in SNPs derived from restriction‐enzyme‐based methods remain largely unknown. Here, we estimated genotyping error rates in SNPs genotyped with double digest RAD sequencing from Mendelian incompatibilities in known mother–offspring dyads of Hoffmans two‐toed sloth (Choloepus hoffmanni) across a range of coverage and sequence quality criteria, for both reference‐aligned and de novo‐assembled data sets. Genotyping error rates were more sensitive to coverage than sequence quality and low coverage yielded high error rates, particularly in de novo‐assembled data sets. For example, coverage ≥5 yielded median genotyping error rates of ≥0.03 and ≥0.11 in reference‐aligned and de novo‐assembled data sets, respectively. Genotyping error rates declined to ≤0.01 in reference‐aligned data sets with a coverage ≥30, but remained ≥0.04 in the de novo‐assembled data sets. We observed approximately 10‐ and 13‐fold declines in the number of loci sampled in the reference‐aligned and de novo‐assembled data sets when coverage was increased from ≥5 to ≥30 at quality score ≥30, respectively. Finally, we assessed the effects of genotyping coverage on a common population genetic application, parentage assignments, and showed that the proportion of incorrectly assigned maternities was relatively high at low coverage. Overall, our results suggest that the trade‐off between sample size and genotyping error rates be considered prior to building sequencing libraries, reporting genotyping error rates become standard practice, and that effects of genotyping errors on inference be evaluated in restriction‐enzyme‐based SNP studies.