Clinal variation in PHY (PAS domain) and CRY (CCT domain)-Signs of local adaptation to light quality in Norway spruce.

Detection of the genomic basis of local adaptation to environmental conditions is challenging in forest trees. Phytochromes (PHY) and cryptochromes (CRY) perceive the red (R)/far-red (FR) and blue light respectively, thus playing a fundamental role in regulating plant growth and development. PHYO and PHYP from conifers are the equivalents of PHYA/PHYC and PHYB in angiosperms, respectively. Norway spruce shows an adaptive latitudinal cline for shade (low R:FR or FR-enriched light) tolerance and requirement of FR light for its growth. We analyzed the exome capture data that included a uniquely large data set of 1654 Norway spruce trees sampled across many latitudes in Sweden to capture the natural clines for photoperiod and FR light exposure during the growth season. Statistically significant clinal variation was detected in allele and genotype frequencies of missense mutations in coding regions belonging to well-defined functional domains of PHYO (PAS-B), PHYP2 (PAS fold-2), CRY1 (CCT1) and CRY2 (CCT2) that strongly correlates with the latitudinal gradient in response to variable light quality in Norway spruce. The missense SNP in PHYO resulting in Asn835Ser, displayed the steepest cline among all other polymorphisms. We propose that these variations in the photoreceptors represent signs of local adaptation to light quality.

Divergent selection predating the Last Glacial Maximum mainly acted on macro-phenotypes in Norway spruce.

The current distribution and population structure of many species were, to a large extent, shaped by cycles of isolation in glacial refugia and subsequent population expansions. Isolation in and postglacial expansion through heterogeneous environments led to either neutral or adaptive divergence. Norway spruce is no exception, and its current distribution is the consequence of a constant interplay between evolutionary and demographic processes. We investigated population differentiation and adaptation of Norway spruce for juvenile growth, diameter of the stem, wood density, and tracheid traits at breast height. Data from 4461 phenotyped and genotyped Norway spruce from 396 half-sib families in two progeny tests were used to test for divergent selection in the framework of Q ST vs. F ST. We show that the macroscopic resultant trait (stem diameter), unlike its microscopic components (tracheid dimensions) and juvenile growth, was under divergent selection that predated the Last Glacial Maximum. Altogether, the current variation in these phenotypic traits in Norway spruce is better explained by local adaptation to ancestral environments than to current ones, where populations were partly preadapted, mainly through growth-related traits.

Analysis of the ASR and LP3 homologous gene families reveal positive selection acting on LP3-3 gene.

Drought has long been established as a major environmental stress for plants which have in turn developed several coping strategies, ranging from physiological to molecular mechanisms. LP3 that was first discovered in loblolly pine (Pinus taeda L.) is a homolog of the Abscisic Acid, Stress and Ripening (ASR) gene belonging to the ABA/WDS gene family that was first detected in tomato. LP3 has been shown to be present in four different paralogs in loblolly pine called LP3-0, LP3-1, LP3-2 and LP3-3. LP3 in loblolly pine has not been as extensively studied as the ASR in tomato. Similar to ASR, the different LP3 paralogs have been shown to be upregulated in response to water deficit stress and to act as transcription factors for genes likely involved in hexose transport. In the current study, we have investigated the evolutionary history of LP3 gene family, with the aim of relating it to that of ASR from a phylogenetic perspective and comparing the differences in selective pressure and codon usage. Phylogenetic trees revealed that LP3 is less divergent across species than ASR even when the trees were solely based on the different sub-sections of the gene. Phylogenetic, GC content, codon usage and selective pressure analyses suggest that LP3-3 is undergoing positive selection.

Adaptive strategies of Scots pine under shade: Increase in lignin synthesis and ecotypic variation in defense-related gene expression.

