Molecular consequences of PQBP1 deficiency, involved in the X-linked Renpenning syndrome.

Mutations in the PQBP1 gene (polyglutamine-binding protein-1) are responsible for a syndromic X-linked form of neurodevelopmental disorder (XL-NDD) with intellectual disability (ID), named Renpenning syndrome. PQBP1 encodes a protein involved in transcriptional and post-transcriptional regulation of gene expression. To investigate the consequences of PQBP1 loss, we used RNA interference to knock-down (KD) PQBP1 in human neural stem cells (hNSC). We observed a decrease of cell proliferation, as well as the deregulation of the expression of 58 genes, comprising genes encoding proteins associated with neurodegenerative diseases, playing a role in mRNA regulation or involved in innate immunity. We also observed an enrichment of genes involved in other forms of NDD (CELF2, APC2, etc). In particular, we identified an increase of a non-canonical isoform of another XL-NDD gene, UPF3B, an actor of nonsense mRNA mediated decay (NMD). This isoform encodes a shorter protein (UPF3B_S) deprived from the domains binding NMD effectors, however no notable change in NMD was observed after PQBP1-KD in fibroblasts containing a premature termination codon. We showed that short non-canonical and long canonical UPF3B isoforms have different interactomes, suggesting they could play distinct roles. The link between PQBP1 loss and increase of UPF3B_S expression was confirmed in mRNA obtained from patients with pathogenic variants in PQBP1, particularly pronounced for truncating variants and missense variants located in the C-terminal domain. We therefore used it as a molecular marker of Renpenning syndrome, to test the pathogenicity of variants of uncertain clinical significance identified in PQPB1 in individuals with NDD, using patient blood mRNA and HeLa cells expressing wild-type or mutant PQBP1 cDNA. We showed that these different approaches were efficient to prove a functional effect of variants in the C-terminal domain of the protein. In conclusion, our study provided information on the pathological mechanisms involved in Renpenning syndrome, but also allowed the identification of a biomarker of PQBP1 deficiency useful to test variant effect.

Single-cell, whole-embryo phenotyping of mammalian developmental disorders.

Mouse models are a critical tool for studying human diseases, particularly developmental disorders1. However, conventional approaches for phenotyping may fail to detect subtle defects throughout the developing mouse2. Here we set out to establish single-cell RNA sequencing of the whole embryo as a scalable platform for the systematic phenotyping of mouse genetic models. We applied combinatorial indexing-based single-cell RNA sequencing3 to profile 101 embryos of 22 mutant and 4 wild-type genotypes at embryonic day 13.5, altogether profiling more than 1.6 million nuclei. The 22 mutants represent a range of anticipated phenotypic severities, from established multisystem disorders to deletions of individual regulatory regions4,5. We developed and applied several analytical frameworks for detecting differences in composition and/or gene expression across 52 cell types or trajectories. Some mutants exhibit changes in dozens of trajectories whereas others exhibit changes in only a few cell types. We also identify differences between widely used wild-type strains, compare phenotyping of gain- versus loss-of-function mutants and characterize deletions of topological associating domain boundaries. Notably, some changes are shared among mutants, suggesting that developmental pleiotropy might be 'decomposable' through further scaling of this approach. Overall, our findings show how single-cell profiling of whole embryos can enable the systematic molecular and cellular phenotypic characterization of mouse mutants with unprecedented breadth and resolution.

Aberrant phase separation and nucleolar dysfunction in rare genetic diseases.

