Genome-wide association study and network analysis of in vitro transformation in Populus trichocarpa support key roles of diverse phytohormone pathways and cross talk.

Wide variation in amenability to transformation and regeneration (TR) among many plant species and genotypes presents a challenge to the use of genetic engineering in research and breeding. To help understand the causes of this variation, we performed association mapping and network analysis using a population of 1204 wild trees of Populus trichocarpa (black cottonwood). To enable precise and high-throughput phenotyping of callus and shoot TR, we developed a computer vision system that cross-referenced complementary red, green, and blue (RGB) and fluorescent-hyperspectral images. We performed association mapping using single-marker and combined variant methods, followed by statistical tests for epistasis and integration of published multi-omic datasets to identify likely regulatory hubs. We report 409 candidate genes implicated by associations within 5 kb of coding sequences, and epistasis tests implicated 81 of these candidate genes as regulators of one another. Gene ontology terms related to protein-protein interactions and transcriptional regulation are overrepresented, among others. In addition to auxin and cytokinin pathways long established as critical to TR, our results highlight the importance of stress and wounding pathways. Potential regulatory hubs of signaling within and across these pathways include GROWTH REGULATORY FACTOR 1 (GRF1), PHOSPHATIDYLINOSITOL 4-KINASE ß1 (PI-4Kß1), and OBF-BINDING PROTEIN 1 (OBP1).

Overexpression of REDUCED WALL ACETYLATION C increases xylan acetylation and biomass recalcitrance in Populus.

Plant lignocellulosic biomass, i.e. secondary cell walls of plants, is a vital alternative source for bioenergy. However, the acetylation of xylan in secondary cell walls impedes the conversion of biomass to biofuels. Previous studies have shown that REDUCED WALL ACETYLATION (RWA) proteins are directly involved in the acetylation of xylan but the regulatory mechanism of RWAs is not fully understood. In this study, we demonstrate that overexpression of a Populus trichocarpa PtRWA-C gene increases the level of xylan acetylation and increases the lignin content and S/G ratio, ultimately yielding poplar woody biomass with reduced saccharification efficiency. Furthermore, through gene coexpression network and expression quantitative trait loci (eQTL) analysis, we found that PtRWA-C was regulated not only by the secondary cell wall hierarchical regulatory network but also by an AP2 family transcription factor HARDY (HRD). Specifically, HRD activates PtRWA-C expression by directly binding to the PtRWA-C promoter, which is also the cis-eQTL for PtRWA-C. Taken together, our findings provide insights into the functional roles of PtRWA-C in xylan acetylation and consequently saccharification and shed light on synthetic biology approaches to manipulate this gene and alter cell wall properties. These findings have substantial implications for genetic engineering of woody species, which could be used as a sustainable source of biofuels, valuable biochemicals, and biomaterials.

Agave REVEILLE1 regulates the onset and release of seasonal dormancy in Populus.

Deciduous woody plants like poplar (Populus spp.) have seasonal bud dormancy. It has been challenging to simultaneously delay the onset of bud dormancy in the fall and advance bud break in the spring, as bud dormancy, and bud break were thought to be controlled by different genetic factors. Here, we demonstrate that heterologous expression of the REVEILLE1 gene (named AaRVE1) from Agave (Agave americana) not only delays the onset of bud dormancy but also accelerates bud break in poplar in field trials. AaRVE1 heterologous expression increases poplar biomass yield by 166% in the greenhouse. Furthermore, we reveal that heterologous expression of AaRVE1 increases cytokinin contents, represses multiple dormancy-related genes, and up-regulates bud break-related genes, and that AaRVE1 functions as a transcriptional repressor and regulates the activity of the DORMANCY-ASSOCIATED PROTEIN 1 (DRM1) promoter. Our findings demonstrate that AaRVE1 appears to function as a regulator of bud dormancy and bud break, which has important implications for extending the growing season of deciduous trees in frost-free temperate and subtropical regions to increase crop yield.

Expression quantitative trait loci mapping identified PtrXB38 as a key hub gene in adventitious root development in Populus.

