Using reduced representation libraries sequencing methods to identify cpDNA polymorphisms in European beech (Fagus sylvatica L)

Reduced representation genomic libraries (RRLs) are increasingly used to answer diverse questions in evolutionary biology, which remained unresolved otherwise. In the case of European beech (Fagus sylvatica L.), former applications of nuclear and chloroplast genetic markers indicated that most of the beech populations of Central, Eastern, and Northern Europe have a very homogeneous genetic structure. Thus, determining differentiation among populations of F. sylvatica at regional scale in Central Europe has remained a challenge for evolutionary biologists, mainly due to lack of variable genetic markers. In this study, we used the RAD-seq and GBS approaches to identify novel cpDNA polymorphisms that could be helpful to identify population structure of European beech at geographical scale. Our datasets allowed us to indicate a finer population structure and stronger patterns of spatial genetic structure than previous studies. It is expected that the applications of cpDNA-based phylogenetic approaches using whole-genome sequencing are able to provide even more detailed picture of species phylogeography than ever before

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

Comparative Analysis of SNP Discovery and Genotyping in <i>Fagus sylvatica</i> L. and <i>Quercus robur</i> L. Using RADseq, GBS, and ddRAD Methods

Next-generation sequencing of reduced representation genomic libraries (RRL) is capable of providing large numbers of genetic markers for population genetic studies at relatively low costs. However, one major concern of these types of markers is the precision of genotyping, which is related to the common problem of missing data, which appears to be particularly important in association and genomic selection studies. We evaluated three RRL approaches (GBS, RADseq, ddRAD) and different SNP identification methods (de novo or based on a reference genome) to find the best solutions for future population genomics studies in two economically and ecologically important broadleaved tree species, namely F. sylvatica and Q. robur. We found that the use of ddRAD method coupled with SNP calling based on reference genomes provided the largest numbers of markers (28 k and 36 k for beech and oak, respectively), given standard filtering criteria. Using technical replicates of samples, we demonstrated that more than 80% of SNP loci should be considered as reliable markers in GBS and ddRAD, but not in RADseq data. According to the reference genomes’ annotations, more than 30% of the identified ddRAD loci appeared to be related to genes. Our findings provide a solid support for using ddRAD-based SNPs for future population genomics studies in beech and oak

Phylogeny of Maleae (Rosaceae) Based on Complete Chloroplast Genomes Supports the Distinction of <em>Aria</em>, <em>Chamaemespilus</em> and <em>Torminalis</em> as Separate Genera, Different from <em>Sorbus</em> sp.

Several genera formerly contained within the genus Sorbus L. sensu lato have been proposed as separate taxa, including Aria, Chamaemespilus and Torminalis. However, molecular evidence for such distinctions are rather scarce. We assembled the complete chloroplast genome of Sorbus aucuparia, another representative of Sorbus s.s., and performed detailed comparisons with the available genomes of Aria edulis, Chamaemespilus alpina and Torminalis glaberrima. Additionally, using 110 complete chloroplast genomes of the Maleae representatives, we constructed the phylogenetic tree of the tribe using Maximum Likelihood methods. The chloroplast genome of S. aucuparia was found to be similar to other species within Maleae. The phylogenetic tree of the Maleae tribe indicated that A. edulis, C. alpina and T. glaberrima formed a concise group belonging to a different clade (related to Malus) than the one including Sorbus s.s. (related to Pyrus). However, Aria and Chamaemespilus appeared to be more closely related to each other than to Torminalis. Our results provide additional support for considering Aria, Chamaemespilus and Torminalis as separate genera different from Sorbus s.s

Complete chloroplast genomes of Fagus sylvatica L. reveal sequence conservation in the inverted repeat and the presence of allelic variation in NUPTs

Growing amounts of genomic data and more efficient assembly tools advance organelle genomics at an unprecedented scale. Genomic resources are increasingly used for phylogenetic analyses of many plant species, but are less frequently used to investigate within-species variability and phylogeography. In this study, we investigated genetic diversity of Fagus sylvatica, an important broadleaved tree species of European forests, based on complete chloroplast genomes of 18 individuals sampled widely across the species distribution. Our results confirm the hypothesis of a low cpDNA diversity in European beech. The chloroplast genome size was remarkably stable (158,428 ± 37 bp). The polymorphic markers, 12 microsatellites (SSR), four SNPs and one indel, were found only in the single copy regions, while inverted repeat regions were monomorphic both in terms of length and sequence, suggesting highly efficient suppression of mutation. The within-individual analysis of polymorphisms showed >9k of markers which were proportionally present in gene and non-gene areas. However, an investigation of the frequency of alternate alleles revealed that the source of this diversity originated likely from nuclear-encoded plastome remnants (NUPTs). Phylogeographic and Mantel correlation analysis based on the complete chloroplast genomes exhibited clustering of individuals according to geographic distance in the first distance class, suggesting that the novel markers and in particular the cpSSRs could provide a more detailed picture of beech population structure in Central Europe

A circular chloroplast genome of fagus sylvatica reveals high conservation between two individuals from Germany and one individual from Poland and an alternate direction of the small single-copy region

