Mobile interspersed repeats are major structural variants in the human genome.

Characterizing structural variants in the human genome is of great importance, but a genome wide analysis to detect interspersed repeats has not been done. Thus, the degree to which mobile DNAs contribute to genetic diversity, heritable disease, and oncogenesis remains speculative. We perform transposon insertion profiling by microarray (TIP-chip) to map human L1(Ta) retrotransposons (LINE-1 s) genome-wide. This identified numerous novel human L1(Ta) insertional polymorphisms with highly variant allelic frequencies. We also explored TIP-chip's usefulness to identify candidate alleles associated with different phenotypes in clinical cohorts. Our data suggest that the occurrence of new insertions is twice as high as previously estimated, and that these repeats are under-recognized as sources of human genomic and phenotypic diversity. We have just begun to probe the universe of human L1(Ta) polymorphisms, and as TIP-chip is applied to other insertions such as Alu SINEs, it will expand the catalog of genomic variants even further.

Identification and characterization of a missense mutation in the O-linked ß-N-acetylglucosamine (O-GlcNAc) transferase gene that segregates with X-linked intellectual disability.

O-GlcNAc is a regulatory post-translational modification of nucleocytoplasmic proteins that has been implicated in multiple biological processes, including transcription. In humans, single genes encode enzymes for its attachment (O-GlcNAc transferase (OGT)) and removal (O-GlcNAcase (OGA)). An X-chromosome exome screen identified a missense mutation, which encodes an amino acid in the tetratricopeptide repeat, in OGT (759G>T (p.L254F)) that segregates with X-linked intellectual disability (XLID) in an affected family. A decrease in steady-state OGT protein levels was observed in isolated lymphoblastoid cell lines from affected individuals, consistent with molecular modeling experiments. Recombinant expression of L254F-OGT demonstrated that the enzyme is active as both a glycosyltransferase and an HCF-1 protease. Despite the reduction in OGT levels seen in the L254F-OGT individual cells, we observed that steady-state global O-GlcNAc levels remained grossly unaltered. Surprisingly, lymphoblastoids from affected individuals displayed a marked decrease in steady-state OGA protein and mRNA levels. We observed an enrichment of the OGT-containing transcriptional repressor complex mSin3A-HDAC1 at the proximal promoter region of OGA and correspondingly decreased OGA promoter activity in affected cells. Global transcriptome analysis of L254F-OGT lymphoblastoids compared with controls revealed a small subset of genes that are differentially expressed. Thus, we have begun to unravel the molecular consequences of the 759G>T (p.L254F) mutation in OGT that uncovered a compensation mechanism, albeit imperfect, given the phenotype of affected individuals, to maintain steady-state O-GlcNAc levels. Thus, a single amino acid substitution in the regulatory domain (the tetratricopeptide repeat domain) of OGT, which catalyzes the O-GlcNAc post-translational modification of nuclear and cytosolic proteins, appears causal for XLID.

StereoGene: rapid estimation of genome-wide correlation of continuous or interval feature data.

MOTIVATION: Genomics features with similar genome-wide distributions are generally hypothesized to be functionally related, for example, colocalization of histones and transcription start sites indicate chromatin regulation of transcription factor activity. Therefore, statistical algorithms to perform spatial, genome-wide correlation among genomic features are required. RESULTS: Here, we propose a method, StereoGene, that rapidly estimates genome-wide correlation among pairs of genomic features. These features may represent high-throughput data mapped to reference genome or sets of genomic annotations in that reference genome. StereoGene enables correlation of continuous data directly, avoiding the data binarization and subsequent data loss. Correlations are computed among neighboring genomic positions using kernel correlation. Representing the correlation as a function of the genome position, StereoGene outputs the local correlation track as part of the analysis. StereoGene also accounts for confounders such as input DNA by partial correlation. We apply our method to numerous comparisons of ChIP-Seq datasets from the Human Epigenome Atlas and FANTOM CAGE to demonstrate its wide applicability. We observe the changes in the correlation between epigenomic features across developmental trajectories of several tissue types consistent with known biology and find a novel spatial correlation of CAGE clusters with donor splice sites and with poly(A) sites. These analyses provide examples for the broad applicability of StereoGene for regulatory genomics. AVAILABILITY AND IMPLEMENTATION: The StereoGene C ++ source code, program documentation, Galaxy integration scripts and examples are available from the project homepage http://stereogene.bioinf.fbb.msu.ru/. CONTACT: favorov@sensi.org. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

