Genomic diversity,population structure,and migration following rapid range expansion in the Balsam Poplar,Populus balsamifera

Rapid range expansions can cause pervasive changes in the genetic diversity and structure of populations. The postglacial history of the Balsam Poplar, Populus balsamifera, involved the colonization of most of northern North America, an area largely covered by continental ice sheets during the last glacial maximum. To characterize how this expansion shaped genomic diversity within and among populations, we developed 412 SNP markers that we assayed for a range‐wide sample of 474 individuals sampled from 34 populations. We complemented the SNP data set with DNA sequence data from 11 nuclear loci from 94 individuals, and used coalescent analyses to estimate historical population size, demographic growth, and patterns of migration. Bayesian clustering identified three geographically separated demes found in the Northern, Central, and Eastern portions of the species’ range. These demes varied significantly in nucleotide diversity, the abundance of private polymorphisms, and population substructure. Most measures supported the Central deme as descended from the primary refuge of diversity. Both SNPs and sequence data suggested recent population growth, and coalescent analyses of historical migration suggested a massive expansion from the Centre to the North and East. Collectively, these data demonstrate the strong influence that range expansions exert on genomic diversity, both within local populations and across the range. Our results suggest that an in‐depth knowledge of nucleotide diversity following expansion requires sampling within multiple populations, and highlight the utility of combining insights from different data types in population genomic studies.     

Genetic signals of demographic expansion in Downy Woodpecker (Picoides pubescens) after the last North American glacial maximum

The glacial cycles of the Pleistocene have been recognized as important, large-scale historical processes that strongly influenced the demographic patterns and genetic structure of many species. Here we present evidence of a postglacial expansion for the Downy Woodpecker (Picoides pubescens), a common member of the forest bird communities in North America with a continental distribution. DNA sequences from the mitochondrial tRNA-Lys, and ATPase 6 and 8 genes, and microsatellite data from seven variable loci were combined with a species distribution model (SDM) to infer possible historical scenarios for this species after the last glacial maximum. Analyses of Downy Woodpeckers from 23 geographic areas suggested little differentiation, shallow genealogical relationships, and limited population structure across the species' range. Microsatellites, which have higher resolution and are able to detect recent differences, revealed two geographic groups where populations along the eastern edge of the Rocky Mountains (Montana, Utah, Colorado, and southern Alberta) were genetically isolated from the rest of the sampled populations. Mitochondrial DNA, an important marker to detect historical patterns, recovered only one group. However, populations in Idaho and southeast BC contained high haplotype diversity and, in general were characterized by the absence of the most common mtDNA haplotype. The SDM suggested several areas in the southern US as containing suitable Downy Woodpecker habitat during the LGM. The lack of considerable geographic structure and the starburst haplotype network, combined with several population genetic tests, suggest a scenario of demographic expansion during the last part of Pleistocene and early Holocene.     

Insights into the neutral and adaptive processes shaping the spatial distribution of genomic variation in the economically important Moroccan locust (Dociostaurus maroccanus)

Understanding the processes that shape neutral and adaptive genomic variation is a fundamental step to determine the demographic and evolutionary dynamics of pest species. Here, we use genomic data obtained via restriction site‐associated DNA sequencing to investigate the genetic structure of Moroccan locust (Dociostaurus maroccanus) populations from the westernmost portion of the species distribution (Iberian Peninsula and Canary Islands), infer demographic trends, and determine the role of neutral versus selective processes in shaping spatial patterns of genomic variation in this pest species of great economic importance. Our analyses showed that Iberian populations are characterized by high gene flow, whereas the highly isolated Canarian populations have experienced strong genetic drift and loss of genetic diversity. Historical demographic reconstructions revealed that all populations have passed through a substantial genetic bottleneck around the last glacial maximum (~21 ka BP) followed by a sharp demographic expansion at the onset of the Holocene, indicating increased effective population sizes during warm periods as expected from the thermophilic nature of the species. Genome scans and environmental association analyses identified several loci putatively under selection, suggesting that local adaptation processes in certain populations might not be impeded by widespread gene flow. Finally, all analyses showed few differences between outbreak and nonoutbreak populations. Integrated pest management practices should consider high population connectivity and the potential importance of local adaptation processes on population persistence.     

Highest genetic diversity at the northern range limit of the rare orchid Isotria medeoloides

