An alfalfa MYB-like transcriptional factor MsMYBH positively regulates alfalfa seedling drought resistance and undergoes MsWAV3-mediated degradation.

Drought is a major threat to alfalfa (Medicago sativa L.) production. The discovery of important alfalfa genes regulating drought response will facilitate breeding for drought-resistant alfalfa cultivars. Here, we report a genome-wide association study of drought resistance in alfalfa. We identified and functionally characterized an MYB-like transcription factor gene (MsMYBH), which increases the drought resistance in alfalfa. Compared with the wild-types, the biomass and forage quality were enhanced in MsMYBH overexpressed plants. Combined RNA-seq, proteomics and chromatin immunoprecipitation analysis showed that MsMYBH can directly bind to the promoters of MsMCP1, MsMCP2, MsPRX1A and MsCARCAB to improve their expression. The outcomes of such interactions include better water balance, high photosynthetic efficiency and scavenge excess H2O2 in response to drought. Furthermore, an E3 ubiquitin ligase (MsWAV3) was found to induce MsMYBH degradation under long-term drought, via the 26S proteasome pathway. Furthermore, variable-number tandem repeats in MsMYBH promoter were characterized among a collection of germplasms, and the variation is associated with promoter activity. Collectively, our findings shed light on the functions of MsMYBH and provide a pivotal gene that could be leveraged for breeding drought-resistant alfalfa. This discovery also offers new insights into the mechanisms of drought resistance in alfalfa.

Genome and systems biology of Melilotus albus provides insights into coumarins biosynthesis.

Melilotus species are used as green manure and rotation crops worldwide and contain abundant pharmacologically active coumarins. However, there is a paucity of information on its genome and coumarin production and function. Here, we reported a chromosome-scale assembly of Melilotus albus genome with 1.04 Gb in eight chromosomes, containing 71.42% repetitive elements. Long terminal repeat retrotransposon bursts coincided with declining of population sizes during the Quaternary glaciation. Resequencing of 94 accessions enabled insights into genetic diversity, population structure, and introgression. Melilotus officinalis had relatively larger genetic diversity than that of M. albus. The introgression existed between M. officinalis group and M. albus group, and gene flows was from M. albus to M. officinalis. Selection sweep analysis identified candidate genes associated with flower colour and coumarin biosynthesis. Combining genomics, BSA, transcriptomics, metabolomics, and biochemistry, we identified a ß-glucosidase (BGLU) gene cluster contributing to coumarin biosynthesis. MaBGLU1 function was verified by overexpression in M. albus, heterologous expression in Escherichia coli, and substrate feeding, revealing its role in scopoletin (coumarin derivative) production and showing that nonsynonymous variation drives BGLU enzyme activity divergence in Melilotus. Our work will accelerate the understanding of biologically active coumarins and their biosynthetic pathways, and contribute to genomics-enabled Melilotus breeding.

The genome of Cleistogenes songorica provides a blueprint for functional dissection of dimorphic flower differentiation and drought adaptability.

Cleistogenes songorica (2n = 4x = 40) is a desert grass with a unique dimorphic flowering mechanism and an ability to survive extreme drought. Little is known about the genetics underlying drought tolerance and its reproductive adaptability. Here, we sequenced and assembled a high-quality chromosome-level C. songorica genome (contig N50 = 21.28 Mb). Complete assemblies of all telomeres, and of ten chromosomes were derived. C. songorica underwent a recent tetraploidization (~19 million years ago) and four major chromosomal rearrangements. Expanded genes were significantly enriched in fatty acid elongation, phenylpropanoid biosynthesis, starch and sucrose metabolism, and circadian rhythm pathways. By comparative transcriptomic analysis we found that conserved drought tolerance related genes were expanded. Transcription of CsMYB genes was associated with differential development of chasmogamous and cleistogamous flowers, as well as drought tolerance. Furthermore, we found that regulation modules encompassing miRNA, transcription factors and target genes are involved in dimorphic flower development, validated by overexpression of CsAP2_9 and its targeted miR172 in rice. Our findings enable further understanding of the mechanisms of drought tolerance and flowering in C. songorica, and provide new insights into the adaptability of native grass species in evolution, along with potential resources for trait improvement in agronomically important species.