Plant establishment requires the formation and development of an extensive root system with architecture modulated by complex genetic networks. Here, we report the identification of the PtrXB38 gene as an expression quantitative trait loci (eQTL) hotspot, mapped using 390 leaf and 444 xylem Populus trichocarpa transcriptomes. Among predicted targets of this trans-eQTL were genes involved in plant hormone responses and root development. Overexpression of PtrXB38 in Populus led to significant increases in callusing and formation of both stem-born roots and base-born adventitious roots. Omics studies revealed that genes and proteins controlling auxin transport and signaling were involved in PtrXB38-mediated adventitious root formation. Protein-protein interaction assays indicated that PtrXB38 interacts with components of endosomal sorting complexes required for transport machinery, implying that PtrXB38-regulated root development may be mediated by regulating endocytosis pathway. Taken together, this work identified a crucial root development regulator and sheds light on the discovery of other plant developmental regulators through combining eQTL mapping and omics approaches.

De Novo Assembly and Annotation of 11 Diverse Shrub Willow (Salix) Genomes Reveals Novel Gene Organization in Sex-Linked Regions.

Poplar and willow species in the Salicaceae are dioecious, yet have been shown to use different sex determination systems located on different chromosomes. Willows in the subgenus Vetrix are interesting for comparative studies of sex determination systems, yet genomic resources for these species are still quite limited. Only a few annotated reference genome assemblies are available, despite many species in use in breeding programs. Here we present de novo assemblies and annotations of 11 shrub willow genomes from six species. Copy number variation of candidate sex determination genes within each genome was characterized and revealed remarkable differences in putative master regulator gene duplication and deletion. We also analyzed copy number and expression of candidate genes involved in floral secondary metabolism, and identified substantial variation across genotypes, which can be used for parental selection in breeding programs. Lastly, we report on a genotype that produces only female descendants and identified gene presence/absence variation in the mitochondrial genome that may be responsible for this unusual inheritance.

Ecosystem consequences of introducing plant growth promoting rhizobacteria to managed systems and potential legacy effects.

The rapidly growing industry of crop biostimulants leverages the application of plant growth promoting rhizobacteria (PGPR) to promote plant growth and health. However, introducing nonnative rhizobacteria may impact other aspects of ecosystem functioning and have legacy effects; these potential consequences are largely unexplored. Nontarget consequences of PGPR may include changes in resident microbiomes, nutrient cycling, pollinator services, functioning of other herbivores, disease suppression, and organic matter persistence. Importantly, we lack knowledge of whether these ecosystem effects may manifest in adjacent ecosystems. The introduced PGPR can leave a functional legacy whether they persist in the community or not. Legacy effects include shifts in resident microbiomes and their temporal dynamics, horizontal transfer of genes from the PGPR to resident taxa, and changes in resident functional groups and interaction networks. Ecosystem functions may be affected by legacies PGPR leave following niche construction, such as when PGPR alter soil pH that in turn alters biogeochemical cycling rates. Here, we highlight new research directions to elucidate how introduced PGPR impact resident microbiomes and ecosystem functions and their capacity for legacy effects.

Towards engineering ectomycorrhization into switchgrass bioenergy crops via a lectin receptor-like kinase.

Soil-borne microbes can establish compatible relationships with host plants, providing a large variety of nutritive and protective compounds in exchange for photosynthesized sugars. However, the molecular mechanisms mediating the establishment of these beneficial relationships remain unclear. Our previous genetic mapping and whole-genome resequencing studies identified a gene deletion event of a Populus trichocarpa lectin receptor-like kinase gene PtLecRLK1 in Populus deltoides that was associated with poor-root colonization by the ectomycorrhizal fungus Laccaria bicolor. By introducing PtLecRLK1 into a perennial grass known to be a non-host of L. bicolor, switchgrass (Panicum virgatum L.), we found that L. bicolor colonizes ZmUbipro-PtLecRLK1 transgenic switchgrass roots, which illustrates that the introduction of PtLecRLK1 has the potential to convert a non-host to a host of L. bicolor. Furthermore, transcriptomic and proteomic analyses on inoculated-transgenic switchgrass roots revealed genes/proteins overrepresented in the compatible interaction and underrepresented in the pathogenic defence pathway, consistent with the view that pathogenic defence response is down-regulated during compatible interaction. Metabolomic profiling revealed that root colonization in the transgenic switchgrass was associated with an increase in N-containing metabolites and a decrease in organic acids, sugars, and aromatic hydroxycinnamate conjugates, which are often seen in the early steps of establishing compatible interactions. These studies illustrate that PtLecRLK1 is able to render a plant susceptible to colonization by the ectomycorrhizal fungus L. bicolor and shed light on engineering mycorrhizal symbiosis into a non-host to enhance plant productivity and fitness on marginal lands.

Deep neural network improves the estimation of polygenic risk scores for breast cancer.

