A transcriptomic time-series reveals differing trajectories during pre-floral development in the apex and leaf in winter and spring varieties of Brassica napus.

Oilseed rape (Brassica napus) is an important global oil crop, with spring and winter varieties grown commercially. To understand the transcriptomic differences between these varieties, we collected transcriptomes from apex and leaf tissue from a spring variety, Westar, and a winter variety, Tapidor, before, during, and after vernalisation treatment, until the plants flowered. Large transcriptomic differences were noted in both varieties during the vernalisation treatment because of temperature and day length changes. Transcriptomic alignment revealed that the apex transcriptome reflects developmental state, whereas the leaf transcriptome is more closely aligned to the age of the plant. Similar numbers of copies of genes were expressed in both varieties during the time series, although key flowering time genes exhibited expression pattern differences. BnaFLC copies on A2 and A10 are the best candidates for the increased vernalisation requirement of Tapidor. Other BnaFLC copies show tissue-dependent reactivation of expression post-cold, with these dynamics suggesting some copies have retained or acquired a perennial nature. BnaSOC1 genes, also related to the vernalisation pathway, have expression profiles which suggest tissue subfunctionalisation. This understanding may help to breed varieties with more consistent or robust vernalisation responses, of special importance due to the milder winters resulting from climate change.

Characterisation of the Introgression of Brassica villosa Genome Into Broccoli to Enhance Methionine-Derived Glucosinolates and Associated Health Benefits.

Broccoli cultivars that have enhanced accumulation of methionine-derived glucosinolates have been developed through the introgression of a novel allele of the MYB28 transcription factor from the wild species Brassica villosa. Through a novel k-mer approach, we characterised the extent of the introgression of unique B. villosa genome sequences into high glucosinolate broccoli genotypes. RNAseq analyses indicated that the introgression of the B. villosa MYB28 C2 allele resulted in the enhanced expression of the MYB28 transcription factor, and modified expression of genes associated with sulphate absorption and reduction, and methionine and glucosinolate biosynthesis when compared to standard broccoli. A adenine-thymine (AT) short tandem repeat (STR) was identified within the 5' untranslated region (UTR) B. villosa MYB28 allele that was absent from two divergent cultivated forms of Brassica oleracea, which may underpin the enhanced expression of B. villosa MYB28.

Natural Variation in OASC Gene for Mitochondrial O-Acetylserine Thiollyase Affects Sulfate Levels in Arabidopsis.

Sulfur plays a vital role in the primary and secondary metabolism of plants, and carries an important function in a large number of different compounds. Despite this importance, compared to other mineral nutrients, relatively little is known about sulfur sensing and signalling, as well as about the mechanisms controlling sulfur metabolism and homeostasis. Sulfur contents in plants vary largely not only among different species, but also among accessions of the same species. We previously used associative transcriptomics to identify several genes potentially controlling variation in sulfate content in the leaves of Brassica napus, including an OASC gene for mitochondrial O-acetylserine thiollyase (OAS-TL), an enzyme involved in cysteine synthesis. Here, we show that loss of OASC in Arabidopsis thaliana lowers not only sulfate, but also glutathione levels in the leaves. The reduced accumulation is caused by lower sulfate uptake and translocation to the shoots; however, the flux through the pathway is not affected. In addition, we identified a single nucleotide polymorphism in the OASC gene among A. thaliana accessions that is linked to variation in sulfate content. Both genetic and transgenic complementation confirmed that the exchange of arginine at position 81 for lysine in numerous accessions resulted in a less active OASC and a lower sulfate content in the leaves. The mitochondrial isoform of OAS-TL is, thus, after the ATPS1 isoform of sulfurylase and the APR2 form of APS reductase 2, the next metabolic enzyme with a role in regulation of sulfate content in Arabidopsis.

Pan-genome analysis identifies intersecting roles for Pseudomonas specialized metabolites in potato pathogen inhibition.