ZC4H2, an XLID gene, is required for the generation of a specific subset of CNS interneurons.

Miles-Carpenter syndrome (MCS) was described in 1991 as an XLID syndrome with fingertip arches and contractures and mapped to proximal Xq. Patients had microcephaly, short stature, mild spasticity, thoracic scoliosis, hyperextendable MCP joints, rocker-bottom feet, hyperextended elbows and knees. A mutation, p.L66H, in ZC4H2, was identified in a XLID re-sequencing project. Additional screening of linked families and next generation sequencing of XLID families identified three ZC4H2 mutations: p.R18K, p.R213W and p.V75in15aa. The families shared some relevant clinical features. In silico modeling of the mutant proteins indicated all alterations would destabilize the protein. Knockout mutations in zc4h2 were created in zebrafish and homozygous mutant larvae exhibited abnormal swimming, increased twitching, defective eye movement and pectoral fin contractures. Because several of the behavioral defects were consistent with hyperactivity, we examined the underlying neuronal defects and found that sensory neurons and motoneurons appeared normal. However, we observed a striking reduction in GABAergic interneurons. Analysis of cell-type-specific markers showed a specific loss of V2 interneurons in the brain and spinal cord, likely arising from mis-specification of neural progenitors. Injected human wt ZC4H2 rescued the mutant phenotype. Mutant zebrafish injected with human p.L66H or p.R213W mRNA failed to be rescued, while the p.R18K mRNA was able to rescue the interneuron defect. Our findings clearly support ZC4H2 as a novel XLID gene with a required function in interneuron development. Loss of function of ZC4H2 thus likely results in altered connectivity of many brain and spinal circuits.

A distinct X-linked syndrome involving joint contractures, keloids, large optic cup-to-disc ratio, and renal stones results from a filamin A (FLNA) mutation.

We further evaluated a previously reported family with an apparently undescribed X-linked syndrome involving joint contractures, keloids, an increased optic cup-to-disc ratio, and renal stones to elucidate the genetic cause. To do this, we obtained medical histories and performed physical examination on 14 individuals in the family, five of whom are affected males and three are obligate carrier females. Linkage analysis was performed on all but one individual and chromosome X-exome sequencing was done on two affected males. The analysis localized the putative gene to Xq27-qter and chromosome X-exome sequencing revealed a mutation in exon 28 (c.4726G>A) of the filamin A (FLNA) gene, predicting that a conserved glycine had been replaced by arginine at amino acid 1576 (p.G1576R). Segregation analysis demonstrated that all known carrier females tested were heterozygous (G/A), all affected males were hemizygous for the mutation (A allele) and all normal males were hemizygous for the normal G allele. The data and the bioinformatic analysis indicate that the G1576R mutation in the FLNA gene is very likely pathogenic in this family. The syndrome affecting the family shares phenotypic overlap with other syndromes caused by FLNA mutations, but appears to be a distinct phenotype, likely representing a unique genetic syndrome.

Gain-of-function glutamate receptor interacting protein 1 variants alter GluA2 recycling and surface distribution in patients with autism.

