The western redcedar genome reveals low genetic diversity in a self-compatible conifer.

We assembled the 9.8-Gbp genome of western redcedar (WRC; Thuja plicata), an ecologically and economically important conifer species of the Cupressaceae. The genome assembly, derived from a uniquely inbred tree produced through five generations of self-fertilization (selfing), was determined to be 86% complete by BUSCO analysis, one of the most complete genome assemblies for a conifer. Population genomic analysis revealed WRC to be one of the most genetically depauperate wild plant species, with an effective population size of approximately 300 and no significant genetic differentiation across its geographic range. Nucleotide diversity, π, is low for a continuous tree species, with many loci showing zero diversity, and the ratio of π at zero- to fourfold degenerate sites is relatively high (approximately 0.33), suggestive of weak purifying selection. Using an array of genetic lines derived from up to five generations of selfing, we explored the relationship between genetic diversity and mating system. Although overall heterozygosity was found to decline faster than expected during selfing, heterozygosity persisted at many loci, and nearly 100 loci were found to deviate from expectations of genetic drift, suggestive of associative overdominance. Nonreference alleles at such loci often harbor deleterious mutations and are rare in natural populations, implying that balanced polymorphisms are maintained by linkage to dominant beneficial alleles. This may account for how WRC remains responsive to natural and artificial selection, despite low genetic diversity.

Single-nuclei transcriptome analysis of the shoot apex vascular system differentiation in Populus.

Differentiation of stem cells in the plant apex gives rise to aerial tissues and organs. Presently, we lack a lineage map of the shoot apex cells in woody perennials - a crucial gap considering their role in determining primary and secondary growth. Here, we used single-nuclei RNA-sequencing to determine cell type-specific transcriptomes of the Populus vegetative shoot apex. We identified highly heterogeneous cell populations clustered into seven broad groups represented by 18 transcriptionally distinct cell clusters. Next, we established the developmental trajectories of the epidermis, leaf mesophyll and vascular tissue. Motivated by the high similarities between Populus and Arabidopsis cell population in the vegetative apex, we applied a pipeline for interspecific single-cell gene expression data integration. We contrasted the developmental trajectories of primary phloem and xylem formation in both species, establishing the first comparison of vascular development between a model annual herbaceous and a woody perennial plant species. Our results offer a valuable resource for investigating the principles underlying cell division and differentiation conserved between herbaceous and perennial species while also allowing us to examine species-specific differences at single-cell resolution.

Xylem systems genetics analysis reveals a key regulator of lignin biosynthesis in Populus deltoides.

Despite the growing resources and tools for high-throughput characterization and analysis of genomic information, the discovery of the genetic elements that regulate complex traits remains a challenge. Systems genetics is an emerging field that aims to understand the flow of biological information that underlies complex traits from genotype to phenotype. In this study, we used a systems genetics approach to identify and evaluate regulators of the lignin biosynthesis pathway in Populus deltoides by combining genome, transcriptome, and phenotype data from a population of 268 unrelated individuals of P. deltoides The discovery of lignin regulators began with the quantitative genetic analysis of the xylem transcriptome and resulted in the detection of 6706 and 4628 significant local- and distant-eQTL associations, respectively. Among the locally regulated genes, we identified the R2R3-MYB transcription factor MYB125 (Potri.003G114100) as a putative trans-regulator of the majority of genes in the lignin biosynthesis pathway. The expression of MYB125 in a diverse population positively correlated with lignin content. Furthermore, overexpression of MYB125 in transgenic poplar resulted in increased lignin content, as well as altered expression of genes in the lignin biosynthesis pathway. Altogether, our findings indicate that MYB125 is involved in the control of a transcriptional coexpression network of lignin biosynthesis genes during secondary cell wall formation in P. deltoides.

Spatiotemporal cytokinin response imaging and ISOPENTENYLTRANSFERASE 3 function in Medicago nodule development.

