Rare inherited A2BP1 deletion in a proband with autism and developmental hemiparesis

To RNA-binding and beyond: Emerging facets of the role of Rbfox proteins in development and disease

The RNA-binding Fox-1 homologue (Rbfox) proteins represent an ancient family of splicing factors, conserved through evolution. All members share an RNA recognition motif (RRM), and a particular affinity for the GCAUG signature in target RNA molecules. The role of Rbfox, as a splice factor, deciding the tissue-specific inclusion/exclusion of an exon, depending on its binding position on the flanking introns, is well known. Rbfox often acts in concert with other splicing factors, and forms splicing regulatory networks. Apart from this canonical role, recent studies show that Rbfox can also function as a transcription co-factor, and affects mRNA stability and translation. The repertoire of Rbfox targets is vast, including genes involved in the development of tissue lineages, such as neurogenesis, myogenesis, and erythropoeiesis, and molecular processes, including cytoskeletal dynamics, and calcium handling. A second layer of complexity is added by the fact that Rbfox expression itself is regulated by multiple mechanisms, and, in vertebrates, exhibits tissue-specific expression. The optimum dosage of Rbfox is critical, and its misexpression is etiological to various disease conditions. In this review, we discuss the contextual roles played by Rbfox as a tissue-specific regulator for the expression of many important genes with diverse functions, through the lens of the emerging data which highlights its involvement in many human diseases. Furthermore, we explore the mechanistic details provided by studies in model organisms, with emphasis on the work with Drosophila. This article is categorized under: RNA Processing > Splicing Mechanisms RNA Interactions with Proteins and Other Molecules > Protein-RNA Interactions: Functional Implications RNA Turnover and Surveillance > Regulation of RNA Stability RNA Processing > Splicing Regulation/Alternative Splicing.

Influence of Age on Skeletal Muscle Hypertrophy and Atrophy Signaling: Established Paradigms and Unexpected Links

Skeletal muscle atrophy in an inevitable occurrence with advancing age, and a consequence of disease including cancer. Muscle atrophy in the elderly is managed by a regimen of resistance exercise and increased protein intake. Understanding the signaling that regulates muscle mass may identify potential therapeutic targets for the prevention and reversal of muscle atrophy in metabolic and neuromuscular diseases. This review covers the major anabolic and catabolic pathways that regulate skeletal muscle mass, with a focus on recent progress and potential new players.

De Novo Frameshift Variants in the Neuronal Splicing Factor NOVA2 Result in a Common C-Terminal Extension and Cause a Severe Form of Neurodevelopmental Disorder

The neuro-oncological ventral antigen 2 (NOVA2) protein is a major factor regulating neuron-specific alternative splicing (AS), previously associated with an acquired neurologic condition, the paraneoplastic opsoclonus-myoclonus ataxia (POMA). We report here six individuals with de novo frameshift variants in NOVA2 affected with a severe neurodevelopmental disorder characterized by intellectual disability (ID), motor and speech delay, autistic features, hypotonia, feeding difficulties, spasticity or ataxic gait, and abnormal brain MRI. The six variants lead to the same reading frame, adding a common proline rich C-terminal part instead of the last KH RNA binding domain. We detected 41 genes differentially spliced after NOVA2 downregulation in human neural cells. The NOVA2 variant protein shows decreased ability to bind target RNA sequences and to regulate target AS events. It also fails to complement the effect on neurite outgrowth induced by NOVA2 downregulation in vitro and to rescue alterations of retinotectal axonal pathfinding induced by loss of NOVA2 ortholog in zebrafish. Our results suggest a partial loss-of-function mechanism rather than a full heterozygous loss-of-function, although a specific contribution of the novel C-terminal extension cannot be excluded.

Copy number variations associated with fetal congenital kidney malformations

Background

Congenital anomalies of the kidney and urinary tract (CAKUT) constitute 20–30% of all congenital malformations. Within the CAKUT phenotypic spectrum, renal hypodysplasia (RHD) is particularly severe. This study aimed to evaluate the applicability of single-nucleotide polymorphism (SNP) array test in prenatal diagnosis of RHD for improving prenatal genetic counseling and to search for evidence of a possible causative role of copy-number variations (CNVs) in RHD.

Results

We performed a systematic survey of CNV burden in 120 fetuses with RHD: 103 cases were isolated RHD and 17 were non-isolated RHD. Single-nucleotide polymorphism (SNP) array test was performed using the Affymetrix CytoScan HD platform. All annotated CNVs were validated by fluorescence in situ hybridization. We identified abnormal CNVs in 15 (12.5%) cases of RHD; of these CNVs, 11 were pathogenic and 4 were variants of uncertain significance. The detection rate of abnormal CNVs in non-isolated RHD was higher (29.4%, 5/17) than that in isolated RHD (9.7%, 10/103) (P = 0.060). Parents are more inclined to terminate the pregnancy if the fetuses have pathogenic results of the SNP-array test.

Conclusions

The variable phenotypes that abnormal CNVs may cause indicate the genetic counseling is needed for RHD cases.

RBFOX1, encoding a splicing regulator, is a candidate gene for aggressive behavior

