Long-insert sequence capture detects high copy numbers in a defence-related beta-glucosidase gene βglu-1 with large variations in white spruce but not Norway spruce

Conifers are long-lived and slow-evolving, thus requiring effective defences against their fast-evolving insect natural enemies. The copy number variation (CNV) of two key acetophenone biosynthesis genes Ugt5/Ugt5b and βglu-1 may provide a plausible mechanism underlying the constitutively variable defence in white spruce (Picea glauca) against its primary defoliator, spruce budworm. This study develops a long-insert sequence capture probe set (Picea_hung_p1.0) for quantifying copy number of βglu-1-like, Ugt5-like genes and single-copy genes on 38 Norway spruce (Picea abies) and 40 P. glauca individuals from eight and nine provenances across Europe and North America respectively. We developed local assemblies (Piabi_c1.0 and Pigla_c.1.0), full-length transcriptomes (PIAB_v1 and PIGL_v1), and gene models to characterise the diversity of βglu-1 and Ugt5 genes. We observed very large copy numbers of βglu-1, with up to 381 copies in a single P. glauca individual. We observed among-provenance CNV of βglu-1 in P. glauca but not P. abies. Ugt5b was predominantly single-copy in both species. This study generates critical hypotheses for testing the emergence and mechanism of extreme CNV, the dosage effect on phenotype, and the varying copy number of genes with the same pathway. We demonstrate new approaches to overcome experimental challenges in genomic research in conifer defences.

Climatic responses and variability in bark anatomical traits of 23 Picea species

In woody plants, bark is an important protective tissue which can participate in photosynthesis, manage water loss, and transport assimilates. Studying the bark anatomical traits can provide insight into plant environmental adaptation strategies. However, a systematic understanding of the variability in bark anatomical traits and their drivers is lacking in woody plants. In this study, the bark anatomical traits of 23 Picea species were determined in a common garden experiment. We analyzed interspecific differences and interpreted the patterns in bark anatomical traits in relation to phylogenetic relationships and climatic factors of each species according to its global distribution. The results showed that there were interspecific differences in bark anatomical traits of Picea species. Phloem thickness was positively correlated with parenchyma cell size, possibly related to the roles of parenchyma cells in the radial transport of assimilates. Sieve cell size was negatively correlated with the radial diameter of resin ducts, and differences in sieve cells were possibly related to the formation and expansion of resin ducts. There were no significant phylogenetic signals for any bark anatomical trait, except the tangential diameter of resin ducts. Phloem thickness and parenchyma cell size were affected by temperature-related factors of their native range, while sieve cell size was influenced by precipitation-related factors. Bark anatomical traits were not significantly different under wet and dry climates. This study makes an important contribution to our understanding of variability in bark anatomical traits among Picea species and their ecological adaptations.

Morphological disparity of extant and extinct pinaceous ovulate cones: How many cones are enough?

PREMISE

The expansion of Pinaceae during the Cretaceous is exemplified by the numerous ovulate cone taxa found in western Europe and North America. The Belgian Wealden facies deposits have delivered hundreds of exceptionally well-preserved yet isolated pinaceous ovulate cones; these cones were placed by convention within form-genera. Ten species have been described in Belgium, representing about 20% of the known fossil record of this period. However, the validity of these taxa is questionable because their intra- and interspecific variabilities have never been thoroughly studied. Moreover, quantifying the expansion of Pinaceae in terms of morphospace occupation is desirable to reveal the dynamics of this critical radiation.

METHODS

We used linear and geometric morphometry to quantify the shape of the extensive sample of Cretaceous cones of Belgium. These methods were also applied to extant pinaceous species to compare the morphological disparity of Cretaceous assemblages against those of today in selected ecosystems. We used ordination methods (PCA) to visualize morphospace occupation and test for species delineation.

RESULTS

The morphological disparity was not higher in fossil species than in extant species we sampled. Both morphological approaches confirmed that the species Pityostrobus andraei presents high morphological variability. Our resampling tests indicate that ovulate cone morphological variability can be satisfactorily quantified with as few as 15 specimens.

CONCLUSIONS

The methodology used here is relevant for quantifying both the variability and the diversity of many fossil assemblages, paving the way for a more robust evaluation of Cretaceous pinaceous diversity.

Physicochemical properties and homology studies of the floral meristem identity gene LFY in nonflowering and flowering plants

In flowering plants, the LEAFY (LFY) gene controls floral meristem activity. In early land plants such as mosses and ferns, it, however, has a minimum role in cell division and development of diploid sporophyte. Homology modeling, an accurate and efficient protein structure prediction method, was used to construct a 3D model of the LEAFY protein in nonflowering and flowering plants. The present study examines the following species: Charophyte green algae, Physcomitrella, Ceratopteris, Picea, and Arabidopsis, as they are the popularly used model organisms for developmental studies. LEAFY protein sequences from the model organisms were aligned by multiple sequence alignment. 3D models of the LEAFY protein from all the model organisms was constructed using the PHYRE2 program with 100% confidence, and the constructed models were evaluated using the MolProbity tool. On the basis of the conserved regions, Charophyte green algae shared 38-46% sequence similarity with Physcomitrella sp., 37-46% similarity with Ceratopteris sp., 33-41% similarity with Picea sp., and 32-38% similarity with Arabidopsis sp. The Motif Finder server identified the protein family domain FLO_LFY and LFY_SAM, whose function is floral meristem development. Secondary structure prediction analysis indicated that the LEAFY protein belongs to the alpha (α) protein class, which is stable against mutation and thus limits structural changes in the LEAFY protein. The study findings reveal two distinct clusters of the LFY gene from the common ancestor green algae. One cluster is present in nonflowering plants that include mosses, pteridophytes, and gymnosperms, and the other cluster is present in flowering plants that include orchids, monocots, dicots, and angiosperms.

