A de novo paradigm for male infertility

De novo mutations are known to play a prominent role in sporadic disorders with reduced fitness. We hypothesize that de novo mutations play an important role in severe male infertility and explain a portion of the genetic causes of this understudied disorder. To test this hypothesis, we utilize trio-based exome sequencing in a cohort of 185 infertile males and their unaffected parents. Following a systematic analysis, 29 of 145 rare (MAF < 0.1%) protein-altering de novo mutations are classified as possibly causative of the male infertility phenotype. We observed a significant enrichment of loss-of-function de novo mutations in loss-of-function-intolerant genes (p-value = 1.00 × 10^-5) in infertile men compared to controls. Additionally, we detected a significant increase in predicted pathogenic de novo missense mutations affecting missense-intolerant genes (p-value = 5.01 × 10^-4) in contrast to predicted benign de novo mutations. One gene we identify, RBM5, is an essential regulator of male germ cell pre-mRNA splicing and has been previously implicated in male infertility in mice. In a follow-up study, 6 rare pathogenic missense mutations affecting this gene are observed in a cohort of 2,506 infertile patients, whilst we find no such mutations in a cohort of 5,784 fertile men (p-value = 0.03). Our results provide evidence for the role of de novo mutations in severe male infertility and point to new candidate genes affecting fertility.

Genetic variation in genes of inborn errors of immunity in children with unexplained encephalitis

Pediatric encephalitis has significant morbidity and mortality, yet 50% of cases are unexplained. Host genetics plays a role in encephalitis' development; however, the contributing variants are poorly understood. One child with anti-NMDA receptor encephalitis and ten with unexplained encephalitis underwent whole genome sequencing to identify rare candidate variants in genes known to cause monogenic immunologic and neurologic disorders, and polymorphisms associated with increased disease risk. Using the professional Human Genetic Mutation Database (Qiagen), we divided the candidate variants into three categories: monogenic deleterious or potentially deleterious variants (1) in a disease-consistent inheritance pattern; (2) in carrier states; and (3) disease-related polymorphisms. Six patients (55%) had a deleterious or potentially deleterious variant in a disease-consistent inheritance pattern, five (45%) were heterozygous carriers for an autosomal recessive condition, and six (55%) carried a disease-related polymorphism. Finally, seven (64%) had more than one variant, suggesting possible polygenetic risk. Among variants identified were those implicated in atypical hemolytic uremic syndrome, common variable immunodeficiency, hemophagocytic lymphohistiocytosis, and systemic lupus erythematosus. This preliminary study shows genetic variation related to inborn errors of immunity in acute pediatric encephalitis. Future research is needed to determine if these variants play a functional role in the development of unexplained encephalitis.

The genetics of obesity: from discovery to biology

The prevalence of obesity has tripled over the past four decades, imposing an enormous burden on people's health. Polygenic (or common) obesity and rare, severe, early-onset monogenic obesity are often polarized as distinct diseases. However, gene discovery studies for both forms of obesity show that they have shared genetic and biological underpinnings, pointing to a key role for the brain in the control of body weight. Genome-wide association studies (GWAS) with increasing sample sizes and advances in sequencing technology are the main drivers behind a recent flurry of new discoveries. However, it is the post-GWAS, cross-disciplinary collaborations, which combine new omics technologies and analytical approaches, that have started to facilitate translation of genetic loci into meaningful biology and new avenues for treatment.

An Investigation of the Knowledge Overlap between Pharmacogenomics and Disease Genetics

Precision medicine faces many challenges, including the gap of knowledge between disease genetics and pharmacogenomics (PGx). Disease genetics interprets the pathogenicity of genetic variants for diagnostic purposes, while PGx investigates the genetic influences on drug responses. Ideally, the quality of health care would be improved from the point of disease diagnosis to drug prescribing if PGx is integrated with disease genetics in clinical care. However, PGx genes or variants are usually not reported as a secondary finding even if they are included in a clinical genetic test for diagnostic purposes. This happens even though the detection of PGx variants can provide valuable drug prescribing recommendations. One underlying reason is the lack of systematic classification of the knowledge overlap between PGx and disease genetics. Here, we address this issue by analyzing gene and genetic variant annotations from multiple expert-curated knowledge databases, including PharmGKB, CPIC, ClinGen and ClinVar. We further classified genes based on the strength of evidence supporting a gene's pathogenic role or PGx effect as well as the level of clinical actionability of a gene. Twenty-six genes were found to have pathogenic variation associated with germline diseases as well as strong evidence for a PGx association. These genes were classified into four sub-categories based on the distinct connection between the gene's pathogenic role and PGx effect. Moreover, we have also found thirteen RYR1 genetic variants that were annotated as pathogenic and at the same time whose PGx effect was supported by a preponderance of evidence and given drug prescribing recommendations. Overall, we identified a nontrivial number of gene and genetic variant overlaps between disease genetics and PGx, which laid out a foundation for combining PGx and disease genetics to improve clinical care from disease diagnoses to drug prescribing and adherence.

Stabilization of negative activation voltages of Cav1.3 L-Type Ca2+-channels by alternative splicing

-->Low voltage-activated Cav1.3 L-type Ca2+-channels are key regulators of neuronal excitability controlling neuronal development and different types of learning and memory. Their physiological functions are enabled by their negative activation voltage-range, which allows Cav1.3 to be active at subthreshold voltages. Alternative splicing in the C-terminus of their pore-forming α1-subunits gives rise to C-terminal long (Cav1.3L) and short (Cav1.3S) splice variants allowing Cav1.3S to activate at even more negative voltages than Cav1.3L. We discovered that inclusion of exons 8b, 11, and 32 in Cav1.3S further shifts activation (-3 to -4 mV) and inactivation (-4 to -6 mV) to more negative voltages as revealed by functional characterization in tsA-201 cells. We found transcripts of these exons in mouse chromaffin cells, the cochlea, and the brain. Our data further suggest that Cav1.3-containing exons 11 and 32 constitute a significant part of native channels in the brain. We therefore investigated the effect of these splice variants on human disease variants. Splicing did not prevent the gating defects of the previously reported human pathogenic variant S652L, which further shifted the voltage-dependence of activation of exon 11-containing channels by more than -12 mV. In contrast, we found no evidence for gating changes of the CACNA1D missense variant R498L, located in exon 11, which has recently been identified in a patient with an epileptic syndrome. Our data demonstrate that alternative splicing outside the C-terminus involving exons 11 and 32 contributes to channel fine-tuning by stabilizing negative activation and inactivation gating properties of wild-type and mutant Cav1.3 channels.

Incorporating functional priors improves polygenic prediction accuracy in UK Biobank and 23andMe data sets

Polygenic risk prediction is a widely investigated topic because of its promising clinical applications. Genetic variants in functional regions of the genome are enriched for complex trait heritability. Here, we introduce a method for polygenic prediction, LDpred-funct, that leverages trait-specific functional priors to increase prediction accuracy. We fit priors using the recently developed baseline-LD model, including coding, conserved, regulatory, and LD-related annotations. We analytically estimate posterior mean causal effect sizes and then use cross-validation to regularize these estimates, improving prediction accuracy for sparse architectures. We applied LDpred-funct to predict 21 highly heritable traits in the UK Biobank (avg N = 373 K as training data). LDpred-funct attained a +4.6% relative improvement in average prediction accuracy (avg prediction R2 = 0.144; highest R2 = 0.413 for height) compared to SBayesR (the best method that does not incorporate functional information). For height, meta-analyzing training data from UK Biobank and 23andMe cohorts (N = 1107 K) increased prediction R2 to 0.431. Our results show that incorporating functional priors improves polygenic prediction accuracy, consistent with the functional architecture of complex traits.