The RBFOX1 gene (or A2BP1) encodes a splicing factor important for neuronal development that has been related to autism spectrum disorder and other neurodevelopmental phenotypes. Evidence from complementary sources suggests that this gene contributes to aggressive behavior. Suggestive associations with RBFOX1 have been identified in genome-wide association studies (GWAS) of anger, conduct disorder, and aggressive behavior. Nominal association signals in RBFOX1 were also found in an epigenome-wide association study (EWAS) of aggressive behavior. Also, variants in this gene affect temporal lobe volume, a brain area that is altered in several aggression-related phenotypes. In animals, this gene has been shown to modulate aggressive behavior in Drosophila. RBFOX1 has also been associated with canine aggression and is upregulated in mice that show increased aggression after frustration of an expected reward. Associated common genetic variants as well as rare duplications and deletions affecting RBFOX1 have been identified in several psychiatric and neurodevelopmental disorders that are often comorbid with aggressive behaviors. In this paper, we comprehensively review the cumulative evidence linking RBFOX1 to aggression behavior and provide new results implicating RBFOX1 in this phenotype. Most of these studies (genetic and epigenetic analyses in humans, neuroimaging genetics, gene expression and animal models) are hypothesis-free, which strengthens the validity of the findings, although all the evidence is nominal and should therefore be taken with caution. Further studies are required to clarify in detail the role of this gene in this complex phenotype.

Genome-wide association study identifies RBFOX1 locus influencing brain glucose metabolism

BACKGROUND

Fluorodeoxyglucose f18 positron emission tomography (18F-FDG PET) is regarded as the only functional neuroimaging biomarker for degeneration which can be used to increase the certainty of Alzheimer's disease (AD) pathophysiological process in research settings or as an optional clinical tool where available. Although a decline in FDG metabolism was confirmed in some regions known to be associated with AD, there was little known about the genetic association of FDG metabolism in AD cohorts. In this study, we present the first genome-wide association study (GWAS) analysis of brain FDG metabolism.

METHODS

A total of 222 individuals were included from the Alzheimer's Disease Neuroimaging Initiative 1 (ADNI-1) cohort. All subjects were restricted to non-Hispanic Caucasians and met all quality control (QC) criteria. Associations of 18F-FDG with the genetic variants were assessed using PLINK 1.07 under the additive genetic model. Genome-wide associations were visualized using a software program R 3.2.3.

RESULTS

One significant SNP rs12444565 in RNA-binding Fox1 (RBFOX1) was found to have a strong association with 18F-FDG (P=6.06×10-8). Rs235141, rs79037, rs12526331 and rs12529764 were identified as four suggestive loci associated with 18F-FDG.

CONCLUSIONS

Our study results suggest that a genome-wide significant SNP (rs12444565) in the RBFOX1, and four suggestive loci (rs235141, rs79037, rs12526331 and rs12529764) are associated with 18F-FDG.

rbFOX1/MBNL1 competition for CCUG RNA repeats binding contributes to myotonic dystrophy type 1/type 2 differences

Myotonic dystrophy type 1 and type 2 (DM1, DM2) are caused by expansions of CTG and CCTG repeats, respectively. RNAs containing expanded CUG or CCUG repeats interfere with the metabolism of other RNAs through titration of the Muscleblind-like (MBNL) RNA binding proteins. DM2 follows a more favorable clinical course than DM1, suggesting that specific modifiers may modulate DM severity. Here, we report that the rbFOX1 RNA binding protein binds to expanded CCUG RNA repeats, but not to expanded CUG RNA repeats. Interestingly, rbFOX1 competes with MBNL1 for binding to CCUG expanded repeats and overexpression of rbFOX1 partly releases MBNL1 from sequestration within CCUG RNA foci in DM2 muscle cells. Furthermore, expression of rbFOX1 corrects alternative splicing alterations and rescues muscle atrophy, climbing and flying defects caused by expression of expanded CCUG repeats in a Drosophila model of DM2.

Genomics of autism spectrum disorder: approach to therapy

Autism spectrum disorder (ASD) is a highly prevalent neurodevelopmental condition with no current treatment available. Although advances in genetics and genomics have identified hundreds of genes associated with ASD, very little is known about the pathophysiology of ASD and the functional contribution of specific genes to ASD phenotypes. Improved understanding of the biological function of ASD-associated genes and how this heterogeneous group of genetic variants leads to the disease is needed in order to develop therapeutic strategies. Here, we review the current state of ASD research related to gene discovery and examples of emerging molecular mechanisms (protein translation and alternative splicing). In addition, we discuss how patient-derived three-dimensional brain organoids might provide an opportunity to model specific genetic variants in order to define molecular and cellular defects that could be amenable for developing and screening personalized therapies related to ASD.

miRNA targeting and alternative splicing in the stress response - events hosted by membrane-less compartments

Stress can be temporary or chronic, and mild or acute. Depending on its extent and severity, cells either alter their metabolism, and adopt a new state, or die. Fluctuations in environmental conditions occur frequently, and such stress disturbs cellular homeostasis, but in general, stresses are reversible and last only a short time. There is increasing evidence that regulation of gene expression in response to temporal stress happens post-transcriptionally in specialized subcellular membrane-less compartments called ribonucleoprotein (RNP) granules. RNP granules assemble through a concentration-dependent liquid-liquid phase separation of RNA-binding proteins that contain low-complexity sequence domains (LCDs). Interestingly, many factors that regulate microRNA (miRNA) biogenesis and alternative splicing are RNA-binding proteins that contain LCDs and localize to stress-induced liquid-like compartments. Consequently, gene silencing through miRNAs and alternative splicing of pre-mRNAs are emerging as crucial post-transcriptional mechanisms that function on a genome-wide scale to regulate the cellular stress response. In this Review, we describe the interplay between these two post-transcriptional processes that occur in liquid-like compartments as an adaptive cellular response to stress.