Pollen Germination and Pollen Tube Growth in Gymnosperms

Pollen germination and pollen tube growth are common to all seed plants, but these processes first developed in gymnosperms and still serve for their successful sexual reproduction. The main body of data on the reproductive physiology, however, was obtained on flowering plants, and one should be careful to extrapolate the discovered patterns to gymnosperms. In recent years, physiological studies of coniferous pollen have been increasing, and both the features of this group and the similarities with flowering plants have already been identified. The main part of the review is devoted to physiological studies carried out on conifer pollen. The main properties and diversity of pollen grains and pollination strategies in gymnosperms are described.

Probability of Spring Frosts, Not Growing Degree-Days, Drives Onset of Spruce Bud Burst in Plantations at the Boreal-Temperate Forest Ecotone

Climate warming-driven early leaf-out is expected to increase forest productivity but concurrently increases leaf exposure to spring frosts, which could reduce forests' net productivity. We hypothesized that due to their damaging effect on buds, spring frosts exert a stronger control on bud phenology than do growing degree-days. We monitored bud flush phenology of three white spruce seed sources (one local seed source from the boreal mixedwood forest and two seed sources from the temperate forest), one black spruce seed source originating from the boreal mixedwood forest and four nonlocal Norway spruce seed sources in 2016 and 2017 in two plantations located on both sides of the temperate-boreal mixedwood forest ecotone in eastern Canada (Quebec). We aimed to determine inter- and intraspecies variations in bud break timing and sensitivity to air temperature and photoperiod. We expected that bud break timing for boreal species and seed sources would be better synchronized with the decrease in frost probability than for nonlocal species and seed sources. We used mixed binomial regressions and AICc model selection to determine the best environmental variables predicting each transition from one stage of bud phenology to the next. At both plantation sites, white spruce bud flush began and ended earlier compared to black and Norway spruce. Buds of all spruce species were sensitive to frost probability for early phenological stages, whereas growing degree-days controlled the remaining stages. Photoperiod sensitivity was higher for white spruce compared to black and Norway spruce and reached its maximum in the temperate forest. At intraspecies level, the two southern white spruce seed sources opened their buds earlier than the local source and were more sensitive to photoperiod, which increased their exposure to spring frosts. Onset of spruce bud flush is driven by spring frosts and photoperiod, but once started, bud phenology responds to temperature. The high photoperiod sensitivity in white spruces could counterbalance climate warming and limit future premature leaf-out, whereas the low photoperiod sensitivity in black spruce should not restrain leaf-out advancement with climate warming. Our results call for adapting the temperature-driven hypotheses of ecophysiological models predicting leaf-out to include spring frost probability.

Simultaneous Real-Time Measurement of Isoprene and 2-Methyl-3-Buten-2-ol Emissions From Trees Using SIFT-MS

The C5 hemiterpenes isoprene and 2-methyl-3-buten-2-ol (MBO) are important biogenic volatiles emitted from terrestrial vegetation. Isoprene is emitted from many plant groups, especially trees such as Populus, while emission of MBO is restricted to certain North American conifers, including species of Pinus. MBO is also a pheromone emitted by several conifer bark beetles. Both isoprene and MBO have typically been measured by proton-transfer reaction mass spectrometry (PTR-MS), but this method cannot accurately distinguish between them because of their signal overlap. Our study developed a method for using selective ion flow tube mass spectrometry (SIFT-MS) that allows simultaneous on-line measurement of isoprene and MBO by employing different reagent ions. The use of m/z(NO⁺) = 68 u for isoprene and m/z(O₂ ⁺) = 71 u for MBO gave minimal interference between the compounds. We tested the suitability of the method by measuring the emission of young trees of Populus, Picea, and Pinus. Our results largely confirm previous findings that Populus nigra, Picea glauca, and Picea abies emit isoprene and Pinus ponderosa emits MBO, but we also found MBO to be emitted by Picea abies. Thus SIFT-MS provides a reliable, easy to use, on-line measuring tool to distinguish between isoprene and MBO. The method should be of use to atmospheric chemists, tree physiologists and forest entomologists, among others.

Seasonal movement of chloroplasts in mesophyll cells of two Picea species

Cold acclimation in evergreen conifers of temperate zone is associated with seasonal structural changes of mesophyll cells. Photoprotective reactions include the movement of chloroplasts from summer position when they are located along the cell walls to winter arrangement with their aggregation in one part of the cell. Special spatial arrangement of mesophyll in Picea species with chloroplasts located along the two opposite cell walls causes the very specific pattern of chloroplast movement. To reveal the intracellular apparatus involved in the seasonal organelle position changes, 3D reconstruction of mesophyll cell structure was applied. Two Picea species, P. obovata and P. pungens, from two geographic regions were studied in a 3-year course. The involvement of small transparent vacuoles in the development of cytoplasmic strands penetrating through the central vacuole and connecting two opposite lateral sides of the cell was shown. The nucleus and cytoplasmic organelles including chloroplasts move inwards the strand forming the cytoplasmic conglomerate enclosed by the vacuole at the cell center. Two Picea species have distinct differences in spatial organization of winter mesophyll cells and in structural events leading to its formation. Analysis of Picea species from two geographic regions over 3 years of monitoring reveals dependence of seasonal organelle movement on the dynamics of temperature decline in autumn.

