Integrated biological and chemical characterisation of a pair of leonardesque canal lock gates.

The Museo Nazionale della Scienza e della Tecnologia "Leonardo da Vinci" in Milan is exposing two pairs of canal lock gates, used to control the water flow in Milan canal system, whose design appears in the Leonardo's Codex Atlanticus. The wood present in the gates has been deeply characterised by mean of a multidisciplinary investigation involving i) DNA barcoding of wood fragments; ii) microbial community characterisation, and iii) chemical analyses. DNA barcoding revealed that two fragments of the gates belonged to wood species widely used in the middle age: Fagus sylvatica and Picea abies. The chemical characterisations were based on the use of ionic liquid as dissolving medium in order to analyse the entire cell wall material by means of Gel Permeation Chromatography (GPC) and 2D-NMR-HSQC techniques. This multidisciplinary analytical approach was able to highlight the complex nature of the degradation occurred during the gate operation (XVI-XVIII centuries): an intricate interplay between microbial populations (i.e. Shewanella), inorganic factors (i.e. iron from nails), physical factors and the lignocellulosic material.

Different Roles of Auxins in Somatic Embryogenesis Efficiency in Two Picea Species.

The effects of auxins 2,4-D (2,4-dichlorophenoxyacetic acid), NAA (1-naphthaleneacetic acid) or picloram (4-amino-3,5,6-trichloropicolinic acid; 9 µM) and cytokinin BA (benzyloadenine; 4.5 µM) applied in the early stages of somatic embryogenesis (SE) on specific stages of SE in Picea abies and P. omorika were investigated. The highest SE initiation frequency was obtained after 2,4-D application in P. omorika (22.00%) and picloram application in P. abies (10.48%). NAA treatment significantly promoted embryogenic tissue (ET) proliferation in P. abies, while 2,4-D treatment reduced it. This reduction was related to the oxidative stress level, which was lower with the presence of NAA in the proliferation medium and higher with the presence of 2,4-D. The reduced oxidative stress level after NAA treatment suggests that hydrogen peroxide (H2O2) acts as a signalling molecule and promotes ET proliferation. NAA and picloram in the proliferation medium decreased the further production and maturation of P. omorika somatic embryos compared with that under 2,4-D. The quality of the germinated P. abies embryos and their development into plantlets depended on the auxin type and were the highest in NAA-originated embryos. These results show that different auxin types can generate different physiological responses in plant materials during SE in both spruce species.

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.

Phylogenomics disentangles the evolutionary history of spruces (Picea) in the Qinghai-Tibetan Plateau: Implications for the design of population genetic studies and species delimitation of conifers.

A laborious and difficult task in current tree of life reconstruction is to resolve evolutionary relationships of closely related congeneric species that originated from recent radiations. This is particularly difficult for forest species with long generation times and large effective population sizes such as conifers. The Qinghai-Tibetan Plateau (QTP) and adjacent areas are considered a species diversity center of Picea, harboring 11 species (including 5 varieties) of this genus, but evolutionary relationships of these species are far from being resolved due to recent radiations, morphological convergence, and frequent interspecific gene flow. In this study, we use these spruce species to test whether phylotranscriptomic analysis, combined with population genetic analysis, can disentangle their evolutionary relationships, and to explore whether reticulate evolution has occurred among them. Phylogenomic analyses indicate that all spruce species in the QTP and neighboring areas, except P. asperata and P. crassifolia, cluster together, and in particular, nearly all taxa (including varieties) reflect reciprocally monophyletic lineages, although the two species P. likiangensis and P. brachytyla are not monophyletic. We found that, compared to herbaceous plants, many more genes (a minimum of 600 OGs for Picea) are required to resolve interspecific relationships of conifers. Contrary to previous studies, our data do not support a hybrid origin of P. purpurea, but suggests a hybrid origin for P. brachytyla var. brachytyla and P. likiangensis var. rubescens. We emphasize that the species or species complex used for population genetic and phylogeographical studies should be monophyletic.

Phylotranscriptomics resolves interspecific relationships and indicates multiple historical out-of-North America dispersals through the Bering Land Bridge for the genus Picea (Pinaceae).