Most legumes can establish a symbiotic association with soil rhizobia that trigger the development of root nodules. These nodules host the rhizobia and allow them to fix nitrogen efficiently. The perception of bacterial lipo-chitooligosaccharides (LCOs) in the epidermis initiates a signaling cascade that allows rhizobial intracellular infection in the root and de-differentiation and activation of cell division that gives rise to the nodule. Thus, nodule organogenesis and rhizobial infection need to be coupled in space and time for successful nodulation. The plant hormone cytokinin (CK) contributes to the coordination of this process, acting as an essential positive regulator of nodule organogenesis. However, the temporal regulation of tissue-specific CK signaling and biosynthesis in response to LCOs or Sinorhizobium meliloti inoculation in Medicago truncatula remains poorly understood. In this study, using a fluorescence-based CK sensor (pTCSn::nls:tGFP), we performed a high-resolution tissue-specific temporal characterization of the sequential activation of CK response during root infection and nodule development in M. truncatula after inoculation with S. meliloti. Loss-of-function mutants of the CK-biosynthetic gene ISOPENTENYLTRANSFERASE 3 (IPT3) showed impairment of nodulation, suggesting that IPT3 is required for nodule development in M. truncatula. Simultaneous live imaging of pIPT3::nls:tdTOMATO and the CK sensor showed that IPT3 induction in the pericycle at the base of nodule primordium contributes to CK biosynthesis, which in turn promotes expression of positive regulators of nodule organogenesis in M. truncatula.

An LCO-responsive homolog of NODULE INCEPTION positively regulates lateral root formation in Populus sp.

The transcription factor NODULE INCEPTION (NIN) has been studied extensively for its multiple roles in root nodule symbiosis within plants of the nitrogen-fixing clade (NFC) that associate with soil bacteria, such as rhizobia and Frankia. However, NIN homologs are present in plants outside the NFC, suggesting a role in other developmental processes. Here, we show that the biofuel crop Populus sp., which is not part of the NFC, contains eight copies of NIN with diversified protein sequence and expression patterns. Lipo-chitooligosaccharides (LCOs) are produced by rhizobia and a wide range of fungi, including mycorrhizal ones, and act as symbiotic signals that promote lateral root formation. RNAseq analysis of Populus sp. treated with purified LCO showed induction of the PtNIN2 subfamily. Moreover, the expression of PtNIN2b correlated with the formation of lateral roots and was suppressed by cytokinin treatment. Constitutive expression of PtNIN2b overcame the inhibition of lateral root development by cytokinin under high nitrate conditions. Lateral root induction in response to LCOs likely represents an ancestral function of NIN retained and repurposed in nodulating plants, as we demonstrate that the role of NIN in LCO-induced root branching is conserved in both Populus sp. and legumes. We further established a visual marker of LCO perception in Populus sp. roots, the putative sulfotransferase PtSS1 that can be used to study symbiotic interactions with the bacterial and fungal symbionts of Populus sp.

The uncharacterized gene EVE contributes to vessel element dimensions in Populus.

The radiation of angiosperms led to the emergence of the vast majority of today's plant species and all our major food crops. Their extraordinary diversification occurred in conjunction with the evolution of a more efficient vascular system for the transport of water, composed of vessel elements. The physical dimensions of these water-conducting specialized cells have played a critical role in angiosperm evolution; they determine resistance to water flow, influence photosynthesis rate, and contribute to plant stature. However, the genetic factors that determine their dimensions are unclear. Here we show that a previously uncharacterized gene, ENLARGED VESSEL ELEMENT (EVE), contributes to the dimensions of vessel elements in Populus, impacting hydraulic conductivity. Our data suggest that EVE is localized in the plasma membrane and is involved in potassium uptake of differentiating xylem cells during vessel development. In plants, EVE first emerged in streptophyte algae, but expanded dramatically among vessel-containing angiosperms. The phylogeny, structure and composition of EVE indicates that it may have been involved in an ancient horizontal gene-transfer event.

Temporal change in chromatin accessibility predicts regulators of nodulation in Medicago truncatula.

