Recent advances in -omics technologies such as transcriptomics, metabolomics, and proteomics along with genotypic profiling have permitted dissection of the genetics of complex traits represented by molecular phenotypes in nonmodel species. To identify the genetic factors underlying variation in primary metabolism in potato (Solanum tuberosum), we have profiled primary metabolite content in a diploid potato mapping population, derived from crosses between S. tuberosum and wild relatives, using gas chromatography-time of flight-mass spectrometry. In total, 139 polar metabolites were detected, of which we identified metabolite quantitative trait loci for approximately 72% of the detected compounds. In order to obtain an insight into the relationships between metabolic traits and classical phenotypic traits, we also analyzed statistical associations between them. The combined analysis of genetic information through quantitative trait locus coincidence and the application of statistical learning methods provide information on putative indicators associated with the alterations in metabolic networks that affect complex phenotypic traits.

Background: Biomarker identification is one of the major and important goal of functional genomics and translational medicine studies. Large scale –omics data are increasingly being accumulated and can provide vital means for the identification of biomarkers for the early diagnosis of complex disease and/or for advanced patient/diseases stratification. These tasks are clearly interlinked, and it is essential that an unbiased and stable methodology is applied in order to address them. Although, recently, many, primarily machine learning based, biomarker identification approaches have been developed, the exploration of potential associations between biomarker identification and the design of future experiments remains a challenge.Methods: In this study, using both simulated and published experimentally derived datasets, we assessed the performance of several state-of-the-art Random Forest (RF) based decision approaches, namely the Boruta method, the permutation based feature selection without correction method, the permutation based feature selection with correction method, and the backward elimination based feature selection method. Moreover, we conducted a power analysis to estimate the number of samples required for potential future studies.Results: We present a number of different RF based stable feature selection methods and compare their performances using simulated, as well as published, experimentally derived, datasets. Across all of the scenarios considered, we found the Boruta method to be the most stable methodology, whilst the Permutation (Raw) approach offered the largest number of relevant features, when allowed to stabilise over a number of iterations. Finally, we developed and made available a web interface (https://joelarkman.shinyapps.io/PowerTools/) to streamline power calculations thereby aiding the design of potential future studies within a translational medicine context.Conclusions: We developed a RF-based biomarker discovery framework and provide a web interface for our framework, termed PowerTools, that caters the design of appropriate and cost-effective subsequent future omics study.

Although a majority of patients with PSC have colitis [PSC-IBD; primary sclerosing cholangitis-inflammatory bowel disease], this is phenotypically different from ulcerative colitis [UC]. We sought to define further the pathophysiological differences between PSC-IBD and UC, by applying a comparative and integrative approach to colonic gene expression, gut microbiota and immune infiltration data. Colonic biopsies were collected from patients with PSC-IBD [n = 10], UC [n = 10], and healthy controls [HC; n = 10]. Shotgun RNA-sequencing for differentially expressed colonic mucosal genes [DEGs], 16S rRNA analysis for microbial profiling, and immunophenotyping were performed followed by multi-omic integration. The colonic transcriptome differed significantly between groups [p = 0.01]. Colonic transcriptomes from HC were different from both UC [1343 DEGs] and PSC-IBD [4312 DEGs]. Of these genes, only 939 had shared differential gene expression in both UC and PSC-IBD compared with HC. Imputed pathways were predominantly associated with upregulation of immune response and microbial defense in both disease cohorts compared with HC. There were 1692 DEGs between PSC-IBD and UC. Bile acid signalling pathways were upregulated in PSC-IBD compared with UC [p = 0.02]. Microbiota profiles were different between the three groups [p = 0.01]; with inferred function in PSC-IBD also being consistent with dysregulation of bile acid metabolism. Th17 cells and IL17-producing CD4 cells were increased in both PSC-IBD and UC when compared with HC [p < 0.05]. Multi-omic integration revealed networks involved in bile acid homeostasis and cancer regulation in PSC-IBD. Colonic transcriptomic and microbiota analysis in PSC-IBD point toward dysregulation of colonic bile acid homeostasis compared with UC. This highlights important mechanisms and suggests the possibility of novel approaches in treating PSC-IBD.

In order to find genetic and metabolic pathways related to phenotypic traits of interest, we analyzed gene expression data, metabolite data obtained with GC-MS and LC-MS, proteomics data and a selected set of tuber quality phenotypic data from a diploid segregating mapping population of potato. In this study we present an approach to integrate these ~ omics data sets for the purpose of predicting phenotypic traits. This gives us networks of relatively small sets of interrelated ~ omics variables that can predict, with higher accuracy, a quality trait of interest. We used Random Forest regression for integrating multiple ~ omics data for prediction of four quality traits of potato: tuber flesh colour, DSC onset, tuber shape and enzymatic discoloration. For tuber flesh colour beta-carotene hydroxylase and zeaxanthin epoxidase were ranked first and forty-fourth respectively both of which have previously been associated with flesh colour in potato tubers. Combining all the significant genes, LC-peaks, GC-peaks and proteins, the variation explained was 75 %, only slightly more than what gene expression or LC-MS data explain by themselves which indicates that there are correlations among the variables across data sets. For tuber shape regressed on the gene expression, LC-MS, GC-MS and proteomics data sets separately, only gene expression data was found to explain significant variation. For DSC onset, we found 12 significant gene expression, 5 metabolite levels (GC) and 2 proteins that are associated with the trait. Using those 19 significant variables, the variation explained was 45 %. Expression QTL (eQTL) analyses showed many associations with genomic regions in chromosome 2 with also the highest explained variation compared to other chromosomes. Transcriptomics and metabolomics analysis on enzymatic discoloration after 5 min resulted in 420 significant genes and 8 significant LC metabolites, among which two were putatively identified as caffeoylquinic acid methyl ester and tyrosine. In this study, we made a strategy for selecting and integrating multiple ~ omics data using random forest method and selected representative individual peaks for networks based on eQTL, mQTL or pQTL information. Network analysis was done to interpret how a particular trait is associated with gene expression, metabolite and protein data.

Diet is a major contributor to metabolic disease risk, but there is controversy as to whether increased incidences of diseases such as non-alcoholic fatty liver disease arise from consumption of saturated fats or free sugars. Here, we investigate whether a sub-set of triacylglycerols (TAGs) were associated with hepatic steatosis and whether they arise from de novo lipogenesis (DNL) from the consumption of carbohydrates. We conduct direct infusion mass spectrometry of lipids in plasma to study the association between specific TAGs and hepatic steatosis assessed by ultrasound and fatty liver index in volunteers from the UK-based Fenland Study and evaluate clustering of TAGs in the National Survey of Health and Development UK cohort. We find that TAGs containing saturated and monounsaturated fatty acids with 16–18 carbons are specifically associated with hepatic steatosis. These TAGs are additionally associated with higher consumption of carbohydrate and saturated fat, hepatic steatosis, and variations in the gene for protein phosphatase 1, regulatory subunit 3b (PPP1R3B), which in part regulates glycogen synthesis. DNL is measured in hyperphagic ob/ob mice, mice on a western diet (high in fat and free sugar) and in healthy humans using stable isotope techniques following high carbohydrate meals, demonstrating the rate of DNL correlates with increased synthesis of this cluster of TAGs. Furthermore, these TAGs are increased in plasma from patients with biopsy-confirmed steatosis. A subset of TAGs is associated with hepatic steatosis, even when correcting for common confounding factors. We suggest that hepatic steatosis risk in western populations is in part driven by increased DNL following carbohydrate rich meals in addition to the consumption of saturated fat.

Background: Rapid non-destructive measurements to predict cassava root yield over the full growing season through large numbers of germplasm and multiple environments is a huge challenge in Cassava breeding programs. As opposed to waiting until the harvest season, multispectral imagery using unmanned aerial vehicles (UAV) are capable of measuring the canopy metrics and vegetation indices (VIs) traits at different time points of the growth cycle. This resourceful time series aerial image processing with appropriate analytical framework is very important for the automatic extraction of phenotypic features from the image data. Many studies have demonstrated the usefulness of advanced remote sensing technologies coupled with machine learning (ML) approaches for accurate prediction of valuable crop traits. Until now, Cassava has received little to no attention in aerial image-based phenotyping and ML model testing.Results: To accelerate image processing, an automated image-analysis framework called CIAT Pheno-i was developed to extract plot level vegetation indices/canopy metrics. Multiple linear regression models were constructed at different key growth stages of cassava, using ground-truth data and vegetation indices obtained from a multispectral sensor. Henceforth, the spectral indices/features were combined to develop models and predict cassava root yield using different Machine learning techniques. Our results showed that (1) Developed CIAT pheno-i image analysis framework was found to be easier and more rapid than manual methods. (2) The correlation analysis of four phenological stages of cassava revealed that elongation (EL) and late bulking (LBK) were the most useful stages to estimate above-ground biomass (AGB), below-ground biomass (BGB) and canopy height (CH). (3) The multi-temporal analysis revealed that cumulative image feature information of EL + early bulky (EBK) stages showed a higher significant correlation (r = 0.77) for Green Normalized Difference Vegetation indices (GNDVI) with BGB than individual time points. Canopy height measured on the ground correlated well with UAV (CHuav)-based measurements (r = 0.92) at late bulking (LBK) stage. Among different image features, normalized difference red edge index (NDRE) data were found to be consistently highly correlated (r = 0.65 to 0.84) with AGB at LBK stage. (4) Among the four ML algorithms used in this study, k-Nearest Neighbours (kNN), Random Forest (RF) and Support Vector Machine (SVM) showed the best performance for root yield prediction with the highest accuracy of R² = 0.67, 0.66 and 0.64, respectively.Conclusion: UAV platforms, time series image acquisition, automated image analytical framework (CIAT Pheno-i), and key vegetation indices (VIs) to estimate phenotyping traits and root yield described in this work have great potential for use as a selection tool in the modern cassava breeding programs around the world to accelerate germplasm and varietal selection. The image analysis software (CIAT Pheno-i) developed from this study can be widely applicable to any other crop to extract phenotypic information rapidly.