Glutamate receptor interacting protein 1 (GRIP1) is a neuronal scaffolding protein that interacts directly with the C termini of glutamate receptors 2/3 (GluA2/3) via its PDZ domains 4 to 6 (PDZ4-6). We found an association (P<0.05) of a SNP within the PDZ4-6 genomic region with autism by genotyping autistic patients (n=480) and matched controls (n=480). Parallel sequencing identified five rare missense variants within or near PDZ4-6 only in the autism cohort, resulting in a higher cumulative mutation load (P=0.032). Two variants correlated with a more severe deficit in reciprocal social interaction in affected sibling pairs from proband families. These variants were associated with altered interactions with GluA2/3 and faster recycling and increased surface distribution of GluA2 in neurons, suggesting gain-of-function because GRIP1/2 deficiency showed opposite phenotypes. Grip1/2 knockout mice exhibited increased sociability and impaired prepulse inhibition. These results support a role for GRIP in social behavior and implicate GRIP1 variants in modulating autistic phenotype.

Affected Sib-Pair Analyses Identify Signaling Networks Associated With Social Behavioral Deficits in Autism.

Autism spectrum disorders (ASDs) are characterized by deficits in three core behavioral domains: reciprocal social interactions, communication, and restricted interests and/or repetitive behaviors. Several hundreds of risk genes for autism have been identified, however, it remains a challenge to associate these genes with specific core behavioral deficits. In multiplex autism families, affected sibs often show significant differences in severity of individual core phenotypes. We hypothesize that a higher mutation burden contributes to a larger difference in the severity of specific core phenotypes between affected sibs. We tested this hypothesis on social behavioral deficits in autism. We sequenced synaptome genes (n = 1,886) in affected male sib-pairs (n = 274) in families from the Autism Genetics Research Exchange (AGRE) and identified rare (MAF ≤ 1%) and predicted functional variants. We selected affected sib-pairs with a large (≥10; n = 92 pairs) or a small (≤4; n = 108 pairs) difference in total cumulative Autism Diagnostic Interview-Revised (ADI-R) social scores (SOCT_CS). We compared burdens of unshared variants present only in sibs with severe social deficits and found a higher burden in SOCT_CS≥10 compared to SOCT_CS ≤ 4 (SOCT_CS≥10: 705.1 ± 16.2; SOCT_CS ≤ 4, 668.3 ± 9.0; p = 0.025). Unshared SOCT_CS≥10 genes only in sibs with severe social deficits are significantly enriched in the SFARI gene set. Network analyses of these genes using InWeb_IM, molecular signatures database (MSigDB), and GeNetMeta identified enrichment for phosphoinositide 3-kinase (PI3K)-AKT-mammalian target of rapamycin (mTOR) (Enrichment Score [eScore] p value = 3.36E-07; n = 8 genes) and Nerve growth factor (NGF) (eScore p value = 8.94E-07; n = 9 genes) networks. These studies support a key role for these signaling networks in social behavioral deficits and present a novel approach to associate risk genes and signaling networks with core behavioral domains in autism.

Dysregulations of sonic hedgehog signaling in MED12-related X-linked intellectual disability disorders.

BACKGROUND: Mutations in mediator of RNA polymerase II transcription subunit 12 homolog (MED12, OMIM 300188) cause X-linked intellectual disability (XLID) disorders including FG, Lujan, and Ohdo syndromes. The Gli3-dependent Sonic Hedgehog (SHH) signaling pathway has been implicated in the original FG syndrome and Lujan syndrome. How are SHH-signaling defects related to the complex clinical phenotype of MED12-associated XLID syndromes are not fully understood. METHODS: Quantitative RT-PCR was used to study expression levels of three SHH-signaling genes in lymophoblast cell lines carrying four MED12 mutations from four unrelated XLID families. Genotype and phenotype correlation studies were performed on these mutations. RESULTS: Three newly identified and one novel MED12 mutations in six affected males from four unrelated XLID families were studied. Three mutations (c.2692A>G; p.N898D, c.3640C>T; p.R1214C, and c.3884G>A; p.R1295H) are located in the LS domain and one (c.617G>A; p.R206Q) is in the L domain of MED12. These mutations involve highly conserved amino acid residues and segregate with ID and related congenital malformations in respective probands families. Patients with the LS-domain mutations share many features of FG syndrome and some features of Lujan syndrome. The patient with the L-domain mutation presented with ID and predominant neuropsychiatric features but little dysmorphic features of either FG or Lujan syndrome. Transcript levels of three Gli3-dependent SHH-signaling genes, CREB5, BMP4, and NEUROG2, were determined by quantitative RT-PCR and found to be significantly elevated in lymphoblasts from patients with three mutations in the MED12-LS domain. CONCLUSIONS: These results support a critical role of MED12 in regulating Gli3-dependent SHH signaling and in developing ID and related congenital malformations in XLID syndromes. Differences in the expression profile of SHH-signaling genes potentially contribute to variability in clinical phenotypes in patients with MED12-related XLID disorders.

