Megabase Length Hypermutation Accompanies Human Structural Variation at 17p11.2.

DNA rearrangements resulting in human genome structural variants (SVs) are caused by diverse mutational mechanisms. We used long- and short-read sequencing technologies to investigate end products of de novo chromosome 17p11.2 rearrangements and query the molecular mechanisms underlying both recurrent and non-recurrent events. Evidence for an increased rate of clustered single-nucleotide variant (SNV) mutation in cis with non-recurrent rearrangements was found. Indel and SNV formation are associated with both copy-number gains and losses of 17p11.2, occur up to ∼1 Mb away from the breakpoint junctions, and favor C > G transversion substitutions; results suggest that single-stranded DNA is formed during the genesis of the SV and provide compelling support for a microhomology-mediated break-induced replication (MMBIR) mechanism for SV formation. Our data show an additional mutational burden of MMBIR consisting of hypermutation confined to the locus and manifesting as SNVs and indels predominantly within genes.

An Organismal CNV Mutator Phenotype Restricted to Early Human Development.

De novo copy number variants (dnCNVs) arising at multiple loci in a personal genome have usually been considered to reflect cancer somatic genomic instabilities. We describe a multiple dnCNV (MdnCNV) phenomenon in which individuals with genomic disorders carry five to ten constitutional dnCNVs. These CNVs originate from independent formation incidences, are predominantly tandem duplications or complex gains, exhibit breakpoint junction features reminiscent of replicative repair, and show increased de novo point mutations flanking the rearrangement junctions. The active CNV mutation shower appears to be restricted to a transient perizygotic period. We propose that a defect in the CNV formation process is responsible for the "CNV-mutator state," and this state is dampened after early embryogenesis. The constitutional MdnCNV phenomenon resembles chromosomal instability in various cancers. Investigations of this phenomenon may provide unique access to understanding genomic disorders, structural variant mutagenesis, human evolution, and cancer biology.

A drosophila genetic resource of mutants to study mechanisms underlying human genetic diseases.

Invertebrate model systems are powerful tools for studying human disease owing to their genetic tractability and ease of screening. We conducted a mosaic genetic screen of lethal mutations on the Drosophila X chromosome to identify genes required for the development, function, and maintenance of the nervous system. We identified 165 genes, most of whose function has not been studied in vivo. In parallel, we investigated rare variant alleles in 1,929 human exomes from families with unsolved Mendelian disease. Genes that are essential in flies and have multiple human homologs were found to be likely to be associated with human diseases. Merging the human data sets with the fly genes allowed us to identify disease-associated mutations in six families and to provide insights into microcephaly associated with brain dysgenesis. This bidirectional synergism between fly genetics and human genomics facilitates the functional annotation of evolutionarily conserved genes involved in human health.

Human CLP1 mutations alter tRNA biogenesis, affecting both peripheral and central nervous system function.

CLP1 is a RNA kinase involved in tRNA splicing. Recently, CLP1 kinase-dead mice were shown to display a neuromuscular disorder with loss of motor neurons and muscle paralysis. Human genome analyses now identified a CLP1 homozygous missense mutation (p.R140H) in five unrelated families, leading to a loss of CLP1 interaction with the tRNA splicing endonuclease (TSEN) complex, largely reduced pre-tRNA cleavage activity, and accumulation of linear tRNA introns. The affected individuals develop severe motor-sensory defects, cortical dysgenesis, and microcephaly. Mice carrying kinase-dead CLP1 also displayed microcephaly and reduced cortical brain volume due to the enhanced cell death of neuronal progenitors that is associated with reduced numbers of cortical neurons. Our data elucidate a neurological syndrome defined by CLP1 mutations that impair tRNA splicing. Reduction of a founder mutation to homozygosity illustrates the importance of rare variations in disease and supports the clan genomics hypothesis.

TLR7 gain-of-function genetic variation causes human lupus.