Chloroplasts are difficult to assemble because of the presence of large inverted repeats. At the same time, correct assemblies are important, as chloroplast loci are frequently used for biogeography and population genetics studies. In an attempt to elucidate the orientation of the single-copy regions and to find suitable loci for chloroplast single nucleotide polymorphism (SNP)-based studies, circular chloroplast sequences for the ultra-centenary reference individual of European Beech (Fagus sylvatica), Bhaga, and an additional Polish individual (named Jamy) was obtained based on hybrid assemblies. The chloroplast genome of Bhaga was 158,458 bp, and that of Jamy was 158,462 bp long. Using long-read mapping on the configuration inferred in this study and the one suggested in a previous study, we found an inverted orientation of the small single-copy region. The chloroplast genome of Bhaga and of the individual from Poland both have only two mismatches as well as three and two indels as compared to the previously published genome, respectively. The low divergence suggests low seed dispersal but high pollen dispersal. However, once chloroplast genomes become available from Pleistocene refugia, where a high degree of variation has been reported, they might prove useful for tracing the migration history of Fagus sylvatica in the Holocene

A comparison of three circular mitochondrial genomes of fagus sylvatica from Germany and Poland reveals low variation and complete identity of the gene space

Similar to chloroplast loci, mitochondrial markers are frequently used for genotyping, phylogenetic studies, and population genetics, as they are easily amplified due to their multiple copies per cell. In a recent study, it was revealed that the chloroplast offers little variation for this purpose in central European populations of beech. Thus, it was the aim of this study to elucidate, if mitochondrial sequences might offer an alternative, or whether they are similarly conserved in central Europe. For this purpose, a circular mitochondrial genome sequence from the more than 300-year-old beech reference individual Bhaga from the German National Park Kellerwald-Edersee was assembled using long and short reads and compared to an individual from the Jamy Nature Reserve in Poland and a recently published mitochondrial genome from eastern Germany. The mitochondrial genome of Bhaga was 504,730 bp, while the mitochondrial genomes of the other two individuals were 15 bases shorter, due to seven indel locations, with four having more bases in Bhaga and three locations having one base less in Bhaga. In addition, 19 SNP locations were found, none of which were inside genes. In these SNP locations, 17 bases were different in Bhaga, as compared to the other two genomes, while 2 SNP locations had the same base in Bhaga and the Polish individual. While these figures are slightly higher than for the chloroplast genome, the comparison confirms the low degree of genetic divergence in organelle DNA of beech in central Europe, suggesting the colonisation from a common gene pool after the Weichsel Glaciation. The mitochondrial genome might have limited use for population studies in central Europe, but once mitochondrial genomes from glacial refugia become available, it might be suitable to pinpoint the origin of migration for the re-colonising beech population

A reference genome of the European beech (Fagus sylvatica L.)

Background: The European beech is arguably the most important climax broad-leaved tree species in Central Europe, widely planted for its valuable wood. Here, we report the 542 Mb draft genome sequence of an up to 300-year-old individual (Bhaga) from an undisturbed stand in the Kellerwald-Edersee National Park in central Germany. Findings: Using a hybrid assembly approach, Illumina reads with short- and long-insert libraries, coupled with long Pacific Biosciences reads, we obtained an assembled genome size of 542 Mb, in line with flow cytometric genome size estimation. The largest scaffold was of 1.15 Mb, the N50 length was 145 kb, and the L50 count was 983. The assembly contained 0.12% of Ns. A Benchmarking with Universal Single-Copy Orthologs (BUSCO) analysis retrieved 94% complete BUSCO genes, well in the range of other high-quality draft genomes of trees. A total of 62,012 protein-coding genes were predicted, assisted by transcriptome sequencing. In addition, we are reporting an efficient method for extracting high-molecular-weight DNA from dormant buds, by which contamination by environmental bacteria and fungi was kept at a minimum. Conclusions: The assembled genome will be a valuable resource and reference for future population genomics studies on the evolution and past climate change adaptation of beech and will be helpful for identifying genes, e.g., involved in drought tolerance, in order to select and breed individuals to adapt forestry to climate change in Europe. A continuously updated genome browser and download page can be accessed from beechgenome.net, which will include future genome versions of the reference individual Bhaga, as new sequencing approaches develop

Supporting Information from Relatives reproduce in synchrony: kinship and individual condition shape intraspecific variation in masting phenotype

Masting, synchronous and interannually variable seed production, is frequently called a reproductive strategy. Yet, it is unclear whether the reproductive behaviour of individuals has a heritable component. To address this, we used 22 years of annual fruit production data from 110 Sorbus aucuparia L. trees to examine the contributions of genetic factors to the reproductive phenotype of individuals, while controlling for environmental variation. Trees sharing close genetic relationships and experiencing similar habitat conditions exhibited similar levels of reproductive synchrony. Trees of comparable sizes displayed similar levels of year-to-year variation in fruiting, with relatedness contributing to this variation. External factors, such as shading, influenced the time intervals between years with abundant fruit production. The effects of genetic relatedness on the synchrony of reproduction among trees and on interannual variation provide long-awaited evidence that the masting phenotype is heritable, and can respond to natural selection