Delineating the Spectrum of Genetic Variants Associated with Bardet-Biedl Syndrome in Consanguineous Pakistani Pedigrees.

This study aimed to find the molecular basis of Bardet-Biedl syndrome (BBS) in Pakistani consanguineous families. A total of 12 affected families were enrolled. Clinical investigations were performed to access the BBS-associated phenotypes. Whole exome sequencing was conducted on one affected individual from each family. The computational functional analysis predicted the variants' pathogenic effects and modeled the mutated proteins. Whole-exome sequencing revealed 9 pathogenic variants in six genes associated with BBS in 12 families. The BBS6/MKS was the most common BBS causative gene identified in five families (5/12, 41.6%), with one novel (c.1226G>A, p.Gly409Glu) and two reported variants. c.774G>A, Thr259LeuTer21 was the most frequent BBS6/MMKS allele in three families 3/5 (60%). Two variants, c.223C>T, p.Arg75Ter and a novel, c. 252delA, p.Lys85STer39 were detected in the BBS9 gene. A novel 8bp deletion c.387_394delAAATAAAA, p. Asn130GlyfsTer3 was found in BBS3 gene. Three known variants were detected in the BBS1, BBS2, and BBS7 genes. Identification of novel likely pathogenic variants in three genes reaffirms the allelic and genetic heterogeneity of BBS in Pakistani patients. The clinical differences among patients carrying the same pathogenic variant may be due to other factors influencing the phenotype, including variants in other modifier genes.

Limited evidence for blood eQTLs in human sexual dimorphism.

BACKGROUND: The genetic underpinning of sexual dimorphism is very poorly understood. The prevalence of many diseases differs between men and women, which could be in part caused by sex-specific genetic effects. Nevertheless, only a few published genome-wide association studies (GWAS) were performed separately in each sex. The reported enrichment of expression quantitative trait loci (eQTLs) among GWAS-associated SNPs suggests a potential role of sex-specific eQTLs in the sex-specific genetic mechanism underlying complex traits. METHODS: To explore this scenario, we combined sex-specific whole blood RNA-seq eQTL data from 3447 European individuals included in BIOS Consortium and GWAS data from UK Biobank. Next, to test the presence of sex-biased causal effect of gene expression on complex traits, we performed sex-specific transcriptome-wide Mendelian randomization (TWMR) analyses on the two most sexually dimorphic traits, waist-to-hip ratio (WHR) and testosterone levels. Finally, we performed power analysis to calculate the GWAS sample size needed to observe sex-specific trait associations driven by sex-biased eQTLs. RESULTS: Among 9 million SNP-gene pairs showing sex-combined associations, we found 18 genes with significant sex-biased cis-eQTLs (FDR 5%). Our phenome-wide association study of the 18 top sex-biased eQTLs on >700 traits unraveled that these eQTLs do not systematically translate into detectable sex-biased trait-associations. In addition, we observed that sex-specific causal effects of gene expression on complex traits are not driven by sex-specific eQTLs. Power analyses using real eQTL- and causal-effect sizes showed that millions of samples would be necessary to observe sex-biased trait associations that are fully driven by sex-biased cis-eQTLs. Compensatory effects may further hamper their detection. CONCLUSIONS: Our results suggest that sex-specific eQTLs in whole blood do not translate to detectable sex-specific trait associations of complex diseases, and vice versa that the observed sex-specific trait associations cannot be explained by sex-specific eQTLs.

CoverageMaster: comprehensive CNV detection and visualization from NGS short reads for genetic medicine applications.

CoverageMaster (CoM) is a copy number variation (CNV) calling algorithm based on depth-of-coverage maps designed to detect CNVs of any size in exome [whole exome sequencing (WES)] and genome [whole genome sequencing (WGS)] data. The core of the algorithm is the compression of sequencing coverage data in a multiscale Wavelet space and the analysis through an iterative Hidden Markov Model. CoM processes WES and WGS data at nucleotide scale resolution and accurately detects and visualizes full size range CNVs, including single or partial exon deletions and duplications. The results obtained with this approach support the possibility for coverage-based CNV callers to replace probe-based methods such as array comparative genomic hybridization and multiplex ligation-dependent probe amplification in the near future.