Populations in previously glaciated regions are often genetically depauperate in comparison with populations at lower latitudes, due either to bottlenecks experienced in post-glacial colonization or to contemporary genetic drift in small, peripheral populations. Populations of the rare self-fertilizing North American orchid Isotria medeoloides are largest in the previously glaciated region near the northern range limit, allowing us to examine the role of historical versus contemporary processes in determining population genetic diversity and structure. If contemporary processes predominate, genetic diversity should increase with increasing census size. In contrast, if sequential bottlenecks associated with colonization are paramount, diversity should decrease with latitude and be relatively insensitive to census size. We genotyped 299 individuals from 20 populations at four variable microsatellite loci to contrast genetic diversity and structure for populations in previously glaciated regions versus previously unglaciated regions. Populations were highly inbred (F=0.95) and highly differentiated (R(ST)=0.485). Across all sampled populations, genetic diversity decreased and genetic differentiation increased with declining population size. Small southern populations were especially differentiated and genetically depauperate. In the glaciated part of the range, genetic diversity increased as populations approached the northern range limit, demonstrating the centrality of contemporary processes for this post-glacial colonist.     

Is Gene Flow Promoting the Reversal of Pleistocene Divergence in the Mountain Chickadee (Poecile gambeli)?

The Pleistocene glacial cycles left a genetic legacy on taxa throughout the world; however, the persistence of genetic lineages that diverged during these cycles is dependent upon levels of gene flow and introgression. The consequences of secondary contact among taxa may reveal new insights into the history of the Pleistocene’s genetic legacy. Here, we use phylogeographic methods, using 20 nuclear loci from regional populations, to infer the consequences of secondary contact following divergence in the Mountain Chickadee (Poecile gambeli). Analysis of nuclear data identified two geographically-structured genetic groups, largely concordant with results from a previous mitochondrial DNA (mtDNA) study. Additionally, the estimated multilocus divergence times indicate a Pleistocene divergence, and are highly concordant with mtDNA. The previous mtDNA study showed a paucity of sympatry between clades, while nuclear patterns of gene flow show highly varied patterns between populations. The observed pattern of gene flow, from coalescent-based analyses, indicates southern populations in both clades exhibit little gene flow within or between clades, while northern populations are experiencing higher gene flow within and between clades. If this pattern were to persist, it is possible the historical legacy of Pleistocene divergence may be preserved in the southern populations only, and the northern populations would become a genetically diverse hybrid species.     

High-density SNP genotyping detects homogeneity of Spanish and French Basques, and confirms their genomic distinctiveness from other European populations

A recent study reported that Basques do not constitute a genetically distinct population, and that Basques from Spanish and French provinces do not show significant genetic similarity. These conclusions disagree with numerous previous studies, and are not consistent with the historical and linguistic evidence that supports the distinctiveness of Basques. In order to further investigate this controversy, we have genotyped 83 Spanish Basque individuals and used these data to infer population structure based on more than 60,000 single nucleotide polymorphisms of several European populations. Here, we present the first high-throughput analysis including Basques from Spanish and French provinces, and show that all Basques constitute a homogeneous group that can be clearly differentiated from other European populations.     

Assessment of genetic diversity and population structure in cultured Australian Pacific oysters

The recent development of Pacific oyster (Crassostrea gigas) SNP genotyping arrays has allowed detailed characterisation of genetic diversity and population structure within and between oyster populations. It also raises the potential of harnessing genomic selection for genetic improvement in oyster breeding programmes. The aim of this study was to characterise a breeding population of Australian oysters through genotyping and analysis of 18 027 SNPs, followed by comparison with genotypes of oyster sampled from Europe and Asia. This revealed that the Australian populations had similar population diversity (H_E) to oysters from New Zealand, the British Isles, France and Japan. Population divergence was assessed using PCA of genetic distance and revealed that Australian oysters were distinct from all other populations tested. Australian Pacific oysters originate from planned introductions sourced from three Japanese populations. Approximately 95% of these introductions were from geographically, and potentially genetically, distinct populations from the Nagasaki oysters assessed in this study. Finally, in preparation for the application of genomic selection in oyster breeding programmes, the strength of LD was evaluated and subsets of loci were tested for their ability to accurately infer relationships. Weak LD was observed on average; however, SNP subsets were shown to accurately reconstitute a genomic relationship matrix constructed using all loci. This suggests that low‐density SNP panels may have utility in the Australian population tested, and the findings represent an important first step towards the design and implementation of genomic approaches for applied breeding in Pacific oysters.     

林麝微卫星座位的分离和鉴定

利用改进的磁珠富集法从林麝(Moschus berezovskii)的基因组中分离到10个新的(AC)n重复微卫星座位,并分析其在24个林麝个体(取样于中国四川金凤山群体)中的多态性。其中,5个位点在24个林麝个体中具有4—13个等位基因,观察杂合度和期望杂合度分别是0.429—0.957和0.587—0.902,平均多态信息含量是0.730。表明这5个微卫星位点具有高度的多态性,可以用于林麝遗传多样性研究,对林麝的保护具有很重要的意义。     

