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

Genetic architecture of inter-specific and -generic grass hybrids by network analysis on multi-omics data

BACKGROUND

Understanding the mechanisms underlining forage production and its biomass nutritive quality at the omics level is crucial for boosting the output of high-quality dry matter per unit of land. Despite the advent of multiple omics integration for the study of biological systems in major crops, investigations on forage species are still scarce.

RESULTS

Our results identified substantial changes in gene co-expression and metabolite-metabolite network topologies as a result of genetic perturbation by hybridizing L. perenne with another species within the genus (L. multiflorum) relative to across genera (F. pratensis). However, conserved hub genes and hub metabolomic features were detected between pedigree classes, some of which were highly heritable and displayed one or more significant edges with agronomic traits in a weighted omics-phenotype network. In spite of tagging relevant biological molecules as, for example, the light-induced rice 1 (LIR1), hub features were not necessarily better explanatory variables for omics-assisted prediction than features stochastically sampled and all available regressors.

CONCLUSIONS

The utilization of computational techniques for the reconstruction of co-expression networks facilitates the identification of key omic features that serve as central nodes and demonstrate correlation with the manifestation of observed traits. Our results also indicate a robust association between early multi-omic traits measured in a greenhouse setting and phenotypic traits evaluated under field conditions.

Analysis of chemotypes and their markers in leaves of core collections of Eucommia ulmoides using metabolomics

The leaves of Eucommia ulmoides contain various active compunds and nutritional components, and have successively been included as raw materials in the Chinese Pharmacopoeia, the Health Food Raw Material Catalogue, and the Feed Raw Material Catalogue. Core collections of E. ulmoides had been constructed from the conserved germplasm resources basing on molecular markers and morphological traits, however, the metabolite diversity and variation in this core population were little understood. Metabolite profiles of E. ulmoides leaves of 193 core collections were comprehensively characterized by GC-MS and LC-MS/MS based non-targeted metabolomics in present study. Totally 1,100 metabolites were identified and that belonged to 18 categories, and contained 120 active ingredients for traditional Chinese medicine (TCM) and 85 disease-resistant metabolites. Four leaf chemotypes of the core collections were established by integrated uses of unsupervised self-organizing map (SOM), supervised orthogonal partial least squares discriminant analysis (OPLS-DA) and random forest (RF) statistical methods, 30, 23, 43, and 23 chemomarkers were screened corresponding to the four chemotypes, respectively. The morphological markers for the chemotypes were obtained by weighted gene co-expression network analysis (WGCNA) between the chenomarkers and the morphological traits, with leaf length (LL), chlorophyll reference value (CRV), leaf dentate height (LDH), and leaf thickness (LT) corresponding to chemotypes I, II, III, and IV, respectively. Contents of quercetin-3-O-pentosidine, isoquercitrin were closely correlated to LL, leaf area (LA), and leaf perimeter (LP), suggesting the quercetin derivatives might influence the growth and development of E. ulmoides leaf shape.

Metabonomic profiling of clubroot-susceptible and clubroot-resistant radish and the assessment of disease-resistant metabolites

Plasmodiophora brassicae causes a serious threat to cruciferous plants including radish (Raphanus sativus L.). Knowledge on the pathogenic regularity and molecular mechanism of P. brassicae and radish is limited, especially on the metabolism level. In the present study, clubroot-susceptible and clubroot-resistant cultivars were inoculated with P. brassicae Race 4, root hairs initial infection of resting spores (10⁷ CFU/mL) at 24 h post-inoculation and root galls symptom arising at cortex splitting stage were identified on both cultivars. Root samples of cortex splitting stage of two cultivars were collected and used for untargeted metabonomic analysis. We demonstrated changes in metabolite regulation and pathways during the cortex splitting stage of diseased roots between clubroot-susceptible and clubroot-resistant cultivars using untargeted metabonomic analysis. We identified a larger number of differentially regulated metabolites and heavier metabolite profile changes in the susceptible cultivar than in the resistant counterpart. The metabolites that were differentially regulated in both cultivars were mostly lipids and lipid-like molecules. Significantly regulated metabolites and pathways according to the P value and variable important in projection score were identified. Moreover, four compounds, including ethyl α-D-thioglucopyranoside, imipenem, ginsenoside Rg1, and 6-gingerol, were selected, and their anti-P. brassicae ability and effects on seedling growth were verified on the susceptible cultivar. Except for ethyl α-D-thioglucopyranoside, the remaining could inhibit clubroot development of varing degree. The use of 5 mg/L ginsenoside Rg1 + 5 mg/L 6-gingerol resulted in the lowest disease incidence and disease index among all treatments and enhanced seedling growth. The regulation of pathways or metabolites of carbapenem and ginsenoside was further explored. The results provide a preliminary understanding of the interaction between radish and P. brassicae at the metabolism level, as well as the development of measures for preventing clubroot.

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.

An efficient LC-MS method for isomer separation and detection of sugars, phosphorylated sugars, and organic acids

