A chromosome-level genome assembly of a model conifer plant, the Japanese cedar, Cryptomeria japonica D. Don.

BACKGROUND: The Japanese cedar (Cryptomeria japonica D. Don) is one of the most important Japanese forest trees, occupying approximately 44% of artificial forests and planted in East Asia, the Azores Archipelago, and certain islands in the Indian Ocean. Although the huge genome of the species (ca. 9 Gbp) with abundant repeat elements may have represented an obstacle for genetic analysis, this species is easily propagated by cutting, flowered by gibberellic acid, transformed by Agrobacterium, and edited by CRISPR/Cas9. These characteristics of C. japonica recommend it as a model conifer species for which reference genome sequences are necessary. RESULTS: Herein, we report the first chromosome-level assembly of C. japonica (2n = 22) using third-generation selfed progeny (estimated homozygosity rate = 0.96). Young leaf tissue was used to extract high molecular weight DNA (> 50 kb) for HiFi PacBio long-read sequencing and to construct an Hi-C/Omni-C library for Illumina short-read sequencing. The 29× and 26× genome coverage of HiFi and Illumina reads, respectively, for de novo assembly yielded 2,651 contigs (9.1 Gbp, N50 contig size 12.0 Mbp). Hi-C analysis mapped 97% of the nucleotides on 11 chromosomes. The assembly was verified through comparison with a consensus linkage map comprising 7,781 markers. BUSCO analysis identified ∼ 91% conserved genes. CONCLUSIONS: Annotations of genes and comparisons of repeat elements with other Cupressaceae and Pinaceae species provide a fundamental resource for conifer research.

Secondary contact zone and genetic introgression in closely related haplodiploid social spider mites.

How frequently hybridisation and gene flow occur in the contact zones of diverging taxa is important for understanding the speciation process. Stigmaeopsis sabelisi and Stigmaeopsis miscanthi high-aggression form (hereafter, S. miscanthi HG) are haplodiploid, social spider mites that infest the Chinese silver grass, Miscanthus sinensis. These two species are closely related and parapatrically distributed in Japan. In mountainous areas, S. sabelisi and S. miscanthi HG are often found in the highlands and lowlands, respectively, suggesting that they are in contact at intermediate altitudes. It is estimated that they diverged from their common ancestors distributed in subtropical regions (south of Japan) during the last glacial period, expanded their distribution into the Japanese Archipelago, and came to have such a parapatric distribution (secondary contact). As their reproductive isolation is strong but incomplete, hybridisation and genetic introgression are expected at their distributional boundaries. In this study, we investigated their spatial distribution patterns along the elevation on Mt. Amagi using male morphological differences, and investigated their hybridisation status using single-nucleotide polymorphisms by MIG-seq. We found their contact zone at altitudes of 150-430 m, suggesting that their contact zone is prevalent in the parapatric area, which is in line with a previous study. Interspecific mating was predicted based on the sex ratio in the contact zone. No obvious hybrids were found, but genetic introgression was detected although it was extremely low.

Divergent mechanisms of reduced growth performance in Betula ermanii saplings from high-altitude and low-latitude range edges.

The reduced growth performance of individuals from range edges is a common phenomenon in various taxa, and considered to be an evolutionary factor that limits the species' range. However, most studies did not distinguish between two mechanisms that can lead to this reduction: genetic load and adaptive selection to harsh conditions. To address this lack of understanding, we investigated the climatic and genetic factors underlying the growth performance of Betula ermanii saplings transplanted from 11 populations including high-altitude edge and low-latitude edge population. We estimated the climatic position of the populations within the overall B. ermanii's distribution, and the genetic composition and diversity using restriction-site associated DNA sequencing, and measured survival, growth rates and individual size of the saplings. The high-altitude edge population (APW) was located below the 95% significance interval for the mean annual temperature range, but did not show any distinctive genetic characteristics. In contrast, the low-latitude edge population (SHK) exhibited a high level of linkage disequilibrium, low genetic diversity, a distinct genetic composition from the other populations, and a high relatedness coefficient. Both APW and SHK saplings displayed lower survival rates, heights and diameters, while SHK saplings also exhibited lower growth rates than the other populations' saplings. The low heights and diameters of APW saplings was likely the result of adaptive selection to harsh conditions, while the low survival and growth rates of SHK saplings was likely the result of genetic load. Our findings shed light on the mechanisms underlying the reduced growth performance of range-edge populations.

Inferring the demographic history of Japanese cedar, Cryptomeria japonica, using amplicon sequencing.

