Explore the vast range of titles available.

While several studies have investigated general properties of the genetic architecture of natural variation in gene expression, few of these have considered natural, outbreeding populations. In parallel, systems biology has established that a general feature of biological networks is that they are scale-free, rendering them buffered against random mutations. To date, few studies have attempted to examine the relationship between the selective processes acting to maintain natural variation of gene expression and the associated co-expression network structure. Here we utilised RNA-Sequencing to assay gene expression in winter buds undergoing bud flush in a natural population of Populus tremula, an outbreeding forest tree species. We performed expression Quantitative Trait Locus (eQTL) mapping and identified 164,290 significant eQTLs associating 6,241 unique genes (eGenes) with 147,419 unique SNPs (eSNPs). We found approximately four times as many local as distant eQTLs, with local eQTLs having significantly higher effect sizes. eQTLs were primarily located in regulatory regions of genes (UTRs or flanking regions), regardless of whether they were local or distant. We used the gene expression data to infer a co-expression network and investigated the relationship between network topology, the genetic architecture of gene expression and signatures of selection. Within the co-expression network, eGenes were underrepresented in network module cores (hubs) and overrepresented in the periphery of the network, with a negative correlation between eQTL effect size and network connectivity. We additionally found that module core genes have experienced stronger selective constraint on coding and non-coding sequence, with connectivity associated with signatures of selection. Our integrated genetics and genomics results suggest that purifying selection is the primary mechanism underlying the genetic architecture of natural variation in gene expression assayed in flushing leaf buds of P. tremula and that connectivity within the co-expression network is linked to the strength of purifying selection.

CONTENTS: Summary 909 I. Introduction 910 II. Genotyping 910 III. Phenotyping 911 IV. Study designs 912 V. The genetics of the ‘omics' 912 VI. Missing heritability: the dark matter of the genome 913 VII. Gene interactions 914 VIII. Many rare alleles 914 IX. Looking in the wrong place 914 X. Looking but not seeing 915 XI. Needles in a haystack 915 XII. Confounding effects 916 XIII. Replicating and verifying associations 916 XIV. The genetic architecture of quantitative traits in plants 917 XV. Outlook 918 Acknowledgements 919 References 919 SUMMARY: Association mapping is rapidly becoming the main method for dissecting the genetic architecture of complex traits in plants. Currently most association mapping studies in plants are preformed using sets of genes selected to be putative candidates for the trait of interest, but rapid developments in genomics will allow for genome-wide mapping in virtually any plant species in the near future. As the costs for genotyping are decreasing, the focus has shifted towards phenotyping. In plants, clonal replication and/or inbred lines allows for replicated phenotyping under many different environmental conditions. Reduced sequencing costs will increase the number of studies that use RNA sequencing data to perform expression quantitative trait locus (eQTL) mapping, which will increase our knowledge of how gene expression variation contributes to phenotypic variation. Current population sizes used in association mapping studies are modest in size and need to be greatly increased if mutations explaining less than a few per cent of the phenotypic variation are to be detected. Association mapping has started to yield insights into the genetic architecture of complex traits in plants, and future studies with greater genome coverage will help to elucidate how plants have managed to adapt to a wide variety of environmental conditions.

Accessing and exploring large-scale genomics data sets remains a significant challenge to researchers without specialist bioinformatics training. We present the integrated PlantGenIE.org platform for exploration of Populus, conifer and Arabidopsis genomics data, which includes expression networks and associated visualization tools. Standard features of a model organism database are provided, including genome browsers, gene list annotation, BLAST homology searches and gene information pages. Community annotation updating is supported via integration of WebApollo. We have produced an RNA-sequencing (RNA-Seq) expression atlas for Populus tremula and have integrated these data within the expression tools. An updated version of the COMPLEX resource for performing comparative plant expression analyses of gene coexpression network conservation between species has also been integrated. The PlantGenIE.org platform provides intuitive access to large-scale and genome-wide genomics data from model forest tree species, facilitating both community contributions to annotation improvement and tools supporting use of the included data resources to inform biological insight.

Abstract Hybridization and resulting introgression are important processes shaping the tree of life and appear to be far more common than previously thought. However, how the genome evolution was shaped by various genetic and evolutionary forces after hybridization remains unresolved. Here we used whole-genome resequencing data of 227 individuals from multiple widespread Populus species to characterize their contemporary patterns of hybridization and to quantify genomic signatures of past introgression. We observe a high frequency of contemporary hybridization and confirm that multiple previously ambiguous species are in fact F1 hybrids. Seven species were identified, which experienced different demographic histories that resulted in strikingly varied efficacy of selection and burdens of deleterious mutations. Frequent past introgression has been found to be a pervasive feature throughout the speciation of these Populus species. The retained introgressed regions, more generally, tend to contain reduced genetic load and to be located in regions of high recombination. We also find that in pairs of species with substantial differences in effective population size, introgressed regions are inferred to have undergone selective sweeps at greater than expected frequencies in the species with lower effective population size, suggesting that introgression likely have higher potential to provide beneficial variation for species with small populations. Our results, therefore, illustrate that demography and recombination have interplayed with both positive and negative selection in determining the genomic evolution after hybridization.