A robust phylogeny is prerequisite to understand the evolution and biogeography of organisms. However, ancient and recent evolutionary radiations occurred in many plant lineages, which pose great challenges for phylogenetic analysis, especially for conifers characterized by large effective population sizes and long generation times. Picea is an important component of the dark coniferous forests in the Northern Hemisphere. Previous studies improved our understanding of its evolutionary history, but its interspecific relationships and biogeographic history remain largely unresolved. In the present study, we reconstructed a well-resolved phylogeny of Picea by comparative transcriptomic analysis based on a complete species sampling. The phylogenetic analysis, together with molecular dating and ancestral area reconstruction, further supports the North American origin hypothesis for Picea, and indicates that this genus experienced multiple out-of-North America dispersals by the Bering Land Bridge. We also found that spruces in the Japanese Archipelago have multiple origins, and P. morrisonicola from the Taiwan Island has a close relationship with species from the Qinghai-Tibetan Plateau and adjacent regions. Our study provides the first complete phylogeny of Picea at the genomic level, which is important for future studies of this genus.

PICEAdatabase: a web database for Picea omics and phenotypic information.

Picea belongs to the Pinaceae family and is a famous commercial tree species because of its straight trunk and excellent timber traits. Recently, omics have been widely used for fundamental and mechanism studies on Picea plants. To improve the accessibility to omics and phenotypic data and facilitate further studies, we compiled the sequences of 2 chloroplast genomes (Picea crassifolia and Picea asperata) and 32 complete omics data sets, including 20 transcriptomes, 4 proteomes, 2 degradomes and 6 microRNAs from P. crassifolia, P. asperata, Picea balfouriana and Picea abies tissues under different treatments, in PICEAdatabase. In addition, phenotypic data on plant growth and wood property traits were collected from two field trials of P. crassifolia. PICEAdatabase also includes useful analysis tools, such as BLAST, DESeq2 and JBrowse, to assist with analyses.

Genome-wide association study identified novel candidate loci affecting wood formation in Norway spruce.

Norway spruce is a boreal forest tree species of significant ecological and economic importance. Hence there is a strong imperative to dissect the genetics underlying important wood quality traits in the species. We performed a functional genome-wide association study (GWAS) of 17 wood traits in Norway spruce using 178 101 single nucleotide polymorphisms (SNPs) generated from exome genotyping of 517 mother trees. The wood traits were defined using functional modelling of wood properties across annual growth rings. We applied a Least Absolute Shrinkage and Selection Operator (LASSO-based) association mapping method using a functional multilocus mapping approach that utilizes latent traits, with a stability selection probability method as the hypothesis testing approach to determine a significant quantitative trait locus. The analysis provided 52 significant SNPs from 39 candidate genes, including genes previously implicated in wood formation and tree growth in spruce and other species. Our study represents a multilocus GWAS for complex wood traits in Norway spruce. The results advance our understanding of the genetics influencing wood traits and identifies candidate genes for future functional studies.

Genetic evidence for species cohesion, substructure and hybrids in spruce.

The origin and history of species are shaped by various evolutionary dynamics, including their persistence in the face of potential gene flow from related taxa. In this study, we use broad geographical and taxonomic sampling (2,219 individuals) to establish the distribution of species, hybrids and cryptic genetic variation within the conifer genus Picea (spruce) across western North America. We demonstrate that the six species of spruce in this region are distinguishable based on their genetic composition, and that the more closely related Engelmann spruce (P. engelmannii) and white spruce (P. glauca) have generated numerous and widespread hybrids. These hybrids occur in the central Rocky Mountains, well to the south of the well-established region of admixture in Canada. Additionally, we provide evidence for subdivision within Engelmann spruce, manifested as a southern Rocky Mountains form, and a northern Rocky Mountain and Cascade mountains (western) form. In the intervening central Rocky Mountains region (forests in Wyoming and adjacent states) we found primarily individuals with admixed ancestry. Following their origin, these species of spruce have interacted repeatedly and in different geographical contexts. Multiple pairs of species have been shown to hybridize, yet the species persist and retain distinguishable compositions. At the same time, large geographical areas exist where hybrids are pervasive. Consequently, spruce provide a case study for the maintenance of species boundaries, particularly for how widespread hybridization need not lead to the collapse and loss of species.

Population genomic analysis reveals that homoploid hybrid speciation can be a lengthy process.

