Normalization and integration of large-scale metabolomics data using support vector regression

Introduction

Untargeted metabolomics studies for biomarker discovery often have hundreds to thousands of human samples. Data acquisition of large-scale samples has to be divided into several batches and may span from months to as long as several years. The signal drift of metabolites during data acquisition (intra- and inter-batch) is unavoidable and is a major confounding factor for large-scale metabolomics studies.

Objectives

We aim to develop a data normalization method to reduce unwanted variations and integrate multiple batches in large-scale metabolomics studies prior to statistical analyses.

Methods

We developed a machine learning algorithm-based method, support vector regression (SVR), for large-scale metabolomics data normalization and integration. An R package named MetNormalizer was developed and provided for data processing using SVR normalization.

Results

After SVR normalization, the portion of metabolite ion peaks with relative standard deviations (RSDs) less than 30 % increased to more than 90 % of the total peaks, which is much better than other common normalization methods. The reduction of unwanted analytical variations helps to improve the performance of multivariate statistical analyses, both unsupervised and supervised, in terms of classification and prediction accuracy so that subtle metabolic changes in epidemiological studies can be detected.

Conclusion

SVR normalization can effectively remove the unwanted intra- and inter-batch variations, and is much better than other common normalization methods.

     

Characterization of interactions between hepatitis C virus NS5B polymerase,annexin A2 and RNA – effects on NS5B catalysis and allosteric inhibition

Background

Direct acting antivirals (DAAs) provide efficient hepatitis C virus (HCV) therapy and clearance for a majority of patients, but are not available or effective for all patients. They risk developing HCV-induced hepatocellular carcinoma (HCC), for which the mechanism remains obscure and therapy is missing. Annexin A2 (AnxA2) has been reported to co-precipitate with the non-structural (NS) HCV proteins NS5B and NS3/NS4A, indicating a role in HCC tumorigenesis and effect on DAA therapy.

Methods

Surface plasmon resonance biosensor technology was used to characterize direct interactions between AnxA2 and HCV NS5B, NS3/NS4 and RNA, and the subsequent effects on catalysis and inhibition.

Results

No direct interaction between AnxA2 and NS3/NS4A was detected, while AnxA2 formed a slowly dissociating, high affinity (K _D?=?30 nM), complex with NS5B, decreasing its catalytic activity and affinity for the allosteric inhibitor filibuvir. The RNA binding of the two proteins was independent and AnxA2 and NS5B interacted with different RNAs in ternary complexes of AnxA2:NS5B:RNA, indicating specific preferences.

Conclusions

The complex interplay revealed between NS5B, AnxA2, RNA and filibuvir, suggests that AnxA2 may have an important role for the progression and treatment of HCV infections and the development of HCC, which should be considered also when designing new allosteric inhibitors.

     

Metastatic tumor cells – genotypes and phenotypes

Background

Metastasis is the primary cause of mortality in cancer patients. Therefore, elucidating the genetics and epigenetics of metastatic tumor cells and the mechanisms by which tumor cells acquire metastatic properties constitute significant challenges in cancer research.

Objective

To summarize the current understandings of the specific genotype and phenotype of the metastatic tumor cells.

Method and Result

In-depth genetic analysis of tumor cells, especially with advances in the next-generation sequencing, have revealed insights of the genotypes of metastatic tumor cells. Also, studies have shown that the cancer stem cell (CSC) and epithelial to mesenchymal transition (EMT) phenotypes are associated with the metastatic cascade.

Conclusion

In this review, we will discuss recent advances in the field by focusing on the genomic instability and phenotypic dynamics of metastatic tumor cells.

     

Determination of fatty acid and sphingoid base composition of eleven ceramide subclasses in <Emphasis Type="Italic">stratum corneum</Emphasis> by UHPLC/scheduled-MRM

Introduction

Ceramides play a key role in skin barrier function in homeostatic and pathological conditions and can be sampled non-invasively through stratum corneum collection.

Objectives

