Bottleneck effects on the sika deer Cervus nippon population in Hokkaido, revealed by ancient DNA analysis

A historic religious sanctuary may have preserved ancestral genetics of Japanese sika deer (Cervus nippon)

Deer have been a major resource for human populations for thousands of years. Anthropogenic activities, such as hunting, have influenced the genetic structure and distribution of deer populations. In Japan, wild Japanese sika deer (Cervus nippon) have been hunted since ancient times but have also been historically protected as sacred animals in several sanctuaries. Sika deer have been protected for over a thousand years in the religious sanctuary around the Kasuga Taisha Shrine on the Kii Peninsula, located in the center of Japan. Here, we used short sequence repeats (SSR) and mitochondrial DNA (mtDNA) to investigate the genetic diversity, population structure, and demography of Japanese sika deer inhabiting the Kii Peninsula, Japan, and discuss possible anthropogenic influences. Using SSR, three distinct genetic groups were distinguished on the Kii Peninsula: an Eastern genetic group, a Western genetic group, and an isolated genetic group with individuals in the religious sanctuary of Kasuga Taisha Shrine in Nara city. The isolated genetic sanctuary group had only the mtDNA haplotype S4. The SSR genotype data suggested a newer divergence time of the genetic groups of the religious sanctuary than would have occurred as a result of Late Quaternary climate change. This time scale coincided with the establishment of the sanctuary with Kasuga Taisha Shrine. Thus, the religious protection conserved genetic variation over a thousand years.

Does harvest affect genetic diversity in grey wolves?

Harvest can affect vital rates such as reproduction and survival, but also genetic measures of individual and population health. Grey wolves (Canis lupus) live and breed in groups, and effective population size is a small fraction of total abundance. As a result, genetic diversity of wolves may be particularly sensitive to harvest. We evaluated how harvest affected genetic diversity and relatedness in wolves. We hypothesized that harvest would (a) reduce relatedness of individuals within groups in a subpopulation but increase relatedness of individuals between groups due to increased local immigration, (b) increase individual heterozygosity and average allelic richness across groups in subpopulations and (c) add new alleles to a subpopulation and decrease the number of private alleles in subpopulations due to an increase in breeding opportunities for unrelated individuals. We found harvest had no effect on observed heterozygosity of individuals or allelic richness at loci within subpopulations but was associated with a small, biologically insignificant effect on within-group relatedness values in grey wolves. Harvest was, however, positively associated with increased relatedness of individuals between groups and a net gain (+16) of alleles into groups in subpopulations monitored since harvest began, although the number of private alleles in subpopulations overall declined. Harvest likely created opportunities for wolves to immigrate into nearby groups and breed, thereby making groups in subpopulations more related over time. Harvest appears to affect genetic diversity in wolves at the group and population levels, but its effects are less apparent at the individual level given the population sizes we studied.

Genetic population structure of invasive raccoons (Procyon lotor) in Hokkaido, Japan: Unique phenomenon caused by pet escape or abandonment

Phylogeographic studies can resolve relationships between genetic population structure of organisms and geographical distributions. Raccoons have become feral in Japan, and in Hokkaido island, they have been rapidly increasing in number and spreading since the 1970s. We analyzed mitochondrial (mtDNA) and microsatellite DNA to understand the current phylogenetic distribution and invasive founder events. Overall, Hokkaido raccoons maintained high genetic diversity (i.e., the level of heterozygosity was comparable to the original habitat, North America). Based on mtDNA distribution and microsatellite diversity, Hokkaido raccoons were divided into six management units. However, mtDNA haplotype distributions and genetic structures based on microsatellites did not always correspond to each other (e.g., two geographically and genetically separated populations showed similar mtDNA distributions). In addition, a high degree of genetic admixture was observed in every unit, and the degree of genetic differentiation was low even between regions separated by long distances. Compared with other countries in Europe where genetic distribution of introduced raccoons is more clearly structured, the current results represent a unique and complex phenomenon of pet escape/abandonment in Hokkaido: i.e., genetically related colonies were introduced into multiple regions as founder events, resulting in the current state in which raccoons are not clearly genetically differentiated even 40 years after introduction.

Gongylonema infection of wild mammals in Japan and Sardinia (Italy)