Spatiotemporal analysis of gene flow in Chesapeake Bay Diamondback Terrapins (Malaclemys terrapin)

There is widespread concern regarding the impacts of anthropogenic activities on connectivity among populations of plants and animals, and understanding how contemporary and historical processes shape metapopulation dynamics is crucial for setting appropriate conservation targets. We used genetic data to identify population clusters and quantify gene flow over historical and contemporary time frames in the Diamondback Terrapin (Malaclemys terrapin). This species has a long and complicated history with humans, including commercial overharvesting and subsequent translocation events during the early twentieth century. Today, terrapins face threats from habitat loss and mortality in fisheries bycatch. To evaluate population structure and gene flow among Diamondback Terrapin populations in the Chesapeake Bay region, we sampled 617 individuals from 15 localities and screened individuals at 12 polymorphic microsatellite loci. Our goals were to demarcate metapopulation structure, quantify genetic diversity, estimate effective population sizes, and document temporal changes in gene flow. We found that terrapins in the Chesapeake Bay region harbour high levels of genetic diversity and form four populations. Effective population sizes were variable. Among most population comparisons, estimates of historical and contemporary terrapin gene flow were generally low (m ≈ 0.01). However, we detected a substantial increase in contemporary gene flow into Chesapeake Bay from populations outside the bay, as well as between two populations within Chesapeake Bay, possibly as a consequence of translocations during the early twentieth century. Our study shows that inferences across multiple time scales are needed to evaluate population connectivity, especially as recent changes may identify threats to population persistence.     

Contrasting levels of genetic differentiation among populations of wolverines (<Emphasis Type="Italic">Gulo gulo</Emphasis>) from northern Canada revealed by nuclear and mitochondrial loci

Habitat loss, fragmentation, overharvest, and other anthropogenic factors have resulted in population and distribution declines for North American wolverines (Gulo gulo). Currently, wolverines east of the Hudson Bay are endangered and possibly extinct, whereas the status of wolverines throughout the remaining Holarctic is vulnerable. Three previous studies using nuclear loci have detected little to no significant structuring among wolverines sampled across northern Canada. Based on these results it has been suggested that wolverines in northern Canada represent a single, panmictic population. However, as has been shown in numerous studies, in cases of female site fidelity, it is possible to have demographically autonomous populations even with male-biased gene flow. To better assess the genetic structure of wolverines in northern Canada, we examined nine microsatellite loci and DNA sequence variation from a 200 bp fragment of the mitochondrial (mtDNA) control region for 270 wolverines from nine collecting areas representing three regions of northern Canada. In agreement with previous studies, microsatellite analyses revealed a lack of significant population substructure (F _ST=0.0004). However, analysis of molecular variance, comparisons of pairwise F _ST values and nested-clade analysis of the mtDNA data revealed considerable genetic structuring among samples of wolverines from these three regions of northern Canada. These mitochondrial data provide evidence that wolverines in Canada are genetically structured due to female philopatry. The contrasting patterns of genetic differentiation based on nuclear and mitochondrial data highlight the importance of examining both nuclear and mitochondrial loci when attempting to elucidate patterns of genetic structure.     

Isolation and characterization of 13 microsatellite loci in Rhinolophus pusillus (least horseshoe bat) with cross-amplification in five related species

We used the enriched genomic library method to isolate and characterize dinucleotide microsatellite loci in the least horseshoe bat, Rhinolophus pusillus. Seventeen loci were obtained and tested on 31 individuals sampled from Guangxi Province in southern China. Thirteen of these markers were polymorphic with expected heterozygosity ranging from 0.821 to 0.909. A total of 164 alleles were detected and the number of alleles per locus ranged from 9 to 16 (mean 12.6). These polymorphic markers will be used to assess population structure in R. pusillus. In addition, successful cross-amplification in five congeneric bat species suggests most of these markers will also be useful for studying related species.     

Genetic similarities within and between human populations