Polygenic risk scores (PRS) estimate the genetic risk of an individual for a complex disease based on many genetic variants across the whole genome. In this study, we compared a series of computational models for estimation of breast cancer PRS. A deep neural network (DNN) was found to outperform alternative machine learning techniques and established statistical algorithms, including BLUP, BayesA, and LDpred. In the test cohort with 50% prevalence, the Area Under the receiver operating characteristic Curve (AUC) were 67.4% for DNN, 64.2% for BLUP, 64.5% for BayesA, and 62.4% for LDpred. BLUP, BayesA, and LPpred all generated PRS that followed a normal distribution in the case population. However, the PRS generated by DNN in the case population followed a bimodal distribution composed of two normal distributions with distinctly different means. This suggests that DNN was able to separate the case population into a high-genetic-risk case subpopulation with an average PRS significantly higher than the control population and a normal-genetic-risk case subpopulation with an average PRS similar to the control population. This allowed DNN to achieve 18.8% recall at 90% precision in the test cohort with 50% prevalence, which can be extrapolated to 65.4% recall at 20% precision in a general population with 12% prevalence. Interpretation of the DNN model identified salient variants that were assigned insignificant p values by association studies, but were important for DNN prediction. These variants may be associated with the phenotype through nonlinear relationships.

A 5-Enolpyruvylshikimate 3-Phosphate Synthase Functions as a Transcriptional Repressor in Populus.

Long-lived perennial plants, with distinctive habits of inter-annual growth, defense, and physiology, are of great economic and ecological importance. However, some biological mechanisms resulting from genome duplication and functional divergence of genes in these systems remain poorly studied. Here, we discovered an association between a poplar (Populus trichocarpa) 5-enolpyruvylshikimate 3-phosphate synthase gene (PtrEPSP) and lignin biosynthesis. Functional characterization of PtrEPSP revealed that this isoform possesses a helix-turn-helix motif in the N terminus and can function as a transcriptional repressor that regulates expression of genes in the phenylpropanoid pathway in addition to performing its canonical biosynthesis function in the shikimate pathway. We demonstrated that this isoform can localize in the nucleus and specifically binds to the promoter and represses the expression of a SLEEPER-like transcriptional regulator, which itself specifically binds to the promoter and represses the expression of PtrMYB021 (known as MYB46 in Arabidopsis thaliana), a master regulator of the phenylpropanoid pathway and lignin biosynthesis. Analyses of overexpression and RNAi lines targeting PtrEPSP confirmed the predicted changes in PtrMYB021 expression patterns. These results demonstrate that PtrEPSP in its regulatory form and PtrhAT form a transcriptional hierarchy regulating phenylpropanoid pathway and lignin biosynthesis in Populus.

Metabolomic Patterns of Septoria Canker Resistant and Susceptible Populus trichocarpa Genotypes 24 Hours Postinoculation.

Sphaerulina musiva is an economically and ecologically important fungal pathogen that causes Septoria stem canker and leaf spot disease of Populus species. To bridge the gap between genetic markers and structural barriers previously found to be linked to Septoria canker disease resistance in poplar, we used hydrophilic interaction liquid chromatography and tandem mass spectrometry to identify and quantify metabolites involved with signaling and cell wall remodeling. Fluctuations in signaling molecules, organic acids, amino acids, sterols, phenolics, and saccharides in resistant and susceptible P. trichocarpa inoculated with S. musiva were observed. The patterns of 222 metabolites in the resistant host implicate systemic acquired resistance (SAR), cell wall apposition, and lignin deposition as modes of resistance to this hemibiotrophic pathogen. This pattern is consistent with the expected response to the biotrophic phase of S. musiva colonization during the first 24 h postinoculation. The fungal pathogen metabolized key regulatory signals of SAR, other phenolics, and precursors of lignin biosynthesis that were depleted in the susceptible host. This is the first study to characterize metabolites associated with the response to initial colonization by S. musiva between resistant and susceptible hosts.

Association mapping, transcriptomics, and transient expression identify candidate genes mediating plant-pathogen interactions in a tree.