Thousands of genetic variants in protein-coding genes have been linked to disease. However, the functional impact of most variants is unknown as they occur within intrinsically disordered protein regions that have poorly defined functions1-3. Intrinsically disordered regions can mediate phase separation and the formation of biomolecular condensates, such as the nucleolus4,5. This suggests that mutations in disordered proteins may alter condensate properties and function6-8. Here we show that a subset of disease-associated variants in disordered regions alter phase separation, cause mispartitioning into the nucleolus and disrupt nucleolar function. We discover de novo frameshift variants in HMGB1 that cause brachyphalangy, polydactyly and tibial aplasia syndrome, a rare complex malformation syndrome. The frameshifts replace the intrinsically disordered acidic tail of HMGB1 with an arginine-rich basic tail. The mutant tail alters HMGB1 phase separation, enhances its partitioning into the nucleolus and causes nucleolar dysfunction. We built a catalogue of more than 200,000 variants in disordered carboxy-terminal tails and identified more than 600 frameshifts that create arginine-rich basic tails in transcription factors and other proteins. For 12 out of the 13 disease-associated variants tested, the mutation enhanced partitioning into the nucleolus, and several variants altered rRNA biogenesis. These data identify the cause of a rare complex syndrome and suggest that a large number of genetic variants may dysregulate nucleoli and other biomolecular condensates in humans.

Large-Scale Functional Assessment of Genes Involved in Rare Diseases with Intellectual Disabilities Unravels Unique Developmental and Behaviour Profiles in Mouse Models.

Major progress has been made over the last decade in identifying novel genes involved in neurodevelopmental disorders, although the task of elucidating their corresponding molecular and pathophysiological mechanisms, which are an essential prerequisite for developing therapies, has fallen far behind. We selected 45 genes for intellectual disabilities to generate and characterize mouse models. Thirty-nine of them were based on the frequency of pathogenic variants in patients and literature reports, with several corresponding to de novo variants, and six other candidate genes. We used an extensive screen covering the development and adult stages, focusing specifically on behaviour and cognition to assess a wide range of functions and their pathologies, ranging from basic neurological reflexes to cognitive abilities. A heatmap of behaviour phenotypes was established, together with the results of selected mutants. Overall, three main classes of mutant lines were identified based on activity phenotypes, with which other motor or cognitive deficits were associated. These data showed the heterogeneity of phenotypes between mutation types, recapitulating several human features, and emphasizing the importance of such systematic approaches for both deciphering genetic etiological causes of ID and autism spectrum disorders, and for building appropriate therapeutic strategies.

Systematic analysis and prediction of genes associated with monogenic disorders on human chromosome X.

Disease gene discovery on chromosome (chr) X is challenging owing to its unique modes of inheritance. We undertook a systematic analysis of human chrX genes. We observe a higher proportion of disorder-associated genes and an enrichment of genes involved in cognition, language, and seizures on chrX compared to autosomes. We analyze gene constraints, exon and promoter conservation, expression, and paralogues, and report 127 genes sharing one or more attributes with known chrX disorder genes. Using machine learning classifiers trained to distinguish disease-associated from dispensable genes, we classify 247 genes, including 115 of the 127, as having high probability of being disease-associated. We provide evidence of an excess of variants in predicted genes in existing databases. Finally, we report damaging variants in CDK16 and TRPC5 in patients with intellectual disability or autism spectrum disorders. This study predicts large-scale gene-disease associations that could be used for prioritization of X-linked pathogenic variants.

Integration of Hi-C with short and long-read genome sequencing reveals the structure of germline rearranged genomes.

Structural variants are a common cause of disease and contribute to a large extent to inter-individual variability, but their detection and interpretation remain a challenge. Here, we investigate 11 individuals with complex genomic rearrangements including germline chromothripsis by combining short- and long-read genome sequencing (GS) with Hi-C. Large-scale genomic rearrangements are identified in Hi-C interaction maps, allowing for an independent assessment of breakpoint calls derived from the GS methods, resulting in >300 genomic junctions. Based on a comprehensive breakpoint detection and Hi-C, we achieve a reconstruction of whole rearranged chromosomes. Integrating information on the three-dimensional organization of chromatin, we observe that breakpoints occur more frequently than expected in lamina-associated domains (LADs) and that a majority reshuffle topologically associating domains (TADs). By applying phased RNA-seq, we observe an enrichment of genes showing allelic imbalanced expression (AIG) within 100 kb around the breakpoints. Interestingly, the AIGs hit by a breakpoint (19/22) display both up- and downregulation, thereby suggesting different mechanisms at play, such as gene disruption and rearrangements of regulatory information. However, the majority of interpretable genes located 200 kb around a breakpoint do not show significant expression changes. Thus, there is an overall robustness in the genome towards large-scale chromosome rearrangements.

