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.

Efficient selection of a biallelic and nonchimeric gene-edited tree using Oxford Nanopore Technologies sequencing.

The clustered regularly interspaced short palindromic repeats (CRISPR)/CRISPR-associated protein 9 nuclease system is a versatile and essential biotechnological tool in the life sciences that allows efficient genome editing. When generating gene-edited trees, T0-generation plants are often used for subsequent analysis because of the time that is required to obtain the desired mutants via crossing. However, T0-generation plants exhibit various unexpected mutations, which emphasizes the need to identify mutants with expected mutation patterns. The two critical checkpoints in this process are to confirm the expected mutation patterns in both alleles and to exclude somatic chimeric plants. In this study, we generated gene-edited Cryptomeria japonica plants and established a method to determine chimerism and mutation patterns using fragment analysis and Oxford Nanopore Technologies (ONT)-based amplicon sequencing. In the first screening, fragment analysis, i.e., indel detection via amplicon analysis, was used to predict indel mutation patterns in both alleles and to discriminate somatic chimeric plants in 188 candidate mutants. In the second screening, we precisely determined the mutation patterns and chimerism in the mutants using ONT-based amplicon sequencing, where confirmation of both alleles can be achieved using allele-specific markers flanking the single guide RNA target site. In the present study, a bioinformatic analysis procedure was developed and provided for the rapid and accurate determination of DNA mutation patterns using ONT-based amplicon sequencing. As ONT amplicon sequencing has a low running cost compared with other long-read analysis methods, such as PacBio, it is a powerful tool in plant genetics and biotechnology to select gene-edited plants with expected indel patterns in the T0-generation.

A single-nucleotide substitution of CjTKPR1 determines pollen production in the gymnosperm plant Cryptomeria japonica.

Pollinosis, also known as pollen allergy or hay fever, is a global problem caused by pollen produced by various plant species. The wind-pollinated Japanese cedar (Cryptomeria japonica) is the largest contributor to severe pollinosis in Japan, where increasing proportions of people have been affected in recent decades. The MALE STERILITY 4 (MS4) locus of Japanese cedar controls pollen production, and its homozygous mutants (ms4/ms4) show abnormal pollen development after the tetrad stage and produce no mature pollen. In this study, we narrowed down the MS4 locus by fine mapping in Japanese cedar and found TETRAKETIDE α-PYRONE REDUCTASE 1 (TKPR1) gene in this region. Transformation experiments using Arabidopsis thaliana showed that single-nucleotide substitution ("T" to "C" at 244-nt position) of CjTKPR1 determines pollen production. Broad conservation of TKPR1 beyond plant division could lead to the creation of pollen-free plants not only for Japanese cedar but also for broader plant species.

Rapid survey of de novo mutations in naturally growing tree species following the March 2011 disaster in Fukushima: The effect of low-dose-rate radiation.

The impact of low-dose-rate radiation on genetics is largely unknown, particularly in natural environments. The Fukushima Dai-ich Nuclear Power Plant disaster resulted in the creation of contaminated natural lands. In this study, de novo mutations (DNMs) in germ line cells were surveyed from double-digest RADseq fragments in Japanese cedar and flowering cherry trees exposed to ambient dose rates ranging from 0.08 to 6.86 µGy h-1. These two species are among the most widely cultivated Japanese gymnosperm and angiosperm trees for forestry and horticultural purpose, respectively. For Japanese flowering cherry, open crossings were performed to produce seedlings, and only two candidate DNMs were detected from uncontaminated area. For Japanese cedar, the haploid megagametophytes were used as next generation samples. The use of megagametophytes from open crossing for next generation mutation screening had many advantages such as reducing exposure to radiation in contaminated areas because artificial crossings are not needed and the ease of data analysis owing to the haploid nature of megagametophytes. A direct comparison of the nucleotide sequences of parents and megagametophytes revealed an average of 1.4 candidate DNMs per megagametophyte sample (range: 0-40) after filtering procedures were optimized based on the validation of DNMs via Sanger sequencing. There was no relationship between the observed mutations and the ambient dose rate in the growing area or the concentration of 137Cs in cedar branches. The present results also suggest that mutation rates differ among lineages and that the growing environment has a relatively large influence on these mutation rates. These results suggested there was no significant increase in the mutation rate of the germplasm of Japanese cedar and flowering cherry trees growing in the contaminated areas.

