Recurrent de novo missense variants across multiple histone H4 genes underlie a neurodevelopmental syndrome.

Chromatin is essentially an array of nucleosomes, each of which consists of the DNA double-stranded fiber wrapped around a histone octamer. This organization supports cellular processes such as DNA replication, DNA transcription, and DNA repair in all eukaryotes. Human histone H4 is encoded by fourteen canonical histone H4 genes, all differing at the nucleotide level but encoding an invariant protein. Here, we present a cohort of 29 subjects with de novo missense variants in six H4 genes (H4C3, H4C4, H4C5, H4C6, H4C9, and H4C11) identified by whole-exome sequencing and matchmaking. All individuals present with neurodevelopmental features of intellectual disability and motor and/or gross developmental delay, while non-neurological features are more variable. Ten amino acids are affected, six recurrently, and are all located within the H4 core or C-terminal tail. These variants cluster to specific regions of the core H4 globular domain, where protein-protein interactions occur with either other histone subunits or histone chaperones. Functional consequences of the identified variants were evaluated in zebrafish embryos, which displayed abnormal general development, defective head organs, and reduced body axis length, providing compelling evidence for the causality of the reported disorder(s). While multiple developmental syndromes have been linked to chromatin-associated factors, missense-bearing histone variants (e.g., H3 oncohistones) are only recently emerging as a major cause of pathogenicity. Our findings establish a broader involvement of H4 variants in developmental syndromes.

A homozygous duplication of the FGG exon 8-intron 8 junction causes congenital afibrinogenemia. Lessons learned from the study of a large consanguineous Turkish family.

Congenital afibrinogenemia is the most severe congenital fibrinogen disorder, characterized by undetectable fibrinogen in circulation. Causative mutations can be divided into two main classes: null mutations with no protein production at all and missense mutations producing abnormal protein chains that are retained inside the cell. The vast majority of cases are due to single base pair mutations or small insertions or deletions in the coding regions or intron-exon junctions of FGB, FGA and FGG. Only a few large rearrangements have been described, all deletions involving FGA. Here we report the characterization of a 403 bp duplication of the FGG exon 8-intron 8 junction accounting for congenital afibrinogenemia in a large consanguineous family from Turkey. This mutation, which had escaped detection by Sanger sequencing of short polymerase chain reaction (PCR) amplicons of coding sequences and splice sites, was identified by studying multiple alignments of reads obtained from whole exome sequencing of a heterozygous individual followed by PCR amplification and sequencing of a larger portion of FGG. Because the mutation duplicates the donor splice site of intron 8, we predicted that the impact of the mutation would be on FGG transcript splicing. Analysis of mRNA produced by cells transiently transfected with normal or mutant minigene constructs showed that the duplication causes production of several aberrant FGG transcripts generating premature truncating codons.

Molecular characterization of pathogenic OTOA gene conversions in hearing loss patients.

Bi-allelic loss-of-function variants of OTOA are a well-known cause of moderate-to-severe hearing loss. Whereas non-allelic homologous recombination-mediated deletions of the gene are well known, gene conversions to pseudogene OTOAP1 have been reported in the literature but never fully described nor their pathogenicity assessed. Here, we report two unrelated patients with moderate hearing-loss, who were compound heterozygotes for a converted allele and a deletion of OTOA. The conversions were initially detected through sequencing depths anomalies at the OTOA locus after exome sequencing, then confirmed with long range polymerase chain reactions. Both conversions lead to loss-of-function by introducing a premature stop codon in exon 22 (p.Glu787*). Using genomic alignments and long read nanopore sequencing, we found that the two probands carry stretches of converted DNA of widely different lengths (at least 9 kbp and around 900 bp, respectively).

Bi-allelic Loss-of-Function Variants in DNMBP Cause Infantile Cataracts.