The gullet worms, classical Gongylonema pulchrum and newly differentiated Gongylonema nepalensis, are prevalent in various mammals in Japan and Sardinia, Italy, respectively. The former species is cosmopolitan in distribution, dwelling in the mucosa of the upper digestive tract of a variety of domestic and wild mammals, and also humans. At present, the geographical distribution of G. nepalensis is known in Nepal and Sardinia, with the nematode having been recorded from the oesophagus of water buffaloes (Nepal), cattle, sheep, goats and wild mouflon (Sardinia). To clarify their natural transmission cycles among domestic and wild mammals, the present study analysed the ribosomal RNA gene (rDNA) and mitochondrial cytochrome c oxidase subunit 1 gene (cox1) of worms of various origins: G. pulchrum worms from sika deer, wild boars, Japanese macaques, and feral alien Reeves's muntjacs in Japan, and G. nepalensis worms from a red fox and a wild boar in Sardinia. Although the internal transcribed spacer (ITS) regions of rDNA and partial cox1 nucleotide sequences of G. pulchrum from native wild mammals in Japan were distinct from those of the worms in cattle, the worms from feral alien Reeves's muntjacs showed the cattle-type ITS genotype and cox1 cattle-I and II haplotypes. The rDNA and cox1 nucleotide sequences of G. nepalensis from a red fox in Sardinia were almost identical to those of the worms from domestic and wild ruminants on the island. The ecological interaction between domestic and wild mammals and their susceptibility to different Gongylonema spp. must be considered when trying to elucidate this spirurid's transmission dynamics in nature.

Metapopulation inbreeding dynamics, effective size and subpopulation differentiation--A general analytical approach for diploid organisms

Motivated by problems in conservation biology we study genetic dynamics in structured populations of diploid organisms (monoecious or dioecious). Our analysis provides an analytical framework that unifies substantial parts of previous work in terms of exact identity by descent (IBD) and identity by state (IBS) recursions. We provide exact conditions under which two structured haploid and diploid populations are equivalent, and some sufficient conditions under which a dioecious diploid population can be treated as a monoecious diploid one. The IBD recursions are used for computing local and metapopulation inbreeding and coancestry effective population sizes and for predictions of several types of fixation indices over different time horizons.

Major decline in marine and terrestrial animal consumption by brown bears (Ursus arctos)

Human activities have had the strongest impacts on natural ecosystems since the last glacial period, including the alteration of interspecific relationships such as food webs. In this paper, we present a historical record of major alterations of trophic structure by revealing millennium-scale dietary shifts of brown bears (Ursus arctos) on the Hokkaido islands, Japan, using carbon, nitrogen, and sulfur stable isotope analysis. Dietary analysis of brown bears revealed that salmon consumption by bears in the eastern region of Hokkaido significantly decreased from 19% to 8%. In addition, consumption of terrestrial animals decreased from 56% to 5% in western region, and 64% to 8% in eastern region. These dietary shifts are likely to have occurred in the last approximately 100–200 years, which coincides with the beginning of modernisation in this region. Our results suggest that human activities have caused an alteration in the trophic structure of brown bears in the Hokkaido islands. This alteration includes a major decline in the marine-terrestrial linkage in eastern region, and a loss of indirect-interactions between bears and wolves, because the interactions potentially enhanced deer predation by brown bears.

Assessing prehistoric genetic structure and diversity of North American elk (Cervus elaphus) populations in Alberta, Canada

North American elk (Cervus elaphus L., 1758) are an important component of Canada’s natural ecosystems. Overhunting and habitat decline in the 19th century led to the near eradication of Rocky Mountain elk (Cervus elaphus nelsoni Bailey, 1935) and Manitoban elk (Cervus elaphus manitobensis Millais, 1915) within Alberta. Though elk populations have been restored within provincial and national parks, it is unknown to what degree historic population declines affected overall genetic diversity and population structuring of the two subspecies. This study targeted 551 bp of mitochondrial D-loop DNA from 50 elk remains recovered from 16 archaeological sites (2260 BCE (before common era) to 1920 CE (common era)) to examine the former genetic diversity and population structure of Alberta’s historic elk populations. Comparisons of ancient and modern haplotype and nucleotide diversity suggest that historic population declines reduced the mitochondrial diversity of Manitoban elk, while translocation of animals from Yellowstone National Park in the early 20th century served to maintain the diversity of Rocky Mountain populations. Gene flow between the two subspecies was significantly higher in the past than today, suggesting that the two subspecies previously formed a continuous population. These data on precontact genetic diversity and gene flow in Alberta elk provide essential baseline data integral for elk management and conservation in the province.