Differentially expressed genes reflect disease-induced rather than disease-causing changes in the transcriptome.

Comparing transcript levels between healthy and diseased individuals allows the identification of differentially expressed genes, which may be causes, consequences or mere correlates of the disease under scrutiny. We propose a method to decompose the observational correlation between gene expression and phenotypes driven by confounders, forward- and reverse causal effects. The bi-directional causal effects between gene expression and complex traits are obtained by Mendelian Randomization integrating summary-level data from GWAS and whole-blood eQTLs. Applying this approach to complex traits reveals that forward effects have negligible contribution. For example, BMI- and triglycerides-gene expression correlation coefficients robustly correlate with trait-to-expression causal effects (rBMI = 0.11, PBMI = 2.0 × 10-51 and rTG = 0.13, PTG = 1.1 × 10-68), but not detectably with expression-to-trait effects. Our results demonstrate that studies comparing the transcriptome of diseased and healthy subjects are more prone to reveal disease-induced gene expression changes rather than disease causing ones.

Biallelic variants in TMEM222 cause a new autosomal recessive neurodevelopmental disorder.

PURPOSE: To elucidate the novel molecular cause in families with a new autosomal recessive neurodevelopmental disorder. METHODS: A combination of exome sequencing and gene matching tools was used to identify pathogenic variants in 17 individuals. Quantitative reverse transcription polymerase chain reaction (RT-qPCR) and subcellular localization studies were used to characterize gene expression profile and localization. RESULTS: Biallelic variants in the TMEM222 gene were identified in 17 individuals from nine unrelated families, presenting with intellectual disability and variable other features, such as aggressive behavior, shy character, body tremors, decreased muscle mass in the lower extremities, and mild hypotonia. We found relatively high TMEM222 expression levels in the human brain, especially in the parietal and occipital cortex. Additionally, subcellular localization analysis in human neurons derived from induced pluripotent stem cells (iPSCs) revealed that TMEM222 localizes to early endosomes in the synapses of mature iPSC-derived neurons. CONCLUSION: Our findings support a role for TMEM222 in brain development and function and adds variants in the gene TMEM222 as a novel underlying cause of an autosomal recessive neurodevelopmental disorder.

Dominant monoallelic variant in the PAK2 gene causes Knobloch syndrome type 2.

Knobloch syndrome is an autosomal recessive phenotype mainly characterized by retinal detachment and encephalocele caused by biallelic pathogenic variants in the COL18A1 gene. However, there are patients clinically diagnosed as Knobloch syndrome with unknown molecular etiology not linked to COL18A1. We studied an historical pedigree (published in 1998) designated as KNO2 (Knobloch type 2 syndrome with intellectual disability, autistic behavior, retinal degeneration, encephalocele). Whole exome sequencing of the two affected siblings and the normal parents resulted in the identification of a PAK2 non-synonymous substitution p.(Glu435Lys) as a causative variant. The variant was monoallelic and apparently de novo in both siblings indicating a likely germ-line mosaicism in one of the parents; the mosaicism, however, could not be observed after deep sequencing of blood parental DNA. PAK2 encodes a member of a small group of serine/threonine kinases; these P21-activating kinases (PAKs) are essential in signal transduction and cellular regulation (cytoskeletal dynamics, cell motility, death and survival signaling and cell cycle progression). Structural analysis of the PAK2 p.(Glu435Lys) variant that is located in the kinase domain of the protein predicts a possible compromise in the kinase activity. Functional analysis of the p.(Glu435Lys) PAK2 variant in transfected HEK293T cells results in a partial loss of the kinase activity. PAK2 has been previously suggested as an autism-related gene. Our results show that PAK2-induced phenotypic spectrum is broad and not fully understood. We conclude that the KNO2 syndrome in the studied family is dominant and caused by a deleterious variant in the PAK2 gene.

Genetic variation near CXCL12 is associated with susceptibility to HIV-related non-Hodgkin lymphoma.

