Nutri-cereal tissue-specific transcriptome atlas during development: Functional integration of gene expression to identify mineral uptake pathways in little millet (Panicum sumatrense).

Little millet (Panicum sumatrense Roth ex Roem. & Schult.) is an essential minor millet of southeast Asia and Africa's temperate and subtropical regions. The plant is stress-tolerant, has a short life cycle, and has a mineral-rich nutritional profile associated with unique health benefits. We report the developmental gene expression atlas of little millet (genotype JK-8) from ten tissues representing different stages of its life cycle, starting from seed germination and vegetative growth to panicle maturation. The developmental transcriptome atlas led to the identification of 342 827 transcripts. The BUSCO analysis and comparison with the transcriptomes of related species confirm that this study presents high-quality, in-depth coverage of the little millet transcriptome. In addition, the eFP browser generated here has a user-friendly interface, allowing interactive visualizations of tissue-specific gene expression. Using these data, we identified transcripts, the orthologs of which in Arabidopsis and rice are involved in nutrient acquisition, transport, and response pathways. The comparative analysis of the expression levels of these transcripts holds great potential for enhancing the mineral content in crops, particularly zinc and iron, to address the issue of "hidden hunger" and to attain nutritional security, making it a valuable asset for translational research.

Multi-omics atlas of combinatorial abiotic stress responses in wheat.

Field-grown crops rarely experience growth conditions in which yield can be maximized. Environmental stresses occur in combination, with advancements in crop tolerance further complicated by its polygenic nature. Strategic targeting of causal genes is required to meet future crop production needs. Here, we employed a systems biology approach in wheat (Triticum aestivum L.) to investigate physio-metabolic adjustments and transcriptome reprogramming involved in acclimations to heat, drought, salinity and all combinations therein. A significant shift in magnitude and complexity of plant response was evident across stress scenarios based on the agronomic losses, increased proline concentrations and 8.7-fold increase in unique differentially expressed transcripts (DETs) observed under the triple stress condition. Transcriptome data from all stress treatments were assembled into an online, open access eFP browser for visualizing gene expression during abiotic stress. Weighted gene co-expression network analysis revealed 152 hub genes of which 32% contained the ethylene-responsive element binding factor-associated amphiphilic repression (EAR) transcriptional repression motif. Cross-referencing against the 31 DETs common to all stress treatments isolated TaWRKY33 as a leading candidate for greater plant tolerance to combinatorial stresses. Integration of our findings with available literature on gene functional characterization allowed us to further suggest flexible gene combinations for future adaptive gene stacking in wheat. Our approach demonstrates the strength of robust multi-omics-based data resources for gene discovery in complex environmental conditions. Accessibility of such datasets will promote cross-validation of candidate genes across studies and aid in accelerating causal gene validation for crop resiliency.

Drought-induced fiber water release and xylem embolism susceptibility of intact balsam poplar saplings.

Balsam poplar (Populus balsamifera L.) is a widespread tree species in North America with significant ecological and economic value. However, little is known about the susceptibility of saplings to drought-induced embolism and its link to water release from surrounding xylem fibers. Questions remain regarding localized mechanisms that contribute to the survival of saplings in vivo of this species under drought. Using X-ray micro-computed tomography on intact saplings of genotypes Gillam-5 and Carnduff-9, we found that functional vessels are embedded in a matrix of water-filled fibers under well-watered conditions in both genotypes. However, water-depleted fibers started to appear under moderate drought stress while vessels remained water-filled in both genotypes. Drought-induced xylem embolism susceptibility was comparable between genotypes, and a greater frequency of smaller diameter vessels in GIL-5 did not increase embolism resistance in this genotype. Despite having smaller vessels and a total vessel number that was comparable to CAR-9, stomatal conductance was generally higher in GIL-5 compared to CAR-9. In conclusion, our in vivo data on intact saplings indicate that differences in embolism susceptibility are negligible between GIL-5 and CAR-9, and that fiber water release should be considered as a mechanism that contributes to the maintenance of vessel functional status in saplings of balsam poplar experiencing their first drought event.

Establishment of willows using the novel DeValix technique: ecological restoration mats designed for phytotechnologies.