Invasive microbes causing diseases such as sudden oak death negatively affect ecosystems and economies around the world. The deployment of resistant genotypes for combating introduced diseases typically relies on breeding programs that can take decades to complete. To demonstrate how this process can be accelerated, we employed a genome-wide association mapping of ca 1,000 resequenced Populus trichocarpa trees individually challenged with Sphaerulina musiva, an invasive fungal pathogen. Among significant associations, three loci associated with resistance were identified and predicted to encode one putative membrane-bound L-type receptor-like kinase and two receptor-like proteins. A susceptibility-associated locus was predicted to encode a putative G-type D-mannose-binding receptor-like kinase. Multiple lines of evidence, including allele analysis, transcriptomics, binding assays, and overexpression, support the hypothesized function of these candidate genes in the P. trichocarpa response to S. musiva.

Role of Reactive Oxygen Species and Hormones in Plant Responses to Temperature Changes.

Temperature stress is one of the major abiotic stresses that adversely affect agricultural productivity worldwide. Temperatures beyond a plant's physiological optimum can trigger significant physiological and biochemical perturbations, reducing plant growth and tolerance to stress. Improving a plant's tolerance to these temperature fluctuations requires a deep understanding of its responses to environmental change. To adapt to temperature fluctuations, plants tailor their acclimatory signal transduction events, and specifically, cellular redox state, that are governed by plant hormones, reactive oxygen species (ROS) regulatory systems, and other molecular components. The role of ROS in plants as important signaling molecules during stress acclimation has recently been established. Here, hormone-triggered ROS produced by NADPH oxidases, feedback regulation, and integrated signaling events during temperature stress activate stress-response pathways and induce acclimation or defense mechanisms. At the other extreme, excess ROS accumulation, following temperature-induced oxidative stress, can have negative consequences on plant growth and stress acclimation. The excessive ROS is regulated by the ROS scavenging system, which subsequently promotes plant tolerance. All these signaling events, including crosstalk between hormones and ROS, modify the plant's transcriptomic, metabolomic, and biochemical states and promote plant acclimation, tolerance, and survival. Here, we provide a comprehensive review of the ROS, hormones, and their joint role in shaping a plant's responses to high and low temperatures, and we conclude by outlining hormone/ROS-regulated plant responsive strategies for developing stress-tolerant crops to combat temperature changes.

Overexpression of a Prefoldin ß subunit gene reduces biomass recalcitrance in the bioenergy crop Populus.

Prefoldin (PFD) is a group II chaperonin that is ubiquitously present in the eukaryotic kingdom. Six subunits (PFD1-6) form a jellyfish-like heterohexameric PFD complex and function in protein folding and cytoskeleton organization. However, little is known about its function in plant cell wall-related processes. Here, we report the functional characterization of a PFD gene from Populus deltoides, designated as PdPFD2.2. There are two copies of PFD2 in Populus, and PdPFD2.2 was ubiquitously expressed with high transcript abundance in the cambial region. PdPFD2.2 can physically interact with DELLA protein RGA1_8g, and its subcellular localization is affected by the interaction. In P. deltoides transgenic plants overexpressing PdPFD2.2, the lignin syringyl/guaiacyl ratio was increased, but cellulose content and crystallinity index were unchanged. In addition, the total released sugar (glucose and xylose) amounts were increased by 7.6% and 6.1%, respectively, in two transgenic lines. Transcriptomic and metabolomic analyses revealed that secondary metabolic pathways, including lignin and flavonoid biosynthesis, were affected by overexpressing PdPFD2.2. A total of eight hub transcription factors (TFs) were identified based on TF binding sites of differentially expressed genes in Populus transgenic plants overexpressing PdPFD2.2. In addition, several known cell wall-related TFs, such as MYB3, MYB4, MYB7, TT8 and XND1, were affected by overexpression of PdPFD2.2. These results suggest that overexpression of PdPFD2.2 can reduce biomass recalcitrance and PdPFD2.2 is a promising target for genetic engineering to improve feedstock characteristics to enhance biofuel conversion and reduce the cost of lignocellulosic biofuel production.

Arabidopsis C-terminal binding protein ANGUSTIFOLIA modulates transcriptional co-regulation of MYB46 and WRKY33.

