Aberrant phase separation and nucleolar dysfunction in rare genetic diseases

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

A proteome-scale map of the SARS-CoV-2-human contactome

Understanding the mechanisms of coronavirus disease 2019 (COVID-19) disease severity to efficiently design therapies for emerging virus variants remains an urgent challenge of the ongoing pandemic. Infection and immune reactions are mediated by direct contacts between viral molecules and the host proteome, and the vast majority of these virus-host contacts (the 'contactome') have not been identified. Here, we present a systematic contactome map of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) with the human host encompassing more than 200 binary virus-host and intraviral protein-protein interactions. We find that host proteins genetically associated with comorbidities of severe illness and long COVID are enriched in SARS-CoV-2 targeted network communities. Evaluating contactome-derived hypotheses, we demonstrate that viral NSP14 activates nuclear factor κB (NF-κB)-dependent transcription, even in the presence of cytokine signaling. Moreover, for several tested host proteins, genetic knock-down substantially reduces viral replication. Additionally, we show for USP25 that this effect is phenocopied by the small-molecule inhibitor AZ1. Our results connect viral proteins to human genetic architecture for COVID-19 severity and offer potential therapeutic targets.

Mapping the minutiae of the human methylome

A study in Nature reports comprehensive, genome-wide mapping of the human methylome that provides mechanistic insights into gene regulation and offers great potential for analysis of cell-free DNA.

Clustering of Juvenile Canavan disease in an Indian community due to population bottleneck and isolation: genomic signatures of a founder event

Mild/juvenile Canavan disease (M/JCD) is less frequently reported in the literature and little is known about its pathogenetic mechanisms. We report a comprehensive investigation into the pathogenetic mechanism of a novel NM_000049.4(ASPA):c.526G>A variant in two families. The families belong to Telugu Devanga Chettiar community (TDC) from southern India. TDC has a complex history of migration from their historical origin centuries ago with high endogamy. TDC probably has the highest clustering M/JCD recorded historically (around 24 cases). The pathogenic variant was shown to cause non-classical splicing defect resulting in two different transcripts. The splicing aberration, a loss of function mechanism coupled with a milder missense effect can explain the milder phenotype compared to the infantile-onset CD. The high clustering of an extremely rare form of neurodegenerative disorder with reduced fitness, led us to speculate the possibility of a founder event. Genotyping array of TDC and multiple distinct populations of Indian origin for several population genetic parameters was performed. It yielded robust signatures of a founder event in TDC, such as a high fixation index, increased runs of homozygosity and identity-by-descent in the absence of consanguinity; a large haplotype with high linkage disequilibrium among markers comprising the pathogenic variant; a robust population structure; mutation dating, estimating the age of the potential founder of TDC at around 375 years; possibly a high carrier rate in TDC. This study has not only focused its attention on natural history and pathogenetics but also paves way for carrier screening programs in TDC and future therapeutic studies.

Mono- and biallelic variant effects on disease at biobank scale

Identifying causal factors for Mendelian and common diseases is an ongoing challenge in medical genetics¹. Population bottleneck events, such as those that occurred in the history of the Finnish population, enrich some homozygous variants to higher frequencies, which facilitates the identification of variants that cause diseases with recessive inheritance^2,3. Here we examine the homozygous and heterozygous effects of 44,370 coding variants on 2,444 disease phenotypes using data from the nationwide electronic health records of 176,899 Finnish individuals. We find associations for homozygous genotypes across a broad spectrum of phenotypes, including known associations with retinal dystrophy and novel associations with adult-onset cataract and female infertility. Of the recessive disease associations that we identify, 13 out of 20 would have been missed by the additive model that is typically used in genome-wide association studies. We use these results to find many known Mendelian variants whose inheritance cannot be adequately described by a conventional definition of dominant or recessive. In particular, we find variants that are known to cause diseases with recessive inheritance with significant heterozygous phenotypic effects. Similarly, we find presumed benign variants with disease effects. Our results show how biobanks, particularly in founder populations, can broaden our understanding of complex dosage effects of Mendelian variants on disease.

Whole-genome sequencing of multiple related individuals with type 2 diabetes reveals an atypical likely pathogenic mutation in the PAX6 gene

Pathogenic variants in more than 14 genes have been implicated in monogenic diabetes; however, a significant fraction of individuals with young-onset diabetes and a strong family history of diabetes have unknown genetic etiology. To identify novel pathogenic alleles for monogenic diabetes, we performed whole-genome sequencing (WGS) on four related individuals with type 2 diabetes - including one individual diagnosed at the age of 31 years - that were negative for mutations in known monogenic diabetes genes. The individuals were ascertained from a large case-control study and had a multi-generation family history of diabetes. Identity-by-descent (IBD) analysis revealed that the four individuals represent two sib-pairs that are third-degree relatives. A novel missense mutation (p.P81S) in the PAX6 gene was one of eight rare coding variants across the genome shared IBD by all individuals and was inherited from affected mothers in both sib-pairs. The mutation affects a highly conserved amino acid located in the paired-domain of PAX6 - a hotspot for missense mutations that cause aniridia and other eye abnormalities. However, no eye-related phenotype was observed in any individual. The well-established functional role of PAX6 in glucose-induced insulin secretion and the co-segregation of diabetes in families with aniridia provide compelling support for the pathogenicity of this mutation for diabetes. The mutation could be classified as "likely pathogenic" with a posterior probability of 0.975 according to the ACMG/AMP guidelines. This is the first PAX6 missense mutation that is likely pathogenic for autosomal-dominant adult-onset diabetes without eye abnormalities.

A rare genotype of biallelic mosaic variants in BCOR gene causing a bilateral ocular anterior segment dysgenesis and cataracts