A DNA barcode reference library for the native woody seed plants of Japan.

DNA barcode databases are increasingly available for a range of organisms, facilitating the wide application of DNA barcode-based studies. Here we announce the development of a comprehensive DNA barcode reference library of Japanese native woody seed plants representing 43 orders, 99 families, 303 genera and 834 species, and comprising 77.3% of the genera and 72.2% of the species of native woody seed plants in Japan. A total of 6216 plant specimens were collected from 223 sites across the subtropical, temperate, boreal and alpine biomes in Japan with most species represented by multiple accessions. This reference library utilized three chloroplast DNA regions (rbcL, trnH-psbA and matK) and consists of 14,403 barcode sequences. Individual regions varied in their identification rates, with species-level and genus-level rates for rbcL, trnH-psbA and matK based on blast being 57.4%/96.2%, 78.5%/99.1% and 67.8%/98.1%, respectively. Identification rates were higher using region combinations, with total species-level rates for two region combinations (rbcL & trnH-psbA, rbcL & matK and trnH-psbA & matK) ranging between 90.6% and 95.8%, and for all three regions being equal to 98.6%. Genus-level identification rates were even higher, ranging between 99.7% and 100% for two region combinations and being 100% for the three regions. These results indicate that this DNA barcode reference library is an effective resource for investigations of native woody seed plants in Japan using DNA barcodes and provides a useful template for the development of libraries for other components of the Japanese flora.

Genotype-by-environment interaction and genetic dissection of heartwood color in Cryptomeria japonica based on multiple common gardens and quantitative trait loci mapping.

The heartwood color of a major plantation tree Cryptomeria japonica shows high variability among clones and cultivars, and brighter heartwood has higher value in the usage of non-laminated wood such as in traditional construction, which makes heartwood color an important trait in breeding of this species. However, the genetic basis of the interactions between genetics and the environment on heartwood color has been understudied while these are necessary for effective breeding programs in multiple environmental condition. The objectives of the present study were to evaluate the effects of genetics and environments on heartwood color and how they interact in contrasting environments, and to identify genomic regions controlling heartwood color in C. japonica across multiple environments. Heartwood color in terms of L*a*b* color space and spectral reflectance was measured in common gardens established in three contrasting sites. Quantitative trait loci (QTL) that affect heartwood color were identified using previously constructed highly saturated linkage maps. Results found that heartwood color was largely genetically controlled, and genotype-by-environment interaction explained one-third of the total genetic variance of heartwood color. The effect of the environment was small compared to the effect of genetics, whereas environmental effects largely varied among heartwood color traits. QTL analysis identified a large number of QTLs with small to moderate effects (phenotypic variation explained of 6.6% on average). Some of these QTLs were stably expressed in multiple environments or had pleiotropic effects on heartwood color and moisture content. These results indicated that genetic variation in phenotypic plasticity plays an important role in regulating heartwood color and that the identified QTLs would maximize the breeding efficiency of heartwood color in C. japonica in heterogeneous environments.

Mutational effects of chronic gamma radiation throughout the life cycle of Arabidopsis thaliana: Insight into radiosensitivity in the reproductive stage.

Organisms living on Earth have always been exposed to natural sources of ionizing radiation, but following recent nuclear disasters, these background levels have often increased regionally due to the addition of man-made sources of radiation. To assess the mutational effects of ubiquitously present radiation on plants, we performed a whole-genome resequencing analysis of mutations induced by chronic irradiation throughout the life cycle of Arabidopsis thaliana grown under controlled conditions. We obtained resequencing data from 36 second generation post-mutagenesis (M2) progeny derived from 12 first generation (M1) lines grown under gamma-irradiation conditions, ranging from 0.0 to 2.0 Gray per day (Gy/day), to identify de novo mutations, including single base substitutions (SBSs) and small insertions/deletions (INDELs). The relationship between de novo mutation frequency and radiation dose rate from 0.0 to 2.0 Gy/day was assessed by statistical modeling. The increase in de novo mutations in response to irradiation dose fit the negative binomial model, which accounted for the high variability of mutation frequency observed. Among the different types of mutations, SBSs were more prevalent than INDELs, and deletions were more frequent than insertions. Furthermore, we observed that the mutational effects of chronic radiation were greater during the reproductive stage. These results will provide valuable insights into practical strategies for analyzing mutational effects in wild plants growing in environments with various mutagens.

