Ethical considerations when co-analyzing ancient DNA and data from private genetic databases.

Ancient DNA studies have begun to explore the possibility of identifying identical DNA segments shared between historical and living people. This research requires access to large genetic datasets to maximize the likelihood of identifying previously unknown, close genetic connections. Direct-to-consumer genetic testing companies, such as 23andMe, Inc., manage by far the largest and most diverse genetic databases that can be used for this purpose. It is therefore important to think carefully about guidelines for carrying out collaborations between researchers and such companies. Such collaborations require consideration of ethical issues, including policies for sharing ancient DNA datasets, and ensuring reproducibility of research findings when access to privately controlled genetic datasets is limited. At the same time, they introduce unique possibilities for returning results to the research participants whose data are analyzed, including those who are identified as close genetic relatives of historical individuals, thereby enabling ancient DNA research to contribute to the restoration of information about ancestral connections that were lost over time, which can be particularly meaningful for families and groups where such history has not been well documented. We explore these issues by describing our experience designing and carrying out a study searching for genetic connections between 18th- and 19th-century enslaved and free African Americans who labored at Catoctin Furnace, Maryland, and 23andMe research participants. We share our experience in the hope of helping future researchers navigate similar ethical considerations, recognizing that our perspective is part of a larger conversation about best ethical practices.

Genetic analysis and natural history of Parkinson's disease due to the LRRK2 G2019S variant.

The LRRK2 G2019S variant is the most common cause of monogenic Parkinson's disease (PD); however, questions remain regarding the penetrance, clinical phenotype and natural history of carriers. We performed a 3.5-year prospective longitudinal online study in a large number of 1286 genotyped LRRK2 G2019S carriers and 109 154 controls, with and without PD, recruited from the 23andMe Research Cohort. We collected self-reported motor and non-motor symptoms every 6 months, as well as demographics, family histories and environmental risk factors. Incident cases of PD (phenoconverters) were identified at follow-up. We determined lifetime risk of PD using accelerated failure time modelling and explored the impact of polygenic risk on penetrance. We also computed the genetic ancestry of all LRRK2 G2019S carriers in the 23andMe database and identified regions of the world where carrier frequencies are highest. We observed that despite a 1 year longer disease duration (P = 0.016), LRRK2 G2019S carriers with PD had similar burden of motor symptoms, yet significantly fewer non-motor symptoms including cognitive difficulties, REM sleep behaviour disorder (RBD) and hyposmia (all P-values ≤ 0.0002). The cumulative incidence of PD in G2019S carriers by age 80 was 49%. G2019S carriers had a 10-fold risk of developing PD versus non-carriers. This rose to a 27-fold risk in G2019S carriers with a PD polygenic risk score in the top 25% versus non-carriers in the bottom 25%. In addition to identifying ancient founding events in people of North African and Ashkenazi descent, our genetic ancestry analyses infer that the G2019S variant was later introduced to Spanish colonial territories in the Americas. Our results suggest LRRK2 G2019S PD appears to be a slowly progressive predominantly motor subtype of PD with a lower prevalence of hyposmia, RBD and cognitive impairment. This suggests that the current prodromal criteria, which are based on idiopathic PD, may lack sensitivity to detect the early phases of LRRK2 PD in G2019S carriers. We show that polygenic burden may contribute to the development of PD in the LRRK2 G2019S carrier population. Collectively, the results should help support screening programmes and candidate enrichment strategies for upcoming trials of LRRK2 inhibitors in early-stage disease.

Genetic Consequences of the Transatlantic Slave Trade in the Americas.

According to historical records of transatlantic slavery, traders forcibly deported an estimated 12.5 million people from ports along the Atlantic coastline of Africa between the 16th and 19th centuries, with global impacts reaching to the present day, more than a century and a half after slavery's abolition. Such records have fueled a broad understanding of the forced migration from Africa to the Americas yet remain underexplored in concert with genetic data. Here, we analyzed genotype array data from 50,281 research participants, which-combined with historical shipping documents-illustrate that the current genetic landscape of the Americas is largely concordant with expectations derived from documentation of slave voyages. For instance, genetic connections between people in slave trading regions of Africa and disembarkation regions of the Americas generally mirror the proportion of individuals forcibly moved between those regions. While some discordances can be explained by additional records of deportations within the Americas, other discordances yield insights into variable survival rates and timing of arrival of enslaved people from specific regions of Africa. Furthermore, the greater contribution of African women to the gene pool compared to African men varies across the Americas, consistent with literature documenting regional differences in slavery practices. This investigation of the transatlantic slave trade, which is broad in scope in terms of both datasets and analyses, establishes genetic links between individuals in the Americas and populations across Atlantic Africa, yielding a more comprehensive understanding of the African roots of peoples of the Americas.

