An atlas of the Norway spruce needle seasonal transcriptome.

Boreal conifers possess a tremendous ability to survive and remain evergreen during harsh winter conditions and resume growth during summer. This is enabled by coordinated regulation of major cellular functions at the level of gene expression, metabolism, and physiology. Here we present a comprehensive characterization of the annual changes in the global transcriptome of Norway spruce (Picea abies) needles as a resource to understand needle development and acclimation processes throughout the year. In young, growing needles (May 15 until June 30), cell walls, organelles, etc., were formed, and this developmental program heavily influenced the transcriptome, explained by over-represented Gene Ontology (GO) categories. Later changes in gene expression were smaller but four phases were recognized: summer (July-August), autumn (September-October), winter (November-February), and spring (March-April), where over-represented GO categories demonstrated how the needles acclimated to the various seasons. Changes in the seasonal global transcriptome profile were accompanied by differential expression of members of the major transcription factor families. We present a tentative model of how cellular activities are regulated over the year in needles of Norway spruce, which demonstrates the value of mining this dataset, accessible in ConGenIE together with advanced visualization tools.

Adaptive Introgression Facilitates Adaptation to High Latitudes in European Aspen (Populus tremula L.).

Understanding local adaptation has become a key research area given the ongoing climate challenge and the concomitant requirement to conserve genetic resources. Perennial plants, such as forest trees, are good models to study local adaptation given their wide geographic distribution, largely outcrossing mating systems, and demographic histories. We evaluated signatures of local adaptation in European aspen (Populus tremula) across Europe by means of whole-genome resequencing of a collection of 411 individual trees. We dissected admixture patterns between aspen lineages and observed a strong genomic mosaicism in Scandinavian trees, evidencing different colonization trajectories into the peninsula from Russia, Central and Western Europe. As a consequence of the secondary contacts between populations after the last glacial maximum, we detected an adaptive introgression event in a genome region of ∼500 kb in chromosome 10, harboring a large-effect locus that has previously been shown to contribute to adaptation to the short growing seasons characteristic of Northern Scandinavia. Demographic simulations and ancestry inference suggest an Eastern origin-probably Russian-of the adaptive Nordic allele which nowadays is present in a homozygous state at the north of Scandinavia. The strength of introgression and positive selection signatures in this region is a unique feature in the genome. Furthermore, we detected signals of balancing selection, shared across regional populations, that highlight the importance of standing variation as a primary source of alleles that facilitate local adaptation. Our results, therefore, emphasize the importance of migration-selection balance underlying the genetic architecture of key adaptive quantitative traits.

Stem girdling affects the onset of autumn senescence in aspen in interaction with metabolic signals.

Autumn senescence in aspen (Populus tremula) is precisely timed every year to relocate nutrients from leaves to storage organs before winter. Here we demonstrate how stem girdling, which leads to the accumulation of photosynthates in the crown, influences senescence. Girdling resulted in an early onset of senescence, but the chlorophyll degradation was slower and nitrogen more efficiently resorbed than during normal autumn senescence. Girdled stems accumulated or retained anthocyanins potentially providing photoprotection in senescing leaves. Girdling of one stem in a clonal stand sharing the same root stock did not affect senescence in the others, showing that the stems were autonomous in this respect. One girdled stem with unusually high chlorophyll and nitrogen contents maintained low carbon-to-nitrogen (C/N) ratio and did not show early senescence or depleted chlorophyll level unlike the other girdled stems suggesting that the responses depended on the genotype or its carbon and nitrogen status. Metabolite analysis highlighted that the tricarboxylic acid (TCA) cycle, salicylic acid pathway, and redox homeostasis are involved in the regulation of girdling-induced senescence. We propose that disrupted sink-source relation and C/N status can provide cues through the TCA cycle and phytohormone signaling to override the phenological control of autumn senescence in the girdled stems.

Functional and evolutionary genomic inferences in Populus through genome and population sequencing of American and European aspen.

