Andrea L. Harper

Associative transcriptomics of traits in the polyploid crop species Brassica napus

Association genetics can quickly and efficiently delineate regions of the genome that control traits and provide markers to accelerate breeding by marker-assisted selection. But most crops are polyploid, making it difficult to identify the required markers and to assemble a genome sequence to order those markers. To circumvent this difficulty, we developed associative transcriptomics, which uses transcriptome sequencing to identify and score molecular markers representing variation in both gene sequences and gene expression, and correlate this with trait variation. Applying the method in the recently formed tetraploid crop Brassica napus, we identified genomic deletions that underlie two quantitative trait loci for glucosinolate content of seeds. The deleted regions contained orthologs of the transcription factor HAG1 (At5g61420), which controls aliphatic glucosinolate biosynthesis in Arabidopsis thaliana. This approach facilitates the application of association genetics in a broad range of crops, even those with complex genomes.

Identification of the pollen self-incompatibility determinant in Papaver rhoeas

Higher plants produce seed through pollination, using specific interactions between pollen and pistil. Self-incompatibility is an important mechanism used in many species to prevent inbreeding; it is controlled by a multi-allelic S locus. ‘Self’ (incompatible) pollen is discriminated from ‘non-self’ (compatible) pollen by interaction of pollen and pistil S locus components, and is subsequently inhibited. In Papaver rhoeas, the pistil S locus product is a small protein that interacts with incompatible pollen, triggering a Ca2+-dependent signalling network, resulting in pollen inhibition and programmed cell death. Here we have cloned three alleles of a highly polymorphic pollen-expressed gene, PrpS (Papaver rhoeas pollen S), from Papaver and provide evidence that this encodes the pollen S locus determinant. PrpS is a single-copy gene linked to the pistil S gene (currently called S, but referred to hereafter as PrsS for Papaver rhoeas stigma S determinant). Sequence analysis indicates that PrsS and PrpS are equally ancient and probably co-evolved. PrpS encodes a novel ∼20-kDa protein. Consistent with predictions that it is a transmembrane protein, PrpS is associated with the plasma membrane. We show that a predicted extracellular loop segment of PrpS interacts with PrsS and, using PrpS antisense oligonucleotides, we demonstrate that PrpS is involved in S-specific inhibition of incompatible pollen. Identification of PrpS represents a major advance in our understanding of the Papaver self-incompatibility system. As a novel cell–cell recognition determinant it contributes to the available information concerning the origins and evolution of cell–cell recognition systems involved in discrimination between self and non-self, which also include histocompatibility systems in primitive chordates and vertebrates.

Genome-Wide Association Study Dissects the Genetic Architecture of Seed Weight and Seed Quality in Rapeseed (Brassica napus L.)

Association mapping can quickly and efficiently dissect complex agronomic traits. Rapeseed is one of the most economically important polyploid oil crops, although its genome sequence is not yet published. In this study, a recently developed 60K Brassica Infinium® SNP array was used to analyse an association panel with 472 accessions. The single-nucleotide polymorphisms (SNPs) of the array were in silico mapped using ‘pseudomolecules’ representative of the genome of rapeseed to establish their hypothetical order and to perform association mapping of seed weight and seed quality. As a result, two significant associations on A8 and C3 of Brassica napus were detected for erucic acid content, and the peak SNPs were found to be only 233 and 128 kb away from the key genes BnaA.FAE1 and BnaC.FAE1. BnaA.FAE1 was also identified to be significantly associated with the oil content. Orthologues of Arabidopsis thaliana HAG1 were identified close to four clusters of SNPs associated with glucosinolate content on A9, C2, C7 and C9. For seed weight, we detected two association signals on A7 and A9, which were consistent with previous studies of quantitative trait loci mapping. The results indicate that our association mapping approach is suitable for fine mapping of the complex traits in rapeseed.

Divergent Regulation of Terpenoid Metabolism in the Trichomes of Wild and Cultivated Tomato Species