BACKGROUND: Symbiotic associations between bacteria and leguminous plants lead to the formation of root nodules that fix nitrogen needed for sustainable agricultural systems. Symbiosis triggers extensive genome and transcriptome remodeling in the plant, yet an integrated understanding of the extent of chromatin changes and transcriptional networks that functionally regulate gene expression associated with symbiosis remains poorly understood. In particular, analyses of early temporal events driving this symbiosis have only captured correlative relationships between regulators and targets at mRNA level. Here, we characterize changes in transcriptome and chromatin accessibility in the model legume Medicago truncatula, in response to rhizobial signals that trigger the formation of root nodules. RESULTS: We profiled the temporal chromatin accessibility (ATAC-seq) and transcriptome (RNA-seq) dynamics of M. truncatula roots treated with bacterial small molecules called lipo-chitooligosaccharides that trigger host symbiotic pathways of nodule development. Using a novel approach, dynamic regulatory module networks, we integrated ATAC-seq and RNA-seq time courses to predict cis-regulatory elements and transcription factors that most significantly contribute to transcriptomic changes associated with symbiosis. Regulators involved in auxin (IAA4-5, SHY2), ethylene (EIN3, ERF1), and abscisic acid (ABI5) hormone response, as well as histone and DNA methylation (IBM1), emerged among those most predictive of transcriptome dynamics. RNAi-based knockdown of EIN3 and ERF1 reduced nodule number in M. truncatula validating the role of these predicted regulators in symbiosis between legumes and rhizobia. CONCLUSIONS: Our transcriptomic and chromatin accessibility datasets provide a valuable resource to understand the gene regulatory programs controlling the early stages of the dynamic process of symbiosis. The regulators identified provide potential targets for future experimental validation, and the engineering of nodulation in species is unable to establish that symbiosis naturally.

Functional and comparative genomics reveals conserved noncoding sequences in the nitrogen-fixing clade.

Nitrogen is one of the most inaccessible plant nutrients, but certain species have overcome this limitation by establishing symbiotic interactions with nitrogen-fixing bacteria in the root nodule. This root-nodule symbiosis (RNS) is restricted to species within a single clade of angiosperms, suggesting a critical, but undetermined, evolutionary event at the base of this clade. To identify putative regulatory sequences implicated in the evolution of RNS, we evaluated the genomes of 25 species capable of nodulation and identified 3091 conserved noncoding sequences (CNS) in the nitrogen-fixing clade (NFC). We show that the chromatin accessibility of 452 CNS correlates significantly with the regulation of genes responding to lipochitooligosaccharides in Medicago truncatula. These included 38 CNS in proximity to 19 known genes involved in RNS. Five such regions are upstream of MtCRE1, Cytokinin Response Element 1, required to activate a suite of downstream transcription factors necessary for nodulation in M. truncatula. Genetic complementation of an Mtcre1 mutant showed a significant decrease of nodulation in the absence of the five CNS, when they are driving the expression of a functional copy of MtCRE1. CNS identified in the NFC may harbor elements required for the regulation of genes controlling RNS in M. truncatula.

Molecular changes in Mesembryanthemum crystallinum guard cells underlying the C3 to CAM transition.

ABSTARCT: KEY MESSAGE: The timing and transcriptomic changes during the C3 to CAM transition of common ice plant support the notion that guard cells themselves can shift from C3 to CAM. Crassulacean acid metabolism (CAM) is a specialized type of photosynthesis: stomata close during the day, enhancing water conservation, and open at night, allowing CO2 uptake. Mesembryanthemum crystallinum (common ice plant) is a facultative CAM species that can shift from C3 photosynthesis to CAM under salt or drought stresses. However, the molecular mechanisms underlying the stress induced transition from C3 to CAM remain unknown. Here we determined the transition time from C3 to CAM in M. crystallinum under salt stress. In parallel, single-cell-type transcriptomic profiling by 3'-mRNA sequencing was conducted in isolated stomatal guard cells to determine the molecular changes in this key cell type during the transition. In total, 495 transcripts showed differential expression between control and salt-treated samples during the transition, including 285 known guard cell genes, seven CAM-related genes, 18 transcription factors, and 185 other genes previously not found to be expressed in guard cells. PEPC1 and PPCK1, which encode key enzymes of CAM photosynthesis, were up-regulated in guard cells after seven days of salt treatment, indicating that guard cells themselves can shift from C3 to CAM. This study provides important information towards introducing CAM stomatal behavior into C3 crops to enhance water use efficiency.