Comprehensive genotype-phenotype correlation in AP-4 deficiency syndrome; Adding data from a large cohort of Iranian patients.

Mutations in adaptor protein complex-4 (AP-4) genes have first been identified in 2009, causing a phenotype termed as AP-4 deficiency syndrome. Since then several patients with overlapping phenotypes, comprised of intellectual disability (ID) and spastic tetraplegia have been reported. To delineate the genotype-phenotype correlation of the AP-4 deficiency syndrome, we add the data from 30 affected individuals from 12 out of 640 Iranian families with ID in whom we detected disease-causing variants in AP-4 complex subunits, using next-generation sequencing. Furthermore, by comparing genotype-phenotype findings of those affected individuals with previously reported patients, we further refine the genotype-phenotype correlation in this syndrome. The most frequent reported clinical findings in the 101 cases consist of ID and/or global developmental delay (97%), speech disorders (92.1%), inability to walk (90.1%), spasticity (77.2%), and microcephaly (75.2%). Spastic tetraplegia has been reported in 72.3% of the investigated patients. The major brain imaging findings are abnormal corpus callosum morphology (63.4%) followed by ventriculomegaly (44.5%). Our result might suggest the AP-4 deficiency syndrome as a major differential diagnostic for unknown hereditary neurodegenerative disorders.

Novel pathogenic EIF2S3 missense variants causing clinically variable MEHMO syndrome with impaired eIF2γ translational function, and literature review.

Rare pathogenic EIF2S3 missense and terminal deletion variants cause the X-linked intellectual disability (ID) syndrome MEHMO, or a milder phenotype including pancreatic dysfunction and hypopituitarism. We present two unrelated male patients who carry novel EIF2S3 pathogenic missense variants (p.(Thr144Ile) and p.(Ile159Leu)) thereby broadening the limited genetic spectrum and underscoring clinically variable expressivity of MEHMO. While the affected male with p.(Thr144Ile) presented with severe motor delay, severe microcephaly, moderate ID, epileptic seizures responsive to treatments, hypogenitalism, central obesity, facial features, and diabetes, the affected male with p.(Ile159Leu) presented with moderate ID, mild motor delay, microcephaly, epileptic seizures resistant to treatment, central obesity, and mild facial features. Both variants are located in the highly conserved guanine nucleotide binding domain of the EIF2S3 encoded eIF2γ subunit of the heterotrimeric translation initiation factor 2 (eIF2) complex. Further, we investigated both variants in a structural model and in yeast. The reduced growth rates and lowered fidelity of translation with increased initiation at non-AUG codons observed for both mutants in these studies strongly support pathogenicity of the variants.

A missense mutation in the CSTF2 gene that impairs the function of the RNA recognition motif and causes defects in 3' end processing is associated with intellectual disability in humans.

CSTF2 encodes an RNA-binding protein that is essential for mRNA cleavage and polyadenylation (C/P). No disease-associated mutations have been described for this gene. Here, we report a mutation in the RNA recognition motif (RRM) of CSTF2 that changes an aspartic acid at position 50 to alanine (p.D50A), resulting in intellectual disability in male patients. In mice, this mutation was sufficient to alter polyadenylation sites in over 1300 genes critical for brain development. Using a reporter gene assay, we demonstrated that C/P efficiency of CSTF2D50A was lower than wild type. To account for this, we determined that p.D50A changed locations of amino acid side chains altering RNA binding sites in the RRM. The changes modified the electrostatic potential of the RRM leading to a greater affinity for RNA. These results highlight the significance of 3' end mRNA processing in expression of genes important for brain plasticity and neuronal development.

Whole genome sequencing identifies a duplicated region encompassing Xq13.2q13.3 in a large Iranian family with intellectual disability.

