To better understand the molecular control of leaf senescence, we examined transcriptome changes during seasonal leaf senescence in Populus trichocarpa Nisqually-1, the Populus reference genome, growing in its natural habitat. Using monthly (from May to October) transcriptomes for three years (2009, 2015, and 2016), we identified 17,974 differentially expressed genes (DEGs; false discovery rate <0.05; log-fold change cutoff = 0) from 36,007 expressed Populus gene models. A total of 14,415 DEGs were directly related to transitions between four major developmental phases – growth, senescence initiation, reorganization, and senescence termination. These DEGs were significantly ( p < 0.05) enriched in 279 gene ontology (GO) terms, including those related to photosynthesis, metabolic process, catalytic activity, protein phosphorylation, kinase activity, pollination, and transport. Also, there were 881 differentially expressed transcription factor (TF) genes from 54 TF families, notably bHLH, MYB, ERF, MYB-related, NAC, and WRKY. We also examined 28 DEGs known as alternative splicing (AS) factors that regulate AS process, and found evidence for a reduced level of AS activity during leaf senescence. Furthermore, we were able to identify a number of promoter sequence motifs associated with leaf senescence. This work provides a comprehensive resource for identification of genes involved in seasonal leaf senescence in trees, and informs efforts to explore the conservation and divergence of molecular mechanisms underlying leaf senescence between annual and perennial species.