Comparison of in silico strategies to prioritize rare genomic variants impacting RNA splicing for the diagnosis of genomic disorders

The development of computational methods to assess pathogenicity of pre-messenger RNA splicing variants is critical for diagnosis of human disease. We assessed the capability of eight algorithms, and a consensus approach, to prioritize 249 variants of uncertain significance (VUSs) that underwent splicing functional analyses. The capability of algorithms to differentiate VUSs away from the immediate splice site as being 'pathogenic' or 'benign' is likely to have substantial impact on diagnostic testing. We show that SpliceAI is the best single strategy in this regard, but that combined usage of tools using a weighted approach can increase accuracy further. We incorporated prioritization strategies alongside diagnostic testing for rare disorders. We show that 15% of 2783 referred individuals carry rare variants expected to impact splicing that were not initially identified as 'pathogenic' or 'likely pathogenic'; one in five of these cases could lead to new or refined diagnoses.

Biallelic variants in genes previously associated with dominant inheritance: CACNA1A, RET and SLC20A2

A subset of families with co-dominant or recessive inheritance has been described in several genes previously associated with dominant inheritance. Those recessive families displayed similar, more severe, or even completely different phenotypes to their dominant counterparts. We report the first patients harboring homozygous disease-related variants in three genes that were previously associated with dominant inheritance: a loss-of-function variant in the CACNA1A gene and two missense variants in the RET and SLC20A2 genes, respectively. All patients presented with a more severe clinical phenotype than the corresponding typical dominant form. We suggest that co-dominant or recessive inheritance for these three genes could explain the phenotypic differences from those documented in their cognate dominant phenotypes. Our results reinforce that geneticists should be aware of the possible different forms of inheritance in genes when WES variant interpretation is performed. We also evidence the need to refine phenotypes and inheritance patterns associated with genes in order to avoid failures during WES analysis and thus, raising the WES diagnostic capacity in the benefit of patients.

Biallelic P4HTM variants associated with HIDEA syndrome and mitochondrial respiratory chain complex I deficiency

We report a patient with profound congenital hypotonia, central hypoventilation, poor visual behaviour with retinal hypopigmentation, and significantly decreased mitochondrial respiratory chain complex I activity in muscle, who died at 7 months of age having made minimal developmental progress. Biallelic predicted truncating P4HTM variants were identified following trio whole-genome sequencing, consistent with a diagnosis of hypotonia, hypoventilation, intellectual disability, dysautonomia, epilepsy and eye abnormalities (HIDEA) syndrome. Very few patients with HIDEA syndrome have been reported previously and mitochondrial abnormalities were observed in three of four previous cases who had a muscle biopsy, suggesting the possibility that HIDEA syndrome represents a primary mitochondrial disorder. P4HTM encodes a transmembrane prolyl 4-hydroxylase with putative targets including hypoxia inducible factors, RNA polymerase II and activating transcription factor 4, which has been implicated in the integrated stress response observed in cell and animal models of mitochondrial disease, and may explain the mitochondrial dysfunction observed in HIDEA syndrome.

Consolidating biallelic SDHD variants as a cause of mitochondrial complex II deficiency

Isolated mitochondrial complex II deficiency is a rare cause of mitochondrial respiratory chain disease. To date biallelic variants in three genes encoding mitochondrial complex II molecular components have been unequivocally associated with mitochondrial disease (SDHA/SDHB/SDHAF1). Additionally, variants in one further complex II component (SDHD) have been identified as a candidate cause of isolated mitochondrial complex II deficiency in just two unrelated affected individuals with clinical features consistent with mitochondrial disease, including progressive encephalomyopathy and lethal infantile cardiomyopathy. We present clinical and genomic investigations in four individuals from an extended Palestinian family with clinical features consistent with an autosomal recessive mitochondrial complex II deficiency, in which our genomic studies identified a homozygous NM_003002.3:c.[205 G > A];[205 G > A];p.[(Glu69Lys)];[(Glu69Lys)] SDHD variant as the likely cause. Reviewing previously published cases, these findings consolidate disruption of SDHD function as a cause of mitochondrial complex II deficiency and further define the phenotypic spectrum associated with SDHD gene variants.

Analyses of oligodontia phenotypes and genetic etiologies

Oligodontia is the congenital absence of six or more teeth and comprises the more severe forms of tooth agenesis. Many genes have been implicated in the etiology of tooth agenesis, which is highly variable in its clinical presentation. The purpose of this study was to identify associations between genetic mutations and clinical features of oligodontia patients. An online systematic search of papers published from January 1992 to June 2021 identified 381 oligodontia cases meeting the eligibility criteria of causative gene mutation, phenotype description, and radiographic records. Additionally, ten families with oligodontia were recruited and their genetic etiologies were determined by whole-exome sequence analyses. We identified a novel mutation in WNT10A (c.99_105dup) and eight previously reported mutations in WNT10A (c.433 G > A; c.682 T > A; c.318 C > G; c.511.C > T; c.321 C > A), EDAR (c.581 C > T), and LRP6 (c.1003 C > T, c.2747 G > T). Collectively, 20 different causative genes were implicated among those 393 cases with oligodontia. For each causative gene, the mean number of missing teeth per case and the frequency of teeth missing at each position were calculated. Genotype-phenotype correlation analysis indicated that molars agenesis is more likely linked to PAX9 mutations, mandibular first premolar agenesis is least associated with PAX9 mutations. Mandibular incisors and maxillary lateral incisor agenesis are most closely linked to EDA mutations.

Cystin genetic variants cause autosomal recessive polycystic kidney disease associated with altered Myc expression

Mutation of the Cys1 gene underlies the renal cystic disease in the Cys1^cpk/cpk (cpk) mouse that phenocopies human autosomal recessive polycystic kidney disease (ARPKD). Cystin, the protein product of Cys1, is expressed in the primary apical cilia of renal ductal epithelial cells. In previous studies, we showed that cystin regulates Myc expression via interaction with the tumor suppressor, necdin. Here, we demonstrate rescue of the cpk renal phenotype by kidney-specific expression of a cystin-GFP fusion protein encoded by a transgene integrated into the Rosa26 locus. In addition, we show that expression of the cystin-GFP fusion protein in collecting duct cells down-regulates expression of Myc in cpk kidneys. Finally, we report the first human patient with an ARPKD phenotype due to homozygosity for a deleterious splicing variant in CYS1. These findings suggest that mutations in Cys1/CYS1 cause an ARPKD phenotype in mouse and human, respectively, and that the renal cystic phenotype in the mouse is driven by overexpression of the Myc proto-oncogene.

ATG4D is the main ATG8 delipidating enzyme in mammalian cells and protects against cerebellar neurodegeneration

Despite the great advances in autophagy research in the last years, the specific functions of the four mammalian Atg4 proteases (ATG4A-D) remain unclear. In yeast, Atg4 mediates both Atg8 proteolytic activation, and its delipidation. However, it is not clear how these two roles are distributed along the members of the ATG4 family of proteases. We show that these two functions are preferentially carried out by distinct ATG4 proteases, being ATG4D the main delipidating enzyme. In mammalian cells, ATG4D loss results in accumulation of membrane-bound forms of mATG8s, increased cellular autophagosome number and reduced autophagosome average size. In mice, ATG4D loss leads to cerebellar neurodegeneration and impaired motor coordination caused by alterations in trafficking/clustering of GABA_A receptors. We also show that human gene variants of ATG4D associated with neurodegeneration are not able to fully restore ATG4D deficiency, highlighting the neuroprotective role of ATG4D in mammals.

