Heterozygosity in an isolated population of a large mammal founded by four individuals is predicted by an individual-based genetic model

Contrasting genomic consequences of anthropogenic reintroduction and natural recolonization in high-arctic wild reindeer

Anthropogenic reintroduction can supplement natural recolonization in reestablishing a species' distribution and abundance. However, both reintroductions and recolonizations can give rise to founder effects that reduce genetic diversity and increase inbreeding, potentially causing the accumulation of genetic load and reduced fitness. Most current populations of the endemic high-arctic Svalbard reindeer (Rangifer tarandus platyrhynchus) originate from recent reintroductions or recolonizations following regional extirpations due to past overharvesting. We investigated and compared the genomic consequences of these two paths to reestablishment using whole-genome shotgun sequencing of 100 Svalbard reindeer across their range. We found little admixture between reintroduced and natural populations. Two reintroduced populations, each founded by 12 individuals around four decades (i.e. 8 reindeer generations) ago, formed two distinct genetic clusters. Compared to the source population, these populations showed only small decreases in genome-wide heterozygosity and increases in inbreeding and lengths of runs of homozygosity. In contrast, the two naturally recolonized populations without admixture possessed much lower heterozygosity, higher inbreeding and longer runs of homozygosity, possibly caused by serial population founder effects and/or fewer or more genetically related founders than in the reintroduction events. Naturally recolonized populations can thus be more vulnerable to the accumulation of genetic load than reintroduced populations. This suggests that in some organisms even small-scale reintroduction programs based on genetically diverse source populations can be more effective than natural recolonization in establishing genetically diverse populations. These findings warrant particular attention in the conservation and management of populations and species threatened by habitat fragmentation and loss.

Density-habitat relationships of white-tailed deer (Odocoileus virginianus) in Finland

In heterogeneous landscapes, resource selection constitutes a crucial link between landscape and population-level processes such as density. We conducted a non-invasive genetic study of white-tailed deer in southern Finland in 2016 and 2017 using fecal DNA samples to understand factors influencing white-tailed deer density and space use in late summer prior to the hunting season. We estimated deer density as a function of landcover types using a spatial capture-recapture (SCR) model with individual identities established using microsatellite markers. The study revealed second-order habitat selection with highest deer densities in fields and mixed forest, and third-order habitat selection (detection probability) for transitional woodlands (clear-cuts) and closeness to fields. Including landscape heterogeneity improved model fit and increased inferred total density compared with models assuming a homogenous landscape. Our findings underline the importance of including habitat covariates when estimating density and exemplifies that resource selection can be studied using non-invasive methods.

Population growth lags in introduced species

When introduced to new ecosystems, species' populations often grow immediately postrelease. Some introduced species, however, maintain a low population size for years or decades before sudden, rapid population growth is observed. Because exponential population growth always starts slowly, it can be difficult to distinguish species experiencing the early phases of slow exponential population growth (inherent lags) from those with actively delayed growth rates (prolonged lags). Introduced ungulates provide an excellent system in which to examine lags, because some introduced ungulate populations have demonstrated rapid population growth immediately postintroduction, while others have not. Using studies from the literature, we investigated which exotic ungulate species and populations (n = 36) showed prolonged population growth lags by comparing the doubling time of real ungulate populations to those predicted from exponential growth models for theoretical populations. Having identified the specific populations that displayed prolonged lags, we examined the impacts of several environmental and biological variables likely to influence the length of lag period. We found that seventeen populations (47%) showed significant prolonged population growth lags. We could not, however, determine the specific factors that contributed to the length of these lag phases, suggesting that these ungulate populations' growth is idiosyncratic and difficult to predict. Introduced species that exhibit delayed growth should be closely monitored by managers, who must be proactive in controlling their growth to minimize the impact such populations may have on their environment.

Genetic signatures of population bottlenecks, relatedness, and inbreeding highlight recent and novel conservation concerns in the Egyptian vulture