Independent and Joint-GWAS for growth traits in Eucalyptus by assembling genome-wide data for 3373 individuals across four breeding populations.

Genome-wide association studies (GWAS) in plants typically suffer from limited statistical power. An alternative to the logistical and cost challenge of increasing sample sizes is to gain power by meta-analysis using information from independent studies. We carried out GWAS for growth traits with six single-marker models and regional heritability mapping (RHM) in four Eucalyptus breeding populations independently and by Joint-GWAS, using gene and segment-based models, with data for 3373 individuals genotyped with a communal EUChip60KSNP platform. While single-single nucleotide polymorphism (SNP) GWAS hardly detected significant associations at high-stringency in each population, gene-based Joint-GWAS revealed nine genes significantly associated with tree height. Associations detected using single-SNP GWAS, RHM and Joint-GWAS set-based models explained on average 3-20% of the phenotypic variance. Whole-genome regression, conversely, captured 64-89% of the pedigree-based heritability in all populations. Several associations independently detected for the same SNPs in different populations provided unprecedented GWAS validation results in forest trees. Rare and common associations were discovered in eight genes involved in cell wall biosynthesis and lignification. With the increasing adoption of genomic prediction of complex phenotypes using shared SNPs and much larger tree breeding populations, Joint-GWAS approaches should provide increasing power to pinpoint discrete associations potentially useful toward tree breeding and molecular applications.

Genomic prediction in contrast to a genome-wide association study in explaining heritable variation of complex growth traits in breeding populations of Eucalyptus.

BACKGROUND: The advent of high-throughput genotyping technologies coupled to genomic prediction methods established a new paradigm to integrate genomics and breeding. We carried out whole-genome prediction and contrasted it to a genome-wide association study (GWAS) for growth traits in breeding populations of Eucalyptus benthamii (n =505) and Eucalyptus pellita (n =732). Both species are of increasing commercial interest for the development of germplasm adapted to environmental stresses. RESULTS: Predictive ability reached 0.16 in E. benthamii and 0.44 in E. pellita for diameter growth. Predictive abilities using either Genomic BLUP or different Bayesian methods were similar, suggesting that growth adequately fits the infinitesimal model. Genomic prediction models using ~5000-10,000 SNPs provided predictive abilities equivalent to using all 13,787 and 19,506 SNPs genotyped in the E. benthamii and E. pellita populations, respectively. No difference was detected in predictive ability when different sets of SNPs were utilized, based on position (equidistantly genome-wide, inside genes, linkage disequilibrium pruned or on single chromosomes), as long as the total number of SNPs used was above ~5000. Predictive abilities obtained by removing relatedness between training and validation sets fell near zero for E. benthamii and were halved for E. pellita. These results corroborate the current view that relatedness is the main driver of genomic prediction, although some short-range historical linkage disequilibrium (LD) was likely captured for E. pellita. A GWAS identified only one significant association for volume growth in E. pellita, illustrating the fact that while genome-wide regression is able to account for large proportions of the heritability, very little or none of it is captured into significant associations using GWAS in breeding populations of the size evaluated in this study. CONCLUSIONS: This study provides further experimental data supporting positive prospects of using genome-wide data to capture large proportions of trait heritability and predict growth traits in trees with accuracies equal or better than those attainable by phenotypic selection. Additionally, our results document the superiority of the whole-genome regression approach in accounting for large proportions of the heritability of complex traits such as growth in contrast to the limited value of the local GWAS approach toward breeding applications in forest trees.

Genome-wide association study reveals putative regulators of bioenergy traits in Populus deltoides.

Genome-wide association studies (GWAS) have been used extensively to dissect the genetic regulation of complex traits in plants. These studies have focused largely on the analysis of common genetic variants despite the abundance of rare polymorphisms in several species, and their potential role in trait variation. Here, we conducted the first GWAS in Populus deltoides, a genetically diverse keystone forest species in North America and an important short rotation woody crop for the bioenergy industry. We searched for associations between eight growth and wood composition traits, and common and low-frequency single-nucleotide polymorphisms detected by targeted resequencing of 18 153 genes in a population of 391 unrelated individuals. To increase power to detect associations with low-frequency variants, multiple-marker association tests were used in combination with single-marker association tests. Significant associations were discovered for all phenotypes and are indicative that low-frequency polymorphisms contribute to phenotypic variance of several bioenergy traits. Our results suggest that both common and low-frequency variants need to be considered for a comprehensive understanding of the genetic regulation of complex traits, particularly in species that carry large numbers of rare polymorphisms. These polymorphisms may be critical for the development of specialized plant feedstocks for bioenergy.

