The mitochondrial genome sequence of Abies alba Mill. reveals a high structural and combinatorial variation.

BACKGROUND: Plant mitogenomes vary widely in size and genomic architecture. Although hundreds of plant mitogenomes of angiosperm species have already been sequence-characterized, only a few mitogenomes are available from gymnosperms. Silver fir (Abies alba) is an economically important gymnosperm species that is widely distributed in Europe and occupies a large range of environmental conditions. Reference sequences of the nuclear and chloroplast genome of A. alba are available, however, the mitogenome has not yet been assembled and studied. RESULTS: Here, we used paired-end Illumina short reads generated from a single haploid megagametophyte in combination with PacBio long reads from high molecular weight DNA of needles to assemble the first mitogenome sequence of A. alba. Assembly and scaffolding resulted in 11 mitogenome scaffolds, with the largest scaffold being 0.25 Mbp long. Two of the scaffolds displayed a potential circular structure supported by PCR. The total size of the A. alba mitogenome was estimated at 1.43 Mbp, similar to the size (1.33 Mbp) of a draft assembly of the Abies firma mitogenome. In total, 53 distinct genes of known function were annotated in the A. alba mitogenome, comprising 41 protein-coding genes, nine tRNA, and three rRNA genes. The proportion of highly repetitive elements (REs) was 0.168. The mitogenome seems to have a complex and dynamic structure featured by high combinatorial variation, which was specifically confirmed by PCR for the contig with the highest mapping coverage. Comparative analysis of all sequenced mitogenomes of gymnosperms revealed a moderate, but significant positive correlation between mitogenome size and proportion of REs. CONCLUSIONS: The A. alba mitogenome provides a basis for new comparative studies and will allow to answer important structural, phylogenetic and other evolutionary questions. Future long-read sequencing with higher coverage of the A. alba mitogenome will be the key to further resolve its physical structure. The observed positive correlation between mitogenome size and proportion of REs will be further validated once available mitogenomes of gymnosperms would become more numerous. To test whether a higher proportion of REs in a mitogenome leads to an increased recombination and higher structural complexity and variability is a prospective avenue for future research.

Tree Genetic Engineering, Genome Editing and Genomics.

In this Special Issue [...].

European oak chemical diversity - from ecotypes to herbivore resistance.

Climate change is increasing insect pressure and forcing plants to adapt. Although chemotypic differentiation and phenotypic plasticity in spatially separated tree populations are known for decades, understanding their importance in herbivory resistance across forests remains challenging. We studied four oak forest stands in Germany using nontarget metabolomics, elemental analysis, and chemometrics and mapped the leaf metabolome of herbivore-resistant (T-) and herbivore-susceptible (S-) European oaks (Quercus robur) to Tortrix viridana, an oak pest that causes severe forest defoliation. Among the detected metabolites, we identified reliable metabolic biomarkers to distinguish S- and T-oak trees. Chemotypic differentiation resulted in metabolic shifts of primary and secondary leaf metabolism. Across forests, T-oaks allocate resources towards constitutive chemical defense enriched of polyphenolic compounds, e.g. the flavonoids kaempferol, kaempferol and quercetin glucosides, while S-oaks towards growth-promoting substances such as carbohydrates and amino-acid derivatives. This extensive work across natural forests shows that oaks' resistance and susceptibility to herbivory are linked to growth-defense trade-offs of leaf metabolism. The discovery of biomarkers and the developed predictive model pave the way to understand Quercus robur's susceptibility to herbivore attack and to support forest management, contributing to the preservation of oak forests in Europe.

Transcriptome analysis of North American sweet birch (Betula lenta) revealed a higher expression of genes involved in the biosynthesis of secondary metabolites than European silver birch (B. pendula).

The North American Betula lenta L. (sweet birch) has been used for medicinal reasons for centuries by native Americans. Although sophisticated technologies have rapidly been developed, a large information gap has been observed regarding genetic regulators of medicinally important compounds in sweet birch. Very little is known on the different genes involved in secondary metabolic biosynthesis in sweet birch. To gain a deeper insight into genetic factors, we performed a transcriptome analysis of each three biological samples from different independent trees of sweet and European silver birch (B. pendula Roth). This allowed us to precisely quantify the transcripts of about 24,000 expressed genes including 29 prominent candidate genes putatively involved in the biosynthesis of secondary metabolites like terpenoids, and aromatic benzoic acids. A total number of 597 genes were differentially expressed between B. lenta and B. pendula, while 264 and 210 genes showed upregulation in the bark and leaf of B. lenta, respectively. Moreover, we identified 39 transcriptional regulatory elements, involved in secondary metabolite biosynthesis, upregulated in B. lenta. Our study demonstrated the potential of RNA sequencing to identify candidate genes interacting in secondary metabolite biosynthesis in sweet birch. The candidate genes identified in this study could be subjected to genetic engineering to functionally characterize them in sweet birch. This knowledge can be beneficial to the increase of therapeutically important compounds.