Misregulation of an Activity-Dependent Splicing Network as a Common Mechanism Underlying Autism Spectrum Disorders

A key challenge in understanding and ultimately treating autism is to identify common molecular mechanisms underlying this genetically heterogeneous disorder. Transcriptomic profiling of autistic brains has revealed correlated misregulation of the neuronal splicing regulator nSR100/SRRM4 and its target microexon splicing program in more than one-third of analyzed individuals. To investigate whether nSR100 misregulation is causally linked to autism, we generated mutant mice with reduced levels of this protein and its target splicing program. Remarkably, these mice display multiple autistic-like features, including altered social behaviors, synaptic density, and signaling. Moreover, increased neuronal activity, which is often associated with autism, results in a rapid decrease in nSR100 and splicing of microexons that significantly overlap those misregulated in autistic brains. Collectively, our results provide evidence that misregulation of an nSR100-dependent splicing network controlled by changes in neuronal activity is causally linked to a substantial fraction of autism cases.

Leveraging blood serotonin as an endophenotype to identify de novo and rare variants involved in autism

BACKGROUND

Autism spectrum disorder (ASD) is one of the most highly heritable neuropsychiatric disorders, but underlying molecular mechanisms are still unresolved due to extreme locus heterogeneity. Leveraging meaningful endophenotypes or biomarkers may be an effective strategy to reduce heterogeneity to identify novel ASD genes. Numerous lines of evidence suggest a link between hyperserotonemia, i.e., elevated serotonin (5-hydroxytryptamine or 5-HT) in whole blood, and ASD. However, the genetic determinants of blood 5-HT level and their relationship to ASD are largely unknown.

METHODS

In this study, pursuing the hypothesis that de novo variants (DNVs) and rare risk alleles acting in a recessive mode may play an important role in predisposition of hyperserotonemia in people with ASD, we carried out whole exome sequencing (WES) in 116 ASD parent-proband trios with most (107) probands having 5-HT measurements.

RESULTS

Combined with published ASD DNVs, we identified USP15 as having recurrent de novo loss of function mutations and discovered evidence supporting two other known genes with recurrent DNVs (FOXP1 and KDM5B). Genes harboring functional DNVs significantly overlap with functional/disease gene sets known to be involved in ASD etiology, including FMRP targets and synaptic formation and transcriptional regulation genes. We grouped the probands into High-5HT and Normal-5HT groups based on normalized serotonin levels, and used network-based gene set enrichment analysis (NGSEA) to identify novel hyperserotonemia-related ASD genes based on LoF and missense DNVs. We found enrichment in the High-5HT group for a gene network module (DAWN-1) previously implicated in ASD, and this points to the TGF-β pathway and cell junction processes. Through analysis of rare recessively acting variants (RAVs), we also found that rare compound heterozygotes (CHs) in the High-5HT group were enriched for loci in an ASD-associated gene set. Finally, we carried out rare variant group-wise transmission disequilibrium tests (gTDT) and observed significant association of rare variants in genes encoding a subset of the serotonin pathway with ASD.

CONCLUSIONS

Our study identified USP15 as a novel gene implicated in ASD based on recurrent DNVs. It also demonstrates the potential value of 5-HT as an effective endophenotype for gene discovery in ASD, and the effectiveness of this strategy needs to be further explored in studies of larger sample sizes.

Reverse Pathway Genetic Approach Identifies Epistasis in Autism Spectrum Disorders

Although gene-gene interaction, or epistasis, plays a large role in complex traits in model organisms, genome-wide by genome-wide searches for two-way interaction have limited power in human studies. We thus used knowledge of a biological pathway in order to identify a contribution of epistasis to autism spectrum disorders (ASDs) in humans, a reverse-pathway genetic approach. Based on previous observation of increased ASD symptoms in Mendelian disorders of the Ras/MAPK pathway (RASopathies), we showed that common SNPs in RASopathy genes show enrichment for association signal in GWAS (P = 0.02). We then screened genome-wide for interactors with RASopathy gene SNPs and showed strong enrichment in ASD-affected individuals (P < 2.2 x 10-16), with a number of pairwise interactions meeting genome-wide criteria for significance. Finally, we utilized quantitative measures of ASD symptoms in RASopathy-affected individuals to perform modifier mapping via GWAS. One top region overlapped between these independent approaches, and we showed dysregulation of a gene in this region, GPR141, in a RASopathy neural cell line. We thus used orthogonal approaches to provide strong evidence for a contribution of epistasis to ASDs, confirm a role for the Ras/MAPK pathway in idiopathic ASDs, and to identify a convergent candidate gene that may interact with the Ras/MAPK pathway.

Prenatal sonographic diagnosis of urorectal septum malformation sequence and chromosomal microarray analysis: A case report and review of the literature

INTRODUCTION

Urorectal septum malformation sequence (URSMS) is a rare congenital abnormal syndrome that is caused by the incomplete division of the cloaca. Based on whether the cloaca membrane breaks down or not, the URSMS are classified as full and partial forms. The prenatal diagnosis of URSMS remains challenging because of poor recognition to this malformation and the relatively non-specific sonographic features. We report a prenatally sonographic diagnosed case of the partial URSMS, and review the literature to summarize the prenatal features.