Successful willow (Salix spp., hybrids and cultivars) establishment is a major determinant of their effectiveness when grown for phytotechnologies. Vertically-planted hardwood cuttings have been shown to produce adequate willow growth and survival, although site conditions at phytoremediation installations can make vertical planting methods unsuitable. The DeValix willow mat restoration technique was designed and tested as an alternative horizontal planting method that can be installed by hand in a variety of environmental applications. The DeValix technique was evaluated by testing five willow clones ("Millbrook"; "Sherburne"; "SX61"; "SX67"; "Tully Champion") grown at two phytoremediation sites (Ontonagon, MI; Manitowoc, WI) for the 2019 growing season. Differences in survival and growth were tested among sites, genotypes, and their interactions. Stem height, diameter, and number of stems per mat were compared to identify clones with greater establishment success and higher phytoremediation potential. Results demonstrated significant effects of site (p < 0.0001) and clone (p < 0.0001) on shoot number. Additionally, the site × clone interaction significantly affected stem height (p = 0.0045) and diameter (p = 0.0166). Stem density ranged from 95,000 to 212,000 stems per hectare, indicating the DeValix technique is a viable establishment method for environmental applications, including phytoremediation and shoreline stabilization.

This research analyzes the establishment success of the DeValix technique, a novel horizontal planting method for willow cultivars, and evaluates the DeValix technique as an alternative to other horizontal techniques and traditional vertical planting methods currently used in environmental applications. Results from this study add to the current knowledge of planting techniques and assesses the use of the DeValix technique for planting willow cultivars that are currently being tested in several phytotechnology systems.

A meta-analysis of mesophyll conductance to CO2 in relation to major abiotic stresses in poplar species.

Mesophyll conductance (gm) determines the diffusion of CO2 from the substomatal cavities to the site of carboxylation in the chloroplasts and represents a critical component of the diffusive limitation of photosynthesis. In this study, we evaluated the average effect sizes of different environmental constraints on gm in Populus spp., a forest tree model. We collected raw data of 815 A-Ci response curves from 26 datasets to estimate gm, using a single curve-fitting method to alleviate method-related bias. We performed a meta-analysis to assess the effects of different abiotic stresses on gm. We found a significant increase in gm from the bottom to the top of the canopy that was concomitant with the increase of maximum rate of carboxylation and light-saturated photosynthetic rate (Amax). gm was positively associated with increases in soil moisture and nutrient availability, but was insensitive to increasing soil copper concentration and did not vary with atmospheric CO2 concentration. Our results showed that gm was strongly related to Amax and to a lesser extent to stomatal conductance (gs). Moreover, a negative exponential relationship was obtained between gm and specific leaf area, which may be used to scale-up gm within the canopy.

Dissecting the Roles of Cuticular Wax in Plant Resistance to Shoot Dehydration and Low-Temperature Stress in Arabidopsis.

Cuticular waxes are a mixture of hydrophobic very-long-chain fatty acids and their derivatives accumulated in the plant cuticle. Most studies define the role of cuticular wax largely based on reducing nonstomatal water loss. The present study investigated the role of cuticular wax in reducing both low-temperature and dehydration stress in plants using Arabidopsis thaliana mutants and transgenic genotypes altered in the formation of cuticular wax. cer3-6, a known Arabidopsis wax-deficient mutant (with distinct reduction in aldehydes, n-alkanes, secondary n-alcohols, and ketones compared to wild type (WT)), was most sensitive to water loss, while dewax, a known wax overproducer (greater alkanes and ketones compared to WT), was more resistant to dehydration compared to WT. Furthermore, cold-acclimated cer3-6 froze at warmer temperatures, while cold-acclimated dewax displayed freezing exotherms at colder temperatures compared to WT. Gas Chromatography-Mass Spectroscopy (GC-MS) analysis identified a characteristic decrease in the accumulation of certain waxes (e.g., alkanes, alcohols) in Arabidopsis cuticles under cold acclimation, which was additionally reduced in cer3-6. Conversely, the dewax mutant showed a greater ability to accumulate waxes under cold acclimation. Fourier Transform Infrared Spectroscopy (FTIR) also supported observations in cuticular wax deposition under cold acclimation. Our data indicate cuticular alkane waxes along with alcohols and fatty acids can facilitate avoidance of both ice formation and leaf water loss under dehydration stress and are promising genetic targets of interest.