The evolution of a species depends on multiple forces, such as demography and natural selection. To understand the trajectory and driving forces of evolution of a target species, it is first necessary to uncover that species' population history, such as past and present population sizes, subdivision and gene flow, by using appropriate genetic markers. Cryptomeria japonica is a long-lived monoecious conifer species that is distributed in Japan. There are two main lines (omote-sugi and ura-sugi), which are distinguished by apparent differences in morphological traits that may have contributed to their local adaptation. The evolution of these morphological traits seems to be related to past climatic changes in East Asia, but no precise estimate is available for the divergence time of these two lines and the subsequent population dynamics in this species. Here, we analyzed the nucleotide variations at 120 nuclear genes in 94 individuals by using amplicon sequencing in combination with high-throughput sequencing technologies. Our analysis indicated that the population on Yakushima Island, the southern distribution limit of C. japonica in Japan, diverged from the other populations 0.85 million years ago (MYA). The divergence time of the other populations on mainland Japan was estimated to be 0.32 MYA suggesting that the divergence of omote-sugi and ura-sugi might have occurred before the last glacial maximum. Although we found modest levels of gene flow between the present populations, the long-term isolation and environmental heterogeneity caused by climatic changes might have contributed to the differentiation of the lines and their local adaptation.

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.

TodoFirGene: Developing Transcriptome Resources for Genetic Analysis of Abies sachalinensis.

Todo-matsu (Abies sachalinensis) is one of the most important forestry species in Hokkaido, Japan and is distributed from near sea level to the alpine zone. Due to its wide spatial distribution, the species adapts to its environment, displaying phenotypes of ecological relevance. In order to identify candidate genes under natural selection, we collected the transcriptome from the female and male flower, leaf and inner bark. De novo assembly with 34.7 Gb of sequencing reads produced 158,542 transcripts from 69,618 loci, whose estimated coverage reached 95.6% of conserved eukaryotic genes. Homology searches against publicly available databases identified 134,190 (84.6%) transcripts with at least one hit. In total, 28,944 simple sequence repeats (SSRs) and 80,758 single nucleotide variants (SNVs) were detected from 23,570 (14.9%) and 25,366 (16.0%) transcripts, which were valuable for use in genetic analysis of the species. All the annotations were included in a relational database, TodoFirGene, which provides an interface for various queries and homology search, and can be accessed at http://plantomics.mind.meiji.ac.jp/todomatsu/. This database hosts not only the A. sachalinensis transcriptome but also links to the proteomes of 13 other species, allowing a comparative genomic study of plant species.

Large contribution of clonal reproduction to the distribution of deciduous liana species (Wisteria floribunda) in an old-growth cool temperate forest: evidence from genetic analysis.

Background and Aims: Extensive clonal (vegetative) reproduction in lianas is a common and important life history strategy for regeneration and colonization success. However, few studies have evaluated the contribution of clonal reproduction to stand-level distribution of lianas in their natural habitat using genetic tools. The objectives of the present study were to investigate (1) the contribution of clonal reproduction to the distribution of Wisteria floribunda, (2) the size of clonal patches and (3) how the distribution patterns of W. floribunda clones are affected by micro-topography. Methods: The contribution of clonal reproduction to the distribution of the deciduous liana species W. floribunda was evaluated using genetic analysis across a 6-ha plot of an old-growth temperate forest in Japan and preference in landform between clonal ramets and non-clonal ramets was assessed. Key Results: Of the 391 ramets sampled, clonal reproduction contributed to 71 and 62 % of the total abundance and basal area, respectively, or 57 and 31 % when the largest ramet within a genet was excluded. The large contribution of clonal reproduction to the density and basal area of W. floribunda was consistent with previous observational studies. The largest genet included a patch size of 0.47 ha and ranged over 180 m. Preferred landforms of clonal and non-clonal ramets were significantly different when evaluated by both abundance and basal area. Non-clonal ramets distributed more on lower part of the slope than other landforms in comparison with clonal ramets and trees, possibly reflecting the limitation of clonal growth by stolons. Conclusions: Using genetic analysis, the present study found evidence of a large contribution of clonal reproduction on the distribution of W. floribunda in its natural habitat. The results indicate that clonal reproduction plays an important role not only in the formation of populations but also in determining the distribution patterns of liana species.

Inferences of population structure and demographic history for Taxodium distichum, a coniferous tree in North America, based on amplicon sequencing analysis.