Mesomelic dysplasias associated with the HOXD locus are caused by regulatory reallocations

Human families with chromosomal rearrangements at 2q31, where the human HOXD locus maps, display mesomelic dysplasia, a severe shortening and bending of the limb. In mice, the dominant Ulnaless inversion of the HoxD cluster produces a similar phenotype suggesting the same origin for these malformations in humans and mice. Here we engineer 1 Mb inversion including the HoxD gene cluster, which positioned Hoxd13 close to proximal limb enhancers. Using this model, we show that these enhancers contact and activate Hoxd13 in proximal cells, inducing the formation of mesomelic dysplasia. We show that a secondary Hoxd13 null mutation in-cis with the inversion completely rescues the alterations, demonstrating that ectopic HOXD13 is directly responsible for this bone anomaly. Single-cell expression analysis and evaluation of HOXD13 binding sites suggests that the phenotype arises primarily by acting through genes normally controlled by HOXD13 in distal limb cells. Altogether, these results provide a conceptual and mechanistic framework to understand and unify the molecular origins of human mesomelic dysplasia associated with 2q31.

The role of genetics in neurodegenerative dementia: a large cohort study in South China

Neurodegenerative dementias are a group of diseases with highly heterogeneous pathology and complicated etiology. There exist potential genetic component overlaps between different neurodegenerative dementias. Here, 1795 patients with neurodegenerative dementias from South China were enrolled, including 1592 with Alzheimer's disease (AD), 110 with frontotemporal dementia (FTD), and 93 with dementia with Lewy bodies (DLB). Genes targeted sequencing analysis were performed. According to the American College of Medical Genetics (ACMG) guidelines, 39 pathogenic/likely pathogenic (P/LP) variants were identified in 47 unrelated patients in 14 different genes, including PSEN1, PSEN2, APP, MAPT, GRN, CHCHD10, TBK1, VCP, HTRA1, OPTN, SQSTM1, SIGMAR1, and abnormal repeat expansions in C9orf72 and HTT. Overall, 33.3% (13/39) of the variants were novel, the identified P/LP variants were seen in 2.2% (35/1592) and 10.9% (12/110) of AD and FTD cases, respectively. The overall molecular diagnostic rate was 2.6%. Among them, PSEN1 was the most frequently mutated gene (46.8%, 22/47), followed by PSEN2 and APP. Additionally, the age at onset of patients with P/LP variants (51.4 years), ranging from 30 to 83 years, was ~10 years earlier than those without P/LP variants (p < 0.05). This study sheds insight into the genetic spectrum and clinical manifestations of neurodegenerative dementias in South China, further expands the existing repertoire of P/LP variants involved in known dementia-associated genes. It provides a new perspective for basic research on genetic pathogenesis and novel guiding for clinical practice of neurodegenerative dementia.

Clinical delineation of SETBP1 haploinsufficiency disorder