Sequencing of two transgenic early-flowering poplar lines confirmed vector-free single-locus T-DNA integration.

Next-generation sequencing (NGS) approaches are attractive alternatives to the PCR-based characterisation of genetically modified plants for safety assessment and labelling since NGS is highly sensitive to the detection of T-DNA inserts as well as vector backbone sequences in transgenic plants. In this study, two independent transgenic male Populus tremula lines, T193-2 and T195-1, both carrying the FLOWERING LOCUS T gene from Arabidopsis thaliana under control of a heat-inducible promoter (pHSP::AtFT) and the non-transgenic control clone W52, were further characterised by NGS and third-generation sequencing. The results support previous findings that the T-DNA was hemizygously inserted in one genomic locus of each line. However, the T-DNA insertions consist of conglomerations of one or two T-DNA copies together with a small T-DNA fragment without AtFT parts. Based on NGS data, no additional T-DNA splinters or vector backbone sequences could be identified in the genome of the two transgenic lines. Seedlings derived from crosses between the pHSP::AtFT transgenic male parents and female wild type plants are therefore expected to be T-DNA splinter or vector backbone free. Thus, PCR analyses amplifying a partial T-DNA fragment with AtFT-specific primers are sufficient to determine whether the seedlings are transgenic or not. An analysis of 72 second generation-seedlings clearly showed that about 50% of them still reveal the presence of the T-DNA, confirming data already published. To prove if unanticipated genomic changes were induced by T-DNA integration, extended future studies using long-range sequencing technologies are required once a suitable chromosome-level P. tremula reference genome sequence is available.

RNA-seq of eight different poplar clones reveals conserved up-regulation of gene expression in response to insect herbivory.

BACKGROUND: Herbivorous insects can have a profound impact on plant growth performance. In some years, canopy damage in poplar plantations exceeds 50% of the total leaf surface, thereby possibly compromising carbon fixation and biomass yield. To assess the transcriptional response of elite poplar clones to insect feeding and to test whether this response varies between different genotypes, we performed an RNA-sequencing experiment. We deeply sequenced the transcriptomes of eight elite clones belonging to three poplar species (Populus trichocarpa, P. nigra and P. maximowiczii), under Phratora vitellinae feeding and control conditions. This allowed us to precisely quantify transcript levels of about 24,000 expressed genes. RESULTS: Our data reveal a striking overall up-regulation of gene expression under insect attack in all eight poplar clones studied. The up-regulated genes were markedly enriched for the biological process 'regulation of transcription' indicating a highly concerted restructuring of the transcriptome. A search for potential cis-regulatory elements (CREs) that may be involved in this process identified the G-box (CACGTG) as the most significant motif in the promoters of the induced genes. In line with the role of the G-box in jasmonate (JA)-mediated activation of gene expression by MYC2, several genes involved in JA biosynthesis and signaling were up-regulated in our dataset. A co-expression network analysis additionally highlighted WRKY transcription factors. Within the most prominent expression module, WRKYs were strongly overrepresented and occupied several network hubs. Finally, the insect-induced genes comprised several protein families known to be involved in plant defenses, e.g. cytochrome P450s, chitinases and protease inhibitors. CONCLUSIONS: Our data represent a comprehensive characterization of the transcriptional response of selected elite poplar clones to insect herbivory. Our results suggest that the concerted up-regulation of gene expression is controlled by JA signaling and WRKY transcription factors, and activates several defense mechanisms. Our data highlight potential targets of selection and may thus contribute to breeding insect-resistant poplar clones.

Knockdown of PCBER1, a gene of neolignan biosynthesis, resulted in increased poplar growth.