Shade is a stressful condition for plants characterized by low Red:Far-Red (R:FR) ratio. The northern latitudes in Sweden daily receive more hours of FR-enriched light (twilight) or shade-like conditions compared to southern forests during the growing season. Scots pine (Pinus sylvestris L.) is a shade-intolerant species. Yet, it is well adapted to this latitudinal variation in light, which is evident by a northward increase in FR requirement to maintain growth. Shade adversely affects plant growth; it makes the plant weak and, therefore, susceptible to pathogen attack. Lignin is involved in plant protection against pathogen invasion mainly by forming a physical barrier. We studied lignin synthesis and expression of defense-related genes (growth-defense trade-offs) under a low R:FR (shade) ratio in Scots pine. A higher number of immunity/defense-related genes were up-regulated in response to shade in northern populations compared to southern ones, which can be viewed as a local adaptation to light quality for optimal growth and survival. Light quality regulates lignin metabolism; light stimulates lignin synthesis, while shade causes a decrease in lignin synthesis in most angiosperms. In contrast, Scots pine shows an increase in lignin synthesis supported by the higher expression of a few key genes in the lignin biosynthetic pathway, a novel finding reported by our study. These findings can be applied to future breeding strategies in forestry to produce disease-resilient trees.

Teasing apart the joint effect of demography and natural selection in the birth of a contact zone.

Vast population movements induced by recurrent climatic cycles have shaped the genetic structure of plant species. During glacial periods species were confined to low-latitude refugia from which they recolonized higher latitudes as the climate improved. This multipronged recolonization led to many lineages that later met and formed large contact zones. We utilize genomic data from 5000 Picea abies trees to test for the presence of natural selection during recolonization and establishment of a contact zone in Scandinavia. Scandinavian P. abies is today made up of a southern genetic cluster originating from the Baltics, and a northern one originating from Northern Russia. The contact zone delineating them closely matches the limit between two major climatic regions. We show that natural selection contributed to its establishment and maintenance. First, an isolation-with-migration model with genome-wide linked selection fits the data better than a purely neutral one. Second, many loci show signatures of selection or are associated with environmental variables. These loci, regrouped in clusters on chromosomes, are often related to phenology. Altogether, our results illustrate how climatic cycles, recolonization and selection can establish strong local adaptation along contact zones and affect the genetic architecture of adaptive traits.

Genetic architecture behind developmental and seasonal control of tree growth and wood properties in Norway spruce.

Genetic control of tree growth and wood formation varies depending on the age of the tree and the time of the year. Single-locus, multi-locus, and multi-trait genome-wide association studies (GWAS) were conducted on 34 growth and wood property traits in 1,303 Norway spruce individuals using exome capture to cover ~130K single-nucleotide polymorphisms (SNPs). GWAS identified associations to the different wood traits in a total of 85 gene models, and several of these were validated in a progenitor population. A multi-locus GWAS model identified more SNPs associated with the studied traits than single-locus or multivariate models. Changes in tree age and annual season influenced the genetic architecture of growth and wood properties in unique ways, manifested by non-overlapping SNP loci. In addition to completely novel candidate genes, SNPs were located in genes previously associated with wood formation, such as cellulose synthases and a NAC transcription factor, but that have not been earlier linked to seasonal or age-dependent regulation of wood properties. Interestingly, SNPs associated with the width of the year rings were identified in homologs of Arabidopsis thaliana BARELY ANY MERISTEM 1 and rice BIG GRAIN 1, which have been previously shown to control cell division and biomass production. The results provide tools for future Norway spruce breeding and functional studies.

Enhanced lignin synthesis and ecotypic variation in defense-related gene expression in response to shade in Norway spruce.

During the growth season, northern forests in Sweden daily receive more hours of far-red (FR)-enriched light or twilight (shade) as compared to southern forests. Norway spruce (shade-tolerant) are adapted to latitudinal variation in twilight characterized by a northward increase in FR requirement to maintain growth. Shade is a stressful condition that affects plant growth and increases plant's susceptibility to pathogen attack. Lignin plays a central role in plant defense and its metabolism is regulated by light wavelength composition (light quality). In the current work, we studied regulation of lignin synthesis and defense-related genes (growth-defense trade-offs) in response to shade in Norway spruce. In most angiosperms, light promotes lignin synthesis, whereas shade decreases lignin production leading to weaker stem, which may make plants more disease susceptible. In contrast, enhanced lignin synthesis was detected in response to shade in Norway spruce. We detected a higher number of immunity/defense-related genes up-regulated in northern populations as compared to south ones in response to shade. Enhanced lignin synthesis coupled with higher defense-related gene expression can be interpreted as an adaptive strategy for better survival in northern populations. Findings will contribute to ensuring deployment of well-adapted genetic material and identifying tree families with enhanced disease resistance.