To develop a novel UHPLC/Scheduled MRM method for the identification and relative distribution of eleven classes of ceramides, which are separated by UHPLC and determined by their specific retention times. The precise composition of the fatty acid and sphingoid base parts of each individual ceramide is determined via mass fragmentation.

Methods

More than 1000 human and pig ceramides were identified. Three human and minipig ceramide classes, CER[AS], CER[NS] and CER[EOS] have been investigated in depth.

Results

Sphingoid bases were characterized by a prevalence of chain lengths with sizes from C16 to C22, whereas fatty acids were mainly observed in the range of C22–C26. Overall, the ceramide profiles between human and minipig stratum corneum were similar. Differences in the CER[AS] and CER[NS] classes included a more homogeneous distribution of fatty acids (16–30 carbon atoms) in minipig, whereas in human longer fatty acid chains (>?24 carbon atoms) predominated.

Conclusion

The method will be useful for the analysis of healthy and pathological skin in various specie, and the measurement of the relative distribution of ceramides as biomarkers for pharmacodynamic studies.

     

Targeted/exome sequencing identified mutations in ten Chinese patients diagnosed with Noonan syndrome and related disorders

Background

Noonan syndrome (NS) and Noonan syndrome with multiple lentigines (NSML) are autosomal dominant developmental disorders. NS and NSML are caused by abnormalities in genes that encode proteins related to the RAS-MAPK pathway, including PTPN11, RAF1, BRAF, and MAP2K. In this study, we diagnosed ten NS or NSML patients via targeted sequencing or whole exome sequencing (TS/WES).

Methods

TS/WES was performed to identify mutations in ten Chinese patients who exhibited the following manifestations: potential facial dysmorphisms, short stature, congenital heart defects, and developmental delay. Sanger sequencing was used to confirm the suspected pathological variants in the patients and their family members.

Results

TS/WES revealed three mutations in the PTPN11 gene, three mutations in RAF1 gene, and four mutations in BRAF gene in the NS and NSML patients who were previously diagnosed based on the abovementioned clinical features. All the identified mutations were determined to be de novo mutations. However, two patients who carried the same mutation in the RAF1 gene presented different clinical features. One patient with multiple lentigines was diagnosed with NSML, while the other patient without lentigines was diagnosed with NS. In addition, a patient who carried a hotspot mutation in the BRAF gene was diagnosed with NS instead of cardiofaciocutaneous syndrome (CFCS).

Conclusions

TS/WES has emerged as a useful tool for definitive diagnosis and accurate genetic counseling of atypical cases. In this study, we analyzed ten Chinese patients diagnosed with NS and related disorders and identified their correspondingPTPN11, RAF1, and BRAF mutations. Among the target genes, BRAF showed the same degree of correlation with NS incidence as that of PTPN11 or RAF1.

     

Hypocapnia is a biological marker for orthostatic intolerance in some patients with chronic fatigue syndrome

Context

Patients with chronic fatigue syndrome and those with orthostatic intolerance share many symptoms, yet questions exist as to whether CFS patients have physiological evidence of orthostatic intolerance.

Objective

To determine if some CFS patients have increased rates of orthostatic hypotension, hypertension, tachycardia, or hypocapnia relative to age-matched controls.

Design

Assess blood pressure, heart rate, respiratory rate, end tidal CO2 and visual analog scales for orthostatic symptoms when supine and when standing for 8 minutes without moving legs.

Setting

Referral practice and research center.

Participants

60 women and 15 men with CFS and 36 women and 4 men serving as age matched controls with analyses confined to 62 patients and 35 controls showing either normal orthostatic testing or a physiological abnormal test.

Main outcome measures

Orthostatic tachycardia; orthostatic hypotension; orthostatic hypertension; orthostatic hypocapnia or combinations thereof.

Results

CFS patients had higher rates of abnormal tests than controls (53% vs 20%, p < .002), but rates of orthostatic tachycardia, orthostatic hypotension, and orthostatic hypertension did not differ significantly between patients and controls (11.3% vs 5.7%, 6.5% vs 2.9%, 19.4% vs 11.4%, respectively). In contrast, rates of orthostatic hypocapnia were significantly higher in CFS than in controls (20.6% vs 2.9%, p < .02). This CFS group reported significantly more feelings of illness and shortness of breath than either controls or CFS patients with normal physiological tests.

