Evolution of MHC class I genes in Eurasian badgers, genus Meles (Carnivora, Mustelidae).

Because of their role in immune defense against pathogens, major histocompatibility complex (MHC) genes are useful in evolutionary studies on how wild vertebrates adapt to their environments. We investigated the molecular evolution of MHC class I (MHCI) genes in four closely related species of Eurasian badgers, genus Meles. All four species of badgers showed similarly high variation in MHCI sequences compared to other Carnivora. We identified 7-21 putatively functional MHCI sequences in each of the badger species, and 2-7 sequences per individual, indicating the existence of 1-4 loci. MHCI exon 2 and 3 sequences encoding domains α1 and α2 exhibited different clade topologies in phylogenetic networks. Non-synonymous nucleotide substitutions at codons for antigen-binding sites exceeded synonymous substitutions for domain α1 but not for domain α2, suggesting that the domains α1 and α2 likely had different evolutionary histories in these species. Positive selection and recombination seem to have shaped the variation in domain α2, whereas positive selection was dominant in shaping the variation in domain α1. In the separate phylogenetic analyses for exon 2, exon 3, and intron 2, each showed three clades of Meles alleles, with rampant trans-species polymorphism, indicative of the long-term maintenance of ancestral MHCI polymorphism by balancing selection.

Low genetic variation in the MHC class II DRB gene and MHC-linked microsatellites in endangered island populations of the leopard cat (Prionailurus bengalensis) in Japan.

Isolated populations of the leopard cat (Prionailurus bengalensis) on Tsushima and Iriomote islands in Japan are classified as subspecies P. b. euptilurus and P. b. iriomotensis, respectively. Because both populations have decreased to roughly 100, an understanding of their genetic diversity is essential for conservation. We genotyped MHC class II DRB exon 2 and MHC-linked microsatellite loci to evaluate the diversity of MHC genes in the Tsushima and Iriomote cat populations. We detected ten and four DRB alleles in these populations, respectively. A phylogenetic analysis showed DRB alleles from both populations to be closely related to those in other felid DRB lineages, indicating trans-species polymorphism. The MHC-linked microsatellites were more polymorphic in the Tsushima than in the Iriomote population. The MHC diversity of both leopard cat populations is much lower than in the domestic cat populations on these islands, probably due to inbreeding associated with founder effects, geographical isolation, or genetic drift. Our results predict low resistance of the two endangered populations to new pathogens introduced to the islands.

Genetic Diversity of MHC Class II DRB1 Exon 2 in the Red Fox (Vulpes vulpes) on Hokkaido, Japan.

To assess the genetic diversity of the red fox (Vulpes vulpes) population on Hokkaido Island, northern Japan, we examined genotypes of exon 2 of the major histocompatibility complex (MHC) class II DRB1 gene for 232 individuals and identified 17 novel alleles. The subpopulation in the Southern area was genetically differentiated from those in all other areas on Hokkaido, suggesting isolation in southern Hokkaido. In addition, the Southern subpopulation is lower in genetic diversity than the other subpopulations, possibly resulting from purifying selection and/or a recent bottleneck. The non-synonymous substitutions exceeded the synonymous substitutions for codons encoding antigen-binding sites (ABSs) in exon 2, indicating that the red fox DRB1 alleles have evolved under positive selection. In a Bayesian phylogenetic tree, although most of the DRB1 alleles from the Hokkaido red fox were contained within a red fox-like canid (Vulpes) clade, some belonged to another canid clade. This means trans-species polymorphism maintained by balancing selection. Our results showed the local variability and the presence of selection on the MHC gene in this population, which contributes to the understanding of the historical background and the molecular evolution.

Genetic variation of major histocompatibility complex genes in the endangered red-crowned crane.