Although circumstantial evidence supports enhanced Toll-like receptor 7 (TLR7) signalling as a mechanism of human systemic autoimmune disease1-7, evidence of lupus-causing TLR7 gene variants is lacking. Here we describe human systemic lupus erythematosus caused by a TLR7 gain-of-function variant. TLR7 is a sensor of viral RNA8,9 and binds to guanosine10-12. We identified a de novo, previously undescribed missense TLR7Y264H variant in a child with severe lupus and additional variants in other patients with lupus. The TLR7Y264H variant selectively increased sensing of guanosine and 2',3'-cGMP10-12, and was sufficient to cause lupus when introduced into mice. We show that enhanced TLR7 signalling drives aberrant survival of B cell receptor (BCR)-activated B cells, and in a cell-intrinsic manner, accumulation of CD11c+ age-associated B cells and germinal centre B cells. Follicular and extrafollicular helper T cells were also increased but these phenotypes were cell-extrinsic. Deficiency of MyD88 (an adaptor protein downstream of TLR7) rescued autoimmunity, aberrant B cell survival, and all cellular and serological phenotypes. Despite prominent spontaneous germinal-centre formation in Tlr7Y264H mice, autoimmunity was not ameliorated by germinal-centre deficiency, suggesting an extrafollicular origin of pathogenic B cells. We establish the importance of TLR7 and guanosine-containing self-ligands for human lupus pathogenesis, which paves the way for therapeutic TLR7 or MyD88 inhibition.

PSMD11 loss-of-function variants correlate with a neurobehavioral phenotype, obesity, and increased interferon response.

Primary proteasomopathies have recently emerged as a new class of rare early-onset neurodevelopmental disorders (NDDs) caused by pathogenic variants in the PSMB1, PSMC1, PSMC3, or PSMD12 proteasome genes. Proteasomes are large multi-subunit protein complexes that maintain cellular protein homeostasis by clearing ubiquitin-tagged damaged, misfolded, or unnecessary proteins. In this study, we have identified PSMD11 as an additional proteasome gene in which pathogenic variation is associated with an NDD-causing proteasomopathy. PSMD11 loss-of-function variants caused early-onset syndromic intellectual disability and neurodevelopmental delay with recurrent obesity in 10 unrelated children. Our findings demonstrate that the cognitive impairment observed in these individuals could be recapitulated in Drosophila melanogaster with depletion of the PMSD11 ortholog Rpn6, which compromised reversal learning. Our investigations in subject samples further revealed that PSMD11 loss of function resulted in impaired 26S proteasome assembly and the acquisition of a persistent type I interferon (IFN) gene signature, mediated by the integrated stress response (ISR) protein kinase R (PKR). In summary, these data identify PSMD11 as an additional member of the growing family of genes associated with neurodevelopmental proteasomopathies and provide insights into proteasomal biology in human health.

Variants in ZFX are associated with an X-linked neurodevelopmental disorder with recurrent facial gestalt.

Pathogenic variants in multiple genes on the X chromosome have been implicated in syndromic and non-syndromic intellectual disability disorders. ZFX on Xp22.11 encodes a transcription factor that has been linked to diverse processes including oncogenesis and development, but germline variants have not been characterized in association with disease. Here, we present clinical and molecular characterization of 18 individuals with germline ZFX variants. Exome or genome sequencing revealed 11 variants in 18 subjects (14 males and 4 females) from 16 unrelated families. Four missense variants were identified in 11 subjects, with seven truncation variants in the remaining individuals. Clinical findings included developmental delay/intellectual disability, behavioral abnormalities, hypotonia, and congenital anomalies. Overlapping and recurrent facial features were identified in all subjects, including thickening and medial broadening of eyebrows, variations in the shape of the face, external eye abnormalities, smooth and/or long philtrum, and ear abnormalities. Hyperparathyroidism was found in four families with missense variants, and enrichment of different tumor types was observed. In molecular studies, DNA-binding domain variants elicited differential expression of a small set of target genes relative to wild-type ZFX in cultured cells, suggesting a gain or loss of transcriptional activity. Additionally, a zebrafish model of ZFX loss displayed an altered behavioral phenotype, providing additional evidence for the functional significance of ZFX. Our clinical and experimental data support that variants in ZFX are associated with an X-linked intellectual disability syndrome characterized by a recurrent facial gestalt, neurocognitive and behavioral abnormalities, and an increased risk for congenital anomalies and hyperparathyroidism.

Core planar cell polarity genes VANGL1 and VANGL2 in predisposition to congenital vertebral malformations.

Congenital scoliosis (CS), affecting approximately 0.5 to 1 in 1,000 live births, is commonly caused by congenital vertebral malformations (CVMs) arising from aberrant somitogenesis or somite differentiation. While Wnt/ß-catenin signaling has been implicated in somite development, the function of Wnt/planar cell polarity (Wnt/PCP) signaling in this process remains unclear. Here, we investigated the role of Vangl1 and Vangl2 in vertebral development and found that their deletion causes vertebral anomalies resembling human CVMs. Analysis of exome sequencing data from multiethnic CS patients revealed a number of rare and deleterious variants in VANGL1 and VANGL2, many of which exhibited loss-of-function and dominant-negative effects. Zebrafish models confirmed the pathogenicity of these variants. Furthermore, we found that Vangl1 knock-in (p.R258H) mice exhibited vertebral malformations in a Vangl gene dose- and environment-dependent manner. Our findings highlight critical roles for PCP signaling in vertebral development and predisposition to CVMs in CS patients, providing insights into the molecular mechanisms underlying this disorder.