Conclusion

A substantial number of CFS patients have orthostatic intolerance in the form of orthostatic hypocapnia. This allows subgrouping of patients with CFS and thus reduces patient pool heterogeneity engendered by use of a clinical case definition.

     

Screening for internal phosphorus utilisation efficiency: comparison of genotypes at equal shoot P content is critical

Aims

Progress in improving the internal phosphorus utilisation efficiency of crops has been limited, which may be due to poor screening methods that allow differences in P uptake among genotypes grown in soil to mask genotypic differences in shoot biomass produced per unit of shoot P (PUE). We investigated alternative soil and hydroponic screening methods for their capacity to produce a consensus ranking of genotypes with regard to PUE.

Methods

Six rice genotypes previously identified in hydroponic screening studies as being high, intermediate or low in PUE were screened using multi P rate hydroponic and soil-based experiments.

Results

Comparisons made at each rate of soil-P supply produced estimates of PUE strongly biased by P uptake differences among genotypes. Using multiple-rate data to derive response functions per genotype showed that similar P content was achieved at different rates of P supply but that high-PUE genotypes clearly separated from intermediate- and low-PUE genotypes if equal P content was used. Ranking analysis suggested that results obtained from soil agreed well with those from the hydroponic study.

Conclusions

PUE was significantly influenced by genotype and P supply, but there was no significant genotype x P supply interaction. Hence, we conclude that screening genotypes using hydroponics at one or two P supply levels is the most cost- and time effective means to screen large numbers of rice genotypes for PUE.

     

Metabolic response of porcine colon explants to in vitro infection by <Emphasis Type="Italic">Brachyspira hyodysenteriae</Emphasis>: a leap into disease pathophysiology

Introduction

Swine dysentery caused by Brachyspira hyodysenteriae is a production limiting disease in pig farming. Currently antimicrobial therapy is the only treatment and control method available.

Objective

The aim of this study was to characterize the metabolic response of porcine colon explants to infection by B. hyodysenteriae.

Methods

Porcine colon explants exposed to B. hyodysenteriae were analyzed for histopathological, metabolic and pro-inflammatory gene expression changes.

Results

Significant epithelial necrosis, increased levels of l-citrulline and IL-1α were observed on explants infected with B. hyodysenteriae.

Conclusions

The spirochete induces necrosis in vitro likely through an inflammatory process mediated by IL-1α and NO.

     

Optimization of fecal sample preparation for untargeted LC-HRMS based metabolomics

Introduction

Collecting feces is easy. It offers direct outcome to endogenous and microbial metabolites.

Objectives

In a context of lack of consensus about fecal sample preparation, especially in animal species, we developed a robust protocol allowing untargeted LC-HRMS fingerprinting.

Methods

The conditions of extraction (quantity, preparation, solvents, dilutions) were investigated in bovine feces.

Results

A rapid and simple protocol involving feces extraction with methanol (1/3, M/V) followed by centrifugation and a step filtration (10 kDa) was developed.

Conclusion

The workflow generated repeatable and informative fingerprints for robust metabolome characterization.

     

A lost opportunity for science: journals promote data sharing in metabolomics but do not enforce it

Introduction

Data sharing is being increasingly required by journals and has been heralded as a solution to the ‘replication crisis’.

Objectives

(i) Review data sharing policies of journals publishing the most metabolomics papers associated with open data and (ii) compare these journals’ policies to those that publish the most metabolomics papers.

Methods

A PubMed search was used to identify metabolomics papers. Metabolomics data repositories were manually searched for linked publications.

Results

Journals that support data sharing are not necessarily those with the most papers associated to open metabolomics data.

Conclusion

Further efforts are required to improve data sharing in metabolomics.

     

Genetic variation in the TAS2R38 taste receptor contributes to the oral microbiota in North and South European locations: a pilot study

Background