SETBP1 haploinsufficiency disorder (MIM#616078) is caused by haploinsufficiency of SETBP1 on chromosome 18q12.3, but there has not yet been any systematic evaluation of the major features of this monogenic syndrome, assessing penetrance and expressivity. We describe the first comprehensive study to delineate the associated clinical phenotype, with findings from 34 individuals, including 24 novel cases, all of whom have a SETBP1 loss-of-function variant or single (coding) gene deletion, confirmed by molecular diagnostics. The most commonly reported clinical features included mild motor developmental delay, speech impairment, intellectual disability, hypotonia, vision impairment, attention/concentration deficits, and hyperactivity. Although there is a mild overlap in certain facial features, the disorder does not lead to a distinctive recognizable facial gestalt. As well as providing insight into the clinical spectrum of SETBP1 haploinsufficiency disorder, this reports puts forward care recommendations for patient management.

Neutrophil-derived reactive oxygen species promote tumor colonization

A single-nucleotide polymorphism of neutrophil cytosolic factor 1 (Ncf1), leading to an impaired generation of reactive oxygen species (ROS), is a causative genetic factor for autoimmune disease. To study a possible tumor protection effect by the Ncf1 mutation in a manner dependent on cell types, we used experimental mouse models of lung colonization assay by B16F10 melanoma cells. We observed fewer tumor foci in Ncf1 mutant mice, irrespective of αβT, γδT, B-cell deficiencies, or of a functional Ncf1 expression in CD68-positive monocytes/macrophages. The susceptibility to tumor colonization was restored by the human S100A8 (MRP8) promoter directing a functional Ncf1 expression to granulocytes. This effect was associated with an increase of both ROS and interleukin 1 beta (IL-1β) production from lung neutrophils. Moreover, neutrophil depletion by anti-Ly6G antibodies increased tumor colonization in wild type but failed in the Ncf1 mutant mice. In conclusion, tumor colonization is counteracted by ROS-activated and IL-1β-secreting tissue neutrophils.

MCM complex members MCM3 and MCM7 are associated with a phenotypic spectrum from Meier-Gorlin syndrome to lipodystrophy and adrenal insufficiency

The MCM2-7 helicase is a heterohexameric complex with essential roles as part of both the pre-replication and pre-initiation complexes in the early stages of DNA replication. Meier-Gorlin syndrome, a rare primordial dwarfism, is strongly associated with disruption to the pre-replication complex, including a single case described with variants in MCM5. Conversely, a biallelic pathogenic variant in MCM4 underlies immune deficiency with growth retardation, features also seen in individuals with pathogenic variants in other pre-initiation complex encoding genes such as GINS1, MCM10, and POLE. Through exome and chromium genome sequencing, supported by functional studies, we identify biallelic pathogenic variants in MCM7 and a strong candidate biallelic pathogenic variant in MCM3. We confirm variants in MCM7 are deleterious and through interfering with MCM complex formation, impact efficiency of S phase progression. The associated phenotypes are striking; one patient has typical Meier-Gorlin syndrome, whereas the second case has a multi-system disorder with neonatal progeroid appearance, lipodystrophy and adrenal insufficiency. We provide further insight into the developmental complexity of disrupted MCM function, highlighted by two patients with a similar variant profile in MCM7 but disparate clinical features. Our results build on other genetic findings linked to disruption of the pre-replication and pre-initiation complexes, and the replisome, and expand the complex clinical genetics landscape emerging due to disruption of DNA replication.

Broadening INPP5E phenotypic spectrum: detection of rare variants in syndromic and non-syndromic IRD

Pathogenic variants in INPP5E cause Joubert syndrome (JBTS), a ciliopathy with retinal involvement. However, despite sporadic cases in large cohort sequencing studies, a clear association with non-syndromic inherited retinal degenerations (IRDs) has not been made. We validate this association by reporting 16 non-syndromic IRD patients from ten families with bi-allelic mutations in INPP5E. Additional two patients showed early onset IRD with limited JBTS features. Detailed phenotypic description for all probands is presented. We report 14 rare INPP5E variants, 12 of which have not been reported in previous studies. We present tertiary protein modeling and analyze all INPP5E variants for deleteriousness and phenotypic correlation. We observe that the combined impact of INPP5E variants in JBTS and non-syndromic IRD patients does not reveal a clear genotype-phenotype correlation, suggesting the involvement of genetic modifiers. Our study cements the wide phenotypic spectrum of INPP5E disease, adding proof that sequence defects in this gene can lead to early-onset non-syndromic IRD.

Novel variants in the stem cell niche factor WNT2B define the disease phenotype as a congenital enteropathy with ocular dysgenesis

WNT2B is a member of the Wnt family, a group of signal transduction proteins involved in embryologic development and stem cell renewal and maintenance. We recently reported homozygous nonsense variants in WNT2B in three individuals with severe, neonatal-onset diarrhea, and intestinal failure. Here we present a fourth case, from a separate family, with neonatal diarrhea associated with novel compound heterozygous WNT2B variants. One of the two variants was a frameshift variant (c.423del [p.Phe141fs]), while the other was a missense change (c.722 G > A [p.G241D]) that we predict through homology modeling to be deleterious, disrupting post-translational acylation. This patient presented as a neonate with severe diet-induced (osmotic) diarrhea and growth failure resulting in dependence on parenteral nutrition. Her gastrointestinal histology revealed abnormal cellular architecture particularly in the stomach and colon, including oxyntic atrophy, abnormal distribution of enteroendocrine cells, and a paucity of colonic crypt glands. In addition to her gastrointestinal findings, she had bilateral corneal clouding and atypical genital development later identified as a testicular 46,XX difference/disorder of sexual development. Upon review of the previously reported cases, two others also had anterior segment ocular anomalies though none had atypical genital development. This growing case series suggests that variants in WNT2B are associated with an oculo-intestinal (and possibly gonadal) syndrome, due to the protein's putative involvement in multiple developmental and stem cell maintenance pathways.

Independent and cumulative coeliac disease-susceptibility loci are associated with distinct disease phenotypes

The phenotype of coeliac disease varies considerably for incompletely understood reasons. We investigated whether established coeliac disease susceptibility variants (SNPs) are individually or cumulatively associated with distinct phenotypes. We also tested whether a polygenic risk score (PRS) based on genome-wide associated (GWA) data could explain the phenotypic variation. The phenotypic association of 39 non-HLA coeliac disease SNPs was tested in 625 thoroughly phenotyped coeliac disease patients and 1817 controls. To assess their cumulative effects a weighted genetic risk score (wGRS39) was built, and stratified by tertiles. In our PRS model in cases, we took the summary statistics from the largest GWA study in coeliac disease and tested their association at eight P value thresholds (PT) with phenotypes. Altogether ten SNPs were associated with distinct phenotypes after correction for multiple testing (PEMP2 ≤ 0.05). The TLR7/TLR8 locus was associated with disease onset before and the SH2B3/ATXN2, ITGA4/UBE2E3 and IL2/IL21 loci after 7 years of age. The latter three loci were associated with a more severe small bowel mucosal damage and SH2B3/ATXN2 with type 1 diabetes. Patients at the highest wGRS39 tertiles had OR > 1.62 for having coeliac disease-related symptoms during childhood, a more severe small bowel mucosal damage, malabsorption and anaemia. PRS was associated only with dermatitis herpetiformis (PT = 0.2, PEMP2 = 0.02). Independent coeliac disease-susceptibility loci are associated with distinct phenotypes, suggesting that genetic factors play a role in determining the disease presentation. Moreover, the increased number of coeliac disease susceptibility SNPs might predispose to a more severe disease course.

Parental mosaicism in Marfan and Ehlers-Danlos syndromes and related disorders

Marfan syndrome (MFS) is a heritable connective tissue disorder (HCTD) caused by pathogenic variants in FBN1 that frequently occur de novo. Although individuals with somatogonadal mosaicisms have been reported with respect to MFS and other HCTD, the overall frequency of parental mosaicism in this pathology is unknown. In an attempt to estimate this frequency, we reviewed all the 333 patients with a disease-causing variant in FBN1. We then used direct sequencing, combined with High Resolution Melting Analysis, to detect mosaicism in their parents, complemented by NGS when a mosaicism was objectivized. We found that (1) the number of apparently de novo events is much higher than the classically admitted number (around 50% of patients and not 25% as expected for FBN1) and (2) around 5% of the FBN1 disease-causing variants were not actually de novo as anticipated, but inherited in a context of somatogonadal mosaicisms revealed in parents from three families. High Resolution Melting Analysis and NGS were more efficient at detecting and evaluating the level of mosaicism compared to direct Sanger sequencing. We also investigated individuals with a causal variant in another gene identified through our "aortic diseases genes" NGS panel and report, for the first time, on an individual with a somatogonadal mosaicism in COL5A1. Our study shows that parental mosaicism is not that rare in Marfan syndrome and should be investigated with appropriate methods given its implications in patient's management.

The genetic landscape of polycystic kidney disease in Ireland

Polycystic kidney diseases (PKDs) comprise the most common Mendelian forms of renal disease. It is characterised by the development of fluid-filled renal cysts, causing progressive loss of kidney function, culminating in the need for renal replacement therapy or kidney transplant. Ireland represents a valuable region for the genetic study of PKD, as family sizes are traditionally large and the population relatively homogenous. Studying a cohort of 169 patients, we describe the genetic landscape of PKD in Ireland for the first time, compare the clinical features of patients with and without a molecular diagnosis and correlate disease severity with autosomal dominant pathogenic variant type. Using a combination of molecular genetic tools, including targeted next-generation sequencing, we report diagnostic rates of 71-83% in Irish PKD patients, depending on which variant classification guidelines are used (ACMG or Mayo clinic respectively). We have catalogued a spectrum of Irish autosomal dominant PKD pathogenic variants including 36 novel variants. We illustrate how apparently unrelated individuals carrying the same autosomal dominant pathogenic variant are highly likely to have inherited that variant from a common ancestor. We highlight issues surrounding the implementation of the ACMG guidelines for variant pathogenicity interpretation in PKD, which have important implications for clinical genetics.

Germline variants at SOHLH2 influence multiple myeloma risk

Multiple myeloma (MM) is caused by the uncontrolled, clonal expansion of plasma cells. While there is epidemiological evidence for inherited susceptibility, the molecular basis remains incompletely understood. We report a genome-wide association study totalling 5,320 cases and 422,289 controls from four Nordic populations, and find a novel MM risk variant at SOHLH2 at 13q13.3 (risk allele frequency = 3.5%; odds ratio = 1.38; P = 2.2 × 10-14). This gene encodes a transcription factor involved in gametogenesis that is normally only weakly expressed in plasma cells. The association is represented by 14 variants in linkage disequilibrium. Among these, rs75712673 maps to a genomic region with open chromatin in plasma cells, and upregulates SOHLH2 in this cell type. Moreover, rs75712673 influences transcriptional activity in luciferase assays, and shows a chromatin looping interaction with the SOHLH2 promoter. Our work provides novel insight into MM susceptibility.

Deletion of CTCF sites in the SHH locus alters enhancer-promoter interactions and leads to acheiropodia

Acheiropodia, congenital limb truncation, is associated with homozygous deletions in the LMBR1 gene around ZRS, an enhancer regulating SHH during limb development. How these deletions lead to this phenotype is unknown. Using whole-genome sequencing, we fine-mapped the acheiropodia-associated region to 12 kb and show that it does not function as an enhancer. CTCF and RAD21 ChIP-seq together with 4C-seq and DNA FISH identify three CTCF sites within the acheiropodia-deleted region that mediate the interaction between the ZRS and the SHH promoter. This interaction is substituted with other CTCF sites centromeric to the ZRS in the disease state. Mouse knockouts of the orthologous 12 kb sequence have no apparent abnormalities, showcasing the challenges in modelling CTCF alterations in animal models due to inherent motif differences between species. Our results show that alterations in CTCF motifs can lead to a Mendelian condition due to altered enhancer-promoter interactions.

Deficiency of TMEM53 causes a previously unknown sclerosing bone disorder by dysregulation of BMP-SMAD signaling

Bone formation represents a heritable trait regulated by many signals and complex mechanisms. Its abnormalities manifest themselves in various diseases, including sclerosing bone disorder (SBD). Exploration of genes that cause SBD has significantly improved our understanding of the mechanisms that regulate bone formation. Here, we discover a previously unknown type of SBD in four independent families caused by bi-allelic loss-of-function pathogenic variants in TMEM53, which encodes a nuclear envelope transmembrane protein. Tmem53-/- mice recapitulate the human skeletal phenotypes. Analyses of the molecular pathophysiology using the primary cells from the Tmem53-/- mice and the TMEM53 knock-out cell lines indicates that TMEM53 inhibits BMP signaling in osteoblast lineage cells by blocking cytoplasm-nucleus translocation of BMP2-activated Smad proteins. Pathogenic variants in the patients impair the TMEM53-mediated blocking effect, thus leading to overactivated BMP signaling that promotes bone formation and contributes to the SBD phenotype. Our results establish a previously unreported SBD entity (craniotubular dysplasia, Ikegawa type) and contribute to a better understanding of the regulation of BMP signaling and bone formation.

Pathogenic variants causing ABL1 malformation syndrome cluster in a myristoyl-binding pocket and increase tyrosine kinase activity

ABL1 is a proto-oncogene encoding a nonreceptor tyrosine kinase, best known in the somatic BCR-ABL fusion gene associated with chronic myeloid leukaemia. Recently, germline missense variants in ABL1 have been found to cause an autosomal dominant developmental syndrome with congenital heart disease, skeletal malformations and characteristic facies. Here, we describe a series of six new unrelated individuals with heterozygous missense variants in ABL1 (including four novel variants) identified via whole exome sequencing. All the affected individuals in this series recapitulate the phenotype of the ABL1 developmental syndrome and additionally we affirm that hearing impairment is a common feature of the condition. Four of the variants cluster in the myristoyl-binding pocket of ABL1, a region critical for auto-inhibitory regulation of the kinase domain. Bio-informatic analysis of transcript-wide conservation and germline/somatic variation reveals that this pocket region is subject to high missense constraint and evolutionary conservation. Functional work to investigate ABL1 kinase activity in vitro by transient transfection of HEK293T cells with variant ABL1 plasmid constructs revealed increased phosphorylation of ABL1-specific substrates compared to wild-type. The increased tyrosine kinase activity was suppressed by imatinib treatment. This case series of six new patients with germline heterozygous ABL1 missense variants further delineates the phenotypic spectrum of this condition and recognises microcephaly as a common finding. Our analysis supports an ABL1 gain-of-function mechanism due to loss of auto-inhibition, and demonstrates the potential for pharmacological inhibition using imatinib.

Beyond DNA repair and chromosome instability-Fanconi anaemia as a cellular senescence-associated syndrome

Fanconi anaemia (FA) is the most frequent inherited bone marrow failure syndrome, due to mutations in genes encoding proteins involved in replication fork protection, DNA interstrand crosslink repair and replication rescue through inducing double-strand break repair and homologous recombination. Clinically, FA is characterised by aplastic anaemia, congenital defects and cancer predisposition. In in vitro studies, FA cells presented hallmarks defining senescent cells, including p53-p21 axis activation, altered telomere length, mitochondrial dysfunction, chromatin alterations, and a pro-inflammatory status. Senescence is a programme leading to proliferation arrest that is involved in different physiological contexts, such as embryogenesis, tissue remodelling and repair and guarantees tumour suppression activity. However, senescence can become a driving force for developmental abnormalities, aging and cancer. Herein, we summarise the current knowledge in the field to highlight the mutual relationships between FA and senescence that lead us to consider FA not only as a DNA repair and chromosome fragility syndrome but also as a "senescence syndrome".

The orthopedic characterization of cfap298tm304 mutants validate zebrafish to faithfully model human AIS

Cerebrospinal fluid (CSF) circulation relies on the beating of motile cilia projecting in the lumen of the brain and spinal cord cavities Mutations in genes involved in cilia motility disturb cerebrospinal fluid circulation and result in scoliosis-like deformities of the spine in juvenile zebrafish. However, these defects in spine alignment have not been validated with clinical criteria used to diagnose adolescent idiopathic scoliosis (AIS). The aim of this study was to describe, using orthopaedic criteria the spinal deformities of a zebrafish mutant model of AIS targeting a gene involved in cilia polarity and motility, cfap298tm304. The zebrafish mutant line cfap298tm304, exhibiting alteration of CSF flow due to defective cilia motility, was raised to the juvenile stage. The analysis of mutant animals was based on micro-computed tomography (micro-CT), which was conducted in a QUANTUM FX CALIPER, with a 59 µm-30 mm protocol. 63% of the cfap298tm304 zebrafish analyzed presented a three-dimensional deformity of the spine, that was evolutive during the juvenile phase, more frequent in females, with a right convexity, a rotational component and involving at least one dislocation. We confirm here that cfap298tm304 scoliotic individuals display a typical AIS phenotype, with orthopedic criteria mirroring patient's diagnosis.

Expression analysis of selected genes involved in tryptophan metabolic pathways in Egyptian children with Autism Spectrum Disorder and learning disabilities

Autism Spectrum Disorder (ASD) and learning disabilities are neurodevelopmental disabilities characterized by dramatically increasing incidence rates, yet the exact etiology for these disabilities is not identified. Impairment in tryptophan metabolism has been suggested to participate in the pathogenesis of ASD, however, further validation of its involvement is required. Additionally, its role in learning disabilities is still uninvestigated. Our objective was to evaluate some aspects of tryptophan metabolism in ASD children (N = 45) compared to children with learning disabilities (N = 44) and healthy controls (N = 40) by measuring the expression levels of the MAOA, HAAO and AADAT genes using real-time RT-qPCR. We also aimed to correlate the expression patterns of these genes with parental ages at the time of childbirth, levels of serum iron, and vitamin D3 and zinc/copper ratio, as possible risk factors for ASD. Results demonstrated a significant decrease in the expression of the selected genes within ASD children (p < 0.001) relative to children with learning disabilities and healthy controls, which significantly associated with the levels of our targeted risk factors (p < 0.05) and negatively correlated to ASD scoring (p < 0.001). In conclusion, this study suggests that the expression of the MAOA, HAAO and AADAT genes may underpin the pathophysiology of ASD.

The TOMM40 '523' polymorphism in disease risk and age of symptom onset in two independent cohorts of Parkinson's disease

Abnormal mitochondrial function is a key process in the pathogenesis of Parkinson's disease (PD). The central pore-forming protein TOM40 of the mitochondria is encoded by the translocase of outer mitochondrial membrane 40 homologue gene (TOMM40). The highly variant '523' poly-T repeat is associated with age-related cognitive decline and age of onset in Alzheimer's disease, but whether it plays a role in modifying the risk or clinical course of PD it yet to be elucidated. The TOMM40 '523' allele length was determined in 634 people with PD and 422 healthy controls from an Australian cohort and the Parkinson's Progression Markers Initiative (PPMI) cohort, using polymerase chain reaction or whole genome sequencing analysis. Genotype and allele frequencies of TOMM40 '523' and APOE ε did not differ significantly between the cohorts. Analyses revealed TOMM40 '523' allele groups were not associated with disease risk, while considering APOE ε genotype. Regression analyses revealed the TOMM40 S/S genotype was associated with a significantly later age of symptom onset in the PPMI PD cohort, but not after correction for covariates, or in the Australian cohort. Whilst variation in the TOMM40 '523' polymorphism was not associated with PD risk, the possibility that it may be a modifying factor for age of symptom onset warrants further investigation in other PD populations.

Investigating the relationship between the SNCA gene and cognitive abilities in idiopathic Parkinson's disease using machine learning

Cognitive impairments are prevalent in Parkinson's disease (PD), but the underlying mechanisms of their development are unknown. In this study, we aimed to predict global cognition (GC) in PD with machine learning (ML) using structural neuroimaging, genetics and clinical and demographic characteristics. As a post-hoc analysis, we aimed to explore the connection between novel selected features and GC more precisely and to investigate whether this relationship is specific to GC or is driven by specific cognitive domains. 101 idiopathic PD patients had a cognitive assessment, structural MRI and blood draw. ML was performed on 102 input features including demographics, cortical thickness and subcortical measures, and several genetic variants (APOE, MAPT, SNCA, etc.). Using the combination of RRELIEFF and Support Vector Regression, 11 features were found to be predictive of GC including sex, rs894280, Edinburgh Handedness Inventory, UPDRS-III, education, five cortical thickness measures (R-parahippocampal, L-entorhinal, R-rostral anterior cingulate, L-middle temporal, and R-transverse temporal), and R-caudate volume. The rs894280 of SNCA gene was selected as the most novel finding of ML. Post-hoc analysis revealed a robust association between rs894280 and GC, attention, and visuospatial abilities. This variant indicates a potential role for the SNCA gene in cognitive impairments of idiopathic PD.

A generic assay for the identification of splicing variants that induce nonsense-mediated decay in Pompe disease

DNA variants affecting mRNA expression and processing in genetic diseases are often missed or poorly characterized. We previously reported a generic assay to identify variants that affect mRNA expression and splicing in Pompe disease, a monogenic disorder caused by deficiency of acid α-glucosidase (GAA). However, this assay could miss mRNA that is subjected to degradation. Here, we inhibited mRNA degradation using cycloheximide and performed unbiased splicing analysis of all GAA exons using exon flanking RT-PCR and exon internal RT-qPCR. In four patients that were suspected of harboring splicing variants but for which aberrant splicing could not be detected in normally growing cells, we detected a total of 10 novel splicing events in cells treated with cycloheximide. In addition, we found that sequences of GAA introns 6 and 12 were naturally included in a subset of transcripts from patients and healthy controls, indicating inefficient canonical splicing. Identification of aberrant splicing caused by the common Asian variant c.546G>T allowed the development of an antisense oligonucleotide that promoted canonical GAA pre-mRNA splicing and elevated GAA enzymatic activity. Our results indicate that this extended generic splicing assay allows the detection of aberrant splicing in cases of mRNA degradation to enable functional analysis of unknown splicing variants and the development of targeted treatment options.

Accurate imputation of human leukocyte antigens with CookHLA

The recent development of imputation methods enabled the prediction of human leukocyte antigen (HLA) alleles from intergenic SNP data, allowing studies to fine-map HLA for immune phenotypes. Here we report an accurate HLA imputation method, CookHLA, which has superior imputation accuracy compared to previous methods. CookHLA differs from other approaches in that it locally embeds prediction markers into highly polymorphic exons to account for exonic variability, and in that it adaptively learns the genetic map within MHC from the data to facilitate imputation. Our benchmarking with real datasets shows that our method achieves high imputation accuracy in a wide range of scenarios, including situations where the reference panel is small or ethnically unmatched.

Whole exome sequencing uncovered highly penetrant recessive mutations for a spectrum of rare genetic pediatric diseases in Bangladesh

Collectively, rare genetic diseases affect a significant number of individuals worldwide. In this study, we have conducted whole-exome sequencing (WES) and identified underlying pathogenic or likely pathogenic variants in five children with rare genetic diseases. We present evidence for disease-causing autosomal recessive variants in a range of disease-associated genes such as DHH-associated 46,XY gonadal dysgenesis (GD) or 46,XY sex reversal 7, GNPTAB-associated mucolipidosis II alpha/beta (ML II), BBS1-associated Bardet-Biedl Syndrome (BBS), SURF1-associated Leigh Syndrome (LS) and AP4B1-associated spastic paraplegia-47 (SPG47) in unrelated affected members from Bangladesh. Our analysis pipeline detected three homozygous mutations, including a novel c. 863 G > C (p.Pro288Arg) variant in DHH, and two compound heterozygous variants, including two novel variants: c.2972dupT (p.Met991Ilefs*) in GNPTAB and c.229 G > C (p.Gly77Arg) in SURF1. All mutations were validated by Sanger sequencing. Collectively, this study adds to the genetic heterogeneity of rare genetic diseases and is the first report elucidating the genetic profile of (consanguineous and nonconsanguineous) rare genetic diseases in the Bangladesh population.

Impaired eIF5A function causes a Mendelian disorder that is partially rescued in model systems by spermidine

The structure of proline prevents it from adopting an optimal position for rapid protein synthesis. Poly-proline-tract (PPT) associated ribosomal stalling is resolved by highly conserved eIF5A, the only protein to contain the amino acid hypusine. We show that de novo heterozygous EIF5A variants cause a disorder characterized by variable combinations of developmental delay, microcephaly, micrognathia and dysmorphism. Yeast growth assays, polysome profiling, total/hypusinated eIF5A levels and PPT-reporters studies reveal that the variants impair eIF5A function, reduce eIF5A-ribosome interactions and impair the synthesis of PPT-containing proteins. Supplementation with 1 mM spermidine partially corrects the yeast growth defects, improves the polysome profiles and restores expression of PPT reporters. In zebrafish, knockdown eif5a partly recapitulates the human phenotype that can be rescued with 1 µM spermidine supplementation. In summary, we uncover the role of eIF5A in human development and disease, demonstrate the mechanistic complexity of EIF5A-related disorder and raise possibilities for its treatment.

Deubiquitylases in developmental ubiquitin signaling and congenital diseases

Metazoan development from a one-cell zygote to a fully formed organism requires complex cellular differentiation and communication pathways. To coordinate these processes, embryos frequently encode signaling information with the small protein modifier ubiquitin, which is typically attached to lysine residues within substrates. During ubiquitin signaling, a three-step enzymatic cascade modifies specific substrates with topologically unique ubiquitin modifications, which mediate changes in the substrate's stability, activity, localization, or interacting proteins. Ubiquitin signaling is critically regulated by deubiquitylases (DUBs), a class of ~100 human enzymes that oppose the conjugation of ubiquitin. DUBs control many essential cellular functions and various aspects of human physiology and development. Recent genetic studies have identified mutations in several DUBs that cause developmental disorders. Here we review principles controlling DUB activity and substrate recruitment that allow these enzymes to regulate ubiquitin signaling during development. We summarize key mechanisms of how DUBs control embryonic and postnatal differentiation processes, highlight developmental disorders that are caused by mutations in particular DUB members, and describe our current understanding of how these mutations disrupt development. Finally, we discuss how emerging tools from human disease genetics will enable the identification and study of novel congenital disease-causing DUBs.

NAA10 p.(N101K) disrupts N-terminal acetyltransferase complex NatA and is associated with developmental delay and hemihypertrophy

Nearly half of all human proteins are acetylated at their N-termini by the NatA N-terminal acetyltransferase complex. NAA10 is evolutionarily conserved as the catalytic subunit of NatA in complex with NAA15, but may also have NatA-independent functions. Several NAA10 variants are associated with genetic disorders. The phenotypic spectrum includes developmental delay, intellectual disability, and cardiac abnormalities. Here, we have identified the previously undescribed NAA10 c.303C>A and c.303C>G p.(N101K) variants in two unrelated girls. These girls have developmental delay, but they both also display hemihypertrophy a feature normally not observed or registered among these cases. Functional studies revealed that NAA10 p.(N101K) is completely impaired in its ability to bind NAA15 and to form an enzymatically active NatA complex. In contrast, the integrity of NAA10 p.(N101K) as a monomeric acetyltransferase is intact. Thus, this NAA10 variant may represent the best example of the impact of NatA mediated N-terminal acetylation, isolated from other potential NAA10-mediated cellular functions and may provide important insights into the phenotypes observed in individuals expressing pathogenic NAA10 variants.

Further confirmation of the association of SLC12A2 with non-syndromic autosomal-dominant hearing impairment

Congenital hearing impairment (HI) is genetically heterogeneous making its genetic diagnosis challenging. Investigation of novel HI genes and variants will enhance our understanding of the molecular mechanisms and to aid genetic diagnosis. We performed exome sequencing and analysis using DNA samples from affected members of two large families from Ghana and Pakistan, segregating autosomal-dominant (AD) non-syndromic HI (NSHI). Using in silico approaches, we modeled and evaluated the effect of the likely pathogenic variants on protein structure and function. We identified two likely pathogenic variants in SLC12A2, c.2935G>A:p.(E979K) and c.2939A>T:p.(E980V), which segregate with NSHI in a Ghanaian and Pakistani family, respectively. SLC12A2 encodes an ion transporter crucial in the homeostasis of the inner ear endolymph and has recently been reported to be implicated in syndromic and non-syndromic HI. Both variants were mapped to alternatively spliced exon 21 of the SLC12A2 gene. Exon 21 encodes for 17 residues in the cytoplasmatic tail of SLC12A2, is highly conserved between species, and preferentially expressed in cochlear tissues. A review of previous studies and our current data showed that out of ten families with either AD non-syndromic or syndromic HI, eight (80%) had variants within the 17 amino acid residue region of exon 21 (48 bp), suggesting that this alternate domain is critical to the transporter activity in the inner ear. The genotypic spectrum of SLC12A2 was expanded and the involvement of SLC12A2 in ADNSHI was confirmed. These results also demonstrate the role that SLC12A2 plays in ADNSHI in diverse populations including sub-Saharan Africans.

Modelling intestinal inflammation and infection using 'mini-gut' organoids

In 2020, major advances to the understanding of gastrointestinal inflammatory and infectious disease have been made using ‘mini-gut’ organoids. Key findings include the discovery of somatic inflammatory gene mutations in ulcerative colitis epithelium, a unique mutational signature in colorectal cancer caused by genotoxic Escherichia coli, and infection of intestinal organoids by SARS-CoV-2.

Unique roles of rare variants in the genetics of complex diseases in humans

Genome-wide association studies have identified >10,000 genetic variants associated with various phenotypes and diseases. Although the majority are common variants, rare variants with >0.1% of minor allele frequency have been investigated by imputation and using disease-specific custom SNP arrays. Rare variants sequencing analysis mainly revealed have played unique roles in the genetics of complex diseases in humans due to their distinctive features, in contrast to common variants. Unique roles are hypothesis-free evidence for gene causality, a precise target of functional analysis for understanding disease mechanisms, a new favorable target for drug development, and a genetic marker with high disease risk for personalized medicine. As whole-genome sequencing continues to identify more rare variants, the roles associated with rare variants will also increase. However, a better estimation of the functional impact of rare variants across whole genome is needed to enhance their contribution to improvements in human health.

Macular retinal thickness differs markedly in age-related macular degeneration driven by risk polymorphisms on chromosomes 1 and 10

The two most common genetic contributors to age-related macular degeneration (AMD), a leading cause of irreversible vision loss worldwide, are variants associated with CFH-CFHR5 on chromosome 1 (Chr1) and ARMS2/HTRA1 on chromosome 10 (Chr10). We sought to determine if risk and protective variants associated with these two loci drive differences in macular retinal thickness prior and subsequent to the onset of clinically observable signs of AMD. We considered 299 individuals (547 eyes) homozygous for risk variants or haplotypes on Chr1 or Chr10 exclusively (Chr1-risk and Chr10-risk, respectively) or homozygous for a neutral haplotype (Chr1-neu), for the protective I62 tagged haplotype (Chr1-prot-I62) or for the protection conferring CFHR1/3 deletion haplotype (Chr1-prot-del) on Chr1 without any risk alleles on Chr10. Among eyes with no clinically observable signs of AMD, the deletion of CFHR1/3, which is strongly protective against this disease, is associated with significantly thicker retinas in the perifovea. When controlling for age, Chr10-risk eyes with early or intermediate AMD have thinner retinas as compared to eyes from the Chr1-risk group with similar disease severity. Our analysis indicates that this difference likely results from distinct biological and disease initiation and progression events associated with Chr1- and Chr10-directed AMD.

Diphthamide-deficiency syndrome: a novel human developmental disorder and ribosomopathy

We describe a novel type of ribosomopathy that is defined by deficiency in diphthamidylation of translation elongation factor 2. The ribosomopathy was identified by correlating phenotypes and biochemical properties of previously described patients with diphthamide biosynthesis gene 1 (DPH1) deficiencies with a new patient that carried inactivating mutations in both alleles of the human diphthamide biosynthesis gene 2 (DPH2). The human DPH1 syndrome is an autosomal recessive disorder associated with developmental delay, abnormal head circumference (microcephaly or macrocephaly), short stature, and congenital heart disease. It is defined by variants with reduced functionality of the DPH1 gene observed so far predominantly in consanguineous homozygous patients carrying identical mutant alleles of DPH1. Here we report a child with a very similar phenotype carrying biallelic variants of the human DPH2. The gene products DPH1 and DPH2 are components of a heterodimeric enzyme complex that mediates the first step of the posttranslational diphthamide modification on the nonredundant eukaryotic translation elongation factor 2 (eEF2). Diphthamide deficiency was shown to reduce the accuracy of ribosomal protein biosynthesis. Both DPH2 variants described here severely impair diphthamide biosynthesis as demonstrated in human and yeast cells. This is the first report of a patient carrying compound heterozygous DPH2 loss-of-function variants with a DPH1 syndrome-like phenotype and implicates diphthamide deficiency as the root cause of this patient's clinical phenotype as well as of DPH1-syndrome. These findings define "diphthamide-deficiency syndrome" as a special ribosomopathy due to reduced functionality of components of the cellular machinery for eEF2-diphthamide synthesis.

Neurological disorder-associated genetic variants in individuals with psychogenic nonepileptic seizures

Psychogenic nonepileptic seizures (PNES) are diagnosed in approximately 30% of patients referred to tertiary care epilepsy centers. Little is known about the molecular pathology of PNES, much less about possible underlying genetic factors. We generated whole-exome sequencing and whole-genome genotyping data to identify rare, pathogenic (P) or likely pathogenic (LP) variants in 102 individuals with PNES and 448 individuals with focal (FE) or generalized (GE) epilepsy. Variants were classified for all individuals based on the ACMG-AMP 2015 guidelines. For research purposes only, we considered genes associated with neurological or psychiatric disorders as candidate genes for PNES. We observe in this first genetic investigation of PNES that six (5.88%) individuals with PNES without coexistent epilepsy carry P/LP variants (deletions at 10q11.22-q11.23, 10q23.1-q23.2, distal 16p11.2, and 17p13.3, and nonsynonymous variants in NSD1 and GABRA5). Notably, the burden of P/LP variants among the individuals with PNES was similar and not significantly different to the burden observed in the individuals with FE (3.05%) or GE (1.82%) (PNES vs. FE vs. GE (3 × 2 χ2), P = 0.30; PNES vs. epilepsy (2 × 2 χ2), P = 0.14). The presence of variants in genes associated with monogenic forms of neurological and psychiatric disorders in individuals with PNES shows that genetic factors are likely to play a role in PNES or its comorbidities in a subset of individuals. Future large-scale genetic research studies are needed to further corroborate these interesting findings in PNES.

Diagnostic yield of panel-based genetic testing in syndromic inherited retinal disease

Thirty percent of all inherited retinal disease (IRD) is accounted for by conditions with extra-ocular features. This study aimed to establish the genetic diagnostic pick-up rate for IRD patients with one or more extra-ocular features undergoing panel-based screening in a clinical setting. One hundred and six participants, tested on a gene panel which contained both isolated and syndromic IRD genes, were retrospectively ascertained from the Manchester Genomic Diagnostics Laboratory database spanning 6 years (2012-2017). Phenotypic features were extracted from the clinical notes and classified according to Human Phenotype Ontology; all identified genetic variants were interpreted in accordance to the American College of Medical Genetics and Genomics guidelines. Overall, 49% (n = 52) of patients received a probable genetic diagnosis. A further 6% (n = 6) had a single disease-associated variant in an autosomal recessive disease-relevant gene. Fifty-two percent (n = 55) of patients had a clinical diagnosis at the time of testing. Of these, 71% (n = 39) received a probable genetic diagnosis. By contrast, for those without a provisional clinical diagnosis (n = 51), only 25% (n = 13) received a probable genetic diagnosis. The clinical diagnosis of Usher (n = 33) and Bardet-Biedl syndrome (n = 10) was confirmed in 67% (n = 22) and 80% (n = 8), respectively. The testing diagnostic rate in patients with clinically diagnosed multisystemic IRD conditions was significantly higher than those without one (71% versus 25%; p value < 0.001). The lower pick-up rate in patients without a clinical diagnosis suggests that panel-based approaches are unlikely to be the most effective means of achieving a molecular diagnosis for this group. Here, we suggest that genome-wide approaches (whole exome or genome) are more appropriate.

Transcript expression-aware annotation improves rare variant interpretation

The acceleration of DNA sequencing in samples from patients and population studies has resulted in extensive catalogues of human genetic variation, but the interpretation of rare genetic variants remains problematic. A notable example of this challenge is the existence of disruptive variants in dosage-sensitive disease genes, even in apparently healthy individuals. Here, by manual curation of putative loss-of-function (pLoF) variants in haploinsufficient disease genes in the Genome Aggregation Database (gnomAD)1, we show that one explanation for this paradox involves alternative splicing of mRNA, which allows exons of a gene to be expressed at varying levels across different cell types. Currently, no existing annotation tool systematically incorporates information about exon expression into the interpretation of variants. We develop a transcript-level annotation metric known as the 'proportion expressed across transcripts', which quantifies isoform expression for variants. We calculate this metric using 11,706 tissue samples from the Genotype Tissue Expression (GTEx) project2 and show that it can differentiate between weakly and highly evolutionarily conserved exons, a proxy for functional importance. We demonstrate that expression-based annotation selectively filters 22.8% of falsely annotated pLoF variants found in haploinsufficient disease genes in gnomAD, while removing less than 4% of high-confidence pathogenic variants in the same genes. Finally, we apply our expression filter to the analysis of de novo variants in patients with autism spectrum disorder and intellectual disability or developmental disorders to show that pLoF variants in weakly expressed regions have similar effect sizes to those of synonymous variants, whereas pLoF variants in highly expressed exons are most strongly enriched among cases. Our annotation is fast, flexible and generalizable, making it possible for any variant file to be annotated with any isoform expression dataset, and will be valuable for the genetic diagnosis of rare diseases, the analysis of rare variant burden in complex disorders, and the curation and prioritization of variants in recall-by-genotype studies.

Whole genome sequencing for mutation discovery in a single case of lysosomal storage disease (MPS type 1) in the dog

Mucopolysaccharidosis (MPS) is a metabolic storage disorder caused by the deficiency of any lysosomal enzyme required for the breakdown of glycosaminoglycans. A 15-month-old Boston Terrier presented with clinical signs consistent with lysosomal storage disease including corneal opacities, multifocal central nervous system disease and progressively worsening clinical course. Diagnosis was confirmed at necropsy based on histopathologic evaluation of multiple organs demonstrating accumulation of mucopolysaccharides. Whole genome sequencing was used to uncover a frame-shift insertion affecting the alpha-L-iduronidase (IDUA) gene (c.19_20insCGGCCCCC), a mutation confirmed in another Boston Terrier presented 2 years later with a similar clinical picture. Both dogs were homozygous for the IDUA mutation and shared coat colors not recognized as normal for the breed by the American Kennel Club. In contrast, the mutation was not detected in 120 unrelated Boston Terriers as well as 202 dogs from other breeds. Recent inbreeding to select for recessive and unusual coat colors may have concentrated this relatively rare allele in the breed. The identification of the variant enables ante-mortem diagnosis of similar cases and selective breeding to avoid the spread of this disease in the breed. Boston Terriers carrying this variant represent a promising model for MPS I with neurological abnormalities in humans.

A recurrent missense variant in HARS2 results in variable sensorineural hearing loss in three unrelated families

HARS2 encodes mitochondrial histidyl-tRNA synthetase (HARS2), which links histidine to its cognate tRNA in the mitochondrial matrix. Biallelic variants in HARS2 are associated with Perrault syndrome, a rare recessive condition characterized by sensorineural hearing loss in both sexes and primary ovarian insufficiency in 46,XX females. Some individuals with Perrault syndrome have a broader phenotypic spectrum with neurological features, including ataxia and peripheral neuropathy. Here, we report a recurrent variant in HARS2 in association with sensorineural hearing loss. In affected individuals from three unrelated families, the variant HARS2 c.1439G>A p.(Arg480His) is present as a heterozygous variant in trans to a putative pathogenic variant. The low prevalence of the allele HARS2 c.1439G>A p.(Arg480His) in the general population and its presence in three families with hearing loss, confirm the pathogenicity of this variant and illustrate the presentation of Perrault syndrome as nonsyndromic hearing loss in males and prepubertal females.