Tree dynamic response and survival in a category-5 tropical cyclone: The case of super typhoon Trami.

In the future with climate change, we expect more forest and tree damage due to the increasing strength and changing trajectories of tropical cyclones (TCs). However, to date, we have limited information to estimate likely damage levels, and nobody has ever measured exactly how forest trees behave mechanically during a TC. In 2018, a category-5 TC destroyed trees in our ongoing research plots, in which we were measuring tree movement and wind speed in two different tree spacing plots. We found damaged trees in only the wider spaced plot. Here, we present how trees dynamically respond to strong winds during a TC. Sustained strong winds obviously trigger the damage to trees and forests but inter-tree spacing is also a key factor because the level of support from neighboring trees modifies the effective "stiffness" against the wind both at the single tree and whole forest stand level.

An Improved and Simplified Propagation System for Pollen-Free Sugi (Cryptomeria japonica) via Somatic Embryogenesis.

Sugi (Japanese cedar, Cryptomeria japonica) is the most important forestry tree species in Japan, covering 44% of the total artificial forest area. Large amounts of pollen released from these forests each spring cause allergic reactions in approximately 40% of the population, which are a serious social and public health problem in Japan. As a countermeasure, there is an urgent need to reforest using male-sterile plants (MSPs; pollen-free plants); however, the production of MSPs via conventional methods is inefficient, time consuming, and requires considerable resources in terms of labor and space. In the present paper, we described an improved and simplified methodology for the efficient propagation of pollen-free Japanese cedar, combining the use of genetic markers (marker-assisted selection or marker-aided selection) for the early selection of male-sterile genotypes and the use of somatic embryogenesis (SE) for the clonal mass propagation of seedlings. We describe all the stages involved in the production process of somatic seedlings. Our results demonstrated that this methodology easily and efficiently produces MSPs with a discrimination rate of 100% in a short period of time. Production of 243.6 ± 163.6 cotyledonary embryos per plate, somatic embryo germination, and plantlet conversion frequencies of 87.1 ± 11.9% and 84.8 ± 12.6%, respectively, and a 77.6 ± 12.1% survival rate after ex vitro acclimatization was achieved. Moreover, we also describe an easy method for the collection of somatic embryos prior to germination, as well as an efficient and practical method for their storage at 5°C. Finally, a representative schedule for the propagation of pollen-free sugi somatic seedlings is presented as a reference for practical uses. This methodology will definitively help to accelerate the production of C. japonica MSPs across Japan.

Marker-Assisted Selection for Pollen-Free Somatic Plants of Sugi (Japanese Cedar, Cryptomeria japonica): A Simple and Effective Methodology for Selecting Male-Sterile Mutants With ms1-1 and ms1-2.

Pollen allergy caused by sugi (Japanese cedar, Cryptomeria japonica) is a serious problem in Japan. One of the measures against pollinosis is the use of male-sterile plants (MSPs; pollen-free plants). In this context, the development of a novel technique for the efficient production of sugi MSPs, which combines marker-assisted selection (MAS) with somatic embryogenesis (SE), was recently reported by our research group. To improve the efficiency of MSP production, in this paper we report improved MAS for male-sterile individuals from embryogenic cells, cotyledonary embryos, and somatic plants of sugi using a newly developed marker in the form of the causative mutation of MS1 itself, selecting individuals with ms1-1 and ms1-2 male-sterile mutations. We also describe simplified methods for extracting DNA from different plant materials and for MAS using LAMP diagnostics. Finally, we show that MAS can be efficiently performed using the one-step indel genotyping (ING) marker developed in this study and using InstaGene for DNA extraction. The combination of SE and 100% accurate marker selection during the embryogenic cell stage enables the mass production of MS1 male-sterile sugi seedlings.

Factors Influencing Somatic Embryo Maturation in Sugi (Japanese Cedar, Cryptomeria japonica (Thunb. ex L.f.) D. Don).