Human immunodeficiency virus (HIV) infection is associated with an increased risk of non-Hodgkin lymphoma (NHL). Even in the era of suppressive antiretroviral treatment, HIV-infected individuals remain at higher risk of developing NHL compared to the general population. To identify potential genetic risk loci, we performed case-control genome-wide association studies and a meta-analysis across three cohorts of HIV+ patients of European ancestry, including a total of 278 cases and 1924 matched controls. We observed a significant association with NHL susceptibility in the C-X-C motif chemokine ligand 12 (CXCL12) region on chromosome 10. A fine mapping analysis identified rs7919208 as the most likely causal variant (P = 4.77e-11), with the G>A polymorphism creating a new transcription factor binding site for BATF and JUND. These results suggest a modulatory role of CXCL12 regulation in the increased susceptibility to NHL observed in the HIV-infected population.

De Novo and Bi-allelic Pathogenic Variants in NARS1 Cause Neurodevelopmental Delay Due to Toxic Gain-of-Function and Partial Loss-of-Function Effects.

Aminoacyl-tRNA synthetases (ARSs) are ubiquitous, ancient enzymes that charge amino acids to cognate tRNA molecules, the essential first step of protein translation. Here, we describe 32 individuals from 21 families, presenting with microcephaly, neurodevelopmental delay, seizures, peripheral neuropathy, and ataxia, with de novo heterozygous and bi-allelic mutations in asparaginyl-tRNA synthetase (NARS1). We demonstrate a reduction in NARS1 mRNA expression as well as in NARS1 enzyme levels and activity in both individual fibroblasts and induced neural progenitor cells (iNPCs). Molecular modeling of the recessive c.1633C>T (p.Arg545Cys) variant shows weaker spatial positioning and tRNA selectivity. We conclude that de novo and bi-allelic mutations in NARS1 are a significant cause of neurodevelopmental disease, where the mechanism for de novo variants could be toxic gain-of-function and for recessive variants, partial loss-of-function.

Pathogenic mosaic variants in congenital hypogonadotropic hypogonadism.

PURPOSE: Congenital hypogonadotropic hypogonadism (CHH) is a rare disorder resulting in absent puberty and infertility. The genetic architecture is complex with multiple loci involved, variable expressivity, and incomplete penetrance. The majority of cases are sporadic, consistent with a disease affecting fertility. The current study aims to investigate mosaicism as a genetic mechanism for CHH, focusing on de novo rare variants in CHH genes. METHODS: We evaluated 60 trios for de novo rare sequencing variants (RSV) in known CHH genes using exome sequencing. Potential mosaicism was suspected among RSVs with altered allelic ratios and confirmed using customized ultradeep sequencing (UDS) in multiple tissues. RESULTS: Among the 60 trios, 10 probands harbored de novo pathogenic variants in CHH genes. Custom UDS demonstrated that three of these de novo variants were in fact postzygotic mosaicism-two in FGFR1 (p.Leu630Pro and p.Gly348Arg), and one in CHD7 (p.Arg2428*). Statistically significant variation across multiple tissues (DNA from blood, buccal, hair follicle, urine) confirmed their mosaic nature. CONCLUSIONS: We identified a significant number of de novo pathogenic variants in CHH of which a notable number (3/10) exhibited mosaicism. This report of postzygotic mosaicism in CHH patients provides valuable information for accurate genetic counseling.

Biallelic variants in PSMB1 encoding the proteasome subunit ß6 cause impairment of proteasome function, microcephaly, intellectual disability, developmental delay and short stature.