Phenotypic plasticity of stomatal and photosynthetic features of four Picea species in two contrasting common gardens

Global climate change is expected to affect mountain ecosystems significantly. Phenotypic plasticity, the ability of any genotype to produce a variety of phenotypes under different environmental conditions, is critical in determining the ability of species to acclimate to current climatic changes. Here, to simulate the impact of climate change, we compared the physiology of species of the genus Picea from different provenances and climatic conditions and quantified their phenotypic plasticity index (PPI) in two contrasting common gardens (dry vs. wet), and then considered phenotypic plastic effects on their future adaptation. The mean PPI of the photosynthetic features studied was higher than that of the stomatal features. Species grown in the arid and humid common gardens were differentiated: the stomatal length (SL) and width (SW) on the adaxial surface, the transpiration rate (Tr) and leaf mass per area (LMA) were more highly correlated with rainfall than other traits. There were no significant relationships between the observed plasticity and the species' original habitat, except in P. crassifolia (from an arid habitat) and P. asperata (from a humid habitat). Picea crassifolia exhibited enhanced instantaneous efficiency of water use (PPI = 0.52) and the ratio of photosynthesis to respiration (PPI = 0.10) remained constant; this species was, therefore, considered to the one best able to acclimate when faced with the effects of climate change. The other three species exhibited reduced physiological activity when exposed to water limitation. These findings indicate how climate change affects the potential roles of plasticity in determining plant physiology, and provide a basis for future reforestation efforts in China.

LTR Retrotransposons Show Low Levels of Unequal Recombination and High Rates of Intraelement Gene Conversion in Large Plant Genomes

The accumulation and removal of transposable elements (TEs) is a major driver of genome size evolution in eukaryotes. In plants, long terminal repeat (LTR) retrotransposons (LTR-RTs) represent the majority of TEs and form most of the nuclear DNA in large genomes. Unequal recombination (UR) between LTRs leads to removal of intervening sequence and formation of solo-LTRs. UR is a major mechanism of LTR-RT removal in many angiosperms, but our understanding of LTR-RT-associated recombination within the large, LTR-RT-rich genomes of conifers is quite limited. We employ a novel read-based methodology to estimate the relative rates of LTR-RT-associated UR within the genomes of four conifer and seven angiosperm species. We found the lowest rates of UR in the largest genomes studied, conifers and the angiosperm maize. Recombination may also resolve as gene conversion, which does not remove sequence, so we analyzed LTR-RT-associated gene conversion events (GCEs) in Norway spruce and six angiosperms. Opposite the trend for UR, we found the highest rates of GCEs in Norway spruce and maize. Unlike previous work in angiosperms, we found no evidence that rates of UR correlate with retroelement structural features in the conifers, suggesting that another process is suppressing UR in these species. Recent results from diverse eukaryotes indicate that heterochromatin affects the resolution of recombination, by favoring gene conversion over crossing-over, similar to our observation of opposed rates of UR and GCEs. Control of LTR-RT proliferation via formation of heterochromatin would be a likely step toward large genomes in eukaryotes carrying high LTR-RT content.

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.

[Distribution changes and refugia of three spruce taxa since the last interglacial.]

Based on the current distribution information and 19 environmental variables data, we used the maximum entropy model to simulate the suitable distribution of Picea likiangensis var. likiangensis, P. purpurea and P. wilsonii in the last interglacial, Last Glacial Maximum, Mid- Holocene and present. The results from such modelling were validated by pollen data. We analyzed the relationship between species distribution dynamics and climate change, and then speculated the cryptic refugia of those species. Both the areas under the receiver operating characteristic curves and the verification results from pollen data indicated high accuracy of the model results. Results showed that isothermality was the most important factor influencing the distribution of P. likiangensis var. likiangensis and mean temperature of the warmest quarter was the most important for the distribution of both P. purpurea and P. wilsonii. Temperature was more important than precipitation in driving species distributions. Three species expanded their distribution ranges in Last Glacial Maximum due to their cold-adapted ecological habitat and the deep canyon topography feature which benefited their migration. There might be refugia of both P. wilsonii and P. purpurea in last interglacial, and they respectively located in Shennongjia Mountain in Hubei and Erlang Mountain and its nearby mountains in Sichuan. Our results, to some extent, made accurate prediction of the suitable distribution of three spruce species in the key periods since last interglacial, and speculated refugia of P. purpurea and P. willsonii. Our findings provided reference for better understanding of the formation mechanism of the present distribution of Picea and prediction of distribution changes in the future and sustainable management and protection of three spruce species.

Challenges of pheromone-based mating disruption of Cydia strobilella and Dioryctria abietella in spruce seed orchards

