Yoshinari Moriguchi

Development and characterization of microsatellite markers for Cryptomeria japonica D.Don.

Abstract.Thirty four microsatellite markers for Cryptomeria japonica D. Don were developed by searching three types of library: a database of C. japonica cDNA sequences, a standard non-enriched genomic DNA library and a microsatellite-enriched library using magnetic particles. The enrichment of microsatellite sequences using magnetic particles is very efficient compared to the other two methods both in terms of the numbers of markers generated, and in the polymorphism they detect. The microsatellites developed from the genomic DNA library generally have longer repeat sequences and are more polymorphic than those from cDNA. All the developed microsatellite markers in this study showed polymorphism among 28 plus trees selected from locations scattered throughout Japan. The mean number of alleles per locus (MNA) detected in the 28 plus trees ranged from 2 to 21 with an average of 7.5. The Polymorphism Information Content (PIC) ranged from 0.160 to 0.936 with an average of 0.666. Co-dominant segregation of alleles in a three-generation pedigree of C. japonica was demonstrated at 34 microsatellite loci, and the segregation was not distorted from Mendelian expectation for all loci. In 12 out of 34 loci, a null allele was detected. Key relationships between polymorphic parameters, such as MNA and PIC, and the characteristics of microsatellite sequences, such as the longest repeat number, total repeat number and total number of nucleotides, were investigated using rank correlation coefficients, Kendalls τ. A positive correlation was found between repeat lengths and polymorphisms. The markers provide sufficient resolution for investigating gene flow within forests and seed orchards, and for genome mapping.

A second generation framework for the analysis of microsatellites in expressed sequence tags and the development of EST-SSR markers for a conifer, Cryptomeria japonica

BackgroundMicrosatellites or simple sequence repeats (SSRs) in expressed sequence tags (ESTs) are useful resources for genome analysis because of their abundance, functionality and polymorphism. The advent of commercial second generation sequencing machines has lead to new strategies for developing EST-SSR markers, necessitating the development of bioinformatic framework that can keep pace with the increasing quality and quantity of sequence data produced. We describe an open scheme for analyzing ESTs and developing EST-SSR markers from reads collected by Sanger sequencing and pyrosequencing of sugi (Cryptomeria japonica).ResultsWe collected 141,097 sequence reads by Sanger sequencing and 1,333,444 by pyrosequencing. After trimming contaminant and low quality sequences, 118,319 Sanger and 1,201,150 pyrosequencing reads were passed to the MIRA assembler, generating 81,284 contigs that were analysed for SSRs. 4,059 SSRs were found in 3,694 (4.54%) contigs, giving an SSR frequency lower than that in seven other plant species with gene indices (5.4–21.9%). The average GC content of the SSR-containing contigs was 41.55%, compared to 40.23% for all contigs. Tri-SSRs were the most common SSRs; the most common motif was AT, which was found in 655 (46.3%) di-SSRs, followed by the AAG motif, found in 342 (25.9%) tri-SSRs. Most (72.8%) tri-SSRs were in coding regions, but 55.6% of the di-SSRs were in non-coding regions; the AT motif was most abundant in 3′ untranslated regions. Gene ontology (GO) annotations showed that six GO terms were significantly overrepresented within SSR-containing contigs. Forty–four EST-SSR markers were developed from 192 primer pairs using two pipelines: read2Marker and the newly-developed CMiB, which combines several open tools. Markers resulting from both pipelines showed no differences in PCR success rate and polymorphisms, but PCR success and polymorphism were significantly affected by the expected PCR product size and number of SSR repeats, respectively. EST-SSR markers exhibited less polymorphism than genomic SSRs.ConclusionsWe have created a new open pipeline for developing EST-SSR markers and applied it in a comprehensive analysis of EST-SSRs and EST-SSR markers in C. japonica. The results will be useful in genomic analyses of conifers and other non-model species.

Gene flow and mating system in five Cryptomeria japonica D. Don seed orchards as revealed by analysis of microsatellite markers

We investigated gene flow in five Cryptomeria japonica D. Don seed orchards of two different types (common and miniature) at widely spaced locations using microsatellite markers. The quality of a seed crop is determined by many factors, including pollen contamination from outside sources, self-fertilization, and the proportion of contributions from constituent clones. Contamination rates were found to vary among ramets both within seed orchards (10.0–76.7% in the most variable seed orchard) and among seed orchards (35.0–65.8% on average). Among ramets, there were significant negative correlations between pollen contamination rate and their distance from the orchard edge; among seed orchards, there were significant positive correlations between the pollen contamination rate and the C. japonica forest area nearby. Some proportion of the pollen (10.7% of total contamination) also migrated from parts of the orchards that had not been treated with gibberellin to induce flowering. Self-fertilization rates varied among seed orchards (1.4–4.4% on average), and there were significant positive correlations between self-fertilization rate and the number of ramets per clone in both types of seed orchard. Contributions as pollen donors differed significantly among clones in all seed orchards. The distance between planted ramets, flowering phenology, and relative pollen fecundity may also have contributed to observed differences in paternal contribution. The influence of these factors on genetic potential did not differ greatly between the two types of orchards.