Despite the global economic and ecological importance of forest trees, the genomic basis of differential adaptation and speciation in tree species is still poorly understood. Populus tremula and Populus tremuloides are two of the most widespread tree species in the Northern Hemisphere. Using whole-genome re-sequencing data of 24 P. tremula and 22 P. tremuloides individuals, we find that the two species diverged ∼2.2–3.1 million years ago, coinciding with the severing of the Bering land bridge and the onset of dramatic climatic oscillations during the Pleistocene. Both species have experienced substantial population expansions following long-term declines after species divergence. We detect widespread and heterogeneous genomic differentiation between species, and in accordance with the expectation of allopatric speciation, coalescent simulations suggest that neutral evolutionary processes can account for most of the observed patterns of genetic differentiation. However, there is an excess of regions exhibiting extreme differentiation relative to those expected under demographic simulations, which is indicative of the action of natural selection. Overall genetic differentiation is negatively associated with recombination rate in both species, providing strong support for a role of linked selection in generating the heterogeneous genomic landscape of differentiation between species. Finally, we identify a number of candidate regions and genes that may have been subject to positive and/or balancing selection during the speciation process.

Seasonal cues influence several aspects of the secondary growth of tree stems, including cambial activity, wood chemistry, and transition to latewood formation. We investigated seasonal changes in cambial activity, secondary cell wall formation, and tracheid cell death in woody tissues of Norway spruce (Picea abies) throughout one seasonal cycle. RNA sequencing was performed simultaneously in both the xylem and cambium/phloem tissues of the stem. Principal component analysis revealed gradual shifts in the transcriptomes that followed a chronological order throughout the season. A notable remodeling of the transcriptome was observed in the winter, with many genes having maximal expression during the coldest months of the year. A highly coexpressed set of monolignol biosynthesis genes showed high expression during the period of secondary cell wall formation as well as a second peak in midwinter. This midwinter peak in expression did not trigger lignin deposition, as determined by pyrolysis-gas chromatography/mass spectrometry. Coexpression consensus network analyses suggested the involvement of transcription factors belonging to the ASYMMETRIC LEAVES2/LATERAL ORGAN BOUNDARIES and MYELOBLASTOSIS-HOMEOBOX families in the seasonal control of secondary cell wall formation of tracheids. Interestingly, the lifetime of the latewood tracheids stretched beyond the winter dormancy period, correlating with a lack of cell death-related gene expression. Our transcriptomic analyses combined with phylogenetic and microscopic analyses also identified the cellulose and lignin biosynthetic genes and putative regulators for latewood formation and tracheid cell death in Norway spruce, providing a toolbox for further physiological and functional assays of these important phase transitions.

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 PtFT 2 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.

Abstract 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.

A comparative transcriptomic study and a single-cell metabolome analysis were combined to determine whether parenchymal ray cells contribute to the biosynthesis of monolignols in the lignifying xylem of Norway spruce (Picea abies). Ray parenchymal cells may function in the lignification of upright tracheids by supplying monolignols. To test this hypothesis, parenchymal ray cells and upright tracheids were dissected with laser-capture microdissection from tangential cryosections of developing xylem of spruce trees. The transcriptome analysis revealed that among the genes involved in processes typical for vascular tissues, genes encoding cell wall biogenesis-related enzymes were highly expressed in both developing tracheids and ray cells. Interestingly, most of the shikimate and monolignol biosynthesis pathway-related genes were equally expressed in both cell types. Nonetheless, 1,073 differentially expressed genes were detected between developing ray cells and tracheids, among which a set of genes expressed only in ray cells was identified. In situ single cell metabolomics of semi-intact plants by picoliter pressure probe-electrospray ionization-mass spectrometry detected monolignols and their glycoconjugates in both cell types, indicating that the biosynthetic route for monolignols is active in both upright tracheids and parenchymal ray cells. The data strongly support the hypothesis that in developing xylem, ray cells produce monolignols that contribute to lignification of tracheid cell walls.

DNA methylation plays important roles in many biological processes, such as silencing of transposable elements, imprinting, and regulating gene expression. Many studies of DNA methylation have shown its essential roles in angiosperms (flowering plants). However, few studies have examined the roles and patterns of DNA methylation in gymnosperms. Here, we present genome-wide high coverage single-base resolution methylation maps of Norway spruce (Picea abies) from both needles and somatic embryogenesis culture cells via whole genome bisulfite sequencing. On average, DNA methylation levels of CG and CHG of Norway spruce were higher than most other plants studied. CHH methylation was found at a relatively low level; however, at least one copy of most of the RNA-directed DNA methylation pathway genes was found in Norway spruce, and CHH methylation was correlated with levels of siRNAs. In comparison with needles, somatic embryogenesis culture cells that are used for clonally propagating spruce trees showed lower levels of CG and CHG methylation but higher level of CHH methylation, suggesting that like in other species, these culture cells show abnormal methylation patterns.