Populations that have drastically decreased in the past often have low genetic variation, which may increase the risk of extinction. The genes of major histocompatibility complex (MHC) play an important role in the adaptive immune response of jawed vertebrates. Maintenance of adaptive genetic diversity such as that of MHC genes is important for wildlife conservation. Here, we determined genotypes of exon 3 of MHC class IA genes (MHCIA) and exon 2 of MHC class IIB genes (MHCIIB) to evaluate genetic variation of the endangered red-crowned crane population on Hokkaido Island, Japan, which experienced severe population decline in the past. We identified 16 and 6 alleles of MHCIA and MHCIIB, respectively, from 152 individuals. We found evidence of a positive selection at the antigen-binding sites in MHCIA exon 3 and MHCIIB exon 2. The phylogenetic analyses indicated evidence of trans-species polymorphism among the crane MHC genes. The genetic variability in both classes of MHC genes at the population level was low. No geographic structure was found based on the genetic diversity of microsatellite and MHC genes. Our study provides useful data for the optimal management of the red-crowned crane population in Hokkaido and can contribute to future studies on MHC genes of the continental populations of the red-crowned crane and other crane species.

Gene duplication and concerted evolution of mitochondrial DNA in crane species.

The gene duplication in mitochondrial DNA (mtDNA) has been reported in diverse bird taxa so far. Although many phylogenetic and population genetic analyses of cranes were carried out based on mtDNA diversity, whether mtDNA contains duplicated regions is unknown. To address the presence or absence of gene duplication in cranes and investigate the molecular evolutionary features of crane mtDNA, we analyzed the gene organization and the molecular phylogeny of mtDNA from 13 crane species. We found that the mtDNA in 13 crane species shared a tandem duplicated region, which consists of duplicated sequence sets including cytochrome b (Cytb), NADH6, control region (CR) and three genes of tRNA. The gene order in the duplicated region was identical among all the 13 crane species, and the nucleotide sequences found within each individual showed high similarities. In addition, phylogenetic trees based on homologous sequences of CR and Cytb indicated the possibility of concerted evolution among the duplicated genes. The results suggested that the duplication event occurred in the common ancestor of crane species or some older ancestors.

Low Genetic Variation of Red-Crowned Cranes on Hokkaido Island, Japan, Over the Hundred Years.

The red-crowned crane (Grus japonensis) is recognized internationally as an endangered species. Migratory populations breed in eastern Russia and northeastern China, whereas the resident population inhabits the island of Hokkaido, Japan. Although the population inhabiting Hokkaido had experienced a severe bottleneck by the end of the 19th century, the population size has recovered to about 1500 and continues to increase now thanks to conservation efforts. A previous study reported that no marked genetic differences were seen in the island population, and that the genetic variation of the whole population on Hokkaido was lower than that of the continental population. However, the precise genetic structure of the island population in the past or near present remains unclear. To better understand the spatiotemporal changes in the genetic structure of the island population, we performed mitochondrial DNA (mtDNA) analyses using stuffed specimens (years 1878-2001) and tissue or blood samples (years 1970-2014). We found three haplotypes in the island population, one of which was a novel mtDNA haplotype in 1997 and 2007 samples. In addition, there was no clear difference in the haplotype frequency through the time span. These results suggest that the low genetic variation of the island population persisted for the last hundred years. It is thus nearly impossible for the island population to recover its genetic variation in isolation. Conservation plans for this species should therefore include the promotion of genetic exchanges between the continental and island populations, such as through artificial introduction to Hokkaido.

Duplication and Variation in the Major Histocompatibility Complex Genes in Blakiston's Fish Owl, Bubo blakistoni.

The major histocompatibility complex (MHC) includes many genes that are essential for the adaptive immune system, and variation in the antigen binding site (ABS) is related to resistance against pathogens. In the present study, quantitative real-time PCR indicated a larger number of MHC gene copies in the endangered population of Blakiston's fish owl (Bubo blakistoni) than in five other owl species, and massively parallel pyrosequencing detected more MHC class IIß per individual alleles in B. blakistoni than in the other species. A chromosomal fluorescence in situ hybridization (FISH) analysis showed that the MHC class I and class IIß loci are closely linked on a single pair of microchromosomes, indicating that the MHC genes were tandemly duplicated in a limited chromosomal region. Because B. blakistoni has twice as many MHC genes as its sister species, the tawny fish owl (Bubo flavipes), the duplication of MHC genes occurred after these species diverged by speciation. A Bayesian molecular phylogenetic analysis showed that the DAB1 and DAB2 lineages of MHC class IIß alleles from various strigid species each formed a separate clade, indicating that the two allelic lineages preceded the radiation of Strigidae and evolved as paralogs. By contrast, the ABS sequences did not form distinct clades between DAB1 and DAB2 alleles but were intermixed, presumably due to gene conversion. Despite the low diversity of alleles per locus, B. blakistoni had many lineages of MHC class IIß alleles. Gene duplication increases variation in the MHC genes in this species, and could have facilitated adaptation in small populations.

