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

GWAS supported by computer vision identifies large numbers of candidate regulators of in planta regeneration in Populus trichocarpa.

Plant regeneration is an important dimension of plant propagation and a key step in the production of transgenic plants. However, regeneration capacity varies widely among genotypes and species, the molecular basis of which is largely unknown. Association mapping methods such as genome-wide association studies (GWAS) have long demonstrated abilities to help uncover the genetic basis of trait variation in plants; however, the performance of these methods depends on the accuracy and scale of phenotyping. To enable a large-scale GWAS of in planta callus and shoot regeneration in the model tree Populus, we developed a phenomics workflow involving semantic segmentation to quantify regenerating plant tissues over time. We found that the resulting statistics were of highly non-normal distributions, and thus employed transformations or permutations to avoid violating assumptions of linear models used in GWAS. We report over 200 statistically supported quantitative trait loci (QTLs), with genes encompassing or near to top QTLs including regulators of cell adhesion, stress signaling, and hormone signaling pathways, as well as other diverse functions. Our results encourage models of hormonal signaling during plant regeneration to consider keystone roles of stress-related signaling (e.g. involving jasmonates and salicylic acid), in addition to the auxin and cytokinin pathways commonly considered. The putative regulatory genes and biological processes we identified provide new insights into the biological complexity of plant regeneration, and may serve as new reagents for improving regeneration and transformation of recalcitrant genotypes and species.

GWAS identifies candidate genes controlling adventitious rooting in Populus trichocarpa.

Adventitious rooting (AR) is critical to the propagation, breeding, and genetic engineering of trees. The capacity for plants to undergo this process is highly heritable and of a polygenic nature; however, the basis of its genetic variation is largely uncharacterized. To identify genetic regulators of AR, we performed a genome-wide association study (GWAS) using 1148 genotypes of Populus trichocarpa. GWASs are often limited by the abilities of researchers to collect precise phenotype data on a high-throughput scale; to help overcome this limitation, we developed a computer vision system to measure an array of traits related to adventitious root development in poplar, including temporal measures of lateral and basal root length and area. GWAS was performed using multiple methods and significance thresholds to handle non-normal phenotype statistics and to gain statistical power. These analyses yielded a total of 277 unique associations, suggesting that genes that control rooting include regulators of hormone signaling, cell division and structure, reactive oxygen species signaling, and other processes with known roles in root development. Numerous genes with uncharacterized functions and/or cryptic roles were also identified. These candidates provide targets for functional analysis, including physiological and epistatic analyses, to better characterize the complex polygenic regulation of AR.

Variation in floral form of CRISPR knock-outs of the poplar homologs of LEAFY and AGAMOUS after FT heat-induced early flowering.

Plant migration and gene flow from genetically modified or exotic trees to nearby lands or by crossing with wild relatives is a major public and regulatory concern. Many genetic strategies exist to mitigate potential gene flow; however, the long delay in onset of flowering is a severe constraint to research progress. We used heat-induced FT overexpression to speed assessment of the expected floral phenotypes after CRISPR knockout of poplar homologs of the key floral genes, LEAFY and AGAMOUS. We selected events with previously characterized CRISPR-Cas9 induced biallelic changes then re-transformed them with the Arabidopsis thaliana FLOWERING LOCUS T (AtFT) gene under control of either a strong constitutive promoter or a heat-inducible promoter. We successfully obtained flowering in both a male and female clones of poplar, observing a wide range of inflorescence and floral forms among flowers, ramets, and insertion events. Overall, flowers obtained from the selected LFY and AG targeted events were consistent with what would be predicted for loss-of-function of these genes. LFY-targeted events showed small catkins with leaf-like organs, AG-targeted events had nested floral organs consistent with reduction in floral determinacy and absence of well-formed carpels or anthers. These findings demonstrate the great developmental plasticity of Populus flowers during genetically accelerated flowering, which may be of horticultural value. They also provide an informative early view of floral phenotypes and apparent sterility from knockouts of both these gene targets.

Knockout of floral and meiosis genes using CRISPR/Cas9 produces male-sterility in Eucalyptus without impacts on vegetative growth.