Seed orchards function as the primary source of high-quality seeds for reforestation in many European countries, but their seed yields can be severely reduced due to seed- and cone-feeding insects. We evaluated various parameters of pheromone-based mating disruption for control of the moths Cydia strobilella and Dioryctria abietella, which are major pests in European Picea abies seed orchards. We applied different types of pheromone dispensers (rubber septa or wax droplets) at different densities and heights, and with different amounts of active components, covering whole orchards or part of an orchard. The efficacy of the treatment was evaluated by analysing male captures in pheromone-baited assessment traps, and presence of larvae in cones. A dramatic decrease (94-100%) in capture of males in traps occurred in all pheromone-treated plots compared to control plots for both moth species. In contrast, a subsequent reduction in larval numbers in cones was only achieved when wax droplets were used as the dispensing formulation at high density and at the highest pheromone dose tested, and only in half of the trials for each pest species. Electrophysiological recordings using antennae of male C. strobilella indicated elevated pheromone concentrations in a treated plot versus a control plot. Our results show that mating disruption has potential to reduce cone damage in spruce seed orchards caused by C. strobilella and D. abietella, but optimisation of the technique is required to achieve consistent and efficient population suppression of these pests.

In Silico Analysis of Small RNAs Suggest Roles for Novel and Conserved miRNAs in the Formation of Epigenetic Memory in Somatic Embryos of Norway Spruce

Epigenetic memory in Norway spruce affects the timing of bud burst and bud set, vitally important adaptive traits for this long-lived forest species. Epigenetic memory is established in response to the temperature conditions during embryogenesis. Somatic embryogenesis at different epitype inducing (EpI) temperatures closely mimics the natural processes of epigenetic memory formation in seeds, giving rise to epigenetically different clonal plants in a reproducible and predictable manner, with respect to altered bud phenology. MicroRNAs (miRNAs) and other small non-coding RNAs (sRNAs) play an essential role in the regulation of plant gene expression and may affect this epigenetic mechanism. We used NGS sequencing and computational in silico methods to identify and profile conserved and novel miRNAs among small RNAs in embryogenic tissues of Norway spruce at three EpI temperatures (18, 23 and 28°C). We detected three predominant classes of sRNAs related to a length of 24 nt, followed by a 21–22 nt class and a third 31 nt class of sRNAs. More than 2100 different miRNAs within the prevailing length 21–22 nt were identified. Profiling these putative miRNAs allowed identification of 1053 highly expressed miRNAs, including 523 conserved and 530 novels. 654 of these miRNAs were found to be differentially expressed (DEM) depending on EpI temperature. For most DEMs, we defined their putative mRNA targets. The targets represented mostly by transcripts of multiple-repeats proteins, like TIR, NBS-LRR, PPR and TPR repeat, Clathrin/VPS proteins, Myb-like, AP2, etc. Notably, 124 DE miRNAs targeted 203 differentially expressed epigenetic regulators. Developing Norway spruce embryos possess a more complex sRNA structure than that reported for somatic tissues. A variety of the predicted miRNAs showed distinct EpI temperature dependent expression patterns. These putative EpI miRNAs target spruce genes with a wide range of functions, including genes known to be involved in epigenetic regulation, which in turn could provide a feedback process leading to the formation of epigenetic marks. We suggest that TIR, NBS and LRR domain containing proteins could fulfill more general functions for signal transduction from external environmental stimuli and conversion them into molecular response. Fine-tuning of the miRNA production likely participates in both developmental regulation and epigenetic memory formation in Norway spruce.

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.

Identification of Norway Spruce MYB-bHLH-WDR Transcription Factor Complex Members Linked to Regulation of the Flavonoid Pathway

Transcription factors (TFs) forming MYB-bHLH-WDR complexes are known to regulate the biosynthesis of specialized metabolites in angiosperms through an intricate network. These specialized metabolites participate in a wide range of biological processes including plant growth, development, reproduction as well as in plant immunity. Studying the regulation of their biosynthesis is thus essential. While MYB (TFs) have been previously shown to control specialized metabolism (SM) in gymnosperms, the identity of their partners, in particular bHLH or WDR members, has not yet been revealed. To gain knowledge about MYB-bHLH-WDR transcription factor complexes in gymnosperms and their regulation of SW, we identified two bHLH homologs of AtTT8, six homologs of the MYB transcription factor AtTT2 and one WDR ortholog of AtTTG1 in Norway spruce. We investigated the expression levels of these genes in diverse tissues and upon treatments with various stimuli including methyl-salicylate, methyl-jasmonate, wounding or fungal inoculation. In addition, we also identified protein-protein interactions among different homologs of MYB, bHLH and WDR. Finally, we generated transgenic spruce cell lines overexpressing four of the Norway spruce AtTT2 homologs and observed differential regulation of genes in the flavonoid pathway and flavonoid contents.

Temperature Range Shifts for Three European Tree Species over the Last 10,000 Years

We quantified the degree to which the relationship between the geographic distribution of three major European tree species, Abies alba, Fagus sylvatica and Picea abies and January temperature (Tjan) has remained stable over the past 10,000 years. We used an extended data-set of fossil pollen records over Europe to reconstruct spatial variation in Tjan values for each 1000-year time slice between 10,000 and 3000 years BP (before present). We evaluated the relationships between the occurrences of the three species at each time slice and the spatially interpolated Tjan values, and compared these to their modern temperature ranges. Our results reveal that F. sylvatica and P. abies experienced Tjan ranges during the Holocene that differ from those of the present, while A. alba occurred over a Tjan range that is comparable to its modern one. Our data suggest the need for re-evaluation of the assumption of stable climate tolerances at a scale of several thousand years. The temperature range instability in our observed data independently validates similar results based exclusively on modeled Holocene temperatures. Our study complements previous studies that used modeled data by identifying variation in frequencies of occurrence of populations within the limits of suitable climate. However, substantial changes that were observed in the realized thermal niches over the Holocene tend to suggest that predicting future species distributions should not solely be based on modern realized niches, and needs to account for the past variation in the climate variables that drive species ranges.