An increasing number of species are thought to have originated by homoploid hybrid speciation (HHS), but in only a handful of cases are details of the process known. A previous study indicated that Picea purpurea, a conifer in the Qinghai-Tibet Plateau (QTP), originated through HHS from P. likiangensis and P. wilsonii. To investigate this origin in more detail, we analysed transcriptome data for 114 individuals collected from 34 populations of the three Picea species from their core distributions in the QTP. Phylogenetic, principal component and admixture analyses of nuclear SNPs showed the species to be delimited genetically and that P. purpurea was admixed with approximately 60% of its ancestry derived from P. wilsonii and 40% from P. likiangensis. Coalescent simulations revealed the best-fitting model of origin involved formation of an intermediate hybrid lineage between P. likiangensis and P. wilsonii approximately 6 million years ago (mya), which backcrossed to P. wilsonii to form P. purpurea approximately one mya. The intermediate hybrid lineage no longer exists and is referred to as a "ghost" lineage. Our study emphasizes the power of population genomic analysis combined with coalescent analysis for reconstructing the stages involved in the origin of a homoploid hybrid species over an extended period. In contrast to other studies, we show that these stages can in some instances span a relatively long period of evolutionary time.

Phylogeny and evolutionary history of Pinaceae updated by transcriptomic analysis.

Pinaceae comprises 11 genera, and represents the largest family of conifers with an extensive wild distribution in the Northern Hemisphere. Intergeneric relationships of Pinaceae have been investigated using many morphological characters and molecular markers, but phylogenetic positions of four genera, including Cathaya, Cedrus, Nothotsuga and Pseudolarix, remain controversial or have not been completely resolved. To completely resolve the intergeneric relationships of Pinaceae, we conducted a comparative transcriptomic study of 14 species representing all Pinaceae genera. Multiple data sets, containing up to 6,369,681 sites across 4676 loci, were analyzed using concatenation and coalescent methods. Our study generated a robust topology, which divides Pinaceae into two clades, one (pinoid) including Cathaya, Larix, Picea, Pinus, and Pseudotsuga, and the other (abietoid) including Abies, Cedrus, Keteleeria, Nothotsuga, Pseudolarix, and Tsuga. Cathaya and Pinus form a clade sister to Picea; Cedrus is sister to the remaining abietoid genera, and the two genera Nothotsuga and Tsuga form a clade sister to Pseudolarix. The discordant positions of Cathaya, Cedrus and Pseudolarix in different gene trees could be explained by ancient radiation and/or molecular homoplastic evolution. The hybrid origin hypothesis of Nothotsuga is not supported. Based on molecular dating, extant Pinaceae genera diverged since about 206 Mya, earlier than the break-up of Pangea, and the divergence among the pinoid genera occurred earlier than the split among the abietoid genera. Moreover, our study indicates that two radiation events occurred in the evolution of Pinaceae genera, and some important morphological characters evolved multiple times based on ancestral state reconstruction.

Timber isoscapes. A case study in a mountain area in the Italian Alps.

BACKGROUND: Local timber is still one of the main sources of work and income for mountain communities. However, illegal logging is a major cause of deforestation in many countries and has significant impacts on local communities and biodiversity. Techniques for tracing timber would provide a useful tool to protect local timber industries and contribute to the fight against illegal logging. Although considerable progress has been made in food traceability, timber provenance is still a somewhat neglected research area. Stable isotope ratios in plants are known to reflect geographical variations. This study reports accurate spatial distribution of δ18O and δ2H in timber from north-eastern Italy (Trentino) in order to trace geographical origin. METHODOLOGY AND PRINCIPAL FINDINGS: We tested the accuracy of four kriging methods using an annual resolution of δ18O and δ2H measured in Picea abies. Pearson's correlation coefficients revealed altitude to be the most appropriate covariate for the cokriging model, which has ultimately proved to be the best method due to its low estimation error. CONCLUSIONS: We present regional maps of interpolated δ18O and δ2H in Picea abies wood together with the 95% confidence intervals. The strong spatial structure of the data demonstrates the potential of multivariate spatial interpolation, even in a highly heterogeneous area such as the Alps. We believe that this geospatial approach can be successfully applied on a wider scale in order to combat illegal logging.

Overexpression of PaNAC03, a stress induced NAC gene family transcription factor in Norway spruce leads to reduced flavonol biosynthesis and aberrant embryo development.