Eucalyptus spp. are widely cultivated for the production of pulp, energy, essential oils, and as ornamentals. However, their dispersal from plantings, especially when grown as an exotic, can cause ecological disruptions. To provide new tools for prevention of sexual dispersal by pollen as well as to induce male-sterility for hybrid breeding, we studied the clustered regularly interspaced short palindromic repeats (CRISPR)/Cas9-mediated knockout of three floral genes in both FT-expressing (early-flowering) and non-FT genotypes. We report male-sterile phenotypes resulting from knockout of the homologs of all three genes, including one involved in meiosis and two regulating early stages of pollen development. The targeted genes were Eucalyptus homologs of REC8 (EREC8), TAPETAL DEVELOPMENT AND FUNCTION 1 (ETDF1), and HECATE3 (EHEC3-like). The erec8 knockouts yielded abnormal pollen grains and a predominance of inviable pollen, whereas the etdf1 and ehec3-like knockouts produced virtually no pollen. In addition to male-sterility, both erec8 and ehec3-like knockouts may provide complete sterility because the failure of erec8 to undergo meiosis is expected to be independent of sex, and ehec3-like knockouts produce flowers with shortened styles and no visible stigmas. When comparing knockouts to controls in wild-type (non-early-flowering) backgrounds, we did not find visible morphological or statistical differences in vegetative traits, including average single-leaf mass, stem volume, density of oil glands, or chlorophyll in leaves. Loss-of-function mutations in any of these three genes show promise as a means of inducing male- or complete sterility without impacting vegetative development.

Robust High-Throughput Phenotyping with Deep Segmentation Enabled by a Web-Based Annotator.

The abilities of plant biologists and breeders to characterize the genetic basis of physiological traits are limited by their abilities to obtain quantitative data representing precise details of trait variation, and particularly to collect this data at a high-throughput scale with low cost. Although deep learning methods have demonstrated unprecedented potential to automate plant phenotyping, these methods commonly rely on large training sets that can be time-consuming to generate. Intelligent algorithms have therefore been proposed to enhance the productivity of these annotations and reduce human efforts. We propose a high-throughput phenotyping system which features a Graphical User Interface (GUI) and a novel interactive segmentation algorithm: Semantic-Guided Interactive Object Segmentation (SGIOS). By providing a user-friendly interface and intelligent assistance with annotation, this system offers potential to streamline and accelerate the generation of training sets, reducing the effort required by the user. Our evaluation shows that our proposed SGIOS model requires fewer user inputs compared to the state-of-art models for interactive segmentation. As a case study of the use of the GUI applied for genetic discovery in plants, we present an example of results from a preliminary genome-wide association study (GWAS) of in planta regeneration in Populus trichocarpa (poplar). We further demonstrate that the inclusion of a semantic prior map with SGIOS can accelerate the training process for future GWAS, using a sample of a dataset extracted from a poplar GWAS of in vitro regeneration. The capabilities of our phenotyping system surpass those of unassisted humans to rapidly and precisely phenotype our traits of interest. The scalability of this system enables large-scale phenomic screens that would otherwise be time-prohibitive, thereby providing increased power for GWAS, mutant screens, and other studies relying on large sample sizes to characterize the genetic basis of trait variation. Our user-friendly system can be used by researchers lacking a computational background, thus helping to democratize the use of deep segmentation as a tool for plant phenotyping.

Functional Diversification of Populus FLOWERING LOCUS D-LIKE3 Transcription Factor and Two Paralogs in Shoot Ontogeny, Flowering, and Vegetative Phenology.

Both the evolution of tree taxa and whole-genome duplication (WGD) have occurred many times during angiosperm evolution. Transcription factors are preferentially retained following WGD suggesting that functional divergence of duplicates could contribute to traits distinctive to the tree growth habit. We used gain- and loss-of-function transgenics, photoperiod treatments, and circannual expression studies in adult trees to study the diversification of three Populus FLOWERING LOCUS D-LIKE (FDL) genes encoding bZIP transcription factors. Expression patterns and transgenic studies indicate that FDL2.2 promotes flowering and that FDL1 and FDL3 function in different vegetative phenophases. Study of dominant repressor FDL versions indicates that the FDL proteins are partially equivalent in their ability to alter shoot growth. Like its paralogs, FDL3 overexpression delays short day-induced growth cessation, but also induces distinct heterochronic shifts in shoot development-more rapid phytomer initiation and coordinated delay in both leaf expansion and the transition to secondary growth in long days, but not in short days. Our results indicate that both regulatory and protein coding sequence variation contributed to diversification of FDL paralogs that has led to a degree of specialization in multiple developmental processes important for trees and their local adaptation.

Overexpression of SHORT VEGETATIVE PHASE-LIKE (SVL) in Populus delays onset and reduces abundance of flowering in field-grown trees.

