Epigenetic Variability Among Saffron Crocus (Crocus sativus L.) Accessions Characterized by Different Phenotypes.

This work represents the first epigenomic study carried out on saffron crocus. Five accessions of saffron, showing differences in tepal pigmentation, yield of saffron and flowering time, were analyzed at the epigenetic level by applying a methylation-sensitive restriction enzyme-sequencing (MRE-seq) approach. Five accession-specific hypomethylomes plus a reference hypomethylome, generated by combining the sequence data from the single accessions, were obtained. Assembled sequences were annotated against existing online databases. In the absence of the Crocus genome, the rice genome was mainly used as the reference as it is the best annotated genome among monocot plants. Comparison of the hypomethylomes revealed many differentially methylated regions, confirming the high epigenetic variability present among saffron accessions, including sequences encoding for proteins that could be good candidates to explain the accessions' alternative phenotypes. In particular, transcription factors involved in flowering process (MADS-box and TFL) and for the production of pigments (MYB) were detected. Finally, by comparing the generated sequences of the different accessions, a high number of SNPs, likely having arisen as a consequence of the prolonged vegetative propagation, were detected, demonstrating surprisingly high genetic variability. Gene ontology (GO) was performed to map and visualize sequence polymorphisms located within the GOs and to compare their distributions among different accessions. As well as suggesting the possible existence of alternative phenotypes with a genetic basis, a clear difference in polymorphic GO is present among accessions based on their geographic origin, supporting a possible signature of selection in the Indian accession with respect to the Spanish ones.

Elucidating Drought Stress Tolerance in European Oaks Through Cross-Species Transcriptomics.

The impact of climate change that comes with a dramatic increase of long periods of extreme summer drought associated with heat is a fundamental challenge for European forests. As a result, forests are expected to shift their distribution patterns toward north-east, which may lead to a dramatic loss in value of European forest land. Consequently, unraveling key processes that underlie drought stress tolerance is not only of great scientific but also of utmost economic importance for forests to withstand future heat and drought wave scenarios. To reveal drought stress-related molecular patterns we applied cross-species comparative transcriptomics of three major European oak species: the less tolerant deciduous pedunculate oak (Quercus robur), the deciduous but quite tolerant pubescent oak (Q. pubescens), and the very tolerant evergreen holm oak (Q. ilex). We found 415, 79, and 222 differentially expressed genes during drought stress in Q. robur, Q. pubescens, and Q. ilex, respectively, indicating species-specific response mechanisms. Further, by comparative orthologous gene family analysis, 517 orthologous genes could be characterized that may play an important role in drought stress adaptation on the genus level. New regulatory candidate pathways and genes in the context of drought stress response were identified, highlighting the importance of the antioxidant capacity, the mitochondrial respiration machinery, the lignification of the water transport system, and the suppression of drought-induced senescence - providing a valuable knowledge base that could be integrated in breeding programs in the face of climate change.

YPR2 is a regulator of light modulated carbon and secondary metabolism in Trichoderma reesei.

BACKGROUND: Filamentous fungi have evolved to succeed in nature by efficient growth and degradation of substrates, but also due to the production of secondary metabolites including mycotoxins. For Trichoderma reesei, as a biotechnological workhorse for homologous and heterologous protein production, secondary metabolite secretion is of particular importance for industrial application. Recent studies revealed an interconnected regulation of enzyme gene expression and carbon metabolism with secondary metabolism. RESULTS: Here, we investigated gene regulation by YPR2, one out of two transcription factors located within the SOR cluster of T. reesei, which is involved in biosynthesis of sorbicillinoids. Transcriptome analysis showed that YPR2 exerts its major function in constant darkness upon growth on cellulose. Targets (direct and indirect) of YPR2 overlap with induction specific genes as well as with targets of the carbon catabolite repressor CRE1 and a considerable proportion is regulated by photoreceptors as well. Functional category analysis revealed both effects on carbon metabolism and secondary metabolism. Further, we found indications for an involvement of YPR2 in regulation of siderophores. In agreement with transcriptome data, mass spectrometric analyses revealed a broad alteration in metabolite patterns in ∆ypr2. Additionally, YPR2 positively influenced alamethicin levels along with transcript levels of the alamethicin synthase tex1 and is essential for production of orsellinic acid in darkness. CONCLUSIONS: YPR2 is an important regulator balancing secondary metabolism with carbon metabolism in darkness and depending on the carbon source. The function of YPR2 reaches beyond the SOR cluster in which ypr2 is located and happens downstream of carbon catabolite repression mediated by CRE1.