Bi-allelic variants in HMGCR cause an autosomal-recessive progressive limb-girdle muscular dystrophy.

Statins are a mainstay intervention for cardiovascular disease prevention, yet their use can cause rare severe myopathy. HMG-CoA reductase, an essential enzyme in the mevalonate pathway, is the target of statins. We identified nine individuals from five unrelated families with unexplained limb-girdle like muscular dystrophy and bi-allelic variants in HMGCR via clinical and research exome sequencing. The clinical features resembled other genetic causes of muscular dystrophy with incidental high CPK levels (>1,000 U/L), proximal muscle weakness, variable age of onset, and progression leading to impaired ambulation. Muscle biopsies in most affected individuals showed non-specific dystrophic changes with non-diagnostic immunohistochemistry. Molecular modeling analyses revealed variants to be destabilizing and affecting protein oligomerization. Protein activity studies using three variants (p.Asp623Asn, p.Tyr792Cys, and p.Arg443Gln) identified in affected individuals confirmed decreased enzymatic activity and reduced protein stability. In summary, we showed that individuals with bi-allelic amorphic (i.e., null and/or hypomorphic) variants in HMGCR display phenotypes that resemble non-genetic causes of myopathy involving this reductase. This study expands our knowledge regarding the mechanisms leading to muscular dystrophy through dysregulation of the mevalonate pathway, autoimmune myopathy, and statin-induced myopathy.

Bi-allelic variants in the ESAM tight-junction gene cause a neurodevelopmental disorder associated with fetal intracranial hemorrhage.

The blood-brain barrier (BBB) is an essential gatekeeper for the central nervous system and incidence of neurodevelopmental disorders (NDDs) is higher in infants with a history of intracerebral hemorrhage (ICH). We discovered a rare disease trait in thirteen individuals, including four fetuses, from eight unrelated families associated with homozygous loss-of-function variant alleles of ESAM which encodes an endothelial cell adhesion molecule. The c.115del (p.Arg39Glyfs∗33) variant, identified in six individuals from four independent families of Southeastern Anatolia, severely impaired the in vitro tubulogenic process of endothelial colony-forming cells, recapitulating previous evidence in null mice, and caused lack of ESAM expression in the capillary endothelial cells of damaged brain. Affected individuals with bi-allelic ESAM variants showed profound global developmental delay/unspecified intellectual disability, epilepsy, absent or severely delayed speech, varying degrees of spasticity, ventriculomegaly, and ICH/cerebral calcifications, the latter being also observed in the fetuses. Phenotypic traits observed in individuals with bi-allelic ESAM variants overlap very closely with other known conditions characterized by endothelial dysfunction due to mutation of genes encoding tight junction molecules. Our findings emphasize the role of brain endothelial dysfunction in NDDs and contribute to the expansion of an emerging group of diseases that we propose to rename as "tightjunctionopathies."

Bi-allelic SNAPC4 variants dysregulate global alternative splicing and lead to neuroregression and progressive spastic paraparesis.

The vast majority of human genes encode multiple isoforms through alternative splicing, and the temporal and spatial regulation of those isoforms is critical for organismal development and function. The spliceosome, which regulates and executes splicing reactions, is primarily composed of small nuclear ribonucleoproteins (snRNPs) that consist of small nuclear RNAs (snRNAs) and protein subunits. snRNA gene transcription is initiated by the snRNA-activating protein complex (SNAPc). Here, we report ten individuals, from eight families, with bi-allelic, deleterious SNAPC4 variants. SNAPC4 encoded one of the five SNAPc subunits that is critical for DNA binding. Most affected individuals presented with delayed motor development and developmental regression after the first year of life, followed by progressive spasticity that led to gait alterations, paraparesis, and oromotor dysfunction. Most individuals had cerebral, cerebellar, or basal ganglia volume loss by brain MRI. In the available cells from affected individuals, SNAPC4 abundance was decreased compared to unaffected controls, suggesting that the bi-allelic variants affect SNAPC4 accumulation. The depletion of SNAPC4 levels in HeLa cell lines via genomic editing led to decreased snRNA expression and global dysregulation of alternative splicing. Analysis of available fibroblasts from affected individuals showed decreased snRNA expression and global dysregulation of alternative splicing compared to unaffected cells. Altogether, these data suggest that these bi-allelic SNAPC4 variants result in loss of function and underlie the neuroregression and progressive spasticity in these affected individuals.

