Species Identification of Pinus Pollen Found in Belukha Glacier , Russian Altai Mountains , Using a Whole-Genome Amplification Method

Pollen taxa in sediment samples can be identified based on morphology. However, closely related species do not differ substantially in pollen morphology, and accurate identification is generally limited to genera or families. Because many pollen grains in glaciers contain protoplasm, genetic information obtained from pollen grains should enable the identification of plant taxa at the species level. In the present study, species identification of Pinus pollen grains was attempted using whole-genome amplification (WGA). We used pollen grains extracted from surface snow (depth, 1.8–1.9 m) from the Belukha glacier in the summer of 2003. WGA was performed using a single pollen grain. Some regions of the chloroplast genome were amplified by PCR, and the DNA products were sequenced to identify the pollen grain. Pinus includes approximately 111 recognized species in two subgenera, four sections, and 11 subsections. The tree species Pinus sibirica and P. sylvestris are currently found at the periphery of the glacier. We identified the pollen grains from the Belukha glacier to the level of section or subsection to which P. sibirica and P. sylvestris belong. Moreover, we specifically identified two pollen grains as P. sibirica or P. cembra. Fifteen species, including P. sibirica, were candidates for the remaining pollen grain.


Introduction
The natural range of the genus Pinus is confined to the Northern Hemisphere, although some species have been introduced to the Southern Hemisphere.The genus is currently a dominant component in forests [1].The pollen grains of Pinus have two sacs and the winged pollen grains can be transported long distances by wind.In fact, Pinus pollen grains can reach remote areas, such as mountain glaciers in the Northern Hemisphere, Arctic glaciers, and the Greenland ice sheet, and are found in snow and ice as a predominant pollen type [2][3][4][5][6][7][8][9].Many pollen grains in glaciers are expected to contain protoplasm, and their maintenance at below 0 • C is favorable for DNA preservation [10].This property is characteristic of glacial pollen; protoplasm is rarely seen in pollen found in other sediment types, such as peat and lacustrine deposits.Therefore, pollen grains in glaciers are advantageous for obtaining genetic information, which should enable identification to the species level.Modern pollen identification focuses on the morphological characteristics of the pollen wall, but this approach is generally limited to the identification of plant genera or families.In the case of Pinus, identification to the genus level is typically possible, although the haploxylon and diploxylon types are sometimes distinguished for Pinus pollen grains on the basis of vesicle morphology and other characters.Hence, alternative techniques are needed for species identification, such as DNA analyses of pollen grains.If Pinus pollen grains in glaciers are identified to the species level, it may be possible to investigate the provenance and transportation routes from source plants.Genetic data also provide valuable information related to physiological ecology, gene flow, and population dynamics [11,12].In addition, if genetic analysis is applied to pollen from ice cores, which are cylindrical samples of ice drilled from glaciers and are used to reconstruct past climate conditions and the environmental history of a particular area, then the abovementioned provenance and ecological studies can be extended to trace back into the past.
Nakazawa et al. [10] analyzed the DNA of Pinus pollen grains collected from subsurface snow layers on the Belukha glacier in the Altai Mountains of Russia in the summer of 2003.They identified Pinus pollen grains to the section level using a PCR-based method.However, it is difficult to achieve species-level identification using this approach.The sequences provide limited information owing to their short lengths, meaning that there is the potential to obtain more DNA information by improving on PCR-based methods.
To identify the grains at the species level, an optimized whole-genome amplification (WGA) approach combined with multiplex PCR was developed in this study.Multiplex PCR can be used to amplify multiple targets in a single PCR experiment.Additionally, DNA barcoding is used to identify materials from known species based on short DNA sequences of standard genomic regions (i.e., DNA barcodes).In general, chloroplast DNA in land plants has a low nucleotide substitution rate, on the order of 10 −9 per site per year [13].Therefore, few mutations are expected within a short period of time, and the most promising DNA barcoding loci for plants are chloroplast genes.DNA barcoding technologies are being developed for applications in palynology (pollen DNA barcoding) [14], and these studies have demonstrated that both chloroplast and nuclear barcoding markers can be amplified from pollen.Unfortunately, most plastid candidate barcodes lack species-level resolution [15].Additionally, DNA barcoding markers with universal primers used in previous pollen DNA barcoding studies [14] provide insufficient information on the species-level taxonomy of Pinus.Species-specific primers are generally used for the precise identification of samples [16,17].Moreover, PCR amplification using limited amounts of DNA template, such as from a single pollen grain, has high risk of contamination, biased amplification, and product redundancy [18,19].Thus, in this study, we designed primers for species-specific DNA barcoding with high resolution at low taxonomic levels to reduce these risks and increase precision.

