Biallelic missense variants in ZBTB11 can cause intellectual disability in humans.

Exploring genes and pathways underlying intellectual disability (ID) provides insight into brain development and function, clarifying the complex puzzle of how cognition develops. As part of ongoing systematic studies to identify candidate ID genes, linkage analysis and next-generation sequencing revealed Zinc Finger and BTB Domain Containing 11 (ZBTB11) as a novel candidate ID gene. ZBTB11 encodes a little-studied transcription regulator, and the two identified missense variants in this study are predicted to disrupt canonical Zn2+-binding residues of its C2H2 zinc finger domain, leading to possible altered DNA binding. Using HEK293T cells transfected with wild-type and mutant GFP-ZBTB11 constructs, we found the ZBTB11 mutants being excluded from the nucleolus, where the wild-type recombinant protein is predominantly localized. Pathway analysis applied to ChIP-seq data deposited in the ENCODE database supports the localization of ZBTB11 in nucleoli, highlighting associated pathways such as ribosomal RNA synthesis, ribosomal assembly, RNA modification and stress sensing, and provides a direct link between subcellular ZBTB11 location and its function. Furthermore, given the report of prominent brain and spinal cord degeneration in a zebrafish Zbtb11 mutant, we investigated ZBTB11-ortholog knockdown in Drosophila melanogaster brain by targeting RNAi using the UAS/Gal4 system. The observed approximate reduction to a third of the mushroom body size-possibly through neuronal reduction or degeneration-may affect neuronal circuits in the brain that are required for adaptive behavior, specifying the role of this gene in the nervous system. In conclusion, we report two ID families segregating ZBTB11 biallelic mutations disrupting Zn2+-binding motifs and provide functional evidence linking ZBTB11 dysfunction to this phenotype.

Genetics of intellectual disability in consanguineous families.

Autosomal recessive (AR) gene defects are the leading genetic cause of intellectual disability (ID) in countries with frequent parental consanguinity, which account for about 1/7th of the world population. Yet, compared to autosomal dominant de novo mutations, which are the predominant cause of ID in Western countries, the identification of AR-ID genes has lagged behind. Here, we report on whole exome and whole genome sequencing in 404 consanguineous predominantly Iranian families with two or more affected offspring. In 219 of these, we found likely causative variants, involving 77 known and 77 novel AR-ID (candidate) genes, 21 X-linked genes, as well as 9 genes previously implicated in diseases other than ID. This study, the largest of its kind published to date, illustrates that high-throughput DNA sequencing in consanguineous families is a superior strategy for elucidating the thousands of hitherto unknown gene defects underlying AR-ID, and it sheds light on their prevalence.

Just Expect It: Compound Heterozygous Variants of POMT1 in a Consanguineous Family-The Role of Next Generation Sequencing in Neuromuscular Disorders.

Muscular dystrophy-dystroglycanopathies (MDDG) are a group of genetically heterogeneous autosomal recessive disorders characterized by hypoglycosylation of α-dystroglycan. Here, we report on two female patients from a consanguineous Lebanese family that presented in early infancy with generalized muscle hypotonia and primary microcephaly. Brain magnetic resonance imaging (MRI) showed different degrees of hypoplasia of the cerebellar vermis and hypoplasia of corpus callosum. Muscle biopsy analyses revealed a muscular dystrophy with reduced expression of α-dystroglycan and merosin in immunoblot analyses. Homozygosity mapping failed to elucidate the causal mutation due to the accepted notion that, in consanguineous families, homozygote mutations cause disease. However, by applying whole exome sequencing, we identified a novel compound heterozygous POMT1 mutation that segregates with the phenotype and is in line with the clinical presentation. This underscores that a less expected compound heterozygous instead of homozygous mutation in a consanguineous marriage results in a recessive disorder and highlights the growing role of next generation sequencing in neuromuscular disorder diagnostics.

Homozygous ARHGEF2 mutation causes intellectual disability and midbrain-hindbrain malformation.