Infantile and childhood-onset cataracts form a heterogeneous group of disorders; among the many genetic causes, numerous pathogenic variants in additional genes associated with autosomal-recessive infantile cataracts remain to be discovered. We identified three consanguineous families affected by bilateral infantile cataracts. Using exome sequencing, we found homozygous loss-of-function variants in DNMBP: nonsense variant c.811C>T (p.Arg271∗) in large family F385 (nine affected individuals; LOD score = 5.18 at θ = 0), frameshift deletion c.2947_2948del (p.Asp983∗) in family F372 (two affected individuals), and frameshift variant c.2852_2855del (p.Thr951Metfs∗41) in family F3 (one affected individual). The phenotypes of all affected individuals include infantile-onset cataracts. RNAi-mediated knockdown of the Drosophila ortholog still life (sif), enriched in lens-secreting cells, affects the development of these cells as well as the localization of E-cadherin, alters the distribution of septate junctions in adjacent cone cells, and leads to a ∼50% reduction in electroretinography amplitudes in young flies. DNMBP regulates the shape of tight junctions, which correspond to the septate junctions in invertebrates, as well as the assembly pattern of E-cadherin in human epithelial cells. E-cadherin has an important role in lens vesicle separation and lens epithelial cell survival in humans. We therefore conclude that DNMBP loss-of-function variants cause infantile-onset cataracts in humans.

Truncating variants in the SHANK1 gene are associated with a spectrum of neurodevelopmental disorders.

PURPOSE: In this study, we aimed to characterize the clinical phenotype of a SHANK1-related disorder and define the functional consequences of SHANK1 truncating variants. METHODS: Exome sequencing (ES) was performed for six individuals who presented with neurodevelopmental disorders. Individuals were ascertained with the use of GeneMatcher and Database of Chromosomal Imbalance and Phenotype in Humans Using Ensembl Resources (DECIPHER). We evaluated potential nonsense-mediated decay (NMD) of two variants by making knock-in cell lines of endogenous truncated SHANK1, and expressed the truncated SHANK1 complementary DNA (cDNA) in HEK293 cells and cultured hippocampal neurons to examine the proteins. RESULTS: ES detected de novo truncating variants in SHANK1 in six individuals. Evaluation of NMD resulted in stable transcripts, and the truncated SHANK1 completely lost binding with Homer1, a linker protein that binds to the C-terminus of SHANK1. These variants may disrupt protein-protein networks in dendritic spines. Dispersed localization of the truncated SHANK1 variants within the spine and dendritic shaft was also observed when expressed in neurons, indicating impaired synaptic localization of truncated SHANK1. CONCLUSION: This report expands the clinical spectrum of individuals with truncating SHANK1 variants and describes the impact these variants may have on the pathophysiology of neurodevelopmental disorders.

Bi-allelic loss of ERGIC1 causes relatively mild arthrogryposis.

Arthrogryposis describes the presence of multiple joint-contractures. Clinical severity of this phenotype is variable, and more than 400 causative genes have been proposed. Among these, ERGIC1 is a recently reported candidate encoding a putative transmembrane protein of the ER-Golgi interface. Two homozygous missense variants have been reported in patients with relatively mild non-syndromic arthrogryposis. In a consanguineous family with two affected siblings presenting congenital arthrogryposis and some facial dysmorphism we performed prenatal array-CGH, postnatal targeted exome and genome sequencing. Genome sequencing identified a homozygous 22.6 Kb deletion encompassing the promoter and first exon of ERGIC1. mRNA quantification showed the complete absence of ERGIC1 expression in the two affected siblings and a decrease in heterozygous parents. Our observations validate the pathogenic role of ERGIC1 in congenital arthrogryposis and demonstrate that complete loss of function causes a relatively mild phenotype. These findings will contribute to improve genetic counseling of ERGIC1 mutations.

Visual impairment and progressive phthisis bulbi caused by recessive pathogenic variant in MARK3.

Developmental eye defects often severely reduce vision. Despite extensive efforts, for a substantial fraction of these cases the molecular causes are unknown. Recessive eye disorders are frequent in consanguineous populations and such large families with multiple affected individuals provide an opportunity to identify recessive causative genes. We studied a Pakistani consanguineous family with three affected individuals with congenital vision loss and progressive eye degeneration. The family was analyzed by exome sequencing of one affected individual and genotyping of all family members. We have identified a non-synonymous homozygous variant (NM_001128918.2: c.1708C > G: p.Arg570Gly) in the MARK3 gene as the likely cause of the phenotype. Given that MARK3 is highly conserved in flies (I: 55%; S: 67%) we knocked down the MARK3 homologue, par-1, in the eye during development. This leads to a significant reduction in eye size, a severe loss of photoreceptors and loss of vision based on electroretinogram (ERG) recordings. Expression of the par-1 p.Arg792Gly mutation (equivalent to the MARK3 variant found in patients) in developing fly eyes also induces loss of eye tissue and reduces the ERG signals. The data in flies and human indicate that the MARK3 variant corresponds to a loss of function. We conclude that the identified mutation in MARK3 establishes a new gene-disease link, since it likely causes structural abnormalities during eye development and visual impairment in humans, and that the function of MARK3/par-1 is evolutionarily conserved in eye development.