Identification and characterization of the Fasciola hepatica sodium- and chloride-dependent taurine transporter.

The parasitic liver fluke Fasciola hepatica infests mainly ruminants, but it can also cause fasciolosis in people, who ingest the metacercariae encysted on plants. The drug of choice to treat fasciolosis is triclabendazole (TBZ), which has been on the market for several decades. This is also true for the other available drugs. Accordingly, drug-resistant flukes have been emerging at an increasing rate making it desirable to identify alternative drug targets. Here, we focused on the fact that adult F. hepatica persists in the hostile environment of the bile ducts of infected organisms. A common way to render bile acids less toxic is to conjugate them to taurine (2-aminoethanesulfonic acid). We cloned a transporter from the solute carrier-6 (SLC6) family, which was most closely related to the GABA-transporter-2 of other organisms. When heterologously expressed, this F. hepatica transporter supported the high-affinity cellular uptake of taurine (KM = 12.0 ± 0.5 µM) but not of GABA. Substrate uptake was dependent on Na+- and Cl- (calculated stoichiometry 2:1). Consistent with the low chloride concentration in mammalian bile, the F. hepatica transporter had a higher apparent affinity for Cl- (EC50 = 14±3 mM) than the human taurine transporter (EC50 = 55±7 mM). We incubated flukes with unconjugated bile acids in the presence and absence of taurine: taurine promoted survival of flukes; the taurine transporter inhibitor guanidinoethansulfonic acid abolished this protective effect of taurine. Based on these observations, we conclude that the taurine transporter is critical for the survival of liver flukes in the bile. Thus, the taurine transporter represents a candidate drug target.

Drought Sensitivity of Norway Spruce at the Species' Warmest Fringe: Quantitative and Molecular Analysis Reveals High Genetic Variation Among and Within Provenances.

Norway spruce (Picea abies) is by far the most important timber species in Europe, but its outstanding role in future forests is jeopardized by its high sensitivity to drought. We analyzed drought response of Norway spruce at the warmest fringe of its natural range. Based on a 35-year old provenance experiment we tested for genetic variation among and within seed provenances across consecutively occurring strong drought events using dendroclimatic time series. Moreover, we tested for associations between ≈1,700 variable SNPs and traits related to drought response, wood characteristics and climate-growth relationships. We found significant adaptive genetic variation among provenances originating from the species' Alpine, Central and Southeastern European range. Genetic variation between individuals varied significantly among provenances explaining up to 44% of the phenotypic variation in drought response. Varying phenotypic correlations between drought response and wood traits confirmed differences in selection intensity among seed provenances. Significant associations were found between 29 SNPs and traits related to drought, climate-growth relationships and wood properties which explained between 11 and 43% of trait variation, though 12 of them were due to single individuals having extreme phenotypes of the respective trait. The majority of these SNPs are located within exons of genes and the most important ones are preferentially expressed in cambium and xylem expansion layers. Phenotype-genotype associations were stronger if only provenances with significant quantitative genetic variation in drought response were considered. The present study confirms the high adaptive variation of Norway spruce in Central and Southeastern Europe and demonstrates how quantitative genetic, dendroclimatic and genomic data can be linked to understand the genetic basis of adaptation to climate extremes in trees.

Association genetics of phenolic needle compounds in Norway spruce with variable susceptibility to needle bladder rust.