Oculofaciocardiodental (OFCD) syndrome is a rare X-linked dominant syndrome characterized by the involvement of the eyes, face, teeth, and heart with variable expressivity. The syndrome is caused by loss-of-function variants in the BCOR gene located on the X chromosome. OFCD affects only females with presumed embryonic lethality among males. We report a first case of a female with biallelic mosaic variants in BCOR gene, leading to a severe ocular phenotype including anterior segment dysgenesis, cataracts, and retinal involvement. The unique condition of biallelic mosaic loss-of-function mutations leads to a variable expression of an allele with the pathogenic variant, independent of the X-Inactivation pattern. This novel mechanism of co-existent biallelic mosaicism should be suspected in unexplained severe cases of OFCD.

A novel FLNA variant in a fetus with skeletal dysplasia

Otopalatodigital spectrum disorder (OPDSD) is characterized by variable phenotypes, including skeletal dysplasia, and is caused by pathogenic variants in filamin A-encoding FLNA. FLNA variants associated with lethal OPDSD primarily alter the CH2 subdomain of the ABD of FLNA. Herein, we report a novel FLNA mutation in a fetus with severe skeletal dysplasia in a pregnant multigravida female with a history of repeated miscarriages and terminations.

Bayesian model and selection signature analyses reveal risk factors for canine atopic dermatitis

Canine atopic dermatitis is an inflammatory skin disease with clinical similarities to human atopic dermatitis. Several dog breeds are at increased risk for developing this disease but previous genetic associations are poorly defined. To identify additional genetic risk factors for canine atopic dermatitis, we here apply a Bayesian mixture model adapted for mapping complex traits and a cross-population extended haplotype test to search for disease-associated loci and selective sweeps in four dog breeds at risk for atopic dermatitis. We define 15 associated loci and eight candidate regions under selection by comparing cases with controls. One associated locus is syntenic to the major genetic risk locus (Filaggrin locus) in human atopic dermatitis. One selection signal in common type Labrador retriever cases positions across the TBC1D1 gene (body weight) and one signal of selection in working type German shepherd controls overlaps the LRP1B gene (brain), near the KYNU gene (psoriasis). In conclusion, we identify candidate genes, including genes belonging to the same biological pathways across multiple loci, with potential relevance to the pathogenesis of canine atopic dermatitis. The results show genetic similarities between dog and human atopic dermatitis, and future across-species genetic comparisons are hereby further motivated.

Genetic epidemiology of inherited retinal diseases in a large patient cohort followed at a single center in Italy

Inherited retinal diseases (IRDs) are the leading cause of vision loss in the working-age population. We performed a retrospective epidemiological study to determine the genetic basis of IRDs in a large Italian cohort (n = 2790) followed at a single referral center. We provided, mainly by next generation sequencing, potentially conclusive molecular diagnosis for 2036 patients (from 1683 unrelated families). We identified a total of 1319 causative sequence variations in 132 genes, including 353 novel variants, and 866 possibly actionable genotypes for therapeutic approaches. ABCA4 was the most frequently mutated gene (n = 535; 26.3% of solved cases), followed by USH2A (n = 228; 11.2%) and RPGR (n = 102; 5.01%). The other 129 genes had a lower contribution to IRD pathogenesis (e.g. CHM 3.5%, RHO 3.5%; MYO7A 3.4%; CRB1 2.7%; RPE65 2%, RP1 1.8%; GUCY2D 1.7%). Seventy-eight genes were mutated in five patients or less. Mitochondrial DNA variants were responsible for 2.1% of cases. Our analysis confirms the complex genetic etiology of IRDs and reveals the high prevalence of ABCA4 and USH2A mutations. This study also uncovers genetic associations with a spectrum of clinical subgroups and highlights a valuable number of cases potentially eligible for clinical trials and, ultimately, for molecular therapies.

Clinical features and genetic spectrum of NMNAT1-associated retinal degeneration

OBJECTIVES

To systematically analyse the NMNAT1 variant spectrum and frequency, the associated phenotypic characteristics, and potential genotype-phenotype correlations based on our data and literature review.

METHODS

Biallelic potential pathogenic variants (PPV) in NMNAT1 were collected from our in-house exome sequencing data. Whole-genome sequencing was conducted subsequently for patients with only one heterozygous PPV detected in NMNAT1. The clinical data were reviewed and evaluated in detail. Furthermore, the literature was reviewed for reports of NMNAT1 variants and their associated phenotypes.

RESULTS

Eleven NMNAT1 variants, including two novel variants, were detected in 8 families from our cohort. All of the 9 available patients showed generalized tapetoretinal dystrophy at an early age (88.9% in the first decade), and disciform macular atrophy was identified in six patients from five unrelated families. Among a total of 125 patients from 8 families of our cohort and 91 families reported by the available literature, 92.9% patients showed onset of disease in the first year after birth, and 89.0% patients showed visual acuity of 0.05 or lower. All of the 39 patients with fundus photos available presented disciform macular atrophy with generalized tapetoretinal dystrophy. Most (54/80, 67.5%) of causative NMNAT1 variants were missense. The most frequent variants in Caucasian and Asian population are p.E257K and p.R237C, respectively.

CONCLUSIONS

Early-onset age, disciform macular atrophy with generalized tapetoretinal dystrophy, and poor visual acuity are the typical features of NMNAT1-associated retinal degeneration. Different variant hot spots of NMNAT1 were observed in different populations.

Genetic variants associated with glaucomatous visual field loss in primary open-angle glaucoma

