Data and code for: Adaptive evolution of freezing tolerance in oaks is key to their dominance in North America
Creators: name: Cavender-Bares, Jeannine
Publisher: name: Zenodo
Year: 2024
Subjects: freezing tolerance
Quercus
Evolution
Cold acclimation
growth rate
freezing tolerance growth rate trade-off
leaf habit
deciduous or evergreen
Ornstein Uhlenbeck model
phylogeny
Oaks

description: Freezing tolerance plays a pivotal role in shaping the distribution and diversification of organisms across diverse habitats and was likely crucial to the expansion and adaptive radiation of the American oaks. We investigated the dynamics of adaptation to climate and potential trade-offs between stem freezing tolerance and growth rate in 48 Quercus species from five genus sections grown under temperate and tropical conditions. Species from colder regions exhibit higher freezing tolerance, lower growth rates and higher acclimation potential compared to species from warmer climates, suggesting evolutionary adaptations to seasonal climate fluctuations. Deciduous species show marked variability in freezing tolerance across their broad climatic range while evergreen species, confined to warm climates, display low freezing tolerance. While annual growth rates are constrained in deciduous species, we found no direct trade-off between freezing tolerance and growth because deciduous species that evolved in warm latitudes lost freezing tolerance. Despite an evolutionary lag, freezing tolerance in Quercus species is closely aligned with its optimal state. The capacity to withstand cold and adapt to a wide range of thermal environments was likely critical to adaptive radiation and current dominance of the North American oaks. Data and code are provided here.
description: Funding provided by: National Science FoundationROR ID: https://ror.org/021nxhr62Award Number: 2021898
Funding provided by: National Science FoundationROR ID: https://ror.org/021nxhr62Award Number: 1146380
description: We examined 48 species of Quercus from 5 different sections: Cerris (3 species), Lobatae (19), Protobalanus (1), Quercus (26), and Virentes (3 species; Table S1). The species are characterized by different growth forms, leaf phenology (34 deciduous and 14 evergreen species), ecological niches, and distribution ranges. Acorns were collected from one to five populations within each of the 48 oak species in the fall of 2010 or 2011 to capture variation across the range of the species and to maximize the phylogenetic coverage of the oaks within North America. Acorns from the Cerris section were gathered from locations outside North America, as these species are not indigenous to the region. Seeds were collected from California, Florida, Illinois, Michigan, Minnesota, North Carolina, New York, Ohio, Utah, and Wyoming from trees occurring in natural areas within each population. The length and width of each acorn were measured, and the acorns were stored for several months at 4ºC at the University of Minnesota until synchronous planting (25 to 40 seedlings per species) in glasshouse facilities at the Franklinville Experimental Research Station in Franklinville, New York or at the University of Minnesota (UMN).                

Seeds were germinated and grown in two replicated climate treatments in each of the two greenhouse locations for two years in (1) a tropical treatment, in which daytime temperature was maintained between 30 and 35ºC, and night-time temperature between 22º and 26ºC; and (2) a temperate treatment, in which winter growth temperatures (mid-November to March) reached a minimum night-time temperature of 4ºC with a minimum daytime temperature of 15ºC, simulating species' range limits.  Plants in the temperate treatment were acclimated to the winter temperature regime approximately 3 months prior to measurement. In New York, each climate regime was replicated in three independently controlled glasshouse rooms (for a total of six rooms) and in Minnesota, each climate regime was replicated in two independently controlled rooms (total of four). Summer photoperiods were left at ambient levels with a maximum of 15.4 hours in Minnesota and 15.16 in New York. The photoperiod was extended to 12 hours during the winter. These climate treatments were repeated every year for three subsequent years.            

Three months after the start of the winter climate treatment, all individuals were monitored for maximum height (h; root collar to apex of tallest stem), the number of leaves (total leaves on all stems), and stem diameter (d; measured at the root collar). Allometric equations used to estimate the biomass for each species. Above-ground biomass was estimated using the maximum height, total number of leaves, and stem diameter. Acorn mass was estimated using the acorn volume.  Absolute growth rate (AGR; g yr-1) was calculated as AGR = (Vfinal -Vinitial) ⁄ (Tfinal - Tinitial), and Relative growth rate (RGR; mm3 mm-3 yr-1) was calculated as RGR= [loge(Vfinal) - loge(Vinitial)] ⁄ (Tfinal - Tinitial). Where Vfinal and Vinitial were calculated as the volume of a cylinder:  Vfinal/initial ­= (dfinal/initial/2)­2 x P x hfinal/initial. Tinitial is the approximate date plants were measured in March 2012 and Tfinal is the date plants.

We used the most comprehensive phylogenetic tree for the genus Quercus and performed a comparative analysis using the 'SLOUCH' package in R. SLOUCH models traits as evolving under an Ornstein-Uhlenbeck (OU) model of evolution in a maximum likelihood framework. Traits are modeled as evolving toward a primary optimum, which is influenced by a continuous, randomly changing predictor variable, in our case, minimum temperature. Estimates of the phylogenetic signal or phylogenetic inertia are given by the phylogenetic half-life (t1/2) and are interpreted relative to the absolute time associated with the phylogeny. When t1/2 is close to 0, the trait is evolving under strong selection by the predictors included in the model, without any influence from the ancestral trait value, whereas increasing values of t1/2 reflect an increasing effect of the ancestral trait value on present-day trait values.

We used index of injury measured at -15oC as the response variable of the model, while minimum temperature of the coldest month in the native climate was treated as a predictor variable modeled as evolving according to a Wiener (Brownian motion) process. We compared models with and without the predictor variable to test the role of climatic niche in freezing tolerance evolution. For all statistical analyses, we used R v.3.4.2 with base packages R core Team (2018).