A network-based method for associating genes with autism spectrum disorder.

Autism spectrum disorder (ASD) is a highly heritable complex disease that affects 1% of the population, yet its underlying molecular mechanisms are largely unknown. Here we study the problem of predicting causal genes for ASD by combining genome-scale data with a network propagation approach. We construct a predictor that integrates multiple omic data sets that assess genomic, transcriptomic, proteomic, and phosphoproteomic associations with ASD. In cross validation our predictor yields mean area under the ROC curve of 0.87 and area under the precision-recall curve of 0.89. We further show that it outperforms previous gene-level predictors of autism association. Finally, we show that we can use the model to predict genes associated with Schizophrenia which is known to share genetic components with ASD.

Human histone H1 variants impact splicing outcome by controlling RNA polymerase II elongation.

Histones shape chromatin structure and the epigenetic landscape. H1, the most diverse histone in the human genome, has 11 variants. Due to the high structural similarity between the H1s, their unique functions in transferring information from the chromatin to mRNA-processing machineries have remained elusive. Here, we generated human cell lines lacking up to five H1 subtypes, allowing us to characterize the genomic binding profiles of six H1 variants. Most H1s bind to specific sites, and binding depends on multiple factors, including GC content. The highly expressed H1.2 has a high affinity for exons, whereas H1.3 binds intronic sequences. H1s are major splicing regulators, especially of exon skipping and intron retention events, through their effects on the elongation of RNA polymerase II (RNAPII). Thus, H1 variants determine splicing fate by modulating RNAPII elongation.

Ataxia Telangiectasia Mutated Signaling Delays Skin Pigmentation upon UV Exposure by Mediating MITF Function toward DNA Repair Mode.

Skin pigmentation is paused after sun exposure; however, the mechanism behind this pausing is unknown. In this study, we found that the UVB-induced DNA repair system, led by the ataxia telangiectasia mutated (ATM) protein kinase, represses MITF transcriptional activity of pigmentation genes while placing MITF in DNA repair mode, thus directly inhibiting pigment production. Phosphoproteomics analysis revealed ATM to be the most significantly enriched pathway among all UVB-induced DNA repair systems. ATM inhibition in mouse or human skin, either genetically or chemically, induces pigmentation. Upon UVB exposure, MITF transcriptional activation is blocked owing to ATM-dependent phosphorylation of MITF on S414, which modifies MITF activity and interactome toward DNA repair, including binding to TRIM28 and RBBP4. Accordingly, MITF genome occupancy is enriched in sites of high DNA damage that are likely repaired. This suggests that ATM harnesses the pigmentation key activator for the necessary rapid, efficient DNA repair, thus optimizing the chances of the cell surviving. Data are available from ProteomeXchange with the identifier PXD041121.

Food-seeking behavior is triggered by skin ultraviolet exposure in males.

Sexual dimorphisms are responsible for profound metabolic differences in health and behavior. Whether males and females react differently to environmental cues, such as solar ultraviolet (UV) exposure, is unknown. Here we show that solar exposure induces food-seeking behavior, food intake, and food-seeking behavior and food intake in men, but not in women, through epidemiological evidence of approximately 3,000 individuals throughout the year. In mice, UVB exposure leads to increased food-seeking behavior, food intake and weight gain, with a sexual dimorphism towards males. In both mice and human males, increased appetite is correlated with elevated levels of circulating ghrelin. Specifically, UVB irradiation leads to p53 transcriptional activation of ghrelin in skin adipocytes, while a conditional p53-knockout in mice abolishes UVB-induced ghrelin expression and food-seeking behavior. In females, estrogen interferes with the p53-chromatin interaction on the ghrelin promoter, thus blocking ghrelin and food-seeking behavior in response to UVB exposure. These results identify the skin as a major mediator of energy homeostasis and may lead to therapeutic opportunities for sex-based treatments of endocrine-related diseases.

Bi-allelic variants in neuronal cell adhesion molecule cause a neurodevelopmental disorder characterized by developmental delay, hypotonia, neuropathy/spasticity.