The molecular cause of the majority of rare autosomal recessive disorders remains unknown. Consanguinity due to extensive homozygosity unravels many recessive phenotypes and facilitates the detection of novel gene-disease links. Here, we report two siblings with phenotypic signs, including intellectual disability (ID), developmental delay and microcephaly from a Pakistani consanguineous family in which we have identified homozygosity for p(Tyr103His) in the PSMB1 gene (Genbank NM_002793) that segregated with the disease phenotype. PSMB1 encodes a ß-type proteasome subunit (i.e. ß6). Modeling of the p(Tyr103His) variant indicates that this variant weakens the interactions between PSMB1/ß6 and PSMA5/α5 proteasome subunits and thus destabilizes the 20S proteasome complex. Biochemical experiments in human SHSY5Y cells revealed that the p(Tyr103His) variant affects both the processing of PSMB1/ß6 and its incorporation into proteasome, thus impairing proteasome activity. CRISPR/Cas9 mutagenesis or morpholino knock-down of the single psmb1 zebrafish orthologue resulted in microcephaly, microphthalmia and reduced brain size. Genetic evidence in the family and functional experiments in human cells and zebrafish indicates that PSMB1/ß6 pathogenic variants are the cause of a recessive disease with ID, microcephaly and developmental delay due to abnormal proteasome assembly.

Taurine treatment of retinal degeneration and cardiomyopathy in a consanguineous family with SLC6A6 taurine transporter deficiency.

In a consanguineous Pakistani family with two affected individuals, a homozygous variant Gly399Val in the eighth transmembrane domain of the taurine transporter SLC6A6 was identified resulting in a hypomorph transporting capacity of ~15% compared with normal. Three-dimensional modeling of this variant has indicated that it likely causes displacement of the Tyr138 (TM3) side chain, important for transport of taurine. The affected individuals presented with rapidly progressive childhood retinal degeneration, cardiomyopathy and almost undetectable plasma taurine levels. Oral taurine supplementation of 100 mg/kg/day resulted in maintenance of normal blood taurine levels. Following approval by the ethics committee, a long-term supplementation treatment was introduced. Remarkably, after 24-months, the cardiomyopathy was corrected in both affected siblings, and in the 6-years-old, the retinal degeneration was arrested, and the vision was clinically improved. Similar therapeutic approaches could be employed in Mendelian phenotypes caused by the dysfunction of the hundreds of other molecular transporters.

Bi-allelic Variants in IQSEC1 Cause Intellectual Disability, Developmental Delay, and Short Stature.

We report two consanguineous families with probands that exhibit intellectual disability, developmental delay, short stature, aphasia, and hypotonia in which homozygous non-synonymous variants were identified in IQSEC1 (GenBank: NM_001134382.3). In a Pakistani family, the IQSEC1 segregating variant is c.1028C>T (p.Thr343Met), while in a Saudi Arabian family the variant is c.962G>A (p.Arg321Gln). IQSEC1-3 encode guanine nucleotide exchange factors for the small GTPase ARF6 and their loss affects a variety of actin-dependent cellular processes, including AMPA receptor trafficking at synapses. The ortholog of IQSECs in the fly is schizo and its loss affects growth cone guidance at the midline in the CNS, also an actin-dependent process. Overexpression of the reference IQSEC1 cDNA in wild-type flies is lethal, but overexpression of the two variant IQSEC1 cDNAs did not affect viability. Loss of schizo caused embryonic lethality that could be rescued to 2nd instar larvae by moderate expression of the human reference cDNA. However, the p.Arg321Gln and p.Thr343Met variants failed to rescue embryonic lethality. These data indicate that the variants behave as loss-of-function mutations. We also show that schizo in photoreceptors is required for phototransduction. Finally, mice with a conditional Iqsec1 deletion in cortical neurons exhibited an increased density of dendritic spines with an immature morphology. The phenotypic similarity of the affecteds and the functional experiments in flies and mice indicate that IQSEC1 variants are the cause of a recessive disease with intellectual disability, developmental delay, and short stature, and that axonal guidance and dendritic projection defects as well as dendritic spine dysgenesis may underlie disease pathogenesis.

Mendelian randomization integrating GWAS and eQTL data reveals genetic determinants of complex and clinical traits.

