An official website of the United States government
PremiseThe ability to sequence genome‐scale data from herbarium specimens would allow for the economical development of data sets with broad taxonomic and geographic sampling that would otherwise not be possible. Here, we evaluate the utility of a basic double‐digest restriction site–associatedDNAsequencing (ddRADseq) protocol usingDNAs from four genera extracted from both silica‐dried and herbarium tissue. MethodsDNAs fromDraba,Boechera,Solidago, andIlexwere processed with a ddRADseq protocol. The effects ofDNAdegradation, taxon, and specimen age were assessed. ResultsAlthough taxon, preservation method, and specimen age affected data recovery, large phylogenetically informative data sets were obtained from the majority of samples. DiscussionThese results suggest that herbarium samples can be incorporated into ddRADseq project designs, and that specimen age can be used as a rapid on‐site guide for sample choice. The detailed protocol we provide will allow users to pursue herbarium‐based ddRADseq projects that minimize the expenses associated with fieldwork and sample evaluation.
more »
« less
Population genomics of divergence among extreme and intermediate color forms in a polymorphic insect
Abstract Geographic variation in insect coloration is among the most intriguing examples of rapid phenotypic evolution and provides opportunities to study mechanisms of phenotypic change and diversification in closely related lineages. The bumble beeBombus bifariuscomprises two geographically disparate color groups characterized by red‐banded and black‐banded abdominal pigmentation, but with a range of spatially and phenotypically intermediate populations across western North America. Microsatellite analyses have revealed thatB. bifariusin the USA are structured into two major groups concordant with geography and color pattern, but also suggest ongoing gene flow among regional populations. In this study, we better resolve the relationships among major color groups to better understand evolutionary mechanisms promoting and maintaining such polymorphism. We analyze >90,000 and >25,000 single‐nucleotide polymorphisms derived from transcriptome (RNAseq) and double digest restriction site associatedDNAsequencing (ddRAD), respectively, in representative samples from spatial and color pattern extremes inB. bifariusas well as phenotypic and geographic intermediates. Both ddRADandRNAseq data illustrate substantial genome‐wide differentiation of the red‐banded (eastern) color form from both black‐banded (western) and intermediate (central) phenotypes and negligible differentiation among the latter populations, with no obvious admixture among bees from the two major lineages. Results thus indicate much stronger background differentiation amongB. bifariuslineages than expected, highlighting potential challenges for revealing loci underlying color polymorphism from population genetic data alone. These findings will have significance for resolving taxonomic confusion in this species and in future efforts to investigate color‐pattern evolution inB. bifariusand other polymorphic bumble bee species.
more »
« less
- Award ID(s):
- 1457659
- PAR ID:
- 10196916
- Publisher / Repository:
- Wiley Blackwell (John Wiley & Sons)
- Date Published:
- Journal Name:
- Ecology and Evolution
- Volume:
- 6
- Issue:
- 4
- ISSN:
- 2045-7758
- Page Range / eLocation ID:
- p. 1075-1091
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract Identifying drivers of dispersal limitation and genetic differentiation is a key goal in biogeography. We examine patterns of population connectivity and genetic diversity using restriction site‐associatedDNAsequencing (RADseq) in two bumble bee species,Bombus vosnesenskiiandBombus bifarius,across latitude and altitude in mountain ranges from California, Oregon and Washington, U.S.A.Bombus vosnesenskii, which occurs across a broader elevational range at most latitudes, exhibits little population structure whileB. bifarius, which occupies a relatively narrow higher elevation niche across most latitudes, exhibits much stronger population differentiation, although gene flow in both species is best explained by isolation with environmental niche resistance. A relationship between elevational habitat breadth and genetic diversity is also apparent, withB. vosnesenskiiexhibiting relatively consistent levels of genetic diversity across its range, whileB. bifariushas reduced genetic diversity at low latitudes, where it is restricted to high‐elevation habitat. The results of this study highlight the importance of the intersect between elevational range and habitat suitability in influencing population connectivity and suggest that future climate warming will have a fragmenting effect even on populations that are presently well connected, as they track their thermal niches upward in montane systems.more » « less
-
Abstract Next‐generation sequencing technologies now allow researchers of non‐model systems to perform genome‐based studies without the requirement of a (often unavailable) closely related genomic reference. We evaluated the role of restriction endonuclease (RE) selection in double‐digest restriction‐site‐associatedDNAsequencing (ddRADseq) by generating reduced representation genome‐wide data using four differentREcombinations. Our expectation was thatREselections targeting longer, more complex restriction sites would recover fewer loci thanREwith shorter, less complex sites. We sequenced a diverse sample of non‐model arachnids, including five congeneric pairs of harvestmen (Opiliones) and four pairs of spiders (Araneae). Sample pairs consisted of either conspecifics or closely related congeneric taxa, and in total 26 sample pair analyses were tested. Sequence demultiplexing, read clustering and variant calling were performed in thepyRADprogram. The 6‐base pair cutterEcoRIcombined with methylated site‐specific 4‐base pair cutterMspIproduced, on average, the greatest numbers of intra‐individual loci and shared loci per sample pair. As expected, the number of shared loci recovered for a sample pair covaried with the degree of genetic divergence, estimated with cytochrome oxidase I sequences, although this relationship was non‐linear. Our comparative results will prove useful in guiding protocol selection for ddRADseq experiments on many arachnid taxa where reference genomes, even from closely related species, are unavailable.more » « less
-
ABSTRACT Sister taxa that have diverged and persisted in sympatry have likely been exposed to the same general environmental changes throughout their evolutionary history and may thus exhibit similar phylogeographies. Here, we compare the phylogeographic patterns of two sister species of isopods (genusTylos) that have broadly overlapping distributions but distinct habitat preferences in the supralittoral zone of Chile. The dynamic geoclimatic history of this region during the Quaternary has been implicated in shaping the evolutionary histories of other coastal taxa.Tylos spinulosusis found in sandy beaches at latitudes ~27°–30° S, whereasTylos chilensishas been found in rocky shores at ~27°–33° S and at ~39°–42° S. We sampled both species across their ranges (collectively from 20 localities) and obtained sequences from at least one mitochondrial gene for 95 T. chilensisand 41 T. spinulosus. We used phylogenetics and population genetics methods to analyze four single‐gene and one concatenated datasets: 12S rDNA (n = 130); 16S rDNA (n = 31); Cytochrome oxidase subunit I (n = 28); Cytochrome b (n = 24); concatenation of the four genes (n = 24). Both species show high levels of isolation of local populations, consistent with expectations from their limited autonomous dispersal potential. However, they exhibit strikingly different mitochondrial phylogeographic patterns.Tylos chilensisshows evidence of multiple relatively deep divergence events leading to geographically restricted lineages that appear to have persisted over multiple glaciations. Surprisingly, one lineage ofT. chilensiswas found in geographically distant localities, suggesting the possibility of human‐mediated dispersal.Tylos spinulosusappears to have undergone a relatively recent bottleneck followed by a population/range expansion. Differences in life histories and habitat preferences or stochasticity may have contributed to these striking phylogeographic differences. Finally, the high levels of differentiation and isolation among populations indicate that they are highly vulnerable to extirpation. We discuss threats to their persistence and recommendations for their conservation.more » « less
-
ABSTRACT Coalescent modelling of hybrid zones can provide novel insights into the historical demography of populations, including divergence times, population sizes, introgression proportions, migration rates and the timing of hybrid zone formation. We used coalescent analysis to determine whether the hybrid zone between phylogeographic lineages of the Plateau Fence Lizard (Sceloporus tristichus) in Arizona formed recently due to human‐induced landscape changes, or if it originated during Pleistocene climatic shifts. Given the presence of mitochondrial DNA from another species in the hybrid zone (Southwestern Fence Lizard,S. cowlesi), we tested for the presence ofS. cowlesinuclear DNA in the hybrid zone as well as reassessed the species boundary betweenS. tristichusandS. cowlesi. No evidence ofS. cowlesinuclear DNA is found in the hybrid zone, and the paraphyly of both species raises concerns about their taxonomic validity. Introgression analysis placed the divergence time between the parental hybrid zone populations at approximately 140 kya and their secondary contact and hybridization at approximately 11 kya at the end of the Pleistocene. Introgression proportions estimated for hybrid populations are correlated with their geographic distance from parental populations. The multispecies coalescent with migration provided significant support for unidirectional migration moving from south to north, which is consistent with spatial cline analyses that suggest a slow but steady northward shift of the centre of the hybrid zone over the last two decades. When analysing hybrid populations sampled along a linear transect, coalescent methods can provide novel insights into hybrid zone dynamics.more » « less