GRASP1 Regulates Synaptic Plasticity and Learning through Endosomal Recycling of AMPA Receptors.

Learning depends on experience-dependent modification of synaptic efficacy and neuronal connectivity in the brain. We provide direct evidence for physiological roles of the recycling endosome protein GRASP1 in glutamatergic synapse function and animal behavior. Mice lacking GRASP1 showed abnormal excitatory synapse number, synaptic plasticity, and hippocampal-dependent learning and memory due to a failure in learning-induced synaptic AMPAR incorporation. We identified two GRASP1 point mutations from intellectual disability (ID) patients that showed convergent disruptive effects on AMPAR recycling and glutamate uncaging-induced structural and functional plasticity. Wild-type GRASP1, but not ID mutants, rescued spine loss in hippocampal CA1 neurons in Grasp1 knockout mice. Together, these results demonstrate a requirement for normal recycling endosome function in AMPAR-dependent synaptic function and neuronal connectivity in vivo, and suggest a potential role for GRASP1 in the pathophysiology of human cognitive disorders.

Affected kindred analysis of human X chromosome exomes to identify novel X-linked intellectual disability genes.

X-linked Intellectual Disability (XLID) is a group of genetically heterogeneous disorders caused by mutations in genes on the X chromosome. Deleterious mutations in ~10% of X chromosome genes are implicated in causing XLID disorders in ~50% of known and suspected XLID families. The remaining XLID genes are expected to be rare and even private to individual families. To systematically identify these XLID genes, we sequenced the X chromosome exome (X-exome) in 56 well-established XLID families (a single affected male from 30 families and two affected males from 26 families) using an Agilent SureSelect X-exome kit and the Illumina HiSeq 2000 platform. To enrich for disease-causing mutations, we first utilized variant filters based on dbSNP, the male-restricted portions of the 1000 Genomes Project, or the Exome Variant Server datasets. However, these databases present limitations as automatic filters for enrichment of XLID genes. We therefore developed and optimized a strategy that uses a cohort of affected male kindred pairs and an additional small cohort of affected unrelated males to enrich for potentially pathological variants and to remove neutral variants. This strategy, which we refer to as Affected Kindred/Cross-Cohort Analysis, achieves a substantial enrichment for potentially pathological variants in known XLID genes compared to variant filters from public reference databases, and it has identified novel XLID candidate genes. We conclude that Affected Kindred/Cross-Cohort Analysis can effectively enrich for disease-causing genes in rare, Mendelian disorders, and that public reference databases can be used effectively, but cautiously, as automatic filters for X-linked disorders.

Effective detection of rare variants in pooled DNA samples using Cross-pool tailcurve analysis.

Sequencing targeted DNA regions in large samples is necessary to discover the full spectrum of rare variants. We report an effective Illumina sequencing strategy utilizing pooled samples with novel quality (Srfim) and filtering (SERVIC4E) algorithms. We sequenced 24 exons in two cohorts of 480 samples each, identifying 47 coding variants, including 30 present once per cohort. Validation by Sanger sequencing revealed an excellent combination of sensitivity and specificity for variant detection in pooled samples of both cohorts as compared to publicly available algorithms.