A Comparison of Visual and Genetic Techniques for Identifying Japanese Marten Scats - Enabling Diet Examination in Relation to Seasonal Food Availability in a Sub-Alpine Area of Japan.

We compared the reliability of visual diagnostic criteria to DNA diagnostic techniques, including newly designed primers, to discriminate Japanese marten (Martes melampus) feces from those of other sympatric carnivore species. Visual criteria proved > 95% reliable for fresh, odoriferous scats in good condition. Based upon this verification, we then examined if and how Japanese marten diet differs among seasons at high elevation study site (1500-2026 m). We also considered how intra-specific competition with the Japanese red fox (Vulpes vulpes japonica) may shape marten feeding ecology. From 120 Japanese marten fecal samples, high elevation diet comprised (frequency of occurrence) 30.6-66.0% mammals, 41.0-72.2% insects and 10.6-46.2% fruits, subject to seasonal variation, with a Shannon-Weaver index value of 2.77. These findings contrast substantially to seasonal marten diet reported in adjacent lowland regions (700-900 m), particularly in terms of fruit consumption, showing the trophic adaptability of the Japanese marten. We also noted a substantial dietary overlap with the red fox (n = 26 scats) with a Shannon-Weaver index of 2.61, inferring little trophic niche mutual exclusion (trophic niche overlap: 0.95), although some specific seasonal prey selection differences were likely related to relative differences in body size between foxes and martens. This additional information on the feeding ecology of the Japanese marten enables a better assessment of the specific risks populations face in mountainous regions.

Population Structure of the Raccoon Dog on the Grounds of the Imperial Palace, Tokyo, Revealed by Microsatellite Analysis of Fecal DNA.

The raccoon dog (Nyctereutes procyonoides, Canidae, Carnivora) is highly adaptable to urban environments. Populations of carnivorans inhabiting urban areas sometimes differ ecologically and genetically from those in rural areas. However, there is little information on urban raccoon dogs. This study focused on raccoon dog populations in Tokyo, one of the most highly urbanized cities in the world. We examined the genotypes of 10 microsatellites for 101 fecal samples from raccoon dogs inhabiting the grounds of the Imperial Palace, a green space in central Tokyo. We successfully genotyped 58 samples originating from 31 individuals. We also analyzed muscle tissue samples from raccoon dogs from the grounds of the Imperial Palace, the Akasaka Imperial Grounds (a green space close to the Imperial Palace), and the surrounding urban area, and then investigated the genetic structure and diversity of these populations, and the genetic differentiation among them. The population on the grounds of the Imperial Palace was genetically differentiated from that in the Akasaka Imperial Grounds, suggesting that the roads and buildings act as barriers to gene flow. In addition, the population on the grounds of the Imperial Palace showed greater genetic difference from that in the surrounding area than that in the Akasaka Imperial Grounds. We speculate that the moats around the Imperial Palace restrict individual ranges within the palace grounds and limit migration and gene flow to other areas.

Isolation and characterization of major histocompatibility complex class II B genes in cranes.

In this study, we isolated and characterized the major histocompatibility complex (MHC) class II B genes in cranes. Genomic sequences spanning exons 1 to 4 were amplified and determined in 13 crane species and three other species closely related to cranes. In all, 55 unique sequences were identified, and at least two polymorphic MHC class II B loci were found in most species. An analysis of sequence polymorphisms showed the signature of positive selection and recombination. A phylogenetic reconstruction based on exon 2 sequences indicated that trans-species polymorphism has persisted for at least 10 million years, whereas phylogenetic analyses of the sequences flanking exon 2 revealed a pattern of concerted evolution. These results suggest that both balancing selection and recombination play important roles in the crane MHC evolution.

Molecular phylogeography of the brown bear (Ursus arctos) in Northeastern Asia based on analyses of complete mitochondrial DNA sequences.