MAIN CONCLUSION: Poplar trees displayed an increased plant height due to the transgenic knockdown of PCBER1, a gene of lignan biosynthesis. The wood composition was slightly altered in both overexpression and knockdown lines. The gene PHENYLCOUMARAN BENZYLIC ETHER REDUCTASE1 (PCBER1) is well known as an important gene in the synthesis of lignans, a group of diverse phenylpropanoid derivatives. They are widely distributed in the plant kingdom and may have a role in both plant defense and growth regulation. To analyze its role in biomass formation and wood composition in poplar, both overexpression and knockdown approaches have been performed. Transgenic lines were analyzed on genetic and phenotypic levels, and partly in regard to their biomass composition. While the PCBER1 overexpression approach remained unremarkable concerning the plant height, biomass composition of obtained transgenic lines was modified. They had a significantly increased amount of ethanol extractives. The PCBER1 knockdown resulted in significantly deviating plants; after 17 months of greenhouse cultivation, transgenic plants were up to 38% higher compared to non-transgenic wild type. Most examined transgenic lines did not reveal a significantly enhanced stem diameter after three vegetation periods in the greenhouse. Significant changes were not obtained with regard to the three major wood components, lignin, cellulose and hemicelluloses. As a slight but not significant reduction in ethanol extractives was detected, the hypothesis arises that the lignan content could be influenced. Lignans become important in the pharmaceutical industry and clinical studies concerning cancer and other diseases, thus further investigations on lignan formation in poplar and its connection to biomass formation seem promising.

Complete Chloroplast Genome Sequences of Four Meliaceae Species and Comparative Analyses.

The Meliaceae family mainly consists of trees and shrubs with a pantropical distribution. In this study, the complete chloroplast genomes of four Meliaceae species were sequenced and compared with each other and with the previously published Azadirachta indica plastome. The five plastomes are circular and exhibit a quadripartite structure with high conservation of gene content and order. They include 130 genes encoding 85 proteins, 37 tRNAs and 8 rRNAs. Inverted repeat expansion resulted in a duplication of rps19 in the five Meliaceae species, which is consistent with that in many other Sapindales, but different from many other rosids. Compared to Azadirachta indica, the four newly sequenced Meliaceae individuals share several large deletions, which mainly contribute to the decreased genome sizes. A whole-plastome phylogeny supports previous findings that the four species form a monophyletic sister clade to Azadirachta indica within the Meliaceae. SNPs and indels identified in all complete Meliaceae plastomes might be suitable targets for the future development of genetic markers at different taxonomic levels. The extended analysis of SNPs in the matK gene led to the identification of four potential Meliaceae-specific SNPs as a basis for future validation and marker development.

Integrated transcriptomics and metabolomics decipher differences in the resistance of pedunculate oak to the herbivore Tortrix viridana L.

BACKGROUND: The interaction between insect pests and their host plants is a never-ending race of evolutionary adaption. Plants have developed an armament against insect herbivore attacks, and attackers continuously learn how to address it. Using a combined transcriptomic and metabolomic approach, we investigated the molecular and biochemical differences between Quercus robur L. trees that resisted (defined as resistant oak type) or were susceptible (defined as susceptible oak type) to infestation by the major oak pest, Tortrix viridana L. RESULTS: Next generation RNA sequencing revealed hundreds of genes that exhibited constitutive and/or inducible differential expression in the resistant oak compared to the susceptible oak. Distinct differences were found in the transcript levels and the metabolic content with regard to tannins, flavonoids, and terpenoids, which are compounds involved in the defence against insect pests. The results of our transcriptomic and metabolomic analyses are in agreement with those of a previous study in which we showed that female moths prefer susceptible oaks due to their specific profile of herbivore-induced volatiles. These data therefore define two oak genotypes that clearly differ on the transcriptomic and metabolomic levels, as reflected by their specific defensive compound profiles. CONCLUSIONS: We conclude that the resistant oak type seem to prefer a strategy of constitutive defence responses in contrast to more induced defence responses of the susceptible oaks triggered by feeding. These results pave the way for the development of biomarkers for an early determination of potentially green oak leaf roller-resistant genotypes in natural pedunculate oak populations in Europe.

A Small Set of Nuclear Markers for Reliable Differentiation of the Two Closely Related Oak Species Quercus Robur and Q. Petraea.

