Downregulation of RdDM during strawberry fruit ripening

Background

Recently, DNA methylation was proposed to regulate fleshy fruit ripening. Fleshy fruits can be distinguished by their ripening process as climacteric fruits, such as tomatoes, or non-climacteric fruits, such as strawberries. Tomatoes undergo a global decrease in DNA methylation during ripening, due to increased expression of a DNA demethylase gene. The dynamics and biological relevance of DNA methylation during the ripening of non-climacteric fruits are unknown.

Results

Here, we generate single-base resolution maps of the DNA methylome in immature and ripe strawberry. We observe an overall loss of DNA methylation during strawberry fruit ripening. Thus, ripening-induced DNA hypomethylation occurs not only in climacteric fruit, but also in non-climacteric fruit. Application of a DNA methylation inhibitor causes an early ripening phenotype, suggesting that DNA hypomethylation is important for strawberry fruit ripening. The mechanisms underlying DNA hypomethylation during the ripening of tomato and strawberry are distinct. Unlike in tomatoes, DNA demethylase genes are not upregulated during the ripening of strawberries. Instead, genes involved in RNA-directed DNA methylation are downregulated during strawberry ripening. Further, ripening-induced DNA hypomethylation is associated with decreased siRNA levels, consistent with reduced RdDM activity. Therefore, we propose that a downregulation of RdDM contributes to DNA hypomethylation during strawberry ripening.

Conclusions

Our findings provide new insight into the DNA methylation dynamics during the ripening of non-climacteric fruit and suggest a novel function of RdDM in regulating an important process in plant development.

     

Associations of sex hormone-binding globulin and testosterone with genome-wide DNA methylation

Background

Levels of sex hormone-binding globulin (SHBG) and the androgen testosterone have been associated with risk of diseases throughout the lifecourse. Although both SHBG and testosterone have been shown to be highly heritable, only a fraction of that heritability has been explained by genetic studies. Epigenetic modifications such as DNA methylation may explain some of the missing heritability and could potentially inform biological knowledge of endocrine disease mechanisms involved in development of later life disease. Using data from the Avon Longitudinal Study of Parents and Children (ALSPAC), we explored cross-sectional associations of SHBG, total testosterone and bioavailable testosterone in childhood (males only) and adolescence (both males and females) with genome-wide DNA methylation. We also report associations of a SHBG polymorphism (rs12150660) with DNA methylation, which leads to differential levels of SHBG in carriers, as a genetic proxy of circulating SHBG levels.

Results

We identified several novel sites and genomic regions where levels of SHBG, total testosterone, and bioavailable testosterone were associated with DNA methylation, including one region associated with total testosterone in males (annotated to the KLHL31 gene) in both childhood and adolescence and a second region associated with bioavailable testosterone (annotated to the CMYA5 gene) at both time-points. We also identified one region where both SHBG and bioavailable testosterone in males in childhood (annotated to the ZNF718 gene) was associated with DNA methylation.

Conclusion

Our findings have important implications in the understanding of the biological processes of SHBG and testosterone, with the potential for future work to determine the molecular mechanisms that could underpin these associations.

     

A hybrid parameter estimation algorithm for beta mixtures and applications to methylation state classification

Background

Mixtures of beta distributions are a flexible tool for modeling data with values on the unit interval, such as methylation levels. However, maximum likelihood parameter estimation with beta distributions suffers from problems because of singularities in the log-likelihood function if some observations take the values 0 or 1.

Methods

While ad-hoc corrections have been proposed to mitigate this problem, we propose a different approach to parameter estimation for beta mixtures where such problems do not arise in the first place. Our algorithm combines latent variables with the method of moments instead of maximum likelihood, which has computational advantages over the popular EM algorithm.

Results

As an application, we demonstrate that methylation state classification is more accurate when using adaptive thresholds from beta mixtures than non-adaptive thresholds on observed methylation levels. We also demonstrate that we can accurately infer the number of mixture components.

Conclusions

The hybrid algorithm between likelihood-based component un-mixing and moment-based parameter estimation is a robust and efficient method for beta mixture estimation. We provide an implementation of the method (“betamix”) as open source software under the MIT license.

     

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.

     

Sex-related factors influence expression of mood-related genes in the basolateral amygdala differentially depending on age and stress exposure

Background

Women are twice as likely to be diagnosed with major depressive disorder (MDD) compared to men, but the molecular mechanisms underlying this sex difference are unclear. Previous studies in the human postmortem brain suggest dysfunction in basolateral amygdala (BLA) inhibitory gamma-aminobutyric acid (GABA) signaling and brain-derived neurotrophic factor (BDNF) function, specifically in females with MDD.