The diversification of chemical production in glandular trichomes is important in the development of resistance against pathogens and pests in two species of tomato. We have used genetic and genomic approaches to uncover some of the biochemical and molecular mechanisms that underlie the divergence in trichome metabolism between the wild species Solanum habrochaites LA1777 and its cultivated relative, Solanum lycopersicum. LA1777 produces high amounts of insecticidal sesquiterpene carboxylic acids (SCAs), whereas cultivated tomatoes lack SCAs and are more susceptible to pests. We show that trichomes of the two species have nearly opposite terpenoid profiles, consisting mainly of monoterpenes and low levels of sesquiterpenes in S. lycopersicum and mainly of SCAs and very low monoterpene levels in LA1777. The accumulation patterns of these terpenoids are different during development, in contrast to the developmental expression profiles of terpenoid pathway genes, which are similar in the two species, but they do not correlate in either case with terpenoid accumulation. However, our data suggest that the accumulation of monoterpenes in S. lycopersicum and major sesquiterpenes in LA1777 are linked both genetically and biochemically. Metabolite analyses after targeted gene silencing, inhibitor treatments, and precursor feeding all show that sesquiterpene biosynthesis relies mainly on products from the plastidic 2-C-methyl-d-erythritol-4-phosphate pathway in LA1777 but less so in the cultivated species. Furthermore, two classes of sesquiterpenes produced by the wild species may be synthesized from distinct pools of precursors via cytosolic and plastidial cyclases. However, highly trichome-expressed sesquiterpene cyclase-like enzymes were ruled out as being involved in the production of major LA1777 sesquiterpenes.

Transcriptomic and Reverse Genetic Analysesof Branched-Chain Fatty Acid and Acyl Sugar Production in Solanum pennellii and Nicotiana benthamiana

Acyl sugars containing branched-chain fatty acids (BCFAs) are exuded by glandular trichomes of many species in Solanaceae, having an important defensive role against insects. From isotope-feeding studies, two modes of BCFA elongation have been proposed: (1) fatty acid synthase-mediated two-carbon elongation in the high acyl sugar-producing tomato species Solanum pennellii and Datura metel; and (2) α-keto acid elongation-mediated one-carbon increments in several tobacco (Nicotiana) species and a Petunia species. To investigate the molecular mechanisms underlying BCFAs and acyl sugar production in trichomes, we have taken a comparative genomic approach to identify critical enzymatic steps followed by gene silencing and metabolite analysis in S. pennellii and Nicotiana benthamiana. Our study verified the existence of distinct mechanisms of acyl sugar synthesis in Solanaceae. From microarray analyses, genes associated with α-keto acid elongation were found to be among the most strongly expressed in N. benthamiana trichomes only, supporting this model in tobacco species. Genes encoding components of the branched-chain keto-acid dehydrogenase complex were expressed at particularly high levels in trichomes of both species, and we show using virus-induced gene silencing that they are required for BCFA production in both cases and for acyl sugar synthesis in N. benthamiana. Functional analysis by down-regulation of specific KAS I genes and cerulenin inhibition indicated the involvement of the fatty acid synthase complex in BCFA production in S. pennellii. In summary, our study highlights both conserved and divergent mechanisms in the production of important defense compounds in Solanaceae and defines potential targets for engineering acyl sugar production in plants for improved pest tolerance.

Genome sequence and genetic diversity of European ash trees

Ash trees (genus Fraxinus, family Oleaceae) are widespread throughout the Northern Hemisphere, but are being devastated in Europe by the fungus Hymenoscyphus fraxineus, causing ash dieback, and in North America by the herbivorous beetle Agrilus planipennis. Here we sequence the genome of a low-heterozygosity Fraxinus excelsior tree from Gloucestershire, UK, annotating 38,852 protein-coding genes of which 25% appear ash specific when compared with the genomes of ten other plant species. Analyses of paralogous genes suggest a whole-genome duplication shared with olive (Olea europaea, Oleaceae). We also re-sequence 37 F. excelsior trees from Europe, finding evidence for apparent long-term decline in effective population size. Using our reference sequence, we re-analyse association transcriptomic data, yielding improved markers for reduced susceptibility to ash dieback. Surveys of these markers in British populations suggest that reduced susceptibility to ash dieback may be more widespread in Great Britain than in Denmark. We also present evidence that susceptibility of trees to H. fraxineus is associated with their iridoid glycoside levels. This rapid, integrated, multidisciplinary research response to an emerging health threat in a non-model organism opens the way for mitigation of the epidemic.

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.