CASE REPORT AND REVIEW

A 37-year old woman was referred at 24 weeks of gestation for fetal abdominal cyst. Detailed sonographic examination was done and revealed the vesicocolic fistula, distended colon, absence of perianal hypoechoic ring, pyelectasis, and small stomach bubble. The URSMS was suspected.Amniocentesis was done and karyotyping revealed 46,XY. Furthermore, chromosomal microarray analysis (CMA) was performed for the first time in URSMS and an alteration of 111.8Kb deletion was detected in 16p13.3 which was located inside the RBFOX1 gene. Parental studies showed that the deletion was inherited from the father who has nomal clinical phenotype.The woman elected to terminate the pregnancy at 25 weeks gestation and postmortem examination confirmed the diagnosis of partial URSMS.The published studies were reviewed and 28 cases of URSMS with conducted prenatal ultrasonography were collected in this report. The most common sonographic description, as suspicious signs of URSMS, were severe oligohydramnios or anhydramnios, urinary tract anomalies, fetal intra-abdominal cysts, and dilated bowel. Also, enterolithiasis and vesicocolic fistula were relatively infrequent but highly specific feature of URSMS.

CONCLUSIONS

URSMS is difficult to be diagnosed prenatally. However, it has characteristic features that can be detected by fetal ultrasonography, and a precise prenatal sonographic examination is crucial for diagnosing URSMS. Besides, more genomic profiling studies are needed to elucidate the causality.

Developmental regulation of RNA processing by Rbfox proteins

The Rbfox genes encode an ancient family of sequence-specific RNA binding proteins (RBPs) that are critical developmental regulators in multiple tissues including skeletal muscle, cardiac muscle, and brain. The hallmark of Rbfox proteins is a single high-affinity RRM domain, highly conserved from insects to humans, that binds preferentially to UGCAUG motifs at diverse regulatory sites in pre-mRNA introns, mRNA 3'UTRs, and pre-miRNAs hairpin structures. Versatile regulatory circuits operate on Rbfox pre-mRNA and mRNA to ensure proper expression of Rbfox1 protein isoforms, which then act on the broader transcriptome to regulate alternative splicing networks, mRNA stability and translation, and microRNA processing. Complex Rbfox expression is encoded in large genes encompassing multiple promoters and alternative splicing options that govern spatiotemporal expression of structurally distinct and tissue-specific protein isoforms with different classes of RNA targets. Nuclear Rbfox1 is a candidate master regulator that binds intronic UGCAUG elements to impact splicing efficiency of target alternative exons, many in transcripts for other splicing regulators. Tissue-specificity of Rbfox-mediated alternative splicing is executed by combinatorial regulation through the integrated activity of Rbfox proteins and synergistic or antagonistic splicing factors. Studies in animal models show that Rbfox1-related genes are critical for diverse developmental processes including germ cell differentiation and memory in Drosophila, neuronal migration and function in mouse brain, myoblast fusion and skeletal muscle function, and normal heart function. Finally, genetic and biochemical evidence suggest that aberrations in Rbfox-regulated circuitry are risk factors for multiple human disorders, especially neurodevelopmental disorders including epilepsy and autism, and cardiac hypertrophy. WIREs RNA 2017, 8:e1398. doi: 10.1002/wrna.1398 For further resources related to this article, please visit the WIREs website.

Repression of Pumilio Protein Expression by Rbfox1 Promotes Germ Cell Differentiation

RNA-binding Fox (Rbfox) proteins have well-established roles in regulating alternative splicing, but specific Rbfox isoforms lack nuclear localization signals and accumulate in the cytoplasm. The potential splicing-independent functions of these proteins remain unknown. Here we demonstrate that cytoplasmic Drosophila Rbfox1 regulates germ cell development and represses the translation of mRNAs containing (U)GCAUG elements within their 3'UTRs. During germline cyst differentiation, Rbfox1 targets pumilio mRNA for destabilization and translational silencing, thereby promoting germ cell development. Mis-expression of pumilio results in the formation of germline tumors, which contain cysts that break down and dedifferentiate back to single, mitotically active cells. Together, these results reveal that cytoplasmic Rbfox family members regulate the translation of specific target mRNAs. In the Drosophila ovary, this activity provides a genetic barrier that prevents germ cells from reverting back to an earlier developmental state. The finding that Rbfox proteins regulate mRNA translation has implications for Rbfox-related diseases.

New genes for focal epilepsies with speech and language disorders

The last 2 years have seen exciting advances in the genetics of Landau-Kleffner syndrome and related disorders, encompassed within the epilepsy-aphasia spectrum (EAS). The striking finding of mutations in the N-methyl-D-aspartate (NMDA) receptor subunit gene GRIN2A as the first monogenic cause in up to 20% of patients with EAS suggests that excitatory glutamate receptors play a key role in these disorders. Patients with GRIN2A mutations have a recognizable speech and language phenotype that may assist with diagnosis. Other molecules involved in RNA binding and cell adhesion have been implicated in EAS; copy number variations are also found. The emerging picture highlights the overlap between the genetic determinants of EAS with speech and language disorders, intellectual disability, autism spectrum disorders and more complex developmental phenotypes.

MiR-980 Is a Memory Suppressor MicroRNA that Regulates the Autism-Susceptibility Gene A2bp1