Mixed support for gene flow as a constraint to local adaptation and contributor to the limited geographic range of an endemic salamander.

Understanding mechanisms that underlie species range limits is at the core of evolutionary ecology. Asymmetric gene flow between larger core populations and smaller edge populations can swamp local adaptation at the range edge and inhibit further range expansion. However, empirical tests of this theory are exceedingly rare. We tested the hypothesis that asymmetric gene flow can constrain local adaptation and thereby species' range limits in an endemic US salamander (Ambystoma barbouri) by determining if gene flow is asymmetric between the core and peripheries of the species' geographic distribution and testing whether local adaptation is swamped at range edges with a reciprocal transplant experiment. Using putatively neutral loci from populations across three core-to-edge transects that covered nearly the entire species' geographic range, we found evidence for asymmetric, core-to-edge gene flow along western and northern transects, but not along a southern transect. Subsequently, the reciprocal transplant experiment suggested that northern and western edge populations are locally adapted despite experiencing asymmetric gene flow, yet have lower fitness in their respective home regions than those of centre population. Conversely, southern populations exhibit low deme quality, experiencing high mortality regardless of where they were reared, probably due to harsher edge habitat conditions. Consequently, we provide rare species-wide evidence that local adaptation can occur despite asymmetric gene flow, though migration from the core may prohibit range expansion by reducing fitness in edge populations. Further, our multitransect study shows that multiple, nonmutually exclusive mechanisms can lead to range limits within a single species.

Selection at a genomic region of major effect is responsible for evolution of complex life histories in anadromous steelhead.

BACKGROUND: Disparity in the timing of biological events occurs across a variety of systems, yet the understanding of genetic basis underlying diverse phenologies remains limited. Variation in maturation timing occurs in steelhead trout, which has been associated with greb1L, an oestrogen target gene. Previous techniques that identified this gene only accounted for about 0.5-2.0% of the genome and solely investigated coastal populations, leaving uncertainty on the genetic basis of this trait and its prevalence across a larger geographic scale. RESULTS: We used a three-tiered approach to interrogate the genomic basis of complex phenology in anadromous steelhead. First, fine scale mapping with 5.3 million SNPs from resequencing data covering 68% of the genome confirmed a 309-kb region consisting of four genes on chromosome 28, including greb1L, to be the genomic region of major effect for maturation timing. Second, broad-scale characterization of candidate greb1L genotypes across 59 populations revealed unexpected patterns in maturation phenology for inland fish migrating long distances relative to those in coastal streams. Finally, genotypes from 890 PIT-tag tracked steelhead determined associations with early versus late arrival to spawning grounds that were previously unknown. CONCLUSIONS: This study clarifies the genetic bases for disparity in phenology observed in steelhead, determining an unanticipated trait association with premature versus mature arrival to spawning grounds and identifying multiple candidate genes potentially contributing to this variation from a single genomic region of major effect. This illustrates how dense genome mapping and detailed phenotypic characterization can clarify genotype to phenotype associations across geographic ranges of species.

Genomic variation underlying complex life-history traits revealed by genome sequencing in Chinook salmon.

