Chloroplast genome-based genetic resources via genome skimming for the subalpine forests of Japan and adjacent regions.

The Japanese subalpine zone is dominated by an ecologically important forest biome, subalpine coniferous forest, constituting a distinct assemblage of cold-tolerant angiosperm and conifer species. While being relatively intact compared to other forest biomes in Japan, subalpine coniferous forests are under significant threat from deer browsing, global warming and small population size effects. However, there is a severe lack of genetic resources available for this biome's major constituent plant species. This study aimed to develop chloroplast genome-based genetic resources for 12 widespread subalpine tree and shrub species (7 angiosperms and 5 conifers) via genome skimming of whole-genomic DNA using short reads (100-150 bp in length). For 10 species, whole chloroplast genomes were assembled via de novo-based methods from 4 to 10 individuals per species sampled from across their ranges in Japan and, for non-Japanese endemic species, elsewhere in northeast Asia. A total of 566 single nucleotide polymorphisms for Japanese samples and 768 for all samples (varying from 2 to 202 per species) were identified which were distributed in geographically restricted lineages in most species. In addition, between 9 and 58 polymorphic simple sequence repeat regions were identified per species. For two Ericaceae species (Rhododendron brachycarpum and Vaccinium vitis-idaea) characterised by large chloroplast genomes, de novo assembly failed, but single nucleotide polymorphisms could be identified using reference mapping. These data will be useful for genetic studies of species taxonomic relationships, investigating phylogeographic patterns within species, developing chloroplast-based markers for conservation genetic studies and has potential application for studies of environmental and ancient DNA.

Phylogenetic origin of two Japanese Torreya taxa found in two regions with strongly contrasting snow depth.

The Japanese archipelago exhibits a notable difference in snow depth in winter, deep snow on the Sea of Japan side and low snow cover on the Pacific Ocean side. This contrasting pattern has shaped the distribution of infraspecific taxon pairs in a range of woody plants, with taxa found on the Sea of Japan side typically exhibiting a stunted shrub form with multiple decumbent stems. The phylogenetic origin of these taxon pairs is unknown, i.e., whether the two taxa diverged from the same species or if they have different origins. This study aimed to reveal the phylogenetic origin of two varieties of Torreya nucifera (Taxaceae); var. nucifera is a tree found on the Pacific Ocean side, whereas var. radicans is a shrub found on the Sea of Japan side. We examined the phylogenetic relationships of the two varieties and worldwide Torreya taxa using whole chloroplast genomes, chloroplast DNA fragments, and the nuclear ribosomal internal transcribed spacer (ITS). The whole chloroplast genome phylogeny indicated that T. nucifera var. radicans was a sister taxon to Chinese T. grandis, rather than to var. nucifera. In contrast, the nuclear ITS phylogeny indicated that while several haplotypes of T. nucifera var. radicans were closely related to T. grandis, most haplotypes of T. nucifera var. radicans formed a single clade with those of var. nucifera. This implies that the homogenization of the ITS has occurred between the two taxa, while taxon-specific chloroplast DNA haplotypes were retained. These discordant phylogenies suggested that the two taxa have different phylogenetic origins, but have an intricate evolutionary history, involving inter-taxa hybridization and gene flow, possibly when their distributions were confined to sympatric refugia. Given the genetic evidence and distinct difference in growth form, we propose that T. nucifera var. radicans should be taxonomically treated as a distinct species, T. fruticosa.

Sika deer presence affects the host-parasite interface of a Japanese land leech.