Unveiling the Complexity of Red Clover (Trifolium pratense L.) Transcriptome and Transcriptional Regulation of Isoflavonoid Biosynthesis Using Integrated Long- and Short-Read RNAseq.

Red clover (Trifolium pratense L.) is used as forage and contains a high level of isoflavonoids. Although isoflavonoids in red clover were discovered a long time ago, the transcriptional regulation of isoflavonoid biosynthesis is virtually unknown because of the lack of accurate and comprehensive characterization of the transcriptome. Here, we used a combination of long-read (PacBio Iso-Seq) and short-read (Illumina) RNAseq sequencing to develop a more comprehensive full-length transcriptome in four tissues (root, stem, leaf, and flower) and to identify transcription factors possibly involved in isoflavonoid biosynthesis in red clover. Overall, we obtained 50,922 isoforms, including 19,860 known genes and 2817 novel isoforms based on the annotation of RefGen Tp_v2.0. We also found 1843 long non-coding RNAs, 1625 fusion genes, and 34,612 alternatively spliced events, with some transcript isoforms validated experimentally. A total of 16,734 differentially expressed genes were identified in the four tissues, including 43 isoflavonoid-biosynthesis-related genes, such as stem-specific expressed TpPAL, TpC4H, and Tp4CL and root-specific expressed TpCHS, TpCHI1, and TpIFS. Further, weighted gene co-expression network analysis and a targeted compound assay were combined to investigate the association between the isoflavonoid content and the transcription factors expression in the four tissues. Twelve transcription factors were identified as key genes for isoflavonoid biosynthesis. Among these transcription factors, the overexpression of TpMYB30 or TpRSM1-2 significantly increased the isoflavonoid content in tobacco. In particular, the glycitin was increased by 50-100 times in the plants overexpressing TpRSM1-2, in comparison to that in the WT plants. Our study provides a comprehensive and accurate annotation of the red clover transcriptome and candidate genes to improve isoflavonoid biosynthesis and accelerate research into molecular breeding in red clover or other crops.

Low pyrrolizidine alkaloid levels in perennial ryegrass is associated with the absence of a homospermidine synthase gene.

BACKGROUND: Pyrrolizidine alkaloids (PAs) are a class of secondary metabolites that function as feeding deterrents in a range of different plant species. In perennial ryegrass (Lolium perenne L.) the only PAs that have been identified are the thesinine-rhamnoside group, which displays significant genetic variation. Homospermidine synthase (HSS) has evolved from deoxyhypusine synthase (DHS) and catalyses the first step in the PA pathway, making it a key candidate for the investigation of genes influencing observed PA trait variation. RESULTS: During PCR amplification and sequence analysis of DHS we identified two putative HSS genes in perennial ryegrass. One of the genes (LpHSS1) was absent in some perennial ryegrass plants. Thesinine-rhamnoside levels were measured using liquid chromatography coupled with mass spectrometry in a diverse association mapping population, consisting of 693 plants free of fungal endophytic symbionts. Association tests that accounted for population structure identified a significant association of absence of the LpHSS1 gene with lower levels of thesinine-rhamnoside PAs. HSS-like gene sequences were identified for other grass species of the Poaceae, including tall fescue, wheat, maize and sorghum. CONCLUSION: HSS is situated at the crucial first step in the PA pathway making it an important candidate gene for investigation of involvement in PA phenotypic variation. In this study, PA level in perennial ryegrass was strongly associated with the presence or absence of the LpHSS1 gene. A genetic marker, developed for the presence/absence of LpHSS1, may be used for marker-assisted breeding to either lower or increase PAs in breeding populations of perennial or Italian ryegrass to investigate a potential role in the deterrence of herbivore pests. The presence of HSS-like genes in several other Poaceae species suggests that PA biosynthesis may occur in plant family members beyond perennial ryegrass and tall fescue and identifies a potential route for manipulating PA levels.