LARS2-Perrault syndrome: a new case report and literature review.

BACKGROUND: Perrault syndrome is a rare recessive and genetically heterogeneous disorder characterized by sensorineural hearing loss in males and females and gonadal dysgenesis in females. Mutations in seven different genes have been identified: HARS2, HSD17B4, CLLP, C10orf, ERAL1, TWNK and LARS2. To date, 19 variants have been reported in 18 individuals with LARS2-Perrault syndrome. CASE PRESENTATION: Here we describe the case of an 8-year-old girl with compound heterozygous missense mutations in the LARS2 gene. We identified two missense mutations [c.457A > C, p.(Asn153His) and c.1565C > A, p.(Thr522Asn)] and subsequent familial segregation showed that each parent had transmitted a mutation. CONCLUSIONS: These results have implications for genetic counseling and provide insight into the functional role of LARS2. This case highlights the importance of an early diagnosis. Systematic genetic screening of children with hearing loss allows the early identification of a Perrault syndrome in order to ensure specific endocrinological surveillance and management to prevent secondary complications. Clinical data are compared with the other cases reported in the literature.

SCN8A heterozygous variants are associated with anoxic-epileptic seizures.

Anoxic-epileptic seizures (AES) are rare outcomes of common childhood reflex anoxic syncope that trigger a true epileptic seizure. The term AES was coined by Stephenson in 1983, to differentiate these events from convulsive syncopes and the more common reflex anoxic syncopes. A genetic susceptibility for AES has been postulated; but, its molecular basis has up to now been elusive. We report here two illustrative cases and show the association of de novo SCN8A variants and AES. One of them had focal or generalized seizures and autonomic symptoms triggered by orthostatism; the second had breath-holding spells triggered by pain or exercise leading to tonic-clonic seizures; both had repeatedly normal EEGs and a family history of reflex syncope. The data of three additional AES patients further suggest, for the first time, a link between SCN8A pathogenic variants and AES. The neurodevelopment of four patients was abnormal. Four of the five SCN8A mutations observed here were previously described in patients with seizure disorders. Seizures responded particularly well to sodium channel blockers. Our observation enriches the spectrum of seizures linked with SCN8A pathogenic variants.

Complex movement disorder in a patient with heterozygous YY1 mutation (Gabriele-de Vries syndrome).

YY1 mutations cause Gabriele-de Vries syndrome, a recently described condition involving cognitive impairment, facial dysmorphism and intrauterine growth restriction. Movement disorders were reported in 5/10 cases of the original series, but no detailed description was provided. Here we present a 21-year-old woman with a mild intellectual deficit, facial dysmorphism and a complex movement disorder including an action tremor, cerebellar ataxia, dystonia, and partial ocular apraxia as the presenting and most striking feature. Whole-exome sequencing revealed a novel heterozygous de novo mutation in YY1 [NM: 003403.4 (YY1): c.907 T > C; p.(Cys303Arg)], classified as pathogenic according to the ACMG guidelines.

Domains of genome-wide gene expression dysregulation in Down's syndrome.

Trisomy 21 is the most frequent genetic cause of cognitive impairment. To assess the perturbations of gene expression in trisomy 21, and to eliminate the noise of genomic variability, we studied the transcriptome of fetal fibroblasts from a pair of monozygotic twins discordant for trisomy 21. Here we show that the differential expression between the twins is organized in domains along all chromosomes that are either upregulated or downregulated. These gene expression dysregulation domains (GEDDs) can be defined by the expression level of their gene content, and are well conserved in induced pluripotent stem cells derived from the twins' fibroblasts. Comparison of the transcriptome of the Ts65Dn mouse model of Down's syndrome and normal littermate mouse fibroblasts also showed GEDDs along the mouse chromosomes that were syntenic in human. The GEDDs correlate with the lamina-associated (LADs) and replication domains of mammalian cells. The overall position of LADs was not altered in trisomic cells; however, the H3K4me3 profile of the trisomic fibroblasts was modified and accurately followed the GEDD pattern. These results indicate that the nuclear compartments of trisomic cells undergo modifications of the chromatin environment influencing the overall transcriptome, and that GEDDs may therefore contribute to some trisomy 21 phenotypes.

