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.

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.

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.

DNA of centuries-old timber can reveal its origin.

Oak wood was highly appreciated and widely used for construction in past centuries. As population sizes expanded in some regions of Europe, local forests were depleted of high-quality timber. Therefore, regions of soaring economies were importing timber initially from the European market and eventually from other continents. Origin of archaeological or historical timber is usually identified by means of dendroprovenancing, i.e. statistical matching of tree-ring-width (TRW) series of timber of unknown origin with TRW reference datasets. However, this method has pitfalls and limitations and therefore alternative techniques are needed. Here, we used three different DNA analysis methods to investigate the potential of using ancient (a)DNA, extracted from oak timber derived from historical buildings and shipwrecks from a variety of countries. All the material had also been analysed dendrochronologically, so its dating and provenance is demonstrated. We included heartwood samples in this analysis, for which DNA extraction is especially challenging as it contains chemicals that inhibit DNA amplification. We succeeded in amplifying DNA for at least one marker from 56% of samples (including heartwood samples), yielding crucial information that allowed us to identify the potential source area of centuries old timber buildings in Latvia and Denmark and of 750-year-old shipwreck material from Germany. Our results prove the strong potential of DNA analyses for identifying timber origin to the regional scale, but by combining these with the dendrochronological results, we can control the exactitude of the aDNA approach and demonstrate a more nuanced examination of the timber sources for these historic structures.

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.

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.

High Level of Conservation of Mitochondrial RNA Editing Sites Among Four Populus Species.

RNA editing occurs in the endosymbiont organelles of higher plants as C-to-U conversions of defined nucleotides. The availability of large quantities of RNA sequencing data makes it possible to identify RNA editing sites and to quantify their editing extent. We have investigated RNA editing in 34 protein-coding mitochondrial transcripts of four Populus species, a genus noteworthy for its remarkably small number of RNA editing sites compared to other angiosperms. 27 of these transcripts were subject to RNA editing in at least one species. In total, 355 RNA editing sites were identified with high confidence, their editing extents ranging from 10 to 100%. The most heavily edited transcripts were ccmB with the highest density of RNA editing sites (53.7 sites / kb) and ccmFn with the highest number of sites (39 sites). Most of the editing events are at position 1 or 2 of the codons, usually altering the encoded amino acid, and are highly conserved among the species, also with regard to their editing extent. However, one SNP was found in the newly sequenced and annotated mitochondrial genome of P. alba resulting in the loss of an RNA editing site compared to P. tremula and P. davidiana This SNP causes a C-to-T transition and an amino acid exchange from Ser to Phe, highlighting the widely discussed role of RNA editing in compensating mutations.

The complete chloroplast genome sequence of Pinus cembra L. (Pinaceae).

The Swiss pine (Pinus cembra) is a montane tree in Central Europe and, therefore, known for its hardiness against severe winter colds. The seeds are harvested and eaten as pine nuts. We assembled and characterized the complete chloroplast genome of P. cembra to serve as a valuable resource in future genetic studies. The complete plastome sequence is 116,609 bp in length and contains 113 genes including 79 protein-coding genes, 30 tRNA genes, and 4 rRNA genes. A phylogenetic analysis of 34 Pinus plastome sequences shows that Pinus sibirica is the nearest relative to P. cembra and that there is a distinct clustering together with the other members of the section Quinquefoliae.

Development of Molecular Markers for Determining Continental Origin of Wood from White Oaks (Quercus L. sect. Quercus).

To detect and avoid illegal logging of valuable tree species, identification methods for the origin of timber are necessary. We used next-generation sequencing to identify chloroplast genome regions that differentiate the origin of white oaks from the three continents; Asia, Europe, and North America. By using the chloroplast genome of Asian Q. mongolica as a reference, we identified 861 variant sites (672 single nucleotide polymorphisms (SNPs); 189 insertion/deletion (indel) polymorphism) from representative species of three continents (Q. mongolica from Asia; Q. petraea and Q. robur from Europe; Q. alba from North America), and we identified additional chloroplast polymorphisms in pools of 20 individuals each from Q. mongolica (789 variant sites) and Q. robur (346 variant sites). Genome sequences were screened for indels to develop markers that identify continental origin of oak species, and that can be easily evaluated using a variety of detection methods. We identified five indels and one SNP that reliably identify continent-of-origin, based on evaluations of up to 1078 individuals representing 13 white oak species and three continents. Due to the size of length polymorphisms revealed, this marker set can be visualized using capillary electrophoresis or high resolution gel (acrylamide or agarose) electrophoresis. With these markers, we provide the wood trading market with an instrument to comply with the U.S. and European laws that require timber companies to avoid the trade of illegally harvested timber.