Study Area and Pollen Samples
The Belukha glacier (49 • 49 N, 86 • 34 E; 4,110 m a.s.l.) is located on the western side of Mt.Belukha (4500 m a.s.l.) in the Russian Altai Mountains and is situated in the border region between Russia, Mongolia, China, and Kazakhstan (Figure 1).In the summer of 2003, a 4-m-deep pit on the plateau of the glacier was examined (4100 m a.s.l.) [20].Pinus pollen grains in a snow sample were obtained at a depth of 1.8-1.9m in the pit and used for the DNA analysis.The Pinus pollen concentration in the sample was 34,900 grains L −1 .The sample was dated to the summer of 2002 by counting the seasonal distribution of pollen [20].The sample was obtained from the same pit as in our previous study [10], but the pollen grains were previously collected from a depth of 0.4-0.5 m and dated to the summer of 2003.The sample was kept in a frozen state until it was analyzed.
Pinus species surrounding the Belukha glacier are P. sibirica, which is distributed between approximately 1000 and 2000 m a.s.l., and P. sylvestris, which typically occurs below the P. sibirica stands [21].For a detailed description of the vegetation, see Nakazawa et al. [10].The sampling site was located on the western side of Mt.Belukha.The distribution map was compiled based on the maps of Farjon [22].The figure is from our previous study [10].

DNA Extraction from a Single Pollen Grain
DNA was extracted using a modified version of the extraction method described by Nakazawa et al. [10].A flow chart of the experimental procedure is shown in Figure 2. Melted snow and ice samples were first filtered through a hydrophilic PTFE membrane filter with a pore size of 10 μm.Next, pollen grains that showed no structural damage were selected from the filter using a micromanipulator (MM-88; Narishige, Tokyo, Japan) under a microscope.Each selected pollen grain was placed onto a new hydrophilic PTFE membrane filter with a pore size of 5 μm, and was then washed by suction filtration with 1 mL of nuclease-free water (Ambion Life Technologies, Foster City, CA, USA).The filter trapping a pollen grain was transferred to a sterile Petri dish.The washed grain was then transferred to the inner side of the lid of a DNA-free PCR tube containing 0.5 μL of water using a pipette.The sampling site was located on the western side of Mt.Belukha.The distribution map was compiled based on the maps of Farjon [22].The figure is from our previous study [10].

DNA Extraction from a Single Pollen Grain
DNA was extracted using a modified version of the extraction method described by Nakazawa et al. [10].A flow chart of the experimental procedure is shown in Figure 2. Melted snow and ice samples were first filtered through a hydrophilic PTFE membrane filter with a pore size of 10 µm.Next, pollen grains that showed no structural damage were selected from the filter using a micromanipulator (MM-88; Narishige, Tokyo, Japan) under a microscope.Each selected pollen grain was placed onto a new hydrophilic PTFE membrane filter with a pore size of 5 µm, and was then washed by suction filtration with 1 mL of nuclease-free water (Ambion Life Technologies, Foster City, CA, USA).The filter trapping a pollen grain was transferred to a sterile Petri dish.The washed grain was then transferred to the inner side of the lid of a DNA-free PCR tube containing 0.5 µL of water using a pipette.The sampling site was located on the western side of Mt.Belukha.The distribution map was compiled based on the maps of Farjon [22].The figure is from our previous study [10].