Leveraging breeding programs and genomic data in Norway spruce (Picea abies L. Karst) for GWAS analysis.

BACKGROUND: Genome-wide association studies (GWAS) identify loci underlying the variation of complex traits. One of the main limitations of GWAS is the availability of reliable phenotypic data, particularly for long-lived tree species. Although an extensive amount of phenotypic data already exists in breeding programs, accounting for its high heterogeneity is a great challenge. We combine spatial and factor-analytics analyses to standardize the heterogeneous data from 120 field experiments of 483,424 progenies of Norway spruce to implement the largest reported GWAS for trees using 134 605 SNPs from exome sequencing of 5056 parental trees. RESULTS: We identify 55 novel quantitative trait loci (QTLs) that are associated with phenotypic variation. The largest number of QTLs is associated with the budburst stage, followed by diameter at breast height, wood quality, and frost damage. Two QTLs with the largest effect have a pleiotropic effect for budburst stage, frost damage, and diameter and are associated with MAP3K genes. Genotype data called from exome capture, recently developed SNP array and gene expression data indirectly support this discovery. CONCLUSION: Several important QTLs associated with growth and frost damage have been verified in several southern and northern progeny plantations, indicating that these loci can be used in QTL-assisted genomic selection. Our study also demonstrates that existing heterogeneous phenotypic data from breeding programs, collected over several decades, is an important source for GWAS and that such integration into GWAS should be a major area of inquiry in the future.

Molecular signatures of local adaptation to light in Norway spruce.

MAIN CONCLUSION: Transcriptomic and exome capture analysis reveal an adaptive cline for shade tolerance in Norway spruce. Genes involved in the lignin pathway and immunity seem to play a potential role in contributing towards local adaptation to light. The study of natural variation is an efficient method to elucidate how plants adapt to local climatic conditions, a key process for the evolution of a species. Norway spruce is a shade-tolerant conifer in which the requirement of far-red light for growth increases latitudinally northwards. The objective of the study is to characterize the genetic control of local adaptation to light enriched in far-red in Norway spruce, motivated by a latitudinal gradient for the Red:Far-red (R:FR) ratio to which Norway spruce has been proven to be genetically adapted. We have established the genomic signatures of local adaptation by conducting transcriptomic (total RNA-sequencing) and genomic analyses (exome capture), for the identification of genes differentially regulated along the cline. RNA-sequencing revealed 274 differentially expressed genes in response to SHADE (low R:FR light), between the southern and northern natural populations in Sweden. Exome capture included analysis of a uniquely large data set (1654 trees) that revealed missense variations in coding regions of nine differentially expressed candidate genes, which followed a latitudinal cline in allele and genotype frequencies. These genes included five transcription factors involved in vital processes like bud-set/bud-flush, lignin pathway, and cold acclimation and other genes that take part in cell-wall remodeling, secondary cell-wall thickening, response to starvation, and immunity. Based on these results, we suggest that the northern populations might not only be able to adjust their growing season in response to low R:FR light, but they may also be better adapted towards disease resistance by up-regulation of the lignin pathway that is linked to immunity. This forms a concrete basis for local adaptation to light quality in Norway spruce, one of the most economically important conifer tree species in Sweden.

Evaluation of antimicrobial, antioxidant and cytotoxic properties of bioactive compounds produced from endophytic fungi of Himalayan yew ( Taxus wallichiana) in Nepal.