A broad portfolio of phenotypic diversity in natural organisms can buffer against exploitation and increase species persistence in disturbed ecosystems. The study of genomic variation that accounts for ecological and evolutionary adaptation can represent a powerful approach to extend understanding of phenotypic variation in nature. Here we present a chromosome-level reference genome assembly for Chinook salmon (Oncorhynchus tshawytscha; 2.36 Gb) that enabled association mapping of life-history variation and phenotypic traits for this species. Whole-genome re-sequencing of populations with distinct life-history traits provided evidence that divergent selection was extensive throughout the genome within and among phylogenetic lineages, indicating that a broad portfolio of phenotypic diversity exists in this species that is related to local adaptation and life-history variation. Association mapping with millions of genome-wide SNPs revealed that a genomic region of major effect on chromosome 28 was associated with phenotypes for premature and mature arrival to spawning grounds and was consistent across three distinct phylogenetic lineages. Our results demonstrate how genomic resources can enlighten the genetic basis of known phenotypes in exploited species and assist in clarifying phenotypic variation that may be difficult to observe in naturally occurring organisms.

Landscape features along migratory routes influence adaptive genomic variation in anadromous steelhead (Oncorhynchus mykiss).

Organisms typically show evidence of adaptation to features within their local environment. However, many species undergo long-distance dispersal or migration across larger geographic regions that consist of highly heterogeneous habitats. Therefore, selection may influence adaptive genetic variation associated with landscape features at residing sites and along migration routes in migratory species. We tested for genomic adaptation to landscape features at natal spawning sites and along migration paths to the ocean of anadromous steelhead trout (Oncorhynchus mykiss) in the Columbia River Basin. Results from multivariate ordination, gene-environment association and outlier analyses using 24,526 single nucleotide polymorphisms (SNPs) provided evidence that adaptive allele frequencies were more commonly associated with landscape features along migration paths than features at natal sites (91.8% vs. 8.2% of adaptive loci, respectively). Among the 45 landscape variables tested, migration distance to the ocean and mean annual precipitation along migration paths were significantly associated with adaptive genetic variation in three distinct genetic groups. Additionally, variables such as minimum migration water temperature and mean migration slope were significant only in inland stocks of steelhead that migrate up to 1,200 km farther than those near the coast, indicating regional differences in migratory selective pressures. This study provides novel approaches for investigating migratory corridors and some of the first evidence that environment along migration paths can lead to substantial divergent selection. Consequently, our approach to understand genetic adaptation to migration conditions can be applied to other migratory species when migration or dispersal paths are generally known.

An approach for identifying cryptic barriers to gene flow that limit species' geographic ranges.

Species' geographic range limits are most often not demarcated by obvious dispersal barriers. Poor-quality habitat at the edge of a species' range can prevent range expansion by preventing outward migration or through reducing adaptive potential resulting from decreased genetic diversity. We identified habitat variables that constrain gene flow across the entire geographic range of an endemic salamander (Ambystoma barbouri) in the eastern United States, and we tested whether increased resistance resulting from these variables provides cryptic dispersal barriers at the range edges. Using polymorphic microsatellite loci, we first identified three genetic clusters that are separated by the Ohio and Kentucky rivers. Through a combination of landscape genetic analyses and generalized dissimilarity modelling, we then classified variables that (i) restrict gene flow in each of the genetic clusters across the geographic distribution of A. barbouri and (ii) become more common towards the peripheries of the distribution. A decrease in limestone availability and an increase in growing season precipitation were correlated with high resistance to gene flow across the range, and both became more common at the edges of the species' distribution. However, other landscape variables were more important for explaining variation in geneflow rates in different portions of the range, such as increased mean annual temperature and frost-free period in the south vs. growing season precipitation in the north. Taken together, these results suggest that there are both range-wide and regionally specific cryptic habitat barriers preventing geographic range expansion. Species 'geographic range limits are probably governed by a set of ecological and evolutionary factors, and our landscape genetic approach could be applied to gain additional insight into many systems.

A test of the central-marginal hypothesis using population genetics and ecological niche modelling in an endemic salamander (Ambystoma barbouri).

The central-marginal hypothesis (CMH) predicts that population size, genetic diversity and genetic connectivity are highest at the core and decrease near the edges of species' geographic distributions. We provide a test of the CMH using three replicated core-to-edge transects that encompass nearly the entire geographic range of the endemic streamside salamander (Ambystoma barbouri). We confirmed that the mapped core of the distribution was the most suitable habitat using ecological niche modelling (ENM) and via genetic estimates of effective population sizes. As predicted by the CMH, we found statistical support for decreased genetic diversity, effective population size and genetic connectivity from core to edge in western and northern transects, yet not along a southern transect. Based on our niche model, habitat suitability is lower towards the southern range edge, presumably leading to conflicting core-to-edge genetic patterns. These results suggest that multiple processes may influence a species' distribution based on the heterogeneity of habitat across a species' range and that replicated sampling may be needed to accurately test the CMH. Our work also emphasizes the importance of identifying the geographic range core with methods other than using the Euclidean centre on a map, which may help to explain discrepancies among other empirical tests of the CMH. Assessing core-to-edge population genetic patterns across an entire species' range accompanied with ENM can inform our general understanding of the mechanisms leading to species' geographic range limits.