Cell adhesion molecules are membrane-bound proteins predominantly expressed in the central nervous system along principal axonal pathways with key roles in nervous system development, neural cell differentiation and migration, axonal growth and guidance, myelination, and synapse formation. Here, we describe ten affected individuals with bi-allelic variants in the neuronal cell adhesion molecule NRCAM that lead to a neurodevelopmental syndrome of varying severity; the individuals are from eight families. This syndrome is characterized by developmental delay/intellectual disability, hypotonia, peripheral neuropathy, and/or spasticity. Computational analyses of NRCAM variants, many of which cluster in the third fibronectin type III (Fn-III) domain, strongly suggest a deleterious effect on NRCAM structure and function, including possible disruption of its interactions with other proteins. These findings are corroborated by previous in vitro studies of murine Nrcam-deficient cells, revealing abnormal neurite outgrowth, synaptogenesis, and formation of nodes of Ranvier on myelinated axons. Our studies on zebrafish nrcamaΔ mutants lacking the third Fn-III domain revealed that mutant larvae displayed significantly altered swimming behavior compared to wild-type larvae (p < 0.03). Moreover, nrcamaΔ mutants displayed a trend toward increased amounts of α-tubulin fibers in the dorsal telencephalon, demonstrating an alteration in white matter tracts and projections. Taken together, our study provides evidence that NRCAM disruption causes a variable form of a neurodevelopmental disorder and broadens the knowledge on the growing role of the cell adhesion molecule family in the nervous system.

Gene architecture directs splicing outcome in separate nuclear spatial regions.

How the splicing machinery defines exons or introns as the spliced unit has remained a puzzle for 30 years. Here, we demonstrate that peripheral and central regions of the nucleus harbor genes with two distinct exon-intron GC content architectures that differ in the splicing outcome. Genes with low GC content exons, flanked by long introns with lower GC content, are localized in the periphery, and the exons are defined as the spliced unit. Alternative splicing of these genes results in exon skipping. In contrast, the nuclear center contains genes with a high GC content in the exons and short flanking introns. Most splicing of these genes occurs via intron definition, and aberrant splicing leads to intron retention. We demonstrate that the nuclear periphery and center generate different environments for the regulation of alternative splicing and that two sets of splicing factors form discrete regulatory subnetworks for the two gene architectures. Our study connects 3D genome organization and splicing, thus demonstrating that exon and intron definition modes of splicing occur in different nuclear regions.

The upstream 5' splice site remains associated to the transcription machinery during intron synthesis.

In the earliest step of spliceosome assembly, the two splice sites flanking an intron are brought into proximity by U1 snRNP and U2AF along with other proteins. The mechanism that facilitates this intron looping is poorly understood. Using a CRISPR interference-based approach to halt RNA polymerase II transcription in the middle of introns in human cells, we discovered that the nascent 5' splice site base pairs with a U1 snRNA that is tethered to RNA polymerase II during intron synthesis. This association functionally corresponds with splicing outcome, involves bona fide 5' splice sites and cryptic intronic sites, and occurs transcriptome-wide. Overall, our findings reveal that the upstream 5' splice sites remain attached to the transcriptional machinery during intron synthesis and are thus brought into proximity of the 3' splice sites; potentially mediating the rapid splicing of long introns.

DNA methylation directs microRNA biogenesis in mammalian cells.

MicroRNA (miRNA) biogenesis initiates co-transcriptionally, but how the Microprocessor machinery pinpoints the locations of short precursor miRNA sequences within long flanking regions of the transcript is not known. Here we show that miRNA biogenesis depends on DNA methylation. When the regions flanking the miRNA coding sequence are highly methylated, the miRNAs are more highly expressed, have greater sequence conservation, and are more likely to drive cancer-related phenotypes than miRNAs encoded by unmethylated loci. We show that the removal of DNA methylation from miRNA loci leads to their downregulation. Further, we found that MeCP2 binding to methylated miRNA loci halts RNA polymerase II elongation, leading to enhanced processing of the primary miRNA by Drosha. Taken together, our data reveal that DNA methylation directly affects miRNA biogenesis.

Histone H1.5 binds over splice sites in chromatin and regulates alternative splicing.

Chromatin organization and epigenetic markers influence splicing, though the magnitudes of these effects and the mechanisms are largely unknown. Here, we demonstrate that linker histone H1.5 influences mRNA splicing. We observed that linker histone H1.5 binds DNA over splice sites of short exons in human lung fibroblasts (IMR90 cells). We found that association of H1.5 with these splice sites correlated with the level of inclusion of alternatively spliced exons. Exons marked by H1.5 had more RNA polymerase II (RNAP II) stalling near the 3' splice site than did exons not associated with H1.5. In cells depleted of H1.5, we showed that the inclusion of five exons evaluated decreased and that RNAP II levels over these exons were also reduced. Our findings indicate that H1.5 is involved in regulation of splice site selection and alternative splicing, a function not previously demonstrated for linker histones.

Combinatorial treatment increases IKAP levels in human cells generated from Familial Dysautonomia patients.