PLS3 missense variants affecting the actin-binding domains cause X-linked congenital diaphragmatic hernia and body-wall defects.

Congenital diaphragmatic hernia (CDH) is a relatively common and genetically heterogeneous structural birth defect associated with high mortality and morbidity. We describe eight unrelated families with an X-linked condition characterized by diaphragm defects, variable anterior body-wall anomalies, and/or facial dysmorphism. Using linkage analysis and exome or genome sequencing, we found that missense variants in plastin 3 (PLS3), a gene encoding an actin bundling protein, co-segregate with disease in all families. Loss-of-function variants in PLS3 have been previously associated with X-linked osteoporosis (MIM: 300910), so we used in silico protein modeling and a mouse model to address these seemingly disparate clinical phenotypes. The missense variants in individuals with CDH are located within the actin-binding domains of the protein but are not predicted to affect protein structure, whereas the variants in individuals with osteoporosis are predicted to result in loss of function. A mouse knockin model of a variant identified in one of the CDH-affected families, c.1497G>C (p.Trp499Cys), shows partial perinatal lethality and recapitulates the key findings of the human phenotype, including diaphragm and abdominal-wall defects. Both the mouse model and one adult human male with a CDH-associated PLS3 variant were observed to have increased rather than decreased bone mineral density. Together, these clinical and functional data in humans and mice reveal that specific missense variants affecting the actin-binding domains of PLS3 might have a gain-of-function effect and cause a Mendelian congenital disorder.

Bi-allelic TTI1 variants cause an autosomal-recessive neurodevelopmental disorder with microcephaly.

Telomere maintenance 2 (TELO2), Tel2 interacting protein 2 (TTI2), and Tel2 interacting protein 1 (TTI1) are the three components of the conserved Triple T (TTT) complex that modulates activity of phosphatidylinositol 3-kinase-related protein kinases (PIKKs), including mTOR, ATM, and ATR, by regulating the assembly of mTOR complex 1 (mTORC1). The TTT complex is essential for the expression, maturation, and stability of ATM and ATR in response to DNA damage. TELO2- and TTI2-related bi-allelic autosomal-recessive (AR) encephalopathies have been described in individuals with moderate to severe intellectual disability (ID), short stature, postnatal microcephaly, and a movement disorder (in the case of variants within TELO2). We present clinical, genomic, and functional data from 11 individuals in 9 unrelated families with bi-allelic variants in TTI1. All present with ID, and most with microcephaly, short stature, and a movement disorder. Functional studies performed in HEK293T cell lines and fibroblasts and lymphoblastoid cells derived from 4 unrelated individuals showed impairment of the TTT complex and of mTOR pathway activity which is improved by treatment with Rapamycin. Our data delineate a TTI1-related neurodevelopmental disorder and expand the group of disorders related to the TTT complex.

ß-Actin G342D as a Cause of NK Cell Deficiency Impairing Lytic Synapse Termination.

NK cell deficiency (NKD) occurs when an individual's major clinical immunodeficiency derives from abnormal NK cells and is associated with several genetic etiologies. Three categories of ß-actin-related diseases with over 60 ACTB (ß-actin) variants have previously been identified, none with a distinct NK cell phenotype. An individual with mild developmental delay, macrothrombocytopenia, and susceptibility to infections, molluscum contagiosum virus, and EBV-associated lymphoma had functional NKD for over a decade. A de novo ACTB variant encoding G342D ß-actin was identified and was consistent with the individual's developmental and platelet phenotype. This novel variant also was found to have direct impact in NK cells because its expression in the human NK cell line YTS (YTS-NKD) caused increased cell spreading in lytic immune synapses created on activating surfaces. YTS-NKD cells were able to degranulate and perform cytotoxicity, but they demonstrated defective serial killing because of prolonged conjugation to the killed target cell and thus were effectively unable to terminate lytic synapses. G342D ß-actin results in a novel, to our knowledge, mechanism of functional NKD via increased synaptic spreading and defective lytic synapse termination with resulting impaired serial killing, leading to overall reductions in NK cell cytotoxicity.

Clan genomics and the complex architecture of human disease.

Human diseases are caused by alleles that encompass the full range of variant types, from single-nucleotide changes to copy-number variants, and these variations span a broad frequency spectrum, from the very rare to the common. The picture emerging from analysis of whole-genome sequences, the 1000 Genomes Project pilot studies, and targeted genomic sequencing derived from very large sample sizes reveals an abundance of rare and private variants. One implication of this realization is that recent mutation may have a greater influence on disease susceptibility or protection than is conferred by variations that arose in distant ancestors.

