Statistical inference of transmission fidelity of DNA methylation patterns over somatic cell divisions in mammals

We develop Bayesian inference methods for a recently-emerging type of epigenetic data to study the transmission fidelity of DNA methylation patterns over cell divisions. The data consist of parent-daughter double-stranded DNA methylation patterns with each pattern coming from a single cell and represented as an unordered pair of binary strings. The data are technically difficult and time-consuming to collect, putting a premium on an efficient inference method. Our aim is to estimate rates for the maintenance and de novo methylation events that gave rise to the observed patterns, while accounting for measurement error. We model data at multiple sites jointly, thus using whole-strand information, and considerably reduce confounding between parameters. We also adopt a hierarchical structure that allows for variation in rates across sites without an explosion in the effective number of parameters. Our context-specific priors capture the expected stationarity, or near-stationarity, of the stochastic process that generated the data analyzed here. This expected stationarity is shown to greatly increase the precision of the estimation. Applying our model to a data set collected at the human FMR1 locus, we find that measurement errors, generally ignored in similar studies, occur at a nontrivial rate (inappropriate bisulfite conversion error: 1.6$$ with 80$$ CI: 0.9--2.3$$). Accounting for these errors has a substantial impact on estimates of key biological parameters. The estimated average failure of maintenance rate and daughter de novo rate decline from 0.04 to 0.024 and from 0.14 to 0.07, respectively, when errors are accounted for. Our results also provide evidence that de novo events may occur on both parent and daughter strands: the median parent and daughter de novo rates are 0.08 (80$$ CI: 0.04--0.13) and 0.07 (80$$ CI: 0.04--0.11), respectively.Comment: Published in at http://dx.doi.org/10.1214/09-AOAS297 the Annals of Applied Statistics (http://www.imstat.org/aoas/) by the Institute of Mathematical Statistics (http://www.imstat.org

High-level dietary cadmium exposure is associated with global DNA hypermethylation in the gastropod hepatopancreas

5-methylcytosine (5mC) is a key epigenetic mark which influences gene expression and phenotype. In vertebrates, this epigenetic mark is sensitive to Cd exposure, but there is no information linking such an event with changes in global 5mC levels in terrestrial gastropods despite their importance as excellent ecotoxicological bioindicators of metal contamination. Therefore, we first evaluated total 5mC content in DNA of the hepatopancreas of adult Cantareus aspersus with the aim to determine whether this epigenetic mark is responsive to Cd exposure. The experiment was conducted under laboratory conditions and involved a continuous exposure, multiple dose- and time-point (14, 28, and 56 days) study design. Hepatopancreas cadmium levels were measured using Flame Atomic Absorption Spectrometry and the percentage of 5-mC in samples using an ELISA-based colorimetric assay. Snail death rates were also assessed. Our results, for the first time, reveal the presence of 5mC in C. aspersus and provide evidence for Cd-induced changes in global 5mC levels in DNA of gastropods and mollusks. Although less sensitive than tissue accumulation, DNA methylation levels responded in a dose- and time-dependent manner to dietary cadmium, with exposure dose having a much stronger effect than exposure duration. An obvious trend of increasing 5mC levels was observed starting at 28 days of exposure to the second highest dose and this trend persisted at the two highest treatments for close to one month, when the experiment was terminated after 56 days. Moreover, a strong association was identified between Cd concentrations in the hepatopancreas and DNA methylation levels in this organ. These data indicate an overall trend towards DNA hypermethylation with elevated Cd exposure. No consistent lethal effect was observed, irrespective of time point and Cd-dosage. Overall, our findings suggest that the total 5mC content in DNA of the hepatopancreas of land snails is responsive to sublethal Cd exposure and give new insights into invertebrate environmental epigenetics

Acute imidacloprid exposure alters mitochondrial function in bumblebee flight muscle and brain

Peer reviewe

Global DNA methylation profiling of manganese-exposed human neuroblastoma SH-SY5Y cells reveals epigenetic alterations in Parkinson’s disease-associated genes