Familial Dysautonomia (FD) is an autosomal recessive congenital neuropathy that results from a point mutation at the 5' splice site of intron 20 in the IKBKAP gene. This mutation decreases production of the IKAP protein, and treatments that increase the level of the full-length IKBKAP transcript are likely to be of therapeutic value. We previously found that phosphatidylserine (PS), an FDA-approved food supplement, elevates IKAP levels in cells generated from FD patients. Here we demonstrate that combined treatment of cells generated from FD patients with PS and kinetin or PS and the histone deacetylase inhibitor trichostatin A (TSA) resulted in an additive elevation of IKAP compared to each drug alone. This indicates that the compounds influence different pathways. We also found that pridopidine enhances production of IKAP in cells generated from FD patients. Pridopidine has an additive effect on IKAP levels when used in combination with kinetin or TSA, but not with PS; suggesting that PS and pridopidine influence IKBKAP levels through the same mechanism. Indeed, we demonstrate that the effect of PS and pridopidine is through sigma-1 receptor-mediated activation of the BDNF signaling pathway. A combination treatment with any of these drugs with different mechanisms has potential to benefit FD patients.

Genetics of hearing loss in the Arab population of Northern Israel.

For multiple generations, much of the Arab population of Northern Israel has lived in communities with consanguineous marriages and large families. These communities have been particularly cooperative and informative for understanding the genetics of recessive traits. We studied the genetics of hearing loss in this population, evaluating 168 families from 46 different villages. All families were screened for founder variants by Sanger sequencing and 13 families were further evaluated by sequencing all known genes for hearing loss using our targeted gene panel HEar-Seq. Deafness in 34 of 168 families (20%) was explained by founder variants in GJB2, SLC26A4, or OTOF. In 6 of 13 families (46%) evaluated using HEar-Seq, deafness was explained by damaging alleles of SLC26A4, MYO15A, OTOG, LOXHD1, and TBC1D24. In some genes critical to hearing, it is particularly difficult to interpret variants that might affect splicing, because the genes are not expressed in accessible tissue. To address this problem for possible splice-altering variants of MYO15A, we evaluated minigenes transfected into HEK293 cells. Results revealed exon skipping in the message of MYO15A c.9083+6T>A, and intron retention in the message of MYO15A c.8340G>A, in each case leading to a premature stop and consistent with co-segregation of homozygosity for each variant with hearing loss. The profile of genetics of hearing loss in this population reflects the genetic heterogeneity of hearing loss and the usefulness of synthetic technologies to evaluate potentially causal variants in genes not expressed in accessible tissues.

The importance of DNA methylation of exons on alternative splicing.

Alternative splicing (AS) contributes to proteome diversity. As splicing occurs cotranscriptionally, epigenetic determinants such as DNA methylation likely play a part in regulation of AS. Previously, we have shown that DNA methylation marks exons and that a loss of DNA methylation alters splicing patterns in a genome-wide manner. To investigate the influence of DNA methylation on splicing of individual genes, we developed a method to manipulate DNA methylation in vivo in a site-specific manner using the deactivated endonuclease Cas9 fused to enzymes that methylate or demethylate DNA. We used this system to directly change the DNA methylation pattern of selected exons and introns. We demonstrated that changes in the methylation pattern of alternatively spliced exons, but not constitutively spliced exons or introns, altered inclusion levels. This is the first direct demonstration that DNA methylation of exon-encoding regions is directly involved in regulation of AS.

Quantitative mass spectrometry analysis reveals a panel of nine proteins as diagnostic markers for colon adenocarcinomas.

Adenocarcinomas are cancers originating from the gland forming cells of the colon and rectal lining, and are known to be the most common type of colorectal cancers. The current diagnosis strategies for colorectal cancers include biopsy, laboratory tests, and colonoscopy which are time consuming. Identification of protein biomarkers could aid in the detection of colon adenocarcinomas (CACs). In this study, tissue proteome of colon adenocarcinomas (n = 11) was compared with the matched control specimens (n = 11) using isobaric tags for relative and absolute quantitation (iTRAQ) based liquid chromatography-mass spectrometry (LC-MS/MS) approach. A list of 285 significantly altered proteins was identified in colon adenocarcinomas as compared to its matched controls, which are associated with growth and malignancy of the tumors. Protein interaction analysis revealed the association of altered proteins in colon adenocarcinomas with various transcription factors and their targets. A panel of nine proteins was validated using multiple reaction monitoring (MRM). Additionally, S100A9 was also validated using immunoblotting. The identified panel of proteins may serve as potential biomarkers and thereby aid in the detection of colon adenocarcinomas.