Predictive ability of genomic selection models in a multi-population perennial ryegrass training set using genotyping-by-sequencing.

KEY MESSAGE: Genomic prediction models for multi-year dry matter yield, via genotyping-by-sequencing in a composite training set, demonstrate potential for genetic gain improvement through within-half sibling family selection. Perennial ryegrass (Lolium perenne L.) is a key source of nutrition for ruminant livestock in temperate environments worldwide. Higher seasonal and annual yield of herbage dry matter (DMY) is a principal breeding objective but the historical realised rate of genetic gain for DMY is modest. Genomic selection was investigated as a tool to enhance the rate of genetic gain. Genotyping-by-sequencing (GBS) was undertaken in a multi-population (MP) training set of five populations, phenotyped as half-sibling (HS) families in five environments over 2 years for mean herbage accumulation (HA), a measure of DMY potential. GBS using the ApeKI enzyme yielded 1.02 million single-nucleotide polymorphism (SNP) markers from a training set of n = 517. MP-based genomic prediction models for HA were effective in all five populations, cross-validation-predictive ability (PA) ranging from 0.07 to 0.43, by trait and target population, and 0.40-0.52 for days-to-heading. Best linear unbiased predictor (BLUP)-based prediction methods, including GBLUP with either a standard or a recently developed (KGD) relatedness estimation, were marginally superior or equal to ridge regression and random forest computational approaches. PA was principally an outcome of SNP modelling genetic relationships between training and validation sets, which may limit application for long-term genomic selection, due to PA decay. However, simulation using data from the training experiment indicated a twofold increase in genetic gain for HA, when applying a prediction model with moderate PA in a single selection cycle, by combining among-HS family selection, based on phenotype, with within-HS family selection using genomic prediction.

Metabolic changes and associated cytokinin signals in response to nitrate assimilation in roots and shoots of Lolium perenne.

The efficiency of inorganic nitrogen (N) assimilation is a critical component of fertilizer use by plants and of forage production in Lolium perenne, an important pasture species worldwide. We present a spatiotemporal description of nitrate use efficiency in terms of metabolic responses and carbohydrate remobilization, together with components of cytokinin signal transduction following nitrate addition to N-impoverished plants. Perennial ryegrass (L. perenne cv. Grasslands Nui) plants were grown for 10 weeks in unfertilized soil and then treated with nitrate (5 mM) hydroponically. Metabolomic analysis by gas chromatography-mass spectrometry and liquid chromatography-tandem mass spectrometry revealed a dynamic interaction between N and carbon metabolism over a week-long time course represented by the relative abundance of amino acids, tricarboxylic acid intermediates and stored water-soluble carbohydrates (WSCs). The initial response to N addition was characterized by a rapid remobilization of carbon stores from the low-molecular weight WSC, along with an increase in N content and assimilation into free amino acids. Subsequently, the shoot became the main source of carbon through remobilization of a large pool of high-molecular weight WSC. Associated quantification of cytokinin levels and expression profiling of putative cytokinin response regulator genes by quantitative reverse transcription polymerase chain reaction support a role for cytokinin in the mediation of the response to N addition in perennial ryegrass. The presence of high levels of cis-zeatin-type cytokinins is discussed in the context of hormonal homeostasis under the stress of steady-state N deficiency.

A novel family of cyclic oligopeptides derived from ribosomal peptide synthesis of an in planta-induced gene, gigA, in Epichloë endophytes of grasses.

Fungal endophytes belonging to the genus Epichloë form associations with temperate grasses belonging to the sub-family Poöideae that range from mutualistic through to pathogenic. We previously identified a novel endophyte gene (designated gigA for grass induced gene) that is one of the most abundantly expressed fungal transcripts in endophyte-infected grasses and which is distributed and highly expressed in a wide range of Epichloë grass associations. Molecular and biochemical analyses indicate that gigA encodes a small secreted protein containing an imperfect 27 amino acid repeat that includes a kexin protease cleavage site. Kexin processing of GigA liberates within the plant multiple related products, named here as epichloëcyclins, which we have demonstrated by MS/MS to be cyclic peptidic in nature. Gene deletion of gigA leads to the elimination of all epichloëcyclins with no conspicuous phenotypic impact on the host grass, suggesting a possible bioactive role. This is a further example of a ribosomal peptide synthetic (RiPS) pathway operating within the Ascomycetes, and is the first description of such a pathway from a mutualistic symbiotic fungus from this Phylum.