PREMISE OF THE STUDY: Studies of natural genetic variation can elucidate the genetic basis of phenotypic variation and the past population structure of species. Our study species, Taxodium distichum, is a unique conifer that inhabits the flood plains and swamps of North America. Morphological and ecological differences in two varieties, T. distichum var. distichum (bald cypress) and T. distichum var. imbricarium (pond cypress), are well known, but little is known about the level of genetic differentiation between the varieties and the demographic history of local populations. METHODS: We analyzed nucleotide polymorphisms at 47 nuclear loci from 96 individuals collected from the Mississippi River Alluvial Valley (MRAV), and Gulf Coastal populations in Texas, Louisiana, and Florida using high-throughput DNA sequencing. Standard population genetic statistics were calculated, and demographic parameters were estimated using a composite-likelihood approach. KEY RESULTS: Taxodium distichum in North America can be divided into at least three genetic groups, bald cypress in the MRAV and Texas, bald cypress in Florida, and pond cypress in Florida. The levels of genetic differentiation among the groups were low but significant. Several loci showed the signatures of positive selection, which might be responsible for local adaptation or varietal differentiation. CONCLUSIONS: Bald cypress was genetically differentiated into two geographical groups, and the boundary was located between the MRAV and Florida. This differentiation could be explained by population expansion from east to west. Despite the overlap of the two varieties' ranges, they were genetically differentiated in Florida. The estimated demographic parameters suggested that pond cypress split from bald cypress during the late Miocene.

Complex pollination of a tropical Asian rainforest canopy tree by flower-feeding thrips and thrips-feeding predators.

PREMISE OF THE STUDY: In tropical rainforests of Southeast Asia, a highly fecund thrips (Thrips spp.) responds rapidly to the mass flowering at multiple-year intervals characteristic of certain species such as the canopy tree studied here, Shorea acuminata, by feeding on flower resources. However, past DNA analyses of pollen adherent to thrips bodies revealed that the thrips promoted a very high level of self-pollination. Here, we identified the pollinator that contributes to cross-pollination and discuss ways that the pollination system has adapted to mass flowering. METHODS: By comparing the patterns of floral visitation and levels of genetic diversity in adherent pollen loads among floral visitors, we evaluated the contribution of each flower visitor to pollination. KEY RESULTS: The big-eyed bug, Geocoris sp., a major thrips predator, was an inadvertent pollinator, and importantly contributed to cross-pollination. The total outcross pollen adhering to thrips was approximately 30% that on the big-eyed bugs. Similarly, 63% of alleles examined in S. acuminata seeds and seedlings occurred in pollen adhering to big-eyed bugs; about 30% was shared with pollen from thrips. CONCLUSIONS: During mass flowering, big-eyed bugs likely travel among flowering S. acuminata trees, attracted by the abundant thrips. Floral visitation patterns of big-eyed bugs vs. other insects suggest that these bugs can maintain their population size between flowering by preying upon another thrips (Haplothrips sp.) that inhabits stipules of S. acuminata throughout the year and quickly respond to mass flowering. Thus, thrips and big-eyed bugs are essential components in the pollination of S. acuminata.

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.

Cryptic diploid lineage of Betula ermanii at its southern boundary populations in Japan.

Polyploidy is thought to enable species diversification and adaptation to extreme environments. Resolving the ecological differences between a taxon's ploidy levels would therefore provide important insights into local adaptation and speciation. The genus Betula includes many polyploids, but estimates of their phylogenetic relationships and evolutionary history are uncertain because of cryptic lineages and species. As one of the southern boundary populations of Betula ermanii in Japan has been shown to have distinctive genetic characteristics and traits, the differences in ploidy levels between three southern boundary and various other Japanese B. ermanii populations were investigated using flow cytometry. Leaf and seed morphologies were also compared. Apart from individuals in southern boundary populations, all those sampled were tetraploid. Individuals from the southern boundary populations were mostly diploid, apart from a few from lower altitude Shikoku populations, which were tetraploid. Leaf and seed morphologies differed between tetraploids and diploids. Diploid individuals were characterized by leaves with a heart-shaped base and many leaf teeth, and seeds with relatively longer wings. The diploid populations could be considered a cryptic relict lineage of B. ermanii, and there is a possibility that this lineage is a diploid ancestor of B. ermanii and a relict population of the Sohayaki element. Further investigation of the Japanese Betula phylogenetic relationships would enable an informed discussion of taxonomic revisions.

Response of stomatal density and size in Betula ermanii to contrasting climate conditions: The contributions of genetic and environmental factors.

