Research Progress HIGHLIGHTS

Activity 1 – Climate adaptation

Activity 1.1 Candidate genes in Douglas-fir


To identify the candidate genes relating to drought and frost tolerance, we have established 1) a large outdoor common garden trial to test frost hardiness, and 2) a large indoor container experiment to examine response to drought stress.  We are growing seedlings from 73 wild collected (natural) populations and 14 orchard seedlots (selectively bred trees).  Of the 73 wild seedlots, 20 will be used for the case-control analysis for genotype-phenotype association and genotype-environment tests.  Seven of the 14 orchard seedlots are part of a separate study on assisted migration – these seven seedlots are presently growing on 45 different sites (i.e. climates) across western Canada and the US where they are being phenotyped for survival and growth.

Seedlot selection was finalized in December 2016.  The seedlots were chosen to provide adequate climatic and geographic variation across the natural range of coastal and interior Douglas-fir.  They span a range of more the 16° latitude, 14° longitude and from 61 m to 1737 m in elevation.

The seeds were stratified in April 2017, double sown into small cones at the greenhouse at UBC in May 2017 and thinned in June 2017. At the end of their first growing season, needle samples were taken from all individuals (7,200 samples) and stored for future DNA extractions.

In their second year of growth, the seedlings were divided in two experimental groups: frost and drought. In March 2018, seedlings in the frost group were planted into large, deep, outdoor raised beds, where they will be exposed to the natural seasonal patterns of daylight length and temperature. They will be measured throughout the growing season (e.g. growth, bud burst, but set timing), and tested for cold hardiness in the fall of 2018.

The seedlings in the drought group remained in the climate controlled chamber in the greenhouse. In April 2018, they were planted into multi-seedling containers in preparation for the drought tolerance experiment. They will spend 8 weeks in these containers under optimal conditions, actively growing and acclimatizing to the containers before the drought treatment begins (summer of 2018).

For the 20 case-control populations in the drought and frost experiments, we will select the 10 least and the 10 most tolerant individuals to be sequenced. Additionally, for the 87 populations, the DNA from only 40 randomly selected individuals (out of 100) per population will be pooled for sequencing.

In year 3 we will put together the phenotype data with the sequencing data and perform genotype-phenotype association analysis.

For more information, contact Rafael Candido-Ribeiro, Pia Smets, Christine Chourmouzis


Illustration of frost hardiness and drought tolerance testing.



Map of seed source locations for Douglas-fir.


Year 1. Douglas-fir germinates (May 2017) and the discards after thinning the double sowed seedlings (June 2017).


Year 1. Douglas-fir seedlings (July 2017) and tissue sampling for DNA extraction (October 2017).


Year 2. Planting Douglas-fir seedlings in the outdoor raised beds (March 2018) in preparation for frost hardiness testing in the fall of 2018.


Year 2. Planting Douglas-fir case control populations into containers (April 2018) in preparation for the drought tolerance testing in summer 2018.


Activity 1.2 Candidate genes in western larch and jack pine

The comparative genomics for Douglas-fir (Activity 1.1 above) and lodgepole pine (from the AdapTree project) will be used to increase the power for identifying climate-adaptation candidate genes in western larch and jack pine. The genomics of climate adaptation in Douglas-fir will be characterized in comparison with western larch (Pseudotsuga and Larix are sister genera) and jack pine will be compared to lodgepole pine. Additionally, there will be a four-way contrast among all species. The extent of parallel adaptation among conifer species can provide a window into the genomics of complex traits and increase the power to separate false positives from true climate-driven genetic patterns.

A genotype x environment interactions association approach will be used to reveal patterns of local adaption to climate. Sequence capture method will be applied to enrich the exonic regions for each species. Sequence capture method combines the hybridization capture and high-throughput sequencing, allowing the rapid and simultaneous description of numerous and potentially large targeted DNA regions. Capture probes are used to target and hybridize the DNA library. We can estimate SNP frequencies from the sequenced target regions.

Capture probes

For designing capture probes, we began by assembling de novo transcriptomes for jack pine and western larch. The transcriptome is the set of all RNA molecules in one cell or a population of cells. De novo transcriptome assembly is the de novo sequence assembly method of creating a transcriptome without the aid of a reference genome. We extracted RNA from foliage of one jack pine and one western larch seedling and sequenced the RNA. The assembled transcripts were then used to design a ~50 Mb sequence capture probe panel for each species.