Methods

We investigated the effects of sex chromosome complement, developmental gonadal sex, and circulating testosterone on expression of 3 GABA-related and 2 BDNF-related genes in the BLA using three cohorts of four core genotypes (FCG) mice. Cohort 1 included gonadally intact pre-pubertal FCG mice; results were analyzed using two-way ANOVA (sex chromosome complement-by-gonadal sex). We examined the same genes under adult non-stressed (cohort 2) and chronically stressed conditions (cohort 3). The results for cohorts 2 and 3 were analyzed by three-way ANOVA (sex chromosome complement-by-gonadal sex-by-hormone). The use of heatmaps and Spearman correlation of BLA gene expression and anxiety-like behavior provides a global interpretation of gene expression patterns.

Results

In weanlings, we found an effect of sex chromosome complement, with lower expression of GABA/BDNF-related genes in XY mice. Most of these effects did not persist into adulthood, although a number of interesting interactions between organizational and activational effects of hormones emerged. In our adult cohorts, we found that testosterone had different effects depending on stress conditions and/or gonadal sex. Notably, in our chronically stressed adults, we found that the BLA pattern of gene expression for the GABA-related gene, somatostatin (Sst), matched the anxiety-like behavior pattern (i.e., lower Sst and higher anxiety-like behavior in XY mice, while testosterone increased Sst and decreased anxiety-like behavior). Additionally, increased Sst gene expression was correlated with decreased anxiety-like behavior.

Conclusions

Sex chromosome complement is an important factor modulating expression of mood-related genes during pre-pubertal development. The observed sex differences under chronically stressed conditions suggest that different molecular profiles may characterize male and female MDD. Our findings here for Sst are especially interesting, and suggest an underlying XY vulnerability that is typically compensated for by circulating testosterone in “normal” males. Without testosterone, women may have lower SST expression in the amygdala, resulting in increased MDD vulnerability.

     

Integrative pathway-based survival prediction utilizing the interaction between gene expression and DNA methylation in breast cancer

Background

Integrative analysis on multi-omics data has gained much attention recently. To investigate the interactive effect of gene expression and DNA methylation on cancer, we propose a directed random walk-based approach on an integrated gene-gene graph that is guided by pathway information.

Methods

Our approach first extracts a single pathway profile matrix out of the gene expression and DNA methylation data by performing the random walk over the integrated graph. We then apply a denoising autoencoder to the pathway profile to further identify important pathway features and genes. The extracted features are validated in the survival prediction task for breast cancer patients.

Results

The results show that the proposed method substantially improves the survival prediction performance compared to that of other pathway-based prediction methods, revealing that the combined effect of gene expression and methylation data is well reflected in the integrated gene-gene graph combined with pathway information. Furthermore, we show that our joint analysis on the methylation features and gene expression profile identifies cancer-specific pathways with genes related to breast cancer.

Conclusions

In this study, we proposed a DRW-based method on an integrated gene-gene graph with expression and methylation profiles in order to utilize the interactions between them. The results showed that the constructed integrated gene-gene graph can successfully reflect the combined effect of methylation features on gene expression profiles. We also found that the selected features by DA can effectively extract topologically important pathways and genes specifically related to breast cancer.

     

A streamlined method for analysing genome-wide DNA methylation patterns from low amounts of FFPE DNA

Background

Formalin fixed paraffin embedded (FFPE) tumor samples are a major source of DNA from patients in cancer research. However, FFPE is a challenging material to work with due to macromolecular fragmentation and nucleic acid crosslinking. FFPE tissue particularly possesses challenges for methylation analysis and for preparing sequencing-based libraries relying on bisulfite conversion. Successful bisulfite conversion is a key requirement for sequencing-based methylation analysis.

Methods

Here we describe a complete and streamlined workflow for preparing next generation sequencing libraries for methylation analysis from FFPE tissues. This includes, counting cells from FFPE blocks and extracting DNA from FFPE slides, testing bisulfite conversion efficiency with a polymerase chain reaction (PCR) based test, preparing reduced representation bisulfite sequencing libraries and massively parallel sequencing.

Results

The main features and advantages of this protocol are:

• An optimized method for extracting good quality DNA from FFPE tissues.
• An efficient bisulfite conversion and next generation sequencing library preparation protocol that uses 50 ng DNA from FFPE tissue.
• Incorporation of a PCR-based test to assess bisulfite conversion efficiency prior to sequencing.

Conclusions

We provide a complete workflow and an integrated protocol for performing DNA methylation analysis at the genome-scale and we believe this will facilitate clinical epigenetic research that involves the use of FFPE tissue.

     

