Whole exome sequencing has proven to be a powerful tool for understanding the genetic architecture of human disease. Here we apply it to more than 2,500 simplex families, each having a child with an autistic spectrum disorder. By comparing affected to unaffected siblings, we show that 13% of de novo missense mutations and 43% of de novo likely gene-disrupting (LGD) mutations contribute to 12% and 9% of diagnoses, respectively. Including copy number variants, coding de novo mutations contribute to about 30% of all simplex and 45% of female diagnoses. Almost all LGD mutations occur opposite wild-type alleles. LGD targets in affected females significantly overlap the targets in males of lower intelligence quotient (IQ), but neither overlaps significantly with targets in males of higher IQ. We estimate that LGD mutation in about 400 genes can contribute to the joint class of affected females and males of lower IQ, with an overlapping and similar number of genes vulnerable to contributory missense mutation. LGD targets in the joint class overlap with published targets for intellectual disability and schizophrenia, and are enriched for chromatin modifiers, FMRP-associated genes and embryonically expressed genes. Most of the significance for the latter comes from affected females.

SETD3 is a member of the protein lysine methyltransferase (PKMT) family, which catalyzes the addition of methyl group to lysine residues. Accumulating data suggest that PKMTs are involved in the regulation of a broad spectrum of biological processes by targeting histone and non-histone proteins. Using a proteomic approach, we have identified 172 new SETD3 interacting proteins. We show that SETD3 binds and methylates the transcription factor FoxM1, which has been previously shown to be associated with the regulation of VEGF expression. We further demonstrate that under hypoxic conditions SETD3 is down-regulated. Mechanistically, we find that under basal conditions, SETD3 and FoxM1 are enriched on the VEGF promoter. Dissociation of both SETD3 and FoxM1 from the VEGF promoter under hypoxia correlates with elevated expression of VEGF. Taken together, our data reveal a new SETD3-dependent methylation-based signaling pathway at chromatin that regulates VEGF expression under normoxic and hypoxic conditions.

SET domain containing 6 (SETD6) monomethylates the RelA subunit of nuclear factor kappa B (NF-κB). The ankyrin repeats of G9a-like protein (GLP) recognizes RelA monomethylated at Lys310. Adjacent to Lys310 is Ser311, a known phosphorylation site of RelA. Ser311 phosphorylation inhibits Lys310 methylation by SETD6 as well as binding of Lys310me1 by GLP. The structure of SETD6 in complex with RelA peptide containing the methylation site, in the presence of S-adenosyl-L-methionine, reveals a V-like protein structure and suggests a model for NF-κB binding to SETD6. In addition, structural modeling of the GLP ankyrin repeats bound to Lys310me1 peptide provides insight into the molecular basis for inhibition of Lys310me1 binding by Ser311 phosphorylation. Together, these findings provide a structural explanation for a key cellular signaling pathway centered on RelA Lys310 methylation, which is generated by SETD6 and recognized by GLP, and incorporate a methylation-phosphorylation switch of adjacent lysine and serine residues. Finally, SETD6 is structurally similar to the Rubisco large subunit methyltransferase. Given the restriction of Rubisco to plant species, this particular appearance of the protein lysine methyltransferase has been evolutionarily well conserved.

Identification of complex molecular networks underlying common human phenotypes is a major challenge of modern genetics. In this study, we develop a method for network-based analysis of genetic associations (NETBAG). We use NETBAG to identify a large biological network of genes affected by rare de novo CNVs in autism. The genes forming the network are primarily related to synapse development, axon targeting, and neuron motility. The identified network is strongly related to genes previously implicated in autism and intellectual disability phenotypes. Our results are also consistent with the hypothesis that significantly stronger functional perturbations are required to trigger the autistic phenotype in females compared to males. Overall, the presented analysis of de novo variants supports the hypothesis that perturbed synaptogenesis is at the heart of autism. More generally, our study provides proof of the principle that networks underlying complex human phenotypes can be identified by a network-based functional analysis of rare genetic variants.