Primary open-angle glaucoma (POAG) is characterized by a progressive optic neuropathy with visual field loss. To investigate the genetic variants associated with visual field loss in POAG, Japanese POAG patients (n = 426) and control subjects (n = 246) were genotyped for 22 genetic variants predisposing to POAG that can be classified into those associated with intraocular pressure (IOP) elevation (IOP-related genetic variants) and optic nerve vulnerability independent of IOP (optic nerve-related genetic variants). The genetic risk score (GRS) of the 17 IOP-related and five optic nerve-related genetic variants was calculated, and the associations between the GRS and the mean deviation (MD) of automated static perimetry as an indicator of the severity of visual field loss and pattern standard deviation (PSD) as an indicator of the focal disturbance were evaluated. There was a significant association (Beta = - 0.51, P = 0.0012) between the IOP-related GRS and MD. The severity of visual field loss may depend on the magnitude of IOP elevation induced by additive effects of IOP-related genetic variants. A significant association (n = 135, Beta = 0.65, P = 0.0097) was found between the optic nerve-related, but not IOP-related, GRS and PSD. The optic nerve-related (optic nerve vulnerability) and IOP-related (IOP elevation) genetic variants may play an important role in the focal and diffuse visual field loss respectively. To our knowledge, this is the first report to show an association between additive effects of genetic variants predisposing to POAG and glaucomatous visual field loss, including severity and focal/diffuse disturbance of visual field loss, in POAG.

A novel SLC5A2 heterozygous variant in a family with familial renal glucosuria

Familial renal glucosuria (FRG) is characterized by persistent glucosuria despite normal blood glucose levels in the absence of overt tubular dysfunction. SGLT2 is a sodium-glucose cotransporter expressed in the proximal tubule; loss-of-function variants in SLC5A2 are the primary cause of FRG. Heterozygous variants have rarely been reported in Japanese individuals. Here, we identified a novel SLC5A2 heterozygous variant, c.1348G>T: p.Gly450Trp, in a Japanese family comprising two children and their father.

Dominant ARF3 variants disrupt Golgi integrity and cause a neurodevelopmental disorder recapitulated in zebrafish

Vesicle biogenesis, trafficking and signaling via Endoplasmic reticulum-Golgi network support essential developmental processes and their disruption lead to neurodevelopmental disorders and neurodegeneration. We report that de novo missense variants in ARF3, encoding a small GTPase regulating Golgi dynamics, cause a developmental disease in humans impairing nervous system and skeletal formation. Microcephaly-associated ARF3 variants affect residues within the guanine nucleotide binding pocket and variably perturb protein stability and GTP/GDP binding. Functional analysis demonstrates variably disruptive consequences of ARF3 variants on Golgi morphology, vesicles assembly and trafficking. Disease modeling in zebrafish validates further the dominant behavior of the mutants and their differential impact on brain and body plan formation, recapitulating the variable disease expression. In-depth in vivo analyses traces back impaired neural precursors' proliferation and planar cell polarity-dependent cell movements as the earliest detectable effects. Our findings document a key role of ARF3 in Golgi function and demonstrate its pleiotropic impact on development.

Five years' experience of the clinical exome sequencing in a Spanish single center

Nowadays, exome sequencing is a robust and cost-efficient genetic diagnostic tool already implemented in many clinical laboratories. Despite it has undoubtedly improved our diagnostic capacity and has allowed the discovery of many new Mendelian-disease genes, it only provides a molecular diagnosis in up to 25-30% of cases. Here, we comprehensively evaluate the results of a large sample set of 4974 clinical exomes performed in our laboratory over a period of 5 years, showing a global diagnostic rate of 24.62% (1391/4974). For the evaluation we establish different groups of diseases and demonstrate how the diagnostic rate is not only dependent on the analyzed group of diseases (43.12% in ophthalmological cases vs 16.61% in neurological cases) but on the specific disorder (47.49% in retinal dystrophies vs 24.02% in optic atrophy; 18.88% in neuropathies/paraparesias vs 11.43% in dementias). We also detail the most frequent mutated genes within each group of disorders and discuss, on our experience, further investigations and directions needed for the benefit of patients.

Analysis of clinically relevant variants from ancestrally diverse Asian genomes

Asian populations are under-represented in human genomics research. Here, we characterize clinically significant genetic variation in 9051 genomes representing East Asian, South Asian, and severely under-represented Austronesian-speaking Southeast Asian ancestries. We observe disparate genetic risk burden attributable to ancestry-specific recurrent variants and identify individuals with variants specific to ancestries discordant to their self-reported ethnicity, mostly due to cryptic admixture. About 27% of severe recessive disorder genes with appreciable carrier frequencies in Asians are missed by carrier screening panels, and we estimate 0.5% Asian couples at-risk of having an affected child. Prevalence of medically-actionable variant carriers is 3.4% and a further 1.6% harbour variants with potential for pathogenic classification upon additional clinical/experimental evidence. We profile 23 pharmacogenes with high-confidence gene-drug associations and find 22.4% of Asians at-risk of Centers for Disease Control and Prevention Tier 1 genetic conditions concurrently harbour pharmacogenetic variants with actionable phenotypes, highlighting the benefits of pre-emptive pharmacogenomics. Our findings illuminate the diversity in genetic disease epidemiology and opportunities for precision medicine for a large, diverse Asian population.

National Center Biobank Network