The spread of transgenes and exotic germplasm from planted crops into wild or feral species is a difficult problem for public and regulatory acceptance of genetically engineered plants, particularly for wind-pollinated trees such as poplar. We report that overexpression of a poplar homolog of the floral repressor SHORT VEGETATIVE PHASE-LIKE (SVL), a homolog of the Arabidopsis MADS-box repressor SHORT VEGETATIVE PHASE (SVP), delayed the onset of flowering several years in three genotypes of field-grown transgenic poplars. Higher expression of SVL correlated with a delay in flowering onset and lower floral abundance, and did not cause morphologically obvious or statistically significant effects on leaf characteristics, tree form, or stem volume. Overexpression effects on reproductive and vegetative phenology in spring was modest and genotype-specific. Our results suggest that use of SVL and related floral repressors can be useful tools to enable a high level of containment for vegetatively propagated short-rotation woody energy or pulp crops.

Genetic containment in vegetatively propagated forest trees: CRISPR disruption of LEAFY function in Eucalyptus gives sterile indeterminate inflorescences and normal juvenile development.

Eucalyptus is among the most widely planted taxa of forest trees worldwide. However, its spread as an exotic or genetically engineered form can create ecological and social problems. To mitigate gene flow via pollen and seeds, we mutated the Eucalyptus orthologue of LEAFY (LFY) by transforming a Eucalyptus grandis × urophylla wild-type hybrid and two Flowering Locus T (FT) overexpressing (and flowering) lines with CRISPR Cas9 targeting its LFY orthologue, ELFY. We achieved high rates of elfy biallelic knockouts, often approaching 100% of transgene insertion events. Frameshift mutations and deletions removing conserved amino acids caused strong floral alterations, including indeterminacy in floral development and an absence of male and female gametes. These mutants were otherwise visibly normal and did not differ statistically from transgenic controls in juvenile vegetative growth rate or leaf morphology in greenhouse trials. Genes upstream or near to ELFY in the floral development pathway were overexpressed, whereas floral organ identity genes downstream of ELFY were severely depressed. We conclude that disruption of ELFY function appears to be a useful tool for sexual containment, without causing statistically significant or large adverse effects on juvenile vegetative growth or leaf morphology.

Vegetative phase change in Populus tremula × alba.

Plants transition through juvenile and adult phases of vegetative development in a process known as vegetative phase change (VPC). In poplars (genus Populus) the differences between these stages are subtle, making it difficult to determine when this transition occurs. Previous studies of VPC in poplars have relied on plants propagated in vitro, leaving the natural progression of this process unknown. We examined developmental morphology of seed-grown and in vitro derived Populus tremula × alba (clone 717-1B4), and compared the phenotype of these to transgenics with manipulated miR156 expression, the master regulator of VPC. In seed-grown plants, most traits changed from node-to-node during the first 3 months of development but remained constant after node 25. Many traits remained unchanged in clones over-expressing miR156, or were enhanced when miR156 was lowered, demonstrating their natural progression is regulated by the miR156/SPL pathway. The characteristic leaf fluttering of Populus is one of these miR156-regulated traits. Vegetative development in plants grown from culture mirrored that of seed-grown plants, allowing direct comparison between plants often used in research and those found in nature. These results provide a foundation for further research on the role of VPC in the ecology and evolution of this economically important genus.

RNAi of AGAMOUS genes in sweetgum alters reproductive organ identity and decreases fruit persistence.

Sweetgums (Liquidambar), members of the family Altingiaceae (Altingiales), have inflorescences and floral organs that are distinctive in structure compared with other angiosperms in which the roles of floral homeotic genes have been studied. To begin to understand the role of AGAMOUS (AG)-a floral homeotic gene that has a major role in stamen and carpel development-in development of the monosexual flowers of sweetgum, we used RNAi to reduce the expression of two members of the AG subfamily. Because AG suppression should induce floral sterility, RNAi might also provide a tool to mitigate the risks of invasiveness-and to reduce the production of its nuisance fruits or allergenic pollen-when sweetgum is used as an exotic shade or forest tree. We tested 33 independent transgenic events and non-transgenic controls during 10 years in the field. The RNAi-AG sweetgum trees maintained normal growth, phenology, and vivid fall coloration during the 10 years of study, but 8 insertion events had highly modified inflorescence and floral morphology. The modified flowers had anthers and carpels that were converted to flat leaf-like structures lacking pollen grains and ovules, respectively. The female inflorescences developed into dry papery structures that failed to produce seeds. These infructescences were smaller than control infructescences, and lost a greater percentage of biomass in a controlled decay assay. RNAi against AG genes was highly effective at impairing fertility and modifying reproductive development without significant vegetative effects in sweetgum and gave phenotypes distinct from, but similar to, that of AG loss of function in other angiosperms.