P21-activated kinase 4 (PAK4), a member of serine/threonine kinases family is over-expressed in numerous cancer tumors and is associated with oncogenic cell proliferation, migration and invasion. Our recent work demonstrated that the SET-domain containing protein 6 (SETD6) interacts with and methylates PAK4 at chromatin in mammalian cells, leading to activation of the Wnt/β-catenin signaling pathway. In our current work, we identified lysine 473 (K473) on PAK4 as the primary methylation site by SETD6. Methylation of PAK4 at K473 activates β-catenin transcriptional activity and inhibits cell adhesion. Specific methylation of PAK4 at K473 also attenuates paxillin localization to focal adhesions leading to overall reduction in adhesion-related features, such as filopodia and actin structures. The altered adhesion of the PAK4 wild-type cells is accompanied with a decrease in the migrative and invasive characteristics of the cells. Taken together, our results suggest that methylation of PAK4 at K473 plays a vital role in the regulation of cell adhesion and migration.

Gliomas are one of the most common and lethal brain tumors among adults. One process that contributes to glioma progression and recurrence is the epithelial to mesenchymal transition (EMT). EMT is regulated by a set of defined transcription factors which tightly regulate this process, among them is the basic helix-loop-helix family member, TWIST1. Here we show that TWIST1 is methylated on lysine-33 at chromatin by SETD6, a methyltransferase with expression levels correlating with poor survival in glioma patients. RNA-seq analysis in U251 glioma cells suggested that both SETD6 and TWIST1 regulate cell adhesion and migration processes. We further show that TWIST1 methylation attenuates the expression of the long-non-coding RNA, LINC-PINT, thereby promoting EMT in glioma. Mechanistically, TWIST1 methylation represses the transcription of LINC-PINT by increasing the occupancy of EZH2 and the catalysis of the repressive H3K27me3 mark at the LINC-PINT locus. Under un-methylated conditions, TWIST1 dissociates from the LINC-PINT locus, allowing the expression of LINC-PINT which leads to increased cell adhesion and decreased cell migration. Together, our findings unravel a new mechanistic dimension for selective expression of LINC-PINT mediated by TWIST1 methylation.

The genesis of the headache phase in migraine with aura is thought to be mediated by cortical spreading depression (CSD) and the subsequent activation and sensitization of primary afferent neurons that innervate the intracranial meninges and their related large vessels. Yet, the exact mechanisms underlying this peripheral meningeal nociceptive response remain poorly understood. We investigated the relative contribution of cortical astrocytes to CSD-evoked meningeal nociception using extracellular single-unit recording of meningeal afferent activity and 2-photon imaging of cortical astrocyte calcium activity, in combination with 2 pharmacological approaches to inhibit astrocytic function. We found that fluoroacetate and l-α-aminoadipate, which inhibit astrocytes through distinct mechanisms, suppressed CSD-evoked afferent mechanical sensitization, but did not affect afferent activation. Pharmacological inhibition of astrocytic function, which ameliorated meningeal afferents' sensitization, reduced basal astrocyte calcium activity but had a minimal effect on the astrocytic calcium wave during CSD. We propose that calcium-independent signaling in cortical astrocytes plays an important role in driving the sensitization of meningeal afferents and the ensuing intracranial mechanical hypersensitivity in migraine with aura.

Migraine headache is hypothesized to involve the activation and sensitization of trigeminal sensory afferents that innervate the cranial meninges. To better understand migraine pathophysiology and improve clinical translation, we used two-photon calcium imaging via a closed cranial window in awake mice to investigate changes in the responses of meningeal afferent fibers using a preclinical model of migraine involving cortical spreading depolarization (CSD). A single CSD episode caused a seconds-long wave of calcium activation that propagated across afferents and along the length of individual afferents. Surprisingly, unlike previous studies in anesthetized animals with exposed meninges, only a very small afferent population was persistently activated in our awake mouse preparation, questioning the relevance of this neuronal response to the onset of migraine pain. In contrast, we identified a larger subset of meningeal afferents that developed augmented responses to acute three-dimensional meningeal deformations that occur in response to locomotion bouts. We observed increased responsiveness in a subset of afferents that were already somewhat sensitive to meningeal deformation before CSD. Furthermore, another subset of previously insensitive afferents also became sensitive to meningeal deformation following CSD. Our data provides new insights into the mechanisms underlying migraine, including the emergence of enhanced meningeal afferent responses to movement-related meningeal deformations as a potential neural substrate underlying the worsening of migraine headache during physical activity.

