Development of EST-SSR markers from an inner bark cDNA library of Fagus crenata (Fagaceae)

A novel synthesis of two decades of microsatellite studies on European beech reveals decreasing genetic diversity from glacial refugia

Genetic diversity influences the evolutionary potential of forest trees under changing environmental conditions, thus indirectly the ecosystem services that forests provide. European beech (Fagus sylvatica L.) is a dominant European forest tree species that increasingly suffers from climate change-related die-back. Here, we conducted a systematic literature review of neutral genetic diversity in European beech and created a meta-data set of expected heterozygosity (He) from all past studies providing nuclear microsatellite data. We propose a novel approach, based on population genetic theory and a min-max scaling to make past studies comparable. Using a new microsatellite data set with unprecedented geographic coverage and various re-sampling schemes to mimic common sampling biases, we show the potential and limitations of the scaling approach. The scaled meta-dataset reveals the expected trend of decreasing genetic diversity from glacial refugia across the species range and also supports the hypothesis that different lineages met and admixed north of the European mountain ranges. As a result, we present a map of genetic diversity across the range of European beech which could help to identify seed source populations harboring greater diversity and guide sampling strategies for future genome-wide and functional investigations of genetic variation. Our approach illustrates how to combine information from several nuclear microsatellite data sets to describe patterns of genetic diversity extending beyond the geographic scale or mean number of loci used in each individual study, and thus is a proof-of-concept for synthesizing knowledge from existing studies also in other species.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11295-022-01577-4.

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

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

Comparative in silico analysis of EST-SSRs in angiosperm and gymnosperm tree genera

BACKGROUND

Simple Sequence Repeats (SSRs) derived from Expressed Sequence Tags (ESTs) belong to the expressed fraction of the genome and are important for gene regulation, recombination, DNA replication, cell cycle and mismatch repair. Here, we present a comparative analysis of the SSR motif distribution in the 5'UTR, ORF and 3'UTR fractions of ESTs across selected genera of woody trees representing gymnosperms (17 species from seven genera) and angiosperms (40 species from eight genera).

RESULTS

Our analysis supports a modest contribution of EST-SSR length to genome size in gymnosperms, while EST-SSR density was not associated with genome size in neither angiosperms nor gymnosperms. Multiple factors seem to have contributed to the lower abundance of EST-SSRs in gymnosperms that has resulted in a non-linear relationship with genome size diversity. The AG/CT motif was found to be the most abundant in SSRs of both angiosperms and gymnosperms, with a relative increase in AT/AT in the latter. Our data also reveals a higher abundance of hexamers across the gymnosperm genera.

CONCLUSIONS

Our analysis provides the foundation for future comparative studies at the species level to unravel the evolutionary processes that control the SSR genesis and divergence between angiosperm and gymnosperm tree species.

Genetic differentiation and genetic diversity of Castanopsis (Fagaceae), the dominant tree species in Japanese broadleaved evergreen forests, revealed by analysis of EST-associated microsatellites

The broadleaved evergreen forests of the East Asian warm temperate zone are characterised by their high biodiversity and endemism, and there is therefore a need to extend our understanding of its genetic diversity and phylogeographic patterns. Castanopsis (Fagaceae) is one of the dominant tree species in the broadleaved evergreen forests of Japan. In this study we investigate the genetic diversity, genetic structure and leaf epidermal morphology of 63 natural populations of C. sieboldii and C. cuspidata, using 32 Expressed Sequence Tag associated microsatellites. The overall genetic differentiation between populations was low (G_ST = 0.069 in C. sieboldii and G_ST = 0.057 in C. cuspidata). Neighbor-joining tree and Bayesian clustering analyses revealed that the populations of C. sieboldii and C. cuspidata were genetically clearly differentiated, a result which is consistent with the morphology of their epidermal cell layers. This suggests that C. sieboldii and C. cuspidata should be treated as independent species, although intermediate morphologies are often observed, especially at sites where the two species coexist. The higher level of genetic diversity observed in the Kyushu region (for both species) and the Ryukyu Islands (for C. sieboldii) is consistent with the available fossil pollen data for Castanopsis-type broadleaved evergreen trees during the Last Glacial Maximum and suggests the existence of refugia for Castanopsis forests in southern Japan. Within the C. sieboldii populations, Bayesian clustering analyses detected three clusters, in the western and eastern parts of the main islands and in the Ryukyu Islands. The west-east genetic differentiation observed for this species in the main islands, a pattern which is also found in several plant and animal species inhabiting Castanopsis forests in Japan, suggests that they have been isolated from each other in the western and eastern populations for an extended period of time, and may imply the existence of eastern refugia.

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.

Phylogeography of Phytophagous Weevils and Plant Species in Broadleaved Evergreen Forests: A Congruent Genetic Gap between Western and Eastern Parts of Japan

The Quaternary climate cycles played an important role in shaping the distribution of biodiversity among current populations, even in warm-temperate zones, where land was not covered by ice sheets. We focused on the Castanopsis-type broadleaved evergreen forest community in Japan, which characterizes the biodiversity and endemism of the warm-temperate zone. A comparison of the phylogeographic patterns of three types of phytophagous weevils associated with Castanopsis (a host-specific seed predator, a generalist seed predator, and a host-specific leaf miner) and several other plant species inhabiting the forests revealed largely congruent patterns of genetic differentiation between western and eastern parts of the main islands of Japan. A genetic gap was detected in the Kii Peninsula to Chugoku-Shikoku region, around the Seto Inland Sea. The patterns of western-eastern differentiation suggest past fragmentation of broadleaved evergreen forests into at least two separate refugia consisting of the southern parts of Kyushu to Shikoku and of Kii to Boso Peninsula. Moreover, the congruent phylogeographic patterns observed in Castanopsis and the phytophagous insect species imply that the plant-herbivore relationship has been largely maintained since the last glacial periods. These results reinforce the robustness of the deduced glacial and postglacial histories of Castanopsis-associated organisms.