The assessment of temporal variation in genetic features can be particularly informative on the factors behind demography and viability of wildlife populations and species. We used molecular methods to evaluate neutral genetic variation, relatedness, bottlenecks, and inbreeding in a declining population of Egyptian vulture (Neophron percnopterus) in central Spain. The results show that the genetic diversity remained relatively stable over a period of twelve years despite the decline in census and effective population sizes in the last decades. A relatively high proportion of nestlings from different and distant territories showed high relatedness in each study year. We also found support for an increasing impact of severe recent (contemporary) rather than distant (historical) past demographic bottlenecks, and the first evidence of inbred mating between full siblings coinciding with lethal malformations in offspring. The inbred nestling with feather malformations was positive to beak and feather disease virus recorded for the first time in this species. These results alert on recent and novel threats potentially affecting health and reducing the adaptive potential of individuals in this threatened species.

An isolated white-tailed deer (Odocoileus virginianus) population on St. John, US Virgin Islands shows low inbreeding and comparable heterozygosity to other larger populations

This is the first study to document the genetic diversity of the white-tailed deer population on St. John, US Virgin Islands. The island population was founded by a small number of animals, has very limited hunting or predation, and recently experienced a reduction in size following an extended drought in 2015. DNA samples were collected from hair from 23 anesthetized adult deer (13 males, 10 females) ranging in age from 1 to 8 years (3.36 ± 1.9 years) and also from fecal DNA samples, for a total of 42 individuals analyzed for genetic diversity. The St. John deer data set averaged 4.19 alleles per marker and demonstrates the second lowest number of alleles (A) when compared to other populations of Odocoileus virginianus (4.19). Heterozygosity was similar to the other studies (0.54) with little evidence of inbreeding. To explain the level of heterozygosity and level of inbreeding within the St. John population, three hypotheses are proposed, including the effect of intrinsic biological traits within the population, a recent infusion of highly heterogeneous loci from North American populations, and a consistent level of immigration from a nearby island. Additional work is needed to further understand the genetic history of the St. John and regional deer populations.

Fine-scale genetic structure in the critically endangered red-fronted macaw in the absence of geographic and ecological barriers

Behavioural and socio-cultural traits are recognized in the restriction of gene flow in species with high cognitive capacity and complex societies. This isolation by social barriers has been generally overlooked in threatened species by assuming disrupted gene flow due to population fragmentation and decline. We examine the genetic structure and ecology of the global population of the Critically Endangered red-fronted macaw (Ara rubrogenys), an endemic species to the inter-Andean valleys of Bolivia. We found a fine-scale genetic structuring in four genetic clusters. Genetic diversity was higher in wild compared to captive-bred macaws, but similar to that of captive wild-caught macaws. We found no clear evidence of severe genetic erosion in the population in recent decades, but it was patent in historic times, overlapping with drastic human habitat transformation and macaw persecution over millennia. We found no evidence of geographical and ecological barriers, owing to the high dispersal ability, nesting and foraging habits between genetic clusters. The lack of genetic intermixing despite long-distance foraging and seasonal movements suggests recruitment in natal colonies and other social factors reinforcing philopatry-related genetic structure. Conservation efforts should be specifically focussed on major threats in each genetic cluster as independent conservation units, and also considered in ex-situ management.

Phylogeny and biogeography of the Japanese rhinoceros beetle, Trypoxylus dichotomus (Coleoptera: Scarabaeidae) based on SNP markers