BACKGROUND: The NAC family of transcription factors is one of the largest gene families of transcription factors in plants and the conifer NAC gene family is at least as large, or possibly larger, as in Arabidopsis. These transcription factors control both developmental and stress induced processes in plants. Yet, conifer NACs controlling stress induced processes has received relatively little attention. This study investigates NAC family transcription factors involved in the responses to the pathogen Heterobasidion annosum (Fr.) Bref. sensu lato. RESULTS: The phylogeny and domain structure in the NAC proteins can be used to organize functional specificities, several well characterized stress-related NAC proteins are found in III-3 in Arabidopsis (Jensen et al. Biochem J 426:183-196, 2010). The Norway spruce genome contain seven genes with similarity to subgroup III-3 NACs. Based on the expression pattern PaNAC03 was selected for detailed analyses. Norway spruce lines overexpressing PaNAC03 exhibited aberrant embryo development in response to maturation initiation and 482 misregulated genes were identified in proliferating cultures. Three key genes in the flavonoid biosynthesis pathway: a CHS, a F3'H and PaLAR3 were consistently down regulated in the overexpression lines. In accordance, the overexpression lines showed reduced levels of specific flavonoids, suggesting that PaNAC03 act as a repressor of this pathway, possibly by directly interacting with the promoter of the repressed genes. However, transactivation studies of PaNAC03 and PaLAR3 in Nicotiana benthamiana showed that PaNAC03 activated PaLAR3A, suggesting that PaNAC03 does not act as an independent negative regulator of flavan-3-ol production through direct interaction with the target flavonoid biosynthetic genes. CONCLUSIONS: PaNAC03 and its orthologs form a sister group to well characterized stress-related angiosperm NAC genes and at least PaNAC03 is responsive to biotic stress and appear to act in the control of defence associated secondary metabolite production.

Assembly of the Complete Sitka Spruce Chloroplast Genome Using 10X Genomics' GemCode Sequencing Data.

The linked read sequencing library preparation platform by 10X Genomics produces barcoded sequencing libraries, which are subsequently sequenced using the Illumina short read sequencing technology. In this new approach, long fragments of DNA are partitioned into separate micro-reactions, where the same index sequence is incorporated into each of the sequencing fragment inserts derived from a given long fragment. In this study, we exploited this property by using reads from index sequences associated with a large number of reads, to assemble the chloroplast genome of the Sitka spruce tree (Picea sitchensis). Here we report on the first Sitka spruce chloroplast genome assembled exclusively from P. sitchensis genomic libraries prepared using the 10X Genomics protocol. We show that the resulting 124,049 base pair long genome shares high sequence similarity with the related white spruce and Norway spruce chloroplast genomes, but diverges substantially from a previously published P. sitchensis- P. thunbergii chimeric genome. The use of reads from high-frequency indices enabled separation of the nuclear genome reads from that of the chloroplast, which resulted in the simplification of the de Bruijn graphs used at the various stages of assembly.

Identifying Genetic Signatures of Natural Selection Using Pooled Population Sequencing in Picea abies.

The joint inference of selection and past demography remain a costly and demanding task. We used next generation sequencing of two pools of 48 Norway spruce mother trees, one corresponding to the Fennoscandian domain, and the other to the Alpine domain, to assess nucleotide polymorphism at 88 nuclear genes. These genes are candidate genes for phenological traits, and most belong to the photoperiod pathway. Estimates of population genetic summary statistics from the pooled data are similar to previous estimates, suggesting that pooled sequencing is reliable. The nonsynonymous SNPs tended to have both lower frequency differences and lower FST values between the two domains than silent ones. These results suggest the presence of purifying selection. The divergence between the two domains based on synonymous changes was around 5 million yr, a time similar to a recent phylogenetic estimate of 6 million yr, but much larger than earlier estimates based on isozymes. Two approaches, one of them novel and that considers both FST and difference in allele frequencies between the two domains, were used to identify SNPs potentially under diversifying selection. SNPs from around 20 genes were detected, including genes previously identified as main target for selection, such as PaPRR3 and PaGI.

The extent and meaning of hybridization and introgression between Siberian spruce (Picea obovata) and Norway spruce (Picea abies): cryptic refugia as stepping stones to the west?