Background: Endophytic fungi are largely underexplored in the discovery of natural bioactive products though being rich sources of novel compounds with promising pharmaceutical potential. In this study, Taxus wallichiana, which has huge medicinal value, was investigated for its endophytic diversity and capability to produce bioactive secondary metabolites by analyzing antioxidant, antimicrobial and cytotoxic properties. Methods: The endophytes were identified by ITS-PCR using genomic DNA samples. The secondary metabolites were extracted by solvent extraction method using ethyl acetate. The antioxidant activity was analyzed by Thin Layer Chromatography, Total Phenol Content (TPC), Total Flavonoid Content (TFC) and DPPH assay, and the antimicrobial activity was analyzed by agar-well diffusion method. Brine shrimp lethality assay was used to analyze the cytotoxicity of the fungal extracts. Results: Out of 16 different Taxus trees sampled from different locations of Dhorpatan, 13 distinctive endophytic fungi were isolated and grouped into 9 different genera: Bjerkandera, Trichoderma, Preussia, Botrytis, Arthrinium, Alternaria, Cladosporium, Sporormiella and Daldinia. The ethyl acetate extracts isolated from three endophytic fungi: Alternaria alternata, Cladosporium cladosporioides and Alternaria brassicae showed significant TPC values of 204±6.144, 312.3±2.147 and 152.7±4.958µg GAE/mg of dry extract, respectively, and TFC values of 177.9±2.911, 644.1±4.202 and 96.38±3.851µg RE/mg of dry extract, respectively. Furthermore, these three extracts showed a dose dependent radical scavenging activity with IC 50 concentration of 22.85, 22.15 and 23.001 µg/ml, respectively. The extracts of C. cladosporioides and A. brassicae also showed promising antimicrobial activity against Escherichia coli, Staphylococcus aureus and Bacillus subtilis with a minimum inhibitory concentration of 250µg/ml for all bacteria. Both the samples showed cytotoxic property against shrimp nauplii with LC 50 of 104.2 and 125.9µg/ml, respectively. Conclusions: The crude fungal extracts obtained from endophytes: A. alternata, C. cladosporioides and A. brassicae upon purification and further identification of the bioactive compounds can be a fascinating source for novel pharmaceutical agents.

Genetic control of tracheid properties in Norway spruce wood.

Through the use of genome-wide association studies (GWAS) mapping it is possible to establish the genetic basis of phenotypic trait variation. Our GWAS study presents the first such effort in Norway spruce (Picea abies (L). Karst.) for the traits related to wood tracheid characteristics. The study employed an exome capture genotyping approach that generated 178 101 Single Nucleotide Polymorphisms (SNPs) from 40 018 probes within a population of 517 Norway spruce mother trees. We applied a least absolute shrinkage and selection operator (LASSO) based association mapping method using a functional multi-locus mapping approach, with a stability selection probability method as the hypothesis testing approach to determine significant Quantitative Trait Loci (QTLs). The analysis has provided 30 significant associations, the majority of which show specific expression in wood-forming tissues or high ubiquitous expression, potentially controlling tracheids dimensions, their cell wall thickness and microfibril angle. Among the most promising candidates based on our results and prior information for other species are: Picea abies BIG GRAIN 2 (PabBG2) with a predicted function in auxin transport and sensitivity, and MA_373300g0010 encoding a protein similar to wall-associated receptor kinases, which were both associated with cell wall thickness. The results demonstrate feasibility of GWAS to identify novel candidate genes controlling industrially-relevant tracheid traits in Norway spruce.

Effect of number of annual rings and tree ages on genomic predictive ability for solid wood properties of Norway spruce.