Genome-wide association studies (GWAS) have identified thousands of variants associated with complex traits, but their biological interpretation often remains unclear. Most of these variants overlap with expression QTLs, indicating their potential involvement in regulation of gene expression. Here, we propose a transcriptome-wide summary statistics-based Mendelian Randomization approach (TWMR) that uses multiple SNPs as instruments and multiple gene expression traits as exposures, simultaneously. Applied to 43 human phenotypes, it uncovers 3,913 putatively causal gene-trait associations, 36% of which have no genome-wide significant SNP nearby in previous GWAS. Using independent association summary statistics, we find that the majority of these loci were missed by GWAS due to power issues. Noteworthy among these links is educational attainment-associated BSCL2, known to carry mutations leading to a Mendelian form of encephalopathy. We also find pleiotropic causal effects suggestive of mechanistic connections. TWMR better accounts for pleiotropy and has the potential to identify biological mechanisms underlying complex traits.

Bi-allelic Variants in DYNC1I2 Cause Syndromic Microcephaly with Intellectual Disability, Cerebral Malformations, and Dysmorphic Facial Features.

Cargo transport along the cytoplasmic microtubular network is essential for neuronal function, and cytoplasmic dynein-1 is an established molecular motor that is critical for neurogenesis and homeostasis. We performed whole-exome sequencing, homozygosity mapping, and chromosomal microarray studies in five individuals from three independent pedigrees and identified likely-pathogenic variants in DYNC1I2 (Dynein Cytoplasmic 1 Intermediate Chain 2), encoding a component of the cytoplasmic dynein 1 complex. In a consanguineous Pakistani family with three affected individuals presenting with microcephaly, severe intellectual disability, simplification of cerebral gyration, corpus callosum hypoplasia, and dysmorphic facial features, we identified a homozygous splice donor site variant (GenBank: NM_001378.2:c.607+1G>A). We report two additional individuals who have similar neurodevelopmental deficits and craniofacial features and harbor deleterious variants; one individual bears a c.740A>G (p.Tyr247Cys) change in trans with a 374 kb deletion encompassing DYNC1I2, and an unrelated individual harbors the compound-heterozygous variants c.868C>T (p.Gln290∗) and c.740A>G (p.Tyr247Cys). Zebrafish larvae subjected to CRISPR-Cas9 gene disruption or transient suppression of dync1i2a displayed significantly altered craniofacial patterning with concomitant reduction in head size. We monitored cell death and cell cycle progression in dync1i2a zebrafish models and observed significantly increased apoptosis, likely due to prolonged mitosis caused by abnormal spindle morphology, and this finding offers initial insights into the cellular basis of microcephaly. Additionally, complementation studies in zebrafish demonstrate that p.Tyr247Cys attenuates gene function, consistent with protein structural analysis. Our genetic and functional data indicate that DYNC1I2 dysfunction probably causes an autosomal-recessive microcephaly syndrome and highlight further the critical roles of the dynein-1 complex in neurodevelopment.

Biallelic variants in FBXL3 cause intellectual disability, delayed motor development and short stature.

FBXL3 (F-Box and Leucine Rich Repeat Protein 3) encodes a protein that contains an F-box and several tandem leucine-rich repeats (LRR) domains. FBXL3 is part of the SCF (Skp1-Cullin-F box protein) ubiquitin ligase complex that binds and leads to phosphorylation-dependent degradation of the central clock protein cryptochromes (CRY1 and CRY2) by the proteasome and its absence causes circadian phenotypes in mice and behavioral problems. No FBXL3-related phenotypes have been described in humans. By a combination of exome sequencing and homozygosity mapping, we analyzed two consanguineous families with intellectual disability and identified homozygous loss-of-function (LoF) variants in FBXL3. In the first family, from Pakistan, an FBXL3 frameshift variant [NM_012158.2:c.885delT:p.(Leu295Phefs*25)] was the onlysegregating variant in five affected individuals in two family loops (LOD score: 3.12). In the second family, from Lebanon, we identified a nonsense variant [NM_012158.2:c.445C>T:p.(Arg149*)]. In a third patient from Italy, a likely deleterious non-synonymous variant [NM_012158.2:c.1072T>C:p.(Cys358Arg)] was identified in homozygosity. Protein 3D modeling predicted that the Cys358Arg change influences the binding with CRY2 by destabilizing the structure of the FBXL3, suggesting that this variant is also likely to be LoF. The eight affected individuals from the three families presented with a similar phenotype that included intellectual disability, developmental delay, short stature and mild facial dysmorphism, mainly large nose with a bulbous tip. The phenotypic similarity and the segregation analysis suggest that FBXL3 biallelic, LoF variants link this gene with syndromic autosomal recessive developmental delay/intellectual disability.