Phosphatidylserine improves axonal transport by inhibition of HDAC and has potential in treatment of neurodegenerative diseases.

Familial dysautonomia (FD) is a rare children neurodegenerative disease caused due to a point mutation in the IKBKAP gene that results in decreased IKK complex-associated protein (IKAP) protein production. The disease affects mostly the dorsal root ganglion (DRG) and the sympathetic ganglion. Recently, we found that the molecular mechanisms underlying neurodegeneration in FD patients are defects in axonal transport of nerve growth factors and microtubule stability in the DRG. Neurons are highly polarized cells with very long axons. In order to survive and maintain proper function, neurons depend on transport of proteins and other cellular components from the neuronal body along the axons. We further demonstrated that IKAP is necessary for axon maintenance and showed that phosphatidylserine acts as an HDAC6 inhibitor to rescue neuronal function in FD cells. In this review, we will highlight our latest research findings.

Calpain 12 Function Revealed through the Study of an Atypical Case of Autosomal Recessive Congenital Ichthyosis.

Congenital erythroderma is a rare and often life-threatening condition, which has been shown to result from mutations in several genes encoding important components of the epidermal differentiation program. Using whole exome sequencing, we identified in a child with congenital exfoliative erythroderma, hypotrichosis, severe nail dystrophy and failure to thrive, two heterozygous mutations in ABCA12 (c.2956C>T, p.R986W; c.5778+2T>C, p. G1900Mfs*16), a gene known to be associated with two forms of ichthyosis, autosomal recessive congenital ichthyosis, and harlequin ichthyosis. Because the patient displayed an atypical phenotype, including severe hair and nail manifestations, we scrutinized the exome sequencing data for additional potentially deleterious genetic variations in genes of relevance to the cornification process. Two mutations were identified in CAPN12, encoding a member of the calpain proteases: a paternal missense mutation (c.1511C>A; p.P504Q) and a maternal deletion due to activation of a cryptic splice site in exon 9 of the gene (c.1090_1129del; p.Val364Lysfs*11). The calpain 12 protein was found to be expressed in both the epidermis and hair follicle of normal skin, but its expression was dramatically reduced in the patient's skin. The downregulation of capn12 expression in zebrafish was associated with abnormal epidermal morphogenesis. Small interfering RNA knockdown of CAPN12 in three-dimensional human skin models was associated with acanthosis, disorganized epidermal architecture, and downregulation of several differentiation markers, including filaggrin. Accordingly, filaggrin expression was almost absent in the patient skin. Using ex vivo live imaging, small interfering RNA knockdown of calpain 12 in skin from K14-H2B GFP mice led to significant hair follicle catagen transformation compared with controls. In summary, our results indicate that calpain 12 plays an essential role during epidermal ontogenesis and normal hair follicle cycling and that its absence may aggravate the clinical manifestations of ABCA12 mutations.

Phosphatidylserine Ameliorates Neurodegenerative Symptoms and Enhances Axonal Transport in a Mouse Model of Familial Dysautonomia.

Familial Dysautonomia (FD) is a neurodegenerative disease in which aberrant tissue-specific splicing of IKBKAP exon 20 leads to reduction of IKAP protein levels in neuronal tissues. Here we generated a conditional knockout (CKO) mouse in which exon 20 of IKBKAP is deleted in the nervous system. The CKO FD mice exhibit developmental delays, sensory abnormalities, and less organized dorsal root ganglia (DRGs) with attenuated axons compared to wild-type mice. Furthermore, the CKO FD DRGs show elevated HDAC6 levels, reduced acetylated α-tubulin, unstable microtubules, and impairment of axonal retrograde transport of nerve growth factor (NGF). These abnormalities in DRG properties underlie neuronal degeneration and FD symptoms. Phosphatidylserine treatment decreased HDAC6 levels and thus increased acetylation of α-tubulin. Further PS treatment resulted in recovery of axonal outgrowth and enhanced retrograde axonal transport by decreasing histone deacetylase 6 (HDAC6) levels and thus increasing acetylation of α-tubulin levels. Thus, we have identified the molecular pathway that leads to neurodegeneration in FD and have demonstrated that phosphatidylserine treatment has the potential to slow progression of neurodegeneration.

How Are Short Exons Flanked by Long Introns Defined and Committed to Splicing?