Associative Transcriptomics Study Dissects the Genetic Architecture of Seed Glucosinolate Content in Brassica napus

Breeding new varieties with low seed glucosinolate (GS) concentrations has long been a prime target in Brassica napus. In this study, a novel association mapping methodology termed ‘associative transcriptomics’ (AT) was applied to a panel of 101 B. napus lines to define genetic regions and also candidate genes controlling total seed GS contents. Over 100,000 informative single-nucleotide polymorphisms (SNPs) and gene expression markers (GEMs) were developed for AT analysis, which led to the identification of 10 SNP and 7 GEM association peaks. Within these peaks, 26 genes were inferred to be involved in GS biosynthesis. A weighted gene co-expression network analysis provided additional 40 candidate genes. The transcript abundance in leaves of two candidate genes, BnaA.GTR2a located on chromosome A2 and BnaC.HAG3b on C9, was correlated with seed GS content, explaining 18.8 and 16.8% of phenotypic variation, respectively. Resequencing of genomic regions revealed six new SNPs in BnaA.GTR2a and four insertions or deletions in BnaC.HAG3b. These deletion polymorphisms were then successfully converted into polymerase chain reaction–based diagnostic markers that can, due to high linkage disequilibrium observed in these regions of the genome, be used for marker-assisted breeding for low seed GS lines.

DYNAMICS OF DRIFT, GENE FLOW, AND SELECTION DURING SPECIATION IN SILENE

The mechanics of speciation with gene flow are still unclear. Disparity among genes in population differentiation (FST) between diverging species is often interpreted as evidence for semipermeable species boundaries, with selection preventing “key” genes from introgressing despite ongoing gene flow. However, FST can remain high before it reaches equilibrium between the lineage sorting of species divergence and the homogenizing effects of gene flow (via secondary contact). Thus, when interpreting FST, the dynamics of drift, gene flow, and selection need to be taken into account. We illustrate this view with a multigenic analyses of gene flow and selection in three closely related Silene species, S. latifolia, S. dioica, and S. diclinis. We report that although S. diclinis appears to have evolved in allopatry, isolation with (bidirectional) gene flow between S. latifolia and S. dioica is likely, perhaps as a result of parapatric speciation followed by more extensive sympatry. Interestingly, we detected the signatures of apparently independent instances of positive selection at the same locus in S. latifolia and S. dioica. Despite gene flow between the species, the adaptive alleles have not crossed the species boundary, suggesting that this gene has independently undergone species‐specific (diversifying or parallel) selection.

Variation in the chemical composition of wheat straw: the role of tissue ratio and composition

BackgroundWheat straw is an attractive substrate for second generation ethanol production because it will complement and augment wheat production rather than competing with food production. However, like other sources of lignocellulosic biomass, even from a single species, it is heterogeneous in nature due to the different tissues and cell types, and this has implications for saccharification efficiency. The aim of this study has been to use Fourier transform infrared (FTIR) spectroscopy and Partial least squares (PLS) modelling to rapidly screen wheat cultivars for the levels of component tissues, the carbohydrate composition and lignin content, and the levels of simple cross-linking phenolics such as ferulic and diferulic acids.ResultsFTIR spectroscopy and PLS modelling was used to analyze the tissue and chemical composition of wheat straw biomass. Predictive models were developed to evaluate the variability in the concentrations of the cell wall sugars, cell wall phenolics and acid-insoluble lignin. Models for the main sugars, phenolics and lignin were validated and then used to evaluate the variation in total biomass composition across 90 cultivars of wheat grown over two seasons.ConclusionsWhilst carbohydrate and lignin components varied across the varieties, this mainly reflected differences in the ratios of the component tissues rather than differences in the composition of those tissues. Further analysis indicated that on a mol% basis, relative levels of sugars within the tissues varied to only a small degree. There were no clear associations between simple phenolics and tissues. The results provide a basis for improving biomass quality for biofuels production through selection of cultivars with appropriate tissue ratios.