DNA Extraction from a Single Pollen Grain
DNA was extracted using a modified version of the extraction method described by Nakazawa et al. [10].A flow chart of the experimental procedure is shown in Figure 2. Melted snow and ice samples were first filtered through a hydrophilic PTFE membrane filter with a pore size of 10 μm.Next, pollen grains that showed no structural damage were selected from the filter using a micromanipulator (MM-88; Narishige, Tokyo, Japan) under a microscope.Each selected pollen grain was placed onto a new hydrophilic PTFE membrane filter with a pore size of 5 μm, and was then washed by suction filtration with 1 mL of nuclease-free water (Ambion Life Technologies, Foster City, CA, USA).The filter trapping a pollen grain was transferred to a sterile Petri dish.The washed grain was then transferred to the inner side of the lid of a DNA-free PCR tube containing 0.5 μL of water using a pipette.The pollen grain in the lid was treated with endonuclease and exonuclease to eliminate potential contaminants from DNA fragments attached to the surface of the grain.One microliter of reaction mixture containing 0.48 µL of water, 0.2 µL of 10× Exonuclease III Buffer (TaKaRa Biotechnology Co., Ltd., Dalian, Japan), 0.1 µL of Exonuclease I (20 U/µL, Epicentre, Madison, WI, USA), 0.02 µL of Exonuclease III (200 U/µL, TaKaRa Biotechnology), and 0.2 µL of DNase I (1 U/µL, Sigma-Aldrich Co., St. Louis, MO, USA) was added to the lid.The mixture was incubated at 37 • C for 3-5 h, then at 98 • C for 10 min.
The treated grain was crushed directly in the lid of the tube using a sterile plastic pipette tip.For crushing, we used an electric toothbrush with an attached tip-rounded pipette tip, which was made by heating with a gas burner in advance.The vibration of the electronic toothbrush facilitated crushing.One microliter of extraction mixture containing 0.7 µL of Tris-HCl (pH 8.0, 20 mM, Nacalai Tesque, Kyoto, Japan), 0.2 µL of proteinase K (1 µg/µL, TaKaRa Biotechnology), and 0.1 µL of cellulase (Sigma-Aldrich) was added to the lid of each sample and spun down for collection at the bottom of the tube.The mixture was incubated at 50 • C for 6 h, at 95 • C for 10 min, and then used as a template.

Whole-Genome Amplification
WGA from a pollen grain was performed using the GenomePlex Single Cell Whole Genome Amplification Kit (WGA4; Sigma-Aldrich) according to the manufacturer's instructions, with slight modifications.After the lysis procedure, 1.3 µL of fragmentation solution, including 0.3 µL of fragmentation buffer and 1 µL of water, was added and heated in a thermal cycler (Bioer Technology Co. Ltd., Hangzhou, China) at 99 • C for 3 min.The PCR tubes were then cooled on a cooling rack (Nippon Genetics Co., Ltd., Tokyo, Japan).For library preparation, 1 µL of the reaction solution including 0.7 µL of Library Preparation Buffer and 0.3 µL of Library Stabilization Solution was added to each sample and was placed in the thermal cycler at 95 • C for 2 min.The samples were cooled on the cooling rack.Next, 1.0 µL of enzyme solution including 0.7 µL of water and 0.3 µL of the Library Preparation Enzyme solution was added to each sample.The samples were placed in the thermal cycler.The reaction time and temperature were based on the instructions provided.
For the amplification, 14.6 µL of the reaction mixture, including 11.0 µL of water, 2.0 µL of 10× Amplification Master Mix, 1.3 µL of WGA DNA Polymerase, 0.2 µL of uracil-N-glycosylase (UNG; TaKaRa Biotechnology), and 0.1 µL of BIOTAQ HS DNA polymerase (Bioline, London, UK), was added to the sample.UNG was used to degrade uracil-containing PCR contaminants from previous PCR prior to the amplification reaction.This is explained in more detail in the next section.Note that the PCR products included dUTP, instead of dTTP.Therefore, the treatment with UNG should allow the selective removal of carryover PCR products.Samples were first incubated at 25 • C for 10 min, and then amplified using an initial denaturation of 95 • C for 10 min followed by 40 cycles each consisting of a denaturation step at 94 • C for 30 s, an annealing step at 52 • C for 1 min, and an extension step at 72 • C for 1 min, and a final extension at 72 • C for 7 min.The WGA DNA in the reaction mixture was stored at −20 • C until further use, without a DNA purification step.