The splice sites (SSs) delimiting an intron are brought together in the earliest step of spliceosome assembly yet it remains obscure how SS pairing occurs, especially when introns are thousands of nucleotides long. Splicing occurs in vivo in mammals within minutes regardless of intron length, implying that SS pairing can instantly follow transcription. Also, factors required for SS pairing, such as the U1 small nuclear ribonucleoprotein (snRNP) and U2AF65, associate with RNA polymerase II (RNAPII), while nucleosomes preferentially bind exonic sequences and associate with U2 snRNP. Based on recent publications, we assume that the 5' SS-bound U1 snRNP can remain tethered to RNAPII until complete synthesis of the downstream intron and exon. An additional U1 snRNP then binds the downstream 5' SS, whereas the RNAPII-associated U2AF65 binds the upstream 3' SS to facilitate SS pairing along with exon definition. Next, the nucleosome-associated U2 snRNP binds the branch site to advance splicing complex assembly. This may explain how RNAPII and chromatin are involved in spliceosome assembly and how introns lengthened during evolution with a relatively minimal compromise in splicing.

A network-based analysis of colon cancer splicing changes reveals a tumorigenesis-favoring regulatory pathway emanating from ELK1.

Splicing aberrations are prominent drivers of cancer, yet the regulatory pathways controlling them are mostly unknown. Here we develop a method that integrates physical interaction, gene expression, and alternative splicing data to construct the largest map of transcriptomic and proteomic interactions leading to cancerous splicing aberrations defined to date, and identify driver pathways therein. We apply our method to colon adenocarcinoma and non-small-cell lung carcinoma. By focusing on colon cancer, we reveal a novel tumor-favoring regulatory pathway involving the induction of the transcription factor MYC by the transcription factor ELK1, as well as the subsequent induction of the alternative splicing factor PTBP1 by both. We show that PTBP1 promotes specific RAC1,NUMB, and PKM splicing isoforms that are major triggers of colon tumorigenesis. By testing the pathway's activity in patient tumor samples, we find ELK1,MYC, and PTBP1 to be overexpressed in conjunction with oncogenic KRAS mutations, and show that these mutations increase ELK1 levels via the RAS-MAPK pathway. We thus illuminate, for the first time, a full regulatory pathway connecting prevalent cancerous mutations to functional tumor-inducing splicing aberrations. Our results demonstrate our method is applicable to different cancers to reveal regulatory pathways promoting splicing aberrations.

Phosphatidylserine enhances IKBKAP transcription by activating the MAPK/ERK signaling pathway.

Familial dysautonomia (FD) is a genetic disorder manifested due to abnormal development and progressive degeneration of the sensory and autonomic nervous system. FD is caused by a point mutation in the IKBKAP gene encoding the IKAP protein, resulting in decreased protein levels. A promising potential treatment for FD is phosphatidylserine (PS); however, the manner by which PS elevates IKAP levels has yet to be identified. Analysis of ChIP-seq results of the IKBKAP promoter region revealed binding of the transcription factors CREB and ELK1, which are regulated by the mitogen-activated protein kinase (MAPK)/extracellular-regulated kinase (ERK) signaling pathway. We show that PS treatment enhanced ERK phosphorylation in cells derived from FD patients. ERK activation resulted in elevated IKBKAP transcription and IKAP protein levels, whereas pretreatment with the MAPK inhibitor U0126 blocked elevation of the IKAP protein level. Overexpression of either ELK1 or CREB activated the IKBKAP promoter, whereas downregulation of these transcription factors resulted in a decrease of the IKAP protein. Additionally, we show that PS improves cell migration, known to be enhanced by MAPK/ERK activation and abrogated in FD cells. In conclusion, our results demonstrate that PS activates the MAPK/ERK signaling pathway, resulting in activation of transcription factors that bind the promoter region of IKBKAP and thus enhancing its transcription. Therefore, compounds that activate the MAPK/ERK signaling pathway could constitute potential treatments for FD.

Regulation of alternative splicing through coupling with transcription and chromatin structure.

Alternative precursor messenger RNA (pre-mRNA) splicing plays a pivotal role in the flow of genetic information from DNA to proteins by expanding the coding capacity of genomes. Regulation of alternative splicing is as important as regulation of transcription to determine cell- and tissue-specific features, normal cell functioning, and responses of eukaryotic cells to external cues. Its importance is confirmed by the evolutionary conservation and diversification of alternative splicing and the fact that its deregulation causes hereditary disease and cancer. This review discusses the multiple layers of cotranscriptional regulation of alternative splicing in which chromatin structure, DNA methylation, histone marks, and nucleosome positioning play a fundamental role in providing a dynamic scaffold for interactions between the splicing and transcription machineries. We focus on evidence for how the kinetics of RNA polymerase II (RNAPII) elongation and the recruitment of splicing factors and adaptor proteins to chromatin components act in coordination to regulate alternative splicing.