Functional DNA methylation signatures for autism spectrum disorder genomic risk loci: 16p11.2 deletions and CHD8 variants

Background: Autism spectrum disorder (ASD) is a common and etiologically heterogeneous neurodevelopmental disorder. Although many genetic causes have been identified (\u3e 200 ASD-risk genes), no single gene variant accounts for \u3e 1% of all ASD cases. A role for epigenetic mechanisms in ASD etiology is supported by the fact that many ASD-risk genes function as epigenetic regulators and evidence that epigenetic dysregulation can interrupt normal brain development. Gene-specific DNAm profiles have been shown to assist in the interpretation of variants of unknown significance. Therefore, we investigated the epigenome in patients with ASD or two of the most common genomic variants conferring increased risk for ASD. Genome-wide DNA methylation (DNAm) was assessed using the Illumina Infinium HumanMethylation450 and MethylationEPIC arrays in blood from individuals with ASD of heterogeneous, undefined etiology (n = 52), and individuals with 16p11.2 deletions (16p11.2del, n = 9) or pathogenic variants in the chromatin modifier CHD8 (CHD8 +/-, n = 7). Results: DNAm patterns did not clearly distinguish heterogeneous ASD cases from controls. However, the homogeneous genetically-defined 16p11.2del and CHD8 +/- subgroups each exhibited unique DNAm signatures that distinguished 16p11.2del or CHD8 +/- individuals from each other and from heterogeneous ASD and control groups with high sensitivity and specificity. These signatures also classified additional 16p11.2del (n = 9) and CHD8 (n = 13) variants as pathogenic or benign. Our findings that DNAm alterations in each signature target unique genes in relevant biological pathways including neural development support their functional relevance. Furthermore, genes identified in our CHD8 +/- DNAm signature in blood overlapped differentially expressed genes in CHD8 +/- human-induced pluripotent cell-derived neurons and cerebral organoids from independent studies. Conclusions: DNAm signatures can provide clinical utility complementary to next-generation sequencing in the interpretation of variants of unknown significance. Our study constitutes a novel approach for ASD risk-associated molecular classification that elucidates the vital cross-talk between genetics and epigenetics in the etiology of ASD

Deterministic mechanical model of T-killer cell polarization reproduces the wandering of aim between simultaneously engaged targets

T-killer cells of the immune system eliminate virus-infected and tumorous cells through direct cell-cell interactions. Reorientation of the killing apparatus inside the T cell to the T-cell interface with the target cell ensures specificity of the immune response. The killing apparatus can also oscillate next to the cell-cell interface. When two target cells are engaged by the T cell simultaneously, the killing apparatus can oscillate between the two interface areas. This oscillation is one of the most striking examples of cell movements that give the microscopist an unmechanistic impression of the cell's fidgety indecision. We have constructed a three-dimensional, numerical biomechanical model of the molecular-motor-driven microtubule cytoskeleton that positions the killing apparatus. The model demonstrates that the cortical pulling mechanism is indeed capable of orienting the killing apparatus into the functional position under a range of conditions. The model also predicts experimentally testable limitations of this commonly hypothesized mechanism of T-cell polarization. After the reorientation, the numerical solution exhibits complex, multidirectional, multiperiodic, and sustained oscillations in the absence of any external guidance or stochasticity. These computational results demonstrate that the strikingly animate wandering of aim in T-killer cells has a purely mechanical and deterministic explanation. © 2009 Kim, Maly

Chromosome Driven Spatial Patterning of Proteins in Bacteria

The spatial patterning of proteins in bacteria plays an important role in many processes, from cell division to chemotaxis. In the asymmetrically dividing bacteria Caulobacter crescentus, a scaffolding protein, PopZ, localizes to both poles and aids the differential patterning of proteins between mother and daughter cells during division. Polar patterning of misfolded proteins in Escherechia coli has also been shown, and likely plays an important role in cellular ageing. Recent experiments on both of the above systems suggest that the presence of chromosome free regions along with protein multimerization may be a mechanism for driving the polar localization of proteins. We have developed a simple physical model for protein localization using only these two driving mechanisms. Our model reproduces all the observed patterns of PopZ and misfolded protein localization - from diffuse, unipolar, and bipolar patterns and can also account for the observed patterns in a variety of mutants. The model also suggests new experiments to further test the role of the chromosome in driving protein patterning, and whether such a mechanism is responsible for helping to drive the differentiation of the cell poles