PoolParty2: An integrated pipeline for analysing pooled or indexed low-coverage whole-genome sequencing data to discover the genetic basis of diversity.

Whole-genome sequencing data allow survey of variation from across the genome, reducing the constraint of balancing genome sub-sampling with estimating recombination rates and linkage between sampled markers and target loci. As sequencing costs decrease, low-coverage whole-genome sequencing of pooled or indexed-individual samples is commonly utilized to identify loci associated with phenotypes or environmental axes in non-model organisms. There are, however, relatively few publicly available bioinformatic pipelines designed explicitly to analyse these types of data, and fewer still that process the raw sequencing data, provide useful metrics of quality control and then execute analyses. Here, we present an updated version of a bioinformatics pipeline called PoolParty2 that can effectively handle either pooled or indexed DNA samples and includes new features to improve computational efficiency. Using simulated data, we demonstrate the ability of our pipeline to recover segregating variants, estimate their allele frequencies accurately, and identify genomic regions harbouring loci under selection. Based on the simulated data set, we benchmark the efficacy of our pipeline with another bioinformatic suite, angsd, and illustrate the compatibility and complementarity of these suites using angsd to generate genotype likelihoods as input for identifying linkage outlier regions using alignment files and variants provided by PoolParty2. Finally, we apply our updated pipeline to an empirical dataset of low-coverage whole genomic data from population samples of Columbia River steelhead trout (Oncorhynchus mykiss), results from which demonstrate the genomic impacts of decades of artificial selection in a prominent hatchery stock. Thus, we not only demonstrate the utility of PoolParty2 for genomic studies that combine sequencing data from multiple individuals, but also illustrate how it compliments other bioinformatics resources such as angsd.

The genetic legacy of African Americans from Catoctin Furnace.

Few African Americans have been able to trace family lineages back to ancestors who died before the 1870 United States Census, the first in which all Black people were listed by name. We analyzed 27 individuals from Maryland's Catoctin Furnace African American Cemetery (1774-1850), identifying 41,799 genetic relatives among consenting research participants in 23andMe, Inc.'s genetic database. One of the highest concentrations of close relatives is in Maryland, suggesting that descendants of the Catoctin individuals remain in the area. We find that many of the Catoctin individuals derived African ancestry from the Wolof or Kongo groups and European ancestry from Great Britain and Ireland. This study demonstrates the power of joint analysis of historical DNA and large datasets generated through direct-to-consumer ancestry testing.

Landscape genomics of the streamside salamander: Implications for species management in the face of environmental change.

Understanding spatial patterns of genetic differentiation and local adaptation is critical in a period of rapid environmental change. Climate change and anthropogenic development have led to population declines and shifting geographic distributions in numerous species. The streamside salamander, Ambystoma barbouri, is an endemic amphibian with a small geographic range that predominantly inhabits small, ephemeral streams. As A. barbouri is listed as near-threatened by the IUCN, we describe range-wide patterns of genetic differentiation and adaptation to assess the species' potential to respond to environmental change. We use outlier scans and genetic-environment association analyses to identify genomic variation putatively underlying local adaptation across the species' geographic range. We find evidence for adaptation with a polygenic architecture and a set of candidate SNPs that identify genes putatively contributing to local adaptation. Our results build on earlier work that suggests that some A. barbouri populations are locally adapted despite evidence for asymmetric gene flow between the range core and periphery. Taken together, the body of work describing the evolutionary genetics of range limits in A. barbouri suggests that the species may be unlikely to respond naturally to environmental challenges through a range shift or in situ adaptation. We suggest that management efforts such as assisted migration may be necessary in future.

