De novo missense variants in exon 9 of SEPHS1 cause a neurodevelopmental condition with developmental delay, poor growth, hypotonia, and dysmorphic features.

Selenophosphate synthetase (SEPHS) plays an essential role in selenium metabolism. Two mammalian SEPHS paralogues, SEPHS1 and SEPHS2, share high sequence identity and structural homology with SEPHS. Here, we report nine individuals from eight families with developmental delay, growth and feeding problems, hypotonia, and dysmorphic features, all with heterozygous missense variants in SEPHS1. Eight of these individuals had a recurrent variant at amino acid position 371 of SEPHS1 (p.Arg371Trp, p.Arg371Gln, and p.Arg371Gly); seven of these variants were known to be de novo. Structural modeling and biochemical assays were used to understand the effect of these variants on SEPHS1 function. We found that a variant at residue Trp352 results in local structural changes of the C-terminal region of SEPHS1 that decrease the overall thermal stability of the enzyme. In contrast, variants of a solvent-exposed residue Arg371 do not impact enzyme stability and folding but could modulate direct protein-protein interactions of SEPSH1 with cellular factors in promoting cell proliferation and development. In neuronal SH-SY5Y cells, we assessed the impact of SEPHS1 variants on cell proliferation and ROS production and investigated the mRNA expression levels of genes encoding stress-related selenoproteins. Our findings provided evidence that the identified SEPHS1 variants enhance cell proliferation by modulating ROS homeostasis. Our study supports the hypothesis that SEPHS1 plays a critical role during human development and provides a basis for further investigation into the molecular mechanisms employed by SEPHS1. Furthermore, our data suggest that variants in SEPHS1 are associated with a neurodevelopmental disorder.

A syndromic neurodevelopmental disorder caused by rare variants in PPFIA3.

PPFIA3 encodes the protein-tyrosine phosphatase, receptor-type, F-polypeptide-interacting-protein-alpha-3 (PPFIA3), which is a member of the LAR-protein-tyrosine phosphatase-interacting-protein (liprin) family involved in synapse formation and function, synaptic vesicle transport, and presynaptic active zone assembly. The protein structure and function are evolutionarily well conserved, but human diseases related to PPFIA3 dysfunction are not yet reported in OMIM. Here, we report 20 individuals with rare PPFIA3 variants (19 heterozygous and 1 compound heterozygous) presenting with developmental delay, intellectual disability, hypotonia, dysmorphisms, microcephaly or macrocephaly, autistic features, and epilepsy with reduced penetrance. Seventeen unique PPFIA3 variants were detected in 18 families. To determine the pathogenicity of PPFIA3 variants in vivo, we generated transgenic fruit flies producing either human wild-type (WT) PPFIA3 or five missense variants using GAL4-UAS targeted gene expression systems. In the fly overexpression assays, we found that the PPFIA3 variants in the region encoding the N-terminal coiled-coil domain exhibited stronger phenotypes compared to those affecting the C-terminal region. In the loss-of-function fly assay, we show that the homozygous loss of fly Liprin-α leads to embryonic lethality. This lethality is partially rescued by the expression of human PPFIA3 WT, suggesting human PPFIA3 function is partially conserved in the fly. However, two of the tested variants failed to rescue the lethality at the larval stage and one variant failed to rescue lethality at the adult stage. Altogether, the human and fruit fly data reveal that the rare PPFIA3 variants are dominant-negative loss-of-function alleles that perturb multiple developmental processes and synapse formation.

Magnetically powered microwheel thrombolysis of occlusive thrombi in zebrafish.

Tissue plasminogen activator (tPA) is the only FDA-approved treatment for ischemic stroke but carries significant risks, including major hemorrhage. Additional options are needed, especially in small vessel thrombi which account for ~25% of ischemic strokes. We have previously shown that tPA-functionalized colloidal microparticles can be assembled into microwheels (µwheels) and manipulated under the control of applied magnetic fields to enable rapid thrombolysis of fibrin gels in microfluidic models of thrombosis. Transparent zebrafish larvae have a highly conserved coagulation cascade that enables studies of hemostasis and thrombosis in the context of intact vasculature, clotting factors, and blood cells. Here, we show that tPA-functionalized µwheels can perform rapid and targeted recanalization in vivo. This effect requires both tPA and µwheels, as minimal to no recanalization is achieved with tPA alone, µwheels alone, or tPA-functionalized microparticles in the absence of a magnetic field. We evaluated tPA-functionalized µwheels in CRISPR-generated plasminogen (plg) heterozygous and homozygous mutants and confirmed that tPA-functionalized µwheels are dose-dependent on plasminogen for lysis. We have found that magnetically powered µwheels as a targeted tPA delivery system are dramatically more efficient at plasmin-mediated thrombolysis than systemic delivery in vivo. Further development of this system in fish and mammalian models could enable a less invasive strategy for alleviating ischemia that is safer than directed thrombectomy or systemic infusion of tPA.