Effect of testosterone replacement on feeding behaviors after acute and chronic stress in gonadectomized male NMRI mice

Background

We study the role of gonadectomy on the response to unavoidable stress and the role of testosterone replacement on gonadectomy in the male Naval Medical Research Institute mice (30±5 g) were studied. For this purpose, the hormonal and metabolic changes were investigated.

Methods

In the experimental group, the gonads were surgically removed, and a cannula was inserted into the left lateral ventricle. For acute and chronic stress induction, animals were placed in the communication box for 30 min for one day and four consecutive days, respectively. The animals received different doses of intraventricular (ICV) testosterone (0.01, 0.05, 0.1 μg/mouse) 5 minutes or intraperitoneal (IP) testosterone (0.05, 0.01, 0.1 mg/kg) 30 minutes before the stress induction.

Results

The results showed that acute and chronic stress increases plasma cortisol concentration. IP testosterone injections of testosterone did not decrease cortisol concentrations in response to acute stress, whereas ICV injections did reduce cortisol concentrations. The stress reduced anorexia time, while the administration of testosterone increased anorexia time. In addition, acute stress reduced food intake in the gonadectomized mice. IP testosterone at 0.01 and 0.05 mg/kg increased food intake. Additionally, stress in gonadectomized mice reduced water intake, while the IP injection of testosterone in chronic stress further reduced water intake. Also, stress reduced the animals’ brain/adrenal volumes, while the IP and ICVinjection of testosterone at 0.01 mg/kg inhibited this effect.

Conclusion

The results showed that the IP (0.05, 0.01, 0.1 mg/kg) and ICV (0.01, 0.05, 0.1 μg/mouse) administration of testosterone in the gonadectomized mice can modulate hormonal and metabolic changes induced by stress.

     

Quality control for Illumina 450K methylation data in the absence of iDat files using correlation structure in pedigrees and repeated measures

Background

An important feature in many genomic studies is quality control and normalization. This is particularly important when analyzing epigenetic data, where the process of obtaining measurements can be bias prone. The GAW20 data was from the Genetics of Lipid Lowering Drugs and Diet Network (GOLDN), a study with multigeneration families, where DNA cytosine-phosphate-guanine (CpG) methylation was measured pre- and posttreatment with fenofibrate. We performed quality control assessment of the GAW20 DNA methylation data, including normalization, assessment of batch effects and detection of sample swaps.

Results

We show that even after normalization, the GOLDN methylation data has systematic differences pre- and posttreatment. Through investigation of (a) CpGs sites containing a single nucleotide polymorphism, (b) the stability of breeding values for methylation across time points, and (c) autosomal gender-associated CpGs, 13 sample swaps were detected, 11 of which were posttreatment.

Conclusions

This paper demonstrates several ways to perform quality control of methylation data in the absence of raw data files and highlights the importance of normalization and quality control of the GAW20 methylation data from the GOLDN study.

     

Pharmacometabolomics applied to zonisamide pharmacokinetic parameter prediction

Introduction

Zonisamide is a new-generation anticonvulsant antiepileptic drug metabolized primarily in the liver, with subsequent elimination via the renal route.

Objectives

Our objective was to evaluate the utility of pharmacometabolomics in the detection of zonisamide metabolites that could be related to its disposition and therefore, to its efficacy and toxicity.

Methods

This study was nested to a bioequivalence clinical trial with 28 healthy volunteers. Each participant received a single dose of zonisamide on two separate occasions (period 1 and period 2), with a washout period between them. Blood samples of zonisamide were obtained from all patients at baseline for each period, before volunteers were administered any medication, for metabolomics analysis.

Results

After a Lasso regression was applied, age, height, branched-chain amino acids, steroids, triacylglycerols, diacyl glycerophosphoethanolamine, glycerophospholipids susceptible to methylation, phosphatidylcholines with 20:4 FA (arachidonic acid) and cholesterol ester and lysophosphatidylcholine were obtained in both periods.

Conclusion

To our knowledge, this is the only research study to date that has attempted to link basal metabolomic status with pharmacokinetic parameters of zonisamide.

     

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.

     

Pseudo current density maps of electrophysiological heart,nerve or brain function and their physical basis

Background

In recent years the visualization of biomagnetic measurement data by so-called pseudo current density maps or Hosaka-Cohen (HC) transformations became popular.

Methods

The physical basis of these intuitive maps is clarified by means of analytically solvable problems.

Results

Examples in magnetocardiography, magnetoencephalography and magnetoneurography demonstrate the usefulness of this method.