UNG Treatment and Multiplex PCR Amplification
The quality of the WGA amplification was evaluated using a multiplex PCR approach.The DNA specimens generated by WGA were subjected to various analyses for chloroplast DNA.Instead of single PCR, a multiplex PCR step was introduced in this study to use the specimens effectively.To preclude carryover contamination of amplification products from previous PCR, the reaction with UNG, an enzyme that degrades uracil-containing DNA, was carried out.Each multiplex PCR amplification was performed with dUTP instead of dTTP.Thus, the PCR products should be obtained from only thymine-containing templates amplified by WGA.
For the multiplex PCR assay, 2-7 primer pairs were designed.PCR was carried out with a 10-µL reaction mixture containing 1 µL of template WGA products.The primers are listed in Table 1.
To identify samples at the section or subsection level, a total of 10 primer pairs were used.Seven or 14 primer pairs were selected to narrow the candidates within each subsection.To evaluate the primer sets, the performance of primers was examined using single pollen grains of P. resinosa belonging to subsection Pinus and needles collected from P. pumila, P. strobus, P. taeda, P. jeffreyi, and P. monophylla, which belong to subsections Strobus, Strobus, Australes, Ponderosae and Cembroides, and P. resinosa, respectively.From the results, some primers were screened out.The primers presented in Table 1 were effective for sequencing.In addition, the number of cycles of multiplex PCR in this study was minimized to avoid the introduction of significant PCR bias.

Nested PCR
The secondary amplification for each strand, which was run with a nested set of primers, was performed using a 0.5-µL aliquot of the purified first PCR products.A new PCR mixture of 10 µL was prepared, containing 0.5 µL of the template, 4.75 µL of water, 2.0 µL of 5× Green GoTaq Flexi Buffer (Promega Co., Madison, WI, USA), 0.9 µL of 25 mM MgCl 2 , 0.8 µL of dU plus dNTP Mixture, 1.0 µL of 5.0 µM primers, and 0.05 µL of GoTaq Hot Start Polymerase (Promega).The primers used for nested PCR are listed in Table 1.The amplification was performed using a thermal cycler (GeneAmp PCR System 9700; Applied Biosystems, Foster City, CA, USA) under the following conditions: initial activation at 95 • C for 2 min, 20 cycles of denaturation at 95 • C for 30 s, annealing at 52 • C for 30 s, and extension at 72 • C for 30 s, followed by a final incubation at 72 • C for 5 min.Amplified PCR products were then sequenced using the BigDye Terminator v.3.1 Sequencing Kit (Applied Biosystems, Foster City, CA, USA) and an ABI 3130xl Genetic Analyzer (Applied Biosystems).

Identification of Individual Pinus Pollen Grains
Each Pinus pollen grain was identified based on parsimony-informative characters in various regions of the chloroplast genome, as shown in Table 1.The genus Pinus has approximately 111 recognized species in two subgenera, four sections, and 17 subsections (Table 2).Sequence data were collected for these regions from GenBank, which are available for almost all Pinus species.In addition, the aligned DNA sequences for each locus were used to determine parsimony-informative sites in MEGA ver. 5 [23] (see Supplementary data).Additionally, a small insertion of 3 bp was considered for the Matsu5 (ycf3) and Matsu7 (rpoB) regions.Genus-level classification was based on Gernandt et al. [24].Their classification, based on a chloroplast DNA phylogeny, was a modification of (1) the influential classification of Little and Critchfield [25], which is based primarily on morphological characters and data from interspecific crosses, and (2) the classification of Price et al. [26], which incorporates additional recently described species.

Table 2.
Classification system for the genus Pinus.The classification system is based on Gernandt et al. [24].The colors indicating various taxa correspond to those in Table 3.