Agricultural soil harbors a diverse microbiome that can form beneficial relationships with plants, including the inhibition of plant pathogens. Pseudomonas spp. are one of the most abundant bacterial genera in the soil and rhizosphere and play important roles in promoting plant health. However, the genetic determinants of this beneficial activity are only partially understood. Here, we genetically and phenotypically characterize the Pseudomonas fluorescens population in a commercial potato field, where we identify strong correlations between specialized metabolite biosynthesis and antagonism of the potato pathogens Streptomyces scabies and Phytophthora infestans. Genetic and chemical analyses identified hydrogen cyanide and cyclic lipopeptides as key specialized metabolites associated with S. scabies inhibition, which was supported by in planta biocontrol experiments. We show that a single potato field contains a hugely diverse and dynamic population of Pseudomonas bacteria, whose capacity to produce specialized metabolites is shaped both by plant colonization and defined environmental inputs.

Potato scab and blight are two major diseases which can cause heavy crop losses. They are caused, respectively, by the bacterium Streptomyces scabies and an oomycete (a fungus-like organism) known as Phytophthora infestans. Fighting these disease-causing microorganisms can involve crop management techniques ­ for example, ensuring that a field is well irrigated helps to keep S. scabies at bay. Harnessing biological control agents can also offer ways to control disease while respecting the environment. Biocontrol bacteria, such as Pseudomonas, can produce compounds that keep S. scabies and P. infestans in check. However, the identity of these molecules and how irrigation can influence Pseudomonas population remains unknown. To examine these questions, Pacheco-Moreno et al. sampled and isolated hundreds of Pseudomonas strains from a commercial potato field, closely examining the genomes of 69 of these. Comparing the genetic information of strains based on whether they could control the growth of S. scabies revealed that compounds known as cyclic lipopeptides are key to controlling the growth of S. scabies and P. infestans. Whether the field was irrigated also had a large impact on the strains forming the Pseudomonas population. Working out how Pseudomonas bacteria block disease could speed up the search for biological control agents. The approach developed by Pacheco-Moreno et al. could help to predict which strains might be most effective based on their genetic features. Similar experiments could also work for other combinations of plants and diseases.

Restriction of cytosolic sucrose hydrolysis profoundly alters development, metabolism, and gene expression in Arabidopsis roots.

Plant roots depend on sucrose imported from leaves as the substrate for metabolism and growth. Sucrose and hexoses derived from it are also signalling molecules that modulate growth and development, but the importance for signalling of endogenous changes in sugar levels is poorly understood. We report that reduced activity of cytosolic invertase, which converts sucrose to hexoses, leads to pronounced metabolic, growth, and developmental defects in roots of Arabidopsis (Arabidopsis thaliana) seedlings. In addition to altered sugar and downstream metabolite levels, roots of cinv1 cinv2 mutants have reduced elongation rates, cell and meristem size, abnormal meristematic cell division patterns, and altered expression of thousands of genes of diverse functions. Provision of exogenous glucose to mutant roots repairs relatively few of the defects. The extensive transcriptional differences between mutant and wild-type roots have hallmarks of both high sucrose and low hexose signalling. We conclude that the mutant phenotype reflects both low carbon availability for metabolism and growth and complex sugar signals derived from elevated sucrose and depressed hexose levels in the cytosol of mutant roots. Such reciprocal changes in endogenous sucrose and hexose levels potentially provide rich information about sugar status that translates into flexible adjustments of growth and development.

QTL-seq identifies BnaFT.A02 and BnaFLC.A02 as candidates for variation in vernalization requirement and response in winter oilseed rape (Brassica napus).