The long-range transport of Pinaceae pollen: an example in Kraków (southern Poland)

High Pinaceae pollen concentrations in the air and on the surface of puddles before the main pollen season started were observed in Kraków (southern Poland) in May 2013. The paper presents the results of detailed studies of the composition and source of the "yellow rain" in 2013, and as a comparison, the Pinaceae pollen concentrations and samples collected from the ground surface in 2014 were considered. The air samples were collected using the volumetric method (Hirst-type device), while pollen grains sampled from the ground surface were processed using a modified Erdtman acetolysis method. Finally, all samples were studied using a light microscope. In 2013, the period of higher Abies, Picea and Pinus pollen concentrations was observed from the 5 to 12 of May, earlier than the main pollen season occurred. The presence of rainfall on the 12 and 13 of May 2013 caused the pollen deposition on the ground surface, where the prevalence of Pinaceae pollen was found. The synoptic situation and the analysis of the back-trajectories and air mass advection at the beginning of May 2013 indicated that Pinaceae pollen grains could have been transported from Ukraine, Romania, Hungary and Slovakia. In contrast, Pinaceae pollen grains deposited on the ground surface as a "yellow" film in May 2014, originated from local sources.

Late Pleistocene climate change promoted divergence between Picea asperata and P. crassifolia on the Qinghai-Tibet Plateau through recent bottlenecks

Divergence during the early stage of speciation can be driven by a population bottleneck via reduced gene flow and enhanced lineage sorting. In this study, we aimed to examine whether such bottlenecks occurred during the initial speciation of two closely related spruce species Picea asperata and P. crassifolia occurring on the Qinghai–Tibet Plateau (QTP). We analyzed sequences of three chloroplast, two mitochondrial DNA fragments and a further 13 nuclear loci from 216 individuals of the two species. Both species showed a low level of genetic diversity in contrast to other congeners occurring in the QTP and adjacent regions. The estimated population sizes of P. asperata and P. crassifolia are less than the ancestral population size before splitting. These results together with multiple statistical tests (Tajima's D, Fu and Li's D* and F*) suggest that these two species underwent recent bottlenecks. Based on approximate Bayesian computation (ABC), we also determined that the period of the population shrinkage was consistent with the interspecific divergence during the late Pleistocene. The reduced population sizes and the divergent selection may together have triggered the initial divergence under high gene flow between these two species. Our results therefore highlight the importance of climatic oscillations during the late Pleistocene in promoting speciation through changing demographic sizes of the ancestral species on the QTP and in adjacent regions.

Morphometric traits capture the climatically driven species turnover of 10 spruce taxa across China

This study explored the relative roles of climate and phylogenetic background in driving morphometric trait variation in 10 spruce taxa in China. The study further addressed the hypothesis that these variations are consistent with species turnover on climatic gradients. Nine morphometric traits of leaves, seed cones, and seeds for the 10 studied spruce taxa were measured at 504 sites. These data were analyzed in combination with species DNA sequences from NCBI GenBank. We detected the effects of phylogeny and climate through trait-variation-based K statistics and phylogenetic eigenvector regression (PVR) analyses. Multivariate analyses were performed to detect trait variation along climatic gradients with species replacement. The estimated K-values for the nine studied morphometric traits ranged from 0.19 to 0.68, and the studied environmental variables explained 39-83% of the total trait variation. Trait variation tended to be determined largely by a temperature gradient varying from wet-cool climates to dry-warm summers and, additionally, by a moisture gradient. As the climate became wetter and cooler, spruce species tended to be replaced by other spruces with smaller needle leaves and seeds but larger cones and seed scales. A regression analysis showed that spruce species tended to be successively replaced by other species, along the gradient, although the trends observed within species were not necessarily consistent with the overall trend. The climatically driven replacement of the spruces in question could be well indicated by the between-species variation in morphometric traits that carry lower phylogenetic signal. Between-species variation in these traits is driven primarily by climatic factors. These species demonstrate a narrower ecological amplitude in temperature but wider ranges on the moisture gradient.

The Ty1-copia LTR retroelement family PARTC is highly conserved in conifers over 200 MY of evolution

Long Terminal Repeat retroelements (LTR-RTs) are a major component of many plant genomes. Although well studied and described in angiosperms, their features and dynamics are poorly understood in gymnosperms. Representative complete copies of a Ty1-copia element isolate in Picea abies and named PARTC were identified in six other conifer species (Picea glauca, Pinus sylvestris, Pinus taeda, Abies sibirica, Taxus baccata and Juniperus communis) covering more than 200 million years of evolution. Here we characterized the structure of this element, assessed its abundance across conifers, studied the modes and timing of its amplification, and evaluated the degree of conservation of its extant copies at nucleotide level over distant species. We demonstrated that the element is ancient, abundant, widespread and its paralogous copies are present in the genera Picea, Pinus and Abies as an LTR-RT family. The amplification leading to the extant copies of PARTC occurred over long evolutionary times spanning 10s of MY and mostly took place after the speciation of the conifers analyzed. The level of conservation of PARTC is striking and may be explained by low substitution rates and limited removal mechanisms for LTR-RTs. These PARTC features and dynamics are representative of a more general scenario for LTR-RTs in gymnosperms quite different from that characterizing the vast majority of LTR-RT elements in angiosperms.

MADS goes genomic in conifers: towards determining the ancestral set of MADS-box genes in seed plants