The construction of a high-density linkage map for identifying SNP markers that are tightly linked to a nuclear-recessive major gene for male sterility in Cryptomeria japonica D. Don

BackgroundHigh-density linkage maps facilitate the mapping of target genes and the construction of partial linkage maps around target loci to develop markers for marker-assisted selection (MAS). MAS is quite challenging in conifers because of their large, complex, and poorly-characterized genomes. Our goal was to construct a high-density linkage map to facilitate the identification of markers that are tightly linked to a major recessive male-sterile gene (ms1) for MAS in C. japonica, a species that is important in Japanese afforestation but which causes serious social pollinosis problems.ResultsWe constructed a high-density saturated genetic linkage map for C. japonica using expressed sequence-derived co-dominant single nucleotide polymorphism (SNP) markers, most of which were genotyped using the GoldenGate genotyping assay. A total of 1261 markers were assigned to 11 linkage groups with an observed map length of 1405.2 cM and a mean distance between two adjacent markers of 1.1 cM; the number of linkage groups matched the basic chromosome number in C. japonica. Using this map, we located ms1 on the 9th linkage group and constructed a partial linkage map around the ms1 locus. This enabled us to identify a marker (hrmSNP970_sf) that is closely linked to the ms1 gene, being separated from it by only 0.5 cM.ConclusionsUsing the high-density map, we located the ms1 gene on the 9th linkage group and constructed a partial linkage map around the ms1 locus. The map distance between the ms1 gene and the tightly linked marker was only 0.5 cM. The identification of markers that are tightly linked to the ms1 gene will facilitate the early selection of male-sterile trees, which should expedite C. japonica breeding programs aimed at alleviating pollinosis problems without harming productivity.

Evidence for cryptic northern refugia in the last glacial period in Cryptomeria japonica

BACKGROUND AND AIMS Distribution shifts and natural selection during past climatic changes are important factors in determining the genetic structure of forest species. In particular, climatic fluctuations during the Quaternary appear to have caused changes in the distribution ranges of plants, and thus strongly affected their genetic structure. This study was undertaken to identify the responses of the conifer Cryptomeria japonica, endemic to the Japanese Archipelago, to past climatic changes using a combination of phylogeography and species distribution modelling (SDM) methods. Specifically, this study focused on the locations of refugia during the last glacial maximum (LGM). METHODS Genetic diversity and structure were examined using 20 microsatellite markers in 37 populations of C. japonica. The locations of glacial refugia were assessed using STRUCTURE analysis, and potential habitats under current and past climate conditions were predicted using SDM. The process of genetic divergence was also examined using the approximate Bayesian computation procedure (ABC) in DIY ABC to test the divergence time between the gene pools detected by the STRUCTURE analysis. KEY RESULTS STRUCTURE analysis identified four gene pools: northern Tohoku district; from Chubu to Chugoku district; from Tohoku to Shikoku district on the Pacific Ocean side of the Archipelago; and Yakushima Island. DIY ABC analysis indicated that the four gene pools diverged at the same time before the LGM. SDM also indicated potential northern cryptic refugia. CONCLUSIONS The combined evidence from microsatellites and SDM clearly indicates that climatic changes have shaped the genetic structure of C. japonica. The gene pool detected in northern Tohoku district is likely to have been established by cryptic northern refugia on the coast of the Japan Sea to the west of the Archipelago. The gene pool in Yakushima Island can probably be explained simply by long-term isolation from the other gene pools since the LGM. These results are supported by those of SDM and the predicted divergence time determined using ABC analysis.

Demonstration of genome-wide association studies for identifying markers for wood property and male strobili traits in Cryptomeria japonica.

Genome-wide association studies (GWAS) are an alternative to bi-parental QTL mapping in long-lived perennials. In the present study, we examined the potential of GWAS in conifers using 367 unrelated plus trees of Cryptomeria japonica D. Don, which is the most widely planted and commercially important tree species in Japan, and tried to detect significant associations between wood property traits and quantity of male strobili on the one hand, and 1,032 single nucleotide polymorphisms (SNPs) assigned to 1,032 genes on the other. Association analysis was performed with the mixed linear model taking into account kinship relationships and subpopulation structure. In total, 6 SNPs were found to have significant associations with the variations in phenotype. These SNPs were not associated with the positions of known genes and QTLs that have been reported to date, thus they may identify novel QTLs. These 6 SNPs were all found in sequences showing similarities with known genes, although further analysis is required to dissect the ways in which they affect wood property traits and abundance of male strobili. These presumptive QTL loci provide opportunities for improvement of C. japonica, based on a marker approach. The results suggest that GWAS has potential for use in future breeding programs in C. japonica.