Chromosome catastrophes involve replication mechanisms generating complex genomic rearrangements.

Complex genomic rearrangements (CGRs) consisting of two or more breakpoint junctions have been observed in genomic disorders. Recently, a chromosome catastrophe phenomenon termed chromothripsis, in which numerous genomic rearrangements are apparently acquired in one single catastrophic event, was described in multiple cancers. Here, we show that constitutionally acquired CGRs share similarities with cancer chromothripsis. In the 17 CGR cases investigated, we observed localization and multiple copy number changes including deletions, duplications, and/or triplications, as well as extensive translocations and inversions. Genomic rearrangements involved varied in size and complexities; in one case, array comparative genomic hybridization revealed 18 copy number changes. Breakpoint sequencing identified characteristic features, including small templated insertions at breakpoints and microhomology at breakpoint junctions, which have been attributed to replicative processes. The resemblance between CGR and chromothripsis suggests similar mechanistic underpinnings. Such chromosome catastrophic events appear to reflect basic DNA metabolism operative throughout an organism's life cycle.

HMZDupFinder: a robust computational approach for detecting intragenic homozygous duplications from exome sequencing data.

Homozygous duplications contribute to genetic disease by altering gene dosage or disrupting gene regulation and can be more deleterious to organismal biology than heterozygous duplications. Intragenic exonic duplications can result in loss-of-function (LoF) or gain-of-function (GoF) alleles that when homozygosed, i.e. brought to homozygous state at a locus by identity by descent or state, could potentially result in autosomal recessive (AR) rare disease traits. However, the detection and functional interpretation of homozygous duplications from exome sequencing data remains a challenge. We developed a framework algorithm, HMZDupFinder, that is designed to detect exonic homozygous duplications from exome sequencing (ES) data. The HMZDupFinder algorithm can efficiently process large datasets and accurately identifies small intragenic duplications, including those associated with rare disease traits. HMZDupFinder called 965 homozygous duplications with three or less exons from 8,707 ES with a recall rate of 70.9% and a precision of 16.1%. We experimentally confirmed 8/10 rare homozygous duplications. Pathogenicity assessment of these copy number variant alleles allowed clinical genomics contextualization for three homozygous duplications alleles, including two affecting known OMIM disease genes EDAR (MIM# 224900), TNNT1(MIM# 605355), and one variant in a novel candidate disease gene: PAAF1.

Hemizygous variants in protein phosphatase 1 regulatory subunit 3F (PPP1R3F) are associated with a neurodevelopmental disorder characterized by developmental delay, intellectual disability and autistic features.

Protein phosphatase 1 regulatory subunit 3F (PPP1R3F) is a member of the glycogen targeting subunits (GTSs), which belong to the large group of regulatory subunits of protein phosphatase 1 (PP1), a major eukaryotic serine/threonine protein phosphatase that regulates diverse cellular processes. Here, we describe the identification of hemizygous variants in PPP1R3F associated with a novel X-linked recessive neurodevelopmental disorder in 13 unrelated individuals. This disorder is characterized by developmental delay, mild intellectual disability, neurobehavioral issues such as autism spectrum disorder, seizures and other neurological findings including tone, gait and cerebellar abnormalities. PPP1R3F variants segregated with disease in affected hemizygous males that inherited the variants from their heterozygous carrier mothers. We show that PPP1R3F is predominantly expressed in brain astrocytes and localizes to the endoplasmic reticulum in cells. Glycogen content in PPP1R3F knockout astrocytoma cells appears to be more sensitive to fluxes in extracellular glucose levels than in wild-type cells, suggesting that PPP1R3F functions in maintaining steady brain glycogen levels under changing glucose conditions. We performed functional studies on nine of the identified variants and observed defects in PP1 binding, protein stability, subcellular localization and regulation of glycogen metabolism in most of them. Collectively, the genetic and molecular data indicate that deleterious variants in PPP1R3F are associated with a new X-linked disorder of glycogen metabolism, highlighting the critical role of GTSs in neurological development. This research expands our understanding of neurodevelopmental disorders and the role of PP1 in brain development and proper function.

Biology in balance: human diploid genome integrity, gene dosage, and genomic medicine.

The path to completion of the functional annotation of the haploid human genome reference build, exploration of the clan genomics hypothesis, understanding human gene and genome functional biology, and gene genome and organismal evolution, is in reach.