There are six national centers (6NCs) for advanced and specialized medicine in Japan that conduct basic and clinical research on major diseases that have a substantial impact on national health. Disease-specific bioresources and information collected by each NC are stored in a separate biobank. The National Center Biobank Network (NCBN) was established in 2011 and coordinates the biobanks and researchers of the 6NCs via an open-access database (Catalogue Database: http://www2.ncbiobank.org/Index_en ) as an efficient means of providing registered biological resources and data for use in research communities. The NCBN resources are characterized by their high-quality and rich medical information and are available for life science research and for the development of novel testing methodologies (biomarkers), new treatments, and drugs for future health care in the scope of personalized medicine through a deeper understanding of disease pathogenesis. Here, we explain the activities of the NCBN and the characteristics of the NCBN Catalogue Database.

Nuclear-embedded mitochondrial DNA sequences in 66,083 human genomes

DNA transfer from cytoplasmic organelles to the cell nucleus is a legacy of the endosymbiotic event-the majority of nuclear-mitochondrial segments (NUMTs) are thought to be ancient, preceding human speciation^1-3. Here we analyse whole-genome sequences from 66,083 people-including 12,509 people with cancer-and demonstrate the ongoing transfer of mitochondrial DNA into the nucleus, contributing to a complex NUMT landscape. More than 99% of individuals had at least one of 1,637 different NUMTs, with 1 in 8 individuals having an ultra-rare NUMT that is present in less than 0.1% of the population. More than 90% of the extant NUMTs that we evaluated inserted into the nuclear genome after humans diverged from apes. Once embedded, the sequences were no longer under the evolutionary constraint seen within the mitochondrion, and NUMT-specific mutations had a different mutational signature to mitochondrial DNA. De novo NUMTs were observed in the germline once in every 10⁴ births and once in every 10³ cancers. NUMTs preferentially involved non-coding mitochondrial DNA, linking transcription and replication to their origin, with nuclear insertion involving multiple mechanisms including double-strand break repair associated with PR domain zinc-finger protein 9 (PRDM9) binding. The frequency of tumour-specific NUMTs differed between cancers, including a probably causal insertion in a myxoid liposarcoma. We found evidence of selection against NUMTs on the basis of size and genomic location, shaping a highly heterogenous and dynamic human NUMT landscape.

Gene therapy: perspectives from young adults with Leber's congenital amaurosis

AIMS/PURPOSE

To investigate Leber congenital amaurosis (LCA) patients' expectations, decision-making processes and gene therapy-related concerns.

METHODS

Using a qualitative approach, we explored perceptions of gene therapy and clinical trials among individuals with LCA. Young adults with a clinical diagnosis of LCA were recruited through the Ocular Genetics Programme at the Hospital for Sick Children. Semi-structured interviews were conducted with ten patients and analysed following the principles of qualitative description.

RESULTS

Study participants were aware of ongoing gene therapy research trials and actively sought information regarding advances in ophthalmology and vision restoration. The majority of participants would enrol or were enrolled in a gene-replacement therapy trial, while a minority was ambivalent or would not enrol if provided an opportunity. Participants attributed different values to clinical trials, which influenced their willingness to participate. Intrinsic factors related to coping, adaptation to vision loss and resilience also influenced decision-making.

DISCUSSION

This study highlights the complex factors involved in gene-therapy-related decision-making and acts as a proponent for adopting patient-centred care strategies when counselling individuals considering gene therapy or clinical trial participation.

Intragenic inversions in NF1 gene as pathogenic mechanism in neurofibromatosis type 1

Neurofibromatosis type 1 (NF1), an autosomal dominant disorder characterized by skin pigmentary lesions and multiple cutaneous neurofibromas, is caused by neurofibromin 1 (NF1) loss of function variants. Currently, a molecular diagnosis is frequently established using a multistep protocol based on cDNA and gDNA sequence analysis and/or Multiplex Ligation-dependent Probe Amplification (MLPA) assay on genomic DNA, providing an overall detection rate of about 95-97%. The small proportion of clinically diagnosed patients, which at present do not obtain a molecular confirmation likely are mosaic, as their pathogenic variant may remain undetected due to low sensitivity of low coverage NGS approaches, or they may carry a type of pathogenic variant refractory to currently used technologies. Here, we report two unrelated patients presenting with two different inversions that disrupt the NF1 coding sequence, resulting in an NF1 phenotype. In one subject, the inversion was associated with microdeletions spanning a few NF1 exons at both breakpoints, while in the other the rearrangement did not cause exon loss, thus testing negative by MLPA assay. Considering the high proportion of repeated regions within the NF1 sequence, we propose that intragenic structural rearrangements should be considered as possible pathogenic mechanisms in patients fulfilling the NIH diagnostic criteria of NF1 but lacking of molecular confirmation and in patients with NF1 intragenic microdeletions.

Systematic Mendelian randomization using the human plasma proteome to discover potential therapeutic targets for stroke

Stroke is the second leading cause of death with substantial unmet therapeutic needs. To identify potential stroke therapeutic targets, we estimate the causal effects of 308 plasma proteins on stroke outcomes in a two-sample Mendelian randomization framework and assess mediation effects by stroke risk factors. We find associations between genetically predicted plasma levels of six proteins and stroke (P ≤ 1.62 × 10-4). The genetic associations with stroke colocalize (Posterior Probability >0.7) with the genetic associations of four proteins (TFPI, TMPRSS5, CD6, CD40). Mendelian randomization supports atrial fibrillation, body mass index, smoking, blood pressure, white matter hyperintensities and type 2 diabetes as stroke risk factors (P ≤ 0.0071). Body mass index, white matter hyperintensity and atrial fibrillation appear to mediate the TFPI, IL6RA, TMPRSS5 associations with stroke. Furthermore, thirty-six proteins are associated with one or more of these risk factors using Mendelian randomization. Our results highlight causal pathways and potential therapeutic targets for stroke.

Variable expressivity in a four-generation ACDMPV family with a non-coding hypermorphic SNV in trans to the frameshifting FOXF1 variant

Heterozygous single nucleotide variants (SNVs) or copy-number variant deletions involving FOXF1 or its distant lung-specific enhancer on chromosome 16q24.1 have been identified in 80-90% of patients with Alveolar capillary dysplasia with misalignment of pulmonary veins (ACDMPV), a lethal neonatal lung developmental disorder. We describe a four-generation family with a deceased ACDMPV neonate, her sibling from the electively terminated pregnancy, healthy mother with a history of pulmonary arterial hypertension (PAH), an unaffected aunt, an aunt deceased due to findings consistent with ACDMPV, and a reportedly unaffected grandmother, all with the frameshifting variant c.881_902dup (p.Gly302Profs*46) in FOXF1, and a deceased great-grandmother with a history of PAH. Genome sequencing analyses in the proband's unaffected mother revealed a non-coding putative regulatory SNV rs560517434-A within the lung-specific distant FOXF1 enhancer in trans to the FOXF1 frameshift mutation. Functional testing of this variant using an in vitro luciferase reporter assay showed that it increased FOXF1 promoter activity 10-fold. Our studies further demonstrate that non-coding SNVs in the FOXF1 enhancer region can rescue the lethal ACDMPV phenotype and support the compound inheritance gene dosage model.

Single nuclei transcriptomics of muscle reveals intra-muscular cell dynamics linked to dystrophin loss and rescue

In Duchenne muscular dystrophy, dystrophin loss leads to chronic muscle damage, dysregulation of repair, fibro-fatty replacement, and weakness. We develop methodology to efficiently isolate individual nuclei from minute quantities of frozen skeletal muscle, allowing single nuclei sequencing of irreplaceable archival samples and from very small samples. We apply this method to identify cell and gene expression dynamics within human DMD and mdx mouse muscle, characterizing effects of dystrophin rescue by exon skipping therapy at single nuclei resolution. DMD exon 23 skipping events are directly observed and increased in myonuclei from treated mice. We describe partial rescue of type IIa and IIx myofibers, expansion of an MDSC-like myeloid population, recovery of repair/remodeling M2-macrophage, and repression of inflammatory POSTN1 + fibroblasts in response to exon skipping and partial dystrophin restoration. Use of this method enables exploration of cellular and transcriptomic mechanisms of dystrophin loss and repair within an intact muscle environment. Our initial findings will scaffold our future work to more directly examine muscular dystrophies and putative recovery pathways.

The ATRX splicing variant c.21-1G>A is asymptomatic

The ATRX variant c.21-1G>A was detected by an exome analysis of a patient with Cockayne syndrome without alpha thalassemia X-linked intellectual disability syndrome (ATR-XS). In addition, variants in ERCC6 were detected. ATRX c.21-1G>A is localized at the splicing acceptor site of intron 1. This splicing event, NM_000489.6: c.21_133del p.S7Rfs*1, induces exon 2 deletion and early termination. The start codon in exon 3 of ATRX is presumed to produce a slightly shorter but functional ATRX protein.

A genomic deletion encompassing CRYBB2-CRYBB2P1 is responsible for autosomal recessive congenital cataracts

Here we report a consanguineous Pakistani family with multiple affected individuals with autosomal recessive congenital cataract (arCC). Exclusion analysis established linkage to chromosome 22q, and Sanger sequencing coupled with PCR-based chromosome walking identified a large homozygous genomic deletion. Our data suggest that this deletion leads to CRYBB2-CRYBB2P1 fusion, consisting of exons 1-5 of CRYBB2 and exon 6 of CRYBB2P1, the latter of which harbors the c.463 C > T (p.Gln155*) mutation, and is responsible for arCC.

The human genetic epidemiology of COVID-19

Human genetics can inform the biology and epidemiology of coronavirus disease 2019 (COVID-19) by pinpointing causal mechanisms that explain why some individuals become more severely affected by the disease upon infection by the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) virus. Large-scale genetic association studies, encompassing both rare and common genetic variants, have used different study designs and multiple disease phenotype definitions to identify several genomic regions associated with COVID-19. Along with a multitude of follow-up studies, these findings have increased our understanding of disease aetiology and provided routes for management of COVID-19. Important emergent opportunities include the clinical translatability of genetic risk prediction, the repurposing of existing drugs, exploration of variable host effects of different viral strains, study of inter-individual variability in vaccination response and understanding the long-term consequences of SARS-CoV-2 infection. Beyond the current pandemic, these transferrable opportunities are likely to affect the study of many infectious diseases.