Pathogenic Variants in PIGG Cause Intellectual Disability with Seizures and Hypotonia.

Glycosylphosphatidylinositol (GPI) is a glycolipid that anchors >150 various proteins to the cell surface. At least 27 genes are involved in biosynthesis and transport of GPI-anchored proteins (GPI-APs). To date, mutations in 13 of these genes are known to cause inherited GPI deficiencies (IGDs), and all are inherited as recessive traits. IGDs mainly manifest as intellectual disability, epilepsy, coarse facial features, and multiple organ anomalies. These symptoms are caused by the decreased surface expression of GPI-APs or by structural abnormalities of GPI. Here, we present five affected individuals (from two consanguineous families from Egypt and Pakistan and one non-consanguineous family from Japan) who show intellectual disability, hypotonia, and early-onset seizures. We identified pathogenic variants in PIGG, a gene in the GPI pathway. In the consanguineous families, homozygous variants c.928C>T (p.Gln310(∗)) and c.2261+1G>C were found, whereas the Japanese individual was compound heterozygous for c.2005C>T (p.Arg669Cys) and a 2.4 Mb deletion involving PIGG. PIGG is the enzyme that modifies the second mannose with ethanolamine phosphate, which is removed soon after GPI is attached to the protein. Physiological significance of this transient modification has been unclear. Using B lymphoblasts from affected individuals of the Egyptian and Japanese families, we revealed that PIGG activity was almost completely abolished; however, the GPI-APs had normal surface levels and normal structure, indicating that the pathogenesis of PIGG deficiency is not yet fully understood. The discovery of pathogenic variants in PIGG expands the spectrum of IGDs and further enhances our understanding of this etiopathogenic class of intellectual disability.

Exome sequencing identifies a de novo FOXA2 variant in a patient with syndromic diabetes.

OBJECTIVE: When diabetes is associated with congenital malformations, without autoimmune antibodies, a genetic cause is suspected. Here, we aimed to identify a defective gene that led to diabetes. RESEARCH DESIGN AND METHODS: We performed an exome analysis of an index case and his healthy parents. RESULTS: The child presented with childhood-onset diabetes, congenital hypopituitarism, cardiac malformation, and anal atresia. A DNA analysis revealed a heterozygous de novo pathogenic variant in the developmental transcription factor, forkhead box A2 (FOXA2). The mutation resided in the DNA-binding domain, which is highly conserved among species. Tridimensional molecular dynamics simulation modeling predicted an altered interaction between the mutated protein and DNA. CONCLUSIONS: A defect in the FOXA2 DNA-binding domain was associated with childhood-onset diabetes and multiple congenital anomalies, which reflected the pleiotropic nature of the gene. This report extends the recently described phenotype of neonatal hypoglycemia to later-onset diabetes. We suggest to include FOXA2 analysis for neonatal hypoglycemia and to implement a long-term follow-up, particularly for the risk of diabetes.

Biased allelic expression in human primary fibroblast single cells.

The study of gene expression in mammalian single cells via genomic technologies now provides the possibility to investigate the patterns of allelic gene expression. We used single-cell RNA sequencing to detect the allele-specific mRNA level in 203 single human primary fibroblasts over 133,633 unique heterozygous single-nucleotide variants (hetSNVs). We observed that at the snapshot of analyses, each cell contained mostly transcripts from one allele from the majority of genes; indeed, 76.4% of the hetSNVs displayed stochastic monoallelic expression in single cells. Remarkably, adjacent hetSNVs exhibited a haplotype-consistent allelic ratio; in contrast, distant sites located in two different genes were independent of the haplotype structure. Moreover, the allele-specific expression in single cells correlated with the abundance of the cellular transcript. We observed that genes expressing both alleles in the majority of the single cells at a given time point were rare and enriched with highly expressed genes. The relative abundance of each allele in a cell was controlled by some regulatory mechanisms given that we observed related single-cell allelic profiles according to genes. Overall, these results have direct implications in cellular phenotypic variability.