This paper presents the results of several experiments identifying basal salts (BS) contained in maturation medium, polyethylene glycol (PEG) concentration, abscisic acid (ABA) concentration, additional supplementation with potassium chloride (KCl), amino acid (AA) concentration, and proliferation culture medium (PCM) as the main culture factors affecting somatic embryo maturation in sugi (Japanese cedar, Cryptomeria japonica, Cupressaceae). Highly efficient embryo maturation was achieved when embryogenic cell lines (ECLs) were cultured on media supplemented with a combination of PEG, ABA, and AAs. More than 1000 embryos per gram of fresh weight (FW) can be produced on EM maturation medium supplemented with 175 g L-1 PEG, 100 µM ABA, 2 g L-1 glutamine, 1 g L-1 asparagine, and 0.5 g L-1 arginine.

Somatic Embryogenesis Initiation in Sugi (Japanese Cedar, Cryptomeria japonica D. Don): Responses from Male-Fertile, Male-Sterile, and Polycross-Pollinated-Derived Seed Explants.

This study aimed to obtain information from several embryogenic cell (EC) genotypes analyzing the factors that affect somatic embryogenesis (SE) initiation in sugi (Cryptomeria japonica, Cupressaceae) to apply them in the improvement of protocols for efficient induction of embryogenic cell lines (ECLs). The results of several years of experiments including studies on the influence of initial explant, seed collection time, and explant genotype as the main factors affecting SE initiation from male-fertile, male-sterile, and polycross-pollinated-derived seeds are described. Initiation frequencies depending on the plant genotype varied from 1.35 to 57.06%. The best induction efficiency was achieved when seeds were collected on mid-July using the entire megagametophyte as initial explants. The extrusion of ECs started approximately after 2 weeks of culture, and the establishment of ECLs was observed mostly 4 weeks after extrusion on media with or without plant growth regulators (PGRs). Subsequently, induced ECLs were maintained and proliferated on media with PGRs by 2-3-week-interval subculture routines. Although, the initial explant, collection time, and culture condition played important roles in ECL induction, the genotype of the plant material of sugi was the most influential factor in SE initiation.

Construction of a reference transcriptome for the analysis of male sterility in sugi (Cryptomeria japonica D. Don) focusing on MALE STERILITY 1 (MS1).

Sugi (Cryptomeria japonica D. Don) is an important conifer used for afforestation in Japan. As the genome of this species is 11 Gbps, it is too large to assemble within a short timeframe. Transcriptomics is one approach that can address this deficiency. Here we designed a workflow consisting of three stages to de novo assemble transcriptome using Oases and Trinity. The three transcriptomic stage used were independent assembly, automatic and semi-manual integration, and refinement by filtering out potential contamination. We identified a set of 49,795 cDNA and an equal number of translated proteins. According to the benchmark set by BUSCO, 87.01% of cDNAs identified were complete genes, and 78.47% were complete and single-copy genes. Compared to other full-length cDNA resources collected by Sanger and PacBio sequencers, the extent of the coverage in our dataset was the highest, indicating that these data can be safely used for further studies. When two tissue-specific libraries were compared, there were significant expression differences between male strobili and leaf and bark sets. Moreover, subtle expression difference between male-fertile and sterile libraries were detected. Orthologous genes from other model plants and conifer species were identified. We demonstrated that our transcriptome assembly output (CJ3006NRE) can serve as a reference transcriptome for future functional genomics and evolutionary biology studies.

Identification and genetic diversity analysis of a male-sterile gene (MS1) in Japanese cedar (Cryptomeria japonica D. Don).

Identifying causative genes for a target trait in conifer reproduction is challenging for species lacking whole-genome sequences. In this study, we searched for the male-sterility gene (MS1) in Cryptomeria japonica, aiming to promote marker-assisted selection (MAS) of male-sterile C. japonica to reduce the pollinosis caused by pollen dispersal from artificial C. japonica forests in Japan. We searched for mRNA sequences expressed in male strobili and found the gene CJt020762, coding for a lipid transfer protein containing a 4-bp deletion specific to male-sterile individuals. We also found a 30-bp deletion by sequencing the entire gene of another individual with the ms1. All nine breeding materials with the allele ms1 had either a 4-bp or 30-bp deletion in gene CJt020762, both of which are expected to result in faulty gene transcription and function. Furthermore, the 30-bp deletion was detected from three of five individuals in the Ishinomaki natural forest. From our findings, CJt020762 was considered to be the causative gene of MS1. Thus, by performing MAS using two deletion mutations as a DNA marker, it will be possible to find novel breeding materials of C. japonica with the allele ms1 adapted to the unique environment of each region of the Japanese archipelago.