To further elucidate the migration history of the brown bears (Ursus arctos) on Hokkaido Island, Japan, we analyzed the complete mitochondrial DNA (mtDNA) sequences of 35 brown bears from Hokkaido, the southern Kuril Islands (Etorofu and Kunashiri), Sakhalin Island, and the Eurasian Continent (continental Russia, Bulgaria, and Tibet), and those of four polar bears. Based on these sequences, we reconstructed the maternal phylogeny of the brown bear and estimated divergence times to investigate the timing of brown bear migrations, especially in northeastern Eurasia. Our gene tree showed the mtDNA haplotypes of all 73 brown and polar bears to be divided into eight divergent lineages. The brown bear on Hokkaido was divided into three lineages (central, eastern, and southern). The Sakhalin brown bear grouped with eastern European and western Alaskan brown bears. Etorofu and Kunashiri brown bears were closely related to eastern Hokkaido brown bears and could have diverged from the eastern Hokkaido lineage after formation of the channel between Hokkaido and the southern Kuril Islands. Tibetan brown bears diverged early in the eastern lineage. Southern Hokkaido brown bears were closely related to North American brown bears.

Limited phylogenetic distribution of a long tandem-repeat cluster in the mitochondrial control region in Bubo (Aves, Strigidae) and cluster variation in Blakiston's fish owl (Bubo blakistoni).

To investigate the phylogenetic position of Blakiston's fish owl (Bubo blakistoni), we sequenced the mitochondrial (mt) DNA control region and cytochrome b (cyt b) for nine Bubo species. Maximum-likelihood analyses of combined control region and cyt b sequences, and cyt b sequences alone, showed that species formerly placed in genus Ketupa comprise a monophyletic group. Unexpectedly, we discovered a long cluster of 20-25 tandem repeat units 77 or 78bp long in the third control region domain in four of the nine Bubo species for which the control region was sequenced (B. blakistoni, B. flavipes, and B. ketupu in the Ketupa clade; B. lacteus), leading to overall control region lengths of 3.0-3.8kpb estimated from agarose gel electrophoresis. The control region in B. lacteus is the longest (3.8kbp) reported to date in vertebrates. Sequencing of eight repeat units at each end of the cluster in 20 B. blakistoni individuals detected several types of repeat units 77 or 78bp long, and six patterns in the order of unit types. The occurrence of a repeat cluster in all three species examined in the Ketupa clade suggests their common ancestor also had a cluster, whereas a maximum parsimony tree showed repeat-unit types grouping by species, rather than by paralog groups, suggesting independent origins of the clusters. We reconcile these results with a turnover model, in which the range in cluster-length variation and unit types at the 5' end are hypothetically functionally constrained by the protein-binding function of the control region, but otherwise there is a continual turnover of units in evolutionary time, with new unit types arising through mutations, proliferating by duplication of single and double repeat blocks, and being lost through deletion. Estimated free energies for reconstructed secondary structures of single and especially pairs of repeat units were higher than for homologous single-unit blocks in species lacking a repeat cluster, supporting slipped-strand mispairing as the mechanism of cluster turnover.

Diversity of mitochondrial D-loop haplotypes from ancient Thracian horses in Bulgaria.

The domestication of the horse began possibly more than 5000 years ago in the western part of the Eurasian steppe, and according to the leading hypothesis, horses first spread from the Steppe toward the region of the Thracian culture, starting in the second half of the 2nd millennium BCE and flourished from the fifth to first centuries BCE, mainly located in present-day Bulgaria. We analyzed 17 horse bone remains excavated from Thracian archaeological sites (fourth to first centuries BCE) in Bulgaria and successfully identified 17 sequences representing 14 different haplotypes of the mitochondrial D-loop. Compared with the mtDNA haplotypes of modern horses around the world, ancient Thracian horses in Bulgaria are thought to be more closely related to modern horses of Southern Europe and less related to those of Central Asia. In addition, the haplotypes we obtained represented 11 previously reported modern horse mtDNA haplogroups: A, B, D, E, G, H, I, L, N, P, and Q. All the haplogroups contain modern and regionally predominant haplotypes occurring in Europe, the Middle East, and Central Asia. Our results indicate that Thracian horses in Bulgaria have had relatively high genetic diversity and are closely related to modern horse breeds.