Plant Reactome (http://plantreactome.gramene.org/) is a free, open-source, curated plant pathway database portal, provided as part of the Gramene project. The database provides intuitive bioinformatics tools for the visualization, analysis and interpretation of pathway knowledge to support genome annotation, genome analysis, modeling, systems biology, basic research and education. Plant Reactome employs the structural framework of a plant cell to show metabolic, transport, genetic, developmental and signaling pathways. We manually curate molecular details of pathways in these domains for reference species Oryza sativa (rice) supported by published literature and annotation of well-characterized genes. Two hundred twenty-two rice pathways, 1025 reactions associated with 1173 proteins, 907 small molecules and 256 literature references have been curated to date. These reference annotations were used to project pathways for 62 model, crop and evolutionarily significant plant species based on gene homology. Database users can search and browse various components of the database, visualize curated baseline expression of pathway-associated genes provided by the Expression Atlas and upload and analyze their Omics datasets. The database also offers data access via Application Programming Interfaces (APIs) and in various standardized pathway formats, such as SBML and BioPAX.

Gene flow from recombinant-DNA-modified (GMO) trees is a major barrier to their public acceptance and regulatory approval. Because many intensively grown trees are vegetatively propagated, complete sexual sterility could be a powerful means to mitigate or prevent gene flow. We tested four pairs of zinc-finger nucleases (ZFNs) as mutagenic agents against the LEAFY and AGAMOUS orthologs in poplar that are expected to be required for sexual fertility. To reduce the potential for pleiotropic effects from mutagenesis, each of the pairs was functionally linked to a heat shock promoter to provide inducible ZFN expression. Using Agrobacterium tumefaciens, we transformed more than 21,000 total explants compromised of both male and female hybrid poplar. The rate of transformation for the ZFN constructs (2 %) was generally reduced compared to the transgenic control (8 %). We produced 391 ZFN transgenic shoots of which only two developed into plants with mutations in a target gene; both were 7-bp deletions in one allele of the PtAG2 locus. No mutations were observed in the PtAG1 or PtLFY loci. Our results indicate a mutation rate of zero to 0.3 % per explant per allele, among the lowest reported for ZFN mutagenesis in plants. The combined effects of low recovery of transgenic plants, a modest mutation frequency, and much higher reported rates of directed mutation for other gene editing methods suggest that the efficient use of ZFNs in poplar requires further technical improvements.

Gramene (http://www.gramene.org) is an online resource for comparative functional genomics in crops and model plant species. Its two main frameworks are genomes (collaboration with Ensembl Plants) and pathways (The Plant Reactome and archival BioCyc databases). Since our last NAR update, the database website adopted a new Drupal management platform. The genomes section features 39 fully assembled reference genomes that are integrated using ontology-based annotation and comparative analyses, and accessed through both visual and programmatic interfaces. Additional community data, such as genetic variation, expression and methylation, are also mapped for a subset of genomes. The Plant Reactome pathway portal (http://plantreactome.gramene.org) provides a reference resource for analyzing plant metabolic and regulatory pathways. In addition to ∼200 curated rice reference pathways, the portal hosts gene homology-based pathway projections for 33 plant species. Both the genome and pathway browsers interface with the EMBL-EBI's Expression Atlas to enable the projection of baseline and differential expression data from curated expression studies in plants. Gramene's archive website (http://archive.gramene.org) continues to provide previously reported resources on comparative maps, markers and QTL. To further aid our users, we have also introduced a live monthly educational webinar series and a Gramene YouTube channel carrying video tutorials.

As a step toward functional annotation of genes required for floral initiation and development within the Eucalyptus genome, we used short read sequencing to analyze transcriptomes of floral buds from early and late developmental stages, and compared these with transcriptomes of diverse vegetative tissues, including leaves, roots, and stems. A subset of 4807 genes (13% of protein-coding genes) were differentially expressed between floral buds of either stage and vegetative tissues. A similar proportion of genes were differentially expressed among all tissues. A total of 479 genes were differentially expressed between early and late stages of floral development. Gene function enrichment identified 158 gene ontology classes that were overrepresented in floral tissues, including 'pollen development' and 'aromatic compound biosynthetic process'. At least 40 floral-dominant genes lacked functional annotations and thus may be novel floral transcripts. We analyzed several genes and gene families in depth, including 49 putative biomarkers of floral development, the MADS-box transcription factors, 'S-domain'-receptor-like kinases, and selected gene family members with phosphatidylethanolamine-binding protein domains. Expanded MADS-box gene subfamilies in Eucalyptus grandis included SUPPRESSOR OF OVEREXPRESSION OF CO 1 (SOC1), SEPALLATA (SEP) and SHORT VEGETATIVE PHASE (SVP) Arabidopsis thaliana homologs. These data provide a rich resource for functional and evolutionary analysis of genes controlling eucalypt floral development, and new tools for breeding and biotechnology.

Eucalypts are the world's most widely planted hardwood trees. Their outstanding diversity, adaptability and growth have made them a global renewable resource of fibre and energy. We sequenced and assembled >94% of the 640-megabase genome of Eucalyptus grandis. Of 36,376 predicted protein-coding genes, 34% occur in tandem duplications, the largest proportion thus far in plant genomes. Eucalyptus also shows the highest diversity of genes for specialized metabolites such as terpenes that act as chemical defence and provide unique pharmaceutical oils. Genome sequencing of the E. grandis sister species E. globulus and a set of inbred E. grandis tree genomes reveals dynamic genome evolution and hotspots of inbreeding depression. The E. grandis genome is the first reference for the eudicot order Myrtales and is placed here sister to the eurosids. This resource expands our understanding of the unique biology of large woody perennials and provides a powerful tool to accelerate comparative biology, breeding and biotechnology.

Allelic variation at the Cu–Zn superoxide dismutase (SOD1) locus has been shown to be associated with resistance of the snail, Biomphalaria glabrata, to infection by the trematode parasite, Schistosoma mansoni. SOD1 catalyses the production of hydrogen peroxide, a known cytotoxic component of the oxidative burst used in defence against pathogens. In our laboratory population of B. glabrata, the most resistant allele at SOD1 is over-expressed relative to the other two alleles. Because hydrogen peroxide also causes oxidative stress on host tissues, we hypothesised that over-expression of SOD1 might be compensated by epistatic interactions with other loci involved in oxidation–reduction (redox) pathways. Catalase, peroxiredoxins and glutathione peroxidases all degrade hydrogen peroxide. We tested whether alleles at each of these loci were in linkage disequilibrium with SOD1 in our population, as might be expected given strong epistatic selection. We found that SOD1, catalase (CAT) and a peroxiredoxin locus (PRX4) are in strong linkage disequilibrium in our population. We also found that these loci are tightly linked, within 1–2cM of each other, which explains the high linkage disequilibrium. This result raises the possibility that there is a linked cluster of redox genes, and perhaps other defence-relevant genes, in the B. glabrata genome. Whether epistatic interactions for fitness actually exist among these loci still needs to be tested. However the close physical linkage among SOD1, PRX4 and CAT, and subsequent high disequilibrium, makes such interactions a plausible hypothesis.

Background Functional annotations of large plant genome projects mostly provide information on gene function and gene families based on the presence of protein domains and gene homology, but not necessarily in association with gene expression or metabolic and regulatory networks. These additional annotations are necessary to understand the physiology, development and adaptation of a plant and its interaction with the environment. Results RiceCyc is a metabolic pathway networks database for rice. It is a snapshot of the substrates, metabolites, enzymes, reactions and pathways of primary and intermediary metabolism in rice. RiceCyc version 3.3 features 316 pathways and 6,643 peptide-coding genes mapped to 2,103 enzyme-catalyzed and 87 protein-mediated transport reactions. The initial functional annotations of rice genes with InterPro, Gene Ontology, MetaCyc, and Enzyme Commission (EC) numbers were enriched with annotations provided by KEGG and Gramene databases. The pathway inferences and the network diagrams were first predicted based on MetaCyc reference networks and plant pathways from the Plant Metabolic Network, using the Pathologic module of Pathway Tools. This was enriched by manually adding metabolic pathways and gene functions specifically reported for rice. The RiceCyc database is hierarchically browsable from pathway diagrams to the associated genes, metabolites and chemical structures. Through the integrated tool OMICs Viewer, users can upload transcriptomic, proteomic and metabolomic data to visualize expression patterns in a virtual cell. RiceCyc, along with additional species-specific pathway databases hosted in the Gramene project, facilitates comparative pathway analysis. Conclusions Here we describe the RiceCyc network development and discuss its contribution to rice genome annotations. As a case study to demonstrate the use of RiceCyc network as a discovery environment we carried out an integrated bioinformatic analysis of rice metabolic genes that are differentially regulated under diurnal photoperiod and biotic stress treatments. The analysis of publicly available rice transcriptome datasets led to the hypothesis that the complete tryptophan biosynthesis and its dependent metabolic pathways including serotonin biosynthesis are induced by taxonomically diverse pathogens while also being under diurnal regulation. The RiceCyc database is available online for free access at http://www.gramene.org/pathway/ .

Gramene (http://www.gramene.org) is a curated online resource for comparative functional genomics in crops and model plant species, currently hosting 27 fully and 10 partially sequenced reference genomes in its build number 38. Its strength derives from the application of a phylogenetic framework for genome comparison and the use of ontologies to integrate structural and functional annotation data. Whole-genome alignments complemented by phylogenetic gene family trees help infer syntenic and orthologous relationships. Genetic variation data, sequences and genome mappings available for 10 species, including Arabidopsis, rice and maize, help infer putative variant effects on genes and transcripts. The pathways section also hosts 10 species-specific metabolic pathways databases developed in-house or by our collaborators using Pathway Tools software, which facilitates searches for pathway, reaction and metabolite annotations, and allows analyses of user-defined expression datasets. Recently, we released a Plant Reactome portal featuring 133 curated rice pathways. This portal will be expanded for Arabidopsis, maize and other plant species. We continue to provide genetic and QTL maps and marker datasets developed by crop researchers. The project provides a unique community platform to support scientific research in plant genomics including studies in evolution, genetics, plant breeding, molecular biology, biochemistry and systems biology.

A framework for understanding the synthesis and catalysis of metabolites and other biochemicals by proteins is crucial for unraveling the physiology of cells. To create such a framework for Zea mays L. subsp. mays (maize), we developed MaizeCyc, a metabolic network of enzyme catalysts, proteins, carbohydrates, lipids, amino acids, secondary plant products, and other metabolites by annotating the genes identified in the maize reference genome sequenced from the B73 variety. MaizeCyc version 2.0.2 is a collection of 391 maize pathways involving 8889 enzyme mapped to 2110 reactions and 1468 metabolites. We used MaizeCyc to describe the development and function of maize organs including leaf, root, anther, embryo, and endosperm by exploring the recently published microarray‐based maize gene expression atlas. We found that 1062 differentially expressed metabolic genes mapped to 524 unique enzymatic reactions associated with 310 pathways. The MaizeCyc pathway database was created by running a library of evidences collected from the maize genome annotation, gene‐based phylogeny trees, and comparison to known genes and pathways from rice (Oryza sativa L.) and Arabidopsis thaliana (L.) Heynh. against the PathoLogic module of Pathway Tools. The network and the database that were also developed as a community resource are freely accessible online at http://maizecyc.maizegdb.org to facilitate analysis and promote studies on metabolic genes in maize.

When salmonid fish that have been raised in hatcheries spawn in the wild, they often produce fewer surviving adult offspring than wild fish. Recent data from steelhead (Oncorhynchus mykiss) in the Hood River (Oregon, USA) show that even one or two generations of hatchery culture can result in dramatic declines in fitness. Although intense domestication selection could cause such declines, it is worth considering alternative explanations. One possibility is heritable epigenetic changes induced by the hatchery environment. Here, we show, using methylation-sensitive amplified fragment length polymorphism, that hatchery and wild adult steelhead from the Hood River do not appear to differ substantially in overall levels of genomic methylation. Thus, although altered methylation of specific DNA sites or other epigenetic processes could still be important, the hatchery environment does not appear to cause a global hypo- or hypermethylation of the genome or create a large number of sites that are differentially methylated.

The microsatellite motifs AG, AC, and ATG were found to be the most abundant in Douglas-fir (Pseudotsuga menziesii (Mirb.) Franco) and several other conifer tree species among di-, tri-, and tetra-nucleotide simple sequence repeats (SSR). Colonies containing AG, AC, and ATG repeats were selected from enriched genomic libraries of Douglas-fir, and 603 were sequenced. Polymerase chain reaction (PCR) primers were designed from flanking sequences in 102 of the SSR clones, of which 50 primer pairs (for 10 AC-repeat microsatellites and 40 AG-repeat microsatellites) produced robust amplification products. Variability was confirmed with 24 unrelated Douglas-fir trees and Medelian segregation with 33-66 progeny from 3 full-sib populations. Forty-eight of the 50 loci were polymorphic, with a mean of 7.5 alleles per locus. Allele sizes ranged from 73 to 292 base pairs. Allele frequencies for the 48 polymorphic loci varied from 0.017 to 0.906 with mean allele frequency of 0.250. Expected heterozygosities among the polymorphic loci varied from 0.174 to 0.926, with a mean of 0.673. Additional, high molecular weight PCR products were amplified by some of the primer pairs, but they did not interfere with the scoring of alleles. Most of the Douglas-fir primer pairs also amplified SSR-containing loci in other conifer species.

Polymorphic microsatellite loci were developed for Western Hemlock [Tsuga heterophylla (Raf.) Sarg], a prominent forest tree species in Western North America. Microsatellite-enriched libraries were screened for (CA)(n) dinucleotide repeats from which 33 positive clones were sequenced. Polymerase chain reaction (PCR) primers for 16 microsatellite loci were prepared and tested against DNA from unrelated Western Hemlock trees. The 12 most informative microsatellite loci are reported here. From four to 22 alleles per locus were observed, with an average expected heterozygousity of 0.799.

Douglas-fir [Pseudotsuga menziesii(Mirb.) France] is one of the world's most valuable timber species and is widely used in reforestation. Therefore there is considerable interest in understanding its genetic structure, In conifers approximately 75% of the genome is composed of repetitive DNA. Thus for conifers characterization of repetitive DNA is a significant part of genome analysis. We have characterized the organization of 5S rRNA genes in Douglas-fir at both the molecular and chromosome levels, 5S DNA repeat units containing the coding sequence for 5S rRNA and the nontranscribed spacer (NTS) were cloned using PCR. Sequencing and Southern hybridization revealed repeat units of 888 and 871 bp in length, the latter with a 17 bp deletion in the NTS, The coding region showed high homology with other eukaryotic 5S rRNA genes. A 35 bp region of the NTS immediately upstream of the 5' end of the coding region showed high similarity to other conifers but not to other published plant 5S rDNA sequences. Physical mapping of 5S rDNA by fluorescent in situ hybridization using a biotinilated homologous probe revealed a single subtelomeric site on one pair of large metacentric chromosomes.

Simple sequence repeats have emerged as the current genetic marker of choice due to their co-dominant, multiallelic, highly polymorphic nature and reproducibility. Therefore as the initial component in a program designed to produce genetic markers for selection of economically important quantitative traits in Douglas-fir and western hemlock we have undertaken an intensive effort to characterize a large number of SSR loci. Hybridization of all 47 possible di, tri and tetra nucleotide SSRs against eight economically important conifer species including the two species of interest revealed AG and ATG motifs to be the most abundant in most of the conifer species studied. Thus far we have developed 95 primer pairs for Douglas-fir and western hemlock and have sequence characterized another 400 SSR containing clones. Most of our primer pairs amplified a product of expected size range. Amplified products were verified by sequence characterization and Southern hybridization.

Putrescine, spermidine, and spermine levels during somatic embryogenesis of interior spruce (Picea glauca x Picea engelmannii complex) were quantified. On abscisic acid supplemented growth medium putrescine and spermidine levels increased two-fold coinciding with maturation of the early somatic embryos to globular embryos. Polyclonal antibodies raised against Escherichia coli arginine decarboxylase (ADC) and ornithine decarboxylase (ODC), following affinity purification specifically recognized spruce ADC and ODC, which corresponded to 85kD and 65kD bands on western blots of total protein ex tracts from embryogenic masses. Immunoassays using these antibodies showed increased ADC levels corresponding to embryo maturation while ODC levels remained the same. From these results it is concluded that polyamines are involved in the maturation of somatic embryos of interior spruce.

Polyamines were localized in embryogenic callus cells of Picea glauca by the fluorescent dye o-phthalaldehyde and polyclonal antibodies raised against putrescine, spermidine, and spermine. Localization by both methods showed higher levels of polyamines in the nuclei and nucleoli than in the cytoplasm. In contrast with previous reports on cytochemical localization of polyamines in animal cells, stain was always excluded from the condensed chromatin of spruce somatic embryo cells. Key words: polyamines, cytochemistry, immunolocalization, white spruce.

We have localized the major ribosomal DNA (rDNA) loci on metaphase chromosomes and in interphase nuclei of white spruce (2n = 24) by fluorescence in situ hybridization. Hybridization sites.of the biotin-labelled rDNA probe were detected using antibody-fluorochrome conjugates and a confocal laser scanning microscope. White spruce has at least 12, and possibly as many as 14, rDNA sites, 1 site present on each of seven separate chromosome pairs. This is one of the highest numbers of rDNA loci yet reported among plant species. The position of the rDNA loci together with secondary constriction patterns permit, for the first time, all homologous pairs of white spruce chromosomes to be distinguished. We discuss the application of molecular cytogenetics in studies relating to the organization and evolution of DNA sequences within conifer genomes.

Trichome morphology of 136 species of Cuphea was studied using light and electron microscopy. Seven different trichome types were recognized, and three were found to have distinct taxonomic value at sectional and subsectional levels. Biarmed trichomes characterize the sections Heteranthus and Amazoniana. Flattened hyaline hairs define a group of species in section Melvilla subsection Pachycalyx. Globose, sessile, glandular trichomes are restricted to sections Archocuphea and Melicyathium. Other types are widely distributed in the genus and can be useful diagnostic characters at the species level. TEM studies of the glandular trichomes revealed subcellular structures thought to be involved in the secretion of resin.