NSD1 mutations generate a genome-wide DNA methylation signature

published_or_final_versio

A new multicompartmental reaction-diffusion modeling method links transient membrane attachment of E. coli MinE to E-ring formation

Many important cellular processes are regulated by reaction-diffusion (RD) of molecules that takes place both in the cytoplasm and on the membrane. To model and analyze such multicompartmental processes, we developed a lattice-based Monte Carlo method, Spatiocyte that supports RD in volume and surface compartments at single molecule resolution. Stochasticity in RD and the excluded volume effect brought by intracellular molecular crowding, both of which can significantly affect RD and thus, cellular processes, are also supported. We verified the method by comparing simulation results of diffusion, irreversible and reversible reactions with the predicted analytical and best available numerical solutions. Moreover, to directly compare the localization patterns of molecules in fluorescence microscopy images with simulation, we devised a visualization method that mimics the microphotography process by showing the trajectory of simulated molecules averaged according to the camera exposure time. In the rod-shaped bacterium _Escherichia coli_, the division site is suppressed at the cell poles by periodic pole-to-pole oscillations of the Min proteins (MinC, MinD and MinE) arising from carefully orchestrated RD in both cytoplasm and membrane compartments. Using Spatiocyte we could model and reproduce the _in vivo_ MinDE localization dynamics by accounting for the established properties of MinE. Our results suggest that the MinE ring, which is essential in preventing polar septation, is largely composed of MinE that is transiently attached to the membrane independently after recruited by MinD. Overall, Spatiocyte allows simulation and visualization of complex spatial and reaction-diffusion mediated cellular processes in volumes and surfaces. As we showed, it can potentially provide mechanistic insights otherwise difficult to obtain experimentally

DNA Methylation Signature for EZH2 Functionally Classifies Sequence Variants in Three PRC2 Complex Genes.

Weaver syndrome (WS), an overgrowth/intellectual disability syndrome (OGID), is caused by pathogenic variants in the histone methyltransferase EZH2, which encodes a core component of the Polycomb repressive complex-2 (PRC2). Using genome-wide DNA methylation (DNAm) data for 187 individuals with OGID and 969 control subjects, we show that pathogenic variants in EZH2 generate a highly specific and sensitive DNAm signature reflecting the phenotype of WS. This signature can be used to distinguish loss-of-function from gain-of-function missense variants and to detect somatic mosaicism. We also show that the signature can accurately classify sequence variants in EED and SUZ12, which encode two other core components of PRC2, and predict the presence of pathogenic variants in undiagnosed individuals with OGID. The discovery of a functionally relevant signature with utility for diagnostic classification of sequence variants in EZH2, EED, and SUZ12 supports the emerging paradigm shift for implementation of DNAm signatures into diagnostics and translational research

Equilibria of Idealized Confined Astral Microtubules and Coupled Spindle Poles

Positioning of the mitotic spindle through the interaction of astral microtubules with the cell boundary often determines whether the cell division will be symmetric or asymmetric. This process plays a crucial role in development. In this paper, a numerical model is presented that deals with the force exerted on the spindle by astral microtubules that are bent by virtue of their confinement within the cell boundary. It is found that depending on parameters, the symmetric position of the spindle can be stable or unstable. Asymmetric stable equilibria also exist, and two or more stable positions can exist simultaneously. The theory poses new types of questions for experimental research. Regarding the cases of symmetric spindle positioning, it is necessary to ask whether the microtubule parameters are controlled by the cell so that the bending mechanics favors symmetry. If they are not, then it is necessary to ask what forces external to the microtubule cytoskeleton counteract the bending effects sufficiently to actively establish symmetry. Conversely, regarding the cases with asymmetry, it is now necessary to investigate whether the cell controls the microtubule parameters so that the bending favors asymmetry apart from any forces that are external to the microtubule cytoskeleton

Uncovering the Dynamics of Cardiac Systems Using Stochastic Pacing and Frequency Domain Analyses

