Evidence for the additivity of rare and common variant burden throughout the spectrum of intellectual disability.

Intellectual disability (ID) is a common disorder, yet there is a wide spectrum of impairment from mild to profoundly affected individuals. Mild ID is seen as the low extreme of the general distribution of intelligence, while severe ID is often seen as a monogenic disorder caused by rare, pathogenic, highly penetrant variants. To investigate the genetic factors influencing mild and severe ID, we evaluated rare and common variation in the Northern Finland Intellectual Disability cohort (n = 1096 ID patients), a cohort with a high percentage of mild ID (n = 550) and from a population bottleneck enriched in rare, damaging variation. Despite this enrichment, we found only a small percentage of ID was due to recessive Finnish-enriched variants (0.5%). A larger proportion was linked to dominant variation, with a significant burden of rare, damaging variation in both mild and severe ID. This rare variant burden was enriched in more severe ID (p = 2.4e-4), patients without a relative with ID (p = 4.76e-4), and in those with features associated with monogenic disorders. We also found a significant burden of common variants associated with decreased cognitive function, with no difference between mild and more severe ID. When we included common and rare variants in a joint model, the rare and common variants had additive effects in both mild and severe ID. A multimodel inference approach also found that common and rare variants together best explained ID status (ΔAIC = 16.8, ΔBIC = 10.2). Overall, we report evidence for the additivity of rare and common variant burden throughout the spectrum of intellectual disability.

Novel patients with NHLRC2 variants expand the phenotypic spectrum of FINCA disease.

Purpose: FINCA disease (Fibrosis, Neurodegeneration and Cerebral Angiomatosis, OMIM 618278) is an infantile-onset neurodevelopmental and multiorgan disease. Since our initial report in 2018, additional patients have been described. FINCA is the first human disease caused by recessive variants in the highly conserved NHLRC2 gene. Our previous studies have shown that Nhlrc2-null mouse embryos die during gastrulation, indicating the essential role of the protein in embryonic development. Defect in NHLRC2 leads to cerebral neurodegeneration and severe pulmonary, hepatic and cardiac fibrosis. Despite having a structure suggestive of an enzymatic role and the clinical importance of NHLRC2 in multiple organs, the specific physiological role of the protein is unknown. Methods: The clinical histories of five novel FINCA patients diagnosed with whole exome sequencing were reviewed. Segregation analysis of the biallelic, potentially pathogenic NHLRC2 variants was performed using Sanger sequencing. Studies on neuropathology and NHLRC2 expression in different brain regions were performed on autopsy samples of three previously described deceased FINCA patients. Results: One patient was homozygous for the pathogenic variant c.442G > T, while the other four were compound heterozygous for this variant and two other pathogenic NHLRC2 gene variants. All five patients presented with multiorgan dysfunction with neurodevelopmental delay, recurrent infections and macrocytic anemia as key features. Interstitial lung disease was pronounced in infancy but often stabilized. Autopsy samples revealed widespread, albeit at a lower intensity than the control, NHLRC2 expression in the brain. Conclusion: This report expands on the characteristic clinical features of FINCA disease. Presentation is typically in infancy, and although patients can live to late adulthood, the key clinical and histopathological features are fibrosis, infection susceptibility/immunodeficiency/intellectual disability, neurodevelopmental disorder/neurodegeneration and chronic anemia/cerebral angiomatosis (hence the acronym FINCA) that enable an early diagnosis confirmed by genetic investigations.

A novel variant in SMG9 causes intellectual disability, confirming a role for nonsense-mediated decay components in neurocognitive development.