BACKGROUND AND AIMS

MADS-box genes comprise a gene family coding for transcription factors. This gene family expanded greatly during land plant evolution such that the number of MADS-box genes ranges from one or two in green algae to around 100 in angiosperms. Given the crucial functions of MADS-box genes for nearly all aspects of plant development, the expansion of this gene family probably contributed to the increasing complexity of plants. However, the expansion of MADS-box genes during one important step of land plant evolution, namely the origin of seed plants, remains poorly understood due to the previous lack of whole-genome data for gymnosperms.

METHODS

The newly available genome sequences of Picea abies, Picea glauca and Pinus taeda were used to identify the complete set of MADS-box genes in these conifers. In addition, MADS-box genes were identified in the growing number of transcriptomes available for gymnosperms. With these datasets, phylogenies were constructed to determine the ancestral set of MADS-box genes of seed plants and to infer the ancestral functions of these genes.

KEY RESULTS

Type I MADS-box genes are under-represented in gymnosperms and only a minimum of two Type I MADS-box genes have been present in the most recent common ancestor (MRCA) of seed plants. In contrast, a large number of Type II MADS-box genes were found in gymnosperms. The MRCA of extant seed plants probably possessed at least 11-14 Type II MADS-box genes. In gymnosperms two duplications of Type II MADS-box genes were found, such that the MRCA of extant gymnosperms had at least 14-16 Type II MADS-box genes.

CONCLUSIONS

The implied ancestral set of MADS-box genes for seed plants shows simplicity for Type I MADS-box genes and remarkable complexity for Type II MADS-box genes in terms of phylogeny and putative functions. The analysis of transcriptome data reveals that gymnosperm MADS-box genes are expressed in a great variety of tissues, indicating diverse roles of MADS-box genes for the development of gymnosperms. This study is the first that provides a comprehensive overview of MADS-box genes in conifers and thus will provide a framework for future work on MADS-box genes in seed plants.

Pinus flexilis and Picea engelmannii share a simple and consistent needle endophyte microbiota with a potential role in nitrogen fixation

Conifers predominantly occur on soils or in climates that are suboptimal for plant growth. This is generally attributed to symbioses with mycorrhizal fungi and to conifer adaptations, but recent experiments suggest that aboveground endophytic bacteria in conifers fix nitrogen (N) and affect host shoot tissue growth. Because most bacteria cannot be grown in the laboratory very little is known about conifer–endophyte associations in the wild. Pinus flexilis (limber pine) and Picea engelmannii (Engelmann spruce) growing in a subalpine, nutrient-limited environment are potential candidates for hosting endophytes with roles in N₂ fixation and abiotic stress tolerance. We used 16S rRNA pyrosequencing to ask whether these conifers host a core of bacterial species that are consistently associated with conifer individuals and therefore potential mutualists. We found that while overall the endophyte communities clustered according to host species, both conifers were consistently dominated by the same phylotype, which made up 19–53% and 14–39% of the sequences in P. flexilis and P. engelmannii, respectively. This phylotype is related to Gluconacetobacter diazotrophicus and other N₂ fixing acetic acid bacterial endophytes. The pattern observed for the P. flexilis and P. engelmannii needle microbiota—a small number of major species that are consistently associated with the host across individuals and species—is unprecedented for an endophyte community, and suggests a specialized beneficial endophyte function. One possibility is endophytic N fixation, which could help explain how conifers can grow in severely nitrogen-limited soil, and why some forest ecosystems accumulate more N than can be accounted for by known nitrogen input pathways.

Increased water-use efficiency does not lead to enhanced tree growth under xeric and mesic conditions

Higher atmospheric CO2 concentrations (c(a)) can under certain conditions increase tree growth by enhancing photosynthesis, resulting in an increase of intrinsic water-use efficiency (i WUE) in trees. However, the magnitude of these effects and their interactions with changing climatic conditions are still poorly understood under xeric and mesic conditions. We combined radial growth analysis with intra- and interannual δ(13)C and δ(18)O measurements to investigate growth and physiological responses of Larix decidua, Picea abies, Pinus sylvestris, Pinus nigra and Pseudotsuga menziesii in relation to rising c(a) and changing climate at a xeric site in the dry inner Alps and at a mesic site in the Swiss lowlands. (i)WUE increased significantly over the last 50 yr by 8-29% and varied depending on species, site water availability, and seasons. Regardless of species and increased (i)WUE, radial growth has significantly declined under xeric conditions, whereas growth has not increased as expected under mesic conditions. Overall, drought-induced stomatal closure has reduced transpiration at the cost of reduced carbon uptake and growth. Our results indicate that, even under mesic conditions, the temperature-induced drought stress has overridden the potential CO2 'fertilization' on tree growth, hence challenging today's predictions of improved forest productivity of temperate forests.

The landscape of nucleotide polymorphism among 13,500 genes of the conifer picea glauca, relationships with functions, and comparison with medicago truncatula