MicroRNAs have been associated with many different biological functions, but little is known about their roles in conditioned behavior. We demonstrate that Drosophila miR-980 is a memory suppressor gene functioning in multiple regions of the adult brain. Memory acquisition and stability were both increased by miR-980 inhibition. Whole cell recordings and functional imaging experiments indicated that miR-980 regulates neuronal excitability. We identified the autism susceptibility gene, A2bp1, as an mRNA target for miR-980. A2bp1 levels varied inversely with miR-980 expression; memory performance was directly related to A2bp1 levels. In addition, A2bp1 knockdown reversed the memory gains produced by miR-980 inhibition, consistent with A2bp1 being a downstream target of miR-980 responsible for the memory phenotypes. Our results indicate that miR-980 represses A2bp1 expression to tune the excitable state of neurons, and the overall state of excitability translates to memory impairment or improvement.

Targeted massively parallel sequencing of autism spectrum disorder-associated genes in a case control cohort reveals rare loss-of-function risk variants

Background

Autism spectrum disorder (ASD) is highly heritable, yet genome-wide association studies (GWAS), copy number variation screens, and candidate gene association studies have found no single factor accounting for a large percentage of genetic risk. ASD trio exome sequencing studies have revealed genes with recurrent de novo loss-of-function variants as strong risk factors, but there are relatively few recurrently affected genes while as many as 1000 genes are predicted to play a role. As such, it is critical to identify the remaining rare and low-frequency variants contributing to ASD.

Methods

We have utilized an approach of prioritization of genes by GWAS and follow-up with massively parallel sequencing in a case-control cohort. Using a previously reported ASD noise reduction GWAS analyses, we prioritized 837 RefSeq genes for custom targeting and sequencing. We sequenced the coding regions of those genes in 2071 ASD cases and 904 controls of European white ancestry. We applied comprehensive annotation to identify single variants which could confer ASD risk and also gene-based association analysis to identify sets of rare variants associated with ASD.

Results

We identified a significant over-representation of rare loss-of-function variants in genes previously associated with ASD, including a de novo premature stop variant in the well-established ASD candidate gene RBFOX1. Furthermore, ASD cases were more likely to have two damaging missense variants in candidate genes than controls. Finally, gene-based rare variant association implicates genes functioning in excitatory neurotransmission and neurite outgrowth and guidance pathways including CACNAD2, KCNH7, and NRXN1.

Conclusions

We find suggestive evidence that rare variants in synaptic genes are associated with ASD and that loss-of-function mutations in ASD candidate genes are a major risk factor, and we implicate damaging mutations in glutamate signaling receptors and neuronal adhesion and guidance molecules. Furthermore, the role of de novo mutations in ASD remains to be fully investigated as we identified the first reported protein-truncating variant in RBFOX1 in ASD. Overall, this work, combined with others in the field, suggests a convergence of genes and molecular pathways underlying ASD etiology.

Electronic supplementary material

The online version of this article (doi:10.1186/s13229-015-0034-z) contains supplementary material, which is available to authorized users.

ZNF804A Transcriptional Networks in Differentiating Neurons Derived from Induced Pluripotent Stem Cells of Human Origin

ZNF804A (Zinc Finger Protein 804A) has been identified as a candidate gene for schizophrenia (SZ), autism spectrum disorders (ASD), and bipolar disorder (BD) in replicated genome wide association studies (GWAS) and by copy number variation (CNV) analysis. Although its function has not been well-characterized, ZNF804A contains a C2H2-type zinc-finger domain, suggesting that it has DNA binding properties, and consequently, a role in regulating gene expression. To further explore the role of ZNF804A on gene expression and its downstream targets, we used a gene knockdown (KD) approach to reduce its expression in neural progenitor cells (NPCs) derived from induced pluripotent stem cells (iPSCs). KD was accomplished by RNA interference (RNAi) using lentiviral particles containing shRNAs that target ZNF804A mRNA. Stable transduced NPC lines were generated after puromycin selection. A control cell line expressing a random (scrambled) shRNA was also generated. Neuronal differentiation was induced, RNA was harvested after 14 days and transcriptome analysis was carried out using RNA-seq. 1815 genes were found to be differentially expressed at a nominally significant level (p<0.05); 809 decreased in expression in the KD samples, while 1106 increased. Of these, 370 achieved genome wide significance (FDR<0.05); 125 were lower in the KD samples, 245 were higher. Pathway analysis showed that genes involved in interferon-signaling were enriched among those that were down-regulated in the KD samples. Correspondingly, ZNF804A KD was found to affect interferon-alpha 2 (IFNA2)-mediated gene expression. The findings suggest that ZNF804A may affect a differentiating neuron’s response to inflammatory cytokines, which is consistent with models of SZ and ASD that support a role for infectious disease, and/or autoimmunity in a subgroup of patients.

The RNA-binding protein Rbfox1 regulates splicing required for skeletal muscle structure and function