Temporal changes of genetic population structure and diversity in the endangered Blakiston's fish owl (Bubo blakistoni) on Hokkaido Island, Japan, revealed by microsatellite analysis.

The Blakiston's fish owl (Bubo blakistoni) population on Hokkaido Island, Japan, decreased to less than one hundred individuals over the last century due to habitat disruption by human activity. Although the ongoing conservation management has slightly restored the population, it remains endangered. In order to assess the genetic variation and population structure of the Blakiston's fish owl in Hokkaido, we genotyped eight microsatellite loci on 120 individuals sampled over the past three decades. The genotype data set showed low levels of genetic variation and gene flow among the geographically isolated five subpopulations. Comparative analysis of past and current populations indicated that some alleles shared by past individuals had been lost, and that genetic variation had declined over the last three decades. The result suggests that the genetic decline may have resulted from inbreeding and/or genetic drift due to bottlenecks in the Hokkaido population. The present study provides invaluable genetic information for the conservation and management of the endangered Blakiston's fish owl in Hokkaido.

Genetic distinctness and variation in the Tsushima Islands population of the Japanese marten, Martes melampus (Carnivora: Mustelidae), revealed by microsatellite analysis.

A carnivoran mammal endemic to Japan, the Japanese marten (Martes melampus) is native in forested regions on Honshu, Shikoku, Kyushu (main islands of Japan), and the Tsushima Islands. The Tsushima population is classified as a different subspecies (M. m. tsuensis) from populations on the main islands (M. m. melampus). To elucidate the genetic structure of the Tsushima population, we genotyped 101 individuals from the Tsushima Islands and 43 individuals from Honshu and Kyushu using 10 microsatellite loci, and performed population genetic analyses on the genotype data. Genetic diversity was lower in the Tsushima population than in three geographic populations on the main islands: heterozygosity was 0.189-0.364 in the former, compared to 0.457-0.747 in the latter. In addition, high pairwise Fst values (0.485-0.682) and Nei's standard distance (0.550-1.183) between the Tsushima and main-island populations indicated a high degree of genetic differentiation. Finally, a Bayesian clustering analysis showed that the Tsushima population is apparently differentiated from the main-island populations and comprises two genetic clusters. A factorial correspondence analysis corroborated these results. Our results suggest that restricted gene flow or inbreeding may have reduced genetic diversity in the Tsushima population, which has been geographically isolated from the main-island populations since the formation of Tsushima Strait.

Comparative molecular phylogeny and evolution of sex chromosome DNA sequences in the family Canidae (Mammalia: Carnivora).

To investigate the molecular phylogeny and evolution of the family Canidae, nucleotide sequences of the zinc-finger-protein gene on the Y chromosome (ZFY, 924-1146 bp) and its homologous gene on the X chromosome (ZFX, 834-839 bp) for twelve canid species were determined. The phylogenetic relationships among species reconstructed by the paternal ZFY sequences closely agreed with those by mtDNA and autosomal DNA trees in previous reports, and strongly supported the phylogenetic affinity between the wolf-like canids clade and the South American canids clade. However, the branching order of some species differed between phylogenies of ZFY and ZFX genes: Cuon alpinus and Canis mesomelas were included in the wolf-like canid clades in the ZFY tree, whereas both species were clustered in a group of Chrysocyon brachyurus and Speothos venaticus in the ZFX tree. The topology difference between ZFY and ZFX trees may have resulted from the two-times higher substitution rate of the former than the latter, which was clarified in the present study. In addition, two types of transposable element sequence (SINE-I and SINE-II) were found to occur in the ZFY final intron of the twelve canid species examined. Because the SINE-I sequences were shared by all the species, they may have been inserted into the ZFY of the common ancestor before species radiation in Canidae. By contract, SINE-II found in only Canis aureus could have been inserted into ZFY independently after the speciation. The molecular diversity of SINE sequences of Canidae reflects evolutionary history of the species radiation.

Diversity of the MHC class II DRB gene in the wolverine (Carnivora: Mustelidae: Gulo gulo) in Finland.