Non-targeted multimodal metabolomics data from ovine rumen fluid fractions.

From an animal health perspective, our understanding of the metabolites in rumen fluid across different host species is poorly understood. Here, we present a metabolomic data set generated using hydrophilic interaction liquid chromatography and semi-polar (C18) chromatography methods coupled to high-resolution mass spectrometry of fractionated ovine rumen samples.

Molecular insights into quality and authentication of sheep meat from proteomics and metabolomics.

Sheep meat (encompassing lamb, hogget and mutton) is an important source of animal protein in many countries, with a unique flavour and sensory profile compared to other red meats. Flavour, colour and texture are the key quality attributes contributing to consumer liking of sheep meat. Over the last decades, various factors from 'farm to fork', including production system (e.g., age, breed, feeding regimes, sex, pre-slaughter stress, and carcass suspension), post-mortem manipulation and processing (e.g., electrical stimulation, ageing, packaging types, and chilled and frozen storage) have been identified as influencing different aspects of sheep meat quality. However conventional meat-quality assessment tools are not able to elucidate the underlying mechanisms and pathways for quality variations. Advances in broad-based analytical techniques have offered opportunities to obtain deeper insights into the molecular changes of sheep meat which may become biomarkers for specific variations in quality traits and meat authenticity. This review provides an overview on how omics techniques, especially proteomics (including peptidomics) and metabolomics (including lipidomics and volatilomics) are applied to elucidate the variations in sheep meat quality, mainly in loin muscles, focusing on colour, texture and flavour, and as tools for authentication. SIGNIFICANCE: From this review, we observed that attempts have been made to utilise proteomics and metabolomics techniques on sheep meat products for elucidating pathways of quality variations due to various factors. For instance, the improvement of colour stability and tenderness could be associated with the changes to glycolysis, energy metabolism and endogenous antioxidant capacity. Several studies identify proteolysis as being important, but potentially conflicting for quality as the enhanced proteolysis improves tenderness and flavour, while reducing colour stability. The use of multiple analytical methods e.g., lipidomics, metabolomics, and volatilomics, detects a wider range of flavour precursors (including both water and lipid soluble compounds) that underlie the possible pathways for sheep meat flavour evolution. The technological advancement in omics (e.g., direct analysis-mass spectrometry) could make analysis of the proteins, lipids and metabolites in sheep meat routine, as well as enhance the confidence in quality determination and molecular-based assurance of meat authenticity.

Untargeted Multimodal Metabolomics Investigation of the Haemonchus contortus Exsheathment Secretome.