The apparent antagonism between salicylic acid (SA) and jasmonic acid (JA)/ethylene (ET) signalling resulting in trade-offs between defence against (hemi)biotrophic and necrotrophic pathogens has been widely described across multiple plant species. However, the underlying mechanism remains to be fully established. The molecular and cellular functions of ANGUSTIFOLIA (AN) were characterised, and its role in regulating the pathogenic response was studied in Arabidopsis. We demonstrated that AN, a plant homologue of mammalian C-TERMINAL BINDING PROTEIN (CtBP), antagonistically regulates plant resistance to the hemibiotrophic pathogen Pseudomonas syringae and the necrotrophic pathogen Botrytis cinerea. Consistent with phenotypic observations, transcription of genes involved in SA and JA/ET pathways was antagonistically regulated by AN. By interacting with another nuclear protein TYROSYL-DNA PHOSPHODIESTERASE1 (TDP1), AN imposes transcriptional repression on MYB46, encoding a transcriptional activator of PHENYLALANINE AMMONIA-LYASE (PAL) genes which are required for SA biosynthesis, while releasing TDP1-imposed transcriptional repression on WRKY33, a master regulator of the JA/ET signalling pathway. These findings demonstrate that transcriptional co-regulation of MYB46 and WRKY33 by AN mediates the coordination of SA and JA/ET pathways to optimise defences against (hemi)biotrophic and necrotrophic pathogens.

The Physcomitrella patens chromosome-scale assembly reveals moss genome structure and evolution.

The draft genome of the moss model, Physcomitrella patens, comprised approximately 2000 unordered scaffolds. In order to enable analyses of genome structure and evolution we generated a chromosome-scale genome assembly using genetic linkage as well as (end) sequencing of long DNA fragments. We find that 57% of the genome comprises transposable elements (TEs), some of which may be actively transposing during the life cycle. Unlike in flowering plant genomes, gene- and TE-rich regions show an overall even distribution along the chromosomes. However, the chromosomes are mono-centric with peaks of a class of Copia elements potentially coinciding with centromeres. Gene body methylation is evident in 5.7% of the protein-coding genes, typically coinciding with low GC and low expression. Some giant virus insertions are transcriptionally active and might protect gametes from viral infection via siRNA mediated silencing. Structure-based detection methods show that the genome evolved via two rounds of whole genome duplications (WGDs), apparently common in mosses but not in liverworts and hornworts. Several hundred genes are present in colinear regions conserved since the last common ancestor of plants. These syntenic regions are enriched for functions related to plant-specific cell growth and tissue organization. The P. patens genome lacks the TE-rich pericentromeric and gene-rich distal regions typical for most flowering plant genomes. More non-seed plant genomes are needed to unravel how plant genomes evolve, and to understand whether the P. patens genome structure is typical for mosses or bryophytes.

Multitrait genome-wide association analysis of Populus trichocarpa identifies key polymorphisms controlling morphological and physiological traits.

Genome-wide association studies (GWAS) have great promise for identifying the loci that contribute to adaptive variation, but the complex genetic architecture of many quantitative traits presents a substantial challenge. We measured 14 morphological and physiological traits and identified single nucleotide polymorphism (SNP)-phenotype associations in a Populus trichocarpa population distributed from California, USA to British Columbia, Canada. We used whole-genome resequencing data of 882 trees with more than 6.78 million SNPs, coupled with multitrait association to detect polymorphisms with potentially pleiotropic effects. Candidate genes were validated with functional data. Broad-sense heritability (H2 ) ranged from 0.30 to 0.56 for morphological traits and 0.08 to 0.36 for physiological traits. In total, 4 and 20 gene models were detected using the single-trait and multitrait association methods, respectively. Several of these associations were corroborated by additional lines of evidence, including co-expression networks, metabolite analyses, and direct confirmation of gene function through RNAi. Multitrait association identified many more significant associations than single-trait association, potentially revealing pleiotropic effects of individual genes. This approach can be particularly useful for challenging physiological traits such as water-use efficiency or complex traits such as leaf morphology, for which we were able to identify credible candidate genes by combining multitrait association with gene co-expression and co-methylation data.

Genetic Analysis of Physcomitrella patens Identifies ABSCISIC ACID NON-RESPONSIVE, a Regulator of ABA Responses Unique to Basal Land Plants and Required for Desiccation Tolerance.

The anatomically simple plants that first colonized land must have acquired molecular and biochemical adaptations to drought stress. Abscisic acid (ABA) coordinates responses leading to desiccation tolerance in all land plants. We identified ABA nonresponsive mutants in the model bryophyte Physcomitrella patens and genotyped a segregating population to map and identify the ABA NON-RESPONSIVE (ANR) gene encoding a modular protein kinase comprising an N-terminal PAS domain, a central EDR domain, and a C-terminal MAPKKK-like domain. anr mutants fail to accumulate dehydration tolerance-associated gene products in response to drought, ABA, or osmotic stress and do not acquire ABA-dependent desiccation tolerance. The crystal structure of the PAS domain, determined to 1.7-Å resolution, shows a conserved PAS-fold that dimerizes through a weak dimerization interface. Targeted mutagenesis of a conserved tryptophan residue within the PAS domain generates plants with ABA nonresponsive growth and strongly attenuated ABA-responsive gene expression, whereas deleting this domain retains a fully ABA-responsive phenotype. ANR orthologs are found in early-diverging land plant lineages and aquatic algae but are absent from more recently diverged vascular plants. We propose that ANR genes represent an ancestral adaptation that enabled drought stress survival of the first terrestrial colonizers but were lost during land plant evolution.