Genetic variation of Gongylonema pulchrum from wild animals and cattle in Japan based on ribosomal RNA and mitochondrial cytochrome c oxidase subunit I genes

The gullet worm (Gongylonema pulchrum) has been recorded from a variety of mammals worldwide, including monkeys and humans. Due to its wide host range, it has been suggested that the worm may be transmitted locally to any mammalian host by chance. To investigate this notion, the ribosomal RNA gene (rDNA), mainly regions of the internal transcribed spacers (ITS) 1 and 2, and a cytochrome c oxidase subunit I (COI) region of mitochondrial DNA of G. pulchrum were characterized using parasites from the following hosts located in Japan: cattle, sika deer, wild boars, Japanese macaques, a feral Reeves's muntjac and captive squirrel monkeys. The rDNA nucleotide sequences of G. pulchrum were generally well conserved regardless of their host origin. However, a few insertions/deletions of nucleotides along with a few base substitutions in the ITS1 and ITS2 regions were observed in G. pulchrum from sika deer, wild boars and Japanese macaques, and those differed from G. pulchrum in cattle, the feral Reeves's muntjac and captive squirrel monkeys. The COI sequences of G. pulchrum were further divided into multiple haplotypes and two groups of haplotypes, i.e. those from a majority of sika deer, wild boars and Japanese macaques and those from cattle and zoo animals, were clearly differentiated. Our findings indicate that domestic and sylvatic transmission cycles of the gullet worm are currently present, at least in Japan.

Whole genome SNP discovery and analysis of genetic diversity in Turkey (Meleagris gallopavo)

Background

The turkey (Meleagris gallopavo) is an important agricultural species and the second largest contributor to the world’s poultry meat production. Genetic improvement is attributed largely to selective breeding programs that rely on highly heritable phenotypic traits, such as body size and breast muscle development. Commercial breeding with small effective population sizes and epistasis can result in loss of genetic diversity, which in turn can lead to reduced individual fitness and reduced response to selection. The presence of genomic diversity in domestic livestock species therefore, is of great importance and a prerequisite for rapid and accurate genetic improvement of selected breeds in various environments, as well as to facilitate rapid adaptation to potential changes in breeding goals. Genomic selection requires a large number of genetic markers such as e.g. single nucleotide polymorphisms (SNPs) the most abundant source of genetic variation within the genome.

Results

Alignment of next generation sequencing data of 32 individual turkeys from different populations was used for the discovery of 5.49 million SNPs, which subsequently were used for the analysis of genetic diversity among the different populations. All of the commercial lines branched from a single node relative to the heritage varieties and the South Mexican turkey population. Heterozygosity of all individuals from the different turkey populations ranged from 0.17-2.73 SNPs/Kb, while heterozygosity of populations ranged from 0.73-1.64 SNPs/Kb. The average frequency of heterozygous SNPs in individual turkeys was 1.07 SNPs/Kb. Five genomic regions with very low nucleotide variation were identified in domestic turkeys that showed state of fixation towards alleles different than wild alleles.

Conclusion

The turkey genome is much less diverse with a relatively low frequency of heterozygous SNPs as compared to other livestock species like chicken and pig. The whole genome SNP discovery study in turkey resulted in the detection of 5.49 million putative SNPs compared to the reference genome. All commercial lines appear to share a common origin. Presence of different alleles/haplotypes in the SM population highlights that specific haplotypes have been selected in the modern domesticated turkey.

Genetic variation and population structure of the Japanese sika deer (Cervus nippon) in the Tohoku District based on mitochondrial D-loop sequences

The sika deer (Cervus nippon) once inhabited the entire Tohoku District, the northeastern part of the main island of Japan. Currently, they are isolated as three discontinuous populations on Mt. Goyo, the Oshika Peninsula, and Kinkazan Island. To assess the genetic diversity and relationships among the sika deer populations in the Tohoku District, we analyzed the mitochondrial DNA D-loop sequences from 177 individuals. We detected a total of five haplotypes. Three haplotypes were present in the population from Mt. Goyo at a haplotype diversity of 0.235 ± 0.061, two haplotypes in the population from the Oshika Peninsula at 0.171 ± 0.064, and only one haplotype was detected in the population from the Kinkazan Island. A significant genetic differentiation was observed among all population pairs. Collectively, our data supports the observed population bottlenecks in the past. Four of the five haplotypes were specific to one of the three populations, whereas only one haplotype was shared between the Mt. Goyo and the Oshika Peninsula populations. This common haplotype may indicate a common ancestral population in the Tohoku District. Conversely, the D-loop haplotypes were completely different among the Kinkazan Island and Oshika Peninsula populations. The lack of a shared haplotype indicates that female gene flow between the two populations is very limited and that the 0.6 km strait acts as a strong barrier.