Since the 1990s, increasing populations of a blood feeding land leech (Haemadipsa japonica) have become a serious issue in several Japanese prefectures, and it may be caused by the increases in sika deer (Cervus nippon) populations seen over the last quarter of the century. Therefore, this study aimed to reveal the host animal species of H. japonica using iDNA (vertebrate DNA isolated from invertebrates) and to test the hypothesis that the increasingly widespread distribution of sika deer results in increased H. japonica populations through changes to the host-parasite interface. We amplified mitochondrial DNA 16S ribosome RNA fragments from iDNA isolated from the blood clots of H. japonica collected across Japan. We identified 17 host animal species, including four orders of Mammalia (Carnivora, Artiodactyla, Rodentia, and Lagomorpha) and two orders of Amphibia (Caudata and Anura). The sika deer was the dominant host species of H. japonica. Additionally, the host animal species composition of H. japonica differed according to the presence or absence of sika deer. In the sites where sika deer were not found, Anura (frog) species were the most commonly identified hosts of H. japonica. These results suggest that the increases in H. japonica populations might have occurred via a change in host preference to sika deer. This change might be driven by the increases in sika deer populations and subsequent increase in the frequency that H. japonica uses the sika deer as easy prey, as well as by sika deer providing more reproductive energy per blood meal than blood meal from frog species. The present study suggests that a more widespread distribution of sika deer resulted in an increase in H. japonica through a change in the host-parasite interface. Therefore, management that focuses on decreasing sika deer populations would likely be an effective method for the reduction of H. japonica populations.

Molecular monitoring of Neoseiulus californicus released from sheltered slow-release sachets for spider mite control in a Japanese pear greenhouse.

A novel system for spider mite control was developed with a slow-release sachet containing Neoseiulus californicus (McGregor) (Acari: Phytoseiidae) protected by a waterproof shelter. Monitoring the efficacy of the predator release system for spider mite control at a Japanese pear greenhouse requires discrimination of N. californicus from other indigenous phytoseiid mite species inhabiting the study site and subsequent identification of the released N. californicus. The report of our earlier study described a PCR-based method for discrimination of N. californicus species. For the present study, we first examined phytoseiid mite species composition in the greenhouse. Subsequently, we developed microsatellite markers to identify the released N. californicus. Finally, we installed the predator release system in the greenhouse and conducted a population survey of phytoseiid and spider mites. Results demonstrated that approximately 1 month is necessary for distribution of the released N. californicus on the leaves.

Nuclear microsatellite and mitochondrial DNA analyses reveal the regional genetic structure and phylogeographical history of a sanguivorous land leech, Haemadipsa japonica, in Japan.