Biallelic variants in LINGO1 are associated with autosomal recessive intellectual disability, microcephaly, speech and motor delay.

PURPOSE: To elucidate the novel molecular cause in two unrelated consanguineous families with autosomal recessive intellectual disability. METHODS: A combination of homozygosity mapping and exome sequencing was used to locate the plausible genetic defect in family F162, while only exome sequencing was followed in the family PKMR65. The protein 3D structure was visualized with the University of California-San Francisco Chimera software. RESULTS: All five patients from both families presented with severe intellectual disability, aggressive behavior, and speech and motor delay. Four of the five patients had microcephaly. We identified homozygous missense variants in LINGO1, p.(Arg290His) in family F162 and p.(Tyr288Cys) in family PKMR65. Both variants were predicted to be pathogenic, and segregated with the phenotype in the respective families. Molecular modeling of LINGO1 suggests that both variants interfere with the glycosylation of the protein. CONCLUSION: LINGO1 is a transmembrane receptor, predominantly found in the central nervous system. Published loss-of-function studies in mouse and zebrafish have established a crucial role of LINGO1 in normal neuronal development and central nervous system myelination by negatively regulating oligodendrocyte differentiation and neuronal survival. Taken together, our results indicate that biallelic LINGO1 missense variants cause autosomal recessive intellectual disability in humans.

Novel NEXMIF pathogenic variant in a boy with severe autistic features, intellectual disability, and epilepsy, and his mildly affected mother.

Intellectual disability (ID) and autism spectrum disorders are complex neurodevelopmental disorders occurring among all ethnic and socioeconomic groups. Pathogenic variants in the neurite extension and migration factor (NEXMIF) gene (formerly named KIAA2022) on the X chromosome are responsible for ID, autistic behavior, epilepsy, or dysmorphic features in males. Most affected females described had a milder phenotype or were asymptomatic obligate carriers. We report here for the first time mother-to-son transmission of a novel NEXMIF truncating variant without X-inactivation skewing in the blood. Truncating gene variant leads to symptomatic mother to severely affected son transmission. Our findings emphasize that NEXMIF sequencing should be strongly considered in patients with unexplained autism spectrum disorder, ID, and epilepsy, irrespective of gender. Such testing could increase our knowledge of the pathogenicity of NEXMIF variants and improve genetic counseling.

Tissue-specific effects of genetic and epigenetic variation on gene regulation and splicing.

Understanding how genetic variation affects distinct cellular phenotypes, such as gene expression levels, alternative splicing and DNA methylation levels, is essential for better understanding of complex diseases and traits. Furthermore, how inter-individual variation of DNA methylation is associated to gene expression is just starting to be studied. In this study, we use the GenCord cohort of 204 newborn Europeans' lymphoblastoid cell lines, T-cells and fibroblasts derived from umbilical cords. The samples were previously genotyped for 2.5 million SNPs, mRNA-sequenced, and assayed for methylation levels in 482,421 CpG sites. We observe that methylation sites associated to expression levels are enriched in enhancers, gene bodies and CpG island shores. We show that while the correlation between DNA methylation and gene expression can be positive or negative, it is very consistent across cell-types. However, this epigenetic association to gene expression appears more tissue-specific than the genetic effects on gene expression or DNA methylation (observed in both sharing estimations based on P-values and effect size correlations between cell-types). This predominance of genetic effects can also be reflected by the observation that allele specific expression differences between individuals dominate over tissue-specific effects. Additionally, we discover genetic effects on alternative splicing and interestingly, a large amount of DNA methylation correlating to alternative splicing, both in a tissue-specific manner. The locations of the SNPs and methylation sites involved in these associations highlight the participation of promoter proximal and distant regulatory regions on alternative splicing. Overall, our results provide high-resolution analyses showing how genome sequence variation has a broad effect on cellular phenotypes across cell-types, whereas epigenetic factors provide a secondary layer of variation that is more tissue-specific. Furthermore, the details of how this tissue-specificity may vary across inter-relations of molecular traits, and where these are occurring, can yield further insights into gene regulation and cellular biology as a whole.