Microbial communities are influenced by environmental factors including host genetics. We investigated the relationship between host bitter taste receptor genotype hTAS2R38 and oral microbiota, together with the influence of geographical location.

Methods

hTAS2R38 polymorphisms and 16S bacterial gene sequencing from oral samples were analyzed from a total of 45 healthy volunteers from different geographical locations.

Results

Genetic variation in the bitter taste receptor TAS2R38 reflected in the microbial composition of oral mucosa in Finnish and Spanish subjects. Multivariate analysis showed significant differences in the microbial composition between country and also dependent on taste genotype. Oral microbiota was shown to be more stable to the geographical location impact among AVI-homozygotes than PAV-homozygotes or heterozygotes (PAV/AVI).

Conclusion

Geographical location and genetic variation in the hTAS2R38 taste receptor impact oral mucosa microbial composition. These findings provide an advance in the knowledge regarding the interactions between taste receptor genes and oral microbiota. This study suggests the role of host-microbiota interactions on the food taste perception in food choices, nutrition, and eating behavior.

     

A methodology for elucidating regulatory mechanisms leading to changes in lipid profiles

Introduction

It is difficult to elucidate the metabolic and regulatory factors causing lipidome perturbations.

Objectives

This work simplifies this process.

Methods

A method has been developed to query an online holistic lipid metabolic network (of 7923 metabolites) to extract the pathways that connect the input list of lipids.

Results

The output enables pathway visualisation and the querying of other databases to identify potential regulators. When used to a study a plasma lipidome dataset of polycystic ovary syndrome, 14 enzymes were identified, of which 3 are linked to ELAVL1—an mRNA stabiliser.

Conclusion

This method provides a simplified approach to identifying potential regulators causing lipid-profile perturbations.

     

Targeted and global pharmacometabolomics in everolimus-based immunosuppression: association of co-medication and lysophosphatidylcholines with dose requirement

Introduction

The immunosuppressive therapy with everolimus (ERL) after heart transplantation is characterized by a narrow therapeutic window and a substantial variability in dose requirement. Factors explaining this variability are largely unknown.

Objectives

Our aim was to evaluate factors affecting ERL metabolism and to identify novel metabolites associated with the individual ERL dose requirement to elucidate mechanisms underlying ERL dose response variability.

Method

We used liquid chromatography coupled with mass spectrometry for quantification of ERL metabolites in 41 heart transplant patients and evaluated the effect of clinical and genetic factors on ERL pharmacokinetics. Non-targeted plasma metabolic profiling by ultra-performance liquid chromatography and high resolution quadrupole-time-of-flight mass spectrometry was used to identify novel metabolites associated with ERL dose requirement.

Results

The determination of ERL metabolites revealed differences in metabolite patterns that were independent from clinical or genetic factors. Whereas higher ERL dose requirement was associated with co-administration of sodium-mycophenolic acid and the CYP3A5 expressor genotype, lower dose was required for patients receiving vitamin K antagonists. Global metabolic profiling revealed several novel metabolites associated with ERL dose requirement. One of them was identified as lysophosphatidylcholine (lysoPC) (16:0/0:0). Subsequent targeted analysis revealed that high levels of several lysoPCs were significantly associated with higher ERL dose requirement.

Conclusion

For the first time, this study describes distinct ERL metabolite patterns in heart transplant patients and detected potentially new drug–drug interactions. The global metabolic profiling facilitated the discovery of novel metabolites associated with ERL dose requirement that might represent new clinically valuable biomarkers to guide ERL therapy.

     

The NS1 protein of influenza a virus interacts with heat shock protein Hsp90 in human alveolar basal epithelial cells: Implication for virus-induced apoptosis

Background

Our previous study showed that the NS1 protein of highly pathogenic avian influenza A virus H5N1 induced caspase-dependent apoptosis in human alveolar basal epithelial cells (A549), supporting its function as a proapoptotic factor during viral infection, but the mechanism is still unknown.

Results