Alternans of cardiac action potential duration (APD) is a well-known arrhythmogenic mechanism which results from dynamical instabilities. The propensity to alternans is classically investigated by examining APD restitution and by deriving APD restitution slopes as predictive markers. However, experiments have shown that such markers are not always accurate for the prediction of alternans. Using a mathematical ventricular cell model known to exhibit unstable dynamics of both membrane potential and Ca2+ cycling, we demonstrate that an accurate marker can be obtained by pacing at cycle lengths (CLs) varying randomly around a basic CL (BCL) and by evaluating the transfer function between the time series of CLs and APDs using an autoregressive-moving-average (ARMA) model. The first pole of this transfer function corresponds to the eigenvalue (λalt) of the dominant eigenmode of the cardiac system, which predicts that alternans occurs when λalt≤−1. For different BCLs, control values of λalt were obtained using eigenmode analysis and compared to the first pole of the transfer function estimated using ARMA model fitting in simulations of random pacing protocols. In all versions of the cell model, this pole provided an accurate estimation of λalt. Furthermore, during slow ramp decreases of BCL or simulated drug application, this approach predicted the onset of alternans by extrapolating the time course of the estimated λalt. In conclusion, stochastic pacing and ARMA model identification represents a novel approach to predict alternans without making any assumptions about its ionic mechanisms. It should therefore be applicable experimentally for any type of myocardial cell

An Experimental and Computational Study of Effects of Microtubule Stabilization on T-Cell Polarity

T-killer cells eliminate infected and cancerous cells with precision by positioning their centrosome near the interface (immunological synapse) with the target cell. The mechanism of centrosome positioning has remained controversial, in particular the role of microtubule dynamics in it. We re-examined the issue in the experimental model of Jurkat cells presented with a T cell receptor-binding artificial substrate, which permits controlled stimulation and reproducible measurements. Neither 1-µM taxol nor 100-nM nocodazole inhibited the centrosome positioning at the “synapse” with the biomimetic substrate. At the same time, in micromolar taxol but not in nanomolar nocodazole the centrosome adopted a distinct peripheral rather than the normally central position within the synapse. This effect was reproduced in a computational energy-minimization model that assumed no microtubule dynamics, but only a taxol-induced increase in the length of the microtubules. Together, the experimental and computational results indicate that microtubule dynamics are not essential for the centrosome positioning, but that the fit of the microtubule array in the deformed body of the conjugated T cell is a major factor. The possibility of modulating the T-cell centrosome position with well-studied drugs and of predicting their effects in silico appears attractive for designing anti-cancer and antiviral therapies

Truncating SRCAP variants outside the Floating-Harbor syndrome locus cause a distinct neurodevelopmental disorder with a specific DNA methylation signature

Truncating variants in exons 33 and 34 of the SNF2-related CREBBP activator protein (SRCAP) gene cause the neurodevelopmental disorder (NDD) Floating-Harbor syndrome (FLHS), characterized by short stature, speech delay, and facial dysmorphism. Here, we present a cohort of 33 individuals with clinical features distinct from FLHS and truncating (mostly de novo) SRCAP variants either proximal (n = 28) or distal (n = 5) to the FLHS locus. Detailed clinical characterization of the proximal SRCAP individuals identified shared characteristics: developmental delay with or without intellectual disability, behavioral and psychiatric problems, non-specific facial features, musculoskeletal issues, and hypotonia. Because FLHS is known to be associated with a unique set of DNA methylation (DNAm) changes in blood, a DNAm signature, we investigated whether there was a distinct signature associated with our affected individuals. A machine-learning model, based on the FLHS DNAm signature, negatively classified all our tested subjects. Comparing proximal variants with typically developing controls, we identified a DNAm signature distinct from the FLHS signature. Based on the DNAm and clinical data, we refer to the condition as "non-FLHS SRCAP-related NDD.'' All five distal variants classified negatively using the FLHS DNAm model while two classified positively using the proximal model. This suggests divergent pathogenicity of these variants, though clinically the distal group presented with NDD, similar to the proximal SRCAP group. In summary, for SRCAP, there is a clear relationship between variant location, DNAm profile, and clinical phenotype. These results highlight the power of combined epigenetic, molecular, and clinical studies to identify and characterize genotype-epigenotype-phenotype correlations