In nematodes that invade the gastro-intestinal tract of the ruminant, the process of larval exsheathment marks the transition from the free-living to the parasitic stages of these parasites. To investigate the secretome associated with larval exsheathment, a closed in vitro system that effectively reproduces the two basic components of an anaerobic rumen environment (CO2 and 39 °C) was developed to trigger exsheathment in one of the most pathogenic and model gastrointestinal parasitic nematodes, Haemonchus contortus (barber's pole worm). This study reports the use of multimodal untargeted metabolomics and lipidomics methodologies to identify the metabolic signatures and compounds secreted during in vitro larval exsheathment in the H. contortus infective third-stage larva (iL3). A combination of statistical and chemoinformatic analyses using three analytical platforms revealed a panel of metabolites detected post exsheathment and associated with amino acids, purines, as well as select organic compounds. The major lipid classes identified by the non-targeted lipidomics method applied were lysophosphatidylglycerols, diglycerides, fatty acyls, glycerophospholipids, and a triglyceride. The identified metabolites may serve as metabolic signatures to improve tractability of parasitic nematodes for characterizing small molecule host-parasite interactions related to pathogenesis, vaccine and drug design, as well as the discovery of metabolic biomarkers.

Identification of urinary biomarkers of colon inflammation in IL10-/- mice using Short-Column LCMS metabolomics.

The interleukin-10-deficient (IL10(-/-)) mouse develops colon inflammation in response to normal intestinal microflora and has been used as a model of Crohn's disease. Short-Column LCMS metabolite profiling of urine from IL10(-/-) and wild-type (WT) mice was used, in two independent experiments, to identify mass spectral ions differing in intensity between these two genotypes. Three differential metabolites were identified as xanthurenic acid and as the glucuronides of xanthurenic acid and of α-CEHC (2,5,7,8-tetramethyl-2-(2'-carboxyethyl)-6-hydroxychroman). The significance of several differential metabolites as potential biomarkers of colon inflammation was evaluated in an experiment which compared metabolite concentrations in IL10(-/-) and WT mice housed, either under conventional conditions and dosed with intestinal microflora, or maintained under specific pathogen-free (SPF) conditions. Concentrations of xanthurenic acid, α-CEHC glucuronide, and an unidentified metabolite m/z 495(-)/497(+) were associated with the degree of inflammation in IL10(-/-) mice and may prove useful as biomarkers of colon inflammation.

Relationship between virulence gene profiles of atypical enteropathogenic Escherichia coli and Shiga toxin-producing E. coli isolates from cattle and sheep in New Zealand.

Virulence gene profiles of atypical enteropathogenic Escherichia coli (aEPEC) and Shiga toxin-producing E. coli (STEC) from cattle, sheep, and humans were examined to determine the relationship between pathotypes. Shared virulence factors (intimin, EHEC hemolysin, serine protease, and a type II secretion system) were identified, suggesting a dynamic evolutionary relationship between aEPEC and STEC.

Multi-Trait Genomic Prediction Improves Predictive Ability for Dry Matter Yield and Water-Soluble Carbohydrates in Perennial Ryegrass.

In perennial ryegrass (Lolium perenne L), annual and seasonal dry matter yield (DMY) and nutritive quality of herbage are high-priority traits targeted for improvement through selective breeding. Genomic prediction (GP) has proven to be a valuable tool for improving complex traits and may be further enhanced through the use of multi-trait (MT) prediction models. In this study, we evaluated the relative performance of MT prediction models to improve predictive ability for DMY and key nutritive quality traits, using two different training populations (TP1, n = 463 and TP2, n = 517) phenotyped at multiple locations. MT models outperformed single-trait (ST) models by 24% to 59% for DMY and 67% to 105% for nutritive quality traits, such as low, high, and total WSC, when a correlated secondary trait was included in both the training and test set (MT-CV2) or in the test set alone (MT-CV3) (trait-assisted genomic selection). However, when a secondary trait was included in training set and not the test set (MT-CV1), the predictive ability was not statistically significant (p > 0.05) compared to the ST model. We evaluated the impact of training set size when using a MT-CV2 model. Using a highly correlated trait (rg = 0.88) as the secondary trait in the MT-CV2 model, there was no loss in predictive ability for DMY even when the training set was reduced to 50% of its original size. In contrast, using a weakly correlated secondary trait (rg = 0.56) in the MT-CV2 model, predictive ability began to decline when the training set size was reduced by only 11% from its original size. Using a ST model, genomic predictive ability in a population unrelated to the training set was poor (rp = -0.06). However, when using an MT-CV2 model, the predictive ability was positive and high (rp = 0.76) for the same population. Our results demonstrate the first assessment of MT models in forage species and illustrate the prospects of using MT genomic selection in forages, and other outcrossing plant species, to accelerate genetic gains for complex agronomical traits, such as DMY and nutritive quality characteristics.

Genomic Predictive Ability for Foliar Nutritive Traits in Perennial Ryegrass.