Genome-wide association studies and expression-based quantitative trait loci analyses reveal roles of HCT2 in caffeoylquinic acid biosynthesis and its regulation by defense-responsive transcription factors in Populus.

3-O-caffeoylquinic acid, also known as chlorogenic acid (CGA), functions as an intermediate in lignin biosynthesis in the phenylpropanoid pathway. It is widely distributed among numerous plant species and acts as an antioxidant in both plants and animals. Using GC-MS, we discovered consistent and extreme variation in CGA content across a population of 739 4-yr-old Populus trichocarpa accessions. We performed genome-wide association studies (GWAS) from 917 P. trichocarpa accessions and expression-based quantitative trait loci (eQTL) analyses to identify key regulators. The GWAS and eQTL analyses resolved an overlapped interval encompassing a hydroxycinnamoyl-CoA:shikimate hydroxycinnamoyl transferase 2 (PtHCT2) that was significantly associated with CGA and partially characterized metabolite abundances. PtHCT2 leaf expression was significantly correlated with CGA abundance and it was regulated by cis-eQTLs containing W-box for WRKY binding. Among all nine PtHCT homologs, PtHCT2 is the only one that responds to infection by the fungal pathogen Sphaerulina musiva (a Populus pathogen). Validation using protoplast-based transient expression system suggests that PtHCT2 is regulated by the defense-responsive WRKY. These results are consistent with reports of CGA functioning as an antioxidant in response to biotic stress. This study provides insights into data-driven and omics-based inference of gene function in woody species.

The Sphagnome Project: enabling ecological and evolutionary insights through a genus-level sequencing project.

Considerable progress has been made in ecological and evolutionary genetics with studies demonstrating how genes underlying plant and microbial traits can influence adaptation and even 'extend' to influence community structure and ecosystem level processes. Progress in this area is limited to model systems with deep genetic and genomic resources that often have negligible ecological impact or interest. Thus, important linkages between genetic adaptations and their consequences at organismal and ecological scales are often lacking. Here we introduce the Sphagnome Project, which incorporates genomics into a long-running history of Sphagnum research that has documented unparalleled contributions to peatland ecology, carbon sequestration, biogeochemistry, microbiome research, niche construction, and ecosystem engineering. The Sphagnome Project encompasses a genus-level sequencing effort that represents a new type of model system driven not only by genetic tractability, but by ecologically relevant questions and hypotheses.

Quantitative trait locus mapping of Populus bark features and stem diameter.

BACKGROUND: Bark plays important roles in photosynthate transport and storage, along with physical and chemical protection. Bark texture varies extensively among species, from smooth to fissured to deeply furrowed, but its genetic control is unknown. This study sought to determine the main genomic regions associated with natural variation in bark features and stem diameter. Quantitative trait loci (QTL) were mapped using an interspecific pseudo-backcross pedigree (Populus trichocarpa x P. deltoides and P. deltoides) for bark texture, bark thickness and diameter collected across three years, two sites and three biological replicates per site. RESULTS: QTL specific to bark texture were highly reproducible in shared intervals across sites, years and replicates. Significant positive correlations and co-localization between trait QTL suggest pleiotropic regulators or closely linked genes. A list of candidate genes with related putative function, location close to QTL maxima and with the highest expression level in the phloem, xylem and cambium was identified. CONCLUSION: Candidate genes for bark texture included an ortholog of Arabidopsis ANAC104 (PopNAC128), which plays a role in lignified fiber cell and ray development, as well as Pinin and Fasciclin (PopFLA) genes with a role in cell adhesion, cell shape and migration. The results presented in this study provide a basis for future genomic characterization of genes found within the QTL for bark texture, bark thickness and diameter in order to better understand stem and bark development in Populus and other woody perennial plants. The QTL mapping approach identified a list of prime candidate genes for further validation using functional genomics or forward genetics approaches.