Factors Influencing Male Reproductive Success in a Cryptomeria japonica Seed Orchard Revealed by Microsatellite Marker Analysis

Abstract We investigated the influence of male flower production, floral synchrony and inter-tree distances on male reproductive success in a miniature seed orchard of Cryptomeria japonica. We used six microsatellite markers to determine the paternity of each seed. In the seed orchard, the average pollen contamination and clonal self-fertilization rates were 38.7% and 1.7%, respectively. The level of male reproductive success of constituent clones varied from 0.0 to 15.7%. Five clones showing the highest male reproductive success contributed ca. 30% of all analyzed seeds as a pollen donor after excluding contamination by external sources of pollen. The statistical analyses showed that male reproductive success was strongly influenced by male flower production of each clone and, possibly, by their distance to the mother trees. The linear regression which included male flower production and floral synchrony as independent variables, however, accounted for only 14.7% of variation of male reproductive success, suggesting that other factors such as pollen competition might also influence male reproductive success. Since we found no significant correlation between male reproductive and female reproductive successes, it may be better to equalize male and female reproductive successes independently

Genetic Differentiation and Evolutionary Adaptation in Cryptomeria japonica

Local adaptation of plant species is a central issue for survival during global climate change, especially for long-lived forest trees, with their lengthy regeneration time and spatially limited gene flow. Identification of loci and/or genomic regions associated with local adaptation is necessary for knowledge of both evolution and molecular breeding for climate change. Cryptomeria japonica is an important species for forestry in Japan; it has a broad natural distribution and can survive in a range of different environments. The genetic structure of 14 natural populations of this species was investigated using 3930 SNP markers. Populations on the Pacific Ocean side of Japan are clearly different from those on the Japan Sea side, as discussed in previous studies. Structure analysis and population network trees show that peripheral populations, including the most northerly and southerly ones, have unique features. We found that the genetic differentiation coefficient is low, FST = 0.05, although it must account for the presence of important genes associated with adaptation to specific environments. In total, 208 outlier loci were detected, of which 43 were associated with environmental variables. Four clumped regions of outlier loci were detected in the genome by linkage analysis. Linkage disequilibrium (LD) was quite high in these clumps of outlier loci, which were found in linkage groups (LGs) 2, 7, 10, and 11, especially between populations of two varieties, and when interchromosomal LD was also detected. The LG7 region is characteristic of the Yakushima population, which is a large, isolated, peripheral population occupying a specific environment resulting from isolation combined with volcanic activity in the region. The detected LD may provide strong evidence for selection between varieties.

Characterization of EST–SSRs from Cryptomeria japonica

Simple sequence repeat (SSR) markers were developed for Cryptomeria japonica from 55,530 expressed sequence tags (ESTs). Of the 219 designed primer pairs, 176 showed clear PCR amplification products and 27 of these revealed polymorphisms in six individuals sampled at sites across Japan. When these 27 loci were examined in 48 individuals from a single population, the number of alleles and the Polymorphism Information Content (PIC) ranged from 2 to 17 and 0.09 to 0.86, respectively. The average PIC value obtained from EST–SSR markers (0.40) was lower than that from genomic SSR markers (0.67 and 0.62). We observed a significant departure from Hardy–Weinberg equilibrium in five of the 27 loci. These EST–SSR markers will be useful for the evaluation of genetic variation in C. japonica and for genetic mapping.

Prominent Differences in Leaf Fatty Acid Composition in the F1 Hybrid Compared with Parent Trees Larix gmelinii var. japonica and L. kaempferi

Fatty acid (FA) compositions in leaves were investigated for two families of F1 hybrids of Larix gmelinii var. japonica × L. kaempferi (F1) and their parent clones. Twenty-one FAs, from C12 to C32, were found in the leaves of both adult trees and seedlings. The levels of 18:1/(18:2 + 18:3) increased in the order L. kaempferi, F1, and L. gmelinii var. japonica, with significant differences between L. gmelinii var. japonica and F1 in adult trees, but these differences were not found in the seedlings. Moreover, in the adult trees, the 18:1/(18:2 + 18:3) levels in the neutral phospholipid fraction and the ΣC18/ΣC16, especially in the glycolipid fraction, showed significant differences among the three species. These characteristics are discussed from the viewpoint of lipid synthesis in the endoplasmic reticulum and chloroplasts, and of the activities and substrate specificities in sequential FA desaturation. Linear discriminant analysis suggested that the FA compositions are useable as an index in the identification of hybrid seedlings.