Genotype and transcriptome effects on somatic embryogenesis in Cryptomeria japonica.

Somatic embryogenesis (SE), which is in vitro regeneration of plant bodies from somatic cells, represents a useful means of clonal propagation and genetic engineering of forest trees. While protocols to obtain calluses and induce regeneration in somatic embryos have been reported for many tree species, the knowledge of molecular mechanisms of SE development is still insufficient to achieve an efficient supply of somatic embryos required for the industrial application. Cryptomeria japonica, a conifer species widely used for plantation forestry in Japan, is one of the tree species waiting for a secure SE protocol; the probability of normal embryo development appears to depend on genotype. To discriminate the embryogenic potential of embryonal masses (EMs) and efficiently obtain normal somatic embryos of C. japonica, we investigated the effects of genotype and transcriptome on the variation in embryogenic potential. Using an induction experiment with 12 EMs each from six genotypes, we showed that embryogenic potential differs between/within genotypes. Comparisons of gene expression profiles among EMs with different embryogenic potentials revealed that 742 differently expressed genes were mainly associated with pattern forming and metabolism. Thus, we suggest that not only genotype but also gene expression profiles can determine success in SE development. Consistent with previous findings for other conifer species, genes encoding leafy cotyledon, wuschel, germin-like proteins, and glutathione-S-transferases are likely to be involved in SE development in C. japonica and indeed highly expressed in EMs with high-embryogenic potential; therefore, these proteins represent candidate markers for distinguishing embryogenic potential.

Contrasted patterns of local adaptation to climate change across the range of an evergreen oak, Quercus aquifolioides.

Long-lived tree species are genetically differentiated and locally adapted with respect to fitness-related traits, but the genetic basis of local adaptation remains largely unresolved. Recent advances in population genetics and landscape genomic analyses enable identification of putative adaptive loci and specific selective pressures acting on local adaptation. Here, we sampled 60 evergreen oak (Quercus aquifolioides) populations throughout the species' range and pool-sequenced 587 individuals at drought-stress candidate genes. We analyzed patterns of genetic diversity and differentiation for 381 single nucleotide polymorphisms (SNPs) from 65 candidate genes and eight microsatellites. Outlier loci were identified by genetic differentiation analysis and genome-environment associations. The response pattern of genetic variation to environmental gradient was assessed by linear isolation-by-distance/environment tests, redundancy analysis, and nonlinear methods. SNPs and microsatellites revealed two genetic lineages: Tibet and Hengduan Mountains-Western Sichuan Plateau (HDM-WSP), with reduced genetic diversity in Tibet lineage. More outlier loci were detected in HDM-WSP lineage than Tibet lineage. Among these, three SNPs in two genes responded to dry season precipitation in the HDM-WSP lineage but not in Tibet. By contrast, genetic variation in the Tibet lineage was related to geographic distance instead of the environment. Furthermore, risk of nonadaptedness (RONA) analyses suggested HDM-WSP lineage will have a better capacity to adapt in the predicted future climate compared with the Tibet lineage. We detected genetic imprints consistent with natural selection and molecular adaptation to drought on the Qinghai-Tibet Plateau (QTP) over a range of long-lived and widely distributed oak species in a changing environment. Our results suggest that different within-species adaptation processes occur in species occurring in heterogeneous environments.

Development of diagnostic PCR and LAMP markers for MALE STERILITY 1 (MS1) in Cryptomeria japonica D. Don.

OBJECTIVE: Due to the allergic nature of the pollen of Cryptomeria japonica, the most important Japanese forestry conifer, a pollen-free cultivar is preferred. Mutant trees detected in nature have been used to produce a pollen-free cultivar. In order to reduce the time and cost needed for production and breeding, we aimed to develop simple diagnostic molecular markers for mutant alleles of the causative gene MALE STERILITY 1 (MS1) in C. japonica to rapidly identify pollen-free mutants. RESULTS: We developed PCR and LAMP markers to detect mutant alleles and to present experimental options depending on available laboratory equipment. LAMP markers were developed for field stations, where PCR machines are unavailable. The LAMP method only needs heat-blocks or a water bath to perform the isothermal amplification and assay results can be read by the naked eye. Because the causative mutations were deletions, we developed two kinds of PCR markers, amplified length polymorphism (ALP) and allele specific PCR (ASP) markers. These assays can be visualized using capillary or agarose gel electrophoresis.