The recurrent de novo c.2011C>T missense variant in MTSS2 causes syndromic intellectual disability.

MTSS2, also known as MTSS1L, binds to plasma membranes and modulates their bending. MTSS2 is highly expressed in the central nervous system (CNS) and appears to be involved in activity-dependent synaptic plasticity. Variants in MTSS2 have not yet been associated with a human phenotype in OMIM. Here we report five individuals with the same heterozygous de novo variant in MTSS2 (GenBank: NM_138383.2: c.2011C>T [p.Arg671Trp]) identified by exome sequencing. The individuals present with global developmental delay, mild intellectual disability, ophthalmological anomalies, microcephaly or relative microcephaly, and shared mild facial dysmorphisms. Immunoblots of fibroblasts from two affected individuals revealed that the variant does not significantly alter MTSS2 levels. We modeled the variant in Drosophila and showed that the fly ortholog missing-in-metastasis (mim) was widely expressed in most neurons and a subset of glia of the CNS. Loss of mim led to a reduction in lifespan, impaired locomotor behavior, and reduced synaptic transmission in adult flies. Expression of the human MTSS2 reference cDNA rescued the mim loss-of-function (LoF) phenotypes, whereas the c.2011C>T variant had decreased rescue ability compared to the reference, suggesting it is a partial LoF allele. However, elevated expression of the variant, but not the reference MTSS2 cDNA, led to similar defects as observed by mim LoF, suggesting that the variant is toxic and may act as a dominant-negative allele when expressed in flies. In summary, our findings support that mim is important for appropriate neural function, and that the MTSS2 c.2011C>T variant causes a syndromic form of intellectual disability.

Bridging clinical care and research in Ontario, Canada: Maximizing diagnoses from reanalysis of clinical exome sequencing data.

We examined the utility of clinical and research processes in the reanalysis of publicly-funded clinical exome sequencing data in Ontario, Canada. In partnership with eight sites, we recruited 287 families with suspected rare genetic diseases tested between 2014 and 2020. Data from seven laboratories was reanalyzed with the referring clinicians. Reanalysis of clinically relevant genes identified diagnoses in 4% (13/287); four were missed by clinical testing. Translational research methods, including analysis of novel candidate genes, identified candidates in 21% (61/287). Of these, 24 families have additional evidence through data sharing to support likely diagnoses (8% of cohort). This study indicates few diagnoses are missed by clinical laboratories, the incremental gain from reanalysis of clinically-relevant genes is modest, and the highest yield comes from validation of novel disease-gene associations. Future implementation of translational research methods, including continued reporting of compelling genes of uncertain significance by clinical laboratories, should be considered to maximize diagnoses.

Biallelic Variants in the Ectonucleotidase ENTPD1 Cause a Complex Neurodevelopmental Disorder with Intellectual Disability, Distinct White Matter Abnormalities, and Spastic Paraplegia.