BACKGROUND: Genomic selection (GS) or genomic prediction is considered as a promising approach to accelerate tree breeding and increase genetic gain by shortening breeding cycle, but the efforts to develop routines for operational breeding are so far limited. We investigated the predictive ability (PA) of GS based on 484 progeny trees from 62 half-sib families in Norway spruce (Picea abies (L.) Karst.) for wood density, modulus of elasticity (MOE) and microfibril angle (MFA) measured with SilviScan, as well as for measurements on standing trees by Pilodyn and Hitman instruments. RESULTS: GS predictive abilities were comparable with those based on pedigree-based prediction. Marker-based PAs were generally 25-30% higher for traits density, MFA and MOE measured with SilviScan than for their respective standing tree-based method which measured with Pilodyn and Hitman. Prediction accuracy (PC) of the standing tree-based methods were similar or even higher than increment core-based method. 78-95% of the maximal PAs of density, MFA and MOE obtained from coring to the pith at high age were reached by using data possible to obtain by drilling 3-5 rings towards the pith at tree age 10-12. CONCLUSIONS: This study indicates standing tree-based measurements is a cost-effective alternative method for GS. PA of GS methods were comparable with those pedigree-based prediction. The highest PAs were reached with at least 80-90% of the dataset used as training set. Selection for trait density could be conducted at an earlier age than for MFA and MOE. Operational breeding can also be optimized by training the model at an earlier age or using 3 to 5 outermost rings at tree age 10 to 12 years, thereby shortening the cycle and reducing the impact on the tree.

Association genetics identifies a specifically regulated Norway spruce laccase gene, PaLAC5, linked to Heterobasidion parviporum resistance.

It is important to improve the understanding of the interactions between the trees and pathogens and integrate this knowledge about disease resistance into tree breeding programs. The conifer Norway spruce (Picea abies) is an important species for the forest industry in Europe. Its major pathogen is Heterobasidion parviporum, causing stem and root rot. In this study, we identified 11 Norway spruce QTLs (Quantitative trait loci) that correlate with variation in resistance to H. parviporum in a population of 466 trees by association genetics. Individual QTLs explained between 2.1 and 5.2% of the phenotypic variance. The expression of candidate genes associated with the QTLs was analysed in silico and in response to H. parviporum hypothesizing that (a) candidate genes linked to control of fungal sapwood growth are more commonly expressed in sapwood, and; (b) candidate genes associated with induced defences are respond to H. parviporum inoculation. The Norway spruce laccase PaLAC5 associated with control of lesion length development is likely to be involved in the induced defences. Expression analyses showed that PaLAC5 responds specifically and strongly in close proximity to the H. parviporum inoculation. Thus, PaLAC5 may be associated with the lignosuberized boundary zone formation in bark adjacent to the inoculation site.

Transcriptome analysis of shade avoidance and shade tolerance in conifers.

MAIN CONCLUSION: Gymnosperms respond differently to light intensity and R:FR; although some aspects of shade response appear conserved, yet underlying mechanisms seem to be diverse in gymnosperms as compared to angiosperms. Shade avoidance syndrome (SAS) is well-characterized in the shade intolerant model species Arabidopsis thaliana whereas much less is known about shade tolerance response (STR), yet regulation of SAS and STR with reference to conifers remains poorly understood. We conducted a comparative study of two conifer species with contrasting responses to shade, Scots pine (shade-intolerant) and Norway spruce (shade-tolerant), with the aim to understand mechanisms behind SAS and STR in conifers. Pine and spruce seedlings were grown under controlled light and shade conditions, and hypocotyl and seedling elongation following different light treatments were determined in both species as indicators of shade responses. Red to far-red light ratio (R:FR) was shown to trigger the shade response in Norway spruce. In Scots pine, we observed an interaction between R:FR and light intensity. RNA sequencing (RNA-Seq) data revealed that SAS and STR responses included changes in expression of genes involved primarily in hormone signalling and pigment biosynthesis. From the RNA-Seq analysis, we propose that although some aspects of shade response appear to be conserved in angiosperms and gymnosperms, yet the underlying mechanisms may be different in gymnosperms that warrants further research.

Analysis of phenotypic- and Estimated Breeding Values (EBV) to dissect the genetic architecture of complex traits in a Scots pine three-generation pedigree design.