The Japanese rhinoceros beetle Trypoxylus dichotomus is one of the largest beetle species in the world and is commonly used in traditional Chinese medicine. Ten subspecies of T. dichotomus and a related Trypoxylus species (T. kanamorii) have been described throughout Asia, but their taxonomic delimitations remain problematic. To clarify issues such as taxonomy, and the degree of genetic differentiation of Trypoxylus populations, we investigated the genetic structure, genetic variability, and phylogeography of 53 specimens of Trypoxylus species from 44 locations in five Asian countries (China, Japan, Korea, Thailand, and Myanmar). Using specific-locus amplified fragment sequencing (SLAF-seq) techniques, we developed 330,799 SLAFs over 114.16M reads, in turn yielding 46,939 high-resolution single nucleotide polymorphisms (SNPs) for genotyping. Phylogenetic analysis of SNPs indicated the presence of three distinct genetic groups, suggesting that the various subspecies could be treated as three groups of populations. PCA and ADMIXTURE analysis also identified three genetic clusters (North, South, West), which corresponded to their locations, suggesting that geographic factors were important in maintaining within population homogeneity and between population divergence. Analyses of SNP data confirmed the monophyly of certain subspecies on islands, while other subspecies (e.g., T. d. septentrionalis) were found to be polyphyletic and nested in more than one lineage. AMOVA demonstrated high level of differentiation among populations/groups. Also, pairwise F ST values revealed high differentiation, particularly between South and West, as well as between North and South. Despite the differentiation, measurable gene flow was inferred between genetic clusters but at varying rates and directions. Our study demonstrated that SLAF-seq derived markers outperformed 16S and COII sequences and provided improved resolution of the genetic differentiation of rhinoceros beetle populations from a large part of the species' range.

Genetic diversity and demographic history of introduced sika deer on the Delmarva Peninsula

The introduction of non-native species can have long-term effects on native plant and animal communities. Introduced populations are occasionally not well understood and offer opportunities to evaluate changes in genetic structure through time and major population changes such as bottleneck and or founder events. Invasive species can often evolve rapidly in new and novel environments, which could be essential to their long-term success. Sika deer are native to East Asia, and their introduction and establishment to the Delmarva Peninsula, USA, is poorly documented, but probably involved ≥1 founder and/or bottleneck events. We quantified neutral genetic diversity in the introduced population and compared genetic differentiation and diversity to the presumed source population from Yakushima Island, Japan, and a captive population of sika deer in Harrington, Delaware, USA. Based on the data from 10 microsatellite DNA loci, we observed reduced genetic variation attributable to founder events, support for historic hybridization events, and evidence that the population did originate from Yakushima Island stocks. Estimates of population structure through Bayesian clustering and demographic history derived from approximate Bayesian computation (ABC), were consistent with the hypothesized founder history of the introduced population in both timing and effective population size (approximately five effective breeding individuals, an estimated 36 generations ago). Our ABC results further supported a single introduction into the wild happening before sika deer spread throughout the Delmarva. We conclude that free-ranging sika deer on Delmarva are descended from ca. five individuals introduced about 100 years ago from captive stocks of deer maintained in the United Kingdom. Free-ranging sika deer on Delmarva have lost neutral diversity due to founder and bottleneck events, yet populations have expanded in recent decades and show no evidence of abnormalities associated with inbreeding. We suggest management practices including increasing harvest areas and specifically managing sika deer outside of Maryland.

Genetic diversity and population structure of the Sapsaree, a native Korean dog breed

Background

The Sapsaree is a breed of dog (Canis familiaris) native to Korea, which became perilously close to extinction in the mid-1980s. However, with systematic genetic conservation and restoration efforts, this breed was rescued from extinction and population sizes have been gradually increasing over the past few decades. The aim of this study was to ascertain novel information about the genetic diversity, population structure, and demographic history of the Sapsaree breed using genome-wide single nucleotide polymorphism data. We characterized the genetic profile of the Sapsaree breed by comparison with seven foreign dog breeds with similar morphologies to estimate genetic differentiation within and among these breeds.

Results