BACKGROUND: The X chromosome has historically been one of the most thoroughly investigated chromosomes regarding intellectual disability (ID), whose etiology is attributed to many factors including copy number variations (CNVs). Duplications of the long arm of the X chromosome have been reported in patients with ID, short stature, facial anomalies, and in many cases hypoplastic genitalia and/or behavioral abnormalities. METHODS: Here, we report on a large Iranian family with X-linked ID caused by a duplication on the X chromosome identified by whole genome sequencing in combination with linkage analysis. RESULTS: Seven affected males in different branches of the family presented with ID, short stature, seizures, facial anomalies, behavioral abnormalities (aggressiveness, self-injury, anxiety, impaired social interactions, and shyness), speech impairment, and micropenis. The duplication of the region Xq13.2q13.3, which is ~1.8 Mb in size, includes seven protein-coding OMIM genes. Three of these genes, namely SLC16A2, RLIM, and NEXMIF, if impaired, can lead to syndromes presenting with ID. Of note, this duplicated region was located within a linkage interval with a LOD score >3. CONCLUSION: Our report indicates that CNVs should be considered in multi-affected families where no candidate gene defect has been identified in sequencing data analysis.

Hi-C Identifies Complex Genomic Rearrangements and TAD-Shuffling in Developmental Diseases.

Genome-wide analysis methods, such as array comparative genomic hybridization (CGH) and whole-genome sequencing (WGS), have greatly advanced the identification of structural variants (SVs) in the human genome. However, even with standard high-throughput sequencing techniques, complex rearrangements with multiple breakpoints are often difficult to resolve, and predicting their effects on gene expression and phenotype remains a challenge. Here, we address these problems by using high-throughput chromosome conformation capture (Hi-C) generated from cultured cells of nine individuals with developmental disorders (DDs). Three individuals had previously been identified as harboring duplications at the SOX9 locus and six had been identified with translocations. Hi-C resolved the positions of the duplications and was instructive in interpreting their distinct pathogenic effects, including the formation of new topologically associating domains (neo-TADs). Hi-C was very sensitive in detecting translocations, and it revealed previously unrecognized complex rearrangements at the breakpoints. In several cases, we observed the formation of fused-TADs promoting ectopic enhancer-promoter interactions that were likely to be involved in the disease pathology. In summary, we show that Hi-C is a sensible method for the detection of complex SVs in a clinical setting. The results help interpret the possible pathogenic effects of the SVs in individuals with DDs.

TAF1, associated with intellectual disability in humans, is essential for embryogenesis and regulates neurodevelopmental processes in zebrafish.

The TATA-box binding protein associated factor 1 (TAF1) protein is a key unit of the transcription factor II D complex that serves a vital function during transcription initiation. Variants of TAF1 have been associated with neurodevelopmental disorders, but TAF1's molecular functions remain elusive. In this study, we present a five-generation family affected with X-linked intellectual disability that co-segregated with a TAF1 c.3568C>T, p.(Arg1190Cys) variant. All affected males presented with intellectual disability and dysmorphic features, while heterozygous females were asymptomatic and had completely skewed X-chromosome inactivation. We investigated the role of TAF1 and its association to neurodevelopment by creating the first complete knockout model of the TAF1 orthologue in zebrafish. A crucial function of human TAF1 during embryogenesis can be inferred from the model, demonstrating that intact taf1 is essential for embryonic development. Transcriptome analysis of taf1 zebrafish knockout revealed enrichment for genes associated with neurodevelopmental processes. In conclusion, we propose that functional TAF1 is essential for embryonic development and specifically neurodevelopmental processes.

Deleterious de novo variants of X-linked ZC4H2 in females cause a variable phenotype with neurogenic arthrogryposis multiplex congenita.