In forest tree breeding, family-based Quantitative Trait Loci (QTL) studies are valuable as methods to dissect the complexity of a trait and as a source of candidate genes. In the field of conifer research, our study contributes to the evaluation of phenotypic and predicted breeding values for the identification of QTL linked to complex traits in a three-generation pedigree population in Scots pine (Pinus sylvestris L.). A total of 11 470 open pollinated F2-progeny trees established at three different locations, were measured for growth and adaptive traits. Breeding values were predicted for their 360 mothers, originating from a single cross of two grand-parents. A multilevel LASSO association analysis was conducted to detect QTL using genotypes of the mothers with the corresponding phenotypes and Estimated Breeding Values (EBV). Different levels of genotype-by-environment (G × E) effects among sites at different years, were detected for survival and height. Moderate-to-low narrow sense heritabilities and EBV accuracies were found for all traits and all sites. We identified 18 AFLPs and 12 SNPs to be associated with QTL for one or more traits. 62 QTL were significant with percentages of variance explained ranging from 1.7 to 18.9%. In those cases where the same marker was associated to a phenotypic or an ebvQTL, the ebvQTL always explained higher proportion of the variance, maybe due to the more accurate nature of Estimated Breeding Values (EBV). Two SNP-QTL showed pleiotropic effects for traits related with hardiness, seed, cone and flower production. Furthermore, we detected several QTL with significant effects across multiple ages, which could be considered as strong candidate loci for early selection. The lack of reproducibility of some QTL detected across sites may be due to environmental heterogeneity reflected by the genotype- and QTL-by-environment effects.

Accuracy of genomic selection for growth and wood quality traits in two control-pollinated progeny trials using exome capture as the genotyping platform in Norway spruce.

BACKGROUND: Genomic selection (GS) can increase genetic gain by reducing the length of breeding cycle in forest trees. Here we genotyped 1370 control-pollinated progeny trees from 128 full-sib families in Norway spruce (Picea abies (L.) Karst.), using exome capture as genotyping platform. We used 116,765 high-quality SNPs to develop genomic prediction models for tree height and wood quality traits. We assessed the impact of different genomic prediction methods, genotype-by-environment interaction (G × E), genetic composition, size of the training and validation set, relatedness, and number of SNPs on accuracy and predictive ability (PA) of GS. RESULTS: Using G matrix slightly altered heritability estimates relative to pedigree-based method. GS accuracies were about 11-14% lower than those based on pedigree-based selection. The efficiency of GS per year varied from 1.71 to 1.78, compared to that of the pedigree-based model if breeding cycle length was halved using GS. Height GS accuracy decreased to more than 30% while using one site as training for GS prediction and using this model to predict the second site, indicating that G × E for tree height should be accommodated in model fitting. Using a half-sib family structure instead of full-sib structure led to a significant reduction in GS accuracy and PA. The full-sib family structure needed only 750 markers to reach similar accuracy and PA, as compared to 100,000 markers required for the half-sib family, indicating that maintaining the high relatedness in the model improves accuracy and PA. Using 4000-8000 markers in full-sib family structure was sufficient to obtain GS model accuracy and PA for tree height and wood quality traits, almost equivalent to that obtained with all markers. CONCLUSIONS: The study indicates that GS would be efficient in reducing generation time of breeding cycle in conifer tree breeding program that requires long-term progeny testing. The sufficient number of trees within-family (16 for growth and 12 for wood quality traits) and number of SNPs (8000) are required for GS with full-sib family relationship. GS methods had little impact on GS efficiency for growth and wood quality traits. GS model should incorporate G × E effect when a strong G × E is detected.

Non-functional plastid ndh gene fragments are present in the nuclear genome of Norway spruce (Picea abies L. Karsch): insights from in silico analysis of nuclear and organellar genomes.

