Drought stress-induced Picea abies transcriptome changes in the context of functional interactions

Molecular responses to drought stress have been mainly studied in deciduous tree species although conifers dominate boreal forests. Here, we analysed the transcriptional response of Picea abies (L.) H. Karst. needles after exposure to severe drought by quantitative RNA-sequencing. In total, 2,402 differentially expressed genes (DEGs) were identified, of which 1,186 were up- and 1,216 downregulated. The upregulated DEGs are mainly involved in responses to stress, nitrogen compound, water deprivation, and abscisic acid as well as in channel activity. Although only one bZIP was identified among the DEGs, several other transcription factors involved in ABA-dependent pathways such as MYB, bHLH and WRKY showed differential expression. AP2/EREBP transcription factors related to ABA-independent pathways were also identified as DEGs. A functional interaction network of the 40 most connected Arabidopsis thaliana homologs of all Picea abies DEGs placed the two top-hubs P5CS1 and P5CS2 in the center. P5CS1 is the key enzyme in the biosynthesis of proline known to be accumulated in plants under abiotic stress. Lignin synthesis and DNA-related processes, among others, are overrepresented in this network. Our data highlight interesting gene targets for functional studies and natural genetic variation analyses to support the future identification and selection of potential drought tolerant trees

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

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 indifferent 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 otherhand, 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

Identification of full-sibling families from natural single-tree ash progenies based on SSR markers and genome-wide SNPs

Common ash, Fraxinus excelsior, is threatened by the invasive pathogen Hymenoscyphus fraxineus, which causes ash dieback. The pathogen is rapidly spreading throughout Europe with severe ecological and economic consequences. Multiple studies have presented evidence for the existence of a small fraction of genotypes with low susceptibility. Such genotypes can be targets for natural and artificial selection to conserve F. excelsior and associated ecosystems. To resolve the genetic architecture of variation in susceptibility it is necessary to analyze segregating populations. Here we employed about 1000 individuals of each of four single-tree progenies from potentially tolerant mother trees to identify full-sibling (full-sib) families. To this end, we first genotyped all 4000 individuals and the four mothers with eight SSR markers. We then used the program COLONY to predict full-sibs without knowledge of the paternal genotypes. For each single-tree progeny, COLONY predicted dozens of full-sib families, ranging from 3–166 individuals. In the next step, 910 individuals assigned to full-sib families with more than 28 individuals were subjected to high-resolution genotyping using over one million genome-wide SNPs which were identified with Illumina low-coverage resequencing. Using these SNP genotyping data in principal component analyses we were able to assign individuals to full-sib families with high confidence. Together the analyses revealed five large families with 73–212 individuals. These can be used to generate genetic linkage maps and to perform quantitative trait locus analyses for ash dieback susceptibility. The elucidation of the genetic basis of natural variation in ash may support breeding and conservation efforts and may contribute to more robust forest ecosystems

Haplotype inference from unphased SNP data in heterozygous polyploids based on SAT

<p>Abstract</p> <p>Background</p> <p>Haplotype inference based on unphased SNP markers is an important task in population genetics. Although there are different approaches to the inference of haplotypes in diploid species, the existing software is not suitable for inferring haplotypes from unphased SNP data in polyploid species, such as the cultivated potato (<it>Solanum tuberosum</it>). Potato species are tetraploid and highly heterozygous.</p> <p>Results</p> <p>Here we present the software SATlotyper which is able to handle polyploid and polyallelic data. SATlo-typer uses the Boolean satisfiability problem to formulate Haplotype Inference by Pure Parsimony. The software excludes existing haplotype inferences, thus allowing for calculation of alternative inferences. As it is not known which of the multiple haplotype inferences are best supported by the given unphased data set, we use a bootstrapping procedure that allows for scoring of alternative inferences. Finally, by means of the bootstrapping scores, it is possible to optimise the phased genotypes belonging to a given haplotype inference. The program is evaluated with simulated and experimental SNP data generated for heterozygous tetraploid populations of potato. We show that, instead of taking the first haplotype inference reported by the program, we can significantly improve the quality of the final result by applying additional methods that include scoring of the alternative haplotype inferences and genotype optimisation. For a sub-population of nineteen individuals, the predicted results computed by SATlotyper were directly compared with results obtained by experimental haplotype inference via sequencing of cloned amplicons. Prediction and experiment gave similar results regarding the inferred haplotypes and phased genotypes.</p> <p>Conclusion</p> <p>Our results suggest that Haplotype Inference by Pure Parsimony can be solved efficiently by the SAT approach, even for data sets of unphased SNP from heterozygous polyploids. SATlotyper is freeware and is distributed as a Java JAR file. The software can be downloaded from the webpage of the GABI Primary Database at <url>http://www.gabipd.org/projects/satlotyper/</url>. The application of SATlotyper will provide haplotype information, which can be used in haplotype association mapping studies of polyploid plants.</p

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

Coming from dry regions Norway spruce seedlings suffer less under drought

Norway spruce seedlings from dry lowland regions as well as from southern provenance are more drought tolerant than those from German highland provenances. • Early stress detection using biomarkers indicates the later mortality dynamics of seedlings under extreme drought well. • Higher drought tolerance is, however, not linked to higher resistance of adult trees of the same provenance against biotic damages

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.Data Availability Statement: DNA sequences from Sanger sequencing of the indel containing fragments derived from this study have been deposited in GenBank (Accession numbers: KU201020- KU201034). Next-generation sequences supporting this study are available from the NCBI GenBank as BioProject PRJNA269970 (SRA Accessions SRS954648, SRS954649, and SRS954650, SRX1050444-1050539). SRA accession numbers of the next-generation sequencing data of the two DNApools are SRR3624658 (Q. robur) and SRR3624680 (Q. mongolica). Draft chloroplast genome references will be submitted to Dryad soon. The open source software Variant Tools is available on  https://github. com/ThuenenFG/varianttools

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