Boreal species were repeatedly exposed to ice ages and went through cycles of contraction and expansion while sister species alternated periods of contact and isolation. The resulting genetic structure is consequently complex, and demographic inferences are intrinsically challenging. The range of Norway spruce (Picea abies) and Siberian spruce (Picea obovata) covers most of northern Eurasia; yet their geographical limits and histories remain poorly understood. To delineate the hybrid zone between the two species and reconstruct their joint demographic history, we analysed variation at nuclear SSR and mitochondrial DNA in 102 and 88 populations, respectively. The dynamics of the hybrid zone was analysed with approximate Bayesian computation (ABC) followed by posterior predictive structure plot reconstruction and the presence of barriers across the range tested with estimated effective migration surfaces. To estimate the divergence time between the two species, nuclear sequences from two well-separated populations of each species were analysed with ABC. Two main barriers divide the range of the two species: one corresponds to the hybrid zone between them, and the other separates the southern and northern domains of Norway spruce. The hybrid zone is centred on the Urals, but the genetic impact of Siberian spruce extends further west. The joint distribution of mitochondrial and nuclear variation indicates an introgression of mitochondrial DNA from Norway spruce into Siberian spruce. Overall, our data reveal a demographic history where the two species interacted frequently and where migrants originating from the Urals and the West Siberian Plain recolonized northern Russia and Scandinavia using scattered refugial populations of Norway spruce as stepping stones towards the west.

Implementation of the Realized Genomic Relationship Matrix to Open-Pollinated White Spruce Family Testing for Disentangling Additive from Nonadditive Genetic Effects.

The open-pollinated (OP) family testing combines the simplest known progeny evaluation and quantitative genetics analyses as candidates' offspring are assumed to represent independent half-sib families. The accuracy of genetic parameter estimates is often questioned as the assumption of "half-sibling" in OP families may often be violated. We compared the pedigree- vs. marker-based genetic models by analysing 22-yr height and 30-yr wood density for 214 white spruce [Picea glauca (Moench) Voss] OP families represented by 1694 individuals growing on one site in Quebec, Canada. Assuming half-sibling, the pedigree-based model was limited to estimating the additive genetic variances which, in turn, were grossly overestimated as they were confounded by very minor dominance and major additive-by-additive epistatic genetic variances. In contrast, the implemented genomic pairwise realized relationship models allowed the disentanglement of additive from all nonadditive factors through genetic variance decomposition. The marker-based models produced more realistic narrow-sense heritability estimates and, for the first time, allowed estimating the dominance and epistatic genetic variances from OP testing. In addition, the genomic models showed better prediction accuracies compared to pedigree models and were able to predict individual breeding values for new individuals from untested families, which was not possible using the pedigree-based model. Clearly, the use of marker-based relationship approach is effective in estimating the quantitative genetic parameters of complex traits even under simple and shallow pedigree structure.

Development of highly reliable in silico SNP resource and genotyping assay from exome capture and sequencing: an example from black spruce (Picea mariana).

Picea mariana is a widely distributed boreal conifer across Canada and the subject of advanced breeding programmes for which population genomics and genomic selection approaches are being developed. Targeted sequencing was achieved after capturing P. mariana exome with probes designed from the sequenced transcriptome of Picea glauca, a distant relative. A high capture efficiency of 75.9% was reached although spruce has a complex and large genome including gene sequences interspersed by some long introns. The results confirmed the relevance of using probes from congeneric species to perform successfully interspecific exome capture in the genus Picea. A bioinformatics pipeline was developed including stringent criteria that helped detect a set of 97,075 highly reliable in silico SNPs. These SNPs were distributed across 14,909 genes. Part of an Infinium iSelect array was used to estimate the rate of true positives by validating 4267 of the predicted in silico SNPs by genotyping trees from P. mariana populations. The true positive rate was 96.2% for in silico SNPs, compared to a genotyping success rate of 96.7% for a set 1115 P. mariana control SNPs recycled from previous genotyping arrays. These results indicate the high success rate of the genotyping array and the relevance of the selection criteria used to delineate the new P. mariana in silico SNP resource. Furthermore, in silico SNPs were generally of medium to high frequency in natural populations, thus providing high informative value for future population genomics applications.

Heterotrimeric G-proteins in Picea abies and their regulation in response to Heterobasidion annosum s.l. infection.