Many genes have been lost from the prokaryote plastidial genome during the early events of endosymbiosis in eukaryotes. Some of them were definitively lost, but others were relocated and functionally integrated to the host nuclear genomes through serial events of gene transfer during plant evolution. In gymnosperms, plastid genome sequencing has revealed the loss of ndh genes from several species of Gnetales and Pinaceae, including Norway spruce (Picea abies). This study aims to trace the ndh genes in the nuclear and organellar Norway spruce genomes. The plastid genomes of higher plants contain 11 ndh genes which are homologues of mitochondrial genes encoding subunits of the proton-pumping NADH-dehydrogenase (nicotinamide adenine dinucleotide dehydrogenase) or complex I (electron transport chain). Ndh genes encode 11 NDH polypeptides forming the Ndh complex (analogous to complex I) which seems to be primarily involved in chloro-respiration processes. We considered ndh genes from the plastidial genome of four gymnosperms (Cryptomeria japonica, Cycas revoluta, Ginkgo biloba, Podocarpus totara) and a single angiosperm species (Arabidopsis thaliana) to trace putative homologs in the nuclear and organellar Norway spruce genomes using tBLASTn to assess the evolutionary fate of ndh genes in Norway spruce and to address their genomic location(s), structure, integrity and functionality. The results obtained from tBLASTn were subsequently analyzed by performing homology search for finding ndh specific conserved domains using conserved domain search. We report the presence of non-functional plastid ndh gene fragments, excepting ndhE and ndhG genes, in the nuclear genome of Norway spruce. Regulatory transcriptional elements like promoters, TATA boxes and enhancers were detected in the upstream regions of some ndh fragments. We also found transposable elements in the flanking regions of few ndh fragments suggesting nuclear rearrangements in those regions. These evidences support the hypothesis that, at least in Picea, ndh translocations from the plastid to the nuclear genome have occurred, and that there might have been a functional machinery at some time during evolution to accommodate them within a nuclear-encoded environment, or attempts to form it.

Comparative in silico analysis of SSRs in coding regions of high confidence predicted genes in Norway spruce (Picea abies) and Loblolly pine (Pinus taeda).

BACKGROUND: Microsatellites or simple sequence repeats (SSRs) are DNA sequences consisting of 1-6 bp tandem repeat motifs present in the genome. SSRs are considered to be one of the most powerful tools in genetic studies. We carried out a comparative study of perfect SSR loci belonging to class I (≥20) and class II (≥12 and <20 bp) types located in coding regions of high confidence genes in Picea abies and Pinus taeda. SSRLocator was used to retrieve SSRs from the full length CDS of predicted genes in both species. RESULTS: Trimers were the most abundant motifs in class I followed by hexamers in Picea abies, while trimers and hexamers were equally abundant in Pinus taeda class I SSRs. Hexamers were most frequent within class II SSRs followed by trimers, in both species. Although the frequency of genes containing SSRs was slightly higher in Pinus taeda, SSR counts per Mbp for class I was similar in both species (P-value = 0.22); while for class II SSRs, it was significantly higher in Picea abies (P-value = 0.00009). AT-rich motifs were higher in abundance than the GC-rich motifs, within class II SSRs in both the species (P-values = 10(-9) and 0). With reference to class I SSRs, AT-rich and GC-rich motifs were detected with equal frequency in Pinus taeda (P-value = 0.24); while in Picea abies, GC-rich motifs were detected with higher frequency than the AT-rich motifs (P-value = 0.0005). CONCLUSIONS: Our study gives a comparative overview of the genome SSRs composition based on high confidence genes in the two recently sequenced and economically important conifers and, also provides information on functional molecular markers that can be applied in genetic studies in Pinus and Picea species.

Evolutionary aspects of functional and pseudogene members of the phytochrome gene family in Scots pine.

According to the neutral theory of evolution, mutation and genetic drift are the only forces that shape unconstrained, neutral, gene evolution. Thus, pseudogenes (which often evolve neutrally) provide opportunities to obtain direct estimates of mutation rates that are not biased by selection, and gene families comprising functional and pseudogene members provide useful material for both estimating neutral mutation rates and identifying sites that appear to be under positive or negative selection pressures. Conifers could be very useful for such analyses since they have large and complex genomes. There is evidence that pseudogenes make significant contributions to the size and complexity of gene families in pines, although few studies have examined the composition and evolution of gene families in conifers. In this work, I examine the complexity and rates of mutation of the phytochrome gene family in Pinus sylvestris and show that it includes not only functional genes but also pseudogenes. As expected, the functional PHYO does not appear to have evolved neutrally, while phytochrome pseudogenes show signs of unconstrained evolution.