Mid-hindbrain malformations can occur during embryogenesis through a disturbance of transient and localized gene expression patterns within these distinct brain structures. Rho guanine nucleotide exchange factor (ARHGEF) family members are key for controlling the spatiotemporal activation of Rho GTPase, to modulate cytoskeleton dynamics, cell division, and cell migration. We identified, by means of whole exome sequencing, a homozygous frameshift mutation in the ARHGEF2 as a cause of intellectual disability, a midbrain-hindbrain malformation, and mild microcephaly in a consanguineous pedigree of Kurdish-Turkish descent. We show that loss of ARHGEF2 perturbs progenitor cell differentiation and that this is associated with a shift of mitotic spindle plane orientation, putatively favoring more symmetric divisions. The ARHGEF2 mutation leads to reduction in the activation of the RhoA/ROCK/MLC pathway crucial for cell migration. We demonstrate that the human brain malformation is recapitulated in Arhgef2 mutant mice and identify an aberrant migration of distinct components of the precerebellar system as a pathomechanism underlying the midbrain-hindbrain phenotype. Our results highlight the crucial function of ARHGEF2 in human brain development and identify a mutation in ARHGEF2 as novel cause of a neurodevelopmental disorder.

Effect of inbreeding on intellectual disability revisited by trio sequencing.

In outbred Western populations, most individuals with intellectual disability (ID) are sporadic cases, dominant de novo mutations (DNM) are frequent, and autosomal recessive ID (ARID) is very rare. Because of the high rate of parental consanguinity, which raises the risk for ARID and other recessive disorders, the prevalence of ID is significantly higher in near- and middle-east countries. Indeed, homozygosity mapping and sequencing in consanguineous families have already identified a plethora of ARID genes, but because of the design of these studies, DNMs could not be systematically assessed, and the proportion of cases that are potentially preventable by avoiding consanguineous marriages or through carrier testing is hitherto unknown. This prompted us to perform whole-exome sequencing in 100 sporadic ID patients from Iran and their healthy consanguineous parents. In 61 patients, we identified apparently causative changes in known ID genes. Of these, 44 were homozygous recessive and 17 dominant DNMs. Assuming that the DNM rate is stable, these results suggest that parental consanguinity raises the ID risk about 3.6-fold, and about 4.1 to 4.25-fold for children of first-cousin unions. These results do not rhyme with recent opinions that consanguinity-related health risks are generally small and have been "overstated" in the past.

Homozygous SLC6A17 mutations cause autosomal-recessive intellectual disability with progressive tremor, speech impairment, and behavioral problems.

We report on Dutch and Iranian families with affected individuals who present with moderate to severe intellectual disability and additional phenotypes including progressive tremor, speech impairment, and behavioral problems in certain individuals. A combination of exome sequencing and homozygosity mapping revealed homozygous mutations c.484G>A (p.Gly162Arg) and c.1898C>G (p.Pro633Arg) in SLC6A17. SLC6A17 is predominantly expressed in the brain, encodes a synaptic vesicular transporter of neutral amino acids and glutamate, and plays an important role in the regulation of glutamatergic synapses. Prediction programs and 3D modeling suggest that the identified mutations are deleterious to protein function. To directly test the functional consequences, we investigated the neuronal subcellular localization of overexpressed wild-type and mutant variants in mouse primary hippocampal neuronal cells. Wild-type protein was present in soma, axons, dendrites, and dendritic spines. p.Pro633Arg altered SLC6A17 was found in soma and proximal dendrites but did not reach spines. p.Gly162Arg altered SLC6A17 showed a normal subcellular distribution but was associated with an abnormal neuronal morphology mainly characterized by the loss of dendritic spines. In summary, our genetic findings implicate homozygous SLC6A17 mutations in autosomal-recessive intellectual disability, and their pathogenic role is strengthened by genetic evidence and in silico and in vitro functional analyses.

Schizophrenia risk polymorphisms in the TCF4 gene interact with smoking in the modulation of auditory sensory gating.