Learning from methylomes: epigenomic correlates of Populus balsamifera traits based on deep learning models of natural DNA methylation.

Epigenomes have remarkable potential for the estimation of plant traits. This study tested the hypothesis that natural variation in DNA methylation can be used to estimate industrially important traits in a genetically diverse population of Populus balsamifera L. (balsam poplar) trees grown at two common garden sites. Statistical learning experiments enabled by deep learning models revealed that plant traits in novel genotypes can be modelled transparently using small numbers of methylated DNA predictors. Using this approach, tissue type, a nonheritable attribute, from which DNA methylomes were derived was assigned, and provenance, a purely heritable trait and an element of population structure, was determined. Significant proportions of phenotypic variance in quantitative wood traits, including total biomass (57.5%), wood density (40.9%), soluble lignin (25.3%) and cell wall carbohydrate (mannose: 44.8%) contents, were also explained from natural variation in DNA methylation. Modelling plant traits using DNA methylation can capture tissue-specific epigenetic mechanisms underlying plant phenotypes in natural environments. DNA methylation-based models offer new insight into natural epigenetic influence on plants and can be used as a strategy to validate the identity, provenance or quality of agroforestry products.

Latitudinal clines in bud flush phenology reflect genetic variation in chilling requirements in balsam poplar, Populus balsamifera.

PREMISE: Boreal and northern temperate forest trees possess finely tuned mechanisms of dormancy, which match bud phenology with local seasonality. After winter dormancy, the accumulation of chilling degree days (CDD) required for rest completion before the accumulation of growing degree days (GDD) during quiescence is an important step in the transition to spring bud flush. While bud flush timing is known to be genetically variable within species, few studies have investigated variation among genotypes from different climates in response to variable chilling duration. METHODS: We performed a controlled environment study using dormant cuttings from 10 genotypes of Populus balsamifera, representing a broad latitudinal gradient (43-58°N). We exposed cuttings to varying amounts of chilling (0-10 weeks) and monitored subsequent GDD to bud flush at a constant forcing temperature. RESULTS: Chilling duration strongly accelerated bud flush timing, with increasing CDD resulting in fewer GDD to flush. Genotypic variation for bud flush was significant and stratified by latitude, with southern genotypes requiring more GDD to flush than northern genotypes. The latitudinal cline was pronounced under minimal chilling, whereas genotypic variation in GDD to bud flush converged as CDD increased. CONCLUSIONS: We demonstrate that increased chilling lessens GDD to bud flush in a genotype-specific manner. Our results emphasize that latitudinal clines in bud flush reflect a critical genotype-by-environment interaction, whereby differences in bud flush between southern vs. northern genotypes depend on chilling. Our results suggest selection has shaped chilling requirements and depth of rest as an adaptive strategy to avoid precocious flush in climates with midwinter warming.

A Cost-Effective and Portable Optical Sensor System to Estimate Leaf Nitrogen and Water Contents in Crops.

Non-invasive determination of leaf nitrogen (N) and water contents is essential for ensuring the healthy growth of the plants. However, most of the existing methods to measure them are expensive. In this paper, a low-cost, portable multispectral sensor system is proposed to determine N and water contents in the leaves, non-invasively. Four different species of plants-canola, corn, soybean, and wheat-are used as test plants to investigate the utility of the proposed device. The sensor system comprises two multispectral sensors, visible (VIS) and near-infrared (NIR), detecting reflectance at 12 wavelengths (six from each sensor). Two separate experiments were performed in a controlled greenhouse environment, including N and water experiments. Spectral data were collected from 307 leaves (121 for N and 186 for water experiment), and the rational quadratic Gaussian process regression (GPR) algorithm was applied to correlate the reflectance data with actual N and water content. By performing five-fold cross-validation, the N estimation showed a coefficient of determination () of 63.91% for canola, 80.05% for corn, 82.29% for soybean, and 63.21% for wheat. For water content estimation, canola showed an of 18.02%, corn showed an of 68.41%, soybean showed an of 46.38%, and wheat showed an of 64.58%. The result reveals that the proposed low-cost sensor with an appropriate regression model can be used to determine N content. However, further investigation is needed to improve the water estimation results using the proposed device.

