Measurement(s)Cone length • Cone weight • Cone width • Height • Viable seed number • Viable seed weight • stem diameter • canopy width • needle length • number of buds • budburst timing • height increment • stem diameter increment • growth cessation • seed viabilityTechnology Type(s)sliding caliper • balance • ruler • natural number • measuring tape • budburst scaleFactor Type(s)field site location • tree family • tree populationSample Characteristic - OrganismPinus sylvestrisSample Characteristic - LocationScotland
Multisite common garden experiments, exposing common pools of genetic diversity to a range of environments, allow quantification of plastic and genetic components of trait variation. For tree species, such studies must be long term as they typically only express mature traits after many years. As well as evaluating standing genetic diversity, these experiments provide an ongoing test of genetic variation against changing environmental conditions and form a vital resource for understanding how species respond to abiotic and biotic variation. Finally, quantitative assessments of phenotypic variation are essential to pair with rapidly accumulating genomic data to advance understanding of the genetic basis of trait variation, and its interaction with climatic change. We describe a multisite, population-progeny, common garden experiment of the economically and ecologically important tree species, Scots pine, collected from across its native range in Scotland and grown in three contrasting environments. Phenotypic traits, including height, stem diameter and budburst were measured over 14 growing seasons from nursery to field site. The datasets presented have a wide range of applications.
Background & SummaryThe need for comprehensive empirical assessments of genetic variation in tree species has never been greater. There is great interest around the world in growing more trees1, for their carbon sequestration abilities in the race for ‘net zero’ carbon emissions, to arrest biodiversity loss and forest decline2, and to manage watersheds3, as well as for the products they provide. However, there is considerable uncertainty in the accuracy of forecasts of future climate and the responses of tree species to those changes4. Intraspecific genetic variation and phenotypic plasticity will play key roles in determining how resilient existing and new tree populations are to the challenges ahead5. We are in urgent need of robust empirical data to calibrate the relationships between genetic variation within species and the environmental variables that will define future climates6. In parallel, the genomic revolution has provided a dramatic increase in the accessibility and scale of molecular data, and trees can now be genotyped faster, more cheaply and in greater number than ever before. However, there is a limit to what these new data can tell us without objective evaluation of the associated phenotypes, for which common garden and reciprocal transplant approaches remain key experimental tools7,8. We aimed to link these genomic and common garden approaches to better understand the genetic basis of phenotypic variation in trees, and to improve forecasting of how tree species will respond to climate change. Here we describe the research platform that we established to conduct these urgently needed long-term studies.Scots pine (Pinus sylvestris L.) is globally very widely distributed, occurring predominantly from eastern Europe and Scandinavia to eastern Siberia, but with substantial populations in Scotland, southern Europe, Turkey and the Caucasus9. In Scotland, the remnant native populations, known locally as the Caledonian pinewoods, are typically small and highly fragmented and are distributed across a highly heterogeneous landscape that varies from oceanic (mild, very wet) environments in the west to more continental (drier, colder) in the east (Fig. 1). This steep gradient of variation on a short spatial scale bears comparison to environmental gradients over much wider spatial scales across Europe (see Fig. 4 in Metzger et al.10). To further the conservation of the pinewoods, specific action plans have been developed, including management of seed movements through a series of seven seed zones (Fig. 211), the context and motivation for which are covered by Salmela et al.12.Fig. 1Ordination of 21 source populations in environmental space. Gray dots: global Scots pine distribution (distribution - EU-Forest database18 and associated climate data (CHELSA database19, https://chelsa-climate.org/). Blue dots: Scots pine distribution in the United Kingdom. Orange dots: populations sampled for the experiment; letter codes match those in Table 1. Yellow dots: trial sites.Fig. 2Location of the 21 source populations sampled for the experiment; letter codes match those in Table 1. Also shown are James Hutton Institute (JHI), location of eastern nursery (NE) and glasshouse (NG); three field sites - Inverewe (FW, also location of western nursery: NW), Glensaugh (FE) and Yair (FS).A key knowledge gap remained in the early 21st century, namely the extent to which the current seed zone system for managing seed sourcing for the native populations of Scottish Scots pine reflects genuine genetic differences, and whether the zones conform to patterns of adaptation which are of relevance to ecological conservation and forestry when assessed using traits such as growth, mortality, phenology and resistance to pests and pathogens. To address this gap, as well as to evaluate the likely responses of indigenous Scots pine to ongoing climate change, an extensive, long-term assessment of genetic variation in Scots pine under contrasting environmental conditions, was founded by the Macaulay Institute (now the James Hutton Institute). A collaboration with the UK Centre for Ecology & Hydrology and Forest Research was subsequently created to develop the studies of adaptation to climate and disease resistance, to apply and advance emerging methods for assessing genetic variation, and to secure the continuation of the study. At the outset, an explicit objective was to sample genetic diversity widely without favouring any particular form or trait, and to undertake a wide ranging evaluation of traits covering different life history dimensions as far as practically possible. The common environment study was also intended to provide a long term experimental platform to facilitate future studies of the basis for variation in the extended phenotype of Scots pine including associated assemblages of organisms and community function.The following describes the origins, design and initial measurements of a multi-site experiment in Scots pine, including protocols adopted during both the initial nursery phase and the final field experiment.MethodsSeed sampling and germinationSeed from ten trees from each of 21 native Scottish Scots pine populations (Table 1) were collected in March 2007 and germinated at the James Hutton Institute, Aberdeen (latitude 57.133214, longitude −2.158764) in June 2007. Populations were chosen to represent the species’ native range in Scotland and to include three populations from each of the seven seed zones (Fig. 2). There was no selection of seed-trees on the basis of any traits except for the possession of cones on the date of sampling. Ten seed trees were sampled from each population according to a spatial protocol designed to cover a circle of approximately 1 km in diameter located around a central tree. The sampling strategy identified nine points each in a pre-determined random direction from the central point, whilst stratifying the number sampled with increasing distance from the central point in the ratio 1: 3: 5. This strategy avoids over-sampling the areas close to the centre point. For smaller fragments of woodland, or where only a few trees with cones were present, then the directions of the sampled trees from the central tree were maintained to give a wide coverage of the woodland area, but the distances between trees varied but were never closer than 50 m. To break dormancy, seeds were soaked for 24 hours on the benchtop at room temperature, after which they were stored in wet paper towels and refrigerated in darkness at 3–5 °C for approximately 4 weeks. Seeds were kept moist and transferred to room temperature until germination began (approx. 5–7 days), then transplanted to 8 cm × 8 cm × 9 cm, 0.4 L pots filled with Levington’s C2a compost and 1.5 g of Osmocote Exact 16–18 months slow release fertiliser and kept in an unheated glasshouse. The compost was covered with a layer of grit to reduce moss and liverwort growth. Seedlings from the same mother tree are described as a family and are assumed to be half-siblings.Table 1Locations and basic environmental data for the populations sampled for seed to establish the trial. See the maternal traits dataset15 for precise data for each mother tree sampled.PopulationLat. (N)Long. (W)Alt. (m)Area (ha)GSLGDDFMT (°C)JMT (°C)Abernethy (AB)57.213.61311–37024522119901.112.7Allt Cul (AC)57.043.35435–51213145513−1.010.4Amat (AM)57.874.6039–2011812148921.212.3Ballochbuie (BB)56.983.30421–531775116446−1.79.5Beinn Eighe BE)57.635.4017–9118228313293.714.2Black Wood of Rannoch (BW)56.684.37250–321101125411382.113.5Coille Coire Chuilc (CC)56.424.71222–311672269281.612.3Conaglen (CG)56.795.3389–1931892468872.211.7Crannach (CR)56.584.68258–3387023110191.812.6Glen Affric (GA)57.264.92205–29315322107690.911.6Glen Cannich (GC)57.354.95182–3813012127781.011.7Glen Derry (GD)57.033.58426–493235168593−0.511.3Glen Einig (GE)57.964.7645–922724210892.213.2Glen Loy (GL)56.915.13136–219741915410.59.8Glen Tanar (GT)57.022.86289–422156423511052.213.6Loch Clair (LC)57.565.3698–16612627712533.413.7Meggernie (MG)56.584.35254–3852772239161.112.0Rhidorroch (RD)57.894.98138–2201032218401.511.6Rothiemurchus (RM)57.153.77295–329174422410871.413.1Shieldaig (SD)57.505.6344–13210327310933.212.8Strath Oykel (SO)57.984.6135–1601425712762.714.0Alt: altitudinal range sampled within each population; Area: core pinewood area16. Average (1961–2000) climate variables from UK Met Office17 – GSL: growing season length (period bounded by daily mean temperature >5 °C for >5 consecutive days and daily mean temperature <5 °C for >5 consecutive days (after 1 July); GDD: growing degree days (the mean number of degrees by which the air temperature has gone above 5 °C calculated day by day and summed over the year); FMT and JMT: February and July mean temperatures (°C).Experimental design: nurseriesThe full collection consisted of 210 families (10 families from each of 21 populations) each consisting of 24 half sibling progeny (total 5,040 individuals); needle tissue was sampled from each seedling and preserved for long term storage, one needle on silica gel, 2–5 needles at −20 °C. After transfer into pots, 8 seedlings per family were moved to one of three nurseries (total 1,680 seedlings per nursery): outdoors at Inverewe Gardens in western Scotland (nursery in the west of Scotland: coded NW, latitude 57.775714, longitude −5.597181, Fig. 2); outdoors in a fruit cage (to minimise browsing) at the James Hutton Institute in Aberdeen (nursery in the east of Scotland: NE); in an unheated glasshouse at the James Hutton Institute in Aberdeen (nursery in a glasshouse: NG). Trees were arranged in 40 randomised trays (blocks) in each nursery. Each block contained two trees per population (total 42 trees). Watering was automatic in NG, and manually as required for NE and NW. No artificial light was used in any of the nurseries. In May 2010 the seedlings from NG were moved outdoors to Glensaugh in Aberdeenshire (latitude 56.893567, longitude −2.535736). In 2010 all plants were repotted into 19 cm (3 L) pots containing Levingtons CNSE Ericaceous compost with added Osmocote STD 16–18 month slow release fertilizer.Experimental design: field sitesIn 2012 the trees were transplanted to one of three field sites: Yair in the Scottish Borders (field site in the south of Scotland: FS, latitude 55.603625, longitude −2.893025); Glensaugh (field site in the east of Scotland: FE); and Inverewe (field site in the west of Scotland: FW). All trees transplanted to FS were raised in the NG and all but four of the trees transplanted to FE were raised locally in the NE (the remainder were grown in NG). In contrast, following mortality and ‘beating up’ (filling gaps where saplings had died), the FW experiment ultimately contained cohorts of trees raised in each of the three nurseries as follows: 290 grown locally in the NW; 132 were grown in the NG; and 82 were grown in the NE.Site historiesThe Yair site (FS) had previously been used for growing Noble fir (Abies procera) for Christmas trees and Lodgepole pine (Pinus contorta), a section of the former were felled and chipped to create a clear area prior to planting. The planting site is also adjacent to a large block of commercial Sitka spruce (Picea sitchensis) forestry, and the Glenkinnon Burn Site of Special Scientific Interest (SSSI NatureScot site code 736; EU site code 135445), an area of mixed broadleaf woodland. Prior to planting, major areas of tall weeds were strimmed. The site was protected by a deer fence. The experiment was planted 8–11 October 2012. The Glensaugh site (FE) is in Forestry Compartment 3 of the Glensaugh Research Station, adjacent to Cleek Loch. It is thought to have been cleared of Scots pine and Larch (Larix decidua) around 1917, after which it reverted to rough grazing. An attempt to reseed part of the site in the 1980s was unsuccessful and it quickly reverted to rough grazing for a second time. The whole site (within which the experimental area is embedded) was deer fenced and re-planted under the Scottish Rural Development Programme (SRDP) in 2012. The experimental plot was planted up 7–9 March 2012. The Inverewe site (FW) had previously been a Sitka spruce and Lodgepole pine plantation (50:50 mix) that had been clear-felled in 2010 following substantial windthrow. The experimental site was deer fenced in early 2012, and the experiment was planted 12–16 March 2012, followed by beating up on 27–28 March 2013 and 22–24 October 2013. There had been minimal preparation of the site in line with current practice for restocking sites. The experimental site is included in the Inverewe Forest Plan, which included deer fencing of a larger area (2014) around the experimental site. Planting of this area was completed in early 2015, funded by NTS (National Trust for Scotland), although natural regeneration is also taking place.At each site, trees were planted in randomised blocks at 3 m × 3 m spacing. There are four randomised blocks in both FS and FE and three in FW. A guard row of Scots pine trees was planted around the periphery of the blocks and between blocks B and C at FS. Each block comprised one individual from each of eight (of the 10 sampled) families per 21 populations (168 trees). Although most families (N = 159) were represented at each of the three sites, families with insufficient trees (N = 9) were replaced in one site (FS) with a different family from the same population. Each experimental site was designed with redundancy such that, if thinning becomes necessary as the trees mature, then the systematic removal of trees (i.e. trees 1,3,5,7, etc of row 1, and 2,4,6,8, etc of row 2, 1,3,5,7,etc of row 3) will maintain a balanced design of the experiment, with sufficient family and population representation to provide an ongoing experiment with full geographic coverage.The field sites generally experience different climates, with FW typically warmer and wetter and with more growing degree days per year and a much longer growing season than both FE and FS (Table 2). The coldest site with the shortest growing season is generally FE.Table 2Average climatic variables at field sites Glensaugh (FE), Inverewe (FW) and Yair (FS) from planting in 2012 until 2020. Climatic variables are derived from data provided by the Met Office (daily mean, minimum and maximum temperatures and monthly rainfall).MeasureFEFWFSAverage Maximum Daily Air Temperature (deg C)11.4312.7312.33Average Minimum Daily Air Temperature (deg C)4.856.505.03Average Mean Daily Temperature (deg C)8.149.628.68Growing Degree Days (day degrees)1146618301664Growing Season Length (days)2258332280Average Total Monthly Rainfall (mm)91.42138.2382.481Growing degree days: the mean number of degrees by which the daily mean air temperature has gone above 5 °C calculated day by day and summed over the year.2Growing season length: period bounded by daily mean air temperature >5 °C for >5 consecutive days and daily mean air temperature <5 °C for >5 consecutive days (after 1 July).Phenotype assessmentsMaternal traitsFollowing seed collection, a range of traits were measured in the mother trees in order to control for maternal effects in subsequent measurements of their progeny (Table 3). For each mother tree, measurements of height and diameter at breast height (DBH) were taken, and ten cones were collected and assessed in detail. Cone width and length were measured prior to drying the cones (when they were still closed). Cone weight was measured post-drying. Seed removed from each cone was assessed for total weight (after wings had been removed) and for the count and percentage of seeds which were classed as viable (viable seed were those which had both a wing and an obvious seed). No further seed sorting was applied.Table 3Traits assessed in mother trees, cones, seeds (dataset: Maternal), nursery seedlings (dataset: Nursery) and field trials (dataset: Field). Within the datasets, traits are recorded in a single column for each year using the format Code-Year (e.g. absolute height in 2008 = HA08) except for the maternal traits datasets which were all measured in 2007.TraitCodeUnitDatasetYear(s)Cone lengthLemmMaternal2007Cone weightWeGMaternal2007Cone widthWimmMaternal2007Height: absoluteHAMmmmmMaternalNurseryField20072007–2011A2013–2020BHeight: incrementHImmmmNurseryField2008–2011AViable seed numberSNCountMaternal2007Seed viabilitySV%Maternal2007Viable seed weightSWGMaternal2007Stem diameter: absoluteCDAMmmmmMaternalNurseryField20072008–2011A2013–2020DStem diameter: incrementDImmmmNurseryField2009–2011ACanopy widthCWmmNursery2008Growth cessationGCDays since 10 SeptemberNursery2008Needle lengthNLmmNursery2007–2009Number of budsBuCountNursery2008Budburst durationEBDDaysField2015–2019Budburst timingFGBBBTDays since 31 MarchNurseryField20082015–2019Where multiple measurements are made in a single year (i.e. for canopy width and needle length) the suffix “_X” is added to the column header, where X is the measurement number (e.g. canopy width measured in 2007, second measurement = CW07_2). Where multiple stages are recorded in a single year (i.e. for budburst timing) the suffix “_Y” is added to the column header, where Y is the budburst stage (see Fig. 3 for budburst stages).A2011 measurement only in FE and FW.B2013 measurement only in FE and FW.CStem diameter measured at breast height for maternal traits dataset and at base for nursery and field traits datasets.DStem diameter only measured at FS in 2020.EBudburst duration calculated as the number of days between budburst stages 4 and 6.FBudburst timing (as scored using stages described in Fig. 3) is limited to stages 4 to 6 in this dataset, due to the high proportion of missing data for other stages.GSee Fig. 3 for budburst stages.Nursery traitsSeedling phenotype assessments were performed annually from 2007–2010 for three different trait types: phenology (budburst and growth cessation); form (total number of buds, needle length); cumulative growth (stem diameter and height, canopy width). Measurements of tree form and cumulative growth traits were taken after the end of each growing season. Phenology was assessed weekly during the spring and autumn of 2008 for budburst and growth cessation, respectively. Budburst was defined as the number of days from 31 March 2008 to the time when newly emerged green needles were observed (budburst stage 6: Fig. 3). In some seedlings in 2008, a secondary flush of growth occurred from terminal buds that had formed during the summer of that year. Therefore, growth cessation was defined retrospectively as the number of days from 10 September 2008 to the date when a terminal bud had formed on the leading shoot of the seedling, providing no further growth was observed either on the stem below that bud, nor from the bud itself. Canopy width (widest point) was measured at two perpendicular points in the horizontal plane. Needle length was measured for three needles per tree. Mortality was recorded each year from 2007 to 2010.Fig. 3Phenological stages of bud burst in Pinus sylvestris assessed in field trials. Inset numbers correspond to budburst stage. Budburst stage 1: bud dormant; 2: bud swelling and showing signs of linear expansion; 3: scales open at base revealing green tissue. Remaining bud remains encased by smooth bud scales; 4: scales open along length of shoot revealing green tissue and partially visible needles; 5: white tipped needles visible along length of the shoot; 6: green needles emerging away from the shoot (bottle brush appearance) along its entire length; 7: Needles have separated and next year’s terminal bud is usually formed and clearly visible.Field traitsTree height was measured in the field in the winter after each growing season from 2013 at FE and FW, and from 2014 to 2020 at all sites. Height was taken as the vertical measurement in cm from top bud straight to the ground. Basal stem diameter was measured at the end of the growing season for trees growing at FE and FW from 2014 to 2020 and for FS in 2020.Phenology assessments were performed in spring at each site from 2015 to 2019. Seven distinct stages of budburst (assessed on the terminal bud) were defined (Fig. 3) although only stages 4 to 6 are included in the dataset and considered for analysis due to high proportions of missing data for the early and late stages. Each tree was assessed for budburst stage at weekly intervals from early spring until budburst was complete. In order to allow comparisons within and among sites and years, the date at which each stage of budburst occurred was considered relative to 31 March of that year. For example, 25 May 2019 is recorded as 55 days since 31 March 2019. The duration of budburst (time taken to reach stage 6 from stage 4) was also estimated.When trees progressed through budburst stages rapidly, skipping a stage between assessments, a mean value was taken from the two assessment dates. For example, if a tree was at stage 4 on day 55 and was recorded as stage 6 at the next assessment on day 62, it is assumed to have reached stage 5 at day 58.5.Data RecordsData are deposited with the Environmental Information Data Centre (http:/eidc.ac.uk)13–15.In all datasets, the first two columns are:1: Population code (code for forest of origin, 21 total. See Table 1 for full population names)2: Family (unique mother tree code: progeny described in nursery traits dataset and field trait dataset from the same family are putative half-siblings)There are two maternal traits datasets: one for traits relating directly to the mother tree (MotherTraits.txt) in which each row represents one tree, and a second for traits relating to cones and seed collected from each mother tree in which each row represents one cone collected from each tree. Columns are defined as follows:Maternal dataset15: MotherTraits.txt3: Population (name of forest of origin, 21 total)4: Seed zone5–13. Location reference and immediate environment for each tree: Latitude [decimal degree]; Longitude [decimal degree]; Aspect; Slope; Altitude [m]; Regeneration [1–4]; soil depth [cm]; soil moisture [1–5]; mean distance to nearest three trees [M]14–15. Mother tree traits (absolute height [M]; diameter at breast height [M])Maternal dataset15: ConeSeedTraits.txt3. Cone number (1–10)4–6. Cone traits (Width [mm], Length [mm], Weight [g]). Missing data = NA7–9. Seed traits (Viable seeds [count], Percentage viable seeds (%), Viable seeds weight [g]). Missing data = NANursery traits dataset13:3. Seedling id4. Nursery site [NE; NW; NG]5–6. Nursery block code: from 2007 to 2010 [1–40]; from 2010 to 2012 [1–98]7. Field code if transplanted [4 digit code; not transplanted = NA]8–40. Nursery traits (see Table 3 for list of traits and format of column header). Missing data = NA41–44. Status [Alive; Dead] for each year 2007–2010 (inclusive)Field traits dataset14:3. Field site [FE; FS; FW]4. Field code [4 digit code]5. Block code [A; B; C; D]6–7. Position of every tree in each site: rows [FE: 1–21; FS: 1–9; FW: 1–21] and columns [FE: 1–32; FS: 1–86; FW: 1–24]. Note: in FS guard rows are present between blocks B and C at columns 43 and 44; in FW tree 6424 was planted between rows 19 and 20 in error, and has therefore been assigned row 19.58–58. Field traits (see Table 3 for list of traits and format of column header). Missing data = NATechnical ValidationMeasurements repeated annually are performed at the same time of year to ensure consistency of the method, e.g. height is measured during winter to avoid the possibility of active growth occurring after the trees are measured. Data were checked after each survey and inconsistent values (i.e. where height was less than the previous year) were re-measured. Where height increment was found to be less than 0 mm (due to an error in measurement or an effective loss in height due to damage) values were removed and classed as missing (NA). Annual stem diameter increment was estimated as the increase in stem diameter from the end of one growing season to the end of the next. Where stem diameter increment was estimated between 0 and −2 mm, the error was assumed to be caused by differences in orientation of measurement between years and increment values were adjusted to 0. Where stem diameter increment was found to be less than −2 mm, the increment and most recent stem diameter measurements were both classed as missing (NA).We used boxplots to visualize data range and data distribution for each trait in each year over all nursery and field sites and all populations (Fig. 4).Fig. 4Box and whisker plots of trait values for each trait in each year over all nursery/field sites and all populations. Trait value units are listed in Table 3. Years are abbreviated (e.g. 2015 = ’15). Solid black lines indicate the median trait value. The bottom and top of boxes indicate the first and third quartile. The upper and lower whiskers extend to the highest and lowest values within 1.5 times the interquartile range. Individual points indicate outliers. Traits derived from the maternal datasets are green, those from the nursery dataset are blue and those from the field dataset are red.