Sex determination of sika deer (Cervus nippon yesoensis) using nested PCR from feces collected in the field

We describe a method for determining the sex of sika deer (Cervus nippon yesoensis) from feces collected in the field. Using a nested polymerase chain reaction (nested PCR), partial sequences of the sex determination region of the Y chromosome (SRY) gene and X zinc finger protein (ZFX) gene were amplified. In 19 individuals with sex information, the correct sex was successfully detected and sequences of target amplicons were completely matched between muscle and feces from the rectum. Among 75 fecal samples collected noninvasively in the field, 68-71 samples (90.7-94.7%) yielded successful sex determinations. Using this technique, feces collected in the field would enhance the utility of genetic analysis. As a result, instead of biomaterials, these samples can serve as new or alternative materials. Finally, it can be expected that this technique will contribute to reveal in advanced detail of the population dynamics and genetic diversity that needed to carry out effective population control and to reduce the extinction risk of sika deer.

Genetic structure of the asiatic black bear in Japan using mitochondrial DNA analysis

The genetic structure of the Asiatic black bear (Ursus thibetanus) in Japan was studied to understand the events that occurred during its evolution. The left domain of the mitochondrial control region (about 240 bp) was sequenced, defining 27 haplotypes that consisted of 23 haplotypes from 333 bears in Japan and 22 bears in the Asian continent. The network tree of the control region indicated that the Japanese population formed a distinct clade from the continental population. The phylogeographic analysis of the haplotypes indicated that the Shikoku and Kii Hanto populations had diverged during the initial phase from the ancestral population. After the 3 dominant haplotypes were rapidly distributed throughout Japan in the early stage of the population dispersal, the Japanese population diverged into eastern and western populations. Using the entire mitochondrial cytochrome b sequence, divergence time between the Japanese and the Continental populations suggested that the Japanese population might have colonized into Japan through the land bridge from the Korean Peninsula around 500 ka, which is consistent with paleontological evidence. Our finding that bears in western Japan exhibit lower genetic diversity and higher levels of genetic differentiation than bears in eastern Japan provides a vital contribution to conservation policy for these isolated populations.

Mitochondrial phylogeography of the Eurasian beaver Castor fiber L

Nucleotide variation in an approximately 490 bp fragment of the mitochondrial DNA control region (mtDNA CR) was used to describe the genetic variation and phylogeographical pattern in the Eurasian beaver (Castor fiber) over its entire range. The sampling effort was focused on the relict populations that survived a drastic population bottleneck, caused by overhunting, at the end of the 19th century. A total of 152 individuals grouped into eight populations representing all currently recognized subspecies were studied. Sixteen haplotypes were detected, none of them shared among populations. Intrapopulation sequence variation was very low, most likely a result of the severe bottleneck. Extreme genetic structure could result from human-mediated extinction of intermediate populations, but it could also be an effect of prior substantial structuring of the beaver populations with watersheds of major Eurasian rivers acting as barriers to gene flow. Phylogenetic analysis revealed the presence of two mtDNA lineages: eastern (Poland, Lithuania, Russia and Mongolia) and western (Germany, Norway and France), the former comprising more divergent haplotypes. The low level of sequence divergence of the entire cytochrome b gene among six individuals representing six subspecies suggests differentiation during the last glacial period and existence of multiple glacial refugia. At least two evolutionary significant units (ESU) can be identified, the western and the eastern haplogroup. The individual relict populations should be regarded as management units, the eastern subspecies possibly also as ESUs. Guidelines for future translocations and reintroductions are proposed.

Surveillance of Chronic Wasting Disease in Sika Deer, Cervus nippon, from Tokachi District in Hokkaido

Surveillance of chronic wasting disease (CWD) was conducted by performing Western blot analysis of tissue samples from 136 sika deer (Cervus nippon) killed by hunters in the Tokachi district of Hokkaido Island. No prion protein (PrPSc) associated with CWD was detected in any of the samples. To assess amino acid polymorphisms of the sika deer PrP gene, nucleotide sequencing of the PrP gene was performed. The only amino acid polymorphisms detected were 3 silent mutations at nucleotide positions 63, 225 and 408. These results suggest that sika deer in the Tokachi district are genetically homogeneous, and are not infected with CWD.