Biallelic loss-of-function variants in the SMG9 gene, encoding a regulatory subunit of the mRNA nonsense-mediated decay (NMD) machinery, are reported to cause heart and brain malformation syndrome. Here we report five patients from three unrelated families with intellectual disability (ID) and a novel pathogenic SMG9 c.551 T > C p.(Val184Ala) homozygous missense variant, identified using exome sequencing. Sanger sequencing confirmed recessive segregation in each family. SMG9 c.551T > C p.(Val184Ala) is most likely an autozygous variant identical by descent. Characteristic clinical findings in patients were mild to moderate ID, intention tremor, pyramidal signs, dyspraxia, and ocular manifestations. We used RNA sequencing of patients and age- and sex-matched healthy controls to assess the effect of the variant. RNA sequencing revealed that the SMG9 c.551T > C variant did not affect the splicing or expression level of SMG9 gene products, and allele-specific expression analysis did not provide evidence that the nonsense mRNA-induced NMD was affected. Differential gene expression analysis identified prevalent upregulation of genes in patients, including the genes SMOX, OSBP2, GPX3, and ZNF155. These findings suggest that normal SMG9 function may be involved in transcriptional regulation without affecting nonsense mRNA-induced NMD. In conclusion, we demonstrate that the SMG9 c.551T > C missense variant causes a neurodevelopmental disorder and impacts gene expression. NMD components have roles beyond aberrant mRNA degradation that are crucial for neurocognitive development.

Correction: Phenotypic spectrum associated with a CRADD founder variant underlying frontotemporal predominant pachygyria in the Finnish population.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Phenotypic spectrum associated with a CRADD founder variant underlying frontotemporal predominant pachygyria in the Finnish population.

Intellectual disability (ID), megalencephaly, frontal predominant pachygyria, and seizures, previously called "thin" lissencephaly, are reported to be caused by recessive variants in CRADD. Among five families of different ethnicities identified, one homozygous missense variant, c.509G>A p.(Arg170His), was of Finnish ancestry. Here we report on the phenotypic variability associated for this potential CRADD founder variant in 22 Finnish individuals. Exome sequencing was used to identify candidate genes in Finnish patients presenting with ID. Targeted Sanger sequencing and restriction enzyme analysis were applied to screen for the c.509G>A CRADD variant in cohorts from Finland. Detailed phenotyping and genealogical studies were performed. Twenty two patients were identified with the c.509G>A p.(Arg170His) homozygous variant in CRADD. The majority of the ancestors originated from Northeastern Finland indicating a founder effect. The hallmark of the disease is frontotemporal predominant pachygyria with mild cortical thickening. All patients show ID of variable severity. Aggressive behavior was found in nearly half of the patients, EEG abnormalities in five patients and megalencephaly in three patients. This study provides detailed data about the phenotypic spectrum of patients with lissencephaly due to a CRADD variant that affects function. High inter- and intrafamilial phenotypic heterogeneity was identified in patients with pachygyria caused by the homozygous CRADD founder variant. The phenotype variability suggests that additional genetic and/or environmental factors play a role in the clinical presentation. Since frontotemporal pachygyria is the hallmark of the disease, brain imaging studies are essential to support the molecular diagnosis for individuals with ID and a CRADD variant.

Contribution of rare and common variants to intellectual disability in a sub-isolate of Northern Finland.

The contribution of de novo variants in severe intellectual disability (ID) has been extensively studied whereas the genetics of mild ID has been less characterized. To elucidate the genetics of milder ID we studied 442 ID patients enriched for mild ID (>50%) from a population isolate of Finland. Using exome sequencing, we show that rare damaging variants in known ID genes are observed significantly more often in severe (27%) than in mild ID (13%) patients. We further observe a significant enrichment of functional variants in genes not yet associated with ID (OR: 2.1). We show that a common variant polygenic risk significantly contributes to ID. The heritability explained by polygenic risk score is the highest for educational attainment (EDU) in mild ID (2.2%) but lower for more severe ID (0.6%). Finally, we identify a Finland enriched homozygote variant in the CRADD ID associated gene.

Quantifying the Impact of Rare and Ultra-rare Coding Variation across the Phenotypic Spectrum.

There is a limited understanding about the impact of rare protein-truncating variants across multiple phenotypes. We explore the impact of this class of variants on 13 quantitative traits and 10 diseases using whole-exome sequencing data from 100,296 individuals. Protein-truncating variants in genes intolerant to this class of mutations increased risk of autism, schizophrenia, bipolar disorder, intellectual disability, and ADHD. In individuals without these disorders, there was an association with shorter height, lower education, increased hospitalization, and reduced age at enrollment. Gene sets implicated from GWASs did not show a significant protein-truncating variants burden beyond what was captured by established Mendelian genes. In conclusion, we provide a thorough investigation of the impact of rare deleterious coding variants on complex traits, suggesting widespread pleiotropic risk.