BACKGROUND: Heterotrimeric G-proteins are important signalling switches, present in all eukaryotic kingdoms. In plants they regulate several developmental functions and play an important role in plant-microbe interactions. The current knowledge on plant G-proteins is mostly based on model angiosperms and little is known about the G-protein repertoire and function in other lineages. In this study we investigate the heterotrimeric G-protein subunit repertoire in Pinaceae, including phylogenetic relationships, radiation and sequence diversity levels in relation to other plant linages. We also investigate functional diversification of the G-protein complex in Picea abies by analysing transcriptional regulation of the G-protein subunits in different tissues and in response to pathogen infection. RESULTS: A full repertoire of G-protein subunits in several conifer species were identified in silico. The full-length P. abies coding regions of one Gα-, one Gß- and four Gγ-subunits were cloned and sequenced. The phylogenetic analysis of the Gγ-subunits showed that PaGG1 clustered with A-type-like subunits, PaGG3 and PaGG4 clustered with C-type-like subunits, while PaGG2 and its orthologs represented a novel conifer-specific putative Gγ-subunit type. Gene expression analyses by quantitative PCR of P. abies G-protein subunits showed specific up-regulation of the Gα-subunit gene PaGPA1 and the Gγ-subunit gene PaGG1 in response to Heterobasidion annosum sensu lato infection. CONCLUSIONS: Conifers possess a full repertoire of G-protein subunits. The differential regulation of PaGPA1 and PaGG1 indicates that the heterotrimeric G-protein complex represents a critical linchpin in Heterobasidion annosum s.l. perception and downstream signaling in P. abies.

Legacy of Pre-Disturbance Spatial Pattern Determines Early Structural Diversity following Severe Disturbance in Montane Spruce Forests.

BACKGROUND: Severe canopy-removing disturbances are native to many temperate forests and radically alter stand structure, but biotic legacies (surviving elements or patterns) can lend continuity to ecosystem function after such events. Poorly understood is the degree to which the structural complexity of an old-growth forest carries over to the next stand. We asked how pre-disturbance spatial pattern acts as a legacy to influence post-disturbance stand structure, and how this legacy influences the structural diversity within the early-seral stand. METHODS: Two stem-mapped one-hectare forest plots in the Czech Republic experienced a severe bark beetle outbreak, thus providing before-and-after data on spatial patterns in live and dead trees, crown projections, down logs, and herb cover. RESULTS: Post-disturbance stands were dominated by an advanced regeneration layer present before the disturbance. Both major species, Norway spruce (Picea abies) and rowan (Sorbus aucuparia), were strongly self-aggregated and also clustered to former canopy trees, pre-disturbance snags, stumps and logs, suggesting positive overstory to understory neighbourhood effects. Thus, although the disturbance dramatically reduced the stand's height profile with ~100% mortality of the canopy layer, the spatial structure of post-disturbance stands still closely reflected the pre-disturbance structure. The former upper tree layer influenced advanced regeneration through microsite and light limitation. Under formerly dense canopies, regeneration density was high but relatively homogeneous in height; while in former small gaps with greater herb cover, regeneration density was lower but with greater heterogeneity in heights. CONCLUSION: These findings suggest that pre-disturbance spatial patterns of forests can persist through severe canopy-removing disturbance, and determine the spatial structure of the succeeding stand. Such patterns constitute a subtle but key legacy effect, promoting structural complexity in early-seral forests as well as variable successional pathways and rates. This influence suggests a continuity in spatial ecosystem structure that may well persist through multiple forest generations.

Tracking the progression of speciation: variable patterns of introgression across the genome provide insights on the species delimitation between progenitor-derivative spruces (Picea mariana × P. rubens).

The genic species concept implies that while most of the genome can be exchanged somewhat freely between species through introgression, some genomic regions remain impermeable to interspecific gene flow. Hence, interspecific differences can be maintained despite ongoing gene exchange within contact zones. This study assessed the heterogeneous patterns of introgression at gene loci across the hybrid zone of an incipient progenitor-derivative species pair, Picea mariana (black spruce) and Picea rubens (red spruce). The spruce taxa likely diverged in geographic isolation during the Pleistocene and came into secondary contact during late Holocene. A total of 300 SNPs distributed across the 12 linkage groups (LG) of black spruce were genotyped for 385 individual trees from 33 populations distributed across the allopatric zone of each species and within the zone of sympatry. An integrative framework combining three population genomic approaches was used to scan the genomes, revealing heterogeneous patterns of introgression. A total of 23 SNPs scattered over 10 LG were considered impermeable to introgression and putatively under diverging selection. These loci revealed the existence of impermeable genomic regions forming the species boundary and are thus indicative of ongoing speciation between these two genetic lineages. Another 238 SNPs reflected selectively neutral diffusion across the porous species barrier. Finally, 39 highly permeable SNPs suggested ancestral polymorphism along with balancing selection. The heterogeneous patterns of introgression across the genome indicated that the speciation process between black spruce and red spruce is young and incomplete, albeit some interspecific differences are maintained, allowing ongoing species divergence even in sympatry. The approach developed in this study can be used to track the progression of ongoing speciation processes.