Identifying and testing marker-trait associations for growth and phenology in three pine species: Implications for genomic prediction.

In tree species, genomic prediction offers the potential to forecast mature trait values in early growth stages, if robust marker-trait associations can be identified. Here we apply a novel multispecies approach using genotypes from a new genotyping array, based on 20,795 single nucleotide polymorphisms (SNPs) from three closely related pine species (Pinus sylvestris, Pinus uncinata and Pinus mugo), to test for associations with growth and phenology data from a common garden study. Predictive models constructed using significantly associated SNPs were then tested and applied to an independent multisite field trial of P. sylvestris and the capability to predict trait values was evaluated. One hundred and eighteen SNPs showed significant associations with the traits in the pine species. Common SNPs (MAF > 0.05) associated with bud set were only found in genes putatively involved in growth and development, whereas those associated with growth and budburst were also located in genes putatively involved in response to environment and, to a lesser extent, reproduction. At one of the two independent sites, the model we developed produced highly significant correlations between predicted values and observed height data (YA, height 2020: r = 0.376, p < 0.001). Predicted values estimated with our budburst model were weakly but positively correlated with duration of budburst at one of the sites (GS, 2015: r = 0.204, p = 0.034; 2018: r = 0.205, p = 0.034-0.037) and negatively associated with budburst timing at the other (YA: r = -0.202, p = 0.046). Genomic prediction resulted in the selection of sets of trees whose mean height was taller than the average for each site. Our results provide tentative support for the capability of prediction models to forecast trait values in trees, while highlighting the need for caution in applying them to trees grown in different environments.

Taming the massive genome of Scots pine with PiSy50k, a new genotyping array for conifer research.

Pinus sylvestris (Scots pine) is the most widespread coniferous tree in the boreal forests of Eurasia, with major economic and ecological importance. However, its large and repetitive genome presents a challenge for conducting genome-wide analyses such as association studies, genetic mapping and genomic selection. We present a new 50K single-nucleotide polymorphism (SNP) genotyping array for Scots pine research, breeding and other applications. To select the SNP set, we first genotyped 480 Scots pine samples on a 407 540 SNP screening array and identified 47 712 high-quality SNPs for the final array (called 'PiSy50k'). Here, we provide details of the design and testing, as well as allele frequency estimates from the discovery panel, functional annotation, tissue-specific expression patterns and expression level information for the SNPs or corresponding genes, when available. We validated the performance of the PiSy50k array using samples from Finland and Scotland. Overall, 39 678 (83.2%) SNPs showed low error rates (mean = 0.9%). Relatedness estimates based on array genotypes were consistent with the expected pedigrees, and the level of Mendelian error was negligible. In addition, array genotypes successfully discriminate between Scots pine populations of Finnish and Scottish origins. The PiSy50k SNP array will be a valuable tool for a wide variety of future genetic studies and forestry applications.

Candidate Genes for the High-Altitude Adaptations of Two Mountain Pine Taxa.

Mountain plants, challenged by vegetation time contractions and dynamic changes in environmental conditions, developed adaptations that help them to balance their growth, reproduction, survival, and regeneration. However, knowledge regarding the genetic basis of species adaptation to higher altitudes remain scarce for most plant species. Here, we attempted to identify such corresponding genomic regions of high evolutionary importance in two closely related European pines, Pinus mugo and P. uncinata, contrasting them with a reference lowland relative-P. sylvestris. We genotyped 438 samples at thousands of single nucleotide polymorphism (SNP) markers, tested their genetic differentiation and population structure followed by outlier detection and gene ontology annotations. Markers clearly differentiated the species and uncovered patterns of population structure in two of them. In P. uncinata three Pyrenean sites were grouped together, while two outlying populations constituted a separate cluster. In P. sylvestris, Spanish population appeared distinct from the remaining four European sites. Between mountain pines and the reference species, 35 candidate genes for altitude-dependent selection were identified, including such encoding proteins responsible for photosynthesis, photorespiration and cell redox homeostasis, regulation of transcription, and mRNA processing. In comparison between two mountain pines, 75 outlier SNPs were found in proteins involved mainly in the gene expression and metabolism.

Hybridization and introgression of native and foreign Sorbus tree species in unique environments of protected mountainous areas.