Chilling-responsive DEMETER-LIKE DNA demethylase mediates in poplar bud break.

Annual dormancy-growth cycle is a developmental and physiological process essential for the survival of deciduous trees in temperate and boreal forests. Seasonal control of shoot growth in woody perennials requires specific genetic programmes responding to environmental signals. The environmental-controlled mechanisms that regulate the shift between winter dormancy and the growth-promoting genetic programmes are still unknown. Here, we show that dynamics in genomic DNA methylation levels are involved in the regulation of dormancy-growth cycle in poplar. The reactivation of growth in the apical shoot during bud break process in spring is preceded by a progressive reduction of genomic DNA methylation in apex tissue. The induction in apex tissue of a chilling-dependent poplar DEMETER-LIKE 10 (PtaDML10) DNA demethylase precedes shoot growth reactivation. Transgenic poplars showing downregulation of PtaDML8/10 caused delayed bud break. Genome-wide transcriptome and methylome analysis and data mining revealed that the gene targets of DEMETER-LIKE-dependent DNA demethylation are genetically associated with bud break. These data point to a chilling-dependent DEMETER-like DNA demethylase mechanisms being involved in the shift from winter dormancy to a condition that precedes shoot apical vegetative growth in poplar.

Overexpression of DEMETER, a DNA demethylase, promotes early apical bud maturation in poplar.

The transition from active growth to dormancy is critical for the survival of perennial plants. We identified a DEMETER-like (CsDML) cDNA from a winter-enriched cDNA subtractive library in chestnut (Castanea sativa Mill.), an economically and ecologically important species. Next, we characterized this DNA demethylase and its putative ortholog in the more experimentally tractable hybrid poplar (Populus tremula × alba), under the signals that trigger bud dormancy in trees. We performed phylogenetic and protein sequence analysis, gene expression profiling, and 5-methyl-cytosine methylation immunodetection studies to evaluate the role of CsDML and its homolog in poplar, PtaDML6. Transgenic hybrid poplars overexpressing CsDML were produced and analysed. Short days and cold temperatures induced CsDML and PtaDML6. Overexpression of CsDML accelerated short-day-induced bud formation, specifically from Stages 1 to 0. Buds acquired a red-brown coloration earlier than wild-type plants, alongside with the up-regulation of flavonoid biosynthesis enzymes and accumulation of flavonoids in the shoot apical meristem and bud scales. Our data show that the CsDML gene induces bud formation needed for the survival of the apical meristem under the harsh conditions of winter.

Integration of genetic, genomic and transcriptomic information identifies putative regulators of adventitious root formation in Populus.