Manganese (Mn) is an essential trace element required for optimal functioning of cellular biochemical pathways in the central nervous system. Elevated exposure to Mn through environmental and occupational exposure can cause neurotoxic effects resulting in manganism, a condition with clinical symptoms identical to idiopathic Parkinson’s disease. Epigenetics is now recognized as a biological mechanism involved in the etiology of various diseases. Here, we investigated the role of DNA methylation alterations induced by chronic Mn (100 µM) exposure in human neuroblastoma (SH-SY5Y) cells in relevance to Parkinson’s disease. A combined analysis of DNA methylation and gene expression data for Parkinson’s disease-associated genes was carried out. Whole-genome bisulfite conversion and sequencing indicate epigenetic perturbation of key genes involved in biological processes associated with neuronal cell health. Integration of DNA methylation data with gene expression reveals epigenetic alterations to PINK1, PARK2 and TH genes that play critical roles in the onset of Parkinsonism. The present study suggests that Mn-induced alteration of DNA methylation of PINK1–PARK2 may influence mitochondrial function and promote Parkinsonism. Our findings provide a basis to further explore and validate the epigenetic basis of Mn-induced neurotoxicity

Epigenetic regulation of 5α reductase-1 underlies adaptive plasticity of reproductive function and pubertal timing

Women facing increased energetic demands in childhood commonly have altered adult ovarian activity and shorter reproductive lifespan, possibly comprising a strategy to optimize reproductive success. Here we sought to understand the mechanisms of early-life programming of reproductive function, by integrating analysis of reproductive tissues in an appropriate mouse model with methylation analysis of proxy tissue DNA in a well-characterized population of Bangladeshi migrants in the UK. Bangladeshi women whose childhood was in Bangladesh were found to have later pubertal onset and lower age-matched ovarian reserve than Bangladeshi women who grew-up in England. Subsequently we aimed to explore the potential relevance to the altered reproductive phenotype of one of the genes that emerged from the screens. Results: Of the genes associated with differential methylation in the Bangladeshi women whose childhood was in Bangladesh as compared to Bangladeshi women who grew up in the UK, 13 correlated with altered expression of the orthologous gene in the mouse model ovaries. These mice had delayed pubertal onset and a smaller ovarian reserve compared to controls. The most relevant of these genes for reproductive function appeared to be SRD5A1, which encodes the steroidogenic enzyme 5α reductase-1. SRD5A1 was more methylated at the same transcriptional enhancer in mice ovaries as in the women’s buccal DNA, and its expression was lower in the hypothalamus of the mice as well, suggesting a possible role in the central control of reproduction. The expression of Kiss1 and Gnrh was also lower in these mice compared to controls, and inhibition of 5α reductase-1 reduced Kiss1 and Gnrh mRNA levels and blocked GnRH release in GnRH neuronal cell cultures. Crucially, we show that inhibition of this enzyme in female mice in vivo delayed pubertal onset. Conclusions: SRD5A1/5α reductase-1 responds epigenetically to the environment and its down-regulation appears to alter the reproductive phenotype. These findings help to explain diversity in reproductive characteristics and how they are shaped by early-life environment, and reveal novel pathways that might be targeted to mitigate health issues caused by life-history trade-offs

Manganese exposure: linking down-regulation of miRNA-7 and miRNA-433 with α-synuclein overexpression and risk of idiopathic Parkinson's disease

Manganese is an essential trace element however elevated environmental and occupational exposure to this element has been correlated with neurotoxicity symptoms clinically identical to idiopathic Parkinson's disease. In the present study we chronically exposed human neuroblastoma SH-SY5Y cells to manganese (100 μM) and carried out expression profiling of miRNAs known to modulate neuronal differentiation and neurodegeneration. The miRNA PCR array results reveal alterations in expression levels of miRNAs, which have previously been associated with the regulation of synaptic transmission and apoptosis. The expressions of miR-7 and miR-433 significantly reduced upon manganese exposure. By in silico homology analysis we identified SNCA and FGF-20as targets of miR-7 and miR-433. We demonstrate an inverse correlation in expression levels where reduction in these two miRNAs causes increases in SNCA and FGF-20. Transient transfection of SH-SY5Y cells with miR-7 and miR-433 mimics resulted in down regulation of SNCA and FGF-20 mRNA levels. Our study is the first to uncover the potential link between manganese exposure, altered miRNA expression and parkinsonism: manganese exposure causes overexpression of SNCA and FGF-20 by diminishing miR-7 and miR-433 levels. These miRNAs may be considered critical for protection from manganese induced neurotoxic mechanism and hence as potential therapeutic targets