Gene families differ in composition, expression, and chromosomal organization between conifers and angiosperms, but little is known regarding nucleotide polymorphism. Using various sequencing strategies, an atlas of 212k high-confidence single nucleotide polymorphisms (SNPs) with a validation rate of more than 92% was developed for the conifer white spruce (Picea glauca). Nonsynonymous and synonymous SNPs were annotated over the corresponding 13,498 white spruce genes representative of 2,457 known gene families. Patterns of nucleotide polymorphisms were analyzed by estimating the ratio of nonsynonymous to synonymous numbers of substitutions per site (A/S). A general excess of synonymous SNPs was expected and observed. However, the analysis from several perspectives enabled to identify groups of genes harboring an excess of nonsynonymous SNPs, thus potentially under positive selection. Four known gene families harbored such an excess: dehydrins, ankyrin-repeats, AP2/DREB, and leucine-rich repeat. Conifer-specific sequences were also generally associated with the highest A/S ratios. A/S values were also distributed asymmetrically across genes specifically expressed in megagametophytes, roots, or in both, harboring on average an excess of nonsynonymous SNPs. These patterns confirm that the breadth of gene expression is a contributing factor to the evolution of nucleotide polymorphism. The A/S ratios of Medicago truncatula genes were also analyzed: several gene families shared between P. glauca and M. truncatula data sets had similar excess of synonymous or nonsynonymous SNPs. However, a number of families with high A/S ratios were found specific to P. glauca, suggesting cases of divergent evolution at the functional level.

Evolutionary history of Purple cone spruce (Picea purpurea) in the Qinghai-Tibet Plateau: homoploid hybrid origin and Pleistocene expansion

Hybridization and introgression can play an important role in speciation. Here, we examine their roles in the origin and evolution of Picea purpurea, a diploid spruce species occurring on the Qinghai-Tibet Plateau (QTP). Phylogenetic relationships and ecological differences between this species and its relatives, P. schrenkiana, P. likiangensis and P. wilsonii, are unclear. To clarify them, we surveyed sequence variation within and between them for 11 nuclear loci, three chloroplast (cp) and two mitochondrial (mt) DNA fragments, and examined their ecological requirements using ecological niche modelling. Initial analyses based on 11 nuclear loci rejected a close relationship between P. schrenkiana and P. purpurea. BP&P tests and ecological niche modelling indicated substantial divergence between the remaining three species and supported the species status of P. purpurea, which contained many private alleles as expected for a well-established species. Sequence variation for cpDNA and mtDNA suggested a close relationship between P. purpurea and P. wilsonii, while variation at the nuclear se1364 gene suggested P. purpurea was more closely related to P. likiangensis. Analyses of genetic divergence, Bayesian clustering and model comparison using approximate Bayesian computation (ABC) of nuclear (nr) DNA variation all supported the hypothesis that P. purpurea originated by homoploid hybrid speciation from P. wilsonii and P. likiangensis. The ABC analysis dated the origin of P. purpurea at the Pleistocene, and the estimated hybrid parameter indicated that 69% of its nuclear composition was contributed by P. likiangensis and 31% by P. wilsonii. Our results further suggested that during or immediately following its formation, P. purpurea was subject to organelle DNA introgression from P. wilsonii such that it came to possess both mtDNA and cpDNA of P. wilsonii. The estimated parameters indicated that following its origin, P. purpurea underwent an expansion during/after the largest Pleistocene glaciation recorded for the QTP.

Review of the Parasa undulata (Cai, 1983) species group with the first conifer-feeding larva for Limacodidae and descriptions of two new species from China and Taiwan (Lepidoptera, Limacodidae)

Although the caterpillars are well-known for the stings and magnificent coloration, the systematics of Limacodidae is historically neglected and chaotic due to the difficulty in matching the larval with adult stages as well as the very conservative and convergent adult morphology. One of the biggest taxonomic problems surrounds a collective group from Southeastern Asia, termed the “green limacodid moths”, which harbours at least 90 species placed in the genus Parasa Walker, 1859 and 14 “subunits”. The P. undulata group was previously composed of 3 species from China and Taiwan, and characterized only by wing pattern. This species group is extensively studied herein with two new species described, i.e. P. viridiflammasp. n. (Taiwan) and P. minwangisp. n. (S. China), and discovery of female genitalia of three species, presenting new phylogenetic insights in this potentially paraphyletic genus. In addition, one limacodid larva was found to be feeding exclusively on Picea (Pinaceae) in Taiwan. Its identity, Parasa pygmy Solovyev, 2010 in P. undulata group, is confirmed through matching its COI sequence to the adult. This discovery is also biologically significant because the previous known host breadth of Parasa was of polyphagy on various angiosperm plant families. This case, therefore, represents the first record of conifer-feeding behavior in this family as well as the first of specialized herbivory in the genus. Meanwhile, the background match between Picea leaves and larval coloration is shared with other Picea-feeding insects. This phenomenon is worth of further investigation in the aspect of convergent evolution of crypsis associated with a particular plant.

Drought response of five conifer species under contrasting water availability suggests high vulnerability of Norway spruce and European larch

The ability of tree species to cope with anticipated decrease in water availability is still poorly understood. We evaluated the potential of Norway spruce, Scots pine, European larch, black pine, and Douglas-fir to withstand drought in a drier future climate by analyzing their past growth and physiological responses at a xeric and a mesic site in Central Europe using dendroecological methods. Earlywood, latewood, and total ring width, as well as the δ(13) C and δ(18) O in early- and latewood were measured and statistically related to a multiscalar soil water deficit index from 1961 to 2009. At the xeric site, δ(13) C values of all species were strongly linked to water deficits that lasted longer than 11 months, indicating a long-term cumulative effect on the carbon pool. Trees at the xeric site were particularly sensitive to soil water recharge in the preceding autumn and early spring. The native species European larch and Norway spruce, growing close to their dry distribution limit at the xeric site, were found to be the most vulnerable species to soil water deficits. At the mesic site, summer water availability was critical for all species, whereas water availability prior to the growing season was less important. Trees at the mesic were more vulnerable to water deficits of shorter duration than the xeric site. We conclude that if summers become drier, trees growing on mesic sites will undergo significant growth reductions, whereas at their dry distribution limit in the Alps, tree growth of the highly sensitive spruce and larch may collapse, likely inducing dieback and compromising the provision of ecosystem services. However, the magnitude of these changes will be mediated strongly by soil water recharge in winter and thus water availability at the beginning of the growing season.