The Rbfox family of RNA-binding proteins is highly conserved with established roles in alternative splicing (AS) regulation. High-throughput studies aimed at understanding transcriptome remodeling have revealed skeletal muscle as displaying one of the largest number of AS events. This finding is consistent with requirements for tissue-specific protein isoforms needed to sustain muscle-specific functions. Rbfox1 is abundant in vertebrate brain, heart and skeletal muscle. Genome-wide genetic approaches have linked the Rbfox1 gene to autism, and a brain-specific knockout mouse revealed a critical role for this splicing regulator in neuronal function. Moreover, a Caenorhabditis elegans Rbfox1 homolog regulates muscle-specific splicing. To determine the role of Rbfox1 in muscle function, we developed a conditional knockout mouse model to specifically delete Rbfox1 in adult tissue. We show that Rbfox1 is required for muscle function but a >70% loss of Rbfox1 in satellite cells does not disrupt muscle regeneration. Deep sequencing identified aberrant splicing of multiple genes including those encoding myofibrillar and cytoskeletal proteins, and proteins that regulate calcium handling. Ultrastructure analysis of Rbfox1(-/-) muscle by electron microscopy revealed abundant tubular aggregates. Immunostaining showed mislocalization of the sarcoplasmic reticulum proteins Serca1 and Ryr1 in a pattern indicative of colocalization with the tubular aggregates. Consistent with mislocalization of Serca1 and Ryr1, calcium handling was drastically altered in Rbfox1(-/-) muscle. Moreover, muscle function was significantly impaired in Rbfox1(-/-) muscle as indicated by decreased force generation. These results demonstrate that Rbfox1 regulates a network of AS events required to maintain multiple aspects of muscle physiology.

Cryptic and complex chromosomal aberrations in early-onset neuropsychiatric disorders

Structural variation (SV) is a significant component of the genetic etiology of both neurodevelopmental and psychiatric disorders; however, routine guidelines for clinical genetic screening have been established only in the former category. Genome-wide chromosomal microarray (CMA) can detect genomic imbalances such as copy-number variants (CNVs), but balanced chromosomal abnormalities (BCAs) still require karyotyping for clinical detection. Moreover, submicroscopic BCAs and subarray threshold CNVs are intractable, or cryptic, to both CMA and karyotyping. Here, we performed whole-genome sequencing using large-insert jumping libraries to delineate both cytogenetically visible and cryptic SVs in a single test among 30 clinically referred youth representing a range of severe neuropsychiatric conditions. We detected 96 SVs per person on average that passed filtering criteria above our highest-confidence resolution (6,305 bp) and an additional 111 SVs per genome below this resolution. These SVs rearranged 3.8 Mb of genomic sequence and resulted in 42 putative loss-of-function (LoF) or gain-of-function mutations per person. We estimate that 80% of the LoF variants were cryptic to clinical CMA. We found myriad complex and cryptic rearrangements, including a "paired" duplication (360 kb, 169 kb) that flanks a 5.25 Mb inversion that appears in 7 additional cases from clinical CNV data among 47,562 individuals. Following convergent genomic profiling of these independent clinical CNV data, we interpreted three SVs to be of potential clinical significance. These data indicate that sequence-based delineation of the full SV mutational spectrum warrants exploration in youth referred for neuropsychiatric evaluation and clinical diagnostic SV screening more broadly.

GRASP: analysis of genotype-phenotype results from 1390 genome-wide association studies and corresponding open access database

SUMMARY

We created a deeply extracted and annotated database of genome-wide association studies (GWAS) results. GRASP v1.0 contains >6.2 million SNP-phenotype association from among 1390 GWAS studies. We re-annotated GWAS results with 16 annotation sources including some rarely compared to GWAS results (e.g. RNAediting sites, lincRNAs, PTMs).

MOTIVATION

To create a high-quality resource to facilitate further use and interpretation of human GWAS results in order to address important scientific questions.

RESULTS

GWAS have grown exponentially, with increases in sample sizes and markers tested, and continuing bias toward European ancestry samples. GRASP contains >100 000 phenotypes, roughly: eQTLs (71.5%), metabolite QTLs (21.2%), methylation QTLs (4.4%) and diseases, biomarkers and other traits (2.8%). cis-eQTLs, meQTLs, mQTLs and MHC region SNPs are highly enriched among significant results. After removing these categories, GRASP still contains a greater proportion of studies and results than comparable GWAS catalogs. Cardiovascular disease and related risk factors pre-dominate remaining GWAS results, followed by immunological, neurological and cancer traits. Significant results in GWAS display a highly gene-centric tendency. Sex chromosome X (OR = 0.18[0.16-0.20]) and Y (OR = 0.003[0.001-0.01]) genes are depleted for GWAS results. Gene length is correlated with GWAS results at nominal significance (P ≤ 0.05) levels. We show this gene-length correlation decays at increasingly more stringent P-value thresholds. Potential pleotropic genes and SNPs enriched for multi-phenotype association in GWAS are identified. However, we note possible population stratification at some of these loci. Finally, via re-annotation we identify compelling functional hypotheses at GWAS loci, in some cases unrealized in studies to date.

CONCLUSION

Pooling summary-level GWAS results and re-annotating with bioinformatics predictions and molecular features provides a good platform for new insights.

AVAILABILITY

The GRASP database is available at http://apps.nhlbi.nih.gov/grasp.

Orchestration of neurodevelopmental programs by RBFOX1: implications for autism spectrum disorder