Pinus Pollen Identification at the Section or Subsection Level
We analyzed 21 pollen grains, and five samples showed positive amplification from at least one locus in the multiplex PCR with the primer sets Matsu1-10 (Table 3).DNA fragmentation and degradation, particularly in ancient samples, make it difficult to amplify long fragments from a single pollen grain in sediment samples [27,28].This problem is alleviated by the higher amplification efficiency of short fragments (<200 bp) [29].To increase the probability of amplification, our primer sets were designed to yield fragments of around 200 bp, even though our samples were not ancient.Although previous pollen DNA barcoding studies have focused on fragments of longer than 300 bp [14], we believe that fragments of around 200 bp may be sufficient to identify pollen at the section or subsection level when the target pollen type is specified and sequence data are available, as demonstrated in the present study.
The controls in our PCRs followed the experimental methodology of Parducci et al. [29].A positive control was not used owing to the high contamination risk.As a negative control, contamination by exogenous chloroplast DNA from Pinus species in the reagents was monitored using a PCR blank that included all reagents.If no band was visible by agarose gel electrophoresis after multiplex PCRs, we concluded that the samples were not contaminated.Samples with bands on the agarose gel were selected for an additional nested PCR amplification step with each single primer set.
Based on the sequence data collected for the clpP, rpoA, atpB, rpoC1, ycf3, ycf4, and rpoB regions by nested PCR, the pollen grains were identified to the section or subsection level (Table 3).We identified four out of five pollen grains as members of the subsection Strobus.One belonged to the section Quinquefoliae, and the remaining grain was part of the section Pinus.
For the four grains in the subsection Strobus, we successfully obtained sequence data for multiple loci (Table 3).Based on these sequence data, we identified each locus as belonging to the same phylogenetic clade, namely, subsection Strobus of the section Quinquefoliae in the subgenus Strobus (Table 2).This consistency confirms that the WGA reactions could increase the amount of accurate DNA data that can be obtained from a single pollen grain.

Species Identification of Pinus Pollen Grains
To further narrow down the candidate species in the subsection Strobus, we used the WGA products for multiplex PCRs with primers to amplify short fragments on the plastid gene ycf1.We designed various primer sets targeting regions in the ycf1 gene (Table 1), which is a highly variable locus that can be used for Pinus species identification [30].The multiplex PCR products were subjected to an additional round of nested PCR.
We successfully amplified target regions for three out of four pollen grains in the subsection Strobus.Regions that were successfully amplified are shown in Table S1.In the sequence analysis based on variable base positions, we identified two pollen grains as P. cembra or P. sibirica (samples No. 1 and 3).Although we examined various regions in the chloroplast genome, to positively identify pollen grains as P. sibirica or P. cembra, it was sufficient to examine only three regions, namely, Strbs3, Strbs7, and Matsu6 (or Matsu8).For the remaining grain, we determined the following 15 candidate species: P. albicaulis, P. armandii, P. bhutanica, P. cembra, P. dalatensis, P. fenzeliana, P. lambertiana, P. monticola, P. morrisonicola, P. parviflora, P. peuce, P. pumila, P. sibirica, P. wallichiana, and P. wangii.
We inferred that the three grains were P. sibirica.Gugerli et al. [31] observed highly similar chloroplast and mitochondrial DNA sequences between P. cembra and P. sibirica.This similarity suggests a relatively recent evolutionary separation of the species, despite their currently disjunct distributions (Figure 1).P. sibirica appears in the Urals and Siberia of Russia, in eastern Kazakhstan, in northern Mongolia, and in Xingjiang, Nei Mongo, and Heilongjiang of China.P. cembra occurs in the Swiss Alps, in the Tirol of Austria, in the High Tatra between Poland and the Slovak Republic, and in the eastern Carpathians of Romania and Ukrine [22].Heinze and Holzer [32] verified that a nearly complete P. cembra chloroplast genome sequence in GenBank (Accession No. FJ899574) is identical to a P. sibirica sequence (Accession No. FJ899558).Although the FJ899574 sequence lacks part of the ycf1 region targeted in this study, the ycf1 sequence identities for both species were validated by a comparative analysis using FJ899558 and KP128626 of P. cembra.Accordingly, the two species likely cannot be discriminated based on comparative sequence analyses of the chloroplast genome.P. sibirica is an extant species and is the only member of the subsection Strobus that is found near the glacier.This species was a major candidate in our study, suggesting that the pollen grains in the glacier originated from P. sibirica trees found in the immediate surroundings.This consistency between the results of our genetic analysis and species distribution data suggests that our method to identify pollen species was reliable.
We were not able to obtain sequence data for three grains.For the pollen grain identified as belonging to the section Pinus, the WGA products were used as templates in the multiplex PCRs, with primer sets for the ycf1, rbcL, and rpl20-rps18 chloroplast regions (Table 1).We could not identify those grains at a lower taxonomic level.However, the grains appeared to be P. sibirica in subsection Strobus, section Quinquefoliae, and P. sylvestris, which belongs to subsection Pinus in section Pinus based on the consistent results at the subsection or section levels.To obtain accurate sequences and facilitate a more detailed taxonomic identification, additional PCRs with newly designed primer sets may be effective.Although we were able to generate sufficient DNA specimens from single pollen grains for additional PCRs using the WGA technique, we did not perform subsequent analyses owing to a lack of time and resources and an expectation of unremarkable results.