The proportion of human genetic variation due to differences between populations is modest, and individuals from different populations can be genetically more similar than individuals from the same population. Yet sufficient genetic data can permit accurate classification of individuals into populations. Both findings can be obtained from the same data set, using the same number of polymorphic loci. This article explains why. Our analysis focuses on the frequency, omega, with which a pair of random individuals from two different populations is genetically more similar than a pair of individuals randomly selected from any single population. We compare omega to the error rates of several classification methods, using data sets that vary in number of loci, average allele frequency, populations sampled, and polymorphism ascertainment strategy. We demonstrate that classification methods achieve higher discriminatory power than omega because of their use of aggregate properties of populations. The number of loci analyzed is the most critical variable: with 100 polymorphisms, accurate classification is possible, but omega remains sizable, even when using populations as distinct as sub-Saharan Africans and Europeans. Phenotypes controlled by a dozen or fewer loci can therefore be expected to show substantial overlap between human populations. This provides empirical justification for caution when using population labels in biomedical settings, with broad implications for personalized medicine, pharmacogenetics, and the meaning of race.     

Multi-locus phylogeographic and population genetic analysis of Anolis carolinensis: historical demography of a genomic model species

The green anole (Anolis carolinensis) has been widely used as an animal model in physiology and neurobiology but has recently emerged as an important genomic model. The recent sequencing of its genome has shed new light on the evolution of vertebrate genomes and on the process that govern species diversification. Surprisingly, the patterns of genetic diversity within natural populations of this widespread and abundant North American lizard remain relatively unknown. In the present study, we use 10 novel nuclear DNA sequence loci (N?=?62 to 152) and one mitochondrial locus (N?=?226) to delimit green anole populations and infer their historical demography. We uncovered four evolutionarily distinct and geographically restricted lineages of green anoles using phylogenetics, bayesian clustering, and genetic distance methods. Molecular dating indicates that these lineages last shared a common ancestor ～2 million years ago. Summary statistics and analysis of the frequency distributions of DNA polymorphisms strongly suggest range-wide expansions in population size. Using Bayesian Skyline Plots, we inferred the timing of population size expansions, which differ across lineages, and found evidence for a relatively recent and rapid westward expansion of green anoles across the Gulf Coastal Plain during the mid-Pleistocene. One surprising result is that the distribution of genetic diversity is not consistent with a latitudinal shift caused by climatic oscillations as is observed for many co-distributed taxa. This suggests that the most recent Pleistocene glacial cycles had a limited impact on the geographic distribution of the green anole at the northern limits of its range.     

Conservation Genetics of Jacquemontia reclinata (Convolvulaceae), an Endangered Species from Southern Florida: Implications for Restoration Management

Guidelines designed to aid in the restoration of rare species have been previously proposed using two primary strategies to select individuals for augmentation and reintroduction: mixing progeny from different populations or separating individuals from different populations. Understanding the genetic structure and diversity of an endangered species can offer insights into conservation management strategies. We used random amplified polymorphic DNA markers to assess the genetic structure and diversity of Jacquemontia reclinata , a federally endangered species endemic to Southeastern Florida. We sampled 20 percent of total number of individuals from eight of the ten known wild populations. Across individuals high levels of polymorphic loci (94.7%) were found and larger populations had greater genetic diversity. Cluster and ordination analyses found that one population was genetically differentiated from all the others; this population grows in a unique habitat. Most genetic variation (77.5%) was found within populations, and genetic distances between populations were not explained by their geographic distances. We recommend the use of two management units in restoration programs for J. reclinata , one consisting of the genetically differentiated population and the second consisting of the other seven populations sampled.     

ABC inference of multi‐population divergence with admixture from unphased population genomic data

Rapidly developing sequencing technologies and declining costs have made it possible to collect genome‐scale data from population‐level samples in nonmodel systems. Inferential tools for historical demography given these data sets are, at present, underdeveloped. In particular, approximate Bayesian computation (ABC) has yet to be widely embraced by researchers generating these data. Here, we demonstrate the promise of ABC for analysis of the large data sets that are now attainable from nonmodel taxa through current genomic sequencing technologies. We develop and test an ABC framework for model selection and parameter estimation, given histories of three‐population divergence with admixture. We then explore different sampling regimes to illustrate how sampling more loci, longer loci or more individuals affects the quality of model selection and parameter estimation in this ABC framework. Our results show that inferences improved substantially with increases in the number and/or length of sequenced loci, while less benefit was gained by sampling large numbers of individuals. Optimal sampling strategies given our inferential models included at least 2000 loci, each approximately 2 kb in length, sampled from five diploid individuals per population, although specific strategies are model and question dependent. We tested our ABC approach through simulation‐based cross‐validations and illustrate its application using previously analysed data from the oak gall wasp, Biorhiza pallida.     

Evidence of genetic distinction and long-term population decline in wolves (Canis lupus) in the Italian Apennines