Neurodevelopmental and neuropsychiatric disorders result from complex interactions between critical genetic factors and as-yet-unknown environmental components. To gain clinical insight, it is critical to develop a comprehensive understanding of these genetic components. RBFOX1, an RNA splicing factor, regulates expression of large genetic networks during early neuronal development, and haploinsufficiency causes severe neurodevelopmental phenotypes including autism spectrum disorder (ASD), intellectual disability, and epilepsy. Genomic testing in individuals and large patient cohorts has identified phenotypically similar cases possessing copy number variations in RBFOX1, implicating the gene as an important cause of neurodevelopmental disease. However, a significant proportion of the observed structural variation is inherited from phenotypically normal individuals, raising questions regarding overall pathogenicity of variation at the RBFOX1 locus. In this chapter, we discuss the molecular, cellular, and clinical evidence supporting the role of RBFOX1 in neurodevelopment and present a comprehensive model for the contribution of structural variation in RBFOX1 to ASD.

Analysis of RBFOX1 gene expression in lymphoblastoid cell lines of Italian discordant autism spectrum disorders sib-pairs

Several lines of evidence suggest that RBFOX1 is a key regulator of transcriptional and splicing programs in neural cells during development, and that it is expressed in a neuronal module enriched for known autism susceptibility genes. We have investigated its expression by semiquantitative RT-PCR in accessible nonbrain resources in eighteen autism spectrum disorder sib-pairs belonging to the Italian Autism Network cohort. RBFOX1 gene expression was detected in lymphoblastoid cell lines but not in lymphocytes. No significant differences between autism spectrum disorders and non-affected brothers were found. We were not able to replicate in lymphoblastoid cell lines the previously reported RBFOX1 gene downregulation in autism, even if a trend was observed. This might be due to less pronounced transcription level differences in RBFOX1 gene expression in lymphoblastoid cell lines than in brain samples.

Copy number variant study of bipolar disorder in Canadian and UK populations implicates synaptic genes

Genome-wide single nucleotide polymorphism (SNP) data from 936 bipolar disorder (BD) individuals and 940 psychiatrically healthy comparison individuals of North European descent were analyzed for copy number variation (CNV). Using multiple CNV calling algorithms, and validating using in vitro molecular analyses, we identified CNVs implicating several candidate genes that encode synaptic proteins, such as DLG1, DLG2, DPP6, NRXN1, NRXN2, NRXN3, SHANK2, and EPHA5, as well as the neuronal splicing regulator RBFOX1 (A2BP1), and neuronal cell adhesion molecule CHL1. We have also identified recurrent CNVs on 15q13.3 and 16p11.2-regions previously reported as risk loci for neuropsychiatric disorders. In addition, we performed CNV analysis of individuals from 215 BD trios and identified de novo CNVs involving the NRXN1 and DRD5 genes. Our study provides further evidence of the occasional involvement of genomic mutations in the etiology of BD, however, there is no evidence of an increased burden of CNVs in BD. Further, the identification of CNVs at multiple members of the neurexin gene family in BD individuals, supports the role of synaptic disruption in the etiology of BD.

Outfoxed by RBFOX1-a caution about ascertainment bias

We report on two patients with intragenic deletions of RBFOX1 and one patient with an intragenic duplication of RBFOX1. These patients, by report, all had autism spectrum disorder and/or developmental delay and had strong family histories of these conditions. We initially hypothesized that RBFOX1 was another susceptibility locus for autism spectrum disorder or developmental delay. However, epidemiological evidence examining large numbers of individuals did not support this hypothesis and the data presented here suggests that RBFOX1 intragenic copy number variants are not pathogenic. This contradicts previous reports that examined smaller numbers of patients and controls. © 2014 Wiley Periodicals, Inc.

Rbfox proteins regulate alternative mRNA splicing through evolutionarily conserved RNA bridges

Alternative splicing (AS) enables programmed diversity of gene expression across tissues and development. We show here that binding in distal intronic regions (>500 nucleotides (nt) from any exon) by Rbfox splicing factors important in development is extensive and is an active mode of splicing regulation. Similarly to exon-proximal sites, distal sites contain evolutionarily conserved GCATG sequences and are associated with AS activation and repression upon modulation of Rbfox abundance in human and mouse experimental systems. As a proof of principle, we validated the activity of two specific Rbfox enhancers in KIF21A and ENAH distal introns and showed that a conserved long-range RNA-RNA base-pairing interaction (an RNA bridge) is necessary for Rbfox-mediated exon inclusion in the ENAH gene. Thus we demonstrate a previously unknown RNA-mediated mechanism for AS control by distally bound RNA-binding proteins.

RBFOX1 and RBFOX3 mutations in rolandic epilepsy

Partial deletions of the gene encoding the neuronal splicing regulator RBFOX1 have been reported in a range of neurodevelopmental diseases, including idiopathic generalized epilepsy. The RBFOX1 protein and its homologues (RBFOX2 and RBFOX3) regulate alternative splicing of many neuronal transcripts involved in the homeostatic control of neuronal excitability. In this study, we explored if structural microdeletions and exonic sequence variations in RBFOX1, RBFOX2, RBFOX3 confer susceptibility to rolandic epilepsy (RE), a common idiopathic focal childhood epilepsy. By high-density SNP array screening of 289 unrelated RE patients, we identified two hemizygous deletions, a 365 kb deletion affecting two untranslated 5′-terminal exons of RBFOX1 and a 43 kb deletion spanning exon 3 of RBFOX3. Exome sequencing of 242 RE patients revealed two novel probably deleterious variants in RBFOX1, a frameshift mutation (p.A233Vfs*74) and a hexanucleotide deletion (p.A299_A300del), and a novel nonsense mutation in RBFOX3 (p.Y287*). Although the three variants were inherited from unaffected parents, they were present in all family members exhibiting the RE trait clinically or electroencephalographically with only one exception. In contrast, no deleterious mutations of RBFOX1 and RBFOX3 were found in the exomes of 6503 non-RE subjects deposited in the Exome Variant Server database. The observed RBFOX3 exon 3 deletion and nonsense mutation suggest that RBFOX3 represents a novel risk factor for RE, indicating that exon deletions and truncating mutations of RBFOX1 and RBFOX3 contribute to the genetic variance of partial and generalized idiopathic epilepsy syndromes.