As plant distribution and performance are determined by both environmental and genetic factors, clarifying the contribution of these two factors is a key for understanding plant adaptation and predicting their distribution under ongoing global warming. Betula ermanii is an ideal species for such research because of its wide distribution across diverse environments. Stomatal density and size are crucial traits that plants undergo changes in to adapt to different environments as these traits directly influence plant photosynthesis and transpiration. In this study, we conducted a multi-location common garden experiment using B. ermanii to (1) clarify the contribution of both environmental and genetic factors to the variation in stomatal density and size of B. ermanii, (2) demonstrate the differences in the plasticity of stomatal density and size among B. ermanii populations, and (3) understand how stomatal density and size of B. ermanii would respond to increased temperature and changing precipitation patterns. Genetic factors played a more significant role in stomatal size than environmental factors, suggesting that B. ermanii struggles to adjust its stomatal size in response to a changing environment. Our results also revealed a positive correlation between stomatal size plasticity and original habitat suitability, indicating that in B. ermanii populations in harsh environments exhibit lower adaptability to environmental shifts. Although stomatal density and size of B. ermanii showed the significant responses to increased temperature and shifting precipitation patterns, the response ranges of stomatal density and size to the environmental factors varied among populations. Our findings highlighted the interplay between genetic and environmental factors in determining the intraspecific variation in stomatal density and size in B. ermanii. This indicated that certain populations of B. ermanii exhibit limited stomatal plasticity and adaptability, which could directly affect photosynthesis and transpiration, suggesting potential population-specific fitness implications for B. ermanii under future climate change.

Climate, not Aboriginal landscape burning, controlled the historical demography and distribution of fire-sensitive conifer populations across Australia.

Climate and fire are the key environmental factors that shape the distribution and demography of plant populations in Australia. Because of limited palaeoecological records in this arid continent, however, it is unclear as to which factor impacted vegetation more strongly, and what were the roles of fire regime changes owing to human activity and megafaunal extinction (since ca 50 kya). To address these questions, we analysed historical genetic, demographic and distributional changes in a widespread conifer species complex that paradoxically grows in fire-prone regions, yet is very sensitive to fire. Genetic demographic analysis showed that the arid populations experienced strong bottlenecks, consistent with range contractions during the Last Glacial Maximum (ca 20 kya) predicted by species distribution models. In southern temperate regions, the population sizes were estimated to have been mostly stable, followed by some expansion coinciding with climate amelioration at the end of the last glacial period. By contrast, in the flammable tropical savannahs, where fire risk is the highest, demographic analysis failed to detect significant population bottlenecks. Collectively, these results suggest that the impact of climate change overwhelmed any modifications to fire regimes by Aboriginal landscape burning and megafaunal extinction, a finding that probably also applies to other fire-prone vegetation across Australia.

Nuclear and chloroplast DNA phylogeography reveals Pleistocene divergence and subsequent secondary contact of two genetic lineages of the tropical rainforest tree species Shorea leprosula (Dipterocarpaceae) in South-East Asia.

Tropical rainforests in South-East Asia have been affected by climatic fluctuations during past glacial eras. To examine how the accompanying changes in land areas and temperature have affected the genetic properties of rainforest trees in the region, we investigated the phylogeographic patterns of a widespread dipterocarp species, Shorea leprosula. Two types of DNA markers were used: expressed sequence tag-based simple sequence repeats and chloroplast DNA (cpDNA) sequence variations. Both sets of markers revealed clear genetic differentiation between populations in Borneo and those in the Malay Peninsula and Sumatra (Malay/Sumatra). However, in the south-western part of Borneo, genetic admixture of the lineages was observed in the two marker types. Coalescent simulation based on cpDNA sequence variation suggested that the two lineages arose 0.28-0.09 million years before present and that following their divergence migration from Malay/Sumatra to Borneo strongly exceeded migration in the opposite direction. We conclude that the genetic structure of S. leprosula was largely formed during the middle Pleistocene and was subsequently modified by eastward migration across the subaerially exposed Sunda Shelf.

Comparing modeling methods of genomic prediction for growth traits of a tropical timber species, Shorea macrophylla.