A large genomic deletion leads to enhancer adoption by the lamin B1 gene: a second path to autosomal dominant adult-onset demyelinating leukodystrophy (ADLD).

Chromosomal rearrangements with duplication of the lamin B1 (LMNB1) gene underlie autosomal dominant adult-onset demyelinating leukodystrophy (ADLD), a rare neurological disorder in which overexpression of LMNB1 causes progressive central nervous system demyelination. However, we previously reported an ADLD family (ADLD-1-TO) without evidence of duplication or other mutation in LMNB1 despite linkage to the LMNB1 locus and lamin B1 overexpression. By custom array-CGH, we further investigated this family and report here that patients carry a large (∼660 kb) heterozygous deletion that begins 66 kb upstream of the LMNB1 promoter. Lamin B1 overexpression was confirmed in further ADLD-1-TO tissues and in a postmortem brain sample, where lamin B1 was increased in the frontal lobe. Through parallel studies, we investigated both loss of genetic material and chromosomal rearrangement as possible causes of LMNB1 overexpression, and found that ADLD-1-TO plausibly results from an enhancer adoption mechanism. The deletion eliminates a genome topological domain boundary, allowing normally forbidden interactions between at least three forebrain-directed enhancers and the LMNB1 promoter, in line with the observed mainly cerebral localization of lamin B1 overexpression and myelin degeneration. This second route to LMNB1 overexpression and ADLD is a new example of the relevance of regulatory landscape modifications in determining Mendelian phenotypes.

Galanin pathogenic mutations in temporal lobe epilepsy.

Temporal lobe epilepsy (TLE) is a common epilepsy syndrome with a complex etiology. Despite evidence for the participation of genetic factors, the genetic basis of TLE remains largely unknown. A role for the galanin neuropeptide in the regulation of epileptic seizures has been established in animal models more than two decades ago. However, until now there was no report of pathogenic mutations in GAL, the galanin-encoding gene, and therefore its role in human epilepsy was not established. Here, we studied a family with a pair of monozygotic twins affected by TLE and two unaffected siblings born to healthy parents. Exome sequencing revealed that both twins carried a novel de novo mutation (p.A39E) in the GAL gene. Functional analysis revealed that the p.A39E mutant showed antagonistic activity against galanin receptor 1 (GalR1)-mediated response, and decreased binding affinity and reduced agonist properties for GalR2. These findings suggest that the p.A39E mutant could impair galanin signaling in the hippocampus, leading to increased glutamatergic excitation and ultimately to TLE. In a cohort of 582 cases, we did not observe any pathogenic mutations indicating that mutations in GAL are a rare cause of TLE. The identification of a novel de novo mutation in a biologically-relevant candidate gene, coupled with functional evidence that the mutant protein disrupts galanin signaling, strongly supports GAL as the causal gene for the TLE in this family. Given the availability of galanin agonists which inhibit seizures, our findings could potentially have direct implications for the development of anti-epileptic treatment.

Simultaneous identification and prioritization of variants in familial, de novo, and somatic genetic disorders with VariantMaster.

There is increasing interest in clinical genetics pertaining to the utilization of high-throughput sequencing data for accurate diagnoses of monogenic diseases. Moreover, massive whole-exome sequencing of tumors has provided significant advances in the understanding of cancer development through the recognition of somatic driver variants. To improve the identification of the variants from HTS, we developed VariantMaster, an original program that accurately and efficiently extracts causative variants in familial and sporadic genetic diseases. The algorithm takes into account predicted variants (SNPs and indels) in affected individuals or tumor samples and utilizes the row (BAM) data to robustly estimate the conditional probability of segregation in a family, as well as the probability of it being de novo or somatic. In familial cases, various modes of inheritance are considered: X-linked, autosomal dominant, and recessive (homozygosity or compound heterozygosity). Moreover, VariantMaster integrates phenotypes and genotypes, and employs Annovar to produce additional information such as allelic frequencies in the general population and damaging scores to further reduce the number of putative variants. As a proof of concept, we successfully applied VariantMaster to identify (1) de novo mutations in a previously described data set, (2) causative variants in a rare Mendelian genetic disease, and (3) known and new "driver" mutations in previously reported cancer data sets. Our results demonstrate that VariantMaster is considerably more accurate in terms of precision and sensitivity compared with previously published algorithms.