An eQTL mapping approach reveals that rare variants in the SEMA5A regulatory network impact autism risk

To date, genome-wide single nucleotide polymorphism (SNP) and copy number variant (CNV) association studies of autism spectrum disorders (ASDs) have led to promising signals but not to easily interpretable or translatable results. Our own genome-wide association study (GWAS) showed significant association to an intergenic SNP near Semaphorin 5A (SEMA5A) and provided evidence for reduced expression of the same gene. In a novel GWAS follow-up approach, we map an expression regulatory pathway for a GWAS candidate gene, SEMA5A, in silico by using population expression and genotype data sets. We find that the SEMA5A regulatory network significantly overlaps rare autism-specific CNVs. The SEMA5A regulatory network includes previous autism candidate genes and regions, including MACROD2, A2BP1, MCPH1, MAST4, CDH8, CADM1, FOXP1, AUTS2, MBD5, 7q21, 20p, USH2A, KIRREL3, DBF4B and RELN, among others. Our results provide: (i) a novel data-derived network implicated in autism, (ii) evidence that the same pathway seeded by an initial SNP association shows association with rare genetic variation in ASDs, (iii) a potential mechanism of action and interpretation for the previous autism candidate genes and genetic variants that fall in this network, and (iv) a novel approach that can be applied to other candidate genes for complex genetic disorders. We take a step towards better understanding of the significance of SEMA5A pathways in autism that can guide interpretation of many other genetic results in ASDs.

Intragenic deletion of RBFOX1 associated with neurodevelopmental/neuropsychiatric disorders and possibly other clinical presentations

BACKGROUND

RBFOX1 is an important splicing factor regulating developmental and tissue-specific alternative splicing in heart, muscle, and neuronal tissues. Constitutional genetic defects in RBFOX1 are implicated in multiple medical conditions.

RESULTS

We identified 14 copy number variants (CNV) involving RBFOX1 from 2,124 consecutive pediatric patients referred for chromosomal microarray analysis (CMA), including 13 intragenic deletions and a single intragenic duplication. The clinical significances of the intragenic deletions of RBFOX1 were evaluated.

CONCLUSIONS

Our data strongly supports the associations of intragenic deletions of RBFOX1 with a diversity of neurodevelopmental and neuropsychiatric disorders, and possibly other clinical features.

Aberrant splicing in neurological diseases

Splicing of precursor messenger RNA (pre-mRNA) removes the intervening sequences (introns) and joins the expressed regions (exons) in the nucleus, before an intron-containing eukaryotic mRNA transcript can be exported and translated into proteins in the cytoplasm. While some sequences are always included or excluded (constitutive splicing), others can be selectively used (alternative splicing) in this process. Particularly by alternative splicing, up to tens of thousands of variant transcripts can be produced from a single gene, which contributes greatly to the proteomic diversity for such complex cellular functions as 'wiring' neurons in the nervous system. Disruption of this process leads to aberrant splicing, which accounts for the defects of up to 50% of mutations that cause certain human genetic diseases. In this review, we describe the different mechanisms of aberrant splicing that cause or have been associated with neurological diseases.

Rbfox1 downregulation and altered calpain 3 splicing by FRG1 in a mouse model of Facioscapulohumeral muscular dystrophy (FSHD)

Facioscapulohumeral muscular dystrophy (FSHD) is a common muscle disease whose molecular pathogenesis remains largely unknown. Over-expression of FSHD region gene 1 (FRG1) in mice, frogs, and worms perturbs muscle development and causes FSHD-like phenotypes. FRG1 has been implicated in splicing, and we asked how splicing might be involved in FSHD by conducting a genome-wide analysis in FRG1 mice. We find that splicing perturbations parallel the responses of different muscles to FRG1 over-expression and disease progression. Interestingly, binding sites for the Rbfox family of splicing factors are over-represented in a subset of FRG1-affected splicing events. Rbfox1 knockdown, over-expression, and RNA-IP confirm that these are direct Rbfox1 targets. We find that FRG1 is associated to the Rbfox1 RNA and decreases its stability. Consistent with this, Rbfox1 expression is down-regulated in mice and cells over-expressing FRG1 as well as in FSHD patients. Among the genes affected is Calpain 3, which is mutated in limb girdle muscular dystrophy, a disease phenotypically similar to FSHD. In FRG1 mice and FSHD patients, the Calpain 3 isoform lacking exon 6 (Capn3 E6-) is increased. Finally, Rbfox1 knockdown and over-expression of Capn3 E6- inhibit muscle differentiation. Collectively, our results suggest that a component of FSHD pathogenesis may arise by over-expression of FRG1, reducing Rbfox1 levels and leading to aberrant expression of an altered Calpain 3 protein through dysregulated splicing.