Migraine headache is hypothesized to involve the activation and sensitization of trigeminal sensory afferents that innervate the cranial meninges. To better understand migraine pathophysiology and improve clinical translation, we used two-photon calcium imaging via a closed cranial window in awake mice to investigate changes in the responses of meningeal afferent fibers using a preclinical model of migraine involving cortical spreading depolarization (CSD). A single CSD episode caused a seconds-long wave of calcium activation that propagated across afferents and along the length of individual afferents. Surprisingly, unlike previous studies in anesthetized animals with exposed meninges, only a very small afferent population was persistently activated in our awake mouse preparation, questioning the relevance of this neuronal response to the onset of migraine pain. In contrast, we identified a larger subset of meningeal afferents that developed augmented responses to acute three-dimensional meningeal deformations that occur in response to locomotion bouts. We observed increased responsiveness in a subset of afferents that were already somewhat sensitive to meningeal deformation before CSD. Furthermore, another subset of previously insensitive afferents also became sensitive to meningeal deformation following CSD. Our data provides new insights into the mechanisms underlying migraine, including the emergence of enhanced meningeal afferent responses to movement-related meningeal deformations as a potential neural substrate underlying the worsening of migraine headache during physical activity.

Signaling via protein lysine methylation has been proposed to play a central role in the regulation of many physiologic and pathologic programs. In contrast to other post-translational modifications such as phosphorylation, proteome-wide approaches to investigate lysine methylation networks do not exist.

High levels of cancer aneuploidy are frequently associated with poor prognosis. To examine the relationship between aneuploidy and cancer progression, we analyzed a series of congenic cell lines that harbor single extra chromosomes. We found that across 13 different trisomic cell lines, 12 trisomies suppressed invasiveness or were largely neutral, while a single trisomy increased metastatic behavior by triggering a partial epithelial-mesenchymal transition. In contrast, we discovered that chromosomal instability activates cGAS/STING signaling but strongly suppresses invasiveness. By analyzing patient copy-number data, we demonstrate that specific aneuploidies are associated with distinct outcomes, and the acquisition of certain aneuploidies is in fact linked with a favorable prognosis. Thus, aneuploidy is not a uniform driver of malignancy, and different aneuploidies can uniquely influence tumor progression. At the same time, the gain of a single chromosome is capable of inducing a profound cell state transition, thereby linking genomic plasticity, phenotypic plasticity, and metastasis.

We show the use of 5'-Acrydite oligonucleotides to copolymerize single-cell DNA or RNA into balls of acrylamide gel (BAGs). Combining this step with split-and-pool techniques for creating barcodes yields a method with advantages in cost and scalability, depth of coverage, ease of operation, minimal cross-contamination, and efficient use of samples. We perform DNA copy number profiling on mixtures of cell lines, nuclei from frozen prostate tumors, and biopsy washes. As applied to RNA, the method has high capture efficiency of transcripts and sufficient consistency to clearly distinguish the expression patterns of cell lines and individual nuclei from neurons dissected from the mouse brain. By using varietal tags (UMIs) to achieve sequence error correction, we show extremely low levels of cross-contamination by tracking source-specific SNVs. The method is readily modifiable, and we will discuss its adaptability and diverse applications.