The Populus genus is one of the major plant model systems, but genomic resources have thus far primarily been available for poplar species, and primarily Populus trichocarpa (Torr. & Gray), which was the first tree with a whole-genome assembly. To further advance evolutionary and functional genomic analyses in Populus, we produced genome assemblies and population genetics resources of two aspen species, Populus tremula L. and Populus tremuloides Michx. The two aspen species have distributions spanning the Northern Hemisphere, where they are keystone species supporting a wide variety of dependent communities and produce a diverse array of secondary metabolites. Our analyses show that the two aspens share a similar genome structure and a highly conserved gene content with P. trichocarpa but display substantially higher levels of heterozygosity. Based on population resequencing data, we observed widespread positive and negative selection acting on both coding and noncoding regions. Furthermore, patterns of genetic diversity and molecular evolution in aspen are influenced by a number of features, such as expression level, coexpression network connectivity, and regulatory variation. To maximize the community utility of these resources, we have integrated all presented data within the PopGenIE web resource (PopGenIE.org).

Storage lipid accumulation is controlled by photoperiodic signal acting via regulators of growth cessation and dormancy in hybrid aspen.

The signalling pathways that control seasonal modulation of carbon metabolism in perennial plants are poorly understood. Using genetic, metabolic and natural variation approaches, we identify factors mediating photoperiodic control of storage lipid accumulation in the model tree hybrid aspen (Populus tremula × tremuloides). We characterized lipid accumulation in transgenic hybrid aspen with impaired photoperiodic and hormonal responses. Genome-wide association mapping was performed in Swedish aspen (P. tremula) genotypes to determine genetic loci associated with genotype variation in lipid content. Our data show that the storage lipid triacylglycerol (TAG) accumulates in cambial meristem and pith rays of aspen in response to photoperiodic signal controlling growth cessation and dormancy induction. We show that photoperiodic control of TAG accumulation is mediated by the FLOWERING LOCUS T/CONSTANS module, which also controls the induction of growth cessation. Hormonal and chromatin remodelling pathways also contribute to TAG accumulation by photoperiodic signal. Natural variation exists in lipid accumulation that is controlled by input from multiple loci. Our data shed light on how the control of storage metabolism is temporally coordinated with growth cessation and dormancy by photoperiodic signal, and reveals that storage lipid accumulation between seeds and perennating organs of trees may involve distinct regulatory circuits.

Autumn senescence in aspen is not triggered by day length.

Autumn senescence in mature aspens, grown under natural conditions, is initiated at almost the same date every year. The mechanism of such precise timing is not understood but we have previously shown that the signal must be derived from light. We studied variation in bud set and autumn senescence in a collection of 116 natural Eurasian aspen (Populus tremula) genotypes, from 12 populations in Sweden and planted in one northern and one southern common garden, to test the hypothesis that onset of autumn senescence is triggered by day length. We confirmed that, although bud set seemed to be triggered by a critical photoperiod/day length, other factors may influence it. The data on initiation of autumn senescence, on the other hand, were incompatible with the trigger being the day length per se, hence the trigger must be some other light-dependent factor.

Both plant genotype and herbivory shape aspen endophyte communities.

Salicinoid phenolic glycosides are common defence substances in salicaceous trees and specialist leaf beetles use these compounds for their own defence against predators. Salicinoids vary qualitatively and qualitatively in aspen (Populus tremula) and this variation has a genetic basis. The foliar endophyte mycobiome is plentiful and we hypothesised that it is related to plant genotype, potentially mediated by salicinoid composition, and that interactions with the leaf beetle Chrysomela tremula may alter this relationship. We studied these three-way interactions in controlled greenhouse experiments. Endophytic fungi were isolated from sterilised leaf tissues with and without beetle damage, and from beetles. We confirmed that endophyte composition was influenced by host genotype. Beetle activity added generalist morphs to the mycobiome that overrode the initial host association. Yeast-like genera (Cryptococcus and Rhodotorula) were isolated only from beetle-damaged tissues and from beetles, whereas fast-growing filamentous fungi dominated beetle-free control plants. Competition experiments between filamentous fungi of plant origin and beetle-related yeasts suggested interaction of both stimulating and inhibiting modes of action amongst the fungi. As a result, we detected examples of amensalism, commensalism, parasitism and competition between the morphs tested, but we found no evidence of mutualism, and consequently no co-evolutionary relationship could be demonstrated, between yeasts carried by beetles, host genotype and associated filamentous morphs. Endophyte studies are method-dependent and high-throughput sequencing technology best define the fungal mycobiome, culturing however continues to be a cheap way to provide fundamental ecological insights and it is also required for experimental studies.