Potential Use of Pollen Grains as a Tracer for Emission Sources
We believe that our method for Pinus pollen identification is suitable for further work aimed at a more detailed characterization of the provenance of aerosols, particularly in Arctic glaciers and the Greenland ice sheet.Pollen is a type of bio-aerosol; Pinus pollen as well as other types of pollen travel long distances.This pollen is regarded as exotic, and many studies have investigated its source area and long-distance transportation by trajectory analyses [33][34][35][36][37].For aerosols reaching ice sheets and mountain glaciers, dust has been used as a tracer for emission sources as well as large-scale atmospheric circulation [38][39][40].Clay mineralogy and Sr-Nd isotopic and elemental compositions have suggested that East Asia (i.e., the Gobi Desert, northern Chinese deserts, and the Taklamakan Desert) is the main source of dust arriving in Greenland, both at present and during the last glacial period [38].Generally, dust seems to originate mainly from arid regions.In contrast, pollen sources are restricted to vegetated areas.Therefore, pollen can be used as another tracer, and a combination of both of these tracers should lead to a better understanding of the provenance of solid aerosols and the materials cycle.This approach has not been used, although some pollen grains, including Pinus pollen grains, have been found in Greenland [2,35,36].
As mentioned in Section 3.2, our analysis strongly suggested that two Pinus pollen grains found in the Belukha glacier are P. sibirica, in consonance with the surrounding Pinus vegetation.Therefore, we can assume that the provenance of the Pinus pollen is the region that extends from the northwest to the east of the glacier, as shown in Figure 1.A Pinus pollen grain is around 50 µm in size; the size of a pollen grain typically ranges from 10 to 200 µm, and the most common size is between 20 and 60 µm [41].Despite the relatively large size of Pinus pollen, the grains are well dispersed in a vesiculate form with two prominent sacci [42].In addition, Pinus is characterized by high pollen grain production.For those reasons, Pinus pollen in palaeoecological records are frequently disregarded as evidence for presence or absence in the arctic region with low local production of pollen [1].However, these properties are favorable for analyses of Pinus pollen grains from Arctic glaciers and the Greenland ice sheet by our method to investigate geographic provenances and the materials cycle.Additionally, their large size seems to be beneficial for the treatment of pollen grains in a laboratory.

Improvement in the Amplification Success Rate
The success rate for obtaining sequence data reported by Nakazawa et al. [10] was 7.6% (n = 105).This was higher than the success rates observed in previous DNA analyses of pollen collected from peat or lacustrine deposits.Suyama et al. [28,43] and Parducci et al. [29] observed success rates ranging from 0 to 3.2%.However, these rates are still insufficient to develop a new field of palynological research based on genetic information.Hence, improving the success rate is a particularly important issue.
In this study, we subjected 21 pollen grains to WGA, and we observed positive amplification from at least one locus for a total of five grains.The success rate for sequence amplification in this study was 24% and exceeded that of Nakazawa et al. [10] who used pollen samples from the same glacier collected from the upper layer of the pit.To improve the success rate of DNA analyses, we introduced multiplex PCR; amplifying multiple loci simultaneously in a single reaction improved the probability that at least one locus was amplified.In addition, the WGA technique enabled an increased quantity of DNA to be obtained from samples with limited DNA content.Since we collected samples from a different layer from that of previous studies, we cannot make a simple comparison of the success rates between studies.However, we were able to obtain sequence data from multiple loci, and this method appeared to be more effective with respect to success rate.