Brain trauma was clinically associated with increased osteogenesis in the appendicular skeleton. We showed previously in C57BL/6J mice that mild traumatic brain injury (mTBI) transiently induced bone formation in the femur via the cannabinoid-1 (CB1) receptor. Here, we subjected ICR mice to mTBI and examined the bone response in the skull using microCT. We also measured mast cell degranulation (MCD)72 h post-injury. Finally, we measured brain and calvarial endocannabinoids levels post-mTBI. mTBI led to decreased bone porosity on the contralateral (untouched) side. This effect was apparent both in young and mature mice. Administration of rimonabant (CB1 inverse agonist) completely abrogated the effect of mTBI on calvarial porosity and significantly reduced MCD, compared with vehicle-treated controls. We also found that mTBI resulted in elevated levels of anandamide, but not 2-arachidonoylglycerol, in the contralateral calvarial bone, whereas brain levels remained unchanged. In C57BL/6J CB1 knockout mice, mTBI did not reduce porosity but in general the porosity was significantly lower than in WT controls. Our findings suggest that mTBI induces a strain-specific CB1-dependent bone anabolic response in the skull, probably mediated by anandamide, but seemingly unrelated to inflammation. The endocannabinoid system is therefore a plausible target in management of bone response following head trauma.

The protein lysine methyltransferase SET domain-containing protein 6 (SETD6) has been shown to influence different cellular activities and to be critically involved in the regulation of diverse developmental and pathological processes. However, the upstream signals that regulate the mRNA expression of SETD6 are not known. Bioinformatic analysis revealed that the SETD6 promoter has a binding site for the transcription factor E2F1. Using various experimental approaches, we show that E2F1 binds to the SETD6 promoter and regulates SETD6 mRNA expression. Our further observation that this phenomenon is SETD6 dependent suggested that SETD6 and E2F1 are linked. We next demonstrate that SETD6 monomethylates E2F1 specifically at K117 in vitro and in cells. Finally, we show that E2F1 methylation at K117 positively regulates the expression level of SETD6 mRNA. Depletion of SETD6 or overexpression of E2F1 K117R mutant, which cannot be methylated by SETD6, reverses the effect. Taken together, our data provide evidence for a positive feedback mechanism, which regulates the expression of SETD6 by E2F1 in a SETD6 methylation-dependent manner, and highlight the importance of protein lysine methyltransferases and lysine methylation signaling in the regulation of gene transcription.

Signaling via lysine methylation by protein lysine methyltransferases (PKMTs), has been linked to diverse biological and disease processes. The mono-methyltransferase SETD6 (SET-domain-containing protein 6) is a member of the PKMT family and was previously shown to regulate essential cellular processes such as the NF-κB, WNT and the oxidative stress pathways. However, on the biochemical level, little is known about the enzymatic mode of action of SETD6. Here we provide evidence that SETD6 forms high-molecular-weight structures. Specifically, we demonstrate that SETD6 monomeric, dimeric and trimeric forms are stabilized by the methyl donor, S-adenosyl-l-methionine. We then show that SETD6 has auto-methylation activity at K39 and K179, which serves as the major auto-methylation sites with a moderate auto-methylation activity toward K372. A point mutation at K179 but not at K39 and K372, located at the SET domain of SETD6, impaired SETD6 ability to form a trimer, strongly implying a link between the auto-methylation and the oligomerization state. Finally, by radioactive in vitro methylation experiments and biochemical kinetics analysis, we show that the auto-methylation at K39 and K179 increases the catalytic rate of SETD6. Collectively, our data support a model by which SETD6 auto-methylation and self-interaction positively regulate its enzymatic activity in vitro and may suggest that other PKMTs are regulated in the same manner.