Populus tremula (European aspen) shows no evidence of sexual dimorphism.

BACKGROUND: Evolutionary theory suggests that males and females may evolve sexually dimorphic phenotypic and biochemical traits concordant with each sex having different optimal strategies of resource investment to maximise reproductive success and fitness. Such sexual dimorphism would result in sex biased gene expression patterns in non-floral organs for autosomal genes associated with the control and development of such phenotypic traits. RESULTS: We examined morphological, biochemical and herbivory traits to test for sexually dimorphic resource allocation strategies within collections of sexually mature and immature Populus tremula (European aspen) trees. In addition we profiled gene expression in mature leaves of sexually mature wild trees using whole-genome oligonucleotide microarrays and RNA-Sequencing. CONCLUSIONS: We found no evidence of sexual dimorphism or differential resource investment strategies between males and females in either sexually immature or mature trees. Similarly, single-gene differential expression and machine learning approaches revealed no evidence of large-scale sex biased gene expression. However, two significantly differentially expressed genes were identified from the RNA-Seq data, one of which is a robust diagnostic marker of sex in P. tremula.

Geographic structure in metabolome and herbivore community co-occurs with genetic structure in plant defence genes.

Plant-herbivore interactions vary across the landscape and have been hypothesised to promote local adaption in plants to the prevailing herbivore regime. Herbivores that feed on European aspen (Populus tremula) change across regional scales and selection on host defence genes may thus change at comparable scales. We have previously observed strong population differentiation in a set of inducible defence genes in Swedish P. tremula. Here, we study the geographic patterns of abundance and diversity of herbivorous insects, the untargeted metabolome of the foliage and genetic variation in a set of wound-induced genes and show that the geographic structure co-occurs in all three data sets. In response to this structure, we observe local maladaptation of herbivores, with fewer herbivores on local trees than on trees originated from more distant localities. Finally, we also identify 28 significant associations between single nucleotide polymorphisms SNPs from defence genes and a number of the herbivore traits and metabolic profiles.

Genetic markers and tree properties predicting wood biorefining potential in aspen (Populus tremula) bioenergy feedstock.