Several polymorphisms of the transcription factor 4 (TCF4) have been shown to increase the risk for schizophrenia, particularly TCF4 rs9960767. This polymorphism is associated with impaired sensorimotor gating measured by prepulse inhibition--an established endophenotype of schizophrenia. We therefore investigated whether TCF4 polymorphisms also affect another proposed endophenotype of schizophrenia, namely sensory gating assessed by P50 suppression of the auditory evoked potential. Although sensorimotor gating and sensory gating are not identical, recent data suggest that they share genetic fundamentals. In a multicenter study at six academic institutions throughout Germany, we applied an auditory P50 suppression paradigm to 1,821 subjects (1,023 never-smokers, 798 smokers) randomly selected from the general population. Samples were genotyped for 21 TCF4 polymorphisms. Given that smoking is highly prevalent in schizophrenia and affects sensory gating, we also assessed smoking behavior, cotinine plasma concentrations, exhaled carbon monoxide, and the Fagerström Test (FTND). P50 suppression was significantly decreased in carriers of schizophrenia risk alleles of the TCF4 polymorphisms rs9960767, rs10401120rs, rs17597926, and 17512836 (P < 0.0002-0.00005). These gene effects were modulated by smoking behavior as indicated by significant interactions of TCF4 genotype and smoking status; heavy smokers (FTND score ≥ 4) showed stronger gene effects on P50 suppression than light smokers and never-smokers. Our finding suggests that sensory gating is modulated by an interaction of TCF4 genotype with smoking, and both factors may play a role in early information processing deficits also in schizophrenia. Consequently, considering smoking behavior may facilitate the search for genetic risk factors for schizophrenia.

Integrated sequence analysis pipeline provides one-stop solution for identifying disease-causing mutations.

Next-generation sequencing has greatly accelerated the search for disease-causing defects, but even for experts the data analysis can be a major challenge. To facilitate the data processing in a clinical setting, we have developed a novel medical resequencing analysis pipeline (MERAP). MERAP assesses the quality of sequencing, and has optimized capacity for calling variants, including single-nucleotide variants, insertions and deletions, copy-number variation, and other structural variants. MERAP identifies polymorphic and known causal variants by filtering against public domain databases, and flags nonsynonymous and splice-site changes. MERAP uses a logistic model to estimate the causal likelihood of a given missense variant. MERAP considers the relevant information such as phenotype and interaction with known disease-causing genes. MERAP compares favorably with GATK, one of the widely used tools, because of its higher sensitivity for detecting indels, its easy installation, and its economical use of computational resources. Upon testing more than 1,200 individuals with mutations in known and novel disease genes, MERAP proved highly reliable, as illustrated here for five families with disease-causing variants. We believe that the clinical implementation of MERAP will expedite the diagnostic process of many disease-causing defects.

Analysis of quantitative trait loci in mice suggests a role of Enoph1 in stress reactivity.

Significant progress in elucidating the genetic etiology of anxiety and depression has been made during the last decade through a combination of human and animal studies. In this study, we aimed to discover genetic loci linked with anxiety as well as depression in order to reveal new candidate genes. Therefore, we initially tested the behavioral sensitivity of 543 F2 animals derived from an intercross of C57BL/6J and C3H/HeJ mice in paradigms for anxiety and depression. Next, all animals were genotyped with 269 microsatellite markers with a mean distance of 5.56 cM. Finally, a Quantitative Trait Loci (QTL) analysis was carried out, followed by selection of candidate genes. The QTL analysis revealed several new QTL on chromosome 5 with a common core interval of 19 Mb. We further narrowed this interval by comparative genomics to a region of 15 Mb. A database search and gene prioritization revealed Enoph1 as the most significant candidate gene on the prioritization list for anxiety and also for depression fulfilling our selection criteria. The Enoph1 gene, which is involved in polyamine biosynthesis, is differently expressed in parental strains, which have different brain spermidine levels and show distinct anxiety and depression-related phenotype. Our result suggests a significant role in polyamines in anxiety and depression-related behaviors.

NDST1 missense mutations in autosomal recessive intellectual disability.

NDST1 was recently proposed as a candidate gene for autosomal recessive intellectual disability in two families. It encodes a bifunctional GlcNAc N-deacetylase/N-sulfotransferase with important functions in heparan sulfate biosynthesis. In mice, Ndst1 is crucial for embryonic development and homozygous null mutations are perinatally lethal. We now report on two additional unrelated families with homozygous missense NDST1 mutations. All mutations described to date predict the substitution of conserved amino acids in the sulfotransferase domain, and mutation modeling predicts drastic alterations in the local protein conformation. Comparing the four families, we noticed significant overlap in the clinical features, including both demonstrated and apparent intellectual disability, muscular hypotonia, epilepsy, and postnatal growth deficiency. Furthermore, in Drosophila, knockdown of sulfateless, the NDST ortholog, impairs long-term memory, highlighting its function in cognition. Our data confirm NDST1 mutations as a cause of autosomal recessive intellectual disability with a distinctive phenotype, and support an important function of NDST1 in human development.