OBJECTIVE: Human genomics established that pathogenic variation in diverse genes can underlie a single disorder. For example, hereditary spastic paraplegia is associated with >80 genes, with frequently only few affected individuals described for each gene. Herein, we characterize a large cohort of individuals with biallelic variation in ENTPD1, a gene previously linked to spastic paraplegia 64 (Mendelian Inheritance in Man # 615683). METHODS: Individuals with biallelic ENTPD1 variants were recruited worldwide. Deep phenotyping and molecular characterization were performed. RESULTS: A total of 27 individuals from 17 unrelated families were studied; additional phenotypic information was collected from published cases. Twelve novel pathogenic ENTPD1 variants are described (NM 001776.6): c.398_399delinsAA; p.(Gly133Glu), c.540del; p.(Thr181Leufs*18), c.640del; p.(Gly216Glufs*75), c.185 T > G; p.(Leu62*), c.1531 T > C; p.(*511Glnext*100), c.967C > T; p.(Gln323*), c.414-2_414-1del, and c.146 A > G; p.(Tyr49Cys) including 4 recurrent variants c.1109 T > A; p.(Leu370*), c.574-6_574-3del, c.770_771del; p.(Gly257Glufs*18), and c.1041del; p.(Ile348Phefs*19). Shared disease traits include childhood onset, progressive spastic paraplegia, intellectual disability (ID), dysarthria, and white matter abnormalities. In vitro assays demonstrate that ENTPD1 expression and function are impaired and that c.574-6_574-3del causes exon skipping. Global metabolomics demonstrate ENTPD1 deficiency leads to impaired nucleotide, lipid, and energy metabolism. INTERPRETATION: The ENTPD1 locus trait consists of childhood disease onset, ID, progressive spastic paraparesis, dysarthria, dysmorphisms, and white matter abnormalities, with some individuals showing neurocognitive regression. Investigation of an allelic series of ENTPD1 (1) expands previously described features of ENTPD1-related neurological disease, (2) highlights the importance of genotype-driven deep phenotyping, (3) documents the need for global collaborative efforts to characterize rare autosomal recessive disease traits, and (4) provides insights into disease trait neurobiology. ANN NEUROL 2022;92:304-321.

The implementation of an enhanced clinical model to improve the diagnostic yield of exome sequencing for patients with a rare genetic disease: A Canadian experience.

The introduction of clinical exome sequencing (ES) has provided a unique opportunity to decrease the diagnostic odyssey for patients living with a rare genetic disease (RGD). ES has been shown to provide a diagnosis in 29%-57% of patients with a suspected RGD, with as many as 70% remaining undiagnosed. There is a need to advance the clinical model of care by more formally integrating approaches that were previously considered research into an enhanced diagnostic workflow. We developed an Exome Clinic, which set out to evaluate a workflow for improving the diagnostic yield of ES for patients with an undiagnosed RGD. Here, we report the outcomes of 47 families who underwent clinical ES in the first year of the clinic. The diagnostic yield from clinical ES was 40% (19/47). Families who remained undiagnosed after ES had the opportunity for follow-up studies that included phenotyping and candidate variant segregation in relatives, genomic matchmaking, and ES reanalysis. This enhanced diagnostic workflow increased the diagnostic yield to 55% (26/47), predominantly through the resolution of variants and genes of uncertain significance. We advocate that this approach be integrated into mainstream clinical practice and highlight the importance of a coordinated translational approach for patients with RGD.

PhenomeCentral: 7 years of rare disease matchmaking.

A major challenge in validating genetic causes for patients with rare diseases (RDs) is the difficulty in identifying other RD patients with overlapping phenotypes and variants in the same candidate gene. This process, known as matchmaking, requires robust data sharing solutions to be effective. In 2014 we launched PhenomeCentral, a RD data repository capable of collecting computer-readable genotypic and phenotypic data for the purposes of RD matchmaking. Over the past 7 years PhenomeCentral's features have been expanded and its data set has consistently grown. There are currently 1615 users registered on PhenomeCentral, which have contributed over 12,000 patient cases. Most of these cases contain detailed phenotypic terms, with a significant portion also providing genomic sequence data or other forms of clinical information. Matchmaking within PhenomeCentral, and with connections to other data repositories in the Matchmaker Exchange, have collectively resulted in over 60,000 matches, which have facilitated multiple gene discoveries. The collection of deep phenotypic and genotypic data has also positioned PhenomeCentral well to support next generation of matchmaking initiatives that utilize genome sequencing data, ensuring that PhenomeCentral will remain a useful tool in solving undiagnosed RD cases in the years to come.

Genomics4RD: An integrated platform to share Canadian deep-phenotype and multiomic data for international rare disease gene discovery.

Despite recent progress in the understanding of the genetic etiologies of rare diseases (RDs), a significant number remain intractable to diagnostic and discovery efforts. Broad data collection and sharing of information among RD researchers is therefore critical. In 2018, the Care4Rare Canada Consortium launched the project C4R-SOLVE, a subaim of which was to collect, harmonize, and share both retrospective and prospective Canadian clinical and multiomic data. Here, we introduce Genomics4RD, an integrated web-accessible platform to share Canadian phenotypic and multiomic data between researchers, both within Canada and internationally, for the purpose of discovering the mechanisms that cause RDs. Genomics4RD has been designed to standardize data collection and processing, and to help users systematically collect, prioritize, and visualize participant information. Data storage, authorization, and access procedures have been developed in collaboration with policy experts and stakeholders to ensure the trusted and secure access of data by external researchers. The breadth and standardization of data offered by Genomics4RD allows researchers to compare candidate disease genes and variants between participants (i.e., matchmaking) for discovery purposes, while facilitating the development of computational approaches for multiomic data analyses and enabling clinical translation efforts for new genetic technologies in the future.

Outcome of over 1500 matches through the Matchmaker Exchange for rare disease gene discovery: The 2-year experience of Care4Rare Canada.

PURPOSE: Matchmaking has emerged as a useful strategy for building evidence toward causality of novel disease genes in patients with undiagnosed rare diseases. The Matchmaker Exchange (MME) is a collaborative initiative that facilitates international data sharing for matchmaking purposes; however, data on user experience is limited. METHODS: Patients enrolled as part of the Finding of Rare Disease Genes in Canada (FORGE) and Care4Rare Canada research programs had their exome sequencing data reanalyzed by a multidisciplinary research team over a 2-year period. Compelling variants in genes not previously associated with a human phenotype were submitted through the MME node PhenomeCentral, and outcomes were collected. RESULTS: In this study, 194 novel candidate genes were submitted to the MME, resulting in 1514 matches, and 15% of the genes submitted resulted in collaborations. Most submissions resulted in at least 1 match, and most matches were with GeneMatcher (82%), where additional email exchange was required to evaluate the match because of the lack of phenotypic or inheritance information. CONCLUSION: Matchmaking through the MME is an effective way to investigate novel candidate genes; however, it is a labor-intensive process. Engagement from the community to contribute phenotypic, genotypic, and inheritance data will ensure that matchmaking continues to be a useful approach in the future.

Causative Mutations and Mechanism of Androgenetic Hydatidiform Moles.

Androgenetic complete hydatidiform moles are human pregnancies with no embryos and affect 1 in every 1,400 pregnancies. They have mostly androgenetic monospermic genomes with all the chromosomes originating from a haploid sperm and no maternal chromosomes. Androgenetic complete hydatidiform moles were described in 1977, but how they occur has remained an open question. We identified bi-allelic deleterious mutations in MEI1, TOP6BL/C11orf80, and REC114, with roles in meiotic double-strand breaks formation in women with recurrent androgenetic complete hydatidiform moles. We investigated the occurrence of androgenesis in Mei1-deficient female mice and discovered that 8% of their oocytes lose all their chromosomes by extruding them with the spindles into the first polar body. We demonstrate that Mei1-/- oocytes are capable of fertilization and 5% produce androgenetic zygotes. Thus, we uncover a meiotic abnormality in mammals and a mechanism for the genesis of androgenetic zygotes that is the extrusion of all maternal chromosomes and their spindles into the first polar body.

Is Evolution in Response to Extreme Events Good for Population Persistence?

AbstractClimate change is predicted to increase the severity of environmental perturbations, including storms and droughts, which act as strong selective agents. These extreme events are often of finite duration (pulse disturbances). Hence, while evolution during an extreme event may be adaptive, the resulting phenotypic changes may become maladaptive when the event ends. Using individual-based models and analytic approximations that fuse quantitative genetics and demography, we explore how heritability and phenotypic variance affect population size and extinction risk in finite populations under an extreme event of fixed duration. Since more evolution leads to greater maladaptation and slower population recovery following an extreme event, greater heritability can increase extinction risk when the extreme event is short. Alternatively, when an extreme event is sufficiently long, heritability often helps a population persist. We also find that when events are severe, the buffering effect of phenotypic variance can outweigh the increased load it causes.

Heterozygous loss-of-function variants significantly expand the phenotypes associated with loss of GDF11.

PURPOSE: Growth differentiation factor 11 (GDF11) is a key signaling protein required for proper development of many organ systems. Only one prior study has associated an inherited GDF11 variant with a dominant human disease in a family with variable craniofacial and vertebral abnormalities. Here, we expand the phenotypic spectrum associated with GDF11 variants and document the nature of the variants. METHODS: We present a cohort of six probands with de novo and inherited nonsense/frameshift (4/6 patients) and missense (2/6) variants in GDF11. We generated gdf11 mutant zebrafish to model loss of gdf11 phenotypes and used an overexpression screen in Drosophila to test variant functionality. RESULTS: Patients with variants in GDF11 presented with craniofacial (5/6), vertebral (5/6), neurological (6/6), visual (4/6), cardiac (3/6), auditory (3/6), and connective tissue abnormalities (3/6). gdf11 mutant zebrafish show craniofacial abnormalities and body segmentation defects that match some patient phenotypes. Expression of the patients' variants in the fly showed that one nonsense variant in GDF11 is a severe loss-of-function (LOF) allele whereas the missense variants in our cohort are partial LOF variants. CONCLUSION: GDF11 is needed for human development, particularly neuronal development, and LOF GDF11 alleles can affect the development of numerous organs and tissues.

UBA2 variants underlie a recognizable syndrome with variable aplasia cutis congenita and ectrodactyly.

PURPOSE: The human chromosome 19q13.11 deletion syndrome is associated with a variable phenotype that includes aplasia cutis congenita (ACC) and ectrodactyly as specific features. UBA2 (ubiquitin-like modifier-activating enzyme 2) lies adjacent to the minimal deletion overlap region. We aimed to define the UBA2-related phenotypic spectrum in humans and zebrafish due to sequence variants and to establish the mechanism of disease. METHODS: Exome sequencing was used to detect UBA2 sequence variants in 16 subjects in 7 unrelated families. uba2 loss of function was modeled in zebrafish. Effects of human missense variants were assessed in zebrafish rescue experiments. RESULTS: Seven human UBA2 loss-of-function and missense sequence variants were detected. UBA2-phenotypes included ACC, ectrodactyly, neurodevelopmental abnormalities, ectodermal, skeletal, craniofacial, cardiac, renal, and genital anomalies. uba2 was expressed in zebrafish eye, brain, and pectoral fins; uba2-null fish showed deficient growth, microcephaly, microphthalmia, mandibular hypoplasia, and abnormal fins. uba2-mRNAs with human missense variants failed to rescue nullizygous zebrafish phenotypes. CONCLUSION: UBA2 variants cause a recognizable syndrome with a wide phenotypic spectrum. Our data suggest that loss of UBA2 function underlies the human UBA2 monogenic disorder and highlights the importance of SUMOylation in the development of affected tissues.

Haploid selection, sex ratio bias, and transitions between sex-determining systems.

Sex determination is remarkably dynamic; many taxa display shifts in the location of sex-determining loci or the evolution of entirely new sex-determining systems. Predominant theories for why we observe such transitions generally conclude that novel sex-determining systems are favoured by selection if they equalise the sex ratio or increase linkage with a locus that experiences different selection in males versus females. We use population genetic models to extend these theories in two ways: (1) We consider the dynamics of loci very tightly linked to the ancestral sex-determining loci, e.g., within the nonrecombining region of the ancestral sex chromosomes. Variation at such loci can favour the spread of new sex-determining systems in which the heterogametic sex changes (XY to ZW or ZW to XY) and the new sex-determining region is less closely linked (or even unlinked) to the locus under selection. (2) We consider selection upon haploid genotypes either during gametic competition (e.g., pollen competition) or meiosis (i.e., nonmendelian segregation), which can cause the zygotic sex ratio to become biased. Haploid selection can drive transitions between sex-determining systems without requiring selection to act differently in diploid males versus females. With haploid selection, we find that transitions between male and female heterogamety can evolve so that linkage with the sex-determining locus is either strengthened or weakened. Furthermore, we find that sex ratio biases may increase or decrease with the spread of new sex chromosomes, which implies that transitions between sex-determining systems cannot be simply predicted by selection to equalise the sex ratio. In fact, under many conditions, we find that transitions in sex determination are favoured equally strongly in cases in which the sex ratio bias increases or decreases. Overall, our models predict that sex determination systems should be highly dynamic, particularly when haploid selection is present, consistent with the evolutionary lability of this trait in many taxa.

PDXK mutations cause polyneuropathy responsive to pyridoxal 5'-phosphate supplementation.

OBJECTIVE: To identify disease-causing variants in autosomal recessive axonal polyneuropathy with optic atrophy and provide targeted replacement therapy. METHODS: We performed genome-wide sequencing, homozygosity mapping, and segregation analysis for novel disease-causing gene discovery. We used circular dichroism to show secondary structure changes and isothermal titration calorimetry to investigate the impact of variants on adenosine triphosphate (ATP) binding. Pathogenicity was further supported by enzymatic assays and mass spectroscopy on recombinant protein, patient-derived fibroblasts, plasma, and erythrocytes. Response to supplementation was measured with clinical validated rating scales, electrophysiology, and biochemical quantification. RESULTS: We identified biallelic mutations in PDXK in 5 individuals from 2 unrelated families with primary axonal polyneuropathy and optic atrophy. The natural history of this disorder suggests that untreated, affected individuals become wheelchair-bound and blind. We identified conformational rearrangement in the mutant enzyme around the ATP-binding pocket. Low PDXK ATP binding resulted in decreased erythrocyte PDXK activity and low pyridoxal 5'-phosphate (PLP) concentrations. We rescued the clinical and biochemical profile with PLP supplementation in 1 family, improvement in power, pain, and fatigue contributing to patients regaining their ability to walk independently during the first year of PLP normalization. INTERPRETATION: We show that mutations in PDXK cause autosomal recessive axonal peripheral polyneuropathy leading to disease via reduced PDXK enzymatic activity and low PLP. We show that the biochemical profile can be rescued with PLP supplementation associated with clinical improvement. As B6 is a cofactor in diverse essential biological pathways, our findings may have direct implications for neuropathies of unknown etiology characterized by reduced PLP levels. ANN NEUROL 2019;86:225-240.

Human trabecular meshwork cell behavior is influenced by collagen scaffold pore architecture and glycosaminoglycan composition.

Glaucoma is a degenerative eye disease in which damage to the optic nerve leads to a characteristic loss of vision. The primary risk factor for glaucoma is an increased intraocular pressure that is caused by an imbalance of aqueous humor generation and subsequent drainage through the trabecular meshwork (TM) drainage system. The small size, donor tissue limitations, and high complexity of the TM make it difficult to research the relationship between the TM cells and their immediate environment. Thus, a biomaterial-based approach may be more appropriate for research manipulations and in vitro drug development platforms. In this work, human TM (hTM) cells were cultured on various collagen scaffolds containing different glycosaminoglycans (GAGs) and different pore architectures to better understand how hTM cells respond to changes in their extracellular environment. Cellular response was measured by quantifying cellular proliferation and expression of an important extracellular matrix protein, fibronectin. The pore architecture of the scaffolds was altered using freeze-casting technique to make both large and small pores that were aligned or with a non-aligned random structure. The composition of the scaffolds was altered with the addition of chondroitin sulfate and/or hyaluronic acid. It was found that the hTM cells grown on large pore scaffolds proliferate more than those grown on small pores. There was an increase in the fibronectin expression with the incorporation of GAGs, and its morphology was changed by the underlying pore architecture. This work will help provide an insight into the behavior of hTM cells when introducing changes in their microenvironment.

Imaging and Analysis of Cellular Locations in Three-Dimensional Tissue Models.

The absence of quantitative in vitro cell-extracellular matrix models represents an important bottleneck for basic research and human health. Randomness of cellular distributions provides an opportunity for the development of a quantitative in vitro model. However, quantification of the randomness of random cell distributions is still lacking. In this paper, we have imaged cellular distributions in an alginate matrix using a multiview light sheet microscope and developed quantification metrics of randomness by modeling it as a Poisson process, a process that has constant probability of occurring in space or time. We imaged fluorescently labeled human mesenchymal stem cells embedded in an alginate matrix of thickness greater than 5 mm with axial resolution, the mean full width at half maximum of the axial intensity profiles of fluorescent particles. Simulated randomness agrees well with the experiments. Quantification of distributions and validation by simulations will enable quantitative study of cell-matrix interactions in tissue models.