Biotechnological strategies for improved photosynthesis in a future of elevated atmospheric CO2.

MAIN CONCLUSION: The improvement of photosynthesis using biotechnological approaches has been the focus of much research. It is now vital that these strategies be assessed under future atmospheric conditions. The demand for crop products is expanding at an alarming rate due to population growth, enhanced affluence, increased per capita calorie consumption, and an escalating need for plant-based bioproducts. While solving this issue will undoubtedly involve a multifaceted approach, improving crop productivity will almost certainly provide one piece of the puzzle. The improvement of photosynthetic efficiency has been a long-standing goal of plant biotechnologists as possibly one of the last remaining means of achieving higher yielding crops. However, the vast majority of these studies have not taken into consideration possible outcomes when these plants are grown long-term under the elevated CO2 concentrations (e[CO2]) that will be evident in the not too distant future. Due to the considerable effect that CO2 levels have on the photosynthetic process, these assessments should become commonplace as a means of ensuring that research in this field focuses on the most effective approaches for our future climate scenarios. In this review, we discuss the main biotechnological research strategies that are currently underway with the aim of improving photosynthetic efficiency and biomass production/yields in the context of a future of e[CO2], as well as alternative approaches that may provide further photosynthetic benefits under these conditions.

Gene copy number variations involved in balsam poplar (Populus balsamifera L.) adaptive variations.

Gene copy number variations (CNVs) involved in phenotypic variations have already been shown in plants, but genomewide testing of CNVs for adaptive variation was not doable until recent technological developments. Thus, reports of the genomic architecture of adaptation involving CNVs remain scarce to date. Here, we investigated F1 progenies of an intraprovenance cross (north-north cross, 58th parallel) and an interprovenances cross (north-south cross, 58th/49th parallels) for CNVs using comparative genomic hybridization on arrays of probes targeting gene sequences in balsam poplar (Populus balsamifera L.), a widespread North American forest tree. A total of 1,721 genes were found in varying copy numbers over the set of 19,823 tested genes. These gene CNVs presented an estimated average size of 8.3 kb and were distributed over poplar's 19 chromosomes including 22 hotspot regions. Gene CNVs number was higher for the interprovenance progeny in accordance with an expected higher genetic diversity related to the composite origin of this family. Regression analyses between gene CNVs and seven adaptive trait variations resulted in 23 significant links; among these adaptive gene CNVs, 30% were located in hotspots. One-to-five gene CNVs were found related to each of the measured adaptive traits and annotated for both biotic and abiotic stress responses. These annotations can be related to the occurrence of a higher pathogenic pressure in the southern parts of balsam poplar's distribution, and higher photosynthetic assimilation rates and water-use efficiency at high latitudes. Overall, our findings suggest that gene CNVs typically having higher mutation rates than SNPs may in fact represent efficient adaptive variations against fast-evolving pathogens.

The adaptive potential of Populus balsamifera L. to phenology requirements in a warmer global climate.

The manner in which organisms adapt to climate change informs a broader understanding of the evolution of biodiversity as well as conservation and mitigation plans. We apply common garden and association mapping approaches to quantify genetic variance and identify loci affecting bud flush and bud set, traits that define a tree's season for height growth, in the boreal forest tree Populus balsamifera L. (balsam poplar). Using data from 478 genotypes grown in each of two common gardens, one near the southern edge and another near the northern edge of P. balsamifera's range, we found that broad-sense heritability for bud flush and bud set was generally high (H(2) > 0.5 in most cases), suggesting that abundant genetic variation exists for phenological response to changes in the length of the growing season. To identify the molecular genetic basis of this variation, we genotyped trees for 346 candidate single nucleotide polymorphisms (SNPs) from 27 candidate genes for the CO/FT pathway in poplar. Mixed-model analyses of variance identified SNPs in 10 genes to be associated with variation in either bud flush or bud set. Multiple SNPs within FRIGIDA were associated with bud flush, whereas multiple SNPs in LEAFY and GIGANTEA 5 were associated with bud set. Although there was strong population structure in stem phenology, the geographic distribution of multilocus association SNP genotypes was widespread except at the most northern populations, indicating that geographic regions may harbour sufficient diversity in functional genes to facilitate adaption to future climatic conditions in many sites.