Winter, spring and biennial varieties of Brassica napus that vary in vernalization requirement are grown for vegetable and oil production. Here, we show that the obligate or facultative nature of the vernalization requirement in European winter oilseed rape is determined by allelic variation at a 10 Mbp region on chromosome A02. This region includes orthologues of the key floral regulators FLOWERING LOCUS C (BnaFLC.A02) and FLOWERING LOCUS T (BnaFT.A02). Polymorphism at BnaFLC.A02 and BnaFT.A02, mostly in cis-regulatory regions, results in distinct gene expression dynamics in response to vernalization treatment. Our data suggest allelic variation at BnaFT.A02 is associated with flowering time in the absence of vernalization, while variation at BnaFLC.A02 is associated with flowering time under vernalizing conditions. We hypothesize selection for BnaFLC.A02 and BnaFT.A02 gene expression variation has facilitated the generation of European winter oilseed rape varieties that are adapted to different winter climates. This knowledge will allow for the selection of alleles of flowering time regulators that alter the vernalization requirement of oilseed rape, informing the generation of new varieties with adapted flowering times and improved yields.

A carbohydrate-binding protein, B-GRANULE CONTENT 1, influences starch granule size distribution in a dose-dependent manner in polyploid wheat.

In Triticeae endosperm (e.g. wheat and barley), starch granules have a bimodal size distribution (with A- and B-type granules) whereas in other grasses the endosperm contains starch granules with a unimodal size distribution. Here, we identify the gene, BGC1 (B-GRANULE CONTENT 1), responsible for B-type starch granule content in Aegilops and wheat. Orthologues of this gene are known to influence starch synthesis in diploids such as rice, Arabidopsis, and barley. However, using polyploid Triticeae species, we uncovered a more complex biological role for BGC1 in starch granule initiation: BGC1 represses the initiation of A-granules in early grain development but promotes the initiation of B-granules in mid grain development. We provide evidence that the influence of BGC1 on starch synthesis is dose dependent and show that three very different starch phenotypes are conditioned by the gene dose of BGC1 in polyploid wheat: normal bimodal starch granule morphology; A-granules with few or no B-granules; or polymorphous starch with few normal A- or B-granules. We conclude from this work that BGC1 participates in controlling B-type starch granule initiation in Triticeae endosperm and that its precise effect on granule size and number varies with gene dose and stage of development.

RNA-seq, de novo transcriptome assembly and flavonoid gene analysis in 13 wild and cultivated berry fruit species with high content of phenolics.