Rare loss-of-function variants in SETD1A are associated with schizophrenia and developmental disorders.

By analyzing the whole-exome sequences of 4,264 schizophrenia cases, 9,343 controls and 1,077 trios, we identified a genome-wide significant association between rare loss-of-function (LoF) variants in SETD1A and risk for schizophrenia (P = 3.3 × 10(-9)). We found only two heterozygous LoF variants in 45,376 exomes from individuals without a neuropsychiatric diagnosis, indicating that SETD1A is substantially depleted of LoF variants in the general population. Seven of the ten individuals with schizophrenia carrying SETD1A LoF variants also had learning difficulties. We further identified four SETD1A LoF carriers among 4,281 children with severe developmental disorders and two more carriers in an independent sample of 5,720 Finnish exomes, both with notable neuropsychiatric phenotypes. Together, our observations indicate that LoF variants in SETD1A cause a range of neurodevelopmental disorders, including schizophrenia. Combining these data with previous common variant evidence, we suggest that epigenetic dysregulation, specifically in the histone H3K4 methylation pathway, is an important mechanism in the pathogenesis of schizophrenia.

Consortium for osteogenesis imperfecta mutations in the helical domain of type I collagen: regions rich in lethal mutations align with collagen binding sites for integrins and proteoglycans.

Osteogenesis imperfecta (OI) is a generalized disorder of connective tissue characterized by fragile bones and easy susceptibility to fracture. Most cases of OI are caused by mutations in type I collagen. We have identified and assembled structural mutations in type I collagen genes (COL1A1 and COL1A2, encoding the proalpha1(I) and proalpha2(I) chains, respectively) that result in OI. Quantitative defects causing type I OI were not included. Of these 832 independent mutations, 682 result in substitution for glycine residues in the triple helical domain of the encoded protein and 150 alter splice sites. Distinct genotype-phenotype relationships emerge for each chain. One-third of the mutations that result in glycine substitutions in alpha1(I) are lethal, especially when the substituting residues are charged or have a branched side chain. Substitutions in the first 200 residues are nonlethal and have variable outcome thereafter, unrelated to folding or helix stability domains. Two exclusively lethal regions (helix positions 691-823 and 910-964) align with major ligand binding regions (MLBRs), suggesting crucial interactions of collagen monomers or fibrils with integrins, matrix metalloproteinases (MMPs), fibronectin, and cartilage oligomeric matrix protein (COMP). Mutations in COL1A2 are predominantly nonlethal (80%). Lethal substitutions are located in eight regularly spaced clusters along the chain, supporting a regional model. The lethal regions align with proteoglycan binding sites along the fibril, suggesting a role in fibril-matrix interactions. Recurrences at the same site in alpha2(I) are generally concordant for outcome, unlike alpha1(I). Splice site mutations comprise 20% of helical mutations identified in OI patients, and may lead to exon skipping, intron inclusion, or the activation of cryptic splice sites. Splice site mutations in COL1A1 are rarely lethal; they often lead to frameshifts and the mild type I phenotype. In alpha2(I), lethal exon skipping events are located in the carboxyl half of the chain. Our data on genotype-phenotype relationships indicate that the two collagen chains play very different roles in matrix integrity and that phenotype depends on intracellular and extracellular events.

Univariate and bivariate variance component linkage analysis of a whole-genome scan for loci contributing to bone mineral density.

Osteoporosis is a common condition characterized by reduced skeletal strength and increased susceptibility to fracture. The single major risk factor for osteoporosis is low bone mineral density (BMD) and strong evidence exists that genetic factors are in part responsible for an individual's BMD. A cohort of 40 multiplex Caucasian families selected through a proband with osteoporosis was genotyped for microsatellite markers spaced at an average of 10 cM, and linkage to femoral neck (FN), lumbar spine (LS) and trochanter (TR) BMD was analyzed using univariate and bivariate variance component linkage analysis. Maximum univariate multipoint lod-scores were 2.87 on chromosome 1p36 for FN BMD, 1.89 on 6q27 for TR BMD, and 2.15 on 7p15 for LS BMD. Results of bivariate linkage analysis were highly correlated with those of the univariate analysis, although generally less significant, suggesting the possibility that some of these susceptibility loci may exert pleiotropic effects on multiple skeletal sites.