A neurodevelopmental disorder caused by a novel de novo SVA insertion in exon 13 of the SRCAP gene

Pathogenic variants in the SRCAP (SNF2-related CREBBP activator protein) gene, which encodes a chromatin-remodeling ATPase, cause neurodevelopmental disorders including Floating Harbor syndrome (FLHS). Here, we report the discovery of a de novo transposon insertion in SRCAP exon 13 from trio genome sequencing in a 28-year-old female with failure to thrive, developmental delay, mood disorder and seizure disorder. The insertion was a full-length (~2.8 kb), antisense-oriented SVA insertion relative to the SRCAP transcript, bearing a 5' transduction and hallmarks of target-primed reverse transcription. The 20-bp 5' transduction allowed us to trace the source SVA element to an intron of a long non-coding RNA on chromosome 12, which is highly expressed in testis. RNA sequencing and qRT-PCR confirmed significant depletion of SRCAP expression and low-level exon skipping in the proband. This case highlights a novel disease-causing structural variant and the importance of transposon analysis in a clinical diagnostic setting.

A novel COL4A5 splicing variant causing X-linked Alport syndrome: A case report

Alport syndrome is a hereditary disorder characterized by renal impairment, hearing loss, and ocular symptoms and is caused by COL4A3, COL4A4, and COL4A5 mutations. Here, we report the case of 3-year-old boy with isolated hematuria detected in routine preventative urinary screening conducted in 3-year-old children. He carried a novel variant, NM_033380.3:c. 1032 + 1 G > A, which caused a splicing abnormality in COL4A5. He was diagnosed with X-linked Alport syndrome.

Arrhythmogenic right ventricular cardiomyopathy in a Japanese patient with a homozygous founder variant of DSG2 in the East Asian population

Arrhythmogenic right ventricular cardiomyopathy (ARVC) is a hereditary cardiomyopathy that results in fatal arrhythmias and heart failure. Herein, we report a Japanese patient with ARVC whose parents were blood relatives. Genetic testing identified a homozygous rare variant, c.1592T > G (p.Phe531Cys), of DSG2 that is presumed to be a founder variant among East Asians. Genetic counseling sessions with precise risk assessment and appropriate follow-up programs were provided to the patient and family members.

Common genetic variation associated with Mendelian disease severity revealed through cryptic phenotype analysis