Somatic Embryogenesis and Plant Regeneration from Sugi (Japanese Cedar, Cryptomeria japonica D. Don, Cupressaceae) Seed Families by Marker Assisted Selection for the Male Sterility Allele ms1.

One of the possible countermeasures for pollinosis caused by sugi (Cryptomeria japonica), a serious public health problem in Japan, is the use of male sterile plants (MSPs; pollen-free plants). However, the production efficiencies of MSPs raised by conventional methods are extremely poor, time consuming, and resulting in a high seedling cost. Here, we report the development of a novel technique for efficient production of MSPs, which combines marker-assisted selection (MAS) and somatic embryogenesis (SE). SE from four full sib seed families of sugi, carrying the male sterility gene MS1, was initiated using megagametophyte explants that originated from four seed collections taken at one-week intervals during the month of July 2017. Embryogenic cell lines (ECLs) were achieved in all families, with initiation rates varying from 0.6% to 59%. Somatic embryos were produced from genetic marker-selected male sterile ECLs on medium containing maltose, abscisic acid (ABA), polyethylene glycol (PEG), and activated charcoal (AC). Subsequently, high frequencies of germination and plant conversion (≥76%) were obtained on plant growth regulator-free medium. Regenerated plantlets were acclimatized successfully, and the initial growth of male sterile somatic plants was monitored in the field.

Spatial variation in bird pollination and its mitigating effects on the genetic diversity of pollen pools accepted by Camellia japonica trees within a population at a landscape level.

Bird pollination can vary spatially in response to spatial fluctuations in flowering even within plant populations. In this study, we examined the hypothesis that the spatial variation in bird pollination may induce mitigating effects, which maintains or increases genetic diversity of pollen pools at local sites with low flowering densities. To test this hypothesis, we analyzed the landscape-level genetic effects within a population of Camellia japonica on the pollen pools accepted by individuals in two reproductive years by using genotypes at eight microsatellite loci of 1323 seeds from 19 seed parents. Regression analyses using the quadratic models of correlated paternity between pollen pools against spatial distances between the seed-parent pairs revealed not only local pollination but also some amount of long-distance pollen dispersal. The genetic diversity of pollen pools accepted by seed parents tended to be negatively related to the densities of flowering individuals near the seed parents during winter (when the effective pollination of C. japonica is mediated mostly by Zosterops japonica). We show that the low density of flowering individuals may induce the expansion of the foraging areas of Z. japonica and consequently increase the genetic diversity of pollen pools. This spatial variation in bird pollination may induce the mitigating effects on the C. japonica population. The comparisons between the two study years indicate that the overall pattern of bird pollination and the genetic effects described here, including the mitigating effects, may be stable over time.

Scanning RNA-Seq and RAD-Seq approach to develop SNP markers closely linked to MALE STERILITY 1 (MS1) in Cryptomeria japonica D. Don.

Cryptomeria japonica is a major forestry tree species in Japan. Male sterility of the species is caused by a recessive gene, which shows dysfunction of pollen development and results in no dispersed pollen. Because the pollen of C. japonica induces pollinosis, breeding of pollen-free C. japonica is desired. In this study, single nucleotide polymorphism (SNP) markers located at 1.78 and 0.58 cM to a male sterility locus (MS1) were identified from an analysis of RNA-Seq and RAD-Seq, respectively. SNPs closely linked to MS1 were first scanned by a method similar to MutMap, where a type of index was calculated to measure the strength of the linkage between a marker sequence and MS1. Linkage analysis of selected SNP markers confirmed a higher efficiency of the current method to construct a partial map around MS1. Allele-specific PCR primer pair for the most closely linked SNP with MS1 was developed as a codominant marker, and visualization of the PCR products on an agarose gel enabled rapid screening of male sterile C. japonica. The allele-specific primers developed in this study would be useful for establishing the selection of male sterile C. japonica.