Childhood-onset osteoarthritis, tall stature, and sensorineural hearing loss associated with Arg75-Cys mutation in procollagen type II gene (COL2A1).

OBJECTIVE: To define the clinical, radiologic, and molecular genetic characteristics of a family with early progressive osteoarthritis mimicking childhood rheumatoid arthritis, Scheuermann-like changes of the spine, tall stature, short 3 and 4 metatarsals, and moderate sensorineural hearing loss. METHODS: We describe a 22-year-old woman and her 54-year-old mother with early progressive osteoarthritis mimicking childhood rheumatoid arthritis. The index case, her mother, and 3 other family members underwent a physical examination, anthropometric measurements, and radiologic studies. Their DNA was sequenced for the procollagen type II (COL2A1) gene. RESULTS: Mild scoliosis was noticed in the proband at the age of 6 years, and at the age of 7 years large Schmorl's nodes were found in the vertebrae L1-2. At the age of 11 years, changes resembling Scheuermann's disease were seen, mostly in the thoracic vertebrae. At the same age, she began to have arthralgia in the weight-bearing joints and osteoarthritis progressed fast, necessitating a hip prosthesis at the age of 18 years. The proband and her mother had bilateral sensorineural hearing loss of moderate degree. Both mother and daughter had an Arg75-Cys mutation in the COL2A1 gene. CONCLUSION: This family is the fourth example of the Arg75-Cys mutation in the COL2A1 gene, which appears to lead to a clearly recognizable phenotype. The finding suggests that sensorineural hearing loss may be a part of this syndrome.

Lack of correlation between the type of COL1A1 or COL1A2 mutation and hearing loss in osteogenesis imperfecta patients.

Osteogenesis imperfecta (OI) is caused by mutations in COL1A1 and COL1A2 that code for the alpha1 and alpha2 chains of type I collagen. Phenotypes correlate with the mutation types in that COL1A1 null mutations lead to OI type I, and structural mutations in alpha1(I) or alpha2(I) lead to more severe OI types (II-IV). However, correlative analysis between mutation types and OI associated hearing loss has not been previously performed. A total of 54 Finnish OI patients with previously diagnosed hearing loss or age 35 or more years were analyzed here for mutations in COL1A1 or COL1A2. Altogether 49 mutations were identified, of which 41 were novel. The 49 mutations represented the molecular genetic background of 41.1% of the Finnish OI population. A total of 38 mutations were in COL1A1 and 11 were in COL1A2. Of these, 16 were glycine substitutions and 16 were splicing mutations in alpha1(I) or alpha2(I). In addition, 17 null allele mutations were detected in COL1A1. A total of 32 patients (65.3%) with a mutation had hearing loss. That is slightly more than in our previous population study on Finnish adults with OI (57.9%). The association between the mutation types and OI type was statistically evident. Patients with COL1A1 mutations more frequently had blue scleras than those with COL1A2 mutations. In addition, patients with COL1A2 mutations tended to be shorter than those with COL1A1 mutations. However, no correlation was found between the mutated gene or mutation type and hearing pattern. These results suggest that the basis of hearing loss in OI is complex, and it is a result of multifactorial, still unknown genetic effects.

Autosomal dominant pseudohypoparathyroidism type Ib is associated with a heterozygous microdeletion that likely disrupts a putative imprinting control element of GNAS.