The Ectomycorrhizal Fungus Laccaria bicolor Produces Lipochitooligosaccharides and Uses the Common Symbiosis Pathway to Colonize Populus Roots.

Mycorrhizal fungi form mutualistic associations with the roots of most land plants and provide them with mineral nutrients from the soil in exchange for fixed carbon derived from photosynthesis. The common symbiosis pathway (CSP) is a conserved molecular signaling pathway in all plants capable of associating with arbuscular mycorrhizal fungi. It is required not only for arbuscular mycorrhizal symbiosis but also for rhizobia-legume and actinorhizal symbioses. Given its role in such diverse symbiotic associations, we hypothesized that the CSP also plays a role in ectomycorrhizal associations. We showed that the ectomycorrhizal fungus Laccaria bicolor produces an array of lipochitooligosaccharides (LCOs) that can trigger both root hair branching in legumes and, most importantly, calcium spiking in the host plant Populus in a CASTOR/POLLUX-dependent manner. Nonsulfated LCOs enhanced lateral root development in Populus in a calcium/calmodulin-dependent protein kinase (CCaMK)-dependent manner, and sulfated LCOs enhanced the colonization of Populus by L. bicolor Compared with the wild-type Populus, the colonization of CASTOR/POLLUX and CCaMK RNA interference lines by L. bicolor was reduced. Our work demonstrates that similar to other root symbioses, L. bicolor uses the CSP for the full establishment of its mutualistic association with Populus.

RNA interference suppression of AGAMOUS and SEEDSTICK alters floral organ identity and impairs floral organ determinacy, ovule differentiation, and seed-hair development in Populus.

The role of the floral homeotic gene AGAMOUS (AG) and its close homologues in development of anemophilous, unisexual catkins has not previously been studied. We transformed two RNA interference (RNAi) constructs, PTG and its matrix-attachment-region flanked version MPG, into the early-flowering female poplar clone 6K10 (Populus alba) to suppress the expression of its two duplicate AG orthologues. By early 2018, six out of 22 flowering PTG events and 11 out of 12 flowering MPG events showed modified floral phenotypes in a field trial in Oregon, USA. Flowers in catkins from modified events had 'carpel-inside-carpel' phenotypes. Complete disruption of seed production was observed in seven events, and sterile anther-like organs in 10 events. Events with strong co-suppression of both the two AG and two SEEDSTICK (STK) paralogues lacked both seeds and associated seed hairs. Alterations in all of the modified floral phenotypes were stable over 4 yr of study. Trees from floral-modified events did not differ significantly (P < 0.05) from nonmodified transgenic or nontransgenic controls in biomass growth or leaf morphology. AG and STK genes show strong conservation of gene function during poplar catkin development and are promising targets for genetic containment of exotic or genetically engineered trees.

Phenotypic Expression and Stability in a Large-Scale Field Study of Genetically Engineered Poplars Containing Sexual Containment Transgenes.

Genetic engineering (GE) has the potential to help meet demand for forest products and ecological services. However, high research and development costs, market restrictions, and regulatory obstacles to performing field tests have severely limited the extent and duration of field research. There is a notable paucity of field studies of flowering GE trees due to the time frame required and regulatory constraints. Here we summarize our findings from field testing over 3,300 GE poplar trees and 948 transformation events in a single, 3.6 hectare field trial for seven growing seasons; this trial appears to be the largest field-based scientific study of GE forest trees in the world. The goal was to assess a diversity of approaches for obtaining bisexual sterility by modifying RNA expression or protein function of floral regulatory genes, including LEAFY, AGAMOUS, APETALA1, SHORT VEGETATIVE PHASE, and FLOWERING LOCUS T. Two female and one male clone were transformed with up to 23 different genetic constructs designed to obtain sterile flowers or delay onset of flowering. To prevent gene flow by pollen and facilitate regulatory approval, the test genotypes chosen were incompatible with native poplars in the area. We monitored tree survival, growth, floral onset, floral abundance, pollen production, seed formation and seed viability. Tree survival was above 95%, and variation in site conditions generally had a larger impact on vegetative performance and onset of flowering than did genetic constructs. Floral traits, when modified, were stable over three to five flowering seasons, and we successfully identified RNAi or overexpression constructs that either postponed floral onset or led to sterile flowers. There was an absence of detectable somaclonal variation; no trees were identified that showed vegetative or floral modifications that did not appear to be related to the transgene added. Surveys for seedling and sucker establishment both within and around the plantation identified small numbers of vegetative shoots (root sprouts) but no seedlings, indicative of a lack of establishment of trees via seeds in the area. Overall, this long term study showed that GE containment traits can be obtained which are effective, stable, and not associated with vegetative abnormalities or somaclonal variation.