Diffusional conductances to CO2 as a target for increasing photosynthesis and photosynthetic water-use efficiency.

A key objective for sustainable agriculture and forestry is to breed plants with both high carbon gain and water-use efficiency (WUE). At the level of leaf physiology, this implies increasing net photosynthesis (A N) relative to stomatal conductance (g s). Here, we review evidence for CO2 diffusional constraints on photosynthesis and WUE. Analyzing past observations for an extensive pool of crop and wild plant species that vary widely in mesophyll conductance to CO2 (g m), g s, and foliage A N, it was shown that both g s and g m limit A N, although the relative importance of each of the two conductances depends on species and conditions. Based on Fick's law of diffusion, intrinsic WUE (the ratio A N/g s) should correlate on the ratio g m/g s, and not g m itself. Such a correlation is indeed often observed in the data. However, since besides diffusion A N also depends on photosynthetic capacity (i.e., V c,max), this relationship is not always sustained. It was shown that only in a very few cases, genotype selection has resulted in simultaneous increases of both A N and WUE. In fact, such a response has never been observed in genetically modified plants specifically engineered for either reduced g s or enhanced g m. Although increasing g m alone would result in increasing photosynthesis, and potentially increasing WUE, in practice, higher WUE seems to be only achieved when there are no parallel changes in g s. We conclude that for simultaneous improvement of A N and WUE, genetic manipulation of g m should avoid parallel changes in g s, and we suggest that the appropriate trait for selection for enhanced WUE is increased g m/g s.

Timing of photoperiodic competency causes phenological mismatch in balsam poplar (Populus balsamifera L.).

Plant phenology is expected to be sensitive to climate warming. In boreal trees, spring flush is primarily temperature driven, whereas height growth cessation and autumn leaf senescence are predominantly controlled by photoperiod. Cuttings of 525 genotypes from the full range of balsam poplar were planted into two common gardens (Vancouver and Indian Head, Canada) at similar latitudes, but with differing winter temperatures and growing seasons. There was clinal variation in spring and, particularly, summer and fall phenology. Bud flush and, despite milder climate, bud set and leaf drop were earlier at Vancouver than at Indian Head by 44, 28 and 7 d, respectively. Although newly flushed growth is insensitive to photoperiod, many genotypes at both sites became competent before the summer solstice. At Vancouver, high-latitude genotypes set dormant terminal buds in mid-spring. Most other genotypes grew until midsummer or set bud temporarily and then experienced a second flush. In both gardens and in a growth chamber experiment, earlier bud set was associated with reduced height growth and higher root/shoot ratios. Shoots attained competency ~5 weeks after flushing, which would normally prevent dormancy induction before the solstice, but may be insufficient if spring advances by more than a few weeks.

Comparative nucleotide diversity across North American and European populus species.