Patients with pseudohypoparathyroidism type Ib (PHP-Ib) have hypocalcemia and hyperphosphatemia due to renal parathyroid hormone (PTH) resistance, but lack physical features of Albright hereditary osteodystrophy. PHP-Ib is thus distinct from PHP-Ia, which is caused by mutations in the GNAS exons encoding the G protein alpha subunit. However, an imprinted autosomal dominant form of PHP-Ib (AD-PHP-Ib) has been mapped to a region of chromosome 20q13.3 containing GNAS. Furthermore, loss of methylation at a differentially methylated region (DMR) of this locus, exon A/B, has been observed thus far in all investigated sporadic PHP-Ib cases and the affected members of multiple AD-PHP-Ib kindreds. We now report that affected members and obligate gene carriers of 12 unrelated AD-PHP-Ib kindreds and four apparently sporadic PHP-Ib patients, but not healthy controls, have a heterozygous approximately 3-kb microdeletion located approximately 220 kb centromeric of GNAS exon A/B. The deleted region, which is flanked by two direct repeats, includes three exons of STX16, the gene encoding syntaxin-16, for which no evidence of imprinting could be found. Affected individuals carrying the microdeletion show loss of exon A/B methylation but no epigenetic abnormalities at other GNAS DMRs. We therefore postulate that this microdeletion disrupts a putative cis-acting element required for methylation at exon A/B, and that this genetic defect underlies the renal PTH resistance in AD-PHP-Ib.

Posterior chorioretinal atrophy and vitreous phenotype in a family with Stickler syndrome from a mutation in the COL2A1 gene.

PURPOSE: To report posterior chorioretinal atrophy (PCRA) and correlate the vitreous phenotype with inheritance of the disease mutation in a family with vitreoretinal dystrophy. DESIGN: Prospective observational case series. METHODS: Twenty-four members of a family with 14 affected individuals were examined, and genetic linkage analysis was performed at the COL2A1, COL11A1, and Wagner disease loci. The vitreous phenotype was prospectively graded as optically empty with retrolenticular membrane, fibrillar, or normal. Ocular ultrasonography and optical coherence tomography (OCT) were performed on selected individuals to study the vitreous structure and vitreoretinal interface. RESULTS: The 6-year-old proband had PCRA and optically empty vitreous without systemic features, suggestive of Wagner disease. The family history was negative for systemic disease, except for one cousin with cleft palate. However, when examined, clinical features of the 14 affected subjects included 5 with small chin, 4 with at least submucosal cleft palate, and 9 with a myopic refractive error greater than 5 diopters. Lens opacity or previous cataract extraction was found in 13 family members. All affected individuals in whom the vitreous could be examined had an optically empty vitreous with retrolental membrane. Posterior chorioretinal atrophy was found in eight of the affected subjects. The finding was not limited to highly myopic subjects, nor did all the high myopes have PCRA. Ultrasonography and OCT revealed vitreous adherent to the retina, but without apparent retinal distortion or edema of the macula. Significant linkage was established to the COL2A1 locus; the other loci were excluded. A single nucleotide insertion mutation (c.2012 2013insC) was identified in exon 34, leading to a downstream premature stop codon in the COL2A1 gene. CONCLUSIONS: Although posterior chorioretinal atrophy and vitreoretinal degeneration have been classically associated with Wagner disease, we demonstrate its presence in a family with typical Stickler syndrome. On the basis of clinical, ultrasonographic, and OCT studies, the etiology of PCRA in this family does not seem to be attributable to vitreomacular traction or myopia. The vitreous findings in this large family confirm reports that mutations in the COL2A1 gene lead to the optically empty vitreous with retrolenticular membrane phenotype.

Radial perivascular retinal degeneration: a key to the clinical diagnosis of an ocular variant of Stickler syndrome with minimal or no systemic manifestations.

PURPOSE: To identify the genetic defect and present the ocular and extraocular findings in a large pedigree of predominantly ocular Stickler syndrome. DESIGN: Observational case series. METHODS: An eight-generation pedigree with hereditary retinal detachments was retrospectively and prospectively studied. Clinical information was obtained by medical records, telephone interviews, medical questionnaires, detailed ophthalmologic examinations, physical examinations, and personal observations. Linkage analysis of the COL2A1 gene was performed on 21 family members, and mutation analysis was performed on three family members. RESULTS: The pedigree consisted of 100 affected individuals. The ocular findings, frequently bilateral, consisted of radial perivascular retinal degeneration (RPRD) (100%), vitreous syneresis (100%), high myopia (76%), retinal detachment (65%), presenile cataract development (occurring before 50 years of age; 78%), and glaucoma (18%). Most (70%) of the retinal detachments occurred between 4 and 18 years of age. Extraocular manifestations, characteristic for Stickler syndrome, were detected in only four of 100 (4%) affected individuals. Linkage analysis with COL2A1 flanking markers showed evidence for linkage to the COL2A1 locus. The COL2A1 gene analysis identified a mutation converting a codon TGC for cysteine(86) to a premature termination codon in the alternatively spliced exon 2. CONCLUSIONS: A variant of Stickler syndrome, caused by mutations in exon 2 of COL2A1, may present in families with all of the ocular findings and no clinically identifiable extraocular findings associated with Stickler syndrome. The predominant ocular findings are a congenitally abnormal vitreous and an acquired radial perivascular retinal degeneration that may lead to complicated childhood and adult retinal detachment.