A newly recognized autosomal recessive syndrome affecting neurologic function and vision.

Genetic factors represent an important etiologic group in the causation of intellectual disability. We describe a Saudi Arabian family with closley related parents in which four of six children were affected by a congenital cognitive disturbance. The four individuals (aged 18, 16, 13, and 2 years when last examined) had motor and cognitive delay with seizures in early childhood, and three of the four (sparing only the youngest child) had progressive, severe cognitive decline with spasticity. Two affected children had ocular malformations, and the three older children had progressive visual loss. The youngest had normal globes with good functional vision when last examined but exhibited the oculodigital sign, which may signify a subclinical visual deficit. A potentially deleterious nucleotide change (c.1A>G; p.Met1Val) in the C12orf57 gene was homozygous in all affected individuals, heterozygous in the parents, and absent in an unaffected sibling and >350 normal individuals. This gene has no known function. This family manifests a autosomal recessive syndrome with some phenotypic variability that includes abnormal development of brain and eyes, delayed cognitive and motor milestones, seizures, and a severe cognitive and visual decline that is associated with a homozygous variant in a newly identified gene.

High levels of brain-type creatine kinase activity in human platelets and leukocytes: a genetic anomaly with autosomal dominant inheritance.

The ectopic expression in peripheral blood cells of the brain-type creatine kinase (CKB) is an autosomal dominant inherited anomaly named CKBE (MIM ID 123270). Here, we characterized the CK activity in serum, platelets (PLT) and leukocytes (WBC) of 22 probands (from 8 unrelated families) and 10 controls. CK activity was measured by standard UV-photometry. Expression of the CKB gene was analyzed by real-time PCR and Western blotting. DNA sequencing including bisulfite treatment was used for molecular analysis of the CKB gene. Serum CK levels were comparable between probands and controls. CKBE probands revealed significantly higher CK activity in PLT (3.7 ± 2.7 versus 179.2 ± 83.0 U/10(12) PLT; p<0.001) and WBC (0.4 ± 0.3 versus 2.6 ± 2.1 U/10(9) WBC; p=0.004). Inhibitory anti-CKM antibodies did not affect CK activity indicating that the CK activity is generated exclusively by the CK-BB isoenzyme. CKB mRNA and protein levels were significantly higher in PLT and WBC from probands compared to controls. Re-sequencing of the entire CKB gene and methylation analysis of a CpG island revealed no alteration in CKBE probands. The genetic basis of CKBE remains unclear, however, we propose that a de-methylated CKB gene is inherited that leads to high CKB expression levels in myeloic precursor cells in the bone marrow.

Evaluating SKI as a candidate gene for non-syndromic cleft lip with or without cleft palate.

Non-syndromic cleft lip with or without cleft palate (NSCL/P) is one of the most common of all congenital malformations and has a multifactorial etiology. Findings in mice suggest that the v-ski sarcoma viral oncogene homolog (SKI) gene is a candidate gene for orofacial clefting. In humans, a significant association between rs2843159 within SKI and NSCL/P has been reported in patients from the Philippines and South America. In the South American patients, the association was driven by the subgroup of patients with non-syndromic cleft lip only (NSCLO). Here we investigated the association with rs2843159 in a Mayan Mesoamerican population (172 NSCL/P patients and 366 controls). In addition, we analyzed the phenotypic subgroups NSCLO and non-syndromic cleft of lip and palate (NSCLP). A trend towards association between rs2843159 and NSCL/P was observed in the Mayan cohort (P = 0.097), and we found a stronger association in the NSCLP subgroup (P = 0.072) despite a limited sample size. To investigate whether other common variants within the SKI gene contribute to NSCL/P susceptibility in European and Asian populations, we also analyzed genotypic data from two recent genome-wide association studies using set-based statistical approaches. These analyses detected a trend toward association in the European population. Our data provide limited support for the hypothesis that common SKI variants are susceptibility factors for NSCL/P.

Genes that determine immunology and inflammation modify the basic defect of impaired ion conductance in cystic fibrosis epithelia.