Biallelic loss of function variants in PPP1R21 cause a neurodevelopmental syndrome with impaired endocytic function.

Next-generation sequencing (NGS) has been instrumental in solving the genetic basis of rare inherited diseases, especially neurodevelopmental syndromes. However, functional workup is essential for precise phenotype definition and to understand the underlying disease mechanisms. Using whole exome (WES) and whole genome sequencing (WGS) in four independent families with hypotonia, neurodevelopmental delay, facial dysmorphism, loss of white matter, and thinning of the corpus callosum, we identified four previously unreported homozygous truncating PPP1R21 alleles: c.347delT p.(Ile116Lysfs*25), c.2170_2171insGGTA p.(Ile724Argfs*8), c.1607dupT p.(Leu536Phefs*7), c.2063delA p.(Lys688Serfs*26) and found that PPP1R21 was absent in fibroblasts of an affected individual, supporting the allele's loss of function effect. PPP1R21 function had not been studied except that a large scale affinity proteomics approach suggested an interaction with PIBF1 defective in Joubert syndrome. Our co-immunoprecipitation studies did not confirm this but in contrast defined the localization of PPP1R21 to the early endosome. Consistent with the subcellular expression pattern and the clinical phenotype exhibiting features of storage diseases, we found patient fibroblasts exhibited a delay in clearance of transferrin-488 while uptake was normal. In summary, we delineate a novel neurodevelopmental syndrome caused by biallelic PPP1R21 loss of function variants, and suggest a role of PPP1R21 within the endosomal sorting process or endosome maturation pathway.

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.

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.

Congenital hypogonadotropic hypogonadism and constitutional delay of growth and puberty have distinct genetic architectures.

OBJECTIVE: Congenital hypogonadotropic hypogonadism (CHH) and constitutional delay of growth and puberty (CDGP) represent rare and common forms of GnRH deficiency, respectively. Both CDGP and CHH present with delayed puberty, and the distinction between these two entities during early adolescence is challenging. More than 30 genes have been implicated in CHH, while the genetic basis of CDGP is poorly understood. DESIGN: We characterized and compared the genetic architectures of CHH and CDGP, to test the hypothesis of a shared genetic basis between these disorders. METHODS: Exome sequencing data were used to identify rare variants in known genes in CHH (n = 116), CDGP (n = 72) and control cohorts (n = 36 874 ExAC and n = 405 CoLaus). RESULTS: Mutations in at least one CHH gene were found in 51% of CHH probands, which is significantly higher than in CDGP (7%, P = 7.6 × 10-11) or controls (18%, P = 5.5 × 10-12). Similarly, oligogenicity (defined as mutations in more than one gene) was common in CHH patients (15%) relative to CDGP (1.4%, P = 0.002) and controls (2%, P = 6.4 × 10-7). CONCLUSIONS: Our data suggest that CDGP and CHH have distinct genetic profiles, and this finding may facilitate the differential diagnosis in patients presenting with delayed puberty.

Biallelic variants in KIF14 cause intellectual disability with microcephaly.

Kinesin proteins are critical for various cellular functions such as intracellular transport and cell division, and many members of the family have been linked to monogenic disorders and cancer. We report eight individuals with intellectual disability and microcephaly from four unrelated families with parental consanguinity. In the affected individuals of each family, homozygosity for likely pathogenic variants in KIF14 were detected; two loss-of-function (p.Asn83Ilefs*3 and p.Ser1478fs), and two missense substitutions (p.Ser841Phe and p.Gly459Arg). KIF14 is a mitotic motor protein that is required for spindle localization of the mitotic citron rho-interacting kinase, CIT, also mutated in microcephaly. Our results demonstrate the involvement of KIF14 in development and reveal a wide phenotypic variability ranging from fetal lethality to moderate developmental delay and microcephaly.