Genome Sequences of Populus tremula Chloroplast and Mitochondrion: Implications for Holistic Poplar Breeding.

Complete Populus genome sequences are available for the nucleus (P. trichocarpa; section Tacamahaca) and for chloroplasts (seven species), but not for mitochondria. Here, we provide the complete genome sequences of the chloroplast and the mitochondrion for the clones P. tremula W52 and P. tremula x P. alba 717-1B4 (section Populus). The organization of the chloroplast genomes of both Populus clones is described. A phylogenetic tree constructed from all available complete chloroplast DNA sequences of Populus was not congruent with the assignment of the related species to different Populus sections. In total, 3,024 variable nucleotide positions were identified among all compared Populus chloroplast DNA sequences. The 5-prime part of the LSC from trnH to atpA showed the highest frequency of variations. The variable positions included 163 positions with SNPs allowing for differentiating the two clones with P. tremula chloroplast genomes (W52, 717-1B4) from the other seven Populus individuals. These potential P. tremula-specific SNPs were displayed as a whole-plastome barcode on the P. tremula W52 chloroplast DNA sequence. Three of these SNPs and one InDel in the trnH-psbA linker were successfully validated by Sanger sequencing in an extended set of Populus individuals. The complete mitochondrial genome sequence of P. tremula is the first in the family of Salicaceae. The mitochondrial genomes of the two clones are 783,442 bp (W52) and 783,513 bp (717-1B4) in size, structurally very similar and organized as single circles. DNA sequence regions with high similarity to the W52 chloroplast sequence account for about 2% of the W52 mitochondrial genome. The mean SNP frequency was found to be nearly six fold higher in the chloroplast than in the mitochondrial genome when comparing 717-1B4 with W52. The availability of the genomic information of all three DNA-containing cell organelles will allow a holistic approach in poplar molecular breeding in the future.

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.

Function of defensive volatiles in pedunculate oak (Quercus robur) is tricked by the moth Tortrix viridana.

The indirect defences of plants are comprised of herbivore-induced plant volatiles (HIPVs) that among other things attract the natural enemies of insects. However, the actual extent of the benefits of HIPV emissions in complex co-evolved plant-herbivore systems is only poorly understood. The observation that a few Quercus robur L. trees constantly tolerated (T-oaks) infestation by a major pest of oaks (Tortrix viridana L.), compared with heavily defoliated trees (susceptible: S-oaks), lead us to a combined biochemical and behavioural study. We used these evidently different phenotypes to analyse whether the resistance of T-oaks to the herbivore was dependent on the amount and scent of HIPVs and/or differences in non-volatile polyphenolic leaf constituents (as quercetin-, kaempferol- and flavonol glycosides). In addition to non-volatile metabolic differences, typically defensive HIPV emissions differed between S-oaks and T-oaks. Female moths were attracted by the blend of HIPVs from S-oaks, showing significantly higher amounts of (E)-4,8-dimethyl-1,3,7-nonatriene (DMNT) and (E)-ß-ocimene and avoid T-oaks with relative high fraction of the sesquiterpenes α-farnesene and germacrene D. Hence, the strategy of T-oaks exhibiting directly herbivore-repellent HIPV emissions instead of high emissions of predator-attracting HIPVs of the S-oaks appears to be the better mechanism for avoiding defoliation.

A very small and isolated population of the green Oak Leaf Roller, Tortrix viridana L., with high genetic diversity--how does this work?

A geographically isolated population of an herbivorous insect (Tortrix viridana, Lepidoptera: Tortricidae) in the Bashkir Transural region and 5 further populations were investigated for genetic variation using eight microsatellite markers. The sample size per population was between 48 and 62 individuals. The genetic variation was higher within the isolated population than within populations in the center of the distribution area. No bottleneck effects were discovered during analyses that could have formatted the gene pool of this population. Balancing or directed selection toward preservation of specific alleles or promotion of heterozygous individuals could be an explanation for the unexpected high genetic diversity within this small and isolated population.

Insects' colonization of human corpses in warm and cold season.

The succession of insects is commonly used for estimation of the postmortem interval (PMI). In the warm season there are many insect species, which enable us to estimate the time of death, but in the cold season there are less species. For the estimation of a seasonal pattern, we examined 117 human corpses found in domestic environment during the year 2001. We found three typical 'summer species', active from May until October, three further species active from April until November. One blowfly species (Calliphora vicina) and the family of scuttle flies were active all over the year. These insects may be useful for PMI estimation in winter; however, further research concerning the developmental rate of scuttle flies is required. The distinct seasonal pattern of the early colonizing blowfly species allows to allocate the time of death to a particular time interval even in cases of long PMIs if empty puparia of these insects are found.