Clinical heterogeneity is common in Mendelian disease, but small sample sizes make it difficult to identify specific contributing factors. However, if a disease represents the severely affected extreme of a spectrum of phenotypic variation, then modifier effects may be apparent within a larger subset of the population. Analyses that take advantage of this full spectrum could have substantially increased power. To test this, we developed cryptic phenotype analysis, a model-based approach that infers quantitative traits that capture disease-related phenotypic variability using qualitative symptom data. By applying this approach to 50 Mendelian diseases in two cohorts, we identify traits that reliably quantify disease severity. We then conduct genome-wide association analyses for five of the inferred cryptic phenotypes, uncovering common variation that is predictive of Mendelian disease-related diagnoses and outcomes. Overall, this study highlights the utility of computationally-derived phenotypes and biobank-scale cohorts for investigating the complex genetic architecture of Mendelian diseases.

Variants in ASPH cause exertional heat illness and are associated with malignant hyperthermia susceptibility

Exertional heat illness (EHI) and malignant hyperthermia (MH) are life threatening conditions associated with muscle breakdown in the setting of triggering factors including volatile anesthetics, exercise, and high environmental temperature. To identify new genetic variants that predispose to EHI and/or MH, we performed genomic sequencing on a cohort with EHI/MH and/or abnormal caffeine-halothane contracture test. In five individuals, we identified rare, pathogenic heterozygous variants in ASPH, a gene encoding junctin, a regulator of excitation-contraction coupling. We validated the pathogenicity of these variants using orthogonal pre-clinical models, CRISPR-edited C2C12 myotubes and transgenic zebrafish. In total, we demonstrate that ASPH variants represent a new cause of EHI and MH susceptibility.

A COL4A1 variant in a neonate with multiple intracranial hemorrhages and congenital cataracts

A 2-day-old neonate presented with seizures, multiple intracranial hemorrhages, and bilateral congenital cataracts. Targeted next-generation sequencing of the collagen type IV alpha 1 chain (COL4A1) gene revealed a heterozygous de novo missense variant (NM_001845.6:c.2291G>A/p.Gly764Asp). This missense variant adds to the compendium of COL4A1 variants and is associated with a COL4A1-related disorder.

Dynamic insights into the effects of nonsynonymous polymorphisms (nsSNPs) on loss of TREM2 function

Single nucleotide variations in Triggering Receptor Expressed on Myeloid Cells 2 (TREM2) are associated with many neurodegenerative diseases, including Nasu-Hakola disease (NHD), frontotemporal dementia (FTD), and late-onset Alzheimer's disease because they disrupt ligand binding to the extracellular domain of TREM2. However, the effects of nonsynonymous single nucleotide polymorphisms (nsSNPs) in TREM2 on disease progression remain unknown. In this study, we identified several high-risk nsSNPs in the TREM2 gene using various deleterious SNP predicting algorithms and analyzed their destabilizing effects on the ligand recognizing region of the TREM2 immunoglobulin (Ig) domain by molecular dynamics (MD) simulation. Cumulative prediction by all tools employed suggested the three most deleterious nsSNPs involved in loss of TREM2 function are rs549402254 (W50S), rs749358844 (R52C), and rs1409131974 (D104G). MD simulation showed that these three variants cause substantial structural alterations and conformational remodeling of the apical loops of the TREM2 Ig domain, which is responsible for ligand recognition. Detailed analysis revealed that these variants substantially increased distances between apical loops and induced conformation remodeling by changing inter-loop nonbonded contacts. Moreover, all nsSNPs changed the electrostatic potentials near the putative ligand-interacting region (PLIR), which suggested they might reduce specificity or loss of binding affinity for TREM2 ligands. Overall, this study identifies three potential high-risk nsSNPs in the TREM2 gene. We propose further studies on the molecular mechanisms responsible for loss of TREM2 function and the associations between TREM2 nsSNPs and neurodegenerative diseases.

A novel TP63 variant in a patient with ankyloblepharon-ectodermal defect-cleft lip/palate syndrome and Rapp-Hodgkin syndrome-like ectodermal dysplasia

Ankyloblepharon-ectodermal defect-cleft lip/palate syndrome and Rapp-Hodgkin syndrome are well-known TP63-related autosomal-dominant genetic disorders with various similar ectodermal dysplasias. In this study, whole-exome sequencing revealed a novel, potentially pathogenic TP63 nonsense variant (NM_001114980.2:c.25 C > T: p.Gln9Ter) in a patient with an atypical clinical phenotype. This variant was detected near translation initiation sites and has an effect only on ΔNp63α, the short isoform protein product of the TP63 gene.

Novel COL2A1 variants in Japanese patients with spondyloepiphyseal dysplasia congenita