Forage nutritive value impacts animal nutrition, which underpins livestock productivity, reproduction and health. Genetic improvement for nutritive traits in perennial ryegrass has been limited, as they are typically expensive and time-consuming to measure through conventional methods. Genomic selection is appropriate for such complex and expensive traits, enabling cost-effective prediction of breeding values using genome-wide markers. The aims of the present study were to assess the potential of genomic selection for a range of nutritive traits in a multi-population training set, and to quantify contributions of family, location and family-by-location variance components to trait variation and heritability for nutritive traits. The training set consisted of a total of 517 half-sibling (half-sib) families, from five advanced breeding populations, evaluated in two distinct New Zealand grazing environments. Autumn-harvested samples were analyzed for 18 nutritive traits and maternal parents of the half-sib families were genotyped using genotyping-by-sequencing. Significant (P < 0.05) family variance was detected for all nutritive traits and genomic heritability (h2g ) was moderate to high (0.20 to 0.74). Family-by-location interactions were significant and particularly large for water soluble carbohydrate (WSC), crude fat, phosphorus (P) and crude protein. GBLUP, KGD-GBLUP and BayesCπ genomic prediction models displayed similar predictive ability, estimated by 10-fold cross validation, for all nutritive traits with values ranging from r = 0.16 to 0.45 using phenotypes from across two locations. High predictive ability was observed for the mineral traits sulfur (0.44), sodium (0.45) and magnesium (0.45) and the lowest values were observed for P (0.16), digestibility (0.22) and high molecular weight WSC (0.23). Predictive ability estimates for most nutritive traits were retained when marker number was reduced from one million to as few as 50,000. The moderate to high predictive abilities observed suggests implementation of genomic selection is feasible for most of the nutritive traits examined.

Environmental regulation of leaf colour in red 35S:PAP1 Arabidopsis thaliana.

* High-temperature, low-light (HTLL) treatment of 35S:PAP1 Arabidopsis thaliana over-expressing the PAP1 (Production of Anthocyanin Pigment 1) gene results in reversible reduction of red colouration, suggesting the action of additional anthocyanin regulators. High-performance liquid chromatography (HPLC), liquid chromatography mass spectrometry (LCMS) and Affimetrix-based microarrays were used to measure changes in anthocyanin, flavonoids, and gene expression in response to HTLL. * HTLL treatment of control and 35S:PAP1 A. thaliana resulted in a reversible reduction in the concentrations of major anthocyanins despite ongoing over-expression of the PAP1 MYB transcription factor. Twenty-one anthocyanins including eight cis-coumaryl esters were identified by LCMS. The concentrations of nine anthocyanins were reduced and those of three were increased, consistent with a sequential process of anthocyanin degradation. Analysis of gene expression showed down-regulation of flavonol and anthocyanin biosynthesis and of transport-related genes within 24 h of HTLL treatment. No catabolic genes up-regulated by HTLL were found. * Reductions in the concentrations of anthocyanins and down-regulation of the genes of anthocyanin biosynthesis were achieved by environmental manipulation, despite ongoing over-expression of PAP1. Quantitative PCR showed reduced expression of three genes (TT8, TTG1 and EGL3) of the PAP1 transcriptional complex, and increased expression of the potential transcriptional repressors AtMYB3, AtMYB6 and AtMYBL2 coincided with HTLL-induced down-regulation of anthocyanin biosynthesis. * HTLL treatment offers a model system with which to explore anthocyanin catabolism and to discover novel genes involved in the environmental control of anthocyanins.

Semi-quantitative and structural metabolic phenotyping by direct infusion ion trap mass spectrometry and its application in genetical metabolomics.