Hybridization and introgression are important processes influencing the genetic diversity and evolution of species. These processes are of particular importance in protected areas, where they can lead to the formation of hybrids between native and foreign species and may ultimately result in the loss of parental species from their natural range. Despite their importance, the contribution of hybridization and introgression to genetic diversity in Sorbus genus remains not fully recognized. We analysed the genetic and morphological variability of several Sorbus species including native (Sorbus aria), foreign (S. intermedia) and potentially hybrid (S. carpatica) individuals from the Polish Carpathian range. Patterns of variation at 13 nuclear microsatellite loci show hybridization between the tested species and confirm the existence of the hybrid form S. carpatica. Biometric analysis on leaves, based of 10 metric features and three parameters, identified several characters for preliminary taxonomic classification; however, none of them could be used as a fully diagnostic marker for faultless annotation of S. intermedia and S. carpatica. The genetic structure analysis indicated complex patterns of population differentiation and its diverse origin. The results allow assessment of genetic variation and identification of parental species participating in hybridization. This knowledge will advance the management of genetic diversity and development of conservation strategies for efficient maintenance of the unique protected ecosystem.

Development of a single nucleotide polymorphism array for population genomic studies in four European pine species.

Pines are some of the most ecologically and economically important tree species in the world, and many have enormous natural distributions or have been extensively planted. However, a lack of rapid genotyping capability is hampering progress in understanding the molecular basis of genetic variation in these species. Here, we deliver an efficient tool for genotyping thousands of single nucleotide polymorphism (SNP) markers across the genome that can be applied to genetic studies in pines. Polymorphisms from resequenced candidate genes and transcriptome sequences of P. sylvestris, P. mugo, P. uncinata, P. uliginosa and P. radiata were used to design a 49,829 SNP array (Axiom_PineGAP, Thermo Fisher). Over a third (34.68%) of the unigenes identified from the P. sylvestris transcriptome were represented on the array, which was used to screen samples of four pine species. The conversion rate for the array on all samples was 42% (N = 20,795 SNPs) and was similar for SNPs sourced from resequenced candidate gene and transcriptome sequences. The broad representation of gene ontology terms by unigenes containing converted SNPs reflected their coverage across the full transcriptome. Over a quarter of successfully converted SNPs were polymorphic among all species, and the data were successful in discriminating among the species and some individual populations. The SNP array provides a valuable new tool to advance genetic studies in these species and demonstrates the effectiveness of the technology for rapid genotyping in species with large and complex genomes.

Patterns of mtDNA variation reveal complex evolutionary history of relict and endangered peat bog pine (Pinus uliginosa).

Estimates of genetic differentiation at intra- and interspecific level are often hindered by the lack of suitable molecular markers. Low phylogeographic resolution limits development of appropriate conservation strategies especially in case of endangered forest tree species with small and disjunct distribution. In this study, we assessed fine-scale genetic structure of relict and endangered peat bog pine (Pinus uliginosa) and two other closely related European pine species (Pinus mugo and Pinus uncinata) using a set of 15 newly developed maternally inherited and seed-mediated mitochondrial DNA (mtDNA) markers and two previously known polymorphic mtDNA regions (nad1, nad7). Three main groups, corresponding in general to three investigated species were revealed in the haplotype network analysis. However, only P. uncinata was clearly distinct at all levels of analysis, whereas great genetic similarity and haplotype sharing was observed between P. uliginosa and P. mugo. Strong phylogeographic structure was found in P. uliginosa that showed high differentiation at relatively short geographical distance among populations and the existence of mitochondrial lineages of different evolutionary history. Hybridization with other pine species has likely contributed to genetic differentiation of P. uliginosa as indicated by contemporary distribution of mtDNA haplotypes. The research emphasizes the importance of accurate assessments of genetic structure of endangered species with complex evolutionary history for development of efficient conservation strategies.

Early phenology and growth trait variation in closely related European pine species.

Closely related taxa occupying different environments are valuable systems for studying evolution. In this study, we examined differences in early phenology (bud set, bud burst) and early growth in a common garden trial of closely related pine species: Pinus sylvestris, P. mugo, and P. uncinata. Seeds for the trial were sourced from populations across the ranges of each species in Europe. Over first 4 years of development, clear differences were observed between species, while the most significant intraspecific differentiation was observed among plants from P. sylvestris populations from continental European locations. Trait differences within P. sylvestris were highly correlated with altitude and latitude of the site of origin. Meanwhile, P. mugo populations from the Carpathians had the earliest bud set and bud flush compared to other populations of the species. Overall, populations from the P. mugo complex from heterogeneous mountain environments and P. sylvestris from the Scottish Highlands showed the highest within-population variation for the focal traits. Although the three species have been shown to be genetically highly similar, this study reveals large differences in key adaptive traits both among and within species.