The results suggest that Sapsarees have higher genetic variance compared with the other breeds analyzed. The majority of the Sapsarees in this study share a discrete genetic pattern, although some individuals were slightly different, possibly as a consequence of the recent restoration process. Concordant results from analyses of linkage disequilibrium, effective population size, genetic diversity, and population structural analyses illustrate a relationship among the Sapsaree and the Tibetan breeds Tibetan terrier and Lhasa Apso, and a small genetic introgression from European breeds. The effective population size of the Sapsaree has contracted dramatically over the past generations, and is currently insufficient to maintain long-term viability of the breed’s genetic diversity.

Conclusions

This study provides novel insights regarding the genetic diversity and population structure of the native Korean dog breed Sapsaree. Our results suggest the importance of a strategic and systematic approach to ensure the genetic diversity and the authenticity of the Sapsaree breed.

Electronic supplementary material

The online version of this article (10.1186/s12863-019-0757-5) contains supplementary material, which is available to authorized users.

Crowdsourcing-based nationwide tick collection reveals the distribution of Ixodes ricinus and I. persulcatus and associated pathogens in Finland

A national crowdsourcing-based tick collection campaign was organized in 2015 with the objective of producing novel data on tick distribution and tick-borne pathogens in Finland. Nearly 20 000 Ixodes ticks were collected. The collected material revealed the nationwide distribution of I. persulcatus for the first time and a shift northwards in the distribution of I. ricinus in Finland. A subset of 2038 tick samples containing both species was screened for Borrelia burgdorferi sensu lato (the prevalence was 14.2% for I. ricinus and 19.8% for I. persulcatus), B. miyamotoi (0.2% and 0.4%, respectively) and tick-borne encephalitis virus (TBEV; 0.2% and 3.0%, respectively). We also report new risk areas for TBEV in Finland and, for the first time, the presence of B. miyamotoi in ticks from mainland Finland. Most importantly, our study demonstrates the overwhelming power of citizen science in accomplishing a collection effort that would have been impossible with the scientific community alone.

Rethinking the influence of hydroelectric development on gene flow in a long-lived fish, the Lake Sturgeon Acipenser fulvescens

Many hydroelectric dams have been in place for 50 - >100 years, which for most fish species means that enough generations have passed for fragmentation induced divergence to have accumulated. However, for long-lived species such as Lake Sturgeon, Acipenser fulvescens, it should be possible to discriminate between historical population structuring and contemporary gene flow and improve the broader understanding of anthropogenic influence. On the Winnipeg River, Manitoba, two hypotheses were tested: 1) Measureable quantities of former reservoir dwelling Lake Sturgeon now reside downstream of the Slave Falls Generating Station, and 2) genetically differentiated populations of Lake Sturgeon occur upstream and downstream, a result of historical structuring. Genetic methods based on ten microsatellite markers were employed, and simulations were conducted to provide context. With regards to contemporary upstream to downstream contributions, the inclusion of length-at-age data proved informative. Both pairwise relatedness and Bayesian clustering analysis substantiated that fast-growing outliers, apparently entrained after residing in the upstream reservoir for several years, accounted for ~15% of the Lake Sturgeon 525–750 mm fork length captured downstream. With regards to historical structuring, upstream and downstream populations were found to be differentiated (F_ST = 0.011, and 0.013–0.014 when fast-growing outliers were excluded), and heterozygosity metrics were higher for downstream versus upstream juveniles. Historical asymmetric (downstream) gene flow in the vicinity of the generating station was the most logical explanation for the observed genetic structuring. In this section of the Winnipeg River, construction of a major dam does not appear to have fragmented a previously panmictic Lake Sturgeon population, but alterations to habitat may be influencing upstream to downstream contributions in unexpected ways.

Variance in disease risk: rural populations and genetic diversity