BACKGROUND: Adventitious roots (AR) develop from tissues other than the primary root, in a process physiologically regulated by phytohormones. Adventitious roots provide structural support and contribute to water and nutrient absorption, and are critical for commercial vegetative propagation of several crops. Here we quantified the number of AR, root architectural traits and root biomass in cuttings from a pseudo-backcross population of Populus deltoides and Populus trichocarpa. Quantitative trait loci (QTL) mapping and whole-transcriptome analysis of individuals with alternative QTL alleles for AR number were used to identify putative regulators of AR development. RESULTS: Parental individuals and progeny showed extensive segregation for AR developmental traits. Quantitative trait loci for number of AR mapped consistently in the same interval of linkage group (LG) II and LG XIV, explaining 7-10 % of the phenotypic variation. A time series transcriptome analysis identified 26,121 genes differentially expressed during AR development, particularly during the first 24 h after cuttings were harvested. Of those, 1929 genes were differentially regulated between individuals carrying alternative alleles for the two QTL for number of AR, in one or more time point. Eighty-one of these genes were physically located within the QTL intervals for number of AR, including putative homologs of the Arabidopsis genes SUPERROOT2 (SUR2) and TRYPTOPHAN SYNTHASE ALPHA CHAIN (TSA1), both of which are involved in the auxin indole-3-acetic acid (IAA) biosynthesis pathway. CONCLUSIONS: This study suggests the involvement of two genes of the tryptophan-dependent auxin biosynthesis pathway, SUR2 and TSA1, in the regulation of a critical trait for the clonal propagation of woody species. A possible model for this regulation is that poplar individuals that have poor AR formation synthesize auxin indole-3-acetic acid (IAA) primarily through the tryptophan (Trp) pathway. Much of the Trp pathway flux appears to be directed to the synthesis of indole glucosinolates (IG), as suggested by the over-expression of SUR2. Individuals that are efficient in AR formation may utilize alternative (non-Trp) pathways to synthesize IAA, based on the observation that they down-regulate the expression of TSA1, one of the critical steps in the synthesis of tryptophan.

Metabolic regulation of triacylglycerol accumulation in the green algae: identification of potential targets for engineering to improve oil yield.

The great need for more sustainable alternatives to fossil fuels has increased our research interests in algal biofuels. Microalgal cells, characterized by high photosynthetic efficiency and rapid cell division, are an excellent source of neutral lipids as potential fuel stocks. Various stress factors, especially nutrient-starvation conditions, induce an increased formation of lipid bodies filled with triacylglycerol in these cells. Here we review our knowledge base on glycerolipid synthesis in the green algae with an emphasis on recent studies on carbon flux, redistribution of lipids under nutrient-limiting conditions and its regulation. We discuss the contributions and limitations of classical and novel approaches used to elucidate the algal triacylglycerol biosynthetic pathway and its regulatory network in green algae. Also discussed are gaps in knowledge and suggestions for much needed research both on the biology of triacylglycerol accumulation and possible avenues to engineer improved algal strains.

Spatial genetic structure, genetic diversity and pollen dispersal in a harvested population of Astrocaryum aculeatum in the Brazilian Amazon.

BACKGROUND: Astrocaryum aculeatum is a palm tree species native to the tropical regions of South America, exploited commercially by local farmers for the pulp extracted from its fruits. The objective of this research was to compare the genetic diversity between adult plants and seedlings from open-pollinated seeds, quantify the pollen flow and dispersal, the spatial genetic structure, and the effective size of a population that has been continuously harvested for its fruits. The study was carried out in a natural population of A. aculeatum distributed over approximately 8 ha in the State of Amazonas (Brazil), separated by 400 m from the closest neighboring population. In total, 112 potential pollen donors, 12 mother plants and 120 offspring were mapped and genotyped. RESULTS: Genetic diversity was high for parents and the offspring. The fixation indexes for adults (F = -0.035) and offspring (F = -0.060) were negative and not significant. A significant spatial genetic structure was detected for the adult plants (up to the distance of 45 m) indicating short-distance seed dispersal. Paternity analysis detected 9.2 % of pollen immigration and the average distance of pollination within the population was 81 m. The average effective pollination neighborhood area between plants was 1.51 ha. CONCLUSIONS: Our results indicate that substantial introduction of new alleles has occurred in the population through pollen immigration, contributing to the maintenance of genetic diversity. Conservation efforts aimed at maintaining the gene pool of the current population or establishing new populations should utilize offspring from mother plants selected to be spaced by at least 50 m to prevent collecting seeds from relatives.

Aluminum tolerance in maize is associated with higher MATE1 gene copy number.