Thiamethoxam exposure deregulates short ORF gene expression in the honey bee and compromises immune response to bacteria

© 2021, The Author(s). Maximizing crop yields relies on the use of agrochemicals to control insect pests. One of the most widely used classes of insecticides are neonicotinoids that interfere with signalling of the neurotransmitter acetylcholine, but these can also disrupt crop-pollination services provided by bees. Here, we analysed whether chronic low dose long-term exposure to the neonicotinoid thiamethoxam alters gene expression and alternative splicing in brains of Africanized honey bees, Apis mellifera, as adaptation to altered neuronal signalling. We find differentially regulated genes that show concentration-dependent responses to thiamethoxam, but no changes in alternative splicing. Most differentially expressed genes have no annotated function but encode short Open Reading Frames, a characteristic feature of anti-microbial peptides. As this suggested that immune responses may be compromised by thiamethoxam exposure, we tested the impact of thiamethoxam on bee immunity by injecting bacteria. We show that intrinsically sub-lethal thiamethoxam exposure makes bees more vulnerable to normally non-pathogenic bacteria. Our findings imply a synergistic mechanism for the observed bee population declines that concern agriculturists, conservation ecologists and the public

Testing the FMR1 Promoter for Mosaicism in DNA Methylation among CpG Sites, Strands, and Cells in FMR1-Expressing Males with Fragile X Syndrome

Variability among individuals in the severity of fragile X syndrome (FXS) is influenced by epigenetic methylation mosaicism, which may also be common in other complex disorders. The epigenetic signal of dense promoter DNA methylation is usually associated with gene silencing, as was initially reported for FMR1 alleles in individuals with FXS. A paradox arose when significant levels of FMR1 mRNA were reported for some males with FXS who had been reported to have predominately methylated alleles. We have used hairpin-bisufite PCR, validated with molecular batch-stamps and barcodes, to collect and assess double-stranded DNA methylation patterns from these previously studied males. These patterns enable us to distinguish among three possible forms of methylation mosaicism, any one of which could explain FMR1 expression in these males. Our data indicate that cryptic inter-cell mosaicism in DNA methylation can account for the presence of FMR1 mRNA in some individuals with FXS

Statistical Inference of In Vivo Properties of Human DNA Methyltransferases from Double-Stranded Methylation Patterns

DNA methyltransferases establish methylation patterns in cells and transmit these patterns over cell generations, thereby influencing each cell's epigenetic states. Three primary DNA methyltransferases have been identified in mammals: DNMT1, DNMT3A and DNMT3B. Extensive in vitro studies have investigated key properties of these enzymes, namely their substrate specificity and processivity. Here we study these properties in vivo, by applying novel statistical analysis methods to double-stranded DNA methylation patterns collected using hairpin-bisulfite PCR. Our analysis fits a novel Hidden Markov Model (HMM) to the observed data, allowing for potential bisulfite conversion errors, and yields statistical estimates of parameters that quantify enzyme processivity and substrate specificity. We apply this model to methylation patterns established in vivo at three loci in humans: two densely methylated inactive X (Xi)-linked loci ( and ), and an autosomal locus (), where methylation densities are tissue-specific but moderate. We find strong evidence for a high level of processivity of DNMT1 at and , with the mean association tract length being a few hundred base pairs. Regardless of tissue types, methylation patterns at are dominated by DNMT1 maintenance events, similar to the two Xi-linked loci, but are insufficiently informative regarding processivity to draw any conclusions about processivity at that locus. At all three loci we find that DNMT1 shows a strong preference for adding methyl groups to hemi-methylated CpG sites over unmethylated sites. The data at all three loci also suggest low (possibly 0) association of the de novo methyltransferases, the DNMT3s, and are consequently uninformative about processivity or preference of these enzymes. We also extend our HMM to reanalyze published data on mouse DNMT1 activities in vitro. The results suggest shorter association tracts (and hence weaker processivity), and much longer non-association tracts than human DNMT1 in vivo