Molecular signatures of divergence and selection in closely related pine taxa.

Efforts to detect loci under selection in plants have mostly focussed on single species. However, assuming that intraspecific divergence may lead to speciation, comparisons of genetic variation within and among recently diverged taxa can help to locate such genes. In this study, coalescent and outlier detection methods were used to assess nucleotide polymorphism and divergence at 79 nuclear gene fragments (1212 SNPs) in 16 populations (153 individuals) of the closely related, but phenotypically and ecologically distinct, pine taxa Pinus mugo, P. uliginosa and P. uncinata across their European distributions. Simultaneously, mitochondrial DNA markers, which are maternally inherited in pines and distributed by seeds at short geographic distance, were used to assess genetic relationships of the focal populations and taxa. The majority of nuclear loci showed homogenous patterns of variation between the taxa due to a high number of shared SNPs and haplotypes, similar levels of polymorphism, and low net divergence. However, against this common genetic background and an overall low population structure within taxa at mitochondrial markers, we identified several genes showing signatures of selection, accompanied by significant intra- and interspecific divergence. Our results indicate that loci involved in species divergence may be involved in intraspecific local adaptation.

Reconstructing the plant mitochondrial genome for marker discovery: a case study using Pinus.

Whole-genome-shotgun (WGS) sequencing of total genomic DNA was used to recover ~1 Mbp of novel mitochondrial (mtDNA) sequence from Pinus sylvestris (L.) and three members of the closely related Pinus mugo species complex. DNA was extracted from megagametophyte tissue from six mother trees from locations across Europe, and 100-bp paired-end sequencing was performed on the Illumina HiSeq platform. Candidate mtDNA sequences were identified by their size and coverage characteristics, and by comparison with published plant mitochondrial genomes. Novel variants were identified, and primers targeting these loci were trialled on a set of 28 individuals from across Europe. In total, 31 SNP loci were successfully resequenced, characterizing 15 unique haplotypes. This approach offers a cost-effective means of developing marker resources for mitochondrial genomes in other plant species where reference sequences are unavailable.

Contrasting patterns of genetic variation in core and peripheral populations of highly outcrossing and wind pollinated forest tree species.

Gene flow tends to have a homogenising effect on a species' background genetic variation over large geographical areas. However, it is usually unknown to what extent the genetic structure of populations is influenced by gene exchange between core and peripheral populations that may represent stands of different evolutionary and demographic history. In this study, we looked at the patterns of population differentiation in Scots pine-a highly outcrossing and wind pollinated conifer species that forms large ecosystems of great ecological and economic importance in Europe and Asia. A set of 13 polymorphic nuclear microsatellite loci was analysed to infer the genetic relationships among 24 populations (676 individuals) from Europe and Asia Minor. The study included specimens from the primary continuous range and from isolated, marginal stands that are considered to be autochthonous populations representative of the species' putative refugial areas. Despite their presumably different histories, a similar level of genetic variation and no evidence of a population bottleneck was found across the populations. Differentiation among populations was relatively low (average FST = 0.035); however, the population structure was not homogenous, which was clearly evident from the allelic frequency spectra and Bayesian assignment analysis. Significant differentiation over short geographical distances was observed between isolated populations within the Iberian and Anatolian Peninsulas (Asia Minor), which contrasted with the absence of genetic differentiation observed between distant populations e.g., between central and northern Europe. The analysed populations were assigned to several groups that corresponded to the geographical regions of their occurrence. These results will be useful in genetics studies in Scots pine that aim to link nucleotide and phenotypic variation across the species distribution range and for development of sustainable breeding and management programs.

Comparative transcriptomics of a complex of four European pine species.