Over 19% of the US population resides in rural areas, where studies of disease risk and disease outcomes are difficult to assess due to smaller populations and lower incidence. While some studies suggest rural disparities for different chronic diseases, the data are inconsistent across geography and definitions of rurality. We reviewed the literature to examine if local variations in population genomic diversity may plausibly explain inconsistencies in estimating disease risk. Many rural communities were founded over 150 years ago by small groups of ethnically and ancestrally similar families. These have since endured relative geographical isolation, similar to groups in other industrialized nations, perhaps resulting in founder effects impacting local disease susceptibility. Studies in Europe and Asia have found that observably different phenotypes may appear in isolated communities within 100 years, and that genomic variation can significantly vary over small geographical scales. Epidemiological studies utilizing common “rural” definitions may miss significant disease differences due to assumptions of risk homogeneity and misinterpretation of administrative definitions of rurality. Local genomic heterogeneity should be an important aspect of chronic disease epidemiology in rural areas, and it is important to consider for designing studies and interpreting results, enabling a better understanding of the heritable components of complex diseases.

MtDNA and Y-chromosomal diversity in the Chachapoya, a population from the northeast Peruvian Andes-Amazon divide

OBJECTIVES

The ancient Chachapoya were an aggregate of several ethnic groups that shared a common language, religion, and material culture. They inhabited a territory at the juncture of the Andes and the Amazon basin. Their position between those ecozones most likely influenced their genetic composition. We attempted to better understand their population history by assessing the contemporary genetic diversity in the Chachapoya and three of their immediate neighbors (Huancas, Jivaro, and Cajamarca). We inferred signatures of demographic history and genetic affinities, and contrasted the findings with data from other populations on local and continental scales.

METHODS

We studied mitochondrial DNA (mtDNA; hypervariable segment [HVSI and HVSII]) and Y chromosome (23 short tandem repeats (STRs)) marker data in 382 modern individuals. We used Sanger sequencing for mtDNA and a commercially available kit for Y-chromosomal STR typing.

RESULTS

The Chachapoya had affinities with various populations of Andean and Amazonian origin. When examining the Native American component, the Chachapoya displayed high levels of genetic diversity. Together with other parameters, for example, large Tajima's D and Fu's Fs, the data indicated no drastic reduction of the population size in the past.

CONCLUSION

The high level of diversity in the Chachapoya, the lack of evidence of drift in the past, and genetic affinities with a broad range of populations in the Americas reflects an intricate population history in the region. The new genetic data from the Chachapoya indeed seems to point to a genetic complexity that is not yet resolved but beginning to be elucidated. Am. J. Hum. Biol. 28:857-867, 2016. © 2016Wiley Periodicals, Inc.

Serological evidence of tick-borne encephalitis virus infection in moose and deer in Finland: sentinels for virus circulation

Background

The incidence of tick-borne encephalitis (TBE) in humans has increased in Finland, and the disease has emerged in new foci. These foci have been investigated to determine the circulating virus subtype, the tick host species and the ecological parameters, but countrywide epidemiological information on the distribution of TBEV has been limited.

Methods

In this study, we screened sera from hunter-harvested wild cervids for the presence of antibodies against tick-borne encephalitis virus (TBEV) with a hemagglutination inhibition test. The positive results were confirmed by a neutralisation assay.

Results

Nine (0.74 %) of 1213 moose, one (0.74 %) of 135 white-tailed deer, and none of the 17 roe deer were found seropositive for TBEV. A close geographical congruence between seropositive cervids and recently reported human TBE cases was observed: nine of the ten seropositive animals were from known endemic areas.

Conclusions

Our results confirm the local circulation of TBEV in several known endemic areas. One seropositive moose had been shot in an area where human TBE cases have not been reported, suggesting a possible new focus. Moose appear to be a useful sentinel animal for the presence of TBEV in the taiga region.

Using heterozygosity-fitness correlations to study inbreeding depression in an isolated population of white-tailed deer founded by few individuals