Conclusions
We described an initial attempt to identify pollen grains from a glacier at the species level based on chloroplast DNA sequences.For precise identification, we applied an optimized WGA technique for single pollen grains.We subjected WGA products to an additional round of multiplex PCRs.The combined DNA sequence data obtained from a single pollen grain suggested that identification at or near the species level is possible.
We analyzed 21 pollen grains, of which five exhibited positive amplification for at least one locus in the multiplex PCRs.One grain appeared to belong to the section Quinquefoliae, four grains were in the subsection Strobus of section Quinquefoliae, and the remaining grain was identified as a member of the section Pinus.In addition, for three out of the four grains in the subsection Strobus, the candidates were narrowed down to two species; 15 candidates remained for the other grain.Owing to the identical P. cembra and P. sibirica chloroplast genome sequences, it was not possible to differentiate between the species using sequence data.However, we inferred that both grains were P. sibirica because it is an extant species that is currently distributed around the glacier.Meanwhile, P. cembra is distributed in some high mountains in Switzerland, Austria, Poland, the Slovak Republic, Romania, and Ukraine [22].Moreover, we could assume the grains traveled from the region that expands from the northwest to the east of the glacier based on the forest distribution of P. sibirica.
Similarly, the remaining grains appeared to be P. sibirica or P. sylvestris, which is also found around the glacier.Our pollen identifications based on DNA sequence data are supported by the vegetation in the region, suggesting that the method established in this study enables reliable identification at the species level for pollen grains.Additionally, the method should be useful for future studies on the provenance of solid aerosols and the materials cycle in the polar region and high mountain regions in the Northern Hemisphere, where glaciers exist.
The DNA amplification success rate in this study exceeded that of our previous study of samples from the same snow pit.However, the samples were collected from a different layer, preventing a simple comparison of success rates between studies.However, we demonstrated that our method could be used to obtain sequence data from multiple loci and effectively increase the success rate.
Further investigations using older samples from ice cores are necessary to extend the applications of these methods and to accumulate data in the field.The rarity of suitable, well-preserved pollen samples in sediments has limited the broad utility of DNA studies for the taxonomic identification of pollen.However, owing to low-temperature conditions, pollen grains in glaciers are not strongly affected by diagenesis and their DNA is therefore likely to be preserved.Accordingly, pollen samples from glaciers have broad utility for studies of taxonomy, past vegetation, population genetics, and climate and environmental conditions.Our method based on WGA and multiplex PCR techniques may enable the generation of DNA specimens from single pollen grains that can be analyzed by multiple molecular techniques for a range of applications.

Figure 1 .
Figure 1.Location of the Belukha glacier in Russia's Altai Republic, and distributions of extant P. cembra and two Pinus species (P.sibirica and P. sylvestris) found in the region surrounding the glacier.The sampling site was located on the western side of Mt.Belukha.The distribution map was compiled based on the maps of Farjon[22].The figure is from our previous study[10].

Figure 1 .
Figure 1.Location of the Belukha glacier in Russia's Altai Republic, and distributions of extant P. cembra and two Pinus species (P.sibirica and P. sylvestris) found in the region surrounding the glacier.The sampling site was located on the western side of Mt.Belukha.The distribution map was compiled based on the maps of Farjon[22].The figure is from our previous study[10].

Figure 1 .
Figure 1.Location of the Belukha glacier in Russia's Altai Republic, and distributions of extant P. cembra and two Pinus species (P.sibirica and P. sylvestris) found in the region surrounding the glacier.The sampling site was located on the western side of Mt.Belukha.The distribution map was compiled based on the maps of Farjon[22].The figure is from our previous study[10].

Figure 2 .
Figure 2. Flow chart of the experimental procedure.

Table 1 .
Primers for the target fragments used in the present study.The primers for nested PCR for Pns1-7 are not shown because the multiplex PCR products showed no visible band by agarose gel electrophoresis, and an additional round of nested PCR was not performed.

Table 3 .
Identification of pollen grains by multiplex PCR.Each pollen grain was identified at the subgenus, section, or subsection level based on regions that were positively amplified.The S, Q, P, and St indicate the subgenus Strobus, section Quinquefoliae, section Pinus, and subsection Strobus, respectively.Dashes indicate a lack of amplification.For the sequence data and identification procedures, see Supplementary data.