KEY MESSAGE: Accumulation of phenolic needle metabolites in Norway spruce is regulated by many genes with small and additive effects and is correlated with the susceptibility against fungal attack. Norway spruce accumulates high foliar concentrations of secondary phenolic metabolites, with important functions for pathogen defence responses. However, the molecular genetic basis underlying the quantitative variation of phenolic compounds and their role in enhanced resistance of spruce to infection by needle bladder rust are unknown. To address these questions, a set of 1035 genome-wide single nucleotide polymorphisms (SNPs) was associated to the quantitative variation of four simple phenylpropanoids, eight stilbenes, nine flavonoids, six related arithmetic parameters and the susceptibility to infection by Chrysomyxa rhododendri in an unstructured natural population of Norway spruce. Thirty-one significant genetic associations for the flavonoids gallocatechin, kaempferol 3-glucoside and quercetin 3-glucoside and the stilbenes resveratrol, piceatannol, astringin and isorhapontin were discovered, explaining 22-59% of phenotypic variation, and indicating a regulation of phenolic accumulation by many genes with small and additive effects. The phenolics profile differed between trees with high and low susceptibility to the fungus, underlining the importance of phenolic compounds in the defence mechanisms of Norway spruce to C. rhododendri. Results highlight the utility of association studies in non-model tree species and may enable marker-assisted selection of Norway spruce adapted to severe pathogen attack.

Selecting Hypomethylated Genomic Regions Using MRE-Seq.

Here, we describe a method capable of filtering the hypomethylated part of plant genomes, the so-called hypomethylome. The principle of the method is based on the filtration and sequence analysis of small DNA fragments generated by methylation-sensitive four-cutter restriction endonucleases, possessing ((5me))CpG motifs in their recognition sites. The majority of these fragments represent genes and their flanking regions containing also regulatory elements-the gene space of the genome. Besides the enrichment of the gene space, another advantage of the method is the simultaneous depletion of repetitive elements due to their methylated nature and its easy application on complex and large plant genomes. Additionally to the wet lab procedure, we describe how to analyze the data using bioinformatics methods and how to apply the method to comparative studies.

How to Isolate a Plant's Hypomethylome in One Shot.

Genome assembly remains a challenge for large and/or complex plant genomes due to their abundant repetitive regions resulting in studies focusing on gene space instead of the whole genome. Thus, DNA enrichment strategies facilitate the assembly by increasing the coverage and simultaneously reducing the complexity of the whole genome. In this paper we provide an easy, fast, and cost-effective variant of MRE-seq to obtain a plant's hypomethylome by an optimized methyl filtration protocol followed by next generation sequencing. The method is demonstrated on three plant species with knowingly large and/or complex (polyploid) genomes: Oryza sativa, Picea abies, and Crocus sativus. The identified hypomethylomes show clear enrichment for genes and their flanking regions and clear reduction of transposable elements. Additionally, genomic sequences around genes are captured including regulatory elements in introns and up- and downstream flanks. High similarity of the results obtained by a de novo assembly approach with a reference based mapping in rice supports the applicability for studying and understanding the genomes of nonmodel organisms. Hence we show the high potential of MRE-seq in a wide range of scenarios for the direct analysis of methylation differences, for example, between ecotypes, individuals, within or across species harbouring large, and complex genomes.

Dissecting plant genomes with the PLAZA comparative genomics platform.

With the arrival of low-cost, next-generation sequencing, a multitude of new plant genomes are being publicly released, providing unseen opportunities and challenges for comparative genomics studies. Here, we present PLAZA 2.5, a user-friendly online research environment to explore genomic information from different plants. This new release features updates to previous genome annotations and a substantial number of newly available plant genomes as well as various new interactive tools and visualizations. Currently, PLAZA hosts 25 organisms covering a broad taxonomic range, including 13 eudicots, five monocots, one lycopod, one moss, and five algae. The available data consist of structural and functional gene annotations, homologous gene families, multiple sequence alignments, phylogenetic trees, and colinear regions within and between species. A new Integrative Orthology Viewer, combining information from different orthology prediction methodologies, was developed to efficiently investigate complex orthology relationships. Cross-species expression analysis revealed that the integration of complementary data types extended the scope of complex orthology relationships, especially between more distantly related species. Finally, based on phylogenetic profiling, we propose a set of core gene families within the green plant lineage that will be instrumental to assess the gene space of draft or newly sequenced plant genomes during the assembly or annotation phase.