BACKGROUND: Pinus sylvestris, P. mugo, P. uliginosa and P. uncinata are closely related but phenotypically and ecologically very distinct European pine species providing an excellent study system for analysis of the genetic basis of adaptive variation and speciation. For comparative genomic analysis of the species, transcriptome sequence was generated for 17 samples collected across the European distribution range using Illumina paired-end sequencing technology. RESULTS: De novo transcriptome assembly of a reference sample of P. sylvestris contained 40968 unigenes, of which fewer than 0.5% were identified as putative retrotransposon sequences. Based on gene annotation approaches, 19659 contigs were identified and assigned to unique genes covering a broad range of gene ontology categories. About 80% of the reads from each sample were successfully mapped to the reference transcriptome of P. sylvestris. Single nucleotide polymorphisms were identified in 22041-24096 of the unigenes providing a set of ~220-262 k SNPs identified for each species. Very similar levels of nucleotide polymorphism were observed across species (π=0.0044-0.0053) and highest pairwise nucleotide divergence (0.006) was found between P. mugo and P. sylvestris at a common set of unigenes. CONCLUSIONS: The study provides whole transcriptome sequence and a large set of SNPs to advance population and association genetic studies in pines. Our study demonstrates that transcriptome sequencing can be a very useful approach for development of novel genomic resources in species with large and complex genomes.

Speciation history of three closely related pines Pinus mugo (T.), P. uliginosa (N.) and P. sylvestris (L.).

Nucleotide polymorphisms at genomic regions including 17 nuclear loci, two chloroplast and one mitochondrial DNA fragments were used to study the speciation history of three pine species: dwarf mountain pine (Pinus mugo), peat-bog pine (P. uliginosa) and Scots pine (P. sylvestris). We set out to investigate three specific speciation scenarios: (I) P. uliginosa is a homoploid hybrid between the other two, (II) the species have evolved without gene flow after divergence and (III) there has been substantial gene flow between the species since their divergence. Overall, the genetic data suggest that P. mugo and P. uliginosa share the same gene pool (average net divergence of 0.0001) and that the phenotypic differences (e.g. growth form) are most likely due to very limited areas of the genome. P. mugo and P. uliginosa are more diverged from P. sylvestris than from each other (average net divergence of 0.0027 and 0.0026, respectively). The nucleotide patterns can best be explained by the divergence with migration speciation scenario, although the hybrid speciation scenario with small genomic contribution from P. sylvestris cannot be completely ruled out. We suggest that the large amount of shared polymorphisms between the pine taxa and the lack of monophyly at all loci studied between P. sylvestris and P. mugo-P. uliginosa can largely be explained by relatively recent speciation history and large effective population sizes but also by interspecific gene flow. These closely related pine taxa form an excellent system for searching for loci involved in adaptive variation as they are differentiated in phenotype and ecology but have very similar genetic background.

Demographic history has influenced nucleotide diversity in European Pinus sylvestris populations.

To infer the role of natural selection in shaping standing genetic diversity, it is necessary to assess the genomewide impact of demographic history on nucleotide diversity. In this study we analyzed sequence diversity of 16 nuclear loci in eight Pinus sylvestris populations. Populations were divided into four geographical groups on the basis of their current location and the geographical history of the region: northern Europe, central Europe, Spain, and Turkey. There were no among-group differences in the level of silent nucleotide diversity, which was approximately 0.005/bp in all groups. There was some evidence that linkage disequilibrium extended further in northern Europe than in central Europe: the estimates of the population recombination rate parameter, rho, were 0.0064 and 0.0294, respectively. The summary statistics of nucleotide diversity in central and northern European populations were compatible with an ancient bottleneck rather than the standard neutral model.

Reciprocal controlled crosses between Pinus sylvestris and P. mugo verified by a species-specific cpDNA marker.

A species-specific marker of cpDNA (paternally inherited in pines) was used to verify the hybrid origin of seedlings from controlled reciprocal crosses between Pinus sylvestris and P. mugo. A very low degree of compatibility between those two species has been revealed. In the three consecutive years of experiments, no filled seeds were obtained in the combination with P. mugo as the seed parent. From P. sylvestris as the seed parent and P. mugo as the pollen donor, we succeeded to obtain four filled seeds (about 1 %), but only in one year. The seedling obtained from the seeds had cpDNA haplotypes specific to P. mugo, which proves their hybrid origin. This method enables verification of the result of controlled crosses. The importance of the results has been discussed in the aspect of postulated natural hybridisation in sympatric populations of the two species.