The identification of quantitative trait loci (QTL) for plant metabolites requires the quantitation of these metabolites across a large range of progeny. We developed a rapid metabolic profiling method using both untargeted and targeted direct infusion tandem mass spectrometry (DIMSMS) with a linear ion trap mass spectrometer yielding sufficient precision and accuracy for the quantification of a large number of metabolites in a high-throughput environment. The untargeted DIMSMS method uses top-down data-dependent fragmentation yielding MS(2) and MS(3) spectra. We have developed software tools to assess the structural homogeneity of the MS(2) and MS(3) spectra hence their utility for phenotyping and genetical metabolomics. In addition we used a targeted DIMS(MS) method for rapid quantitation of specific compounds. This method was compared with targeted LC/MS/MS methods for these compounds. The DIMSMS methods showed sufficient precision and accuracy for QTL discovery. We phenotyped 200 individual Lolium perenne genotypes from a mapping population harvested in two consecutive years. Computational and statistical analyses identified 246 nominal m/z bins with sufficient precision and homogeneity for QTL discovery. Comparison of the data for specific metabolites obtained by DIMSMS with the results from targeted LC/MS/MS analysis showed that quantitation by this metabolic profiling method is reasonably accurate. Of the top 100 MS(1) bins, 22 ions gave one or more reproducible QTL across the 2 years.

A large-scale metabolomics study to harness chemical diversity and explore biochemical mechanisms in ryegrass.

Perennial ryegrass (Lolium perenne) is integral to temperate pastoral agriculture, which contributes most of the milk and meat production worldwide. Chemical profiles and diversity of ryegrass offer several opportunities to harness specific traits and elucidate underlying biological mechanisms for forage improvement. We conducted a large-scale metabolomics study of perennial ryegrass comprising 715 genotypes, representing 118 populations from 21 countries. Liquid/gas chromatography-mass spectrometry based targeted and non-targeted techniques were used to analyse fructan oligosaccharides, lipids, fatty acid methyl esters, polar and semi-polar compounds. Fructan diversity across all genotypes was evaluated, high- and low-sugar groups identified, and fructan accumulation mechanisms explored. Metabolites differentiating the two groups were characterised, modules and pathways they represent deduced, and finally, visualisation and interpretation provided in a biological context. We also demonstrate a workflow for large-scale metabolomics studies from raw data through to statistical and pathway analysis. Raw files and metadata are available at the MetaboLights database.

Phenotypic characterization of transposon-inserted mutants of Clostridium proteoclasticum B316T using extracellular metabolomics.

The development of high-throughput DNA sequencing techniques has enabled the sequencing of several hundred bacterial genomes. However, the major step towards understanding the molecular basis of an organism will be the determination of all gene functions in its genome. Current gene assignments by sequence homology generate numerous hints to putative or unknown functions. Even hits with good homology are often not specific enough to describe the in vivo biochemical functions and the underlying biological roles. In this work we applied metabolic footprinting analysis to characterize Tn916-inserted mutants of a hemicellulose-degrading rumen bacterium grown on complex culture medium. Interestingly, the most distinctive phenotypic difference was observed in a mutant with a transposon insertion in a non-coding region of the genome, while disruption of a gene with high homology to a known alpha-glucosidase/xylosidase showed no distinctive phenotypic effect. Our results demonstrate that extracellular metabolomics data coupled to genome information is a powerful and low-cost approach to rapidly screen and characterize microbial mutants with single gene deletions. However, metabolomics as a stand-alone technique is unlikely to give a complete answer to define gene functions, and, therefore, is an approach to be used to generate hypotheses and direct new experiments to confirm gene function.

E/Z-Thesinine-O-4'-alpha-rhamnoside, pyrrolizidine conjugates produced by grasses (Poaceae).

Based on direct infusion mass spectrometry we identified a novel alkaloid as a major component of perennial ryegrass (Lolium perenne). Initial mass spectral data suggested it to be a pyrrolizidine conjugate. As this class of alkaloids has not been described before from grasses, we isolated it to elucidate its structure. The isolated alkaloid proved to be a mixture of two stereoisomers. The structures of the two compounds as determined by 1D and 2D NMR spectroscopy, were E-thesinine-O-4'-alpha-rhamnoside (1) and Z-thesinine-O-4'-alpha-rhamnoside (2). These identifications were supported by the characterisation by GC-MS and optical rotation of (+)-isoretronecanol as the necine base released on alkaline hydrolysis of these alkaloids. 1 and 2 together with the aglycone and a hexoside were also detected in tall fescue (Festuca arundinacea). This is the first report of pyrrolizidine alkaloids produced by grasses (Poaceae).