Historical information suggests the occurrence of an extensive human-caused contraction in the distribution range of wolves (Canis lupus) during the last few centuries in Europe. Wolves disappeared from the Alps in the 1920s, and thereafter continued to decline in peninsular Italy until the 1970s, when approximately 100 individuals survived, isolated in the central Apennines. In this study we performed a coalescent analysis of multilocus DNA markers to infer patterns and timing of historical population changes in wolves surviving in the Apennines. This population showed a unique mitochondrial DNA control-region haplotype, the absence of private alleles and lower heterozygosity at microsatellite loci, as compared to other wolf populations. Multivariate, clustering and Bayesian assignment procedures consistently assigned all the wolf genotypes sampled in Italy to a single group, supporting their genetic distinction. Bottleneck tests showed evidences of population decline in the Italian wolves, but not in other populations. Results of a Bayesian coalescent model indicate that wolves in Italy underwent a 100- to 1000-fold population contraction over the past 2000-10,000 years. The population decline was stronger and longer in peninsular Italy than elsewhere in Europe, suggesting that wolves have apparently been genetically isolated for thousands of generations south of the Alps. Ice caps covering the Alps at the Last Glacial Maximum (c. 18,000 years before present), and the wide expansion of the Po River, which cut the alluvial plains throughout the Holocene, might have provided effective geographical barriers to wolf dispersal. More recently, the admixture of Alpine and Apennine wolf populations could have been prevented by deforestation, which was already widespread in the fifteenth century in northern Italy. This study suggests that, despite the high potential rates of dispersal and gene flow, local wolf populations may not have mixed for long periods of time.     

Evidence for recent population bottlenecks in northern spotted owls (<Emphasis Type="Italic">Strix occidentalis caurina</Emphasis>)

The northern spotted owl (Strix occidentalis caurina) is one of the most controversial threatened subspecies ever listed under the US Endangered Species Act. Despite protection of its remaining forest habitat, recent field studies show continued declines of northern spotted owls. One potential threat to northern spotted owls which has not yet been shown is loss of genetic variation from population bottlenecks. Bottlenecks can increase the probability of mating among related individuals, potentially causing inbreeding depression, and can decrease adaptive potential. Here we report evidence for recent bottlenecks in northern spotted owls using a large genetic dataset (352 individuals and 10 microsatellite loci). The signature of bottlenecks was strongest in the Washington Cascade Mountains, in agreement with field data. Our results provide independent evidence that northern spotted owls have recently declined, and suggest that loss of genetic variation is an emerging threat to the subspecies’ persistence. Reduced effective population size (N _e) shown here in addition to field evidence for demographic decline highlights the increasing vulnerability of this bird to extinction.     

Phylogeography of spruce beetles (Dendroctonus rufipennis Kirby) (Curculionidae: Scolytinae) in North America

Tree-feeding insects that are widespread in north temperate regions are excellent models for studying how past glaciations have impacted differentiation and speciation. We used mitochondrial DNA (mtDNA) sequences and allele frequencies at nine microsatellite loci to examine genetic population structure across the current range of the spruce beetle (Dendroctonus rufipennis), an economically important insect in North America. Two major haplotype groups occur across northern North America, from Newfoundland to Alaska, on white spruce (Picea glauca), and a third distinctive haplotype group occurs throughout the Rocky Mountains on Engelmann spruce (Picea engelmannii). The two mtDNA lineages found in northern populations are 3-4% divergent from each other and from the lineages found in the Rocky Mountains. Analyses of microsatellite data also suggest the existence of major population groupings associated with different geographical regions. In the Pacific Northwest, concordant contact zones for genetically distinct populations of spruce beetles and their principal hosts appear to reflect recent secondary contact. Although we could detect no evidence of historical mtDNA gene flow between allopatric population groups, patterns of variation in the Pacific Northwest suggest recent hybridization and introgression. Together with the pollen record for spruce, they also suggest that beetles have spread from at least three glacial refugia. A minimum estimate of divergence time between the Rocky Mountain and northern populations was 1.7 Myr (million years), presumably reflecting the combined effects of isolation during multiple glacial cycles.     

Isolation and characterization of microsatellite loci in Pratt’s Leaf-nosed Bat (<Emphasis Type="Italic">Hipposideros pratti</Emphasis>) and cross-species amplification in closely related taxa

We developed nine microsatellite loci using an enriched library method from the genomic DNA of Pratt’s leaf-nosed bat (Hipposideros pratti). These loci were tested on 96 individuals sampled from Sichuan Province, China. The number of alleles per locus varied from 4 to 14. The expected and observed heterozygosity values ranged from 0.513 to 0.886 and from 0.375 to 0.966, respectively. Three microsatellite loci departed significantly from Hardy–Weinberg equilibrium (HWE). No linkage disequilibrium was found. These microsatellite loci will be used in studies of conservation genetics in this species.