SETD3 is a member of the protein lysine methyltransferase (PKMT) family, which catalyzes the addition of methyl group to lysine residues. However, the protein network and the signaling pathways in which SETD3 is involved remain largely unexplored. In the current study, we show that SETD3 is a positive regulator of DNA-damage-induced apoptosis in colon cancer cells. Our data indicate that depletion of SETD3 from HCT-116 cells results in a significant inhibition of apoptosis after doxorubicin treatment. Our results imply that the positive regulation is sustained by methylation, though the substrate remains unknown. We present a functional cross-talk between SETD3 and the tumor suppressor p53. SETD3 binds p53 in cells in response to doxorubicin treatment and positively regulates p53 target genes activation under these conditions. Mechanistically, we provide evidence that the presence of SETD3 and its catalytic activity is required for the recruitment of p53 to its target genes. Finally, Kaplan-Meier survival analysis, of two-independent cohorts of colon cancer patients, revealed that low expression of SETD3 is a reliable predictor of poor survival in these patients, which correlates with our findings. Together, our data uncover a new role of the PKMT SETD3 in the regulation of p53-dependent activation of apoptosis in response to DNA damage.

The protein methyltransferase SETD6 is a key regulator of proliferation and inflammatory processes. However, the role of SETD6 in the regulation of additional cell signaling pathways has not been well studied. Here we show that SETD6 is a negative regulator of the oxidative stress response. Depletion of SETD6 from cells results in elevated Nrf2 levels and a significant increase in Nrf2 antioxidant target gene expression. Using proteomic tools, we uncovered a novel interaction between SETD6 and the oxidative stress sensor DJ1, a protein required for Nrf2-dependent transcription of antioxidant target genes. We show that SETD6 binds DJ1 both in-vitro and in cells but does not methylate DJ1. Under basal conditions, SETD6 and DJ1 are associated at chromatin. Through this interaction, SETD6 inhibits DJ1 activity, which in turn leads to the repression of Nrf2-dependent transcription. In response to oxidative stress, the transcription of Nrf2 antioxidant genes increases. We here show that under this condition, SETD6 mRNA and protein levels are reduced, leading to elevation in Nrf2 expression level and to a weaken interaction between SETD6 and DJ1 at chromatin. Taken together, these findings demonstrate that SETD6 negatively regulates the Nrf2-mediated oxidative stress response through a physical and catalytically independent interaction with DJ1 at chromatin.

Although heritable factors are an important determinant of risk of early-onset cancer, the majority of these malignancies appear to occur sporadically without identifiable risk factors. Germline de novo copy-number variations (CNVs) have been observed in sporadic neurocognitive and cardiovascular disorders. We explored this mechanism in 382 genomes of 116 early-onset cancer case-parent trios and unaffected siblings. Unique de novo germline CNVs were not observed in 107 breast or colon cancer trios or controls but were indeed found in 7% of 43 testicular germ cell tumor trios; this percentage exceeds background CNV rates and suggests a rare de novo genetic paradigm for susceptibility to some human malignancies.

Lysine methylation is a reversible post-translational modification that affects protein function. Lysine methylation is involved in regulating the function of both histone and non-histone proteins, thereby influencing both cellular transcription and the activation of signaling pathways. To date, only a few lysine methyltransferases have been studied in depth. Here, we study the Drosophila homolog of the human lysine methyltransferase SETD3, CG32732/dSETD3. Since mammalian SETD3 is involved in cell proliferation, we tested the effect of dSETD3 on proliferation and growth of Drosophila S2 cells and whole flies. Knockdown of dSETD3 did not alter mTORC1 activity nor proliferation rate of S2 cells. Complete knock-out of dSETD3 in Drosophila flies did not affect their weight, growth rate or fertility. dSETD3 KO flies showed normal responses to starvation and hypoxia. In sum, we could not identify any clear phenotypes for SETD3 knockout animals, indicating that additional work will be required to elucidate the molecular and physiological function of this highly conserved enzyme.

Cortical spreading depression (CSD) is believed to promote migraine headache by enhancing the activity and mechanosensitivity of trigeminal intracranial meningeal afferents. One putative mechanism underlying this afferent response involves an acute excitation of meningeal afferents by cortical efflux of K+ and the ensuing antidromic release of proinflammatory sensory neuropeptides, such as calcitonin gene-related peptide (CGRP).