To characterize the mechanism of NS1-induced apoptosis, we used a two-hybrid system to isolate the potential NS1-interacting partners in A549 cells. We found that heat shock protein 90 (Hsp90) was able to interact with the NS1 proteins derived from both H5N1 and H3N2 viruses, which was verified by co-immunoprecitation assays. Significantly, the NS1 expression in the A549 cells dramatically weakened the interaction between Apaf-1 and Hsp90 but enhanced its interaction with cytochrome c (Cyt c), suggesting that the competitive binding of NS1 to Hsp90 might promote the Apaf-1 to associate with Cyt c and thus facilitate the activation of caspase 9 and caspase 3.

Conclusions

The present results demonstrate that NS1 protein of Influenza A Virus interacts with heat hock protein Hsp90 and meidates the apoptosis induced by influenza A virus through the caspase cascade.

     

Hyper response to ovarian stimulation affects the follicular fluid metabolomic profile of women undergoing IVF similarly to polycystic ovary syndrome

Introduction

During in vitro fertilization (IVF), the hyper response to controlled ovarian stimulation (COS) is a common characteristic among patients diagnosed with polycystic ovary syndrome (PCOS), although non-diagnosed patients may also demonstrate this response.

Objectives

In an effort to investigate follicular metabolic characteristics associated with hyper response to COS, the present study analyzed follicular fluid (FF) samples from patients undergoing IVF.

Methods

FF samples were obtained from patients with PCOS and hyper response during IVF (PCOS group, N?=?15), patients without PCOS but with hyper response during IVF (HR group, N?=?44), and normo-responder patients receiving IVF (control group, N?=?22). FF samples underwent Bligh and Dyer extraction, followed by metabolomic analysis by ultra-performance liquid chromatography mass spectrometry, considering two technical replicates. Clinical data was analyzed by ANOVA and chi-square tests. The metabolomic dataset was analyzed by multivariate statistics, and the significance of biomarkers was confirmed by ANOVA.

Results

Clinical data showed differences regarding follicles production, oocyte and embryo quality. From the 15 proposed biomarkers, 14 were of increased abundance in the control group and attributed as fatty acids, diacylglycerol, triacylglycerol, ceramide, ceramide-phosphate, phosphatidylcholine, and sphingomyelin. The PCOS patients showed increased abundance of a metabolite of m/z 144.0023 that was not attributed to a class.

Conclusion

The clinical and metabolic similarities observed in the FF of hyper responders with and without PCOS diagnosis indicate common biomarkers that could assist on the development of accessory tools for assessment of IVF parameters.

     

Characterization of ankylosing spondylitis and rheumatoid arthritis using <Superscript>1</Superscript>H NMR-based metabolomics of human fecal extracts

Introduction

The differences in fecal metabolome between ankylosing spondylitis (AS)/rheumatoid arthritis (RA) patients and healthy individuals could be the reason for an autoimmune disorder.

Objectives

The study explored the fecal metabolome difference between AS/RA patients and healthy controls to clarify human immune disturbance.

Methods

Fecal samples from 109 individuals (healthy controls 34, AS 40, and RA 35) were analyzed by ¹H NMR spectroscopy. Data were analyzed with principal component analysis (PCA) and orthogonal projection to latent structure discriminant (OPLS-DA) analysis.

Results

Significant differences in the fecal metabolic profiles could distinguish AS/RA patients from healthy controls but could not distinguish between AS and RA patients. The significantly decreased metabolites in AS/RA patients were butyrate, propionate, methionine, and hypoxanthine. Significantly increased metabolites in AS/RA patients were taurine, methanol, fumarate, and tryptophan.

Conclusion

The metabolome variations in feces indicated AS and RA were two homologous diseases that could not be distinguished by ¹H NMR metabolomics.

     

Factor XIII Val34Leu polymorphism and risk of recurrent pregnancy loss in Iranian population: a case control study

BACKGROUND

Recurrent pregnancy loss (RPL) is a heterogeneous condition and thrombophilias have been considered as a probable cause.

OBJECTIVE

The aim of this study was to investigate the prevalence of the coagulation factor XIII Val34Leu polymorphism among women with unexplained RPL.