Assessing central carbon metabolism in plants can be challenging due to the dynamic range in pool sizes, with low levels of important phosphorylated sugars relative to more abundant sugars and organic acids. Here, we report a sensitive liquid chromatography-mass spectrometry method for analysing central metabolites on a hybrid column, where both anion-exchange and hydrophilic interaction chromatography (HILIC) ligands are embedded in the stationary phase. The liquid chromatography method was developed for enhanced selectivity of 27 central metabolites in a single run with sensitivity at femtomole levels observed for most phosphorylated sugars. The method resolved phosphorylated hexose, pentose, and triose isomers that are otherwise challenging. Compared with a standard HILIC approach, these metabolites had improved peak areas using our approach due to ion enhancement or low ion suppression in the biological sample matrix. The approach was applied to investigate metabolism in high lipid-producing tobacco leaves that exhibited increased levels of acetyl-CoA, a precursor for oil biosynthesis. The application of the method to isotopologue detection and quantification was considered through evaluating 13C-labeled seeds from Camelina sativa. The method provides a means to analyse intermediates more comprehensively in central metabolism of plant tissues.

Development of RIKEN Plant Metabolome MetaDatabase

The advancement of metabolomics in terms of techniques for measuring small molecules has enabled the rapid detection and quantification of numerous cellular metabolites. Metabolomic data provide new opportunities to gain a deeper understanding of plant metabolism that can improve the health of both plants and humans that consume them. Although major public repositories for general metabolomic data have been established, the community still has shortcomings related to data sharing, especially in terms of data reanalysis, reusability and reproducibility. To address these issues, we developed the RIKEN Plant Metabolome MetaDatabase (RIKEN PMM, http://metabobank.riken.jp/pmm/db/plantMetabolomics), which stores mass spectrometry-based (e.g. gas chromatography-MS-based) metabolite profiling data of plants together with their detailed, structured experimental metadata, including sampling and experimental procedures. Our metadata are described as Linked Open Data based on the Resource Description Framework using standardized and controlled vocabularies, such as the Metabolomics Standards Initiative Ontology, which are to be integrated with various life and biomedical science data using the World Wide Web. RIKEN PMM implements intuitive and interactive operations for plant metabolome data, including raw data (netCDF format), mass spectra (NIST MSP format) and metabolite annotations. The feature is suitable not only for biologists who are interested in metabolomic phenotypes, but also for researchers who would like to investigate life science in general through plant metabolomic approaches.

DBnorm as an R package for the comparison and selection of appropriate statistical methods for batch effect correction in metabolomic studies

As a powerful phenotyping technology, metabolomics provides new opportunities in biomarker discovery through metabolome-wide association studies (MWAS) and the identification of metabolites having a regulatory effect in various biological processes. While mass spectrometry-based (MS) metabolomics assays are endowed with high throughput and sensitivity, MWAS are doomed to long-term data acquisition generating an overtime-analytical signal drift that can hinder the uncovering of real biologically relevant changes. We developed “dbnorm”, a package in the R environment, which allows for an easy comparison of the model performance of advanced statistical tools commonly used in metabolomics to remove batch effects from large metabolomics datasets. “dbnorm” integrates advanced statistical tools to inspect the dataset structure not only at the macroscopic (sample batches) scale, but also at the microscopic (metabolic features) level. To compare the model performance on data correction, “dbnorm” assigns a score that help users identify the best fitting model for each dataset. In this study, we applied “dbnorm” to two large-scale metabolomics datasets as a proof of concept. We demonstrate that “dbnorm” allows for the accurate selection of the most appropriate statistical tool to efficiently remove the overtime signal drift and to focus on the relevant biological components of complex datasets.

Predicting the quality of ryegrass using hyperspectral imaging

BACKGROUND

The quality of forage plants is a crucial component of animal performance and a limiting factor in pasture based production systems. Key forage attributes that may require improvement include the sugar, lipid, protein and energy contents of the vegetative parts of these plants. The aim of this study was to evaluate the potential capacity of hyperspectral imaging (HSI) for non-invasive assessment of forage chemical composition. Hyperspectral image data within the visible near-infrared range into the extended near-infrared covering 550-1700 nm wavelengths were obtained from 185 accessions of ryegrass (Lolium perenne), which were also analysed for 13 forage quality attributes.

RESULTS

Medium to high predictive power was observed for the HSI models of total sugars (R² validation of 0.58), high molecular weight sugars (R² validation of 0.63), %Ash (R² validation of 0.50) and %Nitrogen (R² validation of 0.70). Significant HSI models with low R² validation of 0.1-0.5 were also obtained for low molecular weight sugars, NDF (%), ADF (%), DOMD (% DM), ME (MJ/kg DM), DM (%), Ca (mg/g) and OM (%). We also observed significant differences in the chemical composition between the pseudostems and leaves of the plants for each accession. The power of HSI for prediction of these differences within plants was also demonstrated.

CONCLUSION

This study paves the way for the HSI technology to be used for in-field estimation of forage composition attributes in perennial ryegrass. This will allow more rapid genetic-based selection and breeding for a trait that is normally expensive to measure providing a cheaper, non-destructive and high throughput screening tool.