Spondyloepiphyseal dysplasia congenita (SEDC) is a multisystemic skeletal disorder caused by pathogenic variants in COL2A1. Here, we report the genotype-phenotype correlations in five Japanese patients with SEDC based on their clinical and radiological findings. All five patients had novel missense variants resulting in glycine substitutions (G474V, G543E, G567S, G594R, and G1170R). Genetic testing is important for early intervention for the extraskeletal complications of SEDC. Spondyloepiphyseal dysplasia congenita (SEDC) (OMIM#183900) is an autosomal dominant chondrodysplasia characterized by disproportionate short stature, abnormal epiphyses, flattened vertebral bodies (skeletal abnormalities), and extraskeletal features, including myopia, retinal degeneration with retinal detachment, and cleft palate. SEDC is caused by a heterozygous variant in the collagen II alpha 1 (COL2A1) gene.

Association of central serous chorioretinopathy with single nucleotide polymorphisms in complement factor H gene in Iranian population

OBJECTIVES

To investigate the association of two different single nucleotide polymorphisms (SNPs) in the complement factor H (CFH) gene with central serous chorioretinopathy (CSCR) in the Iranian population.

METHODS

This is a case-control study with 95 participants in each group who were stratified according to their various ethnical variations. Primers for rs1329428 and rs3753394 polymorphisms were synthesized. DNA was extracted from peripheral blood leukocytes and underwent PCR and high-resolution melt analysis.

RESULTS

The frequency of tt, ct, and cc genotypes for rs1329428 polymorphism was 22 (26.5%), 46 (55.4%), and 15 (18.1%) in acute CSCR and 5 (41.7%), 5 (41.7%), and 2 (16.7%) in chronic CSCR respectively with no significant difference between case and control groups. The frequency of tt, ct, and cc genotypes for rs3753394 polymorphism was 31 (37.3%), 14 (16.9%), and 38 (45.8%) in acute CSCR and 4 (33.3%), 3 (25%), and 5 (41.7%) in chronic CSCR respectively. There was a significant difference between patients of Persian descent and controls in rs3753394 polymorphism (P = 0.00, chi-square test). There was no statistical difference in the frequency of polymorphism between acute and chronic patients (P = 0.64 and P = 0.79 respectively, chi-square test).

CONCLUSIONS

The rs3753394 polymorphism is probably associated with CSCR in Persian ethnicity. Further studies are required to validate the implications of this finding in clinical practice.

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

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

Neonatal-lethal dilated cardiomyopathy due to a homozygous LMOD2 donor splice-site variant

Dilated cardiomyopathy (DCM) is characterized by cardiac enlargement and impaired ventricular contractility leading to heart failure. A single report identified variants in leiomodin-2 (LMOD2) as a cause of neonatally-lethal DCM. Here, we describe two siblings with DCM who died shortly after birth due to heart failure. Exome sequencing identified a homozygous LMOD2 variant in both siblings, (GRCh38)chr7:g.123656237G > A; NM_207163.2:c.273 + 1G > A, ablating the donor 5' splice-site of intron-1. Pre-mRNA splicing studies and western blot analysis on cDNA derived from proband cardiac tissue, MyoD-transduced proband skin fibroblasts and HEK293 cells transfected with LMOD2 gene constructs established variant-associated absence of canonically spliced LMOD2 mRNA and full-length LMOD2 protein. Immunostaining of proband heart tissue unveiled abnormally short actin-thin filaments. Our data are consistent with LMOD2 c.273 + 1G > A abolishing/reducing LMOD2 transcript expression by: (1) variant-associated perturbation in initiation of transcription due to ablation of the intron-1 donor; and/or (2) degradation of aberrant LMOD2 transcripts (resulting from use of alternative transcription start-sites or cryptic splice-sites) by nonsense-mediated decay. LMOD2 expression is critical for life and the absence of LMOD2 is associated with thin filament shortening and severe cardiac contractile dysfunction. This study describes the first splice-site variant in LMOD2 and confirms the role of LMOD2 variants in DCM.

Prevalence and natural history of schwannomas in neurofibromatosis type 2 (NF2): the influence of pathogenic variants

This study explores the natural history of vestibular, trigeminal and lower cranial nerve schwannomas (VS, TS, LCNS) in patients with Neurofibromatosis type 2 (NF2), to understand how pathogenic variants (PVs) of the NF2 gene affect tumour burden and growth rate, via a retrospective analysis of a UK NF2 centre database and imaging. VS, TS and LCNS location and size were measured in accordance with a standardised protocol. PVs were categorised in accordance with the UK NF2 Genetic Severity Score (GSS). 153 patients (age 5-82) had 458 schwannomas, of which 362 were previously untreated comprising: 204 VS, 93 TS, and 65 LCNS (IX, X, XI). 322 schwannomas had sequential imaging allowing growth rate analysis with a mean follow-up of 45 months. VS were universally present, and bilateral in 146/153 cases. 65% of tumours grew >2 mm during the study period at mean rate 2.0 mm/year. Significant association was found between increasing GSS and growth rate. TS occurred in 66/153 patients (bilateral in 27/153); 31% of tumours showed growth (mean 1.8 mm/yr). Significant increase in tumour prevalence was noted with increasing GSS. LCNS were found in 47/153 patients (bilateral in 19/153); 27% of tumours showed growth (mean 1.9 mm/yr). The trend for increased prevalence with increasing GSS did not reach significance. VS growth rate was significantly influenced by GSS and they were much more likely to grow than TS and LCNS. TS prevalence also correlated with increasing GSS.

Singleton exome sequencing of 90 fetuses with ultrasound anomalies revealing novel disease-causing variants and genotype-phenotype correlations

Exome sequencing has been increasingly implemented in prenatal genetic testing for fetuses with morphological abnormalities but normal rapid aneuploidy detection and microarray analysis. We present a retrospective study of 90 fetuses with different abnormal ultrasound findings, in which we employed the singleton exome sequencing (sES; 75 fetuses) or to a lesser extent (15 fetuses) a multigene panel analysis of 6713 genes as a primary tool for the detection of monogenic diseases. The detection rate of pathogenic or likely pathogenic variants in this study was 34.4%. The highest diagnostic rate of 56% was in fetuses with multiple anomalies, followed by cases with skeletal or renal abnormalities (diagnostic rate of 50%, respectively). We report 20 novel disease-causing variants in different known disease-associated genes and new genotype-phenotype associations for the genes KMT2D, MN1, CDK10, and EXOC3L2. Based on our data, we postulate that sES of fetal index cases with a concurrent sampling of parental probes for targeted testing of the origin of detected fetal variants could be a suitable tool to obtain reliable and rapid prenatal results, particularly in situations where a trio analysis is not possible.

Additional findings of tibial dysplasia in a male with orofaciodigital syndrome type XVI

We describe the case of a male patient with orofaciodigital (OFD) syndrome type XVI with a homozygous variant of TMEM107 (p.Phe106del) and the additional findings of tibial dysplasia, which is a pivotal finding of OFD syndrome type IV. His family history included two fetuses with anencephaly with or without cleft lip/palate and polydactyly with no genetic information. Careful attention should be given to the interpretation of this rare pattern.

NeuroSCORE is a genome-wide omics-based model that identifies candidate disease genes of the central nervous system

To identify candidate disease genes of central nervous system (CNS) phenotypes, we created the Neurogenetic Systematic Correlation of Omics-Related Evidence (NeuroSCORE). We identified five genome-wide metrics highly associated with CNS phenotypes to score 19,601 protein-coding genes. Genes scored one point per metric (range: 0-5), identifying 8298 scored genes (scores ≥ 1) and 1601 "high scoring" genes (scores ≥ 3). Using logistic regression, we determined the odds ratio that genes with a NeuroSCORE from 1 to 5 would be associated with known CNS-related phenotypes compared to genes that scored zero. We tested NeuroSCORE using microarray copy number variants (CNVs) in case-control cohorts and aggregate mouse model data. High scoring genes are associated with CNS phenotypes (OR = 5.5, p < 2E-16), enriched in case CNVs, and mouse ortholog genes that cause behavioral and nervous system abnormalities. We identified 1058 high scoring genes with no disease association in OMIM. Transforming the logistic regression results indicates high scoring genes have an 84-92% chance of being associated with a CNS phenotype. Top scoring genes include GRIA1, MAP4K4, SF1, TNPO2, and ZSWIM8. Finally, we interrogated CNVs in the Clinical Genome Resource, finding the majority of clinically significant CNVs contain high scoring genes. These findings can direct future research and improve molecular diagnostics.

Two novel heterozygous variants in ATP1A3 cause movement disorders

Variants in ATP1A3 cause neuropsychiatric disorders, especially those characterized by movement disorders. In this study, we performed whole exome sequencing for two patients with movement disorders and identified two novel heterozygous ATP1A3 variants, a missense c.2408G>A variant and an indel c.2672_2688+10delinsCAG variant. The unique indel variant occurred at the exon-intron boundary at the 3' end of exon 19, and mRNA analysis revealed that this variant caused in-frame indel alteration at the Ser891_Trp896 residue.

A neomorphic variant in SP7 alters sequence specificity and causes a high-turnover bone disorder

SP7/Osterix is a transcription factor critical for osteoblast maturation and bone formation. Homozygous loss-of-function mutations in SP7 cause osteogenesis imperfecta type XII, but neomorphic (gain-of-new-function) mutations of SP7 have not been reported in humans. Here we describe a de novo dominant neomorphic missense variant (c.926 C > G:p.S309W) in SP7 in a patient with craniosynostosis, cranial hyperostosis, and long bone fragility. Histomorphometry shows increased osteoblasts but decreased bone mineralization. Mice with the corresponding variant also show a complex skeletal phenotype distinct from that of Sp7-null mice. The mutation alters the binding specificity of SP7 from AT-rich motifs to a GC-consensus sequence (typical of other SP family members) and produces an aberrant gene expression profile, including increased expression of Col1a1 and endogenous Sp7, but decreased expression of genes involved in matrix mineralization. Our study identifies a pathogenic mechanism in which a mutation in a transcription factor shifts DNA binding specificity and provides important in vivo evidence that the affinity of SP7 for AT-rich motifs, unique among SP proteins, is critical for normal osteoblast differentiation.

Alkaptonuria in Russia

Alkaptonuria is characterized by the accumulation of homogentisic acid (HGA), part of which is excreted in the urine but the excess HGA forms a dark brown ochronotic pigment that deposits in the connective tissue (ochronosis), eventually leading to early-onset severe arthropathy. We analyzed a cohort of 48 Russian AKU families by sequencing all 14 exons (including flanking intronic sequences) of the homogentisate 1,2-dioxygenase gene (HGD) and Multiplex Ligation-dependent Probe Amplification (MLPA) analysis. Nine novel likely pathogenic HGD variants were identified, which have not been reported previously in any other country. Recently, Bychkov et al. [1] reported on the variant spectrum in another cohort of 49 Russian AKU patients. Here we summarize complete data from both cohorts that include 82 Russian AKU families. Taken together, 31 different HGD variants were found in these patients, of which 14 are novel and found only in Russia. The most common variant was c.481G>A (p.(Gly161Arg)), present in almost 54% of all AKU alleles.

Performance comparison of different classification algorithms applied to the diagnosis of familial hypercholesterolemia in paediatric subjects

Familial Hypercholesterolemia (FH) is an inherited disorder of lipid metabolism, characterized by increased low density lipoprotein cholesterol (LDLc) levels. The main purpose of the current work was to explore alternative classification methods to traditional clinical criteria for FH diagnosis, based on several biochemical and biological indicators. Logistic regression (LR), decision tree (DT), random forest (RF) and naive Bayes (NB) algorithms were developed for this purpose, and thresholds were optimized by maximization of Youden index (YI). All models presented similar accuracy (Acc), specificity (Spec) and positive predictive values (PPV). Sensitivity (Sens) and G-mean values were significantly higher in LR and RF models, compared to the DT. When compared to Simon Broome (SB) biochemical criteria for FH diagnosis, all models presented significantly higher Acc, Spec and G-mean values (p < 0.01), and lower negative predictive value (NPV, p < 0.05). Moreover, LR and RF models presented comparable Sens values. Adjustment of the cut-off point by maximizing YI significantly increased Sens values, with no significant loss in Acc. The obtained results suggest such classification algorithms can be a viable alternative to be used as a widespread screening method. An online application has been developed to assess the performance of the LR model in a wider population.

MCAD deficiency caused by compound heterozygous pathogenic variants in ACADM

Medium-chain acyl-coenzyme A dehydrogenase (MCAD) deficiency is an autosomal recessive disease caused by biallelic pathogenic ACADM variants. We report a case of an asymptomatic Japanese girl with MCAD deficiency caused by compound heterozygous pathogenic variants (NM_000016.5:c.1040G > T (p.Gly347Val) and c.449_452delCTGA (p.Thr150ArgfsTer4)). Because the MCAD residual activity in lymphocytes of the patient was below the limit of quantification, both variants are likely to cause complete loss of MCAD enzymatic activity.