METHODS

A total of 140 women with a history of unexplained RPL and 100 age-matched healthy fertile women were recruited. The presence of FXIII Val34Leu polymorphism among the cases and controls was investigated using PCR-RFLP method.

RESULTS

Genotype analyses of the subjects revealed that the patients had a significantly higher prevalence of V/L and L/L than the controls (P<0.05): 33.5% vs. 15%, and 9.2% vs. 2%, respectively.

CONCLUSION

These results indicate a significant association between FXIII Val34Leu polymorphism and unexplained RPL in the Iranian patient.

     

Corneal Collagen Cross-Linking for the Management of Mycotic Keratitis

Purpose

To evaluate the efficiency of corneal collagen cross-linking (CXL) in addition to topical voriconazole in cases with mycotic keratitis.

Design

Retrospective case series in a tertiary university hospital.

Participants

CXL was performed on 13 patients with mycotic keratitis who presented poor or no response to topical voriconazole treatment.

Methods

The clinical features, symptoms, treatment results and complications were recorded retrospectively. The corneal infection was graded according to the depth of infection into the stroma (from grade 1 to grade 3). The visual analogue scale was used to calculate the pain score before and 2 days after surgery.

Main Outcome Measures

Grade of the corneal infection.

Results

Mean age of 13 patients (6 female and 7 male) was 42.4 ± 17.7 years (20–74 years). Fungus was demonstrated in culture (eight patients) or cytological examination (five patients). Seven of the 13 patients (54%) were healed with topical voriconazole and CXL adjuvant treatment in 26 ± 10 days (15–40 days). The remaining six patients did not respond to CXL treatment; they initially presented with higher grade ulcers. Pre- and post-operative pain score values were 8 ± 0.8 and 3.5 ± 1, respectively (p < 0.05).

Conclusions

The current study suggests that adjunctive CXL treatment is effective in patients with small and superficial mycotic ulcers. These observations require further research by large randomized clinical trials.

     

Metabolomics approach for predicting response to neoadjuvant chemotherapy for colorectal cancer

Introduction

Colorectal cancer (CRC) is a clinically heterogeneous disease, which necessitates a variety of treatments and leads to different outcomes. Only some CRC patients will benefit from neoadjuvant chemotherapy (NACT).

Objectives

An accurate prediction of response to NACT in CRC patients would greatly facilitate optimal personalized management, which could improve their long-term survival and clinical outcomes.

Methods

In this study, plasma metabolite profiling was performed to identify potential biomarker candidates that can predict response to NACT for CRC. Metabolic profiles of plasma from non-response (n?=?30) and response (n?=?27) patients to NACT were studied using UHPLC–quadruple time-of-flight)/mass spectrometry analyses and statistical analysis methods.

Results

The concentrations of nine metabolites were significantly different when comparing response to NACT. The area under the receiver operating characteristic curve value of the potential biomarkers was up to 0.83 discriminating the non-response and response group to NACT, superior to the clinical parameters (carcinoembryonic antigen and carbohydrate antigen 199).

Conclusion

These results show promise for larger studies that could result in more personalized treatment protocols for CRC patients.

     

Identification of a urine metabolite constellation characteristic for kidney allograft rejection

Introduction

Allograft rejection is still an important complication after kidney transplantation. Currently, monitoring of these patients mostly relies on the measurement of serum creatinine and clinical evaluation. The gold standard for diagnosing allograft rejection, i.e. performing a renal biopsy is invasive and expensive. So far no adequate biomarkers are available for routine use.

Objectives

We aimed to develop a urine metabolite constellation that is characteristic for acute renal allograft rejection.

Methods

NMR-Spectroscopy was applied to a training cohort of transplant recipients with and without acute rejection.

Results

We obtained a metabolite constellation of four metabolites that shows promising performance to detect renal allograft rejection in the cohorts used (AUC of 0.72 and 0.74, respectively).

Conclusion

A metabolite constellation was defined with the potential for further development of an in-vitro diagnostic test that can support physicians in their clinical assessment of a kidney transplant patient.