BACKGROUND: Wood represents the majority of the biomass on land and constitutes a renewable source of biofuels and other bioproducts. However, wood is recalcitrant to bioconversion, raising a need for feedstock improvement in production of, for instance, biofuels. We investigated the properties of wood that affect bioconversion, as well as the underlying genetics, to help identify superior tree feedstocks for biorefining. RESULTS: We recorded 65 wood-related and growth traits in a population of 113 natural aspen genotypes from Sweden ( https://doi.org/10.5061/dryad.gtht76hrd ). These traits included three growth and field performance traits, 20 traits for wood chemical composition, 17 traits for wood anatomy and structure, and 25 wood saccharification traits as indicators of bioconversion potential. Glucose release after saccharification with acidic pretreatment correlated positively with tree stem height and diameter and the carbohydrate content of the wood, and negatively with the content of lignin and the hemicellulose sugar units. Most of these traits displayed extensive natural variation within the aspen population and high broad-sense heritability, supporting their potential in genetic improvement of feedstocks towards improved bioconversion. Finally, a genome-wide association study (GWAS) revealed 13 genetic loci for saccharification yield (on a whole-tree-biomass basis), with six of them intersecting with associations for either height or stem diameter of the trees. CONCLUSIONS: The simple growth traits of stem height and diameter were identified as good predictors of wood saccharification yield in aspen trees. GWAS elucidated the underlying genetics, revealing putative genetic markers for bioconversion of bioenergy tree feedstocks.

Variation in non-target traits in genetically modified hybrid aspens does not exceed natural variation.

Genetically modified hybrid aspens (Populus tremula L. x P. tremuloides Michx.), selected for increased growth under controlled conditions, have been grown in highly replicated field trials to evaluate how the target trait (growth) translated to natural conditions. Moreover, the variation was compared among genotypes of ecologically important non-target traits: number of shoots, bud set, pathogen infection, amount of insect herbivory, composition of the insect herbivore community and flower bud induction. This variation was compared with the variation in a population of randomly selected natural accessions of P. tremula grown in common garden trials, to estimate how the "unintended variation" present in transgenic trees, which in the future may be commercialized, compares with natural variation. The natural variation in the traits was found to be typically significantly greater. The data suggest that when authorities evaluate the potential risks associated with a field experiment or commercial introduction of transgenic trees, risk evaluation should focus on target traits and that unintentional variation in non-target traits is of less concern.

Mechanisms of enhanced beta-adrenergic reserve from cardiac resynchronization therapy.

BACKGROUND: Cardiac resynchronization therapy (CRT) is the first clinical heart failure treatment that improves chamber systolic function in both the short-term and long-term yet also reduces mortality. The mechanical impact of CRT is immediate and well documented, yet its long-term influences on myocyte function and adrenergic modulation that may contribute to its sustained benefits are largely unknown. METHODS AND RESULTS: We used a canine model of dyssynchronous heart failure (DHF; left bundle ablation, atrial tachypacing for 6 weeks) and CRT (DHF for 3 weeks, biventricular tachypacing for subsequent 3 weeks), contrasting both to nonfailing controls. CRT restored contractile synchrony and improved systolic function compared with DHF. Myocyte sarcomere shortening and calcium transients were markedly depressed at rest and after isoproterenol stimulation in DHF (both anterior and lateral walls), and CRT substantially improved both. In addition, beta(1) and beta(2) stimulation was enhanced, coupled to increased beta(1) receptor abundance but no change in binding affinity. CRT also augmented adenylate cyclase activity over DHF. Inhibitory G-protein (Galpha(i)) suppression of beta-adrenergic stimulation was greater in DHF and reversed by CRT. Galpha(i) expression itself was unaltered; however, expression of negative regulators of Galpha(i) signaling (particularly RGS3) rose uniquely with CRT over DHF and controls. CRT blunted elevated myocardial catecholamines in DHF, restoring levels toward control. CONCLUSIONS: CRT improves rest and beta-adrenergic-stimulated myocyte function and calcium handling, upregulating beta(1) receptors and adenylate cyclase activity and suppressing G(i)-coupled signaling associated with novel RGS upregulation. The result is greater rest and sympathetic reserve despite reduced myocardial neurostimulation as components underlying its net benefit.

Inferring the Genomic Landscape of Recombination Rate Variation in European Aspen (Populus tremula).

The rate of meiotic recombination is one of the central factors determining genome-wide levels of linkage disequilibrium which has important consequences for the efficiency of natural selection and for the dissection of quantitative traits. Here we present a new, high-resolution linkage map for Populus tremula that we use to anchor approximately two thirds of the P. tremula draft genome assembly on to the expected 19 chromosomes, providing us with the first chromosome-scale assembly for P. tremula (Table 2). We then use this resource to estimate variation in recombination rates across the P. tremula genome and compare these results to recombination rates based on linkage disequilibrium in a large number of unrelated individuals. We also assess how variation in recombination rates is associated with a number of genomic features, such as gene density, repeat density and methylation levels. We find that recombination rates obtained from the two methods largely agree, although the LD-based method identifies a number of genomic regions with very high recombination rates that the map-based method fails to detect. Linkage map and LD-based estimates of recombination rates are positively correlated and show similar correlations with other genomic features, showing that both methods can accurately infer recombination rate variation across the genome. Recombination rates are positively correlated with gene density and negatively correlated with repeat density and methylation levels, suggesting that recombination is largely directed toward gene regions in P. tremula.

Leaf shape in Populus tremula is a complex, omnigenic trait.

Leaf shape is a defining feature of how we recognize and classify plant species. Although there is extensive variation in leaf shape within many species, few studies have disentangled the underlying genetic architecture. We characterized the genetic architecture of leaf shape variation in Eurasian aspen (Populus tremula L.) by performing genome-wide association study (GWAS) for physiognomy traits. To ascertain the roles of identified GWAS candidate genes within the leaf development transcriptional program, we generated RNA-Seq data that we used to perform gene co-expression network analyses from a developmental series, which is publicly available within the PlantGenIE resource. We additionally used existing gene expression measurements across the population to analyze GWAS candidate genes in the context of a population-wide co-expression network and to identify genes that were differentially expressed between groups of individuals with contrasting leaf shapes. These data were integrated with expression GWAS (eQTL) results to define a set of candidate genes associated with leaf shape variation. Our results identified no clear adaptive link to leaf shape variation and indicate that leaf shape traits are genetically complex, likely determined by numerous small-effect variations in gene expression. Genes associated with shape variation were peripheral within the population-wide co-expression network, were not highly connected within the leaf development co-expression network, and exhibited signatures of relaxed selection. As such, our results are consistent with the omnigenic model.

A single gene underlies the dynamic evolution of poplar sex determination.

Although hundreds of plant lineages have independently evolved dioecy (that is, separation of the sexes), the underlying genetic basis remains largely elusive1. Here we show that diverse poplar species carry partial duplicates of the ARABIDOPSIS RESPONSE REGULATOR 17 (ARR17) orthologue in the male-specific region of the Y chromosome. These duplicates give rise to small RNAs apparently causing male-specific DNA methylation and silencing of the ARR17 gene. CRISPR-Cas9-induced mutations demonstrate that ARR17 functions as a sex switch, triggering female development when on and male development when off. Despite repeated turnover events, including a transition from the XY system to a ZW system, the sex-specific regulation of ARR17 is conserved across the poplar genus and probably beyond. Our data reveal how a single-gene-based mechanism of dioecy can enable highly dynamic sex-linked regions and contribute to maintaining recombination and integrity of sex chromosomes.

Acetaminophen is cardioprotective against H2O2-induced injury in vivo.

Here we report our ongoing investigation of the cardiovascular effects of acetaminophen, with emphasis on oxidation-induced canine myocardial dysfunction. The objective of the current study was to investigate whether acetaminophen could attenuate exogenous H(2)O(2)-mediated myocardial dysfunction in vivo. Respiratory, metabolic, and hemodynamic indices such as left ventricular function (LVDP and +/-dP/dt(max)), and percent ectopy were measured in anesthetized, open-chest dogs during intravenous administration of 0.88 mM, 2.2 mM, 6.6 mM H(2)O(2). Following 6.6 mM H(2)O(2), tissue from the left ventricle was harvested for electron microscopy. Left ventricular function did not vary significantly between vehicle and acetaminophen groups under baseline conditions. Acetaminophen-treated dogs regained a significantly greater fraction of baseline function after high concentrations of H(2)O(2) than vehicle-treated dogs. Moreover, the incidence of H(2)O(2)-induced ventricular arrhythmias was significantly reduced in the acetaminophen-treated group. Percent ectopy following 6.6 mM concentrations of H(2)O(2) was 1 +/- 0.3 vs. 0.3 +/- 0.1 (P < 0.05) for vehicle- and acetaminophen-treated dogs, respectively. Additionally, electron micrograph images of left ventricular tissue confirmed preservation of tissue ultrastructure in acetaminophen-treated hearts when compared to vehicle. We conclude that, in the canine myocardium, acetaminophen is both functionally cardioprotective and antiarrhythmic against H(2)O(2)-induced oxidative injury.

A major locus controls local adaptation and adaptive life history variation in a perennial plant.

BACKGROUND: The initiation of growth cessation and dormancy represent critical life-history trade-offs between survival and growth and have important fitness effects in perennial plants. Such adaptive life-history traits often show strong local adaptation along environmental gradients but, despite their importance, the genetic architecture of these traits remains poorly understood. RESULTS: We integrate whole genome re-sequencing with environmental and phenotypic data from common garden experiments to investigate the genomic basis of local adaptation across a latitudinal gradient in European aspen (Populus tremula). A single genomic region containing the PtFT2 gene mediates local adaptation in the timing of bud set and explains 65% of the observed genetic variation in bud set. This locus is the likely target of a recent selective sweep that originated right before or during colonization of northern Scandinavia following the last glaciation. Field and greenhouse experiments confirm that variation in PtFT2 gene expression affects the phenotypic variation in bud set that we observe in wild natural populations. CONCLUSIONS: Our results reveal a major effect locus that determines the timing of bud set and that has facilitated rapid adaptation to shorter growing seasons and colder climates in European aspen. The discovery of a single locus explaining a substantial fraction of the variation in a key life-history trait is remarkable, given that such traits are generally considered to be highly polygenic. These findings provide a dramatic illustration of how loci of large-effect for adaptive traits can arise and be maintained over large geographical scales in natural populations.

Relative impacts of environmental variation and evolutionary history on the nestedness and modularity of tree-herbivore networks.

Nestedness and modularity are measures of ecological networks whose causative effects are little understood. We analyzed antagonistic plant-herbivore bipartite networks using common gardens in two contrasting environments comprised of aspen trees with differing evolutionary histories of defence against herbivores. These networks were tightly connected owing to a high level of specialization of arthropod herbivores that spend a large proportion of the life cycle on aspen. The gardens were separated by ten degrees of latitude with resultant differences in abiotic conditions. We evaluated network metrics and reported similar connectance between gardens but greater numbers of links per species in the northern common garden. Interaction matrices revealed clear nestedness, indicating subsetting of the bipartite interactions into specialist divisions, in both the environmental and evolutionary aspen groups, although nestedness values were only significant in the northern garden. Variation in plant vulnerability, measured as the frequency of herbivore specialization in the aspen population, was significantly partitioned by environment (common garden) but not by evolutionary origin of the aspens. Significant values of modularity were observed in all network matrices. Trait-matching indicated that growth traits, leaf morphology, and phenolic metabolites affected modular structure in both the garden and evolutionary groups, whereas extra-floral nectaries had little influence. Further examination of module configuration revealed that plant vulnerability explained considerable variance in web structure. The contrasting conditions between the two gardens resulted in bottom-up effects of the environment, which most strongly influenced the overall network architecture, however, the aspen groups with dissimilar evolutionary history also showed contrasting degrees of nestedness and modularity. Our research therefore shows that, while evolution does affect the structure of aspen-herbivore bipartite networks, the role of environmental variations is a dominant constraint.

Comparative physiology of allopatric Populus species: geographic clines in photosynthesis, height growth, and carbon isotope discrimination in common gardens.

Populus species with wide geographic ranges display strong adaptation to local environments. We studied the clinal patterns in phenology and ecophysiology in allopatric Populus species adapted to similar environments on different continents under common garden settings. As a result of climatic adaptation, both Populus tremula L. and Populus balsamifera L. display latitudinal clines in photosynthetic rates (A), whereby high-latitude trees of P. tremula had higher A compared to low-latitude trees and nearly so in P. balsamifera (p = 0.06). Stomatal conductance (g s) and chlorophyll content index (CCI) follow similar latitudinal trends. However, foliar nitrogen was positively correlated with latitude in P. balsamifera and negatively correlated in P. tremula. No significant trends in carbon isotope composition of the leaf tissue (δ(13)C) were observed for both species; but, intrinsic water-use efficiency (WUEi) was negatively correlated with the latitude of origin in P. balsamifera. In spite of intrinsically higher A, high-latitude trees in both common gardens accomplished less height gain as a result of early bud set. Thus, shoot biomass was determined by height elongation duration (HED), which was well approximated by the number of days available for free growth between bud flush and bud set. We highlight the shortcoming of unreplicated outdoor common gardens for tree improvement and the crucial role of photoperiod in limiting height growth, further complicating interpretation of other secondary effects.