Introduction: Shorea macrophylla is a commercially important tropical tree species grown for timber and oil. It is amenable to plantation forestry due to its fast initial growth. Genomic selection (GS) has been used in tree breeding studies to shorten long breeding cycles but has not previously been applied to S. macrophylla. Methods: To build genomic prediction models for GS, leaves and growth trait data were collected from a half-sib progeny population of S. macrophylla in Sari Bumi Kusuma forest concession, central Kalimantan, Indonesia. 18037 SNP markers were identified in two ddRAD-seq libraries. Genomic prediction models based on these SNPs were then generated for diameter at breast height and total height in the 7th year from planting (D7 and H7). Results and discussion: These traits were chosen because of their relatively high narrow-sense genomic heritability and because seven years was considered long enough to assess initial growth. Genomic prediction models were built using 6 methods and their derivatives with the full set of identified SNPs and subsets of 48, 96, and 192 SNPs selected based on the results of a genome-wide association study (GWAS). The GBLUP and RKHS methods gave the highest predictive ability for D7 and H7 with the sets of selected SNPs and showed that D7 has an additive genetic architecture while H7 has an epistatic genetic architecture. LightGBM and CNN1D also achieved high predictive abilities for D7 with 48 and 96 selected SNPs, and for H7 with 96 and 192 selected SNPs, showing that gradient boosting decision trees and deep learning can be useful in genomic prediction. Predictive abilities were higher in H7 when smaller number of SNP subsets selected by GWAS p-value was used, However, D7 showed the contrary tendency, which might have originated from the difference in genetic architecture between primary and secondary growth of the species. This study suggests that GS with GWAS-based SNP selection can be used in breeding for non-cultivated tree species to improve initial growth and reduce genotyping costs for next-generation seedlings.

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.

BACKGROUND: High-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. RESULTS: We 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. CONCLUSIONS: Using 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.

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.

BACKGROUND: Microsatellites 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). RESULTS: We 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. CONCLUSIONS: We 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.

Genetic Structure of the Liriope muscari Polyploid Complex and the Possibility of Its Genetic Disturbance in Japan.

Anthropogenic activities, such as the movement of plants through greening, can result in genetic disturbance that can interfere with local adaptation in wild populations. Although research is underway to prevent genetic disturbance associated with greening, genetic disturbance of intraspecific polyploidy, which is estimated to be present in 24% of vascular plants, has not been well studied. Liriope muscari is a polyploid complex with known diploid (2n = 36), tetraploid (2n = 72), and hexaploid (2n = 108) forms. The plants of this species tolerate dry and hot conditions and are therefore frequently used for greening and gardening. However, the distribution of this polyploid in Japan, its genetic structure, and genetic disturbance are not known. In this study, we investigated the polyploidy distribution and genetic structure in naturally distributed L. muscari in Japan using chloroplast DNA (cpDNA) haplotypes and nuclear DNA (nDNA). Commercially produced individuals were also studied and compared with natural populations to assess any genetic disturbance of the ploidy complex in this species. Chromosome counts, cpDNA, and nDNA results showed three genetically and cytologically distinct groups in Japan: first, a tetraploid group in mainland Japan; second, a hexaploid group in the Ryukyu Islands; and third, a diploid and tetraploid group in the Ryukyu Islands. Significant isolation by distance was also detected within the three groups (p = 0.001). Genetic disturbance due to greening and gardening should be avoided among the three groups. Genetic disturbance can be reduced by using individuals derived from natural populations that are close to the sites used for greening and gardening. For commercially produced individuals, genetic disturbance is unlikely in the Kanto region, an area of high usage, while genetic disturbance is thought possible in the Ryukyu Islands.

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.

QTL analysis of heterostyly in Primula sieboldii and its application for morph identification in wild populations.

BACKGROUND AND AIMS: Primula sieboldii is a perennial clonal herb that is distributed around the Sea of Japan and is endangered in Japan. Its breeding system is characterized by heteromorphic self-incompatibility, and the morph ratio within a population is very important for reproductive success. The aims of this study were to construct a linkage map, map the S locus as a qualitative trait and quantitative trait loci (QTLs) for floral morphological traits related to heterostyly, and predict the morph type in wild populations by using molecular markers for devising a conservation strategy. METHODS: A linkage map was constructed with 126 markers. The QTLs for four floral traits and the S locus were mapped. Using the genotypes of loci that were located near both the S locus and the QTLs with large effects, morphs of 59 wild genets were predicted. KEY RESULTS: The linkage map consisted of 14 linkage groups (LGs). The S locus was mapped to LG 7. Major QTLs for stigma and anther heights were detected in the same region as the S locus. These QTLs exhibited high logarithm of the odds scores and explained a high percentage of the phenotypic variance (>85 %). By analysing these two traits within each morph, additional QTLs for each trait were detected. Using the four loci linked to the S locus, the morphs of 43 genets in three wild populations could be predicted. CONCLUSIONS: This is the first report of a linkage map and QTL analysis for floral morphology related to heterostyly in P. sieboldii. Floral morphologies related to heterostyly are controlled by the S locus in LG 7 and by several QTLs in other LGs. Additionally, this study showed that molecular markers are effective tools for investigating morph ratios in a population containing the non-flowering individuals or during the non-flowering seasons.