Recent molecular studies have indicated that phylogeographical history of Japanese biota is likely shaped by geohistory along with biological events, such as distribution shifts, isolation, and divergence of populations. However, the genetic structure and phylogeographical history of terrestrial Annelida species, including leech species, are poorly understood. Therefore, we aimed to understand the genetic structure and phylogeographical history across the natural range of Haemadipsa japonica, a sanguivorous land leech species endemic to Japan, by using nine polymorphic nuclear microsatellites (nSSR) and cytochrome oxidase subunit one (COI) sequences of mitochondrial DNA (mtDNA). Analyses using nSSR revealed that H. japonica exhibited a stronger regional genetic differentiation among populations (G'ST = 0.77) than other animal species, probably because of the low mobility of land leech. Analyses using mtDNA indicated that H. japonica exhibited two distinct lineages (A and B), which were estimated to have diverged in the middle Pleistocene and probably because of range fragmentation resulting from climatic change and glacial and interglacial cycles. Lineage A was widely distributed across Japan, and lineage B was found in southwestern Japan. Analyses using nSSR revealed that lineage A was roughly divided into two population groups (i.e., northeastern and southwestern Japan); these analyses also revealed a gradual decrease in genetic diversity with increasing latitude in lineage A and a strong genetic drift in populations of northeastern Japan. Combined with the largely unresolved shallow polytomies from the mtDNA phylogeny, these results implied that lineage A may have undergone a rapid northward migration, probably during the Holocene. Then, the regional genetic structure with local unique gene pools may have been formed within each lineage because of the low mobility of this leech species.

An admixture of Quercus dentata in the coastal ecotype of Q. mongolica var. crispula in northern Hokkaido and genetic and environmental effects on their traits.

In northern Japan, coastal oak forests consist of Quercus dentata (Qd) on the coastal side and Q. mongolica var. crispula (Qc) on the inland side. In the forests of northern Hokkaido, Qd is rare, and a coastal ecotype of Qc with some Qd-like traits grows on the coastal side. To reveal the genetic background of this ecotype, nuclear microsatellite genotypes in closely related oak taxa were obtained from the Eurasian continent, Sakhalin, and Hokkaido. The clustering of these genotypes suggests an admixture of Qd in the coastal ecotype of Qc. Next, we evaluated the effects of admixture and coastal stress on the leaf and shoot traits of Qc and Qd along coastal-inland gradients in northern Hokkaido. The admixture of Qd in Qc was quantified by the Qd ancestry proportions. Coastal stress causes bud mortality in the upper parts of shoots and was quantified by the survival patterns of buds in shoots. The genetic and environmental effects on the traits at Qd-abundant and Qd-rare sites were estimated using linear mixed models. The genetic effect was detected in all traits. Both genetic and environmental effects were detected in most traits. Some traits differed between Qd-abundant and Qd-rare sites in addition to these effects, indicating more Qd-like traits at Qd-rare sites. The findings suggest that an admixture of Qd characterizes the genetic background of the coastal ecotype of Qc and that not only the coastal stress but also the genetic background is responsible for the leaf and shoot traits of Qc and Qd in northern Hokkaido.

Characterization of 13 polymorphic microsatellite loci in the Japanese land leech.

We developed 13 polymorphic microsatellite loci of the Japanese land leech (Haemadipsa japonica; Haemadipsidea) using an Illumina MiSeq sequencing approach. A total of 42,064 nuclear DNA contigs were filtered for microsatellite motifs, among which 30,873 simple sequence repeat loci were identified. From these sequences, we selected 30 primer sets, and 13 of these loci were successfully amplified. Polymorphism of the 13 loci was tested using 16 individuals sampled from sixteen populations across Japan. The number of alleles and polymorphism information content varied from 5 to 17 and 0.335 to 0.883, respectively, and observed and expected heterozygosity values ranged from 0.143 to 0.875 and 0.349 to 0.893, respectively, indicating that these loci are polymorphic. Furthermore, we established useful multiplex PCR using these loci. The 13 microsatellite loci described in this paper are the first nuclear microsatellite markers for a land leech species.

Genetic structure of Sakhalin spruce (Picea glehnii) in northern Japan and adjacent regions revealed by nuclear microsatellites and mitochondrial gene sequences.

The genetic structure of Sakhalin spruce (Picea glehnii) was studied across the natural range of the species, including two small isolated populations in south Sakhalin and Hayachine, by using six microsatellite loci and maternally inherited mitochondrial gene sequences. We also analyzed P. jezoensis, a sympatric spruce in the range. Genetic diversity of P. glehnii was higher in central Hokkaido and the lowest in the Hayachine. Bayesian clustering and principal coordinate analysis by using the microsatellites indicated that the Hayachine was clearly distinct from other populations, implying that it had undergone strong genetic drift since the last glacial period. P. glehnii harbored four mitochondrial haplotypes, two of which were shared with P. jezoensis. One of the two was observed without geographical concentration, suggesting its derivation from ancestral polymorphism. Another was observed in south Sakhalin and in P. jezoensis across Sakhalin. The Bayesian clustering--by using four microsatellite loci, including P. jezoensis populations--indicated unambiguous species delimitation, but with possible admixture of P. jezoensis genes into P. glehnii in south Sakhalin, where P. glehnii is abundantly overwhelmed by P. jezoensis; this might explain the occurrence of introgression of the haplotype of P. jezoensis into P. glehnii.

Phylogeography of the Korean pine (Pinus koraiensis) in northeast Asia: inferences from organelle gene sequences.

Range-wide genetic variation of Korean pine (Pinus koraiensis) was assessed using maternally inherited mtDNA and paternally inherited cpDNA for 16 natural populations throughout northeast Asia in order to study its phylogeographical history during the Quaternary. The cpDNA variation indicated that there was no difference between populations on the Asian continent and those in the Japanese archipelago. In contrast, the mtDNA variation indicated that there was significant difference between the populations from the two regions, with each region having a different lineage. The continental populations exhibited no diversity in the mtDNA examined despite the species' current extensive range and large populations. Conversely, while the Korean pine is rare in Japan, the Japanese populations exhibited greater levels of mtDNA diversity (H (T) = 0.502). The higher mtDNA diversity and evidence from numerous Korean pine macrofossil remains dated to the Pleistocene and recovered various sites in Japan suggest that the Japanese archipelago once served as a refugium to a much larger Korean pine population with a more extensive range than is the case today. The presence of the single mtDNA haplotype across the Asian continent suggests that the present widespread populations could have expanded from a single refugium population after the last glacial periods.