A heterozygosity–fitness correlations (HFCs) may reflect inbreeding depression, but the extent to which they do so is debated. HFCs are particularly likely to occur after demographic disturbances such as population bottleneck or admixture. We here study HFC in an introduced and isolated ungulate population of white-tailed deer Odocoileus virginianus in Finland founded in 1934 by four individuals. A total of 422 ≥ 1-year-old white-tailed deer were collected in the 2012 hunting season in southern Finland and genotyped for 14 microsatellite loci. We find significant identity disequilibrium as estimated by g₂. Heterozygosity was positively associated with size- and age-corrected body mass, but not with jaw size or (in males) antler score. Because of the relatively high identity disequilibrium, heterozygosity of the marker panel explained 51% of variation in inbreeding. Inbreeding explained approximately 4% of the variation in body mass and is thus a minor, although significant source of variation in body mass in this population. The study of HFC is attractive for game- and conservation-oriented wildlife management because it presents an affordable and readily used approach for genetic monitoring that allowing identification of fitness costs associated with genetic substructuring in what may seem like a homogeneous population.

Demographic histories and genetic diversities of Fennoscandian marine and landlocked ringed seal subspecies

Island populations are on average smaller, genetically less diverse, and at a higher risk to go extinct than mainland populations. Low genetic diversity may elevate extinction probability, but the genetic component of the risk can be affected by the mode of diversity loss, which, in turn, is connected to the demographic history of the population. Here, we examined the history of genetic erosion in three Fennoscandian ringed seal subspecies, of which one inhabits the Baltic Sea 'mainland' and two the 'aquatic islands' composed of Lake Saimaa in Finland and Lake Ladoga in Russia. Both lakes were colonized by marine seals after their formation c. 9500 years ago, but Lake Ladoga is larger and more contiguous than Lake Saimaa. All three populations suffered dramatic declines during the 20th century, but the bottleneck was particularly severe in Lake Saimaa. Data from 17 microsatellite loci and mitochondrial control-region sequences show that Saimaa ringed seals have lost most of the genetic diversity present in their Baltic ancestors, while the Ladoga population has experienced only minor reductions. Using Approximate Bayesian computing analyses, we show that the genetic uniformity of the Saimaa subspecies derives from an extended founder event and subsequent slow erosion, rather than from the recent bottleneck. This suggests that the population has persisted for nearly 10,000 years despite having low genetic variation. The relatively high diversity of the Ladoga population appears to result from a high number of initial colonizers and a high post-colonization population size, but possibly also by a shorter isolation period and/or occasional gene flow from the Baltic Sea.

Causes and consequences of fine-scale population structure in a critically endangered freshwater seal

Background

Small, genetically uniform populations may face an elevated risk of extinction due to reduced environmental adaptability and individual fitness. Fragmentation can intensify these genetic adversities and, therefore, dispersal and gene flow among subpopulations within an isolated population is often essential for maintaining its viability. Using microsatellite and mtDNA data, we examined genetic diversity, spatial differentiation, interregional gene flow, and effective population sizes in the critically endangered Saimaa ringed seal (Phoca hispida saimensis), which is endemic to the large but highly fragmented Lake Saimaa in southeastern Finland.

Results

Microsatellite diversity within the subspecies (H_E = 0.36) ranks among the lowest thus far recorded within the order Pinnipedia, with signs of ongoing loss of individual heterozygosity, reflecting very low effective subpopulation sizes. Bayesian assignment analyses of the microsatellite data revealed clear genetic differentiation among the main breeding areas, but interregional structuring was substantially weaker in biparentally inherited microsatellites (F_ST = 0.107) than in maternally inherited mtDNA (F_ST = 0.444), indicating a sevenfold difference in the gene flow mediated by males versus females.

Conclusions

Genetic structuring in the population appears to arise from the joint effects of multiple factors, including small effective subpopulation sizes, a fragmented lacustrine habitat, and behavioural dispersal limitation. The fine-scale differentiation found in the landlocked Saimaa ringed seal is especially surprising when contrasted with marine ringed seals, which often exhibit near-panmixia among subpopulations separated by hundreds or even thousands of kilometres. Our results demonstrate that population structures of endangered animals cannot be predicted based on data on even closely related species or subspecies.