Functional annotation of expressed sequence tags as a tool to understand the molecular mechanism controlling flower bud development in tree peony.

Tree peony (Paeonia suffruticosaAndrews) is an important medicinal and ornamental plant. In China, its root bark is an important ingredient for traditional Chinese medicine. It is valued as an ornamental plant because its flower shows a wide variation in shape and color. We used flower buds at different developmental stages to construct the first cDNA library for this organism. A total of 2241 raw expressed sequence tags (ESTs) were obtained after unidirectional sequencing. After processing and assembly, they resulted in a total of 1300 unigenes [363 contigs with an average size of 3.5 ESTs (ranging from 2 up to 36) and 937 singletons]. Gene Ontology categories were assigned and further summarized into 13 broad families with biological roles according to similar functional characteristics or cellular roles. A total of 185 single nucleotide polymorphisms were detected for all contigs. We were able to detect open reading frames in the consensus sequences of 1268 unigenes. 97.5% of the ESTs showed significant similarity to sequences present in public databases. One hundred and sixty-seven short sequence repeats were obtained in the whole data set. All the analysis and information will be valuable resources for a better understanding of this important plant and also can be used for functional study in Paeoniaceae.

Transcriptional profiling of mature Arabidopsis trichomes reveals that NOECK encodes the MIXTA-like transcriptional regulator MYB106.

Leaf hairs (trichomes) of Arabidopsis (Arabidopsis thaliana) have been extensively used as a model to address general questions in cell and developmental biology. Here, we lay the foundation for a systems-level understanding of the biology of this model cell type by performing genome-wide gene expression analyses. We have identified 3,231 genes that are up-regulated in mature trichomes relative to leaves without trichomes, and we compared wild-type trichomes with two mutants, glabra3 and triptychon, that affect trichome morphology and physiology in contrasting ways. We found that cell wall-related transcripts were particularly overrepresented in trichomes, consistent with their highly elaborated structure. In addition, trichome expression maps revealed high activities of anthocyanin, flavonoid, and glucosinolate pathways, indicative of the roles of trichomes in the biosynthesis of secondary compounds and defense. Interspecies comparisons revealed that Arabidopsis trichomes share many expressed genes with cotton (Gossypium hirsutum) fibers, making them an attractive model to study industrially important fibers. In addition to identifying physiological processes involved in the development of a specific cell type, we also demonstrated the utility of transcript profiling for identifying and analyzing regulatory gene function. One of the genes that are differentially expressed in fibers is the MYB transcription factor GhMYB25. A combination of transcript profiling and map-based cloning revealed that the NOECK gene of Arabidopsis encodes AtMYB106, a MIXTA-like transcription factor and homolog of cotton GhMYB25. However, in contrast to Antirrhinum, in which MIXTA promotes epidermal cell outgrowth, AtMYB106 appears to function as a repressor of cell outgrowth in Arabidopsis.

A role for AtWRKY23 in feeding site establishment of plant-parasitic nematodes.

During the interaction between sedentary plant-parasitic nematodes and their host, complex morphological and physiological changes occur in the infected plant tissue, finally resulting in the establishment of a nematode feeding site. This cellular transformation is the result of altered plant gene expression most likely induced by proteins injected in the plant cell by the nematode. Here, we report on the identification of a WRKY transcription factor expressed during nematode infection. Using both promoter-reporter gene fusions and in situ reverse transcription-polymerase chain reaction, we could show that AtWRKY23 is expressed during the early stages of feeding site establishment. Knocking down the expression of WRKY23 resulted in lower infection of the cyst nematode Heterodera schachtii. WRKY23 is an auxin-inducible gene and in uninfected plants WRKY23 acts downstream of the Aux/IAA protein SLR/IAA14. Although auxin is known to be involved in feeding site formation, our results suggest that, during early stages, auxin-independent signals might be at play to activate the initial expression of WRKY23.