Variation in Mutation Spectra Among CRISPR/Cas9 Mutagenized Poplars.

In an effort to produce reliably contained transgenic trees, we used the CRISPR/Cas9 system to alter three genes expected to be required for normal flowering in poplar (genus Populus). We designed synthetic guide RNAs (sgRNAs) to target the poplar homolog of the floral meristem identity gene, LEAFY (LFY), and the two poplar orthologs of the floral organ identity gene AGAMOUS (AG). We generated 557 transgenic events with sgRNA(s) and the Cas9 transgene and 49 events with Cas9 but no sgRNA, and analyzed all events by Sanger Sequencing of both alleles. Out of the 684 amplicons from events with sgRNAs, 474 had mutations in both alleles (77.5%). We sequenced both AG paralogs for 71 events in INRA clone 717-1B4 and 22 events in INRA clone 353-53, and found that 67 (94.4%) and 21 (95.5%) were double locus knockouts. Due partly to a single nucleotide polymorphism (SNP) present in the target region, one sgRNA targeting the AG paralogs was found to be completely inactive by itself (0%) but showed some activity in generating deletions when used in a construct with a second sgRNA (10.3-24.5%). Small insertion/deletion (indel) mutations were prevalent among mutated alleles of events with only one sgRNA (ranging from 94.3 to 99.1%), while large deletions were prevalent among alleles with two active sgRNAs (mean proportion of mutated alleles was 22.6% for small indels vs. 77.4% for large indels). For both LFY and AG, each individual sgRNA-gene combination had a unique mutation spectrum (p < 0.001). An AG-sgRNA construct with two sgRNAs had similar mutation spectra among two poplar clones (p > 0.05), however, a LFY-sgRNA construct with a single sgRNA gave significantly different mutation spectra among the same two clones (p < 0.001). The 49 empty vector control events had no mutations in either allele, and 310 potential "off-target" sequences also had no mutations in 58 transgenic events studied. CRISPR/Cas9 is a very powerful and precise system for generating loss-of-function mutations in poplars, and should be effective for generating reliably infertile trees that may promote regulatory, market, or public acceptance of genetic engineering technology.

BIG LEAF is a regulator of organ size and adventitious root formation in poplar.

Here we report the discovery through activation tagging and subsequent characterization of the BIG LEAF (BL) gene from poplar. In poplar, BL regulates leaf size via positively affecting cell proliferation. Up and downregulation of the gene led to increased and decreased leaf size, respectively, and these phenotypes corresponded to increased and decreased cell numbers. BL function encompasses the early stages of leaf development as native BL expression was specific to the shoot apical meristem and leaf primordia and was absent from the later stages of leaf development and other organs. Consistently, BL downregulation reduced leaf size at the earliest stages of leaf development. Ectopic expression in mature leaves resulted in continued growth most probably via sustained cell proliferation and thus the increased leaf size. In contrast to the positive effect on leaf growth, ectopic BL expression in stems interfered with and significantly reduced stem thickening, suggesting that BL is a highly specific activator of growth. In addition, stem cuttings from BL overexpressing plants developed roots, whereas the wild type was difficult to root, demonstrating that BL is a positive regulator of adventitious rooting. Large transcriptomic changes in plants that overexpressed BL indicated that BL may have a broad integrative role, encompassing many genes linked to organ growth. We conclude that BL plays a fundamental role in control of leaf size and thus may be a useful tool for modifying plant biomass productivity and adventitious rooting.

FT overexpression induces precocious flowering and normal reproductive development in Eucalyptus.