Pathogenic variants in the X-linked gene ZC4H2, which encodes a zinc-finger protein, cause an infrequently described syndromic form of arthrogryposis multiplex congenita (AMC) with central and peripheral nervous system involvement. We present genetic and detailed phenotypic information on 23 newly identified families and simplex cases that include 19 affected females from 18 families and 14 affected males from nine families. Of note, the 15 females with deleterious de novo ZC4H2 variants presented with phenotypes ranging from mild to severe, and their clinical features overlapped with those seen in affected males. By contrast, of the nine carrier females with inherited ZC4H2 missense variants that were deleterious in affected male relatives, four were symptomatic. We also compared clinical phenotypes with previously published cases of both sexes and provide an overview on 48 males and 57 females from 42 families. The spectrum of ZC4H2 defects comprises novel and recurrent mostly inherited missense variants in affected males, and de novo splicing, frameshift, nonsense, and partial ZC4H2 deletions in affected females. Pathogenicity of two newly identified missense variants was further supported by studies in zebrafish. We propose ZC4H2 as a good candidate for early genetic testing of males and females with a clinical suspicion of fetal hypo-/akinesia and/or (neurogenic) AMC.

Identification of disease-causing variants in the EXOSC gene family underlying autosomal recessive intellectual disability in Iranian families.

Neurodevelopmental delay and intellectual disability (ID) can arise from numerous genetic defects. To date, variants in the EXOSC gene family have been associated with such disorders. Using next-generation sequencing (NGS), known and novel variants in this gene family causing autosomal recessive ID (ARID) have been identified in five Iranian families. By collecting clinical information on these families and comparing their phenotypes with previously reported patients, we further describe the clinical variability of ARID resulting from alterations in the EXOSC gene family, and emphasize the role of RNA processing dysregulation in ARID.

Multigenic truncation of the semaphorin-plexin pathway by a germline chromothriptic rearrangement associated with Moebius syndrome.

Moebius syndrome (MBS) is a congenital disorder caused by paralysis of the facial and abducens nerves. Although a number of candidate genes have been suspected, so far only mutations in PLXND1 and REV3L are confirmed to cause MBS. Here, we fine mapped the breakpoints of a complex chromosomal rearrangement (CCR) 46,XY,t(7;8;11;13) in a patient with MBS, which revealed 41 clustered breakpoints with typical hallmarks of chromothripsis. Among 12 truncated protein-coding genes, SEMA3A is known to bind to the MBS-associated PLXND1. Intriguingly, the CCR also truncated PIK3CG, which in silico interacts with REVL3 encoded by the other known MBS-gene REV3L, and with the SEMA3A/PLXND1 complex via FLT1. Additional studies of other complex rearrangements may reveal whether the multiple breakpoints in germline chromothripsis may predispose to complex multigenic disorders.

Mutation p.R356Q in the Collybistin Phosphoinositide Binding Site Is Associated With Mild Intellectual Disability.

The recruitment of inhibitory GABAA receptors to neuronal synapses requires a complex interplay between receptors, neuroligins, the scaffolding protein gephyrin and the GDP-GTP exchange factor collybistin (CB). Collybistin is regulated by protein-protein interactions at the N-terminal SH3 domain, which can bind neuroligins 2/4 and the GABAAR α2 subunit. Collybistin also harbors a RhoGEF domain which mediates interactions with gephyrin and catalyzes GDP-GTP exchange on Cdc42. Lastly, collybistin has a pleckstrin homology (PH) domain, which binds phosphoinositides, such as phosphatidylinositol 3-phosphate (PI3P/PtdIns3P) and phosphatidylinositol 4-monophosphate (PI4P/PtdIns4P). PI3P located in early/sorting endosomes has recently been shown to regulate the postsynaptic clustering of gephyrin and GABAA receptors and consequently the strength of inhibitory synapses in cultured hippocampal neurons. This process is disrupted by mutations in the collybistin gene (ARHGEF9), which cause X-linked intellectual disability (XLID) by a variety of mechanisms converging on disrupted gephyrin and GABAA receptor clustering at central synapses. Here we report a novel missense mutation (chrX:62875607C>T, p.R356Q) in ARHGEF9 that affects one of the two paired arginine residues in the PH domain that were predicted to be vital for binding phosphoinositides. Functional assays revealed that recombinant collybistin CB3SH3- R356Q was deficient in PI3P binding and was not able to translocate EGFP-gephyrin to submembrane microaggregates in an in vitro clustering assay. Expression of the PI3P-binding mutants CB3SH3- R356Q and CB3SH3- R356N/R357N in cultured hippocampal neurones revealed that the mutant proteins did not accumulate at inhibitory synapses, but instead resulted in a clear decrease in the overall number of synaptic gephyrin clusters compared to controls. Molecular dynamics simulations suggest that the p.R356Q substitution influences PI3P binding by altering the range of structural conformations adopted by collybistin. Taken together, these results suggest that the p.R356Q mutation in ARHGEF9 is the underlying cause of XLID in the probands, disrupting gephyrin clustering at inhibitory GABAergic synapses via loss of collybistin PH domain phosphoinositide binding.