Conclusion

Hardware realizations of the HC-transformation and some similar transformations are discussed which could advantageously support cross-platform comparability of biomagnetic measurements.

     

An optimized rapid bisulfite conversion method with high recovery of cell-free DNA

Background

Methylation analysis of cell-free DNA is a encouraging tool for tumor diagnosis, monitoring and prognosis. Sensitivity of methylation analysis is a very important matter due to the tiny amounts of cell-free DNA available in plasma. Most current methods of DNA methylation analysis are based on the difference of bisulfite-mediated deamination of cytosine between cytosine and 5-methylcytosine. However, the recovery of bisulfite-converted DNA based on current methods is very poor for the methylation analysis of cell-free DNA.

Results

We optimized a rapid method for the crucial steps of bisulfite conversion with high recovery of cell-free DNA. A rapid deamination step and alkaline desulfonation was combined with the purification of DNA on a silica column. The conversion efficiency and recovery of bisulfite-treated DNA was investigated by the droplet digital PCR. The optimization of the reaction results in complete cytosine conversion in 30 min at 70 °C and about 65% of recovery of bisulfite-treated cell-free DNA, which is higher than current methods.

Conclusions

The method allows high recovery from low levels of bisulfite-treated cell-free DNA, enhancing the analysis sensitivity of methylation detection from cell-free DNA.

     

Untargeted metabolomics suffers from incomplete raw data processing

Introduction

Untargeted metabolomics is a powerful tool for biological discoveries. To analyze the complex raw data, significant advances in computational approaches have been made, yet it is not clear how exhaustive and reliable the data analysis results are.

Objectives

Assessment of the quality of raw data processing in untargeted metabolomics.

Methods

Five published untargeted metabolomics studies, were reanalyzed.

Results

Omissions of at least 50 relevant compounds from the original results as well as examples of representative mistakes were reported for each study.

Conclusion

Incomplete raw data processing shows unexplored potential of current and legacy data.

     

Metabolomic profiling of maternal hair suggests rapid development of intrahepatic cholestasis of pregnancy

Introduction

Intrahepatic cholestasis of pregnancy (ICP) is a common maternal liver disease; development can result in devastating consequences, including sudden fetal death and stillbirth. Currently, recognition of ICP only occurs following onset of clinical symptoms.

Objective

Investigate the maternal hair metabolome for predictive biomarkers of ICP.

Methods

The maternal hair metabolome (gestational age of sampling between 17 and 41 weeks) of 38 Chinese women with ICP and 46 pregnant controls was analysed using gas chromatography–mass spectrometry.

Results

Of 105 metabolites detected in hair, none were significantly associated with ICP.

Conclusion

Hair samples represent accumulative environmental exposure over time. Samples collected at the onset of ICP did not reveal any metabolic shifts, suggesting rapid development of the disease.

     

How the chemical features of molecules may have addressed the settlement of metabolic steps

Introduction

While the evolutionary adaptation of enzymes to their own substrates is a well assessed and rationalized field, how molecules have been originally selected in order to initiate and assemble convenient metabolic pathways is a fascinating, but still debated argument.

Objectives

Aim of the present study is to give a rationale for the preferential selection of specific molecules to generate metabolic pathways.

Methods

The comparison of structural features of molecules, through an inductive methodological approach, offer a reading key to cautiously propose a determining factor for their metabolic recruitment.

Results

Starting with some commonplaces occurring in the structural representation of relevant carbohydrates, such as glucose, fructose and ribose, arguments are presented in associating stable structural determinants of these molecules and their peculiar occurrence in metabolic pathways.

Conclusions

Among other possible factors, the reliability of the structural asset of a molecule may be relevant or its selection among structurally and, a priori, functionally similar molecules.

     

Quantitative analysis of tetrahydrofolate metabolites from <Emphasis Type="Italic">clostridium autoethanogenum</Emphasis>

Introduction

Quantification of tetrahydrofolates (THFs), important metabolites in the Wood–Ljungdahl pathway (WLP) of acetogens, is challenging given their sensitivity to oxygen.

Objective

To develop a simple anaerobic protocol to enable reliable THFs quantification from bioreactors.

Methods

Anaerobic cultures were mixed with anaerobic acetonitrile for extraction. Targeted LC–MS/MS was used for quantification.

Results

Tetrahydrofolates can only be quantified if sampled anaerobically. THF levels showed a strong correlation to acetyl-CoA, the end product of the WLP.

Conclusion

Our method is useful for relative quantification of THFs across different growth conditions. Absolute quantification of THFs requires the use of labelled standards.