Utility of pooled sequencing for association mapping in nonmodel organisms.

High-density genome-wide sequencing increases the likelihood of discovering genes of major effect and genomic structural variation in organisms. While there is an increasing availability of reference genomes across broad taxa, the greatest limitation to whole-genome sequencing of multiple individuals continues to be the costs associated with sequencing. To alleviate excessive costs, pooling multiple individuals with similar phenotypes and sequencing the homogenized DNA (Pool-Seq) can achieve high genome coverage, but at the loss of individual genotypes. Although Pool-Seq has been an effective method for association mapping in model organisms, it has not been frequently utilized in natural populations. To extend bioinformatic tools for rapid implementation of Pool-Seq data in nonmodel organisms, we developed a pipeline called PoolParty and illustrate its effectiveness in genetic association mapping. Alignment expectations based on five pooled Chinook salmon (Oncorhynchus tshawytscha) libraries showed that approximately 48% genome coverage per library could be achieved with reasonable sequencing effort. We additionally examined male and female O. tshawytscha libraries to illustrate how Pool-Seq techniques can successfully map known genes associated with functional differences among sexes such as growth hormone 2. Finally, we compared pools of individuals of different spawning ages for each sex to discover novel genes involved with age at maturity in O. tshawytscha such as opsin4 and transmembrane protein19. While not appropriate for every system, Pool-Seq data processed by the PoolParty pipeline is a practical method for identifying genes of major effect in nonmodel organisms when high genome coverage is necessary and cost is a limiting factor.

How spatio-temporal habitat connectivity affects amphibian genetic structure.

Heterogeneous landscapes and fluctuating environmental conditions can affect species dispersal, population genetics, and genetic structure, yet understanding how biotic and abiotic factors affect population dynamics in a fluctuating environment is critical for species management. We evaluated how spatio-temporal habitat connectivity influences dispersal and genetic structure in a population of boreal chorus frogs (Pseudacris maculata) using a landscape genetics approach. We developed gravity models to assess the contribution of various factors to the observed genetic distance as a measure of functional connectivity. We selected (a) wetland (within-site) and (b) landscape matrix (between-site) characteristics; and (c) wetland connectivity metrics using a unique methodology. Specifically, we developed three networks that quantify wetland connectivity based on: (i) P. maculata dispersal ability, (ii) temporal variation in wetland quality, and (iii) contribution of wetland stepping-stones to frog dispersal. We examined 18 wetlands in Colorado, and quantified 12 microsatellite loci from 322 individual frogs. We found that genetic connectivity was related to topographic complexity, within- and between-wetland differences in moisture, and wetland functional connectivity as contributed by stepping-stone wetlands. Our results highlight the role that dynamic environmental factors have on dispersal-limited species and illustrate how complex asynchronous interactions contribute to the structure of spatially-explicit metapopulations.

Adaptive color polymorphism and unusually high local genetic diversity in the side-blotched lizard, Uta stansburiana.

Recently, studies of adaptive color variation have become popular as models for examining the genetics of natural selection. We examined color pattern polymorphism and genetic variation in a population of side-blotched lizards (Uta stansburiana) that is found in habitats with both dark (lava) and light colored (granite) substrates. We conducted a limited experiment for adult phenotypic plasticity in laboratory conditions. We recorded both substrate and lizard color patterns in the field to determine whether lizards tended to match their substrate. Finally we examined genetic variation in a gene (melanocortin 1 receptor) that has been shown to affect lizard color in other species and in a presumably neutral gene (mitochondrial cytochrome b). Populations were sampled in the immediate area of the lava flows as well as from a more distant site to examine the role of population structure. Our captive Uta did not change color to match their background. We show that side-blotched lizards tend to match the substrate on which it was caught in the field and that variation in the melanocortin 1 receptor gene does not correlate well with color pattern in this population. Perhaps the most remarkable result is that this population of side-blotched lizards shows extremely high levels of variation at both genetic markers, in the sense of allele numbers, with relatively low levels of between-allele sequence variation. Genetic variation across this small region was as great or greater than that seen in samples of pelagic fish species collected worldwide. Statistical analysis of genetic variation suggests rapid population expansion may be responsible for the high levels of variation.