Eucalyptus trees are among the most important species for industrial forestry worldwide. However, as with most forest trees, flowering does not begin for one to several years after planting which can limit the rate of conventional and molecular breeding. To speed flowering, we transformed a Eucalyptus grandis × urophylla hybrid (SP7) with a variety of constructs that enable overexpression of FLOWERING LOCUS T (FT). We found that FT expression led to very early flowering, with events showing floral buds within 1-5 months of transplanting to the glasshouse. The most rapid flowering was observed when the cauliflower mosaic virus 35S promoter was used to drive the Arabidopsis thaliana FT gene (AtFT). Early flowering was also observed with AtFT overexpression from a 409S ubiquitin promoter and under heat induction conditions with Populus trichocarpa FT1 (PtFT1) under control of a heat-shock promoter. Early flowering trees grew robustly, but exhibited a highly branched phenotype compared to the strong apical dominance of nonflowering transgenic and control trees. AtFT-induced flowers were morphologically normal and produced viable pollen grains and viable self- and cross-pollinated seeds. Many self-seedlings inherited AtFT and flowered early. FT overexpression-induced flowering in Eucalyptus may be a valuable means for accelerating breeding and genetic studies as the transgene can be easily segregated away in progeny, restoring normal growth and form.

EARLY BUD-BREAK 1 (EBB1) is a regulator of release from seasonal dormancy in poplar trees.

Trees from temperate latitudes transition between growth and dormancy to survive dehydration and freezing stress during winter months. We used activation tagging to isolate a dominant mutation affecting release from dormancy and identified the corresponding gene EARLY BUD-BREAK 1 (EBB1). We demonstrate through positioning of the tag, expression analysis, and retransformation experiments that EBB1 encodes a putative APETALA2/Ethylene responsive factor transcription factor. Transgenic up-regulation of the gene caused early bud-flush, whereas down-regulation delayed bud-break. Native EBB1 expression was highest in actively growing apices, undetectable during the dormancy period, but rapidly increased before bud-break. The EBB1 transcript was localized in the L1/L2 layers of the shoot meristem and leaf primordia. EBB1-overexpressing transgenic plants displayed enlarged shoot meristems, open and poorly differentiated buds, and a higher rate of cell division in the apex. Transcriptome analyses of the EBB1 transgenics identified 971 differentially expressed genes whose expression correlated with the EBB1 expression changes in the transgenic plants. Promoter analysis among the differentially expressed genes for the presence of a canonical EBB1-binding site identified 65 putative target genes, indicative of a broad regulatory context of EBB1 function. Our results suggest that EBB1 has a major and integrative role in reactivation of meristem activity after winter dormancy.

Methylome reorganization during in vitro dedifferentiation and regeneration of Populus trichocarpa.

BACKGROUND: Cytosine DNA methylation (5mC) is an epigenetic modification that is important to genome stability and regulation of gene expression. Perturbations of 5mC have been implicated as a cause of phenotypic variation among plants regenerated through in vitro culture systems. However, the pattern of change in 5mC and its functional role with respect to gene expression, are poorly understood at the genome scale. A fuller understanding of how 5mC changes during in vitro manipulation may aid the development of methods for reducing or amplifying the mutagenic and epigenetic effects of in vitro culture and plant transformation. RESULTS: We investigated the in vitro methylome of the model tree species Populus trichocarpa in a system that mimics routine methods for regeneration and plant transformation in the genus Populus (poplar). Using methylated DNA immunoprecipitation followed by high-throughput sequencing (MeDIP-seq), we compared the methylomes of internode stem segments from micropropagated explants, dedifferentiated calli, and internodes from regenerated plants. We found that more than half (56%) of the methylated portion of the genome appeared to be differentially methylated among the three tissue types. Surprisingly, gene promoter methylation varied little among tissues, however, the percentage of body-methylated genes increased from 9% to 14% between explants and callus tissue, then decreased to 8% in regenerated internodes. Forty-five percent of differentially-methylated genes underwent transient methylation, becoming methylated in calli, and demethylated in regenerants. These genes were more frequent in chromosomal regions with higher gene density. Comparisons with an expression microarray dataset showed that genes methylated at both promoters and gene bodies had lower expression than genes that were unmethylated or only promoter-methylated in all three tissues. Four types of abundant transposable elements showed their highest levels of 5mC in regenerated internodes. CONCLUSIONS: DNA methylation varies in a highly gene- and chromosome-differential manner during in vitro differentiation and regeneration. 5mC in redifferentiated tissues was not reset to that in original explants during the study period. Hypermethylation of gene bodies in dedifferentiated cells did not interfere with transcription, and may serve a protective role against activation of abundant transposable elements.