Quercus robur and Q. petraea are, in addition to Fagus sylvatica, the main economically used deciduous tree species in Europe. Identification of these two species is crucial because they differ in their ecological demands. Because of a changing climate, foresters must know more than ever which species will perform better under given environmental conditions. The search for differentiating molecular markers between these two species has already lasted for decades. Until now, differentiation has only been possible in approaches with a combination of several molecular markers and a subsequent statistical analysis to calculate the probability of being one or the other species. Here, we used MiSeq Illumina data from pools of Q. robur and Q. petraea specimens and identified nuclear SNPs and small InDels versus the Q. robur reference genome. Selected sequence variants with 100% allele frequency difference between the two pools were further validated in an extended set of Q. robur and Q. petraea specimens, and then the number of markers was deliberately reduced to the smallest possible set for species differentiation. A combination of six markers from four nuclear regions is enough to identify Q. robur, Q. petraea or hybrids between these two species quite well and represents a marker set that is cost-efficient and useable in every laboratory.

PlnTFDB: updated content and new features of the plant transcription factor database.

The Plant Transcription Factor Database (PlnTFDB; http://plntfdb.bio.uni-potsdam.de/v3.0/) is an integrative database that provides putatively complete sets of transcription factors (TFs) and other transcriptional regulators (TRs) in plant species (sensu lato) whose genomes have been completely sequenced and annotated. The complete sets of 84 families of TFs and TRs from 19 species ranging from unicellular red and green algae to angiosperms are included in PlnTFDB, representing >1.6 billion years of evolution of gene regulatory networks. For each gene family, a basic description is provided that is complemented by literature references, and multiple sequence alignments of protein domains. TF or TR gene entries include information of expressed sequence tags, 3D protein structures of homologous proteins, domain architecture and cross-links to other computational resources online. Moreover, the different species in PlnTFDB are linked to each other by means of orthologous genes facilitating cross-species comparisons.

PhosPhAt: the Arabidopsis thaliana phosphorylation site database. An update.

The PhosPhAt database of Arabidopsis phosphorylation sites was initially launched in August 2007. Since then, along with 10-fold increase in database entries, functionality of PhosPhAt (phosphat.mpimp-golm.mpg.de) has been considerably upgraded and re-designed. PhosPhAt is now more of a web application with the inclusion of advanced search functions allowing combinatorial searches by Boolean terms. The results output now includes interactive visualization of annotated fragmentation spectra and the ability to export spectra and peptide sequences as text files for use in other applications. We have also implemented dynamic links to other web resources thus augmenting PhosPhAt-specific information with external protein-related data. For experimental phosphorylation sites with information about dynamic behavior in response to external stimuli, we display simple time-resolved diagrams. We have included predictions for pT and pY sites and updated pS predictions. Access to prediction algorithm now allows 'on-the-fly' prediction of phosphorylation of any user-uploaded protein sequence. Protein Pfam domain structures are now mapped onto the protein sequence display next to experimental and predicted phosphorylation sites. Finally, we have implemented functional annotation of proteins using MAPMAN ontology. These new developments make the PhosPhAt resource a useful and powerful tool for the scientific community as a whole beyond the plant sciences.

ARR17 controls dioecy in Populus by repressing B-class MADS-box gene expression.

The number of dioecious species for which the genetic basis of sex determination has been resolved is rapidly increasing. Nevertheless, the molecular mechanisms downstream of the sex determinants remain largely elusive. Here, by RNA-sequencing early-flowering isogenic aspen (Populus tremula) lines differing exclusively for the sex switch gene ARR17, we show that a narrowly defined genetic network controls differential development of female and male flowers. Although ARR17 encodes a type-A response regulator supposedly involved in cytokinin (CK) hormone signalling, clustered regularly interspaced short palindromic repeats (CRISPR)-Cas9-mediated arr17 knockout only affected the expression of a strikingly small number of genes, indicating a specific role in the regulation of floral development rather than a generic function in hormone signalling. Notably, the UNUSUAL FLORAL ORGANS (UFO) gene, encoding an F-box protein acting as a transcriptional cofactor with LEAFY (LFY) to activate B-class MADS-box gene expression, and the B-class gene PISTILLATA (PI), necessary for male floral organ development, were strongly de-repressed in the arr17 CRISPR mutants. Our data highlight a CK-independent role of the poplar response regulator ARR17 and further emphasize the minimal differences between female and male individuals. This article is part of the theme issue 'Sex determination and sex chromosome evolution in land plants'.

GabiPD: the GABI primary database--a plant integrative 'omics' database.

The GABI Primary Database, GabiPD (http://www.gabipd.org/), was established in the frame of the German initiative for Genome Analysis of the Plant Biological System (GABI). The goal of GabiPD is to collect, integrate, analyze and visualize primary information from GABI projects. GabiPD constitutes a repository and analysis platform for a wide array of heterogeneous data from high-throughput experiments in several plant species. Data from different 'omics' fronts are incorporated (i.e. genomics, transcriptomics, proteomics and metabolomics), originating from 14 different model or crop species. We have developed the concept of GreenCards for text-based retrieval of all data types in GabiPD (e.g. clones, genes, mutant lines). All data types point to a central Gene GreenCard, where gene information is integrated from genome projects or NCBI UniGene sets. The centralized Gene GreenCard allows visualizing ESTs aligned to annotated transcripts as well as displaying identified protein domains and gene structure. Moreover, GabiPD makes available interactive genetic maps from potato and barley, and protein 2DE gels from Arabidopsis thaliana and Brassica napus. Gene expression and metabolic-profiling data can be visualized through MapManWeb. By the integration of complex data in a framework of existing knowledge, GabiPD provides new insights and allows for new interpretations of the data.

Proteome-wide survey of phosphorylation patterns affected by nuclear DNA polymorphisms in Arabidopsis thaliana.

BACKGROUND: Protein phosphorylation is an important post-translational modification influencing many aspects of dynamic cellular behavior. Site-specific phosphorylation of amino acid residues serine, threonine, and tyrosine can have profound effects on protein structure, activity, stability, and interaction with other biomolecules. Phosphorylation sites can be affected in diverse ways in members of any species, one such way is through single nucleotide polymorphisms (SNPs). The availability of large numbers of experimentally identified phosphorylation sites, and of natural variation datasets in Arabidopsis thaliana prompted us to analyze the effect of non-synonymous SNPs (nsSNPs) onto phosphorylation sites. RESULTS: From the analyses of 7,178 experimentally identified phosphorylation sites we found that: (i) Proteins with multiple phosphorylation sites occur more often than expected by chance. (ii) Phosphorylation hotspots show a preference to be located outside conserved domains. (iii) nsSNPs affected experimental phosphorylation sites as much as the corresponding non-phosphorylated amino acid residues. (iv) Losses of experimental phosphorylation sites by nsSNPs were identified in 86 A. thaliana proteins, among them receptor proteins were overrepresented.These results were confirmed by similar analyses of predicted phosphorylation sites in A. thaliana. In addition, predicted threonine phosphorylation sites showed a significant enrichment of nsSNPs towards asparagines and a significant depletion of the synonymous substitution. Proteins in which predicted phosphorylation sites were affected by nsSNPs (loss and gain), were determined to be mainly receptor proteins, stress response proteins and proteins involved in nucleotide and protein binding. Proteins involved in metabolism, catalytic activity and biosynthesis were less affected. CONCLUSIONS: We analyzed more than 7,100 experimentally identified phosphorylation sites in almost 4,300 protein-coding loci in silico, thus constituting the largest phosphoproteomics dataset for A. thaliana available to date. Our findings suggest a relatively high variability in the presence or absence of phosphorylation sites between different natural accessions in receptor and other proteins involved in signal transduction. Elucidating the effect of phosphorylation sites affected by nsSNPs on adaptive responses represents an exciting research goal for the future.

Single nucleotide polymorphisms in the allene oxide synthase 2 gene are associated with field resistance to late blight in populations of tetraploid potato cultivars.

The oomycete Phytophthora infestans causes late blight, the most relevant disease of potato (Solanum tuberosum) worldwide. Field resistance to late blight is a complex trait. When potatoes are cultivated under long day conditions in temperate climates, this resistance is correlated with late plant maturity, an undesirable characteristic. Identification of natural gene variation underlying late blight resistance not compromised by late maturity will facilitate the selection of resistant cultivars and give new insight in the mechanisms controlling quantitative pathogen resistance. We tested 24 candidate loci for association with field resistance to late blight and plant maturity in a population of 184 tetraploid potato individuals. The individuals were genotyped for 230 single nucleotide polymorphisms (SNPs) and 166 microsatellite alleles. For association analysis we used a mixed model, taking into account population structure, kinship, allele substitution and interaction effects of the marker alleles at a locus with four allele doses. Nine SNPs were associated with maturity corrected resistance (P < 0.001), which collectively explained 50% of the genetic variance of this trait. A major association was found at the StAOS2 locus encoding allene oxide synthase 2, a key enzyme in the biosynthesis of jasmonates, plant hormones that function in defense signaling. This finding supports StAOS2 as being one of the factors controlling natural variation of pathogen resistance.

The Diversity and Dynamics of Sex Determination in Dioecious Plants.

The diversity of inflorescences among flowering plants is captivating. Such charm is not only due to the variety of sizes, shapes, colors, and flowers displayed, but also to the range of reproductive systems. For instance, hermaphrodites occur abundantly throughout the plant kingdom with both stamens and carpels within the same flower. Nevertheless, 10% of flowering plants have separate unisexual flowers, either in different locations of the same individual (monoecy) or on different individuals (dioecy). Despite their rarity, dioecious plants provide an excellent opportunity to investigate the mechanisms involved in sex expression and the evolution of sex-determining regions (SDRs) and sex chromosomes. The SDRs and the evolution of dioecy have been studied in many species ranging from Ginkgo to important fruit crops. Some of these studies, for example in asparagus or kiwifruit, identified two sex-determining genes within the non-recombining SDR and may thus be consistent with the classical model for the evolution of dioecy from hermaphroditism via gynodioecy, that predicts two successive mutations, the first one affecting male and the second one female function, becoming linked in a region of suppressed recombination. On the other hand, aided by genome sequencing and gene editing, single factor sex determination has emerged in other species, such as persimmon or poplar. Despite the diversity of sex-determining mechanisms, a tentative comparative analysis of the known sex-determining genes and candidates in different species suggests that similar genes and pathways may be employed repeatedly for the evolution of dioecy. The cytokinin signaling pathway appears important for sex determination in several species regardless of the underlying genetic system. Additionally, tapetum-related genes often seem to act as male-promoting factors when sex is determined via two genes. We present a unified model that synthesizes the genetic networks of sex determination in monoecious and dioecious plants and will support the generation of hypothesis regarding candidate sex determinants in future studies.

Mitochondrial Genome of Fagus sylvatica L. as a Source for Taxonomic Marker Development in the Fagales.

European beech, Fagus sylvatica L., is one of the most important and widespread deciduous tree species in Central Europe and is widely managed for its hard wood. The complete DNA sequence of the mitochondrial genome of Fagus sylvatica L. was assembled and annotated based on Illumina MiSeq reads and validated using long reads from nanopore MinION sequencing. The genome assembled into a single DNA sequence of 504,715 bp in length containing 58 genes with predicted function, including 35 protein-coding, 20 tRNA and three rRNA genes. Additionally, 23 putative protein-coding genes were predicted supported by RNA-Seq data. Aiming at the development of taxon-specific mitochondrial genetic markers, the tool SNPtax was developed and applied to select genic SNPs potentially specific for different taxa within the Fagales. Further validation of a small SNP set resulted in the development of four CAPS markers specific for Fagus, Fagaceae, or Fagales, respectively, when considering over 100 individuals from a total of 69 species of deciduous trees and conifers from up to 15 families included in the marker validation. The CAPS marker set is suitable to identify the genus Fagus in DNA samples from tree tissues or wood products, including wood composite products.

Correction to: Long-term study of a subdioecious Populus ×canescens family reveals sex lability of females and reproduction behaviour of cosexual plants.

Table 4 in the original publication reports incomplete genotype names in the column "Cross" and wrong codes in the column "Generation".

Long-term study of a subdioecious Populus ×canescens family reveals sex lability of females and reproduction behaviour of cosexual plants.

KEY MESSAGE: Cosexual Populus ×canescens plants are inconstant females with life course plasticity of sex phenotype and can reproduce by selfing. Populus species are dioecious, but deviations from dioecy are reported in some cases. The objectives of this study were to investigate the phenotypic expression and the inheritance of subdioecy in a Populus ×canescens pedigree. The F1 progeny was monitored for sex during 14 years. Thirty per cent of individuals expressed deviations from dioecy and long-term plasticity of sex. Some plants started flowering as male, then became cosexual, and finally turned female. Two cosexual individuals were self-pollinated and generated a selfed progeny markedly impaired by inbreeding depression, but able to reproduce by outcrossing. Sex segregation of the F1 progeny statistically fitted the expected ratio 1:2:1 (female:male:cosexual). By analysis of DNA markers, the cosexual individuals were genetically clustered with the females. The segregation ratio and the genetic profile indicated that cosexual plants were female with altered sex phenotype. Linkage analysis identified a putative sex-determining region with suppressed recombination on chromosome 19 of the male Populus tremula parent. The male sex trait was linked to the pericentromeric region of the P. tremula chromosome 19, whereas the cosexual trait was linked to chromosome 19 of the female Populus alba parent. A genetic model is proposed to explain inheritance and phenotypic expression of sex.