Genetic and morphological structure of a spruce hybrid (Picea sitchensis x P. glauca) zone along a climatic gradient

PREMISE OF THE STUDY

Historic colonization and contemporary evolutionary processes contribute to patterns of genetic variation and differentiation among populations. However, separating the respective influences of these processes remains a challenge, particularly for natural hybrid zones, where standing genetic variation may result from evolutionary processes both preceding and following contact, influencing the evolutionary trajectory of hybrid populations. Where adaptation to novel environments may be facilitated by interspecific hybridization, teasing apart these processes will have practical implications for forest management in changing environments.

METHODS

We evaluated the neutral genetic architecture of the Picea sitchensis (Sitka spruce) × P. glauca (white spruce) hybrid zone along the Nass and Skeena river valleys in northwestern British Columbia using chloroplast, mitochondrial, and nuclear microsatellite markers, in combination with cone morphological traits.

KEY RESULTS

Sitka spruce mitotype "capture", evidenced by this species dominating the maternal lineage, is consistent with earlier colonization of the region by Sitka spruce. This "capture" differs from the spatial distribution of chloroplast haplotypes, indicating pollen dispersal and its contribution to geographic structure. Genetic ancestry, based on nuclear markers, was strongly influenced by climate and geography. Highly parallel results for replicate transects along environmental gradients provide support for the bounded hybrid superiority model of hybrid zone maintenance. •

CONCLUSIONS

This broad-scale analysis of neutral genetic structure indicates the importance of historic and contemporary gene flow, environmental selection, and their interaction in shaping neutral genetic variation within this hybrid zone, informative to seed transfer development and reforestation for future climates.

Estimating photosynthetic rate and annual carbon gain in conifers from specific leaf weight and leaf biomass

Canopy photosynthesis is difficult to measure directly or to predict with complex models demanding knowledge of seasonal variation in environmental and physiological properties of the canopy. Trees in particular offer a challenge with their large, aerodynamically rough and seasonally-changing canopy properties. In this paper we assess the possibility of using specific leaf weight to predict seasonal and annual net photosynthetic rate in deciduous (Larix sp.) and evergreen (Picea abies) conifers.Annual photosynthetic rate and specific leaf weight of different positions of the crown in both species were highly correlated (r ²=0.930). Annual carbon uptake by different segments in a mature P. abies crown was closely related to leaf biomass. The relationship was improved by adjusting the leaf biomass of each segment in regard to its specific leaf weight relative to the maximum found in the canopy. The adjustment accounted for associated differences in photosynthetic activity. This combined structural index (leaf biomassxrelative specific leaf weight) could, when calibrated, predict the total annual carbon uptake by different parts of the crown. If direct measurements of photosynthesis are not available, the combined structural index may still serve as a comparative estimator of annual carbon uptake.

Carbon dioxide assimilation and growth of red spruce (Picea rubens Sarg.) seedlings in response to ozone, precipitation chemistry, and soil type

The influence of ozone, mist chemistry, rain chemistry, and soil type on CO₂ assimilation and growth of red spruce (Picea rubens Sarg.) seedlings was investigated over a 4-month period under controlled laboratory and glasshouse conditions. Growth was evaluated through interval estimates of aboveground relative growth rates (RGR) and the partitioning of biomass components at harvest to root, stem, and needle fractions. Precipitation chemistry treatments and O₃ exposure dynamics were based on reported characteristics of air chemistry and/or deposition in high-elevation forests of eastern North America. The two soils were collected from Camels Hump in the Green Mountains of Vermont and Acadia National Park on the Maine coast. Soil from Acadia had higher organic content, higher levels of extractable base cations, and lower levels of extractable aluminum and heavy metals. The only treatment variables that consistently influenced the growth of P. rubens were soil type and rain chemistry. In comparison with seedlings grown in soil from Acadia National Park, those grown in Camels Hump soil had significantly less needle (27%), stem (33%), and root (26%) biomass at harvest and statistically lower aboveground RGR within 2 months after initiation of the treatments. Seedlings grown in Camels Hump soil had significantly higher levels of aluminum (6.5X), copper (1.4X), and nickel (2.7X) in new needle tissue. The only influence of precipitation chemistry on the growth of P. rubens was a pattern of greater root and shoot biomass in seedlings experiencing the more acidic rain treatments. Interactive effects among the main treatment variables (e.g., acidic mist and O₃, acidic rain and soil type) on seedling growth were not notable. Rates of CO₂ assimilation and transpiration on a per gram needle dry weight basis [mol·g^-1·s^-1] were not influenced by any of the main treatment variables or their interaction. Because neither soil type nor precipitation chemistry influenced the efficiency of CO₂ assimilation per gram dry weight of needle tissue, the physiological mechanism underlying the growth response of P. rubens is attributed to a change in either whole-plant allocation of carbon resources or a direct toxic effect in the rhizosphere on root growth.