BACKGROUND: Flavonoids are produced in all flowering plants in a wide range of tissues including in berry fruits. These compounds are of considerable interest for their biological activities, health benefits and potential pharmacological applications. However, transcriptomic and genomic resources for wild and cultivated berry fruit species are often limited, despite their value in underpinning the in-depth study of metabolic pathways, fruit ripening as well as in the identification of genotypes rich in bioactive compounds. RESULTS: To access the genetic diversity of wild and cultivated berry fruit species that accumulate high levels of phenolic compounds in their fleshy berry(-like) fruits, we selected 13 species from Europe, South America and Asia representing eight genera, seven families and seven orders within three clades of the kingdom Plantae. RNA from either ripe fruits (ten species) or three ripening stages (two species) as well as leaf RNA (one species) were used to construct, assemble and analyse de novo transcriptomes. The transcriptome sequences are deposited in the BacHBerryGEN database (http://jicbio.nbi.ac.uk/berries) and were used, as a proof of concept, via its BLAST portal (http://jicbio.nbi.ac.uk/berries/blast.html) to identify candidate genes involved in the biosynthesis of phenylpropanoid compounds. Genes encoding regulatory proteins of the anthocyanin biosynthetic pathway (MYB and basic helix-loop-helix (bHLH) transcription factors and WD40 repeat proteins) were isolated using the transcriptomic resources of wild blackberry (Rubus genevieri) and cultivated red raspberry (Rubus idaeus cv. Prestige) and were shown to activate anthocyanin synthesis in Nicotiana benthamiana. Expression patterns of candidate flavonoid gene transcripts were also studied across three fruit developmental stages via the BacHBerryEXP gene expression browser (http://www.bachberryexp.com) in R. genevieri and R. idaeus cv. Prestige. CONCLUSIONS: We report a transcriptome resource that includes data for a wide range of berry(-like) fruit species that has been developed for gene identification and functional analysis to assist in berry fruit improvement. These resources will enable investigations of metabolic processes in berries beyond the phenylpropanoid biosynthetic pathway analysed in this study. The RNA-seq data will be useful for studies of berry fruit development and to select wild plant species useful for plant breeding purposes.

Variation in Expression of the HECT E3 Ligase UPL3 Modulates LEC2 Levels, Seed Size, and Crop Yields in Brassica napus.

Identifying genetic variation that increases crop yields is a primary objective in plant breeding. We used association analyses of oilseed rape/canola (Brassica napus) accessions to identify genetic variation that influences seed size, lipid content, and final crop yield. Variation in the promoter region of the HECT E3 ligase gene BnaUPL3 C03 made a major contribution to variation in seed weight per pod, with accessions exhibiting high seed weight per pod having lower levels of BnaUPL3 C03 expression. We defined a mechanism in which UPL3 mediated the proteasomal degradation of LEC2, a master transcriptional regulator of seed maturation. Accessions with reduced UPL3 expression had increased LEC2 protein levels, larger seeds, and prolonged expression of lipid biosynthetic genes during seed maturation. Natural variation in BnaUPL3 C03 expression appears not to have been exploited in current B napus breeding lines and could therefore be used as a new approach to maximize future yields in this important oil crop.

NIN Acts as a Network Hub Controlling a Growth Module Required for Rhizobial Infection.

The symbiotic infection of root cells by nitrogen-fixing rhizobia during nodulation requires the transcription factor Nodule Inception (NIN). Our root hair transcriptomic study extends NIN's regulon to include Rhizobium Polar Growth and genes involved in cell wall modification, gibberellin biosynthesis, and a comprehensive group of nutrient (N, P, and S) uptake and assimilation genes, suggesting that NIN's recruitment to nodulation was based on its role as a growth module, a role shared with other NIN-Like Proteins. The expression of jasmonic acid genes in nin suggests the involvement of NIN in the resolution of growth versus defense outcomes. We find that the regulation of the growth module component Nodulation Pectate Lyase by NIN, and its function in rhizobial infection, are conserved in hologalegina legumes, highlighting its recruitment as a major event in the evolution of nodulation. We find that Nodulation Pectate Lyase is secreted to the infection chamber and the lumen of the infection thread. Gene network analysis using the transcription factor mutants for ERF Required for Nodulation1 and Nuclear Factor-Y Subunit A1 confirms hierarchical control of NIN over Nuclear Factor-Y Subunit A1 and shows that ERF Required for Nodulation1 acts independently to control infection. We conclude that while NIN shares functions with other NIN-Like Proteins, the conscription of key infection genes to NIN's control has made it a central regulatory hub for rhizobial infection.

Spatio-temporal expression dynamics differ between homologues of flowering time genes in the allopolyploid Brassica napus.

Polyploidy is a recurrent feature of eukaryotic evolution and has been linked to increases in complexity, adaptive radiation and speciation. Within angiosperms such events have occurred repeatedly in many plant lineages. Here we investigate the retention and spatio-temporal expression dynamics of duplicated genes predicted to regulate the floral transition in Brassica napus (oilseed rape, OSR). We show that flowering time genes are preferentially retained relative to other genes in the OSR genome. Using a transcriptome time series in two tissues (leaf and shoot apex) across development we show that 67% of these retained flowering time genes are expressed. Furthermore, between 64% (leaf) and 74% (shoot apex) of the retained gene homologues show diverged expression patterns relative to each other across development, suggesting neo- or subfunctionalization. A case study of homologues of the shoot meristem identity gene TFL1 reveals differences in cis-regulatory elements that could explain this divergence. Such differences in the expression dynamics of duplicated genes highlight the challenges involved in translating gene regulatory networks from diploid model systems to more complex polyploid crop species.

Discovery and Biosynthesis of the Antibiotic Bicyclomycin in Distantly Related Bacterial Classes.

Bicyclomycin (BCM) is a clinically promising antibiotic that is biosynthesized by Streptomyces cinnamoneus DSM 41675. BCM is structurally characterized by a core cyclo(l-Ile-l-Leu) 2,5-diketopiperazine (DKP) that is extensively oxidized. Here, we identify the BCM biosynthetic gene cluster, which shows that the core of BCM is biosynthesized by a cyclodipeptide synthase, and the oxidative modifications are introduced by five 2-oxoglutarate-dependent dioxygenases and one cytochrome P450 monooxygenase. The discovery of the gene cluster enabled the identification of BCM pathways encoded by the genomes of hundreds of Pseudomonas aeruginosa isolates distributed globally, and heterologous expression of the pathway from P. aeruginosa SCV20265 demonstrated that the product is chemically identical to BCM produced by S. cinnamoneus Overall, putative BCM gene clusters have been found in at least seven genera spanning Actinobacteria and Proteobacteria (Alphaproteobacteria, Betaproteobacteria, and Gammaproteobacteria). This represents a rare example of horizontal gene transfer of an intact biosynthetic gene cluster across such distantly related bacteria, and we show that these gene clusters are almost always associated with mobile genetic elements.IMPORTANCE Bicyclomycin is the only natural product antibiotic that selectively inhibits the transcription termination factor Rho. This mechanism of action, combined with its proven biological safety and its activity against clinically relevant Gram-negative bacterial pathogens, makes it a very promising antibiotic candidate. Here, we report the identification of the bicyclomycin biosynthetic gene cluster in the known bicyclomycin-producing organism Streptomyces cinnamoneus, which will enable the engineered production of new bicyclomycin derivatives. The identification of this gene cluster also led to the discovery of hundreds of bicyclomycin pathways encoded in highly diverse bacteria, including in the opportunistic pathogen Pseudomonas aeruginosa This wide distribution of a complex biosynthetic pathway is very unusual and provides an insight into how a pathway for an antibiotic can be transferred between diverse bacteria.

Dissecting the complex regulation of lodging resistance in Brassica napus.

Lodging continues to be a major cause of yield loss in important crop species such as Brassica napus. Understanding the genetic regulation of lodging resistance is therefore of key interest to breeders worldwide. Current strategies aimed at minimising lodging risk involve the incorporation of dwarfing genes or the application of plant growth regulators. However, despite these efforts, lodging continues to be a persistent problem and it is therefore of high interest that novel, complimentary strategies for lodging control are implemented. One approach would be to focus on understanding the genetic properties underlying stem mechanical strength. With this in mind, we screened a training genetic diversity panel of B. napus accession for variation in stem mechanical strength and related traits. Using Associative Transcriptomics, we identified molecular markers for a suite of valuable traits. Using an independent test genetic diversity panel, we show that the methods employed are robust for identification of predictive markers. Furthermore, based on conserved synteny with Arabidopsis thaliana, we are able to provide a biological context to the marker associations detected and provide evidence for a role in pectin methylesterification in contributing to stem mechanical strength in Brassicaceae.

Auxin-Induced Modulation of ETTIN Activity Orchestrates Gene Expression in Arabidopsis.

The phytohormone auxin governs crucial developmental decisions throughout the plant life cycle. Auxin signaling is effectuated by auxin response factors (ARFs) whose activity is repressed by Aux/IAA proteins under low auxin levels, but relieved from repression when cellular auxin concentrations increase. ARF3/ETTIN (ETT) is a conserved noncanonical Arabidopsis thaliana ARF that adopts an alternative auxin-sensing mode of translating auxin levels into multiple transcriptional outcomes. However, a mechanistic model for how this auxin-dependent modulation of ETT activity regulates gene expression has not yet been elucidated. Here, we take a genome-wide approach to show how ETT controls developmental processes in the Arabidopsis shoot through its auxin-sensing property. Moreover, analysis of direct ETT targets suggests that ETT functions as a central node in coordinating auxin dynamics and plant development and reveals tight feedback regulation at both the transcriptional and protein-interaction levels. Finally, we present an example to demonstrate how auxin sensitivity of ETT-protein interactions can shape the composition of downstream transcriptomes to ensure specific developmental outcomes. These results show that direct effects of auxin on protein factors, such as ETT-TF complexes, comprise an important part of auxin biology and likely contribute to the vast number of biological processes affected by this simple molecule.

Awake1, an ABC-Type Transporter, Reveals an Essential Role for Suberin in the Control of Seed Dormancy.

The mother plant plays an important dynamic role in the control of dormancy of her progeny seed in response to environmental signals. In order to further understand the mechanisms by which this dormancy control takes place in Arabidopsis (Arabidopsis thaliana), we conducted a forward genetic screen to isolate mutants that fail to enter dormancy in response to variation in temperature during seed set. We show that, for the first of these mutants, designated awake1, the maternal allele is required for entry into strongly dormant states and that awake1 mutants show seed phenotypes shown previously to be associated with the loss of suberin in the seed. We identify awake1 as an allele of ABCG20, an ATP-binding cassette transporter-encoding gene required for the transport of fatty acids during suberin deposition, and show that further suberin-deficient mutants have seed dormancy defects. Seed coat suberin composition is affected by temperature during seed maturation, but this response appears to be independent of ABCG20. We conclude that seed coat suberin is essential for seed dormancy imposition by low temperature and that the exclusion of oxygen and water from the seed by the suberin and tannin layers is important for dormancy imposition.

Elucidation of the genetic basis of variation for stem strength characteristics in bread wheat by Associative Transcriptomics.

BACKGROUND: The current approach to reducing the tendency for wheat grown under high fertilizer conditions to collapse (lodge) under the weight of its grain is based on reducing stem height via the introduction of Rht genes. However, these reduce the yield of straw (itself an important commodity) and introduce other undesirable characteristics. Identification of alternative height-control loci is therefore of key interest. In addition, the improvement of stem mechanical strength provides a further way through which lodging can be reduced. RESULTS: To investigate the prospects for genetic alternatives to Rht, we assessed variation for plant height and stem strength properties in a training genetic diversity panel of 100 wheat accessions fixed for Rht. Using mRNAseq data derived from RNA purified from leaves, functional genotypes were developed for the panel comprising 42,066 Single Nucleotide Polymorphism (SNP) markers and 94,060 Gene Expression Markers (GEMs). In the first application in wheat of the recently-developed method of Associative Transcriptomics, we identified associations between trait variation and both SNPs and GEMs. Analysis of marker-trait associations revealed candidates for the causative genes underlying the trait variation, implicating xylan acetylation and the COP9 signalosome as contributing to stem strength and auxin in the control of the observed variation for plant height. Predictive capabilities of key markers for stem strength were validated using a test genetic diversity panel of 30 further wheat accessions. CONCLUSIONS: This work illustrates the power of Associative Transcriptomics for the exploration of complex traits of high agronomic importance in wheat. The careful selection of genotypes included in the analysis, allowed for high resolution mapping of novel trait-controlling loci in this staple crop. The use of Gene Expression markers coupled with the more traditional sequence-based markers, provides the power required to understand the biological context of the marker-trait associations observed. This not only adds to the wealth of knowledge that we strive to accumulate regarding gene function and plant adaptation, but also provides breeders with the information required to make more informed decisions regarding the potential consequences of incorporating the use of particular markers into future breeding programmes.

The Starch Granule-Associated Protein EARLY STARVATION1 Is Required for the Control of Starch Degradation in Arabidopsis thaliana Leaves.

To uncover components of the mechanism that adjusts the rate of leaf starch degradation to the length of the night, we devised a screen for mutant Arabidopsis thaliana plants in which starch reserves are prematurely exhausted. The mutation in one such mutant, named early starvation1 (esv1), eliminates a previously uncharacterized protein. Starch in mutant leaves is degraded rapidly and in a nonlinear fashion, so that reserves are exhausted 2 h prior to dawn. The ESV1 protein and a similar uncharacterized Arabidopsis protein (named Like ESV1 [LESV]) are located in the chloroplast stroma and are also bound into starch granules. The region of highest similarity between the two proteins contains a series of near-repeated motifs rich in tryptophan. Both proteins are conserved throughout starch-synthesizing organisms, from angiosperms and monocots to green algae. Analysis of transgenic plants lacking or overexpressing ESV1 or LESV, and of double mutants lacking ESV1 and another protein necessary for starch degradation, leads us to propose that these proteins function in the organization of the starch granule matrix. We argue that their misexpression affects starch degradation indirectly, by altering matrix organization and, thus, accessibility of starch polymers to starch-degrading enzymes.

Delineation of metabolic gene clusters in plant genomes by chromatin signatures.

Plants are a tremendous source of diverse chemicals, including many natural product-derived drugs. It has recently become apparent that the genes for the biosynthesis of numerous different types of plant natural products are organized as metabolic gene clusters, thereby unveiling a highly unusual form of plant genome architecture and offering novel avenues for discovery and exploitation of plant specialized metabolism. Here we show that these clustered pathways are characterized by distinct chromatin signatures of histone 3 lysine trimethylation (H3K27me3) and histone 2 variant H2A.Z, associated with cluster repression and activation, respectively, and represent discrete windows of co-regulation in the genome. We further demonstrate that knowledge of these chromatin signatures along with chromatin mutants can be used to mine genomes for cluster discovery. The roles of H3K27me3 and H2A.Z in repression and activation of single genes in plants are well known. However, our discovery of highly localized operon-like co-regulated regions of chromatin modification is unprecedented in plants. Our findings raise intriguing parallels with groups of physically linked multi-gene complexes in animals and with clustered pathways for specialized metabolism in filamentous fungi.

Molecular markers for tolerance of European ash (Fraxinus excelsior) to dieback disease identified using Associative Transcriptomics.

Tree disease epidemics are a global problem, impacting food security, biodiversity and national economies. The potential for conservation and breeding in trees is hampered by complex genomes and long lifecycles, with most species lacking genomic resources. The European Ash tree Fraxinus excelsior is being devastated by the fungal pathogen Hymenoscyphus fraxineus, which causes ash dieback disease. Taking this system as an example and utilizing Associative Transcriptomics for the first time in a plant pathology study, we discovered gene sequence and gene expression variants across a genetic diversity panel scored for disease symptoms and identified markers strongly associated with canopy damage in infected trees. Using these markers we predicted phenotypes in a test panel of unrelated trees, successfully identifying individuals with a low level of susceptibility to the disease. Co-expression analysis suggested that pre-priming of defence responses may underlie reduced susceptibility to ash dieback.

Genome distribution of differential homoeologue contributions to leaf gene expression in bread wheat.

Using a combination of de novo transcriptome assembly, a newly developed 9495-marker transcriptome SNP genetic linkage map and comparative genomics approaches, we developed an ordered set of nonredundant transcripts for each of the subgenomes of hexaploid wheat: A (47 160 unigenes), B (59 663 unigenes) and D (40 588 unigenes). We used these as reference sequences against which to map Illumina mRNA-Seq reads derived from young leaf tissue. Transcript abundance was quantified for each unigene. Using a three-way reciprocal BLAST approach, 15 527 triplet sets of homoeologues (one from each genome) were identified. Differential expression (P < 0.05) was identified for 5248 unigenes, with 2906 represented at greater abundance than their two homoeologues and 2342 represented at lower abundance than their two homoeologues. Analysis of gene ontology terms revealed no biases between homoeologues. There was no evidence of genomewide dominance effects, rather the more highly transcribed individual genes were distributed throughout all three genomes. Transcriptome display tile plot, a visualization approach based on CMYK colour space, was developed and used to assess the genome for regions of skewed homoeologue transcript abundance. Extensive striation was revealed, indicative of many small regions of genome dominance (transcripts of homoeologues from one genome more abundant than the others) and many larger regions of genome repression (transcripts of homoeologues from one genome less abundant than the others).