Pathogenic variants in E3 ubiquitin ligase RLIM/RNF12 lead to a syndromic X-linked intellectual disability and behavior disorder.

RLIM, also known as RNF12, is an X-linked E3 ubiquitin ligase acting as a negative regulator of LIM-domain containing transcription factors and participates in X-chromosome inactivation (XCI) in mice. We report the genetic and clinical findings of 84 individuals from nine unrelated families, eight of whom who have pathogenic variants in RLIM (RING finger LIM domain-interacting protein). A total of 40 affected males have X-linked intellectual disability (XLID) and variable behavioral anomalies with or without congenital malformations. In contrast, 44 heterozygous female carriers have normal cognition and behavior, but eight showed mild physical features. All RLIM variants identified are missense changes co-segregating with the phenotype and predicted to affect protein function. Eight of the nine altered amino acids are conserved and lie either within a domain essential for binding interacting proteins or in the C-terminal RING finger catalytic domain. In vitro experiments revealed that these amino acid changes in the RLIM RING finger impaired RLIM ubiquitin ligase activity. In vivo experiments in rlim mutant zebrafish showed that wild type RLIM rescued the zebrafish rlim phenotype, whereas the patient-specific missense RLIM variants failed to rescue the phenotype and thus represent likely severe loss-of-function mutations. In summary, we identified a spectrum of RLIM missense variants causing syndromic XLID and affecting the ubiquitin ligase activity of RLIM, suggesting that enzymatic activity of RLIM is required for normal development, cognition and behavior.

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.

A recurrent missense variant in SLC9A7 causes nonsyndromic X-linked intellectual disability with alteration of Golgi acidification and aberrant glycosylation.

We report two unrelated families with multigenerational nonsyndromic intellectual disability (ID) segregating with a recurrent de novo missense variant (c.1543C>T:p.Leu515Phe) in the alkali cation/proton exchanger gene SLC9A7 (also commonly referred to as NHE7). SLC9A7 is located on human X chromosome at Xp11.3 and has not yet been associated with a human phenotype. The gene is widely transcribed, but especially abundant in brain, skeletal muscle and various secretory tissues. Within cells, SLC9A7 resides in the Golgi apparatus, with prominent enrichment in the trans-Golgi network (TGN) and post-Golgi vesicles. In transfected Chinese hamster ovary AP-1 cells, the Leu515Phe mutant protein was correctly targeted to the TGN/post-Golgi vesicles, but its N-linked oligosaccharide maturation as well as that of a co-transfected secretory membrane glycoprotein, vesicular stomatitis virus G (VSVG) glycoprotein, was reduced compared to cells co-expressing SLC9A7 wild-type and VSVG. This correlated with alkalinization of the TGN/post-Golgi compartments, suggestive of a gain-of-function. Membrane trafficking of glycosylation-deficient Leu515Phe and co-transfected VSVG to the cell surface, however, was relatively unaffected. Mass spectrometry analysis of patient sera also revealed an abnormal N-glycosylation profile for transferrin, a clinical diagnostic marker for congenital disorders of glycosylation. These data implicate a crucial role for SLC9A7 in the regulation of TGN/post-Golgi pH homeostasis and glycosylation of exported cargo, which may underlie the cellular pathophysiology and neurodevelopmental deficits associated with this particular nonsyndromic form of X-linked ID.

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.