The wolverine (Gulo gulo) in Finland has undergone significant population declines in the past. Since major histocompatibility complex (MHC) genes encode proteins involved in pathogen recognition, the diversity of these genes provides insights into the immunological fitness of regional populations. We sequenced 862 amplicons (242 bp) of MHC class II DRB exon 2 from 32 Finnish wolverines and identified 11 functional alleles and three pseudogenes. A molecular phylogenetic analysis indicated trans-species polymorphism, and PAML and MEME analyses indicated positive selection, suggesting that the Finnish wolverine DRB genes have evolved under balancing and positive selection. In contrast to DRB gene analyses in other species, allele frequencies in the Finnish wolverines clearly indicated the existence of two regional subpopulations, congruent with previous studies based on neutral genetic markers. In the Finnish wolverine, rapid population declines in the past have promoted genetic drift, resulting in a lower genetic diversity of DRB loci, including fewer alleles and positively selected sites, than other mustelid species analyzed previously. Our data suggest that the MHC region in the Finnish wolverine population was likely affected by a recent bottleneck.

Genetic features of ancient West Siberian people of the Middle Ages, revealed by mitochondrial DNA haplogroup analysis.

In order to investigate the genetic features of ancient West Siberian people of the Middle Ages, we studied ancient DNA from bone remains excavated from two archeological sites in West Siberia: Saigatinsky 6 (eighth to eleventh centuries) and Zeleny Yar (thirteenth century). Polymerase chain reaction amplification and nucleotide sequencing of mitochondrial DNA (mtDNA) succeeded for 9 of 67 specimens examined, and the sequences were assigned to mtDNA haplogroups B4, C4, G2, H and U. This distribution pattern of mtDNA haplogroups in medieval West Siberian people was similar to those previously reported in modern populations living in West Siberia, such as the Mansi, Ket and Nganasan. Exact tests of population differentiation showed no significant differences between the medieval people and modern populations in West Siberia. The findings suggest that some medieval West Siberian people analyzed in the present study are included in direct ancestral lineages of modern populations native to West Siberia.

Phylogeographic sympatry and isolation of the Eurasian badgers (Meles, Mustelidae, Carnivora): Implications for an alternative analysis using maternally as well as paternally inherited genes.

In the present study, to further understand the phylogenetic relationships among the Eurasian badgers (Meles, Mustelidae, Carnivora), which are distributed widely in the Palearctic, partial sequences of the mitochondrial DNA (mtDNA) control region (539-545 base-pairs) as a maternal genetic marker, and the sex-determining region on the Y-chromosome gene (SRY: 1052-1058 base-pairs), as a paternal genetic marker, were examined. The present study revealed ten SRY haplotypes from 47 males of 112 individuals of the Eurasian Continent and Japan. In addition, 39 mtDNA haplotypes were identified from those animals. From the phylogeography of both the uniparentally inherited genes, four lineages were recognized as Japanese, eastern Eurasian, Caucasian, and western Eurasian. The distribution patterns of the mtDNA lineages showed the existence of a sympatric zone between the eastern and western Eurasian lineages around the Volga River in western Russia. Furthermore, the present study suggested that in the Japanese badgers, the larger genetic differentiation of the Shikoku population was attributable to geographic history in the Japanese islands.

Population structures of the red fox (Vulpes vulpes) on the Hokkaido Island, Japan, revealed by microsatellite analysis.

In order to examine the population structures of the red fox (Vulpes vulpes) on the Hokkaido Island in Japan, we conducted analysis on 250 foxes from all over the island for 12 microsatellite loci. Assignment tests using the genotype data set showed that they were divided into 6 subpopulations. Of the 6, one was geographically isolated in the southern region and considered definitive subpopulation, whereas the other 5 were not. The slight differences among the latter 5 subpopulations were explained by the high adaptability and long dispersal of the red fox on the Hokkaido Island. Although there are few ecological data to explain the genetic differentiation of the southern population, we have proposed some hypotheses from the present ecological and geohistorical viewpoints. One convincing reason from the ecological viewpoint is the restriction of gene flow to southern Hokkaido from other areas due to geographical isolation resulting from the land shape. The other explanation is the geohistorical division of southern Hokkaido from other regions on the island during the last interglacial age, resulting in the isolation of the fox population.