BACKGROUND: The cystic fibrosis (CF) basic defect, caused by dysfunction of the apical chloride channel CFTR in the gastrointestinal and respiratory tract epithelia, has not been employed so far to support the role of CF modifier genes. METHODS: Patients were selected from 101 families with a total of 171 F508del-CFTR homozygous CF patients to identify CF modifying genes. A candidate gene based association study of 52 genes on 16 different chromosomes with a total of 182 genetic markers was performed. Differences in haplotype and/or diplotype distribution between case and reference CF subpopulations were analysed. RESULTS: Variants at immunologically relevant genes were associated with the manifestation of the CF basic defect (0.01

Genome-wide SNP-based linkage scan identifies a locus on 8q24 for an age-related hearing impairment trait.

Age-related hearing impairment (ARHI), or presbycusis, is a very common multifactorial disorder. Despite the knowledge that genetics play an important role in the etiology of human ARHI as revealed by heritability studies, to date, its precise genetic determinants remain elusive. Here we report the results of a cross-sectional family-based genetic study employing audiometric data. By using principal component analysis, we were able to reduce the dimensionality of this multivariate phenotype while capturing most of the variation and retaining biologically important features of the audiograms. We conducted a genome-wide association as well as a linkage scan with high-density SNP microarrays. Because of the presence of genetic population substructure, association testing was stratified after which evidence was combined by meta-analysis. No association signals reaching genome-wide significance were detected. Linkage analysis identified a linkage peak on 8q24.13-q24.22 for a trait correlated to audiogram shape. The signal reached genome-wide significance, as assessed by simulations. This finding represents the first locus for an ARHI trait.

Feasible and successful: genome-wide interaction analysis involving all 1.9 x 10(11) pair-wise interaction tests.

The Genome-Wide Association Study (GWAS) is the study design of choice for detecting common genetic risk factors for multifactorial diseases. The performance of full Genome-Wide Interaction Analyses (GWIA) has always been considered computationally challenging. Two-stage strategies to reduce the amount of numerical analysis require the detection of single marker effects or prior pathophysiological hypotheses before the analysis of interaction. This prevents the detection of pure epistatic effects. Our case-control study in idiopathic generalized epilepsy demonstrates that a full GWIA is feasible through use of data compression, specific data representation, interleaved data organization, and parallelization of the analysis on a multi-processor system. Following extensive quality control of the genotypes, our final list of top interaction hits contains only pairs of interacting SNPs with negligible marginal effects. The TOP HIT interaction was between a SNP-pair intragenic to gene DNER (chr 2) and gene CTNNA3 (chr 10). Both of these genes are functionally involved in neuronal migration, synaptogenesis, and the formation of neuronal circuits. Our results therefore indicate a possible interaction between these two genes in epileptogenesis. Results from GWAS are beginning to reveal a 'missing heritability' in complex traits and diseases. Systematic, hypothesis-free analysis of epistatic interaction (GWIA) may help to close this increasingly recognized gap in heritability.

The grainyhead like 2 gene (GRHL2), alias TFCP2L3, is associated with age-related hearing impairment.

Age-related hearing impairment (ARHI) is the most prevalent sensory impairment in the elderly. ARHI is a complex disease caused by an interaction between environmental and genetic factors. The contribution of various environmental factors has been relatively extensively studied. In contrast, investigations to identify the genetic risk factors have only recently been initiated. In this paper we describe the results of an association study performed on 2418 ARHI samples derived from nine centers from seven European countries. In 70 candidate genes, a total of 768 tag single nucleotide polymorphisms (SNPs) were selected based on HAPMAP data. These genes were chosen among the monogenic hearing loss genes identified in mice and men in addition to several strong functional candidates. After genotyping and data polishing, statistical analysis of all samples combined resulted in a P-value that survived correction for multiple testing for one SNP in the GRHL2 gene. Other SNPs in this gene were also associated, albeit to a lesser degree. Subsequently, an analysis of the most significant GRHL2 SNP was performed separately for each center. The direction of the association was identical in all nine centers. Two centers showed significant associations and a third center showed a trend towards significance. Subsequent fine mapping of this locus demonstrated that the majority of the associated SNPs reside in intron 1. We hypothesize that the causative variant may change the expression levels of a GRHL2 isoform.

Hierarchical fine mapping of the cystic fibrosis modifier locus on 19q13 identifies an association with two elements near the genes CEACAM3 and CEACAM6.

On 19q13, TGFB1 and the cystic fibrosis modifier 1 locus (CFM1) have been identified as modifiers of the course of the monogenic disease cystic fibrosis (CF). Recently, we have described a transmission disequilibrium at the microsatellite D19S197, localized between TGFB1 and CFM1. To map the corresponding molecular variants, we have selected informative SNP markers within a 600-kb area and compared two-marker-haplotype-distributions between phenotypically contrasting sib pair groups, intending to type only phylogenetically old markers by aiming for close-to-maximal polymorphism information content of the SNPs. Starting with a seed set of five SNPs that cover intermarker distances of up to 50 kb, we have iteratively added more SNPs to the map, until we could identify two genomic fragments of 3,289 and 2,052 bp for which pairs with contrasting phenotypes showed different haplotype distributions on the final 17-SNP-map (P(raw) = 0.0002, P(corr17SNPs) = 0.0106 and P(raw) = 0.0008, P(corr17SNPs) = 0.0469, respectively). Resequencing of these fragments of four unrelated individuals for each element showed that the mildly and severely affected pairs differ in seven SNPs and concordant pairs differ from discordant pairs in five SNPs. Annotation of these variants indicate that CEACAM6 and a regulatory element near the 3' end of CEACAM3 are associated with CF disease severity and intrapair discordance, respectively. While our approach was only guided by the markers' position, the involvement of genes from the CEACAM family in host defense and innate immunity designates these proteins as likely modifiers of the multi-organ disease cystic fibrosis which is known for its cytokine imbalance and pro-inflammatory phenotype.

Quantitative trait loci contributing to physiological and behavioural ethanol responses after acute and chronic treatment.

The aim of the present study was the identification of gene loci that contribute to the development and manifestation of behaviours related to acute and chronic alcohol exposure, as well as to alcohol withdrawal. For this purpose, we performed a serial behavioural phenotyping of 534 animals from the second filial (F2) generation of a C57BL/6J and C3H/HeJ mice intercross in paradigms with relevance to alcohol dependence. First, ethanol-induced hypothermia was determined in ethanol-naive animals. The mice then received an ethanol solution for several weeks as their only fluid source. Ethanol tolerance, locomotor activity and anxiety-related behaviours were evaluated. The ethanol was next withdrawn and the withdrawal severity was assessed. The ethanol-experienced animals were finally analysed in a two-bottle choice paradigm to determine ethanol preference and stress-induced changes in ethanol preference. The genotypes of these mice were subsequently assessed by microsatellite marker mapping. We genotyped 264 markers with an average marker distance of 5.56 cM, which represents a high-density whole genome coverage. Quantitative trait loci (QTL) were subsequently identified using univariate analysis performed with the R/qtl tool, which is an extensible, interactive environment for mapping QTL in experimental crosses. We found QTL that have already been published, thus validating the serial phenotyping protocol, and identified several novel loci. Our analysis demonstrates that the various responses to ethanol are regulated by independent groups of genes.

A new sex-specific genetic map of the human pseudoautosomal regions (PAR1 and PAR2).

OBJECTIVES: Accurate sex-specific genetic maps are essential for gene mapping projects when using multipoint methods, especially for the pseudoautosomal regions (PARs). This paper describes a new sex-specific genetic map for the human PARs. METHODS: We used multipoint linkage analysis, three-generation pedigrees and to date, the largest set of genetic markers to estimate recombination events in the PARs. The new map combines genotype data from 28 CEPH pedigrees and 29 genetic markers. All genetic markers have been physically mapped, thus providing the opportunity to compare genetic to physical distance through all the intervals. RESULTS: Recombination activity in PAR1 is on average 7-fold higher in male than in female meiosis, resulting in a much larger genetic map for males. This new map shows how genetically different in size is the X chromosome between males and females. CONCLUSIONS: Genetic distances vary extremely between males and females in the PARs. These differences should be considered in genetic studies using accurate sex-specific maps instead of sex-averaged maps. The new map provides a significantly different length compared to some previously published maps, since a much higher number of genetic markers is used. Therefore, we conclude that our map is the most accurate resource for obtaining genetic map information for the PARs.