Mapping the ligand-binding sites and disease-associated mutations on the most abundant protein in the human, type I collagen.

Type I collagen is the most abundant protein in humans, and it helps to maintain the integrity of many tissues via its interactions with cell surfaces, other extracellular matrix molecules, and growth and differentiation factors. Nearly 50 molecules have been found to interact with type I collagen, and for about half of them, binding sites on this collagen have been elucidated. In addition, over 300 mutations in type I collagen associated with human connective tissue disorders have been described. However, the spatial relationships between the known ligand-binding sites and mutation positions have not been examined. To this end, here we have created a map of type I collagen that includes all of its ligand-binding sites and mutations. The map reveals the existence of several hot spots for ligand interactions on type I collagen and that most of the binding sites locate to its C-terminal half. Moreover, on the collagen fibril some potentially relevant relationships between binding sites were observed including the following: fibronectin- and certain integrin-binding regions are near neighbors, which may mechanistically relate to fibronectin-dependent cell-collagen attachment; proteoglycan binding may potentially impact upon collagen fibrillogenesis, cell-collagen attachment, and collagen glycation seen in diabetes and aging; and mutations associated with osteogenesis imperfecta and other disorders show apparently nonrandom distribution patterns within both the monomer and fibril, implying that mutation positions correlate with disease phenotype. These and other observations presented here may provide novel insights into evaluating type I collagen functions and the relationships between its binding partners and mutations.

Procollagen with skipping of alpha 1(I) exon 41 has lower binding affinity for alpha 1(I) C-telopeptide, impaired in vitro fibrillogenesis, and altered fibril morphology.

Previous in vitro data on type I collagen self-assembly into fibrils suggested that the amino acid 776-796 region of the alpha1(I) chain is crucial for fibril formation because it serves as the recognition site for the telopeptide of a docking collagen monomer. We used a natural collagen mutation with a deletion of amino acids 766-801 to confirm the importance of this region for collagen fibril formation. The proband has type III osteogenesis imperfecta and is heterozygous for a COL1A1 IVS 41 A(+4) --> C substitution. The intronic mutation causes splicing of exon 41, confirmed by sequencing of normal and shorter reverse transcriptase-PCR products. Reverse transcriptase-PCR using RNA from proband dermal fibroblasts and clonal cell lines showed the mutant cDNA was about 15% of total alpha1(I) cDNA. The mutant transcript is translated; structurally abnormal alpha chains are demonstrated in the cell layer of proband fibroblasts by SDS-urea-PAGE. The proportion of mutant chains in the secreted procollagen was determined to be 10% by resistance to digestion with MMP-1, since chains lacking exon 41 are missing the vertebral collagenase cleavage site. Secreted proband collagen was used for analysis of kinetics of binding of alpha1(I) C-telopeptide using an optical biosensor. Telopeptide had slower association and faster dissociation from proband than from normal collagen. Purified proband pC-collagen was used to study fibril formation. The presence of the mutant molecules decreases the rate of fibril formation. The fibrils formed in the presence of 10-15% mutant molecules have strikingly increased length compared with normal collagen, but are well organized, as demonstrated by D-periodicity. These results suggest that some collagen molecules containing the mutant chain are incorporated into fibrils and that the absence of the telopeptide binding region from even a small portion of the monomers interferes with fibril growth. Both abnormal fibrils and slower remodeling may contribute to the severe phenotype.in