RNA extraction

RNA extractions on the actively growing jack pine and larch were carried out in the lab at UBC in March 2017.  RNA sequencing and transciptome assembly will be conducted by the Genome Quebec Innovation Centre at McGill.

To prepare the material for extraction the needle tissue was ground in liquid nitrogen to a fine powder using a mortar and pestle. Total RNA was extracted using Spectrum™ Plant Total RNA kit (Sigma, USA) and cleaned up by the On-Column DNase I Digestion Set (Sigma, USA). The purity of the total RNA preparations were assessed with a NanoDrop spectrophotometer (Thermo Scientific, USA). The concentration was determined by Qubit® 2.0 Fluorometer (Invitrogen, USA). The integrity and size distribution of purified RNAs were evaluated with a 2100 Bioanalyzer (Agilent Technologies, USA) equipped with an RNA Nano chip.   For more information on the RNA extraction procedures contact Dragana Obreht Vidakovic.

RNA extraction for jack pine and western larch in the lab at UBC (March 2017).


Western larch shipment arriving at UBC (December 2017) and actively growing in the climate chamber (January 2017).


Jack pine breaking bud (February 2017) and actively growing in the climate chamber (March 2017).

Sequencing libraries

For constructing sequencing libraries, we sampled foliage/bark from existing range-wide western larch and jack pine provenance trials and seed collections. For each species, we selected 40 populations, with 20 individuals, that span the full geographic and climatic range of each species. The DNA is presently being extracted from each individual with representation equalized in pools of all individuals for each population. In spring 2018, we will construct capture libraries for each population by hybridizing the probes to the pooled DNA libraries. After the capture libraries are sequenced, we will estimate SNP frequencies for each population from reads in sequenced population pools and apply the analytical approaches developed in the previous project for genotype-environment association.

We will use sequence capture probes (~50 Mb sequence panels) to capture genomic regions of interest in jack pine and western larch.


Collecting western larch needles in the Lamb Creek BC provenance trial (June 2017).


Bark sampling in jack pine provenance trials (Sept 2016 and July 2017).

For more information, contact Mengmeng Lu and Sam Yeaman.

Activity 1.4 Comparative transciptomics of Douglas-fir

How do trees adapt to different climates? In this activity, we are comparing the physiological responses of four Douglas-fir provenances from different home environments (two coastal and two interior) to simulated future dryer and warmer climates. We will zoom in on gene expression (RNA, the transcriptome) to identify climate-responsive genes. The physiological mechanisms being studied include photo-protection, energy dissipation, water-use efficiency, and photosynthetic gas exchange. This study will result in a more complete causal description of abiotic stress response in Douglas-fir.

We have grown two batches of Douglas-fir seedlings. The first batch has gone through an accelerated growing cycle and are now being used to study the effects of warm autumn to the development of cold hardiness (started in December 2017). The second batch of seedlings are now in their second growing season and will be used to study the response to heatwave (set for March 2018).  Using transcriptome data and detailed phenotyping of Douglas-fir responses to extreme climates, we will select a few climate-responsive genes. These candidate genes will be functionally validated using model plants.

Warm AutumnEffects of a warmer climate on late season performance and the development of cold hardiness.

Douglas-fir seedlings exposed to a short photoperiod (needed to initiate budset) and summer temperature will be compared to seedlings exposed to normal autumn temperature. The difference in photosynthesis and chlorophyll fluorescence will be measured using a portable gas exchange fluorescence system and several pigments will be quantified to understand photoinhibition in these plants. Both sets of plants will be exposed to below freezing conditions to determine their level of freezing tolerance.  The phenotypic data will be complemented by RNA-seq quantification of gene expression in needles, buds and roots.

HeatwaveDifferences in the response of Douglas-fir provenances to heatwave.

In this experiment, we will subject Douglas-fir plants to seven days of no water supply, followed by up to two weeks of extremely hot days. We will compare the responses of heatwave treated Douglas-fir to plants grown under extreme hot and normal summer conditions. Water potential and soil water will be monitored regularly and photosynthesis and photoinhibition will be measured weekly during the treatment. Samples will be collected for quantification of pigments, carbon stable isotope and gene expression profiling.

Douglas fir seedlings growing in the research greenhouse at University of Toronto Mississauga. These seedlings are in their second growing season and will be used to study response to heatwave in March 2018.

For more information, contact Vera Velasco and Ingo Ensminger.