Nucleotide polymorphisms in two North American balsam poplars (Populus trichocarpa Torr. & Gray and P. balsamifera L.; section Tacamahaca), and one Eurasian aspen (P. tremula L.; section Populus) were compared using nine loci involved in defense, stress response, photoperiodism, freezing tolerance, and housekeeping. Nucleotide diversity varied among species and was highest for P. tremula (θ(w) = 0.005, π(T) = 0.007) as compared to P. balsamifera (θ(w) = 0.004, π(T) = 0.005) or P. trichocarpa (θ(w) = 0.002, π(T) = 0.003). Across species, the defense and the stress response loci accounted for the majority of the observed level of nucleotide diversity. In general, the studied loci did not deviate from neutral expectation either at the individual locus (non-significant normalized Fay and Wu's H) or at the multi-locus level (non-significant HKA test). Using molecular clock analysis, section Tacamahaca probably shared a common ancestor with section Populus approximately 4.5 million year ago. Divergence between the two closely related balsam poplars was about 0.8 million years ago, a pattern consistent with an isolation-with-migration (IM) model. As expected, P. tremula showed a five-fold higher substitution rate (2 × 10(-8) substitution/site/year) compared to the North American species (0.4 × 10(-8) substitution/site/year), probably reflecting its complex demographic history. Linkage disequilibrium (LD) varied among species with a more rapid decay in the North American species (<400 bp) in comparison to P. tremula (≫400 bp). The similarities in nucleotide diversity pattern and LD decay of the two balsam poplar species likely reflects the recent time of their divergence.

Enlightening the Pathway of Phytoremediation: Ecophysiology and X-ray Fluorescence Visualization of Two Chilean Hardwoods Exposed to Excess Copper.

In the present climate emergency due to global warming, we are urged to move away from fossil fuels and pursue a speedy conversion to renewable energy systems. Consequently, copper (Cu) will remain in high demand because it is a highly efficient conductor used in clean energy systems to generate power from solar, hydro, thermal and wind energy across the world. Chile is the global leader in copper production, but this position has resulted in Chile having several hundred tailing deposits. We grew two Chilean native hardwood species, quillay (Quillaja saponaria Molina) and espino (Vachellia caven (Molina) Seigler & Ebinger, under three increasing Cu levels (0, 50, and 100 µM) for 6 months in a greenhouse setting. We measured growth, photosynthetic performance and elemental contents of leaves and roots to further evaluate their potential for phytoremediation. Growth of quillay was unaffected by Cu treatment but growth of espino was enhanced, as was its photosynthetic performance, indicating that espino may have an unusually high requirement for copper. Excess Cu was mostly restricted to the roots of both species, where X-ray fluorescence (XRF) mapping indicated some tendency for Cu to accumulate in tissues outside the periderm. Calcium oxalate crystals were prominently visible in XRF images of both species. Nickel (but not Cu) showed a concurrent distribution pattern with these crystals.

Nitrogen isotope discrimination in open-pollinated and hybrid canola suggests indirect selection for enhanced ammonium utilization.

Nitrogen isotope discrimination (Δ15N) may have utility as an indicator of nitrogen use in plants. A simple Δ15N-based isotope mass balance (IMB) model has been proposed to provide estimates of efflux/influx (E/I) ratios across root plasma membranes, the proportion of inorganic nitrogen assimilation in roots (P root) and translocation of inorganic nitrogen to shoots (Ti/Tt) under steady-state conditions. We used the IMB model to investigate whether direct selection for yield in canola (Brassica napus L.) has resulted in indirect selection in traits related to nitrogen use. We selected 23 canola lines developed from 1942 to 2017, including open-pollinated (OP) lines developed prior to 2005 as well as more recent commercial hybrids (CH), and in three separate experiments grew them under hydroponic conditions in a greenhouse with either 0.5 mM ammonium, 0.5 mM nitrate, or 5 mM nitrate. Across all lines, E/I, Proot and Ti/Tt averaged 0.09±0.03, 0.82±0.05 and 0.23±0.06 in the low nitrate experiment, and 0.31±0.06, 0.71±0.07 and 0.42±0.12 in the high nitrate experiment, respectively. In contrast, in the ammonium experiment average E/I was 0.40±0.05 while Ti/Tt averaged 0.07±0.04 and Proot averaged 0.97±0.02. Although there were few consistent differences between OP and CH under nitrate nutrition, commercial hybrids were collectively better able to utilize ammonium as their sole nitrogen source, demonstrating significantly greater overall biomass and a lower Proot and a higher Ti/Tt, suggesting a somewhat greater flux of ammonium to the shoot. Average root and whole-plant Δ15N were also slightly higher in CH lines, suggesting a small increase in E/I. An increased ability to tolerate and/or utilize ammonium in modern canola hybrids may have arisen under intensive mono-cropping.

Genotypic variation in C and N isotope discrimination suggests local adaptation of heart-leaved willow.

Plants acquire multiple resources from the environment and may need to adjust and/or balance their respective resource-use efficiencies to maximize grow and survival, in a locally adaptive manner. In this study, tissue and whole-plant carbon (C) isotopic composition (δ13C) and carbon/nitrogen (C/N) ratios provided long-term measures of use efficiencies for water (WUE) and nitrogen (NUE), and a nitrogen (N) isotopic composition (δ15N)-based mass balance model was used to estimate traits related to N uptake and assimilation in heart-leaved willow (Salix eriocephala Michx.). In an initial common garden experiment consisting of 34 populations, we found population-level variation in δ13C, C/N ratio and δ15N, indicating different patterns in WUE, NUE and N uptake and assimilation. Although there was no relationship between foliar δ13C and C/N ratios among populations, there was a significant negative correlation between these measures across all individuals, implying a genetic and/or plastic trade-off between WUE and NUE not associated with local adaptation. To eliminate any environmental effect, we grew a subset of 21 genotypes hydroponically with nitrate as the sole N source and detected significant variation in δ13C, δ15N and C/N ratios. Variation in δ15N was mainly due to genotypic differences in the nitrate efflux/influx ratio (E/I) at the root. Both experiments suggested clinal variation in δ15N (and thus N uptake efficiency) with latitude of origin, which may relate to water availability and could contribute to global patterns in ecosystem δ15N. There was a tendency for genotypes with higher WUE to come from more water-replete sites with shorter and cooler growing seasons. We found that δ13C, C/N ratio and E/I were not inter-correlated, suggesting that the selection of growth, WUE, NUE and N uptake efficiency can occur without trade-off.

Unweaving the population structure and genetic diversity of Canadian shrub willow.

Perennial shrub willow are increasingly being promoted in short-rotation coppice systems as biomass feedstocks, for phytoremediation applications, and for the diverse ecosystem services that can accrue. This renewed interest has led to widespread willow cultivation, particularly of non-native varieties. However, Canadian willow species have not been widely adopted and their inherent diversity has not yet been thoroughly investigated. In this study, 324 genotypes of Salix famelica and Salix eriocephala collected from 33 sites of origin were analyzed using 26,016 single nucleotide polymorphisms to reveal patterns of population structure and genetic diversity. Analyses by Bayesian methods and principal component analysis detected five main clusters that appeared to be largely shaped by geoclimatic variables including mean annual precipitation and the number of frost-free days. The overall observed (HO) and expected (HE) heterozygosity were 0.126 and 0.179, respectively. An analysis of molecular variance revealed that the highest genetic variation occurred within genotypes (69%), while 8% of the variation existed among clusters and 23% between genotypes within clusters. These findings provide new insights into the extent of genetic variation that exists within native shrub willow species which could be leveraged in pan-Canadian willow breeding programs.

Climate-driven local adaptation of ecophysiology and phenology in balsam poplar, Populus balsamifera L. (Salicaceae).

PREMISE OF THE STUDY: During past episodes of climate change, many plant species experienced large-scale range expansions. Expanding populations likely encountered strong selection as they colonized new environments. In this study we examine the extent to which populations of the widespread forest tree Populus balsamifera L. have become locally adapted as the species expanded into its current range since the last glaciation. METHODS: We tested for adaptive variation in 13 ecophysiology and phenology traits on clonally propagated genotypes originating from a range-wide sample of 20 subpopulations. The hypothesis of local adaption was tested by comparing among-population variation at ecologically important traits (Q(ST)) to expected variation based on demographic history (F(ST)) estimated from a large set of nuclear single nucleotide polymorphism loci. KEY RESULTS: Evidence for divergence in excess of neutral expectations was present for eight of 13 traits. Bud phenology, petiole length, and leaf nitrogen showed the greatest divergence (all Q(ST) > 0.6), whereas traits related to leaf water usage showed the least (all Q(ST) ≤ 0.30) and were not different from neutrality. Strong correlations were present between traits, geography, and climate, and they revealed a general pattern of northern subpopulations adapted to shorter, drier growing seasons compared with populations in the center or eastern regions of the range. CONCLUSIONS: Our study demonstrates pronounced adaptive variation in ecophysiology and phenology among balsam poplar populations. These results suggest that as this widespread forest tree species expanded its range since the end of the last glacial maximum, it evolved rapidly in response to geographically variable selection.