Genome structure variation, including copy number variation and presence/absence variation, comprises a large extent of maize genetic diversity; however, its effect on phenotypes remains largely unexplored. Here, we describe how copy number variation underlies a rare allele that contributes to maize aluminum (Al) tolerance. Al toxicity is the primary limitation for crop production on acid soils, which make up 50% of the world's potentially arable lands. In a recombinant inbred line mapping population, copy number variation of the Al tolerance gene multidrug and toxic compound extrusion 1 (MATE1) is the basis for the quantitative trait locus of largest effect on phenotypic variation. This expansion in MATE1 copy number is associated with higher MATE1 expression, which in turn results in superior Al tolerance. The three MATE1 copies are identical and are part of a tandem triplication. Only three maize inbred lines carrying the three-copy allele were identified from maize and teosinte diversity panels, indicating that copy number variation for MATE1 is a rare, and quite likely recent, event. These maize lines with higher MATE1 copy number are also Al-tolerant, have high MATE1 expression, and originate from regions of highly acidic soils. Our findings show a role for copy number variation in the adaptation of maize to acidic soils in the tropics and suggest that genome structural changes may be a rapid evolutionary response to new environments.

Genetical genomics of Populus leaf shape variation.

BACKGROUND: Leaf morphology varies extensively among plant species and is under strong genetic control. Mutagenic screens in model systems have identified genes and established molecular mechanisms regulating leaf initiation, development, and shape. However, it is not known whether this diversity across plant species is related to naturally occurring variation at these genes. Quantitative trait locus (QTL) analysis has revealed a polygenic control for leaf shape variation in different species suggesting that loci discovered by mutagenesis may only explain part of the naturally occurring variation in leaf shape. Here we undertook a genetical genomics study in a poplar intersectional pseudo-backcross pedigree to identify genetic factors controlling leaf shape. The approach combined QTL discovery in a genetic linkage map anchored to the Populus trichocarpa reference genome sequence and transcriptome analysis. RESULTS: A major QTL for leaf lamina width and length:width ratio was identified in multiple experiments that confirmed its stability. A transcriptome analysis of expanding leaf tissue contrasted gene expression between individuals with alternative QTL alleles, and identified an ADP-ribosylation factor (ARF) GTPase (PtARF1) as a candidate gene for regulating leaf morphology in this pedigree. ARF GTPases are critical elements in the vesicular trafficking machinery. Disruption of the vesicular trafficking function of ARF by the pharmacological agent Brefeldin A (BFA) altered leaf lateral growth in the narrow-leaf P. trichocarpa suggesting a molecular mechanism of leaf shape determination. Inhibition of the vesicular trafficking processes by BFA interferes with cycling of PIN proteins and causes their accumulation in intercellular compartments abolishing polar localization and disrupting normal auxin flux with potential effects on leaf expansion. CONCLUSIONS: In other model systems, ARF proteins have been shown to control the localization of auxin efflux carriers, which function to establish auxin gradients and apical-basal cell polarity in developing plant organs. Our results support a model where PtARF1 transcript abundance changes the dynamics of endocytosis-mediated PIN localization in leaf cells, thus affecting lateral auxin flux and subsequently lamina leaf expansion. This suggests that evolution of differential cellular polarity plays a significant role in leaf morphological variation observed in subgenera of genus Populus.

Discovering candidate genes that regulate resin canal number in Pinus taeda stems by integrating genetic analysis across environments, ages, and populations.

Genetically improving constitutive resin canal development in Pinus stems may enhance the capacity to synthesize terpenes for bark beetle resistance, chemical feedstocks, and biofuels. To discover genes that potentially regulate axial resin canal number (RCN), single nucleotide polymorphisms (SNPs) in 4027 genes were tested for association with RCN in two growth rings and three environments in a complex pedigree of 520 Pinus taeda individuals (CCLONES). The map locations of associated genes were compared with RCN quantitative trait loci (QTLs) in a (P. taeda × Pinus elliottii) × P. elliottii pseudo-backcross of 345 full-sibs (BC1). Resin canal number was heritable (h(2) ˜ 0.12-0.21) and positively genetically correlated with xylem growth (rg ˜ 0.32-0.72) and oleoresin flow (rg ˜ 0.15-0.51). Sixteen well-supported candidate regulators of RCN were discovered in CCLONES, including genes associated across sites and ages, unidirectionally associated with oleoresin flow and xylem growth, and mapped to RCN QTLs in BC1. Breeding is predicted to increase RCN 11% in one generation and could be accelerated with genomic selection at accuracies of 0.45-0.52 across environments. There is significant genetic variation for RCN in loblolly pine, which can be exploited in breeding for elevated terpene content.